Systematic analysis of exonic germline and postzygotic de novo mutations in bipolar disorder.

Bipolar disorder is a severe mental illness characterized by recurrent manic and depressive episodes. To better understand its genetic architecture, we analyze ultra-rare de novo mutations in 354 trios with bipolar disorder. For germline de novo mutations, we find significant enrichment of loss-of-function mutations in constrained genes (corrected-P = 0.0410) and deleterious mutations in presynaptic active zone genes (FDR = 0.0415). An analysis integrating single-cell RNA-sequencing data identifies a subset of excitatory neurons preferentially expressing the genes hit by deleterious mutations, which are also characterized by high expression of developmental disorder genes. In the analysis of postzygotic mutations, we observe significant enrichment of deleterious ones in developmental disorder genes (P = 0.00135), including the SRCAP gene mutated in two unrelated probands. These data collectively indicate the contributions of both germline and postzygotic mutations to the risk of bipolar disorder, supporting the hypothesis that postzygotic mutations of developmental disorder genes may contribute to bipolar disorder.

Determination of novel CYP2D6 haplotype using the targeted sequencing followed by the long-read sequencing and the functional characterization in the Japanese population.

Next-generation sequencing (NGS) has identified variations in cytochrome P450 (CYP) 2D6 associated with drug responses. However, determination of novel haplotypes is difficult because of the short reads generated by NGS. We aimed to identify novel CYP2D6 variants in the Japanese population and predict the CYP2D6 phenotype based on in vitro metabolic studies. Using a targeted NGS panel (PKSeq), 990 Japanese genomes were sequenced, and then novel CYP2D6 haplotypes were determined. Km, Vmax, and intrinsic clearance (Vmax/Km) of N-desmethyl-tamoxifen 4-hydroxylation were calculated by in vitro metabolic studies using cDNA-expressed CYP2D6 proteins. After determination of the CYP2D6 diplotypes, phenotypes of the individuals were predicted based on the in vitro metabolic activities. Targeted NGS identified 14 CYP2D6 variants not registered in the Pharmacogene Variation Consortium (PharmVar) database. Ten novel haplotypes were registered as CYP2D6*128 to *137 alleles in the PharmVar database. Based on the Vmax/Km value of each allele, *128, *129, *130, *131, *132, and *133 were predicted to be nonfunctional alleles. According to the results of the present study, six normal metabolizers (NM) and one intermediate (IM) metabolizers were designated as IM and poor metabolizers (PM), respectively. Our findings provide important insights into novel haplotypes and haplotypes of CYP2D6 and the effects on in vitro metabolic activities.

Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases.

The overwhelming majority of participants in current genetic studies are of European ancestry. To elucidate disease biology in the East Asian population, we conducted a genome-wide association study (GWAS) with 212,453 Japanese individuals across 42 diseases. We detected 320 independent signals in 276 loci for 27 diseases, with 25 novel loci (P < 9.58 × 10-9). East Asian-specific missense variants were identified as candidate causal variants for three novel loci, and we successfully replicated two of them by analyzing independent Japanese cohorts; p.R220W of ATG16L2 (associated with coronary artery disease) and p.V326A of POT1 (associated with lung cancer). We further investigated enrichment of heritability within 2,868 annotations of genome-wide transcription factor occupancy, and identified 378 significant enrichments across nine diseases (false discovery rate < 0.05) (for example, NKX3-1 for prostate cancer). This large-scale GWAS in a Japanese population provides insights into the etiology of complex diseases and highlights the importance of performing GWAS in non-European populations.

Global DNA hypomethylation and its correlation to the betaine level in peripheral blood of patients with schizophrenia.

