Combined proteomic/transcriptomic signature of recurrence post-liver transplantation for hepatocellular carcinoma beyond Milan

Background and aims Liver transplantation (LT) can be offered to patients with Hepatocellular carcinoma (HCC) beyond Milan criteria. However, there are currently limited molecular markers on HCC explant histology to predict recurrence, which arises in up to 20% of LT recipients. The goal of our study was to derive a combined proteomic/transcriptomic signature on HCC explant predictive of recurrence post-transplant using unbiased, high-throughput approaches. Methods Patients who received a LT for HCC beyond Milan criteria in the context of hepatitis B cirrhosis were identified. Tumor explants from patients with post-transplant HCC recurrence (N = 7) versus those without recurrence (N = 4) were analyzed by mass spectrometry and gene expression array. Univariate analysis was used to generate a combined proteomic/transcriptomic signature linked to recurrence. Significantly predictive genes and proteins were verified and internally validated by immunoblotting and immunohistochemistry. Results Seventy-nine proteins and 636 genes were significantly differentially expressed in HCC tumors with subsequent recurrence (p < 0.05). Univariate survival analysis identified Aldehyde Dehydrogenase 1 Family Member A1 (ALDH1A1) gene (HR = 0.084, 95%CI 0.01–0.68, p = 0.0152), ALDH1A1 protein (HR = 0.039, 95%CI 0.16–0.91, p = 0.03), Galectin 3 Binding Protein (LGALS3BP) gene (HR = 7.14, 95%CI 1.20–432.96, p = 0.03), LGALS3BP protein (HR = 2.6, 95%CI 1.1–6.1, p = 0.036), Galectin 3 (LGALS3) gene (HR = 2.89, 95%CI 1.01–8.3, p = 0.049) and LGALS3 protein (HR = 2.6, 95%CI 1.2–5.5, p = 0.015) as key dysregulated analytes in recurrent HCC. In concordance with our proteome findings, HCC recurrence was linked to decreased ALDH1A1 and increased LGALS3 protein expression by Western Blot. LGALS3BP protein expression was validated in 29 independent HCC samples. Conclusions Significantly increased LGALS3 and LGALS3BP gene and protein expression on explant were associated with post-transplant recurrence, whereas increased ALDH1A1 was associated with absence of recurrence in patients transplanted for HCC beyond Milan criteria. This combined proteomic/transcriptomic signature could help in predicting HCC recurrence risk and guide post-transplant surveillance. Supplementary Information The online version contains supplementary material available at 10.1186/s12014-021-09333-x.


Introduction
Hepatocellular carcinoma (HCC) is estimated to represent the fourth most common cause of cancer-related deaths globally [1]. Selected patients with early-stage HCC are candidates for liver transplantation (LT) as a curative therapy. HCC has even become the leading indication for LT worldwide in recent years [1,2]. Currently, the Milan criteria (one lesion less than 5 cm in diameter or three or fewer lesions, all below 3 cm) remain the gold standard for defining the optimum tumor burden that determines transplant eligibility [3]. Nonetheless, it has been recognized that expanded selection criteria can afford patients excellent post-LT outcomes even when traditional morphologic tumor standards are exceeded [4]. Criteria have evolved by expanding the allowable tumor size and number to the incorporation of alphafetoprotein (AFP) and histologic differentiation [4][5][6][7][8][9]. The safe expansion of selection criteria for LT in HCC hinges on an improved understanding of tumor biology through surrogate markers [1]. Consequently, such expanded criteria may allow extending the use of LT as a curative-intent treatment option to a greater number of patients, especially considering that HCC biology and progression are patient-dependent.
Numerous studies have tried to identify molecular biomarkers that predict HCC recurrence after LT or resection [10][11][12][13][14][15][16][17][18]. In this context, proteomic analysis of liver tissue from HCC patients has recently identified some proteins such as melanoma-associated antigen genes [19] and cytokeratin-19 [20]; their expression correlated with early recurrence of HCC post-hepatectomy and allowed stratification into various subtypes with distinct clinical outcomes [11]. Whether such proteomic signatures can stratify patients in a diverse HCC LT population that exceeds Milan criteria is unclear. Identifying novel molecular signatures may enable the selection of patients exceeding traditional standards who may still benefit from LT. Also, signature proteins related to distinct tumor subtypes may potentially serve as actionable targets for individualized prevention and therapy.
We hypothesized that genes and proteins shared among HCC tumors on explant could be used to predict recurrence. Therefore, the goal of the present study was to define a combined proteomic/transcriptomic signature predictive of recurrence post-transplant in patients transplanted for HCC beyond Milan criteria, based on the molecular profiling of the dominant tumour on explant.

