Rapid and accurate SNP genotyping of clonal bacterial pathogens with BioHansel

Seamless, rapid, and accurate analyses of outbreak genomic data using split k-mer analysis

Sequence variation observed in populations of pathogens can be used for important public health and evolutionary genomic analyses, especially outbreak analysis and transmission reconstruction. Identifying this variation is typically achieved by aligning sequence reads to a reference genome, but this approach is susceptible to reference biases and requires careful filtering of called genotypes. There is a need for tools that can process this growing volume of bacterial genome data, providing rapid results, but that remain simple so they can be used without highly trained bioinformaticians, expensive data analysis, and long-term storage and processing of large files. Here we describe split k-mer analysis (SKA2), a method that supports both reference-free and reference-based mapping to quickly and accurately genotype populations of bacteria using sequencing reads or genome assemblies. SKA2 is highly accurate for closely related samples, and in outbreak simulations, we show superior variant recall compared with reference-based methods, with no false positives. SKA2 can also accurately map variants to a reference and be used with recombination detection methods to rapidly reconstruct vertical evolutionary history. SKA2 is many times faster than comparable methods and can be used to add new genomes to an existing call set, allowing sequential use without the need to reanalyze entire collections. With an inherent absence of reference bias, high accuracy, and a robust implementation, SKA2 has the potential to become the tool of choice for genotyping bacteria. SKA2 is implemented in Rust and is freely available as open-source software.

Average Nucleotide Identity and Digital DNA-DNA Hybridization Analysis Following PromethION Nanopore-Based Whole Genome Sequencing Allows for Accurate Prokaryotic Typing

Whole-genome sequencing (WGS) is revolutionizing clinical bacteriology. However, bacterial typing remains investigated by reference techniques with inherent limitations. This stresses the need for alternative methods providing robust and accurate sequence type (ST) classification. This study optimized and evaluated a GridION nanopore sequencing protocol, adapted for the PromethION platform. Forty-eight Escherichia coli clinical isolates with diverse STs were sequenced to assess two alternative typing methods and resistance profiling applications. Multi-locus sequence typing (MLST) was used as the reference typing method. Genomic relatedness was assessed using Average Nucleotide Identity (ANI) and digital DNA-DNA Hybridization (DDH), and cut-offs for discriminative strain resolution were evaluated. WGS-based antibiotic resistance prediction was compared to reference Minimum Inhibitory Concentration (MIC) assays. We found ANI and DDH cut-offs of 99.3% and 94.1%, respectively, which correlated well with MLST classifications and demonstrated potentially higher discriminative resolution than MLST. WGS-based antibiotic resistance prediction showed categorical agreements of ≥ 93% with MIC assays for amoxicillin, ceftazidime, amikacin, tobramycin, and trimethoprim-sulfamethoxazole. Performance was suboptimal (68.8-81.3%) for amoxicillin-clavulanic acid, cefepime, aztreonam, and ciprofloxacin. A minimal sequencing coverage of 12× was required to maintain essential genomic features and typing accuracy. Our protocol allows the integration of PromethION technology in clinical laboratories, with ANI and DDH proving to be accurate and robust alternative typing methods, potentially offering superior resolution. WGS-based antibiotic resistance prediction holds promise for specific antibiotic classes.

Whole genome sequencing-based identification of human tuberculosis caused by animal-lineage Mycobacterium orygis

A recently described member of the Mycobacterium tuberculosis complex (MTBC) is Mycobacterium orygis, which can cause disease primarily in animals but also in humans. Although M. orygis has been reported from different geographic regions around the world, due to a lack of proper identification techniques, the contribution of this emerging pathogen to the global burden of zoonotic tuberculosis is not fully understood. In the present work, we report single nucleotide polymorphism (SNP) analysis using whole genome sequencing (WGS) that can accurately identify M. orygis and differentiate it from other members of the MTBC species. WGS-based SNP analysis was performed for 61 isolates from different provinces in Canada that were identified as M. orygis. A total of 56 M. orygis sequences from the public databases were also included in the analysis. Several unique SNPs in the gyrB, PPE55, Rv2042c, leuS, mmpL6, and mmpS6 genes were used to determine their effectiveness as genetic markers for the identification of M. orygis. To the best of our knowledge, five of these SNPs, viz., gyrB 277 (A→G), gyrB 1478 (T→C), leuS 1064 (A→T), mmpL6 486 (T→C), and mmpS6 334 (C→G), are reported for the first time in this study. Our results also revealed several SNPs specific to other species within MTBC. The phylogenetic analysis shows that the studied genomes were genetically diverse and clustered with M. orygis sequences of human and animal origin reported from different geographic locations. Therefore, the present study provides a new insight into the high-confidence identification of M. orygis from MTBC species based on WGS data, which can be useful for reference and diagnostic laboratories.

