Whole-exome sequencing of individuals from an isolated population implicates rare risk variants in bipolar disorder

Increased burden of rare protein-truncating variants in constrained, brain-specific and synaptic genes in extremely impulsively violent males with antisocial personality disorder

The genetic correlates of extreme impulsive violence are poorly understood, and there have been few studies that have characterized a large group of affected individuals both clinically and genetically. We performed whole exome sequencing (WES) in 290 males with the life-course-persistent, extremely impulsively violent form of antisocial personality disorder (APD) and analyzed the spectrum of rare protein-truncating variants (rPTVs). Comparisons were made with 314 male controls and publicly available genotype data. Functional annotation tools were used for biological interpretation. Participants were significantly more likely to harbor rPTVs in genes that are intolerant to loss-of-function variants (odds ratio [OR] 2.06; p < 0.001), specifically expressed in brain (OR 2.80; p = 0.036) and enriched for those involved in neurotransmitter transport and synaptic processes. In 60 individuals (20%), we identified rPTVs that we classified as clinically relevant based on their clinical associations, biological function and gene expression patterns. Of these, 37 individuals harbored rPTVs in 23 genes that are associated with a monogenic neurological disorder, and 23 individuals harbored rPTVs in 20 genes reportedly intolerant to loss-of-function variants. The analysis presents evidence in support of a model where presence of either one or several private, functionally relevant mutations contribute significantly to individual risk of life-course-persistent APD and reveals multiple individuals who could be affected by clinically unrecognized neuropsychiatric Mendelian disease. Thus, Mendelian diseases and increased rPTV burden may represent important factors for the development of extremely impulsive violent life-course-persistent forms of APD irrespective of their clinical presentation.

Bipolar Disorder

This review focuses on recent advances made towards understanding the neurobiology of bipolar disorder (BD), a chronic neuropsychiatric illness characterized by altered mood and energy states. The past few years have seen the completion of the largest genetic studies by far, which have emphasized the polygenic nature of BD as well as it's connection to other psychiatric illnesses. Furthermore, the use of inducible pluripotent stem cells has rapidly expanded. These studies support previous work that implicates dysregulation of neurodevelopment, mitochondria, and calcium homeostasis, while also allowing for investigation into the underlying mechanisms of individual responsivity to lithium. Sleep and circadian rhythms have also been heavily implicated in BD, from disruptions in activity patterns to molecular abnormalities.

SNP heterozygosity, relatedness and inbreeding of whole genomes from the isolated population of the Faroe Islands

BACKGROUND

The population of the Faroe Islands is an isolated population but very little is known about it from whole genome sequencing. The population of about 50000 people has a high incidence of rare diseases e.g., 1:300 for Primary Carnitine Deficiency. A screening programme was implemented, and eleven persons were also whole genome sequenced at x37 coverage for diagnostic purposes of those cases that were not affected by the known mutations. The purpose of our study is to utilize the high coverage data to explore the genomic variation and the ancestral history of the population. We study the SNP heterozygosity, the pairwise relatedness from kinship, the inbreeding from runs of homozygosity ROH, and we find the minor allele frequency distribution. We estimate the population ancestry and the timing of the founding event by using the whole genomes from eight consenting individuals.

RESULTS

We find the number of SNPs and the heterozygosity for the eight individual samples, and for merged samples, for which we also study the relatedness. We find close relatedness between the supposedly unrelated individuals. From ROH, we interpret the high relatedness as an ancient property of the isolated population. A bottleneck event is estimated starting between years [Formula: see text] with a maximum consanguineous population in year [Formula: see text] and similarly consanguineous between years [Formula: see text]. The ancestry analysis shows the population descends from founders of [Formula: see text] European and [Formula: see text] Admixed American ancestry. A distinct clustering near the central European and British populations of the 1000 Genome Project is likely the result of the population isolation and genetic drift. The minor allele frequency distribution suggests many rare variants.

CONCLUSIONS

The ancestry is mainly European while the inbreeding is higher compared to European populations and population isolates. The Faroese population has inbreeding more like ancient Europeans. We discovered a bottlenecked and consanguineous population event and estimated it starting in the 1st-4th century as compared to the oldest archaeological findings from the 4th-6th century.

Genetics of bipolar disorder: insights into its complex architecture and biology from common and rare variants

