Comparison of RNA isolation methods on RNA-Seq: implications for differential expression and meta-analyses.

BACKGROUND: The increasing number of transcriptomic datasets has allowed for meta-analyses, which can be valuable due to their increased statistical power. However, meta-analyses can be confounded by so-called "batch effects," where technical variation across different batches of RNA-seq experiments can clearly produce spurious signals of differential expression and reduce our power to detect true differences. While batch effects can sometimes be accounted for, albeit with caveats, a better strategy is to understand their sources to better avoid them. In this study, we examined the effects of RNA isolation method as a possible source of batch effects in RNA-seq design. RESULTS: Based on the different chemistries of "classic" hot phenol extraction of RNA compared to common commercial RNA isolation kits, we hypothesized that specific mRNAs may be preferentially extracted depending upon method, which could masquerade as differential expression in downstream RNA-seq analyses. We tested this hypothesis using the Saccharomyces cerevisiae heat shock response as a well-validated environmental response. Comparing technical replicates that only differed in RNA isolation method, we found over one thousand transcripts that appeared "differentially" expressed when comparing hot phenol extraction with the two kits. Strikingly, transcripts with higher abundance in the phenol-extracted samples were enriched for membrane proteins, suggesting that indeed the chemistry of hot phenol extraction better solubilizes those species of mRNA. CONCLUSIONS: Within a self-contained experimental batch (e.g. control versus treatment), the method of RNA isolation had little effect on the ability to identify differentially expressed transcripts. However, we suggest that researchers performing meta-analyses across different experimental batches strongly consider the RNA isolation methods for each experiment.

Transcriptional profiling analysis of Penicillium digitatum, the causal agent of citrus green mold, unravels an inhibited ergosterol biosynthesis pathway in response to citral.

BACKGROUND: Green mold caused by Penicillium digitatum is the most damaging postharvest diseases of citrus fruit. Previously, we have observed that citral dose-dependently inhibited the mycelial growth of P. digitatum, with the minimum inhibitory concentration (MIC) of 1.78 mg/mL, but the underlying molecular mechanism is barely understood. RESULTS: In this study, the transcriptional profiling of the control and 1/2MIC-citral treated P. digitatum mycelia after 30 min of exposure were analyzed by RNA-Seq. A total of 6355 genes, including 2322 up-regulated and 4033 down-regulated genes, were found to be responsive to citral. These genes were mapped to 155 KEGG pathways, mainly concerning mRNA surveillance, RNA polymerase, RNA transport, aminoacyl-tRNA biosynthesis, ABC transporter, glycolysis/gluconeogenesis, citrate cycle, oxidative phosphorylation, sulfur metabolism, nitrogen metabolism, inositol phosphate metabolism, fatty acid biosynthesis, unsaturated fatty acids biosynthesis, fatty acid metabolism, and steroid biosynthesis. Particularly, citral exposure affected the expression levels of five ergosterol biosynthetic genes (e.g. ERG7, ERG11, ERG6, ERG3 and ERG5), which corresponds well with the GC-MS results, the reduction in ergosterol content, and accumulation of massive lanosterol. In addition, ERG11, the gene responsible for lanosterol 14α-demethylase, was observed to be the key down-regulated gene in response to citral. CONCLUSION: Our present finding suggests that citral could exhibit its antifungal activity against P. digitatum by the down-regulation of ergosterol biosynthesis.

LbNrt RNA silencing in the mycorrhizal symbiont Laccaria bicolor reveals a nitrate-independent regulatory role for a eukaryotic NRT2-type nitrate transporter.

Fungal nitrogen metabolism plays a fundamental role in function of mycorrhizal symbiosis and consequently in nutrient cycling of terrestrial ecosystems. Despite its global ecological relevance the information on control and molecular regulation of nitrogen utilization in mycorrhizal fungi is very limited. We have extended the nitrate utilization RNA silencing studies of the model mycorrhizal basidiomycete, Laccaria bicolor, by altering the expression of LbNrt, the sole nitrate transporter-encoding gene of the fungus. Here we report the first nutrient transporter mutants for mycorrhizal fungi. Silencing of LbNrt results in fungal strains with minimal detectable LbNrt transcript levels, significantly reduced growth capacity on nitrate and altered symbiotic interaction with poplar. Transporter silencing also creates marked co-downregulation of whole Laccaria fHANT-AC (fungal high-affinity nitrate assimilation cluster). Most importantly, this effect on the nitrate utilization pathway appears independent of extracellular nitrate or nitrogen status of the fungus. Our results indicate a novel and central nitrate uptake-independent regulatory role for a eukaryotic nitrate transporter. The possible cellular mechanisms behind this regulation mode are discussed in the light of current knowledge on NRT2-type nitrate transporters in different eukaryotes.

