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

Genomic Characterization of IS6110 Insertions in Mycobacterium orygis

Mycobacterium orygis, a subspecies of the Mycobacterium tuberculosis complex (MTBC), has emerged as a significant concern in the context of One Health, with implications for zoonosis or zooanthroponosis or both. MTBC strains are characterized by the unique insertion element IS6110, which is widely used as a diagnostic marker. IS6110 transposition drives genetic modifications in MTBC, imparting genome plasticity and profound biological consequences. While IS6110 insertions are customarily found in the MTBC genomes, the evolutionary trajectory of strains seems to correlate with the number of IS6110 copies, indicating enhanced adaptability with increasing copy numbers. Here, we present a comprehensive analysis of IS6110 insertions in the M. orygis genome, utilizing ISMapper, and elucidate their genetic consequences in promoting successful host adaptation. Our study encompasses a panel of 67 paired-end reads, comprising 11 isolates from our laboratory and 56 sequences downloaded from public databases. Among these sequences, 91% exhibited high-copy, 4.5% low-copy, and 4.5% lacked IS6110 insertions. We identified 255 insertion loci, including 141 intragenic and 114 intergenic insertions. Most of these loci were either unique or shared among a limited number of isolates, potentially influencing strain behavior. Furthermore, we conducted gene ontology and pathway analysis, using eggNOG-mapper 5.0, on the protein sequences disrupted by IS6110 insertions, revealing 63 genes involved in diverse functions of Gene Ontology and 45 genes participating in various KEGG pathways. Our findings offer novel insights into IS6110 insertions, their preferential insertion regions, and their impact on metabolic processes and pathways, providing valuable knowledge on the genetic changes underpinning IS6110 transposition in M. orygis.

Transcriptomic profile of the most successful Mycobacterium tuberculosis strain in Aragon, the MtZ strain, during exponential and stationary growth phases

IMPORTANCE

Aragon Community suffered, during the first years of the beginning of this century, a large outbreak caused by the MtZ strain, producing more than 240 cases to date. MtZ strain and the outbreak have been previously studied from an epidemiological and molecular point of view. In this work, we analyzed the transcriptomic profile of the strain for better understanding of its success among our population. We have discovered that MtZ has some upregulated virulence pathways, such as the ESX-1 system, the cholesterol degradation pathway or the peptidoglycan biosynthesis. Interestingly, MtZ has downregulated the uptake of iron. Another special feature of MtZ strain is the interruption of desA3 gene, essential for producing oleic acid. Although the strain takes a long time to grow in the initial culture media, eventually it is able to reach normal optical densities, suggestive of the presence of another route for obtaining oleic acid in Mycobacterium tuberculosis.

IS6110 Copy Number in Multi-Host Mycobacterium bovis Strains Circulating in Bovine Tuberculosis Endemic French Regions

IS6110 is an insertion sequence found in the Mycobacterium tuberculosis complex, to which Mycobacterium bovis belongs, which can play a role in genome plasticity and in bacterial evolution. In this study, the abundance and location of IS6110 on M. bovis genomic data of French animal field strains were studied. A first analysis was performed on a panel of 81 strains that reflect the national M. bovis population’s genetic diversity. The results show that more than one-third of them are IS6110 multicopy and that 10% have IS6110 in a high copy number (more than 6 copies). Multicopy strains are those circulating in the regions where prevalence was above the national average. Further study of 93 such strains, with an IS6110 copy number of 10-12, showed stability of IS6110 copy number and genome location over time and between host species. The correlation between M. bovis multicopy strains and high bovine tuberculosis (bTB) prevalence leads us to consider whether their epidemiological success could be partly due to genetic changes originated by IS6110 transposition.

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

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

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

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

Proline-Glutamate/Proline-Proline-Glutamate (PE/PPE) proteins of Mycobacterium tuberculosis: The multifaceted immune-modulators

The PE/PPE proteins encoded by seven percent (7%) of Mycobacterium tuberculosis (Mtb) genome are the chief constituents to pathogen's virulence reservoir. The fact that these genes have evolved along ESX secretory system in pathogenic Mtb strains make their investigation very intriguing. There is lot of speculation about the prominent role of these proteins at host pathogen interface and in disease pathogenesis. Nevertheless, the exact function of PE/PPE proteins still remains a mystery which calls for further research targeting these proteins. This article is an effort to document all the facts known so far with regard to these unique proteins which involves their origin, evolution, transcriptional control, and most important their role as host immune-modulators. Our understanding strongly points towards the versatile nature of these PE/PPE proteins as Mtb's host immune sensors and as decisive factors in shaping the outcome of infection. Further investigation on these proteins will surely pave way for newer and effective vaccines and therapeutics to control Tuberculosis (TB).

Aptamer-Based Diagnostic Systems for the Rapid Screening of TB at the Point-of-Care

The transmission of Tuberculosis (TB) is very rapid and the burden it places on health care systems is felt globally. The effective management and prevention of this disease requires that it is detected early. Current TB diagnostic approaches, such as the culture, sputum smear, skin tuberculin, and molecular tests are time-consuming, and some are unaffordable for low-income countries. Rapid tests for disease biomarker detection are mostly based on immunological assays that use antibodies which are costly to produce, have low sensitivity and stability. Aptamers can replace antibodies in these diagnostic tests for the development of new rapid tests that are more cost effective; more stable at high temperatures and therefore have a better shelf life; do not have batch-to-batch variations, and thus more consistently bind to a specific target with similar or higher specificity and selectivity and are therefore more reliable. Advancements in TB research, in particular the application of proteomics to identify TB specific biomarkers, led to the identification of a number of biomarker proteins, that can be used to develop aptamer-based diagnostic assays able to screen individuals at the point-of-care (POC) more efficiently in resource-limited settings.

Rapidly Correcting Frameshift Mutations in the Mycobacterium tuberculosis orn Gene Produce Reversible Ethambutol Resistance and Small-Colony-Variant Morphology

