Identifying pathways modulating sleep duration: from genomics to transcriptomics

Recognizing that insights into the modulation of sleep duration can emerge by exploring the functional relationships among genes, we used this strategy to explore the genome-wide association results for this trait. We detected two major signalling pathways (ion channels and the ERBB signalling family of tyrosine kinases) that could be replicated across independent GWA studies meta-analyses. To investigate the significance of these pathways for sleep modulation, we performed transcriptome analyses of short sleeping flies’ heads (knockdown for the ABCC9 gene homolog; dSur). We found significant alterations in gene-expression in the short sleeping knockdowns versus controls flies, which correspond to pathways associated with sleep duration in our human studies. Most notably, the expression of Rho and EGFR (members of the ERBB signalling pathway) genes was down- and up-regulated, respectively, consistently with the established role of these genes for sleep consolidation in Drosophila. Using a disease multifactorial interaction network, we showed that many of the genes of the pathways indicated to be relevant for sleep duration had functional evidence of their involvement with sleep regulation, circadian rhythms, insulin secretion, gluconeogenesis and lipogenesis.

flies and mammals 9 and support the notion that studies in Drosophila give insights into the molecular basis of sleep in more complex organisms. Therefore, we used a systems biology approach to compare the results of gene expression in dSur KD Drosophila with significant pathways from our GWAS (current and past) 5 meta-analyses results (for study design see Supplementary Material and Methods, Fig. 1). Remarkably, the two main pathways associated with sleep duration in humans were ranking among the top significant pathways in our gene expression analysis in Drosophila. We also confirmed that the expression of several genes, reported to be involved in sleep regulation in model organisms, was down-regulated in the short-sleeping flies. The molecular pathways of these genes are also associated with human sleep duration. Moreover, several genes of the top raking GWAS and Drosophila transcriptome pathways are known to be involved in the development of metabolic dysfunction, as shown in Figure 1.