Patient population
Patients who underwent LT for hepatitis B-induced HCC beyond Milan criteria between 2004 and 2015 were included in our study. This protocol was approved by the University Health Network (UHN) Institutional Review Board (REB#15-9989). Those with recurrence were matched for age and sex to those without recurrence by at least 2 years post-transplant. The first set of patients (n = 11) had snap-frozen HCC explant samples in the biobank, meeting inclusion criteria, and was therefore eligible for proteome and transcriptome analysis. An additional set of 29 samples meeting the above inclusion criteria were from patients with only archived formalinfixed, paraffin-embedded samples available, on which only immunohistochemistry (IHC) validation was feasible. Characteristics of the tumors on explant, including size and number of tumors, microvascular invasion and associated AFP at the time of transplant were documented. Patients could have received bridging therapy but needed to have a viable tumor on the explant for the purposes of this analysis.

Hepatocellular carcinoma sample preparation for mass spectrometry
Frozen HCC tissue samples obtained from the UHN biobank were processed and analyzed in a blinded fashion. These represented the largest, viable tumor in the explant. Rapigest detergent (0.2%) was added to each piece of tissue (~ 100 mg in weight). Samples were then homogenized using microbeads in a tissue dissociator. Ammonium bicarbonate was added, and samples were subjected to sonication (10 s, 3 times) on ice. Samples were then centrifuged at 15,000g at 4 °C for 20 min. Supernatant was collected and vortexed. Total protein concentration was determined using the micro-BCA protein assay kit (Thermo). Each sample was then normalized to 250 µg of total protein. Samples were denatured at 80 °C for 15 min, followed by reduction in 10 mM Conclusions: Significantly increased LGALS3 and LGALS3BP gene and protein expression on explant were associated with post-transplant recurrence, whereas increased ALDH1A1 was associated with absence of recurrence in patients transplanted for HCC beyond Milan criteria. This combined proteomic/transcriptomic signature could help in predicting HCC recurrence risk and guide post-transplant surveillance.
Keywords: Hepatocellular carcinoma, Liver transplantation, Milan criteria, Proteome, Transcriptome, Combined signature, Tumor explant, Predictors of recurrence dithiothreitol (final concentration) for 15 min at 65 °C and alkylation in 20 mM (final concentration) iodoacetamide for 40 min in the dark at room temperature. Finally, trypsin (Promega) was added at 1:50 w/w and incubated overnight at 37 °C. Digested samples were acidified with trifluoroacetic acid (1% v/v) and vortexed for 1 min, then left at room temperature for 5 min. Samples were centrifuged at 15,000g for 10 min to remove Rapigest, then transferred into the new tubes and frozen at -20 °C until further analysis [21,22]. Following digestion and removal of Rapigest, strong cation exchange (SCX) chromatography and fractionation were performed on an HPLC system (Agilent 1100) using a 60-min two-step gradient. The resulting fractions corresponding to chromatographic peaks of eluting peptides were pooled into 10 fractions. Peptides from each fraction were extracted, desalted and diluted to 41 μL with 0.1% v/v formic acid in MS-grade pure water.

Protein identification and quantification
For protein identification, the RAW files of each MS run were generated by XCalibur software v3.0.63 (Thermo Scientific). Raw data were analyzed by MaxQuant software (version 1.5.3.28) and searched against the human Uniprot FASTA database (HUMAN5640_sProt-072016, update of July 20, 2016). Proteins and peptides were identified with a false discovery rate of 1%. A minimum length of 6 amino acids was selected. The false positive rate was determined using reversed mode. Trypsin/P was selected as digestion enzyme, and a maximum of 2 missed cleavages was enabled. While cysteine carbamidomethylation was selected as a fixed modification, methionine oxidation and N-terminal acetylation were set as variable modifications. The initial peptide tolerance against a 'human-first-search' database was set to 20 ppm. The main search peptide mass tolerance was 40 ppm, and the fragment mass MS/MS tolerance was 0.5 Da. Matching between runs was selected. Normalized label-free quantification (LFQ) of proteins was derived from extracted ion current information from razor and unique peptides with a minimum ratio count of 1. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [24] partner repository with the dataset identifier PXD022881 (Reviewer account details: Username: reviewer_pxd022881@ebi. ac.uk; Password: 2GmXJTJo).
Proteomics data were analyzed using Perseus software (version 1.5.2.6). Reverse hits and contaminants were manually checked and removed. Distribution of log2transformed LFQ intensity values of all quantified proteins was examined for each sample (Additional file 1: Fig.  S1A). Following exclusion of proteins identified in < 50% of the samples, we subjected the zero-value intensities to imputation (assuming that low abundance values were missing), keeping a normal distribution, with a downshift of 1.8 standard deviations, and a width of 0.5 for each sample (Additional file 1: Fig. S1B). After imputation, we determined the differentially expressed proteins between recurrent and non-recurrent HCC samples by comparing their mean log2-transformed LFQ intensities using the two-tailed independent t-test (P < 0.05), followed by Benjamini-Hochberg adjustment. Principal component analysis was performed in Perseus. Two components that explained the most variability in samples were selected. Finally, pathway analysis was performed on significantly differentially expressed proteins, using pathDIP (http:// ophid. utoro nto. ca/ pathD IP/) [25].

Gene expression analysis
RNA was extracted from snap-frozen HCC explant tumor specimens and purified using RNeasy Mini Kit (Qiagen; Hilden, Germany). RNA quality was verified by Nanodrop spectrophotometer (VWR; Radnor, PA) and Bioanalyzer (Agilent; Santa Clara, CA). Gene expression microarray analysis was performed using the Affymetrix Human Gene 2.0 ST platform (Thermo Fisher, Waltham, MA). The initial cohort included eleven samples, but two of them did not pass the quality control filter required for microarray processing, therefore the remaining nine were used for further analysis. We compared gene expression of the explant HCC specimens between those who developed recurrence (n = 6) versus those who did not develop recurrence (n = 3).
Microarray data were pre-processed (background subtraction and quantile normalization) by Robust Multi-array Average (RMA) utilizing the oligo package in R version 3.6 [26]. Gene annotation was done using pd.hugene.2.0.st annotation files [27]. Statistical analysis was performed with limma package [28]. After having fit the model with lmFit function (linear model), the differential gene expression was calculated using eBayes function (moderated t-test, p-value, B stats). A gene was considered differentially expressed between the two groups if p-value < 0.05. In addition, filtering criteria by fold change (FC) was selected: FC ≥ 2 (upregulation) or FC ≤ 0.5 (downregulation).
For identification of protein-protein interactions and pathways in explant HCC that were predictive of HCC recurrence, we analyzed the common significant genes/proteins separately for the two groups using STRING software, version 10.5 (https:// string-db. org). The transcriptomic data have been deposited to Gene Expression Ominibus (GEO) repository with the dataset identifier GSE164368.

Protein expression analysis by Western Blot
For validation purposes, we selected only those proteins in common with corresponding differentially expressed genes, namely Galectin 3 (LGALS3), Galectin 3 Binding Protein (LGALS3BP), and Aldehyde Dehydrogenase 1 Family Member A1 (ALDH1A1), that were significantly associated with the time of recurrence by univariate analysis, as described in the statistical analysis section.
We measured protein expression by Western blot in the same human HCC tissue samples used for proteome analysis. Protein concentration was determined using a micro-BCA protein assay kit (Thermo). To verify changes in LGALS3, LGALS3B and ALDH1A1 protein expression, 10 µg of total protein were loaded onto 10% acrylamide gels, separated by SDS-PAGE, and transferred to a PVDF membrane (Millipore). Membranes were then blocked with 5% milk and incubated with mouse monoclonal anti-LGALS3 (1:4000; ab2785, Abcam, previously characterized [29]) or rabbit polyclonal anti-ALDH1A1 (1:2000; ab227948, Abcam [30]). Incubation with rabbit polyclonal anti-LGALS3BP (HPA000554, Atlas Antibodies) did not result in a specific band. Control for protein loading was performed by reblotting membranes using a mouse monoclonal antibody for GAPDH (CB1001, Sigma). The secondary antibodies were HRP-conjugated anti-rabbit (A0545, Sigma) and anti-mouse (P0447, Dako). Following detection in a Gel-Imaging System (Bio-rad), bands were quantified by densitometry using Image J software.

Immunohistochemistry
Protein expression of LGALS3BP was examined by immunohistochemistry on the first set of patients (n = 11), with snap-frozen samples in the biobank meeting inclusion criteria (Table 1). LGALS3BP protein expression was also studied in an additional, independent set of 29 patients (Additional file 8: Table S6). Sections (4 µm) from formalin-fixed paraffin-embedded samples were treated for antigen retrieval following manufacturer's instructions. Rabbit anti-LGALS3BP (Human Protein Atlas, HPA000554) was used at a dilution of 1:500 as primary antibody. After incubation with an HRP-conjugated anti-rabbit secondary antibody, the antibody-HRP complex was visualized with hydrogen peroxide substrate and 3,3′-diaminobenzidine tetrahydrochloride (DAB) chromogen. Slides were then counterstained with Harris Hematoxylin. Entire slides were then digitally scanned in an Aperio ScanScope CS scanning system and analyzed by Aperio Image Scope Viewer software (Leica Biosystems Imaging, Inc, CA) using the Positive Pixel Count v9 algorithm. The intensity of the strongly positive pixels was quantified in each area and normalized to the µm 2 of area analyzed. The strongly positive intensity/µm 2 across all areas were then averaged to determine LGALS3BP protein expression in each HCC sample.

Examination of top candidates in an independent HCC dataset
Several molecules that significantly differentiated recurrent from non-recurrent samples at both gene and protein level (ALDH1A1, LGALS3, LGALS3B) were selected for further investigation. In order to examine their prognostic potential, we analyzed HCC tumors from The Cancer Genome Atlas (TCGA) using KMplotter, a webbased tool that enables survival analysis across multiple cancers and datasets. These tumors represented hepatectomy specimens and not HCC explant samples. Patient samples were split into two groups according to software cut-off for ALDH1A1 (RNAseq probe#216), LGALS3 (RNAseq probe#3958), and LGALS3BP (RNAseq probe#3959). We ran multivariate survival analysis based on the high versus low expression of each of the three genes in tumors. The two groups were compared by a Kaplan-Meier survival plot, and the hazard ratio with 95% confidence intervals and log-rank p-value were calculated.

Statistical analysis
Distributions of Western Blot ratios (LGALS3/GAPDH, ALDH1A1/GAPDH, and LGALS3/ALDH1A1) and LGALS3BP positive staining values were examined using the Shapiro-Wilks normality test. We assessed differences between groups using the independent t-test for variables following a normal distribution, and the Wilcoxon-Mann-Whitney non-parametric test for variables not following a normal distribution. Pearson correlation coefficients were calculated between Western Blot and log2-transformed LFQ protein intensity values of a particular protein, and between log2-transformed LFQ protein intensity values of different proteins of interest, using GraphPad Prism 8.0. Immunohistochemistry results were analyzed using unpaired t-test. P < 0.05 was considered significant. Data are reported as median values with interquartile range. We investigated the association of selected significant genes, proteins or clinical characteristics with the time to recurrence, at the univariate level, with the coxph function of the survival package [31] in R version 3.6. Using the same package, we plotted Kaplan-Meier survival curves, the statistical difference between the curves being calculated with survdiff function.

Proteome analysis of recurrent and non-recurrent HCC cases
Patient characteristics are summarized in Table 1. The study included 4 HCC cases with no recurrence and 7 cases with recurrence post-LT. Tumors showed moderate differentiation and the number of lesions was comparable in the two groups (n = 4 for non-recurrent, and n = 7 for recurrent). The non-recurrent group showed a lower frequency of microvascular invasion (25% of patients) than the recurrent group (57% of patients). Overall survival differed between the two groups, with a median of 6 years in non-recurrent patients and 3 years in recurrent patients. We first analyzed the proteome of HCC explants (Fig. 1A). We identified and quantified a total of 6382 proteins (FDR < 0.01) from 11 HCC samples, using unbiased LC-MS/MS followed by label-free quantification (Fig. 1B). After removal of false-positive (reverse) hits and contaminants, 6277 proteins were further analyzed. In both study groups, the total numbers of identified proteins per sample were comparable and evenly distributed (Additional file 5: Table S1, Additional file 1: Fig. S1A, B). We next excluded from further analysis those proteins identified in < 6 samples, which resulted in 4505 proteins identified and quantified in at least 6 samples. The median number of proteins identified per sample was 4073, after filtering (Additional file 5: Table S1). Our 4505 proteins demonstrated bell-shaped distribution (Additional file 1: Fig. S1A) in each sample. We thus imputed missing protein values from a normal distribution and assumed low abundance (Additional file 1: Fig. S1B). Principal component analysis using all 4505 proteins illustrated that while most recurrent and non-recurrent samples were correctly separated using protein expression data, sample 10 was the most distinct sample (Additional file 1: Fig. S1C). We finally identified 79 proteins that were significantly differentially expressed between recurrent and non-recurrent samples (p < 0.05), with 44 proteins increased and 35 proteins decreased in recurrent samples (Additional file 1: Fig. 1C, Additional file 6: Table S2). We focused on proteins LGALS3BP (recurrent/non-recurrent fold change = 5.8, p = 0.025), member of the beta-galactosidebinding proteins that modulate cell-cell and cell-matrix interactions and LGALS3, a galectin/carbohydrate binding protein (recurrent/non-recurrent fold change = 2.3, p = 0.042), among the proteins significantly increased in recurrent samples. We also focused on ALDH1A1 (recurrent/non-recurrent = 0.2, p = 0.048), which was significantly decreased in recurrent samples. ALDH1A1 is an enzyme member of the aldehyde dehydrogenase family involved in alcohol metabolism. These proteins were selected because of their concomitant significant upregulation at the gene expression level [32][33][34]. Additionally, we noted that these proteins were identified confidently, with many unique peptides, and their expression levels correlated. Levels of ALDH1A1 negatively and significantly correlated with the levels of LGALS3 (Additional file 2: Fig. S2A), while the protein levels of LGALS3 and LGALS3BP showed a positive and significant correlation with each other (Additional file 2: Fig. S2B). Among the 79 proteins differentially expressed in recurrent vs nonrecurrent cases, 21 pathways were significantly enriched (adj. p < 0.05) (Additional file 7: Table S3). Most proteins, including LGALS3BP and LGALS3, belonged to the "immune system" and "innate immune system" based on pathway analysis.

Transcriptome analysis of recurrent and non-recurrent HCC cases
A total of 636 genes were significantly differentially expressed between the same recurrent compared to nonrecurrent HCC samples. We first examined those genes that were also significantly increased at the protein level in recurrent HCC. Sixty-eight upregulated genes and 44 upregulated proteins were identified in HCC explants associated with post-transplant recurrence, with five molecules in alignment with our proteome findings (Fig. 2). The common genes/proteins increased in recurrent HCC were significantly associated with cancer signaling pathways including PI3K-Akt signaling pathway (p = 0.002) and Ras signaling pathway (p = 0.010) (Additional file 8: Table S4), as well as TGF-β signaling pathway (p = 0.0099) and extracellular matrix (ECM) receptor interaction (p = 0.0099). Metabolic pathways  Table S5, Fig. 3).

ALDH1A1, LGALS3, and LGALS3BP are predictors of HCC recurrence
Molecules that significantly differentiated recurrent and non-recurrent samples at both gene and protein level, together with the selected clinical variables (Table 2) were investigated at univariate level. High levels of the ALDH1A1 gene and protein were negatively associated with the HCC recurrence. On the other hand, both LGALS3BP and LGALS3 gene and protein levels were positively associated with HCC recurrence (Fig. 4). Except for age, no other clinical variable, including the RETREAT score that was developed for the HCC within Milan patients, was predictive for HCC recurrence. We noted that the RETREAT score was not significantly different between patients with recurrent HCC vs nonrecurrent (p-value = 0.068). In addition, we plotted survival curves using the mean value as cut-off (Additional file 3: Fig. S3). The results show that ALDH1A1 protein, LGALS3 protein and LGALS3BP gene were still significantly associated with HCC non-recurrence. Due to small sample size, we did not perform multivariable analysis.

Verification and validation of ALDH1A1, LGALS3 and LGALS3BP changes in recurrent and non-recurrent HCC
We employed Western Blot to confirm the observed changes in the protein expression of ALDH1A1, LGALS3 and LGALS3BP. Unfortunately, the antibody against LGALS3BP did not generate a single clear band, forcing us to focus on the first two proteins. In keeping with our proteomics findings, ALDH1A1 protein expression was numerically increased in non-recurrent samples (P = 0.085; Fig. 5A, B), while LGALS3 expression was significantly increased in recurrent samples (P = 0.027; Fig. 5A, C). Reassuringly, there was a strong, significant, and direct correlation between protein intensity levels of ALDH1A1 and LGALS3 measured by mass spectrometry, and their corresponding protein expression by Western Blot (R = 0.837, P = 0.001 for ALDH1A1; R = 0.745, P = 0.008 for LGALS3; Fig. 5D). Moreover, mass spectrometry intensity levels of ALDH1A1 negatively and significantly correlated with mass spectrometry levels of LGALS3 (R = − 0.838, P = 0.001; Additional file 2:    S2A). Since ALDH1A1 and LGALS3 were differentially expressed in opposite directions between recurrent and non-recurrent samples, we evaluated the LGALS3/ ALDH1A1 ratio. Importantly, we observed a > 20-fold significant increase in the LGALS3/ALDH1A1 Western Blot ratio in recurrent, as compared to non-recurrent samples (P = 0.012; Fig. 5E). We next employed immunohistochemistry to confirm the observed changes in the protein expression of LGALS3BP. In concordance with our observations at the proteome and transcriptome level, LGALS3BP positive staining was numerically increased in recurrent HCC tumor explants, compared to non-recurrent, when examining the 11 samples from our discovery cohort (Fig. 6A). Reassuringly, increased LGALS3BP staining was also observed among recurrent cases in an independent, internal cohort of 29 LT recipients with HCC beyond Milan criteria (Fig. 6B).

Examination of top candidates in an independent HCC dataset
The impact of ALDH1A1, LGALS3 and LGALS3BP gene expression on overall survival was analyzed using the TCGA HCC dataset of 364 patients with HCC. ALDH1A1 was significantly predictive of overall survival in HCC patients, with a hazard ratio of 0.66 (95%CI 0.46-0.95, logrank p = 0.024) being associated with high expression level (Additional file 4: Fig. S4A). In contrast, increased LGALS3 expression (Additional file 4: Fig. S4B) was associated with an increased likelihood of death, with a hazard ratio of 1.59 (95%CI 1.09-2.3, logrank p = 0.014). Similarly, high expression of LGALS3BP (Additional file 4: Fig. S4C) was associated with increased likelihood of death, with a hazard ratio of 1.36 (95%CI 0.94-1.97), although below significance threshold (logrank p = 0.1). These findings in the independent dataset suggest that the expression of ALDH1A1, LGALS3 and LGALS3BP may be indicative of the aggressiveness of HCC, and thus its inherent propensity to recur post-transplant.

Discussion
We have delineated a set of proteins/genes on liver explant predictive of recurrence post-transplant, in patients transplanted for HCC beyond Milan. We found that expression levels of 3 specific proteins and their  transcripts (ALDH1A1, LGALS3 and LGALS3BP) distinguish those patients with recurrence post-transplant from those who did not recur. These findings were orthogonally validated in an independent cohort of recipients transplanted for HBV-HCC beyond Milan. These analytes represent a potential prognostic immunohistochemical signature in HCC beyond Milan criteria. Our combined proteomic/transcriptomic signature adds to the current knowledge on predicting recurrence post-LT in HCC beyond Milan criteria. Current predictors of HCC recurrence after LT have been limited to tumor burden on the explant, presence of vascular invasion, and elevated AFP, together comprising scores such as the RETREAT score, used to determine the frequency of surveillance for post-LT recurrence [35]. However, this cannot be used in pre-transplant decision making, as microvascular invasion is only confirmed by posttransplant explant histological assessment. This denotes the limitations of the current image-based only selection algorithm of HCC for LT. Additionally, the ability to predict recurrence is critical for a few reasons: first, those beyond Milan HCC with favourable molecular signatures could be considered for LT. Secondly, the risk-tobenefit ratio might not be in favour of transplant in those with poor molecular features. Those with a higher risk of recurrence would need to be considered for post-transplant surveillance.
We found that high expression of ALDH1A1 protein and transcript level in the primary HCC tumor was associated with lower rates of recurrence post-LT. The expression of ALDH1A1 has been reported as a favorable prognostic factor in HCC, being significantly associated in previous literature with low serum levels of alphafetoprotein, well-differentiated pathology and a favorable clinical outcome [36][37][38]. Similarly, ALDH1A1 had significantly downregulated phosphorylation in HCC compared to adjacent tissue [11]. Consistent with our results, in a gene ontology annotation study of the ALDH1 family [39], HBV-related HCC patients who showed high ALDH1L1 gene expression had superior clinical outcome, with a 57-month recurrence-free survival. High ALDH1B1 expression was protective for HCCs with multiple nodules and high serum alpha-fetoprotein level [39]. Guo et al. [40] have previously reported that long noncoding RNA MACC1-AS1 promoted the stemness of HCC cells by antagonizing ALDH1 [39]. Our study in the LT context is thus in keeping with previous literature on ALDH1A1 in non-transplant HCC.
ALDH1A1 catalyzes the oxidation of the retinol metabolite retinaldehyde to retinoic acid [33]. Thus, lower ALDH1A1 expression in our recurrent cases could represent reduced production of retinoic acid. In concordance, experimental evidence shows that retinoid activity is reduced in HCC cell lines. Moreover, patients with HCC display decreased hepatic retinoid stores and altered retinoid signaling [41]. Intriguingly, retinoic acid treatment inhibited LGALS3 expression in carcinoma cells [42], suggesting that reduced levels of retinoic acid could promote LGALS3 upregulation. In turn, increased levels of LGALS3 have been associated with a higher ability of the tumors to neutralize the immune system [32]. ALDH1A1 and LGALS3 were differentially expressed in opposite directions between recurrent and non-recurrent samples, and their expression strongly and negatively correlated with each other. In addition, we demonstrated a significant increase in the LGALS3/ALDH1A1 ratio in recurrent samples, as compared to non-recurrent samples. This finding may suggest that ALDH1A1 and LGALS3 are functionally related and contribute to HCC invasiveness, particularly in tumours with high LGALS3 relative to ALDH1A1 expression. LGALS3 expression levels were directly associated with a higher risk of HCC recurrence post-transplant in our cohort.
LGALS3 was reported to be critical for Ras signaling and thus supports mitogen activated protein kinase (MAPK) and AKT cascades. Increased levels of LGALS3 have been associated with a higher ability of the tumors to neutralize the immune system [43]. Along the same lines, LGALS3 protein expression was increased in HCC tumours, compared to adjacent tissue [11]. The same study demonstrated that both LGALS3 and LGALS3BP were significantly increased in more aggressive HCC tumours. Our survival investigation on an external HCC dataset from TCGA showed that, despite these samples representing a clinical context different from our own, the results were supportive, by revealing worse prognosis of tumors expressing high LGALS3 and LGALS3BP and low ALDH1A1. LGALS3BP interacts with LGALS3 and was found to be upregulated and highly correlated with LGALS3 in our study. Prior studies indicated that LGALS3BP can be used along with alpha-fetoprotein to improve the screening sensitivity for HCC [43].
LGALS3 and LGALS3BP are thus candidate proteins for identifying HCC tumours with increased likelihood for recurrence post-LT, which might represent potential therapeutic targets to diminish HCC recurrence post-LT [45].
Other studies have examined transcriptional markers of HCC recurrence post-LT. Miltiadous et al. [46] conducted a pivotal study examining gene expression profiling of HCC beyond the Milan Criteria and identified the progenitor phenotype (determined by either CK-19 expression or the Hoshida S2 signature) as predictive of post-LT recurrence. The S2 HCC subclass was linked to high levels of alpha-fetoprotein expression, activation of mTOR and IGF signaling, even at the early stages of the disease. This study certainly pinpointed potential gene signatures predictive of recurrence. Importantly, our study has extended this approach by generating a combined transcriptomic and proteomic signature predictive of HCC recurrence post-LT, which permits cost-effective immunohistochemical evaluation of a proteomic signature.
Our study had limitations. First, our sample size was limited, and it was not possible to adjust for significant clinical and pathological covariates such as tumor burden on explant, presence of microvascular invasion and AFP at time of transplant. Nonetheless, we included only samples from Hepatitis B-induced HCC, the most common reason for HCC liver transplant beyond Milan criteria, in order to minimize heterogeneity. This is a pilot study and additional studies with larger sample sizes and etiologies beyond hepatitis B will be needed to validate our results. Additionally, with long pre-LT waiting times LGALS3, and GAPDH in recurrent (n = 7) and non-recurrent (n = 4) HCC patients. Cropped images of the blots were used in order to improve the clarity and conciseness of the presentation. The protein expression of ALDH1A1 (B) and LGALS3 (C) were measured by densitometry and normalized to GAPDH. The Pearson correlations between Western Blot and mass spectrometry log2-transformed LFQ protein intensity values of ALDH1A1 and LGALS3 were evaluated (D). The Western blot LGALS3/ALDH1A1 ratio was also calculated (E). *P < 0.05. Data are reported as mean ± standard error. LFQ, label-free quantification; LGALS3, galectin-3; ALDH1A1, retinal dehydrogenase 1.
LGALS3, galectin-3; LGALS3BP, galectin-3-binding protein; ALDH1A1, retinal dehydrogenase 1 before transplantation and bridging therapy, the stability of tumor gene/protein expression profiles over time must be considered, should predictive signatures be incorporated into pre-LT patient selection. Despite these shortcomings, our observations provide an initial comprehensive and combined proteomic and transcriptomic profile of HCC explants for patients identified pre-LT beyond Milan criteria. Integration of high-throughput omics profiles is particularly powerful to decipher the key molecules involved in carcinogenesis [47][48][49][50]. Because expression profiling could be performed on liver needle biopsies in a pretransplant setting [51], the positivity of the proposed signatures might also assist in the decision making for selecting patients for transplantation.
In conclusion, we have derived an HCC explant protein signature comprised of ALDH1A1, LGALS3 and LGALS3BP to predict risk of HCC recurrence posttransplant. This short list of proteins was identified using both high-throughput proteomics and transcriptomics, and validated using immunohistochemistry and immunoblotting on HCC explant specimens. Staining of explant specimens for these 3 proteins could provide further guidance regarding screening protocols for patients transplanted for HCC beyond Milan. These findings support the possibility that molecular analysis of HCC performed prior to LT may enable expansion of the current imaging-based selection criteria to include patients beyond Milan Criteria with favorable A B Fig. 6 Validation of LGALS3B protein expression changes in HCC tumour explants. Increased protein expression of LGALS3BP in HCC tumour explant samples from patients that developed recurrence post-LT (n = 7), compared to non-recurrent cases (n = 4), were orthogonally verified by immunohistochemistry in the discovery cohort (A). The same trend was observed in an independent, internal validation cohort of 16 recurrent and 13 non-recurrent HCC cases. Magnification: 20x. Scale bar: 200 μm. HCC, hepatocellular carcinoma; LGALS3BP, galectin-3-binding protein