We have identified a previously unknown mechanism of reversible high-level ethambutol (EMB) resistance in Mycobacterium tuberculosis that is caused by a reversible frameshift mutation in the M. tuberculosis orn gene. A frameshift mutation in orn produces the small-colony-variant (SCV) phenotype, but this mutation does not change the MICs of any drug for wild-type M. tuberculosis However, the same orn mutation in a low-level EMB-resistant double embB-aftA mutant (MIC = 8 μg/ml) produces an SCV with an EMB MIC of 32 μg/ml. Reversible resistance is indistinguishable from a drug-persistent phenotype, because further culture of these orn-embB-aftA SCV mutants results in rapid reversion of the orn frameshifts, reestablishing the correct orn open reading frame, returning the culture to normal colony size, and reversing the EMB MIC back to that (8 μg/ml) of the parental strain. Transcriptomic analysis of orn-embB-aftA mutants compared to wild-type M. tuberculosis identified a 27-fold relative increase in the expression of embC, which is a cellular target for EMB. Expression of embC in orn-embB-aftA mutants was also increased 5-fold compared to that in the parental embB-aftA mutant, whereas large-colony orn frameshift revertants of the orn-embB-aftA mutant had levels of embC expression similar to that of the parental embB-aftA strain. Reversible frameshift mutants may contribute to a reversible form of microbiological drug resistance in human tuberculosis.

Phase variation in Mycobacterium tuberculosis glpK produces transiently heritable drug tolerance

The length and complexity of tuberculosis (TB) therapy, as well as the propensity of Mycobacterium tuberculosis to develop drug resistance, are major barriers to global TB control efforts. M. tuberculosis is known to have the ability to enter into a drug-tolerant state, which may explain many of these impediments to TB treatment. We have identified a mechanism of genetically encoded but rapidly reversible drug tolerance in M. tuberculosis caused by transient frameshift mutations in a homopolymeric tract (HT) of 7 cytosines (7C) in the glpK gene. Inactivating frameshift mutations associated with the 7C HT in glpK produce small colonies that exhibit heritable multidrug increases in minimal inhibitory concentrations and decreases in drug-dependent killing; however, reversion back to a fully drug-susceptible large-colony phenotype occurs rapidly through the introduction of additional insertions or deletions in the same glpK HT region. These reversible frameshift mutations in the 7C HT of M. tuberculosis glpK occur in clinical isolates, accumulate in M. tuberculosis-infected mice with further accumulation during drug treatment, and exhibit a reversible transcriptional profile including induction of dosR and sigH and repression of kstR regulons, similar to that observed in other in vitro models of M. tuberculosis tolerance. These results suggest that GlpK phase variation may contribute to drug tolerance, treatment failure, and relapse in human TB. Drugs effective against phase-variant M. tuberculosis may hasten TB treatment and improve cure rates.

A Mycobacterium tuberculosis surface protein recruits ubiquitin to trigger host xenophagy

Ubiquitin-mediated xenophagy, a type of selective autophagy, plays crucial roles in host defense against intracellular pathogens including Mycobacterium tuberculosis (Mtb). However, the exact mechanism by which host ubiquitin targets invaded microbes to trigger xenophagy remains obscure. Here we show that ubiquitin could recognize Mtb surface protein Rv1468c, a previously unidentified ubiquitin-binding protein containing a eukaryotic-like ubiquitin-associated (UBA) domain. The UBA-mediated direct binding of ubiquitin to, but not E3 ubiquitin ligases-mediated ubiquitination of, Rv1468c recruits autophagy receptor p62 to deliver mycobacteria into LC3-associated autophagosomes. Disruption of Rv1468c-ubiquitin interaction attenuates xenophagic clearance of Mtb in macrophages, and increases bacterial loads in mice with elevated inflammatory responses. Together, our findings reveal a unique mechanism of host xenophagy triggered by direct binding of ubiquitin to the pathogen surface protein, and indicate a diplomatic strategy adopted by Mtb to benefit its persistent intracellular infection through controlling intracellular bacterial loads and restricting host inflammatory responses.

Ubiquitin (Ub)-mediated xenophagy is important in defense against Mycobacterium tuberculosis (Mtb). Here, Chai et al. describe autophagy triggering by Ub binding to the Mtb surface protein Rv1468c, and show that its deletion leads to increased bacterial loads and hyperinflammatory responses in mice.

Primary transcriptome analysis reveals importance of IS elements for the shaping of the transcriptional landscape of Bordetella pertussis

Bordetella pertussis is the causative agent of whooping cough, a respiratory disease still considered as a major public health threat and for which recent re-emergence has been observed. Constant reshuffling of Bordetella pertussis genome organization was observed during evolution. These rearrangements are essentially mediated by Insertion Sequences (IS), a mobile genetic elements present in more than 230 copies in the genome, which are supposed to be one of the driving forces enabling the pathogen to escape from vaccine-induced immunity. Here we use high-throughput sequencing approaches (RNA-seq and differential RNA-seq), to decipher Bordetella pertussis transcriptome characteristics and to evaluate the impact of IS elements on transcriptome architecture. Transcriptional organization was determined by identification of transcription start sites and revealed also a large variety of non-coding RNAs including sRNAs, leaderless mRNAs or long 3' and 5'UTR including seven riboswitches. Unusual topological organizations, such as overlapping 5'- or 3'-extremities between oppositely orientated mRNA were also unveiled. The pivotal role of IS elements in the transcriptome architecture and their effect on the transcription of neighboring genes was examined. This effect is mediated by the introduction of IS harbored promoters or by emergence of hybrid promoters. This study revealed that in addition to their impact on genome rearrangements, most of the IS also impact on the expression of their flanking genes. Furthermore, the transcripts produced by IS are strain-specific due to the strain to strain variation in IS copy number and genomic context.

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

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

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

