Mycobacterium Tuberculosis Genome Research Articles

TB drug susceptibility is more than MIC.

Genome-wide association studies and genetic analyses have identified a clinically prevalent alteration in the Mycobacterium tuberculosis genome that rewires bacterial propionate metabolism, conditionally reduces antibiotic killing without affecting minimum inhibitory concentrations, and may drive emergence of drug resistance.

Recombination Signal in Mycobacterium tuberculosis Stems from Reference-guided Assemblies and Alignment Artefacts.

DNA acquisition via genetic recombination is considered advantageous as it has the potential to bring together beneficial mutations that emerge independently within a population. Furthermore, recombination is considered to contribute to the maintenance of genome stability by purging slightly deleterious mutations. The prevalence of recombination differs among prokaryotic species and depends on the accessibility of DNA transfer mechanisms. An exceptional example is the human pathogen Mycobacterium tuberculosis (MTB) where no clear transfer mechanisms have been so far characterized and the presence of recombination is questioned. Here, we analyze completely assembled MTB genomes in search for evidence of recombination. We find that putative recombination events are enriched in strains reconstructed by reference-guided assembly and in regions with unreliable alignments. In addition, assembly and alignment artefacts introduce phylogenetic signals that are conflicting the established MTB phylogeny. Our results reveal that the so far reported recombination events in MTB are likely to stem from methodological artefacts. We conclude that no reliable signal of recombination is observed in the currently available MTB genomes. Moreover, our study demonstrates the limitations of reference-guided genome assembly for phylogenetic reconstructions. Rigorously de novo assembled genomes of high quality are mandatory in order to distinguish true evolutionary signal from noise, in particular for low diversity species such as MTB.

Comparing Two Mycobacterium tuberculosis Genomes from Chinese Immigrants with Native Genomes Using Mauve Alignments

Comparing Two Mycobacterium tuberculosis Genomes from Chinese Immigrants with Native Genomes Using Mauve Alignments

The PGRS Domain of Mycobacterium tuberculosis PE_PGRS Protein Rv0297 Is Involved in Endoplasmic Reticulum Stress-Mediated Apoptosis through Toll-Like Receptor 4.

ABSTRACTThe genome of Mycobacterium tuberculosis, the causal organism of tuberculosis (TB), encodes a unique protein family known as the PE/PPE/PGRS family, present exclusively in the genus Mycobacterium and nowhere else in the living kingdom, with largely unexplored functions. We describe the functional significance of the PGRS domain of Rv0297, a member of this family. In silico analyses revealed the presence of intrinsically disordered stretches and putative endoplasmic reticulum (ER) localization signals in the PGRS domain of Rv0297 (Rv0297PGRS). The PGRS domain aids in ER localization, which was shown by infecting macrophage cells with M. tuberculosis and by overexpressing the protein by transfection in macrophage cells followed by activation of the unfolded protein response, as evident from increased expression of GRP78/GRP94 and CHOP/ATF4, leading to disruption of intracellular Ca2+ homeostasis and increased nitric oxide (NO) and reactive oxygen species (ROS) production. The consequent activation of the effector caspase-8 resulted in apoptosis of macrophages, which was Toll-like receptor 4 (TLR4) dependent. Administration of recombinant Rv0297PGRS (rRv0297PGRS) also exhibited similar effects. These results implicate a hitherto-unknown role of the PGRS domain of the PE_PGRS protein family in ER stress-mediated cell death through TLR4. Since this protein is already known to be present at later stages of infection in human granulomas it points to the possibility of it being employed by M. tuberculosis for its dissemination via an apoptotic mechanism.

An integrated computational approach of molecular dynamics simulations, receptor binding studies and pharmacophore mapping analysis in search of potent inhibitors against tuberculosis

