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.

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.

Comparative Microbiome Analysis of Three Epidemiologically Important Tick Species in Latvia.

(1) Background: Amplicon-based 16S rRNA profiling is widely used to study whole communities of prokaryotes in many niches. Here, we comparatively examined the microbial composition of three tick species, Ixodes ricinus, Ixodes persulcatus and Dermacentor reticulatus, which were field-collected in Latvia. (2) Methods: Tick DNA samples were used for microbiome analysis targeting bacterial 16S rDNA using next-generation sequencing (NGS). (3) Results: The results showed significant differences in microbial species diversity and composition by tick species and life stage. A close similarity between microbiomes of I. ricinus and I. persulcatus ticks was observed, while the D. reticulatus microbiome composition appeared to be more distinct. Significant differences in alpha and beta microbial diversity were observed between Ixodes tick life stages and sexes, with lower taxa richness indexes obtained for female ticks. The Francisella genus was closely associated with D. reticulatus ticks, while endosymbionts Candidatus Midichlorii and Candidatus Lariskella were associated with I. ricinus and I. persulcatus females, respectively. In I. ricinus females, the endosymbiont load negatively correlated with the presence of the Rickettsia genus. (4) Conclusions: The results of this study revealed important associations between ticks and their microbial community and highlighted the microbiome features of three tick species in Latvia.

Assessment of Amikacin- and Capreomycin-Related Adverse Drug Reactions in Patients with Multidrug-Resistant Tuberculosis and Exploring the Role of Genetic Factors.

Following the introduction of all-oral treatment regimens for patients with drug-resistant tuberculosis (TB), second-line injectable drug applications have been reduced in the last few years. However, they are still important for anti-TB therapy. This study aims to analyze the occurrence of amikacin- and capreomycin-related adverse drug reactions (ADR) in patients with multidrug-resistant tuberculosis (MDR-TB) and evaluate the role of multiple patient-, disease-, and therapy-related factors on the frequency of the observed adverse events. In addition, the possible role of genetic risk factors was studied by full-length mitochondrial DNA sequencing. Toward this aim, we retrospectively evaluated 47 patients with MDR-TB who received amikacin and/or capreomycin. In total, 16 (34.0%) patients developed ototoxicity and 13 (27.7%) developed nephrotoxicity, including 3 (6.4%) patients who experienced both adverse events. Ototoxicity development was more common in patients who received amikacin. No other factors showed a significant impact. Nephrotoxicity was likely associated with previous renal health impairment. Full mitochondrial genome sequencing did not reveal any specific ADR-associated variants, and results showed no differences in adverse event occurrence for any specific variants, mutation count, or mitochondrial haplogroup. The absence of the previously reported ototoxicity-related mtDNA variants in our patients with ototoxicity and nephrotoxicity highlighted the complex nature of the ADR occurrence.

Oral microbiome variations related to ageing: possible implications beyond oral health.

The global population is getting older due to a combination of longer life expectancy and declining birth rates. Growing evidence suggests that the oral microbiota composition and distribution may have a profound effect on how well we age. The purpose of this study was to investigate age-related oral microbiome variations of supragingival plaque and buccal mucosa samples in the general population in Latvia. Our results indicated significant difference between supragingival plaque bacterial profiles of three age groups (20-40; 40-60; 60 + years). Within supragingival plaque samples, age group 20-40 showed the highest bacterial diversity with a decline during the 40-60 age period and uprise again after the age of 60. Among other differences, the important oral commensal Neisseria had declined after the age of 40. Additionally, prevalence of two well-documented opportunistic pathogens Streptococcus anginosus and Gemella sanguinis gradually rose with age within our samples. Furthermore, supragingival plaque and buccal mucosa samples significantly differed in overall bacterial composition.

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.

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.