Accumulating evidence suggests that aberrant epigenetic regulation is involved in the pathophysiology of major psychiatric disorders such as schizophrenia (SZ) and bipolar disorder (BD). We previously showed that the plasma level of betaine (N,N,N-trimethylglycine), a methyl-group donor, was significantly decreased in patients with first episode schizophrenia (FESZ). In this study, we identified decrease of global DNA methylation level in FESZ (N = 24 patients vs N = 42 controls), and found that global DNA methylation level was inversely correlated with scores on the global assessment of functioning (GAF) scale, and positively correlated with plasma betaine level. Notably, correlations between levels of betaine and its metabolites (N,N-dimethylglycine and sarcosine, N-methylglycine) were lower or lost in FESZ plasma, but remained high in controls. We further examined global DNA methylation levels in patients with chronic SZ (N = 388) and BD (N = 414) as well as controls (N = 430), and confirmed significant hypomethylation and decreased betaine level in SZ. We also found that patients with BD type I, but not those with BD type II, showed significant global hypomethylation. These results suggest that global hypomethylation associated with decreased betaine level in blood cells is common to SZ and BD, and may reflect common pathophysiology such as psychotic symptoms.

GWAS of smoking behaviour in 165,436 Japanese people reveals seven new loci and shared genetic architecture.

Cigarette smoking is a risk factor for a wide range of human diseases1. To investigate the genetic components associated with smoking behaviours in the Japanese population, we conducted a genome-wide association study of four smoking-related traits using up to 165,436 individuals. In total, we identified seven new loci, including three loci associated with the number of cigarettes per day (EPHX2-CLU, RET and CUX2-ALDH2), three loci associated with smoking initiation (DLC1, CXCL12-TMEM72-AS1 and GALR1-SALL3) and LINC01793-MIR4432HG, associated with the age of smoking initiation. Of these, three loci (LINC01793-MIR4432HG, CXCL12-TMEM72-AS1 and GALR1-SALL3) were found by conducting an additional sex-stratified genome-wide association study. This additional analysis showed heterogeneity of effects between sexes. The cross-sex linkage disequilibrium score regression2,3 analysis also indicated that the genetic component of smoking initiation was significantly different between the sexes. Cross-trait linkage disequilibrium score regression analysis and trait-relevant tissue analysis showed that the number of cigarettes per day has a specific genetic background distinct from those of the other three smoking behaviours. We also report 11 diseases that share genetic basis with smoking behaviours. Although the current study should be carefully considered owing to the lack of replication samples, our findings characterized the genetic architecture of smoking behaviours. Further studies in East Asian populations are warranted to confirm our findings.

Decreased serum pyridoxal levels in schizophrenia: meta-analysis and Mendelian randomization analysis.

BACKGROUND: Alterations in one-carbon metabolism have been associated with schizophrenia, and vitamin B6 is one of the key components in this pathway. METHODS: We first conducted a case-control study of serum pyridoxal levels and schizophrenia in a large Japanese cohort (n = 1276). Subsequently, we conducted a meta-analysis of association studies (n = 2125). Second, we investigated whether rs4654748, which was identified in a genome-wide association study as a vitamin B6-related single nucleotide polymorphism, was genetically implicated in patients with schizophrenia in the Japanese population (n = 10 689). Finally, we assessed the effect of serum pyridoxal levels on schizophrenia risk using a Mendelian randomization (MR) approach. RESULTS: Serum pyridoxal levels were significantly lower in patients with schizophrenia than in controls, not only in our cohort, but also in the pooled data set of the meta-analysis of association studies (standardized mean difference -0.48, 95% confidence interval [CI] -0.57 to -0.39, p = 9.8 × 10-24). We failed to find a significant association between rs4654748 and schizophrenia. Furthermore, an MR analysis failed to find a causal relationship between pyridoxal levels and schizophrenia risk (odds ratio 0.99, 95% CI 0.65-1.51, p = 0.96). LIMITATIONS: Food consumption and medications may have affected serum pyridoxal levels in our cross-sectional study. Sample size, number of instrumental variables and substantial heterogeneity among patients with schizophrenia are limitations of an MR analysis. CONCLUSION: We found decreased serum pyridoxal levels in patients with schizophrenia in this observational study. However, we failed to obtain data supporting a causal relationship between pyridoxal levels and schizophrenia risk using the MR approach.

Characterization of Gene Expression Patterns among Artificially Developed Cancer Stem Cells Using Spherical Self-Organizing Map.

