Molecular evidence for independent occurrence of IS6110 insertions at the same sites of the genome of Mycobacterium tuberculosis in different clinical isolates

Clustured regularly interspersed short palindromic repeats (CRISPR) genetic diversity studies as a mean to reconstruct the evolution of the Mycobacterium tuberculosis complex

The natural history of tuberculosis may be tackled by various means, among which the record of molecular scars that have been registered by the Mycobacterium tuberculosis complex (MTBC) genomes transmitted from patient to patient for tens of thousands years and possibly more. Recently discovered polymorphic loci, the CRISPR sequences, are indirect witnesses of the historical phage-bacteria struggle, and may be related to the time when the ancestor of today's tubercle bacilli were environmental bacteria, i.e. before becoming intracellular parasites. In this article, we present what are CRISPRs and try to summarize almost 20 years of research results obtained using the genetic diversity of the CRISPR loci in MTBC as a perspective for studying new models. We show that the study of the diversity of CRISPR sequences, thanks to «spoligotyping», has played a great role in our global understanding of the population structure of MTBC.

Mapping IS6110 in high-copy number Mycobacterium tuberculosis strains shows specific insertion points in the Beijing genotype

BACKGROUND

Mycobacterium tuberculosis Beijing strains are characterized by a large number of IS6110 copies, suggesting the potential implication of this element in the virulence and capacity for rapid dissemination characteristic of this family. This work studies the insetion points of IS6110 in high-copy clinical isolates specifically focusing on the Beijing genotype.

RESULTS

In the present work we mapped the insertion points of IS6110 in all the Beijing strains available in the literature and in the DNA sequence databases. We generated a representative primer collection of the IS6110 locations, which was used to analyse 61 high-copy clinical isolates. A total of 440 points of insertion were identified and analysis of their flanking regions determined the exact location, the direct repeats (DRs), the orientation and the distance to neighboring genes of each copy of IS6110. We identified specific points of insertion in Beijing strains that enabled us to obtain a dendrogram that groups the Beijing genotype.

CONCLUSIONS

This work presents a detailed analysis of locations of IS6110 in high-copy clinical isolates, showing points of insertion present with high frequency in the Beijing family and absent in other strains.

Epidemiologic consequences of microvariation in Mycobacterium tuberculosis

BACKGROUND

Evidence from genotype-phenotype studies suggests that genetic diversity in pathogens have clinically relevant manifestations that can impact outcome of infection and epidemiologic success. We studied 5 closely related Mycobacterium tuberculosis strains that collectively caused extensive disease (n = 862), particularly among US-born tuberculosis patients.

METHODS

Representative isolates were selected using population-based genotyping data from New York City and New Jersey. Growth and cytokine/chemokine response were measured in infected human monocytes. Survival was determined in aerosol-infected guinea pigs.

RESULTS

Multiple genotyping methods and phylogenetically informative synonymous single nucleotide polymorphisms showed that all strains were related by descent. In axenic culture, all strains grew similarly. However, infection of monocytes revealed 2 growth phenotypes, slower (doubling ∼55 hours) and faster (∼25 hours). The faster growing strains elicited more tumor necrosis factor α and interleukin 1β than the slower growing strains, even after heat killing, and caused accelerated death of infected guinea pigs (∼9 weeks vs 24 weeks) associated with increased lung inflammation/pathology. Epidemiologically, the faster growing strains were associated with human immunodeficiency virus and more limited in spread, possibly related to their inherent ability to induce a strong protective innate immune response in immune competent hosts.

CONCLUSIONS

Natural variation, with detectable phenotypic changes, among closely related clinical isolates of M. tuberculosis may alter epidemiologic patterns in human populations.

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.

Polymorphic exact tandem repeat A (PETRA): a newly defined lineage of mycobacterium tuberculosis in israel originating predominantly in Sub-Saharan Africa

