Comparative Analysis of the Full Genome of Helicobacter pylori Isolate Sahul64 Identifies Genes of High Divergence

Isolates of Helicobacter pylori can be classified phylogeographically. High genetic diversity and rapid microevolution are a hallmark of H. pylori genomes, a phenomenon that is proposed to play a functional role in persistence and colonization of diverse human populations. To provide further genomic evidence in the lineage of H. pylori and to further characterize diverse strains of this pathogen in different human populations, we report the finished genome sequence of Sahul64, an H. pylori strain isolated from an indigenous Australian. Our analysis identified genes that were highly divergent compared to the 38 publically available genomes, which include genes involved in the biosynthesis and modification of lipopolysaccharide, putative prophage genes, restriction modification components, and hypothetical genes. Furthermore, the virulence-associated vacA locus is a pseudogene and the cag pathogenicity island (cagPAI) is not present. However, the genome does contain a gene cluster associated with pathogenicity, including dupA. Our analysis found that with the addition of Sahul64 to the 38 genomes, the core genome content of H. pylori is reduced by approximately 14% (∼170 genes) and the pan-genome has expanded from 2,070 to 2,238 genes. We have identified three putative horizontally acquired regions, including one that is likely to have been acquired from the closely related Helicobacter cetorum prior to speciation. Our results suggest that Sahul64, with the absence of cagPAI, highly divergent cell envelope proteins, and a predicted nontransportable VacA protein, could be more highly adapted to ancient indigenous Australian people but with lower virulence potential compared to other sequenced and cagPAI-positive H. pylori strains.     

Gain and loss of multiple genes during the evolution of Helicobacter pylori

Sequence diversity and gene content distinguish most isolates of Helicobacter pylori. Even greater sequence differences differentiate distinct populations of H. pylori from different continents, but it was not clear whether these populations also differ in gene content. To address this question, we tested 56 globally representative strains of H. pylori and four strains of Helicobacter acinonychis with whole genome microarrays. Of the weighted average of 1,531 genes present in the two sequenced genomes, 25% are absent in at least one strain of H. pylori and 21% were absent or variable in H. acinonychis. We extrapolate that the core genome present in all isolates of H. pylori contains 1,111 genes. Variable genes tend to be small and possess unusual GC content; many of them have probably been imported by horizontal gene transfer. Phylogenetic trees based on the microarray data differ from those based on sequences of seven genes from the core genome. These discrepancies are due to homoplasies resulting from independent gene loss by deletion or recombination in multiple strains, which distort phylogenetic patterns. The patterns of these discrepancies versus population structure allow a reconstruction of the timing of the acquisition of variable genes within this species. Variable genes that are located within the cag pathogenicity island were apparently first acquired en bloc after speciation. In contrast, most other variable genes are of unknown function or encode restriction/modification enzymes, transposases, or outer membrane proteins. These seem to have been acquired prior to speciation of H. pylori and were subsequently lost by convergent evolution within individual strains. Thus, the use of microarrays can reveal patterns of gene gain or loss when examined within a phylogenetic context that is based on sequences of core genes.     

Comparative Genome Analysis of Helicobacter pylori Strains

DNA macroarrays were used to characterize 17 Helicobacter pylori strains isolated in four geographic regions of Russia (Moscow, St. Petersburg, Kazan, and Novosibirsk). Of all genes, 1272 (81%) proved to occur in all strains and to constitute a functional core of the genome, and 293 (18.7%) were strain-specific and greatly varied among the H. pylori strains. Most (71%) of the latter had unknown functions; the remainder included restriction–modification genes (3–9%), transposition genes (2–4%), and genes coding for outer membrane proteins (2–4%). The Russian H. pylori strains did not differ in genome organization or in the number and distribution of strain-specific genes from strains isolated in other countries.     

Strain-specific genes of Helicobacter pylori: genome evolution driven by a novel type IV secretion system and genomic island transfer

