IS6110 mediates increased transcription of the phoP virulence gene in a multidrug-resistant clinical isolate responsible for tuberculosis outbreaks

Delineating the Acquired Genetic Diversity and Multidrug Resistance in Alcaligenes from Poultry Farms and Nearby Soil

Alcaligenes faecalis is one of the most important and clinically significant environmental pathogens, increasing in importance due to its isolation from soil and nosocomial environments. The Gram-negative soil bacterium is associated with skin endocarditis, bacteremia, dysentery, meningitis, endophthalmitis, urinary tract infections, and pneumonia in patients. With emerging antibiotic resistance in A. faecalis, it has become crucial to understand the origin of such resistance genes within this clinically significant environmental and gut bacterium. In this research, we studied the impact of antibiotic overuse in poultry and its effect on developing resistance in A. faecalis. We sampled soil and faecal materials from five poultry farms, performed whole genome sequencing & analysis and identified four strains of A. faecalis. Furthermore, we characterized the genes in the genomic islands of A. faecalis isolates. We found four multidrug-resistant A. faecalis strains that showed resistance against vancomycin (MIC >1000 μg/ml), ceftazidime (50 μg/ml), colistin (50 μg/ml) and ciprofloxacin (50 μg/ml). From whole genome comparative analysis, we found more than 180 resistance genes compared to the reference sequence. Parts of our assembled contigs were found to be similar to different bacteria which included pbp1A and pbp2 imparting resistance to amoxicillin originally a part of Helicobacter and Bordetella pertussis. We also found the Mycobacterial insertion element IS6110 in the genomic islands of all four genomes. This prominent insertion element can be transferred and induce resistance to other bacterial genomes. The results thus are crucial in understanding the transfer of resistance genes in the environment and can help in developing regimes for antibiotic use in the food and poultry industry.

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.

Characterisation of IS1311 in Mycobacterium avium subspecies paratuberculosis genomes: Typing, continental clustering, microbial evolution and host adaptation

Johne's disease (JD), caused by Mycobacterium avium subspecies paratuberculosis (MAP) is a global burden for livestock producers and has an association with Crohn's disease in humans. Within MAP there are two major lineages, S/Type I/TypeIII and C/Type II, that vary in phenotype including culturability, host preference and virulence. These lineages have been identified using the IS1311 element, which contains a conserved, single nucleotide polymorphism. IS1311 and the closely related IS1245 element belong to the IS256 family of insertion sequences, are dispersed throughout M. avium taxa but remain poorly characterised. To investigate the distribution and diversity of IS1311 in MAP, 805 MAP genomes were collated from public databases. IS1245 was absent, while IS1311 sequence, copy number and insertion loci were conserved between MAP S lineages and varied within the MAP C lineage. One locus was specific to the S strains, which contained nine IS1311 copies. In contrast, C strains contained either seven or eight IS1311 loci. Most insertion loci were associated with the boundaries of homologous regions that had undergone genome rearrangement between the MAP lineages, suggesting that this sequence may be a driver of recombination. Phylogenomic geographic clustering of MAP subtypes was demonstrated for the first time, at continental scale, and indicated that there may have been recent MAP transmission between Europe and North America, in contrast to Australia where importation of live ruminants is generally prohibited. This investigation confirmed the utility of IS1311 typing in epidemiological studies and resolved anomalies in past studies. The results shed light on potential mechanisms of niche/host adaptation, virulence of MAP and global transmission dynamics.

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.

Recent progress in the genotyping of bovine tuberculosis and its rapid diagnosis via nanoparticle-based electrochemical biosensors

Bovine tuberculosis (bTB) is considered a worldwide infectious zoonotic disease. Mycobacterium bovis causes bTB disease. It is one of the Mycobacterium tuberculosis complex (MTBC) members. MTBC is a clonal complex of close relatives with approximately 99.95% similarity. M. bovis is a spillover pathogen that can transmit from animals to humans and rarely from humans to animals with contact. Genotyping techniques are important to discriminate and differentiate between MTBC species. Spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) are widely used but they have some limitations. As an alternative, whole genome sequencing approaches have been utilized due to their high-resolution power. They are employed in typing M. bovis and explain the evolutionary and phylogenetic relationships between isolates. The control of bTB disease has attracted a large amount of attention. Rapid and proper diagnosis is necessary for monitoring the disease as an initial step for its control and treatment. Nanotechnology has a potential impact on the rapid diagnosis and treatment of bTB through the use of nanocarrier and metal nanoparticles (NPs). Special attention has been paid to voltammetric and impedimetric electrochemical strategies as facile, sensitive, and selective methods for the efficient detection of tuberculosis. The efficacy of these sensors is enhanced in the presence of NPs, which act as recognition and/or redox probes. Gold, silver, copper, cobalt, graphene, and magnetic NPs, as well as polypyrrole nanowires and multiwalled carbon nanotubes have been employed for detecting tuberculosis. Overall, NP-based electrochemical sensors represent a promising tool for the diagnosis of bTB.

