Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis

Scanning of genetic diversity of evolutionarily sequential Mycobacterium tuberculosis Beijing family strains based on genome wide analysis

The Beijing family is an endemic lineage of Mycobacterium tuberculosis in eastern Asia. In Japan, five evolutionarily sequential sublineages composing the lineage have predominated. Comparative genomic sequencing based on a microarray technique was conducted for five representative strains of those respective sublineages. Results revealed approximately 200 point mutations specific to each strain. Subsequently, to investigate the genetic diversity of each sublineage, we analysed the phylogenetic divergence of 103 domestic strains belonging to them using genetic markers derived from the mutation information. Results show that the five sublineages have comprised smaller lineages which had diverged at various points. The smaller sub-sublineages have emerged with respective bottlenecks, which are reflected in the excessive monophyletic evolution of the species. Our data provide necessary information to grasp a comprehensive picture of genetic diversity of the lineage constructed in its evolution.

Effect of mutation and genetic background on drug resistance in Mycobacterium tuberculosis

Bacterial factors may contribute to the global emergence and spread of drug-resistant tuberculosis (TB). Only a few studies have reported on the interactions between different bacterial factors. We studied drug-resistant Mycobacterium tuberculosis isolates from a nationwide study conducted from 2000 to 2008 in Switzerland. We determined quantitative drug resistance levels of first-line drugs by using Bactec MGIT-960 and drug resistance genotypes by sequencing the hot-spot regions of the relevant genes. We determined recent transmission by molecular methods and collected clinical data. Overall, we analyzed 158 isolates that were resistant to isoniazid, rifampin, or ethambutol, 48 (30.4%) of which were multidrug resistant. Among 154 isoniazid-resistant strains, katG mutations were associated with high-level and inhA promoter mutations with low-level drug resistance. Only katG(S315T) (65.6% of all isoniazid-resistant strains) and inhA promoter -15C/T (22.7%) were found in molecular clusters. M. tuberculosis lineage 2 (includes Beijing genotype) was associated with any drug resistance (adjusted odds ratio [OR], 3.0; 95% confidence interval [CI], 1.7 to 5.6; P < 0.0001). Lineage 1 was associated with inhA promoter -15C/T mutations (OR, 6.4; 95% CI, 2.0 to 20.7; P = 0.002). We found that the genetic strain background influences the level of isoniazid resistance conveyed by particular mutations (interaction tests of drug resistance mutations across all lineages; P < 0.0001). In conclusion, M. tuberculosis drug resistance mutations were associated with various levels of drug resistance and transmission, and M. tuberculosis lineages were associated with particular drug resistance-conferring mutations and phenotypic drug resistance. Our study also supports a role for epistatic interactions between different drug resistance mutations and strain genetic backgrounds in M. tuberculosis drug resistance.

Multidrug-resistant tuberculosis, People's Republic of China, 2007-2009

Early detection, effective treatment, and infection control measures are needed to reduce transmission.

We conducted a case–control study to investigate risk factors for multidrug-resistant tuberculosis (MDR TB) in the People’s Republic of China. Genotyping analysis was used to estimate the percentage of cases from recent transmission among 100 MDR TB case-patients hospitalized during April 2007–July 2009. Molecular subtyping of isolates showed that 41% of MDR TB strains clustered. Beijing genotype was found in 94% of the MDR TB isolates and 79% of the pan-susceptible isolates. In multivariate analysis, MDR TB was independently associated with Beijing genotype, retreatment for TB, symptoms lasting >3 months before first evaluation at the hospital, lack of health insurance, and being a farmer (vs. being a student). MDR TB was associated with Beijing genotype and lower socioeconomic status. A large percentage of MDR TB cases seemed to result from recent transmission. Early detection, effective treatment, and infection control measures for MDR TB are needed to reduce transmission.

Combination of single nucleotide polymorphism and variable-number tandem repeats for genotyping a homogenous population of Mycobacterium tuberculosis Beijing strains in China

The standard 15- and 24-locus variable-number tandem repeat (VNTR) genotyping methods have demonstrated adequate discriminatory power and a small homoplasy effect for tracing tuberculosis (TB) transmission and predicting Mycobacterium tuberculosis lineages in European and North American countries. However, its validity for the definition of transmission in homogenous M. tuberculosis populations in settings with high TB burdens has been questioned. Here, we genotyped a population-based collection of 191 Beijing strains based on standard 15-locus VNTR (VNTR-15) and 8 single nucleotide polymorphisms (SNPs) in Shanghai, China. Limited discriminatory power and high rates of VNTR homoplasy were observed in the homogenous population of evolutionarily "modern" Beijing strains. Additional typing of three hypervariable loci (VNTR3820, VNTR4120, and VNTR3232) was performed for VNTR-15-based clusters. High variations of hypervariable alleles were observed in clusters with inconsistent SNP sublineages. We concluded that SNPs and hypervariable VNTR loci are helpful to enhance the discriminatory power and decrease the VNTR homoplasy effect for defining clusters. We recommend the combination of standard VNTR-15 and SNPs as first-line typing methods and the hypervariable loci for second-line typing of clustered strains for molecular epidemiology studies of homogenous M. tuberculosis populations.

Pre-existing isoniazid resistance, but not the genotype of Mycobacterium tuberculosis drives rifampicin resistance codon preference in vitro

Both the probability of a mutation occurring and the ability of the mutant to persist will influence the distribution of mutants that arise in a population. We studied the interaction of these factors for the in vitro selection of rifampicin (RIF)-resistant mutants of Mycobacterium tuberculosis. We characterised two series of spontaneous RIF-resistant in vitro mutants from isoniazid (INH)-sensitive and -resistant laboratory strains and clinical isolates, representing various M. tuberculosis genotypes. The first series were selected from multiple parallel 1 ml cultures and the second from single 10 ml cultures. RIF-resistant mutants were screened by Multiplex Ligation-dependent Probe Amplification (MLPA) or by sequencing the rpoB gene. For all strains the mutation rate for RIF resistance was determined with a fluctuation assay. The most striking observation was a shift towards rpoB-S531L (TCG→TTG) mutations in a panel of laboratory-generated INH-resistant mutants selected from the 10-ml cultures (p<0.001). All tested strains showed similar mutation rates (1.33×10⁻⁸ to 2.49×10⁻⁷) except one of the laboratory-generated INH mutants with a mutation rate measured at 5.71×10⁻⁷, more than 10 times higher than that of the INH susceptible parental strain (5.46–7.44×10⁻⁸). No significant, systematic difference in the spectrum of rpoB-mutations between strains of different genotypes was observed. The dramatic shift towards rpoB-S531L in our INH-resistant laboratory mutants suggests that the relative fitness of resistant mutants can dramatically impact the distribution of (subsequent) mutations that accumulate in a M. tuberculosis population, at least in vitro. We conclude that, against specific genetic backgrounds, certain resistance mutations are particularly likely to spread. Molecular screening for these (combinations of) mutations in clinical isolates could rapidly identify these particular pathogenic strains. We therefore recommend that isolates are screened for the distribution of resistance mutations, especially in regions that are highly endemic for (multi)drug resistant tuberculosis.

Investigation on Mycobacterium tuberculosis diversity in China and the origin of the Beijing clade

Background

