Identification of region-specific gene isoforms in the human brain using long-read transcriptome sequencing.

In neurological and neuropsychiatric diseases, different brain regions are affected, and differences in gene expression patterns could potentially explain this mechanism. However, limited studies have precisely explored gene expression in different regions of the human brain. In this study, we performed long-read RNA sequencing on three different brain regions of the same individuals: the cerebellum, hypothalamus, and temporal cortex. Despite stringent filtering criteria excluding isoforms predicted to be artifacts, over half of the isoforms expressed in multiple samples across multiple regions were found to be unregistered in the GENCODE reference. We then especially focused on genes with different major isoforms in each brain region, even with similar overall expression levels, and identified that many of such genes including GAS7 might have distinct roles in dendritic spine and neuronal formation in each region. We also found that DNA methylation might, in part, drive different isoform expressions in different regions. These findings highlight the significance of analyzing isoforms expressed in disease-relevant sites.

Narcolepsy type I-associated DNA methylation and gene expression changes in the human leukocyte antigen region.

Narcolepsy type 1 (NT1) is caused by a loss of hypothalamic orexin-producing cells, and autoreactive CD4+ and CD8+ T cells have been suggested to play a role in the autoimmune mechanism. Although NT1 showed a strong association with human leukocyte antigen (HLA)-DQB1*06:02, the responsible antigens remain unidentified. We analyzed array-based DNA methylation and gene expression data for the HLA region in CD4+ and CD8+ T cells that were separated from the peripheral blood mononuclear cells of Japanese subjects (NT1, N = 42; control, N = 42). As the large number of SNPs in the HLA region might interfere with the affinity of the array probes, we conducted a comprehensive assessment of the reliability of each probe. The criteria were based on a previous study reporting that the presence of frequent SNPs, especially on the 3' side of the probe, makes the probe unreliable. We confirmed that 90.3% of the probes after general filtering in the HLA region do not include frequent SNPs, and are thus suitable for analysis, particularly in Japanese subjects. We then performed an association analysis, and found that several CpG sites in the HLA class II region of the patients were significantly hypomethylated in CD4+ and CD8+ T cells. This association was not detected when the effect of HLA-DQB1*06:02 was considered, suggesting that the hypomethylation was possibly derived from HLA-DQB1*06:02. Further RNA sequencing revealed reduced expression levels of HLA-DQB1 alleles other than HLA-DQB1*06:02 in the patients with NT1. Our results suggest the involvement of epigenetic and expressional changes in HLA-DQB1 in the pathogenesis of NT1.

Low carnitine palmitoyltransferase 1 activity is a risk factor for narcolepsy type 1 and other hypersomnia.

STUDY OBJECTIVES: Narcolepsy type 1 (NT1) is associated with metabolic abnormalities but their etiology remains largely unknown. The gene for carnitine palmitoyltransferase 1B (CPT1B) and abnormally low serum acylcarnitine levels have been linked to NT1. To elucidate the details of altered fatty acid metabolism, we determined levels of individual acylcarnitines and evaluated CPT1 activity in patients with NT1 and other hypersomnia. METHODS: Blood samples from 57 NT1, 51 other hypersomnia patients, and 61 healthy controls were analyzed. The levels of 25 major individual acylcarnitines were determined and the C0/(t[C16] + t[C18]) ratio was used as a CPT1 activity marker. We further performed transcriptome analysis using independent blood samples from 42 NT1 and 42 healthy controls to study the relevance of fatty acid metabolism. NT1-specific changes in CPT1 activity and in expression of related genes were investigated. RESULTS: CPT1 activity was lower in patients with NT1 (p = 0.00064) and other hypersomnia (p = 0.0014) than in controls. Regression analysis revealed that CPT1 activity was an independent risk factor for NT1 (OR: 1.68; p = 0.0031) and for other hypersomnia (OR: 1.64; p = 0.0042). There was a significant interaction between obesity (BMI <25, ≥25) and the SNP rs5770917 status such that nonobese NT1 patients without risk allele had better CPT1 activity (p = 0.0089). The expression levels of carnitine-acylcarnitine translocase (CACT) and CPT2 in carnitine shuttle were lower in NT1 (p = 0.000051 and p = 0.00014, respectively). CONCLUSIONS: These results provide evidences that abnormal fatty acid metabolism is involved in the pathophysiology of NT1 and other hypersomnia.

