Insight into pathogenomics and phylogeography of hypervirulent and highly-lethal Mycobacterium tuberculosis strain cluster

BACKGROUND

. The Mycobacterium tuberculosis Beijing genotype is globally spread lineage with important medical properties that however vary among its subtypes. M. tuberculosis Beijing 14717-15-cluster was recently discovered as both multidrug-resistant, hypervirulent, and highly-lethal strain circulating in the Far Eastern region of Russia. Here, we aimed to analyze its pathogenomic features and phylogeographic pattern.

RESULTS

. The study collection included M. tuberculosis DNA collected between 1996 and 2020 in different world regions. The bacterial DNA was subjected to genotyping and whole genome sequencing followed by bioinformatics and phylogenetic analysis. The PCR-based assay to detect specific SNPs of the Beijing 14717-15-cluster was developed and used for its screening in the global collections. Phylogenomic and phylogeographic analysis confirmed endemic prevalence of the Beijing 14717-15-cluster in the Asian part of Russia, and distant common ancestor with isolates from Korea (> 115 SNPs). The Beijing 14717-15-cluster isolates had two common resistance mutations RpsL Lys88Arg and KatG Ser315Thr and belonged to spoligotype SIT269. The Russian isolates of this cluster were from the Asian Russia while 4 isolates were from the Netherlands and Spain. The cluster-specific SNPs that significantly affect the protein function were identified in silico in genes within different categories (lipid metabolism, regulatory proteins, intermediary metabolism and respiration, PE/PPE, cell wall and cell processes).

CONCLUSIONS

. We developed a simple method based on real-time PCR to detect clinically significant MDR and hypervirulent Beijing 14717-15-cluster. Most of the identified cluster-specific mutations were previously unreported and could potentially be associated with increased pathogenic properties of this hypervirulent M. tuberculosis strain. Further experimental study to assess the pathobiological role of these mutations is warranted.

Insight into Population Structure and Drug Resistance of Pediatric Tuberculosis Strains from China and Russia Gained through Whole-Genome Sequencing

This study aimed to determine phenotypic and genotypic drug resistance patterns of Mycobacterium tuberculosis strains from children with tuberculosis (TB) in China and Russia, two high-burden countries for multi/extensively-drug resistant (MDR/XDR) TB. Whole-genome sequencing data of M. tuberculosis isolates from China (n = 137) and Russia (n = 60) were analyzed for phylogenetic markers and drug-resistance mutations, followed by comparison with phenotypic susceptibility data. The Beijing genotype was detected in 126 Chinese and 50 Russian isolates. The Euro-American lineage was detected in 10 Russian and 11 Chinese isolates. In the Russian collection, the Beijing genotype and Beijing B0/W148-cluster were dominated by MDR strains (68% and 94%, respectively). Ninety percent of B0/W148 strains were phenotypically pre-XDR. In the Chinese collection, neither of the Beijing sublineages was associated with MDR/pre-XDR status. MDR was mostly caused by low fitness cost mutations (rpoB S450L, katG S315T, rpsL K43R). Chinese rifampicin-resistant strains demonstrated a higher diversity of resistance mutations than Russian isolates (p = 0.003). The rifampicin and isoniazid resistance compensatory mutations were detected in some MDR strains, but they were not widespread. The molecular mechanisms of M. tuberculosis adaptation to anti-TB treatment are not unique to the pediatric strains, but they reflect the general situation with TB in Russia and China.

Molecular detection and phylogeography of the hypervirulent subtype of Mycobacterium tuberculosis Beijing 14717-15

