The DNA-binding network of Mycobacterium tuberculosis

Constitutive activation of two-component systems reveals regulatory network interactions in Streptococcus agalactiae

Bacterial two-component systems (TCSs) are signaling modules that control physiology, adaptation, and host interactions. A typical TCS consists of a histidine kinase (HK) that activates a response regulator via phosphorylation in response to environmental signals. Here, we systematically test the effect of inactivating the conserved phosphatase activity of HKs to activate TCS signaling pathways. Transcriptome analyses of 14 HK mutants in Streptococcus agalactiae, the leading cause of neonatal meningitis, validate the conserved HK phosphatase mechanism and its role in the inhibition of TCS activity in vivo. Constitutive TCS activation, independent of environmental signals, enables high-resolution mapping of the regulons for several TCSs (e.g., SaeRS, BceRS, VncRS, DltRS, HK11030, HK02290) and reveals the functional diversity of TCS signaling pathways, ranging from highly specialized to interconnected global regulatory networks. Targeted analysis shows that the SaeRS-regulated PbsP adhesin acts as a signaling molecule to activate CovRS signaling, thereby linking the major regulators of host-pathogen interactions. Furthermore, constitutive BceRS activation reveals drug-independent activity, suggesting a role in cell envelope homeostasis beyond antimicrobial resistance. This study highlights the versatility of constitutive TCS activation, via phosphatase-deficient HKs, to uncover regulatory networks and biological processes.

Predicting bacterial fitness in Mycobacterium tuberculosis with transcriptional regulatory network-informed interpretable machine learning

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis disease, the greatest source of global mortality by a bacterial pathogen. Mtb adapts and responds to diverse stresses such as antibiotics by inducing transcriptional stress-response regulatory programs. Understanding how and when these mycobacterial regulatory programs are activated could enable novel treatment strategies for potentiating the efficacy of new and existing drugs. Here we sought to define and analyze Mtb regulatory programs that modulate bacterial fitness. We assembled a large Mtb RNA expression compendium and applied these to infer a comprehensive Mtb transcriptional regulatory network and compute condition-specific transcription factor activity profiles. We utilized transcriptomic and functional genomics data to train an interpretable machine learning model that can predict Mtb fitness from transcription factor activity profiles. We demonstrated that this transcription factor activity-based model can successfully predict Mtb growth arrest and growth resumption under hypoxia and reaeration using only RNA-seq expression data as a starting point. These integrative network modeling and machine learning analyses thus enable the prediction of mycobacterial fitness under different environmental and genetic contexts. We envision these models can potentially inform the future design of prognostic assays and therapeutic intervention that can cripple Mtb growth and survival to cure tuberculosis disease.

Rhizobium etli CFN42 and Sinorhizobium meliloti 1021 bioinformatic transcriptional regulatory networks from culture and symbiosis

