Comparative genomic and proteomic analyses of PE/PPE multigene family of Mycobacterium tuberculosis H₃₇Rv and H₃₇Ra reveal novel and interesting differences with implications in virulence

PE/PPE mutations in the transmission of Mycobacterium tuberculosis in China revealed by whole genome sequencing

OBJECTIVE

This study aims to examine the impact of PE/PPE gene mutations on the transmission of Mycobacterium tuberculosis (M. tuberculosis) in China.

METHODS

We collected the whole genome sequencing (WGS) data of 3202 M. tuberculosis isolates in China from 2007 to 2018 and investigated the clustering of strains from different lineages. To evaluate the potential role of PE/PPE gene mutations in the dissemination of the pathogen, we employed homoplastic analysis to detect homoplastic single nucleotide polymorphisms (SNPs) within these gene regions. Subsequently, logistic regression analysis was conducted to analyze the statistical association.

RESULTS

Based on nationwide M. tuberculosis WGS data, it has been observed that the majority of the M. tuberculosis burden in China is caused by lineage 2 strains, followed by lineage 4. Lineage 2 exhibited a higher number of transmission clusters, totaling 446 clusters, of which 77 were cross-regional clusters. Conversely, there were only 52 transmission clusters in lineage 4, of which 9 were cross-regional clusters. In the analysis of lineage 2 isolates, regression results showed that 4 specific gene mutations, PE4 (position 190,394; c.46G > A), PE_PGRS10 (839,194; c.744 A > G), PE16 (1,607,005; c.620T > G) and PE_PGRS44 (2,921,883; c.333 C > A), were significantly associated with the transmission of M. tuberculosis. Mutations of PE_PGRS10 (839,334; c.884 A > G), PE_PGRS11 (847,613; c.1455G > C), PE_PGRS47 (3,054,724; c.811 A > G) and PPE66 (4,189,930; c.303G > C) exhibited significant associations with the cross-regional clusters. A total of 13 mutation positions showed a positive correlation with clustering size, indicating a positive association. For lineage 4 strains, no mutations were found to enhance transmission, but 2 mutation sites were identified as risk factors for cross-regional clusters. These included PE_PGRS4 (338,100; c.974 A > G) and PPE13 (976,897; c.1307 A > C).

CONCLUSION

Our results indicate that some PE/PPE gene mutations can increase the risk of M. tuberculosis transmission, which might provide a basis for controlling the spread of tuberculosis.

Expression of a unique M. tuberculosis DNA MTase Rv1509 in M. smegmatis alters the gene expression pattern and enhances virulence

Mycobacterium tuberculosis (M. tb) genome encompasses 4,173 genes, about a quarter of which remain uncharacterized and hypothetical. Considering the current limitations associated with the diagnosis and treatment of tuberculosis, it is imperative to comprehend the pathomechanism of the disease and host-pathogen interactions to identify new drug targets for intervention strategies. Using in-silico comparative genome analysis, we identified one of the M. tb genes, Rv1509, as a signature protein exclusively present in M. tb. To explore the role of Rv1509, a likely methyl transferase, we constructed a knock-in Mycobacterium smegmatis (M. smegmatis) constitutively expressing Rv1509 (Ms_Rv1509). The Ms_Rv1509 led to differential expression of many transcriptional regulator genes as assessed by RNA-seq analysis. Further, in-vitro and in-vivo studies demonstrated an enhanced survival of Ms_Rv1509 inside the host macrophages. Ms_Rv1509 also promoted phagolysosomal escape inside macrophages to boost bacterial replication and dissemination. In-vivo infection studies revealed that Ms_Rv1509 survives better than BCG and causes pathological manifestations in the pancreas after intraperitoneal infection. Long-time survival of Ms_Rv1509 resulted in lymphocyte migration, increased T regulatory cells, giant cell formation, and likely granuloma formation in the pancreas, pointing toward the role of Rv1509 in M. tb pathogenesis.

Pangenome Reconstruction of Mycobacterium tuberculosis as a Guide to Reveal Genomic Features Associated with Strain Clinical Phenotype

Tuberculosis (TB) is one of the leading causes of human deaths worldwide caused by infectious diseases. TB infection by Mycobacterium tuberculosis can occur in the lungs, causing pulmonary tuberculosis (PTB), or in any other organ of the body, resulting in extrapulmonary tuberculosis (EPTB). There is no consensus on the genetic determinants of this pathogen that may contribute to EPTB. In this study, we constructed the M. tuberculosis pangenome and used it as a tool to seek genomic signatures associated with the clinical presentation of TB based on its accessory genome differences. The analysis carried out in the present study includes the raw reads of 490 M. tuberculosis genomes (PTB n = 245, EPTB n = 245) retrieved from public databases that were assembled, as well as ten genomes from Mexican strains (PTB n = 5, EPTB n = 5) that were sequenced and assembled. All genomes were annotated and then used to construct the pangenome with Roary and Panaroo. The pangenome obtained using Roary consisted of 2231 core genes and 3729 accessory genes. On the other hand, the pangenome resulting from Panaroo consisted of 2130 core genes and 5598 accessory genes. Associations between the distribution of accessory genes and the PTB/EPTB phenotypes were examined using the Scoary and Pyseer tools. Both tools found a significant association between the hspR, plcD, Rv2550c, pe_pgrs5, pe_pgrs25, and pe_pgrs57 genes and the PTB genotype. In contrast, the deletion of the aceA, esxR, plcA, and ppe50 genes was significantly associated with the EPTB phenotype. Rv1759c and Rv3740 were found to be associated with the PTB phenotype according to Scoary; however, these associations were not observed when using Pyseer. The robustness of the constructed pangenome and the gene-phenotype associations is supported by several factors, including the analysis of a large number of genomes, the inclusion of the same number of PTB/EPTB genomes, and the reproducibility of results thanks to the different bioinformatic tools used. Such characteristics surpass most of previous M. tuberculosis pangenomes. Thus, it can be inferred that the deletion of these genes can lead to changes in the processes involved in stress response and fatty acid metabolism, conferring phenotypic advantages associated with pulmonary or extrapulmonary presentation of TB. This study represents the first attempt to use the pangenome to seek gene-phenotype associations in M. tuberculosis.

Novel role of PE_PGRS47 in the alteration of mycobacterial cell wall integrity and drug resistance

The proline-glutamic acid and proline-proline-glutamic acid (PE/PPE) family of proteins is widespread in pathogenic mycobacteria and plays different roles in mycobacterial physiology. While several PE/PPE family proteins have been studied, the exact function of most PE/PPE proteins in the physiology of Mycobacterium tuberculosis (Mtb) remains unknown. PE_PGRS47 belongs to the PE/PPE family of proteins reported to help Mtb evade protective host immune responses. In this study, we demonstrate a novel role of PE_PGRS47. Heterologous expression of the pe_pgrs47 gene in a non-pathogenic Mycobacterium smegmatis, intrinsically deficient of PE_PGRS protein, exhibits modulated colony morphology and cell wall lipid profile leading to a marked susceptibility to multiple antibiotics and environmental stressors. Using ethidium bromide/Nile red uptake assays, Mycobacterium smegmatis expressing PE_PGRS47 showed higher cell wall permeability than the control strain. Overall, these data suggested that PE_PGRS47 is cell surface exposed and influences cell wall integrity and the formation of mycobacterial colonies, ultimately potentiating the efficacy of lethal stresses against mycobacteria.