We performed gene expression microarray analysis coupled with spherical self-organizing map (sSOM) for artificially developed cancer stem cells (CSCs). The CSCs were developed from human induced pluripotent stem cells (hiPSCs) with the conditioned media of cancer cell lines, whereas the CSCs were induced from primary cell culture of human cancer tissues with defined factors (OCT3/4, SOX2, and KLF4). These cells commonly expressed human embryonic stem cell (hESC)/hiPSC-specific genes (POU5F1, SOX2, NANOG, LIN28, and SALL4) at a level equivalent to those of control hiPSC 201B7. The sSOM with unsupervised method demonstrated that the CSCs could be divided into three groups based on their culture conditions and original cancer tissues. Furthermore, with supervised method, sSOM nominated TMED9, RNASE1, NGFR, ST3GAL1, TNS4, BTG2, SLC16A3, CD177, CES1, GDF15, STMN2, FAM20A, NPPB, CD99, MYL7, PRSS23, AHNAK, and LOC152573 genes commonly upregulating among the CSCs compared to hiPSC, suggesting the gene signature of the CSCs.

A Unique Procedure to Identify Cell Surface Markers Through a Spherical Self-Organizing Map Applied to DNA Microarray Analysis.

To identify cell-specific markers, we designed a DNA microarray platform with oligonucleotide probes for human membrane-anchored proteins. Human glioma cell lines were analyzed using microarray and compared with normal and fetal brain tissues. For the microarray analysis, we employed a spherical self-organizing map, which is a clustering method suitable for the conversion of multidimensional data into two-dimensional data and displays the relationship on a spherical surface. Based on the gene expression profile, the cell surface characteristics were successfully mirrored onto the spherical surface, thereby distinguishing normal brain tissue from the disease model based on the strength of gene expression. The clustered glioma-specific genes were further analyzed by polymerase chain reaction procedure and immunocytochemical staining of glioma cells. Our platform and the following procedure were successfully demonstrated to categorize the genes coding for cell surface proteins that are specific to glioma cells. Our assessment demonstrates that a spherical self-organizing map is a valuable tool for distinguishing cell surface markers and can be employed in marker discovery studies for the treatment of cancer.

Individual risk alleles of susceptibility to schizophrenia are associated with poor clinical and social outcomes.

Many patients with schizophrenia have poor clinical and social outcomes. Some risk alleles closely related to the onset of schizophrenia have been reported to be associated with their clinical phenotypes, but the direct relationship between genetic vulnerability to schizophrenia and clinical/social outcomes of schizophrenia, as evaluated by both practical clinical scales and 'real-world' function, has not been investigated. We evaluated the clinical and social outcomes of 455 Japanese patients with schizophrenia by severity of illness according to the Clinical Global Impression-Severity Scale (CGI-S) and social outcomes by social adjustment/maladjustment at 5 years after the first visit. We examined whether 46 single nucleotide polymorphisms (SNPs) selected from a Japanese genome-wide association study of susceptibility to schizophrenia were associated with clinical and social outcomes. We also investigated the polygenic risk scores of 46 SNPs. Allele-wise association analysis detected three SNPs, including rs2623659 in the CUB and Sushi multiple domains-1 (CSMD1) gene, associated with severity of illness at end point. The severity of illness at end point was associated with treatment response, but not with the severity of illness at baseline. Three SNPs, including rs2294424 in the C6orf105 gene, were associated with social outcomes. Point estimates of odds ratios showed positive relationships between polygenic risk scores and clinical/social outcomes; however, the results were not statistically significant. Because these results are exploratory, we need to replicate them with a larger sample in a future study.

Association study of BCL9 gene polymorphism rs583583 with schizophrenia and negative symptoms in Japanese population.

B-cell CLL/lymphoma 9 (BCL9) is located within the schizophrenia (SCZ) suspected locus chr1q21.1. A recent study reported that a single nucleotide polyphormism (SNP) within BCL9 (rs583583) is associated with negative symptoms of Schizophrenia, as measured by the Positive and Negative Syndrome Scale (PANSS), in the Caucasian population. We therefore investigated genetic association of rs583583, and its effect on negative symptoms in the Japanese patients. For association analysis, we used a Japanese sample set comprising 1089 SCZ and 950 controls (CON). Analysis of the effect of rs586586 on negative symptoms as examined by PANSS was investigated using 280 SCZ. Furthermore, for analysis of cognitive performance, we investigated 90 SCZ and 51 CON using the Continuous Performance Test (CPT-IP) and the Wisconsin Card Sorting Test (WCST) Keio version. We did not detect association between rs583583 and SCZ. Furthermore, rs583583 was not associated with PANSS negative scores or with CPT-IT or WCST cognitive tests. Considering the results of our previous study, combined with the results of the current study of rs583583, we argue that BCL9 most likely does not harbor a common genetic variant that can increase the risk for SCZ in the Japanese population.

