CENTROIDFOLD: a web server for RNA secondary structure prediction

MicroRNA target prediction: theory and practice

The present study is one of the few that includes tissue samples in the evaluation of target prediction algorithms designed to detect microRNA (miRNA) sequences that might interact with particular messenger RNA (mRNA) sequences. Twelve different target prediction tools were used to find miRNA sequences that might interact with CCL20 gene expression. Different algorithms predicted controversial miRNA sequences for CCL20 regulation due to a different weighting of parameters. Hsa-miR-21 and hsa-miR-145 suggested by four or more programs were chosen for further investigation. Possible real interaction of these miRNA sequences with CCL20 gene expression was monitored using luciferase assays and expression analyses of tissue samples of colorectal adenocarcinoma by either qRT-PCR or ELISA. Folding status of seed-binding sites in complete mRNA and 3'UTR of CCL20 was predicted. Prediction of miRNA expression was attempted based on CCL20 expression data. Eight of the target prediction tools forecasted a role for hsa-miR-21 and four mentioned hsa-miR-145 in CCL20 gene regulation. Laboratory experimentation showed that CCL20 may serve as a target of hsa-miR-21 but not hsa-miR-145. Expression of the molecules resulted in no clear assertion. Folding of seed-binding sites was predicted to be relatively constant for the complete mRNA and 3'UTR. Predicting miRNA expression based on target gene expression was impossible. This might be attributable to the fact that effects of miRNA activity may oscillate between gene product repression and activation. Additional systematic studies are needed to address this issue.

Noncoding RNA in Mycobacteria

Efforts to understand the molecular basis of mycobacterial gene regulation are dominated by a protein-centric view. However, there is a growing appreciation that noncoding RNA, i.e., RNA that is not translated, plays a role in a wide variety of molecular mechanisms. Noncoding RNA comprises rRNA, tRNA, 4.5S RNA, RnpB, and transfer-messenger RNA, as well as a vast population of regulatory RNA, often dubbed "the dark matter of gene regulation." The regulatory RNA species comprise 5' and 3' untranslated regions and a rapidly expanding category of transcripts with the ability to base-pair with mRNAs or to interact with proteins. Regulatory RNA plays a central role in the bacterium's response to changes in the environment, and in this article we review emerging information on the presence and abundance of different types of noncoding RNA in mycobacteria.

Long non-coding RNA and alternative splicing modulations in Parkinson's leukocytes identified by RNA sequencing

The continuously prolonged human lifespan is accompanied by increase in neurodegenerative diseases incidence, calling for the development of inexpensive blood-based diagnostics. Analyzing blood cell transcripts by RNA-Seq is a robust means to identify novel biomarkers that rapidly becomes a commonplace. However, there is lack of tools to discover novel exons, junctions and splicing events and to precisely and sensitively assess differential splicing through RNA-Seq data analysis and across RNA-Seq platforms. Here, we present a new and comprehensive computational workflow for whole-transcriptome RNA-Seq analysis, using an updated version of the software AltAnalyze, to identify both known and novel high-confidence alternative splicing events, and to integrate them with both protein-domains and microRNA binding annotations. We applied the novel workflow on RNA-Seq data from Parkinson's disease (PD) patients' leukocytes pre- and post- Deep Brain Stimulation (DBS) treatment and compared to healthy controls. Disease-mediated changes included decreased usage of alternative promoters and N-termini, 5′-end variations and mutually-exclusive exons. The PD regulated FUS and HNRNP A/B included prion-like domains regulated regions. We also present here a workflow to identify and analyze long non-coding RNAs (lncRNAs) via RNA-Seq data. We identified reduced lncRNA expression and selective PD-induced changes in 13 of over 6,000 detected leukocyte lncRNAs, four of which were inversely altered post-DBS. These included the U1 spliceosomal lncRNA and RP11-462G22.1, each entailing sequence complementarity to numerous microRNAs. Analysis of RNA-Seq from PD and unaffected controls brains revealed over 7,000 brain-expressed lncRNAs, of which 3,495 were co-expressed in the leukocytes including U1, which showed both leukocyte and brain increases. Furthermore, qRT-PCR validations confirmed these co-increases in PD leukocytes and two brain regions, the amygdala and substantia-nigra, compared to controls. This novel workflow allows deep multi-level inspection of RNA-Seq datasets and provides a comprehensive new resource for understanding disease transcriptome modifications in PD and other neurodegenerative diseases.

Long non-coding RNAs (lncRNAs) comprise a novel, fascinating class of RNAs with largely unknown biological functions. Parkinson's-disease (PD) is the most frequent motor disorder, and Deep-brain-stimulation (DBS) treatment alleviates the symptoms, but early disease biomarkers are still unknown and new future genetic interference targets are urgently needed. Using RNA-sequencing technology and a novel computational workflow for in-depth exploration of whole-transcriptome RNA-seq datasets, we detected and analyzed lncRNAs in sequenced libraries from PD patients' leukocytes pre and post-treatment and the brain, adding this full profile resource of over 7,000 lncRNAs to the few human tissues-derived lncRNA datasets that are currently available. Our study includes sample-specific database construction, detecting disease-derived changes in known and novel lncRNAs, exons and junctions and predicting corresponding changes in Polyadenylation choices, protein domains and miRNA binding sites. We report widespread transcript structure variations at the splice junction and exons levels, including novel exons and junctions and alteration of lncRNAs followed by experimental validation in PD leukocytes and two PD brain regions compared with controls. Our results suggest lncRNAs involvement in neurodegenerative diseases, and specifically PD. This comprehensive workflow will be of use to the increasing number of laboratories producing RNA-Seq data in a wide range of biomedical studies.

Getting RIDD of RNA: IRE1 in cell fate regulation

Inositol-requiring enzyme 1 (IRE1) is the most conserved transducer of the unfolded protein response (UPR), a homeostatic response that preserves proteostasis. Intriguingly, via its endoribonuclease activity, IRE1 produces either adaptive or death signals. This occurs through both unconventional splicing of XBP1 mRNA and regulated IRE1-dependent decay of mRNA (RIDD). Whereas XBP1 mRNA splicing is cytoprotective in response to endoplasmic reticulum (ER) stress, RIDD has revealed many unexpected features. For instance, RIDD cleaves RNA at an XBP1-like consensus site but with an activity divergent from XBP1 mRNA splicing and can either preserve ER homeostasis or induce cell death. Here we review recent findings on RIDD and propose a model of how IRE1 RNase activity might control cell fate decisions.

Diversification of 13S globulins, allergenic seed storage proteins, of common buckwheat

