Discrimination between healthy participants and people with panic disorder based on polygenic scores for psychiatric disorders and for intermediate phenotypes using machine learning.

OBJECTIVE: Panic disorder is a modestly heritable condition. Currently, diagnosis is based only on clinical symptoms; identifying objective biomarkers and a more reliable diagnostic procedure is desirable. We investigated whether people with panic disorder can be reliably diagnosed utilizing combinations of multiple polygenic scores for psychiatric disorders and their intermediate phenotypes, compared with single polygenic score approaches, by applying specific machine learning techniques. METHODS: Polygenic scores for 48 psychiatric disorders and intermediate phenotypes based on large-scale genome-wide association studies (n = 7556-1,131,881) were calculated for people with panic disorder (n = 718) and healthy controls (n = 1717). Discrimination between people with panic disorder and healthy controls was based on the 48 polygenic scores using five methods for classification: logistic regression, neural networks, quadratic discriminant analysis, random forests and a support vector machine. Differences in discrimination accuracy (area under the curve) due to an increased number of polygenic score combinations and differences in the accuracy across five classifiers were investigated. RESULTS: All five classifiers performed relatively well for distinguishing people with panic disorder from healthy controls by increasing the number of polygenic scores. Of the 48 polygenic scores, the polygenic score for anxiety UK Biobank was the most useful for discrimination by the classifiers. In combinations of two or three polygenic scores, the polygenic score for anxiety UK Biobank was included as one of polygenic scores in all classifiers. When all 48 polygenic scores were used in combination, the greatest areas under the curve significantly differed among the five classifiers. Support vector machine and logistic regression had higher accuracy than quadratic discriminant analysis and random forests. For each classifier, the greatest area under the curve was 0.600 ± 0.030 for logistic regression (polygenic score combinations N = 14), 0.591 ± 0.039 for neural networks (N = 9), 0.603 ± 0.033 for quadratic discriminant analysis (N = 10), 0.572 ± 0.039 for random forests (N = 25) and 0.617 ± 0.041 for support vector machine (N = 11). The greatest areas under the curve at the best polygenic score combination significantly differed among the five classifiers. Random forests had the lowest accuracy among classifiers. Support vector machine had higher accuracy than neural networks. CONCLUSIONS: These findings suggest that increasing the number of polygenic score combinations up to approximately 10 effectively improved the discrimination accuracy and that support vector machine exhibited greater accuracy among classifiers. However, the discrimination accuracy for panic disorder, when based solely on polygenic score combinations, was found to be modest.

A loss-of-function variant in SUV39H2 identified in autism-spectrum disorder causes altered H3K9 trimethylation and dysregulation of protocadherin ß-cluster genes in the developing brain.

Recent evidence has documented the potential roles of histone-modifying enzymes in autism-spectrum disorder (ASD). Aberrant histone H3 lysine 9 (H3K9) dimethylation resulting from genetic variants in histone methyltransferases is known for neurodevelopmental and behavioral anomalies. However, a systematic examination of H3K9 methylation dynamics in ASD is lacking. Here we resequenced nine genes for histone methyltransferases and demethylases involved in H3K9 methylation in individuals with ASD and healthy controls using targeted next-generation sequencing. We identified a novel rare variant (A211S) in the SUV39H2, which was predicted to be deleterious. The variant showed strongly reduced histone methyltransferase activity in vitro. In silico analysis showed that the variant destabilizes the hydrophobic core and allosterically affects the enzyme activity. The Suv39h2-KO mice displayed hyperactivity and reduced behavioral flexibility in learning the tasks that required complex behavioral adaptation, which is relevant for ASD. The Suv39h2 deficit evoked an elevated expression of a subset of protocadherin ß (Pcdhb) cluster genes in the embryonic brain, which is attributable to the loss of H3K9 trimethylation (me3) at the gene promoters. Reduced H3K9me3 persisted in the cerebellum of Suv39h2-deficient mice to an adult stage. Congruently, reduced expression of SUV39H1 and SUV39H2 in the postmortem brain samples of ASD individuals was observed, underscoring the role of H3K9me3 deficiency in ASD etiology. The present study provides direct evidence for the role of SUV39H2 in ASD and suggests a molecular cascade of SUV39H2 dysfunction leading to H3K9me3 deficiency followed by an untimely, elevated expression of Pcdhb cluster genes during early neurodevelopment.