Rhizobium etli CFN42 proteome-transcriptome mixed data of exponential growth and nitrogen-fixing bacteroids, as well as Sinorhizobium meliloti 1021 transcriptome data of growth and nitrogen-fixing bacteroids, were integrated into transcriptional regulatory networks (TRNs). The one-step construction network consisted of a matrix-clustering analysis of matrices of the gene profile and all matrices of the transcription factors (TFs) of their genome. The networks were constructed with the prediction of regulatory network application of the RhizoBindingSites database (http://rhizobindingsites.ccg.unam.mx/). The deduced free-living Rhizobium etli network contained 1,146 genes, including 380 TFs and 12 sigma factors. In addition, the bacteroid R. etli CFN42 network contained 884 genes, where 364 were TFs, and 12 were sigma factors, whereas the deduced free-living Sinorhizobium meliloti 1021 network contained 643 genes, where 259 were TFs and seven were sigma factors, and the bacteroid Sinorhizobium meliloti 1021 network contained 357 genes, where 210 were TFs and six were sigma factors. The similarity of these deduced condition-dependent networks and the biological E. coli and B. subtilis independent condition networks segregates from the random Erdös-Rényi networks. Deduced networks showed a low average clustering coefficient. They were not scale-free, showing a gradually diminishing hierarchy of TFs in contrast to the hierarchy role of the sigma factor rpoD in the E. coli K12 network. For rhizobia networks, partitioning the genome in the chromosome, chromids, and plasmids, where essential genes are distributed, and the symbiotic ability that is mostly coded in plasmids, may alter the structure of these deduced condition-dependent networks. It provides potential TF gen-target relationship data for constructing regulons, which are the basic units of a TRN.

Nucleoid-associated proteins of mycobacteria come with a distinctive flavor

In every bacterium, nucleoid-associated proteins (NAPs) play crucial roles in chromosome organization, replication, repair, gene expression, and other DNA transactions. Their central role in controlling the chromatin dynamics and transcription has been well-appreciated in several well-studied organisms. Here, we review the diversity, distribution, structure, and function of NAPs from the genus Mycobacterium. We highlight the progress made in our understanding of the effects of these proteins on various processes and in responding to environmental stimuli and stress of mycobacteria in their free-living as well as during distinctive intracellular lifestyles. We project them as potential drug targets and discuss future studies to bridge the information gap with NAPs from well-studied systems.

Regulatory role of Mycobacterium tuberculosis MtrA on dormancy/resuscitation revealed by a novel target gene-mining strategy

Introduction

The unique dormancy of Mycobacterium tuberculosis plays a significant role in the major clinical treatment challenge of tuberculosis, such as its long treatment cycle, antibiotic resistance, immune escape, and high latent infection rate.

Methods

To determine the function of MtrA, the only essential response regulator, one strategy was developed to establish its regulatory network according to high-quality genome-wide binding sites.

Results and discussion

The complex modulation mechanisms were implied by the strong bias distribution of MtrA binding sites in the noncoding regions, and 32.7% of the binding sites were located inside the target genes. The functions of 288 potential MtrA target genes predicted according to 294 confirmed binding sites were highly diverse, and DNA replication and damage repair, lipid metabolism, cell wall component biosynthesis, cell wall assembly, and cell division were the predominant pathways. Among the 53 pathways shared between dormancy/resuscitation and persistence, which accounted for 81.5% and 93.0% of the total number of pathways, respectively, MtrA regulatory genes were identified not only in 73.6% of their mutual pathways, but also in 75.4% of the pathways related to dormancy/resuscitation and persistence respectively. These results suggested the pivotal roles of MtrA in regulating dormancy/resuscitation and the apparent relationship between dormancy/resuscitation and persistence. Furthermore, the finding that 32.6% of the MtrA regulons were essential in vivo and/or in vitro for M. tuberculosis provided new insight into its indispensability. The findings mentioned above indicated that MtrA is a novel promising therapeutic target for tuberculosis treatment since the crucial function of MtrA may be a point of weakness for M. tuberculosis.

Structural insights into the regulation mechanism of Mycobacterium tuberculosis MftR

Mycobacterium tuberculosis, the pathogen of the deadly disease tuberculosis, depends on the redox cofactor mycofactocin (MFT) to adapt to and survive under hypoxic conditions. MftR is a TetR family transcription regulator that binds upstream of the MFT gene cluster and controls MFT synthesis. To elucidate the structural basis underlying MftR regulation, we determined the crystal structure of Mycobacterium tuberculosis MftR (TB-MftR). The structure revealed an interconnected hydrogen bond network in the α1-α2-α3 helices of helix-turn-helix (HTH) DNA-binding domain that is essential for nucleic acid interactions. The ligand-binding domain contains a hydrophobic cavity enclosing long-chain fatty acyl-CoAs like the key regulatory ligand oleoyl-CoA. Despite variations in ligand-binding modes, comparative analyses suggest regulatory mechanisms are largely conserved across TetR family acyl-CoA sensors. By elucidating the intricate structural mechanisms governing DNA and ligand binding by TB-MftR, our study enhances understanding of the regulatory roles of this transcription factor under hypoxic conditions, providing insights that could inform future research into Mycobacterium tuberculosis pathogenesis.

Predictive Biophysical Neural Network Modeling of a Compendium of in vivo Transcription Factor DNA Binding Profiles for Escherichia coli

The DNA binding of most Escherichia coli Transcription Factors (TFs) has not been comprehensively mapped, and few have models that can quantitatively predict binding affinity. We report the global mapping of in vivo DNA binding for 139 E. coli TFs using ChIP-Seq. We used these data to train BoltzNet, a novel neural network that predicts TF binding energy from DNA sequence. BoltzNet mirrors a quantitative biophysical model and provides directly interpretable predictions genome-wide at nucleotide resolution. We used BoltzNet to quantitatively design novel binding sites, which we validated with biophysical experiments on purified protein. We have generated models for 125 TFs that provide insight into global features of TF binding, including clustering of sites, the role of accessory bases, the relevance of weak sites, and the background affinity of the genome. Our paper provides new paradigms for studying TF-DNA binding and for the development of biophysically motivated neural networks.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Structural and biophysical characterization of PadR family protein Rv1176c of Mycobacterium tuberculosis H37Rv

Rv1176c of Mycobacterium tuberculosis H37Rv belongs to the PadR-s1 subfamily of the PadR family of protein. Rv1176c forms a stable dimer in solution. Its stability is characterized by a thermal melting transition temperature (Tm) of 39.4 °C. The crystal structure of Rv1176c was determined at a resolution of 2.94 Å, with two monomers in the asymmetric unit. Each monomer has a characteristic N-terminal winged-helix-turn-helix DNA-binding domain. Rv1176c C-terminal is a coiled-coil dimerization domain formed of α-helices α5 to α7. In the Rv1176c dimer, there is domain-swapping of the C-terminal domain in comparison to other PadR homologs. In the dimer, there is a long inter-subunit tunnel in which different ligands can bind. Rv1176c was found to bind to the promoter region of its own gene with high specificity. M. smegmatis MC2 155 genome lacks homolog of Rv1176c. Therefore, it was used as a surrogate to characterize the functional role of Rv1176c. Expression of Rv1176c in M. smegmatis MC2 155 cells imparted enhanced tolerance towards oxidative stress. Rv1176c expressing M. smegmatis MC2 155 cells exhibited enhanced intracellular survival in J774A.1 murine macrophage cells. Overall, our studies demonstrate Rv1176c to be a PadR-s1 subfamily transcription factor that can moderate the effect of oxidative stress.

MLDSR, the transcriptional regulator of the major lipid droplets protein MLDS, is controlled by long-chain fatty acids and contributes to the lipid-accumulating phenotype in oleaginous Rhodococcus strains

Microorganism lipid droplet small regulator (MLDSR) is a transcriptional regulator of the major lipid droplet (LD)-associated protein MLDS in Rhodococcus jostii RHA1 and Rhodococcus opacus PD630. In this study, we investigated the role of MLDSR on lipid metabolism and triacylglycerol (TAG) accumulation in R. jostii RHA1 at physiological and molecular levels. MLDSR gene deletion promoted a significant decrease of TAG accumulation, whereas inhibition of de novo fatty acid biosynthesis by the addition of cerulenin significantly repressed the expression of the mldsr-mlds cluster under nitrogen-limiting conditions. In vitro and in vivo approaches revealed that MLDSR-DNA binding is inhibited by fatty acids and acyl-CoA residues through changes in the oligomeric or conformational state of the protein. RNAseq analysis indicated that MLDSR not only controls the expression of its own gene cluster but also of several genes involved in central, lipid, and redox metabolism, among others. We also identified putative MLDSR-binding sites on the upstream regions of genes coding for lipid catabolic enzymes and validated them by EMSA assays. Overexpression of mldsr gene under nitrogen-rich conditions promoted an increase of TAG accumulation, and further cell lysis with TAG release to the culture medium. Our results suggested that MLDSR is a fatty acid-responsive regulator that plays a dual role in cells by repression or activation of several metabolic genes in R. jostii RHA1. MLDSR seems to play an important role in the fine-tuning regulation of TAG accumulation, LD formation, and cellular lipid homeostasis, contributing to the oleaginous phenotype of R. jostii RHA1 and R. opacus PD630.

Harnessing regulatory networks in Actinobacteria for natural product discovery

Microbes typically live in complex habitats where they need to rapidly adapt to continuously changing growth conditions. To do so, they produce an astonishing array of natural products with diverse structures and functions. Actinobacteria stand out for their prolific production of bioactive molecules, including antibiotics, anticancer agents, antifungals, and immunosuppressants. Attention has been directed especially towards the identification of the compounds they produce and the mining of the large diversity of biosynthetic gene clusters (BGCs) in their genomes. However, the current return on investment in random screening for bioactive compounds is low, while it is hard to predict which of the millions of BGCs should be prioritized. Moreover, many of the BGCs for yet undiscovered natural products are silent or cryptic under laboratory growth conditions. To identify ways to prioritize and activate these BGCs, knowledge regarding the way their expression is controlled is crucial. Intricate regulatory networks control global gene expression in Actinobacteria, governed by a staggering number of up to 1000 transcription factors per strain. This review highlights recent advances in experimental and computational methods for characterizing and predicting transcription factor binding sites and their applications to guide natural product discovery. We propose that regulation-guided genome mining approaches will open new avenues toward eliciting the expression of BGCs, as well as prioritizing subsets of BGCs for expression using synthetic biology approaches.

ONE-SENTENCE SUMMARY

This review provides insights into advances in experimental and computational methods aimed at predicting transcription factor binding sites and their applications to guide natural product discovery.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Transcription factor mce3R modulates antibiotics and disease persistence in Mycobacteriumtuberculosis

Transcription factors (TFs) of Mycobacterium tuberculosis (Mtb), an etiological agent of tuberculosis, regulate a network of pathways that help prolong the survival of Mtb inside the host. In this study, we have characterized a transcription repressor gene (mce3R) from the TetR family, that encodes for Mce3R protein in Mtb. We demonstrated that the mce3R gene is dispensable for the growth of Mtb on cholesterol. Gene expression analysis suggests that the transcription of genes belonging to the mce3R regulon is independent of the carbon source. We found that, in comparison to the wild type, the mce3R deleted strain (Δmce3R) generated more intracellular ROS and demonstrated reduced susceptibility to oxidative stress. Total lipid analysis suggests that mce3R regulon encoded proteins modulate the biosynthesis of cell wall lipids in Mtb. Interestingly, the absence of Mce3R increased the frequency of generation of antibiotic persisters in Mtb and imparted in-vivo growth advantage phenotype in guinea pigs. In conclusion, genes belonging to the mce3R regulon modulate the frequency of generation of persisters in Mtb. Hence, targeting mce3R regulon encoded proteins could potentiate the current regimen by eliminating persisters during Mtb infection.

MtrA modulates Mycobacterium tuberculosis cell division in host microenvironments to mediate intrinsic resistance and drug tolerance

The success of Mycobacterium tuberculosis (Mtb) is largely attributed to its ability to physiologically adapt and withstand diverse localized stresses within host microenvironments. Here, we present a data-driven model (EGRIN 2.0) that captures the dynamic interplay of environmental cues and genome-encoded regulatory programs in Mtb. Analysis of EGRIN 2.0 shows how modulation of the MtrAB two-component signaling system tunes Mtb growth in response to related host microenvironmental cues. Disruption of MtrAB by tunable CRISPR interference confirms that the signaling system regulates multiple peptidoglycan hydrolases, among other targets, that are important for cell division. Further, MtrA decreases the effectiveness of antibiotics by mechanisms of both intrinsic resistance and drug tolerance. Together, the model-enabled dissection of complex MtrA regulation highlights its importance as a drug target and illustrates how EGRIN 2.0 facilitates discovery and mechanistic characterization of Mtb adaptation to specific host microenvironments within the host.

Biosensor-integrated transposon mutagenesis reveals rv0158 as a coordinator of redox homeostasis in Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is evolutionarily equipped to resist exogenous reactive oxygen species (ROS) but shows vulnerability to an increase in endogenous ROS (eROS). Since eROS is an unavoidable consequence of aerobic metabolism, understanding how Mtb manages eROS levels is essential yet needs to be characterized. By combining the Mrx1-roGFP2 redox biosensor with transposon mutagenesis, we identified 368 genes (redoxosome) responsible for maintaining homeostatic levels of eROS in Mtb. Integrating redoxosome with a global network of transcriptional regulators revealed a hypothetical protein (Rv0158) as a critical node managing eROS in Mtb. Disruption of rv0158 (rv0158 KO) impaired growth, redox balance, respiration, and metabolism of Mtb on glucose but not on fatty acids. Importantly, rv0158 KO exhibited enhanced growth on propionate, and the Rv0158 protein directly binds to methylmalonyl-CoA, a key intermediate in propionate catabolism. Metabolite profiling, ChIP-Seq, and gene-expression analyses indicate that Rv0158 manages metabolic neutralization of propionate toxicity by regulating the methylcitrate cycle. Disruption of rv0158 enhanced the sensitivity of Mtb to oxidative stress, nitric oxide, and anti-TB drugs. Lastly, rv0158 KO showed poor survival in macrophages and persistence defect in mice. Our results suggest that Rv0158 is a metabolic integrator for carbon metabolism and redox balance in Mtb.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Mycobacterium tuberculosis Rv0494 Protein Contributes to Mycobacterial Persistence

Purpose

Fatty acid metabolism plays an important role in the survival and pathogenesis of Mycobacterium tuberculosis. During dormancy, lipids are considered to be the main source of energy. A previous study found that Rv0494 is a starvation-inducible, lipid-responsive transcriptional regulator. However, the role of Rv0494 in bacterial persister survival has not been studied.

Methods

We constructed a Rv0494 deletion mutant strain of Mycobacterium tuberculosis H37Rv and evaluated the susceptibility of the mutant strain to antibiotics using a persistence test.

Results

We found that mutations in Rv0494 lead to survival defects of persisters, which reflected in increased sensitivity to isoniazid.

Conclusion

We conclude that Rv0494 is important for persister survival and may serve as a good target for developing new antibiotics that kill persister bacteria for improved treatment of persistent bacterial infections.

ResR/McdR-regulated protein translation machinery contributes to drug resilience in Mycobacterium tuberculosis

Survival response of the human tuberculosis pathogen, Mycobacterium tuberculosis (Mtb) to a diverse environmental cues is governed through its versatile transcription regulatory mechanisms with the help of a large pool of transcription regulators (TRs). Rv1830 is one such conserved TR, which remains uncharacterized in Mtb. It was named as McdR based on an effect on cell division upon its overexpression in Mycobacterium smegmatis. Recently, it has been implicated in antibiotic resilience in Mtb and reannotated as ResR. While Rv1830 affects cell division by modulating the expression of M. smegmatis whiB2, the underlying cause of its essentiality and regulation of drug resilience in Mtb is yet to be deciphered. Here we show that ResR/McdR, encoded by ERDMAN_2020 in virulent Mtb Erdman, is pivotal for bacterial proliferation and crucial metabolic activities. Importantly, ResR/McdR directly regulates ribosomal gene expression and protein synthesis, requiring distinct disordered N-terminal sequence. Compared to control, bacteria depleted with resR/mcdR exhibit delayed recovery post-antibiotic treatment. A similar effect upon knockdown of rplN operon genes further implicates ResR/McdR-regulated protein translation machinery in attributing drug resilience in Mtb. Overall, findings from this study suggest that chemical inhibitors of ResR/McdR may be proven effective as adjunctive therapy for shortening the duration of TB treatment.

DNA structural properties of DNA binding sites for 21 transcription factors in the mycobacterial genome

Mycobacterium tuberculosis, the causative agent of tuberculosis, has evolved over time into a multidrug resistance strain that poses a serious global pandemic health threat. The ability to survive and remain dormant within the host macrophage relies on multiple transcription factors contributing to virulence. To date, very limited structural insights from crystallographic and NMR studies are available for TFs and TF-DNA binding events. Understanding the role of DNA structure in TF binding is critical to deciphering MTB pathogenicity and has yet to be resolved at the genome scale. In this work, we analyzed the compositional and conformational preference of 21 mycobacterial TFs, evident at their DNA binding sites, in local and global scales. Results suggest that most TFs prefer binding to genomic regions characterized by unique DNA structural signatures, namely, high electrostatic potential, narrow minor grooves, high propeller twist, helical twist, intrinsic curvature, and DNA rigidity compared to the flanking sequences. Additionally, preference for specific trinucleotide motifs, with clear periodic signals of tetranucleotide motifs, are observed in the vicinity of the TF-DNA interactions. Altogether, our study reports nuanced DNA shape and structural preferences of 21 TFs.

Transcription regulation by CarD in mycobacteria is guided by basal promoter kinetics

Bacterial pathogens like Mycobacterium tuberculosis (Mtb) employ transcription factors to adapt their physiology to the diverse environments within their host. CarD is a conserved bacterial transcription factor that is essential for viability in Mtb. Unlike classical transcription factors that recognize promoters by binding to specific DNA sequence motifs, CarD binds directly to the RNA polymerase to stabilize the open complex intermediate (RPo) during transcription initiation. We previously showed using RNA-sequencing that CarD is capable of both activating and repressing transcription in vivo. However, it is unknown how CarD achieves promoter-specific regulatory outcomes in Mtb despite binding indiscriminate of DNA sequence. We propose a model where CarD's regulatory outcome depends on the promoter's basal RPo stability and test this model using in vitro transcription from a panel of promoters with varying levels of RPo stability. We show that CarD directly activates full-length transcript production from the Mtb ribosomal RNA promoter rrnAP3 (AP3) and that the degree of transcription activation by CarD is negatively correlated with RPo stability. Using targeted mutations in the extended -10 and discriminator region of AP3, we show that CarD directly represses transcription from promoters that form relatively stable RPo. DNA supercoiling also influenced RPo stability and affected the direction of CarD regulation, indicating that the outcome of CarD activity can be regulated by factors beyond promoter sequence. Our results provide experimental evidence for how RNA polymerase-binding transcription factors like CarD can exert specific regulatory outcomes based on the kinetic properties of a promoter.

Transcription regulation by CarD in mycobacteria is guided by basal promoter kinetics

Bacterial pathogens like Mycobacterium tuberculosis ( Mtb ) employ transcription factors to adapt their physiology to the diverse environments within their host. CarD is a conserved bacterial transcription factor that is essential for viability in Mtb . Unlike classical transcription factors that recognize promoters by binding to specific DNA sequence motifs, CarD binds directly to the RNA polymerase (RNAP) to stabilize the open complex intermediate (RP _o ) during transcription initiation. We previously showed using RNA-sequencing that CarD is capable of both activating and repressing transcription in vivo . However, it is unknown how CarD achieves promoter specific regulatory outcomes in Mtb despite binding indiscriminate of DNA sequence. We propose a model where CarD's regulatory outcome depends on the promoter's basal RP _o stability and test this model using in vitro transcription from a panel of promoters with varying levels of RP _o stability. We show that CarD directly activates full-length transcript production from the Mtb ribosomal RNA promoter rrnA P3 (AP3) and that the degree of transcription activation by CarD is negatively correlated with RP _o stability. Using targeted mutations in the extended -10 and discriminator region of AP3, we show that CarD directly represses transcription from promoters that form relatively stable RP _o . DNA supercoiling also influenced RP _o stability and affected the direction of CarD regulation, indicating that the outcome of CarD activity can be regulated by factors beyond promoter sequence. Our results provide experimental evidence for how RNAP-binding transcription factors like CarD can exert specific regulatory outcomes based on the kinetic properties of a promoter.

Mycobacterium tuberculosis transcriptional regulator Rv1019 is upregulated in hypoxia, and negatively regulates Rv3230c-Rv3229c operon encoding enzymes in the oleic acid biosynthetic pathway

The main obstacle in eradicating tuberculosis is the ability of Mycobacterium tuberculosis to remain dormant in the host, and then to get reactivated even years later under immunocompromised conditions. Transcriptional regulation in intracellular pathogens plays an important role in their adapting to the challenging environment inside the host cells. Previously, we demonstrated that Rv1019, a putative transcriptional regulator of M. tuberculosis H37Rv, is an autorepressor. We showed that Rv1019 is cotranscribed with Rv1020 (mfd) and Rv1021 (mazG) which encode DNA repair proteins and negatively regulates the expression of these genes. In the present study, we show that Rv1019 regulates the expression of the genes Rv3230c and Rv3229c (desA3) also which form a two-gene operon in M. tuberculosis. Overexpression of Rv1019 in M. tuberculosis significantly downregulated the expression of these genes. Employing Wayne's hypoxia-induced dormancy model of M. tuberculosis, we show that Rv1019 is upregulated three-fold under hypoxia. Finally, by reporter assay, using Mycobacterium smegmatis as a model, we validate that Rv1019 is recruited to the promoter of Rv3230c-Rv3229c during hypoxia, and negatively regulates this operon which is involved in the biosynthesis of oleic acid.

A systematic computational analysis of Mycobacterium tuberculosis H37Rv and human CD34+ genomic expression reveals crucial molecular entities involved in infection progression

The co-evolution of Mycobacterium tuberculosis H37Rv along with its host systems enables the pathogenic bacterium to emerge as a multi-drug resistant form. This creates challenges for a more efficacious treatment strategy that can mitigate the infection. Working towards the same, our study followed a mathematical and statistical approach proposing that mycobacterial transcription factors regulating virulence and adaptation, host cell cytoplasmic component metabolism, oxidoreductase activity and respiratory ETC would be targets for antibiotics against Mycobacterium tuberculosis. Simultaneously, extending the statistical study on Mycobacterium-infected human cord blood CD34+ cells revealed that the human CD34+ genes, S100A8 and FGR (tyrosine-protein kinase, Src2), might be affected in the infection pathogenesis by Mycobacterium. Further, the deduced Mycobacterium-human gene interaction network proposed that mycobacterial coregulators Rv0452 (MarR family regulator) and Rv3862c (WhiB6) triggered genes controlling bacterial metabolism, which influences human immunological pathways involving TLR2 and CXCL8/MAPK8.Communicated by Ramaswamy H. Sarma.

Exploring the mono-/bistability range of positively autoregulated signaling systems in the presence of competing transcription factor binding sites

Binding of transcription factor (TF) proteins to regulatory DNA sites is key to accurate control of gene expression in response to environmental stimuli. Theoretical modeling of transcription regulation is often focused on a limited set of genes of interest, while binding of the TF to other genomic sites is seldom considered. The total number of TF binding sites (TFBSs) affects the availability of TF protein molecules and sequestration of a TF by TFBSs can promote bistability. For many signaling systems where a graded response is desirable for continuous control over the input range, biochemical parameters of the regulatory proteins need be tuned to avoid bistability. Here we analyze the mono-/bistable parameter range for positively autoregulated two-component systems (TCSs) in the presence of different numbers of competing TFBSs. TCS signaling, one of the major bacterial signaling strategies, couples signal perception with output responses via protein phosphorylation. For bistability, competition for TF proteins by TFBSs lowers the requirement for high fold change of the autoregulated transcription but demands high phosphorylation activities of TCS proteins. We show that bistability can be avoided with a low phosphorylation capacity of TCSs, a high TF affinity for the autoregulated promoter or a low fold change in signaling protein levels upon induction. These may represent general design rules for TCSs to ensure uniform graded responses. Examining the mono-/bistability parameter range allows qualitative prediction of steady-state responses, which are experimentally validated in the E. coli CusRS system.

Bacterial Transcription Factors Bind to Coding Regions and Regulate Internal Cryptic Promoters

Transcription factors (TFs) regulate transcription by binding to the specific sequences at the promoter region. However, the mechanisms and functions of TFs binding within the coding sequences (CDS) remain largely elusive in prokaryotes. To this end, we collected 409 data sets for bacterial TFs, including 104 chromatin immunoprecipitation sequencing (ChIP-seq) assays and 305 data sets from the systematic evolution of ligands by exponential enrichment (SELEX) in seven model bacteria. Interestingly, these TFs displayed the same binding capabilities for both coding and intergenic regions. Subsequent biochemical and genetic experiments demonstrated that several TFs bound to the coding regions and regulated the transcription of the binding or adjacent genes. Strand-specific RNA sequencing revealed that these CDS-binding TFs regulated the activity of the cryptic promoters, resulting in the altered transcription of the corresponding antisense RNA. TF RhpR hindered the transcriptional elongation of a subgenic transcript within a CDS. A ChIP-seq and Ribo-seq coanalysis revealed that RhpR influenced the translational efficiency of binding genes. Taken together, the present study reveals three regulatory mechanisms of CDS-bound TFs within individual genes, operons, and antisense RNAs, which demonstrate the variability of the regulatory mechanisms of TFs and expand upon the complexity of bacterial transcriptomes.

Iron deprivation enhances transcriptional responses to in vitro growth arrest of Mycobacterium tuberculosis

The establishment of Mycobacterium tuberculosis (Mtb) long-term infection in vivo depends on several factors, one of which is the availability of key nutrients such as iron (Fe). The relation between Fe deprivation inside and outside the granuloma, and the capacity of Mtb to accumulate lipids and persist in the absence of growth is not well understood. In this context, current knowledge of how Mtb modifies its lipid composition in response to growth arrest, depending on iron availability, is scarce. To shed light on these matters, in this work we compare genome-wide transcriptomic and lipidomic profiles of Mtb at exponential and stationary growth phases using cultures with glycerol as a carbon source, in the presence or absence of iron. As a result, we found that transcriptomic responses to growth arrest, considered as the transition from exponential to stationary phase, are iron dependent for as many as 714 genes (iron-growth interaction contrast, FDR <0.05), and that, in a majority of these genes, iron deprivation enhances the magnitude of the transcriptional responses to growth arrest in either direction. On the one hand, genes whose upregulation upon growth arrest is enhanced by iron deprivation were enriched in functional terms related to homeostasis of ion metals, and responses to several stressful cues considered cardinal features of the intracellular environment. On the other hand, genes showing negative responses to growth arrest that are stronger in iron-poor medium were enriched in energy production processes (TCA cycle, NADH dehydrogenation and cellular respiration), and key controllers of ribosomal activity shut-down, such as the T/A system mazE6/F6. Despite of these findings, a main component of the cell envelope, lipid phthiocerol dimycocerosate (PDIM), was not detected in the stationary phase regardless of iron availability, suggesting that lipid changes during Mtb adaptation to non-dividing phenotypes appear to be iron-independent. Taken together, our results indicate that environmental iron levels act as a key modulator of the intensity of the transcriptional adaptations that take place in the bacterium upon its transition between dividing and dormant-like phenotypes in vitro.

M. tuberculosis AlkX Encoded by rv3249c Regulates a Conserved Alkane Hydroxylase System That Is Important for Replication in Macrophages and Biofilm Formation

Mycobacterium tuberculosis is a highly specialized human pathogen. The success of M. tuberculosis is due to its ability to replicate within host macrophages, resist host immune responses, and ultimately enter a persistent state during a latent tuberculosis infection. Understanding how M. tuberculosis adapts to and replicates in the intracellular environment of the host is crucial for the development of novel, targeted therapeutics. We report the characterization of an M. tuberculosis mutant lacking Rv3249c, a TetR transcriptional regulator. We show that Rv3249c directly represses the adjacent alkB-rubA-rubB operon encoding an alkane hydroxylase/rubredoxin system. For consistency with related systems, we have named the rv3249c gene alkX. The alkX mutant survived better than wild-type M. tuberculosis inside macrophages. This could be phenocopied by overexpression of the alkB-rubA-rubB locus. We hypothesized that the improved intracellular survival phenotype is a result of increased fitness of the mutant; however, we found that the alkX mutant had a defect when grown on some host-associated carbon sources in vitro. We also found that the alkX mutant had a defect in biofilm formation, also linked to the overexpression of the alkB-rubAB genes. Combined, these results define the primary role of AlkX as a transcriptional repressor of the alkB-rubAB operon and suggest the operon contributes to intracellular survival of the pathogen. IMPORTANCE Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), is the leading cause of death worldwide due to a single infectious agent. It is important to understand how M. tuberculosis adapts to and replicates in the intracellular environment of the host. In this study, we characterized the TetR transcriptional regulator Rv3249c and show that it regulates a highly conserved alkane hydroxylase/rubredoxin system. Our data demonstrate that the AlkBRubAB system contributes to the success of the bacterium in host macrophages.

