Association between genotype and drug resistance profiles of Mycobacterium tuberculosis strains circulating in China in a national drug resistance survey

Drug-Resistant Characteristics, Genetic Diversity, and Transmission Dynamics of Multidrug-Resistant Mycobacterium tuberculosis in Jiangxi, China

Purpose

In this study, we aimed to determine the transmission pattern of multidrug-resistant tuberculosis (MDR-TB) isolates circulating in Jiangxi Province with whole-genome sequencing (WGS). In addition, we also sought to describe mutational resistome of MDR-TB isolates.

Patients and Methods

A total of 115 MDR-TB isolates determined by the phenotypic proportion method of drug susceptibility testing between January 2018 and December 2022 from provincial drug surveillance (DRS) in Jiangxi were included in our analysis. The demographic data and treatment history were extracted from the National TB Registry System. WGS was used to analyze the genotypic characteristics of drug resistance and transmissions.

Results

About 62.6% of MDR-TB strains were isolated from cases that received previous anti-tuberculosis treatment. According to the WGS results, 96.5% were genotypic MDR-TB, and more than half of MDR-TB isolates tested were also resistant to streptomycin (59.1%), ethambutol (56.5%), and fluroquinolones (53.0%), while resistance to cycloserine and linezolid was lowest, only in two (1.7%) and one (0.9%) isolate, respectively. Ser450Leu in rpoB (57.9%), Ser315Thr in katG (74.1%), Met306Val in embB (40.0%), Lys43Arg in rpsL (75.0%), Ala90Val in gyrA (32.8%) were predominant mutant types among the rifampin-, isoniazid-, ethambutol-, streptomycin-, fluoroquinolones-resistant isolates, respectively. Lineage 2 (East Asian genotype) occurred at the highest frequency with 97 cases (84.3%), followed by lineage 4 (Euro-American genotype) with 18 cases (15.7%). Additionally, 5 clusters consisting of 10 isolates were identified in the present study, demonstrating a clustering rate of 8.7%.

Conclusion

MDR/Rifampicin-Resistant (RR)-TB epidemic in this region is driven by lineage 2 clade that also show higher resistance to other anti-tuberculosis drugs. Lower cluster rates compared with a relatively higher proportion of new MDR-TB cases indicate that a considerable number of MDR-TB cases remain undiagnosed.

Advancing tuberculosis management: the role of predictive, preventive, and personalized medicine

Tuberculosis is a major global health issue, with approximately 10 million people falling ill and 1.4 million dying yearly. One of the most significant challenges to public health is the emergence of drug-resistant tuberculosis. For the last half-century, treating tuberculosis has adhered to a uniform management strategy in most patients. However, treatment ineffectiveness in some individuals with pulmonary tuberculosis presents a major challenge to the global tuberculosis control initiative. Unfavorable outcomes of tuberculosis treatment (including mortality, treatment failure, loss of follow-up, and unevaluated cases) may result in increased transmission of tuberculosis and the emergence of drug-resistant strains. Treatment failure may occur due to drug-resistant strains, non-adherence to medication, inadequate absorption of drugs, or low-quality healthcare. Identifying the underlying cause and adjusting the treatment accordingly to address treatment failure is important. This is where approaches such as artificial intelligence, genetic screening, and whole genome sequencing can play a critical role. In this review, we suggest a set of particular clinical applications of these approaches, which might have the potential to influence decisions regarding the clinical management of tuberculosis patients.

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.

Genotyping of Mycobacterium tuberculosis complex isolated from humans and animals in northeastern Iran

