mRNA and microRNA analysis reveals modulation of biochemical pathways related to addiction in the ventral tegmental area of methamphetamine self-administering rats

Background Methamphetamine is a highly addictive central nervous system stimulant with increasing levels of abuse worldwide. Alterations to mRNA and miRNA expression within the mesolimbic system can affect addiction-like behaviors and thus play a role in the development of drug addiction. While many studies have investigated the effects of high-dose methamphetamine, and identified neurotoxic effects, few have looked at the role that persistent changes in gene regulation play following methamphetamine self-administration. Therefore, the aim of this study was to identify RNA changes in the ventral tegmental area following methamphetamine self-administration. We performed microarray analyses on RNA extracted from the ventral tegmental area of Sprague–Dawley rats following methamphetamine self-administration training (2 h/day) and 14 days of abstinence. Results We identified 78 miRNA and 150 mRNA transcripts that were differentially expressed (fdr adjusted p < 0.05, absolute log2 fold change >0.5); these included genes not previously associated with addiction (miR-125a-5p, miR-145 and Foxa1), loci encoding receptors related to drug addiction behaviors and genes with previously recognized roles in addiction such as miR-124, miR-181a, DAT and Ret. Conclusion This study provides insight into the effects of methamphetamine on RNA expression in a key brain region associated with addiction, highlighting the possibility that persistent changes in the expression of genes with both known and previously unknown roles in addiction occur.

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.

Unstructured N Terminus of the RNA Polymerase II Subunit Rpb4 Contributes to the Interaction of Rpb4·Rpb7 Subcomplex with the Core RNA Polymerase II of Saccharomyces cerevisiae

Two subunits of eukaryotic RNA polymerase II, Rpb7 and Rpb4, form a subcomplex that has counterparts in RNA polymerases I and III. Although a medium resolution structure has been solved for the 12-subunit RNA polymerase II, the relative contributions of the contact regions between the subcomplex and the core polymerase and the consequences of disrupting them have not been studied in detail. We have identified mutations in the N-terminal ribonucleoprotein-like domain of Saccharomyces cerevisiae Rpb7 that affect its role in certain stress responses, such as growth at high temperature and sporulation. These mutations increase the dependence of Rpb7 on Rpb4 for interaction with the rest of the polymerase. Complementation analysis and RNA polymerase pulldown assays reveal that the Rpb4 center dot Rbp7 subcomplex associates with the rest of the core RNA polymerase II through two crucial interaction points: one at the N-terminal ribonucleoprotein-like domain of Rpb7 and the other at the partially ordered N-terminal region of Rpb4. These findings are in agreement with the crystal structure of the 12-subunit polymerase. We show here that the weak interaction predicted for the N-terminal region of Rpb4 with Rpb2 in the crystal structure actually plays a significant role in interaction of the subcomplex with the core in vivo. Our mutant analysis also suggests that Rpb7 plays an essential role in the cell through its ability to interact with the rest of the polymerase.

Unusual RNA plant virus integration in the soybean genome leads to the production of small RNAs

Horizontal gene transfer (HGT) is known to be a major force in genome evolution. The acquisition of genes from viruses by eukaryotic genomes is a well-studied example of HGT, including rare cases of non-retroviral RNA virus integration. The present study describes the integration of cucumber mosaic virus RNA-1 into soybean genome. After an initial metatranscriptomic analysis of small RNAs derived from soybean, the de novo assembly resulted a 3029-nt contig homologous to RNA-1. The integration of this sequence in the soybean genome was confirmed by DNA deep sequencing. The locus where the integration occurred harbors the full RNA-1 sequence followed by the partial sequence of an endogenous mRNA and another sequence of RNA-1 as an inverted repeat and allowing the formation of a hairpin structure. This region recombined into a retrotransposon located inside an exon of a soybean gene. The nucleotide similarity of the integrated sequence compared to other Cucumber mosaic virus sequences indicates that the integration event occurred recently. We described a rare event of non-retroviral RNA virus integration in soybean that leads to the production of a double-stranded RNA in a similar fashion to virus resistance RNAi plants.

Primer-independent initiation of RNA synthesis by SeMV recombinant RNA-dependent RNA polymerase

