Kleine Levin syndrome is associated with birth difficulties and genetic variants in the TRANK1 gene loci

Aditya Ambati, Ryan Hillary, Smaranda Leu-Semenescu, Hanna M. Ollila, Ling Lin, Emmanuel During, Neal Farber, Thomas J Rico, Juliette Faraco, Eileen Leary, Andrea Goldstein-Piekarski, Yu-Shu Huang, Fang Han, Yakov Sivan, Michel Lecendreux, Pauline Dodet, Makoto Honda, Natan Gadoth, Sona Nevsimalova, Fabio Pizza, Takashi Kanbayashi, Rosa Peraita Adrados, Guy Leschziner, Rosa Hasan, Francesca Canellas, Kazuhiko Kume, Makrina Daniilidou, Patrice Bourgin, David Rye, José L Vicario, Birgit Högl, Seung Chul Hong, Guiseppe Plazzi, Geert Mayer, Anne Marie Landtblom, Yves Dauvilliers, Isabelle Arnulf, Emmanuel Mignot
doi: https://doi.org/10.1101/2021.01.08.20249006

Abstract

Kleine-Levin Syndrome (KLS) is a rare disorder characterized by severe episodic hypersomnia, with cognitive impairment accompanied by apathy or disinhibition. Pathophysiology is unknown, although imaging studies indicate decreased activity in hypothalamic/thalamic areas during episodes. Familial occurrence is increased, and risk is associated with reports of a difficult birth. We conducted a worldwide case-control genome wide association study in 673 KLS cases collected over 14 years, and ethnically matched 15,341 control individuals. We found a strong genome-wide significant association (OR=1.48,rs71947865,p=8.6×10−9) with 20 single nucleotide polymorphisms encompassing a 35kb region located in the 3’ region of TRANK1 gene, previously associated with bipolar disorder and schizophrenia. Strikingly, KLS cases with TRANK1 rs71947865 variant had significantly increased reports of a difficult birth. As perinatal outcomes have dramatically improved over the last 40 years, we further stratified our sample by birth years and found that recent cases had a significantly reduced TRANK1 rs71947865 association. While theTRANK1 rs71947865 association did not replicate in the entire follow-up sample of 171 KLS cases, the TRANK1 rs71947865 was significantly associated with KLS in the subset follow-up sample of 59 KLS cases who reported birth difficulties (OR=1.54;p=0.01). Genetic liability of KLS as explained by polygenic risk scores was increased (pseudo r2=0.15;p<2.0×10−22 at p=0.5 threshold) in the follow-up sample. Pathway analysis of genetic associations identified enrichment of circadian regulation pathway genes in KLS cases. Our results suggest links between KLS, behavioral rhythmicity, and bipolar disorder, and indicates that the TRANK1 polymorphisms in conjunction with reported birth difficulties may predispose to KLS.

Significance Statement Genetic markers in TRANK1 gene and its vicinity have been weakly associated with bipolar disorder and schizophrenia (10% increased risk). We found that the same polymorphisms are associated with Kleine-Levin Syndrome (50% increased risk), a rare sleep disorder characterized by recurrent episodes of severe hypersomnia and cognitive abnormalities. Response to lithium treatment are suggestive of a pathophysiological overlap between KLS and bipolar disorder. The study also shows that variants in the TRANK1 gene region may predispose to KLS when patients have had a difficult birth, suggesting that TRANK1 gene region modulate newborns’ response to brain injury, with consequences for mental and neurological health in adulthood. Another possibility may be that the polymorphism impact birth and KLS.

Introduction

Kleine Levin Syndrome (KLS) was identified as a unique disease entity almost a century ago (1, 2). It affects primarily male teenagers and is characterized by unpredictable periods (median length ∼10 days) of intense hypersomnia (≥18 hours of sleep per day), occasionally associated with megaphagia and sexual disinhibition. More recent descriptions in large case series have confirmed a stereotypic appearance of abrupt onset episodes in mostly male (66-75%) patients and association of severe hypersomniac episodes with cognitive and behavioral disturbances (mental slowness, confusion, apathy, derealization, altered perception, and disinhibited behavior). Remarkably, there is a complete absence of the symptoms between episodes, and a generally favorable evolution, with spontaneous disappearance of episodes after 1-2 decades (3-6). KLS episodes are dramatic in appearance and patients are fully incapacitated. EEG studies occasionally reveal non-specific EEG slowing during episodes. Further, lithium has been shown to reduce episodes in about one-third of cases (4, 7). Episodes may be associated with flu-like symptoms at the onset, with seasonal winter occurrence suggested by some (8).

Although most cases of KLS are sporadic, familial occurrence (3-8% of relatives) has been reported in studies (5, 9) and there are reports of concordant twin pairs (9, 10)which, considering the rarity of the condition (a few cases per million individuals) (11), suggests genetic or shared environmental effects. Early studies had suggested an HLA association suggestive of an autoimmune mechanism (12), but this was not confirmed in larger samples (3, 8, 13). Further, karyotyping was normal in 112 patients, except for one patient with sporadic KLS who had a duplication in Xp22.31 (13).

This study is the first report of genetic associations in a discovery cohort of 673 KLS cases with attempted replication in 171 cases across the world. Additionally, targeted exome sequencing in 32 cases derived from 16 multiplex families (9) was conducted. A KLS Genome-Wide Association Study (GWAS) found a strong signal in the TRANK1 region (rs71947865), similar to that reported in bipolar disorder and schizophrenia. In these latter disorders the TRANK1 association variably maps to rs9834970 in Caucasians or to rs3732386 in multi-ethnic populations, both SNPs located immediately 3’ of the TRANK1 gene (14-16). Strikingly, the TRANK1 signal was mostly present in KLS patients reporting birth difficulties and decreased over time (defined by birth year), paralleling well-established improved perinatal care. Birth difficulties (and maternal infections) have also long been associated with increased risk of bipolar disorder and schizophrenia (17-29). These results indicate possible gene x environment interactions in KLS, with possibly similar effects in other TRANK1-associated disorders such as bipolar disorder and schizophrenia.

Results

Genome-Wide Association Studies (GWAS) of KLS

A worldwide multiethnic cohort of 673 KLS cases, the largest sample of KLS cases ever to be analyzed, was collected between the years of 2003 and 2017 and compared to 15,341 matched controls (the “discovery cohort”; see Dataset S1 for sample construction, and SI Appendix, fig. S1 for corresponding Principal Component Analysis-PCA plots). Genetic association analyses (genomic inflation λGC = 1.06) found a genome-wide significant effect (Figure 1) with no heterogeneity across ethnic groups, with more than 20 tightly linked single nucleotide polymorphisms (SNPs) encompassing a 35kb region of human chromosome 3 in the 3’ region of the tetratricopeptide repeat and ankyrin repeat-containing protein (TRANK1) gene. The TRANK1 gene was observed to have the highest expression in the fallopian tube, closely followed by the brain-cerebellum (GTEx v8) (see Table 1, Figure 2 for zoom plot and Dataset S2 where all variants in this locus are listed). The most significant association was observed at rs71947865 (OR=1.47, p=8.6×10−9, presence of G versus GGGAGCCA), although the highest OR was observed at rs6789885 (OR=1.53, p=8.9×10−8, T versus C, Dataset S2). This TRANK1 region SNP rs71947865 is in absolute Linkage Disequilibrium (LD) with rs3732386 (r2=1), a SNP associated with Bipolar Disorder and Schizophrenia in multiethnic samples (Dataset S3).

View this table:
Table 1:

Top variants in TRANK1 region associated with KLS in worldwide discovery cohort GWAS meta-analysis in 673 KLS cases and 15341 control individuals and in a follow-up replication cohort of 171 KLS cases and 1956 controls.Effect allele is B, genomic coordinates in hg19/GRCh37.

Figure 1:
Figure 1:

Multiethnic GWAS in 673 KLS cases and 15431 control individuals (A) Manhattan plot of the meta-analysis association results from SNP based GWAS study. The X-axis shows genomic location by chromosome and the y-axis shows -log10 p-values. Red horizontal line indicates genome-wide significant p-value threshold of 5×10−8. Individual loci are annotated. (B) Manhattan plot of gene-based association test association results, SNPs were mapped to 18247 protein coding genes (distance 0). Genome wide significance (red dashed line in the plot) was defined at P = 0.05/18247 = 2.74×10−6.

Figure 2:
Figure 2:

Regional association plot of the GWA significance of chromosome 3 flanked by rs76652637 and rs11706780 in the 3’’ region of TRANK1. Highest significance was observed at rs71947865 (OR=1.47, p=8.6×10−9, presence of G versus GGGAGCCA). The variants are color coded by the magnitude of LD with the rs71947865.

The Linkage Disequilibrium (LD) block contains several eQTL (expression quantitative trait loci) effects affecting numerous genes such as TRANK1, LRRFIP2, MLH1, GOLGA4, DCLK3, STAC and long non-coding RNA RP11-640L9.2. The SNP rs11129735 (KLS discovery GWAS p= 9.9×10−6; LD r2=0.4 with our top GWAS significant variants) had the highest coloc (CLPP) probability (>0.95 causal probability) using eCAVIAR, with the expression of LRRFIP2 in the thyroid, DCLK3 in skin and brain frontal cortex, C3orf35 in testis, MLH1 in breast tissue, TRANK1 in the spinal cord and EPM2AIP1 in the hypothalamus (GTEx v7). Querying the GTEx v8 database (30), we found our lead SNP rs71947865 variant to be a strong eQTL for LRRFIP2 gene expression in the thyroid (p=8.2×10−20) and TRANK1 expression in breast-mammary tissue and whole blood. rs71947865 also had multiple effects on LRRFIP2 expression in various tissues (p<10−6, GTEx portal v8). LRRFIP2 together with MLH1 and GOLGA4 are highly expressed in placenta and differentially expressed in placentae from complicated compared to normal pregnancies (31, 32). The tissue expression for most of these genes is widespread in the brain and poorly characterized. For the reason of convenience, we further called this effect “TRANK1-region polymorphisms” or “TRANK1 effect”, although data to date does not allow genetic attribution of the observed effects to secondary changes in TRANK1 gene expression.