Tuberculosis is an infectious chronic disease caused by obligate pathogen Mycobacterium tuberculosis that affects millions of people worldwide. Although many first and second line drugs are available for its treatment, but their irrational use has adversely lead to the emerging cases of multiple drug resistant and extensively drug-resistant tuberculosis. Therefore, there is an intense need to develop novel potent analogues for its treatment. This has prompted us to develop potent analogues against TB. The Mycobacterium tuberculosis genome provides us with number of validated targets to combat against TB. Study of Mtb genome disclosed six epoxide hydrolases (A to F) which convert harmful epoxide into diols and act as a potential drug target for rational drug design. Our current strategy is to develop such analogues which inhibits epoxide hydrolase enzyme present in Mtb genome. To achieve this, we adopted an integrated computational approach involving QSAR, pharmacophore mapping, molecular docking and molecular dynamics simulation studies. The approach envisaged vital information about the role of molecular descriptors, essential pharmacophoric features and binding energy for compounds to bind into the active site of epoxide hydrolase. Molecular docking analysis revealed that analogues exhibited significant binding to Mtb epoxide hydrolase. Further, three docked complexes 2s, 37s and 15s with high, moderate and low docking scores respectively were selected for molecular dynamics simulation studies. RMSD analysis revealed that all complexes are stable with average RMSD below 2 Å throughout the 10 ns simulations. The B-factor analysis showed that the active site residues of epoxide hydrolase are flexible enough to interact with inhibitor. Moreover, to confirm the binding of these urea derivatives, MM-GBSA binding energy analysis were performed. The calculations showed that 37s has more binding affinity (ΔGtotal = −52.24 kcal/mol) towards epoxide hydrolase compared to 2s (ΔGtotal = −51.70 kcal/mol) and 15s (ΔGtotal = −49.97 kcal/mol). The structural features inferred in our study may provide the future directions to the scientists towards the discovery of new chemical entity exhibiting anti-TB property.

System-Wide Analysis Unravels the Differential Regulation and In Vivo Essentiality of Virulence-Associated Proteins B and C Toxin-Antitoxin Systems of Mycobacterium tuberculosis

Toxin-antitoxin (TA) systems are bicistronic genetic modules that are ubiquitously present in bacterial genomes. The Mycobacterium tuberculosis genome encodes 90 putative TA systems, and these are considered to be associated with maintenance of bacterial genomic stability or bacterial survival under unfavorable environmental conditions. The majority of these in M. tuberculosis have been annotated as belonging to the virulence-associated protein B and C (VapBC) family. However, their precise role in bacterial physiology has not been elucidated. Here, we functionally characterized VapC toxins from M. tuberculosis and show that overexpression of some homologs inhibits growth of Mycobacterium bovis bacillus Calmette-Guérin in a bacteriostatic manner. Expression profiling of messenger RNA revealed that these VapC toxins were differentially induced upon exposure of M. tuberculosis to stress conditions. We also unraveled that transcriptional cross-activation exists between TA systems in M. tuberculosis. This study provides the first evidence for the essentiality of VapBC3 and VapBC4 systems in M. tuberculosis virulence.

Outcome of polymerase chain reaction (PCR) analysis in 100 suspected cases of infectious uveitis

BackgroundPolymerase chain reaction (PCR) analysis is an important tool in the diagnosis of infectious uveitis. A retrospective, interventional study of PCR analysis of ocular fluid in suspected infectious uveitis cases between January 2014 to July 2016 was done. Nested, real-time and broad range PCR was performed for detection of the genome of Mycobacterium tuberculosis, herpes virus family, Chikungunya virus, Toxoplasma gondii, fungus, eubacterium and propionibacterium acne.ResultsTotal of 100 cases included, mean age was 39.2 ± 15.4 years. Uveitis was unilateral in 82% and granulomatous in 40%. Mean visual acuity at the initial visit and final visit was 0.73 logMar and 0.63 logMar respectively. PCR analysis confirmed the clinical diagnosis in 70.1% patients. The sensitivity, specificity, positive predictive value and negative predictive value of PCR analysis was 90.2%, 93.9%, 93.9% and 90.2% respectively. The quantitative value of real-time M. tb. Positive PCR ranged from 32c/ml to 2722 c/ml.ConclusionsPCR assay is an accurate technique with high sensitivity and specificity to diagnose the DNA genome in infectious uveitis.

Comparing Two Mycobacterium tuberculosis Genomes from Chinese Immigrants with Native Genomes Using Mauve Alignments.

BackgroundThe number of immigrants with tuberculosis (TB) increases each year in South Korea. Determining the transmission dynamics based on whole genome sequencing (WGS) to cluster the strains has been challenging.MethodsWGS, annotation refinement, and orthology assignment for the GenBank accession number acquisition were performed on two clinical isolates from Chinese immigrants. In addition, the genomes of the two isolates were compared with the genomes of Mycobacterium tuberculosis isolates, from two native Korean and five native Chinese individuals using a phylogenetic topology tree based on the Multiple Alignment of Conserved Genomic Sequence with Rearrangements (Mauve) package.ResultsThe newly assigned accession numbers for two clinical isolates were CP020381.2 (a Korean-Chinese from Yanbian Province) and CP022014.1 (a Chinese from Shandong Province), respectively. Mauve alignment classified all nine TB isolates into a discriminative collinear set with matched regions. The phylogenetic analysis revealed a rooted phylogenetic tree grouping the nine strains into two lineages: strains from Chinese individuals and strains from Korean individuals.ConclusionPhylogenetic trees based on the Mauve alignments were supposed to be useful in revealing the dynamics of TB transmission from immigrants in South Korea, which can provide valuable information for scaling up the TB screening policy for immigrants.