The α polypeptide of the 13S globulin subunit of common buckwheat is the counterpart of the major allergenic β polypeptide. Trypsin digestibility varies between variants of the α polypeptide with and without a tandem repeat insert. To evaluate the intra-species diversity of 13S globulin, the comprehensive screening of a genomic DNA library was performed, resulting in the isolation of 14 and 3 genes for Met-poor and Met-rich subunits, respectively. Although most tandem repeat units were 45 bp in length, the two-repeat gene Glb2B and all one-repeat genes contained an additional 3 bp. Secondary structure predictions and polyacrylamide gel electrophoresis demonstrated that the sense strand of Glb2B-CCG, the additional 3 bp-deletion clone of Glb2B, formed a more rigid secondary structure than that of the wild-type. Thus, the large intra-species variation of 13S globulin revealed in this study and its diversification might be attributable to the unique nature of the tandem repeat sequences.

Mutation of FVS1, encoding a protein with a sterile alpha motif domain, affects asexual reproduction in the fungal plant pathogen Fusarium oxysporum

Fusarium oxysporum produces three kinds of asexual spores: microconidia, macroconidia and chlamydospores. We previously analysed expressed sequence tags during vegetative growth and conidiation in F. oxysporum and found 42 genes that were markedly upregulated during conidiation compared to vegetative growth. One of the genes, FVS1, encodes a protein with a sterile alpha motif (SAM) domain, which functions in protein-protein interactions that are involved in transcriptional or post-transcriptional regulation and signal transduction. Here, we made FVS1-disrupted mutants from the melon wilt pathogen F. oxysporum f. sp. melonis. Although the mutants produced all three kinds of asexual spores with normal morphology, they formed markedly fewer microconidia and macroconidia than the wild type. The mutants appeared to have a defect in the development of the conidiogenesis cells, conidiophores and phialides, required for the formation of microconidia and macroconidia. In contrast, chlamydospore formation was dramatically promoted in the mutants. The growth rates of the mutants on media were slightly reduced, indicating that FVS1 is also involved in, but not essential for, vegetative growth. We also observed that mutation of FVS1 caused defects in conidial germination and virulence, suggesting that the Fvs1 has pleiotropic functions in F. oxysporum.

ISRNA: an integrative online toolkit for short reads from high-throughput sequencing data

Integrative Short Reads NAvigator (ISRNA) is an online toolkit for analyzing high-throughput small RNA sequencing data. Besides the high-speed genome mapping function, ISRNA provides statistics for genomic location, length distribution and nucleotide composition bias analysis of sequence reads. Number of reads mapped to known microRNAs and other classes of short non-coding RNAs, coverage of short reads on genes, expression abundance of sequence reads as well as some other analysis functions are also supported. The versatile search functions enable users to select sequence reads according to their sub-sequences, expression abundance, genomic location, relationship to genes, etc. A specialized genome browser is integrated to visualize the genomic distribution of short reads. ISRNA also supports management and comparison among multiple datasets.

AVAILABILITY

ISRNA is implemented in Java/C++/Perl/MySQL and can be freely accessed at http://omicslab.genetics.ac.cn/ISRNA/.

Aptamers as theranostic agents: modifications, serum stability and functionalisation

Aptamers, and the selection process known as Systematic Evolution of Ligands by Exponential Enrichment (SELEX) used to generate them, were first described more than twenty years ago. Since then, there have been numerous modifications to the selection procedures. This review discusses the use of modified bases as a means of enhancing serum stability and producing effective therapeutic tools, as well as functionalising these nucleic acids to be used as potential diagnostic agents.

Novel role for molecular transporter importin 9 in posttranscriptional regulation of IFN-ε expression

IFN-ε is a unique type I IFN whose constitutive expression in lung, brain, small intestine, and reproductive tissues is only partially understood. Our previous observation that posttranscriptional events participate in the regulation of IFN-ε mRNA expression led us to investigate whether the 5' and/or 3' untranslated regions (UTR) have regulatory functions. Surprisingly, we found that full-length IFN-ε 5'UTR markedly suppressed mRNA expression under basal conditions. Analysis of the secondary structure of this region predicted formation of two stable stem-loop structures, loops 1 and 2. Studies using luciferase constructs harboring various stretches of IFN-ε 5'UTR and mutant constructs in which the conformation of loop structures was disrupted showed that loop 1 is essential for regulation of mRNA expression. Incubation of HeLa cell extracts with agarose-bound RNAs harboring IFN-ε loop structures identified importin 9 (IPO9), a molecular transporter and chaperone, as a candidate that associates with these regions of the 5'UTR. IPO9 overexpression decreased, and IPO9 silencing increased basal IFN-ε expression. Our studies uncover a previously undescribed function for IPO9 as a specific, and negative, posttranscriptional regulator of IFN-ε expression, and they identify key roles for IFN-ε stem-loop structure 1 in this process. IPO9-mediated effects on 5'UTRs appear to extend to additional mRNAs, including hypoxia-inducible factor-1α, that can form specific loop structures.

Comprehensive investigation of microRNAs enhanced by analysis of sequence variants, expression patterns, ARGONAUTE loading, and target cleavage

MicroRNAs (miRNAs) are a class of small RNAs that typically function by guiding the cleavage of target messenger RNAs. They have been shown to play major roles in a variety of plant processes, including development, and responses to pathogens and environmental stresses. To identify new miRNAs and regulation in Arabidopsis (Arabidopsis thaliana), 27 small RNA libraries were constructed and sequenced from various tissues, stresses, and small RNA biogenesis mutants, resulting in 95 million genome-matched sequences. The use of rdr2 to enrich the miRNA population greatly enhanced this analysis and led to the discovery of new miRNAs arising from both known and new precursors, increasing the total number of Arabidopsis miRNAs by about 10%. Parallel Analysis of RNA Ends data provide evidence that the majority guide target cleavage. Many libraries represented novel stress/tissue conditions, such as submergence-stressed flowers, which enabled the identification of new stress regulation of both miRNAs and their targets, all of which were validated in wild-type plants. By combining small RNA expression analysis with ARGONAUTE immunoprecipitation data and global target cleavage data from Parallel Analysis of RNA Ends, a much more complete picture of Arabidopsis miRNAs was obtained. In particular, the discovery of ARGONAUTE loading and target cleavage biases gave important insights into tissue-specific expression patterns, pathogen responses, and the role of sequence variation among closely related miRNA family members that would not be evident without this combinatorial approach.

Fighting against uncertainty: an essential issue in bioinformatics

Many bioinformatics problems, such as sequence alignment, gene prediction, phylogenetic tree estimation and RNA secondary structure prediction, are often affected by the 'uncertainty' of a solution, that is, the probability of the solution is extremely small. This situation arises for estimation problems on high-dimensional discrete spaces in which the number of possible discrete solutions is immense. In the analysis of biological data or the development of prediction algorithms, this uncertainty should be handled carefully and appropriately. In this review, I will explain several methods to combat this uncertainty, presenting a number of examples in bioinformatics. The methods include (i) avoiding point estimation, (ii) maximum expected accuracy (MEA) estimations and (iii) several strategies to design a pipeline involving several prediction methods. I believe that the basic concepts and ideas described in this review will be generally useful for estimation problems in various areas of bioinformatics.