Does cognitive behavioral therapy for anxiety disorders assist the discontinuation of benzodiazepines among patients with anxiety disorders? A systematic review and meta-analysis.

Long-term use of benzodiazepines (BZD) is not recommended for the treatment of anxiety disorders. Cognitive behavioral therapy (CBT) is an effective treatment option for discontinuation of BZD in patients with anxiety disorders. This systematic review and meta-analysis sought to clarify whether CBT is effective for discontinuing BZD anxiolytics in patients with anxiety disorders. This study was preregistered with PROSPERO (registration number: CRD42019125263). A literature search of major electronic databases was conducted in December 2018. Three randomized controlled trials were included in this review, and meta-analyses were performed. The proportion of discontinuing BZD anxiolytics was significantly higher in the CBT plus gradual tapering group than in the gradual tapering alone group, both in the short term (3 months after allocation; number needed to treat: 3.2, 95% confidence interval [CI]: 2.1 to 7.1; risk ratio: 1.96, 95%CI: 1.29 to 2.98, P = 0.002, three studies) and long term (6 to 12 months after allocation; number needed to treat: 2.8, 95%CI: 1.9 to 5.3; risk ratio: 2.16, 95%CI: 1.41 to 3.32, P = 0.0004, three studies). CBT may be effective for discontinuing BZD anxiolytics, both in the short term and in the long term after the allocation. Further studies with larger sample sizes are necessary to draw definitive conclusions regarding the efficacy and safety of CBT for discontinuing BZD anxiolytics in patients with anxiety disorders.

Shared transethnic genetic basis of panic disorder and psychiatric and related intermediate phenotypes.

Panic disorder (PD), a common anxiety disorder, is modestly heritable. The genetic basis of anxiety disorders overlaps with that of other psychiatric disorders and their intermediate phenotypes in individuals of European ancestry. Here, we investigated the transethnic polygenetic features shared between Japanese PD patients and European patients with psychiatric disorders and their intermediate phenotypes by conducting polygenic risk score (PRS) analyses. Large-scale European genome-wide association study (GWAS) datasets (n = 7,556-1,131,881) for ten psychiatric disorders and seven intermediate phenotypes were utilized as discovery samples. PRSs derived from these GWASs were calculated for Japanese target subjects [718 PD patients and 1,717 healthy controls (HCs)]. The effects of these PRSs from European GWASs on the risk of PD in Japanese patients were investigated. The PRSs from European studies of anxiety disorders were marginally higher in Japanese PD patients than in HCs (p = 0.013). Regarding other psychiatric disorders, the PRSs for depression in European patients were significantly higher in Japanese PD patients than in HCs (p = 2.31×10-4), while the PRSs for attention-deficit/hyperactivity disorder in European patients were nominally lower in Japanese PD patients than in HCs (p = 0.024). Regarding health-related, personality-based and cognitive intermediate phenotypes, the PRSs for loneliness (especially isolation) in European individuals were significantly higher in Japanese PD patients than in HCs (p = 9.02×10-4). Furthermore, Japanese PD patients scored nominally higher than HCs in PRSs for neuroticism in European people (p = 3.37×10-3), while Japanese PD patients scored nominally lower than HCs in PRSs for tiredness (p = 0.025), educational attainment (p = 0.035) and cognitive function (p = 9.63×10-3). Our findings suggest that PD shares transethnic genetic etiologies with other psychiatric disorders and related intermediate phenotypes.

