Next Generation Sequencing in Clinical and Public Health Microbiology

Helicobacter pylori Infection, Its Laboratory Diagnosis, and Antimicrobial Resistance: a Perspective of Clinical Relevance

Despite the recent decrease in overall prevalence of Helicobacter pylori infection, morbidity and mortality rates associated with gastric cancer remain high. The antimicrobial resistance developments and treatment failure are fueling the global burden of H. pylori-associated gastric complications. Accurate diagnosis remains the opening move for treatment and eradication of infections caused by microorganisms. Although several reports have been published on diagnostic approaches for H. pylori infection, most lack the data regarding diagnosis from a clinical perspective. Therefore, we provide an intensive, comprehensive, and updated description of the currently available diagnostic methods that can help clinicians, infection diagnosis professionals, and H. pylori researchers working on infection epidemiology to broaden their understanding and to select appropriate diagnostic methods. We also emphasize appropriate diagnostic approaches based on clinical settings (either clinical diagnosis or mass screening), patient factors (either age or other predisposing factors), and clinical factors (either upper gastrointestinal bleeding or partial gastrectomy) and appropriate methods to be considered for evaluating eradication efficacy. Furthermore, to cope with the increasing trend of antimicrobial resistance, a better understanding of its emergence and current diagnostic approaches for resistance detection remain inevitable.

Helicobacter pylori infection and antibiotic resistance - from biology to clinical implications

Helicobacter pylori is a major human pathogen for which increasing antibiotic resistance constitutes a serious threat to human health. Molecular mechanisms underlying this resistance have been intensively studied and are discussed in this Review. Three profiles of resistance - single drug resistance, multidrug resistance and heteroresistance - seem to occur, probably with overlapping fundamental mechanisms and clinical implications. The mechanisms that have been most studied are related to mutational changes encoded chromosomally and disrupt the cellular activity of antibiotics through target-mediated mechanisms. Other biological attributes driving drug resistance in H. pylori have been less explored and this could imply more complex physiological changes (such as impaired regulation of drug uptake and/or efflux, or biofilm and coccoid formation) that remain largely elusive. Resistance-related attributes deployed by the pathogen cause treatment failures, diagnostic difficulties and ambiguity in clinical interpretation of therapeutic outcomes. Subsequent to the increasing antibiotic resistance, a substantial drop in H. pylori treatment efficacy has been noted globally. In the absence of an efficient vaccine, enhanced efforts are needed for setting new treatment strategies and for a better understanding of the emergence and spread of drug-resistant bacteria, as well as for improving diagnostic tools that can help optimize current antimicrobial regimens.

Feasibility and potential significance of rapid in vitro qualitative phenotypic antimicrobial susceptibility testing of gram-negative bacilli with the ProMax system

The emergence and evolution of antibiotic resistance has been accelerated due to the widespread use of antibiotics and a lack of timely diagnostic tests that guide therapeutic treatment with adequate sensitivity, specificity, and antimicrobial susceptibility testing (AST) accuracy. Automated AST instruments are extensively used in clinical microbiology labs and provide a streamlined workflow, simplifying susceptibility testing for pathogenic bacteria isolated from clinical samples. Although currently used commercial systems such as the Vitek2 and BD Phoenix can deliver results in substantially less time than conventional methods, their dependence on traditional AST inoculum concentrations and optical detection limit their speed somewhat. Herein, we describe the GeneFluidics ProMax lab automation system intended for a rapid 3.5-hour molecular AST from clinical isolates. The detection method described utilizes a higher starting inoculum concentration and automated molecular quantification of species-specific 16S rRNA through the use of an electrochemical sensor to assess microbiological responses to antibiotic exposure. A panel of clinical isolates consisting of species of gram-negative rods from the CDC AR bank and two hospitals, New York-Presbyterian Queens and Medical College of Wisconsin, were evaluated against ciprofloxacin, gentamicin, and meropenem in a series of reproducibility and clinical studies. The categorical agreement and reproducibility for Citrobacter freundii, Enterobacter cloacae, Escherichia coli, Klebsiella aerogenes, Klebsiella oxytoca, Klebsiella pneumoniae, and Pseudomonas aeruginosa were 100% and 100% for ciprofloxacin, 98.7% and 100% for gentamicin and 98.5% and 98.5% for meropenem, respectively.