Human alveolar macrophage metabolism is compromised during Mycobacterium tuberculosis infection

Pulmonary macrophages have two distinct ontogenies: long-lived embryonically-seeded alveolar macrophages (AM) and bone marrow-derived macrophages (BMDM). Here, we show that after infection with a virulent strain of Mycobacterium tuberculosis (H37Rv), primary murine AM exhibit a unique transcriptomic signature characterized by metabolic reprogramming distinct from conventional BMDM. In contrast to BMDM, AM failed to shift from oxidative phosphorylation (OXPHOS) to glycolysis and consequently were unable to control infection with an avirulent strain (H37Ra). Importantly, healthy human AM infected with H37Ra equally demonstrated diminished energetics, recapitulating our observation in the murine model system. However, the results from seahorse showed that the shift towards glycolysis in both AM and BMDM was inhibited by H37Rv. We further demonstrated that pharmacological (e.g. metformin or the iron chelator desferrioxamine) reprogramming of AM towards glycolysis reduced necrosis and enhanced AM capacity to control H37Rv growth. Together, our results indicate that the unique bioenergetics of AM renders these cells a perfect target for Mtb survival and that metabolic reprogramming may be a viable host targeted therapy against TB.

The exploitation of host autophagy and ubiquitin machinery by Mycobacterium tuberculosis in shaping immune responses and host defense during infection

Intracellular pathogens have evolved various efficient molecular armaments to subvert innate defenses. Cellular ubiquitination, a normal physiological process to maintain homeostasis, is emerging one such exploited mechanism. Ubiquitin (Ub), a small protein modifier, is conjugated to diverse protein substrates to regulate many functions. Structurally diverse linkages of poly-Ub to target proteins allow enormous functional diversity with specificity being governed by evolutionarily conserved enzymes (E3-Ub ligases). The Ub-binding domain (UBD) and LC3-interacting region (LIR) are critical features of macroautophagy/autophagy receptors that recognize Ub-conjugated on protein substrates. Emerging evidence suggests that E3-Ub ligases unexpectedly protect against intracellular pathogens by tagging poly-Ub on their surfaces and targeting them to phagophores. Two E3-Ub ligases, PRKN and SMURF1, provide immunity against Mycobacterium tuberculosis (M. tb). Both enzymes conjugate K63 and K48-linked poly-Ub to M. tb for successful delivery to phagophores. Intriguingly, M. tb exploits virulence factors to effectively dampen host-directed autophagy utilizing diverse mechanisms. Autophagy receptors contain LIR-motifs that interact with conserved Atg8-family proteins to modulate phagophore biogenesis and fusion to the lysosome. Intracellular pathogens have evolved a vast repertoire of virulence effectors to subdue host-immunity via hijacking the host ubiquitination process. This review highlights the xenophagy-mediated clearance of M. tb involving host E3-Ub ligases and counter-strategy of autophagy inhibition by M. tb using virulence factors. The role of Ub-binding receptors and their mode of autophagy regulation is also explained. We also discuss the co-opting and utilization of the host Ub system by M. tb for its survival and virulence.Abbreviations: APC: anaphase promoting complex/cyclosome; ATG5: autophagy related 5; BCG: bacille Calmette-Guerin; C2: Ca²⁺-binding motif; CALCOCO2: calcium binding and coiled-coil domain 2; CUE: coupling of ubiquitin conjugation to ER degradation domains; DUB: deubiquitinating enzyme; GABARAP: GABA type A receptor-associated protein; HECT: homologous to the E6-AP carboxyl terminus; IBR: in-between-ring fingers; IFN: interferon; IL1B: interleukin 1 beta; KEAP1: kelch like ECH associated protein 1; LAMP1: lysosomal associated membrane protein 1; LGALS: galectin; LIR: LC3-interacting region; MAPK11/p38: mitogen-activated protein kinase 11; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MAP3K7/TAK1: mitogen-activated protein kinase kinase kinase 7; MAPK8/JNK: mitogen-activated protein kinase 8; MHC-II: major histocompatibility complex-II; MTOR: mechanistic target of rapamycin kinase; NBR1: NBR1 autophagy cargo receptor; NFKB1/p50: nuclear factor kappa B subunit 1; OPTN: optineurin; PB1: phox and bem 1; PE/PPE: proline-glutamic acid/proline-proline-glutamic acid; PknG: serine/threonine-protein kinase PknG; PRKN: parkin RBR E3 ubiquitin protein ligase; RBR: RING-in between RING; RING: really interesting new gene; RNF166: RING finger protein 166; ROS: reactive oxygen species; SMURF1: SMAD specific E3 ubiquitin protein ligase 1; SQSTM1: sequestosome 1; STING1: stimulator of interferon response cGAMP interactor 1; TAX1BP1: Tax1 binding protein 1; TBK1: TANK binding kinase 1; TNF: tumor necrosis factor; TRAF6: TNF receptor associated factor 6; Ub: ubiquitin; UBA: ubiquitin-associated; UBAN: ubiquitin-binding domain in ABIN proteins and NEMO; UBD: ubiquitin-binding domain; UBL: ubiquitin-like; ULK1: unc-51 like autophagy activating kinase 1.

Comprehensive sequence and structure analysis of algal lipid catabolic enzyme Triacylglycerol lipase: an in silico study to vitalize the development of optimum engineered strains with high lipid productivity

Microalgae as an alternative renewable resource for biofuel production have captured much significance. Nonetheless, its economic viability is a field of major concern for researchers. Unraveling the lipid catabolic pathway and gaining insights into the sequence-structural features of its primary functioning enzyme, Triacylglycerol lipase, will impart valuable information to target microalgae for augmented lipid content. In the present study, a genome-wide comparative study on putative Triacylglycerol lipase (TAGL) enzyme from algal species belonging to varied phylogenetic lineages was performed. The comprehensive sequence analysis revealed that TAGL comprises of three distinct conserved domains, such as, Patatin, Class III Lipase, and Abhydro_lipase, and also confirmed the ubiquitous presence of GXSXG motif in the sequences analyzed. In the absence of a crystal structure of algal TAGL till date, we developed the first 3D model of patatin domain of TAGL from an oleaginous microalga, Phaedactylum tricornutum, employing homology modeling, docking and molecular dynamic simulations methods. The domain-substrate complex having the low-ranking docking score revealed the binding of palmitic acid to the TAGL patatin domain surface with strong hydrogen bond interactions. The simulation results implied that the substrate-complexed patatin domain and the free enzyme adopted a more stable conformation after 40 ns. This is the first ever attempt to provide in-silico insights into the structural and dynamical insights on catalytic mechanism of the TAGL patatin domain. Subsequently, these findings aided our understanding on their structural stability, folding mechanism and protein-substrate interactions, which could be further utilized to design site-specific mutagenic experiments for engineering microalgal strains.Communicated by Ramaswamy H. Sarma.

