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.     

Lipopolysaccharide structures of Helicobacter pylori genomic strains 26695 and J99, mouse model H. pylori Sydney strain, H. pylori P466 carrying sialyl Lewis X, and H. pylori UA915 expressing Lewis B classification of H. pylori lipopolysaccharides into glycotype families.

This study describes the molecular makeup of the cell-wall lipopolysaccharides (LPSs) (O-chain polysaccharide-->core oligosaccharide-->lipid A) from five Helicobacter pylori strains: H. pylori 26695 and J99, the complete genome sequences of which have been published, the established mouse model Sydney strain (SS1), and the symptomatic strains P466 and UA915. All chemical and serological experiments were performed on the intact LPSs. H. pylori 26695 and SS1 possessed either a low-Mr semi-rough-form LPS carrying mostly a single Ley type-2 blood-group determinant in the O-chain region covalently attached to the core oligosaccharide or a high-Mr smooth-form LPS, as did strain J99, with an elongated partially fucosylated type-2 N-acetyllactosamine (polyLacNAc) O-chain polymer, terminated mainly by a Lex blood-group determinant, connected to the core oligosaccharide. In the midst of semi-rough-form LPS glycoforms, H. pylori 26695 and SS1 also expressed in the O-chain region a difucosylated antigen, alpha-L-Fucp(1-3)-alpha-L-Fucp(1-4)-beta-D-GlcpNAc, and the cancer-cell-related type-1 or type-2 linear B-blood-group antigen, alpha-D-Galp(1-3)-beta-D-Galp(1-3 or 4)-beta-D-GlcpNAc. The LPS of H. pylori strain P466 carried the cancer-associated type-2 sialyl Lex blood-group antigen, and the LPS from strain UA915 expressed a type-1 Leb blood-group unit. These findings should aid investigations that focus on identifying and characterizing genes responsible for LPS biosynthesis in genomic strains 26695 and J99, and in understanding the role of H. pylori LPS in animal model studies. The LPSs from the H. pylori strains studied to date were grouped into specific glycotype families.     

Isolation and characterization of a HpyC1I restriction-modification system in Helicobacter pylori

Using transposon shuttle mutagenesis, we identified six Helicobacter pylori mutants from the NTUH-C1 strain that exhibited decreased adherence and cell elongation. Inverse polymerase chain reaction and DNA sequencing revealed that the same locus was interrupted in these six mutants. Nucleotide and amino acid sequences showed no homologies with H. pylori 26695 and J99 strains. This novel open reading frame contained 1617 base pairs. The amino acid sequence shared 24% identity with a putative nicking enzyme in Bacillus halodurans and 23 and 20% identity with type IIS restriction endonucleases PleI and MlyI, respectively. The purified protein, HpyC1I, showed endonuclease activity with the recognition and cleavage site 5'-CCATC(4/5)-3'. Two open reading frames were located upstream of the gene encoding HpyC1I. Together, HpyC1I and these two putative methyltransferases (M1.HpyC1I and M2.HpyC1I) function as a restriction-modification (R-M) system. The HpyC1I R-M genes were found in 9 of the 15 H. pylori strains tested. When compared with the full genome, significantly lower G + C content of HpyC1I R-M genes implied that these genes might have been acquired by horizontal gene transfer. Plasmid DNA transformation efficiencies and chromosomal DNA digestion assays demonstrated protection from HpyC1I digestion by the R-M system. In conclusion, we have identified a novel R-M system present in approximately 60% of H. pylori strains. Disruption of this R-M system results in cell elongation and susceptibility to HpyC1I digestion.     

应用生物信息学方法筛选幽门螺杆菌疫苗候选抗原

目的：应用生物信息学分析方法筛选幽门螺杆菌新的疫苗候选抗原。方法：从TIGRCMR下载幽门螺杆菌26695和J99株全基因组序列,应用生物信息学SignalP、PredTMBB、LipoP、TMHMM、Phobius、PSORT-B和SubLoc等分析软件,筛选幽门螺杆菌新的外膜蛋白和分泌蛋白疫苗候选抗原。结果：从幽门螺杆菌26695株筛选得到54个编码β-桶型跨膜蛋白、脂蛋白或分泌表达蛋白的疫苗候选蛋白抗原,从幽门螺杆菌J99株得到61个呈现上述表达方式的疫苗候选蛋白抗原;且这2株细菌的疫苗候选蛋白呈现良好的交集状况,即有43个候选疫苗蛋白是相同的。结论：用生物信息学分析方法可以从全基因组范围内快速筛选到保守的分泌或表面暴露的疫苗候选抗原,为疫苗抗原的快速筛选与鉴定奠定了基础。     