Investigation of the genetic diversity of Mycobacterium tuberculosis in China has shown that Beijing genotype strains play a dominant role in the tuberculosis (TB) epidemic. In order to examine the strain diversity in the whole country, and to study the evolutionary development of Beijing strains, we sought to genotype a large collection of isolates using different methods.

Methodology/Principal Findings

We applied a 15-loci VNTR typing analysis on 1,586 isolates from the Beijing municipality and 12 Chinese provinces or autonomous regions. The data was compared to that of 900 isolates from various other worldwide geographic regions outside of China. A total of 1,162/1,586 (73.2%) of the isolates, distributed into 472 VNTR types, were found to belong to the Beijing genotype family and this represented 56 to 94% of the isolates in each of the localizations. VNTR typing revealed that the majority of the non-Beijing isolates fall into two genotype families, which represented 17% of the total number of isolates, and seem largely restricted to China. A small number of East African Indian genotype strains was also observed in this collection. Ancient Beijing strains with an intact region of difference (RD) 181, as well as strains presumably resembling ancestors of the whole Beijing genotype family, were mainly found in the Guangxi autonomous region.

Conclusions/Significance

This is the largest M. tuberculosis VNTR-based genotyping study performed in China to date. The high percentage of Beijing isolates in the whole country and the presence in the South of strains representing early branching points may be an indication that the Beijing lineage originated from China, probably in the Guangxi region. Two modern lineages are shown here to represent the majority of non-Beijing Chinese isolates. The observed geographic distribution of the different lineages within China suggests that natural frontiers are major factors in their diffusion.

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.

DNA fingerprinting of Mycobacterium tuberculosis: from phage typing to whole-genome sequencing

Current typing methods for Mycobacterium tuberculosis complex evolved from simple phenotypic approaches like phage typing and drug susceptibility profiling to DNA-based strain typing methods, such as IS6110-restriction fragment length polymorphisms (RFLP) and variable number of tandem repeats (VNTR) typing. Examples of the usefulness of molecular typing are source case finding and epidemiological linkage of tuberculosis (TB) cases, international transmission of MDR/XDR-TB, the discrimination between endogenous reactivation and exogenous re-infection as a cause of relapses after curative treatment of tuberculosis, the evidence of multiple M. tuberculosis infections, and the disclosure of laboratory cross-contaminations. Simultaneously, phylogenetic analyses were developed based on single nucleotide polymorphisms (SNPs), genomic deletions usually referred to as regions of difference (RDs) and spoligotyping which served both strain typing and phylogenetic analysis. National and international initiatives that rely on the application of these typing methods have brought significant insight into the molecular epidemiology of tuberculosis. However, current DNA fingerprinting methods have important limitations. They can often not distinguish between genetically closely related strains and the turn-over of these markers is variable. Moreover, the suitability of most DNA typing methods for phylogenetic reconstruction is limited as they show a high propensity of convergent evolution or misinfer genetic distances. In order to fully explore the possibilities of genotyping in the molecular epidemiology of tuberculosis and to study the phylogeny of the causative bacteria reliably, the application of whole-genome sequencing (WGS) analysis for all M. tuberculosis isolates is the optimal, although currently still a costly solution. In the last years WGS for typing of pathogens has been explored and yielded important additional information on strain diversity in comparison to the classical DNA typing methods. With the ongoing cost reduction of DNA sequencing it is possible that WGS will become the sole diagnostic tool in the secondary laboratory diagnosis of tuberculosis for identification, drug susceptibility testing and genetic characterization.

"Pseudo-Beijing": evidence for convergent evolution in the direct repeat region of Mycobacterium tuberculosis

Background

Mycobacterium tuberculosis has a global population structure consisting of six main phylogenetic lineages associated with specific geographic regions and human populations. One particular M. tuberculosis genotype known as “Beijing” has repeatedly been associated with drug resistance and has been emerging in some parts of the world. “Beijing” strains are traditionally defined based on a characteristic spoligotyping pattern. We used three alternative genotyping techniques to revisit the phylogenetic classification of M. tuberculosis complex (MTBC) strains exhibiting the typical “Beijing” spoligotyping pattern.

Methods and Findings

MTBC strains were obtained from an ongoing molecular epidemiological study in Switzerland and Nepal. MTBC genotyping was performed based on SNPs, genomic deletions, and 24-loci MIRU-VNTR. We identified three MTBC strains from patients originating from Tibet, Portugal and Nepal which exhibited a spoligotyping patterns identical to the classical Beijing signature. However, based on three alternative molecular markers, these strains were assigned to Lineage 3 (also known as Delhi/CAS) rather than to Lineage 2 (also known as East-Asian lineage). Sequencing of the RD207 in one of these strains showed that the deletion responsible for this “Pseudo-Beijing” spoligotype was about 1,000 base pairs smaller than the usual deletion of RD207 in classical “Beijing” strains, which is consistent with an evolutionarily independent deletion event in the direct repeat (DR) region of MTBC.

Conclusions

We provide an example of convergent evolution in the DR locus of MTBC, and highlight the limitation of using spoligotypes for strain classification. Our results indicate that a proportion of “Beijing” strains may have been misclassified in the past. Markers that are more phylogenetically robust should be used when exploring strain-specific differences in experimental or clinical phenotypes.

Perspective on sequence evolution of microsatellite locus (CCG)n in Rv0050 gene from Mycobacterium tuberculosis

Background

The mycobacterial genome is inclined to polymerase slippage and a high mutation rate in microsatellite regions due to high GC content and absence of a mismatch repair system. However, the exact molecular mechanisms underlying microsatellite variation have not been fully elucidated. Here, we investigated mutation events in the hyper-variable trinucleotide microsatellite locus MML0050 located in the Rv0050 gene of W-Beijing and non-W-Beijing Mycobacterium tuberculosis strains in order to gain insight into the genomic structure and activity of repeated regions.

Results

Size analysis indicated the presence of five alleles that differed in length by three base pairs. Moreover, nucleotide gains occurred more frequently than loses in this trinucleotide microsatellite. Mutation frequency was not completely related with the total length, though the relative frequency in the longest allele was remarkably higher than that in the shortest. Sequence analysis was able to detect seven alleles and revealed that point mutations enhanced the level of locus variation. Introduction of an interruptive motif correlated with the total allele length and genetic lineage, rather than the length of the longest stretch of perfect repeats. Finally, the level of locus variation was drastically different between the two genetic lineages.

Conclusion

The Rv0050 locus encodes the bifunctional penicillin-binding protein ponA1 and is essential to mycobacterial survival. Our investigations of this particularly dynamic genomic region provide insights into the overall mode of microsatellite evolution. Specifically, replication slippage was implicated in the mutational process of this microsatellite and a sequence-based genetic analysis was necessary to determine that point mutation events acted to maintain microsatellite size integrity while providing genomic diversity.

Presence of RD149 deletions in M. tuberculosis Central Asian Strain 1 isolates affect growth and TNFα induction in THP-1 monocytes