As part of the Israel National Program for Prevention and Control of Tuberculosis, the molecular epidemiology of new tuberculosis cases is monitored. Prospective screening showed that about 20% of all new cases of culture-positive tuberculosis (43 of 222) in Israel in the year 2008 were caused by certain Mycobacterium tuberculosis strains of the central Asian (CAS) spoligotype lineage. The identity and similarity of these strains by mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing form a lineage we call PETRA for polymorphic at locus ETR A. The name PETRA was given to 79 strains we have found since the year 2000, because the largest number of strains with MIRU-VNTR profiles identical other than at locus A formed three groups, including 5 of 10 strains that had deleted the ETR A region from their genomes. No PETRA strain was found to be multiple drug resistant (resistant to both isoniazid and rifampin [rifampicin]). Most patients (75% [58 of 77 patients of known origin]) infected with PETRA were of sub-Saharan African origins. The genotypes associated with the 79 PETRA lineage strains presented in this paper suggest that the PETRA lineage is a large, major contributor to new tuberculosis cases in Israel.

Genome analysis of multi- and extensively-drug-resistant tuberculosis from KwaZulu-Natal, South Africa

The KZN strain family of Mycobacterium tuberculosis is a highly virulent strain endemic to the KwaZulu-Natal region of South Africa, which has recently experienced an outbreak of extensively-drug resistant tuberculosis. To investigate the causes and evolution of drug-resistance, we determined the DNA sequences of several clinical isolates--one drug-susceptible, one multi-drug resistant, and nine extensively drug-resistant--using whole-genome sequencing. Analysis of polymorphisms among the strains is consistent with the drug-susceptibility profiles, in that well-known mutations are observed that are correlated with resistance to isoniazid, rifampicin, kanamycin, ofloxacin, ethambutol, and pyrazinamide. However, the mutations responsible for rifampicin resistance in rpoB and pyrazinamide in pncA are in different nucleotide positions in the multi-drug-resistant and extensively drug-resistant strains, clearly showing that they acquired these mutations independently, and that the XDR strain could not have evolved directly from the MDR strain (though it could have arisen from another similar MDR strain). Sequencing of eight additional XDR strains from other areas of KwaZulu-Natal shows that they have identical drug resistant mutations to the first one sequenced, including the same polymorphisms at sites associated with drug resistance, supporting the theory that this represents a case of clonal expansion.

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.

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.

Molecular typing of Mycobacterium tuberculosis by using nine novel variable-number tandem repeats across the Beijing family and low-copy-number IS6110 isolates

Molecular epidemiological tools for genotyping clinical isolates of Mycobacterium tuberculosis have been developed and used to help track and contain transmission of tuberculosis. We identified 87 short sequence repeat loci within the genome of the M. tuberculosis H37Rv strain. Nine tandem repeats were found to be variable (variable-number tandem repeats [VNTRs]) in a set of 91 isolates. Fifty-seven of the isolates had only four IS6110 bands. The other 34 isolates were members of the Beijing strain family. The number of alleles of each these nine VNTRs was determined by examining each isolate. Six of the loci (Mtb-v1, -v4, -v10, -v15, -v18, and -v20) were able to differentiate the Beijing spoligotype identical isolates into seven distinct genotypes. Five of the loci (Mtb-v3, -v5, -v6, -v10, and -v15) were informative in discriminating the four-band IS6110 restriction fragment length polymorphism isolates from each other. The Nei's diversity values of each marker ranged from 0.02 to 0.59, with the number of alleles ranging from two to eight across the entire strain set. These nine loci provide a useful, discriminatory extension of VNTR typing methods for application to molecular epidemiologic studies of M. tuberculosis.

PCR-based methodology for detecting multidrug-resistant strains of Mycobacterium tuberculosis Beijing family circulating in Russia

