Microarray analysis of CircRNAs
The circRNA microarray analysis (CapitalBio Technology, Beijing, China) consisted of specifical probes targeting back-splicing sites of human circRNAs was performed. The circRNAs in 5 pairs of CRC tissues were detected. The data normalization, summarization, and quality control were analyzed by the GeneSpring software. The differential circRNAs are listed in Table S1. All the raw data of the microarray analysis can be found in the Gene Expression Omnibus with the number GSE172229.
Patient samples and cell culture
One hundred CRC patients' tumor samples and paired adjacent normal tissues were collected between 2014 and 2018 from the department of general surgery, First Affiliated Hospital of Nanjing Medical University. Fresh tissues were acquired after surgery and immediately frozen at -80°C. Patients who had received neoadjuvant chemoradiotherapy were excluded. This study was approved by the Human Ethics Committee of First Affiliated Hospital of Nanjing Medical University. Each patient was informed about the study and consented to the specimen donation before the surgery. The raw data of circ-GALNT16 expression in 100 pairs of CRC tissues and relative expression in different pathological subgroups are listed in Table S2.
LoVo, DLD-1, SW-480, Caco-2, RKO, HCT 116, HT-29, NCM460 were purchased from Cell Bank of Type Culture Collection of the Chinese Academy of Science (Shanghai, China) and cultured under recommended medium with 10% fetal bovine serum at 37°C in a 5% CO2 humidified incubator.
RNA extraction and qRT-PCR
Total RNA was extracted from patients’ specimens and cells by TRIzol reagent (Invitrogen, USA), and the procedures were applied as previously described (20). Then, a PrimeScript RT reagent kit (TaKara, Dalian, China) was used for reverse transcription. An SYBR Premix Ex Taq Kit (TaKara) was used for the qRT-PCR experiment on Applied Biosystems 7500 sequence Detection System. GAPDH was used as an internal control. The primers’ sequences are listed in Table S3.
Cell Transfection
The lentivirus containing shRNAs and full-length targeting circ-GALNT16, including corresponding negative control, were synthesized by Obio (Shanghai, China). Serpine1 overexpression plasmid was obtained from Obio. Obio synthesized the truncated hnRNPK plasmids of hnRNPK with a C-terminus 3× Flag tag. The small interfering RNA (siRNA) oligonucleotides targeting hnRNPK, SENP2, and Serpine1 were synthesized by RiboBio (Guangzhou, China). Procedures were described previously, using Lipofectamine 3000 (Invitrogen) (21). The shRNAs and siRNAs sequences are listed in Table S3.
Cell proliferation assays
Cell proliferation assays, including Cell Counting Kit-8 and colony formation assay, were performed as described previously (22). The cell proliferation ability was measured via the 5-Ethynyl-2′-deoxyuridine assay using an EdU kit (Beyotime, Shanghai, China).
Transwell and scratch wound healing assay
Transwell assay was performed with Millipore cell culture inserts (24-well insert, 8-µm pore size). 3×104 cells resuspended in 200µL serum-free medium per well were seeded into the upper chamber of Transwell membrane while the lower layer was filled with 700µL culture medium containing 10% serum for migration assay. After 24h (36 or 48h based on different cell lines), we stained the cells on the underside of the membrane using crystal violet solution and wiped out the cells above the membrane with swabs. Five random views were counted in each well under a microscope. For invasion assay, Matrigel (BD Biosciences, Franklin Lakes, NJ, USA) was spread on the upper layer, and the rest of procedures were followed the steps described above.
A 200µL pipette tip was used to produce consistent length lesions in six-well plates (8×105 cells per well ) for wound healing assay. An inversion microscope was used to take each wound image at 0 and 48 h. We used ImageJ to quantitatively evaluate the gap distance.
Flow cytometry assay of cell cycle and apoptosis
Treated cells were fixed in 75% alcohol overnight at -20°C for cell cycle assay. The cells were washed three times, and a Cell Cycle Analysis Kit (Beyotime, Shanghai, China) was used for propidium iodide (PI) staining.
For apoptosis analysis, all cells were treated with 0.5 mM of H2O2 for 4 h to stimulate apoptosis. An Annexin VAPC/7-AAD Apoptosis Detection Kit (KeyGEN, Jiangsu, China) was used for Annexin V-APC and 7-AAD staining according to the protocol. Finally, the percentage of cell cycle distribution and apoptotic rate were analyzed using BD FACSCanto II (BD Biosciences, San Jose, CA, USA).
Animal models
Four weeks old male BALB/c nude mice were purchased from the Animal Center of Nanjing Medical University (Nanjing, China) for subcutaneous tumor formation and liver metastasis model. For the xenograft model, 106 DLD-1 cells transfected with sh-circ-GALNT16#1 and RKO cells transfected with circ-GALNT16, as well as corresponding control cells, were resuspended in 1mL PBS (phosphate-buffered saline), and respectively subcutaneous injected into the left and right armpits of the mice. Each tumor volume was calculated every five days since the tumor became macroscopic. All the mice were sacrificed 30 days after the subcutaneous injection, and the xenograft tumors were weighed and obtained for IHC.
106 cells suspended in 20µL PBS were injected into the distal tip of the mice's spleen to establish the metastatic model. The liver tissues were dissected and embedded in paraffin for H&E staining 6 weeks later. All these animal experiments were ratified by the Committee on the Ethics of Animal Experiments of Nanjing Medical University.
Immunohistochemistry
Immunohistochemistry (IHC) was performed as previously described (22). All the antibodies used are listed in Table S4.
RNA and protein isolation of nuclear and cytoplasmic fractions
Nuclear and cytoplasmic fractions of CRC cells were separately isolated with a PARIS™ kit (AM1556, Thermo Fisher Scientific, Waltham, USA). RNA was isolated according to the protocol and analyzed through qRT-PCR. U6, 18S, and GAPDH were used as internal controls. The isolated nuclear protein was prepared for subsequent pulldown assay.
Fluorescence in situ hybridization and Immunofluorescence
Circ-GALNT16 specifical Cy3-labelled probe was applied to detect the subcellular localization of circ-GALNT16 in DLD-1 and LoVo using a FISH Kit (RiboBio). Briefly, after cells were fixed with paraformaldehyde for 10min and permeabilized for 5min using PBS with 0.5% Triton X-100, FISH probe hybridized with preheated hybridization buffer was mixed with cells at 37°C overnight. Then, the cells were washed with 4×sodium citrate buffer containing 0.1% Tween-20 for 5 min and 1× SSC for 5 min. The cell nucleus was stained with 4,6-diamidino-2-phenylindole (DAPI). The images were obtained using a confocal fluorescence microscope.
For dual RNA-FISH and immunofluorescence assay, an immunostaining blocking solution (Beyotime) was used for cell blocking for 1 h after incubated with the FISH probe as described above. Then cells were incubated with hnRNPK antibody overnight and labeled with fluorescent secondary antibody for 1h under dark condition. Finally, DAPI was used for nuclear staining.
RNA pulldown assay and mass spectrometry
Biotin labeling pulldown probe targeting circ-GALNT16 and control probe were designed and synthesized by Ribobio. RNA pulldown assay was constructed using Pierce Magnetic RNA-Protein Pull-Down Kit according to the protocol (#20164, Thermo). Mass spectrometry analysis was performed in the elution protein extracted from RNA pulldown assay. The differential proteins identified by mass spectrometry and the RNA binding proteins list (RBPs; http://www.ablife.cc) are shown in Table S5.
Western blot
Western blot (WB) was carried out as previously reported (20). The antibodies are shown in Table S4.
RIP and Coimmunoprecipitation assay
RNA immunoprecipitation (RIP) assay was performed using a RIP Kit (Millipore, Burlington, MA, USA). In Brief, 5µg anti-hnRNPK or anti-FLAG antibody and magnetic beads were mixed and incubated with cells lysis lysed by RIP lysis buffer supplemented with protease and RNase inhibitors overnight at 4°C. The immunoprecipitated RNA was obtained for qRT-PCR after being digested with proteinase K buffer.
Coimmunoprecipitation assay was constructed with an IP/Co-IP Kit (#88828, Thermo) to explore the interactions between hnRNPK with p53 and SENP2. The detailed procedures were performed as previously described (23).
SUMOylation modification analysis
The treated cell lysate was diluted 20-fold with lysis buffer consisted of 20mM Tris-HCl (PH = 8.0), 150mM NaCl, 2 mM NEM, 1 × protease inhibitor cocktail, and 0.2% Triton X-100. The lysate was incubated with anti-hnRNPK for 4 hours at 4°C. Then, the lysate was mixed and incubated with protein A/G-Sepharose (sc-2003; Santa Cruz, USA) overnight at 4°C. Resins were washed with lysis buffer with 1% Triton X-100 for three times and boiled in 50µL SDS sample buffer for 10 min. Finally, WB was performed with supernatants using anti-hnRNPK and anti-SUMO1.
RNA sequencing assay
RNA-seq libraries were prepared after RNA-seq was constructed in RKO cells transfected with circ-GALNT16 overexpression lentivirus and control cells. Firstly, RNA samples that had passed quality inspection were constructed with a starting amount of 1µg and performed RNA Integrity checking by Agilent 4200 TapeStation. The key steps were as follows: RNA purification by polyA Oligo magnetic beads; first-strand cDNA synthesis with random hexamer primers; RNA degradation with RNase H and second-strand cDNA synthesis using DNA polymerase I; double-stranded cDNA fragments end repair and the addition of a single 'A' base at the 3'-end of each strand; special sequencing adapters ligation; PCR amplification. Finally, a HiSeq 2000 system on Pair End (Illumina, San Diego, CA, USA) was used to sequence the purified cDNA. The mRNA sequencing results are shown in Table S6.
Chromatin immunoprecipitation assays
A ChIP Kit (CST, #56383, Danvers, MA, USA) was used for ChIP assays following the protocol. Briefly, chromatin fragments were sheared to be concentrated on 200 to 1000 bp through ChIP sonication lysis buffers. Then beads-antibody complexes were incubated with sheared crosslinked chromatin. After the DNA was purified, qRT-PCR was performed using ChIP primers. The sequences of ChIP primer are listed in Table S3.
Statistical analysis
Each experiment was repeated at least three times. All the statistical analyses were carried out with GraphPad Prism software (La Jolla, CA, USA) and SPSS 13.0 software (Chicago, IL, USA). Student's t-test was performed to analyze the difference between two samples, while ANOVA was used for tests among more than two groups. Pearson correlation analysis was performed to estimate the correlation between circ-GALNT16 and Serpine1. Chi-square test was used to analyze the correlation between circ-GALNT16 expression in tumor tissues and relative clinicopathologic data. Overall survival (OS) rates were estimated using Kaplan-Meier analyses. The significance threshold of each test was set at 0.05.