Whole-genome sequencing is increasingly used in clinical diagnosis of tuberculosis and study of Mycobacterium tuberculosis complex (MTC). MTC consists of several genetically homogenous mycobacteria species which can cause tuberculosis in humans and animals. Regions of difference (RDs) are commonly regarded as gold standard genetic markers for MTC classification. We develop RD-Analyzer, a tool that can accurately infer the species and lineage of MTC isolates from sequence reads based on the presence and absence of a set of 31 RDs. Applied on a publicly available diverse set of 377 sequenced MTC isolates from known major species and lineages, RD-Analyzer achieved an accuracy of 98.14 % (370/377) in species prediction and a concordance of 98.47 % (257/261) in Mycobacterium tuberculosis lineage prediction compared to predictions based on single nucleotide polymorphism markers. By comparing respective sequencing read depths on each genomic position between isolates of different sublineages, we were able to identify the known RD markers in different sublineages of Lineage 4 and provide support for six potential delineating markers having high sensitivities and specificities for sublineage prediction. An extended version of RD-Analyzer was thus developed to allow user-defined RDs for lineage prediction. RD-Analyzer is a useful and accurate tool for species, lineage and sublineage prediction using known RDs of MTC from sequence reads and is extendable to accepting user-defined RDs for analysis. RD-Analyzer is written in Python and is freely available at https://github.com/xiaeryu/RD-Analyzer .

Mycobacterium abscessus is an important infectious agent highly associated with drug resistance and treatment failure. We investigated the drug resistance situation of M. abscessus in Northeast Thailand and the possible genetic basis for this. Sixty-eight M. abscessus clinical isolates were obtained from 26 patients at Srinagarind Hospital during 2012–2016. Drug susceptibility tests and sequencing of erm(41), rrl and rrs genes were performed. Mycobacterium abscessus was resistant to 11/15 antibiotics (nearly 100% resistance in each case). Partial susceptibility to four antibiotics was found (amikacin, tigecycline, clarithromycin and linezolid). Non-massiliense subspecies were significantly associated with clarithromycin resistance (p<0.0001) whereas massiliense subspecies were associated with tigecycline resistance (p = 0.028). Inducible clarithromycin resistance was seen in 22/68 (32.35%) isolates: 21 of these isolates (95.45%) belonged to non-massiliense subspecies and resistance was explicable by the T28C mutation in erm(41). Inducible clarithromycin resistance was found in one isolate of the massiliense subspecies. Acquired clarithromycin resistance explicable by the A2271G/C mutation of rrl was seen in only 7/16 (43.75%) of strains. Inducible and acquired resistance mechanisms can be interchangeable during the course of infection. Rrs mutations were not associated with amikacin resistance in our study. Antibiotic resistance in subspecies of M. abscessus was reported from Northeast Thailand. Known resistance-associated mutations cannot explain all of the resistance patterns observed.

Two-thirds of the world’s population is thought to be infected by Helicobacter pylori. Although most people infected with H. pylori are asymptomatic, this pathogen is associated with several gastric pathologies including cancer. The risk factors for colonization are still unclear and the genetic diversity within individual hosts has never been clearly investigated. This study determined the prevalence of, and explored risk factors for, H. pylori infection directly from paired saliva (n = 110) and stool (n = 110) samples from asymptomatic persons in Northeast Thailand. Samples were subjected to indirect immunofluorescence assay (IFA), 16S rRNA-based real-time PCR and vacA-based semi-nested PCR. Partial vacA gene sequences of H. pylori were compared between saliva and stool samples. The overall prevalence of H. pylori infection in our asymptomatic study population was 64%. Age, gender, occupation and frequency of brushing teeth were not found to be associated with H. pylori colonization. The vacA gene was successfully sequenced from both saliva and stool samples of 12 individuals. For seven of these individuals, saliva and stool sequences fell into different clusters on a phylogenetic tree, indicating intra-host genetic variation of H. pylori. This study reports a high prevalence of H. pylori infection in asymptomatic persons in this region of Thailand and demonstrates that genotypes (vacA gene sequences) of H. pylori may differ between the oral cavity and intestinal tract.

Tuberculous meningitis (TBM) is a severe form of tuberculosis with a high mortality rate. The factors associated with TBM pathogenesis are still unclear. Using comparative whole-genome sequence analysis we compared Mycobacterium tuberculosis (Mtb) isolates from cerebrospinal fluid of TBM cases (n = 73) with those from sputum of pulmonary tuberculosis (PulTB) patients (n = 220) from Thailand. The aim of this study was to seek genetic variants of Mtb associated with TBM. Regardless of Mtb lineage, we found 242 variants that were common to all TBM isolates. Among these variants, 28 were missense SNPs occurring mainly in the pks genes (involving polyketide synthesis) and the PE/PPE gene. Six lineage-independent SNPs were commonly found in TBM isolates, two of which were missense SNPs in Rv0532 (PE_PGRS6). Structural variant analysis revealed that PulTB isolates had 14 genomic regions containing 2–3-fold greater read depth, indicating higher copy number variants and half of these genes belonged to the PE/PPE gene family. Phylogenetic analysis revealed only two small clusters of TBM clonal isolates without support from epidemiological data. This study reported genetic variants of Mtb commonly found in TBM patients compared to PulTB patients. Variants associated with TBM disease warrant further investigation.