ESAT-6 a Major Virulence Factor of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb), the causative agent of human tuberculosis (TB), is one of the most successfully adapted human pathogens. Human-to-human transmission occurs at high rates through aerosols containing bacteria, but the pathogen evolved prior to the establishment of crowded populations. Mtb has developed a particular strategy to ensure persistence in the host until an opportunity for transmission arises. It has refined its lifestyle to obviate the need for virulence factors such as capsules, flagella, pili, or toxins to circumvent mucosal barriers. Instead, the pathogen uses host macrophages, where it establishes intracellular niches for its migration into the lung parenchyma and other tissues and for the induction of long-lived latency in granulomas. Finally, at the end of the infection cycle, Mtb induces necrotic cell death in macrophages to escape to the extracellular milieu and instructs a strong inflammatory response that is required for the progression from latency to disease and transmission. Common to all these events is ESAT-6, one of the major virulence factors secreted by the pathogen. This narrative review highlights the recent advances in understanding the role of ESAT-6 in hijacking macrophage function to establish successful infection and transmission and its use as a target for the development of diagnostic tools and vaccines.

Effects of UV-B radiation on epiphytic bacterial communities on male and female Sargassum thunbergii

The effects of increased UV-B radiation on macroalgae have been widely studied, but knowledge concerning the response of communities of algal epiphytic bacteria to increased UV-B radiation and differences between male and female algae is still lacking. Via 16S rDNA high-throughput sequencing technology, changes in the epiphytic bacterial communities on male and female S. thunbergii under increased UV-B radiation were studied in the lab. Under different UV-B radiation intensities, although the α diversity and community composition of epiphytic bacteria changed little, the β diversity indicated that the community structure of bacteria on S. thunbergii was obviously clustered, and the relative abundance of dominant bacteria and indicator species changed considerably. There were unique bacteria in each experimental group, and the bacteria whose abundance obviously changed were members of groups related to environmental resistance or adaptability. The variation in the abundance of epiphytic bacteria was different in male and female S. thunbergii, and the bacteria whose abundance greatly changed were mainly related to algal growth and metabolism. The abundance of genes with predicted functions related to metabolism, genetic information processing, environmental adaptation and infectious diseases changed with increased UV-B radiation, and those variations differed between epiphytic bacteria on male and female S. thunbergii. This study found that the algal epiphytic bacteria were influenced by the increase in UV-B radiation and underwent certain adaptations through adjustments to community structure and function, and this response was also affected by the sex of the macroalgae. These results are expected to serve as experimental basis and provide reference for further understanding of the response of algae epiphytic bacteria to enhanced UV-B radiation caused by the thinning of the ozone layer and the resulting changes in the relationship between algae and bacteria, which may change the community of the marine ecosystem and affect important marine ecological process.

Infectious and Inflammatory Pathways to Cough

Coughing is a dynamic physiological process resulting from input of vagal sensory neurons innervating the airways and perceived airway irritation. Although cough serves to protect and clear the airways, it can also be exploited by respiratory pathogens to facilitate disease transmission. Microbial components or infection-induced inflammatory mediators can directly interact with sensory nerve receptors to induce a cough response. Analysis of cough-generated aerosols and transmission studies have further demonstrated how infectious disease is spread through coughing. This review summarizes the neurophysiology of cough, cough induction by respiratory pathogens and inflammation, and cough-mediated disease transmission.

Features of Mycobacterium bovis Complete Genomes Belonging to 5 Different Lineages

Mammalian tuberculosis (TB) is a zoonotic disease mainly due to Mycobacterium bovis (M. bovis). A current challenge for its eradication is understanding its transmission within multi-host systems. Improvements in long-read sequencing technologies have made it possible to obtain complete bacterial genomes that provide a comprehensive view of species-specific genomic features. In the context of TB, new genomic references based on complete genomes genetically close to field strains are also essential to perform precise field molecular epidemiological studies. A total of 10 M. bovis strains representing each genetic lineage identified in France and in other countries were selected for performing complete assembly of their genomes. Pangenome analysis revealed a "closed" pangenome composed of 3900 core genes and only 96 accessory genes. Whole genomes-based alignment using progressive Mauve showed remarkable conservation of the genomic synteny except that the genomes have a variable number of copies of IS6110. Characteristic genomic traits of each lineage were identified through the discovery of specific indels. Altogether, these results provide new genetic features that improve the description of M. bovis lineages. The availability of new complete representative genomes of M. bovis will be useful to epidemiological studies and better understand the transmission of this clonal-evolving pathogen.

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.

Comparative Genomics of Disease and Carriage Serotype 1 Pneumococci

The isolation of Streptococcus pneumoniae serotypes in systemic tissues of patients with invasive disease versus the nasopharynx of healthy individuals with asymptomatic carriage varies widely. Some serotypes are hyper-invasive, particularly serotype 1, but the underlying genetics remain poorly understood due to the rarity of carriage isolates, reducing the power of comparison with invasive isolates. Here, we use a well-controlled genome-wide association study to search for genetic variation associated with invasiveness of serotype 1 pneumococci from a serotype 1 endemic setting in Africa. We found no consensus evidence that certain genomic variation is overrepresented among isolates from patients with invasive disease than asymptomatic carriage. Overall, the genomic variation explained negligible phenotypic variability, suggesting a minimal effect on the disease status. Furthermore, changes in lineage distribution were seen with lineages replacing each other over time, highlighting the importance of continued pathogen surveillance. Our findings suggest that the hyper-invasiveness is an intrinsic property of the serotype 1 strains, not specific for a "disease-associated" subpopulation disproportionately harboring unique genomic variation.

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.

