Selective extracellular vesicle-mediated export of an overlapping set of microRNAs from multiple cell types
BMC Genomics volume 13, Article number: 357 (2012)
MicroRNAs (miRNAs) are a class of small RNA molecules that regulate expression of specific mRNA targets. They can be released from cells, often encapsulated within extracellular vesicles (EVs), and therefore have the potential to mediate intercellular communication. It has been suggested that certain miRNAs may be selectively exported, although the mechanism has yet to be identified. Manipulation of the miRNA content of EVs will be important for future therapeutic applications. We therefore wished to assess which endogenous miRNAs are enriched in EVs and how effectively an overexpressed miRNA would be exported.
Small RNA libraries from HEK293T cells and vesicles before or after transfection with a vector for miR-146a overexpression were analysed by deep sequencing. A subset of miRNAs was found to be enriched in EVs; pathway analysis of their predicted target genes suggests a potential role in regulation of endocytosis. RT-qPCR in additional cell types and analysis of publicly available data revealed that many of these miRNAs tend to be widely preferentially exported. Whilst overexpressed miR-146a was highly enriched both in transfected cells and their EVs, the cellular:EV ratios of endogenous miRNAs were not grossly altered. MiR-451 was consistently the most highly exported miRNA in many different cell types. Intriguingly, Argonaute2 (Ago2) is required for miR-451 maturation and knock out of Ago2 has been shown to decrease expression of other preferentially exported miRNAs (eg miR-150 and miR-142-3p).
The global expression data provided by deep sequencing confirms that specific miRNAs are enriched in EVs released by HEK293T cells. Observation of similar patterns in a range of cell types suggests that a common mechanism for selective miRNA export may exist.
Membrane-bound extracellular vesicles (EVs), are referred to by various terms in the literature, including microvesicles or microparticles. However, they may be divided into two main types formed either by budding of the plasma membrane (shedding vesicles, ectosomes; approximately 100 nm −1 μm in diameter) or by exocytosis of multivesicular bodies (MVB) (exosomes; approximately 50–100 nm) [1–3]. While it is important to recognize that these two main populations of EVs have distinct biogenesis pathways, it remains technically challenging to isolate their respective pure populations due to their overlapping size distributions and most preparations are therefore heterogeneous . EVs are targeted to specific recipient cells by surface receptors which mediate their internalisation and uptake of their content, comprising lipids, proteins, mRNAs and microRNAs (miRNAs) [2, 4]. A significant physiological role for EVs is supported by their implication in cancer progression, angiogenesis, cell reprogramming and differentiation [5–9].
MiRNAs are a class of small RNA molecules that regulate gene expression by binding to specific mRNA targets and triggering their degradation or inhibiting their translation. Primary miRNAs (pri-miRNAs), transcribed either as independent genes or as introns of coding genes, are processed to form precursor hairpin loops (pre-miRNAs) and then mature miRNAs by the RNAse III enzymes Drosha and Dicer respectively. One strand of the mature miRNA is loaded into an RNA induced silencing complex (RISC), which mediates the interaction between the miRNA and its target mRNA molecules. Argonaute (Ago) proteins are part of the RISC and mediate cleavage of target mRNAs [10, 11].
Recently, miRNAs have been found extracellularly, being encapsulated within vesicles (EVs), or associated with Ago2 or NPM1 proteins or high-density lipoproteins (HDL) [12–19]. EV preparations isolated by ultracentrifugation are heterogeneous and include both exosomes and other vesicles. Extracellular miRNAs have emerged as prognostic markers for various diseases, such as cancer and cardiovascular disorders [13, 20–23]. Recent studies have demonstrated that miRNAs can be shuttled between cells in EVs and modulate gene expression [6, 24–26] and control different physiological functions in the recipients [27–30].
It has been shown that miRNA export into exosomes occurs via the ceramide-dependent secretory pathway and independently of the endosomal sorting complex required for transport (ESCRT) . Although the mechanism of release of shedding vesicles is less well known, the role of P4-ATPase and ESCRT has been recognized recently . Several studies have suggested that the miRNA content of EVs does not simply reflect the miRNA repertoire of the cells of origin and that some miRNAs are selectively exported or retained within the cell [25, 27, 32–34]. However, the mechanism whereby specific miRNAs are selectively exported is unknown.
The potential for exploiting extracellular vesicles as vehicles for delivery of therapeutic molecules has been appreciated [35–38]. For example, EVs prepared from a patient’s own cells could potentially be enriched with specific siRNAs and/or miRNAs and used as a tool for therapeutic delivery of these regulatory RNA molecules. Alternatively, the therapeutic effect could be achieved by transplanting cells engineered to over-express critical miRNAs, which are delivered to the target cells via EVs[33, 36]. To most effectively harness EVs as therapeutic tools, it will be crucial to understand and control how they are loaded with regulatory RNAs.
As a step towards understanding the mechanism of selective export of miRNAs to EVs we have applied deep sequencing to characterize the global expression pattern of small RNAs in HEK293T cells and the EVs that they release. Although the enrichment of overexpressed miRNA in EVs has been shown by RT-qPCR in HEK293T cells, mesenchymal stem cells, macrophages and immune cells [24, 25, 27, 28, 37], it is not known to what extent overexpressed miRNAs are exported into EVs and how they affect the export of endogenous miRNAs. To address these questions, small RNAs from cells transfected with miR-146a expressing plasmid and their EVs were also subjected to deep sequencing. Ectopic overexpression of miR-146a in HEK293T exosomes was previously shown to exert a functional effect upon recipient cells, including downregulation of target gene expression .
Comparison of our data with other publicly available data sets of miRNA expression in cells and their EVs revealed a subset of miRNAs that tend to be preferentially exported in various cell types, suggesting a common function in intercellular communication.
EVs from HEK293T cells were obtained by ultracentifugation of conditioned media at 100000 g. The presence of spherical membrane bound structures of variable diameter up to approximately 1000 nm in these preparations was confirmed by transmission electron microscopy (Figure 1A). RNA extracted from cells and EVs was analyzed on a Bioanalyzer (Agilent Technologies). The shape of the RNA profiles observed when comparing RNA from pSM30-miR-146a transfected cells and EVs and their non-transfected counterparts were similar (Figure 1B). Total RNA from EVs was enriched in small RNAs and lacked ribosomal RNA peaks typical for cell RNA profiles (although a low broad peak was observed at the level of 18 S RNA). Analysis of the small RNA profiles of EV samples revealed enrichment of sequences of the size expected for miRNAs, with a peak at approximately 20 nt.
To compare the global miRNA profile of non-transfected and pSM30-miR-146a transfected cells and their EVs, small RNA libraries were prepared and sequenced on an Illumina Genome Analyzer. Raw reads were subjected to adaptor removal and aligned to miRBase (Release 17, http://www.mirbase.org/) . Both raw and processed data has been deposited in Gene Expression Omnibus (GEO: http://www.ncbi.nlm.nih.gov/geo/), with accession number GSE38916. The total numbers of reads >15 bp matching each miRNA were calculated. These values were then normalized to number of reads per million mapped (RPMM) in each library, to enable direct comparison of expression levels. As a measure of the relative levels of expression of individual miRNAs in cells and EVs, the log (base two) of the ratio of the read number in EVs over that in cells was calculated (log2 (EV/cell)). The miRNA expression data for cells and EVs is presented in Additional file 1: Table S1.
The enrichment of miR-146a in cells and EVs following transfection with pSM30-miR-146a was striking: Mir-146a was enriched from 10 RPMM (0.1% of all reads) to 317,863 RPMM (32%) in cells and from 1,321 RPMM (1.3%) to 440,496 RPMM (44%) in EVs. Sequencing revealed that overexpressed miR-146a was correctly processed, with the most prevalent isomiR being the same in all cell and EV libraries (Figure 2A) and represented by approximately half of the reads. RT-qPCR from independent pSM30-miR-146a transfections confirmed consistently high overexpression of miR-146a in both cells and EVs (Figure 2B).
To evaluate how the overexpression of miR-146a affected the proportions of other, endogenous miRNAs, we excluded reads for miR-146a from transfected cells and EVs and re-normalized the expression values of all other miRNAs to reads per million mapped (RPMM). This enabled direct comparison between read numbers in transfected and non-transfected cells and EVs. This re-normalisation revealed extremely high correlation between miRNA expression values in non-transfected and transfected cells (R2 = 0.9) (Figure 3A). Moreover, high correlation was also observed between EVs from non-transfected and transfected cells (R2 = 0.8) (Figure 3B). However, correlation between both transfected and non-transfected cells and their respective EVs was much lower (R2 = 0.6) (Figure 3C,D). From all miRNAs that had been detected in both non-transfected cells and EVs, only 19.5% had log2 (EV/cell) between −1 and +1, therefore being represented by similar number of reads in cells and EVs.
The miRNAs preferentially retained within the cell or released in EVs are highly consistent even after over-expression; a scatter plot of log2 ratios of miRNAs represented by more than 1000 reads either in cells or EVs revealed relatively high correlation (R2 = 0.748)(Figure 4A). The most retained and most released miRNAs which are represented by more than 1000 reads either in cells or in EVs are listed in Tables 1 and 2 respectively. MiR-218 was found to be the most retained and miR-451 the most exported miRNA. If miRNAs are incorporated passively into EVs as a random sample of the cellular population one would expect the ratios of EV:cell expression values to be approximately 1:1 when both values are normalized to million reads mapped (i.e. a distribution of log2 (EV/cell) centred on zero). However, we observed a skewed distribution (Figure 4B), with more miRNAs represented by higher RPMM values in EVs than cells (i.e. with a positive log2 (EV/cell) ratio). The values of the positive ratios representing the most released miRNAs, were more extreme than the most negative values. However, several very highly expressed miRNAs had negative values (i.e. preferentially retained within the cell).
To confirm sequencing data, we performed RT-qPCR on RNA samples isolated from biological replicates of non-transfected and pSM30miR-146a transfected cells and their EVs. MiR-26a was chosen for normalization of the RT-qPCR data because this miRNA had relatively high expression (~2000 RPMM) and the numbers of reads in cells and EVs were similar (EV:Cell ratio ~1). The miRNA EV:cell ratios obtained by PCR correlated with sequencing data and confirmed retained and exported miRNAs (Figure 5A).
The most highly exported miRNA in EVs in our data, miR-451, was also reported previously to be highly exported from breast cancer cells lines and mesenchymal stem cells [25, 33]. This was intriguing, and it prompted us to speculate that the same miRNAs may be preferentially exported or retained in various cell types. We therefore tested two other cell types, namely a human microvascular endothelial cell line, (HMEC-1) and primary outgrowth endothelial progenitor cells (OECs) with the same panel of primers used to confirm sequencing data in HEK293T cells. For the 12 miRNAs expressed in both endothelial cells and EVs the proportions retained or released followed a similar pattern to that observed in HEK293T cells (Figure 5B). Two miRNAs, miR-1246 and miR-9 were not detected and miR-7 was not detected in endothelial EVs.
To obtain a broader view of miRNA export in various cells types, we retrieved publicly available data (microarray and RT-qPCR arrays) on miRNA expression in various cell types and their EVs (microvesicles, exosomes). The ratios (or delta Ct for RT-qPCR data) between the expression levels in EVs and cells were calculated and the miRNAs ranked accordingly. Comparison of the ranks of the miRNAs present in our data with their ranks in the other datasets revealed correlations between the different cell types. This was visualized by assigning a color for each miRNA according to its rank in a particular data set: Figure 6A shows the ranks of the top and bottom 100 miRNAs sorted according to the data from HEK293T non-transfected cells (all miRNAs are listed in Additional file 2: Table S2). MiR-451, the most highly enriched miRNA in HEK293T EVs, was among the top 6 exported miRNAs in 6 out of 8 of the other datasets.
MiR-451 is the only reported miRNA that is not processed by Dicer, but instead uses an Ago2-mediated pathway for its maturation [40–42]. It is tempting to speculate that this unique processing pathway might be connected with the high extracellular export of miR-451. In this case one might expect the most exported miRNAs to also be the most affected by loss of Ago2 activity. We therefore ranked miRNAs according to the reduction in their expression in bone marrow of Ago2 ko mice compared to the wild type . Remarkably, in addition to miR-451, several other miRNAs among the top 6 affected by the Ago2 deletion (miR150, miR-486, and miR-144) were also found to be highly exported in several data sets (see column (K) in Figure 6).
As an initial step towards understanding the function of the preferentially exported miRNAs we attempted to determine the genes which are being targeted by these miRNAs. It is difficult to predict these interactions with certainty, but the algorithms available are improving. We employed microT-CDS, which has recently been shown to perform well in terms of sensitivity and specificity compared to other available algorithms . The 2760 genes predicted to be targeted by one or more of the ten most selectively exported miRNAs are listed in Additional file 3: Table S3. Functional analysis of these genes revealed that the most enriched pathway was endocytosis (Kegg hsa04144) (Benjamini-corrected p-value = 0.05) which contained 43 predicted target genes (Figure 7). All the enriched pathways and gene ontologies are listed in Additional file 4: Table S4 as well as clusters of functionally related terms, which include transcriptional regulation, intracellular transport and organelles. Recently developed techniques based upon crosslinking of RISC proteins with their target mRNAs have facilitated experimental determination of miRNA:target interactions [45, 46]. Therefore, to complement the in silico target prediction, the genes with target sites for the same 10 most exported miRNAs predicted by any algorithm that are supported experimentally by CLIP-Seq reads available in Starbase  were downloaded (Additional file 3: Table S3). These genes were also significantly enriched in the endocytosis pathway (Benjamini-corrected p-value = 0.0009) and because only three of the predicted target genes in this pathway are common (Figure 7) the in silico and experimental predictions provide largely independent associations with endocytosis. From the total of 43 microT-CDS and 39 Starbase predicted target genes 7 were common.
MiRNAs are involved in regulation of a wide array of biological processes, from development to immunity. Extracellular and circulating miRNAs have been found both in the non-vesicular fraction and encapsulated in EVs. They have been shown to associate with lipids (HDL) and proteins (Ago2 and NPM1). While the role of EV miRNAs in cell-to-cell communication has now been widely accepted  and HDL-mediated miRNA transfer has been shown to regulate recipient cell gene expression , there is still no direct evidence that miRNAs associated with proteins act in a paracrine fashion [16–18].
Realisation of the potential importance of miRNAs in cell-to-cell communication has triggered research into the mechanisms of miRNA release and their effects upon recipient cells . Evidence is accumulating that miRNA export into EVs is selective, but it is not known how this export is regulated. Several studies have compared the miRNA content of cells and EVs using microarrays or RT-qPCR. Deep sequencing is an alternative approach for measuring miRNA expression which has the advantage of providing a digital measure that does not require an internal standard for normalisation. This is particularly important when comparing miRNA levels between cells and EVs, where no common internal standards have been identified. Although several previous studies have applied deep sequencing in similar scenarios, they involved fewer reads and the expression of only selected miRNAs was reported [48, 49]. Here we report a global comparison of the miRNA pools in HEK293T cells and their EVs. Our analysis confirmed previous observations that some miRNAs are preferentially exported and found to be enriched in EVs.
Enrichment of a miRNA in EVs after its overexpression in cells has been measured previously by RT-qPCR. However, it is important to reveal how such overexpression compares with and affects expression and export of endogenous miRNAs. We found overexpressed miR-146a to contribute almost a third of all reads in cells and even more in EVs, representing almost half of their cargo. It is very important for the miRNA therapeutic field to understand the consequence of miRNA overexpression. Our data suggest that if a miRNA is overexpressed in a specific cell type the possibility that a significant amount of this miRNA will be exported and potentially affect other cells should be considered. Although we cannot assess the effect of miR-146a overexpression upon the absolute expression of other miRNAs, this effect is likely to be similar in both cells and EVs. However, our data demonstrates that the relative expression levels of endogenous miRNAs (and their EV:cell miRNA ratios) remain very stable and are largely unaffected by the miRNA overexpression. The fact that selective miRNA export is tightly regulated suggests that it plays an important role in cell physiology. This will be an important consideration not only for further studies into the mechanism of miRNA export but also for the development of methods for enrichment of EVs with miRNA for therapies. It is possible that the overexpressed miRNA and endogenous miRNAs use different export routes and therefore do not affect each other. Since our EV preparation does not distinguish between exosomes and shedding vesicles we cannot exclude that one of these is preferentially used to carry overexpressed miRNA. Fractionation of EVs will be necessary to further examine this possibility. Nevertheless, there may be an overlap between the two pathways, as a recent study has suggested that shedding vesicles and exosomes require the same GTPase (RAB-11) and ESCRT proteins for their biogenesis and therefore may be regulated by the same mechanisms, regardless of their distinct cellular localization .
Comparison of our data with publicly available expression sets from other cell types revealed that some miRNAs tend to be universally exported. It is possible that this reflects release of the same type of vesicles (eg exosome) by all the cell types analysed. The commonly exported miRNAs possibly fulfill a common function in intercellular communication. MiR-451 was found to be the miRNA with the highest ratio between expression levels in EVs and cells. In relation to its maturation miR-451 is a non-canonical miRNA, processed by the Ago2 protein, in contrast to other miRNAs that commonly require Dicer for this process [40, 41, 43]. Ago2 has been shown previously to be associated with both EV and non-vesicular miRNAs [17, 18, 49]. Comparison of miRNA expression between cells and their EV counterparts from other published data sets reveals miR-451 to be one of the most highly exported miRNAs in multiple cell types. Deletion of Ago2 in bone marrow leads to reduced expression of several other miRNAs in addition to miR-451 , notably some of the miRNAs that are most preferentially exported from various different cell types. These data suggest that there may be a connection between Ago2 processing and the selection of miRNAs for export. Recent studies have demonstrated the physical colocalisation of RISC and MVBs and their functional interactions in regulating miRNA loading and RNA silencing [49, 50]. Being a key component of RISC, Ago2 may potentially be involved in directing miRNA towards exosomal secretion. However it has not been investigated whether Ago2 mediated miR-451 processing also occurs within RISC and in colocalisation with MVB membrane, or in another subcellular compartment. Future studies are required to characterize miRNA export into specific EV populations and must also address the challenge of dissecting the function of Ago2 in miRNA processing from its potential role in miRNA export.
Specific endogenous miRNAs, most notably miR-451 are selectively released from HEK293T cells in EVs by a process that may involve Ago2. Conservation amongst the miRNAs exported by different cells suggests that this process may be widely used to regulate similar target cell functions. Enrichment in genes involved in endocytosis amongst the targets of the most selectively exported miRNAs suggests that these functions may include intercellular vesicular transport itself. The demonstration that EV miRNA content can be manipulated by cellular overexpression of a specific miRNA bodes well for future therapeutic use of this system to deliver miRNAs/siRNAs to target cells.
HEK293T cells from the American Type Culture Collection (ATCC) were cultured in DMEM (Life technologies, Paisley, UK) supplemented with 10% fetal calf serum (FCS) and 100 μg/ml Primocin™ (InvivoGen). Human microvascular endothelial cells (HMEC-1)  were cultured in MCDB 131 medium (Life technologies) containing 10% FCS, L-glutamine (2 mM), epidermal growth factor (10 ng/ml, Life technologies), hydrocortisone (1 μg/ml, Life technologies), and 100 μg/ml Primocin. As described previously , OECs were obtained from the mononuclear cell fraction of umbilical cord blood, obtained under full ethical approval from The Office for Research Ethics Committees Northern Ireland (ORECNI). Cells were resuspended in EGM-2 medium (Lonza) supplemented with 10% FCS and seeded on 24 well culture plates pre-coated with rat tail collagen type 1 (BD Biosciences) at a density of 1×107/ml All cells were cultured under standard conditions at 5% CO2, 37°C, and 95% humidity.
The mature miRNA sequences for miR-146a, were inserted into pSM30 vector according to the published protocol . Briefly, the oligos containing the miR-146a mature sequence and the miR-30 loop sequence were designed as follows: top strand miR-146a oligo: 5'- AGCGGACCCATGGAATTCAGTTCTCATAGTGAAGCCACAGATGTATGAGAACTGAATTCCATGGGTT-3' and the bottom strand miR-146a oligo: 5'-GGCAAACCCATGGAATTCAGTTCTCATACATCTGTGGCTTCACTATGAGAACTGAATTCCATGGGTC-3'). Oligos were synthesized and HPLC purified by Eurogentec (Seraing, Belgium). Five μl of each oligo (100 μM stock solution) were mixed with 8 μl water and 2 μl of 10 x restriction buffer 3 (New England Biolabs), heated at 95°C for 4 minutes and left at room temperature for 15 minutes to anneal. Annealed oligos were diluted 2500 x to get 10 nM for ligation. Annealed oligos were ligated into pSM30 vector using Rapid Ligation Kit (Roche) with 15:1 insert:vector ratio. pSM30 empty vector was digested with BSMBI restriction enzyme and linearized plasmid was gel purified. DH5α competent cells were transformed with the ligation reaction and positive colonies were selected on agar plates supplemented with kanamycin.
HEK293T cells were seeded at 3×104 cells/ml into T75 flasks and transfected after 24 hours with Turbofect in vitro transfection reagent (Fermentas, York, UK) according to the manufacturer’s instructions. After 16 hours media were replaced with media containing 10% serum depleted of EVs by centrifugation for 2 hours at 100000 g. EVs were collected from the conditioned media after 48 h hours.
Isolation and characterisation of EVs
For Electron microscopy HEK293T EVs pellets were fixed in 2.5% glutaraldehyde in 0.1 M sodium cacodylate buffer at 4°C overnight. After post-fixation in 1% osmium tetroxide, pellets were dehydrated with increasing alcohol concentrations and embedded in Spurr resin. Ultrathin sections were prepared using an ultramicrotome, placed on copper grids (Agar Scientific Ltd, UK), stained with uranyl acetate and lead citrate, and examined using a JEOL 100 CX transmission electron microscope.
RNA was extracted using a miRNeasy kit (Qiagen, Crawley, UK) or High pure miRNA isolation kit (Roche, Burgess Hill, UK) according to the manufacturers’ protocols. RNA purity and quantity was determined using a Nanodrop spectrophotometer (Thermo Scientific) or Qubit fluorimeter (Life Technologies). RNA integrity was assessed using RNA 2000 and small RNA chips on a Bioanalyzer (Agilent).
Small RNA libraries were prepared using a sample prep kit v1.5 kit (Illumina) following the manufacturer’s protocol. Cluster generation and sequencing on a Genome Analyzer II was performed at the Trinity Genome Sequencing Laboratory, Dublin (http://www.medicine.tcd.ie/sequencing). Sequencing data were analyzed using Genomics workbench software (CLCbio, Aarhus, Denmark). After removal of adapter sequences reads >15 bp were aligned to miRBase (Release 17) for annotation, allowing no more than 2 mismatches. Reads were normalized to reads per million mapped (RPMM) to enable comparison between different libraries. If a miRNA was found in both EV data sets but was detected in only one data set from cells (due to its low expression value) it was assigned one read in the other cell data set for the purpose of calculating an EV:cells ratio.
Target genes were predicted for miRNAs of interest using the microT-CDS algorithm available at http://www.microrna.gr/microT-CDS with a strict miTG Score Threshold of 0.8 . Lists were integrated using Microsoft Access and pathway analysis of the KEGG database performed using the The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 available at http://david.abcc.ncifcrf.gov/. Predicted target sites supported by CLIP-Seq reads were downloaded from Starbase (http://starbase.sysu.edu.cn/) .
Specific miRNAs were amplified using a modified version of the method described by Shi and Chiang . 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. RNA was then reverse transcribed with 200 U of SuperScript III reverse transcriptase (Life technologies, Paisley, UK) and 0.5 μg poly (T) adapter. Primers for specific miRNAs were designed from miRNA sequences obtained from miRBase and our sequencing data and synthesized by Eurogentec (Seraing, Liège, Belgium) (Additional file 5: Table S5). The reverse primer was the 3' adapter primer (3' RACE outer primer in the FirstChoice RLM-RACE kit; Ambion).
RT-qPCR was performed using Maxima SYBR Green qPCR mastermix (Fermentas) in 10 μl reactions containing 2 μl of 1:15 cDNA dilution. Reactions were performed on a LightCycler 480 (Roche), with the initial preincubation at 50°C for 2 min and activation at 95°C for 10 min, followed by 40 cycles at 95°C for 15 s, and 60°C for 60 s, with fluorescence acquired after 15 s of the 60°C step. Primer efficiencies were determined using serial dilutions of pooled cDNA and the reactions were optimized to achieve efficiencies of ~2 for all primers used. The RNA from three independent experiments was analyzed and all PCR reactions were performed in triplicates. Gene expression data were normalized to miR-26a. The relative expression was determined as 2−ΔCt, where ΔCt = Ct(miRNA)-Ct(miR-26a).
Analysis of publicly available miRNA expression data
Publicly available data were included in analysis if they had expression data for EVs (exosomes, microvesicles) and corresponding cells. Microarray data were downloaded from GEO database [27, 28, 33] or retrieved from supplementary data [58, 59]. Difference between expression in cells and EVs was determined as the ratio of expression EV:cells for microarrays or as delta Ct for RT-qPCR arrays and miRNAs were ranked accordingly. Using Microsoft Access miRNAs that are present in our data were pulled out from all the other data sets with their own ranks kept. Colours were assigned to miRNAs according to their ranks using Microsoft Excel (conditional formatting).
Crosslinking immunoprecipitation sequencing
Human microvascular endothelial cell line
Outgrowth endothelial progenitor cells
RNA-Induced Silencing Complex
Reads per million mapped.
Camussi G, Deregibus MC, Bruno S, Grange C, Fonsato V, Tetta C: Exosome/microvesicle-mediated epigenetic reprogramming of cells. Am J Cancer Res. 2011, 1 (1): 98-110.
Cocucci E, Racchetti G, Meldolesi J: Shedding microvesicles: artefacts no more. Trends Cell Biol. 2009, 19 (2): 43-51. 10.1016/j.tcb.2008.11.003.
Gyorgy B, Szabo TG, Pasztoi M, Pal Z, Misjak P, Aradi B, Laszlo V, Pallinger E, Pap E, Kittel A, Nagy G, Falus A, Buzas EI: Membrane vesicles, current state-of-the-art: emerging role of extracellular vesicles. Cell Mol Life Sci. 2011, 68 (16): 2667-2688. 10.1007/s00018-011-0689-3.
Chen X, Liang H, Zhang J, Zen K, Zhang CY: Secreted microRNAs: a new form of intercellular communication. Trends Cell Biol. 2012, 22 (3): 125-132. 10.1016/j.tcb.2011.12.001.
Anderson HC, Mulhall D, Garimella R: Role of extracellular membrane vesicles in the pathogenesis of various diseases, including cancer, renal diseases, atherosclerosis, and arthritis. Lab Invest. 2010, 90 (11): 1549-1557. 10.1038/labinvest.2010.152.
Zhang Y, Liu D, Chen X, Li J, Li L, Bian Z, Sun F, Lu J, Yin Y, Cai X, Sun Q, Wang K, Ba Y, Wang Q, Wang D, Yang J, Liu P, Xu T, Yan Q, Zhang J, Zen K, Zhang CY: Secreted monocytic miR-150 enhances targeted endothelial cell migration. Mol Cell. 2010, 39 (1): 133-144. 10.1016/j.molcel.2010.06.010.
Quesenberry PJ, Dooner MS, Aliotta JM: Stem cell plasticity revisited: the continuum marrow model and phenotypic changes mediated by microvesicles. Exp Hematol. 2010, 38 (7): 581-592. 10.1016/j.exphem.2010.03.021.
Muralidharan-Chari V, Clancy JW, Sedgwick A, D'Souza-Schorey C: Microvesicles: mediators of extracellular communication during cancer progression. J Cell Sci. 2010, 123 (Pt 10): 1603-1611.
Yuan A, Farber EL, Rapoport AL, Tejada D, Deniskin R, Akhmedov NB, Farber DB: Transfer of microRNAs by embryonic stem cell microvesicles. PLoS One. 2009, 4 (3): e4722-10.1371/journal.pone.0004722.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
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.
Hata T, Murakami K, Nakatani H, Yamamoto Y, Matsuda T, Aoki N: Isolation of bovine milk-derived microvesicles carrying mRNAs and microRNAs. Biochem Biophys Res Commun. 2010, 396 (2): 528-533. 10.1016/j.bbrc.2010.04.135.
Kosaka N, Iguchi H, Ochiya T: Circulating microRNA in body fluid: a new potential biomarker for cancer diagnosis and prognosis. Cancer Sci. 2010, 101 (10): 2087-2092. 10.1111/j.1349-7006.2010.01650.x.
Rabinowits G, Gercel-Taylor C, Day JM, Taylor DD, Kloecker GH: Exosomal microRNA: a diagnostic marker for lung cancer. Clin Lung Cancer. 2009, 10 (1): 42-46. 10.3816/CLC.2009.n.006.
Skog J, Wurdinger T, van Rijn S, Meijer DH, Gainche L, Sena-Esteves M, Curry WT, Carter BS, Krichevsky AM, Breakefield XO: Glioblastoma microvesicles transport RNA and proteins that promote tumour growth and provide diagnostic biomarkers. Nat Cell Biol. 2008, 10 (12): 1470-1476. 10.1038/ncb1800.
Wang K, Zhang S, Weber J, Baxter D, Galas DJ: Export of microRNAs and microRNA-protective protein by mammalian cells. Nucleic Acids Res. 2010, 38 (20): 7248-7259. 10.1093/nar/gkq601.
Turchinovich A, Weiz L, Langheinz A, Burwinkel B: Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011, 39 (16): 7223-7233. 10.1093/nar/gkr254.
Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, Mitchell PS, Bennett CF, Pogosova-Agadjanyan EL, Stirewalt DL, Tait JF, Tewari M: Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci U S A. 2011, 108 (12): 5003-5008. 10.1073/pnas.1019055108.
Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT: MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011, 13 (4): 423-433. 10.1038/ncb2210.
Taylor DD, Gercel-Taylor C: MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol Oncol. 2008, 110 (1): 13-21. 10.1016/j.ygyno.2008.04.033.
Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, Guo J, Zhang Y, Chen J, Guo X, Li Q, Li X, Wang W, Zhang Y, Wang J, Jiang X, Xiang Y, Xu C, Zheng P, Zhang J, Li R, Zhang H, Shang X, Gong T, Ning G, Wang J, Zen K, Zhang J, Zhang CY: Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008, 18 (10): 997-1006. 10.1038/cr.2008.282.
Contu R, Latronico MV, Condorelli G: Circulating microRNAs as potential biomarkers of coronary artery disease: a promise to be fulfilled?. Circ Res. 2010, 107 (5): 573-574. 10.1161/CIRCRESAHA.110.227983.
Gupta SK, Bang C, Thum T: Circulating microRNAs as biomarkers and potential paracrine mediators of cardiovascular disease. Circ Cardiovasc Genet. 2010, 3 (5): 484-488. 10.1161/CIRCGENETICS.110.958363.
Kosaka N, Iguchi H, Yoshioka Y, Takeshita F, Matsuki Y, Ochiya T: Secretory mechanisms and intercellular transfer of microRNAs in living cells. J Biol Chem. 2010, 285 (23): 17442-17452. 10.1074/jbc.M110.107821.
Collino F, Deregibus MC, Bruno S, Sterpone L, Aghemo G, Viltono L, Tetta C, Camussi G: Microvesicles derived from adult human bone marrow and tissue specific mesenchymal stem cells shuttle selected pattern of miRNAs. PLoS One. 2010, 5 (7): e11803-10.1371/journal.pone.0011803.
Iguchi H, Kosaka N, Ochiya T: Secretory microRNAs as a versatile communication tool. Commun Integr Biol. 2010, 3 (5): 478-481. 10.4161/cib.3.5.12693.
Mittelbrunn M, Gutierrez-Vazquez C, Villarroya-Beltri C, Gonzalez S, Sanchez-Cabo F, Gonzalez MA, Bernad A, Sanchez-Madrid F: Unidirectional transfer of microRNA-loaded exosomes from T cells to antigen-presenting cells. Nat Commun. 2011, 2: 282-
Montecalvo A, Larregina AT, Shufesky WJ, Beer Stolz D, Sullivan ML, Karlsson JM, Baty CJ, Gibson GA, Erdos G, Wang Z, Milosevic J, Tkacheva OA, Divito SJ, Jordan R, Lyons-Weiler J, Watkins SC, Morelli AE: Mechanism of transfer of functional microRNAs between mouse dendritic cells via exosomes. Blood. 2012, 119 (3): 756-766. 10.1182/blood-2011-02-338004.
Hergenreider E, Heydt S, Treguer K, Boettger T, Horrevoets AJ, Zeiher AM, Scheffer MP, Frangakis AS, Yin X, Mayr M, Braun T, Urbich C, Boon RA, Dimmeler S: Atheroprotective communication between endothelial cells and smooth muscle cells through miRNAs. Nat Cell Biol. 2012, 14 (3): 249-256. 10.1038/ncb2441.
Cantaluppi V, Gatti S, Medica D, Figliolini F, Bruno S, Deregibus MC, Sordi A, Biancone L, Tetta C, Camussi G: Microvesicles derived from endothelial progenitor cells protect the kidney from ischemia-reperfusion injury by microRNA-dependent reprogramming of resident renal cells. Kidney Int. 2012, 82 (4): 412-27. 10.1038/ki.2012.105.
Wehman AM, Poggioli C, Schweinsberg P, Grant BD, Nance J: The P4-ATPase TAT-5 inhibits the budding of extracellular vesicles in C. elegans embryos. Curr Biol. 2011, 21 (23): 1951-1959. 10.1016/j.cub.2011.10.040.
Chen TS, Lai RC, Lee MM, Choo AB, Lee CN, Lim SK: Mesenchymal stem cell secretes microparticles enriched in pre-microRNAs. Nucleic Acids Res. 2010, 38 (1): 215-224. 10.1093/nar/gkp857.
Pigati L, Yaddanapudi SC, Iyengar R, Kim DJ, Hearn SA, Danforth D, Hastings ML, Duelli DM: Selective release of microRNA species from normal and malignant mammary epithelial cells. PLoS One. 2010, 5 (10): e13515-10.1371/journal.pone.0013515.
Ohshima K, Inoue K, Fujiwara A, Hatakeyama K, Kanto K, Watanabe Y, Muramatsu K, Fukuda Y, Ogura S, Yamaguchi K, Mochizuki T: Let-7 microRNA family is selectively secreted into the extracellular environment via exosomes in a metastatic gastric cancer cell line. PLoS One. 2010, 5 (10): e13247-10.1371/journal.pone.0013247.
Lai RC, Chen TS, Lim SK: Mesenchymal stem cell exosome: a novel stem cell-based therapy for cardiovascular disease. Regen Med. 2011, 6 (4): 481-492. 10.2217/rme.11.35.
Alvarez-Erviti L, Seow Y, Yin H, Betts C, Lakhal S, Wood MJ: Delivery of siRNA to the mouse brain by systemic injection of targeted exosomes. Nat Biotechnol. 2011, 29 (4): 341-345. 10.1038/nbt.1807.
Akao Y, Iio A, Itoh T, Noguchi S, Itoh Y, Ohtsuki Y, Naoe T: Microvesicle-mediated RNA molecule delivery system using monocytes/macrophages. Mol Ther. 2011, 19 (2): 395-399. 10.1038/mt.2010.254.
van Dommelen SM, Vader P, Lakhal S, Kooijmans SA, van Solinge WW, Wood MJ, Schiffelers RM: Microvesicles and exosomes: Opportunities for cell-derived membrane vesicles in drug delivery. J Control Release. 2011, 161 (2): 635-44.
Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database issue): D152-D157.
Cheloufi S, Dos Santos CO, Chong MM, Hannon GJ: A dicer-independent miRNA biogenesis pathway that requires Ago catalysis. Nature. 2010, 465 (7298): 584-589. 10.1038/nature09092.
Cifuentes D, Xue H, Taylor DW, Patnode H, Mishima Y, Cheloufi S, Ma E, Mane S, Hannon GJ, Lawson ND, Wolfe SA, Giraldez AJ: A novel miRNA processing pathway independent of Dicer requires Argonaute2 catalytic activity. Science. 2010, 328 (5986): 1694-1698. 10.1126/science.1190809.
Yang JS, Lai EC: Dicer-independent, Ago2-mediated microRNA biogenesis in vertebrates. Cell Cycle. 2010, 9 (22): 4455-4460. 10.4161/cc.9.22.13958.
Yang JS, Maurin T, Robine N, Rasmussen KD, Jeffrey KL, Chandwani R, Papapetrou EP, Sadelain M, O'Carroll D, Lai EC: Conserved vertebrate mir-451 provides a platform for Dicer-independent, Ago2-mediated microRNA biogenesis. Proc Natl Acad Sci U S A. 2010, 107 (34): 15163-15168. 10.1073/pnas.1006432107.
Reczko M, Maragkakis M, Alexiou P, Grosse I, Hatzigeorgiou AG: Functional microRNA targets in protein coding sequences. Bioinformatics. 2012, 28 (6): 771-776. 10.1093/bioinformatics/bts043.
Hafner M, Landthaler M, Burger L, Khorshid M, Hausser J, Berninger P, Rothballer A, Ascano M, Jungkamp AC, Munschauer M, Ulrich A, Wardle GS, Dewell S, Zavolan M, Tuschl T: Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell. 2010, 14 (1): 129-141.
Chi SW, Zang JB, Mele A, Darnell RB: Argonaute HITS-CLIP decodes microRNA-mRNA interaction maps. Nature. 2009, 460 (7254): 479-486.
Yang JH, Li JH, Shao P, Zhou H, Chen YQ, Qu LH: starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 2011, 39 (Database issue): D202-D209.
Koh W, Sheng CT, Tan B, Lee QY, Kuznetsov V, Kiang LS, Tanavde V: Analysis of deep sequencing microRNA expression profile from human embryonic stem cells derived mesenchymal stem cells reveals possible role of let-7 microRNA family in downstream targeting of hepatic nuclear factor 4 alpha. BMC Genomics. 2010, 11 (Suppl 1): S6-10.1186/1471-2164-11-S1-S6.
Gibbings DJ, Ciaudo C, Erhardt M, Voinnet O: Multivesicular bodies associate with components of miRNA effector complexes and modulate miRNA activity. Nat Cell Biol. 2009, 11 (9): 1143-1149. 10.1038/ncb1929.
Lee YS, Pressman S, Andress AP, Kim K, White JL, Cassidy JJ, Li X, Lubell K, do Lim H, Cho IS, Nakahara K, Preall JB, Bellare P, Sontheimer EJ, Carthew RW: Silencing by small RNAs is linked to endosomal trafficking. Nat Cell Biol. 2009, 11 (9): 1150-1156. 10.1038/ncb1930.
Ades EW, Candal FJ, Swerlick RA, George VG, Summers S, Bosse DC, Lawley TJ: HMEC-1: establishment of an immortalized human microvascular endothelial cell line. J Invest Dermatol. 1992, 99 (6): 683-690. 10.1111/1523-1747.ep12613748.
Medina RJ, O'Neill CL, Humphreys MW, Gardiner TA, Stitt AW: Outgrowth endothelial cells: characterization and their potential for reversing ischemic retinopathy. Invest Ophthalmol Vis Sci. 2010, 51 (11): 5906-5913. 10.1167/iovs.09-4951.
Du G, Yonekubo J, Zeng Y, Osisami M, Frohman MA: Design of expression vectors for RNA interference based on miRNAs and RNA splicing. FEBS J. 2006, 273 (23): 5421-5427. 10.1111/j.1742-4658.2006.05534.x.
Deregibus MC, Cantaluppi V, Calogero R, Lo Iacono M, Tetta C, Biancone L, Bruno S, Bussolati B, Camussi G: Endothelial progenitor cell derived microvesicles activate an angiogenic program in endothelial cells by a horizontal transfer of mRNA. Blood. 2007, 110 (7): 2440-2448. 10.1182/blood-2007-03-078709.
Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38 (Database issue): D355-D360.
da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.
Shi R, Chiang VL: Facile means for quantifying microRNA expression by real-time PCR. Biotechniques. 2005, 39 (4): 519-525. 10.2144/000112010.
Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO: Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007, 9 (6): 654-659. 10.1038/ncb1596.
Hunter MP, Ismail N, Zhang X, Aguda BD, Lee EJ, Yu L, Xiao T, Schafer J, Lee ML, Schmittgen TD, Nana-Sinkam SP, Jarjoura D, Marsh CB: Detection of microRNA expression in human peripheral blood microvesicles. PLoS One. 2008, 3 (11): e3694-10.1371/journal.pone.0003694.
This work was supported by grant number BB/H005498/1from the Biotechnology and Biological Sciences Research Council (BBSRC). The authors thank Gerry Mahon and Liza Colhoun for assistance with electron microscopy and Dr. G. Du (University of Texas Health Science Center at Houston, TX, USA) for kindly providing pSM30. The Trinity Genome Sequencing Laboratory is a core facility funded by Science Foundation Ireland and supported by the Trinity Centre for High Performance Computing.
The authors declare that they have no competing interests.
JG-F contributed to the design of the study, carried out library preparation, sample preparation and RT-qPCR, participated in the analysis of sequence data and drafted the manuscript. AOC constructed the miR-146a expression vector. BC performed RT-qPCR expression analyses. CLO isolated OECs. RJM isolated OECs and assisted with manuscript preparation. DAS conceived, designed and co-ordinated the study, analysed the sequence data and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Table S2. Full list of ranked miRNAs shown in Figure 6. (XLSX 27 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
Guduric-Fuchs, J., O’Connor, A., Camp, B. et al. Selective extracellular vesicle-mediated export of an overlapping set of microRNAs from multiple cell types. BMC Genomics 13, 357 (2012). https://doi.org/10.1186/1471-2164-13-357