- Research article
- Open Access
Prediction of microRNAs affecting mRNA expression during retinal development
BMC Developmental Biology volume 10, Article number: 1 (2010)
MicroRNAs (miRNAs) are small RNA molecules (~22 nucleotides) which have been shown to play an important role both in development and in maintenance of adult tissue. Conditional inactivation of miRNAs in the eye causes loss of visual function and progressive retinal degeneration. In addition to inhibiting translation, miRNAs can mediate degradation of targeted mRNAs. We have previously shown that candidate miRNAs affecting transcript levels in a tissue can be deduced from mRNA microarray expression profiles. The purpose of this study was to predict miRNAs which affect mRNA levels in developing and adult retinal tissue and to confirm their expression.
Microarray expression data from ciliary epithelial retinal stem cells (CE-RSCs), developing and adult mouse retina were generated or downloaded from public repositories. Analysis of gene expression profiles detected the effects of multiple miRNAs in CE-RSCs and retina. The expression of 20 selected miRNAs was confirmed by RT-PCR and the cellular distribution of representative candidates analyzed by in situ hybridization. The expression levels of miRNAs correlated with the significance of their predicted effects upon mRNA expression. Highly expressed miRNAs included miR-124, miR-125a, miR-125b, miR-204 and miR-9. Over-expression of three miRNAs with significant predicted effects upon global mRNA levels resulted in a decrease in mRNA expression of five out of six individual predicted target genes assayed.
This study has detected the effect of miRNAs upon mRNA expression in immature and adult retinal tissue and cells. The validity of these observations is supported by the experimental confirmation of candidate miRNA expression and the regulation of predicted target genes following miRNA over-expression. Identified miRNAs are likely to be important in retinal development and function. Misregulation of these miRNAs might contribute to retinal degeneration and disease. Conversely, manipulation of their expression could potentially be used as a therapeutic tool in the future.
miRNAs are small RNA molecules of 18-22 nucleotides in length which regulate gene expression. Precursor RNAs transcribed either as independent genes or within introns of coding genes are processed by the RNAseIII enzymes Drosha and Dicer to form the mature miRNAs . The first miRNAs were discovered as a result of their role in the timing of adult cell fate determination in Caenorhabditis elegans. Both let-7 and lin-4 are temporally regulated and cause a decrease in expression of target proteins which control developmental stage-specific events [2–4]. The significance of miRNAs in developmental processes has now been shown in many plants and animals .
There is growing evidence that miRNAs play an important role in regulating eye development and function. In the mouse, prevention of miRNA maturation by conditional knockout in retinal tissue of the Dicer protein resulted in formation of photoreceptor rosettes and retinal degeneration and had a profound effect upon visual function, as measured by electroretinogram (ERG) . Li and Carthew reported reciprocal negative regulation between Yan protein and miR-7 in Drosophila retinal cells: Yan directly represses the expression of miR-7 in undifferentiated progenitor cells and miR-7 directly represses the expression of Yan in differentiated photoreceptor cells . The switch between these mutually exclusive expression states occurs when Yan degradation is triggered by epidermal growth factor receptor (EGFR) signaling.
Expression of various miRNAs has been demonstrated in the retina [8–10], notably a sensory organ-specific polycistronic cluster comprising miR-96, miR-182, and miR-183 . The profile of miRNAs expressed in the retina has been shown to change during degeneration in several mouse models [12, 13]. Many genes with retinal-specific functions are predicted to be targeted by miRNAs and functional interactions have been demonstrated .
miRNAs were initially thought to act through inhibition of protein translation . Whilst this mechanism has been widely confirmed  and the global impact of miRNAs on protein levels demonstrated [16, 17], it has become apparent that miRNAs also affect mRNA levels. This was first clearly demonstrated by the observation that those mRNAs with reduced expression following over-expression of a specific miRNA were enriched for predicted target sites of that miRNA . Many miRNAs display tissue-specific expression  and inhibit expression of genes which are found at lower levels relative to other tissues [18, 20]. The action of miRNAs upon mRNA levels has been widely confirmed in subsequent studies [21, 22], which have also improved understanding of the properties of functional predicted target sites [16, 23].
Inhibition of gene expression by miRNAs is mediated by base pairing between the miRNA and complementary sequences in the 3'UTR of target genes. However, unlike siRNAs, which direct cleavage of perfectly matched sequences, miRNAs have imperfect complementarity with their target sites. The features required for a functional target site are not fully understood, but include a stretch of perfect complementarity with a ~7 nucleotide 'seed' region at the 5' end of the miRNA . A range of target prediction algorithms have been developed and each miRNA is thought to target approximately 200 transcripts [25, 26]. Over 550 human miRNAs have been defined to date (miRBase v12 [27, 28]) and have been suggested to regulate 30% of genes .
We have shown recently that the effects of miRNAs upon mRNA levels are such that the consequences of expression of a specific miRNA can be detected within a single gene expression profile . Here this approach is applied to three experimental samples chosen to indicate candidate miRNAs affecting mRNA levels in ciliary epithelium-derived retinal stem cells from the adult eye (CE-RSCs), the developing retina during rod photoreceptor differentiation (postnatal murine retina) and in the mature retina (adult murine retina). This analysis may be indicative of those miRNAs with greatest effects upon target gene expression, or may reveal a subset with targets which are more susceptible to mRNA degradation. In either case, demonstration of the expression of these miRNAs confirms them as candidates for future manipulation aimed at directing differentiation or maintaining the mature retina.
In silico prediction of candidate miRNAs affecting mRNA levels
Each miRNA was mapped to its predicted target genes in each mRNA expression profile and the average levels of the transcripts of their target genes were calculated. The ranked average target gene expression values in the P4 retina are shown in Figure 1A. The probability that each miRNA has a significant effect upon gene expression was calculated using the Wilcoxon ranked sum test; the results for P4 retina are illustrated in Figure 1B. The predictions from all gene expression datasets are listed in Additional file 1: Table S1 and the significant combined probabilities are shown in Table 1 (p < 0.05). Multiple mRNA expression profiles were analysed for murine CE-RSCs, P4 and adult retina. Significant predicted effects upon mRNA expression were detected for multiple miRNAs in the P4 and adult retina and for a smaller number in CE-RSCs. At P4 many miRNAs had highly significant effects, with miR-124, miR-125 and miR-9 being particularly significant (Figure 1, Table 1 and Additional file 1: Table S1). In the adult retina miR-124 and miR-125 were again prominent, but others, including miR-24, miR-326, miR-370, miR-96 and let-7 also had highly significant predicted effects. The greater number of genes assayed for the adult retina may have contributed to the detection of several more miRNAs with significant predicted effects and the smaller p-values than at P4. Whilst most of the candidate miRNAs predicted to affect mRNA levels (P < 0.05) at P4 were also in the adult, only miR-125, miR-378 and miR-24 were detected in all the different tissue types and developmental stages. Gene expression data from human CE-RSCs from adult eyes were also analysed and although many different miRNAs were predicted to be affecting mRNA levels, it was notable that the two miRNAs with the most significant effects in murine CE-RSCs, miR-378 and miR-485-3p, were also the 3rd and 4th most significant in human CE-RSCs (Additional file 1: Table S1).
We hypothesized that miRNAs with detectable effects upon mRNA expression are more highly expressed than those with no effect. We tested this hypothesis in the adult mouse retina using miRNA expression data published by Loscher et al. . A list was compiled of miRNAs whose predicted target genes exhibited significantly low expression. Expression of miRNAs within this group was significantly higher than that of those miRNAs for which no effect on target mRNAs was predicted (Figure 2).
Experimental analysis of miRNA expression
The following miRNAs had highly significant predicted effects on target mRNA levels and were selected for analysis by RT-PCR: miR-124, miR-125, miR-9, and miR-24. A range of other candidates with p-values < 0.05 in at least one experimental condition (miR-128, miR-150, miR-204, miR-25, miR-27, miR-326 miR-34, miR-370, miR-378 and miR-485-5p) were selected to represent different predicted patterns of activity or for their lack of previous association with neural tissue. In addition, miR-122, which had p < 0.1 in CE-RSCs and P4 and adult retina, was of interest because it has been previously characterised as a highly expressed, liver-specific miRNA [19, 22, 30]; the possibility of an alternative role in the retina is intriguing. Let7d has previously been shown to be highly expressed in the retina  and was used as a positive control.
The pattern of relative expression in the three samples considered (porcine CE-RSCs, mouse P4 and adult retina) varied widely amongst miRNAs (Figure 3A). The expression of each individual miRNA corresponded broadly with its predicted effects upon target gene expression (e.g. miR-25 highest at P4 and miR-124 absent in CE-RSCs but present in P4 and adult mouse retina). However, whilst this analysis indicates the sample in which the relative expression of a single miRNA is highest, it gives no indication of the inter-miRNA variation in expression in that sample. It is also important to consider the relative absolute expression of different miRNAs within a single sample when considering the impact of each upon target mRNA expression. Therefore, starting template copy numbers were estimated for each miRNA based on threshold cycle (Ct) and amplification efficiency or by the linear regression efficiency method , which were broadly in agreement (Additional file 2: Figure S1). The miRNAs measured in each sample are shown in Figure 3B ranked by their starting template copy number in the RT-PCR assays. Those miRNAs with the most significant predicted effects upon mRNA expression are highly expressed. Notably, miR-125 has amongst the most significant predicted effects upon mRNA expression at all stages (Table 1) and miR-125a and miR-125b are the most highly expressed miRNAs in all samples (Figure 3B). Furthermore, miR-124 is highly expressed in P4 and adult retina in accordance with its predicted effects. To test the expected relationship between predicted effects and miRNA expression, the estimated miRNA copy numbers were plotted against predicted probability (Figure 4). This revealed a significant correlation (Spearman rank correlation, p = 4.0E-04), with the majority of most highly expressed miRNAs having highly significant predicted effects.
In situ hybridization (ISH) was performed to localise expression of selected miRNAs (Figure 5). Expression in the adult murine retina of the well-characterised neural miRNA miR-124 concurred with published reports . In P4 retina, miR-124 positive signal was mainly colocalised with ganglion and amacrine cells. Expression of miR-34a, which has not previously been described in the retina, was detected in the porcine CE-RSC neurospheres and throughout the mouse P4 retina. In the adult mouse retina the strongest signal was observed in the inner nuclear layer (INL) and in some cells in the ganglion cell layer (GCL). miR-128 was expressed in all samples, notably in CE-RSC neurospheres and the inner plexiform layer of the retina at P4. miR-125b was also expressed at all stages with strong signal in the inner plexiform layer at P4. In the adult retina miR-125b was expressed in the inner and outer nuclear layers. No expression was observed for miR-122 and miR-378 (data not shown).
Predicted targets of candidate miRNAs affecting mRNA levels
In order to identify the genes most likely regulated by miRNAs in CE-RSCs, P4 and adult retina, the target predictions for all the miRNAs shown to have a significant effect upon gene expression were combined. This revealed that some genes were predicted to be targeted by almost half of the miRNAs identified in each sample (Table 2 and Additional file 3: Table S2). These genes are the most likely to be real miRNA targets because if some predicted sites are not functional others may be active. Also, the co-operative action of multiple miRNAs may result in greater effects . Six genes targeted by three miRNAs with the highest expression and greatest predicted effects (miR-124; miR-125 and miR-9.) were selected for validation: ACCN2; ETS1; KLF13; LIN28B; NFIB and SH2B3. Following transfection of HEK293 cells with a pool of miR-124, miR-125 and miR-9 miRNA mimics the mRNA expression of 5 of these 6 genes was significantly reduced (Figure 6).
Comparison of the predicted miRNA target genes with those known to be involved in retinal function and disease  highlighted potential genes of interest (Table 3 and Additional file 3: Table S2). There was considerable overlap between known 'retinal genes' targeted by multiple miRNAs in P4 and adult retina. For example, a gene involved in the elongation of very long chain fatty acids (ELOVL4) is a predicted target of miRNAs with the most significant predicted effects on gene expression and the highest expression in RT-PCR assays. ELOVL4 has been implicated in Stargardts disease and macular dystrophy  and is expressed in photoreceptor inner segments , in agreement with the ISH localisation of miR-124 (Figure 5), its predicted regulator. Likewise, Neurocalcin delta (NCALD) is expressed in retinal amacrine and ganglion cells  in agreement with the expression of miR-204, one of its potential regulators .
To establish the cellular functions which are regulated by the miRNAs shown to affect mRNA expression, enrichment of functional annotation terms within the predicted target genes was assessed . The genes with target sites for those miRNAs identified by our analyses were selected from the set of all genes with a conserved miRNA target site predicted by TargetScan; this list was therefore used as the background to assess enrichment (various categories, including many regulatory functions, are already enriched in this group in comparison to the whole genome; data not shown). In all the different tissue types at different developmental stages, genes with two or more predicted sites for candidate miRNAs affecting mRNA levels were enriched for annotations predominately relating to regulatory functions, notably gene transcription (Additional file 4: Table S3).
miRNAs regulate mRNA expression in the retina
Gene expression is regulated by many factors at both transcriptional and posttranscriptional levels. The effect of miRNAs within these influences is detectable through analysis of microarray expression profiles; groups of mRNAs defined by presence of a putative target site for specific miRNAs were expressed at a significantly lower level than the group of non-target genes. These effects were observed in developing and adult retina and to a lesser extent in CE-RSCs. Not surprisingly, these effects were related to the expression level of the miRNA; those with extremely significant effects, such as miR-125, miR-124 and miR-9 were amongst the most highly expressed in the P4 and adult murine retina. The miRNA expression patterns determined in this study were consistent with previous reports. For example, miR-125, miR-124 and miR-9 have all been independently reported to be highly expressed in the retina [9–11]. The association between the expression level of a miRNA and its predicted effects on mRNA expression is robust because it was demonstrated for both our qRT-PCR analysis of a relatively small number of miRNAs and for global miRNA expression profiles determined by other investigators .
Genes predicted to be targeted by multiple candidate miRNAs affecting mRNA levels are the best candidates for in vivo regulation by this mechanism. It is difficult to test this assertion in retinal cells, however the predicted effects on mRNA expression were observed for five of six genes in an in vitro model. This provides some support to the predictions of individual genes affected by miRNAs, at least for those targeted by multiple miRNAs with highly significant overall effects on mRNA expression,
The pattern of miRNA activity suggests a role in retinal differentiation
The predicted effects of miRNAs were less significant in mouse CE-RSCs than in the P4 and adult retina. miRNA expression was correspondingly lower, with only 4 of 20 miRNAs displaying highest expression in the CE-RSCs (although this may have been partly due to the porcine origin of the CE-RSCs used for RT-PCR analyses). The absolute levels of these miRNAs were also lower than those highly expressed in P4 and adult retina. This is consistent with reports that miRNA expression increases from very low levels at embryonic day 10, when the retina is comprised of retinal progenitor cells with similar properties and gene expression patterns to CE-RSCs , to higher levels in the P1 and adult retina . This pattern of miRNA expression and predicted effects supports a role for miRNAs in retinal differentiation, maturation and maintenance of the adult retina. A similar observation was noted by Wienholds and co-authors, when they were able to detect only a few miRNAs during early zebrafish development .
Although generally lower than in retina, miRNA expression was described for the first time in this study in CE-RSCs. Human and porcine eyes are similar in terms of morphology and physiology. We have previously characterised porcine CE-RSCs  and therefore, although the availability of mouse expression data for CE-RSCs, P4 and adult retina made it the choice for in silico analyses, porcine CE-RSCs were used to confirm expression of predicted miRNAs. CE-RSCs from various species have been shown to have very similar properties and gene expression patterns [38–41] Of the 4 miRNAs expressed most highly in CE-RSCs according to RT-PCR, miR-24 was also predicted to have a significant effect (p < 0.05) upon mRNA expression and miR-122 had p = 0.08. It has recently been shown that miR-24a represses apoptosis and is required for proper retinal development . These miRNAs could also play a role in maintaining the progenitor cell state, as has been shown in other tissues; over-expression of miR-24 causes a delay in maturation of hematopoietic progenitor cells  and over-expression of miR-122 delays differentiation of human embryonic stem cells .
Many of the specific miRNAs shown to affect mRNA expression in the P4 and adult retina have previously been implicated in neuronal differentiation. For example, miR-124 has been associated with the transition from neural progenitor to differentiated neuron in the zebrafish brain . In mouse brain miR-124 expression is restricted to mature neurons  and it promotes neuronal differentiation by triggering brain-specific alternative pre-mRNA splicing . The reported absence of miR-124 in neural and retinal stem/progenitor cells in vivo [9, 48, 49] concurs with the lack of detectable expression in the CE-RSC spheres. Expression of miR-124 in the immature retina (P4) was lower than that in the adult retina and the ISH signal at P4 mainly corresponded to early differentiating neurons, retinal ganglion and amacrine cells. Our ISH for miR-124 in the adult retina confirmed previous reports  and showed miR-124 expression in all retinal layers, with particularly high levels in photoreceptor inner segments.
miR-9 promotes progression of neurogenesis in the zebrafish brain . Expression of miR-9 was not detected in CE-RSCs, but in contrast to miR-124 it peaked in the P4 retina. This pattern is in agreement with Xu et al , who reported a very low level of miR-9 expression at E10 which gradually increased and peaked at P10, being lower in the adult retina. Notably, the transcription factor Hes1, which maintains retinal progenitor pools during development is a predicted target for miR-9 . Similarly to retinal progenitor cells, high levels of Hes1 mRNA are found in CE-RSC spheres (Yanagi et al 2006, our unpublished data). Decrease of miR-9 expression past P10 correlates with the completion of retinal cell differentiation. miR-125 is another miRNA with highly significant predicted effects upon mRNA expression and the high expression levels for miR-125a and miR-125b detected by RT-PCR are in agreement with previous miRNA array data [8, 11]. Together with miR-124 and miR-9, miR-125b is also induced during neural differentiation of embryonic stem cells .
miRNAs differ in their effects upon mRNA turnover and translational repression
Our results suggest that, although many miRNAs are present in the retina, only a few affect mRNA levels sufficiently to have a detectable influence on target gene expression. Despite conclusive evidence of the high expression of miRNAs such as miR-181, -182 and -183 in the retina [8, 9] the genes targeted by these miRNAs did not show significant down regulation at the mRNA level in our analysis. Conversely miR-122, which has previously been shown to direct mRNA cleavage , was expressed at a low level, but had a significant predicted effect upon target gene expression in several samples (combined p = 0.08). Interestingly, miR-125b, one of the miRNAs with the strongest predicted effect upon mRNA levels, has been shown to accelerate deadenylation leading to rapid mRNA decay . It is now accepted that miRNAs can mediate both translational repression and accelerated mRNA turnover . One explanation for the detection of the effect of some highly expressed retinal miRNAs upon mRNA levels and not others could be that they are acting preferentially at the level of mRNA turnover.
The evidence presented in this study suggests that many miRNAs affect the expression levels of their target mRNAs during retinal development. The enrichment of the predicted target genes of these miRNAs for regulatory functions is in keeping with the proposed regulatory role of miRNAs. The lack of a mRNA 'signature' for other miRNAs which are also known to be highly expressed in the retina, suggests that certain miRNAs are less active at the mRNA level and perhaps act predominately upon translation. The identity and expression pattern of those miRNAs which were detected by analysis of target gene expression, such as miR-124, miR-125 and miR-9, provides further evidence that miRNAs play a central role in neuronal differentiation during retinal development.
All animal procedures were performed in compliance with the UK Animals (Scientific Procedures) Act 1986. Adult mice were euthanised in a CO2 chamber and P4 neonatal mice were lethally anesthetised by intra-peritoneal (IP) injection of an overdose of pentobarbitone. Eyes were enucleated and either processed for ISH or the retinas dissected and snap frozen in liquid nitrogen for RNA isolation.
Isolation and culture of CE-RSCs
Mixed sex White Landrace pigs were obtained from the Department of Agriculture and Rural Development Northern Ireland, Hillsborough, UK. CE-RSCs were isolated and grown as described previously . Briefly, 1-2 week old pigs were lethally anesthetised, eyes were enucleated and placed into oxygenated artificial cerebral spinal fluid [aCSF: 124 mM NaCl,5 mM KCl,1.3 mM MgCl2, 26 mM NaHCO3, and 10 mM D-glucose,(Sigma-Aldrich, Poole, UK)]. A strip of ciliary body was dissected and enzymatically digested in Hanks' Balanced Salt Solution (HBSS) containing 2 mg/ml dispase (Sigma-Aldrich, Poole, UK) for 20 minutes at 37°C, followed by digestion in EBSS containing 1.33 mg/ml trypsin, 0.67 mg/ml hyaluronidase and 78 units/ml collagenase (Sigma-Aldrich, Poole, UK) for 20 min at 37°C. The supernatant was decanted and replaced with serum-free medium (SFM, DMEM/F12 (1:1) containing 0.6% (w/v) glucose, 2 mM glutamine, 5 mM HEPES buffer, 2% (v/v) B27, 100 units/ml penicillin and 100 units/ml streptomycin) with 1 mg/ml trypsin inhibitor (Invitrogen, Paisley UK). The cellular debris was gently triturated and dissociated into single cells with a fire-polished pipette. Cells were pelleted at 800 rpm for 10 minutes, resuspended in SFM and passed through a 40 μm cell strainer (BD Biosciences, USA). Cells were counted and plated at a density of 3 × 104 cells/ml in SFM supplemented with 20 ng/ml of epidermal growth factor (EGF) and 10 ng/ml basic fibroblast growth factor (bFGF) (Invitrogen, Paisley, UK). After 7 days newly formed sphere colonies were collected, pelleted at 800 rpm for 10 minutes and digested in Accumax cell counting solution (ICT, San Diego, USA) for 20 min at RT and mechanically dissociated into single cells by pipetting. Cells were washed once in SFM and plated at a density of 3 × 104 cells/ml. After 7 days secondary spheres were collected.
miRNA target confirmation
HEK293-T cells were transfected at 80% confluency with miRNA mimics (Qiagen, Crawley, UK) using Turbofect transfection reagent (Fermentas, York, UK). Cells were collected for RNA isolation at 48 hours post transfection.
For mRNA analyses RNA was extracted using an RNeasy Micro Kit with on column DNAse digestion (Qiagen, West Sussex, UK) according to the manufacturer's protocol. For miRNA analyses total RNA was extracted using a mirVana extraction kit (Ambion, Foster City, CA) according to the manufacturer's protocol.
mRNA microarray analysis and gene expression datasets
RNA samples from retinas of three P4 mice were labeled and hybridized to GeneChip® 430A 2.0 mouse microarrays (Affymetrix, Santa Clara, Ca.) according to the manufacturer's protocols. Raw data were analysed using the Affymetrix MAS 5.0 algorithm and the results deposited in NCBI's Gene Expression Omnibus (GEO) [54, 55]. Publicly available mRNA expression data were downloaded from GEO. All the data were generated from Affymetrix genechips. Samples analysed on the MOE430A platform were available for CE-RSC and are directly comparable with the Mouse 430A 2.0 platform used for the P4 retina samples (~14,000 well characterised transcripts). Data from the MOE430B platform (~22,000 transcripts) was also available for several CE-RSC samples. The Mouse430_2 platform used for adult retinal expression interrogated all the transcripts in the A and B platforms. The tissues and GEO sample accession numbers are listed in Table 4.
miRNA target predictions
Predictions from TargetScan [23, 24] were used because we have previously demonstrated the ability of this algorithm to predict known tissue-specific miRNAs . TargetScan's requirement for a perfect match to the seed region and cross-species conservation reduces the false-positive rate [24, 25, 56]. By minimising background, use of this algorithm maximises the ability to detect effects on expression of real miRNA target genes. Complete miRNA datasets published by Lewis et al.  were downloaded from the TargetScan website (version 3.1) . For controls, 5 random datasets were generated for each of the gene expression datasets. For each of the random sets, the same numbers of target genes as predicted by TargetScan were randomly assigned to each of the miRNAs. The relationships between candidate miRNAs predicted to affect mRNA levels in the retina and their target genes were analysed using Microsoft Access and 'targets of conserved families' from Targetscan version 4.2. Functional characterisation of predicted target genes was performed using the DAVID bioinformatics database .
Bioinformatic analysis of target mRNA levels
The relative expression of the predicted target gene sets for each miRNA were compared as described previously using stand alone Java and R  programs in conjunction with MS Access and MS Excel . Only single gene-specific Affymetrix probesets (suffix_at) were considered and only those expressed in each dataset analysed (designated by a 'present' call if available, or alternatively by a signal strength greater than the median value). Probeset IDs were converted to their cognate Gene Symbols, using the Biomart online suite . For each mRNA expression dataset (Table 4), lists of expressed genes targeted by each miRNA family and their respective expression levels were compiled. The mean expression level was calculated for genes represented by > 1 probeset to provide the 'Average Target Gene Signal'. Only miRNAs with > 50 predicted targets were considered. The expression values of the predicted targets of a specific miRNA in each individual dataset were compared with the set of expression values in that dataset of all the genes predicted to be targeted by miRNAs. To determine whether there was a significant difference between the medians of the ranked gene expression values of the two sets, the nonparametric one-sided Wilcoxon rank sum test was employed. Fisher's combined probability test was used to calculate a composite test statistic from the individual p-values associated with the multiple expression datasets from CE-RSCs, P4 or adult retina. For every miRNA with a combined p < 0.05 we considered that the expression of the set of genes targeted by that miRNA in that tissue was significantly less than the average expression of all the targeted genes.
RT for detection of mRNAs was performed with 1 μg of RNA using random hexamers and SuperscriptIII (Invitrogen, Paisley, UK). Specific miRNAs were detected using a modified version of the method described by Shi and Chiang , in which mature miRNAs are polyadenylated and target sequences for a reverse primer are subsequently incorporated into cDNA by use of an oligo dT adapter. One microgram of RNA was polyadenylated using poly(A) polymerase (PAP; Ambion, Foster City, CA) at 37°C for 1 hour in a 25-μl reaction mixture. RNAs were then reverse transcribed with 200 U reverse transcriptase (SuperScript III; Invitrogen, Paisley, UK) and 0.5 μg poly (T) adapter (3' rapid amplification of complementary DNA ends (RACE) adapter in the FirstChoice RLM-RACE kit; Ambion). Primers for specific miRNAs were based on miRNA sequences obtained from miRBase [62, 63] (Table 5). The reverse primer was the 3' adapter primer (3' RACE outer primer in the FirstChoice RLM-RACE kit; Ambion).
PCR was performed for 45 cycles with denaturation at 94°C for 30 seconds, annealing at 55°C for 30 seconds, and extension at 72°C for 30 seconds (LightCycler 480: Roche, Mannheim, Germany). PCR products were analyzed by polyacrylamide gel electrophoresis (20%; Invitrogen) to confirm the predicted size (approximately 60 bp, including mature miRNA and adapter sequences).
Relative RNA expression was determined by measuring the concentration of template in each sample at the threshold cycle (Ct) from a standard curve (Log concentration against Ct) generated from dilutions of a pooled cDNA sample. Relative expression determined by an alternative method employing delta Ct values and efficiency values calculated from a standard curve gave very similar results (Additional file 2: Figure S1). Starting template copy number was estimated from Ct values using the equation: R0 = RCt /(1 + E)Ct where Ct is the threshold cycle, RCt is the fluorescence at this cycle, E is amplification efficiency (calculated from a standard curve) and R0 is starting fluorescence, which is proportional to the starting template quantity. Starting template copy numbers estimated using the linear regression efficiency method described by Rutledge et al  were in broad agreement with the above method (Additional file 2: Figure S1) and mean values from the two methods were calculated.
In situ hybridization (ISH)
Mice were euthanised as described above. Eyes were enucleated, fixed in 4% PFA for 1 h at 4°C and cryoprotected trough series of 10, 20 and 30% sucrose, for 1-3 h, embedded in optimal cutting temperature compound (OCT, Sakura, Japan) and snap frozen in an isopentane bath on dry ice. CE-RSC spheres were collected and fixed in 4% PFA for 30 min at 4°C, cryoprotected in 30% sucrose over night, embedded in OCT and snap frozen on dry ice. Locked nucleic acid probes (miRCURY™ LNA Detection, Exiqon) for the candidate miRNAs and a scrambled miRNA control probe (Table 5), were pre-labelled with digoxygenin (DIG) at the 5' end. 10 μm cryosections were fixed in 4% PFA for 10 minutes, rinsed in 1× PBS and 100 μl of hybridisation mix (50% formamide, 0.3 M NaCl, 20 mM Tris-HCl pH 8, 10% dextran sulphate, 1× Denhardt's solution, 1 mg/ml yeast rRNA) containing 2 μl of probe/ml was applied to sections. Slides were covered with glass cover slips and hybridisation was performed overnight at 20°C below the Tm. Slides were washed twice in 1 × SSC, 50% formamide, 0.1% Tween 20 for 30 min, and once in 1 × SSC, 0.1% Tween 20 and 0.2 × SSC 0.1% Tween 20 for 30 minutes at 65°C. Two washes in 1× MABT (100 mM maleic acid, 150 mM NaCl, 0.1% Tween 20, pH 7.5) were performed at RT for 30 minutes. After incubation in blocking solution containing 1× MABT, 2% blocking reagent (Roche, Mannheim, Germany) and 20% heat inactivated sheep serum (Sigma-Aldrich, Poole, UK) for 2 hours at RT, anti-digoxygenin Fab fragments (Roche, Mannheim, Germany) were applied at 1:5,000 dilution overnight at RT. Detection staining was performed by incubation with the alkaline phosphatase chromogenic substrate 5-bromo-4-chloro-3-indolylphosphate/nitro-blue tetrazolium (BCIP/NBT) (Roche, Mannheim, Germany).
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854. 10.1016/0092-8674(93)90529-Y.
Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000, 403 (6772): 901-906. 10.1038/35002607.
Wightman B, Ha I, Ruvkun G: Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell. 1993, 75 (5): 855-862. 10.1016/0092-8674(93)90530-4.
Stefani G, Slack FJ: Small non-coding RNAs in animal development. Nat Rev Mol Cell Biol. 2008, 9 (3): 219-230. 10.1038/nrm2347.
Damiani D, Alexander JJ, O'Rourke JR, McManus M, Jadhav AP, Cepko CL, Hauswirth WW, Harfe BD, Strettoi E: Dicer inactivation leads to progressive functional and structural degeneration of the mouse retina. J Neurosci. 2008, 28 (19): 4878-4887. 10.1523/JNEUROSCI.0828-08.2008.
Li X, Carthew RW: A microRNA mediates EGF receptor signaling and promotes photoreceptor differentiation in the Drosophila eye. Cell. 2005, 123 (7): 1267-1277. 10.1016/j.cell.2005.10.040.
Ryan DG, Oliveira-Fernandes M, Lavker RM: MicroRNAs of the mammalian eye display distinct and overlapping tissue specificity. Mol Vis. 2006, 12: 1175-1184.
Karali M, Peluso I, Marigo V, Banfi S: Identification and characterization of microRNAs expressed in the mouse eye. Invest Ophthalmol Vis Sci. 2007, 48 (2): 509-515. 10.1167/iovs.06-0866.
Makarev E, Spence JR, Del Rio-Tsonis K, Tsonis PA: Identification of microRNAs and other small RNAs from the adult newt eye. Mol Vis. 2006, 12: 1386-1391.
Xu S, Witmer PD, Lumayag S, Kovacs B, Valle D: MicroRNA (miRNA) Transcriptome of Mouse Retina and Identification of a Sensory Organ-specific miRNA Cluster. J Biol Chem. 2007, 282 (34): 25053-25066. 10.1074/jbc.M700501200.
Loscher CJ, Hokamp K, Wilson JH, Li T, Humphries P, Farrar GJ, Palfi A: A common microRNA signature in mouse models of retinal degeneration. Exp Eye Res. 2008, 87 (6): 529-534. 10.1016/j.exer.2008.08.016.
Loscher CJ, Hokamp K, Kenna PF, Ivens AC, Humphries P, Palfi A, Farrar GJ: Altered retinal microRNA expression profile in a mouse model of retinitis pigmentosa. Genome Biol. 2007, 8 (11): R248-10.1186/gb-2007-8-11-r248.
Arora A, McKay GJ, Simpson DA: Prediction and Verification of miRNA Expression in Human and Rat Retinas. Invest Ophthalmol Vis Sci. 2007, 48 (9): 3962-3967. 10.1167/iovs.06-1221.
Pillai RS, Bhattacharyya SN, Filipowicz W: Repression of protein synthesis by miRNAs: how many mechanisms?. Trends Cell Biol. 2007, 17 (3): 118-26. 10.1016/j.tcb.2006.12.007.
Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455 (7209): 64-71. 10.1038/nature07242.
Selbach M, Schwanhausser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 455 (7209): 58-63. 10.1038/nature07228.
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. 10.1038/nature03315.
Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002, 12 (9): 735-739. 10.1016/S0960-9822(02)00809-6.
Sood P, Krek A, Zavolan M, Macino G, Rajewsky N: Cell-type-specific signatures of microRNAs on target mRNA expression. Proc Natl Acad Sci USA. 2006, 103 (8): 2746-2751. 10.1073/pnas.0511045103.
Yu Z, Raabe T, Hecht NB: MicroRNA Mirn122a reduces expression of the posttranscriptionally regulated germ cell transition protein 2(Tnp2) messenger RNA (mRNA) by mRNA cleavage. Biol Reprod. 2005, 73 (3): 427-433. 10.1095/biolreprod.105.040998.
Krutzfeldt J, Rajewsky N, Braich R, Rajeev KG, Tuschl T, Manoharan M, Stoffel M: Silencing of microRNAs in vivo with 'antagomirs'. Nature. 2005, 438 (7068): 685-689. 10.1038/nature04303.
Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27 (1): 91-105. 10.1016/j.molcel.2007.06.017.
Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.
Maziere P, Enright AJ: Prediction of microRNA targets. Drug Discov Today. 2007, 12 (11-12): 452-458. 10.1016/j.drudis.2007.04.002.
Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microRNA target predictions. Nat Genet. 2005, 37 (5): 495-500. 10.1038/ng1536.
Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, D154-8. 36 Database
Arora A, Simpson DA: Individual mRNA expression profiles reveal the effects of specific microRNAs. Genome Biol. 2008, 9 (5): R82-10.1186/gb-2008-9-5-r82.
Elmen J, Lindow M, Silahtaroglu A, Bak M, Christensen M, Lind-Thomsen A, Hedtjarn M, Hansen JB, Hansen HF, Straarup EM, McCullagh K, Kearney P, Kauppinen S: Antagonism of microRNA-122 in mice by systemically administered LNA-antimiR leads to up-regulation of a large set of predicted target mRNAs in the liver. Nucleic Acids Res. 2008, 36 (4): 1153-62. 10.1093/nar/gkm1113.
Rutledge RG, Stewart D: A kinetic-based sigmoidal model for the polymerase chain reaction and its application to high-capacity absolute quantitative real-time PCR. BMC Biotechnol. 2008, 8: 47-10.1186/1472-6750-8-47.
Zhang K, Kniazeva M, Han M, Li W, Yu Z, Yang Z, Li Y, Metzker ML, Allikmets R, Zack DJ, Kakuk LE, Lagali PS, Wong PW, MacDonald IM, Sieving PA, Figueroa DJ, Austin CP, Gould RJ, Ayyagari R, Petrukhin K: A 5-bp deletion in ELOVL4 is associated with two related forms of autosomal dominant macular dystrophy. Nat Genet. 2001, 27 (1): 89-93. 10.1038/83817.
Lagali PS, Liu J, Ambasudhan R, Kakuk LE, Bernstein SL, Seigel GM, Wong PW, Ayyagari R: Evolutionarily conserved ELOVL4 gene expression in the vertebrate retina. Invest Ophthalmol Vis Sci. 2003, 44 (7): 2841-2850. 10.1167/iovs.02-0991.
Krishnan A, Venkataraman V, Fik-Rymarkiewicz E, Duda T, Sharma RK: Structural, biochemical, and functional characterization of the calcium sensor neurocalcin delta in the inner retinal neurons and its linkage with the rod outer segment membrane guanylate cyclase transduction system. Biochemistry. 2004, 43 (10): 2708-2723. 10.1021/bi035631v.
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57. 10.1038/nprot.2008.211.
Das AV, James J, Rahnenfuhrer J, Thoreson WB, Bhattacharya S, Zhao X, Ahmad I: Retinal properties and potential of the adult mammalian ciliary epithelium stem cells. Vision Res. 2005, 45 (13): 1653-1666. 10.1016/j.visres.2004.12.017.
Wienholds E, Plasterk RH: MicroRNA function in animal development. FEBS Lett. 2005, 579 (26): 5911-5922. 10.1016/j.febslet.2005.07.070.
Gu P, Harwood LJ, Zhang X, Wylie M, Curry WJ, Cogliati T: Isolation of retinal progenitor and stem cells from the porcine eye. Mol Vis. 2007, 13: 1045-1057.
Coles BL, Angenieux B, Inoue T, Del Rio-Tsonis K, Spence JR, McInnes RR, Arsenijevic Y, Kooy van der D: Facile isolation and the characterization of human retinal stem cells. Proc Natl Acad Sci USA. 2004, 101 (44): 15772-15777. 10.1073/pnas.0401596101.
Inoue Y, Yanagi Y, Tamaki Y, Uchida S, Kawase Y, Araie M, Okochi H: Clonogenic analysis of ciliary epithelial derived retinal progenitor cells in rabbits. Exp Eye Res. 2005, 81 (4): 437-445. 10.1016/j.exer.2005.03.003.
Xu S, Sunderland ME, Coles BL, Kam A, Holowacz T, Ashery-Padan R, Marquardt T, McInnes RR, Kooy van der D: The proliferation and expansion of retinal stem cells require functional Pax6. Dev Biol. 2007, 304 (2): 713-721. 10.1016/j.ydbio.2007.01.021.
Walker JC, Harland RM: microRNA-24a is required to repress apoptosis in the developing neural retina. Genes Dev. 2009, 23 (9): 1046-1051. 10.1101/gad.1777709.
Wang Q, Huang Z, Xue H, Jin C, Ju XL, Han JD, Chen YG: MicroRNA miR-24 inhibits erythropoiesis by targeting activin type I receptor ALK4. Blood. 2008, 111 (2): 588-595. 10.1182/blood-2007-05-092718.
Tzur G, Levy A, Meiri E, Barad O, Spector Y, Bentwich Z, Mizrahi L, Katzenellenbogen M, Ben-Shushan E, Reubinoff BE, Galun E: MicroRNA expression patterns and function in endodermal differentiation of human embryonic stem cells. PLoS ONE. 2008, 3 (11): e3726-10.1371/journal.pone.0003726.
Kapsimali M, Kloosterman WP, de Bruijn E, Rosa F, Plasterk RH, Wilson SW: MicroRNAs show a wide diversity of expression profiles in the developing and mature central nervous system. Genome Biol. 2007, 8 (8): R173-10.1186/gb-2007-8-8-r173.
Smirnova L, Grafe A, Seiler A, Schumacher S, Nitsch R, Wulczyn FG: Regulation of miRNA expression during neural cell specification. Eur J Neurosci. 2005, 21 (6): 1469-1477. 10.1111/j.1460-9568.2005.03978.x.
Makeyev EV, Zhang J, Carrasco MA, Maniatis T: The MicroRNA miR-124 promotes neuronal differentiation by triggering brain-specific alternative pre-mRNA splicing. Mol Cell. 2007, 27 (3): 435-448. 10.1016/j.molcel.2007.07.015.
Bak M, Silahtaroglu A, Moller M, Christensen M, Rath MF, Skryabin B, Tommerup N, Kauppinen S: MicroRNA expression in the adult mouse central nervous system. RNA. 2008, 14 (3): 432-44. 10.1261/rna.783108.
Visvanathan J, Lee S, Lee B, Lee JW, Lee SK: The microRNA miR-124 antagonizes the anti-neural REST/SCP1 pathway during embryonic CNS development. Genes Dev. 2007, 21 (7): 744-749. 10.1101/gad.1519107.
Leucht C, Stigloher C, Wizenmann A, Klafke R, Folchert A, Bally-Cuif L: MicroRNA-9 directs late organizer activity of the midbrain-hindbrain boundary. Nat Neurosci. 2008, 11 (6): 641-648. 10.1038/nn.2115.
Krichevsky AM, Sonntag KC, Isacson O, Kosik KS: Specific microRNAs modulate embryonic stem cell-derived neurogenesis. Stem Cells. 2006, 24 (4): 857-864. 10.1634/stemcells.2005-0441.
Wu L, Fan J, Belasco JG: MicroRNAs direct rapid deadenylation of mRNA. Proc Natl Acad Sci USA. 2006, 103 (11): 4034-4039. 10.1073/pnas.0510928103.
Wu L, Belasco JG: Let me count the ways: mechanisms of gene regulation by miRNAs and siRNAs. Mol Cell. 2008, 29 (1): 1-7. 10.1016/j.molcel.2007.12.010.
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 2007, D760-5. 10.1093/nar/gkl887. 35 Database
Gene Expression Omnibus (GEO). [http://www.ncbi.nlm.nih.gov/projects/geo/]
Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115 (7): 787-798. 10.1016/S0092-8674(03)01018-3.
TargetScan: Prediction of microRNA targets. [http://www.targetscan.org/]
The Comprehensive R Archive Network. [http://cran.r-project.org/]
BioMart Project. [http://www.biomart.org/]
Fisher RA: Statistical Methods for Research Workers. 1932, Edinburgh: Oliver and Boyd
Shi R, Chiang VL: Facile means for quantifying microRNA expression by real-time PCR. BioTechniques. 2005, 39 (4): 519-525. 10.2144/000112010.
miRBase: the microRNA database. [http://www.mirbase.org/]
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, D140-4. 10.1093/nar/gkj112. 34 Database
Paper W, Kroeber M, Heersink S, Stephan DA, Fuchshofer R, Russell P, Tamm ER: Elevated amounts of myocilin in the aqueous humor of transgenic mice cause significant changes in ocular gene expression. Exp Eye Res. 2008, 87 (3): 257-267. 10.1016/j.exer.2008.06.006.
Steele MR, Inman DM, Calkins DJ, Horner PJ, Vetter ML: Microarray analysis of retinal gene expression in the DBA/2J model of glaucoma. Invest Ophthalmol Vis Sci. 2006, 47 (3): 977-985. 10.1167/iovs.05-0865.
Akimoto M, Cheng H, Zhu D, Brzezinski JA, Khanna R, Filippova E, Oh EC, Jing Y, Linares JL, Brooks M, Zareparsi S, Mears AJ, Hero A, Glaser T, Swaroop A: Targeting of GFP to newborn rods by Nrl promoter and temporal expression profiling of flow-sorted photoreceptors. Proc Natl Acad Sci USA. 2006, 103 (10): 3890-3895. 10.1073/pnas.0508214103.
The authors would like to acknowledge the contribution of all the researchers who have made their gene expression data publicly available and without whom this study would not have been possible.
AA performed bioinformatic analyses and some laboratory work; JG-F performed RT-PCR and ISH and contributed to manuscript preparation. LH performed miRNA transfections and mRNA RT-PCR. MD performed microarray analysis with P4 mice; TC contributed to supervision and manuscript preparation; DAS conceived and supervised the study, performed bioinformatic analyses and prepared the manuscript. All authors read and approved the final manuscript.
Amit Arora, Jasenka Guduric-Fuchs contributed equally to this work.
Electronic supplementary material
Additional file 1: Table S1. Significance of predicted effects of miRNAs upon mRNA expression in all gene expression datasets (CE-RSCs, P4 and adult retina) and p-values calculated using Fisher's combined probability test. (XLS 86 KB)
Additional file 2: . miRNA expression detected by RT-PCR. (A) Comparison of relative miRNA expression calculated using a standard curve (as shown in Figure 3A) with the values calculated using ΔCt and efficiency; for each gene the sample with highest expression was normalized to a value of 100. Each point represents the expression of one miRNA in a single sample as determined by both methods, which yielded very similar results (R2 = 0.983). (B) Starting template copy numbers estimated for each miRNA based on threshold cycle (Ct) and amplification efficiency plotted against the values calculated by a linear regression efficiency method (R2 = 0.661). (PDF 143 KB)
Additional file 3: Table S2. All genes predicted to be targeted by candidate miRNAs affecting mRNA levels in each tissue. Genes assessed experimentally for miRNA interactions are highlighted in bold. (XLS 478 KB)
Additional file 4: . The categories of functional annotations enriched within those genes predicted to be targeted by two or more candidate miRNAs affecting mRNA levels within CE-RSCs, PE or adult retina. Analyses were performed using the DAVID bioinformatics database ; categories with p < 0.05 after correction for multiple comparisons (Benjamini) are listed. (XLS 32 KB)
Authors’ original submitted files for images
About this article
Cite this article
Arora, A., Guduric-Fuchs, J., Harwood, L. et al. Prediction of microRNAs affecting mRNA expression during retinal development. BMC Dev Biol 10, 1 (2010). https://doi.org/10.1186/1471-213X-10-1
- miRNA Expression
- Predict Target Gene
- Retinal Development
- Retinal Progenitor Cell