IS6110 transposition and evolutionary scenario of the direct repeat locus in a group of closely related Mycobacterium tuberculosis strains

Development, Evaluation, and Implementation of a House-Made Targeted Next-Generation Sequencing Spoligotyping in a French Laboratory

Epidemiological studies investigating transmission chains of tuberculosis are undertaken worldwide to tackle its spread. CRISPR locus diversity, called spoligotyping, is a widely used genotyping assay for Mycobacterium tuberculosis complex (MTBC) characterization. Herein, we developed a house-made targeted next-generation sequencing (tNGS) spoligotyping, and compared its outputs with those of membrane-based spoligotyping. A total of 144 clinical MTBC strains were retrospectively selected to be representative of the local epidemiology. Data analysis of a training set allowed for the setting of “presence”/“absence” thresholds for each spacer to maximize the sensibility and specificity related to the membrane-based spoligotyping. The thresholds above, in which the spacer was considered present, were 50 read per millions for spacers 10 and 14, 20,000 for spacers 20, 21, and 31, and 1000 for the other spacers. The confirmation of these thresholds was performed using a validation set. The overall agreement on the training and validation sets was 97.5% and 93.8%, respectively. The discrepancies concerned six strains: Two for spacer 14, two for spacer 31, and two for spacer 32. The tNGS spoligotyping, whose thresholds were finely-tuned during a careful bioinformatics pipeline development process, appears be a technique that is reliable, inexpensive, free of handling errors, and automatable through automatic transfer into the laboratory computer system.

Consistency of Mycobacterium tuberculosis Complex Spoligotyping between the Membrane-Based Method and In Silico Approach

To tackle the spread of tuberculosis (TB), epidemiological studies are undertaken worldwide to investigate TB transmission chains. Clustered regulatory interspaced short palindromic repeats (CRISPR) locus diversity, also called spoligotyping, is a widely used genotyping assay for the characterization of Mycobacterium tuberculosis complex (MTBC). We compared herein the spoligotyping of MTBC clinical isolates using a membrane-based method (following an initial PCR step) and whole-genome sequencing (WGS)-based method (i.e., in silico spoligotyping). All MTBC strains isolated at the Lyon University Hospital, France, between November 2016 and December 2020 were included (n = 597). Spoligotyping profiles were also used for species identification among the MTBC. Outputs of both methods were analyzed, and discrepant results were investigated thanks to CRISPRbuilder-TB. The overall agreement was 85.7%. Spacer discrepancies observed between the methods were due to the insertion of IS6110 within the direct repeat (DR) sequence upstream or downstream of spacers, mutated DR sequences, or truncated spacers. Discrepancies did not impact species identification. Although spoligotyping-based species identification was inconclusive for 29 isolates, SNP-based phylogeny conducted after WGS allowed the identification of 23 M. tuberculosis (Mtb), 2 M. canettii, and 4 mixed MTBC infections. WGS yielded very few discrepancies compared to membrane-based spoligotyping. Overall agreement was significantly improved (92.4%) by the CRISPR locus reconstruction using CRISPRbuilder-TB for the MTBC isolates with the shared international type 53 in silico spoligotyping. A smooth transition from the membrane-based to the in silico-based genotyping of M. tuberculosis isolates is, therefore, possible for TB diagnosis and epidemiologic survey. IMPORTANCE Whole-genome sequencing (WGS) has profoundly transformed the perspectives of tuberculosis (TB) diagnosis, providing a better discriminatory power to determine relatedness between Mycobacterium tuberculosis complex (MTBC) isolates. Previous genotyping approaches, such as spoligotyping consisting of an initial PCR step followed by reverse dot hybridization, are currently being replaced by WGS. Several pipelines have been developed to extract a spoligotype from WGS data (in silico spoligotyping) allowing for the continuity of MTBC molecular surveys before and after WGS implementation. The present study found very good overall agreement between hybridization to membrane-based spoligotyping and in silico spoligotyping, indicating the possibility of a smooth transition from the traditional to the in silico-based genotyping of MTBC isolates for TB diagnosis and epidemiological survey.

The role of IS6110 in micro- and macroevolution of Mycobacterium tuberculosis lineage 2

The insertion sequence 6110 (IS6110) is the most studied transposable element in the Mycobacterium tuberculosis complex species. The element plays a significant role in genome plasticity of this important human pathogen, but still many causes and consequences of its transposition have not been fully studied. Here, we analyzed insertion sites for 902 Mycobacterium tuberculosis lineage 2 strains using whole-genome sequencing data. In total, 17,972 insertions were found, corresponding to 827 independent positions in the genome of the reference strain H37Rv. To trace the history of IS6110 expansion since proto-Beijing strains up to modern sublineages, we looked at the distribution of IS6110 across the genome-wide SNP-based phylogenetic tree. This analysis demonstrated a stepwise transposition of IS6110 that occurs by «copy-and-paste» mechanism. Additionally, we detected evolutionary-scale and sublineage-specific integration sites, which can be used for typing and for understanding the reasons for the success of the lineage. A significant part of such insertions affected the genes that are essential for the pathogen. Finally, we identified and confirmed deletions that occurred between differently oriented elements, which is uncommon for this family of insertion elements and appears to be another mechanism of genome variability.

Long reads reveal the diversification and dynamics of CRISPR reservoir in microbiomes

Background

Sequencing of microbiomes has accelerated the characterization of the diversity of CRISPR-Cas immune systems. However, the utilization of next generation short read sequences for the characterization of CRISPR-Cas dynamics remains limited due to the repetitive nature of CRISPR arrays. CRISPR arrays are comprised of short spacer segments (derived from invaders’ genomes) interspaced between flanking repeat sequences. The repetitive structure of CRISPR arrays poses a computational challenge for the accurate assembly of CRISPR arrays from short reads. In this paper we evaluate the use of long read sequences for the analysis of CRISPR-Cas system dynamics in microbiomes.

Results