The objective of this study was to genotype Mycobacterium tuberculosis complex isolated from humans and cattle in northern Iran. Over the course of one year, a total of 120 human and 21 cattle isolates were tested using region of difference (RD)-based polymerase chain reaction (PCR) and mycobacterial interspersed repetitive unites-variable number tandem repeats (MIRU-VNTR). In M. tuberculosis, out of 120 isolates investigated, the most common genotype detected was NEW-1 (53.3%), followed by CAS/ Delhi (24.1%), Haarlem (5%), Beijing (4.16%), Uganda I (4.16%), S (3.3%), Ural (0.83%), TUR (0.83%), Uganda II (0.83%), Lam (0.83%) and Cameroon (0.83%). The HGDI rate was 0.9981 and the clustering rate was 10.83. Of the isolates, QUB26 had the highest allele diversity (h: 0.76), while the loci Mtub29 and MIRU24 had the lowest (h: 0). In M. Bovis, out of 123 collected tissue samples, 21 (17%) grew on culture media. The HGDI rate was 0.71 and clustering rate was 85.7%. The locus ETRC had the highest allele diversity (h: 0.45). The findings of this study suggest that there is high genetic diversity among M. tuberculosis isolates in Khorasan Razavi Province, which is consistent with similar results from other studies in other provinces in Iran and neighboring countries. This indicates that the prevalent genotypes in this study are spreading in the Middle East region. Furthermore, considering that M. Bovis isolates were identified in two clusters, it seems that all of them have a common origin and are circulating among the livestock farms in the province.

Drug Resistance of Mycobacterium tuberculosis Based on Whole-Genome Sequencing in the Yi Ethnic Group, Sichuan Province, China

This study investigated drug-resistant tuberculosis (DR-TB) in the Yi ethnic group. The study was designed to identify risk factors for DR-TB and its relationship with HIV/AIDS. To establish the resistance to antituberculosis drugs, whole-genome sequencing (WGS) was performed using culture-positive Mycobacterium tuberculosis samples collected from people of the Yi ethnic group from March 2019 to March 2021. Baseline characteristics were obtained from China's tuberculosis surveillance system. A total of 116 M. tuberculosis strains were included in the final analysis. Lineage 2.2 (75.86%) was the dominant sublineage, followed by lineage 4.5 (18.97%) and lineage 4.4 (5.17%). The rates of rifampicin-resistant (RR-TB), multidrug-resistant (MDR-TB), and preextensively drug-resistant TB (pre-XDR-TB) were 18.97%, 10.34%, and 6.03%, respectively. Drug-resistant strains were not found in the elderly (age ≥ 65 years). The proportions of RR/MDR-TB and pre-XDR-TB cases among re-treatment patients were higher than those among new patients (χ ² = 12.155, P = 0.003; χ ² = 22.495, P = 0.001, respectively). The pre-XDR-TB case proportions were higher among female patients than among males and higher among referred patients (χ ² = 5.456, P = 0.032; χ ² = 15.134, P = 0.002, respectively). The rates of RR/MDR-TB and pre-XDR-TB did not differ appreciably among groups with different HIV infection statuses nor lineage populations. DR-TB poses a serious challenge to the Yi ethnic group. Re-treatment patients, women, and referred patients were at high risk of MDR/RR-TB or pre-XDR-TB while HIV and lineage 2 had negligible association with drug resistance. Whole-genome sequencing should be used to guide the design of treatment regimens and to tailor public interventions.

Recent Transmission and Prevalent Characterization of the Beijing Family Mycobacterium tuberculosis in Jiangxi, China

