Chun-Qing Li
Zhi-Qin Liu
Sha-Sha Liu
Gao-Tao Zhang
Li Jiang1
Du-Qiang Luo- 1Key Laboratory of Microbial Diversity Research and Application of Hebei Province, College of Life Science, Hebei University, Baoding, China
- 2Key Laboratory of Pharmaceutical Quality Control of Hebei Province, College of Pharmaceutical Science, Hebei University, Baoding, China
- 3College of Science and Technology, Hebei Agricultural University, Huanghua, China
Metformin is a kind of widely used antidiabetic drug that regulates glucose homeostasis by inhibiting liver glucose production and increasing muscle glucose uptake. Recently, some studies showed that metformin exhibits anticancer properties in a variety of cancers. Although several antitumor mechanisms have been proposed for metformin action, its mode of action in human liver cancer remains not elucidated. In our study, we investigated the underlying molecular mechanisms of metformin's antitumor effect on Huh-7 cells of hepatocellular carcinoma (HCC) in vitro. RNA sequencing was performed to explore the effect of metformin on the transcriptome of Huh-7 cells. The results revealed that 4,518 genes (with log2 fold change > 1 or < −1, adjusted p-value < 0.05) were differentially expressed in Huh-7 cells with treatment of 25-mM metformin compared with 0-mM metformin, including 1,812 upregulated and 2,706 downregulated genes. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses identified 54 classical pathways that were significantly enriched, and 16 pathways are closely associated with cancer, such as cell cycle, DNA replication, extracellular matrix–receptor interaction, and so on. We selected 11 differentially expressed genes, which are closely associated with HCC, to validate their differential expressions through a quantitative real-time reverse transcription-polymerase chain reaction. The result exhibited that the genes of fatty acid synthase, mini-chromosome maintenance complex components 6 and 5, myristoylated alanine-rich C-kinase substrate, fatty acid desaturase 2, C-X-C motif chemokine ligand 1, bone morphogenetic protein 4, S-phase kinase-associated protein 2, kininogen 1, and proliferating cell nuclear antigen were downregulated, and Dual-specificity phosphatase-1 is significantly upregulated in Huh-7 cells with treatment of 25-mM metformin. These differentially expressed genes and pathways might play a crucial part in the antitumor effect of metformin and might be potential targets of metformin treating HCC. Further investigations are required to evaluate the metformin mechanisms of anticancer action in vivo.
Introduction
Liver cancer is the most common malignant tumor and the third leading cause of cancer-related mortality worldwide, just next to lung and colorectal cancers. It is estimated that there are approximately 850,000 new cases of liver cancer each year, with a mortality rate of 8.2% of all cancers worldwide (Llovet et al., 2016; Bray et al., 2018). Liver cancer occurs in association with chronic liver injury and hepatocellular death caused by viral hepatitis, chronic alcohol abuse, drugs or chemical insults, obesity, and aflatoxin-B1-contaminated food (Farazi and DePinho, 2006; El-Serag and Rudolph, 2007). Hepatocellular carcinoma (HCC) is the most prevalent type of liver cancer, representing over 80% of all liver cancer cases (Yi et al., 2018). As a highly metastatic and recurrent cancer, HCC tends to be limited with the treatment choices of chemotherapy, surgical resection, liver transplantation, and local ablation in its early stages. Sorafenib is a multikinase inhibitor and has a poor outcome despite being a widely used agent for advanced HCC (Llovet et al., 2008). Similarly, low therapeutic effects were reported for cabozantinib, lenvatinib, and regorafenib (Bruix et al., 2017; Abou-Alfa et al., 2018; Kudo et al., 2018). Therefore, it is urgent to develop new effective drugs to target the HCC cells.
Metformin, as an oral hypoglycemic agent, has been used as the first-line treatment of type 2 diabetes mellitus with a safety profile. Recently, it has been reported that metformin can prevent and control the development of a variety of cancers (Evans et al., 2005; Morales and Morris, 2015).
Abundant epidemiological data have indicated that metformin used in patients with type 2 diabetes mellitus can significantly reduce the incidence of HCC (Lee et al., 2011; Nkontchou et al., 2011; Lee et al., 2012). Moreover, other recent epidemiological studies have demonstrated that metformin treatment can reduce tumor recurrence or improve the prognosis of patients who underwent hepatic resection for HCC (Donadon et al., 2010; Li et al., 2009; Lee et al., 2012; Eisenreich and Leppert, 2017). Several laboratory studies have also confirmed that metformin can inhibit the growth of liver cancer cells by inducing cell arrest and enhancing apoptosis (Buzzai et al., 2007; Ben et al., 2011). Although several mechanisms have been proposed for metformin action, the explicit molecular basis of its antitumor activity in HCC remained incompletely understood. In recent years, with the rapid development of next-generation sequencing technologies, RNA sequencing technology has developed into an effective tool to identify functional genes after treatment with various drugs (Sun et al., 2015). Recently, numerous studies were carried out on the transcriptome profiling of tumor cells for functional gene identification. To gain an intensive insight into the molecular mechanism of metformin treatment on HCC cell Huh-7, we performed a transcriptome screening analysis to identify gene targets and signaling pathways through which metformin might have influenced HCC. Meanwhile, this research is expected to provide an experimental basis for the application of metformin in HCC treatment and help to suggest directions for future investigation.
Materials and Methods
Chemicals and Reagents
Dulbecco's modified Eagle's medium (DMEM, 1 g/L D-glucose) was purchased from Gibco. The TRIzol reagent and 0.25% trypsin– ethylenediaminetetraacetic acid (1×) were purchased from Invitrogen (Carlsbad, CA, United States). Fetal bovine serum (900-108) was from GEMINI. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT, M8180) was obtained from Solarbio (Beijing, China). 1,1-Dimethylbiguanide hydrochloride (Metformin, D150959) was purchased from Sigma-Aldrich (Saint Louis, United States). SYBR Green polymerase chain reaction (PCR) Master Mix was purchased from Vazyme (Nanjing, China). TransScript All-in-One First-Strand complementary DNA (cDNA) Synthesis SuperMix was purchased from TRANSGEN BIOTECH (Beijing, China).
Cell Lines and Cell Culture
HCC cell line Huh-7 was purchased from the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China) and cultured in DMEM medium supplemented with 10% fetal bovine serum at 37°C in a culture incubator with 5% CO2. The cell culture medium was refreshed every other day. All cell culture experiments were carried out within 10 passages after being thawed.
3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide Assay
The exponentially growing Huh-7 cells were grown in 96-well plates with 5,000 cells per well with a 100-µl medium for 24 h and then starved for 12 h at 37°C in a humidified atmosphere containing 5% CO2. Subsequently, the Huh-7 cells were treated with the concentrations of 0- and 20-mM metformin dissolved in a complete DMEM. After culturing for 24, 48, and 72 h in a 37°C humidified incubator, 20 µl of 5 mg/ml MTT reagent was added to each well, and the plates were incubated for 4 h. A total of 100 µl of 0.01-M SDS-HCl was then added into each plate and incubated at 37°C overnight, and the optical density value at 570 nm was examined using a microplate reader (Thermo Fisher Scientifc, Inc., Waltham, MA, United States).
Real-Time Cell Analysis of Metformin
Huh-7 cells in the logarithmic growth phase were planted in 16-well plates (E-plate 16 ACEA Biosciences Inc., San Diego, United States) according to the xCELLigence Real-Time Cell Analyzer (RTCA) DP instrument manual as provided by the manufacturer. To be specific, Huh-7cells were harvested by trypsinization using 0.25% trypsin–ethylenediaminetetraacetic acid and then planted into 16-well plates with 5,000 cells per well with a 100-µl medium. Before planting, cell concentration, aggregation rate, and viability were evaluated by automated Trypan blue exclusion using the Cedex XS Analyzer (Roche, Cedex), and background impedance was measured in 50-μl cell culture medium per well. After planting, impedance was detected in the following schedule: 1-min intervals for 2 h, 5-min intervals for 2 h, and 15-min intervals until metformin treatment. After being starved for 12 h, Huh-7 cells seeded in each plate were treated with the concentrations of 0- and 20-mM metformin dissolved in a complete DMEM. The impedance values were examined in 1-min intervals for 2 h, 5-min intervals for 2 h, and 15-min intervals for the remaining time. The cell index values based on impedance were normalized to the time point of metformin treatment and quantified by the RTCA software program.
Cell Cycle Analysis
Huh-7 cells were planted into six-well plates at a concentration of 1.5 × 105 cells per well with a 2-ml medium and maintained at 37°C and 5% CO2 for 24 h in a humid environment. After being starved in serum-free medium for 12 h, Huh-7 cells were treated with 0- and 20-mM concentrations of metformin for 48 h in a 37°C incubator and then harvested, washed by refrigerated phosphate-buffered saline, fixed in 70% refrigerated ethanol for 4 h at 4°C, and centrifuged at 1,500 rpm for 5 min. Subsequently, the fixed cells were washed by refrigerated phosphate-buffered saline, centrifuged at 1,500 rpm for 5 min and incubated with a mixture of 10-µl RNase, 25-µl propidium iodine (PI), and 500-µl dye buffer in the dark for 30 min at 37°C. Then, cell cycle distribution was analyzed by a fluorescence-activated cell sorting flow cytometer (BD Biosciences, NJ, United States) in G0/G1 phase, S phase, and G2/M phase.
Total RNA Extraction
Huh-7 cells were planted in six-well plates at 1.5 × 105 cells per well and treated with 0- and 20-mM metformin for 48 h after being starved. Total RNA was pooled from cells and extracted with TRIzol reagent following the detailed procedure provided by the manufacturer. The purity of RNA was assessed with a NanoPhotometer spectrophotometer (IMPLEN, CA, United States). The concentration of RNA was tested with Qubit 3.0 Fluorometer (Life Technologies, CA, United States). The integrity of RNA was confirmed with Agilent 2100 RNA Nano 6000 Assay Kit (Agilent Technologies, CA, United States). All samples in this study passed quality control with a minimum RNA integrity number of 8.0.
Library Construction and RNA Sequencing
A transcriptome sequencing trial was carried out by ANOROAD (Beijing, China). Six cDNA libraries were constructed, i.e., three for the 20-mM metformin group and three for the control group. After finishing libraries construction, qualified libraries were sequenced on an Illumina Hiseq platform.
Data Analysis
Raw reads were obtained by converting original sequencing images to sequential reads. Clean reads were gained by completing data processing such as filtering low-quality reads, trimming adapter polluted reads, removing reads rich in ploy-N, and so on. The succeeding analyses were on the basis of high-quality clean reads after filtering. The genome version we referred to was Homo_sapiens.GRCh38.87.chr. Cufflinks v2.1.1 was used to analyze differentially expressed genes (DEGs) between the metformin and control groups. Transcript abundance representing the level of gene expression was examined by fragments per kilobase of exon model per million mapped fragments. Cuffdiff software was used to analyze DEGs according to the following three criteria: 1) log2 fold change value > 1 or < −1; 2) q-value false discovery rate < 0.05; and 3) fragments per kilobase of exon model per million mapped fragments value > 1.
Functional Enrichment Analysis
Gene ontology (GO) functional enrichment analysis was carried out to interpret their biological functions to figure out the DEGs with metformin treatment. We used the GOseq R package to map all DEGs and to find significantly enriched GO terms by comparison with the genomic background and defined the GO terms with a corrected p-value of less than 0.05 to significantly enrich in DEGs. We used KOBAS software to perform Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs compared with transcriptome background through a hypergeometric test. A hypergeometric p < 0.05 was considered to be significantly enriched.
Quantitative Real-Time Polymerase Chain Reaction on Gene Expressions
The expression levels of total RNA from Huh-7 cells with metformin treatment for 48 h were detected by quantitative real-time (RT) PCR. Glyceraldehyde-3-phosphate dehydrogenase was used as a control gene. The primer sequences of the target and the control genes are shown in Table 1. Specific primers were selected from Gene Bank and synthesized by GENEWIZ (Suzhou Genomics Institute, Suzhou, China). Total RNA (1 μg) was extracted from treated Huh-7 cells using Trizol and reversely transcribed using the RT reagent kit. The real-time PCR reaction was carried out in a 10-μl buffer mixture composed of 5-μl AceQ qPCR SYBR Green Master Mix, 0.2 μM of each primer, 0.31-μl template cDNA, and 4.29-μl water free of nuclease. Thermal cycling was carried out following the instructions below: 95°C for 5 min, 40 cycles at 95°C for 10 s, and 60°C for 30 s. Quantitative RT-PCR assay was done in a 96-well optical plate (Roche, Light Cycler 96). The relative quantity of RNA transcript was worked out with the formula 2ΔCt (comparative cycling-threshold) that represents the difference between target gene expression in treated and control groups. The fold change of target gene cDNA relative to control gene was calculated: FC = 2−ΔΔCt.
TABLE 1. List of gene primers.
Statistical Analysis
All data of each group were analyzed by mean ± standard deviation for at least three independent trials using SPSS software version 16.0. Unpaired Student's t-test or Mann–Whitney U test was applied to compare the difference between two groups. One-way analysis of variance and Tukey's post-hoc test were applied to compare the multiple groups' differences. The difference has a statistical significance when the p-value is less than 0.05.
Results
Metformin Inhibits the Proliferation of Human Liver Cancer Huh-7 Cells
Previous studies have shown that metformin has an inhibiting effect in various cancer cells (Hsieh et al., 2014; Shi et al., 2017; Yu et al., 2017). To inspect the effect of metformin on the proliferation of Huh-7 cells, an MTT assay was utilized to examine cell viability. The results revealed that metformin inhibited the proliferation of Huh-7 cells in a time-dependent manner treated with 0- and 20-mM metformin for 24, 48, and 72 h (Figure 1A).
FIGURE 1. Effect of metformin (Met) on cell growth in Huh-7 cells. (A) Cell proliferation was tested through MTT assay. Huh-7 Cells were handled with 20-mM Met for 24, 48, and 72 h, and cell proliferation was tested by MTT reduction assay. Data are presented by mean ± standard deviation of five independent experiments performed in triplicates. ***p < 0.001, **p < 0.01, *p < 0.05 versus CON without metformin. (B) Real-time monitoring of proliferation of Huh-7 cells on RTCES platform. Cells suspensions were added into E-Plates and laid on RTCES reader to real-timely monitor every 30 min during assay. Arrow denotes point-in-time of adding metformin.
Time-Dependent Cell Response Analysis of Metformin by Real-Time Cell Analyzer
RTCA is applied to monitor the effect of small molecule compounds and approved drugs on a cell culture based on impedance readout of time-dependent cell response profiles (TCRPs). The utility of this approach is to test how a bioactive compound affects cell proliferation, migration, infiltration, and intercellular interaction. The research paper described that compounds possessed similar activity if they had displayed similar TCRPs based on impedance, and the mechanism of small molecule compounds might be predicted through the TCRP method (Abassi et al., 2009). Our study showed that a bell-shaped TCRP was tested after metformin treatment and the cell index first increased and then decreased compared with the control group, which reflected a final cytotoxic reaction (Figure 1B). In Abassi et al.'s study, bell-shaped TCRP suggested that metformin might interfere with DNA unwinding, replication, synthesis, transcription, and translation, which is in accordance with our transcriptome outcome.
Metformin Inhibits Cell Cycle Progression in Human Liver Cancer Cells
The cell cycle is composed of a range of tightly controlled events that take place in a cell and drive the replication of DNA and cell division. It was reported that the abnormal control of the cell cycle might be closely associated with the occurrence and development of tumors (Kastan and Bartek, 2004; Ingham and Schwartz, 2017). Therefore, to investigate the effect of pharmacological intervention on the cell cycle is developing into a potential pathway for tumor treatment. To further explore the mechanisms of the cell growth inhibition induced by metformin, we carried out a flow cytometry analysis to inspect alterations of the cell cycle after metformin treatment on Huh-7 cells for 48 h. As shown in Figure 2, metformin could have a significant effect on the cell cycle distribution. Compared with the CON group, 20-mM metformin increased the G1 proportion of Huh-7 cells from 45.18 to 61. 48% and decreased S-phase proportion of cell cycle from 43.16 to 22.00%. These results suggested that metformin inhibited Huh-7 cells growth through arresting cells at the G1 phase.
FIGURE 2. Effect of metformin on cell cycle in Huh-7 cells. Cells were treated with 20-mM Met for 48 h and analyzed for propidium iodide-stained DNA content using flow cytometry. Column diagram shows result of statistical analysis. ***p < 0.001, **p < 0.01, *p < 0.05 versus CON without metformin.
RNA Sequencing Transcriptome Analysis of Metformin Treatment
To further investigate the antitumor biological mechanisms involved with Huh-7 cells treated with metformin, RNA sequencing was performed to detect the effect of metformin on the transcriptome with or without treatment of 20-mM metformin for 48 h in triplicate. A total of 37.24-Gb clean bases were obtained, of which 18.53 Gb from the control group (0-mM metformin, three replicates) and 18.71 Gb from the treatment group (20-mM metformin, three replicates). The clean data of all samples reached more than 6.52 Gb, and the percentage of Q30 base was more than 96%. Then, these clean reads of the two groups were mapped to reference genomes version GRCh38 from Homo_sapiens for quantitative measurement of gene expression, and the mapped rate ranged from 98.15 to 98.56%, indicating high levels of gene expression in all samples. A total of 4,518 DEGs were obtained based on screening conditions of adjusted p-value < 0.05 and |log2FoldChange| > 1 with or without metformin treatment in Huh-7 cells. Among those genes, 1,812 were upregulated, and 2,706 were downregulated (Figures 3A–C). Some significant genes are shown in Tables 2, 3. Those genes were related to such signaling pathways as cell cycle, cell communication, cell adhesion, blood vessel morphogenesis, extracellular matrix organization, DNA replication, positive regulation of phosphorylation, and so on.
FIGURE 3. Differentially expressed genes in Huh-7 cells treated with 20-mM metformin for 48 h. (A) Correlation coefficient plot. Right and bottom of plot represent sample names, and left and top of plot represent sample clustering. Higher value of correlation coefficient, greater correlation between two samples. (B) Volcano map. x-Axis indicates log fold change of differential expression genes between treatment and control groups, and y-axis indicates −log 10 of adjusted p-value. Each point represents a gene. Red points present significantly upregulated genes, green points present significantly downregulated genes, and blue dots were genes with not significantly different expression, which met screening condition of FC > 1, adjusted p-value (adj P) < 0.05. (C) Differential genes cluster heat map. It was a hierarchical clustering analysis of up- and downregulated genes. Red and purple, respectively, indicate up- and downregulated genes.
TABLE 2. Twenty significantly upregulated genes linked to cancer in Huh-7 cells treated with metformin.
TABLE 3. Twenty significantly downregulated genes linked to cancer in Huh-7 cells treated with metformin.
Gene Ontology Analysis of Differentially Expressed Genes
To illustrate the functional changes of gene expression associated with metformin treatment in Huh-7 cells, a GO enrichment analysis of 4,518 DEGs from three libraries was carried out. We derived 575 enriched GO terms (q < 0.05), 430 categories in biological processes, 58 categories in cellular components, and 87 categories in molecular functions. The significant enrichment biological processes included cell cycle, cellular response to chemical stimulus, regulation of signaling, biological adhesion, positive regulation of phosphorylation, blood vessel morphogenesis, cell population proliferation, regulation of cell adhesion, and so on (Figures 4A,B and Table 4). Meanwhile, we performed GO enrichment analysis of upregulated and downregulated DEGs. The downregulated genes were significantly enriched in 212 GO terms (q < 0.05), which are 146 categories in biological process, 29 categories in a cellular component, and 37 categories in molecular function. The top of the ranking categories in the biological process included DNA replication, fatty acid (FA) biosynthetic process, steroid metabolic process, mitotic cell cycle, and so on (Figures 4C,D and Table 5). The upregulated genes were significantly enriched in 618 GO terms (q < 0.05), among which 469 in biological process, 44 in a cellular component, and 105 in molecular function. The top of the ranking list in biological process with upregulated genes included regulation of positive regulation of intracellular signal transduction, cell surface receptor signaling pathway, response to cytokine, and so on (Figures 4E,F and Table 6). The enrichment degree of membrane transport and molecular regulation was the highest, and the enrichment rate of protein binding was the highest.
FIGURE 4. GO enrichment analysis of DEGs in metformin-treated Huh-7. (A) Gene ontology (GO) functional classification on all DEGs between control group and metformin treatment group. (B) Bubble chart of all DEG GO analyses. (C) GO functional classification on upregulated DEGs. (D) Bubble chart of upregulated DEG GO analyses. (E) GO functional classification on downregulated DEGs. (F) Bubble chart of downregulated DEG GO analyses. In column charts, upper x-axis indicates quantity of genes corresponding to different points on fold line; lower x-axis indicates negative logarithm of adjusted p-value of GO terms representing significance enrichment level corresponding to height of column; y-axis indicates GO terms. In bubble charts, y-axis represents GO terms; x-axis indicates Rich factor, which is ratio of sample quantity enriched in GO term to background quantity of annotated genes. Bigger Rich factor, higher enrichment degree. Dot dimension represents quantity of genes of GO terms, and dot color corresponds to different p-value ranges.
TABLE 4. List of 20 biological process entries for all differentially expressed genes.
TABLE 5. List of 20 biological process entries for upregulating differentially expressed genes.
TABLE 6. List of 20 biological process entries for downregulating differentially expressed genes.
Genes in an organism mutually coordinate to perform biological functions, and pathway analysis is conducive to further understanding them. To identify molecular mechanisms of metformin acted on Huh-7 cells, the 4,518 DEGs were subjected to KEGG enrichment analysis. As a consequence, 54 KEGG pathways (p < 0.05) (Figures 5A,B and Table 7) were significantly differentially enriched, which included microRNAs in cancer, AMPK signaling pathway, cytokine–cytokine receptor interaction, complement and coagulation cascades, TNF signaling pathway, pathways in cancer, FA biosynthesis, basal cell carcinoma, cGMP-PKG signaling pathway, PPAR signaling pathway, and so on. The 1,812 upregulated genes were significantly enriched in 63 KEGG pathways (p < 0.05) (Figures 5C,D and Table 8), which included MAPK signaling pathway, AMPK signaling pathway, Ras signaling pathway, phosphatidylinositol signaling system, ErbB signaling pathway, inositol phosphate metabolism, TNF signaling pathway, Hippo signaling pathway—multiple species, FoxO signaling pathway, aldosterone synthesis and secretion, Jak-STAT signaling pathway, PI3K-Akt signaling pathway, and so on. The 2,706 downregulated genes were significantly enriched in 41 KEGG pathways (p < 0.05) (Figures 5E,F and Table 9), which included DNA replication, cell cycle, FA biosynthesis, valine, leucine and isoleucine degradation, extracellular matrix (ECM)–receptor interaction, transforming growth factor (TGF)-β signaling pathway, Fanconi anemia pathway, basal cell carcinoma, synthesis and degradation of ketone bodies, and so on. In these pathways, we paid most attention to these pathways, such as DNA replication, cell cycle, and ECM–receptor interaction.
FIGURE 5. KEGG enrichment analysis of DEGs. (A) KEGG enrichment analysis of all DEGs. (B) Bubble chart of all DEGs about KEGG enrichment analysis. (C) KEGG enrichment analysis of upregulated DEGs. (D) Bubble chart of upregulated DEGs about KEGG enrichment analysis. (E) KEGG enrichment analysis of downregulated DEGs. (F) Bubble chart of downregulated DEGs about KEGG enrichment analysis. In column charts, lower x-axis represents number of genes compared with this pathway, corresponding to different dots on fold line; upper x-axis indicates significance enrichment level, corresponding to height of column. Smaller FDR is, bigger −log10 (q value) is, and more significantly enriched KEGG pathway is. In bubble charts, y-axis represents KEGG pathways and their names; x-axis indicates Rich factor, which is ratio of sample quantity enriched in this pathway to background quantity. Bigger Rich factor, higher enrichment degree. Dot dimension represents gene quantity, and dot color corresponds to different p-value ranges.
TABLE 7. List of 20 KEGG analyses of all differentially expressed genes.
TABLE 8. List of 20 KEGG analyses of upregulated differentially expressed genes.
TABLE 9. List of 20 KEGG analyses of downregulated differentially expressed genes.
Quantitative Real-Time Polymerase Chain Reaction Validation
To verify the transcriptome sequencing outcome and to further investigate gene expression pattern associated with cancer, 11 representative genes were selected for quantitative RT-PCR from KEGG data, including microRNAs in cancer pathway [myristoylated alanine-rich C-kinase substrate (MARCKS)], DNA replication pathway [mini-chromosome maintenance complex component 6 (MCM6), MCM5, and proliferating cell nuclear antigen (PCNA)], pathway in cancer [bone morphogenetic protein 4 (BMP4) and S-phase kinase-associated protein 2 (SKP2)], complement and coagulation cascades [kininogen 1 (KNG1)], FA biosynthesis [fatty acid synthase (FASN)], MAPK pathway [dual-specificity phosphatase-1 (DUSP1)], PPAR signaling pathway [fatty acid desaturase 2 (FADS2)], and TNF signaling pathway [C-X-C motif chemokine ligand 1 (CXCL1)]. The RT-PCR results showed that Huh-7 cells treated with metformin compared with the control group had much higher expression of DUSP1 (p < 0.05) but lower expressions of FASN, MCM6, MCM5, MARCKS, FADS2, CXCL1, BMP4, SKP2, KNG1, and PCNA (p < 0.05), which was consistent with the transcriptome sequencing outcome (Figure 6).
FIGURE 6. Quantitative real-time polymerase chain reaction validation of DEGs. Total RNA was isolated from Huh-7 cells treated with 20-mM metformin and carried out a quantitative real-time polymerase chain reaction trial in triplicate. Glyceraldehyde-3-phosphate dehydrogenase was a reference gene. Use a comparative threshold cycle (Ct) method to analyze data. *p < 0.05,**p < 0.01, ***p < 0.001 versus CON without metformin.
Discussion
Human liver cancer cell lines are commonly used for the detection of anticancer drugs and the research on physiological characteristics of liver cells because they are derived from the tissues of hepatocarcinoma patients, have a set of complete genes, and are able to be relatively stably passaged. The widely used liver cell lines contain Huh-7, HepG-2, SMMC-7721, Bel-7402, and so on. Huh-7 cell line possessing the properties of being positive for AFP (Nakabayashi et al., 1985) and negative for HBV and susceptible to HCV (Sainz et al., 2009) compared with other liver cancer lines can be applied to study the regulation mechanism of gene expression, metabolism, xenotransplantation animal models, and so on.
Transcriptome analysis has been used to identify a potential mechanism of drug action on cancer from the perspective of functional genes and pathway networks (Abaan et al., 2013). An investigative transcriptome analysis was carried out to elucidate an underlying action mechanism of metformin in liver cancer cells. We discovered probable gene targets of metformin through researching differential gene expression and described underlying signaling pathways induced by metformin. Although further investigation is needed to verify these findings, this study may still supply a new perspective for functional validation of succeeding trials. The occurrence and development of cancer are complicated processes involving gene changes. The analysis of participation of a variety of signal molecules in the biological process helps to clarify metformin's antitumor mechanism. The transcriptome analysis revealed that loads of genes were significantly altered in Huh-7 cells with treatment of metformin, implying that these significantly changed genes might play a significant role in the antitumor effect of metformin in vitro. These genes involved 10 downregulated ones and 1 upregulated one in Huh-7 cells treated with metformin.
DUSP1 is a kind of phosphatase and can negatively regulate the activation of JNKs, ERKs, and p38 MAPK, which plays an important role in the regulation of the human cell growth cycle and tumor genesis and development through signaling pathways such as the MAPK pathway (Ueda et al., 2003). In this research, the expression of DUSP1 was upregulated in Huh-7 cells after metformin treatment. Previous studies showed that the DUSP1 gene, as a transcriptional target of the p53 tumor suppressor, is highly expressed in prostate cancer (Hou et al., 2012) and played a vital role to inhibit tumor progression in these tumors such as prostate cancer, bladder cancer, and colon cancer (Slattery et al., 2012).
In addition, some genes were significantly downregulated in Huh-7 cells treated with metformin, including FASN, MCM5, MCM6, PCNA, MARCKS, FADS2, CXCL1, BMP4, SKP2, and KNG1.
FASN is a critical enzyme of synthetizing long-chain FAs from malonyl-CoA and acetyl-coenzyme A (CoA) and has been identified at high levels in several human cancers such as prostate cancer (Souchek et al., 2017), pancreatic cancer (Alo et al., 2007), breast cancer (Gonzalez-Guerrico et al., 2016), and colorectal cancer (Zaytseva et al., 2015). FASN overexpression in cancers had frequently been related to biological aggressiveness and a poor prognosis (Gansler et al., 1997). Some studies revealed that FASN overexpression was closely associated with the PI3K-AKT pathway and SREBP 1c transcriptional regulation. Therefore, FASN was considered one of the more than 600 most promising therapeutic drug targets in human cancers (Behan et al., 2019). Proteins of the mini-chromosome maintenance family, including MCM2–7, are collectively referred to as MCM, which mainly affect the stability of microchromosome mitosis. MCM5 and MCM6, as two members of the MCM family, are associated with the initiation and extension of DNA replication and cell proliferation (Semple and Duncker, 2004). MCM5 overexpression has been discovered in some human cancers containing cervical cancer, bladder cancer, and so on (Alison et al., 2002). Therefore, MCM proteins could always be identified as a biomarker for tumor cells that might improve diagnostic accuracy. PCNA is a kind of nuclear protein that ubiquitously exists in normal proliferating cells and tumor cells. PCNA synthesized during all stages of the cell cycle has a crucial effect on replicating and repairing DNA through supplying DNA polymerases, and loads of protein partners participated in some significant biological processes such as DNA repair and cell cycle (Naryzhny, 2008). PCNA is an indispensable constituent of the ultimate mutual pathway shared by all mitogenic signals and is always in a central position for numerous other cell-signaling pathways. PCNA also has a high level of expression in numerous tumor cells. Therefore, PCNA is identified as a potentially valuable biomarker for cancer diagnosis and prognostic therapeutics (Naryzhny and Lee, 2007). MARCKS, a non-helical, nonspherical, and nonstructural acidic protein, is mainly a kind of PKC substrates taking part in various biological processes involving cytoskeletal control, cell cycle regulation, chemotactic effect, inflammatory reaction, and so on (Lee et al., 2015). In different physiological environments, MARCKS interacts with signals of downstream pathways. Abnormal MARCKS signals may participate in malignant deformation, uncontrolled proliferation, migration, and invasion of malignant tumor cells. Recently, some studies showed that as a proto-oncogene, MARCKS takes part in the tumorigenesis and development of several cancers such as lung cancer (Rohrbach et al., 2015), breast cancer (Chen et al., 2015), renal cell carcinoma (Chen et al., 2017), pancreatic cancer (Brandi et al., 2016), and HCC (Naboulsi et al., 2016) and plays an important role in participation in signaling pathways of many cancers. In addition, MARCKS also has an effect on the sensitivity of cancer cells to chemotherapy drugs and targeted radiation therapy. For these reasons, MARCKS might be an underlying drug target (Chen et al., 2014). FADS2 is one of the FA desaturase gene families.
FADS2 (Ralston et al., 2015), encoding a membrane-binding protein composed of 444 amino acids that cross the cytomembrane four times (Marquardt et al., 2000), has a vital effect on the biosynthesis process of FA (Nakamura and Nara, 2004). The main function of FADS2 is to regulate the synthesis of polyunsaturated fatty acids by introducing a double bond in the hydrocarbon chain of FAs. FADS2 is the first rate-limiting enzyme in polyunsaturated fatty acid synthesis (Lee and Park, 2014). Recent studies showed that FADS2 has a high expression level in breast cancer tissues and esophageal squamous cell carcinoma, which plays a vital role in the occurrence and development of these cancers (Lagarde and Nicolaou, 2015). FADS2 has become a promising drug target in the treatment of cancers. CXCL1, a small molecular cytokine, belongs to the CXC chemokine family, also known as growth regulating oncogene a (GRO a). It produces normal physiological and immune responses by recruiting specific cell groups to infect the specific site of the tumor (Yapa et al., 2017). CXCL1 not only plays a significant role in inflammatory response but also participates in a variety of biological processes, such as angiogenesis tumor genesis (He et al., 2012), atherosclerosis (Boisvert et al., 1998), wound healing, and so on (Ridiandries et al., 2017). Some studies showed that several liver diseases could induce the high expression of CXCL1 in hepatocytes and hepatic astrocytes and eventually lead to tumor genesis, development, invasion, metastasis, and poor prognosis (Miyake et al., 2013). It has also been reported that in the microenvironment of liver cancer, hepatic astrocytes could secrete CXCL1, which helps to induce epithelial–mesenchymal transformation of hepatocarcinoma cells, resistance to regulation, migration, and invasion, and degradation of extracellular matrix by activating the PI3K/AKT pathway (Sobolik et al., 2014). Because of playing an important role in tumor development and potential clinical significance, as a tumor marker, CXCLl might be considered as a therapeutic target of cancer. BMP4, belonging to the BMP subfamily, is a kind of multifunctional signaling molecule that plays a crucial role in regulating important biological functions such as cell growth and development, tissue differentiation, organ formation, apoptosis, and so on (Watanabe et al., 2001). As an upstream signaling molecule, BMP4 carries on its biological functions mainly by activating the intracellular downstream signaling pathway. Some studies showed that in addition to Smad signaling pathway in response to BMPs, BMP4 could also activate downstream p38, JNK, ERK1/2, and other kinases through the interaction between the receptor and MAPKKK, and then p38, JNK, ERK1/2, and other signaling pathways are activated to regulate the transcription of downstream target genes (Selvamurugan et al., 2004). BMP4 can act on tumor growth directly or indirectly by interacting with cytokines or affecting downstream gene expression. On the one hand, BMP4 promotes tumor growth by directly stimulating the expression of genes related to tumor cell proliferation and migration. For instance, in prolactin tumors, BMP4 facilitates the expression of proliferation-related genes such as C-MYC to provoke tumor cell proliferation (Paez-Pereda et al., 2003). On the other hand, BMP4 inhibits tumor growth by interacting with other cytokines. For example, in myeloma cells, BMP4 inhibits DNA synthesis depending on interleukin-6 by downregulating tyrosine phosphorylation of Stat3 induced by interleukin-6 and then inhibits tumor cell proliferation (Hjertner et al., 2001). SKP2, also known as FBXL 1, is an E3 ubiquitin ligase in the ubiquitin degradation system, forming Skp2-SCF ubiquitin ligase complex with Cullin-1, Skp1, and Ring-box 1 for degrading protein substrate. The function of SKP2 as a direct regulator of cyclin-dependent kinases and gene transcription regulators has been established (Frescas and Pagano, 2008). Accumulating studies reported that SKP2 could promote the degradation of several tumor suppressor proteins, including P27, P21, P57, and P130 and also degrade Forkhead box protein O1 by ubiquitination, a transcription factor in the FoxO family with the function of inducing cell cycle arrest, thus to promote the occurrence and development of tumors. It is detected that the protein level of SKP2 is upregulated in various types of human cancers containing liver cancer (Shin et al., 2011), lymphoma, and prostate cancer (Yang et al., 2002). Accumulating studies have reported that SKP2 has a high expression in HCC tissues, and its overexpression might indicate a poor prognosis. KNG1 gene located on chromosome three contains 11 exons expressing high molecular weight kininogene (HK) and low molecular weight kininogene (LK). KNG1 is mainly synthesized in the liver belonging to a plasma protein and a cysteine proteinase inhibitor cleaved into six subchains. KNG1 has been revealed to exert blood coagulation and important anti-inflammatory functions and inhibit the proliferation of endothelial cells and glioma cells (Xu et al., 2018). A recent study demonstrates that KNG1 is connected with the coronavirus disease 2019 infection by reducing ACE2–bradykinin axis related to lung injury and inflammation (Roche and Roche, 2020). Several studies showed that KNG1 had higher content in urine and serum from different cancers, and the content of KNG1 in serum was lower before operation than after operation in colorectal cancer patients (Navaneethan et al., 2014). Recently, KNG1 was also considered a serum biomarker for the early detection of colorectal cancer and advanced colorectal adenoma (Wang et al., 2013). It is noteworthy that the significant decrease in FASN (log2FC= −1.64), MCM5 (log2FC = −1.95), MCM6 (log2FC = −1.90), PCNA (log2FC = −2.31), MARCKS (log2FC = −2.00), FADS2 (log2FC = −3.22), CXCL1 (log2FC = −6.27), BMP4 (log2FC = −2.12), SKP2 (log2FC = −2.03), and KNG1 (log2FC = −3.12), all of which have been related to tumor cell growth, invasion, and metastasis.
Carcinogenesis is a complicated cellular process involved in the dysfunction of various biological pathways (Wang et al., 2013). Through transcriptome analysis, we determined many potential molecular pathways that are significant for carcinogenesis and tumor progression in Huh-7 cells with the treatment of metformin, such as DNA replication, cell cycle, FA biosynthesis, valine, leucine and isoleucine degradation, ECM–receptor interaction, TGF-β signaling pathway, Fanconi anemia pathway, basal cell carcinoma, synthesis and degradation of ketone bodies, microRNAs in cancer, AMPK signaling pathway, cytokine–cytokine receptor interaction, complement and coagulation cascades, TNF signaling pathway, pathways in cancer cGMP-PKG signaling pathway, and PPAR signaling pathway.
In this paper, the pathway of the cell cycle was significantly enriched in downregulated genes of Huh-7 cells treated with metformin (Supplementary Figure S1A). The cell cycle is a pathway associated with cell division mainly driven by a series of genes and other pathways. The regulation of the cell cycle is essential for cell proliferation, differentiation, and apoptosis, especially the compounds synthetized by cyclins and cyclin-dependent kinases. G1/S phase, called Restriction point (R point) in mammals, is the most important checkpoint and controls the transition of a cell from the G1 phase to the S phase. The R point mainly examines whether DNA is damaged or the extracellular environment is suitable, or the cell volume is large enough; therefore, the accurate transition of G1/S is important to the regulation of cell proliferation, and its misregulation may accelerate oncogenesis. SKP2, PCNA, and MCM families (MCM5, MCM6, and MCM7) are significant in the regulation of the G1/S transition (Bertoli et al., 2013) and were significantly downregulated in Huh-7 cells with metformin treatment. The outcomes showed that the arrest of the cell cycle induced by metformin might be a kind of mechanism that causes proliferative inhibition of Huh-7 cells.
DNA and its related substance replication are the basis of cell division, and the perfect mechanism of DNA repair in cells is beneficial to maintaining the integrity of individual genomes. DNA damage caused by intracellular or extracellular stimulation may lead to genomic instability, resulting in cell apoptosis and even malignant transformation of normal cells into tumor ones or other irreversible adverse consequences if it cannot be repaired in time. In our study, the DNA replication pathway in Huh-7 cells treated with metformin was enriched in downregulated genes such as PCNA and MCM families (MCM2, MCM4, MCM5, and MCM6) (Supplementary Figure S1B). The outcomes confirmed that the G1/S phase of cell cycle arrest caused by metformin might attribute to DNA replication.
ECM, as a main constitution of the extracellular microenvironment, plays an important role in the regulation of cell growth, metabolism, migration, proliferation, and differentiation by directly interacting with cells through integrins or other cell surface receptors and is also a significant tissue barrier for tumor metastasis (Wang et al., 2017). A study showed that the formation of ECM tumor fiber formed a cross-linked network structure, which not only helps to prop up and provide nutrition for tumor cells but also plays an important role in tumor growth and invasion. TGF-β produced by tumor cells is a kind of important regulatory factor of activating and promoting fibroblast formation and mediates tumor metastasis by promoting ECM fibrogenesis (Chakravarthy et al., 2018). The basement membrane is a significant part of ECM and acts in the process of ECM promoting tumor cell metastasis. Collagen and laminin are the main components of the basement membrane and play a crucial role in the invasion and development of cancer cells (Kim et al., 2010). In our study, the genes related to the ECM–receptor interaction pathways (Supplementary Figure S1C) are downregulated in Met-treated Huh-7 cells, such as collagen family members (COL1A1, COL1A2, COL2A1, COL9A2, COL9A3, and COL4A1), laminin family members (LAMB3, LAMC2, and LAMA3), and integrin family members (ITGA7 and ITGA5), which are highly consistent with the previous outcomes of migration and invasion assays. The outcomes demonstrated that the inhibition of migration and invasion in Huh-7 cells with metformin treatment might be worked through influencing ECM–receptor interaction pathway.
Conclusion
In summary, this study showed that metformin had a significant inhibitory influence on the proliferation, migration, and invasion of Huh-7 in vitro, which was an outcome of mutual effect among a load of pathways such as DNA replication, cell cycle, ECM–receptor interaction, FA biosynthesis, AMPK signaling pathway, TNF signaling pathway, valine, leucine and isoleucine degradation, TGF-β signaling pathway, pathways in cancer, basal cell carcinoma, PPAR signaling pathway, and multiple of genes related with cell proliferation, cell adhesion, cell cycle, and microRNAs in cancer, substance metabolism, such as MCM family members and DUSP1, FASN, PCNA, SKP2, KNG1, CXCL1, FADS2, MARCKS, and BMP4. In conclusion, this work suggested that metformin might provide some therapeutic effects in the Huh-7 cell line. However, further investigation will be needed to confirm these findings in other HCC cell lines and even in vivo in the future study in the hope of providing a new perspective for future liver cancer therapeutic developments and developing cancer treatment drugs.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/, GSE190076.
Author Contributions
C-QL—Designed and performed the research, analyzed the data, and wrote the manuscript. D-QL and Z-QL—Conceived and designed the research. S-SL and LJ—Collected the literature. G-TZ and CC—Revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by National Natural Science Foundation of China (No. 31672070), Foundation of One Province and One University Project (No. 801260201186), and Changjiang Scholars and Innovative Research Team in University (grant number IRT_15R16).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.822023/full#supplementary-material
References
Abaan, O. D., Polley, E. C., Davis, S. R., Zhu, Y. J., Bilke, S., Walker, R. L., et al. (2013). The Exomes of the NCI-60 Panel: a Genomic Resource for Cancer Biology and Systems Pharmacology. Cancer Res. 73 (14), 4372–4382. doi:10.1158/0008-5472.CAN-12-3342
Abassi, Y. A., Xi, B., Zhang, W., Ye, P., Kirstein, S. L., Gaylord, M. R., et al. (2009). Kinetic Cell-Based Morphological Screening: Prediction of Mechanism of Compound Action and Off-Target Effects. Chem. Biol. 16 (7), 712–723. doi:10.1016/j.chembiol.2009.05.011
Abou-Alfa, G. K., Meyer, T., Cheng, A. L., El-Khoueiry, A. B., Rimassa, L., Ryoo, B. Y., et al. (2018). Cabozantinib in Patients with Advanced and Progressing Hepatocellular Carcinoma. N. Engl. J. Med. 379, 54–63. doi:10.1056/NEJMoa1717002
Alison, M. R., Hunt, T., and Forbes, S. J. (2002). Minichromosome Maintenance (MCM) Proteins May Be Pre-cancer Markers. Gut 50 (3), 290–291. doi:10.1136/gut.50.3.290
Alo, P. L., Amini, M., Piro, F., Pizzuti, L., Sebastiani, V., Botti, C., et al. (2007). Immunohistochemical Expression and Prognostic Significance of Fatty Acid Synthase in Pancreatic Carcinoma. Anticancer Res. 27 (4B), 2523–2527.
Behan, F. M., Iorio, F., Picco, G., Gonçalves, E., Beaver, C. M., Migliardi, G., et al. (2019). Prioritization of Cancer Therapeutic Targets Using CRISPR-Cas9 Screens. Nature 568 (7753), 511–516. doi:10.1038/s41586-019-1103-9
Ben Sahra, I., Regazzetti, C., Robert, G., Laurent, K., Le Marchand-Brustel, Y., Auberger, P., et al. (2011). Metformin, Independent of AMPK, Induces mTOR Inhibition and Cell-Cycle Arrest through REDD1. Cancer Res. 71 (13), 4366–4372. doi:10.1158/0008-5472.CAN-10-1769
Bertoli, C., Skotheim, J. M., and de Bruin, R. A. (2013). Control of Cell Cycle Transcription during G1 and S Phases. Nat. Rev. Mol. Cel Biol. 14 (8), 518–528. doi:10.1038/nrm3629
Boisvert, W. A., Santiago, R., Curtiss, L. K., and Terkeltaub, R. A. (1998). A Leukocyte Homologue of the IL-8 Receptor CXCR-2 Mediates the Accumulation of Macrophages in Atherosclerotic Lesions of LDL Receptor-Deficient Mice. J. Clin. Invest. 101 (2), 353–363. doi:10.1172/JCI1195
Brandi, J., Dalla Pozza, E., Dando, I., Biondani, G., Robotti, E., Jenkins, R., et al. (2016). Secretome Protein Signature of Human Pancreatic Cancer Stem-like Cells. J. Proteomics 136, 1–12. doi:10.1016/j.jprot.2016.01.017
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 68, 394–424. doi:10.3322/caac.21492
Bruix, J., Qin, S., Merle, P., Granito, A., Huang, Y. H., and Bodoky, G. (2017). Regorafenib for Patients with Hepatocellular Carcinoma Who Progressed on Sorafenib Treatment (RESORCE): a Randomised, Double-Blind, Placebo-Controlled, Phase 3 Trial. Lancet 389, 56–66. doi:10.1016/S0140-6736(16)32453-9
Buzzai, M., Jones, R. G., Amaravadi, R. K., Lum, J. J., DeBerardinis, R. J., Zhao, F., et al. (2007). Systemic Treatment with the Antidiabetic Drug Metformin Selectively Impairs P53-Deficient Tumor Cell Growth. Cancer Res. 67 (14), 6745–6752. doi:10.1158/0008-5472.CAN-06-4447
Chakravarthy, A., Khan, L., Bensler, N. P., Bose, P., and Carvalho, D. D., (2018). TGF-β-associated Extracellular Matrix Genes Link Cancer-Associated Fibroblasts to Immune Evasion and Immunotherapy Failure. Nat. Commun. 9 (1), 4692. doi:10.1038/s41467-018-06654-8
Chen, C.-H., Cheng, C.-T., Yuan, Y., Zhai, J., Arif, M., Fong, L. W., et al. (2015). Elevated MARCKS Phosphorylation Contributes to Unresponsiveness of Breast Cancer to Paclitaxel Treatment. Oncotarget 6 (17), 15194–15208. doi:10.18632/oncotarget.3827
Chen, C.-H., Fong, L. W. R., Yu, E., Wu, R., Trott, J. F., and Weiss, R. H. (2017). Upregulation of MARCKS in Kidney Cancer and its Potential as a Therapeutic Target. Oncogene 36 (25), 3588–3598. doi:10.1038/onc.2016.510
Chen, C.-H., Statt, S., Chiu, C.-L., Thai, P., Arif, M., Adler, K. B., et al. (2014). Targeting Myristoylated Alanine-Rich C Kinase Substrate Phosphorylation Site Domain in Lung Cancer. Mechanisms and Therapeutic Implications. Am. J. Respir. Crit. Care Med. 190 (10), 1127–1138. doi:10.1164/rccm.201408-1505OC
Donadon, V., Balbi, M., Mas, M. D., Casarin, P., and Zanette, G. (2010). Metformin and Reduced Risk of Hepatocellular Carcinoma in Diabetic Patients with Chronic Liver Disease. Liver Int. 30 (5), 750–758. doi:10.1111/j.1478-3231.2010.02223.x
Eisenreich, A., and Leppert, U. (2017). Update on the Protective Renal Effects of Metformin in Diabetic Nephropathy. Curr. Medl Chem. 24, 3397–3412. doi:10.2174/0929867324666170404143102
El-Serag, H. B., and Rudolph, K. L. (2007). Hepatocellular Carcinoma: Epidemiology and Molecular Carcinogenesis. Gastroenterology 132, 2557–2576. doi:10.1053/j.gastro.2007.04.061
Evans, J. M. M., Donnelly, L. A., Emslie-Smith, A. M., Alessi, D. R., and Morris, A. D. (2005). Metformin and Reduced Risk of Cancer in Diabetic Patients. BMJ 330, 1304–1305. doi:10.1136/bmj.38415.708634.F7
Farazi, P. A., and DePinho, R. A. (2006). Hepatocellular Carcinoma Pathogenesis: from Genes to Environment. Nat. Rev. Cancer 6, 674–687. doi:10.1038/nrc1934
Frescas, D., and Pagano, M. (2008). Deregulated Proteolysis by the F-Box Proteins SKP2 and Beta-TrCP: Tipping the Scales of Cancer. Nat. Rev. Cancer 8 (6), 438–449. doi:10.1038/nrc2396
Gansler, T. S., Hardman, W., Hunt, D. A., Schaffel, S., Hennigar, R. A., and Hennigar, R. A. (1997). Increased Expression of Fatty Acid Synthase (OA-519) in Ovarian Neoplasms Predicts Shorter Survival. Hum. Pathol. 28 (6), 686–692. doi:10.1016/s0046-8177(97)90177-5
Gaston, K. E., and Grossman, H. B. (2010). Proteomic Assays for the Detection of Urothelial Cancer. Methods Mol. Biol. 641, 303–323. doi:10.1007/978-1-60761-711-2_17
Gonzalez-Guerrico, A. M., Espinoza, I., Schroeder, B., Park, C. H., Kvp, C. M., Khurana, A., et al. (2016). Suppression of Endogenous Lipogenesis Induces Reversion of the Malignant Phenotype and Normalized Differentiation in Breast Cancer. Oncotarget 7 (44), 71151–71168. doi:10.18632/oncotarget.9463
Hjertner, O., Hjorth-Hansen, H., Börset, M., Seidel, C., Waage, A., and Sundan, A. (2001). Bone Morphogenetic Protein-4 Inhibits Proliferation and Induces Apoptosis of Multiple Myeloma Cells. Blood 97 (2), 516–522. doi:10.1182/blood.v97.2.516
Hou, M.-F., Chang, C.-W., Chen, F.-M., Wang, S.-N., Yang, S.-F., Chen, P.-H., et al. (2012). Decreased Total MKP-1 Protein Levels Predict Poor Prognosis in Breast Cancer. World J. Surg. 36 (8), 1922–1932. doi:10.1007/s00268-012-1608-y
Hsieh, S.-C., Tsai, J.-P., Yang, S.-F., Tang, M.-J., and Hsieh, Y.-H. (2014). Metformin Inhibits the Invasion of Human Hepatocellular Carcinoma Cells and Enhances the Chemosensitivity to Sorafenib through a Downregulation of the ERK/JNK-mediated NF-κb-dependent Pathway that Reduces uPA and MMP 9 Expression. Amino Acids 46 (12), 2809–2822. doi:10.1007/s00726-014-1838-4
Ingham, M., and Schwartz, G. K. (2017). Cell-Cycle Therapeutics Come of Age. J. Clin. Oncol. 35 (25), 2949–2959. doi:10.1200/JCO.2016.69.0032
Kastan, M. B., and Bartek, J. (2004). Cell-cycle Checkpoints and Cancer. Nature 432 (18), 316–323. doi:10.1038/nature03097
Kim, B. G., An, H. J., Kang, S., Choi, Y. P., Gao, M. Q., Park, H., et al. (2010). Laminin-332-rich Tumor Microenvironment for Tumor Invasion in the Interface Zone of Breast Cancer. Am. J. Pathol. 178 (1), 373–381. doi:10.1016/j.ajpath.2010.11.028
Kudo, M., Finn, R. S., Qin, S., Han, K. H., Ikeda, K., Piscaglia, F., et al. (2018). Lenvatinib versus Sorafenib in First-Line Treatment of Patients with Unresectable Hepatocellular Carcinoma: a Randomised Phase 3 Non-inferiority Trial. Lancet 391, 1163–1173. doi:10.1016/S0140-6736(18)30207-1
Lagarde, M., and Nicolaou, A. (2015). Oxygenated Metabolism of PUFA: Analysis and Biological Relevance. Biochim. Biophys. Acta 1851 (4), 307. doi:10.1016/j.bbalip.2015.01.005
Lee, H., and Park, W. J. (2014). Unsaturated Fatty Acids, Desaturases, and Human Health. J. Med. Food 17 (2), 189–197. doi:10.1089/jmf.2013.2917
Lee, J. H., Kim, T. I., Jeon, S. M., Hong, S. P., Cheon, J. H., and Kim, W. H. (2012). The Effects of Metformin on the Survival of Colorectal Cancer Patients with Diabetes Mellitus. Int. J. Cancer 131, 752–759. doi:10.1002/ijc.26421
Lee, M.-S., Hsu, C.-C., Wahlqvist, M. L., Tsai, H. N., Chang, Y. H., and Huang, Y. C. (2011). Type 2 Diabetes Increases and Metformin Reduces Total, Colorectal, Liver and Pancreatic Cancer Incidences in Taiwanese: a Representative Population Prospective Cohort Study of 800,000 Individuals. BMC Cancer 11, 20. doi:10.1186/1471-2407-11-20
Lee, S.-M., Suk, K., and Lee, W.-H. (2015). Myristoylated Alanine-Rich C Kinase Substrate (MARCKS) Regulates the Expression of Proinflammatory Cytokines in Macrophages through Activation of P38/JNK MAPK and NF-xB. Cell Immunol 296 (2), 115–121. doi:10.1016/j.cellimm.2015.04.004
Li, D., Yeung, S. C., Hassan, M. M., Konopleva, M., and Abbruzzese, J. L. (2009). Antidiabetic Therapies Affect Risk of Pancreatic Cancer. Gastroenterology 137, 482–488. doi:10.1053/j.gastro.2009.04.013
Liu, Z., Tian, Y., Chen, Q., Zhang, G., Li, C., and Luo, D. Q. (2020). Transcriptome Analysis of MDA-MB-231 Cells Treated with Fumosorinone Isolated from Insect Pathogenic Fungi. Anticancer Agents Med. Chem. 20 (4), 417–428. doi:10.2174/1871520619666191212150322
Llovet, J. M., Ricci, S., Mazzaferro, V., Hilgard, P., Gane, E., Blanc, J. F., et al. (2008). Sorafenib in Advanced Hepatocellular Carcinoma. N. Engl. J. Med. 359, 378–390. doi:10.1056/NEJMoa0708857
Llovet, J. M., Zucman-Rossi, J., Pikarsky, E., Sangro, B., Schwartz, M., Sherman, M., et al. (2016). Hepatocellular Carcinoma. Nat. Rev. Dis. Primers 2, 1–23. doi:10.1038/nrdp.2016.18
Marquardt, A., Stöhr, H., White, K., and Weber, B. H. (2000). cDNA Cloning, Genomic Structure, and Chromosomal Localization of Three Members of the Human Fatty Acid Desaturase Family. Genomics 66 (2), 175–183. doi:10.1006/geno.2000.6196
Miyake, M., Goodison, S., Urquidi, V., Giacoia, E., and Rosser, C. J. (2013). Expression of CXCL1 in Human Endothelial Cells Induces Angiogenesis through the CXCR2 Receptor and the ERK1/2 and EGF Pathways. Lab. Invest. 93 (7), 768–778. doi:10.1038/labinvest.2013.71
Morales, D. R., and Morris, A. D. (2015). Metformin in Cancer Treatment and Prevention. Annu. Rev. Med. 66, 41–413. doi:10.1146/annurev-med-062613-093128
Naboulsi, W., Megger, D. A., Bracht, T., Kohl, M., Turewicz, M., Eisenacher, M., et al. (2016). Quantitative Tissue Proteomics Analysis Reveals Versican as Potential Biomarker for Early-Stage Hepatocellular Carcinoma. J. Proteome Res. 15 (1), 38–47. doi:10.1021/acs.jproteome.5b00420
Nakabayashi, H., Taketa, K., Yamane, T., Oda, M., and Sato, J. (1985). Hormonal Control of Alpha-Fetoprotein Secretion in Human Hepatoma Cell Lines Proliferating in Chemically Defined Medium. Cancer Res. 45 (12 Pt 1), 6379–6383.
Nakamura, M. T., and Nara, T. Y. (2004). Structure, Function, and Dietary Regulation of delta6, delta5, and delta9 Desaturases. Annu. Rev. Nutr. 24, 345–376. doi:10.1146/annurev.nutr.24.121803.063211
Naryzhny, S. N. (2008). Proliferating Cell Nuclear Antigen: a Proteomics View. Cell Mol Life Sci. 65 (23), 3789–3808. doi:10.1007/s00018-008-8305-x
Naryzhny, S. N., and Lee, H. (2007). Characterization of Proliferating Cell Nuclear Antigen (PCNA) Isoforms in normal and Cancer Cells: There Is No Cancer-Associated Form of PCNA. FEBS Lett. 581 (25), 4917–4920. doi:10.1016/j.febslet.2007.09.022
Navaneethan, U., Lourdusamy, V., Venkatesh, P., Willard, B., Sanaka, M. R., and Parsi, M. A. (2014). Bile Proteomics for Differentiation of Malignant from Benign Biliary Strictures: a Pilot Study. Gastroenterol. Rep. (Oxf) 3 (2), 136–143. doi:10.1093/gastro/gou066
Nkontchou, G., Cosson, E., Aout, M., Mahmoudi, A., Bourcier, V., Charif, I., et al. (2011). Impact of Metformin on the Prognosis of Cirrhosis Induced by Viral Hepatitis C in Diabetic Patients. J. Clin. Endocrinol. Metab. 96 (8), 2601–2608. doi:10.1210/jc.2010-2415
Paez-Pereda, M., Giacomini, D., Refojo, D., Nagashima, A. C., Hopfner, U., Grubler, Y., et al. (2003). Involvement of Bone Morphogenetic Protein 4 (BMP-4) in Pituitary Prolactinoma Pathogenesis through a Smad/estrogen Receptor Crosstalk. Proc. Natl. Acad. Sci. U S A. 100 (3), 1034–1039. doi:10.1073/pnas.0237312100
Ralston, J. C., Matravadia, S., Gaudio, N., Holloway, G. P., and Mutch, D. M. (2015). Polyunsaturated Fatty Acid Regulation of Adipocyte FADS1 and FADS2 Expression and Function. Obesity (Silver Spring) 23 (4), 725–728. doi:10.1002/oby.21035
Ridiandries, A., Bursill, C., and Tan, J. (2017). Broad-Spectrum Inhibition of the CC-Chemokine Class Improves Wound Healing and Wound Angiogenesis. Int. J. Mol. Sci. 18 (1), 155. doi:10.3390/ijms18010155
Roche, J. A., and Roche, R. (2020). A Hypothesized Role for Dysregulated Bradykinin Signaling in COVID-19 Respiratory Complications. FASEB J. 34 (6), 7265–7269. doi:10.1096/fj.202000967
Rohrbach, T. D., Jarboe, J. S., Anderson, J. C., Trummell, H. Q., Hicks, P. H., Weaver, A. N., et al. (2015). Targeting the Effector Domain of the Myristoylated Alanine Rich C-Kinase Substrate Enhances Lung Cancer Radiation Sensitivity. Int. J. Oncol. 46 (3), 1079–1088. doi:10.3892/ijo.2014.2799
Sainz, B., Barretto, N., and Uprichard, S. L. (2009). Hepatitis C Virus Infection in Phenotypically Distinct Huh7 Cell Lines. PLoS One 4 (8), e6561. doi:10.1371/journal.pone.0006561
Selvamurugan, N., Kwok, S., Alliston, T., Reiss, M., and Partridge, N. C. (2004). Transforming Growth Factor-Beta 1 Regulation of Collagenase-3 Expression in Osteoblastic Cells by Cross-Talk between the Smad and MAPK Signaling Pathways and Their Components, Smad2 and Runx2. J. Biol. Chem. 279 (18), 19327–19334. doi:10.1074/jbc.M314048200
Semple, J. W., and Duncker, B. P. (2004). ORC-associated Replication Factors as Biomarkers for Cancer. Biotechnol. Adv. 22 (8), 621–631. doi:10.1016/j.biotechadv.2004.06.001
Shi, P., Liu, W., Tala, H., Li, F., Zhang, H., Wang, C., et al. (2017). Metformin Suppresses Triple-Negative Breast Cancer Stem Cells by Targeting KLF5 for Degradation. Cell Discov 3, 17010. doi:10.1038/celldisc.2017.10
Shin, E., Kim, S.-H., Jeong, H.-Y., Jang, J.-J., and Lee, K. (2011). Nuclear Expression of S-phase Kinase-Associated Protein 2 Predicts Poor Prognosis of Hepatocellular Carcinoma. APMIS 120 (5), 349–357. doi:10.1111/j.1600-0463.2011.02838.x
Slattery, M. L., Lundgreen, A., and Wolff, R. K. (2012). MAP Kinase Genes and colon and Rectal Cancer. Carcinogenesis 33 (12), 2398–2408. doi:10.1093/carcin/bgs305
Sobolik, T., Su, Y.-J., Wells, S., Ayers, G. D., Cook, R. S., and Richmond, A. (2014). CXCR4 Drives the Metastatic Phenotype in Breast Cancer through Induction of CXCR2 and Activation of MEK and PI3K Pathways. Mol. Biol. Cel. 25 (5), 566–582. doi:10.1091/mbc.E13-07-0360
Souchek, J. J., Davis, A. L., Hill, T. K., Holmes, M. B., Qi, B., Singh, P. K., et al. (2017). Combination Treatment with Orlistat-Containing Nanoparticles and Taxanes Is Synergistic and Enhances Microtubule Stability in Taxane-Resistant Prostate Cancer Cells. Mol. Cancer Ther. 16 (9), 1819–1830. doi:10.1158/1535-7163.MCT-17-0013
Sun, Q.-L., Zhao, C.-P., Wang, T.-Y., Hao, X.-B., Wang, X.-Y., Zhang, X., et al. (2015). Expression Profile Analysis of Long Non-coding RNA Associated with Vincristine Resistance in colon Cancer Cells by Next-Generation Sequencing. Gene 572, 79–86. doi:10.1016/j.gene.2015.06.087
Ueda, K., Arakawa, H., and Nakamura, Y. (2003). Dual-specifcity Phosphatase 5 (DUSP5) as a Direct Transcriptional Target of Tumor Suppressor P53. Oncogene 22 (36), 5586–5591. doi:10.1038/sj.onc.1206845
Wang, J., Wang, X., Lin, S., Chen, C., Wang, C., Ma, Q., et al. (2013). Identification of Kininogen-1 as a Serum Biomarker for the Early Detection of Advanced Colorectal Adenoma and Colorectal Cancer. PLoS One 8 (7), e70519. doi:10.1371/journal.pone.0070519
Wang, M., Zhao, J., Zhang, L., Wei, F., Lian, Y., Wu, Y., et al. (2017). Role of Tumor Microenvironment in Tumorigenesis. J. Cancer 8 (5), 761–773. doi:10.7150/jca.17648
Watanabe, H., de Caestecker, M. P., and Yamada, Y. (2001). Transcriptional Cross-Talk between Smad, ERK1/2, and P38 Mitogen-Activated Protein Kinase Pathways Regulates Transforming Growth Factor-Beta-Induced Aggrecan Gene Expression in Chondrogenic ATDC5 Cells. J. Biol. Chem. 276 (17), 14466–14473. doi:10.1074/jbc.M005724200
Xu, J., Fang, J., Cheng, Z., Fan, L., Hu, W., Zhou, F., et al. (2018). Overexpression of the Kininogen-1 Inhibits Proliferation and Induces Apoptosis of Glioma Cells. J. Exp. Clin. Cancer Res. 37 (1), 180. doi:10.1186/s13046-018-0833-0
Yang, G., Ayala, G., De Marzo, A., Tian, W., Frolov, A., Wheeler, T. M., et al. (2002). Thompson TC, Harper JW. Elevated Skp2 Protein Expression in Human Prostate Cancer: Association with Loss of the Cyclin-dependent Kinase Inhibitor P27 and PTEN and with Reduced Recurrence-free Survival. Clin. Cancer Res. 8 (11), 3419–3426.
Yapa, S., Mulla, O., Green, V., England, J., and Greenman, J. (2017). The Role of Chemokines in Thyroid Carcinoma. Thyroid 27 (11), 1347–1359. doi:10.1089/thy.2016.0660
Yi, S.-W., Choi, J.-S., Yi, J.-J., Lee, Y.-h., and Han, K. J. (2018). Risk Factors for Hepatocellular Carcinoma by Age, Sex, and Liver Disorder Status: A Prospective Cohort Study in Korea. Cancer 124 (13), 2748–2757. doi:10.1002/cncr.31406
Yu, T., Wang, C., Yang, J., Guo, Y., Wu, Y., and Li, X. (2017). Metformin Inhibits SUV39H1-Mediated Migration of Prostate Cancer Cells. Oncogenesis 6 (5), e324. doi:10.1038/oncsis.2017.28
Zaytseva, Y. Y., Harris, J. W., Mitov, M. I., Kim, J. T., Butterfield, D. A., Lee, E. Y., et al. (2015). Increased Expression of Fatty Acid Synthase Provides a Survival Advantage to Colorectal Cancer Cells via Upregulation of Cellular Respiration. Oncotarget 6 (22), 18891–18904. doi:10.18632/oncotarget.3783
Keywords: metformin, transcriptome analysis, hepatocellular carcinoma, signaling pathway, Huh-7
Citation: Li C-Q, Liu Z-Q, Liu S-S, Zhang G-T, Jiang L, Chen C and Luo D-Q (2022) Transcriptome Analysis of Liver Cancer Cell Huh-7 Treated With Metformin. Front. Pharmacol. 13:822023. doi: 10.3389/fphar.2022.822023
Received: 25 November 2021; Accepted: 18 February 2022;
Published: 23 March 2022.
Edited by:
Madhavi Bandi, Tango Therapeutics, United StatesReviewed by:
Luis Enrique Gomez-Quiroz, Autonomous Metropolitan University, MexicoWunchana Seubwai, Khon Kaen University, Thailand
Copyright © 2022 Li, Liu, Liu, Zhang, Jiang, Chen and Luo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhi-Qin Liu, emhpcWlubGl1OTk5QDEyNi5jb20=; Du-Qiang Luo, ZHUtcWlhbmdsdW9AMTYzLmNvbQ==