Unraveling tuberculosis patient cluster transmission chains: integrating WGS-based network with clinical and epidemiological insights.

Background: Tuberculosis remains a global health threat, and the World Health Organization reports a limited reduction in disease incidence rates, including both new and relapse cases. Therefore, studies targeting tuberculosis transmission chains and recurrent episodes are crucial for developing the most effective control measures. Herein, multiple tuberculosis clusters were retrospectively investigated by integrating patients' epidemiological and clinical information with median-joining networks recreated based on whole genome sequencing (WGS) data of Mycobacterium tuberculosis isolates. Methods: Epidemiologically linked tuberculosis patient clusters were identified during the source case investigation for pediatric tuberculosis patients. Only M. tuberculosis isolate DNA samples with previously determined spoligotypes identical within clusters were subjected to WGS and further median-joining network recreation. Relevant clinical and epidemiological data were obtained from patient medical records. Results: We investigated 18 clusters comprising 100 active tuberculosis patients 29 of whom were children at the time of diagnosis; nine patients experienced recurrent episodes. M. tuberculosis isolates of studied clusters belonged to Lineages 2 (sub-lineage 2.2.1) and 4 (sub-lineages 4.3.3, 4.1.2.1, 4.8, and 4.2.1), while sub-lineage 4.3.3 (LAM) was the most abundant. Isolates of six clusters were drug-resistant. Within clusters, the maximum genetic distance between closely related isolates was only 5-11 single nucleotide variants (SNVs). Recreated median-joining networks, integrated with patients' diagnoses, specimen collection dates, sputum smear microscopy, and epidemiological investigation results indicated transmission directions within clusters and long periods of latent infection. It also facilitated the identification of potential infection sources for pediatric patients and recurrent active tuberculosis episodes refuting the reactivation possibility despite the small genetic distance of ≤5 SNVs between isolates. However, unidentified active tuberculosis cases within the cluster, the variable mycobacterial mutation rate in dormant and active states, and low M. tuberculosis genetic variability inferred precise transmission chain delineation. In some cases, heterozygous SNVs with an allelic frequency of 10-73% proved valuable in identifying direct transmission events. Conclusion: The complex approach of integrating tuberculosis cluster WGS-data-based median-joining networks with relevant epidemiological and clinical data proved valuable in delineating epidemiologically linked patient transmission chains and deciphering causes of recurrent tuberculosis episodes within clusters.

Effect of NAT2, GSTM1 and CYP2E1 genetic polymorphisms on plasma concentration of isoniazid and its metabolites in patients with tuberculosis, and the assessment of exposure-response relationships.

Objectives: Isoniazid is a key drug in the chemotherapy of tuberculosis (TB), however, interindividual variability in pharmacokinetic parameters and drug plasma levels may affect drug responses including drug induced hepatotoxicity. The current study investigated the relationships between isoniazid exposure and isoniazid metabolism-related genetic factors in the context of occurrence of drug induced hepatotoxicity and TB treatment outcomes. Methods: Demographic characteristics and clinical information were collected in a prospective TB cohort study in Latvia (N = 34). Time to sputum culture conversion (tSCC) was used as a treatment response marker. Blood plasma concentrations of isoniazid (INH) and its metabolites acetylisoniazid (AcINH) and isonicotinic acid (INA) were determined at three time points (pre-dose (0 h), 2 h and 6 h after drug intake) using liquid chromatography-tandem mass spectrometry. Genetic variations of three key INH-metabolizing enzymes (NAT2, CYP2E1, and GSTM1) were investigated by application PCR- and Next-generation sequencing-based methods. Depending on variables, group comparisons were performed by Student's t-test, one-way ANOVA, Mann-Whitney-Wilcoxon, and Kruskal-Wallis tests. Pearson correlation coefficient was calculated for the pairs of normally distributed variables; model with rank transformations were used for non-normally distributed variables. Time-to-event analysis was performed to analyze the tSCC data. The cumulative probability of tSCC was obtained using Kaplan-Meier estimators. Cox proportional hazards models were fitted to estimate hazard rate ratios of successful tSCC. Results: High TB treatment success rate (94.1%) was achieved despite the variability in INH exposure. Clinical and demographic factors were not associated with either tSCC, hepatotoxicity, or INH pharmacokinetics parameters. Correlations between plasma concentrations of INH and its metabolites were NAT2 phenotype-dependent, while GSTM1 genetic variants did not showed any effects. CYP2E1*6 (T > A) allelic variant was associated with INH pharmacokinetic parameters. Decreased level of AcINH was associated with hepatotoxicity, while decreased values of INA/INH and AcINH/INH were associated with month two sputum culture positivity. Conclusion: Our findings suggest that CYP2E1, but not GSTM1, significantly affects the INH pharmacokinetics along with NAT2. AcINH plasma level could serve as a biomarker for INH-related hepatotoxicity, and the inclusion of INH metabolite screening in TB therapeutic drug monitoring could be beneficial in clinical studies for determination of optimal dosing strategies.