A rare genetic variant in the cleavage site of prepro-orexin is associated with idiopathic hypersomnia.

Idiopathic hypersomnia (IH) is a rare, heterogeneous sleep disorder characterized by excessive daytime sleepiness. In contrast to narcolepsy type 1, which is a well-defined type of central disorders of hypersomnolence, the etiology of IH is poorly understood. No susceptibility loci associated with IH have been clearly identified, despite the tendency for familial aggregation of IH. We performed a variation screening of the prepro-orexin/hypocretin and orexin receptors genes and an association study for IH in a Japanese population, with replication (598 patients and 9826 controls). We identified a rare missense variant (g.42184347T>C; p.Lys68Arg; rs537376938) in the cleavage site of prepro-orexin that was associated with IH (minor allele frequency of 1.67% in cases versus 0.32% in controls, P = 2.7 × 10-8, odds ratio = 5.36). Two forms of orexin (orexin-A and -B) are generated from cleavage of one precursor peptide, prepro-orexin. The difference in cleavage efficiency between wild-type (Gly-Lys-Arg; GKR) and mutant (Gly-Arg-Arg; GRR) peptides was examined by assays using proprotein convertase subtilisin/kexin (PCSK) type 1 and PCSK type 2. In both PCSK1 and PCSK2 assays, the cleavage efficiency of the mutant peptide was lower than that of the wild-type peptide. We also confirmed that the prepro-orexin peptides themselves transmitted less signaling through orexin receptors than mature orexin-A and orexin-B peptides. These results indicate that a subgroup of IH is associated with decreased orexin signaling, which is believed to be a hallmark of narcolepsy type 1.

A variant in orexin receptor-2 is associated with self-reported daytime sleepiness in the Japanese population.

Excessive daytime sleepiness is characterized by a persistent feeling of having trouble staying awake, typically with inappropriate sleep episodes. Orexin (hypocretin) is a neuropeptide that regulates sleep-wake cycles and rapid eye movement sleep. Several large-scale genome-wide association studies (GWASs) in European populations have found genetic variants in orexin receptor-1 (OX1R) and -2 (OX2R) that are associated with sleep traits including daytime sleepiness. To identify genetic variants associated with daytime sleepiness, we performed an association study of genetic variants in prepro-orexin, OX1R, and OX2R in 14,329 Japanese individuals from the Tohoku Medical Megabank Project cohort. A genetic variant in OX2R was significantly associated with self-reported daytime sleepiness after Bonferroni correction (rs188018846: P = 8.4E-05). In addition, a missense variant in OX2R identified by the European GWASs showed a nominally significant association with daytime sleepiness in a Japanese population (p.Ile308Val, rs2653349: P = 0.044). Multiple genetic variants in OX2R can affect daytime sleepiness in general populations.

Genome-wide association study of idiopathic hypersomnia in a Japanese population.

