A Novel Approach to Helicobacter pylori Pan-Genome Analysis for Identification of Genomic Islands

High-quality pan-genome of Escherichia coli generated by excluding confounding and highly similar strains reveals an association between unique gene clusters and genomic islands

The pan-genome analysis of bacteria provides detailed insight into the diversity and evolution of a bacterial population. However, the genomes involved in the pan-genome analysis should be checked carefully, as the inclusion of confounding strains would have unfavorable effects on the identification of core genes, and the highly similar strains could bias the results of the pan-genome state (open versus closed). In this study, we found that the inclusion of highly similar strains also affects the results of unique genes in pan-genome analysis, which leads to a significant underestimation of the number of unique genes in the pan-genome. Therefore, these strains should be excluded from pan-genome analysis at the early stage of data processing. Currently, tens of thousands of genomes have been sequenced for Escherichia coli, which provides an unprecedented opportunity as well as a challenge for pan-genome analysis of this classical model organism. Using the proposed strategies, a high-quality E. coli pan-genome was obtained, and the unique genes was extracted and analyzed, revealing an association between the unique gene clusters and genomic islands from a pan-genome perspective, which may facilitate the identification of genomic islands.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

A Phylogeny-Informed Proteomics Approach for Species Identification within the Burkholderia cepacia Complex

Ancestral genetic exchange between members of many important bacterial pathogen groups has resulted in phylogenetic relationships better described as networks than as bifurcating trees. In certain cases, these reticulated phylogenies have resulted in phenotypic and molecular overlap that challenges the construction of practical approaches for species identification in the clinical microbiology laboratory. Burkholderia cepacia complex (Bcc), a betaproteobacteria species group responsible for significant morbidity in persons with cystic fibrosis and chronic granulomatous disease, represents one such group where network-structured phylogeny has hampered the development of diagnostic methods for species-level discrimination. Here, we present a phylogeny-informed proteomics approach to facilitate diagnostic classification of pathogen groups with reticulated phylogenies, using Bcc as an example. Starting with a set of more than 800 Bcc and Burkholderia gladioli whole-genome assemblies, we constructed phylogenies with explicit representation of inferred interspecies recombination. Sixteen highly discriminatory peptides were chosen to distinguish B. cepacia, Burkholderia cenocepacia, Burkholderia multivorans, and B. gladioli and multiplexed into a single, rapid liquid chromatography-tandem mass spectrometry multiple reaction monitoring (LC-MS/MS MRM) assay. Testing of a blinded set of isolates containing these four Burkholderia species demonstrated 50/50 correct automatic negative calls (100% accuracy with a 95% confidence interval [CI] of 92.9 to 100%), and 70/70 correct automatic species-level positive identifications (100% accuracy with 95% CI 94.9 to 100%) after accounting for a single initial incorrect identification due to a preanalytic error, correctly identified on retesting. The approach to analysis described here is applicable to other pathogen groups for which development of diagnostic classification methods is complicated by interspecies recombination.

Integrative and Conjugative Elements of Helicobacter pylori Are Hypothetical Virulence Factors Associated With Gastric Cancer

Helicobacter pylori is a bacteria with high genome plasticity that has been associated with diverse gastric pathologies. The genetic diversity of this bacteria has limited the characterization of virulence factors associated with gastric cancer (GC). To identify potentially helpful disease biomarkers, we compared 38 complete genomes and 108 draft genomes of H. pylori isolated worldwide from patients with diverse gastric pathologies and 53 draft genomes of H. pylori isolated from Mexican patients with GC, intestinal metaplasia, gastritis, peptic ulcer, and dyspepsia. H. pylori strains isolated from GC were 3-11 times more likely to harbor any of seven genes encoded within an integrative and conjugative element (ICE) than H. pylori isolated from subjects with other gastric pathologies. We tested the cytopathic effects on AGS cells of selected H. pylori strains with known cytotoxin-associated gene pathogenicity island (cag-PAI) and ICE status (H. pylori strains 29CaP, 29CaCe, 62A9, 7C, 8822, and 26695) and the histopathological damage of H. pylori 29CaP and 62A9 in a mouse model. H. pylori 29CaP, which harbors a complete ICEHptfs3 but lacks cag-PAI, elicited distinctive morphology changes and higher histopathological scores compared with other H. pylori strains carrying cag-PAI and hybrid ICE with incomplete TFSS. The presence of intact segments of ICE regions might be a risk factor to develop GC that needs to be addressed in future studies.

Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori

Epigenetic DNA base methylation plays important roles in gene expression regulation. We here describe a gene expression regulation network consisting of many DNA methyltransferases each frequently changing its target sequence-specificity. Our object Helicobacter pylori, a bacterium responsible for most incidence of stomach cancer, carries a large and variable repertoire of sequence-specific DNA methyltransferases. By creating a dozen of single-gene knockout strains for the methyltransferases, we revealed that they form a network controlling methylome, transcriptome and adaptive phenotype sets. The methyltransferases interact with each other in a hierarchical way, sometimes regulated positively by one methyltransferase but negatively with another. Motility, oxidative stress tolerance and DNA damage repair are likewise regulated by multiple methyltransferases. Their regulation sometimes involves translation start and stop codons suggesting coupling of methylation, transcription and translation. The methyltransferases frequently change their sequence-specificity through gene conversion of their target recognition domain and switch their target sets to remodel the network. The emerging picture of a metamorphosing gene regulation network, or firework, consisting of epigenetic systems ever-changing their specificity in search for adaptation, provides a new paradigm in understanding global gene regulation and adaptive evolution.

Analysis of genetic recombination and the pan-genome of a highly recombinogenic bacteriophage species

Bacteriophages are the most prevalent biological entities impacting on the ecosystem and are characterized by their extensive diversity. However, there are two aspects of phages that have remained largely unexplored: genetic flux by recombination between phage populations and characterization of specific phages in terms of the pan-genome. Here, we examined the recombination and pan-genome in Helicobacter pylori prophages at both the genome and gene level. In the genome-level analysis, we applied, for the first time, chromosome painting and fineSTRUCTURE algorithms to a phage species, and showed novel trends in inter-population genetic flux. Notably, hpEastAsia is a phage population that imported a higher proportion of DNA fragments from other phages, whereas the hpSWEurope phages showed weaker signatures of inter-population recombination, suggesting genetic isolation. The gene-level analysis showed that, after parameter tuning of the prokaryote pan-genome analysis program, H. pylori phages have a pan-genome consisting of 75 genes and a soft-core genome of 10 genes, which includes genes involved in the lytic and lysogenic life cycles. Quantitative analysis of recombination events of the soft-core genes showed no substantial variation in the intensity of recombination across the genes, but rather equally frequent recombination among housekeeping genes that were previously reported to be less prone to recombination. The signature of frequent recombination appears to reflect the host–phage evolutionary arms race, either by contributing to escape from bacterial immunity or by protecting the host by producing defective phages.

Prediction of new vaccine targets in the core genome of Corynebacterium pseudotuberculosis through omics approaches and reverse vaccinology

Corynebacterium pseudotuberculosis is the etiologic agent of veterinary relevance diseases, such as caseous lymphadenitis, affecting different animal species causing damage to the global agribusiness. So far, there are no completely effective treatment methods to overcome the impacts caused by this pathogen. Several genomes of the species are deposited on public databases, allowing the execution of studies related to the pan-genomic approach. In this study, we used an integrated in silico workflow to prospect novel putative targets using the core genome, a set of shared genes among 65 C. pseudotuberculosis strains. Subsequently, through RNA-Seq data of the same abiotic stresses in two strains, we selected only induced genes to compose the reverse vaccinology workflow based in two different strategies. Our results predicted six probable antigens in both analysis, which indicates that they have a strong potential to be used in further studies as vaccine targets against this bacterium.

Genome Mining and Comparative Analysis of Streptococcus intermedius Causing Brain Abscess in a Child

Streptococcus intermedius (SI) is associated with prolonged hospitalization and low survival rates. The genetic mechanisms involved in brain abscess development and genome evolution in comparison to other members of the Streptococcus anginosus group are understudied. We performed a whole-genome comparative analysis of an SI isolate, LAU_SINT, associated with brain abscess following sinusitis with all SI genomes in addition to S. constellatus and S. anginosus. Selective pressure on virulence factors, phages, pan-genome evolution and single-nucleotide polymorphism analysis were assessed. The structural details of the type seven secretion system (T7SS) was elucidated and compared with different organisms. ily and nanA were both abundant and conserved. Nisin resistance determinants were found in 47% of the isolates. Pan-genome and SNPs-based analysis didn’t reveal significant geo-patterns. Our results showed that two SC isolates were misidentified as SI. We propose the presence of four T7SS modules (I–IV) located on various genomic islands. We detected a variety of factors linked to metal ions binding on the GIs carrying T7SS. This is the first detailed report characterizing the T7SS and its link to nisin resistance and metal ions binding in SI. These and yet uncharacterized T7SS transmembrane proteins merit further studies and could represent potential therapeutic targets.

The role of integrating conjugative elements in Helicobacter pylori: a review