Impact of Genetic Diversity on the Biology of Mycobacterium tuberculosis Complex Strains

Tuberculosis (TB) remains the most deadly bacterial infectious disease worldwide. Its treatment and control are threatened by increasing numbers of multidrug-resistant (MDR) or nearly untreatable extensively drug-resistant (XDR) strains. New concepts are therefore urgently needed to understand the factors driving the TB epidemics and the spread of different strain populations, especially in association with drug resistance. Classical genotyping and, more recently, whole-genome sequencing (WGS) revealed that the world population of tubercle bacilli is more diverse than previously thought. Several major phylogenetic lineages can be distinguished, which are associated with their sympatric host population. Distinct clonal (sub)populations can even coexist within infected patients. WGS is now used as the ultimate approach for differentiating clinical isolates and for linking phenotypic to genomic variation from lineage to strain levels. Multiple lines of evidence indicate that the genetic diversity of TB strains translates into pathobiological consequences, and key molecular mechanisms probably involved in differential pathoadaptation of some main lineages have recently been identified. Evidence also accumulates on molecular mechanisms putatively fostering the emergence and rapid expansion of particular MDR and XDR strain groups in some world regions. However, further integrative studies will be needed for complete elucidation of the mechanisms that allow the pathogen to infect its host, acquire multidrug resistance, and transmit so efficiently. Such knowledge will be key for the development of the most effective new diagnostics, drugs, and vaccination strategies.

The impact of insertion sequences on bacterial genome plasticity and adaptability

Transposable elements (TE), small mobile genetic elements unable to exist independently of the host genome, were initially believed to be exclusively deleterious genomic parasites. However, it is now clear that they play an important role as bacterial mutagenic agents, enabling the host to adapt to new environmental challenges and to colonize new niches. This review focuses on the impact of insertion sequences (IS), arguably the smallest TE, on bacterial genome plasticity and concomitant adaptability of phenotypic traits, including resistance to antibacterial agents, virulence, pathogenicity and catabolism. The direct consequence of IS transposition is the insertion of one DNA sequence into another. This event can result in gene inactivation as well as in modulation of neighbouring gene expression. The latter is usually mediated by de-repression or by the introduction of a complete or partial promoter located within the element. Furthermore, transcription and transposition of IS are affected by host factors and in some cases by environmental signals offering the host an adaptive strategy and promoting genetic variability to withstand the environmental challenges.

Biological and Epidemiological Consequences of MTBC Diversity

Tuberculosis is caused by different groups of bacteria belonging to the Mycobacterium tuberculosis complex (MTBC). The combined action of human factors, environmental conditions and bacterial virulence determine the extent and form of human disease. MTBC virulence is a composite of different clinical phenotypes such as transmission rate and disease severity among others. Clinical phenotypes are also influenced by cellular and immunological phenotypes. MTBC phenotypes are determined by the genotype, therefore finding genotypes responsible for clinical phenotypes would allow discovering MTBC virulence factors. Different MTBC strains display different cellular and clinical phenotypes. Strains from Lineage 5 and Lineage 6 are metabolically different, grow slower, and are less virulent. Also, at least certain groups of Lineage 2 and Lineage 4 strains are more virulent in terms of disease severity and human-to-human transmission. Because phenotypic differences are ultimately caused by genotypic differences, different genomic loci have been related to various cellular and clinical phenotypes. However, defining the impact of specific bacterial genomic loci on virulence when other bacterial determinants, human and environmental factors are also impacting the phenotype would contribute to a better knowledge of tuberculosis virulence and ultimately benefit tuberculosis control.

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.

Geographic Differences in the Contribution of ubiA Mutations to High-Level Ethambutol Resistance in Mycobacterium tuberculosis

Ethambutol (EMB) resistance can evolve through a multistep process, and mutations in the ubiA (Rv3806c) gene appear to be responsible for high-level EMB resistance in Mycobacterium tuberculosis We evaluated the prevalence of ubiA and embB (Rv3795) mutations in EMB-resistant strains originating from Africa and South Korea. No differences in embB mutation frequencies were observed between strains from both origins. However, ubiA mutations were present in 45.5% ± 6.5% of the African EMB-resistant isolates but in only 9.5% ± 1.5% of the South Korean EMB-resistant isolates. The ubiA mutations associated with EMB resistance were localized to regions encoding the transmembrane domains of the protein, whereas the embB mutations were localized to regions encoding the extramembrane domains. Larger studies are needed to investigate the causes of increased ubiA mutations as a pathway to high-level EMB resistance in African countries, such as extended EMB usage during tuberculosis treatment.

Evaluation of multiplex PCR using MPB64 and IS6110 primers for rapid diagnosis of tuberculous meningitis

Tuberculous meningitis (TBM) is one of those most serious manifestations of extra-pulmonary tuberculosis and prompt diagnosis and treatment is required for better clinical outcome. It is difficult to diagnose due to lack of rapid, sensitive, and specific tests. Newer methods, which are easy and reliable, are required to diagnose TBM at an early stage. Thus our aim was to evaluate the Multiplex polymerase chain reaction (PCR) technique, using primers directed against the insertion sequence IS6110 and MPB64 gene for the detection of Mycobacterium tuberculosis in Cerebrospinal fluid (CSF), for rapid diagnosis of TBM patients. 102 CSF samples were analyzed from patients suspected with TBM along with a control group of 10 patients having other neurological disorders. CSF sediments were analyzed individually for M. tuberculosis DNA by Multiplex PCR using two set of primers targeting insertion sequence IS6110 and gene MBp64, which is very specific for MTBC. Out of 37 patients diagnosed with TBM clinically, MPB64 PCR was positive in 22, IS6110 PCR was positive in 28, both PCR using Multiplex were positive in 34 and Microscopy was positive in one. Thus Sensitivity of MPB64 PCR, IS6110 PCR, Multiplex PCR and Microscopy were found to be 62.3%, 75.4%, 91.8% and 2.7% respectively. In non TBM group PCR was negative in all cases hence, the specificity was 100%. Multiplex PCR system using primers targeting IS6110 and MPB64, for the detection of M. tuberculosis DNA in CSF samples, has high sensitivity than any one of them alone, and could be used for the early detection of TBM in CSF samples.

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

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

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

