GWAS-identified bipolar disorder risk allele in the FADS1/2 gene region links mood episodes and unsaturated fatty acid metabolism in mutant mice.

Large-scale genome-wide association studies (GWASs) on bipolar disorder (BD) have implicated the involvement of the fatty acid desaturase (FADS) locus. These enzymes (FADS1 and FADS2) are involved in the metabolism of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), which are thought to potentially benefit patients with mood disorders. To model reductions in the activity of FADS1/2 affected by the susceptibility alleles, we generated mutant mice heterozygously lacking both Fads1/2 genes. We measured wheel-running activity over six months and observed bipolar swings in activity, including hyperactivity and hypoactivity. The hyperactivity episodes, in which activity was far above the norm, usually lasted half a day; mice manifested significantly shorter immobility times on the behavioral despair test performed during these episodes. The hypoactivity episodes, which lasted for several weeks, were accompanied by abnormal circadian rhythms and a marked decrease in wheel running, a spontaneous behavior associated with motivation and reward systems. We comprehensively examined lipid composition in the brain and found that levels of certain lipids were significantly altered between wild-type and the heterozygous mutant mice, but no changes were consistent with both sexes and either DHA or EPA was not altered. However, supplementation with DHA or a mixture of DHA and EPA prevented these episodic behavioral changes. Here we propose that heterozygous Fads1/2 knockout mice are a model of BD with robust constitutive, face, and predictive validity, as administration of the mood stabilizer lithium was also effective. This GWAS-based model helps to clarify how lipids and their metabolisms are involved in the pathogenesis and treatment of BD.

Copy number variations in RNF216 and postsynaptic membrane-associated genes are associated with bipolar disorder: a case-control study in the Japanese population.

AIM: Bipolar disorder (BD) is a common psychiatric disorder characterized by alterations between manic/hypomanic and depressive states. Rare pathogenic copy number variations (CNVs) that overlap with exons of synaptic genes have been associated with BD. However, no study has comprehensively explored CNVs in synaptic genes associated with BD. Here, we evaluated the relationship between BD and rare CNVs that overlap with synaptic genes, not limited to exons, in the Japanese population. METHODS: Using array comparative genome hybridization, we detected CNVs in 1839 patients with BD and 2760 controls. We used the Synaptic Gene Ontology database to identify rare CNVs that overlap with synaptic genes. Using gene-based analysis, we compared their frequencies between the BD and control groups. We also searched for synaptic gene sets related to BD. The significance level was set to a false discovery rate of 10%. RESULTS: The RNF216 gene was significantly associated with BD (odds ratio, 4.51 [95% confidence interval, 1.66-14.89], false discovery rate < 10%). The BD-associated CNV that corresponded with RNF216 also partially overlapped with the minimal critical region of the 7p22.1 microduplication syndrome. The integral component of the postsynaptic membrane (Gene Ontology:0099055) was significantly associated with BD. The CNV overlapping with the intron region of GRM5 in this gene set showed a nominal significant association between cases and controls (P < 0.05). CONCLUSION: We provide evidence that CNVs in RNF216 and postsynaptic membrane-related genes confer a risk of BD, contributing to a better understanding of the pathogenesis of BD.

Genetic mapping of simply inherited categorical traits, including anthocyanin accumulation profiles and fruit appearance, in eggplant (Solanum melongena).

BACKGROUND: The appearance quality of the eggplant (Solanum melongena L.) fruit is an important trait that influences its commercial value. It is known that quality traits such as anthocyanin composition and fruit surface pattern are categorical and are inherited simply. However, research examples of gene mapping for the composition (anthocyanin accumulation profile) and the surface pattern in eggplant fruit are limited. METHODS AND RESULTS: To map loci for these traits including the accumulation profiles of two anthocyanins, a widely spreading anthocyanin, delphinidin 3-(p-coumaroyl) rutinoside-5-glucoside (nasunin), and the relatively rare delphinidin 3-glucoside (D3G), we used two F2 intracrossed populations (LWF2 and N28F2). For the LWF2 population, mapping was achieved by reconstructing the linkage map created by Fukuoka et al. [1]. In the case of the N28F2 population, we constructed a linkage map consisting of 13 linkage groups using 238 simple sequence repeats, 75 single-nucleotide polymorphisms. Using the two F2 populations, the nasunin accumulating profile, the striped pattern on the fruit surface, the colors of flowers, fruit, and calyxes, and the D3G accumulating profile were genetically mapped. Furthermore, by utilizing the eggplant reference genome information, mutations in the causative candidate genes for those loci were identified. CONCLUSION: Overall, the results of this study suggest that inactivation of key enzymes of anthocyanin metabolism and the gene orthologous to the tomato u gene are potential causes of observed variety in eggplant appearance traits.