PrrA modulates Mycobacterium tuberculosis response to multiple environmental cues and is critically regulated by serine/threonine protein kinases

The ability of Mycobacterium tuberculosis (Mtb) to adapt to its surrounding environment is critical for the bacterium to successfully colonize its host. Transcriptional changes are a vital mechanism by which Mtb responds to key environmental signals experienced, such as pH, chloride (Cl^-), nitric oxide (NO), and hypoxia. However, much remains unknown regarding how Mtb coordinates its response to the disparate signals seen during infection. Utilizing a transcription factor (TF) overexpression plasmid library in combination with a pH/Cl^--responsive luciferase reporter, we identified the essential TF, PrrA, part of the PrrAB two-component system, as a TF involved in modulation of Mtb response to pH and Cl^-. Further studies revealed that PrrA also affected Mtb response to NO and hypoxia, with prrA overexpression dampening induction of NO and hypoxia-responsive genes. PrrA is phosphorylated not just by its cognate sensor histidine kinase PrrB, but also by serine/threonine protein kinases (STPKs) at a second distinct site. Strikingly, a STPK-phosphoablative PrrA variant was significantly dampened in its response to NO versus wild type Mtb, disrupted in its ability to adaptively enter a non-replicative state upon extended NO exposure, and attenuated for in vivo colonization. Together, our results reveal PrrA as an important regulator of Mtb response to multiple environmental signals, and uncover a critical role of STPK regulation of PrrA in its function.

