Implications of high level pseudogene transcription in Mycobacterium leprae

Background: The Mycobacterium leprae genome has less than 50% coding capacity and 1,133 pseudogenes. Preliminary evidence suggests that some pseudogenes are expressed. Therefore, defining pseudogene transcriptional and translational potentials of this genome should increase our understanding of their impact on M. leprae physiology. Results: Gene expression analysis identified transcripts from 49% of all M. leprae genes including 57% of all ORFs and 43% of all pseudogenes in the genome. Transcribed pseudogenes were randomly distributed throughout the chromosome. Factors resulting in pseudogene transcription included: 1) co-orientation of transcribed pseudogenes with transcribed ORFs within or exclusive of operon-like structures; 2) the paucity of intrinsic stem-loop transcriptional terminators between transcribed ORFs and downstream pseudogenes; and 3) predicted pseudogene promoters. Mechanisms for translational ``silencing'' of pseudogene transcripts included the lack of both translational start codons and strong Shine-Dalgarno (SD) sequences. Transcribed pseudogenes also contained multiple ``in-frame'' stop codons and high Ka/Ks ratios, compared to that of homologs in M. tuberculosis and ORFs in M. leprae. A pseudogene transcript containing an active promoter, strong SD site, a start codon, but containing two in frame stop codons yielded a protein product when expressed in E. coli. Conclusion: Approximately half of M. leprae's transcriptome consists of inactive gene products consuming energy and resources without potential benefit to M. leprae. Presently it is unclear what additional detrimental affect(s) this large number of inactive mRNAs has on the functional capability of this organism. Translation of these pseudogenes may play an important role in overall energy consumption and resultant pathophysiological characteristics of M. leprae. However, this study also demonstrated that multiple translational ``silencing'' mechanisms are present, reducing additional energy and resource expenditure required for protein production from the vast majority of these transcripts.

Involvement of Synthesis and Phosphorylation of Nuclear Protein Factors That Bind to the Positivecis-Acting Element in the Transcriptional Activation of the CYP2B1/B2 Gene by Phenobarbitonein Vivo

The synthesis and phosphorylation of protein factor(s) that bind to the positivecis-acting element (−69 to −98 nt) of the CYP2B1/B2 gene have been examinedin vivoin the rat. Treatment of rats with cycloheximide, a protein synthetic inhibitor, suppresses basal as well as phenobarbitone-induced levels of CYP2B1/B2 mRNA and its run-on transcription. Under these conditions, complex formation of the nuclear extract with the positive element is also inhibited, as judged by gel shift assays. Treatment of rats with 2-aminopurine, a general protein kinase inhibitor, blocks the phenobarbitone-mediated increase in CYP2B1/B2 mRNA, cell-free transcription of a minigene construct containing the positive element, pP450e179DNA, and binding of nuclear proteins to the positive element. Treatment of rats with okadaic acid, a protein phosphatase inhibitor, mimics the effects of phenobarbitone, but only partially. Thus, both phenobarbitone and okadaic acid individually enhance binding of the nuclear protein(s) to the positive element, cell-free transcription of the minigene construct, and phosphorylation of the not, vert, similar26- and 94-kDa proteins binding to the positive element. But unlike phenobarbitone, okadaic acid is not an inducer of CYP2B1/B2 mRNA or its run-on transcription. Thus, phenobarbitone-responsive positive element interactions constitute only a minimal requirement, and okadaic acid is perhaps not able to bring about the total requirement for activation of CYP2B1/B2 gene transcription that should include interaction between the minimal promoter and further upstream elements. An intriguing feature is the antagonistic effect of okadaic acid on phenobarbitone-mediated effects on CYP2B1/B2 mRNA levels, cell-free and run-on transcription, and nuclear protein binding to the positive element. The reason for this antagonism is not clear. It is concluded that phenobarbitone treatment enhancesin vivothe synthesis and phosphorylation of protein factors binding to the positive element and these constitute a minimal requirement for the transcriptional activation of the CYP2B1/B2 gene.

Expression of Human Growth Hormone in Silkworm Larvae through RecombinantBombyx moriNuclear Polyhedrosis Virus

