A postgenomic method for predicting essential genes at subsaturation levels of mutagenesis: application to Mycobacterium tuberculosis

Mycobacterium marinum as a model for understanding principles of mycobacterial pathogenesis

Mycobacterium marinum is a fish pathogen that has become a powerful and well-established model that has accelerated our understanding of the mechanisms of mycobacterial disease. M. marinum is a versatile surrogate for understanding the closely related human pathogen M. tuberculosis, which causes tuberculosis in humans. M. marinum has defined key mechanisms of pathogenesis, both shared with M. tuberculosis and unique to this species. In this review, we discuss the discovery of M. marinum as an occasional human pathogen, the shared aspects of pathogenesis with M. tuberculosis, and how M. marinum has been exploited as a model to define the molecular mechanisms of mycobacterial pathogenesis across several phases of infection.

Domain architecture of the Mycobacterium tuberculosis MabR (Rv2242), a member of the PucR transcription factor family

MabR (Rv2242), a PucR-type transcription factor, plays a crucial role in regulating mycolic acid biosynthesis in Mycobacterium tuberculosis. To understand its regulatory mechanisms, we determined the crystal structures of its N-terminal and C-terminal domains. The N-terminal domain adopts a globin-like fold, while the C-terminal domain comprises an α/β GGDEF domain and an all-α effector domain with a helix-turn-helix DNA-binding motif. This unique domain combination is specific to Actinomycetes. Biochemical and computational studies suggest that full-length MabR forms both dimeric and tetrameric assemblies in solution. Structural analysis revealed two distinct dimerization interfaces within the N- and C-terminal domains, further supporting a tetrameric organization. These findings provide valuable insights into the domain architecture, oligomeric state, and potential regulatory mechanisms of MabR.

Identification of a proteolysis regulator for an essential enzyme in Mycobacterium tuberculosis

In Mycobacterium tuberculosis proteins that are post-translationally modified with Pup, a prokaryotic ubiquitin-like protein, can be degraded by proteasomes. While pupylation is reversible, mechanisms regulating substrate specificity have not been identified. Here, we identify the first depupylation regulators: CoaX, a pseudokinase, and pantothenate, an essential, central metabolite. In a Δ coaX mutant, pantothenate synthesis enzymes were more abundant, including PanB, a substrate of the Pup-proteasome system. Media supplementation with pantothenate decreased PanB levels in a coaX and Pup-proteasome-dependent manner. In vitro , CoaX accelerated depupylation of Pup∼PanB, while addition of pantothenate inhibited this reaction. Collectively, we propose CoaX contributes to proteasomal degradation of PanB by modulating depupylation of Pup∼PanB in response to pantothenate levels.

One Sentence Summary

A pseudo-pantothenate kinase regulates proteasomal degradation of a pantothenate synthesis enzyme in M. tuberculosis .

Identification of genes associated with persistence in Mycobacterium smegmatis

The prevalence of bacterial persisters is related to their phenotypic diversity and is responsible for the relapse of chronic infections. Tolerance to antibiotic therapy is the hallmark of bacterial persistence. In this study, we have screened a transposon library of Mycobacterium smegmatis mc2155 strain using antibiotic tolerance, survival in mouse macrophages, and biofilm-forming ability of the mutants. Out of 10 thousand clones screened, we selected ten mutants defective in all the three phenotypes. Six mutants showed significantly lower persister abundance under different stress conditions. Insertions in three genes belonging to the pathways of oxidative phosphorylation msmeg_3233 (cydA), biotin metabolism msmeg_3194 (bioB), and oxidative metabolism msmeg_0719, a flavoprotein monooxygenase, significantly reduced the number of live cells, suggesting their role in pathways promoting long-term survival. Another group that displayed a moderate reduction in CFU included a glycosyltransferase, msmeg_0392, a hydrogenase subunit, msmeg_2263 (hybC), and a DNA binding protein, msmeg_2211. The study has revealed potential candidates likely to facilitate the long-term survival of M. smegmatis. The findings offer new targets to develop antibiotics against persisters. Further, investigating the corresponding genes in M. tuberculosis may provide valuable leads in improving the treatment of chronic and persistent tuberculosis infections.

The major role of Listeria monocytogenes folic acid metabolism during infection is the generation of N-formylmethionine

IMPORTANCE

Folic acid is an essential vitamin for bacteria, plants, and animals. The lack of folic acid leads to various consequences such as a shortage of amino acids and nucleotides that are fundamental building blocks for life. Though antifolate drugs are widely used for antimicrobial treatments, the underlying mechanism of bacterial folate deficiency during infection is unclear. This study compares the requirements of different folic acid end-products during the infection of Listeria monocytogenes, a facultative intracellular pathogen of animals and humans. The results reveal the critical importance of N-formylmethionine, the amino acid used by bacteria to initiate protein synthesis. This work extends the current understanding of folic acid metabolism in pathogens and potentially provides new insights into antifolate drug development in the future.

Role of succinyl substituents in the mannose-capping of lipoarabinomannan and control of inflammation in Mycobacterium tuberculosis infection

The covalent modification of bacterial (lipo)polysaccharides with discrete substituents may impact their biosynthesis, export and/or biological activity. Whether mycobacteria use a similar strategy to control the biogenesis of its cell envelope polysaccharides and modulate their interaction with the host during infection is unknown despite the report of a number of tailoring substituents modifying the structure of these glycans. Here, we show that discrete succinyl substituents strategically positioned on Mycobacterium tuberculosis (Mtb) lipoarabinomannan govern the mannose-capping of this lipoglycan and, thus, much of the biological activity of the entire molecule. We further show that the absence of succinyl substituents on the two main cell envelope glycans of Mtb, arabinogalactan and lipoarabinomannan, leads to a significant increase of pro-inflammatory cytokines and chemokines in infected murine and human macrophages. Collectively, our results validate polysaccharide succinylation as a critical mechanism by which Mtb controls inflammation.

Inhibitors of riboflavin biosynthetic pathway enzymes as potential antibacterial drugs

Multiple drug resistance is the main obstacle in the treatment of bacterial diseases. Resistance against antibiotics demands the exploration of new antimicrobial drug targets. A variety of in silico and genetic approaches show that the enzymes of the riboflavin biosynthetic pathway are crucial for the survival of bacteria. This pathway is absent in humans thus enzymes of the riboflavin biosynthetic pathway are emerging drug targets for resistant pathogenic bacterial strains. Exploring the structural details, their mechanism of action, intermediate elucidation, and interaction analysis would help in designing suitable inhibitors of these enzymes. The riboflavin biosynthetic pathway consists of seven distinct enzymes, namely, 3,4-dihydroxy-2-butanone 4-phosphate synthase, GTP cyclohydrolase II, pyrimidine deaminase/reductase, phosphatase, lumazine synthase, and riboflavin synthase. The present review summarizes the research work that has been carried out on these enzymes in terms of their structures, active site architectures, and molecular mechanism of catalysis. This review also walks through small molecule inhibitors that have been developed against several of these enzymes.

The human proton pump inhibitors inhibit Mycobacterium tuberculosis rifampicin efflux and macrophage-induced rifampicin tolerance

Tuberculosis treatment requires months-long combination chemotherapy with multiple drugs, with shorter treatments leading to relapses. A major impediment to shortening treatment is that Mycobacterium tuberculosis becomes tolerant to the administered drugs, starting early after infection and within days of infecting macrophages. Multiple lines of evidence suggest that macrophage-induced drug tolerance is mediated by mycobacterial drug efflux pumps. Here, using assays to directly measure drug efflux, we find that M. tuberculosis transports the first-line antitubercular drug rifampicin through a proton gradient-dependent mechanism. We show that verapamil, a known efflux pump inhibitor, which inhibits macrophage-induced rifampicin tolerance, also inhibits M.tuberculosis rifampicin efflux. As with macrophage-induced tolerance, the calcium channel-inhibiting property of verapamil is not required for its inhibition of rifampicin efflux. By testing verapamil analogs, we show that verapamil directly inhibits M. tuberculosis drug efflux pumps through its human P-glycoprotein (PGP)-like inhibitory activity. Screening commonly used drugs with incidental PGP inhibitory activity, we find many inhibit rifampicin efflux, including the proton pump inhibitors (PPIs) such as omeprazole. Like verapamil, the PPIs inhibit macrophage-induced rifampicin tolerance as well as intramacrophage growth, which has also been linked to mycobacterial efflux pump activity. Our assays provide a facile screening platform for M. tuberculosis efflux pump inhibitors that inhibit in vivo drug tolerance and growth.

