N 6-methyladenosine enhances post-transcriptional gene regulation by microRNAs.

MOTIVATION: N 6-methyladenosine (m6A) is the most prevalent modification in eukaryotic messenger RNAs. MicroRNAs (miRNAs) are abundant post-transcriptional regulators of gene expression. Correlation between m6A and miRNA-targeting sites has been reported to suggest possible involvement of m6A in miRNA-mediated gene regulation. However, it is unknown what the regulatory effects might be. In this study, we performed comprehensive analyses of high-throughput data on m6A and miRNA target binding and regulation. RESULTS: We found that the level of miRNA-mediated target suppression is significantly enhanced when m6A is present on target mRNAs. The evolutionary conservation for miRNA-binding sites with m6A modification is significantly higher than that for miRNA-binding sites without modification. These findings suggest functional significance of m6A modification in post-transcriptional gene regulation by miRNAs. We also found that methylated targets have more stable structure than non-methylated targets, as indicated by significantly higher GC content. Furthermore, miRNA-binding sites that can be potentially methylated are significantly less accessible without methylation than those that do not possess potential methylation sites. Since either RNA-binding proteins or m6A modification by itself can destabilize RNA structure, we propose a model in which m6A alters local target secondary structure to increase accessibility for efficient binding by Argonaute proteins, leading to enhanced miRNA-mediated regulation. AVAILABILITY AND IMPLEMENTATION: N/A.

Sfold Tools for MicroRNA Target Prediction.

Computational prediction of miRNA binding sites on target mRNAs facilitates experimental investigation of miRNA functions. In this chapter, we describe STarMir and STarMirDB, two application modules of the Sfold RNA package. STarMir is a Web server for performing miRNA binding site predictions for mRNA and target sequences submitted by users. STarMirDB is a database of precomputed transcriptome-scale predictions. Both STarMir and STarMirDB provide comprehensive sequence, thermodynamic, and target structure features, a logistic probability as a measure of confidence for each predicted site, and a publication-quality diagram of the predicted miRNA-target hybrid. In addition, STarMir now offers a new quantitative score to address combined regulatory effects of multiple seed and seedless sites. This score provides a quantitative measure of the overall regulatory effects of both seed and seedless sites on the target. STarMir and STarMirDB are freely available to all through the Sfold Web application server at http://sfold.wadsworth.org .

Novel determinants of mammalian primary microRNA processing revealed by systematic evaluation of hairpin-containing transcripts and human genetic variation.

Mature microRNAs (miRNAs) are processed from hairpin-containing primary miRNAs (pri-miRNAs). However, rules that distinguish pri-miRNAs from other hairpin-containing transcripts in the genome are incompletely understood. By developing a computational pipeline to systematically evaluate 30 structural and sequence features of mammalian RNA hairpins, we report several new rules that are preferentially utilized in miRNA hairpins and govern efficient pri-miRNA processing. We propose that a hairpin stem length of 36 ± 3 nt is optimal for pri-miRNA processing. We identify two bulge-depleted regions on the miRNA stem, located ∼16-21 nt and ∼28-32 nt from the base of the stem, that are less tolerant of unpaired bases. We further show that the CNNC primary sequence motif selectively enhances the processing of optimal-length hairpins. We predict that a small but significant fraction of human single-nucleotide polymorphisms (SNPs) alter pri-miRNA processing, and confirm several predictions experimentally including a disease-causing mutation. Our study enhances the rules governing mammalian pri-miRNA processing and suggests a diverse impact of human genetic variation on miRNA biogenesis.

STarMir Tools for Prediction of microRNA Binding Sites.

MicroRNAs (miRNAs) are a class of endogenous short noncoding RNAs that regulate gene expression by targeting messenger RNAs (mRNAs), which results in translational repression and/or mRNA degradation. As regulatory molecules, miRNAs are involved in many mammalian biological processes and also in the manifestation of certain human diseases. As miRNAs play central role in the regulation of gene expression, understanding miRNA-binding patterns is essential to gain an insight of miRNA mediated gene regulation and also holds promise for therapeutic applications. Computational prediction of miRNA binding sites on target mRNAs facilitates experimental investigation of miRNA functions. This chapter provides protocols for using the STarMir web server for improved predictions of miRNA binding sites on a target mRNA. As an application module of the Sfold RNA package, the current version of STarMir is an implementation of logistic prediction models developed with high-throughput miRNA binding data from cross-linking immunoprecipitation (CLIP) studies. The models incorporated comprehensive thermodynamic, structural, and sequence features, and were found to make improved predictions of both seed and seedless sites, in comparison to the established algorithms (Liu et al., Nucleic Acids Res 41:e138, 2013). Their broad applicability was indicated by their good performance in cross-species validation. STarMir is freely available at http://sfold.wadsworth.org/starmir.html .

STarMirDB: A database of microRNA binding sites.