Small molecule microarrays of RNA-focused peptoids help identify inhibitors of a pathogenic group I intron.

Peptoids that inhibit the group I intron RNA from Candida albicans, an opportunistic pathogen that kills immunocompromised hosts, have been identified using microarrays. The arrayed peptoid library was constructed using submonomers with moieties similar to ones found in small molecules known to bind RNA. Library members that passed quality control analysis were spotted onto a microarray and screened for binding to the C. albicans group I intron ribozyme. Each ligand binder identified from microarray-based screening inhibited self-splicing in the presence of 1 mM nucleotide concentration of bulk yeast tRNA with IC(50)'s between 150 and 2200 microM. The binding signals and the corresponding IC(50)'s were used to identify features in the peptoids that predispose them for RNA binding. After statistical analysis of the peptoids' structures that bind, a second generation of inhibitors was constructed using these important features; all second generation inhibitors have improved potencies with IC(50)'s of <100 microM. The most potent inhibitor is composed of one phenylguanidine and three tryptamine submonomers and has an IC(50) of 31 microM. This compound is 6-fold more potent than pentamidine, a clinically used drug that inhibits self-splicing. These results show that (i) modulators of RNA function can be identified by designing RNA-focused chemical libraries and screening them via microarray; (ii) statistical analysis of ligand binders can identify features in leads that predispose them for binding to their targets; and (iii) features can then be programmed into second generation inhibitors to design ligands with improved potencies.

Stereospecificity of aminoglycoside-ribosomal interactions.

Aminoglycoside antibiotics bind to the A-site decoding region of bacterial rRNA causing mistranslation and/or premature message termination. Aminoglycoside binding to A-site RNA decoding region constructs is established here to be only weakly stereospecific. Mirror-image prokaryotic A-site decoding region constructs were prepared in the natural D-series and the enantiomeric L-series and tested for binding to a series of aminoglycosides. In general, aminoglycosides bind to the D-series decoding region constructs with 2-3-fold higher affinities than they bind to the enantiomeric L-series. Moreover, L-neamine, the enantiomer of naturally occurring D-neamine, was prepared and shown to bind approximately 2-fold more weakly than D-neamine to the natural series decoding region construct, a result consistent with weakly stereospecific binding. The binding of naturally occurring D-neamine and its synthetic L-enantiomer was further evaluated with respect to binding to prokaryotic and eukaryotic ribosomes. Here, weak stereospecifcity was again observed with L-neamine being the more potent binder by a factor of approximately 2. However, on a functional level, unnatural L-neamine proved to inhibit in vitro translation with significantly lower potency (approximately 5-fold) than D-neamine. In addition, both L- and D-neamine are bacteriocidal toward Gram-(-) bacteria. L-Neamine inhibits the growth of E. coli and P. aeruginosa with 8- and 3-fold higher MIC than D-neamine. Interestingly, L-neamine also inhibits the growth of aminoglycoside-resistant E. coli, which expresses a kinase able to phosphorylate and detoxify aminoglycosides of the D-series. These observations suggest that mirror-image aminoglycosides may avoid certain forms of enzyme-mediated resistance.

SKIP is an indispensable factor for Caenorhabditis elegans development.

SKI-binding protein (SKIP) is a transcription cofactor present in all eukaryotes. Here we show that SKIP is a unique protein that is required for Caenorhabditis elegans viability and development. Expression of CeSKIP (skp-1) assayed by RT-PCR and by GFP fluorescence in transgenic lines starts in embryos and continues to adulthood. Loss of CeSKIP activity by RNA-mediated inhibition results in early embryonic arrest similar to that seen following inhibition of RNA polymerase II. RNA polymerase II phosphorylation appears normal early in CeSKIP RNA-mediated inhibition treated embryos although the expression of several embryonic GFP reporter genes is severely restricted or absent. Our data suggest that CeSKIP is an essential component of many RNA polymerase II transcription complexes and is indispensable for C. elegans development.

Anticodon domain methylated nucleosides of yeast tRNA(Phe) are significant recognition determinants in the binding of a phage display selected peptide.

