Significance: Mycobacterium tuberculosis Beijing family is a group of globally emerging bacterial strains that are responsible for more than a quarter of the global tuberculosis epidemic. Here, we combine whole-genome sequencing and large-scale genotyping to map the temporal and spatial changes of the genetic diversity within this strain family. We reveal a southern East Asia origin and a parallel evolution of this bacterial genotype with modern humans in East Asia during the last 30,000 years. The recently globally emerged Beijing strains mainly belong to a hypervirulent subtype that most likely has initially been selected for adaption to increased population densities during the agricultural transition in northern China.

Isoniazid (INH) and rifampicin (RIF) are the two most effective drugs in tuberculosis therapy. Understanding the molecular mechanisms of resistance to these two drugs is essential to quickly diagnose multidrug-resistant (MDR) tuberculosis and extensive drug-resistant tuberculosis. Nine clinical Mycobacterium tuberculosis isolates resistant to only INH and RIF and 10 clinical pan-sensitive isolates were included to evaluate the expression of 20 putative drug efflux pump genes and sequence mutations in rpoB (RIF), katG (INH), the inhA promoter (INH), and oxyR-ahpC (INH). Nine and three MDR isolates were induced to overexpress efflux pump genes by INH and RIF, respectively. Eight and two efflux pump genes were induced to overexpress by INH and RIF in MDR isolates, respectively. drrA, drrB, efpA, jefA (Rv2459), mmr, Rv0849, Rv1634, and Rv1250 were overexpressed under INH or RIF stress. Most efflux pump genes were overexpressed under INH stress in a MDR isolates that carried the wild-type katG, inhA, and oxyR-ahpC associated with INH resistance than in those that carried mutations. The expression levels of 11 genes (efpA, Rv0849, Rv1250, P55 (Rv1410c), Rv1634, Rv2994, stp, Rv2459, pstB, drrA, and drrB) without drug inducement were significantly higher (P < 0.05) in nine MDR isolates than in 10 pan-sensitive isolates. In conclusion, efflux pumps may play an important role in INH acquired resistance in MDR M. tuberculosis, especially in those strains having no mutations in genes associated with INH resistance; basal expression levels of some efflux pump genes are higher in MDR isolates than in pan-sensitive isolates and the basal expressional differences may be helpful to diagnose and treat resistant tuberculosis.

Objective: To investigate the mutations within the whole rpoB gene of Mycobacterium tuberculosis and analyze their effects on rifampin (RIF) resistance based on crystal structure. Methods: We sequenced the entire rpoB gene in 175 tuberculosis isolates and quantified their minimum inhibitory concentrations using microplate-based assays. Additionally, the structural interactions between wild-type/mutant RpoB and RIF were also analyzed. Results: Results revealed that a total of 34 mutations distributed across 17 different sites within the whole rpoB gene were identified. Of the 34 mutations, 25 could alter the structural interaction between RpoB and RIF and contribute to RIF resistance. Statistical analysis showed that S450L, H445D, H445Y and H445R mutations were associated with high-level RIF resistance, while D435V was associated with moderate-level RIF resistance. Conclusion: Some mutations within the rpoB gene could affect the interaction between RpoB and RIF and thus are associated with RIF resistance. These findings could be helpful to design new antibiotics and develop novel diagnostic tools for drug resistance in TB.

Although the Bacillus Calmette-Guérin (BCG) vaccine against tuberculosis (TB) has been available for more than 75 years, one third of the world's population is still infected with Mycobacterium tuberculosis and approximately 2 million people die of TB every year. To reduce this immense TB burden, a clearer understanding of the functional genes underlying the action of BCG and the development of new vaccines are urgently needed. Comparative genomic analysis of 19 M. tuberculosis complex strains showed that BCG strains underwent repeated human manipulation, had higher region of deletion rates than those of natural M. tuberculosis strains, and lost several essential components such as T-cell epitopes. A total of 188 BCG strain T-cell epitopes were lost to various degrees. The non-virulent BCG Tokyo strain, which has the largest number of T-cell epitopes (359), lost 124. Here we propose that BCG strain protection variability results from different epitopes. This study is the first to present BCG as a model organism for genetics research. BCG strains have a very well-documented history and now detailed genome information. Genome comparison revealed the selection process of BCG strains under human manipulation (1908–1966). Our results revealed the cause of BCG vaccine strain protection variability at the genome level and supported the hypothesis that the restoration of lost BCG Tokyo epitopes is a useful future vaccine development strategy. Furthermore, these detailed BCG vaccine genome investigation results will be useful in microbial genetics, microbial engineering and other research fields.

Ethambutol (EMB) plays a pivotal role in the chemotherapy of drug-resistant tuberculosis (TB), including multidrug-resistant tuberculosis (MDR-TB). Resistance to EMB is considered to be caused by mutations in the embCAB operon ( embC , embA , and embB ). In this study, we analyzed the embCAB mutations among 139 MDR-TB isolates from China and found a possible association between embCAB operon mutation and EMB resistance. Our data indicate that 56.8% of MDR-TB isolates are resistant to EMB, and 82.2% of EMB-resistant isolates belong to the Beijing family. Overall, 110 (79.1%) MDR-TB isolates had at least one mutation in the embCAB operon. The majority of mutations were present in the embB gene and the embA upstream region, which also displayed significant correlations with EMB resistance. The most common mutations occurred at codon 306 in embB ( embB 306), followed by embB 406, embA (−16), and embB 497. Mutations at embB 306 were associated with EMB resistance. DNA sequencing of embB 306–497 was the best strategy for detecting EMB resistance, with 89.9% sensitivity, 58.3% specificity, and 76.3% accuracy. Additionally, embB 306 had limited value as a candidate predictor for EMB resistance among MDR-TB infections in China.

We genotyped 102 Borrelia burgdorferi sensu lato strains isolated from ticks, animals, and patients in 11 provinces in China by PCR–restriction fragment length polymorphism (PCR-RFLP) amplification of 5S ( rrf )-23S ( rrl ) rRNA gene spacer amplicons and multilocus sequence analysis (MLSA). The results showed that Borrelia garinii was the main genotype in China (65/102) and that it was distributed mainly in northern China. Borrelia afzelii was the second most frequently found species (22/102), and it was distributed in both northern and southern China. All Borrelia valaisiana strains were isolated from Guizhou Province. Additionally, one B. burgdorferi strain was isolated from Hunan Province. Our results show the diversity and wide distribution of B. burgdorferi sensu lato in China.