Comparison of Mycobacterium tuberculosis genomes reveals frequent deletions in a 20 kb variable region in clinical isolates

Analysis of the twenty-six largest outbreaks of tuberculosis in Aragon using whole-genome sequencing for surveillance purposes

The incidence of tuberculosis in Aragon, Spain, is around ten cases per 100,000 inhabitants. Since 2004, a molecular surveillance protocol has been carried out; therefore, all M. tuberculosis strains are genotyped. Recently, whole-genome sequencing has been implemented for relevant isolates. The aim of this work is to characterise at the molecular level the causative strains of the 26 largest outbreaks of the community (including ten or more cases), genotyped by IS6110-RFLP and causing 26% of tuberculosis cases. To achieve this objective, two or three isolates of each IS6110-cluster belonging to different years were selected for sequencing. We found that strains of lineages L4.8, L4.3 and L4.1.2 were the most frequent. The threshold of 12 SNPs as the maximum distance for confirming the belonging to an outbreak was met for 18 of the 26 IS6110-clusters. Four pairs of isolates with more than 90 SNPs were identified as not belonging to the same strain, and four other pairs were kept in doubt as the number of SNPs was close to 12, between 14 and 35. The study of Regions of Difference revealed that they are lineage conserved. Moreover, we could analyse the IS6110 locations for all genome-sequenced isolates, finding some frequent locations in isolates belonging to the same lineage and certain IS6110 movements between the paired isolates. In the vast majority, these movements were not captured by the IS6110-RFLP pattern. After classifying the genes containing SNP by their functional category, we could confirm that the number of SNPs detected in genes considered as virulence factors and the number of cases the strain produced were not related, suggesting that a particular SNP is more relevant than the number. The characteristics found in the most successful strains in our community could be useful for other researchers in epidemiology, virulence and pathogenesis.

Identification of gene fusion events in Mycobacterium tuberculosis that encode chimeric proteins

Mycobacterium tuberculosis is a facultative intracellular pathogen responsible for causing tuberculosis. The harsh environment in which M. tuberculosis survives requires this pathogen to continuously adapt in order to maintain an evolutionary advantage. However, the apparent absence of horizontal gene transfer in M. tuberculosis imposes restrictions in the ways by which evolution can occur. Large-scale changes in the genome can be introduced through genome reduction, recombination events and structural variation. Here, we identify a functional chimeric protein in the ppe38-71 locus, the absence of which is known to have an impact on protein secretion and virulence. To examine whether this approach was used more often by this pathogen, we further develop software that detects potential gene fusion events from multigene deletions using whole genome sequencing data. With this software we could identify a number of other putative gene fusion events within the genomes of M. tuberculosis isolates. We were able to demonstrate the expression of one of these gene fusions at the protein level using mass spectrometry. Therefore, gene fusions may provide an additional means of evolution for M. tuberculosis in its natural environment whereby novel chimeric proteins and functions can arise.

Investigation of a rapidly spreading tuberculosis outbreak using whole-genome sequencing

This paper describes the application of whole-genome sequencing (WGS) to investigate an outbreak of Mycobacterium tuberculosis occurring in Aragon, Spain, where strains have been submitted to genotyping since 2004. The responsible outbreak strain appeared in our region first in 2014 and it spread to 14 patients in the following three years. WGS found low variability between the isolates with none of the SNPs differences detected more than once, all of which were attributed to a recent transmission. Although two ambiguous bases linked two cases with those who presented the SNP in the same position, the establishment of a definitive transmission route was not possible. The epidemiological data supported the existence of a super-spreader, probably responsible for the majority of the cases involved since there was a two-year delay in diagnoses among cases. This fact would also help explaining the low variability found. The index case was not identified, possibly because it was not diagnosed in Aragon. In addition WGS characterised the strain as a Linage 4.3.3/LAM family and corroborated the susceptibility to anti-tuberculosis drugs observed by the clinical laboratories. This work shows the need to have epidemiological data to support the genomic data in order to clarify the evolution of tuberculosis outbreaks.

Dispersal of Mycobacterium tuberculosis Driven by Historical European Trade in the South Pacific

Mycobacterium tuberculosis (Mtb) is a globally distributed bacterial pathogen whose population structure has largely been shaped by the activities of its obligate human host. Oceania was the last major global region to be reached by Europeans and is the last region for which the dispersal and evolution of Mtb remains largely unexplored. Here, we investigated the evolutionary history of the Euro-American L4.4 sublineage and its dispersal to the South Pacific. Using a phylodynamics approach and a dataset of 236 global Mtb L4.4 genomes we have traced the origins and dispersal of L4.4 strains to New Zealand. These strains are predominantly found in indigenous Māori and Pacific people and we identify a clade of European, likely French, origin that is prevalent in indigenous populations in both New Zealand and Canada. Molecular dating suggests the expansion of European trade networks in the early 19th century drove the dispersal of this clade to the South Pacific. We also identify historical and social factors within the region that have contributed to the local spread and expansion of these strains, including recent Pacific migrations to New Zealand and the rapid urbanization of Māori in the 20th century. Our results offer new insight into the expansion and dispersal of Mtb in the South Pacific and provide a striking example of the role of historical European migrations in the global dispersal of Mtb.

RD5-mediated lack of PE_PGRS and PPE-MPTR export in BCG vaccine strains results in strong reduction of antigenic repertoire but little impact on protection

Tuberculosis is the deadliest infectious disease worldwide. Although the BCG vaccine is widely used, it does not efficiently protect against pulmonary tuberculosis and an improved tuberculosis vaccine is therefore urgently needed. Mycobacterium tuberculosis uses different ESX/Type VII secretion (T7S) systems to transport proteins important for virulence and host immune responses. We recently reported that secretion of T7S substrates belonging to the mycobacteria-specific Pro-Glu (PE) and Pro-Pro-Glu (PPE) proteins of the PGRS (polymorphic GC-rich sequences) and MPTR (major polymorphic tandem repeat) subfamilies required both a functional ESX-5 system and a functional PPE38/71 protein for secretion. Inactivation of ppe38/71 and the resulting loss of PE_PGRS/PPE-MPTR secretion were linked to increased virulence of M. tuberculosis strains. Here, we show that a predicted total of 89 PE_PGRS/PPE-MPTR surface proteins are not exported by certain animal-adapted strains of the M. tuberculosis complex including M. bovis. This Δppe38/71-associated secretion defect therefore also occurs in the M. bovis-derived tuberculosis vaccine BCG and could be partially restored by introduction of the M. tuberculosis ppe38-locus. Epitope mapping of the PPE-MPTR protein PPE10, further allowed us to monitor T-cell responses in splenocytes from BCG/M. tuberculosis immunized mice, confirming the dependence of PPE10-specific immune-induction on ESX-5/PPE38-mediated secretion. Restoration of PE_PGRS/PPE-MPTR secretion in recombinant BCG neither altered global antigenic presentation or activation of innate immune cells, nor protective efficacy in two different mouse vaccination-infection models. This unexpected finding stimulates a reassessment of the immunomodulatory properties of PE_PGRS/PPE-MPTR proteins, some of which are contained in vaccine formulations currently in clinical evaluation.