Identification of 2-Aryl-Quinolone Inhibitors of Cytochrome bd and Chemical Validation of Combination Strategies for Respiratory Inhibitors against Mycobacterium tuberculosis

Mycobacterium tuberculosis cytochrome bd quinol oxidase (cyt bd), the alternative terminal oxidase of the respiratory chain, has been identified as playing a key role during chronic infection and presents a putative target for the development of novel antitubercular agents. Here, we report confirmation of successful heterologous expression of M. tuberculosis cytochrome bd. The heterologous M. tuberculosis cytochrome bd expression system was used to identify a chemical series of inhibitors based on the 2-aryl-quinolone pharmacophore. Cytochrome bd inhibitors displayed modest efficacy in M. tuberculosis growth suppression assays together with a bacteriostatic phenotype in time-kill curve assays. Significantly, however, inhibitor combinations containing our front-runner cyt bd inhibitor CK-2-63 with either cyt bcc-aa3 inhibitors (e.g., Q203) and/or adenosine triphosphate (ATP) synthase inhibitors (e.g., bedaquiline) displayed enhanced efficacy with respect to the reduction of mycobacterium oxygen consumption, growth suppression, and in vitro sterilization kinetics. In vivo combinations of Q203 and CK-2-63 resulted in a modest lowering of lung burden compared to treatment with Q203 alone. The reduced efficacy in the in vivo experiments compared to in vitro experiments was shown to be a result of high plasma protein binding and a low unbound drug exposure at the target site. While further development is required to improve the tractability of cyt bd inhibitors for clinical evaluation, these data support the approach of using small-molecule inhibitors to target multiple components of the branched respiratory chain of M. tuberculosis as a combination strategy to improve therapeutic and pharmacokinetic/pharmacodynamic (PK/PD) indices related to efficacy.

Biosynthesis of Isonitrile Lipopeptide Metallophores from Pathogenic Mycobacteria

Isonitrile lipopeptides (INLPs) are known to be related to the virulence of pathogenic mycobacteria by mediating metal transport, but their biosynthesis remains obscure. In this work, we use in vitro biochemical assays, site-directed mutagenesis, chemical synthesis, and spectroscopy techniques to scrutinize the activity of core enzymes required for INLP biosynthesis in mycobacteria. Compared to environmental Streptomyces, pathogenic Mycobacterium employ a similar chemical logic and enzymatic machinery in INLP biosynthesis, differing mainly in the fatty-acyl chain length, which is controlled by multiple enzymes in the pathway. Our in-depth study on the non-heme iron(II) and α-ketoglutarate-dependent dioxygenase for isonitrile generation, including Rv0097 from Mycobacterium tuberculosis (Mtb), demonstrates that it recognizes a free-standing small molecule substrate, different from the recent hypothesis that a carrier protein is required for Rv0097 in Mtb. A key residue in Rv0097 is further identified to dictate the varied fatty-acyl chain length specificity between Streptomyces and Mycobacterium.

Mycobacterium tuberculosis senses host Interferon-γ via the membrane protein MmpL10

Mycobacterium tuberculosis (Mtb) is one of the most successful human pathogens. Several cytokines are known to increase virulence of bacterial pathogens, leading us to investigate whether Interferon-γ (IFN-γ), a central regulator of the immune defense against Mtb, has a direct effect on the bacteria. We found that recombinant and T-cell derived IFN-γ rapidly induced a dose-dependent increase in the oxygen consumption rate (OCR) of Mtb, consistent with increased bacterial respiration. This was not observed in attenuated Bacillus Calmette-Guérin (BCG), and did not occur for other cytokines tested, including TNF-α. IFN-γ binds to the cell surface of intact Mtb, but not BCG. Mass spectrometry identified mycobacterial membrane protein large 10 (MmpL10) as the transmembrane binding partner of IFN-γ, supported by molecular modelling studies. IFN-γ binding and the OCR response was absent in Mtb Δmmpl10 strain and restored by complementation with wildtype mmpl10. RNA-sequencing and RT-PCR of Mtb exposed to IFN-γ revealed a distinct transcriptional profile, including genes involved in virulence. In a 3D granuloma model, IFN-γ promoted Mtb growth, which was lost in the Mtb Δmmpl10 strain and restored by complementation, supporting the involvement of MmpL10 in the response to IFN-γ. Finally, IFN-γ addition resulted in sterilization of Mtb cultures treated with isoniazid, indicating clearance of phenotypically resistant bacteria that persist in the presence of drug alone. Together our data are the first description of a mechanism allowing Mtb to respond to host immune activation that may be important in the immunopathogenesis of TB and have use in novel eradication strategies.

Identification of essential genes in Mycobacterium avium subsp. paratuberculosis genome for persistence in dairy calves

To cause disease Mycobacterium avium subsp. paratuberculosis needs to enter mammalian cells, arrest phagosomal maturation and manipulate the host immune system. The genetic basis of the bacterial capacity to achieve these outcomes remains largely unknown. Identifying these genes would allow us to gain a deeper understanding of MAP's pathogenesis and potentially develop a live attenuated Johne's disease vaccine by knocking out these genes. MAP genes demonstrated to be essential for colonization in the natural host, ruminants, are unknown. Genome-wide transposon mutagenesis and high-throughput sequencing were combined to evaluate the essentiality of each coding region in the bacterial genome to survive in dairy calves. A saturated library of 3,852 MAP Tn mutants, with insertions in 56% of TA sites, interrupting 88% of genes, was created using a MycoMarT7 phagemid containing a mariner transposon. Six calves were inoculated with a high dose of a library of MAP mutants, 1011 CFUs, (input) at 2 weeks of age. Following 2 months of incubation, MAP cells were isolated from the ileum, jejunum, and their associated lymph nodes of calves, resulting in approximately 100,000 colonies grown on solid media across 6 animals (output). Targeted next-generation sequencing was used to identify the disrupted genes in all the mutants in the input pool and the output pool recovered from the tissues to identify in vivo essential genes. Statistical analysis for the determination of essential genes was performed by a Hidden Markov Model (HMM), categorizing genes into essential genes that are devoid of insertions and growth-defect genes whose disruption impairs the growth of the organism. Sequence analysis identified 430 in vivo essential and 260 in vivo growth-defect genes. Gene ontology enrichment analysis of the in vivo essential and growth-defect genes with the highest reduction in the tissues revealed a high representation of genes involved in metabolism and respiration, cell wall and cell processing, virulence, and information pathway processes. This study has systematically identified essential genes for the growth and persistence of MAP in the natural host body.

Iron–Sulfur Clusters toward Stresses: Implication for Understanding and Fighting Tuberculosis

Tuberculosis (TB) remains the leading cause of death due to a single pathogen, accounting for 1.5 million deaths annually on the global level. Mycobacterium tuberculosis, the causative agent of TB, is persistently exposed to stresses such as reactive oxygen species (ROS), reactive nitrogen species (RNS), acidic conditions, starvation, and hypoxic conditions, all contributing toward inhibiting bacterial proliferation and survival. Iron–sulfur (Fe-S) clusters, which are among the most ancient protein prosthetic groups, are good targets for ROS and RNS, and are susceptible to Fe starvation. Mtb holds Fe-S containing proteins involved in essential biological process for Mtb. Fe-S cluster assembly is achieved via complex protein machineries. Many organisms contain several Fe-S assembly systems, while the SUF system is the only one in some pathogens such as Mtb. The essentiality of the SUF machinery and its functionality under the stress conditions encountered by Mtb underlines how it constitutes an attractive target for the development of novel anti-TB.