Bipolar disorder (BD) is a common mental disorder characterized by recurrent mood episodes, which causes major socioeconomic burdens globally. Though its disease pathogenesis is largely unknown, the high heritability of BD indicates strong contributions from genetic factors. In this review, we summarize the recent achievements in the genetics of BD, particularly those from genome-wide association study (GWAS) of common variants and next-generation sequencing analysis of rare variants. These include the identification of dozens of robust disease-associated loci, deepening of our understanding of the biology of BD, objective description of correlations with other psychiatric disorders and behavioral traits, formulation of methods for predicting disease risk and drug response, and the discovery of a single gene associated with bipolar disorder and schizophrenia spectrum with a large effect size. On the other hand, the findings to date have not yet made a clear contribution to the improvement of clinical psychiatry of BD. We overview the remaining challenges as well as possible paths to resolve them, referring to studies of other major neuropsychiatric disorders.

dbBIP: a comprehensive bipolar disorder database for genetic research

Bipolar disorder (BIP) is one of the most common hereditary psychiatric disorders worldwide. Elucidating the genetic basis of BIP will play a pivotal role in mechanistic delineation. Genome-wide association studies (GWAS) have successfully reported multiple susceptibility loci conferring BIP risk, thus providing insight into the effects of its underlying pathobiology. However, difficulties remain in the extrication of important and biologically relevant data from genetic discoveries related to psychiatric disorders such as BIP. There is an urgent need for an integrated and comprehensive online database with unified access to genetic and multi-omics data for in-depth data mining. Here, we developed the dbBIP, a database for BIP genetic research based on published data. The dbBIP consists of several modules, i.e.: (i) single nucleotide polymorphism (SNP) module, containing large-scale GWAS genetic summary statistics and functional annotation information relevant to risk variants; (ii) gene module, containing BIP-related candidate risk genes from various sources and (iii) analysis module, providing a simple and user-friendly interface to analyze one’s own data. We also conducted extensive analyses, including functional SNP annotation, integration (including summary-data-based Mendelian randomization and transcriptome-wide association studies), co-expression, gene expression, tissue expression, protein–protein interaction and brain expression quantitative trait loci analyses, thus shedding light on the genetic causes of BIP. Finally, we developed a graphical browser with powerful search tools to facilitate data navigation and access. The dbBIP provides a comprehensive resource for BIP genetic research as well as an integrated analysis platform for researchers and can be accessed online at http://dbbip.xialab.info.

Database URL:  http://dbbip.xialab.info

Genomic and neuroimaging approaches to bipolar disorder

BACKGROUND

To date, besides genome-wide association studies, a variety of other genetic analyses (e.g. polygenic risk scores, whole-exome sequencing and whole-genome sequencing) have been conducted, and a large amount of data has been gathered for investigating the involvement of common, rare and very rare types of DNA sequence variants in bipolar disorder. Also, non-invasive neuroimaging methods can be used to quantify changes in brain structure and function in patients with bipolar disorder.

AIMS

To provide a comprehensive assessment of genetic findings associated with bipolar disorder, based on the evaluation of different genomic approaches and neuroimaging studies.

METHOD

We conducted a PubMed search of all relevant literatures from the beginning to the present, by querying related search strings.

RESULTS

ANK3, CACNA1C, SYNE1, ODZ4 and TRANK1 are five genes that have been replicated as key gene candidates in bipolar disorder pathophysiology, through the investigated studies. The percentage of phenotypic variance explained by the identified variants is small (approximately 4.7%). Bipolar disorder polygenic risk scores are associated with other psychiatric phenotypes. The ENIGMA-BD studies show a replicable pattern of lower cortical thickness, altered white matter integrity and smaller subcortical volumes in bipolar disorder.

CONCLUSIONS

The low amount of explained phenotypic variance highlights the need for further large-scale investigations, especially among non-European populations, to achieve a more complete understanding of the genetic architecture of bipolar disorder and the missing heritability. Combining neuroimaging data with genetic data in large-scale studies might help researchers acquire a better knowledge of the engaged brain regions in bipolar disorder.

Heterozygous De Novo Truncating Mutation of Nucleolin in an ASD Individual Disrupts Its Nucleolar Localization

Nucleolin (NCL/C23; OMIM: 164035) is a major nucleolar protein that plays a critical role in multiple processes, including ribosome assembly and maturation, chromatin decondensation, and pre-rRNA transcription. Due to its diverse functions, nucleolin has frequently been implicated in pathological processes, including cancer and viral infection. We recently identified a de novo frameshifting indel mutation of NCL, p.Gly664Glufs*70, through whole-exome sequencing of autism spectrum disorder trios. Through the transfection of constructs encoding either a wild-type human nucleolin or a mutant nucleolin with the same C-terminal sequence predicted for the autism proband, and by using co-localization with the nucleophosmin (NPM; B23) protein, we have shown that the nucleolin mutation leads to mislocalization of the NCL protein from the nucleolus to the nucleoplasm. Moreover, a construct with a nonsense mutation at the same residue, p.Gly664*, shows a very similar effect on the location of the NCL protein, thus confirming the presence of a predicted nucleolar location signal in this region of the NCL protein. Real-time fluorescence recovery experiments show significant changes in the kinetics and mobility of mutant NCL protein in the nucleoplasm of HEK293Tcells. Several other studies also report de novoNCL mutations in ASD or neurodevelopmental disorders. The altered mislocalization and dynamics of mutant NCL (p.G664Glufs*70/p.G664*) may have relevance to the etiopathlogy of NCL-related ASD and other neurodevelopmental phenotypes.

Emerging perspectives on multidomain phosphatidylinositol transfer proteins

The phosphatidylinositol transfer protein domain (PITPd) is an evolutionarily conserved protein that is able to transfer phosphatidylinositol between membranes in vitro and in vivo. However some animal genomes also include genes that encode proteins where the PITPd is found in cis with a number of additional domains and recent large scale genome sequencing efforts indicate that this type of multidomain architecture is widespread in the animal kingdom. In Drosophila photoreceptors, the multidomain phosphatidylinositol transfer protein RDGB is required to regulate phosphoinositide turnover during G-protein activated phospholipase C signalling. Recent studies in flies and mammalian cell culture models have begun to elucidate functions for the non-PITPd of RDGB and its vertebrate orthologs. We review emerging evidence on the genomics, functional and cell biological perspectives of these multi-domain PITPd containing proteins.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

Integration analysis of methylation quantitative trait loci and GWAS identify three schizophrenia risk variants

Genome-wide association studies (GWAS) have identified hundreds of genetic variants associated with schizophrenia (SCZ). However, prioritizing risk variants and regulatory elements for follow-up functional studies remains a major challenge. Therefore, we performed an integrated analysis to identify variants who affect methylation levels of nearby genes and contribute to the risk of SCZ, and to explore the potential role of these variants in SCZ pathogenesis. First, we used the Summary data-based Mendelian Randomization (SMR) method to integrate GWAS and methylation quantitative trait loci data. Then, the SNP-methylation combinations as associated with SCZ were replicated across multiple samples. Totally, we identified and replicated 14 and one SNP-methylation combinations in blood and brain tissues, respectively, that significantly associated with SCZ. Furthermore, our expression quantitative trait loci analysis, differential methylation analysis, neuroimaging genetics, and cognitive genetics analysis consistently supported the potential roles of these 15 SNPs in the pathogenesis of SCZ. Finally, using the convergent functional genomics method, we prioritized three risk SNPs, including rs3765971 (RERE, P_SMR = 3.87 × 10^-8), rs55742290 (ARL6IP4, P_SMR = 1.50 × 10^-7), and rs7293091 (CENPM, P_SMR = 5.09 × 10^-7), may represent promising risk variants in SCZ. These convergent lines of evidence suggest that three risk variants may be involved in the pathogenesis of SCZ. Further investigation of the roles of these variants in the pathogenesis of SCZ is warranted.

The black sheep of the family- whole-exome sequencing in family of lithium response discordant bipolar monozygotic twins