One of the major findings of the pioneering Mycobacterium tuberculosis H37Rv genome sequencing project was the identification of the highly abundant PE and PPE proteins, named after their N-terminal motifs Pro–Glu (PE) or Pro–Pro–Glu (PPE). Within the 20 years of research since then, many claims were made that PE/PPE proteins, including the two large subgroups encoded by repetitive sequences with very high GC content (PE_PGRS and PPE-MPTR families), are exported to the bacterial surface or beyond, and show broad immunomodulatory impact on host-pathogen interaction. We thus screened strains from different branches of the M. tuberculosis complex, including the attenuated Mycobacterium bovis BCG vaccine strains, for their capacity to export PE_PGRS/PPE-MPTR proteins. Strikingly, we found that BCG strains were unable to export the plethora of PE_PGRS/PPE-MPTR proteins due to the absence of the region of difference RD5, which in M. tuberculosis encodes PPE38, required for PE_PGRS/PPE-MPTR export. Surprisingly, the restoration of PE_PGRS/PPE-MPTR export by genetic complementation in recombinant BCG did not result in immunomodulatory changes or altered protection in mouse models. Our results thus put into perspective the numerous reports on virulence-associated immunomodulatory impact of individual PE_PGRS and PPE-MPTR proteins and open novel questions on their biological function(s).

New insights into the transposition mechanisms of IS6110 and its dynamic distribution between Mycobacterium tuberculosis Complex lineages

The insertion Sequence IS6110, only present in the pathogens of the Mycobacterium tuberculosis Complex (MTBC), has been the gold-standard epidemiological marker for TB for more than 25 years, but biological implications of IS6110 transposition during MTBC adaptation to humans remain elusive. By studying 2,236 clinical isolates typed by IS6110-RFLP and covering the MTBC, we remarked a lineage-specific content of IS6110 being higher in modern globally distributed strains. Once observed the IS6110 distribution in the MTBC, we selected representative isolates and found a correlation between the normalized expression of IS6110 and its abundance in MTBC chromosomes. We also studied the molecular regulation of IS6110 transposition and we found a synergistic action of two post-transcriptional mechanisms: a -1 ribosomal frameshift and a RNA pseudoknot which interferes translation. The construction of a transcriptionally active transposase resulted in 20-fold increase of the transposition frequency. Finally, we examined transposition in M. bovis and M. tuberculosis during laboratory starvation and in a mouse infection model of TB. Our results shown a higher transposition in M. tuberculosis, that preferably happens during TB infection in mice and after one year of laboratory culture, suggesting that IS6110 transposition is dynamically adapted to the host and to adverse growth conditions.

Since the pioneering discovery of transposition by Barbara McClintock in eukaryotes and later in prokaryotes by Robert W. Hedges and Alan E. Jacob, it has become clear the key role of mobile genetics elements in chromosome remodelling, microbial evolution and host adaptation. The insertion sequence IS6110 is widely recognized for its utility in TB diagnosis and epidemiology because it is only present in the M. tuberculosis Complex (MTBC) and its transposition provides an excellent chromosomal polymorphic variability allowing the study of recent TB transmission. This inherent feature of IS6110 leads us to hypothesize that IS6110 plays a crucial role during the TB infectious cycle. However, the biological significance of IS6110 has been hindered by its almost exclusive use as an epidemiological marker. Here, we study the regulatory mechanisms and the distribution of IS6110 in the different MTBC lineages. We discuss the potential biological implications of IS6110, that is much more than an excellent TB epidemiological tool. Since IS6110 could play an important role in the adaptation of MTBC to the host, this study opens new avenues to decipher the biological roles of IS6110 in TB pathogenesis.

Revisiting the role of phospholipases C in virulence and the lifecycle of Mycobacterium tuberculosis

Mycobacterium tuberculosis, the agent of human tuberculosis has developed different virulence mechanisms and virulence-associated tools during its evolution to survive and multiply inside the host. Based on previous reports and by analogy with other bacteria, phospholipases C (PLC) of M. tuberculosis were thought to be among these tools. To get deeper insights into the function of PLCs, we investigated their putative involvement in the intracellular lifestyle of M. tuberculosis, with emphasis on phagosomal rupture and virulence, thereby re-visiting a research theme of longstanding interest. Through the construction and use of an M. tuberculosis H37Rv PLC-null mutant (ΔPLC) and control strains, we found that PLCs of M. tuberculosis were not required for induction of phagosomal rupture and only showed marginal, if any, impact on virulence of M. tuberculosis in the cellular and mouse infection models used in this study. In contrast, we found that PLC-encoding genes were strongly upregulated under phosphate starvation and that PLC-proficient M. tuberculosis strains survived better than ΔPLC mutants under conditions where phosphatidylcholine served as sole phosphate source, opening new perspectives for studies on the role of PLCs in the lifecycle of M. tuberculosis.

Genetic markers, genotyping methods & next generation sequencing in Mycobacterium tuberculosis

Molecular epidemiology (ME) is one of the main areas in tuberculosis research which is widely used to study the transmission epidemics and outbreaks of tubercle bacilli. It exploits the presence of various polymorphisms in the genome of the bacteria that can be widely used as genetic markers. Many DNA typing methods apply these genetic markers to differentiate various strains and to study the evolutionary relationships between them. The three widely used genotyping tools to differentiate Mycobacterium tuberculosis strains are IS6110 restriction fragment length polymorphism (RFLP), spacer oligotyping (Spoligotyping), and mycobacterial interspersed repeat units - variable number of tandem repeats (MIRU-VNTR). A new prospect towards ME was introduced with the development of whole genome sequencing (WGS) and the next generation sequencing (NGS) methods, where the entire genome is sequenced that not only helps in pointing out minute differences between the various sequences but also saves time and the cost. NGS is also found to be useful in identifying single nucleotide polymorphisms (SNPs), comparative genomics and also various aspects about transmission dynamics. These techniques enable the identification of mycobacterial strains and also facilitate the study of their phylogenetic and evolutionary traits.

Marked microevolution of a unique Mycobacterium tuberculosis strain in 17 years of ongoing transmission in a high risk population

The transmission and persistence of Mycobacterium tuberculosis within high risk populations is a threat to tuberculosis (TB) control. In the current study, we used whole genome sequencing (WGS) to decipher the transmission dynamics and microevolution of M. tuberculosis ON-A, an endemic strain responsible for an ongoing outbreak of TB in an urban homeless/under-housed population. Sixty-one M. tuberculosis isolates representing 57 TB cases from 1997 to 2013 were subjected to WGS. Sequencing data was integrated with available epidemiological information and analyzed to determine how the M. tuberculosis ON-A strain has evolved during almost two decades of active transmission. WGS offers higher discriminatory power than traditional genotyping techniques, dividing the M. tuberculosis ON-A strain into 6 sub-clusters, each defined by unique single nucleotide polymorphism profiles. One sub-cluster, designated ON-A^NM (Natural Mutant; 26 isolates from 24 cases) was also defined by a large, 15 kb genomic deletion. WGS analysis reveals the existence of multiple transmission chains within the same population/setting. Our results help validate the utility of WGS as a powerful tool for identifying genomic changes and adaptation of M. tuberculosis.

After the bottleneck: Genome-wide diversification of the Mycobacterium tuberculosis complex by mutation, recombination, and natural selection

Many of the most virulent bacterial pathogens show low genetic diversity and sexual isolation. Accordingly, Mycobacterium tuberculosis, the deadliest human pathogen, is thought to be clonal and evolve by genetic drift. Yet, its genome shows few of the concomitant signs of genome degradation. We analyzed 24 genomes and found an excess of genetic diversity in regions encoding key adaptive functions including the type VII secretion system and the ancient horizontally transferred virulence-related regions. Four different approaches showed evident signs of recombination in M. tuberculosis. Recombination tracts add a high density of polymorphisms, and many are thus predicted to arise from outside the clade. Some of these tracts match Mycobacterium canettii sequences. Recombination introduced an excess of non-synonymous diversity in general and even more in genes expected to be under positive or diversifying selection, e.g., cell wall component genes. Mutations leading to non-synonymous SNPs are effectively purged in MTBC, which shows dominance of purifying selection. MTBC mutation bias toward AT nucleotides is not compensated by biased gene conversion, suggesting the action of natural selection also on synonymous changes. Together, all of these observations point to a strong imprint of recombination and selection in the genome affecting both non-synonymous and synonymous positions. Hence, contrary to some other pathogens and previous proposals concerning M. tuberculosis, this lineage may have come out of its ancestral bottleneck as a very successful pathogen that is rapidly diversifying by the action of mutation, recombination, and natural selection.