Genome history in the symbiotic hybrid Euglena gracilis

Euglena gracilis has a chimeric gene collection in which some genes were inherited from its heterotrophic host and others were acquired from a photoautotrophic endosymbiont during secondary endosymbiosis. The evolutionary reconstruction of such a hybrid genome poses a challenge for standard phylogenetic tools that produce bifurcating trees because genome evolution by endosymbiotic gene transfer is a non tree-like process. We sequenced 2770 ESTs from E. gracilis, of which 841 have homologues in a sample of other eukaryotes. Most of these homologues are found in all of the eukaryotes in our sample, but 117 of them are specific to photoautotrophic eukaryotes. A phylogenetic tree fails to account for this observation but the distribution of homologues and a phylogenetic network clearly show the common origin of E. gracilis from both kinetoplastid and photoautotrophic ancestors.     

Helicobacter Pylori's Plasticity Zones Are Novel Transposable Elements

Background

Genes present in only certain strains of a bacterial species can strongly affect cellular phenotypes and evolutionary potentials. One segment that seemed particularly rich in strain-specific genes was found by comparing the first two sequenced Helicobacter pylori genomes (strains 26695 and J99) and was named a “plasticity zone”.

Principal Findings

We studied the nature and evolution of plasticity zones by sequencing them in five more Helicobacter strains, determining their locations in additional strains, and identifying them in recently released genome sequences. They occurred as discrete units, inserted at numerous chromosomal sites, and were usually flanked by direct repeats of 5′AAGAATG, a sequence generally also present in one copy at unoccupied sites in other strains. This showed that plasticity zones are transposable elements, to be called TnPZs. Each full length TnPZ contained a cluster of type IV protein secretion genes (tfs3), a tyrosine recombinase family gene (“xerT”), and a large (≥2800 codon) orf encoding a protein with helicase and DNA methylase domains, plus additional orfs with no homology to genes of known function. Several TnPZ types were found that differed in gene arrangement or DNA sequence. Our analysis also indicated that the first-identified plasticity zones (in strains 26695 and J99) are complex mosaics of TnPZ remnants, formed by multiple TnPZ insertions, and spontaneous and transposable element mediated deletions. Tests using laboratory-generated deletions showed that TnPZs are not essential for viability, but identified one TnPZ that contributed quantitatively to bacterial growth during mouse infection and another that affected synthesis of proinflammatory cytokines in cell culture.

Conclusions

We propose that plasticity zone genes are contained in conjugative transposons (TnPZs) or remnants of them, that TnPZ insertion is mediated by XerT recombinase, and that some TnPZ genes affect bacterial phenotypes and fitness.     

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.     

Functional analysis of the roles of FliQ and FlhB in flagellar expression in Helicobacter pylori

Expression of the two Helicobacter pylori flagellin proteins FlaA and FlaB is required for full motility and persistent infection of the gastric mucosa. The mechanisms and regulation of the biosynthesis and export of flagella in H. pylori are still poorly understood. Scrutiny of the H. pylori 26695 genome sequence revealed homologues of FliQ and FlhB. The roles of the fliQ and flhB genes in H. pylori were investigated by the construction and characterisation of defined isogenic mutants. The results indicate that these genes are involved in the flagellar expression, adhesion to and colonisation of the gastric mucosa.     

Systematic identification of selective essential genes in Helicobacter pylori by genome prioritization and allelic replacement mutagenesis

A comparative genomic approach was used to identify Helicobacter pylori 26695 open reading frames (ORFs) which are conserved in H. pylori J99 but highly diverged in other eubacteria. A survey of selected pathways of central intermediary metabolism was also carried out, and genes with a potentially selective role in H. pylori were identified. Forty-five ORFs identified in these two analyses were screened using a rapid vector-free allelic replacement mutagenesis technique, and 33 were shown to be essential in vitro. Notably, 13 ORFs gave essentiality results which are unexpected in view of their known or proposed functions, and phylogenetic analysis was used to investigate the annotation of 7 such ORFs which are highly diverged. We propose that the products of a number of these H. pylori-specific essential genes may be suitable targets for novel anti-H. pylori therapies.     

