A predicted operon map for Mycobacterium tuberculosis

The prediction of operons in Mycobacterium tuberculosis (MTB) is a first step toward understanding the regulatory network of this pathogen. Here we apply a statistical model using logistic regression to predict operons in MTB. As predictors, our model incorporates intergenic distance and the correlation of gene expression calculated for adjacent gene pairs from over 474 microarray experiments with MTB RNA. We validate our findings with known examples from the literature and experimentation. From this model, we rank each potential operon pair by the strength of evidence for cotranscription, choose a classification threshold with a true positive rate of over 90% at a false positive rate of 9.1%, and use it to construct an operon map for the MTB genome

The katG mRNA of Mycobacterium tuberculosis and Mycobacterium smegmatis is processed at its 5' end and is stabilized by both a polypurine sequence and translation initiation

<p>Abstract</p> <p>Background</p> <p>In <it>Mycobacterium tuberculosis </it>and in <it>Mycobacterium smegmatis </it>the <it>furA</it>-<it>katG </it>loci, encoding the FurA regulatory protein and the KatG catalase-peroxidase, are highly conserved. In <it>M. tuberculosis furA-katG </it>constitute a single operon, whereas in <it>M. smegmatis </it>a single mRNA covering both genes could not be found. In both species, specific 5' ends have been identified: the first one, located upstream of the <it>furA </it>gene, corresponds to transcription initiation from the <it>furA </it>promoter; the second one is the <it>katG </it>mRNA 5' end, located in the terminal part of <it>furA</it>.</p> <p>Results</p> <p>In this work we demonstrate by in vitro transcription and by RNA polymerase Chromatin immunoprecipitation that no promoter is present in the <it>M. smegmatis </it>region covering the latter 5' end, suggesting that it is produced by specific processing of longer transcripts. Several DNA fragments of <it>M. tuberculosis </it>and <it>M. smegmatis </it>were inserted in a plasmid between the <it>sigA </it>promoter and the <it>lacZ </it>reporter gene, and expression of the reporter gene was measured. A polypurine sequence, located four bp upstream of the <it>katG </it>translation start codon, increased beta-galactosidase activity and stabilized the <it>lacZ </it>transcript. Mutagenesis of this sequence led to destabilization of the mRNA. Analysis of constructs, in which the polypurine sequence of <it>M. smegmatis </it>was followed by an increasing number of <it>katG </it>codons, demonstrated that mRNA stability requires translation of at least 20 amino acids. In order to define the requirements for the 5' processing of the <it>katG </it>transcript, we created several mutations in this region and analyzed the 5' ends of the transcripts: the distance from the polypurine sequence does not seem to influence the processing, neither the sequence around the cutting point. Only mutations which create a double stranded region around the processing site prevented RNA processing.</p> <p>Conclusion</p> <p>This is the first reported case in mycobacteria, in which both a polypurine sequence and translation initiation are shown to contribute to mRNA stability. The <it>furA-katG </it>mRNA is transcribed from the <it>furA </it>promoter and immediately processed; this processing is prevented by a double stranded RNA at the cutting site, suggesting that the endoribonuclease responsible for the cleavage cuts single stranded RNA.</p

A Regulatory Role for NBS1 in Strand-Specific Mutagenesis during Somatic Hypermutation

Activation-induced cytidine deaminase (AID) is believed to initiate somatic hypermutation (SHM) by deamination of deoxycytidines to deoxyuridines within the immunoglobulin variable regions genes. The deaminated bases can subsequently be replicated over, processed by base excision repair or mismatch repair, leading to introduction of different types of point mutations (G/C transitions, G/C transversions and A/T mutations). It is evident that the base excision repair pathway is largely dependent on uracil-DNA glycosylase (UNG) through its uracil excision activity. It is not known, however, which endonuclease acts in the step immediately downstream of UNG, i.e. that cleaves at the abasic sites generated by the latter. Two candidates have been proposed, an apurinic/apyrimidinic endonuclease (APE) and the Mre11-Rad50-NBS1 complex. The latter is intriguing as this might explain how the mutagenic pathway is primed during SHM. We have investigated the latter possibility by studying the in vivo SHM pattern in B cells from ataxia-telangiectasia-like disorder (Mre11 deficient) and Nijmegen breakage syndrome (NBS1 deficient) patients. Our results show that, although the pattern of mutations in the variable heavy chain (VH) genes was altered in NBS1 deficient patients, with a significantly increased number of G (but not C) transversions occurring in the SHM and/or AID targeting hotspots, the general pattern of mutations in the VH genes in Mre11 deficient patients was only slightly altered, with an increased frequency of A to C transversions. The Mre11-Rad50-NBS1 complex is thus unlikely to be the major nuclease involved in cleavage of the abasic sites during SHM, whereas NBS1 might have a specific role in regulating the strand-biased repair during phase Ib mutagenesis

