Taku Miyagawa

Narcolepsy is strongly associated with the T-cell receptor alpha locus

Narcolepsy with cataplexy, characterized by sleepiness and rapid onset into REM sleep, affects 1 in 2,000 individuals. Narcolepsy was first shown to be tightly associated with HLA-DR2 (ref. 3) and later sublocalized to DQB1*0602 (ref. 4). Following studies in dogs and mice, a 95% loss of hypocretin-producing cells in postmortem hypothalami from narcoleptic individuals was reported. Using genome-wide association (GWA) in Caucasians with replication in three ethnic groups, we found association between narcolepsy and polymorphisms in the TRA@ (T-cell receptor alpha) locus, with highest significance at rs1154155 (average allelic odds ratio 1.69, genotypic odds ratios 1.94 and 2.55, P < 10−21, 1,830 cases, 2,164 controls). This is the first documented genetic involvement of the TRA@ locus, encoding the major receptor for HLA-peptide presentation, in any disease. It is still unclear how specific HLA alleles confer susceptibility to over 100 HLA-associated disorders; thus, narcolepsy will provide new insights on how HLA–TCR interactions contribute to organ-specific autoimmune targeting and may serve as a model for over 100 other HLA-associated disorders.

Common variants in P2RY11 are associated with narcolepsy

Growing evidence supports the hypothesis that narcolepsy with cataplexy is an autoimmune disease. We here report genome-wide association analyses for narcolepsy with replication and fine mapping across three ethnic groups (3,406 individuals of European ancestry, 2,414 Asians and 302 African Americans). We identify a SNP in the 3′ untranslated region of P2RY11, the purinergic receptor subtype P2Y11 gene, which is associated with narcolepsy (rs2305795, combined P = 6.1 × 10−10, odds ratio = 1.28, 95% CI 1.19–1.39, n = 5689). The disease-associated allele is correlated with reduced expression of P2RY11 in CD8+ T lymphocytes (72% reduced, P = 0.003) and natural killer (NK) cells (70% reduced, P = 0.031), but not in other peripheral blood mononuclear cell types. The low expression variant is also associated with reduced P2RY11-mediated resistance to ATP-induced cell death in T lymphocytes (P = 0.0007) and natural killer cells (P = 0.001). These results identify P2RY11 as an important regulator of immune-cell survival, with possible implications in narcolepsy and other autoimmune diseases.

Association of the oxytocin receptor (OXTR) gene polymorphisms with autism spectrum disorder (ASD) in the Japanese population.

The oxytocin receptor (OXTR) gene, which is located on chromosome 3p25.3, has been implicated as a candidate gene for susceptibility of autism spectrum disorder (ASD). Positive associations between OXTR and ASD have been reported in earlier studies. However, the results were inconsistent and demand further studies. In this study, we investigated the associations between OXTR and ASD in a Japanese population by analyzing 11 single-nucleotide polymorphisms (SNPs) using both family-based association test (FBAT) and population-based case–control test. No significant signal was detected in the FBAT test. However, significant differences were observed in allelic frequencies of four SNPs, including rs2254298 between patients and controls. The risk allele of rs2254298 was ‘A’, which was consistent with the previous study in Chinese, and not with the observations in Caucasian. The difference in the risk allele of this SNP in previous studies might be attributable to an ethnic difference in the linkage disequilibrium structure between the Asians and Caucasians. In addition, haplotype analysis exhibits a significant association between a five-SNP haplotype and ASD, including rs22542898. In conclusion, our study might support that OXTR has a significant role in conferring the risk of ASD in the Japanese population.

Variant between CPT1B and CHKB associated with susceptibility to narcolepsy

Narcolepsy (hypocretin deficiency), a sleep disorder characterized by sleepiness, cataplexy and rapid eye movement (REM) sleep abnormalities, is tightly associated with HLA-DRB1*1501 (M17378) and HLA-DQB1*0602 (M20432). Susceptibility genes other than those in the HLA region are also likely involved. We conducted a genome-wide association study using 500K SNP microarrays in 222 Japanese individuals with narcolepsy and 389 Japanese controls, with replication of top hits in 159 Japanese individuals with narcolepsy and 190 Japanese controls, followed by the testing of 424 Koreans, 785 individuals of European descent and 184 African Americans. rs5770917, a SNP located between CPT1B and CHKB, was associated with narcolepsy in Japanese (rs5770917[C], odds ratio (OR) = 1.79, combined P = 4.4 × 10−7) and other ancestry groups (OR = 1.40, P = 0.02). Real-time quantitative PCR assays in white blood cells indicated decreased CPT1B and CHKB expression in subjects with the C allele, suggesting that a genetic variant regulating CPT1B or CHKB expression is associated with narcolepsy. Either of these genes is a plausible candidate, as CPT1B regulates β-oxidation, a pathway involved in regulating theta frequency during REM sleep, and CHKB is an enzyme involved in the metabolism of choline, a precursor of the REM- and wake-regulating neurotransmitter acetylcholine.

Genome-wide association study of panic disorder in the Japanese population.

Panic disorder (PD) is an anxiety disorder characterized by panic attacks and anticipatory anxiety. Although a number of association studies have been conducted, no gene has been identified as a susceptibility locus. In this study, we conducted a genome-wide association study of PD in 200 Japanese patients and the same number of controls, using the GeneChip Human Mapping 500 K Array Set. Genotypes were determined using the Bayesian Robust Linear Model with Mahalanobis (BRLMM) genotype calling algorithm. The genotype data were data-cleaned using criteria for SNP call rate (⩾95%), Hardy–Weinberg equilibrium (P⩾0.1%) and minor allele frequency (⩾5%). The significance level of the allele P-value was set at 1.0 × 10−6, to make false discovery rate (FDR) <0.05. As a result, seven SNPs were significantly associated with PD, which were located in or adjacent to genes including PKP1, PLEKHG1, TMEM16B, CALCOCO1, SDK2 and CLU (or APO-J). Studies with other samples are required to confirm the results.