Association of Serotonin Transporter Gene (5-HTTLPR/rs25531) Polymorphism with Comorbidities of Panic Disorder.

INTRODUCTION: Panic disorder (PD) has many comorbidities such as depression, bipolar disorder (BPD), and agoraphobia (AG). PD is a moderately heritable anxiety disorder whose pathogenesis is not well understood. Recently, a tri-allelic serotonin transporter (5-HTTLPR/rs25531) polymorphism was reported to be more sensitive to personality traits compared to the bi-allelic 5-HTTLPR polymorphism. We hypothesized that the 5-HTTLPR/rs25531 polymorphism may lead to a pathological anxious state depending on the presence or absence of a comorbidity in PD. METHODS: In this study, we investigated the relationship between comorbidities in PD and tri-allelic 5-HTTLPR polymorphisms. A total of 515 patients with PD (148 males, 367 females) were genotyped, and the Revised NEO Personality Inventory as well as anxiety-related psychological tests were administered. Depression, BPD, and AG were diagnosed as comorbidities. RESULTS: For the tri-allele 5-HTTLPR genotype, a significant interaction effect was found between openness to experience and comorbid depression. Examination of the interaction between AG and the tri-allelic 5-HTTLPR genotype revealed that L' allele carriers are associated with higher trait anxiety than the S'S' genotype group in PD without AG. CONCLUSION: Some anxiety and personality traits can be characterized by the tri-allelic gene effect of 5-HTTLPR. These results suggest that tri-allelic 5-HTTLPR genotypes have genetic effects on the presence of comorbidities of PD.

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.

Shared genetic etiology between anxiety disorders and psychiatric and related intermediate phenotypes.

BACKGROUND: Psychiatric disorders and related intermediate phenotypes are highly heritable and have a complex, overlapping polygenic architecture. A large-scale genome-wide association study (GWAS) of anxiety disorders identified genetic variants that are significant on a genome-wide. The current study investigated the genetic etiological overlaps between anxiety disorders and frequently cooccurring psychiatric disorders and intermediate phenotypes. METHODS: Using case-control and factor score models, we investigated the genetic correlations of anxiety disorders with eight psychiatric disorders and intermediate phenotypes [the volumes of seven subcortical brain regions, childhood cognition, general cognitive ability and personality traits (subjective well-being, loneliness, neuroticism and extraversion)] from large-scale GWASs (n = 7556-298 420) by linkage disequilibrium score regression. RESULTS: Among psychiatric disorders, the risk of anxiety disorders was positively genetically correlated with the risks of major depressive disorder (MDD) (rg ± standard error = 0.83 ± 0.16, p = 1.97 × 10-7), schizophrenia (SCZ) (0.28 ± 0.09, p = 1.10 × 10-3) and attention-deficit/hyperactivity disorder (ADHD) (0.34 ± 0.13, p = 8.40 × 10-3). Among intermediate phenotypes, significant genetic correlations existed between the risk of anxiety disorders and neuroticism (0.81 ± 0.17, p = 1.30 × 10-6), subjective well-being (-0.73 ± 0.18, p = 4.89 × 10-5), general cognitive ability (-0.23 ± 0.08, p = 4.70 × 10-3) and putamen volume (-0.50 ± 0.18, p = 5.00 × 10-3). No other significant genetic correlations between anxiety disorders and psychiatric or intermediate phenotypes were observed (p > 0.05). The case-control model yielded stronger genetic effect sizes than the factor score model. CONCLUSIONS: Our findings suggest that common genetic variants underlying the risk of anxiety disorders contribute to elevated risks of MDD, SCZ, ADHD and neuroticism and reduced quality of life, putamen volume and cognitive performance. We suggest that the comorbidity of anxiety disorders is partly explained by common genetic variants.

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.

Whole-exome sequencing and gene-based rare variant association tests suggest that PLA2G4E might be a risk gene for panic disorder.