The Beijing genotype is the most common type of tuberculosis in Jiangxi Province, China. The association of population characteristics and their prevalence in the development of recent transmission is still unclear. 1,433 isolates were subjected to drug-resistant tests and MIRU-VNTR analysis. We compared differences in demographic characteristics and drug resistance patterns between the Beijing and non-Beijing family strains. We also explored the association of the clustering rate with the Beijing genotype of Mycobacterium tuberculosis. The Beijing genotype was dominant (78.16%). The results of MIRU-VNTR showed that 775 of 1,433 strains have unique patterns, and the remaining gather into 103 clusters. A recent transmission rate was 31.54% (452/1,433). The Beijing genotype strains were more likely to spread among the recurrent population (p = 0.004), people less than 50 years of age (p = 0.02 or 0.003), and the personnel in the northern regions (p = 0.03). Drug resistance patterns did not show significant differences between Beijing and non-Beijing genotype isolates. Furthermore, we found that HIV-positive cases had a lower clustering rate (p = 0.001). Our results indicated that the recurrent population and people under 50 years of age were more likely to be infected with the Beijing genotype of M. tuberculosis. The strains from the Beijing family were easier to cluster compared to strains isolated from the non-Beijing family. Social activity and AIDS substantially impacted the clustering rate of the Beijing genotype of M. tuberculosis. Multidrug resistant M. tuberculosis affected Beijing genotype transmission.

Distribution and identification of Mycobacterium tuberculosis lineage in Kashgar prefecture

Objectives

Kashgar prefecture is an important transportation and trade hub with a high incidence of tuberculosis. The following study analyzed the composition and differences in Mycobacterium tuberculosis (M.tb) lineage and specific tags to distinguish the lineage of the M.tb in Kashgar prefecture, thus providing a basis for the classification and diagnosis of tuberculosis in this area.

Methods

Whole-genome sequencing (WGS) of 161 M.tb clinical strains was performed. The phylogenetic tree was constructed using Maximum Likelihood (ML) based on single nucleotide polymorphisms (SNPs) and verified through principal component analysis (PCA). The composition structure of M.tb in different regions was analyzed by combining geographic information.

Results

M.tb clinical strains were composed of lineage 2 (73/161, 45.34%), lineage 3 (52/161, 32.30%) and lineage 4 (36/161, 22.36%). Moreover, the 3 lineages were subdivided into 11 sublineages, among which lineage 2 included lineage 2.2.2/Asia Ancestral 1 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 2 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 3 (18/73, 24.66%), and lineage 2.2.1-Modern Beijing (39/73, 53.42%). Lineage 3 included lineage 3.2 (14/52, 26.92%) and lineage 3.3 (38/52, 73.08%), while lineage 4 included lineage 4.1 (3/36, 8.33%), lineage 4.2 (2/36, 5.66%), lineage 4.4.2 (1/36, 2.78%), lineage 4.5 (28/36, 77.78%) and lineage 4.8 (2/36, 5.66%), all of which were consistent with the PCA results. One hundred thirty-six markers were proposed for discriminating known circulating strains. Reconstruction of a phylogenetic tree using the 136 SNPs resulted in a tree with the same number of delineated clades. Based on geographical location analysis, the composition of Lineage 2 in Kashgar prefecture (45.34%) was lower compared to other regions in China (54.35%-90.27%), while the composition of Lineage 3 (32.30%) was much higher than in other regions of China (0.92%-2.01%), but lower compared to the bordering Pakistan (70.40%).

Conclusion

Three lineages were identified in M.tb clinical strains from Kashgar prefecture, with 136 branch-specific SNP. Kashgar borders with countries that have a high incidence of tuberculosis, such as Pakistan and India, which results in a large difference between the M.tb lineage and sublineage distribution in this region and other provinces of China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07307-4.

The Neglected Contribution of Streptomycin to the Tuberculosis Drug Resistance Problem

The airborne pathogen Mycobacterium tuberculosis is responsible for a present major public health problem worsened by the emergence of drug resistance. M. tuberculosis has acquired and developed streptomycin (STR) resistance mechanisms that have been maintained and transmitted in the population over the last decades. Indeed, STR resistant mutations are frequently identified across the main M. tuberculosis lineages that cause tuberculosis outbreaks worldwide. The spread of STR resistance is likely related to the low impact of the most frequent underlying mutations on the fitness of the bacteria. The withdrawal of STR from the first-line treatment of tuberculosis potentially lowered the importance of studying STR resistance. However, the prevalence of STR resistance remains very high, could be underestimated by current genotypic methods, and was found in outbreaks of multi-drug (MDR) and extensively drug (XDR) strains in different geographic regions. Therefore, the contribution of STR resistance to the problem of tuberculosis drug resistance should not be neglected. Here, we review the impact of STR resistance and detail well-known and novel candidate STR resistance mechanisms, genes, and mutations. In addition, we aim to provide insights into the possible role of STR resistance in the development of multi-drug resistant tuberculosis.

A Prospective Cohort Study on the Prevalent and Recurrent Tuberculosis Isolates Using the MIRU-VNTR Typing

The study aims to describe the clustering characteristics of Mycobacterium tuberculosis (M.tb) strains circulating in eastern China and determine the ratio of relapse and reinfection in recurrent patients. We recruited sputum smear-positive pulmonary tuberculosis cases from five cities of Jiangsu Province, China, during August 2013 and December 2015. Patients were followed for the treatment outcomes and recurrence based on a cohort design. M.tb strains were isolated and genotyped using the 12-locus MIRU-VNTR. The Beijing family was identified by the extended Region of Difference (RD) analysis. The Hunter-Gaston Discriminatory Index (HGDI) was used to judge the resolution ability of MIRU-VNTR. The odds ratio (OR) together with 95% confidence interval (CI) were used to estimate the strength of association. We performed a cluster analysis on 2098 M.tb isolates and classified them into 545 genotypes and five categories (I, 0.19%; II, 0.43%; III, 3.34%; IV, 77.46%; V, 18.59%). After adjusting for potential confounders, the Beijing family genotype (OR = 118.63, 95% CI: 79.61–176.79, P = 0.001) was significantly related to the dominant strain infections. Patients infected with non-dominant strains had a higher risk of the pulmonary cavity (OR = 1.39, 95% CI: 1.01–1.91, P = 0.046). Among 37 paired recurrent cases, 22 (59.46%) were determined as endogenous reactivation, and 15 (40.54%) were exogenous reinfection. The type of M.tb strains prevalent in Jiangsu Province is relatively single. Beijing family strains infection is dominant in local tuberculosis cases. Endogenous reactivation appears to be a major cause of recurrent tuberculosis in Eastern China. This finding emphasizes the importance of case follow-up and monitoring after the completion of antituberculosis treatment.

Molecular Characteristics and Drug Resistance of Mycobacterium tuberculosis Isolate Circulating in Shaanxi Province, Northwestern China

Objective: Shaanxi is the most highly populated province with high burdens of tuberculosis in northwestern China. The aim of this study was to investigate the molecular characteristics and drug resistance of Mycobacterium tuberculosis isolates from Shaanxi province of China in 2018. Methods: Phenotypic drug susceptibility testing and spoligotyping methods were performed on 518 M. tuberculosis isolates; drug-resistant isolates were sequenced in 11 drug loci, including katG, inhA, oxyR-ahpC, rpoB, embB, rpsL, rrs₁ (nucleotides 388-1084), gyrA, gyrB, rrs₂ (nucleotides 1158-1674), and eis. Results: The prevalences of isoniazid, rifampicin, ethambutol, streptomycin, ofloxacin, and kanamycin resistance were 22.0%, 19.3%, 7.9%, 23.8%, 10.4%, and 3.3%, respectively. The Beijing family (82.8%) was the predominant genotype, followed by the T (9.3%), H (0.6%), CAS (0.4%), LAM (0.4%), and U (0.4%) families. The percentage of Beijing genotype in a central area (88.1%) was higher than in the south (77.3%) and the north area (80.1%) (p < 0.05), while the sex, age, and treatment history between Beijing and non-Beijing family were not statistically different. Mutation analysis found that the most prevalent mutations were katG315, rpoB531, embB306, rpsL43, gyrA94, and rrs1401; the Beijing family exhibited a high rate of isoniazid-resistant isolates carrying katG315 mutations (p < 0.05). Furthermore, compared with the phenotypic data, the sensitivities of isoniazid, rifampicin, ethambutol, streptomycin, ofloxacin, and kanamycin resistance by sequencing base on 11 loci were 85.1%, 94.0%, 53.7%, 74.8%, 77.8%, and 64.7%, respectively. Conclusions: Shaanxi has a serious epidemic of drug-resistant tuberculosis, Beijing family is the predominant genotype, and the distribution showed geographic diversity. The prevalence of Beijing genotypes has a tendency to promote the transmission of high-level isoniazid-resistant M. tuberculosis. Besides, the hot spot regions localized in the embB, rrs₂, and eis gene appear not to serve as excellent biomarkers for predicting ethambutol and kanamycin resistance in Shaanxi.