We analyzed a dataset of Illumina’s TruSeq Synthetic Long-Reads (SLR) derived from a gut microbiome. We showed that long reads captured CRISPR spacers at a high degree of redundancy, which highlights the spacer conservation of spacer sharing CRISPR variants, enabling the study of CRISPR array dynamics in ways difficult to achieve though short read sequences. We introduce compressed spacer graphs, a visual abstraction of spacer sharing CRISPR arrays, to provide a simplified view of complex organizational structures present within CRISPR array dynamics. Utilizing compressed spacer graphs, several key defining characteristics of CRISPR-Cas system dynamics were observed including spacer acquisition and loss events, conservation of the trailer end spacers, and CRISPR arrays’ directionality (transcription orientation). Other result highlights include the observation of intense array contraction and expansion events, and reconstruction of a full-length genome for a potential invader (Faecalibacterium phage) based on identified spacers.

Conclusion

We demonstrate in an in silico system that long reads provide the necessary context for characterizing the organization of CRISPR arrays in a microbiome, and reveal dynamic and evolutionary features of CRISPR-Cas systems in a microbial population.

Structure and variation of CRISPR and CRISPR-flanking regions in deleted-direct repeat region Mycobacterium tuberculosis complex strains

Background

CRISPR and CRISPR-flanking genomic regions are important for molecular epidemiology of Mycobacterium tuberculosis complex (MTBC) strains, and potentially for adaptive immunity to phage and plasmid DNA, and endogenous roles in the bacterium. Genotyping in the Israel National Mycobacterium Reference Center Tel-Aviv of over 1500 MTBC strains from 2008–2013 showed three strains with validated negative 43-spacer spoligotypes, that is, with putatively deleted direct repeat regions (deleted-DR/CRISPR regions). Two isolates of each of three negative spoligotype MTBC (a total of 6 isolates) were subjected to Next Generation Sequencing (NGS). As positive controls, NGS was performed for three intact-DR isolates belonging to T3_Eth, the largest multiple-drug-resistant (MDR)-containing African-origin cluster in Israel. Other controls consisted of NGS reads and complete whole genome sequences from GenBank for 20 intact-DR MTBC and for 1 deleted-DR MTBC strain recognized as CAS by its defining RD deletion.

Results

NGS reads from negative spoligotype MTBC mapped to reference H37Rv NC_000962.3 suggested that the DR/CRISPR regions were completely deleted except for retention of the middle IS6110 mobile element. Clonally specific deletion of CRISPR-flanking genes also was observed, including deletion of at least cas2 and cas1 genes. Genomic RD deletions defined lineages corresponding to the major spoligotype families Beijing, EAI, and Haarlem, consistent with 24 loci MIRU-VNTR profiles. Analysis of NGS reads, and analysis of contigs obtained by manual PCR confirmed that all 43 gold standard DR/CRISPR spacers were missing in the deleted-DR genomes.

Conclusions

Although many negative spoligotype strains are recorded as spoligotype-international-type (SIT) 2669 in the SITVIT international database, this is the first time to our knowledge that it has been shown that negative spoligotype strains are found in at least 4 different 24 loci MIRU-VNTR and RD deletion families. We report for the first time negative spoligotype-associated total loss of CRISPR region spacers and repeats, with accompanying clonally specific loss of flanking genes, including at least CRISPR-associated genes cas2 and cas1. Since cas1 deleted E.coli shows increased sensitivity to DNA damage and impaired chromosomal segregation, we discussed the possibility of a similar phenotype in the deleted-DR strains and Beijing family strains as both lack the cas1 gene.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3560-6) contains supplementary material, which is available to authorized users.

Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria

Molecular typing has revolutionized epidemiological studies of infectious diseases, including those of a mycobacterial etiology. With the advent of fingerprinting techniques, many traditional concepts regarding transmission, infectivity, or pathogenicity of mycobacterial bacilli have been revisited, and their conventional interpretations have been challenged. Since the mid-1990s, when the first typing methods were introduced, a plethora of other modalities have been proposed. So-called molecular epidemiology has become an essential subdiscipline of modern mycobacteriology. It serves as a resource for understanding the key issues in the epidemiology of tuberculosis and other mycobacterial diseases. Among these issues are disclosing sources of infection, quantifying recent transmission, identifying transmission links, discerning reinfection from relapse, tracking the geographic distribution and clonal expansion of specific strains, and exploring the genetic mechanisms underlying specific phenotypic traits, including virulence, organ tropism, transmissibility, or drug resistance. Since genotyping continues to unravel the biology of mycobacteria, it offers enormous promise in the fight against and prevention of the diseases caused by these pathogens. In this review, molecular typing methods for Mycobacterium tuberculosis and nontuberculous mycobacteria elaborated over the last 2 decades are summarized. The relevance of these methods to the epidemiological investigation, diagnosis, evolution, and control of mycobacterial diseases is discussed.

The Evolutionary History, Demography, and Spread of the Mycobacterium tuberculosis Complex

With the advent of next-generation sequencing technology, the genotyping of clinical Mycobacterium tuberculosis strains went through a major breakup that dramatically improved the field of molecular epidemiology but also revolutionized our deep understanding of the M. tuberculosis complex evolutionary history. The intricate paths of the pathogen and its human host are reflected by a common geographical origin in Africa and strong biogeographical associations that largely reflect the past migration waves out of Africa. This long coevolutionary history is cardinal for our understanding of the host-pathogen dynamic, including past and ongoing demographic components, strains’ genetic background, as well as the immune system genetic architecture of the host. Coalescent- and Bayesian-based analyses allowed us to reconstruct population size changes of M. tuberculosis through time, to date the most recent common ancestor and the several phylogenetic lineages. This information will ultimately help us to understand the spread of the Beijing lineage, the rise of multidrug-resistant sublineages, or the fall of others in the light of socioeconomic events, antibiotic programs, or host population densities. If we leave the present and go through the looking glass, thanks to our ability to handle small degraded molecules combined with targeted capture, paleomicrobiology covering the Pleistocene era will possibly unravel lineage replacements, dig out extinct ones, and eventually ask for major revisions of the current model.

Spoligotyping of Mycobacterium tuberculosis Complex Isolates by Use of Ligation-Based Amplification and Melting Curve Analysis

We report here a ligation-based spoligotyping that can identify unamplified spacers in membrane-based spoligotyping due to asymmetric insertion of IS6110 in the direct repeat locus. Our typing yielded 84.4% (411/487) concordance with traditional typing and 100% (487/487) accuracy when confirmed by DNA sequencing.