The availability of multiple bacterial genome sequences has revealed a surprising extent of variability among strains of the same species. The human gastric pathogen Helicobacter pylori is known as one of the most genetically diverse species. We have compared the genome sequence of the duodenal ulcer strain P12 and six other H. pylori genomes to elucidate the genetic repertoire and genome evolution mechanisms of this species. In agreement with previous findings, we estimate that the core genome comprises about 1200 genes and that H. pylori possesses an open pan-genome. Strain-specific genes are preferentially located at potential genome rearrangement sites or in distinct plasticity zones, suggesting two different mechanisms of genome evolution. The P12 genome contains three plasticity zones, two of which encode type IV secretion systems and have typical features of genomic islands. We demonstrate for the first time that one of these islands is capable of self-excision and horizontal transfer by a conjugative process. We also show that excision is mediated by a protein of the XerD family of tyrosine recombinases. Thus, in addition to its natural transformation competence, conjugative transfer of genomic islands has to be considered as an important source of genetic diversity in H. pylori.     

A 500-year tale of co-evolution,adaptation, and virulence: Helicobacter pylori in the Americas

Helicobacter pylori is a common component of the human stomach microbiota, possibly dating back to the speciation of Homo sapiens. A history of pathogen evolution in allopatry has led to the development of genetically distinct H. pylori subpopulations, associated with different human populations, and more recent admixture among H. pylori subpopulations can provide information about human migrations. However, little is known about the degree to which some H. pylori genes are conserved in the face of admixture, potentially indicating host adaptation, or how virulence genes spread among different populations. We analyzed H. pylori genomes from 14 countries in the Americas, strains from the Iberian Peninsula, and public genomes from Europe, Africa, and Asia, to investigate how admixture varies across different regions and gene families. Whole-genome analyses of 723 H. pylori strains from around the world showed evidence of frequent admixture in the American strains with a complex mosaic of contributions from H. pylori populations originating in the Americas as well as other continents. Despite the complex admixture, distinctive genomic fingerprints were identified for each region, revealing novel American H. pylori subpopulations. A pan-genome Fst analysis showed that variation in virulence genes had the strongest fixation in America, compared with non-American populations, and that much of the variation constituted non-synonymous substitutions in functional domains. Network analyses suggest that these virulence genes have followed unique evolutionary paths in the American populations, spreading into different genetic backgrounds, potentially contributing to the high risk of gastric cancer in the region.Subject terms: Population genetics, Microbial genetics     

Comparative Genomics of a Helicobacter pylori Isolate from a Chinese Yunnan Naxi Ethnic Aborigine Suggests High Genetic Divergence and Phage Insertion

Helicobacter pylori is a common pathogen correlated with several severe digestive diseases. It has been reported that isolates associated with different geographic areas, different diseases and different individuals might have variable genomic features. Here, we describe draft genomic sequences of H. pylori strains YN4-84 and YN1-91 isolated from patients with gastritis from the Naxi and Han populations of Yunnan, China, respectively. The draft sequences were compared to 45 other publically available genomes, and a total of 1059 core genes were identified. Genes involved in restriction modification systems, type four secretion system three (TFS3) and type four secretion system four (TFS4), were identified as highly divergent. Both YN4-84 and YN1-91 harbor intact cag pathogenicity island (cagPAI) and have EPIYA-A/B/D type at the carboxyl terminal of cagA. The vacA gene type is s1m2i1. Another major finding was a 32.5-kb prophage integrated in the YN4-84 genome. The prophage shares most of its genes (30/33) with Helicobacter pylori prophage KHP30. Moreover, a 1,886 bp transposable sequence (IS605) was found in the prophage. Our results imply that the Naxi ethnic minority isolate YN4-84 and Han isolate YN1-91 belong to the hspEAsia subgroup and have diverse genome structure. The genome has been extensively modified in several regions involved in horizontal DNA transfer. The important roles played by phages in the ecology and microevolution of H. pylori were further emphasized. The current data will provide valuable information regarding the H. pylori genome based on historic human migrations and population structure.     

Bioinformatic analysis of ESTs collected by Sanger and pyrosequencing methods for a keystone forest tree species: oak

Background

Helicobacter mustelae causes gastritis, ulcers and gastric cancer in ferrets and other mustelids. H. mustelae remains the only helicobacter other than H. pylori that causes gastric ulceration and cancer in its natural host. To improve understanding of H. mustelae pathogenesis, and the ulcerogenic and carcinogenic potential of helicobacters in general, we sequenced the H. mustelae genome, and identified 425 expressed proteins in the envelope and cytosolic proteome.

