Phenotypic drug-susceptibility profiles and genetic analysis based on whole-genome sequencing of Mycobacterium avium complex isolates in Thailand.

Mycobacterium avium complex (MAC) infections are a significant clinical challenge. Determining drug-susceptibility profiles and the genetic basis of drug resistance is crucial for guiding effective treatment strategies. This study aimed to determine the drug-susceptibility profiles of MAC clinical isolates and to investigate the genetic basis conferring drug resistance using whole-genome sequencing (WGS) analysis. Drug-susceptibility profiles based on minimum inhibitory concentration (MIC) assays were determined for 38 MAC clinical isolates (12 Mycobacterium avium and 26 Mycobacterium intracellulare). Mutations associated with drug resistance were identified through genome analysis of these isolates, and their phylogenetic relationships were also examined. Drug resistance, based on MIC values, was most commonly observed for moxifloxacin (81.6%), followed by linezolid (78.9%), clarithromycin (44.7%) and amikacin (36.8%). We identified specific mutations associated with resistance to amikacin. These include the rrs mutation at C464T in amikacin intermediate-resistance M. avium, and two mutations at T250A and G1453T in amikacin non-susceptible M. intracellulare. Mutations in rrl at A2058G, A2059C and A2059G were potentially linked to clarithromycin resistance. MAC clinical isolates not susceptible to linezolid exhibited mutations in rplC at G237C and C459T, as well as two rplD mutations at G443A and A489G. GyrB substitution Thr521Ala (T521A) was identified in moxifloxacin non-susceptible isolates, which may contribute to this resistance. A phylogeny of our MAC isolates revealed high levels of genetic diversity. Our findings suggest that the standard treatment regimen for MAC infections using moxifloxacin, linezolid, clarithromycin and amikacin may be driving development of resistance, potentially due to specific mutations. The combination of phenotypic and genotypic susceptibility testing can be valuable in guiding the clinical use of drugs for the treatment of MAC infections.

Whole-Genome Sequencing and Drug-Susceptibility Analysis of Serial Mycobacterium abscessus Isolates from Thai Patients.

Mycobacterium abscessus is an important pathogen that can cause serious human diseases and is difficult to treat due to antibiotic resistance. In this study, we analyzed, using whole-genome sequence (WGS) data, M. abscessus strains serially isolated from patients at various time intervals. We undertook genetic diversity analysis between subspecies, mutation-rate estimation and identification of drug-resistant mutations with minimum inhibitory concentration (MIC) analysis. Clonal isolates of M. abscessus:­subsp. abscessus (MAB) and subsp. massiliense (MMAS)­causing persistent infection through time, differed by 0−7 and 0−14 SNPs, respectively, despite being isolated 1 to 659 days apart. Two cases caused by MMAS differed by ≥102 SNPs at 350 days apart and were regarded as examples of reinfection. Isolates collected ≤7 days apart exhibited a high mutation rate (133.83 ± 0.00 SNPs/genome (5 Mb)/year for MMAS and 127.75 SNPs/genome (5 Mb)/year for MAB). Mutation rates declined in a time-dependent manner in both subspecies. Based on isolates collected > 180 days apart, MMAS had a significantly higher average mutation rate than MAB (2.89 ± 1.02 versus 0.82 ± 0.83 SNPs/genome (5 Mb)/year, (p = 0.01), respectively). All well-known drug-resistance mutations were found to be strongly associated with high MIC levels for clarithromycin and ciprofloxacin. No known mutations were identified for strains resistant to linezolid and amikacin. MAB strains in the study were susceptible to amikacin, while most MMAS strains were susceptible to clarithromycin, amikacin and linezolid. No hetero-resistance was found in the strains analyzed. Our study reports the genetic diversity and mutation rate of M. abscessus between the two major subspecies and confirms the drug resistance-associated mutations. Information about drug-resistance and associated mutations can be applied in diagnosis and patient management.

High clustering rate and genotypic drug-susceptibility screening for the newly recommended anti-tuberculosis drugs among global extensively drug-resistant Mycobacterium tuberculosis isolates.