Overview and phylogeny of Mycobacterium tuberculosis complex organisms: implications for diagnostics and legislation of bovine tuberculosis

Members of the Mycobacterium tuberculosis complex (MTBC) cause a serious disease with similar pathology, tuberculosis; in this review, bovine tuberculosis will be considered as disease caused by any member of the MTBC in bovids. Bovine tuberculosis is responsible for significant economic loss due to costly eradication programs and trade limitations and poses a threat to both endangered and protected species as well as to public health. We here give an overview on all members of the MTBC, focusing on their isolation from different animal hosts. We also review the recent advances made in elucidating the evolutionary and phylogenetic relationships of members of the MTBC. Because the nomenclature of the MTBC is controversial, its members have been considered species, subspecies or ecotypes, this review discusses the possible implications for diagnostics and the legal consequences of naming of new species.

MLVA based classification of Mycobacterium tuberculosis complex lineages for a robust phylogeographic snapshot of its worldwide molecular diversity

Multiple-locus variable-number tandem repeat analysis (MLVA) is useful to establish transmission routes and sources of infections for various microorganisms including Mycobacterium tuberculosis complex (MTC). The recently released SITVITWEB database contains 12-loci Mycobacterial Interspersed Repetitive Units – Variable Number of Tandem DNA Repeats (MIRU-VNTR) profiles and spoligotype patterns for thousands of MTC strains; it uses MIRU International Types (MIT) and Spoligotype International Types (SIT) to designate clustered patterns worldwide. Considering existing doubts on the ability of spoligotyping alone to reveal exact phylogenetic relationships between MTC strains, we developed a MLVA based classification for MTC genotypic lineages. We studied 6 different subsets of MTC isolates encompassing 7793 strains worldwide. Minimum spanning trees (MST) were constructed to identify major lineages, and the most common representative located as a central node was taken as the prototype defining different phylogenetic groups. A total of 7 major lineages with their respective prototypes were identified: Indo-Oceanic/MIT57, East Asian and African Indian/MIT17, Euro American/MIT116, West African-I/MIT934, West African-II/MIT664, M. bovis/MIT49, M.canettii/MIT60. Further MST subdivision identified an additional 34 sublineage MIT prototypes. The phylogenetic relationships among the 37 newly defined MIRU-VNTR lineages were inferred using a classification algorithm based on a bayesian approach. This information was used to construct an updated phylogenetic and phylogeographic snapshot of worldwide MTC diversity studied both at the regional, sub-regional, and country level according to the United Nations specifications. We also looked for IS6110 insertional events that are known to modify the results of the spoligotyping in specific circumstances, and showed that a fair portion of convergence leading to the currently observed bias in phylogenetic classification of strains may be traced back to the presence of IS6110. These results shed new light on the evolutionary history of the pathogen in relation to the history of peopling and human migration.

Massive gene acquisitions in Mycobacterium indicus pranii provide a perspective on mycobacterial evolution

Understanding the evolutionary and genomic mechanisms responsible for turning the soil-derived saprophytic mycobacteria into lethal intracellular pathogens is a critical step towards the development of strategies for the control of mycobacterial diseases. In this context, Mycobacterium indicus pranii (MIP) is of specific interest because of its unique immunological and evolutionary significance. Evolutionarily, it is the progenitor of opportunistic pathogens belonging to M. avium complex and is endowed with features that place it between saprophytic and pathogenic species. Herein, we have sequenced the complete MIP genome to understand its unique life style, basis of immunomodulation and habitat diversification in mycobacteria. As a case of massive gene acquisitions, 50.5% of MIP open reading frames (ORFs) are laterally acquired. We show, for the first time for Mycobacterium, that MIP genome has mosaic architecture. These gene acquisitions have led to the enrichment of selected gene families critical to MIP physiology. Comparative genomic analysis indicates a higher antigenic potential of MIP imparting it a unique ability for immunomodulation. Besides, it also suggests an important role of genomic fluidity in habitat diversification within mycobacteria and provides a unique view of evolutionary divergence and putative bottlenecks that might have eventually led to intracellular survival and pathogenic attributes in mycobacteria.

Use of mouse models to study the variability in virulence associated with specific genotypic lineages of Mycobacterium tuberculosis

The host response against Mycobacterium tuberculosis show a wide spectrum of clinical manifestations in those patients who fail to control the infection. The course of the infection and its epidemiological consequences depend upon a complex interplay of host, environmental and bacterial factors. Experimental animal models have helped to define the influence of bacterial genetic diversity on virulence and on the immune response that is induced. For this purpose, experimental animals such as mice, guinea pigs and rabbits have been infected with selected clinical isolates obtained from outbreaks or from clinical epidemiology settings. Here we review the contribution of mouse models to defining the variability in virulence and immune response in relation to mycobacterial genetic diversity. Low dose aerosol infection in C57Bl mice or high dose intratracheal infection in BALB/c mice have demonstrated wide variability in virulence and immune responses induced by different bacterial genotypes, and each genotype has different phenotypes, with high and low virulence variants. In general, these studies have shown that high prevalent strains from big clusters are more virulent than low prevalent sporadic clinical isolates, and highly virulent strains induce non-protective immune responses with some correlation with clinical-epidemiological data. In the future selected strains from these types of studies should be analyzed with molecular technologies. These kinds of study will contribute to the identification of mycobacterial genes associated with virulence and immunogenicity.

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.

European 2--a clonal complex of Mycobacterium bovis dominant in the Iberian Peninsula

Mycobacterium bovis isolates from the Iberian Peninsula are dominated by strains with spoligotype patterns deleted for spacer 21. Whole-genome sequencing of three Spanish strains with spacer 21 missing in their spoligotype pattern revealed a series of SNPs and subsequent screening of a selection of these SNPs identified one in gene guaA that is specific to these strains. This group of strains from the Iberian Peninsula missing spoligotype spacer 21 represents a new clonal complex of M. bovis, defined by the SNP profile with a distinct spoligotype signature. We have named this clonal complex European 2 (Eu2) and found that it was present at low frequency in both France and Italy and absent from the British Isles.