Frequency, composition and mobility of Escherichia coli-derived transposable elements in holdings of plasmid repositories

By providing the scientific community with uniform and standardized resources of consistent quality, plasmid repositories play an important role in enabling scientific reproducibility. Plasmids containing insertion sequence elements (IS elements) represent a challenge from this perspective, as they can change the plasmid structure and function. In this study, we conducted a systematic analysis of a subset of plasmid stocks distributed by plasmid repositories (The Arabidopsis Biological Resource Center and Addgene) which carry unintended integrations of bacterial mobile genetic elements. The integration of insertion sequences was most often found in, but not limited to, pBR322-derived vectors, and did not affect the function of the specific plasmids. In certain cases, the entire stock was affected, but the majority of the stocks tested contained a mixture of the wild-type and the mutated plasmids, suggesting that the acquisition of IS elements likely occurred after the plasmids were acquired by the repositories. However, comparison of the sequencing results of the original samples revealed that some plasmids already carried insertion mutations at the time of donation. While an extensive BLAST analysis of 47 877 plasmids sequenced from the Addgene repository uncovered IS elements in only 1.12%, suggesting that IS contamination is not widespread, further tests showed that plasmid integration of IS elements can propagate in conventional Escherichia coli hosts over a few tens of generations. Use of IS-free E. coli hosts prevented the emergence of IS insertions as well as that of small indels, suggesting that the use of IS-free hosts by donors and repositories could help limit unexpected and unwanted IS integrations into plasmids.

A Low-Prevalence Single-Nucleotide Polymorphism in the Sensor Kinase PhoR in Mycobacterium tuberculosis Suppresses Its Autophosphatase Activity and Reduces Pathogenic Fitness: Implications in Evolutionary Selection

The genome sequencing of Mycobacterium tuberculosis, the causative organism of tuberculosis, has significantly improved our understanding of the mechanisms that drive the establishment of infection and disease progression. Several clinical strains of M. tuberculosis exhibit single-nucleotide polymorphisms (SNPs), the implications of which are only beginning to be understood. Here, we examined the impact of a specific polymorphism in PhoR, the sensor kinase of the PhoPR two-component system. Biochemical analysis revealed reduced autophosphatase/ATPase activity, which led to enhanced downstream gene expression. We complemented M. tuberculosis H37Ra with the wild-type and mutant phoPR genes and characterized the strains in a cell line infection model. We provide an explanation for the low prevalence of the SNP in clinical strains (∼1%), as the mutation causes a survival disadvantage in the host cells. The study provides a rare example of selection of a signaling node under competing evolutionary forces, wherein a biochemically superior mutation aids bacterial adaptation within-host but has low fitness for infection and hence is not selected. Our study highlights the importance of accounting for such SNPs to test therapeutic and co-therapeutic methods to combat TB.

Intragenic Distribution of IS6110 in Clinical Mycobacterium tuberculosis Strains: Bioinformatic Evidence for Gene Disruption Leading to Underdiagnosed Antibiotic Resistance

Antibiotic resistance is a global challenge for tuberculosis control, and accelerating its diagnosis is critical for therapy decisions and controlling transmission. Genotype-based molecular diagnostics now play an increasing role in accelerating the detection of such antibiotic resistance, but their accuracy depends on the instructed detection of genetic variations. Genetic mobile elements such as IS6110 are established sources of genetic variation in Mycobacterium tuberculosis, but their implication in clinical antibiotic resistance has thus far been unclear. Here, we describe the discovery of an intragenic IS6110 insertion into Rv0678 that caused antibiotic resistance in an in vitro-selected M. tuberculosis isolate. The subsequent development of bioinformatics scripts allowed genome-wide analysis of intragenic IS6110 insertions causing gene disruptions in 6,426 clinical M. tuberculosis strains. This analysis identified 10,070 intragenic IS6110 insertions distributed among 333 different genes. Focusing on genes whose disruption leads to antibiotic resistance, 12 clinical isolates were identified with high confidence to be resistant to bedaquiline, clofazimine, pyrazinamide, ethionamide, and para-aminosalicylic acid because of an IS6110-mediated gene disruption event. A number of these IS6110-mediated resistant strains had identical genomic distributions of IS6110 elements and likely represent transmission events of a single resistant isolate. These data provide strong evidence that IS6110-mediated gene disruption is a clinically relevant mechanism of antibiotic resistance in M. tuberculosis that should be considered for molecular diagnostics. Concomitantly, this analysis provides a list of 333 IS6110-disrupted genes in clinical tuberculosis isolates that can be deemed nonessential for human infection.

