Patients and tissue
This study has been approved by the Ethics Committee of Xinhua Hospital (Shanghai, China), and the study was informed in accordance with Declaration of Helsinki. Written informed consent was obtained from the patients before the study began.
Tumor tissues and adjacent normal tissues were obtained from 26 patients underwent radical resection between 2015 and 2017 at Xinhua hospital. Each tissue sample was snap-frozen in liquid nitrogen for further analysis. All of the patients in this study belonged to the same ethnic group. The patients were selected according to the criteria: (1) All clinicopathological diagnoses were confirmed by two pathologists. (2) None of the patients received any treatments before surgery. (4) Availability of complete follow-up data and not lost follow-up. (5) No death in the perioperative period. (6) No history of other synchronous malignancies.
RNA sequencing
Total RNA from tissues and cells were isolated using Hipure Total RNA Mini Kit (Magen) according to the protocol. RNAs were eluted with 50μl of RNase-free water and then re-purified the concentration and integrity of the extracted total RNA was estimated by Qubit 3.0 Fluorometer (Invitrogen, Carlsbad, California), and Agilent 2100 Bioanalyzer (Applied Biosystems, Carlsbad, CA), respectively. RNA samples with a RIN value of at least 7.0 or higher was used for further processing. RNA-seq library was prepared with approximately 1μg of total RNA using KAPA Stranded RNA-Seq Kit with RiboErase (HMR) for Illumina Platforms (Kapa Biosystems, Inc., Woburn, MA). Briefly, ribosomal RNA was removed from total RNA. Next, first strand and directional second strand synthesis were performed. Then the A tailing and adapter ligation were performed with the purified cDNA. Finally, the purified, adapter-ligated DNA was amplified. The library quality and concentration was assessed by utilizing a DNA 1000 chip on an Agilent 2100 Bioanalyzer. Accurate quantification for sequencing applications was determined using the qPCR-based KAPA Biosystems Library Quantification kit (Kapa Biosystems, Inc., Woburn, MA). Each library was diluted to a final concentration of 10 nM and pooled equimolar prior to clustering. Paired-End (PE) sequencing was performed on all samples. For circRNA expression analysis, the read was to mapped genome using the STAR and DCC was used to identify the circRNAs and to estimate the circRNAs expression. TMM (trimmed mean of M-values) was used to normalize the gene expression. Differentially expressed genes were identified using the edgeR program.
Cell cultures
Human cholangiocarcinoma cell lines (RBE, HCCC-9810 and QBC-939), and human normal bile duct epithelium cell line (H69) were obtained from the Cell Bank of the Shanghai Cell Bank of the Chinese Academy of Sciences (Shanghai, China). Cells were maintained at 37℃ in a 5% CO2 humidified incubator, and cultured in RPMI-1640 (Gibco) (RBE and QBC-939) or DMEM (Gibco) (HCCC-9810) supplemented with 10% fetal bovine serum (Gibco) and 1% antibiotics(Gibco). Cells have not been in culture for longer than 2 months.
RNAi and plasmid construction
SiRNAs duplexes were synthesized by GenePharma (Suzhou, China) and transfected into cells using lipofectamine 2000 (Invitrogen) according to the manufacturer’s protocol. Lentivirus for knockdown of circMBOAT2 and plasmid for circMBOAT2 overexpression were obtained from GeneChem (shanghai, China). The target sequences for constructing lentiviral shRNAs and siRNAs are listed in Additional file 2.
RNA extraction and qRT-PCR analysis
Total RNA derived from ICC tissues and cells was extracted using TRIzol reagent (TaKaRa, Dalian, China) according to the manufacturer’s instructions. RNA was reverse transcribed into cDNA using a Primer-Script one step RT-PCR kit (TaKaRa, Dalian, China). Hieff UNICON® qPCR SYBR® Green Master Mix (Yeasen, Shanghai, China) was used for qRT-PCR. The circRNA and mRNA levels were normalized by GAPDH, β-actin or U3. The fold change in relative expression of RNAs was calculated using the 2−ΔΔCt methods. Oligonucleotides sequences are listed in Additional file 2.
RNase R treatment
Two micrograms of total RNA was incubated 30 min at 37 °C in the absence or presence of 5 U/μg RNase R (Geneseed, Guangzhou, China) and the resulting RNA was subsequently purified by RNeasy Mini Kit (Qiagen, Germany), and then analyzed by qRT-PCR.
RNA fluorescence in situ hybridization (FISH)
Oligonucleotide-modified probe sequence for circMBOAT2 and FASN was synthesized from Sangon Biotech (Shanghai, China). Fixed cells were washed in PBS. The cell suspension was pipetted onto autoclaved glass slides, followed by dehydration with 70, 80 and 100% ethanol. Then hybridization was performed at 37 °C overnight in a dark moist chamber. After being washed twice in 50% formamide/2 × SSC for 5 min, the slices were incubated with the regents in Alexa FluorTM 488 Tyramide SuperBoost™ Kits (Thermo Fisher Scientific, Waltham, USA) for 30 min and sealed with parafilm containing DAPI. The images were acquired using a fluorescence microscopy (OLYMPUS, Japan). The probe sequences were shown as below:
circMBOAT2: 5’- Cy3- CACTACAAAGTTGACTTGTGCATGTTCTCCACT- 3’
FASN: 5’- digoxin- GCGTAGGATGGAATCTCGGAAGCGGTC- 3’
In vitro cell phenotypic assays
For CCK-8 proliferation assay, 2×103 cells were seeded in 100 μl complete culture media in 96-well plates for various time periods. Cell Counting Kit-8 assay was performed to measure cell viability according to the manufacturer’s instructions. For 5-Ethynyl-2’-deoxyuridine (EdU) proliferation assay, BeyoClickTM EdU proliferation assay (Beyotime, Shanghai, China) was used according to the manufacturer's protocol. For migration assays, transwell filter champers (8-μm pore size, Corning, NY) were used according to the manufacturer’s instructions. Cells migrated through the membrane were fixed by 4% paraformaldehyde, stained with 0.1% crystal violet, and counted under a light microscope. Cells were incubated with 10μM EdU for two hours, stained with DAB then visualized under a light microscope. For scratch experiment, cells were plated into a six-well plate after 48-hours transfection then scratched by pipette tip. After 24-hours incubation with serum-free medium, the wound width was checked.
Flow cytometric analysis
For cell cycle analysis, after 48 h of incubation, transfected ICC cells were washed with cold phosphate-buffered saline (PBS) and incubated in ice-cold 70% ethanol at 4 °C overnight. Then, cells were incubated with propidium iodide for 30 min and analyzed for cell cycle distribution using a flow cytometer (FACS Calibur, BD Biosciences, USA). Data were analyzed by FlowJo 10.6.2 software and presented as the percentage of cell phase distribution including G0/G1, S and G2/M phases. For cell apoptosis analysis, the cultured cells were stained using annexin V-fluorescein isothiocyanate and propidium iodide (BD Biosciences, USA).
Immunofluorescence assays
Tissues were fixed using 4% paraformaldehyde and embedded in paraffin. The samples were then incubated with primary antibody with PTBP1 at 4 °C overnight and then incubated with goat anti-rabbit IgG with a red or green fluorescent label (Invitrogen, Carlsbad, CA).
Western blot
Protein was extracted from transfected cells with RIPA lysis buffer. Equal amounts of protein samples were loaded and separated by SDS-PAGE and transferred onto PVDF membranes (Merck Millipore, Germany). The membranes were blocked with 5% skim milk in TBST for 1 hour at room temperature. Then, membranes were incubated with diluted primary antibodies anti-PTBP1 for western blot (1:1000 dilutions, 12582-1-AP, Proteintech Group, USA), anti-FASN for western blot (1:1000 dilutions, 10624-2-AP, Proteintech Group, USA) at 4℃ overnight. Then, membranes were washed with TBST three times, 15 minutes/time, followed by incubating with secondary antibody for an hour, washed again with TBST. Finally, the protein bands were visualized by using Gel Doc 2000 (Bio-Rad) and the grey values were measured by Image J software.
Transcription in vitro and RNA pull-down assay
For transcription in vitro, plasmid containing double T7 promoters was digested by single restriction endonuclease. T7 High Yield RNA Transcription Kit (Vazyme, Nanjing, China) was used to transcript forward and reverse linear DNA template to RNA (sense and antisense probe), which was subsequently purified by RNeasy Mini Kit (Qiagen, Germany). Then, Pierce™ RNA 3' End Desthiobiotinylation Kit (Thermo Fisher Scientific, Rockford, USA) was used to label biotin to the 3' End of RNA. For RNA pull-down assay, 1×107 cells were washed in ice-cold PBS, lysed in 500 μl co-IP buffer (Thermo Scientific) supplemented with a cocktail of proteinase inhibitors, phosphatase inhibitors, and RNase inhibitor (Invitrogen), then incubated with 3 μg biotinylated DNA oligo probes against sense or antisense for 2 h at room temperature. A total of 50 μl washed Streptavidin magnetic beads (Thermo Fisher Scientific, Rockford, USA) were added to each binding reaction and further incubated for another hour at room temperature. The beads were washed briefly with elution buffer for five times. Finally, the retrieved proteins were used for mass spectrometry or western blot analysis.
Silver staining and mass spectrometry analysis
Silver staining was performed using the Fast Silver Stain Kit (Beyotime, Shanghai, China) as the protocol described, while MS was done by BGI Genomics (Shenzhen, China). The protein identification uses experimental MS/MS data and aligns them with theoretical MS/MS data from database to obtain results. The whole process starts from converting raw MS data into a peak list and then searching matches in the database. The search results are subject to strict filtering and quality control, and possible protein identifications are produced. Finally, from the final protein identification list, functional annotation analysis such as GO, COG/KOG, and Pathway analysis are performed.
RNA immunoprecipitation
RIP experiments were performed with a RNA Immunoprecipitation Kit (Geneseed, Guangzhou, China) according to the manufacturer’s instructions. Co-precipitated RNA was detected by qRT-PCR.
Lipidomics
Lipid extraction and mass spectrometry-based lipid detection of cell pellet were followed by experimental operation of Huang et al25. A quality control sample was prepared by mixing equal parts of all samples. During the whole experiments, we processed with the same parameters as the analytical sample to evaluate the stability of analytical performance and reliability of data. Ultra-high-performance liquid chromatograph-mass spectrometry analysis was performed on Q Exactive Plus high-resolution mass spectrometer (Thermo Scientific, USA) equipped with an Ultimate 3000 UHPLC system (Thermo Scientific, USA). Lipid identification (structural identification) and a peak table containing the retention time, m/z, peak area (peak matching) were assessed with the MS-DIAL software. The data combination of positive and negative ion modes was defined as relative lipid abundance for subsequent statistical analysis.
Co-Immunoprecipitation
To detect protein–protein interactions, cells were lysed in 500 μl co-IP buffer supplemented with a cocktail of proteinase inhibitors, phosphatase inhibitors, and RNase inhibitor. The lysates were centrifuged at 12,000g for 30 min, and the supernatant was used for immunoprecipitation with agarose beads, which were preincubated with the corresponding antibodies. After incubation at 4 °C overnight, beads were washed three times with PBS. SDS sample buffer was added to the agarose beads and the immunoprecipitates were used for western blot analysis.
Statistical analysis
Results were shown as the mean ± SD. SPSS 22.0 (IBM Corp., Armonk, NY, USA), R 4.0.5 and GraphPad Prism 8.0 were used for comparison analysis. Student t-test was used for comparison between two groups. The Chi-squared test was used for the association of the expression of circMBOAT2 with patients’ clinic pathological parameters. P < 0.05 were considered statistically significant.