Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Resource
Transparent Process
Open Access

Small RNA-Seq reveals novel miRNAs shaping the transcriptomic identity of rat brain structures

Anaïs Soula, Mélissa Valere, María-José López-González, Vicky Ury-Thiery, Alexis Groppi, View ORCID ProfileMarc Landry, Macha Nikolski, View ORCID ProfileAlexandre Favereaux  Correspondence email
Anaïs Soula
1University of Bordeaux, Bordeaux, France
2Centre Nationale de la Recherche Scientifique (CNRS), Unité Mixte de Recherche 5297, Interdisciplinary Institute of Neuroscience, Bordeaux, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mélissa Valere
1University of Bordeaux, Bordeaux, France
2Centre Nationale de la Recherche Scientifique (CNRS), Unité Mixte de Recherche 5297, Interdisciplinary Institute of Neuroscience, Bordeaux, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
María-José López-González
1University of Bordeaux, Bordeaux, France
2Centre Nationale de la Recherche Scientifique (CNRS), Unité Mixte de Recherche 5297, Interdisciplinary Institute of Neuroscience, Bordeaux, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Vicky Ury-Thiery
1University of Bordeaux, Bordeaux, France
2Centre Nationale de la Recherche Scientifique (CNRS), Unité Mixte de Recherche 5297, Interdisciplinary Institute of Neuroscience, Bordeaux, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexis Groppi
3Centre de Bioinformatique de Bordeaux, University of Bordeaux, Bordeaux, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Marc Landry
1University of Bordeaux, Bordeaux, France
2Centre Nationale de la Recherche Scientifique (CNRS), Unité Mixte de Recherche 5297, Interdisciplinary Institute of Neuroscience, Bordeaux, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Marc Landry
Macha Nikolski
3Centre de Bioinformatique de Bordeaux, University of Bordeaux, Bordeaux, France
4CNRS/Laboratoire Bordelais de Recherche en Informatique, University of Bordeaux, Talence, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexandre Favereaux
1University of Bordeaux, Bordeaux, France
2Centre Nationale de la Recherche Scientifique (CNRS), Unité Mixte de Recherche 5297, Interdisciplinary Institute of Neuroscience, Bordeaux, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Alexandre Favereaux
  • For correspondence: alexandre.favereaux@u-bordeaux.fr
Published 30 October 2018. DOI: 10.26508/lsa.201800018
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

In the central nervous system (CNS), miRNAs are involved in key functions, such as neurogenesis and synaptic plasticity. Moreover, they are essential to define specific transcriptomes in tissues and cells. However, few studies were performed to determine the miRNome of the different structures of the rat CNS, although a major model in neuroscience. Here, we determined by small RNA-Seq, the miRNome of the olfactory bulb, the hippocampus, the cortex, the striatum, and the spinal cord and showed the expression of 365 known miRNAs and 90 novel miRNAs. Differential expression analysis showed that several miRNAs were specifically enriched/depleted in these CNS structures. Transcriptome analysis by mRNA-Seq and correlation based on miRNA target predictions suggest that the specifically enriched/depleted miRNAs have a strong impact on the transcriptomic identity of the CNS structures. Altogether, these results suggest the critical role played by these enriched/depleted miRNAs, in particular the novel miRNAs, in the functional identities of CNS structures.

Introduction

miRNAs comprise one of the most abundant classes of gene regulatory molecules in the organism. They are small (22 nt) endogenous noncoding RNAs operating a negative translational regulation on mRNAs (Bartel, 2004). miRNAs are involved in the regulation of most signaling pathways, both in physiology and in pathology, in a broad range of organisms from viruses to animals (Carrington & Ambros, 2003). Their biogenesis mostly starts by the transcription of the miRNA gene by RNA pol II, generating an imperfect stem-loop structure, called the pri-miRNA (Bartel, 2004; Cai et al, 2004; Lee et al, 2004). The pri-miRNA is further processed several times and transported into the cytoplasm. Briefly, the enzyme Drosha cleaves the single-stranded primary miRNA transcripts to produce a stem-loop secondary structure: the pre-miRNA. Then, the pre-miRNA is cleaved by the enzyme Dicer to produce a mature miRNA duplex (Bernstein et al, 2001). Finally, one of the two strands of the mature miRNA duplex, called the guide strand, is incorporated into a proteic complex called the RNA-induced silencing complex (RISC) complex. From this point, this miRNA-RISC complex, referred to as the miRISC, will exert a translation inhibition on target mRNA (Lee et al, 2002). The target recognition by the miRNA is mainly based on a Watson–Crick pairing between a specific region of the miRNA called the seed region (nucleotides two to eight) and a complementary sequence on the 3′UTR of the target mRNA (Bartel, 2009; Wang, 2014). In addition to the seed region binding, the remainder of the miRNA can be engaged in an imperfect or perfect base pairing with the target mRNA, resulting respectively in translational inhibition and/or degradation of the target mRNAs. Eventually, both interaction modes induce a decrease in target mRNA translation into protein. Because of the imperfect binding between miRNAs and target mRNAs, miRNAs are able to precisely regulate the expression of dozens, if not hundreds, of genes. Indeed, a single miRNA can bind to numerous different mRNAs, whereas a single mRNA can present binding sites for multiple miRNA species (Friedman et al, 2009). Consequently, miRNAs are considered key players in the regulation of cellular gene expression.

The importance of miRNAs in the regulation of the main brain functions is now well established. Thus, miRNAs are involved in many important physiological processes, such as the development and the maturation of the nervous system (Díaz et al, 2014), synaptic plasticity, learning, and memory (Bredy et al, 2011; Aksoy-Aksel et al, 2014; Letellier et al, 2014; Rajman et al, 2017). In addition, miRNAs are involved in numerous neurodegenerative diseases, such as Alzheimer (Sarkar et al, 2016), Parkinson, Huntington, and amyotrophic lateral sclerosis (Nelson et al, 2008; Maciotta et al, 2013; Tan et al, 2015).

For example, the role of two brain-specific miRNAs, miR-134 and miR-124, has been well described. miR-134 regulates the size of dendritic spines in rat hippocampal neurons by the translational repression of the Lim-domain-containing protein kinase 1 (LimK1), an activator of actin polymerization (Schratt et al, 2006; Bernard, 2007). Thus, miR-134 regulates the morphology of most excitatory synapses with an impact on synaptic plasticity–dependent functions such as long-term potentiation and memory formation (Gao et al, 2010). miR-124, another brain-specific miRNA, plays a critical role in the transition of progenitor neuronal cells to adult neurons by inhibiting signaling pathways of non-neuronal genes, thereby facilitating the expression of the neuronal transcriptomic identity (Conaco et al, 2006; Makeyev et al, 2007).

Interestingly, miRNAs show different roles depending on the region in which they are expressed. For example, miRNAs are particularly studied in the striatum in the context of the addiction to drugs or alcohol (Bali & Kenny, 2013) because this region serves as a key area for habit learning, one of the major mechanisms for addiction behavior (Volkow et al, 2006). In the spinal cord, miRNAs are mainly studied in a context of spinal cord injury and pain (Favereaux et al, 2011; Elramah et al, 2014; López-González et al, 2017). In the olfactory bulb, in which the regeneration of neurons occurs frequently, the role of miRNAs in neurogenesis mechanisms has been particularly studied. It has been demonstrated that miR-132 plays a critical role to form the basis of a structural plasticity program (Pathania et al, 2012). Another work shows that miR-124 target genes are involved in morphogenesis and neuron differentiation in the olfactory bulb (Li & Ling, 2017).