Vital to successful host colonization by Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis, is the bacterium’s ability to respond and adapt to changes in its local environment during infection. Here, we discover that the essential transcription factor PrrA, part of the PrrAB two-component system (TCS), modulates Mtb response to four important environmental cues encountered within the host: pH, chloride, nitric oxide, and hypoxia. PrrA acts as a rheostat, adjusting the amplitude of Mtb gene expression changes upon bacterial exposure to each of the four environmental signals. Further, we reveal a critical impact of serine/threonine protein kinases (STPKs) on PrrA function, with prevention of STPK phosphorylation of PrrA disrupting adaptive response of Mtb to growth-inhibiting cues and attenuating the bacterium’s ability to colonize its host. Our work uncovers PrrA as a regulator with broad impact across environmental signals, and highlights how two regulatory systems, TCSs and STPKs, critically interact in coordinating Mtb response to environmental cues.

Insight Into Novel Anti-tuberculosis Vaccines by Using Immunoinformatics Approaches

Tuberculosis (TB), an infectious disease, has been a leading cause of morbidity and mortality for decades. The causative agent of TB is the Mycobacterium tuberculosis (Mtb) which can infects various parts of the body, mainly the lungs in pulmonary TB cases. Mycobacterium bovis Bacillus Calmette–Guerin (BCG) is the only approved vaccine for TB, but its efficiency to combat pulmonary TB is limited. Multidrug-resistant (MDR) TB and extensive drug-resistant (XDR) TB requires the evolution of more potent vaccines. Therefore, this research aims to generate a universal TB subunit vaccine using advanced immunoinformatics techniques. In generating a novel multiepitope subunit vaccine, we selected the conserved and experimentally confirmed antigens Rv0058, Rv0101, and Rv3343. After a rigorous evaluation, the top candidates from predicted Helper T-lymphocytes (HTL), Cytotoxic T-lymphocytes (CTL), and B-cell epitopes were considered potential vaccine candidates. Immunogenicity was enhanced by the addition of an adjuvant to the ultimate construct of the vaccine. B-cell epitopes predictions guaranteed the eventual induction of a humoral response. Thereafter, dynamics simulations and molecular docking validated the vaccine-receptor complex’s stability and high affinity for the immune receptor TLR-3. Also, immune simulations revealed the significantly elevated levels of immunoglobulins such as IgM, cytokines such as interleukin-2, helper T (Th) cells, and cytotoxic T-cell populations. These results agreed with the actual inflammatory response and showed rapid antigen clearance after manifold exposure. Finally, the E. coli K12 strain was confirmed via in-silico cloning for quality expression. Nevertheless, in vivo experiments should be performed to validate the safety of the proposed vaccine and its inherent ability to prevent TB infection.

Investigating a putative transcriptional regulatory protein encoded by Rv1719 gene of Mycobacterium tuberculosis

Mycobacterium tuberculosis, the causative agent of tuberculosis, demonstrates immense plasticity with which it adapts to a highly dynamic and hostile host environment. This is facilitated by a web of signalling pathways constantly modulated by a multitude of proteins that regulate the flow of genetic information inside the pathogen. Transcription factors (TFs) belongs to one such family of proteins that modulate the signalling by regulating the abundance of proteins at the transcript level. In the current study, we have characterized the putative transcriptional regulatory protein encoded by the Rv1719 gene of Mycobacterium tuberculosis. This TF belongs to the IclR family of proteins with orthologs found in both bacterial and archaeal species. We cloned the Rv1719 gene into the pET28a expression vector and performed heterologous expression of the recombinant protein with E coli as the host. Further, optimization of the purification protocol by affinity chromatography and characterization of proteins for their functional viability has been demonstrated using various biochemical and/or biophysical approaches. Scale-up of purification yielded approximately 30 mg of ~ 28 kDa protein per litre of culture. In-silico protein domain analysis of Rv1719 protein predicted the presence of the helix-turn-helix (HTH) domain suggesting its ability to bind DNA sequence and modulate transcription; a hallmark of a transcriptional regulatory protein. Further, by performing electrophoretic mobility shift assay (EMSA) we demonstrated that the protein binds to a specific DNA fragment harboring the probable binding site of one of the predicted promoters.

Mycobacterium tuberculosis Transcription Factor EmbR Regulates the Expression of Key Virulence Factors That Aid in Ex Vivo and In Vivo Survival

Mycobacterium tuberculosis encodes ~200 transcription factors that modulate gene expression under different microenvironments in the host. Even though high-throughput chromatin immunoprecipitation sequencing and transcriptome sequencing (RNA-seq) studies have identified the regulatory network for ~80% of transcription factors, many transcription factors remain uncharacterized. EmbR is one such transcription factor whose in vivo regulon and biological function are yet to be elucidated. Previous in vitro studies suggested that phosphorylation of EmbR by PknH upregulates the embCAB operon. Using a gene replacement mutant of embR, we investigated its role in modulating cellular morphology, antibiotic resistance, and survival in the host. Contrary to the prevailing hypothesis, under normal growth conditions, EmbR is neither phosphorylated nor impacted by ethambutol resistance through the regulation of the embCAB operon. The embR deletion mutant displayed attenuated M. tuberculosis survival in vivo. RNA-seq analysis suggested that EmbR regulates operons involved in the secretion pathway, lipid metabolism, virulence, and hypoxia, including well-known hypoxia-inducible genes devS and hspX. Lipidome analysis revealed that EmbR modulates levels of all lysophospholipids, several phospholipids, and M. tuberculosis-specific lipids, which is more pronounced under hypoxic conditions. We found that the EmbR mutant is hypersusceptible to hypoxic stress, and RNA sequencing performed under hypoxic conditions indicated that EmbR majorly regulates genes involved in response to acidic pH, hypoxia, and fatty acid metabolism. We observed condition-specific phosphorylation of EmbR, which contributes to EmbR-mediated transcription of several essential genes, ensuring enhanced survival. Collectively, the study establishes EmbR as a key modulator of hypoxic response that facilitates mycobacterial survival in the host.

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.

Machine Learning of All Mycobacterium tuberculosis H37Rv RNA-seq Data Reveals a Structured Interplay between Metabolism, Stress Response, and Infection

Mycobacterium tuberculosis is one of the most consequential human bacterial pathogens, posing a serious challenge to 21st century medicine. A key feature of its pathogenicity is its ability to adapt its transcriptional response to environmental stresses through its transcriptional regulatory network (TRN). While many studies have sought to characterize specific portions of the M. tuberculosis TRN, and some studies have performed system-level analysis, few have been able to provide a network-based model of the TRN that also provides the relative shifts in transcriptional regulator activity triggered by changing environments. Here, we compiled a compendium of nearly 650 publicly available, high quality M. tuberculosis RNA-sequencing data sets and applied an unsupervised machine learning method to obtain a quantitative, top-down TRN. It consists of 80 independently modulated gene sets known as “iModulons,” 41 of which correspond to known regulons. These iModulons explain 61% of the variance in the organism’s transcriptional response. We show that iModulons (i) reveal the function of poorly characterized regulons, (ii) describe the transcriptional shifts that occur during environmental changes such as shifting carbon sources, oxidative stress, and infection events, and (iii) identify intrinsic clusters of regulons that link several important metabolic systems, including lipid, cholesterol, and sulfur metabolism. This transcriptome-wide analysis of the M. tuberculosis TRN informs future research on effective ways to study and manipulate its transcriptional regulation and presents a knowledge-enhanced database of all published high-quality RNA-seq data for this organism to date.

IMPORTANCEMycobacterium tuberculosis H37Rv is one of the world's most impactful pathogens, and a large part of the success of the organism relies on the differential expression of its genes to adapt to its environment. The expression of the organism's genes is driven primarily by its transcriptional regulatory network, and most research on the TRN focuses on identifying and quantifying clusters of coregulated genes known as regulons. While previous studies have relied on molecular measurements, in the manuscript we utilized an alternative technique that performs machine learning to a large data set of transcriptomic data. This approach is less reliant on hypotheses about the role of specific regulatory systems and allows for the discovery of new biological findings for already collected data. A better understanding of the structure of the M. tuberculosis TRN will have important implications in the design of improved therapeutic approaches.

Multiple genetic paths including massive gene amplification allow Mycobacterium tuberculosis to overcome loss of ESX-3 secretion system substrates

The Mycobacterium tuberculosis (Mtb) ESX-3 type VII secretion system plays a critical role in iron acquisition. Infection of mice with highly attenuated Mtb deletion mutants lacking esxG or esxH, genes encoding key ESX-3 substrates, unexpectedly yielded suppressor mutants with restored capacity to grow in vivo and in vitro in the absence of iron supplementation. Whole-genome sequencing identified two mechanisms of suppression, the disruption of a transcriptional repressor that regulates expression of an ESX-3 paralogous region encoding EsxR and EsxS, and a massive 38- to 60-fold gene amplification of this same region. These data are significant because they reveal a previously unrecognized iron acquisition regulon and inform mechanisms of Mtb chromosome evolution.

Mycobacterium tuberculosis (Mtb) possesses five type VII secretion systems (T7SS), virulence determinants that include the secretion apparatus and associated secretion substrates. Mtb strains deleted for the genes encoding substrates of the ESX-3 T7SS, esxG or esxH, require iron supplementation for in vitro growth and are highly attenuated in vivo. In a subset of infected mice, suppressor mutants of esxG or esxH deletions were isolated, which enabled growth to high titers or restored virulence. Suppression was conferred by mechanisms that cause overexpression of an ESX-3 paralogous region that lacks genes for the secretion apparatus but encodes EsxR and EsxS, apparent ESX-3 orphan substrates that functionally compensate for the lack of EsxG or EsxH. The mechanisms include the disruption of a transcriptional repressor and a massive 38- to 60-fold gene amplification. These data identify an iron acquisition regulon, provide insight into T7SS, and reveal a mechanism of Mtb chromosome evolution involving “accordion-type” amplification.

MadR mediates acyl CoA-dependent regulation of mycolic acid desaturation in mycobacteria

Mycobacterium tuberculosis has a lipid-rich cell envelope that is remodeled throughout infection to enable adaptation within the host. Few transcriptional regulators have been characterized that coordinate synthesis of mycolic acids, the major cell wall lipids of mycobacteria. Here, we show that the mycolic acid desaturase regulator (MadR), a transcriptional repressor of the mycolate desaturase genes desA1 and desA2, controls mycolic acid desaturation and biosynthesis in response to cell envelope stress. A madR-null mutant of M. smegmatis exhibited traits of an impaired cell wall with an altered outer mycomembrane, accumulation of a desaturated α-mycolate, susceptibility to antimycobacterials, and cell surface disruption. Transcriptomic profiling showed that enriched lipid metabolism genes that were significantly down-regulated upon madR deletion included acyl-coenzyme A (aceyl-CoA) dehydrogenases, implicating it in the indirect control of β-oxidation pathways. Electromobility shift assays and binding affinities suggest a unique acyl-CoA pool-sensing mechanism, whereby MadR is able to bind a range of acyl-CoAs, including those with unsaturated as well as saturated acyl chains. MadR repression of desA1/desA2 is relieved upon binding of saturated acyl-CoAs of chain length C16 to C24, while no impact is observed upon binding of shorter chain and unsaturated acyl-CoAs. We propose this mechanism of regulation as distinct to other mycolic acid and fatty acid synthesis regulators and place MadR as the key regulatory checkpoint that coordinates mycolic acid remodeling during infection in response to host-derived cell surface perturbation.

Rv0500A is a transcription factor that links Mycobacterium tuberculosis environmental response with division and impacts host colonization

For Mycobacterium tuberculosis (Mtb) to successfully infect a host, it must be able to adapt to changes in its microenvironment, including to variations in ionic signals such as pH and chloride (Cl^- ), and link these responses to its growth. Transcriptional changes are a key mechanism for Mtb environmental adaptation, and we identify here Rv0500A as a novel transcriptional regulator that links Mtb environmental response and division processes. Global transcriptional profiling revealed that Rv0500A acts as a repressor and influences the expression of genes related to division, with the magnitude of its effect modulated by pH and Cl^- . Rv0500A can directly bind the promoters of several of these target genes, and we identify key residues required for its DNA-binding ability and biological effect. Overexpression of rv0500A disrupted Mtb growth morphology, resulting in filamentation that was exacerbated by high environmental Cl^- levels and acidic pH. Finally, we show that perturbation of rv0500A leads to attenuation of the ability of Mtb to colonize its host in vivo. Our work highlights the important link between Mtb environmental response and growth characteristics, and uncovers a new transcription factor involved in this critical facet of Mtb biology.

Host-pathogen genetic interactions underlie tuberculosis susceptibility in genetically diverse mice

The outcome of an encounter with Mycobacterium tuberculosis (Mtb) depends on the pathogen’s ability to adapt to the variable immune pressures exerted by the host. Understanding this interplay has proven difficult, largely because experimentally tractable animal models do not recapitulate the heterogeneity of tuberculosis disease. We leveraged the genetically diverse Collaborative Cross (CC) mouse panel in conjunction with a library of Mtb mutants to create a resource for associating bacterial genetic requirements with host genetics and immunity. We report that CC strains vary dramatically in their susceptibility to infection and produce qualitatively distinct immune states. Global analysis of Mtb transposon mutant fitness (TnSeq) across the CC panel revealed that many virulence pathways are only required in specific host microenvironments, identifying a large fraction of the pathogen’s genome that has been maintained to ensure fitness in a diverse population. Both immunological and bacterial traits can be associated with genetic variants distributed across the mouse genome, making the CC a unique population for identifying specific host-pathogen genetic interactions that influence pathogenesis.

Proteome Profiling of Mycobacterium tuberculosis Cells Exposed to Nitrosative Stress

Reactive nitrogen species (RNS) are secreted by human cells in response to infection by Mycobacterium tuberculosis (Mtb). Although RNS can kill Mtb under some circumstances, Mtb can adapt and survive in the presence of RNS by a process that involves modulation of gene expression. Previous studies focused primarily on stress-related changes in the Mtb transcriptome. This study unveils changes in the Mtb proteome in response to a sub-lethal dose of nitric oxide (NO) over several hours of exposure. Proteins were identified using liquid chromatography coupled with electrospray ionization mass spectrometry (LC-MS/MS). A total of 2911 Mtb proteins were identified, of which 581 were differentially abundant (DA) after exposure to NO in at least one of the four time points (30 min, 2 h, 6 h, and 20 h). The proteomic response to NO was marked by two phases, with few DA proteins in the early phase and a multitude of DA proteins in the later phase. The efflux pump Rv1687 stood out as being the only protein more abundant at all the time points and might play a role in the early protection of Mtb against nitrosative stress. These changes appeared to be compensatory in nature, contributing to iron homeostasis, energy metabolism, and other stress responses. This study thereby provides new insights into the response of Mtb to NO at the level of proteomics.

Biosynthesis of the redox cofactor mycofactocin is controlled by the transcriptional regulator MftR and induced by long-chain acyl-CoA species