Web tools for molecular epidemiology of tuberculosis

In this study we explore publicly available web tools designed to use molecular epidemiological data to extract information that can be employed for the effective tracking and control of tuberculosis (TB). The application of molecular methods for the epidemiology of TB complement traditional approaches used in public health. DNA fingerprinting methods are now routinely employed in TB surveillance programs and are primarily used to detect recent transmissions and in outbreak investigations. Here we present web tools that facilitate systematic analysis of Mycobacterium tuberculosis complex (MTBC) genotype information and provide a view of the genetic diversity in the MTBC population. These tools help answer questions about the characteristics of MTBC strains, such as their pathogenicity, virulence, immunogenicity, transmissibility, drug-resistance profiles and host-pathogen associativity. They provide an integrated platform for researchers to use molecular epidemiological data to address current challenges in the understanding of TB dynamics and the characteristics of MTBC.

A fragment of 21 ORFs around the direct repeat (DR) region of Mycobacterium tuberculosis is absent from the other sequenced mycobacterial genomes: implications for the evolution of the DR region

The direct repeat (DR) region is a singular locus of the Mycobacterium tuberculosis complex genome. This region consists of 36 bp repetitive sequences separated by non-repetitive unique spacer sequences. Around this region there are several genes coding for proteins of unknown function. To determine whether the M. smegmatis, M. avium, M. marinum and M. leprae genomes contain sequences and ORFs similar to those of the DR locus of the M. tuberculosis complex, we analysed the corresponding regions in these species. As a first step, some conserved genes that flank the DR genes [Rv2785c (rpsO), Rv2786c (ribF), Rv2790c (ltp1 ), Rv2793c (truB), Rv2800, Rv2825, Rv2828, Rv2831 (echA16 ), Rv2838 (rbfA) and Rv2845 (proS )] were used as markers to locate the corresponding orthologues in M. smegmatis, M. avium, M. marinum and M. leprae in silico. Most of these M. tuberculosis marker genes have highly similar orthologues located in the same order and orientation in the other mycobacteria. In contrast, no orthologues were found for ORFs Rv2801–Rv2824, suggesting that these genes are unique to M. tuberculosis within the genus Mycobacterium.We observed that in M. smegmatis and M. avium, Rv2800 and Rv2825 are adjacent. This observation was experimentally confirmed by PCR. In conclusion, as the DR locus and the ORFs around it are absent in M. smegmatis and M. avium and, as it is possible that these species are older than M. tuberculosis, we postulated that the DR locus was acquired by the M. tuberculosis complex species or by an ancestor bacterium.

A review of the contemporary knowledge of bovine tuberculosis and government policy in Iran

In 1931 Carpantier reported bovine TB (BTB) in Iranian cattle. Some eighty years on with a national test-and-slaughter programme in place for over four decades, the efforts to vanquish Mycobacterium bovis (M. bovis) infection in cattle have been in vain as the vast majority of the 30 Iranian provinces still have reports of BTB in their cattle herds every year. This paper reviews the present epidemiology of BTB in Iran and in the region and evaluates the success of government policy in controlling this disease.

Preferential deletion events in the direct repeat locus of Mycobacterium tuberculosis

The "Harlingen" IS6110 restriction fragment length polymorphism (RFLP) cluster has linked over 100 tuberculosis cases in The Netherlands since 1993. Four Mycobacterium tuberculosis isolates that were epidemiologically linked to this cluster had different spoligotype patterns, as well as slightly divergent IS6110 profiles, compared to the majority of the isolates. Sequencing of the direct repeat (DR) locus revealed sequence polymorphisms at the putative deletion sites. These deletion footprints provided evidence for independent deletions of the central region of the DR locus in three isolates, while the different genotype of the fourth isolate was explained by transmission. Our finding suggests that convergent deletions in the DR locus occur frequently. However, deletion footprints are not suitable to detect convergent deletions in the DR because they seem to be exceptional. Deletion footprints in the DR were not described previously, and we did not observe them in any public M. tuberculosis complex sequences. We conclude that preferential deletions in the DR loci of closely related strains are usually an unnoted event that interferes with clustering of closely related strains.

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

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

Replacing reverse line blot hybridization spoligotyping of the Mycobacterium tuberculosis complex

Spoligotyping is a tool for the molecular characterization/typing of Mycobacterium tuberculosis complex (MTBC) strains based on target sequences (spacers) in the direct repeat (DR) region (14). The standard spoligotyping assay involves the hybridization of amplified sample DNA to nylon membrane-immobilized oligonucleotides whose sequences are representative of 43 spacer regions. Variations in the number of spacers as a result of deletions of adjacent blocks of repetitive units allow the differentiation of clinical isolates. In the present study, we developed a new multiplexed primer extension-based spoligotyping assay using automated matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) that improves the classical reverse line blot hybridization assay with respect to reproducibility, throughput, process flow, ease of use, and data analysis. Validation of the MALDI-TOF MS-based spoligotyping assay with two sample sets with a total of 326 samples resulted in 96.6% concordance (315/326) when the full spoligotype patterns were compared with the results of standard spoligotyping and 99.9% concordance when the results were compared with those of individual primer extension assays. Ten strains (including two Mycobacterium canettii strains) showed discordant results with one or two spacer differences from the membrane-based spoligotyping result. Most discordant samples were identified to be the result of ambiguities in the interpretation of weak hybridization signals in the reverse line blot assay and sequence variations in the spacer regions. We established a new automated primer extension assay and successfully validated it for use for the routine typing of MTBC strains in the research and public health laboratory environments. The present multiplex levels of up to 30 are extendable and allow the additional incorporation of controls and antibiotic resistance markers.

The CRISPR system: small RNA-guided defense in bacteria and archaea

All cellular systems evolve ways to combat predators and genomic parasites. In bacteria and archaea, numerous resistance mechanisms have developed against phage. Our understanding of this defensive repertoire has recently been expanded to include the CRISPR system of clustered, regularly interspaced short palindromic repeats. In this remarkable pathway, short sequence tags from invading genetic elements are actively incorporated into the host's CRISPR locus to be transcribed and processed into a set of small RNAs that guide the destruction of foreign genetic material. Here we review the inner workings of this adaptable and heritable immune system and draw comparisons to small RNA-guided defense mechanisms in eukaryotic cells.