Central Asian Strain 1 (CAS1) is the prevalent Mycobacterium tuberculosis genogroup in South Asia. CAS1 strains carry deletions in RD149 and RD152 regions. Significance of these deletions is as yet unknown. We compared CAS1 strains with RD149 and concurrent RD149-RD152 deletions with CAS1 strains without deletions and with the laboratory reference strain, M. tuberculosis H37Rv for growth and for induction of TNFα, IL6, CCL2 and IL10 in THP-1 cells. Growth of CAS1 strains with deletions was slower in broth (RD149; p = 0.024 and RD149-RD152; p = 0.025) than that of strains without deletions. CAS1 strains with RD149 deletion strains further showed reduced intracellular growth (p = 0.013) in THP-1 cells as compared with strains without deletions, and also as compared with H37Rv (p = 0.007) and with CAS1 RD149-RD152 deletion strains (p = 0.029). All CAS1 strains induced higher levels of TNFα and IL10 secretion in THP-1 cells than H37Rv. Additionally, CAS1 strains with RD149 deletions induced more TNFα secretion than those without deletions (p = 0.013). CAS1 RD149 deletion strains from extrapulmonary sources showed more rapid growth and induced lower levels of TNFα and IL6 secretion in THP-1 cells than isolates from pulmonary sources. This data suggests that presence of RD149 reduces growth and increases the induction of TNFα in host cells by CAS1 strains. Differences observed for extrapulmonary strains may indicate an adaptation which increases potential for dissemination and tropism outside the lung. Overall, we hypothesise that RD149 deletions generate genetic diversity within strains and impact interactions of CAS1 strains with host cells with important clinical consequences.

Mycobacterium bovis BCG-mediated protection against W-Beijing strains of Mycobacterium tuberculosis is diminished concomitant with the emergence of regulatory T cells

Despite issues relating to variable efficacy in the past, the Mycobacterium bovis BCG vaccine remains the basis for new-generation recombinant vaccines currently in clinical trials. To date, vaccines have been tested mostly against laboratory strains and not against the newly emerging clinical strains. In this study, we evaluated the ability of BCG Pasteur to protect mice from aerosol infections with two highly virulent W-Beijing clinical strains, HN878 and SA161. In a conventional 30-day protection assay, BCG was highly protective against both strains, but by day 60 of the assay, this protection was diminished. Histological examination of the lungs of vaccinated animals showed reduced lung consolidation and smaller and more-organized granulomas in the vaccinated mice after 30 days, but in both cases, these tissues demonstrated worsening pathology over time. Effector T cell responses were increased in the vaccinated mice infected with HN878, but these diminished in number after day 30 of the infections concomitant with increased CD4(+) Foxp3(+) T cells in the lungs, draining lymph nodes, and the spleen. Given the concomitant decrease in effector immunity and continued expansion of regulatory Foxp3(+) cells observed here, it is reasonable to hypothesize that downregulation of effector immunity by these cells may be a serious impediment to the efficacy of BCG-based vaccines.

A novel method for the rapid and prospective identification of Beijing Mycobacterium tuberculosis strains by high-resolution melting analysis

Genotypic analysis of Mycobacterium tuberculosis (MTB) has enabled the definition of several lineages. The Beijing family, which is considered highly virulent and transmissible, has been associated with resistance in certain settings and involved in severe outbreaks, making it one of the most closely-monitored lineages. Therefore, rapid prospective identification of Beijing MTB strains could be relevant. In the present study, we evaluate a real-time PCR followed by high-resolution melting (HRM) based on the identification of a single nucleotide polymorphism (SNP) in the Rv2629 gene which defines Beijing lineage (A191C for Beijing genotype and A191A for non-Beijing genotype). This combined methodology efficiently differentiated Beijing and non-Beijing strains in 100% of the isolates from a collection of reference strains without requiring specific DNA probes. Additionally, HRM was able to assign a Beijing/non-Beijing genotype in 90.9% of the respiratory specimens assayed. Its applicability was tested on a Peruvian sample of circulating MTB strains, in which it identified 10.7% as belonging to the Beijing genotype; this proportion reached 20% in the North Lima area. HRM analysis of the A191C SNP is a rapid, reliable, and sensitive method for the efficient prospective survey of high-risk Beijing MTB strains, even in developing settings where MTB culture is often not available.

Models to understand the population-level impact of mixed strain M. tuberculosis infections

Over the past decade, numerous studies have identified tuberculosis patients in whom more than one distinct strain of Mycobacterium tuberculosis is present. While it has been shown that these mixed strain infections can reduce the probability of treatment success for individuals simultaneously harboring both drug-sensitive and drug-resistant strains, it is not yet known if and how this phenomenon impacts the long-term dynamics for tuberculosis within communities. Strain-specific differences in immunogenicity and associations with drug resistance suggest that a better understanding of how strains compete within hosts will be necessary to project the effects of mixed strain infections on the future burden of drug-sensitive and drug-resistant tuberculosis. In this paper, we develop a modeling framework that allows us to investigate mechanisms of strain competition within hosts and to assess the long-term effects of such competition on the ecology of strains in a population. These models permit us to systematically evaluate the importance of unknown parameters and to suggest priority areas for future experimental research. Despite the current scarcity of data to inform the values of several model parameters, we are able to draw important qualitative conclusions from this work. We find that mixed strain infections may promote the coexistence of drug-sensitive and drug-resistant strains in two ways. First, mixed strain infections allow a strain with a lower basic reproductive number to persist in a population where it would otherwise be outcompeted if has competitive advantages within a co-infected host. Second, some individuals progressing to phenotypically drug-sensitive tuberculosis from a state of mixed drug-sensitive and drug-resistant infection may retain small subpopulations of drug-resistant bacteria that can flourish once the host is treated with antibiotics. We propose that these types of mixed infections, by increasing the ability of low fitness drug-resistant strains to persist, may provide opportunities for compensatory mutations to accumulate and for relatively fit, highly drug-resistant strains of M. tuberculosis to emerge.

Genetic diversity of the Mycobacterium tuberculosis Beijing family based on IS6110, SNP, LSP and VNTR profiles from Thailand

The Beijing family of Mycobacterium tuberculosis (MTB) has been shown to be highly virulent and is associated with tuberculosis outbreaks and drug resistance. To explore the evolutionary pathway of this genotype, 165 epidemiologically unlinked Beijing MTB isolates from Thailand were characterized based on IS6110, variable number of tandem repeats (VNTRs), region of differences (RDs) and single nucleotide polymorphisms (SNPs). Analysis based on VNTRs showed that 7 VNTRs were phylogenetically informative loci. Three of the seven VNTRs, 4156, 1982 and 1955, were promising markers for defining ancestral and modern Beijing sublineages. The evolutionary history of Beijing sublineages based on VNTR profiles showed dynamic changes in the repeat number, which also correlated with specific SNP STs. We constructed a refined, comprehensive phylogenetic tree based on multiple genetic markers from the Beijing strains from Thailand and combined the data with previous studies that were performed in South Africa and Japan. The phylogenetic tree showed how the modern Beijing and W-Beijing strains evolved from ancestral Beijing strains. The points at which the RD150 deletion occurred repeatedly were also demonstrated. The evolutionary scheme showed that most events were concordant with the SNP STs. Nevertheless, some SNPs used for sub-typing the Beijing family of MTB may not be irreversible, unique events.

Molecular modeling and in silico characterization of Mycobacterium tuberculosis TlyA: possible misannotation of this tubercle bacilli-hemolysin

Background