Evaluation of the frequency of mutation genes in multidrug-resistant tuberculosis (MDR-TB) strains in Beijing, China

The aim of this study was to explore the frequency and distribution of gene mutations that are related to isoniazid (INH) and rifampin (RIF)-resistance in the strains of the multidrug-resistant tuberculosis (MDR-TB) Mycobacterium tuberculosis (M.tb) in Beijing, China. In this retrospective study, the genotypes of 173 MDR-TB strains were analysed by spoligotyping. The katG, inhA genes and the promoter region of inhA, in which genetic mutations confer INH resistance; and the rpoB gene, in which genetic mutations confer RIF resistance, were sequenced. The percentage of resistance-associated nucleotide alterations among the strains of different genotypes was also analysed. In total, 90.8% (157/173) of the MDR strains belonged to the Beijing genotype. Population characteristics were not significantly different among the strains of different genotypes. In total, 50.3% (87/173) strains had mutations at codon S315T of katG; 16.8% (29/173) of strains had mutations in the inhA promoter region; of them, 5.5% (15/173) had point mutations at -15 base (C→T) of the inhA promoter region. In total, 86.7% (150/173) strains had mutations at rpoB gene; of them, 40% (69/173) strains had mutations at codon S531L of rpoB. The frequency of mutations was not significantly higher in Beijing genotypic MDR strains than in non-Beijing genotypes. Beijing genotypic MDR-TB strains were spreading in Beijing and present a major challenge to TB control in this region. A high prevalence of katG Ser315Thr, inhA promoter region (-15C→T) and rpoB (S531L) mutations was observed. Molecular diagnostics based on gene mutations was a useful method for rapid detection of MDR-TB in Beijing, China.

Primary Drug-Resistance Pattern and Trend in Elderly Tuberculosis Patients in Shandong, China, from 2004 to 2019

BACKGROUND

With an aging population, China is facing a huge burden of elderly patients with drug resistant tuberculosis (DR-TB), which has become a significant obstacle for the global TB control. There is still little study on DR-TB in the elderly in China so far. Thus, more research on the epidemiological characteristics and trend of primary DR-TB among the elderly will be necessary.

METHODS

A retrospective study was conducted in Shandong, China from 2004 to 2019. We collected 12,661 primary TB cases, of which 4368 elderly (≥60 years) primary TB cases were involved. Clinical characteristics including age, sex, cavity, smoking, drinking, comorbidity and drug susceptibility data were collected from 36 TB prevention and control institutions of Shandong Province. Sputum samples were collected by each surveillance site, and examined in the TB Reference Laboratory of SPCH. Descriptive statistical analysis, chi-square and linear regression were used for analyzing.

RESULTS

Among 4368 elderly patients with primary TB, the DR-TB and multi-resistant tuberculosis (MDR-TB) accounted for 17.19% and 2.29%, respectively. During 2004-2019, the proportions of MDR-TB, polydrug resistant tuberculosis (PDR-TB), rifampin (RFP)-resistance increased by 160.00%, 18.18%, 231.82%, respectively and the rate of DR-TB among elderly patients with primary cavitary TB increased by 255%. Among the elderly with primary DR-TB during 2004-2019, the proportion of male (from 85.19 to 89.06), cavity (from 7.41 to 46.88), RFP-resistance (from 3.70 to 21.88), and streptomycin (SM)-resistance (from 37.04 to 62.5) increased significantly (P<0.05). And the proportion of female (from 14.81 to 10.94), non-cavity (from 92.59 to 32.81), INH-resistance (from 66.67 to 57.81) decreased significantly (P<0.05).

CONCLUSION

Among the elderly, the proportions of MDR-TB, PDR-TB, RFP-resistance and cavitary DR-TB increased significantly. The pattern of DR-TB changed from female, non-cavity and INH-resistant groups to male, cavity, RFP or SM-resistant groups. For a better control on the elderly DR-TB in the future, we should pay more attention to male, smoking, drinking, chronic obstructive pulmonary disease (COPD) and diabetes subgroups and take targeted measures to control these subgroups.

Molecular epidemiology of clinical Mycobacterium tuberculosis complex isolates in South Omo, Southern Ethiopia

BACKGROUND

Tuberculosis (TB) is caused by Mycobacterium tuberculosis complex (MTBC). Mapping the genetic diversity of MTBC in high TB burden country like Ethiopia is important to understand principles of the disease transmission and to strengthen the regional TB control program. The aim of this study was to investigate the genetic diversity of Mycobacterium tuberculosis complex (MTBC) isolates circulating in the South Omo, southern Ethiopia.

METHODS

MTBC isolates (N = 156) were genetically analyzed using spacer oligotyping (spoligotyping) and mycobacterial interspersed repetitive unit-variable number of tandem repeat (MIRU-VNTR) typing. Major lineages and lineages were identified using MTBC databases. Logistic regression was used to correlate patient characteristics with strain clustering.

RESULTS

The study identified Euro-American (EA), East-African-Indian (EAI), Indo-Oceanic (IO), Lineage_7/Aethiops vertus, Mycobacterium bovis and Mycobacterium africanum major lineages in proportions of 67.3% (105/156), 22.4% (35/156), 6.4% (10/156), 1.9% (3/156), 1.3% (2/156) and 0.6% (1/156), respectively. Lineages identified were Delhi/CAS 23.9% (37/155), Ethiopia_2 20.6% (32/155), Haarlem 14.2% (22/155), URAL 14.2%(22/155), Ethiopia_3 8.4% (13/155), TUR 6.5% (10/155), Lineage_7/Aethiops vertus 1.9% (3/155), Bovis 1.3% (2/155), LAM 1.3% (2/155), EAI 0.6% (1/155), X 0.6% (1/155) and Ethiopia H₃₇Rv-like strain 0.6% (1/155). Of the genotyped isolates 5.8% (9/155) remained unassigned. The recent transmission index (RTI) was 3.9%. Orphan strains compared to shared types (AOR: 0.09, 95% CI: 0.04-0.25) were associated with reduced odds of clustering. The dominant TB lineage in pastoral areas was EAI and in non-pastoral areas was EA.

CONCLUSION

The epidemiological data, highly diverse MTBC strains and a low RTI in South Omo, provide information contributing to the TB Control Program of the country.

Resistance of Mycobacterium tuberculosis strains to Rifampicin: A systematic review and meta-analysis

Introduction

Antitubercular drug resistance strain is a horrifying barrier to effective TB treatment and prevention. The present study aimed to determine the prevalence and geographical distribution of rifampicin-resistance M. tuberculosis (MTB) strains.

Methods

We searched two electronic databases, PubMed and EMBASE, until 26 March 2017 and updated our search on 27 April 2018 and accessed all prevalence studies of MTB strain and their drug susceptibility patterns to rifampicin. The pooled prevalence estimate was determined using random effects model.

Results

We identified 23 studies satisfying the inclusion criteria. The proportion of rifampicin resistance strains was diverged depending on the type of strains, country and Regions. The pooled estimate of rifampicin-resistance strains of MTB for the included studies was 4% (95% CI: 3–5%). In subgroup analysis based on World Health Organization (WHO) Regions, the pooled estimate of rifampicin-resistance strains of MTB was 11% (95% CI: 9–13%) with the Western Pacific Region 24%, Europian Region 10%, South-East Asian Region 6%, African Region 3% and Region of American 1%. Beijing family was the most dominant strain resistance to rifampicin with pooled prevalence of 14% (95% CI: 10–18%). The pooled prevalence of other families, i.e. EAI, T, CAS, MANU, Haarlem, LAM and Ural, was ≤2% for each.

Conclusion

High burden of rifampicin resistance MTB strains was identified in the Western Pacific Region. Of these, Beijing family was predominantly resistance to rifampicin in Western Pacific Region and South-East Asian Region and also spread to European Region and Region of American.

Characterization of Mycobacterium tuberculosis strains in Beijing, China: drug susceptibility phenotypes and Beijing genotype family transmission

Background

The most prevalent strains of Mycobacterium tuberculosis (M.tb) in Beijing belong to the Beijing genotype family. The influence of Beijing genotype prevalence on the development of drug resistance, and the association of infection with Beijing genotype M.tb with population characteristics, in Beijing, however, are still unclear.

Methods

In this retrospective study, 1189 isolates were subjected to drug susceptibility testing (DST) and molecular epidemiological analysis, and differences in the percentage of drug resistance between Beijing and non-Beijing genotype strains were compared. The association between the occurrence of drug resistance and the prevalence of Beijing genotype M.tb was analyzed using statistical methods.

Results

The Beijing genotype family was the dominant genotype (83.3%) among the 1189 M.tb isolates. Beijing genotype M.tb strains were more likely to spread among males [p = 0.018, OR (95% CI):1.127(1.004–1.264)] and people in the 45–64 age group [p = 0.016, OR (95% CI): 1.438 (1.027–2.015)]. On the contrary, non-Beijing genotype M.tb strains were more probably disseminated among the over 65 [p = 0.005, OR (95% CI):0.653 (0.474–0.9)] and non-resident population [p = 0.035, OR (95% CI):1.185(0.985–1.427)]. DST results showed that 849 (71.4%) strains were fully sensitive to first-line drugs, while 340 (28.6%) strains were resistant to at least one drug, and 9% (107/1189) were MDR-TB. The frequency of INH-resistance among Beijing genotype strains was significantly lower than that among non-Beijing genotype strains (p = 0.032). In addition, the Beijing genotype family readily formed clusters.

Conclusions

Our findings indicate that male and middle-aged people were more probably be infected by Beijing genotype M.tb, older people and non-residents were more probably be infected by non-Beijing genotype M.tb. The high percentage of resistance to INH occurring in non-Beijing genotype strains suggested that non-Beijing genotype strains should be given much more interest in Beijing.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3578-7) contains supplementary material, which is available to authorized users.

Wavelet-domain elastic net for clustering on genomes strains

We propose to evaluate genome similarity by combining discrete non-decimated wavelet transform (NDWT) and elastic net. The wavelets represent a signal with levels of detail, that is, hidden components are detected by means of the decomposition of this signal, where each level provides a different characteristic. The main feature of the elastic net is the grouping of correlated variables where the number of predictors is greater than the number of observations. The combination of these two methodologies applied in the clustering analysis of the Mycobacterium tuberculosis genome strains proved very effective, being able to identify clusters at each level of decomposition.

Multifidelity Analysis for Predicting Rare Events in Stochastic Computational Models of Complex Biological Systems

