Associations between the neuron-specific glucocorticoid receptor (NR3C1) Bcl-1 polymorphisms and suicide in cancer patients within the first year of diagnosis
- Subin Park†1,
- Jin Pyo Hong†2,
- Jong-Keuk Lee3,
- Young-Mi Park3,
- Yangsoon Park4,
- Juri Jeon2,
- Myeong Hee Ahn5 and
- Se Chang Yoon2Email author
© The Author(s) 2016
Received: 16 February 2016
Accepted: 15 June 2016
Published: 11 July 2016
Cancer diagnosis is associated with an increased suicide risk, particularly within the first 1 year after diagnosis of cancer. Abnormal function of the hypothalamic–pituitary–adrenal axis has been implicated in the pathophysiology of depression and suicide. We examined genetic associations of the functional Bcl-1 polymorphism of (rs41423247) neuron-specific glucocorticoid receptor (NR3C1) gene, with death by suicide in cancer patients. Suicides occurring within a year of cancer diagnosis (‘early suicide’) were considered separately from those suicides during the second or subsequent year (‘late suicide’) after cancer diagnosis.
The subjects consisted of 343 cancer patients admitted to a general hospital in Seoul, South Korea from 1996 to 2009, of which 182 had died by suicide and 161 were alive on December 31, 2009. Genomic DNA was extracted from formalin-fixed paraffin-embedded tissue sample of patients with cancer. We conducted a case-control association analysis of Bcl-1 polymorphism of NR3C1 gene.
Subjects carrying the GG genotype of Bcl-1 polymorphism were at increased risk of early suicide when compared to those carrying the CC genotype (OR 3.80, 95 % CI 1.02–14.16, p = .047). Similarly, those individuals carrying the GG genotype (recessive mode) had an increased risk of early suicide relative to the CC or CG genotype (OR 3.71, 95 % CI 1.03–13.43, p = .045). However, there were no differences in the genotype distributions of the NR3C1 Bcl-1 polymorphism between late suicide cases and controls.
Our findings suggest that the NR3C1 Bcl-1 polymorphisms may be involved in the susceptibility to suicide within the first year after cancer diagnosis among cancer patients in Korean population.
KeywordsCancer Suicide NR3C1 Genetics
Cancer diagnosis is associated with an increased suicide risk [1–6]. Suicide risk is the highest within the first 1 year after diagnosis of cancer [1, 2, 7–9]. There was a significant decrease in the relative suicide risk over decades, and 5 years after diagnosis, the risk should be about the same as for the general population . This highest suicide risk in early period after cancer diagnosis may be associated with emotional distress induced by diagnosis of cancer, progressive and life-threatening illness [11, 12]. Based on diathesis-stress model of suicidal behavior , among patients who were exposed to major psychological stress (i.e. cancer diagnosis), some individuals may be more vulnerable than others to develop suicidal behavior, and there may be genetic and biological differences that may make someone more or less vulnerable to suicide. To our best knowledge, there has been no genetic association study of suicides in cancer patients.
Depression and stress are among the major risk factors for suicidal behavior. Abnormalities of the hypothalamic–pituitary–adrenal (HPA) axis have been consistently reported in patients with major depressive disorder (MDD)  and in individuals who completed suicide [15, 16]. Whereas about 50 % of patients with MDD have an abnormal dexamethasone suppression test (DST), almost all those MDD patients who eventually completed suicide had an abnormal DST . This suggests a significant association between an abnormal HPA axis and suicide, independently of MDD. The cortisol stress response has been identified as one of the most prominent candidate suicide endophenotypes [15–17].
The dysregulation of HPA axis function observed in depression and suicide is believed to be in part due to a disturbed feedback inhibition by endogenous corticoids [14, 17–19]. In humans, cortisol serves as a negative feedback regulator of HPA axis by its action on two types of receptors, glucocorticoid (GR) and mineralocorticoid (MR) receptors in the brain and the pituitary [20, 21]. It is believed that this impaired HPA axis negative feedback in depression and suicide may be related to altered GR and/or MR function. The most widely studied of polymorphisms associated with GR receptor expression is the Bcl-1 restriction fragment length polymorphism of the nuclear receptor subfamily 3, group C, member 1 (NR3C1) gene, located on chromosome 5q31. A biallelic polymorphism (Bcl-1, C > G) was identified downstream of the exon 2–intron 2 junction of the NR3C1 gene. The polymorphism might affect processing of the primary GR transcripts and cause increased GR sensitivity [22, 23], which affects strength of inhibitory feedback within the HPA axis . Bcl-1 high-function (G) allele carriers have been observed to display reduced cortisol responses following psychosocial stressors  and greater suppression of cortisol after dexamethasone administration . Such findings lead to the suggestion that psychosocial stressors, such as a cancer diagnosis, might have more pronounced effects in homozygous carriers of the G allele. Consistent with this, increased methylation that would cause reduced hippocampal NR3C1 expression has been observed in victims of suicide with a history of childhood abuse .
In the present study, we assessed genetic associations of the functional Bcl-1 polymorphism (rs41423247) of NR3C1 gene, with death by suicide in cancer patients, in particular suicide within the first year of diagnosis. To the best of our knowledge, this is the first association study of suicides in cancer patients, a homogenous stress-exposure group.
The sampling pool for this study consisted of 164,497 adult patients (age ≥20 years) who had been diagnosed with cancer at the University Hospital in Seoul, South Korea, between 1996 and 2009 according to the International Classification of Disease for Oncology tenth revision (ICD-10). This nationwide registry accounts for approximately 10 % of newly diagnosed cancer cases . Data were censored either on the date of death or on December 31, 2009. Person-years were calculated for each patient, and information on patient status was obtained from the database of the National Statistical Office (NSO), which compiles all death notices in South Korea.
The cause of death was established by linking to the NSO database that was matched with the hospital records based on the unique national identification number assigned to each Korean citizen.
Selection of cases and controls
Suicide death was defined by ICD-10 codes X60–X84 (intentional self-harm). Among the 164,497 patients with cancer, 373 had died by suicide by December 31, 2009. Among the 373 cancer patients who had died by suicide, 185 patients had tissue for genetic analysis. One control patient was randomly selected from the study cohort using the statistical analysis system (SAS) program (SAS Institute, Cary, NC, USA) and matched by sex, age (±3 years), and the time of diagnosis (±2 years) with each patient who had died by suicide. Information about demographic factors and anatomic site of cancer was obtained from electronic medical records by two medical doctors. Cancer sites were organized into seven diagnostic groups. These included stomach cancer (C16.x), colorectal cancer (C18.x–20.x), biliary-pancreatic cancer (C22.1, C23.9, and C24.x–25.x), breast cancer (C50.x), liver cancer (C22.0), lung cancer (C34.x), and cancers of other sites (i.e., diagnoses not corresponding to the previously mentioned codes). Suicides occurring within a year of cancer diagnosis (‘early suicide’) were considered separately from those suicides during the second or subsequent year (‘late suicide’) after cancer diagnosis.
DNA extraction and genotyping
Genomic DNA was extracted from formalin-fixed paraffin-embedded (FFPE) tissue sample of patients with cancer using a QIAamp® DNA FFPE Tissue Kit and the QIAGEN Deparaffinization solution (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Approximately 50 ng of genomic DNA was used for TaqMan® SNP genotyping. Among 370 samples, 23 samples below .1 µg were excluded from genetic analyses. Finally, the remaining 347 samples were included in the analysis. We genotyped NR3C1 Bcl-1 polymorphism with a TaqMan® SNP Genotyping Assays on a 7900HT Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). There were four missing genotypes and finally, the 343 samples were included in statistical analyses.
The distribution of the general demographic and clinical features between early suicide cases, late suicide cases, and controls was evaluated by using the 1-way analysis of variance test and the χ2 test for continuous and categorical variables, respectively. Agreement with Hardy–Weinberg equilibrium (HWE) for each SNP was tested using a goodness-of-fit χ2 test. Genotype distribution for each polymorphism was in agreement with the expected values of the HWE (p > .05). To investigate the association between the evaluated genotypes and death by suicide, binary logistic regression analyses were performed. For the main analyses, 3 set of dichotomized outcomes were defined as the all suicide cases versus controls, the early suicide cases versus controls, and the late suicide cases versus controls. A 2-tailed p value of <.05 was considered statistically significant. Data analyses were performed using the Statistical Package for Social Sciences (Version 21.0 for Windows; SPSS, Inc., Chicago, IL, USA). Power calculation was performed using the Quanto (http://www.hydra.usc.edu/gxe/) software to compute the statistical power of our study.
Characteristics of the study population
Controls (N = 161)
Early suicide (N = 47)
Late suicide (N = 135)
Age at diagnosis, mean (SD)
Days to death
Genotype distributions and allele frequencies of Bcl-1 polymorphism of NR3C1 gene between suicide cases and controls
Controls (N = 161)
Suicide (N = 182)
Early suicide (N = 47)
Late suicide (N = 135)
Suicide vs. controls
Early suicide vs. controls
Late suicide vs. controls
OR (95 % CI)
OR (95 % CI)
OR (95 % CI)
CG + GG
CC + CG
This is the first genetic association study on death by suicide of homogenous group exposed to severe psychological stressor, a cancer diagnosis. We found significant differences in genotypic and allelic frequencies of the NR3C1 Bcl-1 polymorphism between early suicide cases and non-suicide controls.
HPA-axis hyperactivity, assessed by means of the DST or 24 h urinary cortisol samples, has been associated with completed suicide among depressed patients. Several studies have reported that suicidal individuals or first-degree relatives of suicide decedents fail to mount a proper HPA axis response to stress, and the cortisol stress response has been identified as one of the most prominent candidate suicide endophenotypes [15–17]. NR3C1 Bcl-1 high function (G) allele carriers have been observed to display reduced cortisol responses following psychosocial stressors [28, 29] and greater suppression of cortisol after dexamethasone . Such findings lead to the conjecture that traumatic experiences might have more pronounced effects in individuals whose genetic propensities code for lesser modulatory feedback within the HPA axis.
In the present study, relative to cancer patients with the CC genotype, those with the GG genotype of Bcl-1 polymorphism had a 3.8-fold higher risk of early suicides, although there were no significant differences in genotype or allele frequencies distribution between total suicide cases and controls. The psychological stress induced by the cancer diagnosis itself may give rise to early suicides, whereas the disease burden, fatigue, exhaustion, and aggravation of preexisting psychiatric disorder may be associated with late suicide . Our results suggest that the Bcl-1 G allele, which reduces cortical response to psychosocial stressor (in this case, the cancer diagnosis), is related to early suicide in cancer patients.
Because of the limited availability of post-mortem brain tissues, there have been only a few studies examining protein and gene expression of GR in suicide completers [26, 30]. Pandey et al.  found that the mRNA and protein expression of GR is decreased in the prefrontal cortex and amygdale of teenage suicide victims compared with control subjects. In a postmortem brain study by McGowan et al. , increased methylation in the exon 1F NR3C1 promoter, that would cause reduced hippocampal NR3C1 expression, was observed in suicide completers with a history of childhood abuse. This result is consistent with our finding of contribution of NR3C1 gene in suicide completion of individuals exposed to severe psychological stressor.
A limitation of the study was the relatively small sample size. On calculating the statistical power of the study, we found that approximately 524 case samples are required to perform a matched case–control study with 80 % power to detect a relative risk of 2.0 in an estimated prevalence of suicide in cancer patients of 60/100,000 . However, in our study, we used a smaller sample size (161 controls and 182 cases). That may be one of the reasons why we could not find any significance in the association between total suicide cases and controls. Furthermore, based on our results shown in Table 2, we found that our study has 65.4 % power with a relative risk of 3.71. As our study has limited power, further study is necessary in an independent sample set with a large sample size. Additionally, we evaluated only one SNP on the NR3C1 gene. Although the NR3C1 gene may affect suicides in cancer patients, the influence of, and interactions between, many target genes, should be considered. Finally, as this was a retrospective chart-review study, we could not evaluate psychiatric symptoms after a cancer diagnosis, which are known to have an important influence on suicide in cancer patients.
Our findings suggest that the NR3C1 Bcl-1 polymorphisms may be involved in the susceptibility to suicide within the first year after cancer diagnosis among cancer patients in Korean population. The study must be replicated in studies using larger samples and with using more SNPs.
abnormal function of the hypothalamic–pituitary–adrenal
the functional Bcl-1 polymorphism
neuron-specific glucocorticoid receptor
major depressive disorder
SP contributed to the data analysis, and wrote the first drafts of the manuscript. JPH contributed to the design of the study, data collection process and revisions of the manuscript; JKL and YMP contributed to the DNA analysis; YP contributed to recruitment of samples; JJ was involved in the data analysis; SCY contributed to the design of the study and revisions of the manuscript. All authors read and approved the final manuscript.
This work was supported by the Korea Healthcare Technology R&D project, Ministry of Health and Welfare, Republic of Korea. [Grant Number A120051].
The authors declare that they have no competing interests. The study was approved by the institutional review board for human subjects at the Asan Medical Center, Ulsan University College of Medicine.
Availability of data and materials
Because this is a small dataset of single SNP, this data do not deserve to be deposited in the repositories. But dataset can be shared with any scientist who is interested.
Open AccessThis 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.
- Ahn E, Shin DW, Cho SI, Park S, Won YJ, Yun YH. Suicide rates and risk factors among Korean cancer patients, 1993–2005. Cancer Epidemiol Biomarkers Prev. 2010;19:2097–105.View ArticlePubMedGoogle Scholar
- Hem E, Loge JH, Haldorsen T, Ekeberg O. Suicide risk in cancer patients from 1960 to 1999. J Clin Oncol. 2004;22:4209–16.View ArticlePubMedGoogle Scholar
- Mann JJ. A current perspective of suicide and attempted suicide. Ann Intern Med. 2002;136:302–11.View ArticlePubMedGoogle Scholar
- Misono S, Weiss NS, Fann JR, Redman M, Yueh B. Incidence of suicide in persons with cancer. J Clin Oncol. 2008;26:4731–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Crocetti E, Arniani S, Acciai S, Barchielli A, Buiatti E. High suicide mortality soon after diagnosis among cancer patients in central Italy. Br J Cancer. 1998;77:1194–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Nasseri K, Mills PK, Mirshahidi HR, Moulton LH. Suicide in cancer patients in California, 1997–2006. Arch Suicide Res. 2012;16:324–33.View ArticlePubMedGoogle Scholar
- Fang F, Fall K, Mittleman MA, Sparen P, Ye W, Adami HO, et al. Suicide and cardiovascular death after a cancer diagnosis. N Engl J Med. 2012;366:1310–8.View ArticlePubMedGoogle Scholar
- Henriksson MM, Isometsa ET, Hietanen PS, Aro HM, Lonnqvist JK. Mental disorders in cancer suicides. J Affect Disord. 1995;36:11–20.View ArticlePubMedGoogle Scholar
- Robinson D, Renshaw C, Okello C, Moller H, Davies EA. Suicide in cancer patients in South East England from 1996 to 2005: a population-based study. Br J Cancer. 2009;101:198–201.View ArticlePubMedPubMed CentralGoogle Scholar
- Allebeck P, Bolund C, Ringback G. Increased suicide rate in cancer patients. A cohort study based on the Swedish Cancer-Environment Register. J Clin Epidemiol. 1989;42:611–6.View ArticlePubMedGoogle Scholar
- Iwamitsu Y, Shimoda K, Abe H, Tani T, Okawa M, Buck R. Anxiety, emotional suppression, and psychological distress before and after breast cancer diagnosis. Psychosomatics. 2005;46:19–24.View ArticlePubMedGoogle Scholar
- Ahn MH, Park S, Lee HB, Ramsey CM, Na R, Kim SO, et al. Suicide in cancer patients within the first year of diagnosis. Psychooncology. 2015;24:601–7.View ArticlePubMedGoogle Scholar
- Mann JJ, Arango V. Integration of neurobiology and psychopathology in a unified model of suicidal behavior. J Clin Psychopharmacol. 1992;12:2S–7S.View ArticlePubMedGoogle Scholar
- Holsboer F. The corticosteroid receptor hypothesis of depression. Neuropsychopharmacology. 2000;23:477–501.View ArticlePubMedGoogle Scholar
- Mann JJ, Arango VA, Avenevoli S, Brent DA, Champagne FA, Clayton P, et al. Candidate endophenotypes for genetic studies of suicidal behavior. Biol Psychiatry. 2009;65:556–63.View ArticlePubMedPubMed CentralGoogle Scholar
- McGirr A, Diaconu G, Berlim MT, Pruessner JC, Sable R, Cabot S, et al. Dysregulation of the sympathetic nervous system, hypothalamic-pituitary-adrenal axis and executive function in individuals at risk for suicide. J Psychiatry Neurosci. 2010;35:399–408.View ArticlePubMedPubMed CentralGoogle Scholar
- Coryell W, Schlesser M. The dexamethasone suppression test and suicide prediction. Am J Psychiatry. 2001;158:748–53.View ArticlePubMedGoogle Scholar
- Nemeroff CB. The corticotropin-releasing factor (CRF) hypothesis of depression: new findings and new directions. Mol Psychiatry. 1996;1:336–42.PubMedGoogle Scholar
- Pariante CM, Miller AH. Glucocorticoid receptors in major depression: relevance to pathophysiology and treatment. Biol Psychiatry. 2001;49:391–404.View ArticlePubMedGoogle Scholar
- Bao AM, Swaab DF. Corticotropin-releasing hormone and arginine vasopressin in depression focus on the human postmortem hypothalamus. Vitam Horm. 2010;82:339–65.View ArticlePubMedGoogle Scholar
- Jacobson L. Hypothalamic-pituitary-adrenocortical axis regulation. Endocrinol Metab Clin North Am. 2005;34:271–92 (vii).View ArticlePubMedGoogle Scholar
- Fleury I, Beaulieu P, Primeau M, Labuda D, Sinnett D, Krajinovic M. Characterization of the Bcl-1 polymorphism in the glucocorticoid receptor gene. Clin Chem. 2003;49:1528–31.View ArticlePubMedGoogle Scholar
- van Rossum EF, Koper JW, van den Beld AW, Uitterlinden AG, Arp P, Ester W, et al. Identification of the Bcl-1 polymorphism in the glucocorticoid receptor gene: association with sensitivity to glucocorticoids in vivo and body mass index. Clin Endocrinol (Oxf). 2003;59:585–92.View ArticleGoogle Scholar
- DeRijk RH, Schaaf M, de Kloet ER. Glucocorticoid receptor variants: clinical implications. J Steroid Biochem Mol Biol. 2002;81:103–22.View ArticlePubMedGoogle Scholar
- Wust S, Federenko IS, van Rossum EF, Koper JW, Kumsta R, Entringer S, et al. A psychobiological perspective on genetic determinants of hypothalamus-pituitary-adrenal axis activity. Ann NY Acad Sci. 2004;1032:52–62.View ArticlePubMedGoogle Scholar
- McGowan PO, Sasaki A, D’Alessio AC, Dymov S, Labonte B, Szyf M, et al. Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse. Nat Neurosci. 2009;12:342–8.View ArticlePubMedPubMed CentralGoogle Scholar
- The Korea Central Cancer Registry [KCCR]. Annual report of cancer statistics in Korea in 2009. Seoul: The Korea Central Cancer Registry; 2011.Google Scholar
- Wust S, Van Rossum EF, Federenko IS, Koper JW, Kumsta R, Hellhammer DH. Common polymorphisms in the glucocorticoid receptor gene are associated with adrenocortical responses to psychosocial stress. J Clin Endocrinol Metab. 2004;89:565–73.View ArticlePubMedGoogle Scholar
- Kumsta R, Entringer S, Koper JW, van Rossum EF, Hellhammer DH, Wust S. Sex specific associations between common glucocorticoid receptor gene variants and hypothalamus-pituitary-adrenal axis responses to psychosocial stress. Biol Psychiatry. 2007;62:863–9.View ArticlePubMedGoogle Scholar
- Pandey GN, Rizavi HS, Tripathi M, Ren X. Region-specific dysregulation of glycogen synthase kinase-3beta and beta-catenin in the postmortem brains of subjects with bipolar disorder and schizophrenia. Bipolar Disord. 2015;17:160–71.View ArticlePubMedGoogle Scholar