Benchmarking the empirical accuracy of short-read sequencing across the M. tuberculosis genome.

MOTIVATION: Short-read whole-genome sequencing (WGS) is a vital tool for clinical applications and basic research. Genetic divergence from the reference genome, repetitive sequences and sequencing bias reduces the performance of variant calling using short-read alignment, but the loss in recall and specificity has not been adequately characterized. To benchmark short-read variant calling, we used 36 diverse clinical Mycobacterium tuberculosis (Mtb) isolates dually sequenced with Illumina short-reads and PacBio long-reads. We systematically studied the short-read variant calling accuracy and the influence of sequence uniqueness, reference bias and GC content. RESULTS: Reference-based Illumina variant calling demonstrated a maximum recall of 89.0% and minimum precision of 98.5% across parameters evaluated. The approach that maximized variant recall while still maintaining high precision (<99%) was tuning the mapping quality filtering threshold, i.e. confidence of the read mapping (recall = 85.8%, precision = 99.1%, MQ ≥ 40). Additional masking of repetitive sequence content is an alternative conservative approach to variant calling that increases precision at cost to recall (recall = 70.2%, precision = 99.6%, MQ ≥ 40). Of the genomic positions typically excluded for Mtb, 68% are accurately called using Illumina WGS including 52/168 PE/PPE genes (34.5%). From these results, we present a refined list of low confidence regions across the Mtb genome, which we found to frequently overlap with regions with structural variation, low sequence uniqueness and low sequencing coverage. Our benchmarking results have broad implications for the use of WGS in the study of Mtb biology, inference of transmission in public health surveillance systems and more generally for WGS applications in other organisms. AVAILABILITY AND IMPLEMENTATION: All relevant code is available at https://github.com/farhat-lab/mtb-illumina-wgs-evaluation. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Intraspecies plasmid and genomic variation of Mycobacterium kubicae revealed by the complete genome sequences of two clinical isolates.

Mycobacterium kubicae is 1 of nearly 200 species of nontuberculous mycobacteria (NTM), environmental micro-organisms that in some situations can infect humans and cause severe lung, skin and soft tissue infections. Although numerous studies have investigated the genetic variation among prevalent clinical NTM species, including Mycobacterium abscessus and Mycobacterium avium, many of the less common but clinically relevant NTM species, including M. kubicae, still lack complete genomes to serve as a comparative reference. Well-characterized representative genomes for each NTM species are important both for investigating the pathogenic potential of NTM, as well as for use in diagnostic methods, even for species that less frequently cause human disease. Here, we report the complete genomes of two M. kubicae strains, isolated from two unrelated patients. Hybrid short-read and long-read sequencing and assembly, using sequence reads from Illumina and Oxford Nanopore Technologies platforms, were utilized to resolve the chromosome and plasmid sequences of each isolate. The genome of NJH_MKUB1 had 5135 coding sequences (CDSs), a circular chromosome of length 5.3 Mb and two plasmids. The genome of NJH_MKUB2 had 5957 CDSs, a circular chromosome of 6.0 Mb and five plasmids. We compared our completed genomic assemblies to four recently released draft genomes of M. kubicae in order to better understand intraspecies genomic conservation and variability. We also identified genes implicated in drug resistance, virulence and persistence in the M. kubicae chromosome and plasmids. Virulence factors encoded in the genome and in the plasmids of M. kubicae provide a foundation for investigating how opportunistic environmental NTM may cause disease.

In-host population dynamics of Mycobacterium tuberculosis complex during active disease.

Tuberculosis (TB) is a leading cause of death globally. Understanding the population dynamics of TB's causative agent Mycobacterium tuberculosis complex (Mtbc) in-host is vital for understanding the efficacy of antibiotic treatment. We use longitudinally collected clinical Mtbc isolates that underwent Whole-Genome Sequencing from the sputa of 200 patients to investigate Mtbc diversity during the course of active TB disease after excluding 107 cases suspected of reinfection, mixed infection or contamination. Of the 178/200 patients with persistent clonal infection >2 months, 27 developed new resistance mutations between sampling with 20/27 occurring in patients with pre-existing resistance. Low abundance resistance variants at a purity of ≥19% in the first isolate predict fixation in the subsequent sample. We identify significant in-host variation in 27 genes, including antibiotic resistance genes, metabolic genes and genes known to modulate host innate immunity and confirm several to be under positive selection by assessing phylogenetic convergence across a genetically diverse sample of 20,352 isolates.

Formulation and pharmacokinetics of self-assembled rifampicin nanoparticle systems for pulmonary delivery.

PURPOSE: To formulate rifampicin, an anti-tuberculosis antibiotic, for aerosol delivery in a dry powder 'porous nanoparticle-aggregate particle' (PNAP) form suited for shelf stability, effective dispersibility and extended release with local lung and systemic drug delivery. METHODS: Rifampicin was encapsulated in PLGA nanoparticles by a solvent evaporation process, spray dried into PNAPs containing varying amounts of nanoparticles, and characterized for physical and aerosol properties. Pharmacokinetic studies were performed with formulations delivered to guinea pigs by intratracheal insufflation and compared to oral and intravenous delivery of rifampicin. RESULTS: The PNAP formulations possessed properties suitable for efficient deposition in the lungs. In vitro release showed an initial burst of rifampicin, with the remainder available for release beyond eight hours. PNAPs delivered to guinea pigs by insufflation achieved systemic levels of rifampicin detected for six to eight hours. Moreover, rifampicin concentrations remained detectable in lung tissue and cells up to and beyond eight hours. Conversely, after pulmonary delivery of an aerosol without nanoparticles, rifampicin could not be detected in the lungs at eight hours. CONCLUSIONS: Our results indicate that rifampicin can be formulated into an aggregated nanoparticle form that, once delivered to animals, achieves systemic exposure and extends levels of drug in the lungs.

Evaluation of a Novel MALDI Biotyper Algorithm to Distinguish Mycobacterium intracellulare From Mycobacterium chimaera.

Accurate and timely mycobacterial species identification is imperative for successful diagnosis, treatment, and management of disease caused by nontuberculous mycobacteria (NTM). The current most widely utilized method for NTM species identification is Sanger sequencing of one or more genomic loci, followed by BLAST sequence analysis. MALDI-TOF MS offers a less expensive and increasingly accurate alternative to sequencing, but the commercially available assays used in clinical mycobacteriology cannot differentiate between Mycobacterium intracellulare and Mycobacterium chimaera, two closely related potentially pathogenic species of NTM that are members of the Mycobacterium avium complex (MAC). Because this differentiation of MAC species is challenging in a diagnostic setting, Bruker has developed an improved spectral interpretation algorithm to differentiate M. chimaera and M. intracellulare based on differential spectral peak signatures. Here, we utilize a set of 185 MAC isolates that have been characterized using rpoB locus sequencing followed by whole genome sequencing in some cases, to test the accuracy of the Bruker subtyper software to identify M. chimaera (n = 49) and M. intracellulare (n = 55). 100% of the M. intracellulare and 82% of the M. chimaera isolates were accurately identified using the MALDI Biotyper algorithm. This subtyper module is available with the MALDI Biotyper Compass software and offers a promising mechanism for rapid and inexpensive species determination for M. chimaera and M. intracellulare.