Transfer of miR-100 and miR-125b increases 3D growth and invasiveness in recipient cancer cells
Abstract
Aim: Extracellular communication via the transfer of vesicles and nanoparticles is now recognized to play an important role in tumor microenvironment interactions. Cancer cells upregulate and secrete abundant levels of miR-100 and miR-125b that can alter gene expression in donor and recipient cells. In this study, we sought to identify targets of miR-100 and miR-125b and conclusively demonstrate that microRNAs (miRNAs) can be functionally transferred from donor to recipient cells.
Methods: To identify targets of miR-100 and miR-125b, we used bioinformatic approaches comparing multiple colorectal cancer (CRC) cell lines, including knockout lines lacking one or both of these miRNAs. We also used spheroid and 3D growth conditions in collagen to test colony growth and invasiveness. We also used Transwell co-culture systems to demonstrate functional miRNA transfer.
Results: From an initial list of 96 potential mRNA targets, we identified and tested 15 targets, with 8 showing significant downregulation in the presence of miR-100 and miR-125b. Among these, cingulin (CGN) and protein tyrosine phosphatase receptor type-R (PTPRR) are downregulated in multiple cancers, consistent with regulation by increased levels of miR-100 and miR-125b. We also show that increased cellular levels of miR-100 and miR-125b enhance 3D growth and invasiveness in CRC and glioblastoma cell lines. Lastly, we demonstrate that extracellular transfer of miR-100 and miR-125b can silence both reporter and endogenous mRNA targets in recipient cells and also increase the invasiveness of recipient spheroid colonies when grown under 3D conditions in type I collagen.
Conclusion:miR-100 and miR-125b target multiple mRNAs that can regulate 3D cell-autonomous growth and invasiveness. By extracellular transfer, miR-100 and miR-125b can also increase colony growth and invasiveness in recipient cells through non-cell-autonomous mechanisms.
Keywords
INTRODUCTION
MicroRNAs (miRNAs) are small (22-23nt) RNAs that base pair with target mRNAs to regulate gene expression by decreasing translation and promoting RNA degradation[1,2]. In cancer, miRNA expression patterns are frequently altered, contributing to all aspects of tumor progression from proliferation to metastasis[3]. We and others have focused on the role of miR-100 and miR-125b in cancer[4-9]. These miRNAs are encoded in the third intron of MIR100HG, a 3kb lncRNA that is transcribed from the large 330 kb MIR100HG locus[10]. Increased expression of miR-100, miR-125b, and MIR100HG can promote cetuximab (anti-EGFR) resistance in colorectal cancer (CRC) cells, drive the progression of pancreatic cancer, and promote the epithelial to mesenchymal transition (EMT)[7,8,11]. These effects are driven in part by regulating the expression of target mRNAs including mTOR[4], negative regulators of Wnt signaling[8], cystic fibrosis transmembrane regulator (CFTR)[9], cingulin (CGN)[9], and checkpoint kinases (CHK)[11]. Beyond these mRNAs, TargetScan predicts a large number of additional potential targets for miR-100 (59 targets) and miR-125b (931 targets) (https://www.targetscan.org/vert_80/). Base pairing between miRNAs and target
Extracellular vesicles (EVs) and nanoparticles participate in cell-cell communication through the transfer of RNA, protein, and lipids[15,16]. Among RNA cargo, the transfer of miRNA is the most well-characterized due to their small size and the ability to experimentally monitor the silencing of target mRNAs in recipient cells using reporter constructs[17]. We previously showed that miR-100 and miR-125b are secreted from CRC cells, with Transwell co-culture experiments supporting the extracellular transfer of these miRNAs to recipient cells[4]. However, since the recipient cells in those transfer experiments expressed endogenous
To identify additional targets of miR-100 and miR-125b and to definitively demonstrate extracellular transfer, we created cell lines with CRISPR/Cas9-mediated knockouts of miR-100, miR-125b, or both. Through extensive bioinformatic analyses of RNAseq data, we identified and tested 15 candidate targets for miR-100 and miR-125b, with the two most statistically significant being CGN and protein tyrosine phosphatase receptor type R (PTPRR). We also show that miR-100 and miR-125b contribute to enhanced 3D growth and invasiveness of CRC and glioblastoma cells. Importantly, we demonstrate the functional transfer of miR-100 and miR-125b from donor to recipient cells, leading to the silencing of both reporter constructs and endogenous genes, as well as increased invasiveness of spheroid cultures grown in type I collagen.
MATERIALS AND METHODS
Cell lines
DLD-1 cells are a heterozygous KRAS colorectal cancer cell line that was originally isolated by Daniel L. Dexter. The Sasazuki lab later generated isogenic KRAS lines (DKO-1 and Dks-8) from the parental DLD-1 cells and generously provided them for our use. These cells are also available from ATCC (CCL-221) and from Horizon Discovery. The LN-18 glioblastoma line was generously provided by the King lab at Vanderbilt and is also available from ATCC (CRL-2610).
Cell culture
The CRC cell line HCA-7 was plated in 3D culture of type I collagen to generate cystic colonies (CC) cells[8]. CC cells were then subjected to iterative selection in 2D (cetuximab resistant) and 3D (cetuximab sensitive) conditions until a resistant line (CC-CR) was generated that is capable of growth in 3D in the presence of cetuximab[8]. The LN-18 glioblastoma line (ATCC CRL-2610) was a kind gift from the lab of Dr. Michael King at Vanderbilt. All cell lines were confirmed to be free of mycoplasma contamination [VenorTMGeM Mycoplasma Detection Kit, polymerase chain reaction (PCR)-based, Sigma Aldrich Cat #MP0025]. Cells were grown in Dulbecco’s Modified Eagle’s Medium (DMEM, Corning) supplemented with 10% fetal bovine serum (FBS, Corning), 1% L-glutamine (200 nM, Gibco), 1% MEM nonessential amino acids (Sigma Aldrich), and 1% penicillin-streptomycin (10,000 U/mL, Gibco) (termed complete DMEM) in 5% CO2 at
Generation of knockout cell lines
We used CRISPR/Cas9 multiplex genome engineering[19] to delete miR-100, miR-125b, or both. Guide RNAs (gRNAs) flanking either miR-100 or miR-125b or both were designed using the Benchling CRISPR gRNA Design Tool (https://www.benchling.com/crispr) [Supplementary Figure 1 and Supplementary Table 1]. Oligonucleotides containing gRNA sequences were first cloned into sgRNA expression vectors (Addgene #53186-53189) and then multiple gRNA cassettes were cloned into lentiviral expression vectors [pLV hUbC-Cas9-T2A-green fluorescent protein (GFP); Addgene, #53190] using Golden Gate cloning technology. Lentiviral stocks were prepared in 293T cells and used to transduce CC-CR cells. After transduction, single GFP+ cells were sorted into 96 cell plates and expanded. PCR amplification of genomic DNA from clonal GFP+ cells was performed using primers flanking the gRNA target sequences to identify deletion clones.
EV isolation
Cells were plated at 9.0 × 106 cells in T-175 flasks (Corning), with each EV collection using at least 3 flasks. Cells were grown to 80% confluency and washed 3 times with 1X phosphate-buffered saline (PBS) before media were replaced with DMEM lacking FBS and grown for 48 h. Culture media were then collected and adherent cells were trypsinized and counted. Conditioned media were then centrifuged at increasing speeds and time: 1,000 × g for 10 min (room temperature), 2,000 × g for 20 min (4 °C), and 10,000 × g for 30 min
Quantitative reverse transcripiotn polymerase chain reaction
Total RNA from both whole cells and EVs was isolated using TRIzol (Life Technologies). Taqman small RNA assays (Life Technologies) were used to quantify miRNA levels. A total of 10 ng of RNA was utilized for each reverse transcription (RT) reaction; 0.67 µL of cDNA was used in each 10 µL quantitative polymerase chain reaction (qPCR) reaction. qPCR was completed in either 96-well or 384-well plates using Biorad CFX96 or CFX384 instruments. Fold changes were calculated as previously described[4].
Transfer assays
For Transwell co-culture assays, cells were plated in Transwell dishes with 0.4 µm filters (Corning, 3460), with 50,000 cells in the donor well and 300,000 cells in the recipient well. Transwell co-cultures were grown for five days, changing media every 72 h. Cells were washed 3 times with 2 mL of 1X PBS; RNA was then collected from either donor or recipient cells.
AGO2 immunoprecipitation
For immunoprecipitation, anti-AGO2 (clone 11A9 antibody; MABE253, Sigma-Aldrich) or control IgG antibodies were used. Ten 150 cm dishes of either CC or CC-CR cells were grown to 95% confluency. Cells were then scraped in ice-cold PBS and centrifuged to collect cell pellets in ice-cold lysis buffer (20 mM Tris HCl pH 7.5, 150 mM KCl, 0.5% NP40, 2 mM EDTA, 1 mM NaF, 0.5 mM DTT). After centrifugation of insoluble cell debris, cell lysates were incubated with Dynabeads. Prior to cell lysis, DynabeadsTM Protein G (ThermoFisher, 10003D) were conjugated with antibodies by incubating overnight with rotation at 4 °C and washed two times with ice-cold lysis buffer. Cell lysates were then added to the conjugated beads and incubated at 4 °C for 6 h. Beads were gently washed with fresh lysis buffer and proteins were eluted. For RNA isolation, beads were incubated with DNase and Proteinase K for 10 min at room temperature. Trizol was then added and incubated at room temperature for 10 min and RNA was isolated.
Small RNA sequencing
RNA libraries were generated using the NEXTFlex Small RNA Library Preparation Kits v3 (Perkin) with the following modifications: (1) 3’- and 5’-adaptors were diluted 1:8; (2) 3’-adaptor ligations were performed overnight in multiple steps: 25 °C for 2 h, 20 °C for 4 h, and 16 °C overnight; (3) following cDNA synthesis and before barcoding PCR, clean up was performed using the No Size Selection protocol (Perkin); and (4) PCR amplification was 20 cycles. Following PCR amplification, individual libraries were size-selected (136-200 bp product) using Pippin Prep (Sage Sciences). Size-selected libraries were quantified using a Qubit Fluorometer. Libraries were checked for quality and sequenced using Illumina short-read technology. Libraries were pooled and paired-end sequencing (PE-150) (equimolar multiplexed libraries) was performed on the NovaSeq6000 platform using the VANTAGE core (Vanderbilt University). Demultiplexing and bioinformatic analyses were performed using the TIGER pipeline[21]. Briefly, Cutadapt (v1.16) was used to trim 3’ adaptors and all reads with < 16 nucleotides (nts) were removed[22]. Quality control on both raw reads and adaptor-trimmed reads was performed using FastQC (v0.11.9) (www.bioinformatics.babraham.ac.uk/projects/fastqc). The adaptor-trimmed reads were mapped to the hg19 genome, with additional rRNA and tRNA reference sequences using Bowtie1 (v1.1.2), allowing only one mismatch[23].
Total RNA sequencing
Bulk RNA sequencing libraries were prepared using Universal RNAseq kits (Tecan). Libraries were cleaned (Zymo), checked for quality using the Agilent bioanalyzer, quantified (Qubit), and pooled based on equimolar concentrations. Pooled libraries were sequenced using Illumina short-read technology based on PE150 on the NovaSeq6000 (Illumina). After sequencing, samples (libraries) were demultiplexed and analyzed. Briefly, adapter sequences were removed using Cutadapt v2.10)[22]. Quality control of both raw reads and adaptor-trimmed reads was performed using FastQC (v0.11.9). After adaptor trimming, reads were aligned to the Gencode GRCh38.p13 genome using STAR (v2.7.8a)[24]. FeatureCounts (v2.0.2)[25] was used to count the number of reads mapped to each gene. Heatmap3[26] was used for cluster analysis and visualization. Differential expression was analyzed by DESeq2 (v1.30.1)[27]. Significant differentially expressed genes were determined with fold change > 2 or < 0.5, and adjusted P values (padj) < 0.05. Genome Ontology and KEGG pathway over-representation analyses were performed using the WebGestaltR package (NULL)[28]. Gene set enrichment analysis (GSEA) was performed using GSEA package (v4.3.2)[29] on database v2022.1.Hs.
Gene expression analysis
The Gene Expression database of Normal and Tumor Tissues (GENT2) was used to compare levels of expression between 72 paired normal and tumor tissues utilizing the GPL570 platform. Statistical significance was determined using two-sample t-tests [GPL570 platform (HG-U133_Plus_2)].
Transfection of luciferase plasmids
Cells were grown to approximately 50%-60% confluency with transfections performed using TransIT®-2020 (Mirus, MIR5404). The Renilla luciferase construct (pClneo-RL, plasmid #115366, Addgene) was transfected at 0.5 μg/mL and the firefly luciferase construct was transfected at 1μg/mL. DNA and TransIT-2020 were incubated for 30 min at room temperature in Opti-MEM (Gibco). Transfection complexes were diluted 1:10 in DMEM and added to cells.
Luciferase reporter assays
3’-UTRs from mRNAs of interest were cloned using RT/PCR into the pMIR-REPORTTM miRNA Expression Reporter Vector System (Firefly Luciferase) using NEBuilder HiFi DNA Assembly (New England BioLabs). PCR amplification primers are listed in Supplementary Table 1. All plasmids were verified by sequencing (Genewiz). Transfected cells were cultured alone or in co-cultures using Transwells with 0.4 mm membranes for 72 h, followed by washing with 1X PBS. Luciferase levels were then detected using the Dual-Glo® Luciferase Assay System (Promega). Changes in firefly luciferase levels were normalized to Renilla luciferase levels and compared to luminescence from the empty (no 3’UTR) vector.
Immunofluorescent staining and imaging of cells
For cell images, cells were plated on coverslips and grown for three days, after which coverslips were washed three times in 1X PBS and fixed in 4% paraformaldehyde for 30 min at room temperature. Samples were then washed and allowed to block for 1 h at room temperature in blocking buffer (1X PBS with 0.001% Triton-X and 0.03% donkey serum). Primary antibody incubation (CGN, ab244406, Abcam) was overnight at 4 °C, followed by three 1X PBS wash steps, and incubation with secondary antibodies (Cy3 anti-Rabbit antibody, Jackson ImmunoResearch and Alexa Fluor 488 conjugated phalloidin, ThermoFisher) for 2 h at room temperature. Samples were washed three times with 1X PBS and mounted on glass slides using Vectashield with 4’,6-diamidino-2-phenylindole (DAPI) and imaged using a Zeiss LSM 880 microscope (Vanderbilt University, CISR). All images were analyzed using ImageJ[30].
For quantification of immunofluorescence in cells, cells on glass coverslips were fixed in neutral-buffered formalin and incubated in a 0.1% Triton X-100 solution (Sigma-Aldrich X100) for 10 min. Cells were blocked for 1 h at room temperature with 10% normal goat serum (Vector Labs S-1000) and 22.5 mg/mL glycine and then incubated overnight at 4 °C in a humidified chamber with anti-CGN (1:500) in 1% bovine serum albumin (Sigma A1470) and 0.1% TWEEN 20 (Sigma-Aldrich P1379). Coverslips were mounted onto slides using ProLong Glass Antifade Mountant with NucBlue (Invitrogen P36981) following secondary antibody incubation with goat anti-rabbit IgG AlexaFluor 594 (1:500, Invitrogen A-11012). A corresponding no primary antibody control was performed for all conditions to confirm specificity. Stained samples were imaged using a Leica DMi8 inverted microscope with Leica DFC9000GT sCMOS fluorescence digital camera. For fluorescence microscopy, a Lumencor mercury-free SOLA light engine was used for the illumination source. The microscope was fitted with DAPI, GFP, Texas Red (TXR), and Y5 filter cubes. Images (8 bit) were captured using LASX imaging software. Images of CGN were captured as Z-stacks with 0.5 μm spacing between slices. All images were analyzed using Fiji software[31]. The area of positive staining was set based on the no primary control images and thresholded accordingly. Z-stacks were merged using maximum intensity projections. For CGN analysis, the area of positive expression was normalized by nuclei number.
Spheroid cultures
Cells (10,000) were plated in Nunclon Sphera low adherent plates (Thermo Scientific) and centrifuged at 1,000 × g for 10 min at room temperature, as previously described[32] in complete DMEM. Cells were grown for three days before transfer into 2 mg/mL type I collagen (Advanced Biomatrix). Once collagen was solidified, complete DMEM was added and changed every two days. Colonies were imaged using Muvicyte Live Cell imager for five days. Images were processed and analyzed using ImageJ[30].
Overexpression of CGN
The CGN coding region was amplified by RT/PCR from RNA isolated from CC cells using primers as shown in Supplementary Table 1. PCR products were cloned into a Sleeping Beauty-based Tet-On vector (pSB829) using NEBuilder HiFi DNA assembly. The resulting construct was confirmed by sequencing. For transfections, 0.5 μg plasmid DNA (pCDNA3.1-SB100X) and 0.5 μg of overexpression vector were incubated with 2 μg of Lipofectamine 2000 (ThermoFisher) in Opti-MEM (Gibco) for 30 min at room temperature. Transfection complexes were diluted in 1:10 in complete DMEM and added to LN-18 cells and incubated for 24 h. Cells were then selected using 300 μg/mL of hygromycin. CGN expression was induced by adding 300 ng/mL of doxycycline (DOX).
Protein collection and western blotting
Proteins were collected using 1X-RIPA buffer (Life Technologies) and quantified using BCA assays (ThermoFisher). 40 μg of the total protein was loaded onto pre-cast SDS gels (4%-20% Mini-PROTEAN TGX 50 μL pre-cast gels; BIO-RAD). Separated proteins were transferred to nylon membranes using the Trans-Blot Turbo Transfer System (BIO-RAD). Membranes were blocked with Intercept Blocking Buffer (IBB) (LI-COR) for 1 h at room temperature. Primary antibodies [CGN, Abcam and glyceraldehyde
RESULTS
To facilitate the identification of mRNA targets for miR-100 and miR-125b and to develop cell lines to test for functional miRNA transfer, we used CRISPR/Cas9 technology to generate deletions of these two miRNAs within the MIR100HG locus [Figure 1A]. A single lentiviral vector expressing multiple gRNAs, Cas9, and GFP was created and lentiviruses were transduced into cetuximab-resistant CC-CR cells[8]. The gRNAs were designed to base pair with regions flanking miR-100, miR-125b, or both [Supplementary Figure 1]. Targeting both miR-100 and miR-125b also results in the deletion of let-7a and the BLID open reading frame, but previous work did not support a role for these genes in cetuximab resistance[8]. Individual GFP+ clones were sorted, targeted deletions were identified by PCR, and verified by sequencing. While CC-CR cells express high levels of miR-100 and miR-125b, the parental, cetuximab-sensitive CC cell line expresses ~50-fold less of both miRNAs[8] [Figure 1B]. As expected, the knockout lines showed undetectable levels of miR-100 or miR-125b or both [Figure 1C].
Figure 1. Expression levels of miR-100 and miR-125b in CC, CC-CR and ΔmiR-100/miR-125b cell lines. (A) One of multiple transcription start sites within the 330 kb MIR100HG locus encodes a 115 kb transcript that is spliced (black exons) to generate a 3 kb lncRNA
Identification of mRNA targets of miR-100 and miR-125b
RNAseq was performed on CC, CC-CR, and knockout cell lines (ΔmiR-100, ΔmiR-125b, and
Figure 2. Identification of targets for miR-100 and miR-125b in CRC cells. (A) Differentially expressed genes were identified by mining RNAseq data comparing CC, CC-CR, ΔmiR-100, ΔmiR-125b, ΔmiR-100/miR-125b cells. Venn diagram showing the overlap between genes upregulated in CC and knockout cell lines, predicted mRNA targets of miR-100 and miR-125b using three different algorithms, and association with AGO2 in CC-CR cells; (B) Heatmaps were generated from the differential RNAseq data to identify potential targets upregulated in CC cells (> 2 fold) and knockout cells (> 1.5 fold) compared to CC-CR cells; (C) Gene Ontology analysis of the targets of miR-100 and miR-125b showed enrichment in genes associated with cell migration and motility. Figure was created using ShinyGO[37]. CRC: Colorectal cancer; CC: cystic colonies.
Silencing of mRNA targets of miR-100 and miR-125b
From the targets generated above, we examined miRNA:mRNA pairing [Supplementary Figure 2] and selected 15 candidate mRNAs to test for the ability of miR-100 (3 candidates), miR-125b (7 candidates), or both (5 candidates) to silence luciferase reporter constructs. We designed three different categories of luciferase vectors in which the firefly luciferase open reading frame was fused to (1) no 3’-UTR (empty vector); (2) 3’-UTRs containing three perfect binding sites for miR-100 or miR-125b (100/125 vector); or (3) 3’-UTRs from candidate targets (experimental vectors) [Figure 3A]. These vectors were transfected into either CC-CR or ΔmiR-100/miR-125b cells and luciferase levels were determined. Bona fide targets of
Figure 3. Silencing of reporter constructs containing target 3’-UTRs for miR-100 and miR-125b. (A) Reporter constructs were generated in which the firefly luciferase open reading frame was fused to either plasmid 3’-UTR sequences (empty vector), 3’-UTR sequences containing three perfect binding sites for miR-100 or miR-125b (100/125 vector), or 3’-UTR sequences from candidate mRNA targets of miR-100 or miR-125b (experimental vector); (B) Quantification of relative luciferase levels after transfection of the indicated constructs into CC-CR cells or ΔmiR-100/miR-125b cells. Relative firefly luciferase levels were calculated by normalizing to a co-transfected internal control vector expressing Renilla luciferase. Data (mean ± SEM, n = 3) were analyzed using two-way ANOVA tests, with * indicating P < 0.05, ** indicating P < 0.01, *** indicating P < 0.001, and **** indicating P < 0.0001. SEM: Scanning electron microscopy; ANOVA: analysis of variance.
CGN and PTPRR are downregulated in CRC
miR-100 and miR-125b display elevated gene expression and target multiple mRNAs in CRC and other cancer types[4,8,9,11,38-40]. CC cells grown in 3D in type I collagen form hollow cysts with a central lumen lined by polarized cells, whereas CC-CR cells form solid, disorganized colonies, indicative of aggressive behavior[8]. We thus sought to test if decreased expression of CGN and PTPRR due to increased expression of miR-100 and miR-125b might be observed in CRC and other cancers. We used the GENT2 database to compare the levels of CGN and PTPRR between normal and cancer tissues [Supplementary Figure 3A and B] and found that CGN is downregulated in 19/35 cancers including CRC (red brackets) while PTPRR is downregulated in 16/35 including CRC. There are likely additional mechanisms that can lead to decreased levels of CGN and PTPRR because not all cancers show increased expression of miR-100[38] or miR-125b[39], but the data are in agreement that decreased levels of CGN and PTPRR correlate with cancer.
To test whether our CRC cell lines match the GENT2 data, we used antibodies to perform western blots and immunostaining of CGN protein levels in CC, CC-CR, and knockout cell lines. Even though we observed silencing of PTPRR reporter constructs [Figure 3], we did not detect reduced protein levels by western blots, so we decided to focus hereafter on CGN. CGN connects tight junctions to the cytoskeleton and localizes to apical junctions joining adjacent epithelial cells[9,41]. In cells that express low to undetectable levels of miR-100 or miR-125b (CC and ΔmiR-100/miR-125b cells), CGN was readily detected by western blots [Figure 4A and B] and localized to cell-cell junctions by immunofluorescence [Figure 4C]. Quantification of fluorescent images demonstrated that cells with low or no miR-100 and miR-125b had increased levels of CGN per cell compared to CC-CR cells with high expression of miR-100 and miR-125b [Figure 4D]. Immunofluorescent signals were more readily visualized in CC cells compared to ΔmiR-100/miR-125b cells, potentially due to the regulation of other factors by miR-100 or miR-125b that also seem to affect overall morphology, as indicated by differences in F-actin staining and changes in CGN localization. In cells expressing high levels of miR-100 and miR-125b (CC-CR cells), CGN detection by both western blots and immunofluorescence was dramatically reduced [Figure 4]. These results are consistent with targeting of CGN by miR-100 and miR-125b.
Figure 4. Decreased Expression of CGN in Cells Expressing High Levels of miR-100 and miR-125b. (A and B) Western blots and quantitation of CGN levels in the indicated cell lines; (C) The indicated cells were fixed and immunostained with DAPI to stain nuclei (blue), Alexa Fluor 488-phalloidin to stain F-actin (green), or antibodies against CGN (magenta). Merged images are shown at the far right. Scale bar indicates 10 mm; (D) Quantification of CGN/cell for the indicated cell lines. For (B), data (mean ± SEM, n = 3) were analyzed using Student’s t test. For (C), data (mean ± SEM, n = 3) were analyzed using one-way ANOVA. For significance, * indicates P < 0.05, ** indicates P < 0.01, *** indicates P < 0.001. CGN: Cingulin; DAPI: 4’,6-diamidino-2-phenylindole; SEM: scanning electron microscopy; ANOVA: analysis of variance.
miR-100 and miR-125b contribute to enhanced 3D growth and invasiveness in CRC cells
Gene Ontology analysis of targets of miR-100 and miR-125b showed enrichment of genes associated with cell migration and motility [Figure 2C]. This is consistent with the downregulation of CGN, which leads to increased invasiveness, metastatic potential, and EMT[9,42-44]. To test the effect of miR-100 and miR-125b and their targeting of CGN on 3D growth, we cultured CC, CC-CR, and ΔmiR-100/miR-125b cells under spheroid growth conditions[32] in low adherence plates for 3 days before transfer to growth in type I collagen. Under these conditions at both day 0 and day 5 in collagen, the morphology of ΔmiR-100/miR-125b colonies was more similar to CC colonies than to colonies derived from parental CC-CR cells
Figure 5. miR-100 and miR-125b contribute to enhanced 3D growth and invasiveness. (A) The indicated cell types were plated in low adherent, spheroid cultures for 3 days before transfer to growth in type I collagen and live imaging over 5 days. Scale bar indicates
Quantification of overall spheroid growth revealed that colonies derived from CC-CR cells showed a significant increase in total spheroid area compared to colonies derived from CC or ΔmiR-100/miR-125b cells [Figure 5B]. We also quantified the number and size of projections into the collagen as a measure of invasiveness by determining the increase in projection surface area after growth in collagen for five days. CC-CR colonies showed a significantly greater increase in projection area compared to CC and ΔmiR-100/miR-125b colonies [Figure 5C]. While invasiveness is likely due to multiple changes, the decreased projections we observe in CC and ΔmiR-100/miR-125b colonies are consistent with increased expression of CGN due to lower levels of miR-100 and miR-125b.
miR-100 and miR-125b contribute to enhanced 3D growth and invasiveness in glioblastoma cells
To test whether targeting of CGN by miR-100 and miR-125b can alter the growth characteristics of non-CRC cells, we used the LN-18 glioblastoma line[45]. To our surprise, these cells express even higher levels of miR-100 and miR-125b than CC-CR cells, consistent with undetectable levels of CGN [Figure 6A-D], consistent with the GENT2 analysis in Supplementary Figure 3A which shows decreased expression in brain cancer. Consistent with our data in LN-18 cells, CGN is not highly expressed in glioma cells within The Cancer Genome Atlas (TCGA). When we plated LN-18 cells under spheroid conditions and transferred the colonies to growth in 3D type I collagen over 5 days, we observed an increase in colony growth [Figure 6E]. However, when we overexpressed CGN levels in LN-18 cells, we detected a significant decrease in colony growth [Figure 6E]. This demonstrates that CGN levels can affect colony growth and that the effect of CGN is not restricted to CRC cells.
Figure 6. Increased CGN expression in glioblastoma contributes to decreased 3D growth. (A and B) Relative levels of miR-100 and miR-125b in CC, CC-CR, and LN-18 (glioblastoma) cell lines. Data represent mean ± SEM with ** indicating P < 0.01 and **** indicating P < 0.0001. Significance was determined using Student’s t test; n = 3; (C) CGN levels were determined using western blots in CC, CC-CR, LN-18, and doxycyclin inducible LN-18 cell lines. Lane loading was monitored with blots against GAPDH; (D) Quantification of blots as in (C) showing fold changes in CGN in the indicated cell lines normalized to GAPDH. Data represent mean ± SEM with n = 3. Significance was determined using Student’s t test with **** indicating P < 0.0001; ns indicating not significant; (E) Increase in spheroid area was measured after transfer to collagen from 0-5 days. Surface area tracing was performed using ImageJ. Significance was determined using two-way ANOVA tests, n = 6. Data represent mean ± SEM with ** indicating P < 0.01, *** indicating P < 0.001, **** indicating P < 0.0001; ns indicating not significant. CGN: Cingulin; CC: cystic colonies; SEM: scanning electron microscopy; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; ANOVA: analysis of variance.
Extracellular transfer of miR-100 and miR-125b
Previous work identified CGN and other targets for miR-100 and miR-125b and showed that cellular changes in these targets can drive migration and invasion of CRC cells[9]. Here, we sought to extend those experiments to focus on the extracellular transfer of miR-100 and miR-125b because we had previously shown that these miRNAs can be secreted from CRC cells[4,46]. However, definitive evidence of extracellular transfer of these miRNAs was lacking because the experiments used recipient cells that express miR-100 and miR-125b which can lead to indirect silencing of reporter constructs and potential false positives[18]. Our knockout cell lines provide an ideal experimental system to conclusively test for extracellular transfer of miR-100 and miR-125b between donor and recipient cells. As donor cells, we used CC-CR cells because they not only express high levels of miR-100 and miR-125b, but they also secrete both miRNAs [Figure 7A]. Thus, we set up Transwell assays with CC-CR donor cells and miR-100/miR-125b recipient cells and co-cultured the cells for 5 days before isolation of RNA from the recipient cells and analysis of miR-100 and miR-125b levels. Compared to control transfer experiments in which ΔmiR-100/miR-125b cells were used as both the donor and recipient cells, we observed dramatic increases in the levels of miR-100 and miR-125b in ΔmiR-100/miR-125b recipient cells (~150- and 60-fold, respectively) when CC-CR cells were used as the donor [Figure 7B].
Figure 7. Functional transfer of miR-100 and miR-125b. (A) RNA was collected from EVs from CC and CC-CR cells and the levels of miR-100 and miR-125b were quantified by RT-qPCR. The left panel shows the fold-change enrichment of these miRNAs in EVs compared to cellular levels. The right panel shows the fold change in EV levels of miR-100 and miR-125b when comparing CC-CR and CC cells; (B) Transwell co-culture experiments were performed with ΔmiR-100/miR-125b recipient cells and the donor cells as indicated on the X-axis. RNA was collected from the recipient cells after 5 days in culture and fold changes in miR-100 and miR-125b levels were determined by RT-qPCR; (C) Schematic of Transwell co-culture experiments to test silencing of CGN and PTPRR luciferase reporter constructs in miR-100/miR125b recipient cells; (D) ΔmiR-100/miR-125b recipient cells were transfected with luciferase reporters and co-cultured with either ΔmiR-100/miR-125b (grey) or CC-CR (black) donor cells. Relative luciferase values (y-axis) were calculated in recipient cells, normalizing as in Figure 3. For (A and B), data (mean ± SEM, n = 3) were analyzed using one-way ANOVA tests while the data in panel (D) (mean ± SEM, n = 3) were analyzed using two-way ANOVA tests. For significance, * indicates P < 0.05, ** indicates P < 0.01, *** indicates P < 0.001, and **** indicates P < 0.0001; ns is not significant. EVs: Extracellular vesicles; CC: cystic colonies; RT-qPCR: reverse transcription-quantitative polymerase chain reaction; CGN: cingulin; PTPRR: protein tyrosine phosphatase receptor type-R; SEM: scanning electron microscopy; ANOVA: analysis of variance.
Knowing that miR-100 and miR-125b can be successfully transferred between cells, we tested whether their transfer from CC-CR cells to ΔmiR-100/miR-125b recipient cells can result in functional miRNA transfer. For this, we used Transwell co-culture assays in which the ΔmiR-100/miR-125b recipient cells were transfected with vectors in which the firefly luciferase open reading frame was fused to the 3’-UTRs from either CGN or PTPRR [Figure 7C]. Compared to controls, we observed a significant reduction in luciferase levels when CC-CR cells were the donor cells and ΔmiR-100/miR-125b were the recipient cells [Figure 7D]. These data demonstrate that miR-100 and miR-125b can be functionally transferred between cells to silence reporter mRNAs. To determine whether miR-100 and miR-125b are transferred as part of EVs, we purified crude small EVs (sEVs) from either CC-CR or ΔmiR-100/miR-125b cells and incubated them directly with
To test whether transfer of these miRNAs can silence endogenous CGN in recipient cells, we used Transwell co-culture experiments with ΔmiR-100/miR-125b recipient cells and either CC-CR or ΔmiR-100/miR-125b donor cells followed by immunofluorescent staining of CGN in recipient cells. With ΔmiR-100/miR-125b donor cells, CGN immunostaining was readily apparent in ΔmiR-100/miR-125b recipient cells [Figure 8A] at levels similar to those observed in Figure 4C. However, when CC-CR cells were used as the donor cells, CGN was undetectable in ΔmiR-100/miR-125b recipient cells [Figure 8A]. These results demonstrate that endogenous CGN can be targeted by functional transfer of miR-100 and miR-125b, leading to decreased endogenous target protein levels.
Figure 8. Transfer of miR-100 and miR-125b decreases CGN immunofluorescence in recipient cells. (A) Transwell co-culture experiments were performed with ΔmiR-100/miR-125b recipient cells and the indicated donor cells. After 5 days, recipient cells were immunostained as in Figure 5C. Scale bar indicates 10 mm; (B) Schematic of Transwell co-culture experiments to test growth changes in spheroids due to transfer of material from donor 2-D cells to recipient spheroids grown in 3D; (C) The increase in spheroid colony area was determined by measuring the surface area of spheroids from day 0 to day 5 after exposure to donor cells. Significance was determined using Student’s t test; n = 6. Data represent mean ± SEM with ** indicating P < 0.01. CGN: Cingulin; SEM: scanning electron microscopy.
As a final assay to test for functional extracellular transfer of miR-100 and miR-125b, we used Transwell experiments in which spheroid cultures were plated in 3D type I collagen and co-cultured with either CC-CR or ΔmiR-100/miR-125b cells. Exposure of spheroid colonies to CC-CR cells which express high levels of miR-100 and miR-125b led to an increase in the area of the resultant colonies compared to co-culture with
DISCUSSION
While many reports have suggested that miRNAs can undergo extracellular transfer between cells[4,49-54], the recipient cells in nearly all of these experiments express the same miRNAs, complicating the interpretation of target mRNA silencing after transfer[18]. Silencing under these conditions could be due to indirect effects including unexpected increased expression of endogenous miRNAs. Our knockout cell lines provide an ideal experimental system to obviate such concerns. Using donor cells that express high levels of miR-100 and miR-125b and recipient cells with targeted deletions of these genes, we identified mRNA targets and showed functional transfer of miR-100 and miR-125b to silence both reporter constructs and endogenous genes in recipient cells.
All cells release a heterogeneous mixture of membrane-bound vesicles and nanoparticles with current work devoted to understanding biogenesis and cargo loading mechanisms[16]. Our experiments conclusively demonstrate extracellular trafficking of miR-100 and miR-125b but whether they are transferred as RNA-protein complexes, nanoparticles (exomeres and supermeres), or membrane-encased small or large vesicles remains to be determined[47,55].
Tumor microenvironment implications
Invasiveness and metastasis are indicative of differing stages of cancer and can also be indicative of patient prognosis[56,57]. EMT is a pivotal stage of metastasis involving the breakdown of cell-cell adhesions[58]. Located on apical regions of adjacent epithelial cells, tight junctions play a crucial role in cell-cell adhesion and tissue integrity. Disruption of tight junctions can allow cells to become invasive and metastatic[59]. CGN serves as a tether between the cytoskeleton and tight junction components, interacting with ZO-1 at apical junctions[41,60]. Loss of CGN in epithelial cells disrupts tight junction formation and cell polarity[42-44]. Coupled with our finding that miR-100 and miR-125b are enriched in EVs from CRC cells and that miR-100 and miR-125b can be transferred between cells, our data support a model whereby cell-cell communication within the tumor microenvironment can alter tight junctions, leading to altered growth and increased metastatic potential in recipient cells.
Disruption of tight junctions can also alter the blood-brain barrier in glioblastoma, referred to as the blood-tumor barrier[61]. In this case, tight junctions and adherens junctions in endothelial cells are disrupted, leading to impaired barrier function, development of EMT, and increased invasiveness[62]. Our finding that glioblastoma cells express high levels of miR-100 and miR-125b and correspondingly low levels of CGN are consistent with the increased invasiveness we observe in 3D colony growth and further support the finding that disruption of tight junctions through decreased CGN levels can facilitate tumor progression through both cell-autonomous and non-cell-autonomous mechanisms.
Our experiments provide a simplified model of cell-cell communication within the tumor microenvironment using Transwell assays with donor and recipient cells. In vivo, such communication is extensive, with both tumor and normal (stromal) cells secreting vesicles and nanoparticles that can engage in two-way exchange. The release of EVs and nanoparticles (EVPs) by normal cells, fibroblasts, endothelial cells, and immune cells can inhibit or restrain cancer cell growth, counteracted by the release of EVPs from tumor cells to evade such effects[63]. EVPs released by tumor cells can also travel to distant organs and generate the metastatic niche which can remodel local gene expression patterns to enhance metastatic cell growth[64]. Both miRNA and protein cargo transfer can mediate this cross-communication, as exemplified by findings that immune modulators such as PD-L1 are carried on tumor EVs that function to suppress immune rejection[65]. Beyond cross-communication, it is also possible that the transfer of EVs between tumor cells themselves can enhance the overall progression of events such as EMT by transfer of miRNAs such as miR-100 and miR-125b between cells that express high levels of these miRNAs (and other cargo) compared to lower expressing tumor cells.
MIR100HG, tumor growth, invasiveness and EMT
Beyond demonstrating miRNA transfer between cells, our knockout lines provided a model to test the effects of 3D growth in the presence and absence of miR-100 and miR-125b. In 3D growth in collagen, CC cells form hollow cysts with a central lumen lined by polarized cells, whereas CC-CR cells grow into the central lumen with disorganized, solid colonies[8]. A prior deletion of exon 4 within the MIR100HG locus supported the idea that MIR100HG contributes to both cetuximab resistance and metastasis[7]. Our precise deletions of miR-100 and miR-125b allowed us to test the effect of these miRNAs on 3D growth and showed that the loss of these miRNAs results in slower overall growth with decreased edge dynamics and fewer invasive projections. Since the expression of MIR100HG is unaffected in our knockout lines, our data agree with previous results that MIR100HG can induce EMT, but also that miR-100 and miR-125b play a similar role through the downregulation of CGN[9]. The data are also consistent with our GO analysis of the targets of miR-100 and miR-125b being enriched in cell motility and migration [Figure 2C].
In summary, our data support findings that increased expression of miR-100 and miR-125b in cancer promotes disruption of tight junctions, activation of EMT, and increased metastatic potential. Given that these miRNAs are abundant in EVs secreted from CRC and glioblastoma cells, the data support both cell-autonomous and non-cell-autonomous roles for miR-100 and miR-125b. Given that miRNAs are abundant in both vesicular and non-vesicular extracellular particles, a limitation of our studies is that we have not yet defined which particle(s) is most efficient at transferring miR-100 and miR-125b and we have not defined the precise mechanisms as to how these particles are taken up by recipient cells.
DECLARATIONS
Acknowledgments
The authors would like to thank all members of the Coffey, Weaver, and Patton labs for advice and suggestions. The authors would also like to thank the Coffey and King labs for the generous gift of the cell lines. We acknowledge the use of The Cancer Genome Atlas (TCGA) and the TCGA Research Network.
Authors’ contributions
Generated all data in the manuscript: Nelson HM, Qu S, Huang L, Chapman SN, Luthcke NL, Schuster SA, Turnbull LA, Shameer M, Guy LL, Stamaris TD, Michell DL, Semler EM, Corn KC, Liu X
Provided scientific advice and suggestions and edited the manuscript: Vickers KC, Liu Q, Franklin JL, Weaver AM, Rafat M, Coffey RJ
Wrote the paper: Nelson HM, Qu S, Franklin JL, Rafat M, Patton JG
Availability of data and materials
All RNAseq data are freely available on the GEO website (GSE259371, GSE259370, GSE259373, and GSE259374). All constructs and cell lines are available upon request.
Financial support and sponsorship
This work was supported by PO1 CA229123 to Coffey RJ, Weaver AM, Vickers KC, Liu Q, and Patton JG and an American Cancer Society Research Scholar grant to Rafat M.
Conflicts of interest
All authors declared that there are no conflicts of interest.
Ethical approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Copyright
© The Author(s) 2024.
Supplementary Materials
REFERENCES
1. He L, Hannon GJ. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet 2004;5:522-31.
2. Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet 2010;11:597-610.
3. Peng Y, Croce CM. The role of MicroRNAs in human cancer. Signal Transduct Target Ther 2016;1:15004.
4. Cha DJ, Franklin JL, Dou Y, et al. KRAS-dependent sorting of miRNA to exosomes. Elife 2015;4:e07197.
5. Qin C, Huang RY, Wang ZX. Potential role of miR-100 in cancer diagnosis, prognosis, and therapy. Tumour Biol 2015;36:1403-9.
7. Liu H, Li D, Sun L, et al. Interaction of lncRNA MIR100HG with hnRNPA2B1 facilitates m6A-dependent stabilization of TCF7L2 mRNA and colorectal cancer progression. Mol Cancer 2022;21:74.
8. Lu Y, Zhao X, Liu Q, et al. lncRNA MIR100HG-derived miR-100 and miR-125b mediate cetuximab resistance via Wnt/β-catenin signaling. Nat Med 2017;23:1331-41.
9. Zhang X, Li T, Han YN, et al. miR-125b promotes colorectal cancer migration and invasion by dual-targeting CFTR and CGN. Cancers 2021;13:5710.
10. Emmrich S, Streltsov A, Schmidt F, Thangapandi VR, Reinhardt D, Klusmann JH. LincRNAs MONC and MIR100HG act as oncogenes in acute megakaryoblastic leukemia. Mol Cancer 2014;13:171.
11. Ottaviani S, Stebbing J, Frampton AE, et al. TGF-β induces miR-100 and miR-125b but blocks let-7a through LIN28B controlling PDAC progression. Nat Commun 2018;9:1845.
13. Li N, Flynt AS, Kim HR, Solnica-Krezel L, Patton JG. Dispatched Homolog 2 is targeted by miR-214 through a combination of three weak microRNA recognition sites. Nucleic Acids Res 2008;36:4277-85.
14. Riolo G, Cantara S, Marzocchi C, Ricci C. miRNA targets: from prediction tools to experimental validation. Methods Protoc 2020;4:1.
15. Couch Y, Buzàs EI, Di Vizio D, et al. A brief history of nearly EV-erything - The rise and rise of extracellular vesicles. J Extracell Vesicles 2021;10:e12144.
16. Dixson AC, Dawson TR, Di Vizio D, Weaver AM. Context-specific regulation of extracellular vesicle biogenesis and cargo selection. Nat Rev Mol Cell Biol 2023;24:454-76.
17. O'Brien K, Breyne K, Ughetto S, Laurent LC, Breakefield XO. RNA delivery by extracellular vesicles in mammalian cells and its applications. Nat Rev Mol Cell Biol 2020;21:585-606.
18. Gruner HN, McManus MT. Examining the evidence for extracellular RNA function in mammals. Nat Rev Genet 2021;22:448-58.
19. Kabadi AM, Ousterout DG, Hilton IB, Gersbach CA. Multiplex CRISPR/Cas9-based genome engineering from a single lentiviral vector. Nucleic Acids Res 2014;42:e147.
20. Jeppesen DK, Fenix AM, Franklin JL, et al. Reassessment of exosome composition. Cell 2019;177:428-45.e18.
21. Allen RM, Zhao S, Ramirez Solano MA, et al. Bioinformatic analysis of endogenous and exogenous small RNAs on lipoproteins. J Extracell Vesicles 2018;7:1506198.
22. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet j 2011;17:10.
24. Dobin A, Davis CA, Schlesinger F, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013;29:15-21.
25. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014;30:923-30.
26. Zhao S, Guo Y, Sheng Q, Shyr Y. Advanced heat map and clustering analysis using heatmap3. Biomed Res Int 2014;2014:986048.
27. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 2014;15:550.
28. Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res 2017;45:W130-7.
29. Subramanian A, Tamayo P, Mootha VK, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 2005;102:15545-50.
30. Aebi M, Hornig H, Padgett RA, Reiser J, Weissmann C. Sequence requirements for splicing of higher eukaryotic nuclear pre-mRNA. Cell 1986;47:555-65.
31. Schindelin J, Arganda-Carreras I, Frise E, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods 2012;9:676-82.
32. Carey SP, Starchenko A, McGregor AL, Reinhart-King CA. Leading malignant cells initiate collective epithelial cell invasion in a three-dimensional heterotypic tumor spheroid model. Clin Exp Metastasis 2013;30:615-30.
33. Müller M, Fazi F, Ciaudo C. Argonaute proteins: from structure to function in development and pathological cell fate determination. Front Cell Dev Biol 2019;7:360.
35. Farazi TA, Juranek SA, Tuschl T. The growing catalog of small RNAs and their association with distinct Argonaute/Piwi family members. Development 2008;135:1201-14.
36. Liu J, Carmell MA, Rivas FV, et al. Argonaute2 is the catalytic engine of mammalian RNAi. Science 2004;305:1437-41.
37. Ge SX, Jung D, Yao R. ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics 2020;36:2628-9.
38. Li C, Gao Y, Zhang K, et al. Multiple roles of microRNA-100 in human cancer and its therapeutic potential. Cell Physiol Biochem 2015;37:2143-59.
39. Wang Y, Zeng G, Jiang Y. The emerging roles of miR-125b in cancers. Cancer Manag Res 2020;12:1079-88.
40. Vallacchi V, Rodolfo M. miR-100 and miR-125b regulate epithelial-mesenchymal transition and drug resistance in tumors. Non-coding RNA Investig 2018;2:57.
41. Citi S, Sabanay H, Jakes R, Geiger B, Kendrick-Jones J. Cingulin, a new peripheral component of tight junctions. Nature 1988;333:272-6.
42. Guillemot L, Citi S. Cingulin regulates claudin-2 expression and cell proliferation through the small GTPase RhoA. Mol Biol Cell 2006;17:3569-77.
43. Mangan AJ, Sietsema DV, Li D, Moore JK, Citi S, Prekeris R. Cingulin and actin mediate midbody-dependent apical lumen formation during polarization of epithelial cells. Nat Commun 2016;7:12426.
44. Schossleitner K, Rauscher S, Gröger M, et al. Evidence that cingulin regulates endothelial barrier function in vitro and in vivo. Arterioscler Thromb Vasc Biol 2016;36:647-54.
45. Diserens AC, de Tribolet N, Martin-Achard A, Gaide AC, Schnegg JF, Carrel S. Characterization of an established human malignant glioma cell line: LN-18. Acta Neuropathol 1981;53:21-8.
46. Abner JJ, Franklin JL, Clement MA, et al. Depletion of METTL3 alters cellular and extracellular levels of miRNAs containing m6A consensus sequences. Heliyon 2021;7:e08519.
47. Zhang Q, Higginbotham JN, Jeppesen DK, et al. Transfer of functional cargo in exomeres. Cell Rep 2019;27:940-54.e6.
48. Zhang Q, Jeppesen DK, Higginbotham JN, et al. Supermeres are functional extracellular nanoparticles replete with disease biomarkers and therapeutic targets. Nat Cell Biol 2021;23:1240-54.
49. Mittelbrunn M, Gutiérrez-Vázquez C, Villarroya-Beltri C, et al. Unidirectional transfer of microRNA-loaded exosomes from T cells to antigen-presenting cells. Nat Commun 2011;2:282.
50. Valadi H, Ekström K, Bossios A, Sjöstrand M, Lee JJ, Lötvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol 2007;9:654-9.
51. Skog J, Würdinger T, van Rijn S, et al. Glioblastoma microvesicles transport RNA and proteins that promote tumour growth and provide diagnostic biomarkers. Nat Cell Biol 2008;10:1470-6.
52. Boelens MC, Wu TJ, Nabet BY, et al. Exosome transfer from stromal to breast cancer cells regulates therapy resistance pathways. Cell 2014;159:499-513.
53. Squadrito ML, Baer C, Burdet F, et al. Endogenous RNAs modulate microRNA sorting to exosomes and transfer to acceptor cells. Cell Rep 2014;8:1432-46.
54. Barman B, Sung BH, Krystofiak E, et al. VAP-A and its binding partner CERT drive biogenesis of RNA-containing extracellular vesicles at ER membrane contact sites. Dev Cell 2022;57:974-94.e8.
55. Jeppesen DK, Zhang Q, Franklin JL, Coffey RJ. Extracellular vesicles and nanoparticles: emerging complexities. Trends Cell Biol 2023;33:667-81.
57. Zeineddine FA, Zeineddine MA, Yousef A, et al. Survival improvement for patients with metastatic colorectal cancer over twenty years. NPJ Precis Oncol 2023;7:16.
58. Pretzsch E, Bösch F, Neumann J, et al. Mechanisms of metastasis in colorectal cancer and metastatic organotropism: hematogenous versus peritoneal spread. J Oncol 2019;2019:7407190.
59. Martin TA. The role of tight junctions in cancer metastasis. Semin Cell Dev Biol 2014;36:224-31.
60. Zihni C, Mills C, Matter K, Balda MS. Tight junctions: from simple barriers to multifunctional molecular gates. Nat Rev Mol Cell Biol 2016;17:564-80.
61. Steeg PS. The blood-tumour barrier in cancer biology and therapy. Nat Rev Clin Oncol 2021;18:696-714.
62. Iwadate Y. Epithelial-mesenchymal transition in glioblastoma progression. Oncol Lett 2016;11:1615-20.
63. Reale A, Khong T, Spencer A. Extracellular vesicles and their roles in the tumor immune microenvironment. J Clin Med 2022;11:6892.
64. Wortzel I, Dror S, Kenific CM, Lyden D. Exosome-mediated metastasis: communication from a distance. Dev Cell 2019;49:347-60.
Cite This Article
How to Cite
Nelson, H. M.; Qu S.; Huang L.; Shameer M.; Corn K. C.; Chapman S. N.; Luthcke N. L.; Schuster S. A.; Stamaris T. D.; Turnbull L. A.; Guy L. L.; Liu X.; Michell D. L.; Semler E. M.; Vickers K. C.; Liu Q.; Franklin J. L.; Weaver A. M.; Rafat M.; Coffey R. J.; Patton J. G. Transfer of miR-100 and miR-125b increases 3D growth and invasiveness in recipient cancer cells. Extracell. Vesicles. Circ. Nucleic. Acids. 2024, 5, 497-516. http://dx.doi.org/10.20517/evcna.2024.43
Download Citation
Export Citation File:
Type of Import
Tips on Downloading Citation
Citation Manager File Format
Type of Import
Direct Import: When the Direct Import option is selected (the default state), a dialogue box will give you the option to Save or Open the downloaded citation data. Choosing Open will either launch your citation manager or give you a choice of applications with which to use the metadata. The Save option saves the file locally for later use.
Indirect Import: When the Indirect Import option is selected, the metadata is displayed and may be copied and pasted as needed.
Comments
Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.