Stem Cell Rev and Rep (2017) 13:226–243 DOI 10.1007/s12015-016-9713-1
Exosome and Microvesicle-Enriched Fractions Isolated
from Mesenchymal Stem Cells by Gradient Separation Showed Different Molecular Signatures and Functions on Renal Tubular Epithelial Cells
Federica Collino1,2 & Margherita Pomatto2 & Stefania Bruno3 & Rafael Soares Lindoso1 & Marta Tapparo2 & Wen Sicheng4 & Peter Quesenberry4 & Giovanni Camussi2
Published online: 9 January 2017
The Author(s) 2017. This article is published with open access at Springerlink.com
Abstract Several studies have suggested that extracellular vesicles (EVs) released from mesenchymal stem cells (MSCs) may mediate MSC paracrine action on kidney regen- eration. This activity has been, at least in part, ascribed to the transfer of proteins/transcription factors and different RNA species. Information on the RNA/protein content of different MSC EV subpopulations and the correlation with their bio- logical activity is currently incomplete. The aim of this study was to evaluate the molecular composition and the functional properties on renal target cells of MSC EV sub-populations separated by gradient floatation. The results demonstrated het- erogeneity in quantity and composition of MSC EVs. Two peaks of diameter were observed (90–110 and 170–190 nm). The distribution of exosomal markers and miRNAs evaluated in the twelve gradient fractions showed an enrichment in frac- tions with a flotation density of 1.08–1.14 g/mL. Based on this observation, we evaluated the biological activity on renal cell proliferation and apoptosis resistance of low (CF1), medium
Electronic supplementary material The online version of this article (doi:10.1007/s12015-016-9713-1) contains supplementary material, which is available to authorized users.
(CF2) and high (CF3) floatation density fractions. EVs de- rived from all fractions, were internalized by renal cells, CF1 and CF2 but not CF3 fraction stimulated signifi- cant cell proliferation. CF2 also inhibited apoptosis on renal tubular cells submitted to ischemia-reperfusion in- jury. Comparative miRNomic and proteomic profiles reveal a cluster of miRNAs and proteins common to all three fractions and an enrichment of selected mole- cules related to renal regeneration in CF2 fraction. In conclu- sion, the CF2 fraction enriched in exosomal markers was the most active on renal tubular cell proliferation and protection from apoptosis.
Keywords Extracellular vesicles . Mesenchymal stem cells . Exosomes . Microvesicles . Acute kidney injury . Kidney regeneration
Extracellular vesicles (EVs) are small extracellular membrane fragments heterogeneous for their origin, dimensions and con- tent. EVs are mainly composed by exosomes, small homoge- neous vesicles (50–150 nm) formed from the endosomal cell compartment, and by microvesicles, with a more heteroge- neous dimension profile, produced by the direct extrusion of the cell plasma membrane. Ratajczak et al.  first demon- strated that embryonic stem cells may release EVs capable to reprogram hematopoietic progenitors. EVs derived from adult stem/progenitor cells have been described to retain features resembling their cells of origin. In particular, EVs derived from mesenchymal stem cells (MSCs) can mediate the para- crine effects of MSCs in different models of tissue regenera- tion [2–5]. Selective molecules such as proteins, RNAs and
Giovanni Camussi firstname.lastname@example.org
Carlos Chagas Filho Institute of Biophysics, Federal University of Rio de Janeiro, Rio de Janeiro, RJ, Brazil
Department of Medical Sciences and 2i3T, University of Torino, Torino, Italy
Department of Molecular Biotechnology and Healthy Science, Molecular Biotechnology Center, University of Torino, Torino, Italy
Division of Hematology/Oncology, Rhode Island Hospital, Brown University, Providence, RI, USA
Stem Cell Rev and Rep (2017) 13:226–243
small non coding RNAs are compartmentalized inside EVs [1, 6–10], and may be transferred to target cells via EVs. Numerous studies reported the heterogeneity in EV popula- tions. Exosomes showed the presence of subsets of cellular proteins, such as the tetraspanin family members (CD63, CD81 and CD9), heat-shock proteins , and endosomal sorting complexes proteins, such as Alix and TSG101 . Moreover, specific subsets of small RNAs have been de- scribed as selectively incorporated in different EV populations [13, 14]. Functional studies have been usually performed using a heterogeneous EV combination or selective sub- fractions [15, 16]. Numerous protocols to isolate EV sub- populations have been described , but to identify a method to separate a pure population still remains challenging . EVs from MSCs (MSC EVs) have been extensively used for the treatment of acute kidney injury (AKI), showing regener- ative properties. These were ascribed, at least in part, to the transfer of RNA species, such as mRNAs and miRNAs [2, 8, 19, 20]. Comprehensive information on the RNA/protein con- tent of different subpopulation of MSC EVs is currently incomplete.
The aim of this study was to analyze the composition and activity of different EV sub-populations released from MSCs, defining the molecular profile associated with their activity on renal tubular epithelial cells. For this purpose, we fractioned by density gradient the conditioned medium (CM) of MSCs and investigated the different EV subpopulations obtained. Moreover, we combined different gradient fractions (CFs), on the bases of differential expression of exosomal markers and density, and then we characterized their biological activity on renal target cells. The signaling pathway (s) associated with the protein and miRNA cargo content of different EV subpop- ulations have also been investigated in relation with their bi- ological activity.
Materials and Methods
Bone marrow MSCs were purchased by Lonza (Basel, Switzerland) and cultured as previously described . Cells were used within the seven passages. MSCs characterization was performed by cytofluorimetric analysis for the expression of the typical mesenchymal markers as described .
Murine tubular epithelial cells (mTEC) were isolated from the kidneys of healthy C57 mice. Cells were cultured and characterized as previously described by Bruno et al. . mTEC were positive for classical epithelial markers such as: cytokeratin, actin, alkaline phosphatase, aminopeptidase A, and megalin, and negative for CD45, von Willebrand factor, desmin and nephrin.
Preparation of Extracellular Vesicles-Containing Conditioned Medium
For the preparation of the conditioned medium (CM), MSCs (passage 3–7) were cultured in the presence of their expansion medium until 80% of confluence. Conditioned medium was obtained from supernatants of 8.95 ± 0.46 × 103/cm2 MSCs maintained in RPMI medium supplemented with 0.1% BSA for 16 h. The viability of MSCs after starvation was about 86 ± 0.5% as detected by the Muse® Count &Viability Assay Kit (CTRL normal medium 88,9% of vitality) (Millipore, MA, USA). Supernatant was first centrifuged at 1500 g for 20 min, to remove debris and apoptotic bodies and then concentrated at 4 °C, approximately 200-fold, using ul- trafiltration units (Amicon Ultra-PL 3, Millipore) with a 3 kDa molecular weight cut-off as previous described . After the concentration CM-containing EVs in 1% dimethyl sulfoxide was kept at −80 °C until use.
Density Gradient Separation of EVs
A discontinuous iodixanol gradient was prepared as described by Tauro et al. . Solutions of 5, 10, 20 and 40% iodixanol were obtained by mixing an OptiPrepTM (60% w/v aqueous iodixanol solution) (Sigma-Aldrich, St. Louis, MO) with the appropriate amounts of homogenization buffer (0.25 M sucrose/ 10 mM EDTA/ 10 mMTris-HCL, pH 7.4). The gra- dient was formed by layering 3 ml of 40%, 3 ml of 20%, 3 ml of 10% and 2.5 ml of 5% solutions in a 13 ml open top polyallomer tube (Beckman Coulter). 500 μl of CM- containing EVs (cCM-EVs) were overlaid onto the top of the gradient  or at the bottom  and centrifuged for 18 h at 100,000 g at 4 °C (SW 40 Ti rotor, Beckman Coulter Optima L-90 K ultracentrifuge, Indianapolis, IN). Preliminary experiments demonstrated the same vesicles distribution using the protocols listened above (not shown) and the loading from the top of the gradient  was used in the following experiments.
Twelve gradient fractions of 1 ml, were collected from the top of the gradient, diluted with 10 ml PBS for washing and centrifuged at 100,000 g for 2 h at 4 °C. The pellets were re- suspended in 100 μl PBS or medium based on the following use. The density of each fraction was measured by weighing a fixed volume . In some experiments, fractions were divid- ed into three groups (combined fractions, CFs), based on their density and on the expression of specific surface markers (Fractions 1–4 low density, CF1; Fractions 5–8, medium den- sity CF2; Fractions 9–12, high density CF3).
To trace EVs by fluorescence microscopy, MSCs were labeled with Vybrant Cell Tracers Dil and Syto-RNA (Life
Stem Cell Rev and Rep (2017) 13:226–243
Technologies, Carlsbad, CA) as previous described . EVs obtained from labeled cells were concentrated and subjected to gradient separation as described above. For EV incorpora- tion, mTEC were plated in 24-well plated and treated with different doses of labeled cCM-EVs (50,000, 150,000, 300,000 or 600,000 EVs/cell) for 24 h. Quantitative analysis of the EV uptake was conducted by FACS. To determine the incorporation of each fraction, mTEC were seeded into 24- well plates and incubated with labeled EVs (150,000 EVs/ cell) from the different combined fractions for 24 h. The up- take of EVs was analyzed by microscope analysis using the ApoTome system (Carl Zeiss, Oberkochen, Germany). Hoechst 33,258 dye (Sigma-Aldrich) was added for nuclear staining.
In Vitro Models
mTEC were seeded at 1500 cells/well into 96-well plates and cultured in serum free low-glucose DMEM in the absence (vehicle, CTR-) or presence of EVs (1 × 107 EVs/ml, 1 × 108 EVs/ml or 1 × 109 EVs/ml). Cells maintained in low-glucose DMEM (Sigma-Aldrich) plus 10% FCS were used as positive control (CTR+). In selected experiments, to mimic the ischemia damage on renal cells, mTEC cultured in serum free DMEM were placed in hypoxic chambers with 1% O2 for 48 h. The re-oxygenation step was conducted for 24 h, in the absence (vehicle, HY/CTR-) or presence of EVs (1 × 107 EVs/ml). DMEM plus 10% FCS (HY/CTR+) or EGF (10 ng/ml, Sigma) (HY/EGF) in the re-oxygenation step were used as controls. For cell proliferation, DNA synthesis was detected as incorporation of 5-bromo-2′-deoxy-uridine (BrdU) into the cellular DNA at 48 h (Roche Applied Science, Mannheim, Germany). Apoptosis/necrosis was mea- sured by MuseTM Caspase-3/7 Kit (Millipore) following the instructions. Percentage for live, apoptotic, and necrotic cells was measured.
Gradient isolated EVs were analyzed by nanoparticle tracking analysis (NTA), using the NanoSight LM10 system (NanoSight Ltd., Amesbury, UK), equipped with a 405 nm laser and with the NTA 2.3 analytic software, to define their dimension and profile. Camera levels were for all the acquisi- tion at 16 and for each sample, five videos of 30 s duration were recorded. Briefly, cCM-EVs or gradient-separated EVs were diluted (1: 1000 and 1: 10, respectively) in 1 ml vesicle- free physiologic solution (Fresenius Kabi, Runcorn, UK). NTA post-acquisition settings were optimized and maintained constant among fractions, and each video was then analyzed to measure EV mean, distribution and concentration.
EVs from the twelve fractions were characterized by cytofluorimetric analysis using the Guava easyCyte Flow
Cytometer (Millipore) with InCyte software. The following FITC or APC conjugated antibodies were used: CD107, CD81 and CD63. FITC or APC mouse non-immune isotypic IgG (Miltenyi Biotec, Bergisch Gladbach, Germany) were used as controls as previously described . Briefly, imme- diately after labelling for 15 min at 4 °C with antibodies, EVs (1.5 × 108 particles) diluted 1 to 3 were acquired.
Western Blot Analysis
For protein analysis, EVs from different fractions and MSCs were lysed at 4 °C for 30 min in RIPA buffer (20 nM Tris-HCl, 150 nM NaCl, 1% deoxycholate, 0.1% SDS, 1% Triton X-100, pH 7.8) supplemented with protease and phosphatase inhibitors cocktail (Sigma-Aldrich). Protein content in EV fractions were quantified by BCA Protein Assay Kit (Pierce, Thermo Fisher Scientific, Waltham, MA). Ten μg of proteins were then separated by 4% to 15% gradient sodium dodecyl sulfate–polyacrylamide gel electrophoresis. The proteins were transferred onto a PVDF membrane by the iBlotTM Dry Blotting System (Life Technology) and then immunoblotted with the following antibodies: CD63 and ANXA2 (Santa Cruz Biotechnology, Santa Cruz CA), HLA1 (Abcam, Cambridge, United Kingdom), CD29 (Thermo Fisher Scientific) and Integrin alpha-5 (Millipore). The protein bands were visual- ized using a ChemiDocTM XRS + (BioRad) with an enhanced chemiluminescence detection kit (ECL) (GE healthcare, Amersham, Buckinghamshire, UK).
Transmission electron microscopy of EVs was performed by loading EVs from different CF onto 200 mesh nickel formvar carbon coated grids (Electron Microscopy Science, Hatfield, PA) for 20 min. EVs were then fixed with a solution contain- ing 2.5% glutaraldehyde and 2% sucrose and after repeated washings in distilled water, samples were negatively stained with NanoVan (Nanoprobes, Yaphank, NK, USA) and exam- ined by Jeol JEM 1010 electron microscope.
Total RNA was isolated from gradient separated EVs using the mirVana RNA isolation kit (Applied Biosystem) according to the manufacturer’s protocol. RNA from all the twelve fractions or from the three CFs was quantified (Nanodrop ND-1000, Wilmington DE) and the small RNA composition of different CFs was assessed by capillary electrophoresis on an Agilent 2100 Bioanalyzer using the small RNAs kit (Agilent Technologies, Inc., Santa Clara, CA).
Stem Cell Rev and Rep (2017) 13:226–243
GA), and 15 μg of EV proteins from all the CFs, were used for RayBio Label-based (L-Series) Human Antibody Array 1000 (RayBiotech) according to the manufacturer instructions. The arrays were performed in duplicates, using a pull of EVs de- rived from two different MSC preparations/array. The array provides detection of 1000 proteins.
Data analysis was conducted after background signal sub- traction and normalization to positive controls (Mean back- ground +4 standard deviations, accuracy ≈ 99%). Comparison of signal intensities among array images was used to define relative differences in expression levels of each protein among the CF fractions. Differential expression analysis was con- ducted using CF2 as reference fraction. Proteins were consid- ered co-expressed when they showed a Fold change (FC) = 0.65 ≤ FC ≤ 1.5 in all sample tested. Differential ex- pression among CFs was considered when FC distribution: <0.65 or >1.5 in all the sample tested.
Pathway and Gene Ontology Analysis of EV Content
For CF2 miRNA target prediction and biological pathway enrichment analysis, the web-based program DIANA- mirPath  was used. The algorithm microT-CDS was cho- sen to predict EV-derived miRNA targets using default thresh- old (microT = 0.8). Only biological pathways showing P value < 0.01 to all known Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were considered as signifi- cantly enriched.
For protein class analysis and pathway classification Panther classification system was used (http://pantherdb. org/). Gene Ontology (GO) analysis was conducted using David functional annotation tool (http://david-d.ncifcrf.gov.). Functional sorting of the proteins differentially expressed between CF2 and the other fractions was done using Funrich analysis tool .
Data were analyzed using the GraphPad Prism 6.0 Demo program. Statistical analyses were conducted using One-way ANOVA with Dunnett’s or Turkey’s post-tests, where appropriated. Statistical significance was established at P < 0.05.
Heterogeneity of EVs Released from MSCs.
To isolate a total EV population from MSCs, we pre-purified their CM as described in Material and Methods. The CM was then subjected to concentration and the number of isolated EVs was counted by NTA. We observed that the mean number
Purified CF-derived EVs isolated from three different MSCs preparations were analyzed for their miRNA content by quan- titative real time (qRT) PCR using the Applied Biosystems TaqManH MicroRNA Assay Human Panel Early Access kit (Life Technologies), able to profile 754 human mature miRNAs by sequential steps of reverse transcription (Megaplex RT Pools; Life Technologies) using an Applied Biosystems 7900H qRT-PCR instrument as previously de- scribed . Briefly, single stranded cDNA was generated from total RNA sample (80 ng) by reverse transcription using a mixture of looped primers (Megaplex RT kit, Life Technologies) following manufacturer’s protocol. The pre- amplification reaction for each sample was performed using a TaqMan® PreAmp Master Mix 2X (Life Technologies) mixed with specific MegaplexTM PreAmp Primers (10X) (Life Technologies). Pre-amplified products were then diluted, loaded in the TaqMan MicroRNA Array and qRT-PCR exper- iments were performed.
Raw Ct values, automatic baseline and threshold were cal- culated using the SDS software version 2.3. Comparison of miRNA expression was conducted using the Expression Suite software (Life Technologies). Fold change (Rq) in miRNA expression among the three fractions was calculated as 2-ΔΔCt using one of the fraction (CF2) as control and normal- izing the data using global normalization . Confirmation of the expression of specific miRNAs in the three fractions was conducted using the miScript SYBR Green PCR Kit (Qiagen, Valencia, CA, USA). Briefly, 50 ng of input RNA were reverse transcribed using the miScript Reverse Transcription Kit and the cDNA was then used to detect and quantify miRNAs of interest. Experiments were run in tripli- cate using 3 ng of cDNA for each reaction as described by the manufacturer’s protocol (Qiagen). The following miRNAs were screened in all the CF: miR-100, miR-21, miR-24, miR-214, miR-34a, miR-127, miR-30c, miR-29a, miR- 125b, miR-10b, let-7c, miR-99a, miR-17 and miR-20a.
Analysis of miR-451 distribution in the twelve fractions was carried on using the miRCURY LNATM Universal RT microRNA PCR kit (Exiqon, A/S, Vedbaek, Denmark). Fifty pg of reverse transcription reaction products were then com- bined with SYBR Green Master Mix (Exiqon) and LNATM PCR primer mix and analyzed as described by the manufac- turer’s protocol. All the qRT-PCR data were normalized using the UniSp6 and the UniSp2 RNA Spike-in templates respec- tively as cDNA synthesis and RNA extraction controls (Exiqon).
Purified CF-derived EVs isolated from different MSC prep- arations were lysed in 2× Cell Lysis Buffer (RayBiotech, Inc.,
of EVs present in the concentrated CM was around 12,750 ± 3187 particles/cell. The cCM-EVs represented a mix population with different diameter, ranging from 50 to 390 nm (Fig. 1a). Two peaks of diameter were defined around 90–110 and 170–190 nm (Fig. 1a) with a mean diameter in the total population of 136 ± 16 nm and a mode of 119 ± 21 nm. The identification of a heterogeneous population prompts us to define its molecular content and activity. For this purpose, we fractionated the cCM-EVs by density gradient separation, isolating twelve different EV fractions. Density of each fraction was measured by weighing a fixed volume (Fig. 1b) . The analysis of EV distribution by NTA, showed an enrichment of vesicles in fractions 1, 4 and 8 (Fig. 1c).
Fig. 1 Characterization of EVs isolated from the conditioned medium (CM) of human MSCs and subjected to OptiPrep gradient separation. a Nanoparticle Tracking Analysis (NTA) profiles of particle size and distribution of EVs present in the concentrated CM of MSCs (cCM- EVs). cCM-EVs represent a extracellular vesicle population containing exosomes and microvesicles. Two peaks around 90–110 and 170–190 nm were detected. b Weight of the twelve different EV fractions generated by gradient separation of cCM-EVs. c NTA analysis of the percentage of EVs distributed in different fractions after gradient separation. Data reported are mean ± SEM of four different experiments
Stem Cell Rev and Rep (2017) 13:226–243 Characterization of EVs in Different Fractions.
By FACS analysis the expression of CD63 endosome-derived exosome marker was mainly observed in fractions 5–8, with a flotation density of 1.08–1.14 g/mL (Fig. 2a). The other tetraspanin family member, CD81 and the lysosomal- associated membrane protein 1, CD107, showed almost the same sinusoidal pattern (Fig. 2a). As shown in Fig. 2b, the biological activity on mTEC is widely distributed in different fractions, except for the densest fractions, but is significantly higher in fraction 6.
Total RNA was extracted from all the twelve fractions and the major concentration of RNA was isolated from fractions 5, 8, 9 and 11 (not shown). Since it has been suggested a selec- tive export of specific miRNAs from different cell compart- ments , we evaluated the distribution of a group of miRNAs that we previously described in MSC-derived EVs . miR-21 was expressed in almost all the EV fractions. However, as shown in Fig. 2c, the Raw Ct of miR-21 indicat- ed its enrichment in the central fractions (Fig. 2c). Similar enrichment pattern was observed for miR-100, −99a and −24 (data not shown). Conversely, miR-451 was detected only in central fractions enriched in the tetraspanin family markers (Fig. 2C) and this distribution was compatible with the exosomal derivation of this miRNA. To exclude the possibil- ity of confounding technical factors on the obtained results, we performed miR-451 quantification adding a synthetic spike-in (UniSp6) during the cDNA synthesis. We obtained the same results, demonstrating the effective compartmental- ization of miR-451 mainly in the CD63 positive fractions (Fig. 2d).
Biological Effect and Characterization of EV Combined Fractions
In order to pair the biological activity with the molecular con- tent of the EV subpopulations, we combined the gradient frac- tions into three groups (Fractions 1–4, low density CF1; Fractions 5–8, medium density CF2; Fractions 9–12, high density CF3), based on their density and on the expression of exosome surface markers and miRNAs. We first analyzed the up-take of labelled cCM-EVs by mTEC. Figure 3a shows a dose dependent up-take of cCM-EVs. We then tested the uptake of EVs from the different fractions using confocal mi- croscopy (Fig. 3b-d). Using EVs labelled with Vybrant Dil and with Syto-RNA, we observed the concomitant presence within the cells of the two dyes, suggesting the EV delivery of RNAs inside mTEC.
To discriminate the effects of the different EV fractions on cell proliferation and apoptosis, in vitro experiments were per- formed. Despite all EV fractions slightly promoted cell prolif- eration on mTEC, only the low density CF1 and medium density CF2 fractions induced a statistically significant
Stem Cell Rev and Rep (2017) 13:226–243
Fig. 2 Characterization of specific exosome markers and activity of the twelve EV fractions isolated by gradient separation. a FACS analysis of the expression of the tetraspanin members, CD63 and CD81 and the lysosomal-associated protein, CD107 in different EV fractions. The exosome markers showed a relevant co-localization in fractions 5–8, characterized by a flotation density of 1.08–1.14 g/mL. b Evaluation of the effects of different EV fractions on mTEC proliferation after 48 h of stimulation, assessed by BrdU up-take (1 × 107 EVs/ml) in respect to control cells (CTR-, DMEM no FCS). Cells cultured in DMEM plus 10% FCS were used as positive control (CTR + =1). Data are expressed as Ratio means ±SEM. ANOVA with Dunnett’s multicomparison test.
proliferation at 1 × 108 and 1 × 109 particles/ml concentration (Fig. 3e). We next evaluated the effects of the combined frac- tions in a hypoxia/reperfusion model that mimics renal tubule damage during the ischemia reperfusion injury (IRI). For this purpose, mTEC were treated with CFs during the re- oxygenation phase of IRI. Only the CF2 medium-density ves- icles together with the EGF treatment during the reperfusion phase were anti-apoptotic (Fig. 3f). On the contrary, CF1 low- density vesicles or CF3 high density EVs were not significant- ly protective on mTEC submitted to IRI.
*P < 0.05 vs CTR-. c Representive qRT-PCR expression of miR-21 and
miR-451 in the twelve fractions. Raw data analysis of miR-21( shows its expression in almost all the EV fractions. On the contrary, miR-
d Representative qRT-PCR analysis showing the relative quantity of miR- 451 in respect to the synthetic spike-in (UniSp6) used as normalizer. Normalized data showed the distribution of miR-451 inside the exosome- enriched EV fractions. Three experiments were conducted with similar results
To define the molecular composition of the EV gradient fractions, accountable for their different biological potential, further characterization studies were conducted. Transmission electron microscopy analysis performed on purified EVs showed their spheroid morphology (Fig. 4a). The size of EVs present in CF1 was more heterogeneous of the CF2 ones. NTA quantification of the particle distribution among the CF demonstrated an enrichment of vesicles in CF1 and CF2 in respect to CF3 fraction (Fig. 4b). As shown in Fig. 4c, NTA measure of the percentage of size distribution indicated that
) shows selective compartmentalization in the central fractions.
Stem Cell Rev and Rep (2017) 13:226–243
Fig. 3 Activity of the combined gradient fractions on tubular cell proliferation and apoptosis. a Quantitative FACS analysis of the up-take of labelled cCM-EVs by mTEC. A dose dependent up-take of cCM-EVs (50,000–600,000 EVs/cell) was observed after incubation for 24 h. b-d Representative micrograph of the internalization of EVs from CF1 (b), CF2 (c) and CF3 (d) by mTEC after 24 h, observed by confocal microscopy. EVs were collected from MSCs double-stained with Syto- RNA (green) and Vybrant Dil (red). Three experiments were performed with similar results. Nuclei were counterstained with Hoechst dye. Original magnification: ×630. (E-F) Evaluation of the effects of the combined fractions on mTEC proliferation (e) and protection from apoptosis (f). e Absorbance ratio of the BrdU up-take by mTEC incubated for 48 h with EVs from different CFs (1 × 107 EVs/ml,
CF1 and CF2 contain smaller EVs (50–150 nm) whereas the CF3 is also enriched in EVs of larger size (>200 nm). The mean diameter of the different populations was respectively: 134.7 ± 19,7 for CF1, 147 ± 27,5 for CF2 and 169 ± 39 for CF3 (Fig. 4c, lower panel).
The expression of some classical vesicular markers (CD63, ANXA2 and CD29) was compared among the three CFs and their cells of origin by Western blot analysis (Fig. 4d-e). The integrin alpha-5 and HLA-1 were also screened in the different fractions (Fig. 4e). CD63 was almost absent in low-density vesicles of CF1, being enriched in the medium-density EVs of CF2 fraction and into a lesser extent in the high-density CF3 fraction (Fig. 4d). The integrin β1 (CD29), showed the same distribution of the tetraspanin CD63, being enriched in CF2 and CF3 and expressed at low levels in CF1 fraction (Fig. 4e). ANXA2, involved in recruitment of miRNAs in EVs, resulted mainly present in fraction CF2 and CF3 (Fig. 4e), supporting the relevant enrichment of miRNAs detected in these two fractions (Fig. 4f). The HLA-class I and α5-integrin (ITGA5) were selectively expressed by the medium-density fraction CF2 (exosome enriched fraction) and almost absent in
1 × 108 EVs/ml or 1 × 109 EVs/ml) in respect to control cells (CTR-, DMEM no FCS). Cells cultured in DMEM plus 10% FCS were used as positive control (CTR+). f Cell death analysis on mTEC subjected to hypoxia/reperfusion was measured by MuseTM Caspase-3/7 Kit. Cells were subjected to hypoxia for 48 h, then EVs from CFs were added during the 24 h of reperfusion (1 × 107 EVs/ml). Cells maintained in the absence of serum were used as negative control (HY/CTR-). Cells cultured in DMEM plus 10% FCS in the reperfusion phase were used as positive control (HY/CTR+). Black bars indicate apoptosis and grey bars represent necrosis. Data are expressed as means ±SEM. ANOVA with Dunnett’s multicomparison test. *P < 0.05, **P < 0.01, ***P < 0.001 vs CTR- or HYP/CTR-, respectively
CF1 low-density fraction (Fig. 4e). Low expression of these two markers was observed in the CF3 high-density fraction (microvesicle enriched fraction).
miRNA Compartmentalization inside CF Fractions
RNA was extracted by different CFs and spectrometrically quantified. The bioanalyzer profile showed the enrichment of RNA of the size of miRNAs in the medium-density CF2 and high-density CF3 fractions (Fig. 4f). No relevant differ- ences in the total RNA isolated in respect to the EV quantity among the three CFs was observed (CF1: 2.62 ± 0.97, CF2: 2.75 ± 0.78, CF3: 2.94 ± 0.88 × 10−9 ng/particle) (Fig. 5a). The expression of 754 human miRNAs was measured by qRT-PCR. Analyzing the normalized Cycle threshold distri- bution (ΔCt, based on Median calculation) of all expressed miRNAs (Ct < 40), the medium-density EVs of CF2 fraction showed a low correlation with both CF3 and CF1 EV populations (Pearson correlation: 0.50 and 0.69, respectively) (Fig. 5b).
Stem Cell Rev and Rep (2017) 13:226–243
Fig. 4 Characterization of morphology, protein surface and miRNA content of combined fraction EVs. a Transmission electron microscopy (original magnification × 75,000) of EVs from different CFs. All the CFs contain vesicles with a typical cup-shaped morphology (scale bar 100 nm). Particles with a different electron density were observed in the CF3 fraction. b Percentage of EVs isolated for each combined fraction detected by NTA. The CF3 resulted the fraction with less EVs in respect to the low-density CF1 and medium-density CF2 fractions. c Percentage of size distribution of EVs in different CF fractions measured by NTA (upper panel). The mean diameter (nm) of each CF populations was also
Using a cut-off ≤35 Ct value in miRNA expression, we found 162 miRNAs present in all the EV fractions (Fig. 5c; Table 1). The analysis showed a selective package of miRNA subsets in medium-density CF2 EV fraction which were unde- tected in the other fractions (Table 2) (Fig. 5c). Among these miRNAs, we observed 19 miRNA* sequences detected only in the CF2 fraction, as previously described for tumor derived exosome population .
By analyzing the fold change distribution of the co- expressed miRNAs among all the fractions (miRNA intersec- tion), we observed a relevant enrichment of miRNAs in CF2 fraction containing exosome enriched-EVs. Among the 162 miRNAs expressed by all the fractions, we found that 112 miRNAs were significantly up-regulated in the CF2 fraction in respect to CF1 (Fig. 5d) and only 12 miRNAs were down- regulated. The same trend was observed between CF2 and CF3 (89 miRNAs up-regulated and 15 down-regulated in CF2 in respect to CF3) (Fig. 5e).
In CF2, we detected 106 miRNA families with an enrichment of the following families: miR-10, miR-17, miR-154, miR-30, miR-27, miR-15, let-7, miR-379, miR-26, miR-34 and miR- 548 families (not shown).
measured (lower panel). Four different gradients were tested in triplicate. (d-e) Representative Western blot analysis of CD63 (d), integrin β1 (CD29), α5-integrin (ITGA5), annexin A2 (ANXA2) and HLA-class I (e) on MSCs and CF derived EVs. Three different experiments were performed with similar results. (f) Representative bioanalyzer profile of small RNAs performed on EVs from the CFs, showing a relevant enrichment of total miRNAs in the medium density CF2 and high density CF3 fractions in respect to CF1. Three different samples tested in triplicate with similar results
We confirmed by qRT-PCR the compartmentalization of miRNAs highly enriched in MSC EVs and/or connected with kidney regeneration in the CF2 fraction (miR-100, −21, −24, −214, −34a, −127, −30c, −29a, −125b, −10b, −let-7c, −99a, −17 and miR-20a) [8, 27–32] (Fig. 5f). Some miRNAs, such as miR-21, −24, −214, and −17 were also enriched in the dense fraction CF3, possiblly being a common signature of different EV sub-fractions. Reduced miRNA content was de- tected in the low-density EV CF1 fraction (Fig. 5f).
Analysis of the pathways over-represented by the predicted targets of enriched/selective CF2 miRNAs was performed with the DIANA mirPath software, as previously described . We detected an enrichment of 14 KEGG biological pro- cesses for the CF2 enriched miRNAs and 11 KEGG processed for the CF2 selective miRNAs (P < 0.01, FDR corrected) (S1 Table and S2 Table). Sixty-nine miRNAs from the two groups showed strong correlation with the same pathways assembled in metabolic, stem cell associated- and migration/ inflammation processes.
The following metabolic-related pathways were over-rep- resented: fatty acid biosynthesis and metabolism, biosynthesis of unsaturated fatty acids, mucin type O-glycan biosynthesis
Stem Cell Rev and Rep (2017) 13:226–243
Fig. 5 Analysis of miRNA compartmentalization inside the CF fractions. a Quantification of the total RNA isolated from the different CFs. Data are expressed as means ±SEM (ng RNA/EV) of four different experiments. No significant differences were observed among the three fractions. (b-e) qRT-PCR profile of 754 mature miRNAs in the CF fractions. b Scatter plot of normalized Cycle threshold distribution (ΔCt) of all expressed miRNAs (Ct < 40) between CF2 and the other fractions. Pearson correlation was calculated (CF2 vs CF3: 0.50 and vs CF1: 0.69). c Venn diagram showing the miRNAs present in all the EV fractions (n = 162). A subset of miRNAs was specific of the medium- density CF2 EVs and undetected in the others. Fold change distribution of
and lysine degradation (Table 3). Pathways related to stem cells were: Hippo signaling, Wnt signaling and pluripotent stem cells regulated-pathways (Table 3). Migration- inflammation related processes such as ECM-receptor interac- tion, TGF-β signaling pathway, glioma and proteoglycans in cancer were also over-represented (Table 3). Heatmap repre- sentation of the most significantly enriched pathways poten- tially modulated by the 69 miRNAs is shown in Fig. 6.
Protein Characterization in CF Fractions.
Protein composition of different CFs could also define their origin and activity. For this reason, proteins were isolated from
the co-expressed miRNAs (miRNA intersection) between CF1 and CF2 fractions (d) and between CF3 and CF2 fractions (e), showing a general enrichment of miRNAs in fraction CF2. (f) Validation of the different compartmentalization of specific miRNAs in CFs. The expression of miRNAs enriched in MSC EVs and/or connected with kidney regenera- tion was analyzed by qRT-PCR. All the miRNAs tested were reduced in CF1, demonstrating less ability of this fraction to compartmentalize miRNAs. The relative quantity of each miRNA (RQ) was measured using the synthetic spike-in (UniSp2) as normalizer. Three different samples tested in triplicate with similar results. Data are expressed as means ±SEM. ANOVA with Dunnett’s multicomparison test. *P < 0.05 vs CF2
the three CFs and quantitative proteomic analysis of 1000 proteins was then conducted. Raw data was reported as sup- plementary information (S3 Table). We observed first a differ- ent distribution of total proteins among the three fractions with an enrichment of isolated proteins in respect to the EV quan- tity in CF3 high-density fraction (Fig. 7a). Less amount of proteins was isolated in CF1 in respect to the high quantity of EVs detected in this fraction.
Proteomic profile of all CFs showed the presence of 655 proteins in the CF2 fraction over the 1000 analyzed. 435 pro- teins were compartmentalized in the CF3 high-density frac- tion and 581 in the CF1 low-density fraction. A group of 413 proteins was detected in all the EV populations. Panther
Stem Cell Rev and Rep (2017) 13:226–243 235
List of the miRNAs detected in all the CF fractions. The miRNAs confirmed by qRT-PCR were labeled in red
List of miRNAs present in all the CFs
hsa-miR-1243 hsa-miR-20a hsa-miR-15b
hsa-miR-15a hsa-miR-302c hsa-miR-720 hsa-miR-99a hsa-miR-19b hsa-miR-367 hsa-miR-21 hsa-miR-574-3p hsa-miR-197 hsa-miR-16 hsa-miR-409-3p
hsa-miR-886-5p hsa-miR-92a hsa-miR-30e-3p hsa-miR-152 hsa-miR-31# hsa-miR-335 hsa-miR-127 hsa-miR-331 hsa-miR-181a hsa-miR-1290 hsa-miR-125a-5p
hsa-miR-27b hsa-miR-218 hsa-miR-34a hsa-miR-625# hsa-miR-539 hsa-miR-146a hsa-miR-130a hsa-miR-411 hsa-miR-126 mmu-miR-374-5p hsa-miR-590-5p
hsa-miR-532-3p hsa-miR-214# hsa-miR-320B hsa-miR-125b-1# hsa-miR-629 hsa-miR-571 hsa-miR-22# hsa-miR-139-5p hsa-miR-302a hsa-miR-136# hsa-miR-27a# hsa-miR-655 hsa-miR-455 hsa-miR-638 hsa-miR-204 hsa-miR-193a-3p hsa-miR-628-3p hsa-miR-202 hsa-miR-572
hsa-miR-1274B hsa-miR-376c hsa-miR-25
hsa-miR-1274A hsa-miR-106b hsa-miR-766
hsa-miR-24 hsa-miR-26a hsa-miR-532
hsa-miR-100 hsa-miR-1260 hsa-let-7g
hsa-miR-125b hsa-miR-382 hsa-miR-133a
hsa-miR-191 hsa-miR-365 hsa-miR-345
hsa-miR-214 hsa-miR-132 hsa-miR-195
hsa-miR-193b hsa-miR-328 hsa-miR-148a
hsa-miR-29a hsa-miR-222# hsa-miR-483-5p
hsa-miR-31 hsa-miR-224 hsa-miR-30d
hsa-miR-30c hsa-miR-28-3p hsa-miR-26b
hsa-let-7e mmu-miR-140 hsa-miR-1271
hsa-miR-145 rno-miR-7# hsa-miR-660
hsa-miR-320 hsa-miR-484 hsa-miR-30b hsa-miR-10b hsa-miR-422a hsa-miR-657 hsa-miR-99b hsa-miR-30a-3p
hsa-miR-27a hsa-miR-212 hsa-miR-143 mmu-miR-93 hsa-miR-486 hsa-miR-10b# hsa-miR-146b hsa-miR-376a
hsa-miR-34b hsa-miR-196b hsa-miR-339-3p hsa-miR-335# hsa-miR-296 mmu-miR-137 hsa-miR-130b hsa-miR-301
hsa-miR-17 hsa-miR-374 hsa-miR-20b
hsa-miR-221 hsa-miR-10a hsa-miR-93#
hsa-miR-199a-3p hsa-miR-425-5p hsa-miR-432
hsa-miR-106a hsa-miR-34a# hsa-miR-494
hsa-miR-138 hsa-miR-150 hsa-miR-193b#
hsa-let-7b mmu-miR-134 hsa-miR-28
hsa-miR-886-3p hsa-miR-193a-5p hsa-miR-744
hsa-miR-30a-5p hsa-let-7c mmu-miR-495
Stem Cell Rev and Rep (2017) 13:226–243
List of the miRNAs selectively detected inside the CF2 Selective miRNA in medium-density CF2 EVs
both fractions. Selective proteins were up-regulated in CF1 and CF3 in respect to CF2 fractions (FC CF2 < 0.65) (CF1: TGFBR2 and SFRP1 and CF3: NTRK2, EPHA5, SRMS, LTK, TNFRSF1B, TGFBR1).
To investigate the differences at biological levels among fractions, we performed functional enrichment analyses of the proteins upregulated in CF2 in respect to CF3 and CF1. Enriched GO molecular function (GO-MF) terms over- represented by the proteins upregulated in CF2 in respect to CF3, were mainly incorporated in: transmembrane receptor and receptor-binding activities, G-protein couple receptor ac- tivity, cytokine and hormone activities and cell adhesion mol- ecules (Fig. 8a, outer chart). The GO-MF overrepresented for proteins upregulated in CF2 in respect to the low-density frac- tion CF1 were largely associated with cell adhesion mole- cules, receptor binding and activity (Fig. 8a, inner chart). Selective molecular functions were metallopeptidases and protein tyrosine-kinase activity (Fig. 8a, inner chart).
Enriched GO biological processes (GO-BP) of proteins upregulated in CF2 in respect to CF3 were chemotaxis, signal transduction via receptor interaction, cell proliferation, re- sponse to wounding, cell-cell signaling, development and in- flammation (Fig. 8b). Enriched GO-BP were also observed, for proteins down-regulated in CF1 in respect to CF2. These include common processes such as signal transduction via receptor interaction, cell proliferation, response to wounding and cell-cell signaling. Selective processes such as cell adhe- sion, regulation of cell differentiation and apoptosis were also detected (Fig. 8c).
The results of the present study demonstrated that EVs isolat- ed from the conditioned medium of MSCs, using the gradient separation technique, are heterogeneous in their quantity and composition. Based on the differential expression of specific exosome enriched-markers and different density among the twelve fractions, we grouped them in three combined frac- tions. The combined fractions, displayed differential pro- proliferative and anti-apoptotic activities on renal tubular ep- ithelial cells. Comparative miRNome and proteomic profiles, revealed a cluster of miRNAs and proteins common to the three vesicle fractions and fraction specific subsets of RNAs and proteins.
EVs have been described as important players of the MSC secretome. In fact, the beneficial effect of MSC treatment on different models of acute and chronic damages was mimic by a vesicles-based therapy [2, 16, 34–36]. EVs are composed of several subpopulations and most of their characteristics are not exclusive of a specific class of EVs . For this reason, the possibility that they act as distinct biological entities is now under evaluation .
hsa-miR-107 hsa-miR-122 hsa-miR-135b hsa-miR-190 hsa-miR-198 hsa-miR-199b hsa-miR-205 hsa-miR-211 hsa-miR-301b hsa-miR-302b hsa-miR-32 hsa-miR-331-5p hsa-miR-34c hsa-miR-369-3p hsa-miR-375 hsa-miR-381 hsa-miR-383 hsa-miR-449 hsa-miR-450a hsa-miR-492 hsa-miR-500 hsa-miR-502 hsa-miR-518a-3p hsa-miR-542-3p
hsa-miR-542-5p hsa-miR-545 hsa-miR-570 hsa-miR-579 hsa-miR-616 hsa-miR-95 hsa-miR-98 dme-miR-7 hsa-let-7f-2# hsa-let-7 g# hsa-miR-106b# hsa-miR-1179 hsa-miR-1201 hsa-miR-1226# hsa-miR-1227 hsa-miR-1254 hsa-miR-1270 hsa-miR-1276 hsa-miR-1291 hsa-miR-130b# hsa-miR-144# hsa-miR-15a# hsa-miR-15b# hsa-miR-16-1#
hsa-miR-191# hsa-miR-206 hsa-miR-20a# hsa-miR-223# hsa-miR-24-1# hsa-miR-24-2# hsa-miR-25# hsa-miR-27b# hsa-miR-302d hsa-miR-361-3p hsa-miR-378 hsa-miR-380-5p hsa-miR-411# hsa-miR-550 hsa-miR-580 hsa-miR-603 hsa-miR-744# hsa-miR-770-5p hsa-miR-941 hsa-miR-942 hsa-miR-944 hsa-miR-99a#
pathway analysis on the CF proteome (413 proteins) found high representation of the following pathways: interleukin mediated signaling, inflammation mediated by chemokine and cytokine, TGF-β signaling, gonadotropin release hor- mone receptor, angiogenesis and Wnt signaling pathways (Fig. 7b). CFs contained proteins which originate from differ- ent cellular compartments such as extracellular region, plasma membrane, extracellular matrix, vesicle lumen and secretory granules (Table 4). Moreover, proteins of the CF proteome were mainly incorporated in the following classes: signaling molecules, receptors, defense/immunity proteins, enzyme modulators and cell adhesion molecules (Fig. 7c). The top ten GO biological processes overrepresented by the CF prote- ome were mainly associated with regulation of cell prolifera- tion, response to wounding, enzyme linked receptor protein signaling pathway, inflammatory response, receptor linked signal transduction, immune response and regulation of phos- phorylation (Fig. 7d).
Differential expression analysis reveals that of the 655 pro- teins expressed in the medium-density EVs from CF2 fraction, 75 proteins were significantly up-regulated in respect to CF3 and 92 proteins resulted up-regulated in respect to CF1 (Table 5). Of them, 27 were commonly down-regulated in
Stem Cell Rev and Rep (2017) 13:226–243 Table3 Biologicalpathways
4 5 12 7 40 10 24 16 89 18 48 19 59 22
100 24 128 30 9 8 17 12 43 5 2 2 59 10 9 10 163 10 19 6 58 5 54 3 40 8 24 4 74 3 38 6 61 10 26 6 16 3
over-represented by the 69 enriched/selected miRNAs of the medium-density EVs in CF2 fraction (P < 0,01; corrected FDR)
Fatty acid biosynthesis
Fatty acid metabolism
Hippo signaling pathway
TGF-beta signaling pathway
Signaling pathways regulating pluripotency of stem cells Proteoglycans in cancer
Biosynthesis of unsaturated fatty acids Mucin type O-Glycan biosynthesis Wnt signaling pathway
Thyroid hormone signaling pathway
Glycosaminoglycan biosynthesis – heparan sulfate / heparin Pathways in cancer
Estrogen signaling pathway
FoxO signaling pathway
PI3K-Akt signaling pathway
<1e-325 <1e-325 <1e-325 <1e-325 <1e-325 <1e-325 <1e-325 <1e-325 <1e-325 8.88E-16 8.65E-11 8.54E-10 1.99E-09 1.56E-08 3.91E-07 9.90E-07 2.44E-06 4.67E-06 8.59E-06 1.84E-05 1.90E-05 0.0004 0.0005 0.0015 0.0026 0.0031
EVs released by MSCs, were a mixed population with dif- ferent diameters and with two prominent peaks around 100 and 180 nm, as detected by NTA. EVs isolated by discontinuous density gradient separation showed a heterogeneity in quantity
and expression of the classical exosomal markers. The exosomal markers were mainly detected in the medium- density gradient fractions (1.08–1.14 g/mL density) which co- expressed the CD63 with classical mesenchymal vesicular
Fig. 6 Heatmap representation of the most significantly enriched pathways potentially modulated by the selective/enriched miRNAs in CF2 fraction. Enrichment analysis of the pathways over-represented by the predicted targets of enriched/selective CF2 miRNAs was conducted using the software DIANA mirPath. miRNAs with similar patterns in
targeting significant pathways clustered together (69 miRNAs). Strong correlation with the pathways assembled in metabolic, stem cell associated- and migration/inflammation processes was observed. Only pathways targeted by the selected miRNAs with a P-value < 0,01 (FDR corrected) were considered for the analysis
Stem Cell Rev and Rep (2017) 13:226–243
Fig. 7 Analysis of CF proteome. The protein composition of the CFs defines their origin and potential activity. (a) Quantification of the total proteins isolated from EVs derived from different CFs, demonstrating an enrichment of isolated proteins/EV in CF3 in respect to the other CFs. Data are expressed as means ±SEM (ng protein/EV) of four different experiments. ANOVA with Turkey’s multicomparison test. *P < 0.05 vs CF1. (b) Panther pathway analysis on the CF proteome shows abundance of interleukin, inflammation, TGF-β, gonadotropin release hormone receptor, angiogenesis and Wnt signaling pathways. (c) Distribution of markers such as CD29 and ITGA5 , and ANXA2 an impor- tant effector of miRNA recruitment in EVs . The selective expression of miR-451, in the central frac- tions of the gradient EVs, further supports the enrichment in exosomes of this fraction, since previous studies in tumor cell lines showed that this miRNA is preferentially sorted into exosomes . According to Xu et al. , the denser fractions contained microvesicle-enriched populations. CF3 denser fraction, de- spite the presence of CD63, does not or barely express the HLA-1 and ITGA5 suggesting that is a different population from CF2. Moreover, CF3 is enriched in EVs of larger size (>200 nm) and showed the highest mean diameter among the
the 413 proteins compartmentalized in all the CFs in different protein classes. The CF proteome contained the following class of proteins: signaling molecules, receptors, defense/immunity proteins, enzymes and cell adhesion molecules. d David GO-BB overrepresented by the CF proteome (P < 0,001; FDR 1%). The top ten processes were associated with regulation of cell proliferation, response to wounding, enzyme linked receptor protein signaling pathway, inflammatory response, receptor linked signal transduction, immune response, and regulation of phosphorylation
different EV populations, supporting the enrichment in microvesicles of this fraction. The low-density CF1 fraction contained an undefined population of EVs, that expressed low levels of CD29 and ANXA2 and were negative for all the other markers tested.
Xu et al. demonstrated that exosomes and microvesicles isolated by differential centrifugation from human colon can- cer cells, have distinct biological activities, being able to pro- mote invasiveness at different rating . Moreover, Aliotta et al. showed that exosomes-enriched and microvesicles- enriched populations from mouse MSCs, have different ef- fects when infused into mice, in a model of monocrotaline induced pulmonary hypertension . The same biological
Stem Cell Rev and Rep (2017) 13:226–243 Table4 DAVIDanalysisto
5.26E-91 1.45E-86 5.40E-86 7.31E-15 1.84E-14 4.34E-12 4.95E-08 1.01E-06 8.08E-06 1.13E-04 5.17E-04
8.07E-04 0.002 0.003 0.003 0.005 0.006
determine GO term enrichment (GOTERM_CC_FAT) for genes coding for proteins of the CF proteome (n = 413) (P < 0,01; FDR 1%)
GO:0005615 ~ extracellular space GO:0044421 ~ extracellular region part GO:0005576 ~ extracellular region GO:0031226 ~ intrinsic to plasma membrane GO:0005887 ~ integral to plasma membrane GO:0009986 ~ cell surface
GO:0044459 ~ plasma membrane part GO:0009897 ~ external side of plasma membrane GO:0043235 ~ receptor complex
GO:0031012 ~ extracellular matrix
GO:0060205 ~ cytoplasmic membrane-bounded vesicle lumen
GO:0031983 ~ vesicle lumen
GO:0031093 ~ platelet alpha granule lumen GO:0005578 ~ proteinaceous extracellular matrix GO:0005886 ~ plasma membrane
GO:0031091 ~ platelet alpha granule GO:0030141 ~ secretory granule
150 38.66 167 43.04 226 58.25
92 23.71 90 23.20 43 11.08
116 29.90 24 6.19 19 4.90 31 7.99 11 2.84
11 2.84 10 2.58 27 6.96
152 39.18 11 2.84 19 4.90
4.03 E-94 1.11E-89 4.14E-89 5.60E-18 1.41E-17 3.29E-15 3.79E-11 7.76E-10 6.19E-09 8.64E-08 3.96E-07
6.18E-07 2.09E-06 2.32E-06 2.62E-06 4.17E-06 5.13E-06
1.04E-91 2.85E-87 1.06E-86 1.44E-15 3.63E-15 8.56E-13 9.75E-09 2.00E-07 1.59E-06 2.22E-05 1.02E-04
1.59E-04 5.37E-04 5.97E-04 6.74E-04 0.001 0.001
differences in mouse and human MSC EVs sub-populations were demonstrated by Wen et al., using a model of marrow radiation damage . We previously showed the ability of MSC EVs, to promote proliferation and to protect from apo- ptosis murine and human renal epithelial tubular cells [2, 43]. These effects were ascribed to a heterogeneous population of EVs containing both microvesicles and exosomes. We here demonstrated that different EV fractions separated by discon- tinuous iodixanol gradient based on their density, display dif- ferential pro-proliferative and anti-apoptotic activities on re- cipient renal tubular cells despite they were equally internal- ized by mTEC. EVs from CF1 and CF2 low and medium- density fractions induced proliferation on renal tubular cells, cultured in the absence of serum. However, the CF2 enriched- exosome population was more effective than the others com- bined fractions in the protection from apoptosis induced by hypoxia/reperfusion injury on renal target cells. High-density EVs in fraction CF3, conceivable containing a microvesicle- enriched population, did not induce significant protection from apoptosis and proliferation of renal tubular cells, also after enhancing the dose of EVs administered.
The activity of EVs has been at least in part ascribed to their miRNA content [42–45]. A subset of miRNA families, was specific signature of the biologically active CF2 medium- density fraction. Some of these miRNAs, such as miR-17-5p and miR-106a have been described together with miR-21 as differentially expressed during different phases of renal injury, being potential biomarkers for AKI . Interestingly, differ- ent miRNA* sequences were only detected in CF2 fraction, supporting the potential of exosomes to transfer miRNA* out- side the cells .
Several other miRNAs enriched in CF2 fraction, were re- ported as protective in the context of AKI including miR-34 , miR-125b , miR-199a-3p, miR-214  and miR- 127 . miR-21 and miR-29, enhanced in CF2, have been linked to the modulation of kidney fibrosis . miR-451 se- lectively present in the central fractions of the gradient, has been demonstrated to act together with miR-144, in the pro- tection against IRI in the heart . All the miRNAs analyzed were reduced in the CF1, supporting a less enrichment of miRNAs in this fraction.
Enrichment analysis of the pathways over-represented by the predicted targets of enriched/selective CF2 miRNAs showed strong correlation with metabolic, stem cell associated- and migration/inflammation processes. High rep- resentation of fatty acid biosynthesis and metabolism, Wnt signaling and pluripotent stem cells regulated-pathways, ECM-receptor interaction and TGF-β signaling pathway have been detected. Interestingly, downregulation of fatty acid ox- idation in tubular cells has been described as a key component of the pathogenesis of AKI . Moreover, stem cell and Wnt/β-catenin signaling pathways have been also correlated with repair processes after ischemic AKI .
Proteins are together with RNA, important effectors for the EV activity. By protein array, we identified a protein cargo of 413 proteins detected in all the different EV fractions. Cross match with Vesiclepedia database  identified 205 classical EV proteins detected in all the EV fractions. Proteins of the CF proteome were mainly incorporated as signaling molecules, receptors and cell adhesion molecules, recently connected with the potential therapeutic effects of the microvesicles from MSCs . Defense/immunity proteins and enzyme
Stem Cell Rev and Rep
Table 5 List of proteins over-expressed in CF2 (FC CF2 > 1.5) List of proteins over-expressed in CF2 in respect to CF3
hCG alpha FGFBP
GDNF Thrombospondin4 PLUNC
Thymidine Kinase1 SPINK1
IL3 R alpha Trypsin 1 OSM CXCR3 ENA78 MCSF IL17R CXCR2 TFPI CCR5 TOPORS
BDNF Frizzled6 PlGF IL17RC HGF
CCR8 CCR6 IGFBP1 TNFRSF1A CD56
IL12 p40 proGlucagon Chem R23 NT4
ICAM5 NeurokininA CD14 ALCAM CKMB HADHA TNFRSF17 TNFRSF5 CCR5 Frizzled5 BDNF
APJ IFNgamma R1 ApoC1 CD117
CCL11 VEGF R3 11bHSD1 TSLP R Neuropilin2 FGF R4 VWF S100A6 S100A10 TNFRSF11B IL28A Tissue Factor Gastrin BMP2 TNFSF8
SRMS BCAM CCR6 CCL16 CD56 HBEGF proGlucagon MMP25 ErbB2 CCL11 Protein p65 FAP
HSP10 CD80 TNFSF3 ErbB3
Activin B betaCatenin TNFRSF12 GPNMB
Leptin (OB) CCR1
IFNalpha / beta R1 LRP6
DANCE NeuroD1 CCR2 ApoE3
IL15 R alpha BMP2 TNFSF8 Activin B betaCatenin IL1 R3 ICAM2 LRP6 DANCE NeuroD1 IL17C
FGF9 CCR2 GLO1 PDX1 ApoE3
List of proteins over-expressed in CF2 in respect to CF1
Thrombospondin1 Fibronectin ApoA4
Factor XIII A Clusterin TYRO10
ALPP Ubiquitin + 1 Btk TNFRSF21 BAI1 Activin A SERPINA6 Beta IGH3 CA 125 ABL1 TIMP2 ALK4
Plasminogen hCG alpha PTHLP CNDP1 APC FGFBP
CA 153 ACTH
PAI1 Cardiotrophin1 VEGFC TRKB
Thymidine Kinase1 C9
modulators class of proteins were also detected in the CF proteome. Moreover, they were classified as extracellular, plasma membrane, matrix and vesicle proteins. Pathway anal- ysis of the CF proteome reveals high representation of path- ways such as interleukin mediated signaling and chemokine/ cytokine mediated inflammation, Wnt signaling, TGF-β and angiogenesis pathways. MSC CF proteome was relevantly enriched in anti- and pro-inflammatory cytokines.
Moreover, EVs from all CF fractions expressed the chemokines receptors of cell of origin, such as CXCR1, CXCR6, CXCR4, CXCR3, CCR3 and CCR7, known to be involved in the MSC migration to the sites of inflammation
. Eph A4 and different subclasses of the Eph B receptors, regulators of MSC attachment/migration  and inhibitors of T cell proliferation , were also present in all EV fractions. Moreover, numerous pro-angiogenic and pro-migratory mol- ecules such as the soluble factors VEGF, TGF-β, ΙL-8 and PDGF and PDGFRα/β were compartmentalized in all the CF fractions. Interestingly, increasing angiogenesis together with renal blood flow have been recently observed in IRI rats treated with adipose mesenchymal stem cell exosomes and associated with the reversion of the kidney damage . The catalytic subunit PI3-Kinase p85-beta, detected among the 50 major expressed proteins of CF proteome, has been
Stem Cell Rev and Rep (2017) 13:226–243
Fig. 8 Differential expression analysis of proteins in EVs from different CF fractions. a David GO-MF terms over-represented by the proteins upregulated in the medium density CF2 in respect to the CF3 high density fraction (n = 75, FC CF2 > 1.5, outer chart). The same analysis was conducted for proteins upregulated in CF2 in respect to the low- density fraction CF1 (n = 92, FC CF2 > 1.5, inner chart). (B-C)
recently correlated with cell cycle re-entry and proliferation via PI3K/Akt signaling pathway activation .
The protein associated biological processes overrepresent- ed in CF2, which was the most biologically active fraction, included positive regulators of proliferation, response to wounding and cell-cell signaling possibly accounting for the pro-proliferative activity of CF EVs. Migration/proliferation processes were also over-represented by the enriched miRNAs in the medium-density CF2 fraction, showing a col- laborative pattern of proteins and miRNAs shuttled by MSCs in promoting regenerative processes. Interestingly, TGF-β pathway crucial in AKI progression to chronic kidney disease  was overrepresented both by miRNAs and proteins com- partmentalized in the CF2 fraction, suggesting the regulation of this pathway by MSC EVs treatment.
Of interest, the major cytokine enriched in CF2 fraction was the IL-13 that together with IL-10 and IL-4, mediates differentiation of monocytes in non-inflammatory (M2)
Distribution analysis of the proteins enriched in CF2 in respect to CF3 b and to CF1 c, in different biological processes (GO-BB). The main representative classes were: signal transduction via receptor interaction, cell proliferation, response to wounding, cell-cell signaling and inflammatory response. (P-value < 0,01, FDR 1%)
macrophages . Recently, MSCs have shown the ability to educate macrophages to acquire an anti-inflammatory M2 phenotype, promoting kidney repair in rhabdomyolysis- induced acute kidney injury . Further, CF2 EVs contained high levels of ApoC3 and ApoA4, recently detected inside the EVs produced by umbilical cord blood-derived MSCs with other family members and involved in inflammation and tis- sue repair .
In conclusion, this study demonstrated that EVs derived from MSCs are heterogeneous with specific signatures ac- counting for the biological activity of different EV popula- tions. Proteins and miRNAs shuttled by MSCs showed a col- laborative pattern in crucial processes activated after injury, such as metabolic, stem cell, inflammation/migration and an- giogenic related processes. Moreover, the medium-density CF2 fraction containing exosome-enriched population of EVs showed the best activity in promoting renal protection from injury in vitro. This fraction was enriched in miRNAs
Stem Cell Rev and Rep (2017) 13:226–243
and proteins associated with biological processes fundamental in kidney regeneration.
Acknowledgements This work was supported by Associazione Italiana per la Ricerca sul Cancro (AIRC) and National Center For Advancing Translational Sciences of the National Institutes of Health under Award number UH2TR000880, UH3TR000880-03S1. The con- tent is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Compliance with Ethical Standards
RNA from immune cell-derived vesicles uncovers the selective incorporation of small non-coding RNA biotypes with potential regulatory functions. Nucleic Acids Research, 40(18), 9272–9285.
- Konala, V. B., Mamidi, M. K., Bhonde, R., Das, A. K., Pochampally, R., & Pal, R. (2016). The current landscape of the mesenchymal stromal cell secretome: a new paradigm for cell-free regeneration. Cytotherapy, 18(1), 13–24.
- Chiasserini, D., Mazzoni, M., Bordi, F., et al. (2015). Identification and partial characterization of two populations of Prostasomes by a combination of dynamic light scattering and proteomic analysis. The Journal of Membrane Biology, 248(6), 991–1004.
- Guduric-Fuchs, J., O’Connor, A., Camp, B., O’Neill, C. L., Medina, R. J., & Simpson, D. A. (2012). Selective extracellular vesicle- mediated export of an overlapping set of microRNAs from multiple cell types. BMC Genomics, 13, 357. doi:10.1186/1471-2164-13- 357.
- Ji, H., Chen, M., Greening, D. W., et al. (2014). Deep sequencing of RNA from three different extracellular vesicle (EV) subtypes re- leased from the human LIM1863 colon cancer cell line uncovers distinct miRNA-enrichment signatures. PloS One, 9(10), e110314.
- Lenassi, M., Cagney, G., Liao, M., et al. (2010). HIV Nef is secreted in exosomes and triggers apoptosis in bystander CD4+ T cells. Traffic, 11(1), 110–122.
- Lai, R. C., Arslan, F., Lee, M. M., et al. (2010). Exosome secreted by MSCs reduces myocardial ischemia/reperfusion injury. Stem Cell Research, 4(3), 214–222.
- Greening, D. W., Xu, R., Ji, H., Tauro, B. J., & Simpson, R. J. (2015). A protocol for exosome isolation and characterization: eval- uation of ultracentrifugation, density-gradient separation, and im- munoaffinity capture methods. Methods in Molecular Biology, 1295, 179–209.
- Taylor, D. D. (2015). Isolation and molecular characterization of extracellular vesicles. Methods, 87, 1–2. doi:10.1016/j. ymeth.2015.08.006.
- Tomasoni, S., Longaretti, L., Rota, C., et al. (2013). Transfer of growth factor receptor mRNA via exosomes unravels the regener- ative effect of mesenchymal stem cells. Stem Cells and Development, 22(5), 772–780.
- Collino, F., Bruno, S., Incarnato, D., et al. (2015). AKI recovery induced by mesenchymal stromal cell-derived extracellular vesicles carrying MicroRNAs. Journal of the American Society of Nephrology, 26(10), 2349–2360.
- Herrera, M. B., Fonsato, V., Bruno, S., et al. (2013). Human liver stem cells improve liver injury in a model of fulminant liver failure. Hepatology, 57(1), 311–319.
- Tauro, B. J., Greening, D. W., Mathias, R. A., et al. (2012). Comparison of ultracentrifugation, density gradient separation, and immunoaffinity capture methods for isolating human colon cancer cell line LIM1863-derived exosomes. Methods, 56(2), 293–304.
- Kowal, J., Arras, G., Colombo, M., et al. (2016). Proteomic com- parison defines novel markers to characterize heterogeneous popu- lations of extracellular vesicle subtypes. Proceedings of the National Academy of Sciences of the United States of America, 113(8), E968–E977.
- Mestdagh, P., Van Vlierberghe, P., De Weer, A., et al. (2009). A novel and universal method for microRNA RT-qPCR data normal- ization. Genome Biology, 10(6), R64. doi:10.1186/gb-2009-10-6-
The authors indicate no potential conflicts of interest.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
- Ratajczak, J., Miekus, K., Kucia, M., et al. (2006). Embryonic stem cell-derived microvesicles reprogram hematopoietic progenitors: evidence for horizontal transfer of mRNA and protein delivery. Leukemia, 20(5), 847–856.
- Bruno, S., Grange, C., Deregibus, M. C., et al. (2009). Mesenchymal stem cell-derived microvesicles protect against acute tubular injury. Journal of the American Society of Nephrology, 20(5), 1053–1067.
- Tan, C. Y., Lai, R. C., Wong, W., Dan, Y. Y., Lim, S. K., & Ho, H. K. (2014). Mesenchymal stem cell-derived exosomes promote he- patic regeneration in drug-induced liver injury models. Stem Cell Research & Therapy, 5(3), 76–89.
- Zhu, Y. J., Sun, J., & Du, T. (2015). Microvesicles derived from human umbilical cord mesenchymal stem cells facilitate tubular epithelial cell dedifferentiation and growth via hepatocyte growth factor induction. PloS One, 10(3), e0121534.
- Nakamura, Y., Miyaki, S., Ishitobi, H., et al. (2015). Mesenchymal- stem-cell-derived exosomes accelerate skeletal muscle regenera- tion. FEBS Letters, 589(11), 1257–1265.
- Deregibus, M. C., Cantaluppi, V., Calogero, R., et al. (2007). Endothelial progenitor cell derived microvesicles activate an angio- genic program in endothelial cells by a horizontal transfer of mRNA. Blood, 110(7), 2440–2448.
- Valadi, H., Ekström, K., Bossios, A., Sjöstrand, M., Lee, J. J., & Lötvall, J. O. (2007). Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nature Cell Biology, 9(6), 654–659.
- Collino, F., Deregibus, M. C., Bruno, S., et al. (2010).
Microvesicles derived from adult human bone marrow and tissue
specific mesenchymal stem cells shuttle selected pattern of r64. miRNAs. PloS One, 5(7), e11803.
- Lai, R. C., Tan, S. S., Teh, B. J., et al. (2012). Proteolytic potential of the MSCS exosome proteome: implications for an exosome- mediated delivery of therapeutic proteasome. The International Journal of Prosthodontics. doi:10.1155/2012/971907.
- Nolte-‘t Hoen, E.N., Buermans, H.P., Waasdorp, M., Stoorvogel, W., Wauben, M.H., ‘t Hoen, P.A. (2012) Deep sequencing of
- Vlachos, I. S., Zagganas, K., Paraskevopoulou, M. D., et al. (2015). DIANA-miRPath v3.0: deciphering microRNA function with ex- perimental support. Nucleic Acids Research, 43(W1), W460– W466.
- Pathan, M., Keerthikumar, S., Ang, C. S., et al. (2015). FunRich: a standalone tool for functional enrichment analysis. Proteomics, 15, 2597–2601.
Stem Cell Rev and Rep (2017) 13:226–243 243
- Bhatt, K., Zhou, L., Mi, Q. S., Huang, S., She, J. X., & Dong, Z. (2010). MicroRNA-34a is induced via p53 during cisplatin neph- rotoxicity and contributes to cell survival. Molecular Medicine, 16, 409–416.
- Godwin, J. G., Geb, X., Stephana, K., Jurischb, A., Tulliusb, S. G., & Iacominia, J. (2010). Identification of a microRNA signature of renal ischemia reperfusion injury. Proceedings of the National Academy of Sciences, 107(32), 14339–14344.
- Aguado-Fraile, E., Ramos, E., Saenz-Morales, D., et al. (2012). miR-127 protects proximal tubule cells against ischemia/reperfu- sion: identification of kinesin family member 3B as miR-127 target. PloS One, 7, e44305.
- Patel, V., & Noureddine, L. (2012). MicroRNAs and fibrosis. Current Opinion in Nephrology and Hypertension, 21(4), 410–416.
- Joo, M. S., Lee, C. G., Koo, J. H., & Kim, S. G. (2013). miR-125b transcriptionally increased by Nrf2 inhibits AhR repressor, which protects kidney from cisplatin-induced injury. Cell Death &
Disease, 4, e899.
- Kaucsár, T., Révész, C., & Godó, M. (2013). Activation of the miR-
17 family and miR-21 during murine kidney ischemia-reperfusion
injury. Nucleic Acid Therapeutics, 23(5), 344–354.
- Chiabotto, G., Bruno, S., Collino, F., & Camussi, G. (2016). Mesenchymal stromal cells epithelial transition induced by renal tubular cells-derived extracellular vesicles. PloS One, 11(7),
- Gatti, S., Bruno, S., Deregibus, M. C., et al. (2011). Microvesicles
derived from human adult mesenchymal stem cells protect against ischaemia-reperfusion-induced acute and chronic kidney injury. Nephrology, Dialysis, Transplantation, 26(5), 1474–1483.
- Zhang, Y., Chopp, M., Meng, Y., et al. (2015). Effect of exosomes derived from multipluripotent mesenchymal stromal cells on func- tional recovery and neurovascular plasticity in rats after traumatic brain injury. Journal of Neurosurgery, 122(4), 856–867.
- Mitsialis, S. A., & Kourembanas, S. (2016). Stem cell-based ther- apies for the newborn lung and brain: possibilities and challenges. Seminars in Perinatology, 40(3), 138–151.
- Lai, R. C., Yeo, R. W., & Lim, S. K. (2015). Mesenchymal stem cell exosomes. Seminars in Cell & Developmental Biology, 40, 82–88.
- Hagiwara, K., Katsuda, T., Gailhouste, L., Kosaka, N., & Ochiya, T. (2015). Commitment of annexin A2 in recruitment of microRNAs into extracellular vesicles. FEBS Letters, 589(24 Pt B), 4071–4078.
- Palma, J., Yaddanapudi, S. C., Pigati, L., et al. (2012). MicroRNAs are exported from malignant cells in customized particles. Nucleic
Acids Research, 40(18), 9125–9138.
- Xu, R., Greening, D. W., Rai, A., Ji, H., & Simpson, R. J. (2015).
Highly-purified exosomes and shed microvesicles isolated from the human colon cancer cell line LIM1863 by sequential centrifugal ultrafiltration are biochemically and functionally distinct. Methods, 87, 11–25.
- Aliotta, J. M., Pereira, M., Wen, S., et al. (2016). Exosomes induce and reverse monocrotaline-induced pulmonary hypertension in mice. Cardiovascular Research, 110(3), 319–330.
- Wen, S., Dooner, M., Cheng, Y., et al. (2016). Mesenchymal stro- mal cell-derived extracellular vesicles rescue radiation damage to murine marrow hematopoietic cells. Leukemia. doi:10.1038 /leu.2016.107.
- Lindoso, R. S., Collino, F., Bruno, S., et al. (2014). Extracellular vesicles released from mesenchymal stromal cells modulate miRNA in renal tubular cells and inhibit ATP depletion injury. Stem Cells and Development, 23(15), 1809–1819.
- Momen-Heravi, F., Bala, S., Kodys, K., & Szabo, G. (2015). Exosomes derived from alcohol-treated hepatocytes horizontally transfer liver specific miRNA-122 and sensitize monocytes to LPS. Scientific Reports, 5, 9991.
- Wang, B., Yao, K., & Huuskes, B. M. (2016). Mesenchymal stem cells deliver exogenous MicroRNA-let7c via exosomes to attenuate renal fibrosis. Molecular Therapy. doi:10.1038/mt.2016.90.
- Zhang, X., Wang, X., & Zhu, H. (2010). Synergistic effects of the GATA-4-mediated miR-144/451 cluster in protection against sim- ulated ischemia/reperfusion-induced cardiomyocyte death. Journal of Molecular and Cellular Cardiology, 49(5), 841–850.
- Simon, N., & Hertig, A. (2015). Alteration of fatty acid oxidation in tubular epithelial cells: from acute kidney injury to renal Fibrogenesis. Frontiers in Medicine (Lausanne), 2, 52. doi:10.3389/fmed.2015.00052.
- Zhou, D., Tan, R. J., Fu, H., & Liu, Y. (2016). Wnt/β-catenin signaling in kidney injury and repair: a double-edged sword. Laboratory Investigation, 96(2), 156–167.
- Kalra, H., Simpson, R. J., Ji, H., et al. (2012). Vesiclepedia: a com- pendium for extracellular vesicles with continuous community an- notation. PLoS Biology, 10(12), e1001450.
- Kim, H. S., Choi, D. Y., Yun, S. J., et al. (2012). Proteomic analysis of microvesicles derived from human mesenchymal stem cells. Journal of Proteome Research, 11(2), 839–849.
- Du, T., & Zhu, Y. J. (2014). The regulation of inflammatory medi- ators in acute kidney injury via exogenous mesenchymal stem cells. Mediators of Inflammation, 2014, 261697. doi:10.1155/2014 /261697.
- Arthur, A., Zannettino, A., Panagopoulos, R., et al. (2011). EphB/ ephrin-B interactions mediate human MSCS attachment, migration and osteochondral differentiation. Bone, 48(3), 533–542.
- Nguyen, T. M., Arthur, A., Hayball, J. D., & Gronthos, S. (2013). EphB and ephrin-B interactions mediate human mesenchymal stem cell suppression of activated T-cells. Stem Cells and Development, 22(20), 2751–2764.
- Lin, K. C., Yip, H. K., & Shao, P. L. (2016). Combination of adipose-derived mesenchymal stem cells (ADMSCS) and ADMSCS-derived exosomes for protecting kidney from acute ischemia-reperfusion injury. International Journal of Cardiology, 216, 173–185.
- Gao, J., Zhou, X. L., Kong, R. N., Ji, L. M., He, L. L., & Zhao, D. B. (2016). microRNA-126 targeting PIK3R2 promotes rheumatoid arthritis synovial fibroblasts proliferation and resistance to apopto- sis by regulating PI3K/AKT pathway. Experimental and Molecular Pathology, 100(1), 192–198.
- Tanaka, T. (2016) A mechanistic link between renal ischemia and fibrosis. Medical Molecular Morphology, 2016 Jul 20. [Epub ahead of print].
- Ricardo, S. D., van Goor, H., & Eddy, A. A. (2008). Macrophage diversity in renal injury and repair. Journal of Clinical Investigation, 118(11), 3522–3530.
- Geng, Y., Zhang, L., & Fu, B. (2014). Mesenchymal stem cells ameliorate rhabdomyolysis-induced acute kidney injury via the ac- tivation of M2 macrophages. Stem Cell Research & Therapy, 5(3), 80–94.
- Kilpinen, L., Impola, U., & Sankkila, L. (2013). Extracellular mem- brane vesicles from umbilical cord blood-derived MSCS protect against ischemic acute kidney injury, a feature that is lost after inflammatory conditioning. Journal of Extracellular Vesicles. doi:10.3402/jev.v2i0.21927.