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

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

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

Several characteristics of Mycobacterium tuberculosis (e.g., conserved genome and low growth rate) have severely restricted the study of the microorganism. The discovery of IS6110 raised hopes of overcoming these obstacles. However, our knowledge of this IS element is relatively limited; even its two basic characteristics (transposition mechanism and target site selection) are far from well understood. In this study, IS6110 insertions in ipl loci (iplA and iplB) in two collections of clinical isolates of M. tuberculosis from different geographic locations, one from Scotland and the other from Thailand, were investigated. Five different IS6110 insertions in the loci were identified: ipl-4::IS6110, ipl-5::IS6110, ipl-11::IS6110, ipl-12::IS6110, and ipl-13::IS6110. An attempt to establish the phylogenetic relationship of the isolates containing these insertions was unsuccessful, suggesting that some of these insertions may have arisen from more than one event. This possibility is further supported by the observation that IS6110 copies existed in the same site but with different orientations in different isolates, and the insertion site of ipl-1::IS6110 harbored IS6110 copies in both iplA and iplB in different strains. All these suggest the independent occurrence of IS6110 insertions at the same sites of the genome of M. tuberculosis in different clinical isolates. The implications of this finding are discussed.     

IS6110-Mediated Deletions of Wild-Type Chromosomes of Mycobacterium tuberculosis

The ipl locus is a site for the preferential insertion of IS6110 and has been identified as an insertion sequence, IS1547, in its own right. Various deletions around the ipl locus of clinical isolates of Mycobacterium tuberculosis were identified, and these deletions ranged in length from several hundred base pairs up to several kilobase pairs. The most obvious feature shared by these deletions was the presence of an IS6110 copy at the deletion sites, which suggested two possible mechanisms for their occurrence, IS6110 transposition and homologous recombination. To clarify the mechanism, an investigation was conducted; the results suggest that although deletion transpositionally mediated by IS6110 was a possibility, homologous recombination was a more likely one. The implications of such chromosomal rearrangements for the evolution of M. tuberculosis, for IS6110-mediated mutagenesis, and for the development of genetic tools are discussed. The deletion of genomic DNA in isolates of M. tuberculosis has previously been noted at only a few sites. This study examined the deletional loss of genetic material at a new site and suggests that such losses may occur elsewhere too and may be more prevalent than was previously thought. Distinct from the study of laboratory-induced mutations, the detailed analysis of clinical isolates, in combination with knowledge of their evolutionary relationships to each other, gives us the opportunity to study mutational diversity in isolates that have survived in the human host and therefore offers a different perspective on the importance of particular genetic markers in pathogenesis.     

The control of copy number of IS6110 in Mycobacterium tuberculosis

Insertion sequence (IS) elements are bacterial genes that are able to transpose to different locations in the genome. These elements are often used in molecular epidemiology as genetic markers that track the spread of pathogens. Transposable elements have frequently been described as "selfish DNA" because they facilitate their own transposition, causing damage when they insert into coding regions, while contributing little if anything to the bacterial host. According to this hypothesis, the expansion of copy number of insertion sequences is opposed by negative selection against high copy numbers. From an alternative point of view, we might expect IS elements to intrinsically regulate transposition within cells, thereby limiting damage to their bacterial host. Here, we report evidence that the copy number of IS6110 in Mycobacterium tuberculosis is controlled by selection against the element. We first construct 12 different models of marker change resulting from a combination of possible transposition functions and selective regimes. We then compute the Akaike Information Criterion for each model to identify the models that best explain data consisting of serial isolates of M. tuberculosis genotyped with IS6110. We find that the best performing models all include selection against the accumulation of copies. Specifically, our analysis points to the interaction of separate copies of the element causing lethal effects. We discuss the implications of these findings for genome evolution and molecular epidemiology.     

IS6110 restriction fragment length polymorphism typing of Mycobacterium tuberculosis isolates from East Azerbaijan Province of Iran

To investigate the genetic variation among Mycobacterium tuberculosis isolates in the East Azerbaijan Province of Iran and to evaluate the level of and risk factors for recent transmission of tuberculosis (TB), we performed IS6110-based restriction fragment length polymorphism analysis of strains, isolated from 105 patients during the period of September 2002 to March 2003 in TB centers and university hospitals of the province. Among 105 isolates, 81 different IS6110 patterns were found, of which 70 were observed only once and 11 were shared by two to eight isolates. Ninety-six isolates (91.4%) were found to have more than five copies of IS6110 and together with high patterns polymorphism, shows that IS6110-RFLP typing could be useful for studying the epidemiology of TB in Azerbaijan. The minimum estimated rate of recent transmission was 23%, suggesting that the degree of recent transmission in East Azerbaijan Province is relatively low. Clustering was not associated with age, sex or site of infection of TB but drug-resistant isolates were less likely to be clustered than sensitive isolates (p < 0.05).     

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.     