A novel DNA chromatography method to discriminate Mycobacterium abscessus subspecies and macrolide susceptibility

Background

The clinical impact of infection with Mycobacterium (M.) abscessus complex (MABC), a group of emerging non-tuberculosis mycobacteria (NTM), is increasing. M. abscessus subsp. abscessus/bolletii frequently shows natural resistance to macrolide antibiotics, whereas M. abscessus subsp. massiliense is generally susceptible. Therefore, rapid and accurate discrimination of macrolide-susceptible MABC subgroups is required for effective clinical decisions about macrolide treatments for MABC infection. We aimed to develop a simple and rapid diagnostic that can identify MABC isolates showing macrolide susceptibility.

Methods

Whole genome sequencing (WGS) was performed for 148 clinical or environmental MABC isolates from Japan to identify genetic markers that can discriminate three MABC subspecies and the macrolide-susceptible erm(41) T28C sequevar. Using the identified genetic markers, we established PCR based- or DNA chromatography-based assays. Validation testing was performed using MABC isolates from Taiwan.

Finding

We identified unique sequence regions that could be used to differentiate the three subspecies. Our WGS-based phylogenetic analysis indicated that M. abscessus carrying the macrolide-susceptible erm(41) T28C sequevar were tightly clustered, and identified 11 genes that were significantly associated with the lineage for use as genetic markers. To detect these genetic markers and the erm(41) locus, we developed a DNA chromatography method that identified three subspecies, the erm(41) T28C sequevar and intact erm(41) for MABC in a single assay within one hour. The agreement rate between the DNA chromatography-based and WGS-based identification was 99·7%.

Interpretation

We developed a novel, rapid and simple DNA chromatography method for identification of MABC macrolide susceptibility with high accuracy.

Funding

AMED, JSPS KAKENHI

Next-Generation Sequencing of the Whole Bacterial Genome for Tracking Molecular Insight into the Broad-Spectrum Antimicrobial Resistance of Helicobacter pylori Clinical Isolates from the Democratic Republic of Congo

Antimicrobial susceptibility testing (AST) is increasingly needed to guide the Helicobacter pylori (H. pylori) treatment but remains laborious and unavailable in most African countries. To assess the clinical relevance of bacterial whole genome sequencing (WGS)-based methods for predicting drug susceptibility in African H. pylori, 102 strains isolated from the Democratic Republic of Congo were subjected to the phenotypic AST and next-generation sequencing (NGS). WGS was used to screen for the occurrence of genotypes encoding antimicrobial resistance (AMR). We noted the broad-spectrum AMR of H. pylori (rates from 23.5 to 90.0%). A WGS-based method validated for variant discovery in AMR-related genes (discovery rates of 100%) helped in identifying mutations of key genes statistically related to the phenotypic AMR. These included mutations often reported in Western and Asian populations and, interestingly, several putative AMR-related new genotypes in the pbp1A (e.g., T558S, F366L), gyrA (e.g., A92T, A129T), gyrB (e.g., R579C), and rdxA (e.g., R131_K166del) genes. WGS showed high performance for predicting AST phenotypes, especially for amoxicillin, clarithromycin, and levofloxacin (Youden's index and Cohen's Kappa > 0.80). Therefore, WGS is an accurate alternative to the phenotypic AST that provides substantial decision-making information for public health policy makers and clinicians in Africa, while providing insight into AMR mechanisms for researchers.

Genomics for Molecular Epidemiology and Detecting Transmission of Carbapenemase-Producing Enterobacterales in Victoria, Australia, 2012 to 2016

Carbapenemase-producing Enterobacterales (CPE) are being increasingly reported in Australia, and integrated clinical and genomic surveillance is critical to effectively manage this threat. We sought to systematically characterize CPE in Victoria, Australia, from 2012 to 2016.