Evolution of the chloroplast genome

We discuss the suggestion that differences in the nucleotide composition between plastid and nuclear genomes may provide a selective advantage in the transposition of genes from plastid to nucleus. We show that in the adenine, thymine (AT)-rich genome of Borrelia burgdorferi several genes have an AT-content lower than the average for the genome as a whole. However, genes whose plant homologues have moved from plastid to nucleus are no less AT-rich than genes whose plant homologues have remained in the plastid, indicating that both classes of gene are able to support a high AT-content. We describe the anomalous organization of dinoflagellate plastid genes. These are located on small circles of 2-3 kbp, in contrast to the usual plastid genome organization of a single large circle of 100-200 kbp. Most circles contain a single gene. Some circles contain two genes and some contain none. Dinoflagellate plastids have retained far fewer genes than other plastids. We discuss a similarity between the dinoflagellate minicircles and the bacterial integron system.     

Role of Helicobacter pylori rfaJ genes (HP0159 and HP1416) in lipopolysaccharide synthesis

The genome of Helicobacter pylori 26695 has been sequenced and the lipopolysaccharide (LPS) O sidechain of this strain has been shown to express both Lewis x and Lewis y units. To determine the role of HP0159 and HP1416, genes recognized as rfaJ homologs and implicated in LPS synthesis, isogenic mutants of H. pylori 26695 were generated. The LPS of mutant 26695::HP0159Kan did not express either Lewis epitope as detected by immunoblotting, whereas the control strain and 26695::HP1416Kan produced both epitopes. Structural analysis of the LPS of the mutants showed that HP0159 encodes an alpha(1,2/3)-glucosyltransferase whereas HP1416 encodes an alpha(1,2/4)-glucosyltransferase.     

Complementation system for Helicobacter pylori

Previously Langford et al. (2006) developed the pIR203C04 complementation system for Helicobacter pylori, which can be used to complement and restore phenotypic effects in H. pylori mutant, and furthermore they used the complementation system in vivo experiments to animals without altering the ability of strain SSI to colonize mice. In their previous study, the pIR203C04 was able to transform 26695, SSI, J99, and 43504 H. pylori strains by an electroporation method. However, in the present study using a natural transformation the pIR203C04 transformed only 26695 H. pylori but not SSI, J99, 7.13, and G27 H. pylori strains. Since the useful complementation system has a limitation of narrow selection among H. pylori strains, we redesigned the complementation system for the improvement. The same intergenic chromosomal site between hp0203 and hp0204 was utilized for the new complementation system because the insertion at the intergenic site didn’t show any polar effects and disruption of other H. pylori genes. The genome sequence analysis showed that the intergenic regions among H. pylori strains may have too low homology to each others to do a homologous recombination. Thus, in addition to the short intergenic region, the fragments of the new complementation system included 3′ conserved parts of hp0203 and hp0204 coding regions. Between the fragments there are a chloramphenicol acetyltransferase cassette and multicloning sites, resulting in pKJMSH. DNA fragment of the interest can be cloned into the multicloning sites of pKJMSH and the fragment can be integrated at the intergenic region of H. pylori chromosome by the homologous recombination. Indeed, by the natural transformation, pKJMSH was able to transform all five H. pylori strains of 26695, SSI, J99, 7.13, and G27, which are common for the investigation of molecular pathogenesis. Thus, the new pKJMSH complementation system is applicable to most H. pylori wild-type stains.     

Comparative proteome analysis of Helicobacter pylori