Multidrug-resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB) make TB difficult to control. Global susceptibility data for six newly recommended anti-TB drugs against M/XDR-TB are still limited. Using publicly available whole-genome sequences, we determined the proportion of 513 phenotypically XDR-TB isolates that carried mutations associated with resistance against these drugs (bedaquiline, clofazimine, linezolid, delamanid, pretomanid and cycloserine). Mutations of Rv0678 and Rv1979c were detected in 69/513 isolates (13.5%) for bedaquiline resistance and 79/513 isolates (15.4%) for clofazimine resistance with additional mmpL5 mutations. Mutations conferring resistance to delamanid were detected in fbiB and ddn genes for 11/513 isolates (2.1%). For pretomanid, a mutation was detected in the ddn gene for 3/513 isolates (0.6%). Nineteen mutations of pykA, cycA, ald, and alr genes, conferring resistance to cycloserine, were found in 153/513 isolates (29.8%). No known mutations associated with linezolid resistance were detected. Cluster analysis showed that 408/513 isolates fell within 99 clusters and that 354 of these isolates were possible primary drug-resistant TB (292 XDR-TB, 57 pre-XDR-TB and 5 MDR-TB). Clonal transmission of primary XDR isolates might contribute significantly to the high prevalence of DR-TB globally.

Optimized analysis parameters of variant calling for whole genome-based phylogeny of Mycobacteroides abscessus.

Whole-genome sequence (WGS) analysis provides the best resolution for reconstructing bacterial phylogeny. However, the resulting tree could vary according to parameters used in the WGS pipeline, making it difficult to compare results across multiple studies. This study compares effects on phylogenies when applying different parameter stringencies. We used as the study model to optimize parameters strains of Mycobacteroides abscessus serially isolated at various intervals, isolates known to represent persistent infection (PI) cases or re-infection (RI) cases and isolates from different subspecies. Un-optimized parameters with low stringency provided an excessive number of SNPs (823) compared to the optimized setting (3 SNPs) between paired strains isolated 1 day apart from PI cases, discordant tree topology and misclassification of subspecies and of instances of RI. We demonstrated that using high-quality variants provides more accuracy for recognizing serial isolates of the same clone versus different clones and for phylogenetic analysis of M. abscessus. Our approach might be used as a model for analyses requiring phylogenetic reconstruction of other bacteria.

High rate of reinfection and possible transmission of Mycobacterium avium complex in Northeast Thailand.

The Mycobacterium avium complex (MAC) includes two main species of non-tuberculous mycobacteria (NTM), M. avium and Mycobacterium intracellulare. These can cause serious disease, especially in immunocompromised patients. Little information is available concerning genetic diversity of NTM. We used multilocus sequence typing (MLST) based on a highly discriminative gene set to analyze MAC serially isolated from patients to determine the rate of MAC reinfection. Genomic DNA was sequenced from 49 MAC isolates (15 cases comprised of 11 true infections and 4 instances of colonization). More than half of the MAC isolates tested were found to be multidrug resistant. The discriminatory power was assessed of 24 house-keeping genes (fusA, atpD, pheT, glnA, topA, secA, argH, glpK, murC, cya, pta, rrl, rrs, hsp65, rpoB, 16S-23S rRNA ITS, recF, lipT, pepB, gnd, aspB, groEL, sodA and est) previously used for genotyping of MAC and other NTM. Seven genes (fusA, secA, rpoB, hsp65, 16S rRNA, 23S rRNA, 16S-23S rRNA ITS) had a discriminatory power index higher than 0.9 and were included in the optimized set that we used. This set was significantly better for genotyping and diagnosis of MAC than previously used 4-gene, 5-gene and 9-gene sets. MLST using our 7-gene set indicated that the rate of reinfection was 54.55% (6/11 cases). Persistent infections (n = 5 cases, 45.45%) were found. A changing of clone in the same patient was found in 1/4 (25%) of the colonization cases. Two small clusters of possible MAC transmission between humans were found. Our study demonstrated that the high frequency of apparent treatment failure of MAC might be artefactual, as a consequence of a high rate of MAC reinfection in Thai population. Our useful highly discriminative gene set for MAC species and clonal strain analysis could be further applied for the diagnosis and patient management.

Assessment of antimycobacterial activities of pure compounds extracted from Thai medicinal plants against clarithromycin-resistant Mycobacterium abscessus.

BACKGROUND: Infection with Mycobacterium abscessus is usually chronic and is associated with clarithromycin resistance. Increasing drug resistance is a major public-health problem and has led to the search for new antimycobacterial agents. We evaluated the antimycobacterial activity, toxicity, and synergistic effects of several plant secondary metabolites against M. abscessus. METHODS: Twenty-three compounds were evaluated for antimycobacterial activity against thirty M. abscessus clinical isolates by broth microdilution to determine their minimum inhibitory concentration (MIC) values. Toxicity was evaluated using red and white blood cells (RBCs and WBCs). The compounds were used in combination with clarithromycin to investigate the possibility of synergistic activity. RESULTS: Five out of twenty-three compounds (RL008, RL009, RL011, RL012 and RL013) exhibited interesting antimycobacterial activity against M. abscessus, with MIC values ranging from <1 to >128 µg/mL. These extracts did not induce hemolytic effect on RBCs and displayed low toxicity against WBCs. The five least-toxic compounds were tested for synergism with clarithromycin against seven isolates with inducible clarithromycin resistance and seven with acquired clarithromycin resistance. The best synergistic results against these isolates were observed for RL008 and RL009 (8/14 isolates; 57%). CONCLUSIONS: This study demonstrated antimycobacterial and synergistic activities of pure compounds extracted from medicinal plants against clarithromycin-resistant M. abscessus. This synergistic action, together with clarithromycin, may be effective for treating infections and should be further studied for the development of novel antimicrobial agents.