Sesbania mosaic virus (SeMV),a single-strand positive-sense RNA plant virus, belongs to the genus Sobemoviruses. Mechanism of replication in Sobemoviruses is poorly understood. In the present study, SeMV RNA-dependent RNA polymerase (RdRp) was overexpressed and purified as a thioredoxin-tagged protein. The recombinant SeMV RdRp could synthesize RNA from genomic or subgenomic RNA templates, even in the absence ofthe protein primer, VPg. Analysis of the product indicated that it was double-stranded and that the mode of initiation was de novo. Mutational analysis of the 3' UTR of subgenomic RNA revealed that a stem-loop structure at the 3' end was important. Further, analysis of this stem-loop showed that the SeMV RdRp was capable of recognizing stem-loop structures of various lengths and forms. These results demonstrate that the SeMV RdRp is capable of primer-independent RNAsynthesis in vitro. (C) 2010 Elsevier Inc. All rights reserved.

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

Transcriptional analysis of persistent Chlamydia pneumoniae infection in vitro

Chlamydia pneumoniae can cause acute respiratory infections including pneumonia. Repeated and persistent Chlamydia infections occur and persistent C. pneumoniae infection may have a role in the pathogenesis of atherosclerosis and coronary heart disease and may also contribute to the development of chronic inflammatory lung diseases like chronic obstructive pulmonary disease (COPD) and asthma. In this thesis in vitro models for persistent C. pneumonia infection were established in epithelial and monocyte/macrophage cell lines. Expression of host cell genes in the persistent C. pneumoniae infection model of epithelial cells was studied by microarray and RT-PCR. In the monocyte/macrophage infection model expression of selected C. pneumoniae genes were studied by RT-PCR and immunofluorescence microscopy. Chlamydia is able to modulate host cell gene expression and apoptosis of host cells, which may assist Chlamydia to evade the host cells' immune responses. This, in turn, may lead to extended survival of the organism inside epithelial cells and promote the development of persistent infection. To simulate persistent C. pneumoniae infection in vivo, we set up a persistent infection model exposing the HL cell cultures to IFN-gamma. When HL cell cultures were treated with moderate concentration of IFN-gamma, the replication of C. pneumoniae DNA was unaffected while differentiation into infectious elementary bodies (EB) was strongly inhibited. By transmission electron microscopy small atypical inclusions were identified in IFN-gamma treated cultures. No second cycle of infection was observed in cells exposed to IFN-gamma , whereas C. pneumoniae was able to undergo a second cycle of infection in unexposed HL cells. Although monocytic cells can naturally restrict chlamydial growth, IFN-gamma further reduced production of infectious C. pneumoniae in Mono Mac 6 cells. Under both studied conditions no second cycle of infection could be detected in monocytic cell line suggesting persistent infection in these cells. As a step toward understanding the role of host genes in the development and pathogenesis of persistent C. pneumoniae infection, modulation of host cell gene expression during IFN-gamma induced persistent infection was examined and compared to that seen during active C. pneumoniae infection or IFN-gamma treatment. Total RNA was collected at 6 to 150 h after infection of an epithelial cell line (HL) and analyzed by a cDNA array (available at that time) representing approximately 4000 human transcripts. In initial analysis 250 of the 4000 genes were identified as differentially expressed upon active and persistent chlamydial infection and IFN-gamma treatment. In persistent infection more potent up-regulation of many genes was observed in IFN-gamma induced persistent infection than in active infection or in IFN-gamma treated cell cultures. Also sustained up-regulation was observed for some genes. In addition, we could identify nine host cell genes whose transcription was specifically altered during the IFN-gamma induced persistent C. pneumoniae infection. Strongest up-regulation in persistent infection in relation to controls was identified for insulin like growth factor binding protein 6, interferon-stimulated protein 15 kDa, cyclin D1 and interleukin 7 receptor. These results suggest that during persistent infection, C. pneumoniae reprograms the host transcriptional machinery regulating a variety of cellular processes including adhesion, cell cycle regulation, growth and inflammatory response, all of which may play important roles in the pathogenesis of persistent C. pneumoniae infection. C. pneumoniae DNA can be detected in peripheral blood mononuclear cells indicating that the bacterium can also infect monocytic cells in vivo and thereby monocytes can assist the spread of infection from the lungs to other anatomical sites. Persistent infection established at these sites could promote inflammation and enhance pathology. Thus, the mononuclear cells are in a strategic position in the development of persistent infection. To investigate the intracellular replication and fate of C. pneumoniae in mononuclear cells we analyzed the transcription of 11 C. pneumoniae genes in Mono Mac 6 cells during infection by real time RT-PCR. Our results suggest that the transcriptional profile of the studied genes in monocytes is different from that seen in epithelial cells and that IFN-gamma has a less significant effect on C. pneumoniae transcription in monocytes. Furthermore, our study shows that type III secretion system (T3SS) related genes are transcribed and that Chlamydia possesses a functional T3SS during infection in monocytes. Since C. pneumoniae infection in monocytes has been implicated to have reduced antibiotic susceptibility, this creates opportunities for novel therapeutics targeting T3SS in the management of chlamydial infection in monocytes.

