Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition

The rapid increase in genome sequence information has necessitated the annotation of their functional elements, particularly those occurring in the non-coding regions, in the genomic context. Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. In this analysis, we demonstrate that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. We have therefore obtained a set of free energy threshold values, for genomic DNA with varying GC content and used them as generic criteria for predicting promoter regions in several microbial genomes, using an in-house developed tool `PromPredict'. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (50.8% GC) and B. subtilis (43.5% GC) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained.

Kinetic characterization of rat brain type IIA sodium channel alpha-subunit stably expressed in a somatic cell line

1. The rat brain type IIA Na+ channel alpha-subunit was stably expressed in Chinese hamster ovary (CHO) cells. Current through the expressed Na+ channels was studied using the whole-cell configuration of the patch clamp technique. The transient Na+ current was sensitive to TTX and showed a bell-shaped peak current vs. membrane potential relation. 2. Na+ current inactivation was better described by the sum of two exponentials in the potential range -30 to +40 mV, with. a dominating fast component and a small slower component. 3. The steady-state inactivation, h(infinity), was related to potential by a Boltzmann distribution, underlying thr ee states of the inactivation gate. 4. Recovery of the channels from inactivation at different potentials in the range -70 to -120 mV were characterized by al? initial delay which decreased with hyperpolarization. The time course was well fitted by the sum of two exponentials. In this case the slower exponential was the major component, and both time constants decreased with hyperpolarization. 5. For a working description of the Na+ channel inactivation in this preparation, with a minimal deviation from the Hodgkin-Huxley model, a three-state scheme of the form O reversible arrow I-1 reversible arrow I-2 was proposed, replacing the original two-state scheme of the Hodgkin-Huxley model, and the rate constants are reported. 6. The instantaneous current-voltage relationship showed marked deviation from linearity and was satisfactorily fitted by the constant-field equation. 7. The time course of activation was described by an m(x) model. However, the best-fitted value of x varied with the membrane potential and had a mean value of 2. 8. Effective gating charge was determined to be 4.7e from the slope of the activation plot, plotted on a logarithmic scale. 9. The rate constants of activation, alpha(m) and beta(m), were determined. Their functional dependence on the membrane potential was investigated.

CsrA interacting small RNAs in Haemophilus spp genomes: a theoretical analysis

The csrA is a carbon storage regulator gene that encodes a protein with multiple RNA interaction sites. Bacterial non-coding small RNAs like csrB, csrC and their counterparts in diverse bacterial genus are identified to control the regulatory activities of CsrA and its orthologs. An attempt has been made in this study to identify 'novel' non-coding small RNAs that are involved in the regulatory activities of csrA gene. All CsrA-interacting small RNAs are computationally fingerprinted to have multiple occurrence of 7-nucleotide CsrA interacting repeats [CAGGA(U/A/C)G] along with a 18-nucleotide upstream binding site. However, in several of the genomes like Haemophilus spp, the upstream binding site is not identified. The current methodology overcomes this difficulty by identifying small RNA-specific orphan transcriptional units within the intergenic regions of the genome. The results could identify all known CsrA-interacting small RNAs in E. coli, Vibrio cholerae and Pseudomonas aeruginosa genomes and additionally has picked six new possible CsrA-interacting small RNA regions in E. coli. Our computational analysis indicates that known rygD and rprA sRNAs in E. coli could possibly interact with CsrA proteins. The study is extended to three of the Haemophilus genomes that could identify seven new possible CsrA interacting small RNAs.

Transcription of Human Resistin Gene Involves an Interaction of Sp1 with Peroxisome Proliferator-Activating Receptor Gamma (PPAR gamma)

Background: Resistin is a cysteine rich protein, mainly expressed and secreted by circulating human mononuclear cells. While several factors responsible for transcription of mouse resistin gene have been identified, not much is known about the factors responsible for the differential expression of human resistin.Methodology/Principal Finding: We show that the minimal promoter of human resistin lies within similar to 80 bp sequence upstream of the transcriptional start site (-240) whereas binding sites for cRel, CCAAT enhancer binding protein alpha (C/EBP-alpha), activating transcription factor 2 (ATF-2) and activator protein 1 (AP-1) transcription factors, important for induced expression, are present within sequences up to -619. Specificity Protein 1(Sp1) binding site (-276 to -295) is also present and an interaction of Sp1 with peroxisome proliferator activating receptor gamma (PPAR gamma) is necessary for constitutive expression in U937 cells. Indeed co-immunoprecipitation assay demonstrated a direct physical interaction of Sp1 with PPAR gamma in whole cell extracts of U937 cells. Phorbol myristate acetate (PMA) upregulated the expression of resistin mRNA in U937 cells by increasing the recruitment of Sp1, ATF-2 and PPAR gamma on the resistin gene promoter. Furthermore, PMA stimulation of U937 cells resulted in the disruption of Sp1 and PPAR gamma interaction. Chromatin immunoprecipitation (ChIP) assay confirmed the recruitment of transcription factors phospho ATF-2, Sp1, Sp3, PPAR gamma, chromatin modifier histone deacetylase 1 (HDAC1) and the acetylated form of histone H3 but not cRel, C/EBP-alpha and phospho c-Jun during resistingene transcription.Conclusion: Our findings suggest a complex interplay of Sp1 and PPAR gamma along with other transcription factors that drives the expression of resistin in human monocytic U937 cells.

Methylenetetrahydrofolate reductase gene polymorphisms and risk of acute lymphoblastic leukemia in children

Methylenetetrahydrofolate reductase (MTHFR) is a critical enzyme in folate metabolism and is involved in DNA synthesis, DNA repair and DNA methylation. Genetic polymorphisms of this enzyme have been shown to impact several diseases, including cancer. Leukemias are malignancies arising from rapidly proliferating hematopoietic cells having great requirement of DNA synthesis. This case-control study was undertaken to analyze the association of the MTHFR gene polymorphisms 677 C"T and 1298 A"C and the risk of acute lymphoblastic leukemia in children. Materials and Methods: Eighty-six patients aged below 15 years with a confirmed diagnosis of acute lymphoblastic leukemia (ALL) and 99 matched controls were taken for this study. Analysis of the polymorphisms was done using the polymerase chain reaction -restriction fragment length polymorphism (PCR-RFLP) method. Results: Frequency of MTHFR 677 CC and CT were 85.9% and 14.1% in the controls, and 84.9% and 15.1% in the cases. The 'T' allele frequency was 7% and 7.5% in cases and controls respectively. The frequency of MTHFR 1298 AA, AC, and CC were 28.3%, 55.6% and 16.1% for controls and 23.3%, 59.3% and 17.4% for cases respectively. The 'C' allele frequency for 1298 A→C was 43.9% and 47% respectively for controls and cases. The odds ratio (OR) for C677T was 1.08 (95% CI 0.48- 2.45, p = 0.851) and OR for A1298C was 1.29(95% CI 0.65-2.29, p = 0.46) and OR for 1298 CC was 1.31 (95% CI 0.53-3.26, p =0.56). The OR for the combined heterozygous status (677 CT and 1298 AC) was 1.94 (95% CI 0.58 -6.52, p = 0.286). Conclusion: The prevalence of 'T' allele for 677 MTHFR polymorphism was low in the population studied. There was no association between MTHFR 677 C→T and 1298 A→C gene polymorphisms and risk of ALL, which may be due to the small sample size.

Trishanku, a novel regulator of cell-type stability and morphogenesis in Dictyostelium discoideum

We have identified a novel gene, trishanku (triA), by random insertional mutagenesis of Dictyostelium discoideum. TriA is a Broad complex Tramtrack bric-a-brac domain-containing protein that is expressed strongly during the late G2 phase of cell cycle and in presumptive spore (prespore (psp)) cells. Disrupting triA destabilizes cell fate and reduces aggregate size; the fruiting body has a thick stalk, a lowered spore: stalk ratio, a sub-terminal spore mass and small, rounded spores. These changes revert when the wild-type triA gene is re-expressed under a constitutive or a psp-specific promoter. By using short- and long-lived reporter proteins, we show that in triA(-) slugs the prestalk (pst)/psp proportion is normal, but that there is inappropriate transdifferentiation between the two cell types. During culmination, regardless of their current fate, all cells with a history of pst gene expression contribute to the stalk, which could account for the altered cell-type proportion in the mutant.

Comparison of promoter-specific events during transcription initiation in mycobacteria

DNA protein interactions that occur during transcription initiation play an important role in regulating gene expression. To initiate transcription, RNA polymerase (RNAP) binds to promoters in a sequence-specific fashion. This is followed by a series of steps governed by the equilibrium binding and kinetic rate constants, which in turn determine the overall efficiency of the transcription process. We present here the first detailed kinetic analysis of promoter RNAP interactions during transcription initiation in the sigma(A)-dependent promoters P-rrnAPCL1, P-rrnB and P-gyr of Mycobacterium smegmatis. The promoters show comparable equilibrium binding affinity but differ significantly in open complex formation, kinetics of isomerization and promoter clearance. Furthermore, the two rrn promoters exhibit varied kinetic properties during transcription initiation and appear to be subjected to different modes of regulation. In addition to distinct kinetic patterns, each one of the housekeeping promoters studied has its own rate-limiting step in the initiation pathway, indicating the differences in their regulation.

sRNAscanner: A Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes

Background:Bacterial non-coding small RNAs (sRNAs) have attracted considerable attention due to their ubiquitous nature and contribution to numerous cellular processes including survival, adaptation and pathogenesis. Existing computational approaches for identifying bacterial sRNAs demonstrate varying levels of success and there remains considerable room for improvement. Methodology/Principal Findings: Here we have proposed a transcriptional signal-based computational method to identify intergenic sRNA transcriptional units (TUs) in completely sequenced bacterial genomes. Our sRNAscanner tool uses position weight matrices derived from experimentally defined E. coli K-12 MG1655 sRNA promoter and rho-independent terminator signals to identify intergenic sRNA TUs through sliding window based genome scans. Analysis of genomes representative of twelve species suggested that sRNAscanner demonstrated equivalent sensitivity to sRNAPredict2, the best performing bioinformatics tool available presently. However, each algorithm yielded substantial numbers of known and uncharacterized hits that were unique to one or the other tool only. sRNAscanner identified 118 novel putative intergenic sRNA genes in Salmonella enterica Typhimurium LT2, none of which were flagged by sRNAPredict2. Candidate sRNA locations were compared with available deep sequencing libraries derived from Hfq-co-immunoprecipitated RNA purified from a second Typhimurium strain (Sittka et al. (2008) PLoS Genetics 4: e1000163). Sixteen potential novel sRNAs computationally predicted and detected in deep sequencing libraries were selected for experimental validation by Northern analysis using total RNA isolated from bacteria grown under eleven different growth conditions. RNA bands of expected sizes were detected in Northern blots for six of the examined candidates. Furthermore, the 5'-ends of these six Northern-supported sRNA candidates were successfully mapped using 5'-RACE analysis. Conclusions/Significance: We have developed, computationally examined and experimentally validated the sRNAscanner algorithm. Data derived from this study has successfully identified six novel S. Typhimurium sRNA genes. In addition, the computational specificity analysis we have undertaken suggests that similar to 40% of sRNAscanner hits with high cumulative sum of scores represent genuine, undiscovered sRNA genes. Collectively, these data strongly support the utility of sRNAscanner and offer a glimpse of its potential to reveal large numbers of sRNA genes that have to date defied identification. sRNAscanner is available from: http://bicmku.in:8081/sRNAscanner or http://cluster.physics.iisc.ernet.in/sRNAscanner/.

DNA Clasping by Mycobacterial HU: The C-Terminal Region of HupB Mediates Increased Specificity of DNA Binding

Background: HU a small, basic, histone like protein is a major component of the bacterial nucleoid. E. coli has two subunits of HU coded by hupA and hupB genes whereas Mycobacterium tuberculosis (Mtb) has only one subunit of HU coded by ORF Rv2986c (hupB gene). One noticeable feature regarding Mtb HupB, based on sequence alignment of HU orthologs from different bacteria, was that HupB(Mtb) bears at its C-terminal end, a highly basic extension and this prompted an examination of its role in Mtb HupB function. Methodology/Principal Findings: With this objective two clones of Mtb HupB were generated; one expressing full length HupB protein (HupB(Mtb)) and another which expresses only the N terminal region (first 95 amino acid) of hupB (HupB(MtbN)). Gel retardation assays revealed that HupBMtbN is almost like E. coli HU (heat stable nucleoid protein) in terms of its DNA binding, with a binding constant (K-d) for linear dsDNA greater than 1000 nM, a value comparable to that obtained for the HU alpha alpha and HU alpha beta forms. However CTR (C-terminal Region) of HupB(Mtb) imparts greater specificity in DNA binding. HupB(Mtb) protein binds more strongly to supercoiled plasmid DNA than to linear DNA, also this binding is very stable as it provides DNase I protection even up to 5 minutes. Similar results were obtained when the abilities of both proteins to mediate protection against DNA strand cleavage by hydroxyl radicals generated by the Fenton's reaction, were compared. It was also observed that both the proteins have DNA binding preference for A: T rich DNA which may occur at the regulatory regions of ORFs and the oriC region of Mtb. Conclusions/Significance: These data thus point that HupB(Mtb) may participate in chromosome organization in-vivo, it may also play a passive, possibly an architectural role.

Identifying N60D mutation in omega subunit of Escherichia coli RNA polymerase by bottom-up proteomic approach

Escherichia coli RNA polymerase is a multi-subunit enzyme containing alpha(2)beta beta'omega sigma, which transcribes DNA template to intermediate RNA product in a sequence specific manner. Although most of the subunits are essential for its function, the smallest subunit omega (average molecular mass similar to 10,105 Da) can be deleted without affecting bacterial growth. Creating a mutant of the omega subunit can aid in improving the understanding of its role. Sequencing of rpoZ gene that codes for omega subunit from a mutant variant suggested a substitution mutation at position 60 of the protein: asparagine (N) -> aspartic acid (D). This mutation was verified at the protein level by following a typical mass spectrometry (MS) based bottom-up proteomic approach. Characterization of in-gel trypsin digested samples by reverse phase liquid chromatography (LC) coupled to electrospray ionization (ESI)-tandem mass spectrometry (MS/MS) enabled in ascertaining this mutation. Electron transfer dissociation (ETD) of triply charged (M + 3H)(3+)] tryptic peptides (residues 53-67]), EIEEGLINNQILDVR from wild-type and EIEEGLIDNQILDVR from mutant, facilitated in unambiguously determining the site of mutation at residue 60.

Functional Identification of Regulatory Elements within the Promoter Region of Platelet-Derived Growth Factor 2

Human platelet-derived growth factor (PDGF) is composed of two polypeptide chains, PDGF-1 and PDGF-2,the human homolog of the v-sis oncogene. Deregulation of PDGF-2 expression can confer a growth advantage to cells possessing the cognate receptor and, thus, may contribute to the malignant phenotype. We investigated the regulation of PDGF-2 mRNA expression during megakaryocytic differentiation of K562 cells. Induction by 12-O-tetradecanoylphorbol-13-acetate (TPA) led to a greater than 200-fold increase in PDGF-2 transcript levels in these cells. Induction was dependent on protein synthesis and was not enhanced by cycloheximide exposure.In our initial investigation of the PDGF-2 promoter, a minimal promoter region, which included sequences extending only 42 base pairs upstream of the TATA signal, was found to be as efficient as 4 kilobase pairs upstream of the TATA signal in driving expression of a reporter gene in uninduced K562 cells. We also functionally identified different regulatory sequence elements of the PDGF-2 promoter in TPA-induced K562 cells. One region acted as a transcriptional silencer, while another region was necessary for maximal activity of the promoter in megakaryoblasts. This region was shown to bind nuclear factors and was the target for trans-activation in normal and tumor cells. In one tumor cell line, which expressed high PDGF-2 mRNA levels, the presence of the positive regulatory region resulted in a 30-fold increase in promoter activity. However, the ability of the minimal PDGF-2 promoter to drive reporter gene expression in uninduced K562 cells and normal fibroblasts, which contained no detectable PDGF-2 transcripts, implies the existence of other negative control mechanisms beyond the regulation of promoter activity.

High-quality annotation of promoter regions for 913 bacterial genomes

Motivation: The number of bacterial genomes being sequenced is increasing very rapidly and hence, it is crucial to have procedures for rapid and reliable annotation of their functional elements such as promoter regions, which control the expression of each gene or each transcription unit of the genome. The present work addresses this requirement and presents a generic method applicable across organisms. Results: Relative stability of the DNA double helical sequences has been used to discriminate promoter regions from non-promoter regions. Based on the difference in stability between neighboring regions, an algorithm has been implemented to predict promoter regions on a large scale over 913 microbial genome sequences. The average free energy values for the promoter regions as well as their downstream regions are found to differ, depending on their GC content. Threshold values to identify promoter regions have been derived using sequences flanking a subset of translation start sites from all microbial genomes and then used to predict promoters over the complete genome sequences. An average recall value of 72% (which indicates the percentage of protein and RNA coding genes with predicted promoter regions assigned to them) and precision of 56% is achieved over the 913 microbial genome dataset.

Specificity in the regulation of eukaryotic gene transcription

The regulation of eukaryotic gene transcription poses major challenges in terms of the innumerable protein factors required to ensure tissue or cell-type specificity. While this specificity is sought to be explained by the interaction of cis-acting DNA elements and thetrans-acting protein factor(s), considerable amount of degeneracy has been observed in this interaction. Immunoglobulin heavy chain gene expression in B cells and liver-specific gene expression are discussed as examples of this complexity in this article. Heterodimerization and post-translational modification of transcription factors and the organization of composite promoter elements are strategies by which diverse sets of genes can be regulated in a specific manner using a finite number of protein factors

High level stable expression of rat brain type IIA sodium channel ?-subunit in CHO cells

The neuronal sodium channels are responsible for the rising phase of action potential and are composed of three subunits, of which the alpha-subunit has been shown to be adequate for most of its functional properties. We have stably expressed the rat brain type IIA sodium channel alpha-subunit in CHO cell tine using a CMV promoter-based vector. The expression was confirmed by detecting a 6.5 kb RNA corresponding to sodium channel alpha-subunit using Northern hybridization. The cells stably expressing the alpha-subunit, yield isolated sodium currents of amplitudes greater than 4nA when studied in whole-cell configuration of the patch-clamp technique. The sodium currents are characterized by activation and inactivation properties similar to neuronal sodium channels, and are blocked by the voltage gated sodium channel blocker tetrodotoxin (TTX).

Cardiac remodeling in Drosophila arises from changes in actin gene expression and from a contribution of lymph gland-like cells to the heart musculature

Understanding the basis of normal heart remodeling can provide insight into the plasticity of the cardiac state, and into the potential for treating diseased tissue. In Drosophila, the adult heart arises during metamorphosis from a series of events, that include the remodeling of an existing cardiac tube, the elaboration of new inflow tracts, and the addition of a layer of longitudinal muscle fibers. We have identified genes active in all these three processes, and studied their expression in order to characterize in greater detail normal cardiac remodeling. Using a Transglutaminase-lacZ transgenic line, that is expressed in the inflow tracts of the larval and adult heart, we confirm the existence of five inflow tracts in the adult structure. In addition, expression of the Actin87E actin gene is initiated in the remodeling cardiac tube, but not in the longitudinal fibers, and we have identified an Act87E promoter fragment that recapitulates this switch in expression. We also establish that the longitudinal fibers are multinucleated, characterizing these cells as specialized skeletal muscles. Furthermore, we have defined the origin of the longitudinal fibers, as a subset of lymph gland cells associated with the larval dorsal vessel. These studies underline the myriad contributors to the formation of the adult Drosophila heart, and provide new molecular insights into the development of this complex organ. (C) 2011 Elsevier Ireland Ltd. All rights reserved.