The Beijing genotype of Mycobacterium tuberculosis has been identified in 40-50% of the clinical isolates studied in Russia during the last decade. This genotype has been reported to be associated with multiple drug resistance and possesses some significant pathogenic properties. Therefore, early identification of such strains is of extreme importance in the timely detection of drug resistance. The present study was performed on 354 strains isolated in Russia from 1996 to 2002 and previously characterised by IS 6110-restriction fragment length polymorphism (RFLP) typing and spoligotyping. These strains included 198 Beijing family strains and 156 strains of other genotypes (IS 6110-RFLP profiles). A subsequent polymerase chain reaction (PCR) analysis with IS 6110-derived outwardly oriented primers (IS 6110-PCR) easily discriminated the Beijing strains from non-Beijing strains. The multiplex allele-specific (MAS)-PCR assays were further used to detect mutations in katG315 and rpoB531, associated with resistance to isoniazid and rifampin, respectively. The katG315 and rpoB531 mutations were found to be more prevalent among Beijing (96.8% and 77.3%) than among non-Beijing strains (85.7% and 28%). Consequently, we propose a two-step methodology based on routine PCR and simple agarose gel electrophoresis in order to detect (i) a Beijing family strain using IS 6110-PCR, and, (ii) its possible resistance to the major anti-tuberculosis drugs using specific MAS-PCR assays.

Evolutionary relationships among strains of Mycobacterium tuberculosis with few copies of IS6110

Molecular typing of Mycobacterium tuberculosis by using IS6110 shows low discrimination when there are fewer than five copies of the insertion sequence. Using a collection of such isolates from a study of the epidemiology of tuberculosis in London, we have shown a substantial degree of congruence between IS6110 patterns and both spoligotype and PGRS type. This indicates that the IS6110 types mainly represent distinct families of strains rather than arising through the convergent insertion of IS6110 into favored positions. This is supported by identification of the genomic sites of the insertion of IS6110 in these strains. The combined data enable identification of the putative evolutionary relationships of these strains, comprising three lineages broadly associated with patients born in South Asia (India and Pakistan), Africa, and Europe, respectively. These lineages appear to be quite distinct from M. tuberculosis isolates with multiple copies of IS6110.

Phylogenetic reconstruction within Mycobacterium tuberculosis Beijing genotype in northwestern Russia

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

Determining the genomic locations of repetitive DNA sequences with a whole-genome microarray: IS6110 in Mycobacterium tuberculosis

The mycobacterial insertion sequence IS6110 has been exploited extensively as a clonal marker in molecular epidemiologic studies of tuberculosis. In addition, it has been hypothesized that this element is an important driving force behind genotypic variability that may have phenotypic consequences. We present here a novel, DNA microarray-based methodology, designated SiteMapping, that simultaneously maps the locations and orientations of multiple copies of IS6110 within the genome. To investigate the sensitivity, accuracy, and limitations of the technique, it was applied to eight Mycobacterium tuberculosis strains for which complete or partial IS6110 insertion site information had been determined previously. SiteMapping correctly located 64% (38 of 59) of the IS6110 copies predicted by restriction fragment length polymorphism analysis. The technique is highly specific; 97% of the predicted insertion sites were true insertions. Eight previously unknown insertions were identified and confirmed by PCR or sequencing. The performance could be improved by modifications in the experimental protocol and in the approach to data analysis. SiteMapping has general applicability and demonstrates an expansion in the applications of microarrays that complements conventional approaches in the study of genome architecture.

Novel IS6110 insertion sites in the direct repeat locus of Mycobacterium tuberculosis clinical strains from the St. Petersburg area of Russia and evolutionary and epidemiological considerations

A modification of spoligotyping with primers derived from the direct repeat (DR) and IS6110 sequences was used to identify IS6110 insertions in the DR locus of Mycobacterium tuberculosis clinical strains from the St. Petersburg area of Russia. Novel IS6110 insertions were identified: (i) in two epidemiologically unlinked Beijing family strains, an asymmetrical direct insertion in DR37; (ii) in a non-Beijing strain, an asymmetrical insertion in the opposite orientation in DR38; (iii) in another non-Beijing strain, a direct insertion in DR38 and one in the opposite orientation in DR14 (DR numbering is according to standard spoligotyping). Our results strengthen an observation that the DR locus structure is extremely conserved in the Beijing genotype. Asymmetrical insertions prevented detection of the adjacent spacer by standard spoligotyping. This, therefore, should be taken into consideration when similar spoligoprofiles that differ only in signals 37 and 38 are interpreted.