Choosing a benchtop sequencing machine to characterise Helicobacter pylori genomes

Genomic analysis of Helicobacter pylori in Australia: Antimicrobial resistance, phylogenetic patterns, and virulence factors

BACKGROUND AND AIM

Rates of antimicrobial-resistant Helicobacter pylori infection are rising globally, but little is known about contemporary resistance patterns, virulence factors, and phylogenetic patterns of isolates within Australia. We aimed to characterize antimicrobial resistance and genetic mutations associated with adverse clinical outcomes.

METHODS

Whole genome sequencing, culturing, and antibiotic sensitivity data for refractory H. pylori isolates at Australian centers were collected between 2013 and 2022. Phylogenetic origins, antibiotic resistance mutations, and virulence factors were examined with phenotypic resistance profiles.

RESULTS

One hundred thirty-five isolates underwent culture, with 109 of these undergoing whole genome sequencing. Forty-three isolates were isolated from patients in South Australia and 66 from Western Australia. Isolates originated primarily from hpEurope (59.6%), hpEastAsia (25.7%), and hpNEAfrica (6.4%). Antimicrobial resistance to clarithromycin was seen in 85% of isolates, metronidazole in 52%, levofloxacin in 18%, rifampicin in 14%, and amoxicillin in 9%. Most isolates (59%) were multi-drug resistant. Resistance concordance between genetically determined resistance and phenotypic resistance was 92% for clarithromycin and 94% for levofloxacin. Analysis of virulence factors demonstrated cag pathogenicity island (cagPAI) in 67% of isolates and cagA in 61%, correlating with isolate genetic origin. The most virulent s1m1 vacuolating cytotoxin A genotype was present in 26% of isolates.

CONCLUSION

Refractory H. pylori isolates in Australia emanate from multiple global origins. Strong concordance between genetic and phenotypic antibiotic resistance profiles raises the possibility of utilizing genetic profiling in clinical practice. The dynamic landscape of H. pylori in Australia warrants the establishment of a national database to monitor H. pylori resistance and evolving virulence.

Cancer-related genetic variants of Helicobacter pylori strains determined using gastric wash-based whole-genome analysis with single-molecule real-time technology

Helicobacter pylori (H. pylori) are a primary factor in the pathogenesis of gastric cancer (GC); GC ranks third among cancer-related mortality. A clear understanding of the H. pylori genome factors underlying GC is necessary to develop more effective methods to prevent GC. A single-molecule real-time DNA sequencing-based H. pylori genome-wide association study analysis was performed using the H. pylori genome present in five early-stage GC (EGC) and five non-GC clinical DNA samples recovered from gastric washes. A total of 275 genes with 702 nucleotide variants (NVs) were found to be common to three or more patients with EGC but no non-GC patients (single-NV: 654/702, 93.2%; multi-NV: 40/702, 5.7%; deletion: 3/702, 0.4%; insertion: 3/702, 0.7%). Gene ontology analysis of H. pylori revealed that genes involved in the mitochondrial electron transport system, glycolytic processes and the TCA cycle were highly enriched. Cancer-related NVs were most frequently found in a member of the Helicobacter outer membrane protein family, hopL. In particular, one of the NVs in hopL was a novel six-nucleotide insertion (1159095̂1159096, TACTTC); this mutant was detected more frequently in a validation set of 50 additional EGC samples (22/50, 44.0%) than in 18 non-GC samples (3/18, 16.7%, P = .04). These results suggest that the hopL variant is associated with the development of GC and may serve as a genetic biomarker of H. pylori virulence and GC risk. Our assay can serve as a potent tool to expand our understanding of bacteria-associated tumorigenesis.

Biomimetic self-assembly of subcellular structures

This article summarizes recent progress on biomimetic subcellular structures and discusses integration of these isolated systems.

Optimizing sequencing protocols for leaderboard metagenomics by combining long and short reads

As metagenomic studies move to increasing numbers of samples, communities like the human gut may benefit more from the assembly of abundant microbes in many samples, rather than the exhaustive assembly of fewer samples. We term this approach leaderboard metagenome sequencing. To explore protocol optimization for leaderboard metagenomics in real samples, we introduce a benchmark of library prep and sequencing using internal references generated by synthetic long-read technology, allowing us to evaluate high-throughput library preparation methods against gold-standard reference genomes derived from the samples themselves. We introduce a low-cost protocol for high-throughput library preparation and sequencing.

Review of current diagnostic methods and advances in Helicobacter pylori diagnostics in the era of next generation sequencing

