MiR-103a-3p targets the 5\u27 UTR of GPRC5A in pancreatic cells.

MicroRNAs (miRNAs) are short noncoding RNAs that regulate the expression of their targets in a sequence-dependent manner. For protein-coding transcripts, miRNAs regulate expression levels through binding sites in either the 3\u27 untranslated region (3\u27 UTR) or the amino acid coding sequence (CDS) of the targeted messenger RNA (mRNA). Currently, for the 5\u27 untranslated region (5\u27 UTR) of mRNAs, very few naturally occurring examples exist whereby the targeting miRNA down-regulates the expression of the corresponding mRNA in a seed-dependent manner. Here we describe and characterize two miR-103a-3p target sites in the 5\u27 UTR of GPRC5A, a gene that acts as a tumor suppressor in some cancer contexts and as an ongocene in other cancer contexts. In particular, we show that the interaction of miR-103a-3p with each of these two 5\u27 UTR targets reduces the expression levels of both GPRC5A mRNA and GPRC5A protein in one normal epithelial and two pancreatic cancer cell lines. By ectopically expressing sponges that contain instances of the wild-type 5\u27 UTR targets we also show that we can reduce miR-103a-3p levels and increase GPRC5A mRNA and protein levels. These findings provide some first knowledge on the post-transcriptional regulation of this tumor suppressor/oncogene and present additional evidence for the participation of 5\u27 UTRs in miRNA driven post-transcriptional regulatory control

Pattern-based phylogenetic distance estimation and tree reconstruction

