Next Article in Journal
Optimization of ZnAl/Chitosan Supra-Nano Hybrid Preparation as Efficient Antibacterial Material
Next Article in Special Issue
Molecular Interactions between Pathogens and the Circadian Clock
Previous Article in Journal
rAAV Engineering for Capsid-Protein Enzyme Insertions and Mosaicism Reveals Resilience to Mutational, Structural and Thermal Perturbations
Previous Article in Special Issue
The Mechanisms Underlying the Cytotoxic Effects of Copper Via Differentiated Embryonic Chondrocyte Gene 1
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Circadian Gene Polymorphisms Associated with Breast Cancer Susceptibility

1
Department of Molecular Genetics and Epigenetics, Nofer Institute of Occupational Medicine, 91-348 Lodz, Poland
2
Department of Environmental Epidemiology, Nofer Institute of Occupational Medicine, 91-348 Lodz, Poland
3
Department of Biological and Environmental Monitoring, Nofer Institute of Occupational Medicine, 91-348 Lodz, Poland
4
Department of Medical Laboratory Diagnostics and Bank of Frozen Tissues and Genetic Specimens, Medical University of Gdansk, 80-211 Gdansk, Poland
5
Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.PL), 80-211 Gdansk, Poland
6
Department of Surgical Oncology, Medical University of Gdansk, 80-214 Gdansk, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(22), 5704; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20225704
Submission received: 21 October 2019 / Revised: 10 November 2019 / Accepted: 11 November 2019 / Published: 14 November 2019
(This article belongs to the Special Issue Crosstalk between Circadian Rhythm and Diseases)

Abstract

:
Breast cancer (BC) is a major problem for civilization, manifested by continuously increasing morbidity and mortality among women worldwide. Core circadian genes may play an important role in cancer development and progression. To evaluate the effects of single nucleotide polymorphism (SNP) in circadian genes in BC risk, 16 functional SNPs were genotyped in 321 BC patients and 364 healthy women using the TaqMan fluorescence-labelled probes or High-Resolution Melt Curve technique in the Real-Time PCR system. The selected SNPs were analyzed for the risk of BC, progression, and the influence on gene expression in BC tissue pairs to demonstrate the functionality of genetic variants. The study showed a relationship between an increased BC risk under the dominant genetic model of CRY2 rs10838524, PER2 rs934945, and recessive genetic model of PER1 rs2735611. A protective effect of BMAL1 rs2279287 was observed among carriers with at least one variant allele. Moreover, we found an increased risk of estrogen-/progesterone-positive tumors under the dominant genetic model of PER2 rs934945 and estrogen negative tumors under the variant genotype of CRY2 rs10838524, PER1 rs2735611. We demonstrated significantly altered gene expression of BMAL1, CRY2, PER1, PER2, PER3 according to particular genotypes in the BC tissue pairs. Our findings support the hypothesized role of circadian genes in breast carcinogenesis and indicate probable biomarkers for breast cancer susceptibility.

1. Introduction

Despite significant progress in early diagnostics and modern therapies of BC in the last decades, breast cancer (BC) still constitutes the most common malignancy among women worldwide. According to the latest report of Global Cancer Statistics 2018 (GLOBOCAN 2018), breast cancer is ranked first place as a deadly disease among female cancer patients [1]. As the scientific environment is trying to find a new approach to elucidate the mechanism of mammary tumorigenesis, new biomarkers are needed to extend the diagnostics.
Single-nucleotide polymorphism (SNP) is the most common genetic variant in the human genome, which is considered a stable biomarker of genetic background to predict risk, progression, and treatment response to various diseases [2]. Researchers have revealed that several SNPs in circadian genes are associated with breast cancer susceptibility [3,4]. SNPs are located in different gene bodies including promoters, exons, introns as well as 5′- and 3′ UTRs [5]. Therefore, changes in gene expression and their cancer predisposition differ depending on the location of SNPs. The location of SNPs may affect gene expression by altering promoter activity, binding transcription factors, as well as DNA CpG sites methylation [2,6,7]. Additionally, cancer risk may depend on exonal SNPs by suppressing gene transcriptions and translations. SNPs in introns region are also not indifferent to gene functions. Such a location may generate splice variants of transcripts and promote or disrupt binding and function of long non-coding RNAs. SNPs in the 5′UTR affect translation, whereas SNPs in the 3′UTR affect binding of microRNA [2].
So far, there has not been enough research in relation to the incidence of breast cancer and the polymorphic variants of the circadian genes. Polymorphisms of positive BMAL1, CLOCK, NPAS2, and negative regulators of circadian rhythms PERIOD (1,2,3), CRYPTOCHROME (1,2), as well as TIMELESS have been most commonly studied. Only in individual epidemiological studies, statistically significant correlations between circadian gene variants and breast cancer have been demonstrated. Some of the analyzed SNPs have also been linked to other cancers and diseases [3]. Therefore, we conducted an association study among 321 newly diagnosed BC patients and 364 healthy women, living in the same longitude of Poland, in Gdansk and Lodz. The main objective of our project was to investigate an association between genetic variants of crucial circadian genes CLOCK, NPAS2, BMAL1, PER1, 2, 3, CRY1, 2, and TIMELESS and the risk of BC, progression, and the influence of gene expression on BC tissue pairs to demonstrate their functional significance in the process of carcinogenesis in the mammary gland.

2. Results

The analysis was based on the cases of primary breast cancer patients and healthy volunteers, all of Caucasian ancestry. The selected characteristics of the study population are presented in Table 1. Among breast cancer patients, there were more non-smoking and pre-menopausal women compared to the control group. The women from both groups were at a similar age. The mean age for the breast cancer patients was 58.85 (SD 11.50) and for the healthy subject it was 60.80 (SD 7.11), p = 0.007. Other demographic characteristics did not differ significantly between the cases and controls. A total of 16 preselected SNPs in nine core circadian genes were genotyped and most of them did not departure from the Hardy–Weinberg equilibrium (HWE) except for three SNPs, which were not in HWE, namely: CLOCK rs12505266, PER1 rs3027178, PER3 rs2640909 in the control group (Table 2).

2.1. Selected Circadian Gene Polymorphism Is Associated with Breast Cancer Risk