Background

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

Methodology/Principal Findings

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

Conclusions/Significance

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

In Vivo IS6110 Profile Changes in a Mycobacterium tuberculosis Strain as Determined by Tracking over 14 Years

Transposition and homologous recombination of IS6110 appear in Mycobacterium tuberculosis along in vivo sequential infections. These events were checked in different clones of a successful strain, M. tuberculosis Zaragoza, with the focus on a variant in which integration of a copy of IS6110 in the origin of replication (oriC) region occurred.

Phylogeny to function: PE/PPE protein evolution and impact on Mycobacterium tuberculosis pathogenicity

The pe/ppe genes represent one of the most intriguing aspects of the Mycobacterium tuberculosis genome. These genes are especially abundant in pathogenic mycobacteria, with more than 160 members in M. tuberculosis. Despite being discovered over 15 years ago, their function remains unclear, although various lines of evidence implicate selected family members in mycobacterial virulence. In this review, we use PE/PPE phylogeny as a framework within which we examine the diversity and putative functions of these proteins. We report on the evolution and diversity of the respective gene families, as well as the implications thereof for function and host immune recognition. We summarize recent findings on pe/ppe gene regulation, also placing this in the context of PE/PPE phylogeny. We collate data from several large proteomics datasets, providing an overview of PE/PPE localization, and discuss the implications this may have for host responses. Assessment of the current knowledge of PE/PPE diversity suggests that these proteins are not variable antigens as has been so widely speculated; however, they do clearly play important roles in virulence. Viewing the growing body of pe/ppe literature through the lens of phylogeny reveals trends in features and function that may be associated with the evolution of mycobacterial pathogenicity.

Consequences of genomic diversity in Mycobacterium tuberculosis

The causative agent of human tuberculosis, Mycobacterium tuberculosis complex (MTBC), comprises seven phylogenetically distinct lineages associated with different geographical regions. Here we review the latest findings on the nature and amount of genomic diversity within and between MTBC lineages. We then review recent evidence for the effect of this genomic diversity on mycobacterial phenotypes measured experimentally and in clinical settings. We conclude that overall, the most geographically widespread Lineage 2 (includes Beijing) and Lineage 4 (also known as Euro-American) are more virulent than other lineages that are more geographically restricted. This increased virulence is associated with delayed or reduced pro-inflammatory host immune responses, greater severity of disease, and enhanced transmission. Future work should focus on the interaction between MTBC and human genetic diversity, as well as on the environmental factors that modulate these interactions.

Unraveling Mycobacterium tuberculosis genomic diversity and evolution in Lisbon, Portugal, a highly drug resistant setting

Background

Multidrug- (MDR) and extensively drug resistant (XDR) tuberculosis (TB) presents a challenge to disease control and elimination goals. In Lisbon, Portugal, specific and successful XDR-TB strains have been found in circulation for almost two decades.

Results

In the present study we have genotyped and sequenced the genomes of 56 Mycobacterium tuberculosis isolates recovered mostly from Lisbon. The genotyping data revealed three major clusters associated with MDR-TB, two of which are associated with XDR-TB. Whilst the genomic data contributed to elucidate the phylogenetic positioning of circulating MDR-TB strains, showing a high predominance of a single SNP cluster group 5. Furthermore, a genome-wide phylogeny analysis from these strains, together with 19 publicly available genomes of Mycobacterium tuberculosis clinical isolates, revealed two major clades responsible for M/XDR-TB in the region: Lisboa3 and Q1 (LAM).

The data presented by this study yielded insights on microevolution and identification of novel compensatory mutations associated with rifampicin resistance in rpoB and rpoC. The screening for other structural variations revealed putative clade-defining variants. One deletion in PPE41, found among Lisboa3 isolates, is proposed to contribute to immune evasion and as a selective advantage. Insertion sequence (IS) mapping has also demonstrated the role of IS6110 as a major driver in mycobacterial evolution by affecting gene integrity and regulation.

Conclusions

Globally, this study contributes with novel genome-wide phylogenetic data and has led to the identification of new genomic variants that support the notion of a growing genomic diversity facing both setting and host adaptation.

Electronic supplementary material

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

Evolutionary history of tuberculosis shaped by conserved mutations in the PhoPR virulence regulator

Although the bovine tuberculosis (TB) agent, Mycobacterium bovis, may infect humans and cause disease, long-term epidemiological data indicate that humans represent a spill-over host in which infection with M. bovis is not self-maintaining. Indeed, human-to-human transmission of M. bovis strains and other members of the animal lineage of the tubercle bacilli is very rare. Here, we report on three mutations affecting the two-component virulence regulation system PhoP/PhoR (PhoPR) in M. bovis and in the closely linked Mycobacterium africanum lineage 6 (L6) that likely account for this discrepancy. Genetic transfer of these mutations into the human TB agent, Mycobacterium tuberculosis, resulted in down-regulation of the PhoP regulon, with loss of biologically active lipids, reduced secretion of the 6-kDa early antigenic target (ESAT-6), and lower virulence. Remarkably, the deleterious effects of the phoPR mutations were partly compensated by a deletion, specific to the animal-adapted and M. africanum L6 lineages, that restores ESAT-6 secretion by a PhoPR-independent mechanism. Similarly, we also observed that insertion of an IS6110 element upstream of the phoPR locus may completely revert the phoPR-bovis-associated fitness loss, which is the case for an exceptional M. bovis human outbreak strain from Spain. Our findings ultimately explain the long-term epidemiological data, suggesting that M. bovis and related phoPR-mutated strains pose a lower risk for progression to overt human TB, with major impact on the evolutionary history of TB.