IMPORTANCE To help control the spread of drug-resistant tuberculosis and to guide treatment choices, it is important that rapid and accurate molecular diagnostic tools are used. Current molecular diagnostic tools detect the most common antibiotic-resistance-conferring mutations in the form of single nucleotide changes, small deletions, or insertions. Mobile genetic elements, named IS6110, are also known to move within the M. tuberculosis genome and cause significant genetic variations, although the role of this variation in clinical drug resistance remains unclear. In this work, we show that both in vitro and in data analyzed from 6,426 clinical M. tuberculosis strains, IS6110 elements are found that disrupt specific genes essential for the function of a number of pivotal antituberculosis drugs. By providing ample evidence of clinically relevant IS6110-mediated drug resistance, we believe that this shows that this form of genetic variation must not be overlooked in molecular diagnostics of drug resistance.

Analysis of Mycobacterium africanum in the last 17 years in Aragon identifies a specific location of IS6110 in Lineage 6

The purpose of this study was to increase our knowledge about Mycobacterium africanum and report the incidence and characteristics of tuberculosis (TB) due to their lineages in Aragon, Spain, over the period 2003–2019. The study includes all the cases in our region, where all the M. tuberculosis complex isolates are systematically characterised. We detected 31 cases of M. africanum among 2598 cases of TB in the period studied. TB caused by M. africanum is rare (1.19%) in our population, and it affects mainly men of economically productive age coming from West African countries. Among the isolates, Lineage (L) 6 was more frequent than L5. The genotyping of these strains identified five clusters and 13 strains with a unique pattern. The isolates’ characterisation identified a copy of IS6110 within the moaX gene, which turned out to be specific for L6. It will allow the differentiation of this lineage from the rest of MTBC with a simple PCR reaction. It remains to be established whether this polymorphism may limit M. africanum transmission. Furthermore, a mutation in the mutT2 promoter was found as specific for L6 strains, which could be related to the high variability found for L6 compared to L5.

The role of Mycobacterium tuberculosis complex species on apoptosis and necroptosis state of macrophages derived from active pulmonary tuberculosis patients

Objective

The role of Mycobacterium tuberculosis complex (MTBC) species in tuberculosis (TB) infection in human is still questioned. The aim of this study was to determine whether M. tuberculosis and M. bovis is associated with apoptosis and necroptosis by measuring the expression of specific signaling pathways components (Fas-associated protein with death domain (FADD) and receptor interacting protein 3 (RIP3)), and the level of apoptosis.

Results

We recruited 30 patients with pulmonary TB; 24 patients were infected with M. tuberculosis Beijing strain and six patients with M. bovis BCG strain. M. tuberculosis-infected patients were more likely to have severe lung damage compared to those infected with M. bovis (odds ratio [OR] 7.60; 95% confidence interval [CI] 1.07–54.09). M. tuberculosis infection was associated with lower expression of FADD and lower apoptosis level of macrophages compared to M. bovis. No significant different of RIP3 between MTBC species groups. In conclusion, M. tuberculosis Beijing strain was associated with severe pulmonary damage, inhibited FADD expression and reduced apoptosis level of macrophages derived from pulmonary TB patients. This suggests that the M. tuberculosis Beijing strain is potentially to be used as determinant of disease progressivity and tissue damage in TB cases.

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

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

Comparing the Metabolic Capabilities of Bacteria in the Mycobacterium tuberculosis Complex

Pathogenic mycobacteria are known for their ability to maintain persistent infections in various mammals. The canonical pathogen in this genus is Mycobacterium tuberculosis and this bacterium is particularly successful at surviving and replicating within macrophages. Here, we will highlight the metabolic processes that M. tuberculosis employs during infection in macrophages and compare these findings with what is understood for other pathogens in the M. tuberculosis complex.

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

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

RD5-mediated lack of PE_PGRS and PPE-MPTR export in BCG vaccine strains results in strong reduction of antigenic repertoire but little impact on protection

Tuberculosis is the deadliest infectious disease worldwide. Although the BCG vaccine is widely used, it does not efficiently protect against pulmonary tuberculosis and an improved tuberculosis vaccine is therefore urgently needed. Mycobacterium tuberculosis uses different ESX/Type VII secretion (T7S) systems to transport proteins important for virulence and host immune responses. We recently reported that secretion of T7S substrates belonging to the mycobacteria-specific Pro-Glu (PE) and Pro-Pro-Glu (PPE) proteins of the PGRS (polymorphic GC-rich sequences) and MPTR (major polymorphic tandem repeat) subfamilies required both a functional ESX-5 system and a functional PPE38/71 protein for secretion. Inactivation of ppe38/71 and the resulting loss of PE_PGRS/PPE-MPTR secretion were linked to increased virulence of M. tuberculosis strains. Here, we show that a predicted total of 89 PE_PGRS/PPE-MPTR surface proteins are not exported by certain animal-adapted strains of the M. tuberculosis complex including M. bovis. This Δppe38/71-associated secretion defect therefore also occurs in the M. bovis-derived tuberculosis vaccine BCG and could be partially restored by introduction of the M. tuberculosis ppe38-locus. Epitope mapping of the PPE-MPTR protein PPE10, further allowed us to monitor T-cell responses in splenocytes from BCG/M. tuberculosis immunized mice, confirming the dependence of PPE10-specific immune-induction on ESX-5/PPE38-mediated secretion. Restoration of PE_PGRS/PPE-MPTR secretion in recombinant BCG neither altered global antigenic presentation or activation of innate immune cells, nor protective efficacy in two different mouse vaccination-infection models. This unexpected finding stimulates a reassessment of the immunomodulatory properties of PE_PGRS/PPE-MPTR proteins, some of which are contained in vaccine formulations currently in clinical evaluation.