To determine whether the particular genotypes of core circadian genes could influence breast cancer risk, we performed genotyping analysis using TaqMan probes or High-Resolution Melting technique on DNA of 107 matched pairs tissues and 214 blood samples from breast cancer patients.
Of the 16 investigated SNPs five were found to be significantly associated with breast cancer risk. Rs2279287 is located in the 5′ UTR region of BMAL1, according to SNP info NIH platform (https://snpinfo.niehs.nih.gov/snpfunc.htm, Table 2). This polymorphic site is the probable point of attachment of transcription factors. The carriers of the minor allele (T) rs2279287 had a reduced predisposition to breast cancer assuming a recessive homozygous genotype OR = 0.47 (0.27–0.80) p = 0.02 and a dominant genetic model OR = 0.69 (0.50–0.95) p = 0.02. Additionally, we also found a potential protective effect (at marginal statistical significance) of minor allele (G) rs3027178 for PER1 assuming a recessive genotype OR = 0.54 (0.28–1.08); p = 0.09 and a recessive genetic model OR = 0.55 (0.28–1.08); p = 0.08 (Table 3).
The women who had one or more protective alleles (derived at least one gene BMAL1 rs2279287 and/or PER1 rs3027178) had a significantly reduced breast cancer risk. OR = 0.49 (0.32–0.77) p = 0.002 (Table 4).
Rs10838524 is located on an intron of the CRY2 locus. The present analysis suggested that minor allele (G) is associated at marginal statistical significance with increased breast cancer susceptibility of 65% under a recessive genotype OR = 1.65 (1.05–2.58), p = 0.07 and significantly of 45% under a dominant genetic model OR = 1.45 (1.00–2.10), p = 0.05.
Intronic rs12505266 of CLOCK was significantly associated with an increased predisposition to breast cancer among recessive homozygous genotypes OR = 1.22 (0.83–1.79), p = 0.057 and at marginal significance under a recessive genetic model OR = 1.38 (0.97–1.96), p = 0.07. Similar associations were observed between PER2 missense variant rs934945 under a dominant genetic model and an increased breast cancer risk OR = 1.56 (1.09–2.23), p = 0.01. Potential significance of an increased predisposition to breast cancer was demonstrated for a recessive genetic model of PER2 rs11894491 OR = 1.53 (0.93–2.50) p = 0.09. We found a significant association between a heterozygous genotype of TIMELESS rs2279665 and breast cancer risk OR = 0.69 (0.481.00) p = 0.02. For another SNP we found marginal significance including a recessive homozygous genotype and a recessive genetic model for PER1 rs3027178 OR = 0.54 (0.28–1.08) p = 0.09; OR = 0.55 (0.28–1.08) p = 0.08, respectively (Table 3).
The patients having at least one or more risk alleles (among three significant SNPs CRY2 rs10838524; PER1 rs2735611; PER2 rs934945) were significantly associated with an increased breast cancer risk OR = 1.66 (1.17–2.35) p = 0.005 (Table 4).

2.2. Circadian Gene Variants Are Associated with an Estrogen and Progesterone Receptor Status

In addition, to the main effects model, the control group was simultaneously compared to the breast cancer patients stratified into two groups according to their hormonal receptor status-estrogen-positive (ER+)/progesterone-positive (PR+) tumors and estrogen-negative (ER−)/progesterone-negative (PR−) tumors. The genotyping analysis were performed on DNA samples derived from 107 matched tissue pairs and 214 blood samples from breast cancer patients. Interestingly, four of 16 were significantly associated with a breast cancer risk, including two SNPs with ER+/PR+ tumors and two SNPs with both types of tumors. Recessive homozygotes of BMAL1 rs2279287 had a significant protective effect of ER+/PR+ tumors OR = 0.45 (0.24–0.86) p = 0.04 in comparison to an increased risk effect of PER1 rs2735611 OR = 3.76 (1.38–10.24) p = 0.01 and PER2 rs934945 OR = 1.75 (1.16–2.63) p = 0.01. While recessive homozygotes of CRY2 rs10838524 OR = 2.26 (1.03–4.94) p = 0.03 and PER1 OR = 6.08 (1.80–20.56) p = 0.003 were significantly associated with ER-/PR- tumors Table 5.

2.3. A Putative Functional Effect of SNP on Circadian Gene Expression

To determine whether the genotypes of SNP could alter mRNA transcription we compared data obtained from the gene expression analysis using the quantitative real-time PCR among 107 matched pairs tissues from breast cancer patients to their genotype. As shown in Table 6, we found that mRNA expression level was significantly changed according to particular gene variants, indicating that SNPs may affect transcription of these genes. The carriers who had at least one variant allele of BMAL1 rs2279297 (p = 0.05) or were a recessive homozygote of CRY2 rs3824872 (p = 0.0004), had a significantly decreased gene expression level in breast cancer tissue compared to the women with a dominant homozygous genotype. In contrast, the carriers of a dominant genotype of PER3 rs10462020 had a significantly reduced level of mRNA expression in breast cancer tissue p = 0.02. We also observed significant differences in a transcript level in adjacent non-tumor tissues. The carriers with a variant allele had reduced gene expression of PER1 rs2735611, (p = 0.04), PER1 rs3027178 (p = 0.002), or an increased transcript level of PER2 rs11894491 (p = 0.03). The bioinformatic analysis (GTex eQTL calculator, www.gtexportal.org) was employed to explore potential biological effects of selected SNPs of core circadian genes on their mRNA expression in normal breast-mammary tissue [8]. As shown in Figure 1, we observed that three CLOCK rs12505266, p = 0.000048; NPAS2 rs2305160, p = 0.017; PER2 rs11894491, p = 0.036 of the investigated genes had a significantly different gene expression pattern in breast mammary tissue according to different genotypes and two SNPs including CLOCK rs1801260 p = 0.065 and CRY1 rs8192440 p = 0.069 had a different expression pattern at the marginal statistical significance (Figure 1). Only one SNP-rs11894491-showed a significant effect on a PER2 transcript level in mammary cancer-free tissues, based on the results from our dataset and the GTex eQTL calculator. However, we identified higher gene expression in the carriers with a minor genotype in comparison to GTex platform, where the carriers with a major genotype had a higher mRNA level than the carriers with other gene variants.

3. Discussion

Our analysis indicates that five SNPs, BMAL1 rs2279287, CLOCK rs12505266, CRY2 rs10838524, PER1 rs2735611, PER2 rs934945, are significantly associated with the breast cancer risk in the Polish association study. Additionally, we observed differences in the gene expression level according to particular gene variants: BMAL1 rs2279287, CRY2 rs10838524, PER1 rs2735611, rs3027178, PER2 rs11894491, PER3 rs10462020.
There are several studies demonstrating that alterations of circadian rhythm together with disruptions caused by shift work and exposure to light at night (LAN) may promote breast cancer development [9,10]. There are a few studies on circadian gene variants in shift workers and susceptibility to breast carcinogenesis [11,12]. Molecular epidemiological studies reveal a significant association of the polymorphic variants of circadian genes and the risk of various cancer types such as breast cancer [13,14,15,16,17], colorectal cancer [18,19,20], glioma [21], prostate cancer [22,23], hepatocellular carcinoma [24,25], lung cancer [26], non-Hodgkin’s lymphoma [27], as well as mood [28,29] and metabolic disorders [30].
Genetic variants such as SNPs play a substantial role in gene expression regulation, mRNA translation, and degradation protein structures all of which may affect gene functions and human phenotype [2]. Circadian genes play essential roles in regulation of gene expression, including cell proliferation, DNA damage repair, cell cycle control, and apoptosis. Due to their physiological role, alteration of these genes may lead to cancer development [3,4]. Considering this fact, it is sensible to assume that SNPs of circadian genes may affect cancer susceptibility, cancer cell proliferation, invasion, as well as a response to treatment and even patients’ survival [31,32,33]. Consequently, explaining the mechanism of circadian genes in cancer development will be useful in better diagnostics and improving patient’s healthcare.
In the present study, we described findings of the first approach, to the best of our knowledge, which has investigated the relationship between various polymorphic sites of circadian genes and the transcript level among breast cancer patients. The results were based on genotyping from 364 healthy women and 321 breast cancer patients. Additionally, genotyping results were analyzed according to gene expression data obtained from 107 breast cancer tissue pairs. We hypothesized that circadian gene SNPs may influence their expression level and susceptibility to breast cancer. Our results support this hypothesis. In fact, five of the 16 analyzed genetic variants were statistically significantly associated with a predisposition to breast cancer employing a dominant or recessive genetic model of inheritance. Up to now, some significant genetic variants of circadian genes have been recognized as associated with breast cancer risk [13,14,16,34,35,36,37,38].
CLOCK rs1801260 located on the 3′ UTR region, is one of the most common variants studied by many authors with regard to malignancies as well as other disorders [3,4]. Regarding investigations of rs1801260 and breast cancer predisposition, Hoffman et al. have indicated that under the additive genetic model, there is an increased risk of ER/PR negative breast cancer OR = 2.57 (1.14–5.82) [34]. Subsequently, CECILE study has shown that rs1801260 is linked with breast cancer risk among postmenopausal women [16]. According to SNPinfo Web Server (http://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html), rs1801260 is predicted as a microRNA binding site. Further, a performed bioinformatic analysis microRNA.org—Targets and Expression (http://www.microrna.org/microrna/home.do) has revealed that rs1801260 is within the miRNA hsa-miR-141 that was identified as a probable diagnostic marker of breast cancer [39]. However, our analysis did not show a significant association between CLOCK and breast cancer risk. What is more, we did not demonstrate significant changes at a mRNA level under this CLOCK variant. We selected another rs12505266 CLOCK for the analysis. We found, at marginal significance, that the carriers of a minor allele T demonstrated an increased breast cancer risk in comparison to the healthy population OR = 1.21 (0.82–1.78) p = 0.057. At a marginal significant level, the carriers with TT genotype had a higher risk of developing ER+ tumors than ER-, OR = 1.28 (0.82–2.01). This is an intronic polymorphic site, but according to SNPinfo Web Server (http://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html), rs12505266 is predicted as a transcription factor binding site of CLOCK, which may regulate expression of CLOCK gene. Furthermore, the in-silico analysis using the GTEx eQTL calculator (https://gtexportal.org/home/testyourown) [8], has shown a significantly elevated gene expression level in breast mammary tissues carrying variant-containing TT genotypes. Nevertheless, our gene expression analysis showed changes at marginal statistical significance on a mRNA level under rs12505266 CLOCK in adjacent non-tumor tissues, where tissues with a minor genotype had an increased transcript level of CLOCK. This inconsistency may be due to a smaller size population used in our study. Other epidemiological studies have indicated a negative effect of intronic rs3805151 [37], rs11932595, rs7698022, and 3′UTR SNPs rs1048004, rs1801260 [16] CLOCK polymorphism as well as a positive effect under the dominant genetic model of CLOCK rs6850524, rs13102385, rs11133391 on breast cancer susceptibility [34].
Of particular interest, the SNP rs2279287 within the BMAL1 gene was found to be associated with a reduced breast cancer risk, among the carriers with the variant containing CT/TT genotypes than among those with the homozygous wild-type genotype (CC). A reduced risk of breast tumorigenesis was also observed in ER+ tumors compared to ER- in carriers of at least one variant allele. The polymorphic site rs2279287 as predicted on SNPinfo Web Server http://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html is the probable point of transcription factor attachment, which may have an influence on BMAL1 gene expression. Indeed, we found a significant relationship between BMAL1 genetic polymorphism (rs2279287) and mRNA transcription. Breast cancer tissues carrying a variant containing CT/TT genotypes had a reduced level of BMAL1 expression. In the literature, we found one publication with significant results where rs2279287 was associated with a seasonal affective disorder in Indian families [40].
Moreover, there are a few SNPs which were distinguished in Period and Cryptochrome genes. Among circadian genes, PERs are recognized as tumor suppressor genes, which is relevant for the process of carcinogenesis [41]. Studies on animals have shown that PER2 knockout rodents are more prone to cancer initiation after gamma radiation in comparison to wild type mice [42,43]. Moreover, overexpression of PER2 in MCF-7 leads to significant inhibition of breast cancer cell proliferation [41]. Therefore, most of the investigated SNPs referred to Period genes (PER1 rs2735611, PER1 rs3027178, PER2 rs11894491, PER2 rs2304672, PER2 rs934945, PER3 rs10462020, PER3 rs2640909) and were incorporated into the analysis. Selected SNPs according to SNPinfo Web Server (http://snpinfo.niehs.nih.gov/snpinfo/snpfunc.html) are predicted as an exonic splicing enhancer or silencer (ESE or ESS), which may have a role in gene expression regulation. Rs934945 is located on the last exon of PER2 and has a missense functional effect leading to the submission of Glycine-Glutamic Acid. According to NCBI Protein database, the rs934945 site is also a CRY-binding domain. Carriers who had at least one minor allele (T) had an increased breast cancer risk employing a dominant genetic model. Similar observations have been found by Dai et al., i.e., patients carrying both the PER2 (rs934945) TT and CLOCK (rs3805151) CC genotypes had an increased breast cancer risk [37]. Opposite results have been demonstrated by Benna et al. showing that patients with recessive genotype PER2 (rs934945) TT had a reduced predisposition to sarcoma and liposarcoma [32]. Additionally, other research has not indicated associations between PER2 (rs934945) and gastric cancer [31]. Interestingly, our research partially confirmed previous findings of Benna et al. showing that PER1 rs3027178 had a protective effect and reduced risk of soft tissue sarcoma [32]. We observed a similar effect at the statistical significance margin under the recessive genetic model. The carriers with at least one G allele may have a reduced breast cancer risk of 45%. Opposite findings have been presented by Zhang et al. where carriers with at least one minor allele (PER1 rs3027178) among hepatocellular carcinoma patients, had a worse prognosis of recurrence-free survival [24]. A later study of Qu et al. has shown that PER1 rs3027178 was also associated with worse overall survival of gastric cancer patients under a dominant genetic model [31]. In our study, the strongest identified association was rs2735611 in the exon 18 of PER1. The carriers of recessive genotype (AA) had an increased risk of breast cancer OR = 4.37 (1.78–10.69) p = 0.001. According to the results of gene expression, the carriers with at least one variant allele had decreased mRNA expression of PER1. The decreased expression level of PER1 was observed in breast cancer tissue and associated with worse breast cancer type ER-. Among gastric cancer patients this association has not been demonstrated [31].
Together with Period genes, CRY1 and CRY2 form a negative arm of the circadian feedback loop and are essential for the proper functioning of circadian machinery. These genes have potential to influence expression of many biological pathways directly or indirectly [44]. For instance, knockdown of CRY2 in MCF-7 cells increased accumulation of mutagen-induced DNA damage and altered expression of genes involved in the DNA damage response and cell cycle regulatory pathways including cyclin-dependent kinase inhibitor p21 (CDKN1A) and induction of cyclin D1 (CCND1)—an oncogene that is often overexpressed in breast cancer cases [45]. There are a few studies which focus on CRY1 and CRY2 polymorphism, but there are some significant gene variants associated with breast cancer risk [12,37]. In our analysis we found only one significant association among Cryptochrome genes—CRY2 rs10838524. The patients with at least one variant allele had an increased breast cancer risk. Moreover, a minor homozygous genotype was significantly associated with breast cancer risk in the women with ER- tumors. There are no other publications where this genotype is associated with an increased cancer predisposition, but CRY2 rs10838524 was associated with winter depression among Finnish and Swedish populations [46]. Other cryptochrome gene variants have previously been shown to be associated with a decreased breast cancer risk including GT genotype of rs1056560 in CRY1 (OR = 0.84 95% CI = 0.71–0.99) and CC genotype of rs1401417 in CRY2 (OR = 0.24 95% CI = 0.08–0.73) in premenopausal women [37]. However, a significant association with breast cancer risk in postmenopausal patients for three SNPs in CRY2 has been found. GC or CC genotypes of rs1401417 (OR = 0.44 95% CI = 0.21–0.92), AG or GG genotypes of rs11038689, (OR = 0.71 95% CI = 0.51–0.99) and AA genotype of rs7123390, (OR = 0.44 95% CI = 0.22–0.86) [35]. Among Han Chinese gastric cancer patients, SNP rs1056560 in the CRY1, under an additive model had a protective effect on the overall survival of the patients with a HR of 0.72 (95% CI 0.58–0.88, p = 0.021) [31].
Another polymorphic gene with potential relevance in breast cancer development is NPAS2. A first case-control study on non-synonymous polymorphisms rs2305160 Ala394Thr of NPAS2, among 431 breast cancer women and 476 controls (most of them were Caucasians) has revealed that rs2305160 Ala/Thr genotype had a significant association with a decreased breast cancer risk in comparison to more frequent Ala/Ala genotype (OR = 0.61 95% CI = 0.46–0.81) [13]. The interesting issue is that considering the menopausal status, both groups, pre- and post-menopausal patients, had a significant relationship between heterozygous genotype and breast cancer risk, but stronger association has been observed among premenopausal women (OR = 0.44 95% CI = 0.25–0.77, p = 0.004) than in postmenopausal women (OR = 0.65 95% CI = 0.46–0.91) [13]. Further studies on Ala394Thr on 348 breast cancer tissue using the TaqMan allelic discrimination assay have shown a borderline, but non-significant association between homozygous Thr/Thr genotype and a poor prognosis in breast cancer survival [14]. A later study of Yuan et al. has demonstrated a significant association between rs2305160 in NPAS2 and overall death risk in hepatocellular carcinoma patients under genetic dominant model (HR = 1.63 (1.29–2.07), p < 0.001) [25]. Our data analysis did not show a significant effect of this polymorphism on breast cancer among the study group. Only the bioinformatic analysis (based on GTex eQTL calculator) revealed significant differences in mRNA level according to the genotypes of rs2305160.
Among genetic association studies on circadian genes and their impact on cancer development, only one publication has demonstrated a significant relationship between polymorphic variants of TIMELESS (rs7302060, rs2291738) and breast cancer risk [36]. This study has reported a significant association between both tagging SNPs rs7302060, rs2291738 and breast cancer, showing that the C allele ssof rs7302060 has been linked to reduction of breast cancer risk (OR = 0.54 95% CI = 0.54–0.99). A similar observation was made for two genotypes CC of rs730260 (OR = 0.35 95% CI = 0.16–0.78) and GG of rs2291738 (OR = 0.45 95% CI = 0.217–0.97). Both had an association with a decreased breast cancer risk among ER/PR positive tumors [36]. Our analysis revealed that a heterozygous genotype of rs2279665 is linked to a reduced predisposition to breast cancer risk also in the case of ER+/PR+ tumors.
Moreover, the recent GWAS study conducted by Mocellin et al. has shown a highly significant association between the breast cancer risk and circadian pathway variation (p-value 1.9 × 10−6) The top genes and number of SNPs involved in cancer development were RORA and RORB (eight SNPs and five SNPs, respectively); other significant genes with single SNPs were PER1 (two SNPs), ARNTL, CRY2, CLOCK, CRY1, and RORC, in order of significance. A similar observation has been also found for ER- breast tumors [27].
A potential limitation of our study was the relatively small size of the study group. Additionally, we were not able to perform prognostic analysis stratified by circadian SNPs, according to patients’ overall survival because of the lack of relevant data. Therefore, associations between circadian gene polymorphisms and BC should be treated with caution and validated on a larger population.

4. Methods and Materials

4.1. Study Population and Biological Materials

A total of 321 Polish women, all of Caucasian origin, with a newly and histologically diagnosed BC and without previous chemotherapy were included in this association study. The patients were derived from two medical centers in Poland: Clinic of Oncologic Surgery at the Medical University of Gdansk and the Copernicus Memorial Hospital, Lodz, Poland. The breast cancer patients were hospitalized over the years 2006–2015. The study group consisted of breast cancer patients who were consecutively admitted to the Clinics and agreed to donate the biological material for scientific purposes. Biological material consisted of 107 BC tissue pairs tumor tissue and tumor-adjacent (cancer-free) tissue (0.6 × 0.6 cm from the same patient) and 214 blood samples (7.5 mL). Matched BC tissue samples were derived from a group of patients hospitalized at Clinic of Oncologic Surgery at the Medical University of Gdansk while blood samples were derived from a group of patients hospitalized at the Copernicus Memorial Hospital. All the tissue specimens were collected during the tumor removal surgery and evaluated by pathologists to confirm the diagnosis of both tumor and adjacent tissue (the second one was collected close to the resection margin). Then the samples of BC tissue pairs were frozen at −70 °C and deposited in Bank of Frozen Tissues and Genetic Specimens at Medical University of Gdansk for a further molecular analysis. The blood samples (7.5 mL) were collected from the patients into heparinized test tubes and separated by centrifugation into a buffy coat for DNA isolation. Each fraction was stored at −20 °C until a molecular analysis at NIOM.
The control group consisted of 364 of healthy women who were gathered for projects carried out at NIOM between 2008 and 2014. Eligible controls at above 50 years of age and with lack of shift work history were selected from general population, who voluntarily agreed to donate biological material for scientific purposes. Blood samples (7.5 mL) were collected from each healthy volunteer into heparinized test tubes and separated by centrifugation into the buffy coat for DNA isolation. Each fraction was stored at −20 °C until a molecular analysis at NIOM. The women from the control group who were reported to work in shifts (at the time of recruitment) were not included in the study, because shiftwork that involves circadian disruption has been classified by the International Agency for Research on Cancer as a probable human carcinogen group 2A [47,48].
Basic epidemiological characteristics (age, BMI, smoking, and menopausal status) were collected using individual questionnaires (from the BC patients and the control group), whereas clinical data (receptor status, tumor stage and grade, histological type of tumor) were obtained from medical records of both groups of the BC patients (Table 1).

4.2. Ethics Declaration

This study was conducted in compliance with the Declaration of Helsinki and was approved by the relevant Local Ethics Committees (Independent Ethics Committee at Medical University of Gdansk, resolution No. NKEBN/781/2005, on 15 December 2005 and Ethical Institutional Review Board at the Nofer Institute of Occupational Medicine, resolution No. 4/2019; on 4 April 2019. We obtained a written consent of the participants to deposit their biological material in the Bank of Frozen Tissues and Genetic Specimens of the Medical University of Gdansk (tissue specimens) and in the NIOM (blood specimens) to use them in further analyses.

4.3. DNA Isolation Procedures

Total genomic DNA from the tissue pairs (adjacent non-tumor breast tissue and breast cancer tissue) and buffy coat were isolated by the DNeasy Blood and Tissue KIT (Qiagen, Hilden, Germany) according to the manufacturer’s protocol in order to evaluate SNPs in circadian genes. The amount of the total extracted DNA was determined in microplates using the spectrophotometric method (MultiscanTM Go Microplate Spectrophotometer, Thermo Fisher; Waltham, MA, USA).

4.4. SNP Selection, Genotyping, and Gene Expression.

Potential functional SNPs that were included had to meet the following criteria: SNPs of core circadian genes were selected based on their position within the gene, predicted function according to SNPinfo Web Server (https://snpinfo.niehs.nih.gov/snpfunc.htm); validated SNPs with minor allele >10% in the Caucasian population according to the SNP base in the National Center for Biotechnology Information, https://www.ncbi.nlm.nih.gov, and on the basis of relevant literature [3,4]. Finally, a total of 16 functional SNPs in nine core circadian genes were selected including non-synonymous variation, potential transcription factor binding site, potential splice site, or potential miRNA- binding site (Table 2).
Allelic discrimination for CRY2 (rs10838524), PER3 (rs10462020), PER1 (rs2735611), PER2 (rs934945), PER2 (rs2304672), NPAS2 (rs2305160), TIMELESS (rs2279665) was conducted using the TaqMan SNP Genotyping Assays (Thermo Fisher Scientific, Waltham, MA, USA) and for genes CLOCK (rs12505266), CRY2 (rs3824872) PER1 (rs3027178) PER2 (rs11894491) PER3 (rs2640909) TIMELESS (rs774047) we used the TaqMan probes and primers (Eurofins Genomics, Ebersberg, Germany). Genotyping was performed using the Real-Time PCR method and the Light Cycler 96 Roche (Roche, Basel Switzerland) and FastStart DNA probes master mix (Roche, Basel Switzerland). PCR reactions were carried out with 10 ng of DNA in a final volume of 10 μL. For genes CLOCK (rs1801260), BMAL1 (rs2279287), CRY1 (rs8192440) we used the High-Resolution Melt Curve technique and restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR) for genotype verification. Primers were designed with Beacon Designer 7.01 (PREMIER Biosoft International, Palo alto, CA, USA) according to the GenBank® genetic sequence database. HRM assay was performed on Light Cycler 96 Roche (Roche, Basel, Switzerland) using 10 ng of DNA in a final volume of 20 μL and SsoFastTM EvaGreen Supermix (Bio-Rad Laboratories, Inc, Hercules, CA, USA). Genotyping was performed using the control DNA samples with known genotypes and a negative control.

4.5. Gene Expression

Gene expression data obtained from the 107 tissue pairs of breast cancer patients were compared with the results of allelic discrimination in order to check the influence of selected SNPs on the transcriptome of particular circadian genes. Gene expression assay has been described in detail in our previous work [47] Our results from allelic discrimination were also compared with The Genotype-Tissue Expression (GTEx) using the eQTL calculator (https://gtexportal.org/home/testyourown) [8]. Results of the in silico analysis were obtained using the following data source: GTEx Analysis Release V8 (dbGaP Accession phs000424.v8.p2) on 14 October, 2019. The clinical material used in GTEx platform consisted of 396 Breast Mammary Tissues from the women between 20 and 79 years of age and of mixed ethnicity. GTEx platform does not provide eQTL analysis for cancer tissues, therefore we could not perform in silico analysis for BC tissues.

4.6. Statistics Analysis

Means with standard deviation and frequencies of the basic features were calculated. Normality of the data was evaluated with the Shapiro–Wilk test. The Chi-square test was performed to determine any discrepancies of distribution from the Hardy–Weinberg Equilibrium (HWE). Odds ratios (ORs) and their corresponding 95% confidence intervals (CI) were calculated by the logistic regression analysis in order to measure the strength of association of selected circadian genes polymorphisms and their interactions with breast cancer risk. The p-value was adjusted by confounding variables (age, BMI, smoking status, menopausal status). Additionally, two genetic models: dominant genetic model (genotypes with a presence of at least one copy of the minor allele vs. dominant homozygous genotype) and recessive genetic model (recessive homozygous genotype vs. genotypes with a presence of at least one copy of the major allele) were considered. Mean Normalized Expression (MNE) was calculated previously [49]. The obtained data from gene expression were analyzed according to the best fitting genetic models using the Student’s t-test. All the performed analyses were calculated using STATISTICA 12 software package (Statsoft, Tulsa, OK, USA).

5. Conclusions

Our association study indicated that circadian gene variants, including BMAL1 rs2279287, CLOCK rs12505266, CRY2 rs10838524, PER1 rs2735611, PER2 rs934945, rs11894491, may be associated with breast carcinogenesis. This relationship seems to have a reflection in a biological effect because polymorphisms may influence gene expression, protein function like protein-protein interactions, or mRNA stability. This is consistent with the results obtained in our previous studies [49,50]. However, due to the limited size of the group, the association between functional SNPs in circadian genes and breast cancer susceptibility needs to be interpreted with caution and validated in a larger population. Finally, the majority of breast cancer association studies and circadian genes SNPs were conducted among a Caucasian population. Importantly, it should be remembered, that variability of alleles exists in populations, so allele variants of specific genes are frequently observed in specific populations like Caucasian or Asian populations.

Author Contributions

Conceptualization, M.L. and E.R.; Formal analysis, M.L.; Investigation, M.L.; Methodology, M.L.; Project administration, M.L. and E.R.; Resources, B.P., J.G., B.S., L.K. and J.S.; Supervision, E.R.; Writing–original draft, M.L.; Writing–review and editing, M.L., E.J., E.W. and E.R.

Funding

This work was supported by the National Science Center (Grant No. 2014/15/N/NZ5/01671) (to M.L.) and partially supported by Ministry of Science and Higher Education, Poland DIR/WK/2017/01 (to L.K.), Internal grant for statutory activity at Nofer Institute of Occupational Medicine IMP 14.5/2019 (to M.L., E.R., E.W.).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BCBreast cancer
SNPSingle Nucleotide Polymorphism
EREstrogen Receptor
PRProgesterone Receptor
HER2Human Epidermal Growth Factor Receptor 2
GWASGenome-Wide Association Study
OROdds Ratio
SDStandard Deviation
CIConfidence Intervals
HWEHardy–Weinberg Equilibrium
MNEMean Normalized Expression
GTexThe Genotype-Tissue Expression
eQTLExpression quantitative trait loci
HRMHigh Resolution Melt
TFBSTranscription factor binding site
ESSexonic splicing silencer
ESEexonic splicing enhancer
LANlight at night

References

  1. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed]
  2. Deng, N.; Zhou, H.; Fan, H.; Yuan, Y. Single nucleotide polymorphisms and cancer susceptibility. Oncotarget 2017, 8, 110635–110649. [Google Scholar] [CrossRef]
  3. Valenzuela, F.J.; Vera, J.; Venegas, C.; Muñoz, S.; Oyarce, S.; Muñoz, K.; Lagunas, C. Evidences of Polymorphism Associated with Circadian System and Risk of Pathologies: A Review of the Literature. Int. J. Endocrinol. 2016, 2016, 2746909. [Google Scholar] [CrossRef] [PubMed]
  4. Reszka, E.; Przybek, M.; Muurlink, O.; Peplonska, B. Circadian gene variants and breast cancer. Cancer Lett. 2017, 390, 137–145. [Google Scholar] [CrossRef]
  5. The 1000 Genomes Project Consortium. A global reference for human genetic variation. Nature 2015, 526, 68–74. [Google Scholar] [CrossRef] [PubMed]
  6. Wu, H.; Zhang, K.; Gong, P.; Qiao, F.; Wang, L.; Cui, H.; Sui, X.; Gao, J.; Fan, H. A Novel Functional TagSNP Rs7560488 in the DNMT3A1 Promoter Is Associated with Susceptibility to Gastric Cancer by Modulating Promoter Activity. PLoS ONE 2014, 9, e92911. [Google Scholar] [CrossRef] [PubMed]
  7. Huang, Y.T.; VanderWeele, T.J.; Lin, X. Joint Analysis of SNP and Gene Expression data in genetic association studies of complex diseases. Ann. Appl Stat. 2014, 8, 352–376. [Google Scholar] [CrossRef] [PubMed]
  8. Consortium, G.T. Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans. Science 2015, 348, 648–660. [Google Scholar] [CrossRef] [PubMed]
  9. Stevens, R.G.; Brainard, G.C.; Blask, D.E.; Lockley, S.W.; Motta, M.E. Breast cancer and circadian disruption from electric lighting in the modern world. CA Cancer J. Clin. 2014, 64, 207–218. [Google Scholar] [CrossRef]
  10. Bonde, J.P.; Hansen, J.; Kolstad, H.A.; Mikkelsen, S.; Olsen, J.H.; Blask, D.E.; Härmä, M.; Kjuus, H.; de Koning, H.J.; Olsen, J.; et al. Work at night and breast cancer--report on evidence-based options for preventive actions. Scand. J. Work Environ. Health. 2012, 38, 380–390. [Google Scholar] [CrossRef]
  11. Monsees, G.M.; Kraft, P.; Hankinson, S.E.; Hunter, D.J.; Schernhammer, E.S. Circadian genes and breast cancer susceptibility in rotating shift workers. Int. J. Cancer 2012, 131, 2547–2552. [Google Scholar] [CrossRef]
  12. Zienolddiny, S.; Haugen, A.; Lie, J.A.; Kjuus, H.; Anmarkrud, K.H.; Kjærheim, K. Analysis of polymorphisms in the circadian-related genes and breast cancer risk in Norwegian nurses working night shifts. Breast Cancer Res. 2013, 15, R53. [Google Scholar] [CrossRef] [PubMed]
  13. Zhu, Y.; Stevens, R.G.; Leaderer, D.; Hoffman, A.; Holford, T.; Zhang, Y.; Brown, H.N.; Zheng, T. Non-synonymous polymorphisms in the circadian gene NPAS2 and breast cancer risk. Breast Cancer Res. Treat. 2008, 107, 421–425. [Google Scholar] [CrossRef] [PubMed]
  14. Yi, C.; Mu, L.; de la Longrais, I.A.; Sochirca, O.; Arisio, R.; Yu, H.; Hoffman, A.E.; Zhu, Y.; Katsaro, D. The circadian gene NPAS2 is a novel prognostic biomarker for breast cancer. Breast Cancer Res. Treat. 2010, 120, 663–669. [Google Scholar] [CrossRef] [PubMed]
  15. Wirth, M.D.; Burch, J.B.; Hebert, J.R.; Kowtal, P.; Mehrotra-Kapoor, A.; Steck, S.E.; Hurley, T.G.; Gupta, P.C.; Pednekar, M.S.; Youngstedt, S.D.; et al. Case-control study of breast cancer in India: Role of PERIOD3 clock gene length polymorphism and chronotype. Cancer Investig. 2014, 32, 321–329. [Google Scholar] [CrossRef] [PubMed]
  16. Truong, T.; Liquet, B.; Menegaux, F.; Plancoulaine, S.; Laurent-Puig, P.; Mulot, C.; Cordina-Duverger, E.; Sanchez, M.; Arveux, P.; Kerbrat, P.; et al. Breast cancer risk, nightwork, and circadian clock gene polymorphisms. Endocr. Relat. Cancer 2014, 21, 629–638. [Google Scholar] [CrossRef]
  17. Rabstein, S.; Harth, V.; Justenhoven, C.; Pesch, B.; Plöttner, S.; Heinze, E.; Lotz, A.; Baisch, C.; Schiffermann, M.; Brauch, H.; et al. Polymorphisms in circadian genes, night work and breast cancer: Results from the GENICA study. Chronobiol. Int. 2014, 31, 1115–1122. [Google Scholar] [CrossRef]
  18. Zhou, F.; He, X.; Liu, H.; Zhu, Y.; Jin, T.; Chen, C.; Qu, F.; Li, Y.; Bao, G.; Chen, Z.; et al. Functional polymorphisms of circadian positive feedback regulation genes and clinical outcome of Chinese patients with resected colorectal cancer. Cancer 2012, 118, 937–946. [Google Scholar] [CrossRef]
  19. Karantanos, T.; Theodoropoulos, G.; Gazouli, M.; Vaiopoulou, A.; Karantanou, C.; Stravopodis, D.J.; Bramis, K.; Lymperi, M.; Pektasidis, D. Association of the clock genes polymorphisms with colorectal cancer susceptibility. J. Surg. Oncol. 2013, 108, 563–567. [Google Scholar] [CrossRef]
  20. Alexander, M.; Burch, J.B.; Steck, S.E.; Chen, C.F.; Hurley, T.G.; Cavicchia, P.; Ray, M.; Shivappa, N.; Guess, J.; Zhang, H.; et al. Case-control study of the PERIOD3 clock gene length polymorphism and colorectal adenoma formation. Oncol. Rep. 2015, 33, 935–941. [Google Scholar] [CrossRef]
  21. Madden, M.H.; Anic, G.M.; Thompson, R.C.; Nabors, L.B.; Olson, J.J.; Browning, J.E.; Monteiro, A.N.; Egan, K.M. Circadian pathway genes in relation to glioma risk and outcome. Cancer Causes Control. 2014, 25, 25–32. [Google Scholar] [CrossRef] [PubMed]
  22. Markt, S.C.; Valdimarsdottir, U.A.; Shui, I.M.; Sigurdardottir, L.G.; Rider, J.R.; Tamimi, R.M.; Batista, J.L.; Haneuse, S.; Flynn-Evans, E.; Lockley, S.W.; et al. Circadian clock genes and risk of fatal prostate cancer. Cancer Causes Control. 2015, 26, 25–33. [Google Scholar] [CrossRef] [PubMed]
  23. Chu, L.W.; Zhu, Y.; Yu, K.; Zheng, T.; Yu, H.; Zhang, Y.; Sesterhenn, I.; Chokkalingam, A.P.; Danforth, K.N.; Shen, M.C.; et al. Variants in circadian genes and prostate cancer risk: A population-based study in China. Prostate Cancer Prostatic. Dis. 2008, 11, 342–348. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, Z.; Ma, F.; Zhou, F.; Chen, Y.; Wang, X.; Zhang, H.; Zhu, Y.; Bi, J.; Zhang, Y. Functional polymorphisms of circadian negative feedback regulation genes are associated with clinical outcome in hepatocellular carcinoma patients receiving radical resection. Med. Oncol. 2014, 31, 179. [Google Scholar] [CrossRef] [PubMed]
  25. Yuan, P.; Wang, S.; Zhou, F.; Wan, S.; Yang, Y.; Huang, X.; Zhang, Z.; Zhu, Y.; Zhang, H.; Xing, J. Functional polymorphisms in the NPAS2 gene are associated with overall survival in transcatheter arterial chemoembolization-treated hepatocellular carcinoma patients. Cancer Sci. 2014, 105, 825–832. [Google Scholar] [CrossRef] [PubMed]
  26. Couto, P.; Miranda, D.; Vieira, R.; Vilhena, A.; De Marco, L.; Bastos-Rodrigues, L. Association between CLOCK, PER3 and CCRN4L with non-small cell lung cancer in Brazilian patients. Mol. Med. Rep. 2014, 10, 435–440. [Google Scholar] [CrossRef]
  27. Zhu, Y.; Zheng, T. Clock-cancer connection in non-Hodgkin’s lymphoma. Med. Hypotheses 2008, 70, 788–792. [Google Scholar] [CrossRef]
  28. Benedetti, F.; Riccaboni, R.; Dallaspezia, S.; Locatelli, C.; Smeraldi, E.; Colombo, C. Effects of CLOCK gene variants and early stress on hopelessness and suicide in bipolar depression. Chronobiol. Int. 2015, 32, 1156–1161. [Google Scholar] [CrossRef]
  29. Kim, H.I.; Lee, H.J.; Cho, C.H.; Kang, S.G.; Yoon, H.K.; Park, Y.M.; Lee, S.H.; Moon, J.H.; Song, H.M.; Lee, E.; et al. Association of CLOCK, ARNTL, and NPAS2 gene polymorphisms and seasonal variations in mood and behavior. Chronobiol. Int. 2015, 32, 785–791. [Google Scholar] [CrossRef]
  30. Corella, D.; Asensio, E.M.; Coltell, O.; Sorlí, J.V.; Estruch, R.; Martínez-González, M.; Salas-Salvadó, J.; Castañer, O.; Arós, F.; Lapetra, J.; et al. CLOCK gene variation is associated with incidence of type-2 diabetes and cardiovascular diseases in type-2 diabetic subjects: Dietary modulation in the PREDIMED randomized trial. Cardiovasc. Diabetol. 2016, 15, 4. [Google Scholar] [CrossRef]
  31. Qu, F.; Qiao, Q.; Wang, N.; Ji, G.; Zhao, H.; He, L.; Wang, H.; Bao, G. Genetic polymorphisms in circadian negative feedback regulation genes predict overall survival and response to chemotherapy in gastric cancer patients. Sci. Rep. 2016, 6, 22424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Benna, C.; Rajendran, S.; Spiro, G.; Tropea, S.; Del Fiore, P.; Rossi, C.R.; Mocellin, S. Associations of clock genes polymorphisms with soft tissue sarcoma susceptibility and prognosis. J. Transl. Med. 2018, 16, 338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chen, Y.; Wang, D.; Song, Y.; Zhang, X.; Jiao, Z.; Dong, J.; Lu, L.; Zou, Z.; Du, W.; Qu, F. Functional polymorphisms in circadian positive feedback loop genes predict postsurgical prognosis of gastric cancer. Cancer Med. 2019, 8, 1919–1929. [Google Scholar] [CrossRef] [PubMed]
  34. Hoffman, A.E.; Yi, C.H.; Zheng, T.; Stevens, R.G.; Leaderer, D.; Zhang, Y.; Holford, T.R.; Hansen, J.; Paulson, J.; Zhu, Y. CLOCK in breast tumorigenesis: Genetic, epigenetic, and transcriptional profiling analyses. Cancer Res. 2010, 70, 1459–1468. [Google Scholar] [CrossRef] [Green Version]
  35. Hoffman, A.E.; Zheng, T.; Yi, C.H.; Stevens, R.G.; Ba, Y.; Zhang, Y.; Leaderer, D.; Holford, T.; Hansen, J.; Zhu, Y. The core circadian gene Cryptochrome 2 influences breast cancer risk, possibly by mediating hormone signaling. Cancer Prev. Res. 2010, 3, 539–548. [Google Scholar] [CrossRef] [Green Version]
  36. Fu, A.; Leaderer, D.; Zheng, T.; Hoffman, A.E.; Stevens, R.G.; Zhu, Y. Genetic and epigenetic associations of circadian gene TIMELESS and breast cancer risk. Mol. Carcinog. 2012, 51, 923–929. [Google Scholar] [CrossRef]
  37. Dai, H.; Zhang, L.; Cao, M.; Song, F.; Zheng, H.; Zhu, X.; Wei, Q.; Zhang, W.; Chen, K. The role of polymorphisms in circadian pathway genes in breast tumorigenesis. Breast Cancer Res. Treat. 2011, 127, 531–540. [Google Scholar] [CrossRef]
  38. Grundy, A.; Schuetz, J.M.; Lai, A.S.; Janoo-Gilani, R.; Leach, S.; Burstyn, I.; Richardson, H.; Brooks-Wilson, A.; Spinelli, J.J.; Aronson, K.J. Shift work, circadian gene variants and risk of breast cancer. Cancer Epidemiol. 2013, 37, 606–612. [Google Scholar] [CrossRef]
  39. Xiong, D.D.; Lv, J.; Wei, K.L.; Feng, Z.B.; Chen, J.T.; Liu, K.C.; Chen, G.; Luo, D.Z. A nine-miRNA signature as a potential diagnostic marker for breast carcinoma: An integrated study of 1,110 cases. Oncol. Rep. 2017, 37, 3297–3304. [Google Scholar] [CrossRef]
  40. Rajendran, B.; Janakarajan, V.N. Circadian clock gene aryl hydrocarbon receptor nuclear translocator-like polymorphisms are associated with seasonal affective disorder: An Indian family study. Indian J. Psychiatry 2016, 58, 57–60. [Google Scholar] [CrossRef] [Green Version]
  41. Xiang, S.; Coffelt, S.B.; Mao, L.; Yuan, L.; Cheng, Q.; Hill, S.M. Period-2: A tumor suppressor gene in breast cancer. J. Circadian Rhythms 2008, 6, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. 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, 41–50. [Google Scholar] [CrossRef] [Green Version]
  43. Lee, S.; Donehower, L.A.; Herron, A.J.; Moore, D.D.; Fu, L. Disrupting circadian homeostasis of sympathetic signaling promotes tumor development in mice. PLoS ONE 2010, 5, e10995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Takahashi, J.S.; Hong, H.K.; Ko, C.H.; McDearmon, E.L. The genetics of mammalian circadian order and disorder: Implications for physiology and disease. Nat. Rev. Genet. 2008, 9, 764–775. [Google Scholar] [CrossRef]
  45. Hoffman, A.E.; Zheng, T.; Ba, Y.; Stevens, R.G.; Yi, C.H.; Leaderer, D.; Zhu, Y. Phenotypic effects of the circadian gene Cryptochrome 2 on cancer-related pathways. BMC Cancer 2010, 10, 110. [Google Scholar] [CrossRef] [Green Version]
  46. Lavebratt, C.; Sjoholm, L.K.; Soronen, P.; Paunio, T.; Vawter, M.P.; Bunney, W.E.; Adolfsson, R.; Forsell, Y.; Wu, J.C.; Kelsoe, J.R.; et al. CRY2 is associated with depression. PLoS ONE 2010, 5, e9407. [Google Scholar] [CrossRef] [Green Version]
  47. Straif, K.; Baan, R.; Grosse, Y.; Secretan, B.; El Ghissassi, F.; Bouvard, V.; Altieri, A.; Benbrahim-Tallaa, L.; Cogliano, V. Carcinogenicity of shift-work, painting, and fire-fighting. Lancet Oncol. 2007, 8, 1065–1066. [Google Scholar] [CrossRef]
  48. Stevens, R.G.; Hansen, J.; Costa, G.; Haus, E.; Kauppinen, T.; Aronson, K.J.; Castano-Vinyals, G.; Davis, S.; Frings-Dresen, M.H.; Fritschi, L.; et al. Considerations of circadian impact for defining ‘shift work’ in cancer studies: IARC Working Group Report. Occup. Environ. Med. 2011, 68, 154–162. [Google Scholar] [CrossRef]
  49. Lesicka, M.; Jablonska, E.; Wieczorek, E.; Seroczynska, B.; Siekierzycka, A.; Skokowski, J.; Kalinowski, L.; Wasowicz, W.; Reszka, E. Altered circadian genes expression in breast cancer tissue according to the clinical characteristics. PLoS ONE 2018, 13, 15. [Google Scholar] [CrossRef] [Green Version]
  50. Lesicka, M.; Jablonska, E.; Wieczorek, E.; Seroczynska, B.; Kalinowski, L.; Skokowski, J.; Reszka, E. A different methylation profile of circadian genes promoter in breast cancer patients according to clinicopathological features. Chronobiol. Int. 2019, 36, 1103–1114. [Google Scholar] [CrossRef]
Figure 1. The bioinformatic analysis gene expression patterns in breast mammary tissue according to the selected genotypes of particular circadian genes using eQTL calculator on the GTex platform. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 1. The bioinformatic analysis gene expression patterns in breast mammary tissue according to the selected genotypes of particular circadian genes using eQTL calculator on the GTex platform. * p < 0.05; ** p < 0.01; *** p < 0.001.
Ijms 20 05704 g001
Table 1. Selected demographic and clinical characteristics of the breast cancer patients and healthy population in the association study.
Table 1. Selected demographic and clinical characteristics of the breast cancer patients and healthy population in the association study.
CharacteristicsBreast Cancer Patients
(n = 321)
Healthy Population
(n = 364)
p-Value
Age (years), mean (SD)58.85 (11.50)60.80 (7.11)0.007 1
BMI (kg/m2); mean (SD)27.13 (4.67)27.08 (4.62)0.90 1
Menopausal status
Pre-menopause77370.00001 2
Post-menopause233325
Unknown112
Smoking status
Never smokers2492250.00001 2
Past/Current smokers69138
Unknown31
Tumor stage
I157
II83
III58
IV9
Unknown14
Nodal status
Yes123
No183
Unknown15
Histological type
Carcinoma ductale174
Carcinoma lobulare30
Mix type of carcinoma115
Unknown2
Differentiation Grade (G)
Well27
Moderate151
Poor104
Unknown39
Estrogen receptor status (ER)
Positive227
Negative80
Unknown14
Progesterone Receptor status (PR)
Positive192
Negative115
Unknown14
HER2 receptor status
positive139
negative161
Unknown21
Notes: significant p-values are marked in bold. 1 The p values were calculated using the Student’s t-test. 2 The p values were calculated using the Pearson Chi-Square test.
Table 2. Selected circadian genes’ single nucleotide polymorphisms subject to analysis and their predicted functions.
Table 2. Selected circadian genes’ single nucleotide polymorphisms subject to analysis and their predicted functions.
GeneSNP IDChromosome PositionAllele
Major/Minor
RegionMAF GlobalMAF
Control Group
Probable FunctionGenotyping
BMAL1rs2279287chr11:13298485C/T5′UTR0.430.37TFBSHRM
CLOCKrs12505266chr4: 56412169C/Tintron0.330.43 1TFBSTaqMan probes
CLOCKrs1801260chr4:56301369A/G3′UTR0.250.33probable miRNA binding site hsa-miR-141HRM
CRY1rs8192440chr12:107395106C/Texon 50.210.36Splicing (ESE or ESS) Gly/GlyHRM
CRY2rs3824872chr11:45905605C/Anear gene0.460.24TFBSTaqMan probes
CRY2rs10838524chr11:45870177A/Gintron0.330.47N/ATaqMan probes
NPAS2rs2305160chr2:101591304G/Aintron0.20.36Splicing (ESE or ESS) Thr/AlaTaqMan probes
PER1rs2735611chr17:8048283G/Aexon 180.420.13Splicing (ESE or ESS) Gly/GlyTaqMan probes
PER1rs3027178chr17: 8053085T/Gexon 50.290.25 1Splicing (ESE or ESS) Thr/ThrTaqMan probes
PER2rs11894491chr2:239198325G/Aintron0.240.34TFBSTaqMan probes
PER2rs2304672chr2 239186589G/C5′UTR0.090.13Splicing (ESE or ESS)TaqMan probes
PER2rs934945chr2:239155053C/Texon 230.220.13CRY binding domain Gly/GluTaqMan probes
PER3rs10462020chr1:7880683T/Gexon 150.120.19Splicing (ESE or ESS) Val/GlyTaqMan probes
PER3rs2640909chr1:7830057T/Cexon 180.180.30 1Met/ThrTaqMan probes
TIMELESSrs2279665chr12:56827694G/Cexon 30.460.43Splicing (ESE or ESS) Leu/LeuTaqMan probes
TIMELESSrs774047chr12:56815922C/Texon 200.490.43Splicing (ESE or ESS) Gln/ArgTaqMan probes
1 Not in HWE; HRM High Resolution Melt; TFBS Transcription factor binding site; ESE Exonic splicing enhancer; ESS Exonic splicing silencer; N/A not applicable.
Table 3. Association between the selected single nucleotide polymorphisms (SNPs) of circadian genes with breast cancer risk.
Table 3. Association between the selected single nucleotide polymorphisms (SNPs) of circadian genes with breast cancer risk.
GeneGenotypeCtrlCasesOR (95%CI)p-Value 3
BMAL1
rs2279287
CC152159Ref.
CT1551270.76 (0.54–1.07)0.56
TT57280.47 (0.27–0.80)0.02
CT + TT 12121550.69 (0.50–0.95)0.02
TT 257280.53 (0.32–0.89)0.02
unknown07
CLOCK
rs1801260
AA160146Ref.
AG1681250.79 (0.56–1.12)0.18
GG34401.03 (0.60–1.78)0.57
AG + GG 12021650.84 (0.60–1.16)0.28
GG 234401.16 (0.69–1.94)0.58
unknown210
CLOCK
rs12505266
CC145134Ref.
CT125880.75 (0.51–1.10)0.03
TT94981.22 (0.83–1.79)0.057
TT + CT 12191860.95 (0.69–1.31)0.75
TT 294981.38 (0.97–1.96)0.07
unknown01
CRY1
rs8192440
CC148105Ref.
CT1651180.98 (0.68–1.41)0.94
TT49340.98 (0.58–1.66)0.97
CT + TT 12141520.98 (0.69–1.38)0.89
TT 249340.99 (0.60–1.62)0.97
unknown264
CRY2
rs10838524
AA10469Ref.
AG1771571.36 (0.91–2.01)0.73
GG83911.65 (1.05–2.58)0.07
AG + GG 12602481.45 (1.00–2.10)0.05
GG 283911.35 (0.93–1.95)0.11
unknown04
CRY2
rs3824872
CC210195Ref.
CA1331080.95 (0.68–1.34)0.73
AA21171.06 (0.53–2.13)0.82
CA + AA 11541250.97 (0.70–1.34)0.84
AA 221171.08 (0.54–2.14)0.83
unknown01
PER1
rs2735611
GG278222Ref.
GA79690.97 (0.66–1.44)0.01
AA7224.34 (1.77–10.67)0.001
GA + AA 186911.23 (0.85–1.77)0.26
AA 27224.37 (1.78–10.69)0.001
unknown08
PER1
rs3027178
TT214180Ref.
TG1201200.96 (0.68–1.35)0.24
GG30180.54 (0.28–1.08)0.09
TG + GG 11501380.88 (0.63–1.22)0.43
GG 230180.55 (0.28–1.08)0.08
unknown03
PER2
rs934945
CC275215Ref.
CT82911.52 (1.05–2.19)0.87
TT7112.08 (0.74–5.83)0.32
CT + TT 1891021.56 (1.09–2.23)0.01
TT 27111.85 (0.66–5.17)0.24
unknown04
PER2
rs2304672
GC280254Ref.
GC77600.86 (0.57–1.29)0.31
CC730.3 (0.06–1.46)0.16
GC + CC184630.81 (0.54–1.20)0.28
CC 2730.31 (0.06–1.50)0.14
unknown04
PER2
rs11894491
GG156139Ref.
GA1701360.94 (0.67–1.32)0.15
AA38451.48 (0.88–2.50)0.09
AA + GA 1 1.03 (0.75–1.42)0.86
AA 238451.53 (0.93–2.50)0.09
unknown01
PER3
rs10462020
TT241206Ref.
TG110951.01 (0.71–1.43)0.19
GG13191.98 (0.88–4.44)0.10
TG + GG 11231141.09 (0.78–1.53)0.60
GG 213191.97 (0.89–4.39)0.10
unknown01
PER3
rs2640909
TT193162Ref.
TC1251060.98 (0.69–1.40)0.33
CC45521.38 (0.85–2.23)0.16
TC + CC 11701581.08 (0.79–1.49)0.63
CC 245521.39 (0.88–2.20)0.16
unknown11
NPAS2
rs2305160
GG150131Ref.
GA1661360.92 (0.65–1.31)0.38
AA48491.16 (0.71–1.89)0.41
GA + AA 12141850.98 (0.71–1.35)0.89
AA 248491.21 (0.77–1.90)0.41
unknown05
TIMELESS
rs774047
CC11899Ref.
CT1831550.98 (0.68–1.41)0.35
TT63661.32 (0.83–2.10)0.18
CT + TT 12462211.06 (0.76–1.50)0.72
TT 263661.33 (0.88–2.01)0.17
unknown01
TIMELESS
rs2279665
GG114112Ref.
GC1881330.69 (0.48–1.00)0.02
CC62621.04 (0.66–1.66)0.29
GC + CC 12501950.78 (0.55–1.10)0.15
CC 262621.29 (0.86–1.95)0.22
unknown014
The significant values are marked in bold (p ≤ 0.05). OR, odds ratio; CI, confidence interval. 1 A dominant genetic model. 2 A recessive genetic model. 3 p-value adjusted by age, BMI, smoking status, menopausal status.
Table 4. Association between risk or protective SNPs of circadian genes with breast cancer predisposition.
Table 4. Association between risk or protective SNPs of circadian genes with breast cancer predisposition.
GeneCtrlCasesOR (95%CI)p-Value 3
One risk allele 1901001.53 (1.07–2.20)0.020
Two risk alleles 12128.21 (1.71–39.42)0.009
One or two risk alleles 1921121.67 (1.17–2.37)0.004
One or more risk 1 alleles931121.66 (1.17–2.35)0.005
One protective allele 277410.48 (0.30–0.77)0.002
Two protective alleles 2530.71 (0.16–3.08)0.650
One or two protective alleles 282440.49 (0.32–0.77)0.002
The significant values are marked in bold (p ≤ 0.05). OR, odds ratio; CI, confidence interval. 1 risk allele: CRY2 rs10838524; PER1 rs2735611; PER2 rs934945. 2 protective alleles: BMAL1 rs2279287; PER1 rs3027178. 3 p-value adjusted by age, BMI, smoking status, menopausal status.
Table 5. Selected SNPs of circadian genes and breast cancer risk according to the hormonal receptor status.
Table 5. Selected SNPs of circadian genes and breast cancer risk according to the hormonal receptor status.
ER+/PR+ER-/PR-
GeneGenotypesCtrlCasesOR (95%CI)p-Value 3CasesOR (95%CI)p-Value 3
BMAL1
rs2279287
CC15294Ref. 37Ref.
CT155760.74 (0.50–1.10)0.67300.83 (0.46–1.50)0.46
TT57150.45 (0.24–0.86)0.0460.42 (0.15–1.19)0.13
CT + TT 1212910.66 (0.46–0.97)0.03360.73 (0.42–1.28)0.27
TT 257150.52 (0.28–0.97)0.0460.46 (0.17–1.25)0.13
CLOCK
rs1801260
AA16088Ref. 33Ref.
AG168760.76 (0.51–1.13)0.52260.68 (0.37–1.26)0.07
GG34180.77 (0.40–1.5)0.70141.46 (0.64–3.35)0.16
AG + GG 1202940.76 (0.52–1.11)0.16400.83 (0.47–1.44)0.50
GG 234180.88 (0.47–1.66)0.70141.74 (0.80–3.82)0.17
CLOCK
rs12505266
CC14577Ref. 32Ref.
CT125520.81 (0.52–1.27)0.11240.80 (0.41–1.56)0.48
TT94591.28 (0.82–2.01)0.09201.00 (0.51–1.95)0.73
TT + CT 12191111.01 (0.69–1.48)0.95440.89 (0.51–1.55)0.68
TT 294591.40 (0.93–2.11)0.10201.09 (0.59–2.02)0.78
CRY1
rs8192440
CC14869Ref. 22Ref.
CT165610.83 (0.54–1.28)0.24301.13 (0.59–2.17)0.51
TT49241.16 (0.64–2.09)0.3960.82 (0.30–2.26)0.59
CT + TT 1214850.91 (0.61–1.36)0.63361.06 (0.57–1.97)0.86
TT 249241.27 (0.73–2.20)0.3960.76 (0.29–1.98)0.58
CRY2
rs10838524
AA10438Ref. 14Ref.
AG177961.52 (0.96–2.43)0.39371.32 (0.64–2.73)0.66
GG83511.67 (0.98–2.84)0.18242.26 (1.03–4.94)0.03
AG + GG 12601471.57 (1.01–2.44)0.05611.60 (0.81–3.16)0.17
GG 283511.26 (0.82–1.93)0.29241.87 (1.02–3.42)0.04
CRY2
rs3824872
CC210112Ref. 50Ref.
CA133681.07 (0.72–1.58)0.68240.66 (0.36–1.21)0.90
AA2180.91 (0.38–2.18)0.7820.50 (0.10–2.36)0.53
CA + AA 1154761.05 (0.72–1.53)0.81260.64 (0.36–1.15)0.14
AA 22180.89 (0.38–2.09)0.7920.57 (0.12–2.69)0.48
PER1
rs2735611
GG278131Ref. 52Ref.
GA79400.93 (0.58–1.47)0.03160.89 (0.45–1.77)0.02
AA7123.76 (1.38–10.24)0.0166.08 (1.80–20.56)0.003
GA+AA 186521.14 (0.75–1.75)0.54221.24 (0.67–2.28)0.50
AA 27123.82 (1.41–10.35)0.0166.23 (1.86–20.9)0.003
PER1
rs3027178
TT214113Ref. 45Ref.
TG120620.87 (0.58–1.30)0.79270.75 (0.41–1.38)0.57
GG30130.65 (0.31–1.40)0.3540.35 (0.10–1.31)0.17
TG + GG 1150750.82 (0.56–1.21)0.32310.67 (0.37–1.19)0.17
GG 230130.69 (0.33–1.45)0.3340.40 (0.11–1.43)0.16
PER2
rs934945
CC275122Ref. 52Ref.
CT82571.73 (1.13–2.63)0.59191.44 (0.76–2.71)0.64
TT762.02 (0.6–6.76)0.4873.16 (0.74–13.42)0.19
CT + TT 189631.75 (1.16–2.63)0.012021.57 (0.86–2.87)0.14
TT 2761.71 (0.52–5.68)0.3872.88 (0.69–12.08)0.15
PER2
rs2304672
GC280154Ref. 57Ref.
GC77300.71 (0.43–1.16)0.52171.01 (0.52–1.95)0.70
CC710.24 (0.03–1.98)0.2410.63 (0.07–5.37)0.67
GC + CC 184310.66 (0.41–1.07)0.10180.97 (0.51–1.85)0.92
CC 2710.25 (0.03–2.11)0.2010.63 (0.07–5.33)0.67
PER2
rs11894491
GG15682Ref. 30Ref.
GA170790.94 (0.63–1.4)0.15371.11 (0.62–1.98)0.95
AA38271.61 (0.88–2.94)0.0891.28 (0.49–3.30)0.66
AA + GA 12081061.05 (0.72–1.53)0.80461.14 (0.65–1.98)0.65
AA 238271.66 (0.94–2.93)0.0891.21 (0.49–2.97)0.68
PER3
rs10462020
TT241121Ref. 48Ref.
TG110580.99 (0.66–1.49)0.37231.15 (0.63–2.08)0.54
GG1381.69 (0.63–4.50)0.2952.2 (0.58–8.44)0.29
TG + GG 1123661.05 (0.71–1.55)0.82281.23 (0.70–2.18)0.47
GG 21381.69 (0.64–4.47)0.2952.10 (0.56–7.9)0.27
PER3
rs2640909
TT19397Ref. 40Ref.
TC125610.93 (0.62–1.41)0.51230.79 (0.42–1.5)0.11
CC45301.16 (0.66–2.03)0.50131.79 (0.83–3.86)0.06
TC + CC 1170910.99 (0.69–1.44)0.98361.03 (0.59–1.79)0.92
CC 245301.19 (0.7–2.04)0.52131.95 (0.94–4.07)0.07
NPAS2
rs2305160
GG15076Ref. 30Ref.
GA166810.97 (0.65–1.46)0.60331.05 (0.57–1.92)0.93
AA48281.17 (0.67–2.05)0.52111.04 (0.44–2.46)0.96
GA + AA 12141091.02 (0.7–1.49)0.93441.05 (0.59–1.85)0.87
AA 248281.19 (0.70–2.00)0.52111.02 (0.46–2.26)0.96
TIMELESS
rs774047
CC11859Ref. 22Ref.
CT183941.02 (0.67–1.55)0.81381.05 (0.56–1.99)0.65
TT63351.13 (0.66–1.96)0.64161.44 (0.65–3.17)0.34
CT + TT 12461291.05 (0.70–1.56)0.83541.15 (0.63–2.10)0.65
TT 263351.12 (0.69–1.83)0.64161.39 (0.70–2.76)0.34
TIMELESS
rs2279665
GG11471Ref. 20Ref.
GC188750.64 (0.42–0.98)0.08381.13 (0.58–2.20)0.51
CC62320.82 (0.48–1.42)0.91161.86 (0.83–4.17)0.11
GC + CC 12501070.69 (0.46–1.02)0.06541.30 (0.69–2.43)0.42
CC 262321.06 (0.65–1.74)0.81161.72 (0.87–3.40)0.12
The significant values are marked in bold (p ≤ 0.05). OR, odds ratio; CI, confidence interval. 1 A dominant genetic model. 2 A recessive genetic model. 3 p-value adjusted by age, BMI, smoking status, menopausal status.
Table 6. Gene expression of the selected circadian genes according to gene variants in the breast cancer tissues and adjacent non-tumor tissues.
Table 6. Gene expression of the selected circadian genes according to gene variants in the breast cancer tissues and adjacent non-tumor tissues.
GeneAdjacent Non-Tumor TissuesBreast Cancer Tissues
Mean ± SDN 1Mean ± SDN 1p-ValueMean ± SDN1Mean ± SDN1p-Value 2
BMAL1
rs2279287
CCCT + TT CCCT + TT
1.78 ± 0.37411.8 ± 0.29380.821.88 ± 0.34501.75 ± 0.30450.05
CLOCK
rs1801260
AA + AGGG AA + AGGG
1.82 ± 0.34721.97 ± 0.24110.172.03 ± 0.41822.00 ± 0.40150.79
CLOCK
rs12505266
CC + CTTT CC + CTTT
1.79 ± 0.37571.91 ± 0.27340.092.03 ± 0.42672.01 ± 0.4380.80
CRY1
rs8192440
CC + CTTT CC + CTTT
2.1 ± 0.38702.05 ± 0.2580.742.13 ± 0.38751.91 ± 0.56100.10
CRY2
rs10838524
CC + CAAA CC + CAAA
2.37 ± 0.39912.16 ± 0.6650.261.98 ± 0.45991.27 ± 0.4660.0004
CRY2
rs3824872
AG + GGAA AG + GGAA
2.38 ± 0.42752.27 ± 0.38170.341.97 ± 0.48831.81 ± 0.42190.18
PER1
rs3027178
TTTG + GG TTTG + GG
2.12 ± 0.45581.77 ± 0.37400.0021.63 ± 0.56641.66 ± 0.94410.88
PER1
rs2735611
GA + GGAA GA + GGAA
2.14 ± 0.45831.81 ± 0.4790.041.67 ± 0.56881.41 ± 0.67120.14
PER2
rs934945
CT + TTCC CT + TTCC
2.2 ± 0.45332.12 ± 0.31600.302.04 ± 0.4381.97 ± 0.46630.45
PER2
rs11894491
GA + GGAA GA + GGAA
2.11 ± 0.36862.42 ± 0.23110.0071.99 ± 0.44922.06 ± 0.49120.61
PER2
rs2304672
GGGC + CC CCCG + GG
2.16 ± 0.36752.1 ± 0.39180.611.99 ± 0.45822.05 ± 0.4190.57
PER3
rs10462020
TT TG + GG TT TG + GG
2.42 ± 0.34702.5 ± 0.38260.332.14 ± 0.52762.4 ± 0.55280.02
PER3
rs2640909
TC + TTCC TC + TTCC
2.44 ± 0.36652.46 ± 0.35310.782.22 ± 0.51712.18 ± 0.61330.76
NPAS2
rs2305160
GG + GAAA GG + GAAA
1.29 ± 0.54661.56 ± 0.66150.061.44 ± 0.64701.55 ± 0.65170.53
TIMELESS
rs2279665
GGGC + CC GGGC + CC
1.34 ± 0.48261.33 ± 0.57590.971.66 ± 0.64301.66 ± 0.54720.98
TIMELESS
rs774047
CC + CTTT CC + CTTT
1.31 ± 0.53731.44 ± 0.59150.381.69 ± 0.57801.58 ± 0.58250.39
The significant values are marked in bold (p ≤ 0.05); 1 N reflects the number of tissues with complete gene expression and genotyping results, 2 p-value obtained by means of the Student’s t-test.

Share and Cite

MDPI and ACS Style

Lesicka, M.; Jabłońska, E.; Wieczorek, E.; Pepłońska, B.; Gromadzińska, J.; Seroczyńska, B.; Kalinowski, L.; Skokowski, J.; Reszka, E. Circadian Gene Polymorphisms Associated with Breast Cancer Susceptibility. Int. J. Mol. Sci. 2019, 20, 5704. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20225704

AMA Style

Lesicka M, Jabłońska E, Wieczorek E, Pepłońska B, Gromadzińska J, Seroczyńska B, Kalinowski L, Skokowski J, Reszka E. Circadian Gene Polymorphisms Associated with Breast Cancer Susceptibility. International Journal of Molecular Sciences. 2019; 20(22):5704. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20225704

Chicago/Turabian Style

Lesicka, Monika, Ewa Jabłońska, Edyta Wieczorek, Beata Pepłońska, Jolanta Gromadzińska, Barbara Seroczyńska, Leszek Kalinowski, Jarosław Skokowski, and Edyta Reszka. 2019. "Circadian Gene Polymorphisms Associated with Breast Cancer Susceptibility" International Journal of Molecular Sciences 20, no. 22: 5704. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20225704

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop