Natural transposon mutagenesis of clinical isolates of Mycobacterium tuberculosis: how many genes does a pathogen need?

Functional genetic variation in pe/ppe genes contributes to diversity in Mycobacterium tuberculosis lineages and potential interactions with the human host

Introduction

Around 10% of the coding potential of Mycobacterium tuberculosisis constituted by two poorly understood gene families, the pe and ppe loci, thought to be involved in host-pathogen interactions. Their repetitive nature and high GC content have hindered sequence analysis, leading to exclusion from whole-genome studies. Understanding the genetic diversity of pe/ppe families is essential to facilitate their potential translation into tools for tuberculosis prevention and treatment.

Methods

To investigate the genetic diversity of the 169 pe/ppe genes, we performed a sequence analysis across 73 long-read assemblies representing seven different lineages of M. tuberculosis and M. bovis BCG. Individual pe/ppe gene alignments were extracted and diversity and conservation across the different lineages studied.

Results

The pe/ppe genes were classified into three groups based on the level of protein sequence conservation relative to H37Rv, finding that >50% were conserved, with indels in pe_pgrs and ppe_mptr sub-families being major drivers of structural variation. Gene rearrangements, such as duplications and gene fusions, were observed between pe and pe_pgrs genes. Inter-lineage diversity revealed lineage-specific SNPs and indels.

Discussion

The high level of pe/ppe genes conservation, together with the lineage-specific findings, suggest their phylogenetic informativeness. However, structural variants and gene rearrangements differing from the reference were also identified, with potential implications for pathogenicity. Overall, improving our knowledge of these complex gene families may have insights into pathogenicity and inform the development of much-needed tools for tuberculosis control.

The MtZ Strain: Molecular Characteristics and Outbreak Investigation of the Most Successful Mycobacterium tuberculosis Strain in Aragon Using Whole-Genome Sequencing

Since 2004, a tuberculosis surveillance protocol has been carried out in Aragon, thereby managing to detect all tuberculosis outbreaks that take place in the community. The largest outbreak was caused by a strain named Mycobacterium tuberculosis Zaragoza (MtZ), causing 242 cases as of 2020. The main objective of this work was to analyze this outbreak and the molecular characteristics of this successful strain that could be related to its greater transmission. To do this, we first applied whole-genome sequencing to 57 of the isolates. This revealed two principal transmission clusters and six subclusters arising from them. The MtZ strain belongs to L4.8 and had eight specific single nucleotide polymorphisms (SNPs) in genes considered to be virulence factors [ptpA, mc3D, mc3F, VapB41, pks15 (two SNPs), virS, and VapC50]. Second, a transcriptomic study was carried out to better understand the multiple IS6110 copies present in its genome. This allowed us to observe three effects of IS6110: the disruption of the gene in which the IS6110 is inserted (desA3), the overexpression of a gene (ppe38), and the absence of transcription of genes (cut1:Rv1765c) due to the recombination of two IS6110 copies. Finally, because of the disruption of ppe38 and ppe71 genes by an IS6110, a study of PE_PGRS secretion was carried out, showing that MtZ secretes these factors in higher amounts than the reference strain, thereby differing from the hypervirulent phenotype described for the Beijing strains. In conclusion, MtZ consists of several SNPs in genes related to virulence, pathogenesis, and survival, as well as other genomic polymorphisms, which may be implicated in its success among our population.

Defining Discriminatory Antibody Fingerprints in Active and Latent Tuberculosis