Materials and Methods
Study Populations. Data from 3 GWAS consisting of 1564 phenotyped and genotyped individuals (all of European ancestry and informative for sleep duration) were meta-analysed in the discovery phase (Table 1) Table S1 and Supplementary Text). The data of the GEC cohort was obtained through an online-based platform and biological samples were obtained by mailing participants self-collecting saliva sample-kits (Oragene; DNA Genotek, Inc., Ottawa, Ontario, Canada) for DNA extraction. The consistence of the method was tested comparing intra-individual blood and saliva DNA of a sub-sample of individuals ( Supplementary Fig. S2). To replicate the results of our pathway analyses, we made use of independent meta-analysis results on sleep duration GWAS that we had published 5 . A description of the 7 cohorts included in that study (KORA -Cooperative health research in the Region of Augsburg, KORCULA -the Korcula study in Croatia, ERF -Erasmus Rucphen Family, EGCUT -Estonian Genome Center, University of Tartu, MICROS -the Microisolates in South Tyrol Study, NESDA -the Netherlands Study of Depression and Anxiety and the ORCADES Figure 1. Human gene vs. disease multifactorial interaction network. The network shows interactions between genes from significant pathways (ERBB signaling family of tyrosine kinases and ion channels) identified based on the GWAS datasets (Meta3, green nodes; Meta7, light-blue nodes; overlapping genes between Meta3 and Meta7, orange nodes) and the Drosophila transcriptome GSEA. Only human homologs of the respective Drosophila genes with altered expression in the dSur KD flies in relation to controls (Fig. 2) were included. Relationships between the respective genes and biological processes (diabetes, carbohydrate and lipid metabolism, orange; cardiovascular diseases, beige), protein complexes and phenotypes (abnormal sleep and circadian behaviour, blue) are also shown (decreased activity/expression, grey edges with cross bar; increased activity/expression, green colored arrows; modulated activity/expression, green edges with open diamond). Red edges indicate gene expression for flies pooled every 3 h of the 24 hours period (decreasing expression, cross bar; increased expression in relation to wild-type controls, arrows; ratio RNAi/wt). Protein-protein interactions are displayed as black edges; interlinking genes are shown as grey nodes. A list of interactions with literature references is available in the Supplementary Table S8. Genotyping and Imputation. Genotyping of BREC, EGCUT and GEC were accomplished all with the same platform (Illumina Human OmniExpress 700 K) at the Helmholtz Centre Munich or Estonian Biocentre Genotyping Facilities. The overall quality-control criteria excluded individuals with low call rates or deviation from Hardy-Weinberg equilibrium, excess heterozygosity, gender mismatch and bad clustering in the stratification analysis ( Supplementary Fig. S3). Non-genotyped SNPs were imputed by using the 1000 Genomes reference panel and the IMPUTE2 12 software. Cohorts included in the meta-analysed published GWAS 5 were genotyped on a variety of platforms. The quality control and imputation details can be found in Supplementary Materials and Methods, Figure 4. Differentially expressed genes in dSur KD flies vs. controls. Heat map representing a selection of deregulated transcripts homologs to human genes (indicated in brackets) associated with sleep duration in the meta-analysis results used for the GSEA in (a) "Pooled" (Drosophila pooled every 3 h of the 24 hours period) and (b) "Night" (3 h into the night) conditions. A color-coded scale for the normalized expression values is used as follows: yellow and blue represent high and low expression levels in dSur KD with respect to control, respectively. The expression level of each transcript was calculated as the log2 (dSur KD/control), and the complete lists of differentially expressed genes identified by LIMMA software are provided in the Supplementary Table S7. Statistical Analysis. Genotypes consisting of both directly typed and imputed SNPs (N ≈ 4.9-9 Mio, as specified in Supplementary Table S3) entered the GWA analyses. To avoid over-inflation of test statistics due to population structure or relatedness, we applied genomic controls for the independent studies and meta-analysis (Supplementary Materials and Methods, Supplementary Fig. S4 and Supplementary Table S3). Linear regression (PLINK) for associations with normalized sleep duration was performed under an additive model, with SNP allele dosage as predictor and with age, age2, gender, normalized MSFsc (Midsleep on free days, corrected for sleep-depth accumulated during the workweek), season of assessment (dichotomized based on time of the year; day-light savings time -DST or standard zone time assessments) and BMI as covariates ( Supplementary  Fig. S5). A fixed-effects meta-analysis was conducted using the inverse-variance-weighted method in PLINK ( Supplementary Fig. S6). All SNPs with low MAF (<0.01) and low imputation quality (Rsq/proper_info <0.3) were dropped from the meta-analysis. Corresponding to Bonferroni adjustment for one million independent tests 13 , we specified a threshold of P < 5 × 10 −8 for genome-wide significance. To assess the number of independent loci associated with sleep duration at the P < 10 −5 (Table 1 and Supplementary Table S4), correlated SNPs were grouped using a LD-based result clumping procedure (PLINK).
Gene set enrichment analysis for GWAS. To identify pathways/gene sets associated with sleep duration, we used an improved gene set enrichment analysis for GWAS (i-GSEA4GWAS) 8 . This method overcomes the limitations of single SNP analysis in GWAS, using a gene set enrichment approach to yield pathways/gene sets associated with the trait of interest. Genes were mapped within 100 kb up-and downstream of the SNP (threshold P < 0.05) region. We searched for gene sets relevant for canonical pathways from a variety of online resources and GO biological processes, molecular function and cellular components. Pathways/gene sets with FDR < 0.25 were regarded as to be possibly associated with traits, whereas those with FDR < 0.05 were regarded as high confidence or with statistical significance (Supplementary Tables S5 and S6).

Gene versus diseases network analyses.
To integrate information of genes/proteins belonging to the most significant pathways with their validated relevance in several biological processes, we generated a multifactorial interaction network (Fig. 1). We used CIDeR, a publicly available database that integrates interactions between heterogeneous factors associated with human diseases 14 . All interactions shown in Figure 1 are based on experimental evidence from scientific literature (Supplementary Table S8: CIDeR interactions 14 ). Our analyses covered disease-related interactions from sleep disorders, as well as the metabolic dysfunction (obesity and type I and II diabetes).
Fly strains and Behavioural Assays. Flies were raised at 23 °C on yeast cornmeal agar food. An elav-Gal4 strain was crossed to the RNAi UAS-Sur line (V104241) from the Vienna Drosophila RNAi Centre. The KD effect was quantified and found to be ~40% (Supplementary Fig. S7 and Supplementary Table S9). Control experiments were performed with V10424 and elav-Gal4 flies. Drosophila Activity Monitoring System devices (DAMS) from Trikinetics (Waltham, MA) were used to collect locomotor activity and sleep data in 1-min bins. Two-day-old males (32 flies per genotype) were acclimatised in 12 h:12 h light/dark cycles at 23 °C for two days, and activity was recorded under these conditions for 3 further days. Sleep deprivation was calculated from the locomotor activity data by using a Microsoft Excel script in which sleep was defined as 5 min of consecutive inactivity of the   flies 15 . All animal procedures were approved by the local ethical committee at University of Padova and followed international guidelines.

Microarray labelling and hybridization. Gene expression profiling was carried out on dSur KD
(elav-Gal4/RNAi-104241 knockdown) and parental control strain (RNAi-104241) fruit flies using the Drosophila 1.0 custom platform (Agilent Technologies) in two experimental conditions: "Pooled" (pooled samples collected every three hours over a 24 h period, including also ZT15) and "Night" (collected in the dark phase at ZT15 as representative of the time-point of pronounced sleep shortening in the mutant flies). Flies were collected from the incubator during darkness at the night time. Red light was used for a few seconds only to transfer the flies into vials prior to freezing. For each condition, "Pooled" and "Night", total RNA was extracted from the head of 30 flies for each genotype, each biological replicate, and each time point. Four biological replicates were analysed for dSur KD and control samples respectively for each time point for a total of 16 microarray experiments. Details on total RNA isolation procedure, microarray design and microarray labelling are presented in Supplementary Materials and Methods. Gene expression data are available in the GEO database with the accession number: GSE52764.

Statistical analysis of gene expression data.
Inter-array normalization of expression levels was performed with quantile method 16  Functional enrichment analysis of the transcriptomic's results. Functional enrichment analysis of differentially expressed genes in the dSur KD relative to expression observed in the control flies in the "Night" (flies collected at ZT15) and "Pooled" (flies collected every 3 hours over a 24 hours period) experiments was performed using the Graphite web tool 19 , which is based on Reactome pathway database 20 . The hypergeometric test (Fisher Exact test), which estimates the chance probability of observing a given number of genes from a pathway among the selected differentially expressed genes, was applied. Raw p-values were adjusted using Benjamini-Hockberg method.

Validation of relative gene expression by quantitative RT-PCR. Quantitative RT-PCR was used
to validate the expression values of four differentially expressed genes obtained from microarray experiments (Fig. 2). Experimental details can be found in the Supplementary

Results
Here, we report for the first time the GWAS data of our current meta-analysis on sleep duration (GEC, BREC and EGCUT; see Supplementary Figs S3-S6, S8 and Supplementary Tables S3 and S4). Following the traditional genomic analyses, no population stratification was observed after genomic control in any of the three investigated cohorts (inflation factors λ, 1.00 to 1.03 for the independent studies and 0.99 for the overall genome-wide meta-analysis), as shown on the quantile-quantile plots ( Supplementary Fig. S4). Meta-analysis of all 3 studies (N = 1564) did not yield a genome-wide significant signal ( Supplementary Fig. S6); while independent genome-wide significant hits were observed in the GEC cohort ( Supplementary Fig. S5). A summary of the best ranking loci based on P values and number of correlated SNPs (see Materials and Methods for details) is listed in Table 1 (see also Supplementary Table S4 and Supplementary Fig. S8). We had shown earlier 5 that the effect of ABCC9 gene variant on epidemiological variation in human sleep duration is influenced by inter-individual differences both in chronotype (phase of entrainment) and seasonal differences in entrainment of the biological clock. We have therefore adjusted our analysis for season of assessment and chronotype. Studies not considering these confounding factors, may not be able 21 to replicate our findings.
To get insights into the modulation of sleep duration by exploring the functional relationships among genes, we further analysed our current and past 5 GWAS meta-analyses results to identify overrepresented gene-sets or pathways associated with sleep duration. Using a gene set enrichment analysis (GSEA) for the GWAS meta-analyses, we discovered a number of significant pathways (Table 2), including Circadian Exercise, a circadian pathway. This pathway was significant in both GSEA for the GWAS (Meta 7) and transcriptomic data ("Pooled" Drosophila), as shown in Table 2. Only two types of signalling pathways (ion channel genes and the ERBB signalling family of tyrosine kinases) associated with sleep duration in both GWAS meta-analyses and in the transcriptomic experiments; Table 2). We succeeded thereby to replicate these pathways, proving that this method can increase power for uncovering genetic factors relevant for sleep duration. As to be expected, only a few genes in these pathways were overlapping across studies (Fig. 1). To get insights into the functional relevance of these gene sets/pathways for sleep duration, we combined these findings to the results of transcriptome analysis from a short sleeping Drosophila model, knockdown for the ABCC9 gene homolog (dSur), versus control flies (Figs 1 and 2, experimental validation is presented in Supplementary Fig. S9). We performed the gene expression profiling of "Pooled" and "Night" dSur KD since in these flies the quantity of sleep was affected both during the 24 h (22% of reduction) and at ZT15 (36% of reduction) as shown in Supplementary Fig. S9 LIMMA software with default settings (FDR ≤ 5%), we identified 4941 differentially expressed genes (DEGs), of which 2220 were up-regulated (45%) and 2721 were down-regulated (55%) and, 2885 DEG of which 1343 were up-regulated (47%) and 1542 were down-regulated (53%) in RNAi "Pooled" and "Night" samples with respect to controls (Supplementary Table S7). Functional enrichment analysis of DEGs in the dSur KD relative to expression observed in the control flies in the "Night" and "Pooled" experiments were conducted by using the Graphite web tool. The hypergeometric test (Fisher Exact test) was used to rank biological significance of differentially expressed genes in dSur KD versus control flies in the "Night" and "Pooled" conditions. A list of all significant pathways can be found in the Supplementary Tables S5 and 6. Using this strategy, we confirmed that pathways replicated across GWAS meta-analyses were highly significant in the transcriptomic experiments ( Table 2). Pathways that were significant only for one of the GWAS meta-analyses (in processes involving cell death, immuneresponse, signal transduction and metabolic processes), but had several corresponding significant pathways in the Drosophila, gene expression profiles comparisons are also shown in Table 2.
To get a better idea of the network interactions among genes belonging to the replicated pathways in the human studies, we focused specifically on their respective Drosophila gene homologs (Table 2, Fig. 2). Here, we performed a multifactorial interaction networks analysis (Fig. 1) to connect these datasets to: (i) our data on differential gene expression in the dSur KD vs. control flies (Fig. 2), (ii) information from the literature on gene-gene, and gene-protein interactions, and, (iii) published information on the relevance of these genes in the development of human diseases (Fig. 1, Supplementary Table S8). This approach yielded surprising evidences of the relevance of the GWAS associated pathways, which were significantly different in the Drosophila microarray experiments of the KD vs. control flies, i.e., the ERBB signalling pathway ( Table 2 and Figs 2 and 3). Additionally, potassium channel genes from top ranking GWAS pathways were down regulated in the Drosophila short sleeping vs. control flies (Figs 1 and 2). Other top ranking pathways in humans and flies were related to insulin signalling and metabolism, with a down regulation of the genes involved in this pathway in the dSur KD flies (Table 2, Fig. 2 and Supplementary Fig. S10). Differential expression of genes Egfr (Epidermal growth factor receptor), InaC (inactivation no afterpotential C), Hk (Hyperkinetic), and Sh-isoF (Shaker isoform F) was validated by qRT-PCR in both of the experimental conditions (night and pooled flies), as presented in Supplementary Materials and Methods, Supplementary Fig. S11 and Supplementary Table S9. Several other genes associated with T2D, involved in insulin secretion and lipid metabolism had differential expression in the dSur KD vs. control flies and were nominally associated with sleep duration in humans (Figs 1 and 2).