Twin studies are among the most promising strategies for studying heritable disorders, including bipolar disorder (BD). The aim of the present study was to identify distinguishing genes between monozygotic (MZ) twins with different BD phenotype and compare them to their non-affected siblings. Whole-exome sequencing (WES) can identify rare and structural variants that could detect the polygenetic burden of complex disorders. WES was performed on a family composed of two MZ twins with BD, their unaffected brother and unaffected parents. The twins have a discordant response to lithium and distinct course of illness. Following WES, six genes of particular interest emerged: Neurofibromin type 1 (NF1), Biorientation of chromosomes in cell division 1 (BOD1), Golgi-associated gamma adaptin ear-containing ARF binding protein 3 (GGA3), Disrupted in schizophrenia 1 (DISC1), Neuromedin U receptor 2 (NMUR2), and Huntingtin interacting protein 1-related (HIP1R). Interestingly, many of these influence glutamatergic pathways and thus the findings may have therapeutical implications. These results may provide important insights to unveil genetic underpinnings of BD and the response to lithium.

Contribution of common and rare damaging variants in familial forms of bipolar disorder and phenotypic outcome

Genome-wide association studies on bipolar disorders (BD) have revealed an additive polygenic contribution of common single-nucleotide polymorphisms (SNPs). However, these SNPs explain only 25% of the overall genetic variance and suggest a role of rare variants in BD vulnerability. Here, we combined high-throughput genotyping data and whole-exome sequencing in cohorts of individuals with BD as well as in multiplex families with a high density of affected individuals in order to determine the contribution of both common and rare variants to BD genetic vulnerability. Using polygenic risk scores (PRS), we showed a strong contribution of common polymorphisms previously associated with BD and schizophrenia (SZ) and noticed that those specifically associated with SZ contributed more in familial forms of BD than in non-familial ones. The analysis of rare damaging variants shared by affected individuals in multiplex families with BD revealed a single interaction network enriched in neuronal and developmental biological pathways, as well as in the regulation of gene expression. We identified four genes with a higher mutation rate in individuals with BD than in the general population and showed that mutations in two of them were associated with specific clinical manifestations. In addition, we showed a significant negative correlation between PRS and the number of rare damaging variants specifically in unaffected individuals of multiplex families. Altogether, our results suggest that common and rare genetic variants both contribute to the familial aggregation of BD and this genetic architecture may explain the heterogeneity of clinical manifestations in multiplex families.

Association between CANCA1C gene rs1034936 polymorphism and alcohol dependence in bipolar disorder

INTRODUCTION

Bipolar disorder (BD) is a highly heritable and disabling mental illness, commonly associated with substance abuse, being alcohol abuse the most frequent. Comorbid BD and substance abuse disorders are often associated with high levels of health service utilization and destabilization of the course of illness resulting in poor treatment outcomes. Although recent genome-wide association studies have detected a number of risk genes for BD, the data is still sparse and inconclusive for those genes that may contribute to the increased risk of comorbid alcohol abuse (AA) in BD. The primary aim of the present study was to investigate the effects of 46 single-nucleotide polymorphisms (SNPs) within eight genes on different phenotypes of BD patients, such as comorbid alcohol abuse. We further assessed clinical variables associated with AA.

METHODS

One-hundred fifty-eight BD I and II patients were enrolled in a naturalistic cohort study. Genomic DNA of 92 patients was extracted from whole blood using standard procedures and 46 tag SNPs in eight genes of interest (ANK, CACNA1C, CACNB2, FKBP5, GRM7, ITIH3, SYNE1 and TCF4) were genotyped.

RESULTS

Seventy-one patients out of 158 (45%) satisfied diagnostic criteria for comorbid AA. Among 46 SNPs analyzed, the only SNP associated with comorbid AA was rs1034936 polymorphism in the CANCA1C gene. This polymorphism was also associated with lifetime cocaine abuse, manic switch and current atypical antipsychotics.

CONCLUSIONS

Our findings suggest a role of rs1034936 CACNA1C gene variant in BD-AA group. Despite their preliminary nature, the present results may provide new insight on mechanisms underlying AA in BD.

Whole-exome sequencing identifies variants associated with structural MRI markers in patients with bipolar disorders

BACKGROUND

Bipolar disorder (BD) is one of the most heritable psychiatric disorders. A growing number of whole-exome sequencing (WES) studies for BD has been performed, however, no research has examined the association between single nucleotide variants (SNVs) from WES and structural magnetic resonance imaging (MRI) data.

METHODS

We sequenced whole-exomes in 53 patients with BD and 82 healthy control participants at an initial discovery stage and investigated the impacts of SNVs in risk genes from WES analysis on the cortical gray-matter thickness and integrity of white matter tracts and in the following stage. Cortical thickness and white matter integrity were investigated using the FreeSurfer and TRACULA (Tracts Constrained by UnderLying Anatomy).

