Transcription factors are targeted by differentially expressed miRNAs in primates

Identification of microRNAs Derived from Transposable Elements in the Macaca mulatta (Rhesus Monkey) Genome

Transposable elements (TEs) are mobile DNA entities that can move within the host genome. Over long periods of evolutionary time, TEs are typically silenced via the accumulation of mutations in the genome, ultimately resulting in their immobilization. However, they still play an important role in the host genome by acting as regulatory elements. They influence host transcription in various ways, one of which as the origin of the generation of microRNAs (miRNAs), which are so-called miRNAs derived from TEs (MDTEs). miRNAs are small non-coding RNAs that are involved in many biological processes by regulating gene expression at the post-transcriptional level. Here, we identified MDTEs in the Macaca mulatta (rhesus monkey) genome, which is phylogenetically close species to humans, based on the genome coordinates of miRNAs and TEs. The expression of 5 out of 17 MDTEs that were exclusively registered in M. mulatta from the miRBase database (v22) was examined via quantitative polymerase chain reaction (qPCR). Moreover, Gene Ontology analysis was performed to examine the functional implications of the putative target genes of the five MDTEs.

Analysis of conserved miRNAs in cynomolgus macaque genome using small RNA sequencing and homology searching

MicroRNAs (miRNAs) are important regulators that fine-tune diverse cellular activities. Cynomolgus macaques (Macaca fascicularis) are used extensively in biomedical and pharmaceutical research; however, substantially fewer miRNAs have been identified in this species than in humans. Consequently, we investigated conserved miRNA profiles in cynomolgus macaques by homology searching and small RNA sequencing. In total, 1,455 high-confidence miRNA gene loci were identified, 408 of which were also confirmed by RNA sequencing, including 73 new miRNA loci reported in cynomolgus macaques for the first time. Comparing miRNA expression with age, we found a positive correlation between sequence conservation and expression levels during miRNA evolution. Additionally, we found that the miRNA gene locations in cynomolgus macaque genome were very flexible. Most were embedded in intergenic spaces or introns and clustered together. Several miRNAs were found in certain gene locations, including 64 exon-resident miRNAs, six splice-site-overlapping miRNAs (SO-miRNAs), and two pairs of distinct mirror miRNAs. We also identified 78 miRNA clusters, 68 of which were conserved in the human genome, including 10 large miRNA clusters predicted to regulate diverse developmental and cellular processes in cynomolgus macaque. Thus, this study not only expands the number of identified miRNAs in cynomolgus macaques but also provides clues for future research on the differences in miRNA repertoire between macaques and humans.

Effects of Cardiac Sympathetic Neurodegeneration and PPARγ Activation on Rhesus Macaque Whole Blood miRNA and mRNA Expression Profiles

Degeneration of sympathetic innervation of the heart occurs in numerous diseases, including diabetes, idiopathic REM sleep disorder, and Parkinson's disease (PD). In PD, cardiac sympathetic denervation occurs in 80-90% of patients and can begin before the onset of motor symptoms. Today, there are no disease-modifying therapies for cardiac sympathetic neurodegeneration, and biomarkers are limited to radioimaging techniques. Analysis of expression levels of coding mRNA and noncoding RNAs, such as microRNAs (miRNAs), can uncover pathways involved in disease, leading to the discovery of biomarkers, pathological mechanisms, and potential drug targets. Whole blood in particular is a clinically relevant source of biomarkers, as blood sampling is inexpensive and simple to perform. Our research group has previously developed a nonhuman primate model of cardiac sympathetic denervation by intravenous administration of the catecholaminergic neurotoxin 6-hydroxydopamine (6-OHDA). In this rhesus macaque (Macaca mulatta) model, imaging with positron emission tomography showed that oral administration of the peroxisome proliferator-activated receptor gamma (PPARγ) agonist pioglitazone (n = 5; 5 mg/kg daily) significantly decreased cardiac inflammation and oxidative stress compared to placebo (n = 5). Here, we report our analysis of miRNA and mRNA expression levels over time in the whole blood of these monkeys. Differential expression of three miRNAs was induced by 6-OHDA (mml-miR-16-2-3p, mml-miR-133d-3p, and mml-miR-1262-5p) and two miRNAs by pioglitazone (mml-miR-204-5p and mml-miR-146b-5p) at 12 weeks posttoxin, while expression of mRNAs involved in inflammatory cytokines and receptors was not significantly affected. Overall, this study contributes to the characterization of rhesus coding and noncoding RNA profiles in normal and disease-like conditions, which may facilitate the identification and clinical translation of biomarkers of cardiac neurodegeneration and neuroprotection.

