- Open Access
Genome-wide mutation analysis of Helicobacter pylori after inoculation to Mongolian gerbils
Gut Pathogens volume 11, Article number: 45 (2019)
Helicobacter pylori is a pathogenic bacterium that causes various gastrointestinal diseases in the human stomach. H. pylori is well adapted to the human stomach but does not easily infect other animals. As a model animal, Mongolian gerbils are often used, however, the genome of the inoculated H. pylori may accumulate mutations to adapt to the new host. To investigate mutations occurring in H. pylori after infection in Mongolian gerbils, we compared the whole genome sequence of TN2 wild type strain (TN2wt) and next generation sequencing data of retrieved strains from the animals after different lengths of infection.
We identified mutations in 21 loci of 17 genes of the post-inoculation strains. Of the 17 genes, five were outer membrane proteins that potentially influence on the colonization and inflammation. Missense and nonsense mutations were observed in 15 and 6 loci, respectively. Multiple mutations were observed in three genes. Mutated genes included babA, tlpB, and gltS, which are known to be associated with adaptation to murine. Other mutations were involved with chemoreceptor, pH regulator, and outer membrane proteins, which also have potential to influence on the adaptation to the new host.
We confirmed mutations in genes previously reported to be associated with adaptation to Mongolian gerbils. We also listed up genes that mutated during the infection to the gerbils, though it needs experiments to prove the influence on adaptation.
Helicobacter pylori (H. pylori) is known to a risk factor of various gastrointestinal diseases [1,2,3,4]. Previous studies investigated genetic diversification of H. pylori in the time course of chronic infection or transmission and revealed that the mutation rate of this bacterium is high [5,6,7,8].
However, H. pylori is well adapted to the human stomach but does not easily infect other animals. In search of a good animal model, experimental infection was attempted in Rhesus monkeys , mice , and Mongolian gerbils [11,12,13,14,15,16]. Genetic diversification of H. pylori in the infected animals was also studied [9, 10, 16,17,18,19].
Model animals are expected to respond to the stimulation in the similar manner to humans and be maintained on reasonable cost and handling efforts. Small rodent Mongolian gerbils develop similar symptoms to human by H. pylori infection as gastric inflammation, ulceration and cancer [13, 15, 20, 21]. Thus, they work as the good animal model.
We also used Mongolian gerbils as the model animal and discovered that babA expression in H. pylori initially increased upon infection but reduced over time, then lost after 6 months  and that infection with oipA or babA mutants resulted in significantly reduced cytokine levels but alpAB mutant did not infect Mongolian gerbils .
Earlier studies used PCR to investigate changes in genes during animal infection. However, DNA sequencing advancements enabled the extensive exploration of mutations by sequencing bacterial genomes before and after infection [16, 19]. Here, we used the whole genome sequence of TN2 wild type (TN2wt) as a reference and sequenced short reads from three derivative strains to identify genomic mutations during infection in Mongolian gerbils. We detected mutations in agreement with previous studies and identified new mutations that may be associated with adaptation of the bacteria to different hosts.
Inoculation, euthanasia and isolation of H. pylori
We inoculated TN2 wild type strain  (TNwt) to Mongolian gerbils as described in our previous paper . Six-week-old male Mongolian gerbils (MGS/Sea; Harlan Sprague Dawley) were orogastrically inoculated 3 times (days 0, 1, 2) with 1.0 mL of H. pylori (109 colon-forming units/mL) or sterile brain–heart infusion (BHI) broth using gastric intubation needles after 16 h of fasting . No specific pretreatments were administered prior to orogastric H. pylori inoculation. Inoculated Mongolian gerbils were sacrificed after 1 month (TN2-1M), 3 months (TN2-3M), and 6 months (TN2-6M). At necropsy, an ~ 1-mm2 piece of gastric mucosa from the antrum was collected for culturing of H. pylori and subsequent DNA extraction.
Bacterial culture and DNA extraction
Helicobacter pylori were cultured on confluent plates expanded from a single colony under microaerobic conditions (12% CO2) at 37 °C. Bacterial DNA was extracted from the plates using a commercially available kit (QIAGEN Inc., Valencia, CA, USA).
Sequencing of the genomic DNA
The whole genome sequence of TN2wt was provided by our collaborator at the Okinawa Institute of Advanced Sciences. The whole-genome sequencing of TN2wt was carried out using the PacBio RS II (Pacific Biosciences, Menlo Park, CA) platform. De-novo assembly was performed using the hierarchical genome assembly process (HGAP) workflow , including consensus polishing with Quiver v. 2.3.3. By this workflow, the complete genome sequence of TN2wt was obtained. Annotation was performed by MiGap service provided by National Institute of Genetics. The genome DNA of H. pylori strains retrieved from the Mongolian gerbils were sequenced by HiSeq2000 (paired end, 2 × 100 bp). DNA was quantified by Qubit fluorometric method (Thermo Fisher Scientific). DNA purity was assessed by the UV absorbance ratio at 260/280 with 1.8–2.0. Finally, 500 ng of DNA input was used for DNA library preparation. The numbers of reads obtained were 13,574,248, 14,583,596, and 13,938,018 for TN2-1M, TN2-3M, and TN2-6M, respectively; 99.69%, 99.74%, and 99.75% of the reads mapped to the reference TN2wt genome, resulting in average mapping depths of 758.8, 815.7, and 779.6 for TN2-1M, TN2-3M, and TN2-6M, respectively. The coverage of the reference genome was 100% in the all strains.
Short read data of genomic DNA from the retrieved strains (TN2-1M, TN2-3M, and TN2-6M) were mapped to the complete genome sequence of TN2wt using Genomics Workbench v. 7.0.4 (CLC QIAGEN) with default parameter setting. We also attempted de-novo assembly, but the assembly produced around 30 contigs and the total length was shorter than the original genome. Therefore, we used the reference mapping results for the analysis. We selected non-synonymous mutations that were identified in more than 90% of the mapped reads. If available, protein structure data were downloaded from PDB (https://www.rcsb.org/) [24, 25] and the location of the mutated locus was visualized by Chimera v. 1.10.2 .
Results and discussion
Non-synonymous mutations in the retrieved strains
Compared with the original TN2wt genome, strains TN2-1M, TN2-3M, and TN2-6M had 6, 9, and 6 non-synonymous mutations, respectively (Table 1, Fig. 1). These mutations were resided in 17 genes. In accordance with our previous report , 5 of the 17 genes were outer membrane proteins that potentially influence on colonization and inflammation.
Some genes had multiple mutations. TN2-1M had two missense mutations in kefB and single missense mutation in other three genes. A nucleotide insertion in hofH of TN2-1M (1290th nucleotide in the gene) caused frameshift, however, it did not cause a premature stop codon. Instead, the frameshift delayed the occurrence of a stop codon and elongated the gene 15 bp. Consequently, mutations observed in TN2-1M were all missense. KefB is a component of potassium ion (K+) transportation system that regulates cytoplasmic pH and influence on bacterial growth and survival . UreI is a pH-gated urea channel that enable H. pylori to colonize in acidic environment [28,29,30]. Missense mutations in these genes might change reactivity to pH fluctuation. GltS is a Glu-specific transporter and known also to be essential for colonization of H. pylori in Mongolian gerbils [31, 32].
TN2-3M contained seven missense and two nonsense mutations. Nucleotide deletion in oppA that leads to the premature stop codon was observed both in TN2-3M and TN2-6M. OppA is one of the ABC-type transporter genes for oligopeptide transport. Previous in-vitro study reported that disruption of oppA did not significantly change the growth of the mutant from the wild type . This may suggest that the nonsense mutation in oppA was allowed because this gene is not essential for growth. Another possibility is that loss of oppA is neutral in vitro or in the originated human stomach but rather advantageous in the Mongolian gerbil stomach. Considering that the nonsense mutation of oppA was observed both in TN2-3M and TN2-6M, the latter hypothesis is also probable.
TN2-6M contained two missense and four nonsense mutations. In this strain, babA, oppA, tlpB, and outer membrane protein had nonsense mutations. As for tlpB, two missense mutations were also observed in TN2-3M. TlpB and babA are known to be involved with H. pylori adaptation to Mongolian gerbils. Our previous study revealed that infection with mutated babA reduced cytokine levels and inflammatory cell infiltrations of the host  and that babA expression disappeared 6 months after inoculation to Mongolian gerbils . TlpB is a chemoreceptor that detect acidity and urea [34, 35]. Similar to babA, mutants lacking tlpB colonized as good as wild type but caused less inflammations in the stomach of mice and Mongolian gerbils [36, 37]. TlpB accepts posttranslational regulation by small RNA that targets guanin repeat (G-repeat) upstream of the gene . Because expression of tlpB is affected by the G-repeat length, we counted the G-repeat length of our strain. The lengths were 12 for TN2wt, TN2-1M, and TN2-6M and 11 for TN2-3M, which are associated with low level of tlpB expression .
Mutations in oppA and tlpB have also been reported  (Table 2), but the inoculated animal in this study was a mouse. There were no genes in common with another genome study using the Mongolian gerbil as a model . Another research group compared the H. pylori genome before (PMSS1) and after (SS1) inoculation . They reported that oppA was disrupted in the original strain; we also observed disruption of this gene in the derived strains. The authors also reported a change at the 443rd amino acid in tlpB. Although the details of the mutations were different, these genes may be associated with the host change, since they were observed in independent studies, which occurs rarely by chance.
We previously performed a PCR-based study  wherein we examined 20 samples of Mongolian gerbils inoculated with H. pylori. TM2-6M is one of the strains used in the study. Although the disruption of babA by nucleotide deletion/insertion was observed in half of the samples, the deletion/insertion locations and lengths were different. The frequency of disrupted babA increased over time after inoculation. This suggested a possible advantage to losing babA.
Apart from babA, increasing number of nonsense mutations were observed in the current study. The frequencies of the nonsense mutations were 0/6, 2/9, and 3/6 in TN2-1M, TN2-3M, and TN2-6M. Disruption of a gene will not be desirable for the bacteria in its native environment, but it may be selected for if it is advantageous in a new environment. Gene disruption also occurs more easily than gain of a new function by substitution because genes can be broken in various ways, like in babA.
Mutated loci on the protein structure
Protein structure data were available for ureI (3UX4)  and virB11 (1NLZ) . We downloaded the data and marked the mutated loci on the structure.
UreI channel consists of six protomers that form a hexametric ring. Figure 2 shows the half of the hexametric ring and the location of H131R in each protomer. H131 is located in periplasmic loop 2 (PL2). Previous study substituted amino acids of various loci in PL2 and reported that H131R hampered urea transportation in Xenopus laevis oocytes . Figure 3 shows the location of H314Y in VirB11. VirB11 also form a hexametric assembly. H314Y is located in a b-sheet near the end of the protomer, however, no function is reported about this locus.
Structure data of TlpB was also available but G26W and G275W were outside of the analyzed region. According to protein domain information, G26W is contained in the tm1 (transmembrane helices 1) and G275W is in HAMP (histidine kinase, adenylyl cyclase, methyl-binding protein, phosphatase) domain. Tm1 mediates signal transmission across the membrane by piston-like motion of tm2 relative to tm1. HAMP domain is supposed to constitutes a switch region that translates the piston-like motion into a different type of transition within the distal portions . Therefore, mutations G26W and G275W may influence on the function of the chemoreceptor for acidity and urea.
We compared H. pylori genomes between original TN2wt and three strains retrieved after inoculation to Mongolian gerbils. We identified mutations in 21 loci of 17 genes of the post-inoculation strains. Mutated genes included babA, tlpB, and gltS, which is known to be associated with adaptation to murine. Other mutations were involved with chemoreceptor, pH regulator, and outer membrane proteins, which also have potential to influence on the adaptation to the new host.
Availability of data and materials
Genome sequence data of TN2wt is available from GenBank under the Accession number AP019730.
Cover TL, Blaser MJ. Helicobacter pylori in health and disease. Gastroenterology. 2009;136(6):1863–73.
Ernst PB, Gold BD. The disease spectrum of Helicobacter pylori: the immunopathogenesis of gastroduodenal ulcer and gastric cancer. Annu Rev Microbiol. 2000;54:615–40.
Basso D, Plebani M, Kusters JG. Pathogenesis of Helicobacter pylori infection. Helicobacter. 2010;15(Suppl 1):14–20.
Yamaoka Y, Graham DY. Helicobacter pylori virulence and cancer pathogenesis. Future Oncol. 2014;10(8):1487–500.
Linz B, Windsor HM, Gajewski JP, Hake CM, Drautz DI, Schuster SC, Marshall BJ. Helicobacter pylori genomic microevolution during naturally occurring transmission between adults. PLoS ONE. 2013;8(12):e82187.
Linz B, Windsor HM, McGraw JJ, Hansen LM, Gajewski JP, Tomsho LP, Hake CM, Solnick JV, Schuster SC, Marshall BJ. A mutation burst during the acute phase of Helicobacter pylori infection in humans and rhesus macaques. Nat Commun. 2014;5:4165.
Kennemann L, Didelot X, Aebischer T, Kuhn S, Drescher B, Droege M, Reinhardt R, Correa P, Meyer TF, Josenhans C, et al. Helicobacter pylori genome evolution during human infection. Proc Natl Acad Sci USA. 2011;108(12):5033–8.
Didelot X, Nell S, Yang I, Woltemate S, van der Merwe S, Suerbaum S. Genomic evolution and transmission of Helicobacter pylori in two South African families. Proc Natl Acad Sci USA. 2013;110(34):13880–5.
Solnick JV, Hansen LM, Salama NR, Boonjakuakul JK, Syvanen M. Modification of Helicobacter pylori outer membrane protein expression during experimental infection of rhesus macaques. Proc Natl Acad Sci USA. 2004;101(7):2106–11.
Yamaoka Y, Kita M, Kodama T, Imamura S, Ohno T, Sawai N, Ishimaru A, Imanishi J, Graham DY. Helicobacter pylori infection in mice: role of outer membrane proteins in colonization and inflammation. Gastroenterology. 2002;123(6):1992–2004.
Styer CM, Hansen LM, Cooke CL, Gundersen AM, Choi SS, Berg DE, Benghezal M, Marshall BJ, Peek RM Jr, Boren T, et al. Expression of the BabA adhesin during experimental infection with Helicobacter pylori. Infect Immun. 2010;78(4):1593–600.
Ohno T, Vallstrom A, Rugge M, Ota H, Graham DY, Arnqvist A, Yamaoka Y. Effects of blood group antigen-binding adhesin expression during Helicobacter pylori infection of Mongolian gerbils. J Infect Dis. 2011;203(5):726–35.
Hirayama F, Takagi S, Yokoyama Y, Iwao E, Ikeda Y. Establishment of gastric Helicobacter pylori infection in Mongolian gerbils. J Gastroenterol. 1996;31(Suppl 9):24–8.
Ikeno T, Ota H, Sugiyama A, Ishida K, Katsuyama T, Genta RM, Kawasaki S. Helicobacter pylori-induced chronic active gastritis, intestinal metaplasia, and gastric ulcer in Mongolian gerbils. Am J Pathol. 1999;154(3):951–60.
Ogura K, Maeda S, Nakao M, Watanabe T, Tada M, Kyutoku T, Yoshida H, Shiratori Y, Omata M. Virulence factors of Helicobacter pylori responsible for gastric diseases in Mongolian gerbil. J Exp Med. 2000;192(11):1601–10.
Beckett AC, Loh JT, Chopra A, Leary S, Lin AS, McDonnell WJ, Dixon B, Noto JM, Israel DA, Peek RM Jr, et al. Helicobacter pylori genetic diversification in the Mongolian gerbil model. PeerJ. 2018;6:e4803.
Behrens W, Schweinitzer T, Bal J, Dorsch M, Bleich A, Kops F, Brenneke B, Didelot X, Suerbaum S, Josenhans C. Role of energy sensor TlpD of Helicobacter pylori in gerbil colonization and genome analyses after adaptation in the gerbil. Infect Immun. 2013;81(10):3534–51.
Harris AG, Wilson JE, Danon SJ, Dixon MF, Donegan K, Hazell SL. Catalase (KatA) and KatA-associated protein (KapA) are essential to persistent colonization in the Helicobacter pylori SS1 mouse model. Microbiology. 2003;149(Pt 3):665–72.
Draper JL, Hansen LM, Bernick DL, Abedrabbo S, Underwood JG, Kong N, Huang BC, Weis AM, Weimer BC, van Vliet AH, et al. Fallacy of the unique genome: sequence diversity within single Helicobacter pylori strains. MBio. 2017. 10.1128/mBio.02321-16.
Noto JM, Gaddy JA, Lee JY, Piazuelo MB, Friedman DB, Colvin DC, Romero-Gallo J, Suarez G, Loh J, Slaughter JC, et al. Iron deficiency accelerates Helicobacter pylori-induced carcinogenesis in rodents and humans. J Clin Invest. 2013;123(1):479–92.
Watanabe T, Tada M, Nagai H, Sasaki S, Nakao M. Helicobacter pylori infection induces gastric cancer in mongolian gerbils. Gastroenterology. 1998;115(3):642–8.
Sugimoto M, Ohno T, Graham DY, Yamaoka Y. Helicobacter pylori outer membrane proteins on gastric mucosal interleukin 6 and 11 expression in Mongolian gerbils. J Gastroenterol Hepatol. 2011;26(11):1677–84.
Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, Clum A, Copeland A, Huddleston J, Eichler EE, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10(6):563–9.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE. The Protein Data Bank. Nucleic Acids Res. 2000;28(1):235–42.
Burley SK, Berman HM, Bhikadiya C, Bi C, Chen L, Di Costanzo L, Christie C, Dalenberg K, Duarte JM, Dutta S, et al. RCSB Protein Data Bank: biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy. Nucleic Acids Res. 2019;47(D1):D464–D474474.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.
Healy J, Ekkerman S, Pliotas C, Richard M, Bartlett W, Grayer SC, Morris GM, Miller S, Booth IR, Conway SJ, et al. Understanding the structural requirements for activators of the Kef bacterial potassium efflux system. Biochemistry. 2014;53(12):1982–92.
Skouloubris S, Thiberge JM, Labigne A, De Reuse H. The Helicobacter pylori UreI protein is not involved in urease activity but is essential for bacterial survival in vivo. Infect Immun. 1998;66(9):4517–21.
Scott D, Weeks D, Melchers K, Sachs G. UreI-mediated urea transport in Helicobacter pylori: an open and shut case? Trends Microbiol. 2000;8(8):348–9.
Bury-Mone S, Skouloubris S, Labigne A, De Reuse H. The Helicobacter pylori UreI protein: role in adaptation to acidity and identification of residues essential for its activity and for acid activation. Mol Microbiol. 2001;42(4):1021–34.
Kavermann H, Burns BP, Angermuller K, Odenbreit S, Fischer W, Melchers K, Haas R. Identification and characterization of Helicobacter pylori genes essential for gastric colonization. J Exp Med. 2003;197(7):813–22.
Leduc D, Gallaud J, Stingl K, de Reuse H. Coupled amino acid deamidase-transport systems essential for Helicobacter pylori colonization. Infect Immun. 2010;78(6):2782–92.
Weinberg MV, Maier RJ. Peptide transport in Helicobacter pylori: roles of dpp and opp systems and evidence for additional peptide transporters. J Bacteriol. 2007;189(9):3392–402.
Huang JY, Goers Sweeney E, Guillemin K, Amieva MR. Multiple acid sensors control Helicobacter pylori colonization of the stomach. PLoS Pathog. 2017;13(1):e1006118.
Huang JY, Sweeney EG, Sigal M, Zhang HC, Remington SJ, Cantrell MA, Kuo CJ, Guillemin K, Amieva MR. Chemodetection and destruction of host urea allows Helicobacter pylori to locate the epithelium. Cell Host Microbe. 2015;18(2):147–56.
Williams SM, Chen YT, Andermann TM, Carter JE, McGee DJ, Ottemann KM. Helicobacter pylori chemotaxis modulates inflammation and bacterium-gastric epithelium interactions in infected mice. Infect Immun. 2007;75(8):3747–57.
McGee DJ, Langford ML, Watson EL, Carter JE, Chen YT, Ottemann KM. Colonization and inflammation deficiencies in Mongolian gerbils infected by Helicobacter pylori chemotaxis mutants. Infect Immun. 2005;73(3):1820–7.
Pernitzsch SR, Tirier SM, Beier D, Sharma CM. A variable homopolymeric G-repeat defines small RNA-mediated posttranscriptional regulation of a chemotaxis receptor in Helicobacter pylori. Proc Natl Acad Sci USA. 2014;111(4):E501–510.
Strugatsky D, McNulty R, Munson K, Chen CK, Soltis SM, Sachs G, Luecke H. Structure of the proton-gated urea channel from the gastric pathogen Helicobacter pylori. Nature. 2013;493(7431):255–8.
Savvides SN, Yeo HJ, Beck MR, Blaesing F, Lurz R, Lanka E, Buhrdorf R, Fischer W, Haas R, Waksman G. VirB11 ATPases are dynamic hexameric assemblies: new insights into bacterial type IV secretion. EMBO J. 2003;22(9):1969–80.
Weeks DL, Sachs G. Sites of pH regulation of the urea channel of Helicobacter pylori. Mol Microbiol. 2001;40(6):1249–59.
Goers Sweeney E, Henderson JN, Goers J, Wreden C, Hicks KG, Foster JK, Parthasarathy R, Remington SJ, Guillemin K. Structure and proposed mechanism for the pH-sensing Helicobacter pylori chemoreceptor TlpB. Structure. 2012;20(7):1177–88.
This work was supported by the Okinawa Prefectural Government. Molecular graphics and analyses performed with UCSF Chimera, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311.
Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan (15H02657, 16H06279, 16H05191 and 18KK0266) (Yamaoka Y) and 17K08834 (Suzuki R) and National Institutes of Health Grants DK62813 (Yamaoka Y).
Ethics approval and consent to participate
The animals used in this study were cared for in accordance with our institutional guidelines. Gerbils had free access to food and drinking water throughout the experiment. The experimental protocol was approved by the Animal Care Committee of the Michael E. DeBakey Veterans Affairs Medical Center, Houston, Texas.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Suzuki, R., Satou, K., Shiroma, A. et al. Genome-wide mutation analysis of Helicobacter pylori after inoculation to Mongolian gerbils. Gut Pathog 11, 45 (2019). https://doi.org/10.1186/s13099-019-0326-5
- Helicobacter pylori
- Mongolian gerbil
- Animal model
- Genome comparison
- Adaptive mutation
- Protein structure