Meta-analysis of genome-wide association studies for panic disorder in the Japanese population

Panic disorder (PD) is a moderately heritable anxiety disorder whose pathogenesis is not well understood. Due to the lack of power in previous association studies, genes that are truly associated with PD might not be detected. In this study, we conducted a genome-wide association study (GWAS) in two independent data sets using the Affymetrix Mapping 500K Array or Genome-Wide Human SNP Array 6.0. We obtained imputed genotypes for each GWAS and performed a meta-analysis of two GWAS data sets (718 cases and 1717 controls). For follow-up, 12 single-nucleotide polymorphisms (SNPs) were tested in 329 cases and 861 controls. Gene ontology enrichment and candidate gene analyses were conducted using the GWAS or meta-analysis results. We also applied the polygenic score analysis to our two GWAS samples to test the hypothesis of polygenic components contributing to PD. Although genome-wide significant SNPs were not detected in either of the GWAS nor the meta-analysis, suggestive associations were observed in several loci such as BDKRB2 (P=1.3 × 10−5, odds ratio=1.31). Among previous candidate genes, supportive evidence for association of NPY5R with PD was obtained (gene-wise corrected P=6.4 × 10−4). Polygenic scores calculated from weakly associated SNPs (P<0.3 and 0.4) in the discovery sample were significantly associated with PD status in the target sample in both directions (sample I to sample II and vice versa) (P<0.05). Our findings suggest that large sets of common variants of small effects collectively account for risk of PD.

Replication of a genome-wide association study of panic disorder in a Japanese population

Panic disorder (PD) is an anxiety disorder characterized by recurrent and unexpected panic attacks, subsequent worry and phobic avoidance. Although a number of association and linkage studies have been conducted, no gene has been identified as a susceptibility locus. We previously conducted a genome-wide association analysis of PD in 200 Japanese patients and the same number of controls, using a 500 K single nucleotide polymorphisms (SNPs) chip. In this study, we report a replication analysis of PD using the DigTag2 assay. The second stage sample consisted of 558 Japanese patients and 566 controls. Thirty-two markers were tested in a replication sample. As a result, no significant association was found after correction for multiple testing. However, the difference was observed at the nominal allele P-value <0.05 for two SNPs (rs6733840 and rs132617). We also conducted haplotype analyses of SNPs in the APOL3 and CLU genes. Our results failed to show any significant association with PD in these genes. Further studies on these variants with a larger sample size may be worth testing to confirm the results.

Appropriate data cleaning methods for genome-wide association study

AbstractGenome-wide association studies (GWAS) using a large number of single nucleotide polymorphisms (SNPs) have successfully been applied to identify genetic variants of common diseases. However, genotyping using the new array technologies is often associated with spurious results that could unfavorably affect analyses of GWAS. Consequently, data cleaning is of paramount importance in excluding spurious genotyping results. In this study, we investigated the criteria required for the appropriate cleaning of 389 unrelated healthy Japanese samples analyzed using the GeneChip Human Mapping 500K Array Set for GWAS. The samples were randomly subdivided into two groups, and the allele frequencies in the groups were compared for individual SNPs as a quasi-case-control study. Then, observed results were filtered by four parameters (SNP call rate, confidence score obtained using the Bayesian Robust Linear Model with Mahalanobis genotype-calling algorithm, Hardy-Weinberg equilibrium, and minor allele frequency) and assessed for deviation from the null hypothesis. We found that appropriate data cleaning could be achieved using these four parameters. Our findings offer an avenue for obtaining appropriate data from GWAS.

Polymorphisms of promoter and coding regions of the arylamine N-acetyltransferase 2 (NAT2) gene in the Indonesian population: proposal for a new nomenclature

AbstractPolymorphisms of arylamine N-acetyltransferase 2 (NAT2) are reportedly associated with the risk of drug toxicities and development of various diseases. The present study examined NAT2 polymorphisms in both promoter and coding regions in the Indonesian population using PCR direct sequencing. The promoter and coding regions of NAT2 displayed 23 polymorphisms/variations, including eight new ones. Seven haplotypes in the promoter region and six haplotypes in the coding region were inferred. The haplotypes in promoter and coding regions showed limited combinations, and 13 combined haplotypes were inferred. The most frequent haplotypes were U1 (38.9%), U2 (33.5%) in the promoter region and NAT2*4 (37.3%), NAT2*6A (36.8%) in the coding region. When converted to predicted phenotypes, the studied population comprised 65.4% rapid acetylators and 35.6% slow acetylators according to bimodal distribution. According to trimodal distribution, frequencies of predicted phenotypes were 13.6, 50.8 and 35.6% for rapid, intermediate and slow acetylators, respectively. Frequencies of NAT2 alleles for the Indonesian population resembled those of other Southeast Asian populations. We also propose a new NAT2 nomenclature composed of haplotypes in the promoter region and conventional NAT2 haplotypes in the coding region, symbolized by NAT2*4.U1, NAT2*4.U2, NAT2*4.U3, NAT2*4.U5, NAT2*4.U6, NAT2*4.U7, NAT2*6A.U1, NAT2*7B.U2, NAT2*7B.U3, NAT2*5B.U1, NAT2*5B.U4, NAT2*12A.U4 and NAT2*13.U1.

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.