Tuberculosis (TB) is among the leading causes of death worldwide from a single infectious agent, second only to COVID-19 in 2020. TB is caused by infection with Mycobacterium tuberculosis (Mtb), that results either in a latent or active form of disease, the latter associated with Mtb spread. In the absence of an effective vaccine, epidemiologic modeling suggests that aggressive treatment of individuals with active TB (ATB) may curb spread. Yet, clinical discrimination between latent (LTB) and ATB remains a challenge. While antibodies are widely used to diagnose many infections, the utility of antibody-based tests to diagnose ATB has only regained significant traction recently. Specifically, recent interest in the humoral immune response to TB has pointed to potential differences in both targeted antigens and antibody features that can discriminate latent and active TB. Here we aimed to integrate these observations and broadly profile the humoral immune response across individuals with LTB or ATB, with and without HIV co-infection, to define the most discriminatory humoral properties and diagnose TB disease more easily. Using 209 Mtb antigens, striking differences in antigen-recognition were observed across latently and actively infected individuals that was modulated by HIV serostatus. However, ATB and LTB could be discriminated, irrespective of HIV-status, based on a combination of both antibody levels and Fc receptor-binding characteristics targeting both well characterized (like lipoarabinomannan, 38 kDa or antigen 85) but also novel Mtb antigens (including Rv1792, Rv1528, Rv2435C or Rv1508). These data reveal new Mtb-specific immunologic markers that can improve the classification of ATB versus LTB.

In-depth Analysis of IS6110 Genomic Variability in the Mycobacterium tuberculosis Complex

The insertion sequence (IS) 6110 is a repetitive mobile element specific for the Mycobacterium tuberculosis complex (MTBC) used for years to diagnose and genotype this pathogen. It contains the overlapping reading frames orfA and orfB that encode a transposase. Its genetic variability is difficult to study because multiple copies are present in the genome. IS6110 is randomly located, nevertheless some preferential locations have been reported, which could be related to the behaviour of the strains. The aim of this work was to determine the intra- and inter-strain genetic conservation of this element in the MTBC. For this purpose, we analysed 158 sequences of IS6110 copies from 55 strains. Eighty-four copies were from 17 strains for which we knew all the locations in their genome. In addition, we studied 74 IS6110 copies in 38 different MTBC strains in which the location was characteristic of different families including Haarlem, LAM, S, and L6 strains. We observed mutation in 13.3% of the copies studied and we found 10 IS6110 variants in 21 copies belonging to 16 strains. The high copy number strains showed 6.2% of their IS6110 copies mutated, in contrast with the 31.1% in the low-copy-number strains. The apparently more ancient copy localised in the DR region was that with more variant copies, probably because this was the most studied location. Notably, all Haarlem and X family strains studied have an IS6110 in Rv0403c, suggesting a common origin for both families. Nevertheless, we detected a variant specific for the X family that would have occurred in this location after the phylogenetic separation. This variant does not prevent transposition although it may occur at a lower frequency, as X strains remain with low copy number (LCN) of IS6110.

Label-Free Comparative Proteomics of Differentially Expressed Mycobacterium tuberculosis Protein in Rifampicin-Related Drug-Resistant Strains

Rifampicin (RIF) is one of the most important first-line anti-tuberculosis (TB) drugs, and more than 90% of RIF-resistant (RR) Mycobacterium tuberculosis clinical isolates belong to multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB. In order to identify specific candidate target proteins as diagnostic markers or drug targets, differential protein expression between drug-sensitive (DS) and drug-resistant (DR) strains remains to be investigated. In the present study, a label-free, quantitative proteomics technique was performed to compare the proteome of DS, RR, MDR, and XDR clinical strains. We found iniC, Rv2141c, folB, and Rv2561 were up-regulated in both RR and MDR strains, while fadE9, espB, espL, esxK, and Rv3175 were down-regulated in the three DR strains when compared to the DS strain. In addition, lprF, mce2R, mce2B, and Rv2627c were specifically expressed in the three DR strains, and 41 proteins were not detected in the DS strain. Functional category showed that these differentially expressed proteins were mainly involved in the cell wall and cell processes. When compared to the RR strain, Rv2272, smtB, lpqB, icd1, and folK were up-regulated, while esxK, PPE19, Rv1534, rpmI, ureA, tpx, mpt64, frr, Rv3678c, esxB, esxA, and espL were down-regulated in both MDR and XDR strains. Additionally, nrp, PPE3, mntH, Rv1188, Rv1473, nadB, PPE36, and sseA were specifically expressed in both MDR and XDR strains, whereas 292 proteins were not identified when compared to the RR strain. When compared between MDR and XDR strains, 52 proteins were up-regulated, while 45 proteins were down-regulated in the XDR strain. 316 proteins were especially expressed in the XDR strain, while 92 proteins were especially detected in the MDR strain. Protein interaction networks further revealed the mechanism of their involvement in virulence and drug resistance. Therefore, these differentially expressed proteins are of great significance for exploring effective control strategies of DR-TB.