Helicobacter pylori, the causative agent of gastritis, ulcer and stomach carcinoma, infects approximately half of the worlds population. After sequencing the complete genome of two strains, 26695 and J99, we have approached the demanding task of investigating the functional part of the genetic information containing macromolecules, the proteome. The proteins of three strains of H. pylori, 26695 and J99, and a prominent strain used in animal models SS1, were separated by a high-resolution two-dimensional electrophoresis technique with a resolution power of 5000 protein spots. Up to 1800 protein species were separated from H. pylori which had been cultivated for 5 days on agar plates. Using matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) peptide mass fingerprinting we have identified 152 proteins, including nine known virulence factors and 28 antigens. The three strains investigated had only a few protein spots in common. We observe that proteins with an amino acid exchange resulting in a net change of only one charge are shifted in the two-dimensional electrophoresis (2-DE) pattern. The expression of 27 predicted conserved hypothetical open reading frames (ORFs) and six unknown ORFs were confirmed. The growth conditions of the bacteria were shown to have an effect on the presence of certain proteins. A preliminary immunoblotting study using human sera revealed that this approach is ideal for identifying proteins of diagnostic or therapeutic value. H. pylori 2-DE patterns with their identified protein species were added to the dynamic 2D-PAGE database (http://www.mpiib-berlin.mpg.de/2D-PAGE/). This basic knowledge of the proteome in the public domain will be an effective instrument for the identification of new virulence or pathogenic factors, and antigens of potentially diagnostic or curative value against H. pylori.     

A Variable Gene in a Conserved Region of the Helicobacter pylori Genome: Isotopic Gene Replacement or Rapid Evolution?

The present study concerns the identification of a novel coding sequence in a region of the Helicobacter pylori genome, located between JHP1069/HP1141 and JHP1071/HP1143 according to the numbering of the J99 and 26695 reference strains, respectively, and spanning three different coding DNA sequences (CDSs). The CDSs located at the centre of this locus were highly polymorphic, as determined by the analysis of 24 European isolates, 3 Asian, and 3 African isolates. Phylogenetic and molecular evolutionary analyses showed that the CDSs were not restricted to the geographical origin of the strains. Despite a very high variability observed in the deduced protein sequences, significant similarity was observed, always with the same protein families, i.e. ATPase and bacteriophage receptor/invasion proteins. Although this variability could be explained by isotopic gene replacement via horizontal transfer of a gene with the same function but coming from a variety of sources, it seems more likely that the very high sequence variation observed at this locus is the result of a strong selection pressure exerted on the corresponding gene product. The CDSs identified in the present study could be used as strain specific markers.     

Nucleotide correspondence between protein-coding sequences of Helicobacter pylori 26695 and J99 strains

Comparison of open-reading frames (ORFs) H. pylori 26695 and J99 strains has been revealed prevalence of nucleotide replacements as transitions (more than 3%) above transversions (less than 1%). Prevalence of nucleotide transitions is caused by high speed of C : G to T : A transitions in a coding strand of DNA (3.5-5.3%) and not coding strand (2.9-3.9%). The correspondence rate of transversion (A --> C, A --> T, C --> A, C --> G, G --> C, G --> T, T --> A and T --> G) did not exceed 0.84%. The highest correspondence frequency between C and T was detected in ACGT-ATGT (28.3%) - the site of methylation by active methyltransferase M.Hpy99XI in H. pylori 26695 and J99. Thus one can speculate that predominant transition taking place in H. pylori is mutation of C into T, which is realized through cytosine methylation-deamination mechanism.     

Finding pathogenicity islands and gene transfer events in genome data

MOTIVATION: There is a growing literature on wavelet theory and wavelet methods showing improvements on more classical techniques, especially in the contexts of smoothing and extraction of fundamental components of signals. G+C patterns occur at different lengths (scales) and, for this reason, G+C plots are usually difficult to interpret. Current methods for genome analysis choose a window size and compute a chi(2) statistics of the average value for each window with respect to the whole genome. RESULTS: Firstly, wavelets are used to smooth G+C profiles to locate characteristic patterns in genome sequences. The method we use is based on performing a chi(2) statistics on the wavelet coefficients of a profile; thus we do not need to choose a fixed window size, in that the smoothing occurs at a set of different scales. Secondly, a wavelet scalogram is used as a measure for sequence profile comparison; this tool is very general and can be applied to other sequence profiles commonly used in genome analysis. We show applications to the analysis of Deinococcus radiodurans chromosome I, of two strains of Helicobacter pylori (26695, J99) and two of Neisseria meningitidis (serogroup B strain MC58 and serogroup A strain Z2491). We report a list of loci that have different G+C content with respect to the nearby regions; the analysis of N. meningitidis serogroup B shows two new large regions with low G+C content that are putative pathogenicity islands. AVAILABILITY: Software and numerical results (profiles, scalograms, high and low frequency components) for all the genome sequences analyzed are available upon request from the authors.