Hiromi Toyoda

High-accuracy imputation for HLA class I and II genes based on high-resolution SNP data of population-specific references.

Statistical imputation of classical human leukocyte antigen (HLA) alleles is becoming an indispensable tool for fine-mappings of disease association signals from case–control genome-wide association studies. However, most currently available HLA imputation tools are based on European reference populations and are not suitable for direct application to non-European populations. Among the HLA imputation tools, The HIBAG R package is a flexible HLA imputation tool that is equipped with a wide range of population-based classifiers; moreover, HIBAG R enables individual researchers to build custom classifiers. Here, two data sets, each comprising data from healthy Japanese individuals of difference sample sizes, were used to build custom classifiers. HLA imputation accuracy in five HLA classes (HLA-A, HLA-B, HLA-DRB1, HLA-DQB1 and HLA-DPB1) increased from the 82.5–98.8% obtained with the original HIBAG references to 95.2–99.5% with our custom classifiers. A call threshold (CT) of 0.4 is recommended for our Japanese classifiers; in contrast, HIBAG references recommend a CT of 0.5. Finally, our classifiers could be used to identify the risk haplotypes for Japanese narcolepsy with cataplexy, HLA-DRB1*15:01 and HLA-DQB1*06:02, with 100% and 99.7% accuracy, respectively; therefore, these classifiers can be used to supplement the current lack of HLA genotyping data in widely available genome-wide association study data sets.

New susceptibility variants to narcolepsy identified in HLA class II region

Narcolepsy, a sleep disorder characterized by excessive daytime sleepiness, cataplexy and rapid eye movement sleep abnormalities, is tightly associated with human leukocyte antigen HLA-DQB1*06:02. DQB1*06:02 is common in the general population (10-30%); therefore, additional genetic factors are needed for the development of narcolepsy. In the present study, HLA-DQB1 in 664 Japanese narcoleptic subjects and 3131 Japanese control subjects was examined to determine whether HLA-DQB1 alleles located in trans of DQB1*06:02 are associated with narcolepsy. The strongest association was with DQB1*06:01 (P = 1.4 × 10(-10), odds ratio, OR = 0.39), as reported in previous studies. Additional predisposing effects of DQB1*03:02 were also found (P = 2.5 × 10(-9), OR = 1.97). A comparison between DQB1*06:02 heterozygous cases and controls revealed dominant protective effects of DQB1*06:01 and DQB1*05:01. In addition, a single-nucleotide polymorphism-based conditional analysis controlling for the effect of HLA-DQB1 was performed to determine whether there were other independent HLA associations outside of HLA-DQB1. This analysis revealed associations at HLA-DPB1 in the HLA class II region (rs3117242, P = 4.1 × 10(-5), OR = 2.45; DPB1*05:01, P = 8.1 × 10(-3), OR = 1.39). These results indicate that complex HLA class II associations contribute to the genetic predisposition to narcolepsy.

A polymorphism in CCR1/CCR3 is associated with narcolepsy.

Etiology of narcolepsy-cataplexy involves multiple genetic and environmental factors. While the human leukocyte antigen (HLA)-DRB1*15:01-DQB1*06:02 haplotype is strongly associated with narcolepsy, it is not sufficient for disease development. To identify additional, non-HLA susceptibility genes, we conducted a genome-wide association study (GWAS) using Japanese samples. An initial sample set comprising 409 cases and 1562 controls was used for the GWAS of 525,196 single nucleotide polymorphisms (SNPs) located outside the HLA region. An independent sample set comprising 240 cases and 869 controls was then genotyped at 37 SNPs identified in the GWAS. We found that narcolepsy was associated with a SNP in the promoter region of chemokine (C-C motif) receptor 1 (CCR1) (rs3181077, P=1.6×10(-5), odds ratio [OR]=1.86). This rs3181077 association was replicated with the independent sample set (P=0.032, OR=1.36). We measured mRNA levels of candidate genes in peripheral blood samples of 38 cases and 37 controls. CCR1 and CCR3 mRNA levels were significantly lower in patients than in healthy controls, and CCR1 mRNA levels were associated with rs3181077 genotypes. In vitro chemotaxis assays were also performed to measure monocyte migration. We observed that monocytes from carriers of the rs3181077 risk allele had lower migration indices with a CCR1 ligand. CCR1 and CCR3 are newly discovered susceptibility genes for narcolepsy. These results highlight the potential role of CCR genes in narcolepsy and support the hypothesis that patients with narcolepsy have impaired immune function.