Carbapenemase-producing Enterobacterales (CPE) are being increasingly reported in Australia, and integrated clinical and genomic surveillance is critical to effectively manage this threat. We sought to systematically characterize CPE in Victoria, Australia, from 2012 to 2016. Suspected CPE were referred to the state public health laboratory in Victoria, Australia, from 2012 to 2016 and examined using phenotypic, multiplex PCR and whole-genome sequencing (WGS) methods and compared with epidemiological metadata. Carbapenemase genes were detected in 361 isolates from 291 patients (30.8% of suspected CPE isolates), mostly from urine (42.1%) or screening samples (34.8%). IMP-4 (28.0% of patients), KPC-2 (25.3%), NDM (24.1%), and OXA carbapenemases (22.0%) were most common. Klebsiella pneumoniae (48.8% of patients) and Escherichia coli (26.1%) were the dominant species. Carbapenemase-inactivation method (CIM) testing reliably detected carbapenemase-positive isolates (100% sensitivity, 96.9% specificity), identifying an additional five CPE among 159 PCR-negative isolates (IMI and SME carbapenemases). When epidemiologic investigations were performed, all pairs of patients designated “highly likely” or “possible” local transmission had ≤23 pairwise single-nucleotide polymorphisms (SNPs) by genomic transmission analysis; conversely, all patient pairs designated “highly unlikely” local transmission had ≥26 pairwise SNPs. Using this proposed threshold, possible local transmission was identified involving a further 16 patients for whom epidemiologic data were unavailable. Systematic application of genomics has uncovered the emergence of polyclonal CPE as a significant threat in Australia, providing important insights to inform local public health guidelines and interventions. Using our workflow, pairwise SNP distances between CPE isolates of ≤23 SNPs suggest local transmission.

Status and potential of bacterial genomics for public health practice: a scoping review

BACKGROUND

Next-generation sequencing (NGS) is increasingly being translated into routine public health practice, affecting the surveillance and control of many pathogens. The purpose of this scoping review is to identify and characterize the recent literature concerning the application of bacterial pathogen genomics for public health practice and to assess the added value, challenges, and needs related to its implementation from an epidemiologist's perspective.

METHODS

In this scoping review, a systematic PubMed search with forward and backward snowballing was performed to identify manuscripts in English published between January 2015 and September 2018. Included studies had to describe the application of NGS on bacterial isolates within a public health setting. The studied pathogen, year of publication, country, number of isolates, sampling fraction, setting, public health application, study aim, level of implementation, time orientation of the NGS analyses, and key findings were extracted from each study. Due to a large heterogeneity of settings, applications, pathogens, and study measurements, a descriptive narrative synthesis of the eligible studies was performed.

RESULTS

Out of the 275 included articles, 164 were outbreak investigations, 70 focused on strategy-oriented surveillance, and 41 on control-oriented surveillance. Main applications included the use of whole-genome sequencing (WGS) data for (1) source tracing, (2) early outbreak detection, (3) unraveling transmission dynamics, (4) monitoring drug resistance, (5) detecting cross-border transmission events, (6) identifying the emergence of strains with enhanced virulence or zoonotic potential, and (7) assessing the impact of prevention and control programs. The superior resolution over conventional typing methods to infer transmission routes was reported as an added value, as well as the ability to simultaneously characterize the resistome and virulome of the studied pathogen. However, the full potential of pathogen genomics can only be reached through its integration with high-quality contextual data.

CONCLUSIONS

For several pathogens, it is time for a shift from proof-of-concept studies to routine use of WGS during outbreak investigations and surveillance activities. However, some implementation challenges from the epidemiologist's perspective remain, such as data integration, quality of contextual data, sampling strategies, and meaningful interpretations. Interdisciplinary, inter-sectoral, and international collaborations are key for an appropriate genomics-informed surveillance.

Inguinal microbiome in patients undergoing an endovascular aneurysm repair: Application of next-generation sequencing of the 16S-23S rRNA regions

BACKGROUND

Surgical site infection (SSI) remains a hazardous complication after vascular surgery. In this pilot study we investigated the inguinal microbiome in skin biopsies using histology and 16S-23S rDNA Next Generation Sequencing (NGS). Our hypothesis was that causative microorganisms of SSI are present in the inguinal microbiome.

METHODS

Data on surgical site infections and skin samples from the Percutaneous in Endovascular Repair versus Open (PiERO) trail were evaluated. Two patients with SSI were matched for age and comorbidity to eight matching patients of the PiERO trial. All patients were treated for an abdominal aortic aneurysm with endovascular repair. Nasal and perineal cultures were taken preoperatively to detect Staphylococcus aureus carriage. After disinfection with chlorhexidine, groin biopsies were taken to identify bacteria in deeper skin layers. All samples were subjected to histological analysis and culture-free 16S-23S rDNA NGS.

RESULTS

Staphylococcus aureus species were cultured in 5 out of 20 preoperative nasal and perineal swaps. Histology detected only a few bacteria. NGS of the 16S-23S rRNA regions identified DNA of bacterial species in all biopsies (20/20). Most identified genera and species proved to be known skin flora bacteria. No relation was found between SSIs and the preoperative microbiome.

CONCLUSION

In this pilot study, an innovative analysis of the preoperative microbiome using 16S-23S rDNA NGS did not show a relation with the occurrence of a surgical site infection. No pathogenic bacterial species were present in the inguinal skin after disinfection with chlorhexidine.

Identification and characterization of putative Aeromonas spp. T3SS effectors

The genetic determinants of bacterial pathogenicity are highly variable between species and strains. However, a factor that is commonly associated with virulent Gram-negative bacteria, including many Aeromonas spp., is the type 3 secretion system (T3SS), which is used to inject effector proteins into target eukaryotic cells. In this study, we developed a bioinformatics pipeline to identify T3SS effector proteins, applied this approach to the genomes of 105 Aeromonas strains isolated from environmental, mutualistic, or pathogenic contexts and evaluated the cytotoxicity of the identified effectors through their heterologous expression in yeast. The developed pipeline uses a two-step approach, where candidate Aeromonas gene families are initially selected using Hidden Markov Model (HMM) profile searches against the Virulence Factors DataBase (VFDB), followed by strict comparisons against positive and negative control datasets, greatly reducing the number of false positives. This approach identified 21 Aeromonas T3SS likely effector families, of which 8 represent known or characterized effectors, while the remaining 13 have not previously been described in Aeromonas. We experimentally validated our in silico findings by assessing the cytotoxicity of representative effectors in Saccharomyces cerevisiae BY4741, with 15 out of 21 assayed proteins eliciting a cytotoxic effect in yeast. The results of this study demonstrate the utility of our approach, combining a novel in silico search method with in vivo experimental validation, and will be useful in future research aimed at identifying and authenticating bacterial effector proteins from other genera.

A Comparison of Three Different Bioinformatics Analyses of the 16S-23S rRNA Encoding Region for Bacterial Identification

Rapid and reliable identification of bacterial pathogens directly from patient samples is required for optimizing antimicrobial therapy. Although Sanger sequencing of the 16S ribosomal RNA (rRNA) gene is used as a molecular method, species identification and discrimination is not always achievable for bacteria as their 16S rRNA genes have sometimes high sequence homology. Recently, next generation sequencing (NGS) of the 16S–23S rRNA encoding region has been proposed for reliable identification of pathogens directly from patient samples. However, data analysis is laborious and time-consuming and a database for the complete 16S–23S rRNA encoding region is not available. Therefore, a better, faster, and stronger approach is needed for NGS data analysis of the 16S–23S rRNA encoding region. We compared speed and diagnostic accuracy of different data analysis approaches: de novo assembly followed by Basic Local Alignment Search Tool (BLAST), operational taxonomic unit (OTU) clustering, or mapping using an in-house developed 16S–23S rRNA encoding region database for the identification of bacterial species. De novo assembly followed by BLAST using the in-house database was superior to the other methods, resulting in the shortest turnaround time (2 h and 5 min), approximately 2 h less than OTU clustering and 4.5 h less than mapping, and a sensitivity of 80%. Mapping was the slowest and most laborious data analysis approach with a sensitivity of 60%, whereas OTU clustering was the least laborious approach with 70% sensitivity. Although the in-house database requires more sequence entries to improve the sensitivity, the combination of de novo assembly and BLAST currently appears to be the optimal approach for data analysis.