Helicobacter pylori (H. pylori) infection is highly prevalent in the human population and may lead to severe gastrointestinal pathology including gastric and duodenal ulcers, mucosa associated tissue lymphoma and gastric adenocarcinoma. In recent years, an alarming increase in antimicrobial resistance and subsequently failing empiric H. pylori eradication therapies have been noted worldwide, also in many European countries. Therefore, rapid and accurate determination of H. pylori’s antibiotic susceptibility prior to the administration of eradication regimens becomes ever more important. Traditionally, detection of H. pylori and its antimicrobial resistance is done by culture and phenotypic drug susceptibility testing that are cumbersome with a long turn-around-time. Recent advances in diagnostics provide new tools, like real-time polymerase chain reaction (PCR) and line probe assays, to diagnose H. pylori infection and antimicrobial resistance to certain antibiotics, directly from clinical specimens. Moreover, high-throughput whole genome sequencing technologies allow the rapid analysis of the pathogen’s genome, thereby allowing identification of resistance mutations and associated antibiotic resistance. In the first part of this review, we will give an overview on currently available diagnostic methods for detection of H. pylori and its drug resistance and their implementation in H. pylori management. The second part of the review focusses on the use of next generation sequencing technology in H. pylori research. To this end, we conducted a literature search for original research articles in English using the terms “Helicobacter”, “transcriptomic”, “transcriptome”, “next generation sequencing” and “whole genome sequencing”. This review is aimed to bridge the gap between current diagnostic practice (histology, rapid urease test, H. pylori culture, PCR and line probe assays) and new sequencing technologies and their potential implementation in diagnostic laboratory settings in order to complement the currently recommended H. pylori management guidelines and subsequently improve public health.

Increased sequencing depth does not increase captured diversity of arbuscular mycorrhizal fungi

The arrival of 454 sequencing represented a major breakthrough by allowing deeper sequencing of environmental samples than was possible with existing Sanger approaches. Illumina MiSeq provides a further increase in sequencing depth but shorter read length compared with 454 sequencing. We explored whether Illumina sequencing improves estimates of arbuscular mycorrhizal (AM) fungal richness in plant root samples, compared with 454 sequencing. We identified AM fungi in root samples by sequencing amplicons of the SSU rRNA gene with 454 and Illumina MiSeq paired-end sequencing. In addition, we sequenced metagenomic DNA without prior PCR amplification. Amplicon-based Illumina sequencing yielded two orders of magnitude higher sequencing depth per sample than 454 sequencing. Initial analysis with minimal quality control recorded five times higher AM fungal richness per sample with Illumina sequencing. Additional quality control of Illumina samples, including restriction of the marker region to the most variable amplicon fragment, revealed AM fungal richness values close to those produced by 454 sequencing. Furthermore, AM fungal richness estimates were not correlated with sequencing depth between 300 and 30,000 reads per sample, suggesting that the lower end of this range is sufficient for adequate description of AM fungal communities. By contrast, metagenomic Illumina sequencing yielded very few AM fungal reads and taxa and was dominated by plant DNA, suggesting that AM fungal DNA is present at prohibitively low abundance in colonised root samples. In conclusion, Illumina MiSeq sequencing yielded higher sequencing depth, but similar richness of AM fungi in root samples, compared with 454 sequencing.

Quantum changes in Helicobacter pylori gene expression accompany host-adaptation

Helicobacter pylori is a highly successful gastric pathogen. High genomic plasticity allows its adaptation to changing host environments. Complete genomes of H. pylori clinical isolate UM032 and its mice-adapted serial derivatives 298 and 299, generated using both PacBio RS and Illumina MiSeq sequencing technologies, were compared to identify novel elements responsible for host-adaptation. The acquisition of a jhp0562-like allele, which encodes for a galactosyltransferase, was identified in the mice-adapted strains. Our analysis implies a new β-1,4-galactosyltransferase role for this enzyme, essential for Le^y antigen expression. Intragenomic recombination between babA and babB genes was also observed. Further, we expanded on the list of candidate genes whose expression patterns have been mediated by upstream homopolymer-length alterations to facilitate host adaption. Importantly, greater than four-fold reduction of mRNA levels was demonstrated in five genes. Among the down-regulated genes, three encode for outer membrane proteins, including BabA, BabB and HopD. As expected, a substantial reduction in BabA protein abundance was detected in mice-adapted strains 298 and 299 via Western analysis. Our results suggest that the expression of Le^y antigen and reduced outer membrane protein expressions may facilitate H. pylori colonisation of mouse gastric epithelium.

Comparative Genomics Revealed Multiple Helicobacter pylori Genes Associated with Biofilm Formation In Vitro

Background

Biofilm formation by Helicobacter pylori may be one of the factors influencing eradication outcome. However, genetic differences between good and poor biofilm forming strains have not been studied.