Independent large scale duplications in multiple M. tuberculosis lineages overlapping the same genomic region

Mycobacterium tuberculosis, the causative agent of most human tuberculosis, infects one third of the world's population and kills an estimated 1.7 million people a year. With the world-wide emergence of drug resistance, and the finding of more functional genetic diversity than previously expected, there is a renewed interest in understanding the forces driving genome evolution of this important pathogen. Genetic diversity in M. tuberculosis is dominated by single nucleotide polymorphisms and small scale gene deletion, with little or no evidence for large scale genome rearrangements seen in other bacteria. Recently, a single report described a large scale genome duplication that was suggested to be specific to the Beijing lineage. We report here multiple independent large-scale duplications of the same genomic region of M. tuberculosis detected through whole-genome sequencing. The duplications occur in strains belonging to both M. tuberculosis lineage 2 and 4, and are thus not limited to Beijing strains. The duplications occur in both drug-resistant and drug susceptible strains. The duplicated regions also have substantially different boundaries in different strains, indicating different originating duplication events. We further identify a smaller segmental duplication of a different genomic region of a lab strain of H37Rv. The presence of multiple independent duplications of the same genomic region suggests either instability in this region, a selective advantage conferred by the duplication, or both. The identified duplications suggest that large-scale gene duplication may be more common in M. tuberculosis than previously considered.

Systematic survey of clonal complexity in tuberculosis at a populational level and detailed characterization of the isolates involved

Clonally complex infections by Mycobacterium tuberculosis are progressively more accepted. Studies of their dimension in epidemiological scenarios where the infective pressure is not high are scarce. Our study systematically searched for clonally complex infections (mixed infections by more than one strain and simultaneous presence of clonal variants) by applying mycobacterial interspersed repetitive-unit (MIRU)-variable-number tandem-repeat (VNTR) analysis to M. tuberculosis isolates from two population-based samples of respiratory (703 cases) and respiratory-extrapulmonary (R+E) tuberculosis (TB) cases (71 cases) in a context of moderate TB incidence. Clonally complex infections were found in 11 (1.6%) of the respiratory TB cases and in 10 (14.1%) of those with R+E TB. Among the 21 cases with clonally complex TB, 9 were infected by 2 independent strains and the remaining 12 showed the simultaneous presence of 2 to 3 clonal variants. For the 10 R+E TB cases with clonally complex infections, compartmentalization (different compositions of strains/clonal variants in independent infected sites) was found in 9 of them. All the strains/clonal variants were also genotyped by IS6110-based restriction fragment length polymorphism analysis, which split two MIRU-defined clonal variants, although in general, it showed a lower discriminatory power to identify the clonal heterogeneity revealed by MIRU-VNTR analysis. The comparative analysis of IS6110 insertion sites between coinfecting clonal variants showed differences in the genes coding for a cutinase, a PPE family protein, and two conserved hypothetical proteins. Diagnostic delay, existence of previous TB, risk for overexposure, and clustered/orphan status of the involved strains were analyzed to propose possible explanations for the cases with clonally complex infections. Our study characterizes in detail all the clonally complex infections by M. tuberculosis found in a systematic survey and contributes to the characterization that these phenomena can be found to an extent higher than expected, even in an unselected population-based sample lacking high infective pressure.

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.

Tuberculosis vaccine strain Mycobacterium bovis BCG Russia is a natural recA mutant

Background

The current tuberculosis vaccine is a live vaccine derived from Mycobacterium bovis and attenuated by serial in vitro passaging. All vaccine substrains in use stem from one source, strain Bacille Calmette-Guérin. However, they differ in regions of genomic deletions, antigen expression levels, immunogenicity, and protective efficacy.

Results

As a RecA phenotype increases genetic stability and may contribute restricting the ongoing evolution of the various BCG substrains while maintaining their protective efficacy, we aimed to inactivate recA by allelic replacement in BCG vaccine strains representing different phylogenetic lineages (Pasteur, Frappier, Denmark, Russia). Homologous gene replacement was achieved successfully in three out of four strains. However, only illegitimate recombination was observed in BCG substrain Russia. Sequence analyses of recA revealed that a single nucleotide insertion in the 5' part of recA led to a translational frameshift with an early stop codon making BCG Russia a natural recA mutant. At the protein level BCG Russia failed to express RecA.

Conclusion

According to phylogenetic analyses BCG Russia is an ancient vaccine strain most closely related to the parental M. bovis. We hypothesize that recA inactivation in BCG Russia occurred early and is in part responsible for its high degree of genomic stability, resulting in a substrain that has less genetic alterations than other vaccine substrains with respect to M. bovis AF2122/97 wild-type.

Genetic basis of virulence attenuation revealed by comparative genomic analysis of Mycobacterium tuberculosis strain H37Ra versus H37Rv

Tuberculosis, caused by Mycobacterium tuberculosis, remains a leading infectious disease despite the availability of chemotherapy and BCG vaccine. The commonly used avirulent M. tuberculosis strain H37Ra was derived from virulent strain H37 in 1935 but the basis of virulence attenuation has remained obscure despite numerous studies. We determined the complete genomic sequence of H37Ra ATCC25177 and compared that with its virulent counterpart H37Rv and a clinical isolate CDC1551. The H37Ra genome is highly similar to that of H37Rv with respect to gene content and order but is 8,445 bp larger as a result of 53 insertions and 21 deletions in H37Ra relative to H37Rv. Variations in repetitive sequences such as IS6110 and PE/PPE/PE-PGRS family genes are responsible for most of the gross genetic changes. A total of 198 single nucleotide variations (SNVs) that are different between H37Ra and H37Rv were identified, yet 119 of them are identical between H37Ra and CDC1551 and 3 are due to H37Rv strain variation, leaving only 76 H37Ra-specific SNVs that affect only 32 genes. The biological impact of missense mutations in protein coding sequences was analyzed in silico while nucleotide variations in potential promoter regions of several important genes were verified by quantitative RT-PCR. Mutations affecting transcription factors and/or global metabolic regulations related to in vitro survival under aging stress, and mutations affecting cell envelope, primary metabolism, in vivo growth as well as variations in the PE/PPE/PE-PGRS family genes, may underlie the basis of virulence attenuation. These findings have implications not only for improved understanding of pathogenesis of M. tuberculosis but also for development of new vaccines and new therapeutic agents.

Stabilization of the genome of the mismatch repair deficient Mycobacterium tuberculosis by context-dependent codon choice

BACKGROUND

The rate at which a stretch of DNA mutates is determined by the cellular systems for DNA replication and repair, and by the nucleotide sequence of the stretch itself. One sequence feature with a particularly strong influence on the mutation rate are nucleotide repeats. Some microbial pathogens use nucleotide repeats in their genome to stochastically vary phenotypic traits and thereby evade host defense. However, such unstable sequences also come at a cost, as mutations are often deleterious. Here, we analyzed how these opposing forces shaped genome stability in the human pathogen Mycobacterium tuberculosis. M. tuberculosis lacks a mismatch repair system, and this renders nucleotide repeats particularly unstable.

RESULTS

We found that proteins of M. tuberculosis are encoded by using codons in a context-dependent manner that prevents the emergence of nucleotide repeats. This context-dependent codon choice leads to a strong decrease in the estimated frame-shift mutation rate and thus to an increase in genome stability.

CONCLUSION

These results indicate that a context-specific codon choice can partially compensate for the lack of a mismatch repair system, and helps to maintain genome integrity in this pathogen.