Characterization of foot-and-mouth disease virus types Ο and Asia 1 RNA

Foot-and-mouth disease is an acute and highly contagious febrile disease affecting cloven-footed animals. Identification of the foot-and-mouth disease virus (FMDV), the causative agent of the disease, posed problems because of the occurrence of many types and subtypes of the virus. A molecular approach based on oligonucleotide mapping of FMDV RNA has been used for the identification and characterization of virus isolates obtained in a disease outbreak (King et al., 1981). One-dimensional oligonucleotide mapping was used for rapid analysis of FMDV RNA (LaTorre et al., 1982). FMDV types Ο and Asia 1 of Indian origin are being routinely used for vaccine production in India. This report presents the differences between FMDV types Ο and Asia 1 at molecular level based on one-dimensional oligonucleotide mapping of virus-induced poly (A) RNA.

Dexamethasone negatively regulates phenobarbitone-activated transcription but synergistically enhances cytoplasmic levels of cytochrome P-450b/e messenger RNA

Dexamethasone has a potentiating effect on phenobarbitone mediated induction of cytochrome P-450b + e mRNAs in adult rat liver. However, the glucocorticoid inhibits phenobarbitone-activated transcription of cytochrome P-450b + e mRNAs by 60-70%. This inhibitory effect is evident in run-off transcription of the endogenous genes as well as in the transcription of an added cloned gene fragment. Dexamethasone inhibits the phenobarbitone-mediated increase in the binding of a transcription factor(s) to the upstream region of the gene as evidenced by gel retardation and Southwestern blot analysis. The glucocorticoid does not stabilize the phenobarbitone-induced polyribosomal cytochrome P-450b + e mRNAs but appears to stabilize the nuclear transcripts. It is proposed that a negative element may mediate the action of dexamethasone at the level of nuclear transcription and stabilization of the nuclear transcript may account for the potentiating effect of the glucocorticoid on phenobarbitone-mediated increase in cytochrome P-450b + e mRNAs in the cytoplasm of the adult rat liver. However, the cytochrome P-450b protein levels are slightly lower in phenobarbitone + dexamethasone treatment than in phenobarbitone-treated liver microsomes.

Analysis of sequence dependent variations in secondary and tertiary structure of tRNA molecules

The double helical regions of the five tRNA(Phe) and two tRNA(Asp) crystal structures have been analyzed using the local basepair step parameters. The sequence dependent effects in the mini double helices of tRNA are very similar to those observed in the crystal structures of oligonucleotides in the A-form, the purine-pyrimidine and purine-purine steps have small roll angles when compared to the fiber models of A-DNA as well as A-RNA, while the pyrimidine-purine doublet steps have large roll angles. The orientation of the basepairs in the D-stem is unusual and invariant i.e. they are different from the other three stems but are very similar in all the five tRNA(Phe) crystal structures, presumably due to tertiary interaction of the Watson-Crick basepairs with other bases, with all bases being highly conserved. The origin of the differences between the tertiary structures of tRNA(Phe) and tRNA(Asp) from yeast has also been investigated. It is found that even though the angle between the acceptor arm and the D-stem is very similar in the two structures, the angle subtended by the acceptor arm and the anticodon arm is smaller in the tRNA(Phe) structure (by more than 10 degrees). This is due to differences in the orientation of the two mini helices constituting the anticodon arm, which are inclined to each other by approximately 25 degrees in tRNA(Phe) and 16 degrees in tRNA(Asp). In addition, the acceptor arm, the D-stem and the anticodon stem are nearly coplanar in tRNA(Phe), while in tRNA(Asp) the anticodon stem projects out of the plane defined by the acceptor arm and the anticodon stem. These two features together lead to a larger separation between the acceptor and anticodon ends in tRNA(Asp) and indicate that the junction between the D-stem and the anticodon stem is quite variable, with features characteristic of a ball-and-socket type joint and determined for each tRNA molecule by the base sequence at the junction.

Purification and characterization of DNA-dependent RNA polymerase II from Candida utilis