One of the major findings of the pioneering Mycobacterium tuberculosis H37Rv genome sequencing project was the identification of the highly abundant PE and PPE proteins, named after their N-terminal motifs Pro–Glu (PE) or Pro–Pro–Glu (PPE). Within the 20 years of research since then, many claims were made that PE/PPE proteins, including the two large subgroups encoded by repetitive sequences with very high GC content (PE_PGRS and PPE-MPTR families), are exported to the bacterial surface or beyond, and show broad immunomodulatory impact on host-pathogen interaction. We thus screened strains from different branches of the M. tuberculosis complex, including the attenuated Mycobacterium bovis BCG vaccine strains, for their capacity to export PE_PGRS/PPE-MPTR proteins. Strikingly, we found that BCG strains were unable to export the plethora of PE_PGRS/PPE-MPTR proteins due to the absence of the region of difference RD5, which in M. tuberculosis encodes PPE38, required for PE_PGRS/PPE-MPTR export. Surprisingly, the restoration of PE_PGRS/PPE-MPTR export by genetic complementation in recombinant BCG did not result in immunomodulatory changes or altered protection in mouse models. Our results thus put into perspective the numerous reports on virulence-associated immunomodulatory impact of individual PE_PGRS and PPE-MPTR proteins and open novel questions on their biological function(s).

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.

Mycobacterium tuberculosis Complex Members Adapted to Wild and Domestic Animals

The Mycobacterium tuberculosis complex (MTBC) is composed of several highly genetically related species that can be broadly classified into those that are human-host adapted and those that possess the ability to propagate and transmit in a variety of wild and domesticated animals. Since the initial description of the bovine tubercle bacillus, now known as Mycobacterium bovis, by Theobald Smith in the late 1800's, isolates originating from a wide range of animal hosts have been identified and characterized as M. microti, M. pinnipedii, the Dassie bacillus, M. mungi, M. caprae, M. orygis and M. suricattae. This chapter outlines the events resulting in the identification of each of these animal-adapted species, their close genetic relationships, and how genome-based phylogenetic analyses of species-specific variation amongst MTBC members is beginning to unravel the events that resulted in the evolution of the MTBC and the observed host tropism between the human- and animal-adapted member species.

Comparative 'omics analyses differentiate Mycobacterium tuberculosis and Mycobacterium bovis and reveal distinct macrophage responses to infection with the human and bovine tubercle bacilli

Members of the Mycobacterium tuberculosis complex (MTBC) are the causative agents of tuberculosis in a range of mammals, including humans. A key feature of MTBC pathogens is their high degree of genetic identity yet distinct host tropism. Notably, while Mycobacterium bovis is highly virulent and pathogenic for cattle, the human pathogen M. tuberculosis is attenuated in cattle. Previous research also suggests that host preference amongst MTBC members has a basis in host innate immune responses. To explore MTBC host tropism, we present in-depth profiling of the MTBC reference strains M. bovis AF2122/97 and M. tuberculosis H37Rv at both the global transcriptional and the translational level via RNA-sequencing and SWATH MS. Furthermore, a bovine alveolar macrophage infection time course model was used to investigate the shared and divergent host transcriptomic response to infection with M. tuberculosis H37Rv or M. bovis AF2122/97. Significant differential expression of virulence-associated pathways between the two bacilli was revealed, including the ESX-1 secretion system. A divergent transcriptional response was observed between M. tuberculosis H37Rv and M. bovis AF2122/97 infection of bovine alveolar macrophages, in particular cytosolic DNA-sensing pathways at 48 h post-infection, and highlights a distinct engagement of M. bovis with the bovine innate immune system. The work presented here therefore provides a basis for the identification of host innate immune mechanisms subverted by virulent host-adapted mycobacteria to promote their survival during the early stages of infection.

Mycobacterium tuberculosis Acquires Limited Genetic Diversity in Prolonged Infections, Reactivations and Transmissions Involving Multiple Hosts