The genome of Helicobacter pylori contains many putative genes, including a genetic region known as the Integrating Conjugative Elements of H. pylori type four secretion system (ICEHptfs). This genetic regions were originally termed as "plasticity zones/regions" due to the great genetic diversity between the original two H. pylori whole genome sequences. Upon analysis of additional genome sequences, the regions were reported to be extremely common within the genome of H. pylori. Moreover, these regions were also considered conserved rather than genetically plastic and were believed to act as mobile genetic elements transferred via conjugation. Although ICEHptfs(s) are highly conserved, these regions display great allele diversity, especially on ICEHptfs4, with three different subtypes: ICEHptfs4a, 4b, and 4c. ICEHptfs were also reported to contain a novel type 4 secretion system (T4SS) with both epidemiological and in vitro infection model studies highlighting that this novel T4SS functions primarily as a virulence factor. However, there is currently no information regarding the structure, the genes responsible for forming the T4SS, and the interaction between this T4SS and other virulence genes. Unlike the cag pathogenicity island (PAI), which contains CagA, a gene found to be essential for H. pylori virulence, these novel T4SSs have not yet been reported to contain genes that contribute significant effects to the entire system. This notion prompted the hypothesis that these novel T4SSs may have different mechanisms involving cag PAI.

Genomics of Helicobacter pylori

As Helicobacter pylori infects half the world's population and displays an extensive intraspecies diversity, genomics is a powerful tool to understand evolution and disease, to identify factors that confer higher risk of severe sequelae, and to find new approaches for therapy both among bacterial and host targets. In line with these objectives, this review article summarizes the major findings in Helicobacter genomics in papers published between April 2016 and March 2017.

Fallacy of the Unique Genome: Sequence Diversity within Single Helicobacter pylori Strains

Many bacterial genomes are highly variable but nonetheless are typically published as a single assembled genome. Experiments tracking bacterial genome evolution have not looked at the variation present at a given point in time. Here, we analyzed the mouse-passaged Helicobacter pylori strain SS1 and its parent PMSS1 to assess intra- and intergenomic variability. Using high sequence coverage depth and experimental validation, we detected extensive genome plasticity within these H. pylori isolates, including movement of the transposable element IS607, large and small inversions, multiple single nucleotide polymorphisms, and variation in cagA copy number. The cagA gene was found as 1 to 4 tandem copies located off the cag island in both SS1 and PMSS1; this copy number variation correlated with protein expression. To gain insight into the changes that occurred during mouse adaptation, we also compared SS1 and PMSS1 and observed 46 differences that were distinct from the within-genome variation. The most substantial was an insertion in cagY, which encodes a protein required for a type IV secretion system function. We detected modifications in genes coding for two proteins known to affect mouse colonization, the HpaA neuraminyllactose-binding protein and the FutB α-1,3 lipopolysaccharide (LPS) fucosyltransferase, as well as genes predicted to modulate diverse properties. In sum, our work suggests that data from consensus genome assemblies from single colonies may be misleading by failing to represent the variability present. Furthermore, we show that high-depth genomic sequencing data of a population can be analyzed to gain insight into the normal variation within bacterial strains.IMPORTANCE Although it is well known that many bacterial genomes are highly variable, it is nonetheless traditional to refer to, analyze, and publish "the genome" of a bacterial strain. Variability is usually reduced ("only sequence from a single colony"), ignored ("just publish the consensus"), or placed in the "too-hard" basket ("analysis of raw read data is more robust"). Now that whole-genome sequences are regularly used to assess virulence and track outbreaks, a better understanding of the baseline genomic variation present within single strains is needed. Here, we describe the variability seen in typical working stocks and colonies of pathogen Helicobacter pylori model strains SS1 and PMSS1 as revealed by use of high-coverage mate pair next-generation sequencing (NGS) and confirmed by traditional laboratory techniques. This work demonstrates that reliance on a consensus assembly as "the genome" of a bacterial strain may be misleading.

Use of pan-genome analysis for the identification of lineage-specific genes of Helicobacter pylori

The human bacterial pathogen Helicobacter pylori has a highly variable genome, with significant allelic and sequence diversity between isolates and even within well-characterised strains, hampering comparative genomics of H. pylori In this study, pan-genome analysis has been used to identify lineage-specific genes of H. pylori A total of 346 H. pylori genomes spanning the hpAfrica1, hpAfrica2, hpAsia2, hpEurope, hspAmerind and hspEAsia multilocus sequence typing (MLST) lineages were searched for genes specifically overrepresented or underrepresented in MLST lineages or associated with the cag pathogenicity island. The only genes overrepresented in cag-positive genomes were the cag pathogenicity island genes themselves. In contrast, a total of 125 genes were either overrepresented or underrepresented in one or more MLST lineages. Of these 125 genes, alcohol/aldehyde-reducing enzymes linked with acid resistance and production of toxic aldehydes were found to be overrepresented in African lineages. Conversely, the FecA2 ferric citrate receptor was missing from hspAmerind genomes, but present in all other lineages. This work shows the applicability of pan-genome analysis for identification of lineage-specific genes of H. pylori, facilitating further investigation to allow linkage of differential distribution of genes with disease outcome or virulence of H. pylori.