Association of symptoms with polymorphisms

Clinical data collected from 592 of 673 KLS patients from the discovery cohort were next scrutinized for genotype effects. Our a priori hypothesis for testing was to explore whether TRANK1 polymorphisms were associated with either more depressive symptoms, considering lithium responsiveness in the syndrome, and/or a differential effect on birth difficulties, a robust finding in our prior study (3). Other variables were compared as exploratory analyses. TRANK1 rs71947865 (G) genotype was most significantly associated with reports of birth difficulties (OR=1.65; p<0.005 Fisher’s exact test of allelic counts) and developmental delays. In addition, patients with TRANK1 rs71947865 (G) have more severe cognitive symptoms, a longer duration of 1st episode, a more frequent family history of KLS, an older age of onset, and more behavioral disinhibition and altered speech during episodes, generally suggesting more severe disease (Table 2, full list of variables tested in Dataset S4). Types of birth difficulties are summarized in Dataset S5. In a further analysis of the association of clinical variables with the genotype dose of the TRANK1 deletion, we adjusted for reports of birth difficulties. We noted that cognitive signs, age of onset and longer 1st episode duration were significantly (p<0.05) increased in KLS cases with TRANK1 deletion allele, while there was a trend to association of TRANK1 deletion allele with current depression in KLS cases (Table 3).

View this table:
Table 2:

TRANK1 deletion (rs71947865) and its univariate association with KLS clinical variables. We coded the minor allele additively and used linear regression for continuous variables and logistic regression for categorical variables : Model1 - the minor allele genotype was fit as a function of the clinical variable, Model2 – was similar to Model1 except for including Ethnicity, Gender and Missing data proportion as three covariates. A Fisher’s exact test of allelic counts was further performed for categorical variables. In these analyses 592 KLS cases were included (Frequency % | mean ± sem was computed). BB-allele G/G, AB- GGGAGCCA/G, AA - GGGAGCCA/ GGGAGCCA, *Minor Allele G (deletion) was the effect allele.

View this table:
Table 3:

Association of clinical variables with genotype dose of TRANK1 deletion (rs71947865) after adjusting for medical history of birth difficulties, gender, referrer/race & missing data proportion in 592 KLS cases (linear regression for continuous variables and logistic regression for categorical variables). (Frequency % | mean ± sem was computed). BB- allele G, AB- GGGAGCCA/G, AA - GGGAGCCA/ GGGAGCCA, *Minor Allele G (deletion) was the effect allele.

View this table:
Table 4:

A gene set enrichment analysis was carried out on the combined discovery and follow-up cohort of KLS cases using the MAGMA module. A total of 5346 sets including curated gene sets, KEGG pathways, Reactome pathways and GO biological processes derived from msigdb6.1 were used. (5% FDR P value was computed). Beta is the direction of association of the pathway/geneset with the KLS case control status.

TRANK1 genotyping was conducted in 195 controls with birth difficulty data available. These controls were from the KLS case control study reported by Arnulf et al. (3), with further addition of controls (friend or spouse of KLS patient) over more recent years. In the control sample, 22 subjects reported birth difficulties and, unlike in the KLS cases we found no significant association of birth difficulties with TRANK1 rs71947865 genotype (p=0.5, OR=1.05). To formally test for an interaction effect of TRANK1 genotype on birth difficulties in KLS versus controls, we used the combined KLS cases-control samples matched by PCA. Whereas both birth difficulties and TRANK1 had a significant effect on KLS risk (p<0.05), the interaction term of TRANK1 genotype and presence of birth difficulties did not reach statistical significance. Note that the sample of controls with clinical (birth difficulty) data available is small and, although it was also used in the GWAS as a control, most of the controls used for GWAS were drawn from the Genetic Epidemiology Research on Aging (GERA)(33) cohort and/or other controls obtained through other studies from the same countries, which come from similar countries and collaborators.

The TRANK1 association is birth year dependent

As indicated above, our analysis revealed a significant association of the TRANK1 variant with birth difficulties in KLS patients but not controls. As perinatal care has improved over the last decades (34-39); https://www.cdc.gov/nchs/data/nvsr/nvsr56/nvsr56_10.pdf; https://mchb.hrsa.gov/chusa14/health-status-behaviors/infants/infant-mortality.html), most notably in the 1990-2000s, we asked if the year of birth was an essential factor correlated to KLS in conjunction with TRANK1 rs71947865 variant. A possible birth year dependent TRANK1 interaction is also supported by the recent largest bipolar disorder GWAS ever published, where it was noteworthy that the strong TRANK1 SNP association found in the discovery sample in bipolar disorder cases did not replicate in a new, presumably more recent, bipolar GWAS replication set (15), a caveat being that birth year data was not made available in this study (15) for all but one replication cohort where the birth year range was 1981-2005. Further, later reported TRANK1 ORs for bipolar disorder decreased in size from earlier studies in Brazil (40) and other countries (41-43) to the most recent metanalysis (15), even independently of the non-significant replication sample of this most recent metanalysis (15). This is illustrated in Dataset S3, where published TRANK1 association results are reported, with the caveat that the same samples are typically used cumulatively across studies.

In a subset of the full cohort (discovery and replication sets) when a valid birth year was available (n=650), we performed TRANK1 SNP association by birth year stratified into three approximately equal-sized groups; 1987 or earlier (n=220, median year of birth: 1978), 1988-1994 (n=202, median year of birth: 1992) and 1995-2000’s (n=228, median year of birth: 1999). Using genetic principle component matched controls for each stratum mainly derived from the GERA cohort (birth year range ∼1928 to 1989), we observed a decrease in TRANK1 SNP rs71947865 association with time (Figure 3A and Dataset S6A). Notably, while individuals born in 1987 or earlier had the strongest association TRANK1 association (p=8.3×10−4; OR=1.41), this effect decreased in KLS cases born 1988-1994 (p=7.48.3×10−3; OR=1.34) to finally disappear in the 1995-2000’s group (p=0.31; OR=1.11). In addition, within KLS cases, a Fisher’s test of allelic counts was performed, birth year bin 1987 or earlier was compared to 1988-1994 and 1995-2000’s in 2×2 contingency table. We observed highest association of TRANK1 SNP rs71947865 (p=0.002) in year bin 1987 or earlier when compared to the minor allele count of TRANK1 SNP rs71947865 in the birth year bin 1995-2000’s (see Dataset S6B).

Figure 3:
Figure 3:

(A)The TRANK1 region SNP associations in KLS cases stratified by birth years, the y-axis represents the odds ratio of the effect allele while the x-axis is the birth year. Inset P value of the association is labelled. (B) The TRANK1 region SNP associations in KLS cases stratified by birth difficulties, the y-axis represents the odds ratio of the effect allele while the x-axis is the presence of birth difficulties.

GWAS Analysis of follow up KLS cases with and without reports of birth difficulties

To replicate these findings, we collected and genotyped 171 new KLS cases collected between 2017-2019, collecting information on birth difficulties whenever possible. In the resulting sample set, birth difficulties data were available for 151 of the 171 cases. In the overall replication sample (n=171), generally born more recently (median year of birth: 1998, Range: 1960 to 2015), TRANK1 (rs71947865) was not significantly associated with KLS (p=0.1, OR=1.1) (Table 1 and Dataset S2). Strikingly, however, when we stratified this sample into cases with versus without birth difficulties, then a significant TRANK1 association was seen only in the first group (p=0.01; OR=1.54 in n=59 cases reporting birth difficulties, Dataset S7). In KLS cases with no birth difficulties, little effect was observed (p=0.4; OR=1.12 in n=88 cases). No birth information was available for n=24 patients. A combined TRANK1 SNP meta-analysis (OR=1.4;p=3.5×10−8) of both discovery and replication cohorts is presented in Table 1

Our prior analyses of the discovery cohort had revealed significant association of TRANK1 variant rs71947865 with birth difficulties, we therefore determined the association of TRANK1 rs71947865 variants in KLS cases with birth difficulties when compared to controls individuals in both discovery and replication cohorts. We observed that TRANK1 rs71947865 variants were significantly associated birth difficulties in KLS cases from the both discovery cohort (p=0.002; OR=1.61, see Dataset S8A) and the follow-up sample of KLS cases (p=0.03; OR=1.53). A combined analysis of KLS cases with and without birth difficulties from combining both discovery and replication cohort compared to genetic PCA matched controls also revealed significantly increased TRANK1 variant rs71947865 frequencies in KLS cases with birth difficulties (Fisher’s test, OR=1.54; p=4.83×10−4, Figure 3B & Dataset S8B) while in KLS cases without birth difficulties, the effect was weaker (Fisher’s test, OR=1.20;P=0.04).

Fine mapping of the KLS signal across ethnicity

GWAS studies in schizophrenia and bipolar disorders have mapped the TRANK1 signal at SNP rs9834970 (r2=0.74 with KLS SNP rs71947865) in Caucasians (15, 42-44), while when a mixed ethnic sample that included Asians and Caucasian was used, the top SNP is rs3732386 (16, 45-47), a signal entirely linked (r2=1) with the KLS associated SNP rs71947865 (see Dataset S3; SI Appendix, fig. S2A-sC).

Since our KLS sample is derived from multiple ethnicities, as an additional validation, we sought to determine TRANK1 loci association KLS cases of European/Caucasian origin only. In the overall European KLS meta-analysis, we observed that the TRANK1 region SNP rs71947865 (OR=1.28; p=1.9×10−4) was retained as the top association (see Dataset S9 and SI Appendix, fig. S2A-C) with no evidence of heterogeneity. Our KLS European/Caucasians only analysis (discovery plus replication cohort) revealed that rs9834970 was among the top associations in the loci (SI Appendix, fig. S2D and Dataset S9). The rs3732386 and rs71947865 were the top associations in the TRANK1 region and were dependent on the existence of associated birth difficulties (SI Appendix, fig. S2). Further, on comparing the associations in the TRANK1 loci in discovery KLS versus combined schizophrenia and bipolar disorder GWAS, there is remarkable overlap of signals with only the magnitude of ORs being 1.3x times increased in KLS (see SI Appendix, fig. S2D).

