Skip to main content

Alterations in MicroRNA gene expression profile in liver transplant patients with hepatocellular carcinoma



Hepatocellular carcinoma (HCC) can lead to liver failure which renders to liver transplant. miRNAs might be detected as biomarkers in subclinical stage of several hepatobiliary disorders like HCC. Therefore, in the present study, alterations in miRNAs as biomarkers were detected in LT patients with HCC.


Fourteen tissue samples composed of 5 rejected and 9 non-rejected ones were used for studying the miRNAs expression pattern using LNA-array probe assay and the result was evaluated by in house SYBR Green Real-time PCR protocols on 30 other tissue samples composed of 10 rejected and 20 non-rejected ones for the selected miRNAs. All samples were collected from liver transplanted patients with HCC.


The study results revealed that in rejected patients compared to non-rejected ones, hsa-miR-3158-5p, -4449, -4511, and -4633-5p were up-regulated and hsa-miR-122-3p, -194-5p, 548as-3p, and -4284 were down-regulated. ROC curve analysis also confirmed that miR194-5p and -548as-3p in up-regulated and also, miR-3158-5p, -4449 in down-regulated microRNAs are significantly important molecules in rejection.


Finally, the tissue levels of specific miRNAs (especially hsa-miR-3158-5p, -4449, -194-5p and -548as-3p) significantly correlated with the development of HCC, which can be present as biomarkers after further completing studies.

Peer Review reports


The main treatment approach in many end stage liver failures is liver transplant (LT) [1]. Ranked as the fifth most common diagnosed cancer worldwide, hepatocellular carcinoma (HCC) refers to a multistage process, involving many genes [2], which can render to hepatocyte turnover, oxidative DNA damage, and inflammation [3]. Even though innovated therapeutic strategies are being established, handling advanced HCC is poorly efficacious and merely at the point of diagnosis. Apparently, transplantation seems to be the only known effective treatment approach for this disease [4].

A considerable part of human genome is consisted of non-coding regions, including non-coding RNAs, known as microRNAs (miRNAs) [1]. These regions encode dynamic miRNAs, regulating around 20–30% of mammalian genes that comprise nearly 5% of transcriptome [1, 5]. MiRNAs are evolutionary conserved, made up of 18–25 nucleotides. They are also involved in various biological processes e.g. cell growth, apoptosis, hematopoietic lineage differentiation, and gene regulation [1, 6, 7].

There are also many miRNAs with specific expression patterns in the liver at different stages [8]. Investigations have demonstrated that some miRNAs are involved in regulating many metabolic pathways. Changes in miRNA expression levels might reflect upon the underlying sources of inflammation [9]. The research which was conducted by Murakami et al., in 2006 was one of the pioneer studies on investigation the role of miRNAs in HCC development [10]. After that a number of studies have been investigated the miRNA changes in HCC [11,12,13,14,15]. As an illustration, it has been reported that miR-17-92, miR-21, miR-221, miR-222, and miR-224 are frequently up-regulated in HCC. Furthermore, the importance of miRNAs in liver function has been broadly documented in the related literature [16]. Studies have also revealed that miR-122 as a blood biomarker changes in viral, alcoholic, and chemical injuries induced in the liver [8]. The importance of miRNAs is also, studied in LT rejection and it seems that some of miRNAs are playing specific role in this regard. Ectopic expression of serum miRNA signatures is known to have prognostic, diagnostic, and biologic impact in liver allograft rejection [17]. Other sources showed that in rat model of LT rejection, miR-146a, 15b, 223, 23a, 27a, 34a and 451 were upregulated [18].

Furthermore, different studies have been used the LNA-array technique for detection of miRNAs’ level in different samples and discussed the advantages of this method over other related methods. The ability of this method to produce high affinity hybridizations render to accurate discriminating signals, have been announced in earlier reports [5, 19,20,21].

Based on invasiveness and time-consuming nature of the pathological tests, mainly biopsy, the necessity for a faster and less invasive method seems crucial. Therefore, LNA-array (locked nucleic acid-array) probe and real-time polymerase chain reaction (RT-PCR) was utilized to detect alterations in miRNA expression in the tissues of LT patients with HCC. This study is tried to elucidate the role of miRNAs in regulating the rejection of liver transplant in HCC patients.


Patients and samples

This study was fulfilled in the Pathology Ward of Namazi Hospital, Shiraz, Iran (2016–2018). To conduct this study, a total number of 44 formalin-fixed paraffin-embedded (FFPE) samples were collected from LT patients with HCC. Macro-dissection was also performed. In brief, there were attempts to adopt a standard procedure to dissect the middle region of the tumor to harvest a proper sample, surrounded by the cancerous cells. The harvesting procedure was then confirmed microscopically to ensure the nonexistence of non-cancerous cells in each sample. All of the participants were asked to completed a questionnaire on demographic characteristics, history of other cancers except HCC, alcohol use, and tobacco use. Those with a history of other cancers, alcohol use, or tobacco use were excluded from the study.