We have developed an alignment-free method that calculates phylogenetic distances using a maximum likelihood approach for a model of sequence change on patterns that are discovered in unaligned sequences. To evaluate the phylogenetic accuracy of our method, and to conduct a comprehensive comparison of existing alignment-free methods (freely available as Python package decaf+py at http://www.bioinformatics.org.au), we have created a dataset of reference trees covering a wide range of phylogenetic distances. Amino acid sequences were evolved along the trees and input to the tested methods; from their calculated distances we infered trees whose topologies we compared to the reference trees. We find our pattern-based method statistically superior to all other tested alignment-free methods on this dataset. We also demonstrate the general advantage of alignment-free methods over an approach based on automated alignments when sequences violate the assumption of collinearity. Similarly, we compare methods on empirical data from an existing alignment benchmark set that we used to derive reference distances and trees. Our pattern-based approach yields distances that show a linear relationship to reference distances over a substantially longer range than other alignment-free methods. The pattern-based approach outperforms alignment-free methods and its phylogenetic accuracy is statistically indistinguishable from alignment-based distances.Comment: 21 pages, 3 figures, 2 table

A new computational method for the detection of horizontal gene transfer events

In recent years, the increase in the amounts of available genomic data has made it easier to appreciate the extent by which organisms increase their genetic diversity through horizontally transferred genetic material. Such transfers have the potential to give rise to extremely dynamic genomes where a significant proportion of their coding DNA has been contributed by external sources. Because of the impact of these horizontal transfers on the ecological and pathogenic character of the recipient organisms, methods are continuously sought that are able to computationally determine which of the genes of a given genome are products of transfer events. In this paper, we introduce and discuss a novel computational method for identifying horizontal transfers that relies on a gene's nucleotide composition and obviates the need for knowledge of codon boundaries. In addition to being applicable to individual genes, the method can be easily extended to the case of clusters of horizontally transferred genes. With the help of an extensive and carefully designed set of experiments on 123 archaeal and bacterial genomes, we demonstrate that the new method exhibits significant improvement in sensitivity when compared to previously published approaches. In fact, it achieves an average relative improvement across genomes of between 11 and 41% compared to the Codon Adaptation Index method in distinguishing native from foreign genes. Our method's horizontal gene transfer predictions for 123 microbial genomes are available online at

A sensitive, support-vector-machine method for the detection of horizontal gene transfers in viral, archaeal and bacterial genomes

In earlier work, we introduced and discussed a generalized computational framework for identifying horizontal transfers. This framework relied on a gene's nucleotide composition, obviated the need for knowledge of codon boundaries and database searches, and was shown to perform very well across a wide range of archaeal and bacterial genomes when compared with previously published approaches, such as Codon Adaptation Index and C + G content. Nonetheless, two considerations remained outstanding: we wanted to further increase the sensitivity of detecting horizontal transfers and also to be able to apply the method to increasingly smaller genomes. In the discussion that follows, we present such a method, Wn-SVM, and show that it exhibits a very significant improvement in sensitivity compared with earlier approaches. Wn-SVM uses a one-class support-vector machine and can learn using rather small training sets. This property makes Wn-SVM particularly suitable for studying small-size genomes, similar to those of viruses, as well as the typically larger archaeal and bacterial genomes. We show experimentally that the new method results in a superior performance across a wide range of organisms and that it improves even upon our own earlier method by an average of 10% across all examined genomes. As a small-genome case study, we analyze the genome of the human cytomegalovirus and demonstrate that Wn-SVM correctly identifies regions that are known to be conserved and prototypical of all beta-herpesvirinae, regions that are known to have been acquired horizontally from the human host and, finally, regions that had not up to now been suspected to be horizontally transferred. Atypical region predictions for many eukaryotic viruses, including the α-, β- and γ-herpesvirinae, and 123 archaeal and bacterial genomes, have been made available online at

YAMAT-seq: an efficient method for high-throughput sequencing of mature transfer RNAs.

Besides translation, transfer RNAs (tRNAs) play many non-canonical roles in various biological pathways and exhibit highly variable expression profiles. To unravel the emerging complexities of tRNA biology and molecular mechanisms underlying them, an efficient tRNA sequencing method is required. However, the rigid structure of tRNA has been presenting a challenge to the development of such methods. We report the development of Y-shaped Adapter-ligated MAture TRNA sequencing (YAMAT-seq), an efficient and convenient method for high-throughput sequencing of mature tRNAs. YAMAT-seq circumvents the issue of inefficient adapter ligation, a characteristic of conventional RNA sequencing methods for mature tRNAs, by employing the efficient and specific ligation of Y-shaped adapter to mature tRNAs using T4 RNA Ligase 2. Subsequent cDNA amplification and next-generation sequencing successfully yield numerous mature tRNA sequences. YAMAT-seq has high specificity for mature tRNAs and high sensitivity to detect most isoacceptors from minute amount of total RNA. Moreover, YAMAT-seq shows quantitative capability to estimate expression levels of mature tRNAs, and has high reproducibility and broad applicability for various cell lines. YAMAT-seq thus provides high-throughput technique for identifying tRNA profiles and their regulations in various transcriptomes, which could play important regulatory roles in translation and other biological processes

Knowledge about the presence or absence of miRNA isoforms (isomiRs) can successfully discriminate amongst 32 TCGA cancer types.

Isoforms of human miRNAs (isomiRs) are constitutively expressed with tissue- and disease-subtype-dependencies. We studied 10 271 tumor datasets from The Cancer Genome Atlas (TCGA) to evaluate whether isomiRs can distinguish amongst 32 TCGA cancers. Unlike previous approaches, we built a classifier that relied solely on \u27binarized\u27 isomiR profiles: each isomiR is simply labeled as \u27present\u27 or \u27absent\u27. The resulting classifier successfully labeled tumor datasets with an average sensitivity of 90% and a false discovery rate (FDR) of 3%, surpassing the performance of expression-based classification. The classifier maintained its power even after a 15× reduction in the number of isomiRs that were used for training. Notably, the classifier could correctly predict the cancer type in non-TCGA datasets from diverse platforms. Our analysis revealed that the most discriminatory isomiRs happen to also be differentially expressed between normal tissue and cancer. Even so, we find that these highly discriminating isomiRs have not been attracting the most research attention in the literature. Given their ability to successfully classify datasets from 32 cancers, isomiRs and our resulting \u27Pan-cancer Atlas\u27 of isomiR expression could serve as a suitable framework to explore novel cancer biomarkers

Aging-associated changes in microRNA expression profile of internal anal sphincter smooth muscle: Role of microRNA-133a.

A comprehensive genomic and proteomic, computational, and physiological approach was employed to examine the (previously unexplored) role of microRNAs (miRNAs) as regulators of internal anal sphincter (IAS) smooth muscle contractile phenotype and basal tone. miRNA profiling, genome-wide expression, validation, and network analyses were employed to assess changes in mRNA and miRNA expression in IAS smooth muscles from young vs. aging rats. Multiple miRNAs, including rno-miR-1, rno-miR-340-5p, rno-miR-185, rno-miR-199a-3p, rno-miR-200c, rno-miR-200b, rno-miR-31, rno-miR-133a, and rno-miR-206, were found to be upregulated in aging IAS. qPCR confirmed the upregulated expression of these miRNAs and downregulation of multiple, predicted targets (Eln, Col3a1, Col1a1, Zeb2, Myocd, Srf, Smad1, Smad2, Rhoa/Rock2, Fn1, Tagln v2, Klf4, and Acta2) involved in regulation of smooth muscle contractility. Subsequent studies demonstrated an aging-associated increase in the expression of miR-133a, corresponding decreases in RhoA, ROCK2, MYOCD, SRF, and SM22α protein expression, RhoA-signaling, and a decrease in basal and agonist [U-46619 (thromboxane A2 analog)]-induced increase in the IAS tone. Moreover, in vitro transfection of miR-133a caused a dose-dependent increase of IAS tone in strips, which was reversed by anti-miR-133a. Last, in vivo perianal injection of anti-miR-133a reversed the loss of IAS tone associated with age. This work establishes the important regulatory effect of miRNA-133a on basal and agonist-stimulated IAS tone. Moreover, reversal of age-associated loss of tone via anti-miR delivery strongly implicates miR dysregulation as a causal factor in the aging-associated decrease in IAS tone and suggests that miR-133a is a feasible therapeutic target in aging-associated rectoanal incontinence

Increasing cell density globally enhances the biogenesis of Piwi-interacting RNAs in Bombyx mori germ cells.

Piwi proteins and their bound Piwi-interacting RNAs (piRNAs) are predominantly expressed in the germline and play crucial roles in germline development by silencing transposons and other targets. Bombyx mori BmN4 cells are culturable germ cells that equip the piRNA pathway. Because of the scarcity of piRNA-expressing culturable cells, BmN4 cells are being utilized for the analyses of piRNA biogenesis. We here report that the piRNA biogenesis in BmN4 cells is regulated by cell density. As cell density increased, the abundance of Piwi proteins and piRNA biogenesis factors was commonly upregulated, resulting in an increased number of perinuclear nuage-like granules where Piwi proteins localize. Along with these phenomena, the abundance of mature piRNAs also globally increased, whereas levels of long piRNA precursor and transposons decreased, suggesting that increasing cell density promotes piRNA biogenesis pathway and that the resultant accumulation of mature piRNAs is functionally significant for transposon silencing. Our study reveals a previously uncharacterized link between cell density and piRNA biogenesis, designates cell density as a critical variable in piRNA studies using BmN4 cell system, and suggests the alteration of cell density as a useful tool to monitor piRNA biogenesis and function

IsomiR Expression Profiles in Human Lymphoblastoid Cell Lines Exhibit Population and Gender Dependencies.

For many years it was believed that each mature microRNA (miRNA) existed as a single entity with fixed endpoints and a \u27static\u27 and unchangeable primary sequence. However, recent evidence suggests that mature miRNAs are more \u27dynamic\u27 and that each miRNA precursor arm gives rise to multiple isoforms, the isomiRs. Here we report on our identification of numerous and abundant isomiRs in the lymphoblastoid cell lines (LCLs) of 452 men and women from five different population groups. Unexpectedly, we find that these isomiRs exhibit an expression profile that is population-dependent and gender-dependent. This is important as it indicates that the LCLs of each gender/population combination have their own unique collection of mature miRNA transcripts. Moreover, each identified isomiR has its own characteristic abundance that remains consistent across biological replicates indicating that these are not degradation products. The primary sequences of identified isomiRs differ from the known miRBase miRNA either at their 5´-endpoint (leads to a different \u27seed\u27 sequence and suggests a different targetome), their 3´-endpoint, or both simultaneously. Our analysis of Argonaute PAR-CLIP data from LCLs supports the association of many of these newly identified isomiRs with the Argonaute silencing complex and thus their functional roles through participation in the RNA interference pathway