The Mycobacterium tuberculosis PE_PGRS Protein Family Acts as an Immunological Decoy to Subvert Host Immune Response

Mycobacterium tuberculosis (M.tb) is a successful pathogen that can reside within the alveolar macrophages of the host and can survive in a latent stage. The pathogen has evolved and developed multiple strategies to resist the host immune responses. M.tb escapes from host macrophage through evasion or subversion of immune effector functions. M.tb genome codes for PE/PPE/PE_PGRS proteins, which are intrinsically disordered, redundant and antigenic in nature. These proteins perform multiple functions that intensify the virulence competence of M.tb majorly by modulating immune responses, thereby affecting immune mediated clearance of the pathogen. The highly repetitive, redundant and antigenic nature of PE/PPE/PE_PGRS proteins provide a critical edge over other M.tb proteins in terms of imparting a higher level of virulence and also as a decoy molecule that masks the effect of effector molecules, thereby modulating immuno-surveillance. An understanding of how these proteins subvert the host immunological machinery may add to the current knowledge about M.tb virulence and pathogenesis. This can help in redirecting our strategies for tackling M.tb infections.

PGRS Domain of Rv0297 of Mycobacterium tuberculosis Functions in A Calcium Dependent Manner

Mycobacterium tuberculosis (M.tb), the pathogen causing tuberculosis, is a major threat to human health worldwide. Nearly 10% of M.tb genome encodes for a unique family of PE/PPE/PGRS proteins present exclusively in the genus Mycobacterium. The functions of most of these proteins are yet unexplored. The PGRS domains of these proteins have been hypothesized to consist of Ca²⁺ binding motifs that help these intrinsically disordered proteins to modulate the host cellular responses. Ca²⁺ is an important secondary messenger that is involved in the pathogenesis of tuberculosis in diverse ways. This study presents the calcium-dependent function of the PGRS domain of Rv0297 (PE_PGRS5) in M.tb virulence and pathogenesis. Tandem repeat search revealed the presence of repetitive Ca²⁺ binding motifs in the PGRS domain of the Rv0297 protein (Rv0297PGRS). Molecular Dynamics simulations and fluorescence spectroscopy revealed Ca²⁺ dependent stabilization of the Rv0297PGRS protein. Calcium stabilized Rv0297PGRS enhances the interaction of Rv0297PGRS with surface localized Toll like receptor 4 (TLR4) of macrophages. The Ca²⁺ stabilized binding of Rv0297PGRS with the surface receptor of macrophages enhances its downstream consequences in terms of Nitric Oxide (NO) production and cytokine release. Thus, this study points to hitherto unidentified roles of calcium-modulated PE_PGRS proteins in the virulence of M.tb. Understanding the pathogenic potential of Ca²⁺ dependent PE_PGRS proteins can aid in targeting these proteins for therapeutic interventions.

Teleological cooption of Mycobacterium tuberculosis PE/PPE proteins as porins: Role in molecular immigration and emigration

Permeation through bacterial cells for exchange or uptake of biomolecules and ions invariably depend upon the existence of pore-forming proteins (porins) in their outer membrane. Mycobacterium tuberculosis (M. tb) harbours one of the most rigid cell envelopes across bacterial genera and is devoid of the classical porins for solute transport across the cell membrane. Though canonical porins are incompatible with the evolution of permeability barrier, porin like activity has been reported from membrane preparations of pathogenic mycobacteria. This suggests a sophisticated transport mechanism that has been elusive until now, along with the protein family responsible for it. Recent evidence suggests that these slow-growing mycobacteria have co-opted some of PE/PPE family proteins as molecular transport channels, in place of porins, to facilitate uptake of nutrients required to thrive in the restrictive host environment. These reports advocate that PE/PPE proteins, due to their structural ability, have a potential role in importing small molecules to the cell's interior. This mechanism unveils how a successful pathogen overcomes its restrictive membrane's transport limitations for selective uptake of nutrients. If extrapolated to have a role in drug transport, these channels could help understand the emergence of drug resistance. Further, as these proteins are associated with the export of virulence factors, they can be exploited as novel drug targets. There remains, however, an interesting question that as the PE/PPE proteins can allow the 'import' of molecules from outside the cell, is the reverse transport also possible across the M. tb membrane. In this review, we have discussed recent evidence supporting PE/PPE's role as a specific transport channel for selective uptake of small molecule nutrients and, as possible molecular export machinery of M. tb. This newly discovered role as transmembrane channels demands further research on this enigmatic family of proteins to comprehend the pathomechanism of this very smart pathogen.

Mycobacterium tuberculosis RipA Dampens TLR4-Mediated Host Protective Response Using a Multi-Pronged Approach Involving Autophagy, Apoptosis, Metabolic Repurposing, and Immune Modulation

Reductive evolution has endowed Mycobacterium tuberculosis (M. tb) with moonlighting in protein functions. We demonstrate that RipA (Rv1477), a peptidoglycan hydrolase, activates the NFκB signaling pathway and elicits the production of pro-inflammatory cytokines, TNF-α, IL-6, and IL-12, through the activation of an innate immune-receptor, toll-like receptor (TLR)4. RipA also induces an enhanced expression of macrophage activation markers MHC-II, CD80, and CD86, suggestive of M1 polarization. RipA harbors LC3 (Microtubule-associated protein 1A/1B-light chain 3) motifs known to be involved in autophagy regulation and indeed alters the levels of autophagy markers LC3BII and P62/SQSTM1 (Sequestosome-1), along with an increase in the ratio of P62/Beclin1, a hallmark of autophagy inhibition. The use of pharmacological agents, rapamycin and bafilomycin A1, reveals that RipA activates PI3K-AKT-mTORC1 signaling cascade that ultimately culminates in the inhibition of autophagy initiating kinase ULK1 (Unc-51 like autophagy activating kinase). This inhibition of autophagy translates into efficient intracellular survival, within macrophages, of recombinant Mycobacterium smegmatis expressing M. tb RipA. RipA, which also localizes into mitochondria, inhibits the production of oxidative phosphorylation enzymes to promote a Warburg-like phenotype in macrophages that favors bacterial replication. Furthermore, RipA also inhibited caspase-dependent programed cell death in macrophages, thus hindering an efficient innate antibacterial response. Collectively, our results highlight the role of an endopeptidase to create a permissive replication niche in host cells by inducing the repression of autophagy and apoptosis, along with metabolic reprogramming, and pointing to the role of RipA in disease pathogenesis.

Distribution of sigma factors delineates segregation of virulent and avirulent Mycobacterium

The genus Mycobacterium includes a wide range of species of both slow and rapid growth under major pathogens, opportunists, and saprophytes. The number and combination of sigma factors are extremely diversified among various species of Mycobacterium. The comparative genome analysis illustrates that SigC, SigD, SigG, SigH, SigK and SigI are dominant among the pathogens. Evolutionary analysis using Bayesian inference on 16S rRNA and MLST-based phylogeny using 14 housekeeping genes distinctly differentiate the slow-growing Mycobacterium from fast growers and segregate pathogens from opportunists and saprophytes. Based on the similarity coefficient upon the allotment of sigma factors in mycobacterial species through UPGMA dendrogram analysis, it is apparent that the pathogens are grouped separately following the similar trend observed from the evolutionary approach. Predominance of a set of sigma factors particularly the pathogenic Mycobacterium co-exists with the distribution of six well-known virulence factors of Mycobacterium (PhoP, PcaA, FbpA, Mce1B, KatG and PE_PGRS30). The pathogenicity responsible sigma factors elicit close resemblance with few notable characters of the known virulence factors. Thus the analysis renders that the distribution of sigma factors of different species of Mycobacterium can be a potential tool to predict their pathogenicity index.

Comprehensive analysis of protein acetyltransferases of human pathogen Mycobacterium tuberculosis

Tuberculosis (TB), a leading infectious disease caused by Mycobacterium tuberculosis strain, takes four human lives every minute globally. Paucity of knowledge on M. tuberculosis virulence and antibiotic resistance is the major challenge for tuberculosis control. We have identified 47 acetyltransferases in the M. tuberculosis, which use diverse substrates including antibiotic, amino acids, and other chemical molecules. Through comparative analysis of the protein file of the virulent M. tuberculosis H37Rv strain and the avirulent M. tuberculosis H37Ra strain, we identified one acetyltransferase that shows significant variations with N-terminal deletion, possibly influencing its physicochemical properties. We also found that one acetyltransferase has three types of post-translation modifications (lysine acetylation, succinylation, and glutarylation). The genome context analysis showed that many acetyltransferases with their neighboring genes belong to one operon. By data mining from published transcriptional profiles of M. tuberculosis exposed to diverse treatments, we revealed that several acetyltransferases may be functional during M. tuberculosis infection. Insights obtained from the present study can potentially provide clues for developing novel TB therapeutic interventions.

PGRS Domain of Rv0297 of Mycobacterium tuberculosis Is Involved in Modulation of Macrophage Functions to Favor Bacterial Persistence

Mycobacterium tuberculosis (M. tb) Rv0297-encoded PE_PGRS5 has been known to be expressed at the later stages of infection and in acidified phagosomes during transcriptome and proteomic studies. The possible role of Rv0297 in the modulation of phagosomal maturation and in providing protection against a microbicidal environment has been hypothesized. We show that Rv0297PGRS is involved in modulating the calcium homeostasis of macrophages followed by impedance of the phagolysosomal acidification process. This is evident from the downregulation of the late endosomal markers (Rab7 and cathepsin D) in the macrophages infected with recombinant Mycobacterium smegmatis (rM.smeg)—M.smeg_Rv0297 and M.smeg_Rv0297PGRS—or treated with recombinant Rv0297PGRS protein. Macrophages infected with rM.smeg expressing Rv0297 produce nitric oxide and undergo apoptosis, which may aid in the dissemination of pathogen in the later stages of infection. Rv0297 was also found to be involved in rescuing the bacterium from oxidative and hypoxic stress employed by macrophages and augmented the survivability of the recombinant bacterium. These results attribute to the functional significance of this protein in M.tb virulence mechanism. The fact that this protein gets expressed at the later stages of lung granulomas during M.tb infection suggests that the bacterium possibly employs Rv0297 as its dissemination and survival strategy.

Mucosal-Associated Invariant T Cells Develop an Innate-Like Transcriptomic Program in Anti-mycobacterial Responses

Conventional T cells exhibit a delayed response to the initial priming of peptide antigens presented by major histocompatibility complex (MHC) proteins. Unlike conventional T cells, mucosal-associated invariant T (MAIT) cells quickly respond to non-peptidic metabolite antigens presented by MHC-related protein 1 (MR1). To elucidate the MR1-dependent activation program of MAIT cells in response to mycobacterial infections, we determined the surface markers, transcriptomic profiles, and effector responses of activated human MAIT cells. Results revealed that mycobacterial-incubated antigen-presenting cells stimulated abundant human CD8⁺ MAIT cells to upregulate the co-expression of CD69 and CD26, as a combinatorial activation marker. Further transcriptomic analyses demonstrated that CD69⁺CD26⁺⁺ CD8⁺MAIT cells highly expressed numerous genes for mediating anti-mycobacterial immune responses, including pro-inflammatory cytokines, cytolytic molecules, NK cell receptors, and transcription factors, in contrast to inactivated counterparts CD69^+/−CD26^+/− CD8⁺MAIT cells. Gene co-expression, enrichment, and pathway analyses yielded high statistical significance to strongly support that activated CD8⁺ MAIT cells shared gene expression and numerous pathways with NK and CD8⁺ T cells in activation, cytokine production, cytokine signaling, and effector functions. Flow cytometry detected that activated CD8⁺MAIT cells produced TNFα, IFNγ, and granulysin to inhibit mycobacterial growth and fight mycobacterial infection. Together, results strongly support that the combinatorial activation marker CD69⁺CD26⁺⁺ labels the activated CD8⁺MAIT cells that develop an innate-like activation program in anti-mycobacterial immune responses. We speculate that the rapid production of anti-mycobacterial effector molecules facilitates MAIT cells to fight early mycobacterial infection in humans.

Disorder-to-order transition in PE-PPE proteins of Mycobacterium tuberculosis augments the pro-pathogen immune response

A growing body of evidence supports the hypothesis that intrinsically disordered proteins often mediate host–pathogen interactions and modulate host functions for pathogen survival and virulence. Mycobacterium tuberculosis (M.tb) has evolved largely through reductive evolution, with a few exceptions such as the glycine–alanine‐rich PE–PPE/PGRS protein family, which has been expanding in pathogenic mycobacteria. Here, our analyses of the M.tb proteome and secretome revealed that the PE–PGRS subfamily is enriched for disordered regions and disordered binding sites, pointing to their importance in host–pathogen interactions. As a case study, the secondary structure of PE35–PPE68 and PE32–PPE65 of the pathogenesis‐related RD1 and RD8 regions was analyzed through Fourier‐transform infrared spectroscopy. These disordered proteins displayed a considerable structural shift from disordered to ordered while engaged in the formation of complexes. While these proteins are immunogenic individually and enhance the pro‐pathogen response, their corresponding complexes enhanced the responses manifold as displayed here by PE35 and PPE68. It is likely that M.tb exploits such disorder–order structural dynamics as a strategy to mount a pro‐pathogen response and subvert host defense for productive infection. This functional gain also serves as a means to compensate genomic content loss due to reductive evolution.

Expression and regulatory networks of Mycobacterium tuberculosis PE/PPE family antigens