Materials and Methods

Biofilm yield of 32 Helicobacter pylori strains (standard strain and 31 clinical strains) were determined by crystal-violet assay and grouped into poor, moderate and good biofilm forming groups. Whole genome sequencing of these 32 clinical strains was performed on the Illumina MiSeq platform. Annotation and comparison of the differences between the genomic sequences were carried out using RAST (Rapid Annotation using Subsystem Technology) and SEED viewer. Genes identified were confirmed using PCR.

Results

Genes identified to be associated with biofilm formation in H. pylori includes alpha (1,3)-fucosyltransferase, flagellar protein, 3 hypothetical proteins, outer membrane protein and a cag pathogenicity island protein. These genes play a role in bacterial motility, lipopolysaccharide (LPS) synthesis, Lewis antigen synthesis, adhesion and/or the type-IV secretion system (T4SS). Deletion of cagA and cagPAI confirmed that CagA and T4SS were involved in H. pylori biofilm formation.

Conclusions

Results from this study suggest that biofilm formation in H. pylori might be genetically determined and might be influenced by multiple genes. Good, moderate and poor biofilm forming strain might differ during the initiation of biofilm formation.

Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

BACKGROUND

Across species, diversity at the Major Histocompatibility Complex (MHC) is critical to disease resistance and population health; however, use of MHC diversity to quantify the genetic health of populations has been hampered by the extreme variation found in MHC genes. Next generation sequencing (NGS) technology generates sufficient data to genotype even the most diverse species, but workflows for distinguishing artifacts from alleles are still under development. We used NGS to evaluate the MHC diversity of over 300 captive and wild ring-tailed lemurs (Lemur catta: Primates: Mammalia). We modified a published workflow to address errors that arise from deep sequencing individuals and tested for evidence of selection at the most diverse MHC genes.

RESULTS

In addition to evaluating the accuracy of 454 Titanium and Ion Torrent PGM for genotyping large populations at hypervariable genes, we suggested modifications to improve current methods of allele calling. Using these modifications, we genotyped 302 out of 319 individuals, obtaining an average sequencing depth of over 1000 reads per amplicon. We identified 55 MHC-DRB alleles, 51 of which were previously undescribed, and provide the first sequences of five additional MHC genes: DOA, DOB, DPA, DQA, and DRA. The additional five MHC genes had one or two alleles each with little sequence variation; however, the 55 MHC-DRB alleles showed a high dN/dS ratio and trans-species polymorphism, indicating a history of positive selection. Because each individual possessed 1-7 MHC-DRB alleles, we suggest that ring-tailed lemurs have four, putatively functional, MHC-DRB copies.

CONCLUSIONS

In the future, accurate genotyping methods for NGS data will be critical to assessing genetic variation in non-model species. We recommend that future NGS studies increase the proportion of replicated samples, both within and across platforms, particularly for hypervariable genes like the MHC. Quantifying MHC diversity within non-model species is the first step to assessing the relationship of genetic diversity at functional loci to individual fitness and population viability. Owing to MHC-DRB diversity and copy number, ring-tailed lemurs may serve as an ideal model for estimating the interaction between genetic diversity, fitness, and environment, especially regarding endangered species.

Divergence between the Highly Virulent Zoonotic Pathogen Helicobacter heilmannii and Its Closest Relative, the Low-Virulence "Helicobacter ailurogastricus" sp. nov

Helicobacter heilmannii naturally colonizes the stomachs of dogs and cats and has been associated with gastric disorders in humans. Nine feline Helicobacter strains, classified as H. heilmannii based on ureAB and 16S rRNA gene sequences, were divided into a highly virulent and a low-virulence group. The genomes of these strains were sequenced to investigate their phylogenetic relationships, to define their gene content and diversity, and to determine if the differences in pathogenicity were associated with the presence or absence of potential virulence genes. The capacities of these helicobacters to bind to the gastric mucosa were investigated as well. Our analyses revealed that the low-virulence strains do not belong to the species H. heilmannii but to a novel, closely related species for which we propose the name Helicobacter ailurogastricus. Several homologs of H. pylori virulence factors, such as IceA1, HrgA, and jhp0562-like glycosyltransferase, are present in H. heilmannii but absent in H. ailurogastricus. Both species contain a VacA-like autotransporter, for which the passenger domain is remarkably larger in H. ailurogastricus than in H. heilmannii. In addition, H. ailurogastricus shows clear differences in binding to the gastric mucosa compared to H. heilmannii. These findings highlight the low-virulence character of this novel Helicobacter species.