microRNAs (miRNAs) are an abundant class of small endogenous non-coding RNAs (ncRNAs) of ∼22 nucleotides (nts) in length. These small regulatory molecules are involved in diverse developmental, physiological and pathological processes. miRNAs target mRNAs (mRNAs) for translational repression and/or mRNA degradation. Predictions of miRNA binding sites facilitate experimental validation of miRNA targets. Models developed with data from CLIP studies have been used for predictions of miRNA binding sites in the whole transcriptomes of human, mouse and worm. The prediction results have been assembled into STarMirDB, a new database of miRNA binding sites available at http://sfold.wadsworth.org/starmirDB.php . STarMirDB can be searched by miRNAs or mRNAs separately or in combination. The search results are categorized into seed and seedless sites in 3' UTR, CDS and 5' UTR. For each predicted site, STarMirDB provides a comprehensive list of sequence, thermodynamic and target structural features that are known to influence miRNA: target interaction. A high resolution PDF diagram of the conformation of the miRNA:target hybrid is also available for visualization and publication. The results of a database search are available through both an interactive viewer and downloadable text files.

Effects of genetic variations on microRNA: target interactions.

Genetic variations within microRNA (miRNA) binding sites can affect miRNA-mediated gene regulation, which may lead to phenotypes and diseases. We perform a transcriptome-scale analysis of genetic variants and miRNA:target interactions identified by CLASH. This analysis reveals that rare variants tend to reside in CDSs, whereas common variants tend to reside in the 3' UTRs. miRNA binding sites are more likely to reside within those targets in the transcriptome with lower variant densities, especially target regions in which nucleotides have low mutation frequencies. Furthermore, an overwhelming majority of genetic variants within or near miRNA binding sites can alter not only the potential of miRNA:target hybridization but also the structural accessibility of the binding sites and flanking regions. These suggest an interpretation for certain associations between genetic variants and diseases, i.e. modulation of miRNA-mediated gene regulation by common or rare variants within or near miRNA binding sites, likely through target structure alterations. Our data will be valuable for discovering new associations among miRNAs, genetic variations and human diseases.

STarMir: a web server for prediction of microRNA binding sites.

STarMir web server predicts microRNA (miRNA) binding sites on a target ribonucleic acid (RNA). STarMir is an implementation of logistic prediction models developed with miRNA binding data from crosslinking immunoprecipitation (CLIP) studies (Liu,C., Mallick, B., Long, D., Rennie, W.A., Wolenc, A., Carmack, C.S. and Ding, Y. (2013). CLIP-based prediction of mammalian microRNA binding sites. Nucleic Acids Res., 41(14), e138). In both intra-dataset and inter-dataset validations, the models showed major improvements over established algorithms in predictions of both seed and seedless sites. General applicability of the models was indicated by good performance in cross-species validations. The input data for STarMir is processed by the web server to perform prediction of miRNA binding sites, compute comprehensive sequence, thermodynamic and target structure features and a logistic probability as a measure of confidence for each predicted site. For each of seed and seedless sites and for all three regions of a mRNA (3' UTR, CDS and 5' UTR), STarMir output includes the computed binding site features, the logistic probability and a publication-quality diagram of the predicted miRNA:target hybrid. The prediction results are available through both an interactive viewer and downloadable text files. As an application module of the Sfold RNA package (http://sfold.wadsworth.org), STarMir is freely available to all at http://sfold.wadsworth.org/starmir.html.

MicroRNA binding sites in C. elegans 3' UTRs.

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression. Since the discovery of lin-4, the founding member of the miRNA family, over 360 miRNAs have been identified for Caenorhabditis elegans (C. elegans). Prediction and validation of targets are essential for elucidation of regulatory functions of these miRNAs. For C. elegans, crosslinking immunoprecipitation (CLIP) has been successfully performed for the identification of target mRNA sequences bound by Argonaute protein ALG-1. In addition, reliable annotation of the 3' untranslated regions (3' UTRs) as well as developmental stage-specific expression profiles for both miRNAs and 3' UTR isoforms are available. By utilizing these data, we developed statistical models and bioinformatics tools for both transcriptome-scale and developmental stage-specific predictions of miRNA binding sites in C. elegans 3' UTRs. In performance evaluation via cross validation on the ALG-1 CLIP data, the models were found to offer major improvements over established algorithms for predicting both seed sites and seedless sites. In particular, our top-ranked predictions have a substantially higher true positive rate, suggesting a much higher likelihood of positive experimental validation. A gene ontology analysis of stage-specific predictions suggests that miRNAs are involved in dynamic regulation of biological functions during C. elegans development. In particular, miRNAs preferentially target genes related to development, cell cycle, trafficking, and cell signaling processes. A database for both transcriptome-scale and stage-specific predictions and software for implementing the prediction models are available through the Sfold web server at http://sfold.wadsworth.org.

A study on the influence of different promoter and 5'UTR (URM) cassettes from Arabidopsis thaliana on the expression level of the reporter gene ß glucuronidase in tobacco and cotton.