Idiopathic hypersomnia (IH) is a rare sleep disorder characterized by excessive daytime sleepiness, great difficulty upon awakening, and prolonged sleep time. In contrast to narcolepsy type 1, which is a well-recognized hypersomnia, the etiology of IH remains poorly understood. No susceptibility loci for IH have been identified, although familial aggregations have been observed among patients with IH. Narcolepsy type 1 is strongly associated with human leukocyte antigen (HLA)-DQB1*06:02; however, no significant associations between IH and HLA alleles have been reported. To identify genetic variants that affect susceptibility to IH, we performed a genome-wide association study (GWAS) and two replication studies involving a total of 414 Japanese patients with IH and 6587 healthy Japanese individuals. A meta-analysis of the three studies found no single-nucleotide polymorphisms (SNPs) that reached the genome-wide significance level. However, we identified several candidate SNPs for IH. For instance, a common genetic variant (rs2250870) within an intron of PDE9A was suggestively associated with IH. rs2250870 was significantly associated with expression levels of PDE9A in not only whole blood but also brain tissues. The leading SNP in the PDE9A region was the same in associations with both IH and PDE9A expression. PDE9A is a potential target in the treatment of several brain diseases, such as depression, schizophrenia, and Alzheimer's disease. It will be necessary to examine whether PDE9A inhibitors that have demonstrated effects on neurophysiologic and cognitive function can contribute to the development of new treatments for IH, as higher expression levels of PDE9A were observed with regard to the risk allele of rs2250870. The present study constitutes the first GWAS of genetic variants associated with IH. A larger replication study will be required to confirm these associations. Supplementary Information: The online version contains supplementary material available at 10.1007/s41105-021-00349-2.

Metabolome analysis using cerebrospinal fluid from narcolepsy type 1 patients.

Narcolepsy type 1 (NT1) is a hypersomnia characterized by excessive daytime sleepiness and cataplexy. Inappropriate regulation of fatty acid metabolism has been suggested to be involved in the pathophysiology of NT1, but the detailed mechanisms remain uncertain. Here we performed a metabolomic analysis of cerebrospinal fluid samples from 14 NT1 and 17 control subjects using a novel capillary electrophoresis coupled with Fourier transform mass spectrometry. A total of 268 metabolites were identified and the amount of histidine was the most significantly increased in NT1 patients (p = 4.0 × 10-4). Validation analysis using high-performance liquid chromatography (HPLC) including independent replication samples also identified the association of histidine (p = 2.02 × 10-3). Further, levels of histamine, which is synthesized from histidine, were also examined using HPLC and were found to be significantly decreased in NT1 patients (p = 6.12 × 10-4). Pathway analysis with nominally significant metabolites identified several pathways related to the metabolism of glycogenic amino acids, suggesting that glycogenesis is enhanced in NT1 as a compensatory mechanism for fatty acid metabolism. We performed further exploratory analysis, searching for metabolites associated with sleep variables from polysomnography and the multiple sleep latency test. As a result, 5'-deoxy-5'-methylthioadenosine showed a significant association with apnea-hypopnea index (p = 2.66 ×10-6). Moreover, gamma aminobutyric acid displayed a negative correlation with rapid eye movement sleep latency (REML), and thus might represent an intriguing target for future studies to elucidate how the controlling circuit of REM sleep is associated with abnormally short REML in NT1.

Case-control association study of rare nonsynonymous variants of SCN1A and KCNQ2 in acute encephalopathy with biphasic seizures and late reduced diffusion.

PURPOSE: Acute encephalopathy with biphasic seizures and late reduced diffusion (AESD) is characterized by prolonged febrile seizures at onset and subsequent damage to the cerebral cortex of infants and children. The pathogenesis is suspected to be excitotoxicity leading to neuronal death. SCN1A and KCNQ2 are causative genes of genetic epilepsy including Dravet syndrome and Ohtahara syndrome. Here we conducted a case-control rare-variant association study of the two genes in AESD. METHODS: The coding regions of SCN1A and KCNQ2 were sequenced by the Sanger method for 175 and 111 patients, respectively, with AESD. As control subjects, we used genetic data from 3554 subjects provided by the Integrative Japanese Genome Variation Database (iJGVD). Then we performed a case-control association study of rare missense and splice region variants (minor allele frequency < 0.005) of each gene with AESD using Weighted Sum Statistics (WSS) and Sequence Kernel Association Test (SKAT). RESULTS: SCN1A rare variants had a significant association with AESD after correction for multiple tests (WSS, permutated p value 4.00 × 10-3: SKAT, p value 2.51 × 10-4). The association was more significant when we focused on deleterious variants (WSS, permutated p = 9.00 × 10-4; SKAT, p = 4.99 × 10-5). Although KCNQ2 rare nonsynonymous variants tended to be more frequent in patients than in controls, there was no significant difference. CONCLUSION: Our study provided statistical evidence of an association between SCN1A and AESD for the first time, and established SCN1A as one of the susceptibility genes for AESD.