Replication of Polygenic risk score of KLS

We found significant evidence that the KLS discovery cohort polygenic risk scores (PRS) were indicative of increased genetic liability to KLS phenotype in a follow-up replication cohort (n=171) under various p-value thresholds using PRSice2 (48). The most predictive p-value threshold was 0.5 (pseudo R2 p=2×10−22; Figure 4) using 81,570 variants. P-value thresholds ranging from 0.05 to 0.4 were also indicative of increased genetic liability to KLS (p<1×10−20, mean pseudo R2 =0.14, see Dataset S10). Further, as an additional control, we also used only Caucasian KLS sample discovery meta-analysis (n=299) to build PRS with validation only in Caucasian KLS replication sample (n=160), the PRSs were indicative again of increased genetic liability to KLS (pseudo R2 =0.094) with the best score at p=0.4 threshold (see SI Appendix, fig. S3).

Figure 4:
Figure 4:

Variance explained by PRS (polygenic risk scores) built on the KLS discovery cohort when used to predict on the replication cohort, the y-axis is the % variance explained by the PRS while the x-axis represents the p-value thresholds.

Considering the results above, we further built PRSs for Bipolar disorder (15), Schizophrenia (49), circadian chronotype (50), sleep duration (51), and excessive daytime sleepiness (52) using previously published GWAS studies, and evaluated their predictive performance on our KLS follow up cohort. Shared KLS genetic liability (pseudo R2 =0.06) was found using bipolar disorder and schizophrenia GWAS summary statistics as the training set (see SI Appendix, fig. S3 and Dataset S10). LD score regression analysis estimated the total observed scale heritability h2=0.09 (intercept=1.044) and also revealed significant genetic overlap of KLS with bipolar disorder (rg=0.41;p<0.05) and depressive symptoms (rg=0.32;p<0.05, Dataset S11).

Gene Pathway analysis

Biological processes pathway enrichment analysis was next conducted using the MAGMA (53). The MAGMA module gene-based association analysis was first performed in the combined discovery-replication sample summary statistics mapping SNPs to genes with 1kb window upstream and downstream. This first step entails collapsing mean individual SNP statistics into a per -gene level statistic. Then a gene p-value is obtained by using a known approximation of the sampling distribution as detailed in MAGMA (53). The gene-based p values converted to z-scores are then used to conduct gene-set enrichment analyses on gene sets derived from MsigDBv6.1 (5346 curated gene sets). At 5% FDR (false discovery rate), pathways “Reactome-RORA activates circadian expression”(54) and “Reactome-Circadian Repression of Expression By Rev ERBA”(55) were observed as significantly enriched in KLS genetic association analyses (see Table 4). Notably these 2 pathways shared identical gene members except for single gene. Gene based association tests within these enriched pathways are presented in Dataset S12. Both ROR Alpha (retinoic acid receptor-related alpha) and Rev ERBA are critical transcriptional regulators of periodic expression of BMAl1 gene, directly regulating circadian pathways (56-59). Other circadian pathways of interest included “Reactome-Bmal1 Clock NPAS2 Activates Circadian Expression” but did not reach 5% FDR threshold.

KLS Exome sequencing analysis

Targeted exome sequencing was performed in 32 individuals with KLS derived from 16 multiplex families (see Dataset S13A). The variant frequencies and variant by variant analyses compared to 286 European controls from 1000 genomes phase 3 were largely unremarkable. Damaging or probably damaging variants in each gene were collapsed and mutational burden in KLS probands and/or shared by at least 3 families were identified and presented below.

LMOD3 gene variants are not associated with KLS

Following up on a recent study suggesting increased mutation load in the LMOD3 gene in familial KLS cases (60), we analyzed the LMOD3 burden in our cohort of KLS probands and sporadic cases (n=17) compared to controls (n=286). We observed no significant differences in mutational burden in LMOD3 gene between cases and controls (p=0.13; see Dataset S13B). Variant level analyses in KLS probands and sporadic cases (n=17) revealed a previously reported single site rs35740823 (chr3:69171290; p.R83H) to be enriched in KLS cases (allelic OR=3.8; p=0.06) as compared to controls (60) (see Dataset S13B). We therefore imputed this LMOD3 variant rs35740823 and analyzed for associations in 151 KLS cases and 478 controls derived from replication cohort and predominantly of European origin but did not find any remarkable differences (OR=1.1; p=0.77; Dataset S13B).

Gene mutational burden test in KLS

We also compared the overall frequency of modifier variants tagged using SnpEff (61) as moderate or high effects in all genes in all familial KLS cases versus 1000 genome controls (n=286) running gene based mutational burden analyses (Dataset S14A), matched for ethnicity and SNP principal components. Although significant effects were found, none were found in genes within the TRANK1 loci. We however observed increased mutational burden in the rare alleles in genes NRXN2, GABRG2, CHORDC1, AATK and CASP7 genes (FDR p<1×10−4) among probands (n=17) relative to 1000 genomes reference samples (Dataset S14A). Genes with excessive rare alleles across each KLS family were next analyzed. Notably, GABRG2, CASP7 and ANKRD31 among other genes were shared across at least 3 KLS families (see Dataset S14B).

Further, pathway burden test from the list of 21 genes belonging to the circadian regulation pathway identified through the MAGMA analysis revealed significantly increased minor allele burden (see Dataset S15A) in the whole KLS cohort (n=32) compared to European origin controls. We also used pathway sets from msigDB6.1 to profile mutational burden in each pathway in KLS probands (n=14). Notably, we found increased mutational burden in GO GAMMA AMINOBUTYRIC ACID SIGNALING PATHWAY (FDR p=0.007) in KLS probands among other pathways that included many neuronal pathways (see Dataset S15B). This agrees with results of the gene burden test where we observed significantly increased burden in GABRG2 in KLS probands.

Discussion

Our study is the first to explore the interaction of genetic and environmental contributions to the pathophysiology of KLS. Results of GWAS revealed that TRANK1 gene variants were predisposing and associated with birth difficulties in KLS. PRS constructed on discovery GWAS analysis found strongly replicable PRS that strikingly overlapped with bipolar disorder and schizophrenia. Pathway analysis also revealed the involvement of behavioral and circadian rhythmicity genes.

This large study did not confirm findings made in other smaller genetic studies of KLS. For example, a 20-year-old study had suggested an HLA-DQ2 association (12) in KLS, but this was not confirmed either in other studies (3, 8, 13) nor in our current analyses. A recent study (60) focused on a large Saudi pedigree KLS affected family where a rs111848977 p.E142D mutation in LMOD3 was observed to cosegregate. This study also reported an increased burden of damaging mutations, including p.R23H in LMOD3 in KLS, a gene involved in nemaline myopathy 10 (62). In our sample of 32 familial cases, we could not find any overall increase in the frequency of LMOD3 mutations and, although p.R83H was nominally increased in familial cases, we did not find this association in a larger cohort of 151 KLS cases (see Dataset S13B).

Although the phenomenology of KLS is quite different from that of bipolar disorder or schizophrenia, these syndromes have some clinical similarities. Most notably, in both KLS and bipolar disorder, patients are asymptomatic between episodes. Similarly, exacerbations are usually episodic in some schizophrenia. In KLS, however, onset/offset of episodes is rapid (typically within one day), while it is generally slower in most bipolar patients, except for “rapid cyclers” (63). Further, although not a dominant feature, depression, anxiety, or hallucinations/delusions can occur during KLS episodes (3, 4, 6, 64) and at the end of episodes, some patients may have brief (1-2 days) period of elation, insomnia, and logorrhea which can resemble a maniac mini-episode.

Most interestingly, hypersomnia is a common feature of depressive episodes in bipolar type 1 disorder (65-67), unlike in other depressions where patients most often exhibit insomnia. Like KLS, bipolar disorder is also associated with cognitive deficiencies, although these are milder (68, 69) as compared to KLS patients. This, together with the fact that lithium is the only treatment reported to have an effect in preventing KLS episodes (7), suggests possible shared pathophysiology with bipolar disorder. Pathophysiological overlap among KLS, bipolar and schizophrenia is also in line with genetic studies in the United Kingdom (UK) Biobank cohort that have shown that the genetic architecture of long self-reported sleep duration overlaps with that of bipolar disorder, in contrast with genetic architecture of insomnia that shares similarity with anxiety and Restless Leg Syndrome (70). Most strikingly, strong correlations between PRS scores of long sleep duration and bipolar disorder type 1 in the UK Biobank sample was found (71). This suggests that hypersomnia, periodic or not, could share some commonality with bipolar disorder and schizophrenia.

One possible hypothesis that KLS, like bipolar disorder and other psychosis, is characterized by periodic instability of specific brain networks involving, in the case of KLS, reduced oxygen flow in cortical and subcortical areas, as reported by imaging studies(72). In this syndrome, cognitive and derealization symptoms may correlate with decreased activity in cortical areas, whereas decreased activity in thalamus and hypothalamus may be involved in the functional hypersomnia observed in disorder (72-74). Involvement of circadian or behavioral rhythmicity genes may be essential for maintaining hypersomnia episodes once an episode is initiated, for example, if subjects are suddenly spending much time in bed, unexposed to light, a behavior which could result in difficulties entraining in genetically susceptible individuals (75). This hypothesis could also explain why patients with KLS can exhibit Sleep Onset REM Periods during MSLTs (the circadian clock strongly drives occurrence of REM sleep) (76, 77), a feature also associated with shift work (78). Other investigators have found abnormal circadian entrainment during KLS episodes using actigraphy (79). A role for circadian abnormalities has also been suggested in bipolar disorder type 1 (67, 80, 81). Subjects with bipolar type 1 also reported increased eveningness tendencies (71), although this was not found in patients with KLS (3).