Mycofactocin (MFT) is a ribosomally synthesized and post-translationally-modified redox cofactor found in pathogenic mycobacteria. While MFT biosynthetic proteins have been extensively characterized, the physiological conditions under which MFT biosynthesis is required are not well understood. To gain insights into the mechanisms of regulation of MFT expression in Mycobacterium smegmatis mc²155, we investigated the DNA-binding and ligand-binding activities of the putative TetR-like transcription regulator, MftR. In this study, we demonstrated that MftR binds to the mft promoter region. We used DNase I footprinting to identify the 27 bp palindromic operator located 5′ to mftA and found it to be highly conserved in Mycobacterium tuberculosis, Mycobacterium bovis, Mycobacterium ulcerans, and Mycobacterium marinum. To determine under which conditions the mft biosynthetic gene cluster (BGC) is induced, we screened for effectors of MftR. As a result, we found that MftR binds to long-chain acyl-CoAs with low micromolar affinities. To demonstrate that oleoyl-CoA induces the mft BGC in vivo, we re-engineered a fluorescent protein reporter system to express an MftA–mCherry fusion protein. Using this mCherry fluorescent readout, we show that the mft BGC is upregulated in M. smegmatis mc²155 when oleic acid is supplemented to the media. These results suggest that MftR controls expression of the mft BGC and that MFT production is induced by long-chain acyl-CoAs. Since MFT-dependent dehydrogenases are known to colocalize with acyl carrier protein/CoA-modifying enzymes, these results suggest that MFT might be critical for fatty acid metabolism or cell wall reorganization.

Quantitative prediction of conditional vulnerabilities in regulatory and metabolic networks using PRIME

The ability of Mycobacterium tuberculosis (Mtb) to adopt heterogeneous physiological states underlies its success in evading the immune system and tolerating antibiotic killing. Drug tolerant phenotypes are a major reason why the tuberculosis (TB) mortality rate is so high, with over 1.8 million deaths annually. To develop new TB therapeutics that better treat the infection (faster and more completely), a systems-level approach is needed to reveal the complexity of network-based adaptations of Mtb. Here, we report a new predictive model called PRIME (Phenotype of Regulatory influences Integrated with Metabolism and Environment) to uncover environment-specific vulnerabilities within the regulatory and metabolic networks of Mtb. Through extensive performance evaluations using genome-wide fitness screens, we demonstrate that PRIME makes mechanistically accurate predictions of context-specific vulnerabilities within the integrated regulatory and metabolic networks of Mtb, accurately rank-ordering targets for potentiating treatment with frontline drugs.

Methylation Motifs in Promoter Sequences May Contribute to the Maintenance of a Conserved m5C Methyltransferase in Helicobacter pylori

DNA methylomes of Helicobacter pylori strains are complex due to the large number of DNA methyltransferases (MTases) they possess. H. pylori J99 M.Hpy99III is a 5-methylcytosine (^m5C) MTase that converts GCGC motifs to G^m5CGC. Homologs of M.Hpy99III are found in essentially all H. pylori strains. Most of these homologs are orphan MTases that lack a cognate restriction endonuclease, and their retention in H. pylori strains suggest they have roles in gene regulation. To address this hypothesis, green fluorescent protein (GFP) reporter genes were constructed with six putative promoters that had a GCGC motif in the extended −10 region, and the expression of the reporter genes was compared in wild-type H. pylori G27 and a mutant lacking the M.Hpy99III homolog (M.HpyGIII). The expression of three of the GFP reporter genes was decreased significantly in the mutant lacking M.HpyGIII. In addition, the growth rate of the H. pylori G27 mutant lacking M.HpyGIII was reduced markedly compared to that of the wild type. These findings suggest that the methylation of the GCGC motif in many H. pylori GCGC-containing promoters is required for the robust expression of genes controlled by these promoters, which may account for the universal retention of M.Hpy99III homologs in H. pylori strains.

A balancing act in transcription regulation by response regulators: titration of transcription factor activity by decoy DNA binding sites

Studies of transcription regulation are often focused on binding of transcription factors (TFs) to a small number of promoters of interest. It is often assumed that TFs are in great excess to their binding sites (TFBSs) and competition for TFs between DNA sites is seldom considered. With increasing evidence that TFBSs are exceedingly abundant for many TFs and significant variations in TF and TFBS numbers occur during growth, the interplay between a TF and all TFBSs should not be ignored. Here, we use additional decoy DNA sites to quantitatively analyze how the relative abundance of a TF to its TFBSs impacts the steady-state level and onset time of gene expression for the auto-activated Escherichia coli PhoB response regulator. We show that increasing numbers of decoy sites progressively delayed transcription activation and lowered promoter activities. Perturbation of transcription regulation by additional TFBSs did not require extreme numbers of decoys, suggesting that PhoB is approximately at capacity for its DNA sites. Addition of decoys also converted a graded response to a bi-modal response. We developed a binding competition model that captures the major features of experimental observations, providing a quantitative framework to assess how variations in TFs and TFBSs influence transcriptional responses.

OsnR is an autoregulatory negative transcription factor controlling redox-dependent stress responses in Corynebacterium glutamicum

Background

Corynebacterium glutamicum is used in the industrial production of amino acids and nucleotides. During the course of fermentation, C. glutamicum cells face various stresses and employ multiple regulatory genes to cope with the oxidative stress. The osnR gene plays a negative regulatory role in redox-dependent oxidative-stress responses, but the underlying mechanism is not known yet.

Results

Overexpression of the osnR gene in C. glutamicum affected the expression of genes involved in the mycothiol metabolism. ChIP-seq analysis revealed that OsnR binds to the promoter region of multiple genes, including osnR and cg0026, which seems to function in the membrane-associated redox metabolism. Studies on the role of the osnR gene involving in vitro assays employing purified OsnR proteins and in vivo physiological analyses have identified that OsnR inhibits the transcription of its own gene. Further, oxidant diamide stimulates OsnR-binding to the promoter region of the osnR gene. The genes affected by the overexpression of osnR have been found to be under the control of σ^H. In the osnR-overexpressing strain, the transcription of sigH is significantly decreased and the stimulation of sigH transcription by external stress is lost, suggesting that osnR and sigH form an intimate regulatory network.

Conclusions

Our study suggests that OsnR not only functions as a transcriptional repressor of its own gene and of those involved in redox-dependent stress responses but also participates in the global transcriptional regulation by controlling the transcription of other master regulators, such as sigH.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01693-1.

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacterium tuberculosis

Nucleoid-associated proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacterium tuberculosis. Results suggest narrower minor groove width (MGW) and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM have heterogeneous DNA structural predilections. In contrast, WhiB4-DNA-binding sites were characterized by wider MGW, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

Expression Dysregulation as a Mediator of Fitness Costs in Antibiotic Resistance

Antimicrobial resistance (AMR) poses a threat to global health and the economy. Rifampicin-resistant Mycobacterium tuberculosis accounts for a third of the global AMR burden. Gaining the upper hand on AMR requires a deeper understanding of the physiology of resistance. AMR often results in a fitness cost in the absence of drug. Identifying the molecular mechanisms underpinning this cost could help strengthen future treatment regimens. Here, we used a collection of M. tuberculosis strains that provide an evolutionary and phylogenetic snapshot of rifampicin resistance and subjected them to genome-wide transcriptomic and proteomic profiling to identify key perturbations of normal physiology. We found that the clinically most common rifampicin resistance-conferring mutation, RpoB Ser450Leu, imparts considerable gene expression changes, many of which are mitigated by the compensatory mutation in RpoC Leu516Pro. However, our data also provide evidence for pervasive epistasis—the same resistance mutation imposed a different fitness cost and functionally distinct changes to gene expression in genetically unrelated clinical strains. Finally, we report a likely posttranscriptional modulation of gene expression that is shared in most of the tested strains carrying RpoB Ser450Leu, resulting in an increased abundance of proteins involved in central carbon metabolism. These changes contribute to a more general trend in which the disruption of the composition of the proteome correlates with the fitness cost of the RpoB Ser450Leu mutation in different strains.

Nucleoid-associated proteins shape chromatin structure and transcriptional regulation across the bacterial kingdom