Sensitivity to gene dosage and gene expression affects genes with copy number variants observed among neuropsychiatric diseases.

BACKGROUND: Copy number variants (CNVs) have been reported to be associated with diseases, traits, and evolution. However, it is hard to determine which gene should have priority as a target for further functional experiments if a CNV is rare or a singleton. In this study, we attempted to overcome this issue by using two approaches: by assessing the influences of gene dosage sensitivity and gene expression sensitivity. Dosage sensitive genes derived from two-round whole-genome duplication in previous studies. In addition, we proposed a cross-sectional omics approach that utilizes open data from GTEx to assess the effect of whole-genome CNVs on gene expression. METHODS: Affymetrix Genome-Wide SNP Array 6.0 was used to detect CNVs by PennCNV and CNV Workshop. After quality controls for population stratification, family relationship and CNV detection, 287 patients with narcolepsy, 133 patients with essential hypersomnia, 380 patients with panic disorders, 164 patients with autism, 784 patients with Alzheimer disease and 1280 healthy individuals remained for the enrichment analysis. RESULTS: Overall, significant enrichment of dosage sensitive genes was found across patients with narcolepsy, panic disorders and autism. Particularly, significant enrichment of dosage-sensitive genes in duplications was observed across all diseases except for Alzheimer disease. For deletions, less or no enrichment of dosage-sensitive genes with deletions was seen in the patients when compared to the healthy individuals. Interestingly, significant enrichments of genes with expression sensitivity in brain were observed in patients with panic disorder and autism. While duplications presented a higher burden, deletions did not cause significant differences when compared to the healthy individuals. When we assess the effect of sensitivity to genome dosage and gene expression at the same time, the highest ratio of enrichment was observed in the group including dosage-sensitive genes and genes with expression sensitivity only in brain. In addition, shared CNV regions among the five neuropsychiatric diseases were also investigated. CONCLUSIONS: This study contributed the evidence that dosage-sensitive genes are associated with CNVs among neuropsychiatric diseases. In addition, we utilized open data from GTEx to assess the effect of whole-genome CNVs on gene expression. We also investigated shared CNV region among neuropsychiatric diseases.

Epigenome-wide association study of narcolepsy-affected lateral hypothalamic brains, and overlapping DNA methylation profiles between narcolepsy and multiple sclerosis.

Narcolepsy with cataplexy is a sleep disorder caused by a deficiency in hypocretin neurons in the lateral hypothalamus (LH). Here we performed an epigenome-wide association study (EWAS) of DNA methylation for narcolepsy and replication analyses using DNA samples extracted from two brain regions: LH (Cases: N = 4; Controls: N = 4) and temporal cortex (Cases: N = 7; Controls: N = 7). Seventy-seven differentially methylated regions (DMRs) were identified in the LH analysis, with the top association of a DMR in the myelin basic protein (MBP) region. Only five DMRs were detected in the temporal cortex analysis. Genes annotated to LH DMRs were significantly associated with pathways related to fatty acid response or metabolism. Two additional analyses applying the EWAS data were performed: (1) investigation of methylation profiles shared between narcolepsy and other disorders and (2) an integrative analysis of DNA methylation data and a genome-wide association study for narcolepsy. The results of the two approaches, which included significant overlap of methylated positions associated with narcolepsy and multiple sclerosis, indicated that the two diseases may partly share their pathogenesis. In conclusion, DNA methylation in LH where loss of orexin-producing neurons occurs may play a role in the pathophysiology of the disease.