The most striking finding of this study is the observation of a strong genome-wide significant association of KLS with TRANK1-region variants [OR∼1.50] that are strongly associated with bipolar disorder [OR∼1.12] (43, 82, 83) and schizophrenia [OR∼1.08] (46, 84-86), but not mood instability (87). Although it remains possible this effect could be due to a rarer variant linked with this haplotype in KLS, this is more likely mediated by a similar effect of these common variants across these diseases. Of notable importance is the fact the TRANK1 effect size in this study is 4-5 orders of magnitude higher (OR∼1.5) than in studies in these other pathologies (OR∼1.1), explaining detection of association in a small sample of less than 1,000 KLS patients. This association with KLS is thus unlikely to be explained by contamination of our current KLS cohort with bipolar cases, although the opposite, contamination of bipolar or schizophrenia samples with KLS patients could explain weaker associations in these other pathologies.

The reason for the weaker than expected TRANK1 replication in the KLS replication cohort became evident after stratifying our discovery sample by birth decade, finding that the TRANK1 association was highest in subjects born before the 1980s and decreased thereafter. Further, KLS subjects with TRANK1 polymorphism have reported increased birth difficulties (summarized in Dataset S5). We hypothesize that improved perinatal outcomes that have occurred in the past 50 years may explain this tendency (34-39); https://www.cdc.gov/nchs/data/nvsr/nvsr56/nvsr56_10.pdf; https://mchb.hrsa.gov/chusa14/health-status-behaviors/infants/infant-mortality.html). The TRANK1 birth difficulties relationship was also confirmed in the independent replication sample, where we found that only the subset of KLS patient reporting birth difficulties had a significant TRANK1 association. Overall, patients who reported birth difficulties had the highest association (Fisher’s test, OR=1.57; p=2.5×10−4) both in the discovery and replication cohorts (Dataset S9A-B). TRANK1 association also increases disease severity during episodes independently of birth difficulty effects (Table 3).

Interestingly, although the association with TRANK1 was the highest genome significant finding in bipolar disorder in a recent metanalysis of 32 cohorts from 14 countries in Europe, North America and Australia, totaling 20,352 cases and 31,358 controls of European descent (rs9834970; OR=1.10, p=5.5 × 10−14) (15), but, as in our current KLS sample, replication in follow-up bipolar disorder samples totaling 9,412 cases and 137,760 controls of European ancestry did not replicate (OR=1.02, p=0.15). This finding deviated more than expected by chance (p=0.01, although not significant if Bonferroni corrected for all loci detected in the study). We hypothesize a similar time dependent TRANK1 birth difficulty interaction may play a role in bipolar disorder and schizophrenia. Indeed, these disorders are known to be strongly associated with early life complications (17-28, 88, 89), and data in the UK Biobank cohort indicates a dose-dependent increased prevalence of major depressive disorders and bipolar disorders in association with reported low birth weight (90), even when controlled for a history of the physical disorder known to be associated with low birth weight.

Two possibilities, which are not necessarily mutually exclusive, could explain this data. First, TRANK1 region genes could affect placental physiology, modulating the risk of birth difficulties and long-term fetal damage that could manifest as KLS or other disorders. Ursini et al. (31) found strong significant genetic interaction between the PRS of schizophrenia and the occurrence of early life complications in multiple case-control samples. Further, MLH1, LRRFIP2, and GOLGA4 associated with these overall genetic components were highly expressed in the placenta and modulated by pathological placental processes such as preeclampsia and intrauterine growth restriction but not in other disease processes (31). Problematically, however, recent papers did not observe that the KLS associated SNPs identified in this study had substantial effects on placental expression of these genes (31, 32, 91, 92). Further, TRANK1-region polymorphisms are not known maternal or fetal genetic factors associated with low birth weight (93-95), prematurity (96), or birth complications (91, 97) in large samples, and were not associated with birth difficulties in a small sample of 198 controls similarly ascertained to our KLS subjects for birth difficulties. Finally, this hypothesis would not explain the time dependency/cohort effect of the TRANK1 association we observed with KLS, which is likely due to improved or changes in perinatal care with time.

Another possible explanation, which we favor as it better fits our results and the lack of effects of TRANK1 on birth difficulties in previously published studies, is that newborns carrying TRANK1-region polymorphisms are more susceptible to injury in the presence of birth difficulties (notably hypoxia) in a way that produces behavioral instability such as KLS or other psychiatric conditions such as bipolar disorder in adulthood. This would be a true gene x environment interaction, although our study lacked sufficient control to formally test this hypothesis statistically. The mechanism underlying this phenomenon, which could have evolutionary advantages in times of stress (98), remains to be investigated. Interestingly, rs9834970, an SNP linked with the TRANK1 region susceptibility haplotype is associated with decreased expression of TRANK1 in iPSCs and neural progenitor cells (99); such effect has not been reported in the adult brain, but, as discussed, TRANK1 associated polymorphism reduces expression in various dividing tissues such as blood, esophageal mucosa, and B cells. The iPSC/neuronal precursor cell effect was rescued by valproic acid but not lithium, and the TRANK1 gene knockdown was found to modulate the expression of many genes, suggesting developmental effects (99).

Because the susceptibility haplotype modulates expression of many other genes in a tissue-specific manner, however, we believe the study of other genes in this region across tissues and development with the susceptibility haplotype will be necessary before concluding that TRANK1 is the gene behind this genetic association which is a limitation of our analyses. Additional studies aiming at studying more precisely which types of birth difficulties are associated with TRANK1 in these disorders is also needed. Our results offer new clues on the developmental window in which the TRANK1-susceptibility haplotype could affect risk for KLS and mental disorders in adulthood and mandates a more careful analysis of what and how the reported birth difficulties interact with specific TRANK1 genetic variants.

Our data support a genetic predisposition in the presence of TRANK1 polymorphisms and links KLS, behavioral rhythmicity, and bipolar disorder. Polygenic risk scores derived from discovery KLS samples are strongly predictive of the KLS phenotype in the replication sample. Further, PRSs support a genetic overlap with bipolar disorder and schizophrenia with the caveat being that these associations in various p thresholds were unadjusted for multiple comparisons and likely are underpowered since only replication follow up KLS sample was used. Strikingly, the temporal effect of TRANK1 polymorphism was observed in KLS cases in relation to the year of birth. In the presence of birth difficulties, the TRANK1 polymorphism may predispose to pathophysiological consequences such as KLS, bipolar disorder, and schizophrenia.

Materials and Methods

Inclusion criteria for KLS and origin of cases

Patients with KLS were included if they met International Classification of Sleep Disorders (ICSD3) criteria for the disorder including the presence of recurrent episodes (≥2) of severe hypersomnia (>18/24h) lasting days to weeks, associated with cognitive disturbances (confusion, amnesia, derealization), behavioral abnormalities (major apathy, megaphagia, and hypersexuality). These episodes had to be interspersed by long periods of normal sleep, cognition, behavior, and mood. A sleep-trained physician-diagnosed all patients. They were also asked to answer a detailed questionnaire, or a systematic interview as described (3, 6). Informed consent, in accordance with governing institutions, was obtained from all subjects. The research protocol was approved by Institutional Review Board (IRB) panels on Medical Human Subjects at Stanford University, and by respective IRB panels in each country providing samples for the study.

KLS cases

The discovery sample of 673 KLS cases includes 261 cases recruited in the United States (US), with the help of the Kleine Levin Syndrome Foundation (KLSF). US patients were diagnosed by an outside physician or referred by the KLSF (Steve Maier, Neal Farber), and the diagnosis confirmed at Stanford, where questionnaire and blood sample was collected. Other major recruitment sites include Pitié-Salpêtrière, Paris, France (Dr. Isabelle Arnulf, n=140), University of Peking, Beijing, People’s Republic of China, (Dr. Fang Han, n=62), Taiwan University, Taipei, Republic of China, Taiwan (Dr. Yu-Shu Huang, n=58), Guy de Chauliac Hospital, Montpellier, France (Dr. Yves Dauvilliers, n=35), Safra Children’s Hospital, Sheba Medical Center, Tel Aviv University, Tel Aviv, Israel (Dr. Yakov Sivan, n=34), Hepata Clinic, Germany (Dr. Gert Mayer, n=24), Uppsala University, Sweden (Dr. Anne-Marie Landtblom, n= 15), Hospital Robert Debre, Paris, France(Dr. Michel Lecendreux, n=7), Tokyo Metropolitan Institute of Medical Science Tokyo, Japan (Dr. Makoto Honda, n=7),University of Bologna, Italy (Drs Fabio Pizza and Guiseppe Plazzi, n=6), Maynei Hayeshua Medical Center, Bnei Barak, and The Sackler Faculty of Medicine, Tel-Aviv University (Dr. Natan Gadoth, n=6), Akita University, Akita, Japan (Dr. Kanbayashi, n=2), St. Vincent’s Hospital,The Catholic University of Korea (Dr. Seung Chul Hong, n=2), 1st Faculty of Medicine and General Teaching Hospital, Prague, Czech Republic (Dr. Sona Nevsimalova, n=2), Bambino Gesù Children Hospital, Rome, Italy (Dr. Marina Macchialio, n=2) and other centers who contributed single sample (n=10).

