OBJECTIVE

Common circadian-related gene variants associate with increased risk for metabolic alterations including type 2 diabetes. However, little is known about whether diet and sleep could modify associations between circadian-related variants (CLOCK-rs1801260, CRY2-rs11605924, MTNR1B-rs1387153, MTNR1B-rs10830963, NR1D1-rs2314339) and cardiometabolic traits (fasting glucose [FG], HOMA-insulin resistance, BMI, waist circumference, and HDL-cholesterol) to facilitate personalized recommendations.

RESEARCH DESIGN AND METHODS

We conducted inverse-variance weighted, fixed-effect meta-analyses of results of adjusted associations and interactions between dietary intake/sleep duration and selected variants on cardiometabolic traits from 15 cohort studies including up to 28,190 participants of European descent from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium.

RESULTS

We observed significant associations between relative macronutrient intakes and glycemic traits and short sleep duration (<7 h) and higher FG and replicated known MTNR1B associations with glycemic traits. No interactions were evident after accounting for multiple comparisons. However, we observed nominally significant interactions (all P < 0.01) between carbohydrate intake and MTNR1B-rs1387153 for FG with a 0.003 mmol/L higher FG with each additional 1% carbohydrate intake in the presence of the T allele, between sleep duration and CRY2-rs11605924 for HDL-cholesterol with a 0.010 mmol/L higher HDL-cholesterol with each additional hour of sleep in the presence of the A allele, and between long sleep duration (≥9 h) and MTNR1B-rs1387153 for BMI with a 0.60 kg/m2 higher BMI with long sleep duration in the presence of the T allele relative to normal sleep duration (≥7 to <9 h).

CONCLUSIONS

Our results suggest that lower carbohydrate intake and normal sleep duration may ameliorate cardiometabolic abnormalities conferred by common circadian-related genetic variants. Until further mechanistic examination of the nominally significant interactions is conducted, recommendations applicable to the general population regarding diet—specifically higher carbohydrate and lower fat composition—and normal sleep duration should continue to be emphasized among individuals with the investigated circadian-related gene variants.

Gene-environment interactions can identify potential opportunities for personalized health interventions for individuals who are genetically susceptible to type 2 diabetes and related chronic diseases (1). During the past decade, researchers have examined lifestyle interventions, particularly those related to diet, physical activity, and sleep, for individuals at increased genetic risk for metabolic alterations such as type 2 diabetes (2,3). For example, dietary changes in carbohydrate (CHO) and fat intake have been shown to attenuate a genetic predisposition to elevated fasting glucose (FG) (4), insulin resistance (5), and type 2 diabetes (6). Sleep duration has also been assessed as a modifying factor in the associations between genetics and type 2 diabetes because of the association of sleep with diet and chronic disease and its effect on the genetic risk for obesity (7). Identifying optimal and personalized therapies for the primary prevention of type 2 diabetes through gene-environment investigations is critical to public health for a disease that affects an estimated 29.1 million Americans (9.3% of the U.S. population) (1,8). In addition, because only 10% of the total heritability of type 2 diabetes is accounted for by genetic variants (9), gene-environment investigations may also reveal novel biological pathways and genetic loci pertinent to this disease.

Glucose homeostasis and insulin secretion are among several biological processes that are controlled by the circadian biological clock, which is maintained endogenously through a transcription-translation feedback loop composed of clock genes (10). Glycemic control is mediated through multiple processes, including circadian regulation of hepatic glucose metabolism (11,12); secretion of adipokines, such as leptin and adiponectin (13,14); and the pancreatic secretion of insulin and glucagon (15). Experiments in clock mutant mice showing disrupted glucose homeostasis, insulin secretion and sensitivity, and other metabolic processes, along with circadian disruption in humans (16), emphasize the importance of circadian control in metabolic control (15,17).

In support of the link between the circadian system and glycemic control in humans are results from genome-wide association studies (GWAS) of FG (18,19) and type 2 diabetes (9,20) that have reported associations with clock gene CRY2, encoding cryptochrome 2, and the circadian-related melatonin receptor 1B gene MTNR1B. In addition, the circadian locomotor output of clock genes cycles kaput CLOCK and, more recently, the nuclear receptor rev-erb-α NR1D1 have been associated with related metabolic traits, including lower circulating concentrations of HDL-cholesterol (HDL-C) and elevated central adiposity (2123). Because metabolic traits are important predictors of type 2 diabetes, these loci may also be relevant to the pathogenesis of type 2 diabetes (24). Thus, investigating whether lifestyle modifications—particularly diet, for its potent role in entraining circadian clocks in metabolic tissues (25), and sleep, for its putative effect on disease risk (7)—attenuate circadian-related genetic predispositions to metabolic disruptions may facilitate the development of personalized recommendations to improve type 2 diabetes prevention strategies.

In cross-sectional meta-analyses of large population-based cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium, we tested whether dietary intake (total CHO, total fat, polyunsaturated fatty acid [PUFA], monounsaturated fatty acid [MUFA], and saturated fatty acid [SFA]) and sleep duration (continuous and categorical) modify the associations between five common circadian-related gene variants (CLOCK-rs1801260, CRY2-rs11605924, MTNR1B-rs1387153, MTNR1B-rs10830963, and NR1D1-rs2314339) and the two glycemic traits of FG and HOMA-insulin resistance (HOMA-IR), as well as related anthropometric (BMI and waist circumference) and lipid (HDL-C) traits. These outcomes are related to cardiometabolic disease and have previously been shown to associate with the selected genetic variants.

Cohorts

The present cross-sectional meta-analyses include up to 28,190 participants of European descent from the following 15 cohort studies of the CHARGE Consortium Nutrition Working Group (Supplementary Table 1): Coronary Artery Risk Development in Young Adults Study (CARDIA); Corogene Controls; Cardiovascular Health Study (CHS); Dietary, Lifestyle, and Genetic determinants of Obesity and Metabolic syndrome (DILGOM); Family Heart Study (FamHS); Framingham Offspring Study (FOS); Genetics of Lipid Lowering Drugs and Diet Network (GOLDN); GOYA MALE; Helsinki Birth Cohort Study (HBCS); Invecchiare in Chianti (aging in the Chianti area, InCHIANTI); Inter99; Multi-Ethnic Study of Atherosclerosis (MESA); The Rotterdam Study; The Hellenic Study of Interactions between SNPs and Eating in Atherosclerosis Susceptibility (THISEAS); and the Cardiovascular Risk in Young Finns Study (YFS). Participants provided written informed consent, and the study protocols were approved by local institutional review boards and/or oversight committees.

Dietary Assessment

Dietary data were collected via validated food-frequency questionnaires (13 cohorts), dietary recall (1 cohort), and food record (1 cohort) (Supplementary Table 2). The type of food-frequency questionnaire used in each cohort differed slightly to capture the dietary habits of the population of interest. The present analysis quantified total CHO intake and total fat intake as percentages of total energy intake. Additional analyses used percentage of energy from specific fatty acids, including PUFA, MUFA, and SFA. Total energy intake from protein was not included in the present analysis because of the lack of evidence supporting protein intake in gene-environment interactions (3).

Sleep Assessment

Data on habitual weekday/workday nighttime sleep duration in hours per night were obtained from self-reported responses to questions such as, “How many hours of sleep do you usually get at night?” or were calculated from self-reported weekday/workday bed and rise times (Supplementary Table 3). Responses were analyzed as continuous and categorical variables. Commonly accepted cutoffs were used to create three sleep duration categories: short (<7 h), normal (≥7 to <9 h), and long (≥9 h).

Outcome Measurements

Cohort-specific assessment methods for FG (mmol/L), BMI (kg/m2), waist circumference (cm), and HDL-C (mmol/L) are described in detail in Supplementary Table 3. HOMA-IR was estimated from fasting insulin and FG concentrations, using the previously validated equation [HOMA-IR = FG (mmol/L) × fasting insulin (mU/L)/22.5], and was natural log-transformed to reduce skew before data analysis.

Genotyping

We selected five single nucleotide polymorphisms (SNPs) in circadian-related genes based on previous reports from GWAS meta-analysis (CRY2-rs11605924, MTNR1B-rs10830963), replicated candidate gene association studies (NR1D1-rs2314339), gene-environment interaction studies, or a combination of these findings (CLOCK-rs1801260, MTNR1B-rs1387153) that showed associations with type 2 diabetes, FG, and/or BMI (4,9,1821). SNPs and/or SNPs in linkage disequilibrium (r2 > 0.80; HapMap III release 2 data set) were previously directly genotyped or imputed by participating cohorts before inclusion in this analysis (Supplementary Table 4). SNPs were assessed for quality control: genotyped SNPs were excluded on the basis of low call rate (<95%) and departure from Hardy-Weinberg equilibrium (<1E-06), and imputed SNPs were removed on the basis of low imputation quality (MACH: R2 < 0.3 or IMPUTE: proper info <0.4). Not all SNPs were available in all participating cohorts (Supplementary Table 5), and therefore, total sample sizes for analyses varied.

Cohort-Specific Analyses

All participating cohort-specific statistical analyses followed a uniform analysis plan. First, main associations between dietary intake or sleep duration and all outcomes were estimated with adjustment for age, sex, BMI (except for BMI outcome), and study site (in CARDIA, CHS, FamHS, GOLDN, InCHIANTI, and MESA) using linear regression models. Second, main associations between selected SNPs and all outcomes were investigated by using linear regression models and an additive genetic model adjusted for the aforementioned covariates in addition to family or population structure (in Corogene Controls, DILGOM, FamHS, FOS, GOLDN, MESA, The Rotterdam Study, and YFS), and genotype batch (in FamHS). Third, for our primary analysis of interest, 175 interactions (7 environmental variables × 5 SNPs × 5 outcomes) between dietary intake or sleep duration (continuous and categorical) and the selected SNPs on all outcomes were tested by using intake/sleep duration × SNP cross-product terms and including main-effect terms in linear regression models adjusted for the aforementioned covariates. Participants within each cohort were excluded from the analysis if they were shift workers, on sleep medications or antidepressant medications, reported bedtimes after 5 a.m. or before 6 p.m., and/or reported sleep duration <3 or ≥16 h. For glycemic outcomes, participants with type 2 diabetes within each cohort were excluded from main association and interaction analyses.

Meta-analyses

We conducted inverse-variance weighted, fixed-effect meta-analyses using METAL (version released 2011-03-25) for 1) main associations of dietary intake/sleep duration on the outcomes, 2) main associations of the selected SNPs on the outcomes, and 3) interactions between SNPs and dietary intake/sleep duration on the outcomes.

