- Research article
- Open Access
Integrating transcriptomics and proteomics to show that tanshinone IIA suppresses cell growth by blocking glucose metabolism in gastric cancer cells
BMC Genomics volume 16, Article number: 41 (2015)
Tanshinone IIA (TIIA) is a diterpene quinone extracted from the plant Danshen (Salvia miltiorrhiza) used in traditional Chinese herbal medicine. It has been reported to have anti-tumor potential against several kinds of cancer, including gastric cancer. In most solid tumors, a metabolic switch to glucose is a hallmark of cancer cells, which do this to provide nutrients for cell proliferation. However, the mechanism associated with glucose metabolism by which TIIA acts on gastric cancer cells remains to be elucidated.
We found that TIIA treatment is able to significantly inhibit cell growth and the proliferation of gastric cancer in a dose-dependent manner. Using next-generation sequencing-based RNA-seq transcriptomics and quantitative proteomics-isobaric tags for relative and absolute quantification (iTRAQ), we characterized the mechanism of TIIA regulation in gastric cancer cell line AGS. In total, 16,603 unique transcripts and 102 proteins were identified. After enrichment analysis, we found that TIIA regulated genes are involved in carbohydrate metabolism, the cell cycle, apoptosis, DNA damage and cytoskeleton reorganization. Our proteomics data revealed the downregulation of intracellular ATP levels, glucose-6-phosphate isomerase and L-lactate dehydrogenase B chains by TIIA, which might work with disorders of glucose metabolism and extracellular lactate levels to suppress cell proliferation. The up-regulation of p53 and down-regulation of AKT was shown in TIIA- treated cells, which indicates the transformation of oncogenes. Severe DNA damage, cell cycle arrest at the G2/M transition and apoptosis with cytoskeleton reorganization were detected in TIIA-treated gastric cancer cells.
Combining transcriptomics and proteomics results, we propose that TIIA treatment could lead cell stresses, including nutrient deficiency and DNA damage, by inhibiting the glucose metabolism of cancer cells. This study provides an insight into how the TIIA regulatory metabolism in gastric cancer cells suppresses cell growth, and may help improve the development of cancer therapy.
Gastric cancer is one of the most notorious cancers worldwide. It is the fourth most frequently occurring cancer, and the second most common cause of death for both sexes among all cancers, claiming over 736,000 lives worldwide in 2008 . However, anti-cancer drug development of gastric cancer today remains slow and costly, and drug-resistance remains a potential obstacle .
Tanshinone IIA (TIIA) is a diterpene quinone extracted from the plant Danshen (Salvia miltiorrhiza), which is used in traditional Chinese medicine . The first usage of Danshen extractions as herbal medicine can be traced back to more than 1800 years ago . Its applications include prevention of cardiac diseases , protection of the nervous system  and hepatocytes , and inhibition of osteoporosis . Drug repositioning is considered to be a promising and valuable method for the reduction of the side effects and cost of anti-cancer drug research and development. The functions of many chemicals extracted from Danshen are therefore explored for their anti-cancer potential, and TIIA is the most abundant and structurally representative of these . TIIA has been recently reported to have anti-cancer potential against several cancers including breast cancer , prostate cancer , colorectal cancer , lung cancer , liver cancer , leukemia , and gastric cancer . However, the regulatory mechanism of TIIA in gastric cancer cells remains unclear.
Metabolic transformation accompanying nutritional imbalance is one of the leading causes of cancer progression. Glucose is a primary source for the pentose phosphate pathway, which makes RNA and DNA. Glycolytic intermediates can be used to assist lipid biosynthesis to produce ATP and non-essential amino acids, such as alanine, for growth . Proliferating cells are commonly maintained by enhancing aerobic glycolysis, also called the Warburg effect, which causes lactate accumulation and contributes to the development of malignancies. In oncogenic pathways, phosphatidylinositide 3-kinases (PI3K)/protein kinase B (AKT) has been shown to enhance glycolysis, while the tumor suppressor p53 inhibits it, suggesting a glycolytic switch intrinsically associated with oncogenic transformation . Through the inference of oncogene expression and glucose metabolism, the supply of nutrients to cancer cells could be blocked, which would stunt their proliferative potential. This may be an important finding for drug discovery.
High-throughput data can be used to provide a comprehensive inventory of all the biological processes of cells, display their complexity, and increase data accuracy. An accurate picture of the differential expression of experimental samples is important for defining precise targets and networks. Here we use two types of high-throughput data to uncover the regulatory mechanism of TIIA in gastric cancer cells: transcription levels from next-generation sequencing (NGS) data, and isobaric tags for relative and absolute quantification (iTRAQ)-based quantitative proteomics analysis. NGS is one of newer transcriptome sequencing approaches and can perform high-throughput sequencing by generating thousands or millions of sequences in parallel . It has better sensitivity and lower background noise than microarray analysis, and more researchers are using it to investigate the mechanisms of anti-cancer drugs [19,20]. iTRAQ is a promising new technique for quantitative proteomics which can quantify proteins from different sources using their different tags . Because of its sensitivity in labeling peptides, it is considered a more powerful quantitative proteomic technique than others, including the two frequently used quantitative proteomics techniques 2D-DIGE (difference gel electrophoresis) and cICAT (cleavable isotope-coded affinity tags) [22,23].
In this study, we identified a total of 16,603 unique transcripts and 102 TIIA-regulated proteins that are involved in glucose metabolic process. We further investigated and found that TIIA blocks glycolysis and gluconeogenesis in cancer cells by altering protein expression, and causes DNA damage, cell cycle arrest, cytoskeleton reorganization and apoptosis.
TIIA reduces the growth rate and suppresses proliferation of gastric cancer AGS cells
To determine whether TIIA can affect gastric cancer cell survival, we treated AGS cells with different concentrations to detect real-time cell growth rates by an RTCA (Real-Time Cell Analyzer; xCELLigence™; Roche Applied Science & ACEA Biosciences) system. The IC50 of TIIA at 48 hr was calculated to be 5.3 μM. Figure 1A shows that 0.1 μM TIIA was only slightly effective, whereas 10 μM TIIA killed almost all cells. Each curve is significantly different from the other (p < 0.0001, Wilcoxon Signed-Rank Test). Our results demonstrate that the cell growth rate was significantly decreased under TIIA treatment conditions compared with controls, suggesting that AGS cell growth was significantly inhibited by TIIA in a dose-dependent manner.
To explore whether TIIA could affect gastric cancer cell proliferation rates, we treated AGS cells at different concentrations of TIIA and counted the number of cell colonies stained with crystal violet. The results show that 0.53 μM (1/10 IC50) TIIA reduces colonies to about 60%, whereas 5.3 μM (IC50) TIIA results in almost no colonies. It indicates that the proliferation rate of AGS cells was significantly suppressed by TIIA in a dose-dependent manner (*p < 0.05, **p < 0.01, Student’s t-test) (Figure 1B and C).
Functional annotation enrichment of TIIA -regulated genes
To uncover the TIIA regulatory mechanism in AGS cells, we performed RNA-seq analysis to profile the transcriptomes of gastric cancer cells when treated with dimethyl sulfoxide (DMSO) (control) or with TIIA. In order to characterize the gene expression profile to response of biological functions, our initial analysis performed functional enrichment of all identified genes onto MetaCore pathway analysis. Major bio-functional networks were significantly enriched and showed in Figure 2A. Of six bio-functional networks, “Catabolic process” is top-ranked. Previous studies have shown that dysregulation of metabolism is an important indicator for tumorgenesis . Thus, we forced on metabolic networks to enrich all identified genes by MetaCore. As indicated in Figure 2B, a large number of glucose networks were affected, including Phosphatidylinositol-4,5-diphosphate pathway, Phosphatidylinositol-3,4,5-triphosphate pathway, Ceramide pathway, Pentose phosphate pathways and transport, Glucose pathway, Glycolysis, Gluconeogenesis and glucose transport, and Sucrose metabolism and transport. It has been known that PI3K, an key enzyme of phosphatidylinositol-4,5-diphosphate and phosphatidylinositol-3,4,5-triphosphate pathways, increases glycolysis in cancer cells . Based on these results, we propose that TIIA-regulated pathways involving in the enriched Glycolysis, Gluconeogenesis and glucose transport network (Figure 2C).
Moreover, a total of 2,761 differentially expressed genes (DEGs) were considered that have different expression between control and TIIA treatment samples (FDR ≤ 0.001 and |log2(fold-change)| ≥ 1) (Additional files 1 and 2). Among these, 1,010 genes were up-regulated and 1,751 genes were down-regulated by TIIA. MetaCore networks analysis of DEGs reveals significantly enriched biological processes that frequently occur in gastric cancer cells with TIIA treatment. As shown in Table 1 and Additional file 3, major bio-functional networks were in correspondence with the enrichment analysis of integrative transcriptome (Figure 2A), such as Cell cycle_G2-M, DNA damage_DBS repair, Apoptotic nucleus, Cytoskeleton_Intermediate filaments.
Proteomic expression profiling of TIIA-treated gastric cells
To further elucidate cellular mechanism and molecular function, we performed iTRAQ-based proteomics analysis to assess the protein expression profiles in AGS cells with TIIA treatment. Equal amounts of peptides collected from control samples or TIIA-treated samples were labeled with iTRAQ reagents, respectively, and used for replication (control samples labeled by 114 or 115; TIIA-treated samples labeled by 116 or 117). All labeled peptides were mixed and analyzed by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). A total of 102 differentially expressed proteins were identified with a false discovery rate (FDR) of 3.94%; annotated MS/MS spectra are collected in Table 2 and Additional file 4. There were 100% of identified peptides that were labeled with iTRAQ reagents. The intensity levels of iTRAQ signals distribution plots show that positive correlation among our treatments (Additional file 5), suggesting that labeled samples with a high reproducibility were present in this study.
In order to validate iTRAQ data, we evaluated certain identified protein expression patterns (fold change > 1.5 or < 0.65) by Western blot analysis, including proteasome subunit β type-3 (PSMB3), 40S ribosomal protein S2 (RS2) and glucose-6-phosphate isomerase (G6PI). In our MS/MS spectra, PSMB3 (Figure 3A), RS2 (Figure 3B) and G6PI (Figure 3C) displayed a high intensity across different samples. After western blotting analysis, PSMB3 can be detected to be significantly up-regulated (fold change = 1.45, p < 0.05) whereas RS2 (fold change = 0.49, p < 0.01) and G6PI (fold change = 0.1, p < 0.01) were significantly down-regulated (Figure 3D). As our results, expression patterns of PSMB3, RS2 and G6PI correspond with our iTRAQ data for those proteins, suggesting the consistence between iTRAQ and western blotting data.
TIIA suppresses glucose metabolism of gastric cancer cells
To delineate the mechanism of TIIA in gastric cancer cells, we compared our quantitative proteomics and transcriptomics data. We found TIIA regulated not only genes but also proteins involved in glycolysis, the cell cycle, apoptosis, DNA damage and cytoskeleton rearrangement. As shown in Figure 2C, five proteins (PGAM4, ENO1, LDHA, LDHB and G6PI) were identified at both the protein and transcription level. After western blotting analysis, we found that the protein expression of G6PI and LDHB was downregulated by TIIA, and ENO1 was not changed (Figures 3D and 4A). Because G6PI and LDHB are important enzymes in glycolysis and gluconeogenesis, respectively, we further examined if TIIA can regulate other proteins related to glucose metabolism, such as aldolase C (ALDOC), malate dehydrogenase 1 (MDH1), phosphoenolpyruvate carboxykinase 2 (PCK2) and phosphoglycerate kinase 1 (PGK1), which are involved in gluconeogenesis. ALDOC and PCK2 were also identified in our transcriptomics data (Additional file 2). MDH1, which reduces oxaloacetate to malate in the mitochondria, was downregulated by TIIA. PCK2, which converts oxaloacetate to phosphoenolpyruvate, was up-regulated by TIIA (Figure 4B). The transformation of oncogenes, such as p53 and AKT, is also involved in the glucose metabolism switch in cancer cells . We found that p53 increased and AKT decreased in gastric cancer cells following TIIA treatment (Figure 4C). Additionally, TIIA treatment significantly decreased the intracellular ATP levels in AGS cells compared with the control sample (Figure 4D).
TIIA arrests the cell cycle at the G2/M phase transition
The ability to monitor response to regulation of the cell cycle is an enriched function from our transcriptomics data (Figure 2A and Table 1). We treated AGS cells with different concentrations of TIIA and measured DNA distributions by using flow cytometry to detect the cell cycle distribution of a population of cancer cells. The percentage of AGS cells in the G2/M phase increased up to 13.68% above control levels after 5.3 μM TIIA treatment, showing that TIIA induces cell cycle arrest of AGS at G2/M in a dosage-dependent manner (Figure 5A).
As indicated in Table 1, CDK1, cyclin B1 and Cdc25C are associated with the cell cycle at the G2/M phase, and were identified in our transcriptomics data. CDK1 activation can regulate the progression of the cell cycle from the G2 to the M phase, which is dependent on coordination with cyclin B [24,25]. The activation of the CDK1/cyclin-B complex is maintained through phosphorylation of Thr161 and dephosphorylation of Thr14 and Tyr15 in CDK1. Dephosphorylation of Thr14 and Tyr15 in CDK1 is catalyzed by phosphatase Cdc25C, which is considered a rate-limiting step for the G2 to M phase transition [24,25]. Previous reports suggest monitoring the alteration of CDK1, cyclin B1, Cdc25C, and phospho-CDK1 (CDK1-Thr161) protein expressions is a useful way to validate the occurrence of cell cycle arrest at the G2/M transition . For these reasons, to confirm whether TIIA induces cell cycle arrest at G2/M in gastric cancer cells, we treated AGS cells with TIIA at a concentration of 5.3 μM (IC50), and then measured protein expression levels using western blotting analysis. Levels of phospho-CDK1 (CDK1-Thr161), cyclin B1, and Cdc25C were all reduced in cells treated with TIIA (Figure 5B). Our results indicate that TIIA treatment induced characteristic cell cycle arrest at G2/M in AGS cells by altering cyclin B1 and Cdc25C expression as well as the phosphorylation of CDK1.
TIIA treatment causes apoptosis and reorganization of actin filaments and microtubules
Our cell cycle analysis also indicated that levels of cells in the sub-G1 phase were increased 5.5% above control levels under 5 μM TIIA treatment conditions (Figure 5A). A significant increase of cells in the sub-G1 phase is widely accepted as a sign of apoptosis induction . As our transcriptomics and proteomics data (Figure 2A and Table 1) show, apoptosis-related genes could be induced by TIIA. We treated AGS cells at different concentrations of TIIA and detected the proportions of cells undergoing apoptosis or necrosis using flow cytometry. Our results show that the proportion of cells undergoing apoptosis (including early- and late-phase apoptosis) significantly increased by 15.3% above control levels under 5.3 μM TIIA treatment conditions, showing that TIIA induced apoptosis of AGS cells in a dosage-dependent manner (Figure 6A).
It is widely understood that reorganization of the cytoskeleton, including actin filaments and microtubules, plays a crucial role in apoptosis [28,29]. Links between this process and TIIA treatment can be seen in our transcriptomics and proteomics data (Figure 2A and Table 1). To detect whether the cytoskeletons of AGS cells undergo reorganization after TIIA exposure, we treated AGS cells with TIIA at different concentrations, then examined the consequent distribution of actin filaments and microtubules using immunofluorescence staining. Many cells were seen to manifest shrinking morphology after TIIA treatment (Figure 6B). Actin filaments under TIIA treatment became more condensed, especially at the cell periphery, and underwent crumbling. On the other hand, microtubules aggregated to become thick bundles, and were distributed along nuclear fragmentation sites with condensed chromatin (Figure 6B). These kinds of cytoskeletal reorganizations, combined with nuclear fragmentation, are all characteristic of apoptosis [2,28-30], showing that TIIA induced cytoskeleton reorganization arising from apoptosis in AGS cells.
TIIA triggers γ-H2AX nuclear foci in response to DNA double strand breaks
Based on our previous results, TIIA could induce DNA damage in gastric cancer cells (Figure 2A and Table 1). DNA damage, including double strand breaks (DSB), often leads to genetic instability; proper cellular responses to DNA damage are crucial for cell function and survival [31,32]. Previous studies have shown that phosphorylation of the histone variant H2AX, producing γ-H2AX at nuclear foci, plays an important role in the DNA damage response triggered by DSB [33,34]. The change in H2AX levels was also apparent in our data (Table 1). To examine whether TIIA triggers DNA damage in gastric cancer cells, we treated AGS cells with different levels of TIIA to examine the subsequent localization of γ-H2AX using immunofluorescence staining. Numerous γ-H2AX foci were localized in the nuclei of TIIA treated cells, while γ-H2AX was only represented in a few foci in control cells (Figure 7A). Increased protein expression of γ-H2AX was also detected in TIIA-treated cells (Figure 7B). These results suggest that TIIA triggers DSB, triggering a DNA damage response in AGS cells.
Botanical herbs have been used for disease treatment and prevention, and as alternative and complementary therapies . For instance, paclitaxel isolated from one kind of botanical herb, Taxus brevifolia, has been used as second-line chemotherapy for advanced or recurrent gastric cancer . In this study, we treated gastric cancer with TIIA, which can be used to improve heart function by limiting apoptosis or oxidative damage . Recently, many studies have shown that TIIA exhibits anti-tumor potential [10-16], and that on cancers such as breast cancer, prostate cancer, and leukemia, it can induce the mitochondria-dependent apoptosis pathway [38-40]. In MKN45 gastric cancer cells, TIIA was reported to inhibit cell growth, induce G2/M cell cycle arrest and apoptosis [16,41]. However, few of them actually constructed possible pathways or mechanisms in gastric cancer.
Genes for glycolysis are overexpressed in 70% of all human cancers worldwide, including gastric cancer . By comparing our transcriptomics and proteomics data, we are the first to propose that the transformation of glucose metabolism, glycolysis and gluconeogenesis in gastric cancer cells is a major biological process that is regulated by TIIA treatment. An enzyme involved in glycolysis, G6PI, catalyzes the reversible isomerization of glucose-6-phosphate to fructose-6-phosphate, and is downregulated by TIIA . Increased synthesis of G6PI, known as autocrine motility factor (AMF), is considered to be a unique feature of cancer cells, which stimulates cell growth and contributes to cancer metastasis and malignancy [44-47]. AMF can also down-regulate caspase-9 and Apaf-1, making cancer cells more resistant to mitochondria-dependent apoptosis . Previous studies have reported that TIIA can induce mitochondria-dependent apoptosis by regulating caspase-9 and Apaf-1 in several types of cancer [38-40]. In gastric cancer, we also found that TIIA may induce apoptosis (Figure 6A). Based on our results, we suggest that TIIA might induce the occurrence of apoptosis by suppressing G6PI expression, which decreases glucose consumption and inhibits glycolysis in cancer cells.
In gluconeogenesis, expression of LDHB, one of the subunits of lactate dehydrogenase (LDH), which converts lactate to pyruvate, was shown to decrease in AGS cells after TIIA treatment. Pyruvate is the end product of glycolysis and contributes to gluconeogenesis, acetyl-CoA which enters the Krebs (TCA) cycle, lipid synthesis and nonessential amino acid synthesis for proliferative responses in tumor cells . Decreased LDHB expression might enhance negative regulation of the Warburg effect, which mediates pyruvate into lactate and back in tumor cells . In mouse immortalized cell lines, LDHB is critical in the mechanistic target of rapamycin (mTOR) pathway to induce tumor formation . LDHB is also considered as a tumor marker that increases in many cancers because it facilitates tumor growth and cell proliferation [51-53]. Serum LDH is also considered to be predictors for overall survival of advanced nasopharyngeal carcinoma patients . In gastric cancer cells, we suggest that TIIA treatment might block nutrient supply, which reduces cell survival and proliferation by causing lactate dehydrogenase deficiency.
MDH1 and PCK2 in mitochondria are used to convert pyruvate to phosphoenolpyruvate. In colon cancer, increased MDH levels could enhance glycolysis and lead to cell proliferation and tumorigenesis . Silencing PCK2 in colon cancer cells could reduce susceptibility to preoperative 5-fluorouracil-based radiation therapy  and responses to high glucose levels . In this study, we found that TIIA treatment can decrease MDH1 expression and increase PCK2 expression (Figure 4B). Increased PGK1 in gastric cancer is considered as a biomarker of advanced gastric cancer and is followed by increasing intracellular ATP levels . However, TIIA reduced intracellular ATP levels and slightly up-regulated PGK1, which indicates that some functions of TIIA are still unclear. Moreover, the oncogene AKT was down-regulated by TIIA, which can stimulate aerobic glycolysis in cancer cells , suggesting that TIIA has the potential to be used as chemotherapy for gastric cancer because of its effectiveness in transforming glucose metabolism (Figure 8).
The glucose metabolic switch is associated with the cell cycle, apoptosis and DNA repair in tumor cells . Our transcriptomics analyses and proteomics data also showed that TIIA can regulate these functions and their related proteins, such as PSMB3 and RS2. Previous studies have shown that the expression of PSMB3 and RS2 is involved the development of cancer [60,61]. PSMB3, which belongs to the proteasome B-type family, was up-regulated by TIIA. Proteasomes are necessary for the DNA damage response, and accumulate at double strand break sites to recruit other repair proteins . TIIA treatment increased the likelihood of the occurrence of DNA damage in gastric cancer cells (Figure 7), which might stimulate cancer cells to increase PSMB3 expression.
Moreover, RS2, as a cancer-related ribosomal protein, was down-regulated by TIIA. Increased RS2 is present when p53 mutates, which occurs in cell cycle progression or apoptosis in response to DNA damage . Elevated expression of p53, cell cycle arrest in the G2/M phase, apoptosis and cytoskeleton reorganization were observed in AGS cells following TIIA treatment. DNA damage is one of the major causes of cell cycle arrest in the G2/M phase . Cytoskeleton reorganization, including the redistribution of actin and microtubules, is a characteristic of chromatin remodeling during apoptosis . This has also been observed in cancer cells treated with clinical anti-cancer drugs, such as etoposide, doxorubicin and taxol .
We integrate transcriptomics and proteomics data to uncover TIIA-regulated mechanisms in gastric cancer. In glycolysis and gluconeogenesis, TIIA reduces glucose consumption and the production of pyruvate by regulating G6PI, LDHB, MDH1, PCK2 and PGK1 expression. AKT decreases and p53 increases in response to apoptosis and DNA damage in gastric cancer cells following TIIA treatment, while the glucose metabolism switch induced by the transformation of oncogenes is destroyed (Figure 8). We suggest that TIIA treatment could cause cell stress, nutrient deficiency and DNA damage, to induce cycle arrest at the G2/M transition, apoptosis and cytoskeleton reorganization, all of which would inhibit cell proliferation. In this study, we provide new insight into the TIIA regulatory mechanism in gastric cancer cells, so that future cancer research can improve therapeutic strategies.
Human stomach adenocarcinoma AGS (CRL-1739; ATCC) cells were grown in 90% RPMI 1640 medium (Biological Industries, Beit Haemek, Israel) supplemented with 10% fetal bovine serum (Biological Industries). Cells were cultured at 37°C in an incubator with controlled humidified atmosphere containing 5% CO2. Cells were disaggregated for subculturing with trypsin plus ethylenediaminetetraacetic acid (EDTA), then spun down by centrifugation at 1200 rpm for 5 min, after which the pellet was re-suspended in culture medium.
Tanshinone IIA powder (T4952; Sigma-Aldrich) was dissolved in dimethyl sulfoxide (DMSO) (BioShop), which was diluted to applicable concentrations for treatment.
Growth curves and IC50 measurement
Growth curves of AGS cells were recorded by an RTCA DP® system (xCELLigence™; Roche Applied Science & ACEA Biosciences). Cells were treated with TIIA or 0.1% DMSO as control for 72 hr after 24 hr of seeding (5,000 cells/well in 96-well E-plates). TIIA treatment concentrations were 0.1 μM, 1 μM, 5.3 μM, and 10 μM; 0.1% DMSO dissolved in culture medium was used as control. Each condition was performed in triplicate. The IC50 of TIIA at 48 hr was calculated with RTCA software 1.2 (xCELLigence™) by using AGS cell growth curves under TIIA treatments of 0.1 μM, 1 μM, 5.3 μM, and 10 μM.
Colony formation assay
AGS cells were treated with TIIA or 0.1% DMSO control medium for 5 days after 24 hr of seeding (2,500 cells/well in 6-well plates). TIIA treatment concentrations were 0.53 μM, 1.06 μM, and 5.3 μM; 0.1% DMSO dissolved in culture medium was used as control. After treatment, cells were fixed with methanol for 10 min, stained with 200 μg/mL crystal violet (Sigma-Aldrich) for 10 min, and then colonies were counted using AlphaView SA 3.4.0 (ProteinSimple). Each condition was performed in triplicate.
RNA sequencing and data analysis
Total RNA was extracted from AGS cells treated with DMSO (control) or TIIA, using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s instructions. The quantity and quality of RNA were checked using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA), and were found to have an RNA Integrity Number (RIN) value of more than 9. Poly(A) mRNA was isolated with oligo(dT)-bound magnetic beads and incubated with fragmentation buffer to form short RNA fragments. Reverse transcriptase and random hexamer primers were used to synthesize the first-strand cDNA, and DNA polymerase I and RNaseH were used to synthesize the second-strand cDNA. Double stranded cDNA was end-repaired with T4 DNA polymerase, Klenow fragments, and T4 Polynucleotide Kinase, after which a single “A” base was added, and the entire sequence was ligated to Illumina sequencing adapters (Illumina, San Diego, CA). PCR was performed to amplify the fragments. Using Illumina HiSeqTM 2000 (Illumina), the cDNA library was sequenced on a flow cell after validation on an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System.
After the base calling, the reads with adaptor sequences and low quality scores were removed. The reads with low quality scores are defined as the reads with greater than 10% of unknown bases (N) or greater than 50% of low quality bases which quality values are less than 5. The reads with high quality scores were mapped to the human reference genome hg19 assembly using SOAPaligner in SOAP2 with 2 mismatch allowance  and annotated based on the GENCODE . Read counts for individual GENCODE genes were subsequently determined using HTSeq-count (http://www-huber.embl.de/users/anders/HTSeq), by considering only uniquely mapped reads. Expression of each individual gene was quantified by using the Reads Per Kilobase per Million mapped reads (RPKM) method. The differentially expressed genes were identified by Poisson distribution model . The sequences reported in this study have been deposited in the Sequence Read Archive database with accession number SRP049450.
Functional enrichment analysis
We used MetaCore version 6.18 (GeneGo, St. Joseph, MI, USA) to perform a gene function analysis and to enrich the functional networks of identified genes collected from RNA-seq. There were 16,110 of 16,603 identified genes eligible for network enrichment via GO Processes and Metabolic Networks analysis. For analysis of differentially expressed genes (|log2(fold–change)| > 1), 2,717 of 2,761 DEGs were eligible for GO Processes enrichment analysis, and displayed classified genes (Table 1). All enrichment analysis was tested using the p-value threshold p < 0.0001 for the data inputs.
iTRAQ labeling and LC-MS/MS analysis
AGS cells were treated with 5.3 μM (IC50) TIIA or 0.1% DMSO control medium for 48 hr after 24 hr of seeding (8 × 104 cells/well in 6-well plates). Cells were harvested with trypsin/EDTA and then total proteins were extracted by using lysis buffer (1% SDS (Bioman), 50 mM Tris–HCl (pH 6.8), 10% glycerol), and 1× protease inhibitor (Bioman) and sonication. Proteins were processed through reduction, alkylation, and gel-assisted trypsin digestion overnight to yield peptides, which were later extracted from the gels as previous methods . Equal amounts of peptides from control and TIIA-treated samples were labeled by different iTRAQ reagents (AB SCIEX; control samples labeled by 114 or 115; TIIA-treated samples labeled by 116 or 117) and incubated at room temperature for 1 hr. Peptides were combined together and dried with a centrifugal evaporator (CVE-2000; EYELA). After iTRAQ-labeling, samples were desalted and analyzed by using a LC-ESI-Q-TOF mass spectrometer (Waters SYNAPT® G2 HDMS; Waters Corp.). Samples were injected into a 180 mm × 2 cm capillary trap column and separated by a 75 mm × 25 cm nanoACQUITY UPLC™ 1.7 mm Ethylene Bridged Hybrid C18 column using a nanoACQUITY Ultra Performance LC™ System (Waters Corp.). The mass spectrometer (MS) was operated in electrospray ionization sensitivity mode, and calibrated with a synthetic human [Glu1]-Fibrinopeptide B solution (1 pmol/ml; Sigma-Aldrich) delivered through a NanoLockSpray™ source, which was used for accurate mass measurements. Data was acquired in the data directed analysis (DDA) mode, which included one full MS scan (m/z 350–1700, 1 s) and three sequential MS/MS scans (m/z 100–1990. 1.5 s for each scan) on the three most intense ions present in the full scan mass spectrum.
The peak list resulting from MS/MS spectra was generated by Mascot Distiller v2.3.2 (Matrix Science, London, United Kingdom). Data files were searched against the sequence database (containing 536,789 sequences entries) of the Swiss-Prot human database, using Mascot search engine v2.3.02 (Matrix Science, London, United Kingdom). Both the precursor peptide and fragment ion tolerances were set to ±0.1 Da. The search parameter settings were as follows: allowances for two missed cleavages from trypsin digestion and variable modifications of deamidation (NQ), oxidation (M), iTRAQ (N terminal), iTRAQ (K), and methylthio (C). The peptide charge was set to Mr, the instrument was set to ESI-QUAD-TOF, and the decoy database was searched. Mascot search results were filtered using “Significance threshold” set at p < 0.05 and “Ions score or expect cut-off” set at 0.05. To evaluate the false discovery rate (FDR), we compared a decoy database search against a randomized decoy database created by Mascot using identical search parameters and validation criteria. FDR was calculated as D/R × 100%, where D and R are the number of matches above identity threshold using the decoy and real databases, respectively. The mass spectrometry proteomics data have been uploaded to the ProteomeXchange Consortium  via the PRIDE partner repository with the data set identifier PXD000998 and DOI 10.6019/PXD000998.
For protein quantitation, signature ions (m/z = 114, 115, 116 and 117) and peptides were detected and analyzed using Multi-Q software (v22.214.171.124) . Peptides that satisfied the following four criteria were subjected to further analysis. Firstly, the peptide is labeled with iTRAQ tags; secondly, the peptide has an ion score higher than the Mascot identity score (p < 0.05); thirdly, the peptide is nondegenerate (unique); fourthly, the iTRAQ signature ion peak intensity (ion count) of the peptide is within the dynamic range (the peak intensity of each iTRAQ signature ion must be > 0, and the average of the peak intensities of all iTRAQ signature ions must be ≥ 30). Before quantitation of the expression of each protein, the peak intensity of the iTRAQ signature ion was normalized, as “Method 1” of our previous study . To determine the expression ratio of identified proteins in AGS cells from both the control and the TIIA treatment, the normalized peptide iTRAQ signal of each identified protein was summarized, to calculate protein ratios (TIIA treatment/control).
Western blot analysis
AGS cells were treated with 5.3 μM (IC50) TIIA or 0.1% DMSO control medium for 48 hr after 24 hr of seeding (8 × 104 cells/well in 6-well plates). Cells were harvested with trypsin/EDTA and total proteins were extracted. Then, proteins from control and TIIA-treated samples were separated in 12% SDS-PAGE gels, and transferred onto 0.45 μm PVDF membranes (Millipore) in a Trans-Blot® SD Semi-Dry Transfer Cell (Bio-Rad) for 50 min at 400 mA. The membrane was blocked for 1 hr at room temperature in 5% non-fat milk powder/PBS-T (1× PBS, 0.1% Tween 20 (Sigma-Aldrich)) and incubated overnight at 4°C with blocking buffer containing rabbit monoclonal antibodies to human RS2 (GeneTex; 1:1,000), PSMB3 (GeneTex; 1:1,000), phospho-CDK1 (Santa Cruz; 1:100), CDK1 (Santa Cruz; 1:100), Cyclin B1 (GeneTax; 1:500), Cdc25C (GeneTex; 1:1,000), G6PI (GeneTex; 1:1000), ENO1 (GeneTex; 1:2000), MDH1 (GeneTex; 1:500), PGK1 (GeneTex; 1:500), ALDOC (GeneTex; 1:250), PCK2 (GeneTex; 1:1000), LDH-B (GeneTex; 1:100), p53 (Santa Cruz; 1:500) or AKT (Santa Cruz; 1:1000). The membrane was washed with PBS-T, incubated 1 hr with 5% non-fat milk powder/PBS-T containing anti-rabbit IgG antibodies (1:10,000) (Sigma-Aldrich) or anti-mouse IgG antibodies (Sigma-Aldrich, 1:10,000), washed and imaged with enhanced chemiluminescence (PerkinElmer). The membrane image was then analyzed by an AutoChemi Image System (UVP) or exposed to Fuji medical X-ray film, followed by quantification with AlphaView SA 3.4.0 (ProteinSimple).
Intracellular ATP generation assay
Cells were seeded onto 6-well plates at 8 × 104 cells/well, and then incubated for 24 h. For the control, 0.1% DMSO was added to the medium, and for the treatment, 5.3 μM TIIA was added. After 48 h of drug exposure, the medium was removed, and then cells were washed twice with PBS. The levels of intracellular ATP were determined using a bioluminescent somatic cell assay kit (Sigma-Aldrich), according to the manufacturer’s instructions, and normalized to protein concentrations. Luminescence was detected using a FlexStation III (Molecular Devices). The ATP content of each sample was calculated as the average of the relative light readings and based on the ATP standard curve.
For cell cycle analysis, AGS cells were treated with TIIA or DMSO as control for 48 hr. TIIA treatment concentrations were 0.625 μM, 1.25 μM, 2.5 μM, and 5.3 μM; 0.1% DMSO dissolved in culture medium was used as control. After treatment, cells were harvested with trypsin/EDTA, fixed with 70% ethanol, then spun down, after which ethanol was removed. Then each sample was mixed with RNase A (100 μg/mL), incubated at 37°C for 1 hr, and stained with propidium iodide (PI) (Santa Cruz Biotechnology, Inc.) at a concentration of 100 μg/mL in the dark at room temperature for 15 min. For apoptosis analysis, AGS cells were treated with TIIA or DMSO control medium for 48 hr after 24 hr of seeding (3.5 × 105 cells in 10-cm plates). TIIA treatment concentrations were 1.25 μM, and 5.3 μM; 0.1% DMSO dissolved in culture medium was used as control. After treatment, cells were harvested with trypsin/EDTA, suspended, and counted. Then each sample was adjusted to a concentration of 106 cells/tube and stained with Annexin V-FITC (Santa Cruz Biotechnology, Inc.) and PI (Santa Cruz Biotechnology, Inc.) dissolved in binding buffer (Santa Cruz Biotechnology, Inc.) in the dark at room temperature for 15 min. Both cell cycle distribution and apoptotic cells proportion were then analyzed with a BD FACSCanto II flow cytometer (BD Biosciences) and FCS Express 4 (BD Biosciences).
AGS cells were treated with TIIA or DMSO control medium for 48 hr after 24 hr of seeding (4 × 104 cells in 6-well plates). TIIA treatment concentrations were 1.25 μM and 5.3 μM, and 0.1% DMSO dissolved in culture medium was used as control. Cells were washed, fixed with 4% paraformaldehyde (Sigma-Aldrich) in PBS for 20 min at 37°C, then permeabilized with 0.25% Triton X-100 (Sigma-Aldrich) in PBS for 10 min at room temperature. Then cells were incubated with 1% BSA (Bioshop) in PBS as blocking buffer for 30 min at room temperature, and labeled with mouse monoclonal antibodies to human α-tubulin (Millipore; 1:500) or γ-H2AX (Abcam; 1:500) dissolved in blocking buffer at 4°C overnight. After being washed with PBS three times, cells were labeled with anti-mouse FITC-IgG (Millipore; 1:100) and TRITC-conjugated Phalloidin (Millipore; 1:2000) dissolved in blocking buffer for 1 hr in the dark at room temperature. Then cells were washed with PBS three times and mounted with ProLong® Gold reagent with DAPI (Invitrogen). Images were acquired by using a fluorescence microscope with Leica HCX FL PLAN 1006/1.25 oil objective, a SPOT camera (Diagnostic Instruments), and were analyzed with SPOT Advanced software (Diagnostic Instruments).
Data were expressed as mean ± standard deviation (SD) and analyzed using two-tailed Student’s t-tests. In the cell proliferation assays, data were analyzed using the Wilcoxon Signed-Rank test. A P-value of less than 0.05 was taken to indicate statistical significance.
Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010;127(12):2893–917.
Carden CP, Sarker D, Postel-Vinay S, Yap TA, Attard G, Banerji U, et al. Can molecular biomarker-based patient selection in phase I trials accelerate anticancer drug development? Drug Discov Today. 2010;15(3–4):88–97.
Lin R, Wang WR, Liu JT, Yang GD, Han CJ. Protective effect of tanshinone IIA on human umbilical vein endothelial cell injured by hydrogen peroxide and its mechanism. J Ethnopharmacol. 2006;108(2):217–22.
Chen J, Tina C. Chinese medical herbology and pharmacology. CA: Art of Medicine Press Inc; 2001.
Ding M, Zhao GR, Yuan YJ, Guo ZX. Aqueous extract of Salvia miltiorrhoza regulates adhesion molecule expression of tumor necrosis factor alpha-induced endothelial cells by blocking activation of nuclear factor kappaB. J Cardiovasc Pharmacol. 2005;45(6):516–24.
Ren Y, Houghton PJ, Hider RC, Howes MJ. Novel diterpenoid acetylcholinesterase inhibitors from Salvia miltiorhiza. Planta Med. 2004;70(3):201–4.
Wasser S, Ho JM, Ang HK, Tan CE. Salvia miltiorrhiza reduces experimentally-induced hepatic fibrosis in rats. J Hepatol. 1998;29(5):760–71.
Kim HH, Kim JH, Kwak HB, Huang H, Han SH, Ha H, et al. Inhibition of osteoclast differentiation and bone resorption by tanshinone IIA isolated from Salvia miltiorrhiza Bunge. Biochem Pharmacol. 2004;67(9):1647–56.
Tian G, Zhang Y, Zhang T, Yang F, Ito Y. Separation of tanshinones from Salvia miltiorrhiza Bunge by high-speed counter-current chromatography using stepwise elution. J Chromatogr A. 2000;904(1):107–11.
Wang X, Wei Y, Yuan S, Liu G, Lu Y, Zhang J, et al. Potential anticancer activity of tanshinone IIA against human breast cancer. Int J Cancer. 2005;116(5):799–807.
Won SH, Lee HJ, Jeong SJ, Lu J, Kim SH. Activation of p53 signaling and inhibition of androgen receptor mediate tanshinone IIA induced G1 arrest in LNCaP prostate cancer cells. Phytother Res. 2012;26(5):669–74.
Shan YF, Shen X, Xie YK, Chen JC, Shi HQ, Yu ZP, et al. Inhibitory effects of tanshinone II-A on invasion and metastasis of human colon carcinoma cells. Acta Pharmacol Sin. 2009;30(11):1537–42.
Chiu TL, Su CC. Tanshinone IIA induces apoptosis in human lung cancer A549 cells through the induction of reactive oxygen species and decreasing the mitochondrial membrane potential. Int J Mol Med. 2010;25(2):231–6.
Liu J, Shen HM, Ong CN. Salvia miltiorrhiza inhibits cell growth and induces apoptosis in human hepatoma HepG(2) cells. Cancer Lett. 2000;153(1–2):85–93.
Yoon Y, Kim YO, Jeon WK, Park HJ, Sung HJ. Tanshinone IIA isolated from Salvia miltiorrhiza BUNGE induced apoptosis in HL60 human premyelocytic leukemia cell line. J Ethnopharmacol. 1999;68(1–3):121–7.
Chen J, Shi DY, Liu SL, Zhong L. Tanshinone IIA induces growth inhibition and apoptosis in gastric cancer in vitro and in vivo. Oncol Rep. 2012;27(2):523–8.
Munoz-Pinedo C, El Mjiyad N, Ricci JE. Cancer metabolism: current perspectives and future directions. Cell Death Dis. 2012;3:e248.
Gottlob K, Majewski N, Kennedy S, Kandel E, Robey RB, Hay N. Inhibition of early apoptotic events by Akt/PKB is dependent on the first committed step of glycolysis and mitochondrial hexokinase. Genes Dev. 2001;15(11):1406–18.
Metzker ML. Sequencing technologies - the next generation. Nat Rev Genet. 2010;11(1):31–46.
Lin LL, Huang HC, Juan HF, for the Taida Cancer Systems Biology Study G. Deciphering molecular determinants of chemotherapy in gastrointestinal malignancy using systems biology approaches. Drug Discov Today. 2014;19(9):1402–9.
Wiese S, Reidegeld KA, Meyer HE, Warscheid B. Protein labeling by iTRAQ: a new tool for quantitative mass spectrometry in proteome research. Proteomics. 2007;7(3):340–50.
Wu WW, Wang G, Baek SJ, Shen RF. Comparative study of three proteomic quantitative methods, DIGE, cICAT, and iTRAQ, using 2D gel- or LC-MALDI TOF/TOF. J Proteome Res. 2006;5(3):651–8.
DeSouza L, Diehl G, Rodrigues MJ, Guo J, Romaschin AD, Colgan TJ, et al. Search for cancer markers from endometrial tissues using differentially labeled tags iTRAQ and cICAT with multidimensional liquid chromatography and tandem mass spectrometry. J Proteome Res. 2005;4(2):377–86.
Malumbres M, Barbacid M. Cell cycle, CDKs and cancer: a changing paradigm. Nat Rev Cancer. 2009;9(3):153–66.
Lapenna S, Giordano A. Cell cycle kinases as therapeutic targets for cancer. Nat Rev Drug Discov. 2009;8(7):547–66.
Chang HY, Shih MH, Huang HC, Tsai SR, Juan HF. Middle infrared radiation induces G2/M cell cycle arrest in A549 lung cancer cells. PLoS One. 2013;8:e54117.
Luo J, Su F, Chen D, Shiloh A, Gu W. Deacetylation of p53 modulates its effect on cell growth and apoptosis. Nature. 2000;408(6810):377–81.
Soldani C, Croce AC, Bottone MG, Fraschini A, Biggiogera M, Bottiroli G, et al. Apoptosis in tumour cells photosensitized with Rose Bengal acetate is induced by multiple organelle photodamage. Histochem Cell Biol. 2007;128(5):485–95.
Bottone M, Santin G, Aredia F, Bernocchi G, Pellicciari C. Morphological features of organelles during apoptosis: an overview. Cells. 2013;2:294–305.
Moss DK, Betin VM, Malesinski SD, Lane JD. A novel role for microtubules in apoptotic chromatin dynamics and cellular fragmentation. J Cell Sci. 2006;119(Pt 11):2362–74.
Khanna KK, Jackson SP. DNA double-strand breaks: signaling, repair and the cancer connection. Nat Genet. 2001;27(3):247–54.
Rossi DJ, Bryder D, Seita J, Nussenzweig A, Hoeijmakers J, Weissman IL. Deficiencies in DNA damage repair limit the function of haematopoietic stem cells with age. Nature. 2007;447(7145):725–9.
Celeste A, Petersen S, Romanienko PJ, Fernandez-Capetillo O, Chen HT, Sedelnikova OA, et al. Genomic instability in mice lacking histone H2AX. Science. 2002;296(5569):922–7.
Podhorecka M, Skladanowski A, Bozko P. H2AX phosphorylation: its role in DNA damage response and cancer therapy. J Nucleic Acids. 2010;2010:920161.
Lin LL, Wang YH, Lai CY, Chau CL, Su GC, Yang CY, et al. Systems biology of meridians, acupoints, and chinese herbs in disease. Evid Based Complement Alternat Med. 2012;2012:372670.
Hironaka S, Zenda S, Boku N, Fukutomi A, Yoshino T, Onozawa Y. Weekly paclitaxel as second-line chemotherapy for advanced or recurrent gastric cancer. Gastric Cancer. 2006;9(1):14–8.
Pang H, Han B, Yu T, Peng Z. The complex regulation of tanshinone IIA in rats with hypertension-induced left ventricular hypertrophy. PLoS One. 2014;9(3):e92216.
Tian HL, Yu T, Xu NN, Feng C, Zhou LY, Luo HW, et al. A novel compound modified from tanshinone inhibits tumor growth in vivo via activation of the intrinsic apoptotic pathway. Cancer Lett. 2010;297(1):18–30.
Won SH, Lee HJ, Jeong SJ, Lee HJ, Lee EO, Jung DB, et al. Tanshinone IIA induces mitochondria dependent apoptosis in prostate cancer cells in association with an inhibition of phosphoinositide 3-kinase/AKT pathway. Biol Pharm Bull. 2010;33(11):1828–34.
Liu JJ, Lin DJ, Liu PQ, Huang M, Li XD, Huang RW. Induction of apoptosis and inhibition of cell adhesive and invasive effects by tanshinone IIA in acute promyelocytic leukemia cells in vitro. J Biomed Sci. 2006;13(6):813–23.
Dong X, Dong J, Peng G. Growth-inhibiting and apoptosis-inducing effects of tanshinone II A on human gastric carcinoma cells. J Huazhong Univ Sci Technolog Med Sci. 2007;27(6):706–9.
Altenberg B, Greulich KO. Genes of glycolysis are ubiquitously overexpressed in 24 cancer classes. Genomics. 2004;84(6):1014–20.
Gaitonde MK, Murray E, Cunningham VJ. Effect of 6-phosphogluconate on phosphoglucose isomerase in rat brain in vitro and in vivo. J Neurochem. 1989;52(5):1348–52.
Liotta LA, Mandler R, Murano G, Katz DA, Gordon RK, Chiang PK, et al. Tumor cell autocrine motility factor. Proc Natl Acad Sci U S A. 1986;83(10):3302–6.
Silletti S, Raz A. Autocrine motility factor is a growth factor. Biochem Biophys Res Commun. 1993;194(1):446–57.
Funasaka T, Haga A, Raz A, Nagase H. Tumor autocrine motility factor induces hyperpermeability of endothelial and mesothelial cells leading to accumulation of ascites fluid. Biochem Biophys Res Commun. 2002;293(1):192–200.
Jiang WG, Raz A, Douglas-Jones A, Mansel RE. Expression of autocrine motility factor (AMF) and its receptor, AMFR, in human breast cancer. J Histochem Cytochem. 2006;54(2):231–41.
Haga A, Funasaka T, Niinaka Y, Raz A, Nagase H. Autocrine motility factor signaling induces tumor apoptotic resistance by regulations Apaf-1 and Caspase-9 apoptosome expression. Int J Cancer. 2003;107(5):707–14.
Feron O. Pyruvate into lactate and back: from the warburg effect to symbiotic energy fuel exchange in cancer cells. Radiother Oncol. 2009;92(3):329–33.
Zha X, Wang F, Wang Y, He S, Jing Y, Wu X, et al. Lactate dehydrogenase B is critical for hyperactive mTOR-mediated tumorigenesis. Cancer Res. 2011;71(1):13–8.
Dennison JB, Molina JR, Mitra S, Gonzalez-Angulo AM, Balko JM, Kuba MG, et al. Lactate dehydrogenase B: a metabolic marker of response to neoadjuvant chemotherapy in breast cancer. Clin Cancer Res. 2013;19(13):3703–13.
McCleland ML, Adler AS, Shang Y, Hunsaker T, Truong T, Peterson D, et al. An integrated genomic screen identifies LDHB as an essential gene for triple-negative breast cancer. Cancer Res. 2012;72(22):5812–23.
Chen Y, Zhang H, Xu A, Li N, Liu J, Liu C, et al. Elevation of serum l-lactate dehydrogenase B correlated with the clinical stage of lung cancer. Lung Cancer. 2006;54(1):95–102.
Wan XB, Wei L, Li H, Dong M, Lin Q, Ma XK, et al. High pretreatment serum lactate dehydrogenase level correlates with disease relapse and predicts an inferior outcome in locally advanced nasopharyngeal carcinoma. Eur J Cancer. 2013;49(10):2356–64.
Bi X, Lin Q, Foo TW, Joshi S, You T, Shen HM, et al. Proteomic analysis of colorectal cancer reveals alterations in metabolic pathways: mechanism of tumorigenesis. Mol Cell Proteomics. 2006;5(6):1119–30.
Park JW, Kim SC, Kim WK, Hong JP, Kim KH, Yeo HY, et al. Expression of phosphoenolpyruvate carboxykinase linked to chemoradiation susceptibility of human colon cancer cells. BMC Cancer. 2014;14:160.
Zhao S, Xu W, Jiang W, Yu W, Lin Y, Zhang T, et al. Regulation of cellular metabolism by protein lysine acetylation. Science. 2010;327(5968):1000–4.
Zieker D, Konigsrainer I, Traub F, Nieselt K, Knapp B, Schillinger C, et al. PGK1 a potential marker for peritoneal dissemination in gastric cancer. Cell Physiol Biochem. 2008;21(5–6):429–36.
Elstrom RL, Bauer DE, Buzzai M, Karnauskas R, Harris MH, Plas DR, et al. Akt stimulates aerobic glycolysis in cancer cells. Cancer Res. 2004;64(11):3892–9.
Loging WT, Reisman D. Elevated expression of ribosomal protein genes L37, RPP-1, and S2 in the presence of mutant p53. Cancer Epidemiol Biomarkers Prev. 1999;8(11):1011–6.
Yan MY, Chien SY, Kuo SJ, Chen DR, Su CC. Tanshinone IIA inhibits BT-20 human breast cancer cell proliferation through increasing caspase 12, GADD153 and phospho-p38 protein expression. Int J Mol Med. 2012;29(5):855–63.
Jacquemont C, Taniguchi T. Proteasome function is required for DNA damage response and fanconi anemia pathway activation. Cancer Res. 2007;67(15):7395–405.
Yarden RI, Pardo-Reoyo S, Sgagias M, Cowan KH, Brody LC. BRCA1 regulates the G2/M checkpoint by activating Chk1 kinase upon DNA damage. Nat Genet. 2002;30(3):285–9.
Elmore S. Apoptosis: a review of programmed cell death. Toxicol Pathol. 2007;35(4):495–516.
Grzanka A, Grzanka D, Orlikowska M. Cytoskeletal reorganization during process of apoptosis induced by cytostatic drugs in K-562 and HL-60 leukemia cell lines. Biochem Pharmacol. 2003;66(8):1611–7.
Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, et al. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009;25(15):1966–7.
Harrow J, Frankish A, Gonzalez JM, Tapanari E, Diekhans M, Kokocinski F, et al. GENCODE: the reference human genome annotation for the ENCODE project. Genome Res. 2012;22(9):1760–74.
Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95.
Han CL, Chien CW, Chen WC, Chen YR, Wu CP, Li H, et al. A multiplexed quantitative strategy for membrane proteomics: opportunities for mining therapeutic targets for autosomal dominant polycystic kidney disease. Mol Cell Proteomics. 2008;7(10):1983–97.
Vizcaino JA, Deutsch EW, Wang R, Csordas A, Reisinger F, Rios D, et al. ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nat Biotechnol. 2014;32(3):223–6.
Lin WT, Hung WN, Yian YH, Wu KP, Han CL, Chen YR, et al. Multi-Q: a fully automated tool for multiplexed protein quantitation. J Proteome Res. 2006;5(9):2328–38.
Wu YH, Hu CW, Chien CW, Chen YJ, Huang HC, Juan HF. Quantitative proteomic analysis of human lung tumor xenografts treated with the ectopic ATP synthase inhibitor citreoviridin. PLoS One. 2013;8(8):e70642.
We acknowledged Technology Commons in College of Life Science, National Taiwan University for technical assistance with the Mascot search engine and flow cytometer. This work was supported by the National Science Council of Taiwan (NSC 102-2628-B-002-041-MY3, 102-2627-B-002-002, and 102-2311-B-010-004) and the National Taiwan University Cutting-Edge Steering Research Project (103R7602C3) and Translational Medicine Piloted Model Cooperation Projects of National Taiwan University Medical Campus and College of Life Science (103-L004).
The authors declare that they have no competing interests.
LLL and CRH designed and performed experiments. LLL, CRH, CLH, HCH and HFJ analyzed data. HCH and HFJ contributed reagents and guided all aspects of the project. LLL, CRH, HCH and HFJ co-wrote the paper. All authors read and approved the final manuscript.
The scatter plot represents gene expression levels (mRNA-RPKM) in the control and TIIA treatment samples. Colored dots represent DEGs (FDR ≤ 0.001 and |log2(fold–change)| ≥ 1).
List of RNA-seq identified genes and their expression ratios (TIIA/Control).
The top 50 enrichment functions of RNA-seq-identified DEGs analyzed by MetaCore software.
Detailed information on iTRAQ identified proteins and peptide ratios (TIIA/Control).
High reproducibility of iTRAQ-based proteomics analysis. (A) Control sample peptides were equally divided into two parts, C1 and C2, and labeled with iTRAQ 114 and 115 separately. The distribution plot of iTRAQ ion intensity in control samples (711 peptides identified) shows a high correlation (R2 = 0.975) between the two labels. (B) Treatment (IC50 dose of TIIA for 48 hr) sample peptides were equally divided into two parts, T1 and T2, and labeled with iTRAQ 116 and 117 separately. The distribution plot of iTRAQ ion intensity in treatment samples (711 peptides identified) shows a high correlation (R2 = 0.984) between the two labels.
About this article
Cite this article
Lin, LL., Hsia, CR., Hsu, CL. et al. Integrating transcriptomics and proteomics to show that tanshinone IIA suppresses cell growth by blocking glucose metabolism in gastric cancer cells. BMC Genomics 16, 41 (2015). https://doi.org/10.1186/s12864-015-1230-0
- Tanshinone IIA
- Gastric cancer
- Isobaric tags for relative and absolute quantification
- Next generation sequencing