A missense variant in PER2 is associated with delayed sleep-wake phase disorder in a Japanese population.

Delayed sleep-wake phase disorder (DSWPD) is a subtype of circadian rhythm sleep-wake disorders, and is characterized by an inability to fall asleep until late at night and wake up at a socially acceptable time in the morning. The study aim was to identify low-frequency nonsense and missense variants that are associated with DSWPD. Candidate variants in circadian rhythm-related genes were extracted by integration of genetic variation databases and in silico assessment. We narrowed down the candidates to six variants. To examine whether the six variants are associated with DSWPD, we performed an association study in 236 Japanese patients with DSWPD and 1436 controls. A low-frequency missense variant (p.Val1205Met) in PER2 showed a significant association with DSWPD (2.5% in cases and 1.1% in controls, P = 0.026, odds ratio (OR) = 2.32). The variant was also associated with idiopathic hypersomnia known to have a tendency toward phase delay (P = 0.038, OR = 2.07). PER2 forms a heterodimer with CRY, and the heterodimer plays an important role in the regulation of circadian rhythms. Val1205 is located in the CRY-binding domain of PER2 and was hypothesized to interact with CRY. The p.Val1205Met substitution could be a potential genetic marker for DSWPD.

Thermolabile polymorphism of carnitine palmitoyltransferase 2: A genetic risk factor of overall acute encephalopathy.

OBJECTIVES: Acute encephalopathy is an acute brain dysfunction after preceding infection, consisting of multiple syndromes. Some syndromes, such as acute encephalopathy with biphasic seizures and late reduced diffusion (AESD), are severe with poor outcome, whereas others, such as clinically mild encephalitis/encephalopathy with reversible splenial lesion (MERS), are mild with favorable outcome. Previous study reported the association of the thermolabile polymorphism in Carnitine Palmitoyltransferase 2 (CPT2) gene and severe syndromes of acute encephalopathy. To further explore the pathogenetic role of CPT2 in acute encephalopathy, we conducted a case-control association study of a typical thermolabile CPT2 polymorphism, rs2229291, in 416 patients of acute encephalopathy, including both severe and mild syndromes. METHODS: The case cohort consisted of 416 patients, including AESD, MERS, and other syndromes. The control subjects were 100 healthy Japanese. rs2229291 was genotyped by Sanger sequencing. Genetic distribution was compared between the patients and controls using Cochran-Armitage trend test. RESULTS: Minor allele frequency of rs2229291 was significantly higher in AESD (p = 0.044), MERS (p = 0.015) and entire acute encephalopathy (p = 0.044) compared to the controls. The polymorphism showed no significant association with influenza virus, or with outcome. CONCLUSIONS: This study provided evidence that CPT2 is a susceptibility gene for overall acute encephalopathy, including both severe and mild syndromes, and suggested that impairment of mitochondrial metabolism is common to various syndromes of acute encephalopathy.

Genetics of narcolepsy.

Narcolepsy is a term that was initially coined by Gélineáu in 1880 and is a chronic neurological sleep disorder that manifests as a difficulty in maintaining wakefulness and sleep for long periods. Currently, narcolepsy is subdivided into two types according to the International Classification of Sleep Disorders, 3rd edition: narcolepsy type 1 (NT1) and narcolepsy type 2 (NT2). NT1 is characterized by excessive daytime sleepiness, cataplexy, hypnagogic hallucinations, and sleep paralysis and is caused by a marked reduction in neurons in the hypothalamus that produce orexin (hypocretin), which is a wakefulness-associated neuropeptide. Except for cataplexy, NT2 exhibits most of the same symptoms as NT1. NT1 is a multifactorial disease, and genetic variations at multiple loci are associated with NT1. Almost all patients with NT1 carry the specific human leukocyte antigen (HLA) allele HLA-DQB1 * 06:02. Genome-wide association studies have uncovered >10 genomic variations associated with NT1. Rare variants associated with NT1 have also been identified by DNA genome sequencing. NT2 is also a complex disorder, but its underlying genetic architecture is poorly understood. However, several studies have revealed loci that increase susceptibility to NT2. The currently identified loci cannot explain the heritability of narcolepsy (NT1 and NT2). We expect that future genomic research will provide important contributions to our understanding of the genetic basis and pathogenesis of narcolepsy.