Use of the insertion element IS6110 for DNA fingerprinting of Mycobacterium tuberculosis isolates presenting various profiles of drug susceptibility

Abstract IS 6100 is an insertion sequence of the IS3 family and it is present in multiple copies in the chromosome of Mycobacterium tuberculosis . Four to 15 copies are present in various strains of M. tuberculosis . In this study, the value of IS 6110 as an epidemiological marker of tuberculosis was examined. Unrelated clinical strains from Greek patients presented, in restriction fragment length polymorphism analysis, a high degree of polymorphism, whereas patterns of related clinical strains from familial outbreaks were identical. Since RFLP analysis with acetylaminofluorene labeled IS 6110 as the probe gave satisfactory results, it is suggested that this non-radioactive probe can be used in hospitals and health centres for the epidemiological survey of M. tuberculosis infections.     

Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome

Mycobacterial interspersed repetitive units (MIRUs) are 40-100 bp DNA elements often found as tandem repeats and dispersed in intergenic regions of the Mycobacterium tuberculosis complex genomes. The M. tuberculosis H37Rv chromosome contains 41 MIRU loci. After polymerase chain reaction (PCR) and sequence analyses of these loci in 31 M. tuberculosis complex strains, 12 of them were found to display variations in tandem repeat copy numbers and, in most cases, sequence variations between repeat units as well. These features are reminiscent of those of certain human variable minisatellites. Of the 12 variable loci, only one was found to vary among genealogically distant BCG substrains, suggesting that these interspersed bacterial minisatellite-like structures evolve slowly in mycobacterial populations.     

Biological properties of clinical isolates of Mycobacterium tuberculosis

The invasiveness of 2 grud sensitive and 8 holy-resistant M. tuberculosis clinical isolated was evaluated in experiments on BALB/c mice. Mycobacterial suspension was injected into the caudal vein of experimental animals. The results were evaluated by the degree of contamination of lungs and spleen of infected animals euthanized at different periods on time. The study revealed high variability in the degree of contamination of the organs of the animals infected with M. tuberculosis drug-resistant clinical isolates.     

IS6110, a Mycobacterium tuberculosis complex-specific insertion sequence, is also present in the genome of Mycobacterium smegmatis, suggestive of lateral gene transfer among mycobacterial species

IS6110 is an insertion element found exclusively within the members of the Mycobacterium tuberculosis complex (MTBC), and because of this exclusivity, it has become an important diagnostic tool in the identification of MTBC species. The restriction of IS6110 to the MTBC is hypothesized to arise from the inability of these bacteria to exchange DNA. We have identified an IS6110-related element in a strain of Mycobacterium smegmatis. The presence of IS6110 indicates that lateral gene transfer has occurred among mycobacterial species, suggesting that the mycobacterial gene pool is larger than previously suspected.     

Comparative genomic analysis of Mycobacterium tuberculosis clinical isolates

Background

Due to excessive antibiotic use, drug-resistant Mycobacterium tuberculosis has become a serious public health threat and a major obstacle to disease control in many countries. To better understand the evolution of drug-resistant M. tuberculosis strains, we performed whole genome sequencing for 7 M. tuberculosis clinical isolates with different antibiotic resistance profiles and conducted comparative genomic analysis of gene variations among them.

Results

We observed that all 7 M. tuberculosis clinical isolates with different levels of drug resistance harbored similar numbers of SNPs, ranging from 1409–1464. The numbers of insertion/deletions (Indels) identified in the 7 isolates were also similar, ranging from 56 to 101. A total of 39 types of mutations were identified in drug resistance-associated loci, including 14 previously reported ones and 25 newly identified ones. Sixteen of the identified large Indels spanned PE-PPE-PGRS genes, which represents a major source of antigenic variability. Aside from SNPs and Indels, a CRISPR locus with varied spacers was observed in all 7 clinical isolates, suggesting that they might play an important role in plasticity of the M. tuberculosis genome. The nucleotide diversity (Л value) and selection intensity (dN/dS value) of the whole genome sequences of the 7 isolates were similar. The dN/dS values were less than 1 for all 7 isolates (range from 0.608885 to 0.637365), supporting the notion that M. tuberculosis genomes undergo purifying selection. The Л values and dN/dS values were comparable between drug-susceptible and drug-resistant strains.

