The Role of Peroxisome Proliferator-Activated Receptors (PPARs) in Disease and Targeted TreatmentsView this Special Issue
Research Article | Open Access
Rui Kong, Nan Wang, Wei Han, Wen Bao, Jie Lu, "Fenofibrate Exerts Antitumor Effects in Colon Cancer via Regulation of DNMT1 and CDKN2A", PPAR Research, vol. 2021, Article ID 6663782, 20 pages, 2021. https://doi.org/10.1155/2021/6663782
Fenofibrate Exerts Antitumor Effects in Colon Cancer via Regulation of DNMT1 and CDKN2A
Peroxisome proliferator-activated receptor alpha (PPARA) is the molecular target of fibrates commonly used to treat dyslipidemia and diabetes. Recently, the potential role of PPARA in other pathological conditions, such as cancers, has been recognized. Here, using bioinformatics analysis, we found that PPARA was expressed at relatively low levels in pancancers, and Kaplan-Meier analyses revealed that high PPARA protein expression was correlated with better survival of patients with colon cancer. In vitro experiments showed that fenofibrate regulated cell cycle distribution, promoted apoptosis, and suppressed cell proliferation and epithelial mesenchymal transition by activating PPARA. PPARA activation inhibited DNMT1 activity and abolished methylation-mediated CDKN2A repression. Downregulation of cyclin-CDK complexes led to the restoration of CDKN2A, which caused cell cycle arrest in the G1 phase via regulation of the CDKN2A/RB/E2F pathway. Finally, we demonstrated that fenofibrate administration inhibited tumor growth and DNMT1 activity in vivo. The PPARA agonist, fenofibrate, might serve as an applicable agent for epigenetic therapy of colon cancer patients.
Colorectal cancer (CRC) ranks third in terms of morbidity and fourth in terms of mortality. CRC is also the most common type of cancer worldwide, with almost 900,000 deaths each year [1, 2]. Although new treatment options, including immunotherapy and neoadjuvant chemotherapy, have significantly improved patient prognosis, the 5-year survival rate of CRC remains below 15% . Thus, investigations of the molecular mechanisms involved in cancer initiation are necessary to develop new therapeutic strategies.
In addition to genetic alterations (i.e., deletion, amplification, and translocation), epigenetic modifications play an important role in malignant progression of tumors. CpG islands (CGI) are contiguous groups of dinucleotides mainly located at the 5 end of a gene and are characterized by high GC content . Most CGIs in gene promoters are unmethylated, allowing active transcription . CGI methylation changes are hallmarks of many human cancers and lead to concomitant gene inactivation [6–8]. DNA methyltransferase 1 (DNMT1) is a major DNA methyltransferase responsible for methylation maintenance during DNA replication, and inactivation of DNMT1 in mice results in early embryonic lethality . DNMT3a and DNMT3b mainly act as de novo methyltransferases .
CRCs are characterized by lower levels of absolute genomic methylation compared with normal tissues, a characteristic that contributes to high genomic instability and results in cancer development . In addition, promoter hypermethylation of specific genes has been identified in CRCs, as well as methylated CGIs that are associated with gene silencing. Hypermethylation in several tumor suppressor genes such as RASSF1, PTEN, and CDKN2A is associated with abnormal cellular activities, including aberrant cell aging, proliferation, and death [12–14].
Peroxisome proliferator-activated receptor alpha (PPARA) is a ligand-activated transcription factor that belongs to the nuclear hormone receptor superfamily [15, 16]. Studies have demonstrated that PPARA plays a critical role in lipid metabolism and the inflammatory response [17, 18]. Fenofibrate is a selective PPARA agonist that regulates lipid transport and metabolism, and is widely used in the treatment of hyperlipidemia . In addition to the metabolic efficacy, recent studies have revealed the antitumor function of PPARs [20–24]. Studies have also confirmed the oncosuppressive effect of fenofibrate in various human cancer cell lines through different signaling pathways [25–27]. Reports also revealed the inhibitory role of PPARA on DNMT1 activity in mouse models . In the present study, we demonstrated that the PPARA agonist, fenofibrate, inhibited DNMT1-mediated methylation of CDKN2A and exerted anticancer effects by promoting cell apoptosis, inhibiting cell migration, and suppressing cell proliferation via the CDKN2A/RB/E2F pathway.
2. Materials and Methods
2.1. Cell Culture
The colorectal cancer cell line HCT116, the colon cancer cell lines SW480 and Caco-2, and the two normal intestine epithelial cell lines NCM460 and HIEC were obtained from Chinese Academy of Sciences Committee Type Culture Collection Cell Bank (Beijing, China). HCT116 was cultured in complete DMEM (Thermo Fisher Scientific, Waltham, MA), and the other cell lines were maintained in RPMI 1640 (Thermo Fisher Scientific, Waltham, MA) supplemented with 10% fetal bovine serum (FBS, Thermo Fisher Scientific, Waltham, MA) and 1% penicillin-streptomycin.
Inhibition of DNMT1 activity was conducted using 5 μM 5-azacytidine (MedChemExpress, Shanghai, China), and the treatment time was 24 h. The pcDNA DNMT1 was transfected into HCT116 and SW480 cells for upregulation of DNMT1 expression. The pcDNA DNMT1 were synthesized by GenePharma Co. Ltd. (Shanghai, China).
Cells were exposed to different concentrations of fenofibrate (Topscience Co., Ltd., Shanghai, China) for 48 hours, and morphology changes were observed using a phase-contrast microscope and imaged (200x).
2.2. Cell Viability Assay
Colon cancer cell lines SW480 and the colorectal cancer cell line HCT116 were plated at a density of 5000 cells/mL in 96-well plates (100 μL medium per well) with three replicates. Cells were treated with fenofibrate in the pharmacologic concentration range 0–260 μM for 24 h and 48 h. Cell viability was detected using the cell counting kit (YEASEN, Shanghai, China) according to the instructions.
2.3. Colony Formation Assay
Cells were trypsinized, counted, and seeded in a 6-well plate at 700 cells per well, treated with different concentrations of fenofibrate. After 14 days, the visible colonies were counted and photographed.
2.4. Wound Healing Assay
HCT116 and SW480 cells were grown on 6-well plates. Scratch was made using a 200 μL pipette tip when the cell confluence reached 80%-90%. After scratching, cells were washed with phosphate-buffered saline and then cultured in serum-free medium. The healing rate was quantified with measurements of the gap size after the culture using ImageJ software.
2.5. Cell Immunofluorescence
HCT116 and SW480 cells were seeded in 6-well culture plates plated with cell climbing slices. After being treated with fenofibrate for 24 h, cells were fixed with 4% paraformaldehyde and permeabilized with 0.1% Triton X-100. Then, cells were incubated with primary antibodies (E-cadherin, vimentin, PCNA). The chromosomes were counter-stained with DAPI (Beyotime, Shanghai, China). Images were viewed with a fluorescent microscope.
Apoptosis was analyzed by means of TUNEL assay using the EdUTP TUNEL cell detection kit (Beyotime, Shanghai, China) according to the manufacturer’s instructions. Photomicrographs were taken under confocal microscopy.
2.6. EDU Staining
5-Ethynyl-2-deoxyuridine (EdU) staining assay was carried out on fenofibrate-treated cells utilizing an EdU immunofluorescence staining kit (Ribobio, China) according to the manufacturer’s instructions . The results were observed using an inverted fluorescence microscope (200x).
2.7. Flow Cytometry Analysis
HCT116 and SW480 cell lines were treated with various doses of fenofibrate for 24 h, then stained with annexin V FITC and propidium iodide (PI) (BD Biosciences, San Jose, CA). Cell apoptosis was analyzed by flow cytometry (Cytomics FC500; Beckman Coulter, Fullerton, CA).
The distribution of cell cycle phases was analyzed by the cell cycle detection kit (Beyotime, China). Samples were fixed with 75% ethanol at −20°C for 24 h. The fixed cells were treated with RNaseA and stained with propidium iodide following the manufacturer’s instruction . The cell cycle analysis was studied by flow cytometry.
2.8. Animal Experiments
Animal experiments were performed according to the National Institutes of Health Guidelines for the Care and Use of Laboratory Animals and were approved by the Animal Care and Use Committee of Shanghai Tongji University, China. The nude mice were randomly divided into 2 groups (), and each mouse was subcutaneously injected with HCT116 cells in the right axilla. Fenofibrate was suspended in saline and intragastrically administered at 200 mg/kg per mouse once a day. Mice were anesthetized and sacrificed three weeks after fenofibrate administration. Tumor size and mouse weight were recorded during the experiment. Tissues were harvested for further analysis.
Xenografted tissues were fixed with 4% paraformaldehyde and embedded in paraffin, sectioned, and stained with hematoxylin and eosin. Histopathological changes were observed by microscopy (200x).
2.9. DNMT1 Content
The measurement of DNMT1 content in fenofibrate-treated cells or tumor tissues was conducted using Human or Mouse DNMT1 ELISA kit (HZBIO, Shanghai, China). Samples were added to each well in enzyme-labeled coated plates and incubated at 37°C for 30 minutes. After washing, coloration was developed using chromogen reagents A and B; then, the reaction was terminated by stop solution. The absorbance value was measured at 450 nm wavelength using a microplate reader (BioTek microplate reader).
2.10. Methylation-Specific PCR
Genomic DNA was extracted using a genomic DNA extraction kit (TIANGEN). Eluted DNA (20 μL) was subjected to bisulfite modification using the EpiTect Fast Bisulfite Conversion Kit (Qiagen). Methylation status of CDKN2A promoter was analyzed using the methylation-specific primers (M) and nonmethylation-specific primers (U). The amplification products were separated by agarose gel electrophoresis and visualized by SYBR Green staining under UV light. Primers (Table S1) used in methylation-specific PCR (MSP) were the same as the sequences described by Herman et al. .
2.11. Bioinformatics Analysis of Human Tumor Samples from TCGA Dataset
RNA sequencing profiles and relevant clinical information of pancancer samples were retrieved from The Cancer Genome Atlas (TCGA) data portal in October, 2020. RNA-seq data was normalized by fragments per kilobase per million (FPKM) using log2 scale. Transcription levels of DNMT1, DNMT3a, DNMT3b, PPARs, ACOX1, and CDKN2A were analyzed. Survival probabilities were computed by the Kaplan-Meier method. In addition, Cox regression analysis was conducted to calculate the hazard ratios of genes of interest in multiple cancer types. Moreover, the correlation between clinical stage and gene expression was evaluated via “ggpubr” package in R.
2.12. Western Blot Analysis
Total protein from cells or tissues was extracted using RIPA lysate (Invitrogen, USA). Equal amount of protein samples (40 μg for cell samples per lane and 80 μg for tissue samples per lane) was run on sodium dodecylsulfate-polyacrylamide (SDS) gel and then transferred to polyvinylidene difluoride (PVDF) membranes. Membranes were blocked with 5% nonfat milk for 1 h and then incubated overnight at 4°C with primary antibodies. Anti-DNMT1 (1 : 500, Abcam, Cambridge, UK), anti-PPARA (1 : 1000, Abcam, Cambridge, UK), anti-E2F1 (1 : 500, Cell Signaling Technology, Danvers, MA), anti-pRb (1 : 500, Cell Signaling Technology, Danvers, MA), anti-CDKN2A (1 : 500, Cell Signaling Technology, Danvers, MA), anti-CyclinD1 (1 : 1000, Cell Signaling Technology, Danvers, MA), anti-CDK4 (1 : 1000, Cell Signaling Technology, Danvers, MA), anti-CDK6 (1 : 1000, Cell Signaling Technology, Danvers, MA), and anti-RB (1 : 200, Santa Cruz Biotechnology, CA) antibodies were used in this experiment. Then, the PVDF membranes were incubated with corresponding secondary antibodies at room temperature for 1 h. Primary and secondary antibody diluent for WB was used in this experiment (YEASEN, Shanghai, China). Signals were detected using the Odyssey Two-color Infrared Laser Imaging System (Li-Cor, Lincoln, NE).
2.13. Quantitative Real-Time PCR Analysis
Total RNA from cells or tissues was extracted using the TRIzol reagent according to the manufacturer’s instructions. Approximately 500 ng of extracted RNA was used to synthesize cDNA using the reverse transcription kit (TaKaRa Biotechnology, Dalian, China). Quantitative PCR was carried out in the Applied Biosystems 7500 Real-Time PCR System using 50 ng of cDNA and a SYBR Green PCR master mix (YEASEN, Shanghai, China). Relative gene expression was calculated based on 2△△CT algorithm. All the primers were designed using the principle of span exons to avoid genomic DNA contamination. Primer sequences utilized in this study are provided in Table S2.
2.14. Statistical Analyses
Two-group comparisons were analyzed by Student’s -test; multigroup comparisons were analyzed via one-way ANOVA. Spearman’s correlation analysis was performed to evaluate expression correlation. Kaplan-Meier analysis was carried out to analyze overall survival. Multivariate analysis was performed using the Cox multivariate regression analysis model. values < 0.05 were considered statistically significant. Statistical analysis was performed with SPSS 17.0 software (IBM).
3.1. Analysis of Gene Expression Patterns of PPARA and DNA Methyltransferase in Human Colon Cancer
First, we performed pancancer analyses to identify the mean expression levels of PPARs in different types of tumors. The results indicated that PPARA was lowly expressed in tumorous tissues, as shown in Figure 1(a). Subsequently, the Cox proportional hazards regression analysis was used to evaluate PPARs as prognostic markers in various tumors. Genes with a hazard were significantly correlated with patient outcome. The forest plot in Figure 1(b) shows that in most cancer types, PPARs, including PPARA, may serve as a prognostic indicator of digestive tract cancers. Subsequently, we focused on the clinical significance of PPARA expression. PPARA was downregulated in colon cancer and correlated with TNM stage in the TCGA COAD dataset (Figure 1(c)). Furthermore, following the evaluation of PPARA protein levels, the prognosis of patients with high PPARA or low PPARA expression revealed that higher expression exhibited a better prognosis (Figure 1(d)).
Recent reports have demonstrated that DNA hypermethylation of tumor suppressor genes contributes to cancer progression. In the current study, the mean expression levels of the DNA methyltransferases DNMT1 and DNMT3a and the protein arginine methyltransferase, PRMT6, were evaluated in pancancers using bioinformatics analyses. High levels of several methyltransferases were observed in various tumor tissues (Figure 1(a)). The forest graph indicated the prognostic value of DNMT1 in most cancers, with an (Figure 1(b)). In addition, colon cancer cell lines with high endogenous expression of methyltransferases and low endogenous expression of PPARA relative to normal NCM460 and HIEC cells were selected for further in vitro assays (Figure 1(e)).
3.2. The PPARA Agonist, Fenofibrate, Attenuated Cell Viability and Proliferation
The effect of the PPAR agonist on colon cancer cell proliferation was assessed by treating 5000 cells with various concentrations of fenofibrate (0–300 μmol/L) for 24 and 48 h. Cell survival was determined using the CCK8 kit. As shown in Figure 2(a), fenofibrate exerted an inhibitory effect on cell proliferation in HCT116 cells with an IC50 range of 180–200 μmol/L (24 h) and 100–120 μmol/L (48 h); fenofibrate also exerted an inhibitory effect on cell proliferation in SW480 cells, with an IC50 range of 160–180 μmol/L (24 h) and 80–100 μmol/L (48 h).
We selected suitable concentrations of fenofibrate (180 and 200 μmol/L for HCT116 cells and 160 and 180 μmol/L for SW480 cells) for subsequent treatments. After incubation with fenofibrate for 24 h, cells were stained with EDU (green) and DAPI (blue). Treatment of fenofibrate decreased the proportion of cells with green fluorescence in a dose-dependent manner () (Figure 2(b)). To further demonstrate the inhibitory effect of fenofibrate on cell growth, the colony formation assay was employed. Treatment of HCT116 and SW480 cells with different doses of fenofibrate reduced clone numbers, especially in the high-dose group (Figure 2(c)). The expression level of PCNA, an endogenous marker of mitogenesis, was detected between the low-dose, high-dose, and vehicle-treated groups using qRT-PCR and immunofluorescence staining. Reduced PCNA mRNA level and protein level upon fenofibrate treatment suggested the suppressive regulation of cell proliferation by fenofibrate (Figures 2(d) and 2(e)).
3.3. The PPAR Agonist, Fenofibrate, Promoted Cell Apoptosis In Vitro
To investigate whether PPARA activation promoted cell apoptosis in vitro, colon cancer cells with or without fenofibrate treatment were assessed by flow cytometry. The data suggested that fenofibrate increased cell apoptosis and, in particular, late apoptosis rates (Figure 3(a)). Subsequently, we confirmed the proapoptotic effects of fenofibrate using TUNEL staining, which detects the DNA breaks in apoptotic cells. Figure 3(b) shows that fenofibrate administration increased the number of positively stained cells (green fluorescence). Moreover, the changes in cell shape after incubation with fenofibrate was examined using light microscopy. The number of shedding cells was increased, accompanied with morphological deformation and crumpled appearances. Such morphological changes were observed in a time- and dose-dependent manner (Figure 3(c)). The expression level of Bax, a typical apoptosis-related protein which promotes apoptosis was increased by treatment with fenofibrate, while prosurvival protein Bcl-2 was decreased as shown by western blotting (Figure 3(d)).
3.4. The PPAR Agonist, Fenofibrate, Inhibited Cell Migration and Epithelial Mesenchymal Transition
Cell migration was measured using the wound healing assay in fenofibrate-treated cells. Data from the scratch healing assay are shown in Figure 4(a). Fenofibrate significantly reduced the migratory capacity of CRC cells, which exhibited delays in the closure of scratches. Next, we examined whether there was a repressive effect of PPARA activation on epithelial mesenchymal transition (EMT). The expression of several EMT biomarkers of EMT was tested using qRT-PCR (Figure 4(b)). The mRNA levels of vimentin and MMP9 in the fenofibrate-treated group were reduced, while E-cadherin levels were upregulated. To further confirm the findings, expression changes in the EMT-associated markers induced by fenofibrate were examined using immunofluorescence, and similar results were observed (Figure 4(c)). Taken together, PPARA activation could exert an antitumor effect by restraining cell migration and EMT.
3.5. The PPAR Agonist, Fenofibrate, Enhanced Tumor Suppressor Gene Expression and Repressed DNMT1 Content
The research of Luo et al.  substantiated that the loss of PPARA resulted in abnormal expression of several methyltransferases and promoted CRC progression in a mouse model. Therefore, we hypothesized that the PPARA agonist, fenofibrate, could reduce the content of DNMT1 and rescue the expression of methylation-silenced tumor suppressor genes. First, we found using qRT-PCR that the DNMT1 mRNA level decreased when PPARA was activated following fenofibrate treatment (Figure 5(a)).
Next, the expression of several tumor suppressor genes was determined (Figure 5(b)). It has been documented that the expression of silenced genes is correlated with promoter hypermethylation. Studies showed that p21, p27, CDKN2A, MLH1, and RASSF1A were reactivated in two cancer cell lines after fenofibrate treatment. We further discovered that the expression of the upstream regulators of DNMT1, Oct4, Nanog, and Sox9 was decreased in fenofibrate-treated cells and were correlated with stemness. In addition, we measured the DNMT1 enzyme content in colon cancer cells and found a significant reduction in enzyme concentration upon high-dose treatment of fenofibrate (Figure 5(c)).
3.6. Fenofibrate Recovered the Expression of CDKN2A via Downregulation of DNMT1
To determine the mechanisms of fenofibrate treatment on the demethylation of the CDKN2A promoter, MSP was carried out to evaluate the methylation status of the CDKN2A promoter. In the Figure 6(a), M and U referred to the PCR products of methylated and unmethylated alleles, respectively. Analyses showed that in the fenofibrate-free group, a methylation product was observed, while no band was observed in the fenofibrate-treated group. Thus, we inferred that the methylation status of the promoter was abolished by fenofibrate treatment (Figure 6(a)).
5-Azacytidine is an effective inhibitor of DNMT1. To investigate the effects of the overexpression or downexpression of DNMT1 on CDKN2A, colon cancer cells were treated with 5-azacytidine, DNMT1 overexpression plasmid, and fenofibrate. The results showed that DNMT1 downexpression resulted in increased CDKN2A mRNA and protein levels, while DNMT1 overexpression caused a decreased level of CDKN2A (Figures 6(b) and 6(c)). The expression changes of CDKN2A induced by 5-azacytidine could be enhanced by fenofibrate treatment. In addition, fenofibrate could reverse the low expression of CDKN2A caused by DNMT1 overexpression plasmid. Taken together, these findings indicated that the PPARA agonist, fenofibrate, could upregulate the expression of CDKN2A by inhibiting gene hypermethylation mediated by DNMT1.
3.7. Fenofibrate Regulates the Cell Cycle via the CDKN2A/RB1/E2F1 Pathway
DNMT1 is responsible for maintaining DNA methylation after each round of the cell cycle. CDKN2A, a CDK inhibitor, acts as a negative regulator of cell cycle process. In the current study, fenofibrate administration suppressed the expression of DNMT1 and relieved the DNMT1-mediated silencing of CDKN2A. Therefore, we investigated the molecular mechanism by which fenofibrate regulated cell cycle distribution. PI staining was conducted on control and drug-treated cells, and the results of cell cycle analysis are shown in Figure 7(a). Fenofibrate treatment suppressed the G1→S transition, induced G0/G1 phase cell arrest, and blocked S phase entry. Since the activity of the cyclin D/CDK4/CDK6 complex is essential to the G1/S transition and can be inhibited by CDK inhibitors, we detected the changes in expression of cyclinD1, CDK4, and CDK6 in cells following fenofibrate treatment. The results from western blot and qRT-PCR analyses suggested that fenofibrate induced the activation of CDKN2A, resulting in low expression levels of CDKs (Figures 7(b) and 7(c)). During the G1 phase, the level of E2F1 was upregulated by the activated cyclin-CDK complex and was released from the RB/E2F contigs in response to abnormal growth stimulation.
To further validate our findings, the expression of key factors involved in the RB/E2F signaling pathway was examined. The expression of total RB was upregulated, while the levels of phosphorylated RB and E2F1 were reduced following fenofibrate treatment. These results demonstrated that fenofibrate treatment decreased the protein and transcript levels of pRB and E2F through the downregulation of cyclin-CDKs by activated CDKN2A.
3.8. The PPAR Agonist, Fenofibrate, Inhibited Tumor Growth in an Animal Model
We established an animal model of HCT116 cell-bearing nude mice. When tumor tissues reached an average volume of 1 mm3, mice were randomized into two groups. Mice in the treatment group received 200 mg/kg fenofibrate suspended in 200 μL saline by gavage every day. The control group was gavaged with an equal volume of saline. All mice bearing tumors survived during the experiment. Tumor size was significantly smaller in the drug-treated group compared to the cancer group.
The weight of each mouse was recorded twice weekly, and the data are shown in Figure 8(a). Tumor xenografts were sectioned and stained with H&E to observe the pathological changes. Slides from fenofibrate-treated tissues revealed necrotic lesions (Figure 8(b)). DMNT1 content in tissue homogenates was measured by an ELISA assay. Following fenofibrate treatment, DNMT1 enzyme level in tissues from the treatment group was lower than that from the vehicle group (Figure 8(c)). Additionally, fenofibrate significantly downregulated the mRNA and protein levels of DNMT1 and CDK4 and increased the expression levels of PPARA and CDKN2A (Figures 8(d) and 8(e)). Taken together, PPARA activation suppressed tumors in vivo by upregulating the expression of CDKN2A.
The role of PPARA in tumor initiation and development remains controversial. Several studies have elucidated the antitumor effect of PPARA in various cancers, including breast cancer, prostate cancer, and ovarian cancer [32–34]. Some studies have drawn contradictory conclusions, suggesting that prolonged administration of PPARA agonists might cause hepatocarcinogenesis; however, the detailed mechanism remains unclear . In the present study, we analyzed the expression levels of PPARs in pancancers and noted that PPARA was expressed at low levels in several types of tumors, including colon cancer. Colon cancer samples with high PPARA protein expression were observed to have a better prognosis than those with low PPARA levels. These results suggested that PPARA might serve as a tumor suppressor gene in colon cancer. Most recently, PPARA-specific agonists were reported to exhibit anticancer effects in a variety of tumors. Fenofibrate is a PPARA activator that belongs to the fibrate class of drugs. An increasing number of studies have revealed its potential role as an antitumor agent that affects multiple biological pathways [36–39].
DNMT1 is responsible for maintaining global methylation and aberrant CGI methylation in human cancer cells, whereas DNMT3a and DNMT3b are believed to act as maintenance and de novo methyltransferases. The elevated expression of DNMT1 has been reported in colon adenocarcinomas, hepatocarcinomas, and lung cancer [40–42].
Hypermethylation of gene promoter regions leads to transcriptional repression. Tse et al.  showed that promoter methylation of tumor suppressor genes promoted carcinogenesis of colon cancer. In this study, we observed that fenofibrate treatment increased PPARA expression and decreased DNMT1 activity, accompanied with the elevated expression of a series of established tumor suppressor genes, including RASSF1A, MLH1, p21, and p27. CDKN2A mRNA and protein levels were upregulated in both HCT116 and SW480 fenofibrate-treated cells compared to controls. To confirm the hypothesis that fenofibrate abrogated the hypermethylation of CDKN2A, we detected the methylation status of its promoter using methylation-specific PCR. In reactions using methylation-specific primers, no band for methylated CDKN2A was observed in the fenofibrate-treated group. Furthermore, the expression levels of DNMT1 and CDKN2A were measured in colon cancer cells following treatment with the DNMT1 inhibitor, 5-azacytidine, and DNMT1 overexpression plasmid. The results indicated that fenofibrate functions as a repressor, similar to a methyltransferase inhibitor.
There are several reports suggesting that PPARA activation inhibits cell proliferation by targeting the cyclin-dependent kinase inhibitor, CDKN2A [43, 44]. In the present study, we found that cells treated with fenofibrate were arrested in the G1 phase, and the number of G2/M cells was reduced. Moreover, fenofibrate elevated the level of CDKN2A by suppressing DNMT1 expression, reduced the activation of cyclin D1/CDK complexes, and phosphorylated RB. Cyclin D1/CDK complexes are cell cycle-related molecules that facilitate the G1/S transition . In quiescent cells, hypophosphorylated RB protein interacts with E2F and inhibits its transcription activity. Under growth stimulation or cancerous states, cyclin-CDK complexes are activated to induce RB phosphorylation. Phosphorylated RB then releases the E2F transcription factor, which triggers the transition of the cell cycle from the G1 phase to the S phase, thereby, enabling uncontrolled cell proliferation . In summary, the RB/E2F pathway was involved in fenofibrate-mediated epigenetic changes on CDKN2A, which resulted in the alterations in cell cycle distribution.
Some studies have shown that DNMT1 cooperated directly with E2F1 and HDAC to accelerate aberrant methylation in tumors . The free E2F1 is released by phosphorylated RB and binds to its cognate sites on the DNMT1 promoter region, which played a positive role on DNMT1 expression in cell cycle process . Elevated expression of DNMT1 induced DNA hypermethylation of several tumor suppressor genes, including CDKN2A. However, conflicting results have emerged regarding the relationship between DNMT1 and E2F1. In mesenchymal stem cells, the expression of E2F1 was not correlated with that of DNMT1. Complete cell cycle arrest by serum starvation did not affect the expression of DNMT1, while E2F1 expression was decreased . The above findings revealed that DNMT1 may not be a responsive target of E2F1 during cell cycle arrest.
In the present study, we found that fenofibrate may act in a similar manner as a methylation transferase inhibitor. It reduced DNMT1 activity and E2F1 expression. The mechanistic details behind PPAR agonist on DNMT1 inhibition have not been determined. Whether DNMT1 downregulation caused decreases in free E2F1 requires further investigation.
In addition, we demonstrated that fenofibrate inhibited tumor progression by regulating cell apoptosis and migration. The results from flow cytometry analysis and TUNEL assays showed that fenofibrate caused an increase in late apoptosis in a dose-dependent manner. However, the underlying mechanism for such findings was not determined in our study.
Cellular plasticity mediated by EMT regulatory circuits enhances the invasive properties of cancer cells . Transcriptional repression of E-cadherin is frequently observed in malignant tumor cells. Some studies have verified that DNMT1 caused the suppression of E-cadherin through hypermethylation of its promoter region . We found that E-cadherin expression was increased following fenofibrate treatment in vitro; however, whether the upregulation was correlated with reduced DNMT1 activity and promoter demethylation requires further investigation.
Finally, we carried out tumor xenograft experiments using HCT116 cells to investigate the antitumor efficacy of fenofibrate in vivo. Fenofibrate decreased the tumor volume significantly compared to the vehicle-treated mice. The necrotic area was identified in H&E-stained samples from fenofibrate-treated mice. The expression of DNMT1, CDK4, and CDKN2A was effectively reduced by fenofibrate treatment compared with control cells. These results demonstrated that fenofibrate could ablate tumors and retard tumor growth.
These findings, coupled with the reversibility of DNA methylation, support the possibility of fenofibrate as a potential epigenetic treatment in colon cancer patients.
In conclusion, the present work illustrated that activation of PPARA by fenofibrate administration protected against colon cancer progression through epigenetic modifications. Fenofibrate weakened DNMT1 activity and restored the expression of the tumor suppressor gene, CDKN2A, which suppressed cell proliferation by blocking the G1 to S transition through the RB/E2F pathway. In addition, fenofibrate inhibited cancer cell invasion by regulating EMT. Therefore, we conclude that fenofibrate could act as an adjuvant agent in colon cancer treatment.
|BLCA:||Bladder urothelial carcinoma|
|BRCA:||Breast invasive carcinoma|
|CESC:||Cervical squamous cell carcinoma and endocervical adenocarcinoma|
|DLBC:||Lymphoid neoplasm diffuse large B-cell lymphoma|
|HNSC:||Head and neck squamous cell carcinoma|
|KIRC:||Kidney renal clear cell carcinoma|
|KIRP:||Kidney renal papillary cell carcinoma|
|AML:||Acute myeloid leukemia|
|LGG:||Brain lower grade glioma|
|LIHC:||Liver hepatocellular carcinoma|
|LUSC:||Lung squamous cell carcinoma|
|OV:||Ovarian serous cystadenocarcinoma|
|PCPG:||Pheochromocytoma and paraganglioma|
|SKCM:||Skin cutaneous melanoma|
|TGCT:||Testicular germ cell tumors|
|UCEC:||Uterine corpus endometrial carcinoma|
|FPKM:||Fragments per kilobase per million|
|TCGA:||The Cancer Genome Atlas|
|PPARA:||Peroxisome proliferator-activated receptor alpha|
|DNMT1:||DNA methyltransferase 1|
|EMT:||Epithelial mesenchymal transition|
|qRT-PCR:||Quantitative reverse transcription PCR|
TCGA gene expression data and clinical data were obtained from TCGA data portal (https://portal.gdc.cancer.gov).
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Rui Kong provided the main contribution to this work.
This study received support from National Natural Science Foundation of China (81300340, 81471537).
Table S1: list of primers used in the qRT-PCR assay. Table S2: list of primers used in the methylation-specific PCR assay. (Supplementary Materials)
- F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
- M. Arnold, M. S. Sierra, M. Laversanne, I. Soerjomataram, A. Jemal, and F. Bray, “Global patterns and trends in colorectal cancer incidence and mortality,” Gut, vol. 66, no. 4, pp. 683–691, 2017.
- E. Dekker, P. J. Tanis, J. L. A. Vleugels, P. M. Kasi, and M. B. Wallace, “Colorectal cancer,” Lancet, vol. 394, no. 10207, pp. 1467–1480, 2019.
- J. W. T. Tse, L. J. Jenkins, F. Chionh, and J. M. Mariadason, “Aberrant DNA methylation in colorectal cancer: what should we target?” Trends in Cancer, vol. 3, no. 10, pp. 698–712, 2017.
- V. V. Lao and W. M. Grady, “Epigenetics and colorectal cancer,” Nature Reviews Gastroenterology & Hepatology, vol. 8, no. 12, pp. 686–700, 2011.
- C. Ili, K. Buchegger, H. Demond et al., “Landscape of genome-wide DNA methylation of colorectal cancer metastasis,” Cancers, vol. 12, no. 9, p. 2710, 2020.
- P. Liu, L. Wu, H. Chand, C. Li, X. Hu, and Y. Li, “Silencing of miR-152 contributes to DNMT1-mediated CpG methylation of the PTEN promoter in bladder cancer,” Life Sciences, vol. 261, p. 118311, 2020.
- M. Mamo, I. C. Ye, J. W. DiGiacomo, J. Y. Park, B. Downs, and D. M. Gilkes, “Hypoxia alters the response to anti-EGFR therapy by regulating EGFR expression and downstream signaling in a DNA methylation-specific and HIF-dependent manner,” Cancer Research, vol. 80, no. 22, pp. 4998–5010, 2020.
- M. Kulis and M. Esteller, “DNA methylation and cancer,” Advances in Genetics, vol. 70, pp. 27–56, 2010.
- M. Laisné, N. Gupta, O. Kirsh, S. Pradhan, and P.-A. Defossez, “Mechanisms of DNA methyltransferase recruitment in mammals,” Genes, vol. 9, no. 12, p. 617, 2018.
- L. Gao, M. Emperle, Y. Guo et al., “Comprehensive structure-function characterization of DNMT3B and DNMT3A reveals distinctive de novo DNA methylation mechanisms,” Nature Communications, vol. 11, no. 1, article 3355, 2020.
- T. Hinoue, D. J. Weisenberger, C. P. E. Lange et al., “Genome-scale analysis of aberrant DNA methylation in colorectal cancer,” Genome Research, vol. 22, no. 2, pp. 271–282, 2012.
- M. S. Fernandes, F. Carneiro, C. Oliveira, and R. Seruca, “Colorectal cancer and RASSF family–a special emphasis on RASSF1A,” International Journal of Cancer, vol. 132, no. 2, pp. 251–258, 2013.
- D. S. Byun, N. Ahmed, S. Nasser et al., “Intestinal epithelial-specific PTEN inactivation results in tumor formation,” American Journal of Physiology-Gastrointestinal and Liver Physiology, vol. 301, no. 5, pp. G856–G864, 2011.
- K. Huang, M. Du, X. Tan et al., “PARP1-mediated PPARα poly(ADP-ribosyl)ation suppresses fatty acid oxidation in non-alcoholic fatty liver disease,” Journal of Hepatology, vol. 66, no. 5, pp. 962–977, 2017.
- L. J. Holm, M. Ø. Monsted, M. Haupt-Jorgensen, and K. Buschard, “PPARs and the development of type 1 diabetes,” PPAR Research, vol. 2020, Article ID 6198628, 11 pages, 2020.
- M. Pawlak, P. Lefebvre, and B. Staels, “Molecular mechanism of PPARα action and its impact on lipid metabolism, inflammation and fibrosis in non-alcoholic fatty liver disease,” Journal of Hepatology, vol. 62, no. 3, pp. 720–733, 2015.
- N. Bougarne, B. Weyers, S. J. Desmet et al., “Molecular actions of PPARα in lipid metabolism and inflammation,” Endocrine Reviews, vol. 39, no. 5, pp. 760–802, 2018.
- L. Wu, J. Li, J. Feng et al., “Crosstalk between PPARs and gut microbiota in NAFLD,” Biomedicine & Pharmacotherapy, vol. 136, article 111255, 2021.
- T. Shigeto, Y. Yokoyama, B. Xin, and H. Mizunuma, “Peroxisome proliferator-activated receptor alpha and gamma ligands inhibit the growth of human ovarian cancer,” Oncology Reports, vol. 18, no. 4, pp. 833–840, 2007.
- Y. Gao, X. Chen, Q. He et al., “Adipocytes promote breast tumorigenesis through TAZ-dependent secretion of resistin,” Proceedings of the National Academy of Sciences of the United States of America, vol. 117, no. 52, pp. 33295–33304, 2020.
- S. Shi, G. Yu, B. Huang, Y. Mi, Y. Kang, and J. P. Simon, “PPARG could work as a valid therapeutic strategy for the treatment of lung squamous cell carcinoma,” PPAR Research, vol. 2020, Article ID 2510951, 9 pages, 2020.
- Y. Xu, X. Li, Y. Han et al., “A new prognostic risk model based on PPAR pathway-related genes in kidney renal clear cell carcinoma,” PPAR Research, vol. 2020, Article ID 6937475, 13 pages, 2020.
- S. Du, N. Wagner, and K.-D. Wagner, “The emerging role of PPAR beta/delta in tumor angiogenesis,” PPAR Research, vol. 2020, Article ID 3608315, 16 pages, 2020.
- T. Li, Q. Zhang, J. Zhang et al., “Fenofibrate induces apoptosis of triple-negative breast cancer cells via activation of NF-κB pathway,” BMC Cancer, vol. 14, no. 1, p. 96, 2014.
- N. W. Chang, M. H. Tsai, C. Lin et al., “Fenofibrate exhibits a high potential to suppress the formation of squamous cell carcinoma in an oral-specific 4-nitroquinoline 1-oxide/arecoline mouse model,” Biochimica et Biophysica Acta, vol. 1812, no. 4, pp. 558–564, 2011.
- G. Augimeri, L. Gelsomino, P. Plastina et al., “Natural and synthetic PPARγ ligands in tumor microenvironment: a new potential strategy against breast cancer,” International Journal of Molecular Sciences, vol. 21, no. 24, article 9721, 2020.
- Y. Luo, C. Xie, C. N. Brocker et al., “Intestinal PPARα Protects Against Colon Carcinogenesis via Regulation of Methyltransferases DNMT1 and PRMT6,” Gastroenterology, vol. 157, no. 3, pp. 744–759.e4, 2019.
- Q. Xu, Y. Xiang, Q. Wang et al., “SETD2 regulates the maternal epigenome, genomic imprinting and embryonic development,” Nature Genetics, vol. 51, no. 5, pp. 844–856, 2019.
- H. R. Liu, L. Y. Meng, Z. Y. Lin, Y. Shen, Y. Q. Yu, and Y. Z. Zhu, “Cochinchina momordica seed extract induces apoptosis and cell cycle arrest in human gastric cancer cells via PARP and p53 signal pathways,” Nutrition and Cancer, vol. 64, no. 7, pp. 1070–1077, 2012.
- J. G. Herman, J. R. Graff, S. Myohanen, B. D. Nelkin, and S. B. Baylin, “Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 18, pp. 9821–9826, 1996.
- Y. Segawa, R. Yoshimura, T. Hase et al., “Expression of peroxisome proliferator-activated receptor (PPAR) in human prostate cancer,” Prostate, vol. 51, no. 2, pp. 108–116, 2002.
- S. C. Kwong, A. H. A. Jamil, A. Rhodes, N. A. Taib, and I. Chung, “Metabolic role of fatty acid binding protein 7 in mediating triple-negative breast cancer cell death via PPAR-α signaling,” Journal of Lipid Research, vol. 60, no. 11, pp. 1807–1817, 2019.
- Y. Yokoyama, B. Xin, T. Shigeto et al., “Clofibric acid, a peroxisome proliferator-activated receptor α ligand, inhibits growth of human ovarian cancer,” Molecular Cancer Therapeutics, vol. 6, no. 4, pp. 1379–1386, 2007.
- S. Feng, H. Wang, Y. Wang et al., “Apatinib induces 3-hydroxybutyric acid production in the liver of mice by peroxisome proliferator-activated receptor α activation to aid its antitumor effect,” Cancer Science, vol. 110, no. 10, pp. 3328–3339, 2019.
- D. Panigrahy, A. Kaipainen, S. Huang et al., “PPARalpha agonist fenofibrate suppresses tumor growth through direct and indirect angiogenesis inhibition,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 3, pp. 985–990, 2008.
- D. Yamasaki, N. Kawabe, H. Nakamura et al., “Fenofibrate suppresses growth of the human hepatocellular carcinoma cell via PPARα-independent mechanisms,” European Journal of Cell Biology, vol. 90, no. 8, pp. 657–664, 2011.
- E. Binello, E. Mormone, L. Emdad, H. Kothari, and I. M. Germano, “Characterization of fenofibrate-mediated anti-proliferative pro-apoptotic effects on high-grade gliomas and anti-invasive effects on glioma stem cells,” Journal of Neuro-Oncology, vol. 117, no. 2, pp. 225–234, 2014.
- T. Morinishi, Y. Tokuhara, H. Ohsaki, E. Ibuki, K. Kadota, and E. Hirakawa, “Activation and expression of peroxisome proliferator-activated receptor alpha are associated with tumorigenesis in colorectal carcinoma,” PPAR Research, vol. 2019, Article ID 7486727, 9 pages, 2019.
- L. Chen, J. Peng, Y. Wang et al., “Fenofibrate-induced mitochondrial dysfunction and metabolic reprogramming reversal: the anti-tumor effects in gastric carcinoma cells mediated by the PPAR pathway,” American Journal of Translational Research, vol. 12, no. 2, pp. 428–446, 2020.
- Y. Bai, L. Lang, W. Zhao, and R. Niu, “Long non-coding RNA HOXA11-AS promotes non-small cell lung cancer tumorigenesis through microRNA-148a-3p/DNMT1 regulatory axis,” Oncotargets and Therapy, vol. 12, pp. 11195–11206, 2019.
- J. M. Serrano-Morales, M. D. Vázquez-Carretero, M. J. Peral, A. A. Ilundáin, and P. García-Miranda, “Reelin-Dab1 signaling system in human colorectal cancer,” Molecular Carcinogenesis, vol. 56, no. 2, pp. 712–721, 2017.
- T. L. Kruer, S. M. Dougherty, L. Reynolds et al., “Expression of the lncRNA maternally expressed gene 3 (MEG3) contributes to the control of lung cancer cell proliferation by the Rb pathway,” PLoS One, vol. 11, no. 11, article e0166363, 2016.
- F. Gizard, T. Nomiyama, Y. Zhao et al., “The PPARalpha/p16INK4a pathway inhibits vascular smooth muscle cell proliferation by repressing cell cycle-dependent telomerase activation,” Circulation Research, vol. 103, no. 10, pp. 1155–1163, 2008.
- Y. Majeed, R. Upadhyay, S. Alhousseiny et al., “Potent and PPARα-independent anti-proliferative action of the hypolipidemic drug fenofibrate in VEGF-dependent angiosarcomas in vitro,” Scientific Reports, vol. 9, no. 1, article 6316, 2019.
- A. N. Srinivas, D. Suresh, F. Mirshahi, P. K. Santhekadur, A. J. Sanyal, and D. P. Kumar, “Emerging roles of AATF: checkpoint signaling and beyond,” Journal of Cellular Physiology, vol. 236, no. 5, pp. 3383–3395, 2021.
- K. D. Robertson, S. Ait-Si-Ali, T. Yokochi, P. A. Wade, P. L. Jones, and A. P. Wolffe, “DNMT1 forms a complex with Rb, E2F1 and HDAC1 and represses transcription from E2F-responsive promoters,” Nature Genetics, vol. 25, no. 3, pp. 338–342, 2000.
- F. Leng, J. Yu, C. Zhang et al., “Methylated DNMT1 and E2F1 are targeted for proteolysis by L3MBTL3 and CRL4DCAF5 ubiquitin ligase,” Nature Communications, vol. 9, no. 1, article 1641, 2018.
- C. C. Tsai, P. F. Su, Y. F. Huang, T. L. Yew, and S. C. Hung, “Oct4 and Nanog directly regulate Dnmt1 to maintain self-renewal and undifferentiated state in mesenchymal stem cells,” Molecular Cell, vol. 47, no. 2, pp. 169–182, 2012.
- M. Saitoh, “Involvement of partial EMT in cancer progression,” Journal of Biochemistry, vol. 164, no. 4, pp. 257–264, 2018.
- A. Fukagawa, H. Ishii, K. Miyazawa, and M. Saitoh, “δEF1 associates with DNMT1 and maintains DNA methylation of the E-cadherin promoter in breast cancer cells,” Cancer Medicine, vol. 4, no. 1, pp. 125–135, 2015.
Copyright © 2021 Rui Kong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.