- Methodology article
- Open Access
- Published:
A next generation sequencing based approach to identify extracellular vesicle mediated mRNA transfers between cells
BMC Genomics volume 18, Article number: 987 (2017)
Abstract
Background
Exosomes and other extracellular vesicles (EVs) have emerged as an important mechanism of cell-to-cell communication. However, previous studies either did not fully resolve what genetic materials were shuttled by exosomes or only focused on a specific set of miRNAs and mRNAs. A more systematic method is required to identify the genetic materials that are potentially transferred during cell-to-cell communication through EVs in an unbiased manner.
Results
In this work, we present a novel next generation of sequencing (NGS) based approach to identify EV mediated mRNA exchanges between co-cultured adipocyte and macrophage cells. We performed molecular and genomic profiling and jointly considered data from RNA sequencing (RNA-seq) and genotyping to track the “sequence varying mRNAs” transferred between cells. We identified 8 mRNAs being transferred from macrophages to adipocytes and 21 mRNAs being transferred in the opposite direction. These mRNAs represented biological functions including extracellular matrix, cell adhesion, glycoprotein, and signal peptides.
Conclusions
Our study sheds new light on EV mediated RNA communications between adipocyte and macrophage cells, which may play a significant role in developing insulin resistance in diabetic patients. This work establishes a new method that is applicable to examining genetic material exchanges in many cellular systems and has the potential to be extended to in vivo studies as well.
Background
Cell-to-cell communication plays a key role in maintaining the integrity of multicellular systems. The most studied cell-to-cell communication mechanisms include chemical or hormone-mediated signaling and direct cell-to-cell contacts. In late 1990’s, exchange of cellular information was also demonstrated to occur via release of intracellular contents packaged in lipid bilayer vesicles called extracellular vesicles (EVs) or exosomes as reviewed by Colombo et al. [1]. Various types of EVs exist and are generally classified according to their sub-cellular origins. Microvesicles (MVs) are EVs formed and released by budding from cells plasma membrane and display a diverse range of sizes (100–1000 nm in diameter). Exosomes, in contrast, are of endosomal origin and are released as a consequence of multivesicular endosomes fusing with the plasma membrane and are generally small in size (30–150 nm). Importantly, EV secretion appears to be conserved throughout evolution and is a characteristic of most cell types including adipocytes, macrophages, hematopoietic, neuronal, fibroblastic, and various tumour cells [2, 3]. Furthermore, these EVs have been found to contain various types of cargo, such as mRNA, microRNA, proteins, lipids, and DNA [4,5,6], representing diverse ways of mediating cell-to-cell communications. As EVs contain surface molecules that can be recognized by recipient cells, these molecules can be readily shuttled from one cell to another and thereby influence the biological state of the recipient cell in multiple ways [1].
With the potential of EVs as critical mediators of cell-to-cell communication, they have become a key focus of research for numerous pathological settings as biomarkers as well as mediators of disease. A role for EVs has been subsequently identified in various disease settings, including diabetes, cardiovascular disease, inflammation and pain, degenerative brain disorders and cancer, to name a few. For example, Deng et al. demonstrated that exosome-like vesicles obtained from obese mice was sufficient to induce insulin resistance when injected into wild-type C57BL/6 mice [2]. Ibrahim et al. pinpointed exosomes secreted by human cardio sphere-derived cells (CDCs) as critical agents of regeneration of injured heart muscle and cardio protection. They found that the injection of exosomes into injured mouse hearts recapitulated the regenerative and functional effects produced by CDC transplantation, whereas inhibition of exosome production by CDCs blocks these benefits [7]. With respect to cancers, exosomal microRNAs and RNA are in fact being used as diagnostic biomarkers for several cancers such as ovarian cancer [8]. Moreover, exosomes have been adopted as carriers to load MHC class I and class II peptides for vaccinating metastatic melanoma patients [9] and to deliver antitumor microRNA let-7a to treat breast cancer in mice [10]. The readers are referred to recent review articles for detailed roles of exosomes in disease processes [11, 12].
Despite the seemingly extensive characterization of EVs to date, however, limitations exist. First, though several studies have identified many genetic materials within EVs, they were not detailed as necessarily being transported and released into recipient cells [13]. Secondly, there is a lack of estimation of the transferring rate, e.g., the proportion of total genetic materials (in the form of mRNA, miRNA, protein, etc.) being transferred from one cell to another. Mittelbrumn et al. demonstrated the existence of antigen-driven unidirectional miRNA transfer from the T cell to the antigen-presenting cells [14]. However, it is not clear if the unidirectional transfer is true between other cell types. For these reasons, a more systematic method is needed to identify the genetic materials that are potentially transferred during cell-to-cell communication through EVs in an unbiased manner.
In this study, we jointly consider data from RNA sequencing and genotype array to systematically discover mRNA exchanges between two co-cultured cell lines of different genetic backgrounds. We relate those exchanges to potential EV mediated transfer by verifying them with mRNA sequencing of purified EVs. Given the precedent for potential crosstalk between macrophages and adipocytes in adipose tissue under obese conditions leading to insulin resistance [15], we chose to apply our novel methodology to co-cultures of human differentiated adipocytes and macrophages.
Methods
Experimental design and data generation
Experimental design
Our main experiment was performed on an in vitro co-culture cellular system, in which two types of human cells, namely, the adipocytes and macrophages were cultured in transwell plates with porous membrane inserts (pore size of 0.4 μm) to prevent them from being mixed together (Corning, Inc. Costar, NY, USA, see Additional file 1). The porous membrane allows small particles (size less than 0.4 μm) to pass through, making EV mediated mRNA exchanges between the two cell lines possible, which has been demonstrated by Garcia et al. [16] and Zheng et al. [17] (Fig. 1a and Additional file 1: Figure S1). The two cell lines were derived from donors with no known familial relationships. Therefore, a large number of single nucleotide polymorphic (SNP) markers between the two cell lines were expected which would allow us to easily distinguish cell origins. With this experimental setting, we planned to address: if there are mRNA exchanges detectable between the two cell lines; what the genes being transferred and their functions are; and finally if any transferred mRNAs are likely mediated by EVs in the media.
The experimental workflow used to identify mRNA transfers between adipocytes and macrophages in a co-culture system: (a) A schema illustrating the experimental design for cell culture and sample collection including (a) adipocyte cells cultured alone, from which two cell pellet samples ADaloneN1 and ADaloneN2 were retrieved for RNA-seq, adipocyte B1 and B2 were technical duplicates of genotyping using Illumina Omni2.5 SNP array, and an exosome extraction ADexosome was prepared for RNA-seq; (b) same as (a) but for macrophage cells; (c) co-culture of two cell lines, from which three adipocyte cell pellet samples ADcoN1-N3, three macrophage cell pellet samples MOcoN1-N3, and an exosome sample ADMOexosome were profiled by RNA-seq. b The analytical pipeline including the pre-processing of (a) genotype data and (b) RNA-seq data, (c) the Bayesian model to call mRNA transfers between two cell lines, and (d) the final output of this pipeline
Data generation
We performed genotype profiling using the Illumina HumanOmni2.5Exome-8 BeadChip for the two cell lines with two technical replicates for each cell line (Adipocyte (AD) B1/B2 and Macrophage (MO) B1/B2). We also performed RNA-seq profiling for both single cultured and co-cultured cells. For single cultured cells, we generated RNA-seq data from two replicate cell pellet samples (ADaloneN1/N2, and MOaloneN1/N2). For the co-culture cells, we obtained the RNA-seq data from triplicate cell pellet samples in the adipocyte layer (ADcoN1-N3) and triplicate cell pellet samples in the macrophage layer (MOcoN1-N3). EVs were isolated using sequential ultracentrifugation and western blotting demonstrated enrichment for CD9 expression, a marker of endosomal origin. Considering the precedence set in the field, where a mixture of exosomes and MV is likely, we are being conservative and referring to what we isolated in the media broadly as EVs also referred to as ‘exosome-like vesicles’. We performed RNA-sequencing on isolated EVs from media of single cultured and co-cultured cells and profiled the data of ADexosome, MOexosome, and ADMOexosome (Fig. 1a). The readers are referred to Additional file 1 for details regarding genotyping and RNA sequencing on both cells and EVs.
Data processing and quality control
NGS data processing
We first conducted a multiple-step quality control (QC). Briefly, we filtered out ribosomal RNAs (rRNAs) using SortMeRNA [18] and trimmed Illumina adaptors using Trimmomatic [19] on raw pair-ended RNA-seq reads (of length 100 bps) for each sample. 11~90 M pairs of reads were obtained after these two QC steps for all samples except ADexosome which generated much fewer reads (528 K pairs reads). These reads were then mapped to reference human genome (hg38) using an annotation file (GENCODE V21) by STAR (2.4.0.1) [20]. The mapping rates were over 90% for all cell samples and over 80% for exosomes (Additional file 1: Table S1). The uniquely mapped reads were further processed by several steps, such as marking duplicates, Split’N’Trim, reassigning mapping quality, base call recalibration, etc., similar to the protocol used in [21] (Fig. 1b). SAMtools mpileup [22] was used to retrieve the base profile at each locus (see Additional file 1 for more details). We disregarded mapped indel, placeholder, and reference skip bases and only considered remaining bases for our analysis.
We adopted TopHat, HTSeq, and edgeR/DEseq protocols [23] to identify differentially expressed genes between single-cultured and co-cultured samples. Cufflinks protocol [24] was used to quantify gene expression, and GATK [21] to call variants in each sample. For cross sample quality control, we applied Principal Component Analysis (PCA) and the results based on the first two PCs for gene expressions were plotted in Additional file 1: Figure S2. We also showed the consistency of genotypes on common SNPs across 10 cell samples and 2 exosome samples in Additional file 1: Figure S3. Both gene expression levels and variants were consistent with sample annotations, suggesting no sample mislabelling occurred. To maximize the power of detecting mRNA transfer, we also merged RNA-seq reads from the same cell lines and denoted them as ADalone, ADco, MOalone, and MOco, respectively.
Analytical procedures for identifying transferred mRNAs
At a SNP locus, five typical scenarios of nucleotide base composition exist as illustrated in Fig. 1b (c). We define a SNP as a “tag” SNP if it satisfies the following two criteria: (1) shows a polymorphism between the two cell lines; and (2) it has a homozygous genotype in the recipient cell. Tag SNPs are informative to us to infer if donor cells have provided some copies of their mRNAs to the recipient cells using RNA-seq data. Two candidate tag SNPs are highlighted by grey background shading. Using the first highlighted case as an example, the genotype at this locus is “CT” for macrophages and “CC” for adipocytes. For adipocytes under co-culture condition, the mapped reads data contain both “C” and “T” bases, in which the “T” bases may come from macrophages during in vitro co-culture. Alternatively, it is possible that the “T”s in the adipocytes may originate from other mechanisms such as sequencing error, mapping error, and RNA editing. We developed a Bayesian model to distinguish mRNA transfers from the alternative mechanisms. Considering adipocytes co-cultured with macrophages, for a given locus, the Bayesian model takes the following as inputs: (1) the genotype information of donor (macrophage) and recipient (adipocyte) cell lines; (2) the counts of the nucleotides at the locus, and (3) base qualities of the reads mapped to the locus.
Specifically as shown in Additional file 1: Figure S4, we denote the mapped RNA-seq read data at a particular genome position from four profiles by A I , A C , M I , and M C , corresponding to adipocyte alone, adipocyte co-cultured, macrophage alone, and macrophage co-cultured cells, respectively. Let G d and G r be the genotypes of donor cells and receptor cell respectively. When examining the sequence data of adipocyte co-cultured with macrophage, G d would be the genotype of macrophage (donor cell) and G r would be the genotype of adipocyte (receptor cell). We assume that the read depth of A C is N and denote t r (0 ≤ t r ≤ N) as the number of reads in A C that were transferred from donor cells. Obviously, a genetic material transfer happens at the position of consideration if t r ≥ 1.
We frame the question as a Bayesian inference. We calculate the Bayes factor r to measure the confidence ratio of observing the data under two hypotheses, i.e., there exists at least one transfer vs. there is no transfer (null hypothesis).
We reject null hypothesis, and claim a positive genetic material transfer if ris greater than a predefined threshold. We assume that the true genotype of the donor and receptor cells are known, i.e., G r and G d are predefined constants (based on genotype array data). We next calculate P[t r = i| A C ] with 0 ≤ i ≤ N, that is, the posterior probability of exactly i nucleotides being transferred given A C .
To calculate P[t r = i| A C , G r , G d ], based on Bayesian rule,
where the denominator P[A C , G r , G d ] is a constant and is often omitted in Bayesian calculation. The P[t r = i, G r , G d ] = P[t r = i]P[G r , G d ] holds as we assume that the number of reads being transferred is independent of the genotype of donor or acceptor cells.
It is of note that we derive G r and G d based on genotype array data, but given we also have the RNA-seq data, an alternative approach is to estimate P[G r ] and P[G d ] by P[G r | A I ] and P[G d | M I ], which could be calculated in the same way as McKenna et al. [25]. This alternative approach can be useful when genotype information is not available.
P[t r = i] is the prior “belief” that i reads were from transfer at the position under consideration. Without much prior knowledge of transfer, we assume a uniform prior, so that, \( P\left[{t}_r=i\right]=\left\{\begin{array}{l}\frac{1}{2}, if\ i=0\\ {}\frac{1}{2N}, if\ i\ne 0\end{array}\right. \). That is to assume equal probabilities of having genetic material transfer and not having a transfer, and further assume that the probability of transferring i nucleotides to be equal for any i with 1 ≤ i ≤ N.
To calculateP[A C | t r = i, G r , G d ], we assume that reads are independent to each other, thus \( P\left[{A}_C|{t}_r=i,{G}_r,{G}_d\right]=\prod \limits_{j=1}^NP\left[{b}_j|{t}_r=i,{G}_r,{G}_d\right] \), where b j is the nucleotide observed for read j. For N reads in A C , if i of them were from transfer, each read has a probability of \( \frac{i}{N} \) for coming from a transfer and a probability of \( \frac{N-i}{N} \) for not from transfer. Thus, \( P\left[{b}_j|{t}_r=i,{G}_r,{G}_d\right]=\frac{i}{N}P\left[{b}_j|{G}_d\right]+\frac{N-i}{N}P\left[{b}_j|{G}_r\right] \). In addition, for those transferred reads, they are from either one of the parental chromosomes in the donor cell. By assuming that the probability of transferring from maternal chromosome is equal to the probability from paternal chromosome, (i.e., we ignore the situation of allele specific expression), we get\( P\left[{b}_j|{G}_d\right]=P\left[{b}_j|\left\{{A}_1,{A}_2\right\}\right]=\frac{1}{2}P\left[{b}_j|{A}_1\right]+\frac{1}{2}P\left[{b}_j|{A}_2\right] \), where the genotype G d = {A 1, A 2} is decomposed into its two alleles. The probability of observing a base given an allele is
Where Q is the phred scaled recalibrated quality score at the base. Similarly, we can calculateP[b j | G r ].
Based on the method described above, we calculate the Bayesian factor for each locus. We call genetic material transfer occurred at a locus if the Bayes factor is greater than a predefined threshold β (β = 20). We performed the analysis on all ~605,000 tag SNPs (Fig. 2) with read depth equal to or greater than 10 in the co-cultured samples. Finally, we evaluated the possibility of mRNA transfers being mediated through exosomes by cross-referencing to the identity of mRNAs in the EVs (Fig. 1b).
Genotyping quality control and comparison between adipocyte and macrophage cell lines: There are two cell lines each with two technical replicates for their genotyping profiling (Adipocyte B1/B2 and Macrophage B1/B2). The number in each eclipse denotes the number of SNPs kept at that step for the corresponding sample. There are 604,540 “tag” SNPs showing polymorphisms between the two cell lines
It is of note that a very small fraction of loci showed inconsistency between the genotyping data and the RNA-seq reads profile. For example, at chr1:145,746,933, the genotype of adipocyte cell line is “GG”, however all the 56 reads and all the 38 reads mapped to this locus are “C”s in ADaloneN1 and ADaloneN2, respectively (Additional file 2: Dataset S1). Such inconsistent loci were likely caused by genotyping errors, alignment errors, errors in the reference genome, etc., and could be problematic for our down-stream analysis. Therefore, we filtered out genotypes if in the single cultured recipient cells, the proportion of reads different from genotype was larger than a predefined threshold γ. With γ = 0.005, 9861 (2.26%) and 14,028 (3.22%) loci were filtered out for adipocytes and macrophages, respectively (Additional file 2: Dataset S1).
Estimating FDR
We calculated the Bayesian factors for triplicate co-cultured samples (i.e., ADcoN1-N3: adipocytes co-cultured with macrophages, and MOcoN1-N3: macrophages co-cultured with adipocytes) and ranked target loci accordingly (loci with larger Bayesian factor rank to the top) in each sample. We then adopted a robust rank aggregation method [26] to identify loci that were ranked consistently better than expected under null hypothesis of uncorrelated inputs and assigned a significance score (p-value) to each locus. Finally, we adjusted the p-values for multiple comparisons to estimate the false discovery rate (FDR).
Results
Identify mRNA transfers potentially mediated by EVs in co-cultured adipocytes and macrophages
We explored mRNA transfers in two directions, i.e., from macrophage to adipocyte and from adipocyte to macrophage, respectively.
From macrophage to adipocyte
Ten thousand ninety-five loci survived the filtering steps in the analytical procedures. We further required the potential transfer loci to (1) have Bayesian factor larger than or equal to 20 in at least one sample (i.e., ADcoN1-N3), and (2) contain at least 2 reads that are potentially transferred from donor cells (e.g., “T” allele for the example in analytical procedures in ADco (the merged data from ADcoN1-N3)). Three hundred twenty-one loci satisfied such requirements (Additional file 3: Dataset S2). Among these loci, we identified 8 (corresponding to 8 unique genes) that are putative mRNA transfers from macrophage to adipocytes with high confidence (FDR < 0.1, Table 1). All these 8 genes were likely mediated by EVs as they were all expressed in both MOexosome and ADMOexosome (Table 1). It is of note that the number of total transcripts in GENCODE V21 is 60,566, among which 13,697 have FPKM larger than 0 and only 1966 larger than 1 in ADMOexosome RNA expression. The Fisher’s exact test indicates that the overlap between inferred transferred genes and those expressed in ADMOexosome are very significant (p-value < 2.2E-16) (Table 2). We also used Integrative Genomics Viewer (IGV) [27] to visually verify the 8 genes and showed the alignment (of MOaloneN1-N2, ADcoN1-N3, and ADaloneN1-N2) at chr19:10,286,547 (ICAM1) in Additional file 1: Figure S5. As shown in Additional file 1: Figure S5, the read counts summarized by IGV were consistent with the read counts we calculated and listed in Table 1.
A close inspection of the individual genes showed that several of them belonged to the extracellular space (GO:0005615; PAPPA, ICAM1, CTSZ and COL6A3) and lysosome (GO:0005764 NPC1 and CTSZ) GO categories. One unifying theme identified from the top three transfers was related to impacting pathways associated with insulin resistance. For example, our top hit ICAM1 (Intercellular Adhesion Molecule 1) from macrophages to adipocytes is an endothelial- and leukocyte-associated transmembrane protein long known for its importance in stabilizing cell-cell interactions and facilitating leukocyte endothelial transmigration [28]. It has also been demonstrated to associate with insulin resistance and diabetic retinopathy in type 2 diabetes (T2D) mellitus [29, 30]. Interestingly, although ICAM-1 is often annotated as a transmembrane protein, two types of extracellular ICAM-1 have also been detected outside of cells or in serum including a soluble form and a membranous form associated with exosomes.
In addition to inflammatory mediators like ICAM-1, factors related to extracellular matrix (ECM) components of the adipose tissue have recently emerged as important mediators in obesity-related pathogenesis. In particular one of the most abundantly expressed collagens in the adipose tissue forming part of the ECM structure is COL6 and its alpha 3 chain, COL6A3 has been associated with adipose tissue inflammation and fibrosis. In collagen VI knockout (KO) mouse for example on an ob/ob background, adipocytes of the knockout mice were larger than wildtype and blood glucose was normalized suggesting that elements of the ECM restrict expansion of adipocyte during obese insults. Relevance has also been observed in obese humans, with elevated levels of collagen VI being detected as well as significant correlations with macrophage infiltration. Observing mRNA for COL6A3 in exosomes from macrophages suggests another level by which the ECM may be influenced by the cells of the adipose tissue depots [31, 32].
Finally, rounding out the top three hits, is Pregnancy-associated plasma protein-A (PAPP-A), a secreted metalloproteinase. PAPP-A cleaves insulin-like growth factor binding proteins (IGFBPs), thereby functioning as a growth-promoting enzyme by releasing bioactive IGF in close proximity to the IGF receptor. As PAPP-A has demonstrated to have fat depot-specific expression in humans and mice, and as the IGF signalling is known to regulate various adipose tissue processes in part through influencing the insulin/insulin receptor signalling axis, exosomal derived PAPP-A mRNA may serve as another mechanism whereby PAPP-A elicits its autocrine/paracrine actions [33].
From adipocyte to macrophage
Thirty-seven thousand eight hundred fifty-three loci survived the filtering steps in the analytical procedures, among which 599 are potentially involved in the mRNA transfer from adipocyte to macrophage (Additional file 4: Dataset S3). There are 21 high confidence loci with FDR < 0.1, among which 17 are likely mediated by EVs as listed in Table 3. The Fisher’s exact test indicates that the overlap between our inferred genes and those expressed in ADMOexosome are also very significant (p-value <2.2E-16) (Table 2). Furthermore, GO ontology analysis of these genes indicated that several of these genes were annotated under the extracellular vesicular exosome (GO:0070062) category (FCER2, UBL3, AHNAK, DSTN, SOD2, HSPG2 and LAMC1). Similarly, we used Integrative Genomics Viewer [27] to visually check the 17 genes and an example at chr6:75,666,836 (SENP6) is shown in Additional file 1: Fig. S6.
A search in PubMed for a role for our top 3 mRNA transfers, namely, PIEZO1, RMRP and LAMC1 indicate their relevance in cellular communication processes.PIEZO1 is an ion channel mediating mechanosensory transduction in mammalian cells [34]; RMRP is a lncRNA with multiple RNA targets [35] and LAMC1 is a member of laminins(subunit gamma 1), which are a family of extracellular matrix glycoproteins [36]. Moreover, LAMC1 is also a candidate gene for diabetic nephropathy [37]. Several other significant mRNA transfers include HSPG2 and SPTLC2. SPTLC2 encodes an enzyme involved in sphingolipid synthesis, which has been shown upon heterozygous deficiency to protect mice from insulin resistance [38]. Finally, the neuroblast differentiation-associated protein AHNAK is highlighted for its important role in the regulation of thermogenesis and lipolysis in WAT via b-adrenergic signalling. AHNAK−/− mice under a high-fat diet had enhanced insulin sensitivity and browning of the WAT depot [39]. Interestingly, AHNAK is a large plasma membrane protein of various functions including plasma membrane support, calcium signalling as well as regulated exocytosis via its key membership of a specialized vesicle, called enlargeosomes. Enlargeosomes are non-secretory, cytoplasmic vesicles competent of regulated exocytosis after rising intracellular calcium and contribute to plasma membrane repair and vesicle shedding [40, 41]. Finding that AHNAK mRNA is transferred from adipocyte to macrophages in exosomes suggests an additional level of complexity to the role of AHNAK in exocytosis.
Although it remains to be ascertained whether these mRNAs become functional proteins, evidence from others suggest that mRNAs transferred via exosomes do become functional proteins [4]. Thus these 25 genes (from both directions) are viable candidates for further studies with respect to communication between adipocytes and macrophages.
Estimate the mRNA transfer rate between the two cell lines
To estimate the transfer rate and test the performance of our model at different transfer rates, we simulated mRNA transfers at various rates and performed analyses based on simulated data.
mRNA transfer simulation
The general work flow of the simulation process is shown in Fig. 3. Using simulated mRNA transfer from macrophages to adipocytes as an example, we randomly selected a predefined fraction of reads from MOalone and merged them with ADalone data. We then ran our previously described pipeline on the merged ADalone data to see how many loci with “transferred” reads could be correctly identified. For each predefined transfer rate, we performed 10 simulation runs (denoted as \( {AD}_{co}^{Si} \) for 1 ≤ i ≤ 10with S indicates simulated data). Similarly, we constructed \( {MO}_{co}^{Si} \) for 1 ≤ i ≤ 10 to study genetic transfer from adipocytes to macrophages. The simulation study has its unique advantage since we know exactly if a base mapped to a locus is from the donor cells or not, this would allow us to assess the performance of our pipeline of detecting mRNA transfer at different transfer rates.
A pipeline to generate simulated data for adipocytes co-cultured with macrophages: (a) Steps used to generate simulated data. b A schematic to illustrate how we sampled reads from one type of cells and merged them with the reads from the other type of cells to simulate a sample with known mRNA transfers
Calling genetic exchange on simulated data and calculating accuracy
We calculated the Bayesian factor at each “target” locus for the simulated samples. By merging the results from these samples and using Bayesian factor as a binary classifier, we created receiver operating characteristic (ROC) curves and calculated the area under curves (AUCs) of our method at 4 different transfer rates, i.e., 0.0001, 0.001, 0.01, and 0.1 (Fig. 4). As shown in Fig. 4, our method performs very well on identifying mRNA transfers at high transfer rates. At transfer rate 0.1, we achieved AUCs of 0.95 for transfer from macrophages to adipocytes and 0.90 for transfer in the opposite direction. The performance of our method decreases as transfer rate decreases. This is not surprising since when transfer rate decreases, the number of “alien” reads at each locus decreases, making signal to noise ratio to drop.
Performance of the Bayesian framework on simulated data with different transfer rates: (a) The ROC curves and AUCs for our method on simulated data from macrophages to adipocytes at 4 different transfer rates: 0.1, 0.01, 0.001, and 0.0001. b The ROC curves and AUCs for our method on simulated data from adipocytes to macrophages
Estimation of mRNA transfer rate in the co-cultured samples
We provided a rough estimation of the transfer rate in our co-culture system by examining the difference in the overall distribution of Bayesian factor statistics obtained from the real co-culture data and the ones from simulated samples. The rationale underlying this analysis is that when the simulated transfer rate is similar to the real transfer rate, the difference between the two distributions should be minimal. It is of note that the numbers of reads are different in simulated samples and real co-cultured samples; this may cause the distribution of Bayesian factor statistics to change even when the underlying transfer rates are the same. To adjust for variation in the sample total read numbers to achieve a fair comparison, we down-sampled the reads data in ADco (MOco) such that it had roughly the same number of reads as the simulated co-cultured samples. We performed 10 runs of down-sampling, and generated 10 reads profiles denoted by \( {AD}_{co}^{Ri} \) (\( {MO}_{co}^{Ri} \)) (1 ≤ i ≤ 10 with R denotes real data), respectively for each transfer rate, i.e., 0.0001, 0.001, 0.01, and 0.1. We then calculated the Bayesian factor at each “target” locus for the down-sampled data. At a specific transfer rate, the overall distribution of Bayesian factor statistics in real co-culture is estimated by merging all the Bayesian factors from \( {AD}_{co}^{Ri} \) (\( {MO}_{co}^{Ri} \)) (1 ≤ i ≤ 10), and similar procedures were performed for simulated data. The Kolmogorov–Smirnov (KS) test was applied to compare the two distributions and we estimated the transfer rate to be the one with the lowest KS statistics. For both transfer directions, the minimal KS-statistics were achieved at transfer rate close to 0.001 (Table 4), indicating that the transfer rate in our co-culture system under quiescent condition was around 0.001.
Differentially expressed genes between cell lines cultured alone and co-cultured
To evaluate the impact of co-culturing to the transcriptome, we used cell lines cultured alone as control, and identified differentially expressed genes (DEGs) in the co-cultured system. We adopted two commonly used methods DESeq and EdgeR [23, 42] and summarized the results in Table 5. As can be seen in Table 5, the DEGs identified by the two methods were consistent even though DESeq inferred more differentially expressed genes than EdgeR. There were 575 and 142 DEGs (FDR ≤ 0.05) identified by both methods for adipocytes and macrophages respectively (Additional file 5: Dataset S4).
We performed function enrichment analysis of the DEGs using David tools (Version 6.7) [43] and the full results were provided in supplementary materials (Additional file 1: Table S2 for adipocytes and Table S3 for macrophages). The DEGs for adipocytes were mostly enriched in type 1 repeats in thrombospondin family, which are multimeric multidomain glycoproteins that function at cell surfaces and in the extracellular matrix milieu. THBS1, which encodes thrombosponin 1 was one of the genes found within this pathway significantly down-regulated in the co-cultured cells relative to the adipocytes cultured alone. THBS1 is interesting with respect to type 2 diabetes as it has been shown to be elevated in the circulation in obese and insulin-resistance individuals [44] and loss of THBS1 in mice, protects them from diet-induced weight gain and adipocyte hypertrophy [45]. Interestingly, the top 1, 3 and 6 genes identified in our transfer study, i.e., ICAM1, PAPPA, and NPC1 are also significantly differentially expressed between ADalone and ADco with FDRs 1.02E-7, 1.36E-7, and 1.68E-5, respectively. ICAM1 expression has been shown to relate with obesity and insulin resistance [46]. NPC1 haploinsufficiency also promotes weight gain and metabolic features associated with insulin resistance [47].
In contrast to the adipose tissue, the DEGs identified by comparing MOco vs. MOalone are mostly enriched for transmembrane proteins, immune response pathways such as leukocyte activation, and EGF-like domains. Most occurrences of the EGF-like domain are found in the extracellular domain of membrane-bound proteins or in proteins known to be secreted. The EGF receptor family is in part involved in Notch signaling, which controls cell-cell communication [48]. It is of note that a lot of studies have shown that ECM modulates epidermal growth factor receptor activation and leukocytes [49, 50]. Thus, the differential expressed genes could be some down-stream effects of the transferred transcripts from adipocytes to macrophages.
Discussion
There is accumulating evidence that exosomes via horizontal transfer (i.e., the movement of genetic material between cells) of genetic information can play a key role in cell-to-cell communications [4,5,6, 51, 52]. Numerous studies have thus focused on providing a comprehensive characterization of the content of EVs and these efforts have led to the creation of databases, such as EVpedia and Vesiclepedia [53, 54], which record molecules (proteins, mRNAs, microRNAs or lipids) observed within these vesicles. However, being identified in exosomes is not necessarily indicating that the RNAs or proteins will be transferred into other cells. Thus in the present study, we complement these efforts by providing a computational approach to identify genetic material that has been transferred between two in vitro co-cultured cells mediated by EVs. Next generation sequencing of cellular genetic material which differed between the co-cultured cell types was used as the finger print to place donor-derived mRNAs at the scene of the recipient’s cellular RNA pool. In comparison to other labelling technologies, using DNA sequence polymorphism as marker has its advantages since they are naturally occurring and introduce no artificial modifications to the biological systems. On the other hand, we are limited to loci with polymorphisms for mRNA transfer detection, and therefore this approach while at genome-scale provides semi-whole genome coverage.
In comparison of our identified 25 genetic transfers (in both directions) to those catalogued in the exosome database, ExoCarta (~1700 distinct human mRNAs across various tissues sources) [55], we identified 7 mRNAs in common (ITGB1, NFATC3, SOD2, FOS, AHNAK, LAMC1, and SPTLC2), all of which are in the direction from adipocytes to macrophages (with p-value 2.45E-3). The overlap is significant despite the diversity seen in exosome cargo depending on the cell type under study as well as the cellular state. Nonetheless this serves to further underlie the importance of characterization of EV contents. While our method has general applicability, as described below, we perceptively chose to apply this to co-cultures of macrophage and adipocyte given the precedence in the literature that there is crosstalk between adipocytes and immune cells in the adipose tissue that contributes to metabolic dysregulation and obesity [56]. The function of the mRNAs transferred range from protein coding to transcriptional regulators and thus the impact on the recipient cell could vary greatly if they are translated into functional units. Although in our study we did not assess this, evidence from other studies suggest the transferred mRNA can be functional in the recipient cell. From this perspective, it is of great interest to see that one of the mRNAs we identify as transferred from adipocyte to macrophage, as well as confirmed in ExoCarta, is AHNAK. Given the most recent identification that AHNAK in the regulation of thermogenesis and lipolysis in WAT via β-adrenergic signalling makes this observation of AHNAK mRNA transfer in adipocyte exosomes to macrophages very relevant to the new field of extracellular vesicle biology and the possible impact for new therapeutics for T2D. It also highlights the value of this computational approach in generating novel hypothesis.
Currently our method is optimized for detecting nucleic acid exchanges between two cells lines at loci with known polymorphisms. It is known that exosome RNA contains not only mRNA, but also non-coding RNA species such as small microRNAs. For example, by analysing miRNA expression levels in a variety of cell lines and their derived exosomes, Guduric-Fuchs et al. found that miRNAs like miR-150, miR-142-3p, and miR-451 preferentially enter exosomes [57]. Huang et al. characterized the human plasma-derived exosomal RNAs by deep sequencing [58]. Although we did sequence miRNAs from the cells alongside total mRNA, we could however, not find sufficient genetic diversity amongst the miRNAs between the adipocyte and macrophages in order to allow for the identification of transfers via genetic differences mainly due to short length and relatively small number of miRNA species that could be detected.
Although accumulating evidence supports the important role of exosomes/EVs for mediating the cell-to-cell communication and exchange of genetic information, exosomes or lipid vesicles in general may not be the only mechanism for RNAs to transfer between cells. For example, it has been demonstrated that miRNA can be protected in the extracellular environment by forming complex with high-density lipoproteins (HDL) and RNA-binding proteins [59, 60]. Since our experimental design did not exclude other mechanisms of genetic material transfer besides EVs, we consider the identified transfers are likely mediated by exosomes but could also be contributed by other mechanisms. There are a few possibilities that can lead to false discoveries. First, reads from RNA-seq contain errors due to the technical limitation of the next generation sequencing technology [61]. For example, the raw quality scores of Illumina sequencing are calculated from the signal intensity, which do not always accurately represent the true error rate. Our study is particularly sensitive to this error rate since the transfer rate of transcripts is also quite low in this study as indicated by our simulation study. When the transfer rate is low (as seen in our experiment), the low signal/noise ratio could be the major factor for the high FDRs. Second, all the aligners may have some mapping errors. A comparison study with various sequence aligners shows that STAR-2pass with annotation has the best alignment performance but is still not completely error free [62]. Third, there are also genotyping error by SNP-arrays [63]. Finally, there are a few biological mechanisms such as RNA editing which will cause false discovery in our method.
Our simulation study has shown that our method has its ideal performance when the transfer rates between the two in vitro co-cultured cell lines are high. Inducing exosome secretion could be done via chemical means, such as altering cellular ceramide levels [64]. Our methods also require a good level of genetic diversity considering the donor and recipient systems. Optimization of genetic diversity can be done using cultures of cells that are knowingly from different genetic individuals such as in our experimental design. Other experimentally relevant systems to optimize genetic diversity, include investigations surrounding human derived exosomes (e.g., from plasma) and their function by injecting into the mouse [65]. In this case, one could survey the mouse tissues in order to identify tissues which have been impacted by the nucleic acid cargo of the donor exosome. Another naturally expressing genetic diversity is seen in cancers. It has been known that cancer cells secrete exosomes including during cell migration [66] and invasion [67]. Importantly, tumour derived microvesicles have been identified to contain mutated and amplified oncogenic DNA sequences and potentially have a role in genetic communication between cells as well as provide a potential source of tumour biomarkers. Thus our approach, we predict, would be highly amenable to identifying the genetic materials transferred between cancer and surrounding or distant normal cells [52] and further highlights the overall potential impact of this methodology.
Conclusions
In this study we present a novel systematic framework to call genes involved in the process of mRNA exchange between two co-cultured cell lines of different genetic backgrounds and investigate the role of exosomes as a vehicle in mediating the exchange. The systematic framework includes a protocol to perform quality control, alignment, mapping, and base call recalibration on the raw SNP array and RNA sequencing reads data, a Bayesian model to evaluate the significance of genotypic variation of a cell line under in vitro co-culture, and a method to estimate the rate of false discovered loci involving in the transfer process. By applying the framework to a co-culture between adipocyte and macrophage cell lines, we identified with high confidence, 8 mRNAs being transferred from macrophages to adipocytes and 21 mRNAs transferred in the opposite direction. These mRNAs represent biological functions including extracellular matrix, cell adhesion, glycoprotein, and signal peptides. We also estimate the transfer rate to be 0.001 in both directions. Our work provides a novel solution to studying EV mediated mRNA transfers between cells and this work is able to be extended to in vivo studies as well.
Abbreviations
- AD:
-
Adipocyte
- AUCs:
-
Area under curves
- CDCs:
-
Cardio sphere-derived cells
- DEGs:
-
Differentially expressed genes
- EVs:
-
Extracellular vesicles
- exRNA:
-
Extracellular RNA
- FDR:
-
False discovery rate
- FPKM:
-
Fragments per kilobase of transcript per million
- GO:
-
Gene ontology
- HDL:
-
High-density lipoproteins
- IGFBPs:
-
Insulin-like growth factor binding proteins
- IGV:
-
Integrative genomics viewer
- KO:
-
Knockout
- KS:
-
Kolmogorov–Smirnov
- MO:
-
Macrophage
- MVs:
-
Microvesicles
- NGS:
-
Next generation of sequencing
- PAPP-A:
-
Pregnancy-associated plasma protein-A
- PCA:
-
Principal component analysis
- QC:
-
Quality control
- RNA-seq:
-
RNA sequencing
- rRNAs:
-
Ribosomal RNAs
- SNP:
-
Single nucleotide polymorphism
- T2D:
-
Type 2 diabetes
References
Colombo M, Raposo G, Thery C. Biogenesis, secretion, and intercellular interactions of exosomes and other extracellular vesicles. Annu Rev Cell Dev Biol. 2014;30:255–89.
Deng Z-b, Poliakov A, Hardy RW, Clements R, Liu C, Liu Y, Wang J, Xiang X, Zhang S, Zhuang X, et al. Adipose tissue exosome-like vesicles mediate activation of macrophage-induced insulin resistance. Diabetes. 2009;58(11):2498–505.
Raposo G, Stoorvogel W. Extracellular vesicles: Exosomes, microvesicles, and friends. J Cell Biol. 2013;200(4):373–83.
Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007;9(6):654–9.
Skog J, Wurdinger T, van Rijn S, Meijer DH, Gainche L, Curry WT, Carter BS, Krichevsky AM, Breakefield XO. Glioblastoma microvesicles transport RNA and proteins that promote tumour growth and provide diagnostic biomarkers. Nat Cell Biol. 2008;10(12):1470–6.
Mathivanan S, Fahner CJ, Reid GE, Simpson RJ. ExoCarta 2012: database of exosomal proteins, RNA and lipids. Nucleic Acids Res. 2012;40(Database issue):D1241–4.
Ibrahim AGE, Cheng K, Marban E. Exosomes as critical agents of cardiac regeneration triggered by cell therapy. Stem Cell Rep. 2014;2(5):606–19.
Taylor DD, Gercel-Taylor C. MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol Oncol. 2008;110(1):13–21.
Escudier B, Dorval T, Chaput N, Andre F, Caby MP, Novault S, Flament C, Leboulaire C, Borg C, Amigorena S, et al. Vaccination of metastatic melanoma patients with autologous dendritic cell (DC) derived-exosomes: results of thefirst phase I clinical trial. J Transl Med. 2005;3(1):10.
Ohno S, Takanashi M, Sudo K, Ueda S, Ishikawa A, Matsuyama N, Fujita K, Mizutani T, Ohgi T, Ochiya T, et al. Systemically injected exosomes targeted to EGFR deliver antitumor microRNA to breast cancer cells. Molecular therapy : the journal of the American Society of Gene Therapy. 2013;21(1):185–91.
Vlassov AV, Magdaleno S, Setterquist R, Conrad R. Exosomes: current knowledge of their composition, biological functions, and diagnostic and therapeutic potentials. Biochim Biophys Acta. 2012;1820(7):940–8.
Lawson C, Vicencio JM, Yellon DM, Davidson SM. Microvesicles and exosomes: new players in metabolic and cardiovascular disease. J Endocrinol. 2016;228(2):R57–71.
Budnik V, Ruiz-Canada C, Wendler F. Extracellular vesicles round off communication in the nervous system. Nat Rev Neurosci. 2016;17(3):160–72.
Mittelbrunn M, Gutierrez-Vazquez C, Villarroya-Beltri C, Gonzalez S, Sanchez-Cabo F, Gonzalez MA, Bernad A, Sanchez-Madrid F. Unidirectional transfer of microRNA-loaded exosomes from T cells to antigen-presenting cells. Nat Commun. 2011;2:282.
Suganami T, Ogawa Y. Adipose tissue macrophages: their role in adipose tissue remodeling. J Leukoc Biol. 2010;88(1):33–9.
Garcia NA, Moncayo-Arlandi J, Sepulveda P, Diez-Juan A. Cardiomyocyte exosomes regulate glycolytic flux in endothelium by direct transfer of GLUT transporters and glycolytic enzymes. Cardiovasc Res. 2016;109(3):397–408.
Zheng PM, Chen L, Yuan XL, Luo Q, Liu Y, Xie GH, Ma YH, Shen LS. Exosomal transfer of tumor-associated macrophage-derived miR-21 confers cisplatin resistance in gastric cancer cells. J Exp Clin Cancer Res. 2017;36(1):53.
Kopylova E, Noe L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics (Oxford, England). 2012;28(24):3211–7.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics (Oxford, England). 2014;30(15):2114–20.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics (Oxford, England). 2013;29(1):15–21.
Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J et al: From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Current protocols in bioinformatics / editoral board, Andreas D Baxevanis [et al] 2013, 11(1110):11-10.
Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011;27(21):2987–93.
Anders S, McCarthy DJ, Chen Y, Okoniewski M, Smyth GK, Huber W, Robinson MD. Count-based differential expression analysis of RNA sequencing data using R and bioconductor. Nat Protoc. 2013;8(9):1765–86.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7(3):562–78.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Kolde R, Laur S, Adler P, Vilo J. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics (Oxford, England). 2012;28(4):573–80.
Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6.
Segura E, Nicco C, Lombard B, Veron P, Raposo G, Batteux F, Amigorena S, Thery C. ICAM-1 on exosomes from mature dendritic cells is critical for efficient naive T-cell priming. Blood. 2005;106(1):216–23.
Leinonen E, Hurt-Camejo E, Wiklund O, Hulten LM, Hiukka A, Taskinen M-R. Insulin resistance and adiposity correlate with acute-phase reaction and soluble cell adhesion molecules in type 2 diabetes. Atherosclerosis. 2003;166(2):387–94.
Kamiuchi K, Hasegawa G, Obayashi H, Kitamura A, Ishii M, Yano M, Kanatsuna T, Yoshikawa T, Nakamura N. Intercellular adhesion molecule-1 (ICAM-1) polymorphism is associated with diabetic retinopathy in type 2 diabetes mellitus. Diabet Med. 2002;19(5):371–6.
Khan T, Muise ES, Iyengar P, Wang ZV, Chandalia M, Abate N, Zhang BB, Bonaldo P, Chua S, Scherer PE. Metabolic Dysregulation and adipose tissue fibrosis: role of collagen VI. Mol Cell Biol. 2009;29(6):1575–91.
Dankel SN, Svard J, Mattha S, Claussnitzer M, Kloting N, Glunk V, Fandalyuk Z, Grytten E, Solsvik MH, Nielsen H-J, et al. COL6A3 expression in adipocytes associates with insulin resistance and depends on PPARgamma and adipocyte size. Obesity (Silver Spring). 2014;22(8):1807–13.
Conover CA, Harstad SL, Tchkonia T, Kirkland JL. Preferential impact of pregnancy-associated plasma protein-a deficiency on visceral fat in mice on high-fat diet. Am J Physiol-Endoc M. 2013;305(9):E1145–53.
Coste B, Mathur J, Schmidt M, Earley TJ, Ranade S, Petrus MJ, Dubin AE, Patapoutian A. Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels. Science. 2010;330(6000):55–60.
Maida Y, Yasukawa M, Furuuchi M, Lassmann T, Possemato R, Okamoto N, Kasim V, Hayashizaki Y, Hahn WC, Masutomi K. An RNA-dependent RNA polymerase formed by TERT and the RMRP RNA. Nature. 2009;461(7261):230–U104.
Davis GE, Senger DR. Endothelial extracellular matrix - biosynthesis, remodeling, and functions during vascular morphogenesis and neovessel stabilization. Circ Res. 2005;97(11):1093–107.
Ewens KG, George RA, Sharma K, Ziyadeh FN, Spielman RS. Assessment of 115 candidate genes for diabetic nephropathy by transmission/disequilibrium test. Diabetes. 2005;54(11):3305–18.
Li Z, Zhang H, Liu J, Liang C-P, Li Y, Li Y, Teitelman G, Beyer T, Bui HH, Peake DA, et al. Reducing plasma membrane sphingomyelin increases insulin sensitivity. Mol Cell Biol. 2011;31(20):4205–18.
Shin JH, Lee SH, Kim YN, Kim IY, Kim YJ, Kyeong DS, Lim HJ, Cho SY, Choi J, Wi YJ, et al. AHNAK deficiency promotes browning and lipolysis in mice via increased responsiveness to beta-adrenergic signalling. Sci Rep-Uk. 2016;6:23426.
Haase H. Ahnak, a new player in beta-adrenergic regulation of the cardiac L-type Ca2+ channel. Cardiovasc Res. 2007;73(1):19–25.
Borgonovo B, Cocucci E, Racchetti G, Podini P, Bachi A, Meldolesi J. Regulated exocytosis: a novel, widely expressed system. Nat Cell Biol. 2002;4(12):955–62.
Seyednasrollah F, Laiho A, Elo LL. Comparison of software packages for detecting differential expression in RNA-seq studies. Brief Bioinform. 2015;16(1):59–70.
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.
Matsuo Y, Tanaka M, Yamakage H, Sasaki Y, Muranaka K, Hata H, Ikai I, Shimatsu A, Inoue M, Chun T-H, et al. Thrombospondin 1 as a novel biological marker of obesity and metabolic syndrome. Metabolism. 2015;64(11):1490–9.
Inoue M, Jiang Y, Barnes RH 2nd, Tokunaga M, Martinez-Santibanez G, Geletka L, Lumeng CN, Buchner DA, Chun T-H. Thrombospondin 1 mediates high-fat diet-induced muscle fibrosis and insulin resistance in male mice. Endocrinology. 2013;154(12):4548–59.
Kent JW Jr, Comuzzie AG, Mahaney MC, Almasy L, Rainwater DL, VandeBerg JL, MacCluer JW, Blangero J. Intercellular adhesion molecule-1 concentration is genetically correlated with insulin resistance, obesity, and HDL concentration in Mexican Americans. Diabetes. 2004;53(10):2691–5.
Jelinek D, Millward V, Birdi A, Trouard TP, Heidenreich RA, Garver WS. Npc1 haploinsufficiency promotes weight gain and metabolic features associated with insulin resistance. Hum Mol Genet. 2011;20(2):312–21.
Lai EC. Notch signaling: control of cell communication and cell fate. Development. 2004;131(5):965–73.
Sorokin L. The impact of the extracellular matrix on inflammation. Nat Rev Immunol. 2010;10(10):712–23.
Yarwood SJ, Woodgett JR. Extracellular matrix composition determines the transcriptional response to epidermal growth factor receptor activation. Proc Natl Acad Sci U S A. 2001;98(8):4472–7.
Gibbings DJ, Ciaudo C, Erhardt M, Voinnet O. Multivesicular bodies associate with components of miRNA effector complexes and modulate miRNA activity. Nat Cell Biol. 2009;11(9):1143–9.
Balaj L, Lessard R, Dai L, Cho Y-J, Pomeroy SL, Breakefield XO, Skog J. Tumour microvesicles contain retrotransposon elements and amplified oncogene sequences. Nat Commun. 2011;2:180.
Kim D-K, Lee J, Simpson RJ, Lotvall J, Gho YS. EVpedia: a community web resource for prokaryotic and eukaryotic extracellular vesicles research. Semin Cell Dev Biol. 2015;40:4–7.
Kalra H, Simpson RJ, Ji H, Aikawa E, Altevogt P, Askenase P, Bond VC, Borras FE, Breakefield X, Budnik V, et al. Vesiclepedia: a compendium for extracellular vesicles with continuous community annotation. PLoS Biol. 2012;10(12):e1001450.
Keerthikumar S, Chisanga D, Ariyaratne D, Al Saffar H, Anand S, Zhao K, Samuel M, Pathan M, Jois M, Chilamkurti N, et al. ExoCarta: a web-based compendium of Exosomal cargo. J Mol Biol. 2016;428(4):688–92.
Hill AA, Reid Bolus W, Hasty AH. A decade of progress in adipose tissue macrophage biology. Immunol Rev. 2014;262(1):134–52.
Guduric-Fuchs J, O'Connor A, Camp B, O'Neill CL, Medina RJ, Simpson DA. Selective extracellular vesicle-mediated export of an overlapping set of microRNAs from multiple cell types. BMC Genomics. 2012;13:357.
Huang X, Yuan T, Tschannen M, Sun Z, Jacob H, Du M, Liang M, Dittmar RL, Liu Y, Liang M, et al. Characterization of human plasma-derived exosomal RNAs by deep sequencing. BMC Genomics. 2013;14:319.
Etheridge A, Gomes CP, Pereira RW, Galas D, Wang K. The complexity, function and applications of RNA in circulation. Front Genet. 2013;4:115.
Wang K, Zhang S, Weber J, Baxter D, Galas DJ. Export of microRNAs and microRNA-protective protein by mammalian cells. Nucleic Acids Res. 2010;38(20):7248–59.
Le H-S, Schulz MH, McCauley BM, Hinman VF, Bar-Joseph Z. Probabilistic error correction for RNA sequencing. Nucleic Acids Res. 2013;41(10):e109.
Engstrom PG, Steijger T, Sipos B, Grant GR, Kahles A, Ratsch G, Goldman N, Hubbard TJ, Harrow J, Guigo R, et al. Systematic evaluation of spliced alignment programs for RNA-seq data. Nat Methods. 2013;10(12):1185–91.
LaFramboise T. Single nucleotide polymorphism arrays: a decade of biological, computational and technological advances. Nucleic Acids Res. 2009;37(13):4181–93.
Trajkovic K, Hsu C, Chiantia S, Rajendran L, Wenzel D, Wieland F, Schwille P, Brugger B, Simons M. Ceramide triggers budding of exosome vesicles into multivesicular endosomes. Science. 2008;319(5867):1244–7.
Geiger A, Walker A, Nissen E. Human fibrocyte-derived exosomes accelerate wound healing in genetically diabetic mice. Biochem Biophys Res Commun. 2015;467(2):303–9.
Luga V, Zhang L, Viloria-Petit AM, Ogunjimi AA, Inanlou MR, Chiu E, Buchanan M, Hosein AN, Basik M, Wrana JL. Exosomes mediate stromal mobilization of autocrine Wnt-PCP signaling in breast cancer cell migration. Cell. 2012;151(7):1542–56.
Higginbotham JN, Demory Beckler M, Gephart JD, Franklin JL, Bogatcheva G, Kremers G-J, Piston DW, Ayers GD, McConnell RE, Tyska MJ, et al. Amphiregulin exosomes increase cancer cell invasion. Curr Biol. 2011;21(9):779–86.
Funding
JY is supported through Berg postdoc fellowship. ZT receives financial support from Berg Pharma as a consultant. This work is funded by a joint collaboration program between Berg Pharma and Icahn School of Medicine at Mount Sinai.
Availability of data and materials
All the data were uploaded to GEO with access https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=ozyfywuchjsnzwx&acc=GSE94155. The programs are available for download at http://research.mssm.edu/tulab/software/EVtransfer.html.
Author information
Authors and Affiliations
Contributions
ES conceived the overall experiment plan. CA, ZT, BZ, JZ, SH, AK, VV, VA, RS and NN designed the experiment. JY, KG, KA, SS JR, SG, BS and MM performed experiments. JY, ZT, and FP analyzed the data. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable
Consent for publication
Not applicable
Competing interests
The authors declare that they have no competing interests.
Additional files
Additional file 1:
A next generation sequencing based approach to identify extracellular vesicle mediated mRNA transfers between cells. (DOCX 1282 kb)
Additional file 2: Dataset S1.
Informaiton of filted loci from adipocyte to macrophage and from macrophage to adipocyte. (XLS 3896 kb)
Additional file 3: Dataset S2.
Details of mRNA transfer analysis from macrophage to adipocyte. (XLSX 1903 kb)
Additional file 4: Dataset S3.
Details of mRNA transfer analysis from adipocyte to macrophage. (XLSX 6308 kb)
Additional file 5: Dataset S4.
Differential genes between adipocyte cultured alone and co-cultured with macrophage and between macrophage alone and co-cultured with adipocyte. (XLSX 96 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Yang, J., Hagen, J., Guntur, K.V. et al. A next generation sequencing based approach to identify extracellular vesicle mediated mRNA transfers between cells. BMC Genomics 18, 987 (2017). https://doi.org/10.1186/s12864-017-4359-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-017-4359-1
Keywords
- Extracellular RNA (exRNA)
- Exosome
- Cell-to-cell communication
- RNA-seq
- Macrophages
- Adipocytes
- Bayesian method