Original Research ARTICLE
Melanoma-Derived Exosomes Induce PD-1 Overexpression and Tumor Progression via Mesenchymal Stem Cell Oncogenic Reprogramming
- 1Laboratory of Microscopic Image Analysis and Machine Learning, Institute of Biochemistry, Biological Research Centre of Hungarian Academy of Sciences, Szeged, Hungary
- 2Doctoral School of Interdisciplinary Sciences, Faculty of Medicine, University of Szeged, Szeged, Hungary
- 3Department of Dermatology and Allergology, University of Szeged, Szeged, Hungary
- 4Department of Immunology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
- 5Laboratory of Functional Genomics, Institute of Genetics, Biological Research Centre of Hungarian Academy of Sciences, Szeged, Hungary
- 6Laboratory of Proteomics Research, Institute of Biochemistry, Biological Research Centre of Hungarian Academy of Sciences, Szeged, Hungary
- 7Artificial Chromosome and Stem Cell Research Laboratory, Institute of Genetics, Biological Research Centre of Hungarian Academy of Sciences, Szeged, Hungary
- 8Sequencing Platform, Institute of Biochemistry, Biological Research Centre of Hungarian Academy of Sciences, Szeged, Hungary
- 9Department of Pharmacology and Pharmacotherapy, Faculty of Medicine, University of Szeged, Szeged, Hungary
- 10Department of Biochemistry and Molecular Biology, Faculty of Science and Informatics, University of Szeged, Szeged, Hungary
- 11Advanced Optical Imaging Group, Department of Optics and Quantum Electronics, Faculty of Science and Informatics, University of Szeged, Szeged, Hungary
- 12Atomic Force Microscope Laboratory, Institute of Biophysics, Biological Research Centre of Hungarian Academy of Sciences, Szeged, Hungary
- 13MTA-SE Immuno-proteogenomics Extracellular Vesicle Research Group, Department of Genetics, Cell- and Immunobiology, Faculty of Medicine, Semmelweis University, Budapest, Hungary
- 14Department of Oral Biology and Experimental Dental Research, Faculty of Dentistry, University of Szeged, Szeged, Hungary
Recently, it has been described that programmed cell death protein 1 (PD-1) overexpressing melanoma cells are highly aggressive. However, until now it has not been defined which factors lead to the generation of PD-1 overexpressing subpopulations. Here, we present that melanoma-derived exosomes, conveying oncogenic molecular reprogramming, induce the formation of a melanoma-like, PD-1 overexpressing cell population (mMSCPD-1+) from naïve mesenchymal stem cells (MSCs). Exosomes and mMSCPD-1+ cells induce tumor progression and expression of oncogenic factors in vivo. Finally, we revealed a characteristic, tumorigenic signaling network combining the upregulated molecules (e.g., PD-1, MET, RAF1, BCL2, MTOR) and their upstream exosomal regulating proteins and miRNAs. Our study highlights the complexity of exosomal communication during tumor progression and contributes to the detailed understanding of metastatic processes.
Heterogeneous tumor tissue is comprised of a wide variety of collocated cells. Their spatiotemporal co-existence facilitates direct communication between them. Cancer cells contribute to tumor niche formation not only by soluble factor production and receptor-ligand interactions (1), but also by releasing vesicles whose molecular contents add up to a complex information package. Previous studies demonstrated that cultured human tumor cells release extracellular vesicles such as exosomes of 20–120 nm diameters (2). Among others, exosomes carry structural and signaling proteins, MHC molecules, cell surface molecules typically associated with apoptosis, and mRNAs/miRNAs with multiple functions (3). Therefore, these exosomal-molecular-patterns, as unique entities of the complex intercellular communication, are not independent of the quality or state of the mother cell (4, 5). Exosomes have been recognized long ago, but their identification, characterization, and isolation are still under intense investigation. Further, whereas the definition of exosomes is based on the pathway of biogenesis, oncosomes form a functional class of extracellular vesicles. Indeed, oncosomes are suggested being capable of carrying and conveying tumor-related information (6) and might have a significant role in formation of tumor microenvironment (7).
Mesenchymal stem cells (MSCs) were first described as stromal cells of the bone marrow with multipotent differentiation potential and characteristic immunomodulatory effects (8). In relation to their immunological and differentiation properties, there is a debate about the role of MSCs in tumor progression (9). Indeed, the cellular fate could depend on the cancer type and also on the status of the affected MSCs. Activated MSCs can secrete pro-angiogenic soluble factors and are able to differentiate to vessel wall associated cells (10) or to cancer associated fibroblasts (CAFs) (11). Furthermore, Baglio et al. showed that tumor secreted extracellular vesicles promote osteosarcoma progression via TGFβ signaling induced IL-6 production by MSCs (12). Moreover, Peinado et al. demonstrated that exosomes from highly metastatic melanomas increase the metastatic behavior of primary tumors by permanently educating bone marrow progenitors through the receptor tyrosine kinase MET (13).
Malignant melanoma is a highly invasive, metastatic cancer with poor prognosis and survival rate. Since melanoma has been long recognized as a highly immunogenic tumor, therapeutic approaches target different immunological mechanisms to treat patients with this aggressive skin cancer. In the past decade, besides IL-2 and CTLA-4 therapies, the PD-1:PD-L1 blockade proved to be an effective treatment in metastatic melanoma (14). Programmed cell death protein 1 (PD-1), mainly expressed by immune cells such as activated T cells, dendritic cells and macrophages is a cell surface receptor with a central role in modulation of T cell responses. PD-1 binding to its ligand, PD-L1, can result in apoptosis and anergy of immune cells.
Notably, PD-1 expression is not a default property of melanoma cells. Of further importance, recently Kleffel et al. demonstrated that PD-1 overexpressing melanoma cell subpopulations are especially aggressive, and that melanoma PD-1:PD-L1 interactions modulate mTOR signaling, whereby they interfere with programmed cell death (15). This work is remarkably instrumental with respect to the consequences of PD-1 overexpression and experimental inhibition of PD-1 functions (Pdcd1-shRNA or anti-PD-1 antibody). However, it does not define the factors the activation of which lead to the generation of PD-1 overexpressing subpopulations. This question, though, is an especially intriguing one, if one considers that one of the most effective immunotherapies available today is based on the blockade of the aforementioned PD-1:PD-L1 interaction (16).
Based on the above, in the current study, we aimed at specifically and systematically addressing the following questions:
1. Can we unambiguously define the cellular and molecular signs of melanoma-derived exosome-induced, intercellular communication-mediated malignant transformation of MSC cultures?
2. Can we detect the melanoma-derived exosome-induced tumor progression in vivo?
3. What alteration can be detected in the expression of PD-1, a melanoma progression marker and therapeutic target, upon exposure to melanoma-derived exosomes?
Here, we present the first piece of evidence that melanoma exosomes, by conveying a complex oncogenic molecular reprogramming, induce the formation of a PD-1 overexpressing cell population (melanoma-like MSCPD-1+; mMSCPD-1+) from naïve MSCs. These mMSCPD-1+ cells represent a new entity with melanoma-like gene expression profile and phenotypic properties. Further, exosomes and exosome-activated mMSCPD-1+ cells induce rapid tumor progression in vivo, due to their strong expression of oncogenic dominance factors and decreased susceptibility to programmed cell death.
Materials and Methods
Cell Culture and Cell Line
B16F1 mouse melanoma cell line was obtained from ECACC and cultured in Dulbecco's modified Eagle's medium (DMEM, Lonza, Basel, Switzerland) containing 10% fetal bovine serum (FBS) (Euroclone, Milano, Italy), 1% sodium pyruvate, 1% MEM non-essential amino acids, 1% MEM vitamin solution, 2 mM L-glutamine, and 1% Penicillin-Streptomycin-Amphotericin B Mixture (all from Lonza).
Murine MSCs were isolated from C57BL/6N (Charles River Laboratories, Sulzfeld) 8 week-old male mice (17). Abdominal inguinal fat pads were excised, rinsed with RPMI 1640 medium, transferred to sterile tissue culture dishes and mechanically dissociated. The dissociated tissue was resuspended in RPMI 1640 containing 100 μg/ml collagenase (Sigma-Aldrich, St. Louis, MO) and incubated at 37°C for 50 min. Collagenase was neutralized with a growth medium containing 10% FBS. After centrifugation at 470 g for 15 min, cell pellets were resuspended and washed in the culture medium. After the centrifugation, cell pellets were resuspended in complete MesenCult medium (Stemcell Technologies, Vancouver, British Columbia, Canada) and filtered through a 100 μm cell strainer (BD Biosciences, Franklin Lakes, NJ) to tissue culture dishes, and cultured for 48 h at 37°C in 5% CO2 and 90% humidity. Unattached cells and debris were then removed and fresh medium was added to the adherent cells. The cells were cultured to 80% confluence before being released with trypsin–EDTA and sub-cultured.
Purity of MSC cultures were checked by flow cytometry analysis using the Mouse Multipotent Mesenchymal Stromal Cell Marker Antibody Panel (R&D Systems, Minneapolis, MN) according to the manufacturer's instructions.
Isolation and Characterization of the Melanoma-Derived Exosomes
Exosome Isolation by Filtration and Differential Centrifugation
Exosomes were isolated by adapting the protocol of Peinado et al. (13) from melanoma cell culture supernatant. B16F1 supernatants were harvested, supplemented by complete protease inhibitor cocktail (Roche, Mannheim, Germany) and centrifuged at 780 g for 5 min at 4°C to remove intact cells. Then, the supernatants were centrifuged at 3,900 g for 15 min at 4°C and filtered through a 0.2 μm membrane (Millipore, Billerica, MA) to remove larger cell fragments and microvesicles. Exosomes were pelleted by ultracentrifugation at 150,000 g (T-1270 rotor at 40,500 rpm) for 1 h at 4°C. The pellet was washed twice and resuspended in DPBS and stored at −80°C.
The concentration of exosomal proteins was determined using a Pierce BCA Protein assay kit (Thermo Fisher Scientific, Waltham, MA) and a Benchmark Microplate Reader (Bio-Rad, Hercules, CA) according to the manufacturer's instructions.
Determination of Shape and Size of the Isolated Vesicles
Scanning electron microscopy (SEM)
Five microliters of isolated exosomes in DPBS were spotted on glass coverslips coated with 0.01% Poly-L-Lysine (Sigma-Aldrich) and incubated overnight at 4°C. Using sterile forceps, coverslips were transferred into a plastic plate. Exosomes were fixed by gently adding 2% paraformaldehyde buffer pH 7.2 diluted in DPBS for 30 min. The coverslips were washed twice with DPBS and dehydrated with a graded ethanol series (40, 60, 80, and 100% ethanol, each for 5 min). The samples were dried with a critical point dryer (QUORUM K850, Quorum Technologies Ltd, Laughton, UK) and the coverslips were mounted onto a microscope stub at a time using carbon tape, followed by 7 nm gold coating (QUORUM Q150, Quorum Technologies Ltd) and observed under a field-emission scanning electron microscope (JEOL JSM-7100F/LV).
Atomic force microscopy (AFM)
Exosomes were let to adsorb to the freshly cleaved muscovite mica (SPI-ChemTM Mica Sheets, Structure Probe, Inc., West Chester, PA) surface directly from DPBS. All AFM measurements were carried out with an Asylum MFP-3D head and Molecular Force Probe Controller (Asylum Research, Santa Barbara, CA), using Asylum Research MFP-3D program (version 15.09.112) written in Igor Pro software (version 6.37, Wavemetrics, Lake Oswego, OR). Image procession and data calculation were made using the same software. Silicon rectangular cantilevers (OMLC-AC240TS, Olympus Optical Co., Ltd., Tokyo, Japan) with the typical spring constant of 2 N/m were used in dry condition and silicon nitride rectangular cantilevers with “V” shaped tips (Bio-Lever BL-RC150VB, Olympus Optical Co. Ltd.) with the typical spring constant of 0.03 N/m in buffer solution. Typically, 512 × 512 points was taken at 1 line/s scan rate. The measurements presented here are amplitude images and height profile.
Detection of Exosomal Markers by Western Blotting
Protein samples were resuspended in 4x sample buffer (NuPAGE LDS Sample Buffer (4X), NuPAGE Sample Reducing Agent (10X), Thermo Fisher Scientific), boiled at 96°C for 10 min, and immediately cooled on ice. Electrophoresis of the proteins was performed using 4–12% Bis-Tris Protein Gels (NuPAGE Novex, Thermo Fisher Scientific), ProSieve Color Protein Markers (Lonza), and electrophoresis buffer (NuPAGE MOPS SDS Running Buffer, Thermo Fisher Scientific) at 200 V and 0.03 A for 40 min.
The proteins of electrophoresed gels were transferred to Immobilon transfer membrane (Millipore, Darmstadt, Germany) using transfer buffer [NuPAGE Transfer Buffer (20X)] at 30 V and 170–110 mA for 60 min. Membranes were blocked in TBST buffer (25 mM Tris-HCl, 150 mM NaCl, 0.05% Tween-20, pH 7.2) containing 5% non-fat milk for 60 min at room temperature (RT). After blocking, membranes were incubated with each primary antibody in TBST buffer containing 1% non-fat milk overnight at 4°C. Membranes were washed three times for 10 min with TBST buffer and incubated for 60 min at room temperature with secondary antibody in TBST buffer containing 1% non-fat milk. Membranes were washed three times for 10 min with TBST buffer.
For exosomal marker identification western blot analyses were performed with an anti-CD63 polyclonal antibody (1:250 dilution; Biorbyt, Cambridge, UK) anti-CD9 monoclonal antibody (clone EPR2949, 1:500 dilution, LifeSpan Biosciences, Seattle, WA), anti-CD81 monoclonal antibody (clone: EAT2, 1:1,000 dilution, LifeSpan Biosciences, Seattle, WA), and anti-HSP70 monoclonal antibody (clone: C92F3A-5, 1:8,000 dilution; Enzo Life Sciences, Farmingdale, NY), anti-rabbit IgG HRP-Conjugated antibody (1:1,000 dilution, R&D Systems), anti-Hamster IgG HRP-Conjugated antibody (1:30,000 dilution, Thermo Fisher Scientific). Bound antibodies were visualized by chemiluminescence using an ECL Plus Western Blotting detection system (Advansta, Menlo Park, CA). Immunoreactive signals were detected by using LI-COR ODYSSEY Fc (Dual-mode imaging system) imager followed by analysis with Odyssey v1.2, Image Studio Lite v5.2.
Identification of Exosomal Proteins by Mass Spectrometry
Exosomal proteins (24 μg) were separated in 4–12% Bis-Tris Protein Gels and stained with Coomassie blue (0.1% Coomassie Brilliant Blue R-250, 50% methanol and 10% acetic acid).
Each lane was cut to 12–12 equal bands and subjected to in-gel digestion. Gel bands were diced to smaller pieces, and the SDS and CBB dyes were washed out with 3 × 50 μl 25 mM ammonium-bicarbonate (ABC)/50% Acetonitrile (AcN). After reduction with DTT (1,4-dithiothreitol, Sigma-Aldrich; 20 μl, 10 mM DTT in 25 mM ABC) at 56°C for 30 min, and alkylation with IAM (iodoacetamide, Sigma- Aldrich; 20 μl, 55 mM IAM in 25 mM ABC) at RT in the dark for 30 min, the gel samples were dried in a vacuum centrifuge and after that rehydrated in 20 μl of trypsin (Sequencing Grade Modified Trypsin, Promega, Fitchburg, WI; 5 ng/μl in 25 mM ABC) and incubated at 37°C. The digestion was stopped after 4 h by lowering the pH of the buffer below 3, by adding 2 μl of 10% formic acid. Tryptic peptides were extracted from the gel with 3 × 50 μl of 2% formic acid in 50% AcN and dried. Prior mass spectrometric analysis, all samples were redissolved in 50–50 μl of 0.1% formic acid (FA).
Samples were analyzed on an LTQ-Orbitrap Elite (Thermo Fisher Scientific) mass spectrometer on-line coupled with a nanoHPLC (nanoAcquity, Waters, Milford, MA) system. 5–5 μl of the in-gel digests were loaded (for 3 min at 8 μl/min flow, using 0.1% FA in 3% Acetonitrile−97% water) onto a reversed phase trap column (Waters, Symmetry C18, 0.180 × 20 mm) and separated on a C18 reversed phase (Waters, Milford, MA, BEH300C18 1.7 μm) nanocolumn (0.075 × 200 mm). The flow rate was 330 nl/min and a linear gradient was used from 3 to 40% B in 37 min (solvent A was 0.1% FA in water and solvent B was 0.1% FA in Acetonitrile).
The high voltage (1.2 kV) was applied through liquid junction between the chromatographic column and the non-coated silica nanospray emitter (NewObjective, Woburn, MA, 10 μm tip ID). The mass spectrometer operated in data-dependent mode: the survey mass spectra were detected in the orbitrap with high resolution (R = 60 k @ m/z: 400, mass range m/z: 380–1,400) and the most abundant multiply charged 20 peaks were selected for ion-trap fragmentation (NCE: 35%; activation q: 0.25; activation time: 10 ms; minimum signal intensity: 5,000 counts). The MSMS spectra were detected in the ion trap. Dynamic exclusion was used, the precursors were excluded for 15 s after the first fragmentation event.
Data analysis: searchable peaklists (mgf format) were extracted using ProteomeDiscoverer (ver:1.4 Thermo Fisher Scientific) and subjected to database search on our in-house ProteinProspector (ver: 5.14.1) search engine using the following parameters: parent ion tolerant: 5 ppm; fragment ion tolerant: 0.6 Da; Cys carbamidomethylation was set as constant and Met oxidation, cyclisation of peptide N-terminal Glu to pyroglutamic acid, protein N terminal acetylation were set as variable modifications. Only fully tryptic peptides were considered with maximum of 2 missed cleavage sites. The Mus musculus and Bos taurus protein sequences of the Uniprot (UniProtKB.06.11.2014) database completed with human keratins and pig trypsin, altogether 106,330 protein sequences were searched. For the false discovery rate (FDR) estimation, the searches were performed on the database concatenated with the randomized sequences. Protein identification was accepted if the ProteinProspector expectation value was <0.01 and the protein was identified with at least 2 unique peptides (expectation value <0.05 and score higher than 15). FDR values were <1% in all cases.
For functional validation, the resulted protein list was analyzed by the “Core Analysis” function included in Ingenuity Pathway Analysis (IPA, Quiagen Bioinformatics) software.
Isolation and Detection of Exosomal miRNAs
miRNA sequencing was performed using SOLiD Total RNA-Seq lit for Small RNA Libraries (Applied Biosystems now part of Thermo Fisher Scientific) according to the manufacturer's instructions. Purification was performed on 10% TBE-Urea gels stained with Sybr Gold nucleic acid gel stain (both from Invitrogen now part of Thermo Fisher Scientific). Final purification was performed using PureLink PCR Micro Kit (Invitrogen). Final libraries were quality checked using High Sense DNA kit on Bioanalyzer (all from Agilent, Santa Clara, CA). Concentration of each library was determined using the SOLID Library TaqMan Quantitation Kit (Life Technologies now part of Thermo Fisher Scientific). Each library was clonally amplified on SOLiD P1 DNA Beads by emulsion PCR (ePCR). Emulsions were broken with butanol, and ePCR beads enriched for template-positive beads by hybridization with magnetic enrichment beads. Template-enriched beads were extended at the 3′ end in the presence of terminal transferase and 3′ bead linker. Beads with the clonally amplified DNA were deposited onto SOLiD sequencing slide and sequenced on SOLiD 5500 Instrument using the 50-base sequencing chemistry.
Bioinformatic Analysis Raw data quality assessment, read trimming read mapping and miRNA expression profiling were carried out in CLC Genomics Workbench tool version 8.0.2 (CLC Bio now part of Qiagen, Venlo, Netherlands) using annotated Mus musculus miRNA sequences according to the miRBase release 21 as a mapping reference.
In vitro Experiments
6 × 104 cell/ml passage 2 MSCs were plated in cell culture dishes (1.5 × 104/cm2). After 24 h incubation, MSC cultures were exposed to B16F1-derived exosomes (40 μg/ml exosomal proteins; 1.5 × 1011 exosomes) at every 24 h. Samples were exposed to exosomes for 24, 48, 72, and 96 h and then harvested in method-competent buffers.
Visualization of Labeled Exosome Internalization in MSCs
To examine the uptake of exosomes by MSCs, cells were plated to black 24-well Visiplates (1 × 104 cells/well) and incubated for 24 h. The exosomes were labeled with Dil dye (1,1′-dioctadecyl-3,3,3′,3′-tetramethylindocarbocyanine perchlorate, PromoKine, Heidelberg, Germany) and the MSC cultures were labeled with DiO dye (3,3′-dioctadecyloxacarbo-cyanine perchlorate, PromoKine) according to the manufacturer's instructions. Dil-labeled exosomes were washed in DPBS by ultracentrifugation (at 150,000 g for 1 h at 4°C). Forty micrograms per milliliter DiL-labeled exosomes were added to DiO-labeled MSC cultures and the exosome uptake was followed for 24 h in the Celldiscoverer 7 automated live cell imaging system (Zeiss, Oberkochen, Germany). After 24 h, the cells were fixed with 4% paraformaldehyde solution and a nucleus staining was performed using DAPI (Life Technologies now part of Thermo Fisher Scientific). Then, 5 image z-stacks were acquired for both channels by Operetta High Content Screening System (Perkin Elmer, Waltham, MA). The stacks were maximum intensity projected and then analyzed automatically using a customized version of CellProfiler (18). Nuclei were detected with Otsu-adaptive threshold combined with diameter based filtering, then cytoplasms were identified with propagation method seeded from the nuclei and using the exosome channel. Exosomes were located with a customized version of A-trous wavelet transform based spot detection (19). Several wavelet levels were used to ensure the detection of exosomes with various size and then the overlaps were removed based on circularity measures. Finally, the exosome numbers per cell were identified using MATLAB programming, the diagrams were created in Microsoft Excel.
After 72 h incubation, exosome-exposed and control MSC cultures were dissociated with trypsin from the culture surface. Cells were washed in medium and counted in a Bürker chamber and a cell counter (Bio-Rad, TC10 Automated Cell Counter).
Detection of Apoptosis
Exosome-exposed MSCs and control cells were treated with 100 ng/ml mouse TNFα (R&D Systems). After 24 h incubation, cell death was determined by the Annexin V Apoptosis Detection Kit with PI (Biolegend, San Diego, CA) according to the manufacturer's recommendations. Samples were measured by FACS Calibur flow cytometer (BD Biosciences), data were analyzed by Flowing Software (Cell Imaging Core, Turku Center for Biotechnology, Finland) where percent of positive cells was determined by relative fluorescence intensity and the results were expressed as mean of percentage of positive cells (%) ±SD. Cells that are annexin-V/PI double positive show the sign of late apoptosis, while cells that are annexin-V positive and PI negative indicate early apoptosis. Annexin-V negative and PI positive cells are necrotic, viable cells are both annexin-V, and propidium iodide (PI) negative.
RNA Preparation, Melan-A, and Mitf Quantitative Real-Time PCR (QRT-PCR)
Total RNA of biological samples was purified using the Quick-RNA MiniPrep isolation kit of #R1054S (Zymo Research Irvine, CA). All the preparation steps were carried out according the manufacturer's instructions. RNA samples were stored at −80°C in the presence 30 U of Prime RNAse inhibitor (Fermentas, part of Thermo Fisher Scientific) for further analysis. The quantity of isolated RNA samples was checked by spectrophotometry (NanoDrop 3.1.0, Rockland, DE).
QRT-PCR was performed on a RotorGene 3000 instrument (Corbett Research, Sydney, Australia) with gene-specific primers and SybrGreen protocol to monitor gene expression. One microgram of total RNA was reverse transcribed with random primers using the High-Capacity cDNA Archive Kit (Applied Biosystems) according to the manufacturer's instructions in final volume of 30 μl. The temperature profile of the reverse transcription was the following: 10 min at room temperature, 2 h at 37°C, 5 min on ice, and finally 10 min at 75°C for enzyme inactivation. These steps were carried out in a Thermal Cycler machine (MJ Research, Marshall Scientific, Hampton, NH). After dilution with 30 μl of water, 1 μl of the diluted reaction mix was used as template in the QRT-PCR. Reactions were done with FastStart SYBR Green Master mix (Roche) according to the manufacturer's instructions at a final primer concentration of 250 nM under the following conditions: 15 min at 95°C, 40 cycles of 95°C for 15 s, 60°C for 25 s, and 72°C for 25 s. The fluorescence intensity of SybrGreen dye was detected after each amplification step. Melting temperature analysis was done after each reaction to check the quality of the products. Primers were designed using the online Roche Universal Probe Library Assay Design Center. The quality of the primers was verified by MS analysis provided by Bioneer (Daejeon, Republic of Korea). Relative expression ratios were calculated as normalized ratios to MmRpl27 (Mus musculus ribosomal protein L27) gene. Non-template control sample was used for each PCR run to check the primer-dimer formation. The final relative gene expression ratios were calculated as delta-delta Ct values. Information about the genes and the primers is collected in Table S1.
RNA Isolation, Reverse Transcription, and TaqMan Panel
RNAs were isolated from MSC cells with Qiagen RNeasy mini kit (Qiagen) based on the manufacturer's instruction. RNA concentrations were measured by NanoDrop spectrophotometer (NanoDrop). cDNAs were reverse transcribed with TaqMan Reverse Transcription Reagents (Thermo Fisher Scientific) following the manufacturer's instruction. For the TaqMan QRT-PCR panel, cDNA mixture of three parallel samples was used in case of each condition: control, 6, 24, and 72 h exosome treatment applied on MSC cells derived from four mice. Eighty nanogram cDNA and TaqMan Gene Expression Master Mix (Thermo Fisher Scientific) were used for the qPCR experiment. Expression of the examined 44 genes was calculated by ΔΔCt method and normalized to the average Ct values of 4 internal controls (PPIA, 18S RNA, ACTB, and GAPDH).
PD-1 detection by Western blotting
The cells were washed three times in PBS, then were lysed in TENT Buffer (50 mM Tris-HCl, 2 mM EDTA, 150 mM NaCl, 1% TritonX-100, completed with 1x protease inhibitor cocktail (Roche). Protein samples were separated and blotted as described above.
For PD-1 protein level detection, anti-PD-1 (clone: RMO1-14, Biolegend) antibody was used in 1:1,000 dilution. Anti-rat horseradish peroxidase (HRP)-conjugated secondary antibody (R&D Systems) was used in 1:500 dilution.
PD-1 detection by Immunofluorescence microscopy
MSC cultures were fixed with 4% paraformaldehyde for 10 min, and blocked in PBS buffer containing 5% BSA (Sigma-Aldrich) for 60 min at room temperature. After blocking, cells were incubated with primary antibody in DPBS buffer containing 1.2% BSA overnight at 4°C. As primary antibodies, 1:200 anti-PD-1 (clone: RMO1-14, Biolegend,) and 1:200 anti-MLANA (Byorbit) were used. After incubation, cells were washed three times in PBS, and incubated for 1 h at room temperature with 1:100 Alexa Fluor 647-conjugated anti-rat antibody (Jackson ImmunoResearch Laboratories, Baltimore, PA) or 1:500 Alexa Fluor 555-conjugated anti-rabbit antibody (Thermo Fisher Scientific). Tubulin was stained with 1:500 Alexa Fluor 488-conjugated anti-tubulin-α antibody (clone: 10D8, Biolegend). Nucleus staining was performed using DAPI (Life Technologies). Slides were washed three times with DPBS between each step. Images were obtained at 60x magnification using an Olympus confocal laser scanning microscope.
PD-1 detection by STORM super-resolution imaging
All dSTORM super-resolution experiments were performed on a custom-made inverted microscope based on a Nikon Eclipse Ti-E frame. After being conditioned (through spatial filtering via fiber coupling and beam expansion) the applied laser beams were focused into the back focal plane of the microscope objective (Nikon CFI Apo 100x, NA = 1.49), which produced a collimated beam on the sample. The angle of illumination was set via a tilting mirror mounted into a motorized gimbal holder and placed into the conjugate plane of the sample. All the dSTORM images were captured under EPI illumination at an excitation wavelength of 647 nm (Nikon: 647 nm, 300 mW). The laser power, controlled via an acousto-optic tunable filter (AOTF), was set to 4 kW/cm2 on the sample plane. An additional laser (Nichia, Anan, Tokushima, Japan, 405 nm, 60 mW) was used for both reactivation and reference measurements. Images were captured by an Andor iXon3 897 BV EMCCD digital camera (512 × 512 pixels with pixel size of 16 μm). The size of the illuminated region of the sample was matched to the size of the detector, which determined the field of view (FOV = 80 × 80 μm2). Frame stacks for dSTORM super-resolution imaging were typically captured at a reduced image size (crop mode), when only the central 128 × 128 pixel region was selected. A fluorescence filter set (Semrock, Rochester, NY, LF405/488/561/635-A-000) was used to select and separate the excitation and emission lights in the microscope. Additional emission filters (Semrock, BLP01-647R-25) were used in the detection path to further clean the fluorescent light spectrally for the reduction of spectral crosstalk.
During the measurements, the perfect focus system of the microscope was used to keep the sample in focus with a precision of <30 nm. The storage buffer on the sample was replaced with a special switching buffer (20). Typically, 10,000 frames were captured with an exposure time of 30 ms. Reference images with full size FOV were captured at low intensity when the majority of fluorescent molecules were still active and the overall structure of the labeled sample could be visualized.
The captured and stored image stacks were evaluated and analyzed by rainSTORM localization software (21). The individual images of single molecules were fitted with a Gaussian point spread function and their center positions were associated with the position of the fluorescent molecule. Localizations were filtered via their intensity, ellipticity and standard deviation values. Localizations with precisions of <45 nm were only used to form the final image. The estimated mean precision of the accepted localizations was 19 nm. Mechanical drift introduced by either the mechanical movement of the sample or thermal effects was analyzed and reduced by means of a blind drift correction algorithm. Spatial coordinates of the localized molecules were stored and the final super-resolved image was visualized. The multicolor merged images were generated by ImageJ software.
In vivo Experiments
B16F1 melanoma cells were administrated intravenously (1 × 105 cell/100 μl) to 6–8 week old female C57BL/6N mice (Charles River Laboratories). One week later, tumor bearing mice were randomized and divided into 3 groups (n = 10/group). Mice were injected intravenously with control buffer (100 μl), exosome exposed MSCs (1 × 105 cell/100 μl) or exosomes (40 μg/100 μl) 7, 8, 9, 10, 11 days after the tumor injection (Table S2). One week after the first MSC administration, 3 animals/group were euthanized, their lungs were removed, photographed and stored at −80°C for further protein, mRNA and histological analyses. The remained mice were observed for 10 more days. At the end point, the animals were euthanized, and the tumor metastases were investigated not only in their lung, but also in their entire body and removed for histological analysis. Experiments were repeated 3 times. All animal experiments were performed in accordance with national (1998. XXVIII; 40/2013) and European (2010/63/EU) animal ethics guidelines. The experimental protocols were approved by the Animal Experimentation and Ethics Committee of the Biological Research Centre of the Hungarian Academy of Sciences and the Hungarian National Animal Experimentation and Ethics Board (clearance number: XVI./03521/2011 and XV./78/2018).
Tumor coverage of lungs was determined by the analysis of acquired images using the ImageJ software. The area of tumors and the healthy regions was measured and the mean percentage, SD, and p values were calculated in Microsoft Excel.
Quantification of Metastases Associated Vessel Diameters
From native animal lungs 4 μm criostat sections were made on silanized slides, than fixed and retrieved by Fix and Perm A-B solution (Thermo Fisher Scientific, USA) for 20 min. Sections were counterstained by conventional hematoxilin for 30 min than washed in tap water and coverslipped. Sections were digitalized by automatic slide-scanner (3DHistech, Hungary), using the software 3DHISTECH Pannoramic Viewer (3DHistech, Hungary). Strict tumor border was carefully marked then vessel diameters were measured.
Cytokine and Chemokine Array
Lung samples were lysed in NP40 cell lysis buffer (Thermo Fisher Scientific, USA) and protein content was measured by the Pierce BCA Protein Assay kit (Thermo Fisher Scientific, USA). Expression levels of different cytokines in pooled lung specimens were determined using Mouse Cytokine Array Panel A (R&D Systems, Cat. no. ARY006), according to the manufacturer's instructions. Immunoreactive signals were detected by using LI-COR ODYSSEY Fc (Dual-mode imaging system) imager followed by analysis with Image Studio Lite v5.2.
Custom TaqMan Array Panels
QRT-PCR was performed on an ABI Prism 7000 sequence detection system (Applied Biosystems) using specific pre-designed customized 96-well TaqMan Array oncogene panels (Thermo Fisher Scientific) containing 44 specifically selected primers and probes according to extensive literature work. Total RNA was isolated using TRIzol (Thermo Fisher Scientific), DNase treatment was performed according to the manufacturer's protocol, and then total RNA was transcribed into cDNA using High Capacity cDNA Kit (Thermo Fisher Scientific). PCR amplification was performed using TaqMan primers and probes and thermal cycle conditions were set as follows: 2 min at 50°C, 10 min at 95°C, and 40 cycles of 15 s at 95°C and 1 min at 60°C. As internal control transcripts of ACTB (Assay ID: Mm00607939_s1), GAPDH (Assay ID: Mm99999915_g1), and PPIA (Assay ID: Mm02342430_g1) were determined. The amount of the transcripts was normalized to those of the housekeeping genes using the ΔCt method. Finally the results were normalized to the expression of the vehicle control (ΔΔCt method).
For relationship discovery Hierarchical Cluster Analysis was performed by R software. In detail, the “bottom up” agglomerative hierarchical clustering strategy was used and results represented in a tree-based dendrogram [For refer the R software: (22)].
Network Representation by Ingenuity Pathway Analysis
According to literature data, we established a protein network from the in vivo overexpressed genes. A custom graphical representation of this network was generated using the Path Explorer tool of the IPA Path Designer. Genes are represented as red nodes, using various shapes that represent the functional class of the gene product. In a few cases proteins are substituted with the complex, which they are involved in. To identify the potential exosomal factors, which can induce the activation of the network, a list was generated from the exosomal proteins detected by LC-MS/MS and the exosomal miRNAs identified by SOLiD sequencing. The Grow tool in the IPA Path Designer revealed significant interactions between the network and the generated list (the interacting exosomal factors are listed in gray boxes and the types of relationships are indicated in parentheses). During construction in the IPA, the significance level was set to “experimentally observed” data to avoid the representation of predicted, unproven interactions.
All of the data are presented as the mean ± SD or SEM and represent minimum of three independent experiments. Statistical parameters including statistical analysis, statistical significance, and n value are reported in the figure legends. For in vivo experiments n = number of animals. For statistical comparison, we performed two-tailed Student's t-test. A value of p < 0.05 was considered significant [represented as *p < 0.05, not significant (n.s.)].
Isolated Vesicles Show Exosomal Properties
First, we isolated extracellular vesicles from B16F1 mouse melanoma cells. As it was shown by SEM and AFM (Figures 1A,B), the isolated fraction indeed contained exosomes as the particles were cap-shaped, and their size was within the 40–120 nm range. Then, presence of molecules (CD9, CD63, CD81, and HSP70), characteristic for exosomes (4, 17) was assessed by Western blotting (Figure 1C).
Figure 1. Characterization of B16F1 melanoma cell culture-derived exosomes. (A) Scanning electron micrograph of melanoma exosomes. (B) Atomic force microscopy images of exosomes. Left and middle pictures represent the shape and surface topography of vesicles, right graph represents height profile. (C) Western blot analysis of common exosome markers (CD9, CD63, CD81, and HSP70). (D,E) Top 5 molecular and cellular functions identified by Ingenuity Pathway Analysis of exosomal proteins and miRNAs.
Exosomes were then subjected to large-scale analysis to determine their protein and miRNA profiles. Whole proteome analysis (using LC-MS/MS) and bioinformatics tools (used to identify elements of the exosome's proteomics spectrum in the UniProt database, and to compare the findings with entries of the ExoCarta database) revealed that 95 distinct proteins were identified in melanoma exosomes (Table 1). These proteins exhibited 86.3% overlap with molecules listed as characteristic exosome protein markers in ExoCarta. Further, miRNA sequencing (using SOLiD 5500xl technology) identified 168 known miRNA elements (Table 2) which, similarly to the proteomics data, exhibited a large (93.5%) overlap with molecules listed as characteristic exosome miRNA markers in ExoCarta.
To uncover the functional significance of the proteomics and miRNA sequencing data, IPA was employed. This analysis has shown that the identified proteins most probably participated in cellular and molecular processes such as “Cell Death and Survival,” “Cellular Movement,” “Cell-to-Cell Signaling and Interaction,” “Cellular Growth and Proliferation,” and “Cell Morphology” (Prange = 7.53 × 10−15-9.32 × 10−4 significance range) (Figure 1D). Very similarly to these data, functions of the identified miRNAs were suggested to be linked to mechanisms of “Cellular Development,” “Cellular Growth and Proliferation,” “Cellular Movement,” “Cell Cycle,” and “Cell Death and Proliferation” (Prange = 1.25 × 10−12-4.88 × 10−2 significance range) (Figure 1E).
Tumor Exosome Exposure Resulted in Oncogenic Reprogramming of MSCs in vitro
Then, we investigated the effect of exosomes on biological processes (e.g., proliferation, survival, malignant transformation, etc.) of MSCs, which are generally considered as proper in vitro models of tumor stroma (12). For these experiments, MSC cultures were initiated from mouse abdominal adipose tissue (17) and were subjected to melanoma exosome treatment.
First, we assessed whether exosomes were internalized by MSCs. High-throughput microscopy showed that MSCs [labeled green by the DiOC18(3) lipid dye] indeed took up exosomes [labeled red by the DilC18(3) lipid dye] as early as 1–2 h after application (Figure S1). Importantly, after 24 h, the majority of MSCs were loaded by exosomes (Figures 2A,B). Indeed, image analysis and statistical evaluation revealed a 91% internalization efficacy. This suggests that functional alterations demonstrated by this study were due to exosome-induced cell-population, and not individual cell level effects.
Figure 2. Internalized melanoma exosomes induce malignant transformation of the recipient MSCs. (A) Fluorescent images of exosome uptake by MSCs. DiO (green lipid dye)-labeled cells were exposed to Dil (red lipid dye)-labeled exosomes for 24 h and fixed in 4% PFA. Control cells were treated in the same manner without exosomes. Pictures were acquired by an Operetta high content screening system (Perkin Elmer). Lower two images represent an experiment, where nuclei were counterstained with DAPI. Exosomes were detected with a customized version of A-trous wavelet transform, and were highlighted with yellow. (B) Quantitative analysis of exosome uptake. The graph shows the number of detected exosomes (Y axis) in each cell (X axis) in three exosome-exposed (24 h) and three control cell cultures (ctrl). (C) Cell proliferation assay of exosome-exposed MSCs. Cells were plated at 1 × 104 cell/cm2 density in control and exosome-exposed cultures, and 72 h after the exosome treatment, the cell number was determined by manual counting and by an automated cell counter. Both methods showed a significantly increased cell proliferation of exosome-treated cells. Results are presented as mean ± SD (n = 3). (D) Apoptosis analysis of exosome-exposed MSCs by flow cytometry. Exosome-pre-exposed cells were treated by 100 ng/ml TNFα for 24 h, stained with Annexin V-FITC and propidium iodide. They were analyzed by flow cytometry both in the case of the TNFα-treated and the untreated control cell cultures. The graph represents the percentage of early apoptotic, late apoptotic and necrotic cells. Both TNFα-induced total cell death (Σ) and necrosis were significantly lower (Σp = 0.03 and p = 0.02, respectively) in exosome-pretreated cell cultures compared to the corresponding control cells. Results are presented as mean + SD (n = 3). (E) QRT-PCR analysis of Mlana and Mitf in MSC cultures treated by exosomes in every 24 h indicated by arrows in the graph. The expression of both mRNAs increased after exosome exposure, but they showed different kinetics. (F) Fluorescent immunocytochemistry of MLANA in exosome-exposed MSC cultures using a primary rabbit antibody to MLANA and a secondary AlexaFluor555-conjugated antibody to rabbit IgG (red). α-tubulin network of cells was directly labeled by an AlexaFluor488-conjugated antibody (green) and the nuclei were stained with DAPI (blue).
Then, we determined whether the internalized exosomes could induce a melanoma-like malignant transformation of the MSCs. By employing two complementary cell-counting methods, we found that proliferation rate of the MSCs significantly accelerated 72 h after exosome exposure (Figure 2C). By flow cytometry, we also showed that the exosome-treated MSCs exhibited a partial resistance to the cell death-inducing effects of 100 ng/ml tumor necrosis factor-α (TNFα) as the fraction of the dead cells was significantly decreased in these cultures (Figure 2D).
Since the exosomes were isolated from melanoma cells, we were then intrigued to uncover whether the above alterations (which all argue for the malignant transformation of the susceptible cells) also resulted in de novo appearance of melanoma-specific features in the transformed MSCs. To answer this question, expression of the melanoma-specific markers MLANA and MITF were investigated. By QRT-PCR, we found that mRNA transcript levels of both markers elevated markedly in MSCs upon exosome treatment (Figure 2E and Figure S2), albeit the kinetics of elevation of the two molecules were slightly different. Moreover, in good accordance with the mRNA data, immunofluorescence labeling showed that exosome exposure markedly increased the expression of MLANA at the protein level as well (Figure 2F).
Next, we assessed whether the above effects of exosomes inducing a malignant-like transformation of the MSCs was accompanied by a cellular-molecular oncogenic reprogramming of the target cells. Naïve MSCs were exposed to a standardized volume of exosomes for various time intervals (to avoid the experimental fluctuations, we collected pooled samples from multiple independent in vitro experiments). Then, samples were subjected to QRT-PCR analysis using a self-designed panel of 40 oncogenes and tumor suppressor genes which were previously suggested to play a role in melanoma progression (Table 3).
As shown in Figure 3A, gene expression pattern of MSCs exposed to melanoma exosomes, exhibited a clear oncogenic dominance (compared to the non-exposed cells). This was verified by statistical analysis of the mean relative gene expression levels of all molecules investigated. Statistically higher values were obtained in the case of exosome-treated cells (p = 1.9 × 10−5, p = 0.031 and p = 2.3 × 10−8 for the 6, 24, and 72 h time points, respectively).
Figure 3. Exosome re-educated MSCs show oncogene dominance and PD-1 expression. (A) QRT-PCR analysis of 40 tumor-related genes in exosome-exposed MSCs using a self-designed panel. The graph shows the relative expression values for each gene after 6, 24, and 72 h of exosome exposure. The trend line of altered gene expression pattern (indicated by thick black line) shows an increasing tendency over time (mean ± SD). (B) QRT-PCR analysis of PD-1 in MSCs after 2 and 24 h of exosome exposure. The graph represents mean + SEM (n = 3). (C) Representative immunoblot of PD-1 protein expression in control and exosome-exposed MSCs after 24 h of exosome treatment. (D,E) Fluorescent immunocytochemistry of PD-1 in 24 h exosome-exposed MSC cultures, using a primary rat antibody to PD-1 and a secondary AlexaFluor647-conjugated antibody to rat IgG (red). Nuclei were stained with DAPI. (D) α-tubulin network of cells was directly labeled by an AlexaFluor488-conjugated antibody (green). Images were acquired by confocal microscopy. (E) Images were taken by STORM super-resolution microscopy. STORM super-resolution imaging of PD-1 revealed that PD-1 was localized mostly around the nucleus, which were blurred and were not resolvable using diffraction-limited confocal microscopy (D) of the same region.
As mentioned above, Kleffel et al. has recently shown that melanoma cell subpopulations which overexpress PD-1, quite intriguingly exhibit remarkably increased invasiveness and aggressive growth properties (15). However, the authors did not define the factor(s) which induced the above PD-1 overexpression. Since the above finding strongly suggested an “MSC re-education” capacity of melanoma exosomes to induce malignant-like behavior, we next assessed the expression of PD-1 in MSC cultures.
As expected, only insignificant PD-1 expression (both at the mRNA and protein levels) could be identified in control, non-treated MSCs. In contrast, a marked, significant, and time-dependent elevation of PD-1 expression was detected upon exosome treatment by QRT-PCR (Figure 3B), Western blot and an immunocytochemical analysis (Figures 3C,D). Further, by employing super-resolution microscopy, we were able to identify a dramatic upregulation of PD-1 at the single molecular level in exosome-treated MSCs (Figure 3E).
Importantly, since proteomics analysis did not identify the presence of PD-1 in exosomes, these data suggest that the high PD-1 protein content in exosome-exposed MSCs was a result of de novo induction and not of exosome-mediated molecular transfer. Our findings therefore suggest that melanoma exosome-mediated “re-education” of the cells resulted in a novel MSC population which could be identified as MSCPD-1+.
B16F1 Exosomes Augment in vivo Tumorigenesis and Tumor Progression
After presenting evidence for the in vitro tumorigenic induction potential of exosomes on cultured MSCs, we hypothesized that this phenomenon could be identified in vivo as well. To probe this assumption, we employed the well-known animal model, routinely used in our laboratories (46), in which tumors, developed mostly in the lungs, are induced in mice by intravenous administration of mouse B16F1 melanoma cells (to the tail vein). Then, tumor-bearing mice received buffer, or exosomes isolated from the same B16F1 melanoma cells, or exosome-induced MSCPD-1+ cells.
Notably, the exosome-related groups (i.e., exosome, MSCPD-1+) were characterized by a markedly increased size of tumor-covered lung tissues (the increase proved to be significant in the MSCPD-1+ group) (Figure 4A). Of further importance, we also found that in both exosome-related groups, the numbers of distant metastases were significantly elevated when compared to the control (Figure 4B), and diameter of metastases associated blood vessels significantly increased in exosome and MSCPD-1+ injected groups of tumor bearing mice (Figure S3). Using of cytokine and chemokine arrays, tumor supportive cytokine and chemokine levels were elevated in exosome related groups (Figure S4). As we have seen during our previous studies (46), these metastases were mostly localized to the ovaries and kidneys (and very rarely to the lymph nodes) of control tumor-bearing animals. However, besides these sites, the presence of exosomes resulted in frequent metastases in the lymph nodes and, as a new location, in the liver. Interestingly, in MSCPD-1+ treated mice, exosome-transformed MSCs could be identified in the para-aortic lymph nodes by FISH (Figure 4C) verifying the successful in vivo adherence of MSCPD-1+ cells.
Figure 4. Melanoma exosomes promote tumor progression and metastasis formation in vivo. (A) Tumor coverage of lungs on day 15 in different animal groups. The graph represents mean + SD (n = 3). (B) Number of distant metastases on day 25. One dot represents one animal in each group, red lines show the average number of metastases per animal. (C) FISH analysis of a paraortic lymph node metastases, which showed the presence of exosome-exposed MSC. Y chromosome (red dot) of the male mouse-derived MSC was detected in the metastases of a female mouse. (D) Heatmap and cluster analysis of gene expression pattern, which show protooncogenic dominance in exosome or MSCPD-1+ groups. Robust hierarchical clustering based on fold changes in the gene expression data between selected groups divide the treated groups into several upper classes. Groups treated by exosomes or MSCPD-1+ were compared to untreated controls. Exosomes MSCPD-1+ treatments were closely related. (E) QRT-PCR analysis of PD-1 in the lung samples using TaqMan probes (n = 3). (F) The Venn diagram shows possible relations between a finite collections of different sets of genes measured during gene expression profiling (D).
Lung tissues of the different groups were then subjected to in-depth expressional profiling 14 days after injection of exosomes or MSCPD-1+ cells. Namely, QRT-PCR analysis was performed using a self-designed panel of 40 genes; besides housekeeping genes, we assessed expression of (i) proto-oncogenes; (ii) genes reportedly involved in malignant transformation; and/or (iii) genes found to be involved in melanoma development and progression.
Hierarchical cluster analysis (HCA) of the gene expression patterns clearly showed a robust proto-oncogenic dominance in lung samples of the exosome-related groups when compared to the control tissues. Indeed, we identified significant de novo induction of 26 and 23 genes: Alcam1, Eng, Flot2, Itga4, Itga6, Kit, Pecam1, Prom1, Thy, Cdc42, Tiam1, Bcl2, Bax, Casp9, Casp8, Ets1, Hgf, Jak2, Met, Myb, Map2k1, Map2k2, Mapk1 in both group and Elk1, Rb1, Itga2 in the exosome group. An additional six and three genes were identified with at least 10-fold increase (Cd44, Itgb1, Muc1, Pik3ca, Akt1, Rac1 in the exosome group and Cd44, Itgb1, Rac1 in the MSCPD-1+ group). Furthermore, six and eight genes showed at least two-fold mRNA level elevations (Bmi1, Cdk4, Stat3, Kitl, Raf1, Mtor and Bmi1, Muc1, Cdk4, Stat3, Kitl, Pik3ca, Mtor, Akt1, respectively). Figure 4D presents a heatmap of indicative data of overexpressed genes induced by exosome exposure.
Furthermore, we showed that, besides the above genes, expression of PD-1 was also significantly increased in both exosome-related groups (Figure 4E). Notably, although mRNA transcript level of PD-1 was close to two-fold in MSCPD-1+ lung tissues in comparison to samples in the exosome group, the difference was not significant (most probably due to the large inter-animal variability and standard error).
We constructed a Venn-diagram (Figure 4F) to show all possible logical connections between the various gene expression patterns presented in Figure 4D. Importantly, according to cluster analysis, we could not identify a single gene which was missing from the exosome-related groups in comparison to the control tumor-bearing mice. In other words, whereas these exosome-related groups did exhibit gene expression profiles which were characteristic only for them, such individual profiles could not be detected in the control group.
Specifically, the following gene expression patterns were defined:
• three genes (Elk1, Rb1, Igta2) were exclusively induced only in the exosome-treated group;
• twenty-three genes (Prom1, Tiam1, Bcl2, Bax, Casp9, Hgf, Jak2, Met, Mapk2k2, Alcam1, Eng, Flot2, Itga4, Itga6, Kit, Pecam1, Thy, Cdc42, Casp8, Ets1, Myb, Map2k1, Mapk1) were found to be upregulated both in the exosome and MSCPD-1+ cell treated groups;
• fourteen genes (Bmi1, Cd44, Itgb1, Muc1, Cdk4, Nras, Sat3, Kitl, Raf1, Mtor, Akt1, Mapk3, Pik3ca, Rac1) were found to be upregulated in all three groups.
Finally, it should be noted that the dramatic gene expression alterations seen in the exosome-related groups were exclusively due to the presence of the exosomes as the “MSCPD-1+ only” cluster contained no genes (Figure 4F).
Melanoma-Derived Exosomes Promote Tumorigenic and Cell Survival Signaling Pathway(s)
Figure 5 shows the map of interaction pathways based on overexpressed molecules (red symbols) in in vivo experiments. We detected overexpressed elements of three main pathways which do participate in tumor progression and metastasis formation (47). Using the IPA Path Designer Grow tool, we generated an interaction map which contains proteins (or their established complexes) encoded by the overexpressed genes, and exosomal miRNAs and proteins (gray boxes) which were previously shown to control or affect the marked signaling molecules.
Figure 5. Integrated associative network of the in vivo overexpressed genes supplemented with the interacting exosomal factors. The network of overexpressed genes (red symbols) was conceived by us based on literature data. The relationships between molecules were supported also by the IPA knowledge base. Network visualization was performed using the Path Explorer tool of the IPA Path Designer. The exosomal proteins and miRNAs (gray boxes) were connected to elements of the network by the Grow tool of the IPA Path Designer based on experimental data of the IPA knowledge base. Activation of the established network by exosomal components may support the survival, migration, growth, and metastasis of tumor cells.
To demonstrate these results and to explore the potential causative exosomal factors, we again performed an IPA. The analysis allowed further refinement of the underlying molecular processes and pathways involved in the summarized effect of this gene expression profile. The previously described 40-gene QRT-PCR panel was enriched in genes related to cellular movement and migration, cell survival and connective tissue development and function, involving tumor cells. Additional processes included growth and proliferation of tumor cells as well as the PD-1:PD-L1 interaction (and its consequences). Focusing on PD-1 and mTOR pathway, a network was built from the overexpressed genes, and their complexes known to be related to tumor progression using the Path Explorer tool in the IPA Path Designer. Then a list was generated that contained the exosomal proteins detected by MS and the exosomal miRNAs identified by SOLiD sequencing. The Grow tool in the IPA Path Designer revealed significant interactions between the exosome-induced expression network and the molecular content of these vesicles. The resulting interaction network shows both direct and indirect interactions, but it contains only experimentally observed relationships based on the Ingenuity Knowledge Base. This network demonstrates that 61 exosomal molecules (Table 4) may affect tumor progression through pathways controlled by key components including MET, Ras, RAF1, Mek, ERK1/2, MITF, BCL2, PI3K, Akt, mTOR, PD-1, KIT, JAK STAT3, or ETS1.
Table 4. List of exosomal proteins and miRNAs connected to the integrated associative network by IPA Path Designer.
Taken together, these findings demonstrate that interaction between exosomes and MSCs induces a tumor-like phenotype with PD-1 overexpression of naïve MSC in vitro and a fast tumor progression in vivo.
Considering the fact that metastatic complications are responsible for almost 90% of cancer-driven mortality (48), a deep understanding of tumor metastatic processes is one of the most important challenges in both fundamental and applied research directions in tumor biology. Whilst certain organotropic aspects of tumor cells were addressed by Paget's “seed-and-soil” hypothesis in the early nineteenth century (49), the detailed description of the spatiotemporal characteristics of the metastatic processes requires an extended model. MSCs, as multipotent stromal cells, are important players in dynamic tumor microenvironment (12). MSCs could stimulate angiogenesis via paracrine signaling and participate in metastasis formation to promote the epithelial-mesenchymal transition (12). Consequently, in the last decade, numerous studies suggested the importance of soluble and vesicle-like components of the tumor microenvironment in epithelial-mesenchymal transition as crucial factors in the metastatic establishment (50). Among these, recent findings of intercellular communication mechanisms unquestionably highlighted the crucial role of EVs (e.g., exosomes) serving as “information packages” within the coordinated cell-to-cell signaling (51).
In the current study, we demonstrate that MSCs, widely abundant in solid tumors as well as in healthy tissues, undergo a marked re-education process upon communication with metastatic cancer cells via exosome-mediated information transfer. This transformation process results in characteristic response-patterns corresponding to a given cancerous cell lines (12).
First, we intended to precisely define the EV population which we later employed via in vitro and in vivo experimental approaches. Besides applying the experimental requirements of the International Society for Extracellular Vesicles (ISEV) (i.e., determination of shape, size distribution and protein markers) (52), we thoroughly defined the miRNA and protein content of melanoma exosomes. These systematic assays concluded that the isolated vesicles contained factors, which may regulate signaling pathways related to cell death and survival, cellular movement, development and cell proliferation. Given the appropriate uptake of these vesicles by the targeted, multipotent recipient cells, it is highly probable that most signaling pathways mentioned above were significantly affected. In accordance with previous results in the literature, we showed successful internalization of exosomes by recipient MSCs. As expected, following 72 h of exposure, MSCs displayed accelerated proliferation presumably initiated by the transforming exosome cargo. Moreover, melanoma exosomes were able to hinder the induction of the cell death response commonly triggered by TNFα. TNFα is a frequently identified, soluble, tumor-suppressing factor in cancerous microenvironments (53). This inhibition strongly suggested that exosomes initialize and sustain malignant transformation of MSCs.
Further, we described numerous, exosome-induced phenotypic manifestation of the above malignant re-education process in the recipient cells. Apart from identifying tumor markers, the transformation of exposed cells is characterized by markedly modified gene expression patterns. As discussed previously, specific gene expression alterations are often associated with specific steps of cancerous transformation. Malignant cell lines possess altered expression of genes responsible for apoptotic regulation (Bax, Bcl2, caspases), regulating metastasis (integrins, ITGA2-6, ITGB1, KIT) or various growth factors (ETS1, HGF, MET, etc.). Importantly, we successfully identified two melanoma-specific markers, i.e., MITF (Microphthalmia-associated transcription factor) and MLANA (Melanoma antigen recognized by T-cells) in the exosome-exposed MSCs. MLANA and its regulator MITF play a fundamental role in melanocyte development, tumor progression and are overexpressed in melanoma cells (54). The presence of these markers in exosome-exposed MSCs shows effective transfer of the encoded tumorigenic information, and may have a pivotal role in disease progression. Besides, we also identified a characteristic set of genes, belonging mostly to proto-oncogenes, whose expression was significantly elevated upon exosome treatment. Among these molecules, we emphasize the increased expression of CD44 which was previously assigned as cancer stem cell marker and of mTOR, a characteristic indicator of the engagement of the PD-1:PD-L1 pathway (15). Indeed, activation of mTOR signaling is particularly important. As it was previously demonstrated by Kleffel et al., overexpression of PD-1 in exceptionally aggressive melanoma subpopulations facilitated tumor progression by activating the mTOR signaling pathway (actually, this response was highly unexpected compared to that of the classic PD-1:PD-L1 dependent T-cell anergy, a crucial aspect of immunotherapy in the clinics) (15). Quite strikingly, in exosome treated MSCs, we detected a significant induction of PD-1 and mTOR. These novel data strongly support that these vesicles might trigger the formation of an aggressive, melanoma-like subpopulation of re-educated recipient MSCs.
Our findings therefore propose that the re-educated melanoma-like MSCPD-1+ subpopulation, originating from multipotent tissue-derived MSCs, is an autonomous entity. It is characterized by typical properties of cancerous transformation such as hyperproliferation, resistance to apoptosis, expression of melanoma markers, proto-oncogenic gene expression patterns, and PD-1 over-expression. Indeed, those melanoma exosomes that induced in vitro generation of the melanoma-like MSCPD-1+ cells, also facilitated in vivo metastasis formation in the lungs or other distant organs of treated animals. Actually, our findings that melanoma exosomes or exosome re-educated MSCPD-1+ promote tumor progression, are in good accordance with previous finding of Kleffel et al. describing that shRNA or antibody mediated inhibition of PD-1 signaling inhibits metastasis formation and tumor progression in experimental animals (15). Assessment of the expression patterns of genes related to tumor progression in the exposed animal tissues revealed numerous cases of exosome-induced overexpression; this supports our hypothesis that melanoma exosomes should be considered as oncosomes (6).
Our results also demonstrated that melanoma exosomes also generated a characteristic signaling pattern in MSCPD-1+ recipient population. As it was shown in numerous studies PI3K/Akt, Ras/MAPK, and STAT3 pathways are commonly activated in tumor cells (47, 55). In this study, based on our in vitro and in vivo results and literature data, we introduce a novel, comprehensive network of common tumor-related proteins, such as PD-1, MET, RAF1, STAT3, BCL2, or mTOR. The network highlights not only the relationship of these elements, but also contains upstream exosomal regulating factors which may contribute to the activation of tumorigenic signaling, and hence, fast tumor progression. According to our best knowledge, this associative network between overexpressed genes and the potential exosomal inducers, is the first tumor progression signaling pattern which connected experimental response-patterns with experimentally detected exosomal-molecular-patterns. Further, our data also indicate that the above complexity of exosomal communication requires system-level approaches.
Finally, it is important to note that our conclusions are not based on in silico predictions exclusively, but rather on carefully designed and systematically executed in vitro and in vivo experiments, all suggesting that the specific, robust, molecular content of the isolated exosomes can indeed generate a unique intercellular niche responsible for the re-education of neighboring cells via oncogenic transformation. These re-educated melanoma-like MSCPD-1+ cells, in turn, facilitate metastatic disease progression via ignition of a complex series of subsequent events both locally and systematically. Thus, based on our results and recently published, additional evidences about exosomes, there is an urgent demand to supplement and extend the outdated “seed-and-soil” hypothesis with the oncosome-driven re-education process.
Data Availability Statement
All datasets generated for this study are included in the manuscript/Supplementary Files.
All animal experiments were performed in accordance with national (1998. XXVIII; 40/2013) and European (2010/63/EU) animal ethics guidelines. The experimental protocols were approved by the Animal Experimentation and Ethics Committee of the Biological Research Centre of the Hungarian Academy of Sciences and the Hungarian National Animal Experimentation and Ethics Board (clearance number: XVI./03521/2011 and XVI./78/2018).
EG-S: acquisition, analysis and interpretation of data (in vitro and in vivo experiments), drafting or revising the article. MH and GD: acquisition, analysis and interpretation of data, drafting or revising the article. IBN: histology. JM: acquisition, analysis and interpretation of data-in vivo QRT-PCR. ÁZ: acquisition, analysis and interpretation of data (MLANA, MITF). ÉH-G: acquisition, analysis and interpretation of Proteomics data. RK: providing and supervising of stem cell protocols. IN: acquisition, analysis and interpretation of miRNA data. PH: supervising of machine learning. ÁB: high content screening. ÁS: acquisition, analysis and interpretation of machine learning data. MK: acquisition, analysis and interpretation of confocal data. TP: drafting or revising the article, PD-1 QRT-PCR. BB: acquisition, analysis of data (in vitro QRT-PCR). ME: STORM microscopy. ZS: AFM microscopy. ZV: acquisition, analysis and interpretation of apoptosis data. EB: drafting or revising the article. LK: clinical relevance and drafting or revising the article. TB: conception and design and drafting or revising the article. KB: conception and design, acquisition of data, analysis and interpretation of data, drafting or revising the article and founding acquisition.
This work was in part supported by the National Research, Development and Innovation Fund of Hungary, financed under the NKFI-6-K funding scheme (11493 project), GINOP-2.3.2-15-2016-00015, GINOP-2.2.1-15-2017-00052, János Bolyai Research Scholarship of the Hungarian Academy of Sciences ÚNKP-19-4 New National Excellence Program of the Ministry for Innovation and Technology and University of Szeged Open Access Fund. The dSTORM experiments were funded by the Hungarian Brain Research Program (KTIA_13_NAP-A-I/14) and the GINOP-2.3.2-15-2016-00036.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors wish to thank Gábor Braunitzer and Csaba Vizler, for scientific discussion, Annamária Marton for technical issues and adaptation of MSC isolation protocol, Lilla Pintér for technical support.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.02459/full#supplementary-material
3. Braicu C, Tomuleasa C, Monroig P, Cucuianu A, Berindan-Neagoe I, Calin GA. Exosomes as divine messengers: are they the Hermes of modern molecular oncology? Cell Death Differ. (2015) 22:34–45. doi: 10.1038/cdd.2014.130
4. Yáñez-Mó M, Siljander PR, Andreu Z, Zavec AB, Borràs FE, Buzas E, et al. Biological properties of extracellular vesicles and their physiological functions. J Extracell Vesicles. (2015) 4:27066. doi: 10.3402/jev.v4.27066
5. Harmati M, Tarnai Z, Decsi G, Kormondi S, Szegletes Z, Janovak L, et al. Stressors alter intercellular communication and exosome profile of nasopharyngeal carcinoma cells. J Oral Pathol Med. (2017) 46:259–66. doi: 10.1111/jop.12486
7. Minciacchi VR, Freeman MR, Di Vizio D. Extracellular vesicles in cancer: exosomes, microvesicles and the emerging role of large oncosomes. Semin Cell Dev Biol. (2015) 40:41–51. doi: 10.1016/j.semcdb.2015.02.010
10. Suzuki K, Sun R, Origuchi M, Kanehira M, Takahata T, Itoh J, et al. Mesenchymal stromal cells promote tumor growth through the enhancement of neovascularization. Mol Med. (2011) 17:579–87. doi: 10.2119/molmed.2010.00157
12. Baglio SR, Lagerweij T, Pérez-Lanzón M, Ho XD, Léveillé N, Melo SA, et al. Blocking tumor-educated MSC paracrine activity halts osteosarcoma progression. Clin Cancer Res. (2017) 23:3721–33. doi: 10.1158/1078-0432.CCR-16-2726
13. Peinado H, Alečković M, Lavotshkin S, Matei I, Costa-Silva B, Moreno-Bueno G, et al. Melanoma exosomes educate bone marrow progenitor cells toward a pro-metastatic phenotype through MET. Nat Med. (2012) 18:883–91. doi: 10.1038/nm.2753
15. Kleffel S, Posch C, Barthel SR, Mueller H, Schlapbach C, Guenova E, et al. Melanoma cell-intrinsic PD-1 receptor functions promote tumor growth. Cell. (2015) 162:1242–56. doi: 10.1016/j.cell.2015.08.052
16. Bersanelli M, Buti S. From targeting the tumor to targeting the immune system: transversal challenges in oncology with the inhibition of the PD-1/PD-L1 axis. World J Clin Oncol. (2017) 8:37–53. doi: 10.5306/wjco.v8.i1.37
17. Szebeni GJ, Kriston-Pál É, Blazsó P, Katona RL, Novák J, Szabó E, et al. Identification of galectin-1 as a critical factor in function of mouse mesenchymal stromal cell-mediated tumor promotion. PLoS ONE. (2012) 7:411372. doi: 10.1371/journal.pone.0041372
18. Carpenter AE, Jones TR, Lamprecht MR, Clarke C, Kang I, Friman O, et al. CellProfiler: image analysis software for identifying and quantifying cell phenotypes. Genome Biol. (2006) 7:R100. doi: 10.1186/gb-2006-7-10-r100
20. van de Linde S, Löschberger A, Klein T, Heidbreder M, Wolter S, Heilemann M, et al. Direct stochastic optical reconstruction microscopy with standard fluorescent probes. Nat Protoc. (2011) 6:991–1009. doi: 10.1038/nprot.2011.336
23. Donizy P, Zietek M, Halon A, Leskiewicz M, Kozyra C, Matkowski R. Prognostic significance of ALCAM (CD166/MEMD) expression in cutaneous melanoma patients. Diagn Pathol. (2015) 10:86. doi: 10.1186/s13000-015-0331-z
26. Muñoz R, Arias Y, Ferreras JM, Jiménez P, Langa C, Rojo MA, et al. In vitro and in vivo effects of an anti-mouse endoglin (CD105)-immunotoxin on the early stages of mouse B16MEL4A5 melanoma tumours. Cancer Immunol Immunother. (2013) 62:541–51. doi: 10.1007/s00262-012-1357-7
27. Hazarika P, McCarty MF, Prieto VG, George S, Babu D, Koul D, et al. Up-regulation of Flotillin-2 is associated with melanoma progression and modulates expression of the thrombin receptor protease activated receptor 1. Cancer Res. (2004) 64:7361–9. doi: 10.1158/0008-5472.CAN-04-0823
29. Ahmed KBR, Davies AM. New molecular targets for the systemic therapy of melanoma. In: Murph M, editor. Research on Melanoma - A Glimpse into Current Directions and Future Trends. London: Intech (2011). p. 161–80.
31. Madjd Z, Erfani E, Gheytanchi E, Moradi-Lakeh M, Shariftabrizi A, Asadi-Lari M. Expression of CD133 cancer stem cell marker in melanoma: a systematic review and meta-analysis. Int J Biol Markers. (2016) 31:118–25. doi: 10.5301/jbm.5000209
32. Schubert K, Gutknecht D, Köberle M, Anderegg U, Saalbach A. Melanoma cells use Thy-1 (CD90) on endothelial cells for metastasis formation. Am J Pathol. (2013) 182:266–76. doi: 10.1016/j.ajpath.2012.10.003
34. Minard ME, Kim LS, Price JE, Gallick GE. The role of the guanine nucleotide exchange factor Tiam1 in cellular migration, invasion, adhesion and tumor progression. Breast Cancer Res Treat. (2004) 84:21–32. doi: 10.1023/B:BREA.0000018421.31632.e6
35. Raisova M, Hossini AM, Eberle J, Riebeling C, Wieder T, Sturm I, et al. The Bax/Bcl-2 ratio determines the susceptibility of human melanoma cells to CD95/Fas-mediated apoptosis. J Invest Dermatol. (2001) 117:333–40. doi: 10.1046/j.0022-202x.2001.01409.x
37. Chen J, Zhang XD. Nanodelivery of anticancer agents in melanoma: encouraging, but a long way to go. In: Hamblin M, Avci P, Prow T, editors. Nanoscience in Dermatology. Cambridge: Academic Press (2016). p. 189–202. doi: 10.1016/B978-0-12-802926-8.00015-X
38. Hügel R, Muendlein A, Volbeding L, Drexel H, Richtig E, Wehkamp U, et al. Serum levels of hepatocyte growth factor as a potential tumor marker in patients with malignant melanoma. Melanoma Res. (2016) 26:354–60. doi: 10.1097/CMR.0000000000000269
42. Cao HH, Chu JH, Kwan HY, Su T, Yu H, Cheng CY, et al. Inhibition of the STAT3 signaling pathway contributes to apigenin-mediated anti-metastatic effect in melanoma. Sci Rep. (2016) 6:21731. doi: 10.1038/srep21731
43. Roesch A, Becker B, Meyer S, Hafner C, Wild PJ, Landthaler M, et al. Overexpression and hyperphosphorylation of retinoblastoma protein in the progression of malignant melanoma. Mod Pathol. (2005) 18:565–72. doi: 10.1038/modpathol.3800324
46. Buzás K, Marton A, Vizler C, Gyukity-Sebestyén E, Harmati M, Nagy K, et al. Bacterial sepsis increases survival in metastatic melanoma: Chlamydophila pneumoniae induces macrophage polarization and tumor regression. J Invest Dermatol. (2016) 136:862–5. doi: 10.1016/j.jid.2015.12.032
50. Gopal SK, Greening DW, Rai A, Chen M, Xu R, Shafiq A, et al. Extracellular vesicles: their role in cancer biology and epithelial-mesenchymal transition. Biochem J. (2017) 474:21–45. doi: 10.1042/BCJ20160006
52. Lötvall J, Hill AF, Hochberg F, Buzás E, Di Vizio D, Gardiner C, et al. Minimal experimental requirements for definition of extracellular vesicles and their functions: a position statement from the International Society for Extracellular Vesicles. J Extracell Vesicles. (2014) 3:26913. doi: 10.3402/jev.v3.26913
54. Schepsky A, Bruser K, Gunnarsso GJ, Goodall J, Hallsson JH, Goding CR, et al. The microphthalmia-associated transcription factor Mitf interacts with beta-catenin to determine target gene expression. Mol Cell Biol. (2006) 26:8914–27. doi: 10.1128/MCB.02299-05
Keywords: PD-1, exosome, melanoma/tumor progression, stem cell, reprogramming, signalization pattern, metastasis
Citation: Gyukity-Sebestyén E, Harmati M, Dobra G, Németh IB, Mihály J, Zvara Á, Hunyadi-Gulyás É, Katona R, Nagy I, Horváth P, Bálind Á, Szkalisity Á, Kovács M, Pankotai T, Borsos B, Erdélyi M, Szegletes Z, Veréb ZJ, Buzás EI, Kemény L, Bíró T and Buzás K (2019) Melanoma-Derived Exosomes Induce PD-1 Overexpression and Tumor Progression via Mesenchymal Stem Cell Oncogenic Reprogramming. Front. Immunol. 10:2459. doi: 10.3389/fimmu.2019.02459
Received: 21 March 2019; Accepted: 01 October 2019;
Published: 18 October 2019.
Edited by:Fabrizio Mattei, National Institute of Health (ISS), Italy
Reviewed by:Kawaljit Kaur, University of California, Los Angeles, United States
Fatemeh Momen-Heravi, Columbia University, United States
Copyright © 2019 Gyukity-Sebestyén, Harmati, Dobra, Németh, Mihály, Zvara, Hunyadi-Gulyás, Katona, Nagy, Horváth, Bálind, Szkalisity, Kovács, Pankotai, Borsos, Erdélyi, Szegletes, Veréb, Buzás, Kemény, Bíró and Buzás. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Krisztina Buzás, firstname.lastname@example.org