Mycolactone Gene Expression Is Controlled by Strong SigA-Like Promoters with Utility in Studies of Mycobacterium ulcerans and Buruli Ulcer

Mycolactone A/B is a lipophilic macrocyclic polyketide that is the primary virulence factor produced by Mycobacterium ulcerans, a human pathogen and the causative agent of Buruli ulcer. In M. ulcerans strain Agy99 the mycolactone polyketide synthase (PKS) locus spans a 120 kb region of a 174 kb megaplasmid. Here we have identified promoter regions of this PKS locus using GFP reporter assays, in silico analysis, primer extension, and site-directed mutagenesis. Transcription of the large PKS genes mlsA1 (51 kb), mlsA2 (7 kb) and mlsB (42 kb) is driven by a novel and powerful SigA-like promoter sequence situated 533 bp upstream of both the mlsA1 and mlsB initiation codons, which is also functional in Escherichia coli, Mycobacterium smegmatis and Mycobacterium marinum. Promoter regions were also identified upstream of the putative mycolactone accessory genes mup045 and mup053. We transformed M. ulcerans with a GFP-reporter plasmid under the control of the mls promoter to produce a highly green-fluorescent bacterium. The strain remained virulent, producing both GFP and mycolactone and causing ulcerative disease in mice. Mosquitoes have been proposed as a potential vector of M. ulcerans so we utilized M. ulcerans-GFP in microcosm feeding experiments with captured mosquito larvae. M. ulcerans-GFP accumulated within the mouth and midgut of the insect over four instars, whereas the closely related, non-mycolactone-producing species M. marinum harbouring the same GFP reporter system did not. This is the first report to identify M. ulcerans toxin gene promoters, and we have used our findings to develop M. ulcerans-GFP, a strain in which fluorescence and toxin gene expression are linked, thus providing a tool for studying Buruli ulcer pathogenesis and potential transmission to humans

Sequence-Based Analysis Uncovers an Abundance of Non-Coding RNA in the Total Transcriptome of Mycobacterium tuberculosis

RNA sequencing provides a new perspective on the genome of Mycobacterium tuberculosis by revealing an extensive presence of non-coding RNA, including long 5’ and 3’ untranslated regions, antisense transcripts, and intergenic small RNA (sRNA) molecules. More than a quarter of all sequence reads mapping outside of ribosomal RNA genes represent non-coding RNA, and the density of reads mapping to intergenic regions was more than two-fold higher than that mapping to annotated coding sequences. Selected sRNAs were found at increased abundance in stationary phase cultures and accumulated to remarkably high levels in the lungs of chronically infected mice, indicating a potential contribution to pathogenesis. The ability of tubercle bacilli to adapt to changing environments within the host is critical to their ability to cause disease and to persist during drug treatment; it is likely that novel post-transcriptional regulatory networks will play an important role in these adaptive responses

Optimum DNA curvature using a hybrid approach involving an artificial neural network and genetic algorithm

In the present paper, a hybrid technique involving artificial neural network (ANN) and genetic algorithm (GA) has been proposed for performing modeling and optimization of complex biological systems. In this approach, first an ANN approximates (models) the nonlinear relationship(s) existing between its input and output example data sets. Next, the GA, which is a stochastic optimization technique, searches the input space of the ANN with a view to optimize the ANN output. The efficacy of this formalism has been tested by conducting a case study involving optimization of DNA curvature characterized in terms of the RL value. Using the ANN-GA methodology, a number of sequences possessing high RL values have been obtained and analyzed to verify the existence of features known to be responsible for the occurrence of curvature. A couple of sequences have also been tested experimentally. The experimental results validate qualitatively and also near-quantitatively, the solutions obtained using the hybrid formalism. The ANN-GA technique is a useful tool to obtain, ahead of experimentation, sequences that yield high RL values. The methodology is a general one and can be suitably employed for optimizing any other biological feature

Optimum DNA Curvature Using a Hybrid Approach Involving an Artificial Neural Network and Genetic Algorithm

In the present paper, a hybrid technique involving artificial neural network (ANN) and genetic algorithm (GA) has been proposed for performing modeling and optimization of complex biological systems. In this approach, first an ANN approximates (models) the non-linear relationship(s) existing between its input and output example data sets. Next, the GA, which is a stochastic optimization technique, searches the input space of the ANN with a view to optimize the ANN output. The efficacy of this formalism has been tested by conducting a case study involving optimization of DNA curvature characterized in terms of the RL value. Using the ANN-GA methodology, a number of sequences possessing high RL values have been obtained and analyzed to verify the existence of features known to be responsible for the occurrence of curvature. A couple of sequences have also been tested experimentally. The experimental results validate qualitatively and also near-quantitatively, the solutions obtained using the hybrid formalism. The ANN-GA technique is a useful tool to obtain, ahead of experimentation, sequences that yield high RL values. The methodology is a general one and can be suitably employed for optimizing any other biological feature