The TlyA protein has a controversial function as a virulence factor in Mycobacterium tuberculosis (M. tuberculosis). At present, its dual activity as hemolysin and RNA methyltransferase in M. tuberculosis has been indirectly proposed based on in vitro results. There is no evidence however for TlyA relevance in the survival of tubercle bacilli inside host cells or whether both activities are functionally linked. A thorough analysis of structure prediction for this mycobacterial protein in this study shows the need for reevaluating TlyA's function in virulence.

Results

Bioinformatics analysis of TlyA identified a ribosomal protein binding domain (S4 domain), located between residues 5 and 68 as well as an FtsJ-like methyltranferase domain encompassing residues 62 and 247, all of which have been previously described in translation machinery-associated proteins. Subcellular localization prediction showed that TlyA lacks a signal peptide and its hydrophobicity profile showed no evidence of transmembrane helices. These findings suggested that it may not be attached to the membrane, which is consistent with a cytoplasmic localization. Three-dimensional modeling of TlyA showed a consensus structure, having a common core formed by a six-stranded β-sheet between two α-helix layers, which is consistent with an RNA methyltransferase structure. Phylogenetic analyses showed high conservation of the tlyA gene among Mycobacterium species. Additionally, the nucleotide substitution rates suggested purifying selection during tlyA gene evolution and the absence of a common ancestor between TlyA proteins and bacterial pore-forming proteins.

Conclusion

Altogether, our manual in silico curation suggested that TlyA is involved in ribosomal biogenesis and that there is a functional annotation error regarding this protein family in several microbial and plant genomes, including the M. tuberculosis genome.

Genotypes and characteristics of clustering and drug susceptibility of Mycobacterium tuberculosis isolates collected in Heilongjiang Province, China

For the last decade China has occupied second place, after India, among the top five countries with high burdens of tuberculosis (TB). Heilongjiang Province is located in northeastern China. The prevalence of drug-resistant TB in Heilongjiang Province is higher than the average level in China. To determine the transmission characteristics of Mycobacterium tuberculosis strains isolated in this area and their genetic relationships, especially among the Beijing family strains, we investigated their genotypes. From May 2007 to October 2008, 200 M. tuberculosis isolates from patients presenting pulmonary TB were analyzed by molecular typing using PCR-based methods: spacer-oligonucleotide typing (spoligotyping), Beijing family-specific PCR (detection of the deletion of region of difference 105 [RD105]), and mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) analysis. Different combinations of MIRU-VNTR loci were evaluated to define the genotypes and clustering characteristics of the local strains. We found that Beijing family strains represented 89.5% of the isolates studied. However, the rates of multidrug-resistant (MDR) M. tuberculosis among Beijing and non-Beijing family strains were not statistically different. The 15-locus set is considered the optimal MIRU-VNTR locus combination for analyzing the M. tuberculosis strains epidemic in this area, while the 10-locus set is an ideal set for first-line molecular typing. We found that the clustering rate of all the M. tuberculosis isolates analyzed was 10.0% using the 15-locus set typing. We conclude that the Beijing family genotype is predominant and that highly epidemic TB and MDR TB are less likely associated with the active transmission of M. tuberculosis in the study area.

Sublineages of Beijing Strain of Mycobacterium tuberculosis in Sri Lanka

Strains of the Beijing/W genotype of Mycobacterium tuberculosis have been responsible for large outbreaks of tuberculosis around the world, sometimes involving multi-drug resistance. It has been shown that more recently evolved Beijing sublineages are prone to cause outbreaks. Furthermore Beijing is the single predominant cluster in Sri Lanka. The present study identifies that recently evolved sublineages of Beijing strains are present in the study population. The majority of Beijing isolates (92.85%) were pan-susceptible. However, these findings may have important implications for the control and prevention of tuberculosis in Sri Lanka.

African 2, a clonal complex of Mycobacterium bovis epidemiologically important in East Africa

We have identified a clonal complex of Mycobacterium bovis isolated at high frequency from cattle in Uganda, Burundi, Tanzania, and Ethiopia. We have named this related group of M. bovis strains the African 2 (Af2) clonal complex of M. bovis. Af2 strains are defined by a specific chromosomal deletion (RDAf2) and can be identified by the absence of spacers 3 to 7 in their spoligotype patterns. Deletion analysis of M. bovis isolates from Algeria, Mali, Chad, Nigeria, Cameroon, South Africa, and Mozambique did not identify any strains of the Af2 clonal complex, suggesting that this clonal complex of M. bovis is localized in East Africa. The specific spoligotype pattern of the Af2 clonal complex was rarely identified among isolates from outside Africa, and the few isolates that were found and tested were intact at the RDAf2 locus. We conclude that the Af2 clonal complex is localized to cattle in East Africa. We found that strains of the Af2 clonal complex of M. bovis have, in general, four or more copies of the insertion sequence IS6110, in contrast to the majority of M. bovis strains isolated from cattle, which are thought to carry only one or a few copies.

Mutations in the regulatory network underlie the recent clonal expansion of a dominant subclone of the Mycobacterium tuberculosis Beijing genotype

The Beijing genotype family is an epidemiologically important sub-group of Mycobacterium tuberculosis. It has been suggested that the high frequency of the Beijing isolates in some areas could be explained by selective advantages. Some evidence suggests that the emerging and most frequently isolated "Typical Beijing" lineage has the ability to circumvent BCG-induced immunity. To investigate the phylogeny of the Beijing genotype of M. tuberculosis, the genome of six Beijing strains from three different countries was sequenced with next-generation sequencing. The phylogeny of these strains was established using single nucleotide polymorphisms (SNPs). The three Typical Beijing strains clustered very tightly in the Beijing phylogeny suggesting that Typical Beijing strains represent a monophyletic lineage and resulted from recent diversification. Typing of 150 M. tuberculosis strains with a subset of the SNPs and comparison of the IS6110 restriction-fragment length polymorphism (RFLP) patterns of these strains to a database of 1522 Beijing RFLP patterns revealed that about 80% of all Beijing strains belong to the Typical Beijing subclone, which indicates clonal expansion. To identify the genomic changes that are characteristic for all Typical Beijing strains and to reconstruct their most recent common ancestor, the presence of SNPs was assayed in other Beijing strains. We identified 51 SNPs that define the minimal set of polymorphisms for all Typical Beijing strains. Nonsynonymous polymorphisms in genes coding for the regulatory network were over-represented in this set of mutations. We suggest that alterations in the response to environmental signals may have enabled Typical Beijing strains to develop the emerging phenotype.

Phylogeny of Mycobacterium tuberculosis Beijing strains constructed from polymorphisms in genes involved in DNA replication, recombination and repair

Background

The Beijing family is a successful group of M. tuberculosis strains, often associated with drug resistance and widely distributed throughout the world. Polymorphic genetic markers have been used to type particular M. tuberculosis strains. We recently identified a group of polymorphic DNA repair replication and recombination (3R) genes. It was shown that evolution of M. tuberculosis complex strains can be studied using 3R SNPs and a high-resolution tool for strain discrimination was developed. Here we investigated the genetic diversity and propose a phylogeny for Beijing strains by analyzing polymorphisms in 3R genes.

Methodology/Principal Findings