Gene regulatory network stabilized by pervasive weak repressions: microRNA functions revealed by the May-Wigner theory

Food web and gene regulatory networks (GRNs) are large biological networks, both of which can be analyzed using the May-Wigner theory. According to the theory, networks as large as mammalian GRNs would require dedicated gene products for stabilization. We propose that microRNAs (miRNAs) are those products. More than 30% of genes are repressed by miRNAs, but most repressions are too weak to have a phenotypic consequence. The theory shows that (i) weak repressions cumulatively enhance the stability of GRNs, and (ii) broad and weak repressions confer greater stability than a few strong ones. Hence, the diffuse actions of miRNAs in mammalian cells appear to function mainly in stabilizing GRNs. The postulated link between mRNA repression and GRN stability can be seen in a different light in yeast, which do not have miRNAs. Yeast cells rely on non-specific RNA nucleases to strongly degrade mRNAs for GRN stability. The strategy is suited to GRNs of small and rapidly dividing yeast cells, but not the larger mammalian cells. In conclusion, the May-Wigner theory, supplanting the analysis of small motifs, provides a mathematical solution to GRN stability, thus linking miRNAs explicitly to 'developmental canalization'.

Identification and expression analysis of miRNA in hybrid snakehead by deep sequencing approach and their targets prediction

MicroRNAs (miRNAs) play important regulatory roles in numerous biological processes, but there is no report on miRNAs of hybrid snakehead. In this study, four independent small RNA libraries were constructed from the spleen, liver kidney and muscle of hybrid snakehead. These libraries were sequenced using deep sequencing technology, as result, a total of 1,067,172, 1,152,002, 1,653,941 and 970,866 clean reads from these four libraries were obtained. 252 known miRNAs and 63 putative novel miRNAs in these small RNA dataset were identified. The stem-loop RT-qPCR analysis indicated that eight known miRNAs and two novel miRNAs show different expression in eight different kinds of tissues. For better understanding the functions of miRNAs, 95,947 predicated target genes were analyzed by GO and their pathways, the results indicated that these targets of the identified miRNAs are involved in a broad range of physiological functions.

Systems biology approach to study the role of miRNA in promoter targeting during megakaryopoiesis

The distinct process of megakaryopoiesis requires occurrence of endomitosis for polyploidization of the megakaryocytes. Although, Cyclins, CDKs and have been described to regulate endomitosis, the exact mechanism still remains an enigma. miRNA which were otherwise known as post transcriptional gene silencers are now emerging with various non-canonical functions including gene regulation at pre-transcriptional level by miRNA binding at promoter region. Out of the many processes they regulate, miRNA have been manifested to play a role in megakaryocyte differentiation. In this study an attempt has been made to identify miRNA that could regulate cell cycle genes (Cyclins and CDKs) by targeting their promoters, during megakaryopoiesis. A new computational algorithm was implemented using Perl programming to identify putative targets of miRNA in CDK and Cyclin promoters. Perl script was also used to check nuclear localizing miRNA based on the presence of a consensus sequence. Real-time PCR was performed to analyze the expression of miRNA and their predicted targets in Dami vs. PMA treated Dami cells. Putative targets of miRNAs with longest, high complementarity matches in CDK/Cyclin promoters were obtained. We identified two significant miRNA, miR-1273g-3p and miR-619-5p with longest seed sequence matches. We further identified three main targets (CDK10, CDK11, Cyclin F) through which these two miRNA could regulate cell cycle during megakaryopoiesis. Our results reinforce the role of promoting targeting miRNA in regulation of cell cycle through certain CDK/Cyclins to support the process of endomitosis during megakaryopoiesis.

Computational identification of piRNA targets on mouse mRNAs

MOTIVATION

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are highly abundant in the germline. One important role of piRNAs is to defend genome integrity by guiding PIWI proteins to silence transposable elements (TEs), which have a high potential to cause deleterious effects on their host. The mechanism of piRNA-mediated post-transcriptional silencing was also observed to affect mRNAs, suggesting that piRNAs might play a broad role in gene expression regulation. However, there has been no systematic report with regard to how many protein-coding genes might be targeted and regulated by piRNAs.

RESULTS

