- Short report
- Open Access
Phylogenetic analysis of Helicobacter pylori cagA gene of Turkish isolates and the association with gastric pathology
Gut Pathogensvolume 5, Article number: 33 (2013)
The cagA gene is one of the important virulence factors of Helicobacter pylori. The diversity of cagA 5′ conserved region is thought to reflect the phylogenetic relationships between different H. pylori isolates and their association with peptic ulceration. Significant geographical differences among isolates have been reported. The aim of this study is to compare Turkish H. pylori isolates with isolates from different geographical locations and to correlate the association with peptic ulceration.
Total of 52 isolates of which 19 were Turkish and 33 from other geographic locations were studied. Gastric antral biopsies collected from 19 Turkish patients (Gastritis = 12, ulcer = 7) were used to amplify the cagA 5′ region by PCR then followed by DNA sequencing.
The phylogenetic tree displayed 3 groups: A) a mix of 2 sub-groups “Asian” and “African/Anatolian/Asian/European”, B) “Anatolian/European” and C) “American-Indian”. Turkish H. pylori isolates clustered in the mixed sub-group A were mostly from gastritis patients while those clustered in group B were from peptic ulcer patients. A phylogenetic tree constructed for our Turkish isolates detected distinctive features among those from gastritis and ulcer patients. We have found that 2/3 of the gastritis isolates were clustered alone while 1/3 was clustered together with the ulcer isolates. Several amino acids were found to be shared between the later groups but not with the first group of gastritis.
This study provided an additional insight into the profile of our cagA gene which implies a relationship in geographic locations of the isolates.
Helicobacter pylori possess several virulence factors that play important role in the development of gastric diseases. The cag pathogenicity island (cagPAI) a 40-kb locus contains 31 genes among which the cytotoxin-associated gene (cagA) was found to exert a severe damaging effect [1, 2]. Several studies showed that cagA-positive H. pylori strains are more often isolated from patients with gastric ulcers (GU), duodenal ulcers (DU) and gastric cancer (GC) than those with gastritis (G) [3, 4]. Also studies have revealed that the prevalence of cagA-positive strains varies among different populations. In East Asia, nearly all isolated strains were cagA-positive whereas in Western countries this frequency is much lower [5, 6].
The structure of the cagA gene reveals a 5′ highly conserved region and a 3′ variable region. The variation in the size of the CagA protein has been correlated with the varying number of repeat sequences located in the 3′ variable region of the gene that encodes the EPIYA motifs . Based on the types of these motifs the CagA proteins were divided into Western and East Asian types. In a previous study it has been shown that there is less than 53% homology between the CagA repeat sequences of Western and East Asian strains indicating the existence of gene variability that might affect the toxin strength .
H. pylori genetic diversity appears to be a reflection of evolution through thousands of years even before migration out of Africa that lead to geographic spread throughout the world. cagA gene diversity among strains from different geographic regions has been analyzed using different approaches. In one approach EPIYA motifs were used to determine sequence diversity while in others phylogenetic analysis were undertaken on larger portions of the cagA gene . Kawai et al.  recently compared the genome sequence of 4 Japanese H. pylori strains to the available complete genome sequences of Korean, Amerind, European, and West African strains. They reported differences in gene sequences of virulent factors of the Japanese and Korean strains from the others. Camorlinga-Ponce et al.  did phylogenetic analyses of the cagA 3′ region of isolates from various geographic locations and reported that although most H. pylori isolates from indigenous communities of Mexico were of the Western type, a new Amerindian cluster neither Western nor Asian was found. Mane et al.  also studied the cagA gene of a strain isolated from an Amerindian subject and reported that there was a substantial divergence of that Amerindian strain from the Old World strains. Recently Hirai et al.  studied the cagA 3′ region of isolates from asymptomatic healthy Japanese and Thai subjects and compared with those from patients with GC. They suggested that cagA sequence differences between those subjects might be linked to GC incidence in Japan and Thailand. Another study showed that the origin of H. pylori isolates from GC patients were different from those with duodenal ulcer in Japan based on a full length cagA sequences analysis . It has been also shown that differences in the prevalence or expression of cagA gene play an important role in the etiology of these diseases . It is also thought that the diversity of cagA 5′ conserved region reflects the phylogenetic relationships among H. pylori isolates and their association with the clinical outcome [9–11]. Therefore a comparative analysis of H. pylori cagA gene among different isolates may shed light on these concerns. In addition such a study dealing with such analysis of cagA gene is still lacking in Turkey. In this study the diversity of H. pylori cagA gene of isolates from Turkish dyspeptic patients was investigated using DNA sequence analysis and a correlation with the clinical outcome was established. We have also conducted a comparative phylogenetic analysis between the cagA 5′ conserved region of our isolates and isolates from different geographic locations of the world to determine the clustering of these isolates within different regions.
Results and discussion
H. pylori colonies were identified by urease test, gram stain and PCR. A single colony was picked, subcultured and then genomic DNA was isolated. Amplification of the cagA 5′ region revealed a 350 bp band on agarose gel. DNA sequencing of the amplified product was done and the amino acid sequence was deduced. Our constructed phylogenetic tree (Figure 1a) made it possible to distinguish 3 general groups; A, B and C. Group A is a mixed group that contained isolates from nearly all geographic regions and it composed of 2 sub-groups one designated as “Asian” as it contained only isolates from Asia and the other one as “African/Anatolian/Asian/European” as it contained isolates from all those regions. Group B contained isolates from the Anatolian and European region only including the 2 reference strains STR_26695 (gi|306479659|) and STR_J99 (gi|4155035|). Group C is positioned like an out-group with 4 isolates of American-Indian origin. We have found that our Turkish H. pylori isolates from the Anatolian region were clustered in 2 locations on the tree as being part of both group A and B. Most of our isolates from patients with gastritis were clustered in the sub-group A (African/Anatolian/Asian/European) but still others were gathered in sub-group B (Anatolian/European). However isolates from DU and GU patients with one exception (strain 216DU) were clustered only in group B (Anatolian/European). Figure 1b delineate those groups in a more visual solid matrix as an unrooted form of the constructed tree. In this figure the bootstrap values represents the reliability of the tree topology.
We have also constructed a phylogenetic tree only for our Turkish isolates using the same settings as for the whole tree to detect any distinctive features among those from gastritis and ulcer patients (Figure 2). The results showed that isolates from those patients were divergent. We have found that 2/3 of our isolates from gastritis patients were clustered alone while the other 1/3 was clustered together with those isolates from peptic ulcer patients. One interesting finding was the detection of some amino acids that was shared between the later 2 groups (1/3 gastritis and peptic ulcer) but not with the first group of gastritis (longitudinal boxes). Such finding can be further investigated using large number of isolates with known gastric pathology. We also found that some amino acids were unique to the ulcer isolates but not found in the gastritis isolates. Figure 3 depicts a unified origin for the majority of isolates studied that were divergent from the Amerindian isolates which was clearly shown in the alignment set up.
H. pylori possess an enormous genomic diversity that facilitates host adaptation [16, 17]. The high prevalence of such diversity exerted an impact on the clinical outcome. Thus genetic typing of H. pylori isolates can shed the light on the severity of the disease. Since the cagA gene plays a very important role in such diseases, it was the most extensively studied. Comparative sequence analysis of cagA gene of isolates from patients with different gastric pathology, ethnic groups and geographical locations seems to be a method of choice in this regard. Furthermore no such analysis has been reported in Turkey yet. In this study we attempted to assess the genetic relationship between our isolates and those from other locations based on the cagA 5′ region sequences.
The sequence diversities of the cagA gene have been used to deduce the genetic relationships among isolates from Europe, USA, Asia and Africa [18, 19]. Albert et al.  attempted to analyze the cagA genes of 9 isolates from ethnic Arab Kuwaiti patients and reported that isolates were closely related to the Indo-European group but distinct from East Asian strains. Cortes et al.  studied the cagA gene of isolates from 19 Philippinos and reported that the majority of those isolates were of the Western type suggesting Western influence. Rahman et al.  also did phylogenetic analysis using a 219 bp fragment near the 5′ end of cagA gene and showed that Bangladeshi isolates were more closely related to the Indian and Western isolates than those from Chinese and Japanese (East Asian) isolates. van der Ende et al.  also showed earlier that H. pylori cagA-positive isolates from China and Netherlands were distinct.
Based on the assumption of monophyletic origin of cagA sequences we can infer that our H. pylori isolates from the Anatolian region were distributed into two main groups as shown in our constructed phylogenetic tree the mixed group of African/Anatolian/Asian/European and the Anatolian/European sub-groups. The fact that Anatolian location connects the three continents Asia, Africa and Europe made it possible to suggest that H. pylori strains in those regions have undergone genetic modifications through migration events within those continents. Actually this has been shown on our phylogenetic tree by spotting location of Anatolian origin sequence clusters as one of those clusters was closely related to a distinct European group while the other was embedded into the mixed group which is close to the Asian cluster.
Although there is no particular cagA conserved region that could differentiate our Turkish isolates from the other isolates, our constructed phylogenetic tree showed that they were clustered based on the clinical outcome. We have found that our Turkish isolates from gastritis patients were clustered into 2 groups one that included the majority of the isolates and the other included the rest of the gastritis isolates together with isolates from peptic ulcer patients. Isolates from the latter groups were found to share several amino acids that were not found in the first group of gastritis alone. This might suggest that the progression of the disease from gastritis to peptic ulceration might be the result of frequent mutations that occur during the years of infection. The detection of conserved domains specific to those groups could be the reflection of association of those isolates with the clinical disease. In addition some amino acids were found only in peptic ulcer isolates but not in the gastritis isolates which might indicates that such isolates undergo further mutations. Nevertheless further studies comprising additional 5′ region sequences of cagA including isolates from cancer patients can provide more evidence on such correlation. In this study we have shown that our isolates were closely related to the European isolates although some have a tendency to form a separate cluster close to the mixed group. This shows that our isolates have a mixed gene pool while the others were much more scattered. Since the cagA gene undergoes selective pressure over the years a study that utilizes house-keeping genes will further confirm these results.
This study provided an additional insight into the profile of our cagA gene which implies a relationship in geographic locations of the isolates.
Gastric antral biopsies collected from 19 Turkish patients (G = 12, DU = 5, GU = 2) were homogenized, inoculated onto Columbia agar plates with 5% horse blood and incubated at 37°C for 4–5 days under microaerophilic conditions.
Genomic DNA extraction was done using the QIAamp DNA Mini Kit (QIAGEN, Germany) according to the manufacturer′s instructions. Amplification of the cagA 5′ conserved region was done using the primers F1 (5′-GATAACAGGCAAGCTTTTTGAGG-3′) and B1 (5′-CTGCAAAAG ATTGTTTGGCAG-3′) . The amplification steps were done under the following conditions: 94°C for 1 min; 35 cycles of 94°C for 1 min and 57°C for 1 min, and 68°C for 1 min. Final extension was done at 68°C for 5 min. The mixture was stored at 4°C. PCR products were separated by 2% agarose gel electrophoresis and examined under UV illumination. The PCR products were then purified using QIAquick PCR purification kit (QIAGEN, Germany) according to the manufacturer’s instructions.
DNA sequencing was performed by the BigDye Terminator v.1.1 Cycle Sequencing Kit in our lab using ABI PRISM 310 Genetic Analyzer (Applied Biosystems, USA).
The analysis included cagA 5′ region amino acid sequences of 52 isolates of which 19 were Turkish isolates and 33 from other countries including 2 reference strain sequences retrieved from GenBank. To clarify the phylogenetic relationship between our Turkish H. pylori isolates and those from other countries the cagA 5′ region was sequenced and then translated into amino acid sequence using the Open Reading Frame (ORF) Finder program. All possible transcripts of sequences were checked with BLAST program to determine the correct amino acids sequences. The translated 5′ region sequences of Turkish isolates were then aligned with other representative sequences from European, African, Asian and American-Indian isolates using multiple alignments with MUSCLE program . Selection of the conserved blocks and removal of ambiguously aligned divergent blocks from alignment was done with Gblocks 0.91b software [25, 26]. The alignment result was then subjected to PhyML 3.0 program  hosted on Web Server “Phylogeny.fr”  to build phylogenetic tree using Maximum likelihood algorithm . Since the approximate likelihood ratio test (aLRT) and Shimodaira-Hasegawa (SH) provide a compelling alternative to conventional methods , phylogenetic tree was constructed using the aLRT/SH-like. The web based program Boxshade  was used for highlighting sequences according to the similarities. The non-Turkish sequences and their initials with GeneBank accession numbers used in the tree construction were the followings:
hpEurope_1:gi|306480322|,hpEurope_2:gi|306480290|,hpEurope_3:gi|306479848|,hpEurope_4:gi|306479659|,hpEurope_5:gi|306479690|,hpEurope_6:gi|394913|,hpEurope_Swe:gi|37811871|,hpEastAsia_hspEAsia_1:gi|306480475|,hpEastAsia_hspEAsia_2: gi|306480413|,hpEastAsia_hspEAsia_3:gi|306480260|,hpEastAsia_hspEAsia_4:gi|306480140|,hpEastAsia_hspEAsia_5:gi|306479818|,hpEastAsia_hspMaori_1:gi|306480200|,hpEastAsia_hspMaori_2:gi|306479989|,hpEastAsia_hspMaori_3:gi|306479960|,hpEastAsia_hspAmerind_1:gi|306480508|,hpEastAsia_hspAmerind_2:gi|306480494|,hpEastAsia_hspAmerind_3:gi|306479932|,hpEastAsia_hspAmerind_4:gi|306479916|,hpEastAsia_hspAmerind_5:gi|306479895|,hpAfrica1_hspWAfrica_1:gi|306480170|,hpAfrica1_hspWAfrica_2:gi|306480230|,hpAfrica1_hspWAfrica_3:gi|306479787|,hpAfrica1_hspSAfrica_2:gi|306479721|,hpNEAfrica_1:gi|306480445|,hpSahul:gi|306480354|,hpAsiaII_6: gi|306480019|,hpAsiaII_4:gi|306480049|,hpAsiaII_3:gi|306480079|,hpAsiaII_1:gi|306480384|,hpAsiaII_5:gi|306479878|,STR_26695:gi|306479659|, STR_J99:gi|4155035|.
Censini S, Lange C, Xiang Z, Crabtree JE, Ghiara P, Borodovsky M, Rappuoli R: Covacci A: cag, a pathogenicity island of Helicobacter pylori, encodes type I-specific and disease associated virulence factors. Proc Natl Acad Sci U S A. 1996, 93: 14648-14653. 10.1073/pnas.93.25.14648.
Hatakeyama M: Oncogenic mechanism of Helicobacter pylori. Nihon Rinsho Meneki Gakkai Kaishi. 2008, 3: 132-140.
Hatakeyama M: Helicobacter pylori CagA-a bacterial intruder conspiring gastric carcinogenesis. Int J Cancer. 2006, 119: 1217-1223. 10.1002/ijc.21831.
Con SA, Valerin AL, Takeuchi H, Con-Wong R, Con-Chin VG, Con-Chin GR, Yagi-Chaves SN, Mena F, Brenes Pino F, Echandi G, Kobayashi M, Monge-Izaguirre M, Nishioka M, Morimoto N, Sugiura T, Araki K: Helicobacter pylori CagA status associated with gastric cancer incidence rate variability in Costa Rican regions. J Gastroenterol. 2006, 41: 632-637. 10.1007/s00535-006-1812-3.
Pan ZJ, van der Hulst RW, Feller M, Xiao SD, Tytgat GN, Dankert J, van der Ende A: Equally high prevalence of infection with cagA positive Helicobacter pylori in Chinese patients with peptic ulcer disease and those with chronic gastritis associated dyspepsia. J Clin Microbiol. 1997, 35: 1344-1347.
van der Ende A, Pan ZJ, Bart A, van der Hulst RW, Feller M, Xiao SD, Tytgat GN, Dankert J: cagA positive Helicobacter pylori population in China and Netherlands are distinct. Infect Immun. 1998, 66: 1822-1826.
Covacci A, Censini S, Bugnoli M: Molecular characterization of the 128-kDa immunodominant antigen of Helicobacter pylori associated with cytotoxicity and duodenal ulcer. Proc Natl Acad Sci U S A. 1993, 90: 5791-5795. 10.1073/pnas.90.12.5791.
Yamazaki S, Yamakawa A, Okuda T, Ohtani M, Suto H, Ito Y, Yamazaki Y, Keida Y, Higashi H, Hatakeyama M, Azuma T: Distinct diversity of vacA, cagA, and cagE genes of Helicobacter pylori associated with peptic ulcer in Japan. J Clin Microbiol. 2005, 8: 3906-3916.
Duncan SS, Valk PL, Shaffer CL, Bordenstein SR, Cover TL: J-Western forms of Helicobacter pylori cagA constitute a distinct phylogenetic group with a widespread geographic distribution. J Bacteriol. 2012, 194: 1593-1604. 10.1128/JB.06340-11.
Kawai M, Furuta Y, Yahara K, Tsuru T, Oshima K, Handa N, Takahashi N, Yoshida M, Azuma T, Hattori M, Uchiyama I, Kobayashi I: Evolution in an oncogenic bacterial species with extreme genome plasticity: Helicobacter pylori East Asian genomes. BMC Microbiol. 2011, 16: 11-104.
Camorlinga-Ponce M, Perez-Perez G, Gonzalez-Valencia G, Mendoza I, Peñaloza-Espinosa R, Ramos I, Kersulyte D, Reyes-Leon A, Romo C, Granados J, Muñoz L, Berg DE, Torres J: Helicobacter pylori genotyping from American indigenous groups shows novel Amerindian vacA and cagA alleles and Asian. African and European admixture. PLoS One. 2011, 6: e27212-10.1371/journal.pone.0027212.
Mane SP, Dominguez-Bello MG, Blaser MJ, Sobral BW, Hontecillas R, Skoneczka J, Mohapatra SK, Crasta OR, Evans C, Modise T, Shallom S, Shukla M, Varon C, Mégraud F, Maldonado-Contreras AL, Williams KP, Bassaganya-Riera J: Host-interactive genes in Amerindian Helicobacter pylori diverge from their old world homologs and mediate inflammatory responses. J Bacteriol. 2010, 192: 3078-3092. 10.1128/JB.00063-10.
Hirai I, Yoshinaga A, Kimoto A, Sasaki T, Yamamoto Y: Sequence analysis of East Asian cagA of Helicobacter pylori isolated from asymptomatic healthy Japanese and Thai individuals. Curr Microbiol. 2011, 62: 855-860. 10.1007/s00284-010-9797-9.
Satomi S, Yamakawa A, Matsunaga S, Masaki R, Inagaki T, Okuda T, Suto H, Ito Y, Yamazaki Y, Kuriyama M, Keida Y, Kutsumi H, Azuma T: Relationship between the diversity of the cagA gene of Helicobacter pylori and gastric cancer in Okinawa, Japan. J Gastroenterol. 2006, 7: 668-673.
Con SA, Takeuchi H, Valerín AL, Con-Wong R, Con-Chin GR, Con-Chin VG, Nishioka M, Mena F, Brenes F, Yasuda N, Araki K, Sugiura T: Diversity of Helicobacter pylori cagA and vacA genes in Costa Rica: its relationship with atrophic gastritis and gastric cancer. Helicobacter. 2007, 5: 547-552.
Falush D, Wirth T, Linz B, Pritchard JK, Stephens M, Kidd M, Blaser MJ, Graham DY, Vacher S, Perez-Perez GI, Yamaoka Y, Megraud F, Otto K, Reichard U, Katzowitsch E, Wang X, Achtman M, Suerbaum S: Traces of human migrations in Helicobacter pylori populations. Science. 2003, 299: 1582-1585. 10.1126/science.1080857.
Linz B, Balloux F, Moodley Y, Manica A, Liu H, Roumagnac P, Falush D, Stamer C, Prugnolle F, Van Der Merwe SW, Yamaoka Y, Graham DY, Perez-Trallero E, Wadstrom T, Suerbaum S, Achtman M: An African origin for the intimate association between humans and Helicobacter pylori. Nature. 2007, 445: 915-918. 10.1038/nature05562.
Selbach M, Moese S, Hauck CR, Meyer TF, Backert S: Src is the kinase of the Helicobacter pylori CagA protein in vitro and in vivo. J Biol Chem. 2002, 277: 6775-6778. 10.1074/jbc.C100754200.
Higashi H, Tsutsumi R, Muto S, Sugiyama T, Azuma T, Asaka M, Hatakeyama M: SHP-2 tyrosine phosphatase as an intracellular target of Helicobacter pylori CagA protein. Science. 2002, 295: 683-686. 10.1126/science.1067147.
Albert MJ, Al-Akbal HM, Dhar R, De R, Mukhopadhyay AK: Genetic affinities of Helicobacter pylori isolates from ethnic Arabs in Kuwait. Gut Pathog. 2010, 2: 6-10.1186/1757-4749-2-6.
Cortes MC, Yamakawa A, Casingal CR, Fajardo LS, Juan ML, De Guzman BB, Bondoc EM, Mahachai V, Yamazaki Y, Yoshida M, Kutsumi H, Natividad FF, Azuma T, St. Luke's Helicobacter pylori Study Group: Diversity of the cagA gene of Helicobacter pylori strains from patients with gastroduodenal diseases in the Philippines. FEMS Immunol Med Microbiol. 2010, 60: 90-97. 10.1111/j.1574-695X.2010.00722.x.
Rahman M, Mukhopadhyay AK, Nahar S, Datta S, Ahmad MM, Sarker S, Masud IM, Engstrand L, Albert MJ, Nair GB, Berg DE: DNA-level characterization of Helicobacter pylori strains from patients with overt disease and with benign infections in Bangladesh. J Clin Microbiol. 2003, 41: 2008-2014. 10.1128/JCM.41.5.2008-2014.2003.
Tummuru MK, Cover TL, Blaser MJ: Cloning and Expression of a high-molecular-mass major antigen of Helicobacter pylori: evidence of linkage to cytotoxin production. Infect Immun. 1993, 61: 1799-1809.
Edgar RC: MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17: 540-552. 10.1093/oxfordjournals.molbev.a026334.
Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Sys Biol. 2007, 56: 564-577. 10.1080/10635150701472164.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, Claverie JM, Gascuel O: Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36: 465-469. 10.1093/nar/gkn180.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Anisimova M, Gascuel O: Approximate likelihood ratio test for branches: A fast, accurate and powerful alternative. Syst Biol. 2006, 55: 539-552. 10.1080/10635150600755453.
Anisimova M, Gil M, Dufayard JF, Dessimoz C, Gascuel O: Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. Syst Biol. 2011, 60: 685-699. 10.1093/sysbio/syr041.
Pretty printing and shading of multiple-alignment files.http://www.ch.embnet.org/software/BOX_form.html,
This study was supported by the grants of Fatih University (No: P50030903_2) and of the Scientific and Technological Research Council of Turkey (TÜBİTAK) (No: 111 T370).
The authors declare that they have no competing interests.
BAS developed the idea, designed the methodology, analyzed the results and drafted the manuscript. SA carried out endoscopies, biopsy collection and patient evaluation. BKB grew the organism, PCR amplification and DNA sequencing of the gene. MTY carried out detailed analysis of the phylogenetic tree. All authors read and approved the final manuscript.