Comparison of restriction fragment length polymorphism with the polymorphic guanine-cytosine-rich sequence and spoligotyping for differentiation of Mycobacterium tuberculosis isolates with five or fewer copies of IS6110

The use of IS6110 as a marker for molecular epidemiological studies is limited when a Mycobacterium tuberculosis isolate has five or fewer copies of IS6110. Restriction fragment length polymorphism analysis with a highly polymorphic GC-rich repetitive sequence located in the plasmid pTBN12 (PGRS RFLP) and spoligotyping (based on the polymorphism of the DR region) are two frequently used secondary typing methods. The aim of this study was to compare the performance of these two methods in a population-based study in San Francisco. We included all patients with culture-positive tuberculosis from 1999 to 2007 with IS6110 RFLP results presenting five or fewer bands. PGRS RFLP and spoligotyping were performed using standardized methods. We determined the concordance between the two methods regarding cluster status and the risk factors for an isolate to be in a cluster with each of the methods. Our data indicate that both methods had similar discriminatory power and that the risk factors associated with clustering by either method were the same. Although the cluster/unique status was concordant in 84% of the isolates, patients were clustered differently depending on the method. Therefore, the methods are not interchangeable, and the same method should be used for longitudinal studies.

Mycobacterium bovis infection in Holstein Friesian cattle, Iran

To identify strains of Mycobacterium bovis circulating in Iran, we used region of difference, spoligotypes, and variable number tandem repeats to genotype 132 M. bovis isolates from Holstein Friesian cattle. Despite wide geographic origins, the strains were genetically homogeneous. Increased distribution of cattle herds and inadequate control measures may have contributed to strain dispersion.

African 1, an epidemiologically important clonal complex of Mycobacterium bovis dominant in Mali, Nigeria, Cameroon, and Chad

We have identified a clonal complex of Mycobacterium bovis present at high frequency in cattle in population samples from several sub-Saharan west-central African countries. This closely related group of bacteria is defined by a specific chromosomal deletion (RDAf1) and can be identified by the absence of spacer 30 in the standard spoligotype typing scheme. We have named this group of strains the African 1 (Af1) clonal complex and have defined the spoligotype signature of this clonal complex as being the same as the M. bovis BCG vaccine strain but with the deletion of spacer 30. Strains of the Af1 clonal complex were found at high frequency in population samples of M. bovis from cattle in Mali, Cameroon, Nigeria, and Chad, and using a combination of variable-number tandem repeat typing and spoligotyping, we show that the population of M. bovis in each of these countries is distinct, suggesting that the recent mixing of strains between countries is not common in this area of Africa. Strains with the Af1-specific deletion (RDAf1) were not identified in M. bovis isolates from Algeria, Burundi, Ethiopia, Madagascar, Mozambique, South Africa, Tanzania, and Uganda. Furthermore, the spoligotype signature of the Af1 clonal complex has not been identified in population samples of bovine tuberculosis from Europe, Iran, and South America. These observations suggest that the Af1 clonal complex is geographically localized, albeit to several African countries, and we suggest that the dominance of the clonal complex in this region is the result of an original introduction into cows naïve to bovine tuberculosis.

Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes

BACKGROUND

Molecular typing methods are commonly used to study genetic relationships among bacterial isolates. Many of these methods have become standardized and produce portable data. A popular approach for analyzing such data is to construct graphs, including phylogenies. Inferences from graph representations of data assist in understanding the patterns of transmission of bacterial pathogens, and basing these graph constructs on biological models of evolution of the molecular marker helps make these inferences. Spoligotyping is a widely used method for genotyping isolates of Mycobacterium tuberculosis that exploits polymorphism in the direct repeat region. Our goal was to examine a range of models describing the evolution of spoligotypes in order to develop a visualization method to represent likely relationships among M. tuberculosis isolates.

RESULTS

We found that inferred mutations of spoligotypes frequently involve the loss of a single or very few adjacent spacers. Using a second-order variant of Akaike's Information Criterion, we selected the Zipf model as the basis for resolving ambiguities in the ancestry of spoligotypes. We developed a method to construct graphs of spoligotypes (which we call spoligoforests). To demonstrate this method, we applied it to a tuberculosis data set from Cuba and compared the method to some existing methods.

CONCLUSION

We propose a new approach in analyzing relationships of M. tuberculosis isolates using spoligotypes. The spoligoforest recovers a plausible history of transmission and mutation events based on the selected deletion model. The method may be suitable to study markers based on loci of similar structure from other bacteria. The groupings and relationships in the spoligoforest can be analyzed along with the clinical features of strains to provide an understanding of the evolution of spoligotypes.

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.

On-line resources for bacterial micro-evolution studies using MLVA or CRISPR typing

The control of bacterial pathogens requires the development of tools allowing the precise identification of strains at the subspecies level. It is now widely accepted that these tools will need to be DNA-based assays (in contrast to identification at the species level, where biochemical based assays are still widely used, even though very powerful 16S DNA sequence databases exist). Typing assays need to be cheap and amenable to the designing of international databases. The success of such subspecies typing tools will eventually be measured by the size of the associated reference databases accessible over the internet. Three methods have shown some potential in this direction, the so-called spoligotyping assay (Mycobacterium tuberculosis, 40,000 entries database), Multiple Loci Sequence Typing (MLST; up to a few thousands entries for the more than 20 bacterial species), and more recently Multiple Loci VNTR Analysis (MLVA; up to a few hundred entries, assays available for more than 20 pathogens). In the present report we will review the current status of the tools and resources we have developed along the past seven years to help in the setting-up or the use of MLVA assays or lately for analysing Clustered Regularly Interspaced Short Palindromic Repeats called CRISPRs which are the basis for spoligotyping assays.

Rapidly evolving CRISPRs implicated in acquired resistance of microorganisms to viruses