A group of 3R genes was sequenced in a collection of Beijing strains from different geographic origins. Sequence analysis and comparison with the ones of non-Beijing strains identified several SNPs. These SNPs were used to type a larger collection of Beijing strains and allowed identification of 26 different sequence types for which a phylogeny was constructed. Phylogenetic relationships established by sequence types were in agreement with evolutionary pathways suggested by other genetic markers, such as Large Sequence Polymorphisms (LSPs). A recent Beijing genotype (Bmyc10), which included 60% of strains from distinct parts of the world, appeared to be predominant.

Conclusions/Significance

We found SNPs in 3R genes associated with the Beijing family, which enabled discrimination of different groups and the proposal of a phylogeny. The Beijing family can be divided into different groups characterized by particular genetic polymorphisms that may reflect pathogenic features. These SNPs are new, potential genetic markers that may contribute to better understand the success of the Beijing family.

Mixed infection with Beijing and non-Beijing strains in pulmonary tuberculosis in Taiwan: prevalence, risk factors, and dominant strain

Patients with pulmonary tuberculosis (TB) can be simultaneously infected with different strains of Mycobacterium tuberculosis (mixed infection). We investigated the prevalence and risk factors of mixed infection by Beijing and non-Beijing strains in pulmonary TB patients in Taiwan. We developed a quantitative PCR method to simultaneously detect the presence of Beijing and non-Beijing strains. A total of 868 pretreatment samples (from 868 patients), including 563 sputum samples smear-positive for acid-fast bacilli and 305 liquid medium samples culture-positive for mycobacteria, were tested. Medical records of patients with culture-confirmed pulmonary TB were reviewed. The detection limit of our quantitative PCR method was five copies of target sequences. With mycobacterial culture result as the reference standard, the sensitivity and specificity of our quantitative PCR method were 95% and 98%, respectively. M. tuberculosis strains were isolated in 466 samples, of which 231 (49.6%) were infected with a Beijing strain. Another 14 patients (3.0%) had mixed infection, with the Beijing strain being the dominant strain in 13 (93%). Age <25 years with pulmonary cavities was associated with mixed infection. In patients infected with non-Beijing strains, the bacterial load of non-Beijing strains was lower among those with mixed infection than among those without. Our quantitative PCR method was accurate in detecting Beijing and non-Beijing strains in smear-positive sputum and culture-positive liquid medium samples. Mixed infection was present in pulmonary TB patients (3.0%), especially in those aged <25 years with pulmonary cavities. Beijing strains seem to be more dominant than non-Beijing strains in patients with mixed infection.

The non-clonality of drug resistance in Beijing-genotype isolates of Mycobacterium tuberculosis from the Western Cape of South Africa

Background

The Beijing genotype of M. tuberculosis is a virulent strain that is disseminating worldwide and has a strong association with drug resistance. In the Western Cape of South Africa, epidemiological studies have identified the R220 cluster of the Beijing genotype as a major contributor to a recent outbreak of drug-resistant tuberculosis. Although the outbreak is considered to be due to clonal transmission, the relationship among drug resistant isolates has not yet been established.

Results

To better understand the evolution of drug resistance among these strains, 14 drug-resistant clinical isolates of the Beijing genotype were sequenced by whole-genome sequencing, including eight from R220 and six from a more ancestral Beijing cluster, R86, for comparison. While each cluster shares a distinct resistance mutation for isoniazid, mapping of other drug-resistance mutations onto a phylogenetic tree constructed from single nucleotide polymorphisms shows that resistance mutations to many drugs have arisen multiple times independently within each cluster of isolates. Thus, drug resistance among these isolates appears to be acquired, not clonally derived. This observation suggests that, although the Beijing genotype as a whole might have selective advantages enabling its rapid dissemination, the XDR isolates are relatively less fit and do not propagate well. Although it has been hypothesized that the increased frequency of drug resistance in some Beijing lineages might be caused by a mutator phenotype, no significant shift in synonymous substitution patterns is observed in the genomes.

Conclusion

While MDR-TB is spreading by transmission in the Western Cape, our data suggests that further drug resistance (i.e. XDR-TB) at this stage is acquired.

Innate immune response to Mycobacterium tuberculosis Beijing and other genotypes

Background

As a species, Mycobacterium tuberculosis is more diverse than previously thought. In particular, the Beijing family of M. tuberculosis strains is spreading and evoluating throughout the world and this is giving rise to public health concerns. Genetic diversity within this family has recently been delineated further and a specific genotype, called Bmyc10, has been shown to represent over 60% of all Beijing clinical isolates in several parts of the world. How the host immune system senses and responds to various M. tuberculosis strains may profoundly influence clinical outcome and the relative epidemiological success of the different mycobacterial lineages. We hypothesised that the success of the Bmyc10 group may, at least in part, rely upon its ability to alter innate immune responses and the secretion of cytokines and chemokines by host phagocytes.

Methodology/Principal Findings

We infected human macrophages and dendritic cells with a collection of genetically well-defined M. tuberculosis clinical isolates belonging to various mycobacterial families, including Beijing. We analyzed cytokine and chemokine secretion on a semi-global level using antibody arrays allowing the detection of sixty-five immunity-related soluble molecules. Our data indicate that Beijing strains induce significantly less interleukin (IL)-6, tumor necrosis factor (TNF), IL-10 and GRO-α than the H37Rv reference strain, a feature that is variously shared by other modern and ancient M. tuberculosis families and which constitutes a signature of the Beijing family as a whole. However, Beijing strains did not differ relative to each other in their ability to modulate cytokine secretion.

Conclusions/Significance

Our results confirm and expand upon previous reports showing that M. tuberculosis Beijing strains in general are poor in vitro cytokine inducers in human phagocytes. The results suggest that the epidemiological success of the Beijing Bmyc10 is unlikely to rely upon any specific ability of this group of strains to impair anti-mycobacterial innate immunity.

Descriptive proteomic analysis shows protein variability between closely related clinical isolates of Mycobacterium tuberculosis

The use of isobaric tags such as iTRAQ allows the relative and absolute quantification of hundreds of proteins in a single experiment for up to eight different samples. More classical techniques such as 2-DE can offer a complimentary approach for the analysis of complex protein samples. In this study, the proteomes of secreted and cytosolic proteins of genetically closely related strains of Mycobacterium tuberculosis were analyzed. Analysis of 2-D gels afforded 28 spots with variations in protein abundance between strains. These were identified by MS/MS. Meanwhile, a rigorous statistical analysis of iTRAQ data allowed the identification and quantification of 101 and 137 proteins in the secreted and cytosolic fractions, respectively. Interestingly, several differences in protein levels were observed between the closely related strains BE, C28 and H6. Seven proteins related to cell wall and cell processes were more abundant in BE, while enzymes related to metabolic pathways (GltA2, SucC, Gnd1, Eno) presented lower levels in the BE strain. Proteins involved in iron and sulfur acquisition (BfrB, ViuB, TB15.3 and SseC2) were more abundant in C28 and H6. In general, iTRAQ afforded rapid identification of fine differences between protein levels such as those presented between closely related strains. This provides a platform from which the relevance of these differences can be assessed further using complimentary proteomic and biological modeling methods.

Strains of Mycobacterium tuberculosis from western Maharashtra, India, exhibit a high degree of diversity and strain-specific associations with drug resistance, cavitary disease, and treatment failure