Panic disorder (PD) is characterized by recurrent and unexpected panic attacks, subsequent anticipatory anxiety, and phobic avoidance. Recent epidemiological and genetic studies have revealed that genetic factors contribute to the pathogenesis of PD. We performed whole-exome sequencing on one Japanese family, including multiple patients with panic disorder, which identified seven rare protein-altering variants. We then screened these genes in a Japanese PD case-control group (384 sporadic PD patients and 571 controls), resulting in the detection of three novel single nucleotide variants as potential candidates for PD (chr15: 42631993, T>C in GANC; chr15: 42342861, G>T in PLA2G4E; chr20: 3641457, G>C in GFRA4). Statistical analyses of these three genes showed that PLA2G4E yielded the lowest p value in gene-based rare variant association tests by Efficient and Parallelizable Association Container Toolbox algorithms; however, the p value did not reach the significance threshold in the Japanese. Likewise, in a German case-control study (96 sporadic PD patients and 96 controls), PLA2G4E showed the lowest p value but again did not reach the significance threshold. In conclusion, we failed to find any significant variants or genes responsible for the development of PD. Nonetheless, our results still leave open the possibility that rare protein-altering variants in PLA2G4E contribute to the risk of PD, considering the function of this gene.

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.

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.

Oxytocin receptor gene variations predict neural and behavioral response to oxytocin in autism.

Oxytocin appears beneficial for autism spectrum disorder (ASD), and more than 20 single-nucleotide polymorphisms (SNPs) in oxytocin receptor (OXTR) are relevant to ASD. However, neither biological functions of OXTR SNPs in ASD nor critical OXTR SNPs that determine oxytocin's effects on ASD remains known. Here, using a machine-learning algorithm that was designed to evaluate collective effects of multiple SNPs and automatically identify most informative SNPs, we examined relationships between 27 representative OXTR SNPs and six types of behavioral/neural response to oxytocin in ASD individuals. The oxytocin effects were extracted from our previous placebo-controlled within-participant clinical trial administering single-dose intranasal oxytocin to 38 high-functioning adult Japanese ASD males. Consequently, we identified six different SNP sets that could accurately predict the six different oxytocin efficacies, and confirmed the robustness of these SNP selections against variations of the datasets and analysis parameters. Moreover, major alleles of several prominent OXTR SNPs-including rs53576 and rs2254298-were found to have dissociable effects on the oxytocin efficacies. These findings suggest biological functions of the OXTR SNP variants on autistic oxytocin responses, and implied that clinical oxytocin efficacy may be genetically predicted before its actual administration, which would contribute to establishment of future precision medicines for ASD.

The First Pilot Genome-Wide Gene-Environment Study of Depression in the Japanese Population.

Stressful events have been identified as a risk factor for depression. Although gene-environment (G × E) interaction in a limited number of candidate genes has been explored, no genome-wide search has been reported. The aim of the present study is to identify genes that influence the association of stressful events with depression. Therefore, we performed a genome-wide G × E interaction analysis in the Japanese population. A genome-wide screen with 320 subjects was performed using the Affymetrix Genome-Wide Human Array 6.0. Stressful life events were assessed using the Social Readjustment Rating Scale (SRRS) and depression symptoms were assessed with self-rating questionnaires using the Center for Epidemiologic Studies Depression (CES-D) scale. The p values for interactions between single nucleotide polymorphisms (SNPs) and stressful events were calculated using the linear regression model adjusted for sex and age. After quality control of genotype data, a total of 534,848 SNPs on autosomal chromosomes were further analyzed. Although none surpassed the level of the genome-wide significance, a marginal significant association of interaction between SRRS and rs10510057 with depression were found (p = 4.5 × 10-8). The SNP is located on 10q26 near Regulators of G-protein signaling 10 (RGS10), which encodes a regulatory molecule involved in stress response. When we investigated a similar G × E interaction between depression (K6 scale) and work-related stress in an independent sample (n = 439), a significant G × E effect on depression was observed (p = 0.015). Our findings suggest that rs10510057, interacting with stressors, may be involved in depression risk. Incorporating G × E interaction into GWAS can contribute to find susceptibility locus that are potentially missed by conventional GWAS.