Wavelet-domain elastic net for clustering on genomes strains.

We propose to evaluate genome similarity by combining discrete non-decimated wavelet transform (NDWT) and elastic net. The wavelets represent a signal with levels of detail, that is, hidden components are detected by means of the decomposition of this signal, where each level provides a different characteristic. The main feature of the elastic net is the grouping of correlated variables where the number of predictors is greater than the number of observations. The combination of these two methodologies applied in the clustering analysis of the Mycobacterium tuberculosis genome strains proved very effective, being able to identify clusters at each level of decomposition.

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.

Mycobacterium tuberculosis Rv3651 is a triple sensor-domain protein.

The genome of the human pathogen Mycobacterium tuberculosis (Mtb) encodes ∼4,400 proteins, but one third of them have unknown functions. We solved the crystal structure of Rv3651, a hypothetical protein with no discernible similarity to proteins with known function. Rv3651 has a three-domain architecture that combines one cGMP-specific phosphodiesterases, adenylyl cyclases and FhlA (GAF) domain and two Per-ARNT-Sim (PAS) domains. GAF and PAS domains are sensor domains that are typically linked to signaling effector molecules. Unlike these sensor-effector proteins, Rv3651 is an unusual sensor domain-only protein with highly divergent sequence. The structure suggests that Rv3651 integrates multiple different signals and serves as a scaffold to facilitate signal transfer.

Complete Genome Sequencing of Mycobacterium bovis SP38 and Comparative Genomics of Mycobacterium bovis and M. tuberculosis Strains.

Mycobacterium bovis causes bovine tuberculosis and is the main organism responsible for zoonotic tuberculosis in humans. We performed the sequencing, assembly and annotation of a Brazilian strain of M. bovis named SP38, and performed comparative genomics of M. bovis genomes deposited in GenBank. M. bovis SP38 has a traditional tuberculous mycobacterium genome of 4,347,648 bp, with 65.5% GC, and 4,216 genes. The majority of CDSs (2,805, 69.3%) have predictive function, while 1,206 (30.07%) are hypothetical. For comparative analysis, 31 M. bovis, 32 M. bovis BCG, and 23 Mycobacterium tuberculosis genomes available in GenBank were selected. M. bovis RDs (regions of difference) and Clonal Complexes (CC) were identified in silico. Genome dynamics of bacterial groups were analyzed by gene orthology and polymorphic sites identification. M. bovis polymorphic sites were used to construct a phylogenetic tree. Our RD analyses resulted in the exclusion of three genomes, mistakenly annotated as virulent M. bovis. M. bovis SP38 along with strain 35 represent the first report of CC European 2 in Brazil, whereas two other M. bovis strains failed to be classified within current CC. Results of M. bovis orthologous genes analysis suggest a process of genome remodeling through genomic decay and gene duplication. Quantification, pairwise comparisons and distribution analyses of polymorphic sites demonstrate greater genetic variability of M. tuberculosis when compared to M. bovis and M. bovis BCG (p ≤ 0.05), indicating that currently defined M. tuberculosis lineages are more genetically diverse than M. bovis CC and animal-adapted MTC (M. tuberculosis Complex) species. As expected, polymorphic sites annotation shows that M. bovis BCG are subjected to different evolutionary pressures when compared to virulent mycobacteria. Lastly, M. bovis phylogeny indicates that polymorphic sites may be used as markers of M. bovis lineages in association with CC. Our findings highlight the need to better understand host-pathogen co-evolution in genetically homogeneous and/or diverse host populations, considering the fact that M. bovis has a broader host range when compared to M. tuberculosis. Also, the identification of M. bovis genomes not classified within CC indicates that the diversity of M. bovis lineages may be larger than previously thought or that current classification should be reviewed.

METHODS OF DIAGNOSIS OF LATENT TUBERCULOSIS

The decoding of Mycobacterium tuberculosis genome opened new possibilities for the development of new diagnostic methods based on detection of antigen-specific immune response. New methods detecting the immune response to proteins which are mostly produced by actively propogating pathogenic microbacteriummake it possible to assess the risk for development of active tuberculosis in patients with latent tuberculosis as well as to prescribe preventive treatment. Their results may be used as a criterion for complex assessment of active tuberculosis and efficiency of anti-tuberculosis therapy.