We performed spoligotyping of Mycobacterium tuberculosis isolates from 833 systematically sampled pulmonary tuberculosis (TB) patients in urban Mumbai, India (723 patients), and adjacent rural areas in western India (110 patients). The urban cohort consisted of two groups of patients, new cases (646 patients) and first-time treatment failures (77 patients), while only new cases were recruited in the rural areas. The isolates from urban new cases showed 71% clustering, with 168 Manu1, 62 CAS, 22 Beijing, and 30 EAI-5 isolates. The isolates from first-time treatment failures were 69% clustered, with 14 Manu1, 8 CAS, 8 Beijing, and 6 EAI-5 isolates. The proportion of Beijing strains was higher in this group than in urban new cases (odds ratio [OR], 3.29; 95% confidence limit [95% CL], 1.29 to 8.14; P = 0.003). The isolates from rural new cases showed 69% clustering, with 38 Manu1, 7 CAS, and 1 EAI-5 isolate. Beijing was absent in the rural cohort. Manu1 was found to be more common in the rural cohort (OR, 0.67; 95% CL, 0.42 to 1.05; P = 0.06). In total, 71% of isolates were clustered into 58 spoligotypes with 4 predominant strains, Manu1 (26%), CAS (9%), EAI-5 (4%), and Beijing (4%), along with 246 unique spoligotypes. In the isolates from urban new cases, we found Beijing to be associated with multidrug resistance (MDR) (OR, 3.40; 95% CL, 1.20 to 9.62; P = 0.02). CAS was found to be associated with pansensitivity (OR, 1.83; 95% CL, 1.03 to 3.24; P = 0.03) and cavities as seen on chest radiographs (OR, 2.72; 95% CL, 1.34 to 5.53; P = 0.006). We recorded 239 new spoligotypes yet unreported in the global databases, suggesting that the local TB strains exhibit a high degree of diversity.

Massive gene duplication event among clinical isolates of the Mycobacterium tuberculosis W/Beijing family

As part of our effort to uncover the molecular basis for the phenotypic variation among clinical Mycobacterium tuberculosis isolates, we have previously reported that isolates belonging to the W/Beijing lineage constitutively overexpress the DosR-regulated transcriptional program. While generating dosR knockouts in two independent W/Beijing sublineages, we were surprised to discover that they possess two copies of dosR. This dosR amplification is part of a massive genomic duplication spanning 350 kb and encompassing >300 genes. In total, this equates to 8% of the genome being present as two copies. The presence of IS6110 elements at both ends of the region of duplication, and in the novel junction region, suggests that it arose through unequal homologous recombination of sister chromatids at the IS6110 sequences. Analysis of isolates representing the major M. tuberculosis lineages has revealed that the 350-kb duplication is restricted to the most recently evolved sublineages of the W/Beijing family. Within these isolates, the duplication is partly responsible for the constitutive dosR overexpression phenotype. Although the nature of the selection event giving rise to the duplication remains unresolved, its evolution is almost certainly the result of specific selective pressure(s) encountered inside the host. A preliminary in vitro screen has failed to reveal a role of the duplication in conferring resistance to common antitubercular drugs, a trait frequently associated with W/Beijing isolates. Nevertheless, this first description of a genetic remodeling event of this nature for M. tuberculosis further highlights the potential for the evolution of diversity in this important global pathogen.

Phylogeographical particularity of the Mycobacterium tuberculosis Beijing family in South Korea based on international comparison with surrounding countries

To understand the domestic population structure of Mycobacterium tuberculosis clinical isolates in the Republic of Korea, we genotypically analysed 80 isolates obtained from various geographical origins in the country. Of these, 64 (80.0 %) isolates were identified as Beijing family strains. It is particularly interesting that their phylogenetic classification, based on the ancient/modern separation and the presence/absence of the genomic region RD181, revealed a majority of the ancient (RD181+) subfamily in the population. The 15 loci of variable number of tandem repeat(s) of mycobacterial interspersed repetitive units (15-MIRU-VNTR) were also analysed. Combination with the previous VNTR data reported from surrounding countries revealed that the topology of the minimum spanning tree was linked tightly not to the geographical origins of the patients but to the phylogenetic characteristics of the isolates. These results show that the phylogeographical distribution of the M. tuberculosis Beijing family around far-eastern Asia could be estimated using international accumulation and comparison of VNTR genotyping data.

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.

Drug resistant Mycobacterium tuberculosis of the Beijing genotype does not spread in Sweden

BACKGROUND

Drug resistant (DR) and multi-drug resistant (MDR) tuberculosis (TB) is increasing worldwide. In some parts of the world 10% or more of new TB cases are MDR. The Beijing genotype is a distinct genetic lineage of Mycobacterium tuberculosis, which is distributed worldwide, and has caused large outbreaks of MDR-TB. It has been proposed that certain lineages of M. tuberculosis, such as the Beijing lineage, may have specific adaptive advantages. We have investigated the presence and transmission of DR Beijing strains in the Swedish population.

METHODOLOGY/PRINCIPAL FINDINGS

All DR M. tuberculosis complex isolates between 1994 and 2008 were studied. Isolates that were of Beijing genotype were investigated for specific resistance mutations and phylogenetic markers. Seventy (13%) of 536 DR strains were of Beijing genotype. The majority of the patients with Beijing strains were foreign born, and their country of origin reflects the countries where the Beijing genotype is most prevalent. Multidrug-resistance was significantly more common in Beijing strains than in non-Beijing strains. There was a correlation between the Beijing genotype and specific resistance mutations in the katG gene, the mabA-inhA-promotor and the rpoB gene. By a combined use of RD deletions, spoligotyping, IS1547, mutT gene polymorphism and Rv3135 gene analysis the Beijing strains could be divided into 11 genomic sublineages. Of the patients with Beijing strains 28 (41%) were found in altogether 10 clusters (2-5 per cluster), as defined by RFLP IS6110, while 52% of the patients with non-Beijing strains were in clusters. By 24 loci MIRU-VNTR 31 (45%) of the patients with Beijing strains were found in altogether 7 clusters (2-11 per cluster). Contact tracing established possible epidemiological linkage between only two patients with Beijing strains.

CONCLUSIONS/SIGNIFICANCE

Although extensive outbreaks with non-Beijing TB strains have occurred in Sweden, Beijing strains have not taken hold, in spite of the proximity to high prevalence countries such as Russia and the Baltic countries. The Beijing sublineages so far introduced in Sweden may not be adapted to spread in the Scandinavian population.

Characterization of Mycobacterium tuberculosis Beijing isolates from the Mediterranean area

BACKGROUND

The Beijing lineage of Mycobacterium tuberculosis is causing concern due to its global distribution and its involvement in severe outbreaks. Studies focused on this lineage are mainly restricted to geographical settings where its prevalence is high, whereas those in other areas are scarce. In this study, we analyze Beijing isolates in the Mediterranean area, where this lineage is not prevalent and is mainly associated with immigrant cases.

RESULTS