Genetic Diversity in Mycobacterium tuberculosis Clinical Isolates and Resulting Outcomes of Tuberculosis Infection and Disease

Tuberculosis claims more human lives than any other bacterial infectious disease and represents a clear and present danger to global health as new tools for vaccination, treatment, and interruption of transmission have been slow to emerge. Additionally, tuberculosis presents with notable clinical heterogeneity, which complicates diagnosis, treatment, and the establishment of nonrelapsing cure. How this heterogeneity is driven by the diversity ofclinical isolates of the causative agent, Mycobacterium tuberculosis, has recently garnered attention. Herein, we review advances in the understanding of how naturally occurring variation in clinical isolates affects transmissibility, pathogenesis, immune modulation, and drug resistance. We also summarize how specific changes in transcriptional responses can modulate infection or disease outcome, together with strain-specific effects on gene essentiality. Further understanding of how this diversity of M. tuberculosis isolates affects disease and treatment outcomes will enable the development of more effective therapeutic options and vaccines for this dreaded disease.

Isolation and comparative genomics of Mycobacterium tuberculosis isolates from cattle and their attendants in South India

The major human pathogen Mycobacterium tuberculosis is rarely reported to cause disease in other animals. Cases in livestock are thought to occur through contact with infected handlers, but previous studies evaluating putative livestock-human transmission used typing techniques with limited resolution. Here, we undertook cross-sectional surveillance for tuberculosis in 271 livestock handlers and 167 cattle on three farms in Chennai, India and defined the relatedness of cultured isolates using whole genome sequencing. Humans and livestock were screened for active mycobacterial infection, and opportunistic post-mortem examination was performed on comparative intradermal test-positive cattle that died. Four cattle and 6 handlers on two farms were culture-positive for M. tuberculosis; M. bovis was not isolated. All 10 isolates (one from each case) belonged to Lineage 1. Pairwise genome comparisons of single nucleotide polymorphism (SNP) differences ranged from 1 to 600 SNPs, but 3 isolate pairs were less than 5 SNPs different. Two pairs were from handlers and the third pair were from two cattle on the same farm. The minimum pairwise SNP difference between a cattle and human isolate was >250 SNPs. Our study confirms the presence of M. tuberculosis infection in cattle in India, sequencing of which characterised relatedness between human and cattle-derived isolates.

Integrated Multi-Omic Analysis of Mycobacterium tuberculosis H37Ra Redefines Virulence Attributes

H37Ra is a virulence attenuated strain of Mycobacterium tuberculosis widely employed as a model to investigate virulence mechanisms. Comparative high-throughput studies have earlier correlated its avirulence to the presence of specific mutations or absence of certain proteins. However, a recent sequencing study of H37Ra, has disproved several genomic differences earlier reported to be associated with virulence. This warrants further investigations on the H37Ra proteome as well. In this study, we carried out an integrated analysis of the genome, transcriptome, and proteome of H37Ra. In addition to confirming single nucleotide variations (SNVs) and insertion-deletions that were reported earlier, our study provides novel insights into the mutation spectrum in the promoter regions of 7 genes. We also provide transcriptional and proteomic evidence for 3,900 genes representing ~80% of the total predicted gene count including 408 proteins that have not been identified previously. We identified 9 genes whose coding potential was hitherto reported to be absent in H37Ra. These include 2 putative virulence factors belonging to ESAT-6 like family of proteins. Furthermore, proteogenomic analysis enabled us to identify 63 novel proteins coding genes and correct 25 existing gene models in H37Ra genome. A majority of these were found to be conserved in the virulent strain H37Rv as well as in other mycobacterial species suggesting that the differences in the virulent and avirulent strains of M. tuberculosis are not entirely dependent on the expression of certain proteins or their absence but may possibly be ascertained to functional changes.