RNase HI Depletion Strongly Potentiates Cell Killing by Rifampicin in Mycobacteria

Multidrug-resistant (MDR) tuberculosis (TB) is defined by the resistance of Mycobacterium tuberculosis, the causative organism, to the first-line antibiotics rifampicin and isoniazid. Mitigating or reversing resistance to these drugs offers a means of preserving and extending their use in TB treatment. R-loops are RNA/DNA hybrids that are formed in the genome during transcription, and they can be lethal to the cell if not resolved. RNase HI is an enzyme that removes R-loops, and this activity is essential in M. tuberculosis: knockouts of rnhC, the gene encoding RNase HI, are nonviable. This essentiality makes it a candidate target for the development of new antibiotics. In the model organism Mycolicibacterium smegmatis, RNase HI activity is provided by two enzymes, RnhA and RnhC. We show that the partial depletion of RNase HI activity in M. smegmatis, by knocking out either of the genes encoding RnhA or RnhC, led to the accumulation of R-loops. The sensitivity of the knockout strains to the antibiotics moxifloxacin, streptomycin, and rifampicin was increased, the latter by a striking near 100-fold. We also show that R-loop accumulation accompanies partial transcriptional inhibition, suggesting a mechanistic basis for the synergy between RNase HI depletion and rifampicin. A model of how transcriptional inhibition can potentiate R-loop accumulation is presented. Finally, we identified four small molecules that inhibit recombinant RnhC activity and that also potentiated rifampicin activity in whole-cell assays against M. tuberculosis, supporting an on-target mode of action and providing the first step in developing a new class of antimycobacterial drug.

Evaluating the Performance of PPE44, HSPX, ESAT-6 and CFP-10 Factors in Tuberculosis Subunit Vaccines

Mycobacterium tuberculosis (M. tuberculosis) is an intracellular pathogen causing long-term infection in humans that mainly attacks macrophages and can escape from the immune system with the various mechanisms. The only FDA-approved vaccine against M. tuberculosis (MTB) is Mycobacterium bovis bacillus Calmette-Guérin (BCG). The protection of this vaccine typically lasts 10–15 years. Due to the increasing number of people becoming ill with MTB each year worldwide, the need to develop a new effective treatment against the disease has been increased. During the past two decades, the research budget for TB vaccine has quadrupled to over half a billion dollars. Most of these research projects were based on amplifying and stimulating the response of T-cells and developing the subunit vaccines. Additionally, these studies have demonstrated that secretory and immunogenic proteins of MTB play a key role in the pathogenesis of the bacteria. Therefore, these proteins were used to develop the new subunit vaccines. In this review, based on the use of these proteins in the successful new subunit vaccines, the PPE44, HSPX, CFP-10 and ESAT-6 antigens were selected and the role of these antigens in designing and developing new subunit vaccines against TB and for the prevention of TB were investigated.

Glby, Encoded by MAB_3167c, Is Required for In Vivo Growth of Mycobacteroides abscessus and Exhibits Mild β-Lactamase Activity

Mycobacteroides abscessus (Mab; also known as Mycobacterium abscessus) is an emerging opportunistic pathogen. Patients with structural lung conditions such as bronchiectasis, cystic fibrosis, and chronic obstructive pulmonary disease are at high risk of developing pulmonary Mab disease. This disease is often chronic as the current treatment regimens are sub-efficacious. Here, we characterize the phenotype of a Mab strain lacking the MAB_3167c locus, which encodes a protein hereafter referred to as Glby. We demonstrate that the loss of Glby impairs normal planktonic growth in liquid broth, results in longer average cell length, and a melding of surfaces between cells. Glby also exhibits a mild β-lactamase activity. We also present evidence that amino acid substitutions that potentially alter Glby function are not favored. Lastly, we demonstrate that, in a mouse model of pulmonary Mab infection, the mutant lacking Glby was unable to proliferate, gradually cleared, and was undetectable after 3 weeks. These data suggest that an agent that inhibits Glby in vivo may be an efficacious treatment against Mab disease. IMPORTANCE Mycobacteroides abscessus can cause chronic pulmonary infections requiring administration of multiple antibiotics, still resulting in a low cure rate. The incidence of M. abscessus disease is increasing in the United States and the developed regions of the world. We show for the first time that a protein, Glby, affects growth of this bacterium. Using a mouse model of lung M. abscessus disease, we demonstrate that Glby is required for this bacterium to cause disease.

A Feedback Regulatory Loop Containing McdR and WhiB2 Controls Cell Division and DNA Repair in Mycobacteria

Cell division must be coordinated with DNA repair, which is strictly regulated in response to different drugs and environmental stresses in bacteria. However, the mechanisms by which mycobacteria orchestrate these two processes remain largely uncharacterized. Here, we report a regulatory loop between two essential mycobacterial regulators, McdR (Rv1830) and WhiB2, in coordinating the processes of cell division and DNA repair. McdR inhibits cell division-associated whiB2 expression by binding to the AATnACAnnnnTGTnATT motif in the promoter region. Furthermore, McdR overexpression simultaneously activates imuAB and dnaE2 expression to promote error-prone DNA repair, which facilitates genetic adaptation to stress conditions. Through a feedback mechanism, WhiB2 activates mcdR expression by binding to the cGACACGc motif in the promoter region. Importantly, analyses of mutations in clinical Mycobacterium tuberculosis strains indicate that disruption of this McdR-WhiB2 feedback regulatory loop influences expression of both cell growth- and DNA repair-associated genes, which further supports the contribution of McdR-WhiB2 regulatory loop in regulating mycobacterial cell growth and drug resistance. This highly conserved feedback regulatory loop provides fresh insight into the link between mycobacterial cell growth control and stress responses. IMPORTANCE Drug-resistant M. tuberculosis poses a threat to the control and prevention of tuberculosis (TB) worldwide. Thus, there is a need to identify the mechanisms enabling M. tuberculosis to adapt and grow under drug-induced stress. Rv1830 has been shown to be associated with drug resistance in M. tuberculosis, but its mechanisms have not yet been elucidated. Here, we reveal a regulatory role of Rv1830, which coordinates cell division and DNA repair in mycobacteria, and rename it McdR (mycobacterial cell division regulator). An increase in McdR levels represses the expression of cell division-associated whiB2 but activates the DNA repair-associated, error-prone enzymes ImuA/B and DnaE2, which in turn facilitates adaptation to stress responses and drug resistance. Furthermore, WhiB2 activates the transcription of mcdR to form a conserved regulatory loop. These data provide new insights into the mechanisms controlling mycobacterial cell growth and stress responses.

Reprogramming Mycobacterium tuberculosis CRISPR System for Gene Editing and Genome-wide RNA Interference Screening

Mycobacterium tuberculosis is the causative agent of tuberculosis (TB), which is still the leading cause of mortality from a single infectious disease worldwide. The development of novel anti-TB drugs and vaccines is severely hampered by the complicated and time-consuming genetic manipulation techniques for M. tuberculosis. Here, we harnessed an endogenous type III-A CRISPR/Cas10 system of M. tuberculosis for efficient gene editing and RNA interference (RNAi). This simple and easy method only needs to transform a single mini-CRISPR array plasmid, thus avoiding the introduction of exogenous protein and minimizing proteotoxicity. We demonstrated that M. tuberculosis genes can be efficiently and specifically knocked in/out by this system as confirmed by DNA high-throughput sequencing. This system was further applied to single- and multiple-gene RNAi. Moreover, we successfully performed genome-wide RNAi screening to identify M. tuberculosis genes regulating in vitro and intracellular growth. This system can be extensively used for exploring the functional genomics of M. tuberculosis and facilitate the development of novel anti-TB drugs and vaccines.

