Genome-wide gene by environment study of time spent in daylight and chronotype identifies emerging genetic architecture underlying light sensitivity ==================================================================================================================================================== * Angus C. Burns * Andrew J. K. Phillips * Martin K. Rutter * Richa Saxena * Sean W. Cain * Jacqueline M. Lane ## Abstract Light is the primary stimulus for synchronizing the circadian clock in humans. There are very large interindividual differences in the sensitivity of the circadian clock to light. Little is currently known about the genetic basis for these interindividual differences. We performed a genome-wide gene-by-environment interaction study (GWIS) in 280,897 individuals from the UK Biobank cohort to identify genetic variants that moderate the effect of daytime light exposure on chronotype (individual time of day preference), acting as ‘light sensitivity’ variants for the impact of daylight on the circadian system. We identified a genome-wide significant SNP mapped to the *ARL14EP* gene (rs3847634; *p* < 5×10−8), where additional minor alleles were found to enhance the morningness effect of daytime light exposure (*β**GxE* = -.03, SE = 0.005) and were associated with increased gene *ARL14EP* expression in brain and retinal tissues. Gene-property analysis showed light sensitivity loci were enriched for genes in the G protein-coupled glutamate receptor signaling pathway and in *Per2*+ hypothalamic neurons. Linkage disequilibrium score regression identified significant genetic correlations of the light sensitivity GWIS with chronotype and sleep duration, such that greater light sensitivity was associated with later chronotype, greater insomnia symptoms and shorter sleep duration. Greater light sensitivity was also genetically correlated with greater risk for PTSD. This study is the first to assess light as an important exposure in the genomics of chronotype and is a critical first step in uncovering the genetic architecture of human circadian light sensitivity and its links to sleep and mental health. Keywords * Light sensitivity * gene by environment * GWAS * chronotype * daylight * circadian rhythms * sleep * mood disorders * PTSD ## Introduction In humans, daily rhythms in gene expression, physiology, and behavior are regulated by the circadian clock located in the suprachiasmatic nuclei1 (SCN). The circadian clock is primarily synchronized by ocular light exposure2. These rhythms are the overt manifestation of the activity of well-described molecular feedback-loops. Light information transmitted to the SCN from the retina via the retinohypothalamic tract (RHT) is translated into shifts in the phase of molecular biological clocks3 by altering the expression of core clock genes4. There are large interindividual differences in the sensitivity of the human circadian clock to light. The circadian phase shifting response of humans to the same light stimulus can vary by as much as two hours5 and recent work has shown that the melatonin suppression response, a measure of circadian light sensitivity, can vary between individuals by over an order of magnitude6. These interindividual differences play a role in the aetiology of delayed sleep-wake phase disorder (DSWPD)7, major depressive disorder8,9, bipolar disorder 10,11 and other sleep disturbances12. Light sensitivity has also been related to tolerance of shift work, with individuals varying in their ability to adjust their biological clock so that their biological day aligns with the environmental night13-15. Circadian light sensitivity is a heritable trait16, as measured by twin MZ-DZ correlations. Chronotype, a proxy-measure of circadian phase17, has significant SNP-heritability (SNP *h*2 range = .12-.37) in published GWAS and loci for this trait are enriched for circadian clock and glutamate signaling pathways18,19. These findings indicate that both circadian light sensitivity and the downstream traits it impacts are in part determined by genetic variation in the human population. A common variable-number tandem repeat polymorphism in the *PER3* gene has been associated with light sensitivity20 in a small sample; however, no other work has yet examined the genomics of circadian light sensitivity. The timing of light exposure has a differential impact on the timing of circadian rhythms, as described by the human light phase-response curve (PRC). Light exposure in the late night and early day advances the rhythm while light exposure in the evening and early night delays the rhythm21. Consistent with the light PRC, daytime light exposure tends to produce advances in the phase of the circadian clock and thus is robustly associated with earlier chronotype and greater ease of awakening22-24. Here, we leverage the phenotypic association between daytime light exposure and chronotype to conduct a Genome-Wide Interaction Study (GWIS) of loci that dampen or enhance this effect. Using unrelated individuals from the UK Biobank, we test for gene-by-environment (GxE) interactions between SNP alleles and daytime light exposure on chronotype to identify the genetic architecture of light sensitivity. ## Methods ### UK Biobank and genetic data The UK Biobank prospective general population cohort contains more than 502,000 UK residents (aged 37–73 years; ∼54% women) recruited via National Health Service (NHS) patient registers from 2006 to 2010, with the study population described in detail elsewhere25,26. Participants provided extensive demographic, lifestyle, health, mood, and physical information via assessment and touch-screen questionnaires as well as physiological samples for the purpose of genome-wide array genotyping. Participants who accepted the invitation to join the UK Biobank cohort provided written, informed consent and the UK Biobank has generic ethical approval from the North West Multi-Center Research Ethics Committee (ref 11/NW/03820). The current analyses were conducted under UK Biobank application number 6818 (Martin Rutter). In brief, blood, saliva, and urine were collected from participants, and DNA was extracted from the buffy coat samples. Participant DNA was genotyped on two arrays, Affymetrix UK BiLEVE Axiom array (initial 50,000 participants) and the Affymetrix UK Biobank Axiom Array with >95% common content and genotypes for ∼800,000 autosomal SNPs imputed to two reference panels (Haplotype Reference Consortium (HRC) and UK10K haplotype resource). Genotypes were called using Affymetrix Power Tools software. Sample and SNPs for quality control were selected from a set of 489,212 samples across 812,428 unique markers. Sample quality control (QC) was conducted using 605,876 high-quality autosomal markers. Samples were removed for high missingness or heterozygosity (968 samples) and sex chromosome abnormalities (652 samples). Genotypes for 488,377 samples passed sample QC (∼99.9% of total samples). Marker-based QC measures were tested in the European ancestry subset (n = 463,844), which was identified based on principal components of ancestry. SNPs were tested for batch effects (197 SNPs/ batch), plate effects (284 SNPs/batch), Hardy– Weinberg equilibrium (572 SNPs/ batch), sex effects (45 SNPs/batch), array effects (5417 SNPs), and discordance across control replicates (622 on UK BiLEVE Axiom array and 632 UK Biobank Axiom array; P value <10−12 or <95% for all tests). For each batch (106 batches total), markers that failed at least one test were set to missing. Before imputation, 805,426 SNPs passed QC in at least one batch (>99% of the array content). Population structure was captured by principal component analysis on the samples using a subset of high-quality (missingness < 1.5%), high-frequency (>2.5%) SNPs (∼100,000 SNPs) and identified the subsample of white British descent. In addition to the calculated population structure by the UK Biobank, we locally further clustered participants into four ancestry clusters using K-means clustering on the principal components, identifying 453,964 participants of European ancestry. Related individuals were identified by estimating kinship coefficients for all pairs of samples, using only markers weakly informative of ancestral background. For the current analysis, individuals of non-white ethnicity and related individuals were excluded to limit confounding effects. The UK Biobank centrally imputed autosomal SNPs to UK10K haplotype, 1000 Genomes Phase 3, and Haplotype Reference Consortium (HRC). Autosomal SNPs were pre-phased using SHAPEIT3 and imputed using IMPUTE4. In total ∼96 million SNPs were imputed.. For all analyses, we used ∼15.0M HRC imputed variants with an imputation r2 ≥ 0.3, MAF ≥ 0.001 (0.1%) and with a Hardy–Weinberg equilibrium (HWE) P>1×10-12 (chi-squared; 1 degree of freedom). We excluded non-HRC imputed variants. Further details on the UK Biobank genotyping, quality control and imputation procedures can be found elsewhere27. ### Measurement of time spent in daylight, chronotype, and other self-reported phenotypes At the baseline assessment, all participants self-reported how many hours they spent outdoors during the day on a typical day in both summer and winter (UK Biobank data-fields 1050 and 1060). Participants reported an integer using a touch-screen number pad or selected among alternate options including “Less than an hour a day”, “Do not know”, or “Prefer not to answer”. Initial cleaning of the data involved excluding participants that rated “Prefer not to answer” or “Do not know”, re-coding “Less than an hour a day” as zero (*n* = 98,431) and excluding values larger than the typical day length in the UK during summer (16 hours; *n* = 253) and winter (8 hours; *n* = 5,474). Summer and winter reports were strongly correlated (Pearson’s *r* = .65), indicating people who spent more time in outdoor light in winter also tended to do so in summer. Furthermore, a subset of approximately 20,000 participants completed a repeat assessment of this question from 2012-2013. Both summer and winter reports were strongly correlated (both *r* > .6) across repeat assessments, indicating good reliability of the measure. To yield a single measure of time spent in outdoor light, winter and summer reports were averaged within participants. For analysis, extreme bins (9-12 h spent outdoors) with low density (total *n* = 1,016 (0.2%) were collapsed into a single bin. Finally, due to the rightward skew in the distribution of the time spent in outdoor light, a log10 transformation was applied. We extracted self-reported chronotype from the baseline assessment (“Morningness-Eveningness”; data-field 1180) which was measured with the question “Do you consider yourself to be?” with six options: “Definitely a ‘morning’ person”, “More a ‘morning’ than ‘evening’ person”, “More an ‘evening’ than a ‘morning’ person”, “Definitely an ‘evening’ person”, “Do not know” or “Prefer not to answer”. The latter two categories were coded as missing. Participants further self-reported age and sex. Finally, weight and height were measured, and body-mass index (BMI) was calculated as weight (kg)/height2 (m2; UK Biobank data field 21001). ### Genome-wide association studies for light sensitivity in the UK Biobank A genome-wide interaction study (GWIS) for circadian light sensitivity was conducted in unrelated individuals of white British ancestry (*n* = 280,897) to reduce confounding by stratification and kinship using the –gxe function in PLINK28. The logarithm of time spent in daylight (continuous) was used as the moderating exposure for the genotype effect on chronotype (continuous) with the following equation: ![Formula][1] Here, SNP refers to the allelic dosage at a given SNP coded additively, daylight refers to the log-transformed average time spent outdoors during the day and the covariates included in the model were age, sex, genotyping array (UK BiLEVE or UK Biobank Axiom), BMI, and the top ten principal components of ancestry to control for population stratification. The interaction term (*β**GxE*) was used to identify “light sensitivity” variants, or variants that moderate the strength of the effect of daylight exposure on chronotype. Negative values of *β**GxE* weights indicate additional minor alleles for a given SNP enhance the morningness effect of light on chronotype, whereas positive *β**GxE* weights indicate additional minor alleles dampen the morningness effect of light on chronotype. In supplement to the GWIS, we conducted two separate GWAS of chronotype stratified by time spent in daylight to examine whether light exposure modifies the genetic architecture of chronotype. Low (≤1.5 h; *n* = 121,922) and high (≥3.5 h; *n* = 118,429) daytime light exposure groups were derived by applying a quartile split of the time spent in daylight variable, with the two intermediary quartiles excluded. Within these stratified samples, genetic association analysis across the autosomes was performed in European participants with BOLT-LMM29 linear mixed models allowing for relatedness using an additive genetic model adjusting for age, sex, ten principal components of ancestry, genotyping array, and genetic correlation matrix. ![Formula][2] A maximum SNP missingness of 10% and per-sample missingness of 40% was tolerated. A SNP imputation quality threshold of 0.80 and minor allele frequency of 0.001 (0.1%) was used for all genome-wide analyses. A genome-wide significance threshold of 5×10−8 was used to identify significant SNPs. ### Gene, Pathway, Tissue, and Single-Cell Expression Post-GWIS gene based and pathway analyses were carried out with MAGMA30 in FUMA31. Gene based analysis was performed using default MAGMA parameters and a SNP-wide mean model of SNPs mapped to 18,960 protein coding genes (Bonferroni *p* = 0.05/18960 = 2.64×10−6). Pathway analysis was conducted using MAGMA30 gene-set analysis in FUMA, which tests 10,894 gene sets (curated gene sets = 4,728; gene ontology terms = 6,166) for enrichment using the full distribution of SNP *p* values from the GWIS (Bonferroni *p* = 0.05/10,651 = 4.69×10−6). Tissue enrichment analyses (GTEx V8) were conducted using MAGMA gene-property analysis in FUMA using default parameters to test the relationship between tissue-specific gene-expression profiles and the GWIS light sensitivity associations. Cell-type specificity analysis was conducted in the same manner, using single-cell RNA-seq gene-expression profiles as the outcome, on the three circadian hypothalamic cell subtypes (*Per2**+*, *VIP**+* & *NMS**+*) identified in Romanov, Zeisel, Bakker, Girach, Hellysaz, Tomer, Alpár, Mulder, Clotman, Keimpema, Hsueh, Crow, Martens, Schwindling, Calvigioni, Bains, Máté, Szabó, Yanagawa, Zhang, Rendeiro, Farlik, Uhlén, Wulff, Bock, Broberger, Deisseroth, Hökfelt, Linnarsson, Horvath and Harkany 32. Queries of expression, methylation and histone modification QTLs for genome-wide significant loci were done using QTLbase33 and EyeGEx34. ### Genetic Correlations Linkage Disequilibrium (LD) score regression was implemented in the R package *GenomicSEM**35* to examine the genetic correlation (*r*g) between GWIS light sensitivity genetic architecture (indexed by *β**GxE*) and sleep, circadian, psychiatric, and metabolic traits with publicly available GWAS data that have been linked phenotypically to circadian light sensitivity36-39. A second set of genetic correlations was conducted for stratified GWAS of chronotype in high and low daylight exposure groups and the above traits to examine whether the exposure moderated the association of chronotype with these traits. Genetic correlations for high and low daylight exposure chronotype GWAS were compared using a χ2 test of a constrained 1-df model, where the genetic correlation of chronotype in the high and low daylight groups was constrained to be equal against an unconstrained 0-df model in *GenomicSEM*. ## Results ### Phenotypic results The median time spent in daylight among UK Biobank participants was 2.5 h (IQR = 1.5-3.5 hours) and time spent in daylight was correlated with chronotype, such that greater daylight exposure was associated with greater morningness (*r* = -.10, *p* < 5×10−150). This association has been previously reported24 and is independent of lifestyle, demographic, and employment covariates. Demographic characteristics of the GWIS sample are presented in Supplemental Table 1. ### Genome-wide interaction study The Manhattan and Q-Q plots for the GWIS of light sensitivity SNPs using 14,055,103 imputed variants in unrelated white British UK Biobank participants (*n* = 280,897; λGC = 1.04) are presented in Figure 1a-b. A single genome-wide significant hit, rs3847634 (see Supplemental Figure 1 for locus plot), was identified that mapped to the *ARL14EP/FSHB* genomic region on chromosome 11 (*chr11:30,344,598-30,359,774*), *β**GxE* = -.03, SE = 0.005, *Z* = -5.48, *p* = 4.37 ×10−8), such that additional T allelic dosage conferred a stronger morningness effect of daylight (see Figure 1c; see Supplemental Table 2 for top GWIS hits). Sensitivity analysis that removed the adjustment for BMI replicated this finding (*β**GxE* = -.03, *p* = 4.5×10−8). The marginal effect of this SNP in the GWIS analysis was positive and genome-wide significant (*β**G* = .04, SE = .007, *Z* = 6.42, *p* = 1.35 ×10−10), indicating that, independent of daytime light exposure, additional T alleles conferred greater eveningness. Importantly, this marginal effect was not identified in a prior GWAS of chronotype18 (rs3847634 β = 0.01, SE = 0.003, *p* = 0.0004) that did not take into account daytime light exposure and its interaction with genotype. An independent *ARL14EP/MPPED2* locus (rs621421; *r*2 with rs3847634 = .008, 1000 Genomes EUR population) has been associated with chronotype18 in prior GWAS; however, neither the marginal nor interaction effect for this locus reached significance in the GWIS (*p*both > .05). The *p* values for the marginal effect of SNPs on chronotype in the GWIS and the published chronotype GWAS18 effect correlated strongly (*r* = .56) and the genetic correlation between the two indicated agreement (*r*g = 1.01). The GxE effect from the GWIS for all previously reported chronotype loci are reported in Supplemental Table 3. No other novel genome-wide significant marginal effects were identified in the GWIS as compared to the prior chronotype GWAS. ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/30/2022.06.29.22277078/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2022/06/30/2022.06.29.22277078/F1) Figure 1 Genome-wide interaction study of light sensitivity. (A) Manhattan and (B) Q-Q plots of genome-wide interaction study (GWIS) for daylight exposure on chronotype, horizontal red dotted line indicates genome-wide significance threshold p < 5×10−8. The genome-wide significant SNP rs3847634 at the *ARL14EP* locus and SNPs in LD with it are highlighted in red (C) Conditional regression plot demonstrating the interaction between genotype at rs3847634 and daylight on chronotype, such that additional T alleles confer a stronger morningness effect of daylight. ### QTL, Gene, and Pathway Analysis The rs3847634 SNP is an expression quantitative trait locus (eQTL; QTLbase & EyeGEx) for *ARL14EP* in the hippocampus (β = .061, *p* = 1.15×10−7) and retina (β = .44, *p* = 5.49×10−10), such that additional T alleles were associated with increased *ARL14EP* expression in these tissues. It was also associated with increased *FSHB* gene expression in the brain (β = .25, *p* = 5.23×10−12; Supplemental Tables 4 & 5). Analysis of histone and methylation QTLs revealed that the T allele is associated with increased H3K27ac histone modifications (β = .06, *p* = 2.52×10−6) and reduced methylation (cg06241208; chr11:30,344,200; β = -.34, *p* = 3.76×10−16; Supplemental Table 6) at the promoter region of the *ARL14EP* gene. MAGMA gene-based tests combining SNP-level information within a gene did not identify any Bonferroni-significant genes (see Supplemental Figure 2; Supplemental Table 7). However, suggestive signals were identified in the *HOMER1* (*p =* 1.37×10−5), *GOLIM4* (*p* = 3.34×10−5), and *ARL14EP* (*p* = 4.43×10−5) genes. Pathway analysis combining gene-level information identified a single Bonferroni-significant Gene-Ontology pathway, the G Protein-Coupled Glutamate Receptor Signaling Pathway (*p* = 1.87×10−6; see Supplemental Table 8), which contains 13 genes, including *HOMER1* and the glutamate metabotropic receptor 1 gene (*MgluR1/GRM1*; Gene-based test *p* = 0.006; see Supplemental Table 8). Additionally, genes identified in the light sensitivity GWIS by QTL and positional mapping (RGS12 & RIC8B) significantly overlapped with the Gene Ontology *G Protein Alpha Subunit Binding* (*p*adj = .01; see Supplemental Figure 3). ### Tissue and Single-Cell Expression Analysis MAGMA tissue-expression analysis did not identify significant enrichment in any of the 53 superordinate GTEx tissue types (Supplemental Figure 4), nor any of the brain samples across the stages of development in the BrainSpan database (Supplemental Figure 5). Single-cell RNA-seq gene-property analysis of circadian subtypes of hypothalamic neurons32 revealed a Bonferroni-significant enrichment in *Per2*+ hypothalamic cells (*p* = 0.007; Bonferroni *p* = 0.017), but not in *VIP*+ or *NMS*+ circadian subtypes (Supplemental Table 9). ### Genetic Correlations We performed LD-score regression40 (LDSC) analyses to assess the genetic overlap of light sensitivity (i.e., the genetic architecture of loci moderating the effect of daylight on chronotype) and publicly available sleep, circadian, psychiatric, and metabolic traits (see Figure 2 and Supplemental Table 10). Light sensitivity (*h*2 = .2%, *se* = .1%) had a strong negative genetic correlation with chronotype (*r*g = - .62, se = 0.10, *p* = 1.47×10−9), indicating that SNPs which tend to confer greater light sensitivity (i.e., enhancing the morningness effect of daylight on chronotype) tend to also be eveningness SNPs. Sleep duration also significantly correlated with light sensitivity (*r*g = .35, se = 0.11, *p* = 0.002), indicating that SNPs that tend to confer lower light sensitivity tend to also confer longer sleep duration. There was a negative genetic correlation between light sensitivity and PTSD (*r*g = -.71, se = 0.30, *p* = 0.01), indicating that SNPs associated with higher light sensitivity tended to confer greater risk for PTSD; however, the uncertainty for this estimate was large. Similarly, frequent insomnia had a significant genetic correlation with light sensitivity (*r*g = -.26, se = 0.12, *p* = .03), such that greater light sensitivity was associated with greater insomnia symptomology. No other traits were significantly genetically correlated with light sensitivity. ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/30/2022.06.29.22277078/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2022/06/30/2022.06.29.22277078/F2) Figure 2 Light sensitivity genetic correlations. Genetic correlations between the light sensitivity GWIS and sleep, psychiatric and metabolic traits previously linked to light sensitivity. Note: as negative values of the βGxE term indicate an enhancing effect of the SNP on light, negative genetic correlations denote the correlation between greater light sensitivity with a given trait. We then tested if the amount of daytime light exposure moderates the genetic correlations of chronotype with the same set of sleep, circadian, psychiatric, and metabolic traits by performing LDSC on chronotype GWAS stratified by daylight exposure (bottom quartile, low light ≤ 1.5h, *n* = 121,922; top quartile, high light ≥ 3.5h, *n* = 118,429; see Supplemental Figures 6 & 7 for Manhattan and Q-Q plots). The SNP-*h*2 of chronotype was moderated by light exposure as calculated by BOLT-REML, with heritability being higher in the low daylight exposure group (ChronotypeLow Light *h*2 = 0.20, se = 0.003; ChronotypeHigh Light *h*2 = 0.16, se = 0.003). The genetic correlation between chronotype in those exposed to low and high daytime light was strong (*r*g = 0.93, se = 0.05, *p* = 6.92×10−88). The lead SNP of the low light chronotype GWAS was rs12424741, mapped to the ALG10 locus (*p* = 5.29×10−16) and the lead SNP of the high light chronotype was rs3769124 mapped to the ASB1 locus (*p* = 2.6×10−10), both previously reported18. Consistent with the light sensitivity genetic correlations, there was a significant difference between the genetic correlation of chronotype and sleep duration in the high (*r*g = 0.11, se = 0.03, *p* = 1.5×10−4) versus low (*r*g = 0.02, se = 0.03, *p* = 0.49; χ2 = 6.06, *p* = 0.01) light exposure groups, such that eveningness was associated with longer sleep duration in those exposed to high daytime light but not in those exposed to low daytime light (see Figure 3 and Supplemental Table 11). Similarly, the genetic correlations between chronotype and PTSD (χ2 = 4.81, *p* = 0.03) and chronotype and daytime sleepiness (χ2 = 4.12, *p* = 0.04) were significantly different, such that eveningness was correlated with reduced risk for PTSD and daytime sleepiness among those exposed to high levels of daytime light. Eveningness was more strongly associated with major depressive disorder among those with low daylight exposure (*r*g = 0.10, se = 0.03, *p* = .0004) compared to those with high daylight exposure (*r*g = 0.06, se = 0.03, *p* = .03), though this difference was only suggestive (χ2 = 2.99; *p* = .08). ![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/30/2022.06.29.22277078/F3.medium.gif) [Figure 3](http://medrxiv.org/content/early/2022/06/30/2022.06.29.22277078/F3) Figure 3 Chronotype genetic correlations stratified by daylight exposure. Genetic correlations between chronotype (eveningness) and sleep, psychiatric and metabolic outcomes for GWAS split by environmental exposure to daylight (bottom quartile, blue points, low light ≤ 1.5h, *n* = 121,922; top quartile, red points, high light ≥ 3.5h, *n* = 118,429). Grey points depict the genetic correlations of the published chronotype GWAS with sleep, psychiatric and metabolic outcomes for comparison. ## Discussion Here we report the first study of genome-wide gene by environment interactions for daytime light exposure on chronotype and uncover an emerging genetic architecture of human circadian light sensitivity. We identified a novel locus, candidate genes, and pathways for the modulation of the advancing effect of daytime light exposure on chronotype. We also observed enrichment for light sensitivity genes in hypothalamic *Per2*+ cell subtypes. Genetic correlations indicated that the genetic architecture of light sensitivity significantly overlapped with chronotype and sleep duration, while being nominally associated with sleep and psychiatric disturbance. This study provides important insights into the biological basis of human light sensitivity and its link to sleep and psychiatric traits. We found that one locus, rs3847634, mapped to the *ARL14EP* gene on chromosome 11, reached genome-wide significance in the GWIS. Additional T alleles at this locus were associated with a stronger morningness effect of daytime light exposure on chronotype, thus conferring enhanced light sensitivity. In addition to identifying this GxE effect, we observed that the rs3847634 SNP reached genome-wide significance for its marginal effect on chronotype, such that additional T alleles at this locus were associated with greater eveningness, an effect previously unobserved18,19. An independent SNP in the *ARL14EP* gene (rs621421, *r*2= 0.0075 with rs3847634) that has previously been associated with chronotype18, provides additional evidence for this gene playing a role in circadian timing. This SNP did not interact with light in the present study. Replication of this interaction in independent datasets with light exposure, chronotype and genetic data will be important. This result highlights the importance of considering relevant environmental exposures when examining the genomics of a given trait. While potentially counterintuitive, this result is consistent with a variant that enhances the overall effect of light on the timing of the circadian clock, independent of the timing of light exposure. Enhanced sensitivity to light would boost phase-advancing effects of daytime light exposure, thereby strengthening the association between daylight exposure and morning chronotype. However, the same enhancement of light sensitivity would also boost phase-delaying effects of evening light exposure41. Given modern human lighting conditions are characterized by reduced bright light exposure during the day and bright electric light exposure at night42-45 a SNP that enhances light sensitivity would therefore be expected to confer both of these properties. We found that rs3847634 is associated with reduced CpG methylation and increased h3k27ac histone modifications in the promoter region of the *ARL14EP* gene. Such epigenetic modifications are consistent with a more active gene promoter by diminishing the repressive effects of DNA methylation and opening chromatin to transcriptional machinery, respectively46,47. These modifications may therefore explain our observation that rs3847634 is an eQTL for increased *ARL14EP* gene expression in brain and retinal tissues. Together, these findings hint at a plausible molecular mechanism for the effect of rs3847634 on *ARL14EP* gene expression and subsequent phenotypic effects. The *ARL14EP* gene is important for appropriate brain and axonal development. In humans, *ARL14EP* haploinsufficiency is associated with callosal hypoplasia while its knockdown in animal models results in severe axonal arborization deficits and the disruption of transcollosal connectivity48. WAGR syndrome (Wilms’ tumour, Aniridia, Genital abnormalities, and mental Retardation) is caused by an 11p14-12 deletion that includes the *ARL14EP* gene and is characterized by a photosensitivity and intellectual disability phenotype49. Furthermore, an independent SNP at the *ARL14EP* locus is associated with retinal refractive error50 and non-synonymous mutations in the CRD domain of the ARL14EP protein are associated with neurodevelopmental disorders51,52. There appears to be overlap in biological function of *ARL14EP* with other genes that encode for axon development/guidance (via semaphorins, ephrins). The semaphorin Sema6a is hyper-expressed in *ARL14EP*-knockdown neurons and epigenetic editing of Sema6a gene promoters has been shown to rescue axon function in *ARL14EP* knockdown animal models48. Semaphorins and ephrins have been shown to be important in the development of retinal ganglion cells and a SNP at the Sema6a locus is associated with chronotype18,53,54. While indirect, collectively these data in concert with our findings point toward a putative role for *ARL14EP* in axonal development pertinent to the retina and phototransduction – for example, by altering retinohypothalamic tract development and subsequently the strength of retinal inputs to the circadian clock. Pathway analysis identified enrichment of GWIS SNPs in the G protein-coupled glutamatergic signalling pathway. Nominal enrichment at the gene level was also observed for the *HOMER1* and *MgluR1* genes that are constituents of this pathway. These results are consistent with a large body of research that demonstrates the critical role of glutamatergic signalling in circadian photoreception, in particular via metabotropic glutamate receptors 1 and 5 expressed in the SCN55,56. Similarly, there is strong evidence for the role of HOMER1 in circadian phototransduction. HOMER1 participates in activity-dependent control of metabotropic glutamate receptors, including MgluR1 and MgluR5, is rhythmic in the SCN with a peak in the subjective night when the clock is most sensitive to phase shifts, and its expression is induced by light exposure in the subjective night57-59. Furthermore, cell-type specificity analysis focusing on hypothalamic neurons expressing circadian-related genes32 identified significant enrichment in SCN neurons expressing the core clock gene Per2. Together, these findings provide convergent evidence towards a circadian light sensitivity understanding of the genetic architecture uncovered by the present GWIS. Finally, genetic correlation analyses identified significant genetic overlap between light sensitivity and chronotype, such that greater light sensitivity was associated with greater eveningness. This genetic correlation was strong and is likely explained by a similar phenomenon described above, whereby modern human lighting environments are depleted for daytime light43 and enriched with night-time light44,60 – hence resulting in light sensitivity SNPs being associated with eveningness as individuals with a genetic load for greater light sensitivity experience greater delaying phase shifts to evening light and a later chronotype21,61. Importantly, there is not a statistical dependence between the light sensitivity GxE term in the model and the marginal SNP effect. For example, a recent stress-sensitivity GWIS for depression found zero genetic correlation between stress sensitivity and the marginal SNP effect on depression62. Light sensitivity was also found to genetically correlate with sleep duration, such that greater light sensitivity genetic load was associated with lower sleep duration genetic load. This finding is consistent with greater light sensitivity amplifying the phase delaying effect of evening light exposure, resulting in the delayed timing of sleep and, given a constant wake time, shortened sleep12,44. Greater light sensitivity may also shorten sleep via direct effects on arousal63,64, thereby enhancing the disruptive effects of evening light on subsequent sleep65. A novel association between greater light sensitivity and PTSD diagnosis was also observed. This finding may elucidate recent work demonstrating circadian dysrhythmia is linked to symptom severity in PTSD and that melatonin and phase advancing light may be effective treatments for PTSD symptoms66-68. We note that there are limitations of both the sample and the subjective nature of predictors and outcome measures used in this study. Firstly, there are significant changes in the properties of the circadian system that occur with age69, including earlier phase70, dampened amplitude70 and reduced light sensitivity71-73. Older adults also spend a larger portion of their day in bright (>1000lux) light74. Given these changes and the makeup of the UK Biobank sample being largely middle-older aged adults, it is feasible that the genetic architecture of light sensitivity, defined here as the moderating effect of genotype on the daylight-chronotype relationship, may differ in younger adults. Secondly, the present work relies on a self-report measure of daytime light exposure. While this measure may broadly capture time spent in bright light it does not capture objective intensity and duration of exposure at the ocular level. Future work should address these limitations by testing examining genetic light sensitivity across the lifespan using objective measures of light and, if feasible, circadian phase. Finally, it is important to note that selective gene-environment correlation may provide an alternate explanation for our results, such that individuals with a greater genetic propensity to seek out more daylight may also therefore tend to have earlier chronotype. The present work is the first to identify common genetic variation underlying light sensitivity in the human genome; its links to molecular pathways and cell subtypes in the hypothalamus; and genetic overlap with sleep duration, sleep disturbance and PTSD. This study emphasises the relevance of light as a critical environmental exposure when studying the genomics of circadian rhythms traits and provides the basis for future lines of research to improve the methodology beyond subjective measures of light and circadian phase. ## Supporting information Supplementary Tables [[supplements/277078_file02.xlsx]](pending:yes) Supplementary Figures [[supplements/277078_file03.docx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## Author Contributions *Study conceptualisation and design:* A.C.B., J.M.L., R.S., S.W.C., A.J.K.P., M.K.R.; *Analysis:* A.C.B., J.M.L.; *Overseeing project:* A.C.B., J.M.L., R.S., S.W.C., A.J.K.P.; *Writing and editing:* A.C.B., J.M.L., R.S., S.W.C., A.J.K.P., M.K.R.; *Critical revisions of the manuscript:* A.C.B., J.M.L., R.S., S.W.C., A.J.K.P., M.K.R. ## Disclosure Statement The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. AJKP and SWC have received research funding from Delos and Versalux, and they are co-founders and co-directors of Circadian Health Innovations PTY LTD. SWC has also received research funding from Beacon Lighting and has consulted for Dyson. ## Data Availability Statement This work utilized the UK Biobank resource (application 26209; Martin Rutter). ## Acknowledgements This research and Angus C. Burns was supported by a Research Training Program (RTP) scholarship from the Australian Government. * Received June 29, 2022. * Revision received June 29, 2022. * Accepted June 30, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Buijs RM, Kalsbeek A. Hypothalamic integration of central and peripheral clocks. Nature Reviews Neuroscience. 2001;2(7):521–526. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/35081582&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11433377&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169860200020&link_type=ISI) 2. 2.Boivin DB, Duffy JF, Kronauer RE, Czeisler CA. Dose-response relationships for resetting of human circadian clock by light. Nature. 1996;379(6565):540–542. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/379540a0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8596632&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996TU69300050&link_type=ISI) 3. 3.Reppert SM, Weaver DR. Coordination of circadian timing in mammals. Nature. 2002;418(6901):935–941. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature00965&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12198538&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177677500035&link_type=ISI) 4. 4.Shigeyoshi Y, Taguchi K, Yamamoto S, et al. Light-induced resetting of a mammalian circadian clock is associated with rapid induction of the mPer1 transcript. Cell. 1997;91(7):1043–1053. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(00)80494-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9428526&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071281400020&link_type=ISI) 5. 5.Lockley SW, Brainard GC, Czeisler CA. High Sensitivity of the Human Circadian Melatonin Rhythm to Resetting by Short Wavelength Light. The Journal of Clinical Endocrinology & Metabolism. 2003;88(9):4502–4505. doi:10.1210/jc.2003-030570 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2003-030570&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12970330&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185258700073&link_type=ISI) 6. 6.Phillips AJ, Vidafar P, Burns AC, et al. High sensitivity and interindividual variability in the response of the human circadian system to evening light. Proceedings of the National Academy of Sciences. 2019;116(24):12019–12024. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE2LzI0LzEyMDE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMzAvMjAyMi4wNi4yOS4yMjI3NzA3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 7. 7.Watson LA, Phillips AJK, Hosken IT, et al. Increased sensitivity of the circadian system to light in delayed sleep–wake phase disorder. The Journal of Physiology. 2018;()doi:doi:10.1113/JP275917 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1113/JP275917&link_type=DOI) 8. 8.McGlashan EM, Coleman MY, Vidafar P, Phillips AJK, Cain SW. Decreased sensitivity of the circadian system to light in current, but not remitted depression. Journal of Affective Disorders. Sep 1 2019;256:386–392. doi:10.1016/j.jad.2019.05.076 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jad.2019.05.076&link_type=DOI) 9. 9.McGlashan EM, Nandam LS, Vidafar P, Mansfield DR, Rajaratnam SMW, Cain SW. The SSRI citalopram increases the sensitivity of the human circadian system to light in an acute dose. Psychopharmacology. 2018/11/01 2018;235(11):3201–3209. doi:10.1007/s00213-018-5019-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00213-018-5019-0&link_type=DOI) 10. 10.Bullock B, McGlashan EM, Burns AC, Lu BS, Cain SW. Traits related to bipolar disorder are associated with an increased post-illumination pupil response. Psychiatry Research. 2019/08/01/ 2019;278:35–41. doi:[https://doi.org/10.1016/j.psychres.2019.05.025](https://doi.org/10.1016/j.psychres.2019.05.025) 11. 11.Lewy AJ, Wehr TA, Goodwin FK, Newsome DA, Rosenthal NE. Manic-depressive patients may be supersensitive to light. Lancet (London, England). Feb 14 1981;1(8216):383–4. doi:10.1016/s0140-6736(81)91697-4 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(81)91697-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6110010&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 12. 12.Chellappa SL. Individual differences in light sensitivity affect sleep and circadian rhythms. Sleep. 2021;44(2):zsaa214. doi:10.1093/sleep/zsaa214 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/sleep/zsaa214&link_type=DOI) 13. 13.Boivin DB, Tremblay GM, James FO. Working on atypical schedules. Sleep medicine. 2007;8(6):578–589. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.sleep.2007.03.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17481949&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249934000005&link_type=ISI) 14. 14.Postnova S, Robinson PA, Postnov DD. Adaptation to Shift Work: Physiologically Based Modeling of the Effects of Lighting and Shifts’ Start Time. PLOS ONE. 2013;8(1):e53379. doi:10.1371/journal.pone.0053379 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0053379&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23308206&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 15. 15.Martin SK, Eastman CI. Medium-intensity Light Produces Circadian Rhythm Adaptation to Simulated Night-shift Work. Sleep. 1998;21(2):154–165. doi:10.1093/sleep/21.2.154 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/sleep/21.2.154&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9542799&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000072642500006&link_type=ISI) 16. 16.Hallam KT, Olver JS, Chambers V, Begg DP, McGrath C, Norman TR. The heritability of melatonin secretion and sensitivity to bright nocturnal light in twins. Psychoneuroendocrinology. 2006/08/01/ 2006;31(7):867–875. doi:[https://doi.org/10.1016/j.psyneuen.2006.04.004](https://doi.org/10.1016/j.psyneuen.2006.04.004) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16769177&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 17. 17.Baehr EK, Revelle W, Eastman CI. Individual differences in the phase and amplitude of the human circadian temperature rhythm: with an emphasis on morningness–eveningness. Journal of sleep research. 2000;9(2):117–127. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2869.2000.00196.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10849238&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000087863500002&link_type=ISI) 18. 18.Jones SE, Lane JM, Wood AR, et al. Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nature Communications. 2019/01/29 2019;10(1):343. doi:10.1038/s41467-018-08259-7 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-08259-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30696823&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 19. 19.Lane JM, Vlasac I, Anderson SG, et al. Genome-wide association analysis identifies novel loci for chronotype in 100,420 individuals from the UK Biobank. Nature communications. 2016;7(1):1–10. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms10165&link_type=DOI) 20. 20.Chellappa SL, Viola AU, Schmidt C, et al. Human Melatonin and Alerting Response to Blue-Enriched Light Depend on a Polymorphism in the Clock Gene PER3. The Journal of Clinical Endocrinology & Metabolism. 2012;97(3):E433–E437. doi:10.1210/jc.2011-2391 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2011-2391&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22188742&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301229600015&link_type=ISI) 21. 21.Khalsa SBS, Jewett ME, Cajochen C, Czeisler CA. A phase response curve to single bright light pulses in human subjects. The Journal of physiology. 2003;549(3):945–952. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1113/jphysiol.2003.040477&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12717008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183815200026&link_type=ISI) 22. 22.Martin JS, Hébert M, Ledoux É, Gaudreault M, Laberge L. Relationship of Chronotype to Sleep, Light Exposure, and Work-Related Fatigue in Student Workers. Chronobiology International. 2012/04/01 2012;29(3):295–304. doi:10.3109/07420528.2011.653656 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/07420528.2011.653656&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22390242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 23. 23.Roenneberg T, Wirz-Justice A, Merrow M. Life between Clocks: Daily Temporal Patterns of Human Chronotypes. Journal of Biological Rhythms. 2003/02/01 2003;18(1):80–90. doi:10.1177/0748730402239679 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730402239679&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12568247&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180630100009&link_type=ISI) 24. 24.Burns AC, Saxena R, Vetter C, Phillips AJK, Lane JM, Cain SW. Time spent in outdoor light is associated with mood, sleep, and circadian rhythm-related outcomes: A cross-sectional and longitudinal study in over 400,000 UK Biobank participants. Journal of Affective Disorders. 2021/12/01/ 2021;295:347–352. doi:[https://doi.org/10.1016/j.jad.2021.08.056](https://doi.org/10.1016/j.jad.2021.08.056) 25. 25.Sudlow C, Gallacher J, Allen N, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLOS Medicine. 2015;12(3):e1001779. 26. 26.Collins R. What makes UK Biobank special? The Lancet. 2012;9822(379):1173–1174. 27. 27.Bycroft C, Freeman C, Petkova D, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562(7726):203–209. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 28. 28.Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. The American journal of human genetics. 2007;81(3):559–575. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 29. 29.Loh P-R, Kichaev G, Gazal S, Schoech AP, Price AL. Mixed-model association for biobank-scale datasets. Nature Genetics. 2018/07/01 2018;50(7):906–908. doi:10.1038/s41588-018-0144-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0144-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29892013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 30. 30.de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: generalized gene-set analysis of GWAS data. PLoS computational biology. 2015;11(4):e1004219. 31. 31.Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nature Communications. 2017/11/28 2017;8(1):1826. doi:10.1038/s41467-017-01261-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-01261-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 32. 32.Romanov RA, Zeisel A, Bakker J, et al. Molecular interrogation of hypothalamic organization reveals distinct dopamine neuronal subtypes. Nature Neuroscience. 2017/02/01 2017;20(2):176–188. doi:10.1038/nn.4462 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.4462&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27991900&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 33. 33.Zheng Z, Huang D, Wang J, et al. QTLbase: an integrative resource for quantitative trait loci across multiple human molecular phenotypes. Nucleic acids research. 2020;48(D1):D983–D991. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkz888&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 34. 34.Ratnapriya R, Sosina OA, Starostik MR, et al. Retinal transcriptome and eQTL analyses identify genes associated with age-related macular degeneration. Nature Genetics. 2019/04/01 2019;51(4):606–610. doi:10.1038/s41588-019-0351-9 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0351-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30742112&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 35. 35.Grotzinger AD, Rhemtulla M, de Vlaming R, et al. Genomic structural equation modelling provides insights into the multivariate genetic architecture of complex traits. Nature Human Behaviour. 2019/05/01 2019;3(5):513–525. doi:10.1038/s41562-019-0566-x [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41562-019-0566-x&link_type=DOI) 36. 36.Wirz-Justice A. Biological rhythm disturbances in mood disorders. International Clinical Psychopharmacology. 2006;21:S11–S15. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/01.yic.0000195660.37267.cf&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=AMBIGUOUS 16&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 37. 37.Agorastos A, Olff M. Traumatic stress and the circadian system: neurobiology, timing and treatment of posttraumatic chronodisruption. European Journal of Psychotraumatology. 2020;11(1):1833644. 38. 38.Arendt J. Melatonin, circadian rhythms, and sleep. Mass Medical Soc; 2000. p. 1114–1116. 39. 39.Maury E, Ramsey KM, Bass J. Circadian rhythms and metabolic syndrome: from experimental genetics to human disease. Circulation research. 2010;106(3):447–462. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6OToiMTA2LzMvNDQ3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMzAvMjAyMi4wNi4yOS4yMjI3NzA3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 40. 40.Bulik-Sullivan BK, Loh P-R, Finucane HK, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nature genetics. 2015;47(3):291–295. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 41. 41.Stone JE, McGlashan EM, Quin N, et al. The Role of Light Sensitivity and Intrinsic Circadian Period in Predicting Individual Circadian Timing. J Biol Rhythms. Dec 2020;35(6):628–640. doi:10.1177/0748730420962598 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730420962598&link_type=DOI) 42. 42.Goulet G, Mongrain V, Desrosiers C, Paquet J, Dumont M. Daily light exposure in morning-type and evening-type individuals. Journal of Biological Rhythms. 2007;22(2):151–158. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730406297780&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17440216&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245246600006&link_type=ISI) 43. 43.Martinez-Nicolas A, Martinez-Madrid MJ, Almaida-Pagan PF, Bonmati-Carrion M-A, Madrid JA, Rol MA. Assessing Chronotypes by Ambulatory Circadian Monitoring. Front Physiol. 2019;10:1396–1396. doi:10.3389/fphys.2019.01396 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fphys.2019.01396&link_type=DOI) 44. 44.Cain SW, McGlashan EM, Vidafar P, et al. Evening home lighting adversely impacts the circadian system and sleep. Scientific reports. 2020;10(1):1–10. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-61712-w&link_type=DOI) 45. 45.Wright Kenneth P, McHill Andrew W, Birks Brian R, Griffin Brandon R, Rusterholz T, Chinoy Evan D. Entrainment of the Human Circadian Clock to the Natural Light-Dark Cycle. Current Biology. 2013/08/19/ 2013;23(16):1554–1558. doi:[https://doi.org/10.1016/j.cub.2013.06.039](https://doi.org/10.1016/j.cub.2013.06.039) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cub.2013.06.039&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23910656&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 46. 46.Pradeepa MM. Causal role of histone acetylations in enhancer function. Transcription. 2017;8(1):40–47. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/21541264.2016.1253529&link_type=DOI) 47. 47.Moore LD, L. T, Fan G. DNA Methylation and Its Basic Function. Neuropsychopharmacology. 2013/01/01 2013;38(1):23–38. doi:10.1038/npp.2012.112 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2012.112&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22781841&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312099500003&link_type=ISI) 48. 48.Peter CJ, Saito A, Hasegawa Y, et al. In vivo epigenetic editing of Sema6a promoter reverses transcallosal dysconnectivity caused by C11orf46/Arl14ep risk gene. Nat Commun. Sep 11 2019;10(1):4112. doi:10.1038/s41467-019-12013-y [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-12013-y&link_type=DOI) 49. 49.Fischbach BV, Trout KL, Lewis J, Luis CA, Sika M. WAGR syndrome: a clinical review of 54 cases. Pediatrics. 2005;116(4):984–988. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBlZGlhdHJpY3MiO3M6NToicmVzaWQiO3M6OToiMTE2LzQvOTg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMzAvMjAyMi4wNi4yOS4yMjI3NzA3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 50. 50.Hysi PG, Choquet H, Khawaja AP, et al. Meta-analysis of 542,934 subjects of European ancestry identifies new genes and mechanisms predisposing to refractive error and myopia. Nature Genetics. 2020/04/01 2020;52(4):401–407. doi:10.1038/s41588-020-0599-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-0599-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32231278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 51. 51.Najmabadi H, Hu H, Garshasbi M, et al. Deep sequencing reveals 50 novel genes for recessive cognitive disorders. Nature. 2011;478(7367):57–63. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10423&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21937992&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295575400034&link_type=ISI) 52. 52.Alazami AM, Patel N, Shamseldin HE, et al. Accelerating novel candidate gene discovery in neurogenetic disorders via whole-exome sequencing of prescreened multiplex consanguineous families. Cell reports. 2015;10(2):148–161. 53. 53.Hornberger MR, Dütting D, Ciossek T, et al. Modulation of EphA Receptor Function by Coexpressed EphrinA Ligands on Retinal Ganglion Cell Axons. Neuron. 1999/04/01/ 1999;22(4):731–742. doi:[https://doi.org/10.1016/S0896-6273(00)80732-1](https://doi.org/10.1016/S0896-6273(00)80732-1) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0896-6273(00)80732-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10230793&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000079976200016&link_type=ISI) 54. 54.Erskine L, Herrera E. The retinal ganglion cell axon’s journey: insights into molecular mechanisms of axon guidance. Developmental biology. 2007;308(1):1–14. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ydbio.2007.05.013&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17560562&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248589700001&link_type=ISI) 55. 55.Gompf HS, Fuller PM, Hattar S, Saper CB, Lu J. Impaired circadian photosensitivity in mice lacking glutamate transmission from retinal melanopsin cells. Journal of biological rhythms. 2015;30(1):35–41. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730414561545&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25512304&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 56. 56.Golombek DA, Rosenstein RE. Physiology of circadian entrainment. Physiological reviews. 2010;90(3):1063–1102. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.00009.2009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20664079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280404600007&link_type=ISI) 57. 57.Park HT, Kang EK, Bae KW. Light regulates Homer mRNA expression in the rat suprachiasmatic nucleus. Molecular brain research. 1997;52(2):318–322. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0169-328X(97)00292-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9495554&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071730800017&link_type=ISI) 58. 58.Brandstätter JH, Dick O, Boeckers TM. The postsynaptic scaffold proteins ProSAP1/Shank2 and Homer1 are associated with glutamate receptor complexes at rat retinal synapses. Journal of Comparative Neurology. 2004;475(4):551–563. doi:[https://doi.org/10.1002/cne.20194](https://doi.org/10.1002/cne.20194) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cne.20194&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15236236&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222453300007&link_type=ISI) 59. 59.Xiao B, Tu JC, Petralia RS, et al. Homer regulates the association of group 1 metabotropic glutamate receptors with multivalent complexes of homer-related, synaptic proteins. Neuron. 1998;21(4):707–716. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0896-6273(00)80588-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9808458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000076697300015&link_type=ISI) 60. 60.Pilz LK, Levandovski R, Oliveira MA, Hidalgo MP, Roenneberg T. Sleep and light exposure across different levels of urbanisation in Brazilian communities. Scientific reports. 2018;8(1):1–11. 61. 61.Lewy AJ, Bauer VK, Ahmed S, et al. The Human Phase Response Curve (Prc) to Melatonin is About 12 Hours out of Phase with the Prc to Light. Chronobiology International. 1998/01/01 1998;15(1):71–83. doi:10.3109/07420529808998671 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/07420529808998671&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9493716&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071928100008&link_type=ISI) 62. 62.Arnau-Soler A, Adams MJ, Generation S, Major Depressive Disorder Working Group of the Psychiatric Genomics C, Hayward C, Thomson PA. Genome-wide interaction study of a proxy for stress-sensitivity and its prediction of major depressive disorder. PLOS ONE. 2018;13(12):e0209160. doi:10.1371/journal.pone.0209160 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0209160&link_type=DOI) 63. 63.Badia P, Myers B, Boecker M, Culpepper J, Harsh J. Bright light effects on body temperature, alertness, EEG and behavior. Physiology & behavior. 1991;50(3):583–588. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0031-9384(91)90549-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1801013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GK17100019&link_type=ISI) 64. 64.Jung CM, Khalsa SBS, Scheer FA, et al. Acute effects of bright light exposure on cortisol levels. Journal of biological rhythms. 2010;25(3):208–216. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730410368413&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20484692&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277861500006&link_type=ISI) 65. 65.Chellappa SL, Steiner R, Oelhafen P, et al. Acute exposure to evening blue-enriched light impacts on human sleep. J Sleep Res. Oct 2013;22(5):573–80. doi:10.1111/jsr.12050 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jsr.12050&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23509952&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000324335300012&link_type=ISI) 66. 66.Sandahl H, Baandrup L, Vindbjerg E, Jennum P, Carlsson J. Social zeitgebers and circadian dysrhythmia are associated with severity of symptoms of PTSD and depression in trauma-affected refugees. European Archives of Psychiatry and Clinical Neuroscience. 2021;271(7):1319–1329. 67. 67.Agorastos A, Linthorst AC. Potential pleiotropic beneficial effects of adjuvant melatonergic treatment in posttraumatic stress disorder. Journal of pineal research. 2016;61(1):3–26. 68. 68.Youngstedt SD, Kline CE, Reynolds AM, et al. Bright Light Treatment of Combat-related PTSD: A Randomized Controlled Trial. Military Medicine. 2022;187(3-4):e435–e444. doi:10.1093/milmed/usab014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/milmed/usab014&link_type=DOI) 69. 69.Duffy JF, Czeisler CA. Age-related change in the relationship between circadian period, circadian phase, and diurnal preference in humans. Neuroscience letters. 2002;318(3):117–120. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0304-3940(01)02427-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11803113&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173643600002&link_type=ISI) 70. 70.Czeisler CA, Dumont M, Duffy JF, et al. Association of sleep-wake habits in older people with changes in output of circadian pacemaker. Lancet (London, England). Oct 17 1992;340(8825):933–6. doi:10.1016/0140-6736(92)92817-y [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0140-6736(92)92817-Y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1357348&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JT82300003&link_type=ISI) 71. 71.Crowley SJ, Cain SW, Burns AC, Acebo C, Carskadon MA. Increased Sensitivity of the Circadian System to Light in Early/Mid-Puberty. The Journal of Clinical Endocrinology & Metabolism. 2015;100(11):4067–4073. doi:10.1210/jc.2015-2775 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2015-2775&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26301944&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) 72. 72.Herljevic M, Middleton B, Thapan K, Skene DJ. Light-induced melatonin suppression: age-related reduction in response to short wavelength light. Experimental gerontology. Mar 2005;40(3):237–42. doi:10.1016/j.exger.2004.12.001 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.exger.2004.12.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15763401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228158500014&link_type=ISI) 73. 73.Sletten TL, Revell VL, Middleton B, Lederle KA, Skene DJ. Age-related changes in acute and phase-advancing responses to monochromatic light. J Biol Rhythms. Feb 2009;24(1):73–84. doi:10.1177/0748730408328973 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730408328973&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19227580&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262736600008&link_type=ISI) 74. 74.Scheuermaier K, Laffan AM, Duffy JF. Light Exposure Patterns in Healthy Older and Young Adults. Journal of Biological Rhythms. 2010/04/01 2010;25(2):113–122. doi:10.1177/0748730410361916 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0748730410361916&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20348462&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F30%2F2022.06.29.22277078.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276010100005&link_type=ISI) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif