Genomic diversity among Beijing and non-Beijing Mycobacterium tuberculosis isolates from Myanmar

Genome-Wide Study of Drug Resistant Mycobacterium tuberculosis and Its Intra-Host Evolution during Treatment

The emergence of drug resistant Mycobacterium tuberculosis (MTB) strains has become a global public health problem, while, at the same time, there has been development of new antimicrobial agents. The main goals of this study were to determine new variants associated with drug resistance in MTB and to observe which polymorphisms emerge in MTB genomes after anti-tuberculosis treatment. We performed whole-genome sequencing of 152 MTB isolates including 70 isolates as 32 series of pre- and post-treatment MTB. Based on genotypes and phenotypic drug susceptibility, we conducted phylogenetic convergence-based genome-wide association study (GWAS) with streptomycin-, isoniazid-, rifampicin-, ethambutol-, fluoroquinolones-, and aminoglycosides-resistant MTB against susceptible ones. GWAS revealed statistically significant associations of SNPs within Rv2820c, cyp123 and indels in Rv1269c, Rv1907c, Rv1883c, Rv2407, Rv3785 genes with resistant MTB phenotypes. Comparisons of serial isolates showed that treatment induced different patterns of intra-host evolution. We found indels within Rv1435c and ppsA that were not lineage-specific. In addition, Beijing-specific polymorphisms within Rv0036c, Rv0678, Rv3433c, and dop genes were detected in post-treatment isolates. The appearance of Rv3785 frameshift insertion in 2 post-treatment strains compared to pre-treatment was also observed. We propose that the insertion within Rv3785, which was a GWAS hit, might affect cell wall biosynthesis and probably mediates a compensatory mechanism in response to treatment. These results may shed light on the mechanisms of MTB adaptation to chemotherapy and drug resistance formation.

Correlations between drug resistance of Beijing/W lineage clinical isolates of Mycobacterium tuberculosis and sublineages: a 2009-2013 prospective study in Xinjiang province, China

BACKGROUND

The prevalence of drug-resistant tuberculosis (TB) in Xinjiang is higher than in other regions of China, and Beijing/W lineage Mycobacterium tuberculosis (MTB) is the dominant strain of MTB in Xinjiang. However, information on multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB, particularly the correlation between MDR and the Beijing/W lineage and the correlation between drug resistance and the Beijing/W sublineage strains, is limited.

MATERIAL/METHODS

We conducted a prospective study to describe the prevalence of MDR/XDR TB, Beijing/W lineage and sublineage strains in Xinjiang in China from 2009 to 2013. All MTB underwent drug susceptibility testing to the first- and second-line anti-tuberculosis drugs. The Beijing/W lineages and sublineages were detected by large-sequence polymorphisms with polymerase chain reaction.

RESULTS

A total of 410 clinical isolates were identified. The overall percentage of MDR and XDR cases in Xinjiang was 13.2% (54/410) and 13.0% (7/54), respectively. Overall, 9.8% (14/143) of the Beijing lineage MTB were MDR patients, and 15.6% (40/257) of the Non-Beijing lineage MTB were MDR patients. In the 143 Beijing MTB lineages, 11.2% isolates were in sublineage 105, 15.4% isolates were in sublineage 207, 69.2% isolates were in sublineage 181, and 4.2% isolates were in sublineage 150. None of the isolates were detected in sublineage 142. Significant differences between the Beijing/W and non-Beijing/W strains were observed regarding INH and EMB resistance, respectively.

CONCLUSIONS

The prevalence of the MDR TB in Xinjiang remains high and imposes challenges for TB control. Four Beijing/W sublineage isolates were observed in Xinjiang. There was no correlation between MDR and the Beijing/W lineage and no correlation between drug resistance and the Beijing/W sublineage strains. Surveillance of the clinical isolates of MTB is recommended to strengthen the identification of MDR/XDR TB and sublineages of the Beijing/W strains.

Characterisation of pks15/1 in clinical isolates of Mycobacterium tuberculosis from Mexico

Tuberculosis (TB) is an infectocontagious respiratory disease caused by members of the Mycobacterium tuberculosis complex. A 7 base pair (bp) deletion in the locus polyketide synthase (pks)15/1 is described as polymorphic among members of the M. tuberculosis complex, enabling the identification of Euro-American, Indo-Oceanic and Asian lineages. The aim of this study was to characterise this locus in TB isolates from Mexico. One hundred twenty clinical isolates were recovered from the states of Veracruz and Estado de Mexico. We determined the nucleotide sequence of a ± 400 bp fragment of the locus pks15/1, while genotypic characterisation was performed by spoligotyping. One hundred and fifty isolates contained the 7 bp deletion, while five had the wild type locus. Lineages X (22%), LAM (18%) and T (17%) were the most frequent; only three (2%) of the isolates were identified as Beijing and two (1%) EAI-Manila. The wild type pks15/1 locus was observed in all Asian lineage isolates tested. Our results confirm the utility of locus pks15/1 as a molecular marker for identifying Asian lineages of the M. tuberculosis complex. This marker could be of great value in the epidemiological surveillance of TB, especially in countries like Mexico, where the prevalence of such lineages is unknown.

There is no correlation between sublineages and drug resistance of Mycobacterium tuberculosis Beijing/W lineage clinical isolates in Xinjiang, China

The Beijing/W lineage strains are the major prevalent strains in China. The prevalence, mortality and drug-resistant rates of tuberculosis in Xinjiang, Northwestern China are higher than in other parts of the country. Our previous study results showed that the dominant strains of Mycobacterium tuberculosis (MTB) were ‘Beijing/W lineage’ MTB in Xinjiang; those strains had no significant correlation with drug resistance. We investigated whether the prevalence of ‘Beijing/W lineage’ sublineage strains was associated with drug resistance. We collected 478 sputum specimens from patients with pulmonary tuberculosis. Beijing/W strains and their sublineages were identified by distinguishing five specific large sequence polymorphisms, using polymerase chain reaction. All strains were subjected to a drug susceptibility test using the proportion method on Löwenstein–Jensen culture medium. In total, 379 clinical isolates of MTB were isolated and identified, 57·26% of these isolates were identified as Beijing/W strains, of which 11·06% isolates were in sublineage 105, 14·74% isolates in sublineage 207, 69·59% isolates in sublineage 181, and 4·61% isolates in sublineage 150. None of the isolates was in sublineage 142. Our data showed there were four sublineages of Beijing/W isolates in Xinjiang province, China. However, there were no correlations between drug resistance and the sublineages of Beijing/W strains.

Detection of Mycobacterium bovis in bovine and bubaline tissues using nested-PCR for TbD1

In the present study, a nested-PCR system, targeting the TbD1 region, involving the performance of conventional PCR followed by real-time PCR, was developed to detect Mycobacterium bovis in bovine/bubaline tissue homogenates. The sensitivity and specificity of the reactions were assessed with DNA samples extracted from tuberculous and non-tuberculous mycobacteria, as well as other actinomycetales species and DNA samples extracted directly from bovine and bubaline tissue homogenates. In terms of analytical sensitivity, the DNA of M. bovis AN5 was detected up to 1.56 ng with conventional PCR, 97.6 pg with real-time PCR, and 1.53 pg with nested-PCR in the reaction mixture. The nested-PCR exhibited 100% analytical specificity for M. bovis when tested with the DNA of reference strains of environmental mycobacteria and closely-related Actinomycetales. A clinical sensitivity value of 76.0% was detected with tissue samples from animals that exhibited positive results in the comparative intradermal tuberculin test (CITT), as well as from those with lesions compatible with tuberculosis (LCT) that rendered positive cultures. A clinical specificity value of 100% was detected with tissue samples from animals with CITT- results, with no visible lesions (NVL) and negative cultures. No significant differences were found between the nested-PCR and culture in terms of detecting CITT+ animals with LCT or with NVL. No significant differences were recorded in the detection of CITT- animals with NVL. However, nested-PCR detected a significantly higher number of positive animals than the culture in the group of animals exhibiting LCT with no previous records of CITT. The use of the nested-PCR assay to detect M. bovis in tissue homogenates provided a rapid diagnosis of bovine and bubaline tuberculosis.

Somatic maintenance resources in the honeybee worker fat body are distributed to withstand the most life-threatening challenges at each life stage

