Abstract
Objective:
Pancreatic cancer (PC) is a serious disease with poor outcomes, and its prevalence has been increasing steadily. The circadian rhythm (CR) is involved in multiple physiological events and maintains homeostasis. Alterations in the CR elevate the risk of developing cancer. The present case-control research was carried out to estimate the possible association between PERIOD2/PERIOD3 (PER2/PER3) gene variable number tandem repeat polymorphism (VNTR) variants and PC in the Turkish population.
Materials and Methods:
A total of 198 subjects (78 patients with PC and 120 healthy controls) were enrolled in this work. Genomic DNA was collected from peripheral blood mononuclear cells, and genotypic analyses was performed using a polymerase chain reaction (PCR) method. Odds ratio (OR) with a 95% confidence interval (95% CI) was calculated using the χ2 test.
Results:
The frequency of the 4R (4 repeats)/3R (3 repeats), 3R/3R genotypes, and 3R allele of PER2 VNTR in patients with PC was significantly higher than in the control group (p = 0003, p = 0.00004, respectively). PER2 VNTR 4/5 genotype was related to perineural invasion (p = 0.040). The genotype and allele distribution of PER3 VNTR variant did not show any statistical difference between the two groups (p > 0.05). The PER2/PER3 VNTR 4/5-4R/3R combined genotype was increased in the patient group (p = 0.013), while 4/5-4R/4R combined genotype was increased in the control group (p = 0.0001).
Conclusions:
Our work has indicated that PER2 VNTR 3R allele may play a crucial role in the pathogenesis of PC in Turkish patients, which may become a useful marker for predicting the development of PC. Furthermore, the PER2 VNTR genotype seems to be related to perineural invasion in PC.
Introduction
Pancreatic cancer (PC) is a serious disease with poor outcomes, and its prevalence has been increasing steadily. PC is among the major causes of cancer-related mortality, ranking fourth in the United States and fifth in Europe (Sonoyama et al., 2011). Pancreatic adenocarcinoma and its related types account for 90% of all PCs (Feldmann et al., 2007). Although there have been marked advances in this field, PC remains to be one of the most aggressive and therapy-resistant types of malignancy, with the lowest 5-year survival rate of all types of malignancies (Jemal et al., 2009).
The poor outcome is mainly due to the delayed recognition of the disease along with the high invasiveness and drug resistance of the tumor. The current advancements in the investigation of the molecular and genetic changes in PC have not yet been reflected into significant improvement in patient survival or reduced mortality. There is a crucial need for new treatment targets and molecular markers in PC management to attain higher treatment efficiency. The cause of sporadic PC remains unclear. However, growing evidence reveals that pancreatic carcinogenesis involves a complex interaction between genetic mutations, epigenetic changes, and environmental risk factors (Bardeesy and DePinho, 2002).
Some biochemical pathways can result in carcinogenesis. Among these circadian rhythms (CRs) in physiological and behavioral processes take place. Numerous works have proved that circadian clock gene deregulation plays a task in the occurrence of cancer and other diseases. Impairment in CRs enhances tumor progression, and restoring its restoration is expected to improve prognosis.
The CR organizers aim clock-controlled genes to produce physiological clock at all essential cellular metabolism levels in both suprachiasmatic nucleus (SCN) neurons and peripheral organs, leading to a molecular rhythmicity that expresses 3-10% of all mRNAs in tissues, that is due to time-dependent communication between the circadian modulators with certain gene promoter sequences, transcription factors, or transcriptional initiation, elongation, and termination complexes, as well as the main factors regulating chromatin remodeling (Duong et al., 2011; Padmanabhan et al., 2012). Clock-controlled genes that were investigated include the essential regulators of cell proliferation, metabolism, senescence, and DNA damage response (Fu and Kettner, 2013).
Both lower expression and polymorphism of the core circadian genes PER1, PER2, and PER3 are often encountered in the human mammillary gland, uterus, prostate, pancreatic, colorectal, and NSCLC, as well as hepatocellular carcinoma, squamous cell carcinoma of the head and neck, glioma, acute myeloid lymphoma, and chronic myeloid lymphoma compared to normal host tissues (Fu and Kettner, 2013).
The circadian system regulates various biological functions during 24 h, including sleep/activity, body temperature, heart rate, blood glucose level, cortisol synthesis, and oxidative stress (Valenzuela et al., 2016). In mammalian animals, this system involves a central clock found in the SCN of the hypothalamus and in a series of peripheral oscillators like the liver, lungs, adrenal glands, fibroblast cells, and other tissues (Valenzuela et al., 2016). Some clock genes, such as PER1, PER2, PER3, CRY1, CRY2, BMAL1, CLOCK, and casein kinase 1 delta/epsilon (CK1D and CK1E), interact in translational-transcriptional autoregulatory feedback loops to produce CRs (Hirayama and Sassone-Corsi, 2005).
Any alteration in the Clock gene regulation can affect any part of the organism, including the cardiovascular, gastrointestinal, endocrine, central nervous systems (CNSs), sleep-wake cycles, and cellular processes that are very important in carcinogenesis (e.g., cell proliferation, DNA damage response, and apoptosis) (Gery and Koeffler, 2010). The PER (period circadian clock) gene family is deemed to be the core group member of the circadian clock genes family, which possesses the PER1/2/3 genes as subgroups. Changes in PER gene manifestation have been implied in the genesis and prognosis of cancer (Zhao et al., 2016).
PER1/2 expression levels have been stated to be considerably suppressed in PC (Relles et al., 2013), gastric (Zhao et al., 2014), colorectal (Karantanos et al., 2013; Momma et al., 2017), prostate (Chu et al., 2008; Cao et al., 2009), breast cancer (Winter et al., 2007; Dai et al., 2011), head and neck squamous cell carcinoma (HNSCC) (Rahman et al., 2019), nonsmall cell lung cancer (NSCLC) (Liu et al., 2014), chronic lymphocytic leukemia (Rana et al., 2014), and melanoma (De Assis et al., 2018). It was reported that PER3 expression is considerably suppressed in PC, colorectal cancer, HNSCC, NSCLC, and hepatocellular carcinoma (Deng and Yang, 2019).
Because CR alteration is a known risk factor for several cancers, we aimed to determine the relationship of PER2/PER3 gene variable number tandem repeat polymorphism (VNTR) variants with development, clinical, and pathological characteristics of PC in a Turkish population using polymerase chain reaction (PCR) method for PER2/PER3 VNTR variant genotyping (Shumay et al., 2012).
Materials and Methods
Study population
The patient group was composed of 78 unrelated patients (49 males and 29 females; mean age (±SD) = 66.8 ± 9.9 years) with a clinical diagnosis of PC recruited from those who were treated and followed up in the Department of Medical Oncology, Education and Research Hospital, Gaziosmanpasa University, Tokat, Turkey. The diagnosis of PC was confirmed by immuno-histo-pathological examination of biopsy or resected tissue specimens as reported by the World Health Organization classification. A total of 120 unrelated healthy subjects (71 male and 49 female; mean age (±SD) = 65.4 ± 8.4 years) were enlisted as control group. Controls were matched with patients by sex and age. All controls were confirmed to be free of any malignancy. All cases and controls were Turkish.
The demographical and clinical characteristics (age, gender, smoking, age at diagnosis, tumor diameter, tumor localization, tumor subtypes, tumor differentiation, number of metastatic lymph nodes, and perineural invasion) of the patients were obtained using questionnaires and medical records. Written informed consent was obtained from all patients and healthy controls. The study protocols were performed in accordance with the principles of the Declaration of Helsinki. The experimental study protocol and process were assessed and approved by the Local Research Ethics Committee of Tokat Gaziosmanpaşa University, registered under the number 18-KAEK-262.
Genotyping
The genomic DNA was isolated from peripheral lymphocytes using DNA Isolation Kit (PureLink Genomic DNA Mini Kit; Invitrogen). The isolated DNA was used for description of the following polymorphic PER2/PER3 genes. PCR method was used for PER2/PER3 VNTR variant genotyping (Shumay et al., 2012). The studied VNTR region chr2: 239 184 578-239 185 009 has been defined by Shumay et al. (2012) and that rs10462023, rs13391269, and rs2304673 single nucleotide polymorphisms (SNPs) have been defined in regions close to this VNTR.
Statistical analysis
Data were analyzed using the Statistical Package for Social Sciences (SPSS) software version 20.0 for Windows (SPSS, Inc., Chicago, IL). Mean and standard deviation were used for the presentation of continuous quantitative variables. Frequencies and percentages were used for categorical data. The PER2/PER3 VNTR overall genotype distribution was compared by chi-square (χ2) test, and the specific genotype and allele distributions were compared using Fisher's exact test. The odds ratios (ORs) and 95% confidence intervals (CIs) were used to determine the relationships between the PER2/PER3 VNTR allelic and genotypic variants and their occurrence in the patients. A p-value smaller than 0.05 was considered statistically significant.
Sample size estimations (Priori) and power analysis (Post hoc)
The sample size was calculated for Chi-square test, which was used to test the main hypothesis of our study. As a result of the sample size analysis conducted using previous study knowledge (Relles et al., 2013), it was found that minimum 150 individuals needed to be involved in the study in order to reveal the significant differences in the groups using 80% power (1-β = 0.80) and α = 0.05 error (95% confidence interval). Post hoc power analysis was performed to determine the power of the study with a type 1 error value of 0.05 for the basic hypotheses that were found to be statistically significant. The G*power (version 3.1.9.7) package was used for sample size estimation and post hoc power analysis.
Results
A total of 78 PC patients and 120 age- and gender-matched healthy controls were genotyped for PER2/PER3 gene variants. Demographic and clinical characteristics of the study participants are shown in Table 1.
Baseline Characteristics of Pancreatic Cancer Patients and Controls
Data were analyzed by analysis of variance and χ2 test. Mean plus standard deviation values are presented for age, age of diagnosis, and tumor diameter.
PC, pancreatic cancer; IPMN, intraductal papillary mucinous neoplasm.
The distribution of genotypes and allele frequencies of PER2/PER3 VNTR variants in the patient and control groups are shown in Table 2 and Table 3. There was a significant difference in genotype and allele frequencies of PER2 VNTR variant between patients and control subjects. The power for genotypes was 96.8%, and power for alleles was 98.7%. PER2 VNTR 4R (4 repeats)/3R (3 repeats) and 3R/3R genotype were more common in patient group (p = 0003). In addition, 3R allele increased in patient group (p = 0.00004). The genotype and allele distribution of PER3 VNTR variant did not show any statistical difference between patients and controls (p > 0.05).
Genotype and Allele Frequencies of PER2 Variable Number Tandem Repeat Polymorphism Variant in Pancreatic Cancer Patients and Controls
Data were analyzed by χ2. The results that are statistically significant are typed in bold.
PER2, period circadian clock 2; VNTR, variable number tandem repeat polymorphism.
Genotype and Allele Frequencies of PER3 Variable Number Tandem Repeat Polymorphism Variant in Pancreatic Cancer Patients and Controls
Data were analyzed by χ2.
PER3, period circadian clock 3.
Furthermore, we also analyzed whether any differences existed in the clinical and pathological features of patients based on the genotype distributions of PER2/PER3 VNTR variants (Table 4). There was a significant difference between the PER2 variant and a clinical feature. PER2 VNTR 4/5 genotype was related with perineural invasion (p = 0.040).
Clinical and Pathological Characteristics of Pancreatic Cancer Patients Stratified According to PER3 and PER2 Variable Number Tandem Repeat Polymorphism Variants
Data were analyzed by Anova or χ2 tests.
The results that are statistically significant are shown in boldface.
PC, pancreatic cancer; PER2, period circadian clock 2; PER3, period circadian clock 3.
In addition, we studied the risk related to the combined genotypes for the PER2 and PER3 VNTR variants (Table 5). Our results also revealed a significant difference between the combined genotypes in the patient and control groups. PER3/PER2 VNTR 4/5-4R/4R combined genotype was higher in the control group (p = 0.0001). 4/5-4R/3R combined genotype increased in the patient group (p = 0.013).
Comparative Analysis of Combined Genotypes of Patients and Controls
Data were analyzed by χ2 or fisher's exact tests. The results that are statistically significant are typed in bold.
Discussion
In this article, we investigated whether the PER2 VNTR (in intron 3) and PER3 VNTR variants are related to the risk of PC and evaluated its relationship with clinical and pathological findings. In most trials studying the task of the molecular clock in human cancers, disturbance or variations of several or all whole circadian genes are encountered. The PER2 gene belongs to the period family type of genes containing PER1, PER2, and PER3 and is mostly found in the CNS, including the SCN and the peripheral nervous systems. The PER2 is localized on chromosome 2 at 2q37.3 (Toh et al., 2001). The PER2 gene is vital in controlling the CR and plays a crucial function in tumor suppression (Oda et al., 2009).
PER2 mutant mice exhibited a shorter circadian period compared to wild type (WT) mice and diminished PER1 expression in the SCN, implying that PER2 controls PER1 (Kim et al., 2018). Downregulation of PER2 enhances β-catenin protein levels and its target cyclin D, resulting in cell proliferation in colon cancer cell lines and intestinal and colonic polyp development. This indicates that the PER2 gene product suppresses tumorigenesis in the small intestine and colon by downregulation of β-catenin and β-catenin target gene signaling pathways (Wood et al., 2009). Other article stated that overexpression of mPER2 in mouse tumor cells caused apoptotic cell death and the suppression of tumor growth (Hua et al., 2006). Relles et al. (2013) found that expression levels of PER1, PER2, PER3, Cry1, Cry2, Tipin, Tim, CK1e, Bmal-ARNTL, and Clock genes were reduced in the pancreatic ductal adenocarcinoma tissue in comparison with their matched nearby tissue.
The PER2 gene can turn on the c-Myc signaling pathways producing a genomic imbalance and cell growth. PER2 malfunction can also disrupt p53-mediated apoptosis and cause genomic instability and the accumulation of deficient cells (Benna et al., 2018). The PER2 locus has elevated levels of interpopulation genetic differentiation relative to other loci, indicating a task for geographically limited positive selection (Cruciani et al., 2008).
The SNPs in the PER2 gene have been linked with aberrant circadian parameters, chronotypes, depression, and also with increased alcohol intake; however, no correlation with cocaine dependence was encountered in a study that investigated some PER2 SNPs (Shumay et al., 2012). It was declared that PER2 rs934945 variant T allele was associated with reduced predisposition to soft tissue sarcoma (Benna et al., 2018). Although these VNTRs are often classified and neglected as “junk” DNA, the number of repetitions of the sequence can differ within and between individuals, thus making these VNTRs a polymorphic entity.
The PER3 VNTR (rs57875989) consists of four or five copies of a 54-bp sequence encoding 18 amino acids (Nadkarni et al., 2005). The 5-repeat variant supplements some potential phosphorylation motifs to the gene, and PER3's that interoperate with circadian activity could be increased in those people (Wirth et al., 2013). The 5-repeat PER3 allele is linked with a rather penetrant phenotype that has morning circadian preference, overt cognitive decline as a reaction to sleep deficiency, discrepancies in rate or schedule of melatonin or cortisol secretion, and a propensity for depressive manifestations or an earlier commencement of bipolar disorder (ALExANDER et al., 2015).
In a study, it was found that the PER3 5/5 variant carriers were related to increased amount of interleukin-6, B-type natriuretic peptide, and lower vitamin A levels (Lipkova et al., 2014b). Golalipour et al. (2017) reported that PER3 4/4 genotype and 4-repeat allele were connected with multiple sclerosis. No significant difference was observed in allele periodicities of PER3 VNTR between the obese and the nonobese individuals (Bienertova-Vasku et al., 2014). In addition, there was no significant difference between chronic heart failure patients and controls regarding PER3 VNTR allele and genotype frequencies (Lipkova et al., 2014a).
Geng et al. (2015) found a moderate increase in risk of cancer among individuals with the PER3 VNTR 5 allele compared to individuals with the 4 allele in a meta-analysis. Alexander et al. (2015) indicated that the individuals carrying the 5-repeat PER3 VNTR variant may be more susceptible to colorectal adenoma formation. Yeğin et al. (2020) reported that PER3 VNTR variant was not associated with bladder carcinoma. But they found that the survival times of patients decreased in the patient group (progression and cystectomy positive) for PER3 4/4 genotype and (recurrence, progression, and cystectomy positive) for PER3 4/5 genotype.
In this study, the statistical analyses showed that genotype and allele distribution of PER2 were statistically significantly different between the control and patient groups. PER2 VNTR 4R/3R, 3R/3R genotype and 3R allele increased in patient group (p < 0.05) (Table 2). So PER2 VNTR 3R allele seems to be a risk factor in PC formation. Our study was unable to show any significant association of the PER3 VNTR genotype and allele distribution with PC (p > 0.05) (Table 3). Then, we also analyzed if any differences existed in the demographical and clinical characteristics of patients according to allele frequencies. We found that PER2 VNTR 4/5 genotype increased in patients with perineural invasion (p < 0.05) (Table 4). In PER3/PER2 combined analysis, 4/5-4R/3R combined genotype was more prevalent in patient group (p < 0.05) (Table 5).
There are some limitations of this study. First, we investigated only two variants involved in the pathway of PER2/PER3; other regulatory genes in the clock gene family signaling pathway may also play a role in the genesis of PC. Second, even though the PER3 VNTR variant that we examined was not found to be associated with PC, other variants of PER2/PER3 gene have not been examined yet. Third, owing to the rather small sample size, the frequencies of some homozygous variants were low in groups and consequently diminished the statistical power. Finally, absence of assessment of expression levels of PER2/PER3 is also a limitation of this study.
Conclusion
In conclusion, our work has indicated that PER2 VNTR 3R allele may play a significant task in the development of PC in Turkish patients, which may become a useful marker for predicting the development of PC. Furthermore, the PER2 VNTR genotype seems to be associated with perineural invasion in PC. These findings imply that continued research into clock gene variants will be an important source of information in understanding the pathogenesis of PC. Further multicenter, well-designed studies with larger populations that evaluate gene-environment interactions are necessary to endorse our results.
Data Availability Statement
The data sets generated during the current study are available from the corresponding author upon request.
Statement of Licensing Committee and Institution
Ethics committee approval for this study was obtained from the Local Research Ethics Committee registered under the number 18-KAEK-262. All methods were performed in accordance with the relevant guidelines and regulations of the institution. Written informed consent was obtained from subjects and patients who participated in this study.
Footnotes
Authors' Contribution
H.D.: Conceptualization, Software, Validation, Formal analysis, Data Curation, Writing—Review & Editing, Funding acquisition.
S.Y.: Methodology, Software, Validation, Formal analysis, Writing—Review & Editing, Visualization, Supervision.
A.F.N.: Methodology, Software, Validation, Resources, Writing—Original Draft, Supervision, Project administration.
E.D.: Validation, Formal analysis, Data Curation, Writing—Review & Editing, Supervision.
O.G.: Conceptualization, Validation, Formal analysis, Resources, Data Curation, Writing—Review & Editing, Project administration.
Author Disclosure Statement
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this article. The author confirms on behalf of the other authors that there is no conflict of interest to declare.
Funding Information
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