Essential Protein Recognition via Community Significance

Essential protein plays a vital role in understanding the cellular life. With the advance in high-throughput technologies, a number of protein-protein interaction (PPI) networks have been constructed such that essential proteins can be identified from a system biology perspective. Although a series of network-based essential protein discovery methods have been proposed, these existing methods still have some drawbacks. Recently, it has been shown that the significance-based method SigEP is promising on overcoming the defects that are inherent in currently available essential protein identification methods. However, the SigEP method is developed under the unrealistic Erdös-Rényi (E-R) model and its time complexity is very high. Hence, we propose a new significance-based essential protein recognition method named EPCS in which the essential protein discovery problem is formulated as a community significance testing problem. Experimental results on four PPI networks show that EPCS performs better than nine state-of-the-art essential protein identification methods and the only significance-based essential protein identification method SigEP.

Prediction of Essential Genes in Comparison States Using Machine Learning

Identifying essential genes in comparison states (EGS) is vital to understanding cell differentiation, performing drug discovery, and identifying disease causes. Here, we present a machine learning method termed Prediction of Essential Genes in Comparison States (PreEGS). To capture the alteration of the network in comparison states, PreEGS extracts topological and gene expression features of each gene in a five-dimensional vector. PreEGS also recruits a positive sample expansion method to address the problem of unbalanced positive and negative samples, which is often encountered in practical applications. Different classifiers are applied to the simulated datasets, and the PreEGS based on the random forests model (PreEGSRF) was chosen for optimal performance. PreEGSRF was then compared with six other methods, including three machine learning methods, to predict EGS in a specific state. On real datasets with four gene regulatory networks, PreEGSRF predicted five essential genes related to leukemia and five enriched KEGG pathways. Four of the predicted essential genes and all predicted pathways were consistent with previous studies and highly correlated with leukemia. With high prediction accuracy and generalization ability, PreEGSRF is broadly applicable for the discovery of disease-causing genes, driver genes for cell fate decisions, and complex biomarkers of biological systems.

Therapeutic target identification via differential genome analysis of antibiotic resistant Shigella sonnei and inhibitor evaluation against a selected drug target

Shigella sonnei has been implicated in bloody diarrhea (accompanied by abdominal pain and fever) and is an emerging pathogen of concern, especially in developing countries. The major means of transmission is the fecal-oral route while sexual transmission has also been reported. In children, the impact might be stunted growth due to life-threatening illness. Resistance has been reported in this species for several types of antibiotics. In this study, we retrieved the antibiotic-resistant labeled whole genome sequences of the species from the PATRIC database and performed a pan-genome analysis to filter out core genes. Antibiotic resistance was studied in the core, accessory and unique genome. Core genes were utilized as seed substance for essentiality analysis and drug candidate assignment. Product of the gene aroG, i.e. chorismate biosynthetic process 3-deoxy-7-phosphoheptulonate synthase enzyme, responsible for aromatic amino acid family biosynthetic process, was taken for further downstream processing. Natural product libraries of flavonoids (n = 178), ZINC database derived inhibitor compounds of the 3-deoxy-7-phosphoheptulonate synthase enzyme (n = 112), and streptomycin compounds (n = 737) were docked to find out potent inhibitors, followed by dynamics simulation of 50 ns each for top compounds.. Physicochemical and ADMET profiling of the top compounds was done to analyze their safety for consumption. We propose that the top compounds: Phytoene from Streptomycin library and ZINC000036444158 (synonym:1,16-bis[(dihydroxyphosphinyl)oxy]hexadecane) from 3-deoxy-7-phosphoheptulonate synthase inhibitor library of ZINC database (and used as a control in this study) should be tested in vitro against Shigella sonnei, to fully determine their efficacy. This could add to the drying pipeline of potent drug molecules against emerging pathogens.

Genome-Wide Essentiality Analysis of Mycobacterium abscessus by Saturated Transposon Mutagenesis and Deep Sequencing

Mycobacterium abscessus is an emerging opportunistic human pathogen that naturally resists most major classes of antibiotics, making infections difficult to treat. Thus far, little is known about M. abscessus physiology, pathogenesis, and drug resistance. Genome-wide analyses have comprehensively catalogued genes with essential functions in Mycobacterium tuberculosis and Mycobacterium avium subsp. hominissuis (here, M. avium) but not in M. abscessus. By optimizing transduction conditions, we achieved full saturation of TA insertion sites with Himar1 transposon mutagenesis in the M. abscessus ATCC 19977^T genome, as confirmed by deep sequencing prior to essentiality analyses of annotated genes and other genomic features. The overall densities of inserted TA sites (85.7%), unoccupied TA sites (14.3%), and nonpermissive TA sites (8.1%) were similar to results in M. tuberculosis and M. avium. Of the 4,920 annotated genes, 326 were identified as essential, 269 (83%) of which have mutual homology with essential M. tuberculosis genes, while 39 (12%) are homologous to genes that are not essential in M. tuberculosis and M. avium, and 11 (3.4%) only have homologs in M. avium. Interestingly, 7 (2.1%) essential M. abscessus genes have no homologs in either M. tuberculosis or M. avium, two of which were found in phage-like elements. Most essential genes are involved in DNA replication, RNA transcription and translation, and posttranslational events to synthesize important macromolecules. Some essential genes may be involved in M. abscessus pathogenesis and antibiotics response, including certain essential tRNAs and new short open reading frames. Our findings will help to pave the way for better understanding of M. abscessus and benefit development of novel bactericidal drugs against M. abscessus. IMPORTANCE Limited knowledge regarding Mycobacterium abscessus pathogenesis and intrinsic resistance to most classes of antibiotics is a major obstacle to developing more effective strategies to prevent and mitigate disease. Using optimized procedures for Himar1 transposon mutagenesis and deep sequencing, we performed a comprehensive analysis to identify M. abscessus genetic elements essential for in vitro growth and compare them to similar data sets for M. tuberculosis and M. avium subsp. hominissuis. Most essential M. abscessus genes have mutual homology with essential M. tuberculosis genes, providing a foundation for leveraging available knowledge from M. tuberculosis to develop more effective drugs and other interventions against M. abscessus. A small number of essential genes unique to M. abscessus deserve further attention to gain insights into what makes M. abscessus different from other mycobacteria. The essential genes and other genomic features such as short open reading frames and noncoding RNA identified here will provide useful information for future study of M. abscessus pathogenicity and new drug development.

Biostatistical prediction of genes essential for growth of Anaplasma phagocytophilum in a human promyelocytic cell line using a random transposon mutant library

Anaplasma phagocytophilum (Ap), agent of human anaplasmosis, is an intracellular bacterium that causes the second most common tick-borne illness in North America. To address the lack of a genetic system for these pathogens, we used random Himar1 transposon mutagenesis to generate a library of Ap mutants capable of replicating in human promyelocytes (HL-60 cells). Illumina sequencing identified 1195 non-randomly distributed insertions. As the density of mutants was non-saturating, genes without insertions were either essential for Ap, or spared randomly. To resolve this question, we applied a biostatistical method for prediction of essential genes. Since the chances that a transposon was inserted into genomic TA dinucleotide sites should be the same for all loci, we used a Markov chain Monte Carlo model to estimate the probability that a non-mutated gene was essential for Ap. Predicted essential genes included those coding for structural ribosomal proteins, enzymes involved in metabolism, components of the type IV secretion system, antioxidant defense molecules and hypothetical proteins. We have used an in silico post-genomic approach to predict genes with high probability of being essential for replication of Ap in HL-60 cells. These results will help target genes to investigate their role in the pathogenesis of human anaplasmosis.

The ATP-Binding Cassette (ABC) Transport Systems in Mycobacterium tuberculosis: Structure, Function, and Possible Targets for Therapeutics

Simple Summary