Genome architecture has proven to be critical in determining gene regulation across almost all domains of life. While many of the key components and mechanisms of eukaryotic genome organization have been described, the interplay between bacterial DNA organization and gene regulation is only now being fully appreciated. An increasing pool of evidence has demonstrated that the bacterial chromosome can reasonably be thought of as chromatin, and that bacterial chromosomes contain transcriptionally silent and transcriptionally active regions analogous to heterochromatin and euchromatin, respectively. The roles played by histones in eukaryotic systems appear to be shared across a range of nucleoid-associated proteins (NAPs) in bacteria, which function to compact, structure, and regulate large portions of bacterial chromosomes. The broad range of extant NAPs, and the extent to which they differ from species to species, has raised additional challenges in identifying and characterizing their roles in all but a handful of model bacteria. Here we review the regulatory roles played by NAPs in several well-studied bacteria and use the resulting state of knowledge to provide a working definition for NAPs, based on their function, binding pattern, and expression levels. We present a screening procedure which can be applied to any species for which transcriptomic data are available. Finally, we note that NAPs tend to play two major regulatory roles - xenogeneic silencers and developmental regulators - and that many unrecognized potential NAPs exist in each bacterial species examined.

Mycobacterium tuberculosis SufR responds to nitric oxide via its 4Fe-4S cluster and regulates Fe-S cluster biogenesis for persistence in mice

The persistence of Mycobacterium tuberculosis (Mtb) is a major problem in managing tuberculosis (TB). Host-generated nitric oxide (NO) is perceived as one of the signals by Mtb to reprogram metabolism and respiration for persistence. However, the mechanisms involved in NO sensing and reorganizing Mtb's physiology are not fully understood. Since NO damages iron-sulfur (Fe-S) clusters of essential enzymes, the mechanism(s) involved in regulating Fe-S cluster biogenesis could help Mtb persist in host tissues. Here, we show that a transcription factor SufR (Rv1460) senses NO via its 4Fe-4S cluster and promotes persistence of Mtb by mobilizing the Fe-S cluster biogenesis system; suf operon (Rv1460-Rv1466). Analysis of anaerobically purified SufR by UV-visible spectroscopy, circular dichroism, and iron-sulfide estimation confirms the presence of a 4Fe-4S cluster. Atmospheric O₂ and H₂O₂ gradually degrade the 4Fe-4S cluster of SufR. Furthermore, electron paramagnetic resonance (EPR) analysis demonstrates that NO directly targets SufR 4Fe-4S cluster by forming a protein-bound dinitrosyl-iron-dithiol complex. DNase I footprinting, gel-shift, and in vitro transcription assays confirm that SufR directly regulates the expression of the suf operon in response to NO. Consistent with this, RNA-sequencing of MtbΔsufR demonstrates deregulation of the suf operon under NO stress. Strikingly, NO inflicted irreversible damage upon Fe-S clusters to exhaust respiratory and redox buffering capacity of MtbΔsufR. Lastly, MtbΔsufR failed to recover from a NO-induced non-growing state and displayed persistence defect inside immune-activated macrophages and murine lungs in a NO-dependent manner. Data suggest that SufR is a sensor of NO that supports persistence by reprogramming Fe-S cluster metabolism and bioenergetics.

Collateral Sensitivity to β-Lactam Drugs in Drug-Resistant Tuberculosis Is Driven by the Transcriptional Wiring of BlaI Operon Genes

The evolution of resistance to one antimicrobial can result in enhanced sensitivity to another, known as "collateral sensitivity." This underexplored phenomenon opens new therapeutic possibilities for patients infected with pathogens unresponsive to classical treatments. Intrinsic resistance to β-lactams in Mycobacterium tuberculosis (the causative agent of tuberculosis) has traditionally curtailed the use of these low-cost and easy-to-administer drugs for tuberculosis treatment. Recently, β-lactam sensitivity has been reported in strains resistant to classical tuberculosis therapy, resurging the interest in β-lactams for tuberculosis. However, a lack of understanding of the molecular underpinnings of this sensitivity has delayed exploration in the clinic. We performed gene expression and network analyses and in silico knockout simulations of genes associated with β-lactam sensitivity and genes associated with resistance to classical tuberculosis drugs to investigate regulatory interactions and identify key gene mediators. We found activation of the key inhibitor of β-lactam resistance, blaI, following classical drug treatment as well as transcriptional links between genes associated with β-lactam sensitivity and those associated with resistance to classical treatment, suggesting that regulatory links might explain collateral sensitivity to β-lactams. Our results support M. tuberculosis β-lactam sensitivity as a collateral consequence of the evolution of resistance to classical tuberculosis drugs, mediated through changes to transcriptional regulation. These findings support continued exploration of β-lactams for the treatment of patients infected with tuberculosis strains resistant to classical therapies. IMPORTANCE Tuberculosis remains a significant cause of global mortality, with strains resistant to classical drug treatment considered a major health concern by the World Health Organization. Challenging treatment regimens and difficulty accessing drugs in low-income communities have led to a high prevalence of strains resistant to multiple drugs, making the development of alternative therapies a priority. Although Mycobacterium tuberculosis is naturally resistant to β-lactam drugs, previous studies have shown sensitivity in strains resistant to classical drug treatment, but we currently lack understanding of the molecular underpinnings behind this phenomenon. We found that genes involved in β-lactam susceptibility are activated after classical drug treatment resulting from tight regulatory links with genes involved in drug resistance. Our study supports the hypothesis that β-lactam susceptibility observed in drug-resistant strains results from the underlying regulatory network of M. tuberculosis, supporting further exploration of the use of β-lactams for tuberculosis treatment.

Intricate Genetic Programs Controlling Dormancy in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) displays the remarkable ability to transition in and out of dormancy, a hallmark of the pathogen’s capacity to evade the immune system and exploit susceptible individuals. Uncovering the gene regulatory programs that underlie the phenotypic shifts in MTB during disease latency and reactivation has posed a challenge. We develop an experimental system to precisely control dissolved oxygen levels in MTB cultures in order to capture the transcriptional events that unfold as MTB transitions into and out of hypoxia-induced dormancy. Using a comprehensive genome-wide transcription factor binding map and insights from network topology analysis, we identify regulatory circuits that deterministically drive sequential transitions across six transcriptionally and functionally distinct states encompassing more than three-fifths of the MTB genome. The architecture of the genetic programs explains the transcriptional dynamics underlying synchronous entry of cells into a dormant state that is primed to infect the host upon encountering favorable conditions.

Mycobacterium tuberculosis (MTB) persists within the host by counteracting disparate stressors including hypoxia. Peterson et al. report a transcriptional program that coordinates sequential state transitions to drive MTB in and out of hypoxia-induced dormancy. Among varied properties, this program encodes advanced preparedness to infect the host in favorable conditions.

Dissecting the RecA-(In)dependent Response to Mitomycin C in Mycobacterium tuberculosis Using Transcriptional Profiling and Proteomics Analyses

Mycobacteria exploit at least two independent global systems in response to DNA damage: the LexA/RecA-dependent SOS response and the PafBC-regulated pathway. Intracellular pathogens, such as Mycobacterium tuberculosis, are exposed to oxidative and nitrosative stress during the course of infection while residing inside host macrophages. The current understanding of RecA-independent responses to DNA damage is based on the saprophytic model of Mycobacterium smegmatis, a free-living and nonpathogenic mycobacterium. The aim of the present study was to identify elements of RecA-independent responses to DNA damage in pathogenic intracellular mycobacteria. With the help of global transcriptional profiling, we were able to dissect RecA-dependent and RecA-independent pathways. We profiled the DNA damage responses of an M. tuberculosis strain lacking the recA gene, a strain with an undetectable level of the PafBC regulatory system, and a strain with both systems tuned down simultaneously. RNA-Seq profiling was correlated with the evaluation of cell survival in response to DNA damage to estimate the relevance of each system to the overall sensitivity to genotoxic agents. We also carried out whole-cell proteomics analysis of the M. tuberculosis strains in response to mitomycin C. This approach highlighted that LexA, a well-defined key element of the SOS system, is proteolytically inactivated during RecA-dependent DNA repair, which we found to be transcriptionally repressed in response to DNA-damaging agents in the absence of RecA. Proteomics profiling revealed that AlkB was significantly overproduced in the ΔrecA pafBC^{CRISPRi/dCas9} strain and that Holliday junction resolvase RuvX was a DNA damage response factor that was significantly upregulated regardless of the presence of functional RecA and PafBC systems, thus falling into a third category of DNA damage factors: RecA- and PafBC-independent. While invisible to the mass spectrometer, the genes encoding alkA, dnaB, and dnaE2 were significantly overexpressed in the ΔrecA pafBC^{CRISPRi/dCas9} strain at the transcript level.