Identification of operator sites within the upstream region of the putativemce2Rgene from mycobacteria

PcaR, a GntR/FadR Family Transcriptional Repressor Controls the Transcription of Phenazine-1-Carboxylic Acid 1,2-Dioxygenase Gene Cluster in Sphingomonas histidinilytica DS-9

In our previous study, the phenazine-1-carboxylic acid (PCA) 1,2-dioxygenase gene cluster (pcaA1A2A3A4 cluster) in Sphingomonas histidinilytica DS-9 was identified to be responsible for the conversion of PCA to 1,2-dihydroxyphenazine (Ren Y, Zhang M, Gao S, Zhu Q, et al. 2022. Appl Environ Microbiol 88:e00543-22). However, the regulatory mechanism of the pcaA1A2A3A4 cluster has not been elucidated yet. In this study, the pcaA1A2A3A4 cluster was found to be transcribed as two divergent operons: pcaA3-ORF5205 (named A3-5205 operon) and pcaA1A2-ORF5208-pcaA4-ORF5210 (named A1-5210 operon). The promoter regions of the two operons were overlapped. PcaR acts as a transcriptional repressor of the pcaA1A2A3A4 cluster, and it belongs to GntR/FadR family transcriptional regulator. Gene disruption of pcaR can shorten the lag phase of PCA degradation. The results of electrophoretic mobility shift assay and DNase I footprinting showed that PcaR binds to a 25-bp motif in the ORF5205-pcaA1 intergenic promoter region to regulate the expression of two operons. The 25-bp motif covers the -10 region of the promoter of A3-5205 operon and the -35 region and -10 region of the promoter of A1-5210 operon. The TNGT/ANCNA box within the motif was essential for PcaR binding to the two promoters. PCA acted as an effector of PcaR, preventing it from binding to the promoter region and repressing the transcription of the pcaA1A2A3A4 cluster. In addition, PcaR represses its own transcription, and this repression can be relieved by PCA. This study reveals the regulatory mechanism of PCA degradation in strain DS-9, and the identification of PcaR increases the variety of regulatory model of the GntR/FadR-type regulator. IMPORTANCE Sphingomonas histidinilytica DS-9 is a phenazine-1-carboxylic acid (PCA)-degrading strain. The 1,2-dioxygenase gene cluster (pcaA1A2A3A4 cluster, encoding dioxygenase PcaA1A2, reductase PcaA3, and ferredoxin PcaA4) is responsible for the initial degradation step of PCA and widely distributed in Sphingomonads, but its regulatory mechanism has not been investigated yet. In this study, a GntR/FadR-type transcriptional regulator PcaR repressing the transcription of pcaA1A2A3A4 cluster and pcaR gene was identified and characterized. The binding site of PcaR in ORF5205-pcaA1 intergenic promoter region contains a TNGT/ANCNA box, which is important for the binding. These findings enhance our understanding of the molecular mechanism of PCA degradation.

Mycobacterial MCE proteins as transporters that control lipid homeostasis of the cell wall

Mammalian cell entry (mce) genes are not only present in genomes of pathogenic mycobacteria, including Mycobacterium tuberculosis (the causative agent of tuberculosis), but also in saprophytic and opportunistic mycobacterial species. MCE are conserved cell-wall proteins encoded by mce operons, which maintain an identical structure in all mycobacteria: two yrbE genes (A and B) followed by six mce genes (A, B, C, D, E and F). Although these proteins are known to participate in the virulence of pathogenic mycobacteria, the presence of the operons in nonpathogenic mycobacteria and other bacteria indicates that they play another role apart from host cell invasion. In this respect, more recent studies suggest that they are functionally similar to ABC transporters and form part of lipid transporters in Actinobacteria. To date, most reviews on mce operons in the literature discuss their role in virulence. However, according to data from transcriptional studies, mce genes, particularly the mce1 and mce4 operons, modify their expression according to the carbon source and upon hypoxia, starvation, surface stress and oxidative stress; which suggests a role of MCE proteins in the response of Mycobacteria to external stressors. In addition to these data, this review also summarizes the studies demonstrating the role of MCE proteins as lipid transporters as well as the relevance of their transport function in the interaction of pathogenic Mycobacteria with the hosts. Altogether, the evidence to date would indicate that MCE proteins participate in the response to the stress conditions that mycobacteria encounter during infection, by participating in the cell wall remodelling and possibly contributing to lipid homeostasis.

Mce2R/Rv0586 of Mycobacterium tuberculosis is the functional homologue of FadRE. coli