Recent experimental evidence has demonstrated that bacteria acquire resistance to viruses by incorporation of short transcribed nucleotide sequences into regions of clustered regularly interspaced short palindromic repeats (CRISPR). We have analysed community genomic data from acidophilic microbial biofilms and discovered that evolution of the CRISPR regions in two distinct Leptospirillum group II bacteria occurs fast enough to promote individuality in otherwise nearly clonal populations. Comparative genomics strongly indicates very recent lateral transfer of the CRISPR locus between these populations, followed by significant loss of spacer sequences and locus expansion by unidirectional heterogeneous addition of new spacer sequences. Diversification of the CRISPR region is inferred to be a population-level response to the rapidly changing selective pressure of phage predation. Results reinforce the importance of phage-host interactions in shaping microbial ecology and evolution over very short time scales.

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.

Improvement of spoligotyping with additional spacer sequences for characterization of Mycobacterium bovis and M. caprae isolates from Spain

Spoligotyping is a typing tool used worldwide for epidemiological studies on Mycobacterium tuberculosis complex organisms; however it has received little attention regarding improvement of its discriminatory power (DP). In this study, we have evaluated a spoligotyping membrane prepared with 25 novel spacer sequences selected from a previous study [van der Zanden AG, Kremer K, Schouls LM. Improvement of differentiation and interpretability of spoligotyping for Mycobacterium tuberculosis complex isolates by introduction of new spacer oligonucleotides. J Clin Microbiol 2002;40:4628-39] on 308 M. bovis and 88 M. caprae Spanish isolates in comparison with the traditional spoligotyping membrane. The results obtained by combining the two membranes together revealed an improvement of 45 patterns instead of 31. The spacers used in the second membrane were able to distinguish 8 out of the 16 M. bovis types that had more than one isolate. Seven of these types were differentiated into two subtypes with the second-generation membrane, while spb-7, the most prevalent in Spain, was further differentiated into eight subtypes. This second-generation membrane also differentiates M. bovis from M. caprae. A set of 39 spacers (1, 2, 4-8, 10-15, 17-21, 23, 26-32, 37, 44-49, 51-54, 56 and 57) contain all the DP for both M. bovis and M. caprae isolates; and a set of 35 spacers (1, 2, 4-8, 10-15, 17-21, 26-32, 37, 44-48, 52-54 and 57) had all the DP for the M. bovis isolates. Our results show that the research on new spacers and the design of a new membrane may be useful for epidemiological studies of M. bovis and M. caprae isolates.

Molecular evolution of Mycobacterium tuberculosis

Tuberculosis continues to be the main cause of death from a single infectious agent in developing countries. The causative agent, Mycobacterium tuberculosis, is thought to have diverged from its common ancestor as recently as 15,000 years ago. Subsequently, various genetic elements have evolved over time at different rates and can be used to elucidate patterns of infection. When individual elements are studied within genetic families, very low rates of variation are observed for almost every marker. For example, when all M. tuberculosis genetic families are considered, the number of alleles observed at each mycobacterial interspersed repetitive unit (MIRU) locus usually drops when viewed within a single genetic family, indicating that the rate of repeat variation may be low, as each member of that family is a descendant of a single common ancestor. Also, the low level of silent nucleotide variation observed indicates that M. tuberculosis is, in evolutionary terms, very young. Mapping the variation of the different markers used in molecular epidemiology within a genetic framework enables the relative rates of variation of these markers to be determined and, together with a complete chronology, allows the identification of more informative panels of markers tailored to individual genetic families.

Microevolution of the direct repeat locus of Mycobacterium tuberculosis in a strain prevalent in San Francisco

We describe a microevolutionary event of a prevalent strain of Mycobacterium tuberculosis that caused two outbreaks in San Francisco. During the second outbreak, a direct variable repeat was lost. We discuss the mechanisms of this change and the implications of analyzing multiple genetic loci in this context.

Detecting emerging strains of tuberculosis by using spoligotypes

The W-Beijing strain of tuberculosis has been identified in many molecular epidemiological studies as being particularly prevalent. This identification has been made possible through the development of a number of genotyping technologies including spoligotyping. Highly prevalent genotypes associated with outbreaks, such as the W-Beijing strain, are implicitly regarded as fast spreading. Here we present a quantitative method to identify "emerging" strains, those that are spreading faster than the background rate inferred from spoligotype data. The approach uses information about the mutation process specific to spoligotypes, combined with a model of both transmission and mutation. The core principle is that if two comparable strains have the same number of isolates, then the strain with fewer inferred mutation events must have spread faster if the mutation process is common. Applying this method to four different data sets, we find not only the W-Beijing strain, but also a number of other strains, to be emerging in this sense. Importantly, the strains that are identified as emerging are not simply those with the largest number of cases. The use of this method should facilitate the targeting of individual genotypes in intervention programs.

Bottlenecks and broomsticks: the molecular evolution of Mycobacterium bovis

Mycobacterium bovis is the cause of tuberculosis in cattle and is a member of the Mycobacterium tuberculosis complex. In contrast to many other pathogenic bacterial species, there is little evidence for the transfer and recombination of genes between cells. The clonality of this group of organisms indicates that the population structure is dominated by reductions in diversity, caused either by population bottlenecks or selective sweeps as entire chromosomes become fixed in the population. We describe how these forces have shaped not only the phylogeny of this group but also, at a very local level, the population structure of Mycobacterium bovis in the British Isles. We also discuss the practical implications of applying this knowledge to understanding the spread of infection and the development of improved vaccines and diagnostic tests.

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.

Can 15-locus mycobacterial interspersed repetitive unit-variable-number tandem repeat analysis provide insight into the evolution of Mycobacterium tuberculosis?

The phylogeny and evolution of the bacterium Mycobacterium tuberculosis is still poorly understood despite the application of a variety of molecular techniques. We analyzed 469 M. tuberculosis and 49 Mycobacterium bovis isolates to evaluate if the mycobacterial interspersed repetitive units-variable-number tandem repeats (MIRU-VNTR) commonly used for epidemiological studies can define the phylogeny of the M. tuberculosis complex. This population was characterized by previously identified silent single-nucleotide polymorphisms (sSNPs) or by a macroarray based on these sSNPs that was developed in this study. MIRU-VNTR phylogenetic codes capable of differentiating between phylogenetic lineages were identified. Overall, there was 90.9% concordance between the lineages of isolates as defined by the MIRU-VNTR and sSNP analyses. The MIRU-VNTR phylogenetic code was unique to M. bovis and was not observed in any M. tuberculosis isolates. The codes were able to differentiate between different M. tuberculosis strain families such as Beijing, Delhi, and East African-Indian. Discrepant isolates with similar but not identical MIRU-VNTR codes often displayed a stepwise trend suggestive of bidirectional evolution. A lineage-specific panel of MIRU-VNTR can be used to subdivide each lineage for epidemiological purposes. MIRU-VNTR is a valuable tool for phylogenetic studies and could define an evolutionarily uncharacterized population of M. tuberculosis complex organisms.

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.

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