Assessment of in vitro activities of novel modified antimicrobial peptides against clarithromycin resistant Mycobacterium abscessus.

Mycobacterium abscessus (Mab) is one of the most drug resistant bacteria with a high treatment failure rate. Antimicrobial peptides (AMPs) are alternative therapeutic agents against this infection. This study was aimed to assess the in vitro activities of thirteen AMPs (S5, S52, S6, S61, S62, S63, KLK, KLK1, KLK2, Pug-1, Pug-2, Pug-3 and Pug-4) that have never been investigated against drug resistant Mab isolates. Only four novel modified AMPs (S61, S62, S63 and KLK1) provided the lowest minimum inhibitory concentration (MIC) values ranging from 200-400 µg/ml against the Mab ATCC19977 strain. These four potential AMPs were further tested with 16 clinical isolates of clarithromycin resistant Mab. The majority of the tested strains (10/16 isolates, 62.5%) showed ~99% kill by all four AMPs within 24 hours with an MIC <50 µg/ml. Only two isolates (12.5%) with acquired clarithromycin resistance, however, exhibited values <50 µg/ml of four potential AMPs, S61, S62, S63 and KLK1 after 3-days-incubation. At the MICs level, S63 showed the lowest toxicity with 1.50% hemolysis and 100% PBMC viability whereas KLK1 showed the highest hemolysis (10.21%) and lowest PBMC viability (93.52%). S61, S62 and S63 were further tested with clarithromycin-AMP interaction assays and found that 5/10 (50%) of selected isolates exhibited a synergistic interaction with 0.02-0.41 FICI values. This present study demonstrated the potential application of novel AMPs as an adjunctive treatment with clarithromycin against drug resistant Mab infection.

Whole-Genome Sequencing Analysis to Identify Infection with Multiple Species of Nontuberculous Mycobacteria.

Mixed infection with multiple species of nontuberculous mycobacteria (NTM) is difficult to identify and to treat. Current conventional molecular-based methods for identifying mixed infections are limited due to low specificity. Here, we evaluated the utility of whole-genome sequencing (WGS) analysis to detect and identify mixed NTM infections. Analytical tools used included PubMLST, MetaPhlAn3, Kraken2, Mykrobe-Predictor and analysis of heterozygous SNP frequencies. The ability of each to identify mixed infections of NTM species was compared. Sensitivity was tested using 101 samples (sequence sets) including 100 in-silico simulated mixed samples with various proportions of known NTM species and one sample of known mixed NTM species from a public database. Single-species NTM control samples (155 WGS samples from public databases and 15 samples from simulated reads) were tested for specificity. Kraken2 exhibited 100% sensitivity and 98.23% specificity for detection and identification of mixed NTM species with accurate estimation of relative abundance of each species in the mixture. PubMLST (99% and 96.47%) and MetaPhlAn3 (95.04% and 83.52%) had slightly lower sensitivity and specificity. Mykrobe-Predictor had the lowest sensitivity (57.42%). Analysis of read frequencies supporting single nucleotide polymorphisms (SNPs) could not detect mixed NTM samples. Clinical NTM samples (n = 16), suspected on the basis of a 16S-23S rRNA gene sequence-based line-probe assay (LPA) to contain more than one NTM species, were investigated using WGS-analysis tools. This identified only a small proportion (37.5%, 6/16 samples) of the samples as mixed infections and exhibited only partial agreement with LPA results. LPAs seem to be inadequate for detecting mixed NTM species infection. This study demonstrated that WGS-analysis tools can be used for diagnosis of mixed infections with different species of NTM.

Whole-genome sequence analysis and comparisons between drug-resistance mutations and minimum inhibitory concentrations of Mycobacterium tuberculosis isolates causing M/XDR-TB.