Whole-genome and targeted sequencing of drug-resistant Mycobacterium tuberculosis on the iSeq100 and MiSeq: A performance, ease-of-use, and cost evaluation

BACKGROUND

Accurate, comprehensive, and timely detection of drug-resistant tuberculosis (TB) is essential to inform patient treatment and enable public health surveillance. This is crucial for effective control of TB globally. Whole-genome sequencing (WGS) and targeted next-generation sequencing (NGS) approaches have potential as rapid in vitro diagnostics (IVDs), but the complexity of workflows, interpretation of results, high costs, and vulnerability of instrumentation have been barriers to broad uptake outside of reference laboratories, especially in low- and middle-income countries. A new, solid-state, tabletop sequencing instrument, Illumina iSeq100, has the potential to decentralize NGS for individual patient care.

METHODS AND FINDINGS

In this study, we evaluated WGS and targeted NGS for TB on both the new iSeq100 and the widely used MiSeq (both manufactured by Illumina) and compared sequencing performance, costs, and usability. We utilized DNA libraries produced from Mycobacterium tuberculosis clinical isolates for the evaluation. We conducted WGS on three strains and observed equivalent uniform genome coverage with both platforms and found the depth of coverage obtained was consistent with the expected data output. Utilizing the standardized, cloud-based ReSeqTB bioinformatics pipeline for variant analysis, we found the two platforms to have 94.0% (CI 93.1%-94.8%) agreement, in comparison to 97.6% (CI 97%-98.1%) agreement for the same libraries on two MiSeq instruments. For the targeted NGS approach, 46 M. tuberculosis-specific amplicon libraries had 99.6% (CI 98.0%-99.9%) agreement between the iSeq100 and MiSeq data sets in drug resistance-associated SNPs. The upfront capital costs are almost 5-fold lower for the iSeq100 ($19,900 USD) platform in comparison to the MiSeq ($99,000 USD); however, because of difference in the batching capabilities, the price per sample for WGS was higher on the iSeq100. For WGS of M. tuberculosis at the minimum depth of coverage of 30x, the cost per sample on the iSeq100 was $69.44 USD versus $28.21 USD on the MiSeq, assuming a 2 × 150 bp run on a v3 kit. In terms of ease of use, the sequencing workflow of iSeq100 has been optimized to only require 27 minutes total of hands-on time pre- and post-run, and the maintenance is simplified by a single-use cartridge-based fluidic system. As these are the first sequencing attempts on the iSeq100 for M. tuberculosis, the sequencing pool loading concentration still needs optimization, which will affect sequencing error and depth of coverage. Additionally, the costs are based on current equipment and reagent costs, which are subject to change.

CONCLUSIONS

The iSeq100 instrument is capable of running existing TB WGS and targeted NGS library preparations with comparable accuracy to the MiSeq. The iSeq100 has reduced sequencing workflow hands-on time and is able to deliver sequencing results in <24 hours. Reduced capital and maintenance costs and lower-throughput capabilities also give the iSeq100 an advantage over MiSeq in settings of individualized care but not in high-throughput settings such as reference laboratories, where sample batching can be optimized to minimize cost at the expense of workflow complexity and time.

Platforms and Analytical Tools Used in Nucleic Acid Sequence-Based Microbial Genotyping Procedures

