Germline polymorphisms of circadian genes and gastric cancer predisposition
Authors: Senthilkumar Rajendran, Clara Benna, Alberto Marchet, Donato Nitti, Simone Mocellin
Journal: Cancer Communications (2020)

Dear Editor, Gastric cancer represents a remarkable disease burden worldwide, ranking among the first five tumor types in incidence and mortality [1]. Germline DNA variation has been extensively investigated in terms of predisposition to sporadic gastric cancer, which represents more than 90% of all cases [2]. Currently available evidence shows that the fraction of disease burden that can be attributable to known risk polymorphisms is small (< 20%) [2]. Single germline variations of circadian genes (also called clock genes) have been associated with the predisposition of different tumor types [3]. The circadian clock is a time-tracking rhythmic biological system with a periodicity of about 24 hours that enables organisms to anticipate environmental changes and allow them to modify their behavior and physiological functions in the most efficient way. Circadian rhythms are controlled by proteins encoded by circadian genes, which have been discovered in all studied species. Remarkably, the disruption of these rhythms has been linked with risk of different diseases including cancer. In regards to the latter, a growing wealth of evidence supports the potential tumor suppressor role of the biological clock [3, 4]. As the role of circadian gene germline variants has never been explored in the field of gastric cancer susceptibility, with the present work, we intended to test the hypothesis that specific single nucleotide polymorphisms (SNPs) of the circadian genes, such as CLOCK, NPAS2, PER1, PER2, RORA, and TIMELESS, could significantly increase or decrease the predisposition to develop gastric cancer. We considered the 10 SNPs of the above listed 6 circadian genes that are known to be functional or associated with cancer risk or prognosis. The main features of the SNPs are described in our previous study [5]. We conducted a retrospective study based on a total of 1065 subjects comprising of 455 cases of gastric cancer and 610 healthy controls. All of them were of European ancestry. The median age of onset for gastric cancer was 67 years (range, 27-90 years). Among these gastric cancer patients, 249 (54.7%) were males and 206 (45.3%) were females. The median survival was 30.0 months, ranging from 1.0 to 293.0 months. These datasets were already employed in our previous studies [5, 6] and the detailed characteristics of the subjects are summarized in Table 1 and Supplementary Table 1. Genotyping was performed by real-time PCR. Multivariate logistic regression analysis was performed to assess the associations employing four models of inheritance: allelic, recessive, dominant, and co-dominant. The detailed methods are available in Supplementary information. All the preselected SNPs were successfully genotyped, and no departures from Hardy-Weinberg equilibrium were observed (Supplementary Table 2). The average genotyping success rate of selected SNPs in all participants was 98.9% (range, 96.0%-100%). The mean statistical power for this analysis was 61%. Detailed statistical power for each SNP is reported in Supplementary Table 3. Associations between the selected circadian genes genetic variations and gastric cancer predisposition were tested assuming 4 models of inheritance. The results are summarized in Table 2. We used odds ratios (ORs) and their corresponding 95% confidence intervals (CI) to measure the strength of association between each polymorphism and gastric cancer susceptibility. Overall, the genetic variants significantly associated with gastric cancer predisposition were: NPAS2 rs895520, PER1 rs3027178, PER2 rs934945, RORA rs339972. In particular, the present analysis suggested that NPAS2 rs895520 minor allele (A) was associated with an increased susceptibility to gastric cancer of 24% under an additive (per allele OR, 1.24; 95% CI, 1.01-1.52; P = 0.036), recessive (OR, 1.56; 95% CI, 1.09-2.24; P = 0.016) and co-dominant (OR, 1.62; 95% CI, 1.07-2.44; P = 0.022) model of inheritance. PER1 rs3027178, a genetic variant with a synonymous functional effect was associated with a reduced predisposition (per allele OR, 0.80; 95% CI, 0.64-0.99; P = 0.037). PER2 rs934945 (C > T) is located on the last exon of PER2 locus and has a missense functional effect, leading to the substitution of Glycine-Glutamic acid. Carriers of at least one copy of the minor allele had a decreased predisposition to develop gastric cancer (28%) employing a dominant genetic model (OR, 0.72; 95% CI, 0.53-0.98; P = 0.037). Employing a co-dominant model heterozygotes had a 31% risk reduction as compared to homozygotes for the common allele (C) (OR 0.69; 95% CI 0.50-0.94; P = 0.019). RORA rs339972 C allele was associated with a decreased predisposition to develop gastric cancer assuming an additive (per allele OR, 0.78; 95% CI, 0.63-0.98; P = 0.032) or dominant (OR, 0.75; 95% CI, 0.56-1.00; P = 0.049) genetic model. To the best of our knowledge, this is the first scientific work investigating the relations between circadian genes DNA genetic variations and the susceptibility to gastric cancer. Therefore, we could not know a priori the genotype-phenotype relation of these SNPs; as a consequence, we tested 4 genetic models of inheritance: allelic, recessive, dominant and co-dominant. When testing the allelic/recessive/dominant models, for those polymorphisms which were significantly associated with the phenotype in more than one model, the best fitting model was considered the one with the lower P value. Our results indicated that NPAS2 rs895520 best-fitted model for the association with gastric cancer was the recessive model of inheritance, while RORA rs339972 was the allelic model. Interestingly, we found similar results regarding NPAS2 rs895520 in our previous work on associations of circadian genes polymorphisms with soft tissue sarcoma susceptibility [5], while there was no difference in terms of P value for RORA rs339972 comparing the allelic and the dominant model, nevertheless, both were associated with sarcoma susceptibility as it was for gastric cancer. Since the maximum power was reached when the 'true' mode of inheritance of the disease susceptibility loci and the genetic model used in the analysis were concordant [7], it is worth determining the genotype-phenotype relation for each SNP. We tested the co-dominant model as well, for two reasons: its robust method [7] and its application in testing the circadian genes SNPs associations with different neoplasms [8, 9]. Employing the co-dominant model PER2 rs934945 heterozygotes had a decreased predisposition compared to homozygotes for the common allele (C) of 31%. Karantanos et al. [9] found no association of PER2 rs934945 with colorectal cancer neither with the allelic nor with the co-dominant model. Dai and colleagues [8] found no association of PER2 rs934945 with breast cancer in overall analysis while found a significant association in subgroup analysis. Homozygotes for the minor allele (T) had an increased risk of developing breast cancer only in a specific CLOCK rs3805151 background (homozygosis for the common allele C). This was in line with the shared idea that genetic variations have different effects in different neoplasms. In particular, this was recently highlighted for prognosis in an interesting work performed by Chang and Lai [4]. They performed a comprehensive study of circadian genes in 21 cancer types that considered genomic, transcriptomic and phenotypic (clinical prognosis) data and they found that circadian genes were substantially altered by somatically acquired deletions and amplifications. Core circadian genes, PERs, CRY2, CLOCK, NR1D2, RORA and RORB exhibited global patterns of somatic loss and downregulation across multiple tumor types and that loss-of-function of these genes resulted in increased death risks in patients. However, tumor suppressive qualities appeared to be cancer type-specific. Opposite trend was obtained for bladder and stomach cancers as their "low" loss-of-function of putative tumor-suppressive circadian genes were found to be associated with adverse survival outcomes [4]. In our previous study concerning the associations of gastric cancer prognosis and germline variation of circadian genes [6] we had a similar approach. We found that germline polymorphisms in the circadian pathway were associated with the survival of patients with gastric cancer, independently of established prognostic factors such as disease stage and patient age at diagnosis. In particular, combined information deriving from two SNPs (rs3749474 and rs1801260, two variants of the CLOCK gene 3'-UTR) allowed us to classify patients into a high or low CLOCK transcription, with the latter showing a significantly worse prognosis (about 70% increased risk of death). This apparent discrepancy highlights that gastric cancer prognosis and circadian genes relations need further in-depth analysis. Moreover, we could not replicate the data reported by Qu and colleagues [10] on the association between PER variants and prognosis. Different ethnicity (European vs. Asian), sample size (the Asian series was more than two-fold larger) and disease stage composition (only our study included patients with advanced and metastatic gastric cancer) might partly explain this discrepancy. Nevertheless, differences were found by 2 groups studying PER2 expression as a prognostic factor for gastric cancer in patients with Asian ethnicity. Zhao and colleagues [11] found that PER2 expression was downregulated in most gastric cancer tissues, while Hu and colleagues [12] found that it was upregulated. To our knowledge, this is the first analysis investigating the hypothesis of an association between germline genetic variations of the circadian pathway with gastric cancer susceptibility. The power of our study is not optimal, and the present study should be considered as a pilot work that warrants further validation in different datasets. Nevertheless, our results showed that the 4 circadian clock variants were clinically and statistically associated with gastric cancer predisposition. This study was approved by the Ethics Committee of Padova University Hospital (identifier: prot#448). Written informed consent was obtained from all patients. Not applicable. The authors declare that they have no competing interests. This work was supported by the University of Padova (DOR1944742). C.B. and S.M.: analyzed data and co-wrote the manuscript. S.R.: performed experiments. A.M.: managed clinico-pathological data. D.N.: provided critical revision of the manuscript. All authors read and approved the final manuscript. We thank the personnel of the Biobank of the First Surgical Clinic (Padova University Hospital, Padua, Italy) and particular Dr. Andrea Ferron (Padova University Hospital, Padua, Italy) for organizing the sampling activity and Dr. Enrico Lion (Padova University Hospital, Padua, Italy) for organizing the informed consent retrieval. All data generated or analyzed during this study are included in this published article and its additional files. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.