Increased plasma IL-6, IL-8, TNF-alpha, and G-CSF in Japanese narcolepsy

Narcolepsy is a chronic hypersomnia involving excessive daytime sleepiness and cataplexy. Some susceptibility genes and environmental factors suggest that post-infectious immunological alterations underlie its pathophysiology. To investigate the immunological alterations in narcolepsy patients, we examined cytokines. Nine healthy controls and twenty-one narcolepsy patients with cataplexy were studied. All subjects were positive for the HLA-DRB1(∗)1501-DQB1(∗)0602 allele. Age-, sex-, and body mass index -matched healthy controls were selected. Plasma samples were separated using EDTA-2K-coated blood collection tubes. Bioplex Pro Human Cytokine 17-Plex Assays were used to measure plasma cytokines. Elevations of interleukin (IL)-6, IL-8, granulocyte- colony stimulating factor (G-CSF), and tumor necrosis factor-alpha were found in the narcolepsy group compared with healthy controls (p<0.05). G-CSF values were significantly correlated with the disease duration in narcolepsy patients (r=0.426, p<0.05). IL-8 and G-CSF play major roles in neutrophil activation in respiratory diseases. Since environmental factors including infection are reportedly associated with narcolepsy onset, elevated IL-8 and G-CSF may be involved in the pathophysiology of narcolepsy.

Genome-wide association study of HLA-DQB1*06:02 negative essential hypersomnia

Essential hypersomnia (EHS), a sleep disorder characterized by excessive daytime sleepiness, can be divided into two broad classes based on the presence or absence of the HLA-DQB1*06:02 allele. HLA-DQB1*06:02-positive EHS and narcolepsy with cataplexy are associated with the same susceptibility genes. In contrast, there are fewer studies of HLA-DQB1*06:02 negative EHS which, we hypothesized, involves a different pathophysiological pathway than does narcolepsy with cataplexy. In order to identify susceptibility genes associated with HLA-DQB1*06:02 negative EHS, we conducted a genome-wide association study (GWAS) of 125 unrelated Japanese EHS patients lacking the HLA-DQB1*06:02 allele and 562 Japanese healthy controls. A comparative study was also performed on 268 HLA-DQB1*06:02 negative Caucasian hypersomnia patients and 1761 HLA-DQB1*06:02 negative Caucasian healthy controls. We identified three SNPs that each represented a unique locus— rs16826005 (P = 1.02E-07; NCKAP5), rs11854769 (P = 6.69E-07; SPRED1), and rs10988217 (P = 3.43E-06; CRAT) that were associated with an increased risk of EHS in this Japanese population. Interestingly, rs10988217 showed a similar tendency in its association with both HLA-DQB1*06:02 negative EHS and narcolepsy with cataplexy in both Japanese and Caucasian populations. This is the first GWAS of HLA-DQB1*06:02 negative EHS, and the identification of these three new susceptibility loci should provide additional insights to the pathophysiological pathway of this condition.

Transcriptional regulation of the hypocretin/orexin gene by NR6A1.