Exploring Variability in Rifampicin Plasma Exposure and Development of Anti-Tuberculosis Drug-Induced Liver Injury among Patients with Pulmonary Tuberculosis from the Pharmacogenetic Perspective.

Genetic polymorphisms can exert a considerable impact on drug pharmacokinetics (PK) and the development of adverse drug reactions (ADR). However, the effect of genetic polymorphisms on the anti-tuberculosis (anti-TB) drug, and particularly rifampicin (RIF), exposure or anti-TB drug-induced liver injury (DILI) remains uncertain. Here, we evaluated the relationship between single nucleotide polymorphisms (SNPs) detected in the RIF pharmacogenes (AADAC, SLCO1B1, SLCO1B3, ABCB1, and NR1I2) and RIF PK parameters, as well as anti-TB treatment-associated DILI. In total, the study enrolled 46 patients with drug-susceptible pulmonary TB. The RIF plasma concentration was measured using the LC-MS/MS method in the blood samples collected pre-dose and 2 and 6 h post-dose, whilst the DILI status was established using the results from blood biochemical analysis performed before and 10-12 days after treatment onset. The genotyping was conducted using a targeted NGS approach. After adjustment for confounders, the patients carrying the rs3732357 GA/AA genotype of the NR1I2 gene were found to have significantly lower RIF plasma AUC0-6 h in comparison to those with GG genotype, while the difference in RIF plasma Cmax was insignificant. None of the analyzed SNPs was related to DILI. Hence, we are the first to report NR1I2 intronic SNP rs3732357 as the genetic component of variability in RIF exposure. Regarding anti-TB treatment-associated DILI, the other preexisting factors promoting this ADR should be considered.

Genotypic and phenotypic comparison of drug resistance profiles of clinical multidrug-resistant Mycobacterium tuberculosis isolates using whole genome sequencing in Latvia.

BACKGROUND: Multidrug-resistant tuberculosis (MDR-TB) remains a major public health problem in many high tuberculosis (TB) burden countries. Phenotypic drug susceptibility testing (DST) take several weeks or months to result, but line probe assays and Xpert/Rif Ultra assay detect a limited number of resistance conferring gene mutations. Whole genome sequencing (WGS) is an advanced molecular testing method which theoretically can predict the resistance of M. tuberculosis (Mtb) isolates to all anti-TB agents through a single analysis. METHODS: Here, we aimed to identify the level of concordance between the phenotypic and WGS-based genotypic drug susceptibility (DS) patterns of MDR-TB isolates. Overall, data for 12 anti-TB medications were analyzed. RESULTS: In total, 63 MDR-TB Mtb isolates were included in the analysis, representing 27.4% of the total number of MDR-TB cases in Latvia in 2012-2014. Among them, five different sublineages were detected, and 2.2.1 (Beijing group) and 4.3.3 (Latin American-Mediterranean group) were the most abundant. There were 100% agreement between phenotypic and genotypic DS pattern for isoniazid, rifampicin, and linezolid. High concordance rate (> 90%) between phenotypic and genotypic DST results was detected for ofloxacin (93.7%), pyrazinamide (93.7%) and streptomycin (95.4%). Phenotypic and genotypic DS patterns were poorly correlated for ethionamide (agreement 56.4%), ethambutol (85.7%), amikacin (82.5%), capreomycin (81.0%), kanamycin (85.4%), and moxifloxacin (77.8%). For capreomycin, resistance conferring mutations were not identified in several phenotypically resistant isolates, and, in contrary, for ethionamide, ethambutol, amikacin, kanamycin, and moxifloxacin the resistance-related mutations were identified in several phenotypically sensitive isolates. CONCLUSIONS: WGS is a valuable tool for rapid genotypic DST for all anti-TB agents. For isoniazid and rifampicin phenotypic DST potentially can be replaced by genotypic DST based on 100% agreement between the tests. However, discrepant results for other anti-TB agents limit their prescription based solely on WGS data. For clinical decision, at the current level of knowledge, there is a need for combination of genotypic DST with modern, validated phenotypic DST methodologies for those medications which did not showed 100% agreement between the methods.

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.

Next-Generation Sequencing and Bioinformatics-Based Protocol for the Full-Length CYP2E1 Gene Polymorphism Analysis.

