Cell culture and establishment of depleted and increased mtDNA copy number
Human ESCC cell lines, TE8 (RBRC-RCB2098) and TE11 (RBRC-2100), were purchased from the RIKEN BioResource Center (Tsukuba, Ibaraki, Japan). These cells were cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum (Thermo Fisher Scientific, Waltham, MA, USA), penicillin (100 IU/ml), and streptomycin (100 µg/ml) at 37°C in a humidified incubator with 5% CO2. Depletion of mtDNA content was induced by knockdown of mitochondrial transcription factor A (TFAM) expression (27) (28). Silencing of TFAM reduces mtDNA levels, since this factor is important for mtDNA packaging and maintenance (16) (29). A short hairpin RNA designed by Sigma–Aldrich (St. Louis, MO, USA), MISSION TRC-Hs1.0, was applied for knockdown of TFAM expression in TE8 and TE11 cells. The depleted mtDNA cell lines “tfam-sh1” and “tfam-sh2” were established by targeting the following sequences of the TFAM gene: 5’-CGTCGCACAATAAAGAACAA-3’ (TRCN0000016094) and 5’-GCAGATTTAAAGAACAGCTAA-3’ (TRCN0000016097), respectively. In addition, the nontargeted sequence 5’-GGCGCGATAGCGCTAATAATTT-3’ (SHC016, Sigma–Aldrich, defined as ‘control-sh’) was used as the control for comparison (control-sh). The morphology of the cells was evaluated by optical microscopy (BZ-X710, KEYENCE, Osaka, Japan) (16). In addition, an increased mtDNA copy number was established as follows: 1 × 106 cells were plated in 1.25 ml of complete medium without antibiotics in 6-well plates and transfected with 7 mg of a plasmid containing the TFAM-encoding gene (OriGene, RC215488) using 12 µl of Lipofectamine 3000 reagent (L3000015, Thermo Fisher Scientific, Tokyo, Japan) following the manufacturer’s instructions. The protein expression of TFAM was lower in tfam-sh1 and tfam-sh2 cells than in control-sh cells, and tfam-sh1 and tfam-sh2 ESCC cells had an approximately 40–60% reduction in mtDNA copy number compared with control-sh cells (Supplemental Fig. 1A, 1B). Additionally, the proliferation rate was significantly lower in tfam-sh1 and tfam-sh2 cells than in control-sh cells (Supplemental Fig. 1C).
Clinical samples
From November 2006 to December 2011, 88 esophageal squamous cell cancer patients underwent R0 resection after neoadjuvant chemotherapy at Osaka University Hospital (Osaka, Japan). Neoadjuvant chemotherapy followed by surgery at our hospital was performed for patients with cStage I (excluding T1N0), II, III, or IV disease without distant organ metastasis. The neoadjuvant chemotherapy regimen consisted of either ACF (adriamycin 35 mg/m2, cisplatin 70 mg/m2 on Day 1, and continuous fluorouracil 700 mg/m2 infusion for 7 days) every 4 weeks or DCF (docetaxel 70 mg/m2, cisplatin 70 mg/m2 on Day 1, and continuous fluorouracil 700 mg/m2 infusion for 5 days) every 3 weeks, as previously described. Surgery was performed after two cycles of neoadjuvant chemotherapy (30). Clinicopathological findings were classified according to the Union for International Cancer Control (UICC)-TNM classification, seventh edition (31).
Formalin-fixed, paraffin-embedded samples from the patients were collected at our hospital. Cancerous ESCC nests were subjected to DNA extraction after surgery using laser microdissection with a Leica LMD7000 instrument (Leica Microsystems, Wetzlar, Germany). Moreover, of the patients, formalin-fixed, paraffin-embedded samples were collected before and after neoadjuvant therapy from 4 patients with ESCC who had undergone esophagectomy with neoadjuvant therapy. All patients provided written informed consent for the use of the resected specimens. In addition, this study was approved by the ethics committee of Osaka University, Graduate School of Medicine (approval number #15401) and was conducted in accordance with the Declaration of Helsinki.
Measurement of mtDNA copy number and mRNA expression levels
The mtDNA copy number was measured by quantitative real-time PCR (qPCR) using specific primers for the mtDNA-coded cytochrome oxidase I (MTCO-1) gene and normalized to the expression of the nuclear DNA-encoded cytochrome oxidase IV (COX IV) gene (a mitochondrial respiratory chain enzyme), and the mtDNA copy number was adjusted by setting the mtDNA copy number of control-sh or TE11 cells as 1.00. Complex I/IV-related genes, EMT-related genes including epithelial markers (E-cadherin) and mesenchymal markers (N-cadherin, vimentin, and ZEB1), and DNA methylation-related genes such as DNMT1, DNMT3A, and DNMT3B were analyzed by qPCR. In addition, the house-keeping gene GAPDH was analyzed in this study. The primers are shown in Supplemental Table 1.
Chemosensitivity assay (viability assay)
A total of 2.0 ×103 cells per well were seeded in 96-well plates, and cell viability was evaluated using Cell Counting Kit-F (CK06, Dojindo, Japan) at 48 hours after incubation under chemotherapy, cisplatin (CDDP), 5-fluorouracil (5-FU) and docetaxel (DTX) exposure. The chemotherapy resistance rate was determined as the ratio of the proliferation under chemotherapy at each concentration to that under control treatment. The fluorescence intensity was measured by an iMark™ Microplate Absorbance Reader (BIO-RAD, Tokyo, Japan) using a plate reader at an excitation wavelength of 490 nm and emission wavelength of 515 nm.
Apoptosis analysis
For the apoptosis assay, an Annexin V-Cy3 Apoptosis Staining/Detection Kit was used (ab14142, Abcam, Cambridge, UK). A total of 2×105 cells collected by centrifugation were resuspended in binding buffer, and Annexin V-Cy3 was added. Then, the cells were incubated at room temperature for 5 min in the dark. Flow cytometry was carried out using a FACSCanto II flow cytometer (BD Biosciences), the data were analyzed by FlowJo ver 10.3 software (TOMY DIGITAL BIOLOGY, Tokyo, Japan), and mean fluorescence intensities (MFIs) were measured.
Flow cytometry analysis
Harvested cells (1.0 ×105 cells per well) were seeded in 6-well plates and stained using the MitoProbe JC-1 assay (M34152, Invitrogen) according to the manufacturer’s instructions. The mitochondrial membrane potential (MMP) was analyzed by flow cytometry. The MitoProbe JC-1 assay features a conversion of PE-A (red) to FITC-A (green) when the MMP is depolarized. Flow cytometry was carried out using a FACSCanto II flow cytometer (BD Biosciences), the data were analyzed by FlowJo ver 10.3 software (TOMY DIGITAL BIOLOGY, Tokyo, Japan), and the MFI was measured.
Fluorescent immunostaining
A total of 2.0 ×104 cells per well were seeded in multiwell glass bottom dishes (D141400, MATSUNAMI, Osaka, Japan). mtDNA in cells was stained with SYBR Green I at a 1:600,000 dilution for 30 minutes and washed with PBS four times as previously described (32). Then, the MMP was determined via staining with MitoTracker Orange (M7510, Invitrogen, California, UA) for 10 minutes followed by two washes with PBS. Fluorescence immunostaining was evaluated by a Confocal Laser Scanning Microscope FV3000 (Olympus, Tokyo, Japan). This fluorescence was quantified by ImageJ v1.53 software, which is an open source and public domain image processing software.
Complex I and IV activity assays
Complex I activity was investigated using a Complex I Enzyme Activity Microplate Assay Kit (ab14142, Abcam, Cambridge, UK). The protein in cells was added to the microplate wells that had been precoated with a specific capture antibody for complex I. Then, samples were immobilized in the well. In the assay, complex I activity is determined by detecting the oxidation of NADH to NAD + and the simultaneous decrease in dye, which leads to increased absorbance at an optimal density of 450 nm by the SH-9000lab system (Corona Electric, Ibaraki, Japan). Complex IV activity was evaluated using a Complex IV Rodent Enzyme Activity Microplate Assay Kit (ab109911, Abcam, Cambridge, UK). This assay kit was used to determine the activity of cytochrome c oxidase, which participates in complex IV activity, the activity was determined calorimetrically by detecting the oxidation of reduced cytochrome c based on the absorbance change at 550 nm.
Immunoblotting analysis
Relative protein expression levels were investigated by immunoblotting. E-cadherin, N-cadherin, vimentin, and Zeb expression levels were evaluated using a commercially available EMT antibody sample kit (1:1000 dilution, #9782, Cell Signaling Technology). DNMT1, DNMT3A and DNMT3B levels were evaluated with anti-DNMT1 rabbit monoclonal antibody (1:1000 dilution, #5032, Cell Signaling Technology), anti-DNMT3A rabbit monoclonal antibody (1:1000 dilution, ab227823, Abcam, Cambridge, UK) and anti-DNMT3B rabbit monoclonal antibody (1:1000 dilution, #67259, Cell Signaling Technology), respectively. In addition, TFAM expression levels were assessed using TFAM rabbit polyclonal antibody (1:1000 dilution, ab47517, Abcam, Cambridge, UK).
DNA methylation assay
The overall DNA methylation in cells was analyzed using a methylated DNA quantification kit (ab117128, Abcam, Cambridge, UK), which quantifies global DNA methylation by specifically measuring the levels of 5-methylcytosine (5-mC) in a microplate-based format. Genomic DNA was applied to microplates and bound to assay wells. Then, capture antibody, detection antibody, enhancer solution and developing solution were added after washing well. The fluorescence intensity was measured with an SH-9000lab instrument (Corona Electric, Ibaraki, Japan) using a plate reader at 450 nm optimal density.
Cell proliferation assay
A total of 5.0 ×102 cells per well were seeded in 96-well plates, and cell viability was evaluated using Cell Counting Kit-F (CK06, Dojindo). The fluorescence intensity was measured at an excitation wavelength of 490 nm and emission wavelength of 515 nm by an iMark™ Microplate Absorbance Reader (BIO-RAD, Tokyo, Japan).
Immunohistochemistry
Tumor specimens were fixed with 10% formalin, and paraffin-embedded tissue blocks were sectioned into 2.5-µm slices. The sections were deparaffinized with xylene, dehydrated with ethanol in stages and incubated in 10 mM citrate buffer at 110°C using a pressure cooker for 15 min for antigen removal. Endogenous peroxidase activity of the tissue specimens was blocked by incubating the slides in 3% hydrogen peroxide (H2O2) dissolved in methanol for 20 min at room temperature. Subsequently, the sample was treated with 1% horse serum albumin for 30 min at room temperature to block nonspecific reactions, all sections were incubated with primary antibodies at 4°C overnight in a humidified environment. The antibodies used in this study were anti-E-cadherin (CDH1) monoclonal antibody (#3195, dilution 1:100, Cell Signaling Technology), anti-N-cadherin (CDH2) monoclonal antibody (33-3900, dilution 1:100, Invitrogen), and anti-vimentin (VIM) monoclonal antibody (#5741, dilution 1:300, Cell Signaling Technology). After incubation with secondary antibodies for 20 min at room temperature, the reactions were visualized using the VECTASTAIN® Elite® ABC Kit (PK-6100, Vector Laboratories), which stains the targeted antigen brown, and counterstaining with hematoxylin. The degree of CDH1, CDH2 and VIM staining was evaluated based on the extent of membranous staining. Five stained hot spots were evaluated (n = 3), and the average of each stained area was estimated by microscopy (BZ-X 710; Keyence).
TUNEL assay
TUNEL assays were used to investigate apoptosis in formalin-fixed, paraffin-embedded xenograft tumor tissue samples. In brief, paraffin sections (2.5 µm) were deparaffinized in xylene and rehydrated with 70% and 90% alcohol. The TUNEL signal was detected using the ApopTag Fluorescein In Situ Apoptosis Detection Kit (S7110, Sigma–Aldrich). Nuclei were counterstained using VECTASHIELD Mounting Medium with DAPI (Vector Laboratories). Green fluorescence from apoptotic cells was analyzed by fluorescence microscopy (BZ-X 710; Keyence).
Statistical analysis
The Mann–Whitney U test, the χ2 test, or Student’s t test was used to compare patient characteristics. The disease-specific survival (DSS), including cancer or recurrence death, and recurrence-free survival (RFS) rates from November 2006 to December 2011 were calculated from the date of random assignment, validated by the Kaplan–Meier method, and compared with the log-rank test on an intent-to-treat basis; the corresponding HRs were calculated with the 95% CIs. Cox proportional hazards regression models were used to identify variables significantly associated with prognosis. Continuous variables are expressed as the mean ± SD, unless otherwise stated. Statistical significance was indicated at a p value < 0.05. All analyses were performed using JMP® 14 (SAS Institute Inc., Cary, NC, USA)