We trained a support vector machine classifier based on a combination of Miwi CLIP-Seq-derived features and position-derived features to predict the potential targets of piRNAs on mRNAs in the mouse. Reanalysis of a published microarray dataset suggested that the expression level of the 2587 protein-coding genes predicted as piRNA targets showed significant upregulation as a whole after abolishing the slicer activity of Miwi, supporting the conclusion that they are subject to piRNA-mediated regulation.

AVAILABILITY AND IMPLEMENTATION

A web version of the method called pirnaPre as well as our results for browse is available at http://www.regulatoryrna.org/software/piRNA/piRNA_target_mRNA/index.php

CONTACT

crs@sun5.ibp.ac.cn or heshunmin@gmail.com

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of microRNAs in Macaca fascicularis (Cynomolgus Monkey) by Homology Search and Experimental Validation by Small RNA-Seq and RT-qPCR Using Kidney Cortex Tissues

MicroRNAs (miRNAs) present in tissues and biofluids are emerging as sensitive and specific safety biomarkers. MiRNAs have not been thoroughly described in M. fascicularis, an animal model used in pharmaceutical industry especially in drug safety evaluation. Here we investigated the miRNAs in M. fascicularis. For Macaca mulatta, a closely related species of M. fascicularis, 619 stem-loop precursor miRNAs (pre-miRNAs) and 914 mature miRNAs are available in miRBase version 21. Using M. mulatta miRNAs as a reference list and homology search tools, we identified 604 pre-miRNAs and 913 mature miRNAs in the genome of M. fascicularis. In order to validate the miRNAs identified by homology search we attempted to sequence miRNAs expressed in kidney cortex from M. fascicularis. MiRNAs expressed in kidney cortex may indeed be released in urine upon kidney cortex damage and be potentially used to monitor drug induced kidney injury. Hence small RNA sequencing libraries were prepared using kidney cortex tissues obtained from three naive M. fascicularis and sequenced. Analysis of sequencing data indicated that 432 out of 913 mature miRNAs were expressed in kidney cortex tissues. Assigning these 432 miRNAs to pre-miRNAs revealed that 273 were expressed from both the -5p and -3p arms of 150 pre-miRNAs and 159 miRNAs expressed from either the -5p or -3p arm of 176 pre-miRNAs. Mapping sequencing reads to pre-miRNAs also facilitated the detection of twenty-two new miRNAs. To substantiate miRNAs identified by small RNA sequencing, 313 miRNAs were examined by RT-qPCR. Expression of 262 miRNAs in kidney cortex tissues ware confirmed by TaqMan microRNA RT-qPCR assays. Analysis of kidney cortex miRNA targeted genes suggested that they play important role in kidney development and function. Data presented in this study may serve as a valuable resource to assess the renal safety biomarker potential of miRNAs in Cynomolgus monkeys.

Lineage-Specific Changes in Biomarkers in Great Apes and Humans

Although human biomedical and physiological information is readily available, such information for great apes is limited. We analyzed clinical chemical biomarkers in serum samples from 277 wild- and captive-born great apes and from 312 healthy human volunteers as well as from 20 rhesus macaques. For each individual, we determined a maximum of 33 markers of heart, liver, kidney, thyroid and pancreas function, hemoglobin and lipid metabolism and one marker of inflammation. We identified biomarkers that show differences between humans and the great apes in their average level or activity. Using the rhesus macaques as an outgroup, we identified human-specific differences in the levels of bilirubin, cholinesterase and lactate dehydrogenase, and bonobo-specific differences in the level of apolipoprotein A-I. For the remaining twenty-nine biomarkers there was no evidence for lineage-specific differences. In fact, we find that many biomarkers show differences between individuals of the same species in different environments. Of the four lineage-specific biomarkers, only bilirubin showed no differences between wild- and captive-born great apes. We show that the major factor explaining the human-specific difference in bilirubin levels may be genetic. There are human-specific changes in the sequence of the promoter and the protein-coding sequence of uridine diphosphoglucuronosyltransferase 1 (UGT1A1), the enzyme that transforms bilirubin and toxic plant compounds into water-soluble, excretable metabolites. Experimental evidence that UGT1A1 is down-regulated in the human liver suggests that changes in the promoter may be responsible for the human-specific increase in bilirubin. We speculate that since cooking reduces toxic plant compounds, consumption of cooked foods, which is specific to humans, may have resulted in relaxed constraint on UGT1A1 which has in turn led to higher serum levels of bilirubin in humans.