Rare events such as genetic mutations or cell-cell interactions are important contributors to dynamics in complex biological systems, eg, in drug-resistant infections. Computational approaches can help analyze rare events that are difficult to study experimentally. However, analyzing the frequency and dynamics of rare events in computational models can also be challenging due to high computational resource demands, especially for high-fidelity stochastic computational models. To facilitate analysis of rare events in complex biological systems, we present a multifidelity analysis approach that uses medium-fidelity analysis (Monte Carlo simulations) and/or low-fidelity analysis (Markov chain models) to analyze high-fidelity stochastic model results. Medium-fidelity analysis can produce large numbers of possible rare event trajectories for a single high-fidelity model simulation. This allows prediction of both rare event dynamics and probability distributions at much lower frequencies than high-fidelity models. Low-fidelity analysis can calculate probability distributions for rare events over time for any frequency by updating the probabilities of the rare event state space after each discrete event of the high-fidelity model. To validate the approach, we apply multifidelity analysis to a high-fidelity model of tuberculosis disease. We validate the method against high-fidelity model results and illustrate the application of multifidelity analysis in predicting rare event trajectories, performing sensitivity analyses and extrapolating predictions to very low frequencies in complex systems. We believe that our approach will complement ongoing efforts to enable accurate prediction of rare event dynamics in high-fidelity computational models.

Whole genome SNP analysis suggests unique virulence factor differences of the Beijing and Manila families of Mycobacterium tuberculosis found in Hawaii

While tuberculosis (TB) remains a global disease, the WHO estimates that 62% of the incident TB cases in 2016 occurred in the WHO South-East Asia and Western Pacific regions. TB in the Pacific is composed predominantly of two genetic families of Mycobacterium tuberculosis (Mtb): Beijing and Manila. The Manila family is historically under-studied relative to the families that comprise the majority of TB in Europe and North America (e.g. lineage 4), and it remains unclear why this lineage has persisted in Filipino populations despite the predominance of more globally successful Mtb lineages in most of the world. The Beijing family is of particular interest as it is increasingly associated with drug resistance throughout the world. Both of these lineages are important to the State of Hawaii, where they comprise over two-thirds of TB cases. Here, we performed whole genome sequencing on 82 Beijing family, Manila family, and outgroup clinical Mtb isolates from Hawaii to identify lineage-specific SNPs (SNPs found in all isolates from their respective families, and exclusively in those families) in established virulence factor genes. Six non-silent lineage-specific virulence factor SNPs were found in the Beijing family, including mutations in alternative sigma factor sigG and polyketide synthases pks5 and pks7. The Manila family displayed more than eleven non-silent lineage-specific and characteristic virulence factor mutations, including in genes coding for MCE-family protein Mce1B, two mutations in fatty-acid-AMP ligase FadD26, and virulence-regulating transcriptional regulator VirS. This study further identified an ancient clade that shared some virulence factor mutations with the Manila family, and investigated the relationship of those and other “Manila-like” spoligotypes to the Manila family with this SNP dataset. This work identified a set of virulence genes that are worth pursuing to determine potential differences in transmission or virulence displayed by these two Mtb families.

Characteristics of distribution of Mycobacterium tuberculosis lineages in China

The genotyping methods of Mycobacterium tuberculosis would dramatically improve our understanding of the molecular epidemiology of tuberculosis. 3,929 isolates, from a National Survey of Drug-Resistant Tuberculosis in 2007 in China, were successfully genotyped by large sequence polymorphisms and 15 loci variable number tandem repeats. We found that 2,905 (2,905/3,929, 73.9%) cases belonged to Lineage 2, dominated in the east and central regions, 975 cases (975/3,929, 24.8%) were Lineage 4, highly prevailed in the west regions, and 36 and 13 cases were Lineage 3 and Lineage 1, respectively. We also explored the associations between lineages (Lineage 2 vs. Lineage 4) and clinical characteristics by logistic regression. For Lineage 2, the risk factors were Han-ethnicity population and fever. However, for Lineage 4, they were occupation (farmer), and degree of education (non-literate). Fully understanding of the distribution of Mycobacterium tuberculosis lineage and its risk factors would play a critical role in tuberculosis prevention, control, and treatment.