The evolution of mycobacterial pathogenicity: clues from comparative genomics

DiGAlign: Versatile and Interactive Visualization of Sequence Alignment for Comparative Genomics

With the explosion of available genomic information, comparative genomics has become a central approach to understanding microbial ecology and evolution. We developed DiGAlign (https://www.genome.jp/digalign/), a web server that provides versatile functionality for comparative genomics with an intuitive interface. It allows the user to perform the highly customizable visualization of a synteny map by simply uploading nucleotide sequences of interest, ranging from a specific region to the whole genome landscape of microorganisms and viruses. DiGAlign will serve a wide range of biological researchers, particularly experimental biologists, with multifaceted features that allow the rapid characterization of genomic sequences of interest and the generation of a publication-ready figure.

Phylogenetic distribution of malonate semialdehyde decarboxylase (MSAD) genes among strains within the genus Mycobacterium: evidence of MSAD gene loss in the evolution of pathogenic mycobacteria

Despite the great diversity of malonate semialdehyde decarboxylases (MSADs), one of five subgroups of the tautomerase superfamily (TSF) found throughout the biosphere, their distribution among strains within the genus Mycobacterium remains unknown. In this study, we sought to investigate the phylogenetic distribution of MSAD genes of mycobacterial species via genome analysis of 192 different reference Mycobacterium species or subspecies retrieved from NCBI databases. We found that in a total of 87 of 192 strains (45.3%), MSAD-1 and MSAD-2 were distributed in an exclusive manner among Mycobacterium species except for 12 strains, including Mycobacterium chelonae members, with both in their genome. Of note, Mycobacterium strains better adapted to the host and of high virulence potential, such as the Mycobacterium tuberculosis complex, Mycobacterium leprae, Mycobacterium marinum, Mycobacterium ulcerans, and Mycobacterium avium subsp. paratuberculosis, had no orthologs of MSAD in their genome, suggesting MSAD loss during species differentiation in pathogenic slow-growing Mycobacterium. To investigate the MSAD distribution among strains of M. avium subspecies, the genome sequences of a total of 255 reference strains from the four subspecies of M. avium (43 of subspecies avium, 162 of subspecies hominissuis, 49 of subspecies paratuberculosis, and 1 of subspecies silvaticum) were further analyzed. We found that only 121 of 255 strains (47.4%) had MSADs in their genome, with none of the 49 M. avium subsp. paratuberculosis strains having MSAD genes. Even in 13 of 121 M. avium strains with the MSAD-1 gene in their genome, deletion mutations in the 98th codon causing premature termination of MSAD were found, further highlighting the occurrence of MSAD pseudogenization during species or subspecies differentiation of M. avium. In conclusion, our data indicated that there are two distinct types of MSADs, MSAD-1 and MSAD-2, among strains in the Mycobacterium genus, but more than half of the strains, including pathogenic mycobacteria, M. tuberculosis and M. leprae, have no orthologs in their genome, suggesting MSAD loss during host adaptation of pathogenic mycobacteria. In the future, the role of two distinct MSADs, MSAD-1 and MSAD-2, in mycobacterial pathogenesis or evolution should be investigated.

The Antibacterial Type VII Secretion System of Bacillus subtilis: Structure and Interactions of the Pseudokinase YukC/EssB

Type VIIb secretion systems (T7SSb) were recently proposed to mediate different aspects of Firmicutes physiology, including bacterial pathogenicity and competition. However, their architecture and mechanism of action remain largely obscure. Here, we present a detailed analysis of the T7SSb-mediated bacterial competition in Bacillus subtilis, using the effector YxiD as a model for the LXG secreted toxins. By systematically investigating protein-protein interactions, we reveal that the membrane subunit YukC contacts all T7SSb components, including the WXG100 substrate YukE and the LXG effector YxiD. YukC’s crystal structure shows unique features, suggesting an intrinsic flexibility that is required for T7SSb antibacterial activity. Overall, our results shed light on the role and molecular organization of the T7SSb and demonstrate the potential of B. subtilis as a model system for extensive structure-function studies of these secretion machineries.

Comparative Analysis on Proteomics Profiles of Intracellular and Extracellular M.tb and BCG From Infected Human Macrophages

Tuberculosis is the second cause in infectious diseases leading to human death. Understanding the virulence mechanism is inevitable if the disease needs to be fully cured. Therefore, this study aimed to reveal this mechanism by comparing proteomic profiles of intracellular and extracellular virulent strain M.tb and bacille Calmette–Guérin (BCG) from infected THP-1cells. First, M.tb and BCG infected THP-1 at MOI 10:1. Twelve hours postinfection, intracellular bacteria of M.tb and BCG were collected, whereas the two bacilli cultured in 7H9 broth media were used as the control. Then four groups of bacilli were subjected to proteomic analysis, and differential proteomic profiles between M.tb and BCG were comparatively analyzed with bioinformatics tools. As a result, we identified a total of 1,557 proteins. Further, they were divided into four groups for comparison of M.tb versus BCG under 7H9 culture (shorten as out), M.tb in (intracellular) versus M.tb out, BCG in versus BCG out and M.tb in versus BCG in. Between M.tb in versus BCG in, a total of 211 differentially expressed proteins were found. Eight proteins like ESAT-6 distributed in six RDs and some known proteins related to virulence. Besides, five uncharacterized proteins were differentially expressed. Further analysis revealed enriched pathways were associated with glyoxylate and dicarboxylate metabolism pathways. In M.tb out versus BCG out, a total of 144 differential proteins were identified and mainly involved in metabolism pathways. Then, 121 differential proteins in the group of M.tb in versus M.tb out were enriched in ribosome and oxidative phosphorylation related to adaptation to the host environment. The group of BCG in versus BCG out shared the same trend of different pathways to the M.tb in versus M.tb out. Finally, 42 proteins were identified to be up-regulated only in intracellular M.tb including eight RD proteins, whereas 22 up-regulated uniquely in intracellular BCG. Besides, only two proteins (Pks13 and Rv1405c) were commonly up-regulated in intracellular M.tb and BCG. Further, some unknown proteins were uniquely up-regulated in the intracellular M.tb and BCG. These findings provide valuable data for further exploration of molecular mechanism for M.tb virulence and BCG immune response.

Quantitative proteomics reveals that dormancy-related proteins mediate the attenuation in mycobacterium strains

Although members of the Mycobacterium tuberculosis complex (MTBC) exhibit high similarity, they are characterized by differences with respect to virulence, immune response, and transmissibility. To understand the virulence of these bacteria and identify potential novel therapeutic targets, we systemically investigated the total cell protein contents of virulent H37Rv, attenuated H37Ra, and avirulent M. bovis BCG vaccine strains at the log and stationary phases, based on tandem mass tag (TMT) quantitative proteomics. Data analysis revealed that we obtained deep-coverage protein identification and high quantification. Although 272 genetic variations were reported in H37Ra and H37Rv, they showed very little expression difference in log and stationary phase. Quantitative comparison revealed H37Ra and H37Rv had significantly dysregulation in log phase (227) compared with stationary phase (61). While BCG and H37Rv, and BCG and H37Ra showed notable differences in stationary phase (1171 and 1124) with respect to log phase (381 and 414). In the log phase, similar patterns of protein abundance were observed between H37Ra and BCG, whereas a more similar expression pattern was observed between H37Rv and H37Ra in the stationary phase. Bioinformatic analysis revealed that the upregulated proteins detected for H37Rv and H37Ra in log phase were virulence-related factors. In both log and stationary phases, the dysregulated proteins detected for BCG, which have also been identified as M. tuberculosis response proteins under dormancy conditions. We accordingly describe the proteomic profiles of H37Rv, H37Ra, and BCG, which we believe will potentially provide a better understanding of H37Rv pathogenesis, H37Ra attenuation, and BCG immuno protection.

CRISPR-Cas Systems: Prospects for Use in Medicine

CRISPR-Cas systems, widespread in bacteria and archaea, are mainly responsible for adaptive cellular immunity against exogenous DNA (plasmid and phage). However, the latest research shows their involvement in other functions, such as gene expression regulation, DNA repair and virulence. In recent years, they have undergone intensive research as convenient tools for genomic editing, with Cas9 being the most commonly used nuclease. Gene editing may be of interest in biotechnology, medicine (treatment of inherited disorders, cancer, etc.), and in the development of model systems for various genetic diseases. The dCas9 system, based on a modified Cas9 devoid of nuclease activity, called CRISPRi, is widely used to control gene expression in bacteria for new drug biotargets validation and is also promising for therapy of genetic diseases. In addition to direct use for genomic editing in medicine, CRISPR-Cas can also be used in diagnostics, for microorganisms’ genotyping, controlling the spread of drug resistance, or even directly as “smart” antibiotics. This review focuses on the main applications of CRISPR-Cas in medicine, and challenges and perspectives of these approaches.

The Biotin Biosynthetic Pathway in Mycobacterium tuberculosis is a Validated Target for the Development of Antibacterial Agents

Mycobacterium tuberculosis, responsible for Tuberculosis (TB), remains the leading cause of mortality among infectious diseases worldwide from a single infectious agent, with an estimated 1.7 million deaths in 2016. Biotin is an essential cofactor in M. tuberculosis that is required for lipid biosynthesis and gluconeogenesis. M. tuberculosis relies on de novo biotin biosynthesis to obtain this vital cofactor since it cannot scavenge sufficient biotin from a mammalian host. The biotin biosynthetic pathway in M. tuberculosis has been well studied and rigorously genetically validated providing a solid foundation for medicinal chemistry efforts. This review examines the mechanism and structure of the enzymes involved in biotin biosynthesis and ligation, summarizes the reported genetic validation studies of the pathway, and then analyzes the most promising inhibitors and natural products obtained from structure-based drug design and phenotypic screening.

Comparative in-silico proteomic analysis discerns potential granuloma proteins of Yersinia pseudotuberculosis

Yersinia pseudotuberculosis is one of the three pathogenic species of the genus Yersinia. Most studies regarding pathogenesis of Y. pseudotuberculosis are based on the proteins related to Type III secretion system, which is a well-known primary virulence factor in pathogenic Gram-negative bacteria, including Y. pseudotuberculosis. Information related to the factors involved in Y. pseudotuberculosis granuloma formation is scarce. In the present study we have used a computational approach to identify proteins that might be potentially involved in formation of Y. pseudotuberculosis granuloma. A comparative proteome analysis and conserved orthologous protein identification was performed between two different genera of bacteria - Mycobacterium and Yersinia, their only common pathogenic trait being ability to form necrotizing granuloma. Comprehensive analysis of orthologous proteins was performed in proteomes of seven bacterial species. This included M. tuberculosis, M. bovis and M. avium paratuberculosis - the known granuloma forming Mycobacterium species, Y. pestis and Y. frederiksenii - the non-granuloma forming Yersinia species and, Y. enterocolitica - that forms micro-granuloma and, Y. pseudotuberculosis - a prominent granuloma forming Yersinia species. In silico proteome analysis indicated that seven proteins (UniProt id A0A0U1QT64, A0A0U1QTE0, A0A0U1QWK3, A0A0U1R1R0, A0A0U1R1Z2, A0A0U1R2S7, A7FMD4) might play some role in Y. pseudotuberculosis granuloma. Validation of the probable involvement of the seven proposed Y. pseudotuberculosis granuloma proteins was done using transcriptome data analysis and, by mapping on a composite protein-protein interaction map of experimentally proved M. tuberculosis granuloma proteins (RD1 locus proteins, ESAT-6 secretion system proteins and intra-macrophage secreted proteins). Though, additional experiments involving knocking out of each of these seven proteins are required to confirm their role in Y. pseudotuberculosis granuloma our study can serve as a basis for further studies on Y. pseudotuberculosis granuloma.

Accumulation of EBI3 induced by virulent Mycobacterium tuberculosis inhibits apoptosis in murine macrophages

Macrophages are the primary host target cells of Mycobacterium tuberculosis (M. tb). As a subunit of immunoregulatory cytokines IL-27 and IL-35, Epstein–Barr virus-induced gene 3 (EBI3) has typically been explored as the secreted form and assessed in terms of its effects triggered by extracellular EBI3. However, little is known about intracellular EBI3 function. In the current study, we report that EBI3 production by macrophages is elevated in TB patients. We further demonstrate that increased EBI3 accumulates in virulent M. tb-treated murine macrophages. Eukaryotic translation elongation factor 1-alpha 1 (eEF1A1) binds to intracellular EBI3 to reduce Lys48 (K48)-linked ubiquitination of EBI3, leading to EBI3 accumulation. Moreover, the intracellular EBI3 inhibits caspase-3-mediated apoptosis in M. tb-treated macrophages. Herein, we propose a novel mechanism for accumulating intracellular EBI3 and its regulation of macrophage apoptosis in response to virulent M. tb.

Metabolite Profiling: A Tool for the Biochemical Characterisation of Mycobacterium sp

Over the last decades, the prevalence of drug-resistance in Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, has increased. These findings have rekindled interest in elucidating the unique adaptive molecular and biochemistry physiology of Mycobacterium. The use of metabolite profiling independently or in combination with other levels of "-omic" analyses has proven an effective approach to elucidate key physiological/biochemical mechanisms associated with Mtb throughout infection. The following review discusses the use of metabolite profiling in the study of tuberculosis, future approaches, and the technical and logistical limitations of the methodology.

Recovery of Non-tuberculous Mycobacteria from Water is Influenced by Phenotypic Characteristics and Decontamination Methods

Infections related to non-tuberculous mycobacteria (NTM) have recently increased worldwide. The transmission of these microorganisms from the environment has been suggested as the main source for human infections. To elucidate the epidemiological aspects and distribution of these pathogens, many studies have evaluated several decontamination methods and protocols to properly isolate NTM from environmental samples, mainly from water. However, no satisfactory strategy has been found for isolation of most of the NTM species harboring different phenotypic characteristics. Here, we evaluated the susceptibility of 23 NTM strains presenting variable growth rate and pigmentation patterns to eight different methods: oxalic acid (2.5% and 5%), cetylpyridinium chloride (CPC) (0.0025% and 0.005%), sodium hydroxide (NaOH) (2% and 4%), and sodium dodecyl sulfate (SDS) plus NaOH (SDS 1.5%-NaOH 0.5% and SDS 3%-NaOH 1%). It was found that the viability of NTM exposed to different decontamination methods varies according to their phenotypic characteristics and two methods (SDS 1.5% plus NaOH 0.5% and CPC 0.0025%) were necessary for effective isolation of all of the species tested. These findings supply important insights for future studies on the environmental occurrence of mycobacteria and improving the sensibility of traditional strategies.

Pathogenomics of Virulence Traits of Plesiomonas shigelloides That Were Deemed Inconclusive by Traditional Experimental Approaches

One of the major challenges of modern medicine includes the failure of conventional protocols to characterize the pathogenicity of emerging pathogens. This is particularly apparent in the case of Plesiomonas shigelloides. Although a number of infections have been linked to this microorganism, experimental evidence of its virulence factors (VFs), obtained by traditional approaches, is somewhat inconclusive. Hence, it remains unclear whether P. shigelloides is a true or opportunistic one. In the current study, four publicly available whole-genome sequences of P. shigelloides (GN7, NCTC10360, 302-73, and LS1) were profiled using bioinformatics platforms to determine the putative candidate VFs to characterize the bacterial pathogenicity. Overall, 134 unique open reading frames (ORFs) were identified that were homologous or orthologous to virulence genes identified in other pathogens. Of these, 52.24% (70/134) were jointly shared by the strains. The numbers of strain-specific virulence traits were 4 in LS1; 7 in NCTC10360; 10 in 302-73; and 15 in GN7. The pathogenicity islands (PAIs) common to all the strains accounted for 24.07% ORFs. The numbers of PAIs exclusive to each strain were 8 in 302-73; 11 in NCTC10360; 14 in GN7; and 18 in LS1. A PAI encoding Vibrio cholerae ToxR-activated gene d protein was specific to 302-73, GN7, and NCTC10360 strains. Out of 33 antibiotic multi-resistance genes identified, 16 (48.48%) genes were intrinsic to all strains. Further, 17 (22.08%) of 77 antibiotic resistance islands were found in all the strains. Out of 23 identified distinct insertion sequences, 13 were only harbored by strain LS1. The number of intact prophages identified in the strains was 1 in GN7; 2 in 302-73; and 2 in NCTC10360. Further, 1 CRISPR element was identified in LS1; 2 in NCTC10360; and 8 in 302-73. Fifteen (78.95%) of 19 secretion systems and secretion effector variants were identified in all the strains. In conclusion, certain P. shigelloides strains might possess VFs associated with gastroenteritis and extraintestinal infections. However, the role of host factors in the onset of infections should not be undermined.

CD82 hypomethylation is essential for tuberculosis pathogenesis via regulation of RUNX1-Rab5/22

The tumor suppressor gene CD82/KAI1 is a member of the tetraspanin superfamily and organizes various membrane-based processes. Mycobacterium tuberculosis (MTB) persists in host macrophages by interfering with phagolysosome biogenesis and inflammatory responses, but the role of CD82 in controlling the intracellular survival of pathogenic mycobacteria within macrophages remains poorly understood. In this study, we demonstrated that the virulent MTB strain H37Rv (MTB Rv) induced CD82 promoter hypomethylation, resulting in CD82 expression. Targeting of the runt-related transcription factor 1 (RUNX1) by CD82 is essential for phagosome arrest via interacting with Rab5/22. This arrest is required for the intracellular growth of MTB in vitro and in vivo, but not for that of MTB H37Ra (MTB Ra) in macrophages. In addition, knockdown or knockout of CD82 or RUNX1 increased antibacterial host defense via phagolysosome biogenesis, inflammatory cytokine production, and subsequent antimicrobial activity both in vitro and in vivo. Notably, the levels of CD82 and RUNX1 in granulomas were elevated in tuberculosis (TB) patients, indicating that CD82 and RUNX1 have clinical significance in human TB. Our findings identify a previously unrecognized role of CD82 hypomethylation in the regulation of phagosome maturation, enhanced intracellular survival, and the innate host immune response to MTB. Thus, the CD82-RUNX1-Rab5/22 axis may be a previously unrecognized virulence mechanism of MTB pathogenesis.

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

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

The Rate and Spectrum of Spontaneous Mutations in Mycobacterium smegmatis, a Bacterium Naturally Devoid of the Postreplicative Mismatch Repair Pathway

Mycobacterium smegmatis is a bacterium that is naturally devoid of known postreplicative DNA mismatch repair (MMR) homologs, mutS and mutL, providing an opportunity to investigate how the mutation rate and spectrum has evolved in the absence of a highly conserved primary repair pathway. Mutation accumulation experiments of M. smegmatis yielded a base-substitution mutation rate of 5.27 × 10⁻¹⁰ per site per generation, or 0.0036 per genome per generation, which is surprisingly similar to the mutation rate in MMR-functional unicellular organisms. Transitions were found more frequently than transversions, with the A:T→G:C transition rate significantly higher than the G:C→A:T transition rate, opposite to what is observed in most studied bacteria. We also found that the transition-mutation rate of M. smegmatis is significantly lower than that of other naturally MMR-devoid or MMR-knockout organisms. Two possible candidates that could be responsible for maintaining high DNA fidelity in this MMR-deficient organism are the ancestral-like DNA polymerase DnaE1, which contains a highly efficient DNA proofreading histidinol phosphatase (PHP) domain, and/or the existence of a uracil-DNA glycosylase B (UdgB) homolog that might protect the GC-rich M. smegmatis genome against DNA damage arising from oxidation or deamination. Our results suggest that M. smegmatis has a noncanonical Dam (DNA adenine methylase) methylation system, with target motifs differing from those previously reported. The mutation features of M. smegmatis provide further evidence that genomes harbor alternative routes for improving replication fidelity, even in the absence of major repair pathways.

Mycobacterial Infections in Rabbits: From the Wild to the Laboratory

Tuberculous mycobacterial diseases such as leprosy and tuberculosis are ancient diseases that currently continue threatening human health in some countries. Non-tuberculous mycobacterial (NTM) infections cause a series of well-defined pathological entities, as well as some opportunistic diseases that have also increased worldwide, being more common among immunocompromised patients but rising also in immunocompetent individuals. Reports on natural infections by mycobacteria in rabbits are scarce and mainly involve NTM such as Mycobacterium avium subsp. avium in pigmy rabbits in the United States and Mycobacterium avium subsp. paratuberculosis in wild rabbits in Europe. Rabbits have been used as laboratory animals through the years, both to generate immunological reagents and as infection models. Mycobacterial infection models have been developed in this animal species showing different susceptibility patterns to mycobacteria in laboratory conditions. The latent tuberculosis model and the cavitary tuberculosis model have been widely used to elucidate pathogenic mechanisms and to evaluate chemotherapy and vaccination strategies. Rabbits have also been used as bovine paratuberculosis infection models. This review aimed to gather both wildlife and experimental infection data on mycobacteriosis in rabbits to assess their role in the spread of these infections as well as their potential use in the experimental study of mycobacterial pathogenesis and treatment.

Mycobacterial proteomics: analysis of expressed proteomes and post-translational modifications to identify candidate virulence factors

The Mycobacterium tuberculosis bacillus has a number of unique features that make it a particularly effective human pathogen. Although genomic analysis has added to our current understanding of the molecular basis by which M. tuberculosis damages its host, proteomics may be better suited to describe the dynamic interactions between mycobacterial and host systems that underpin this disease. The M. tuberculosis proteome has been investigated using proteomics for over a decade, with increasingly sophisticated mass spectrometry technology and sensitive methods for comparative proteomic profiling. Deeper coverage of the M. tuberculosis proteome has led to the identification of hundreds of putative virulence determinants, as well as an unsurpassed coverage of post-translational modifications. Proteomics is therefore uniquely poised to contribute to our understanding of this pathogen, which may ultimately lead to better management of the disease.

The discovery, function and development of the variable number tandem repeats in different Mycobacterium species

The method of genotyping by variable number tandem repeats (VNTRs) facilitates the epidemiological studies of different Mycobacterium species worldwide. Until now, the VNTR method is not fully understood, for example, its discovery, function and classification. The inconsistent nomenclature and terminology of VNTR is especially confusing. In this review, we first describe in detail the VNTRs in Mycobacterium tuberculosis (M. tuberculosis), as this pathogen resulted in more deaths than any other microbial pathogen as well as for which extensive studies of VNTRs were carried out, and then we outline the recent progress of the VNTR-related epidemiological research in several other Mycobacterium species, such as M. abscessus, M. africanum, M. avium, M. bovis, M. canettii, M. caprae, M. intracellulare, M. leprae, M. marinum, M. microti, M. pinnipedii and M. ulcerans from different countries and regions. This article is aimed mainly at the practical notes of VNTR to help the scientists in better understanding and performing this method.

Genetic markers, genotyping methods & next generation sequencing in Mycobacterium tuberculosis

Molecular epidemiology (ME) is one of the main areas in tuberculosis research which is widely used to study the transmission epidemics and outbreaks of tubercle bacilli. It exploits the presence of various polymorphisms in the genome of the bacteria that can be widely used as genetic markers. Many DNA typing methods apply these genetic markers to differentiate various strains and to study the evolutionary relationships between them. The three widely used genotyping tools to differentiate Mycobacterium tuberculosis strains are IS6110 restriction fragment length polymorphism (RFLP), spacer oligotyping (Spoligotyping), and mycobacterial interspersed repeat units - variable number of tandem repeats (MIRU-VNTR). A new prospect towards ME was introduced with the development of whole genome sequencing (WGS) and the next generation sequencing (NGS) methods, where the entire genome is sequenced that not only helps in pointing out minute differences between the various sequences but also saves time and the cost. NGS is also found to be useful in identifying single nucleotide polymorphisms (SNPs), comparative genomics and also various aspects about transmission dynamics. These techniques enable the identification of mycobacterial strains and also facilitate the study of their phylogenetic and evolutionary traits.

Role of the horizontal gene exchange in evolution of pathogenic Mycobacteria

BACKGROUND

Mycobacterium tuberculosis is one of the most dangerous human pathogens, the causative agent of tuberculosis. While this pathogen is considered as extremely clonal and resistant to horizontal gene exchange, there are many facts supporting the hypothesis that on the early stages of evolution the development of pathogenicity of ancestral Mtb has started with a horizontal acquisition of virulence factors. Episodes of infections caused by non-tuberculosis Mycobacteria reported worldwide may suggest a potential for new pathogens to appear. If so, what is the role of horizontal gene transfer in this process?

RESULTS

Availing of accessibility of complete genomes sequences of multiple pathogenic, conditionally pathogenic and saprophytic Mycobacteria, a genome comparative study was performed to investigate the distribution of genomic islands among bacteria and identify ontological links between these mobile elements. It was shown that the ancient genomic islands from M. tuberculosis still may be rooted to the pool of mobile genetic vectors distributed among Mycobacteria. A frequent exchange of genes was observed between M. marinum and several saprophytic and conditionally pathogenic species. Among them M. avium was the most promiscuous species acquiring genetic materials from diverse origins.

CONCLUSIONS

Recent activation of genetic vectors circulating among Mycobacteria potentially may lead to emergence of new pathogens from environmental and conditionally pathogenic Mycobacteria. The species which require monitoring are M. marinum and M. avium as they eagerly acquire genes from different sources and may become donors of virulence gene cassettes to other micro-organisms.

The Mycobacterium tuberculosis outer membrane channel protein CpnT confers susceptibility to toxic molecules

Mycobacterium tuberculosis, the causative agent of tuberculosis, is protected from toxic solutes by an effective outer membrane permeability barrier. Recently, we showed that the outer membrane channel protein CpnT is required for efficient nutrient uptake by M. tuberculosis and Mycobacterium bovis BCG. In this study, we found that the cpnT mutant of M. bovis BCG is more resistant than the wild type to a large number of drugs and antibiotics, including rifampin, ethambutol, clarithromycin, tetracycline, and ampicillin, by 8- to 32-fold. Furthermore, the cpnT mutant of M. bovis BCG was 100-fold more resistant to nitric oxide, a major bactericidal agent required to control M. tuberculosis infections in mice. Thus, CpnT constitutes the first outer membrane susceptibility factor in slow-growing mycobacteria. The dual functions of CpnT in uptake of nutrients and mediating susceptibility to toxic molecules are reflected in macrophage infection experiments: while loss of CpnT was detrimental for M. bovis BCG in macrophages that enable bacterial replication, presumably due to inadequate nutrient uptake, it conferred a survival advantage in macrophages that mount a strong bactericidal response. Importantly, the cpnT gene showed a significantly higher density of nonsynonymous mutations in drug-resistant clinical M. tuberculosis strains, indicating that CpnT is under selective pressure in human tuberculosis and/or during chemotherapy. Our results indicate that the CpnT channel constitutes an outer membrane gateway controlling the influx of nutrients and toxic molecules into slow-growing mycobacteria. This study revealed that reducing protein-mediated outer membrane permeability might constitute a new drug resistance mechanism in slow-growing mycobacteria.

Whole-genome sequencing for comparative genomics and de novo genome assembly

Next-generation sequencing technologies for whole-genome sequencing of mycobacteria are rapidly becoming an attractive alternative to more traditional sequencing methods. In particular this technology is proving useful for genome-wide identification of mutations in mycobacteria (comparative genomics) as well as for de novo assembly of whole genomes. Next-generation sequencing however generates a vast quantity of data that can only be transformed into a usable and comprehensible form using bioinformatics. Here we describe the methodology one would use to prepare libraries for whole-genome sequencing, and the basic bioinformatics to identify mutations in a genome following Illumina HiSeq or MiSeq sequencing, as well as de novo genome assembly following sequencing using Pacific Biosciences (PacBio).

Comparative functional genomics and the bovine macrophage response to strains of the mycobacterium genus

Mycobacterial infections are major causes of morbidity and mortality in cattle and are also potential zoonotic agents with implications for human health. Despite the implementation of comprehensive animal surveillance programs, many mycobacterial diseases have remained recalcitrant to eradication in several industrialized countries. Two major mycobacterial pathogens of cattle are Mycobacterium bovis and Mycobacterium avium subspecies paratuberculosis (MAP), the causative agents of bovine tuberculosis (BTB) and Johne's disease (JD), respectively. BTB is a chronic, granulomatous disease of the respiratory tract that is spread via aerosol transmission, while JD is a chronic granulomatous disease of the intestines that is transmitted via the fecal-oral route. Although these diseases exhibit differential tissue tropism and distinct complex etiologies, both M. bovis and MAP infect, reside, and replicate in host macrophages - the key host innate immune cell that encounters mycobacterial pathogens after initial exposure and mediates the subsequent immune response. The persistence of M. bovis and MAP in macrophages relies on a diverse series of immunomodulatory mechanisms, including the inhibition of phagosome maturation and apoptosis, generation of cytokine-induced necrosis enabling dissemination of infection through the host, local pathology, and ultimately shedding of the pathogen. Here, we review the bovine macrophage response to infection with M. bovis and MAP. In particular, we describe how recent advances in functional genomics are shedding light on the host macrophage-pathogen interactions that underlie different mycobacterial diseases. To illustrate this, we present new analyses of previously published bovine macrophage transcriptomics data following in vitro infection with virulent M. bovis, the attenuated vaccine strain M. bovis BCG, and MAP, and discuss our findings with respect to the differing etiologies of BTB and JD.

Mycobacterium tuberculosis pili (MTP), a putative biomarker for a tuberculosis diagnostic test

Novel biomarkers are urgently needed for point of care TB diagnostics. In this study, we investigated the potential of the pilin subunit protein encoded by the mtp gene as a diagnostic biomarker. BLAST analysis of the mtp gene on published genome databases, and amplicon sequencing were performed in Mycobacterium tuberculosis Complex (MTBC) strains and other organisms. The protein secondary structure of the amino acid sequences of non-tuberculous Mycobacteria that partially aligned with the mtp sequence was analysed with PredictProtein software. The mtp gene and corresponding amino acid sequence of MTBC were 100% homologous with H37Rv, in contrast to the partial alignment of the non-tuberculous Mycobacteria. The mtp gene was present in all 91 clinical isolates of MTBC. Except for 2 strains with point mutations, the sequence was 100% conserved among the clinical strains. The mtp gene could not be amplified in all non-tuberculous Mycobacteria and respiratory organisms. The predicted MTP protein structure of Mycobacterium avium, Mycobacterium ulcerans and Mycobacterium abscessus differed significantly from that of the M. tuberculosis, which was similar to Mycobacterium marinum. The absence of the mtp gene in non-tuberculous Mycobacteria and other respiratory bacteria suggests that its encoded product, the pilin subunit protein of M. tuberculosis may be a suitable marker for a point of care TB test.

Comparison of the membrane proteome of virulent Mycobacterium tuberculosis and the attenuated Mycobacterium bovis BCG vaccine strain by label-free quantitative proteomics

The Mycobacterium tuberculosis membrane is rich in antigens that are potential targets for diagnostics and the development of new vaccines. To better understand the mechanisms underlying MTB virulence and identify new targets for therapeutic intervention, we investigated the differential composition of membrane proteomes between virulent M. tuberculosis H37Rv (MTB) and the Mycobacterium bovis BCG vaccine strain. To compare the membrane proteomes, we used LC-MS/MS analysis in combination with label-free quantitative proteomics, utilizing the area under the curve of the extracted ion chromatograms of peptides obtained from m/z and retention time alignment of MS1 features. With this approach, we obtained relative abundance ratios for 2203 identified membrane-associated proteins in high confidence. Of these proteins, 294 showed statistically significant differences of at least two fold in relative abundance between MTB and BCG membrane fractions. Our comparative analysis detected several proteins associated with known genomic regions of difference between MTB and BCG as being absent, which validated the accuracy of our approach. In further support of our label-free quantitative data, we verified select protein differences by immunoblotting. To our knowledge, we have generated the first comprehensive and high-coverage profile of comparative membrane proteome changes between virulent MTB and its attenuated relative BCG, which helps elucidate the proteomic basis of the intrinsic virulence of the MTB pathogen.

Genome-level analyses of Mycobacterium bovis lineages reveal the role of SNPs and antisense transcription in differential gene expression

Background

Bovine tuberculosis (bTB) is a disease with major implications for animal welfare and productivity, as well as having the potential for zoonotic transmission. In Great Britain (GB) alone, controlling bTB costs in the region of £100 million annually, with the current control scheme seemingly unable to stop the inexorable spread of infection. One aspect that may be driving the epidemic is evolution of the causative pathogen, Mycobacterium bovis. To understand the underlying genetic changes that may be responsible for this evolution, we performed a comprehensive genome-level analyses of 4 M. bovis strains that encompass the main molecular types of the pathogen circulating in GB.

Results

We have used a combination of genome sequencing, transcriptome analyses, and recombinant DNA technology to define genetic differences across the major M. bovis lineages circulating in GB that may give rise to phenotypic differences of practical importance. The genomes of three M. bovis field isolates were sequenced using Illumina sequencing technology and strain specific differences in gene expression were measured during in vitro growth and in ex vivo bovine alveolar macrophages using a whole genome amplicon microarray and a whole genome tiled oligonucleotide microarray. SNP/small base pair insertion and deletions and gene expression data were overlaid onto the genomic sequence of the fully sequenced strain of M. bovis 2122/97 to link observed strain specific genomic differences with differences in RNA expression.

Conclusions

We show that while these strains show extensive similarities in their genetic make-up and gene expression profiles, they exhibit distinct expression of a subset of genes. We provide genomic, transcriptomic and functional data to show that synonymous point mutations (sSNPs) on the coding strand can lead to the expression of antisense transcripts on the opposing strand, a finding with implications for how we define a 'silent’ nucleotide change. Furthermore, we show that transcriptomic data based solely on amplicon arrays can generate spurious results in terms of gene expression profiles due to hybridisation of antisense transcripts. Overall our data suggest that subtle genetic differences, such as sSNPS, may have important consequences for gene expression and subsequent phenotype.

Replicon-dependent bacterial genome evolution: the case of Sinorhizobium meliloti

Many bacterial species, such as the alphaproteobacterium Sinorhizobium meliloti, are characterized by open pangenomes and contain multipartite genomes consisting of a chromosome and other large-sized replicons, such as chromids, megaplasmids, and plasmids. The evolutionary forces in both functional and structural aspects that shape the pangenome of species with multipartite genomes are still poorly understood. Therefore, we sequenced the genomes of 10 new S. meliloti strains, analyzed with four publicly available additional genomic sequences. Results indicated that the three main replicons present in these strains (a chromosome, a chromid, and a megaplasmid) partly show replicon-specific behaviors related to strain differentiation. In particular, the pSymB chromid was shown to be a hot spot for positively selected genes, and, unexpectedly, genes resident in the pSymB chromid were also found to be more widespread in distant taxa than those located in the other replicons. Moreover, through the exploitation of a DNA proximity network, a series of conserved “DNA backbones” were found to shape the evolution of the genome structure, with the rest of the genome experiencing rearrangements. The presented data allow depicting a scenario where the pSymB chromid has a distinctive role in intraspecies differentiation and in evolution through positive selection, whereas the pSymA megaplasmid mostly contributes to structural fluidity and to the emergence of new functions, indicating a specific evolutionary role for each replicon in the pangenome evolution.

Comparative genomic structures of Mycobacterium CRISPR-Cas

Clustered regularly interspaced short palindromic repeats (CRISPR) are inheritable genetic elements of many archaea and bacteria, conferring acquired immunity against invading nucleic acids. CRISPR might be indicative of the bacterial niche adaptation and evolutionary. Mycobacterium is an important genus occupying diverse niches with profound medical and environmental significance. To present a comparative genomic landscape of the Mycobacterium CRISPR, the feature of mycobacterium CRISPR structures with sequenced complete genomes were bioinformatically analyzed. The results show that CRISPR structures can be found among 14 mycobacteria, and all loci are chromosomally located. Long CRISPRs present in three species, namely M. tuberculosis, M. bovis, and M. avium. Integrated CRISPR-Cas system can only be found in M. tuberculosis and M. bovis, with highly conserved repeat sequences, very short leaders, and promoterless. M. tuberculosis and M. bovis repeat sequences cannot form stable RNA secondary structure, consistent with a Cas6-binding sequence. M. avium repeat sequences can form classical stem-loop structure. A three-step model of M. tuberculosis CRISPR-Cas system action was put forward based on the composition and function of cas genes cluster. M. tuberculosis and M. bovis CRISPRs might interfere with the invading nucleic acids, but have somehow lost the capacity to incorporate new spacers and co-evolve with corresponding mycobacteriophages.

Molecular detection of nontuberculous mycobacteria: advantages and limits of a broad-range sequencing approach

The isolation of nontuberculous mycobacteria (NTM) from clinical specimens has become very common in recent years. Such organisms are typically environmental and occasionally pathogenic for humans and animals. Standard diagnosis of mycobacterial infections relies on direct examination and culture. However, molecular tools are now available which allow quicker and more accurate diagnosis. Detection of NTM can be performed directly from clinical samples, although identification is mostly carried out after isolation. Sequencing of genomic targets (such as 16S rRNA, ITS, rpoB or hsp65) allows accurate and rapid identification, but has some technical limitations. A brief summary of the molecular methods available for NTM identification and a discussion of the problems associated with the use of sequencing analysis together with a description of available algorithms for NTM identification are the major objectives of this review.

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.

Computational genomics-proteomics and Phylogeny analysis of twenty one mycobacterial genomes (Tuberculosis & non Tuberculosis strains)

Background

The genus Mycobacterium comprises different species, among them the most contagious and infectious bacteria. The members of the complex Mycobacterium tuberculosis are the most virulent microorganisms that have killed human and other mammals since millennia. Additionally, with the many different mycobacterial sequences available, there is a crucial need for the visualization and the simplification of their data. In this present study, we aim to highlight a comparative genome, proteome and phylogeny analysis between twenty-one mycobacterial (Tuberculosis and non tuberculosis) strains using a set of computational and bioinformatics tools (Pan and Core genome plotting, BLAST matrix and phylogeny analysis).

Results

Considerably the result of pan and core genome Plotting demonstrated that less than 1250 Mycobacterium gene families are conserved across all species, and a total set of about 20,000 gene families within the Mycobacterium pan-genome of twenty one mycobacterial genomes.

Viewing the BLAST matrix a high similarity was found among the species of the complex Mycobacterium tuberculosis and less conservation is found with other slow growing pathogenic mycobacteria.

Phylogeny analysis based on both protein conservation, as well as rRNA clearly resolve known relationships between slow growing mycobacteria.

Conclusion

Mycobacteria include important pathogenic species for human and animals and the Mycobacterium tuberculosis complex is the most cause of death of the humankind. The comparative genome analysis could provide a new insight for better controlling and preventing these diseases.

Comparative analysis of mycobacterial truncated hemoglobin promoters and the groEL2 promoter in free-living and intracellular mycobacteria

The success of Mycobacterium tuberculosis depends on its ability to withstand and survive the hazardous environment inside the macrophages that are created by reactive oxygen intermediates, reactive nitrogen intermediates, severe hypoxia, low pH, and high CO(2) levels. Therefore, an effective detoxification system is required for the pathogen to persist in vivo. The genome of M. tuberculosis contains a new family of hemoproteins named truncated hemoglobin O (trHbO) and truncated hemoglobin N (trHbN), encoded by the glbO and glbN genes, respectively, important in the survival of M. tuberculosis in macrophages. Mycobacterial heat shock proteins are known to undergo rapid upregulation under stress conditions. The expression profiles of the promoters of these genes were studied by constructing transcriptional fusions with green fluorescent protein and monitoring the promoter activity in both free-living and intracellular milieus at different time points. Whereas glbN showed an early response to the oxidative and nitrosative stresses tested, glbO gave a lasting response to lower concentrations of both stresses. At all time points and under all stress conditions tested, groEL2 showed higher expression than both trHb promoters and expression of both promoters showed an increase while inside the macrophages. Real-time PCR analysis of trHb and groEL2 mRNAs showed an initial upregulation at 24 h postinfection. The presence of the glbO protein imparted an increased survival to M. smegmatis in THP-1 differentiated macrophages compared to that imparted by the glbN and hsp65 proteins. The comparative upregulation shown by both trHb promoters while grown inside macrophages indicates the importance of these promoters for the survival of M. tuberculosis in the hostile environment of the host.

Independent large scale duplications in multiple M. tuberculosis lineages overlapping the same genomic region

Mycobacterium tuberculosis, the causative agent of most human tuberculosis, infects one third of the world's population and kills an estimated 1.7 million people a year. With the world-wide emergence of drug resistance, and the finding of more functional genetic diversity than previously expected, there is a renewed interest in understanding the forces driving genome evolution of this important pathogen. Genetic diversity in M. tuberculosis is dominated by single nucleotide polymorphisms and small scale gene deletion, with little or no evidence for large scale genome rearrangements seen in other bacteria. Recently, a single report described a large scale genome duplication that was suggested to be specific to the Beijing lineage. We report here multiple independent large-scale duplications of the same genomic region of M. tuberculosis detected through whole-genome sequencing. The duplications occur in strains belonging to both M. tuberculosis lineage 2 and 4, and are thus not limited to Beijing strains. The duplications occur in both drug-resistant and drug susceptible strains. The duplicated regions also have substantially different boundaries in different strains, indicating different originating duplication events. We further identify a smaller segmental duplication of a different genomic region of a lab strain of H37Rv. The presence of multiple independent duplications of the same genomic region suggests either instability in this region, a selective advantage conferred by the duplication, or both. The identified duplications suggest that large-scale gene duplication may be more common in M. tuberculosis than previously considered.

Alignment of multiple complete genomes suggests that gene rearrangements may contribute towards the speciation of Mycobacteria

To more accurately define the taxonomic relationships among species belonging to the genus Mycobacterium we have applied and compared three complete genome sequence comparison procedures to existing systems. These included a nucleotide sequence comparison including both coding and no-coding regions of the genome and two genomic-order comparisons using MAUVE and M-GCAT software to provide comparative gene synteny. These methods clearly differentiated a panel of genomes from reference mycobacterial species. Overall, the speciation of bacteria through determination of gene rearrangements were consistent with the gold standard method for species definition in bacteria, DNA-DNA hybridization however within the context of this system, individual components of the Mycobacterium tuberculosis complex (MTBC) did not show sufficient diversity to classify them as a separate species. The high number of gene rearrangements observed between the species tested suggests that gene reorganization of the genome represents an important contributor to speciation within the genus Mycobacterium and other related genera. The absence of rearrangements amongst MTBC supports their consideration as a single genospecies. Some gene rearrangements provided clear internal synteny between genomes of mycobacterial strains belonging to a same species and we suggest these could be used to classify subspecies.

Environmental mycobacteria: a threat to human health?

In many cases, bacterial pathogens are close relatives to nonpathogens. Pathogens seem to be limited lineages within nonpathogenic bacteria. Nonpathogenic isolates are generally more diverse and widespread in the environment and it is generally considered that environmental bacteria do not pose a risk to human health as clinical isolates do; this may not be the case with mycobacteria, but environmental mycobacteria have not been well studied. It is documented that several environmental mycobacteria constitute a source for human infections. Diverse mycobacterial environmental isolates are rarely involved in human disease. Environmental mycobacteria may have a role in degradation of different compounds. Environmental mycobacteria have had a long interaction with humans, maybe as long as the human species, and may have contributed to human evolution.

Proteomic profile of culture filtrate from the Brazilian vaccine strain Mycobacterium bovis BCG Moreau compared to M. bovis BCG Pasteur

Background

Bacille Calmette-Guerin (BCG) is currently the only available vaccine against tuberculosis (TB) and comprises a heterogeneous family of sub-strains with genotypic and phenotypic differences. The World Health Organization (WHO) affirms that the characterization of BCG sub-strains, both on genomic and proteomic levels, is crucial for a better comprehension of the vaccine. In addition, these studies can contribute in the development of a more efficient vaccine against TB. Here, we combine two-dimensional electrophoresis (2DE) and mass spectrometry to analyse the proteomic profile of culture filtrate proteins (CFPs) from M. bovis BCG Moreau, the Brazilian vaccine strain, comparing it to that of BCG Pasteur. CFPs are considered of great importance given their dominant immunogenicity and role in pathogenesis, being available for interaction with host cells since early infection.

Results

The 2DE proteomic map of M. bovis BCG Moreau CFPs in the pH range 3 - 8 allowed the identification of 158 spots corresponding to 101 different proteins, identified by MS/MS. Comparison to BCG Pasteur highlights the great similarity between these BCG strains. However, quantitative analysis shows a higher expression of immunogenic proteins such as Rv1860 (BCG1896, Apa), Rv1926c (BCG1965c, Mpb63) and Rv1886c (BCG1923c, Ag85B) in BCG Moreau when compared to BCG Pasteur, while some heat shock proteins, such as Rv0440 (BCG0479, GroEL2) and Rv0350 (BCG0389, DnaK), show the opposite pattern.

Conclusions

Here we report the detailed 2DE profile of CFPs from M. bovis BCG Moreau and its comparison to BCG Pasteur, identifying differences that may provide relevant information on vaccine efficacy. These findings contribute to the detailed characterization of the Brazilian vaccine strain against TB, revealing aspects that may lead to a better understanding of the factors leading to BCG's variable protective efficacy against TB.

Mycobacterium tuberculosis spoligotypes and drug susceptibility pattern of isolates from tuberculosis patients in South-Western Uganda

Background

Determination of the prevalence and drug susceptibility of the M. tuberculosis strains is important in tuberculosis control. We determined the genetic diversity and susceptibility profiles of mycobacteria isolated from tuberculosis patients in Mbarara, South Western Uganda.

Methods

We enrolled, consecutively; all newly diagnosed and previously treated smear-positive TB patients aged ≥ 18 years. The isolates were characterized using regions of difference (RD) analysis and spoligotyping. Drug resistance against rifampicin and isoniazid were tested using the Genotype^® MDRTBplus assay and the indirect proportion method on Lowenstein-Jensen media. HIV-1 testing was performed using two rapid HIV tests.

Results

A total of 125 isolates from 167 TB suspects (60% males) with a mean age 33.7 years and HIV prevalence of 67.9% (55/81) were analyzed. Majority (92.8%) were new cases while only 7.2% were retreatment cases. All the 125 isolates were identified as M. tuberculosis strict sense with the majority (92.8%) of the isolates being modern strains while seven (7.2%) isolates were ancestral strains. Spoligotyping revealed 79 spoligotype patterns, with an overall diversity of 63.2%. Sixty two (49.6%) of the isolates formed 16 clusters consisting of 2-15 isolates each. A majority (59.2%) of the isolates belong to the Uganda genotype group of strains. The major shared spoligotypes in our sample were SIT 135 (T2-Uganda) with 15 isolates and SIT 128 (T2) with 3 isolates. Sixty nine (87%) of the 79 patterns had not yet been defined in the SpolDB4.0.database. Resistance mutations to either RIF or INH were detected in 6.4% of the isolates. Multidrug resistance, INH and RIF resistance was 1.6%, 3.2% and 4.8%, respectively. The rpoβ gene mutations seen in the sample were D516V, S531L, H526Y H526D and D516V, while one strain had a Δ1 mutation in the wild type probes. There were three strains with katG (codon 315) gene mutations only while one strain showed the inhA promoter gene mutation.

Conclusion

The present study shows that the TB epidemic in Mbarara is caused by modern M. tuberculosis strains mainly belonging to the Uganda genotype and anti-TB drug resistance rate in the region is low.

Modulation of cell death by M. tuberculosis as a strategy for pathogen survival

It has been clearly demonstrated that in vitro, virulent M. tuberculosis can favor necrosis over apoptosis in infected macrophages, and this has been suggested as a mechanism for evading the host immune response. We recently reported that an effect consistent with this hypothesis could be observed in cells from the blood of TB patients, and in this paper, we review what is known about evasion strategies employed by M. tuberculosis and in particular consider the possible interaction of the apoptosis-inhibiting effects of M. tuberculosis infection with another factor (IL-4) whose expression is thought to play a role in the failure to control M. tuberculosis infection. It has been noted that IL-4 may exacerbate TNF-α-induced pathology, though the mechanism remains unexplained. Since pathology in TB typically involves inflammatory aggregates around infected cells, where TNF-α plays an important role, we predicted that IL-4 would inhibit the ability of cells to remove M. tuberculosis by apoptosis of infected cells, through the extrinsic pathway, which is activated by TNF-α. Infection of human monocytic cells with mycobacteria in vitro, in the presence of IL-4, appears to promote necrosis over apoptosis in infected cells—a finding consistent with its suggested role as a factor in pathology during M. tuberculosis infection.

A survey of nucleotide cyclases in actinobacteria: unique domain organization and expansion of the class III cyclase family in Mycobacterium tuberculosis

Cyclic nucleotides are well-known second messengers involved in the regulation of important metabolic pathways or virulence factors. There are six different classes of nucleotide cyclases that can accomplish the task of generating cAMP, and four of these are restricted to the prokaryotes. The role of cAMP has been implicated in the virulence and regulation of secondary metabolites in the phylum Actinobacteria, which contains important pathogens, such as Mycobacterium tuberculosis, M. leprae, M. bovis and Corynebacterium, and industrial organisms from the genus Streptomyces. We have analysed the actinobacterial genome sequences found in current databases for the presence of different classes of nucleotide cyclases, and find that only class III cyclases are present in these organisms. Importantly, prominent members such as M. tuberculosis and M. leprae have 17 and 4 class III cyclases, respectively, encoded in their genomes, some of which display interesting domain fusions seen for the first time. In addition, a pseudogene corresponding to a cyclase from M. avium has been identified as the only cyclase pseudogene in M. tuberculosis and M. bovis. The Corynebacterium and Streptomyces genomes encode only a single adenylyl cyclase each, both of which have corresponding orthologues in M. tuberculosis. A clustering of the cyclase domains in Actinobacteria reveals the presence of typical eukaryote-like, fungi-like and other bacteria-like class III cyclase sequences within this phylum, suggesting that these proteins may have significant roles to play in this important group of organisms.

Comparison of promoter-specific events during transcription initiation in mycobacteria

DNA-protein interactions that occur during transcription initiation play an important role in regulating gene expression. To initiate transcription, RNA polymerase (RNAP) binds to promoters in a sequence-specific fashion. This is followed by a series of steps governed by the equilibrium binding and kinetic rate constants, which in turn determine the overall efficiency of the transcription process. We present here the first detailed kinetic analysis of promoter-RNAP interactions during transcription initiation in the sigma(A)-dependent promoters P(rrnAPCL1), P(rrnB) and P(gyr) of Mycobacterium smegmatis. The promoters show comparable equilibrium binding affinity but differ significantly in open complex formation, kinetics of isomerization and promoter clearance. Furthermore, the two rrn promoters exhibit varied kinetic properties during transcription initiation and appear to be subjected to different modes of regulation. In addition to distinct kinetic patterns, each one of the housekeeping promoters studied has its own rate-limiting step in the initiation pathway, indicating the differences in their regulation.

Genetics of tuberculosis in Irish Holstein-Friesian dairy herds

Information is lacking on genetic parameters for tuberculosis (TB) susceptibility in dairy cattle. Mycobacterium bovis is the principal agent of tuberculosis in cattle. The objective of this study was to quantify the genetic variation present among Irish Holstein-Friesian dairy herds in their susceptibility to M. bovis infection. A total of 15,182 cow and 8,104 heifer single intradermal comparative tuberculin test (SICTT, a test for M. bovis exposure and presumed infection) records from November 1, 2002, to October 31, 2005, were available for inclusion in the analysis. Data on observed carcass TB lesions from abattoirs were also available for inclusion in the analysis. The only animals retained were those present in a herd during episodes in which at least 2 animals showed evidence of infection; this ensured a high likelihood of exposure to M. bovis. Linear animal models, and sire and animal threshold models were used to estimate the variance components for susceptibility to M. bovis-purified protein derivative (PPD) responsiveness and confirmed M. bovis infection. The heritability estimates from the threshold sire models were biased upward because the relatedness between dam-daughter pairs was ignored. The threshold animal model produced heritability estimates of 0.14 in cows and 0.12 in heifers for susceptibility to M. bovis-PPD responsiveness, and 0.18 in cows for confirmed M. bovis infection susceptibility. Therefore, exploitable genetic variation exists among Irish dairy cows for susceptibility to M. bovis infection. Sire rankings from the linear and threshold animal models were similar, indicating that either model could be used for the analysis of susceptibility to M. bovis-PPD responsiveness. A favorable genetic correlation close to unity was observed between susceptibility to confirmed M. bovis infection and M. bovis-PPD responsiveness, indicating that direct selection for resistance to M. bovis-PPD responsiveness will indirectly reduce susceptibility to confirmed M. bovis infection. Data from the national TB eradication program could be used routinely to estimate breeding values for susceptibility to M. bovis infection.

Complete genome sequence of Francisella tularensis subspecies holarctica FTNF002-00

Francisella tularensis subspecies holarctica FTNF002-00 strain was originally obtained from the first known clinical case of bacteremic F. tularensis pneumonia in Southern Europe isolated from an immunocompetent individual. The FTNF002-00 complete genome contains the RD₂₃ deletion and represents a type strain for a clonal population from the first epidemic tularemia outbreak in Spain between 1997–1998. Here, we present the complete sequence analysis of the FTNF002-00 genome. The complete genome sequence of FTNF002-00 revealed several large as well as small genomic differences with respect to two other published complete genome sequences of F. tularensis subsp. holarctica strains, LVS and OSU18. The FTNF002-00 genome shares >99.9% sequence similarity with LVS and OSU18, and is also ∼5 MB smaller by comparison. The overall organization of the FTNF002-00 genome is remarkably identical to those of LVS and OSU18, except for a single 3.9 kb inversion in FTNF002-00. Twelve regions of difference ranging from 0.1–1.5 kb and forty-two small insertions and deletions were identified in a comparative analysis of FTNF002-00, LVS, and OSU18 genomes. Two small deletions appear to inactivate two genes in FTNF002-00 causing them to become pseudogenes; the intact genes encode a protein of unknown function and a drug:H⁺ antiporter. In addition, we identified ninety-nine proteins in FTNF002-00 containing amino acid mutations compared to LVS and OSU18. Several non-conserved amino acid replacements were identified, one of which occurs in the virulence-associated intracellular growth locus subunit D protein. Many of these changes in FTNF002-00 are likely the consequence of direct selection that increases the fitness of this subsp. holarctica clone within its endemic population. Our complete genome sequence analyses lay the foundation for experimental testing of these possibilities.

Single nucleotide polymorphisms on the road to strain differentiation in Mycobacterium ulcerans

The genomic fine-typing of strains of Mycobacterium ulcerans, the causative agent of the emerging human disease Buruli ulcer, is difficult due to the clonal population structure of geographical lineages. Although large sequence polymorphisms (LSPs) resulted in the clustering of patient isolates originating from across the globe, differentiation of strains within continents using conventional typing methods is very limited. In this study, we analyzed M. ulcerans LSP haplotype-specific insertion sequence elements among 83 M. ulcerans strains and identified single nucleotide polymorphisms (SNPs) that differentiate between regional strains. This is the first genetic discrimination based on SNPs of M. ulcerans strains from African countries where Buruli ulcer is endemic, resulting in the highest geographic resolution of genotyping so far. The findings support the concept of genome-wide SNP analyses as tools to study the epidemiology and evolution of M. ulcerans at a local level.