Results

The H. mustelae genome lacks orthologs of major H. pylori virulence factors including CagA, VacA, BabA, SabA and OipA. However, it encodes ten autotransporter surface proteins, seven of which were detected in the expressed proteome, and which, except for the Hsr protein, are of unknown function. There are 26 putative outer membrane proteins in H. mustelae, some of which are most similar to the Hof proteins of H. pylori. Although homologs of putative virulence determinants of H. pylori (NapA, plasminogen adhesin, collagenase) and Campylobacter jejuni (CiaB, Peb4a) are present in the H. mustelae genome, it also includes a distinct complement of virulence-related genes including a haemagglutinin/haemolysin protein, and a glycosyl transferase for producing blood group A/B on its lipopolysaccharide. The most highly expressed 264 proteins in the cytosolic proteome included many corresponding proteins from H. pylori, but the rank profile in H. mustelae was distinctive. Of 27 genes shown to be essential for H. pylori colonization of the gerbil, all but three had orthologs in H. mustelae, identifying a shared set of core proteins for gastric persistence.

Conclusions

The determination of the genome sequence and expressed proteome of the ulcerogenic species H mustelae provides a comparative model for H. pylori to investigate bacterial gastric carcinogenesis in mammals, and to suggest ways whereby cag minus H. pylori strains might cause ulceration and cancer. The genome sequence was deposited in EMBL/GenBank/DDBJ under accession number FN555004.     

Helicobacter pylori Lipopolysaccharide Is Synthesized via a Novel Pathway with an Evolutionary Connection to Protein N-Glycosylation

Lipopolysaccharide (LPS) is a major component on the surface of Gram negative bacteria and is composed of lipid A-core and the O antigen polysaccharide. O polysaccharides of the gastric pathogen Helicobacter pylori contain Lewis antigens, mimicking glycan structures produced by human cells. The interaction of Lewis antigens with human dendritic cells induces a modulation of the immune response, contributing to the H. pylori virulence. The amount and position of Lewis antigens in the LPS varies among H. pylori isolates, indicating an adaptation to the host. In contrast to most bacteria, the genes for H. pylori O antigen biosynthesis are spread throughout the chromosome, which likely contributed to the fact that the LPS assembly pathway remained uncharacterized. In this study, two enzymes typically involved in LPS biosynthesis were found encoded in the H. pylori genome; the initiating glycosyltransferase WecA, and the O antigen ligase WaaL. Fluorescence microscopy and analysis of LPS from H. pylori mutants revealed that WecA and WaaL are involved in LPS production. Activity of WecA was additionally demonstrated with complementation experiments in Escherichia coli. WaaL ligase activity was shown in vitro. Analysis of the H. pylori genome failed to detect a flippase typically involved in O antigen synthesis. Instead, we identified a homolog of a flippase involved in protein N-glycosylation in other bacteria, although this pathway is not present in H. pylori. This flippase named Wzk was essential for O antigen display in H. pylori and was able to transport various glycans in E. coli. Whereas the O antigen mutants showed normal swimming motility and injection of the toxin CagA into host cells, the uptake of DNA seemed to be affected. We conclude that H. pylori uses a novel LPS biosynthetic pathway, evolutionarily connected to bacterial protein N-glycosylation.     

The complex methylome of the human gastric pathogen Helicobacter pylori

The genome of Helicobacter pylori is remarkable for its large number of restriction-modification (R-M) systems, and strain-specific diversity in R-M systems has been suggested to limit natural transformation, the major driving force of genetic diversification in H. pylori. We have determined the comprehensive methylomes of two H. pylori strains at single base resolution, using Single Molecule Real-Time (SMRT®) sequencing. For strains 26695 and J99-R3, 17 and 22 methylated sequence motifs were identified, respectively. For most motifs, almost all sites occurring in the genome were detected as methylated. Twelve novel methylation patterns corresponding to nine recognition sequences were detected (26695, 3; J99-R3, 6). Functional inactivation, correction of frameshifts as well as cloning and expression of candidate methyltransferases (MTases) permitted not only the functional characterization of multiple, yet undescribed, MTases, but also revealed novel features of both Type I and Type II R-M systems, including frameshift-mediated changes of sequence specificity and the interaction of one MTase with two alternative specificity subunits resulting in different methylation patterns. The methylomes of these well-characterized H. pylori strains will provide a valuable resource for future studies investigating the role of H. pylori R-M systems in limiting transformation as well as in gene regulation and host interaction.     