Whole-genome analyses reveal a novel prophage and cgSNPs-derived sublineages of Brachyspira hyodysenteriae ST196

Background

Brachyspira (B.) hyodysenteriae is a fastidious anaerobe spirochete that can cause swine dysentery, a severe mucohaemorragic colitis that affects pig production and animal welfare worldwide. In Switzerland, the population of B. hyodysenteriae is characterized by the predominance of macrolide-lincosamide-resistant B. hyodysenteriae isolates of sequence type (ST) ST196, prompting us to obtain deeper insights into the genomic structure and variability of ST196 using pangenome and whole genome variant analyses.

Results

The draft genome of 14 B. hyodysenteriae isolates of ST196, sampled during a 7-year period from geographically distant pig herds, was obtained by whole-genome sequencing (WGS) and compared to the complete genome of the B. hyodysenteriae isolate Bh743-7 of ST196 used as reference. Variability results revealed the existence of 30 to 52 single nucleotide polymorphisms (SNPs), resulting in eight sublineages of ST196. The pangenome analysis led to the identification of a novel prophage, pphBhCH20, of the Siphoviridae family in a single isolate of ST196, which suggests that horizontal gene transfer events may drive changes in genomic structure.

Conclusions

This study contributes to the catalogue of publicly available genomes and provides relevant bioinformatic tools and information for further comparative genomic analyses for B. hyodysenteriae. It reveals that Swiss B. hyodysenteriae isolates of the same ST may have evolved independently over time by point mutations and acquisition of larger genetic elements. In line with this, the third type of mobile genetic element described so far in B. hyodysenteriae, the novel prophage pphBhCH20, has been identified in a single isolate of B. hyodysenteriae of ST196.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08347-5.

How much do smear-negative patients really contribute to tuberculosis transmissions? Re-examining an old question with new tools

BACKGROUND

Sputum smear microscopy is a common surrogate for tuberculosis infectiousness. Previous estimates that smear-negative patients contribute 13-20% of transmissions and are, on average, 20 to 25% as infectious as smear-positive cases are understood to be high. Herein, we use an ideal real-world setting, a comprehensive dataset, and new high-resolution techniques to more accurately estimate the true transmission risk of smear-negative cases.

METHODS

We treated all adult culture-positive pulmonary TB patients diagnosed in the province of Alberta, Canada from 2003 to 2016 as potential transmitters. The primary data sources were the Alberta TB Registry and the Provincial Laboratory for Public Health. We measured, as primary outcomes, the proportion of transmissions attributable to smear-negative sources and the relative transmission rate. First, we replicated previous studies by using molecular (DNA) fingerprint clustering. Then, using a prospectively collected registry of TB contacts, we defined transmission events as active TB amongst identified contacts who either had a 100% DNA fingerprint match to the source case or a clinical diagnosis. We supplemented our analysis with genome sequencing on temporally and geographically linked DNA fingerprint clusters of cases not identified as contacts.

FINDINGS

There were 1176 cases, 563 smear-negative and 613 smear-positive, and 23,131 contacts. Replicating previous studies, the proportion of transmissions attributable to smear-negative source cases was 16% (95% CI, 12-19%) and the relative transmission rate was 0.19 (95% CI, 0.14-0.26). With our combined approach, the proportion of transmission was 8% (95% CI, 3-14%) and the relative transmission rate became 0.10 (95% CI, 0.05-0.19).

INTERPRETATION

When we examined the same outcomes as in previous studies but refined transmission ascertainment with the addition of conventional epidemiology and genomics, we found that smear-negative cases were ∼50% less infectious than previously thought.

FUNDING

Alberta Innovates Health Solutions.