Background:Mycobacterium tuberculosis (MTB) has limited ability to acquire variability. Analysis of its microevolution might help us to evaluate the pathways followed to acquire greater infective success. Whole-genome sequencing (WGS) in the analysis of the transmission of MTB has elucidated the magnitude of variability in MTB. Analysis of transmission currently depends on the identification of clusters, according to the threshold of variability (<5 SNPs) between isolates. Objective: We evaluated whether the acquisition of variability in MTB, was more frequent in situations which could favor it, namely intrapatient, prolonged infections or reactivations and interpatient transmissions involving multiple sequential hosts. Methods: We used WGS to analyze the accumulation of variability in sequential isolates from prolonged infections or translations from latency to reactivation. We then measured microevolution in transmission clusters with prolonged transmission time, high number of involved cases, simultaneous involvement of latency and active transmission. Results: Intrapatient and interpatient acquisition of variability was limited, within the ranges expected according to the thresholds of variability proposed, even though bursts of variability were observed. Conclusions: The thresholds of variability proposed for MTB seem to be valid in most circumstances, including those theoretically favoring acquisition of variability. Our data point to multifactorial modulation of microevolution, although further studies are necessary to elucidate the factors underlying this modulation.

Comparative Whole-Genomic Analysis of an Ancient L2 Lineage Mycobacterium tuberculosis Reveals a Novel Phylogenetic Clade and Common Genetic Determinants of Hypervirulent Strains

Background: Development of improved therapeutics against tuberculosis (TB) is hindered by an inadequate understanding of the relationship between disease severity and genetic diversity of its causative agent, Mycobacterium tuberculosis. We previously isolated a hypervirulent M. tuberculosis strain H112 from an HIV-negative patient with an aggressive disease progression from pulmonary TB to tuberculous meningitis—the most severe manifestation of tuberculosis. Human macrophage challenge experiment demonstrated that the strain H112 exhibited significantly better intracellular survivability and induced lower level of TNF-α than the reference virulent strain H37Rv and other 123 clinical isolates.

Aim: The present study aimed to identify the potential genetic determinants of mycobacterial virulence that were common to strain H112 and hypervirulent M. tuberculosis strains of the same phylogenetic clade isolated in other global regions.

Methods: A low-virulent M. tuberculosis strain H54 which belonged to the same phylogenetic lineage (L2) as strain H112 was selected from a collection of 115 clinical isolates. Both H112 and H54 were whole-genome-sequenced using PacBio sequencing technology. A comparative genomics approach was adopted to identify mutations present in strain H112 but absent in strain H54. Subsequently, an extensive phylogenetic analysis was conducted by including all publically available M. tuberculosis genomes. Single-nucleotide-polymorphisms (SNPs) and structural variations (SVs) common to hypervirulent strains in the global collection of genomes were considered as potential genetic determinants of hypervirulence.

Results:Sequencing data revealed that both H112 and H54 were identified as members of the same sub-lineage L2.2.1. After excluding the lineage-related mutations shared between H112 and H54, we analyzed the phylogenetic relatedness of H112 with global collection of M. tuberculosis genomes (n = 4,338), and identified a novel phylogenetic clade in which four hypervirulent strains isolated from geographically diverse regions were clustered together. All hypervirulent strains in the clade shared 12 SNPs and 5 SVs with H112, including those affecting key virulence-associated loci, notably, a deleterious SNP (rv0178 p. D150E) within mce1 operon and an intergenic deletion (854259_ 854261delCC) in close-proximity to phoP.

Conclusion: The present study identified common genetic factors in a novel phylogenetic clade of hypervirulent M. tuberculosis. The causative role of these mutations in mycobacterial virulence should be validated in future study.

Comparative Genomics of a Bovine Mycobacterium tuberculosis Isolate and Other Strains Reveals Its Potential Mechanism of Bovine Adaptation

The Mycobacterium tuberculosis complex causes tuberculosis (TB) in humans and other animal species, but Mycobacterium tuberculosis has a distinct host preference to humans. The present study aimed to determine whether a bovine M. tb strain 1458 has evolved some genetic properties in their genome that might be associated with their bovine adaptation. The genome of the M. tb strain 1458 was sequenced and subjected to an extensive comparative genomic analysis. A phylogenetic analysis showed that strain 1458 is most closely related to a Chinese M. tb strain, CCDC5079, of the same Beijing family. Compared with three human M. tb Beijing family strains, the strain 1458 has the fewest unique genes. However, there are most (21) IS6110 insertion sequences in the strain 1458 genome at either intragenic or intergenic sites, resulting in the interruption of 11 genes including three PPE family-encoding genes (PPE16, PPE38, and PPE59). Only the strain 1458 genome has the upstream insertion in esxS and phoP genes. PCR confirmed four upstream insertions and qPCR determined that transcription of esxS, phoP, dnaN, and ctpD genes differed significantly between M. tb strain 1458 and H37Rv or M. bovis. A Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that the genes affected by non-synonymous SNPs are enriched in RNA polymerase. Moreover, 127 of the 133 unique SNPs in strain 1458 are either different to those in the M. bovis genome. In conclusion, some critical genes responsible for bacterial virulence and immunogenicity were interrupted in the genome of bovine M. tb strain 1458 by IS insertions and non-synonymous SNPs, which might contribute to its bovine adaptation, and the modification of its virulence and immunogenicity in cattle.

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.

A High-Throughput Assay for Developing Inhibitors of PhoP, a Virulence Factor of Mycobacterium tuberculosis

BACKGROUND