Analysis of IS6110 insertion sites provide a glimpse into genome evolution of Mycobacterium tuberculosis

Insertion sequence (IS) 6110 is found at multiple sites in the Mycobacterium tuberculosis genome and displays a high degree of polymorphism with respect to copy number and insertion sites. Therefore, IS6110 is considered to be a useful molecular marker for diagnosis and strain typing of M. tuberculosis. Generally IS6110 elements are identified using experimental methods, useful for analysis of a limited number of isolates. Since short read genome sequences generated using next-generation sequencing (NGS) platforms are available for a large number of isolates, a computational pipeline for identification of IS6110 elements from these datasets was developed. This study shows results from analysis of NGS data of 1377 M. tuberculosis isolates. These isolates represent all seven major global lineages of M. tuberculosis. Lineage specific copy number patterns and preferential insertion regions were observed. Intra-lineage differences were further analyzed for identifying spoligotype specific variations. Copy number distribution and preferential locations of IS6110 in different lineages imply independent evolution of IS6110, governed mainly through ancestral insertion, fitness (gene truncation, promoter activity) and recombinational loss of some copies. A phylogenetic tree based on IS6110 insertion data of different isolates was constructed in order to understand genome level variations of different markers across different lineages.

Evolutionary Trends of the Transposase-Encoding Open Reading Frames A and B (orfA and orfB) of the Mycobacterial IS6110 Insertion Sequence

Background

The IS6110 insertion sequence, a member of the IS3 family of insertion sequences, was found to be specific to the Mycobacterium tuberculosis complex (MTBC). Although IS6110 has been extensively characterized as a transposable genetic marker, the evolutionary history of its own transposase-encoding sequence has not, to the best of our knowledge, been investigated.

Methodology/Principal Findings

Here we explored the evolution of the IS6110 sequence by analysing the genetic variability and the selective forces acting on its transposase-encoding open reading frames (ORFs) A and B (orfA and orfB). For this purpose, we used a strain collection consisting of smooth tubercle bacilli (STB), an early branching lineage of the MTBC, and present-day M. tuberculosis strains representing the full breadth of genetic diversity in Tunisia. In each ORF, we found a major haplotype that dominated over a flat distribution of rare descendent haplotypes, consisting mainly of single- and double-nucleotide variant singletons. The predominant haplotypes consisted of both ancestral and present-day strains, suggesting that IS6110 acquisition predated the emergence of the MTBC. There was no evidence of recombination and both ORFs were subjected to strict purifying selection, as demonstrated by their dN/dS ratios (0.29 and 0.51, respectively), as well as their significantly negative Tajima’s D statistics. Strikingly, the purifying selection acting on orfA proved much more stringent, suggesting its critical role in regulating the transpositional process. Maximum likelihood analyses further excluded any possibility of positive selection acting on single amino acid residues.

Conclusions/Significance

Taken together our data fit with an evolutionary scenario according to which the observed variability pattern of the IS6110 transposase-encoding ORFs is generated mainly through random point mutations that accrued on a functionally optimal IS6110 copy, whose acquisition predated the emergence of the MTBC complex. Background selection acting against deleterious mutations led to an excess of low-frequency variants.

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.

Current methods in the molecular typing of Mycobacterium tuberculosis and other mycobacteria

In the epidemiology of tuberculosis (TB) and nontuberculous mycobacterial (NTM) diseases, as in all infectious diseases, the key issue is to define the source of infection and to disclose its routes of transmission and dissemination in the environment. For this to be accomplished, the ability of discerning and tracking individual Mycobacterium strains is of critical importance. Molecular typing methods have greatly improved our understanding of the biology of mycobacteria and provide powerful tools to combat the diseases caused by these pathogens. The utility of various typing methods depends on the Mycobacterium species under investigation as well as on the research question. For tuberculosis, different methods have different roles in phylogenetic analyses and person-to-person transmission studies. In NTM diseases, most investigations involve the search for environmental sources or phylogenetic relationships. Here, too, the type of setting determines which methodology is most suitable. Within this review, we summarize currently available molecular methods for strain typing of M. tuberculosis and some NTM species, most commonly associated with human disease. For the various methods, technical practicalities as well as discriminatory power and accomplishments are reviewed.

Global study of IS6110 in a successful Mycobacterium tuberculosis strain: clues for deciphering its behavior and for its rapid detection

The Mycobacterium tuberculosis insertion sequence IS6110, besides being a very useful tool in molecular epidemiology, seems to have an impact on the biology of bacilli. In the present work, we mapped the 12 points of insertion of IS6110 in the genome of a successful strain named M. tuberculosis Zaragoza (which has been referred to as the MTZ strain). This strain, belonging to principal genetic group 3, caused a large unsuspected tuberculosis outbreak involving 85 patients in Zaragoza, Spain, in 2001 to 2004. The mapping of the points of insertion of IS6110 in the genome of the Zaragoza strain offers clues for a better understanding of the adaptability and virulence of M. tuberculosis. Surprisingly, the presence of one copy of IS6110 was found in Rv2286c, as was recently described for a successful Beijing sublineage. As a result of this analysis, a rapid method for detecting this particular M. tuberculosis strain has been designed.

IS-seq: a novel high throughput survey of in vivo IS6110 transposition in multiple Mycobacterium tuberculosis genomes

Background

The insertion element IS6110 is one of the main sources of genomic variability in Mycobacterium tuberculosis, the etiological agent of human tuberculosis. Although IS 6110 has been used extensively as an epidemiological marker, the identification of the precise chromosomal insertion sites has been limited by technical challenges. Here, we present IS-seq, a novel method that combines high-throughput sequencing using Illumina technology with efficient combinatorial sample multiplexing to simultaneously probe 519 clinical isolates, identifying almost all the flanking regions of the element in a single experiment.

Results

We identified a total of 6,976 IS6110 flanking regions on the different isolates. When validated using reference strains, the method had 100% specificity and 98% positive predictive value. The insertions mapped to both coding and non-coding regions, and in some cases interrupted genes thought to be essential for virulence or in vitro growth. Strains were classified into families using insertion sites, and high agreement with previous studies was observed.

Conclusions

This high-throughput IS-seq method, which can also be used to map insertions in other organisms, extends previous surveys of in vivo interrupted loci and provides a baseline for probing the consequences of disruptions in M. tuberculosis strains.

Comparative analysis of Mycobacterium tuberculosis pe and ppe genes reveals high sequence variation and an apparent absence of selective constraints

Mycobacterium tuberculosis complex (MTBC) genomes contain 2 large gene families termed pe and ppe. The function of pe/ppe proteins remains enigmatic but studies suggest that they are secreted or cell surface associated and are involved in bacterial virulence. Previous studies have also shown that some pe/ppe genes are polymorphic, a finding that suggests involvement in antigenic variation. Using comparative sequence analysis of 18 publicly available MTBC whole genome sequences, we have performed alignments of 33 pe (excluding pe_pgrs) and 66 ppe genes in order to detect the frequency and nature of genetic variation. This work has been supplemented by whole gene sequencing of 14 pe/ppe (including 5 pe_pgrs) genes in a cohort of 40 diverse and well defined clinical isolates covering all the main lineages of the M. tuberculosis phylogenetic tree. We show that nsSNP's in pe (excluding pgrs) and ppe genes are 3.0 and 3.3 times higher than in non-pe/ppe genes respectively and that numerous other mutation types are also present at a high frequency. It has previously been shown that non-pe/ppe M. tuberculosis genes display a remarkably low level of purifying selection. Here, we also show that compared to these genes those of the pe/ppe families show a further reduction of selection pressure that suggests neutral evolution. This is inconsistent with the positive selection pressure of "classical" antigenic variation. Finally, by analyzing such a large number of genes we were able to detect large differences in mutation type and frequency between both individual genes and gene sub-families. The high variation rates and absence of selective constraints provides valuable insights into potential pe/ppe function. Since pe/ppe proteins are highly antigenic and have been studied as potential vaccine components these results should also prove informative for aspects of M. tuberculosis vaccine design.