Drug resistance (DR) remains a major challenge for tuberculosis (TB) control. Whole-genome sequencing (WGS) provides the highest genetic resolution for genotypic drug-susceptibility tests (DST). We compared DST profiles of 60 Mycobacterium tuberculosis isolates which were drug resistant according to agar proportion tests (one poly DR-TB, 34 multidrug-resistant TB and 25 extensively drug-resistant TB). We additionally performed minimum inhibitory concentration (MIC) tests using Sensititre MYCOTBI plates (MYCOTB) and a WGS-based DST. Agreement between WGS-based DST and MYCOTB was high for all drugs except ethambutol (65%) and ethionamide (62%). Isolates harboring the -15 c/t inhA promoter mutation had a significantly lower MIC for isoniazid than did isolates with the katG Ser315Thr mutation (p < 0.001). Similar patterns were seen for ethambutol (embB Gly406Asp vs. embB Met306Ile), streptomycin (gid Gly73Ala vs. rpsL Lys43Arg), moxifloxacin (gyrA Ala90Val vs. gyrA Asp94Gly) and rifabutin (rpoB Asp435Phe/Tyr/Val vs. rpoB Ser450Leu). For genotypic heteroresistance, isolates with lower proportion of mapped read tended to has lower MIC of anti-TB drugs than those with higher proportion. These results emphasize the high applicability of WGS for determination of DR-TB and the association of particular mutations with MIC levels.

Whole-genome sequencing analysis of multidrug-resistant Mycobacterium tuberculosis from Java, Indonesia.

Introduction. Multidrug-resistant tuberculosis (MDR-TB) is a major public health problem globally, including in Indonesia. Whole-genome sequencing (WGS) analysis has rarely been used for the study of TB and MDR-TB in Indonesia.Aim. We evaluated the use of WGS for drug-susceptibility testing (DST) and to investigate the population structure of drug-resistant Mycobacterium tuberculosis in Java, Indonesia.Methodology. Thirty suspected MDR-TB isolates were subjected to MGIT 960 system (MGIT)-based DST and to WGS. Phylogenetic analysis was done using the WGS data. Results obtained using MGIT-based DST and WGS-based DST were compared.Results. Agreement between WGS and MGIT was 93.33 % for rifampicin, 83.33 % for isoniazid and 76.67 % for streptomycin but only 63.33 % for ethambutol. Moderate WGS-MGIT agreement was found for second-line drugs including amikacin, kanamycin and fluoroquinolone (73.33-76.67 %). MDR-TB was more common in isolates of the East Asian Lineage (63.3%). No evidence of clonal transmission of DR-TB was found among members of the tested population.Conclusion. Our study demonstrated the applicability of WGS for DST and molecular epidemiology of DR-TB in Java, Indonesia. We found no transmission of DR-TB in Indonesia.

Comparisons of whole-genome sequencing and phenotypic drug susceptibility testing for Mycobacterium tuberculosis causing MDR-TB and XDR-TB in Thailand.

Drug-resistant tuberculosis (TB) is a major public health problem. There is little information regarding the genotypic-phenotypic association of anti-TB drugs, especially for second-line drugs. This study compared phenotypic drug susceptibility testing (DST) with predictions based on whole-genome sequencing (WGS) data for 266 Mycobacterium tuberculosis isolates. Phenotypic DST used the standard proportional method. Clinical isolates of M. tuberculosis collected in Thailand between 1998 and 2013 comprised 51 drug-sensitive strains, six mono-resistant strains, two multiple-resistant strains, 88 multi-drug-resistant strains, 95 pre-extensively drug-resistant strains and 24 extensively drug-resistant strains. WGS analysis was performed using the computer programs PhyResSE and TB-Profiler. TB-Profiler had higher average concordance with phenotypic DST than PhyResSE for both first-line (91.96% vs. 91.4%) and second-line (79.67% vs. 78.20%) anti-TB drugs. The average sensitivity for all anti-TB drugs was also higher (83.13% vs. 72.08%) with slightly lower specificity (83.50% vs. 86.68%). Regardless of the program used, isoniazid, rifampicin and amikacin had the highest concordance with phenotypic DST (96.2%, 93.5% and 95.6%, respectively). Ethambutol, ethionamide and fluoroquinolones had the lowest concordance (87.34%, 81.44% and 73.85%, respectively). Concordance rates of ofloxacin (a second-generation fluoroquinolone), levofloxacin, moxifloxacin and gatifloxacin (third- and fourth-generation fluoroquinolones) were 91.79%, 76.62%, 72.64% and 57.35%, respectively. Discordance between phenotypic and WGS-based DSTs may be due, in part, to the choice of critical concentration and variable reproducibility of the phenotypic tests. It may also be due to limitations of the mutation databases (especially for the second-line drugs) and the analysis program used. Mutations related to fluoroquinolone resistance, especially the later generations, need to be identified.