Ongoing genome reduction in Mycobacterium ulcerans

M. ulcerans is adapting to a more stable environment.

Elucidation of the transmission, epidemiology, and evolution of Mycobacterium ulcerans, the causative agent of Buruli ulcer, is hampered by the striking lack of genetic diversity of this emerging pathogen. However, by using a prototype plasmid-based microarray that covered 10% of the genome, we found multiple genomic DNA deletions among 30 M. ulcerans clinical isolates of diverse geographic origins. Many of the changes appear to have been mediated by insertion sequence (IS) elements IS2404 and IS2606, which have high copy numbers. Classification of the deleted genes according to their biological functions supports the hypothesis that M. ulcerans has recently evolved from the generalist environmental M. marinum to become a niche-adapted specialist. The substantial genomic diversity, along with a prototype microarray that covered a small portion of the genome, suggests that a genome-wide microarray will make available a genetic fingerprinting method with the high resolution required for microepidemiologic studies.

Discovery of a novel Mycobacterium tuberculosis lineage that is a major cause of tuberculosis in Rio de Janeiro, Brazil

The current study evaluated Mycobacterium tuberculosis isolates from Rio de Janeiro, Brazil, for genomic deletions. One locus in our panel of PCR targets failed to amplify in approximately 30% of strains. A single novel long sequence polymorphism (>26.3 kb) was characterized and designated RD(Rio). Homologous recombination between two similar protein-coding genes is proposed as the mechanism for deleting or modifying 10 genes, including two potentially immunogenic PPE proteins. The flanking regions of the RD(Rio) locus were identical in all strains bearing the deletion. Genetic testing by principal genetic group, spoligotyping, variable-number tandem repeats of mycobacterial interspersed repetitive units (MIRU-VNTR), and IS6110-based restriction fragment length polymorphism analysis cumulatively support the idea that RD(Rio) strains are derived from a common ancestor belonging solely to the Latin American-Mediterranean spoligotype family. The RD(Rio) lineage is therefore the predominant clade causing tuberculosis (TB) in Rio de Janeiro and, as indicated by genotypic clustering in MIRU-VNTR analysis, the most significant source of recent transmission. Limited retrospective reviews of bacteriological and patient records showed a lack of association with multidrug resistance or specific risk factors for TB. However, trends in the data did suggest that RD(Rio) strains may cause a form of TB with a distinct clinical presentation. Overall, the high prevalence of this genotype may be related to enhanced virulence, transmissibility, and/or specific adaptation to a Euro-Latin American host population. The identification of RD(Rio) strains outside of Brazil points to the ongoing intercontinental dissemination of this important genotype. Further studies are needed to determine the differential strain-specific features, pathobiology, and worldwide prevalence of RD(Rio) M. tuberculosis.

Genome analysis shows a common evolutionary origin for the dominant strains of Mycobacterium tuberculosis in a UK South Asian community

We have investigated the Mycobacterium tuberculosis strain types present in the South Asian population of the UK, in which tuberculosis is particularly prevalent. In contrast to the widespread Beijing strains which have the variable number tandem repeats (VNTR) profile 42435, isolates with the VNTR profile 42235, jointly with 02335 or 42234 profiles, appear more frequently in tuberculosis patients of South Asian ethnic origin (SA-strains) in the UK than in any other ethnic group. Using microarray-based comparative genomics to distinguish total or partially deleted genes, we found that three of the common deleted regions in the SA-strains were identical to some deleted genes in the strain CH, which caused an outbreak among South Asian patients in Leicester in 2001 but were different from genomic deletions found in Beijing/W strains. Analysis of some of the deleted regions revealed differences in comparison to the strain CH including the polymorphism in some of the PE/PPE and Esat-6 genes, which may be responsible for the diversity of antigenic variation or differences in the activation of the host immune response. Interrupted genes or the replacement by insertion elements was confirmed in some of the deleted genomic regions. Our results are consistent with the hypothesis that the SA-strains may present common features, implying a common origin for this group of strains.

The role of IS6110 in the evolution of Mycobacterium tuberculosis

Members of the Mycobacterium tuberculosis complex contain the transposable element IS6110 which, due to its high numerical and positional polymorphism, has become a widely used marker in epidemiological studies. Here, we review the evidence that IS6110 is not simply a passive or 'junk' DNA sequence, but that, through its transposable activity, it is able to generate genotypic variation that translates into strain-specific phenotypic variation. We also speculate on the role that this variation has played in the evolution of M. tuberculosis and conclude that the presence of a moderate IS6110 copy number within the genome may provide the pathogen with a selective advantage that has aided its virulence.

Purification and characterization of mycobacterial phospholipase A: an activity associated with mycobacterial cutinase

We describe mycobacterial phospholipase A activity (MPLA) and, using reverse genetics, have associated this activity with putative mycobacterial cutinase. PLAs, which hydrolyze fatty acids on phospholipids, play a significant role in human inflammatory states and disease pathogenesis. In prokaryotes, the recognition of their role in virulence is more recent. Cutinases are serine esterases whose primary substrate is cutin, the waxy exterior layer of plants. Mycobacterium tuberculosis has maintained seven putative cutinases, though it should not encounter cutin; we demonstrate that known cutinases and MPLA cleave phospholipids in a PLA-type manner and also hydrolyze Tween. We analyzed cutinase motifs in mycobacteria and found the motif very prevalent. All mycobacteria tested had MPLA activity. These studies suggest an alternative use for putative cutinases by the M. tuberculosis group that is likely related to MPLA activity and lipid metabolism.

Genetic characterization of Mycobacterium tuberculosis clinical isolates with deletions in the plcA–plcB–plcC locus

SETTING

The basis for Mycobacterium tuberculosis virulence is not completely understood. Analysis of the genomic structure of clinical isolates will give information that can be related to biological activities involved in virulence.

OBJECTIVE

To determine the extension of the deletion in the plcA-plcB-plcC locus of selected M. tuberculosis isolates, as well as other changes in the chromosome.

DESIGN

In the present work we characterized a group of M. tuberculosis isolates devoid of the plcA-plcB-plcC locus by PCR, sequencing and microarrays.

RESULTS

PCR amplification of this region demonstrated a complete lack of plcA and plcB ORF's in all of the isolates. The plcC gene was completely deleted in one of the strains (DR-689) and the other three isolates still conserved part of this ORF. The loss of lateral DNA sequences ranged from 3723 to 7646bp. An IS6110 element was present in all tested strains cases, and some isolates presented the insertion of ORF's coding for proteins homologous to the ESAT-6 and QILSS families. Genomic DNA of all the strains was extracted and analyzed with an in-house microarray system to observe loss of other genes possibly implicated in attenuated virulence. Two of the strains presented novel deletions; the rest of the isolates showed deletions already reported for other M. tuberculosis strains. DR-689, a Beijing type M. tuberculosis strain isolated in Canada, showed an IS6110 RFLP and a genomic deletion pattern similar to a San Francisco family of strains, although completely unrelated epidemiologically.

CONCLUSION

Genomic changes in M. tuberculosis seem to occur in a controlled manner and they are possibly related to changes in its pathogenic properties.

Role of large sequence polymorphisms (LSPs) in generating genomic diversity among clinical isolates of Mycobacterium tuberculosis and the utility of LSPs in phylogenetic analysis