The hypocretin (also known as orexin) neuropeptide system coordinates the regulation of various physiological processes. A reduction in Nr6a1 expression was observed in hypocretin neuron-ablated transgenic mice. To show that prepro-hypocretin transcription is functionally modulated by NR6A1, we performed chromatin immunoprecipitation (ChIP) analysis, double-immunostaining, a luciferase reporter assay, and an in utero electroporation study. ChIP analysis showed that endogenous NR6A1 binds to a putative NR6A1-binding site. Double-immunostaining indicated almost all hypocretin neurons were positive for NR6A1 immunoreactivity. NR6A1 overexpression in SH-SY5Y cells modulated hypocretin promoter activity, an effect that was countered by lacking a putative NR6A1-binding site. Electroporation with Nr6a1 in the foetal hypothalamus promoted hypocretin transcription as compared to GFP-electroporation. These experiments confirmed that NR6A1 works as a regulator for hypocretin transcription.

Evaluation of polygenic risks for narcolepsy and essential hypersomnia

In humans, narcolepsy is a sleep disorder that is characterized by sleepiness, cataplexy and rapid eye movement (REM) sleep abnormalities. Essential hypersomnia (EHS) is another type of sleep disorder that is characterized by excessive daytime sleepiness without cataplexy. A human leukocyte antigen (HLA) class II allele, HLA-DQB1*06:02, is a major genetic factor for narcolepsy. Almost all narcoleptic patients are carriers of this HLA allele, while 30–50% of EHS patients and 12% of all healthy individuals in Japan carry this allele. The pathogenesis of narcolepsy and EHS is thought to be partially shared. To evaluate the contribution of common single-nucleotide polymorphisms (SNPs) to narcolepsy onset and to assess the common genetic background of narcolepsy and EHS, we conducted a polygenic analysis that included 393 narcoleptic patients, 38 EHS patients with HLA-DQB1*06:02, 119 EHS patients without HLA-DQB1*06:02 and 1582 healthy individuals. We also included 376 individuals with panic disorder and 213 individuals with autism to confirm whether the results were biased. Polygenic risks in narcolepsy were estimated to explain 58.1% (PHLA-DQB1*06:02=2.30 × 10−48, Pwhole genome without HLA-DQB1*06:02=6.73 × 10−2) including HLA-DQB1*06:02 effects and 1.3% (Pwhole genome without HLA-DQB1*06:02=2.43 × 10−2) excluding HLA-DQB1*06:02 effects. The results also indicated that small-effect SNPs contributed to the development of narcolepsy. Reported susceptibility SNPs for narcolepsy in the Japanese population, CPT1B (carnitine palmitoyltransferase 1B), TRA@ (T-cell receptor alpha) and P2RY11 (purinergic receptor P2Y, G-protein coupled, 11), were found to explain 0.8% of narcolepsy onset (Pwhole genome without HLA-DQB1*06:02=9.74 × 10−2). EHS patients with HLA-DQB1*06:02 were estimated to have higher shared genetic background to narcoleptic patients than EHS patients without HLA-DQB1*06:02 even when the effects of HLA-DQB1*06:02 were excluded (EHS with HLA-DQB1*06:02: 40.4%, PHLA-DQB1*06:02=7.02 × 10−14, Pwhole genome without HLA-DQB1*06:02=1.34 × 10−1, EHS without HLA-DQB1*06:02: 0.4%, Pwhole genome without HLA-DQB1*06:02=3.06 × 10−1). Meanwhile, the polygenic risks for narcolepsy could not explain the onset of panic disorder and autism, suggesting that our results were reasonable.

An association analysis of HLA-DQB1 with narcolepsy without cataplexy and idiopathic hypersomnia with/without long sleep time in a Japanese population

Narcolepsy without cataplexy (NA w/o CA) (narcolepsy type 2) is a lifelong disorder characterized by excessive daytime sleepiness and rapid eye movement (REM) sleep abnormalities, but no cataplexy. In the present study, we examined the human leukocyte antigen HLA-DQB1 in 160 Japanese patients with NA w/o CA and 1,418 control subjects. Frequencies of DQB1*06:02 were significantly higher in patients with NA w/o CA compared with controls (allele frequency: 16.6 vs. 7.8%, P=1.1×10−7, odds ratio (OR)=2.36; carrier frequency: 31.3 vs. 14.7%, P=7.6×10−8, OR=2.64). Distributions of HLA-DQB1 alleles other than DQB1*06:02 were compared between NA w/o CA and narcolepsy with cataplexy (NA-CA) to assess whether the genetic backgrounds of the two diseases have similarities. The distribution of the HLA-DQB1 alleles in DQB1*06:02-negative NA w/o CA was significantly different from that in NA-CA (P=5.8×10−7). On the other hand, the patterns of the HLA-DQB1 alleles were similar between DQB1*06:02-positive NA w/o CA and NA-CA. HLA-DQB1 analysis was also performed in 186 Japanese patients with idiopathic hypersomnia (IHS) with/without long sleep time, but no significant associations were observed.