Comparative lipidomics of drug sensitive and resistant Mycobacterium tuberculosis reveals altered lipid imprints.

Lipids are most adaptable molecules that acclimatize to the development of multidrug resistance (MDR). The precise molecular mechanism of this acclimatization achieved in Mycobacterium tuberculosis (MTB) remains elusive. Although lipids of MTB have been characterized to some details, a comparable resource does not exist between drug sensitive (DS) and resistant (DR) strains of MTB. Here, by employing high-throughput mass spectrometry-based lipidomic approach, we attempted to analyze the differential lipidome profile of DS and DR MTB clinical isolates. We analyzed three major classes of lipids viz fatty acyls, glycerophospholipids and glycerolipids and their respective subclasses. Notably, we observed differential fatty acyls and glycerophospholipids as evident from increased mycolic acids phosphatidylinositol mannosides, phosphatidylinositol, cardiolipin and triacylglycerides abundance, respectively, which are crucial for MTB virulence and pathogenicity. Considering the fact that 30% of the MTB genome codes for lipid, this comprehensive lipidomic approach unravels extensive lipid alterations in DS and DR that will serve as a resource for identifying biomarkers aimed at disrupting the functions of MTB lipids responsible for MDR acquisition in MTB.

Recombinant expression and biochemical characterization of Mycobacterium tuberculosis 3Fe-4S ferredoxin Rv1786.

Ferredoxins are iron-sulfur protein that mediate electron transfer in cytochrome P450 mono-oxygenase (CYP)-related catalytic reactions in a wide variety of organisms. Rv1786 is a putative ferredoxin, encoded by a gene located downstream of the gene encoding CYP143A1 in the Mycobacterium tuberculosis genome. However, the structure and function of Rv1786 have remained unclear. Here, the recombinant Mtb Rv1786 was expressed, purified as a His-tagged form and characterized with [3Fe-4S] clusters as its cofactors using a series of measurements including SDS-PAGE, western blot, UV/Visible, MALDI-TOF/TOF-MS, and electron paramagnetic resonance spectroscopic analysis. Based on the assessments of surface plasmon resonance (SPR) and steady state kinetic assays, Rv1786 was found to be able to couple with both ferredoxin reductase A (FdrA) and flavoprotein reductase A (FprA) as redox partner, but with a stronger binding to FprA and a better coupling activity to FdrA. Preliminary structural and biochemical characterization of Mtb Rv1786 as a redox partner is presented here.

Complete auxotrophy for unsaturated fatty acids requires deletion of two sets of genes in Mycobacterium smegmatis.

The synthesis of unsaturated fatty acids in Mycobacterium smegmatis is poorly characterized. Bioinformatic analysis revealed four putative fatty acid desaturases in its genome, one of which, MSMEG_1886, is highly homologous to desA3, the only palmitoyl/stearoyl desaturase present in the Mycobacterium tuberculosis genome. A MSMEG_1886 deletion mutant was partially auxotrophic for oleic acid and viable at 37°C and 25°C, although with a long lag phase in liquid medium. Fatty acid analysis suggested that MSMEG_1886 is a palmitoyl/stearoyl desaturase, as the synthesis of palmitoleic acid was abrogated, while oleic acid contents dropped by half in the mutant. Deletion of the operon MSMEG_1741-1743 (highly homologous to a Pseudomonas aeruginosa acyl-CoA desaturase) had little effect on growth of the parental strain; however the double mutant MSMEG_1886-MSMEG_1741-1743 strictly required oleic acid for growth. The ΔMSMEG_1886-ΔMSMEG_1741 double mutant was able to grow (poorly but better than the ΔMSMEG_1886 single mutant) in solid and liquid media devoid of oleic acid, suggesting a repressor role for ΔMSMEG_1741. Fatty acid analysis of the described mutants suggested that MSMEG_1742-43 desaturates C18:0 and C24:0 fatty acids. Thus, although the M. smegmatis desA3 homologue is the major player in unsaturated fatty acid synthesis, a second set of genes is also involved.

A Genomic Signature for Genotyping Mycobacterium tuberculosis.