In a global transcriptome analysis of three natural and three manipulated honeybee worker phenotypes at different ages, we have investigated the distribution of investment in somatic maintenance of the fat body. Gene expression is modulated so that the bees are able to resist the most life-threatening challenges at the actual life stage. Different modes of maintenance and repair are regulated, apparently to meet the environmental challenges most detrimental to survival and reproductive potential for the hive. We observed a broad down-regulation of genomic and cellular maintenance in the short-lived foragers and nurse bees compared to the long-lived winter bees. Our results show that survival and reproduction of the entire hive is given priority over the individual bees, hence supporting the idea of the honeybee society as a superorganism. Our results also fit the disposable soma theory of aging.

Genotypes and drug susceptibility of Mycobacterium tuberculosis Isolates in Shihezi, Xinjiang Province, China

Background

Tuberculosis (TB) remains a major global health problem. To investigate the genotypes of Mycobacterium tuberculosis (MTB) and the distribution of Beijing family strains, molecular epidemiology technologies have been used widely.

Methods

From June 2010 to June 2011, 55 M. tuberculosis isolates from patients with pulmonary TB were studied by Beijing family-specific PCR (detection of the deletion of region of difference 105 [RD105]), and mycobacterial interspersed repetitive units variable number tandem repeat (MIRU-VNTR) analysis. Twenty-four MIRU-VNTR loci defined the genotypes and clustering characteristics of the local strains. All strains were subjected to a drug susceptibility test (DST) by the proportion method on Lowenstein-Jensen (LJ) culture media.

Results

Fifty-five clinical isolates of MTB were collected. Beijing family strains represented 85.5% of the isolates studied. Using 24 loci MIRU-VNTR typing categorized the strains into eight gene groups, 46 genotypes, and seven clusters. 83.6% (46/55) of the isolates belonged to the largest gene group. Thirty-six isolates (65.5%) were susceptible, nineteen (34.5%) were resistant to at least one drug, seven (12.8%) were Multidrug-Resistant Tuberculosis (MDR TB), and two (3.6%) were extremely drug-resistant tuberculosis (XDR-TB).

Conclusion

The results showed there were obvious polymorphisms of VNTRs of MTB clinical strains. Beijing family strains of MTB were predominant in the Shihezi region of Xinjiang province. There was no correlation between the drug-resistance and Beijing family strains of MTB. It is necessary to strengthen the monitoring, treatment, and management of drug-resistance TB in Shihezi region, Xinjiang.

SNP/RD typing of Mycobacterium tuberculosis Beijing strains reveals local and worldwide disseminated clonal complexes

The Beijing strain is one of the most successful genotypes of Mycobacterium tuberculosis worldwide and appears to be highly homogenous according to existing genotyping methods. To type Beijing strains reliably we developed a robust typing scheme using single nucleotide polymorphisms (SNPs) and regions of difference (RDs) derived from whole-genome sequencing data of eight Beijing strains. SNP/RD typing of 259 M. tuberculosis isolates originating from 45 countries worldwide discriminated 27 clonal complexes within the Beijing genotype family. A total of 16 Beijing clonal complexes contained more than one isolate of known origin, of which two clonal complexes were strongly associated with South African origin. The remaining 14 clonal complexes encompassed isolates from different countries. Even highly resolved clonal complexes comprised isolates from distinct geographical sites. Our results suggest that Beijing strains spread globally on multiple occasions and that the tuberculosis epidemic caused by the Beijing genotype is at least partially driven by modern migration patterns. The SNPs and RDs presented in this study will facilitate future molecular epidemiological and phylogenetic studies on Beijing strains.

Concordance of variable-number tandem repeat (VNTR) and large sequence polymorphism (LSP) analyses of Mycobacterium tuberculosis strains