Directed evolution study unveiling key sequence factors that affect translation efficiency in Escherichia coli

Synonymous mutations in protein coding genes significantly impact translation efficiency. We synthesized a pair of genes encoding green fluorescent protein that were separated by 160 synonymous mutations to investigate key factors that affect translation efficiency. One sequence was optimized for Escherichia coli (GFP(Eco)) and the other for Bacillus subtilis (GFP(Bsu)). When the genes were expressed in E. coli, GFP(Eco) fluoresced 12-fold stronger than GFP(Bsu), confirming the suboptimal nature of the GFP(Bsu) gene. We then employed directed evolution to improve the expression of GFP(Bsu). Random mutagenesis and DNA shuffling was used to generate mutant libraries, which were screened for fluorescence. A variant showing 6-fold fluorescence enhancement was identified, which contained a single mutation (G10A) in a rare codon for Gly-4. However, the substitution generated another type of rare codon, AGA, for Arg, suggesting that the improvement was caused by a factor other than the rare codon. We next applied saturation mutagenesis to Gly-4. The darkest variant contained a GGG codon (GFP(Bsu)-G) for Gly-4. Taking the location of the mutation into account, we hypothesized that destabilization of the mRNA secondary structure around the initiation codon improved the expression. We then randomized the nucleotide triplet in 5'-untranslated region (5'UTR) of GFP(Bsu), which is complementary to the Gly-4 codon. A variant showing 6-fold fluorescence enhancement was identified, which exhibited a destabilized secondary structure. When this 5'UTR sequence was combined with GFP(Bsu)-G, 22-fold fluorescent improvement was achieved. Collectively, the stability of the mRNA secondary structure around the initiation codon predominantly affected the translation efficiency.

Structural requirement in Clostridium perfringens collagenase mRNA 5' leader sequence for translational induction through small RNA-mRNA base pairing

