Cell culture and conditioned medium (CM) collection
HCC cell lines MHCC97H and MHCC97H with integrin-β1 knockdown (Liver Cancer Institute of Fudan University, Shanghai, China), Hep3B and HepG2 (ATCC, USA), Huh7 (Japanese Cancer Research Bank) were grown in DMEM with 10% FBS (Gibco) and 1% penicillin/streptomycin. Cell lines were authenticated by short tandem repeat validation analysis during the study period. Primary human hepatic stellate cells (pHSCs) (Sciencell, USA) were maintained in the provided medium and LX2 cells (a gift from S. Friedman) were cultured in DMEM with 2% FBS. All cell cultures were carried out in a 37 °C incubator with a humidified atmosphere in presence of 5% CO2.
As in our previous description [26], conditioned medium was collected from activated HSCs (HSC-CM) and anti-human POSTN antibody (2.5 μg/mL) (Abcam, Cambridge, UK) was added into HSC-CM to neutralize the activity of POSTN.
To obtain conditioned medium from calcipotriol-treated HSCs, pHSCs or LX2 cells were pre-stimulated using 10 ng/mL TGF-β1 and then incubated with 100 nM calcipotriol (Sigma-Aldrich) for 12 h, replenished with fresh medium for another 24 h and the medium was collected for the subsequent experiments.
In vitro heat treatment
As previously described [26], HCC cells were heated at the pre-set temperatures, seeded into 96-well plates and cell viability was measured at 48 h after heat treatment. IT50 was calculated to indicate the temperature of inducing a 50% reduction in cell viability compared with the 37 °C control. The IT50 data was used to simulate in vivo sublethal heat condition.
Cell proliferation
Cell proliferation was measured using the cell viability assay WST-1 or BrdU ELISA kit (Roche, Manheim, Germany) according to the manufacturer’s instructions. Briefly, HCC cells were subjected to sublethal heat treatment and plated in 96-well plates (8 × 103 cells/well). Equivalent amount of HSC-CM or control medium was added into each well and replaced daily to remove cell debris. The absorbance was measured at 450 nm on a multiskan spectrum reader (Thermo Scientific).
Cell co-culture were performed using a transwell system containing the 0.4 μm pore filter insert (Millipore, Switzerland) allowing medium components freely diffusion but restricting cells migration. Briefly, heat-treated HCC cells were plated in the bottom alone or together with HSCs cells in the transwell insert (4:1 ratio). After 4 days, the proliferation of HCC cells in the bottom wells were measured using WST-1 reagent.
Cell migration and invasion assay
As previously described [27], migration and invasion assays were carried out. Briefly, in the migration assay, heat-treated residual HCC cells (1 × 105) were seeded into the upper chamber while 600 μL HSC-CM or control medium was added into the lower chamber, and then incubated for 36 h. The invasion assay was done in a similar manner except transwell insert pre-coated with Matrigel (BD Biosciences) and culture for 72 h. At the end of experiments, the cells in the upper chamber were removed with cotton swabs and the cells on the lower surface of the inserts were fixed, stained and photographed from five random fields (100× magnification) under phase contrast microscope (Leica, Germany).
Colony formation assay
After sublethal heat treatment, HCC cells (2 × 103 cells) were seeded into 6-well plates and equal volume of HSC-CM or control medium was then added into each well for 8 days. Cell colony with a diameter larger than 50 μm was counted using ImageJ software (National Institutes of Health, Bethesda, MD, USA) after fixed with methanol, stained with crystal violet and photographed.
Flow cytometry analysis
Cell apoptosis was detected using annexin V staining kit (Invitrogen, Carlsbad, California, USA) according to the manufacturer’s instruction. In brief, heat-treated HCC cells were incubated with HSC-CM or control medium for 3 days. Then, cells were harvested and resuspended in annexin-binding buffer (1 × 106 cells/mL). Subsequently, appropriate amount of Alexa Fluor 488 annexin V and PI working solution were added. Early and late apoptosis rates were analyzed by FACS Calibur flow cytometer (BD Biosciences, USA) and FlowJo software (Tree Star Inc, Ashland, Ore).
Quantitative reverse transcription-PCR (qRT-PCR)
Total RNA was extracted using TRIzol reagent, transcribed and amplified with the use of RevertAid First Strand cDNA Synthesis kit and the Maxima SYBR Green qPCR Master Mix kit (Thermo Fisher Scientific, Waltham, MA, USA). The relative gene expression was calculated using the equation 2−ΔΔCt where the Ct value of GAPDH or β-actin was used as the normalization. The PCR primers are shown in Additional file 1: Table S1.
Western blot analysis
Western blot was performed as the previously described [26]. Briefly, total protein was extracted using RIPA buffer added with PMSF and phosphatase inhibitors. Protein concentration was determined using BCA protein assay (Millipore, Switzerland). Then, 20 μg of protein was separated by SDS-PAGE and transferred to polyvinylidene difluoride membranes (PVDF, Millipore, USA). After blocked with 5% non-fat milk, the membrane was incubated with primary antibody against PCNA (1:2000), Snail (1:1000), Vimentin (1:1000), E-cadherin (1:1000), N-cadherin (1:1000), integrin β1 (1:1000), Shc (1:1000), phosopho-Shc (1:2000), ERK1/2 (1:1000) or phosopho-ERK1/2 (Thr202/Tyr204) (1:2000) (Cell Signaling Technology, USA), periostin (1:1000), collagen I (1:5000), α-SMA (1:300) or anti-vitamin D receptor antibody (1:1000) (Abcam, USA), GAPDH (1:1000), tubulin (1:1000) or β-actin (1:1000) (Beyotime, China) at 4 °C overnight and the corresponding HRP-conjugated secondary antibody for 1 h at room temperature on the next day. The membrane was developed with enhanced chemiluminescence (ECL; New Cell & Molecular Biotech Co., China).
Immunohistochemical analysis
As previously described [28], immunohistochemistry was conducted using the streptavidin–peroxidase two-step method. Briefly, tissue section was deparaffinised, rehydrated, heated with antigen retrieval, blocked with 3% hydrogen peroxide and then incubated with primary antibody against POSTN (1:200), α-SMA (1:100), PCNA (1:400), cleaved-Caspase-3 (1:100) or E-cadherin (1:100) (Cell Signaling Technology, USA) at 4 °C overnight. The immunoreactivity and nuclear counterstain was performed using EnVision two-step visualization system (GeneTech, Shanghai, China) and hematoxylin. Photographs of three representative fields (200× magnification) were captured by microscope (Leica, Germany).
Dual luciferase assay
The 1284-base pair cDNA containing the coding region of vitamin D receptor (VDR), named GV141-VDR, was purchased from Genechem (Shanghai, China). A luciferase reporter construct plasmid GV238 that contained the 5ʹ-flanking region of the POSTN gene from − 2000 to − 1 relative to the transcription start site (corresponds to − 2000 relative to the first AGA of 5′ UTR, amplified with the primers 5′-GATAGGTACCGCAAAGAACGACTAGGTTAAAATTG-3′ and 5′-AGATCTCGAGGAACTCTTTCCAGGAAGCATCGG-3′), as well as a series of deletion constructs (deleted at intervals of 500 base pairs from the 5ʹ end), were generated by GeneChem (Shanghai, China). Based on the conserved motif of VDR, the JASPAR database (http://jaspar.genereg.net/) was used to predict the binding sites of VDR on POSTN promoter sequence. Two putative sites in the promoter sequence of POSTN were identified (− 1829 to − 1815 and − 1415 to − 1401). To determine functionally relevant VDR binding sites, the putative binding sites for the seed region of VDR were mutated using site-specific sequence mutagenesis (CCCCACGGTTTCCAT). All promoter constructs were sequenced. HEK 293T cells were seeded into 24-well plates at 50–60% confluence per well. Cells were transfected with 1 μg POSTN promoter plasmids containing the firefly luciferase reporter, 1 μg VDR expression vector GV141-VDR and 20 ng Renilla luciferase expressing construct (as an internal control) into HEK 293T cells using Lipofectamine 2000 reagent (Invitrogen, CA, USA) and according to the manufacturer’s protocol, and cells were dealt with 100 nM calcipotriol for 12 h after 12 h of transfection. Twenty-four hours after transfection, the cells were detected for luciferase activity using the Dual-Luciferase Reporter Assay System (Promega, Madison, WI, USA) and multi-plate reader. Firefly luciferase activity is normalized to the Renilla luciferase activity. Each experiment was independently performed at least two times in triplicate.
Gene expression profile
Total RNA was extracted and purified from heat-exposed residual MHCC97H cells treated with or without 100 ng/mL POSTN using miRNeasy Mini kit (QIAGEN, Germany). Genome-wide expression profiling was conducted using Affymetrix GeneChip Human Genome U133 Plus 2.0 Array according to the manufacturer’s instructions. Briefly, total RNA samples were amplified, labelled and purified to obtain biotin-labeled cRNA using GeneChip 3′IVT PLUS reagent Kit (Affymetrix, CA, US). The arrays were hybridized, washed and stained using the GeneChip Hybridization Wash and Stain Kit, and then scanned with a GeneChip Scanner 3000 (Affymetrix, CA, US) for subsequent generation of raw data. The gene expression data were processed and further analyzed using MAS 5.0 algorithm, Affy packages in R. Genes significantly differentially expressed were selected based on fold change > 2 or < 0.5, and T-test with P < 0.05. Gene annotation was based on Gene Ontology analysis. KEGG pathway enrichment analysis was conducted to examine the functional association between differentially expressed genes and to generate the significant gene networks. Differentially expressed genes were mapped to STRING database (http://string-db.org/) and then imported into Cytoscape software (http://www.cytoscape.org/) to acquire protein–protein interaction (PPI) network. The microarray data have been deposited into Gene Expression Omnibus (GEO) database (GSE108853).
Gene expression profiles of HCC cohorts (374 tumoral samples and 50 normal control samples) were downloaded from Genomic Data Commons Data Portal (GDC, https://gdc.cancer.gov/). Data processing and analysis were performed by using R Statistical Software. The expression enrichment of a gene signature was assessed significance for differences between the two groups using unpaired two-tailed Student’s t-test. The strength of associations between two parameters was assessed using Spearman’s correlation test. Survival plots were generated using Kaplan–Meier analyses and splitting the tumor samples into high- and low-expression groups. The median of tumor tissues was token as the threshold and significance for differences of survival between both groups was assessed by a log-rank test.
Animal experiments
All animal experiments were conducted according to the guidelines of the Shanghai Medical Experimental Animal Care Commission and approved by the Ethical Committee on Animal Experiments of Fudan University, Shanghai, China.
Tumorigenicity assay was performed as described [29]. In brief, 2 × 104 heat-treated residual MHCC 97H cells resuspended in 50% growth-factor-reduced Matrigel (BD Biosciences) were subcutaneously injected into the flanks of NOD/SCID mice alone or with 100 ng/mL POSTN. Mice were sacrificed 2 months post injection, at which time tumors were harvested.
As we previously described [26], 4-week old BALB/c nude mice were inoculated subcutaneously with heat-treated residual MHCC 97H cells (2 × 107, n = 6) with pHSCs (ratio, 4:1) into the flanks of mice. Eighteen tumor-bearing mice were randomly assigned into three treatment groups: saline, calcipotriol (60 μg/kg, intraperitoneal, daily), or calcipotriol (60 μg/kg, intraperitoneal, daily) + cisplatin (4 mg/kg, twice a week). After 2 weeks, mice were sacrificed and tumor samples were harvested.
To examine the presence of activated HSCs in the tumors, CFSE labelled pHSCs (5 × 106) with heat-treated residual MHCC 97H cells were inoculated subcutaneously into nude mice. At 2 weeks post injection, mice were euthanized and tumors were cryosectioned to detect the HSCs using fluorescence microscope (Leica, Germany).
Statistical analysis
All statistical calculations were performed using Prism 6.0 (GraphPad Software, Inc., La Jolla, CA). Data were expressed as either mean ± standard deviation (SD) or percentage (%) when appropriate. Continuous data were analyzed by unpaired Student’s t-test or ANOVA between two groups or multiple groups whereas categorical data were compared using Chi-square test or Fisher’s exact test. Differences were considered statistically significant when P < 0.05 (two-side).