Aim. To develop a method for the detection of Mycobacterium tuberculosis strains of the highly lethal, hypervirulent, and multidrug-resistant Beijing 14717-15-cluster and its application for screening retrospective collections of M. tuberculosis from regions of Russia. The collection included M. tuberculosis DNA samples collected within population studies. Spoligotyping and genotyping of 24 variable number of tandem repeats loci were performed according to standard protocols. Phylogenetic analysis was based on the whole genome sequencing data. Polymorphism at the genomic position 2423040 AG, specific for the Beijing 14717-15 cluster, was detected using PCR-RFLP and HhaI restriction enzyme. The bioinformatic and phylogenetic analysis of whole genome sequencing data of 205 strains of the early ancient sublineage of the M. tuberculosis Beijing genotype showed that strains with the VNTR 14717-15 profile (according to MIRU-VNTRplus.org) and those close to it were grouped into one monophyletic cluster of related strains, which we defined as Beijing 14717-15-cluster. Among the SNPs specific for the 14717-15 cluster, we chose a functionally neutral SNP at the genomic position 2423040 AG (Rv2161c Val33Ala) and developed a method for its detection in the PCR-RFLP format. Next, we applied this method to screening DNA collections from the regions of the European and Asian parts of Russia and Asian countries. The results demonstrate the presence of strains of the Beijing 14717-15 cluster only in the Asian part of Russia, with a peak in Buryatia (18%). We have developed a method for the detection of the hypervirulent and highly lethal genetic cluster M. tuberculosis Beijing 14717-15 based on the detection of a specific mutation in the Rv2161c gene using PCR followed by HhaI restriction of the PCR product and agarose gel electrophoresis to discriminate between wild-type and mutant alleles. The advantages of the proposed method are the speed, simple and unambiguous interpretation of the results, the ability to detect contamination and mixed infection, suitability for efficient analysis of large collections of M. tuberculosis strains.

Influence of hypomagnetic field on the heartbeat in zebrafish embryos

The magnetic environment may influence the functioning of the cardiovascular system. It was reported that low-frequency and static magnetic fields affect hemodynamics, heart rate, and heart rate variability in animals and humans. Moreover, recent data suggest that magnetic fields affect the circadian rhythms of physiological processes. The influence of the magnetic environment on heart functionating during early development has been studied insufficiently. We utilized transparent zebrafish embryos to evaluate the effect of the hypomagnetic field on the characteristics of cardiac function using a noninvasive optical approach based on photoplethysmographic microscopic imaging. The embryos were exposed to the geomagnetic and hypomagnetic fields from the second to the 116th hour post fertilization under a 16 h light/8 h dark cycle or constant illumination. The exposure of embryos to the hypomagnetic field in both lighting modes led to increased embryo mortality, the appearance of abnormal phenotypes, and a significant increase in the embryo’s heartbeat rate. The difference between maximal and minimal heartbeat intervals, maximal to minimal heartbeat intervals ratio, and the coefficient of variation of heartbeat rate were increased in the embryos exposed to the hypomagnetic field under constant illumination from 96 to 116 h post fertilization. The dynamics of heartbeat rate changes followed a circadian pattern in all studied groups except zebrafish exposed to the hypomagnetic field under constant illumination. The results demonstrate the importance of natural magnetic background for the early development of zebrafish. The possible mechanisms of observed effects are discussed.

Pharmacokinetic-Pharmacodynamic Determinants of Clinical Outcomes for Rifampin-Resistant Tuberculosis: A Multisite Prospective Cohort Study

BACKGROUND

Rifampin-resistant and/or multidrug-resistant tuberculosis (RR/MDR-TB) treatment requires multiple drugs, and outcomes remain suboptimal. Some drugs are associated with improved outcome. It is unknown whether particular pharmacokinetic-pharmacodynamic relationships predict outcome.

METHODS

Adults with pulmonary RR/MDR-TB in Tanzania, Bangladesh, and the Russian Federation receiving local regimens were enrolled from June 2016 to July 2018. Serum was collected after 2, 4, and 8 weeks for each drug's area under the concentration-time curve over 24 hours (AUC0-24). Quantitative susceptibility of the M. tuberculosis isolate was measured by minimum inhibitory concentrations (MICs). Individual drug AUC0-24/MIC targets were assessed by adjusted odds ratios (ORs) for favorable treatment outcome, and hazard ratios (HRs) for time to sputum culture conversion. K-means clustering algorithm separated the cohort of the most common multidrug regimen into 4 clusters by AUC0-24/MIC exposures.