Physiological impact of transposable elements encoding DDE transposases in the environmental adaptation of Streptococcus agalactiae

We have referenced and described Streptococcus agalactiae transposable elements encoding DDE transposases. These elements belonged to nine families of insertion sequences (ISs) and to a family of conjugative transposons (TnGBSs). An overview of the physiological impact of the insertion of all these elements is provided. DDE-transposable elements affect S. agalactiae in a number of aspects of its capability to adapt to various environments and modulate the expression of several virulence genes, the scpB–lmB genomic region and the genes involved in capsule expression and haemolysin transport being the targets of several different mobile elements. The referenced mobile elements modify S. agalactiae behaviour by transferring new gene(s) to its genome, by modifying the expression of neighbouring genes at the integration site or by promoting genomic rearrangements. Transposition of some of these elements occurs in vivo, suggesting that by dynamically regulating some adaptation and/or virulence genes, they improve the ability of S. agalactiae to reach different niches within its host and ensure the ‘success’ of the infectious process.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

Evolution of high-level ethambutol-resistant tuberculosis through interacting mutations in decaprenylphosphoryl-β-D-arabinose biosynthetic and utilization pathway genes

To study the evolution of drug resistance, we genetically and biochemically characterized Mycobacterium tuberculosis strains selected in vitro for ethambutol resistance. Mutations in decaprenylphosphoryl-β-D-arabinose (DPA) biosynthetic and utilization pathway genes Rv3806c, Rv3792, embB and embC accumulated to produce a wide range of ethambutol minimal inhibitory concentrations (MICs) that depended on mutation type and number. Rv3806c mutations increased DPA synthesis, causing MICs to double from 2 to 4 μg/ml in a wild-type background and to increase from 16 to 32 μg/ml in an embB codon 306 mutant background. Synonymous mutations in Rv3792 increased the expression of downstream embC, an ethambutol target, resulting in MICs of 8 μg/ml. Multistep selection was required for high-level resistance. Mutations in embC or very high embC expression were observed at the highest resistance level. In clinical isolates, Rv3806c mutations were associated with high-level resistance and had multiplicative effects with embB mutations on MICs. Ethambutol resistance is acquired through the acquisition of mutations that interact in complex ways to produce a range of MICs, from those falling below breakpoint values to ones representing high-level resistance.

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

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

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.

Genetic heterogeneity revealed by sequence analysis of Mycobacterium tuberculosis isolates from extra-pulmonary tuberculosis patients

BACKGROUND

Tuberculosis remains a major public health problem. Clinical tuberculosis manifests often as pulmonary and occasionally as extra-pulmonary tuberculosis. The emergence of drug resistant tubercle bacilli and its association with HIV is a formidable challenge to curb the spread of tuberculosis. There have been concerted efforts by whole genome sequencing and bioinformatics analysis to identify genomic patterns and to establish a relationship between the genotype of the organism and clinical manifestation of tuberculosis. Extra-pulmonary TB constitutes 15-20 percent of the total clinical cases of tuberculosis reported among immunocompetent patients, whereas among HIV patients the incidence is more than 50 percent. Genomic analysis of M. tuberculosis isolates from extra pulmonary patients has not been explored.

RESULTS

The genomic DNA of 5 extra-pulmonary clinical isolates of M. tuberculosis derived from cerebrospinal fluid, lymph node fine needle aspirates (FNAC) / biopsies, were sequenced. Next generation sequencing approach (NGS) was employed to identify Single Nucleotide Variations (SNVs) and computational methods used to predict their consequence on functional genes. Analysis of distribution of SNVs led to the finding that there are mixed genotypes in patient isolates and that many SNVs are likely to influence either gene function or their expression. Phylogenetic relationship between the isolates correlated with the origin of the isolates. In addition, insertion sites of IS elements were identified and their distribution revealed a variation in number and position of the element in the 5 extra-pulmonary isolates compared to the reference M. tuberculosis H37Rv strain.

CONCLUSIONS

The results suggest that NGS sequencing is able to identify small variations in genomes of M. tuberculosis isolates including changes in IS element insertion sites. Moreover, variations in isolates of M. tuberculosis from non-pulmonary sites were documented. The analysis of our results indicates genomic heterogeneity in the clinical isolates.

The complex architecture of mycobacterial promoters

The genus Mycobacterium includes a variety of species with differing phenotypic properties, including growth rate, pathogenicity and environment- and host-specificity. Although many mycobacterial species have been extensively studied and their genomes sequenced, the reasons for phenotypic variation between closely related species remain unclear. Variation in gene expression may contribute to these characteristics and enable the bacteria to respond to changing environmental conditions. Gene expression is controlled primarily at the level of transcription, where the main element of regulation is the promoter. Transcriptional regulation and associated promoter sequences have been studied extensively in E. coli. This review describes the complex structure and characteristics of mycobacterial promoters, in comparison to the classical E. coli prokaryotic promoter structure. Some components of mycobacterial promoters are similar to those of E. coli. These include the predominant guanine residue at the transcriptional start point, conserved -10 hexamer, similar interhexameric distances, the use of ATG as a start codon, the guanine- and adenine-rich ribosome binding site and the presence of extended -10 (TGn) motifs in strong promoters. However, these components are much more variable in sequence in mycobacterial promoters and no conserved -35 hexamer sequence (clearly defined in E. coli) can be identified. This may be a result of the high G+C content of mycobacterial genomes, as well as the large number of sigma factors present in mycobacteria, which may recognise different promoter sequences. Mycobacteria possess a complex transcriptional regulatory network. Numerous regulatory motifs have been identified in mycobacterial promoters, predominantly in the interhexameric region. These are bound by specific transcriptional regulators in response to environmental changes. The combination of specific promoter sequences, transcriptional regulators and a variety of sigma factors enables rapid and specific responses to diverse conditions and different stages of infection. This review aims to provide an overview of the complex architecture of mycobacterial transcriptional regulation.