The Gram-positive anaerobic bacterium Clostridium perfringens is pathogenic to humans and animals, and the production of its toxins is strictly regulated during the exponential phase. We recently found that the 5' leader sequence of the colA transcript encoding collagenase, which is a major toxin of this organism, is processed and stabilized in the presence of the small RNA VR-RNA. The primary colA 5'-untranslated region (5'UTR) forms a long stem-loop structure containing an internal bulge and masks its own ribosomal binding site. Here we found that VR-RNA directly regulates colA expression through base pairing with colA mRNA in vivo. However, when the internal bulge structure was closed by point mutations in colA mRNA, translation ceased despite the presence of VR-RNA. In addition, a mutation disrupting the colA stem-loop structure induced mRNA processing and ColA-FLAG translational activation in the absence of VR-RNA, indicating that the stem-loop and internal bulge structure of the colA 5' leader sequence is important for regulation by VR-RNA. On the other hand, processing was required for maximal ColA expression but was not essential for VR-RNA-dependent colA regulation. Finally, colA processing and translational activation were induced at a high temperature without VR-RNA. These results suggest that inhibition of the colA 5' leader structure through base pairing is the primary role of VR-RNA in colA regulation and that the colA 5' leader structure is a possible thermosensor.

Stress-induced lncRNAs evade nuclear degradation and enter the translational machinery

Long noncoding RNAs (lncRNAs) play important roles in the regulation of gene expression. In fission yeast, glucose starvation triggers a transcriptional cascade of polyadenylated lncRNAs in the upstream region of the fructose-1,6-bisphosphatase gene (fbp1(+) ), which is correlated with stepwise chromatin remodeling and necessary for the massive induction of fbp1(+) mRNA. Here, we show that these novel metabolic stress-induced lncRNAs (mlonRNAs) are 5'-capped, less stable than fbp1(+) mRNA and sensitive to a certain extent to the nuclear exosome cofactor Rrp6. However, most mlonRNAs seem to escape nuclear degradation and are exported to the cytoplasm, where they localize to polysomes precisely during glucose starvation-induced global translation inhibition. It is likely that ribosomes tend to accumulate in the upstream region of mlonRNAs. Although mlonRNAs contain an unusual amount of upstream AUGs (uAUGs) and small open reading frames (uORFs), they escape Upf1-mediated targeting to the non-sense-mediated decay (NMD) pathway. The deletion of Upf1 had no effect on mlonRNA stability, but considerably destabilized fbp1(+) mRNA, hinting toward a possible novel role of Upf1. Our findings suggest that the stability of mlonRNAs is distinctly regulated from mRNA and previously described noncoding transcripts.

CompaRNA: a server for continuous benchmarking of automated methods for RNA secondary structure prediction

We present a continuous benchmarking approach for the assessment of RNA secondary structure prediction methods implemented in the CompaRNA web server. As of 3 October 2012, the performance of 28 single-sequence and 13 comparative methods has been evaluated on RNA sequences/structures released weekly by the Protein Data Bank. We also provide a static benchmark generated on RNA 2D structures derived from the RNAstrand database. Benchmarks on both data sets offer insight into the relative performance of RNA secondary structure prediction methods on RNAs of different size and with respect to different types of structure. According to our tests, on the average, the most accurate predictions obtained by a comparative approach are generated by CentroidAlifold, MXScarna, RNAalifold and TurboFold. On the average, the most accurate predictions obtained by single-sequence analyses are generated by CentroidFold, ContextFold and IPknot. The best comparative methods typically outperform the best single-sequence methods if an alignment of homologous RNA sequences is available. This article presents the results of our benchmarks as of 3 October 2012, whereas the rankings presented online are continuously updated. We will gladly include new prediction methods and new measures of accuracy in the new editions of CompaRNA benchmarks.

Admixture mapping in lupus identifies multiple functional variants within IFIH1 associated with apoptosis, inflammation, and autoantibody production

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22–24 (LOD = 6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ∼1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [P_meta = 5.20×10⁻¹⁴; odds ratio, 95% confidence interval = 0.82 (0.78–0.87)], and two missense variants, rs1990760 (Ala946Thr) [P_meta = 3.08×10⁻⁷; 0.88 (0.84–0.93)] and rs10930046 (Arg460His) [P_dom = 1.16×10⁻⁸; 0.70 (0.62–0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

African-Americans (AA) are at increased risk of systemic lupus erythematosus (SLE), but the genetic basis of this risk increase is largely unknown. We used admixture mapping to localize disease-causing genetic variants that differ in frequency across populations. This approach is advantageous for localizing susceptibility genes in recently admixed populations like AA. Our genome-wide admixture scan identified seven admixture signals, and we followed the best signal at 2q22–24 with fine-mapping, imputation-based association analysis and experimental validation. We identified two independent coding variants and a non-coding variant within the IFIH1 gene associated with SLE. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Suppression of microRNAs by dual-targeting and clustered Tough Decoy inhibitors

MicroRNAs (miRNAs) are ubiquitous regulators of gene expression that contribute to almost any cellular process. Methods for managing of miRNA activity are attracting increasing attention in relation to diverse experimental and therapeutic applications. DNA-encoded miRNA inhibitors expressed from plasmid or virus-based vectors provide persistent miRNA suppression and options of tissue-directed micromanaging. In this report, we explore the potential of exploiting short, hairpin-shaped RNAs for simultaneous suppression of two or more miRNAs. Based on the "Tough Decoy" (TuD) design, we create dual-targeting hairpins carrying two miRNA recognition sites and demonstrate potent co-suppression of different pairs of unrelated miRNAs by a single DNA-encoded inhibitor RNA. In addition, enhanced miRNA suppression is achieved by expression of RNA polymerase II-transcribed inhibitors carrying clustered TuD hairpins with up to a total of eight miRNA recognition sites. Notably, by expressing clustered TuD inhibitors harboring a single recognition site for each of a total of six miRNAs, we document robust parallel suppression of multiple miRNAs by inhibitor RNA molecules encoded by a single expression cassette. These findings unveil a new potential of TuD-based miRNA inhibitors and pave the way for standardizing synchronized suppression of families or clusters of miRNAs.

Identification of small RNAs in Mycobacterium smegmatis using heterologous Hfq

Gene regulation by small RNAs (sRNAs) has been extensively studied in various bacteria. However, the presence and roles of sRNAs in mycobacteria remain largely unclear. Immunoprecipitation of RNA chaperone Hfq to enrich for sRNAs is one of the effective methods to isolate sRNAs. However, the lack of an identified mycobacterial hfq restricts the feasibility of this approach. We developed a novel method that takes advantage of the conserved inherent sRNAs-binding capability of heterologous Hfq from Escherichia coli to enrich sRNAs from Mycobacterium smegmatis, a model organism for studying Mycobacterium tuberculosis. We validated 12 trans-encoded and 12 cis-encoded novel sRNAs in M. smegmatis. Many of these sRNAs are differentially expressed at exponential phase compared with stationary phase, suggesting that sRNAs are involved in the growth of mycobacteria. Intriguingly, five of the cis-encoded novel sRNAs target known transposases. Phylogenetic conservation analysis shows that these sRNAs are pathogenicity dependent. We believe that our findings will serve as an important reference for future analysis of sRNAs regulation in mycobacteria and will contribute significantly to the development of sRNAs prediction programs. Moreover, this novel method of using heterologous Hfq for sRNAs enrichment can be of general use for the discovery of bacterial sRNAs in which no endogenous Hfq is identified.

In silico experiment with an-antigen-toll like receptor-5 agonist fusion construct for immunogenic application to Helicobacter pylori

BACKGROUNDS

Helicobacter pylori colonize the gastric mucosa of half of the world's population. Although it is classified as a definitive type I carcinogen by World Health Organization, there is no effective vaccine against this bacterium. H. pylori evade the host immune response by avoiding toll-like detection, such as detection via toll-like receptor-5 (TLR-5). Thus, a chimeric construct consisting of selected epitopes from virulence factors that is incorporated into a TLR-5 ligand (Pseudomonas flagellin) could result in more potent innate and adaptive immune responses.

MATERIALS AND METHODS

Based on the histocompatibility antigens of BALB/c mice, in silico techniques were used to select several fragments from H. pylori virulence factors with a high density of B- and T-cell epitopes.

RESULTS

These segments consist of cytotoxin-associated geneA (residue 162-283), neutrophil activating protein (residue 30-135) and outer inflammatory protein A (residue 155-268). The secondary and tertiary structure of the chimeric constructs and other bioinformatics analyses such as stability, solubility, and antigenicity were performed. The chimeric construct containing antigenic segments of H. pylori proteins was fused with the D3 domain of Pseudomonas flagellin. This recombinant chimeric gene was optimized for expression in Escherichia coli. The in silico results showed that the conserved C- and N-terminal domains of flagellin and the antigenicity of selected fragments were retained.

DISCUSSION

In silico analysis showed that Pseudomonas flagellin is a suitable platform for incorporation of an antigenic construct from H. pylori. This strategy may be an effective tool for the control of H. pylori and other persistent infections.

Class I/II hybrid inhibitory oligodeoxynucleotide exerts Th1 and Th2 double immunosuppression

We designed class I/II hybrid inhibitory oligodeoxynucleotides (iODNs), called iSG, and found that the sequence 5′-TTAGGG-3′, which has a six-base loop head structure, and a 3′-oligo (dG)_3–5 tail sequence are important for potent immunosuppressive activity. Interestingly, splenocytes isolated from ovalbumin (OVA)-immunized mice and treated with iSG3 showed suppression of not only interleukin (IL)-6, IL-12p35, IL-12p40, and interferon (IFN) γ mRNA expression, but also IL-4 and IL-13 mRNA expression. Thus, both Th2 and Th1 immune responses can be strongly suppressed by iODNs in splenocytes from allergen-immunized mice, suggesting usefulness in the treatment of diseases induced by over-active immune activation.

Fibroblast growth factor (FGF) signaling during gastrulation negatively modulates the abundance of microRNAs that regulate proteins required for cell migration and embryo patterning

FGF signaling plays a pivotal role in regulating cell movements and lineage induction during gastrulation. Here we identify 44 microRNAs that are expressed in the primitive streak region of gastrula stage chicken embryos. We show that the primary effect of FGF signaling on microRNA abundance is to negatively regulate the levels of miR-let-7b, -9, -19b, -107, -130b, and -218. LIN28B inhibits microRNA processing and is positively regulated by FGF signaling. Gain- and loss-of-function experiments show that LIN28B negatively regulates the expression of miR-19b, -130b, and let-7b, whereas negative modulation of miR-9, -107, and -218 appears to be independent of LIN28B function. Predicted mRNA targets of the FGF-regulated microRNAs are over-represented in serine/threonine and tyrosine kinase receptors, including ACVR1, ACVR2B, PDGFRA, TGFBR1, and TGFBR3. Luciferase assays show that these and other candidates are targeted by FGF-regulated microRNAs. PDGFRA, a receptor whose activity is required for cell migration through the primitive streak, is a target of miR-130b and -218 in vivo. These results identify a novel mechanism by which FGF signaling regulates gene expression by negatively modulating microRNA abundance through both LIN28B-dependent and LIN28B-independent pathways.

Evolution of an inhibitory RNA aptamer against T7 RNA polymerase

Aptamers are promising gene components that can be used for the construction of synthetic gene circuits. In this study, we isolated an RNA aptamer that specifically inhibits transcription of T7 RNA polymerase (RNAP). The 38-nucleotide aptamer, which was a shortened variant of an initial SELEX isolate, showed moderate inhibitory activity. By stepwise doped-SELEX, we isolated evolved variants with strong inhibitory activity. A 29-nucleotide variant of a doped-SELEX isolate showed 50% inhibitory concentration at 11 nM under typical in vitro transcription conditions. Pull-down experiments revealed that the aptamer inhibited the association of T7 RNAP with T7 promoter DNA.

Rtips: fast and accurate tools for RNA 2D structure prediction using integer programming

We present a web-based tool set Rtips for fast and accurate prediction of RNA 2D complex structures. Rtips comprises two computational tools based on integer programming, IPknot for predicting RNA secondary structures with pseudoknots and RactIP for predicting RNA–RNA interactions with kissing hairpins. Both servers can run much faster than existing services with the same purpose on large data sets as well as being at least comparable in prediction accuracy. The Rtips web server along with the stand-alone programs is freely accessible at http://rna.naist.jp/.

HEXIM1-binding elements on mRNAs identified through transcriptomic SELEX and computational screening

The positive transcription elongation factor b (P-TEFb) is one of the main regulatory factors of the transcription mediated by RNA polymerase II (RNAPII). P-TEFb promotes transcriptional elongation by phosphorylating its targets, which include the C-terminal domain of RNAPII. The activity of P-TEFb is negatively regulated by an RNA-binding protein HEXIM1 in association with 7SK snRNA. To search for other cellular RNAs that bind to HEXIM1, we used systematic evolution of ligands by exponential enrichment (SELEX) with the HeLa cDNA library as the initial pool source. We identified cad mRNA as a HEXIM1-binding RNA and confirmed their association in HeLa cells. In vitro mutational analysis showed that cad mRNA binds to HEXIM1 through its bulged stem structure located in exon 11. In addition, a computational search revealed other RNAs with similar stem structures, including brd4 and tcf3 mRNAs, both of which were shown to be coimmunoprecipitable with anti-HEXIM1 antibody in HeLa cells. Our findings suggest a possible role for HEXIM1 in the regulation of specific gene expressions.

RNA-PAIRS: RNA probabilistic assignment of imino resonance shifts

The significant biological role of RNA has further highlighted the need for improving the accuracy, efficiency and the reach of methods for investigating RNA structure and function. Nuclear magnetic resonance (NMR) spectroscopy is vital to furthering the goals of RNA structural biology because of its distinctive capabilities. However, the dispersion pattern in the NMR spectra of RNA makes automated resonance assignment, a key step in NMR investigation of biomolecules, remarkably challenging. Herein we present RNA Probabilistic Assignment of Imino Resonance Shifts (RNA-PAIRS), a method for the automated assignment of RNA imino resonances with synchronized verification and correction of predicted secondary structure. RNA-PAIRS represents an advance in modeling the assignment paradigm because it seeds the probabilistic network for assignment with experimental NMR data, and predicted RNA secondary structure, simultaneously and from the start. Subsequently, RNA-PAIRS sets in motion a dynamic network that reverberates between predictions and experimental evidence in order to reconcile and rectify resonance assignments and secondary structure information. The procedure is halted when assignments and base-parings are deemed to be most consistent with observed crosspeaks. The current implementation of RNA-PAIRS uses an initial peak list derived from proton-nitrogen heteronuclear multiple quantum correlation ((1)H-(15)N 2D HMQC) and proton-proton nuclear Overhauser enhancement spectroscopy ((1)H-(1)H 2D NOESY) experiments. We have evaluated the performance of RNA-PAIRS by using it to analyze NMR datasets from 26 previously studied RNAs, including a 111-nucleotide complex. For moderately sized RNA molecules, and over a range of comparatively complex structural motifs, the average assignment accuracy exceeds 90%, while the average base pair prediction accuracy exceeded 93%. RNA-PAIRS yielded accurate assignments and base pairings consistent with imino resonances for a majority of the NMR resonances, even when the initial predictions are only modestly accurate. RNA-PAIRS is available as a public web-server at http://pine.nmrfam.wisc.edu/RNA/.

Identification of cis- and trans-acting factors involved in the localization of MALAT-1 noncoding RNA to nuclear speckles

MALAT-1 noncoding RNA is localized to nuclear speckles despite its mRNA-like characteristics. Here, we report the identification of several key factors that promote the localization of MALAT-1 to nuclear speckles and also provide evidence that MALAT-1 is involved in the regulation of gene expression. Heterokaryon assays revealed that MALAT-1 does not shuttle between the nucleus and cytoplasm. RNAi-mediated repression of the nuclear speckle proteins, RNPS1, SRm160, or IBP160, which are well-known mRNA processing factors, resulted in the diffusion of MALAT-1 to the nucleoplasm. We demonstrated that MALAT-1 contains two distinct elements directing transcripts to nuclear speckles, which were also capable of binding to RNPS1 in vitro. Depletion of MALAT-1 represses the expression of several genes. Taken together, our results suggest that RNPS1, SRm160, and IBP160 contribute to the localization of MALAT-1 to nuclear speckles, where MALAT-1 could be involved in regulating gene expression.

Improved prediction of RNA tertiary structure with insights into native state dynamics

The importance of RNA tertiary structure is evident from the growing number of published high resolution NMR and X-ray crystallographic structures of RNA molecules. These structures provide insights into function and create a knowledge base that is leveraged by programs such as Assemble, ModeRNA, RNABuilder, NAST, FARNA, Mc-Sym, RNA2D3D, and iFoldRNA for tertiary structure prediction and design. While these methods sample native-like RNA structures during simulations, all struggle to capture the native RNA conformation after scoring. We propose RSIM, an improved RNA fragment assembly method that preserves RNA global secondary structure while sampling conformations. This approach enhances the quality of predicted RNA tertiary structure, provides insights into the native state dynamics, and generates a powerful visualization of the RNA conformational space. RSIM is available for download from http://www.github.com/jpbida/rsim.

A proposal for a new HIV-1 DLS structural model

The dimer initiation site/dimer linkage sequence (DIS/DLS) region of the human immunodeficiency virus type 1 (HIV-1) RNA genome is suggested to play essential roles at various stages of the viral life cycle. Through a novel assay we had recently developed, we reported on the necessary and sufficient region for RNA dimerization in the HIV-1 virion. Using this system, we performed further detailed mapping of the functional base pairs necessary for HIV-1 DLS structure. Interestingly, the study revealed a previously unnoticed stem formation between two distantly positioned regions. Based on this and other findings on functional base pairing in vivo, we propose new 3D models of the HIV-1 DLS which contain a unique pseudoknot-like conformation. Since this pseudoknot-like conformation appears to be thermodynamically stable, forms a foundational skeleton for the DLS and sterically restricts the spontaneous diversification of DLS conformations, its unique shape may contribute to the viral life cycle and potentially serve as a novel target for anti-HIV-1 therapies.

A potent 2'-O-methylated RNA-based microRNA inhibitor with unique secondary structures

MicroRNAs (miRNAs) are involved in various biological processes and human diseases. The development of strong low-molecular weight inhibitors of specific miRNAs is thus expected to be useful in providing tools for basic research or in generating promising new therapeutic drugs. We have previously described the development of ‘Tough Decoy (TuD) RNA’ molecules, which achieve the long-term suppression of specific miRNA activity in mammalian cells when expressed from a lentivirus vector. In our current study, we describe new synthetic miRNA inhibitors, designated as S-TuD (Synthetic TuD), which are composed of two fully 2′-O-methylated RNA strands. Each of these strands includes a miRNA-binding site. Following the hybridization of paired strands, the resultant S-TuD forms a secondary structure with two stems, which resembles the corresponding TuD RNA molecule. By analyzing the effects of S-TuD against miR-21, miR-200c, miR-16 and miR-106b, we have elucidated the critical design features of S-TuD molecules that will provide optimum inhibitory effects following transfection into human cell lines. We further show that the inhibitory effects of a single transfection of S-TuD-miR200c are quite long-lasting (>7 days) and induce partial EMT, the full establishment of which requires 11 days when using a lentivirus vector that expresses TuD-miR200c continuously.

Genetic and structural analysis of HIV-1 Rev responsive element related to V38A and T18A enfuvirtide resistance mutations

BACKGROUND

For the expression of late viral genes, HIV-1 efficiently exploits the nuclear export by using Rev viral protein, which specifically binds the RNA Rev Responsive Element (RRE). This region is contained within the gp120-gp41 encoding sequence. Enfuvirtide is the first approved HIV-1 fusion-inhibitor, and gp41 codons associated with primary enfuvirtide-resistance (amino-acids 36-45) are localized within the RRE structure. We previously found the co-presence of V38A+T18A resistance mutations in patients failing enfuvirtide.

METHODS

Collecting 476 and 135 HIV-1 B-subtype gp41 sequences from enfuvirtide-naïve and enfuvirtide-treated patients, respectively, two mutations previously found associated with enfuvirtide treatment, T18A and V38A, were analyzed. Moreover, the RNA secondary structure was displayed by CONTRAfold-software and the gp41 evolutionary pathways by a mutagenetic tree.

RESULTS

By modeling the RRE structure, we show that the T18 and V38 codons are base pairing within the RRE-stem-IIA, an important domain involved in Rev binding. While a structural RRE impairment in the presence of V38A alone was found, a restoration of the original RRE structure occurred in co-presence of V38A+T18A. By mutagenetic tree analysis, a compensatory evolution confirming our hypothesis on the structural modification mechanism was observed.

CONCLUSION

We show that enfuvirtide pressure may also affect specific RRE domains involved in Rev binding, thus requiring a compensatory evolution able to preserve the secondary structure of the RRE.

Massive analysis of rice small RNAs: mechanistic implications of regulated microRNAs and variants for differential target RNA cleavage

Small RNAs have a variety of important roles in plant development, stress responses, and other processes. They exert their influence by guiding mRNA cleavage, translational repression, and chromatin modification. To identify previously unknown rice (Oryza sativa) microRNAs (miRNAs) and those regulated by environmental stress, 62 small RNA libraries were constructed from rice plants and used for deep sequencing with Illumina technology. The libraries represent several tissues from control plants and plants subjected to different environmental stress treatments. More than 94 million genome-matched reads were obtained, resulting in more than 16 million distinct small RNA sequences. This allowed an evaluation of ~400 annotated miRNAs with current criteria and the finding that among these, ~150 had small interfering RNA-like characteristics. Seventy-six new miRNAs were found, and miRNAs regulated in response to water stress, nutrient stress, or temperature stress were identified. Among the new examples of miRNA regulation were members of the same miRNA family that were differentially regulated in different organs and had distinct sequences Some of these distinct family members result in differential target cleavage and provide new insight about how an agriculturally important rice phenotype could be regulated in the panicle. This high-resolution analysis of rice miRNAs should be relevant to plant miRNAs in general, particularly in the Poaceae.

A distinct class of internal ribosomal entry site in members of the Kobuvirus and proposed Salivirus and Paraturdivirus genera of the Picornaviridae

The 5'-untranslated regions (5' UTRs) of picornavirus genomes contain an internal ribosomal entry site (IRES) that promotes the end-independent initiation of translation. Picornavirus IRESs are classified into four structurally distinct groups, each with different initiation factor requirements. Here, we identify a fifth IRES class in members of Kobuvirus, Salivirus, and Paraturdivirus genera of Picornaviridae: Aichi virus (AV), bovine kobuvirus (BKV), canine kobuvirus (CKoV), mouse kobuvirus (MKoV), sheep kobuvirus (SKV), salivirus A (SV-A), turdivirus 2 (TV2), and TV3. The 410-nucleotide (nt)-long AV IRES comprises four domains (I to L), including a hairpin (L) that overlaps a Yn-Xm-AUG (pyrimidine tract/spacer/initiation codon) motif. SV-A, CKoV, and MKoV also contain these four domains, whereas BKV, SKV, and TV2/TV3 5' UTRs contain domains that are related to domain I and equivalent to domains J and K but lack an AV-like domain L. These IRESs are located at different relative positions between a conserved 5'-terminal origin of replication and divergent coding sequences. Elements in these IRESs also occur elsewhere: domain J's apical subdomain, which contains a GNRA tetraloop, matches an element in type 1 IRESs, and eIF4G-binding motifs in domain K and in type 2 IRESs are identical. Other elements are unique, and their presence leads to unique initiation factor requirements. In vitro reconstitution experiments showed that like AV, but in contrast to other currently characterized IRESs, SV-A requires the DExH-box protein DHX29 during initiation, which likely ensures that the initiation codon sequestered in domain L is properly accommodated in the ribosomal mRNA-binding cleft.

pre-miRNA profiles obtained through application of locked nucleic acids and deep sequencing reveals complex 5'/3' arm variation including concomitant cleavage and polyuridylation patterns

Recent research hints at an underappreciated complexity in pre-miRNA processing and regulation. Global profiling of pre-miRNA and its potential to increase understanding of the pre-miRNA landscape is impeded by overlap with highly expressed classes of other non coding (nc) RNA. Here, we present a data set excluding these RNA before sequencing through locked nucleic acids (LNA), greatly increasing pre-miRNA sequence counts with no discernable effect on pre-miRNA or mature miRNA sequencing. Analysis of profiles generated in total, nuclear and cytoplasmic cell fractions reveals that pre-miRNAs are subject to a wide range of regulatory processes involving loci-specific 3′- and 5′-end variation entailing complex cleavage patterns with co-occurring polyuridylation. Additionally, examination of nuclear-enriched flanking sequences of pre-miRNA, particularly those derived from polycistronic miRNA transcripts, provides insight into miRNA and miRNA-offset (moRNA) production, specifically identifying novel classes of RNA potentially functioning as moRNA precursors. Our findings point to particularly intricate regulation of the let-7 family in many ways reminiscent of DICER1-independent, pre-mir-451-like processing, introduce novel and unify known forms of pre-miRNA regulation and processing, and shed new light on overlooked products of miRNA processing pathways.

Association of TLR7 variants with AIDS-like disease and AIDS vaccine efficacy in rhesus macaques

In HIV infection, TLR7-triggered IFN-α production exerts a direct antiviral effect through the inhibition of viral replication, but may also be involved in immune pathogenesis leading to AIDS. TLR7 could also be an important mediator of vaccine efficacy. In this study, we analyzed polymorphisms in the X-linked TLR7 gene in the rhesus macaque model of AIDS. Upon resequencing of the TLR7 gene in 36 rhesus macaques of Indian origin, 12 polymorphic sites were detected. Next, we identified three tightly linked single nucleotide polymorphisms (SNP) as being associated with survival time. Genotyping of 119 untreated, simian immunodeficiency virus (SIV)-infected male rhesus macaques, including an 'MHC adjusted' subset, revealed that the three TLR7 SNPs are also significantly associated with set-point viral load. Surprisingly, this effect was not observed in 72 immunized SIV-infected male monkeys. We hypothesize (i) that SNP c.13G>A in the leader peptide is causative for the observed genotype-phenotype association and that (ii) the underlying mechanism is related to RNA secondary structure formation. Therefore, we investigated a fourth SNP (c.-17C>T), located 17 bp upstream of the ATG translation initiation codon, that is also potentially capable of influencing RNA structure. In c.13A carriers, neither set-point viral load nor survival time were related to the c.-17C>T genotype. In c.13G carriers, by contrast, the c.-17C allele was significantly associated with prolonged survival. Again, no such association was detected among immunized SIV-infected macaques. Our results highlight the dual role of TLR7 in immunodeficiency virus infection and vaccination and imply that it may be important to control human AIDS vaccine trials, not only for MHC genotype, but also for TLR7 genotype.

Systematic analysis of mitochondrial genes associated with hearing loss in the Japanese population: dHPLC reveals a new candidate mutation

Background

Variants of mitochondrial DNA (mtDNA) have been evaluated for their association with hearing loss. Although ethnic background affects the spectrum of mtDNA variants, systematic mutational analysis of mtDNA in Japanese patients with hearing loss has not been reported.

Methods

Using denaturing high-performance liquid chromatography combined with direct sequencing and cloning-sequencing, Japanese patients with prelingual (N = 54) or postlingual (N = 80) sensorineural hearing loss not having pathogenic mutations of m.1555A > G and m.3243A > G nor GJB2 were subjected to mutational analysis of mtDNA genes (12S rRNA, tRNA^Leu(UUR), tRNA^Ser(UCN), tRNA^Lys, tRNA^His, tRNA^Ser(AGY), and tRNA^Glu).

Results

We discovered 15 variants in 12S rRNA and one homoplasmic m.7501A > G variant in tRNA^Ser(UCN); no variants were detected in the other genes. Two criteria, namely the low frequency in the controls and the high conservation among animals, selected the m.904C > T and the m.1105T > C variants in 12S rRNA as candidate pathogenic mutations. Alterations in the secondary structures of the two variant transcripts as well as that of m.7501A > G in tRNA^Ser(UCN) were predicted.

Conclusions

The m.904C > T variant was found to be a new candidate mutation associated with hearing loss. The m.1105T > C variant is unlikely to be pathogenic. The pathogenicity of the homoplasmic m.7501T > A variant awaits further study.

The mechanism of translation initiation on Aichivirus RNA mediated by a novel type of picornavirus IRES

Picornavirus mRNAs contain IRESs that sustain their translation during infection, when host protein synthesis is shut off. The major classes of picornavirus IRESs (Types 1 and 2) have distinct structures and sequences, but initiation on both is determined by their specific interaction with eIF4G. We report here that Aichivirus (AV), a member of the Kobuvirus genus of Picornaviridae, contains an IRES that differs structurally from Type 1 and Type 2 IRESs. Its function similarly involves interaction with eIF4G, but its eIF4G-interacting domain is structurally distinct, although it contains an apical eIF4G-interacting motif similar to that in Type 2 IRESs. Like Type 1 and Type 2 IRESs, AV IRES function is enhanced by pyrimidine tract-binding protein (PTB), but the pattern of PTB's interaction with each of these IRESs is distinct. Unlike all known IRESs, the AV IRES is absolutely dependent on DHX29, a requirement imposed by sequestration of its initiation codon in a stable hairpin.

A NEW SPECIES OF VOLVOX SECT. MERRILLOSPHAERA (VOLVOCACEAE, CHLOROPHYCEAE) FROM TEXAS1

Smith (1944) divided the familiar genus Volvox L. into four sections, placing seven species that lacked cytoplasmic bridges between adult cells in the section Merrillosphaera. Herein, we describe a new member of the section Merrillosphaera originating from Texas (USA): Volvox ovalis Pocock ex Nozaki et A. W. Coleman sp. nov. Asexual spheroids of V. ovalis are ovoid or elliptical, with a monolayer of 1,000-2,000 somatic cells that are not linked by cytoplasmic bridges, an expanded anterior region, and 8-12 gonidia in the posterior region. Visibly asymmetric cleavage divisions do not occur in V. ovalis embryos as they do Volvox carteri F. Stein, Volvox obversus (W. Shaw) Printz, and Volvox africanus G. S. West, so the gonidia of the next generation are not yet recognizable in V. ovalis embryos prior to inversion. Molecular phylogenetic analyses of the five chloroplast genes and the internal transcribed spacer (ITS) regions of nuclear rDNA indicated that V. ovalis is closely related to Volvox spermatosphaera Powers (Powers 1908, as "spermatosphara") and/or Volvox tertius Art. Mey.; however, V. ovalis can be distinguished from V. spermatosphaera by its larger gonidia, and from V. tertius by visible differences in gonidial chloroplast morphology.

CentroidHomfold-LAST: accurate prediction of RNA secondary structure using automatically collected homologous sequences

Although secondary structure predictions of an individual RNA sequence have been widely used in a number of sequence analyses of RNAs, accuracy is still limited. Recently, we proposed a method (called 'CentroidHomfold'), which includes information about homologous sequences into the prediction of the secondary structure of the target sequence, and showed that it substantially improved the performance of secondary structure predictions. CentroidHomfold, however, forces users to prepare homologous sequences of the target sequence. We have developed a Web application (CentroidHomfold-LAST) that predicts the secondary structure of the target sequence using automatically collected homologous sequences. LAST, which is a fast and sensitive local aligner, and CentroidHomfold are employed in the Web application. Computational experiments with a commonly-used data set indicated that CentroidHomfold-LAST substantially outperformed conventional secondary structure predictions including CentroidFold and RNAfold.

Uniparental inheritance of chloroplast DNA is strict in the isogamous volvocalean Gonium

Background

A problem has remained unresolved regarding the exceptions to the unilateral inheritance of chloroplast DNA (cpDNA) from MT+/female in Chlamydomonas and other volvocaleans demonstrated by the previous genetic analyses. For identification of the parental types of cpDNA, these studies used parents that have differences in restriction fragment length polymorphisms and exhibit partial sexual incompatibility.

Methodology/Principal Findings

In the present study, we used sexually compatible parents of the isogamous colonial volvocalean Gonium maiaprilis that seemed an ideal species to identify the pattern of cpDNA inheritance based on the length difference in the putative group I intron interrupted in the Rubisco large subunit gene and objective identification of mating types by the presence or absence of the minus-dominance (MID) gene. We examined patterns of inheritance of cpDNA and presence/absence of a MID ortholog (GmMID) in 107 F₁ progeny of G. maiaprilis that were obtained by inducing germination of separated single zygotes. The results demonstrated no exception of the uniparental inheritance of cpDNA from the MT+ parent (lacking GmMID) in sexually compatible or genetically less divergent strains of G. maiaprilis.

Conclusions/Significance

The present data suggest that the uniparental inheritance of cpDNA is likely more strict in crossings of less diverged strains or sexually compatible parental volvocaleans, and some genetic inconsistency between the parents may cause exceptional uniparental inheritance of cpDNA.

Generalized centroid estimators in bioinformatics

In a number of estimation problems in bioinformatics, accuracy measures of the target problem are usually given, and it is important to design estimators that are suitable to those accuracy measures. However, there is often a discrepancy between an employed estimator and a given accuracy measure of the problem. In this study, we introduce a general class of efficient estimators for estimation problems on high-dimensional binary spaces, which represent many fundamental problems in bioinformatics. Theoretical analysis reveals that the proposed estimators generally fit with commonly-used accuracy measures (e.g. sensitivity, PPV, MCC and F-score) as well as it can be computed efficiently in many cases, and cover a wide range of problems in bioinformatics from the viewpoint of the principle of maximum expected accuracy (MEA). It is also shown that some important algorithms in bioinformatics can be interpreted in a unified manner. Not only the concept presented in this paper gives a useful framework to design MEA-based estimators but also it is highly extendable and sheds new light on many problems in bioinformatics.

Two-step cleavage of hairpin RNA with 5' overhangs by human DICER

Background

DICER is an RNase III family endoribonuclease that processes precursor microRNAs (pre-miRNAs) and long double-stranded RNAs, generating microRNA (miRNA) duplexes and short interfering RNA duplexes with 20~23 nucleotides (nts) in length. The typical form of pre-miRNA processed by the Drosha protein is a hairpin RNA with 2-nt 3' overhangs. On the other hand, production of mature miRNA from an endogenous hairpin RNA with 5' overhangs has also been reported, although the mechanism for this process is unknown.

Results

In this study, we show that human recombinant DICER protein (rDICER) processes a hairpin RNA with 5' overhangs in vitro and generates an intermediate duplex with a 29 nt-5' strand and a 23 nt-3' strand, which was eventually cleaved into a canonical miRNA duplex via a two-step cleavage. The previously identified endogenous pre-miRNA with 5' overhangs, pre-mmu-mir-1982 RNA, is also determined to be a substrate of rDICER through the same two-step cleavage.

Conclusions

The two-step cleavage of a hairpin RNA with 5' overhangs shows that DICER releases double-stranded RNAs after the first cleavage and binds them again in the inverse direction for a second cleavage. These findings have implications for how DICER may be able to interact with or process differing precursor structures.

Biochemical identification of A-to-I RNA editing sites by the inosine chemical erasing (ICE) method

Adenosine-to-inosine (A-to-I) RNA editing is a biologically important posttranscriptional processing event involved in the transcriptome diversification. The most conventional method of editing site identification is to compare the cDNA sequence with its corresponding genomic sequence; however, using this method, it is difficult to discriminate between guanosine residue that originated from inosine and errors or noise in the sequencing chromatograms. To address this issue, we developed the inosine chemical erasing (ICE) method to identify inosines in RNA strands utilizing inosine cyanoethylation and reverse transcription PCR. Since this method requires only a limited quantity of total RNA, it can be used in the genome-wide profiling of A-to-I editing sites in tissues and cells from various organisms, including clinical specimens.

RNA structure prediction

The prediction of RNA structure can be a first important step for the functional characterization of novel ncRNAs. Especially for the very meaningful secondary structure, there is a multitude of computational prediction tools. They differ not only in algorithmic details and the underlying models but also in what exactly they are trying to predict. This chapter gives an overview of different programs that aim to predict RNA secondary structure. We will introduce the ViennaRNA software package and web server as a solution that implements most of the varieties of RNA secondary structure prediction that have been published over the years. We focus on algorithms going beyond the mere prediction of a static structure.

Prediction of RNA secondary structure by maximizing pseudo-expected accuracy

Background

Recent studies have revealed the importance of considering the entire distribution of possible secondary structures in RNA secondary structure predictions; therefore, a new type of estimator is proposed including the maximum expected accuracy (MEA) estimator. The MEA-based estimators have been designed to maximize the expected accuracy of the base-pairs and have achieved the highest level of accuracy. Those methods, however, do not give the single best prediction of the structure, but employ parameters to control the trade-off between the sensitivity and the positive predictive value (PPV). It is unclear what parameter value we should use, and even the well-trained default parameter value does not, in general, give the best result in popular accuracy measures to each RNA sequence.

Results

Instead of using the expected values of the popular accuracy measures for RNA secondary structure prediction, which is difficult to be calculated, the pseudo-expected accuracy, which can easily be computed from base-pairing probabilities, is introduced. It is shown that the pseudo-expected accuracy is a good approximation in terms of sensitivity, PPV, MCC, or F-score. The pseudo-expected accuracy can be approximately maximized for each RNA sequence by stochastic sampling. It is also shown that well-balanced secondary structures between sensitivity and PPV can be predicted with a small computational overhead by combining the pseudo-expected accuracy of MCC or F-score with the γ-centroid estimator.

Conclusions

This study gives not only a method for predicting the secondary structure that balances between sensitivity and PPV, but also a general method for approximately maximizing the (pseudo-)expected accuracy with respect to various evaluation measures including MCC and F-score.

Improving the accuracy of predicting secondary structure for aligned RNA sequences

Considerable attention has been focused on predicting the secondary structure for aligned RNA sequences since it is useful not only for improving the limiting accuracy of conventional secondary structure prediction but also for finding non-coding RNAs in genomic sequences. Although there exist many algorithms of predicting secondary structure for aligned RNA sequences, further improvement of the accuracy is still awaited. In this article, toward improving the accuracy, a theoretical classification of state-of-the-art algorithms of predicting secondary structure for aligned RNA sequences is presented. The classification is based on the viewpoint of maximum expected accuracy (MEA), which has been successfully applied in various problems in bioinformatics. The classification reveals several disadvantages of the current algorithms but we propose an improvement of a previously introduced algorithm (CentroidAlifold). Finally, computational experiments strongly support the theoretical classification and indicate that the improved CentroidAlifold substantially outperforms other algorithms.