Introduction: Pharmacogenetics studies provide clinically relevant information on the identified associations between genetic variants and individual variability in drug response, which, in turn, offers great promise for guiding personalized drug therapy and clinical trial design. However, there is a lack of information concerning the evidence-based clinical annotations of specific CYP2E1 genetic variants. Aim: To design and evaluate the next-generation sequencing-based method for full-length CYP2E1 gene polymorphism analysis. Materials and Methods: Seven gene-specific oligonucleotide primer pairs targeting overlapping CYP2E1 gene fragments spanning all nine gene exons with interleaving introns, untranslated (UTR) and intergenic regions were designed. Human DNA samples (n = 3) were used as a training set to check the primer performance and to optimize the PCR conditions. The effectiveness of the developed target amplification and sequencing protocol was evaluated using the test set comprising human DNA samples (n = 3) obtained from tuberculosis patients. Sequencing data analysis was performed on the Galaxy online-based platform. Results: The sequencing data quality was sufficient for the detection of genetic variants dispersed throughout the CYP2E1 gene with a high degree of confidence in fully covered regions achieving optimal reading depth of the targeted fragment with high base call accuracy. Conclusion: Developed protocol can be applied in subpopulation-level association studies to determine whether single nucleotide variants (SNVs) or variant combinations from multiple regions of the CYP2E1 gene are of clinical significance.

LC-MS/MS method for simultaneous quantification of the first-line anti-tuberculosis drugs and six primary metabolites in patient plasma: Implications for therapeutic drug monitoring.

The pharmacokinetic profiling of drug substances and corresponding metabolites in the biological matrix is one of the most informative tools for the treatment efficacy assessment. Therefore, to satisfy the need for comprehensive monitoring of anti-tuberculosis drugs in human plasma, a liquid chromatography-tandem mass spectrometry (LC-MS/MS) method was developed and validated for simultaneous quantification of first-line anti-tuberculosis drugs (ethambutol, isoniazid, pyrazinamide, and rifampicin) along with their six primary metabolites. Simple single-step protein precipitation with methanol was chosen as the most convenient sample pre-treatment method. Chromatographic separation of the ten analyte mixture was achieved within 10 minutes on a reverse-phase C8 column using mobile phase gradient mode. The multiple reaction monitoring mode (MRM) was used for analyte detection and quantification in patient samples. The chosen quantification ranges fully covered expected plasma concentrations. The method exhibited acceptable selectivity; the within- and between-run accuracy ranged from 87.2 to 113.6%, but within- and between-run precision was between 1.6 and 14.9% (at the LLOQ level CV < 20%). Although the response of the isonicotinic acid varied depending on the matrix source (CV 21.8%), validation results proved that such inconsistency does not affect the accuracy and precision of results. If stored at room temperature plasma samples should be processed within 4 h after collection, temporary storage at -20 °C up to 24 h is acceptable due to stability issues of analytes. The developed method was applied for the patient sample analysis (n = 34) receiving anti-tuberculosis treatment with the first-line drugs.

Implementation of a next-generation sequencing-based targeted approach for full-length CYP3A4 gene sequencing.

Aim: To evaluate the application of next-generation sequencing-based targeted protocol for full-length CYP3A4 gene sequencing analysis. Materials & methods: The developed sequencing protocol was applied to analyze human DNA samples (n = 7) obtained from tuberculosis patients admitted to the Riga East University Hospital, Center of Tuberculosis and Lung diseases. Results: The sequencing data quality was sufficient for the detection of already known genetic variants, as well as for identifying rare and novel variants dispersed throughout the CYP3A4 gene with a high degree of confidence. Conclusion: Developed protocol can be applied in subpopulation level association studies to determine whether specific genetic variants or variant combinations from multiple regions of the CYP3A4 gene are of clinical significance.

The prevalence of mitochondrial mutations associated with aminoglycoside-induced deafness in ethnic Latvian population: the appraisal of the evidence.

Aminoglycosides are potent antibiotics which are used to treat severe gram-negative infections, neonatal sepsis, and multidrug-resistant tuberculosis. Ototoxicity is a well-known side effect of aminoglycosides, and a rapid, profound, and irreversible hearing loss can occur in predisposed individuals. MT-RNR1 gene encoding the mitochondrial ribosomal 12S subunit is a hot spot for aminoglycoside-induced hearing loss mutations, however, a variability in the nature and frequency of genetic changes in different populations exists. The objective of this study was to analyze MT-RNR1 gene mutations in a Baltic-speaking Latvian population, and to estimate the prevalence of such genetic changes in the population-specific mitochondrial haplogroups. In the cohort of 191 ethnic non-related Latvians, the presence of two deafness-associated mutations, m.1555A>G and m.827A>G, three potentially pathogenic variations, m.961insC(n), m.961T>G and m.951G>A, and one unknown substitution, m961T>A was detected, and the aggregate frequency of all variants was 7.3%. All genetic changes were detected in samples belonged to the haplogroups H, U, T, and J. The presence of several aminoglycoside ototoxicity-related MT-RNR1 gene mutations in Baltic-speaking Latvian population indicates the necessity to include ototoxicity-related mutation analysis in the future studies in order to determine the feasibility of DNA screening for patients before administration of aminoglycoside therapy.