Variable-number tandem repeat (VNTR) and large sequence polymorphism (LSP) analyses were compared to determine whether VNTR analysis was effective for population genetic analysis of Mycobacterium tuberculosis strains. A total of 682 strains, 510 Beijing genotype and 172 non-Beijing genotype strains, were studied. The number of repeats was investigated for 24 VNTR loci: the 15 loci of "optimized miru", the 8 loci of "Beijing option", and 1 locus for "JATA12". Six loci (miru31, Mtub4, QUB4156c, QUB3232, VNTR3820, and VNTR4120) showed significantly different median numbers of repeats in strains belonging to different lineages defined by LSP (P<0.01, Mann-Whitney U test). When a minimum-spanning tree (MST) was reconstructed using these 6 loci, most strains clustered in the expected branches in the MST branches. However, topology of the MST was not congruent with the evolutional hypothesis of M. tuberculosis, indicating that MST analysis using VNTR data should not use for phylogeny of the organism. When the standardized index of association (sI(A)) was calculated using data for the 6 VNTR loci, the value of sI(A) was significantly different from zero (Monte Carlo simulation with 10,000 resamplings) in every lineage, indicating the linkage disequilibrium in different lineage strains of M. tuberculosis. These results were consistent with the hypothesis that clonal evolution of lineages of the organism has occurred. Therefore, the 6 loci identified in this study would be effective for M. tuberculosis population genetic analysis due to their significantly different median numbers of repeat and linkage disequilibrium though VNTR data was not effective for phylogeny of the organism.

Immunoproteomic identification of human T cell antigens of Mycobacterium tuberculosis that differentiate healthy contacts from tuberculosis patients

Identification of Mycobacterium tuberculosis antigens inducing cellular immune responses is required to improve the diagnosis of and vaccine development against tuberculosis. To identify the antigens of M. tuberculosis that differentiated between tuberculosis (TB) patients and healthy contacts based on T cell reactivity, the culture filtrate of in vitro grown M. tuberculosis was fractionated by two-dimensional liquid phase electrophoresis and tested for the ability to stimulate T cells in a whole blood assay. This approach separated the culture filtrate into 350 fractions with sufficient protein quantity (at least 200 microg of protein) for mass spectrometry and immunological analyses. High levels of interferon-gamma (IFN-gamma) secretion were induced by 105 fractions in healthy contacts compared with TB patients (p < 0.05). Most interesting was the identification of 10 fractions that specifically induced strong IFN-gamma production in the healthy contact population but not in TB patients. Other immunological measurements showed 42 fractions that induced significant lymphocyte proliferative responses in the healthy contact group compared with the TB patients. The tumor necrosis factor-alpha response for most of the fractions did not significantly differ in the tested groups, and the interleukin-4 response was below the detectable range for all fractions and both study groups. Proteomic characterization of the 105 fractions that induced a significant IFN-gamma response in the healthy contacts compared with the TB patients led to the identification of 59 proteins of which 24 represented potentially novel T cell antigens. Likewise, the protein identification in the 10 healthy "contact-specific fractions" revealed 16 proteins that are key candidates as vaccine or diagnostic targets.

Predominance of Mycobacterium tuberculosis EAI and Beijing lineages in Yangon, Myanmar

Isolates of the Mycobacterium tuberculosis Beijing lineage are associated with high rates of transmission, hypervirulence and drug resistance. The Beijing lineage has been shown to dominate the tuberculosis (TB) epidemic in East Asia; however, the diversity and frequency of M. tuberculosis genotypes from Myanmar are unknown. We present the first comprehensive study describing the M. tuberculosis isolates circulating in Yangon, Myanmar. Thus, 310 isolates from pulmonary TB patients from Yangon, Myanmar, were genotyped by spoligotyping and IS6110-based restriction fragment length polymorphism analysis (IS6110 RFLP). The most frequent lineages observed were the East African-Indian (EAI; 48.4%; n = 150) and Beijing (31.9%; n = 99) lineages. Isolates belonging to the most frequent shared types (STs), ST1 (n = 98; Beijing), ST292 (n = 28; EAI), and ST89 (n = 11; EAI), had >or=75% similarity in their IS6110 patterns. Five of 11 Beijing isolates comprising five clusters with identical IS6110 RFLP patterns could be discriminated by mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) analysis. Of the 150 EAI isolates, 40 isolates (26.7%) had only one IS6110 copy, and 17 of these isolates could be discriminated by MIRU-VNTR analysis. The findings from this study suggest that although there is a predominance of the ancient EAI lineage in Yangon, the TB epidemic in Yangon is driven by clonal expansion of the ST1 genotype. The Beijing lineage isolates (21.4%) were more likely (P = 0.009) than EAI lineage isolates to be multidrug resistant (MDR) (1.3%; odds ratio, 3.2, adjusted for the patients' history of exposure to anti-TB drugs), suggesting that the spread of MDR Beijing isolates is a major problem in Yangon.