http://orcid.org/0000-0003-2095-8681
Figures
Abstract
Breast cancer has a multifactorial etiology. One of the supposed and novel mechanisms is an alteration of circadian gene expression. Circadian genes play a crucial role in many physiological processes. These processes, such as genomic stability, DNA repair mechanism and apoptosis, are frequently disrupted in breast tumors. To assess the significance of circadian gene expression in breast cancer, we carried out an analysis of CLOCK, BMAL1, NPAS2, PER1, PER2, PER3 and CRY1, CRY2, TIMELESS, CSNK1E expression by the use of the quantitative Real-Time PCR technique in tumor tissue and non-tumor adjacent normal tissue sampled from 107 women with a newly diagnosed disease. The obtained data were compared to the clinical and histopathological features. PER1, PER2, PER3, CRY2 were found to be significantly down-expressed, while CLOCK, TIMELESS were over-expressed in the studied tumor samples compared to the non-tumor samples. Only gene expression of CRY1 was significantly down-regulated with progression according to the TNM classification. We found significantly decreased expression of CRY2, PER1, PER2 genes in the ER/PR negative breast tumors compared to the ER/PR positive tumors. Additionally, expression of CRY2, NPAS2 genes had a decreased level in the poorly differentiated tumors in comparison with the well and moderately differentiated ones. Our results indicate that circadian gene expression is altered in breast cancer tissue, which confirms previous observations from various animal and in vitro studies.
Citation: Lesicka M, Jabłońska E, Wieczorek E, Seroczyńska B, Siekierzycka A, Skokowski J, et al. (2018) Altered circadian genes expression in breast cancer tissue according to the clinical characteristics. PLoS ONE 13(6): e0199622. https://doi.org/10.1371/journal.pone.0199622
Editor: William B. Coleman, University of North Carolina at Chapel Hill School of Medicine, UNITED STATES
Received: November 22, 2017; Accepted: May 2, 2018; Published: June 29, 2018
Copyright: © 2018 Lesicka et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are only within the paper.
Funding: This work was supported by the National Science Center (Grant No. 2014/15/N/NZ5/01671) (to ML), https://www.ncn.gov.pl/?language=en, and partially supported by Ministry of Science and Higher Education Poland DIR/WK/2017/01 (to LK), http://www.nauka.gov.pl/en/, and Internal grant for statutory activity at Nofer Institute of Occupational Medicine IMP 14.4/2018. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Circadian rhythm is the most important cyclic changes system in all living organisms across the world, generated by the endogenous circadian clock. Circadian rhythm plays a crucial role in the maintenance of genomic stability and homeostasis. Globalization and fast-paced life of the contemporary world contribute to widespread changes in circadian coordination, resulting in civilizational diseases including metabolic syndromes, cardiovascular diseases, ageing, psychological disorders, neurodegenerative diseases as well as many types of cancer [1].
The core pacemaker of circadian rhythm is located in the suprachiasmatic nucleus of the hypothalamus, which synchronizes with both: endogenous and external factors [2]. Peripheral control systems are located in many tissues including breast tissue [3]. The molecular circadian clock at a transcriptional level consists of two groups of genes: positive and negative regulators. Transcriptional factors, i.e., CLOCK, BMAL1, and NPAS2 that form heterodimers CLOCK/BMAL1 and NPAS2/BMAL1 belong to the first group. These heterodimers initiate transcription via E-boxes in promoters of the second group, i.e.: CRY1, CRY2, PER1, PER2, PER3. CRYs and PERs proteins inhibit heterodimers complexes as a negative loop of the circadian clock. The auxiliary loop consisting of REV-ERB/ROR ensures stabilization of this system [4, 5]. The described above physiological processes occur consecutively in each normal cell; however, in cancerous cells, this order is not preserved [6, 7]. Cancer disease with a complex multi-step etiology is manifested by abnormal physiology on a cellular level. Currently, there is more and more evidence, that circadian clock genes expression may play a crucial role in mammary gland pathophysiology. Circadian rhythm disruption may initiate various tumor-related processes [8]. A relevant process for cells, i.e. cell cycle is indirectly regulated by clock genes. Therefore, alterations of this process are manifested by abnormal cell division in many types of cancer [9]. This is probably caused by the disrupted synthesis of melatonin—nocturnal hormone, due to a permanent exposure to artificial light during the night, which particularly relates to night shift workers [10]. Clock genes belong to tumor suppressor genes that can control cell cycle, apoptosis, DNA damage and repair system [11]. The basic genetic regulation of the circadian clock is highly conserved across living organisms. The first evidence that deregulation of circadian rhythm due to environmental factors initiates breast carcinogenesis was hypothesized by Hamilton in 1969 [12]. Disruptions of the circadian gene expression pattern as an early event that leads to breast cancer initiation have been hypothesized by Stevens [13].
Disruption of circadian machinery in the human body may be influenced by three factors such as genetic defects, ageing of the body and lifestyle, also including shift work system [6, 14]. These factors together or individually may lead to dysregulation of the circadian clock, which is manifested by an altered cell cycle, reduced apoptosis, altered cell metabolism and disturbed melatonin synthesis. All the mentioned factors reveal a potential link with breast carcinogenesis [15]. An altered pattern of circadian gene expression is observed in the case of in vitro models, rodent animals as well as in human genetic association studies [16]. Clock genes alterations in mammary epithelium may lead to cell cycle dysregulation, which is caused by aberrant cell divisions, which can be associated with tumor progression and a more aggressive type of breast cancer tumors [17].
The main objective of the study was to present expression pattern of ten core clock genes (CLOCK, NPAS2, BMAL, CRY1, CRY2, PER1, PER2, PER3, TIMELESS, and CSNK1E) in breast cancer tissues in comparison to the paired non-tumor adjacent tissue according to the histological and clinical features, including. estrogen (ER), progesterone (PR) receptor status, HER2 receptor status, molecular type of the breast cancer, Tumor Node Metastasis stage (TNM) classification and grading.
Materials and methods
Ethics declaration
This study was conducted in compliance with the Declaration of Helsinki and was approved by relevant Local Ethics Committees (Independent Ethics Committee at Medical University of Gdansk, resolution No. NKEBN/781/2005, and Ethical Institutional Review Board at the Nofer Institute of Occupational Medicine resolution No. 20/2014). We have obtained the written consent of participants to deposit their tissue samples in the Bank of Frozen Tissues and Genetic Specimens of the Medical University of Gdansk and forward for further analyzes.
Study population
We analyzed 107 Polish women, all Caucasian origin with a newly diagnosed breast cancer and without previous chemotherapy. The patients were admitted to the Clinic of Oncologic Surgery at the Medical University of Gdansk over the years 2006–2015 and divided into four groups according to the TNM classification. During the tumor removal surgery, tumor tissue and tumor-adjacent normal tissue (0.6 cm x 0.6 cm) were collected from each patient, then the pairs of tissue samples were frozen in -70°C and deposited in Bank of Frozen Tissues and Genetic Specimens of Medical University of Gdansk. Finally, they were shipped to the Nofer Institute of Occupational Medicine in Lodz (NIOM) for a further molecular analysis. The patient’s characteristics are presented in Table 1. All the tissue specimens were evaluated by pathologists to confirm the diagnosis of both tumor and adjacent normal tissue that was collected close to the resection margins. Subsequently, samples were specified according to estrogen receptor status (ER), progesterone receptor status (PR), Human Epidermal Growth Factor Receptor 2 (HER2) receptor status. Antigen retrieval and staining were performed by immunohistochemistry using the automatic devices: BenchMark GX (Ventana Medical Systems, Roche; Tucson, Arizona, AZ, USA) for all mentioned receptors and additionally fluorescent in situ hybridization analysis (only for equivocal HER2 positive samples with IHC 2+ score). Following criteria were included: estrogen receptor expression was determined as positive—IHC ranged from 4 to 8 as the intensity of the expression of ER receptor and as negative—IHC ranged from 0 to 3 as the intensity of the expression of ER receptor. Moreover, patients underwent a comprehensive histo-clinicopathological and post-surgical examination of tumor size (T-status), lymph node (N) status, and tumor grade (G1- well differentiated; G2- moderately differentiated; G3- poorly differentiated) histological type of breast cancer (Carcinoma ductale, Carcinoma lobulare). Additionally, molecular type of breast cancer were determined, according to this characteristics patients were divided into four molecular subtypes of breast cancer including Luminal A (estrogen-receptor and/or progesterone-receptor positive, HER2 negative) Luminal B (estrogen-receptor and/or progesterone-receptor positive HER2 positive) Basal-like (estrogen-receptor and progesterone-receptor negative, HER2 negative), HER2-enriched (estrogen-receptor and progesterone-receptor negative, HER2 positive).
RNA isolation procedures
Total RNA was isolated by the use of the RNeasy Lipid Tissue KIT (Qiagen, Hilden, Germany,) according to the manufacturer’s protocol. The amount of the total extracted RNA was determined using the spectrophotometer method—microplate (Multiskan™ GO Microplate Spectrophotometer, Thermo Fisher; Waltham, Massachusetts, MA, USA).
Gene expression analysis
Expression analysis of CLOCK, NPAS2, BMAL, CRY1, CRY2, PER1, PER2, PER3, TIMELESS and CSNK1E in the breast tumor tissue and the adjacent tumor tissue was performed by means of quantitative Real-Time PCR (qRT-PCR). Primers were designed with BEACON DESIGNER 7.01 (PREMIER Biosoft International, Palo Alto, CA, USA). Expression of core circadian genes was quantified with the FastStart Essential DNA Green Master (Roche, Indianapolis, IN, USA) and using glyceraldehydes 3-phosphate dehydrogenase (GAPDH) and 60S acidic ribosomal protein P0 (RPLP0) as endogenous controls which present stable expression levels in RNA isolated from the breast tumor tissue and adjacent normal breast tissue samples. cDNA of the tumor and tumor-adjacent samples were synthesized on 200ng RNA with Transcriptor First Strand cDNA Synth, Kit 200 (Roche, Indianapolis, IN, USA). All the gene expression procedures fulfilled the MIQE guidelines (Minimum Information for Publication of Quantitative Real-Time PCR Experiment). PCR efficiency was calculated using dilutions of five pooled cDNA samples, which were reversely transcribed from the total human RNA isolated from the breast tissue samples. Each of the predesigned primer sets should meet the criterion of 90% qRT-PCR efficiency. Randomly selected samples were simultaneously amplified on the same plate to determine the inter-assay coefficients of variability (CV). The samples were assayed on the PTC-200 thermocycler (MJ Research, Bio-Rad Laboratories, Inc. Hercules, California, CA, USA) (cDNA synthesis) and LightCycler® 96 (Roche, Indianapolis, IN, USA) (qRT-PCR assay). In order to calculate gene expression, we used reference gene-normalized relative quantification with efficiency correction for the target and reference genes (based on Livak method). Values of the Mean Normalized Expression (MNE) were log10 transformed and multiplied by 1000 to obtain a normal distribution of variables.
Statistical analysis
Continuous variables were expressed as means with standard error or medians with interquartile ranges. Gene expression levels in the breast cancer tissue samples were compared with the paired adjacent non-tumorous tissue samples, calculated using the mean normalized expression and log-transformed so as to achieve normal distribution. Normality was assessed using the Shapiro–Wilk W-test. Comparison between the multiple groups was performed using the analysis of variance (ANOVA), followed by the in between-group comparisons with the Tukey HSD test if ANOVA yielded a significant difference. Differences among the paired tissue sets were evaluated by the paired t-test or the Wilcoxon signed rank sum test. Additionally, p-value has adjusted for age, BMI, smoking status, menopausal status by linear regression. Significance was declared for p levels less than 0.05. Statistica12.0 (StatSoft, Tulsa, OK, USA) was used for the statistical analysis.
Results
Subjects characteristics
Table 1 shows clinicopathological features of the patients with primary breast cancer without previous chemotherapy. A total of 107 patients (60.22±13.17 years of age and BMI 26.84 ±4.57) were enrolled into the study. Seventy-four and seventy-seven percent of the patients were non-smokers (n = 80) and 24.20% were smokers (n = 26). One patient had an unknown smoking status (0.93%). Majority of the breast cancer patients were classified as ER, PR positive (n = 76, 71.03%; n = 68, 63.55%, respectively) in comparison to the patients with negative ER, PR status (n = 27, 25.23%; n = 35, 32.71, respectively). This is consistent with the literature, which states that about 70% of breast cancer patients are ER, PR positive. Four women did not specify their hormonal receptor status. Fifty-nine patients were HER2 negative (55.14%), while the HER2 positivity rate was 37.38% (n = 40). TNM was specified for all the patients. We divided them into four groups: T1N0 (n = 30; 28.04%); T1N1-3(n = 24; 22.43%); T3N0 (n = 23; 21.49%); T3N1-3 (n = 30; 28.04%).
Disrupted circadian gene expression in the breast cancer tissue samples compared to the non-tumor tissue samples
In order to determine whether circadian genes expression, including CLOCK, NPAS2, BMAL, CRY1, CRY2, PER1, PER2, PER3, TIMELESS, and CSNK1E, was altered in the breast tumor samples, we performed quantitative Real-Time PCR with intercalating SYBR Green in the 107 pairs of the tumor tissue and non-tumor adjacent normal breast tissue samples. We did not obtain detectable transcript levels for all the tissue samples in the total analysis, thus, the number of tissue specimens for the statistical analysis was limited and it is included in presented tables. (Tables 2 and 3). Our results indicated that the analyzed genes presented different expression patterns compared to the adjacent normal breast tissue samples. The strong statistical significance of the correlation between the investigated genes was observed in the case of the tumor tissue samples in comparison with the correlation analysis between the genes in the adjacent normal tissue samples and between types of the tissues. The transcript level was significantly up-regulated in the breast cancer tissue samples in comparison with the adjacent normal tissue samples for CLOCK (β-coefficient 0.219, p = 0.002) and TIMELESS (β-coefficient 0.298, p<0.0001). A significantly decreased mRNA level in the breast cancer tissue samples was observed for CRY2 (β-coefficient -0.424, p<0.0001), PER1 (β-coefficient -0.371, p<0.0001), PER2 (β-coefficient -0.149, p = 0.037) and PER3 (β-coefficient -0.231, p = 0.001). However, we did not observe a statistically significant correlation in the gene expression of BMAL1, CRY2, NPAS2, CSNK1E in the tissue pairs (Table 2).
Circadian gene expression according to the TNM classification
To investigate whether circadian genes expression is associated with the tumor stage, we analyzed the breast cancer tumor tissue samples according to the TNM classification (Table 3). The analysis included four groups consisting of T1N0–29 samples, T1N1-3–24 samples, T3N0–22 samples and T3N-3–30 samples. A significant difference (p<0.001) was reported only for CRY1. Additionally, the analysis showed marginal statistical significance for NPAS2. In general, some changes in the transcript level between the groups were observed; however, they were not statistically significant.
Decreased circadian gene expression in the negative ER and PR breast cancer tissue samples
Hormone receptor status is an essential clinical factor in the evaluation of breast cancer diagnosis which determines a proper therapy. We carried out an analysis of each circadian gene expression according to the ER/PR status in almost each breast cancer tissue sample. Associations for few genes were found. The ER-negative (ER-) breast cancer tissue samples had significantly lower PER1 gene expression than the ER-positive breast cancer tissue samples (ER+) (β-coefficient -0.291 p = 0.005). A similar marginal significant difference was observed for CRY2 (β-coefficient -0.178 p = 0.088). For other genes, we did not observe any changes in the transcript levels, between ER- and ER+ (Fig 1). A further comparison by PR receptor status showed that the patients with PR-negative tumors (PR-) had significantly lower levels of CRY2, PER,1 PER2 (β-coefficient -0.266; p = 0.008; β-coefficient -0.335, p = 0.001; β-coefficient -0.217; p = 0.034; respectively) and marginally significant lower levels of BMAL1, PER3, NPAS2 gene expression (β-coefficient -0.177 p = 0.09, β-coefficient -0.187; p = 0.068, β-coefficient -0.215 p = 0.057, respectively) than the patients with PR positive tumors (PR+) (Fig 1). This may suggest that tumor suppressing circadian genes (PER1,2,3, NPAS2) are particularly relevant for ER/PR negative tumorigenesis. All the specimens were classified into the molecular groups according to the presence of receptors on their surface: Luminal A, Luminal B, basal-like, HER2 overexpressed. mRNA level of each gene was compared to the receptor status including ER, PR, and HER2 (ER/PR positive, HER2 negative). The data analysis showed a significant difference between the molecular types for CRY2, PER1, PER2, PER3 (p<0.01, p = 0.04, p = 0.03, p = 0.03, respectively) (Fig 1). For other genes, no significant results were obtained (data not shown).
(A) Estrogen receptor status. (B) Progesterone receptor status. (C) Molecular type of breast cancer. Results presented as a mean of the Mean Normalized Expression (MNE)±SEM; statistical significance marked as *—p <0.05; **—p <0.01 obtained by Student’s t-test or One way ANOVA test.
Decreased circadian gene expression in the positive HER-2 breast cancer tissue samples
Expression status of proto-oncogene human epidermal growth factor receptor -2 (HER-2) in cancer cells determines the aggressiveness of breast cancer among the patients. It is used as a prognostic and predictive marker in breast cancer treatment. Our data was correlated with the HER-2 receptor status. Significant associations were observed for CRY2, PER3, PER2 (β-coefficient -0.266 p = 0.009, β-coefficient -0.289 p = 0.005, β-coefficient -0.217 p = 0.039, respectively). mRNA level of those genes was decreased in the HER2 positive tumors (Fig 2).
Results presented as mean of the Mean Normalized Expression (MNE)±SEM statistical significance marked as *—p <0.05; **—p <0.01 obtained by Student’s t-test.
Expression of circadian genes according to the histological type of carcinoma
All the analyzed samples were characterized by histopathological features. We divided the tissues into 3 groups–Carcinoma ductale, Carcinoma lobulare and mix type of carcinomas including Adenocarcinoma mucinosum, Carcinoma adenoides cysticum, Carcinoma ductale partim lobulare, Carcinoma metaplasticum mamame ac diferentiatio planoepitheliele, Carcinoma mucinosum, Carcinoma tubulolobulare, Carcinosarcoma. Significantly decreased transcript levels of CRY1, PER2 were observed in Carcinoma ductale and mix types of carcinomas (Fig 3).
Mix type of carcinomas include: Adenocarcinoma mucinosum, Carcinoma adenoides cysticum, Carcinoma ductale partim lobulare, Carcinoma metaplasticum mamame ac diferentiatio planoepitheliele, Carcinoma mucinosum, Carcinoma tubulolobulare, Carcinosarcoma. Results presented as a mean of the Mean Normalized Expression (MNE)±SEM statistical significance marked as *—p <0.05; **—p <0.01 obtained by One-way ANOVA test.
Expression of circadian genes according to the tumor grading
According to the tumor grading, breast cancer tissue samples were divided into three groups: G1-Well differentiated, G2- moderately differentiated, G3-poorly differentiated. Due to the limited number of tissue samples in the G1 group, this group was combined with the G2 for the statistical analysis. Data analysis showed that genes CRY2, PER1, NPAS2 had a decreased level of mRNA expression in the G3-poorly differentiated tumors in comparison with the G1- well and the G2-moderately differentiated tumors (β coefficient 0.28, p = 0.004; β coefficient 0.39 p = 0.00005; β coefficient -0,3 p = 0.009 respectively) (Fig 4).
Results presented as a mean of the Mean Normalized Expression (MNE)±SEM; statistical significance marked as *—p <0.05; **—p <0.01 obtained by Student’s t-test.
Discussion
In general, circadian genes are deregulated in breast cancer patients. Therefore, they may play an important role in cancer-related processes [16, 18–20]. So far studies suggest that the disruption of the circadian pattern may lead to vast consequences at a molecular level—previous evidence has shown that 2–10% of all mammalian genes display circadian oscillation [21] [22]. Our results support the hypothesis that disruption of clock genes in breast cancer may play an important role in the disease-related pathways. To elucidate the association between breast cancer and circadian rhythm, along with better understanding of breast tumorigenesis, we decided to correlate expression of the core clock genes with the clinical and histological characteristics. Circadian gene expression was detected in the paired tumor and non-tumor adjacent tissue samples obtained from the newly diagnosed breast cancer patients. Expression of some investigated circadian genes was significantly decreased in the cancerous tissue samples in comparison with the normal tissue samples from the breast cancer patients. The present data indicate that expression of TIMELESS, CLOCK in breast cancer tissue was over-expressed unlike the reduced transcript level of CRY2 and PER1, 2, 3. So far, nine studies on circadian genes expression in human breast cancer tissues have been performed [18, 20, 23–29]. The first study on circadian gene expression in human breast cancer has been conducted by Chen et. al. Using the immunohistochemical assay they have found that PER proteins were deregulated in breast cancer tissue [18]. Similar observations have been pointed out by Kuo et. al. [20], who have found that homogeneous expression of PER2 was associated with lymph node metastasis and a poor prognosis [20]. Some studies have indicated similar results showing that genes PER1, PER2, PER3, and CRY2 were down-expressed in breast cancer tissues. On the other hand, CLOCK and TIMELESS were found to be over-expressed in tumor tissue in comparison with normal adjacent breast cancer tissues [16].
Our results comply with the previous in vitro and animal studies related to circadian disturbances in breast cancer. Firstly, the described alterations in the functioning of circadian rhythm have been observed in animal and in vitro models. Downregulation of Per1 and Per2 in rodent models initiates higher cancer cell growth, but on the other hand, overexpressed Per2 leads to apoptosis in mouse mammary carcinoma cell lines ETM6. Per2 has a leading function in the circadian rhythm regulation [30]. Per2 directly controls c-myc—a gene responsible for proper functioning of the cell cycle [30]. Per2 mutant mice have a deregulated c-myc and p53 gene functions, leading to apoptosis reduction and cancer development after γ radiation exposure [30]. Furthermore, PER and CRY genes have been shown to be engaged in cell proliferation, including roles in DNA damage checkpoint control and regulation of important genes of cell cycle progression protein like c-myc, Wee1, Myt1, Cyclin B, Cdc2 [31]. Due to the features specified above PER genes can act as tumor suppressor genes [30]. A few studies have shown a low level of PER and CRY expression in human breast cancer [18, 25, 26, 31–33]. Similar functions have been observed in human studies, where PER1 and PER2 in normal cells promote apoptosis and act as tumor suppressor genes. In the case of breast cancer, a decreased transcript level of these genes is observed. PER1, PER2 indirectly suppress c-myc transcription by inhibiting E-box mediated transactivation by BMAL1/NPAS2. The important role of circadian genes as one essential part of cell cycle machinery is stimulation of ATM and CHK2 gene expression. Deficiency of PER1 influences the altered expression of cell cycle checkpoint genes. It can be one of the early events of cancer initiation because of deregulations of these genes in breast cancer. In vitro studies on breast cancer cell lines show that the induced overexpression of PERs gene elicits cell cycle arrest, growth inhibition of cancerous cells and induces apoptosis. CRYs play a similar role in tumorigenesis activation via the cell cycle. Deficiency of CRY leads to the deregulated expression of Wee-1, Cyclin D1 and reveals disrupted cell cycle regulation. Our experiment showed that expression of CRYs, PERs are statistically decreased in the breast cancer tissue samples compared to the normal tissues samples.
What is more, we compared expression levels of clock genes to the clinical and pathological features such as status of estrogen, progesterone and HER2 receptors, type of histological type of breast cancer, stage of tumor progression and tumor grading. A significantly increased CRY1 expression level was observed in the most advanced stage—T3N1-3 in comparison with T1N0. This fact might be related to the defence mechanisms. According to the molecular type of breast cancer, the statistically significant lowest expression of PER1 was observed in the triple negative breast cancer patients, while in the case of HER2-enriched, the lowest expression was observed for the negative loop of the central clock genes such as CRY2, PER2, PER3. We found a strong statistical association between the core circadian genes and hormone receptor status as well as tumor grading. They displayed significantly lower expression in the ER-negative, PR-negative and poor-differentiated tumors. These results are in agreement with the previously published studies showing how the loss of circadian gene expression may lead to the loss of cell proliferation control. Breast cancer patients with a negative hormone receptor status for ER and PR are less responsive to the treatment and have a worse course of the disease and progression. Thus, a suggestion that clock genes expression (like CRYs and PERs) and their higher activity in ER/PR positive tumors in comparison with ER/PR negative may imply an important relationship between ER, PR—signalling and circadian rhythm. This hypothesis is supported by recent studies conducted by Gery et al., who have found a link between ER status and PER2 expression [33]. Cancer cell line MCF-7 with an ER-positive status treated with E2 had higher expression of PER2 than cell lines without ER [33]. In the PER2 promoter, they have found potential Estrogen Responsive Element (ERE) binding site at 356 bp. The same authors have found a significant relationship between the low level of PER1 expression and the negative ER status [33]. Similar observations for CRY genes have been made by Hoffman et al. explaining why circadian genes expression is decreased in ER-negative tumors.
In addition, it is generally recognized that estrogen stimulation is one of the major initiators of breast tumorigenesis. When the ER status is taken into account, patients with ER-positive breast cancer possess a lower level of CLOCK expression than those with ER-negative tumors. This data imply that abnormal overexpression of CLOCK might be one of the early events in breast tumorigenesis [23]. Further experiments conducted by Xiao et al. have revealed opposite results. As much as 74% of examined ER-positive breast cancer tissues had a higher CLOCK expression than ER-negative samples measured by immunohistochemical assay [34]. A very interesting observation has been recently made by the same group. They have found that CLOCK is required for the proliferation of breast cancer. Those authors suggest that CLOCK, together with ER, could promote cancer cell proliferation. Breast cancer MCF-7 cell line with overexpression of CLOCK and ER showed higher growth than cell lines transfected with empty vector. This discovery indicates that it seems reasonable to assume that induction of CLOCK expression via ER is a core event in the proliferation of ER-positive breast cancer cells [34]. In our experiment, we did not observe similar directions. The transcript level did not differ according to the various tumor features. We only observed higher CLOCK expression in the tumor tissue samples in comparison with the non-tumor adjacent tissue samples.
There are not many studies on TIMELESS expression and its association with breast cancer. Our results are similar to those obtained by Mao et al. that indicate TIMELESS overexpression is different in breast tumor specimens compared to normal tissues. They have also identified a cancer-relevant network of transcripts with altered expression following TIMELESS knockdown, which contained many genes with known functions in cancer development and progression. Additionally, cell proliferation of TIMELESS knockdown MCF7 cell line was significantly decreased [35]. A recent study confirmed that TIMELESS levels in both breast cancer cell lines and tissues were significantly over-expressed. Upregulation of TIMELESS dramatically enhanced, while knockdown of TIMELESS suppressed the self-renewal of cancer stem cells (CSCs), cell invasion and migration abilities of breast cancer cells in vitro [36].
A general observation from this study is a statistically significant diminished expression level of CRY2, and PERs and an increased expression level of CLOCK and TIMELESS in the breast cancer tissue samples in comparison with the adjacent normal tissue samples. Another observation from our research indicates downregulation of CRYs and PERs in the ER/PR negative breast cancer tissue samples and in Carcinoma ductale histological type.
Circadian genes are involved in many physiological functions and their alterations could be manifested in various neoplasms, including breast cancer. According to the existing data regarding CLOCK-related modulation of breast cancer development, we shall focus on transcriptional and epigenetic regulation of the core clock genes assuming that: 1) circadian genes can modulate cancer suppression via caretaker, gatekeeper and landscaper genes and direct modification of chromatin structure; 2) expression of the clock genes seems to be epigenetically modulated, 3) there is cross-talk between the clock genes and ER status.
References
- 1. Kettner NM, Katchy CA, Fu L. Circadian gene variants in cancer. Ann Med. 2014;46(4):208–20. Epub 2014/06/05. pmid:24901356.
- 2. Weaver DR. The suprachiasmatic nucleus: a 25-year retrospective. J Biol Rhythms. 1998;13(2):100–12. pmid:9554572.
- 3. Dibner C, Schibler U, Albrecht U. The mammalian circadian timing system: organization and coordination of central and peripheral clocks. Annu Rev Physiol. 2010;72:517–49. Epub 2010/02/13. pmid:20148687.
- 4. Reppert SM, Weaver DR. Molecular analysis of mammalian circadian rhythms. Annu Rev Physiol. 2001;63:647–76. pmid:11181971.
- 5. Ukai H, Ueda HR. Systems biology of mammalian circadian clocks. Annu Rev Physiol. 2010;72:579–603. pmid:20148689.
- 6. Evans JA, Davidson AJ. Health consequences of circadian disruption in humans and animal models. Prog Mol Biol Transl Sci. 2013;119:283–323. pmid:23899601.
- 7. Lee S, Donehower LA, Herron AJ, Moore DD, Fu L. Disrupting circadian homeostasis of sympathetic signaling promotes tumor development in mice. PLoS One. 2010;5(6):e10995. Epub 2010/06/07. pmid:20539819.
- 8. Kelleher FC, Rao A, Maguire A. Circadian molecular clocks and cancer. Cancer Lett. 2014;342(1):9–18. Epub 2013/10/04. pmid:24099911.
- 9. Filipski E, King VM, Li X, Granda TG, Mormont MC, Liu X, et al. Host circadian clock as a control point in tumor progression. Journal of the National Cancer Institute. 2002;94(9):690–7. pmid:11983758.
- 10. Stevens RG. Artificial lighting in the industrialized world: circadian disruption and breast cancer. Cancer causes & control: CCC. 2006;17(4):501–7. pmid:16596303.
- 11. Fu L, Kettner NM. The circadian clock in cancer development and therapy. Prog Mol Biol Transl Sci. 2013;119:221–82. pmid:23899600.
- 12. Hamilton T. Influence of environmental light and melatonin upon mammary tumour induction. Br J Surg. 1969;56(10):764–6. pmid:5347580.
- 13. Stevens RG, Davis S. The melatonin hypothesis: electric power and breast cancer. Environ Health Perspect. 1996;104 Suppl 1:135–40. Epub 1996/03/01. pmid:8722117.
- 14. Stevens RG, Brainard GC, Blask DE, Lockley SW, Motta ME. Breast cancer and circadian disruption from electric lighting in the modern world. CA Cancer J Clin. 2014;64(3):207–18. pmid:24604162.
- 15. Stevens RG. Circadian disruption and breast cancer: from melatonin to clock genes. Epidemiology. 2005;16(2):254–8. pmid:15703542.
- 16. Reszka E, Przybek M. Circadian Genes in Breast Cancer. Adv Clin Chem. 2016;75:53–70. Epub 2016/04/15. pmid:27346616.
- 17. Blakeman V, Williams JL, Meng QJ, Streuli CH. Circadian clocks and breast cancer. Breast cancer research: BCR. 2016;18(1):89. Epub 2016/09/02. pmid:27590298.
- 18. Chen ST, Choo KB, Hou MF, Yeh KT, Kuo SJ, Chang JG. Deregulated expression of the PER1, PER2 and PER3 genes in breast cancers. Carcinogenesis. 2005;26(7):1241–6. Epub 2005/03/24. pmid:15790588.
- 19. Tarone RE. Circadian disruption and breast cancer. Epidemiology. 2005;16(5):710–1; author reply 1. pmid:16135953.
- 20. Kuo SJ, Chen ST, Yeh KT, Hou MF, Chang YS, Hsu NC, et al. Disturbance of circadian gene expression in breast cancer. Virchows Arch. 2009;454(4):467–74. Epub 2009/03/19. pmid:19296127.
- 21. Storch KF, Lipan O, Leykin I, Viswanathan N, Davis FC, Wong WH, et al. Extensive and divergent circadian gene expression in liver and heart. Nature. 2002;417(6884):78–83. Epub 2002/04/21. pmid:11967526.
- 22. Lowrey PL, Takahashi JS. Mammalian circadian biology: elucidating genome-wide levels of temporal organization. Annu Rev Genomics Hum Genet. 2004;5:407–41. pmid:15485355.
- 23. Hoffman AE, Yi CH, Zheng T, Stevens RG, Leaderer D, Zhang Y, et al. CLOCK in breast tumorigenesis: genetic, epigenetic, and transcriptional profiling analyses. Cancer Res. 2010;70(4):1459–68. Epub 2010/02/02. pmid:20124474.
- 24. Fu A, Leaderer D, Zheng T, Hoffman AE, Stevens RG, Zhu Y. Genetic and epigenetic associations of circadian gene TIMELESS and breast cancer risk. Mol Carcinog. 2012;51(12):923–9. Epub 2011/10/17. pmid:22006848.
- 25. Winter SL, Bosnoyan-Collins L, Pinnaduwage D, Andrulis IL. Expression of the circadian clock genes Per1 and Per2 in sporadic and familial breast tumors. Neoplasia. 2007;9(10):797–800. pmid:17971899.
- 26. Hoffman AE, Zheng T, Yi CH, Stevens RG, Ba Y, Zhang Y, et al. The core circadian gene Cryptochrome 2 influences breast cancer risk, possibly by mediating hormone signaling. Cancer prevention research. 2010;3(4):539–48. Epub 2010/03/16. pmid:20233903.
- 27. Yi C, Mu L, de la Longrais IA, Sochirca O, Arisio R, Yu H, et al. The circadian gene NPAS2 is a novel prognostic biomarker for breast cancer. Breast Cancer Res Treat. 2010;120(3):663–9. Epub 2009/08/01. pmid:19649706.
- 28. Cadenas C, van de Sandt L, Edlund K, Lohr M, Hellwig B, Marchan R, et al. Loss of circadian clock gene expression is associated with tumor progression in breast cancer. Cell Cycle. 2014;13(20):3282–91. Epub 2014/12/09. pmid:25485508.
- 29. Climent J, Perez-Losada J, Quigley DA, Kim IJ, Delrosario R, Jen KY, et al. Deletion of the PER3 gene on chromosome 1p36 in recurrent ER-positive breast cancer. J Clin Oncol. 2010;28(23):3770–8. Epub 2010/07/12. pmid:20625127.
- 30. Fu L, Pelicano H, Liu J, Huang P, Lee C. The circadian gene Period2 plays an important role in tumor suppression and DNA damage response in vivo. Cell. 2002;111(1):41–50. pmid:12372299.
- 31. Gery S, Komatsu N, Baldjyan L, Yu A, Koo D, Koeffler HP. The circadian gene per1 plays an important role in cell growth and DNA damage control in human cancer cells. Mol Cell. 2006;22(3):375–82. pmid:16678109.
- 32. Hua H, Wang Y, Wan C, Liu Y, Zhu B, Yang C, et al. Circadian gene mPer2 overexpression induces cancer cell apoptosis. Cancer Sci. 2006;97(7):589–96. pmid:16827798.
- 33. Gery S, Virk RK, Chumakov K, Yu A, Koeffler HP. The clock gene Per2 links the circadian system to the estrogen receptor. Oncogene. 2007;26(57):7916–20. Epub 2007/06/18. pmid:17599055.
- 34. Xiao L, Chang AK, Zang MX, Bi H, Li S, Wang M, et al. Induction of the CLOCK gene by E2-ERα signaling promotes the proliferation of breast cancer cells. PLoS One. 2014;9(5):e95878. Epub 2014/05/02. pmid:24789043.
- 35. Mao Y, Fu A, Leaderer D, Zheng T, Chen K, Zhu Y. Potential cancer-related role of circadian gene TIMELESS suggested by expression profiling and in vitro analyses. BMC Cancer. 2013;13:498. Epub 2013/10/25. pmid:24161199.
- 36. Chi L, Zou Y, Qin L, Ma W, Hao Y, Tang Y, et al. TIMELESS contributes to the progression of breast cancer through activation of MYC. Breast cancer research: BCR. 2017;19(1):53. Epub 2017/05/04. pmid:28464854.