A highly expressed miR-101 isomiR is a functional silencing small RNA
© Llorens et al.; licensee BioMed Central Ltd. 2013
Received: 15 August 2012
Accepted: 30 January 2013
Published: 15 February 2013
Skip to main content
© Llorens et al.; licensee BioMed Central Ltd. 2013
Received: 15 August 2012
Accepted: 30 January 2013
Published: 15 February 2013
MicroRNAs (miRNAs) are short non-coding regulatory RNAs that control gene expression usually producing translational repression and gene silencing. High-throughput sequencing technologies have revealed heterogeneity at length and sequence level for the majority of mature miRNAs (IsomiRs). Most isomiRs can be explained by variability in either Dicer1 or Drosha cleavage during miRNA biogenesis at 5’ or 3’ of the miRNA (trimming variants). Although isomiRs have been described in different tissues and organisms, their functional validation as modulators of gene expression remains elusive. Here we have characterized the expression and function of a highly abundant miR-101 5’-trimming variant (5’-isomiR-101).
The analysis of small RNA sequencing data in several human tissues and cell lines indicates that 5’-isomiR-101 is ubiquitously detected and a highly abundant, especially in the brain. 5’-isomiR-101 was found in Ago-2 immunocomplexes and complementary approaches showed that 5’-isomiR-101 interacted with different members of the silencing (RISC) complex. In addition, 5’-isomiR-101 decreased the expression of five validated miR-101 targets, suggesting that it is a functional variant. Both the binding to RISC members and the degree of silencing were less efficient for 5’-isomiR-101 compared with miR-101. For some targets, both miR-101 and 5’-isomiR-101 significantly decreased protein expression with no changes in the respective mRNA levels. Although a high number of overlapping predicted targets suggest similar targeted biological pathways, a correlation analysis of the expression profiles of miR-101 variants and predicted mRNA targets in human brains at different ages, suggest specific functions for miR-101- and 5’-isomiR-101.
These results suggest that isomiRs are functional variants and further indicate that for a given miRNA, the different isomiRs may contribute to the overall effect as quantitative and qualitative fine-tuners of gene expression.
miRNAs are small non-coding RNAs that generally act as negative post-transcriptional regulators of gene expression. It is estimated that miRNAs regulate at least 20% of human genes , constituting an important layer of regulation in gene expression networks. The influence of a particular miRNA on the transcriptome is difficult to dissect, since genes can harbor binding sites for several miRNAs, and a specific miRNA can modulate the expression of hundreds of mRNA targets . miRNAs are involved in almost every biological process examined and their deregulation has been associated with a number of pathological conditions, including cancer and neurological disorders [3–6].
Mature miRNAs function as components of the RNA-Induced Silencing Complex (RISC) , guiding the RISC to specific gene targets through base-pairing interactions between the miRNA and the mRNA. Members of the Argonaute (Ago) protein family are central to RISC function. The seed region extends from the second to the seventh or eighth position of mature miRNA and has traditionally been considered as a major determinant of the target mRNA repertoire, through perfect base pairing with the 3’-untranslated regions (3’-UTR) of mRNAs . Other structures of the 3’-end of the miRNA are also involved in duplex stability, providing a supplementary site if the seed is fully paired or a compensatory site when the 5’-end has mismatches or bugles [9, 10]. However, our knowledge of the determinants governing gene targeting is far from complete. In fact, recent findings show that targeting can occur through sites other than the 3’-UTR, and that seed region perfect base pairing is not always required [11–14]. This miRNA/mRNA binding leads either to target mRNA degradation or to translational inhibition, depending, among other factors, on the extent of base-pairing and the type of RISC/Ago proteins .
A novel degree of complexity for miRNA transcriptome has been highlighted using next generation sequencing (NGS) strategies. NGS has revealed post-transcriptional editing processes in miRNAs, mainly consisting in variations in the 3’- and 5’-terminus and to a minor extent, nucleotide substitutions along the miRNA sequence [16–19]. Variations at the miRNA ends may be generated by a shift in Drosha and Dicer cleavage sites during miRNA biogenesis (trimming variants) [19, 20]. Another type of variability is related to pri-miRNA post-transcriptional editing as a consequence of adenosine or cytidine deaminase activity (nucleotide substitution variants). Finally, nucleotide additions at the 3’-end of the miRNA are another source of miRNA variability [21–23]. The resulting sequences differ slightly from the annotated miRNA and have been termed isomiRs .
Several lines of evidence suggest that isomiRs are not experimental artifacts derived from RNA degradation during sample preparation for NGS [24, 25]. In fact the frequencies of nucleotide substitutions are very diferent from sequencing error rates; and the detection of isomiRs in many species, tissues and physiological conditions [26, 27], suggests that they may have a biological role. This notion is further supported by their differential expression across developmental stages [27, 28]. This scenario suggests that, for a given miRNA, posttranscriptional gene expression may depend on different isomiRs. However, their functionality as posttranscriptional gene expression regulators needs to be proven. Our recent studies of miRNA profiling in human brain samples have shown a multitude of miRNA variants . miR-101 is an interesting example, since it presents 5’-trimming isomiRs (5’-isomiR-101) expressed at levels similar to those of the reference sequence. Here, we characterize the ability of this highly expressed human isomiR-101 to work as a functional miRNA.
Considering all the sequencing datasets examined, the percentage of sequences mapping onto miR-101 varied between tissues and cell lines, with the brain presenting the highest levels (Figure 1C). These data suggest that high miR-101 levels may be important for adult brain physiology. In fact, miR-101 expression is regulated under a number of physiological and pathological conditions, including angiogenesis , tumors [31, 32] and neurodegeneration [16, 33, 34]. In the last condition, our earlier study identified miR-101 as an up-regulated gene in HD . However, that study was an overall analysis of all the sequences mapping onto miR-101. To discern the specific expression pattern of miR-101 5’-isomiRs, we evaluated the expression of miR-101 variants harboring each of the two main seeds in control brains and brains of patients with HD (Figure 1D). Both types of variants were up-regulated in the frontal cortex and striatum of HD patients. In addition, we performed a similar analysis in other physiological conditions of the nervous system, including in vitro differentiation of the SH-SY5Y neuronal cell line to a postmitotic dopaminergic phenotype (Figure 1E) and early postnatal development and ageing of the human brain  (Figure 1F). Interestingly, the relative amounts of both types of miR-101 variants in the brain differed according to age. These data suggest that miR-101 is modulated in the development and differentiation processes of the nervous system. Furthermore, the context-dependent variations in the relative amounts of miR-101 and 5’-isomiR-101 suggest that the mechanisms regulating their biogenesis and/or stability may differ depending on the biological process.
Functional studies aiming to elucidate the possible relevance of the more abundant miR-101 isomiRs in the nervous system were performed in the neuronal cell line SH-SY5Y. We used miRNA mimics for two abundant representative sequences containing ACAGUAC or UACAGUA seeds: The miRBase reference miR-101 UACAGUACUGUGAUAACUGAA, and one of the most abundant miR-101 5’-isomiRs presenting the same length as the reference miR-101 GUACAGUACUGUGAUAACUGA (Additional file 1: Figure S1).
We then examined whether transfected miR-101 and 5’-isomiR-101 mIRIDIAN were incorporated into the RISC, using SH-SY5Y cells stably expressing Ago2-FLAG. We performed immunoprecipitation (IP) assays using anti-FLAG antibodies for Ago2-IP or control antibodies, and subsequently RNA in the IP was isolated. We detected a significant enrichment of either miR-101 or 5’-isomiR-101 in Ago2 IP. The amount of miR-101 or 5’-isomiR-101 in cells transfected with the corresponding miRNA mimics was similar (Figure 1A). However, increased amounts of miR-101 were detected in the Ago2 IP, compared with 5’-isomiR-101 (Figure 2B), suggesting that miR-101 was more efficiently loaded into RISC.
The ability of miR-101 and 5’-isomiR-101 to bind to the RISC was further confirmed in complementary experiments designed to evaluate the presence of RISC members in transfected biotinylated miR-101 and 5’-isomiR-101 (Figure 2C). mIRIDIAN 3’-biotin-miR-101 and 3’-biotin-5’-isomiR-101 were transfected in SH-SY5Y cells and then pulled down using streptavidin beads. The samples were blotted against Ago2 and Rck/p54 which represses translation and associates with the RISC complex in order to regulate miRNA mediated gene silencing . Both biotin-miR-101 and biotin-5’-isomiR-101 were able to interact with Ago2 and Rck/p54. Although similar levels of each biotinylated mIRIDIAN were detected in transfected SH-SY5Y cells, increased levels of both proteins were found in biotin-miR-101 pull downs compared with those of biotin-5’-isomiR-101. This result suggests that miR-101 binds with RISC members more efficiently, which is in agreement with results shown in Figure 1B.
To validate the presence of both variants in endogenous RISC, we immunoprecipitated Ago2 from the human frontal cortex. The control experiment included a sample in which the Ago2 antibody had previously been denaturized (d-Ago2) and was therefore unable to interact with Ago2. Total RNA was isolated from Ago2 and d-Ago2 IPs and subsequently we performed RT-qPCR for miR-101, 5’-isomiR-101, miR-29a, SNORD44 and U6. Data were expressed as the fold change obtained from the Ago2 vs d-Ago2 crossing thresholds (CTs) (Figure 2D). miR-101 and 5’-isomiR-101 were enriched in Ago2 IPs compared to control d-Ago2 antibody. miR-29a, a highly abundant miRNA in brain  was also enriched in Ago2 complexes. In contrast, the small nucleolar RNAs U6 and SNORD44, which were not expected to bind Ago2, were detected at insignificant levels. The presence of both types of miR-101 variants in Ago2 immunocomplexes was confirmed in publicly available high-throughput sequencing data of endogenous Ago2-associated sRNAs in HEK293T and mouse NIH-3T3 cells (GEO datasets GSM337571 and GSM849857 respectively).
To evaluate the possibility that 5’-isomiR-101 and miR-101 are sorted into different Agos, we analyzed publicly available small-RNA sequencing datasets from Ago1, Ago2 and Ago3 IP and total cell extracts, in THP-1 human monocytic cells  (Additional file 3: Figure S2). We found all types of miR-101 in each Ago (the trimming, substitution and 3’-additions variants). In addition, the percentage of miR-101 sequences loaded onto each Ago was comparable. Grouping the variants according to the two main seeds we found that in total cell homogenates the ratio miR-101/5’-isomiR-101 in THP-1 cells was close to 2. However, this ratio varied slightly depending on the Ago protein, with the highest levels being found in Ago2 and the lowest in Ago3. This suggests that although in terms of quantity miR-101 exceeded 5’-IsomiR-101 in all Agos, miR-101 species loaded more onto Ago2 than 5’-isomiR-101 variants. This is in agreement with our findings for SH-SY5Y.
Since miR-101 and 5’-isomiR-101 have a different seed region, we evaluated the predicted conserved targets for each variant, using the TargetScan 5.2 release , by searching for the presence of 8mer and 7mer sites that match miR-101 and 5’-isomiR-101 seed regions. TargetScan 5.2 detected that the majority of putative mRNA targets (>70%) were coincidental, suggesting that miR-101 and 5’-isomiR-101 may target highly overlapping pathways. In an attempt to identify possible differentially deregulated genes, SH-SY5Y cells were transfected with miR-101 or 5’-isomiR-101 mimics or a scrambled mimic in three independent experiments, and 48 h later overall gene expression was evaluated using Illumina microarrays (HumanHT Expression BeadChips). Transfection of each miR-101 variant resulted in few significantly deregulated genes that showed few changes in expression (q<5%, fold change < −1,2 or > 1,2; Additional file 4: Table S2–Additional file 5: Table S3). miR-101 induced a significant downregulation of 16 genes, eight being targeted by at least one prediction algorithm. Nine out of the 16 genes were also downregulated by 5’-isomiR-101 (q<10%), while the rest did not show significant variations. Thus, 5’-isomiR-101 seems to induce an overall weaker inhibition. In agreement, transfection of 5’-isomiR-101 resulted in a slight downregulation of two genes, for which miR-101 did not induce significant expression changes. These data indicate that under our experimental conditions modulation of gene expression by miR-101 or 5’-isomiR-101 may not involve mRNA degradation as the main mechanism, and suggest further that miR-101 inhibits gene expression more efficiently than 5’-isomiR-101.
To gain insights into the possible role of 5’isomiR-101 in physiological conditions, we analyzed publicly available mRNA and small-RNA datasets of human brain samples at multiple time-points between early postnatal development and aging . We chose this paradigm for two reasons: i/ miR-101 expression levels increased through aging, suggesting a role in this process (Figure 1F), and ii/ variants carrying each seed are not equally expressed at the different time-points (compare the oldest and the middle-aged individuals, Figure 1F), suggesting a seed-specific role in the modulation of the age-related genes.
To assess the possible seed-specificity in the modulation of the age-related transcriptome, we analyzed the correlation between the expression of miR-101 variants and that of the age-related genes at different ages through life. First, we generated a list of candidate age-related genes through a multiple regression analysis of the gene expression profiles using age as the predicted variable (a total of 3801). We then identified the miR-101 or 5'-isomiR-101 predicted targets among the pool of age-related genes, using the TargetScan 5.2 algorithm (837 for miR-101 and 899 for 5’-isomiR-101, out of the 3801; the majority, 73%-78%, were common to the two miR-101 seeds). Finally, we measured the correlation between the expression levels of age-related genes and those of the most abundant reference seed sequence (miR-101, UACAGUACUGUGAUAACUGAA) or the most abundant 5’-isomiR sequence (5’-isomiR-101, GUACAGUACUGUGAUAACUGA).
The negative correlation between the expression profiles of the mRNAs and those of miR-101 or 5’-isomiR-101 was considered the readout for miR-101 or 5’-isomiR-101 effective/detectable gene targeting. The expression of a total of 64 and 153 age-related genes anti-correlated with that of the ref-miR-101 or 5’-isomiR-101 respectively (r<−0.7), with a significantly lower number of common genes between the two miR-101 seeds (16%-39%, p<0.01). These results indicate that the majority of anti-correlated targets were seed-specific, which supports the notion of a selective function of each miR-101 variant in modulating age-related genes (Additional file 6: Figure S3).
To further confirm 5’-isomiR-101 functionality we determined the effect of miR-101 and 5’-isomiR-101 mimics on the expression of COX-2, MCL-1, APP, EZH2 and MKP-1 at the protein level (Figure 3B). We found that comparable expression levels of exogenously transfected miR-101 or 5’-isomiR-101 in SH-SY5Y cells decreased the expression of all proteins. However, miR-101 produced a significantly higher down-regulation of the proteins studied compared with 5’-isomiR-101 (Figure 3B). In SH-SY5Y cells, miR-101 repressed the expression of COX-2, MCL-1, EZH2, APP and MKP-1 at the same levels as those reported in other cell systems [38–42].
The decreased silencing activity of 5’-isomiR-101 agrees with the reduced binding of this variant to Ago2 complexes. To further confirm that the target mRNAs were less represented in the RISC, we transfected SH-SY5Y cells with 3’biotinylated miR-101 and 5’-isomiR-101 mimics and subsequently pulled down the complexes with streptavidin. Both 3’biotinylated-miR-101 and 3’biotinylated-5’-isomiR-101 interacted with the mRNA of the five genes analyzed in the present study, although higher mRNA levels were found in the 3’biotinylated-miR-101 pulldowns (Figure 3C).
High-throughput sequencing has revealed the existence of miRNA sequence variants for virtually all miRNAs in many species and tissues and in both physiological and pathological conditions [16, 43]. This scenario suggests a functional role for isomiRs that may introduce a new layer of complexity in gene expression regulation. Here, we characterized the functionality of a highly abundant type of miR-101 variant that differs in the seed region. miR-101 has mainly been studied for its tumor suppressive functions, by directly targeting EZH2, Cox-2, Mcl-1 and Fos [31, 39, 41, 44–46]. miR-101 has also been shown to modulate the expression of target genes involved in inflammatory processes , autophagy  and angiogenesis . In addition, our analysis showing the regulated expression of miR-101 in brain development and ageing and in a neurodegenerative process suggests important roles for this miRNA in nervous system physiology and pathology. Reinforcing this idea, genes that are relevant in brain function including APP  and ATX1  are proven targets of miR-101.
For the majority of miRNAs, the most abundant isomiRs include 3’-trimming variants and nucleotide additions in the 3’-end , suggesting that modifications in the 3’ end are more permissive. In contrast, the diversity and abundance of isomiRs affecting the 5’-end of mature miRNAs (5’-trimming variants) is considerably reduced, which may be indicative of an increased impact of these variants in biological processes. In fact, miRNA 5’-trimming variants show a modified seed region (nucleotides 2–8 of the mature sequence) which is one of the major determinants defining target genes [50, 51]. Thus, for a given miRNA, isomiRs expressed in high amounts which differ in the seed region are expected to have some impact on gene expression regulation. Our results show that among all the sequences mapping onto miR-101, 5’-trimming variants affecting a single nucleotide position are highly abundant in human samples and show variable expression depending on the sample. This suggests that the biogenesis and/or stability of miR-101 isomiRs are dependent on the cell context. 5’-trimming isomiRs may be generated by variations of the Dicer-cutting site on the precursor miRNA. In support of this idea, both types of miR-101 variants are strongly and similarly decreased (data not shown) when analyzing publicly available sRNA sequencing datasets in MCF-7 cells knocked down for Dicer . Recent results indicate that Dicer activity can be modulated by a number of proteins, interacting either with Dicer or with miRNA precursors . Thus, whether differential modulation of Dicer activity in diverse cell types and/or physiological conditions is underlying the relative abundance of 5’-trimming miR-101 isomiRs is an open question that deserves specific research.
Another possible mechanism accounting for the differential abundance of all variants mapping onto miR-101 may be related to the stability of the mature forms. It has been shown that the post-transcriptional addition of non-genome-encoded nucleotides to the 3’ end of mature miRNA increases miRNA stability or abundance [24, 54]. However, the proportion of species with 3’-nucleotide additions mapping as miR-101 or 5’-isomiR-101 did not correlate with the differential miR-101/5’-isomiRs-101 relative abundance detected in the samples examined (data not shown). This suggests that nucleotide additions may not be a major explanation for the differential relative abundance.
Several findings indicate that miR-101 5’-trimming variants are functional silencing sRNA: i) both 5’-isomiR-101 mimic and endogenous 5’-isomiR-101 were loaded into Ago2 immunocomplexes, ii) biotin-streptavidin pull-down showed that 5’-isomiR-101 binds to both Ago2 and Rck/p54 RISC components iii) biotinylated 5’-isomiR-101 pull-down assays purified known miR-101 target mRNAs and iv) 5’-isomiR-101 mimic inhibited the expression of known miR-101 targets.
In trying to identify possible candidates that are differentially modulated by miR-101 and 5’-isomiR-101, high-throughput array-based approach failed to show significant regulation of mRNAs by either variant. It may be that miR-101 modulates gene expression mainly at the translational level, producing mild effects on RNA levels that are difficult to verify in a microarray approach. In spite of this, the arrays identified a small number of genes that were regulated by only one type of variant. In addition, the non-overlapping anti-correlation between the expression profiles of miR-101 variants and a number of age-related genes in the human brain agrees with specific miR-101- and 5’-isomiR-101 gene targeting. Although selective functions of the seed variants may add complexity to miR-101 regulatory effect, it has been proposed that isomiRs function cooperatively to target common biological pathways . In agreement with this suggestion, the high number of commonly predicted targets indicates that similar biological pathways are regulated by isomiR-101 and 5’-isomiR-101.
A number of results suggest that 5’-isomiR-101 is less efficient than miR-101 in gene silencing: i) determinations of 5’-isomiR-101 and miR-101 with highly specific RT-PCR assays in IPs of endogenous Ago2 indicate that miR-101 is more efficiently loaded into the RISC, although the analysis of SH-SY5Y small RNA sequencing data suggests that 5’-isomiRs-101 are more abundant than the reference miR-101; ii) increased amounts of miR-101 were also detected in Ago2 IPs of SH-SY5Y cells transfected with equivalent quantities of miR-101 or 5’-isomiR-101 mimics; iii) biotinylated miR-101 bound increased amounts of both Ago2 and Rck/p54 RISC members compared with biotinylated 5’-isomiR-101 and iv) 5’-isomiR-101 showed a decreased ability to inhibit several miR-101 targets. Overall, these data indicate that the silencing activity of miR-101 variants is related to their efficient loading into the RISC and further suggest that the amount of the different miR-101 variants expressed within a cell does not strictly correlate with the quantities detected in the Ago2 complexes. Since Ago2 is the only Ago with mRNA slicing activity, decreased loading of 5’-isomiR-101 onto Ago2 may explain the less efficient gene silencing of at least some target genes showing RNA destabilization. However, the reduced loading of 5’-isomiR-101 in Ago1 and Ago3 detected in monocytic human cells suggests that less efficient binding to other Agos may contribute to the decreased silencing at the protein level. It has been shown that the 5’-terminal residue in a miRNA influences its loading into the RISC, with terminal 5’-U enhancing the Ago2 complex . Thus, the 5’ terminal G in the isomiR may provide a less efficient incorporation into the RISC.
Our analysis suggests that the relative amount of the two main miR-101 seeds was not equivalent in the different Agos in monocytic human cells. It has recently been shown that a subset of miRNAs give rise to distinct isomiRs which associate mainly with distinct Ago proteins, both in plants  and in human samples [58, 59], supporting an isomiR-dependent differential output in gene silencing.
In summary, our data show that 5’-isomiR-101 is a functional silencing small RNA. We propose that in addition to specific gene targeting, 5’-isomiR-101 variants present a subtle effect through regulatory mechanisms involving decreased silencing efficiency. Changes in the relative expression of seeds ratios may provide scenarios for isomiR-selective qualitative and quantitative gene expression regulation. This may add complexity to the fine-tuned nature of miRNA regulation of gene expression.
The mechanisms by which miR-101 variants are generated and a full characterization of the repertoire of targets modulated by miR-101 and 5’-isomiR-101 may shed light on the respective biological function. Extensive study of these aspects will broaden our understanding of the overall biological impact of isomiRs in different tissues and physiological states.
Anti-Ago2 was from Millipore, anti-Rck/p54 was from MBL, anti-Mcl-1 was from Cell Signaling, anti-COX-2 was from Cayman, anti-EZH2 was obtained from Dr. Kristian Helllin (BRIC-University of Copenhagen), anti-MKP-1 was from Santa Cruz, anti-APP was from Abcam and anti-tubulin and anti-FLAG (M2) were from Sigma. miRIDIAN™ hsa-miRNA-101 mimics and siGLO Green (transfection indicator) were from Dharmacon: hsa-miRNA-101 (mature sequence: UACAGUACUGUGAUAACUGAA), hsa-5’-isomiR-101 (mature sequence: GUACAGUACUGUGAUAACUGA), 3’biotinilated hsa-miRNA-101 (mature sequence: UACAGUACUGUGAUAACUGAA-Bio) and 3’biotinilated hsa-5’-isomiR-101 (mature sequence: GUACAGUACUGUGAUAACUGA-Bio).
Brain samples of individuals without clinical neurological disease and no major histopathological lesions were obtained from the Institute of Neuropathology and the University of Barcelona Brain Bank, after the informed consent of the patients or their relatives and the approval of the local ethics committee had been given. Blood samples of healthy individuals without clinical neurological disease were obtained from the Hospital Clinic of Barcelona after written informed consent of the participants and approval of the local Ethics Committee had been given.
We used publicly available small RNA sequencing datasets from the GEO repository (GSE31069, GSM522374, GSE16579, GSE18069). In addition we sequenced small RNA in human blood samples (4 individuals), human brains (amygdala area of 4 individuals without major neurodegenerative histopathological signs) and SH-SY5Y neuroblastoma cells before and after differentiation to a postmitotic phenotype as previously described . In these samples, total RNA was isolated using PAXgene blood RNA Kit (Qiagen) or the miRNAeasy kit (Qiagen). From each sample, 1μg of total RNA was used to prepare libraries according to the “TruSeq Small RNA Sample Preparation Guide” (Part# 15004197 Rev. A). Libraries were subsequently sequenced on one Illumina Genome Analyzer II run using 36 single-read cycles. For processing of the sequencing data (recognition and removal of the adapter) and characterization of miRNA variants we used a stand-alone version of the Seqbuster tool (http://estivill_lab.crg.es/seqbuster). The parameters were set up to detect isomiRs by: a) trimming at the 3’/5’-end (up to 3 nucleotides shifted from the reference dicing site), b) addition at the 3-end (up to 3 nucleotides added) and c) nucleotide substitution (up to 1 mismatch with respect to the reference miRNA sequence).
To identify the genes differently expressed in the brain at different ages we applied the same method described by Somel et al. . Briefly, the intensity of gene expression at 12 different time-points along life was used to correlate to age. Expression was fitted to 8 non-linear functions capable of detect linear, quadratic or more complex correlation between gene expression levels and age. All genes that fitted to one of those 8 models with p-value lower than 0.01 were considered age-related genes.
We considered an age-related gene to be targeted by miR-101 or 5’-isomiR-101 sequences only if: i) there is a pair complementary between the miRNA seed and the gene according to TargetScan Custom prediction algorithm  and ii) there is a negative expression correlation between the miRNA sequence and the gene. Only those genes that have been detected in the gene expression profile of the brain samples, according to the public data from the original work, were used as input for TargetScan Custom algorithm. For each miR-101-sequence-gene pair predicted by TargetScan Custom, a linear regression model was performed by the standard R function (’lm’). To consider a gene regulated by a miR-101 variant, the p-value of the model should be less than 0,01 (F-test) and Pearson’s product moment correlation lower than −0.70. Normalized intensities of microarray data for gene expression, and logarithms of counts for miRNA sequences were used. The numbers of assigned miR-101 or 5’-isomiR-101 predicted targets were calculated according to the first requirement (seed match) or both of them (seed match and anti-correlation). The statistical differences between the numbers of commonly predicted miR-101 and 5’-isomiR-101 targets, considering or not the anti-correlation criterion was determined with the prop.test function in R.
SH-SY5Y and SH-SY5Y cells stably expressing Ago2-FLAG (SH-SY5Y-Ago2-FLAG)  were cultured at 37°C in a 95:5 air/CO2 water saturated atmosphere in 1:1 mixture of DMEM and Ham’s F12 medium and 10% supplemental fetal bovine serum. For cell transfection, SH-SY5Y cells were seeded at 0.2 × 106 and 50 ng of hsa-miRNA mimic or SiGLO Green (Control) were transfected with Lipofectamine Plus according to the manufacturer’s instructions. For western-blotting, 48 h post-transfection the medium was removed and the cells were washed twice with cold phosphate-buffered saline (PBS) and lysed with Lysis Buffer (50 mM Tris–HCl, pH 7.5, 150 mM NaCl, 1.5 mM MgCl2, 1.5 mM EGTA, 1 mM EDTA, 1% Triton X-100, 1 mM DTT, and Protease and Phosphatase Inhibitors). After 15 min on ice, the lysate was removed from the dishes and centrifuged at 12,000 g for 20 min at 4°C. Supernatants were quantified using Bradford method and kept until use. For RNA analysis cells were harvested and treated as described above.
Biotin-tagged miRNA complexes were purified as described before . Briefly, transfected SH-SY5Y cells were harvested 48 h after transfection, washed in PBS and scraped in 1 ml of lysis buffer (20 mM Tris; pH 7.5, 200 mM NaCl, 2.5 mM MgCl2, 0.05% Igepal, 60 U Superase-In/ml (Ambion), 1 mM DTT and Pefabloc (Roche). Lysate was vortexed twice for 20 seconds and cleared by spinning at 12,000g at 4 °C for 20 min in a microcentrifuge. Fifty μl of the sample were saved for input analysis and placed on dry ice. The remaining lysate was incubated for 1 h at 4 °C rotating with 25 μl pre-washed streptavidin-sepharose beads. In order to prevent non-specific binding of RNA and protein complexes, the beads were coated with RNase-free BSA and yeast tRNA (both from Ambion). Beads were washed three times with 500 μl ice-cold lysis buffer. In order to isolate the RNA, TRIzol was added to both input and affinity-washed beads. 40 μl of chloroform was added and after shaking samples were spin at 12,000 g at 4 °C for 15 min in a microcentrifuge. Supernatant was transferred to a fresh tube and RNA was precipitated adding 1 μl of GlycoBlue co-precipitant (Ambion), 6 μl of 5 M NaCl and 360 μl of ice-cold 95% ethanol overnight at −20 °C. Samples were pelleted at 12,000g for 15 min at 4°C, washed in 70% ethanol in RNase-free water and air-dried. RNA was dissolved in 120 μl of RNase-free water (Ambion). In order to isolate the proteins bound to the biotynilated miRNAs, washed beads were resuspended in SDS-buffer and boiled at 100°C.
The microRNA fraction from human frontal cortex lysate was purified using the microRNA isolation kit from Wako based on Ago2 immunoprepipitation complexes, following manufacturer’s instructions.
The protein content of either immunoprecipitations or cellular extracts were run on 8%-12% SDS-PAGE gels, transferred on to PVDF membranes, and incubated with the corresponding antibodies at the following concentrations: Ago-2 1:2000; Rck/p54: 1:1000; Mcl-1: 1:1000; COX-2: 1:500; EZH2: 1:1000; MKP-1: 1:250; APP: 1:1000; Tubulin: 1:10.000. The membranes were developed with the enhanced chemiluminescence method (ECL) following the manufacturer’s instructions.
SH-SY5Y cells were transfected with miR-101 or 5’-isomiR-101 as previously described in three independent experiments, and total RNA was isolated 48 h later. RNA concentration was measured with a Nanodrop™, and RNA quality was assessed by Bioanalyzer™ with RIN (RNA integrity number) ranging between 9 and 10. For each sample, 200 ng of total RNA was reverse transcribed, amplified by in vitro transcription and labeled with biotin-UTP using the Illumina Total Prep RNA amplification kit (IL1791, Applied Biosystem/Ambion, Austin, TX, USA) following the manufacturer’s instructions. 750 ng of biotinylated cRNA was hybridized in a BeadChip Hyb Chamber with rocking for 16 h at 58°C. Bead arrays were washed and blocked in a rocking incubator. Finally, bead arrays were dried by centrifugation followed by scanning in the Illumina Beadstation. The raw data were summarized per probe using BeadStudio software Gene Expression. The SAM (significance analysis of microarrays) two-class unpaired comparison test was applied to detect significant differences in gene expression between treated and control conditions, initially setting the statistical significance at a false discovery rate of 5%, with an arbitrary absolute fold chance cutoff set at 1.2.
Quantitative real time PCR was performed using the miRCURY LNA™ microRNA PCR System (Exiqon) on total RNA extracted from SH-SY5Y with mirVana’s isolation kit (Ambion) following the manufacturer’s instructions. PCR amplification and detection were performed with the ROCHE LightCycler 480 detector, using 2x SYBR GREEN Master Mix. The reaction profile was: Polymerase Activation/Denaturation cycle (95°C for 10 min) followed by 40 amplification cycles (95°C-10”, 60°C-20”). Each sample was amplified in triplicate. miRNA levels were calculated using the LightCycler 480 software. Both problem and calibrator samples were normalized by the relative expression of reference small non-coding RNAs (SNORD44 and SNORD48).
Total RNA was extracted from SH-SY5Y cells with mirVana isolation kit (Ambion). Purified RNAs were treated with RNase-free DNAse (DNA-free, Ambion) and reverse-transcribed, with Superscript II (Invitrogen) to generate the corresponding cDNAs that served as PCR templates for mRNA quantification. The primers used for RT-qPCR were (5’->3’): COX-2; F-CCTGTGCCTGATGATTGC, R-CTGATGCGTGAAGTGCTG, MCL-1; F-AAAGAGGCTGGGATGGGTTT, R-CAAAAGCCAGCAGCACATTC, APP; F-TCAGGTTGACGCCGCTGT, R-TTCGTAGCCGTTCTGCTGC, EZH2; F-AGTGTGACCCTGACCTCTGT, R-AGATGGTGCCAGCAATAGAT, DUSP-1 F-CAGCTGCTGCAGTTTGAGTC, R- AGAGGTCGTAATGGGGCTCT. PCR amplification and detection were performed with the ROCHE LightCycler 480 detector, using 2X SYBR GREEN Master Mix as reagent following the manufacturer’s instructions. The reaction profile was: denaturation/activation cycle (95° for 10 min) followed by 40 cycles of denaturation-annealing-extension (95°-10’, 60°-40”, 72°-10”) and a final melting cycle (95°-5’, 72°-10”, 98°-continuous). Each sample was amplified in triplicate. mRNA levels were calculated using the LightCycler 480 software. Samples were normalized by the relative expression of housekeeping gene (Tubulin).
Statistical analysis was performed with G-STAT software (http://www.gstat.org). The Mann–Whitney non-parametric test was used for the statistical analysis (two-tailed) of the RT-qPCR and western blot assays in transfection experiments; P≤ 0.05 was considered statistically significant.
In Figures 2A and 3C, results are expressed as the mean of the Fold Change ratios (microRNA - Biotin-mRNA / control siGLO - Biotin) ± SEM. In Figure 2B and 2D, results are expressed as the mean of the Fold Change ratios (FLAG-Ago2 Antibodies / control Antibodies) ± SEM. In Figure 2C, a miRNA-101 transfected sample was used to normalize data and assigned a value of 1 (Arbitrary Units). In Figure 3A and 3B, a control sample was used to normalize Fold Change values in each set of experiments and SEM was calculated. This control sample was assigned a value of 1.
RNA induced silencing complex.
The authors would like to thank Dra. Juana Diez (UPF-Barcelona), Dra. Carme Caelles (IRB-Barcelona) and Dr. Kristian Helin (BRIC-University of Copenhagen) for providing Rck/p54, MKP-1 and EZH2 antibodies respectively.
We also thank Michael Maudsley (Language Services – University of Barcelona) for correcting the English manuscript.
This work was supported by the Spanish Government and FEDER (Fondo Europeo de Desarrollo Regional): PN de I+D+I 2008–2011, PI081367 (E. Martí) and PN de I+D+I 2012–2015 PI11/02036 (E. Martí), PI1100968 (I. Ferrer) Instituto Carlos III –ISCIII-, Subdirección General de Evaluación y Fomento de la Investigación, SAF2008-00357 (X. Estivill) Ministerio de Economia y competitividad, BESAD (JA de Río), DEMTEST Joint Programing of Neurodegenerative diseases (EU) (JA del Río), MICINN BFU2009-10848 and BFU2012-32617 (JA del Río), FP7-PRIORITY (JA del Río), Eusko Fundazioa (BIO12/AL/004) (JA del Río) and Generalitat de Catalunya SGR2009-366 (JA del Río); the Sixth Framework Programme of the European Commission through the SIROCCO integrated project LSHG-CT-2006-037900. The Spanish Government supports: M. Bañez-Coronel (Sara Borrell Postdoctoral contract, ISCIII) and partially supports E.Martí (Programa Miguel Servet, ISCIII).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.