As this study was a retrospective one, the samples were collected and divided into acute rejected and nonrejected patients due to their pathology reports. The pathology results were closely examined by expert pathologists and all of tissue samples that were classified as rejected were at stages III and IV based on AJCC (American Joint Committee on Cancer) staging system [22].

The samples were divided into two main groups. The 1st group (test group) was composed of 14 liver transplant samples (5 rejected and 9 nonrejected). The 2nd group (validation group) was composed of 30 samples (10 rejected and 20 nonrejected). For more clarification the sample division is shown in Fig. 1. All the protocols applied in this study were in complete agreement with the Helsinki Declaration and its later amendment. The local Ethics Committee of Shiraz University of Medical Sciences also approved each stage of this study. After explaining the study objectives, written informed consent was further obtained from each patient.

Fig. 1

The study sample division chart. The studied patents samples were divided into two categories. The first category was composed of 14 individuals (test group) which were participated in LNA-array probe test and were composed of 5 rejected and 9 nonrejected samples. The second group was consisted of 30 individuals (validation group). These patient samples were used for validating the results of LNA-array test

De-paraffinization of FFPE samples, RNA isolation and quality control

Using a microtome, 10 cuts of 5 μm of each sample were prepared. All the samples were RNA isolated, using Trizol™ (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. Sample RNA Quality Control was further performed by an Agilent 2100 Bioanalyzer, to provide an electropherogram for each sample. In addition to the measurement of the rRNA ratio in a traditional manner (28S/18S), the Bioanalyzer provided an RNA Integrity Number (RIN) value (ranging from 0 to 10) and a reliable impression of RNA quality. The recommended RIN values were higher than 7 for good array performance. The NanoDrop instrument was also used to measure the concentration (A260), protein contamination (ratio A260/A280), and contamination with buffer components, or organic compounds (ratio A260/A230) accurately.

Sample preparation for LNA-array probe

Labeling, hybridization, and normalization

Using the miRCURY LNA™ miRNA Hi-Power Labeling Kit (Exiqon, Denmark), 400 ng total RNA were labeled by Hy3™ and Hy5™ fluorescent label for labeling the samples and the references, respectively to produce highly efficient and uniform labeling. The samples (labeled with Hy3™) and the references (labeled with Hy5™) were then mixed pairwise and recruited for hybridization. The hybridization procedure was performed using Tecan’s HS 4800™ hybridization station (Tecan, Austria) followed by miRCURY LNA™ miRNA array instruction manual.

Later, the labeled specimens were hybridized by miRCURY LNA™ miRNA array 7th Gen (Exiqon, Denmark) using the manufacturer’s protocol. The quantified signals of the background were also corrected (Normexp with offset value 10) and normalized using the global Locally Weighted Scatterplot Smoothing (LOWESS) regression algorithm.

Principle component analysis (PCA) plot

In order to decrease the proportions of the bulky data and to explore sample classes arising naturally and according to the expression profile, the PCA plot was applied. The top 50 miRNAs that had major variations across all the samples were also included.

Array slide quality control using spike-ins (labeling controls) and normalization

The spike-in controls were added in various concentrations in both Hy3™ and Hy5™ labeling reactions, giving the opportunity to evaluate the labeling reaction, hybridization, and performance of array experiment in a general manner (The level correlation of the signal intensities across all the slides is shown in the Additional file 1: Figure A1).

The LOWESS also allows the correction of systematic deviations in the MA plot resulting in an intensity-dependent adjustment of the MA-data to a straight line. The positive effect of this normalization is illustrated in the MA plots for each slide, showing the plot before and after normalization (see Additional file 1: Figure A2).

Real-time PCR

Poly A polymerization, cDNA synthesis, and real-time PCR

To this end, for study the expression level of selected miRNAs in the validation group of samples, 5 µg of each isolated total RNAs was used for poly A polymerization by means of microRNA cDNA Synthesis Kit (Parsgenome, Iran). Then, 2 µg of the product was applied for the next step, i.e., cDNA synthesis using specific reverse primers for each miRNA. Finally, these miR specific cDNAs were employed as templates for quantitative SYBR Green Real-time analysis (Parsgenome, Iran). The expression level of the selected miRNA molecules was further determined employing SYBR Green Real-time PCR (ABI, USA). In addition, U6 miRNA primers were included in the mentioned kit as an internal control.

The mix used for Real-time PCR for each reaction contained SYBR Green Premix (10 µl of Ex taq, Takara, Japan), Rox reference dye (0.4 µl), forward and reverse primers (1 µl of 10 pM mix of each primer pairs), and template (2 µl of synthesized cDNA, between 0.1 to 100 ng). Nuclease free water also added to this reaction up to 20 µl total volume. The used forward primer sequences are shown in Table 1. The program employed for Real-time PCR was one cycle 95 °C -5 min, followed by 40 cycles of 95 °C- 5 s, 61 °C- 20 s, and 72 °C- 30 s, followed by melting curve analysis for specificity of each reaction.

Table 1 The forward primer sequences used for validating the selected microRNAs

Statistical analysis

The data were coded and imported into the IBM SPSS Statistics software (version 24). The quantitative variables were then compared using the non-parametric Kruskal–Wallis test for a comparison between more than two groups and the Mann–Whitney U test was employed to compare two groups. The calculations for expression comparison in LNA-array tested samples were also executed by the R/Bioconductor software, using the LIMMA package. The fold changes in the validation group of samples were calculated by livak method (2−ΔΔct). The P-value cut-off score for a significant difference was considered less than 0.05.


Descriptive statistics

The test group was composed of 14 LT patients with the mean age of 47.4 ± 7.8 years (age range 31–56) and included 9 (64.3%) males. These samples were used for LNA-array probe. More data is shown in Fig. 1 and Table 2. The validation group consisted of 20 LT patients with no acute rejection of LT and 10 patients showing acute rejection. The mean age of nonrejected ones was 49 ± 15.3 years (age range 14–69) and included 14 (70%) males. The rejected samples had mean age of 49.7 ± 17.1 years (age range 13–68) and included 6 (75%) males.

Table 2 Descriptive statistics related to first group (test group) of patients; the total number of patients participating in the test group and the number and percent of male and female participated in test group is shown in the table

PCA plot analysis

The top 50 miRNAs that had major variations across all the samples were used for PCA plot analysis. The results are depicted in Fig. 2a, demonstrating variations of the biological patterns or technical factors in the test group (n = 14). Accordingly, PC1 and PC2 describe variations related to sample groups or treatments, whereas PC3 and later PC’s might describe underlying or less variable factors like sample preparation conditions, operator, storage time, etc.

Fig. 2

(a) Matrix PCA plot. The PC1 to PC4 analysis was performed on the results of microRNA expression variation using LNA-array probe on test group of samples; due to the color key information, the most differences belong to PC1 and between samples PR and PNR. Normalized log ratio values were used for the analysis. (b) Each row represents a microRNA and each column represents a sample. The color scale illustrates the relative expression level of microRNAs. Red color represents an expression level below the reference, and green color represents expression higher than the reference genes. Pooled sample groups were used in order to reach more exact analyses and the results mentioned here is related to pooled samples. Here, PR means pool of reject samples while PNR means pool of nonrejected samples and R means one reject sample while NR means one nonrejected sample

The unsupervised PCA and the hierarchical clustering could clearly show a tight grouping of the control samples (i.e., the non-rejected), whereas the treated samples (namely, the rejected) differed both from the controls as well as from each other. A simple expression analysis by comparison did detect a subset of differentially expressed miRNAs including four up-regulated and four down-regulated miRNAs selected for further analyses.

Heat map and unsupervised hierarchical clustering

The heat map diagram shows the results of a two-way hierarchical clustering of miRNAs and the samples in the test group (n = 14). The clustering was done using the complete-linkage method together with the Euclidean distance. A small subset of the miRNAs was also excluded from the heat map. As well, the normalized log ratio values were used for the analysis (Fig. 2b).

MiRNA expression analysis for rejected versus non-rejected in LNA-arrayed samples

An expression analysis for the rejected cases versus the non-rejected ones was performed using a simple comparison method in the test group (n = 14). The miRNA profiling accordingly identified a subset of top miRNAs (where the absolute value of the log fold change was larger than 1) out of the total number of miRNAs analyzed (371) by the miRCURY LNA™ miRNA Array, differentially expressed in the rejected and the non-rejected samples. Table 3 provides a list of miRNAs with the most differential expression. The results showed that hsa-miR-122-3p, -548as-3p, -4284, and -194-5p had the highest levels of reduction amongst all the studied miRs, and hsa-miR-4511, -4449, -4633-5p, and -3158-5p were the most increased ones; respectively.

Table 3 The table shows the selected differentially expressed microRNA candidates, ranked according to the absolute value of the log fold change from the test group of samples. Also, the table shows the average Hy3 (fluorescent dye) used for detecting the expression level (number of microRNAs tested) and the fold change of each microRNA in each category (nonrejected and rejected). Finally, logFC (logarithmic fold change) of each microRNA calculated in order to understand the extent of expression of each microRNA in reject samples in comparison to nonrejected group

Selected miRNA expression level comparison between rejected and non-rejected groups of patients after LT

The expression level of the selected miRs was tested in the 2nd groups of patients (validation group) composed of rejected and non-rejected ones. The results of this selection are demonstrated in Fig. 3a for the down-regulated miRs and Fig. 3b for the up-regulated ones based on a fold change. The results also revealed that hsa-miR-548as-3p and -194-5p had significantly reduced and hsa-miR-4449 and -3158-5p had shown a significant rise in the patients with LT rejection.

Fig. 3

The expression level of selected miRs for down regulated miRs (a) and for up regulated ones (b) based on fold change. The results of LNA-array probe test on the test group of samples were further validated using qRt-PCR on the validation group of samples. The selected microRNAs had the most variation among all the tested microRNAs. The fold changes are calculated by livak method (2−ΔΔct)

Calculating sensitivity and specificity of down-regulated miRNAs between nonrejected and rejected groups of patients after LT

The sensitivity and the specificity of the down-regulated miRNAs between the non-rejected and rejected samples of patients in the validation group were calculated (Fig. 4). The area under curve (AUC), p value, and 95% confidence interval (95% CI) of each down-regulated miRNA are demonstrated in Table 4. Two of the miRNAs (i.e., has-miR194-5p and -548as-3p) had a significant difference (p < 0.05), making them candidates as promising biomarkers.

Fig. 4

Calculating the sensitivity and specificity of down-regulated microRNAs between nonrejected and rejected groups of patients after transplantation. In order to further validate the sensitivity and specificity of selected down-regulated microRNAs between nonrejected and rejected groups of patients, the ROC curve analysis performed

Table 4 Sensitivity and specificity of down and up-regulated microRNAs between nonrejected and rejected groups of patients (validation group) after transplantation. The AUC and p value of each microRNA is shown to describe the importance of each studied microRNA

Calculating sensitivity and specificity of up-regulated miRNAs between nonrejected and rejected groups of patients after LT

The sensitivity and the specificity of the up-regulated miRNAs between the non-rejected and rejected samples of patients in the validation group were calculated (Fig. 5). The AUC, p-value, and 95% CI of each up-regulated miRNA are demonstrated in Table 4. Two of the miRNAs (namely, has-miR4449 and -3158-5p) established a significant difference (p < 0.05), making them candidates as promising biomarkers.

Fig. 5

Calculating the sensitivity and specificity of up-regulated microRNAs between nonrejected and rejected groups of patients after transplantation. In order to further validate the sensitivity and specificity of selected up-regulated microRNAs between nonrejected and rejected groups of patients, the ROC curve analysis performed


The functions of miRNAs can provide valuable diagnostic measurements as well as prognostic indicators for several diseases including types of cancer [5]. A healthy liver consists of different cells producing a variety of miRNA expression profiles according to complex intrinsic and extrinsic signals. Under microbial pathogenesis, these cells can also alter the profile of miRNAs [23]. Commonly, HCC is accompanied by loss of liver and the remaining treatment would be transplantation. Although many studied have been accomplished in order to find early detecting markers on the basis of miRNAs, still no effective biomarkers have been introduced in this regard. Knowing the fact that transplantation is the final therapeutic choice for HCC, the maintenance of the grafted liver and prevention from rejection is the priority for these patients. In the present research, it is tried to evaluate the alterations in miRNA expression profile of LT patients with HCC in order to introduce probable miRNA biomarker in early detection of rejection. Although it has been suggested that some of the circulating blood miRNAs have the potential to be used as miRNA-based blood biomarkers in cancer detection [24, 25], it is noteworthy to investigate the possibility of introducing more accurate miRNAs for following rejection outcomes. In a study done by Miyaaki et al., they used the same techniques as we did, for studying the miRNA profile of liver transplanted patients with HCV/HIV co-infection and confirmed that miR-101b and miR-149 are significantly decreased in co-infected group of patients comparing to HCV-infected ones [26].

The findings revealed that miR-4449, miR-4511, miR-3158-5p, and miR-4633-5p had been up-regulated and miR-4284, miR-122-3p, miR-194-5p, and miR-548as-3p down-regulated in patients with LT rejection once compared with non-rejected ones. The results from LT patients were also in accordance with the LNA-array probe results. With reference to the results collected from the PCA plot, there was a significant difference between miRNA expression in the rejected LT patients versus the non-rejected ones.

The miR-122-3p is a liver-specific miRNA, highly expressed, comprised of more than 70% of all cloned liver miRNAs, unlike other tissues [1, 27]. Moreover, this miRNA plays a role in the liver metabolism, for instance, many genes interfering in lipid metabolism regulation can be down-regulated when using an antisense strategy to knock it down [28]. The sequence and the expression pattern of the miR-122-3p in the liver is also highly conserved [29]. In a study, HCC patients had shown a reduction in miR-122-3p expression [30]. MiR-122-3p, miR-26a, and miR-195 have been also detected as tumor suppressors molecules in the liver [31]. The fact that miR-122-3p reduces in HCC patients was accordingly certified in the present study. Other studies on liver transplant recipients with recurrent HCV and acute cellular rejection demonstrated the importance of miR-122 and miR-194 in liver transplant rejection either [24, 32]. Furthermore, the correlation between miR122-3p and -194-5p has been a subject of research studies. Therefore, this miRNA was analyzed between HCC non-rejected and rejected patients, indicating that the given miRNA could significantly (p = 0.041) decrease in the rejected LT cases.

The miR-194-5p expression had been previously detected, reporting that such an miRNA could play a role in activation of stellate cells during liver fibrogenesis [33, 34]. Another study had further profiled the expression of miR-194-5p as a marker of liver HCC metastasis by showing the overexpression of miR-194 in cancerous liver cell lines to down-regulate N-cadherin expression and to suppress migration, invasion, and metastasis [31]. An increase in the expression of miR-194-5p in HCV-infected hepatoma cells has been also well documented although the miR-194-5p expression level had significantly reduced in HCC patients [35]. In addition, Farid et al. had demonstrated that the expression levels of miR-122 could significantly correlate with miR-194 among HCC patients (r = 0.322, p = 0.007). They had also reported no significant differences between non-rejected and rejected groups of patients in the expression of miR-194-5p (p = 0.134). However, the difference in miR-122 expression level between two studied group was significant (p = 0.001) [36]. In addition, miR-194-5p had been down-regulated in the tissues involved in HCC and was able to reduce cell viability and proliferation by inducing G1 arrest and apoptosis in HCC cells [37].

Some miRNAs might be originated from repetitive elements like transposable elements (TEs) [38]. Hsa-mir-548as-3p is also a member of a large human gene family, derived from transposable elements named Made1 [39]. Approximately 69 members are found in almost all human chromosomes, especially chromosomes 6, 8, and X (30.43%). The predicted targets for hsa-miR-548as-3p are also functional molecules in many biological processes like MAPK signaling, cell cycle, p53 signaling pathway, colorectal cancer (CRC), non-small cell lung cancer, B cell receptor signaling pathway, transforming growth factor β (TGF-β) signaling pathway and renal cell carcinoma [40]. TGF-β signaling pathway is also important in post-transplant inflammation. Besides, it is considered as one of the main mediators and inducers of fibrosis. Therefore, miRNA down-regulation during the rejection period renders to overexpressed SMAD4, enhancing the signaling pathway of TGF-β [41]. In addition, it has been reported that Tg737 is a target gene for miR-548as-3p, whose down-regulation facilitates HCC cell proliferation in vivo and in vitro, followed by an increase in colony formation [42]. The present study showed that the given miRNA had significantly (p = 0.041) reduced in rejected LT patients with HCC.

The miRNA-4248 is located on chromosome 7, but there is not much data published about it. One study had used cross-mapping to detect miRNAs and had reported hsa-miR-4284 in cross-mapping events [43]. In addition, another research team had reported that hsa-miR-4284 could be the most down-regulated ones in the tissue samples of patients with ulcerative colitis. They had further claimed that miR-4284 level correlated with the disease activity by regulating CXCL5 mRNA expression 44. Another research had found that miR-4284 and -4484 could be considered as diagnostic biomarkers for diffuse large B-cell lymphoma (DLBL) [45]. The down-regulation of this miRNA was further detected in LT patients with HCC in the present study. Furthermore, the sensitivity and the specificity of analyses of the four down-regulated miRNAs showed that miR194-5p and miR548as-3p might have enough qualities (Table 4) for being biomarker candidates for rejection especially among HCC patients.

MiR-3158-5p is detected to be located on chromosome 10 and has more than 500 predicted targets. The bioinformatic studies here showed that TGF-β2 gene could be observed among its upper 95% target scores (see Additional file 1: Table A1; “”). TGF-β also plays an important role in preparing a suitable microenvironment for tumor cell growth in liver diseases. This factor causes a signaling pathway that can promote progression of HCC in two ways; first, this factor has an intrinsic activity as an autocrine or paracrine growth factor, and second, it has extrinsic activity inducing changes in microenvironments like variations in cancer-associated fibroblasts, regulatory T cells (Treg cells), and inflammatory mediators [46]. The other up-regulated miRNA is miR-4449, located on chromosome 4 with 22 predicted targets in miRDB (see Additional file 1: Table A2). Among the predicted targets, the cyclin-dependent kinase 5 regulatory subunit 2, p39 (CDK5R2) is detected as the most highly scored one. Previously, P39 had been introduced as a potential clinical prognostic marker for HCC [47]. Also, it had been detected that miR-4449 was a potential blood-based marker in multiple myeloma [48]. The present study suggested that miR-3158-5p and miR-4449 had significantly up-regulated in HCC patients with LT rejection. Even in analyzing sensitivity and specificity, these two miRNAs had AUC greater than 0.8 (Table 4). Therefore, they were proposed as a biomarker candidate for rejection in patients with HCC.

The third up-regulated miRNA in HCC patients with LT rejection was miR-4633-5p, located on chromosome 5 with more than 100 predicted targets in miRDB (see Additional file 1: Table A3). In 2015, a research group had further noted that miR-4633-5p could be significantly expressed in human epidermal growth factor receptor 2 (HER2)-positive breast carcinomas among 85 tested miRNAs versus normal population [49]. The last up-regulated studied miRNA in LT rejected patients with HCC is miR-4511, located on chromosome 15 with more than 400 predicted targets in miRDB (see Additional file 1: Table 4A). There are also few studies on the role of miRNA in LT and further studies are warranted to detect the exact role of this miR and its targets. The results of the present study showed the up-regulation of these two miRNAs in HCC patients with rejected LT.

It is of note that more than 97% of liver biopsies are hepatocytes and biliary cells, and the rest are Kupffer cells and endothelial ones. In LT rejected biopsies, some macrophages, plasma cells, and lymphocytes are also added to their population [50]. Given the fact that the selected miRNAs are expressed in the non-rejected group of patients and the majority of cells in biopsies are hepatocytes. Alterations in the expression profile of the selected miRNAs might be due to the rejection process in hepatocytes.

Putting together, noticing the fact that there is un urgent need for detecting biomarkers in diagnosis and prognosis of HCC, the findings of this study might be the first step in this pathway and seems to be preliminary keys for further confirmation of their regulatory role in rejection in LT patients with HCC. For instance, among down-regulated miRNAs, miR194-5p and miR548as-3p, both showed significant difference between rejected and nonrejected samples and the ROC curve analysis of these two miRNAs showed significant high specificity and sensitivity. Also, among up-regulated miRNAs, miR-4449 and miR-3158-5p can be considered to have important roles in causing rejection in such patients. Besides the ROC curve analysis of these two microRNAs showed significant high specificity and sensitivity. Further studies on detecting the important targets of these miRNAs, obviously would be a large finding in this field.

Although finding miRNAs as biomarkers might be very helpful for more accurate prognosis and diagnosis of HCC. This study was suffered from some limitations including: small sample size based on nature and prevalence of HCC, difficulties for follow up of patients and ethical issue related to collecting biopsy samples from normal tissue.


In summary, the identification of miRNAs and small RNA species seems to represent only the tip of the iceberg. In addition, the prediction of an individual miRNA, its target, and even its function is still one of the biggest challenges in research. Therefore, based on these results, two down-regulated (miR194-5p and miR584as-3p) and two up-regulated (namely, miR-4449 and miR-3158-5p) miRNAs were introduced as candidate biomarkers for LT rejection in HCC patients which need to be confirmed by further studies. Even though the biological function of miRNAs and results of this research encourage researches to study more this subject.



alanine transaminase


aspartate transaminase




hepatocellular carcinoma




3′-untranslated regions


locked nucleic acid-array


formalin-fixed paraffin-embedded


transposable elements


transforming growth factor β


large B-cell lymphoma

Treg cells:

regulatory T cells


colorectal cancer


cyclin-dependent kinase 5 regulatory subunit 2, p39


human epidermal growth factor receptor 2


  1. 1.

    Chen XM. MicroRNA signatures in liver diseases. World J Gastroenterol. 2009;15(14):1665–72.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Huang S, He X. The role of microRNAs in liver cancer progression. Br J Cancer. 2011;104(2):235–40.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Dragani TA. Risk of HCC: Genetic heterogeneity and complex genetics. J Hepatol. 2010;52(2):252–7.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Jemal A, Center MM, DeSantis C, Ward EM. Global patterns of cancer incidence and mortality rates and trends. Cancer Epidemiol Biomark Prev. 2010;19(8):1893–907.

    Article  Google Scholar 

  5. 5.

    Castoldi M, Schmidt S, Benes V, et al. A sensitive array for microRNA expression profiling ( miChip ) based on locked nucleic acids ( LNA ) A sensitive array for microRNA expression profiling ( miChip ) based on locked nucleic acids ( LNA ). RNA Soc. 2006;12:913–20.

    CAS  Article  Google Scholar 

  6. 6.

    Barad O, Meiri E, Avniel A, et al. MicroRNA expression detected by oligonucleotide microarrays: system establishment and expression profiling in human tissues. Methods. 2004;14:2486–94.

    CAS  Article  Google Scholar 

  7. 7.

    Song Q, Song J, Wang Q, et al. miR-548d-3p/TP53BP2 axis regulates the proliferation and apoptosis of breast cancer cells. Cancer Med. 2016;5(2):315–24.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Lewis AP, Jopling CL. Regulation and biological function of the liver-specific miR-122. Biochem Soc Trans. 2010;38(6):1553–7.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Hayes C, Chayama K. MicroRNAs as biomarkers for liver disease and hepatocellular carcinoma. Int J Mol Sci. 2016;17(3):280.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Murakami Y, Yasuda T, Saigo K, et al. Comprehensive analysis of microRNA expression patterns in hepatocellular carcinoma and non-tumorous tissues. Oncogene. 2006;25(17):2537–45.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Yang N, Ekanem NR, Sakyi CA, Ray SD. Hepatocellular carcinoma and microRNA: new perspectives on therapeutics and diagnostics. Adv Drug Deliv Rev. 2015;81:62–74.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Gramantieri L, Ferracin M, Fornari F, et al. Cyclin G1 is a target of miR-122a, a MicroRNA frequently down-regulated in human hepatocellular carcinoma. Cancer Res. 2007;67(13):6092–9.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Burchard J, Zhang C, Liu AM, et al. MicroRNA-122 as a regulator of mitochondrial metabolic gene network in hepatocellular carcinoma. Mol Syst Biol. 2010;6:402.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Mizuguchi Y, Mishima T, Yokomuro S, et al. Sequencing and Bioinformatics-Based Analyses of the microRNA Transcriptome in Hepatitis B–Related Hepatocellular Carcinoma. Guan X, ed. PLoS One. 2011;6(1):e15304.

  15. 15.

    Wang Y, Lee ATC, Ma JZI, et al. Profiling microRNA expression in hepatocellular carcinoma reveals microRNA-224 up-regulation and apoptosis inhibitor-5 as a microRNA-224-specific target. J Biol Chem. 2008;283(19):13205–15.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Sui W, Lin H, Peng W, et al. Molecular dysfunctions in acute rejection after renal transplantation revealed by integrated analysis of transcription factor, microRNA and long noncoding RNA. Genomics. 2013;102(4):310–22.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Shaked A, Chang BL, Barnes MR, et al. An ectopically expressed serum miRNA signature is prognostic, diagnostic, and biologically related to liver allograft rejection. Hepatology. 2017;65(1):269–80.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Morita M, Chen J, Fujino M, et al. Identification of microRNAs involved in acute rejection and spontaneous tolerance in murine hepatic allografts. Sci Rep. 2014.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Castoldi M, V B, MW H, MU M. miChip: a microarray platform for expression profiling of microRNAs based on locked nucleic acid (LNA) oligonucleotide capture probes. Methods. 2007.

  20. 20.

    Castoldi M, S S, V B, MW H, MU M. miChip: an array-based method for microRNA expression profiling using locked nucleic acid capture probes. Nat Protoc. 2008;3(2). doi:

  21. 21.

    Li W, K R. MicroRNA detection by microarray. Anal Bioanal Chem. 2009;394(4). doi:

  22. 22.

    American Joint Committee on Cancer. Accessed November 22, 2020.

  23. 23.

    Chen X-M, Splinter PL, O’Hara SP, LaRusso NF. A cellular micro-RNA, let-7i, regulates Toll-like receptor 4 expression and contributes to cholangiocyte immune responses against Cryptosporidium parvum infection. J Biol Chem. 2007;282(39):28929–38.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Asaoka T, Hernandez D, Tryphonopoulos P, et al. Clinical significance of intragraft miR-122 and -155 expression after liver transplantation. Hepatol Res. 2015;45(8):898–905.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Jin Y, Wong YS, Goh BKP, et al. Circulating microRNAs as Potential Diagnostic and Prognostic Biomarkers in Hepatocellular Carcinoma. Sci Rep. 2019;9(1):1–12.

    CAS  Article  Google Scholar 

  26. 26.

    Miyaaki H, Takatsuki M, Ichikawa T, et al. Intrahepatic microRNA profile of liver transplant recipients with hepatitis C virus co-infected with human immunodeficiency virus. Ann Transplant. 2017;22:701–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Chang J, Nicolas E, Marks D, et al. miR-122, a mammalian liver-specific microRNA, is processed from hcr mRNA and may downregulate the high affinity cationic amino acid transporter CAT-1. RNA Biol. 2004;1(2):106–13.

    CAS  Article  Google Scholar 

  28. 28.

    Arora A, Simpson DA. Individual mRNA expression profiles reveal the effects of specific microRNAs. Genome Biol. 2008;9(5):R82.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Wienholds E, Kloosterman WP, Miska E, et al. MicroRNA expression in Zebrafish embryonic development. Science (80-). 2005;309(5732):310–1.

    CAS  Article  Google Scholar 

  30. 30.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Meng Z, Fu X, Chen X, et al. miR-194 is a marker of hepatic epithelial cells and suppresses metastasis of liver cancer cells in mice. Hepatology. 2010;52(6):2148–57.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Farid WRR, Pan Q, van der Meer AJP, et al. Hepatocyte-derived microRNAs as serum biomarkers of hepatic injury and rejection after liver transplantation. Liver Transpl. 2012;18(3):290–7.

    Article  PubMed  Google Scholar 

  33. 33.

    Raymond CK, Roberts BS, Garrett-Engele P, Lim LP, Johnson JM. Simple, quantitative primer-extension PCR assay for direct monitoring of microRNAs and short-interfering RNAs. RNA. 2005;11(11):1737–44.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Venugopal SK, Jiang J, Kim T-H, et al. Liver fibrosis causes downregulation of miRNA-150 and miRNA-194 in hepatic stellate cells, and their overexpression causes decreased stellate cell activation. Am J Physiol Gastrointest Liver Physiol. 2010;298(1):G101–6.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Jung HS, Seo Y-R, Yang YM, et al. Gα12gep oncogene inhibits FOXO1 in hepatocellular carcinoma as a consequence of miR-135b and miR-194 dysregulation. Cell Signal. 2014;26(7):1456–65.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Farid WRR, Pan Q, Van Der Meer AJP, et al. Hepatocyte-derived microRNAs as serum biomarkers of hepatic injury and rejection after liver transplantation. Liver Transplant. 2012;18(3):290–7.

    Article  Google Scholar 

  37. 37.

    Zhao Y, Li F, Zhang X, et al. MicroRNA-194 acts as a prognostic marker and inhibits proliferation in hepatocellular carcinoma by targeting MAP4K4. Int J Clin Exp Pathol. 2015;8(10):12446–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Yuan Z, Sun X, Liu H, Xie J. MicroRNA genes derived from repetitive elements and expanded by segmental duplication events in mammalian genomes. Tjwa M, ed. PLoS One. 2011;6(3):e17666.

  39. 39.

    Piriyapongsa J, Jordan IK. A family of human microRNA genes from miniature inverted-repeat transposable elements. Christoffels A, ed. PLoS One. 2007;2(2):e203. doi:

  40. 40.

    Liang T, Guo L, Liu C. Genome-wide analysis of mir-548 gene family reveals evolutionary and functional implications. J Biomed Biotechnol. 2012;2012:1–8.

    CAS  Article  Google Scholar 

  41. 41.

    Hamdorf M, Kawakita S, Everly M. The potential of MicroRNAs as novel biomarkers for transplant rejection. J Immunol Res. 2017;2017:1–12.

    CAS  Article  Google Scholar 

  42. 42.

    Zhao G, Wang T, Huang Q-K, et al. MicroRNA-548a-5p promotes proliferation and inhibits apoptosis in hepatocellular carcinoma cells by targeting Tg737. World J Gastroenterol. 2016;22(23):5364.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Guo L, Liang T, Gu W, Xu Y, Bai Y, Lu Z. Cross-mapping events in miRNAs reveal potential miRNA-Mimics and evolutionary implications. PLoS ONE. 2011.

    Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Koukos G, Polytarchou C, Kaplan JL, et al. A MicroRNA signature in pediatric ulcerative colitis. Inflamm Bowel Dis. 2015;21(5):996–1005.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Tamaddon G, Geramizadeh B, Karimi MH, Mowla SJ, Abroun S. miR-4284 and miR-4484 as putative biomarkers for diffuse large B-Cell lymphoma. Iran J Med Sci. 2016;41(4):334–9.

    PubMed  Google Scholar 

  46. 46.

    Giannelli G, Villa E, Lahn M. Transforming growth factor-β as a therapeutic target in hepatocellular carcinoma. Cancer Res. 2014;74(7):1890–4.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Lu J-W, Chang J-G, Yeh K-T, Chen R-M, Tsai JJP, Hu R-M. Decreased expression of p39 is associated with a poor prognosis in human hepatocellular carcinoma. Med Oncol. 2011;28(S1):239–45.

    CAS  Article  Google Scholar 

  48. 48.

    Shen X, Ye Y, Qi J, et al. Identification of a novel microRNA, miR-4449, as a potential blood based marker in multiple myeloma. Clin Chem Lab Med. 2017;55(5):748–54.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Vos S, Vesuna F, Raman V, van Diest PJ, van der Groep P. miRNA expression patterns in normal breast tissue and invasive breast cancers of BRCA1 and BRCA2 germ-line mutation carriers. Oncotarget. 2015;6(31):32115–37.

    Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Ross MH, Kaye GI, Pawlina W. Histology : a text and atlas with cell and molecular biology. Lippincott Williams & Wilkins; 2003. Accessed April 9, 2018.

Download references


The authors wish to thank Mr. H. Argasi at the Research Consultation Center (RCC) of Shiraz University of Medical Sciences for his invaluable assistance in editing this manuscript.


This work was supported by grants from Transplant Research Center, Shiraz University of Medical Sciences, with grant No. 91/181. Funding sources had no influence over study design, data collection, analysis or interpretation or manuscript preparation and submission.

Author information




AA performed the majority of experiments and analyzed the data; AA, KM, YR performed the molecular investigations; AA, KM, YR and MJ designed and coordinated the research; AA, KM, YR and MJ wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mohammad Hossein Karimi.

Ethics declarations

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Consent for publication

Written informed consent for publication was obtained.

Ethics approval and consent to participate

All applied protocols in this study, were in complete agreement with the Helsinki Declaration and its later amendment. The local Ethics committee of Shiraz University of Medical Sciences approved every step of this study. After explaining the study objectives, written informed consent was obtained from each patient.

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

. Extra plots and data of evaluated microRNA targets.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Afshari, A., Yaghobi, R., Karimi, M.H. et al. Alterations in MicroRNA gene expression profile in liver transplant patients with hepatocellular carcinoma. BMC Gastroenterol 21, 262 (2021).

Download citation


  • LNA-array
  • microRNA
  • Liver
  • Transplantation
  • Hepatocellular carcinoma
  • RT-PCR