Genome-wide analysis of CNV (copy number variation) and their associations with narcolepsy in a Japanese population.

In humans, narcolepsy with cataplexy (narcolepsy) is a sleep disorder that is characterized by sleepiness, cataplexy and rapid eye movement (REM) sleep abnormalities. Narcolepsy is caused by a reduction in the number of neurons that produce hypocretin (orexin) neuropeptide. Both genetic and environmental factors contribute to the development of narcolepsy.Rare and large copy number variations (CNVs) reportedly play a role in the etiology of a number of neuropsychiatric disorders. Narcolepsy is considered a neurological disorder; therefore, we sought to investigate any possible association between rare and large CNVs and human narcolepsy. We used DNA microarray data and a CNV detection software application, PennCNV-Affy, to detect CNVs in 426 Japanese narcoleptic patients and 562 healthy individuals. Overall, we found a significant enrichment of rare and large CNVs (frequency ⩽1%, size ⩾100 kb) in the patients (case–control ratio of CNV count=1.54, P=5.00 × 10−4). Next, we extended a region-based association analysis by including CNVs with its size ⩾30 kb. Rare and large CNVs in PARK2 region showed a significant association with narcolepsy. Four patients were assessed to carry duplications of the gene region, whereas no controls carried the duplication, which was further confirmed by quantitative PCR assay. This duplication was also found in 2 essential hypersomnia (EHS) patients out of 171 patients. Furthermore, a pathway analysis revealed enrichments of gene disruptions by rare and large CNVs in immune response, acetyltransferase activity, cell cycle regulation and regulation of cell development. This study constitutes the first report on the risk association between multiple rare and large CNVs and the pathogenesis of narcolepsy. In the future, replication studies are needed to confirm the associations.

Narcolepsy susceptibility gene CCR3 modulates sleep-wake patterns in mice

Narcolepsy is caused by the loss of hypocretin (Hcrt) neurons and is associated with multiple genetic and environmental factors. Although abnormalities in immunity are suggested to be involved in the etiology of narcolepsy, no decisive mechanism has been established. We previously reported chemokine (C-C motif) receptor 3 (CCR3) as a novel susceptibility gene for narcolepsy. To understand the role of CCR3 in the development of narcolepsy, we investigated sleep-wake patterns of Ccr3 knockout (KO) mice. Ccr3 KO mice exhibited fragmented sleep patterns in the light phase, whereas the overall sleep structure in the dark phase did not differ between Ccr3 KO mice and wild-type (WT) littermates. Intraperitoneal injection of lipopolysaccharide (LPS) promoted wakefulness and suppressed both REM and NREM sleep in the light phase in both Ccr3 KO and WT mice. Conversely, LPS suppressed wakefulness and promoted NREM sleep in the dark phase in both genotypes. After LPS administration, the proportion of time spent in wakefulness was higher, and the proportion of time spent in NREM sleep was lower in Ccr3 KO compared to WT mice only in the light phase. LPS-induced changes in sleep patterns were larger in Ccr3 KO compared to WT mice. Furthermore, we quantified the number of Hcrt neurons and found that Ccr3 KO mice had fewer Hcrt neurons in the lateral hypothalamus compared to WT mice. We found abnormalities in sleep patterns in the resting phase and in the number of Hcrt neurons in Ccr3 KO mice. These observations suggest a role for CCR3 in sleep-wake regulation in narcolepsy patients.