A variant at 9q34.11 is associated with HLA-DQB1*06:02 negative essential hypersomnia.

Essential hypersomnia (EHS) is a lifelong disorder characterized by excessive daytime sleepiness without cataplexy. EHS is associated with human leukocyte antigen (HLA)-DQB1*06:02, similar to narcolepsy with cataplexy (narcolepsy). Previous studies suggest that DQB1*06:02-positive and -negative EHS are different in terms of their clinical features and follow different pathological pathways. DQB1*06:02-positive EHS and narcolepsy share the same susceptibility genes. In the present study, we report a genome-wide association study with replication for DQB1*06:02-negative EHS (408 patients and 2247 healthy controls, all Japanese). One single-nucleotide polymorphism, rs10988217, which is located 15-kb upstream of carnitine O-acetyltransferase (CRAT), was significantly associated with DQB1*06:02-negative EHS (P = 7.5 × 10-9, odds ratio = 2.63). The risk allele of the disease-associated SNP was correlated with higher expression levels of CRAT in various tissues and cell types, including brain tissue. In addition, the risk allele was associated with levels of succinylcarnitine (P = 1.4 × 10-18) in human blood. The leading SNP in this region was the same in associations with both DQB1*06:02-negative EHS and succinylcarnitine levels. The results suggest that DQB1*06:02-negative EHS may be associated with an underlying dysfunction in energy metabolic pathways.

Epigenome-wide association study of DNA methylation in narcolepsy: an integrated genetic and epigenetic approach.

Narcolepsy with cataplexy, which is a hypersomnia characterized by excessive daytime sleepiness and cataplexy, is a multifactorial disease caused by both genetic and environmental factors. Several genetic factors including HLA-DQB1*06:02 have been identified; however, the disease etiology is still unclear. Epigenetic modifications, such as DNA methylation, have been suggested to play an important role in the pathogenesis of complex diseases. Here, we examined DNA methylation profiles of blood samples from narcolepsy and healthy control individuals and performed an epigenome-wide association study (EWAS) to investigate methylation loci associated with narcolepsy. Moreover, data from the EWAS and a previously performed narcolepsy genome-wide association study were integrated to search for methylation loci with causal links to the disease. We found that (1) genes annotated to the top-ranked differentially methylated positions (DMPs) in narcolepsy were associated with pathways of hormone secretion and monocarboxylic acid metabolism. (2) Top-ranked narcolepsy-associated DMPs were significantly more abundant in non-CpG island regions and more than 95 per cent of such sites were hypomethylated in narcolepsy patients. (3) The integrative analysis identified the CCR3 region where both a single methylation site and multiple single-nucleotide polymorphisms were found to be associated with the disease as a candidate region responsible for narcolepsy. The findings of this study suggest the importance of future replication studies, using methylation technologies with wider genome coverage and/or larger number of samples, to confirm and expand on these results.

An epigenome-wide methylation study of healthy individuals with or without depressive symptoms.

Major depressive disorder is a common psychiatric disorder that is thought to be triggered by both genetic and environmental factors. Depressive symptoms are an important public health problem and contribute to vulnerability to major depression. Although a substantial number of genetic and epigenetic studies have been performed to date, the detailed etiology of depression remains unclear and there are no validated biomarkers. DNA methylation is one of the major epigenetic modifications that play diverse roles in the etiology of complex diseases. In this study, we performed an epigenome-wide association study (EWAS) of DNA methylation on subjects with (N = 20) or without (N = 27) depressive symptoms in order to examine whether different levels of DNA methylation were associated with depressive tendencies. Employing methylation-array technology, a total of 363,887 methylation sites across the genomes were investigated and several candidate CpG sites associated with depressive symptoms were identified, especially annotated to genes linked to a G-protein coupled receptor protein signaling pathway. These data provide a strong impetus for validation studies using a larger cohort and support the possibility that G-protein coupled receptor protein signaling pathways are involved in the pathogenesis of depression.

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.