RESULTS

We identified 122 BD-related genes including KMT2C, AHNAK, CDH23, DCHS1, FRAS1, MACF1 and RYR3 and observed 27 recurrent copy number alteration regions including gain on 8p23.1 and loss on 15q11.1 - q11.2. Among them, single nucleotide polymorphism (SNP) rs4639425 in KMT2C gene, which regulates histone H3 lysine 4 (H3K4) methylation involved in chromatin remodeling, was associated with widespread alterations of white matter integrity including the cingulum, uncinate fasciculus, cortico-spinal tract, and superior longitudinal fasciculus.

LIMITATION

The small sample size of patients with BD in the genome data may cause our study to be underpowered when searching for putative rare mutations.

CONCLUSION

This study first combined a WES approach and neuroimaging findings in psychiatric disorders. We postulate the rs4639425 may be associated with BD-related microstructural changes of white matter tracts.

Nitric oxide signalling and antidepressant action revisited

Studies about the pathogenesis of mood disorders have consistently shown that multiple factors, including genetic and environmental, play a crucial role on their development and neurobiology. Multiple pathological theories have been proposed, of which several ultimately affects or is a consequence of dysfunction in brain neuroplasticity and homeostatic mechanisms. However, current clinical available pharmacological intervention, which is predominantly monoamine-based, suffers from a partial and lacking response even after weeks of continuous treatment. These issues raise the need for better understanding of aetiologies and brain abnormalities in depression, as well as developing novel treatment strategies. Nitric oxide (NO) is a gaseous unconventional neurotransmitter, which regulates and governs several important physiological functions in the central nervous system, including processes, which can be associated with the development of mood disorders. This review will present general aspects of the NO system in depression, highlighting potential targets that may be utilized and further explored as novel therapeutic targets in the future pharmacotherapy of depression. In particular, the review will link the importance of neuroplasticity mechanisms governed by NO to a possible molecular basis for the antidepressant effects.

A preliminary investigation of rare variants associated with genetic risk for PTSD in a natural disaster-exposed adolescent sample

Background: Posttraumatic stress disorder (PTSD) involves a complex interaction of biological, psychological, and social factors. Numerous studies have demonstrated genetic variation associated with the development of PTSD, primarily in adults. However, the contribution of low frequency and rare genetic variants to PTSD is unknown to date. Moreover, there is limited work on genetic risk for PTSD in child and adolescent populations. Objective: This preliminary study aimed to identify the low frequency and rare genetic variation that contributes to PTSD using an exome array. Method: This post-disaster, adolescent sample (n = 707, 51% females, M _age = 14.54) was assessed for PTSD diagnosis and symptom count following tornado exposure. Results: Gene-based models, covarying for ancestry principal components, age, sex, tornado severity, and previous trauma identified variants in four genes associated with diagnosis and 276 genes associated with symptom count (at p _adj < .001). Functional class analyses suggested an association with variants in the nonsense class (nonsynonymous variant that results in truncation of, and usually non-functional, protein) with both outcomes. An exploratory gene network pathway analysis showed a great number of significant genes involved in brain and immune function, illustrating the usefulness of downstream examination of gene-based findings that may point to relevant biological processes. Conclusions: While further investigation in larger samples is warranted, findings align with extant PTSD literature that has identified variants associated with biological conditions such as immune function.

Rare Risk Variants Identification by Identity-by-Descent Mapping and Whole-Exome Sequencing Implicates Neuronal Development Pathways in Schizophrenia and Bipolar Disorder