Low-Cost, High-Throughput Sequencing of DNA Assemblies Using a Highly Multiplexed Nextera Process

In recent years, next-generation sequencing (NGS) technology has greatly reduced the cost of sequencing whole genomes, whereas the cost of sequence verification of plasmids via Sanger sequencing has remained high. Consequently, industrial-scale strain engineers either limit the number of designs or take short cuts in quality control. Here, we show that over 4000 plasmids can be completely sequenced in one Illumina MiSeq run for less than $3 each (15× coverage), which is a 20-fold reduction over using Sanger sequencing (2× coverage). We reduced the volume of the Nextera tagmentation reaction by 100-fold and developed an automated workflow to prepare thousands of samples for sequencing. We also developed software to track the samples and associated sequence data and to rapidly identify correctly assembled constructs having the fewest defects. As DNA synthesis and assembly become a centralized commodity, this NGS quality control (QC) process will be essential to groups operating high-throughput pipelines for DNA construction.

The complete methylome of Helicobacter pylori UM032

BACKGROUND

The genome of the human gastric pathogen Helicobacter pylori encodes a large number of DNA methyltransferases (MTases), some of which are shared among many strains, and others of which are unique to a given strain. The MTases have potential roles in the survival of the bacterium. In this study, we sequenced a Malaysian H. pylori clinical strain, designated UM032, by using a combination of PacBio Single Molecule, Real-Time (SMRT) and Illumina MiSeq next generation sequencing platforms, and used the SMRT data to characterize the set of methylated bases (the methylome).

RESULTS

The N4-methylcytosine and N6-methyladenine modifications detected at single-base resolution using SMRT technology revealed 17 methylated sequence motifs corresponding to one Type I and 16 Type II restriction-modification (R-M) systems. Previously unassigned methylation motifs were now assigned to their respective MTases-coding genes. Furthermore, one gene that appears to be inactive in the H. pylori UM032 genome during normal growth was characterized by cloning.

CONCLUSION

Consistent with previously-studied H. pylori strains, we show that strain UM032 contains a relatively large number of R-M systems, including some MTase activities with novel specificities. Additional studies are underway to further elucidating the biological significance of the R-M systems in the physiology and pathogenesis of H. pylori.

Evaluation of convenient pretreatment protocols for RNA virus metagenomics in serum and tissue samples

Viral metagenomic approaches are increasingly being used for viral discovery. Various strategies are applied to enrich viral sequences, but there is often a lack of knowledge about their effective influence on the viral discovery sensitivity. We evaluate some convenient and widely used approaches for RNA virus discovery in clinical samples in order to reveal their sensitivity and potential bias introduced by the enrichment or amplifications steps. An RNA virus was artificially spiked at a fixed titer in serum and lung tissue, respectively, low and high nucleic acid content matrices. For serum, a simple DNase treatment on the RNA extract gave the maximum gain in proportion of viral sequences (83×), and a subsequent ribosomal RNA removal nearly doubled once more the proportion of viral sequences. For lung tissue, a ribosomal RNA depletion step on the RNA extract had the biggest gain in proportion of viral sequences (32×). We show also that direct sequencing of cDNA is recommended above an extra random PCR amplification step, and a that the virion enrichment strategy (filtration and nuclease treatment) has a beneficial effect for sequencing-based virus discovery. Our findings provide sample-dependent guidelines for targeted virus discovery strategies.

Molecular oncology testing in resource-limited settings

Cancer prevalence and mortality are high in developing nations, where resources for cancer control are inadequate. Nearly one-quarter of cancers in resource-limited nations are infection related, and molecular assays can capitalize on this relationship by detecting pertinent pathogen genomes and human gene variants to identify those at highest risk for progression to cancer, to classify lesions, to predict effective therapy, and to monitor tumor burden over time. Prime examples are human papillomavirus in cervical neoplasia, Helicobacter pylori and Epstein-Barr virus in gastric adenocarcinoma and lymphoma, and hepatitis B or C virus in hepatocellular cancer. Research is underway to engineer devices that overcome social, economic, and technical barriers limiting effective laboratory support. Additional challenges include an educated workforce, infrastructure for quality metrics and record keeping, and funds to sustain molecular test services. The combination of well-designed interfaces, novel and robust electrochemical technology, and telemedicine tools will promote adoption by frontline providers. Fast turnaround is crucial for surmounting loss to follow-up, although increased use of cell phones, even in rural areas, enhances options for patient education and engagement. Links to a broadband network facilitate consultation and centralized storage of medical data. Molecular technology shows promise to address gaps in health care through rapid, user-friendly, and cost-effective devices reflecting clinical priorities in resource-poor areas.

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population