A bacterial genome in transition--an exceptional enrichment of IS elements but lack of evidence for recent transposition in the symbiont Amoebophilus asiaticus

Background

Insertion sequence (IS) elements are important mediators of genome plasticity and are widespread among bacterial and archaeal genomes. The 1.88 Mbp genome of the obligate intracellular amoeba symbiont Amoebophilus asiaticus contains an unusually large number of transposase genes (n = 354; 23% of all genes).

Results

The transposase genes in the A. asiaticus genome can be assigned to 16 different IS elements termed ISCaa1 to ISCaa16, which are represented by 2 to 24 full-length copies, respectively. Despite this high IS element load, the A. asiaticus genome displays a GC skew pattern typical for most bacterial genomes, indicating that no major rearrangements have occurred recently. Additionally, the high sequence divergence of some IS elements, the high number of truncated IS element copies (n = 143), as well as the absence of direct repeats in most IS elements suggest that the IS elements of A. asiaticus are transpositionally inactive. Although we could show transcription of 13 IS elements, we did not find experimental evidence for transpositional activity, corroborating our results from sequence analyses. However, we detected contiguous transcripts between IS elements and their downstream genes at nine loci in the A. asiaticus genome, indicating that some IS elements influence the transcription of downstream genes, some of which might be important for host cell interaction.

Conclusions

Taken together, the IS elements in the A. asiaticus genome are currently in the process of degradation and largely represent reflections of the evolutionary past of A. asiaticus in which its genome was shaped by their activity.

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

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

Relaxed natural selection alone does not permit transposable element expansion within 4,000 generations in Escherichia coli

Insertion sequences (ISs) are transposable genetic elements in bacterial genomes. IS elements are common among bacteria but are generally rare within free-living species, probably because of the negative fitness effects they have on their hosts. Conversely, ISs frequently proliferate in intracellular symbionts and pathogens that recently transitioned from a free-living lifestyle. IS elements can profoundly influence the genomic evolution of their bacterial hosts, although it is unknown why they often expand in intracellular bacteria. We designed a laboratory evolution experiment with Escherichia coli K-12 to test the hypotheses that IS elements often expand in intracellular bacteria because of relaxed natural selection due to (1) their generally small effective population sizes (N (e)) and thus enhanced genetic drift, and (2) their nutrient rich environment, which makes many biosynthetic genes unnecessary and thus selectively neutral territory for IS insertion. We propagated 12 populations under four experimental conditions: large N (e) versus small N (e), and nutrient rich medium versus minimal medium. We found that relaxed selection over 4,000 generations was not sufficient to permit IS element expansion in any experimental population, thus leading us to hypothesize that IS expansion in intracellular symbionts may often be spurred by enhanced transposition rates, possibly due to environmental stress, coupled with relaxed natural selection.

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

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

Intergenic transposable elements are not randomly distributed in bacteria

Insertion sequences (ISs) are mobile genetic elements in bacterial genomes. In general, intergenic IS elements are probably less deleterious for their hosts than intragenic ISs, simply because they have a lower likelihood of disrupting native genes. However, since promoters, Shine–Dalgarno sequences, and transcription factor binding sites are intergenic and upstream of genes, I hypothesized that not all neighboring gene orientations (NGOs) are selectively equivalent for IS insertion. To test this, I analyzed the NGOs of all intergenic ISs in 326 fully sequenced bacterial chromosomes. Of the 116 genomes with enough IS elements for statistical analysis, 68 have significantly more ISs between convergently oriented genes than expected, and 46 have significantly fewer ISs between divergently oriented genes. This suggests that natural selection molds intergenic IS distributions because they are least intrusive between convergent gene pairs and most intrusive between divergent gene pairs.

A major subgroup of Beijing familyMycobacterium tuberculosisis associated with multidrug resistance and increased transmissibility

This study investigated further the association between the Beijing familyMycobacterium tuberculosiscirculating in rural China and anti-tuberculosis (TB) drug resistance. In total, 351M. tuberculosisisolates were collected through a population-based epidemiological study, 223 (63·5%) of which were resistant to at least one anti-TB drug, including 53 (15·1%) multidrug-resistant (MDR) isolates. Spoligotyping found 243 isolates (69·2%) that belonged to the Beijing family. A major subgroup of the Beijing family identified by mycobacterial interspersed repetitive unit (MIRU) genotyping (223325173533), showed significantly higher frequencies of MDR (44·7%vs. 13·7%, OR 6·18, 95% CI 2·68–14·23),katGandrpoBmutations (31·6%vs. 9·3%, OR 4·27, 95% CI 1·86–9·80), and being clustered by IS6110RFLP genotyping (60·5%vs. 21·0%, OR 6·14, 95% CI 2·82–13·37) in comparison with other Beijing family isolates. Our data suggest that MIRU genotype 223325173533 of the Beijing family is associated with MDR and increased transmissibility.

Allelic exchange and mutant selection demonstrate that common clinical embCAB gene mutations only modestly increase resistance to ethambutol in Mycobacterium tuberculosis

