Discordance Between Phenotypic and WGS-Based Drug Susceptibility Testing Results for Some Anti-Tuberculosis Drugs: A Snapshot Study of Paired Mycobacterium tuberculosis Isolates with Small Genetic Distance.

Background: Current tuberculosis treatment regimens primarily rely on phenotypic drug susceptibility testing and rapid molecular assays. Although whole-genome sequencing (WGS) offers a promising alternative, disagreements between phenotypic and molecular testing methods remain. In this retrospective study, we compared the phenotypic and WGS-predicted drug resistance profiles of paired Mycobacterium tuberculosis isolates with small genetic distances (≤10 single nucleotide variants) obtained from patients with longitudinal single-episode or recurrent tuberculosis. Additionally, we investigated the distribution of drug-resistance-conferring variants among the identified M. tuberculosis genotypes. Methods: Paired M. tuberculosis isolates from 46 patients with pulmonary tuberculosis (2002-2019) were analyzed. Spoligotyping was performed for all the isolates. WGS data were processed using TB-Profiler software to genotype the strains and detect variants in M. tuberculosis genes associated with drug resistance. The significance of these variants was evaluated using the M. tuberculosis variant catalog developed by the World Health Organization. Phenotypic drug susceptibility test results were obtained from patients' medical records. Results: Among the 46 isolate pairs, 25 (54.3%) harbored drug-resistance-associated variants, with 20 demonstrating identical WGS-predicted drug resistance profiles. Drug-resistant isolate pairs belonged to Lineages 2 and 4, with the most common sub-lineages being 2.2.1 (SIT1 and SIT190 spoligotypes), and 4.3.3 (SIT42). Agreement between phenotypic and WGS-based drug susceptibility testing was highest (>90%) for rifampicin, isoniazid, ethambutol, fluoroquinolones, streptomycin, and amikacin when calculated for M. tuberculosis isolates or isolate pairs. In most discordant cases, isolate pairs harbored variants that could cause low- or moderate-level resistance or were previously associated with variable minimum inhibitory concentrations. Notably, such discrepancies mostly occurred in one isolate from the pair. In addition, differences in resistance-related variant distributions among M. tuberculosis genotypes were observed for most of the analyzed drugs. Conclusion: The simultaneous performance of phenotypic and WGS-based drug susceptibility testing creates the most accurate drug resistance profile for M. tuberculosis isolates and eliminates important limitations of each method.

Advantages of analysing both pairwise SNV-distance and differing SNVs between Mycobacterium tuberculosis isolates for recurrent tuberculosis cause determination.

Endogenous reactivation and exogenous reinfection are two possible causes of recurrent tuberculosis (TB). However, in some cases, precise cause determination can be challenging. In this study, we used whole genome sequencing to determine pairwise SNV distances and detect differing SNVs in initial and subsequent isolates for recurrent TB cases when the first and second episodes were caused by Mycobacterium tuberculosis (Mtb) strains with an identical spoligotype pattern. In total, 104 Mtb isolates from 36 recurrent TB and 16 single TB episode patients were included in the study. Most isolate pairs belonged to the SIT1 (n=21), SIT42 (n=9), SIT53 (n=9), and SIT254 (n=7) spoligotypes, and in 27 cases, resistance to at least one anti-TB drug was found in either isolate. Drug susceptibility was more common in the recurrent TB patient cohort, and longitudinal single TB episode isolates were more prone to be drug-resistant (p=0.03), while the association between patient cohort and spoligotype was not statistically significant (p=0.07). The pairwise SNV-distance between the longitudinal single TB episode isolates was small (0-7 SNVs). Among the recurrent TB isolates, based on the high SNV-distance (38-273 SNVs), six reinfection cases (16.7%) were identified. This distance was small (<10 SNVs) in the remaining 30 isolate pairs. Further analysis of differing SNVs revealed that 22 (61.1%) cases could be classified as possible reactivation. Notably, despite the small distance of 2-7 SNVs, initial isolates of eight patients (22.2%) had several SNVs that were not found in the second isolates; therefore, these cases were classified as reinfection with a closely related Mtb strain. No statistically significant difference in the time interval between specimen collection in the reactivation and reinfection Mtb sample groups (p=0.13) or an association between recurrence cause and drug resistance status (p=0.62) or spoligotype (p=0.79) could be detected. The mycobacterial median mutation rate of longitudinal single TB episodes and possible reactivation isolate pairs (n=37) was 0.12 SNVs/genome/year (IQR 0-0.39), and in 18 cases (48.6%), it was equal to zero. No statistically significant differences in mutation rate were found between recurrent TB and longitudinal single TB episode isolates (p=0.087), drug-susceptible and resistant isolates (p=0.37) or isolates of Beijing and other genotype families (p=0.33). Furthermore, four cases of fluoroquinolone resistance development through the acquired SNVs in the gyrA gene were identified. To conclude, this study highlighted the complexity of recurrent episode cause determination and showed the usefulness of differing SNV identification in both Mtb isolates in such cases. Expected drug susceptibility was the only discriminative factor for recurrent TB episode-causing mycobacterial strains, while no differences between reactivation and reinfection sample groups could be identified.