Heterogeneity across studies was tested by using Cochran’s Q statistic and quantified using the I2 statistic. All association and interaction analyses with moderate heterogeneity (I2 > 30%) were further assessed for potential sources of heterogeneity by conducting meta-regression and sensitivity analyses. Meta-regression analyses were conducted using the R metafor package (R version 3.1.0) to assess the effect of the following moderator variables on heterogeneity of association/interaction: geographical location (U.S. vs. northern Europe vs. Mediterranean), mean age of cohort (20–64 years vs. 65–80 years), and total energy intake (<2,000 vs. ≥2,000 kcal/day). Sensitivity analyses assessed the influence of a single cohort on the meta-analyzed estimate by repeating analyses removing one cohort study at a time in the association and interaction analyses. Statistical significance for each outcome was defined at a level of 0.003, based on Bonferroni correction for 15 total independent interaction tests (3 independent environmental variables [percent energy from CHO, percent energy from total fat, and sleep duration] × 5 independent SNPs). We performed power calculations using Quanto version 1.2.4 (http://biostats.usc.edu/software) for NR1D1-rs2314339, which represents the SNP with the least possible power (based on lowest minor allele frequency and sample size). At 80% power, our sample size meets the estimated sample size to detect an interaction effect between dietary macronutrient intake (per 1% intake) and SNPs (per effect allele) on FG.

General characteristics of participants are reported in Table 1. Mean ages ranged between 32.6 and 70.2 years, and women comprised 37–62% of participants in each cohort, except for GOYA MALE (all men). Mean BMI exceeded 25 kg/m2 for all cohorts, and there were no substantial variations in waist circumference, FG, HOMA-IR, or HDL-C among cohorts. Total energy intake was generally higher in U.S. cohorts, but relative macronutrient intakes were similar across studies and differed only for MUFA and SFA intakes. The Mediterranean cohorts (InCHIANTI and THISEAS) had higher mean MUFA intake than northern European and U.S. cohorts (P < 0.0001), and The Rotterdam Study reported the highest SFA intake (P < 0.0001). In addition, the average habitual weekday/workday sleep duration was similar across cohorts, ranging between 6.8 and 8.3 h, the prevalence of short sleep duration (<7 h) ranged between 18.1 and 41.6%, and long sleep duration (≥9 h) ranged between 4.0 and 11.9%.

Table 1

General characteristics of participating CHARGE cohorts*

CARDIACorogene ControlsCHSDILGOMFamHSFOSGOLDNGOYA MALEHBCSInCHIANTIInter99MESARotterdamTHISEASYFS
Characteristicsn = 1,073n = 549n = 2,918n = 3,136n = 3,326n = 917n = 810n = 187n = 1,184n = 922n = 6,000n = 2,028n = 2,367n = 931n = 1,730
Age (years) 32.6 ± 3.3 51.2 ± 13.4 71.2 ± 4 49.9 ± 13.3 51.4 ± 13.3 57.5 ± 9 48.4 ± 16.3 44.7 ± 6.5 69 ± 2.9 64.5 ± 14.2 46.2 ± 7.9 61.7 ± 9.7 70.2 ± 5 57 ± 12.4 37.8 ± 5 
Sex (% women) 52.4 54.8 61.9 56.6 53.3 55.4 52.0 60.6 53.5 51.4 50.96 58.0 37.1 55.4 
BMI (kg/m225.7 ± 5.2 26.5 ± 4.5 26.4 ± 4.5 26.7 ± 4.6 27.7 ± 5.5 28 ± 5.2 28.3 ± 5.6 25.6 ± 3.2 27.5 ± 4.6 27.3 ± 4.2 26.3 ± 4.6 27.8 ± 5.1 26.8 ± 3.9 28.1 ± 4.3 25.9 ± 4.6 
Waist circumference (cm) 106.9 ± 36.1 89.4 ± 13.1 93.2 ± 12.9 89.5 ± 12.6 97.3 ± 15.3 97.3 ± 15.3 96.6 ± 16.6 93 ± 10.6 94.4 ± 13 91.5 ± 11.2 86.7 ± 13.4 98 ± 14.5 93.1 ± 11.5 97.3 ± 13.7 88.2 ± 13.3 
FG (mmol/L) 5.0 ± 0.7 5.8 ± 0.5 6 ± 1.6 5.7 ± 0.5 5.5 ± 1.4 5.5 ± 0.8 5.6 ± 1 5.7 ± 0.6 5.7 ± 1.2 5.2 ± 1.5 5.6 ± 1.1 5.1 ± 1.2 5.6 ± 0.9 5.9 ± 1.8 5.3 ± 0.9 
Fasting insulin (pmol/L) 30.4 ± 29.5 42.9 ± 26.2 108.1 ± 125 43.9 ± 29.7 79.4 ± 94.7 NA 95.4 ± 56.7 37.4 ± 26.7 70.9 ± 56 77.2 ± 42.8 42.9 ± 29.6 63.2 ± 38.2 77.9 ± 54.7 NA 61.4 ± 60.1 
HOMA-IR 2.9 ± 4.2 1.6 ± 1.1 4.5 ± 9.2 1.7 ± 1.3 3.1 ± 6.1 NA 3.5 ± 2.5 1.4 ± 1.1 2.6 ± 2.3 2.6 ± 1.9 1.9 ± 1.7 1.4 ± 1.3 2.9 ± 2.5 NA 2.2 ± 3.3 
HDL-C (mmol/L) 1.3 ± 0.4 1.5 ± 0.3 1.4 ± 0.4 1.5 ± 0.3 1.3 ± 0.4 1.3 ± 0.4 1.2 ± 0.3 1.3 ± 0.3 1.6 ± 0.4 1.5 ± 0.4 1.4 ± 0.4 1.4 ± 0.4 1.4 ± 0.4 1.3 ± 0.4 1.3 ± 0.3 
Sleep duration (h) NA 7.2 ± 1 7.2 ± 1.3 7.3 ± 1 NA 7.1 ± 1.1 NA NA 7.4 ± 1.1 6.8 ± 1.4 7.3 ± 0.9 7 ± 1.2 NA 7.1 ± 1.5 7.4 ± 0.8 
 Short <7 h, n (%) NA 99 (18.1) 389 (26.8) 568 (18.2) NA 231 (25.7) NA NA 492 (41.6) NA 851 (20.0) 586 (29.3) NA 198 (36.1) 330 (19.1) 
 Normal ≥7 to <9 h, n (%) NA 406 (74.4) 892 (61.3) 2324 (74.3) NA 613 (68.2) NA NA 651 (55.0) NA 3174 (74.7) 1281 (52.4) NA 287 (52.4) 1330 (76.9) 
 Long ≥9 h, n (%) NA 41 (7.5) 173 (11.9) 234 (7.5) NA 55 (6.1) NA NA 41 (5.3) NA 225 (5.3) 130 (11.5) NA 63 (11.5) 70 (4.0) 
Dietary intake                
 Total energy (kcal/day) 2,587 ± 928 2,302 ± 776 1,999 ± 613 2,538 ± 903 1,755 ± 619 1,878 ± 616 2,115 ± 1,181 2,490 ± 558 2,212 ± 807 2,066 ± 601 2,332 ± 870 1,602 ± 705 2,007 ± 520 1,922 ± 1,046 2,383 ± 761 
 Total fats (% total energy) 34 ± 6.7 31.2 ± 4.9 32.3 ± 6 31.3 ± 4.8 30.8 ± 7.4 29.5 ± 6.4 35.4 ± 6.9 32.4 ± 4.6 32.5 ± 5.2 31 ± 5.2 32.6 ± 7 32.1 ± 7.3 36.0 ± 6 36.1 ± 6 32.9 ± 4.8 
 PUFA (% total energy) 6.9 ± 2 5.5 ± 1.2 7.5 ± 2.2 5.4 ± 1.1 4.5 ± 1.4 5.7 ± 1.6 6.9 ± 2.2 4.1 ± 0.7 5.2 ± 1.2 3.4 ± 0.7 5.3 ± 1.5 6.3 ± 1.8 6.9 ± 2.7 5 ± 1.4 5.3 ± 1.1 
 MUFA (% total energy) 12.5 ± 2.8 11.1 ± 2.2 11.6 ± 2.5 11.1 ± 2.1 12 ± 3.2 11 ± 2.6 13.3 ± 2.8 9.8 ± 1.7 10.8 ± 2.1 15.6 ± 3.2 10.6 ± 2.7 12.3 ± 3 12.3 ± 2.6 16.2 ± 3.6 11.1 ± 2 
 SFA (% total energy) 11.9 ± 3.2 11.3 ± 2.4 10.3 ± 2.2 11.8 ± 2.5 11.3 ± 3.2 10.1 ± 2.8 11.8 ± 2.8 13.5 ± 2.4 11.8 ± 2.5 10.4 ± 2.2 12.3 ± 7.9 10.8 ± 3.4 14.1 ± 3 11.8 ± 3.2 11.8 ± 2.4 
 CHO (% total energy) 50.2 ± 7.8 48.1 ± 6.2 52.2 ± 7.9 49 ± 6.1 50.8 ± 10.1 51.6 ± 8.2 48.9 ± 8.7 41.8 ± 6.6 47 ± 6.5 51.4 ± 6.7 47.5 ± 7.9 50.7 ± 9 43.3 ± 6.9 44.5 ± 8.7 45.9 ± 5.8 
 Proteins (% total energy) 15 ± 2.6 17.6 ± 2.6 19.1 ± 3.2 17.7 ± 2.5 18.3 ± 4.1 17.1 ± 3.1 15.8 ± 2.9 15.2 ± 2.3 17.1 ± 2.5 15.7 ± 2.1 15.1 ± 2.6 15.7 ± 3.3 16.9 ± 3 17.9 ± 2.4 17.5 ± 2.4 
CARDIACorogene ControlsCHSDILGOMFamHSFOSGOLDNGOYA MALEHBCSInCHIANTIInter99MESARotterdamTHISEASYFS
Characteristicsn = 1,073n = 549n = 2,918n = 3,136n = 3,326n = 917n = 810n = 187n = 1,184n = 922n = 6,000n = 2,028n = 2,367n = 931n = 1,730
Age (years) 32.6 ± 3.3 51.2 ± 13.4 71.2 ± 4 49.9 ± 13.3 51.4 ± 13.3 57.5 ± 9 48.4 ± 16.3 44.7 ± 6.5 69 ± 2.9 64.5 ± 14.2 46.2 ± 7.9 61.7 ± 9.7 70.2 ± 5 57 ± 12.4 37.8 ± 5 
Sex (% women) 52.4 54.8 61.9 56.6 53.3 55.4 52.0 60.6 53.5 51.4 50.96 58.0 37.1 55.4 
BMI (kg/m225.7 ± 5.2 26.5 ± 4.5 26.4 ± 4.5 26.7 ± 4.6 27.7 ± 5.5 28 ± 5.2 28.3 ± 5.6 25.6 ± 3.2 27.5 ± 4.6 27.3 ± 4.2 26.3 ± 4.6 27.8 ± 5.1 26.8 ± 3.9 28.1 ± 4.3 25.9 ± 4.6 
Waist circumference (cm) 106.9 ± 36.1 89.4 ± 13.1 93.2 ± 12.9 89.5 ± 12.6 97.3 ± 15.3 97.3 ± 15.3 96.6 ± 16.6 93 ± 10.6 94.4 ± 13 91.5 ± 11.2 86.7 ± 13.4 98 ± 14.5 93.1 ± 11.5 97.3 ± 13.7 88.2 ± 13.3 
FG (mmol/L) 5.0 ± 0.7 5.8 ± 0.5 6 ± 1.6 5.7 ± 0.5 5.5 ± 1.4 5.5 ± 0.8 5.6 ± 1 5.7 ± 0.6 5.7 ± 1.2 5.2 ± 1.5 5.6 ± 1.1 5.1 ± 1.2 5.6 ± 0.9 5.9 ± 1.8 5.3 ± 0.9 
Fasting insulin (pmol/L) 30.4 ± 29.5 42.9 ± 26.2 108.1 ± 125 43.9 ± 29.7 79.4 ± 94.7 NA 95.4 ± 56.7 37.4 ± 26.7 70.9 ± 56 77.2 ± 42.8 42.9 ± 29.6 63.2 ± 38.2 77.9 ± 54.7 NA 61.4 ± 60.1 
HOMA-IR 2.9 ± 4.2 1.6 ± 1.1 4.5 ± 9.2 1.7 ± 1.3 3.1 ± 6.1 NA 3.5 ± 2.5 1.4 ± 1.1 2.6 ± 2.3 2.6 ± 1.9 1.9 ± 1.7 1.4 ± 1.3 2.9 ± 2.5 NA 2.2 ± 3.3 
HDL-C (mmol/L) 1.3 ± 0.4 1.5 ± 0.3 1.4 ± 0.4 1.5 ± 0.3 1.3 ± 0.4 1.3 ± 0.4 1.2 ± 0.3 1.3 ± 0.3 1.6 ± 0.4 1.5 ± 0.4 1.4 ± 0.4 1.4 ± 0.4 1.4 ± 0.4 1.3 ± 0.4 1.3 ± 0.3 
Sleep duration (h) NA 7.2 ± 1 7.2 ± 1.3 7.3 ± 1 NA 7.1 ± 1.1 NA NA 7.4 ± 1.1 6.8 ± 1.4 7.3 ± 0.9 7 ± 1.2 NA 7.1 ± 1.5 7.4 ± 0.8 
 Short <7 h, n (%) NA 99 (18.1) 389 (26.8) 568 (18.2) NA 231 (25.7) NA NA 492 (41.6) NA 851 (20.0) 586 (29.3) NA 198 (36.1) 330 (19.1) 
 Normal ≥7 to <9 h, n (%) NA 406 (74.4) 892 (61.3) 2324 (74.3) NA 613 (68.2) NA NA 651 (55.0) NA 3174 (74.7) 1281 (52.4) NA 287 (52.4) 1330 (76.9) 
 Long ≥9 h, n (%) NA 41 (7.5) 173 (11.9) 234 (7.5) NA 55 (6.1) NA NA 41 (5.3) NA 225 (5.3) 130 (11.5) NA 63 (11.5) 70 (4.0) 
Dietary intake                
 Total energy (kcal/day) 2,587 ± 928 2,302 ± 776 1,999 ± 613 2,538 ± 903 1,755 ± 619 1,878 ± 616 2,115 ± 1,181 2,490 ± 558 2,212 ± 807 2,066 ± 601 2,332 ± 870 1,602 ± 705 2,007 ± 520 1,922 ± 1,046 2,383 ± 761 
 Total fats (% total energy) 34 ± 6.7 31.2 ± 4.9 32.3 ± 6 31.3 ± 4.8 30.8 ± 7.4 29.5 ± 6.4 35.4 ± 6.9 32.4 ± 4.6 32.5 ± 5.2 31 ± 5.2 32.6 ± 7 32.1 ± 7.3 36.0 ± 6 36.1 ± 6 32.9 ± 4.8 
 PUFA (% total energy) 6.9 ± 2 5.5 ± 1.2 7.5 ± 2.2 5.4 ± 1.1 4.5 ± 1.4 5.7 ± 1.6 6.9 ± 2.2 4.1 ± 0.7 5.2 ± 1.2 3.4 ± 0.7 5.3 ± 1.5 6.3 ± 1.8 6.9 ± 2.7 5 ± 1.4 5.3 ± 1.1 
 MUFA (% total energy) 12.5 ± 2.8 11.1 ± 2.2 11.6 ± 2.5 11.1 ± 2.1 12 ± 3.2 11 ± 2.6 13.3 ± 2.8 9.8 ± 1.7 10.8 ± 2.1 15.6 ± 3.2 10.6 ± 2.7 12.3 ± 3 12.3 ± 2.6 16.2 ± 3.6 11.1 ± 2 
 SFA (% total energy) 11.9 ± 3.2 11.3 ± 2.4 10.3 ± 2.2 11.8 ± 2.5 11.3 ± 3.2 10.1 ± 2.8 11.8 ± 2.8 13.5 ± 2.4 11.8 ± 2.5 10.4 ± 2.2 12.3 ± 7.9 10.8 ± 3.4 14.1 ± 3 11.8 ± 3.2 11.8 ± 2.4 
 CHO (% total energy) 50.2 ± 7.8 48.1 ± 6.2 52.2 ± 7.9 49 ± 6.1 50.8 ± 10.1 51.6 ± 8.2 48.9 ± 8.7 41.8 ± 6.6 47 ± 6.5 51.4 ± 6.7 47.5 ± 7.9 50.7 ± 9 43.3 ± 6.9 44.5 ± 8.7 45.9 ± 5.8 
 Proteins (% total energy) 15 ± 2.6 17.6 ± 2.6 19.1 ± 3.2 17.7 ± 2.5 18.3 ± 4.1 17.1 ± 3.1 15.8 ± 2.9 15.2 ± 2.3 17.1 ± 2.5 15.7 ± 2.1 15.1 ± 2.6 15.7 ± 3.3 16.9 ± 3 17.9 ± 2.4 17.5 ± 2.4 

Values are means ± SD or percentages.

NA, not available.

*Cohort study name (study acronym) (country): Coronary Artery Risk Development in Young Adults (CARDIA) (U.S.); Corogene Controls (Finland); Cardiovascular Health Study (CHS) (U.S.); Dietary, Lifestyle, and Genetic determinants of Obesity and Metabolic syndrome (DILGOM) (Finland); Family Heart Study (FamHS) (U.S.); Framingham Offspring Study (FOS) (U.S.); Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) (U.S.); GOYA MALE (Denmark); Helsinki Birth Cohort Study (HBCS) (Finland); Invecchiare in Chianti (aging in the Chianti Area; InCHIANTI) (Italy); Inter99 (Denmark); Multi-Ethnic Study of Atherosclerosis (MESA) (U.S.); The Rotterdam Study (the Netherlands); The Hellenic Study of Interactions between SNPs and Eating in Atherosclerosis Susceptibility (THISEAS) (Greece); Cardiovascular Risk in Young Finns Study (YFS) (Finland).

†Weekday or workday self-reported sleep duration as usual hours of sleep per night.

Associations of Dietary Intake/Sleep Duration With Cardiometabolic Traits

After adjustment for age, sex, BMI, and study site, we identified significant main associations between dietary intake/sleep duration with glycemic traits (Supplementary Fig. 1 and Supplementary Table 6). Each additional 1% of CHO intake was associated with 0.007 mmol/L lower FG (β ± SE = −0.007 ± 0.001 mmol/L, P = 2.7E-29) and a 0.001 lower HOMA-IR (β ± SE = −0.001 ± 0.0004 [ln], P = 0.002). Each additional 1% of total fat intake was associated with a 0.01 mmol/L higher FG (β ± SE = 0.01 ± 0.001 mmol/L, P = 8.36E-13) and a 0.005 higher HOMA-IR (β ± SE = 0.005 ± 0.0005 [ln], P = 2.76E-16). Similar trends were evident for MUFA and SFA intakes: each additional 1% of MUFA intake was associated with a 0.01 mmol/L higher FG (β ± SE = 0.01 ± 0.001 mmol/L, P = 4.23E-12) and a 0.01 higher HOMA-IR (β ± SE = 0.01 ± 0.001 [ln], P = 1.12E-13), whereas each additional 1% of SFA intake was associated with a 0.01 mmol/L higher FG (β ± SE = 0.01 ± 0.001 mmol/L, P = 9.26E-13) and a 0.01 higher HOMA-IR (β ± SE = 0.01 ± 0.001 [ln], P = 1.66E-19). However, no associations were evident for PUFA intake.

When sleep was assessed linearly, we observed a marginal association between sleep duration and FG: each additional hour of sleep was associated with 0.01 mmol/L lower FG (β ± SE = −0.01 ± 0.005 mmol/L, P = 0.05). However, when sleep was categorical, we observed a significant association between short sleep duration (<7 h) and FG: short sleep duration was associated with a 0.03 mmol/L higher FG (β ± SE = 0.03 ± 0.01 mmol/L, P = 0.01) and with a 0.45 kg/m2 higher BMI (β ± SE = 0.45 ± 0.09 kg/m2, P < 0.001) relative to normal sleep duration (≥7 to <9 h). In addition, long sleep duration (≥9 h) was associated with a 0.33 kg/m2 higher BMI (β ± SE = 0.33 ± 0.14 kg/m2, P = 0.02) and a 0.76 cm higher waist circumference (β ± SE = 0.76 ± 0.18 cm, P < 0.001) and was marginally associated with a higher HOMA-IR (β ± SE = 0.03 ± 0.02 [ln], P = 0.05), relative to normal sleep duration (≥7 to <9 h).

Additional main associations of dietary intake and sleep duration on anthropometric and lipid traits are presented in Supplementary Table 6.

Associations of SNPs With Cardiometabolic Traits

Sensitivity analyses indicated that substantial heterogeneity (I2 > 30%) was introduced by one cohort (GOLDN) for glycemic trait outcomes; consequently, GOLDN was excluded from the association and interaction meta-analyses for glycemic traits.

Main associations of selected SNPs on cardiometabolic traits are presented in Supplementary Table 7. We replicated previously published associations between MTNR1B variants and glycemic traits in the present meta-analysis (18,19). In short, MTNR1B-rs1387153 was associated with FG (β ± SE = 0.058 ± 0.007 mmol/L per additional T allele, P = 1.7E-17), whereas MTNR1B-rs10830963 was associated with FG (β ± SE = 0.1 ± 0.008 mmol/L per additional G allele, P = 4.2E-35) and HOMA-IR (β ± SE = 0.016 ± 0.004 [ln] per additional G allele, P = 0.004). We did not replicate previous associations between CRY2-rs11605924 and FG (P = 0.06) (18,19). Consistent with previous findings, no associations were observed for CLOCK-rs1801260 or NR1D1-rs2314339 on glycemic traits. Main associations of selected SNPs on anthropometric and lipid traits are presented in Supplementary Table 7.

Interactions Between Dietary Intake and Selected SNPs on Cardiometabolic Traits

Meta-analyzed estimates of the interactions between dietary intake and selected SNPs on cardiometabolic traits are presented in Table 2. We observed no interactions at the prespecified Bonferroni-corrected significance level of P < 0.003. We observed a nominal interaction (i.e., P < 0.05) between CHO intake and MTNR1B-rs1387153 for FG (β ± SE = 0.003 ± 0.001 mmol/L, P = 0.01), which suggests a 0.003 mmol/L higher FG with each additional 1% of CHO intake in the presence of the effect T allele (Supplementary Fig. 2A and Supplementary Fig. 3A). In other words, although higher CHO intake is associated with lower FG when evaluated independently of genotype, the protective association of a higher CHO intake was 0.003 mmol/L weaker (per 1% difference in CHO intake) in the presence of each additional T allele, implying that the T allele attenuates the inverse association between CHO intake and FG. Another nominal interaction was evident for MUFA intake and the same MTNR1B variant for FG (β ± SE = −0.007 ± 0.003 mmol/L, P = 0.02), which suggests 0.007 mmol/L lower FG with each additional 1% of MUFA intake in the presence of the effect T allele (Supplementary Fig. 2B). In addition, a nominal interaction was present between total fat intake and NR1D1-rs2314339 on HOMA-IR (β ± SE = 0.0024 ± 0.001 [ln], P = 0.04) (Supplementary Fig. 2C), suggesting a 0.0024 higher HOMA-IR with each additional 1% of total fat intake in the presence of the effect T allele.

Table 2

Meta-analyzed interactions between dietary intake and SNPs on cardiometabolic traits*

Glycemic traitsAnthropometric traitsLipid trait
FG (mmol/L)ln-HOMA-IRBMI (kg/m2)Waist circumference (cm)HDL-C (mmol/L)
SNPGeneAllelesnβ ± SEPβ ± SEPβ ± SEPβ ± SEPβ ± SEP
Total fat (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.003 ± 0.002 0.12 0.0013 ± 0.001 0.29 0.0006 ± 0.001 0.33 0.0002 ± 0.0003 0.55 −2.6E-05 ± 2.1E-05 0.14 
 rs11605924 CRY2 A/C 18,389 −0.001 ± 0.001 0.28 0.0012 ± 0.001 0.09§ −0.0005 ± 0.001 0.40§ 0.0001 ± 0.0002 0.93 −1.3E-05 ± 1.8E-05 0.51§ 
 rs1387153 MTNR1B T/C 26,532 −0.002 ± 0.001 0.09 0.0006 ± 0.001 0.46§ −0.0010 ± 0.001 0.09§ −0.0001 ± 0.0003 0.58 5.2E-05 ± 2.1E-05 0.02 
 rs10830963 MTNR1B G/C 27,607 −0.001 ± 0.002 0.72 −0.0005 ± 0.001 0.60 −0.0004 ± 0.001 0.52 −0.0003 ± 0.0003 0.31 4.7E-05 ± 2.1E-05 0.02 
 rs2314339 NR1D1 T/C 25,016 −0.002 ± 0.002 0.39 0.0024 ± 0.001 0.04 0.0010 ± 0.001 0.09 0.0002 ± 0.0003 0.58 1.6E-05 ± 2.6E-05 0.55 
PUFA (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.005 ± 0.005 0.39§ −0.0003 ± 0.004 0.95 0.0022 ± 0.002 0.25 0.0011 ± 0.0010 0.15 −5.2E-05 ± 5.2E-05 0.49§ 
 rs11605924 CRY2 A/C 18,389 −0.003 ± 0.004 0.41 0.0002 ± 0.003 0.94 −0.0014 ± 0.002 0.39§ 0.0001 ± 0.0010 0.94 −4.1E-05 ± 5.2E-05 0.46 
 rs1387153 MTNR1B T/C 26,532 0.001 ± 0.004 0.91§ 0.0013 ± 0.003 0.66§ −0.0009 ± 0.002 0.65 0.0004 ± 0.0010 0.64 1.0E-04 ± 5.2E-05 0.06 
 rs10830963 MTNR1B G/C 27,607 0.003 ± 0.005 0.61 −0.0002 ± 0.003 0.96 0.0001 ± 0.002 0.97 −0.0002 ± 0.0010 0.74 7.8E-05 ± 5.2E-05 0.17 
 rs2314339 NR1D1 T/C 25,016 −0.001 ± 0.008 0.89 0.0047 ± 0.005 0.30 0.0010 ± 0.003 0.59 0.0004 ± 0.0010 0.69 6.5E-05 ± 7.8E-05 0.42 
MUFA (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.006 ± 0.004 0.11 0.0047 ± 0.003 0.09 0.0008 ± 0.002 0.58 0.0002 ± 0.0010 0.79 −6.5E-05 ± 5.2E-05 0.19 
 rs11605924 CRY2 A/C 18,389 −0.001 ± 0.003 0.60 0.0019 ± 0.002 0.31 −0.0016 ± 0.001 0.21 0.0002 ± 0.0010 0.78 −1.6E-05 ± 5.2E-05 0.69§ 
 rs1387153 MTNR1B T/C 26,532 −0.007 ± 0.003 0.02 0.0008 ± 0.002 0.67§ −0.0035 ± 0.002 0.02§ −0.0003 ± 0.0010 0.67 1.0E-04 ± 5.2E-05 0.01 
 rs10830963 MTNR1B G/C 27,607 −0.002 ± 0.003 0.58 −0.0015 ± 0.002 0.50 −0.0020 ± 0.002 0.31 −0.0007 ± 0.0010 0.30 1.0E-04 ± 5.2E-05 0.02 
 rs2314339 NR1D1 T/C 25,016 −0.005 ± 0.006 0.35 0.0038 ± 0.003 0.16 0.0026 ± 0.002 0.17 0.0003 ± 0.0010 0.71 5.2E-05 ± 5.2E-05 0.52 
SFA (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.006 ± 0.004 0.13§ 0.0025 ± 0.003 0.35 0.0020 ± 0.002 0.31 0.0001 ± 0.0010 0.92 −8.5E-05 ± 5.2E-05 0.11§ 
 rs11605924 CRY2 A/C 18,389 −0.003 ± 0.002 0.21 −0.0001 ± 0.002 0.97§ 0.0003 ± 0.001 0.84§ 0.0001 ± 0.0010 0.99 −1.0E-05 ± 5.2E-05 0.82 
 rs1387153 MTNR1B T/C 26,532 −0.003 ± 0.003 0.21 0.0013 ± 0.002 0.45 −0.0030 ± 0.002 0.04§ −0.0009 ± 0.0010 0.17 7.5E-05 ± 5.2E-05 0.15 
 rs10830963 MTNR1B G/C 27,607 −0.001 ± 0.003 0.66 −0.0016 ± 0.002 0.41 −0.001 ± 0.002 0.38 −0.0007 ± 0.0010 0.26§ 9.6E-05 ± 5.2E-05 0.06 
 rs2314339 NR1D1 T/C 25,016 −0.005 ± 0.005 0.29 0.0045 ± 0.003 0.09 0.0050 ± 0.002 0.01 0.0007 ± 0.0010 0.42 1.3E-05 ± 7.8E-05 0.87 
CHO (% total energy)            
 rs1801260 CLOCK C/T 18,401 −0.002 ± 0.001 0.15 −0.0012 ± 0.001 0.19 −0.0003 ± 0.001 0.53 −0.0002 ± 0.0002 0.47 8.0E-06 ± 1.6E-05 0.64 
 rs11605924 CRY2 A/C 18,389 0.001 ± 0.001 0.57 −0.0007 ± 0.001 0.24§ −0.0001 ± 0.001 0.81§ 0.0001 ± 0.0002 0.86 1.0E-05 ± 1.6E-05 0.49 
 rs1387153 MTNR1B T/C 26,532 0.003 ± 0.001 0.008 −0.0001 ± 0.001 0.87 0.0010 ± 0.001 0.18 0.0001 ± 0.0002 0.97§ −2.8E-05 ± 1.6E-05 0.07 
 rs10830963 MTNR1B G/C 27,607 0.002 ± 0.001 0.08 0.0008 ± 0.001 0.24 0.0001 ± 0.001 0.76 0.0001 ± 0.0002 0.80§ −2.6E-05 ± 1.6E-05 0.08 
 rs2314339 NR1D1 T/C 25,016 0.001 ± 0.002 0.59 −0.0017 ± 0.001 0.05§ −0.0010 ± 0.001 0.03 −0.0002 ± 0.0003 0.51 −4.4E-05 ± 2.1E-05 0.03 
Glycemic traitsAnthropometric traitsLipid trait
FG (mmol/L)ln-HOMA-IRBMI (kg/m2)Waist circumference (cm)HDL-C (mmol/L)
SNPGeneAllelesnβ ± SEPβ ± SEPβ ± SEPβ ± SEPβ ± SEP
Total fat (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.003 ± 0.002 0.12 0.0013 ± 0.001 0.29 0.0006 ± 0.001 0.33 0.0002 ± 0.0003 0.55 −2.6E-05 ± 2.1E-05 0.14 
 rs11605924 CRY2 A/C 18,389 −0.001 ± 0.001 0.28 0.0012 ± 0.001 0.09§ −0.0005 ± 0.001 0.40§ 0.0001 ± 0.0002 0.93 −1.3E-05 ± 1.8E-05 0.51§ 
 rs1387153 MTNR1B T/C 26,532 −0.002 ± 0.001 0.09 0.0006 ± 0.001 0.46§ −0.0010 ± 0.001 0.09§ −0.0001 ± 0.0003 0.58 5.2E-05 ± 2.1E-05 0.02 
 rs10830963 MTNR1B G/C 27,607 −0.001 ± 0.002 0.72 −0.0005 ± 0.001 0.60 −0.0004 ± 0.001 0.52 −0.0003 ± 0.0003 0.31 4.7E-05 ± 2.1E-05 0.02 
 rs2314339 NR1D1 T/C 25,016 −0.002 ± 0.002 0.39 0.0024 ± 0.001 0.04 0.0010 ± 0.001 0.09 0.0002 ± 0.0003 0.58 1.6E-05 ± 2.6E-05 0.55 
PUFA (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.005 ± 0.005 0.39§ −0.0003 ± 0.004 0.95 0.0022 ± 0.002 0.25 0.0011 ± 0.0010 0.15 −5.2E-05 ± 5.2E-05 0.49§ 
 rs11605924 CRY2 A/C 18,389 −0.003 ± 0.004 0.41 0.0002 ± 0.003 0.94 −0.0014 ± 0.002 0.39§ 0.0001 ± 0.0010 0.94 −4.1E-05 ± 5.2E-05 0.46 
 rs1387153 MTNR1B T/C 26,532 0.001 ± 0.004 0.91§ 0.0013 ± 0.003 0.66§ −0.0009 ± 0.002 0.65 0.0004 ± 0.0010 0.64 1.0E-04 ± 5.2E-05 0.06 
 rs10830963 MTNR1B G/C 27,607 0.003 ± 0.005 0.61 −0.0002 ± 0.003 0.96 0.0001 ± 0.002 0.97 −0.0002 ± 0.0010 0.74 7.8E-05 ± 5.2E-05 0.17 
 rs2314339 NR1D1 T/C 25,016 −0.001 ± 0.008 0.89 0.0047 ± 0.005 0.30 0.0010 ± 0.003 0.59 0.0004 ± 0.0010 0.69 6.5E-05 ± 7.8E-05 0.42 
MUFA (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.006 ± 0.004 0.11 0.0047 ± 0.003 0.09 0.0008 ± 0.002 0.58 0.0002 ± 0.0010 0.79 −6.5E-05 ± 5.2E-05 0.19 
 rs11605924 CRY2 A/C 18,389 −0.001 ± 0.003 0.60 0.0019 ± 0.002 0.31 −0.0016 ± 0.001 0.21 0.0002 ± 0.0010 0.78 −1.6E-05 ± 5.2E-05 0.69§ 
 rs1387153 MTNR1B T/C 26,532 −0.007 ± 0.003 0.02 0.0008 ± 0.002 0.67§ −0.0035 ± 0.002 0.02§ −0.0003 ± 0.0010 0.67 1.0E-04 ± 5.2E-05 0.01 
 rs10830963 MTNR1B G/C 27,607 −0.002 ± 0.003 0.58 −0.0015 ± 0.002 0.50 −0.0020 ± 0.002 0.31 −0.0007 ± 0.0010 0.30 1.0E-04 ± 5.2E-05 0.02 
 rs2314339 NR1D1 T/C 25,016 −0.005 ± 0.006 0.35 0.0038 ± 0.003 0.16 0.0026 ± 0.002 0.17 0.0003 ± 0.0010 0.71 5.2E-05 ± 5.2E-05 0.52 
SFA (% total energy)            
 rs1801260 CLOCK C/T 18,401 0.006 ± 0.004 0.13§ 0.0025 ± 0.003 0.35 0.0020 ± 0.002 0.31 0.0001 ± 0.0010 0.92 −8.5E-05 ± 5.2E-05 0.11§ 
 rs11605924 CRY2 A/C 18,389 −0.003 ± 0.002 0.21 −0.0001 ± 0.002 0.97§ 0.0003 ± 0.001 0.84§ 0.0001 ± 0.0010 0.99 −1.0E-05 ± 5.2E-05 0.82 
 rs1387153 MTNR1B T/C 26,532 −0.003 ± 0.003 0.21 0.0013 ± 0.002 0.45 −0.0030 ± 0.002 0.04§ −0.0009 ± 0.0010 0.17 7.5E-05 ± 5.2E-05 0.15 
 rs10830963 MTNR1B G/C 27,607 −0.001 ± 0.003 0.66 −0.0016 ± 0.002 0.41 −0.001 ± 0.002 0.38 −0.0007 ± 0.0010 0.26§ 9.6E-05 ± 5.2E-05 0.06 
 rs2314339 NR1D1 T/C 25,016 −0.005 ± 0.005 0.29 0.0045 ± 0.003 0.09 0.0050 ± 0.002 0.01 0.0007 ± 0.0010 0.42 1.3E-05 ± 7.8E-05 0.87 
CHO (% total energy)            
 rs1801260 CLOCK C/T 18,401 −0.002 ± 0.001 0.15 −0.0012 ± 0.001 0.19 −0.0003 ± 0.001 0.53 −0.0002 ± 0.0002 0.47 8.0E-06 ± 1.6E-05 0.64 
 rs11605924 CRY2 A/C 18,389 0.001 ± 0.001 0.57 −0.0007 ± 0.001 0.24§ −0.0001 ± 0.001 0.81§ 0.0001 ± 0.0002 0.86 1.0E-05 ± 1.6E-05 0.49 
 rs1387153 MTNR1B T/C 26,532 0.003 ± 0.001 0.008 −0.0001 ± 0.001 0.87 0.0010 ± 0.001 0.18 0.0001 ± 0.0002 0.97§ −2.8E-05 ± 1.6E-05 0.07 
 rs10830963 MTNR1B G/C 27,607 0.002 ± 0.001 0.08 0.0008 ± 0.001 0.24 0.0001 ± 0.001 0.76 0.0001 ± 0.0002 0.80§ −2.6E-05 ± 1.6E-05 0.08 
 rs2314339 NR1D1 T/C 25,016 0.001 ± 0.002 0.59 −0.0017 ± 0.001 0.05§ −0.0010 ± 0.001 0.03 −0.0002 ± 0.0003 0.51 −4.4E-05 ± 2.1E-05 0.03 

*Additive allele mode, adjusted for age, sex, BMI (except when assessing BMI outcome), study site (in CARDIA; CHS; FamHS; GOLDN; InCHIANTI; MESA), family or population structure (in Corogene Controls; DILGOM; FamHS; FOS; GOLDN; MESA; The Rotterdam Study; YFS), and genotype batch (in FamHS). Interaction coefficients are shown as β ± SE. β represents the direction and magnitude of the change in outcome trait with each additional minor allele, per each additional % of macronutrient intake. Boldface type indicates nominally significant values (P < 0.05).

†Alleles presented as effect/noneffect alleles.

‡The number of independent observations in each interaction analysis. Exact numbers of observations vary depending on availability of phenotype and genotype information and are presented in Supplementary Table 8.

§I2 >30%, where I2 represents the heterogeneity statistic presented as percent. Exact I2 is presented in Supplementary Table 8.

For anthropometric and lipid traits, we observed nominal interactions (i.e., P < 0.05), including that between SFA intake and NR1D1-rs2314339 for BMI (β ± SE = 0.005 ± 0.002 kg/m2, P = 0.01); the interaction suggests a 0.005 kg/m2 higher BMI with each additional 1% of SFA intake in the presence of the effect T allele (Supplementary Fig. 1D).

Interactions Between Sleep Duration and Selected SNPs on Cardiometabolic Traits

Meta-analyzed estimates of the interactions between sleep duration (continuous and categorical) and selected SNPs on cardiometabolic traits are presented in Table 3. We observed no interactions at the prespecified Bonferroni-corrected significance level of P < 0.003. A nominal interaction was evident between sleep duration and CRY2-rs11605924 for HDL-C (β ± SE = 0.010 ± 0.004 mmol/L, P = 0.005), suggesting a 0.010 mmol/L higher HDL-C with each additional hour of sleep in the presence of the effect A allele (Supplementary Fig. 2E and Supplementary Fig. 3B). That is, in the presence of each additional A allele, the protective association of higher sleep duration on HDL-C was 0.01 mmol/L stronger (per 1 h of additional sleep), such that the A allele appears to strengthen the positive association observed with longer sleep duration. No interactions were evident between categories of sleep duration and this variant for HDL-C (short sleep duration, P = 0.15; long sleep duration, P = 0.21).

Table 3

Meta-analyzed interactions between sleep duration (continuous and categorical) and SNPs on cardiometabolic traits*

Glycemic traitsAnthropometric traitsLipid trait
FG (mmol/L)ln-HOMA-IRBMI (kg/m2)Waist circumference (cm)HDL-C (mmol/L)
SNPGeneAllelesnβ ± SEPβ ± SEPβ ± SEPβ ± SEPβ ± SEP
Continuous, h             
 rs1801260 CLOCK C/T 10,707 0.005 ± 0.02 0.67§ −0.004 ± 0.01 0.61 −0.07 ± 0.07 0.32 0.08 ± 0.09 0.38 0.001 ± 0.005 0.82§ 
 rs11605924 CRY2 A/C 10,706 −0.001 ± 0.01 0.85 −0.004 ± 0.01 0.52 0.02 ± 0.05 0.66 −0.07 ± 0.06 0.24 0.010 ± 0.004 0.005 
 rs1387153 MTNR1B T/C 19,911 0.006 ± 0.01 0.42 0.007 ± 0.01 0.23 0.04 ± 0.05 0.35 0.09 ± 0.06 0.14 −0.001 ± 0.004 0.78 
 rs10830963 MTNR1B G/C 19,913 0.002 ± 0.01 0.85§ −0.004 ± 0.01 0.60 −0.03 ± 0.06 0.58 0.06 ± 0.07 0.41 −0.006 ± 0.004 0.13 
 rs2314339 NR1D1 T/C 18,404 0.007 ± 0.02 0.64 0.015 ± 0.01 0.19§ −0.14 ± 0.08 0.10 0.03 ± 0.11 0.78 0.002 ± 0.006 0.79 
Short (<7 h)            
 rs1801260 CLOCK C/T 2,158 −0.010 ± 0.03 0.70§ 0.005 ± 0.23 0.82 −0.06 ± 0.19 0.60 −0.11 ± 0.23 0.65 −0.008 ± 0.01 0.57 
 rs11605924 CRY2 A/C 3,294 −0.002 ± 0.02 0.92§ 0.006 ± 0.01 0.66§ 0.02 ± 0.12 0.89§ −0.08 ± 0.15 0.58 0.010 ± 0.01 0.15 
 rs1387153 MTNR1B T/C 2,158 0.010 ± 0.02 0.49 −0.001 ± 0.01 0.92§ 0.25 ± 0.12 0.04 −0.23 ± 0.15 0.13 0.005 ± 0.01 0.62 
 rs10830963 MTNR1B G/C 3,525 0.020 ± 0.02 0.38§ 0.020 ± 0.02 0.13§ −0.08 ± 0.15 0.18§ −0.11 ± 0.18 0.53 −0.005 ± 0.01 0.62 
 rs2314339 NR1D1 T/C 3,525 0.0004 ± 0.04 0.88 −0.030 ± 0.03 0.26 0.45 ± 0.23 0.05 0.35 ± 0.28 0.22 −0.003 ± 0.02 0.86 
Long (≥9 h)             
 rs1801260 CLOCK C/T 522 −0.006 ± 0.04 0.89 −0.030 ± 0.04 0.48 0.14 ± 0.31 0.65 −0.30 ± 0.41 0.47 −0.020 ± 0.02 0.45 
 rs11605924 CRY2 A/C 906 0.002 ± 0.03 0.95 0.010 ± 0.02 0.68 0.27 ± 0.21 0.19 −0.03 ± 0.27 0.92 −0.010 ± 0.02 0.21 
 rs1387153 MTNR1B T/C 522 0.020 ± 0.03 0.54 0.010 ± 0.02 0.70 0.60 ± 0.20 0.003 −0.40 ± 0.25 0.10§ 0.010 ± 0.03 0.57 
 rs10830963 MTNR1B G/C 961 0.008 ± 0.03 0.79 0.030 ± 0.03 0.30 −0.20 ± 0.24 0.42§ 0.20 ± 0.30 0.50§ −0.020 ± 0.02 0.64 
 rs2314339 NR1D1 T/C 961 0.070 ± 0.06 0.49 −0.040 ± 0.05 0.50§ 0.44 ± 0.40 0.28 0.58 ± 0.50 0.25 0.010 ± 0.02 0.83 
Glycemic traitsAnthropometric traitsLipid trait
FG (mmol/L)ln-HOMA-IRBMI (kg/m2)Waist circumference (cm)HDL-C (mmol/L)
SNPGeneAllelesnβ ± SEPβ ± SEPβ ± SEPβ ± SEPβ ± SEP
Continuous, h             
 rs1801260 CLOCK C/T 10,707 0.005 ± 0.02 0.67§ −0.004 ± 0.01 0.61 −0.07 ± 0.07 0.32 0.08 ± 0.09 0.38 0.001 ± 0.005 0.82§ 
 rs11605924 CRY2 A/C 10,706 −0.001 ± 0.01 0.85 −0.004 ± 0.01 0.52 0.02 ± 0.05 0.66 −0.07 ± 0.06 0.24 0.010 ± 0.004 0.005 
 rs1387153 MTNR1B T/C 19,911 0.006 ± 0.01 0.42 0.007 ± 0.01 0.23 0.04 ± 0.05 0.35 0.09 ± 0.06 0.14 −0.001 ± 0.004 0.78 
 rs10830963 MTNR1B G/C 19,913 0.002 ± 0.01 0.85§ −0.004 ± 0.01 0.60 −0.03 ± 0.06 0.58 0.06 ± 0.07 0.41 −0.006 ± 0.004 0.13 
 rs2314339 NR1D1 T/C 18,404 0.007 ± 0.02 0.64 0.015 ± 0.01 0.19§ −0.14 ± 0.08 0.10 0.03 ± 0.11 0.78 0.002 ± 0.006 0.79 
Short (<7 h)            
 rs1801260 CLOCK C/T 2,158 −0.010 ± 0.03 0.70§ 0.005 ± 0.23 0.82 −0.06 ± 0.19 0.60 −0.11 ± 0.23 0.65 −0.008 ± 0.01 0.57 
 rs11605924 CRY2 A/C 3,294 −0.002 ± 0.02 0.92§ 0.006 ± 0.01 0.66§ 0.02 ± 0.12 0.89§ −0.08 ± 0.15 0.58 0.010 ± 0.01 0.15 
 rs1387153 MTNR1B T/C 2,158 0.010 ± 0.02 0.49 −0.001 ± 0.01 0.92§ 0.25 ± 0.12 0.04 −0.23 ± 0.15 0.13 0.005 ± 0.01 0.62 
 rs10830963 MTNR1B G/C 3,525 0.020 ± 0.02 0.38§ 0.020 ± 0.02 0.13§ −0.08 ± 0.15 0.18§ −0.11 ± 0.18 0.53 −0.005 ± 0.01 0.62 
 rs2314339 NR1D1 T/C 3,525 0.0004 ± 0.04 0.88 −0.030 ± 0.03 0.26 0.45 ± 0.23 0.05 0.35 ± 0.28 0.22 −0.003 ± 0.02 0.86 
Long (≥9 h)             
 rs1801260 CLOCK C/T 522 −0.006 ± 0.04 0.89 −0.030 ± 0.04 0.48 0.14 ± 0.31 0.65 −0.30 ± 0.41 0.47 −0.020 ± 0.02 0.45 
 rs11605924 CRY2 A/C 906 0.002 ± 0.03 0.95 0.010 ± 0.02 0.68 0.27 ± 0.21 0.19 −0.03 ± 0.27 0.92 −0.010 ± 0.02 0.21 
 rs1387153 MTNR1B T/C 522 0.020 ± 0.03 0.54 0.010 ± 0.02 0.70 0.60 ± 0.20 0.003 −0.40 ± 0.25 0.10§ 0.010 ± 0.03 0.57 
 rs10830963 MTNR1B G/C 961 0.008 ± 0.03 0.79 0.030 ± 0.03 0.30 −0.20 ± 0.24 0.42§ 0.20 ± 0.30 0.50§ −0.020 ± 0.02 0.64 
 rs2314339 NR1D1 T/C 961 0.070 ± 0.06 0.49 −0.040 ± 0.05 0.50§ 0.44 ± 0.40 0.28 0.58 ± 0.50 0.25 0.010 ± 0.02 0.83 

*Additive allele mode, adjusted for age, sex, BMI (except when assessing BMI outcome), study site (in CHS; InCHIANTI; MESA), family or population structure (in Corogene Controls; DILGOM; FOS; MESA; YFS), and genotype batch (in FamHS). Interaction coefficients are shown as β ± SE. β represents the direction and magnitude of the change in outcome trait with each additional minor allele, per each additional hour of sleep or compared with the reference sleep group (≥7 to <9 h). Boldface type indicates nominally significant values (P < 0.05).

†Alleles presented as effect/noneffect alleles.

‡The number of independent observations in each interaction analysis. Exact numbers of observations vary depending on availability of phenotype and genotype information and are presented in Supplementary Table 8.

§I2 >30%, where I2 represents the heterogeneity statistic presented as percent. Exact I2 is presented in Supplementary Table 8.

Finally, we observed a nominal interaction between short sleep duration (<7 h) and MTNR1B-rs1387153 for BMI (β ± SE = 0.25 ± 0.12 kg/m2, P = 0.04) (Supplementary Fig. 2F) and a stronger interaction between long sleep duration (≥9 h) and the same variant for BMI (β ± SE = 0.60 ± 0.20 kg/m2, P = 0.003) (Supplementary Fig. 2G); these interactions suggest 0.25 and 0.60 kg/m2 higher BMIs with short and long sleep durations, respectively, in the presence of the effect T allele, relative to normal sleep duration (≥7 to <9 h).

Results from meta-regressions and sensitivity analyses did not substantively (P > 0.05) affect association or interaction results, neither did they reveal any clear sources of heterogeneity between dietary intake/sleep duration on glycemic, anthropometric, or lipid traits (results not shown).

Our meta-analyses of 15 cohorts did not identify statistically significant gene-environment interactions between modifiable lifestyle factors (relative CHO and fat intakes and sleep duration) and common circadian-related gene variants on glycemic and related metabolic traits. However, the observed nominally significant interactions suggest that lower CHO intake and normal sleep duration may ameliorate cardiometabolic abnormalities conferred by the investigated variants. Until further examination of these interactions is conducted, recommendations should continue to emphasize favorable lifestyle behaviors applicable to the general population, namely, maintaining a healthy diet and obtaining sufficient sleep (≥7 to <9 h) for reducing the risk of metabolic disorders associated with these variants.

The dietary intake associations we observed suggest that diets with a higher percentage of total energy intake from CHO and a lower percentage of total energy intake from total fat, MUFA, and SFA are associated with lower FG. The FG-raising associations of MUFA and SFA intakes observed in the current study are twice the magnitude of the FG-lowering effect of CHO intake, emphasizing the importance of limiting fat intake and in line with current dietary recommendations for individuals with diabetes (26). The similar associations of MUFA and SFA on FG may be partly driven by the strong correlation between the two macronutrients in some cohorts rather than the intake of each dietary component independently. Despite not taking CHO quality into account, our observations that diets higher in relative CHO intake are associated with lower FG are consistent with results from meta-analyses for whole grains (27). Furthermore, the association observed between short sleep duration (<7 h) and higher FG supports previously reported associations in single cohorts (28) and supports previously reported associations between short sleep duration and type 2 diabetes (7). As such, lifestyle recommendations should include dietary modifications related to higher CHO and lower fat composition and achieving normal sleep durations (≥7 to <9 h).

Our evaluation of gene-environment interactions suggest novel putative interactions that fell short of our Bonferroni-corrected cut point but are supported by biological plausibility and may be important for understanding the etiology of type 2 diabetes. The strongest nominal interaction for glycemic traits was an interaction of CHO intake and the MTNR1B-rs1387153 variant on FG, which suggests that every 1% increase in CHO intake exacerbates the FG-raising effect of the T allele that interacts similarly with CHO as another clock variant in CRY1 (5). Along with other nominal interactions observed for MTNR1B-rs1387153, these findings suggest lower CHO and higher MUFA intakes for lower FG among those with the effect T allele. The high frequency of the effect T allele among individuals of European descent (minor allele frequency = 0.28) and the consistency of the rs1387153-FG association across different ethnicities (29) warrant further investigation of this interaction and examination of the potential role of CHO quality in the FG-raising effect of the T allele. These nominal findings, in conjunction with confirmed associations between two common MTNR1B variants and FG of effect sizes similar to those reported earlier (18,30), indicate that continuing efforts to identify lifestyle modifications that offset this genetic risk should remain an important area of active research. Consistent with previous findings, no interactions were evident between sleep duration and the selected variants on glycemic traits (7).

Our previous work suggests that variants identified through GWAS or candidate gene association studies for type 2 diabetes may show gene-environment interactions for related metabolic traits beyond glycemic traits (3). We identified three nominal interactions that are supported by previous reports and biological plausibility. The first is a nominal interaction between SFA intake and NR1D1-rs2314339 on BMI. The obesity-associated NR1D1 gene encodes the nuclear receptor rev-erb-α, which plays a critical role in metabolism and was reported to respond to dietary MUFA for the outcome of BMI (21,31,32). The second is a nominal interaction between FG-associated CRY2-rs11605924 and sleep duration for HDL-C. We observed a positive association between HDL-C and sleep only when considered in the context of CRY2, a clock gene that inhibits CLOCK:BMAL1-mediated transcription of genes involved in lipid metabolism (33). Supporting the circadian control of HDL-C are results from the Global Lipids Genetics Consortium GWAS for HDL-C for this CRY2 variant (β ± SE = 0.0004 ± 0.0001 mmol/L per additional A allele, P < 0.001) (34), which suggest marginal associations between CRY2 and HDL-C (Supplementary Fig. 4). The CRY2 variant is in linkage disequilibrium with rs6843722 (r2 = 1.00 in the 1000 Genomes Project data set), a CRY2 intronic variant that was shown to abolish the upregulation of CRY2 expression in human peripheral blood mononuclear cells after sleep restriction and has functional evidence to affect transcriptional regulation of CCCTC-binding factor and glucocorticoid receptor, two transcription factors associated with HDL-C (3537). Therefore, it is possible that short sleep duration results in differences in CRY2 expression, influencing CRY2 control of downstream pathways, namely HDL-C. Finally, we observed nominal interactions between short and long sleep duration, both of which are associated with higher BMI, and MTNR1B-rs1387153 on BMI, suggesting even higher BMI with short and long sleep duration among carriers of the effect T allele. This interaction provides additional support for the potential role of sleep duration in modifying the associations between circadian-related genetic variants and cardiometabolic outcomes (38) and the importance of normal sleep duration (≥7 to <9 h) for optimal health.

The strengths of the present observational study from 15 cohorts include a large sample size necessary to detect gene-environment interactions. Our collaborative approach enabled us to standardize our analytic approach across cohorts, and despite the wide range of cohorts included in the study, we observed little evidence of heterogeneity in our overall analysis. However the present investigation also has limitations. The reported findings are limited to individuals of European descent, and exploring the interaction in other populations is warranted considering replication of the SNP associations in different ethnicities (29,39).

Our use of self-reported dietary intake and sleep duration was susceptible to reporting bias, and the use of different assessment tools across cohorts could have increased measurement error, biasing our results toward the null (40). In addition, we failed to replicate a previously reported significant association between CRY2 variant and FG, although we observed an effect size similar to that of the discovery GWAS of up to 46,000 individuals (18); it is possible our sample size was too small to replicate the significant associations. Although we have selected circadian-related gene variants showing strong associations with metabolic traits from GWAS and candidate-gene studies, it is possible that interactions might be observable for other circadian-related variants.

Lastly, these cross-sectional meta-analyses of observational studies can only lead us toward hypotheses rather than demonstrate the temporal relationships or causal pathways linking clock genes, diet, or sleep, with glycemic, anthropometric, and lipid traits. Other studies are required to establish these mechanistic links, including studies of genetic modification of the effects of experimental changes in diet composition or sleep duration.

Our findings contribute to the understanding of how lifestyle can reduce the risk of type 2 diabetes and cardiometabolic disorders in genetically susceptible individuals. Results from the present large observational study from 15 cohorts suggest the potential presence of selected common circadian-related gene-environment interactions on metabolic traits. The nominal interaction between CHO intake and the MTNR1B variant on FG, suggesting that CHO intake could exacerbate the FG-raising effects of this well-studied MTNR1B variant, the evidence supporting the role for CRY2 in HDL-C control and its responsiveness to sleep duration, and the interaction between long sleep duration and MTNR1B variant on BMI, suggesting that the association between long sleep duration and higher BMI is exacerbated among carriers of the MTNR1B variant, are interesting and merit further study. Mechanistic examinations of the novel nominal interactions and further investigations in larger cohorts are necessary before personalized recommendations are framed for individuals at genetic risk for metabolic disruption. Moreover, the observed associations between diet—specifically, higher CHO and lower fat composition—and normal sleep duration (≥7 to <9 h) on glycemic traits—particularly FG—suggest that emphasis should be placed on these modifiable lifestyle factors to offset the growing prevalence of type 2 diabetes and cardiovascular disease.

Appendix

Members of the CHARGE Nutrition Study Group. Hassan S. Dashti, PhD, Nutrition and Genomics Laboratory, Jean Mayer U.S. Department of Agriculture Human Nutrition Research Center on Aging at Tufts University, Boston, MA; Jack L. Follis, PhD, Department of Mathematics, Computer Science, and Cooperative Engineering, University of St. Thomas, Houston, TX; Caren E. Smith, MS, DVM, Nutrition and Genomics Laboratory, Jean Mayer U.S. Department of Agriculture Human Nutrition Research Center on Aging at Tufts University, Boston, MA; Toshiko Tanaka, PhD, Translational Gerontology Branch, National Institute on Aging, Baltimore, MD; Marta Garaulet, PhD, Department of Physiology, University of Murcia, Murcia, Spain; Daniel J. Gottlieb, MD, MPH, Division of Sleep and Circadian Disorders, Brigham and Women’s Hospital, Boston, MA, Division of Sleep Medicine, Harvard Medical School, Boston, MA, and Sleep Disorders Center, VA Boston Healthcare System, Boston, MA; Adela Hruby, PhD, MPH, Department of Nutrition, Harvard School of Public Health, Boston, MA; Paul F. Jacques, ScD, Nutritional Epidemiology Laboratory, Jean Mayer U.S. Department of Agriculture Human Nutrition Research Center on Aging at Tufts University, Boston, MA; Jessica C. Kiefte-de Jong, RD, PhD, Department of Epidemiology, Erasmus Medical Center, Rotterdam, the Netherlands, and Global Public Health, Leiden University College, The Hague, the Netherlands; Stefania Lamon-Fava, MD, PhD, Cardiovascular Nutrition Laboratory, Jean Mayer U.S. Department of Agriculture Human Nutrition Research Center on Aging at Tufts University, Boston, MA; Frank A.J.L. Scheer, PhD, Division of Sleep and Circadian Disorders, Brigham and Women’s Hospital, Boston, MA, and Division of Sleep Medicine, Harvard Medical School, Boston, MA; Traci M. Bartz, MS, Cardiovascular Health Research Unit, Department of Medicine, University of Washington, Seattle, WA, and Department of Biostatistics, University of Washington, Seattle, WA; Leena Kovanen, MS, Department of Mental Health and Substance Abuse Services, National Institute for Health and Welfare (THL), Helsinki, Finland; Mary K. Wojczynski, PhD, Department of Genetics, Washington University School of Medicine, St. Louis, MO; Alexis C. Frazier-Wood, PhD, U.S. Department of Agriculture/Agricultural Research Service Children's Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX; Tarunveer S. Ahluwalia, PhD, The Novo Nordisk Foundation Centre for Basic Metabolic Research, Section of Metabolic Genetics, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark, and Copenhagen Prospective Studies on Asthma in Childhood, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark, and Danish Pediatric Asthma Centre, Gentofte Hospital, The Capital Region, Copenhagen, Denmark; Mia-Maria Perälä, MSc, Department of Chronic Disease Prevention, National Institute for Health and Welfare (THL), Helsinki, Finland; Anna Jonsson, PhD, The Novo Nordisk Foundation Centre for Basic Metabolic Research, Section of Metabolic Genetics, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark; Taulant Muka, MD, MSc, DSc, Department of Epidemiology, Erasmus Medical Center, Rotterdam, the Netherlands; Ioanna P. Kalafati, MSc, Department of Nutrition and Dietetics, Harokopio University, Athens, Greece; Vera Mikkilä, PhD, Department of Food and Environmental Sciences, Division of Nutrition, University of Helsinki, Helsinki, Finland, and Research Centre of Applied and Preventive Cardiovascular Medicine, University of Turku, Turku, Finland; Rozenn N. Lemaitre, PhD, MPH, Cardiovascular Health Research Unit, Department of Medicine, University of Washington, Seattle, WA; Timo Partonen, MD, PhD, Department of Mental Health and Substance Abuse Services, National Institute for Health and Welfare (THL), Helsinki, Finland; Tapani Ebeling, MD, PhD, Oulu University Hospital, Department of Internal Medicine, Division of Endocrinology, Oulu, Finland; Paul N. Hopkins, MD, MSPH, School of Medicine, University of Utah, Salt Lake City, UT; Lavinia Paternoster, PhD, MRC Integrative Epidemiology Unit, School of Social and Community Medicine, University of Bristol, Bristol, U.K.; Jari Lahti, PhD, Institute of Behavioural Sciences, University of Helsinki, Helsinki, Finland, and Folkhälsan Research Centre, Helsinki, Finland; Dena G. Hernandez, MS, Laboratory of Neurogenetics, National Institute on Aging, Baltimore, MD; Ulla Toft, PhD, Research Centre for Prevention and Health, Glostrup University Hospital, Glostrup, Denmark; Richa Saxena, PhD, Division of Sleep and Circadian Disorders, Brigham and Women’s Hospital, Boston, MA, and Center for Human Genetic Research and Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, MA; Anna Vitezova, PharmD, MSc, Department of Epidemiology, Erasmus Medical Center, Rotterdam, the Netherlands; Stavroula Kanoni, PhD, William Harvey Research Institute, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, London, U.K.; Olli T. Raitakari, PhD, Research Centre of Applied and Preventive Cardiovascular Medicine, University of Turku, Turku, Finland, and Department of Clinical Physiology and Nuclear Medicine, University of Turku and Turku University Hospital, Turku, Finland; Bruce M. Psaty, MD, PhD, Cardiovascular Health Research Unit, Department of Medicine, University of Washington, Seattle, WA, Departments of Epidemiology and Health Services, University of Washington, Seattle, WA, and Group Health Research Institute, Group Health, Seattle, WA; Markus Perola, MD, PhD, Department of Chronic Disease Prevention, National Institute for Health and Welfare (THL), Helsinki, Finland; Satu Männistö, PhD, Department of Chronic Disease Prevention, National Institute for Health and Welfare (THL), Helsinki, Finland; Robert J. Straka, BSc, PharmD, FCCP, Experimental and Clinical Pharmacology, College of Pharmacy, University of Minnesota, MN; Torben Hansen, MD, PhD, The Novo Nordisk Foundation Centre for Basic Metabolic Research, Section of Metabolic Genetics, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark; Katri Räikkönen, PhD, Institute of Behavioural Sciences, University of Helsinki, Helsinki, Finland; Luigi Ferrucci, MD, PhD, Translational Gerontology Branch, National Institute on Aging, Baltimore, MD; Niels Grarup, MD, PhD, The Novo Nordisk Foundation Centre for Basic Metabolic Research, Section of Metabolic Genetics, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark; W. Craig Johnson, MS, Department of Biostatistics, University of Washington, Seattle, WA; Loukianos Rallidis, MD, Second Department of Cardiology, University General Hospital “Attikon,” Athens, Greece; Mika Kähönen, PhD, Department of Clinical Physiology, Tampere University Hospital and University of Tampere, Tampere, Finland; David S. Siscovick, MD, MPH, New York Academy of Medicine, New York, NY; Aki S. Havulinna, DSc, Department of Chronic Disease Prevention, National Institute for Health and Welfare (THL), Helsinki, Finland; Arne Astrup, DMSc, Department of Nutrition, Exercise, and Sports, Faculty of Science, University of Copenhagen, Copenhagen, Denmark; Torben Jørgensen, MD, Research Centre for Prevention and Health, Capital Region, Institute of Public Health, University of Copenhagen, Copenhagen, Denmark, and Faculty of Medicine, University of Aalborg, Aalborg, Denmark; Tzu-An Chen, PhD, U.S. Department of Agriculture/Agricultural Research Service, Children's Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX; Albert Hofman, MD, PhD, Department of Epidemiology, Erasmus Medical Center, Rotterdam, the Netherlands; Panos Deloukas, PhD, William Harvey Research Institute, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, London, U.K., and Princess Al-Jawhara Al-Brahim Centre of Excellence in Research of Hereditary Disorders (PACER-HD), King Abdulaziz University, Jeddah, Saudi Arabia; Jorma S.A. Viikari, PhD, Department of Medicine, University of Turku, Turku, Finland, and Division of Medicine, Turku University Hospital, Turku, Finland; Dariush Mozaffarian, MD, DrPH, Friedman School of Nutrition Science and Policy, Tufts University, Boston, MA; Oluf Pedersen, MD, PhD, The Novo Nordisk Foundation Centre for Basic Metabolic Research, Section of Metabolic Genetics, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark; Jerome I. Rotter, MD, Institute for Translational Genomics and Population Sciences, Los Angeles Biomedical Research Institute at Harbor-UCLA Medical Center, Torrance, CA; André G. Uitterlinden, PhD, Department of Internal Medicine, Erasmus Medical Center, Rotterdam, the Netherlands; Ilkka Seppälä, MSc, Department of Clinical Chemistry, Fimlab Laboratories, University of Tampere School of Medicine, Tampere, Finland; Henning Tiemeier, MD, PhD, Department of Epidemiology, Erasmus Medical Center, Rotterdam, the Netherlands; Veikko Salomaa, MD, PhD, Department of Chronic Disease Prevention, National Institute for Health and Welfare (THL), Helsinki, Finland; Sina A. Gharib, MD, Cardiovascular Health Research Unit, Department of Medicine, University of Washington, Seattle, WA, and Computational Medicine Core, Center for Lung Biology, UW Medicine Sleep Center, University of Washington, Seattle, WA; Ingrid B. Borecki, PhD, Department of Genetics, Washington University School of Medicine, St. Louis, MO; Donna K. Arnett, PhD, Department of Epidemiology, School of Public Health, University of Alabama at Birmingham, Birmingham, AL; Thorklid I.A. Sørensen, DrMedSci, The Novo Nordisk Foundation Centre for Basic Metabolic Research, Section of Metabolic Genetics, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark, and Institute of Preventive Medicine, Bispebjerg and Frederiksberg Hospitals, The Capital Region, Copenhagen, Denmark; Johan G. Eriksson, PhD, MD, Department of Chronic Disease Prevention, National Institute for Health and Welfare (THL), Helsinki, Finland, Folkhälsan Research Centre, Helsinki, Finland, Department of General Practice and Primary Health Care, University of Helsinki, Helsinki, Finland, Helsinki University Central Hospital, Unit of General Practice, Helsinki, Finland, and Vasa Central Hospital, Vasa, Finland; Stefania Bandinelli, MD, Geriatric Unit, Azienda Sanitaria Firenze (ASF), Florence, Italy; Allan Linneberg, MD, PhD, Research Centre for Prevention and Health, Glostrup University Hospital, Glostrup, Denmark, Department of Clinical Experimental Research, Glostrup University Hospital, Glostrup, Denmark, and Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark; Stephen S. Rich, PhD, Center for Public Health Genomics, University of Virginia, Charlottesville, VA; Oscar H. Franco, MD, PhD, Department of Epidemiology, Erasmus Medical Center, Rotterdam, the Netherlands; George Dedoussis, PhD, Department of Nutrition and Dietetics, Harokopio University, Athens, Greece; Terho Lehtimäki, PhD, Department of Clinical Chemistry, Fimlab Laboratories, University of Tampere School of Medicine, Tampere, Finland; José M. Ordovás, PhD, Nutrition and Genomics Laboratory, Jean Mayer U.S. Department of Agriculture Human Nutrition Research Center on Aging at Tufts University, Boston, MA, Department of Epidemiology, Centro Nacional de Investigaciones Cardiovasculares (CNIC), Madrid, Spain, and Instituto Madrileño de Estudios Avanzados en Alimentación (IMDEA-FOOD), Madrid, Spain.

Clinical trial reg. nos. NCT00005130 (CARDIA), NCT00005133 (Cardiovascular Health Study), NCT00005136 (Family Heart Study), NCT00005121 (Framingham Offspring Study), NCT00083369 (Genetic and Environmental Determinants of Triglycerides), NCT01331512 (InCHIANTI Study), NCT00289237 (Inter99), and NCT00005487 (Multi-Ethnic Study of Atherosclerosis), clinicaltrials.gov.

Funding. Infrastructure for the CHARGE Consortium is supported in part by the National Heart, Lung, and Blood Institute grant HHSN268200800007C. Support was provided by the U.S. Department of Agriculture, under agreement no. 58-1950-0-014, American Heart Association (grant 14BGIA18740011), Academy of Finland (grant 117787), Italian Ministry of Health (grant ICS110.1/RF97.71), National Institute of Diabetes and Digestive and Kidney Diseases (grant DK063491), National Center for Advancing Translational Sciences (grant UL1TR000124). C.E.S. is supported by K08-HL-112845-01. F.A.J.L.S. was supported in part by National Institutes of Health grants R21-DK-089378 and R01-HL-094806. Cohort-specific sources of support and acknowledgments are presented in Supplementary Table 1. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the authors and do not necessarily reflect the view of the U.S. Department of Agriculture.

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

Author Contributions. H.S.D., C.E.S., T.T., M.G., D.J.G., F.A.J.L.S., and J.M.O. designed the study. H.S.D., J.L.F., T.T., T.M.B., L.K., M.K.W., A.C.F.-W., T.S.A., M.-M.P., A.J., T.M., I.P.K., and V.M. conducted research and contributed to statistical analyses. H.S.D., C.E.S., M.G., D.J.G., A.H., P.F.J., S.L.-F., F.A.J.L.S., and J.M.O. interpreted data. H.S.D., J.L.F., C.E.S., T.T., M.G., D.J.G., A.H., P.F.J., J.C.K.-d.J., S.L.-F., F.A.J.L.S., and J.M.O. wrote the manuscript. All authors read and approved the final version of the manuscript. H.S.D. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Prior Presentation. Preliminary results were presented at the 78th Scientific Sessions and Annual Meeting of the American Society for Nutrition in conjunction with Experimental Biology 2014, San Diego, CA, 26–30 April 2014, and at the 2014 Society for Research on Biological Rhythms Meeting, Big Sky, MT, 14–18 June 2014.

1.
Franks
PW
,
Pearson
E
,
Florez
JC
.
Gene-environment and gene-treatment interactions in type 2 diabetes: progress, pitfalls, and prospects
.
Diabetes Care
2013
;
36
:
1413
1421
[PubMed]
2.
Ordovas
JM
,
Corella
D
.
Nutritional genomics
.
Annu Rev Genomics Hum Genet
2004
;
5
:
71
118
[PubMed]
3.
Parnell
LD
,
Blokker
BA
,
Dashti
HS
, et al
.
CardioGxE, a catalog of gene-environment interactions for cardiometabolic traits
.
BioData Min
2014
;
7
:
21
[PubMed]
4.
Garaulet
M
,
Lee
YC
,
Shen
J
, et al
.
CLOCK genetic variation and metabolic syndrome risk: modulation by monounsaturated fatty acids
.
Am J Clin Nutr
2009
;
90
:
1466
1475
[PubMed]
5.
Dashti
HS
,
Smith
CE
,
Lee
YC
, et al
.
CRY1 circadian gene variant interacts with carbohydrate intake for insulin resistance in two independent populations: Mediterranean and North American
.
Chronobiol Int
2014
;
31
:
660
667
[PubMed]
6.
Zheng
JS
,
Arnett
DK
,
Parnell
LD
, et al
.
Genetic variants at PSMD3 interact with dietary fat and carbohydrate to modulate insulin resistance
.
J Nutr
2013
;
143
:
354
361
[PubMed]
7.
Tare
A
,
Lane
JM
,
Cade
BE
, et al
.
Sleep duration does not mediate or modify association of common genetic variants with type 2 diabetes
.
Diabetologia
2014
;
57
:
339
346
[PubMed]
8.
American Diabetes Association
.
Diagnosis and classification of diabetes mellitus
.
Diabetes Care
2014
;
37
(
Suppl. 1
):
S81
S90
[PubMed]
9.
Voight
BF
,
Scott
LJ
,
Steinthorsdottir
V
, et al.;
MAGIC investigators
;
GIANT Consortium
.
Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis
.
Nat Genet
2010
;
42
:
579
589
[PubMed]
10.
Kalsbeek
A
,
la Fleur
S
,
Fliers
E
.
Circadian control of glucose metabolism
.
Mol Metab
2014
;
3
:
372
383
[PubMed]
11.
Panda
S
,
Antoch
MP
,
Miller
BH
, et al
.
Coordinated transcription of key pathways in the mouse by the circadian clock
.
Cell
2002
;
109
:
307
320
[PubMed]
12.
Lamia
KA
,
Storch
KF
,
Weitz
CJ
.
Physiological significance of a peripheral tissue circadian clock
.
Proc Natl Acad Sci U S A
2008
;
105
:
15172
15177
[PubMed]
13.
Scheer
FAJL
,
Chan
JL
,
Fargnoli
J
, et al
.
Day/night variations of high-molecular-weight adiponectin and lipocalin-2 in healthy men studied under fed and fasted conditions
.
Diabetologia
2010
;
53
:
2401
2405
[PubMed]
14.
Shea
SA
,
Hilton
MF
,
Orlova
C
,
Ayers
RT
,
Mantzoros
CS
.
Independent circadian and sleep/wake regulation of adipokines and glucose in humans
.
J Clin Endocrinol Metab
2005
;
90
:
2537
2544
[PubMed]
15.
Marcheva
B
,
Ramsey
KM
,
Buhr
ED
, et al
.
Disruption of the clock components CLOCK and BMAL1 leads to hypoinsulinaemia and diabetes
.
Nature
2010
;
466
:
627
631
[PubMed]
16.
Scheer
FAJL
,
Hilton
MF
,
Mantzoros
CS
,
Shea
SA
.
Adverse metabolic and cardiovascular consequences of circadian misalignment
.
Proc Natl Acad Sci U S A
2009
;
106
:
4453
4458
[PubMed]
17.
Turek
FW
,
Joshu
C
,
Kohsaka
A
, et al
.
Obesity and metabolic syndrome in circadian Clock mutant mice
.
Science
2005
;
308
:
1043
1045
[PubMed]
18.
Dupuis
J
,
Langenberg
C
,
Prokopenko
I
, et al.;
DIAGRAM Consortium
;
GIANT Consortium
;
Global BPgen Consortium
;
Anders Hamsten on behalf of Procardis Consortium
;
MAGIC investigators
.
New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk
.
Nat Genet
2010
;
42
:
105
116
[PubMed]
19.
Prokopenko
I
,
Langenberg
C
,
Florez
JC
, et al
.
Variants in MTNR1B influence fasting glucose levels
.
Nat Genet
2009
;
41
:
77
81
[PubMed]
20.
Liu
C
,
Li
H
,
Qi
L
, et al
.
Variants in GLIS3 and CRY2 are associated with type 2 diabetes and impaired fasting glucose in Chinese Hans
.
PLoS One
2011
;
6
:
e21464
[PubMed]
21.
Garaulet
M
,
Smith
CE
,
Gomez-Abellán
P
, et al
.
REV-ERB-ALPHA circadian gene variant associates with obesity in two independent populations: Mediterranean and North American
.
Mol Nutr Food Res
2014
;
58
:
821
829
[PubMed]
22.
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
:
658
662
[PubMed]
23.
Gómez-Abellán
P
,
Madrid
JA
,
Luján
JA
, et al
.
Sexual dimorphism in clock genes expression in human adipose tissue
.
Obes Surg
2012
;
22
:
105
112
[PubMed]
24.
Lorenzo
C
,
Okoloise
M
,
Williams
K
,
Stern
MP
,
Haffner
SM
;
San Antonio Heart Study
.
The metabolic syndrome as predictor of type 2 diabetes: the San Antonio heart study
.
Diabetes Care
2003
;
26
:
3153
3159
[PubMed]
25.
Froy
O
.
The relationship between nutrition and circadian rhythms in mammals
.
Front Neuroendocrinol
2007
;
28
:
61
71
[PubMed]
26.
American Diabetes Association; Bantle JP, Wylie-Rosett J, Albright AL, et al. Nutrition recommendations and interventions for diabetes: a position statement of the American Diabetes Association. Diabetes Care 2008;31(Suppl. 1):S61–S78
27.
Nettleton
JA
,
McKeown
NM
,
Kanoni
S
, et al.;
MAGIC Investigators
.
Interactions of dietary whole-grain intake with fasting glucose- and insulin-related genetic loci in individuals of European descent: a meta-analysis of 14 cohort studies
.
Diabetes Care
2010
;
33
:
2684
2691
[PubMed]
28.
Knutson
KL
,
Van Cauter
E
,
Zee
P
,
Liu
K
,
Lauderdale
DS
.
Cross-sectional associations between measures of sleep and markers of glucose metabolism among subjects with and without diabetes: the Coronary Artery Risk Development in Young Adults (CARDIA) Sleep Study
.
Diabetes Care
2011
;
34
:
1171
1176
[PubMed]
29.
Chambers
JC
,
Zhang
W
,
Zabaneh
D
, et al
.
Common genetic variation near melatonin receptor MTNR1B contributes to raised plasma glucose and increased risk of type 2 diabetes among Indian Asians and European Caucasians
.
Diabetes
2009
;
58
:
2703
2708
[PubMed]
30.
Bouatia-Naji
N
,
Bonnefond
A
,
Cavalcanti-Proença
C
, et al
.
A variant near MTNR1B is associated with increased fasting plasma glucose levels and type 2 diabetes risk
.
Nat Genet
2009
;
41
:
89
94
[PubMed]
31.
Goumidi
L
,
Grechez
A
,
Dumont
J
, et al
.
Impact of REV-ERB alpha gene polymorphisms on obesity phenotypes in adult and adolescent samples
.
Int J Obes (Lond)
2013
;
37
:
666
672
[PubMed]
32.
Yang
X
,
Downes
M
,
Yu
RT
, et al
.
Nuclear receptor expression links the circadian clock to metabolism
.
Cell
2006
;
126
:
801
810
[PubMed]
33.
van der Horst
GT
,
Muijtjens
M
,
Kobayashi
K
, et al
.
Mammalian Cry1 and Cry2 are essential for maintenance of circadian rhythms
.
Nature
1999
;
398
:
627
630
[PubMed]
34.
Willer
CJ
,
Schmidt
EM
,
Sengupta
S
, et al.;
Global Lipids Genetics Consortium
.
Discovery and refinement of loci associated with lipid levels
.
Nat Genet
2013
;
45
:
1274
1283
[PubMed]
35.
Boyle
AP
,
Hong
EL
,
Hariharan
M
, et al
.
Annotation of functional variation in personal genomes using RegulomeDB
.
Genome Res
2012
;
22
:
1790
1797
[PubMed]
36.
Lavebratt
C
,
Sjöholm
LK
,
Soronen
P
, et al
.
CRY2 is associated with depression
.
PLoS One
2010
;
5
:
e9407
[PubMed]
37.
Figtree
GA
,
Grieve
SM
,
Speller
B
, et al
.
A commonly occurring polymorphism upstream of the estrogen receptor alpha alters transcription and is associated with increased HDL
.
Atherosclerosis
2008
;
199
:
354
361
[PubMed]
38.
Dashti
HS
,
Follis
JL
,
Smith
CE
, et al
.
Habitual sleep duration is associated with BMI and macronutrient intake and may be modified by CLOCK genetic variants
.
Am J Clin Nutr
2015
;
101
:
135
143
[PubMed]
39.
Kan
MY
,
Zhou
DZ
,
Zhang
D
, et al
.
Two susceptible diabetogenic variants near/in MTNR1B are associated with fasting plasma glucose in a Han Chinese cohort
.
Diabet Med
2010
;
27
:
598
602
[PubMed]
40.
Freedman
LS
,
Schatzkin
A
,
Midthune
D
,
Kipnis
V
.
Dealing with dietary measurement error in nutritional cohort studies
.
J Natl Cancer Inst
2011
;
103
:
1086
1092
[PubMed]

Supplementary data