Biological and Pharmaceutical Bulletin Advance Publication by J-STAGE Advance Publication DOI:10.1248/bpb.b19-00527 October 9, 2019
Biol. Pharm. Bull.
Regular Article
A Genome-Wide Association Study Identifies Five Novel Genetic Markers for
Trastuzumab-Induced Cardiotoxicity in Japanese Population
Mari Hara Nakano,a,b Chihiro Udagawa,a,c Arata Shimo,b Yasuyuki Kojima,b Reiko Yoshie,b
Hisamitsu Zaha,d Norie Abe,d Tokiwa Motonari,d Mikiko Unesoko,d Kenji Tamura,e
Tatsunori Shimoi,e Masayuki Yoshida,f Teruhiko Yoshida,g Hiromi Sakamoto,g Ken Kato,h
Taisei Mushiroda,i Koichiro Tsugawa,b Hitoshi Zembutsu*a,g
aProject for Development of Liquid Biopsy Diagnosis, Japanese Foundation for Cancer
Research, Research Institute; 3-8-31 Ariake, Koto-ku, Tokyo 135-8550, Japan: bDivision of
Breast and Endocrine Surgery, Department of Surgery, St.Marianna University School of
Medicine; 2-16-1 Sugao, Miyamae-ku, Kawasaki 216-8511, Japan: cNew Business
Development Life Science Group, Toyo Kohan Co., Ltd.; 2-18-1, Higashi-Gotanda,
Shinagawa-ku, Tokyo 141-8260, Japan: dDepartment of Breast Surgery, Nakagami Hospital;
610, Noborikawa, Okinawa 904-2195, Japan: eDepartment of Breast and Medical Oncology,
National Cancer Center Hospital; 5-1-1 Tsukiji, Chuo-ku, Tokyo 104-0045, Japan: f
Department of Pathology and Clinical Laboratories, National Cancer Center Hospital; 5-1-1
Tsukiji, Chuo-ku, Tokyo 104-0045, Japan: gFundamental Innovative Oncology Core, National
Cancer Center Research Institute; 5-1-1 Tsukiji, Chuo-ku, Tokyo 104-0045, Japan:
hDepartment of Gastrointestinal Medical Oncology, National Cancer Center Hospital; 5-1-1
Tsukiji, Chuo-ku, Tokyo 104-0045, Japan: and iRIKEN, Center for Integrative Medical
Sciences; 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan.
Ⓒ 2019 The Pharmaceutical Society of Japan
*Corresponding author.
Phone: +81-3-3570-0509 / FAX: +81-3-3570-0515 / E-mail: hitoshi.zembutsu@ jfcr.or.jp
Biological and Pharmaceutical Bulletin Advance Publication
Summary
Trastuzumab has been administered to patients with HER2-positive cancer, however, the cardiotoxicity is identified as one of the life-threatening toxicities.
Clinically useful biomarker for trastuzumab-induced cardiotoxicity has been expected to be developed. To identify a novel genetic marker(s) determining the risk of trastuzumab-induced cardiotoxicity, we performed a first genome-wide association study (GWAS) in Japanese population. We enrolled 481 patients who had been treated with trastuzumab and carried out a GWAS using 11 cases (with cardiotoxicity) and 257 controls (without cardiotoxicity). Top 100 single nucleotide polymorphisms (SNPs) which revealed the smallest P values in GWAS (P = 7.60 x 10-7 - 2.01 x 10-4) were further examined using replication samples consisted of 14 cases and 199 controls.
The combined analysis of the GWAS and replication study indicated possible association of five loci with trastuzumab-induced cardiotoxicity (rs9316695 on chromosome 13q14.3, rs28415722 on chromosome 15q26.3, rs7406710 on chromosome 17q25.3, rs11932853 on chromosome 4q25, and rs8032978 on
-6 -5 -4 -4 chromosome 15q26.3, Pcombined = 6.00 x 10 , 8.88 x 10 , 1.07 x 10 , 1.42 x 10 , 1.60 x
10-4, respectively). Furthermore, we developed a risk prediction model for trastuzumab-induced cardiotoxicity using the five marker SNPs. The incidence of trastuzumab-induced cardiotoxicity in patients with risk score ≥ 5 was significantly higher (42.5%) compared to that in patients with score ≤ 4 (1.8%) (P = 7.82 x 10-15, odds ratio = 40.0). These findings suggest the potential to improve the ability of physicians to avoid the trastuzumab-induced cardiotoxicity for patients with
HER2-positive cancer.
Biological and Pharmaceutical Bulletin Advance Publication
Key words
cardiotoxicity, human epidermal growth factor receptor 2 (HER2), breast cancer, trastuzumab, genome-wide association study (GWAS)
Biological and Pharmaceutical Bulletin Advance Publication
Introduction
Human epidermal growth factor receptor (EGFR) type2 (HER2) is a member of the EGFR family, which is a transmembrane tyrosine kinase receptor involved in cell differentiation and proliferation through the activation of intracellular signaling pathways in response to homo- and hetero-dimerization of HER2 with EGFR family members. It is reported that HER2 is overexpressed or amplified in approximately
25%-30% of breast cancer patients.1-3) Furthermore, patients with HER2-positive cancer have a higher risk of recurrence and mortality than those with HER2-negative cancer.1,4) Trastuzumab is a recombinant humanized monoclonal antibody and selectively targets the extracellular domain of the HER2 protein.1,3-5) Trastuzumab induces antibody-dependent cellular cytotoxicity (ADCC) and prevents HER2 receptor from its dimerization, resulting in inhibiting downstream signaling pathways: the phosphoinositide 3-kinase (PI3K) pathway, AKT pathway and MAPK pathway, which are associated with cell survival and proliferation.3,5,6) It is reported that patients treated with chemotherapy plus first-line trastuzumab showed a significantly better response than chemotherapy alone, and prolonged survival.7-9) However, it is reported that one of the most serious side effects of trastuzumab is cardiotoxicity, which causes cardiac dysfunction, either asymptomatic or symptomatic.6,10,11) Although trastuzumab is an effective drug for HER2-positive cancer, some patients cannot continue the trastuzumab therapy due to the cardiotoxicity. The mechanism of trastuzumab-induced cardiotoxicity is not yet fully understood, however, inhibition of HER2 signaling in cardiomyocytes has been suggested to be involved in this toxicity.12,13)
Anthracycline-induced cardiotoxicity is thought to be cumulative, dose-dependent, and irreversible, whereas trastuzumab-induced cardiotoxicity is mostly reversible and
Biological and Pharmaceutical Bulletin Advance Publication
dose-independent, suggesting that genetic factors play important role in this toxicity.6,14,15)
The association studies of candidate genes including Erb-B2 receptor tyrosine kinase 2 (ERBB2), which encodes HER2, have been reported, however, these results have not yet been sufficiently validated,16-18) and clinically available predictive marker for trastuzumab-induced cardiotoxicity has not yet been developed. Hence, it has been suggested that the other unknown genetic factors which could contribute to the risk of trastuzumab-induced cardiotoxicity should exist. In this study, to identify common genetic variants determining the risk of trastuzumab-induced cardiotoxicity, we carried out a genome-wide association study (GWAS) and replication study of 481 patients treated with trastuzumab, and identified five loci that were likely to be associated with risk of trastuzumab-induced cardiotoxicity.
Materials and Methods
Patients
To identify a genetic marker(s) for the risk of trastuzumab-induced cardiotoxicity, we selected 268 patients treated with trastuzumab from the samples registered in NCC biobank (http://www.ncc.go.jp/jp/biobank/) from February 2010 to
December 2015 for GWAS stage. In the replication study, 213 patients (14 cases and
199 controls) treated with trastuzumab were collected from St. Marianna University
School of Medicine Hospital and Nakagami Hospital from October 2017 to March 2018, and collected from the samples registered in NCC Biobank from January 2016 to
October 2017. We defined trastuzumab-induced cardiotoxicity as left ventricular ejection fraction (LVEF) < 45% or LVEF < 50% with an absolute decrease of 10% from
Biological and Pharmaceutical Bulletin Advance Publication
baseline according to the criteria in Herceptin Adjuvant (HERA) trial.10,12,15,19)
Transthoracic echocardiography was performed for assessment of LVEF at baseline and during follow-up (every 1.0 - 37.7 months, median; 3.1 months). Echocardiography was performed according to standard procedures by one or two sonographers and analyzed by one cardiologist using the iE33 (Philips, Netherlands), SONOS5500
(Philips, Netherlands) or VividE95 (GE, USA) system. All participants provided written informed consent to this study. The research protocol of the study was reviewed and approved by the Research Ethics Committee of the National Cancer
Center (Tokyo, Japan), St. Marianna University School of Medicine (Kawasaki, Japan),
Nakagami Hospital (Okinawa, Japan) and Japanese Foundation for Cancer Research
(Tokyo, Japan), and the members of the committee included experts in the genetic research and ethics.
Genotyping and quality control
Genomic DNA was extracted from peripheral blood or formalin-fixed paraffin-embedded (FFPE) lymph nodes which did not contain cancer cells. For the
GWAS analysis, 268 patients were genotyped using the Infinium OmniExpressExome-8 v1.4 DNA Analysis Kit (Illumina, San Diego, USA). We applied SNP quality control
(call rate of ≥ 99%) in both cases and controls and a Hardy–Weinberg equilibrium (P >
1.0 x 10-6) in controls and SNPs with minor allele frequency of ≥ 0.01. A total of
543,807 SNPs in autosomal chromosomes passed the quality control filters. We utilized the identity-by-state method to evaluate cryptic relatedness for the samples included in this study. The population stratification was examined by principal component analysis (PCA) using the EIGENSTRAT software 6.0.1
Biological and Pharmaceutical Bulletin Advance Publication
(https://www.hsph.harvard.edu/alkes-price/software/). The PCA was carried out by comparing the distribution of the sample populations with three reference populations from the 1,000 Genomes Project Phase3 database that included Europeans (represented by Caucasian from UTAH, CEU), Africans (represented by Yoruba from Ibadan, YRI) and East Asians (represented by Japanese from Tokyo, JPT, Han Chinese from Beijing,
CHB, Southern Han Chinese, CHS, Chinese Dai in Xishuangbanna, CDX, and Kinh in
Ho Chi Minh City, Vietnam, KHV). The PCA was performed on the basis of the genotype information from the samples included in this study. The quantile-quantile
(Q-Q) plot was generated between observed P values of Fisher's exact test against expected P values and showed that genomic inflation factor was 1.172, representing no significant population stratification (Supplementary Fig.1).
Genotype imputation
Genome-wide imputation was performed for the 268 patients used in GWAS.
Because 43 SNPs located in the 3 genomic regions (chromosome 13q14.3 and 2 independent loci on 15q26.3) showed greater significance than marker SNPs (tag-SNPs) at each locus in the genome-wide imputation analysis, these genomic regions
(chromosome 13q14.3, 54593774-54618139, 15q26.3, 98578726-98646496, and
15q26.3, 101796748-101800094, respectively) were also imputed for 213 patients used in replication study. For imputation, we used the reference panel which was based on the 1,000 Genomes Project Phase 3 integrated release version 5 for individuals of East
Asian descent comprising Japanese from Tokyo, Chinese from Beijing and Chinese from southern China. Missing genotypes were imputed with Minimac3 software. We extracted variants with imputation quality threshold of RSQ > 0.3. For imputation
Biological and Pharmaceutical Bulletin Advance Publication
analysis, we performed SNP quality control by excluding SNPs that had low genotyping rate < 99%, and showed deviations from Hardy-Weinberg equilibrium (P ≤ 1.0 x 10-6) in controls and SNPs with minor allele frequency of < 0.01.
Statistical analysis
In the GWAS and the replication study, a case-control study by Fisher's exact test was applied to three genetic models: an allelic frequency model, a dominant-inheritance model, and a recessive-inheritance model. Odds ratios (ORs) and confidence intervals (CIs) were calculated for the genetic model with smallest P value in those three genetic models, using a non-risk allele or a non-risk genotype as a reference. Genome-wide significance levels after adjustment of Bonferroni correction for multiple testing of three genetic models were P = 3.06 x 10-8 [0.05 / (543,807 x 3)] in the GWAS and P = 1.67 x 10-4 [0.05 / (100 x 3)] in the replication study. For combination analysis, the genotype count of the replication study was added to that of the GWAS. The differences in the distribution of age, primary cancer site, pretreatment with anthracycline, HER2 status and hormone receptor status were evaluated by logistic regression analysis (Supplementary Table 1). For the predictive scoring system of the risk of trastuzumab-induced cardiotoxicity, we assigned a score of
2 to individuals homozygous for risk allele, 1 to individuals heterozygous for risk allele, and 0 to individuals homozygous for non-risk allele at rs9316695, rs11932853 and rs8032978, respectively. We further assigned a score of 2 to individuals homozygous for the risk allele and 0 to individuals with the other genotypes at rs28415722 and rs7406710 and summed up the scores for each SNP to obtain individuals’ scores.
Based on this system, each patient was classified into any of the nine prediction score
Biological and Pharmaceutical Bulletin Advance Publication
groups (score 0, 1, 2, 3, 4, 5, 6, 7, or 8). All the statistical analyses were carried out using R statistical environment version 3.3.1 (http://www.r-project.org/), PLINK version
1.07,20) or the BellCurve for Excel (Social Survey Research Information Co., Ltd.,
Tokyo, Japan). Regional association plots were generated using Locus Zoom
(http://locuszoom.sph.umich.edu/).
Results
Patient characteristics
To identify a genetic marker(s) determining the risk of trastuzumab-induced cardiotoxicity, we recruited 481 patients treated with trastuzumab, including 25 cases
(with trastuzumab-induced cardiotoxicity) and 456 controls (without trastuzumab-induced cardiotoxicity). Table 1 showed the characteristics of these 481 patients receiving trastuzumab treatment. Their median age at the beginning of the trastuzumab treatment was 56 years old in both case and control groups. The distributions of gender (percentage of female patients) were 100 and 76.7% in the case and control groups, respectively, in the GWAS, and 100 and 91.5%, respectively, in the replication samples (Supplementary Table 2). The incidence of cardiotoxicity in gastric and breast cancers were 1.1% (1/90) and 6.2% (24/387), respectively (Table 1).
Among the characteristics listed in Table 1, none of them showed significant association with trastuzumab-induced cardiotoxicity in logistic regression analysis (Supplementary
Table 1).
GWAS and replication analysis
In this study, trastuzumab-induced cardiotoxicity was defined as LVEF <45%
Biological and Pharmaceutical Bulletin Advance Publication
or LVEF <50% with an absolute decrease of 10% from baseline according to the criteria in Herceptin Adjuvant (HERA) trial. 10, 12, 15, 19) We carried out a GWAS of 268 patients who received trastuzumab therapy using the Infinium OmniExpressExome-8 v1.4 DNA Analysis Kit (Illumina). After the standard quality control, an association analysis was carried out for 543,807 SNPs by Fisher’s exact test on the basis of three genetic models: allelic, dominant, and recessive. We generated a quantile-quantile plot and observed the genomic inflation factor of 1.172, indicating no population stratification (Supplementary Fig.1). However, we could not observe the SNP which reached the genome-wide significance level (Supplementary Fig.2). The top 100 SNPs which revealed the smallest P value in GWAS showed possible associations (P =7.60 x
10-7 – 2.01 x 10-4, Supplementary Table 3). To further validate the results of the
GWAS analysis, we carried out a replication study of the 100 SNPs using 213 independent patients, including 14 cases and 199 controls. Although no SNP reached the significance level for three genetic models after the Bonferroni correction (P = 1.67 x 10-4 [0.05/ (100 x 3)]) in the replication study, we identified 16 SNPs with suggestive associations (P = 3.57 x 10-2 - 9.37 x 10-2, Supplementary Table 4). A combined analysis of the GWAS and the replication study for the 16 SNPs indicated that nine
SNPs in five loci showed stronger association (smaller P value) than those in GWAS results (rs9316695, rs9527156, rs12583122 and rs11617903 on chromosome 13q14.3,
-6 -6 -5 -5 Pcombined = 6.00 x 10 , 8.64 x 10 , 1.92 x 10 and 2.22 x 10 , respectively, rs1383149
-4 -4 and rs12372962 on chromosome 15q26.3, Pcombined = 1.01 x 10 and 1.15 x 10 ,
-4 respectively, rs7406710 on chromosome 17q25.3, Pcombined = 1.07 x 10 , rs11932853 on
-4 chromosome 4q25, Pcombined =1.42 x 10 , and rs8032978 on chromosome 15q26.3,
-4 Pcombined = 1.60 x 10 , Supplementary Table 5, Fig. 1), however, none of them reached
Biological and Pharmaceutical Bulletin Advance Publication
the genome-wide significance level in this study (P = 3.06 x 10-8).
Furthermore, we carried out subgroup analysis according to cancer type and hormone receptor status (Supplementary Table 6), and found that rs8032978 on chromosome 15q26.3 showed stronger association in patients with hormone receptor positive breast cancer compared to that in original combined analysis (P = 6.55 x 10-5,
OR=8.53, 95% CI=2.91-23.39, Supplementary Table 6). This finding suggests that rs8032978 might also be a promising marker of trastuzumab-induced cardiotoxicity in patients with hormone receptor positive breast cancer. Association analysis was also carried out by logistic regression analysis with adjusting medication for cardiovascular disease and gender which could be associated with the risk of cardiotoxicity
(Supplementary Table 2), however, none of the above SNPs did not reach genome-wide significance level (Supplementary Table 7).
Imputation analysis
To identify additional candidate loci associated with trastuzumab-induced cardiotoxicity, we examined the associations by using genome-wide imputed genotypes of GWAS samples, however, we could not observe novel candidate locus associated with the trastuzumab-induced cardiotoxicity. Of the five candidate loci identified in
GWAS (chromosome 4q25, 13q14.3, 2 independent loci on 15q26.3, and 17q25.3), imputation analysis identified 43 novel SNPs in three loci (chromosome 13q14.3 and 2 independent loci on 15q26.3) of which significance were equal to or greater than the marker SNPs (tag-SNPs) identified in GWAS (Fig. 1A, B, E, Supplementary Table 8).
We further performed a replication study using 213 independent patients for the 43
SNPs and found that no SNP showed stronger association than each marker SNP in the
Biological and Pharmaceutical Bulletin Advance Publication
above three loci (Supplementary Table 8). However, a combined result of the GWAS and replication stage for the 43 SNPs revealed that the association of an imputed SNP, rs28415722 on chromosome 15q26.3, was the strongest in the region (Pcombined = 8.88 x
10−5, Supplementary Table 5, 8). Hence, we regarded rs28415722 as a marker SNP in the region (Table 2, Fig. 1B)
Predictive scoring system for trastuzumab-induced cardiotoxicity using the five marker SNPs
Because the five SNPs, which showed smallest P value in the combined study at each candidate locus (rs9316695, rs28415722, rs7406710, rs11932853 and rs8032978,
Table 2), were independent predictors of trastuzumab-induced cardiotoxicity by logistic regression analysis (P = 2.82 x 10-4 - 4.15 x 10-3, Supplementary Table 1), we investigated combined effects of the five loci on the risk of cardiotoxicity in patients treated with trastuzumab by using a scoring system. For the predictive scoring system, each patient was scored considering the total number of risk alleles and genotypes. We gave a score of 2 to individuals homozygous for risk allele and 0 to individuals with the other genotypes at rs28415722 and rs7406710 because these SNPs showed smallest P value in recessive-inheritance model. We further gave a score of 2 to individuals homozygous for risk allele, 1 to individuals heterozygous for risk allele and 0 to individuals homozygous for non-risk allele at rs9316695, rs11932853 and rs8032978.
The patients’ scores were obtained by summing up the scores of the five SNPs. The patients used in this study were classified into nine prediction score groups (score 0-8,
Table 3). The proportion of patients with trastuzumab-induced cardiotoxicity was likely to be increased in groups with higher prediction scores; the incidences of the
Biological and Pharmaceutical Bulletin Advance Publication
cardiotoxicity were 1.8% (8/441) in the score 0-4 group, 36.4% (8/22) in the score 5 group, 22.2% (2/9) in the score 6 group, 75.0% (3/4) in the score 7 group and 80.0%
(4/5) in the score 8 group (Table 3). We further classify the patients into two groups; score 0-4 group and 5-8 group, considering the best threshold for the diagnostic test
(sensitivity, 68.0% and specificity, 95.0%). The incidences of trastuzumab-induced cardiotoxicity in score 5-8 group were 45.8% (11/24), 37.5% (6/16) and 42.5% (17/40), whereas those in score 0-4 group were 0% (0/244), 4.1% (8/197) and 1.8% (8/441) in the GWAS, replication and combined stage, respectively, suggesting that the predictive scoring system using the above five SNPs could predict the risk of cardiotoxicity prior
-15 to initiation of trastuzumab therapy (Pcombined = 7.82 x 10 , OR=40.0, 95%
CI=15.6-102.3, Table 3 , Fig. 2). Although we investigated the relationship between the prediction score and cardiotoxicity grade (LVEF) or time-to-event (cardiotoxicity), we could not observe significant relationship between them (Supplementary Fig.3).
Moreover, we performed risk stratification analysis based on medication for cardiovascular disease because it showed significant association with cardiotoxicity (P
= 1.80 x 10-4, Supplementary Table 2). However, this predictive scoring system appeared to be medication-nonspecific (independent of medication) because the associations of cardiotoxicity with prediction scores were comparable between medicated and unmedicated patients (medicated; P = 6.82 x 10-6, unmedicated; P = 1.01 x 10-6, Supplementary Table 9).
Discussion
The safe and effective treatment with trastuzumab based on individual germline and/or somatic genetic information is important to provide personalized medicine for
Biological and Pharmaceutical Bulletin Advance Publication
patients with HER2-positive cancer. Although a GWAS of cardiotoxicity in Adjuvant
Trastuzumab Trial (NCCTG N9831) in Caucasian population has already been reported,21) this study is the first GWAS which attempts to identify common genetic variants associated with trastuzumab-induced cardiotoxicity in Japanese population.
The genome-wide significant SNP was not identified in this study, however, we identified five independent loci showing possible association with trastuzumab-induced cardiotoxicity (Table 2, Fig. 1). Interestingly, of the above loci, proprotein convertase subtilisin/kexin type 6 (PCSK6) and small nuclear ribonucleoprotein polypeptide A
(SNRPA1), which have been suggested to promote the proliferation of breast cancer and colon cancer, respectively, were included in the genomic region of chromosome
15q26.3 (Fig. 1E). 22,23) We observed significantly higher incidence of trastuzumab-induced cardiotoxicity in patients with higher score (≥ 5) compared to those with lower score (≤ 4) in predictive scoring system using the five marker SNPs
(Fig. 2, Table 2, 3). This result suggests the potential to improve the ability of physicians to avoid the risk of trastuzumab-induced cardiotoxicity based on the result of this predictive scoring system for the treatment of patients with HER2-positive cancer.
-6 rs9316695, the marker SNP on chromosome 13q14.3 (Pcombined = 6.00 x 10 ,
OR= 4.46, Table 2), which showed the strongest association in a combined analysis of
GWAS and replication study, was located in intergenic region, and long intergenic non-protein coding RNA 558 (LINC00558) was included in this region (Fig. 1A).
NOTCH3, which encodes the notch receptor 3, has been reported to be one of the target genes of LINC00558,24) and NOTCH3 signaling pathways are significantly activated in
ErbB2-negative cells.25,26) The Notch signaling system is present in the cardiomyocytes and suggested to play a role in myocardial remodeling in heart
Biological and Pharmaceutical Bulletin Advance Publication
failure.27,28) Hence, genetic variation in LINC00558 might contribute to interindividual variability in trastuzumab-induced cardiotoxicity through effects on Notch signaling in the stressed (inactivation of HER2 signaling) cardiomyocyte, although further validation studies will be needed to prove the true association.
rs8032978 on chromosome 15q26.3, which is one of the five marker SNPs in
-4 this study (Pcombined = 1.60 x 10 , OR= 5.83, Table 2), was located 44 kb downstream from the PCSK6 gene (Fig. 1E). PCSK6 activates corin which is a type-II transmembrane serine protease found mainly in the heart and converts it to a two-chain active enzyme through cleavage at Arg-801.29-32) Corin is known to play an essential role in the regulation of water and salt balance by converting natriuretic peptides to their active form,30,33) which promotes vasodilation, diuresis and natriuresis, and prevent cardiac remodeling, in patients with heart failure.29,33) These lines of evidence suggest that impaired activity and/or gene expression of PCSK6 might contribute to the inability to recover the cardiac function in damaged (inactivation of HER2 signaling) heart due to insufficient activation of natriuretic peptides.29-33)
To effectively clarify the mechanism of trastuzumab-induced cardiotoxicity, we performed expression quantitative trait locus (cis-eQTL) analysis for the five candidate loci. We observed cis-regulatory effect of rs7406710, which is located 1.5kb downstream from fascin actin-bundling protein 2, retinal (FSCN2) on chromosome
17q25.3 (Fig. 1C), on gene expression level of FSCN2 (P = 3.20 x 10-10, Supplementary
Fig.4). FSCN2 is actin-binding proteins that cross-link filamentous actin into tightly packed parallel bundles.34,35) FSCN2 can bind with multiple actins,34,36) and has been reported to bind with beta-, gamma-, and alpha-actin which is highly expressed in cardiomyocyte.34,37,38) Although the relationship between HER2 signaling and actin
Biological and Pharmaceutical Bulletin Advance Publication
crosslinking protein, FSCN2, in myocardium is not yet clarified, lower expression of
FSCN2 in cardiomyocyte might impair cardiac contraction under the condition of HER2 inhibition (Supplementary Fig.4). Further studies will be required to clarify the functional associations of FSCN2 with cardiotoxicity in patients treated with trastuzumab.
The candidate gene association studies previously reported that rs1058808
(P1170A) and rs1136201 (I655V) in ERBB2 (HER2) could be candidates that confer susceptibility to trastuzumab-induced cardiotoxicity.16-18) In our GWAS, the two SNPs showed no association with trastuzumab-induced cardiotoxicity (rs1058808; P = 1.21 x
10-1, rs1136201; P = 6.86 x 10-1). Moreover, a first GWAS of trastuzumab-induced cardiotoxicity in Caucasian population suggested that SNPs at six loci, LIM domain binding 2 (LDB2), BMP/retinoic acid inducible neural specific 1 (BRINP1), chr6 intergenic, RAB22A, member RAS oncogene family (RAB22A), transient receptor potential cation channel subfamily C member 6 (TRPC6) and long intergenic non-protein coding RNA 1060 (LINC01060), were associated with a decline in LVEF (P
= 8.93 x 10-8 - 7.73 x 10-6) in patients who received chemotherapy plus trastuzumab.21)
In our GWAS, the SNPs in these candidate genes showed no or weak association with trastuzumab-induced cardiotoxicity (P = 3.33 x 10-3 - 9.99 x 10-1). Although the small sample size in our study might be one of the causes of the discrepant results, interethnic difference in allele frequency of the above SNPs could also affect the discrepancy.
Moreover, mechanism of action involved in the trastuzumab-induced cardiotoxicity might differ between ethnic groups. The P values of the SNPs in ATP binding cassette subfamily B member 1 (ABCB1), solute carrier family 28 member 3 (SLC28A3), carbonyl reductase 3 (CBR3), neutrophil cytosolic factor (NCF4), Rac family small
Biological and Pharmaceutical Bulletin Advance Publication
GTPase 2 (RAC2), retinoic acid receptor gamma (RARG) and UDP glucuronosyltransferase family 1 member A6 (UGT1A6), which were reported to be associated with anthracycline-induced cardiotoxicity,39-42) were from 1.68 x 10-2 to 9.99 x 10-1, indicating no or weak association in our GWAS. The sample size used in our study might not be enough to detect true associations of the above SNPs because the effect sizes of these candidate genes were not large.16,17,21,39-42)
In conclusion, our GWAS using 481 patients treated with trastuzumab identified five novel candidate loci, chromosome 4q25, 13q14.3, 17q25.3, and 2 independent loci on 15q26.3, which were associated with trastuzumab-induced cardiotoxicity. Furthermore, the combined analysis of the five SNP loci revealed that the predictive score based on risk genotypes was significantly associated with trastuzumab-induced cardiotoxicity. These findings provide new insights into personalized trastuzumab therapy for patients with HER2 positive cancer. However, larger-scale replication study and further functional analysis are required to verify our results and to clarify their biological mechanisms which have effects on trastuzumab-induced cardiotoxicity.
Acknowledgments
We express our heartfelt gratitude to all the study participants. We thank Mr.
Satoshi Kasahara, Ms. Aya Imai and Ms. Emi Shigino for their contribution to the collection of clinical information, and all staffs in Nakagami Hospital, St. Marianna
University School of Medicine, and National Cancer Center Biobank for their contribution to the sample collection.
Biological and Pharmaceutical Bulletin Advance Publication
Conflict of Interest
The authors declare no conflict of interest.
Supplementary Materials
The online version of this article contains supplementary materials.
Biological and Pharmaceutical Bulletin Advance Publication
References
1) Bonifazi M, Franchi M, Rossi M, Moja L, Zambelli A, Zambon A, T Corrao G,
La Vecchia C, Zocchetti C, Negri E. Trastuzumab-related cardiotoxicity in early breast cancer: a cohort study. Oncologist., 18, 795-801 (2013).
2) Hsieh AC, Moasser MM. Targeting HER proteins in cancer therapy and the role of the non-target HER3. Br J Cancer., 97, 453-457 (2007).
3) Valabrega G, Montemurro F, Aglietta M. Trastuzumab: mechanism of action, resistance and future perspectives in HER2-overexpressing breast cancer. Ann Oncol.,
18, 977-984 (2007).
4) Ponde NF, Lambertini M, de Azambuja E. Twenty years of anti-HER2 therapy-associated cardiotoxicity. ESMO Open., 1, e000073 (2016).
5) Jones KL, Buzdar AU. Evolving novel anti-HER2 strategies. Lancet Oncol., 10,
1179-1187 (2009).
6) Chen ZI, Ai DI. Cardiotoxicity associated with targeted cancer therapies. Mol
Clin Oncol., 4, 675-681 (2016).
7) Jeyakumar A, Younis T. Trastuzumab for HER2-Positive Metastatic Breast
Cancer: Clinical and Economic Considerations. Clin Med Insights Oncol., 6, 179-187
(2012).
8) Slamon DJ, Leyland-Jones B, Shak S, Fuchs H, Paton V, Bajamonde A,
Fleming T, Eiermann W, Wolter J, Pegram M, Baselga J, Norton L. Use of chemotherapy plus a monoclonal antibody against HER2 for metastatic breast cancer that overexpresses HER2. N Engl J Med., 344, 783-792 (2001).
9) Montemurro F, Donadio M, Clavarezza M, Redana S, Jacomuzzi ME,
Valabrega G, Danese S, Vietti-Ramus G, Durando A, Venturini M, Aglietta M. Outcome
Biological and Pharmaceutical Bulletin Advance Publication
of patients with HER2-positive advanced breast cancer progressing during trastuzumab-based therapy. Oncologist., 11, 318-324 (2006).
10) Sengupta PP, Northfelt DW, Gentile F, Zamorano JL, Khandheria BK.
Trastuzumab-induced cardiotoxicity: heart failure at the crossroads. Mayo Clin Proc.,
83, 197-203 (2008).
11) Advani PP, Ballman KV, Dockter TJ, Colon-Otero G, Perez EA. Long-Term
Cardiac Safety Analysis of NCCTG N9831 (Alliance) Adjuvant Trastuzumab Trial. J
Clin Oncol., 34, 581-587 (2016).
12) Dias A, Claudino W, Sinha R, Perez CA, Jain D. Human epidermal growth factor antagonists and cardiotoxicity-A short review of the problem and preventative measures. Crit Rev Oncol Hematol., 104, 42-51 (2016).
13) Florido R, Smith KL, Cuomo KK, Russell SD. Cardiotoxicity From Human
Epidermal Growth Factor Receptor-2 (HER2) Targeted Therapies. J Am Heart Assoc., 6, e006915 (2017).
14) Wilkinson EL, Sidaway JE, Cross MJ. Cardiotoxic drugs Herceptin and doxorubicin inhibit cardiac microvascular endothelial cell barrier formation resulting in increased drug permeability. Biol Open., 5, 1362-1370 (2016).
15) Bloom MW, Hamo CE, Cardinale D, Ky B, Nohria A, Baer L, Skopicki H,
Lenihan DJ, Gheorghiade M, Lyon AR, Butler J. Cancer Therapy-Related Cardiac
Dysfunction and Heart Failure: Part 1: Definitions, Pathophysiology, Risk Factors, and
Imaging. Circ Heart Fail., 9, e002661 (2016).
16) Stanton SE, Ward MM, Christos P, Sanford R, Lam C, Cobham MV, Donovan
D, Scheff RJ, Cigler T, Moore A, Vahdat LT, Lane ME, Chuang E. Pro1170 Ala polymorphism in HER2-neu is associated with risk of trastuzumab cardiotoxicity. BMC
Biological and Pharmaceutical Bulletin Advance Publication
Cancer., 15, 267 (2015).
17) Beauclair S, Formento P, Fischel JL, Lescaut W, Largillier R, Chamorey E,
Hofman P, Ferrero JM, Pagès G, Milano G. Role of the HER2 [Ile655Val] genetic polymorphism in tumorogenesis and in the risk of trastuzumab-related cardiotoxicity.
Ann Oncol., 18, 1335-1341 (2007).
18) Lemieux J, Diorio C, Cote MA, Provencher L, Barabe F, Jacob S, St-Pierre C,
Demers E, Tremblay-Lemay R, Nadeau-Larochelle C, Michaud A, Laflamme C.
Alcohol and HER2 polymorphisms as risk factor for cardiotoxicity in breast cancer treated with trastuzumab. Anticancer Res., 33, 2569-2576 (2013).
19) Cameron D, Piccart-Gebhart MJ, Gelber RD, Procter M, Goldhirsch A, de
Azambuja E, Castro G Jr, Untch M, Smith I, Gianni L, Baselga J, Al-Sakaff N, Lauer S,
McFadden E, Leyland-Jones B, Bell R, Dowsett M, Jackisch C. Herceptin Adjuvant
(HERA) Trial Study Team. 11 years' follow-up of trastuzumab after adjuvant chemotherapy in HER2-positive early breast cancer: final analysis of the HERceptin
Adjuvant (HERA) trial., Lancet., 389, 1195-1205 (2017).
20) Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller
J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet., 81, 559-575
(2007).
21) Serie DJ, Crook JE, Necela BM, Dockter TJ, Wang X, Asmann YW,
Fairweather D, Bruno KA, Colon-Otero G, Perez EA, Thompson EA, Norton N.
Genome-wide association study of cardiotoxicity in the NCCTG N9831 (Alliance) adjuvant trastuzumab trial. Pharmacogenet Genomics., 27, 378-385 (2017).
22) Wang P, Wang F, Wang L, Pan J. Proprotein convertase subtilisin/kexin type 6
Biological and Pharmaceutical Bulletin Advance Publication
activates the extracellular signal-regulated kinase 1/2 and Wnt family member 3A pathways and promotes in vitro proliferation, migration and invasion of breast cancer
MDA-MB-231 cells. Oncol Lett., 16,145-150 (2018).
23) Zeng Q, Lei F, Chang Y, Gao Z, Wang Y, Gao Q, Niu P, Li Q. An oncogenic gene, SNRPA1, regulates PIK3R1, VEGFC, MKI67, CDK1 and other genes in colorectal cancer. Biomed Pharmacother.,117, 109076 (2019).
24) Chen X, Jung JG, Shajahan-Haq AN, Clarke R, Shih Ie M, Wang Y, Magnani L,
Wang TL, Xuan J. ChIP-BIT: Bayesian inference of target genes using a novel joint probabilistic model of ChIP-seq profiles. Nucleic Acids Res., 44, e65 (2016).
25) Hirose H, Ishii H, Mimori K, Ohta D, Ohkuma M, Tsujii H, Saito T, Sekimoto
M, Doki Y, Mori M. Notch pathway as candidate therapeutic target in Her2/Neu/ErbB2 receptor-negative breast tumors. Oncol Rep., 23, 35-43 (2010).
26) Yamaguchi N, Oyama T, Ito E, Satoh H, Azuma S, Hayashi M, Shimizu K,
Honma R, Yanagisawa Y, Nishikawa A, Kawamura M, Imai J, Ohwada S, Tatsuta K,
Inoue J, Semba K, Watanabe S. NOTCH3 signaling pathway plays crucial roles in the proliferation of ErbB2-negative human breast cancer cells. Cancer Res., 68, 1881-1888
(2008).
27) Oie E, Sandberg WJ, Ahmed MS, Yndestad A, Laerum OD, Attramadal H,
Aukrust P, Eiken HG. Activation of Notch signaling in cardiomyocytes during post-infarction remodeling. Scand Cardiovasc J., 44, 359-366 (2010).
28) Rizzo P, Mele D, Caliceti C, Pannella M, Fortini C, Clementz AG, Morelli MB,
Aquila G, Ameri P, Ferrari R. The role of notch in the cardiovascular system: potential adverse effects of investigational notch inhibitors. Front Oncol., 4, 384 (2014).
29) Khoury EE, Kinaneh S, Aronson D, Amir O, Ghanim D, Volinsky N, Azzam Z,
Biological and Pharmaceutical Bulletin Advance Publication
Abassi Z. Natriuretic peptides system in the pulmonary tissue of rats with heart failure: potential involvement in lung edema and inflammation. Oncotarget., 9, 21715-21730
(2018).
30) Yan W, Wu F, Morser J, Wu Q. Corin, a transmembrane cardiac serine protease, acts as a pro-atrial natriuretic peptide-converting enzyme. Proc Natl Acad Sci U S A., 97,
8525-8529 (2000).
31) Chen S, Cao P, Dong N, Peng J, Zhang C, Wang H, Zhou T, Yang J, Zhang Y,
Martelli EE, Naga Prasad SV, Miller RE, Malfait AM, Zhou Y, Wu Q. PCSK6-mediated corin activation is essential for normal blood pressure. Nat Med., 21, 1048-1053 (2015).
32) Yan W, Sheng N, Seto M, Morser J, Wu Q. Corin, a mosaic transmembrane serine protease encoded by a novel cDNA from human heart. J Biol Chem., 274,
14926-14935 (1999).
33) Wu F, Yan W, Pan J, Morser J, Wu Q. Processing of pro-atrial natriuretic peptide by corin in cardiac myocytes. J Biol Chem., 277, 16900-16905 (2002).
34) Liu X, Zhao M, Xie Y, Li P, Wang O, Zhou B, Yang L, Nie Y, Cheng L, Song X,
Jin C, Han F. Null Mutation of the Fascin2 Gene by TALEN Leading to Progressive
Hearing Loss and Retinal Degeneration in C57BL/6J Mice. G3 (Bethesda)., 8,
3221-3230 (2018).
35) Hashimoto Y, Kim DJ, Adams JC. The roles of fascins in health and disease. J
Pathol., 224, 289-300 (2011).
36) Avenarius MR, Saylor KW, Lundeberg MR, Wilmarth PA, Shin JB, Spinelli KJ,
Pagana JM, Andrade L, Kachar B, Choi D, David LL, Barr-Gillespie PG. Correlation of actin crosslinker and capper expression levels with stereocilia growth phases. Mol Cell
Proteomics., 13, 606-620 (2014).
Biological and Pharmaceutical Bulletin Advance Publication
37) Perrin BJ, Strandjord DM, Narayanan P, Henderson DM, Johnson KR, Ervasti
JM. beta-Actin and fascin-2 cooperate to maintain stereocilia length. J Neurosci., 33,
8114-8121 (2013).
38) Shin JB, Longo-Guess CM, Gagnon LH, Saylor KW, Dumont RA, Spinelli KJ,
Pagana JM, Wilmarth PA, David LL, Gillespie PG, Johnson KR. The R109H variant of fascin-2, a developmentally regulated actin crosslinker in hair-cell stereocilia, underlies early-onset hearing loss of DBA/2J mice. J Neurosci., 30, 9683-9694 (2010).
39) Hertz DL, Caram MV, Kidwell KM, Thibert JN, Gersch C, Seewald NJ,
Smerage J, Rubenfire M, Henry NL, Cooney KA, Leja M, Griggs JJ, Rae JM. Evidence for association of SNPs in ABCB1 and CBR3, but not RAC2, NCF4, SLC28A3 or
TOP2B, with chronic cardiotoxicity in a cohort of breast cancer patients treated with anthracyclines. Pharmacogenomics., 17, 231-240 (2016).
40) Aminkeng F, Bhavsar AP, Visscher H, Rassekh SR, Li Y, Lee JW, Brunham LR,
Caron HN, van Dalen EC, Kremer LC, van der Pal HJ, Amstutz U, Rieder MJ,
Bernstein D, Carleton BC, Hayden MR, Ross CJ; Canadian Pharmacogenomics
Network for Drug Safety Consortium. A coding variant in RARG confers susceptibility to anthracycline-induced cardiotoxicity in childhood cancer. Nat Genet., 47, 1079-1084
(2015).
41) Aminkeng F, Ross CJ, Rassekh SR, Hwang S, Rieder MJ, Bhavsar AP, Smith A,
Sanatani S, Gelmon KA, Bernstein D, Hayden MR, Amstutz U, Carleton BC; CPNDS
Clinical Practice Recommendations Group. Recommendations for genetic testing to reduce the incidence of anthracycline-induced cardiotoxicity. Br J Clin Pharmacol., 82,
683-695 (2016).
42) Visscher H, Ross CJ, Rassekh SR, Barhdadi A, Dube MP, Al-Saloos H, Sandor
Biological and Pharmaceutical Bulletin Advance Publication
GS, Caron HN, van Dalen EC, Kremer LC, van der Pal HJ, Brown AM, Rogers PC,
Phillips MS, Rieder MJ, Carleton BC, Hayden MR; Canadian Pharmacogenomics
Network for Drug Safety Consortium. Pharmacogenomic prediction of anthracycline-induced cardiotoxicity in children. J Clin Oncol., 30, 1422-1428 (2012)
Biological and Pharmaceutical Bulletin Advance Publication
Table 1. Patient demographics and clinical characteristics.
ER, estrogen receptor; PR, progesterone receptor.
Biological and Pharmaceutical Bulletin Advance Publication
Table 2. Summary of association results of GWAS and replication study for five loci.
Case Control P value Chromosomal Alleleb Risk Genotyped/ Positiona SNP Gene Stage Genotype Genotype Odds ratio (95% CI)c region [1/2] allele Imputed RAF RAF Allelic Dominant Recessive 11 12 22 11 12 22 GWAS 3 5 3 0.50 199 55 3 0.12 2.59 x 10-5 8.22 x 10-4 9.71 x 10-4 7.38 (2.77-19.66) 13q14.3 54,615,773 rs9316695 No gene C/A A replication Genotyped 7 6 1 0.29 153 43 3 0.12 3.79 x 10-2 4.80 x 10-2 2.40 x 10-1 2.84 (1.02-7.19) combined 10 11 4 0.38 352 98 6 0.12 6.00 x 10-6 1.22 x 10-4 9.73 x 10-4 4.46 (2.30-8.47) GWAS 1 2 8 0.82 95 119 43 0.40 1.97 x 10-4 1.04 x 10-1 1.10 x 10-4 13.08(2.99-79.58) d 15q26.3 98,609,746 rs28415722 No gene G/A A replication Imputed 6 3 5 0.46 82 85 32 0.37 4.21 x 10-1 1.00 7.31 x 10-2 2.88(0.71-10.34) combined 7 5 13 0.62 177 204 75 0.39 1.62 x 10-3 3.98 x 10-1 8.88 x 10-5 5.48(2.21-13.69) GWAS 11 0 0 1.00 111 119 27 0.66 2.38 x 10-4 6.10 x 10-1 1.35 x 10-4 NA 17q25.3 79,505,624 rs7406710 No gene C/T C replication Genotyped 10 3 1 0.82 90 93 16 0.69 2.01 x 10-1 1.00 9.37 x 10-2 3.01 (0.83-13.61) combined 21 3 1 0.90 201 212 43 0.67 4.55 x 10-4 7.17 x 10-1 1.07 x 10-4 6.64 (2.19-27.01) GWAS 7 3 1 0.77 31 129 97 0.37 2.26 x 10-4 6.04 x 10-2 1.47 x 10-4 12.54 (2.99-61.88) 4q25 112,024,388 rs11932853 No gene T/C T replication Genotyped 4 8 2 0.57 27 99 73 0.38 7.05 x 10-2 1.45 x 10-1 1.28 x 10-1 2.13 (0.92-5.08) combined 11 11 3 0.66 58 228 170 0.38 1.42 x 10-4 9.56 x 10-3 2.10 x 10-4 3.20 (1.70-6.23) GWAS 6 4 1 0.27 239 18 0 0.04 1.95 x 10-4 9.84 x 10-4 4.10 x 10-2 10.22 (2.93-31.90) 15q26.3 101,799,790 rs8032978 No gene A/G G replication Genotyped 11 3 0 0.11 185 13 1 0.04 1.07 x 10-1 8.87 x 10-2 1.00 3.57 (0.57-15.86) combined 17 7 1 0.18 424 31 1 0.04 1.60 x 10-4 4.29 x 10-4 1.01 x 10-1 5.83 (2.30-13.51) a Based on GRCh37 genome assembly. b Reference allele (GRCh37) was defined as allele 1. c Odds ratios were shown for the model with minimum P values. d Genotype of rs28415722 was validated using TaqMan SNP genotyping assay.
GWAS, genome-wide association study; SNP, single nucleotide polymorphism; RAF, risk allele frequency; CI, confidence interval; NA, not available.
Biological and Pharmaceutical Bulletin Advance Publication
Table 3. Prediction scores of trastuzumab-induced cardiotoxicity using rs9316695, rs28415722, rs7406710, rs11932853 and rs8032978.
N (%) Score % of Case P value Odds ratio (95% CI) Case Control 0 0 (0.0) 53 (11.6) 1 0 (0.0) 93 (20.4) 2 3 (12.0) 109 (23.9) 1.8 1.00 (reference) 3 2 (8.0) 111 (24.3) 4 3 (12.0) 67 (14.7) P = 7.82 x 10-15 5 8 (32.0) 14 (3.1) 36.4 6.68 x 10-8 30.9 (10.1-94.3) Odds ratio = 40.0 6 2 (8.0) 7 (1.5) 22.2 1.47 x 10-2 15.5 (2.8-86.4） (95%CI, 15.6-102.3） 7 3 (12.0) 1 (0.2) 75.0 4.46 x 10-5 162.4 (15.2-1734.8) 8 4 (16.0) 1 (0.2) 80.0 1.50 x 10-6 216.5 (21.7-2159.8)
CI, confidence interval.
Biological and Pharmaceutical Bulletin Advance Publication
Biological and Pharmaceutical Bulletin Advance Publication
Fig. 1. Regional association plots of five loci associated with trastuzumab-induced cardiotoxicity. Upper panel; P values of genotyped SNPs (circle) and imputed SNPs (square) are plotted (as -log10P value) against their physical location on chromosome 13q14.3(A), 15q26.3(B), 17q25.3(C), 4q25(D) and 15q26.3(E). The genetic recombination rates estimated from 1,000 Genomes samples (JPT + CHB) are shown with a blue line.
Biological and Pharmaceutical Bulletin Advance Publication
SNP's color indicates LD with rs9316695 (A), rs28415722 (B), rs7406710 (C), rs11932853 (D) and rs8032978 (E) according to a scale from r2 = 0 to 1 based on pair-wise r2 values from ASN data in 1,000 Genomes Project. Lower panel; gene annotations from the University of California Santa Cruz genome browser.
Biological and Pharmaceutical Bulletin Advance Publication
Fig. 2. Predictive risk score and incidence of trastuzumab-induced cardiotoxicity.
The incidence of trastuzumab-induced cardiotoxicity in patients with score ≥ 5 was significantly higher compared to those in patients with score ≤ 4 (Fisher’s exact test; P
= 7.82 x 10-15). Stippled bar: GWAS stage, striped bar: replication stage, and black bar: combined stage.
Biological and Pharmaceutical Bulletin Advance Publication