Tuberculosis (TB) kills over 1.5 million people per year despite the available anti-TB drugs. The long duration needed to treat TB by the current TB drugs, which target the essential cellular activities, inevitably leads to the emergence of drug-resistance. The emergence of drug-resistant TB prompts for an urgent need for new and more effective drugs.

OBJECTIVE

The response regulator PhoP, an essential virulence factor of Mycobacterium tuberculosis (MTB), is an attractive target for developing novel anti- TB drugs. This study aims to develop a robust high-throughput screening assay to identify PhoP inhibitors that disrupt the PhoP-DNA binding.

METHOD

Guided by the crystal structure of the PhoP-DNA complex, we designed and developed an assay based on Foster resonance energy transfer (FRET) by labeling Cy3 on the DNA and Cy5 on PhoP. We screened compound libraries for inhibitors that dissociated the PhoP-DNA complex by detection of the FRET signal. Hits were confirmed for their direct binding to PhoP by thermal shift assays.

RESULTS

From a test screening of ~6,000 bioactive compounds and approved drugs, three active compounds were identified that directly bound to PhoP and inhibited the PhoP-DNA interactions. These three PhoP inhibitors can be further developed to improve potency and are useful to study the mechanism of inhibition.

CONCLUSION

Our results demonstrated that this FRET-based PhoP-DNA binding assay is valid for additional compound library screening to identify new leads for developing novel TB drugs that target the virulence of MTB.

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.

Two-Component Regulatory Systems of Mycobacteria

Two-component regulatory systems (2CRSs) are widely used by bacteria to sense and respond to environmental stimuli with coordinated changes in gene expression. Systems are normally comprised of a sensory kinase protein that activates a transcriptional regulator by phosphorylation. Mycobacteria have few 2CRSs, but they are of key importance for bacterial survival and play important roles in pathogenicity. Mycobacterium tuberculosis has 12 paired two-component regulatory systems (which include a system with two regulators and one sensor, and a split sensor system), as well as four orphan regulators. Several systems are involved in virulence, and disruption of different systems leads to attenuation or hypervirulence. PhoPR plays a major role in regulating cell wall composition, and its inactivation results in sufficient attenuation of M. tuberculosis that deletion strains are live vaccine candidates. MprAB controls the stress response and is required for persistent infections. SenX3-RegX3 is required for control of aerobic respiration and phosphate uptake, and PrrAB is required for adaptation to intracellular infection. MtrAB is an essential system that controls DNA replication and cell division. The remaining systems (KdpDE, NarL, TrcRS, TcrXY, TcrA, PdtaRS, and four orphan regulators) are less well understood. The structure and binding motifs for several regulators have been characterized, revealing variations in function and operation. The sensors are less well characterized, and stimuli for many remain to be confirmed. This chapter reviews our current understanding of the role of two-component systems in mycobacteria, in particular M. tuberculosis.

Evolutionary landscape of the Mycobacterium tuberculosis complex from the viewpoint of PhoPR: implications for virulence regulation and application to vaccine development

Different members of the Mycobacterium genus have evolved to cause tuberculosis in diverse human populations and in a variety of animal species. Our cumulative knowledge of mycobacterial genomes indicates that mutations in the PhoPR two-component virulence system were acquired not only during the natural evolution of mycobacterial species but also during in vitro subculture, which has given rise to the attenuated reference strain H37Ra or to different daughter strains of Mycobacterium bovis BCG. PhoPR is a well-known regulator of pathogenic phenotypes, including secretion of the virulence factor ESAT-6, biosynthesis of acyltrehalose-based lipids, and modulation of antigen export, in members of the Mycobacterium tuberculosis complex (MTBC). Evolutionarily conserved polymorphisms in PhoPR from Mycobacterium africanum, M. bovis, or M. tuberculosis H37Ra result in loss of functional phenotypes. Interestingly, some members of the MTBC have acquired compensatory mutations to counteract these polymorphisms and, probably, to maintain their pathogenic potential. Some of these compensatory mutations include the insertion of the IS6110 element upstream from phoPR in a particular M. bovis strain that is able to transmit between humans or polymorphisms in M. africanum and M. bovis that affect the regulatory region of the espACD operon, allowing PhoPR-independent ESAT-6 secretion. This review highlights the increasing knowledge of the significance of PhoPR in the evolution of the MTBC and its potential application in the construction of new attenuated vaccines based on phoPR inactivation. In this context, the live attenuated vaccine MTBVAC, based on a phoP fadD26 deletion mutant of M. tuberculosis, is the first vaccine of this kind to successfully enter into clinical development, representing a historic milestone in the field of human vaccinology.

ISMapper: identifying transposase insertion sites in bacterial genomes from short read sequence data

Background

Insertion sequences (IS) are small transposable elements, commonly found in bacterial genomes. Identifying the location of IS in bacterial genomes can be useful for a variety of purposes including epidemiological tracking and predicting antibiotic resistance. However IS are commonly present in multiple copies in a single genome, which complicates genome assembly and the identification of IS insertion sites. Here we present ISMapper, a mapping-based tool for identification of the site and orientation of IS insertions in bacterial genomes, directly from paired-end short read data.

Results