PE/PPE family antigens are distributed mainly in pathogenic mycobacteria and serve as potential antituberculosis (TB) vaccine components. Some PE/PPE family antigens can regulate the host innate immune response, interfere with macrophage activation and phagolysosome fusion, and serve as major sources of antigenic variation. PE/PPE antigens have been associated with mycobacteria pathogenesis; pe/ppe genes are mainly found in pathogenic mycobacteria and are differentially expressed between Mtb and Mycobacterium bovis. PE/PPE proteins were essential for the growth of Mtb, and PE/PPE proteins were differentially expressed under a variety of conditions. Multiple mycobacterial-virulence-related transcription factors, sigma factors, the global transcriptional regulation factor Lsr2, MprAB, and PhoPR two-component regulatory systems, and cyclic adenine monophosphate-dependent regulators, regulate the expression of PE/PPE family antigens. Multiple-scale integrative analysis revealed the expression and regulatory networks of PE/PPE family antigens underlying the virulence and pathogenesis of Mtb, providing important clues for the discovery of new anti-TB measures.

Linking common non-coding RNAs of human lung cancer and M. tuberculosis

Lung cancer and pulmonary tuberculosis caused by Mycobacterium are two major causes of deaths worldwide. Tuberculosis linked lung cancer is known. However, the precise molecular mechanism of Mycobacterium associated increased risk of lung cancer is not understood. We report 45 common human miRNAs deregulated in both pulmonary tuberculosis and lung cancer. We show that sRNA_1096 and sRNA_1414 from M. tuberculosis have sequence homology with human mir-21. Hence, the potential role of these three small non-coding RNAs in rifampicin resistance in pulmonary tuberculosis is implied. Further, the linking of sRNA_1096 and sRNA_1414 from M. tuberculosis with the host lung tumorigenesis is inferred. Nonetheless, further analysis and validation is required to associate these three non-coding RNAs with Mycobacterium associated increased risk of lung cancer.

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

The 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 Ca²⁺ 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.

Apoptosis is generally thought to be a defense mechanism in protecting the host against Mycobacterium tuberculosis in early stages of infection. However, apoptosis during later stages in lung granulomas may favor the bacterium in disseminating the disease. ER stress has been found to induce apoptosis in TB granulomas, in zones where apoptotic macrophages accumulate in mice and humans. In this study, we report ER stress-mediated apoptosis of host cells by the Rv0297-encoded PE_PGRS5 protein of M. tuberculosis exceptionally present in the pathogenic Mycobacterium genus. The PGRS domain of Rv0297 aids the protein in localizing to the ER and induces the unfolded protein response followed by apoptosis of macrophages. The effect of the Rv0297PGRS domain was found to be TLR4 dependent. This study presents novel insights on the strategies employed by M. tuberculosis to disseminate the disease.

Cadmium-induced conformational changes in type 2 metallothionein of medicinal plant Coptis japonica: insights from molecular dynamics studies of apo, partially and fully metalated forms

Plants play an important role in the removal of excess heavy metals from soil and water. Medicinal plants can also have non-traditional use in phytoremediation technologies. Among the heavy metals, Cadmium (Cd) is the most abundant and readily taken up by the crop plants. Plant metallothioneins (MTs) are small proteins having cysteine-rich residues and appear to play key roles in metal homoeostasis. Plant metallothionein 2 (MT 2) from Coptis japonica (Gold-thread; CjMT 2) is a typical member of this subfamily and features two cysteine-rich regions containing eight and six cysteine residues, respectively, separated by 42 amino acids long linker region. In-silico analysis of MT 2 protein sequences of C. japonica was performed. In this study, ab initio methods were utilised for the prediction of three-dimensional structure of CjMT 2. After structure validation, heavy metal-binding sites were predicted for the selected modelled structures of CjMT 2. To obtain Cd_i-CjMT 2 (i = 1-7), metalated complex individual docking experiments were performed. The stability of the metalated docked structures was assessed by molecular dynamics (MD) simulation studies. Our study showed that CjMT 2 binds up to 4 Cd²⁺ ions in two distinct domains: a N-terminal β-domain that binds to 2 Cd²⁺ ions and a C-terminal α-domain that binds with 2 Cd²⁺ ions. Our analysis revealed that Cys residues of alpha and beta domain and some residues of spacer region of CjMT 2 protein might be important for the cadmium interaction. MD simulation studies provided insight into metal-induced conformational changes and mechanism of metalation of CjMT 2, an intrinsically disordered protein. This study provides useful insights into mechanism of cadmium-type 2 metallothionein interaction.

Comparative Genomic Analysis of Two Clonally Related Multidrug Resistant Mycobacterium tuberculosis by Single Molecule Real Time Sequencing

Background: Multidrug-resistant tuberculosis (MDR-TB) is posing a major threat to global TB control. In this study, we focused on two consecutive MDR-TB isolated from the same patient before and after the initiation of anti-TB treatment. To better understand the genomic characteristics of MDR-TB, Single Molecule Real-Time (SMRT) Sequencing and comparative genomic analyses was performed to identify mutations that contributed to the stepwise development of drug resistance and growth fitness in MDR-TB under in vivo challenge of anti-TB drugs. Result: Both pre-treatment and post-treatment strain demonstrated concordant phenotypic and genotypic susceptibility profiles toward rifampicin, pyrazinamide, streptomycin, fluoroquinolones, aminoglycosides, cycloserine, ethionamide, and para-aminosalicylic acid. However, although both strains carried identical missense mutations at rpoB S531L, inhA C-15T, and embB M306V, MYCOTB Sensititre assay showed that the post-treatment strain had 16-, 8-, and 4-fold elevation in the minimum inhibitory concentrations (MICs) toward rifabutin, isoniazid, and ethambutol respectively. The results have indicated the presence of additional resistant-related mutations governing the stepwise development of MDR-TB. Further comparative genomic analyses have identified three additional polymorphisms between the clinical isolates. These include a single nucleotide deletion at nucleotide position 360 of rv0888 in pre-treatment strain, and a missense mutation at rv3303c (lpdA) V44I and a 6-bp inframe deletion at codon 67-68 in rv2071c (cobM) in the post-treatment strain. Multiple sequence alignment showed that these mutations were occurring at highly conserved regions among pathogenic mycobacteria. Using structural-based and sequence-based algorithms, we further predicted that the mutations potentially have deleterious effect on protein function. Conclusion: This is the first study that compared the full genomes of two clonally-related MDR-TB clinical isolates during the course of anti-TB treatment. Our work has demonstrated the robustness of SMRT Sequencing in identifying mutations among MDR-TB clinical isolates. Comparative genome analysis also suggested novel mutations at rv0888, lpdA, and cobM that might explain the difference in antibiotic resistance and growth pattern between the two MDR-TB strains.