The replication cohort included 171 additional KLS cases. These included cases from Pitié-Salpêtrière, Paris, France (Dr. Isabelle Arnulf, Dr. Smaranda Leu-Semenescu and Dr. Pauline Dodet, n=66), Uppsala University, Sweden (Dr. Anne-Marie Landtblom, n= 26), Kleine Levin Syndrome Foundation (n=23), University of Bologna, Italy (Drs Fabio Pizza and Guiseppe Plazzi, n=13), London Bridge Hospital, London, UK (Dr. Guy Leschziner n=11), Guy de Chauliac Hospital, Montpellier, France (Dr. Yves Dauvilliers, n=9), Emory University, Atlanta (Dr. David Rye n=6), Hepata Clinic, Germany (Dr. Gert Mayer, n=5), Stanford University, (Dr. Mignot, n=5), Tokyo Metropolitan Institute of Medical Science, Tokyo, Japan (Dr. Makoto Honda n =4) and other centers (n=3).

GWAS patients and matched controls

Because the samples were collected over a long-time frame (2003-2019), genotyping platforms evolved, and hence each platform required its cohort of cases and controls (Dataset S1). As a result, the discovery sample (n=673, 2003-2017) was subdivided into seven subcohorts of cases and controls based on genotyping platform/ethnicity differences. Cases and ethnically matched controls (n=15,339) were genotyped at UCSF genomics core with Affymetrix Axiom® Genome-Wide LAT, EUR, CHB, Axiom® 55066 arrays, Axiom® world arrays, while multiethnic Axiom® PMRA arrays were genotyped at Akesogen inc. These include three European ancestry cohorts, an Asian cohort, a Latino cohort, and Ashkenazi and Non-Ashkenazi Jewish cohorts. The replication sample (n=171, collected between the years 2017-2019), subcohort 8, was treated as a single sample and genotyped solely on Axiom® PMRA arrays. DNA was extracted using Qiagen DNeasy kit, and samples were stored at -20c until genotyping, a minimum of ∼200ng was used for genotyping.

Genotype analysis and imputation

Raw CEL files were imported into Affymetrix Genotyping console, and Affymetrix best practices workflow was followed. Genotype calls were exported as plink ped/map files (100). For CEL files generated from Axiom® world arrays, primarily AffyPipe was used to call the genotypes and exported as plink ped/map files (101). Stringent QC was applied and variants and individuals with >0.05 missingness excluded, related individuals with identity by descent pihat>0.2 excluded and markers with Hardy-Weinberg equilibrium p<1e-6 in controls excluded. HapMap and 1000 genomes phase 3 reference population genotypes(102) were used to guide ethnicity matching in our sample sub cohorts. Population covariates were extracted using plink and visually inspected for unexpected deviations and outliers, after which Euclidean distance was used to accurately and automatically match cases and controls. Post QC, genotypes were phased one chromosome at a time with shapeit (103, 104) using the build37 coordinate system of the human genome. Phased haplotypes were then subject to imputation to 1000 Genomes phase 1 in 1mb chunks for each chromosome using impute2 (105).

Association testing and variant prioritization

Frequentist association tests were then performed on the imputed genotypes by including genetic principal components (PCA) and missing data proportion using snptest (106). Meta-analysis was subsequently performed using cohort specific β estimates using a fixed effects model(106). The genome-wide significance threshold was set at p<5e-8 based on Bonferroni adjustment. Gene-based association tests, including pathway analysis were performed with MAGMA modules (107). Manhattan plots were constructed with R package qqman (108). Additional eQTL annotation was carried out using SNIPA, webserver (109), and GTEx v8 (30). eCAVIAR was used to perform eQTL colocalization (110). LD hub was used to calculate LD score regression and find shared SNP heritability across phenotypes (111). Polygenic risk scores were computed using PRSice2 (48), briefly imputed genotypes in BGEN format from the replication data set were quality controlled (maf >0.05 and impute info score >=0.9) and clumped for linkage disequilibrium (250 kb window; r2 = 0.1) and used as target dataset along with first five genetic PCAs. The common variants remaining after clumping from the base and target dataset were 104626. Apart from using summary statistics from the KLS discovery, we also used summary statistics from various other phenotypes. PRSs at various p-value thresholds (5×10−10,5×10−8,1×10−7,1×10−5, 0.0001, 0.001, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5 and 1) were derived and validated in KLS replication cohort (n=171).

Exome sequencing sample and matched controls

