Minos: variant adjudication and joint genotyping of cohorts of bacterial genomes

Short-read variant calling for bacterial genomics is a mature field, and there are many widely-used software tools. Different underlying approaches (eg pileup, local or global assembly, paired-read use, haplotype use) lend each tool different strengths, especially when considering non-SNP (single nucleotide polymorphism) variation or potentially distant reference genomes. It would therefore be valuable to be able to integrate the results from multiple variant callers, using a robust statistical approach to "adjudicate" at loci where there is disagreement between callers. To this end, we present a tool, Minos, for variant adjudication by mapping reads to a genome graph of variant calls. Minos allows users to combine output from multiple variant callers without loss of precision. Minos also addresses a second problem of joint genotyping SNPs and indels in bacterial cohorts, which can also be framed as an adjudication problem. We benchmark on 62 samples from 3 species (Mycobacterium tuberculosis, Staphylococcus aureus, Klebsiella pneumoniae) and an outbreak of 385 M. tuberculosis samples. Finally, we joint genotype a large M. tuberculosis cohort (N≈15k) for which the rifampicin phenotype is known. We build a map of non-synonymous variants in the RRDR (rifampicin resistance determining region) of the rpoB gene and extend current knowledge relating RRDR SNPs to heterogeneity in rifampicin resistance levels. We replicate this finding in a second M. tuberculosis cohort (N≈13k). Minos is released under the MIT license, available at https://github.com/iqbal-lab-org/minos.

Brucella Genomics: Macro and Micro Evolution

Brucella organisms are responsible for one of the most widespread bacterial zoonoses, named brucellosis. The disease affects several species of animals, including humans. One of the most intriguing aspects of the brucellae is that the various species show a ~97% similarity at the genome level. Still, the distinct Brucella species display different host preferences, zoonotic risk, and virulence. After 133 years of research, there are many aspects of the Brucella biology that remain poorly understood, such as host adaptation and virulence mechanisms. A strategy to understand these characteristics focuses on the relationship between the genomic diversity and host preference of the various Brucella species. Pseudogenization, genome reduction, single nucleotide polymorphism variation, number of tandem repeats, and mobile genetic elements are unveiled markers for host adaptation and virulence. Understanding the mechanisms of genome variability in the Brucella genus is relevant to comprehend the emergence of pathogens.

Use of whole genome sequencing to complement characterisation of a typhoid fever outbreak among a Marshallese community: Oklahoma, 2015

Typhoid fever is an illness caused bySalmonella entericaserotype Typhi. In developing regions, it affects an estimated 20 million people annually, causing 200 000 deaths. Although uncommon, cases occur in the USA each year, predominantly due to international travel. During February 2015, the Oklahoma State Department of Health (OSDH) detected an outbreak of typhoid fever among residents of northwestern Oklahoma. OSDH conducted case-patient interviews to identify the source and symptomatic contacts. Whole genome sequencing (WGS) was performed to characterise the genetic relatedness of isolates among the four outbreak-associated pulsed-field gel electrophoresis (PFGE) patterns. We identified 38 cases, 25 confirmed and 13 probable, in two states. WGS revealed a 0–10 single-nucleotide polymorphism variation between isolates. Although we were unable to determine the source, almost all case-patients were members of the Marshallese community that attended a common event in Oklahoma, or were contacts to a confirmed case. This is the largest outbreak of typhoid fever in the USA since 1989, and first to apply WGS to complement interpretation of PFGE results during a typhoid fever outbreak investigation. This investigation illustrates the potential risk of outbreaks among communities comprised of international populations from regions where typhoid fever remains endemic.