Schizophrenia (SCZ) and bipolar disorder (BPD) are highly heritable disorders with an estimated co-heritability of 68%. Hundreds of common alleles have been implicated, but recently a role for rare, high-penetrant variants has been also suggested in both disorders. This study investigated a familial cohort of SCZ and BPD patients from a closed population sample, where the high recurrence of the disorders and the homogenous genetic background indicate a possible enrichment in rare risk alleles. A total of 230 subjects (161 cases, 22 unaffected relatives, and 47 controls) were genetically investigated through an innovative strategy that integrates identity-by-descent (IBD) mapping and whole-exome sequencing (WES). IBD analysis allowed to track high-risk haplotypes (IBD_risk) shared exclusively by multiple patients from different families and possibly carrying the most penetrant alleles. A total of 444 non-synonymous sequence variants, of which 137 disruptive, were identified in IBD_risk haplotypes by WES. Interestingly, gene sets previously implicated in SCZ (i.e., post-synaptic density (PSD) proteins, voltage-gated calcium channels (VGCCs), and fragile X mental retardation protein (FMRP) targets) were found significantly enriched in genes carrying IBD_risk variants. Further, IBD_risk variants were preferentially affecting genes involved in the extracellular matrix (ECM) biology and axon guidance processes which appeared to be functionally connected in the pathway-derived meta-network analysis. Results thus confirm rare risk variants as key factors in SCZ and BPD pathogenesis and highlight a role for the development of neuronal connectivity in the etiology of both disorders.

NGS Technologies as a Turning Point in Rare Disease Research , Diagnosis and Treatment

Approximately 25-50 million Americans, 30 million Europeans, and 8% of the Australian population have a rare disease. Rare diseases are thus a common problem for clinicians and account for enormous healthcare costs worldwide due to the difficulty of establishing a specific diagnosis. In this article, we review the milestones achieved in our understanding of rare diseases since the emergence of next-generation sequencing (NGS) technologies and analyze how these advances have influenced research and diagnosis. The first half of this review describes how NGS has changed diagnostic workflows and provided an unprecedented, simple way of discovering novel disease-associated genes. We focus particularly on metabolic and neurodevelopmental disorders. NGS has enabled cheap and rapid genetic diagnosis, highlighted the relevance of mosaic and de novo mutations, brought to light the wide phenotypic spectrum of most genes, detected digenic inheritance or the presence of more than one rare disease in the same patient, and paved the way for promising new therapies. In the second part of the review, we look at the limitations and challenges of NGS, including determination of variant causality, the loss of variants in coding and non-coding regions, and the detection of somatic mosaicism variants and epigenetic mutations, and discuss how these can be overcome in the near future.

Weighted functional linear regression models for gene-based association analysis

Functional linear regression models are effectively used in gene-based association analysis of complex traits. These models combine information about individual genetic variants, taking into account their positions and reducing the influence of noise and/or observation errors. To increase the power of methods, where several differently informative components are combined, weights are introduced to give the advantage to more informative components. Allele-specific weights have been introduced to collapsing and kernel-based approaches to gene-based association analysis. Here we have for the first time introduced weights to functional linear regression models adapted for both independent and family samples. Using data simulated on the basis of GAW17 genotypes and weights defined by allele frequencies via the beta distribution, we demonstrated that type I errors correspond to declared values and that increasing the weights of causal variants allows the power of functional linear models to be increased. We applied the new method to real data on blood pressure from the ORCADES sample. Five of the six known genes with P < 0.1 in at least one analysis had lower P values with weighted models. Moreover, we found an association between diastolic blood pressure and the VMP1 gene (P = 8.18×10-6), when we used a weighted functional model. For this gene, the unweighted functional and weighted kernel-based models had P = 0.004 and 0.006, respectively. The new method has been implemented in the program package FREGAT, which is freely available at https://cran.r-project.org/web/packages/FREGAT/index.html.

A Genome-Wide Association Study and Complex Network Identify Four Core Hub Genes in Bipolar Disorder

Bipolar disorder is a common and severe mental illness with unsolved pathophysiology. A genome-wide association study (GWAS) has been used to find a number of risk genes, but it is difficult for a GWAS to find genes indirectly associated with a disease. To find core hub genes, we introduce a network analysis after the GWAS was conducted. Six thousand four hundred fifty eight single nucleotide polymorphisms (SNPs) with p < 0.01 were sifted out from Wellcome Trust Case Control Consortium (WTCCC) dataset and mapped to 2045 genes, which are then compared with the protein–protein network. One hundred twelve genes with a degree >17 were chosen as hub genes from which five significant modules and four core hub genes (FBXL13, WDFY2, bFGF, and MTHFD1L) were found. These core hub genes have not been reported to be directly associated with BD but may function by interacting with genes directly related to BD. Our method engenders new thoughts on finding genes indirectly associated with, but important for, complex diseases.