1. Introduction
The circadian timing system consists of the brain’s central clock in the hypothalamic suprachiasmatic nucleus (SCN) and various clocks in peripheral tissues, including muscle, adipose tissue, and liver, that orchestrate 24-hour cycles at the transcriptional level [1,2]. The SCN receives a 24-hours solar day via a pathway from the retina that synchronize the circadian rhythm of SCN [2]. The circadian clock is implicated in the regulation of numerous cellular, physiological, metabolic, endocrinal, and behavioral systems. The circadian system regulates energy expenditure, food intake, body insulin sensitivity, and the sleep–wake cycle [3,4].
In mammals, the central circadian rhythm is regulated by an autonomous transcription–translation feedback loop of cells that provides a pattern of rhythmic expression to the clock genes. During the day, the Circadian Locomotor Output Cycles Kaput (CLOCK) transcription factors and Neuronal PAS Domain Protein 2 (NPAS2) form complexes with Brain and Muscle ARNT-Like 1-2 (ARNTL or BMAL1 and ARNTL2 or BMAL2) in order to activate the transcription of Period Circadian Clock (PER1/PER2) and Cryptochrome Circadian Clock (CRY1/CRY2). The heterodimerized proteins Period Circadian Regulator (PER) and Cryptochrome Circadian Regulator (CRY) move to the nucleus and suppress the action of CLOCK/BMAL1, thus forming a negative feedback loop [5]. Moreover, recent genome-wide association studies GWASs have associated the Vaccinia-related Kinase (VRK) Serine/Threonine Kinase 2 (VRK2) gene with the circadian cycle and, more specifically, with sleep [6,7] although the mechanism has yet to be explained.
It has, however, been shown that alterations in lifestyle patterns due to social or occupational patterns can contribute to the disruption of the circadian cycle [1,4], with possible effect on metabolic risk factors [8]. Furthermore, genetic differences may also be associated with inter-individual factors associated with the circadian system [9]. Thus, gene expression analysis may be a potential method for identifying molecular pathways associated with cardiometabolic diseases [10,11]. Our group is well-acquainted with the importance of genes involved in the circadian cycle, having observed a significant association between the CLOCK-rs4580704 single-nucleotide polymorphism (SNP) and type 2 diabetes (T2D) incidence in 3671 non-T2D Prevention with Mediterranean Diet (PREDIMED) study participants, suggesting a gene effect modulation induced by the Mediterranean diet [3]. In addition, we recently observed that the association between genetic variations in the melatonin receptor 1B gene (MTNR1B) consisting of the rs10830963 SNP was age-dependent, having an increased effect in young populations in comparison with older populations [12]. However, as far as we know, very few studies have analyzed the association between gene expression and clock DNA polymorphisms or sleep-related variables in humans. Therefore, our aim was to analyze the effect of selected DNA-SNPs in the main clock genes on leukocyte gene expression. In addition, we analyzed the modulation by sex and other lifestyle variables as well as their association with sleep-related variables.
2. Materials and Methods
2.1. Study Subjects
We carried out a cross-sectional study in a subsample of subjects from a general Mediterranean population participating in the Obesity, Nutrition & Information and Communication Technologies (OBENUTIC) study. Briefly, OBENUTIC is a case-controlled study of more than 1000 Caucasian participants aged between 18 and 80 years old from the Valencia Region, Spain. Inclusion criteria have been described previously [13]. We randomly selected 130 participants for the clock gene expression sub-study. Finally, we included a total of 102 participants for gene expression analysis after discarding those that did not meet RNA quality. The study was undertaken at the Department of Preventive Medicine and Public Health, School of Medicine at the University of Valencia, Valencia. Participants provided written informed consent, and study protocol and procedures were approved according to the ethical standards of the Helsinki Declaration and by the Human Research Ethics Committee of the University of Valencia, Valencia (reference number: H1488282121722; 6 April 2017).
2.2. Lifestyle and Clinical Evaluation
Socio-demographic and lifestyle variables (smoking, diet, and physical activity) were obtained through a standardized questionnaire previously used in our studies. In addition, we gathered information on various weekday and weekend sleep-related items. Participants were asked: (1) How many hours do you sleep on weekdays? (2) How many hours do you sleep on Saturdays? (3) How many hours do you sleep on Sunday? In addition, we asked the same questions, but this time related to the time of awakening. From these data, we calculated the following variables: Average sleep duration (h) = ((5 × sleep duration on weekdays) + sleep duration on Saturdays + sleep duration on Sunday)/7; Average waketime= ((5 × awakening time on weekdays) + awakening time on Saturdays + awakening time on Sunday)/7. Average weekend waketime was calculated with the previous formula but using only Saturday and Sunday data.
Participants were measured for height with a standard stadiometer (SECA Mod 220, SECA Deutschland Gmbh & Co, Hamburg, Germany). Weight was determined by bioimpedance using a TANITA-BC-420-S scale (Tanita Corporation 1-14-2 Maeno-Cho, Habashi-ku, Tokyo, Japan). Waist circumference was measured at the midpoint between the lowest rib and the iliac crest after a normal exhalation (ISAK, 2011). Systolic and diastolic blood pressure was obtained using an automatic sphygmomanometer (OMRON HEM-705CP). Two blood pressure readings were taken in sitting position, resting 5 min between each measurement and considering the level of blood pressure as the mean of the readings. Blood samples were obtained after an overnight fast. Fasting glucose and lipids were measured using a Beckman Coulter Olympus AU5400 (Pasadena, CA, USA).
2.3. DNA Extraction and Genotyping
Genomic DNA was isolated from white blood cells. The CLOCK-rs4580704, ARNTL-rs1982350, and MTNR1B-rs10830963 were genotyped with an ABI Prism 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA) using a fluorescent allelic discrimination TaqManTM assay. Genotype frequencies did not deviate from Hardy–Weinberg equilibrium expectations.
Besides the genotyping of these pre-selected polymorphisms in the clock genes, a high-density genotyping at the genome-wide level using the Illumina Infinium OmniExpress-24 BeadChip genotyping array (ver. 1.0 and 1.1) or the Illumina Global Screening Array (GSA) (ver. 1.0) (Illumina Inc., San Diego, CA, USA), was then carried out as previously described [14]. Allele detection and genotype calling were performed in the GenomeStudio genotyping module (ver. 2.0.5) (Illumina, Inc., San Diego, CA, USA), and quality-control procedures were applied. We then selected the main clock genes (CLOCK, ARNTL1, PER1, PER3, and VRK2) and extracted all the genotyped SNPs in these genes using a dedicated Python script. As the SNP coverage of the OmniExpress-24 array is different from that of the GSA array, we only analyzed SNPs that matched in both arrays in order to only use directly genotyped SNPs instead of the imputed ones. We tested the linkage disequilibrium between the available SNPs in the clock genes, and only independent SNPs were used in further analyses.
2.4. RNA Isolation, cDNA Synthesis, and qPCR
Total RNA was extracted from plasma samples using Trizol Reagent, employing the manual method described by Chomczynski P and Sacchi N in 1987 [15]. RNA was quantified by measuring absorbance at 260 and 280 nm using a NanoPhotometerTM spectrophotometer Implen P300 (Implen, Munich, Germany), and integrity was checked using an Agilent 2100 Bioanalyzer (Agilent Technologies, Foster City, CA, USA). Briefly, 500 ng total RNA from all samples were reverse transcribed using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystem, Foster City, CA, USA) and GeneAmp Polymerase Chain Reaction PCR System 9700 RT-PCR (Applied Biosystem, Foster City, CA, USA). Quantitative real-time PCR was performed using an ABI PRISM 7900 HT Sequence Detection System (Applied Biosystem). The genes were amplified using commercially available Taqman probes (Applied Biosystems). The Taqman Gene Expression kit was used for real-time PCR analysis. Duplicate amplification reactions were performed under the following conditions: 50 °C for 2 min, 95 °C for 10 min, and 40 cycles of 95° for 15 s and 60 °C for 1 min. The software used for the analysis of the results was Sequence Detection Systems (SDS) Automation Controller Software ver. 2.3 and Relative Quantitation (RQ) manager ver. 1.2 (Applied Biosystem, Foster City, CA, USA). The relative mRNA level for each transcript was calculated by the methods proposed by Livak KJ et al. for gene expression. Gene expression values at baseline are shown as 2-ΔCt [10,16].
2.5. Statistical Analysis
Anthropometric, blood pressure, and biochemical measurements were expressed as mean ± SD (standard deviation). The results for genes expression are presented as mean ± SEM (standard error of mean). Variables that did not fulfil normal distribution criteria were Ln transformed to improve symmetry for subsequent analyses. A polygenic risk score (PRS) for each clock gene was created by the unweighted sum of the highly significant clock SNP related with circadian phenotypes. Individual participant scores were created by adding up the number of risk alleles at each SNP. To test the differences between clock genes expression and lifestyle or clinical variables, we used linear regression tests. Linear regression coefficients were used to analyze the associations between SNPs genotypes and genes expression. Multivariate models were carried out to adjust for confounding factors. Statistical analyses were performed with IBM SPSS Statistics, version 27.0 (Armonk, NY, USA). All tests were two-tailed, and p-values < 0.05 were considered statistically significant.
3. Results
3.1. Population Characteristics
The 102 participants included in the analysis, 52 being women, had an average age of 43.1 ± 13.8 years and had a Body mass index (BMI) of 27.2 ± 5.1 kg/m2. Clinical and lifestyle characteristics are represented in Table 1. Briefly, 30.4% of the whole population presented obesity, 25.7% were hypertensive, 16.5% were current smokers, 33% had hypercholesterolemia, and 5% had diabetes. Differences per sex were detected in the prevalence of hypertension, obesity, and diabetes, which were higher in men (p ˂ 0.05). Individual clinical parameters also showed the same differences between males and females (Table 1).
3.2. Circadian Clock Gene Expression
Results of the global peripheral blood expression genes are represented in Figure 1. Women showed higher clock genes expression than men in all the genes analyzed, but this difference was only significant for PER1 expression (p = 0.034).
No differences in these gene expressions were found according to age, obesity, diabetes, diet, sleep (duration and waketime), or smoking. However, in this last lifestyle factor, we detected a lower VRK2 expression in smokers, but it did not reach statistical significance (β: −1.30 ± 0.70; p = 0.067), Table S1.
3.3. Impact of the Selected SNPs on Clock Gene Expression
We analyzed the differences in clock genes expression in all clock genes using the Infinium Ommni Express Array and the Infinium Global Screening Array methods. This netted a total of 32 polymorphic sites in the CLOCK, ARNTL, PER1, PER3, and VRK2 genes. We also analyzed 3 more TaqIB SNP in the CLOCK, ARNTL, and MTNR1B genes, with a total of 35 SNPs analyzed. Of these 35 polymorphisms, 1 was a duplicate (rs1982350), and a total of 7 SNPs did not reach the Hardy–Weinberg equilibrium expectations at p < 0.05, so we removed it from our analysis. Of these 28 polymorphisms, 4 polymorphisms (rs1982350, rs7924734, rs10832027, and rs2678902) were associated with CLOCK gene expression. Two more SNPs (rs228642 and rs10127838) in the PER3 gene were associated with PER1 expression, and for the VRK2 gene, three polymorphisms were linked (rs10832027, rs11022778, and rs10830963) with its expression (Table 2).
Furthermore, we tested whether these polymorphisms were related to sleep factors, and we found an association between the rs238666 in the PER3 gene (β: 0.354 ± 0.354, p = 0.008) and the rs4580704 in CLOCK gene (β: −0.360 ± 0.155, p = 0.023) with weekly sleep duration after adjusting for sex and age. Similar results were detected for the two SNPs in the VRK2 gene and waketime on the weekends (rs1051061: β: 0.435 ± 0.187, p = 0.023; and rs848291: β: −0.437 ± 0.187, p = 0.022.
In addition, we created a specific PRS with the most relevant SNPs for each clock gene. Thus, we detected that the three PRS could be used as a predictor of clock genes expression: CLOCK expression = β: 0.516 ± 0.135, p < 0.001; PER1 expression = β: 0.867 ± 0.248, p = 0.001, and VRK2 expression = 0.508 ± 0.127, p < 0.001 (Table 2).
4. Discussion
We analyzed gene expression levels in leukocytes of the main clock-related genes in a Mediterranean general population. To our knowledge, this study provides the first evidence of sex-differences expression of core molecular clock components in human leukocytes from blood. We also found, for the first time, differences in clock gene expressions according to different genotypes from circadian genes (ARNTL, PER3, and VRK2). However, due to our sample size, these findings need further investigation in a larger population and in different cohorts.
There is strong evidence of sex differences in physiology, health, and diseases. These differences have been attributed to hormones, sex chromosomes, genotypes, behavior, or environmental exposures [17]. To date, the underlying mechanisms continue to be unknown. The Genotype-Tissue Expression (GTEx) project has evaluated sex gene regulation in 29 human healthy tissue using 8279 whole-genome expression profiles to understand the biological mechanism basis of sex differences [18]. They found that the number of genes differentially expressed by sex is small (around 64 genes) in some tissues apart from the breast [18]. Our sex-specific analyses suggest that females present a higher clock gene expression in all the genes analyzed. However, statistical significance was achieved only for the PER1 gene. Similar results were detected by Gómez-Abellán in obese women in the expression of PER2, ARNTL1, and CRY1 in adipose tissue [19]. Lim et al. found a timing expression of three clock genes (PER2, PER3, and ARNTL1) earlier in females than in men [20], and this also could explain the differences per sex identified in our sample. Our findings may confirm the growing evidence of sex differences in circadian rhythms, such as intrinsic period, chronotype, or differences in sleep phases [20,21], and according to some authors, estrogens could partly modulate circadian clock genes expression [22]. On the other hand, although we did not find any significant differences in clock gene expression by lifestyle and clinical variables, other authors detected that clock genes expression could be influenced by other factors such as sleep disturbances and chronic pathologies. Emekli et al. noticed a decrease in the expression levels of PER1 and PER2 genes in night shift workers and patients with chronic insomnia [23]. One week of insufficient sleep (5.70 h) was linked with the expression of circadian rhythms, mainly involving the PER2 and PER3 genes, two of the most down-regulated genes after sleep deprivation intervention [24]. Furthermore, studies in mice observed that 28 days of sleep restriction was associated with the disappearance of circadian rhythms in the Cry1 and Per1/2/3 genes [25]. Obesity [26] and diabetes [27,28] have been also related to a reduction in the expression of clock genes in peripheral leucocytes. Since these results raise more questions, future studies including epigenetic variability—essential in the gene regulatory process—are necessary.
There have been very few studies analyzing the gene expression of the main genes related to circadian rhythm and clock polymorphisms, with our results adding new information to the field. In particular, we showed that specific SNPs in the ARNTL (rs7924734, rs10832027, rs11022778), PER3 (rs228642, rs10127838), VRK2 (rs2678902), and MTNR1B (rs10830963) genes presented differences in their clock gene expression according to genotypes. Unexpectedly, we did not notice significant results in the CLOCK’s polymorphisms. A prior study in mice found that the rs1801260 resulted in differences in CLOCK and PER2 expression according to genotype, with the ones that had a greater expression being the carriers of the mutated allele [29]. We could not observe this difference, probably due to our sample size and differences between humans and mice. Changes in the expression of circadian genes may have an implication on physiological processes, as the molecular clock core regulates hormonal, sleep/wake, food intake, or blood pressure rhythms. Thus, clock gene expression has been related to different components of cardiometabolic risk factors. Vieira et al. detected alterations in clock gene expression in visceral adipose from obese subjects and a correlation between CLOCK gene and large-density lipoprotein (LDL) levels [30]. Similar results were found by Gómez-Abellán in morbidly obese women. PER2 expression was inversely correlated to waist circumference, and ARNTL 1, PER1, and CRY1 were also negatively associated with LDL cholesterol and total cholesterol [31].
Genetic variations could explain some phenotypic consequence, such as duration of sleep, waketime, addiction, insomnia, day preferences, etc. In our work, we noticed that the polymorphisms rs238666 in the PER3 and the rs4580704 in the CLOCK gene were associated with sleep duration and the rs1051061 and rs848291 in the VRK2 with waketimes. Heritability of sleep duration has been estimated as 46% of the variability and 44% for sleep quality. However, these values vary substantially with age, being less influential at younger ages [32]. Moreover, in recent years, different GWAS studies related to sleep traits have been carried out [33,34,35]. These results highlight the role of genetic factors in explaining differences in sleep parameters among the population. Nevertheless, to validate this association, further characterization of the molecular consequences of these polymorphisms is required and evaluating clock gene expression could be one approach [36].
5. Conclusions
In short, in these Mediterranean subjects, we observed that some polymorphisms in the circadian rhythm genes are significantly associated with one or more of the main clock genes expressed in leukocytes from blood. Moreover, we detected sex-specific differences in gene expression levels of circadian genes. Changes in gene expression could be a measure in the development of metabolic disturbances and a biomarker for effective treatments of metabolic pathologies.
The authors’ responsibilities were as follows: Conceptualization, C.O.-A., R.B. and D.C.; methodology, C.O.-A., E.M.A., R.F.-C., I.G.-M. and R.B.; software, O.C.; formal analysis, O.C., D.C. and R.B.; investigation, R.B., O.P., E.M.A., C.O.-A., I.G.-M., J.V.S., L.V.V. and D.C.; resources, D.C.; data curation, R.B. and D.C.; writing—original draft preparation, C.O.-A., E.M.A., D.C. and R.B.; writing—review and editing, all the authors. All authors have read and agreed to the published version of the manuscript.
The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Human Research Ethics Committee of Valencia University, Valencia (ethical approval code H1488282121722, 6 April 2017).
Informed consent was obtained from all subjects involved in the study.
Neither the participants’ consent forms nor ethics approval included permission for open access. However, we follow a controlled data-sharing collaboration model, and data for collaborations will be available upon request pending application and approval. Investigators who are interested in this study can contact the corresponding author Dolores Corella (
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. Representation of the clock gene expression by sex. (A) Differences in CLOCK genes between males and females were not statistical significance (p = 0.189). (B) Females showed higher PER1 gene expression than men (p = 0.034). (C) A trend was observed for a higher VRK2 gene expression in females than in males (p = 0.089).
General characteristics of the study participants.
| Total (n = 102) | Male (n = 50) | Female (n = 52) | p | |
|---|---|---|---|---|
| Age (years) | 43.1 ± 13.8 | 44.0 ± 13.5 | 42.3 ± 14.2 | 0.531 |
| BMI (kg/m2) | 27.2 ± 5.1 | 29.0 ± 4.7 | 25.5 ± 4.8 | ˂0.001 |
| Waist circumference (cm) | 92.0 ± 16.2 | 100.5 ± 15.1 | 83.9 ± 12.6 | ˂0.001 |
| SBP (mm Hg) | 123.3 ± 16.4 | 129.1 ± 14.4 | 117.9 ± 16.5 | ˂0.001 |
| DBP (mm Hg) | 77.9 ± 10.5 | 80.5 ± 12.8 | 75.5 ± 7.2 | 0.020 |
| Total cholesterol (mg/dL) | 207.4 ± 40.1 | 206.0 ± 36.0 | 208.6 ± 43.9 | 0.745 |
| LDL cholesterol (mg/dL) | 136.1 ± 32.5 | 138.7 ± 30.2 | 133.6 ± 34.7 | 0.434 |
| HDL cholesterol (mg/dL) | 56.9 ± 14.1 | 49.0 ± 9.8 | 64.5 ± 13.4 | ˂0.001 |
| Triglycerides (mg/dL) | 106.7 ± 63.0 | 127.2 ± 71.6 | 87.3 ± 46.4 | ˂0.001 |
| Fasting glucose (mg/dL) | 95.1 ± 26.2 | 101.8 ± 34.8 | 88.7 ± 10.4 | 0.013 |
| Hypertension (%) | 25.7 | 38.8 | 13.5 | 0.004 |
| Obesity cases (%) | 30.4 | 42.0 | 19.2 | 0.012 |
| Type 2 diabetes (%) | 5 | 10.2 | 0.0 | 0.018 |
| Hypercholesterolemia (%) | 33 | 37.5 | 28.8 | 0.358 |
| Current smoker (%) | 16.5 | 13.3 | 19.2 | 0.435 |
| Sedentary (%) | 26.5 | 30.4 | 23.1 | 0.380 |
| Sleep duration (h) | 7.1 ± 0.9 | 7.1 ± 0.9 | 7.1 ± 0.9 | 0.791 |
| Waketime (hour:min) | 7:45 ± 1:02 | 7:48 ± 1:18 | 7:42 ± 0:42 | 0.573 |
Values are mean ± SE (standar error) for continuous variables and % for categorical variables; p, p-value for the comparisons (means or %) between men and women; SBP, systolic blood pressure; DBP, diastolic blood pressure; LDL, large-density lipoprotein; HDL, high-density lipoprotein obesity, BMI ≥ 30 kg/m2; type 2 diabetes, antidiabetic drug or glucose ≥ 126 mg/dL; hypercholesterolemia, hypolipidemic drugs or LDL cholesterol ≥ 160 mg/dL; hypertension, antihypertensive drug or SBP ≥ 140 mmHg or DBP ≥ 90 mmHg; Waketime was represented as hours:minute. BMI, body mass index.
Relationship between clock gene polymorphisms and clock gene expression.
| CLOCK Gene Expression | PER1 Gene Expression | VRK2 Gene Expression | |||||
|---|---|---|---|---|---|---|---|
| SNP | Gene | β ± SE | p-value | β ± SE | p-value | β ± SE | p-value |
| rs4580704 | CLOCK | 0.296 ± 0.462 | 0.524 | 0.153 ± 0.468 | 0.745 | 0.221 ± 0.426 | 0.606 |
| rs1056547 | CLOCK | −0.112 ±0.434 | 0.797 | −0.425 ± 0.436 | 0.333 | −0.153 ± 0.406 | 0.702 |
| rs1801260 | CLOCK | −0.203 ± 0.443 | 0.648 | −0.326 ± 0.440 | 0.461 | −0.529 ± 0.403 | 0.193 |
| rs1982350 | ARNTL | −0.711 ± 0.414 | 0.090 | −0.94 ± 0.441 | 0.831 | −0.594 ± 0.385 | 0.126 |
| rs10741616 | ARNTL | −0.600 ± 0.442 | 0.179 | 0.316 ± 0.462 | 0.496 | −0.303 ± 0.417 | 0.471 |
| rs7924734 | ARNTL | 1.092 ± 0.484 | 0.027 | 0.259 ± 0.523 | 0.622 | 0.881 ± 0.455 | 0.056 |
| rs6486121 | ARNTL | 0.513 ± 0.400 | 0.204 | −0.353 ± 0.403 | 0.385 | 0.055 ± 0.375 | 0.883 |
| rs10832027 | ARNTL | 0.959 ± 0.469 | 0.044 | 0.302 ± 0.506 | 0.552 | 1.043 ± 0.435 | 0.019 |
| rs7937060 | ARNTL | 0.561 ± 0.390 | 0.155 | −0.222 ± 0.396 | 0.557 | 0.073 ± 0.373 | 0.845 |
| rs11022778 | ARNTL | −0.756 ± 0.442 | 0.091 | −0.623 ± 0.446 | 0.166 | −1.193 ± 0.401 | 0.004 |
| rs969485 | ARNTL | −0.497 ± 0.447 | 0.270 | −0.384 ± 0.447 | 0.393 | −0.708 ± 0.411 | 0.089 |
| rs2518023 | PER1 | 0.749 ± 1.066 | 0.484 | 0.426 ± 1.085 | 0.696 | 0.741 ± 1.000 | 0.461 |
| rs228729 | PER3 | 0.196 ± 0.459 | 0.670 | −0.149 ± 0.458 | 0.746 | −0.182 ± 0.427 | 0.671 |
| rs707467 | PER3 | 0.617 ± 0.631 | 0.332 | −0.436 ± 0.662 | 0.513 | 0.522 ± 0.596 | 0.384 |
| rs228642 | PER3 | −0.392 ± 0.496 | 0.431 | −1.180 ± 0.487 | 0.018 | −0.364 ± 0.468 | 0.439 |
| rs228666 | PER3 | −0.255 ± 0.406 | 0.532 | 0.152 ± 0.410 | 0.711 | −0.492 ± 0.384 | 0.205 |
| rs10127838 | PER3 | 0.568 ± 0.427 | 0.188 | 0.922 ± 0.429 | 0.035 | 0.731 ± 0.398 | 0.071 |
| rs2678902 | VRK2 | 0.987 ± 0.434 | 0.026 | 0.728 ± 0.453 | 0.112 | 0.607 ± 0.417 | 0.150 |
| rs2678880 | VRK2 | −0.093 ± 0.409 | 0.821 | 0.040 ± 0.415 | 0.923 | −0.219 ± 0.387 | 0.574 |
| rs7596775 | VRK2 | 0.382 ± 0.410 | 0.355 | −0.331 ± 0.417 | 0.430 | −0.144 ± 0.399 | 0.718 |
| rs2312147 | VRK2 | 0.109 ± 0.434 | 0.803 | 0.158 ± 0.442 | 0.721 | 0.235 ± 0.413 | 0.570 |
| rs7561688 | VRK2 | 0.749 ± 0.436 | 0.090 * | 0.542 ± 0.449 | 0.231 | 0.427 ± 0.418 | 0.310 |
| rs986512 | VRK2 | −0.346 ± 0.422 | 0.414 | −0.019 ± 0.434 | 0.964 | −0.403 ± 0.398 | 0.314 |
| rs1051061 | VRK2 | −0.060 ± 0.462 | 0.897 | 0.057 ± 0.470 | 0.904 | −0.228 ± 0.437 | 0.604 |
| rs848292 | VRK2 | 0.429 ± 0.462 | 0.356 | −0.057 ± 0.474 | 0.904 | 0.391 ± 0.440 | 0.377 |
| rs848291 | VRK2 | 0.135 ± 0.453 | 0.766 | 0.002 ± 0.461 | 0.997 | 0.314 ± 0.427 | 0.465 |
| rs10830963 | MNTR1B | −0.400 ± 0.485 | 0.412 | −0.586 ± 0.489 | 0.235 | −1.117 ± 0.430 | 0.011 |
| Specific PRS | PRS | 0.516 ± 0.135 | <0.001 | 0.867 ± 0.248 | 0.001 | 0.508 ± 0.127 | <0.001 |
Values are presented as β ± SE (standard error); SNP, single nucleotide polymorphism; PRS, polygenic risk score; p-values were obtained in multivariable linear regressions models adjusted for age, sex, and kcal.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Rijo-Ferreira, F.; Takahashi, J.S. Genomics of circadian rhythms in health and disease. Genome Med.; 2019; 11, 82. [DOI: https://dx.doi.org/10.1186/s13073-019-0704-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31847894]
2. Stenvers, D.J.; Scheer, F.A.J.L.; Schrauwen, P.; la Fleur, S.E.; Kalsbeek, A. Circadian clocks and insulin resistance. Nat. Rev. Endocrinol.; 2019; 15, pp. 75-89. [DOI: https://dx.doi.org/10.1038/s41574-018-0122-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30531917]
3. 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. [DOI: https://dx.doi.org/10.1186/s12933-015-0327-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26739996]
4. Froy, O.; Garaulet, M. The circadian clock in white and brown adipose tissue: Mechanistic, endocrine, and clinical aspects. Endocr. Rev.; 2018; 39, pp. 261-273. [DOI: https://dx.doi.org/10.1210/er.2017-00193]
5. 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. [DOI: https://dx.doi.org/10.1155/2016/2746909]
6. Jones, S.E.; Tyrrell, J.; Wood, A.R.; Beaumont, R.N.; Ruth, K.S.; Tuke, M.A.; Yaghootkar, H.; Hu, Y.; Teder-Laving, M.; Hayward, C. et al. Genome-wide association analyses in 128,266 individuals identifies new morningness and sleep duration loci. PLoS Genet.; 2016; 12, e1006125. [DOI: https://dx.doi.org/10.1371/journal.pgen.1006125]
7. Veatch, O.J.; Keenan, B.T.; Gehrman, P.R.; Malow, B.A.; Pack, A.I. Pleiotropic genetic effects influencing sleep and neurological disorders. Lancet Neurol.; 2017; 16, pp. 158-170. [DOI: https://dx.doi.org/10.1016/S1474-4422(16)30339-8]
8. St-Onge, M.P.; Grandner, M.A.; Brown, D.; Conroy, M.B.; Jean-Louis, G.; Coons, M.; Bhatt, D.L. Sleep duration and quality: Impact on lifestyle behaviors and cardiometabolic health: A scientific statement from the American heart association. Circulation; 2016; 134, pp. e367-e386. [DOI: https://dx.doi.org/10.1161/CIR.0000000000000444]
9. Yamazaki, E.M.; Goel, N. Genetics of circadian and sleep measures in adults: Implications for sleep medicine. Curr. Sleep Med. Rep.; 2020; 6, pp. 32-45. [DOI: https://dx.doi.org/10.1007/s40675-020-00165-z]
10. Ulloa-Martínez, M.; Burguete-García, A.I.; Murugesan, S.; Hoyo-Vadillo, C.; Cruz-Lopez, M.; García-Mena, J. Expression of candidate genes associated with obesity in peripheral white blood cells of Mexican children. Arch. Med. Sci.; 2016; 12, pp. 969-976. [DOI: https://dx.doi.org/10.5114/aoms.2016.58126]
11. Leder, L.; Kolehmainen, M.; Narverud, I.; Dahlman, I.; Myhrstad, M.C.W.; de Mello, V.D.; Paananen, J.; Carlberg, C.; Schwab, U.; Herzig, K.H. et al. Effects of a healthy Nordic diet on gene expression changes in peripheral blood mononuclear cells in response to an oral glucose tolerance test in subjects with metabolic syndrome: A SYSDIET sub-study. Genes Nutr.; 2016; 11, 3. [DOI: https://dx.doi.org/10.1186/s12263-016-0521-4]
12. Sorlí, J.V.; Barragán, R.; Coltell, O.; Portolés, O.; Pascual, E.C.; Ortega-Azorín, C.; González, J.I.; Estruch, R.; Saiz, C.; Pérez-Fidalgo, A. et al. Chronological age interacts with the circadian melatonin receptor 1B gene variation, determining fasting glucose concentrations in mediterranean populations. additional analyses on Type-2 Diabetes risk. Nutrients; 2020; 12, 3323. [DOI: https://dx.doi.org/10.3390/nu12113323]
13. Barragán, R.; Coltell, O.; Portolés, O.; Asensio, E.M.; Sorlí, J.V.; Ortega-Azorín, C.; González, J.I.; Sáiz, C.; Fernández-Carrión, R.; Ordovas, J.M. et al. Bitter, sweet, salty, sour and umami taste perception decreases with age: Sex-specific analysis, modulation by genetic variants and taste-preference associations in 18 to 80 year-old subjects. Nutrients; 2018; 10, 1536. [DOI: https://dx.doi.org/10.3390/nu10101539]
14. Coltell, O.; Ortega-Azorín, C.; Sorlí, J.V.; Portolés, O.; Asensio, E.M.; Saiz, C.; Barragán, R.; Estruch, R.; Corella, D. Circulating adiponectin and its association with metabolic traits and Type 2 Diabetes: Gene-diet interactions focusing on selected gene variants and at the genome-wide level in high-cardiovascular risk Mediterranean subjects. Nutrients; 2021; 13, 541. [DOI: https://dx.doi.org/10.3390/nu13020541]
15. Chomczynski, P.; Sacchi, N. Single-Step Method of RNA Isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal. Biochem.; 1987; 162, pp. 156-159. [DOI: https://dx.doi.org/10.1016/0003-2697(87)90021-2]
16. Schmittgen, T.D.; Livak, K.J. Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc.; 2008; 3, pp. 1101-1108. [DOI: https://dx.doi.org/10.1038/nprot.2008.73]
17. Aguet, F.; Barbeira, A.N.; Bonazzola, R.; Brown, A.; Castel, S.E.; Jo, B.; Kasela, S.; Kim-Hellmuth, S.; Liang, Y.; Oliva, M. et al. The impact of sex on gene expression across human tissues. Science; 2020; 369, eaba3066. [DOI: https://dx.doi.org/10.1126/SCIENCE.ABA3066]
18. Lopes-Ramos, C.M.; Chen, C.Y.; Kuijjer, M.L.; Paulson, J.N.; Sonawane, A.R.; Fagny, M.; Platig, J.; Glass, K.; Quackenbush, J.; DeMeo, D.L. Sex differences in gene expression and regulatory networks across 29 human tissues. Cell Rep.; 2020; 31, 107795. [DOI: https://dx.doi.org/10.1016/j.celrep.2020.107795]
19. Gómez-Abellán, P.; Madrid, J.A.; Luján, J.A.; Frutos, M.D.; González, R.; Martínez-Augustín, O.; de Medina, F.S.; Ordovás, J.M.; Garaulet, M. Sexual dimorphism in clock genes expression in human adipose tissue. Obes. Surg.; 2012; 22, pp. 105-112. [DOI: https://dx.doi.org/10.1007/s11695-011-0539-2]
20. Lim, A.S.P.; Myers, A.J.; Yu, L.; Buchman, A.S.; Duffy, J.F.; de Jager, P.L.; Bennett, D.A. Sex difference in daily rhythms of clock gene expression in the aged human cerebral cortex. J. Biol. Rhythms; 2013; 28, pp. 117-129. [DOI: https://dx.doi.org/10.1177/0748730413478552]
21. Škrlec, I.; Talapko, J.; Juzbašić, M.; Steiner, R. Sex differences in circadian clock genes and myocardial infarction susceptibility. J. Cardiovasc. Dev. Dis.; 2021; 8, 53. [DOI: https://dx.doi.org/10.3390/jcdd8050053]
22. Hatcher, K.M.; Royston, S.E.; Mahoney, M.M. Modulation of circadian rhythms through estrogen receptor signaling. Eur. J. Neurosci.; 2020; 51, pp. 217-228. [DOI: https://dx.doi.org/10.1111/ejn.14184]
23. Emeklİ, R.; İsmaİloğullari, S.; Bayram, A.; Akalin, H.; Tuncel, G.; Dündar, M. Comparing expression levels of PERIOD genes PER1, PER2 and PER3 in chronic insomnia patients and medical staff working in the night shift. Sleep Med.; 2020; 73, pp. 101-105. [DOI: https://dx.doi.org/10.1016/j.sleep.2020.04.027]
24. Möller-Levet, C.S.; Archer, S.N.; Bucca, G.; Laing, E.E.; Slak, A.; Kabiljo, R.; Lo, J.C.; Santhi, N.; von Schantz, M.; Smith, C.P. et al. Effects of insufficient sleep on circadian rhythmicity and expression amplitude of the human blood transcriptome. Proc. Natl. Acad. Sci. USA.; 2013; 19, pp. E1132-E1141. [DOI: https://dx.doi.org/10.1073/pnas.1217154110]
25. Li, W.; Wang, Z.; Cao, J.; Dong, Y.; Chen, Y. Role of sleep restriction in daily rhythms of expression of hypothalamic core clock genes in mice. Curr. Issues Mol. Biol.; 2022; 25, pp. 609-625. [DOI: https://dx.doi.org/10.3390/cimb44020042]
26. Yamaoka, M.; Maeda, N.; Nakamura, S.; Kashine, S.; Nakagawa, Y.; Hiuge-Shimizu, A.; Okita, K.; Imagawa, A.; Matsuzawa, Y.; Matsubara, K. et al. A pilot investigation of visceral fat adiposity and gene expression profile in peripheral blood cells. PLoS ONE; 2012; 7, e47377. [DOI: https://dx.doi.org/10.1371/journal.pone.0047377]
27. Ando, H.; Takamura, T.; Matsuzawa-Nagata, N.; Shima, K.R.; Eto, T.; Misu, H.; Shiramoto, M.; Tsuru, T.; Irie, S.; Fujimura, A. et al. Clock gene expression in peripheral leucocytes of patients with Type 2 Diabetes. Diabetologia; 2009; 52, pp. 329-335. [DOI: https://dx.doi.org/10.1007/s00125-008-1194-6]
28. Yu, R.; Tian, L.; Ding, Y.; Gao, Y.; Li, D.; Tang, Y. Correlation between inflammatory markers and impaired circadian clock gene expression in type 2 diabetes mellitus. Diabetes Res. Clin. Pract.; 2019; 156, 107831. [DOI: https://dx.doi.org/10.1016/j.diabres.2019.107831]
29. Ozburn, A.R.; Purohit, K.; Parekh, P.K.; Kaplan, G.N.; Falcon, E.; Mukherjee, S.; Cates, H.M.; McClung, C.A. Functional implications of the CLOCK 3111T/C single-nucleotide polymorphism. Front. Psychiatry; 2016; 7, 67. [DOI: https://dx.doi.org/10.3389/fpsyt.2016.00067]
30. Vieira, E.; Ruano, E.G.; Figueroa, A.L.C.; Aranda, G.; Momblan, D.; Carmona, F.; Gomis, R.; Vidal, J.; Hanzu, F.A. Altered clock gene expression in obese visceral adipose tissue is associated with metabolic syndrome. PLoS ONE; 2014; 9, e0111678. [DOI: https://dx.doi.org/10.1371/journal.pone.0111678]
31. Gómez-Abellán, P.; Hernández-Morante, J.J.; Luján, J.A.; Madrid, J.A.; Garaulet, M. Clock genes are implicated in the human metabolic syndrome. Int. J. Obes.; 2008; 32, pp. 121-128. [DOI: https://dx.doi.org/10.1038/sj.ijo.0803689] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17653067]
32. Kocevska, D.; Barclay, N.L.; Bramer, W.M.; Gehrman, P.R.; van Someren, E.J.W. Heritability of sleep duration and quality: A systematic review and meta-analysis. Sleep Med. Rev.; 2021; 59, 101448. [DOI: https://dx.doi.org/10.1016/j.smrv.2021.101448] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33636423]
33. Dashti, H.S.; Jones, S.E.; Wood, A.R.; Lane, J.M.; van Hees, V.T.; Wang, H.; Rhodes, J.A.; Song, Y.; Patel, K.; Anderson, S.G. et al. Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates. Nat. Commun.; 2019; 10, 1100. [DOI: https://dx.doi.org/10.1038/s41467-019-08917-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30846698]
34. Jones, S.E.; Lane, J.M.; Wood, A.R.; van Hees, V.T.; Tyrrell, J.; Beaumont, R.N.; Jeffries, A.R.; Dashti, H.S.; Hillsdon, M.; Ruth, K.S. et al. Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nat. Commun.; 2019; 10, 343. [DOI: https://dx.doi.org/10.1038/s41467-018-08259-7]
35. Dashti, H.S.; Daghlas, I.; Lane, J.M.; Huang, Y.; Udler, M.S.; Wang, H.; Ollila, H.M.; Jones, S.E.; Kim, J.; Wood, A.R. et al. Genetic determinants of daytime napping and effects on cardiometabolic health. Nat. Commun.; 2021; 12, 900. [DOI: https://dx.doi.org/10.1038/s41467-020-20585-3]
36. Zhang, L.; Ptáček, L.J.; Fu, Y.H. Diversity of human clock genotypes and consequences. Prog. Mol. Biol. Transl. Sci.; 2013; 119, pp. 51-81. [DOI: https://dx.doi.org/10.1016/B978-0-12-396971-2.00003-8]
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Sorlí, José V 1
; Coltell, Oscar 2
; Gonzalez-Monje, Inmaculada 3 ; Fernández-Carrión, Rebeca 1 ; Villamil, Laura V 4 ; Portolés, Olga 1 ; Corella, Dolores 1
; Ortega-Azorín, Carolina 5
; Asensio, Eva M 1
1 Department of Preventive Medicine and Public Health, School of Medicine, University of Valencia, 46010 Valencia, Spain; CIBER Fisiopatología de la Obesidad y Nutrición, Instituto de Salud Carlos III, 28029 Madrid, Spain
2 CIBER Fisiopatología de la Obesidad y Nutrición, Instituto de Salud Carlos III, 28029 Madrid, Spain; Department of Computer Languages and Systems, Universitat Jaume I, 12071 Castellón, Spain
3 Department of Preventive Medicine and Public Health, School of Medicine, University of Valencia, 46010 Valencia, Spain
4 Department of Physiology, School of Medicine, University Antonio Nariño, 111511 Bogotá, Colombia
5 Department of Preventive Medicine and Public Health, School of Medicine, University of Valencia, 46010 Valencia, Spain; CIBER Fisiopatología de la Obesidad y Nutrición, Instituto de Salud Carlos III, 28029 Madrid, Spain; Department of Computer Languages and Systems, Universitat Jaume I, 12071 Castellón, Spain