Author Response: KRAS-dependent Sorting of Mirna to Exosomes
Diana J Cha,Jeffrey L Franklin,Yongchao Dou,Qi Liu,James N Higginbotham,Michelle Demory Beckler,Alissa M Weaver,Kasey Vickers,Nirpesh Prasad,Shawn Levy,Bing Zhang,Robert J Coffey,James G Patton
DOI: https://doi.org/10.7554/elife.07197.027
2015-01-01
Abstract:Article Figures and data Abstract eLife digest Introduction Results Discussion Materials and methods References Decision letter Author response Article and author information Metrics Abstract Mutant KRAS colorectal cancer (CRC) cells release protein-laden exosomes that can alter the tumor microenvironment. To test whether exosomal RNAs also contribute to changes in gene expression in recipient cells, and whether mutant KRAS might regulate the composition of secreted microRNAs (miRNAs), we compared small RNAs of cells and matched exosomes from isogenic CRC cell lines differing only in KRAS status. We show that exosomal profiles are distinct from cellular profiles, and mutant exosomes cluster separately from wild-type KRAS exosomes. miR-10b was selectively increased in wild-type exosomes, while miR-100 was increased in mutant exosomes. Neutral sphingomyelinase inhibition caused accumulation of miR-100 only in mutant cells, suggesting KRAS-dependent miRNA export. In Transwell co-culture experiments, mutant donor cells conferred miR-100-mediated target repression in wild-type-recipient cells. These findings suggest that extracellular miRNAs can function in target cells and uncover a potential new mode of action for mutant KRAS in CRC. https://doi.org/10.7554/eLife.07197.001 eLife digest Cells use several different methods to control which genes are expressed to produce the proteins and RNA molecules that they need to work efficiently. The first step of gene expression is to transcribe a gene to form an RNA molecule. Protein-coding mRNA molecules can then be translated to make proteins. However, many RNA transcripts do not encode proteins. One example of these non-coding RNAs is a class of small RNAs called microRNAs (miRNAs), which are predicted to target more than 60% of protein-coding genes and can control which proteins are made. It was once thought that miRNAs exist only within the cell where they are synthesized. Recently, however, miRNAs have been found outside the cell bound to lipids and proteins, or encased in extracellular vesicles, such as exosomes. Exosomes are small bubble-like structures used by cells to export material into the space outside of cells. Exosomes containing miRNAs can circulate throughout the body, potentially transferring information between cells to alter gene expression in recipient cells. Many colorectal cancer cells have mutations in a gene that encodes a protein called KRAS. In 2011 and 2013, researchers found that the contents of the exosomes released from these mutant KRAS colorectal cancer cells can influence normal cells in ways that would help a cancer to spread. Furthermore, the exosomes released from the KRAS mutant cells contain different proteins than non-mutant cells. Now, Cha, Franklin et al.—including several researchers who worked on the 2011 and 2013 studies—show that exosomes released by mutant KRAS cells also contain miRNAs, and that these miRNAs are different from the ones exported in exosomes by cells with a normal copy of the KRAS gene. In particular, several miRNAs that suppress cancer growth in a healthy cell are found at lower levels in mutant KRAS cells. Instead, these miRNAs are highly represented in the exosomes that are released by the KRAS mutant cells. When cells with a normal copy of the KRAS gene were exposed to the contents of the exosomes released from KRAS mutant cells, an important gene involved in cell growth was suppressed. This indicates that the miRNAs exported from cancerous cells can influence gene expression in neighboring cells. Getting rid of such cancer-suppressing miRNAs could give cancer cells a growth advantage over normal cells to promote tumor growth. Cha, Franklin et al. also suggest that it might be possible to create a non-invasive test to detect colorectal cancer by monitoring the levels of circulating miRNAs in patients. Potential treatments for the disease could also target these miRNAs. https://doi.org/10.7554/eLife.07197.002 Introduction An emerging paradigm in the study of cell signaling is the potential role for post-transcriptional gene regulation by extracellular RNAs. microRNAs (miRNAs) are perhaps the best characterized class of small noncoding RNAs (ncRNAs) that have been detected in extracellular fluids (Valadi et al., 2007). Mature miRNAs are 21–23 nucleotides in length and bind to target mRNAs to inhibit their expression (Krol et al., 2010). Because miRNAs imperfectly pair with their mRNA targets, they can potentially regulate hundreds of transcripts within a genome (Bartel and Chen, 2004). However, individual miRNAs exhibit exquisite tissue-specific patterns of expression (Wienholds et al., 2005), control cell fate decisions (Alvarez-Garcia and Miska, 2005), and are often aberrantly expressed in human cancers (Thomson et al., 2006), affording possible disease-specific signatures with diagnostic, prognostic, and therapeutic potential (Lu et al., 2005; Volinia et al., 2006). In addition to their intracellular roles, recent experiments have identified miRNAs outside the cell in extracellular vesicles (EVs) including exosomes or larger vesicles (Valadi et al., 2007; Crescitelli et al., 2013), in high-density lipoprotein particles (Vickers et al., 2011), or in smaller complexes with Argonaute 2 protein (Arroyo et al., 2011). Exosomes are small 40–130 nm vesicles of endosomal origin that are secreted by all cells and can fuse and be internalized by recipient cells (Valadi et al., 2007; Kosaka et al., 2010; Higginbotham et al., 2011; Mittelbrunn et al., 2011; Montecalvo et al., 2012). It has been suggested that protein cargo transfer by exosomes between cells is associated with tumor aggressiveness and metastasis (Skog et al., 2008; Higginbotham et al., 2011; Luga et al., 2012; Hoshino et al., 2013; Costa-Silva et al., 2015). With the discovery that miRNAs and other RNAs can also be packaged into EVs, or exported by other extracellular mechanisms, it remains unclear the extent to which RNA cargo is sorted for export and how it is dysregulated in disease conditions, such as cancer. Despite accumulating evidence that exosomes are biologically active, little is known regarding how oncogenic signaling affects the repertoire of miRNAs or proteins that are selected for secretion. Given the potential of cancer-derived secreted RNAs to modulate the tumor microenvironment, elucidation of the potential mechanisms for selective sorting of cargo into exosomes is critical to understanding extracellular signaling by RNA. KRAS mutations occur in approximately 34–45% of colon cancers (Wong and Cunningham, 2008). We have previously shown that exosomes from mutant KRAS colorectal cancer (CRC) cells can be transferred to wild-type cells to induce cell growth and migration (Higginbotham et al., 2011; Demory Beckler et al., 2013). Compared to exosomes derived from isogenically matched wild-type cells, exosomes derived from mutant KRAS cells contain dramatically different protein cargo (Demory Beckler et al., 2013). Here, we show that KRAS status also prominently affects the miRNA profile in cells and their corresponding exosomes. Exosomal miRNA profiles are distinct from cellular miRNA patterns, and exosomal miRNA profiles are better predictors of KRAS status than cellular miRNA profiles. Furthermore, we show that cellular trafficking of miRNAs is sensitive to neutral sphingomyelinase (nSMase) inhibition in mutant, but not wild type, KRAS cells and that transfer of miRNAs between cells can functionally alter gene expression in recipient cells. Results Small ncRNAs are differentially distributed in exosomes Because small RNAs are thought to be sorted at endosomal membranes and since KRAS signaling can also occur on late endosomes (Lu et al., 2005), we hypothesized that oncogenic KRAS signaling could alter RNA export into exosomes. We prepared small RNA libraries from both exosomes and whole cells using isogenically matched CRC cell lines that differ only in KRAS status (Figure 1—source data 1) (Shirasawa et al., 1993). Exosomes were purified using differential centrifugation and consisted of vesicles ranging in size from 40 to 130 nm (Higginbotham et al., 2011; Demory Beckler et al., 2013). These preparations exclude larger microvesicles but contain smaller lipoproteins and probably other small RNA–protein complexes (unpublished observation). Comprehensive sequencing analyses of both cellular and exosomal small RNAs from all three cell lines revealed that more than 85% of the reads from the cellular RNA libraries mapped to the genome, compared to only 50–71% from the exosomal libraries (Figure 1A). The non-mappable reads consisted largely of sequences that contain mismatches to genomic sequences. Figure 1 with 1 supplement see all Download asset Open asset Small RNA sequencing analysis of cellular and exosomal RNAs from CRC cell lines. Shown are (A) total read numbers (y-axis) and the total percentage of mappable reads (red), percentage of unique mappable reads (green), reads that map to multiple genomic locations (dark blue), and those that could not be mapped (cyan). (B) The majority of mappable small RNA reads were derived from noncoding RNAs in cells and exosomes. In cells, the majority of small RNA reads mapped to microRNAs (miRNAs) (miRbase 19), whereas in exosomes, the majority of small RNA reads mapped to repetitive elements. (C) The origin of repetitive reads from exosomal small RNA sequencing is shown. Repeat reads were annotated based on RepeatMasker and Rfam classified into tRNAs, rRNAs, snRNAs, and others. (D) The length distribution of reads mapping to miRNA hairpins was determined for small RNA reads from the three CRC cell lines and their purified exosomes. Colors represent the nucleotides identified for the 5′ base, T (cyan), A (red), G (green), and C (blue). Figure 1—figure supplement 1. https://doi.org/10.7554/eLife.07197.003 Figure 1—source data 1 Colorectal cancer (CRC) cell lines. Small RNA sequencing libraries were prepared from three isogenic CRC cell lines with the indicated alleles of KRAS. Table is based on work done in Demory Beckler et al. (2013). https://doi.org/10.7554/eLife.07197.004 Download elife-07197-fig1-data1-v2.docx The global small RNA profiles identified reads from various classes of RNA, including miRNAs, with differential enrichment of specific RNAs in both the cellular and exosomal fractions. Compared to cellular RNA samples, which displayed an enrichment of miRNA sequences (∼70%), miRNA reads in exosomal samples comprised a smaller percentage of the total small reads (5–18%) compared to other ncRNA classes (e.g., tRNAs, rRNAs, snRNAs) (Figure 1B,C, Supplementary file 1). Most of these reads appear to be the fragments of larger RNAs, both cytoplasmic and nuclear. It is not clear how these RNAs are associated with and/or deposited into exosomes. The size distribution of cellular small RNA matched that expected from miRNA-derived reads (21–23 nucleotides). However, the small RNA read distribution from exosomes was much broader with many reads smaller than 22 nucleotides in length (Figure 1—figure supplement 1). Given that these reads map to RNAs other than known miRNAs, these data suggest that a large proportion of small exosomal RNA reads is derived from processing of other RNAs, in addition to post-transcriptionally modified miRNA reads that are apparently subject to editing, trimming, and/or tailing (Koppers-Lalic et al., 2014). Consistent with this, when read identity was restricted to miRNAs by mapping back to known miRNA hairpin sequences, the length distribution of mappable reads was nearly identical between cells and exosomes (Figure 1D). miRNAs are differentially enriched in exosomes dependent on KRAS status Focusing on mappable reads, we sought to ascertain whether miRNAs might be differentially represented when comparing cells to their secreted exosomes. For this, we quantified the relative abundance of individual miRNAs and made pairwise comparisons between normalized miRNA reads. Spearman correlation analyses demonstrated high correlation between replicates of individual cell lines (r = 0.95–0.96) and between cellular data sets differing only in KRAS status (r = 0.92–0.96) (Figure 2—figure supplements 1–3). In contrast, the miRNA profiles from exosomes compared to their parent cells were less correlated (DKO-1 r = 0.67–0.81, DKs-8 r = 0.64–0.71, DLD-1 r = 0.64–0.69) (Figure 2—figure supplements 1, 2, 4). We next utilized principal component (PC) analysis to determine whether the overall miRNA profiles could distinguish between cells and exosomes and/or between wild-type and mutant KRAS status. The miRNA profiles from the three cell lines all clustered close to one another indicating that overall miRNA expression profiles are fairly similar among the different cell types (Figure 2). In marked contrast, PC analysis revealed that exosomal miRNA profiles clearly segregate according to KRAS status (Figure 2). Relatively, minor differences between cellular miRNA expression profiles become much more prominent when comparing exosomal miRNA patterns (Figure 2—figure supplement 3). This indicates that the presence of a mutant KRAS allele alters sorting of specific miRNAs to exosomes, a finding that has potentially important implications for biomarker development. Figure 2 with 4 supplements see all Download asset Open asset Small RNA composition segregates with KRAS status. Principal component analysis was performed comparing small RNA sequencing data sets from CRC cells and exosomes. The small RNA composition from cells differed significantly from exosomes. Nevertheless, clustering showed that mutant KRAS status could be inferred from small RNA composition. Also see Figure 2—figure supplements 1–4. https://doi.org/10.7554/eLife.07197.006 To gain more insight into the relative abundance of miRNAs in cells vs matched exosomes, we examined the most abundant miRNA species in the various sequencing libraries (determined by mean reads of individual miRNAs). For many miRNAs, exosomal abundance correlated with cellular abundance (Supplementary file 2). However, calculation of fold changes among the three isogenic KRAS cell lines, and exosomes released from these cells, showed that distinct subsets of miRNAs are enriched in either exosomes or cells (Tables 1, 2). For all three cell lines, 25 miRNAs were consistently upregulated in cells and 29 miRNAs were consistently upregulated in exosomes (Figure 3A,B). Additionally, the diversity of miRNAs was substantially greater in mutant KRAS DKO-1 exosomes (94 unique miRNAs) compared to parental DLD-1 or wild-type KRAS DKs-8 exosomes (Figure 3B). A select subset of cell and exosomally targeted miRNAs were validated separately by quantitative reverse-transcription PCR (qRT/PCR) (Figure 3C). Collectively, these data indicate that the miRNA profiles observed in exosomes are distinct from their parental cells with specific miRNAs preferentially overrepresented or underrepresented in exosomes. We observed a mutant KRAS-specific pattern of secreted miRNAs, consistent with the hypothesis that dysregulation of miRNA metabolism is associated with tumorigenesis, a previously unrecognized feature of mutant KRAS. Table 1 Differential expression of miRNAs in colorectal cancer cells* https://doi.org/10.7554/eLife.07197.011 DKO-1 hsa-miR-548uhsa-miR-16-1-3phsa-miR-33a-3phsa-miR-33a-5p hsa-miR-31-5phsa-miR-181b-3phsa-miR-450a-5phsa-miR-424-5p hsa-miR-9-5phsa-miR-219-5phsa-miR-190ahsa-miR-573 hsa-miR-30d-3phsa-miR-204-5phsa-miR-1226-3phsa-miR-499a-5p hsa-miR-450b-5phsa-miR-499b-3phsa-miR-3662hsa-miR-20a-3p hsa-miR-27b-5phsa-miR-5701hsa-miR-4677-3phsa-let-7i-5p hsa-miR-331-3phsa-miR-31-3phsa-miR-651hsa-miR-1306-5p hsa-miR-147bhsa-miR-3611hsa-miR-1305hsa-miR-148a-3p hsa-miR-27b-3phsa-miR-1306-3phsa-miR-374b-3phsa-miR-1260b hsa-miR-3940-3phsa-miR-200c-5phsa-miR-548ar-3pDKs-8 hsa-miR-132-5phsa-miR-484hsa-miR-374a-5phsa-miR-1180 hsa-miR-1307-3phsa-miR-200a-5phsa-miR-548o-3phsa-miR-149-5p hsa-miR-3615hsa-miR-100-5phsa-miR-197-3phsa-miR-378a-5p hsa-let-7a-3pDLD-1 hsa-miR-141-3phsa-miR-26b-5phsa-miR-24-3phsa-miR-3074-5p hsa-miR-15a-5phsa-miR-27a-3phsa-miR-3613-5phsa-miR-30b-5p hsa-miR-29a-3phsa-miR-301a-5phsa-let-7i-3phsa-miR-185-5p hsa-let-7g-5phsa-miR-23b-3phsa-miR-22-3pDKO-1 and DKs-8 hsa-miR-141-5phsa-miR-582-5pDKO-1 and DLD-1 hsa-miR-556-3phsa-miR-374a-3phsa-miR-106b-5phsa-miR-17-3p hsa-miR-24-1-5phsa-miR-340-3pDLD-1 and DKs-8 hsa-miR-24-2-5phsa-miR-106a-5phsa-miR-30e-5phsa-miR-107 hsa-miR-429hsa-miR-98-5phsa-miR-425-5phsa-miR-140-5p hsa-miR-93-5phsa-miR-210hsa-miR-126-3phsa-miR-194-5p hsa-miR-29b-3phsa-miR-15b-5phsa-miR-362-5phsa-miR-27a-5p hsa-miR-454-3phsa-miR-452-5phsa-miR-196b-5pDKO-1, DLD-1 and DKs-8 hsa-miR-32-5phsa-miR-582-3phsa-miR-542-3phsa-miR-96-5p hsa-miR-101-3phsa-miR-18a-5phsa-miR-3529-3phsa-miR-7-5p hsa-miR-19a-3phsa-miR-142-3phsa-miR-20a-5phsa-miR-32-3p hsa-miR-130b-5phsa-miR-1278hsa-miR-7-1-3phsa-miR-590-3p hsa-miR-4473hsa-miR-17-5phsa-miR-103a-3phsa-miR-103b hsa-miR-19b-3phsa-miR-340-5phsa-miR-200a-3phsa-miR-34a-5p hsa-miR-372 * miRNAs differentially enriched in cells when comparing mean reads in exosomes vs cell. miRNAs expression patterns were compared between DKs-8, DKO-1, and DLD-1 cells. miRNAs were identified that were enriched in just one of the three cell types or that overlapped between combinations of cells. For miRNAs, 25 were identified that are expressed in all three cell types, 13 were enriched in DKs-8 cells, 15 in DLD-1 cells, and 39 were unique to DKO-1 cells. Table 2 Differential distribution of miRNAs in exosomes* https://doi.org/10.7554/eLife.07197.012 DKO-1 hsa-miR-139-5phsa-miR-3178hsa-miR-151bhsa-miR-125b-1-3p hsa-miR-193b-3phsa-miR-935hsa-miR-130b-3phsa-miR-628-3p hsa-miR-139-3phsa-let-7d-3phsa-miR-589-3phsa-miR-4532 hsa-miR-451ahsa-miR-6087hsa-miR-151a-5phsa-miR-940 hsa-miR-222-3phsa-miR-766-5phsa-miR-505-5phsa-miR-3187-3p hsa-miR-125a-3phsa-miR-3679-5phsa-miR-4436b-3phsa-miR-4787-3p hsa-miR-2277-3phsa-miR-361-5phsa-miR-1293hsa-miR-3183 hsa-miR-3162-5phsa-miR-642a-3phsa-miR-642b-5phsa-miR-197-5p hsa-miR-324-3phsa-miR-145-3phsa-miR-3182hsa-miR-3127-3p hsa-miR-3127-5phsa-miR-4728-3phsa-miR-3184-5phsa-miR-125b-5p hsa-miR-186-5phsa-miR-1hsa-miR-100-5phsa-miR-423-3p hsa-miR-766-3phsa-miR-4753-5phsa-miR-145-5phsa-miR-4724-5p hsa-miR-373-3phsa-miR-223-5phsa-miR-1307-5phsa-miR-1914-3p hsa-miR-3121-3phsa-miR-3613-3phsa-miR-205-5phsa-miR-98-3p hsa-miR-23a-3phsa-miR-3124-5phsa-miR-3656hsa-miR-3918 hsa-miR-4449hsa-miR-378chsa-miR-3138hsa-miR-1910 hsa-miR-3174hsa-miR-4466hsa-miR-3679-3phsa-miR-3200-5p hsa-miR-6511b-5phsa-miR-1247-5phsa-miR-22-3phsa-miR-877-5p hsa-miR-4687-3phsa-miR-1292-5phsa-miR-181c-5phsa-miR-6131 hsa-miR-6513-5phsa-miR-3661hsa-miR-132-3phsa-miR-214-3p hsa-miR-574-3phsa-miR-3190-3phsa-miR-326hsa-miR-3191-5p hsa-miR-3198hsa-miR-3928hsa-miR-629-3phsa-miR-4489 hsa-miR-4700-5phsa-miR-5006-5phsa-miR-5088hsa-miR-2110 hsa-miR-3911hsa-miR-3146DKs-8 hsa-miR-1224-5phsa-let-7b-5phsa-miR-155-5phsa-let-7c hsa-let-7a-5phsa-miR-146b-5phsa-miR-4647hsa-miR-4494 hsa-miR-711hsa-miR-1263DLD-1 hsa-miR-1226-5phsa-miR-4745-5phsa-miR-4435hsa-miR-939-5p hsa-miR-409-3phsa-miR-1304-3pDKO-1 and DKs-8 hsa-miR-146a-5phsa-miR-4508hsa-miR-224-5phsa-miR-4429 hsa-miR-222-5phsa-miR-629-5phsa-miR-4492hsa-miR-3653 hsa-miR-320ahsa-miR-1290hsa-miR-1262hsa-miR-5010-5p hsa-miR-204-3phsa-miR-4461hsa-miR-5187-5pDKO-1 and DLD-1 hsa-miR-483-5phsa-miR-4658hsa-miR-4758-5phsa-miR-492 hsa-miR-5001-5phsa-miR-371a-5phsa-miR-1323hsa-miR-371b-3p hsa-miR-501-3phsa-miR-4446-3phsa-miR-6511a-5phsa-miR-30a-3p hsa-miR-4727-3pDLD-1 and DKs-8 hsa-miR-28-3phsa-miR-3934-5pDKO-1, DLD-1 and DKs-8 hsa-miR-658hsa-miR-320dhsa-miR-4792hsa-miR-1246 hsa-miR-320ehsa-miR-4516hsa-miR-320bhsa-miR-4488 hsa-miR-1291hsa-miR-320chsa-miR-4634hsa-miR-3605-5p hsa-miR-4741hsa-miR-3591-3phsa-miR-122-5phsa-miR-486-3p hsa-miR-184hsa-miR-223-3phsa-miR-3651hsa-miR-486-5p hsa-miR-3180hsa-miR-3180-3phsa-miR-3168hsa-miR-4497 hsa-miR-423-5phsa-miR-3184-3phsa-miR-150-5phsa-miR-664a-5p hsa-miR-182-5p * miRNAs differentially enriched in exosomes when comparing mean reads in exosomes vs cell. miRNAs expression patterns were compared between exosomes purified from DKs-8, DKO-1, and DLD-1 cells. miRNAs were identified that were enriched in exosomes from just one of the three cell lines or that overlapped between combinations of cell lines. 29 miRNAs were common between exosomes from all three cell lines. 94 were enriched in exosomes from DKO-1 cells, 10 in exosomes from DKs-8 cells, and only 6 in exosomes from DLD-1 cells. Figure 3 Download asset Open asset KRAS-dependent regulation of miRNAs in exosomes and cells. Differentially distributed miRNAs in (A) cells and (B) exosomes from the three CRC cell lines differing in KRAS status. (C) qRT-PCR validation of selected miRs from DKs-8 and DKO-1 cellular and exosomal RNA samples normalized to U6 snRNA. Fold changes were calculated using the ΔΔC(t) method comparing exosomes to cells. Negative fold changes indicate greater enrichment in cells, and positive fold changes indicate greater enrichment in exosomes. Also see Supplementary file 2. https://doi.org/10.7554/eLife.07197.013 KRAS-dependent sorting of miRNAs miR-100 Down regulation of miR-100-5p was observed in mutant KRAS DKO-1 and parental DLD-1 cells compared to wild-type KRAS DKs-8 cells (Table 1). This is consistent with reports that have shown decreased miR-100 expression in metastatic cancers (Petrelli et al., 2012; Gebeshuber and Martinez, 2013). miR-100 has also been shown to negatively regulate migration, invasion, and the epithelial–mesenchymal transition (EMT) (Chen et al., 2014; Wang et al., 2014; Zhou et al., 2015). Interestingly, miR-100 was enriched in exosomes derived from mutant KRAS cells (>eightfold and >threefold enriched in DKO-1 and DLD-1 exosomes, respectively; Supplementary file 2), suggesting that decrease of miR-100 in cells is due to secretion in exosomes. This is in line with findings that circulating levels of miR-100 are upregulated in the plasma of mutant KRAS-expressing mouse pancreatic cancer models and in patients with pancreatic cancer (LaConti et al., 2011). More broadly, the observation that miR-100-5p specifically accumulates in exosomes suggests that there may be sequence-specific requirements for the sorting of certain miRNAs into exosomes. miR-10b Our RNA sequencing data identified miR-10b as preferentially secreted in exosomes isolated from cells harboring a wild-type KRAS allele (>threefold-change and >twofold-change enrichment in DKs-8 and DLD-1 exosomes, respectively) but retained in mutant KRAS DKO-1 cells (∼threefold-change cell enrichment). miR-10b is referred to as an oncomiR because it is frequently upregulated during progression of various cancers, including CRC (Ma, 2010). miR-320 miR-320 is aberrantly expressed in several types of cancer, including colon cancer. It is expressed in the proliferative compartment of normal colonic crypts (Schepeler et al., 2008; Hsieh et al., 2013). miR-320 members (miR-320b, -c, d, and -e) were abundant in both mutant KRAS (DKO-1) and wild-type KRAS (DKs-8) exosomes, but underrepresented in the matched cells, indicating that some miRNAs are transcribed and predominantly exported into exosomes, independent of KRAS status (Table 2, Supplementary file 1A). Of these family members, miR-320a and miR-320b were the most abundant species represented in exosomes by our RNA sequencing analyses (miR-320a in DKO-1 exosomes, and miR-320b in DKs-8 and DKO-1 exosomes). Interestingly, however, we observed the largest enrichment for miR-320d (fold changes >241 in DKs-8 and >229 in DKO-1) in exosomes relative to cells, despite being ∼fourfold less abundant than miR-320b levels. Because the 3′-terminus may be important in regulating miRNA stability and turnover, coupled to the fact that the sequences of miR-320a-d members differ only at their 3′-termini, enrichment of certain miRNAs in exosomes could be due to higher turnover/decay rates in cells. Exosomal secretion and strand selection Because we observed differential export of specific miRNAs, we investigated whether there might be miRNA sequence-specific sorting signals. Previous reports have shown differential accumulation of 5p or 3p strands in exosomes compared to parental cells (Ji et al., 2014). Thus, we analyzed our data sets to test whether exosomes might be preferentially enriched for one strand over the other. We were able to identify individual miRNAs where the two strands differentially sorted between cells and exosomes. For example, the -5p strands of miR-423 were overrepresented in DKO-1 exosomes but in exosomes from DKs-8 cells, both strands were overrepresented compared to cells (data not shown). This indicates that KRAS status may differentially affect selection of passenger or guide strands for sorting into exosomes for select individual miRNAs. Individual miRNAs often exist as populations of variants (isomiRs) that differ in length and/or nucleotide composition generated by template- or non-template-directed variation (Burroughs et al., 2010; Newman et al., 2011; Neilsen et al., 2012). When we analyzed our sequencing data sets, we did not detect differential accumulation of isomers with variable 5′ termini (data not shown). For cellular miRNAs, most reads were full length with a slight enrichment in 3′ non-templated addition of A-tailed miRNAs, regardless of KRAS status (Figure 4; Figure 4—figure supplement 1). For exosomes, we observed a slight enrichment for C residues added to the 3′ ends of miRNAs from wild-type KRAS cells (Figure 4—figure supplement 1). We did not observe this in mutant KRAS exosomes, where instead, we noticed an increase in 3′ trimming of miRNAs (Figure 4, Figure 4—figure supplement 2). Overall, it remains to be determined whether such modifications constitute a global exosomal sorting signal in these cells. Figure 4 with 3 supplements see all Download asset Open asset Comparison of miRNA 3′ trimming and tailing between cells and exosomes. Data from the heat maps shown in Figure 4—figure supplement 2 were pooled to illustrate overall changes in either 3′ nucleotide additions (tailing) or 3′ resection (trimming) compared to full length miRNA sequences (intact). Overall, the patterns between cells and between exosomes are very similar. A comparison of cells to exosomes shows that exosomes display a slight increase in trimmed miRNAs. Also see Figure 4—figure supplement 2. https://doi.org/10.7554/eLife.07197.014 Consistent with published data, we have shown that miRNA expression patterns vary between parental cells and their cognate exosomes (Tables 1, 2, Figure 3A,B) (Valadi et al., 2007; Mittelbrunn et al., 2011; Ekstrom et al., 2012; Montecalvo et al., 2012; Squadrito et al., 2014). Differential export suggests that specific signals must exist to sort distinct miRNAs (Batagov et al., 2011; Villarroya-Beltri et al., 2013). We therefore conducted MEME analysis to attempt to identify sequence motifs that might serve as targeting signals. When we examined all miRNA reads detected in exosomes, we did not find any global enrichment for specific sequences or motifs, including those reported to be bound by hnRNP A2B1 (GGAG or U/CC) (Bolukbasi et al., 2012; Villarroya-Beltri et al., 2013) (Figure 4—figure supplement 3). However, when we analyzed miR-320 because it is preferentially exported to exosomes independent of KRAS status, we were able to identify the GGAG sequence contained within the 3′ end of the mature sequence. Additionally, upon restricting our analysis to reads from the most differentially expressed miRNAs when comparing exosomes to cells, we found a slight enrichment for C residues, possibly alternating C residues in exosomal miRNAs (Figure 4—figure supplement 3). Sphingomyelinase-dependent sorting of miRNAs to exosomes Although little is understood regarding the molecular mechanisms for packaging exosomal miRNAs, recent evidence suggests that the secretion of miRNAs in exosomes is dependent on ceramide via its production by neutral sphingomyelinase 2 (nSMase 2) (Kosaka et al., 2010; Mittelbrunn et al., 2011). Inhibition of de novo ceramide synthesis by treatment with a nSMase inhibitor impaired exosomal miRNA release, apparently due to decreased formation of miRNA-containing exosomes (Kosaka et al., 2010; Mittelbrunn et al., 2011). To test the role of nSMase in miRNA secretion in our system, we treated CRC cells with the nSMase inhibitor, GW4869. We determined the effect of this inhibitor on miR-10b since it is preferentially found in wild-type KRAS DKs-8 exosomes, miR-100 since it is preferentially found in mutant KRAS DKO-1 and DLD-1 exosomes, and miR-320 which sorts into exosomes regardless of KRAS status. For miR-10b, we did not observe significant changes in its cellular levels after treatment with GW4869 in either wild-type KRAS DKs-8 or mutant KRAS DKO-1 cells (Figure 5C). In contrast, inhibition of nSMase caused a ∼threefold increase in intracellular levels of miR-100 in mutant KRAS DKO-1 cells but remained unchanged in wild-type DKs-8 KRAS cells (Figure 5A,B,C). Similarly, miR-320 levels were found to increase (∼2.5 fold) only in GW4869-treated mutant KRAS DKO-1 cells (Figure 5C). These data are most consistent with the hypothesis that impaired ceramide synthesis alters ce