Strain-specific genes of Helicobacter pylori: distribution, function and dynamics

Whole-genome clustering of the two available genome sequences of Helicobacter pylori strains 26695 and J99 allows the detection of 110 and 52 strain-specific genes, respectively. This set of strain-specific genes was compared with the sets obtained with other computational approaches of direct genome comparison as well as experimental data from microarray analysis. A considerable number of novel function assignments is possible using database-driven sequence annotation, although the function of the majority of the identified genes remains unknown. Using whole-genome clustering, it is also possible to detect species-specific genes by comparing the two H.pylori strains against the genome sequence of Campylobacter jejuni. It is interesting that the majority of strain-specific genes appear to be species specific. Finally, we introduce a novel approach to gene position analysis by employing measures from directional statistics. We show that although the two strains exhibit differences with respect to strain-specific gene distributions, this is due to the extensive genome rearrangements. If these are taken into account, a common pattern for the genome dynamics of the two Helicobacter strains emerges, suggestive of certain spatial constraints that may act as control mechanisms of gene flux.     

The Helicobacter pylori HpyAXII restriction-modification system limits exogenous DNA uptake by targeting GTAC sites but shows asymmetric conservation of the DNA methyltransferase and restriction endonuclease components

The naturally competent organism Helicobacter pylori encodes a large number of restriction–modification (R–M) systems that consist of a restriction endonuclease and a DNA methyltransferase. R–M systems are not only believed to limit DNA exchange among bacteria but may also have other cellular functions. We report a previously uncharacterized H. pylori type II R–M system, M.HpyAXII/R.HpyAXII. We show that this system targets GTAC sites, which are rare in the H. pylori chromosome but numerous in ribosomal RNA genes. As predicted, this type II R–M system showed attributes of a selfish element. Deletion of the methyltransferase M.HpyAXII is lethal when associated with an active endonuclease R.HpyAXII unless compensated by adaptive mutation or gene amplification. R.HpyAXII effectively restricted both unmethylated plasmid and chromosomal DNA during natural transformation and was predicted to belong to the novel ‘half pipe’ structural family of endonucleases. Analysis of a panel of clinical isolates revealed that R.HpyAXII was functional in a small number of H. pylori strains (18.9%, n = 37), whereas the activity of M.HpyAXII was highly conserved (92%, n = 50), suggesting that GTAC methylation confers a selective advantage to H. pylori. However, M.HpyAXII activity did not enhance H. pylori fitness during stomach colonization of a mouse infection model.     

New implications on genomic adaptation derived from the Helicobacter pylori genome comparison

Background

Helicobacter pylori has a reduced genome and lives in a tough environment for long-term persistence. It evolved with its particular characteristics for biological adaptation. Because several H. pylori genome sequences are available, comparative analysis could help to better understand genomic adaptation of this particular bacterium.

Principal Findings

We analyzed nine H. pylori genomes with emphasis on microevolution from a different perspective. Inversion was an important factor to shape the genome structure. Illegitimate recombination not only led to genomic inversion but also inverted fragment duplication, both of which contributed to the creation of new genes and gene family, and further, homological recombination contributed to events of inversion. Based on the information of genomic rearrangement, the first genome scaffold structure of H. pylori last common ancestor was produced. The core genome consists of 1186 genes, of which 22 genes could particularly adapt to human stomach niche. H. pylori contains high proportion of pseudogenes whose genesis was principally caused by homopolynucleotide (HPN) mutations. Such mutations are reversible and facilitate the control of gene expression through the change of DNA structure. The reversible mutations and a quasi-panmictic feature could allow such genes or gene fragments frequently transferred within or between populations. Hence, pseudogenes could be a reservoir of adaptation materials and the HPN mutations could be favorable to H. pylori adaptation, leading to HPN accumulation on the genomes, which corresponds to a special feature of Helicobacter species: extremely high HPN composition of genome.

Conclusion