However, analysis of miRBase (Griffiths-Jones, 2006) (http://www.mirbase.org), the most common database referencing miRNAs, shows a large difference between the number of miRNAs referenced in the mouse (Mus musculus, 1,915 mature miRNAs and 1193 precursors), the human (Homo sapiens, 2,558 mature miRNAs and 1,881 precursors), and in the rat (Rattus norvegicus, 765 mature miRNAs and 495 precursors). Thus, it seems that a large percentage of the miRNAs expressed in the rat are not yet uncovered. Rat is a gold standard model in neuroscience, and it is thus of prime importance to get an extensive knowledge of miRNA expression in the central nervous system (CNS).

In this study, we used the small RNA-Seq technology to perform an exhaustive analysis of the miRNome of the rat CNS. Even if previous studies explored the miRNome of the rat brain, either they focused their attention on one specific structure, such as the hippocampus (Shinohara et al, 2011), and the hypothalamus (Amar et al, 2012), or they analyzed the miRNome of a structure over different developmental stages (Yao et al, 2012; Yin et al, 2015). A more extensive study has been performed on five different regions of the brain, but unfortunately using the microarray technology, which limits the number of miRNA quantified and occludes the discovery of novel miRNAs (Olsen et al, 2009). To our knowledge, an extensive analysis of the miRNome of several structures of rat CNS has never been performed.

We focused our attention on the olfactory bulb, the cortex, the hippocampus, the striatum, and the spinal cord and found that of the 495 miRNAs referenced in miRBase, 365 are expressed in these five CNS structures. In addition, we discovered 90 novel miRNAs, some of them having orthologous sequences in the mouse or human that play a crucial role in the regulation of the functions of neurons. We also showed that each CNS structure expresses a particular miRNome, and interestingly, the novel miRNAs seem to play a key role in defining structure-specific miRNomes. In addition, we quantified the expression of all mRNAs using mRNA-Seq on the very same samples used for miRNA quantification. Then, we correlated miRNA and mRNA expression data using target predictions and selected the top 20 regulated mRNAs for further analysis. Gene Ontology (GO) (Ashburner et al, 2000; The Gene Ontology Consortium, 2017) analysis revealed that most of the regulated targets were already known for their role in the development or the function of the nervous system.

Results

Small RNA sequencing of five structures of the CNS

We performed small RNA-Seq to reveal the miRNome of five structures of the CNS: the olfactory bulb, the cortex, the hippocampus, the striatum, and the spinal cord. To perform a statistical analysis of the miRNome, we dissected these structures in triplicate from three male adult Wistar rats. Thus, we synthesized 15 independent libraries that were finally pooled for sequencing. RNA sequencing performed well, and we obtained a total of 454,577,868 reads from which 1,150,734 (0.3% of total) were filtered out because they were too short (size < 15 nt). To assess the quality of the reads, we used FastQC software (version 0.11.5, default parameters), and all sequences were flagged as good quality (Table S1).

Table S1 Small-RNA Seq reads quality assessment with FastQC software.

The CNS expresses known and novel miRNAs

We first checked the number of miRNAs known and referenced in miRBase that were expressed among the five CNS structures. From the 495 miRNAs referenced in rat miRBase, 365 miRNAs were expressed in our samples, corresponding to 73.7% of the miRBase-referenced miRNAs (Fig 1A).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. Pre-miRNAs expressed in the CNS structures as revealed by RNA-Seq.

(A) Pie chart of miRBase-referenced pre-miRNAs, 365 of them are expressed in the CNS. (B) Predictive structures of novel miRNAs were calculated with the RNAfold program. Colors represent the probability of matching between the bases. The predictive structures of novel miRNAs are similar to those of known miRNAs. (C) Histogram of the number of known and novel precursors of miRNAs expressed in each structure. (D) Venn diagramm of the repartition of pre-miRNAs in the CNS structures considering known and novel pre-miRNAs.

To discover novel miRNAs expressed in the brain, we used the miRPro software (Shi et al, 2015). Based on miRDeep2, miRPro analyzes the sequences obtained from small RNA-Seq experiments seeking for novel miRNAs using the following pipeline. First, miRPro aligns sequences of potential novel mature miRNAs on the genome. Then, the algorithm excises the corresponding potential pre-miRNAs from the genomic sequence. Finally, miRPro determines the most relevant pre-miRNA, in particular by the analysis of pre-miRNA folding with the RNAfold tool (Hofacker, 2003). miRPro analysis of the reads obtained from the five CNS structures leaded to the prediction of 8,416 novel pre-miRNAs (Table S2). To discard potential false-predictive novel pre-miRNAs, we refined the analysis by only conserving those expressed at least in two replicated samples from the same CNS structure and with a minimum expression of 1 count per million (cpm). Finally, we identified 107 novel pre-miRNAs (Table S3), which after miRNA processing gave 90 novel mature miRNAs (Table S4), not yet referenced in rat miRBase. The number of pre-miRNAs was higher than the number of mature miRNAs because 8 novel pre-miRNAs were duplicated in the genome. For example, novel-pre-miR-19 was present on both chromosomes 3 and 20, novel-miR-8 was present in 3 copies on chromosome 1, and novel-miR-2 was present on chromosomes 7, 12, 15, and 18 and twice on chromosome 1 (Table 1).

Table S2 miRPro prediction of novel pre-miRNAs from the rat central nervous system.

Table 3 Final list of novel pre-miRNAs from the rat central nervous system.

Table S4 List of novel mature miRNAs from the rat central nervous system.

View this table:
  • View inline
  • View popup
Table 1.

Chromosomic organization of the novel miRNAs.

The length of the novel miRNAs ranged from 18 to 24 nt, with a mean length of 20.59 nt, which is in line with the conventional characteristics of miRNAs. To assess whether these novel miRNAs display the canonical stem-loop structure of miRNAs, we checked the predictive folding of the novel pre-miRNA sequences. Fig 1B shows, as an example, the secondary structure of three novel pre-miRNAs: rno-mir-novel-42, the most expressed in our samples; rno-mir-novel-35, an miRNA with an expression level corresponding to the average expression of all novel miRNAs; and finally rno-mir-novel-23, which had the lowest expression. All these novel miRNAs display a stem-loop structure similar to the one of prototypical miRNAs (e.g., rno-mir-124-1, Fig S1). Thus, the predictive secondary structure of these novel pre-miRNAs is compatible with functional miRNAs.

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Predictive structure of rno-miR-124-1.

Predictive structure of rno-miR-124-1 has been performed with the RNAfold program. Colors represent the probability of matching between the bases.

Pre-miRNA expression in the five CNS structures

We compared pre-miRNA expression across the five CNS structures. Globally, all the studied structures expressed more or approximately the same amount of known pre-miRNA species (from 319 for the striatum and up to 342 for the spinal cord, Fig 1C). In addition, concerning novel pre-miRNAs, four of the five CNS structures under study expressed approximately the same amount of novel pre-miRNAs (49 for the olfactory bulb, 55 for the cortex, 53 for the hippocampus, and 60 for the striatum). In contrast, the spinal cord expressed a higher number of novel pre-miRNAs compared with the other structures, with 68 novel pre-miRNAs, which is more than the half of the total number of novel pre-miRNAs discovered.

Then, we analyzed the repertoire of pre-miRNAs detected in each CNS structure. We used a Venn diagramm representation in which each CNS structure is represented by a color: the cortex in red, the hippocampus in brown, the striatum in yellow, the spinal cord in green, and the olfactory bulb in blue. This representation (Fig 1D) showed that most of the miRNAs, known or novel, were detected in all the five CNS structures studied (319 of 472, 67.58%). However, some pre-miRNAs were detected only in one structure (9, 9, 7, 13, and 33 for the olfactory bulb, the cortex, the hippocampus, the striatum, and the spinal cord, respectively).

Chromosomal repartition of pre-miRNAs

To try to understand why these new miRNAs were not discovered before and to look for possible miRNA clusters, we analyzed the localization of the genes coding for the miRNAs along the 22 chromosomes of the rat genome.

First, we studied the strand of the DNA molecule on which miRNAs were located. Globally, there were more miRNAs, known and novel, along the positive strand compared with the negative strand (Fig 2A). In addition, the proportion of miRNAs located on the positive strand is significantly higher for novel miRNAs than for known miRNAs (Fisher exact test, P < 0.05). However, the biological meaning of this bias toward the positive strand for miRNA genes is still unknown.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. Chromosomic repartition of miRNAs.

(A) Pie charts of the repartition of pre-miRNAs per DNA strand. There are more novel miRNAs on the positive strand than for miRBase-referenced miRNAs, Fisher exact test, P < 0.05. (B) Histogram of the repartition of known and novel pre-miRNAs along the 22 chromosomes of the rat genome. (C) Correlative analysis of the number of miRNAs in function of the length of chromosomes, R2 = 0.3814, P = 0.0022. (D) Pie chart represents genome mapping of the novel pre-miRNAs: most of the new miRNAs are intergenic (57.94%) or intronic (28.04%). Genome browser view of an intergenic pre-miRNA (novel-miR-42, E) and an intronic pre-miRNA (novel-miR-13, F).

The chromosomal repartition of the miRNAs expressed in the five CNS structures combining both miRBase-referenced and novel miRNAs (Fig 2B) was similar to the one of miRBase-referenced miRNAs only (Fig S2). Indeed, chromosomes 1, 6, and X were the largest contributors to miRNAs expression in the CNS (49, 65, and 46, respectively, Fig 2B and Table S5) and chromosomes 1, 7, and 10 contained the highest number of novel miRNAs (12, 8, and 15 respectively, Fig 2B and Table S5).

Table S5 Chromosomal repartition of pre-miRNAs.

Figure S2. Chromosomal repartition of miRBase-referenced miRNAs expressed in the five CNS structures.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Chromosomal repartition of miRBase-referenced miRNAs expressed in the five CNS structures.

Considering rat miRNAs already referenced in miRBase, chromosomes 1, 6, and X contain the highest number of miRNAs (64, 73, and 59, respectively). Analysis of the expression of the miRBase-referenced miRNAs in the five CNS structures studied shows that chromosomes 1, 6, and X host the highest number of miRNAs (37, 61, and 45, respectively).

Linear regression and statistical analysis showed a moderate correlation between the size of the chromosomes and the number of the miRNAs genes (Fig 2C, R2 = 0.3814, P = 0.0022). To go further on this analysis, we looked for miRNAs clusters. Because criteria for miRNA clusters are versatile in the literature (Altuvia et al, 2005), we selected three different ranges for miRNA clusters. We considered two miRNA genes as clustered if they were located within a 1-kb, 2-kb, or 5-kb region of the genomic DNA. According to these criteria, none of the novel miRNA genes were clustered to other known or novel miRNA gene (Table S6).

Table S6 miRNA cluster analysis.

Genome mapping analysis of the novel pre-miRNAs revealed that most of them were located in intergenic regions 62 (57.94%, Fig 2D and E) or introns 30 (28.04%, Fig 2D and F); the remaining were located in the promoter transcription start site: 10 (9.35%); transcription termination site: 10 (9.35%); exon: 1 (0.93%); or 3′UTR: 1 (0.93%).

Novel mature miRNAs

Maturation of the pre-miRNA produces two miRNAs corresponding to the 5′ and 3′ ends of the stem of the hairpin structure. One of the miRNAs, called the “passenger” miRNA, is degraded, whereas the other one, the “guide” or “mature” miRNA, is loaded into the miRISC complex to target mRNAs. To decipher the role of miRNAs in the transcriptomic regulation of CNS, it is thus of prime importance to analyze the expression of mature miRNAs.

First, we analyzed the nucleotidic composition of the novel miRNAs. We found that U is the most representative nucleotide in the 5′ position of the novel miRNA sequences, followed by C and A (Fig 3A).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. Analysis of the sequence of novel miRNAs.

(A) The nucleotidic composition of all novel miRNAs was compiled, and the probability of each nucleotide at each position is plotted; U is the predominant nucleotide at the 5′ position. (B) Proof of concept for the “functional” ortholog miRNA hypothesis, mmu-miR-676-3p, and novel-rno-miR-21-5p share the same seed region on the LCE2D 3′UTR. Luciferase experiment proves that both miRNAs are able to interact with LCE2D 3′UTR and mediate luciferase translation inhibition. Cel-miR-67, a miRNA from C. elegans known to have no target in mammals, is used as control. Data shown are mean ± SD; numbers within the bars indicate biological replicates, one-way ANOVA, followed by the Dunnett post hoc test, ***P < 0.001, *P < 0.05.

Then, we wondered if, of the 90 novel miRNAs we discovered, some were orthologs with a mouse and/or a human miRNA already referenced in miRBase. Using Blast software (Altschul et al, 1990), we defined as ortholog a miRNA sequence with a minimum of 95% similarities with a miRNA from another species. We found 12 novel rat miRNAs that were orthologs of mouse and/or human miRNAs (Table 2). More particularly, there were five miRNAs that were orthologs of mouse miRNAs only, whereas the other seven miRNAs were orthologs of both human and mouse miRNAs. For example, novel-miR-7-5p and novel-miR-8-3p were canonical orthologs of mouse miR-344-5p and miR-344h-3p, respectively, whereas novel-miR-14-3p was a canonical ortholog of both mouse and human miR-330-5p (Table 2).

View this table:
  • View inline
  • View popup
Table 2.

Some novel miRNAs have orthologous sequences in human and/or mouse genomes.

miRNA targets are commonly predicted on the basis of a perfect complementarity between the seed region of the miRNA and the 3′UTR of the target mRNA. Thus, miRNAs sharing the same seed region often share the same targets and are grouped as a family of miRNAs (Ambros et al, 2003). Consequently, in addition to canonical orthologs, we searched for what we called “functional orthologs,” which are miRNAs sharing, at least, the same seed region with an miRNA from another species.

We found in total 23 “functional orthologs” either of human or mouse miRNAs: 8 “functional orthologs” of mouse miRNAs, 12 “functional orthologs” of human miRNAs, and three novel rat miRNAs sharing the seed region with both mouse and human miRNAs (Table 3). For example, novel-miR-78-5p was a “functional ortholog” of mouse miR-3070-5p, novel-miR-73-3p was a “functional ortholog” of human miR-765, and novel-miR-84-5p was a “functional ortholog” of both human and mouse miR-134-3p. To evaluate the ability of these “functional orthologs” to repress mRNA expression the same way as miRNAs sharing the same seed region, we designed a proof-of-concept experiment based on a luciferase assay. Novel-rno-miR-21-5p was predicted to be a “functional ortholog” of mouse mmu-miR-676-3p and thus may regulate the same target genes. We searched TargetScan for the predicted targets of mmu-miR-676-3p, and from the top five predicted targets, we selected LCE2D. LCE2D stands for late cornified envelope protein 2D, a protein involved in the development of skin. We selected LCE2D because the interaction with mmu-miR-676-3p relied only on the hybridization of the seed region, and interestingly, novel-rno-miR-21-5p was also predicted to interact with LCE2D only via the seed region. We cloned the 3′UTR sequence of LCE2D in a luciferase reporter plasmid and tested in COS cells the ability of novel-rno-miR-21-5p to regulate LCE2D. As a control miRNA, we used Cel-miR-67, a miRNA from C. elegans known to have no target in mammals. Fig 3B shows that mmu-miR-676-3p inhibited the expression of the luciferase reporter by 15% compared with control condition, and more interestingly, novel-rno-miR-21-5p inhibited the reporter gene by 22%. The level of regulation was moderate, but this is a common feature of miRNAs that are known as fine-tuner of gene expression. The key point is that both miRNAs were able to regulate the same target, and this can be considered a clue that “functional orthologs” may regulate the same targets.

View this table:
  • View inline
  • View popup
Table 3.

Novel miRNAs with their corresponding functional orthologous sequences in the mouse and/or human.

In conclusion, among the 90 novel mature miRNAs discovered in the rat CNS by small RNA-Seq, only 35 were orthologs, either classical or “functional” orthologs, of mouse or human miRNAs. Consequently, most novel miRNAs that we discovered are totally new, with no orthologs at least in the mouse or human.

Clustering of the CNS structures in regard to their miRNome

To test whether miRNA expression discriminated the different brain structures, we performed hierarchical clustering analysis of the miRNome of the five CNS structures tested. To perform this analysis, we used the count per million value for each miRNA (Table S7), and we represented the results as a heatmap.

Table S7 Expression data for all miRNAs for all biological replicates (count per million).

miRbase-referenced miRNAs

First, we analyzed the expression of miRBase-referenced miRNAs in the five CNS structures studied. Hierarchical clustering based on the expression of known miRNAs, followed by heatmap representation, showed that replicates of the same CNS structure are almost all grouped together (Fig 4A). Thus, except for one cortex replicate that was not clustered with the two others biological replicates, all other biological replicates were grouped together. Particularly, we can notice that most of the miRBase-referenced miRNAs were expressed approximately at the same level across the different samples, and this is in line with the Venn diagramm analysis (Fig 1D). However, specific miRNAs were differentially expressed, giving a miRNA signature for each CNS structure. For instance, the high expression level of miR-96-5p, -200a-3p, -200b-3p, and -200c-3p in the olfactory bulb discriminated this structure from the others.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. Hierarchical clustering of the miRNome of the different CNS structures.

(A, B) Heatmaps representing the expression of each miRNA from all the studied structures. miRNAs are represented at the bottom, biological replicates of the different CNS structures are represented on the right, and statistical dendrogram of clusterization of the samples is represented on the left. Colors represent the level of miRNA expression (log2 of count per million); red: high expression; green: low expression. (A) Considering miRBase-referenced miRNAs, all biological replicates of the same structure are grouped together, except cortex replicate #3. (B) Hierarchical clustering of novel miRNAs shows a perfect clustering of the biological replicates from the same structure. Black arrowheads indicate miRNAs those expression are strongly different in a structure compared with all other samples.

Novel miRNAs

Then, we did the same analysis on the newly discovered miRNAs. This second heatmap showed a perfect clustering of the three replicates from the same structure (Fig 4B). This perfect clustering may result from the enrichment of some structures for specific miRNAs. For example, novel-miR-47-5p had a low expression in all structures except in the olfactory bulb. The other way around, the low expression of some miRNAs may enable the clustering of a given CNS structure. For instance, novel-miR-28-3p was well expressed in all structures but the cortex. These results suggest that the novel miRNAs discovered may be important to the specification of the CNS structures.

miRNAs are differentially expressed in CNS structures

Hierarchical clustering analysis suggests that we could use the expression of miRNAs as a criterion for CNS structures classification. To identify which miRNAs were specifying CNS structures, we performed a statistical analysis of miRNA expression data with the DESeq2 algorithm (Love et al, 2014) (default parameters). Thus, we were able to determine miRNAs that were underexpressed or overexpressed in a CNS structure compared with all the others.

We found differentially expressed miRNAs in each structure but in variable amount. A complete list of all differentially expressed miRNAs, their fold changes to the other CNS structures, and statistical P-values can be found in Table S8. In particular, the olfactory bulb had the highest number of differentially expressed miRNAs with 44 overexpressed and 34 down-regulated (Fig 5A). In the spinal cord, there were 32 overexpressed miRNAs and 33 downregulated (Fig 5E). Four miRNAs were overexpressed in the cortex and seven down-regulated (Fig 5B). In the hippocampus, only two miRNAs were differentially expressed, one up-regulated and one down-regulated (Fig 5C). Finally, in the striatum, we found 20 miRNAs overexpressed and two underexpressed (Fig 5D).

Table S8 Differential expression analysis of miRNAs in each structure determined with the DESeq2 algorithm.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5. Differentially expressed miRNAs between the structures.

All the depicted miRNAs are statistically dysregulated in the structure considered compared with all other structures as assessed with the DESeq2 algorithm and Wald test (P < 0.05). The log2 of the fold change of over- (green) and down-regulated (red) miRNAs is indicated; miRNAs written in bold correspond to novel miRNAs. (A) The olfactory bulb shows the highest number of differentially regulated miRNAs: 44 are up-regulated and 34 are down-regulated. In the cortex (B), four miRNAs are significantly up-regulated and seven are down-regulated, whereas in the hippocampus (C), only one miRNA is overexpressed and one is underexpressed. (D) In the striatum, 20 and 2 miRNAs are up- and down-regulated, respectively. (E) In the spinal cord, 32 and 34 miRNAs are up- and down-regulated, respectively.

Interestingly, among the differentially expressed miRNAs, some were newly discovered miRNAs. For instance, in the olfactory bulb, novel miRNAs rno-novel-miR-47-5p and -56-5p were up-regulated, whereas rno-novel-miR-42-5p and -59-3p were down-regulated. In total, 9 of the 90 newly discovered miRNAs were differentially expressed in CNS structures. These results suggest that the expression of some miRNAs may regulate specifically gene expression in a given structure.

Structure-specific miRNAs and the regulation of their predicted target genes

The rationale was that to have a strong impact on gene expression and finally on the physiological function of the different CNS structures, specific miRNAs should be highly regulated in one structure compared with all other ones. In addition to the statistical criterion used to define differentially expressed miRNAs, we used a second criterion for specificity. We first searched for enriched miRNAs, and we selected as specifically enriched miRNAs, those with a log 2 of the fold change (log2FC) greater than four, meaning that the expression of the given miRNA was at least 16 times higher in this structure than in any other structures.

In addition to specifically enriched miRNAs, the specific depletion of miRNAs in a structure could have a strong impact on CNS physiology by unleashing gene expression of specific targets. Thus, in addition to specifically enriched miRNAs, we also looked for specifically depleted miRNAs. We selected as depleted miRNAs, miRNAs which had a log 2 of the fold change lower than −4, corresponding to, at least, a 16 times decrease in miRNA expression.

Enrichment/depletion analysis of the CNS miRNome showed that except for the striatum, all structures had specific miRNAs. In particular, there were many miRNAs enriched in the olfactory bulb (12) and in the spinal cord (6) but none in the cortex and in the hippocampus (Fig 6). Concerning down-regulated miRNAs, five miRNAs were depleted in the spinal cord but only one in the olfactory bulb, the cortex, and the hippocampus. Interestingly, the only specific miRNA of the cortex was a novel miRNAs, rno-novel-miR-28-3p. In addition, of the five specifically depleted miRNAs in the spinal cord, three were novel miRNAs (Fig 6).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6. Enrichment/depletion of miRNAs in specific structures.

To focus on the miRNAs that could have the most relevant role in transcriptome regulation, we defined as enriched miRNAs those with a fold change >16, and depleted miRNAs those with a fold change <−16. Except for the striatum, all structures express specifically enriched or depleted miRNAs. The olfactory bulb shows the highest number of specific miRNAs (13 miRNAs), the spinal cord expresses 11 specific miRNAs, the cortex and the hippocampus exhibit only one specifically depleted miRNA, and interestingly, the cortex-specific miRNA is a novel miRNA (novel miRNAs are written in bold).

These results suggest that these specific miRNAs should differentially modulate the transcriptome of these CNS structures. To test this hypothesis, we quantified the expression of all mRNAs using mRNA-Seq on the very same samples used for miRNA quantification. We first determined regulated mRNAs by performing a statistical analysis of mRNA expression data with the DESeq2 algorithm (Love et al, 2014) (Table S9, default parameters). Then, we correlated miRNA and mRNA expression data using target predictions as determined by TargetScan (default parameters) for known miRNAs (Agarwal et al, 2015) and miRDB tools for novel miRNAs. This correlation analysis revealed many miRNA-mRNA pairs regulated in the CNS structures studied except for the striatum (Table S10). Next, we assessed, for each structure, if predicted target genes of the specific miRNAs were more regulated than nontarget genes. To do so, we analyzed the cumulative frequency distribution of expression change of all miRNA target genes versus nontarget genes. Then, in an attempt to highlight the most biologically relevant targets for each structure, we selected the mRNAs displaying the most important regulation (either the 20 uppest or 20 lowest regulated mRNAs). From this list, we performed a GO term enrichment analysis using the Metascape tool (Tripathi et al, 2015). Concerning specifically enriched miRNAs in the olfactory bulb, the cumulative frequency distribution of expression changes of all predicted targets was significantly different from nontarget genes (Fig 7A). Although most of the predicted targets showed an expected down-regulation, some were up-regulated. However, gene-by-gene analysis showed a statistical difference in 12 genes that were all down-regulated (Fig 7B). Interestingly, some regulated mRNAs were predicted to be targeted by multiple structure-specific miRNAs. For instance, Col25a1 (log2FC = −2.18), a brain-specific membrane collagen, was predicted to be inhibited by miR-96-5p and miR-182. In addition, Caln1 (log2FC = −2.85), a protein with high similarity to the calcium-binding proteins of the calmodulin family, was predicted to be targeted by miR-183-5p and novel-miR-47-5p. Among the most regulated mRNAs were multiple members of the solute carrier family, such as Slc18a3 (log2FC = −4.31) and Slc6a5 (log2FC = −3.48), which transport acetylcholine and glycine, respectively, two important neurotransmitters in the brain. In the same structure, one miRNA was specifically depleted, miR-544-3p (log2FC = −18.6), and the cumulative frequency distribution of expression changes of all its predicted targets was significantly different from nontarget genes (Fig 7C). From the 40 statistically up-regulated targets of miR-544-3p (Table S10), the 20 most up-regulated were shown (Fig 7D). For instance, the expression of Cdhr1 (log2FC = 6.02), a cadherin-related protein, the transcription factors Sp8 (log2FC = 4.97) and Barhl2 (log2FC = 3.79), the neurogenesis marker Doublecortin (Dcx, log2FC = 3.55), and Wnt5a (log2FC = 3.39) were highly up-regulated in the olfactory bulb compared with all other structures studied. GO term enrichment suggests that the regulated miRNA-mRNA pairs could have a role in the specification of this brain structure (Fig 7E). Indeed terms such as “axon extension”, “regulation of anatomical structure size,” and “chemical synaptic transmission” were statistically enriched.

Table S9 Differential expression analysis of mRNAs in each structure determined with the DESeq2 algorithm.

Table S10 Correlation of miRNA and mRNA expression data according to miRNA target prediction from Tagetscan and miRDB.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7. Olfactory bulb–specific miRNAs and their regulated targets.

(A) Concerning specifically enriched miRNAs, cumulative frequency distribution of expression changes of all predicted targets is significantly different from nontarget genes (blue line versus black line, Kolmogorov–Smirnov test). Distribution of expression changes of statistically regulated target genes (red) shows a down-regulation compared with nontarget genes. (B) Individual gene statistical analysis revealed that 12 target genes of the enriched miRNAs are statistically down-regulated (DESeq2 algorithm and Wald test, P < 0.05). Interestingly, some regulated genes are predicted to be the target of multiple miRNAs (gene names written in italics). Data shown are mean ± SD form three biological replicates. (C) Concerning specifically depleted miRNAs, cumulative frequency distribution of expression changes of all predicted targets is significantly different from nontarget genes (blue line versus black line, Kolmogorov–Smirnov test). Distribution of expression changes of statistically regulated target genes (red) shows an up-regulation compared with nontarget genes. (D) Individual gene statistical analysis revealed that 40 target genes of miR-544-3p are up-regulated (DESeq2 algorithm and Wald test, P < 0.05); only the 20 most regulated are depicted. Data shown are mean ± SD form three biological replicates. (E) GO term enrichment analysis of the selected targets.

In the cortex, only one miRNA, the newly discovered novel-miR-28-3p, was specifically depleted (log2FC = −7.55), and the cumulative frequency distribution of expression changes of all its predicted targets was significantly shifted to the right compared with nontarget genes, revealing a global up-regulation (Fig 8A). Twenty of novel-miR-28-3p predicted targets were statistically up-regulated (Fig 8B), for example, Kcns1 (log2FC = 5.27), a voltage-gated potassium channel subunit, Vip (log2FC = 2), a neurotransmitter and neuromodulator with neurotrophic properties, cerebellin 2 (Cbln2, log2FC = 1.96), a protein forming homohexamers at the presynaptic compartment, and Lynx1 (log2FC = 1.2), a protein associated with the cholinergic regulation. GO term analysis revealed neuroscience-familiar keywords such as “potassium ion transport” and “actin cytoskeleton organization” (Fig 8C).

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8. Cortex-specific miRNA and its regulated targets.

(A) Cumulative frequency distribution of expression changes of all predicted targets is significantly shifted to the right compared with nontarget genes, revealing a global up-regulation (blue line versus black line, Kolmogorov–Smirnov test). Distribution of expression changes of statistically regulated target genes (red) shows an up-regulation compared with nontarget genes. (B) Individual gene statistical analysis revealed that 20 target genes of novel-miR-28-3p are up-regulated (DESeq2 algorithm and Wald test, P < 0.05). Data shown are mean ± SD form three biological replicates. (C) GO term enrichment analysis of the selected targets.

In the hippocampus, the specific depletion of miR-3065-5p (log2FC = −5.85) was associated with the up-regulation of many of its predicted targets. Hence, the cumulative frequency distribution of expression changes of all its predicted targets was significantly shifted to the right compared with nontarget genes, revealing an overall up-regulation (Fig 9A). Statistical analysis showed that 70 target genes of miR-3065-5p were significantly up-regulated in the hippocampus compared with all other CNS structures (Table S10). Among these, many have a known role in neuronal function such as, the LIM homeobox 9 protein (Lhx9, log2FC = 4.53), the interleukin IL-16 (log2FC = 3.6), the gonadotrophin-releasing hormone receptor (Gnrhr, log2FC = 3.55), the protein kinase C γ (Prkcg, log2FC = 2.5), the protocadherin-19 (log2FC = 1.88), and the GluA1 AMPA receptor (Gria1, log2FC = 1.85), Fig 9B. In addition, these enriched mRNAs were associated with relevant GO terms such as “negative regulation of nervous system development” (Fig 9C).

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9. Hippocampus-specific miRNA and its regulated targets.

(A) Cumulative frequency distribution of expression changes of all predicted targets is significantly shifted to the right compared with nontarget genes, revealing a global up-regulation (blue line versus black line, Kolmogorov–Smirnov test). Distribution of expression changes of statistically regulated target genes (red) shows an up-regulation compared with nontarget genes. (B) Individual gene statistical analysis revealed that 70 target genes of novel-miR-28-3p are up-regulated (DESeq2 algorithm and Wald test, P < 0.05); only the 20 most regulated are depicted. Data shown are mean ± SD form three biological replicates. (C) GO term enrichment analysis of the selected targets.

Finally, concerning specifically enriched miRNAs in the spinal cord, the cumulative frequency distribution of expression changes of all predicted targets was significantly different from nontarget genes (Fig 10A). Although most of the predicted targets showed an expected down-regulation, some were up-regulated. Among the predicted targets, 71 were statistically down-regulated (Table S10). Fig 10B shows the 20 most down-regulated mRNAs, such as the transcription factor Bcl11b (also known as Ctip2, log2FC = −4.08) targeted by both miR-10b-5p and miR-615, the vesicular glutamate transporter 1 (Slc17a7, log2FC = −5.44), the transcription factor subunit FosB (log2FC = −4.67), the activity-regulated cytoskeleton-associated protein (Arc, log2FC = −4.25), or the synaptopodin protein (Synpo, log2FC = −3.83). Concerning specifically depleted miRNAs in the spinal cord, the cumulative frequency distribution of expression changes of all predicted targets was significantly different from nontarget genes (Fig 10C). Although most of the predicted targets showed an expected up-regulation, some were down-regulated. Among the predicted targets, 91 were statistically up-regulated (Table S10). Hence, many predicted targets of miR-344 g (log2FC = −8.11) were up-regulated (Fig 10D) such as the transcription factors paired box gene 2 (Pax2, log2FC = 7.0), Lmx1b (log2FC = 4.69), and Lbx1 (log2FC = 4.62). The glycine neurotransmitter transporter Slc6a5 (log2FC = 6.51) and the apelin receptor (Aplnr, log2FC = 4.35) were also up-regulated. miR-551b-5p was another miRNAs strongly depleted (log2FC = −7.74) in the spinal cord, and its predicted targets were largely up-regulated, for instance, the transcription factors Hoxb3 (log2FC = 6.91), Pax8 (log2FC = 6.67), Nkx6-1 (log2FC = 4.65), and Lhx4 (log2FC = 4.31), polycystic kidney disease 2–like 1 protein (Pkd2l1 also known as Trpp3, log2FC = 5.51), a cationic channel, and the glucagon-like peptide 1 receptor (Glp1r, log2FC = 5.06). GO term analysis highlighted many keywords relevant in the neuroscience and neurodevelopment field such as “pattern specification process”, “spinal cord development” and “chemical synaptic transmission” (Fig 10E).

Figure 10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 10. Spinal cord–specific miRNAs and their regulated targets.

(A) Concerning specifically enriched miRNAs, cumulative frequency distribution of expression changes of all predicted targets is significantly different from nontarget genes (blue line versus black line, Kolmogorov–Smirnov test). Distribution of expression changes of statistically regulated target genes (red) shows a down-regulation compared with nontarget genes. (B) Individual gene statistical analysis revealed that 71 target genes of the enriched miRNAs are statistically down-regulated (DESeq2 algorithm and Wald test, P < 0.05); only the 20 most regulated are depicted. Interestingly, some regulated genes are predicted to be the target of multiple miRNAs (gene names written in italics). Data shown are mean ± SD form three biological replicates. (C) Concerning specifically depleted miRNAs, cumulative frequency distribution of expression changes of all predicted targets is significantly different from nontarget genes (blue line versus black line, Kolmogorov–Smirnov test). Distribution of expression changes of statistically regulated target genes (red) shows an up-regulation compared with nontarget genes. (D) Individual gene statistical analysis revealed that 91 target genes of the depleted miRNAs are up-regulated (DESeq2 algorithm and Wald test, P < 0.05); only the 20 most regulated are depicted. Interestingly, some regulated genes are predicted to be the target of multiple miRNAs (gene names written in italics). Data shown are mean ± SD form three biological replicates. (E) GO term enrichment analysis of the selected targets.

Discussion

The aim of this study was to determine the miRNome of five major structures of the CNS, the cortex, the hippocampus, the olfactory bulb, the striatum, and the spinal cord, and to evaluate the impact of enriched/depleted miRNAs on the transcriptome.

Three principal methods are used to measure the expression levels of miRNAs in a tissue: reverse transcription quantitative PCR (RT-qPCR) (Chen et al, 2005; Shi & Chiang, 2005), microarray hybridization (Yin et al, 2008), and small RNA-Seq (Hafner et al, 2008), all of which face unique challenges compared with their use in miRNA profiling. However, RNA-Seq offers competitive advantages, such as enabling discovery of new sequences variants and new sequences (Willenbrock et al, 2009). Consequently, we performed small RNA-Seq to obtain an exhaustive analysis of the miRNome of these five CNS structures, and we found a total of 90 novel miRNAs.

Previous studies used small RNA-Seq to analyze the miRNome of the rat brain. Unfortunately, all these studies analyzed miRNA expression in only one brain structure, occluding any attempt to identify structure-specific miRNAs. For instance, Shinohara et al (2011) focused their attention on the hippocampus and overall found a relatively modest number of novel miRNAs (17), which was restricted to two candidates when they applied a filter based on the abundance. None of the two miRNA candidates identified by Shinohara et al (2011) were found in our data. In addition, Amar et al (2012) analyzed the expression of miRNAs in the hypothalamus only, and the number of reads they obtained for each sample was not compatible with novel miRNA discovery. Yao et al analyzed the miRNome of the developing cerebral cortex from E10 to P28 (Yao et al, 2012) found 101 potential novel miRNAs. Finally, Yin et al used as starting material the whole brain with all structures mixed together and found 171 candidate novel miRNAs. Unfortunately, the novel miRNAs discovered in Yin's study were not yet implemented in miRBase. However, we aligned the sequences of their novel miRNAs against ours and found matches for five miRNAs. Thus, Yin's novel-miR-16, novel-miR-18, novel-miR-20, novel-miR-32, and novel-miR-33 match with ours novel-miR-79-3p, novel-miR-64-5p, novel-miR-38-5p, novel-miR-87-5p, and novel-miR-67-5p, respectively.

All the computational tools developed for identifying miRNAs from small RNA-Seq data sets suffer from high false-positive and false-negative rates and also of a lack of consistency across species (Li et al, 2012; Williamson et al, 2013; Kang & Friedländer, 2015). A strategy to limit the number of false-positive novel miRNAs and to focus on the most biologically relevant ones could be to apply a filter on the predicted novel miRNAs based on the abundance. Thus, we initially obtained more than 8000 novel miRNAs from our samples with the miRPro algorithm (based on miRDeep2), and we reduced this number to 90 by applying two filters: the number of reads (>1 cpm of total reads) and the presence of the novel miRNA in at least two independent samples. The number of novel miRNAs discovered by small RNA-Seq experiments is hard to predict because it depends on many parameters such as the diversity of the nature of samples (culture, tissues, and fluids), the number of reads for each sample, the quality of the library, and finally the bioinformatic tools to predict novel miRNA and the efficiency to discard false-positive novel miRNAs. In this regard, our results are in accordance with the literature. Indeed, using miRDeep2 as miRNA-prediction algorithm, recent studies found from 71 novel miRNAs in the mouse brain after an infection with Japanese encephalitis viruses (Li et al, 2016), up to 171 novel miRNAs for the Yin’s study (Yin et al, 2015).

In comparison with the miRNA database (miRBase) of the mouse (Mus musculus, 1,915 mature miRNAs and 1,193 precursors) and human (Homo sapiens, 2,558 mature miRNAs and 1,881 precursors), the content of the rat miRNA database is low (Rattus norvegicus, 765 mature miRNAs and 495 precursors). This may reflect the fact that the human and the mouse are two species intensively studied from a genetic point of view compared with the rat. Even if we implement our novel miRNAs in miRBase, the rat miRBase would still be more than half less referenced than the mouse one and three times less referenced than the human one. This may be due to the fact that we focused our study only on the CNS and that we may have discovered more miRNAs if we had made a pan-organism study. Nevertheless, our work significantly increases the knowledge on rat brain miRNA expression.

We searched orthologous sequences for our novel miRNAs in the mouse and human genomes because miRNAs that are well conserved across species often play crucial roles. Among the 90 novel miRNAs discovered, 12 have orthologous sequences in mouse and/or human organisms (we determined orthologs by a minimum of 95% of similarities between the sequences). Thus, most of these newly discovered miRNAs in the rat brain have no homologue in the mouse or human. Because mouse and human brains had already been analyzed for new miRNAs, it may suggest that these miRNAs are rat specific, supporting rat-specific gene regulations. Studies on human genome suggest that recently discovered miRNAs are mainly evolutionarily young (Friedländer et al, 2014). In addition, it has been recently shown that within a species, interstrain variation exists with functional significance on the targeted mRNAs (Trontti et al, 2018). However, we cannot exclude that the mouse and human orthologs of these miRNAs have not been found yet due to the bioinformatic tools used for miRNA prediction. For example, both novel-miR-7-5p and novel-miR-8-3p are orthologs of the mouse miR-344-5p and miR-344h-3p, respectively. Interestingly, a work showed by RT-qPCR and in situ hybridization that miR-344-3p is detected in cortical regions surrounding the ventricular system, choroid plexus, glomerular layer of the olfactory bulb, and granular cell layer of the cerebellar cortex (Liu et al, 2014). They showed that miR-344-3p is neural specific during mouse embryonic development and that it may play an important role in CNS morphogenesis. In our analysis, we found that novel-miR-7-5p and novel-miR-8-3p are specifically down-regulated in the spinal cord, but expressed in all the other structures, approximately at the same level. It should be interesting to define more specifically the spatial and temporal distribution of these novel-miRNAs to compare with mouse miR-344 expression and function.

Furthermore, among the 90 novel miRNAs discovered, we found 23 miRNAs that we called “functional orthologs”, meaning that they share the same seed region with other miRNAs. Because miRNA target prediction is based on the complementarity between the seed region of the miRNA and the 3′UTR of the target mRNA, we hypothesized that these “functional orthologs” would have the same targets. Indeed, miRNAs sharing the same seed region often share the same targets and are grouped as a family of miRNAs (Ambros et al, 2003). To test this hypothesis, we designed a proof-of-concept experiment based on a luciferase assay using the novel-rno-miR-21-5p that was predicted to be a “functional ortholog” of mouse mmu-miR-676-3p. We used LCE2D, a predicted target gene of mmu-miR-676-3p sharing with novel-rno-miR-21-5p the same interaction based on the seed region only. Luciferase experiment showed that both mmu-miR-676-3p and novel-rno-miR-21-5p were able to inhibit the expression of the luciferase-LCE2D construct. This result can be considered a clue that “functional orthologs” may regulate the same targets. However, it cannot be assumed from one experiment that they are functionally equivalent and further reporter gene assays have to be performed to confirm the targets of all “functional orthologs”.

The increasing amount of small RNA-Seq data provides a more complete view of the miRNome, with the discovery of novel canonical miRNAs, and many other miRNA-like small RNAs that originate from non-canonical miRNA genes. In addition, the increasing number of novel miRNAs discovered across the literature, sometimes repeatedly such as Yin's work and ours, increases the complexity of miRNA indexation. These arguments enforce to revise miRNA's nomenclature conventions to ease the comparison of miRNA research across species (Desvignes et al, 2015). Characterization of 5′-end sequence of miRNAs is important because this so-called seed region of miRNAs is critical in miRNA-target mRNA binding (Engels & Hutvagner, 2006). Analysis of nucleotide sequences in eukaryotic miRNAs showed a clear bias for U or A at the 5′ position (Frank et al, 2010). Here, we found that U was the most representative nucleotide in the 5′ position of our novel miRNA sequences, followed by C and A. Results from others studies on other organisms are consistent (Zhou et al, 2009; Wei et al, 2011; Ji et al, 2012; Pacholewska et al, 2016). This result reinforces the biological pertinence of our novel miRNAs.

It is well described that tissue- and cell-specific expression patterns of miRNAs are responsible for their particular biological function within a specific system. For example, the most famous brain-enriched miRNA, miR-124, defines brain-specific gene expression in human cells (Lim et al, 2005).

The rationale of our study was to find structure-specific miRNAs that could regulate a specific set of genes in a given CNS structure. We reasoned that to have a strong impact on gene expression and finally on the physiological function of the different CNS structures, specific miRNAs should be highly regulated in one structure compared with all other ones. Thus, we selected as specifically enriched miRNAs, those with a log 2 of the fold-change greater than 4, meaning that the expression of the given miRNA was at least 16 times higher in this structure than in any other structure. In addition to specifically enriched miRNAs, the specific depletion of miRNAs in a structure could have a strong impact on CNS physiology by unleashing gene expression of specific targets. Therefore, we selected as depleted miRNAs, miRNAs which had a log 2 of the fold change lower than −4, corresponding to, at least, a 16 times decrease in miRNA expression. This analysis revealed specific miRNAs for each structure except for the striatum. To evaluate the potential impact of these specific miRNAs, we measured the expression of their predicted target mRNAs. Thus, we quantified the expression of all mRNAs using mRNA-Seq on the very same samples used for miRNA quantification. After statistical analysis of the differentially regulated mRNAs, we correlated miRNA and mRNA expression data using target predictions as determined by TargetScan for known miRNAs and miRDB tools for novel miRNAs. To evaluate whether the predicted target genes were more regulated than nontarget genes, we analyzed the cumulative frequency distribution of expression change of all miRNA target genes versus nontarget genes. This analysis showed that the expression of predicted targets is significantly different from the one of nontargeted genes. Although in most cases the expression of the predicted targets was regulated as expected (opposite to miRNA expression because miRNAs act as inhibitors), on some occasion the predicted targets showed an unexpected regulation. This result is not surprising because miRNAs are not the only regulators of gene expression. Hence, other mechanisms, such as transcription factors, affect mRNA levels, and the net effect on mRNA expression is a combination of all these mechanisms. To try to focus on the most biologically relevant targets for each structure, we selected the mRNAs displaying the most important regulation (either the 20 uppest or 20 lowest regulated mRNAs). GO enrichment analysis of these top 20 mRNAs revealed that most of the regulated targets were genes known for their involvement in the development or the function of the nervous system. In addition, literature review showed that some of these targets specifically regulated in a CNS structure were already known to have a specific role in the given structure. Table 4 summarizes the correlation between our results and the literature.

View this table:
  • View inline
  • View popup
Table 4.

Relevance of the correlation between miRNA and mRNA expression.

We found 12 enriched miRNAs in the olfactory bulb, and among them the expression of all the members of the well-described miR-200 family (Senfter et al, 2016) composed of miR-200a, b, c, and miR-429. Recent studies demonstrate that the miR-200 family has an important role for the regulation of the proliferation and the differentiation of neuronal cells (Pandey et al, 2015; Beclin et al, 2016). More interestingly, a work incriminates the miR-200 family in the regulation of the olfactory neurogenesis (Choi et al, 2008), which is in accordance with our result showing an over-representation of all the members of this family of miRNAs in the olfactory bulb. In addition, we found that miR-544-3p was strongly down-regulated and among its predicted targets several genes are highly up-regulated. For instance, Sp8, a zinc-finger transcription factor known to regulate olfactory bulb interneuron development (Li et al, 2018) Bar H Like Homeobox 2 (Barhl2), another transcription factor, is known to play a key role in the development of region-specific differences in embryonic mouse forebrain through interaction with Pax6 (Parish et al, 2016) and with the diencephalon patterning by regulating Shh (Ding et al, 2017). In addition, Doublecortin (Dcx) is mandatory for proper establishment of the olfactory tract (Belvindrah et al, 2011), and Wnt5a is required for the development of the olfactory bulb (Zaghetto et al, 2007; Paina et al, 2011; Pino et al, 2011).

In the cortex, the novel-miR-28-3p is specifically depleted, and cortex-relevant predicted targets are highly up-regulated such as the vasoactive intestinal peptide (Vip). Vip-expressing interneurons have a key role in cortical circuit development, suggesting a possible contribution to pathophysiology in neurodevelopmental disorders (Batista-Brito et al, 2017). Cerebellins are secreted hexameric proteins that form tripartite complexes with the presynaptic cell-adhesion molecules neurexins and the postsynaptic glutamate receptor–related proteins GluD1 and GluD2. These tripartite complexes are believed to regulate synapses. In particular, cerebellin 2 (Cbln2) is specifically expressed in a subpopulation of excitatory cortical neurons (Seigneur & Südhof, 2017). Finally, Lynx1 is a cholinergic brake involved in the maintenance of the stability of cortical networks (Morishita et al, 2010).

In the hippocampus, miR-3065-5p was strongly down-regulated, and among its putative targets, many were up-regulated and related to hippocampal functions. Hence, the developmental regulatory gene Lhx9 seems to be involved in the development of the hippocampal subdivisions (Abellán et al, 2014). The neuronal interleukin-16 (Il16), a PDZ-containing protein, enriched in the hippocampus is responsible for the modulation of Kv4.2K+ currents, thereby regulating the intrinsic neuronal properties (Fenster et al, 2007). In addition, the gonadotrophin-releasing hormone receptor (Gnrhr) was up-regulated, and a recent study strongly suggests a role of specific Gnrhr activation in neuronal plasticity in the hippocampus (Schang et al, 2011). Plasticity mechanisms such as long-term potentiation are also supported by a regulated target of miR-3065-5p, namely the proteinase C γ (Prkcg) (Gärtner et al, 2006). Non-clustered protocadherins (PCDHs) are calcium-dependent adhesion molecules hypothesized to be involved in neuronal circuit formation and plasticity. Thus, the protocadherin-19 (Pcdh19), a potential target of miR-3065-5p, is regulated in the hippocampus in response to a plasticity paradigm (Kim et al, 2010), and mutations on Pcdh19 are responsible for epilepsy and mental retardation (Dibbens et al, 2008). Finally, the GluA1 AMPA receptor, a key element of homeostatic plasticity in the hippocampus (Sutton et al, 2006), was up-regulated. Interestingly, we previously showed that GluA1 fine-tuning involved miRNA regulations (Letellier et al, 2014).

Finally, in the spinal cord, we found both specifically enriched and depleted miRNAs. We found six miRNAs specifically enriched in the spinal cord. Two families are represented: miR-196 family, represented by miR-196a and b (Chen et al, 2011), and miR-10 family, represented by miR-10a and b (Tehler et al, 2011). Interestingly, some studies demonstrate the role of the miR-196 in the regulation of the motor neuron programs (Asli & Kessel, 2010). Furthermore, our results are in accordance with the literature: miR-10a and miR-10b that are specifically expressed in the spinal cord in our study are also found specific of this region in the mouse (Bak et al, 2008). One of the predicted targets of miR10b-5p was down-regulated and known for its function in the spinal cord. Hence, in loss-of-function experiments, it has been shown that Bcl11b plays a critical role in the development of axonal projections to the spinal cord in vivo (Arlotta et al, 2005). miR-615 was another enriched miRNA, and interestingly, among its regulated targets, many were associated with the pain-processing mechanisms in the spinal cord. Hence, expression of the vesicular glutamate transporter 1 (Vglut1 also known as Slc17a7) defines a subpopulation of neurons involves in pain processing (Landry et al, 2004) and nerve injury induces a decrease in Vglut1 expression in the dorsal horn of the spinal cord (Brumovsky et al, 2007). The expression of the fos-related fosB gene, sharing many properties with c-Fos, is increased during spinal cord development (Redemann-Fibi et al, 1991), and interestingly, fosB is involved in the transcriptomic cascade triggered during nociception (Herdegen et al, 1991). In addition, the expression of the Arc/Arg3.1 protein is specifically enhanced in response to nociceptive stimulations (Hossaini et al, 2010) through the activation of the intracellular mGluR5 receptor (Vincent et al, 2016). Finally, we previously demonstrated that synaptopodin (Synpo), a key element of the spine apparatus, was regulated in the spinal cord in an animal model of bone cancer pain (Elramah et al, 2017).

Concerning depleted miRNAs, down-regulation of miR-344 g and miR-551b-5p was associated with a strong up-regulation of many of their predicted targets involved in the development of the spinal cord; for instance, Pax gene members, Pax2 and Pax8 (Batista & Lewis, 2008; Larsson, 2017), Hoxb3 (Yau et al, 2002), the Nkx homeobox gene member Nkx6-1 (Sander et al, 2000; Vallstedt et al, 2001), and the LIM homeobox transcription factor 1−β (Lmx1b [Ding et al, 2004; Hilinski et al, 2016]). Again, in the regulated targets, several were involved in pain mechanisms such as the homeodomain factor Lbx1 (Gross et al, 2002; Cheng et al, 2017), the apelin receptor (Aplnr [Xiong et al, 2017]), or the glucagon-like peptide-1 receptor (Glp1r [Gong et al, 2014a, 2014b]). Finally, a predicted target of miR-551b-5p, the transient receptor potential channel Pkd2l1, was up-regulated. Interestingly, Pkd2l1 is a specific marker of spinal cerebrospinal fluid–contacting neurons (Djenoune et al, 2014), an evolutionarily conserved population of neurons regulating locomotion by relaying mechanical stimuli to spinal circuits (Böhm et al, 2016).

Although the link between the enriched/depleted miRNAs and the regulated mRNAs in this study is only correlative and based on target predictions, the GO term analysis and the relevance of the literature associated with these targets strongly suggest that the specific miRNAs may play a role in the transcriptomic specificity of the CNS structures. The miRNAs playing the most important role on CNS structure function might be the specifically depleted miRNAs. Thus, by unleashing gene expression of their target mRNAs, they may play an important role in the specification and the function of the specific structures of the rat CNS. In the future, it would be interesting to artificially modify the expression of these miRNAs during development to confirm their role in the specification of the structures of the CNS. In that perspective, our work will help further studies in dissecting biological functions of miRNAs in the CNS.

Materials and Methods

Animals

Three male Wistar rats aged 12 wk were dissected to collect CNS tissues. All experimental procedures followed the ethical guidelines of the University of Bordeaux's ethics committee. Animals were anesthetized with an intraperitoneal injection of pentobarbital and then perfused with aCSF (NaCl 130.5 mM; KCl 2.4 mM, MgSO4 1.3 mM; KH2PO4 1.2 mM; Hepes 1.25 mM; CaCl2 2.4 mM; Glucose 10 mM; NaHCO3 19.5 mM) to wash out blood from the nervous tissue and avoid contaminations with blood cells RNA.

Sample collection and RNA extraction

From each rat, five structures of the CNS (the olfactory bulb, a part of the cortex, the hippocampus, the striatum, and the dorsal spinal cord) were dissected and immediately placed into a 1.5-ml tube containing 700 μl of TRIzol reagent (QIAGEN) on ice. To avoid cross-contamination, we used new disposable dissecting instruments for each sample. Samples were crushed with the Fast Prep 24 instrument (MP Biomedicals). Total RNA was extracted following the manufacturer's instructions, allowing isolation of all RNAs including the small RNAs (miRNeasy Micro Kit; QIAGEN). The quality and the quantity of the RNA samples were determined with a spectrophotometer NanoDrop One (Ozyme) and with an RNA chip (Agilent RNA 6000 Nano kit) into the BioAnalyzer 2100 (Agilent).

Preparation and sequencing of small RNA libraries

The libraries were synthesized from 1 μg of RNA for each sample, following the manufacturer's instructions (NEXTflex Illumina Small RNA Seq Prep; Bioo Scientific Corporation). Briefly, 3′ adenylated adaptaters, followed by 5′ adaptaters, were sequentially added to RNA strands. Then, RNA samples were reverse transcribed using the 3′ adaptaters as the template for the RT primer. Finally, cDNAs were amplified by PCR. Then, the quality and the quantity of cDNA were analyzed with a DNA chip (Agilent).

A size selection of the cDNA libraries (150 pb) was performed with a Pippin Prep (Sage Science). All the samples were pooled in an equimolar manner and purified on column (MinElute PCR Purification Kit; QIAGEN). Finally, the quality and quantity of cDNA were checked with a DNA chip (Agilent) before sequencing. The samples were sent to the functional genomic platform of Nice Sophia Antipolis, France (https://www.france-genomique.org). The next-generation single-read sequencing was performed using the sequencer NextSeq 500, Illumina, using the single-end 75-bp high output sequencing mode. We selected this mode to achieve a minimum of 8 millions of reads per sample to enable accurate detection of novel miRNAs with the miRPro algorithm. The preliminary results are the following:

Total reads: 454,577,868.

Filtered too short reads (size < 15 b): 1,150,734 (0.3% of total).

Not filtered reads:

  • 1) reads ready to be mapped with Bowtie2: 453,427,134 (99.7% of total)

  • 2) unmapped reads, with low quality mm > 2 or contamination reads: 17,632,800 (3.9% of unfiltered reads)

  • 3) mapped reads

    • -ready to be annotated: 435,794,334 (96.1% of unfiltered reads)

    • -not annotated reads, aligned in unannotated regions: 72,073,770 (16.5% of mapped reads)

    • -annotated reads: 363,720,564 (83.5% of mapped reads)

    • -miRNA reads: 363,720,564 (100% of annotated)