The Role of miR-378a in Metabolism, Angiogenesis, and Muscle Biology

MicroRNA-378a (miR-378a, previously known as miR-378) is one of the small noncoding RNA molecules able to regulate gene expression at posttranscriptional level. Its two mature strands, miR-378a-3p and miR-378a-5p, originate from the first intron of the peroxisome proliferator-activated receptor gamma, coactivator 1 beta (ppargc1b) gene encoding PGC-1β. Embedding in the sequence of this transcriptional regulator of oxidative energy metabolism implies involvement of miR-378a in metabolic pathways, mitochondrial energy homeostasis, and related biological processes such as muscle development, differentiation, and regeneration. On the other hand, modulating the expression of proangiogenic factors such as vascular endothelial growth factor, angiopoietin-1, or interleukin-8, influencing inflammatory reaction, and affecting tumor suppressors, such as SuFu and Fus-1, miR-378a is considered as a part of an angiogenic network in tumors. In the latter, miR-378a can evoke broader actions by enhancing cell survival, reducing apoptosis, and promoting cell migration and invasion. This review describes the current knowledge on miR-378a linking oxidative/lipid metabolism, muscle biology, and blood vessel formation.

Evolutionary dynamics of coding and non-coding transcriptomes

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

A non-canonical landscape of the microRNA system

Microribonucleic acids, best known as microRNAs or miRNAs, are small, non-coding RNAs with important regulatory roles in eukaryotic cells. Here, I present a broad review on highly relevant but generally non-depicted features of miRNAs, among which stand out the non-conventional miRNA seed sites, the unusual messenger RNA (mRNA) target regions, the non-canonical miRNA-guided mechanisms of gene expression regulation, and the recently identified new class of miRNA ligands. Furthermore, I address the miRNA uncommon genomic location, transcription, and subcellular localization. Altogether, these unusual features and roles place the miRNA system as a very diverse, complex, and intriguing biological mechanism.

The role of gene regulatory factors in the evolutionary history of humans

Deciphering the molecular basis of how modern human phenotypes have evolved is one of the most fascinating challenges in biology. Here, we will focus on the roles of gene regulatory factors (GRFs), in particular transcription factors (TFs) and long non-coding RNAs (lncRNAs) during human evolution. We will present examples of TFs and lncRNAs that have changed or show signs of positive selection in humans compared to chimpanzees, in modern humans compared to archaic humans, or within modern human populations. On the basis of current knowledge about the functions of these GRF genes, we speculate that they have been involved in speciation as well as in shaping phenotypes such as brain functions, skeletal morphology, and metabolic processes.

In transition: primate genomics at a time of rapid change

The field of nonhuman primate genomics is undergoing rapid change and making impressive progress. Exploiting new technologies for DNA sequencing, researchers have generated new whole-genome sequence assemblies for multiple primate species over the past 6 years. In addition, investigations of within-species genetic variation, gene expression and RNA sequences, conservation of non-protein-coding regions of the genome, and other aspects of comparative genomics are moving at an accelerating speed. This progress is opening a wide array of new research opportunities in the analysis of comparative primate genome content and evolution. It also creates new possibilities for the use of nonhuman primates as model organisms in biomedical research. This transition, based on both new technology and the new information being generated in regard to human genetics, provides an important justification for reevaluating the research goals, strategies, and study designs used in primate genetics and genomics.

Comparative primate genomics: emerging patterns of genome content and dynamics

Advances in genome sequencing technologies have created new opportunities for comparative primate genomics. Genome assemblies have been published for various primate species, and analyses of several others are underway. Whole-genome assemblies for the great apes provide remarkable new information about the evolutionary origins of the human genome and the processes involved. Genomic data for macaques and other non-human primates offer valuable insights into genetic similarities and differences among species that are used as models for disease-related research. This Review summarizes current knowledge regarding primate genome content and dynamics, and proposes a series of goals for the near future.

Human-specific microRNA regulation of FOXO1: implications for microRNA recognition element evolution