Several reports of promoters from plants, viral and artificial origin that confer high constitutive expression are known. Among these the CaMV 35S promoter is used extensively for transgene expression in plants. We identified candidate promoters from Arabidopsis based on their transcript levels (meta-analysis of available microarray control datasets) to test their activity in comparison to the CaMV 35S promoter. A set of 11 candidate genes were identified which showed high transcript levels in the aerial tissue (i.e. leaf, shoot, flower and stem). In the initial part of the study binary vectors were developed wherein the promoter and 5'UTR region of these candidate genes (Upstream Regulatory Module, URM) were cloned upstream to the reporter gene ß glucuronidase (gus). The promoter strengths were tested in transformed callus of Nicotiana tabacum and Gossypium hirsutum. On the basis of the results obtained from the callus, the influence of the URM cassettes on transgene expression was tested in transgenic tobacco. The URM regions of the genes encoding a subunit of photosystem I (PHOTO) and geranyl geranyl reductase (GGR) in A. thaliana genome showed significantly high levels of GUS activity in comparison to the CaMV 35S promoter. Further, when the 5'UTRs of both the genes were placed downstream to the CaMV 35S promoter it led to a substantial increase in GUS activity in transgenic tobacco lines and cotton callus. The enhancement observed was even higher to that observed with the viral leader sequences like Ω and AMV, known translational enhancers. Our results indicate that the two URM cassettes or the 5'UTR regions of PHOTO and GGR when placed downstream to the CaMV 35S promoter can be used to drive high levels of transgene expression in dicotyledons.

A 28 nt long synthetic 5'UTR (synJ) as an enhancer of transgene expression in dicotyledonous plants.

BACKGROUND: A high level of transgene expression is required, in several applications of transgenic technology. While use of strong promoters has been the main focus in such instances, 5'UTRs have also been shown to enhance transgene expression. Here, we present a 28 nt long synthetic 5'UTR (synJ), which enhances gene expression in tobacco and cotton. RESULTS: The influence of synJ on transgene expression was studied in callus cultures of cotton and different tissues of transgenic tobacco plants. The study was based on comparing the expression of reporter gene gus and gfp, with and without synJ as its 5'UTR. Mutations in synJ were also analyzed to identify the region important for enhancement. synJ, enhances gene expression by 10 to 50 fold in tobacco and cotton depending upon the tissue studied. This finding is based on the experiments comparing the expression of gus gene, encoding the synJ as 5'UTR under the control of 35S promoter with expression cassettes based on vectors like pBI121 or pRT100. Further, the enhancement was in most cases equivalent to that observed with the viral leader sequences known to enhance translation like Ω and AMV. In case of transformed cotton callus as well as in the roots of tobacco transgenic plants, the up-regulation mediated by synJ was much higher than that observed in the presence of both Ω as well as AMV. The enhancement mediated by synJ was found to be at the post-transcriptional level. The study also demonstrates the importance of a 5'UTR in realizing the full potential of the promoter strength. synJ has been utilized to design four cloning vectors: pGEN01, pBGEN02, pBGEN02-hpt and pBGEN02-ALSdm each of which can be used for cloning the desired transgene and achieving high level of expression in the resulting transgenic plants. CONCLUSIONS: synJ, a synthetic 5'UTR, can enhance transgene expression under a strong promoter like 35S as well as under a weak promoter like nos in dicotyledonous plants. synJ can be incorporated as the 5'UTR of transgenes, especially in cases where high levels of expression is required. A set of vectors has also been designed to facilitate this process.

Detrimental effect of expression of Bt endotoxin Cry1Ac on in vitro regeneration, in vivo growth and development of tobacco and cotton transgenics.

High levels of expression of the cry1Ac gene from Bacillus thuringiensis cannot be routinely achieved in transgenic plants despite modifications made in the gene to improve its expression. This has been attributed to the instability of the transcript in a few reports. In the present study, based on the genetic transformation of cotton and tobacco, we show that the expression of the Cry1Ac endotoxin has detrimental effects on both the in vitro and in vivo growth and development of transgenic plants. A number of experiments on developing transgenics in cotton with different versions of cry1Ac gene showed that the majority of the plants did not express any Cry1Ac protein. Based on Southern blot analysis, it was also observed that a substantial number of lines did not contain the cry1Ac gene cassette although they contained the marker gene nptII. More significantly, all the lines that showed appreciable levels of expression were found to be phenotypically abnormal. Experiments on transformation of tobacco with different constructs expressing the cry1Ac gene showed that in vitro regeneration was inhibited by the encoded protein. Further, out of a total of 145 independent events generated with the different cry1Ac gene constructs in tobacco, only 21 showed expression of the Cry1Ac protein, confirming observations made in cotton that regenerants that express high levels of the Cry1Ac protein are selected against during regeneration of transformed events. This problem was circumvented by targeting the Cry1Ac protein to the chloroplast, which also significantly improved the expression of the protein.