Analysis of Mycobacterium tuberculosis genetic lineages circulating in Riga and Riga region, Latvia, isolated between 2008 and 2012.

Although the number of new tuberculosis (TB) cases registered per year has decreased by 3-fold between 2001 and 2017 in Latvia, the TB incidence and rates of multidrug resistant TB in this Baltic country remain substantially higher than in most other European countries. Molecular typing methods of Mycobacterium tuberculosis (MTB) play an important role both in clinical studies of the disease and the epidemiological investigations, allowing to describe and characterize the pathogen's population structure and spread of particular genotypes. Aim of this study was to examine the prevalence of MTB lineages in Riga and Riga region of Latvia within a five-year period (2008-2012), and to evaluate the discriminatory power (DP) of spoligotyping, standard 24-locus MIRU-VNTR and IS6110-RFLP methods in this setting. The results showed that the main MTB spoligotype families were Beijing (25.3%) and LAM (24.3%), followed by T (22.1%), Ural (11.2%), Haarlem (6.6%) and X superfamily (3.4%). This distribution remained stable over the five consecutive years. 67.6% of MTB isolates were pan-susceptible, and 32.4% were resistant to any drug; multi-drug resistance was found in 5.8% of MTB strains, and 7.6% of MTB isolates were extensively drug-resistant. Drug resistance was associated with SIT1, SIT283 and SIT42 genotypes, while SIT1 and SIT42 were overrepresented among multi drug-resistant MTB strains. Overall, DP of spoligotyping method alone was 0.8953, while DP of both 24-locus MIRU-VNTR analysis and IS6110 RFLP was higher (DP = 0.9846 and 0.9927, respectively), mainly due to the improvement of the resolution for the Beijing strains. In conclusion, this work represents the first comprehensive molecular epidemiological description of TB in Latvia, highlighting the high genetic diversity of MTB strains circulating in Riga and Riga region. In combination with detailed epidemiological data this approach was helpful for the in-depth understanding of epidemiological processes in settings where the Next-Gen sequencing is not available as a routine method.

Genotypic and phenotypic characteristics of aminoglycoside-resistant Mycobacterium tuberculosis isolates in Latvia.

Mutations causing resistance to aminoglycosides, such as kanamycin (KAN), amikacin (AMK), and streptomycin, are not completely understood. In this study, polymorphisms of aminoglycoside resistance influencing genes such as rrs, eis, rpsL, and gidB in 41 drug-resistant and 17 pan-sensitive Mycobacterium tuberculosis clinical isolates in Latvia were analyzed. Mutation A1400G in rrs gene was detected in 92% isolates with high resistance level to KAN and diverse MIC level to AMK. Mutations in promoter region of eis were detected in 80% isolates with low-level MIC of KAN. The association of K43R mutation in rpsL gene, a mutation in the rrs gene at position 513, and various polymorphisms in gidB gene with distinct genetic lineages of M. tuberculosis was observed. The results of this study suggest that association of different controversial mutations of M. tuberculosis genes to the drug resistance phenotype should be done in respect to genetic lineages.

Characterisation of rpsL, rrs and embB mutations associated with streptomycin and ethambutol resistance in Mycobacterium tuberculosis.

In order to characterise molecular mechanisms of first-line drug resistance in Mycobacterium tuberculosis and to evaluate the use of molecular markers of resistance (gene point mutations), we analysed 66 multi-drug-resistant (MDR) isolates from Latvian tuberculosis patients. They were all resistant to rifampin (RIF), isoniazid (INH) and streptomycin (SM), and 33 were resistant to ethambutol (EMB). Enzymatic digestion by MboII and nucleotide sequencing of the rpsL gene fragment detected a single nucleotide substitution K43R in 40 (61%) of the 66 SM-resistant M. tuberculosis isolates. Of the other 26 SM-resistant isolates, 16 (24%) had mutations at positions 513A-->C and 516C-->T of the rrs gene and 10 (15%) had the wild-type sequence. The single-stranded DNA conformation polymorphism (SSCP) method was used to detect mutations in the embB gene associated with EMB resistance. Substitutions in the embB gene were found by SSCP analysis in 15 (45%) and by sequencing in 17 (52%) of the 33 EMB-resistant isolates. Surprisingly, SSCP revealed a nucleotide mutation at codon M306 in five (15%) of 33 in vitro EMB-susceptible MDR isolates.