Mycobacterium tuberculosis is a bacterium of great medical importance because it causes tuberculosis, a disease that affects millions of people worldwide. Two important features are related to this bacterium: its ability to infect and survive inside the host, minimizing the immune response, and the burden of clinical isolates that are highly resistant to antibiotics treatment. These two phenomena are directly affected by cell envelope proteins, such as proteins from the ATP-Binding Cassette (ABC transporters) superfamily. In this review, we have compiled information on all the M. tuberculosis ABC transporters described so far, both from a functional and structural point of view, and show their relevance for the bacillus and the potential targets for studies aiming to control the microorganism and structural features.

Abstract

Mycobacterium tuberculosis is the etiological agent of tuberculosis (TB), a disease that affects millions of people in the world and that is associated with several human diseases. The bacillus is highly adapted to infect and survive inside the host, mainly because of its cellular envelope plasticity, which can be modulated to adapt to an unfriendly host environment; to manipulate the host immune response; and to resist therapeutic treatment, increasing in this way the drug resistance of TB. The superfamily of ATP-Binding Cassette (ABC) transporters are integral membrane proteins that include both importers and exporters. Both types share a similar structural organization, yet only importers have a periplasmic substrate-binding domain, which is essential for substrate uptake and transport. ABC transporter-type importers play an important role in the bacillus physiology through the transport of several substrates that will interfere with nutrition, pathogenesis, and virulence. Equally relevant, exporters have been involved in cell detoxification, nutrient recycling, and antibiotics and drug efflux, largely affecting the survival and development of multiple drug-resistant strains. Here, we review known ABC transporters from M. tuberculosis, with particular focus on the diversity of their structural features and relevance in infection and drug resistance.

Predicting essential genes of 41 prokaryotes by a semi-supervised method

Essential genes are vitally important to the survival and reproduction of organisms. Many machine learning methods have been widely employed to predict essential genes and have obtained satisfactory results. However, most of these methods are supervised methods and may not obtain the desired result when the labeled data are insufficient. In this paper, we proposed a learning with local and global consistency (LGC) method-based classifier, which was employed to predict the essential genes of 41 prokaryotes. LGC is a graph-based semi-supervised learning method that can construct a prediction model using finite label and constraint information. The performance of the proposed classifier was evaluated by employing intra-organism prediction and leave-one-species-out validation. The average AUC value of 41 organisms in intra-organisms prediction was 0.723 when the labeled sample ratio was 0.5. The results of this study indicate that the proposed method can achieve acceptable prediction performance with limited labeled data. Additionally, the results demonstrate that this method has good universality.

Functional Disassociation Between the Protein Domains of MSMEG_4305 of Mycolicibacterium smegmatis (Mycobacterium smegmatis) in vivo

MSMEG_4305 is a two-domain protein of Mycolicibacterium smegmatis (Mycobacterium smegmatis) (Mycolicibacterium smegmatis). The N-terminal domain of MSMEG_4305 encodes an RNase H type I. The C-terminal domain is a presumed CobC, predicted to be involved in the aerobic synthesis of vitamin B12. Both domains reach their maximum at distinct pH, approximately 8.5 and 4.5, respectively. The presence of the CobC domain influenced RNase activity in vitro in homolog Rv2228c. Here, we analyzed the role of MSMEG_4305 in vitamin B12 synthesis and the functional association between both domains in vivo in M. smegmatis. We used knock-out mutant of M. smegmatis, deficient in MSMEG_4305. Whole-cell lysates of the mutants strain contained a lower concentration of vitamin B12, as it determined with immunoenzimatic assay. We observed growth deficits, related to vitamin B12 production, on media containing sulfamethazine and propionate. Removal of the CobC domain of MSMEG_4305 in ΔrnhA background hardly affected the growth rate of M. smegmatis in vivo. The strain carrying truncation showed no fitness deficit in the competitive assay and it did not show increased level of RNA/DNA hybrids in its genome. We show that homologs of MSMEG_4305 are present only in the Actinomycetales phylogenetic branch (according to the old classification system). The domains of MSMEG_4305 homologs accumulate mutations at a different rate, while the linker region is highly variable. We conclude that MSMEG_4305 is a multidomain protein that most probably was fixed in the phylogenetic tree of life due to genetic drift.

Development of a Counterselectable Transposon To Create Markerless Knockouts from an 18,432-Clone Ordered Mycobacterium bovis Bacillus Calmette-Guérin Mutant Resource

Mutant resources are essential to improve our understanding of the biology of slow-growing mycobacteria, which include the causative agents of tuberculosis in various species, including humans. The generation of deletion mutants in slow-growing mycobacteria in a gene-by-gene approach in order to make genome-wide ordered mutant resources is still a laborious and costly approach, despite the recent development of improved methods. On the other hand, transposon mutagenesis in combination with Cartesian pooling-coordinate sequencing (CP-CSeq) allows the creation of large archived Mycobacterium transposon insertion libraries. However, such mutants contain selection marker genes with a risk of polar gene effects, which are undesired both for research and for use of these mutants as live attenuated vaccines. In this paper, a derivative of the Himar1 transposon is described which allows the generation of clean, markerless knockouts from archived transposon libraries. By incorporating FRT sites for FlpE/FRT-mediated recombination and I-SceI sites for ISceIM-based transposon removal, we enable two thoroughly experimentally validated possibilities to create unmarked mutants from such marked transposon mutants. The FRT approach is highly efficient but leaves an FRT scar in the genome, whereas the I-SceI-mediated approach can create mutants without any heterologous DNA in the genome. The combined use of CP-CSeq and this optimized transposon was applied in the BCG Danish 1331 vaccine strain (WHO reference 07/270), creating the largest ordered, characterized resource of mutants in a member of the Mycobacterium tuberculosis complex (18,432 clones, mutating 83% of the nonessential M. tuberculosis homologues), from which markerless knockouts can be easily generated.IMPORTANCE While speeding up research for many fields of biology (e.g., yeast, plant, and Caenorhabditis elegans), genome-wide ordered mutant collections are still elusive in mycobacterial research. We developed methods to generate such resources in a time- and cost-effective manner and developed a newly engineered transposon from which unmarked mutants can be efficiently generated. Our library in the WHO reference vaccine strain of Mycobacterium bovis BCG Danish targets 83% of all nonessential genes and was made publicly available via the BCCM/ITM Mycobacteria Collection. This resource will speed up Mycobacterium research (e.g., drug resistance research and vaccine development) and paves the way to similar genome-wide mutant collections in other strains of the Mycobacterium tuberculosis complex. The stretch to a full collection of mutants in all nonessential genes is now much shorter, with just 17% remaining genes to be targeted using gene-by-gene approaches, for which highly effective methods have recently also been described.

Comparative Analysis, Structural Insights, and Substrate/Drug Interaction of CYP128A1 in Mycobacterium tuberculosis

Cytochrome P450 monooxygenases (CYPs/P450s) are well known for their role in organisms' primary and secondary metabolism. Among 20 P450s of the tuberculosis-causing Mycobacterium tuberculosis H37Rv, CYP128A1 is particularly important owing to its involvement in synthesizing electron transport molecules such as menaquinone-9 (MK9). This study employs different in silico approaches to understand CYP128 P450 family's distribution and structural aspects. Genome data-mining of 4250 mycobacterial species has revealed the presence of 2674 CYP128 P450s in 2646 mycobacterial species belonging to six different categories. Contrast features were observed in the CYP128 gene distribution, subfamily patterns, and characteristics of the secondary metabolite biosynthetic gene cluster (BGCs) between M. tuberculosis complex (MTBC) and other mycobacterial category species. In all MTBC species (except one) CYP128 P450s belong to subfamily A, whereas subfamily B is predominant in another four mycobacterial category species. Of CYP128 P450s, 78% was a part of BGCs with CYP124A1, or together with CYP124A1 and CYP121A1. The CYP128 family ranked fifth in the conservation ranking. Unique amino acid patterns are present at the EXXR and CXG motifs. Molecular dynamic simulation studies indicate that the CYP128A1 bind to MK9 with the highest affinity compared to the azole drugs analyzed. This study provides comprehensive comparative analysis and structural insights of CYP128A1 in M. tuberculosis.

The geno-spatio analysis of Mycobacterium tuberculosis complex in hot and cold spots of Guangxi, China

BACKGROUND

At present, there are few studies on polymorphism of Mycobacterium tuberculosis (Mtb) gene and how it affects the TB epidemic. This study aimed to document the differences of polymorphisms between tuberculosis hot and cold spot areas of Guangxi Zhuang Autonomous Region, China.

METHODS

The cold and hot spot areas, each with 3 counties, had been pre-identified by TB incidence for 5 years from the surveillance database. Whole genome sequencing analysis was performed on all sputum Mtb isolates from the detected cases during January and June 2018. Single nucleotide polymorphism (SNP) of each isolate compared to the H37Rv strain were called and used for lineage and sub-lineage identification. Pairwise SNP differences between every pair of isolates were computed. Analyses of Molecular Variance (AMOVA) across counties of the same hot or cold spot area and between the two areas were performed.

RESULTS

As a whole, 59.8% (57.7% sub-lineage 2.2 and 2.1% sub-lineage 2.1) and 39.8% (17.8% sub-lineage 4.4, 6.5% sub-lineage 4.2 and 15.5% sub-lineage 4.5) of the Mtb strains were Lineage 2 and Lineage 4 respectively. The percentages of sub-lineage 2.2 (Beijing family strains) are significantly higher in hot spots. Through the MDS dimension reduction, the genomic population structure in the three hot spot counties is significantly different from those three cold spot counties (T-test p = 0.05). The median of SNPs distances among Mtb isolates in cold spots was greater than that in hot spots (897 vs 746, Rank-sum test p < 0.001). Three genomic clusters, each with genomic distance ≤12 SNPs, were identified with 2, 3 and 4 consanguineous strains. Two clusters were from hot spots and one was from cold spots.

CONCLUSION

Narrower genotype diversity in the hot area may indicate higher transmissibility of the Mtb strains in the area compared to those in the cold spot area.

A systematic evaluation of Mycobacterium tuberculosis Genome-Scale Metabolic Networks

Metabolism underpins the pathogenic strategy of the causative agent of TB, Mycobacterium tuberculosis (Mtb), and therefore metabolic pathways have recently re-emerged as attractive drug targets. A powerful approach to study Mtb metabolism as a whole, rather than just individual enzymatic components, is to use a systems biology framework, such as a Genome-Scale Metabolic Network (GSMN) that allows the dynamic interactions of all the components of metabolism to be interrogated together. Several GSMNs networks have been constructed for Mtb and used to study the complex relationship between the Mtb genotype and its phenotype. However, the utility of this approach is hampered by the existence of multiple models, each with varying properties and performances. Here we systematically evaluate eight recently published metabolic models of Mtb-H37Rv to facilitate model choice. The best performing models, sMtb2018 and iEK1011, were refined and improved for use in future studies by the TB research community.

The current state of animal models and genomic approaches towards identifying and validating molecular determinants of Mycobacterium tuberculosis infection and tuberculosis disease

Animal models are important in understanding both the pathogenesis of and immunity to tuberculosis (TB). Unfortunately, we are beginning to understand that no animal model perfectly recapitulates the human TB syndrome, which encompasses numerous different stages. Furthermore, Mycobacterium tuberculosis infection is a very heterogeneous event at both the levels of pathogenesis and immunity. This review seeks to establish the current understanding of TB pathogenesis and immunity, as validated in the animal models of TB in active use today. We especially focus on the use of modern genomic approaches in these models to determine the mechanism and the role of specific molecular pathways. Animal models have significantly enhanced our understanding of TB. Incorporation of contemporary technologies such as single cell transcriptomics, high-parameter flow cytometric immune profiling, proteomics, proteomic flow cytometry and immunocytometry into the animal models in use will further enhance our understanding of TB and facilitate the development of treatment and vaccination strategies.

Network-based methods for predicting essential genes or proteins: a survey

Genes that are thought to be critical for the survival of organisms or cells are called essential genes. The prediction of essential genes and their products (essential proteins) is of great value in exploring the mechanism of complex diseases, the study of the minimal required genome for living cells and the development of new drug targets. As laboratory methods are often complicated, costly and time-consuming, a great many of computational methods have been proposed to identify essential genes/proteins from the perspective of the network level with the in-depth understanding of network biology and the rapid development of biotechnologies. Through analyzing the topological characteristics of essential genes/proteins in protein-protein interaction networks (PINs), integrating biological information and considering the dynamic features of PINs, network-based methods have been proved to be effective in the identification of essential genes/proteins. In this paper, we survey the advanced methods for network-based prediction of essential genes/proteins and present the challenges and directions for future research.

Creating a Library of Random Transposon Mutants in Leptospira

Generation of a random transposon mutant library is advantageous in Leptospira as site-directed mutagenesis remains a challenge, especially in pathogenic species. This procedure is typically completed by transformation of Leptospira with a Himar1 containing plasmid via conjugation with Escherichia coli as a donor cell. Here we describe the methodology to generate random transposon mutants in the saprophyte Leptospira biflexa via conjugation of plasmid pSW29T-TKS2 harbored in E. coli β2163. Determination of transposon insertion site by semi-random nested PCR will also be described. A similar methodology may be employed to generate Tn mutants of pathogenic Leptospira species.

Mycobacterium tuberculosis LipE Has a Lipase/Esterase Activity and Is Important for Intracellular Growth and In Vivo Infection

Mycobacterium tuberculosis Rv3775 (LipE) was annotated as a putative lipase. However, its lipase activity has never been characterized, and its precise role in tuberculosis (TB) pathogenesis has not been thoroughly studied to date. We overexpressed and purified the recombinant LipE (rLipE) protein and demonstrated that LipE has a lipase/esterase activity. rLipE prefers medium-chain ester substrates, with the maximal activity on hexanoate. Its activity is the highest at 40°C and pH 9. We determined that rLipE hydrolyzes trioctanoate. Using site-directed mutagenesis, we confirmed that the predicted putative activity triad residues Ser97, Gly342, and His363 are essential for the lipase activity of rLipE. The expression of the lipE gene was induced under stressed conditions mimicking M. tuberculosis' intracellular niche. The gene-disrupting mutation of lipE led to significantly reduced bacterial growth inside THP-1 cells and human peripheral blood mononuclear cell-derived macrophages and attenuated M. tuberculosis infection in mice (with ∼8-fold bacterial load reduction in mouse lungs). Our data suggest that LipE functions as a lipase and is important for M. tuberculosis intracellular growth and in vivo infection.

Rv2617c and P36 are virulence factors of pathogenic mycobacteria involved in resistance to oxidative stress

In this study, we characterized the role of Rv2617c in the virulence of Mycobacterium tuberculosis. Rv2617c is a protein of unknown function unique to M. tuberculosis complex (MTC) and Mycobacterium leprae. In vitro, this protein interacts with the virulence factor P36 (also named Erp) and KdpF, a protein linked to nitrosative stress. Here, we showed that knockout of the Rv2617c gene in M. tuberculosis CDC1551 reduced the replication of the pathogen in a mouse model of infection and favored the trafficking of mycobacteria to phagolysosomes. We also demonstrated that Rv2617c and P36 are required for resistance to in vitro hydrogen peroxide treatment in M. tuberculosis and Mycobacterium bovis, respectively. These findings indicate Rv2617c and P36 act in concert to prevent bacterial damage upon oxidative stress.

Mycobacterium tuberculosis CysA2 is a dual sulfurtransferase with activity against thiosulfate and 3-mercaptopyruvate and interacts with mammalian cells

Cyanide is a toxic compound that is converted to the non-toxic thiocyanate by a rhodanese enzyme. Rhodaneses belong to the family of transferases (sulfurtransferases), which are largely studied. The sulfur donor defines the subfamily of these enzymes as thiosulfate:cyanide sulfurtransferases or rhodaneses (TSTs) or 3-mercaptopyruvate sulfurtransfeases (MSTs). In Mycobacterium tuberculosis, the causative agent of tuberculosis, the gene Rv0815c encodes the protein CysA2, a putative uncharacterized thiosulfate:cyanide sulfurtransferase that belongs to the essential sulfur assimilation pathway in the bacillus and is secreted during infection. In this work, we characterized the functional and structural properties of CysA2 and its kinetic parameters. The recombinant CysA2 is a α/β protein with two rhodanese-like domains that maintains the functional motifs and a catalytic cysteine. Sulfurtransferase activity was determined using thiosulfate and 3-mercaptopyruvate as sulfur donors. The assays showed Km values of 2.89 mM and 7.02 mM for thiosulfate and 3-mercaptopyruvate, respectively, indicating the protein has dual activity as TST and MST. Immunological assays revealed that CysA2 interacted with pulmonary cells, and it was capable to activate macrophages and dendritic cells, indicating the stimulation of the immune response, which is important for its use as an antigen for vaccine development and immunodiagnostic.

Peptide Deformylase (def) is essential in Mycobacterium smegmatis, but the essentiality is compensated by inactivation of methionine formylation

Background

Co-translational processes in bacteria are attractive drug targets, but while some processes are essential, others are not. The essentiality of Peptide Deformylase (PDF, def) for vitality of mycobacteria was speculated, but never unequivocally proven.

Results

Here we show by targeted deletion experiments that def can only be deleted from M. smegmatis when an additional copy is present; that prior deletion of tRNA^fMet-Formyl Transferase (FMT, encoded by fmt) renders def completely dispensable; and that re-introduction of fmt into a Δdef mutant is not possible – constituting a definitive proof for the essentiality of def in mycobacteria.

Conclusions

Peptide deformylase is essential in M. smegmatis, but the fact that inactivation of fmt renders the gene completely dispensable, and thus any inhibitor of def useless, casts doubt on the usefulness of PDF as a drug-target in mycobacteria.

Inhibition of CorA-Dependent Magnesium Homeostasis Is Cidal in Mycobacterium tuberculosis

Mechanisms of magnesium homeostasis in Mycobacterium tuberculosis are poorly understood. Here, we describe the characterization of a pyrimidinetrione amide scaffold that disrupts magnesium homeostasis in the pathogen by direct binding to the CorA Mg²⁺/Co²⁺ transporter. Mutations in domains of CorA that are predicted to regulate the pore opening in response to Mg²⁺ ions conferred resistance to this scaffold. The pyrimidinetrione amides were cidal against the pathogen under both actively replicating and nonreplicating conditions in vitro and were efficacious against the organism during macrophage infection. However, the compound lacked efficacy in infected mice, possibly due to limited exposure. Our results indicate that inhibition of Mg²⁺ homeostasis by CorA is an attractive target for tuberculosis drug discovery and encourage identification of improved CorA inhibitors.

A guide to Mycobacterium mutagenesis

The genus Mycobacterium includes several pathogens that cause severe disease in humans, like Mycobacterium tuberculosis (M. tb), the infectious agent causing tuberculosis. Genetic tools to engineer mycobacterial genomes, in a targeted or random fashion, have provided opportunities to investigate M. tb infection and pathogenesis. Furthermore, they have allowed the identification and validation of potential targets for the diagnosis, prevention, and treatment of tuberculosis. This review describes the various methods that are available for the generation of mutants in Mycobacterium species, focusing specifically on tools for altering slow-growing mycobacteria from the M. tb complex. Among others, it incorporates the recent new molecular biological technologies (e.g. ORBIT) to rapidly and/or genome-wide comprehensively obtain targeted mutants in mycobacteria. As such, this review can be used as a guide to select the appropriate genetic tools to generate mycobacterial mutants of interest, which can be used as tools to aid understanding of M. tb infection or to help developing TB intervention strategies.

Direct transmission of within-host Mycobacterium tuberculosis diversity to secondary cases can lead to variable between-host heterogeneity without de novo mutation: A genomic investigation

BACKGROUND

Whole genome sequencing (WGS) has enabled the development of new approaches to track Mycobacterium tuberculosis (Mtb) transmission between tuberculosis (TB) cases but its utility may be challenged by the discovery that Mtb diversifies within hosts. Nevertheless, there is limited data on the presence and degree of within-host evolution.

METHODS

We profiled a well-documented Mtb transmission cluster with three pulmonary TB cases to investigate within-host evolution and describe its impact on recent transmission estimates. We used deep sequencing to track minority allele frequencies (<50·0% abundance) during transmission and standard treatment.

FINDINGS

Pre-treatment (n = 3) and serial samples collected over 2 months of antibiotic treatment (n = 16) from all three cases were analysed. Consistent with the epidemiological data, zero fixed SNP separated all genomes. However, we identified six subclones between the three cases with an allele frequency ranging from 35·0% to 100·0% across sampling intervals. Five subclones were identified within the index case pre-treatment and shared with one secondary case, while only the dominant clone was observed in the other secondary case. By tracking the frequency of these heterogeneous alleles over the two-month therapy, we observed distinct signatures of drift and negative selection, but limited evidence for de novo mutations, even under drug pressure.

INTERPRETATION

We document within-host Mtb diversity in an index case, which led to transmission of minority alleles to a secondary case. Incorporating data on heterogeneous alleles may refine our understanding of Mtb transmission dynamics. However, more evidence is needed on the role of transmission bottleneck on observed heterogeneity between cases.

Genomewide Assessment of Mycobacterium tuberculosis Conditionally Essential Metabolic Pathways

A better understanding of essential cellular functions in pathogenic bacteria is important for the development of more effective antimicrobial agents. We performed a comprehensive identification of essential genes in Mycobacterium tuberculosis, the major causative agent of tuberculosis, using a combination of transposon insertion sequencing (Tn-seq) and comparative genomic analysis. To identify conditionally essential genes by Tn-seq, we used media with different nutrient compositions. Although many conditional gene essentialities were affected by the presence of relevant nutrient sources, we also found that the essentiality of genes in a subset of metabolic pathways was unaffected by metabolite availability. Comparative genomic analysis revealed that not all essential genes identified by Tn-seq were fully conserved within the M. tuberculosis complex, including some existing antitubercular drug target genes. In addition, we utilized an available M. tuberculosis genome-scale metabolic model, iSM810, to predict M. tuberculosis gene essentiality in silico Comparing the sets of essential genes experimentally identified by Tn-seq to those predicted in silico reveals the capabilities and limitations of gene essentiality predictions, highlighting the complexity of M. tuberculosis essential metabolic functions. This study provides a promising platform to study essential cellular functions in M. tuberculosis IMPORTANCE Mycobacterium tuberculosis causes 10 million cases of tuberculosis (TB), resulting in over 1 million deaths each year. TB therapy is challenging because it requires a minimum of 6 months of treatment with multiple drugs. Protracted treatment times and the emergent spread of drug-resistant M. tuberculosis necessitate the identification of novel targets for drug discovery to curb this global health threat. Essential functions, defined as those indispensable for growth and/or survival, are potential targets for new antimicrobial drugs. In this study, we aimed to define gene essentialities of M. tuberculosis on a genomewide scale to comprehensively identify potential targets for drug discovery. We utilized a combination of experimental (functional genomics) and in silico approaches (comparative genomics and flux balance analysis). Our functional genomics approach identified sets of genes whose essentiality was affected by nutrient availability. Comparative genomics revealed that not all essential genes were fully conserved within the M. tuberculosis complex. Comparing sets of essential genes identified by functional genomics to those predicted by flux balance analysis highlighted gaps in current knowledge regarding M. tuberculosis metabolic capabilities. Thus, our study identifies numerous potential antitubercular drug targets and provides a comprehensive picture of the complexity of M. tuberculosis essential cellular functions.

High-Resolution Structure of ClpC1-Rufomycin and Ligand Binding Studies Provide a Framework to Design and Optimize Anti-Tuberculosis Leads

Addressing the urgent need to develop novel drugs against drug-resistant Mycobacterium tuberculosis ( M. tb) strains, ecumicin (ECU) and rufomycin I (RUFI) are being explored as promising new leads targeting cellular proteostasis via the caseinolytic protein ClpC1. Details of the binding topology and chemical mode of (inter)action of these cyclopeptides help drive further development of novel potency-optimized entities as tuberculosis drugs. ClpC1 M. tb protein constructs with mutations driving resistance to ECU and RUFI show reduced binding affinity by surface plasmon resonance (SPR). Despite certain structural similarities, ECU and RUFI resistant mutation sites did not overlap in their SPR binding patterns. SPR competition experiments show ECU prevents RUFI binding, whereas RUFI partially inhibits ECU binding. The X-ray structure of the ClpC1-NTD-RUFI complex reveals distinct differences compared to the previously reported ClpC1-NTD-cyclomarin A structure. Surprisingly, the complex structure revealed that the epoxide moiety of RUFI opened and covalently bound to ClpC1-NTD via the sulfur atom of Met1. Furthermore, RUFI analogues indicate that the epoxy group of RUFI is critical for binding and bactericidal activity. The outcomes demonstrate the significance of ClpC1 as a novel target and the importance of SAR analysis of identified macrocyclic peptides for drug discovery.

Inorganic polyphosphate accumulation suppresses the dormancy response and virulence in Mycobacterium tuberculosis

Stringent response pathways involving inorganic polyphosphate (PolyP) play an essential role in bacterial stress adaptation and virulence. The intracellular levels of PolyP are modulated by the activities of polyphosphate kinase-1 (PPK1), polyphosphate kinase-2 (PPK2), and exopolyphosphatases (PPXs). The genome of Mycobacterium tuberculosis encodes two functional PPXs, and simultaneous deletion of ppx1 and ppx2 results in a defect in biofilm formation. We demonstrate here that these PPXs cumulatively contribute to the ability of M. tuberculosis to survive in nutrient-limiting, low-oxygen growth conditions and also in macrophages. Characterization of single (Δppx2) and double knockout (dkppx) strains of M. tuberculosis indicated that PPX-mediated PolyP degradation is essential for establishing bacterial infection in guinea pigs. RNA-Seq-based transcriptional profiling revealed that relative to the parental strain, the expression levels of DosR regulon-regulated dormancy genes were significantly reduced in the dkppx mutant strain. In concordance, we also provide evidence that PolyP inhibits the autophosphorylation activities associated with DosT and DosS sensor kinases. The results in this study uncover that enzymes involved in PolyP homeostasis play a critical role in M. tuberculosis physiology and virulence and are attractive targets for developing more effective therapeutic interventions.

Structural analysis of molybdopterin synthases from two mycobacterial pathogens

The molybdenum cofactor, composed of molybdopterin and molybdenum, is a necessary compound for the catalytic activity of molybdenum enzymes. Molybdenum cofactor biosynthesis is a conserved multi-step process involving several enzymes. Molybdopterin synthase, a hetero-tetrameric enzyme composed of a pair of MoaE-MoaD subunits, catalyzes the generation of the cis-dithiolene group of molybdopterin in the second step of the process. The cis-dithiolene group can covalently bind molybdenum. Most mycobacterial species possess several genes encoding the full pathway of molybdenum cofactor biosynthesis. In M. smegmatis, the moaD2 and moaE2 genes encode the functional molybdopterin synthase. However, M. tuberculosis has genes encoding several molybdopterin synthase subunit homologs, including moaD1, moaD2, moaE1, moaE2, and moaX, which encodes a MoaD-MoaE fusion protein. Previous studies have shown that moaD2 and moaE2 encode functional molybdopterin synthase. Here, we report the crystal structures of two substrate-free molybdopterin synthases from two different mycobacterial pathogens, M. tuberculosis and M. smegmatis, at 2.1 Å and 2.6 Å resolutions, respectively. The overall structure of both molybdopterin synthases was hetero-tetrameric, consisting of a MoaE2 dimer flanked on either side by single MoaD2 subunits. The carboxyl-terminal domain of MoaD2 inserted into MoaE2, forming the active pocket. A comparison with previously reported molybdopterin synthase structures showed that substrate-binding and catalytic residues were conserved, despite low sequence similarity among these enzymes. The low sequence identity at the MoaE-MoaD heterodimer interface may provide the structural basis to explore mycobacterial inhibitors.

Molecular epidemiology and whole genome sequencing analysis of clinical Mycobacterium bovis from Ghana

BACKGROUND

Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a re-emerging problem in both livestock and humans. The association of some M. bovis strains with hyper-virulence, MDR-TB and disseminated disease makes it imperative to understand the biology of the pathogen.

METHODS

Mycobacterium bovis (15) among 1755 M. tuberculosis complex (MTBC) isolated between 2012 and 2014 were characterized and analyzed for associated patient demography and other risk factors. Five of the M. bovis isolates were whole-genome sequenced and comparatively analyzed against a global collection of published M. bovis genomes.

RESULTS

Mycobacterium bovis was isolated from 3/560(0.5%) females and 12/1195(1.0%) males with pulmonary TB. The average age of M. bovis infected cases was 46.8 years (7-72years). TB patients from the Northern region of Ghana (1.9%;4/212) had a higher rate of infection with M. bovis (OR = 2.7,p = 0.0968) compared to those from the Greater Accra region (0.7%;11/1543). Among TB patients with available HIV status, the odds of isolating M. bovis from HIV patients (2/119) was 3.3 higher relative to non-HIV patients (4/774). Direct contact with livestock or their unpasteurized products was significantly associated with bTB (p<0.0001, OR = 124.4,95% CI = 30.1-508.3). Two (13.3%) of the M. bovis isolates were INH resistant due to the S315T mutation in katG whereas one (6.7%) was RIF resistant with Q432P and I1491S mutations in rpoB. M. bovis from Ghana resolved as mono-phyletic branch among mostly M. bovis from Africa irrespective of the host and were closest to the root of the global M. bovis phylogeny. M. bovis-specific amino acid mutations were detected among MTBC core genes such as mce1A, mmpL1, pks6, phoT, pstB, glgP and Rv2955c. Additional mutations P6T in chaA, G187E in mgtC, T35A in Rv1979c, S387A in narK1, L400F in fas and A563T in eccA1 were restricted to the 5 clinical M. bovis from Ghana.

CONCLUSION

Our data indicate potential zoonotic transmission of bTB in Ghana and hence calls for intensified public education on bTB, especially among risk groups.

Function, essentiality, and expression of cytochrome P450 enzymes and their cognate redox partners in Mycobacterium tuberculosis: are they drug targets?

This review covers the current knowledge of the cytochrome P450 enzymes (CYPs) of the human pathogen Mycobacterium tuberculosis (Mtb) and their endogenous redox partners, focusing on their biological function, expression, regulation, involvement in antibiotic resistance, and suitability for exploitation as antitubercular targets. The Mtb genome encodes twenty  CYPs and nine associated redox partners required for CYP catalytic activity. Transposon insertion mutagenesis studies have established the (conditional) essentiality of several of these enzymes for in vitro growth and host infection. Biochemical characterization of a handful of Mtb CYPs has revealed that they have specific physiological functions in bacterial virulence and persistence in the host. Analysis of the transcriptional response of Mtb CYPs and redox partners to external insults and to first-line antibiotics used to treat tuberculosis showed a diverse expression landscape, suggesting for some enzymes a potential role in drug resistance. Combining the knowledge about the physiological roles and expression profiles indicates that, at least five Mtb CYPs, CYP121A1, CYP125A1, CYP139A1, CYP142A1, and CYP143A1, as well as two ferredoxins, FdxA and FdxC, can be considered promising novel therapeutic targets.