RESULTS

Among 290 patients, 62 (21%) experienced treatment failure, including 30 deaths. Moxifloxacin AUC0-24/MIC target of 58 was associated with favorable treatment outcome (OR, 3.75; 95% confidence interval, 1.21-11.56; P = .022); levofloxacin AUC0-24/MIC of 118.3, clofazimine AUC0-24/MIC of 50.5, and pyrazinamide AUC0-24 of 379 mg × h/L were associated with faster culture conversion (HR >1.0, P < .05). Other individual drug exposures were not predictive. Clustering by AUC0-24/MIC revealed that those with the lowest multidrug exposures had the slowest culture conversion.

CONCLUSIONS

Amidst multidrug regimens for RR/MDR-TB, serum pharmacokinetics and M. tuberculosis MICs were variable, yet defined parameters to certain drugs-fluoroquinolones, pyrazinamide, clofazimine-were predictive and should be optimized to improve clinical outcome.

CLINICAL TRIALS REGISTRATION

NCT03559582.

Process evaluation for the adaptation, testing and dissemination of a mobile health platform to support people with HIV and tuberculosis in Irkutsk, Siberia

OBJECTIVES

We developed and tested a mobile health-based programme to enhance integration of HIV and tuberculosis (TB) care and to promote a patient-centred approach in a region of high coinfection burden. Phases of programme development included planning, stakeholder interviews and platform re-build, testing and iteration.

SETTING

In Irkutsk, Siberia, HIV/TB coinfection prevalence is high relative to the rest of the Russian Federation.

PARTICIPANTS

Pilot testing occurred for a cohort of 60 people with HIV and TB.

RESULTS

Key steps emerged to ensure the mobile health-based programme could be operational and adequately adapted for the context, including platform language adaptation, optimisation of server management, iteration of platform features, and organisational practice integration. Pilot testing of the platform rebuild yielded favourable patient perceptions of usability and acceptability at 6 months (n=47 surveyed), with 18 of 20 items showing scores above 4 (on a scale from 1 to 5) on average. Development of this mobile health-based programme for integrated care of infections highlighted the importance of several considerations for tailoring these interventions contextually, including language adaptation and technological capacity, but also, importantly, contextualised patient preferences related to privacy and communication with peers and/or providers, existing regional capacity for care coordination of different comorbidities, and infection severity and treatment requirements.

CONCLUSIONS

Our experience demonstrated that integration of care for TB and HIV can be well served by using multimodal mobile health-based programmes, which can enhance communication and streamline workflow between providers across multiple collaborating institutions and improve continuity between inpatient and outpatient care settings. Further study of programme impact on contextual disease-related stigma and social isolation as well as evaluation of implementation on a broader scale for HIV care is currently under way.

TRIAL REGISTRATION NUMBER

NCT03819374.

Transborder molecular analysis of drug-resistant tuberculosis in Mongolia and Eastern Siberia, Russia

Eastern Siberia (Russia) and Mongolia are borderline regions in Asia with a high incidence of tuberculosis (TB). In this study, we investigated the transborder transmission of Mycobacterium tuberculosis with a focus on endemic and epidemic clones and drug resistance. M. tuberculosis isolates (287 from Mongolia and 754 from Russia) were collected using cross-sectional population-based surveys between 2010 and 2016. The isolates were genotyped using 24 variable number of tandem repeat (VNTR) loci and by testing of the key markers to discriminate within the Beijing genotype. All isolates were divided into 427 MIRU-types that were assigned to Lineage 2 (Beijing) and Lineage 4 (Ural, Haarlem, Latin American-Mediterranean [LAM], S, unclassified). The Beijing genotype was dominant in both countries (69% in Russia, 75% in Mongolia). However, the Beijing isolates differed significantly between the countries, in terms of the identified subtypes. LAM was the most common non-Beijing genotype (11.1% in Mongolia and 14.9% in Russia) and LAM isolates mostly belonged to the LAM-RUS branch in both countries. The MDR rate was higher in Russia than in Mongolia among newly diagnosed patients: 29.4% versus 5.6% (p < 0.001). In Mongolia, the MDR rate was similar in Beijing (29.7%) and non-Beijing (27.5%) genotypes. In Russia, a higher MDR rate was observed in (i) Beijing compared to non-Beijing (48.7% versus 38.3%, p = 0.03) and (ii) Beijing B0/W148 compared to Beijing Central Asian/Russian (63.4% versus 37.3%, p<0.001). In conclusion, the M. tuberculosis population structure in Mongolia was shaped by mainly historical interaction with China (dominance of the Beijing genotype) and Northern Eurasia (presence of the LAM-RUS branch). In contrast, the transborder transmission of M. tuberculosis since the 1990s between Mongolia and its neighbors has been negligible, and the adverse trends of MDR-TB in Russia did not impact the current situation in Mongolia. This article is protected by copyright. All rights reserved.

Extremely lethal and hypervirulent Mycobacterium tuberculosis strain cluster emerging in Far East, Russia

Mycobacterium tuberculosis strains of the early ancient sublineage of the Beijing genotype are mostly drug susceptible and mainly circulate in East Asia. We have recently discovered two clusters of this sublineage emerging in the Asian part of Russia (VNTR-defined 1071-32 and 14717-15 types) and, to our surprise, both were strongly MDR/XDR-associated. Here, we evaluated their pathogenic features. The clinical isolates and reference laboratory strain H37Rv were investigated in the C57BL/6 mouse model to assess their virulence and lethality properties. The BACTEC MGIT 960 system was used to study the in vitro growth characteristics. In the murine model, strains 396 (14717-15-cluster, from Buryatia, Far East) and 6691 (1071-32-cluster, from Omsk, Siberia) demonstrated contrasting properties. The 396-infected group had significantly higher mortality, more weight loss, higher bacterial burden, and more severe lung pathology. Furthermore, compared to the previously published data on other Russian epidemic Beijing strains (B0/W148, CAO, Central Asian Russian), strain 396 demonstrated the highest mortality. Under the in vitro growth experiment, cluster 14717-15 isolates had significantly shorter lag-phase. To conclude, low-virulent MDR strain 6691 belongs to the Beijing 1071-32-cluster widespread across FSU countries but at low prevalence. This corresponds to common expectation that multiple drug resistance mutations reduce fitness and virulence. In contrast, highly lethal and hypervirulent MDR strain 396 represents an intriguing Beijing 14717-15 cluster predominant only in Buryatia, Far East (16%), sporadically found beyond it, but not forming clusters of transmission. Further in-depth study of this most virulent Russian Beijing cluster is warranted.

Implementation of a Mobile Health Strategy to Improve Linkage to and Engagement with HIV Care for People Living with HIV, Tuberculosis, and Substance Use in Irkutsk, Siberia

In Irkutsk, Siberia, there is a high prevalence of HIV and tuberculosis (TB) coinfection. Mobile health (mHealth) strategies have shown promise for increasing linkage to and engagement in care for people living with HIV (PLWH) in other contexts. We evaluated outcomes for a cohort of PLWH, TB, and substance use in Irkutsk after participation in a multi-feature mHealth intervention called MOCT. Sixty patients were enrolled during hospitalization for TB. We evaluated participant app usage, linkage to HIV care postdischarge, perception of self-efficacy related to HIV care, and HIV-related clinical outcomes at 6 months. We also performed an exploratory analysis to compare a subset of 49 patients with a pre-intervention cohort matched for age and gender. Participants demonstrated engagement with app features examined at 6 months. The majority linked to HIV care by 6 months (83%). Self-scoring of confidence in ability to communicate with HIV providers improved from baseline (median score 8, scale 1-10) to 6 months (10, p = 0.004). A higher proportion of the MOCT subset refilled antiretroviral therapy (69% vs. 43% in pre-intervention cohort, p = 0.01), with fewer deaths in the MOCT subset at 6 months (1 death vs. 10 deaths in pre-intervention cohort, p = 0.02) and a decreased likelihood of developing the composite outcome of death/failure to achieve viral suppression at 6 months (adjusted odds ratio = 0.33, p = 0.029). This study demonstrates preliminary intervention uptake and improvement in short-term outcomes for an urban cohort of PLWH, TB, and substance use enrolled in a multi-feature mHealth intervention, a novel strategy for the context. Clinical Trial Registration Number: NCT03819374.