Analysis of the VAV3 as candidate gene for schizophrenia: evidences from voxel-based morphometry and mutation screening.

In recently completed Japanese genome-wide association studies (GWAS) of schizophrenia (JPN_GWAS) one of the top association signals was detected in the region of VAV3, a gene that maps to the chromosome 1p13.3. In order to complement JPN_GWAS findings, we tested the association of rs1410403 with brain structure in healthy individuals and schizophrenic patients and performed exon resequencing of VAV3. We performed voxel-based morphometry (VBM) and mutation screening of VAV3. Four independent samples were used in the present study: (1) for VBM analysis, we used case-control sample comprising 100 patients with schizophrenia and 264 healthy controls, (2) mutation analysis was performed on a total of 321 patients suffering from schizophrenia, and 2 case-control samples (3) 729 unrelated patients with schizophrenia and 564 healthy comparison subjects, and (4) sample comprising 1511 cases and 1517 healthy comparison subjects and were used for genetic association analysis of novel coding variants with schizophrenia. The VBM analysis suggests that rs1410403 might affect the volume of the left superior and middle temporal gyri (P=.011 and P=.013, respectively), which were reduced in patients with schizophrenia compared with healthy subjects. Moreover, 4 rare novel missense variants were detected. The mutations were followed-up in large independent sample, and one of the novel variants (Glu741Gly) was associated with schizophrenia (P=.02). These findings demonstrate that VAV3 can be seen as novel candidate gene for schizophrenia in which both rare and common variants may be related to increased genetic risk for schizophrenia in Japanese population.

Effects of cyclophosphamide on the kaolin consumption (pica behavior) in five strains of adult male rats.

It is known that pica, the consumption of non-nutritive substances such as kaolin, can be induced by administration of toxins or emetic agents in rats. In the present study, we examined the effects of intraperitoneal (i.p.) administration of cyclophosphamide on pica behavior and on the concentration of 5-hydroxyindoleacetic acids (5HIAA) in cerebrospinal fluid (CSF) in the following five strains of adult male rats: Sprague Dawley (SD), Wistar, Fischer 344 (F344), Wistar-Imamichi (WI) and Long Evans (LE). Cyclophosphamide (25 mg or 50 mg/kg) was injected (i.p.) into the rats and kaolin and food intake were measured at 24 hr after injection. The animals were anesthetized with urethane (1 g/kg) at 3 hr after injection of cyclophosphamide, and CSF was collected from the cisterna magna. WI and LE rats clearly showed pica behavior as compared with the other strains. In LE rats, the concentration of 5HIAA in CSF also increased in a dose-dependent manner of cyclophosphamide. The pretreatment with ondansetron (5-HT(3) antagonist) restored both changes (kaolin consumption and 5HIAA levels) induced by cyclophosphamide. These results suggest that the LE rat is sensitive to cyclophosphamide, that pica induced by cyclophosphamide mimics many aspects of emesis including the serotonergic response in the central nervous system and that use of the pica model would be a practical method for evaluating the effects of antiemetic drugs in addition to the mechanism of emesis.

Serotonin 1A receptor gene is associated with Japanese methamphetamine-induced psychosis patients.