Lipid metabolism is critical to Mycobacterium tuberculosis survival and infection. Unlike Escherichia coli, which has a single FadR, the M. tuberculosis genome encodes five proteins of the FadR sub-family. While the role of E. coli FadR as a regulator of fatty acid metabolism is well known, the definitive functions of M. tuberculosis FadR proteins are still under investigation. An interesting question about the M. tuberculosis FadRs remains open: which one of these proteins is the functional homologue of E. coli FadR? To address this, we have applied two different approaches. The first one was the bioinformatics approach and the second one was the classical molecular genetic approach involving complementation studies. Surprisingly, the results of these two approaches did not agree. Among the five M. tuberculosis FadRs, Rv0494 shared the highest sequence similarity with FadR^{E. coli} and Rv0586 was the second best match. However, only Rv0586, but not Rv0494, could complement E. coli ∆fadR, indicating that Rv0586 is the M. tuberculosis functional homologue of FadR^{E. coli}. Further studies showed that both regulators, Rv0494 and Rv0586, show similar responsiveness to LCFA, and have conserved critical residues for DNA binding. However, analysis of the operator site indicated that the inter-palindromic distance required for DNA binding differs for the two regulators. The differences in the binding site selection helped in the success of Rv0586 binding to fadB upstream over Rv0494 and may have played a critical role in complementing E. coli ∆fadR. Further, for the first time, we report the lipid-responsive nature of Rv0586.

ClosIndb: A resource for computationally derived information from clostridial genomes

Over the past few years, several clostridial genomes have been sequenced, and since then new sequencing projects are also under way. Clostridia is one of the most sequenced genera, and presently, complete genome sequences of 49 clostridial species are available in public archives. Unraveling this wealth of genomic information opens up potential avenues in clostridial research. In the present study, we have carried out in silico analysis to decipher the genomic data. Subsequently, a web resource, ClosIndb, has been developed which collates the computationally derived information associated with all clostridial genes. It features various aspects of coding regions as well as non-coding regions, such as putative orthologs, proteins physicochemical properties, operons and cis-regulatory elements. It provides users with comparative details of all clostridial proteins across the firmicutes. ClosIndb is a comprehensive resource for all completely sequenced clostridial genomes and is under constant development. ClosIndb is freely accessible at http://bif.uohyd.ac.in/closindb/.

Genetic heterogeneity revealed by sequence analysis of Mycobacterium tuberculosis isolates from extra-pulmonary tuberculosis patients

BACKGROUND

Tuberculosis remains a major public health problem. Clinical tuberculosis manifests often as pulmonary and occasionally as extra-pulmonary tuberculosis. The emergence of drug resistant tubercle bacilli and its association with HIV is a formidable challenge to curb the spread of tuberculosis. There have been concerted efforts by whole genome sequencing and bioinformatics analysis to identify genomic patterns and to establish a relationship between the genotype of the organism and clinical manifestation of tuberculosis. Extra-pulmonary TB constitutes 15-20 percent of the total clinical cases of tuberculosis reported among immunocompetent patients, whereas among HIV patients the incidence is more than 50 percent. Genomic analysis of M. tuberculosis isolates from extra pulmonary patients has not been explored.

RESULTS

The genomic DNA of 5 extra-pulmonary clinical isolates of M. tuberculosis derived from cerebrospinal fluid, lymph node fine needle aspirates (FNAC) / biopsies, were sequenced. Next generation sequencing approach (NGS) was employed to identify Single Nucleotide Variations (SNVs) and computational methods used to predict their consequence on functional genes. Analysis of distribution of SNVs led to the finding that there are mixed genotypes in patient isolates and that many SNVs are likely to influence either gene function or their expression. Phylogenetic relationship between the isolates correlated with the origin of the isolates. In addition, insertion sites of IS elements were identified and their distribution revealed a variation in number and position of the element in the 5 extra-pulmonary isolates compared to the reference M. tuberculosis H37Rv strain.

CONCLUSIONS

The results suggest that NGS sequencing is able to identify small variations in genomes of M. tuberculosis isolates including changes in IS element insertion sites. Moreover, variations in isolates of M. tuberculosis from non-pulmonary sites were documented. The analysis of our results indicates genomic heterogeneity in the clinical isolates.

Identification and characterization of Rv0494: a fatty acid-responsive protein of the GntR/FadR family from Mycobacterium tuberculosis

Escherichia coli FadR, a member of the GntR family of transcription factors, plays dual roles in fatty acid metabolism. FadR-DNA binding is inhibited by fatty acyl-CoAs, and thus FadR acts as a sensor of the fatty acid level in bacteria. We have identified FadR-binding sites in the upstream regions of genes showing altered expression after the disruption of fatty acid biosynthesis in Mycobacterium tuberculosis. A FadR homologue in M. tuberculosis, Rv0494, was identified, which binds to its operator in the upstream region of the kas operon. We have shown that FadRMt (Rv0494) directly binds to long-chain fatty acyl-CoA and that binding quenches the intrinsic fluorescence of the purified protein. FadR-DNA binding can be impaired by long-chain fatty acyl-CoA compounds. Overexpression of Rv0494 in Mycobacterium bovis BCG reduced the basal level expression of kas operon genes, thereby suggesting the repressor nature of this protein in fatty acid synthase II regulation. This is the first report, to the best of our knowledge, of a GntR/FadR family protein acting as a fatty acid-responsive transcriptional regulator in M. tuberculosis, suggesting a possible role for this protein in mycolic acid biosynthesis.