Structural basis of the PE-PPE protein interaction in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, has developed multiple strategies to adapt to the human host. The five type VII secretion systems, ESX-1-5, direct the export of many virulence-promoting protein effectors across the complex mycobacterial cell wall. One class of ESX substrates is the PE-PPE family of proteins, which is unique to mycobacteria and essential for infection, antigenic variation, and host-pathogen interactions. The genome of Mtb encodes 168 PE-PPE proteins. Many of them are thought to be secreted through ESX-5 secretion system and to function in pairs. However, understanding of the specific pairing of PE-PPE proteins and their structure-function relationship is limited by the challenging purification of many PE-PPE proteins, and our knowledge of the PE-PPE interactions therefore has been restricted to the PE25-PPE41 pair and its complex with the ESX-5 secretion system chaperone EspG5. Here, we report the crystal structure of a new PE-PPE pair, PE8-PPE15, in complex with EspG5. Our structure revealed that the EspG5-binding sites on PPE15 are relatively conserved among Mtb PPE proteins, suggesting that EspG5-PPE15 represents a more typical model for EspG5-PPE interactions than EspG5-PPE41. A structural comparison with the PE25-PPE41 complex disclosed conformational changes in the four-helix bundle structure and a unique binding mode in the PE8-PPE15 pair. Moreover, homology-modeling and mutagenesis studies further delineated the molecular determinants of the specific PE-PPE interactions. These findings help develop an atomic algorithm of ESX-5 substrate recognition and PE-PPE pairing.

SMRT genome assembly corrects reference errors, resolving the genetic basis of virulence in Mycobacterium tuberculosis

BACKGROUND

The genetic basis of virulence in Mycobacterium tuberculosis has been investigated through genome comparisons of virulent (H37Rv) and attenuated (H37Ra) sister strains. Such analysis, however, relies heavily on the accuracy of the sequences. While the H37Rv reference genome has had several corrections to date, that of H37Ra is unmodified since its original publication.

RESULTS

Here, we report the assembly and finishing of the H37Ra genome from single-molecule, real-time (SMRT) sequencing. Our assembly reveals that the number of H37Ra-specific variants is less than half of what the Sanger-based H37Ra reference sequence indicates, undermining and, in some cases, invalidating the conclusions of several studies. PE_PPE family genes, which are intractable to commonly-used sequencing platforms because of their repetitive and GC-rich nature, are overrepresented in the set of genes in which all reported H37Ra-specific variants are contradicted. Further, one of the sequencing errors in H37Ra masks a true variant in common with the clinical strain CDC1551 which, when considered in the context of previous work, corresponds to a sequencing error in the H37Rv reference genome.

CONCLUSIONS

Our results constrain the set of genomic differences possibly affecting virulence by more than half, which focuses laboratory investigation on pertinent targets and demonstrates the power of SMRT sequencing for producing high-quality reference genomes.

Identifying and characterising PPE7 (Rv0354c) high activity binding peptides and their role in inhibiting cell invasion

This study was aimed at characterising the PPE7 protein from the PE/PPE protein family. The presence and transcription of the rv0354c gene in the Mycobacterium tuberculosis complex was determined and the subcellular localisation of the PPE7 protein on mycobacterial membrane was confirmed by immunoelectron microscope. Two peptides were identified as having high binding activity (HABPs) and were tested in vitro regarding the invasion of Mycobacterium tuberculosis H37Rv. HABP 39224 inhibited invasion in A549 epithelial cells and U937 macrophages by more than 50%, whilst HABP 39225 inhibited invasion by 40% in U937 cells. HABP 39224, located in the protein's C-terminal region, has a completely conserved amino acid sequence in M. tuberculosis complex species and could be selected as a base peptide when designing a subunit-based, anti-tuberculosis vaccine.

India's Computational Biology Growth and Challenges

India's computational science is growing swiftly due to the outburst of internet and information technology services. The bioinformatics sector of India has been transforming rapidly by creating a competitive position in global bioinformatics market. Bioinformatics is widely used across India to address a wide range of biological issues. Recently, computational researchers and biologists are collaborating in projects such as database development, sequence analysis, genomic prospects and algorithm generations. In this paper, we have presented the Indian computational biology scenario highlighting bioinformatics-related educational activities, manpower development, internet boom, service industry, research activities, conferences and trainings undertaken by the corporate and government sectors. Nonetheless, this new field of science faces lots of challenges.

Mycobacterium tuberculosis Co-operonic PE32/PPE65 Proteins Alter Host Immune Responses by Hampering Th1 Response

PE/PPE genes, present in cluster with ESAT-6 like genes, are suspected to have a role in antigenic variation and virulence of Mycobacterium tuberculosis. Their roles in immune evasion and immune modulation of host are also well documented. We present evidence that PE32/PPE65 present within the RD8 region are co-operonic, co-transcribed, and co-translated, and play role in modulating host immune responses. Experiments with macrophage cell lines revealed that this protein complex suppresses pro-inflammatory cytokines such as TNF-α and IL-6 whereas also inducing high expression of anti-inflammatory IL-10. Immunization of mice with these recombinant proteins dampens an effective Th1 response as evident from reduced frequency of IFN-γ and IL-2 producing CD4⁺ and CD8⁺ T cells. IgG sub-typing from serum of immunized mice revealed high levels of IgG1 when compared with IgG2a and IgG2b. Further IgG1/IgG2a ratio clearly demonstrated that the protein complex manipulates the host immune response favorable to the pathogen. Our results demonstrate that the co-transcribed and co-translated PE32 and PPE65 antigens are involved specifically in modulating anti-mycobacterial host immune response by hampering Th1 response.

Analyses of methyltransferases across the pathogenicity spectrum of different mycobacterial species point to an extremophile connection

Tuberculosis is a devastating disease, taking one human life every 20 seconds globally. We hypothesize that professional pathogens such as M.tb have acquired specific features that might assist in causing infection, persistence and transmissible pathology in their host. We have identified 121 methyltransferases (MTases) in the M.tb proteome, which use a variety of substrates - DNA, RNA, protein, intermediates of mycolic acid biosynthesis and other fatty acids - that are involved in cellular maintenance within the host. A comparative analysis of the proteome of the virulent strain H37Rv and the avirulent strain H37Ra identified 3 MTases, which displayed significant variations in terms of N-terminal extension/deletion and point mutations, possibly impacting various physicochemical properties. The cross-proteomic comparison of MTases of M.tb H37Rv with 15 different Mycobacterium species revealed the acquisition of novel MTases in a MTB complex as a function of evolution. Phylogenetic analysis revealed that these newly acquired MTases showed common roots with certain extremophiles such as halophilic and acidophilic organisms. Our results establish an evolutionary relationship of M.tb with halotolerant organisms and also the role of MTases of M.tb in withstanding the host osmotic stress, thereby pointing to their likely role in pathogenesis, virulence and niche adaptation.

Interaction of Mycobacterium tuberculosis Virulence Factor RipA with Chaperone MoxR1 Is Required for Transport through the TAT Secretion System