In the decade and a half since the introduction of next-generation sequencing (NGS), the technical feasibility, cost, and overall utility of sequencing have changed dramatically, including applications for infectious disease epidemiology. Massively parallel sequencing technologies have decreased the cost of sequencing by more than 6 orders or magnitude over this time, with a corresponding increase in data generation and complexity. This review provides an overview of the basic principles, chemistry, and operational mechanics of current sequencing technologies, including both conventional Sanger and NGS approaches. As the generation of large amounts of sequence data becomes increasingly routine, the role of bioinformatics in data analysis and reporting becomes all the more critical, and the successful deployment of NGS in public health settings requires careful consideration of changing information technology, bioinformatics, workforce, and regulatory requirements. While there remain important challenges to the sustainable implementation of NGS in public health, in terms of both laboratory and bioinformatics capacity, the impact of these technologies on infectious disease surveillance and outbreak investigations has been nothing short of revolutionary. Understanding the important role that NGS plays in modern public health laboratory practice is critical, as is the need to ensure appropriate workforce, infrastructure, facilities, and funding consideration for routine NGS applications, future innovation, and rapidly scaling NGS-based infectious disease surveillance and outbreak response activities. *This article is part of a curated collection.

HCV adaptation to HIV coinfection

Human immunodeficiency virus (HIV) infection is rising as a leading cause of morbidity and mortality among hepatitis C virus (HCV)-infected patients. Both viruses interact in co-infected hosts, which may affect their intra-host evolution, potentially leading to differing genetic composition of viral populations in co-infected (CIP) and mono-infected (MIP) patients. Here, we investigate genetic differences between intra-host variants of the HCV hypervariable region 1 (HVR1) sampled from CIP and MIP. Nucleotide (nt) sequences of intra-host HCV HVR1 variants (N = 28,622) obtained from CIP (N = 112) and MIP (n = 176) were represented using 148 physical-chemical (PhyChem) indexes of DNA nt dimers. Significant (p < .0001) differences in the means and frequency distributions of 7 PhyChem properties were found between HVR1 variants from both groups. Linear projection analysis of 29 PhyChem features extracted from such PhyChem properties showed that the CIP and MIP HVR1 variants have a distinct distribution in the modeled 2D-space, with only ~1.3% of PhyChem profiles (N = 6782), shared by all HVR1 variants, being found in both groups. Probabilistic neural network (PNN) and naïve Bayesian (NB) classifiers trained on the PhyChem features accurately classified HVR1 variants by the group in cross-validation experiments (AUROC ≥ 0.96). Similarly, both models showed a high accuracy (AUROC ≥ 0.95) when evaluated on a test dataset of HVR1 sequences obtained from 10 patients, data from whom were not used for model building. Both models performed at the expected lower accuracy on randomly labeled datasets in cross-validation experiments (AUROC = 0.50). The random-label trained PNN showed a similar drop in accuracy on the test dataset (AUROC = 0.48), indicating that the detected associations were unlikely due to random correlations. Marked differences in genetic composition of HCV HVR1 variants sampled from CIP and MIP suggest differing intra-host HCV evolution in the presence of HIV infection. PhyChem features identified here may be used for detection of HIV infection from intra-host HCV variants alone in co-infected patients, thus facilitating monitoring for HIV introduction to high-risk populations with high HCV prevalence.

Next-generation sequencing applications in clinical bacteriology

With the rapid advances in next generation sequencing (NGS) technologies, clinical and public health microbiology laboratories are increasingly adopting NGS technology in their workflows into their existing diagnostic cycles. In this bacteriology focused review, we review aspects and considerations for applying NGS in the clinical microbiology settings, and highlight the impact of such implementation on the analytical and post-analytical stages of diagnosis

Right-Sizing Technology in the Era of Consumer-Driven Health Care

Technology for modern clinical and public health microbiology laboratories has evolved at an impressive rate over the last two decades. Contemporary diagnostics can rapidly provide powerful data that can impact patient lives and support infectious disease outbreak investigations. At the same time, dramatic changes to health care delivery are putting new pressures on a system that is now focusing on patient-centric, value-driven, convenient care. For laboratories, balancing all these demands in a cost-contained environment remains a challenge. This article explores the current and future directions of diagnostics in our dynamic health care environment.