Conclusions

In this study, we show that clinical M. tuberculosis isolates exhibit distinct variations in terms of the distribution of SNP, Indels, CRISPR-cas locus, as well as the nucleotide diversity and selection intensity, but there are no generalizable differences between drug-susceptible and drug-resistant isolates on the genomic scale. Our study provides evidence strengthening the notion that the evolution of drug resistance among clinical M. tuberculosis isolates is clearly a complex and diversified process.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-469) contains supplementary material, which is available to authorized users.     

两株超级广泛耐药结核分枝杆菌全基因组比较分析

【目的】结合现有数据,通过对两株临床超级广泛耐药的结核分枝杆菌全基因组的测序和分析,发现其型别相关的突变位点,解释发生广泛耐药的基因组突变机制。【方法】利用Solexa第二代测序技术对两株广泛耐药结核分枝杆菌(FJ05194和GuangZ0019)进行全基因组测序分析。以H37Rv为参考序列得到两株广泛耐药菌株的单核苷酸多态性(SNPs),构建系统发育树鉴定菌株型别,判断突变位点中型别相关和非型别相关的SNPs。定位SNPs所在的基因组区域,对型别相关的突变基因进行KEGG通路的富集分析,对非型别相关的突变基因和间隔区判断是否与耐药相关。【结果】两株广泛耐药菌株分别属于Lineage2和Lineage4型别,两菌株在碱基替换方面存在差异性,Lineage2型别相关的基因功能富集于ABC转运蛋白和核苷酸切除修复的通路。耐药方面,发现了已知的耐药相关基因的突变(rpoB、katG、rpsl、gyrA、gyrB、embB和ethA等),但卷曲霉素和卡那霉素相关的rrs、tlyA和eis启动子区域未发生突变,不足以解释其耐药性的产生。与最新报道的候选耐药基因比较,发现了卷曲霉素和卡那霉素相关的突变(Rv1393c、Rv0265c和narX等)和外排泵相关的pstB、Rv2333c和Rv2687c突变。【结论】结核分枝杆菌Lineage2型别相关的SNPs中含有影响结核分枝杆菌突变率和耐药性的突变。对于两株超级广泛耐药的结核菌,已知的激活药物或药靶相关的单耐药基因突变集合不能完全解释其广泛耐药性,还涉及新候选结核耐药基因、外排泵和补偿等其他潜在机制的相关基因突变。     

Drug resistance and IS6110‐RFLP patterns of Mycobacterium tuberculosis in patients with recurrent tuberculosis in northern Thailand

The emergence of drug resistant Mycobacterium tuberculosis has become a global threat to tuberculosis (TB) prevention and control efforts. This study aimed to determine the drug resistance profiles and DNA fingerprints of M. tuberculosis strains isolated from patients with relapsed or retreatment pulmonary TB in Chiang Rai province in northern Thailand. Significant differences in multidrug resistance (MDR) (P = 0.025) and resistance to isoniazid (P = 0.025) and rifampin (P = 0.046) between first and second registrations of patients with retreatment TB were found. However, there were no significant differences in resistance to any drugs in patients with relapsed TB. The rate of MDR‐TB strains was 12.2% among new patients at first registration, 22.5% among patients with recurrence who had previously undergone treatment at second registration and 12.5% at third registration. Two retreatment patients whose initial treatment had failed had developed MDR‐TB with resistance to all TB drugs tested, including rifampin, isoniazid, streptomycin and ethambutol. IS6110‐RFLP analysis revealed that 66.7% (10/15 isolates) of MDR‐TB belonged to the Beijing family. In most cases, IS6110‐RFLP patterns of isolates from the same patients were identical in relapse and retreatment groups. However, some pairs of isolates from retreatment patients after treatment failure had non‐identical IS6110‐RFLP patterns. These results suggest that, after failure and default treatment, patients with retreatment tuberculosis have a significantly greater risk of MDR‐TB, isoniazid and rifampin resistance than do other patients.     

IS6110 fingerprinting of sensitive and resistant strains (1991-1992) of Mycobacterium tuberculosis in Colombia