Only 1% (N = 26) of the isolates from two population-based studies in Spain corresponded to Beijing strains, most of which were pan-susceptible and from Peruvian and Ecuadorian patients. Restriction fragment length polymorphism typing with the insertion sequence IS6110 identified three small clusters (2-3 cases). Mycobacterial interspersed repetitive unit-variable number tandem repeat typing (MIRU-15) offered low discriminatory power, requiring the introduction of five additional loci. A selection of the Beijing isolates identified in the Spanish sample, together with a sample of Beijing strains from Italy, to broaden the analysis context in the Mediterranean area, were assayed in an infection model with THP-1 cells. A wide range of intracellular growth rates was observed with only two isolates showing an increased intracellular replication, in both cases associated with contained production of TNF-alpha. No correlation was observed between virulence and the Beijing phylogenetic group, clustered/orphan status, or resistance. The Beijing strain responsible for extensive spread on Gran Canaria Island was also identified in Madrid, but did not lead to secondary cases and did not show high infectivity in the infection model.

CONCLUSIONS

The Beijing lineage in our area is a non-homogeneous family, with only certain highly virulent representatives. The specific characterization of Beijing isolates in different settings could help us to accurately identify the virulent representatives before making general assumptions about this lineage.

Differences among sublineages of the East-Asian lineage of Mycobacterium tuberculosis in genotypic clustering

SETTING

The East-Asian lineage of Mycobacterium tuberculosis is composed of five sublineages, and includes the strains from the Beijing spoligotype family. In some studies these strains were highly pathogenic, although other studies did not support this finding.

OBJECTIVE

To determine if the sublineages of the East-Asian lineage of M. tuberculosis differ in their capacity to cause secondary cases, as assessed by genotypic clustering of isolates.

DESIGN

In a population-based study of 545 patients with M. tuberculosis from the East-Asian lineage in San Francisco, we used DNA-based fingerprinting to identify genotypic clustering, which was compared among the different sublineages defined by large sequence polymorphism.

RESULTS

Strains from sublineage 207 had the highest frequency of genotypic clustering. In the multivariate analysis, only patients born in the United States were associated with clustering.

CONCLUSIONS

We found evidence in a univariate analysis that the different East-Asian sublineages of M. tuberculosis have different frequencies of genotypic clustering. The effect size for this difference was unchanged in multivariate analysis, although loss of observations due to missing data resulted in a non-significant P value. It is tantalizing to hypothesize that the different East-Asian sublineages may differ in their capacity to cause secondary cases.

Does M. tuberculosis genomic diversity explain disease diversity?

The outcome of tuberculosis infection and disease is highly variable. This variation has been attributed primarily to host and environmental factors, but better understanding of the global genomic diversity in the M. tuberculosis complex (MTBC) suggests that bacterial factors could also be involved. Review of nearly 100 published reports shows that MTBC strains differ in their virulence and immunogenicity in experimental models, but whether this phenotypic variation plays a role in human disease remains unclear. Given the complex interactions between the host, the pathogen and the environment, linking MTBC genotypic diversity to experimental and clinical phenotypes requires an integrated systems epidemiology approach embedded in a robust evolutionary framework.

Systems biology approaches to understanding mycobacterial survival mechanisms

The advent of high-throughput platforms for the interrogation of biological systems at the cellular and molecular level have allowed living cells to be observed and understood at a hitherto unprecedented level of detail and have enabled the construction of comprehensive, predictive in silico models. Here, we review the application of such high-throughput, systems-biological techniques to mycobacteria-specifically to the pernicious human pathogen Mycobacterium tuberculosis (MTb) and its ability to survive in human hosts. We discuss the development and application of transcriptomic, proteomic, regulomic, and metabolomic techniques for MTb as well as the development and application of genome-scale in silico models. Thus far, systems-biological approaches have largely focused on in vitro models of MTb growth; reliably extending these approaches to in vivo conditions relevant to infection is a significant challenge for the future that holds the ultimate promise of novel chemotherapeutic interventions.

Possible underlying mechanisms for successful emergence of the Mycobacterium tuberculosis Beijing genotype strains

The wide geographic distribution of one clade of Mycobacterium tuberculosis, the Beijing genotype family, and its genetic homogeneity, suggests that strains belonging to this grouping might have a selective advantage over other M tuberculosis strains. This hypothesis was addressed by reviewing molecular-epidemiological, experimental, and clinical studies. Beijing strains represent about 50% of strains in east Asia and at least 13% of strains worldwide. Their emergence might be linked to escape from BCG vaccination, and to multidrug resistance, which is associated with the Beijing genotype in many areas. Different animal models have shown Beijing strains to be more virulent, and to cause more histopathological changes, higher outgrowth, and increased mortality. At a molecular level, Beijing strains have specific properties in terms of protein and lipid structures and their interaction with the human immune system. Finally, the Beijing genotype has been linked to polymorphisms in immune genes, suggesting the possibility of human-mycobacterial co-evolution. The emergence of the Beijing genotype family might represent an evolutionary response of M tuberculosis to vaccination or antibiotic treatment, with an important negative impact on tuberculosis control. More research is needed to further unravel the mechanisms underlying the emergence of M tuberculosis Beijing genotype strains, and examine the implications for future control strategies.

Strains of the East Asian (W/Beijing) lineage of Mycobacterium tuberculosis are DosS/DosT-DosR two-component regulatory system natural mutants

As part of our ongoing efforts to uncover the phenotypic consequences of genetic variability among clinical Mycobacterium tuberculosis isolates, we previously reported that isolates of the "East Asian" or "W/Beijing" lineage constitutively overexpress the coordinately regulated transcriptional program known as the DosR regulon under standard in vitro conditions. This phenotype distinguishes the W/Beijing lineage from all other M. tuberculosis lineages, which normally induce expression of this regulon only once exposed to low oxygen or nitric oxide, both of which result in inhibition of bacterial respiration and replication. Transcription of the DosR regulon is controlled through a two-component regulatory system comprising the transcription factor DosR and two possible cognate histidine sensor kinases, DosS and DosT. Through sequence analysis of a carefully selected set of isolates representing each of the major M. tuberculosis lineages, we describe herein a naturally occurring frameshift mutation in the gene encoding the DosT sensor kinase for isolates of the most recently evolved W/Beijing sublineages. Intriguingly, the occurrence of the frameshift mutation correlates precisely with the appearance of the constitutive DosR regulon phenotype displayed by the same "modern" W/Beijing strains. However, complementation studies have revealed that the mutation in dosT alone is not directly responsible for the constitutive DosR regulon phenotype. Our data serve to highlight the evolutionary pressure that exists among distinct M. tuberculosis lineages to maintain tight control over DosR regulon expression.

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.

Selection of mutations to detect multidrug-resistant Mycobacterium tuberculosis strains in Shanghai, China

