Skip to main content

Advertisement

Construction of possible integrated predictive index based on EGFR and ANXA3 polymorphisms for chemotherapy response in fluoropyrimidine-treated Japanese gastric cancer patients using a bioinformatic method

BMC Cancer volume 15, Article number: 718 (2015) Cite this article

Article metrics

Abstract

Background

Variability in drug response between individual patients is a serious concern in medicine. To identify single-nucleotide polymorphisms (SNPs) related to drug response variability, many genome-wide association studies have been conducted.

Methods

We previously applied a knowledge-based bioinformatic approach to a pharmacogenomics study in which 119 fluoropyrimidine-treated gastric cancer patients were genotyped at 109,365 SNPs using the Illumina Human-1 BeadChip. We identified the SNP rs2293347 in the human epidermal growth factor receptor (EGFR) gene as a novel genetic factor related to chemotherapeutic response. In the present study, we reanalyzed these hypothesis-free genomic data using extended knowledge.

Results

We identified rs2867461 in annexin A3 (ANXA3) gene as another candidate. Using logistic regression, we confirmed that the performance of the rs2867461 + rs2293347 model was superior to those of the single factor models. Furthermore, we propose a novel integrated predictive index (iEA) based on these two polymorphisms in EGFR and ANXA3. The p value for iEA was 1.47 × 10−8 by Fisher’s exact test. Recent studies showed that the mutations in EGFR is associated with high expression of dihydropyrimidine dehydrogenase, which is an inactivating and rate-limiting enzyme for fluoropyrimidine, and suggested that the combination of chemotherapy with fluoropyrimidine and EGFR-targeting agents is effective against EGFR-overexpressing gastric tumors, while ANXA3 overexpression confers resistance to tyrosine kinase inhibitors targeting the EGFR pathway.

Conclusions

These results suggest that the iEA index or a combination of polymorphisms in EGFR and ANXA3 may serve as predictive factors of drug response, and therefore could be useful for optimal selection of chemotherapy regimens.

Open Peer Review reports

Background

Inter-individual variation in drug response is clinically expected, but relatively difficult to predict [1, 2]. Chemotherapy, in particular, is plagued by highly variable response rates as well as significant toxicity [1]. Genetic variation is an important cause of inter-individual variability in drug response. Dihydropyrimidine dehydrogenase (DPD), an enzyme encoded by the DPYD gene, plays a key role in the adverse effects of fluoropyrimidine treatment: it participates in the catabolism of fluoropyrimidines, such as 5-fluorouracil (5-FU) and its prodrugs capecitabine and S-1 (trade name TS-1, the 5-fluorouracil derivative developed by Tetsuhiko Shirasaka). DPD is an inactivating and rate-limiting enzyme for 5-FU, which is used in various chemotherapeutic regimens to treat gastrointestinal, breast, and head/neck cancers [3]. The antitumor effect of 5-FU is due to its intracellular conversion into antiproliferative nucleotides via anabolic pathways. DPD affects 5-FU availability by rapidly degrading it to 5,6-dihydrofluorouracil (DHFU) [4]. 5-FU catabolism occurs in various tissues including tumors, but is most active in the liver [5, 6].

Wide variability in DPD activity (8- to 21-fold) was shown in Caucasians, and 3–5 % of Caucasians had reduced DPD activity [7, 8]. To date, at least 68 variant DPYD alleles exerting various effects on DPD activity have been reported [3, 9–13]. Of these alleles, the splice site polymorphism IVS14 + 1G>A, which causes skipping of exon 14, is occasionally detected in Northern Europeans with an allele frequency of 0.01–0.02 [9]. Of the patients with a 5-FU-associated grade 3 or 4 adverse event, 24–28 % are heterozygous or homozygous for the IVS14 + 1G>A single nucleotide polymorphism (SNP) [9]. This SNP, however, has not been reported in Japanese or African-American populations [3], and therefore this SNP is not predictive of antitumor effect.

A genome-wide association study (GWAS) is an examination of many common genetic variants in different individuals to determine whether a particular variant is associated with a trait. GWAS using hypothesis-free genomic data is a powerful approach to identify common genetic variants between patients. However, multiple testing problems are a limitation of this approach. We addressed this issue in previous reports by proposing a combined method consisting of a knowledge-based algorithm, two stages of screening, and permutation test to identify significant SNPs [14]. The usability of our combined method was confirmed by applying it into another dataset [15]. In general, the objective of statistical or bioinformatics analysis is the enrichment of important information from a large dataset [16–25]. The use of a knowledge-based algorithm is not a novel concept, but is both practical and useful [26–36]. In the previous study, we applied our combined method to data from gastric cancer patients treated with fluoropyrimidine [14]. We found that rs2293347 in the human epidermal growth factor receptor (EGFR) is a candidate SNP related to chemotherapeutic response and antitumor effect. Nonetheless, the comprehensiveness of the method was limited.

