- Research article
- Open access
- Published:
Comprehensive transcriptomic analyses of tissue, serum, and serum exosomes from hepatocellular carcinoma patients
BMC Cancer volume 19, Article number: 1007 (2019)
Abstract
Background
The expression of microRNAs (miRNAs) is a promising prognostic and diagnostic tool in hepatocellular carcinoma (HCC). Here we performed small RNA sequencing (sRNA-seq) of tissue, serum and serum exosomes to investigate changes in miRNA expression between the different sample types and correlated the expression with clinical parameters. We also performed gene expression arrays on tumor and normal tissue.
Results
Paired tissue, serum and serum exosomes sequencing revealed consistent positive correlation of miR-21 between serum exosomes and tumor tissue, indicating that miR-21 could be exported from tissue to circulation via exosomes. We found that let-7 miRNAs are generally upregulated in serum exosomes compared to whole serum, indicating that these miRNAs could be preferentially loaded into exosomes. Comparing serum from HCC patients with serum from healthy individuals revealed a global increase of miRNAs in serum from HCC patients, including an almost 4-fold increase of several miRNAs, including the liver-specific miR-122. When correlating miRNA expression with clinical parameters we detected significant association between hepatitis B virus (HBV) infection and miR-122 in serum as well as several serum and tissue-miRNAs that correlated with surgery type. We found that miR-141 and miR-146 correlated with cirrhosis in tumor tissue and normal tissue, respectively. Finally, high expression of miR-21 in tumors were associated with poor survival. Focusing on gene expression we found several significant messenger RNAs (mRNAs) between tumor and normal tissue and a Gene Ontology (GO) analysis revealed that these changes were mainly related to cell cycle and metabolism. Further, we detected mRNAs that correlated with cirrhosis and HBV infection in tissue. Finally, GO analysis of predicted targets for miRNAs down-regulated in tumor found that these were enriched for functions related to collagen synthesis.
Conclusions
Our combined data point to altered miRNA and mRNA expression contributing to both generally impaired lipid metabolism and increased cell proliferation and a miRNA-driven increase in collagen synthesis in HCC. Our results further indicate a correlation in miRNA expression between exosomes, serum, and tissue samples suggesting export from tumors via exosomes. This correlation could provide a basis for a more tumor-specific miRNA profile in serum.
Background
Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide with more than 780,000 new cases in 2012 (http://globocan.iarc.fr/Pages/fact_sheets_cancer.aspx). Common risk factors of HCC include alcohol consumption, hepatitis B virus (HBV), hepatitis C virus (HCV), and non-alcoholic fatty liver disease. Circulating biomarkers in serum and exosomes, such as microRNAs (miRNAs) and other small non-coding RNAs (ncRNAs) are promising tools for discovering and monitoring HCC [1,2,3,4,5]. For example, levels of miR-21 in serum exosomes are significantly higher in patients with HCC compared to controls [6, 7]. In whole serum, several studies have identified potential diagnostic miRNAs [8,9,10,11,12,13,14]; however, no consensus has been reached on a panel of miRNAs that robustly identifies HCC at an early stage of the disease. In tissue, several miRNAs are found to be dysregulated between tumor and normal tissue [15,16,17,18,19], for instance miR-21, miR-199, and miR-221. Different miRNAs are shown to correlate with the main risk factors of HCC, including liver cirrhosis [20], HCV [21, 22] and non-alcoholic fatty liver disease [23], indicating that miRNAs have a diagnostic potential in detecting HCC in an early phase of the disease. Other small non-coding RNAs, including small nucleolar RNAs (snoRNAs) and transfer RNAs (tRNAs) are abundant in HCC and have been implicated in the tumorigenesis of HCC [24,25,26].
MicroRNAs can be actively sorted into exosomes in complex with Argonaute (Ago) proteins. These exosomes can be released into circulation from apoptosis or from active export mechanisms [27] and may provide more tumor specific miRNA profiles than whole blood and serum. Specific exosomal miRNA profiles have been identified in HCC [2, 4, 28]. Furthermore, exosomes can transport miRNAs between cells, as shown for miR-122 between Huh7 and HepG2 human liver cancer cell lines [29] and for miRNAs between human and mouse liver cells and primary B cells [30]. This exosome-mediated transport can lead to the transported miRNAs targeting messenger RNAs (mRNAs) in the recipient cells. For example, exosomes containing miR-122 can lead to increased sensitivity to chemotherapy due to down-regulation of miR-122 target genes [31].
Here we investigated differences in miRNA expression in tissue, serum, and exosomes of HCC patients by sequencing small RNAs (sRNA-seq) in tumor tissue, normal tissue, serum exosomes, and whole serum from the same patient. We found that miRNAs upregulated in tumor tissue, including miR-21, were upregulated in serum exosomes, indicating that exosomes of tumor origin are enriched among serum exosomes. Further, we correlated miRNA and mRNA expression from HCC tissue and serum with clinical parameters and detected miRNAs and mRNAs associated with Cirrhosis and HBV.
Methods
Sample material
A total of 89 patients with primary HCC who underwent a curative liver resection and 56 patients who underwent liver transplantation at Fundeni Clinical Institute, Bucharest, Romania were included in this study. Liver tumor samples and adjacent non-tumor tissue were collected at the time of surgery in a stabilizing solution RNAlater (Sigma, St. Louis, MO), and stored at − 80 °C until collection of all samples. Plasma and serum aliquots were prepared from blood samples obtained prior to surgery and were stored at − 80 °C.
The study conformed to the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the Ethics Committee of the Fundeni Clinical Institute (30,884/22.10.2014). All patients signed a written informed consent. Follow-up was completed on 25 August 2016. The period of follow up was defined from the date of surgery to the date of patient’s death or the last follow-up point.
RNA and exosome isolation
For the small RNA sequencing, total RNA was isolated from tissue samples by using miRVana RNA isolation kit (ThermoFisher Scientific, Waltham, MA) following the total RNA procedure. RNA was eluted using RNase-free water and stored at − 80 °C. RNA from serum was isolated using the Qiagen Plasma/Serum kit (Qiagen, Hilden, Germany) using 200 μL of serum.
For the miRNA and mRNA microarray experiments, total RNA was isolated with TRIzol reagent according to the manufacturer’s instructions (Invitrogen, Carlsbad, CA).
Exosomes were isolated from serum using size exclusion chromatography following the protocol of Böing et al. [32]. RNA was subsequently isolated using the above mentioned Qiagen Plasma/Serum kit.
RNA quantification and quality assessment of isolated RNA
RNA purity and concentration was measured with NanoDrop™ ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA). For further assessment of RNA quality and relative size, samples were measured using Eukaryote total RNA assay on the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA) to detect the presence of small RNAs and to calculate RIN values. A RIN value > 9 was regarded as high quality and sufficient for sequencing.
Preparation of cDNA library for small RNA sequencing
The tissue validation dataset was prepared using NEBNext® Small RNA Library Prep (#E7330L) and the serum validation data set was prepared using the TruSeq small RNA protocol (Illumina, San Diego, CA), both according to the manufacturer’s instructions. The paired tissue, serum and exosome samples for the discovery dataset and the validation dataset was prepared using NEXTFLEX® Small RNA-Seq Kit (#NOVA-5132-06), according to the manufacturer’s instructions. PCR Amplification was performed using 13 cycles. Calibrator RNAs were added during the 3′ ligation step as previously described [33, 34]. The miRNA fragments were sequenced on the Illumina HiSeq 2500 and HiSeq 4000 systems (Illumina, San Diego, CA) using 50 base pair single read, at the Genomics Core Facility (GCF) in Trondheim, Norway.
Processing of sequence data
The raw data was processed according to Mjelle et al. [34] to identify mature miRNAs and isomiRs.
MicroRNA microarray analysis
The miRNA experiment was performed in Center of Excellence in Translational Medicine on Fundeni Clinical Hospital using One-colour Microarray-based Gene Expression (Agilent Technologies, Santa Clara, CA, USA) according to manufacturer protocol.
Briefly, total RNA (100 ng) was dephosphorylated, labeled with Cyanine 3-pCp using miRNA Complete Labeling and Hyb Kit. The labeled miRNAs were hybridized on G3 Human miRNA Microarray Kit, Release 21, 8 × 60 K. After hybridization and washing, slides were scanned with Agilent Microarray Scanner with SureScan High -Resolution Technology (G2505B).
Messenger RNA microarray analysis
Messenger RNA microarray was performed using the Agilent technology. Specifically, the samples were run using the protocol “One Color Microarray-Based Gene Expression Analysis, Low Input Quick Amp Labelling, version 6.9.1, August 2015”. Data were extracted from the machine using the software “Feature Extraction 12.0.3.1” and the microarray slide was scanned using “Agilent Microarray Scanner G250”.
The acquired data from miRNA and mRNA microarray was analyzed using Scan Control software and further extracted with Feature Extraction 12.0.3.1(Agilent Technologies) in order to obtain fluorescent image slides in TIF format and respectively raw data with signal intensity values.
Differential expression analysis
Differentially expressed miRNAs and isomiRs were identified using the Bioconductor package limma combined with voom transformation. To compare miRNA expression between samples, read counts were normalized using counts per million (cpm) normalization. For miRNAs, ncRNAs and isomiRs, we required expression of 1 cpm in at least 50% of the samples when comparing tumor vs normal. The miRNA microarray data was analyzed in R using the Bioconductor package AgiMicroRna. The mRNA microarray data was analyzed in R using the Bioconductor package limma and the functions normalizeBetweenArrays and backgroundCorrect.
Gene ontology analysis
Gene ontology analysis on mRNAs was performed using the Bioconductor package clusterProfiler and the enrichGO profiler with the parameters p value Cutoff = 0.01; pAdjustMethod = “BH”. Up-regulated and down-regulated transcripts were defined as transcripts with adjusted p-value below 0.05 and log2 fold change (logFC) above 0.5 or below − 0.5, respectively. The “universe” was defined as transcripts detected in at least 50% of the samples. Only probes that could be assigned a unique Entrez ID were included in the analysis. To reduce redundant GO terms we used the “simplify” function in clusterProfiler with default parameters.
Gene ontology analysis for miRNA targets
Sixty-nine miRNAs were differentially expressed in both miRNA sequencing data and in microarray data. TargetScan gave us ~ 8500 unique target genes for these miRNAs. We further filtered the number of miRNA-target gene pairs based on expression levels of miRNA (keeping highly expressed miRNAs having 75% percentile expression across samples > 10 log2 cpm), weighted context scores from TargetScan (less than − 0.4) and selecting the miRNA target gene pairs, if more than 75% of the number of miRNAs targeting one gene are either all up or all downregulated [35]. Further, we selected only the negatively correlated miRNA-target gene pairs. We finally had 31 miRNA and 315 unique target genes. Gene Ontology (GO) enrichment analysis was done for target genes of both up and downregulated miRNAs with gProfileR package in R.
Survival analysis
The survival analyses were performed in R by using the package “survival” and the survival-plots were generated using “ggsurvplot”. “High” and “Low”-expression in the plots indicate the 50% highest and lowest percentile of the expression values. The p-values were calculated form a coxph-model and were adjusted for multiple testing by using Bonferroni Correction.
Results
Overview of the study and data
Our study material consisted of liver tumor and adjacent normal liver samples from HCC patients, serum samples from the same patients and from non-cancer controls, and serum exosomes isolated from the patient serum samples (Additional file 1: Figure S1). We divided the participants into a discovery cohort, consisting of the 19 non-cancer controls and 17 patients (cases) with matched tissue, serum, and exosome samples, and a validation cohort, consisting of the remaining patients.
In the discovery cohort, we used sRNA-seq to characterize miRNAs and other small RNAs in the samples. The aim of this experiment was to identify candidate HCC biomarkers by investigating the distribution of small RNAs in the different biological sample types, identifying expression differences between tumor and normal tissue, patient and control serum, and patient serum and serum exosomes, and investigating correlations in expression between exosome, serum and tissue samples.
We then used the larger validation cohort to confirm the results from the discovery cohort comparisons and to identify RNAs that correlated with clinical parameters. Specifically, we did a second sRNA-seq experiment that included 80 serum samples, 8 matched exosome samples, 78 tumor tissue samples, 80 normal tissue samples (78 matched tumor and normal tissue samples). For additional validation, 73 of the matched tissue samples were analyzed on miRNA microarrays. Finally, we used mRNA microarrays to analyze 59 of the matched tissue samples to identify differentially expressed mRNAs, mRNAs correlating with clinical parameters, and target candidates for the differentially expressed miRNAs. The following sections detail our results.
MicroRNA expression in tissue, serum and serum exosomes
To detect differentially expressed miRNAs between sample classes and to identify potential correlations between the circulating miRNAs and tissue miRNAs, we sequenced small RNAs in tumor and normal tissue, and whole serum and serum exosomes collected from 17 patients (discovery cohort; Additional file 2: Table S1). For comparison, we sequenced small RNAs in whole serum from 19 non-cancer controls. On average, 11,460,874 reads mapped to the human genome. Of these, on average 3,509,859 reads aligned to miRBase and 5,496,217 reads aligned to the RNACentral database of ncRNAs (see Methods) (Additional file 3: Figure S2A). The composition of detected RNAs varied across the sample types; exosome samples were depleted for transfer RNAs (tRNAs) compared to serum and tissue, whereas tissue samples were enriched for small nucleolar RNAs (snoRNAs) compared to serum and exosome samples (Additional file 3: Figure S2B). When analyzing the relative expression of the different ncRNAs within each sample type we observed high expression of small cytoplasmic RNAs (scRNAs), which included Y RNAs, in the serum and exosome samples and high expression of snoRNAs in the tissue samples (Additional file 3: Figure S2C). All samples included a set of 10 calibrator RNAs that were added during sample preparation before sequencing and used for normalization in all analysis (Additional file 4: Table S2).
A principal component analysis (PCA) of mature miRNA expression found that the tissue samples were markedly different from the serum and exosome samples (Fig. 1a; PC1). Within the serum and exosome samples, control and cancer serum formed partially overlapping sub-clusters along the PCA’s second principal component (PC2), with the cancer exosome samples clustering between the two serum groups. In comparison, except for two cancer samples, the cancer and normal tissue samples formed nearly distinct clusters along PC2. Overall, the intra-group variation was lower for the tissue samples than for the serum and exosome samples, and normal tissue had the lowest intra-group variation in miRNA expression.
Differentially expressed miRNAs, isomiRs, and ncRNAs in tissue, exosomes, and serum
To investigate specific expression differences between the sample types for the 17 patients and control samples, we performed differential expression analysis on tumor vs normal tissue, exosomes vs cancer serum, and cancer serum vs control serum. We detected 58 significant miRNAs between exosomes and cancer serum (Additional file 5: Figure S3A), 40 significant miRNAs between cancer serum and control serum (Additional file 5: Figure S3B), and 56 significant miRNAs between tumor and normal tissue (Additional file 5: Figure S3C, Additional file 6: Table S3). In general, we observed different sets of significant miRNAs between the three comparisons, of which tumor tissue vs normal tissue had the highest number of unique significant miRNAs (Fig. 1b). One miRNA, miR-532-5p, was significant across all three comparisons.
Next, we focused on isomiRs (Additional file 7: Figure S4A) and other ncRNAs that did not match any known miRNAs. We detected 176 significant isomiRs between exosomes and cancer serum, 237 significant isomiRs between tumor and normal tissue and 123 significant isomiRs between cancer serum and control serum (Additional file 7: Figure S4B). For the ncRNAs, we detected 20 significant ncRNAs between tumor and normal tissue, 64 significant ncRNAs between exosomes and cancer serum and 48 significant ncRNAs between cancer serum and control serum. We observed an enrichment of down-regulated snoRNAs in tumor tissue vs normal tissue (Additional file 7: Figure S4C). In exosome vs cancer serum, most of the highly significant ncRNAs were down-regulate in exosomes (Additional file 7: Figure S4D). In cancer serum vs controls serum, we observed an enrichment of up-regulated snoRNAs (Additional file 7: Figure S4E).
Expression differences between tumor and normal tissue are partly mirrored between serum exosomes and whole serum
As tumor cells release miRNA-containing exosomes that enter circulation and can have both local and remote functions [36], we reasoned that tumor-associated miRNAs could be enriched in serum exosomes compared with whole serum. Of the 114 miRNAs that were significant in either the exosome vs cancer serum or tumor vs normal tissue comparison, 77 unique miRNAs were expressed in all four sample types and 37 miRNAs had identical signs of their log fold change (logFC) values between the two comparisons (Additional file 8: Figure S5A). Focusing on the six miRNAs that were significant in both comparisons (exosome vs cancer serum and tumor vs normal), five were either up-regulated (logFC> 0; 2 miRNAs) or down-regulated (logFC< 0; 3 miRNAs) in both tumors and exosomes compared with normal tissue and whole serum, respectively (Fig. 1c). The up-regulated miRNAs were the oncogenic miRNAs miR-21 and miR-10b, supporting that miRNA dysregulation in tumors is reflected in circulation.
To validate these results, we performed a second sequencing experiment that included eight patients with matched exosome and cancer serum samples and 75 patients with matched tumor and normal tissue samples of which five patients were matched across all four sample types (validation cohort; Additional file 2: Table S1). Because of the higher number of patients, we detected a higher number of miRNAs with significant expression differences between tumor and normal tissue when considering all 75 patients (Additional file 6: Table S3). All miRNAs that were significant between tumor and normal tissue in the discovery dataset showed consistent expression in the validation dataset, when using all 75 patients, except miR-183-3p and miR-34a-3p which were not detected in the validation dataset (Additional file 8: Figure S5B).
For the correlation analysis, we used data from the validation cohort patients that were matched across all four sample types. The results showed that five of the six miRNAs had similar direction of expression in this new cohort, with miR-21 being positively correlated in both cohorts (Fig. 1c).
As miR-21 reproducibly informed on tumor dysregulation within both tissue and serum samples, we next asked whether this miRNA could inform on patient survival. A survival analysis in the validation dataset found that high expression levels of miR-21 within tumor tissue was significantly associated with poorer survival (Fig. 1d). We found no significant association between patient survival and miRNA expression in serum or normal tissue.
Next, we investigated whether isomiRs and other ncRNAs also showed correlated expression differences between tumor and normal tissue and serum exosomes and whole serum. In the discovery cohort data, we detected 15 isomiRs that were significant in both exosome vs cancer serum and tumor vs normal; 13 of these isomiRs were either up-regulated or down-regulated in both tumors and exosomes compared with normal tissue and whole serum, respectively (Additional file 8: Figure S5C). In the validation cohort, 12 of these 15 isomiRs showed correlated expression differences, including isomiRs for miR-21, miR-122, miR-375 and other miRNAs.
Sixty-four ncRNAs were significant in exosome vs cancer serum, 20 ncRNAs were significant in tumor vs normal, and one ncRNA were significant in both comparisons. Most ncRNAs were up-regulated in tumor vs normal and downregulated in exosomes vs serum and only one of the significant ncRNAs were detected in the validation cohort.
In summary, the miRNA miR-21 showed increased mature and isomiR expression in tumor and serum exosomes compared with normal liver and cancer serum, supporting a model were tumor cells more than normal cells release miR-21-containing exosomes into circulation.
Correlating miRNA expression in serum samples with clinical parameters
Based on the second sRNA-seq experiment on 75 patients we asked if the expression of miRNAs in serum correlated with the patients’ clinical parameters. A PCA plot showed that sample age is one of the main factors explaining the variation in miRNA expression, suggesting that RNA degradation may contribute significantly to miRNA expression profiles (Additional file 9: Figure S6A). All further analysis on the serum samples were therefore adjusted for sample age.
We detected two miRNA, miR-194-5p and miR-122-3p, as significantly upregulated in patients with HBV infection (Table 1). When comparing resection versus transplantation regimens we detected 17 miRNAs that were significantly differentially expressed between the two treatment groups (Table 1).
MicroRNA expression in tissue correlates with surgery type, HBV infection, and cirrhosis
Having analyzed the serum samples, we went on investigating miRNA expression in tumor and normal tissue for the same patients. The tissue samples were analyzed using both sequencing and microarray approaches. First, we investigated the correlation between the two platforms by comparing fold-change values for the differences between tumor and normal. We detected 137 significant miRNAs using microarrays and 234 significant miRNAs using sequencing; 68 miRNAs were significant on both platforms (Additional file 6: Table S3). We observed high correlation between the platforms and when focusing on miRNAs that were detected on both platforms and significant on at least one of the platforms, 116 of 128 miRNAs (91%) were changing in the same direction (Additional file 9: Figure S6B). Similar as in the serum samples, we identified sample age as a major factor affecting miRNA expression profiles in the sequencing data (Additional file 9: Figure S6C-D). The microarray data, however, were not affected by sample age (Additional file 9: Figure S6E-F).
Further, we compared tissue miRNA expression with the clinical parameters, using the microarray data as a validation for the sequencing data. Comparing resection vs transplantation and adjusting for sample age, we identified 8 significant miRNAs in the tumor samples and 18 significant miRNAs in the normal samples (Table 2). For the tumor samples, six of the eight miRNAs were expressed in the same direction in the microarray; however, none were significant in the microarray. For the normal samples, five of the 18 significant miRNAs were expressed in the same direction in the microarray; however, none were significant in the microarray.
Comparing cirrhosis and non-cirrhosis, we identified miR-141-3p to have lower expression in the tumor samples of cirrhotic as compared with non-cirrhotic patients and miR-146b-5p to have higher expression in the normal samples of cirrhotic compared with non-cirrhotic patients (Table 2). The microarray data confirmed the expression differences for miR-146b-5p, whereas miR-141-3p was not detected as expressed in the microarray data.
Gene expression in tissue correlates with surgery type, HBV infection and cirrhosis
Next, we used microarrays to measure gene expression profiles in tumor and normal samples from 59 HCC patients. Of the 50,739 measured probes we detected 6280 differentially expressed transcripts between tumor and normal samples. Focusing on genes with an average log2 expression signal above 10, we observed a trend towards more down-regulated genes in tumor compared to normal, especially for the most significant genes (Additional file 10: Figure S7A). A PCA plot revealed clustering of the normal samples and more distributed tumor samples, indicating tumor samples to be more heterogeneous than normal samples (Additional file 10: Figure S7B). We observed no clustering related to the main clinical parameters (Additional file 10: Figure S7B). We performed gene ontology (GO) analysis on the significantly up-regulated and down-regulated transcripts (see Methods). The GO analysis showed that genes up-regulated in tumors were enriched for terms related to cell division and DNA replication, whereas genes down-regulated in tumors were enriched for terms related to metabolism and lipid (Additional file 11: Figure S8A). Subsequent KEGG analysis on the same two groups revealed enrichment of the terms cell cycle and DNA replication for the up-regulated genes and metabolism for the down-regulated genes (Additional file 11: Figure S8B). In addition, for the up-regulated genes we found enrichment of the terms alcoholism and viral carcinogenesis, two factors that are considered as major risk-factors in HCC.
Next, we compared gene expression with clinical parameters and detected significant genes with respect to cirrhosis, HBV, and treatment. Specifically, we detected 14 significant genes between cirrhosis and non-cirrhosis for the tumor samples; five significant genes between resection and transplantation for the normal samples and 13 significant genes between HBV and non-HBV for the normal samples (Table 3).
Finally, we validated the gene expression differences between tumor and normal samples by analyzing a publicly available HCC dataset. We analyzed the dataset of Grinchuk et al. that contained 115 primary tumor samples and 52 adjacent normal tissue samples [37]. A correlation analysis between the two datasets showed that the logFC values for the comparison tumor vs normal were highly comparable (r = 0.56) (Additional file 12: Figure S9). Two thousand three hundred eighty genes were significant in both datasets and the correlation value was 0.76 when only including these genes. These results indicate that the differences in gene expression observed in the current study is highly comparable with other studies from different cohorts.
Gene ontology analyses for miRNA target genes
Having established that both miRNAs and mRNAs are highly dysregulated in HCC tissue, we wanted to investigate if the differentially expressed miRNAs were targeting a specific set of mRNAs. We used 69 miRNAs that were differentially expressed in tumor vs normal tissue and identified their candidate target mRNAs using TargetScan [38] and additional filtering (see Methods). Gene ontology analysis showed that miRNAs down-regulated in tumor were targeting genes that were enriched for functions related to collagen, whereas miRNAs up-regulated in tumor were targeting genes that were enriched for functions related to regulation of neutrophils, legionellosis, and basal cell carcinoma (Additional file 13: Figure S10).
Discussion
We here performed paired small RNA profiling of tumor and normal tissue and circulating exosomes and serum from HCC patients. Comparing the differences between tumor and normal tissue and the differences between exosomes and serum we detected six miRNAs that differed in both comparisons. Four of the six miRNAs were concomitantly enriched (two miRNAs) or depleted (two miRNAs) within tumor tissue and serum exosomes in both the discovery and validation data. The most striking is miR-21, which is upregulated in tumor tissue and in serum exosomes as well as being associated with overall survival. These findings indicate that the expression levels in serum exosomes of some miRNAs, including miR-21, depend on exosomes released from tumor cells. Indeed, increased exosomal miR-21 levels has been associated with cancer in other studies. In ovarian cancer, colon cancer, and lung cancer, positive correlations were found between miR-21 in the tumor and circulating miR-21 in exosomes [39,40,41]. It has been suggested that exosomal miR-21 could be used as a universal biomarker for cancer in combination with other clinical parameters [42]. In HCC, one study has shown a positive correlation between exosomal miR-21 in cell lines and miR-21 in the cell culture supernatant [43]. Another study showed significantly higher miR-21 levels in exosomes than in exosome-depleted serum or the whole serum [6]. Together, several previous results show that exosomal miR-21 could be a good biomarker for HCC, supporting our findings on patients’ samples.
Comparing miRNA expression in exosomes and whole serum we observed enrichment of several let-7 miRNAs in the exosomes. One study on gastric cancer cell lines showed that let-7 was secreted from cells via exosomes into the extracellular environment and thereby maintain the oncogenic properties of the cells [44]. Potentially more relevant to our results is the finding by Okoye et al. showing that Foxp3+ T regulatory (Treg) cells released exosomes containing let-7, and that these exosomes were further transferred to T helper 1 (Th1) cells, suppressing Th1 cell proliferation. This means that the increased let-7d levels we observed could be due to immune responses that cause exosomal let-7d, and potentially other let-7 members, to be released into circulation [45].
Comparing miRNA expression in serum from cancer patients and controls we observed a strong and highly significant upregulation of miR-122. MicroRNA miR-122 is highly expressed and strongly enriched in hepatocytes and is frequently downregulated in HCC tissue [46], which is shown to result in metastatic properties of hepatocytes [47]. Indeed, we found reduced levels of miR-122 in tumor tissue in our discovery set (p-value = 0.059), and a significant reduction in the validation set and in the microarray data (Additional file 6: Table S3). Reports also show that miR-122 is increased in serum of HCC patients [48, 49], in line with our results. In addition to increased level of miR-122, we generally observed a global upregulation of miRNAs in serum of cancer patients. Also, the miRNAs that were upregulated were affected to a much larger extent than the miRNAs that were downregulated.
Focusing on other small ncRNAs, we observed several significant ncRNAs in the different comparisons. In tissue, we observed a general downregulation of snoRNAs in tumor compared to normal. Studies have shown that snoRNAs are indeed downregulated in tumors [50] and in HCC [51]. In cancer serum vs control serum, snoRNAs were generally up-regulated and tRNAs were generally downregulated. In exosomes vs serum, we observed a similar trend as for miRNAs where most of the significant ncRNAs were downregulated in exosomes compared to serum.
A general theme for the sequencing data was that sample age was a confounding factor in the data, both in serum and in tissues. We observed separate clustering of samples collected recently compared to samples collected several years ago. Interestingly, the miRNA microarray data did not show similar clustering, indicating that the signal from the microarray probes could be less sensitive to RNA degradation or factors related to the samples’ age. In the tissue samples, the sample age bias is less evident than in the tumor samples, but is clearly visible in the normal samples. Notably, sample age largely overlaps with treatment type (resection vs transplantation), meaning that we cannot accurately predict if the bias is due to treatment or sample age. However, no clear separation with regard to treatment is found in the microarray data or in the tumor tissue.
Analyzing the tissue mRNA data, we found a large set of differentially expressed transcripts. The most significant of these mRNAs were regulated in tumor compared to normal tissue. Gene ontology analysis of the differentially expressed mRNAs revealed that cell cycle terms were enriched among genes upregulated in tumors. This is expected, as cell cycle genes often are dysregulated in HCC and cancer in general [52]. Downregulated genes were enriched for terms related to metabolism in general and lipid processes in particular. The liver plays an important role in catabolism of plasma lipoproteins and these processes are often impaired in liver cancer [53]. We also observed downregulation of some inflammation and complement related genes including CXCL14 and CFP. The chemokine CXCL14 is a known tumor suppressor involved in antimicrobial immunity and inflammatory processes [54]. Enriched KEGG pathways show similar biological functions as the gene ontologies. Here, alcoholism is enriched among the upregulated genes, in addition to cell cycle, microRNAs, and viral carcinogenesis. Alcohol and viral infections, especially hepatitis B and hepatitis C virus, are among the most significant risk factors of HCC. When comparing our gene expression results with the dataset of Grinchuk et al., the significant genes were generally changing in the same direction in the two datasets. Interestingly, and in accordance with the results above, the most significantly down-regulated genes common in both datasets were related to immunity. For instance, CLEC1B and CLEC4G are involved in T-cell immune responses and FCN2 is a liver-specific gene related to complement activation.
Finally, we investigated the gene ontologies for the predicted target genes of the dysregulated mRNAs in tissue. Here, several terms related to collagen were enriched among the mRNAs being targeted by miRNAs that were downregulated in tumor. This means that these genes have potentially less miRNA regulation and thereby increased expression in HCC tumors. Indeed, studies have shown that collagen related genes are upregulated in HCC and contribute to cancer progression [55, 56].
Conclusions
This study indicates a correlation in miRNA and isomiR expression between exosomes and tissue, suggesting that these RNAs, including miR-21, are exported from tumors into circulation via exosomes. Further, the study presents a comprehensive profiling of miRNAs, other small ncRNAs, and mRNAs in HCC. Several of these RNAs are associated with clinical parameters such as Cirrhosis and HBV infection. We further found that high tumor expression of miR-21 is associated with poorer survival, pointing towards miR-21 as a potential prognostic and diagnostic biomarker for HCC. Gene ontology analyses of altered miRNA and mRNA expression pointed towards impaired lipid metabolism, increased cell proliferation and a miRNA-driven increase in collagen synthesis in HCC.
Availability of data and materials
All the data obtained and materials analyzed in this research are available with the corresponding author upon reasonable request.
Abbreviations
- Ago:
-
Argonaute
- CPM:
-
Counts per million
- DNA:
-
Deoxy-ribonucleic acid
- GO:
-
Gene Ontology
- HBV:
-
Hepatitis B virus
- HCC:
-
Hepatocellular carcinoma
- HCV:
-
Hepatitis C virus
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- logFC:
-
Log fold change
- miRNA:
-
MicroRNA
- miRs:
-
MicroRNA
- mRNA:
-
Messenger RNA
- ncRNAs:
-
Non-coding RNA
- OR:
-
odds ratio
- PC:
-
Principal component
- PCA:
-
Principal component analysis
- r:
-
Pearson correlation coefficient
- RFA:
-
Radiofrequency ablation
- RNA:
-
Ribonucleic acid
- scRNAs:
-
Small cytoplasmic RNAs
- snoRNAs:
-
Small nucleolar RNA
- sRNA-seq:
-
Small RNA sequencing
- Th1:
-
T helper 1
- tRNA:
-
Transfer RNA
References
Yang J, Li C, Zhang L, Wang X. Extracellular vesicles as carriers of non-coding RNAs in liver diseases. Front Pharmacol. 2018;9:415.
Li LM, Liu ZX, Cheng QY. Exosome plays an important role in the development of hepatocellular carcinoma. Pathol Res Pract. 2019;215(8):152468.
Peng C, Ye Y, Wang Z, Guan L, Bao S, Li B, et al. Circulating microRNAs for the diagnosis of hepatocellular carcinoma. Dig Liver Dis. 2019;51(5):621–31.
Jun L, Yang G, Zhisu L. The utility of serum exosomal microRNAs in hepatocellular carcinoma. Biomed Pharmacother. 2019;111:1221–7.
Sadri Nahand J, Bokharaei-Salim F, Salmaninejad A, Nesaei A, Mohajeri F, Moshtzan A, et al. microRNAs: key players in virus-associated hepatocellular carcinoma. J Cell Physiol. 2019;234(8):12188–225.
Wang H, Hou L, Li A, Duan Y, Gao H, Song X. Expression of serum exosomal microRNA-21 in human hepatocellular carcinoma. Biomed Res Int. 2014;2014:864894.
Xu J, Wu C, Che X, Wang L, Yu D, Zhang T, et al. Circulating microRNAs, miR-21, miR-122, and miR-223, in patients with hepatocellular carcinoma or chronic hepatitis. Mol Carcinog. 2011;50(2):136–42.
Zhang Y, Li T, Qiu Y, Zhang T, Guo P, Ma X, et al. Serum microRNA panel for early diagnosis of the onset of hepatocellular carcinoma. Medicine (Baltimore). 2017;96(2):e5642.
Liao Q, Han P, Huang Y, Wu Z, Chen Q, Li S, et al. Potential role of circulating microRNA-21 for hepatocellular carcinoma diagnosis: a meta-analysis. PLoS One. 2015;10(6):e0130677.
Shen J, Wang A, Wang Q, Gurvich I, Siegel AB, Remotti H, et al. Exploration of genome-wide circulating microRNA in hepatocellular carcinoma: MiR-483-5p as a potential biomarker. Cancer Epidemiol Biomark Prev. 2013;22(12):2364–73.
Dhayat SA, Husing A, Senninger N, Schmidt HH, Haier J, Wolters H, et al. Circulating microRNA-200 family as diagnostic marker in hepatocellular carcinoma. PLoS One. 2015;10(10):e0140066.
El-Abd NE, Fawzy NA, El-Sheikh SM, Soliman ME. Circulating miRNA-122, miRNA-199a, and miRNA-16 as biomarkers for early detection of hepatocellular carcinoma in Egyptian patients with chronic hepatitis C virus infection. Mol Diagn Ther. 2015;19(4):213–20.
Huang YH, Liang KH, Chien RN, Hu TH, Lin KH, Hsu CW, et al. A circulating MicroRNA signature capable of assessing the risk of hepatocellular carcinoma in cirrhotic patients. Sci Rep. 2017;7(1):523.
Zhang EL, Gu J, Zhang ZY, Dong KS, Liang BY, Huang ZY. MicroRNA expression profiling in patients with hepatocellular carcinoma of familial aggregation and hepatitis B virus infection. Oncol Lett. 2017;14(1):971–6.
Yang J, Han S, Huang W, Chen T, Liu Y, Pan S, et al. A meta-analysis of microRNA expression in liver cancer. PLoS One. 2014;9(12):e114533.
Gao B, Ning S, Li J, Liu H, Wei W, Wu F, et al. Integrated analysis of differentially expressed mRNAs and miRNAs between hepatocellular carcinoma and their matched adjacent normal liver tissues. Oncol Rep. 2015;34(1):325–33.
Hou J, Lin L, Zhou W, Wang Z, Ding G, Dong Q, et al. Identification of miRNomes in human liver and hepatocellular carcinoma reveals miR-199a/b-3p as therapeutic target for hepatocellular carcinoma. Cancer Cell. 2011;19(2):232–43.
Meng F, Henson R, Wehbe-Janek H, Ghoshal K, Jacob ST, Patel T. MicroRNA-21 regulates expression of the PTEN tumor suppressor gene in human hepatocellular cancer. Gastroenterology. 2007;133(2):647–58.
Fornari F, Gramantieri L, Ferracin M, Veronese A, Sabbioni S, Calin GA, et al. MiR-221 controls CDKN1C/p57 and CDKN1B/p27 expression in human hepatocellular carcinoma. Oncogene. 2008;27(43):5651–61.
Guo L, Li W, Hu L, Zhou H, Zheng L, Yu L, et al. Diagnostic value of circulating microRNAs for liver cirrhosis: a meta-analysis. Oncotarget. 2018;9(4):5397–405.
Mourad L, El-Ahwany E, Zoheiry M, Abu-Taleb H, Hassan M, Ouf A, et al. Expression analysis of liver-specific circulating microRNAs in HCV-induced hepatocellular carcinoma in Egyptian patients. Cancer Biol Ther. 2018;19(5):400–6.
Amr KS, Atia HA, Elbnhawy RA, Ezzat WM. Early diagnostic evaluation of miR-122 and miR-224 as biomarkers for hepatocellular carcinoma. Genes Dis. 2017;4(4):215–21.
Liu CH, Ampuero J, Gil-Gomez A, Montero-Vallejo R, Rojas A, Munoz-Hernandez R, et al. miRNAs in patients with non-alcoholic fatty liver disease: a systematic review and meta-analysis. J Hepatol. 2018;69(6):1335–48.
Wang H, Ma P, Liu P, Chen B, Liu Z. Small nucleolar RNA U2_19 promotes hepatocellular carcinoma progression by regulating Wnt/beta-catenin signaling. Biochem Biophys Res Commun. 2018;500(2):351–6.
Wu L, Zheng J, Chen P, Liu Q, Yuan Y. Small nucleolar RNA ACA11 promotes proliferation, migration and invasion in hepatocellular carcinoma by targeting the PI3K/AKT signaling pathway. Biomed Pharmacother. 2017;90:705–12.
Selitsky SR, Baran-Gale J, Honda M, Yamane D, Masaki T, Fannin EE, et al. Small tRNA-derived RNAs are increased and more abundant than microRNAs in chronic hepatitis B and C. Sci Rep. 2015;5:7675.
Bell E, Taylor MA. Functional roles for exosomal microRNAs in the tumour microenvironment. Comput Struct Biotechnol J. 2017;15:8–13.
Kogure T, Lin WL, Yan IK, Braconi C, Patel T. Intercellular nanovesicle-mediated microRNA transfer: a mechanism of environmental modulation of hepatocellular cancer cell growth. Hepatology. 2011;54(4):1237–48.
Basu S, Bhattacharyya SN. Insulin-like growth factor-1 prevents miR-122 production in neighbouring cells to curtail its intercellular transfer to ensure proliferation of human hepatoma cells. Nucleic Acids Res. 2014;42(11):7170–85.
Pan Q, Ramakrishnaiah V, Henry S, Fouraschen S, de Ruiter PE, Kwekkeboom J, et al. Hepatic cell-to-cell transmission of small silencing RNA can extend the therapeutic reach of RNA interference (RNAi). Gut. 2012;61(9):1330–9.
Lou G, Song X, Yang F, Wu S, Wang J, Chen Z, et al. Exosomes derived from miR-122-modified adipose tissue-derived MSCs increase chemosensitivity of hepatocellular carcinoma. J Hematol Oncol. 2015;8:122.
WTO. 2018. http://globocan.iarc.fr/Pages/fact_sheets_cancer.aspx. Accessed 1 Oct 2019.
Hafner M, Renwick N, Farazi TA, Mihailovic A, Pena JT, Tuschl T. Barcoded cDNA library preparation for small RNA profiling by next-generation sequencing. Methods. 2012;58(2):164–70.
Mjelle R, Sellaeg K, Saetrom P, Thommesen L, Sjursen W, Hofsli E. Identification of metastasis-associated microRNAs in serum from rectal cancer patients. Oncotarget. 2017;8(52):90077–89.
Olsen LC, O'Reilly KC, Liabakk NB, Witter MP, Saetrom P. MicroRNAs contribute to postnatal development of laminar differences and neuronal subtypes in the rat medial entorhinal cortex. Brain Struct Funct. 2017;222(7):3107–26.
Maia J, Caja S, Strano Moraes MC, Couto N, Costa-Silva B. Exosome-based cell-cell communication in the tumor microenvironment. Front Cell Dev Biol. 2018;6:18.
Grinchuk OV, Yenamandra SP, Iyer R, Singh M, Lee HK, Lim KH, et al. Tumor-adjacent tissue co-expression profile analysis reveals pro-oncogenic ribosomal gene signature for prognosis of resectable hepatocellular carcinoma. Mol Oncol. 2018;12(1):89–113.
Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009;19(1):92–105.
Rabinowits G, Gercel-Taylor C, Day JM, Taylor DD, Kloecker GH. Exosomal microRNA: a diagnostic marker for lung cancer. Clin Lung Cancer. 2009;10(1):42–6.
Taylor DD, Gercel-Taylor C. MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol Oncol. 2008;110(1):13–21.
Tsukamoto M, Iinuma H, Yagi T, Matsuda K, Hashiguchi Y. Circulating exosomal MicroRNA-21 as a biomarker in each tumor stage of colorectal cancer. Oncology. 2017;92(6):360–70.
Shi J. Considering exosomal miR-21 as a biomarker for cancer. J Clin Med. 2016;5(4):42.
Fornari F, Ferracin M, Trere D, Milazzo M, Marinelli S, Galassi M, et al. Circulating microRNAs, miR-939, miR-595, miR-519d and miR-494, identify cirrhotic patients with HCC. PLoS One. 2015;10(10):e0141448.
Ohshima K, Inoue K, Fujiwara A, Hatakeyama K, Kanto K, Watanabe Y, et al. Let-7 microRNA family is selectively secreted into the extracellular environment via exosomes in a metastatic gastric cancer cell line. PLoS One. 2010;5(10):e13247.
Okoye IS, Coomes SM, Pelly VS, Czieso S, Papayannopoulos V, Tolmachova T, et al. MicroRNA-containing T-regulatory-cell-derived Exosomes suppress pathogenic T helper 1 cells. Immunity. 2014;41(3):503.
Kutay H, Bai S, Datta J, Motiwala T, Pogribny I, Frankel W, et al. Downregulation of miR-122 in the rodent and human hepatocellular carcinomas. J Cell Biochem. 2006;99(3):671–8.
Coulouarn C, Factor VM, Andersen JB, Durkin ME, Thorgeirsson SS. Loss of miR-122 expression in liver cancer correlates with suppression of the hepatic phenotype and gain of metastatic properties. Oncogene. 2009;28(40):3526–36.
An Y, Gao S, Zhao WC, Qiu BA, Xia NX, Zhang PJ, et al. Novel serum microRNAs panel on the diagnostic and prognostic implications of hepatocellular carcinoma. World J Gastroenterol. 2018;24(24):2596–604.
Qi P, Cheng SQ, Wang H, Li N, Chen YF, Gao CF. Serum microRNAs as biomarkers for hepatocellular carcinoma in Chinese patients with chronic hepatitis B virus infection. PLoS One. 2011;6(12):e28486.
Mannoor K, Liao J, Jiang F. Small nucleolar RNAs in cancer. Biochim Biophys Acta. 2012;1826(1):121–8.
Baral D, Wu L, Katwal G, Yan X, Wang Y, Ye Q. Clinical significance and biological roles of small nucleolar RNAs in hepatocellular carcinoma. Biomed Rep. 2018;8(4):319–24.
Greenbaum LE. Cell cycle regulation and hepatocarcinogenesis. Cancer Biol Ther. 2004;3(12):1200–7.
Jiang J, Nilsson-Ehle P, Xu N. Influence of liver cancer on lipid and lipoprotein metabolism. Lipids Health Dis. 2006;5:4.
Lu J, Chatterjee M, Schmid H, Beck S, Gawaz M. CXCL14 as an emerging immune and inflammatory modulator. J Inflamm (Lond). 2016;13:1.
Zheng X, Liu W, Xiang J, Liu P, Ke M, Wang B, et al. Collagen I promotes hepatocellular carcinoma cell proliferation by regulating integrin beta1/FAK signaling pathway in nonalcoholic fatty liver. Oncotarget. 2017;8(56):95586–95.
Fang M, Yuan J, Peng C, Li Y. Collagen as a double-edged sword in tumor progression. Tumour Biol. 2014;35(4):2871–82.
Acknowledgements
We thank the Genomics Core Facility at The Norwegian University of Science and Technology for performing sequencing.
Funding
The research leading to these results has received funding from EEA Financial Mechanism 2009–2014 under the project contract no 4SEE/30.06.2014. The Programme Operator of EEA contributed to project dissemination by organizing the “Communication and Publicity Seminar” and “Research Ethics Seminar” that provided tailored information and training during project implementation. In this way, the collection of data was made in line with international and national regulations. Moreover, the program operator constantly ensured the project monitoring and contributed to the financial, scientific and technical report, which improved data collection and study design.
Author information
Authors and Affiliations
Contributions
RM: Preparing the manuscript, data analysis, and producing sequencing data. SOD: Conceived the study design, involved in biological sample collection and patients follow-up. NB: Involved in procurement of biological sample. KC: Data analysis. AS: Performed microarray experiments. DC: Performed exosomes isolation. VH: Involved in histopathological review of the samples. PS: Data analysis and statistics. IP: conceived and directed the study. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Ethical approval number: 30884/22.10.2014. The study is performed in accordance with the Declaration of Helsinki.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1:
Figure S1. Study overview.
Additional file 2: Table S1.
Patient sample overview.
Additional file 3: Figure S2.
Sequencing statistics.
Additional file 4: Table S2.
Calibrator RNAs.
Additional file 5: Figure S3.
Volcano plots miRNAs.
Additional file 6: Table S3.
Differentially expressed miRNAs in tissue.
Additional file 7: Figure S4.
IsomiRs and other ncRNAs.
Additional file 8: Figure S5.
Correlation between exosomes and tissue.
Additional file 9: Figure S6.
PCA plot of miRNAs.
Additional file 10: Figure S7.
Volcano and PCA plot of mRNA data.
Additional file 11: Figure S8.
GO plots of mRNA data.
Additional file 12: Figure S9.
Correlation of mRNAs in Mjelle et al. and Grinchuk et al.
Additional file 13: Figure S10.
GO plot of miRNA target data.
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
Mjelle, R., Dima, S.O., Bacalbasa, N. et al. Comprehensive transcriptomic analyses of tissue, serum, and serum exosomes from hepatocellular carcinoma patients. BMC Cancer 19, 1007 (2019). https://doi.org/10.1186/s12885-019-6249-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12885-019-6249-1