Transposable elements can affect an organism's fitness through the insertional inactivation of genes and can therefore be used to identify genes that are nonessential for growth in vitro or in animal models. However, these models may not adequately represent the genetic requirements during chains of human infection. We have therefore conducted a genome-wide survey of transposon mutations in Mycobacterium tuberculosis isolates from cases of human infection, identifying the precise, base-specific insertion sites of the naturally occurring transposable element IS6110. Of 294 distinct insertions mapped to the strain H37Rv genome, 180 were intragenic, affecting 100 open reading frames. The number of genes carrying IS6110 in clinical isolates, and hence apparently not essential for infection and transmission, is very much lower than the estimates of nonessential genes derived from in vitro studies. This suggests that most genes in M. tuberculosis play a significant role in human infection chains. IS6110 insertions were underrepresented in genes associated with virulence, information pathways, lipid metabolism, and membrane proteins but overrepresented in multicopy genes of the PPE family, genes of unknown function, and intergenic sequences. Population genomic analysis of isolates recovered from an organism's natural habitat is an important tool for determining the significance of genes or classes of genes in the natural biology of an organism.

Ecotypes of the Mycobacterium tuberculosis complex

A phylogeny of the Mycobacterium tuberculosis complex has recently shown that the animal-adapted strains are found in a single lineage marked by the deletion of chromosomal region 9 (RD9) [Brosch et al., 2002. A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc. Natl Acad. Sci. USA 99 (6), 3684-3689]. We have obtained the spoligotype patterns of the RD9 deleted strains used to generate this new evolutionary scenario and we show that the presence of spoligotype spacers 3, 9, 16, 39, and 40-43 is phylogenetically informative in this lineage. We have used the phylogenetically informative spoligotype spacers to screen a database of spoligotype patterns and have identified further members of a group of strains apparently host-adapted to antelopes. The presence of the spoligotype spacers is congruent with the phylogeny generated by chromosomal deletions, suggesting that recombination is rare or absent between strains of this lineage. The phylogenetically informative spacers, in concert with the previously identified single nucleotide mutations and chromosomal deletions, can be used to identify a series of clades in the RD9 deleted lineage each with a separate host preference. Finally, we discuss the application of the ecotype concept to this series of clades and suggest that the M. tuberculosis complex may best be described as a series of host-adapted ecotypes.

Clonal expansion of a globally disseminated lineage of Mycobacterium tuberculosis with low IS6110 copy numbers

Knowledge of the clonal expansion of Mycobacterium tuberculosis and accurate identification of predominant evolutionary lineages in this species remain limited, especially with regard to low-IS6110-copy-number strains. In this study, 170 M. tuberculosis isolates with </=6 IS6110 insertions identified in Cape Town, South Africa, were characterized by principal genetic grouping, restriction fragment length polymorphism analysis, spoligotyping, IS6110 insertion site mapping, and variable-number tandem repeat (VNTR) typing. These analyses indicated that all but one of the isolates analyzed were members of principal genetic group 2 and of the same low-IS6110-copy-number lineage. The remaining isolate was a member of principal genetic group 1 and a different low-IS6110-copy-number lineage. Phylogenetic reconstruction suggests clonal expansion through sequential acquisition of additional IS6110 copies, expansion and contraction of VNTR sequences, and the deletion of specific direct-variable-repeat sequences. Furthermore, comparison of the genotypic data of 91 representative low-IS6110-copy-number isolates from Cape Town, other southern African regions, Europe, and the United States suggests that certain low-IS6110-copy-number strain spoligotypes and IS6110 fingerprints were acquired in the distant past. These clones have subsequently become widely disseminated and now play an important role in the global tuberculosis epidemic.

Molecular characteristics of strains of the cameroon family, the major group of Mycobacterium tuberculosis in a country with a high prevalence of tuberculosis

A preliminary investigation of the genetic biodiversity of Mycobacterium tuberculosis complex strains in Cameroon, a country with a high prevalence of tuberculosis, described a group of closely related M. tuberculosis strains (the Cameroon family) currently responsible for more than 40% of smear-positive pulmonary tuberculosis cases. Here, we used various molecular methods to study the genetic characteristics of this family of strains. Cameroon family M. tuberculosis strains (i) are part of the major genetic group 2 and lack the TbD1 region like other families of epidemic strains, (ii) lack spacers 23, 24, and 25 in their direct repeat (DR) region, (iii) have an identical number of repeats in 8 of 12 variable-number tandem repeats of mycobacterial interspersed repetitive unit (MIRU-VNTR) loci, (iv) have similar IS6110-restriction fragment length polymorphism (RFLP) multiband patterns (10 to 15 copies) with seven common IS6110 bands, (v) do not have an IS6110 element in their DR locus, and (vi) have four IS6110 elements in open reading frames (adenylate cyclase, phospholipase C, moeY, and ATP binding genes). Analysis by spoligotyping, MIRU-VNTR, and IS6110-RFLP typing methods revealed differences not observed in previous studies; polymorphism as assessed by MIRU-VNTR typing was lower than suggested by spoligotyping, and in rare cases, strains with identical IS6110-RFLP patterns had spoligotypes differing by as much as 15 spacers. Our findings confirm the recent expansion of this family in Cameroon and indicate that the interpretation of molecular typing results has to be adapted to the characteristics of the strain population within each setting. The knowledge of this particular genotype, with its large involvement in tuberculosis in Cameroon, allows greater refinement of tuberculosis transmission studies by interpreting data in the context of this geographic area.

Spoligotyping profile change caused by deletion of a direct variable repeat in a Mycobacterium tuberculosis isogenic laboratory strain