In the present study, to achieve a more comprehensive analysis, we applied our combined method based on an extended knowledge to the dataset of the previous study. Using this approach, we identified rs2867461 in annexin A3 (ANXA3) gene related to the chemotherapeutic response as a novel candidate SNP. Based on discovery of this SNP, we proposed an integrated predictive index based on these two polymorphisms in EGFR and ANXA3 and tested performance of this index. Furthermore, we constructed an EGFR and ANXA3 relation model related to fluoropyrimidine resistance, according to the literature.

Methods

Ethics statement

This study was conducted according to the principles expressed in the Declaration of Helsinki. The ethics committees of the National Cancer Center and National Institute of Health Sciences, Japan, approved the study protocol. All patients provided written informed consent.

Preparation of hypothesis-free genomic data on gastric cancer patients treated with fluoropyrimidine

This study was performed within the framework of the Millennium Genome Project in Japan. A total of 128 Japanese fluoropyrimidine-naïve gastric cancer patients at the National Cancer Center Hospital and National Cancer Center Hospital East were included in the study. DNA samples were extracted from peripheral blood mononuclear cells and 109,365 SNPs were genotyped using the Illumina Human-1 BeadChip. We further restricted our analysis to 119 of the 128 patients whose chemotherapeutic responses were evaluated using Response Evaluation Criteria in Solid Tumors (RECIST). Among the 119 gastric cancer patients, 58 patients were treated with S-1, 27 patients were treated with 5-FU/methotrexate (5-FU/MTX), 33 patients were treated with high-dose 5-FU, and 1 patient was treated with low-dose 5-FU. We defined the 58 patients treated with S-1 as the first dataset and the collection of all 119 patients treated with fluoropyrimidine (including S-1, 5-FU/MTX, high-dose 5-FU, and low-dose 5-FU) as the second dataset in the same way as in the previous study [14].

Patient characteristics and clinical parameters

A summary of the patients’ characteristics from the two datasets is shown in Additional file 1: Table S1. The association of genetic or clinical parameters with chemotherapeutic response was examined using Fisher’s exact test. Chemotherapeutic responses (complete response: CR, partial response: PR, no change: NC, progressive disease: PD) were evaluated using RECIST. We defined two groups: “CR + PR” (CR or PR) and “NC + PD” (NC or PD). Grading of clinical test values was defined using National Cancer Institute - Common Toxicity Criteria (NCI-CTC Version 2.0).

Statistical analyses

Patients’ genotype data and clinical parameters were statistically analyzed by R packages (version 3.1.2) (http://www.r-project.org/). Further detailed theories and algorithms are shown in Additional file 2.

Results

Identification of rs2867461 in ANXA3

We reanalyzed hypothesis-free genomic data from gastric cancer patients treated with fluoropyrimidine by applying applied our combined method with extended knowledge as described in our previous study [14], as shown in Fig. 1. Using this approach, we extracted rs2867461 in ANXA3 as another candidate SNP related to chemotherapeutic response. Further detailed analyses and the procedure are shown in Additional file 3.

Fig. 1
figure1

Extraction of candidate SNPs by an extended KB-SNP. We performed extended KB-SNP to identify novel candidate SNPs related to chemotherapy response. a SNPs linked to any PubMed IDs were extracted and the SNPs related to cancer were removed, as we had already analyzed SNPs related to cancer in the previous study. b A total of 1,767 SNPs were extracted from 109,365 SNPs by the extended KB-SNP and the basic filtering in the present study

Full size image

Comparison of the models based on rs2867461 in ANXA3

We analyzed not only an allele model, but also dominant and recessive models of rs2867461 in ANXA3 in the first (S1-treated gastric cancer patients) and second datasets (fluoropyrimidine-treated gastric cancer patients; Fig. 2). Figure 2a shows that in the first dataset the p value of the allele model was the lowest (p = 1.02 × 10−6, OR = 0.084), and the p value of the recessive model (p = 2.50 × 10−5, OR = 0.033) was lower than the p value of the dominant model (p = 3.24 × 10−4, OR = 0). Similarly, Fig. 2b shows that in the second dataset the p value of the allele model was also the lowest (p = 5.75 × 10−5, OR = 0.22), and the p value of the recessive model (p = 3.52 × 10−4, OR = 0.13) was lower than the p value of the dominant model (p = 7.78 × 10−4, OR = 0.15). Therefore, the recessive model is the best model for rs2867461 in ANXA3. To evaluate combination effects of multiple factors, the proportional odds model was used to construct multiple logistic regression models.

Fig. 2
figure2

Contingency tables for rs2867461 in ANXA3 for each model using each dataset. a S-1-treated gastric cancer patients (first dataset). b Fluoropyrimidine (including S-1)-treated gastric cancer patients (second dataset). P values were calculated using Fisher’s exact test. OR: odds ratio, CI: confidence interval, RECIST: Response Evaluation Criteria in Solid Tumors, CR: complete response, PR: partial response, NC: no change, PD: progressive disease

Full size image

Selection of a model based on rs2867461 in ANXA3 and construction of multiple regression models

We compared AICs and AUCs between 10 models: NULL (without parameters), rs2293347 (genotype of rs2293347 in EGFR), Cr (grade of creatinine), Chem (a history of chemotherapy), rs2867461 (the genotype of rs2867461 in ANXA3), rs2867461 + rs2293347, rs2867461 + rs2293347 + Cr, rs2867461 + rs2293347 + Chem, and rs2867461 + rs2293347 + Cr + Chem model (Fig. 3a). ROC curves for the five logistic regression models, Cr + Chem, rs2867461, rs2293347, rs2867461 + rs2293347, and rs2867461 + rs2293347 + Cr, are shown in Fig. 3b. All models performed better than the NULL model, although the Cr + Chem model was better than either Cr or Chem alone, and the rs2293347 and rs2867461 models performed better than the Cr + Chem model, as shown in Fig. 3a and b. Finally, the rs2867461 + rs2293347 + Cr model had the lowest AIC among the 10 models tested. Although the rs2867461 + rs2293347 + Cr model gave the best results, the best cutoff value was at a sensitivity of 68.0 % and specificity of 100.0 %, with performance depending on only rs2867461 + rs2293347, as shown in Fig. 3b. Therefore, we selected the rs2867461 + rs2293347 model as the best model in the present study, and the best cutoff value was found at a sensitivity of 69.0 % and specificity of 100.0 %. The integrated genetic factor consisting of rs2867461 and rs2293347 is a possible predictive factor of efficacy of treatment in fluoropyrimidine-treated gastric cancer patients.

Fig. 3
figure3

Comparison of AIC, AUC, and ROC curves between logistic regression models. a Parameters used for each model. b ROC curves for the following models: rs2293347, rs2867461, Cr + Chem, rs2867461 + rs2293347, and rs2867461 + rs2293347 + Cr. ROC: receiver operating characteristic, AUC: area under the ROC curve, NULL: model without any parameters. Each genetic factor indicates proportional odds model, AIC: Akaike’s information criterion, Sens.: sensitivity (%), Spec: specificity (%), Chem: a history of chemotherapy, Cr: grade of creatinine

Full size image

The integrated predictive index based on two polymorphisms in EGFR and ANXA3

To define a novel predictive factor consisting of two polymorphisms in EGFR and ANXA3, we defined the total number of minor alleles of rs2293347 and rs2867461 as an integrated predictive index based on EGFR and ANXA3 (iEA index). Contingency tables and the ROC curve for this novel predictive factor, iEA, are shown in Fig. 4. This figure shows that the p value of iEA was 2.56 × 10−8 by Fisher’s exact test, and a higher iEA was correlated with a formula for the better response rate (RR): ((CR + PR)/(CR + PR + NC + PD)). For example, RR = 0 % (iEA = 0 or 1), 28.1 % (iEA = 2), 46.2 % (iEA = 3), and 75.0 % (iEA = 4). Figure 4b shows that the ROC curve for the regression model based on iEA is approximately the same as the ROC for the rs2867461 + rs2293347 model. We constructed a 2 × 2 contingency table by combining contingency tables of iEA, as shown in Fig. 4c. Figure 4c shows that the p value of iEA was 1.47 × 10−8 by Fisher’s exact test. These results suggested that iEA may be an important predictive factor of response rate in fluoropyrimidine-treated gastric cancer patients. Nonetheless, clinical utility of iEA needs to be validated in future studies.

Fig. 4
figure4

Contingency tables for integrated predictive index using polymorphisms in EGFR and ANXA3 and ROC curve. a Contingency table for the iEA index. b ROC curve for the iEA index. c The combined contingency table for the iEA index. Abbreviations are the same as defined in Figs. 2 and 3

Full size image

Discussion

In the previous study, we extracted RS numbers (SNP IDs) related to cancer using a combination of National Center for Biotechnology Information (NCBI) dbSNP and NCBI PubMed [14]. In the present study, we extracted all SNP numbers linked to PubMed IDs on the basis of dbSNP but excluded SNPs related to cancer, as we had already analyzed SNPs related to cancer in the previous study. However, among these SNPs not directly related to cancer, the SNPs could still be indirectly related to cancer, as they may be involved in cellular differentiation, apoptosis, drug metabolism, transporter and immune system processes. Thus, this information may be potentially useful. Therefore, we used information of SNPs linked to any function except for cancer in the present study. Furthermore, Illumina Human-1 BeadChip is one of the most preliminary types of arrays; their detectable SNPs are not tag SNPs and it is difficult to reduce multiple comparisons problem by constructing linkage disequilibrium blocks. Therefore, we focused on the combination of dbSNP and PubMed as the most reliable information.

An SNP extracted using the combined method, rs2867461 in ANXA3, was previously reported as a genetic factor associated with rheumatoid arthritis, systemic lupus erythematosus, and Graves’ disease in a Japanese population [37]. Although the relationship between cancer and rs2867461 in ANXA3 has not been reported to date, many studies have recently been published on the association between ANXA3 and drug resistance or chemotherapy response [38]. The annexin family is a well-known multigene family of Ca2+-regulated phospholipid- and membrane-binding proteins [39]. ANXA3 is a member of the annexin family, and important functions of ANXA3 in tumor development, metastasis, and drug resistance have been demonstrated [38]. For example, ANXA3 overexpression was found to correlate with enhanced drug resistance in ovarian cancer, promote the development of colorectal adenocarcinoma and pancreatic carcinoma, and facilitate metastasis of lung adenocarcinoma and hepatocarcinoma. In contrast, decreased ANXA3 expression negatively correlates with the development of prostate and renal carcinoma [38]. To identify drug resistance mechanisms, Pénzváltó et al. tested 45 cancer cell lines for sensitivity to five tyrosine kinase inhibitors targeting the ERBB/RAS pathway: sunitinib, erlotinib, lapatinib, sorafenib, and gefitinib [40]. The authors identified ANXA3 as one of the two significant genes from microarray analysis and this finding was validated by quantitative real-time PCR. To identify key proteins related to multidrug resistance (MDR) of hepatocellular carcinoma, Tong et al. analyzed the 5-FU-resistant BEL7402/5-FU cell line and parental BEL7402 cells [41]. Among the highly expressed proteins in BEL7402/5-FU associated with MDR, only the expression of ANXA3 was verified using an isobaric tag for relative and absolute quantitation-coupled two-dimensional liquid chromatography tandem mass spectrometry. Furthermore, in a recent study that compared EGFR-mutated and EGFR-wild type tumors, ANXA3 was identified as one of only four downregulated genes involved in prostate cancer progression [42]. These and other results suggest that ANXA3 is a tyrosine phosphorylation target of EGFR [43] and expression of EGFR may generally suppress expression of ANXA3 [42]. Therefore, high expression of ANXA3 may confer drug resistance.

According to our previous report, the rs2293347 SNP in EGFR was extracted as a potential predictive factor of chemotherapeutic response in Japanese gastric cancer patients treated with fluoropyrimidine [14]. This study showed that the rs2293347GA/AA genotype was associated with a lower risk of progressive disease compared with the rs2293347GG genotype (OR = 0.048, p = 6.32 × 10−5). Recently, Mochinaga et al. reported that high expression of DPD in lung adenocarcinoma is associated with mutations in EGFR [44]. Several studies have demonstrated that high DPD levels result in low sensitivity to fluoropyrimidine for various cancers, such as gastric cancer [45, 46], colon cancer [47], bladder cancer [48], and breast cancer [49]. Therefore, rs2293347 might affect DPD expression related to sensitivity to fluoropyrimidine.

The rs2293347G>A polymorphism located in exon 25 of EGFR is a synonymous SNP (D994D), while the rs2867461G>A polymorphism is located in intron 7 of ANXA3. These polymorphisms do not change the amino acid sequence of the protein. However, if rs2867461 and rs2293347 have no function, these SNPs are possible predictive factors linked with other functional polymorphisms in ANXA3 and EGFR, respectively. Therefore, rs2867461 in ANXA3 and rs2293347 in EGFR are promising predictive factors that can be used for selection of chemotherapy regimens: for instance, fluoropyrimidine alone or a combination of fluoropyrimidine with EGFR-targeting agents. Further research is needed to elucidate the clinical relevance of these SNPs.

As mentioned above, many studies suggest that the EGFR and ANXA3 genes have relevance to fluoropyrimidine resistance and their polymorphisms have links with biological functions. Because the IntPath database is currently the most powerful tool and also the most comprehensive integrated pathway database, we first conducted pathway analysis using the IntPath database [50] to draw the genetic networks related to EGFR and ANXA3. However, we could not identify pathway information using this database. Therefore, we manually constructed a hypothetical model of relationship between EGFR and ANXA3 (Fig. 5) according to the literature.

Fig. 5
figure5

Hypothetical model of EGFR and ANXA3 to fluoropyrimidine resistance in fluoropyrimidine-treated gastric cancer patients. ANXA3 overexpression confers resistance tyrosine kinase inhibitors targeting ERBB/RAS pathway. High expression of DPD is associated with mutations in EGFR. DPD is an inactivating and rate-limiting enzyme for fluoropyrimidine

Full size image

In this study, we extracted rs2867461 (which showed statistical significance according to p (0.0406) < 0.05) using a combination of two stages of screening and permutation testing of prefiltered SNPs for both of first and second sets. When only the first dataset was used, the q value calculated by the BH method was 0.00159, as shown in Additional file 4: Table S2. This q value is statistical significance.

Using our combined method involving two stages of screening, we identified rs2867461 as a possible genetic predictive factor. We note that our filtering methodology may have also eliminated several interesting regulatory marker SNPs that might be relevant to drug response, as shown in Fig. 1. However, the sample size of this study is not enough to identify all of these marker SNPs without omission. Therefore, we prioritized control of type I error at the cost of statistical power (type II error) in the present study. All statistical information regarding the chemotherapeutic response of gastric cancer patients treated with fluoropyrimidine (p < 0.05) for each SNP is shown in Additional file 5: Table S3, and the data are also provided on the website Genome Medicine Database of Japan (GeMDBJ) [51] (http://gemdbj.ncc.go.jp/omics/). These data will be useful for confirmation studies or meta-analyses in the future.

Conclusions

In the present study, we reanalyzed hypothesis-free genomic data from gastric cancer patients treated with fluoropyrimidine by applying our combined method with extended knowledge. Using this approach, we identified rs2867461 in ANXA3 as a candidate SNP related to response to chemotherapeutic response. The rs2867461 + rs2293347 model has greater predictive performance than clinical parameters, each single SNP (rs2867461/rs2293347), or environmental factors, and the rs2867461 + rs2293347 model had a sensitivity of 69.0 % and specificity of 100.0 %. Furthermore, in the present study, we propose a novel integrated predictive index based on the polymorphisms in EGFR and ANXA3, the iEA index. The p value for iEA is 1.47 × 10−8 by Fisher’s exact test. Collectively, iEA or the combination of rs2867461 and rs2293347 may serve as predictive factors for selecting chemotherapy regimens for the treatment of gastric cancer patients.

Availability of supporting data

The data set supporting the results of this article is available in the Genome Medicine Database of Japan (GeMDBJ) (http://gemdbj.ncc.go.jp/omics/) with the accession number GWAS030.

References

  1. 1.

    Efferth T, Volm M. Pharmacogenetics for individualized cancer chemotherapy. Pharmacol Ther. 2005;107(2):155–76.

  2. 2.

    Evans WE, Relling MV. Pharmacogenomics: translating functional genomics into rational therapeutics. Science. 1999;286(5439):487–91.

  3. 3.

    Maekawa K, Saeki M, Saito Y, Ozawa S, Kurose K, Kaniwa N, et al. Genetic variations and haplotype structures of the DPYD gene encoding dihydropyrimidine dehydrogenase in Japanese and their ethnic differences. J Hum Genet. 2007;52(10):804–19.

  4. 4.

    Heggie GD, Sommadossi JP, Cross DS, Huster WJ, Diasio RB. Clinical pharmacokinetics of 5-fluorouracil and its metabolites in plasma, urine, and bile. Cancer Res. 1987;47(8):2203–6.

  5. 5.

    Naguib FN, el Kouni MH, Cha S. Enzymes of uracil catabolism in normal and neoplastic human tissues. Cancer Res. 1985;45(11 Pt 1):5405–12.

  6. 6.

    Lu Z, Zhang R, Diasio RB. Dihydropyrimidine dehydrogenase activity in human peripheral blood mononuclear cells and liver: population characteristics, newly identified deficient patients, and clinical implication in 5-fluorouracil chemotherapy. Cancer Res. 1993;53(22):5433–8.

  7. 7.

    Etienne MC, Lagrange JL, Dassonville O, Fleming R, Thyss A, Renee N, et al. Population study of dihydropyrimidine dehydrogenase in cancer patients. J Clin Oncol. 1994;12(11):2248–53.

  8. 8.

    Lu Z, Zhang R, Carpenter JT, Diasio RB. Decreased dihydropyrimidine dehydrogenase activity in a population of patients with breast cancer: implication for 5-fluorouracil-based chemotherapy. Clin Cancer Res. 1998;4(2):325–9.

  9. 9.

    van Kuilenburg AB. Dihydropyrimidine dehydrogenase and the efficacy and toxicity of 5-fluorouracil. Eur J Cancer. 2004;40(7):939–50.

  10. 10.

    Gross E, Ullrich T, Seck K, Mueller V, de Wit M, von Schilling C, et al. Detailed analysis of five mutations in dihydropyrimidine dehydrogenase detected in cancer patients with 5-fluorouracil-related side effects. Hum Mutat. 2003;22(6):498.

  11. 11.

    Ogura K, Ohnuma T, Minamide Y, Mizuno A, Nishiyama T, Nagashima S, et al. Dihydropyrimidine dehydrogenase activity in 150 healthy Japanese volunteers and identification of novel mutations. Clin Cancer Res. 2005;11(14):5104–11.

  12. 12.

    Seck K, Riemer S, Kates R, Ullrich T, Lutz V, Harbeck N, et al. Analysis of the DPYD gene implicated in 5-fluorouracil catabolism in a cohort of Caucasian individuals. Clin Cancer Res. 2005;11(16):5886–92.

  13. 13.

    Zhu AX, Puchalski TA, Stanton Jr VP, Ryan DP, Clark JW, Nesbitt S, et al. Dihydropyrimidine dehydrogenase and thymidylate synthase polymorphisms and their association with 5-fluorouracil/leucovorin chemotherapy in colorectal cancer. Clin Colorectal Cancer. 2004;3(4):225–34.

  14. 14.

    Takahashi H, Kaniwa N, Saito Y, Sai K, Hamaguchi T, Shirao K, et al. Identification of a candidate single-nucleotide polymorphism related to chemotherapeutic response through a combination of knowledge-based algorithm and hypothesis-free genomic data. J Biosci Bioeng. 2013;116(6):768–73.

  15. 15.

    Takahashi H, Sai K, Saito Y, Kaniwa N, Matsumura Y, Hamaguchi T, et al. Application of a combination of a knowledge-based algorithm and 2-stage screening to hypothesis-free genomic data on irinotecan-treated patients for identification of a candidate single nucleotide polymorphism related to an adverse effect. PLoS One. 2014;9(8), e105160.

  16. 16.

    Takahashi H, Honda H. Modified signal-to-noise: a new simple and practical gene filtering approach based on the concept of projective adaptive resonance theory (PART) filtering method. Bioinformatics. 2006;22(13):1662–4.

  17. 17.

    Takahashi H, Iwakawa H, Nakao S, Ojio T, Morishita R, Morikawa S, et al. Knowledge-based fuzzy adaptive resonance theory and its application to the analysis of gene expression in plants. J Biosci Bioeng. 2008;106(6):587–93.

  18. 18.

    Takahashi H, Kobayashi T, Honda H. Construction of robust prognostic predictors by using projective adaptive resonance theory as a gene filtering method. Bioinformatics. 2005;21(2):179–86.

  19. 19.

    Takahashi H, Masuda K, Ando T, Kobayashi T, Honda H. Prognostic predictor with multiple fuzzy neural models using expression profiles from DNA microarray for metastases of breast cancer. J Biosci Bioeng. 2004;98(3):193–9.

  20. 20.

    Takahashi H, Murase Y, Kobayashi T, Honda H. New cancer diagnosis modeling using boosting and projective adaptive resonance theory with improved reliable index. Biochem Eng J. 2007;33(2):100–9.

  21. 21.

    Takahashi H, Nakagawa A, Kojima S, Takahashi A, Cha BY, Woo JT, et al. Discovery of novel rules for G-quadruplex-forming sequences in plants by using bioinformatics methods. J Biosci Bioeng. 2012;114(5):570–5.

  22. 22.

    Takahashi H, Nakayama R, Hayashi S, Nemoto T, Murase Y, Nomura K, et al. Macrophage migration inhibitory factor and stearoyl-CoA desaturase 1: potential prognostic markers for soft tissue sarcomas based on bioinformatics analyses. PLoS One. 2013;8(10), e78250.

  23. 23.

    Takahashi H, Takahashi A, Naito S, Onouchi H. BAIUCAS: a novel BLAST-based algorithm for the identification of upstream open reading frames with conserved amino acid sequences and its application to the Arabidopsis thaliana genome. Bioinformatics. 2012;28(17):2231–41.

  24. 24.

    Takahashi H, Tomida S, Kobayashi T, Honda H. Inference of common genetic network using fuzzy adaptive resonance theory associated matrix method. J Biosci Bioeng. 2003;96(2):154–60.

  25. 25.

    Takahashi H, Honda H. A new reliable cancer diagnosis method using boosted fuzzy classifier with a SWEEP operator method. J Chem Eng Jpn. 2005;38(9):763–73.

  26. 26.

    Kawamura T, Takahashi H, Honda H. Proposal of new gene filtering method, BagPART, for gene expression analysis with small sample. J Biosci Bioeng. 2008;105(1):81–4.

  27. 27.

    Takahashi H, Honda H. Lymphoma prognostication from expression profiling using a combination method of boosting and projective adaptive resonance theory. J Chem Eng Jpn. 2006;39(7):767–71.

  28. 28.

    Sano M, Aoyagi K, Takahashi H, Kawamura T, Mabuchi T, Igaki H, et al. Forkhead box A1 transcriptional pathway in KRT7-expressing esophageal squamous cell carcinomas with extensive lymph node metastasis. Int J Oncol. 2010;36(2):321–30.

  29. 29.

    Takahashi H, Iwakawa H, Ishibashi N, Kojima S, Matsumura Y, Prananingrum P, et al. Meta-analyses of microarrays of arabidopsis asymmetric leaves1 (as1), as2 and their modifying mutants reveal a critical role for the ETT pathway in stabilization of adaxial-abaxial patterning and cell division during leaf development. Plant Cell Physiol. 2013;54(3):418–31.

  30. 30.

    Kumasaka MY, Yajima I, Iida M, Takahashi H, Inoue Y, Fukushima S, et al. Correlated expression levels of endothelin receptor B and Plexin C1 in melanoma. Am J Cancer Res. 2015;5(3):1117–23.

  31. 31.

    Takahashi H, Aoyagi K, Nakanishi Y, Sasaki H, Yoshida T, Honda H. Classification of intramural metastases and lymph node metastases of esophageal cancer from gene expression based on boosting and projective adaptive resonance theory. J Biosci Bioeng. 2006;102(1):46–52.

  32. 32.

    Yajima I, Kumasaka MY, Naito Y, Yoshikawa T, Takahashi H, Funasaka Y, et al. Reduced GNG2 expression levels in mouse malignant melanomas and human melanoma cell lines. Am J Cancer Res. 2012;2(3):322–9.

  33. 33.

    Takahashi A, Nakayama R, Ishibashi N, Doi A, Ichinohe R, Ikuyo Y, et al. Analysis of gene expression profiles of soft tissue sarcoma using a combination of knowledge-based filtering with integration of multiple statistics. PLoS One. 2014;9(9), e106801.

  34. 34.

    Takahashi H, Honda H. Prediction of peptide binding to major histocompatibility complex class II molecules through use of boosted fuzzy classifier with SWEEP operator method. J Biosci Bioeng. 2006;101(2):137–41.

  35. 35.

    Iwasaki M, Takahashi H, Iwakawa H, Nakagawa A, Ishikawa T, Tanaka H, et al. Dual regulation of ETTIN (ARF3) gene expression by AS1-AS2, which maintains the DNA methylation level, is involved in stabilization of leaf adaxial-abaxial partitioning in Arabidopsis. Development. 2013;140(9):1958–69.

  36. 36.

    Takahashi H, Nemoto T, Yoshida T, Honda H, Hasegawa T. Cancer diagnosis marker extraction for soft tissue sarcomas based on gene expression profiling data by using projective adaptive resonance theory (PART) filtering method. BMC Bioinformatics. 2006;7:399.

  37. 37.

    Okada Y, Terao C, Ikari K, Kochi Y, Ohmura K, Suzuki A, et al. Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population. Nat Genet. 2012;44(5):511–6.

  38. 38.

    Wu N, Liu S, Guo C, Hou Z, Sun MZ. The role of annexin A3 playing in cancers. Clin Transl Oncol. 2013;15(2):106–10.

  39. 39.

    Moss SE, Morgan RO. The annexins. Genome Biol. 2004;5(4):219.

  40. 40.

    Pénzváltó Z, Tegze B, Szász AM, Sztupinszki Z, Likó I, Szendrői A, et al. Identifying resistance mechanisms against five tyrosine kinase inhibitors targeting the ERBB/RAS pathway in 45 cancer cell lines. PLoS One. 2013;8(3), e59503.

  41. 41.

    Tong SW, Yang YX, Hu HD, An X, Ye F, Hu P, et al. Proteomic investigation of 5-fluorouracil resistance in a human hepatocellular carcinoma cell line. J Cell Biochem. 2012;113(5):1671–80.

  42. 42.

    Peraldo-Neia C, Migliardi G, Mello-Grand M, Montemurro F, Segir R, Pignochino Y, et al. Epidermal Growth Factor Receptor (EGFR) mutation analysis, gene expression profiling and EGFR protein expression in primary prostate cancer. BMC Cancer. 2011;11:31.

  43. 43.

    Ruan L, Wang GL, Chen Y, Yi H, Tang CE, Zhang PF, et al. Identification of tyrosine phosphoproteins in signaling pathway triggered TGF-a by using functional proteomics technology. Med Oncol. 2010;27(4):1407–14.

  44. 44.

    Mochinaga K, Tsuchiya T, Nagasaki T, Arai J, Tominaga T, Yamasaki N, et al. High expression of dihydropyrimidine dehydrogenase in lung adenocarcinoma is associated with mutations in epidermal growth factor receptor: implications for the treatment of non--small-cell lung cancer using 5-fluorouracil. Clin Lung Cancer. 2014;15(2):136–44. e134.

  45. 45.

    Ishikawa Y, Kubota T, Otani Y, Watanabe M, Teramoto T, Kumai K, et al. Dihydropyrimidine dehydrogenase and messenger RNA levels in gastric cancer: possible predictor for sensitivity to 5-fluorouracil. Jpn J Cancer Res. 2000;91(1):105–12.

  46. 46.

    Inada T, Ogata Y, Kubota T, Tomikawa M, Yamamoto S, Andoh J, et al. 5-fluorouracil sensitivity and dihydropyrimidine dehydrogenase activity in advanced gastric cancer. Anticancer Res. 2000;20(4):2457–62.

  47. 47.

    Salonga D, Danenberg KD, Johnson M, Metzger R, Groshen S, Tsao-Wei DD, et al. Colorectal tumors responding to 5-fluorouracil have low gene expression levels of dihydropyrimidine dehydrogenase, thymidylate synthase, and thymidine phosphorylase. Clin Cancer Res. 2000;6(4):1322–7.

  48. 48.

    Mizutani Y, Wada H, Fukushima M, Yoshida O, Ukimura O, Kawauchi A, et al. The significance of dihydropyrimidine dehydrogenase (DPD) activity in bladder cancer. Eur J Cancer. 2001;37(5):569–75.

  49. 49.

    Horiguchi J, Takei H, Koibuchi Y, Iijima K, Ninomiya J, Uchida K, et al. Prognostic significance of dihydropyrimidine dehydrogenase expression in breast cancer. Br J Cancer. 2002;86(2):222–5.

  50. 50.

    Zhou H, Jin J, Zhang H, Yi B, Wozniak M, Wong L. IntPath--an integrated pathway gene relationship database for model organisms and important pathogens. BMC Syst Biol. 2012;6 Suppl 2:S2.

  51. 51.

    Yoshida T, Ono H, Kuchiba A, Saeki N, Sakamoto H. Genome-wide germline analyses on cancer susceptibility and GeMDBJ database: Gastric cancer as an example. Cancer Sci. 2010;101(7):1582–9.

Download references

Acknowledgements

We thank Ms. Sumiko Ohnami for help with SNP genotyping. This work was supported in part by the Ministry of Education, Culture, Sports, Science, and Technology of Japan (MEXT): Grants-in-Aid for Scientific Research for Young Scientists (B) (nos. 21710211 and 24710222 to H.T.) and a Grant-in-Aid for Scientific Research on Innovative Areas (no. 26114703 to H.T.). This work was also supported by the Advanced Research for Medical Products Mining Program of the National Institute of Biomedical Innovation (NIBIO ID10-41), the Research Foundation for the Electrotechnology of Chubu, and the Nakajima Foundation.

Author information

Affiliations

  1. Graduate School of Horticulture, Chiba University, 648 Matsudo, Matsudo, Chiba, 271-8510, Japan

    • Hiro Takahashi

  2. Plant Biology Research Center, Chubu University, Matsumoto-cho 1200, Kasugai, Aichi, 487-8501, Japan

    • Hiro Takahashi
    •  & Anna Takahashi

  3. Division of Genetics, National Cancer Center Research Institute, 5-1-1 Tsukiji, Chuo-ku, Tokyo, 104-0045, Japan

    • Hiro Takahashi
    • , Yoko Odaka
    • , Misuzu Okuyama
    • , Hiromi Sakamoto
    •  & Teruhiko Yoshida

  4. Division of Medicinal Safety Science, National Institute of Health Sciences, 1-18-1 Kamiyoga, Setagaya-ku, Tokyo, 158-8501, Japan

    • Nahoko Kaniwa
    • , Yoshiro Saito
    •  & Kimie Sai

  5. Gastrointestinal Medical Oncology Division, National Cancer Center Hospital, 5-1-1 Tsukiji, Chuo-ku, Tokyo, 104-0045, Japan

    • Tetsuya Hamaguchi
    • , Kuniaki Shirao
    •  & Yasuhiro Shimada

  6. Division of Developmental Therapeutics, Research Center for Innovative Oncology, National Cancer Center Hospital East, 6-5-1, Kashiwanoha, Kashiwa, Chiba, 277-8577, Japan

    • Yasuhiro Matsumura

  7. Department of Gastrointestinal Oncology, National Cancer Center Hospital East, 6-5-1, Kashiwanoha, Kashiwa, Chiba, 277-8577, Japan

    • Atsushi Ohtsu
    • , Takayuki Yoshino
    •  & Toshihiko Doi

  8. Division of Functional Biochemistry and Genomics, National Institute of Health Sciences, 1-18-1 Kamiyoga, Setagaya-ku, Tokyo, 158-8501, Japan

    • Jun-ichi Sawada

  9. Present address: Pharmaceutical and Medical Devices Agency, Shinkasumigaseki-building, 3-3-2 Kasumigaseki, Chiyoda-ku, Tokyo, 100-0013, Japan

    • Jun-ichi Sawada

Authors

  1. Search for Hiro Takahashi in:

  2. Search for Nahoko Kaniwa in:

  3. Search for Yoshiro Saito in:

  4. Search for Kimie Sai in:

  5. Search for Tetsuya Hamaguchi in:

  6. Search for Kuniaki Shirao in:

  7. Search for Yasuhiro Shimada in:

  8. Search for Yasuhiro Matsumura in:

  9. Search for Atsushi Ohtsu in:

  10. Search for Takayuki Yoshino in:

  11. Search for Toshihiko Doi in:

  12. Search for Anna Takahashi in:

  13. Search for Yoko Odaka in:

  14. Search for Misuzu Okuyama in:

  15. Search for Jun-ichi Sawada in:

  16. Search for Hiromi Sakamoto in:

  17. Search for Teruhiko Yoshida in:

Corresponding author

Correspondence to Hiro Takahashi.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HT, KS, JS, HS, and TYoshida conceived and designed the experiments. YSaito, KSai, NK, YO, and MO performed the experiments. HT and AT analyzed the data. HT, TH, KShirao, YShimada, YM, AO, TYoshino, and TD contributed reagents/materials/analysis tools. HT wrote the paper. All authors read and approved the final manuscript.

Additional files

Additional file 1: Table S1.

Selected genetic and clinical parameters of gastric cancer patients. (XLS 33 kb)

Additional file 2:

Supplementary Methods. (PDF 91 kb)

Additional file 3:

Supplementary Results. (PDF 125 kb)

Additional file 4: Table S2.

Extracted SNPs with q < 0.95 for the first dataset. (XLS 35 kb)

Additional file 5: Table S3.

Statistical information about the chemotherapeutic response of gastric cancer patients treated with fluoropyrimidine (p < 0.05). (XLS 772 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Takahashi, H., Kaniwa, N., Saito, Y. et al. Construction of possible integrated predictive index based on EGFR and ANXA3 polymorphisms for chemotherapy response in fluoropyrimidine-treated Japanese gastric cancer patients using a bioinformatic method. BMC Cancer 15, 718 (2015) doi:10.1186/s12885-015-1721-z

Download citation

Keywords

  • Single nucleotide polymorphisms
  • Bioinformatics
  • Gastric cancer
  • Genome-wide association study
  • Fluoropyrimidine]

BMC Cancer

ISSN: 1471-2407

Contact us