BACKGROUND

Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously.

RESULTS

Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%.

CONCLUSIONS

Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

Whole-genome sequencing of clarithromycin resistant Helicobacter pylori characterizes unidentified variants of multidrug resistant efflux pump genes

BACKGROUND

Clarithromycin (CLR) is the key drug in eradication therapy of Helicobacter pylori (H. pylori) infection, and widespread use of CLR has led to an increase in primary CLR-resistant H. pylori. The known mechanism of CLR resistance has been established in A2146G and A2147G mutations in the 23S rRNA gene, but evidence of the involvement of other genetic mechanisms is lacking. Using the MiSeq platform, whole-genome sequencing of the 19 clinical strains and the reference strain ATCC26695 was performed to identify single nucleotide variants (SNVs) of multi-drug resistant efflux pump genes in the CLR-resistant phenotype.

RESULTS

Based on sequencing data of ATCC26695, over one million sequencing reads with over 50-fold coverage were sufficient to detect SNVs, but not indels in the bacterial genome. Sequencing reads of the clinical isolates ranged from 1.82 to 10.8 million, and average coverage ranged from 90.9- to 686.3-fold, which were acceptable criteria for detecting SNVs. Utilizing the conventional approach of allele-specific PCR, point mutations in the 23S rRNA gene were detected in 12 clinical resistant isolates, but not in 7 clinical susceptible isolates. All sequencing reads of CLR-resistant strains had a G mutation in an identical position of the 23S rRNA gene. In addition, genetic variants of four gene clusters (hp0605-hp0607, hp0971-hp0969, hp1327-hp1329, and hp1489-hp1487) of TolC homologues, which have been implicated in multi-drug resistance, were examined. Specific SNVs were dominantly found in resistant strains.

CONCLUSIONS

Gene clusters of TolC homologues are involved in CLR susceptibility profiles in individual H. pylori strains. Whole-genome sequencing has yielded novel understanding of genotype-phenotype relationships.

Analysis of plant microbe interactions in the era of next generation sequencing technologies

Next generation sequencing (NGS) technologies have impressively accelerated research in biological science during the last years by enabling the production of large volumes of sequence data to a drastically lower price per base, compared to traditional sequencing methods. The recent and ongoing developments in the field allow addressing research questions in plant-microbe biology that were not conceivable just a few years ago. The present review provides an overview of NGS technologies and their usefulness for the analysis of microorganisms that live in association with plants. Possible limitations of the different sequencing systems, in particular sources of errors and bias, are critically discussed and methods are disclosed that help to overcome these shortcomings. A focus will be on the application of NGS methods in metagenomic studies, including the analysis of microbial communities by amplicon sequencing, which can be considered as a targeted metagenomic approach. Different applications of NGS technologies are exemplified by selected research articles that address the biology of the plant associated microbiota to demonstrate the worth of the new methods.

Complete Genome Sequences of Eight Helicobacter pylori Strains with Different Virulence Factor Genotypes and Methylation Profiles, Isolated from Patients with Diverse Gastrointestinal Diseases on Okinawa Island, Japan, Determined Using PacBio Single-Molecule Real-Time Technology

We report the complete genome sequences of eight Helicobacter pylori strains isolated from patients with gastrointestinal diseases in Okinawa, Japan. Whole-genome sequencing and DNA methylation detection were performed using the PacBio platform. De novo assembly determined a single, complete contig for each strain. Furthermore, methylation analysis identified virulence factor genotype-dependent motifs.

Live genomics for pathogen monitoring in public health

Whole genome analysis based on next generation sequencing (NGS) now represents an affordable framework in public health systems. Robust analytical pipelines of genomic data provides in short laps of time (hours) information about taxonomy, comparative genomics (pan-genome) and single polymorphisms profiles. Pathogenic organisms of interest can be tracked at the genomic level, allowing monitoring at one-time several variables including: epidemiology, pathogenicity, resistance to antibiotics, virulence, persistence factors, mobile elements and adaptation features. Such information can be obtained not only at large spectra, but also at the "local" level, such as in the event of a recurrent or emergency outbreak. This paper reviews the state of the art in infection diagnostics in the context of modern NGS methodologies. We describe how actuation protocols in a public health environment will benefit from a "streaming approach" (pipeline). Such pipeline would NGS data quality assessment, data mining for comparative analysis, searching differential genetic features, such as virulence, resistance persistence factors and mutation profiles (SNPs and InDels) and formatted "comprehensible" results. Such analytical protocols will enable a quick response to the needs of locally circumscribed outbreaks, providing information on the causes of resistance and genetic tracking elements for rapid detection, and monitoring actuations for present and future occurrences.