The standardized method to study the polymorphism of IS 6110 was used to characterize 53 isolates of Mycobacterium tuberculosis obtained during 1991-1992 from 14 regions in Colombia. In Valle region cluster rate was 25% (4/16). The mean number of IS6110 band was 10 +/- 3. Similarity between strains was of 60% in 81% of strains and this tended to be correlated with geographic origin. For the first time M. tuberculosis without IS6110 bands in restriction fragment length polymorphism analysis was found in Colombia. Additional studies are necessaries in order to best characterize the situation in relation to human immunodeficiency virus epidemic and recent changes in tuberculosis control program.     

Stimulation of transposition of the Mycobacterium tuberculosis insertion sequence IS6110 by exposure to a microaerobic environment

The Mycobacterium tuberculosis-specific insertion sequence IS6110/986 has been widely used as a probe because of the multiple polymorphism observed among different strains. To investigate transposition of IS6110, a series of artificially constructed composite transposons containing IS6110 and a kanamycin resistance marker were constructed. The composite transposons were inserted into a conditionally replicating, thermosensitive, Escherichia coli-mycobacterial shuttle vector and introduced into M. smegmatis mc2155. Lawns of transformants were grown at the permissive temperature on kanamycin-supplemented agar and subsequently prevented from further growth by shifting to the non-permissive temperature. Under normal atmospheric conditions, kanamycin-resistant papillae appeared after only about 5-6 weeks of incubation. However, these events were not associated with transposon mobilization. In contrast, lawns that were exposed to a 48 h microaerobic shock generated kanamycin-resistant papillae after only 6-14 days. These events were generated by conservative transposition of the IS6110 composite transposon into the M. smegmatis chromosome, with loss of the shuttle vector. In common with other IS3 family elements, transposition of IS6110 is thought to be controlled by translational frameshifting. However, we were unable to detect any significant frameshifting within the putative frameshifting site of IS6110, and the level of frameshifting was not affected by microaerobic incubation. The finding that transposition of IS6110 is stimulated by incubation at reduced oxygen tensions may be relevant to transposition of IS6110 in M. tuberculosis harboured within TB lesions.     

Genetic diversity of Iranian clinical isolates of Mycobacterium tuberculosis

In the current study we aimed to execute a rather less complicated molecular tying method, i.e. the random amplification of polymorphic DNA (RAPD) analysis to find the heterogeneity of Iranian strains of Mycobacterium tuberculosis. The isolates comprised a total of 96 strains of M. tuberculosis collected from clinical specimens of patients in Isfahan and Tehran. The isolates were assigned to the species M. tuberculosis by the key conventional and molecular methods. They were then subjected to RAPD analysis by four arbitrary primers, namely, the primers 27F, 1525R, MS- GF and INS-2. They were then evaluated for the number and intensity of the band patterns. The RAPD profiles of the Iranian isolates showed a degree of heterogeneity which varied based on the primer used. However, analysis of the isolates by primer INS-2 revealed the highest degree of diversity yielding 31 distinguishable RAPD types. RAPD analysis provides a rapid and easy means of identifying heterogeneity among the M. tuberculosis isolates. This typing system might be considered a valuable alternative molecular typing for countries with limited resources provided that the reproducibility and reliability of the method is carefully assured.     

Analysis of sequence diversity among IS6110 sequence of Mycobacterium tuberculosis: possible implications for PCR based detection

The IS6110 belongs to the family of insertion sequences (IS) of the IS3 category. This insertion sequence was reported to be specific for Mycobacterium tuberculosis complex and hence is extensively exploited for laboratory detection of the agent of tuberculosis and for epidemiological investigations based on polymerase chain reaction. IS6110 is 1361-bp long and within this sequence different regions have been utilized as targets in the identification of M. tuberculosis by PCR. However, the results are not always consistent, specific and sensitive. In recent years, a few clinical investigations raised concerns over IS6110 specificity and sensitivity in the diagnosis of tuberculosis due to false-positive (homology with other target DNA besides M. tuberculosis) or false negative (due to absence of copies of IS6110) results with IS6110 specific primers. To unravel the variations in IS6110 sequences, an insilico analysis of IS6110 sequence of different strains of M. tuberculosis was carried out. Our results of comparative analysis of IS6110 insertion sequences of M. tuberculosis complex suggests that, IS6110 insertion sequences harbored variations in its sequence, which is evident from the phylogenetic analysis. Importantly, IS6110 sequence has divergence within the copies of same strain and formed different clusters. A list of IS6110 specific primers used in various clinical investigation of tuberculosis was obtained from the literature and their performance scrutinized. Our study emphasizes the need to develop PCR assays (multiplex format) targeting more than one region of the genome of M. tuberculosis.