We have generated a recombinantBombyx morinuclear polyhedrosis virus, vBmhGH, harboring the full-length human growth hormone gene (2.4-kb genomic DNA, with four introns and the signal peptide sequences) under the control of the polyhedrin promoter. BmN cells in culture infected with the recombinant virus showed the presence of RNA corresponding to the authentic growth hormone mRNA as well as its incompletly processed precusor. Electrophoretic analysis and immunoprecipitation of proteins of recombinant virus-infected BmN cells revealed the presence of the growth hormone protein. Infection of silkworm larvae with vBmhGH led to the synthesis and efficient secretion of the protein into hemolymph. The recombinant human growth hormone was biologically active in a radioreceptor competition binding assay. The secreted protein was isolated and purified to homogeneity by a single step immunoaffinity chromatography, to a specific activity of 2.4 × 104U/mg. The recombinant hGH retained the immunological and biolological properties of the native peptide. We conclude that BmNPV vectors can be used successfully for expressing chromosomal genes harboring multiple introns.

Synthesis of leader RNA and editing of P mRNA during transcription by rinderpest virus

Purified rinderpest virus was earlier shown to transcribe in vitro, all virus-specific mRNAs with the promoter-proximal N mRNA being the most abundant. Presently, this transcription system has been shown to synthesize full length monocistronic mRNAs comparable to those made in infected cells. Small quantities of bi- and tricistronic mRNAs are also synthesized. Rinderpest virus synthesizes in vitro, a leader RNA of not, vert, similar 55 nucleotides in length. Purified rinderpest virus also exhibits RNA editing activity during the synthesis of P mRNA as shown by primer extension analysis of the mRNA products.

The Mycobacterial MsDps2 Protein Is a Nucleoid-Forming DNA Binding Protein Regulated by Sigma Factors sigma(A) and sigma(B)

The Dps (DNA-binding protein from starved cells) proteins from Mycobacterium smegmatis MsDps1 and MsDps2 are both DNA-binding proteins with some differences. While MsDps1 has two oligomeric states, with one of them responsible for DNA binding, MsDps2 has only one DNA-binding oligomeric state. Both the proteins however, show iron-binding activity. The MsDps1 protein has been shown previously to be induced under conditions of starvation and osmotic stress and is regulated by the extra cellular sigma factors sigma(H) and sigma(F). We show here, that the second Dps homologue in M. smegmatis, namely MsDps2, is purified in a DNA-bound form and exhibits nucleoid-like structures under the atomic force microscope. It appears that the N-terminal sequence of Dps2 plays a role in nucleoid formation. MsDps2, unlike MsDps1, does not show elevated expression in nutritionally starved or stationary phase conditions; rather its promoter is recognized by RNA polymerase containing sigma(A) or sigma(B), under in vitro conditions. We propose that due to the nucleoid-condensing ability, the expression of MsDps2 is tightly regulated inside the cells.

Cyclic AMP in Mycobacteria: Characterization and Functional Role of the Rv1647 Ortholog in Mycobacterium smegmatis{triangledown}

Mycobacterial genomes are endowed with many eukaryote-like nucleotide cyclase genes encoding proteins that can synthesize 3',5'-cyclic AMP (cAMP). However, the roles of cAMP and the need for such redundancy in terms of adenylyl cyclase genes remain unknown. We measured cAMP levels in Mycobacterium smegmatis during growth and under various stress conditions and report the first biochemical and functional characterization of the MSMEG_3780 adenylyl cyclase, whose orthologs in Mycobacterium tuberculosis (Rv1647) and Mycobacterium leprae (ML1399) have been recently characterized in vitro. MSMEG_3780 was important for producing cAMP levels in the logarithmic phase of growth, since the {Delta}MSMEG_3780 strain showed lower intracellular cAMP levels at this stage of growth. cAMP levels decreased in wild-type M. smegmatis under conditions of acid stress but not in the {Delta}MSMEG_3780 strain. This was correlated with a reduction in MSMEG_3780 promoter activity, indicating that the effect of the reduction in cAMP levels on acid stress was caused by a decrease in the transcription of MSMEG_3780. Complementation of the {Delta}MSMEG_3780 strain with the genomic integration of MSMEG_3780 or the Rv1647 gene could restore cAMP levels during logarithmic growth. The Rv1647 promoter was also acid sensitive, emphasizing the biochemical and functional similarities in these two adenylyl cyclases. This study therefore represents the first detailed biochemical and functional analysis of an adenylyl cyclase that is important for maintaining cAMP levels in mycobacteria and underscores the subtle roles that these genes may play in the physiology of the organism.

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.

Evidence for an association with type 2 diabetes mellitus at the PPARG locus in a South Indian population