Novel tools are urgently needed for the rapid, reliable detection of multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains of Mycobacterium tuberculosis. To develop such tools, we need information about the frequency and distribution of the mycobacterial mutations and genotypes that are associated with phenotypic drug resistance. In a population-based study, we sequenced specific genes of M. tuberculosis that were associated with resistance to rifampin and isoniazid in 242 phenotypically MDR isolates and 50 phenotypically pan-susceptible isolates from tuberculosis (TB) cases in Shanghai, China. We estimated the sensitivity and specificity of the mutations, using the results of conventional, culture-based phenotypic drug susceptibility testing as the standard. We detected mutations within the 81-bp core region of rpoB in 96.3% of phenotypically MDR isolates. Mutations in two structural genes (katG and inhA) and two regulatory regions (the promoter of mabA-inhA and the intergenic region of oxyR-ahpC) were found in 89.3% of the MDR isolates. In total, 88.0% (213/242 strains) of the phenotypic MDR strains were confirmed by mutations in the sequenced regions. Mutations in embB306 were also considered a marker for MDR and significantly increased the sensitivity of the approach. Based on our findings, an approach that prospectively screens for mutations in 11 sites of the M. tuberculosis genome (rpoB531, rpoB526, rpoB516, rpoB533, and rpoB513, katG315, inhA-15, ahpC-10, ahpC-6, and ahpC-12, and embB306) could detect 86.8% of MDR strains in Shanghai. This study lays the foundation for the development of a rapid, reliable molecular genetic test to detect MDR strains of M. tuberculosis in China.

Pyrosequencing identification of Mycobacterium tuberculosis W-Beijing

Background

The worldwide expanding Mycobacterium tuberculosis W-Beijing family is associated with treatment failure and relapse. Its identification currently relies on spoligotyping and conventional sequencing. We developed pyrosequencing as an alternative method for its identification.

Findings

Pyrosequencing found a G/A substitution in the Rv0927c-pstS3 intergenic spacer and a RD105 deletion, identifying 8/104 M. tuberculosis isolates as W-Beijing isolates. In addition, pyrosequencing found a previously unreported TGC deletion in the Rv0927c gene of W-Beijing isolates. Total concordance was found between the pyrosequencing data and conventional sequencing, as well as reference molecular identification. Multispacer Sequence Typing assigned the W-Beijing isolates to the Asian lineage and the 96 non-W-Beijing isolates to the Euro-American lineage (P < 10^-5). The W-Beijing isolates were all susceptible to streptomycin, rifampin, isoniazid, ethambutol, and pyrazinamide; no resistance-associated mutations were detected in these eight W-Beijing isolates. There were no statistically significant differences in the antibiotic susceptibility of W-Beijing and non-W-Beijing isolates (p = 0.2, X² test). Pyrosequencing correctly identified M. tuberculosis organisms in 26/26 sputum specimens exhibiting acid-fast bacilli. Pyrosequencing results were obtained within four hours, incurring an estimated cost of 1.86 €/test.

Conclusion

Pyrosequencing of the Rv0927c gene and adjacent intergenic spacer is an efficient, low-cost technique for the rapid identification of W-Beijing isolates.

Lupus vulgaris occurring in a locus minoris resistentiae

BACKGROUND

The pathogenesis of lupus vulgaris, a form of cutaneous tuberculosis, is not always clear, especially in patients who do not have coexistent extracutaneous tuberculosis and in patients with single lesions.

OBJECTIVES

To report a case of lupus vulgaris in a locus minoris resistentiae (a site of reduced resistance) and to use a unique set of clinical circumstances and laboratory tests to reconstruct the pathogenesis of the lesion and the response to treatment.

CONCLUSION

Lupus vulgaris can occur in a locus minoris resistentiae; local trauma and possibly other factors, such as increased temperature, topical corticosteroids, and the virulence of the infecting strain, may facilitate the growth of Mycobacterium tuberculosis present at a locus minoris resistentiae as a result of a silent bacillemia.

Genotyping of genetically monomorphic bacteria: DNA sequencing in Mycobacterium tuberculosis highlights the limitations of current methodologies

Because genetically monomorphic bacterial pathogens harbour little DNA sequence diversity, most current genotyping techniques used to study the epidemiology of these organisms are based on mobile or repetitive genetic elements. Molecular markers commonly used in these bacteria include Clustered Regulatory Short Palindromic Repeats (CRISPR) and Variable Number Tandem Repeats (VNTR). These methods are also increasingly being applied to phylogenetic and population genetic studies. Using the Mycobacterium tuberculosis complex (MTBC) as a model, we evaluated the phylogenetic accuracy of CRISPR- and VNTR-based genotyping, which in MTBC are known as spoligotyping and Mycobacterial Interspersed Repetitive Units (MIRU)-VNTR-typing, respectively. We used as a gold standard the complete DNA sequences of 89 coding genes from a global strain collection. Our results showed that phylogenetic trees derived from these multilocus sequence data were highly congruent and statistically robust, irrespective of the phylogenetic methods used. By contrast, corresponding phylogenies inferred from spoligotyping or 15-loci-MIRU-VNTR were incongruent with respect to the sequence-based trees. Although 24-loci-MIRU-VNTR performed better, it was still unable to detect all strain lineages. The DNA sequence data showed virtually no homoplasy, but the opposite was true for spoligotyping and MIRU-VNTR, which was consistent with high rates of convergent evolution and the low statistical support obtained for phylogenetic groupings defined by these markers. Our results also revealed that the discriminatory power of the standard 24 MIRU-VNTR loci varied by strain lineage. Taken together, our findings suggest strain lineages in MTBC should be defined based on phylogenetically robust markers such as single nucleotide polymorphisms or large sequence polymorphisms, and that for epidemiological purposes, MIRU-VNTR loci should be used in a lineage-dependent manner. Our findings have implications for strain typing in other genetically monomorphic bacteria.

Genomic diversity among drug sensitive and multidrug resistant isolates of Mycobacterium tuberculosis with identical DNA fingerprints

Background

Mycobacterium tuberculosis complex (MTBC), the causative agent of tuberculosis (TB), is characterized by low sequence diversity making this bacterium one of the classical examples of a genetically monomorphic pathogen. Because of this limited DNA sequence variation, routine genotyping of clinical MTBC isolates for epidemiological purposes relies on highly discriminatory DNA fingerprinting methods based on mobile and repetitive genetic elements. According to the standard view, isolates exhibiting the same fingerprinting pattern are considered direct progeny of the same bacterial clone, and most likely reflect ongoing transmission or disease relapse within individual patients.

Methodology/Principal Findings

Here we further investigated this assumption and used massively parallel whole-genome sequencing to compare one drug-susceptible (K-1) and one multidrug resistant (MDR) isolate (K-2) of a rapidly spreading M. tuberculosis Beijing genotype clone from a high incidence region (Karakalpakstan, Uzbekistan). Both isolates shared the same IS6110 RFLP pattern and the same allele at 23 out of 24 MIRU-VNTR loci.

We generated 23.9 million (K-1) and 33.0 million (K-2) paired 50 bp purity filtered reads corresponding to a mean coverage of 483.5 fold and 656.1 fold respectively. Compared with the laboratory strain H37Rv both Beijing isolates shared 1,209 SNPs. The two Beijing isolates differed by 130 SNPs and one large deletion. The susceptible isolate had 55 specific SNPs, while the MDR variant had 75 specific SNPs, including the five known resistance-conferring mutations.

Conclusions

Our results suggest that M. tuberculosis isolates exhibiting identical DNA fingerprinting patterns can harbour substantial genomic diversity. Because this heterogeneity is not captured by traditional genotyping of MTBC, some aspects of the transmission dynamics of tuberculosis could be missed or misinterpreted. Furthermore, a valid differentiation between disease relapse and exogenous reinfection might be impossible using standard genotyping tools if the overall diversity of circulating clones is limited. These findings have important implications for clinical trials of new anti-tuberculosis drugs.