Mutations within codon 306 of the Mycobacterium tuberculosis embB gene modestly increase ethambutol (EMB) MICs. To identify other causes of EMB resistance and to identify causes of high-level resistance, we generated EMB-resistant M. tuberculosis isolates in vitro and performed allelic exchange studies of embB codon 406 (embB406) and embB497 mutations. In vitro selection produced mutations already identified clinically in embB306, embB397, embB497, embB1024, and embC13, which result in EMB MICs of 8 or 14 microg/ml, 5 microg/ml, 12 microg/ml, 3 microg/ml, and 4 microg/ml, respectively, and mutations at embB320, embB324, and embB445, which have not been identified in clinical M. tuberculosis isolates and which result in EMB MICs of 8 microg/ml, 8 microg/ml, and 2 to 8 microg/ml, respectively. To definitively identify the effect of the common clinical embB497 and embB406 mutations on EMB susceptibility, we created a series of isogenic mutants, exchanging the wild-type embB497 CAG codon in EMB-susceptible M. tuberculosis strain 210 for the embB497 CGG codon and the wild-type embB406 GGC codon for either the embB406 GCC, embB406 TGC, embB406 TCC, or embB406 GAC codon. These new mutants showed 6-fold and 3- to 3.5-fold increases in the EMB MICs, respectively. In contrast to the embB306 mutants, the isogenic embB497 and embB406 mutants did not have preferential growth in the presence of isoniazid or rifampin (rifampicin) at their MICs. These results demonstrate that individual embCAB mutations confer low to moderate increases in EMB MICs. Discrepancies between the EMB MICs of laboratory mutants and clinical M. tuberculosis strains with identical mutations suggest that clinical EMB resistance is multigenic and that high-level EMB resistance requires mutations in currently unknown loci.

Evidence for a rapid rate of molecular evolution at the hypervariable and immunogenic Mycobacterium tuberculosis PPE38 gene region

BACKGROUND

PPE38 (Rv2352c) is a member of the large PPE gene family of Mycobacterium tuberculosis and related mycobacteria. The function of PPE proteins is unknown but evidence suggests that many are cell-surface associated and recognised by the host immune system. Previous studies targeting other PPE gene members suggest that some display high levels of polymorphism and it is thought that this might represent a means of providing antigenic variation. We have analysed the genetic variability of the PPE38 genomic region on a cohort of M. tuberculosis clinical isolates representing all of the major phylogenetic lineages, along with the ancestral M. tuberculosis complex (MTBC) member M. canettii, and supplemented this with analysis of publicly available whole genome sequences representing additional M. tuberculosis clinical isolates, other MTBC members and non tuberculous mycobacteria (NTM). Where possible we have extended this analysis to include the adjacent plcABC and PPE39/40 genomic regions.

RESULTS

We show that the ancestral MTBC PPE38 region comprises 2 homologous PPE genes (PPE38 and PPE71), separated by 2 esat-6 (esx)-like genes and that this structure derives from an esx/esx/PPE duplication in the common ancestor of M. tuberculosis and M. marinum. We also demonstrate that this region of the genome is hypervariable due to frequent IS6110 integration, IS6110-associated recombination, and homologous recombination and gene conversion events between PPE38 and PPE71. These mutations result in combinations of gene deletion, gene truncation and gene disruption in the majority of clinical isolates. These mutations were generally found to be IS6110 strain lineage-specific, although examples of additional within-lineage and even within-cluster mutations were observed. Furthermore, we provide evidence that the published M. tuberculosis H37Rv whole genome sequence is inaccurate regarding this region.

CONCLUSION

Our results show that this antigen-encoding region of the M. tuberculosis genome is hypervariable. The observation that numerous different mutations have become fixed within specific lineages demonstrates that this genomic region is undergoing rapid molecular evolution and that further lineage-specific evolutionary expansion and diversification has occurred subsequent to the lineage-defining mutational events. We predict that functional loss of these genes could aid immune evasion. Finally, we also show that the PPE38 region of the published M. tuberculosis H37Rv whole genome sequence is not representative of the ATCC H37Rv reference strain.

Global transcriptional response to spermine, a component of the intramacrophage environment, reveals regulation of Francisella gene expression through insertion sequence elements

Tularemia is caused by the category A biodefense agent Francisella tularensis. This bacterium is associated with diverse environments and a plethora of arthropod and mammalian hosts. How F. tularensis adapts to these different conditions, particularly the eukaryotic intracellular environment in which it replicates, is poorly understood. Here, we demonstrate that the polyamines spermine and spermidine are environmental signals that alter bacterial stimulation of host cells. Genomewide analysis showed that F. tularensis LVS undergoes considerable changes in gene expression in response to spermine. Unexpectedly, analysis of gene expression showed that multiple members of two classes of Francisella insertion sequence (IS) elements, ISFtu1 and ISFtu2, and the genes adjacent to these elements were induced by spermine. Spermine was sufficient to activate transcription of these IS elements and of nearby genes in broth culture and in macrophages. Importantly, the virulent strain of F. tularensis, Schu S4, exhibited similar phenotypes of cytokine induction and gene regulation in response to spermine. Distinctions in gene expression changes between Schu S4 and LVS at one orthologous locus, however, correlated with differences in IS element location. Our results indicate that spermine and spermidine are novel triggers to alert F. tularensis of its eukaryotic host environment. The results reported here also identify an unexpected mechanism of gene regulation controlled by a spermine-responsive promoter contained within IS elements. Different arrangements of these mobile genetic elements among Francisella strains may contribute to virulence by conveying new expression patterns for genes from different strains.

A novel insertion sequence, IS1642, of Mycobacterium avium, which forms long direct repeats of variable length

A new insertion sequence (IS), IS1642, was identified in a Mycobacterium avium strain isolated from a human patient. IS1642 had a size of 1642 bp and contained a single ORF encoding a probable transposase of 503 amino acid residues homologous (79% identity) to that of IS1549 found in Mycobacterium smegmatis. The IS1642 included imperfect inverted repeats (5'-cctgacttttatca-3', 5'-tgataaaagtcggg-3') on its ends, and was flanked by direct repeats of variable length ranging from 5 to 161 bp. It was suggested that the IS1642 was widely distributed in many M. avium strains of human patients, and the Southern blot profile of IS1642 was very diverse among the strains examined. The transposition event of IS1642 was observed by in vitro repeated passages, showing that the IS1642 is actually a transposable element. In light of these characteristics, IS1642 could be a new useful marker when genotyping with high discrimination is required.

Analysis of the genetic variation in Mycobacterium tuberculosis strains by multiple genome alignments

Background

