Summary
Oral squamous cell carcinoma (OSCC) become a heavy burden of public health, with approximately 300 000 newly diagnosed cases and 145 000 deaths worldwide per year. Nucleotide metabolism fuel DNA replication and RNA synthesis, which is indispensable for cell proliferation. But how tumor cells orchestrate nucleotide metabolic enzymes to support their rapid growth is largely unknown. Here we show that expression of pyrimidine metabolic enzyme dihydroorotate dehydrogenase (DHODH) is upregulated in OSCC tissues, compared to non-cancerous adjacent tissues. Enhanced expression of DHODH is correlated with a shortened patient survival time. Inhibition of DHODH by either shRNA or selective inhibitors impairs proliferation of OSCC cells and growth of tumor xenograft. Further, loss of functional DHODH imped de novo pyrimidine synthesis, and disrupt mitochondrial respiration probably through destabilizing the MICOS complex. Mechanistic study shows that transcriptional factor SOX2 plays an important role in the upregulation of DHODH in OSCC. Our findings add to the knowledge of how cancer cells co-opt nucleotide metabolism to support their rapid growth, and thereby highlight DHODH as a potential prognostic and therapeutic target for OSCC treatment.
Subject terms: Oral cancer, Metabolomics
Introduction
Head and neck cancer squamous cell carcinoma (HNSCC) ranks as the sixth-most common cancer worldwide, with an annual incidence of greater than 650 000 cases and causing 350 000 deaths annually.1 Oral squamous cell carcinoma (OSCC) is the most common type of HNSCC which refers to cancers that exist in the oral cavity, such as the mucosa of the lips, floor of the mouth, tongue, buccal mucosa, lower and upper gingival, hard palate, and retromolar trigone. Approximately 300 000 cases of OSCC occur every year, resulting in nearly 145 000 deaths worldwide.2,3 Development of OSCC is associated with multiple risk factors, among which tobacco use, alcohol consumption, and betel quid chewing are the most noteworthy. Other possible risk factors include oral microbiome, immune status, environmental pollutants, occupational exposures, heritable conditions.2 Though multiple therapeutic strategies, such as surgery, chemotherapy, radiotherapy, immunotherapy, or a combination of these modalities have been applied in OSCC treatment, lymph node metastasis and postoperative recurrence lead to a particularly poor prognosis and low 5-year survival rate.3,4 Therefore, identifying the key factors responsible for OSCC development and progression at molecular, genetic, and epigenetic levels is essential for improving the survival rate of OSCC patients.
Pyrimidines, including cytosine (C), thymine (T), and uracil (U), are essential building blocks for DNA replication and RNA transcription.5 The synthesis of pyrimidines is accomplished via two major pathways: de novo pyrimidine synthesis and salvage pathway. In mammalian cells, the de novo pyrimidine synthesis is the main source for producing pyrimidine nucleotides.6 CO2, glutamine, water, and aspartate are used as substrates to produce dihydroorotate and subsequently orotate. Orotate is then attached to the ribose-5-phosphate ring to form orotidine 5′-monophosphate (OMP), which is the precursor to generate cytosine, thymine, or uracil nucleotides (Fig. 1a).7 In this pathway, DHODH, a ubiquitously expressed mitochondrial inner membrane protein, governs the rate-limiting step by catalyzing the conversion of dihydroorotate to orotate through an oxidation reaction.6,8,9 It is widely established that nucleotide synthesis is frequently enhanced in cancer cells, to meet the increased demand for rapid cell proliferation.10 However, how DHODH is modulated in OSCC, and the role of DHODH in the regulation of OSCC cell proliferation remains elusive.
In this study, we demonstrate that the expression of DHODH was markedly increased in OSCC tissues, and was associated with multiple clinical parameters. Inhibition of DHODH by either shRNAs or selective antagonists reduced the proliferation rate and tumorigenicity of OSCC cell lines. Loss of DHODH dampens pyrimidine nucleotide synthesis and impairs mitochondrial homeostasis. Further, mechanistic study shows that the aberrant expression of DHODH in OSCC is regulated at transcriptional level by transcription factor SOX2.
Results
Expression of DHODH is upregulated in OSCC
We set out by examining the expression of DHODH in OSCC (n = 75) and non-cancerous adjacent oral mucous tissues (n = 30) through immunostaining. As shown in Fig. 1b, DHODH immunoreactivity was mainly observed in the perinuclear region, which was in line with previous reports showing mitochondrial localization of this protein.11 Notably, strong or moderate staining of DHODH was detected in tumor tissues; in contrast, most of non-cancerous tissues showed relatively weak DHODH signal (P < 0.001).
We further examined the correlation of DHODH with a set of clinicopathologic parameters. In contrast to those tumors diagnosed as Grade 1, DHODH staining was found to be more intense in the tumors of Grade 2 or 3 (Grade 2 vs. Grade 1, P = 0.018 1; Grade 3 vs. Grade 1, P = 0.012 9; Fig. 1c). No apparent differences in DHODH expression were found between Grade 2 or 3 tumors (Fig. 1c). Further, DHODH expression was substantially enhanced in OSCC patient with lymph node metastasis (P = 0.046 4, Fig. 1d). Nevertheless, no palpable correlation was observed between DHODH expression and tumor locations, including tongue, gums, pharynx, jaw, or lower lip (Fig. 1e).
To evaluate the impact of DHODH expression on the outcome of OSCC patients, we employed the online engine Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/detail.php) to analyzed 467 cases, whose gene expression profiling data is archived in The Cancer Genome Atlas (TCGA) or Genotype-Tissue Expression (GTEx). The survival curves of patients with high (n = 233) or low (n = 234) DHODH expression were established by Kaplan–Meier method, and statistical comparison of these curves was performed by Log-rank test. As results, those patients harboring a compromised DHODH expression showed a median survival time for approximately 67 months. While a largely shortened survival time was found for those patients with high DHODH expression (P = 0.035, Fig. 1f). These results strongly suggest that the expression of DHODH is upregulated in OSCC, and is corelated with the poor outcome of patients.
Inhibition of DHODH impeded OSCC cell proliferation
To explore the role of DHODH in OSCC development, we examined the impact of DHODH expression on OSCC cell proliferation. In line with the results from clinical samples, DHODH expression in OSCC CAL-27, HSC-3, and HN4 cell lines are much higher than human normal oral keratinocytes (HOK) (Fig. 2a). Since HSC-3 cells showed the highest expression level, we thus knockdown DHODH expression in this cell line by shRNAs (Fig. 2b). As results, loss of DHODH markedly reduced the cell viability by nearly 70%, shown by CCK8 assay and Edu incorporation assay (Fig. 2c–d). Consistently, in colony formation assay, the average number of clones formed by untreated HSC-3 cells was 85.33 ± 3.18, while treatment of shRNA #1 or #2 reduced the clone number to 13 ± 3.46 or 11 ± 3.06, respectively (Fig. 2e). These changes in OSCC cell proliferation were not likely due to the off-target effects, since these two shRNAs targeted distinct sites on DHODH mRNA. Similar results were obtained from HN4 cells (Fig. S1a–c).
Next, we tested whether enzymatic activity is prerequisite for DHODH-dependent OSCC cell proliferation. BAY 2402234 was reported as a selective potent antagonist of DHODH, through binding the ubiquinone binding site of DHODH between the N-terminal helices.12 Indeed, incubation with BAY 2402234 obliterate roughly 90% of DHODH activity in the HSC-3 cell lysates (Fig. 2f). Notably, BAY 2402234 treatment markedly abolished the light absorption values in CCK8 assay (Fig. 2g), the percentage of Edu-positive cells in Edu incorporation assay (Fig. 2h), and the clone numbers in colony formation assay (Fig. 2i). Similar results were observed by the treatment with another DHODH inhibitor Vidofludimus (Fig. 2f–i).13
To estimate the role of DHODH in regulating OSCC tumor growth, we established a mouse xenograft model by subcutaneous injection of HSC-3 cells that were stably expressed with non-target control (NC) shRNA or DHODH shRNAs. In spite of no obvious differences during the first 8 days after injection, the tumors formed by both DHODH shRNA-expressed cells grew much slower than the NC control tumors (Fig. 2j). Accordingly, the number of cells with Ki-67 expression was markedly decreased upon loss of DHODH (Fig. 2k). Consistently, inhibition of DHODH by BAY 2402234 or Vidofludimus treatment largely retarded xenograft growth and abolished Ki-67-positive cells (Fig. 2l–m). These results suggest that DHODH expression was required for OSCC cell proliferation.
Inhibition of DHODH impairs pyrimidine nucleotide synthesis
Considering that cancer cells frequently hijack nucleotide metabolism to boost cell proliferation,14 we examined the role of DHODH in cellular de novo pyrimidine synthesis in OSCC cells. As shown in Fig. 3a, b, loss of DHODH resulted in a compromised level of orotate and an augmented level of dihydroorotate in HSC-3 cells. To further monitor the metabolic flux of this pathway, we incubated HSC-3 cells with 15N-labeled glutamine to trace the metabolites generated from glutamine (Fig. 3c). Indeed, shRNA-mediated silence of DHODH substantially eliminated glutamine-derived orotate in HSC-3 cell lysates, with an obvious accumulation of glutamine-derived dihydroorotate, suggesting a blockage in conversion of dihydroorotate to orotate (Fig. 3d, e). Accordingly, the levels of 15N-labeled downstream metabolites, including uridine monophosphate (UMP), uridine triphosphate (UTP), and cytidine triphosphate (CTP), were dampened in the DHODH shRNAs-expressed cells (Fig. 3f–h, and S1d).
Salvage pathway support cellular pyrimidine pool by recycling nucleosides and free bases generated by DNA and RNA breakdown, and functions as an alternative route when de novo synthesis pathway is retarded.15 However, DHODH-silenced HSC-3 cells exhibited a 60–70% reduction in the level of total cellular UMP, UTP, or CTP, compared to NC shRNA-treated cells (Fig. 3i–k). These data suggest that DHODH is required for pyrimidine nucleotide synthesis in OSCC cells, which cannot be compensated by the salvage pathway.
Inhibition of DHODH destabilizes mitochondrial contact site and cristae organizing system (MICOS) complex and disrupts mitochondrial homeostasis
DHODH protein is mainly located in the mitochondrial inner membrane.11 Strikingly, the amount of mitochondrial (mt)DNA was largely reduced in DHODH-silenced HSC-3 cells (Fig. 4a). Therefore, we examined the potential involvement of DHODH in mitochondrial homeostasis. As shown in Fig. 4b–c, knockdown of DHODH diminished both the ATP content and oxygen consumption rate (OCR) in mitochondria isolated from HSC-3 cells. As normalized to the total mass of mitochondrial protein, the lowered mitochondrial ATP and OCR was supposedly due to mitochondrial dysfunction, rather than the reduced amounts of mitochondria. Additionally, since loss of DHODH caused accumulation of cellular dihydroorotate (Fig. 3b), we hypothesized that the treatment effect of dihydroorotate might recapitulate the effects of DHODH shRNA. Indeed, incubation with dihydroorotate resulted in an obliterate mtDNA as well as a repressed mitochondrial ATP and OCR (Fig. 4d–f).
Mitochondrial ultrastructure analyses by transmission electron microscopy (TEM) revealed a reduced and disturbed structure of mitochondrial cristae in DHODH-silenced cells (Fig. 4g). Mitochondrial cristae are folds of the inner membrane, which provides the indispensable platform for respiratory chain.16 MISCO complex, formed by a set of MIC protein components, plays a pivotal role in maintaining the proper structure of mitochondrial cristae.17 Of our particular interest, we tested the impact of DHODH on MISCO complex assembly. We found that either knockdown of DHODH (Fig. 4h and S1e) or treatment with dihydroorotate (Fig. 4i) partially abolished the amount of MIC27, MIC19, or MIC10 proteins in MIC60 immunoprecipitates, suggesting that loss of DHODH may destabilize the MISCO complex. Further, we knockdown the endogenous DHODH in HSC-3 cells and reconstituted the expression of an enzymatic-dead R135C mutant,18 and found that this mutant similarly attenuated the interactions between MIC protein components (Fig. 4j–k), suggesting the enzymatic activity of DHODH was required for maintaining MISCO complex. Together, these results suggest that inhibition of DHODH causes mitochondrial dysfunction probably through disturbing MISCO complex assembly.
Sex determining region Y-box 2 (SOX2) transcriptionally regulates DHODH expression in OSCC
Next, we turned to explore the mechanism responsible for DHODH upregulation in OSCC. Quantification of DHODH mRNA showed an apparent increased in CAL-27, HSC-3, and HN4 cells compared to HOK cells (Fig. 5a), elucidating that DHODH was likely regulated at transcriptional level. Analyses of DHODH promoter region (−1500 – 0) by PROMO online engine (Version 8.3, http://alggen.lsi.upc.es/cgi-bin/promo_v3/promo/promoinit.cgi?dirDB=TF_8.3) revealed putative binding sites for three OSCC-associated transcriptional factors, c-Jun (ref. 19), SOX2 (ref. 20), and c-Myc (ref. 21) (Fig. 5b). However, only knockdown of SOX2 largely abolished DHODH expression in HSC-3 cells (Fig. 5c–d). Further, transfection of various amounts of DHODH-expressing vector caused an enhanced DHODH expression at a distinct extent in HOK cells (Fig. 5e), suggesting that SOX2, but not c-Jun or c-Myc, was involved in the regulation of DHODH transcription. Consistently, modulation of SOX2 expression in HOK (Fig. 5f) or HSC-3 cells (Fig. 5g) efficiently affected the transcriptional activity of DHODH, shown by a luciferase-based reporter gene assay.
We determined the correlation between the expression of DHODH and SOX2. As shown in Fig. 5h, the expression pattern of SOX2 among HOK and OSCC cell lines was comparable to DHODH expression. Further, immunostaining of SOX2 revealed that most cases with high DHODH expression showed a stronger SOX2 immunoreactivity (Fig. 5i). Quantification of the staining showed that these correlations were significant. Taken together, these results suggest that upregulation of DHODH in OSCC is regulated by SOX2 at transcriptional level.
Discussion
Nucleotide metabolism fuel DNA replication and RNA synthesis, which is indispensable for cell proliferation. But how tumor cells orchestrate nucleotide metabolic enzymes to support their rapid growth is largely unknown.22 Here, we show that pyrimidine metabolic enzyme DHODH was overexpressed in OSCC, and is correlated with poor outcome of OSCC patients. Further, inhibition of DHODH by either shRNA or selective inhibitor counteracted OSCC cell proliferation and xenograft growth. Our data suggest that DHODH might be a prognostic biomarker or potential treatment target for OSCC.
We demonstrate that loss of DHODH blocked pyrimidine synthesis, and abolished cellular pyrimidine pool, which was unlikely replenished by salvage pathway. Tumor cells frequently harbor an enhanced de novo synthesis of pyrimidine nucleotides to support the rapid cell proliferation.23 Pyrimidine nucleotides are essential building blocks for DNA replication and RNA transcription.24 Additionally, CTP is prerequisite for phosphatidylcholine synthesis, which is the major component of cell membrane.25 While UTP has a role in the metabolism of galactose, where the activated form UDP-galactose is converted to UDP-glucose.26 Considering DHODH governs the rate-limiting step for pyrimidine synthesis, it is reasonable to infer that inhibition of DHODH impairs OSCC cell proliferation by disrupting pyrimidine nucleotide supply.
It was an unexpected finding that knockdown of DHODH defaced MICOS complex and perturbed the cristae structure in mitochondria. Mitochondria are essential organelles within all proliferating cells, including tumor cells. Mitochondria are responsible for nutrients transport, energy metabolism, and synthesis of numerous cellular components required for growth and proliferation.24 These functions of mitochondria are largely dependent on cristae formed by the folding of the inner membrane, whose structure increases the surface area of the inner membrane. MICOS protein complexes are crucial for the formation and maintenance of the cristae structure.16 Nine types of MICOS complex subunits have been described so far, including MIC 60 (mitofilin), MIC 10 (MINOS1), MIC 19 (CHCHD3), MIC 25 (CHCHD6), MIC 26 (APOO), MIC 23, MIC 27 (ApooL), MIC 13 (Qil1), and MIC 14 (CHCHD10). Accumulating evidences indicate that MIC60 is the core component of MICOS. MIC 60, mitochondrial outer membrane protein Sam 50 and MIC 19 interact with each other to form bridge-like complex, which mediates the contact between the mitochondrial outer and inner membranes.27,28 Considering that mitochondrial respiration complexes are assembled and located on the inner membrane, it is reasonable to infer that impaired mitochondrial respiration, observed in DHODH-silenced cells, was probably caused by the modulation of the MICOS complex. In addition, since we also show that treatment with dihydroorotate caused similar effects on mitochondria, further work is still needed to decipher the biological consequences for dihydroorotate overload in mitochondrial inner membrane.
Further, our data suggest that SOX2 probably governs DHODH expression in OSCC. We show that both expression or transcription activity of DHODH in OSCC cells could be manipulated by modulating the expression SOX2. SOX2 plays an important role in the regulation of embryonic cell development and the stem-cell nature of various adult stem cell populations.29 However, SOX2 is frequently activated in tumor cells with stem cell-like characteristics, regulating cell self-renewal, maintaining cell pluripotency, and thereby promoting tumor growth. The abnormal expression of SOX2, probably due to increased copy number in genome, is associated with multiple tumor types, including OSCC.30 It is reported that upregulation of SOX2 is associated with a more aggressive phenotype of OSCC cells, and rendered OSCC cells resistant to chemotherapy.31 In this study, our data underscores the role of SOX2 in the reprogramming of pyrimidine synthesis, through depicting its regulation of DHODH transcription. More work is merited to determine whether SOX2 also regulates other steps of nucleotide synthesis, such as the pentose phosphate pathway or salvage pathway.
Materials and Methods
Antibodies recognizing Tubulin, COX IV, SOX2, c-Myc, were obtained from Cell Signaling Technology. Antibodies recognizing DHODH, and MIC19 were obtained from Abcam. Antibody against Flag, anti-Flag M2 agarose beads, dihydroorotate, and bovine serum albumin were purchased from Sigma. Antibody recognizing MIC60 was obtained from NOVUS. Antibody recognizing MIC10 was obtained from PROSPEC. Antibody recognizing MIC27, horseradish peroxidase-conjugated goat anti-mouse, or rabbit secondary antibodies were obtained from Thermo Fisher Scientific. [15N2]-glutamine was purchased from Cambridge Isotope Laboratories. BAY 2402234 and Vidofludimus were obtained from TOPSCIENCE. Lipofectamine 2000 was obtained from Thermo Fisher Scientific.
DNA constructs and mutagenesis
PCR-amplified human DHODH, MIC60, SOX2, and DHODH promoter region (−1 500-0) were subcloned into pcDNA3.1/hygro(+)-Flag, or pGL3 vectors. DHODH R135C mutant was constructed using the QuikChange site-directed mutagenesis kit (Stratagene, La Jolla, CA).
The following pGIPZ shRNAs were used: control shRNA, GCT TCT AAC ACC GGA GGT CTT; DHODH shRNA #1, AAT ATT CAA TGT CCT TGC A (targeting non-coding region); DHODH shRNA #2, TCC ATG ACT TTT TCC TCC T (targeting non-coding region); SOX2 shRNA #1, TCT CAG CTT ATA AAC AAT G; SOX2 shRNA #2, AAA ACA TTT TTT TCG TCG C. Mixed siRNAs targeting c-Myc or c-Jun was obtained from Santa Cruz.
Cell culture
CAL-27 cells were obtained from ATCC, and HSC-3 cells were obtained from JCRB Cell Bank. HOK cells were obtained from ScienCell. HN4 cells were provided by J. Silvio Gutkind. These cells maintained with Dulbecco’s Modified Eagle Medium supplemented with 10% fetal bovine serum. Cells were treated with BAY 2402234 at 10 nM or Vidofludimus at 10 μmol·L−1.
For generating shRNA-depleted stable cell lines, cells were transfected with shRNA plasmids and selected by puromycin. For generating gene-expressing stable cell lines, shRNA-depleted cells were infected with lentivirus carrying WT or mutant gene and selected by hygromycin or/and G418.
Transfection of plasmids were performed by using Lipofectamine 2000, following the manufacture’s instructions.
Immunoprecipitation and immunoblot analysis
Protein was extracted from cultured cells using a lysis buffer (50 mmol·L−1 Tris-HCl [pH 7.5], 0.1% SDS, 1% Triton X-100, 150 mmol·L−1 NaCl, 1 mmol·L−1 dithiothreitol, 0.5 mmol·L−1 EDTA, 100 μmol·L−1 PMSF, 100 μmol·L−1 leupeptin, 1 μmol·L−1 aprotinin, 100 μmol·L−1 sodium orthovanadate, 100 μmol·L−1 sodium pyrophosphate, and 1 mmol·L−1 sodium fluoride). The cell lysates was centrifuged at 13 400 g, and supernatants (2 mg protein per mL) were immunoprecipitated overnight at 4 °C using the indicated antibodies. Following the overnight incubation, protein A or G agarose beads were added and left for an additional 3 h. The immunocomplexes were then washed with lysis buffer 3 times and prepared for immunoblot analyses with corresponding antibodies as described previously.
Analysis of intermediate metabolites
Metabolites were confirmed by high-resolution mass spectrometry (HRMS). For preparation of cell sample, 5 × 107 cells were washed with ice-cold PBS to remove culture medium. The tissue samples (500 mg) of mice intestine were homogenized on ice immediately after sampling. Metabolites were extracted by treatment with 90/9/1 (v/v/v) acetonitrile/water/formic acid. The samples were centrifuged at 117 000 g for 10 min, and then the supernatant was transferred to a clean tube and evaporated to dryness under a nitrogen gas stream. Dried samples were reconstituted in 0.2% ammonium hydroxide in ammonium acetate (10 mmol·L−1), then 10 μL was injected into a Thermo Scientific Vanquish liquid chromatography system containing a Thermo Hypercarb 100 × 3 mm 3 μm HPLC column heated to 35 °C with mobile phase A (MPA) consisting of 0.2% ammonium hydroxide in ammonium acetate (10 mmol·L−1) and mobile phase B (MPB) consisting of 0.2% ammonium hydroxide in acetonitrile. The gradient elution was performed at a flow rate of 0.3 mL·min−1: 0 min (0% MPB)-2.0 min (0% MPB)-15.0 min (30% MPB)-15.1 min (95% MPB)-20.0 min (95% MPB)-20.1 min (0% MPB)-25.0 min (STOP). Thermo Orbitrap Fusion Tribrid Mass Spectrometer in Selected Ion Mode (SIM) electrospray positive mode was used to acquire data. Then, use Thermo TraceFinder software to qualitatively analyze the Peak integration and area.
Cells were incubated with [15N2]-glucose (4 mmol·L−1), for indicated time to track metabolic flux. The metabolites were analyzed by HRMS with similar protocol.
EdU incorporation assay
EdU incorporation assay was performed by using EdU Proliferation Kit (Abcam). Briefly, after treatment, cells were incubated with 20 µmol·L−1 EdU for 3 h. Cells were then sequentially incubated with 1X Fixative Solution for 15 min, 1X Permeabilization Buffer for 20 min, EdU Additive Solution Reaction Buffer for 30 min, and 5 µg·mL−1 Hoechst for 30 min. The EdU- or Hoechst-positive cells were counted in at least 20 fields. The percentage of Edu-positive cells was calculated as the ratio between EdU- and Hoechst-positive cells.
Colony formation assay
Cells were seeded into a six-well plate at a density of 100 cells per well, and cultured for 14 days. Colonies were fixed and stained with Crystal Violet solution (0.5% Crystal violet in 20% ethanol), colonies of more than 50 cells were counted.
CCK8 assay
Cell Counting Kit 8 (WST-8/CCK8) (Abcam) assessed Cell proliferation. Cells were seeded in a 96 well plate (5 000 cells per well) overnight. A total of 10 µL of the cell proliferation reagent WST-8 solution was added to each well of the 96-well plate and then incubated for 3 h at 37 °C in dark. The change in absorbance was measured at 460 nm.
Immunohistochemical staining
Immunohistochemical staining was performed according to manufacturer instructions using the VECTASTAIN ABC kit (VECTOR LABORATORIES, CA). Sections were scored according to the percentage of positive cells and staining intensity. Scoring was established as follows: the percentage of positive cells was scored as 0 if 0% of the tumor cells showed positive staining, 1 if 0% to 10%, 2 if 11% to 30%, 3 if 31% to 70%, and 4 if 71% to 100%, the staining intensity was rated on a scale of 0 to 4: 0, negative; 1, weak; 2, moderate; 3, strong; 4, very strong. We then multiply the proportion and intensity scores to obtain a total score (range: 0–16).
Transmission electron microscopy (TEM)
The cells were fixed with 0.1% glutaraldehyde in 0.1 mol·L−1 sodium iodate and 1% OsO4 successively for 2 h and 1.5 h, respectively. Then the samples were stained in 3% uranyl acetate aqueous solution for 1 h after washing. After that, the samples were rinsed again with water, then dehydrated with gradient alcohol (50%, 75%, and 95%–100% alcohol), and finally embedded in Epon-Araldite resin (Canemco, 034). Ultrathin sections were prepared on Reichert Ultramicrotome, and the sections were counterstained with 0.3% lead citrate, and then observed under a Philips EM420 transmission electron microscope. Use Image Pro Plus version 3 to obtain the value of the area occupied by autophagic vesicles and cytoplasm.
Reverse transcription polymerase chain reaction (RT-PCR)
RT-PCR was performed following a previous report.32 The primers are as follows: DHODH forward, GCT GTC ATT AAC AGG TAT GGA TTT AAC AG; DHODH reverse, GTT GAT AAA TCC CGG AGG GGC TTC.
Measurement of mitochondrial DNA content
QIAamp DNA mini kit (QIAGENE, Germantown, MD) was used to extract total DNA according to the manufacturer’s instructions. The level of mitochondrial DNA D-loop structure was measured by qRT-PCR and standardized to β-actin encoded by genomic DNA. The primers are as follows: D-Loop forward, GAT TTG GGT ACC ACC CAA GTA TTG; D-Loop reverse, GTA CAA TAT TCA TGG TGG CTG GCA; Actin forward, TCA CCC ACA CTG TGC CCA TCT ACG A; Actin reverse, CAG CGG AAC CGC TCA TTG CCA ATG G.
Measurement of DHODH activity
DHODH activity in cell lysates were measured following the previous report.33 A total of 300 μL cell lysate was mixed with in an aqueous solution (total volume, 1.0 mL) containing 500 μmol·L−1 DHO, 200 mmol·L−1 K2CO3-HCI (pH 8.0), 0.2% triton X-100, and 100 μmol·L−1 coenzyme Q10 at 37 °C 30 min. An aliquot (100 μL) of the mixture of enzyme reaction mixture was mixed with 100 μL of 0.5 μmol·L−1orotic acid, 50 μL of H2O, 250 μL of 4.0 mmol·L−1 4-TFMBAO, 250 μL of 8.0 mmol·L−1 K3[Fe(CN)6], and 250 μL of 80 mmol·L−1 K2CO3 and then heated at 80 °C for 4.0 min. The reaction was stopped by cooling in an ice-water bath and the FL intensity was measured with a spectrofluorometer (FP-6300 Jasco, Tokyo, Japan): excitation and emission wavelengths were 340 nm and 460 nm, respectively.
Mitochondria isolation
Mitochondria was isolated by using Mammalian Mitochondria Isolation Kit (Biovision), following the manufacture’s instruction.
Measurement of ATP content
Mitochondrial ATP content was measured in freshly isolated mitochondrial by using ATP Assay Kit (Colorimetric/Fluorometric, Abcam), following the manufacture’s instruction.
Measurement of mitochondrial OCR
Mitochondrial OCR was measured in freshly isolated mitochondrial by using MitoCheck® Mitochondrial OCR Assay Kit, following the manufacture’s instruction.
Mice xenograft
A total of 1 × 106 HSC-3 cells was subcutaneously implanted into 6-week-old female athymic nude mice (n = 7). The tumor volume was calculated by the formula: 0.5 × length × width2. The mice were sacrificed 27 days after implantation, and tumor tissues were subjected to immunohistochemical staining.
Quantification and statistical analysis
Statistical analysis was performed using the two-tailed unpaired Student’s t-test with a P-value of < 0.05 was considered as statistical significance unless specifically indicated. All data were analyzed as the mean ± standard deviation from three independent experiments/samples unless otherwise stated. Significance levels are: *P < 0.05; **P < 0.01; ***P < 0.001.
Supplementary information
Acknowledgements
This work was supported by grant from the CAMS Innovation Fund for Medical Sciences (CIFMS) 2019-I2M-5-004 (Q.C.); National Natural Science Foundation of China grants 81672674 (R.L.), 81872218 (R.L.) and 81872208 (L.J.); Science & Technology Department of Sichuan Province Applied Basic Research Program 2020YJ0451 (T.R.); Fok Ying Tong Education Foundation grant 161036 (R.L.); Young Talent Program of China National Nuclear Corporation CNNC201948 (S.J.); Chengdu Medical College Fund Natural Science General Project CYZ18-17 (S.J.).
Author contributions
Q.C., R.L., and L.J. conceived and designed the study; X.Q., S.J., T.R., Y.X., and Y.H. performed the experiments; and R.L. and X.Q. wrote the paper, with comments from all authors.
Competing interests
The authors declare no competing interests.
Footnotes
These authors contributed equally: Xuemei Qiu, Sheng Jiang
Contributor Information
Tao Ren, Email: rentao509@cmc.edu.cn.
Lu Jiang, Email: jianglu@scu.edu.cn.
Rui Liu, Email: liurui_scu@hotmail.com.
Supplementary information
The online version of this article (10.1038/s41368-020-00109-x) contains supplementary material, which is available to authorized users.
References
- 1.Chandler KB, et al. β-Catenin/CBP inhibition alters epidermal growth factor receptor fucosylation status in oral squamous cell carcinoma. Mol. Omics. 2020;16:195–209. doi: 10.1039/D0MO00009D. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 2.Chi AC, Day TA, Neville BW. Oral cavity and oropharyngeal squamous cell carcinoma–an update. CA Cancer J. Clin. 2015;65:401–421. doi: 10.3322/caac.21293. [DOI] [PubMed] [Google Scholar]
- 3.Kato MG, et al. Update on oral and oropharyngeal cancer staging - International perspectives. World J. Otorhinolaryngol. Head. Neck Surg. 2020;6:66–75. doi: 10.1016/j.wjorl.2019.06.001. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 4.Chen X, et al. Zymosan promotes proliferation, Candida albicans adhesion and IL-1β production of oral squamous cell carcinoma in vitro. Infect. Agent Cancer. 2020;15:51. doi: 10.1186/s13027-020-00315-6. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 5.Abbas H. First-principles study of interaction of serine with nucleobases of DNA and RNA. J. Biol. Phys. 2017;43:105–111. doi: 10.1007/s10867-016-9436-1. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 6.Löffler M, Jöckel J, Schuster G, Becker C. Dihydroorotat-ubiquinone oxidoreductase links mitochondria in the biosynthesis of pyrimidine nucleotides. Mol. Cell Biochem. 1997;174:125–129. doi: 10.1023/A:1006859115450. [DOI] [PubMed] [Google Scholar]
- 7.Zhang Y, et al. Transition state analogues of Plasmodium falciparum and human orotate phosphoribosyltransferases. J. Biol. Chem. 2013;288:34746–34754. doi: 10.1074/jbc.M113.521955. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 8.Löffler, M., Carrey, E. A. & Knecht, W. The pathway to pyrimidines: the essential focus on dihydroorotate dehydrogenase, the mitochondrial enzyme coupled to the respiratory chain. Nucleosides Nucleotides Nucleic Acids39, 1–25 (2020). [DOI] [PubMed]
- 9.Yin S, Kabashima T, Zhu Q, Shibata T, Kai M. Fluorescence assay of dihydroorotate dehydrogenase that may become a cancer biomarker. Sci. Rep. 2017;7:40670. doi: 10.1038/srep40670. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 10.Villa, E., Ali, E. S., Sahu, U. & Ben-Sahra, I. Cancer Cells Tune the Signaling Pathways to Empower de Novo Synthesis of Nucleotides. Cancers (Basel)11, 10.3390/cancers11050688 (2019). [DOI] [PMC free article] [PubMed]
- 11.Oliver JD, et al. F901318 represents a novel class of antifungal drug that inhibits dihydroorotate dehydrogenase. Proc. Natl Acad. Sci. USA. 2016;113:12809–12814. doi: 10.1073/pnas.1608304113. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 12.Christian S, et al. The novel dihydroorotate dehydrogenase (DHODH) inhibitor BAY 2402234 triggers differentiation and is effective in the treatment of myeloid malignancies. Leukemia. 2019;33:2403–2415. doi: 10.1038/s41375-019-0461-5. [DOI] [PubMed] [Google Scholar]
- 13.Muehler A, Kohlhof H, Groeppel M, Vitt D. Safety, Tolerability and Pharmacokinetics of Vidofludimus calcium (IMU-838) After Single and Multiple Ascending Oral Doses in Healthy Male Subjects. Eur. J. Drug Metab. Pharmacokinetics. 2020;45:557–573. doi: 10.1007/s13318-020-00623-7. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 14.Santoro V, et al. SLC25A32 sustains cancer cell proliferation by regulating flavin adenine nucleotide (FAD) metabolism. Oncotarget. 2020;11:801–812. doi: 10.18632/oncotarget.27486. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 15.Squadrito F, et al. Pharmacological activity and clinical use of PDRN. Front Pharm. 2017;8:224. doi: 10.3389/fphar.2017.00224. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 16.Kozjak-Pavlovic V. The MICOS complex of human mitochondria. Cell Tissue Res. 2017;367:83–93. doi: 10.1007/s00441-016-2433-7. [DOI] [PubMed] [Google Scholar]
- 17.Rampelt H, Zerbes RM, van der Laan M, Pfanner N. Role of the mitochondrial contact site and cristae organizing system in membrane architecture and dynamics. Biochim Biophys. Acta Mol. Cell Res. 2017;1864:737–746. doi: 10.1016/j.bbamcr.2016.05.020. [DOI] [PubMed] [Google Scholar]
- 18.Fang J, et al. Dihydro-orotate dehydrogenase is physically associated with the respiratory complex and its loss leads to mitochondrial dysfunction. Biosci. Rep. 2013;33:e00021. doi: 10.1042/BSR20120097. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 19.Turatti, E., da Costa Neves A., Fau - de Magalhães, M. H. C. G., de Magalhães Mh Fau - de Sousa, S. O. M. & de Sousa, S. O. Assessment of c-Jun, c-Fos and cyclin D1 in premalignant and malignant oral lesions. [DOI] [PubMed]
- 20.Kokalj Vokač N, Cizmarević B, Zagorac A, Zagradišnik B, Lanišnik B. An evaluation of SOX2 and hTERC gene amplifications as screening markers in oral and oropharyngeal squamous cell carcinomas. Mol. Cytogenet. 2014;7:5. doi: 10.1186/1755-8166-7-5. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 21.Pallavi N, Nalabolu GRK, Hiremath SKS. Bcl-2 and c-Myc expression in oral dysplasia and oral squamous cell carcinoma: An immunohistochemical study to assess tumor progression. J. Oral. Maxillofac. Pathol. 2018;22:325–331. doi: 10.4103/jomfp.JOMFP_197_18. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 22.Robinson, A. D., Eich, M. L. & Varambally, S. Dysregulation of de novo nucleotide biosynthetic pathway enzymes in cancer and targeting opportunities. [DOI] [PubMed]
- 23.Yeo KL, et al. Synergistic suppression of dengue virus replication using a combination of nucleoside analogs and nucleoside synthesis inhibitors. Antimicrob. Agents Chemother. 2015;59:2086–2093. doi: 10.1128/AAC.04779-14. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 24.Cha YH, Yook JI, Kim HS, Kim NH. Catabolic metabolism during cancer EMT. Arch. Pharm. Res. 2015;38:313–320. doi: 10.1007/s12272-015-0567-x. [DOI] [PubMed] [Google Scholar]
- 25.Lütje S, et al. PSMA ligands for radionuclide imaging and therapy of prostate cancer: clinical status. Theranostics. 2015;5:1388–1401. doi: 10.7150/thno.13348. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Mok WK, Tan YX, Lee J, Kim J, Chen WN. A metabolomic approach to understand the solid-state fermentation of okara using Bacillus subtilis WX-17 for enhanced nutritional profile. AMB Express. 2019;9:60. doi: 10.1186/s13568-019-0786-5. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 27.Tang J, et al. Sam50-Mic19-Mic60 axis determines mitochondrial cristae architecture by mediating mitochondrial outer and inner membrane contact. Cell Death Differ. 2020;27:146–160. doi: 10.1038/s41418-019-0345-2. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28.Stephan T, et al. MICOS assembly controls mitochondrial inner membrane remodeling and crista junction redistribution to mediate cristae formation. EMBO J. 2020;39:e104105. doi: 10.15252/embj.2019104105. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 29.Novak, D. et al. SOX2 in development and cancer biology. Semin Cancer Biol. 10.1016/j.semcancer.2019.08.007 (2019). [DOI] [PubMed]
- 30.Freier K, et al. Recurrent copy number gain of transcription factor SOX2 and corresponding high protein expression in oral squamous cell carcinoma. Genes Chromosomes Cancer. 2010;49:9–16. doi: 10.1002/gcc.20714. [DOI] [PubMed] [Google Scholar]
- 31.Liang X, et al. Combined class I histone deacetylase and mTORC1/C2 inhibition suppresses the initiation and recurrence of oral squamous cell carcinomas by repressing SOX2. Cancer Lett. 2019;454:108–119. doi: 10.1016/j.canlet.2019.04.010. [DOI] [PubMed] [Google Scholar]
- 32.Oh JY, Park MA, Kim YC. Peppermint oil promotes hair growth without toxic signs. Toxicol. Res. 2014;30:297–304. doi: 10.5487/TR.2014.30.4.297. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 33.Yin, S., Kabashima, T., Zhu, Q., Shibata, T. & Kai, M. Fluorescence assay of dihydroorotate dehydrogenase that may become a cancer biomarker. Sci. Rep.7, 40670 (2017). [DOI] [PMC free article] [PubMed]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.