The row data were uploaded in the European Nucleotide Archive under the accession number PRJEB24026.

The sequence quality was assessed with FastQC software (www.bioinformatics.babraham.ac.uk/projects/fastqc/).

Identification of conserved and novel miRNAs by bioinformatic analysis

To obtain clean reads, the adaptator sequences (from the preparation of libraries) and the four nucleotides inserted between the adaptor sequence and the sequence of the miRNAs were removed using the Cutadapt tool. The Rattus norgevicus genome reference was downloaded from UCS Genome Browser assembly ID:m6. All the bioinformatic analysis were performed using the program miRPro (Shi et al, 2015), which uses the main algorithm of the software miRDeep2 (Friedländer et al, 2012). The main advantage of miRPro compared with miRDeep2 is that it proposes consistent and unified names for novel precursors and their mature miRNAs in all libraries of the samples. The reads issued from the sequencing are aligned against the genome of reference using Novoalign (Li & Homer, 2010).

Sequences that matched known rat small RNAs such as rRNA, scRNA, snoRNA, snRNA, and tRNA or degradation fragments of mRNAs were excuded in further analysis, and sequences that perfectly matched the rat genome along their entire lenght and recognized as miRNAs were subjected to subsequent analysis.

The small cleaned RNAs aligned against the miRNA database that cannot be annotated to any previous category were subjected to novel miRNA prediction software miRPro (Shi et al, 2015). Briefly, after the matching of the reads with the genome of reference, the genome is scanned from 5′ to 3′ to excise potential precursor sequences. Then, the software creates an index of potential precursors. Afterwards, the reads are aligned with the potential precursors. Finally, the predictive structures are performed by the software RNAfold.