The recent determination of the complete nucleotide sequence of several Mycobacterium tuberculosis (MTB) genomes allows the use of comparative genomics as a tool for dissecting the nature and consequence of genetic variability within this species. The multiple alignment of the genomes of clinical strains (CDC1551, F11, Haarlem and C), along with the genomes of laboratory strains (H37Rv and H37Ra), provides new insights on the mechanisms of adaptation of this bacterium to the human host.

Findings

The genetic variation found in six M. tuberculosis strains does not involve significant genomic rearrangements. Most of the variation results from deletion and transposition events preferentially associated with insertion sequences and genes of the PE/PPE family but not with genes implicated in virulence. Using a Perl-based software islandsanalyser, which creates a representation of the genetic variation in the genome, we identified differences in the patterns of distribution and frequency of the polymorphisms across the genome. The identification of genes displaying strain-specific polymorphisms and the extrapolation of the number of strain-specific polymorphisms to an unlimited number of genomes indicates that the different strains contain a limited number of unique polymorphisms.

Conclusion

The comparison of multiple genomes demonstrates that the M. tuberculosis genome is currently undergoing an active process of gene decay, analogous to the adaptation process of obligate bacterial symbionts. This observation opens new perspectives into the evolution and the understanding of the pathogenesis of this bacterium.

IS6110 in oriC affects the morphology and growth of Mycobacterium tuberculosis and attenuates virulence in mice

The IS6110 element is widely used in studies of molecular epidemiology of tuberculosis and it is considered the gold standard for genotyping Mycobacterium tuberculosis strains. Because of its high frequency of transposition, IS6110 is probably a major contributor to the evolution of M. tuberculosis. Nevertheless, very few studies of the effect of IS6110 insertions on the virulence of M. tuberculosis have been reported. We analysed two isogenic groups of M. tuberculosis strains isolated from the sputa of two patients. Strains belonging to the same isogenic group differed from one another by one IS6110-oriC hybridising band, but they showed identical spoligo and MIRU-VNTR profiles. Isogenic strains containing the IS6110 element in oriC exhibited a diminished growth rate and average dimensions of the bacilli were modified; moreover, they were less virulent in a mouse model.

Spontaneous transposition of IS1096 or ISMsm3 leads to glycopeptidolipid overproduction and affects surface properties in Mycobacterium smegmatis

Natural modification of the colony appearance is a phenomenon that has not been fully understood in mycobacteria. Here, we show that Mycobacterium smegmatis ATCC607 displays a low-frequency spontaneous morphological variation that correlates with the acquisition of a panel of new phenotypes, such as aggregation, biofilm formation and sliding motility. These variants produce larger amounts of glycopeptidolipid (GPL), a cell-surface component, than did the wild-type strain. This conversion results from the transposition of two types of insertion sequences, IS1096 and ISMsm3, into two loci. One locus is the promoter region of the mps operon, the GPL biosynthesis gene cluster, leading to the overexpression of these genes. The other locus is the lsr2 gene, which encodes a small basic histone-like protein that likely plays a regulatory role at the mps promoter and also controls pigment production. This study demonstrates that insertion sequence mobility play a crucial role in the acquisition of new phenotypes.

IS1096-mediated DNA rearrangements play a key role in genome evolution of Mycobacterium smegmatis

The acquisition of DNA and the loss of genetic information are two important mechanisms that contribute to strain-specific differences in genome content. In this study, comparative genomics has allowed us to infer the roles of genomic rearrangement and changes in both distribution and copy number of the insertion element, IS1096, in the evolution of Mycobacterium smegmatis mc2155 from its progenitor, M. smegmatis ATCC 607. Comparative analysis revealed that the ATCC 607 genome contains only 11 IS1096 elements against the 24 reported in mc2155. As mc2155 evolved, there was a considerable expansion in the copy number of IS1096 (+13) as well as duplication of a 56-kb fragment flanked on both sides by IS1096; concurrently, a single IS1096 element and its flank were deleted. This study demonstrates that insertion sequence (IS) expansion and IS-induced rearrangements such as duplication, deletion and shuffling are major forces driving genomic diversity and evolution.

Transfer of embB codon 306 mutations into clinical Mycobacterium tuberculosis strains alters susceptibility to ethambutol, isoniazid, and rifampin

Implicated as a major mechanism of ethambutol (EMB) resistance in clinical studies of Mycobacterium tuberculosis, mutations in codon 306 of the embB gene (embB306) have also been detected in EMB-susceptible clinical isolates. Other studies have found strong associations between embB306 mutations and multidrug resistance, but not EMB resistance. We performed allelic exchange studies in EMB-susceptible and EMB-resistant clinical M. tuberculosis isolates to identify the role of embB306 mutations in any type of drug resistance. Replacing wild-type embB306 ATG from EMB-susceptible clinical M. tuberculosis strain 210 with embB306 ATA, ATC, CTG, or GTG increased the EMB MIC from 2 microg/ml to 7, 7, 8.5, and 14 microg/ml, respectively. Replacing embB306 ATC or GTG from two high-level EMB-resistant clinical strains with wild-type ATG lowered EMB MICs from 20 microg/ml or 28 microg/ml, respectively, to 3 microg/ml. All parental and isogenic mutant strains had identical isoniazid (INH) and rifampin (RIF) MICs. However, embB306 CTG mutants had growth advantages compared to strain 210 at sub-MICs of INH or RIF in monocultures and at sub-MICs of INH in competition assays. CTG mutants were also more resistant to the additive or synergistic activities of INH, RIF, or EMB used in different combinations. These results demonstrate that embB306 mutations cause an increase in the EMB MIC, a variable degree of EMB resistance, and are necessary but not sufficient for high-level EMB resistance. The unusual growth property of embB306 mutants in other antibiotics suggests that they may be amplified during treatment in humans and that a single mutation may affect antibiotic susceptibility against multiple first-line antibiotics.