Mycobacterium tuberculosis is a leading cause of death worldwide. The M. tuberculosis TAT (twin-arginine translocation) protein secretion system is present at the cytoplasmic membrane of mycobacteria and is known to transport folded proteins. The TAT secretion system is reported to be essential for many important bacterial processes that include cell wall biosynthesis. The M. tuberculosis secretion and invasion protein RipA has endopeptidase activity and interacts with one of the resuscitation antigens (RpfB) that are expressed during pathogen reactivation. MoxR1, a member of the ATPase family that is associated with various cellular activities, was predicted to interact with RipA based on in silico analyses. A bimolecular fluorescence complementation (BiFC) assay confirmed the interaction of these two proteins in HEK293T cells. The overexpression of RipA in Mycobacterium smegmatis and copurification with MoxR1 further validated their interaction in vivo. Recombinant MoxR1 protein, expressed in Escherichia coli, displays ATP-enhanced chaperone activity. Secretion of recombinant RipA (rRipA) protein into the E. coli culture filtrate was not observed in the absence of RipA-MoxR interaction. Inhibition of this export system in M. tuberculosis, including the key players, will prevent localization of peptidoglycan hydrolase and result in sensitivity to existing β-lactam antibiotics, opening up new candidates for drug repurposing.

The virulence mechanism of mycobacteria is very complex. Broadly, the virulence factors can be classified as secretion factors, cell surface components, enzymes involved in cellular metabolism, and transcriptional regulators. The mycobacteria have evolved several mechanisms to secrete its proteins. Here, we have identified one of the virulence proteins of Mycobacterium tuberculosis, RipA, possessing peptidoglycan hydrolase activities secreted by the TAT secretion pathway. We also identified MoxR1 as a protein-protein interaction partner of RipA and demonstrated chaperone activity of this protein. We show that MoxR1-mediated folding is critical for the secretion of RipA within the TAT system. Inhibition of this export system in M. tuberculosis will prevent localization of peptidoglycan hydrolase and result in sensitivity to existing β-lactam antibiotics, opening up new candidates for drug repurposing.

Prediction of Certain Well-Characterized Domains of Known Functions within the PE and PPE Proteins of Mycobacteria

The PE and PPE protein family are unique to mycobacteria. Though the complete genome sequences for over 500 M. tuberculosis strains and mycobacterial species are available, few PE and PPE proteins have been structurally and functionally characterized. We have therefore used bioinformatics tools to characterize the structure and function of these proteins. We selected representative members of the PE and PPE protein family by phylogeny analysis and using structure-based sequence annotation identified ten well-characterized protein domains of known function. Some of these domains were observed to be common to all mycobacterial species and some were species specific.

tRNA is a new target for cleavage by a MazF toxin

Toxin-antitoxin (TA) systems play key roles in bacterial persistence, biofilm formation and stress responses. The MazF toxin from the Escherichia coli mazEF TA system is a sequence- and single-strand-specific endoribonuclease, and many studies have led to the proposal that MazF family members exclusively target mRNA. However, recent data indicate some MazF toxins can cleave specific sites within rRNA in concert with mRNA. In this report, we identified the repertoire of RNAs cleaved by Mycobacterium tuberculosis toxin MazF-mt9 using an RNA-seq-based approach. This analysis revealed that two tRNAs were the principal targets of MazF-mt9, and each was cleaved at a single site in either the tRNA(Pro14) D-loop or within the tRNA(Lys43) anticodon. This highly selective target discrimination occurs through recognition of not only sequence but also structural determinants. Thus, MazF-mt9 represents the only MazF family member known to target tRNA and to require RNA structure for recognition and cleavage. Interestingly, the tRNase activity of MazF-mt9 mirrors basic features of eukaryotic tRNases that also generate stable tRNA-derived fragments that can inhibit translation in response to stress. Our data also suggest a role for tRNA distinct from its canonical adapter function in translation, as cleavage of tRNAs by MazF-mt9 downregulates bacterial growth.

RD-1 encoded EspJ protein gets phosphorylated prior to affect the growth and intracellular survival of mycobacteria

Mycobacterium tuberculosis (MTB) synchronizes a number of processes and controls a series of events to subvert host defense mechanisms for the sake of residing inside macrophages. Besides these, MTB also possesses a wide range of signal enzyme systems, including eleven serine threonine protein kinases (STPKs). The present study describes STPK modulated modification in one of the hypothetical proteins of the RD1 region; EspJ (ESX-1 secretion associated protein), which is predicted to be involved in virulence of MTB. We have employed knock-out MTB, and M. bovis BCG as a surrogate strain to elaborate the consequence of the phosphorylation of EspJ. The molecular and mass spectrometric analyses in this study, confirmed EspJ as one of the substrates of STPKs. The ectopic expression of phosphoablative mutants of espJ in M. bovis BCG also articulated the effect of phosphorylation on the growth and in survival of mycobacteria. Importantly, the level of phosphorylation of EspJ also differed between pathogenic H₃₇ Rv (Rv) and non pathogenic H₃₇ Ra (Ra) strains of MTB. This further suggested that to a certain extent, the STPKs mediated phosphorylation may be accountable, in determining the growth and in intra-cellular survival of mycobacteria.

PE11 (Rv1169c) selectively alters fatty acid components of Mycobacterium smegmatis and host cell interleukin-6 level accompanied with cell death

PE/PPE family proteins, named after their conserved PE (Pro-Glu) and PPE (Pro-Pro-Glu) domains of N-terminal, are most intriguing aspects of pathologic mycobacterial genome. The roles of most members of this family remain unknown, although selected genes of this family are related to the virulence of Mycobacterium tuberculosis. In order to decipher the role of Rv1169c, the Mycobacterium smegmatis strain heterologous expressed this ORF was constructed and identified that Rv1169c was a cell wall associated protein with a novel function in modifying the cell wall fatty acids. The growth of Rv1169c expressing strain was affected under surface stress, acidic condition and antibiotics treatment. M. smegmatis expressing Rv1169c induced necrotic cell death of macrophage after infection and significantly decreased interlukin-6 production compared to controls. In general, these results underscore a proposing role of Rv1169c in virulence of M. tuberculosis, as it's role in the susceptibility of anti-mycobacteria factors caused by modified cell wall fatty acid, and the induced necrotic cell death by Rv1169c is crucial for M. tuberculosis virulence during infection.

Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement

Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3–5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.

Comparative analyses of nonpathogenic, opportunistic, and totally pathogenic mycobacteria reveal genomic and biochemical variabilities and highlight the survival attributes of Mycobacterium tuberculosis