Our research demonstrated that both genome content and structure of H. pylori have been highly adapted to its particular life style.     

Development of a Multiple-Locus Variable number of tandem repeat Analysis (MLVA) for Leptospira interrogans and its application to Leptospira interrogans serovar Australis isolates from Far North Queensland, Australia

The human gastric pathogen Helicobacter pylori causes chronic gastritis, peptic ulcer disease, gastric carcinoma, and mucosa-associated lymphoid tissue (MALT) lymphoma. It infects over 50% of the worlds' population, however, only a small subset of infected people experience H. pylori-associated illnesses. Associations with disease-specific factors remain enigmatic years after the genome sequences were deciphered. Infection with strains of Helicobacter pylori that carry the cytotoxin-associated antigen A (cagA) gene is associated with gastric carcinoma. Recent studies revealed mechanisms through which the cagA protein triggers oncopathogenic activities. Other candidate genes such as some members of the so-called plasticity region cluster are also implicated to be associated with carcinoma of stomach. Study of the evolution of polymorphisms and sequence variation in H. pylori populations on a global basis has provided a window into the history of human population migration and co-evolution of this pathogen with its host. Possible symbiotic relationships were debated since the discovery of this pathogen. The debate has been further intensified as some studies have posed the possibility that H. pylori infection may be beneficial in some humans. This assumption is based on increased incidence of gastro-oesophageal reflux disease (GERD), Barrett's oesophagus and adenocarcinoma of the oesophagus following H. pylori eradication in some countries. The contribution of comparative genomics to our understanding of the genome organisation and diversity of H. pylori and its pathophysiological importance to human healthcare is exemplified in this review.     

A stable shuttle vector system for efficient genetic complementation of Helicobacter pylori strains by transformation and conjugation

A versatile plasmid shuttle vector system was constructed, which is useful for genetic complementation of Helicobacter pylori strains or mutants with cloned genes of homologous or heterologous origin. The individual plasmid vectors consist of the minimal essential genetic elements, including an origin of replication for Escherichia coli, a H. pylori-specific replicon originally identified on a small cryptic H. pylori plasmid, an oriT sequence and a multiple cloning site. Shuttle plasmid pHel2 carries a chloramphenicol resistance cassette (cat _GC) and pHel3 contains a kanamycin resistance gene (aphA-3) as the selectable marker; both are functional in E. coli and H. pylori. The shuttle plasmids were introduced into the H. pylori strain P1 by natural transformation. A efficiency of 7.0?×?10^?7 and 4.7?×?10^?7 transformants per viable recipient was achieved with pHel2 and pHel3, respectively, and both vectors showed stable, autonomous replication in H. pylori. An approximately 100-fold higher H. pylori transformation rate was obtained when the shuttle vectors for transformation were isolated from the homologous H. pylori strain, rather than E. coli, indicating that DNA restriction and modification mechanisms play a crucial role in plasmid transformation. Interestingly, both shuttle vectors could also be mobilized efficiently from E. coli into different H.?pylori recipients, with pHel2 showing an efficiency of 2.0?×?10^?5 transconjugants per viable H. pylori P1 recipient. Thus, DNA restriction seems to be strongly reduced or absent during conjugal transfer. The functional complementation of a recA-deficient H. pylori mutant by the cloned H. pylorirecA ⁺ gene, and the expression of the heterologous green fluorescent protein (GFP) in H.?pylori demonstrate the general usefulness of?this system, which will significantly facilitate the molecular analysis of H. pylori virulence factors in the future.     

Variability of gene order in different Helicobacter pylori strains contributes to genome diversity

Considerable genomic microdiversity has been reported previously among Helicobacter pylori isolates. We have constructed genome maps of four unrelated H. pylori strains (NCTC11637, NCTC11639, UA802 and UA861) using pulsed-field gel electrophoresis (PFGE) with Notl and Nrul, hybridization with extracted PFGE DNA fragments and probing with 17 gene probes. These strains of H. pylori were compared with a fifth unrelated H. pylori strain NCTC11638 mapped previously. Considerable diversity in gene arrangement was evident among the five H, pylori maps, and no consistent gene clustering was found. The association of only four genes, katA (catalase gene), vacA (vacuo-lating cytotoxin gene), hpaA (a putative adhesin gene), and pfr (bacterial ferritin gene) were generally conserved within approximately the same 25% of the genome; however, the order of these genes also varied. Our study demonstrates that macrodiversity, i.e. variability in gene order, in addition to microdiversity, is a characteristic of the H. pylori genome.     