MicroRNAs (miRNAs) have been established as important negative post-transcriptional regulators for gene expression. Within the past decade, miRNAs targeting transcription factors (TFs) has emerged as an important mechanism for gene expression regulation. Here, we tested the hypothesis that in TF 3'UTRs, human-specific single nucleotide change(s) that create novel miRNA recognition elements (MREs) contribute to species-specific differences in TF expression. From several potential human-specific TF MREs, one candidate, a member of the Forkhead Box O (FOXO) subclass in the Forkhead family known as Forkhead Box O1 (FOXO1; FKHR; NM_002015) was tested further. Human FOXO1 contains two sites predicted to confer miR-183-mediated post-transcriptional regulation: one specific to humans and the other conserved. Utilizing dual luciferase expression reporters, we show that only the human FOXO1 3'UTR contains a functional miR-183 site, not found in chimpanzee or mouse 3'untranslated regions (UTRs). Site-directed mutagenesis supports functionality of the human-specific miR-183 site, but not the conserved miR-183 site. Via overexpression and target site protection assays, we show that human FOXO1 is regulated by miR-183, but mouse FOXO1 is not. Finally, FOXO1-regulated cellular phenotypes, including cell invasion and proliferation, are impacted by miR-183 targeting only in human cells. These results provide strong evidence for human-specific gain of TF MREs, a process that may underlie evolutionary differences between phylogenic groups.

Genome-wide characterization of microRNA in foxtail millet (Setaria italica)

BACKGROUND

MicroRNAs (miRNAs) are a class of short non-coding, endogenous RNAs that play key roles in many biological processes in both animals and plants. Although many miRNAs have been identified in a large number of organisms, the miRNAs in foxtail millet (Setaria italica) have, until now, been poorly understood.

RESULTS

In this study, two replicate small RNA libraries from foxtail millet shoots were sequenced, and 40 million reads representing over 10 million unique sequences were generated. We identified 43 known miRNAs, 172 novel miRNAs and 2 mirtron precursor candidates in foxtail millet. Some miRNA*s of the known and novel miRNAs were detected as well. Further, eight novel miRNAs were validated by stem-loop RT-PCR. Potential targets of the foxtail millet miRNAs were predicted based on our strict criteria. Of the predicted target genes, 79% (351) had functional annotations in InterPro and GO analyses, indicating the targets of the miRNAs were involved in a wide range of regulatory functions and some specific biological processes. A total of 69 pairs of syntenic miRNA precursors that were conserved between foxtail millet and sorghum were found. Additionally, stem-loop RT-PCR was conducted to confirm the tissue-specific expression of some miRNAs in the four tissues identified by deep-sequencing.

CONCLUSIONS

We predicted, for the first time, 215 miRNAs and 447 miRNA targets in foxtail millet at a genome-wide level. The precursors, expression levels, miRNA* sequences, target functions, conservation, and evolution of miRNAs we identified were investigated. Some of the novel foxtail millet miRNAs and miRNA targets were validated experimentally.

Species-specific microRNA regulation influences phenotypic variability: perspectives on species-specific microRNA regulation

Phenotypic divergence among animal species may be due in part to species-specific (SS) regulation of gene expression by small, non-coding regulatory RNAs termed "microRNAs". This phenomenon can be modulated by several variables. First, microRNA genes vary by their level of conservation, many of them being SS, or unique to a particular evolutionary lineage. Second, microRNA expression levels vary spatially and temporally in different species. Lastly, while microRNAs bind the 3'UTR of target genes in order to silence their expression, the binding sites themselves are often non-conserved. The variability of the miRNA-target paradigm between different species is thus multifactorial, and this paradigm has only just started to gain attention from researchers in various fields. Here we present and discuss recent findings regarding the characteristics and implications of SS microRNA regulation.

ROCK: a resource for integrative breast cancer data analysis

Given the steady increase in breast cancer rates in both the developed and developing world, there has been a concerted research effort undertaken worldwide to understand the molecular mechanisms underpinning the disease. The data generated from numerous clinical trials and experimental studies shed light on different aspects of the disease. We present a new version of the ROCK database (rock.icr.ac.uk), which integrates such diverse data types allowing unique analyses of published breast cancer experimental data. We have added several new data types and analysis modules to ROCK, which allow the user to interactively query and research the huge amounts of available experimental data and perform complex correlations across studies and data types such as gene expression, genomic copy number aberrations, micro RNA expression, RNA interference, survival analysis, clinical annotation and signalling protein networks. We present the recent and major functional updates and enhancements to the ROCK resource, including new analysis modules and microRNA and NGS data integration, and illustrate how ROCK can be used to confirm known experimental results as well as generate novel leads and new experimental hypotheses using the Wnt signalling cell surface receptor FZD7 and the Myc oncogene. ROCK provides a unique breast cancer analysis platform of integrated experimental datasets at the genomic, transcriptomic and proteomic level. This paper presents how ROCK has transitioned from being simply a database to an interactive resource useful to the broader breast cancer research community in our effort to facilitate research into the underlying molecular mechanisms of breast cancer.