Mycobacterial evolution involves various processes, such as genome reduction, gene cooption, and critical gene acquisition. Our comparative genome size analysis of 44 mycobacterial genomes revealed that the nonpathogenic (NP) genomes were bigger than those of opportunistic (OP) or totally pathogenic (TP) mycobacteria, with the TP genomes being smaller yet variable in size—their genomic plasticity reflected their ability to evolve and survive under various environmental conditions. From the 44 mycobacterial species, 13 species, representing TP, OP, and NP, were selected for genomic-relatedness analyses. Analysis of homologous protein-coding genes shared between Mycobacterium indicus pranii (NP), Mycobacterium intracellulare ATCC 13950 (OP), and Mycobacterium tuberculosis H37Rv (TP) revealed that 4,995 (i.e., ~95%) M. indicaus pranii proteins have homology with M. intracellulare, whereas the homologies among M. indicus pranii, M. intracellulare ATCC 13950, and M. tuberculosis H37Rv were significantly lower. A total of 4,153 (~79%) M. indicus pranii proteins and 4,093 (~79%) M. intracellulare ATCC 13950 proteins exhibited homology with the M. tuberculosis H37Rv proteome, while 3,301 (~82%) and 3,295 (~82%) M. tuberculosis H37Rv proteins showed homology with M. indicus pranii and M. intracellulare ATCC 13950 proteomes, respectively. Comparative metabolic pathway analyses of TP/OP/NP mycobacteria showed enzymatic plasticity between M. indicus pranii (NP) and M. intracellulare ATCC 13950 (OP), Mycobacterium avium 104 (OP), and M. tuberculosis H37Rv (TP). Mycobacterium tuberculosis seems to have acquired novel alternate pathways with possible roles in metabolism, host-pathogen interactions, virulence, and intracellular survival, and by implication some of these could be potential drug targets.

The complete sequence analysis of Mycobacterium indicus pranii, a novel species of Mycobacterium shown earlier to have strong immunomodulatory properties and currently in use for the treatment of leprosy, places it evolutionarily at the point of transition to pathogenicity. With the purpose of establishing the importance of M. indicus pranii in providing insight into the virulence mechanism of tuberculous and nontuberculous mycobacteria, we carried out comparative genomic and proteomic analyses of 44 mycobacterial species representing nonpathogenic (NP), opportunistic (OP), and totally pathogenic (TP) mycobacteria. Our results clearly placed M. indicus pranii as an ancestor of the M. avium complex. Analyses of comparative metabolic pathways between M. indicus pranii (NP), M. tuberculosis (TP), and M. intracellulare (OP) pointed to the presence of novel alternative pathways in M. tuberculosis with implications for pathogenesis and survival in the human host and identification of new drug targets.

Mycobacterium tuberculosis PE25/PPE41 protein complex induces necrosis in macrophages: Role in virulence and disease reactivation?

•

The Mycobacterium secreted protein PE25/PPE41 drives TNF-α secretion.

•

PE25/PPE41 protein induces necrotic cell death, but not apoptosis, in macrophages.

•

Necotic cell death induced by PE25/PPE41 is independent of TNF-α/NFκB/AP-1 pathways.

•

PE25/PPE41 possibly acts as virulence factor, by an ‘immune quorum sensing’ mechanism.

•

Necrotic cell death may help in mycobacterial dissemination and re-activation.

Necrotic cell death during TB infection is an important prerequisite for bacterial dissemination and virulence. The underlying mechanisms and the bacterial factors involved therein are not well understood. The Mycobacterium tuberculosis (M. tuberculosis) co-operonic PE25/PPE41 protein complex, similar to ESAT-6/CFP-10, belonging to the PE/PPE and ESAT-6 families of genes has co-expanded and co-evolved in the genomes of pathogenic mycobacteria. We report a novel role of this highly immunogenic PE25/PPE41 protein complex in inducing necrosis, but not apoptosis, in macrophages. We propose that these protein complexes of M. tuberculosis, secreted by similar/unique transport system (Type VII), have an important role in M. tuberculosis virulence and disease reactivation.

A novel molecular typing method of Mycobacteria based on DNA barcoding visualization

Different subtypes of Mycobacterium tuberculosis (MTB) may induce diverse severe human infections, and some of their symptoms are similar to other pathogenes, e.g. Nontuberculosis mycobacteria (NTM). So determination of mycobacterium subtypes facilitates the effective control of MTB infection and proliferation. This study exploits a novel DNA barcoding visualization method for molecular typing of 17 mycobacteria genomes published in the NCBI prokaryotic genome database. Three mycobacterium genes (Rv0279c, Rv3508 and Rv3514) from the PE/PPE family of MT Band were detected to best represent the inter-strain pathogenetic variations. An accurate and fast MTB substrain typing method was proposed based on the combination of the aforementioned three biomarker genes and the 16S rRNA gene. The protocol of establishing a bacterial substrain typing system used in this study may also be applied to the other pathogenes.

Analysis of differentially expressed proteins in late-stationary growth phase of Mycobacterium tuberculosis H37Rv

Mycobacterium tuberculosis is a causative agent of tuberculosis (TB). The ability of M. tuberculosis to be quiescent in the cell has caused the emergence of latent infection. A comprehensive proteomic analysis of M. tuberculosis H37Rv over three growth phases, namely mid-log (14-day culture), early stationary (28-day culture), and late stationary (50-day culture), was performed in order to study the change in proteome from the mid-log phase to late-stationary phase. Combination methods of two-dimensional electrophoresis (2-DE) and tandem mass spectrometry were used to generate proteome maps of M. tuberculosis at different growth phases. Ten proteins were detected differentially expressed in the late-stationary phase compared with the other two phases. These proteins were SucD, TrpD, and Rv2161c, which belong to metabolic pathway proteins; FadE5, AccD5, DesA1, and Rv1139c are proteins involved in cell wall or lipid biosynthesis, whereas TB21.7 and Rv3224 are conserved hypothetical proteins with unknown function. A surface antigen protein, DesA1, was not detectable in the late-stationary phase, although present in both log and early-stationary phases. The changes in the expression levels of these proteins were in line with the growth environment changes of the bacteria from mid-log phase to late-stationary phase. The information gathered may be valuable in the intervention against latent TB infection.

In search of actionable targets for agrigenomics and microalgal biofuel production: sequence-structural diversity studies on algal and higher plants with a focus on GPAT protein

The triacylglycerol (TAG) pathway provides several targets for genetic engineering to optimize microalgal lipid productivity. GPAT (glycerol-3-phosphate acyltransferase) is a crucial enzyme that catalyzes the initial step of TAG biosynthesis. Despite many recent biochemical studies, a comprehensive sequence-structure analysis of GPAT across diverse lipid-yielding organisms is lacking. Hence, we performed a comparative genomic analysis of plastid-located GPAT proteins from 7 microalgae and 3 higher plants species. The close evolutionary relationship observed between red algae/diatoms and green algae/plant lineages in the phylogenetic tree were further corroborated by motif and gene structure analysis. The predicted molecular weight, amino acid composition, Instability Index, and hydropathicity profile gave an overall representation of the biochemical features of GPAT protein across the species under study. Furthermore, homology models of GPAT from Chlamydomonas reinhardtii, Arabidopsis thaliana, and Glycine max provided deep insights into the protein architecture and substrate binding sites. Despite low sequence identity found between algal and plant GPATs, the developed models exhibited strikingly conserved topology consisting of 14α helices and 9β sheets arranged in two domains. However, subtle variations in amino acids of fatty acyl binding site were identified that might influence the substrate selectivity of GPAT. Together, the results will provide useful resources to understand the functional and evolutionary relationship of GPAT and potentially benefit in development of engineered enzyme for augmenting algal biofuel production.