New epidemic cluster of pre-extensively drug resistant isolates of Mycobacterium tuberculosis Ural family emerging in Eastern Europe

Background

Ural genetic family is a part of the Euro-American lineage of Mycobacterium tuberculosis and is endemic in Northern Eurasia (former Soviet Union [FSU]). These strains were long described as drug susceptible and of low virulence, but recent studies reported an increasing circulation of the multidrug-resistant (MDR) and extensively drug-resistant (XDR) Ural strains. Here, we analyzed all publicly available whole genome sequence data of Ural genotype isolates, in order to elucidate their phylogenomic diversity with a special focus on MDR and potentially epidemic clones.

Results

A total of 149 M. tuberculosis genomes of Ural isolates from FSU countries were mined from the GMTV database and TB-ARC project. We identified 6002 variable amino acid positions that were assessed for functional significance and used to build ML, NJ trees and for Bayesian TMRCA estimation. Three robust monophyletic clades were identified: Clade A (31 isolates from Russia, Belarus, Moldova), Clade B (52 isolates from Russia), and Clade C (37 isolates from Moldova, 2 from Belarus). Clade C was significantly associated with XDR or pre-XDR status compared to the pooled Clades A and B (33/39 versus 5/83, P < 0.0001). Time of origin was estimated for Clade A at 77.7–137 years ago and for Clade B at 56.3–99.2 years ago compared to the significantly more recent origin for Clade C. in silico spoligotyping identified signatures specific of the Clade A (spoligotype SIT35), and Clades B and C (both SIT262).

Conclusions

A genetically compact and evolutionarily young Ural Clade C, likely originated after collapse of the Soviet Union, and reached epidemic proportions in Moldova in the last 20 years. This epidemic pre-XDR clone (mostly rifampin, isoniazid and kanamycin resistant) is characterized by a specific combination of mutations: KatG Ser315Thr, fabG1 -15C > T, RpoB Ser450Leu, RpsL Lys88Arg, eis -12G > A and EmbB Ser297Ala/T > G. Its further dissemination may occur towards both Russia and European Union and should be taken into consideration by health authorities. The identified spoligotyping signatures can serve for rapid preliminary detection and surveillance of the more hazardous pre-XDR associated strains of the Ural family, both in populations from countries of their endemic circulation and migrant communities.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5162-3) contains supplementary material, which is available to authorized users.

Undertreated HIV and drug-resistant tuberculosis at a referral hospital in Irkutsk, Siberia

SETTING

A referral hospital for tuberculosis (TB) in Irkutsk, the Russian Federation.

OBJECTIVE

To describe disease characteristics, treatment and hospital outcomes of TB-HIV (human immunodeficiency virus).

DESIGN

Observational cohort of HIV-infected patients admitted for anti-tuberculosis treatment over 6 months.

RESULTS

A total of 98 patients were enrolled with a median CD4 count of 147 cells/mm(3) and viral load of 205 943 copies/ml. Among patients with drug susceptibility testing (DST) results, 29 (64%) were multidrug-resistant (MDR), including 12 without previous anti-tuberculosis treatment. Nineteen patients were on antiretroviral therapy (ART) at admission, and 10 (13% ART-naïve) were started during hospitalization. Barriers to timely ART initiation included death, in-patient treatment interruption, and patient refusal. Of 96 evaluable patients, 21 (22%) died, 14 (15%) interrupted treatment, and 10 (10%) showed no microbiological or radiographic improvement. Patients with a cavitary chest X-ray (aOR 7.4, 95%CI 2.3-23.7, P = 0.001) or central nervous system disease (aOR 6.5, 95%CI 1.2-36.1, P = 0.03) were more likely to have one of these poor outcomes.