Mycobacterium tuberculosis strains contain different genomic insertions or deletions called large sequence polymorphisms (LSPs). Distinguishing between LSPs that occur one time versus ones that occur repeatedly in a genomic region may provide insights into the biological roles of LSPs and identify useful phylogenetic markers. We analyzed 163 clinical M. tuberculosis isolates for 17 LSPs identified in a genomic comparison of M. tuberculosis strains H37Rv and CDC1551. LSPs were mapped onto a single-nucleotide polymorphism (SNP)-based phylogenetic tree created using nine novel SNP markers that were found to reproduce a 212-SNP-based phylogeny. Four LSPs (group A) mapped to a single SNP tree segment. Two LSPs (group B) and 11 LSPs (group C) were inferred to have arisen independently in the same genomic region either two or more than two times, respectively. None of the group A LSPs but one group B LSP and five group C LSPs were flanked by IS6110 sequences in the references strains. Genes encoding members of the proline-glutamic acid or proline-proline-glutamic acid protein families were present only in group B or C LSPs. SNP- versus LSP-based phylogenies were also compared. We classified each isolate into 58 LSP types by using a separate LSP-based phylogenetic analysis and mapped the LSP types onto the SNP tree. LSPs often assigned isolates to the correct phylogenetic lineage; however, significant mistakes occurred for 6/58 (10%) of the LSP types. In conclusion, most LSPs occur in genomic regions that are prone to repeated insertion/deletion events and were responsible for an unexpectedly high degree of genomic variation in clinical M. tuberculosis. Group B and C LSPs may represent polymorphisms that occur due to selective pressure and affect the phenotype of the organism, while group A LSPs are preferable phylogenetic markers.

Novel genetic polymorphisms that further delineate the phylogeny of the Mycobacterium tuberculosis complex

In a previous report, we described a PCR protocol for the differentiation of the various species of the Mycobacterium tuberculosis complex (MTC) on the basis of genomic deletions (R. C. Huard, L. C. de Oliveira Lazzarini, W. R. Butler, D. van Soolingen, and J. L. Ho, J. Clin. Microbiol. 41:1637-1650, 2003). That report also provided a broad cross-comparison of several previously identified, phylogenetically relevant, long-sequence and single-nucleotide polymorphisms (LSPs and SNPs, respectively). In the present companion report, we expand upon the previous work (i) by continuing the evaluation of known MTC phylogenetic markers in a larger collection of tubercle bacilli (n = 125), (ii) by evaluating additional recently reported MTC species-specific and interspecific polymorphisms, and (iii) by describing the identification and distribution of a number of novel LSPs and SNPs. Notably, new genomic deletions were found in various Mycobacterium tuberculosis strains, new species-specific SNPs were identified for "Mycobacterium canettii," Mycobacterium microti, and Mycobacterium pinnipedii, and, for the first time, intraspecific single-nucleotide DNA differences were discovered for the dassie bacillus, the oryx bacillus, and the two Mycobacterium africanum subtype I variants. Surprisingly, coincident polymorphisms linked one M. africanum subtype I genotype with the dassie bacillus and M. microti with M. pinnipedii, thereby suggesting closer evolutionary ties within each pair of species than had been previously thought. Overall, the presented data add to the genetic definitions of several MTC organisms as well as fine-tune current models for the evolutionary history of the MTC.

Identification of an IS6110 insertion site in plcD, the unique phospholipase C gene of Mycobacterium bovis

The IS6110 repetitive element is present in multiple copies in most Mycobacterium tuberculosis complex bacteria, except for Mycobacterium bovis strains, which usually contain a single copy of IS6110 located on a 1·9 kb PvuII fragment of the direct repeat region. IS6110 transposition can disrupt coding regions and is a major force of genomic variation. In a previous work it was demonstrated that phospholipase C genes are preferential loci for IS6110 transposition in M. tuberculosis clinical strains. Bacterial phospholipase C enzymes participate in pathogenic mechanisms used by different organisms, and have been implicated in intracellular survival, cytolysis and cell-to-cell spread. Four phospholipase C genes (plcA, plcB, plcC and plcD) were detected in the genomes of M. tuberculosis, Mycobacterium africanum, Mycobacterium microti and ‘Mycobacterium canettii’. M. bovis and the vaccine strain M. bovis Bacillus Calmette–Guérin contain only the plcD gene. In the present work, the existence of IS6110 insertions within plcD, the unique phospholipase C gene of M. bovis, has been investigated by PCR, Southern blot hybridization and sequencing analysis. In 18 (7·3 %) of 245 isolates analysed, the plcD gene was interrupted by the insertion of one copy of IS6110, which in all cases was transposed in the same orientation and at the same position, 1 972 894, relative to the genome of M. bovis AF2122/97. These 18 isolates were distributed in 6 different spoligotype patterns and contained 4 to 8 IS6110 copies. In contrast, strains showing an intact plcD gene contained one (87 %), two (9·4 %) or three (2·4 %) IS6110 copies, and only a single isolate (1·2 %) had four IS6110 copies. The implications of plcD gene disruption in M. bovis have not been fully investigated, but no differences in the organ distribution of the disease were detected when animals infected with strains from the same spoligotype patterns bearing plcD : : IS6110 and intact plcD were compared.

The genetic portrait of an outbreak strain

Unique events in the genome of Mycobacterium tuberculosis, including deletions and IS6110 insertions, have been proposed to be responsible for the virulence phenotype of outbreak strains. Based on this premise, we determined ten IS6110 insertion sites in the genome of the M. tuberculosis CH strain, which was responsible for a large outbreak in Leicestershire, England. Together with previous data on genomic deletions, it was found that 16 genes were mutated either by IS6110 insertions or deletions. The likely impact of these genomic events on the phenotype of the CH strain is discussed.

A genome-wide sequence-independent comparative analysis of insertion-deletion polymorphisms in multiple Mycobacterium tuberculosis strains

We applied an enhanced version of subtractive hybridization for comparative analyses of indel differences between genomes of several Mycobacterium tuberculosis strains widespread in Russian regions, and the H37Rv reference strain. A number of differences were detected and partially analyzed, thus demonstrating the practicality of the approach. A majority of the insertions found were shared by all Russian strains, except for strain 1540 that revealed the highest virulence in animal tests. This strain possesses a number of genes absent from other clinical strains. Two of the differential genes were found to encode putative membrane proteins and are presumed to affect mycobacterial interaction with the host cell, thus enhancing virulent properties of the isolate. The method used is of general application, and enables the elaboration of a catalogue of indel polymorphic genomic differences between closely related strains.

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

Beijing/W strains of Mycobacterium tuberculosis are geographically widespread and hypervirulent. To enhance our understanding of their origin and evolution, we sought phylogenetically informative large sequence polymorphisms (LSPs) within the Beijing/W family. Comparative whole-genome hybridization of Beijing/W strains revealed 21 LSPs, 7 of which were previously unreported. We show that some of these LSPs are unique event polymorphisms that can be used to define and subdivide the Beijing/W family. One LSP (RD105) was seen in all Beijing/W strains and thus serves as a useful marker for the identification of this family of strains. Additional LSPs (RD142, RD150, and RD181) further divided this family into four monophyletic subgroups, demonstrating a deeper population structure than previously appreciated. All Beijing/W strains were also observed to have an intact pks15/1 gene that is involved in the biosynthesis of a phenolic glycolipid, a putative virulence factor. A simple PCR assay using these Beijing/W strain-defining deletions will facilitate molecular epidemiological studies and may assist in the identification of the molecular basis of phenotypes associated with this important lineage of M. tuberculosis.

Insertion- and deletion-associated genetic diversity of Mycobacterium tuberculosis phospholipase C-encoding genes among 106 clinical isolates from Turkey