BACKGROUND: Several investigations have reported associations the serotonin 1A (5-HT1A) receptor to schizophrenia and psychotic disorders, making 5-HT1A receptor gene (HTR1A) an adequate candidate gene for the pathophysiology of schizophrenia and methamphetamine (METH)-induced psychosis. Huang and colleagues reported that rs6295 in HTR1A was associated with schizophrenia. The symptoms of methamphetamine (METH)-induced psychosis are similar to those of paranoid type schizophrenia. It may indicate that METH-induced psychosis and schizophrenia have common susceptibility genes. In support of this hypothesis, we reported that the V-act murine thymoma viral oncogene homologue 1 (AKT1) gene was associated with METH-induced psychosis and schizophrenia in the Japanese population. Furthermore, we conducted an analysis of the association of HTR1A with METH-induced psychosis. METHOD: Using one functional SNP (rs6295) and one tagging SNP (rs878567), we conducted a genetic association analysis of case-control samples (197 METH-induced psychosis patients and 337 controls) in the Japanese population. The age and sex of the control subjects did not differ from those of the methamphetamine dependence patients. RESULTS: Rs878567 was associated with METH-induced psychosis patients in the allele/genotype-wise analysis. Moreover, this significance remained after Bonferroni correction. In addition, we detected an association between rs6295 and rs878567 in HTR1A and METH-induced psychosis patients in the haplotype-wise analysis. Although we detected an association between rs6295 and METH-induced psychosis patients, this significance disappeared after Bonferroni correction. CONCLUSION: HTR1A may play an important role in the pathophysiology of METH-induced psychosis in the Japanese population. However, because we did not perform a mutation scan of HTR1A, a replication study using a larger sample may be required for conclusive results.

The different effects of i.c.v. injection of pituitary adenylate cyclase activating polypeptide (PACAP) on prolactin secretion in adult male and lactating rats.

The effects of intracerebroventricular (i.c.v.) administration of pituitary adenylate cyclase activating polypeptide38 (PACAP38) on prolactin (PRL) secretion and the activity of tyrosine hydroxylase (TH) were examined in adult male and lactating rats with or without suckling stimulus. In adult male rats and lactating rats with suckling stimulus, administration of PACAP38 (0.25 or 1 nmol) decreased PRL secretion and increased the activity of TH in the stalk-median eminence. On the other hand, the injection of PACAP38 did not affect PRL secretion and TH activity in lactating rats without sucking stimulus. Administration of PACAP6-38 (4 nmol), a specific receptor antagonist, also had no effect on PRL secretion and TH activity in adult male rats. These results suggest that i.c.v. administration of PACAP inhibits PRL secretion mediated by dopamine neuron within the hypothalamus, but the effects of PACAP differ depending on the physiological condition of animals. These observed effects of PACAP on PRL release may be pharmacological responses rather than physiological responses.

A functional polymorphism in estrogen receptor alpha gene is associated with Japanese methamphetamine induced psychosis.

BACKGROUND: A recent study reported an association between rs2234693, which influences enhancer activity levels in estrogen receptor alpha gene (ESR1), and schizophrenia. This study reported that schizophrenic patients with the CC genotype have significantly lower ESR1 mRNA levels in the prefrontal cortex than patients with other genotypes. The symptoms of methamphetamine induced psychosis are similar to those of paranoid type schizophrenia. Therefore, we conducted an association analysis of rs2234693 with Japanese methamphetamine induced psychosis patients. METHOD: Using rs2234693, we conducted a genetic association analysis of case-control samples (197 methamphetamine induced psychosis patients and 197 healthy controls). The age and sex of the control subjects did not differ from those of the methamphetamine induced psychosis patients. RESULTS: We detected a significant association between ESR1 and methamphetamine induced psychosis patients in allele/genotype-wise analysis. For further interpretation of these associations, we performed single marker analysis of subjects divided by sex. Rs2234693 was associated with male methamphetamine induced psychosis. DISCUSSION: Our results suggest that rs2234693 in ESR1 may play a role in the pathophysiology of Japanese methamphetamine induced psychosis patients.

Glutamate cysteine ligase modifier (GCLM) subunit gene is not associated with methamphetamine-use disorder or schizophrenia in the Japanese population.