CONCLUSION

High rates of MDR-TB, treatment interruption and death were found in an HIV-infected population hospitalized in Irkutsk. There are opportunities for integration of HIV and TB services to overcome barriers to timely ART initiation, increase the use of anti-tuberculosis regimens informed by second-line DST, and strengthen out-patient diagnosis and treatment networks.

Administrative interventions associated with increased initiation on antiretroviral therapy in Irkutsk, Siberia

A bundle of initiatives to integrate human immunodeficiency virus (HIV) and tuberculosis (TB) services was assessed for the impact on antiretroviral therapy (ART) initiation at a TB referral hospital in Irkutsk, Russian Federation, from February 2014 to December 2015. The ART initiation rates in 166 ART-naïve patients undergoing anti-tuberculosis treatment (34.1% with multidrug or extensively drug-resistant TB) increased significantly from 14 (17%) pre-intervention to 44 (54%) post-intervention (P < 0.001). A survey of TB hospital staff identified administrative prioritisation as the most important initiative for increasing ART initiation.

Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography

Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans.

Primary multidrug-resistant Mycobacterium tuberculosis in 2 regions, Eastern Siberia, Russian Federation

Of 235 Mycobacterium tuberculosis isolates from patients who had not received tuberculosis treatment in the Irkutsk oblast and the Sakha Republic (Yakutia), eastern Siberia, 61 (26%) were multidrug resistant. A novel strain, S 256, clustered among these isolates and carried eis-related kanamycin resistance, indicating a need for locally informed diagnosis and treatment strategies.

'Lethal' combination of Mycobacterium tuberculosis Beijing genotype and human CD209 -336G allele in Russian male population

An interaction of different human alleles and endemic bacterial strains may be clinically manifested as different outcome of the disease in different hosts infected with the same genotype. The primary objective of this study was to investigate this issue in the model of Mycobacterium tuberculosis and human DC-SIGN encoding CD209 promoter SNP (rs4804803) in Russian Siberian population. We sought to find a possible combination of M. tuberculosis lineage and human host allele/genotype correlating with unfavorable outcome of the disease. The 101 paired DNA samples from patients with pulmonary TB (human and M. tuberculosis DNA) were studied by 12-loci MIRU-VNTR typing (M. tuberculosis strains) and CD209 -336 A/G typing (human DNA). Ninety autopsy DNA samples as a source of human and mycobacterial nucleic acids from persons who died from TB were also subjected to the same genotyping procedures. A human control group consisted of 177 healthy individuals. The Beijing genotype was more frequently identified in autopsy versus patient samples, in 70.0% and 51.5%, respectively (χ(2)=6.06, P=0.01). Regarding other M. tuberculosis genetic families, no significant difference in LAM family distribution among patient strains and autopsy samples has been found. In contrast, Ural genotype was significantly less frequently detected in the autopsy samples (χ(2)=6.12, P=0.01). Allelic and genotypic frequencies of the CD209 -336A/G did not differ significantly under global comparison when contrasting controls versus patients versus autopsy samples. However intriguing and contrasting significant associations were found in the male subgroup under M. tuberculosis genotype-stratified comparisons. Firstly, male carriers of -336AA genotype were more frequently infected with Beijing genotype (χ(2)=5.2, P=0.02). Secondly and remarkably, this association was inverted in the autopsy sample: male carriers of -336AA genotype died less frequently due to TB caused by a Beijing rather than a non-Beijing strain (χ(2)=5.37, P=0.02). In conclusion, we hypothesize that although carriers of CD209 -336A allele are more sensitive to infection with a Beijing strain, a combination of human CD209 -336G allele and M. tuberculosis Beijing genotype leads more frequently to the lethal outcome in pulmonary TB male patients in Russian (Caucasian) population.