Shared genetic components between metabolic syndrome and schizophrenia: Genetic correlation using multipopulation data sets.

AIM: The genetic relationship between schizophrenia (SCZ) and other nonpsychiatric disorders remains largely unknown. We examined the shared genetic components between these disorders based on multipopulation data sets. METHODS: We used two data sets for East Asian (EAS) and European (EUR) samples. SCZ data was based on the Psychiatric Genomics Consortium Asia with our own genome-wide association study for EAS and Psychiatric Genomics Consortium for EUR. Nonpsychiatric data (20 binary traits [mainly nonpsychiatric complex disorders] and 34 quantitative traits [mainly laboratory examinations and physical characteristics]) were obtained from Biobank Japan and UK Biobank for EAS and EUR samples, respectively. To evaluate genetic correlation, linkage disequilibrium score regression analysis was utilized with further meta-analysis for each result from EAS and EUR samples to obtain robust evidence. Subsequent mendelian randomization analysis was also included to examine the causal effect. RESULTS: A significant genetic correlation between SCZ and several metabolic syndrome (MetS) traits was detected in the combined samples (meta-analysis between EAS and EUR data) (body mass index [rg  = -0.10, q-value = 1.0 × 10-9 ], high-density-lipoprotein cholesterol [rg  = 0.072, q-value = 2.9 × 10-3 ], blood sugar [rg  = -0.068, q-value = 1.4 × 10-2 ], triglycerides [rg  = -0.052, q-value = 2.4 × 10-2 ], systolic blood pressure [rg  = -0.054, q-value = 3.5 × 10-2 ], and C-reactive protein [rg  = -0.076, q-value = 7.8 × 10-5 ]. However, no causal relationship on SCZ susceptibility was detected for these traits based on the mendelian randomization analysis. CONCLUSION: Our results indicate shared genetic components between SCZ and MetS traits and C-reactive protein. Specifically, we found it interesting that the correlation between MetS traits and SCZ was the opposite of that expected from clinical studies: this genetic study suggests that SCZ susceptibility was associated with reduced MetS. This implied that MetS in patients with SCZ was not associated with genetic components but with environmental factors, including antipsychotics, lifestyle changes, poor diet, lack of exercise, and living conditions.

Polygenic risk score as clinical utility in psychiatry: a clinical viewpoint.

Genome-wide association studies (GWASs) have detected many susceptible variants for common diseases, including psychiatric disorders. However, because of the small effect size of each variant, clinical utility that aims for risk prediction and/or diagnostic assistance based on the individual "variants" is difficult to use. Therefore, to improve the statistical power, polygenic risk score (PRS) has been established and applied in the GWAS as a robust analytic tool. Although PRS has potential predictive ability, because of its current "insufficient" discriminative power at the individual level for clinical use, it remains limited solely in the research area, specifically in the psychiatric field. For a better understanding of the PRS, in this review, we (1) introduce the clinical features of psychiatric disorders, (2) summarize the recent GWAS/PRS findings in the psychiatric disorders, (3) evaluate the problems of PRS, and (4) propose its possible utility to apply PRS into the psychiatric clinical setting.

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.

Polygenic risk scores in schizophrenia with clinically significant copy number variants.

AIMS: Recent studies have revealed that the interplay between polygenic risk scores (PRS) and large copy number variants (CNV; >500kb) is essential for the etiology of schizophrenia (SCZ). To replicate previous findings, including those for smaller CNV (>10kb), the PRS between SCZ patients with and without CNV were compared. METHODS: The PRS were calculated for 724 patients with SCZ and 1178 healthy controls (HC), genotyped using array-based comparative genomic hybridization and single nucleotide polymorphisms chips, and comparisons were made between cases and HC, or between subjects with and without 'clinically significant' CNV. RESULTS: First, we replicated the higher PRS in patients with SCZ compared to that in HC (without taking into account the CNV). For clinically significant CNV, as defined by the American College of Medical Genetics ('pathogenic' and 'uncertain clinical significance, likely pathogenic' CNV), 66 patients with SCZ carried clinically significant CNV, whereas 658 SCZ patients had no such CNV. In the comparison of PRS between cases with/without the CNV, despite no significant difference in PRS, significant enrichment of the well-established risk CNV (22q11.2 deletion and 47,XXY/47,XXX) was observed in the lowest decile of PRS in SCZ patients with the CNV. CONCLUSION: Although the present study failed to replicate the significant difference in PRS between SCZ patients with and without clinically significant CNV, SCZ patients with well-established risk CNV tended to have a lower PRS. Therefore, we speculate that the CNV in SCZ patients with lower PRS may contain 'genuine' risk; PRS is a possible tool for prioritizing clinically significant CNV because the power of the CNV association analysis is limited due to their rarity.

Fine mapping of a major locus representing the lack of prickles in eggplant revealed the availability of a 0.5-kb insertion/deletion for marker-assisted selection.

As prickles cause labour inefficiency during cultivation and scratches on the skin of fruits during transportation, they are considered undesirable traits of eggplant (Solanum melongena L.). Because the molecular basis of prickle emergence has not been entirely revealed in plants, we mapped an eggplant semi-dominant Prickle (Pl) gene locus, which causes the absence of prickles, on chromosome 6 of a linkage map of the F2 population derived from crossing the no-prickly cultivar 'Togenashi-senryo-nigo' and the prickly line LS1934. By performing synteny mapping with tomato, the genomic region corresponding to the eggplant Pl locus was identified. Through bacterial artificial chromosome (BAC) screening, positive BAC clones and the contig sequence that harbour the Pl locus in the prickly eggplant genome were revealed. The BAC contig length was 133 kb, and it contained 16 predicted genes. Among them, a characteristic 0.5-kb insertion/deletion was detected. As the 0.5-kb insertion was commonly identified with the prickly phenotype worldwide, a primer pair that amplifies the insertion/deletion could be used for marker-assisted selection of the no-prickly phenotype. Such findings contribute to map-based-cloning of the Pl gene and the understanding of gene function, ultimately providing new insights into the regulatory molecular mechanisms underlying prickle emergence in plants.

Construction of a core collection of eggplant (Solanum melongena L.) based on genome-wide SNP and SSR genotypes.

A core collection of eggplant (Solanum melongena L.) was developed based on a dataset of genome-wide 831 SNP and 50 SSR genotypes analyzed in 893 accessions of eggplant genetic resources collected in the NARO Genebank using the Core Hunter II program. The 893 accessions were collected worldwide, mainly Asia. Genetic variation and population structure among the 893 eggplant accessions were characterized. The genetic diversity of the Asian accessions, especially the South Asian and Southeast Asian accessions, forming the center of diversity in eggplant, was higher than that of the other regions. The resulting core collection, World Eggplant Core (WEC) collection consisted of 100 accessions basically collected from the high genetic diversity countries. Based on the results of the cluster and STRUCTURE analyses with SNP genotypes, the WEC collection was divided into four clusters (S1-S4). Each cluster corresponds to a geographical group as below, S1; the European, American and African countries, S2; the East Asian countries, S3; the Southeast Asian countries, S4; the South Asian and Southeast Asian countries. The genotype and phenotype data of the WEC collection are available from the VegMarks database (https://vegmarks.nivot.affrc.go.jp/resource/), and seed samples are available from the NARO Genebank (https://www.gene.affrc.go.jp/databases-core_collections.php).

Genome-wide association studies of bipolar disorder: A systematic review of recent findings and their clinical implications.

Recent advances in molecular genetics have enabled assessments of the associations among genetic variants (e.g., single-nucleotide polymorphisms) and susceptibility for complex diseases, including psychiatric disorders. Specifically, genome-wide association studies (GWAS), meta-analyses of the GWAS summary statistics, and mega-analyses (which use raw data, not summary statistics) of GWAS have provided revolutionary results and have identified numerous susceptibility genes or single-nucleotide polymorphisms. By using several tens of thousands of subjects, >40 genes have been identified as being associated with susceptibility for bipolar disorder so far. The purpose of this systematic review was to summarize the recent findings of bipolar disorder GWAS and discuss their clinical implications.

DNA methylation analyses of the candidate genes identified by a methylome-wide association study revealed common epigenetic alterations in schizophrenia and bipolar disorder.

AIM: Schizophrenia (SZ) and bipolar disorder (BD) have been known to share genetic and environmental risk factors, and complex gene-environmental interactions may contribute to their pathophysiology. In contrast to high genetic overlap between SZ and BD, as revealed by genome-wide association studies, the extent of epigenetic overlap remains largely unknown. In the present study, we explored whether SZ and BD share epigenetic risk factors in the same manner as they share genetic components. METHODS: We performed DNA methylation analyses of the CpG sites in the top five candidate regions (FAM63B, ARHGAP26, CTAGE11P, TBC1D22A, and intergenic region [IR] on chromosome 16) reported in a previous methylome-wide association study (MWAS) of SZ, using whole blood samples from subjects with BD and controls. RESULTS: Among the five candidate regions, the CpG sites in FAM63B and IR on chromosome 16 were significantly hypomethylated in the samples from subjects with BD as well as those from subjects with SZ. On the other hand, the CpG sites in TBC1D22A were hypermethylated in the samples from subjects with BD, in contrast to hypomethylation in the samples from subjects with SZ. CONCLUSION: Hypomethylation of FAM63B and IR on chromosome 16 could be common epigenetic risk factors for SZ and BD. Further comprehensive epigenetic studies for BD, such as MWAS, will uncover the extent of similarity and uniqueness of epigenetic alterations.

Detailed mapping of a resistance locus against Fusarium wilt in cultivated eggplant (Solanum melongena).

KEY MESSAGE: This is the first report on genetic mapping of a resistance locus against Fusarium wilt caused by the plant pathogen Fusarium oxysporum f. sp. melongenae in cultivated eggplant. ABSTRACT: Fusarium wilt, caused by the plant pathogen Fusarium oxysporum f. sp. melongenae, is a major soil-borne disease threatening stable production in eggplant (Solanum melongena). Although three eggplant germplasms, LS1934, LS174, and LS2436, are known to be highly resistant to the pathogen, their resistance loci have not been mapped. In this study, we performed quantitative trait locus analyses in F2:3 populations and detected a resistance locus, FM1, at the end of chromosome 2, with two alleles, Fm1(L) and Fm1(E), in the F2 populations LWF2 [LS1934 × WCGR112-8 (susceptible)] and EWF2 [EPL-1 (derived from LS174) × WCGR112-8], respectively. The percentage of phenotypic variance explained by Fm1(L) derived from LS1934 was 75.0% [Logarithm of the odds (LOD) = 29.3], and that explained by Fm1(E) derived from EPL-1 was 92.2% (LOD = 65.8). Using backcrossed inbred lines, we mapped FM1 between two simple sequence repeat markers located ~4.881 cM apart from each other. Comparing the location of the above locus to those of previously reported ones, the resistance locus Rfo-sa1 from an eggplant ally (Solanum aethiopicum gr. Gilo) was mapped very close to FM1, whereas another resistance locus, from LS2436, was mapped to the middle of chromosome 4. This is the first report of mapping of a Fusarium resistance locus in cultivated eggplant. The availability of resistance-linked markers will enable the application of marker-assisted selection to overcome problems posed by self-incompatibility and introduction of negative traits because of linkage drag, and will lead to clear understanding of genetic mechanism of Fusarium resistance.

Risk factors and clinical characteristics of the depressive state induced by pegylated interferon therapy in patients with hepatitis C virus infection: A prospective study.

AIM: Pegylated interferon (PegIFN) therapies for hepatitis C virus (HCV) infection often induce a depressive state. This study aimed to identify the risk factors for and clinical characteristics of PegIFN-induced depressive state. METHODS: Sixty-nine subjects with HCV who received PegIFN therapy were enrolled. Before beginning therapy, all subjects were evaluated using the Neuroticism-Extraversion-Openness Five-Factor Inventory and the List of Threatening Events Questionnaire. Beck Depression Inventory (BDI) scores were also evaluated at baseline, 2-4 weeks after initiating therapy, and every 4 weeks thereafter. RESULTS: During the study, 18 subjects (24.3%) developed a depressive state (BDI ≥ 10). A bimodal peak of onset was observed during the early (2-8 weeks) and late (after 20 weeks) therapy phases. Moreover, we observed that baseline BDI scores (odds ratio [OR] = 1.40, P = 0.0104) and neuroticism (OR = 1.14, P = 0.0275) were significant risk factors for developing a depressive state. To determine the specific characteristics of this condition, we compared the BDI subscales between the 'PegIFN-induced' and 'general' depressive state reported previously. We found that the score at 'somatic symptoms' was higher in the 'PegIFN-induced' group. CONCLUSION: Our results indicate the following: (i) PegIFN-induced depressive state most frequently develops during the first 8 weeks of therapy; (ii) baseline BDI and neuroticism scores are risk factors for PegIFN-induced depressive state; and (iii) the core symptoms of PegIFN-induced depressive state are different from those of 'general' depression.

Rare UNC13B variations and risk of schizophrenia: Whole-exome sequencing in a multiplex family and follow-up resequencing and a case-control study.

Rare genomic variations inherited in multiplex schizophrenia families are suggested to play a role in the genetic etiology of the disease. To identify rare variations with large effects on the risk of developing schizophrenia, we performed whole-exome sequencing (WES) in two affected and one unaffected individual of a multiplex family with 10 affected individuals. We also performed follow-up resequencing of the unc-13 homolog B (Caenorhabditis elegans) (UNC13B) gene, a potential risk gene identified by WES, in the multiplex family and undertook a case-control study to investigate association between UNC13B and schizophrenia. UNC13B coding regions (39 exons) from 15 individuals of the multiplex family and 111 affected offspring for whom parental DNA samples were available were resequenced. Rare missense UNC13B variations identified by resequencing were further tested for association with schizophrenia in two independent case-control populations comprising a total of 1,753 patients and 1,602 controls. A rare missense variation (V1525M) in UNC13B was identified by WES in the multiplex family; this variation was present in five of six affected individuals, but not in eight unaffected individuals or one individual of unknown disease status. Resequencing UNC13B coding regions identified five rare missense variations (T103M, M813T, P1349T, I1362T, and V1525M). In the case-control study, there was no significant association between rare missense UNC13B variations and schizophrenia, although single-variant meta-analysis indicated that M813T was nominally associated with schizophrenia. These results do not support a contribution of rare missense UNC13B variations to the genetic etiology of schizophrenia in the Japanese population. © 2016 Wiley Periodicals, Inc.

Genetic association study between the detected risk variants based upon type II diabetes GWAS and psychotic disorders in the Japanese population.

Several epidemiological and genetic studies have suggested that the risk of type II diabetes (T2D) is likely to overlap with the susceptibility to psychotic disorders such as schizophrenia (SCZ) and bipolar disorder (BD). In this study, we aimed to examine the association of single-nucleotide polymorphisms (SNPs) detected in previous T2D genome-wide association studies (GWAS) with SCZ, BD and psychosis (SCZ plus BD). A total of 37 SNPs were selected from the literature. A two-stage analysis was conducted using a first set of screening samples (total N=3037) and a second set of replication samples (N=4950). None of the SNPs showed a significant association to the screening samples after correction for multiple testing. To avoid type II error, we genotyped the top three SNPs in BCL11A, HMG20A and HNF4A showing associations with any of the phenotypes (Puncorrected <0.01) using independent samples to replicate the nominal associations. However, we were unable to find any significant associations based on the screening results (Puncorrected>0.05). Our findings did not support the shared genetic risk between T2D and psychotic disorders in the Japanese population. However, further replication using a larger sample size is required.

Replication and cross-phenotype study based upon schizophrenia GWASs data in the Japanese population: support for association of MHC region with psychosis.

Recent genome-wide association studies (GWASs) of schizophrenia (SCZ) identified several susceptibility genes and suggested shared genetic components between SCZ and bipolar disorder (BD). We conducted a genetic association study of single nucleotide polymorphisms (SNPs) selected according to previous SCZ GWAS targeting psychotic disorders (SCZ and BD) in the Japanese population. Fifty-one SNPs were analyzed in a two-stage design using first-set screening samples (all SNPs: 1,032 SCZ, 1,012 BD, and 993 controls) and second-set replication samples ("significant" SNPs in the first-set screening analysis: 1,808 SCZ, 821 BD, and 2,321 controls). We assessed allelic associations between the selected SNPs and the three phenotypes (SCZ, BD, and "psychosis" [SCZ + BD]). Nine SNPs revealed nominal association signals for all comparisons (P(uncorrected) < 0.05), of which two SNPs located in the major histocompatibility complex region (rs7759855 in zinc finger and SCAN domain containing 31 [ZSCAN31] and rs1736913 in HLA-F antisense RNA1 [HLA-F-AS1]) were further assessed in the second-set replication samples. The associations were confirmed for rs7759855 (P(corrected) = 0.026 for psychosis; P(corrected) = 0.032 for SCZ), although the direction of effect was opposite to that in the original GWAS of the Chinese population. Finally, a meta-analysis was conducted using our two samples and using our data and data from Psychiatric GWAS Consortium (PGC), which have shown the same direction of effect. SNP in ZSCAN31 (rs7759855) had the strongest association with the phenotypes (best P = 6.8 × 10(-5) for psychosis: present plus PGC results). These data support shared risk SNPs between SCZ and BD in the Japanese population and association between MHC and psychosis.