Assessment and characterization of phenotypic heterogeneity of anxiety disorders across five large cohorts.

To achieve sample sizes necessary for effectively conducting genome-wide association studies (GWASs), researchers often combine data from samples possessing multiple potential sources of heterogeneity. This is particularly relevant for psychiatric disorders, where symptom self-report, differing assessment instruments, and diagnostic comorbidity complicates the phenotypes and contribute to difficulties with detecting and replicating genetic association signals. We investigated sources of heterogeneity of anxiety disorders (ADs) across five large cohorts used in a GWAS meta-analysis project using a dimensional structural modeling approach including confirmatory factor analyses (CFAs) and measurement invariance (MI) testing. CFA indicated a single-factor model provided the best fit in each sample with the same pattern of factor loadings. MI testing indicated degrees of failure of metric and scalar invariance which depended on the inclusion of the effects of sex and age in the model. This is the first study to examine the phenotypic structure of psychiatric disorder phenotypes simultaneously across multiple, large cohorts used for GWAS. The analyses provide evidence for higher order invariance but possible break-down at more detailed levels that can be subtly influenced by included covariates, suggesting caution when combining such data. These methods have significance for large-scale collaborative studies that draw on multiple, potentially heterogeneous datasets. Copyright © 2016 John Wiley & Sons, Ltd.

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.

Erratum: Novel rare variations of the oxytocin receptor (OXTR) gene in autism spectrum disorder individuals.

[This corrects the article DOI: 10.1038/hgv.2015.24.][This corrects the article DOI: 10.1038/hgv.2015.24.].

Polymorphisms in the TMEM132D region are associated with panic disorder in HLA-DRB1*13:02-negative individuals of a Japanese population.

We herein report an association between TMEM132D and panic disorder (PD) in a Japanese population, evaluating the effects of HLA-DRB1*13:02, which we previously reported as a susceptibility genetic factor for PD. SNPs in TMEM132D showed significant associations with PD in subjects without HLA-DRB1*13:02 (rs4759997; P=5.02×10(-6), odds ratio=1.50) but not in those with the HLA allele. TMEM132D might have a role in the development of PD in subjects without HLA-DRB1*13:02.

Genome-wide Association Study of Autism Spectrum Disorder in the East Asian Populations.

Autism spectrum disorder is a heterogeneous neurodevelopmental disorder with strong genetic basis. To identify common genetic variations conferring the risk of ASD, we performed a two-stage genome-wide association study using ASD family and healthy control samples obtained from East Asian populations. A total of 166 ASD families (n = 500) and 642 healthy controls from the Japanese population were used as the discovery cohort. Approximately 900,000 single nucleotide polymorphisms (SNPs) were genotyped using Affymetrix Genome-Wide Human SNP array 6.0 chips. In the replication stage, 205 Japanese ASD cases and 184 healthy controls, as well as 418 Chinese Han trios (n = 1,254), were genotyped by TaqMan platform. Case-control analysis, family based association test, and transmission/disequilibrium test (TDT) were then conducted to test the association. In the discovery stage, significant associations were suggested for 14 loci, including 5 known ASD candidate genes: GPC6, JARID2, YTHDC2, CNTN4, and CSMD1. In addition, significant associations were identified for several novel genes with intriguing functions, such as JPH3, PTPRD, CUX1, and RIT2. After a meta-analysis combining the Japanese replication samples, the strongest signal was found at rs16976358 (P = 6.04 × 10(-7)), which is located near the RIT2 gene. In summary, our results provide independent support to known ASD candidate genes and highlight a number of novel genes warranted to be further investigated in a larger sample set in an effort to improve our understanding of the genetic basis of ASD.

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.