A recent study showed a significant association between schizophrenia in European samples and the glutamate cysteine ligase modifier (GCLM) subunit gene, which is the key glutathione (GSH)-synthesizing enzyme. Since the symptoms of methamphetamine (METH)-induced psychosis are similar to those of schizophrenia, the GCLM gene is thought to be a good candidate gene for METH-use disorder or related disorders. To evaluate the association between the GCLM gene and METH-use disorder and schizophrenia, we conducted a case-control study of Japanese subjects (METH-use disorder, 185 cases; schizophrenia, 742 cases; and controls, 819). Four SNPs (2 SNPs from an original report and JSNP database, and 2 "tagging SNPs" from HapMap database) in the GCLM gene were examined in this association analysis; one SNP showed an association with both METH-use disorder and METH-induced psychosis. After Bonferroni's correction for multiple testing, however, this significance disappeared. No significant association was found with schizophrenia. Our findings suggest that a common genetic variation in the GCLM gene might not contribute to the risk of METH-use disorder and schizophrenia in the Japanese population.

Association analysis of AKT1 and schizophrenia in a UK case control sample.

AKT1 (V-akt murine thyoma viral oncogene homolog 1) is involved in intracellular signalling pathways postulated as of aetiological importance in schizophrenia. Markers in the AKT1 gene have also recently been associated with schizophrenia in two samples of European origin and in Japanese and Iranian samples. Aiming to replicate these findings, we examined ten SNPs spanning AKT1 in a UK case-control sample (schizophrenia cases n=673, controls n=716). These included all SNPs previously reported to be associated in European, Japanese and Iranian samples, alone or in haplotypes, as well as additional markers defined by the Haploview Tagger program (pair-wise tagging, minimum r(2)=0.8, minor allele frequency=0.02). We found no association with single markers (min p=0.17). We found weak evidence for association (p=0.04) with a four marker haplotype reported as significant in the original positive European sample of Emamian et al. [Emamian, E.S., Hall, D., Birnbaum, M.J., Karayiorgou, M., Gogos, J.A., 2004. Convergent evidence for impaired AKT1-GSK3beta signaling in schizophrenia. Nat. Genet. 36, 131-137] and also an overlapping three marker haplotype (p=0.016) that had previously been reported as significant in a Japanese sample. Nominal p-values for these haplotypes did not survive correction for multiple testing. Our study provides at best weak support for the hypothesis that AKT1 is a susceptibility gene for schizophrenia. Examination of our own data and those of other groups leads us to conclude that overall, the evidence for association of AKT1 as a susceptibility gene for schizophrenia is weakly positive, but not yet convincing.

Positive association of AKT1 haplotype to Japanese methamphetamine use disorder.

Recent evidence suggests that the AKT1-GSK3beta signalling cascade partially mediates dopamine-dependent behaviours. In relation to the pathophysiology of schizophrenia or methamphetamine (Meth) use disorder, AKT1 is a good candidate gene for such conditions. For schizophrenia, positive associations of SNPs and AKT1 haplotypes were reported in US and Japanese samples. To evaluate the association between AKT1 and Meth-use disorder, we conducted a case-control study of Japanese samples (182 patients and 437 controls). A positive association between a SNP and haplotypes was found, and the 'signal' SNP was the same SNP found to be associated with US schizophrenia, but not with Japanese schizophrenia. Our results indicate that AKT1 may play a possible role in the development of Meth-use disorder. Further investigation of these associations, together with evidence from previous animal studies, may open the way to elucidation of the pathophysiology of this condition.

No association was found between a functional SNP in ZDHHC8 and schizophrenia in a Japanese case-control population.

ZDHHC8 is a new and attractive candidate for a schizophrenia-susceptibility factor. First, several lines of linkage studies showed that 22q11, on which ZDHHC8 is located, is a "hot" region. Second, fine linkage disequilibrium mapping revealed a significant association around ZDHHC8. Moreover, a very recent study reported that one single nucleotide polymorphism (SNP: rs175174) in ZDHHC8 might affect the splicing process, the ZDHHC8 knock-out mice showed the gender-specific phenotype, and the transmission disequilibrium test (TDT) using this SNP also showed significant association with human female schizophrenia. Thus, we attempted a replication study of this SNP using relatively large Japanese case-control samples (561 schizophrenics and 529 controls). No association was found between schizophrenia and controls even after dividing samples by gender. Because our sample size provided quite high power, ZDHHC8 may not play a major role in Japanese schizophrenia. And our results did not support the gender-specific effect of this SNP.