Two panels of plasma microRNAs as non-invasive biomarkers for prediction of recurrence in resectable NSCLC

The diagnosis of non-small cell lung carcinoma (NSCLC) at an early stage, as well as better prediction of outcome remains clinically challenging due to the lack of specific and robust non-invasive markers. The discovery of microRNAs (miRNAs), particularly those found in the bloodstream, has opened up new perspectives for tumor diagnosis and prognosis. The aim of our study was to determine whether expression profiles of specific miRNAs in plasma could accurately discriminate between NSCLC patients and controls, and whether they are able to predict the prognosis of resectable NSCLC patients. We therefore evaluated a series of seventeen NSCLC-related miRNAs by quantitative real-time (qRT)-PCR in plasma from 52 patients with I-IIIA stages NSCLC, 10 patients with chronic obstructive pulmonary disease (COPD) and 20-age, sex and smoking status-matched healthy individuals. We identified an eleven-plasma miRNA panel that could distinguish NSCLC patients from healthy subjects (AUC = 0.879). A six-plasma miRNA panel was able to discriminate between NSCLC patients and COPD patients (AUC = 0.944). Furthermore, we identified a three-miRNA plasma signature (high miR-155-5p, high miR-223-3p, and low miR-126-3p) that significantly associated with a higher risk for progression in adenocarcinoma patients. In addition, a three-miRNA plasma panel (high miR-20a-5p, low miR-152-3p, and low miR-199a-5p) significantly predicted survival of squamous cell carcinoma patients. In conclusion, we identified two plasma miRNA expression profiles that may be useful for predicting the outcome of patients with resectable NSCLC.

MicroRNAs in opioid addiction: elucidating evolution

Three reviews in the Frontiers Research Topic “Non-Coding RNA and Addiction” (He and Wang, 2012; Rodriguez, 2012; Zheng et al., 2012), grouped under the chapter “MicroRNAs and Morphine,” focus on the contribution of microRNAs to opioid abuse. Although animal models have been fundamental to our understanding of addiction pathways, the assumption that microRNAs implicated in opioid tolerance – and their binding sites in mRNAs – are conserved in mammalian evolution was not examined by the authors. Inspired by recent reports which highlight a surprising lack of evolutionary conservation in non-coding RNA genes, in this perspective we use public genome, annotation, and transcriptome datasets to verify microRNA host gene, mature microRNA, and microRNA binding site conservation at key loci functional in opioid addiction. We reveal a complex evolutionary landscape in which certain directional regulatory edges of the microRNA–mRNA hub-and-spoke network lack pan-mammalian conservation.

Annotation of primate miRNAs by high throughput sequencing of small RNA libraries

Background

In addition to genome sequencing, accurate functional annotation of genomes is required in order to carry out comparative and evolutionary analyses between species. Among primates, the human genome is the most extensively annotated. Human miRNA gene annotation is based on multiple lines of evidence including evidence for expression as well as prediction of the characteristic hairpin structure. In contrast, most miRNA genes in non-human primates are annotated based on homology without any expression evidence. We have sequenced small-RNA libraries from chimpanzee, gorilla, orangutan and rhesus macaque from multiple individuals and tissues. Using patterns of miRNA expression in conjunction with a model of miRNA biogenesis we used these high-throughput sequencing data to identify novel miRNAs in non-human primates.

Results

We predicted 47 new miRNAs in chimpanzee, 240 in gorilla, 55 in orangutan and 47 in rhesus macaque. The algorithm we used was able to predict 64% of the previously known miRNAs in chimpanzee, 94% in gorilla, 61% in orangutan and 71% in rhesus macaque. We therefore added evidence for expression in between one and five tissues to miRNAs that were previously annotated based only on homology to human miRNAs. We increased from 60 to 175 the number miRNAs that are located in orthologous regions in humans and the four non-human primate species studied here.

Conclusions

In this study we provide expression evidence for homology-based annotated miRNAs and predict de novo miRNAs in four non-human primate species. We increased the number of annotated miRNA genes and provided evidence for their expression in four non-human primates. Similar approaches using different individuals and tissues would improve annotation in non-human primates and allow for further comparative studies in the future.