Epigenome-wide association study of DNA methylation in panic disorder.

BACKGROUND: Panic disorder (PD) is considered to be a multifactorial disorder emerging from interactions among multiple genetic and environmental factors. To date, although genetic studies reported several susceptibility genes with PD, few of them were replicated and the pathogenesis of PD remains to be clarified. Epigenetics is considered to play an important role in etiology of complex traits and diseases, and DNA methylation is one of the major forms of epigenetic modifications. In this study, we performed an epigenome-wide association study of PD using DNA methylation arrays so as to investigate the possibility that different levels of DNA methylation might be associated with PD. METHODS: The DNA methylation levels of CpG sites across the genome were examined with genomic DNA samples (PD, N = 48, control, N = 48) extracted from peripheral blood. Methylation arrays were used for the analysis. ß values, which represent the levels of DNA methylation, were normalized via an appropriate pipeline. Then, ß values were converted to M values via the logit transformation for epigenome-wide association study. The relationship between each DNA methylation site and PD was assessed by linear regression analysis with adjustments for the effects of leukocyte subsets. RESULTS: Forty CpG sites showed significant association with PD at 5% FDR correction, though the differences of the DNA methylation levels were relatively small. Most of the significant CpG sites (37/40 CpG sites) were located in or around CpG islands. Many of the significant CpG sites (27/40 CpG sites) were located upstream of genes, and all such CpG sites with the exception of two were hypomethylated in PD subjects. A pathway analysis on the genes annotated to the significant CpG sites identified several pathways, including "positive regulation of lymphocyte activation." CONCLUSIONS: Although future studies with larger number of samples are necessary to confirm the small DNA methylation abnormalities associated with PD, there is a possibility that several CpG sites might be associated, together as a group, with PD.

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.

Phenome-wide association study maps new diseases to the human major histocompatibility complex region.

BACKGROUND: Over 160 disease phenotypes have been mapped to the major histocompatibility complex (MHC) region on chromosome 6 by genome-wide association study (GWAS), suggesting that the MHC region as a whole may be involved in the aetiology of many phenotypes, including unstudied diseases. The phenome-wide association study (PheWAS), a powerful and complementary approach to GWAS, has demonstrated its ability to discover and rediscover genetic associations. The objective of this study is to comprehensively investigate the MHC region by PheWAS to identify new phenotypes mapped to this genetically important region. METHODS: In the current study, we systematically explored the MHC region using PheWAS to associate 2692 MHC-linked variants (minor allele frequency ≥0.01) with 6221 phenotypes in a cohort of 7481 subjects from the Marshfield Clinic Personalized Medicine Research Project. RESULTS: Findings showed that expected associations previously identified by GWAS could be identified by PheWAS (eg, psoriasis, ankylosing spondylitis, type I diabetes and coeliac disease) with some having strong cross-phenotype associations potentially driven by pleiotropic effects. Importantly, novel associations with eight diseases not previously assessed by GWAS (eg, lichen planus) were also identified and replicated in an independent population. Many of these associated diseases appear to be immune-related disorders. Further assessment of these diseases in 16 484 Marshfield Clinic twins suggests that some of these diseases, including lichen planus, may have genetic aetiologies. CONCLUSIONS: These results demonstrate that the PheWAS approach is a powerful and novel method to discover SNP-disease associations, and is ideal when characterising cross-phenotype associations, and further emphasise the importance of the MHC region in human health and disease.