Peroxisome proliferator activated receptor-gamma 2 (PPARG2) is a nuclear hormone receptor of ligand-dependent ranscription factor involved in adipogenesis and a molecular target of the insulin sensitizers thiazolidinediones. We addressed the question of whether the 3 variants (-1279G/A, Pro12Ala, and His478His) in the PPARG2 gene are associated with type 2 diabetes mellitus and its related traits in a South Indian population. The study subjects (1000 type 2 diabetes mellitus and 1000 normal glucose-tolerant subjects) were chosen randomly from the Chennai Urban Rural Epidemiology Study, an ongoing population-based study in southern India. The variants were screened by single-stranded conformational variant, direct sequencing, and restriction fragment length polymorphism. Linkage disequilibrium was estimated from the estimates of haplotypic frequencies. The -1279G/A, Pro12Ala, and His478His variants of the PPARG2 gene were not associated with type 2 diabetes mellitus. However, the 2-loci analyses showed that, in the presence of Pro/Pro genotype of the Pro12Ala variant, the -1279G/A promoter variant showed increased susceptibility to type 2 diabetes mellitus (odds ratio, 2.092; 95% confidence interval, 1.22-3.59; P = .008), whereas in the presence of 12Ala allele, the -1279G/A showed a protective effect against type 2 diabetes mellitus (odds ratio, 0.270; 95% confidence interval, 0.15-0.49; P < .0001). The 3-loci haplotype analysis showed that the A-Ala-T (-1279G/A-Pro12Ala-His478His) haplotype was associated with a reduced risk of type 2 diabetes mellitus (P < .0001). Although our data indicate that the PPARG2 gene variants, independently, have no association with type 2 diabetes mellitus, the 2-loci genotype analysis involving -1279G/A and Pro12Ala variants and the 3-loci haplotype analysis have shown a significant association with type 2 diabetes mellitus in this South Indian population. (C) 2010 Elsevier Inc. All rights reserved.

Structural and Biochemical Bases for the Redox Sensitivity of Mycobacterium tuberculosis RslA

An effective transcriptional response to redox stimuli is of particular importance for Mycobacterium tuberculosis, as it adapts to the environment of host alveoli and macrophages. The M. tuberculosis a factor sigma(L) regulates the expression of genes involved in cell-wall and polyketide syntheses. sigma(L) interacts with the cytosolic anti-sigma domain of a membrane-associated protein, RslA. Here we demonstrate that RslA binds Zn2+ and can sequester sigma(L) in a reducing environment. In response to an oxidative stimulus, proximal cysteines in the CXXC motif of RslA form a disulfide bond, releasing bound Zn2+. This results in a substantial rearrangement of the sigma(L)/RslA complex, leading to an 8-fold decrease in the affinity of RslA for sigma(L). The crystal structure of the -35-element recognition domain of sigma(L), sigma(L)(4), bound to RslA reveals that RslA inactivates sigma(L) by sterically occluding promoter DNA and RNpolymerase binding sites. The crystal structure further reveals that the cysteine residues that coordinate Zn2+ in RslA are solvent exposed in the complex, thus providing a structural basis for the redox sensitivity of RslA. The biophysical parameters of sigma(L)/RslA interactions provide a template for understanding how variations in the rate of Zn2+ release and associated conformational changes could regulate the activity of a Zn2+-associated anti-sigma factor. (C) 2010 Elsevier Ltd. All rights reserved.

A novel approach for large-scale polypeptide folding based on elastic networks using continuous optimization

We present a new computationally efficient method for large-scale polypeptide folding using coarse-grained elastic networks and gradient-based continuous optimization techniques. The folding is governed by minimization of energy based on Miyazawa–Jernigan contact potentials. Using this method we are able to substantially reduce the computation time on ordinary desktop computers for simulation of polypeptide folding starting from a fully unfolded state. We compare our results with available native state structures from Protein Data Bank (PDB) for a few de-novo proteins and two natural proteins, Ubiquitin and Lysozyme. Based on our simulations we are able to draw the energy landscape for a small de-novo protein, Chignolin. We also use two well known protein structure prediction software, MODELLER and GROMACS to compare our results. In the end, we show how a modification of normal elastic network model can lead to higher accuracy and lower time required for simulation.

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.

Analysis of the beta-glucoside utilization (bgl) genes of Shigella sonnei: evolutionary implications for their maintenance in a cryptic state