Bacterial phospholipase C has been reported to play a role in the pathogenesis of many bacteria. In order to gain a better understanding of the potential role of Mycobacterium tuberculosis phospholipase C in the pathogenesis of human tuberculosis, we investigated the genetic diversity of the four M. tuberculosis phospholipase C-encoding genes (plcA, plcB, plcC, and plcD) resulting from the IS6110 insertion and associated deletion, among 106 clinical isolates obtained from Turkey, by using PCR, Southern hybridization, and DNA sequencing. Two sequenced M. tuberculosis strains, H37Rv and CDC1551, were used as the references in the comparison. Sixty-six (62.3%) of the 106 isolates had an intact plcD gene, and 40 (37.7%) showed an interruption of the gene. Of the latter 40 isolates, 19 (47.5%) had an IS6110 insertion with no associated deletion in the plcD gene, 2 (5%) had an IS6110 insertion and an associated deletion within the plcD gene, 15 (37.5%) had an IS6110 insertion in the plcD gene that was associated with a partial deletion of the plcD gene and its right forward adjacent region, and 4 (10%) had a complete deletion of the plcD gene. The proportions of the isolates with an interrupted plcA, plcB, or plcC gene were 3.8, 1.9, and 3.8%, respectively. The data indicate that there is a much higher frequency of IS6110 insertion and deletion in the plcD gene than in the plcA, plcB, and plcC genes of M. tuberculosis.

Clinical relevance of Mycobacterium tuberculosis plcD gene mutations

To identify Mycobacterium tuberculosis virulence factors, we integrated comparative genomics and epidemiologic data analysis to investigate the relationship between certain genomic insertions and deletions in the phospholipase-C gene D (plcD) with the clinical presentation of tuberculosis (TB). Four hundred ninety-six well-characterized M. tuberculosis clinical isolates were studied. Approximately 30% (147) of the isolates had an interruption of the plcD gene. Patients infected with the plcD mutant were twice as likely to have extrathoracic disease as those infected by a strain without an interruption (adjusted odds ratio, 2.19; 95% confidence interval, 1.27, 3.76). When we limited the analysis to the 275 isolates with distinct DNA fingerprint patterns, we observed the same association (adjusted odds ratio, 2.74; 95% confidence interval, 1.35, 5.56). Furthermore, the magnitude of the association appeared to differ with the type of extrathoracic TB. Our findings suggest that the plcD gene of M. tuberculosis is potentially involved in the pathogenesis of TB, and the clinical presentation of the disease may be influenced by the genetic variability of the plcD region.

Tuberculosis: from prehistory to Robert Koch, as revealed by ancient DNA

During the past 10 years palaeomicrobiology, a new scientific discipline, has developed. The study of ancient pathogens by direct detection of their DNA has answered several historical questions and shown changes to pathogens over time. However, ancient DNA (aDNA) continues to be controversial and great care is needed to provide valid data. Here we review the most successful application of the technology, which is the study of tuberculosis. This has provided direct support for the current theory of Mycobacterium tuberculosis evolution, and suggests areas of investigation for the interaction of M tuberculosis with its host.

Use of genome level-informed PCR as a new investigational approach for analysis of outbreak-associated Mycobacterium tuberculosis isolates

Mycobacterium tuberculosis strain CH, the index isolate linked to a major tuberculosis outbreak associated with high levels of transmissibility and virulence, was characterized by microarray analysis by use of a PCR product array representative of the genome of M. tuberculosis strain H37Rv. Seven potential genomic deletions were identified in CH, five of which were confirmed by PCR analysis across the predicted deletion points. The panel of five PCRs required to individually interrogate these loci was collectively referred to as the genome level-informed PCR (GLIP) assay. GLIP analysis was performed with CH, 12 other epidemiologically linked isolates, and 43 recent, non-outbreak-associated isolates derived from patients within the local area. All 13 outbreak-linked isolates showed a profile corresponding to the presence of all five deletions. These 13 isolates were also found to share common variable-number tandem repeat and mycobacterial interspersed repetitive unit profiles. None of the 43 non-outbreak-associated isolates exhibited the five-deletion profile. Although three individual deletions were present in upwards of 44% of the non-outbreak-associated isolates, no single-deletion isolates were detected. Interestingly, none of these deletions had been previously recognized, and sequence analysis of the immediate flanking regions in CH failed to identify a likely mechanism of deletion for four of the five loci. The GLIP assay also proved valuable in ongoing surveillance of the outbreak, rapidly identifying a further two outbreak-associated cases months after the initial cluster and, importantly, dismissing a further 12 epidemiologically suspect cases, which allowed the optimum deployment of public health resources.

Lack of mismatch correction facilitates genome evolution in mycobacteria

In silico genome sequence analyses suggested that mycobacteria are devoid of the highly conserved mutLS-based post-replicative mismatch repair system. Here, we present the first biological evidence for the lack of a classical mismatch repair function in mycobacteria. We found that frameshifts, but not general mutation rates are unusually high in Mycobacterium smegmatis. However, despite the absence of mismatch correction, M. smegmatis establishes a strong barrier to recombination between homeologous DNA sequences. We show that 10-12% of DNA sequence heterology restricts initiation of recombination but not extension of heteroduplex DNA intermediates. Together, the lack of mismatch correction and a high stringency of initiation of homologous recombination provide an adequate strategy for mycobacterial genome evolution, which occurs by gene duplication and divergent evolution.

Genomic characterization of an endemic Mycobacterium tuberculosis strain: evolutionary and epidemiologic implications

In a study of 302 Mycobacterium tuberculosis clinical isolates from the low-incidence Canadian-born population of Quebec, we characterized a large endemic strain family by using genomic deletions. The DS6(Quebec) deleted region (11.4 kb) defined a strain family of 143 isolates encompassing two subgroups: one characterized by pyrazinamide (PZA) susceptibility and the other marked by a PZA-monoresistant phenotype. A second deletion (8 bp) in the pncA gene was shared by all 76 isolates with the PZA resistance phenotype, whereas a third DRv0961 deletion (970 bp) defined a further subset of 15 isolates. From their deletion profiles, we derived a most parsimonious evolutionary scenario and compared multiple standard genotyping modalities (using IS6110 restriction fragment length polymorphism [RFLP], spoligotyping, and mycobacterial interspersed repetitive units [MIRU]) across the deletion-based subgroups. The use of a single genotyping modality yielded an unexpectedly high proportion of clustered isolates for a high IS6110 copy strain (27% by IS6110 RFLP, 61% by MIRU, and 77% by spoligotyping). By combining all three modalities, only 14% were genotypically clustered overall, a result more congruent with the epidemiologic profile of reactivation tuberculosis, as suggested by the older age (mean age, 60 years), rural setting, and low proportion of epidemiologic links. These results provide insight into the evolution of genotypes in endemic strains and the potential for false clustering in molecular epidemiologic studies.

IS6110 functions as a mobile, monocyte-activated promoter in Mycobacterium tuberculosis

The mobile insertion sequence, IS6110, is an important marker in tracking of Mycobacterium tuberculosis strains. Here, we demonstrate that IS6110 can upregulate downstream genes through an outward-directed promoter in its 3' end, thus adding to the significance of this element. Promoter activity was orientation dependent and was localized within a 110 bp fragment adjacent to the right terminal inverted repeat. Transcripts from this promoter, named OP6110, begin approximately 85 bp upstream of the 3' end of IS6110. Use of green fluorescent protein (GFP) expression constructs showed that OP6110 was upregulated in M. tuberculosis during growth in human monocytes and in late growth phases in broth. Analysis of natural insertion sites in M. tuberculosis showed that IS6110 upregulated expression of several downstream genes during growth in human monocytes, including Rv2280 in H37Rv and the PE-PGRS gene, Rv1468c, in the clinical strain 210, which is a member of the Beijing family. Transcription between IS6110 and downstream genes was confirmed by reverse transcription polymerase chain reaction. The ability to activate genes during infection suggests that IS6110 has the potential to influence growth characteristics of different strains, and indicates another mechanism by which IS6110 can impact M. tuberculosis evolution.

IS6110-mediated deletion polymorphism in isogenic strains of Mycobacterium tuberculosis

Previous studies have described IS6110-mediated polymorphism as an important driving force in Mycobacterium tuberculosis genome evolution and have provided indirect evidence for IS6110-driven deletion events. This study provides the first description of an IS6110-mediated deletion event in truly isogenic strains. We also provide further support for the hypothesis that the region from Rv1754 to Rv1765 is a hot spot for IS6110 insertion and deletion events.

Genotypic analysis of Mycobacterium tuberculosis in Bangladesh and prevalence of the Beijing strain

Genotypic analysis was performed on 48 Mycobacterium tuberculosis complex strains collected from a hospital in Dhaka city. Deletion analysis showed that the isolates were all M. tuberculosis; 13 of them were found to be of the "ancestral" type, while 35 were of the "modern" type, indicating that both endemic (ancestral type) and epidemic (modern type) strains cause tuberculosis in Bangladesh. Genotyping based on the spoligotype and variable-number tandem repeats (VNTR) of mycobacterial interspersed repetitive units (MIRU) was also done. A total of 34 strains (71%) were grouped by spoligotyping into nine different clusters; the largest comprised 15 isolates of the Beijing genotype, whereas the remaining eight clusters consisted of two to five isolates. MIRU-VNTR typing detected 32 different patterns among 44 tested strains, and the 15 Beijing strains were further discriminated by MIRU-VNTR typing (7 distinct patterns for the 15 isolates). These results indicate that MIRU-VNTR typing, along with spoligotyping and deletion analysis, can be used effectively for molecular epidemiological studies to determine ongoing transmission clusters; to our knowledge, this is the first report about the type of strains prevailing in Bangladesh.

Macro-array and bioinformatic analyses reveal mycobacterial 'core' genes, variation in the ESAT-6 gene family and new phylogenetic markers for the Mycobacterium tuberculosis complex

To better understand the biology and the virulence determinants of the two major mycobacterial human pathogens Mycobacterium tuberculosis and Mycobacterium leprae, their genome sequences have been determined recently. In silico comparisons revealed that among the 1439 genes common to both M. tuberculosis and M. leprae, 219 genes code for proteins that show no similarity with proteins from other organisms. Therefore, the latter 'core' genes could be specific for mycobacteria or even for the intracellular mycobacterial pathogens. To obtain more information as to whether these genes really were mycobacteria-specific, they were included in a focused macro-array, which also contained genes from previously defined regions of difference (RD) known to be absent from Mycobacterium bovis BCG relative to M. tuberculosis. Hybridization of DNA from 40 strains of the M. tuberculosis complex and in silico comparison of these genes with the near-complete genome sequences from Mycobacterium avium, Mycobacterium marinum and Mycobacterium smegmatis were undertaken to answer this question. The results showed that among the 219 conserved genes, very few were not present in all the strains tested. Some of these missing genes code for proteins of the ESAT-6 family, a group of highly immunogenic small proteins whose presence and number is variable among the genomically highly conserved members of the M. tuberculosis complex. Indeed, the results suggest that, with few exceptions, the 'core' genes conserved among M. tuberculosis H37Rv and M. leprae are also highly conserved among other mycobacterial strains, which makes them interesting potential targets for developing new specific anti-mycobacterial drugs. In contrast, the genes from RD regions showed great variability among certain members of the M. tuberculosis complex, and some new specific deletions in Mycobacterium canettii, Mycobacterium microti and seal isolates were identified and further characterized during this study. Together with the distribution of a particular 6 or 7 bp micro-deletion in the gene encoding the polyketide synthase pks15/1, these results confirm and further extend the revised phylogenetic model for the M. tuberculosis complex recently presented.

Analysis of genetic polymorphisms affecting the four phospholipase C (plc) genes in Mycobacterium tuberculosis complex clinical isolates

The Mycobacterium tuberculosis genome contains four highly related genes which present significant similarity to Pseudomonas aeruginosa genes encoding phospholipase C enzymes. Three of these genes, plcA, plcB and plcC, are organized in tandem (locus plcABC). The fourth gene, plcD, is located in a different region. This study investigates variations in plcABC and plcD genes in clinical isolates of M. tuberculosis, Mycobacterium africanum and ‘Mycobacterium canettii’. Genetic polymorphisms were examined by PCR, Southern blot hybridization, sequence analysis and RT-PCR. Seven M. tuberculosis isolates contain insertions of IS6110 elements within plcA, plcC or plcD. In 19 of 25 M. tuberculosis isolates examined, genomic deletions were identified, resulting in loss of parts of genes or complete genes from the plcABC and/or plcD loci. Partial plcD deletion was observed in one M. africanum isolate. In each case, deletions were associated with the presence of a copy of the IS6110 element and in all occurrences IS6110 was transposed in the same orientation. A mechanism of deletion resulting from homologous recombination of two copies of IS6110 was recognized in a group of genetically related M. tuberculosis isolates. Five M. tuberculosis isolates presented major polymorphisms in the plcABC and plcD regions, along with loss of expression competence that affected all four plc genes. Phospholipase C is a well-known bacterial virulence factor. The precise role of phospholipase C in the pathogenicity of M. tuberculosis is unknown, but considering the potential importance that the plc genes may have in the virulence of the tubercle bacillus, the study of isolates cultured from patients with active tuberculosis bearing genetic variations affecting these genes may provide insights into the significance of phospholipase C enzymes for tuberculosis pathogenicity.

Functional and evolutionary genomics of Mycobacterium tuberculosis: insights from genomic deletions in 100 strains

To better understand genome function and evolution in Mycobacterium tuberculosis, the genomes of 100 epidemiologically well characterized clinical isolates were interrogated by DNA microarrays and sequencing. We identified 68 different large-sequence polymorphisms (comprising 186,137 bp, or 4.2% of the genome) that are present in H37Rv, but absent from one or more clinical isolates. A total of 224 genes (5.5%), including genes in all major functional categories, were found to be partially or completely deleted. Deletions are not distributed randomly throughout the genome but instead tend to be aggregated. The distinct deletions in some aggregations appear in closely related isolates, suggesting a genomically disruptive process specific to an individual mycobacterial lineage. Other genomic aggregations include distinct deletions that appear in phylogenetically unrelated isolates, suggesting that a genomic region is vulnerable throughout the species. Although the deletions identified here are evidently inessential to the causation of disease (they are found in active clinical cases), their frequency spectrum suggests that most are weakly deleterious to the pathogen. For some deletions, short-term evolutionary pressure due to the host immune system or antibiotics may favor the elimination of genes, whereas longer-term physiological requirements maintain the genes in the population.

The majority of inducible DNA repair genes in Mycobacterium tuberculosis are induced independently of RecA

In many species of bacteria most inducible DNA repair genes are regulated by LexA homologues and are dependent on RecA for induction. We have shown previously by analysing the induction of recA that two mechanisms for the induction of gene expression following DNA damage exist in Mycobacterium tuberculosis. Whereas one of these depends on RecA and LexA in the classical way, the other mechanism is independent of both of these proteins and induction occurs in the absence of RecA. Here we investigate the generality of each of these mechanisms by analysing the global response to DNA damage in both wild-type M. tuberculosis and a recA deletion strain of M. tuberculosis using microarrays. This revealed that the majority of the genes that were induced remained inducible in the recA mutant stain. Of particular note most of the inducible genes with known or predicted functions in DNA repair did not depend on recA for induction. Amongst these are genes involved in nucleotide excision repair, base excision repair, damage reversal and recombination. Thus, it appears that this novel mechanism of gene regulation is important for DNA repair in M. tuberculosis.

Genomic interrogation of the dassie bacillus reveals it as a unique RD1 mutant within the Mycobacterium tuberculosis complex