The contributions of the natural modified nucleosides to RNA identity in protein/RNA interactions are not understood. We had demonstrated that 15 amino acid long peptides could be selected from a random phage display library using the criterion of binding to a modified, rather than unmodified, anticodon domain of yeast tRNA(Phe) (ASL(Phe)). Affinity and specificity of the selected peptides for the modified ASL(Phe) have been characterized by fluorescence spectroscopy of the peptides' tryptophans. One of the peptides selected, peptide t(F)2, exhibited the highest specificity and most significant affinity for ASL(Phe) modified with 2'-O-methylated cytidine-32 and guanosine-34 (Cm(32) and Gm(34)) and 5-methylated cytidine-40 (m(5)C(40)) (K(d) = 1.3 +/- 0.4 microM) and a doubly modified ASL(Phe)-Gm(34),m(5)C(40) and native yeast tRNA(Phe) (K(d) congruent with 2.3 and 3.8 microM, respectively) in comparison to that for the unmodified ASL(Phe) (K(d) = 70.1 +/- 12.3 microM). Affinity was reduced when a modification altered the ASL loop structure, and binding was negated by modifications that disfavored hairpin formation. Peptide t(F)2's higher affinity for the ASL(Phe)-Cm(32),Gm(34),m(5)C(40) hairpin and fluorescence resonance energy transfer from its tryptophan to the hypermodified wybutosine-37 in the native tRNA(Phe) placed the peptide across the anticodon loop and onto the 3'-side of the stem. Inhibition of purified yeast phenylalanyl-tRNA synthetase (FRS) catalyzed aminoacylation of cognate yeast tRNA(Phe) corroborated the peptide's binding to the anticodon domain. The phage-selected peptide t(F)2 has three of the four amino acids crucial to G(34) recognition by the beta-structure of the anticodon-binding domain of Thermus thermophilus FRS and exhibited circular dichroism spectral properties characteristic of beta-structure. Thus, modifications as simple as methylations contribute identity elements that a selected peptide specifically recognizes in binding synthetic and native tRNA and in inhibiting tRNA aminoacylation.

Binding enhancement by tertiary interactions and suicide inhibition of a Candida albicans group I intron by phosphoramidate and 2'-O-methyl hexanucleotides.

Candida albicans is one of many infectious pathogens that are evolving resistance to current treatments. RNAs provide a large class of targets for new therapeutics for fighting these organisms. One strategy for targeting RNAs uses short oligonucleotides that exhibit binding enhancement by tertiary interactions in addition to Watson-Crick pairing. A potential RNA target in C. albicans is the self-splicing group I intron in the LSU rRNA precursor. The recognition elements that align the 5' exon splice site for a ribozyme derived from this precursor are complex [Disney, M. D., Haidaris, C. G., and Turner, D. H. (2001) Biochemistry 40, 6507-6519]. These recognition elements have been used to guide design of hexanucleotide mimics of the 5' exon that have backbones modified for nuclease stability. These hexanucleotides bind as much as 100000-fold more tightly to a ribozyme derived from the intron than to a hexanucleotide mimic of the intron's internal guide sequence, r(GGAGGC). Several of these oligonucleotides inhibit precursor self-splicing via a suicide inhibition mechanism. The most promising suicide inhibitor is the ribophosphoramidate rn(GCCUC)rU, which forms more trans-spliced than cis-spliced product at oligonucleotide concentrations of >100 nM at 1 mM Mg(2+). The results indicate that short oligonucleotides modified for nuclease stability can target catalytic RNAs when the elements of tertiary interactions are complex.

Discovery of selective, small-molecule inhibitors of RNA complexes--II. Self-splicing group I intron ribozyme.

Self-splicing group I intron RNA was chosen as a potential therapeutic target for small-molecule intervention. High-throughput screening methodologies have been developed to identify small organic molecules that regulate the activities of these catalytic introns. Group introns derived from pathogenic Pneumocystis carinii and phage T4 were used as model systems. Inhibitors identified from a library of approximately equal to 150,000 compounds were shown to regulate biochemical reactions including the two-step intron splicing and an RNA ligation catalyzed by the group I introns. These inhibitors provide a unique opportunity to understand small-molecule recognition of the self-splicing RNA. The methodologies developed for group I introns should be applicable to studies of other RNA systems.

[Effect of cycloheximide and actinomycin D on the phospholipid and neutral lipid content of Cunninghamella japonica]. / Vliianie tsiklogeksimida i aktinomitsina D na soderzhanie fosfolipidov i neitral'nykh lipidov Cunninghamella japonica.

Changes in the composition of neutral lipids and phospholipids were studied during the growth of Cunninghamella japonica under the action of cycloheximide and actinomycin D. The data were used to discuss how the synthesis of enzymes catalysing the formation of individual lipid classes was regulated and whether it would be possible to control the composition of phospholipids and neutral lipids using inhibitors of RNA and protein synthesis. The results are also indicative of a certain correlation between growth phases of the fungus and changes in certain characteristics of membrane lipids.