Sequence Divergence and Conservation in Genomes of Helicobacter cetorum Strains from a Dolphin and a Whale

Background and Objectives

Strains of Helicobacter cetorum have been cultured from several marine mammals and have been found to be closely related in 16 S rDNA sequence to the human gastric pathogen H. pylori, but their genomes were not characterized further.

Methods

The genomes of H. cetorum strains from a dolphin and a whale were sequenced completely using 454 technology and PCR and capillary sequencing.

Results

These genomes are 1.8 and 1.95 mb in size, some 7–26% larger than H. pylori genomes, and differ markedly from one another in gene content, and sequences and arrangements of shared genes. However, each strain is more related overall to H. pylori and its descendant H. acinonychis than to other known species. These H. cetorum strains lack cag pathogenicity islands, but contain novel alleles of the virulence-associated vacuolating cytotoxin (vacA) gene. Of particular note are (i) an extra triplet of vacA genes with ≤50% protein-level identity to each other in the 5′ two-thirds of the gene needed for host factor interaction; (ii) divergent sets of outer membrane protein genes; (iii) several metabolic genes distinct from those of H. pylori; (iv) genes for an iron-cofactored urease related to those of Helicobacter species from terrestrial carnivores, in addition to genes for a nickel co-factored urease; and (v) members of the slr multigene family, some of which modulate host responses to infection and improve Helicobacter growth with mammalian cells.

Conclusions

Our genome sequence data provide a glimpse into the novelty and great genetic diversity of marine helicobacters. These data should aid further analyses of microbial genome diversity and evolution and infection and disease mechanisms in vast and often fragile ocean ecosystems.     

Comparative Genomic Analysis of East Asian and Non-Asian Helicobacter pylori Strains Identifies Rapidly Evolving Genes

Helicobacter pylori infection is a risk factor for the development of gastric adenocarcinoma, a disease that has a high incidence in East Asia. Genes that are highly divergent in East Asian H. pylori strains compared to non-Asian strains are predicted to encode proteins that differ in functional activity and could represent novel determinants of virulence. To identify such proteins, we undertook a comparative analysis of sixteen H. pylori genomes, selected equally from strains classified as East Asian or non-Asian. As expected, the deduced sequences of two known virulence determinants (CagA and VacA) are highly divergent, with 77% and 87% mean amino acid sequence identities between East Asian and non-Asian groups, respectively. In total, we identified 57 protein sequences that are highly divergent between East Asian and non-Asian strains, but relatively conserved within East Asian strains. The most highly represented functional groups are hypothetical proteins, cell envelope proteins and proteins involved in DNA metabolism. Among the divergent genes with known or predicted functions, population genetic analyses indicate that 86% exhibit evidence of positive selection. McDonald-Kreitman tests further indicate that about one third of these highly divergent genes, including cagA and vacA, are under diversifying selection. We conclude that, similar to cagA and vacA, most of the divergent genes identified in this study evolved under positive selection, and represent candidate factors that may account for the disproportionately high incidence of gastric cancer associated with East Asian H. pylori strains. Moreover, these divergent genes represent robust biomarkers that can be used to differentiate East Asian and non-Asian H. pylori strains.     

Infection of less virulent Helicobacter pylori strains in asymptomatic healthy individuals in Thailand as a potential contributing factor to the Asian enigma

In Thailand, gastric cancer incidence is considerably low despite the high prevalence of Helicobacter pylori infection. We investigated the prevalence of H. pylori infection and the genotypes of cagA by using 179 stool specimens obtained from asymptomatic Thai individuals. In this study, the prevalence of H. pylori infection was 43.6%, and the detection rate of cagA-positive strains was 43.5%. In addition, the proportion of the highly virulent East-Asian type of cagA was 7.2%. These results indicate that the low prevalence of cagA-positive H. pylori strain as well as the low prevalence of East-Asian genotype cagA-positive strains may contribute to the low gastric cancer incidence.