Agilent Sure Select V4 system workflow was used to capture regions of interest; the resulting fragments were subject to Illumina paired-end sequencing run with an average depth of 130x and average read length of 100 in the target regions for the KLS family DNA samples. The KLS family DNA consisted of 32 individuals derived from 16 multiplex families (see Dataset S13A). On the other hand, for controls, exome FASTQ files from 286 individuals part of 1000 genomes (predominantly CEU, TSI, GBR, IBS & FIN) were used. The FASTQ files were aligned to UCSC hg19 reference FASTA using BWA (112) and read groups inserted for downstream variant quality score recalibration. GATK best practices workflow (https://software.broadinstitute.org/gatk/best-practices/) was implemented, and joint genotyping performed on both cases and controls simultaneously. Raw variant call sets were subject to variant quality score recalibration, and variants passing the filters were then imported into plinkseq (https://atgu.mgh.harvard.edu/plinkseq/) and used in downstream association analysis. Annovar and SnpEff were used to annotate variants (61, 113). Following which moderate and high modifier variants were collapsed and gene-based burden tests in various sub-cohorts were carried out using plinkseq.

A gene burden test was carried out within each family, probands only and all the 32 individuals. Each analysis used the 1000 genomes European dataset matched by principle components as a control cohort. The burden analyses were focused only in rare variants (minor allele counts 2-20) and restricted to predicted to be moderate or high modifiers as tagged by Annovar and SnpEff. Further, genes part of the rhythmic behavior pathway were collapsed, and SKAT based gene burden tests were performed (114) with resulting p-values adjusted using a 5% false discovery rate (FDR).

Taqman genotyping

To validate the genotype calls from SNP genotyping arrays, we conducted Taqman assays for TRANK1 SNPs rs71947865 and rs6550435 on a random subset of KLS cases and controls typed across multiple platforms (n=122). Allele concordance for rs71947865 was 100%, while for rs6550435, it was 97% in 122 subjects.

Clinical Associations

Demographic, and various clinical features within the KLS cohort (n=592) were analyzed for associations with TRANK1 rs71947865 variant (deletion). Three kinds of analyses were performed, model1; A linear additive model was fit as a function of the clinical variable, Model2; same as model1 except for including covariates such as ethnicity/referrer code (geographical location) and missing clinical data proportion per individual, Fisher’s exact test of allelic counts was performed for all categorical clinical variables.

Birth Year Associations

Two analysis strategies were pursued, first, when birth year was available for KLS cases (n=650), we stratified the KLS cases into three approximately equal-sized groups; 1987 or before (n=220, median year of birth: 1978), 1988-1994 (n=202, median year of birth: 1992) and 1995-2000’s (n=228, median year of birth: 1999). Within each group, we matched control samples based on genetic principle components. Fisher’s exact and chi-square tests on TRANK1 rs71947865 allelic counts were performed to find significant difference in allele frequencies. A Fisher’s test of TRANK1 SNPs allelic counts was performed within KLS cases stratified by the 3 birth year bins, the minor allele counts in birth year bin 1987 or earlier was compared to minor allele counts 1988-1994 and 1995-2000’s in 2×2 contingency table.

Data Availability

All genotyping calls can be requested by contacting the corresponding author.

Data availability

All genotyping calls can be requested by contacting the corresponding author.

Supplementary Information

Fig. S1.
Fig. S1.

A: The first 2 genetic principle components matching of case-controls in KLS discovery and (B) replication sub -cohorts stratified by presence of birth difficulties. (C) QQplot of the total KLS case-controls combined discovery + replication genetic association p-values.

Fig. S2.
Fig. S2.

(A-C): TRANK1 regional plots in the European origin KLS meta-analysis and sub-cohorts with and with our birth difficulties. (D) – TRANK1 loci plots with KLS discovery in the positive fraction contrasted with Bipolar Disorder-Schizophrenia (BD-SCZ) GWAS summary stats in the negative fraction(1) (Ruderfer DM et al. 2018), horizontal lines indicate GWAS significance threshold (p=5e-8).

Fig. S3.
Fig. S3.

PRS (polygenic risk scores) constructed from various other GWAS studies i.e. bipolar disorder, schizophrenia, circadian chronotype, sleep duration and excessive daytime sleepiness, evaluating their predictive performance on the KLS follow up cohort. EurKLS – PRS built on Caucasian discovery and validated in Caucasian replication KLS samples. The y-axis is the % variance explained by the PRS while the x-axis represents the p-value thresholds. Inset are marked the Pvalue of the fit. The * on the top of the bars represent the empirical p value (p<0.05) associated with the test of p-value threshold.

Dataset S1: The total cohort used in the analysis split by sub-cohort and genotyping array used grouped by KLS cases and controls.

Dataset S2: Variants in the TRANK1 loci associated with KLS in worldwide discovery cohort GWAS meta-analysis in 673 KLS cases and 15341 control individuals and in a Follow-up replication cohort of 171 KLS cases and 1956 controls. The effect allele was B*.

Dataset S3: TRANK1 region SNPs correlation in schizophrenia and bi-polar disorder GWAS studies. * indicates Genome wide significant (p<5e-8). PGC is Psychiatric Genomics Consortium.

Dataset S4: TRANK1 deletion (rs71947865) and its univariate association with KLS clinical variables. We used 2 linear additive models Model1-the additively coded minor allele genotype was fit as a function of the clinical variable, Model2 – was similar to model1 except for including ethnicity, Gender and missing data proportion. A Fisher’s exact test of allelic counts was further performed for categorical variables. In these analyses 592 KLS cases were included. (Frequency % | mean ± sem was computed). BB-allele G, AB-GGGAGCCA/G, AA - GGGAGCCA/ GGGAGCCA, *Minor Allele G (deletion) was the effect allele.

Dataset S5: Birth Difficulty and developmental delay summary in the KLS cases.

Dataset S6: (A) TRANK1 SNPs in birth year stratified analysis, Fisher’s exact count test was performed for case -control analyses. Minor *A allele was effect allele. (B) Allelic counts of TRANK1 SNPs in KLS cases compared to each birth year bin, Fisher’s exact count test was performed for case only analyses. Minor A allele was effect allele.

Dataset S7: Variants in the TRANK1 loci associated with KLS in worldwide discovery cohort GWAS meta-analysis in 673 KLS cases and 15339 control individuals and Follow-up replication sub-cohorts based on presence of birth difficulties. The effect allele was B*. genomic coordinates in hg19/GRCh37.

Dataset S8: (A) TRANK1 SNPs association in KLS cases with birth difficulties compared to controls stratified by discovery and the replication cohort, Fisher’s exact count test was performed. Minor *A allele was effect allele. (B) TRANK1 SNPs in birth difficulty stratified analysis in the combined discovery and replication cohort, Fisher’s exact count test was performed. Minor *A allele was effect allele.

Dataset S9: Meta-analysis of Variants in the TRANK1 loci in European origin KLS cases in 2 cohorts stratified by all, presence of birth difficulties and absence of birth difficulties. Effect allele is B*. Filter on Analysis column to see associations by birth difficulties. genomic coordinates in hg19/GRCh37.

Dataset S10: Polygenic risk score analysis of KLS follow up sample in relation to other disorders including Bipolar Disorder, Schizophrenia, Circadian chronotype and Excessive daytime sleepiness.

Dataset S11: LD hub output containing LD score regression and shared SNP heritability with KLS across phenotypes, the table lists the top traits that share heritability with KLS.

Dataset S12: Top 2 enriched pathways in a MAGMA based gene set enrichment analysis was carried out using A total of 5346 sets including curated gene sets, KEGG pathways, Reactome pathways and GO biological processes derived from msigdb6.1 were used. Genomic coordinates in hg19/GRCh37.

Dataset S13: (A) The KLS cases and families (n=32) used in the exome sequencing studies. Controls were derived from 1000 genomes phase 3. (B): Frequency and burden of all LMOD3 gene variants derived from whole exome sequencing in 17 KLS Probands/sporadic cases and 286 European controls. Genomic coordinates in hg19/GRCh37.

Dataset S14: (A) Exome burden test restricted to rare variants with high and moderate effects grouped by reference genes within each families, probands and total cohort are listed, Columns with prefix KLS implies all the 32 cases and column prefix proband includes the proband in each family in addition to the 3 sporadic cases (n=17). genomic coordinates in hg19/GRCh37.(B): Exome burden test restricted to rare variants with high and moderate effects grouped by reference genes shared by at least 3 families, Columns with prefix KLS implies all the 32 cases and column prefix proband includes the proband in each family in addition to the 3 sporadic cases (n=17). genomic coordinates in hg19/GRCh37.

Dataset S15: (A) Top hits from Exome SKAT model with population covariates in 32 KLS cases and 286 controls with a focus on collapsing variants in the genes belonging to REACTOME_RORA_ACTIVATES_CIRCADIAN_EXPRESSION and REACTOME_CIRCADIAN_REPRESSION_OF_EXPRESSION_BY_REV_ERBA both of which were 5% FDR significant in the joint Discovery-Replication KLS GWAS pathway analyses. (B): Top hits from Exome SKAT model with population covariates in 17 KLS probands and 286 controls with a focus on collapsing variants in the genes in pathways derived from a total of 5346 sets including curated gene sets, KEGG pathways, Reactome pathways and GO biological processes derived from msigdb6.1 were used.

Acknowledgments

We thank countless staff members, clinicians, patients, and patient advocates who are not listed here and have provided material and logistic support over decades for these 15 yearlong studies. The study was primarily funded by the Kleine-Levin Syndrome Foundation (KLSF, 2003-2020) and by NIH 1R01MH080957 (2007-2013) to EM. The French Kleine-Levin Syndrome research program is financed by a grant to IA (PHRC 070138) from the French Ministry of Health. HMO has been supported by the Academy of Finland (#309643), Yrjö Jahnsson Foundation and Oskar Öflund Foundation. Computing for this project was performed on the Stanford Genomics cluster supported by NIH 1S10OD023452-01. We would like to thank the Stanford Research Computing Center for providing computational resources and support that contributed to these research results.

References

  1. 1.↵
    D. W. Kleine, Periodische Schlafsucht. pp. 305–320. European Neurology 57, 305–320 (1925).
    OpenUrl
  2. 2.↵
    M. Levin, Periodic somnolence and morbid hunger: A new syndrome. Brain 59, 494–504 (1936).
    OpenUrlCrossRefWeb of Science
  3. 3.↵
    I. Arnulf et al., Kleine-Levin syndrome: a systematic study of 108 patients. Ann Neurol 63, 482–493 (2008).
    OpenUrlCrossRefPubMed
  4. 4.↵
    I. Arnulf, T. J. Rico, E. Mignot, Diagnosis, disease course, and management of patients with Kleine-Levin syndrome. Lancet Neurol 11, 918–928 (2012).
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    M. Billiard, I. Jaussent, Y. Dauvilliers, A. Besset, Recurrent hypersomnia: a review of 339 cases. Sleep Med Rev 15, 247–257 (2011).
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    I. Arnulf, J. M. Zeitzer, J. File, N. Farber, E. Mignot, Kleine-Levin syndrome: a systematic review of 186 cases in the literature. Brain 128, 2763–2776 (2005).
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    S. Leu-Semenescu et al., Lithium therapy in Kleine-Levin syndrome: An open-label, controlled study in 130 patients. Neurology 85, 1655–1662 (2015).
    OpenUrl
  8. 8.↵
    Y. S. Huang et al., Relationship between Kleine-Levin syndrome and upper respiratory infection in Taiwan. Sleep 35, 123–129 (2012).
    OpenUrlPubMedWeb of Science
  9. 9.↵
    Q. T. Nguyen et al., Familial Kleine-Levin Syndrome: A Specific Entity? Sleep 39, 1535–1542 (2016).
    OpenUrl
  10. 10.↵
    R. Peraita-Adrados, J. L. Vicario, M. Tafti, M. Garcia de Leon, M. Billiard, Monozygotic twins affected with Kleine-Levin syndrome. Sleep 35, 595–596 (2012).
    OpenUrl
  11. 11.↵
    O. Habra, R. Heinzer, J. Haba-Rubio, A. O. Rossetti, Prevalence and Mimics of Kleine-Levin Syndrome: A Survey in French-Speaking Switzerland. J Clin Sleep Med 12, 1083–1087 (2016).
    OpenUrl
  12. 12.↵
    Y. Dauvilliers et al., Kleine-Levin syndrome: an autoimmune hypothesis based on clinical and genetic analyses. Neurology 59, 1739–1745 (2002).
    OpenUrlCrossRefPubMed
  13. 13.↵
    S. Lavault et al., Kleine-Levin syndrome in 120 patients: differential diagnosis and long episodes. Ann Neurol 77, 529–540 (2015).
    OpenUrlPubMed
  14. 14.↵
    M. Lam et al., Comparative genetic architectures of schizophrenia in East Asian and European populations. Nat Genet 51, 1670–1678 (2019).
    OpenUrl
  15. 15.↵
    E. A. Stahl et al., Genome-wide association study identifies 30 loci associated with bipolar disorder. Nat Genet 51, 793–803 (2019).
    OpenUrlCrossRefPubMed
  16. 16.↵
    M. Ikeda et al., Genome-Wide Association Study Detected Novel Susceptibility Genes for Schizophrenia and Shared Trans-Populations/Diseases Genetic Effect. Schizophr Bull 45, 824–834 (2019).
    OpenUrl
  17. 17.↵
    S. A. Stilo, R. M. Murray, Non-Genetic Factors in Schizophrenia. Curr Psychiatry Rep 21, 100 (2019).
    OpenUrl
  18. 18.
    J. Soderlund, S. Wicks, L. Jorgensen, C. Dalman, Comparing cohort incidence of schizophrenia with that of bipolar disorder and affective psychosis in individuals born in Stockholm County 1955-1967. Psychol Med 45, 3433–3439 (2015).
    OpenUrl
  19. 19.
    R. Chudal et al., Perinatal factors and the risk of bipolar disorder in Finland. J Affect Disord 155, 75–80 (2014).
    OpenUrlCrossRefPubMed
  20. 20.
    U. K. Haukvik et al., Pre- and perinatal hypoxia associated with hippocampus/amygdala volume in bipolar disorder. Psychol Med 44, 975–985 (2014).
    OpenUrl
  21. 21.
    C. Nosarti et al., Preterm birth and psychiatric disorders in young adult life. Arch Gen Psychiatry 69, E1–8 (2012).
    OpenUrlCrossRefPubMed
  22. 22.
    J. McGrath, J. Scott, Urban birth and risk of schizophrenia: a worrying example of epidemiology where the data are stronger than the hypotheses. Epidemiol Psichiatr Soc 15, 243–246 (2006).
    OpenUrlPubMed
  23. 23.
    M. Cannon, P. B. Jones, R. M. Murray, Obstetric complications and schizophrenia: historical and meta-analytic review. Am J Psychiatry 159, 1080–1092 (2002).
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.
    H. Kunugi, S. Nanko, R. M. Murray, Obstetric complications and schizophrenia: prenatal underdevelopment and subsequent neurodevelopmental impairment. Br J Psychiatry Suppl 40, s25–29 (2001).
    OpenUrl
  25. 25.
    D. K. Kinney, D. A. Yurgelun-Todd, M. Tohen, S. Tramer, Pre- and perinatal complications and risk for bipolar disorder: a retrospective study. J Affect Disord 50, 117–124 (1998).
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.
    J. R. Geddes, S. M. Lawrie, Obstetric complications and schizophrenia: a meta-analysis.Br J Psychiatry 167, 786–793 (1995).
    OpenUrlAbstract/FREE Full Text
  27. 27.
    T. F. McNeil, Perinatal risk factors and schizophrenia: selective review and methodological concerns. Epidemiol Rev 17, 107–112 (1995).
    OpenUrlPubMedWeb of Science
  28. 28.↵
    S. N. Brixey, B. J. Gallagher, 3rd., J. A. McFalls Jr..,, L. F. Parmelee, Gestational and neonatal factors in the etiology of schizophrenia. J Clin Psychol 49, 447–456 (1993).
    OpenUrlCrossRefPubMed
  29. 29.↵
    H. Shimizu, A. Hirose, T. Tatsuno, M. Nakamura, J. Katsube, Pharmacological properties of SM-3997: a new anxioselective anxiolytic candidate. Jpn J Pharmacol 45, 493–500 (1987).
    OpenUrlCrossRefPubMed
  30. 30.↵
    G. T. Consortium, The Genotype-Tissue Expression (GTEx) project. Nat Genet 45, 580–585 (2013).
    OpenUrlCrossRefPubMed
  31. 31.↵
    G. Ursini et al., Convergence of placenta biology and genetic risk for schizophrenia. Nat Med 24, 792–801 (2018).
    OpenUrl
  32. 32.↵
    S. Peng et al., Expression quantitative trait loci (eQTLs) in human placentas suggest developmental origins of complex diseases. Hum Mol Genet 26, 3432–3441 (2017).
    OpenUrlCrossRef
  33. 33.↵
    T. J. Hoffmann et al., Next generation genome-wide association tool: design and coverage of a high-throughput European-optimized SNP array. Genomics 98, 79–89 (2011).
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    A. Sander, R. Wauer, From single-case analysis of neonatal deaths toward a further reduction of the neonatal mortality rate. J Perinat Med 47, 125–133 (2018).
    OpenUrl
  35. 35.
    M. K. Sugai, S. Gilmour, E. Ota, K. Shibuya, Trends in perinatal mortality and its risk factors in Japan: Analysis of vital registration data, 1979-2010. Sci Rep 7, 46681 (2017).
    OpenUrl
  36. 36.
    E. B. Schneider, Fetal health stagnation: Have health conditions in utero improved in the United States and Western and Northern Europe over the past 150 years? Soc Sci Med 179, 18–26 (2017).
    OpenUrl
  37. 37.
    E. Verstraete et al., Healthcare-associated bloodstream infections in a neonatal intensive care unit over a 20-year period (1992-2011): trends in incidence, pathogens, and mortality. Infect Control Hosp Epidemiol 35, 511–518 (2014).
    OpenUrlCrossRefPubMed
  38. 38.
    S. V. Glinianaia et al., Temporal changes in key maternal and fetal factors affecting birth outcomes: a 32-year population-based study in an industrial city. BMC Pregnancy Childbirth 8, 39 (2008).
    OpenUrlCrossRefPubMed
  39. 39.↵
    D. Rush, J. M. Alvir, D. A. Kenny, S. S. Johnson, D. G. Horvitz, The National WIC Evaluation: evaluation of the Special Supplemental Food Program for Women, Infants, and Children. III. Historical study of pregnancy outcomes. Am J Clin Nutr 48, 412–428 (1988).
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    R. Secolin et al., Family-based association study for bipolar affective disorder. Psychiatr Genet 20, 126–129 (2010).
    OpenUrlPubMed
  41. 41.↵
    A. W. Charney et al., Evidence for genetic heterogeneity between clinical subtypes of bipolar disorder. Transl Psychiatry 7, e993 (2017).
    OpenUrl
  42. 42.↵
    L. Hou et al., Genome-wide association study of 40,000 individuals identifies two novel loci associated with bipolar disorder. Hum Mol Genet 25, 3383–3394 (2016).
    OpenUrlCrossRefPubMed
  43. 43.↵
    T. W. Muhleisen et al., Genome-wide association study reveals two new risk loci for bipolar disorder. Nat Commun 5, 3339 (2014).
    OpenUrlCrossRefPubMed
  44. 44.↵
    D. M. Ruderfer et al., Polygenic dissection of diagnosis and clinical dimensions of bipolar disorder and schizophrenia. Mol Psychiatry 19, 1017–1024 (2014).
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Z. Li et al., Genome-wide association analysis identifies 30 new susceptibility loci for schizophrenia. Nat Genet 49, 1576–1583 (2017).
    OpenUrlCrossRefPubMed
  46. 46.↵
    C. Schizophrenia Working Group of the Psychiatric Genomics, Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421–427 (2014).
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    D. T. Chen et al., Genome-wide association study meta-analysis of European and Asian-ancestry samples identifies three novel loci associated with bipolar disorder. Mol Psychiatry 18, 195–205 (2013).
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    S. W. Choi, P.F. O’Reilly, PRSice-2: Polygenic Risk Score software for biobank-scale data. Gigascience 8, giz082 (2019).
    OpenUrl
  49. 49.↵
    D. M. Ruderfer et al., Genomic Dissection of Bipolar Disorder and Schizophrenia, Including 28 Subphenotypes. Cell 173, 1705-+ (2018).
    OpenUrlCrossRefPubMed
  50. 50.↵
    S. E. Jones et al., Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nature Communications 10, 343 (2019).
    OpenUrl
  51. 51.↵
    H. S. Dashti et al., Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates. Nat Commun 10, 1100 (2019).
    OpenUrlCrossRefPubMed
  52. 52.↵
    H. Wang et al., Genome-wide association analysis of self-reported daytime sleepiness identifies 42 loci that suggest biological subtypes. Nat Commun 10, 3503 (2019).
    OpenUrl
  53. 53.↵
    C. A. de Leeuw, J. M. Mooij, T. Heskes, D. Posthuma, MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput Biol 11, e1004219 (2015).
    OpenUrlCrossRefPubMed
  54. 54.↵
    V. Giguere et al., Isoform-specific amino-terminal domains dictate DNA-binding properties of ROR alpha, a novel family of orphan hormone nuclear receptors. Genes Dev 8, 538–553 (1994).
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    P. Griffin et al., Circadian clock protein Rev-erbα regulates neuroinflammation. Proceedings of the National Academy of Sciences 116, 5102–5107 (2019).
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    Y. Zhang et al., Discrete functions of nuclear receptor Rev-erbα couple metabolism to the clock. Science 348, 1488–1492 (2015).
    OpenUrlAbstract/FREE Full Text
  57. 57.
    L. A. Solt, D. J. Kojetin, T. P. Burris, The REV-ERBs and RORs: molecular links between circadian rhythms and lipid homeostasis. Future Med Chem 3, 623–638 (2011).
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.
    F. Guillaumond, H. Dardente, V. Giguere, N. Cermakian, Differential control of Bmal1 circadian transcription by REV-ERB and ROR nuclear receptors. J Biol Rhythms 20, 391–403 (2005).
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    T. K. Sato et al., A functional genomics strategy reveals Rora as a component of the mammalian circadian clock. Neuron 43, 527–537 (2004).
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    S. M. Al Shareef et al., Kleine-Levin syndrome is associated with LMOD3 variants. J Sleep Res 28, e12718 (2019).
    OpenUrl
  61. 61.↵
    P. Cingolani et al., A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 6, 80–92 (2012).
    OpenUrlCrossRef
  62. 62.↵
    U. A. Schatz et al., Evidence of mild founder LMOD3 mutations causing nemaline myopathy 10 in Germany and Austria. Neurology 91, e1690–e1694 (2018).
    OpenUrl
  63. 63.↵
    R. Tillman, B. Geller, Definitions of rapid, ultrarapid, and ultradian cycling and of episode duration in pediatric and adult bipolar disorders: a proposal to distinguish episodes from cycles. J Child Adolesc Psychopharmacol 13, 267–271 (2003).
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    E. Groos et al., Emerging psychiatric disorders in Kleine-Levin syndrome. J Sleep Res 27, e12690 (2018).
    OpenUrl
  65. 65.↵
    K. A. Kaplan et al., Hypersomnia subtypes, sleep and relapse in bipolar disorder. Psychol Med 45, 1751–1763 (2015).
    OpenUrl
  66. 66.
    R. B. Grigolon et al., Hypersomnia and Bipolar Disorder: A systematic review and meta-analysis of proportion. J Affect Disord 246, 659–666 (2019).
    OpenUrl
  67. 67.↵
    L. Pagani et al., Genetic contributions to circadian activity rhythm and sleep pattern phenotypes in pedigrees segregating for severe bipolar disorder. Proc Natl Acad Sci U S A 113, E754–761 (2016).
    OpenUrlAbstract/FREE Full Text
  68. 68.↵
    E. Bora, Neurocognitive features in clinical subgroups of bipolar disorder: A meta-analysis. J Affect Disord 229, 125–134 (2018).
    OpenUrl
  69. 69.↵
    I. M. M. Lima, A. D. Peckham, S. L. Johnson, Cognitive deficits in bipolar disorders: Implications for emotion. Clin Psychol Rev 59, 126–136 (2018).
    OpenUrl
  70. 70.↵
    J. M. Lane et al., Genome-wide association analyses of sleep disturbance traits identify new loci and highlight shared genetics with neuropsychiatric and metabolic traits. Nat Genet 49, 274–281 (2017).
    OpenUrlCrossRefPubMed
  71. 71.↵
    K. J. S. Lewis et al., Comparison of Genetic Liability for Sleep Traits Among Individuals With Bipolar Disorder I or II and Control Participants. JAMA Psychiatry 10.1001/jamapsychiatry.2019.4079 (2019).
  72. 72.↵
    M. Engstrom, F. Latini, A. M. Landtblom, Neuroimaging in the Kleine-Levin Syndrome.Curr Neurol Neurosci Rep 18, 58 (2018).
    OpenUrl
  73. 73.
    S. B. Hong, Neuroimaging of Narcolepsy and Kleine-Levin Syndrome. Sleep Med Clin 12, 359–368 (2017).
    OpenUrl
  74. 74.↵
    A. Kas, S. Lavault, M. O. Habert, I. Arnulf, Feeling unreal: a functional imaging study in patients with Kleine-Levin syndrome. Brain 137, 2077–2087 (2014).
    OpenUrlCrossRefPubMed
  75. 75.↵
    C. Blume, C. Garbazza, M. Spitschan, Effects of light on human circadian rhythms, sleep and mood. Somnologie (Berl) 23, 147–156 (2019).
    OpenUrl
  76. 76.↵
    C. A. Czeisler, J. C. Zimmerman, J. M. Ronda, M. C. Moore-Ede, E. D. Weitzman, Timing of REM sleep is coupled to the circadian rhythm of body temperature in man. Sleep 2, 329–346 (1980).
    OpenUrlCrossRefPubMedWeb of Science
  77. 77.↵
    M. A. Carskadon, W. C. Dement, Sleep studies on a 90-minute day. Electroencephalogr Clin Neurophysiol 39, 145–155 (1975).
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.↵
    A. Goldbart et al., Narcolepsy and predictors of positive MSLTs in the Wisconsin Sleep Cohort. Sleep 37, 1043–1051 (2014).
    OpenUrl
  79. 79.↵
    L. Chen, Y.-S. Huang., Actigraphy study and different Hypersomnia disorders. 2nd Asian Narcolepsy & Hypersomnolence Society Meeting. 2, 39 (2019).
    OpenUrl
  80. 80.↵
    A. K. Gold, G. Kinrys, Treating Circadian Rhythm Disruption in Bipolar Disorder. Curr Psychiatry Rep 21, 14 (2019).
    OpenUrl
  81. 81.↵
    A. Vreeker et al., Genetic analysis of activity, brain and behavioral associations in extended families with heavy genetic loading for bipolar disorder. Psychol Med doi:10.1017/S0033291719003416, 1–9 (2019).
    OpenUrlCrossRef
  82. 82.↵
    M. Ikeda et al., A genome-wide association study identifies two novel susceptibility loci and trans population polygenicity associated with bipolar disorder. Mol Psychiatry 23, 639–647 (2018).
    OpenUrl
  83. 83.↵
    G. C. B. D. W. G. Psychiatric, Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4. Nat Genet 43, 977–983 (2011).
    OpenUrlCrossRefPubMed
  84. 84.↵
    A. J. Forstner et al., Identification of shared risk loci and pathways for bipolar disorder and schizophrenia. PLoS One 12, e0171595 (2017).
    OpenUrl
  85. 85.
    S. Ripke et al., Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet 45, 1150–1159 (2013).
    OpenUrlCrossRefPubMed
  86. 86.↵
    C. Schizophrenia Psychiatric Genome-Wide Association Study, Genome-wide association study identifies five new schizophrenia loci. Nat Genet 43, 969–976 (2011).
    OpenUrlCrossRefPubMed
  87. 87.↵
    J. Ward et al., Genome-wide analysis in UK Biobank identifies four loci associated with mood instability and genetic correlation with major depressive disorder, anxiety disorder and schizophrenia. Transl Psychiatry 7, 1264 (2017).
    OpenUrl
  88. 88.↵
    G. Weber, Investigation of metal-species (Ca, Mg, Zn, Fe, Cu, Pb, Cd, Sn) in urine by HPLC-AAS. J Trace Elem Electrolytes Health Dis 2, 61–65 (1988).
    OpenUrlPubMed
  89. 89.↵
    K. J. Tsuchiya, M. Byrne, P. B. Mortensen, Risk factors in relation to an emergence of bipolar disorder: a systematic review. Bipolar Disord 5, 231–242 (2003).
    OpenUrlCrossRef
  90. 90.↵
    D. M. Lyall et al., Low birth weight and features of neuroticism and mood disorder in 83 545 participants of the UK Biobank cohort. BJPsych Open 2, 38–44 (2016).
    OpenUrl
  91. 91.↵
    T. Kikas, K. Rull, R. N. Beaumont, R. M. Freathy, M. Laan, The Effect of Genetic Variation on the Placental Transcriptome in Humans. Front Genet 10, 550 (2019).
    OpenUrl
  92. 92.↵
    A. Brockenshire, The “mini-mental state”: a handy tool. Perspectives 11, 7–8 (1987).
    OpenUrl
  93. 93.↵
    M. Horikoshi et al., Genome-wide associations for birth weight and correlations with adult disease. Nature 538, 248–252 (2016).
    OpenUrlCrossRefPubMed
  94. 94.
    N. M. Warrington et al., Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors. Nat Genet 51, 804–814 (2019).
    OpenUrlCrossRefPubMed
  95. 95.↵
    R. N. Beaumont et al., Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics. Hum Mol Genet 27, 742–756 (2018).
    OpenUrlCrossRef
  96. 96.↵
    G. Zhang et al., Genetic Associations with Gestational Duration and Spontaneous Preterm Birth. N Engl J Med 377, 1156–1167 (2017).
    OpenUrlCrossRefPubMed
  97. 97.↵
    L. Srinivasan et al., Genome-wide association study of sepsis in extremely premature infants. Arch Dis Child Fetal Neonatal Ed 102, F439–F445 (2017).
    OpenUrlAbstract/FREE Full Text
  98. 98.↵
    J. R. Pasqualini, P. Mercat, N. Giambiagi, Histone acetylation decreased by estradiol in the MCF-7 human mammary cancer cell line. Breast Cancer Res Treat 14, 101–105 (1989).
    OpenUrlCrossRefPubMed
  99. 99.↵
    X. Jiang et al., Sodium valproate rescues expression of TRANK1 in iPSC-derived neural cells that carry a genetic variant associated with serious mental illness. Mol Psychiatry 24, 613–624 (2019).
    OpenUrl
  100. 100.↵
    S. Purcell et al., PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics 81, 559–575 (2007).
    OpenUrlCrossRefPubMed
  101. 101.↵
    E. L. Nicolazzi, D. Iamartino, J. L. Williams, AffyPipe: an open-source pipeline for Affymetrix Axiom genotyping workflow. Bioinformatics 30, 3118–3119 (2014).
    OpenUrlCrossRefPubMed
  102. 102.↵
    C. Genomes Project et al., A global reference for human genetic variation. Nature 526, 68–74 (2015).
    OpenUrlCrossRefPubMed
  103. 103.↵
    O. Delaneau, J. Marchini, J. F. Zagury, A linear complexity phasing method for thousands of genomes. Nat Methods 9, 179–181 (2012).
    OpenUrlCrossRefPubMedWeb of Science
  104. 104.↵
    O. Delaneau, J. Marchini, C. Genomes Project, C. Genomes Project, Integrating sequence and array data to create an improved 1000 Genomes Project haplotype reference panel. Nat Commun 5, 3934 (2014).
    OpenUrlCrossRefPubMed
  105. 105.↵
    B. N. Howie, P. Donnelly, J. Marchini, A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet 5, e1000529 (2009).
    OpenUrlCrossRefPubMed
  106. 106.↵
    J. Marchini, B. Howie, S. Myers, G. McVean, P. Donnelly, A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet 39, 906–913 (2007).
    OpenUrlCrossRefPubMedWeb of Science
  107. 107.↵
    K. Watanabe, E. Taskesen, A. van Bochoven, D. Posthuma, Functional mapping and annotation of genetic associations with FUMA. Nat Commun 8, 1826 (2017).
    OpenUrlCrossRefPubMed
  108. 108.↵
    S. D. Turner, qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. bioRxiv doi:10.1101/005165 (2014).
    OpenUrlAbstract/FREE Full Text
  109. 109.↵
    M. Arnold, J. Raffler, A. Pfeufer, K. Suhre, G. Kastenmuller, SNiPA: an interactive, genetic variant-centered annotation browser. Bioinformatics 31, 1334–1336 (2015).
    OpenUrlCrossRefPubMed
  110. 110.↵
    F. Hormozdiari et al., Colocalization of GWAS and eQTL Signals Detects Target Genes. American journal of human genetics 99, 1245–1260 (2016).
    OpenUrlCrossRefPubMed
  111. 111.↵
    J. Zheng et al., LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics 33, 272–279 (2017).
    OpenUrlCrossRefPubMed
  112. 112.↵
    H. Li, R. Durbin, Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).
    OpenUrlCrossRefPubMedWeb of Science
  113. 113.↵
    K. Wang, M. Li, H. Hakonarson, ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38, e164 (2010).
    OpenUrlCrossRefPubMed
  114. 114.↵
    S. Lee, G. R. Abecasis, M. Boehnke, X. Lin, Rare-variant association analysis: study designs and statistical tests. American journal of human genetics 95, 5–23 (2014).
    OpenUrlCrossRefPubMed

References

  1. 1.
    D. M. Ruderfer et al., Genomic dissection of bipolar disorder and schizophrenia, including 28 subphenotypes. Cell 173, 1705-1715.e1716 (2018).
    OpenUrlCrossRefPubMed
Back to top
Posted January 16, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Kleine Levin syndrome is associated with birth difficulties and genetic variants in the TRANK1 gene loci
Aditya Ambati, Ryan Hillary, Smaranda Leu-Semenescu, Hanna M. Ollila, Ling Lin, Emmanuel During, Neal Farber, Thomas J Rico, Juliette Faraco, Eileen Leary, Andrea Goldstein-Piekarski, Yu-Shu Huang, Fang Han, Yakov Sivan, Michel Lecendreux, Pauline Dodet, Makoto Honda, Natan Gadoth, Sona Nevsimalova, Fabio Pizza, Takashi Kanbayashi, Rosa Peraita Adrados, Guy Leschziner, Rosa Hasan, Francesca Canellas, Kazuhiko Kume, Makrina Daniilidou, Patrice Bourgin, David Rye, José L Vicario, Birgit Högl, Seung Chul Hong, Guiseppe Plazzi, Geert Mayer, Anne Marie Landtblom, Yves Dauvilliers, Isabelle Arnulf, Emmanuel Mignot
medRxiv 2021.01.08.20249006; doi: https://doi.org/10.1101/2021.01.08.20249006
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area