Discussion
Consistent with our previously published GWAS meta-analyses for sleep duration, the top ranking pathways identified both in human and Drosophila studies contained mostly ion channel genes, which are involved in various aspects of neural functions, immune responses, signal transduction (Figs 1 and 3, Table 2). Voltage-gated potassium channels (e.g., KCNA6-Drosophila shaker-related subfamily, member 6) identified in the top ranking pathways of our current meta-analysis in humans, was slightly down regulated in the dSur KD vs. control flies ( Table 2, Supplementary Fig. S11). Voltage-gated potassium channels genes were identified in Drosophila, and studies of genetic inactivation of their closest homologs in mammals proved their involvement in sleep duration also in mice 22 . However, in mammals there is a larger variety of genes encoding these channels, for instance, there is one Shaker gene in Drosophila melanogaster, but at least 16 genes code for α-subunits of voltage-dependent potassium channels in mammals 22 . Additional channels, such as KCNQ (potassium voltage-gated channel subfamily Q), which evolutionarily displays conserved electrophysiology and pharmacology between Drosophila and mammals 23 , had increased expression in the dSur KD flies and 3 of the 5 KCNQ human genes appeared in the best ranking GWAS pathways (Table 2 and Fig. 2). Functional evidence of the relevance of these potassium voltage channels for sleep in humans comes from patients with a rare autoimmune disorder (Morvan's syndrome; see cider multifactorial interaction network in human diseases, Fig. 1), who have marked sleeplessness associated with the presence of autoantibodies against voltage-dependent potassium channels 24 . A potassium ATP dependent channel (KCNJ3, potassium voltage-gated channel subfamily J member 3), which expression is inhibited in bipolar depression and schizophrenia, was down regulated in the Drosophila KD and the gene appeared as part of one of the significant pathways in the human GWAS (Table 2). Also voltage-gated calcium channels CACNA2D1 (auxiliary subunit alpha2delta 1), CACNG5 (auxiliary subunit gamma 5), CACNA1S (subunit alpha1 S), CACNA1E (subunit alpha1 E), CACNA1C (subunit alpha1 C) and CACNA1A (subunit alpha1 A) were among the best ranking gene sets in our pathway analyses in humans. One of them, CACNA1A, which encodes a subunit of Cav2.1, was characterized as a sleep gene in knockout mice 25 , showing reduced expression in our short-sleeping flies (Fig. 2). Related to these calcium channels, mutations in 3 of the genes present in the diagram CACNA1C, CACNA2D1 (calcium voltage-gated channel auxiliary subunit alpha2delta 1) and SCN5A (sodium voltage-gated channel alpha subunit 5) have been shown to cause the Brugada Syndrome, characterized by typical ECG features, ventricular arrhythmias, which may be accompanied by sleep disordered breathing 26 . ERBB signalling pathway and sleep duration. For the first time the ERBB pathway, protein family of tyrosin kinases that consists of four cell surface receptors ErbB1/EGFR/HER1 (epidermal growth factor receptor), ERBB2/HER2 (erb-b2 receptor tyrosine kinase 2), ErbB3/HER3 (erb-b2 receptor tyrosine kinase 3), and ERBB4/ HER4 (erb-b2 receptor tyrosine kinase 4) 27 , has been associated with sleep in humans. This is supported by robust evidence of the role of EGFR signalling pathway modulating sleep in Drosophila 28 and mammals 29 . EGF activates tyrosin kinases acting in intracellular signalling pathways, and in components of intracellular signal-transduction pathways 29 . EGF is expressed in the brain, influences the production of several sleep promoting substances 29 , and its exogenous administration in rabbits enhances non-rapid eye movement sleep. In our experiments, knocking down dSur increased the expression of EGFR (Fig. 2) and shortened sleep duration in these flies in relation to wild Scientific RepoRts | 7: 4555 | DOI:10.1038/s41598-017-04027-7 type controls. In Drosophila, activation of EGF is dependent of Rho (also known as Rhomboid, Rho1, a Rho family GTPases, G protein-coupled receptors), and inhibition of Rho delays sleep consolidation in flies 28 . An increase in EGFR expression in our dSur KD flies ( Fig. 3 and Supplementary Fig. S11) can potentially be influencing sleep consolidation in a similar manner. Rho and EGF are expressed in the pars intercerebralis, a part of the fly brain that is functionally analogous to the hypothalamus in vertebrates, which is a region of the mammalian brain well established to be a regulator of arousal. The role of Rho in the presentation and intramembrane cleavage of TGF-α-like proteins in Drosophila suggests possible similar roles in the presentation of TGF-α-like proteins in vertebrates 30 . A Drosophila Rho homolog in humans (human rhomboid family 1, RHBDF1) was shown to interact with TGF-α family ligands 31 . Stimulation of EGFR signalling in human articular chondrocytes by TGFα results in the activation of RhoA/ROCK (Rho kinase), MEK (MAPK/ERK kinase)/ERK (extracellular-signal-regulated kinase), and PI3K (phosphoinositide 3-kinase) pathways. Our gene set enrichment analyses of genetic and transcriptomic data indicated that EGFR, PI3K and genes of the MAPK/ERK kinase pathway are linked to sleep shortening (Fig. 1). Moreover, Rho down regulation and EGF up-regulation in the dSur KD link EGFR signalling pathway to the ABCC9 K ATP potassium channel. Supporting this evidence, K ATP channels, activated by pinacidil, triggered signalling through Rho kinase in mice 32 . ABCC9 encodes a subunit of the KATP dependent channel which is crucial for cardiovascular and diabetic phenotypes 33 . EGFR inhibition prevents diabetes-induced up-regulation of multiple gene pathways in the mesenteric vasculature 34 . These results suggest that activation of EGFR signalling is a key initiating step that leads to induction of multiple signalling pathways in the development of diabetes-induced vascular dysfunction 34 . Our results support it, showing that knocking down the potassium channel, which leads to a diabetic/cardiomyopathies related phenotype in Drosophila, has downstream effects activating EGF expression. Moreover, they suggest that this is a shortening sleep duration pathway, thereby relating sleep regulation to key components of diabetes and cardiomyopathies canonical pathways. This is supported by the pathway regulating Insulin Secretion, which shows higher significance for flies collected 3 hours into the night (ZT15) than for the pooled flies experiment. Interestingly, mammalian rhodopsin interacts with a Rho family member and activates this small GTPase in a light-dependent fashion 35 . However, no role has been described for transducing light activity through a rhodopsin/small GTPase pathway. Rhodopsin like receptor activity and a GPCRDB class, rhodopsin like, were significant pathways in the current meta-analysis and had similar pathways (GPCR signalling, Table 2) being altered in the transcriptomic study.
The notion that pathways modulating sleep duration could overlap between humans and flies was earlier laid by studies suggesting that the nervous system of both flies and mammals used similar arousal transmitters 9 .
Various neurotransmitters act at cell-surface G protein-coupled receptors (GPCR, pathway which was also significant in our Drosophila gene expression profiles; Table 2) to activate intracellular signal transduction cascades. Activation of the G protein-coupled receptors (such as dopamine receptors) leads to an increase or decrease in adenylate cyclase, modulating cAMP levels. cAMP in turn activates protein kinase A (PKA), which phosphorylates a number of targets, including the transcription factor CREB (cAMP response element binding protein) 9 . The alpha catalytic subunits of PKA (PRKCA in humans, InaC in Drosophila) levels was increased in the dSur KD ( Fig. 2 and Supplementary Fig. S11) and the gene was detected in associated pathways of both sleep duration GWAS meta-analyses (Fig. 1). CREB activity is linked to sleep homeostasis as a CRE reporter gene is up-regulated in response to sleep deprivation and reduced CREB activity results in an elevated sleep rebound 9 . Supporting evidence from the GWAS could be found in our study for CREB5 (cAMP responsive element binding protein 5) and GRM8 (glutamate metabotropic receptor 8), which had the best association signals in the current GWAS meta-analysis (Table 1, Supplementary Fig. S8). GRM8 is also part of the G protein coupled receptor activity pathway ( Table 2), suggesting that although no genome-wide significance could be observed in single SNP analyses, we see the close connection of this to other genes within pathways that significantly associate with sleep duration (Fig. 1). Whereas only CREB has an established role in the transduction pathway, CREB5 increases during sleep and decrease during sleep deprivation in mice ( Supplementary Fig. S12).
Although the relationship between sleep duration and obesity may partially be caused by environmental influences such as voluntary sleep restriction and circadian misalignment, as indicated by previous epidemiological studies 36,37 , the association of sleep duration with pathways related to energy metabolism indicates that genetic factors are central to it. A relation of sleep duration and metabolism was also demonstrated in Drosophila, as cyc 01 (canonical clock protein cycle) mutant flies can sustain periods of waking longer when starved than when sleep deprived, suggesting that the effect of sleep loss may be attenuated during starvation 38 . Gene expression analysis demonstrated that many of the genes with the greatest transcriptional differences between the sleep deprived and starved flies were genes involved with lipid metabolism suggesting a role in sleep homeostasis maintenance. Interestingly, down-regulation of genes involved in lipolysis and up-regulation of genes involved in triglycerides storage resulted in substantially increased sleep rebound following a night of sleep loss. In agreement with these observations, in the dSur KD flies we observed a significant up-regulation of bubblegum (bgm, an acyl-CoA synthetase with homology to the human SLC27a6), heimdall (hll, a long-chain-fatty-acid-CoA ligase bubblegum-like) and Lsd2 (lipid storage droplet-2 39 ) genes. These genes are involved in the control of lipid storage reducing lipolysis and blocking fatty acid release, and their Drosophila knockouts were shown to have sleep rebound alterations following sleep deprivation 38 . Upregulation of these genes in the dSur KD can be thereby a consequence of sleep deprivation in our dSur KD flies.
Our disease multifactorial interaction network shows that many of the genes from the pathways indicated to be relevant for sleep duration (Fig. 1, based on the human GWAS and Drosophila transcriptomics data), had functional evidence of their involvement with sleep regulation, circadian rhythms, insulin secretion, gluconeogenesis and lipogenesis. Moreover, genes that are key for energy metabolism were activated or inhibited in the short sleeping flies, potentially for playing a role on sleep shortening or responding to it.