Characterization of microevolution events in Mycobacterium tuberculosis strains involved in recent transmission clusters

Under certain circumstances, it is possible to identify clonal variants of Mycobacterium tuberculosis infecting a single patient, probably as a result of subtle genetic rearrangements in part of the bacillary population. We systematically searched for these microevolution events in a different context, namely, recent transmission chains. We studied the clustered cases identified using a population-based universal molecular epidemiology strategy over a 5-year period. Clonal variants of the reference strain defining the cluster were found in 9 (12%) of the 74 clusters identified after the genotyping of 612 M. tuberculosis isolates by IS6110 restriction fragment length polymorphism analysis and mycobacterial interspersed repetitive units-variable-number tandem repeat typing. Clusters with microevolution events were epidemiologically supported and involved 4 to 9 cases diagnosed over a 1- to 5-year period. The IS6110 insertion sites from 16 representative isolates of reference and microevolved variants were mapped by ligation-mediated PCR in order to characterize the genetic background involved in microevolution. Both intragenic and intergenic IS6110 locations resulted from these microevolution events. Among those cases of IS6110 locations in intergenic regions which could have an effect on the regulation of adjacent genes, we identified the overexpression of cytochrome P450 in one microevolved variant using quantitative real-time reverse transcription-PCR. Our results help to define the frequency with which microevolution can be expected in M. tuberculosis transmission chains. They provide a snapshot of the genetic background of these subtle rearrangements and identify an event in which IS6110-mediated microevolution in an isogenic background has functional consequences.

Genomic stability over 9 years of an isoniazid resistant Mycobacterium tuberculosis outbreak strain in Sweden

In molecular epidemiological studies of drug resistant Mycobacterium tuberculosis (TB) in Sweden a large outbreak of an isoniazid resistant strain was identified, involving 115 patients, mainly from the Horn of Africa. During the outbreak period, the genomic pattern of the outbreak strain has stayed virtually unchanged with regard to drug resistance, IS6110 restriction fragment length polymorphism and spoligotyping patterns. Here we present the complete genome sequence analyses of the index isolate and two isolates sampled nine years after the index case as well as experimental data on the virulence of this outbreak strain. Even though the strain has been present in the community for nine years and passaged between patients at least five times in-between the isolates, we only found four single nucleotide polymorphisms in one of the later isolates and a small (4 amino acids) deletion in the other compared to the index isolate. In contrast to many other evolutionarily successful outbreak lineages (e.g. the Beijing lineage) this outbreak strain appears to be genetically very stable yet evolutionarily successful in a low endemic country such as Sweden. These findings further illustrate that the rate of genomic variation in TB can be highly strain dependent, something that can have important implications for epidemiological studies as well as development of resistance.

Novel hot spot of IS6110 insertion in Mycobacterium tuberculosis

We describe a hot spot for the insertion of IS6110 in Mycobacterium tuberculosis located in the area of region of difference 724 (RD724). Because RD724 defines sublineage 724 of M. tuberculosis, caution must be exercised when screening for RD724, as different polymorphisms can be observed in this region.

Bacterial genetic signatures of human social phenomena among M. tuberculosis from an Aboriginal Canadian population

