Associations of clock genes polymorphisms with soft tissue sarcoma susceptibility and prognosis
Journal of Translational Medicinevolume 16, Article number: 338 (2018)
Dysfunction of the circadian clock and polymorphisms of some circadian genes have been linked to cancer development and progression. We investigated the relationship between circadian genes germline variation and susceptibility or prognosis of patients with soft tissue sarcoma.
Patients and methods
We considered the 14 single nucleotide polymorphisms (SNPs) of 6 core circadian genes that have a minor allele frequency > 5% and that are known to be associated with cancer risk or prognosis. Genotyping was performed by q-PCR. Peripheral blood and clinic-pathological data were available for 162 patients with liposarcoma or leiomyosarcoma and 610 healthy donors. Associations between the selected clock genes polymorphisms and sarcoma susceptibility or prognosis were tested assuming 3 models of inheritance: additive, recessive and dominant. Subgroup analysis based on sarcoma histotype was performed under the additive genetic model. Multivariate logistic regression and multivariate Cox proportional hazard regression analyses were utilized to assess the association between SNPs with patient susceptibility and survival, respectively. Pathway variation analysis was conducted employing the Adaptive Rank Truncated Product method.
Six out of the 14 analyzed SNPs were statistically significantly associated with susceptibility or prognosis of soft tissue sarcoma (P < 0.05). The present analysis suggested that carriers of the minor allele of the CLOCK polymorphism rs1801260 (C) or of PER2 rs934945 (T) had a reduced predisposition to sarcoma (26% and 35% respectively with the additive model) and liposarcoma (33% and 41% respectively). The minor allele (A) of NPAS2 rs895520 was associated with an increased predisposition to sarcoma of 33% and leiomyosarcoma of 44%. RORA rs339972 C allele was associated with a decreased predisposition to develop sarcoma assuming an additive model (29%) and leiomyosarcoma (36%). PER1 rs3027178 was associated with a reduced predisposition only in liposarcoma subgroup (32%). rs7602358 located upstream PER2 was significantly associated with liposarcoma survival (HR: 1.98; 95% CI 1.02–3.85; P = 0.04). Germline genetic variation in the circadian pathway was associated with the risk of developing soft tissue sarcoma (P = 0.035).
Genetic variation of circadian genes appears to play a role in the determinism of patient susceptibility and prognosis. These findings prompt further studies to fully dissect the molecular mechanisms.
Circadian rhythms are approximately 24 h oscillations of biochemical, physiologic and behavioral processes, present in almost all living organisms, which arise from an ancient adaptation to the rotation of the earth. Disruption of circadian rhythms has been related to different diseases as diabetes, depression, sleep disorders, obesity, heart attack and cancer . Since the eighties it has been hypothesized that the high risk of human breast and prostate cancer in industrialized societies was caused by circadian rhythms deregulation . In 2007 an expert panel assembled by the International Agency for Research on Cancer (IARC) concluded that “shift-work that involves circadian disruption is probably carcinogenic to humans” . The biological clock is a self-sustained mechanism able to maintain and synchronize circadian rhythms via transcription-translation feedback loops, constituted by core circadian clock genes. They can be divided in: positive activators as CLOCK (clock circadian regulator), NPAS2 (neuronal PAS domain protein 2), ARNTL (aryl hydrocarbon receptor nuclear translocator like, also referred to as brain and muscle Arnt-like protein-1, BMAL1), RORA (RAR related orphan receptor A), and NR1D1 (nuclear receptor subfamily 1 group D member 1 also known as Rev-Erb alpha); negative effectors as CRY1 (cryptochrome circadian clock 1), CRY2 (cryptochrome circadian clock 2), PER1 (period circadian clock 1), PER2 (period circadian clock 2), PER3 (period circadian clock 3) and NR1D2 (nuclear receptor subfamily 1 group D member 2 also known as Rev-Erb beta) and modulators as CSNK1E (casein kinase I epsilon). An additional clock-related gene which probably functions as modulator is TIMELESS (timeless circadian clock).
The circadian system is composed by the central clock, located in the suprachiasmatic nucleus of the brain and by peripheral clocks, located in virtually all body tissues. The two components communicate and synchronize with each other and, in particular, the central clock controls the peripheral clocks . Moreover, the central clock modulates the expression of the so-called clock-controlled genes, which are estimated to be approximately 20% of the genes in mammals [5, 6], many of them regulating cancer-related biological pathways such as cell proliferation, apoptosis, DNA damage and repair, carcinogen metabolism and/or detoxification [7, 8].
A different approach was undertaken for the first time in 2005. Zhu and Colleagues studied a structural variant in the circadian gene PER3, which was detected to be significantly associated with increased risk of breast cancer in young women . Subsequently, in a growing number of molecular epidemiological studies, germline variations in clock genes have been associated by several Authors with different type of cancer susceptibility  and in some cases with the prognosis of cancer patients [11,12,13,14,15,16,17]. Recently the discovery rate of susceptibility loci is being greatly accelerated by genome-wide association studies (GWASs) which can test up to one million single nucleotide polymorphisms (SNPs) in thousands of subjects at a time , nevertheless associations between clock genes variations and sarcoma susceptibility has not been explored yet.
Sarcomas are a family of rare malignant tumors arising from bone and soft tissues with more than 50 different histologies accounting for about 1–2% of cancers in adults and 15–20% in children (worldwide incidence: approximately 200,000 cases per year) .
Given the potential role of circadian genes in tumorigenesis, it has been hypothesized that genetic variations in these genes could be associated with an individual’s susceptibility also with sarcoma. In this exploratory analysis, we examined whether 14 common genetic variants in 6 circadian pathway genes are related to soft tissue sarcoma susceptibility or patient outcome in a series of 162 diagnosed soft tissue sarcoma cases and 610 healthy controls.
Materials and methods
We conducted a retrospective study to test the hypothesis that genetic variation (in terms of candidate single nucleotide polymorphisms, SNPs) of the circadian pathway might be associated to the susceptibility and the prognosis of patients affected with sarcoma. To this aim we extracted the clinico-pathological data of patients treated at our institution (University Hospital of Padova, North East of Italy) between 1992 and 2016, using a prospectively maintained database linked to our institutional biobank (Clinica Chirurgica I–Istituto Oncologico Veneto). To be included in the study, each case had to meet the following requirements: (1) histologically confirmed diagnosis of liposarcoma or leiomyosarcoma; (2) pathology-based information on TNM stage; (3) follow-up data (minimum follow up: 6 months); (4) availability of peripheral blood for genotyping purposes.
Patients and healthy donors
We selected retrospectively 162 consecutive sarcoma patients, 93 patients bearing liposarcoma and 69 patients bearing leiomyosarcoma, and 610 healthy controls. Healthy controls selection was both population-based (n = 270 blood donors) and hospital-based (n = 340, healthy subjects who visited the Clinica Chirurgica I ambulatories for routine check-ups). All patients signed an informed consent form explaining the research purposes of the blood withdrawal.
We focused on 5 of the 12 core clock genes, which are CLOCK, NPAS2, PER1, PER2, RORA. Moreover, we added a clock-related gene TIMELESS, associated with cancer risk in several studies [20,21,22,23]. The CLOCK locus contains 28734 SNPs annotated by the National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/) SNP data repository. The Genome Variation Server of the University of Washington (http://gvs.gs.washington.edu/GVS/) and the TagSNP tool of the US National Institute of Environmental Health Sciences (https://snpinfo.niehs.nih.gov/snpinfo/snptag.html) were interrogated to identify the SNPs that tag CLOCK SNPs with a minimum correlation coefficient (r2) of at least 0.80, minimal genotype data coverage of 50%, and minimal allele frequency of 5%, in the Caucasian population (CEU) genotyped by the HapMap Project. This process yielded an initial group of 134 Tag SNPs for genotyping, divided in 11 bins (groups of SNPs in high linkage disequilibrium, D’ greater than 0.8). Excluding those SNPs tagging less than 10 SNPs, we selected from each bin the Tag SNPs with more functional/cancer related information available. They are the following: rs1801260, rs3736544, rs3749474. Moreover, we included rs34897046 which has a missense functional effect.
The TIMELESS locus contains 9200 SNPs annotated by the NCBI dbSNP. With the same procedure described, rs774027 was chosen as TagSNP. In order to enrich our analysis, rs3809125 (3′-UTR) and rs7302060 were added based on literature (see Table 2) .
For NPAS2, PER1, PER2, RORA we decided to follow a different approach. Despite the elevated number of Tag SNPs for each gene (more than 200 in NPAS2 locus) none of them tags more than 9 SNPs and most of them tags 1 SNP. We relied on our previous meta-analysis  or on literature (see Table 2).
NPAS2 rs895520, PER2 rs7602358, RORA rs339972, and rs10519097 had a statistically significant association with cancer risk. Moreover, regarding the NPAS2 locus, we selected the most studied variant, based on number of datasets, with missense functional effect, rs2305160. PER1 rs3027178 and PER2 rs934945 were selected based on literature [24,25,26,27,28].
DNA extraction and genotyping
Genomic DNA was isolated from peripheral whole blood employing the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions and quantified by Nanodrop 1000 spectrophotometer (Thermofisher Scientific, Waltham, MAUSA).
10 to 20 ng of DNA of each patient were used for TaqMan SNP genotyping assays (comprising primers and fluorescent probes) according to the manufacturer’s instructions (Thermofisher Scientific, Waltham, MA USA). Genotyping was performed by real–time PCR using either allelic discrimination in the 7300 RT–PCR System (Thermofisher Scientific, Waltham, MA USA), either using endpoint genotyping in an LightCycler 480II (Roche Molecular Diagnostics Pleasanton, CA, USA). For the 7300 RT–PCR System PCR parameters involved an initial denaturation at 95 °C for 10 min followed by 40 cycles at 95 °C for 15 s and 60 °C for 1 min. Post run data were analysed by 7300 SDS software (Thermofisher Scientific, Waltham, MA USA) and Automatic calls were assigned with approximately 99.8% quality. A call rate > 95% was considered the cutoff to consider genotyping. For LightCycler 480II PCR parameters involved an initial denaturation at 95 °C for 10 min followed by 40 cycles at 95 °C for 10 s, 60 °C for 1 min, and 72 °C for 10 s. Post run data were analysed by LightCycler 480 Endpoint Genotyping software (Roche Molecular Diagnostics Pleasanton, CA, USA). Two blank (water) controls in each 96-well plate were used for the assay quality control.
We considered three genetic models:
Additive (or allelic; presence of 0 vs 1 vs 2 copies of the variant allele): the genotype for each SNP is considered as a continuous variable;
Recessive (presence of 2 copies of the variant allele vs 0 or 1 copies of the variant allele): the genotype for each SNP is considered as a categorical variable with 2 categories (homozygous for the common allele + heterozygous, homozygous for the variant allele);
Dominant (presence of at least one copy of the variant allele vs no copies of the variant allele): the genotype for each SNP is considered as a categorical variable with 2 categories (homozygous for the variant allele + heterozygous, homozygous for the common allele).
For susceptibility of sarcoma assessment multivariate logistic regression was employed for each genetic model. Odds ratios (OR) and 95% confidence intervals were used as a measure of association. In those multivariate models the evaluated outcome was the presence or absence of sarcoma, while the explanatory variables were the single SNPs adjusted for age and gender.
For prognosis of sarcoma assessment multivariate Cox proportional hazard regression was employed. Overall survival was defined as the time from the date of tumor diagnosis to the date of death by any cause or last follow-up visit. Hazard ratios (HR) and 95% confidence intervals were used as a measure of association. In those multivariate models the evaluated event was the patient’s survival, the time to event were the months of survival, and the explanatory variables were the single SNP adjusted for age, gender and sarcoma stage.
For discriminating whether sarcoma subtype interacts with SNP associations with susceptibility or prognosis, the above mentioned statistical analyses were repeated, for the liposarcoma subgroup and for the leiomyosarcoma subgroup, employing the additive model because the additive model can be considered a conservative choice between recessive and dominant models.
Hardy–Weinberg equilibrium (HWE) was tested for both samples (patients and healthy controls) for each SNP employing OEGE—Online Encyclopedia for Genetic Epidemiology studies , http://www.oege.org/software/hwe-mr-calc.shtml. This tool is a HWE calculator for biallelic SNPs based on Chi-square statistic.
Statistical power was calculated for each SNP employing the on-line tool “Power and Sample Size” of the University of Vanderbilt (http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize) . Power was defined as the probability of correctly rejecting the null hypothesis that the relative risk (OR) was equal to 1, given 162 case patients and 610 controls. The type I error probability α was set to 0.05 and ψ (OR considered clinically relevant in our study) was set to 0.80.
For all the analysis conducted in this study, a P value < 0.05 was used as the cut-off for significance, not adjusted for multiple comparisons.
Rcmdr: R Commander. R package version 2.4-2 was employed for the analyses.
Pathway variation analysis
We also investigated the associations of soft tissue sarcoma with genes, observed as combinations of SNPs, and with the circadian rhythm pathway as a whole, observed as combinations of 6 clock genes. This analysis was conducted employing the Adaptive Rank Truncated Product (ARTP) method  as we already described in detail in .
The analysis was based on a total of 162 cases of sarcoma and 610 healthy controls, all of European ancestry. Subjects’ characteristics and the main features of the SNPs we investigated are summarized in Tables 1 and 2, respectively.
A total of 14 preselected SNPs in 6 circadian clock genes were successfully genotyped, and no departures from Hardy–Weinberg equilibrium were observed neither among the controls nor among the patients, as reported in Additional file 1: Table S1 (P > 0.05). The CLOCK rs34897046 SNP showed minor allele frequencies < 5% across the whole study population (3% in controls and 4% in cases, see Additional file 1: Table S1) and was not evaluated further statistically.
Statistical power for each SNP for primary analysis (sarcoma) and subgroup analysis (liposarcoma and leiomyosarcoma) is reported in Additional file 2: Table S2.
Associations between the selected clock genes polymorphisms and sarcoma predisposition were tested assuming 3 models of inheritance: additive, recessive and dominant. We used odds ratios (ORs) and their corresponding 95% confidence intervals (95% CI) to measure the strength of association between each polymorphism and sarcoma risk. The results are reported in Table 3. We examined further subgroup SNP associations according to sarcoma histology: liposarcoma (n = 93) and leiomyosarcoma (n = 69). The results are reported in Table 4.
CLOCK maps on chromosome 4 at 4q12. CLOCK putative protein belongs to the basic helix-loop-helix PAS family of transcription factors and forms heterodimers with ARNTL (BMAL1) to enhance target gene expression. CLOCK is also involved in growth arrest, DNA repair and apoptosis upon genotoxic stress caused by UV radiation, suggesting that this molecule may represent an important “caretaker” promoting cell cycle arrest upon DNA damage . CLOCK has the properties of a histone acetyl transferase and is involved in chromatin remodeling .
rs1801260 is located on 3′-UTR region of CLOCK. Carriers of the minor allele (C) rs1801260 had a reduced predisposition to sarcoma under an additive (per allele OR 0.74; 95% CI 0.55–1.00; P = 0.05) and a recessive (OR 0.40; 95% CI 0.18–0.88; P = 0.02) genetic model. Subgroup analysis suggested that the protective effect of the C allele was significantly associated with liposarcoma under an additive model (per allele OR 0.67; 95% CI 0.46–0.98; P = 0.04). No significant association was found with leiomyosarcoma.
NPAS2 is the largest human clock gene. It maps on chromosome 2 at 2q11.2 and, like its paralogue CLOCK, encodes for a member of the basic helix-loop-helix PAS class of transcription factors . When dimerized with ARNTL, NPAS2 binds to E-box regulatory elements in target promoter regions and enhances target gene expression. Previous studies reported NPAS2 as a putative tumor suppressor . rs895520 is located on an intron of NPAS2 locus. The present analysis suggested that the minor allele (A) was associated with an increased predisposition to sarcoma of 33% under an additive (per allele OR 1.33; 95% CI 1.02–1.73; P = 0.03) and a recessive (OR 1.70; 95% CI 1.08–2.68; P = 0.02) model. In subgroup analysis this association was statistically significant in leiomyosarcoma (per allele OR 1.44; 95% CI 0.99–2.09; P = 0.05).
PERs code for PAS domain-containing key regulators of the circadian clock. PER genes control their own transcription by directly repressing ARNTL heterodimers, their activators . Moreover, it has been suggested that PER1 and PER2 function as tumor suppressors [51,52,53,54,55,56]. PER2 gene can activate c-Myc signaling pathways leading to genomic instability and cell proliferation. PER2 dysfunction can also impair p53-mediated apoptosis and consequently result in genomic instability and the accumulation of damaged cells [55, 57]. PER1 maps on chromosome 17 at 17p13.1 and PER2 on chromosome 2 at 2q37.3.
rs934945 (C>T) is located on the last exon of PER2 locus and has a missense functional effect, leading to the substitution Glycine–Glutamic Acid. Carriers of the minor allele had a decreased predisposition to develop sarcoma (35%) employing an additive genetic model (per allele OR 0.65; 95% CI 0.45–0.94; P = 0.02) or a dominant model and (OR 0.63; 95% CI 0.42–0.95; P = 0.03). This association is statistically significant also in liposarcoma subgroup (per allele OR 0.59; 95% CI 0.37–0.96; P = 0.03). Moreover, in liposarcoma subgroup PER1 rs3027178, a genetic variant with a synonymous functional effect, was associated with a reduced predisposition (per allele OR 0.68; 95% CI 0.47–0.98; P = 0.04).
RORA maps on chromosome 15 at 15q21-q22, spans a 730 kb large genomic region comprised of 15 exons and encodes for one member of the retinoid orphan nuclear receptor subfamily of orphan receptors. RORA has been reported as potential tumor suppressor [58, 59]. ARNTL-CLOCK or ARNTL-NPAS2 heterodimers promote the transcription of RORA, which in turn activates the transcription of ARNTL.
rs339972 C allele was associated with a decreased predisposition to develop sarcoma assuming an additive (per allele OR 0.71;95% CI 0.53–0.96; P = 0.02) or a dominant (OR 0.64; 95% CI 0.44–0.93; P = 0.02) genetic model. This association is statistically significant in leiomyosarcoma subgroup (per allele OR 0.64; 95% CI 0.42–0.98; P = 0.04).
No significant associations were found for the selected TIMELESS SNPs and the predisposition to develop sarcoma.
rs7602358 located upstream PER2 was significantly associated with liposarcoma survival (HR 1.98; 95% CI 1.02–3.85; P = 0.04) employing an additive genetic model. Carriers of the minor allele (G) had an increased risk of mortality of 98%. No further statistically significant associations with prognosis were found, neither in primary analysis nor in subgroup analysis. Results are reported in Additional file 3: Table S3.
Pathway variation analysis
The results of the pathway variation analysis are reported in Table 5. We found a significant association between circadian pathway variation and risk of developing this tumour (circadian pathway P value 0.035). This result was based on the data regarding 12 SNPs located in six genes. The top genes were PER2 (2 SNPs, circadian gene P value 0.036) and RORA (2 SNPs, circadian gene P value 0.050).
In the present article we described the findings of the first study, to our knowledge, investigating the relations between circadian clock genes DNA genetic variations and the susceptibility to soft tissue sarcoma or the outcome of sarcoma bearing patients. The results were based on genotyping data from 772 people enrolled retrospectively with a control-to-case ratio equal to 3.
We hypothesized that circadian clock genes SNPs may influence the susceptibility to soft tissue sarcoma as suggested for many different types of cancer as breast, prostate, colorectal, ovarian, pancreatic, lung, glioma, and non-Hodgkin lymphoma . Our results argue in favor of this hypothesis, in fact 4 of the 13 analyzed genetic variant were statistically significantly associated with the predisposition to sarcoma employing an additive model of inheritance. They are: CLOCK rs1801260, NPAS2 rs895520, PER2 rs934945, and RORA rs339972. In our previous meta-analysis on sarcoma genetic variations  we highlighted the interesting fact that the susceptibility, defined as odds ratio, statistically significantly associated with single variants ranged between higher values (mean approximately 1.70) than those usually observed for other malignancies such as breast , colorectal , and gastric carcinomas , which generally include odds ratios between 1.10 and 1.30. The present study is in line with this observation (mean approximately 1.40) fueling the speculation that germline DNA variation is especially important in the determinism of the predisposition to this family of tumors.
Susceptibility to soft tissue sarcoma
CLOCK rs1801260, located on 3′-UTR region, was studied by many Authors in relation to breast cancer [27, 33, 34], colorectal cancer [15, 25], esophageal carcinoma [12, 13], and gastric cancer . Our previous meta-analysis on genetic variation of clock genes and cancer risk  failed to reveal an association with cancer risk. Primary meta-analysis relied upon 4 datasets, including 2578 cases and 3349 controls, 3 out of 4 were based on breast cancer patients. Stratified analysis by cancer type showed lack of association with breast cancer. The fourth dataset was retrieved by the publication of Karantanos and Colleagues  on colorectal cancer. Karantanos found that the minor allele (C) and the CC genotype of rs1801260 polymorphism conferred an increased risk (78%) for the CRC development. Our findings in the present study showed that carriers of the C allele of rs1801260 had a reduced predisposition to sarcoma (26%).
CLOCK rs1801260 was considered by Authors also in relation to outcome. According to a study by Zhou et al.  who evaluated the association between clock genes polymorphisms and prognosis in patients with colorectal cancer, the presence of T allele and the TT genotype of rs1801260 SNP were related to decreased overall survival and unfavorable prognosis. In gastric cancer our group  suggested that rs1801260 C allele affect patient survival only if combined with the major allele of CLOCK rs3749474. In the present study we found no association with patients’ prognosis.
Considering all these observations one may speculate that the C allele effect is strongly dependent on the tumor site.
NPAS2 rs895520 is located on an intron of NPAS2 locus. Our previous meta-analysis of four datasets (including 19,865 subjects) revealed a highly significant association with an intermediate level of evidence with cancer risk (summary OR 1.08; 95% CI 1.03–1.13; P = 0.001). This result is also congruent with findings from a NCI genome-wide association study on prostate cancer (CGEMS project), released in 2006 and not included in our meta-analysis, that showed some significant associations between variants in circadian genes and prostate cancer susceptibility. The CGEMS project genotyped 550,000 SNPs in 1182 prostate cancer cases and 1174 controls . A total of 104 SNPs in circadian genes were included. Eight of these, located in four genes, NPAS2 (2 SNPs), CSNK1E (3 SNPs), CRY1 (2 SNPs), and CRY2 (1 SNP), were significantly associated with prostate cancer risk. NPAS2 rs895521 was associated with P ≤ 0.01. The present analysis suggested that the minor allele (A) was associated with an increased predisposition to sarcoma of 33%.
PER2 rs934945 has a missense functional effect causing glycine to glutamic acid substitution in PER2 protein and is probably related to a decreased PER2 activity. Dai and Colleagues  suggested that individuals carrying both the CLOCK rs3805151 (not considered in our study) CC and PER2 TT genotype had an increased breast cancer risk (OR 2.28; 95% CI 1.22–4.26). In the current study we evaluated the association between this polymorphism and the predisposition to develop sarcoma and we found that carriers of the minor allele (T) had a decreased risk (35%). As for CLOCK the effect could be attributed to a different susceptibility depending on the cancer type.
RORA rs339972 is located on an intron of RORA locus. In our previous meta-analysis this SNP was found to be associated to the risk of cancer (summary OR 1.08; 95% CI 1.01–1.15; P = 0.02). Primary meta-analysis relied on 2 datasets of breast cancer patients and one of pancreatic cancer. In the current study rs339972 C allele was associated with a decreased predisposition to develop sarcoma (per allele OR 0.71; 95% CI 0.53–0.96; P = 0.02).
Liposarcoma and leiomyosarcoma subgroups
Stratification by sarcoma histotype showed also potentially interesting results. Employing an additive model of inheritance CLOCK rs1801260 and PER2 rs934945 were statistically significantly associated with liposarcoma, while NPAS2 rs895520 and RORA rs339972 were statistically significantly associated with leiomyosarcoma. Moreover, PER1 rs3027178 was found to be associated only with liposarcoma susceptibility. It can be argued that different circadian clock genes influence specifically a particular sarcoma subtype. NPAS2 is a paralog of CLOCK. Both proteins are major regulators of the molecular clock, and act by forming heterodimers with ARNTL and promoting transcription of target genes. However, while CLOCK is mainly expressed in the ‘‘central pacemaker’’ of the circadian system, the suprachiasmatic nucleus of the hypothalamus, NPAS2 is expressed mainly in the forebrain . This suggests that while these 2 genes are functionally analogous, they might be involved in separate circadian-controlled processes. The circadian clock directs nearly all aspects of diurnal physiology, including lipid metabolism and fat cell differentiation [66,67,68]. Recently, CLOCK and PER2 polymorphisms have been linked to obesity and the metabolic syndrome [69, 70]. Since it has been proposed that liposarcoma could arise from the dedifferentiation of fat cells , one could hypothesize a specific role of CLOCK and PER2 genes and this aspect, worth to be investigated, could open the avenue to new therapeutic targets.
Genetic variation in the circadian system has been associated with tumor aggressiveness or patient survival in hepatocellular carcinoma [28, 46], melanoma  as well as prostate , colorectal , gastric [16, 17] and breast  cancer. The present results also suggest a potential impact of circadian clock on liposarcoma survival. In particular, rs7602358 located upstream PER2, was significantly associated with liposarcoma prognosis. rs7602358 was previously considered by four research groups [23, 27, 38, 40] evaluating the associations with prostate cancer, breast cancer, and glioma risk. Zhu and Colleagues  found a difference in risk association of rs7602358 with prostate cancer between less aggressive and more aggressive prostate cancer subgroup, leading to suppose a role of PER2 in malignant cells aggressiveness.
Pathway variation analysis
We used a gene- and pathway-based approach to investigate the overall effect of circadian clock gene germline variations on soft tissue sarcoma risk. This approach is useful to detect the combined effects of genetic polymorphisms that are weakly associated with the disease but may not be detected in single-SNP analyses and may provide additional insights into the mechanisms underlying disease susceptibility. Our results suggest that genetic variation in the circadian rhythm pathway as a whole is related to soft tissue sarcoma susceptibility. This association was mostly driven by PER2 and RORA. To the best of our knowledge, this is the first time that ARTP-based gene and pathway analysis has been applied to the relationship between circadian genes’ germline variation and soft tissue sarcoma susceptibility.
Strength and limitations
Our study is notable for several strengths. To the best of our knowledge, this is the first study investigating the associations of germline circadian clock genes polymorphisms in relation to risk or prognosis of soft tissue sarcoma. Given the rarity of sarcoma, we leveraged the existing resources of our biobank to evaluate a relatively large cohort of patients with available germline DNA. We used several strategies to select the potentially interesting variations: by using tag SNPs, we were able to efficiently interrogate multiple regions of the CLOCK and TIMELESS genes.
Our study is also limited by a number of weaknesses. Due to sample size considerations and to our resources, we limited the number of circadian clock genes selected and SNPs evaluated. We considered three genetic models of inheritance without correcting for multiple testing. Any multiple testing correction would probably nullify statistical significance, nevertheless the strength of those association is noticeably higher than that usually observed for other tumor types. However, we considered this work a starting point and we did not know a priori the best fitted model. The power of most of our comparisons is around 30–40% and a larger sample would be needed to reach the commonly used statistical power of 80%. We pooled together different sarcoma histotypes to assess whether there are associations linked to soft tissue sarcoma in general but, given the complexity of this pathology, this could be considered a pitfall. In this regard we performed subgroup analysis and in fact we found that potentially histotype interacts with different SNPs. The source of controls is mixed, population and hospital based, nevertheless to avoid selection biases for the hospital based fraction we chose patients from different conditions.
In conclusion, the present work represents the first article on association between genetic variation of clock-related genes and soft tissue sarcoma susceptibility or prognosis. Our findings support this relationship and might become a useful starting point for future investigation, which is certainly needed to shed more light in this promising field of cancer predisposition.
Roenneberg T, Merrow M. The circadian clock and human health. CurrBiol. 2016;26(10):R432–43.
Stevens RG. Electric power use and breast cancer: a hypothesis. Am J Epidemiol. 1987;125(4):556–61.
Straif K, Baan R, Grosse Y, Secretan B, El Ghissassi F, Bouvard V, et al. Carcinogenicity of shift-work, painting, and fire-fighting. Lancet Oncol. 2007;8(12):1065–6.
Yamazaki S, Numano R, Abe M, Hida A, Takahashi R, Ueda M, et al. Resetting central and peripheral circadian oscillators in transgenic rats. Science. 2000;288(5466):682–5.
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.
Koike N, Yoo SH, Huang HC, Kumar V, Lee C, Kim TK, et al. Transcriptional architecture and chromatin landscape of the core circadian clock in mammals. Science. 2012;338(6105):349–54.
Fu L, Lee CC. The circadian clock: pacemaker and tumour suppressor. Nat Rev Cancer. 2003;3(5):350–61.
Sahar S, Sassone-Corsi P. Metabolism and cancer: the circadian clock connection. Nat Rev Cancer. 2009;9(12):886–96.
Zhu Y, Brown HN, Zhang Y, Stevens RG, Zheng T. Period3 structural variation: a circadian biomarker associated with breast cancer in young women. Cancer Epidemiol Biomarkers Prev. 2005;14(1):268–70.
Benna C, Helfrich-Forster C, Rajendran S, Monticelli H, Pilati P, Nitti D, et al. Genetic variation of clock genes and cancer risk: a field synopsis and meta-analysis. Oncotarget. 2017;8(14):23978–95.
Lin DW, FitzGerald LM, Fu R, Kwon EM, Zheng SL, Kolb S, et al. Genetic variants in the LEPR, CRY1, RNASEL, IL4, and ARVCF genes are prognostic markers of prostate cancer-specific mortality. Cancer Epidemiol Biomarkers Prev. 2011;20(9):1928–36.
Miki I, Tamura T, Nakamura T, Makimoto H, Hamana N, Uchiyama H, et al. Circadian variability of pharmacokinetics of 5-fluorouracil and CLOCK T3111C genetic polymorphism in patients with esophageal carcinoma. Ther Drug Monit. 2005;27(3):369–74.
Okuno T, Tamura T, Yamamori M, Chayahara N, Yamada T, Miki I, et al. Favorable genetic polymorphisms predictive of clinical outcome of chemoradiotherapy for stage II/III esophageal squamous cell carcinoma in Japanese. Am J ClinOncol. 2007;30(3):252–7.
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.
Zhou F, He X, Liu H, Zhu Y, Jin T, Chen C, et al. Functional polymorphisms of circadian positive feedback regulation genes and clinical outcome of Chinese patients with resected colorectal cancer. Cancer. 2012;118(4):937–46.
Rajendran S, Benna C, Monticelli H, Spiro G, Menin C, Mocellin S. Germline variation of circadian pathway genes and prognosis of gastric cancer patients. Gut. 2018;67(4):779–80.
Qu F, Qiao Q, Wang N, Ji G, Zhao H, He L, et al. Genetic polymorphisms in circadian negative feedback regulation genes predict overall survival and response to chemotherapy in gastric cancer patients. Sci Rep. 2016;1(6):22424.
Visscher PM, Wray NR, Zhang Q, Sklar P, McCarthy MI, Brown MA, et al. 10 years of GWAS discovery: biology, function, and translation. Am J Hum Genet. 2017;101(1):5–22.
Taylor BS, Barretina J, Maki RG, Antonescu CR, Singer S, Ladanyi M. Advances in sarcoma genomics and new therapeutic targets. Nat Rev Cancer. 2011;11(8):541–57.
Fu A, Leaderer D, Zheng T, Hoffman AE, Stevens RG, Zhu Y. Genetic and epigenetic associations of circadian gene TIMELESS and breast cancer risk. MolCarcinog. 2012;51(12):923–9.
Grundy A, Schuetz JM, Lai AS, Janoo-Gilani R, Leach S, Burstyn I, et al. Shift work, circadian gene variants and risk of breast cancer. Cancer Epidemiol. 2013;37(5):606–12.
Jim HS, Lin HY, Tyrer JP, Lawrenson K, Dennis J, Chornokur G, et al. Common genetic variation in circadian rhythm genes and risk of epithelial ovarian cancer (EOC). J Genet Genome Res. 2015;2(2):017 (Epub 2015 Sep 15).
Zhu Y, Stevens RG, Hoffman AE, Fitzgerald LM, Kwon EM, Ostrander EA, et al. Testing the circadian gene hypothesis in prostate cancer: a population-based case-control study. Cancer Res. 2009;69(24):9315–22.
Dai H, Zhang L, Cao M, Song F, Zheng H, Zhu X, et al. The role of polymorphisms in circadian pathway genes in breast tumorigenesis. Breast Cancer Res Treat. 2011;127(2):531–40.
Karantanos T, Theodoropoulos G, Gazouli M, Vaiopoulou A, Karantanou C, Stravopodis DJ, et al. Association of the clock genes polymorphisms with colorectal cancer susceptibility. J SurgOncol. 2013;108(8):563–7.
Li J, Humphreys K, Heikkinen T, Aittomaki K, Blomqvist C, Pharoah PD, et al. A combined analysis of genome-wide association studies in breast cancer. Breast Cancer Res Treat. 2011;126(3):717–27.
Truong T, Liquet B, Menegaux F, Plancoulaine S, Laurent-Puig P, Mulot C, et al. Breast cancer risk, nightwork, and circadian clock gene polymorphisms. EndocrRelat Cancer. 2014;21(4):629–38.
Zhao B, Lu J, Yin J, Liu H, Guo X, Yang Y, et al. A functional polymorphism in PER3 gene is associated with prognosis in hepatocellular carcinoma. Liver Int. 2012;32(9):1451–9.
Rodriguez S, Gaunt TR, Day IN. Hardy-Weinberg equilibrium testing of biological ascertainment for Mendelian randomization studies. Am J Epidemiol. 2009;169(4):505–14.
Dupont WD, Plummer WD Jr. Power and sample size calculations. A review and computer program. Control Clin Trials. 1990;11(2):116–28.
Yu K, Li Q, Bergen AW, Pfeiffer RM, Rosenberg PS, Caporaso N, et al. Pathway analysis by adaptive combination of P-values. Genet Epidemiol. 2009;33(8):700–9.
Mocellin S, Tropea S, Benna C, Rossi CR. Circadian pathway genetic variation and cancer risk: evidence from genome-wide association studies. BMC Med. 2018;16(1):20.
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.
Monsees GM, Kraft P, Hankinson SE, Hunter DJ, Schernhammer ES. Circadian genes and breast cancer susceptibility in rotating shift workers. Int J Cancer. 2012;131(11):2547–52.
Zienolddiny S, Haugen A, Lie JA, Kjuus H, Anmarkrud KH, Kjaerheim K. Analysis of polymorphisms in the circadian-related genes and breast cancer risk in Norwegian nurses working night shifts. Breast Cancer Res. 2013;15(4):R53.
Hunter DJ, Kraft P, Jacobs KB, Cox DG, Yeager M, Hankinson SE, et al. A genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancer. Nat Genet. 2007;39(7):870–4.
Petersen GM, Amundadottir L, Fuchs CS, Kraft P, Stolzenberg-Solomon RZ, Jacobs KB, et al. A genome-wide association study identifies pancreatic cancer susceptibility loci on chromosomes 13q22.1, 1q32.1 and 5p15.33. Nat Genet. 2010;42(3):224–8.
Markt SC, Valdimarsdottir UA, Shui IM, Sigurdardottir LG, Rider JR, Tamimi RM, et al. Circadian clock genes and risk of fatal prostate cancer. Cancer Causes Control. 2015;26(1):25–33.
Chu L, Zhu Y, Yu K, Zheng T, Yu H, Zhang Y, et al. Variants in circadian genes and prostate cancer risk: a population-based study in China. Pros Cancer Pros Dis. 2008;11(4):342–8.
Madden MH, Anic GM, Thompson RC, Nabors LB, Olson JJ, Browning JE, et al. Circadian pathway genes in relation to glioma risk and outcome. Cancer Causes Control. 2014;25(1):25–32.
Rana S, Shahid A, Ullah H, Mahmood S. Lack of association of the NPAS2 gene Ala394Thr polymorphism (rs2305160:G>A) with risk of chronic lymphocytic leukemia. Asian Pac J Cancer Prev. 2014;15(17):7169–74.
Wang F, Hu Z, Yang R, Tang J, Liu Y, Hemminki K, et al. A variant affecting miRNAs binding in the circadian gene Neuronal PAS domain protein 2 (NPAS2) is not associated with breast cancer risk. Breast Cancer Res Treat. 2011;127(3):769–75.
Zhu Y, Leaderer D, Guss C, Brown HN, Zhang Y, Boyle P, et al. Ala394Thr polymorphism in the clock gene NPAS2: a circadian modifier for the risk of non-Hodgkin’s lymphoma. Int J Cancer. 2007;120(2):432–5.
Zhu Y, Stevens RG, Leaderer D, Hoffman A, Holford T, Zhang Y, et al. Non-synonymous polymorphisms in the circadian gene NPAS2 and breast cancer risk. Breast Cancer Res Treat. 2008;107(3):421–5.
Yuan P, Wang S, Zhou F, Wan S, Yang Y, Huang X, et al. Functional polymorphisms in the NPAS2 gene are associated with overall survival in transcatheter arterial chemoembolization-treated hepatocellular carcinoma patients. Cancer Sci. 2014;105(7):825–32.
Zhang Z, Ma F, Zhou F, Chen Y, Wang X, Zhang H, et al. Functional polymorphisms of circadian negative feedback regulation genes are associated with clinical outcome in hepatocellular carcinoma patients receiving radical resection. Med Oncol. 2014;31(12):179 (Epub 2014 Oct 26).
Alhopuro P, Bjorklund M, Sammalkorpi H, Turunen M, Tuupanen S, Bistrom M, et al. Mutations in the circadian gene CLOCK in colorectal cancer. Mol Cancer Res. 2010;8(7):952–60.
Nakahata Y, Kaluzova M, Grimaldi B, Sahar S, Hirayama J, Chen D, et al. The NAD+-dependent deacetylase SIRT1 modulates CLOCK-mediated chromatin remodeling and circadian control. Cell. 2008;134(2):329–40.
Bersten DC, Sullivan AE, Peet DJ, Whitelaw ML. bHLH-PAS proteins in cancer. Nat Rev Cancer. 2013;13(12):827–41.
Hoffman AE, Zheng T, Ba Y, Zhu Y. The circadian gene NPAS2, a putative tumor suppressor, is involved in DNA damage response. Mol Cancer Res. 2008;6(9):1461–8.
Cao Q, Gery S, Dashti A, Yin D, Zhou Y, Gu J, et al. A role for the clock gene per1 in prostate cancer. Cancer Res. 2009;69(19):7619–25.
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.
Fan W, Chen X, Li C, Yongluo, Chen L, Liu P, et al. The analysis of deregulated expression and methylation of the PER2 genes in gliomas. J Cancer Res Ther. 2014;10(3):636–40.
Fonnes S, Donatsky AM, Gogenur I. Expression of core clock genes in colorectal tumour cells compared with normal mucosa: a systematic review of clinical trials. Colorectal Dis. 2015;17(4):290–7.
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.
Mostafaie N, Kallay E, Sauerzapf E, Bonner E, Kriwanek S, Cross HS, et al. Correlated downregulation of estrogen receptor beta and the circadian clock gene Per1 in human colorectal cancer. MolCarcinog. 2009;48(7):642–7.
Fu L, Lee CC. The circadian clock: pacemaker and tumour suppressor. Nat Rev Cancer. 2003;3(5):350–61.
Xiong G, Wang C, Evers BM, Zhou BP, Xu R. RORalpha suppresses breast tumor invasion by inducing SEMA3F expression. Cancer Res. 2012;72(7):1728–39.
Zhu Y, McAvoy S, Kuhn R, Smith DI. RORA, a large common fragile site gene, is involved in cellular stress response. Oncogene. 2006;25(20):2901–8.
Benna C, Simioni A, Pasquali S, De Boni D, Rajendran S, Spiro G, et al. Genetic susceptibility to bone and soft tissue sarcomas: a field synopsis and meta-analysis. Oncotarget. 2018;9(26):18607–26.
Zhang B, Beeghly-Fadiel A, Long J, Zheng W. Genetic variants associated with breast-cancer risk: comprehensive research synopsis, meta-analysis, and epidemiological evidence. Lancet Oncol. 2011;12(5):477–88.
Theodoratou E, Montazeri Z, Hawken S, Allum GC, Gong J, Tait V, et al. Systematic meta-analyses and field synopsis of genetic association studies in colorectal cancer. J Natl Cancer Inst. 2012;104(19):1433–57.
Mocellin S, Verdi D, Pooley KA, Nitti D. Genetic variation and gastric cancer risk: a field synopsis and meta-analysis. Gut. 2015;64(8):1209–19.
Yeager M, Orr N, Hayes RB, Jacobs KB, Kraft P, Wacholder S, et al. Genome-wide association study of prostate cancer identifies a second risk locus at 8q24. Nat Genet. 2007;39(5):645–9.
Reick M, Garcia JA, Dudley C, McKnight SL. NPAS2: an analog of clock operative in the mammalian forebrain. Science. 2001;293(5529):506–9.
Poggiogalle E, Jamshed H, Peterson CM. Circadian regulation of glucose, lipid, and energy metabolism in humans. Metabolism. 2018;84:11–27.
Brown SA. Circadian metabolism: from mechanisms to metabolomics and medicine. Trends EndocrinolMetab. 2016;27(6):415–26.
Bahrami-Nejad Z, Zhao ML, Tholen S, Hunerdosse D, Tkach KE, van Schie S, et al. A transcriptional circuit filters oscillating circadian hormonal inputs to regulate fat cell differentiation. Cell Metab. 2018;27(4):854–868.e8.
Garaulet M, Corbalan-Tutau MD, Madrid JA, Baraza JC, Parnell LD, Lee YC, et al. PERIOD2 variants are associated with abdominal obesity, psycho-behavioral factors, and attrition in the dietary treatment of obesity. J Am Diet Assoc. 2010;110(6):917–21.
Scott EM, Carter AM, Grant PJ. Association between polymorphisms in the Clock gene, obesity and the metabolic syndrome in man. Int J Obes (Lond). 2008;32(4):658–62.
Bill KL, Casadei L, Prudner BC, Iwenofu H, Strohecker AM, Pollock RE. Liposarcoma: molecular targets and therapeutic implications. Cell Mol Life Sci. 2016;73(19):3711–8.
Li B, Wang Y, Xu Y, Liu H, Bloomer W, Zhu D, et al. Genetic variants in RORA and DNMT1 associated with cutaneous melanoma survival. Int J Cancer. 2018;142(11):2303–12.
CB, SM: statistical analysis and manuscript writing; SR, GS: DNA extraction and genotyping; PDF, ST: clinico-pathological data management; CRR: appraisal of manuscript. All authors read and approved the final manuscript.
This research has been conducted using the resources of the BioBank of the 1st Surgical Clinic, University Hospital of Padova-Istituto Oncologico Veneto (IOV-IRCCS), Padova, Italy.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional files.
Consent for publication
Ethics approval and consent to participate
Padova University Hospital Ethic Committee prot#448.
University of Padova, BIRD168075, “Germline polymorphisms of candidate genes as predictor of risk and prognosis in patients with cutaneous melanoma and soft tissue sarcoma”.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1. Hardy-Weimberg equilibrium (HWE) assumption test for all the SNPs analyzed in this study in healthy controls and sarcoma patients.
Additional file 2: Table S2. Statistical power for all the SNPs tested in this study in sarcoma patients.
Additional file 3: Table S3. Upper panel: associations of circadian pathway genes with prognosis of sarcoma under 3 models of inheritance (additive, recessive, dominant); lower panel: association of circadian pathway genes with prognosis of liposarcoma and leiomyosarcoma in subgroups analysis under the additive genetic model.
About this article
- Soft tissue sarcoma
- Circadian clock gene
- Single nucleotide polymorphism (SNP)
- Genetic variation
- Circadian pathway