Mycobacterium tuberculosis (MTB), the causative agent of tuberculosis (TB), has a vast diversity of genotypes including Beijing, CAS, EAI, Haarlem, LAM, X, Ural, T, AFRI1 and AFRI2. However, genotyping can be expensive, time consuming and in some cases, results may vary depending on methodology used. Here, we proposed a new set of 10 SNPs using a total of 249 MTB genomes, and selected by first the inclusion/ exclusion (IE) criteria using spoligotyping and phylogenies, followed by the selection of the nonsynonymous SNPs present in the most conserved cluster of orthologous groups (COG) of each genotype of MTB. Genotype assignment of the new set of 10 SNPs was validated using an additional of 34 MTB genomes and results showed 100% correlation with their known genotypes. Our set of 10 SNPs have not been previously reported and cover the MTB genotypes that are prevalent worldwide. This set of SNPs could be used for molecular epidemiology with drug resistant markers.

Mapping and characterization of G-quadruplexes in Mycobacterium tuberculosis gene promoter regions

Mycobacterium tuberculosis is the causative agent of tuberculosis (TB), one of the top 10 causes of death worldwide in 2015. The recent emergence of strains resistant to all current drugs urges the development of compounds with new mechanisms of action. G-quadruplexes are nucleic acids secondary structures that may form in G-rich regions to epigenetically regulate cellular functions. Here we implemented a computational tool to scan the presence of putative G-quadruplex forming sequences in the genome of Mycobacterium tuberculosis and analyse their association to transcription start sites. We found that the most stable G-quadruplexes were in the promoter region of genes belonging to definite functional categories. Actual G-quadruplex folding of four selected sequences was assessed by biophysical and biomolecular techniques: all molecules formed stable G-quadruplexes, which were further stabilized by two G-quadruplex ligands. These compounds inhibited Mycobacterium tuberculosis growth with minimal inhibitory concentrations in the low micromolar range. These data support formation of Mycobacterium tuberculosis G-quadruplexes in vivo and their potential regulation of gene transcription, and prompt the use of G4 ligands to develop original antitubercular agents.

Whole Genome Sequencing of Mycobacterium tuberculosis Isolates From Extrapulmonary Sites.

Tuberculosis (TB) remains one of the leading causes of morbidity and mortality worldwide. Extrapulmonary tuberculosis (EPTB) constitutes around 15-20% of TB cases in immunocompetent individuals. Extrapulmonary sites that are affected by TB include bones, lymph nodes, meningitis, pleura, and genitourinary tract. Whole genome sequencing has emerged as a powerful tool to map genetic diversity among Mycobacterium tuberculosis (MTB) isolates and identify the genomic signatures associated with drug resistance, pathogenesis, and disease transmission. Several pulmonary isolates of MTB have been sequenced over the years. However, availability of whole genome sequences of MTB isolates from extrapulmonary sites is limited. Some studies suggest that genetic variations in MTB might contribute to disease presentation in extrapulmonary sites. This can be addressed if whole genome sequence data from large number of extrapulmonary isolates becomes available. In this study, we have performed whole genome sequencing of five MTB clinical isolates derived from EPTB sites using next-generation sequencing platform. We identified 1434 nonsynonymous single nucleotide variations (SNVs), 143 insertions and 105 deletions. This includes 279 SNVs that were not reported before in publicly available datasets. We found several mutations that are known to confer resistance to drugs. All the five isolates belonged to East-African-Indian lineage (lineage 3). We identified 9 putative prophage DNA integrations and 14 predicted clustered regularly interspaced short palindromic repeats (CRISPR) in MTB genome. Our analysis indicates that more work is needed to map the genetic diversity of MTB. Whole genome sequencing in conjunction with comprehensive drug susceptibility testing can reveal clinically relevant mutations associated with drug resistance.

The Enigmatic PE/PPE Multigene Family of Mycobacteria and Tuberculosis Vaccination.

ABSTRACTThe genome of Mycobacterium tuberculosis, the bacterium responsible for the disease tuberculosis, contains an unusual family of abundant antigens (PE/PPEs). To date, certain members of this multigene family occur only in mycobacteria that cause disease. It is possible that the numerous proteins encoded by these mycobacterial genes dictate the immune pathogenesis of this bacterial pathogen. There is also evidence that some of these antigens are present at the cell surface and that they affect the pathology and immunology of the organism in many ways. Also, they elicit both antibodies and T cells, they may be involved in antigenic variation, and they may be good candidates for vaccines and drugs. However, since they are plentiful and extremely homologous, these PE/PPEs are very challenging to study, and it is difficult to be certain what role(s) they have in the pathogenesis of tuberculosis. Consequently, how to develop treatments like vaccines using these antigens as candidates is complex.