DNA-dependent RNA polymerase II from Candida utilis has been purified to near homogeneity. The purified enzyme resolved into three subforms, viz. IIO, IIA and IIB. On SDS-PAGE the enzyme showed ten polypeptides with molecular weights in the range of 205 kDa to 14 kDa. By two dimensional electrophoresis (IEF followed by SDS-PAGE) the presence of basic and acidic polypeptides has been demonstrated. The enzyme showed Km values of 5, 5.6 and 8 mu M for GTP, CTP and ATP, respectively, and the activity was inhibited by low levels of oc-amanitin and antibodies raised against bovine RNA polymerase II. By Western blot analysis the enzyme was found to cross-react with antibodies to bovine RNA polymerase II. RNA polymerase II from G. utilis is a phosphoprotein, the subunits RPB1 and RPB10 were found to be phosphorylated. Analysis of carboxy-terminal domain indicated that it was functionally redundant at least in case of nonspecific transcription, implicating its role in other nuclear processes, such as promoter specific initiation or transcription activation or RNA processing.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

DNA Free Energy-Based Promoter Prediction and Comparative Analysis of Arabidopsis and Rice Genomes

The cis-regulatory regions on DNA serve as binding sites for proteins such as transcription factors and RNA polymerase. The combinatorial interaction of these proteins plays a crucial role in transcription initiation, which is an important point of control in the regulation of gene expression. We present here an analysis of the performance of an in silico method for predicting cis-regulatory regions in the plant genomes of Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) on the basis of free energy of DNA melting. For protein-coding genes, we achieve recall and precision of 96% and 42% for Arabidopsis and 97% and 31% for rice, respectively. For noncoding RNA genes, the program gives recall and precision of 94% and 75% for Arabidopsis and 95% and 90% for rice, respectively. Moreover, 96% of the false-positive predictions were located in noncoding regions of primary transcripts, out of which 20% were found in the first intron alone, indicating possible regulatory roles. The predictions for orthologous genes from the two genomes showed a good correlation with respect to prediction scores and promoter organization. Comparison of our results with an existing program for promoter prediction in plant genomes indicates that our method shows improved prediction capability.

Cryptic AUG is important for 48S ribosomal assembly during internal initiation of translation of coxsackievirus B3 RNA

We have investigated the possible role of a conserved cis-acting element, the cryptic AUG, present in the 5' UTR of coxsackievirus B3 (CVB3) RNA. CVB3 5' UTR contains multiple AUG codons upstream of the initiator AUG, which are not used for the initiation of translation. The 48S ribosomal assembly takes place upstream of the cryptic AUG. We show here that mutation in the cryptic AUG results in reduced efficiency of translation mediated by the CVB3 IRES; mutation also reduces the interaction of mutant IRES with a well characterized IRES trans-acting factor, the human La protein. Furthermore, partial silencing of the La gene showed a decrease in IRES activity in the case of both the wild-type and mutant. We have demonstrated here that the interaction of the 48S ribosomal complex with mutant RNA was weaker compared with wild-type RNA by ribosome assembly analysis. We have also investigated by chemical and enzymic modifications the possible alteration in secondary structure in the mutant RNA. Results suggest that the secondary structure of mutant RNA was only marginally altered. Additionally, we have demonstrated by generating compensatory and non-specific mutations the specific function of the cryptic AUG in internal initiation. Results suggest that the effect of the cryptic AUG is specific and translation could not be rescued. However, a possibility of tertiary interaction of the cryptic AUG with other cis-acting elements cannot be ruled out. Taken together, it appears that the integrity of the cryptic AUG is important for efficient translation initiation by the CVB3 IRES RNA.

Investigations on the RNA binding and phosphorylation of groundnut bud necrosis virus nucleocapsid protein

Groundnut bud necrosis virus belongs to the genus Tospovirus, infects a wide range of crop plants and causes severe losses. To understand the role of the nucleocapsid protein in the viral life cycle, the protein was overexpressed in E. coli and purified by Ni-NTA chromatography. The purified N protein was well folded and was predominantly alpha-helical. Deletion analysis revealed that the C-terminal unfolded region of the N protein was involved in RNA binding. Furthermore, the N protein could be phosphorylated in vitro by Nicotiana benthamiana plant sap and by purified recombinant kinases such as protein kinase CK2 and calcium-dependent protein kinase. This is the first report of phoshphorylation of a nucleocapsid protein in the family Bunyaviridae. The possible implications of the present findings for the viral life cycle are discussed.