Spoligotyping is a major tool for molecular typing of Mycobacterium tuberculosis complex organisms. For epidemiological purposes, strains are considered clonal only when their spoligotyping patterns are identical. We report a change in the spoligotyping profiles of truly isogenic strains (a clinical isolate and a subculture derived in the laboratory) caused by deletion of a direct variable repeat. Without the information about the relationship between them, a link between these strains would have gone unnoticed. Evolutionary events should be taken into account in the interpretation of spoligotyping results and in the design of databases.

Correlation of virulence, lung pathology, bacterial load and delayed type hypersensitivity responses after infection with different Mycobacterium tuberculosis genotypes in a BALB/c mouse model

One of the most intriguing aspects of tuberculosis is that the outcome of an infection with M. tuberculosis (TB) is highly variable between individuals. The possibility of differences in virulence between M. tuberculosis strains or genotypes has only recently been studied. There is evidence of multifactorial genetic predisposition in humans that influences the susceptibility to tuberculosis. A better understanding of differences in virulence between M. tuberculosis genotypes could be important with regard to the efforts at TB control and the development of improved antituberculosis vaccines. Survival, lung pathology, bacterial load and delayed type hypersensitivity (DTH) responses of BALB/c mice after intratracheal infection with any of 19 different M. tuberculosis complex strains of 11 major genotype families were studied. The results indicate that among genetically different M. tuberculosis strains a very broad response was present with respect to virulence, pathology, bacterial load and DTH. 'Low'-responders were the H37Rv, Canetti, Beijing-1 strains, while Beijing-2,3, Africa-2 and Somalia-2 strains were 'high'-responders. A severe pathological response correlates with a high mortality and a high CFU counts in lungs, but poorly with the degree of the DTH response.

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.

A new high-throughput AFLP approach for identification of new genetic polymorphism in the genome of the clonal microorganism Mycobacterium tuberculosis

We have here applied high-throughput amplified fragment length polymorphism (htAFLP) analysis to strains belonging to the five classical species of the Mycobacterium tuberculosis complex. Using 20 strains, three enzyme combinations and eight selective amplification primer pairs, 24 AFLP reactions were performed per strain. Overall, this resulted in 480 DNA fingerprints and more than 1200 htAFLP-amplified PCR fragments were visualised per strain. The cumulative dendrogram correctly clustered strains from the various species, albeit within a distance of 6.5% for most of them. The single isolate of Mycobacterium canettii presented separately at 19% distance. All over, 169 fragments (14%) appeared to be polymorphic. Sixty-eight were specific for M. canetti and forty-five for Mycobacterium bovis. For the 10 different M. tuberculosis strains included in the present analysis, 56 polymorphic markers were identified. Upon sequencing 20 of these marker regions and comparisons with the H37Rv genome sequence, 25% appeared to share homology to members of the antigenically variable PE/PPE surface protein encoding gene family confirming previous findings on the genetic heterogeneity within these genes. In addition, homologues for phage genes and insertion element-encoded genes were detected. Forty-five percent of the sequences derived from ORFs with a currently unknown function, which was corroborated by genome sequence comparison for the clinical M. tuberculosis CD 1551 isolate. Sequence variation in M. tuberculosis was assessed in more detail for a subset of these loci by newly designed PCR restriction fragment length polymorphism (RFLP) tests and direct sequencing. Fourteen novel PCR RFLP tests were developed and twelve novel single nucleotide polymorphisms (SNPs) were identified, all suited for epidemiological analysis of M. tuberculosis. The tests allowed for identification of the major Mycobacterium species and M. tuberculosis variants and clones.

Molecular differentiation of Mycobacterium bovis isolates. Review of main techniques and applications

Until recently, none of the Mycobacterium bovis typing techniques permitted a satisfactory differentiation of isolates. During the last 10 years, the genome of pathogenic mycobacteria has been extensively studied, and phylogenetic analyses have shown that all (except Mycobacterium avium) belong to a single genetic species: the Mycobacterium tuberculosis complex. This increase in knowledge about the genome of these bacteria has lead to the discovery of molecular markers that allow us to differentiate isolates. Because of the phylogenetic proximity of the strains, even if most of these markers have been discovered in M. tuberculosis, they could be successfully adapted to the other bacteria of the M. tuberculosis complex, especially M. bovis. The most common markers in use today are the IS6110 insertion sequence, the direct repeat (DR) region, the poly(GC) rich (PGRS) sequences and the variable number tandem repeats (VNTR) sequences. The corresponding typing techniques are briefly described, and current knowledge of polymorphism and marker stability is detailed. If molecular markers are to offer wide perspectives for field studies, these two characteristics (polymorphism and stability) must be taken into account when choosing the marker(s) used in a study. In this context, examples of the application of molecular typing techniques for M. bovis are reviewed, on the one hand with epidemiological studies for which the major problem is the comparison between isolates and, on the other, with more general studies about the population genetics of M. bovis in a given country, and about its history and its phylogeny.

Locating transposable element polymorphisms in bacterial genomes

Although whole-genome sequencing is greatly extending our knowledge of the genetic capacity of those bacterial species, it is only directly informative for the particular strain sequenced. Many bacterial species exhibit more or less genetic polymorphism within their populations and characterising this variety is an extremely important way of elucidating the biology of these species. Often genomic polymorphisms are associated with multicopy elements, particularly transposable elements. We describe a novel method that efficiently characterises the sequences of such polymorphisms. We have optimised heminested inverse PCR (hINVPCR) to assess the diversity of insertional polymorphisms of a transposable element (IS6110) in clinical isolates of Mycobacterium tuberculosis. To increase the yield of information, genomic DNA was digested with different endonucleases (Bsp1286I, HaeII or PvuI), and primers based on both the 5' and 3' ends of IS6110 were used to amplify and determine the genomic sequence upstream (or downstream) of the transposable element. We found that both the choice of restriction enzyme and the use of primers at both ends of the transposable element significantly increased the diversity of the insertion sites identified. Band stabbing was incorporated into the method as an alternative to cloning in order to screen large number of isolates at a sequence level in a rapid and labour-efficient fashion. We describe some of the purposes to which such data can be put.