Despite their remarkable genetic homology, members of the Mycobacterium tuberculosis complex express very different phenotypes, most notably in their spectra of clinical presentation. For example, M. tuberculosis is regarded as pathogenic to humans, whereas members having deleted RD1, such as Mycobacterium microti and Mycobacterium bovis BCG, are not. The dassie bacillus, an infrequent variant of the M. tuberculosis complex characterized as being most similar to M. microti, is the causative agent of tuberculosis (TB) in the dassie (Procavia capensis). Intriguingly, the dassie bacillus is not pathogenic to rabbits or guinea pigs and has never been documented to infect humans. Although it was identified more than a half-century ago, the reasons behind its attenuation are unknown. Because large sequence polymorphisms have presented themselves as the most obvious genomic distinction among members of the M. tuberculosis complex, the DNA content of the dassie bacillus was interrogated by Affymetrix GeneChip to identify regions that are absent from it but present in M. tuberculosis H37Rv. Comparison has led to the identification of nine regions of difference (RD), five of which are shared with M. microti (RDs 3, 7, 8, 9, and 10). Although the dassie bacillus does not share the other documented deletions in M. microti (RD1(mic), RD5(mic), MID1, MID2, and MID3), it has endured unique deletions in the regions of RD1, RD5, N-RD25, and Rv3081-Rv3082c (virS). RD1(das), affecting only Rv3874-Rv3877, is the smallest natural deletion of the RD1 region uncovered and points to genes within this region that are likely implicated in virulence. Newfound deletions from the dassie bacillus are discussed in relation to their evolutionary and biological significance.

Patients with active tuberculosis often have different strains in the same sputum specimen

It is generally accepted that tuberculosis results from a single infection with a single Mycobacterium tuberculosis strain. Such infections are thought to confer protective immunity against exogenous reinfection. In this study, a novel polymerase chain reaction method was developed to specifically identify M. tuberculosis strains belonging to the Beijing and non-Beijing evolutionary lineages in sputum specimens collected from tuberculosis patients resident in an epidemiologic field site in Cape Town, South Africa. The sensitivity and specificity of the polymerase chain reaction-based strain classification method were 100% (95% confidence interval, 85-100%) when compared with DNA fingerprinting and spacer oligotyping (spoligotyping). Application of this method showed that 19% of all patients were simultaneously infected with Beijing and non-Beijing strains, and 57% of patients infected with a Beijing strain were also infected with a non-Beijing strain. Multiple infections were more frequent in retreatment cases (23%) as compared with new cases (17%), but were not associated with sex, age, or smear grading. These results suggest that multiple infections are frequent, implying high reinfection rates and the absence of efficient protective immunity conferred by the initial infection. This finding could influence our understanding of the epidemiology of disease in high-incidence regions and our understanding for vaccine development.

Expression of SA5K, a secretion antigen of Mycobacterium tuberculosis, inside human macrophages and in sputum from tuberculosis patients

An 8.3 kDa protein (SA5K), secreted by Mycobacterium tuberculosis/Mycobacterium bovis bacillus Calmette-Guérin (BCG) in culture filtrate, has been previously described in our laboratory. In the present study, analysis of the distribution of SA5K gene (Rv1174c) among M. tuberculosis strains, isolated from a wide variety of clinical specimens, revealed that the gene is present in all clinical isolates analyzed (29/29). SA5K expression inside human macrophages infected with BCG was demonstrated by reverse transcription-polymerase chain reaction (RT-PCR) on RNA extracted from bacterial cells following 24 and 48 h of infection. In addition, in order to evaluate whether SA5K gene was also expressed at the site of infection in the lung, a nested RT-PCR assay was developed to detect specific mRNA in sputum samples collected from smear positive tuberculosis patients. SA5K mRNA was detected in all the samples containing high numbers of tubercle bacilli demonstrating that the corresponding gene is expressed during the course of clinical infection.

Different strains of Mycobacterium tuberculosis cause various spectrums of disease in the rabbit model of tuberculosis

The rabbit model of tuberculosis has been used historically to differentiate between Mycobacterium tuberculosis and Mycobacterium bovis based on their relative virulence in this animal host. M. tuberculosis infection in market rabbits is cleared over time, whereas infection with M. bovis results in chronic, progressive, cavitary disease leading to death. Because of the innate resistance of commercial rabbits to M. tuberculosis, 320 to 1,890 log-phase, actively growing inhaled bacilli were required to form one grossly visible pulmonary tubercle at 5 weeks. The range of inhaled doses required to make one tubercle allows us to determine the relative pathogenicities of different strains. Fewer inhaled organisms of the M. tuberculosis Erdman strain were required than of M. tuberculosis H37Rv to produce a visible lesion at 5 weeks. Furthermore, with the Erdman strain, only 7 of 15 rabbits had healed lesions at 16 to 18 weeks; among the other animals, two had chronic, progressive cavitary disease, a phenotype usually seen only with M. bovis infection. Genotypic investigation of the Erdman strain with an H37Rv-based microarray identified gene differences in the RD6 region. Southern blot and PCR structural genetic analysis showed significant differences between M. tuberculosis strains in this region. Correlation of the relative pathogenicity, including disease severity, in the rabbit model with the strain genotype may help identify stage-specific M. tuberculosis genes important in human disease.

IS6110-mediated deletion polymorphism in the direct repeat region of clinical isolates of Mycobacterium tuberculosis

This study investigates the phenomenon of IS6110-mediated deletion polymorphism in the direct repeat (DR) region of the genome of Mycobacterium tuberculosis. Clinical isolates and their putative predecessors were compared using a combination of DR region restriction fragment length polymorphism, IS6110 DNA fingerprinting, spoligotyping, and DNA sequencing, which allowed the mapping of chromosome structure and deletion junctions. The data suggest that adjacently situated IS6110 elements mediate genome deletion. However, in contrast to previous reports, deletions appear to be mediated by inversely oriented IS6110 elements. This suggests that these events may occur via mechanisms other than RecA-mediated homologous recombination. The results underscore the important role of IS6110-associated deletion hypervariability in driving M. tuberculosis genome evolution.

Phylogenetic reconstruction within Mycobacterium tuberculosis Beijing genotype in northwestern Russia

A selection of genetic markers was used to study the evolution of Mycobacterium tuberculosis Beijing family strains in northwestern Russia. A total of 221 of 434 epidemiologically unlinked isolates studied in 1996-2001 belonged to the Beijing family as determined by standard spoligotyping (signals 35-43). Ninety-six percent of these Beijing isolates ("typical") were closely related in IS6110-RFLP (D > 0.85) while 9 remaining isolates (2 different profiles, "atypical") were more distant from the rest (D = 0.6-0.7). Further analysis was performed on a selection of 12 typical and both atypical Beijing strains with different IS6110-RFLP profiles (2 isolates each). All 28 Beijing isolates studied had the KatG 463Leu allele, an intact mtp40 fragment of the mpcA gene, and an identical structure of the DR locus (15 DVRs) with an upstream IS6110 copy in opposite orientation. The IS6110-RFLP based neighbor-joining (distance) and quartet-puzzling (maximum-likelihood) trees showed that the branch lengths were considerably longer for atypical Beijing strains. Typical Beijing strains had the 1.02 kb Rv3135 PPE-family gene and two IS1547 copies (iplA and iplB) one of them (iplB) disrupted by IS6110 insertion. Atypical Beijing strains had the 1.97 kb Rv3135 gene and a single intact IS1547/iplA copy. We suggest that the M. tuberculosis Beijing family strains currently circulating in the northwest of Russia are relatively ancient and thus appear to be endemic in this region since evolutionarily distant time. The prevalent typical Beijing strains (96%) are likely to be of monophyletic origin and their ongoing dissemination has started recently: these strains differ in rapidly evolving IS6110-RFLP but have identical structure of other polymorphic genome regions studied. The atypical Beijing strains (4%) are evolutionary older; they probably had a common (unknown) ancestor with typical Beijing strains.