MycoRRdb: a database of computationally identified regulatory regions within intergenic sequences in mycobacterial genomes

The identification of regulatory regions for a gene is an important step towards deciphering the gene regulation. Regulatory regions tend to be conserved under evolution that facilitates the application of comparative genomics to identify such regions. The present study is an attempt to make use of this attribute to identify regulatory regions in the Mycobacterium species followed by the development of a database, MycoRRdb. It consist the regulatory regions identified within the intergenic distances of 25 mycobacterial species. MycoRRdb allows to retrieve the identified intergenic regulatory elements in the mycobacterial genomes. In addition to the predicted motifs, it also allows user to retrieve the Reciprocal Best BLAST Hits across the mycobacterial genomes. It is a useful resource to understand the transcriptional regulatory mechanism of mycobacterial species. This database is first of its kind which specifically addresses cis-regulatory regions and also comprehensive to the mycobacterial species. Database URL: http://mycorrdb.uohbif.in.

Functional analysis of an intergenic non-coding sequence within mce1 operon of M.tuberculosis

Background

The mce operons play an important role in the entry of M. tuberculosis into macrophages and non-phagocytic cells. Their non-redundant function as well as complex regulation is implied by the phenotype of mce mutants. Recently, mce1 operon was found to extend over 13 genes, fadD5 (Rv0166) being the first gene of the operon. The presence of a non-coding sequence of 200 base pairs between Rv0166 and Rv0167 is peculiar to mce1 among the four mce operons of M.tuberculosis. We have examined the function of this region.

Results

We predicted putative promoter activity of the 200 base pairs of non-coding, intergenic region between Rv0166 and Rv0167 in silico using MEME software and designate it as intergenic promoter, IGPr. We demonstrate both promoter activity and a putative negative regulatory function of this fragment by reporter assays carried out in the surrogate host M.smegmatis. We find that the repressive elements not only control the native promoter but also repress a heterologous promoter of M.smegmatis. The higher activity of the intergenic promoter in a clinical isolate in comparison with the wild type sequence from M.tuberculosis H37Rv could be correlated with a point mutation within the negative element. We have mapped two transcription start sites for mce1 operon both of which are utilized in M.tuberculosis H37Rv as well as the clinical isolate VPCI591. Our studies show that the promoter activity in the non-coding region is relevant not only in reporter gene expression but also in the expression of mce1 operon in M. tuberculosis cells grown in synthetic medium.

Conclusion

The mce operon of M.tuberculosis H37Rv potentially can be transcribed from two promoters P1 and P2, former mapping upstream of Rv0166 and the latter in the non-coding intergenic region between Rv0166 and Rv0167. The transcription initiation from P1 results in a transcript with Rv0166 while that from P2 will be without it. The sequences between the translation start site of Rv0167 and the promoter P2 have a negative regulatory role, as point mutation within the sequence leads to enhanced activity of P2 as well as a heterologous promoter from M.smegmatis. The mutation detected in the clinical isolate VPCI591 therefore behaves like a gain-of-function mutation.

Mce2R from Mycobacterium tuberculosis represses the expression of the mce2 operon

The mce2 operon is one of the four mce operons present in Mycobacterium tuberculosis that encode exported proteins with a probable role in the virulence mechanisms of this bacterium. In the present study we demonstrated that Rv0586, which encodes a putative GntR-like regulator, is part of the mce2 operon. By using a promoter-lacZ fusion approach and bioinformatics tools, we found that Rv0586 represses the expression of Mce2 proteins and of a putative endonuclease IV, encoded by end (Rv0670) gene. For this reason, we have re-named the repressor protein Mce2R. By gel-shift experiments Mce2R binding was determined to be located within the mce2 promoter region. In addition, two FadR-like operator motifs were identified within the promoter regions of both the mce2 operon and the end gene. These motifs overlap putative -10 and -35 promoter boxes. M. tuberculosis carrying mce2 and end promoter-lacZ fusions were used to infect J774 macrophage-like cells. Expression of beta-galactosidase was induced after phagocytocis, suggesting that some cellular factor could be a key component of the molecular switch regulation expression of the mce2 operon. In conclusion, these results add novel evidence of the complex regulation of mce operon expression.