ISMapper was validated using three types of short read data: (i) simulated reads from a variety of species, (ii) Illumina reads from 5 isolates for which finished genome sequences were available for comparison, and (iii) Illumina reads from 7 Acinetobacter baumannii isolates for which predicted IS locations were tested using PCR. A total of 20 genomes, including 13 species and 32 distinct IS, were used for validation. ISMapper correctly identified 97 % of known IS insertions in the analysis of simulated reads, and 98 % in real Illumina reads. Subsampling of real Illumina reads to lower depths indicated ISMapper was able to correctly detect insertions for average genome-wide read depths >20x, although read depths >50x were required to obtain confident calls that were highly-supported by evidence from reads. All ISAba1 insertions identified by ISMapper in the A. baumannii genomes were confirmed by PCR. In each A. baumannii genome, ISMapper successfully identified an IS insertion upstream of the ampC beta-lactamase that could explain phenotypic resistance to third-generation cephalosporins. The utility of ISMapper was further demonstrated by profiling genome-wide IS6110 insertions in 138 publicly available Mycobacterium tuberculosis genomes, revealing lineage-specific insertions and multiple insertion hotspots.

Conclusions

ISMapper provides a rapid and robust method for identifying IS insertion sites directly from short read data, with a high degree of accuracy demonstrated across a wide range of bacteria.

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.

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.

Revealing hidden clonal complexity in Mycobacterium tuberculosis infection by qualitative and quantitative improvement of sampling

The analysis of microevolution events, its functional relevance and impact on molecular epidemiology strategies, constitutes one of the most challenging aspects of the study of clonal complexity in infection by Mycobacterium tuberculosis. In this study, we retrospectively evaluated whether two improved sampling schemes could provide access to the clonal complexity that is undetected by the current standards (analysis of one isolate from one sputum). We evaluated in 48 patients the analysis by mycobacterial interspersed repetitive unit-variable number tandem repeat of M. tuberculosis isolates cultured from bronchial aspirate (BAS) or bronchoalveolar lavage (BAL) and, in another 16 cases, the analysis of a higher number of isolates from independent sputum samples. Analysis of the isolates from BAS/BAL specimens revealed clonal complexity in a very high proportion of cases (5/48); in most of these cases, complexity was not detected when the isolates from sputum samples were analysed. Systematic analysis of isolates from multiple sputum samples also improved the detection of clonal complexity. We found coexisting clonal variants in two of 16 cases that would have gone undetected in the analysis of the isolate from a single sputum specimen. Our results suggest that analysis of isolates from BAS/BAL specimens is highly efficient for recording the true clonal composition of M. tuberculosis in the lungs. When these samples are not available, we recommend increasing the number of isolates from independent sputum specimens, because they might not harbour the same pool of bacteria. Our data suggest that the degree of clonal complexity in tuberculosis has been underestimated because of the deficiencies inherent in a simplified procedure.

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.

Novel GMO-Based Vaccines against Tuberculosis: State of the Art and Biosafety Considerations

Novel efficient vaccines are needed to control tuberculosis (TB), a major cause of morbidity and mortality worldwide. Several TB vaccine candidates are currently in clinical and preclinical development. They fall into two categories, the one of candidates designed as a replacement of the Bacille Calmette Guérin (BCG) to be administered to infants and the one of sub-unit vaccines designed as booster vaccines. The latter are designed as vaccines that will be administered to individuals already vaccinated with BCG (or in the future with a BCG replacement vaccine). In this review we provide up to date information on novel tuberculosis (TB) vaccines in development focusing on the risk assessment of candidates composed of genetically modified organisms (GMO) which are currently evaluated in clinical trials. Indeed, these vaccines administered to volunteers raise biosafety concerns with respect to human health and the environment that need to be assessed and managed.

Current knowledge and pending challenges in zoonosis caused by Mycobacterium bovis: a review

Mycobacterium bovis is both the causative agent of bovine tuberculosis (TB) and a zoonotic pathogen. In humans, considerably fewer cases of TB are caused by M. bovis than M. tuberculosis; nevertheless, diagnostic limitations mean that currently available data on prevalence grossly underestimate the true dimension of the problem. The routes of transmission from animals to humans are well known and include direct exposure to infected animals or consumption of contaminated animal products. Application of fingerprinting tools facilitates analysis of the molecular epidemiology of M. bovis in animal-to-human and human-to-human transmission. Apart from cattle and M. bovis, other animal species and members within the M. tuberculosis complex can contribute to the zoonosis. Improvements in diagnostic techniques, application of more advanced discriminatory genotyping tools, and collaboration between veterinary and human health care researchers are key to our understanding of this zoonosis.

Differences in gene expression between clonal variants of Mycobacterium tuberculosis emerging as a result of microevolution

Clonal variants of Mycobacterium tuberculosis can emerge as a result of microevolution in a single host or after sequential infection of different hosts. The significance of subtle genotypic variations is still unknown. In three of the four loci analyzed from clonal variants differing in only one MIRU-VNTR locus, we found that the expression of the adjacent genes was modulated differently. These data highlight the potential advantages that acquisition of subtle variability may have in M. tuberculosis.

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.