The pattern of expression of the genes involved in the utilization of aryl beta-glucosides such as arbutin and salicin is different in the genus Shigella compared to Escherichia coli. The results presented here indicate that the homologue of the cryptic bgl operon of E. coli is conserved in Shigella sonnei and is the primary system involved in beta-glucoside utilization in the organism. The organization of the bgl genes in 5. sonnei is similar to that of E. coli; however there are three major differences in terms of their pattern of expression. (i) The bglB gene, encoding phospho-beta-glucosidase B, is insertionally inactivated in 5. sonnei. As a result, mutational activation of the silent bgl promoter confers an Arbutin-positive (Arb(+)) phenotype to the cells in a single step; however, acquiring a Salicin-positive (Sal(+)) phenotype requires the reversion or suppression of the bglB mutation in addition. (ii) Unlike in E. coli, a majority of the activating mutations (conferring the Arb(+) phenotype) map within the unlinked hns locus, whereas activation of the E. coli bgl operon under the same conditions is predominantly due to insertions within the bglR locus. (iii) Although the bgl promoter is silent in the wild-type strain of 5. sonnei (as in the case of E. coli), transcriptional and functional analyses indicated a higher basal level of transcription of the downstream genes. This was correlated with a 1 bp deletion within the putative Rho-independent terminator present in the leader sequence preceding the homologue of the bglG gene. The possible evolutionary implications of these differences for the maintenance of the genes in the cryptic state are discussed.

Studies on hypomethylation of liver DNA during early stages of chemical carcinogenesis in rat liver

Our finding that the inhibitors of DNA methylation, 5-azacytidine, 5-azadeoxycytidine or adenosine dialdehyde, given after a carcinogen all potentiated initiation suggested that hypomethylation of DNA during repair synthesis of DNA might play a role in the initiation of the carcinogenic process. To examine this aspect further, we have asked the question, do the nodules which develop from initiated cells after promotion with 1% orotic acid exhibit an altered methylation pattern in their DNA? The methylation status of the DNA from nodules has been examined using the restriction endonucleases HpaII/MspI and HhaI which distinguish between methylated and unmethylated cytosines in their nucleotide recognition DNA 5'-CCGG and 5'-GCGC respectively. The proto-oncogenes, c-myc, c-fos and c-Ha-ras, in the DNA were primarily studied in this investigation because of their possible involvement in cell proliferation and/or in cell transformation and tumorigenesis. The results indicate that in the nodule DNA, c-myc and c-fos are hypomethylated in the sequence of CCGG while the c-Ha-ras shows hypomethylation in the alternating GCGC sequence. This methylation pattern seen in the nodule DNA is not found in the DNA of the non-nodular surrounding liver or liver tissue after exposure to promoter or carcinogen alone. It is also not found in the DNA of regenerating liver. It is particularly significant that the methylation patterns in the c-myc and c-Ha-ras regions are similar to those found in several cancer tissues. The results suggest that this methylation pattern is acquired early in the carcinogenic process and raises the question whether it has any bearing on the process.

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/.

Diverse pathways for salicin utilization in Shigella sonnei and Escherichia coli carrying an impaired bgl operon

Utilization of the aryl-beta-glucosides salicin or arbutin in most wild-type strains of E. coli is achieved by a single-step mutational activation of the bgl operon. Shigella sonnei, a branch of the diverse E. coli strain tree, requires two sequential mutational steps for achieving salicin utilization as the bglB gene, encoding the phospho-beta-glucosidase B, harbors an inactivating insertion. We show that in a natural isolate of S. sonnei, transcriptional activation of the gene SSO1595, encoding a phospho-beta-glucosidase, enables salicin utilization with the permease function being provided by the activated bgl operon. SSO1595 is absent in most commensal strains of E. coli, but is present in extra-intestinal pathogens as bgcA, a component of the bgc operon that enables beta-glucoside utilization at low temperature. Salicin utilization in an E. coli bglB laboratory strain also requires a two-step activation process leading to expression of BglF, the PTS-associated permease encoded by the bgl operon and AscB, the phospho-beta-glucosidase B encoded by the silent asc operon. BglF function is needed since AscF is unable to transport beta-glucosides as it lacks the IIA domain involved in phopho-relay. Activation of the asc operon in the Sal(+) mutant is by a promoter-up mutation and the activated operon is subject to induction. The pathway to achieve salicin utilization is therefore diverse in these two evolutionarily related organisms; however, both show cooperation between two silent genetic systems to achieve a new metabolic capability under selection.