Genetic variation in hippocampal microRNA expression differences in C57BL/6 J X DBA/2 J (BXD) recombinant inbred mouse strains
BMC Genomics volume 13, Article number: 476 (2012)
miRNAs are short single-stranded non-coding RNAs involved in post-transcriptional gene regulation that play a major role in normal biological functions and diseases. Little is currently known about how expression of miRNAs is regulated. We surveyed variation in miRNA abundance in the hippocampus of mouse inbred strains, allowing us to take a genetic approach to the study of miRNA regulation, which is novel for miRNAs. The BXD recombinant inbred panel is a very well characterized genetic reference panel which allows quantitative trait locus (QTL) analysis of miRNA abundance and detection of correlates in a large store of brain and behavioural phenotypes.
We found five suggestive trans QTLs for the regulation of miRNAs investigated. Further analysis of these QTLs revealed two genes, Tnik and Phf17, under the miR-212 regulatory QTLs, whose expression levels were significantly correlated with miR-212 expression. We found that miR-212 expression is correlated with cocaine-related behaviour, consistent with a reported role for this miRNA in the control of cocaine consumption. miR-31 is correlated with anxiety and alcohol related behaviours. KEGG pathway analysis of each miRNA’s expression correlates revealed enrichment of pathways including MAP kinase, cancer, long-term potentiation, axonal guidance and WNT signalling.
The BXD reference panel allowed us to establish genetic regulation and characterize biological function of specific miRNAs. QTL analysis enabled detection of genetic loci that regulate the expression of these miRNAs. eQTLs that regulate miRNA abundance are a new mechanism by which genetic variation influences brain and behaviour. Analysis of one of these QTLs revealed a gene, Tnik, which may regulate the expression of a miRNA, a molecular pathway and a behavioural phenotype. Evidence of genetic covariation of miR-212 abundance and cocaine related behaviours is strongly supported by previous functional studies, demonstrating the value of this approach for discovery of new functional roles and downstream processes regulated by miRNA.
Micro-ribonucleic acids (miRNAs) are short single-stranded non-coding RNAs that are involved in the post-transcriptional regulation of genes. MiRNAs are processed from primary transcripts into smaller stem-loop precursor RNAs in the nucleus by the enzyme DROSHA and then further processed by DICER in the cytoplasm into the mature miRNA . Only the mature form of miRNAs cause post-transcriptional gene silencing by imperfect base pairing, as part of a larger molecular complex, with its target sites which in turn can lead to mRNA cleavage or translational repression [2, 3]. A single miRNA can have predicted target sites on hundreds of different mRNAs and thus have wide ranging effects on mRNA expression . Roughly one-third of all genes have at least one predicted miRNA target site . However, the small size of the miRNAs, combined with the complex nature of the miRNA and miRNA target interactions, makes prediction of miRNA target genes difficult .
MiRNAs play a major role in a wide range of functions including development, immune processes, apoptosis and synapse formation [3, 6–11]. MiRNAs have also been shown to play an important role in numerous diseases including cancer, heart disease and mental disorders [12–14]. The associations between miRNAs and these phenotypes and diseases are likely due to the effects of miRNAs on the expression of specific groups of mRNAs via post-transcriptional regulation. Thus, genetic variation that effects miRNA expression will likely affect these phenotypes. Numerous studies have used genetic correlation of gene expression with a variety of phenotypes to discover gene functions in biological processes including brain and behaviour [15, 16]. By exploring patterns in the networks of gene expression genetically correlated with individual miRNAs we may be able to determine each miRNA’s underlying function.
While many miRNAs are widely expressed, the expression of some miRNAs has a spatial specificity [17, 18]. This is the case in adult brain, where miRNAs expression is often exclusive to particular regions  or where there are families of miRNAs that are preferentially expressed in specific brain areas, such as the hippocampus or frontal cortex . This suggests that individual miRNAs may only play a role in defined anatomical areas.
Though much is known about the general function and regulation of miRNAs, much work is still needed to determine the specific functioning of individual miRNAs. A high-throughput functional genomics strategy to investigate the functions of individual miRNAs makes use of the genetic variation across mouse inbred strains and allows one to see how this variation affects miRNA expression. To date, three studies have used this approach to establish variation in miRNA expression across inbred mouse strains, one in the hippocampus  and two in the liver [21, 22]. We previously identified several differentially expressed hippocampal miRNAs in the C57BL/6 J and DBA/2 J strains . Finding the genetic source of this variation and identifying the functional associates requires a genetic reference population that has been broadly characterized for brain and behavioural function.
In this study we further characterized differentially expressed hippocampal miRNAs between the inbred mouse strains C57BL/6 J and DBA/2 J  in a C57BL/6 J x DBA/2 J recombinant inbred panel (BXD). This is the first study to use BXD RI mice to investigate the functioning of specific miRNAs in the brain. The wide array of cumulative data available for these strains allowed us to evaluate associations with a wide range of phenotypes including mRNA expression, neurological and behavioural phenotypes, and perform QTL mapping. This approach allowed us to identify upstream genetic loci that influence functioning of these specific miRNA and downstream gene expression and behavioural correlates which may be influenced by the same loci.
Genetic variation in miRNA expression in the BXD RI lines
The means (± standard deviations) of hippocampal miRNA expression for all the BXD RI strains investigated are listed in Table 1. We replicated the reported expression differences across the C57BL/6 J vs. DBA/2 J strains for the five miRNAs we previously found to be differentially expressed . The η2 values for miRNA gene expression by strain are 0.21, 0.33, 0.28, 0.27 and 0.37 for miR-15b, miR-31, miR-34c, miR-212 and miR-301a, respectively.
MRNA expression correlations of miRNA expression
The correlations among hippocampal miRNA expression and genomewide mRNA expression, independent of the presence of a predicted miRNA target site, are summarized in Table 2 (for a complete list of nominally significant genes (p-value < 0.05) for each miRNA see Additional file 1: Table S1). We controlled the false detection rate for each gene list individually using the Benjamini and Hochberg method, q < 0.2 . Only a single correlation remained significant following FDR multiple testing correction, which was the correlation between miR-15b and probe set 1437110_at (p-value = 2.6e-6, q = 0.04). This probe set corresponds to the gene 2810474O19Rik.
Χ2 analysis was used to determine whether those genes with predicted miRNA binding sites were overrepresented within the list of nominally significant correlations between miRNA and mRNA. We used target sites predicted by the MIRANDA algorithm as this is one of the most widely used and has a higher sensitivity rate than the other commonly used algorithms . As increases in miRNA expression can degrade the corresponding target mRNA, we conducted these analyses first using just the significant negative correlations, then using all significant correlations (summarized on Table 2). None of the miRNAs showed an overrepresentation of genes with predicted miRNA binding sites when either just the negative correlations or all correlations were used.
KEGG pathway analysis
We used the list of all genes for which there were nominally significant correlations between miRNA vs. mRNA expression, independent of the presence of a miRNA target site, to conduct KEGG pathway analyses for each miRNA. We found a number of pathways that were significantly enhanced for these miRNAs (summarized on Table 3), including numerous neurobehavioural pathways. In particular, we found that genes correlated with miR-15b expression were enriched in a number of pathways related to neuronal plasticity, including axon guidance, long-term potentiation, long-term depression and regulation of the actin cytoskeleton.
QTL analysis of miRNA expression
We conducted quantitative trait loci analysis of miRNA expression for each miRNA using R/QTL. We found a total of five suggestive QTLs which are summarized on Table 4. In order to further investigate what genes may be responsible for these QTLs, we conducted correlations of miRNA vs. mRNA gene expression for those genes within the 1.5 LOD confidence intervals for each of these peaks. We controlled the false detection rate for each gene list individually using the Benjamini and Hochberg method, q < 0.2 . Two of these correlations survived multiple testing correction, 1429870_at (gene Tnik) and 1438412_at (gene Phf17), both lying underneath the expression QTL (eQTL) for miR-212 on chromosome 3 (see Figure 1).
We examined whether there were any genes under these eQTL peaks that contained predicted miRNA binding sites. There were genes with predicted miRNA target sites under all of the eQTL peaks, except the miR-301a peak on chr11 (summarized in Additional file 2: Table S2). The gene Acad9, which codes for a mitochondrial Acyl-CoA dehydrogenase, had two predicted miR-212 target sites and was negatively correlated with miR-212 expression (r = -.446, p = .029). Neither of the genes, Phf17 and Tnik, have a predicted miR-212 target site.
Correlations of miRNA expression and phenotype
The correlations of miRNA expression and phenotype revealed a total of 69 nominally significant correlations (p < 0.05), summarized in Additional file 3: Table S3. Two of these correlations survived FDR multiple testing correction (q < 0.2), miR-31 expression and both blood ethanol concentration (following 2.25 mg/kg ethanol) (r = 0.7, p = 0.001) and percent time spent in the centre area of the open field test (r = 0.68, p = 0.002) (measures from ). In order to determine if the number of significant correlations within a behavioural category was overrepresented, we conducted Fishers Exact Tests which showed that behaviour traits in general were overrepresented with miR-31 expression (p < 0.05) and cocaine related traits were overrepresented with both miR-34c expression (p < 0.05) and miR-212 expression (p < 0.1) (see Table 5).
BXD-RI lines have been successfully used to characterize the roles that genetic variation plays in various molecular and functional phenotypes, including gene expression, and thus in complex behaviour and disease [16, 26, 27]. The genetic and phenotypic variation across the RI strains, combined with extensive collections of cumulative data available in extensive web-based databases including gene expression data, genetic mapping panels and phenotypic data  underlies the utility of these reference strains in such characterization.
One of the greatest advantages of using BXD RI, and similar RI populations, is that existing dense genotypic data allows variation of any quantitative phenotype to be readily mapped using QTL mapping approaches, thus allowing for the discovery of putative local (cis) and distant (trans) genetic loci governing this variation. We were able to find five suggestive distant QTLs for the miRNAs investigated. Of particular interest was the QTL for miR-212 expression on chromosome 3, within which we found two expression probe sets that were significantly correlated with miR-212 expression. Phf17, PHD finger protein17 (also know as Jade1), has also been shown to be involved in WNT pathway signalling , particularly playing a role in anteroposterior axis development . Tnik, TRAF2 and NCK-interacting protein kinase, is a member of the Ste20 group of kinases, known to be regulators of MAP kinase cascades . It has previously been shown to be an activator of WNT target genes  and regulation of the cytoskeleton . KEGG pathway analysis for co-expressed transcripts of miR-212 supports a role for miR-212 in the regulation of the cytoskeleton and of MAPK kinase pathways, functions that overlap with those of Tnik and suggesting that the miR-212 and Tnik correlation may underlie a real functional relationship between these genes.
It should be noted that neither of the genes corresponding to these probe sets, Phf17 and Tnik, have a predicted miR-212 binding site or any RNA genes within their transcripts. The existence of such a site is not required to posit a mechanism by which genetic variation in that a functional process can indirectly regulate the abundance of miR-212. Further studies are required to determine which genetic variant potentially underlies this correlation, and what genes and mechanisms this variant acts upon to indirectly regulate miR-212 expression.
We found that the gene Acad9, which lies under the QTL for miR-212 expression on chromosome 3, had two predicted miR-212 target sites. Furthermore its expression was also negatively correlated with miR-212 expression. This suggests that there is the potential for Acad9 to both indirectly regulate miR-212 expression and in turn be directly regulated by miR-212 expression. While we do not expect to find such a mechanism for every QTL locus, this is an additional mechanism by which genetic variation could influence gene expression.
While we did find five suggestive distant QTLs, we did not find any local QTLs, indicating that polymorphisms in the miRNA we evaluated do not directly influence their abundance in the hippocampus. Two previous studies investigating liver miRNA expression found a local QTL for miR-31 in BXD and F2 panels [21, 22]. The failure to replicate this local QTL in the present study may be due to the use of a different tissue (hippocampus) in our study. A study conducted in human samples looking for eQTLs across numerous tissues found that while 30% of all eQTLs were shared across the three tissues investigated, 29% were tissue-specific .
We found significant correlations between both miR-34c and miR-212 expression and cocaine-related behaviours. This is particularly interesting for miR-212, which has previously been shown to control cocaine intake . In particular, it was shown that striatal miR-212 expression is increased following extended cocaine use in rats and that increases in striatal miR-212 expression leads to decreases in cocaine intake following extended access conditions. Furthermore, miR-212 is known to be regulated by MeCP2, an important regulator of neuroplasticity, and to affect BDNF expression and neuroplasticity in postmitotic neurons . As cocaine addiction is widely believed to result in changes in neurocircuitry , miR-212 is an excellent candidate for susceptibility to cocaine addiction, further supported by a previous finding that miR-212 expression affects dendrite growth and arborisation . Our present result augments this finding by demonstrating that genetic polymorphisms can cause phenotypic variation in this process. The detection of this relationship also demonstrates the utility of a systems genetics strategy for the discovery of specific molecular and functional roles of miRNAs.
We found significant correlations (q < 0.2) of miR-31 expression with both blood ethanol concentration and a measure of anxiety (percent time spent in the centre area of the open field test). We further found an overrepresentation of traits in general behaviour with miR-31 expression, with a particular number from open field and light-dark box measures, suggesting that miR-31 may play a role in anxiety. These are the first indications of miR-31 potentially being involved in alcohol or anxiety related traits. A number of miRNAs have previously been shown to be altered in the brain of alcoholics, including miR-15b, miR-34c and miR-301a . This suggests that miR-31 should be investigated further for a role in anxiety and susceptibility to alcohol use, but that it may not be associated with long term alcohol exposure.
The WNT signalling pathway has previously been associated with adult neurogenesis in the hippocampus . Both the MAPK and WNT signalling pathways were enhanced within the KEGG pathway analysis for miR-31 expression. These pathways have been associated with miR-31 in a previous KEGG pathway analysis in human tumour cells . Over expression of miR-31 increased Wnt-5a expression, which adds further support for the involvement of miR-31 in the WNT signalling pathway . Together this suggests that miR-31 may possibly play a role in adult neurogenesis via the WNT signalling pathway.
MiR-301a plays a role in wide range of cancers including lung cancer , breast cancer , pancreatic cancer . KEGG analysis revealed that miR-301a was associated with colorectal cancer. The expression of miR-301a is altered in p53-deficient mice, a known tumor suppressor gene, further supporting its involvement in cancer .
A recent study has shown that miR-34c plays a role in anxiety at the level of the amygdala . More specifically, this gene was upregulated following acute and chronic stress, and lenti-virus mediated overexpression of miR-34c in the amygdala induced anxiolytic behaviour after challenge. We found the miR-34c expression was positively correlated with open arm duration (r = 0.43, p = 0.03), a measure of anxiety. Additionally, miR-34c has been shown to reduce the cellular response to corticotrophin releasing factor receptor type 1 (CRFR1), possibly acting via a miR-34c target site on the CRFR1 mRNA . Together these findings suggest a role for miR-34c in regulating the central stress response. This gene has also been shown to be elevated in the hippocampus of Alzheimer’s disease patients and the corresponding mouse models, and overexpressing miR-34c leads to memory impairment . Together this data suggests a role for this gene in the mechanisms of anxiety and memory which should be further investigated.
MiR-34c has also been shown to be involved in various cancers [49–51]. More specifically, miR-34c is thought to act as a tumor suppressor as part of a negative feedback loop including Myc and Mapkapk5, part of the MAPK signalling pathway . We found that the MAPK signalling pathway is enhanced for miR-34c, and while miR-34c was not significantly associated with Mapkapk5 expression (r = -0.24, p = 0.27), it was with Mapkapk3 (r = 0.55, p = 0.009).
Our KEGG pathway analysis for miR-15b suggests it plays a role in long-term potentiation, long-term depression and axon guidance. A study investigating the localization of miRNAs in sympathetic neurons revealed that miR-15b is more highly abundant in the distal axons compared to the cell bodies . Taken together, these findings suggest a role for miR-15b in neuronal plasticity. Expression studies have also linked miR-15b expression with various cancers [46, 54, 55], including pancreatic cancer . Similarly, we found an enhanced KEGG pathway for pancreatic cancer for this gene.
The only correlation between miRNA and mRNA expression that had a genome-wide significance level was between miR-15b and 2810474O19Rik expression. 2810474O19Rik, which has no predicted miR-15b binding sites, is thought to play a role in development, being expressed in both gonadal  and preimplantation mouse development . This gene also has been suggested to play a role in cell potency, more specifically in pluripotent cell identity via an interaction with Oct4. There is limited evidence that miR-15b plays a role in cell potency, with it being expressed in multipotent cells during osteogenic differentiation . The expression of miR-15b was upregulated in the umbilical vein  and human placenta , but there is no known link between miR-15b and development.
Our study made use of existing mRNA data collected in independent mice. While this enables one to clearly ascribe correlation to genetic factors, each sample came from environmentally distinct mice. Thus, one may conclude that environmental variation in mRNA across these two population samples was sufficient to exceed the genetic variation accounted for by the loci we detected. Nonetheless, relevant co-expression was detected and future studies in which samples are collected from a single population for phenotype, miRNA and gene expression are warranted. This may be done in a genetic reference population to exploit the breadth of existing data or a large experimental cross or mapping population such as the Diversity Outbred to improve mapping power and precision.
In our study we investigated the potential effects of miRNAs on mRNA by correlating their gene expression. For this approach to be successful, we need a significant percentage of miRNAs to regulate their target genes via RNA degradations rather than by blocking translation, as blocking translation would not necessarily change the mRNA levels. Numerous groups have successfully used this approach to show that miRNAs do reduce the expression of a significant fraction of their targets [63–66]. This approach thus allows us to investigate the effects of miRNA on their potential gene targets for a great number of genes.
If a significant fraction of predicted miRNA target sites are real and if miRNAs commonly cause RNA degradation of their target genes then we should see an over-representation of miRNA vs. mRNA correlations for genes with predicted target sites. We failed to find this for the miRNAs that we investigated (both all correlated genes and just the negatively correlated genes). This generally held true in another study investigating liver miRNA expression in inbred and BXD RI strains . It is possible that there is little degradation of mRNA transcripts by miRNA in mammals, and instead miRNAs block translation of the targeted mRNA. The difficulty in accurate prediction of miRNA binding sites due to the small size of the recognition sequences could also account for this result . An alternative interpretation is that additional cellular processes, possibility including indirect mechanisms of miRNA gene regulation, act in an indirect fashion to modulate or mask the regulatory effects of miRNAs in vivo.
Systems genetic analyses make use of naturally occurring genetic polymorphisms to simultaneously map sources of genetic variation and identify relations among biological entities across biological scale (in this case miRNA, mRNA and behaviour). Genetic analysis of miRNA expression and co-expression in the BXD recombinant inbred panel, takes advantage of the extensive cumulative data available for this panel, to further characterize five miRNAs that were previously shown to be differentially expressed in the BXD progenitor strains . This is the first study of brain miRNA abundance in the BXD genetic reference panel, and it revealed that a well-characterized miRNA to behaviour relationship is subject to genetic control.
Using this approach, we conducted QTL analysis that enabled the detection of genetic loci that regulate the expression of these miRNAs. eQTLs that regulate miRNA abundance are a new mechanism by which genetic variation influences brain and behaviour. Analysis of one of these QTLs revealed a gene, Tnik, which may regulate the expression of a miRNA, a molecular pathway and a behavioural phenotype. Furthermore, we found evidence of genetic covariation of miR-212 abundance and cocaine related behaviours that is strongly supported by previous functional studies, demonstrating the value of this approach for discovery of new functional roles and downstream processes regulated by miRNA.
In summary, systems genetic analysis of miRNA abundance is a promising approach to discovery of the functional significance or particular miRNA in phenotypic variation and disease. With a simple profiling study across the original BXD RI strains, robust associations of miRNA to gene expression and behaviour can be detected. Simultaneous estimation of gene expression and miRNA abundance in the full BXD RI panel or high-power, high precision populations such as large experimental crosses and advanced mouse populations recently developed by the Complex Trait Consortium may yield many more such associations to bridge the gap between molecular discovery of miRNA and functional biology.
Hippocampal tissue was collected from 24 strains of C57BL/6 J x DBA/2 J recombinant inbred mice (BXD): BXD1/TyJ, BXD5/TyJ, BXD6/TyJ, BXD8/TyJ, BXD9/TyJ, BXD11/TyJ, BXD12/TyJ, BXD16/TyJ, BXD18/TyJ, BXD19/TyJ, BXD21/TyJ, BXD22/TyJ, BXD23/TyJ, BXD24/TyJ, BXD27/TyJ, BXD28/TyJ, BXD30/TyJ, BXD31/TyJ, BXD32/TyJ, BXD33/TyJ, BXD34/TyJ, BXD39/TyJ, BXD40/TyJ, BXD42/TyJ (n = 4 for all BXD strains except BXD30/TyJ where n = 2). The colony of BXD RI strains was maintained at the Institute of Psychiatry using original stocks purchased fromsbl;lll The Jackson Laboratory (Bar Harbor, ME, USA). All animals were males aged between 100-110 days when killed. After cervical dislocation, bilateral hippocampi were dissected in their entirety from fresh brains within two minutes from the time of death. Any connecting tissue was trimmed off and the hippocampi were immediately snap frozen on dry ice and stored at -800C. All housing and experimental procedures were carried out in accordance with the UK Home Office Animals (Scientific Procedures) Act 1986 under License PPL No. 70/5113.
In order to obtain reliable levels of miRNA, total RNA was isolated from hippocampus using a miRNA isolation kit (Abmbion, Life Technologies, UK). The total RNA concentration was determined using a nanodrop ND-1000 (Thermo Fisher Scientific).
Taqman real-time polymerase chain reactions (RT-PCR)
The mature miRNA expression for the five miRNAs miR-15b, miR-31, miR-34c, miR-212, and miR-301a were quantified in the individual BXD animals using Taqman RT-PCR assays (Applied Biosystems, Life Technologies, Foster City, CA, USA). The reverse transcription (RT) reactions were conducted using the miRNA reverse transcription primer specific for each assay, using the standard TaqMan® micro RNA Reverse Transcription protocol (http://www.appliedbiosystems.com, document #: 4364031). The RT-PCR reactions were performed in triplicate using 0.5 μl of 20x PCR Probe/Primer Mix, 1.5 μl of product from the RT reaction (diluted 1:10), 5 μl of 2x TaqMan Master Mix (No UNG) and 3 μls nuclease-free water. A sample minus reverse transcription buffer was used as a negative control for the RT-PCR reactions. (All of the negative controls used failed to reach threshold by 45 cycles).
Reactions were run on a 7900HT Fast real-Time PCR System (Applied Biosystems, Foster City, CA, USA) in 384 well format. RNU19, miR-9 and miR-99a were used as controls for each of the test assays investigated in the individual samples. These assays were chosen as controls because their expression did not differ across strain and they had low variability within strain.
Taqman RT-PCR statistical analysis
Relative expression was calculated using the standard 2 relative expression method in Microsoft Excel. The relative expression of the individual samples was normalized to the geometric mean of RNU19, miR-9 and miR-99a (previously shown not to differ between C57BL/6 J and DBA 2 J, ). The η2 values were calculated from the genetic and environmental (error-term) sums of squares derived from a one-way ANOVA using the open source statistical program R.
Correlations of behaviour and miRNA expression
Correlations were conducted between miRNA expression and the phenotypic measures from the behavioural batteries previously described in  (n = 40) and  (n = 243). We limited our analysis to those measures for which we had the corresponding phenotype and miRNA expression data in at least 18 BXD RI strains. Correlations were conducted using PASW Statistics 18, Release Version 18.0.0 (SPSS Inc., 2009, Chicago, IL, USA). False Discovery Rate control  was used to account for multiple testing (q < 0.2). For each miRNA, we conducted Fisher’s exact tests (one-sided) to determine if the number of significantly correlated behavioural traits within a given behavioural category was greater than that predicted by chance.
Correlations of mRNA and miRNA expression
Lists of genes with MIRANDA predicted miRNA sites  were downloaded from the online miRNA database, miRBase (http://microrna.sanger.ac.uk/). Hippocampal mRNA expression for the BXD RI mice was obtained from the Hippocampus Expression Consortium . The hippocampus, excluding most of the subiculum, from two to three animals was dissected and pooled for hybridization to a single Affymetrix M430 2.0 array. Raw microarray data was transformed using the PDNN, MAS5 and RMA methods. to 2z + 8, thus yielding a data set with a standard deviation of 2 and an overall mean of 8.
The expression data and predicted target site information were merged and correlations between miRNA and mRNA expression were conducted using the open source statistical program R . False discovery rate  was used to account for multiple testing. KEGG pathway analysis was conducted using the lists of nominally significant genes for each miRNA using WebGestalt (http://bioinfo.vanderbilt.edu/webgestalt). Chi squared tests for independence (two-sided without Yates correction) were calculated to see if there was a relationship between the presence of a target site and a significant correlation between mRNA and miRNA expression (p <0.05).
QTL analysis of miRNA expression
QTL analysis was conducted for miRNA expression for those miRNAs investigated in the BXD RI strains (miR-15b, miR-31, miR-34c, miR-212, and miR- 301a). These analyses were done using the R-QTL package in the open source statistical program R as previously described in . Confidence intervals (1 and 1.5 LOD score drops) were calculated for any suggestive QTL peaks. A series of correlations between miRNA expression and mRNA expression were conducted for all of the genes underlying all the suggestive eQTL peaks using PASW Statistics 18, Release Version 18.0.0 (SPSS Inc., 2009, Chicago, IL, USA). False discovery rate control  was used to account for multiple testing (q < 0.2).
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297.
Yekta S, Shih IH, Bartel DP: MicroRNA-directed cleavage of HOXB8 mRNA. Science. 2004, 304 (5670): 594-596.
Kloosterman WP, Plasterk RH: The diverse functions of microRNAs in animal development and disease. Dev Cell. 2006, 11 (4): 441-450.
Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9 (2): 102-114.
Witkos TM, Koscianska E, Krzyzosiak WJ: Practical Aspects of microRNA Target Prediction. Curr Mol Med. 2011, 11 (2): 93-109.
Sayed D, Abdellatif M: MicroRNAs in development and disease. Physiol Rev. 2011, 91 (3): 827-887.
Chen CZ, Li L, Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004, 303 (5654): 83-86.
Zhang B, Wang Q, Pan X: MicroRNAs and their regulatory roles in animals and plants. J Cell Physiol. 2007, 210 (2): 279-289.
Rodriguez A, Vigorito E, Clare S, Warren MV, Couttet P, Soond DR, van Dongen S, Grocock RJ, Das PP, Miska EA, Vetrie D, Okkenhaug K, Enright AJ, Dougan G, Turner M, Bradley A: Requirement of bic/microRNA-155 for normal immune function. Science. 2007, 316 (5824): 608-611.
Thai TH, Calado DP, Casola S, Ansel KM, Xiao C, Xue Y, Murphy A, Frendewey D, Valenzuela D, Kutok JL, Schmidt-Supprian M, Rajewsky N, Yancopoulos G, Rao A, Rajewsky K: Regulation of the germinal centre response by microRNA-155. Science. 2007, 316 (5824): 604-608.
Fiore R, Schratt G: MicroRNAs in synapse development: tiny molecules to remember. Expert Opin Biol Ther. 2007, 7 (12): 1823-1831.
Abelson JF, Kwan KY, O'Roak BJ, Baek DY, Stillman AA, Morgan TM, Mathews CA, Pauls DL, Rasin MR, Gunel M, Davis NR, Ercan-Sencicek AG, Guez DH, Spertus JA, Leckman JF, Dure LS, Kurlan R, Singer HS, Gilbert DL, Farhi A, Louvi A, Lifton RP, Sestan N, State MW: Sequence variants in SLITRK1 are associated with Tourette's syndrome. Science. 2005, 310 (5746): 317-320.
Gregory RI, Shiekhattar R: MicroRNA biogenesis and cancer. Cancer Res. 2005, 65 (9): 3509-3512.
van Rooij E, Olson EN: MicroRNAs: powerful new regulators of heart disease and provocative therapeutic targets. J Clin Invest. 2007, 117 (9): 2369-2376.
Chesler EJ, Wang J, Lu L, Qu Y, Manly KF, Williams RW: Genetic correlates of gene expression in recombinant inbred strains: a relational model system to explore neurobehavioral phenotypes. Neuroinformatics. 2003, 1 (4): 343-357.
Churchill GA: Recombinant inbred strain panels: a tool for systems genetics. Physiol Genomics. 2007, 31 (2): 174-175.
Juhila J, Sipilä T, Icay K, Nicorici D, Ellonen P, Kallio A, Korpelainen E, Greco D, Hovatta I: MicroRNA expression profiling reveals MiRNA families regulating specific biological pathways in mouse frontal cortex and hippocampus. PLoS One. 2011, 6 (6): e21495-
Liu S, Zhang L, Li Q, Zhao P, Duan J, Cheng D, Xiang Z, Xia Q: MicroRNA expression profiling during the life cycle of the silkworm (Bombyx mori). BMC Genomics. 2009, 10: 455-
Olsen L, Klausen M, Helboe L, Nielsen FC, Werge T: MicroRNAs show mutually exclusive expression patterns in the brain of adult male rats. PLoS One. 2009, 4 (10): e7225-
Parsons MJ, Grimm C, Nietfeld W, Lehrech H, Schalkwyk LC: Differences in miRNA expression across common inbred strains using Taqman RT-PCR. Mamm Genome. 2008, 19 (7-8): 552-560.
Zhao E, Keller MP, Rabaglia ME, Oler AT, Stapleton DS, Schueler KL, Neto EC, Moon JY, Wang P, Wang IM, Lum PY, Ivanovska I, Cleary M, Greenawalt D, Tsang J, Choi YJ, Kleinhanz R, Shang J, Zhou YP, Howard AD, Zhang BB, Kendziorski C, Thornberry NA, Yandell BS, Schadt EE, Attie AD: Obesity and genetics regulate microRNAs in islets, liver, and adipose of diabetic mice. Mamm Genome. 2009, 20 (8): 476-485.
Gatti DM, Lu L, Williams RW, Sun W, Wright FA, Threadgill DW, Rusyn I: MicroRNA expression in the livers of inbred mice. Mutat Res. 2011, 714 (1–2): 126-133.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
Alexiou P, Maragkakis M, Papadopoulos GL, Reczko M, Hatzigeorgiou AG: Lost in translation: an assessment and perspective for computational microRNA target identification. Bioinformatics. 2009, 25 (23): 3049-3055.
Philip VM, Duvvuru S, Gomero B, Ansah TA, Blaha CD, Cook MN, Hamre KM, Lariviere WR, Matthews DB, Mittleman G, Goldowitz D, Chesler EJ: High-throughput behavioral phenotyping in the expanded panel of BXD recombinant inbred strains. Genes Brain Behav. 2010, 9 (2): 129-159.
Chesler EJ, Lu L, Shou S, Qu Y, Gu J, Wang J, Hsu HC, Mountz JD, Baldwin NE, Langston MA, Threadgill DW, Manly KF, Williams RW: Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function. Nat Genet. 2005, 37 (3): 233-242.
Bystrykh L, Weersing E, Dontje B, Sutton S, Pletcher MT, Wiltshire T, Su AI, Vellenga E, Wang J, Manly KF, Lu L, Chesler EJ, Alberts R, Jansen RC, Williams RW, Cooke MP, de Haan G: Uncovering regulatory pathways that affect hematopoietic stem cell function using 'genetical genomics'. Nat Genet. 2005, 37 (3): 225-232.
Chesler EJ, Lu L, Wang J, Williams RW, Manly KF: WebQTL: rapid exploratory analysis of gene expression and genetic networks for brain and behaviour. Nat Neurosci. 2004, 7 (5): 485-486.
Chitalia VC, Foy RL, Bachschmid MM, Zeng L, Panchenko MV, Zhou MI, Bharti A, Seldin DC, Lecker SH, Dominguez I, Cohen HT: Jade-1 inhibits Wnt signalling by ubiquitylating beta-catenin and mediates Wnt pathway inhibition by pVHL. Nat Cell Biol. 2008, 10 (10): 1208-1216.
Tzouanacou E, Tweedie S, Wilson V: Identification of Jade1, a gene encoding a PHD zinc finger protein, in a gene trap mutagenesis screen for genes involved in anteroposterior axis development. Mol Cell Biol. 2003, 23 (23): 8553-8562.
Dan I, Watanabe NM, Kusumi A: The Ste20 group kinases as regulators of MAP kinase cascades. Trends Cell Biol. 2001, 11 (5): 220-230.
Mahmoudi T, Li VS, Ng SS, Taouatas N, Vries RG, Mohammed S, Heck AJ, Clevers H: The kinase TNIK is an essential activator of Wnt target genes. EMBO J. 2009, 28 (21): 3329-3340.
Fu CA, Shen M, Huang BC, Lasaga J, Payan DG, Luo Y: TNIK, a novel member of the germinal centre kinase family that activates the c-Jun N-terminal kinase pathway and regulates the cytoskeleton. J Biol Chem. 1999, 22 (43): 30729-30737. 274
Nica AC, Parts L, Glass D, Nisbet J, Barrett A, Sekowska M, Travers M, Potter S, Grundberg E, Small K, Hedman AK, Bataille V, Tzenova Bell J, Surdulescu G, Dimas AS, Ingle C, Nestle FO, di Meglio P, Min JL, Wilk A, Hammond CJ, Hassanali N, Yang TP, Montgomery SB, O'Rahilly S, Lindgren CM, Zondervan KT, Soranzo N, Barroso I, Durbin R, Ahmadi K, Deloukas P, McCarthy MI, Dermitzakis ET, Spector TD: MuTHER Consortium. The architecture of gene regulatory variation across multiple human tissues: the MuTHER study. PLoS Genet. 2011, 7 (2): e1002003-
Hollander JA, Im HI, Amelio AL, Kocerha J, Bali P, Lu Q, Willoughby D, Wahlestedt C, Conkright MD, Kenny PJ: Striatal microRNA controls cocaine intake through CREB signalling. Nature. 2010, 466 (7303).
Im HI, Hollander JA, Bali P, Kenny PJ: MeCP2 controls BDNF expression and cocaine intake through homeostatic interactions with microRNA-212. Nat Neurosci. 2010, 13 (9): 1120-1127.
Uys JD, Reissner KJ: Glutamatergic neuroplasticity in cocaine addiction. Prog Mol Biol Transl Sci. 2011, 98: 367-400.
Magill ST, Cambronne XA, Luikart BW, Lioy DT, Leighton BH, Westbrook GL, Mandel G, Goodman RH: microRNA-132 regulates dendritic growth and arborization of newborn neurons in the adult hippocampus. Proc Natl Acad Sci USA. 2010, 107 (47): 20382-20387.
Lewohl JM, Nunez YO, Dodd PR, Tiwari GR, Harris RA, Mayfield RD: Up-regulation of MicroRNAs in brain of human alcoholics. Alcohol Clin Exp Res. 2011, 35 (11): 1928-1937.
Lie DC, Colamarino SA, Song HJ, Désiré L, Mira H, Consiglio A, Lein ES, Jessberger S, Lansford H, Dearie AR, Gage FH: Wnt signalling regulates adult hippocampal neurogenesis. Nature. 2005, 437 (7063): 1370-1375.
Necela BM, Carr JM, Asmann YW, Thompson EA: Differential expression of microRNAs in tumors from chronically inflamed or genetic (APC(Min/+)) models of colon cancer. PLoS One. 2011, 6 (4): e18501-
Xi S, Yang M, Tao Y, Xu H, Shan J, Inchauste S, Zhang M, Mercedes L, Hong JA, Rao M, Schrump DS: Cigarette smoke induces C/EBP-β-mediated activation of miR-31 in normal human respiratory epithelia and lung cancer cells. PLoS One. 2010, 5 (10): e13764-
Silva J, García V, Zaballos Á, Provencio M, Lombardía L, Almonacid L, García JM, Domínguez G, Peña C, Diaz R, Herrera M, Varela A, Bonilla F: Vesicle-related microRNAs in plasma of nonsmall cell lung cancer patients and correlation with survival. Eur Respir J. 2011, 37 (3): 617-623.
Shi W, Gerster K, Alajez NM, Tsang J, Waldron L, Pintilie M, Hui AB, Sykes J, P'ng C, Miller N, McCready D, Fyles A, Liu FF: MicroRNA-301 mediates proliferation and invasion in human breast cancer. Cancer Res. 2011, 71 (8): 2926-2937.
Lee EJ, Gusev Y, Jiang J, Nuovo GJ, Lerner MR, Frankel WL, Morgan DL, Postier RG, Brackett DJ, Schmittgen TD: Expression profiling identifies microRNA signature in pancreatic cancer. Int J Cancer. 2007, 120 (5): 1046-1054.
Hosako H, Martin GS, Barrier M, Chen YA, Ivanov IV, Mirkes PE: Gene and microRNA expression in p53-deficient day 8.5 mouse embryos. Birth Defects Res A Clin Mol Terato. 2009, 85 (6): 546-555.
Haramati S, Navon I, Issler O, Ezra-Nevo G, Gil S, Zwang R, Hornstein E, Chen A: MicroRNA as repressors of stress-induced anxiety: the case of amygdalar miR-34. J Neurosci. 2011, 31 (40): 14191-14203.
Zovoilis A, Agbemenyah HY, Agis-Balboa RC, Stilling RM, Edbauer D, Rao P, Farinelli L, Delalle I, Schmitt A, Falkai P, Bahari-Javan S, Burkhardt S, Sananbenesi F, Fischer A: microRNA-34c is a novel target to treat dementias. EMBO J. 2011, 30 (20): 4299-4308.
Suzuki H, Yamamoto E, Nojima M, Kai M, Yamano HO, Yoshikawa K, Kimura T, Kudo T, Harada E, Sugai T, Takamaru H, Niinuma T, Maruyama R, Yamamoto H, Tokino T, Imai K, Toyota M, Shinomura Y: Methylation-associated silencing of microRNA-34b/c in gastric cancer and its involvement in an epigenetic field defect. Carcinogenesis. 2010, 31 (12): 2066-2073.
Hagman Z, Larne O, Edsjö A, Bjartell A, Ehrnström RA, Ulmert D, Lilja H, Ceder Y: miR-34c is downregulated in prostate cancer and exerts tumor suppressive functions. Int J Cancer. 2010, 127 (12): 2768-2776.
Wang Z, Chen Z, Gao Y, Li N, Li B, Tan F, Tan X, Lu N, Sun Y, Sun J, Sun N, He J: DNA hypermethylation of microRNA-34b/c has prognostic value for stage I non-small cell lung cancer. Cancer Biol Ther. 2011, 11 (5): 490-496.
Kress TR, Cannell IG, Brenkman AB, Samans B, Gaestel M, Roepman P, Burgering BM, Bushell M, Rosenwald A, Eilers M: The MK5/PRAK kinase and Myc form a negative feedback loop that is disrupted during colorectal tumorigenesis. Mol Cell. 2011, 41 (4): 445-457.
Natera-Naranjo O, Aschrafi A, Gioio AE, Kaplan BB: Identification and quantitative analyses of microRNAs located in the distal axons of sympathetic neurons. RNA. 2010, 16 (8): 1516-1529.
Chung GE, Yoon JH, Myung SJ, Lee JH, Lee SH, Lee SM, Kim SJ, Hwang SY, Lee HS, Kim CY: High expression of microRNA-15b predicts a low risk of tumor recurrence following curative resection of hepatocellular carcinoma. Oncol Rep. 2010, 23 (1): 113-119.
Satzger I, Mattern A, Kuettler U, Weinspach D, Voelker B, Kapp A, Gutzmer R: MicroRNA-15b represents an independent prognostic parameter and is correlated with tumor cell proliferation and apoptosis in malignant melanoma. Int J Cancer. 2010, 126 (11): 2553-2562.
Zhang Y, Li M, Wang H, Fisher WE, Lin PH, Yao Q, Chen C: Profiling of 95 microRNAs in pancreatic cancer cell lines and surgical specimens by real-time PCR analysis. World J Surg. 2009, 33 (4): 698-709.
Hudson QJ, Smith CA, Sinclair AH: Conserved expression of a novel gene during gonadal development. Dev Dyn. 2005, 233 (3): 1083-1090.
Ko MS, Kitchen JR, Wang X, Threat TA, Wang X, Hasegawa A, Sun T, Grahovac MJ, Kargul GJ, Lim MK, Cui Y, Sano Y, Tanaka T, Liang Y, Mason S, Paonessa PD, Sauls AD, DePalma GE, Sharara R, Rowe LB, Eppig J, Morrell C, Doi H: Large-scale cDNA analysis reveals phased gene expression patterns during preimplantation mouse development. Development. 2000, 127 (8): 1737-1749.
van den Berg DL, Snoek T, Mullin NP, Yates A, Bezstarosti K, Demmers J, Chambers I, Poot RA: An Oct4-centered protein interaction network in embryonic stem cells. Cell Stem Cell. 2010, 6 (4): 369-381.
Gao J, Yang T, Han J, Yan K, Qiu X, Zhou Y, Fan Q, Ma B: MicroRNA expression during osteogenic differentiation of human multipotent mesenchymal stromal cells from bone marrow. J Cell Biochem. 2011, 112 (7): 1844-1856.
Liu XD, Wu X, Yin YL, Liu YQ, Geng MM, Yang HS, Blachier F, Wu GY: Effects of dietary L: -arginine or N-carbamylglutamate supplementation during late gestation of sows on the miR-15b/16, miR-221/222, VEGFA and eNOS expression in umbilical vein. Amino Acids. 2011, -in press
Mayor-Lynn K, Toloubeydokhti T, Cruz AC, Chegini N: Expression profile of microRNAs and mRNAs in human placentas from pregnancies complicated by preeclampsia and preterm labor. Reprod Sci. 2011, 18 (1): 46-56.
Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van Dongen S, Inoue K, Enright AJ, Schier AF: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 312 (5770): 75-79.
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005, 433 (7027): 769-773.
Selbach M, Schwanhäusser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 4 (7209): 58-63. 455
Shahab SW, Matyunina LV, Mezencev R, Walker LD, Bowen NJ, Benigno BB, McDonald JF: Evidence for the complexity of microRNA-mediated regulation in ovarian cancer: a systems approach. PLoS One. 2011, 6 (7): e22508-
Lad HV, Liu L, Paya-Cano JL, Parsons MJ, Kember R, Fernandes C, Schalkwyk LC: Behavioural battery testing: evaluation and behavioural outcomes in 8 inbred mouse strains. Physiol Behav. 2010, 99 (3): 301-316.
Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008, 36 (Database issue): D149-D153.
Overall RW, Kempermann G, Peirce J, Lu L, Goldowitz D, Gage FH, Goodwin S, Smit AB, Airey DC, Rosen GD, Schalkwyk LC, Sutter TR, Nowakowski RS, Whatley S, Williams RW: Genetics of the hippocampal transcriptome in mouse: a systematic survey and online neurogenomics resource. Front Neurosci. 2009, 3: 55-
Ihaka R, Gentleman R: R: a language for data analysis and graphics. J Comput Graph Stat. 1996, 5: 299-314.
Lad HV, Liu L, Payá-Cano JL, Fernandes C, Schalkwyk LC: Quantitative traits for the tail suspension test: automation, optimization, and BXD RI mapping. Mamm Genome. 2007, 18 (6–7): 482-491.
This work was supported by a grant from the European New and Emerging Science and Technology program, contract NEST 028594 (LS and HL).
The authors declare that they have no competing interests.
MJP conducted the real-time PCRs, the correlational analysis for miRNA vs. mRNAs, the KEGG pathway analysis and drafted the manuscript. LS carried out the QTL analysis and helped conceive the study. LL, JLPC and CF assisted with the RNA extractions and the behavioural analysis from . EJC and VMP assisted with the analysis and interpretation of behavioural data from . CG helped with the real-time PCR analysis. WN and HL assisted with the conceiving and planning of the project. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1. Summary of all significant correlations of miRNA expression and mRNA expression (p-values < 0.05). (DOC 3 MB)
Additional file 3: Table S3. Summary of all significant correlations of miRNA expression and phenotype measures (p-values < 0.05). (DOC 102 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Parsons, M.J., Grimm, C., Paya-Cano, J.L. et al. Genetic variation in hippocampal microRNA expression differences in C57BL/6 J X DBA/2 J (BXD) recombinant inbred mouse strains. BMC Genomics 13, 476 (2012). https://doi.org/10.1186/1471-2164-13-476