Despite a widespread global distribution and highly variable disease phenotype, there is little DNA sequence diversity among isolates of Mycobacterium tuberculosis. In addition, many regional population genetic surveys have revealed a stereotypical structure in which a single clone, lineage, or clade makes up the majority of the population. It is often assumed that dominant clones are highly adapted, that is, the overall structure of M. tuberculosis populations is the result of positive selection. In order to test this assumption, we analyzed genetic data from extant populations of bacteria circulating in Aboriginal communities in Saskatchewan, Canada. Demographic parameters of the bacterial population were estimated from archival epidemiological data collected over approximately 130 years since the onset of epidemic tuberculosis in the host communities. Bacterial genetic data were tested against neutral theory expectations and the local evolutionary history of M. tuberculosis investigated by phylogenetic analysis. Our findings are not consistent with positive selection on the bacterial population. Instead, we uncovered founder effects persisting over decades and barriers to gene flow within the bacterial population. Simulation experiments suggested that a combination of these neutral influences could result in the stereotypical structure of M. tuberculosis populations. Some aspects of population structure were suggestive of background selection, and data were on the whole consistent with combined effects of population bottlenecks, subdivision, and background selection. Neutral phenomena, namely, bottlenecks and partitions within populations, are prominent influences on the evolution of M. tuberculosis and likely contribute to restricted genetic diversity observed within this species. Given these influences, a complex evolutionary model will be required to define the relative fitness of different M. tuberculosis lineages and, ultimately, to uncover the genetic basis for its success as a pathogen.

Multiple, independent, identical IS6110 insertions in Mycobacterium tuberculosis PPE genes

IS6110 is a transposable element found in Mycobacterium tuberculosis complex members. Regions of preferential IS6110 integration occur within the M. tuberculosis genome but the element has not previously been shown to exhibit any sequence-specific integration preferences. Here we provide evidence for multiple independent IS6110 insertions into identical, or near-identical, positions within the highly homologous 5' region of 3 PPE genes.

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.

A genomic library-based amplification approach (GL-PCR) for the mapping of multiple IS6110 insertion sites and strain differentiation of Mycobacterium tuberculosis

Evidence suggests that insertion of the IS6110 element is not without consequence to the biology of Mycobacterium tuberculosis complex strains. Thus, mapping of multiple IS6110 insertion sites in the genome of biomedically relevant clinical isolates would result in a better understanding of the role of this mobile element, particularly with regard to transmission, adaptability and virulence. In the present paper, we describe a versatile strategy, referred to as GL-PCR, that amplifies IS6110-flanking sequences based on the construction of a genomic library. M. tuberculosis chromosomal DNA is fully digested with HincII and then ligated into a plasmid vector between T7 and T3 promoter sequences. The ligation reaction product is transformed into Escherichia coli and selective PCR amplification targeting both 5' and 3' IS6110-flanking sequences are performed on the plasmid library DNA. For this purpose, four separate PCR reactions are performed, each combining an outward primer specific for one IS6110 end with either T7 or T3 primer. Determination of the nucleotide sequence of the PCR products generated from a single ligation reaction allowed mapping of 21 out of the 24 IS6110 copies of two 12 banded M. tuberculosis strains, yielding an overall sensitivity of 87,5%. Furthermore, by simply comparing the migration pattern of GL-PCR-generated products, the strategy proved to be as valuable as IS6110 RFLP for molecular typing of M. tuberculosis complex strains. Importantly, GL-PCR was able to discriminate between strains differing by a single IS6110 band.

Seeking completeness in bacterial mutant hunts

The identification of most or all of the genetic functions that are required for a particular biological process could be achieved through phenotypic studies of high genome-coverage mutant collections. Technologies for creating such collections, in the form of mixed populations or individually arrayed sequence-defined mutants, are now available for numerous bacterial species. The analysis of mixed mutant collections using microarray-based detection procedures appears to be particularly effective in identifying functions required for complex processes such as virulence. The phenotypic analysis of sequence-defined mutant libraries provides a virtually complete identification of nonessential genes required for processes for which suitable screens can be devised. Such libraries also serve as a source of individual mutants for examining the biological relevance of gene associations revealed by transcriptional profiling or homology.

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.