Selection of miRNAs

To avoid false-positive miRNAs, we kept miRNAs found at least in two samples of each structure. Then, we kept miRNAs that obtained more than 1 count per million of sequences, in at least two samples per structure.

Orthologs among new miRNAs discovered

Sequences of new miRNAs were aligned against mouse and human sequences (obtained from miRBase release 21 [http://www.mirbase.org/]) (Kozomara & Griffiths-Jones, 2011, 2014) using Blast software (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Ortholog sequences were determined from a sequence homology more than 95% of identity between the two sequences. Moreover, functional orthologs were determined with a perfect homology into the seed region.

Characterization of novel miRNA nucleotidic sequences

We used the software suite MEME (http://meme-suite.org) to statistically analyze the nucleotidic composition of novel miRNA sequences.

Identification of differentially expressed miRNAs and specific miRNAs

We compared the miRNA expression profiles between the nervous system structures, using DESeq2 (Love et al, 2014) under software R. We considered miRNAs as differentially expressed between structures when P < 0.05. We considered an miRNA as specific of a structure when it was mainly overexpressed or underexpressed in one structure compared with all the others, following these criteria: log2 fold change > 4 or < −4 and P < 0.05.

miRNA target prediction

To understand the function of selected known miRNAs, based on their abundance, and/or their ability to be differentially expressed, we used the algorithm of the database Targetscan 7.1 (http://www.targetscan.org/vert_71/) (Shin et al, 2010; Agarwal et al, 2015). For novel miRNAs, we used miRDB (http://www.mirdb.org/). We selected mRNA targets that have particular known functions in the nervous system. The signaling pathway of mRNA targets was analyzed using GO Enrichment analysis (http://www.geneontology.org/).

Preparation and sequencing of mRNA libraries

The libraries were synthesized from 1 μg of RNA for each sample, following the manufacturer's instructions (NEXTflex Illumina Rapid Directional mRNA-Seq, Bioo Scientific Corporation). Briefly, mRNAs were isolated with poly(A) magnetic beads. Then, poly(A) RNA is fragmented, followed by the first and second strand syntheses. 5′ and 3′ adenylated adaptaters are added, and cDNA is amplified by PCR. Finally, the quality and the quantity of libraries were analyzed with a DNA chip (Agilent) and pooled before sequencing. Sequencing was performed by the Genewiz company on a Illumina HiSeq in the 2 × 150 bp configuration. From the 15 samples, we obtained a total of 408,628,260 reads with a mean quality score of >38.44 and more than 92.48% of bases with a Phred quality score of >30. Reads were aligned with STAR (Dobin et al, 2013), and abundance data (gene counts) were generated with the –quantMode option. Raw mRNA-Seq data and gene count numbers were submitted to the Gene Expression Omnibus database and recorded with the accession number GSE119349.

Luciferase experiment

The wild-type 3′UTR LCE2D luciferase reporter was obtained by annealing 50-bp synthesized oligonucleotides containing the putative mmu-miR-676 binding site. Mutated 3′UTR construct was obtained using the same 50-bp synthesized oligonucleotides, although the putative mmu-miR-676 binding site composed of 8 nt was replaced with antisense nucleotides. Annealed oligonucleotides were then ligated into pmiRGLO vector (Promega) downstream of the Firefly luciferase reporter. The pmiRGLO vector also contains a Renilla luciferase cassette that is used as a transfection normalizer. For the assay, LCE2D luciferase reporter in its wild type or mutated form was co-transfected in COS cells with a pcDNA3.1 plasmid expressing the mmu-miR-676, the novel-miR-21, or the Cel-miR-67 as control. Hence, Cel-miR-67 is a miRNA from C. elegans known to have no target in mammals. The activity of both Firefly and Renilla luciferases was assessed 24 h after transfection using the Dual-Luciferase Reporter Assay System kit (Promega). Translation activity was reported as the ratio between Firefly and Renilla luciferase.

Statistical analysis

To evaluate the differential expression of miRNAs in the five structures, we used the DESeq2 algorithm in R software. The Shapiro–Wilk test was realized with R on the sequence count results to evaluate the normal distribution of the sample. Then, in regard to the distribution of the sample, we performed either a Kruskal–Wallis test, followed by a Dunn multiple comparison test, or a one-way analysis of variance, followed by a Bonferroni post hoc test. To evaluate the differential expression of mRNAs in the five structures, we used the DESeq2 algorithm in R software with the Wald test to calculate the P-value. To assess the difference in the cumulative frequency distribution of expression change of all miRNA target genes versus nontarget genes, we used the Kolmogorov–Smirnov test. The significance level was set at P < 0.05.

Acknowledgements

A Soula was funded by the French Ministry for Research and Fondation pour la Recherche Médicale. Part of the computations, reads alignment, novel miRNA prediction, and gene abundance estimation were performed at the Bordeaux Bioinformatics Center (CBiB). Computer resources were provided by the CBiB, Université de Bordeaux.

Author Contributions

  • A Soula: conceptualization, data curation, formal analysis, investigation, methodology, and writing—original draft, review, and editing.

  • M Valere: data curation, investigation, and methodology.

  • MJ López-González: investigation and methodology.

  • V Ury-Thiery: investigation and methodology.

  • A Groppi: resources, software, methodology, and writing—review and editing.

  • M Landry: conceptualization, investigation, and methodology.

  • M Nikolski: conceptualization, resources, software, writing—review and editing.

  • A Favereaux: conceptualization, formal analysis, funding acquisition, investigation, methodology, project administration, and writing—original draft, review, and editing.

Conflict of Interest Statement

The authors declare that they have no conflict of interest.

  • Received January 9, 2018.
  • Revision received October 10, 2018.
  • Accepted October 11, 2018.
  • © 2018 Soula et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Abellán A,
    2. Desfilis E,
    3. Medina L
    (2014) Combinatorial expression of Lef1, Lhx2, Lhx5, Lhx9, Lmo3, Lmo4, and Prox1 helps to identify comparable subdivisions in the developing hippocampal formation of mouse and chicken. Front Neuroanat 8: 59. doi:10.3389/fnana.2014.00059
    OpenUrlCrossRefPubMed
  2. ↵
    1. Agarwal V,
    2. Bell GW,
    3. Nam J-W,
    4. Bartel DP
    (2015) Predicting effective microRNA target sites in mammalian mRNAs. eLife 4. doi:10.7554/elife.05005
    OpenUrlCrossRef
  3. ↵
    1. Aksoy-Aksel A,
    2. Zampa F,
    3. Schratt G
    (2014) MicroRNAs and synaptic plasticity: A mutual relationship. Philos Trans R Soc Lond B Biol Sci 369: 20130515. doi:10.1098/rstb.2013.0515
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Altschul SF,
    2. Gish W,
    3. Miller W,
    4. Myers EW,
    5. Lipman DJ
    (1990) Basic local alignment search tool. J Mol Biol 215: 403–410. doi:10.1006/jmbi.1990.9999
    OpenUrlCrossRefPubMed
  5. ↵
    1. Altuvia Y,
    2. Landgraf P,
    3. Lithwick G,
    4. Elefant N,
    5. Pfeffer S,
    6. Aravin A,
    7. Brownstein MJ,
    8. Tuschl T,
    9. Margalit H
    (2005) Clustering and conservation patterns of human microRNAs. Nucleic Acids Res 33: 2697–2706. doi:10.1093/nar/gki567
    OpenUrlCrossRefPubMed
  6. ↵
    1. Amar L,
    2. Benoit C,
    3. Beaumont G,
    4. Vacher CM,
    5. Crepin D,
    6. Taouis M,
    7. Baroin-Tourancheau A
    (2012) MicroRNA expression profiling of hypothalamic arcuate and paraventricular nuclei from single rats using Illumina sequencing technology. J Neurosci Methods 209: 134–143. doi:10.1016/j.jneumeth.2012.05.033
    OpenUrlCrossRefPubMed
  7. ↵
    1. Ambros V,
    2. Lee RC,
    3. Lavanway A,
    4. Williams PT,
    5. Jewell D
    (2003) MicroRNAs and other tiny endogenous RNAs in C. elegans. Curr Biol CB 13: 807–818. doi:10.1016/s0960-9822(03)00287-2
    OpenUrlCrossRef
  8. ↵
    1. Arlotta P,
    2. Molyneaux BJ,
    3. Chen J,
    4. Inoue J,
    5. Kominami R,
    6. Macklis JD
    (2005) Neuronal subtype-specific genes that control corticospinal motor neuron development in vivo. Neuron 45: 207–221. doi:10.1016/j.neuron.2004.12.036
    OpenUrlCrossRefPubMed
  9. ↵
    1. Ashburner M,
    2. Ball CA,
    3. Blake JA,
    4. Botstein D,
    5. Butler H,
    6. Cherry JM,
    7. Davis AP,
    8. Dolinski K,
    9. Dwight SS,
    10. Eppig JT, et al.
    (2000) Gene ontology: Tool for the unification of biology. The gene Ontology Consortium. Nat Genet 25: 25–29. doi:10.1038/75556
    OpenUrlCrossRefPubMed
  10. ↵
    1. Asli NS,
    2. Kessel M
    (2010) Spatiotemporally restricted regulation of generic motor neuron programs by miR-196-mediated repression of Hoxb8. Dev Biol 344: 857–868. doi:10.1016/j.ydbio.2010.06.003
    OpenUrlCrossRefPubMed
  11. ↵
    1. Bak M,
    2. Silahtaroglu A,
    3. Møller M,
    4. Christensen M,
    5. Rath MF,
    6. Skryabin B,
    7. Tommerup N,
    8. Kauppinen S
    (2008) MicroRNA expression in the adult mouse central nervous system. RNA N Y N 14: 432–444. doi:10.1261/rna.783108
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Bali P,
    2. Kenny PJ
    (2013) MicroRNAs and drug addiction. Front Genet 4: 43. doi:10.3389/fgene.2013.00043
    OpenUrlCrossRef
  13. ↵
    1. Bartel DP
    (2004) MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 116: 281–297. doi:10.1016/s0092-8674(04)00045-5
    OpenUrlCrossRefPubMed
  14. ↵
    1. Bartel DP
    (2009) MicroRNAs: Target recognition and regulatory functions. Cell 136: 215–233. doi:10.1016/j.cell.2009.01.002
    OpenUrlCrossRefPubMed
  15. ↵
    1. Batista MF,
    2. Lewis KE
    (2008) Pax2/8 act redundantly to specify glycinergic and GABAergic fates of multiple spinal interneurons. Dev Biol 323: 88–97. doi:10.1016/j.ydbio.2008.08.009
    OpenUrlCrossRefPubMed
  16. ↵
    1. Batista-Brito R,
    2. Vinck M,
    3. Ferguson KA,
    4. Chang JT,
    5. Laubender D,
    6. Lur G,
    7. Mossner JM,
    8. Hernandez VG,
    9. Ramakrishnan C,
    10. Deisseroth K, et al.
    (2017) Developmental dysfunction of VIP interneurons impairs cortical circuits. Neuron 95: 884–895. e9. doi:10.1016/j.neuron.2017.07.034
    OpenUrlCrossRefPubMed
  17. ↵
    1. Beclin C,
    2. Follert P,
    3. Stappers E,
    4. Barral S,
    5. Nathalie C,
    6. de Chevigny A,
    7. Magnone V,
    8. Lebrigand K,
    9. Bissels U,
    10. Huylebroeck D, et al.
    (2016) miR-200 family controls late steps of postnatal forebrain neurogenesis via Zeb2 inhibition. Sci Rep 6: 35729. doi:10.1038/srep39368
    OpenUrlCrossRef
  18. ↵
    1. Belvindrah R,
    2. Nissant A,
    3. Lledo P-M
    (2011) Abnormal neuronal migration changes the fate of developing neurons in the postnatal olfactory bulb. J Neurosci Off J Soc Neurosci 31: 7551–7562. doi:10.1523/jneurosci.6716-10.2011
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Bernard O
    (2007) Lim kinases, regulators of actin dynamics. Int J Biochem Cel Biol 39: 1071–1076. doi:10.1016/j.biocel.2006.11.011
    OpenUrlCrossRefPubMed
  20. ↵
    1. Bernstein E,
    2. Caudy AA,
    3. Hammond SM,
    4. Hannon GJ
    (2001) Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature 409: 363–366. doi:10.1038/35053110
    OpenUrlCrossRefPubMed
  21. ↵
    1. Böhm UL,
    2. Prendergast A,
    3. Djenoune L,
    4. Nunes Figueiredo S,
    5. Gomez J,
    6. Stokes C,
    7. Kaiser S,
    8. Suster M,
    9. Kawakami K,
    10. Charpentier M, et al.
    (2016) CSF-contacting neurons regulate locomotion by relaying mechanical stimuli to spinal circuits. Nat Commun 7: 10866. doi:10.1038/ncomms10866
    OpenUrlCrossRefPubMed
  22. ↵
    1. Bredy TW,
    2. Lin Q,
    3. Wei W,
    4. Baker-Andresen D,
    5. Mattick JS
    (2011) MicroRNA regulation of neural plasticity and memory. Neurobiol Learn Mem 96: 89–94. doi:10.1016/j.nlm.2011.04.004
    OpenUrlCrossRefPubMed
  23. ↵
    1. Brumovsky P,
    2. Watanabe M,
    3. Hökfelt T
    (2007) Expression of the vesicular glutamate transporters-1 and -2 in adult mouse dorsal root ganglia and spinal cord and their regulation by nerve injury. Neuroscience 147: 469–490. doi:10.1016/j.neuroscience.2007.02.068
    OpenUrlCrossRefPubMed
  24. ↵
    1. Cai X,
    2. Hagedorn CH,
    3. Cullen BR
    (2004) Human microRNAs are processed from capped, polyadenylated transcripts that can also function as mRNAs. RNA N Y N 10: 1957–1966. doi:10.1261/rna.7135204
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Carrington JC,
    2. Ambros V
    (2003) Role of microRNAs in plant and animal development. Science 301: 336–338. doi:10.1126/science.1085242
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Chen C,
    2. Ridzon DA,
    3. Broomer AJ,
    4. Zhou Z,
    5. Lee DH,
    6. Nguyen JT,
    7. Barbisin M,
    8. Xu NL,
    9. Mahuvakar VR,
    10. Andersen MR, et al.
    (2005) Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res 33: e179. doi:10.1093/nar/gni178
    OpenUrlCrossRefPubMed
  27. ↵
    1. Chen C,
    2. Zhang Y,
    3. Zhang L,
    4. Weakley SM,
    5. Yao Q
    (2011) MicroRNA-196: Critical roles and clinical applications in development and cancer. J Cel Mol Med 15: 14–23. doi:10.1111/j.1582-4934.2010.01219.x
    OpenUrlCrossRefPubMed
  28. ↵
    1. Cheng L,
    2. Duan B,
    3. Huang T,
    4. Zhang Y,
    5. Chen Y,
    6. Britz O,
    7. Garcia-Campmany L,
    8. Ren X,
    9. Vong L,
    10. Lowell BB, et al.
    (2017) Identification of spinal circuits involved in touch-evoked dynamic mechanical pain. Nat Neurosci 20: 804–814. doi:10.1038/nn.4549
    OpenUrlCrossRefPubMed
  29. ↵
    1. Choi PS,
    2. Zakhary L,
    3. Choi W-Y,
    4. Caron S,
    5. Alvarez-Saavedra E,
    6. Miska EA,
    7. McManus M,
    8. Harfe B,
    9. Giraldez AJ,
    10. Horvitz HR, et al.
    (2008) Members of the miRNA-200 family regulate olfactory neurogenesis. Neuron 57: 41–55. doi:10.1016/j.neuron.2007.11.018
    OpenUrlCrossRefPubMed
  30. ↵
    1. Conaco C,
    2. Otto S,
    3. Han J-J,
    4. Mandel G
    (2006) Reciprocal actions of REST and a microRNA promote neuronal identity. Proc Natl Acad Sci U S A 103: 2422–2427. doi:10.1073/pnas.0511041103
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Desvignes T,
    2. Batzel P,
    3. Berezikov E,
    4. Eilbeck K,
    5. Eppig JT,
    6. McAndrews MS,
    7. Singer A,
    8. Postlethwait JH
    (2015) miRNA nomenclature: A view incorporating genetic origins, biosynthetic pathways, and sequence variants. Trends Genet TIG 31: 613–626. doi:10.1016/j.tig.2015.09.002
    OpenUrlCrossRefPubMed
  32. ↵
    1. Díaz NF,
    2. Cruz-Reséndiz MS,
    3. Flores-Herrera H,
    4. García-López G,
    5. Molina-Hernández A
    (2014) MicroRNAs in central nervous system development. Rev Neurosci 25: 675–686. doi:10.1515/revneuro-2014-0014
    OpenUrlCrossRefPubMed
  33. ↵
    1. Dibbens LM,
    2. Tarpey PS,
    3. Hynes K,
    4. Bayly MA,
    5. Scheffer IE,
    6. Smith R,
    7. Bomar J,
    8. Sutton E,
    9. Vandeleur L,
    10. Shoubridge C, et al.
    (2008) X-linked protocadherin 19 mutations cause female-limited epilepsy and cognitive impairment. Nat Genet 40: 776–781. doi:10.1038/ng.149
    OpenUrlCrossRefPubMed
  34. ↵
    1. Ding Q,
    2. Balasubramanian R,
    3. Zheng D,
    4. Liang G,
    5. Gan L
    (2017) Barhl2 determines the early patterning of the diencephalon by regulating Shh. Mol Neurobiol 54: 4414–4420. doi:10.1007/s12035-016-0001-5
    OpenUrlCrossRef
  35. ↵
    1. Ding Y-Q,
    2. Yin J,
    3. Kania A,
    4. Zhao Z-Q,
    5. Johnson RL,
    6. Chen Z-F
    (2004) Lmx1b controls the differentiation and migration of the superficial dorsal horn neurons of the spinal cord. Dev Camb Engl 131: 3693–3703. doi:10.1242/dev.01250
    OpenUrlCrossRef
  36. ↵
    1. Djenoune L,
    2. Khabou H,
    3. Joubert F,
    4. Quan FB,
    5. Nunes Figueiredo S,
    6. Bodineau L,
    7. Del Bene F,
    8. Burcklé C,
    9. Tostivint H,
    10. Wyart C
    (2014) Investigation of spinal cerebrospinal fluid-contacting neurons expressing PKD2L1: Evidence for a conserved system from fish to primates. Front Neuroanat 8: 26. doi:10.3389/fnana.2014.00026
    OpenUrlCrossRefPubMed
  37. ↵
    1. Dobin A,
    2. Davis CA,
    3. Schlesinger F,
    4. Drenkow J,
    5. Zaleski C,
    6. Jha S,
    7. Batut P,
    8. Chaisson M,
    9. Gingeras TR
    (2013) STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 29: 15–21. doi:10.1093/bioinformatics/bts635
    OpenUrlCrossRefPubMed
  38. ↵
    1. Elramah S,
    2. Landry M,
    3. Favereaux A
    (2014) MicroRNAs regulate neuronal plasticity and are involved in pain mechanisms. Front Cel Neurosci 8: 31. doi:10.3389/fncel.2014.00031
    OpenUrlCrossRefPubMed
  39. ↵
    1. Elramah S,
    2. López-González MJ,
    3. Bastide M,
    4. Dixmérias F,
    5. Roca-Lapirot O,
    6. Wielanek-Bachelet A-C,
    7. Vital A,
    8. Leste-Lasserre T,
    9. Brochard A,
    10. Landry M, et al.
    (2017) Spinal miRNA-124 regulates synaptopodin and nociception in an animal model of bone cancer pain. Sci Rep 7: 10949. doi:10.1038/s41598-017-10224-1
    OpenUrlCrossRef
  40. ↵
    1. Engels BM,
    2. Hutvagner G
    (2006) Principles and effects of microRNA-mediated post-transcriptional gene regulation. Oncogene 25: 6163–6169. doi:10.1038/sj.onc.1209909
    OpenUrlCrossRefPubMed
  41. ↵
    1. Favereaux A,
    2. Thoumine O,
    3. Bouali-Benazzouz R,
    4. Roques V,
    5. Papon M-A,
    6. Salam SA,
    7. Drutel G,
    8. Léger C,
    9. Calas A,
    10. Nagy F, et al.
    (2011) Bidirectional integrative regulation of Cav1.2 calcium channel by microRNA mir-103: Role in pain. EMBO J 30: 3830–3841. doi:10.1038/emboj.2011.249
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Fenster CP,
    2. Fenster SD,
    3. Leahy HP,
    4. Kurschner C,
    5. Blundon JA
    (2007) Modulation of Kv4.2 K+ currents by neuronal interleukin-16, a PDZ domain-containing protein expressed in the hippocampus and cerebellum. Brain Res 1162: 19–31. doi:10.1016/j.brainres.2007.05.051
    OpenUrlCrossRefPubMed
  43. ↵
    1. Frank F,
    2. Sonenberg N,
    3. Nagar B
    (2010) Structural basis for 5ʹ-nucleotide base-specific recognition of guide RNA by human AGO2. Nature 465: 818–822. doi:10.1038/nature09039
    OpenUrlCrossRefPubMed
  44. ↵
    1. Friedländer MR,
    2. Lizano E,
    3. Houben AJ,
    4. Bezdan D,
    5. Báñez-Coronel M,
    6. Kudla G,
    7. Mateu-Huertas E,
    8. Kagerbauer B,
    9. González J,
    10. Chen KC, et al.
    (2014) Evidence for the biogenesis of more than 1,000 novel human microRNAs. Genome Biol 15: R57. doi:10.1186/gb-2014-15-4-r57
    OpenUrlCrossRefPubMed
  45. ↵
    1. Friedländer MR,
    2. Mackowiak SD,
    3. Li N,
    4. Chen W,
    5. Rajewsky N
    (2012) miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res 40: 37–52. doi:10.1093/nar/gkr688
    OpenUrlCrossRefPubMed
  46. ↵
    1. Friedman RC,
    2. Farh KK-H,
    3. Burge CB,
    4. Bartel DP
    (2009) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 19: 92–105. doi:10.1101/gr.082701.108
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Gao J,
    2. Wang W-Y,
    3. Mao Y-W,
    4. Gräff J,
    5. Guan J-S,
    6. Pan L,
    7. Mak G,
    8. Kim D,
    9. Su SC,
    10. Tsai L-H
    (2010) A novel pathway regulates memory and plasticity via SIRT1 and miR-134. Nature 466: 1105–1109. doi:10.1038/nature09271
    OpenUrlCrossRefPubMed
  48. ↵
    1. Gärtner A,
    2. Polnau DG,
    3. Staiger V,
    4. Sciarretta C,
    5. Minichiello L,
    6. Thoenen H,
    7. Bonhoeffer T,
    8. Korte M
    (2006) Hippocampal long-term potentiation is supported by presynaptic and postsynaptic tyrosine receptor kinase B-mediated phospholipase Cgamma signaling. J Neurosci Off J Soc Neurosci 26: 3496–3504. doi:10.1523/jneurosci.3792-05.2006
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Gong N,
    2. Fan H,
    3. Ma A-N,
    4. Xiao Q,
    5. Wang Y-X
    (2014a) Geniposide and its iridoid analogs exhibit antinociception by acting at the spinal GLP-1 receptors. Neuropharmacology 84: 31–45. doi:10.1016/j.neuropharm.2014.04.007
    OpenUrlCrossRefPubMed
  50. ↵
    1. Gong N,
    2. Xiao Q,
    3. Zhu B,
    4. Zhang C-Y,
    5. Wang Y-C,
    6. Fan H,
    7. Ma A-N,
    8. Wang Y-X
    (2014b) Activation of spinal glucagon-like peptide-1 receptors specifically suppresses pain hypersensitivity. J Neurosci Off J Soc Neurosci 34: 5322–5334. doi:10.1523/jneurosci.4703-13.2014
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Griffiths-Jones S
    (2006) miRBase: The microRNA sequence database. Methods Mol Biol 342: 129–138. doi:10.1385/1-59745-123-1:129
    OpenUrlCrossRefPubMed
  52. ↵
    1. Gross MK,
    2. Dottori M,
    3. Goulding M
    (2002) Lbx1 specifies somatosensory association interneurons in the dorsal spinal cord. Neuron 34: 535–549. doi:10.1016/s0896-6273(02)00690-6
    OpenUrlCrossRefPubMed
  53. ↵
    1. Hafner M,
    2. Landgraf P,
    3. Ludwig J,
    4. Rice A,
    5. Ojo T,
    6. Lin C,
    7. Holoch D,
    8. Lim C,
    9. Tuschl T
    (2008) Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing. Methods San Diego Calif 44: 3–12. doi:10.1016/j.ymeth.2007.09.009
    OpenUrlCrossRefPubMed
  54. ↵
    1. Herdegen T,
    2. Tölle TR,
    3. Bravo R,
    4. Zieglgänsberger W,
    5. Zimmermann M
    (1991) Sequential expression of JUN B, JUN D and FOS B proteins in rat spinal neurons: Cascade of transcriptional operations during nociception. Neurosci Lett 129: 221–224. doi:10.1016/0304-3940(91)90466-7
    OpenUrlCrossRefPubMed
  55. ↵
    1. Hilinski WC,
    2. Bostrom JR,
    3. England SJ,
    4. Juárez-Morales JL,
    5. de Jager S,
    6. Armant O,
    7. Legradi J,
    8. Strähle U,
    9. Link BA,
    10. Lewis KE
    (2016) Lmx1b is required for the glutamatergic fates of a subset of spinal cord neurons. Neural Develop 11: 16. doi:10.1186/s13064-016-0070-1
    OpenUrlCrossRef
  56. ↵
    1. Hofacker IL
    (2003) Vienna RNA secondary structure server. Nucleic Acids Res 31: 3429–3431. doi:10.1093/nar/gkg599
    OpenUrlCrossRefPubMed
  57. ↵
    1. Hossaini M,
    2. Jongen JLM,
    3. Biesheuvel K,
    4. Kuhl D,
    5. Holstege JC
    (2010) Nociceptive stimulation induces expression of Arc/Arg3.1 in the spinal cord with a preference for neurons containing enkephalin. Mol Pain 6: 43. doi:10.1186/1744-8069-6-43
    OpenUrlCrossRefPubMed
  58. ↵
    1. Ji Z,
    2. Wang G,
    3. Xie Z,
    4. Zhang C,
    5. Wang J
    (2012) Identification and characterization of microRNA in the dairy goat (Capra hircus) mammary gland by Solexa deep-sequencing technology. Mol Biol Rep 39: 9361–9371. doi:10.1007/s11033-012-1779-5
    OpenUrlCrossRefPubMed
  59. ↵
    1. Kang W,
    2. Friedländer MR
    (2015) Computational prediction of miRNA genes from small RNA sequencing data. Front Bioeng Biotechnol 3: 7. doi:10.3389/fbioe.2015.00007
    OpenUrlCrossRef
  60. ↵
    1. Kim SY,
    2. Mo JW,
    3. Han S,
    4. Choi SY,
    5. Han SB,
    6. Moon BH,
    7. Rhyu IJ,
    8. Sun W,
    9. Kim H
    (2010) The expression of non-clustered protocadherins in adult rat hippocampal formation and the connecting brain regions. Neuroscience 170: 189–199. doi:10.1016/j.neuroscience.2010.05.027
    OpenUrlCrossRefPubMed
  61. ↵
    1. Kozomara A,
    2. Griffiths-Jones S
    (2011) miRBase: Integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res 39: D152–D157. doi:10.1093/nar/gkq1027
    OpenUrlCrossRefPubMed
  62. ↵
    1. Kozomara A,
    2. Griffiths-Jones S
    (2014) miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res 42: D68–D73. doi:10.1093/nar/gkt1181
    OpenUrlCrossRefPubMed
  63. ↵
    1. Landry M,
    2. Bouali-Benazzouz R,
    3. El Mestikawy S,
    4. Ravassard P,
    5. Nagy F
    (2004) Expression of vesicular glutamate transporters in rat lumbar spinal cord, with a note on dorsal root ganglia. J Comp Neurol 468: 380–394. doi:10.1002/cne.10988
    OpenUrlCrossRefPubMed
  64. ↵
    1. Larsson M
    (2017) Pax2 is persistently expressed by GABAergic neurons throughout the adult rat dorsal horn. Neurosci Lett 638: 96–101. doi:10.1016/j.neulet.2016.12.015
    OpenUrlCrossRefPubMed
  65. ↵
    1. Lee Y,
    2. Jeon K,
    3. Lee J-T,
    4. Kim S,
    5. Kim VN
    (2002) MicroRNA maturation: Stepwise processing and subcellular localization. EMBO J 21: 4663–4670. doi:10.1093/emboj/cdf476
    OpenUrlAbstract
  66. ↵
    1. Lee Y,
    2. Kim M,
    3. Han J,
    4. Yeom K-H,
    5. Lee S,
    6. Baek SH,
    7. Kim VN
    (2004) MicroRNA genes are transcribed by RNA polymerase II. EMBO J 23: 4051–4060. doi:10.1038/sj.emboj.7600385
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Letellier M,
    2. Elramah S,
    3. Mondin M,
    4. Soula A,
    5. Penn A,
    6. Choquet D,
    7. Landry M,
    8. Thoumine O,
    9. Favereaux A
    (2014) miR-92a regulates expression of synaptic GluA1-containing AMPA receptors during homeostatic scaling. Nat Neurosci 17: 1040–1042. doi:10.1038/nn.3762
    OpenUrlCrossRefPubMed
  68. ↵
    1. Li G,
    2. Ling S
    (2017) MiR-124 promotes newborn olfactory bulb neuron dendritic morphogenesis and spine density. J Mol Neurosci MN 61: 159–168. doi:10.1007/s12031-016-0873-x
    OpenUrlCrossRef
  69. ↵
    1. Li H,
    2. Homer N
    (2010) A survey of sequence alignment algorithms for next-generation sequencing. Brief Bioinform 11: 473–483. doi:10.1093/bib/bbq015
    OpenUrlCrossRefPubMed
  70. ↵
    1. Li J,
    2. Wang C,
    3. Zhang Z,
    4. Wen Y,
    5. An L,
    6. Liang Q,
    7. Xu Z,
    8. Wei S,
    9. Li W,
    10. Guo T, et al.
    (2018) Transcription factors Sp8 and Sp9 coordinately regulate olfactory bulb interneuron development. Cereb Cortex 28: 3278–3294. doi:10.1093/cercor/bhx199
    OpenUrlCrossRef
  71. ↵
    1. Li X-F,
    2. Cao R-B,
    3. Luo J,
    4. Fan J-M,
    5. Wang J-M,
    6. Zhang Y-P,
    7. Gu J-Y,
    8. Feng X-L,
    9. Zhou B,
    10. Chen P-Y
    (2016) MicroRNA transcriptome profiling of mice brains infected with Japanese encephalitis virus by RNA sequencing. Infect Genet Evol J Mol Epidemiol Evol Genet Infect Dis 39: 249–257. doi:10.1016/j.meegid.2016.01.028
    OpenUrlCrossRef
  72. ↵
    1. Li Y,
    2. Zhang Z,
    3. Liu F,
    4. Vongsangnak W,
    5. Jing Q,
    6. Shen B
    (2012) Performance comparison and evaluation of software tools for microRNA deep-sequencing data analysis. Nucleic Acids Res 40: 4298–4305. doi:10.1093/nar/gks043
    OpenUrlCrossRefPubMed
  73. ↵
    1. Lim LP,
    2. Lau NC,
    3. Garrett-Engele P,
    4. Grimson A,
    5. Schelter JM,
    6. Castle J,
    7. Bartel DP,
    8. Linsley PS,
    9. Johnson JM
    (2005) Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature 433: 769–773. doi:10.1038/nature03315
    OpenUrlCrossRefPubMed
  74. ↵
    1. Liu Q,
    2. He H,
    3. Zeng T,
    4. Huang Z,
    5. Fan T,
    6. Wu Q
    (2014) Neural-specific expression of miR-344-3p during mouse embryonic development. J Mol Histol 45: 363–372. doi:10.1007/s10735-013-9555-y
    OpenUrlCrossRef
  75. ↵
    1. López-González MJ,
    2. Landry M,
    3. Favereaux A
    (2017) MicroRNA and chronic pain: From mechanisms to therapeutic potential. Pharmacol Ther 180: 1–15. doi:10.1016/j.pharmthera.2017.06.001
    OpenUrlCrossRef
  76. ↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15: 550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  77. ↵
    1. Maciotta S,
    2. Meregalli M,
    3. Torrente Y
    (2013) The involvement of microRNAs in neurodegenerative diseases. Front Cel Neurosci 7: 265. doi:10.3389/fncel.2013.00265
    OpenUrlCrossRefPubMed
  78. ↵
    1. Makeyev EV,
    2. Zhang J,
    3. Carrasco MA,
    4. Maniatis T
    (2007) The MicroRNA miR-124 promotes neuronal differentiation by triggering brain-specific alternative pre-mRNA splicing. Mol Cel 27: 435–448. doi:10.1016/j.molcel.2007.07.015
    OpenUrlCrossRefPubMed
  79. ↵
    1. Morishita H,
    2. Miwa JM,
    3. Heintz N,
    4. Hensch TK
    (2010) Lynx1, a cholinergic brake, limits plasticity in adult visual cortex. Science 330: 1238–1240. doi:10.1126/science.1195320
    OpenUrlAbstract/FREE Full Text
  80. ↵
    1. Nelson PT,
    2. Wang W-X,
    3. Rajeev BW
    (2008) MicroRNAs (miRNAs) in neurodegenerative diseases. Brain Pathol Zurich Switz 18: 130–138. doi:10.1111/j.1750-3639.2007.00120.x
    OpenUrlCrossRef
  81. ↵
    1. Olsen L,
    2. Klausen M,
    3. Helboe L,
    4. Nielsen FC,
    5. Werge T
    (2009) MicroRNAs show mutually exclusive expression patterns in the brain of adult male rats. PLoS One 4: e7225. doi:10.1371/journal.pone.0007225
    OpenUrlCrossRefPubMed
  82. ↵
    1. Pacholewska A,
    2. Mach N,
    3. Mata X,
    4. Vaiman A,
    5. Schibler L,
    6. Barrey E,
    7. Gerber V
    (2016) Novel equine tissue miRNAs and breed-related miRNA expressed in serum. BMC Genomics 17: 831. doi:10.1186/s12864-016-3168-2
    OpenUrlCrossRef
  83. ↵
    1. Paina S,
    2. Garzotto D,
    3. DeMarchis S,
    4. Marino M,
    5. Moiana A,
    6. Conti L,
    7. Cattaneo E,
    8. Perera M,
    9. Corte G,
    10. Calautti E, et al.
    (2011) Wnt5a is a transcriptional target of Dlx homeogenes and promotes differentiation of interneuron progenitors in vitro and in vivo. J Neurosci Off J Soc Neurosci 31: 2675–2687. doi:10.1523/jneurosci.3110-10.2011
    OpenUrlAbstract/FREE Full Text
  84. ↵
    1. Pandey A,
    2. Singh P,
    3. Jauhari A,
    4. Singh T,
    5. Khan F,
    6. Pant AB,
    7. Parmar D,
    8. Yadav S
    (2015) Critical role of the miR-200 family in regulating differentiation and proliferation of neurons. J Neurochem 133: 640–652. doi:10.1111/jnc.13089
    OpenUrlCrossRefPubMed
  85. ↵
    1. Parish EV,
    2. Mason JO,
    3. Price DJ
    (2016) Expression of Barhl2 and its relationship with Pax6 expression in the forebrain of the mouse embryo. BMC Neurosci 17: 76. doi:10.1186/s12868-016-0311-6
    OpenUrlCrossRef
  86. ↵
    1. Pathania M,
    2. Torres-Reveron J,
    3. Yan L,
    4. Kimura T,
    5. Lin TV,
    6. Gordon V,
    7. Teng Z-Q,
    8. Zhao X,
    9. Fulga TA,
    10. Van Vactor D, et al.
    (2012) miR-132 enhances dendritic morphogenesis, spine density, synaptic integration, and survival of newborn olfactory bulb neurons. PLoS One 7: e38174. doi:10.1371/journal.pone.0038174
    OpenUrlCrossRefPubMed
  87. ↵
    1. Pino D,
    2. Choe Y,
    3. Pleasure SJ
    (2011) Wnt5a controls neurite development in olfactory bulb interneurons. ASN Neuro 3: e00059. doi:10.1042/an20100038
    OpenUrlCrossRefPubMed
  88. ↵
    1. Rajman M,
    2. Metge F,
    3. Fiore R,
    4. Khudayberdiev S,
    5. Aksoy-Aksel A,
    6. Bicker S,
    7. Ruedell Reschke C,
    8. Raoof R,
    9. Brennan GP,
    10. Delanty N, et al.
    (2017) A microRNA-129-5p/Rbfox crosstalk coordinates homeostatic downscaling of excitatory synapses. EMBO J 36: 1770–1787. doi:10.15252/embj.201695748
    OpenUrlCrossRef
  89. ↵
    1. Redemann-Fibi B,
    2. Schuermann M,
    3. Müller R
    (1991) Stage and tissue-specific expression of fosB during mouse development. Differ Res Biol Divers 46: 43–49. doi:10.1111/j.1432-0436.1991.tb00864.x
    OpenUrlCrossRefPubMed
  90. ↵
    1. Sander M,
    2. Paydar S,
    3. Ericson J,
    4. Briscoe J,
    5. Berber E,
    6. German M,
    7. Jessell TM,
    8. Rubenstein JL
    (2000) Ventral neural patterning by Nkx homeobox genes: Nkx6.1 controls somatic motor neuron and ventral interneuron fates. Genes Dev 14: 2134–2139. doi:10.1101/gad.820400
    OpenUrlAbstract/FREE Full Text
  91. ↵
    1. Sarkar S,
    2. Jun S,
    3. Rellick S,
    4. Quintana DD,
    5. Cavendish JZ,
    6. Simpkins JW
    (2016) Expression of microRNA-34a in Alzheimer’s disease brain targets genes linked to synaptic plasticity, energy metabolism, and resting state network activity. Brain Res 1646: 139–151. doi:10.1016/j.brainres.2016.05.026
    OpenUrlCrossRef
  92. ↵
    1. Schang A-L,
    2. Ngô-Muller V,
    3. Bleux C,
    4. Granger A,
    5. Chenut M-C,
    6. Loudes C,
    7. Magre S,
    8. Counis R,
    9. Cohen-Tannoudji J,
    10. Laverrière J-N
    (2011) GnRH receptor gene expression in the developing rat hippocampus: Transcriptional regulation and potential roles in neuronal plasticity. Endocrinology 152: 568–580. doi:10.1210/en.2010-0840
    OpenUrlCrossRefPubMed
  93. ↵
    1. Schratt GM,
    2. Tuebing F,
    3. Nigh EA,
    4. Kane CG,
    5. Sabatini ME,
    6. Kiebler M,
    7. Greenberg ME
    (2006) A brain-specific microRNA regulates dendritic spine development. Nature 439: 283–289. doi:10.1038/nature04367
    OpenUrlCrossRefPubMed
  94. ↵
    1. Seigneur E,
    2. Südhof TC
    (2017) Cerebellins are differentially expressed in selective subsets of neurons throughout the brain. J Comp Neurol 525: 3286–3311. doi:10.1002/cne.24278
    OpenUrlCrossRefPubMed
  95. ↵
    1. Senfter D,
    2. Madlener S,
    3. Krupitza G,
    4. Mader RM
    (2016) The microRNA-200 family: Still much to discover. Biomol Concepts 7: 311–319. doi:10.1515/bmc-2016-0020
    OpenUrlCrossRef
  96. ↵
    1. Shi J,
    2. Dong M,
    3. Li L,
    4. Liu L,
    5. Luz-Madrigal A,
    6. Tsonis PA,
    7. Del Rio-Tsonis K,
    8. Liang C
    (2015) mirPRo-a novel standalone program for differential expression and variation analysis of miRNAs. Sci Rep 5: 14617. doi:10.1038/srep14617
    OpenUrlCrossRef
  97. ↵
    1. Shi R,
    2. Chiang VL
    (2005) Facile means for quantifying microRNA expression by real-time PCR. BioTechniques 39: 519–525. doi:10.2144/000112010
    OpenUrlCrossRefPubMed
  98. ↵
    1. Shin C,
    2. Nam J-W,
    3. Farh KK-H,
    4. Chiang HR,
    5. Shkumatava A,
    6. Bartel DP
    (2010) Expanding the microRNA targeting code: Functional sites with centered pairing. Mol Cel 38: 789–802. doi:10.1016/j.molcel.2010.06.005
    OpenUrlCrossRefPubMed
  99. ↵
    1. Shinohara Y,
    2. Yahagi K,
    3. Kawano M,
    4. Nishiyori H,
    5. Kawazu C,
    6. Suzuki N,
    7. Manabe R,
    8. Hirase H
    (2011) miRNA profiling of bilateral rat hippocampal CA3 by deep sequencing. Biochem Biophys Res Commun 409: 293–298. doi:10.1016/j.bbrc.2011.05.004
    OpenUrlCrossRefPubMed
  100. ↵
    1. Sutton MA,
    2. Ito HT,
    3. Cressy P,
    4. Kempf C,
    5. Woo JC,
    6. Schuman EM
    (2006) Miniature neurotransmission stabilizes synaptic function via tonic suppression of local dendritic protein synthesis. Cell 125: 785–799. doi:10.1016/j.cell.2006.03.040
    OpenUrlCrossRefPubMed
  101. ↵
    1. Tan L,
    2. Yu J-T,
    3. Tan L
    (2015) Causes and consequences of MicroRNA dysregulation in neurodegenerative diseases. Mol Neurobiol 51: 1249–1262. doi:10.1007/s12035-014-8803-9
    OpenUrlCrossRefPubMed
  102. ↵
    1. Tehler D,
    2. Høyland-Kroghsbo NM,
    3. Lund AH
    (2011) The miR-10 microRNA precursor family. RNA Biol 8: 728–734. doi:10.4161/rna.8.5.16324
    OpenUrlCrossRefPubMed
  103. ↵
    1. The Gene Ontology Consortium
    (2017) Expansion of the gene ontology knowledgebase and resources. Nucleic Acids Res 45: D331–D338. doi:10.1093/nar/gkw1108
    OpenUrlCrossRefPubMed
  104. ↵
    1. Tripathi S,
    2. Pohl MO,
    3. Zhou Y,
    4. Rodriguez-Frandsen A,
    5. Wang G,
    6. Stein DA,
    7. Moulton HM,
    8. DeJesus P,
    9. Che J,
    10. Mulder LCF, et al.
    (2015) Meta- and orthogonal integration of influenza ‘OMICs’ data defines a role for UBR4 in virus budding. Cell Host Microbe 18: 723–735. doi:10.1016/j.chom.2015.11.002
    OpenUrlCrossRefPubMed
  105. ↵
    1. Trontti K,
    2. Väänänen J,
    3. Sipilä T,
    4. Greco D,
    5. Hovatta I
    (2018) Strong conservation of inbred mouse strain microRNA loci but broad variation in brain microRNAs due to RNA editing and isomiR expression. RNA N Y N 24: 643–655. doi:10.1261/rna.064881.117
    OpenUrlAbstract/FREE Full Text
  106. ↵
    1. Vallstedt A,
    2. Muhr J,
    3. Pattyn A,
    4. Pierani A,
    5. Mendelsohn M,
    6. Sander M,
    7. Jessell TM,
    8. Ericson J
    (2001) Different levels of repressor activity assign redundant and specific roles to Nkx6 genes in motor neuron and interneuron specification. Neuron 31: 743–755. doi:10.1016/s0896-6273(01)00412-3
    OpenUrlCrossRefPubMed
  107. ↵
    1. Vincent K,
    2. Cornea VM,
    3. Jong Y-JI,
    4. Laferrière A,
    5. Kumar N,
    6. Mickeviciute A,
    7. Fung JST,
    8. Bandegi P,
    9. Ribeiro-da-Silva A,
    10. O’Malley KL, et al.
    (2016) Intracellular mGluR5 plays a critical role in neuropathic pain. Nat Commun 7: 10604. doi:10.1038/ncomms10604
    OpenUrlCrossRefPubMed
  108. ↵
    1. Volkow ND,
    2. Wang G-J,
    3. Telang F,
    4. Fowler JS,
    5. Logan J,
    6. Childress A-R,
    7. Jayne M,
    8. Ma Y,
    9. Wong C
    (2006) Cocaine cues and dopamine in dorsal striatum: Mechanism of craving in cocaine addiction. J Neurosci Off J Soc Neurosci 26: 6583–6588. doi:10.1523/jneurosci.1544-06.2006
    OpenUrlAbstract/FREE Full Text
  109. ↵
    1. Wang X
    (2014) Composition of seed sequence is a major determinant of microRNA targeting patterns. Bioinforma Oxf Engl 30: 1377–1383. doi:10.1093/bioinformatics/btu045
    OpenUrlCrossRefPubMed
  110. ↵
    1. Wei LQ,
    2. Yan LF,
    3. Wang T
    (2011) Deep sequencing on genome-wide scale reveals the unique composition and expression patterns of microRNAs in developing pollen of Oryza sativa. Genome Biol 12: R53. doi:10.1186/gb-2011-12-6-r53
    OpenUrlCrossRefPubMed
  111. ↵
    1. Willenbrock H,
    2. Salomon J,
    3. Søkilde R,
    4. Barken KB,
    5. Hansen TN,
    6. Nielsen FC,
    7. Møller S,
    8. Litman T
    (2009) Quantitative miRNA expression analysis: Comparing microarrays with next-generation sequencing. RNA N Y N 15: 2028–2034. doi:10.1261/rna.1699809
    OpenUrlAbstract/FREE Full Text
  112. ↵
    1. Williamson V,
    2. Kim A,
    3. Xie B,
    4. McMichael GO,
    5. Gao Y,
    6. Vladimirov V
    (2013) Detecting miRNAs in deep-sequencing data: A software performance comparison and evaluation. Brief Bioinform 14: 36–45. doi:10.1093/bib/bbs010
    OpenUrlCrossRefPubMed
  113. ↵
    1. Xiong Q,
    2. He W,
    3. Wang H,
    4. Zhou J,
    5. Zhang Y,
    6. He J,
    7. Yang C,
    8. Zhang B
    (2017) Effect of the spinal apelin-APJ system on the pathogenesis of chronic constriction injury-induced neuropathic pain in rats. Mol Med Rep 16: 1223–1231. doi:10.3892/mmr.2017.6734
    OpenUrlCrossRef
  114. ↵
    1. Yao M-J,
    2. Chen G,
    3. Zhao P-P,
    4. Lu M-H,
    5. Jian J,
    6. Liu M-F,
    7. Yuan X-B
    (2012) Transcriptome analysis of microRNAs in developing cerebral cortex of rat. BMC Genomics 13: 232. doi:10.1186/1471-2164-13-232
    OpenUrlCrossRefPubMed
  115. ↵
    1. Yau TO,
    2. Kwan CT,
    3. Jakt LM,
    4. Stallwood N,
    5. Cordes S,
    6. Sham MH
    (2002) Auto/cross-regulation of Hoxb3 expression in posterior hindbrain and spinal cord. Dev Biol 252: 287–300. doi:10.1006/dbio.2002.0849
    OpenUrlCrossRefPubMed
  116. ↵
    1. Yin JQ,
    2. Zhao RC,
    3. Morris KV
    (2008) Profiling microRNA expression with microarrays. Trends Biotechnol 26: 70–76. doi:10.1016/j.tibtech.2007.11.007
    OpenUrlCrossRefPubMed
  117. ↵
    1. Yin L,
    2. Sun Y,
    3. Wu J,
    4. Yan S,
    5. Deng Z,
    6. Wang J,
    7. Liao S,
    8. Yin D,
    9. Li G
    (2015) Discovering novel microRNAs and age-related nonlinear changes in rat brains using deep sequencing. Neurobiol Aging 36: 1037–1044. doi:10.1016/j.neurobiolaging.2014.11.001
    OpenUrlCrossRef
  118. ↵
    1. Zaghetto AA,
    2. Paina S,
    3. Mantero S,
    4. Platonova N,
    5. Peretto P,
    6. Bovetti S,
    7. Puche A,
    8. Piccolo S,
    9. Merlo GR
    (2007) Activation of the Wnt-beta catenin pathway in a cell population on the surface of the forebrain is essential for the establishment of olfactory axon connections. J Neurosci Off J Soc Neurosci 27: 9757–9768. doi:10.1523/jneurosci.0763-07.2007
    OpenUrlAbstract/FREE Full Text
  119. ↵
    1. Zhou M,
    2. Wang Q,
    3. Sun J,
    4. Li X,
    5. Xu L,
    6. Yang H,
    7. Shi H,
    8. Ning S,
    9. Chen L,
    10. Li Y, et al.
    (2009) In silico detection and characteristics of novel microRNA genes in the Equus caballus genome using an integrated ab initio and comparative genomic approach. Genomics 94: 125–131. doi:10.1016/j.ygeno.2009.04.006
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Small RNA-Seq reveals novel miRNAs shaping the transcriptomic identity of rat brain structures
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
miRNAs shape brain transcriptome
Anaïs Soula, Mélissa Valere, María-José López-González, Vicky Ury-Thiery, Alexis Groppi, Marc Landry, Macha Nikolski, Alexandre Favereaux
Life Science Alliance Oct 2018, 1 (5) e201800018; DOI: 10.26508/lsa.201800018

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
miRNAs shape brain transcriptome
Anaïs Soula, Mélissa Valere, María-José López-González, Vicky Ury-Thiery, Alexis Groppi, Marc Landry, Macha Nikolski, Alexandre Favereaux
Life Science Alliance Oct 2018, 1 (5) e201800018; DOI: 10.26508/lsa.201800018
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 1, No. 5
October 2018
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Genomics & Functional Genomics
  • Neuroscience

Related Articles

  • No related articles found.

Cited By...

  • mir124-dependent tagging of glutamatergic synapses by synaptopodin controls non-uniform and input-specific homeostatic synaptic plasticity
  • Google Scholar

More in this TOC Section

  • Proteome profiling in heart valve disease
  • IARA: a curated genetic atlas of human spliceosome
  • Concordance of MERFISH with single-cell RNA-Seq
Show more Resource

Similar Articles

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2023 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.