Identification of Susceptible Loci and Enriched Pathways for Bipolar II Disorder Using Genome-Wide Association Studies

Genetic architecture of the structural connectome

Myelinated axons form long-range connections that enable rapid communication between distant brain regions, but how genetics governs the strength and organization of these connections remains unclear. We perform genome-wide association studies of 206 structural connectivity measures derived from diffusion magnetic resonance imaging tractography of 26,333 UK Biobank participants, each representing the density of myelinated connections within or between a pair of cortical networks, subcortical structures or cortical hemispheres. We identify 30 independent genome-wide significant variants after Bonferroni correction for the number of measures studied (126 variants at nominal genome-wide significance) implicating genes involved in myelination (SEMA3A), neurite elongation and guidance (NUAK1, STRN, DPYSL2, EPHA3, SEMA3A, HGF, SHTN1), neural cell proliferation and differentiation (GMNC, CELF4, HGF), neuronal migration (CCDC88C), cytoskeletal organization (CTTNBP2, MAPT, DAAM1, MYO16, PLEC), and brain metal transport (SLC39A8). These variants have four broad patterns of spatial association with structural connectivity: some have disproportionately strong associations with corticothalamic connectivity, interhemispheric connectivity, or both, while others are more spatially diffuse. Structural connectivity measures are highly polygenic, with a median of 9.1 percent of common variants estimated to have non-zero effects on each measure, and exhibited signatures of negative selection. Structural connectivity measures have significant genetic correlations with a variety of neuropsychiatric and cognitive traits, indicating that connectivity-altering variants tend to influence brain health and cognitive function. Heritability is enriched in regions with increased chromatin accessibility in adult oligodendrocytes (as well as microglia, inhibitory neurons and astrocytes) and multiple fetal cell types, suggesting that genetic control of structural connectivity is partially mediated by effects on myelination and early brain development. Our results indicate pervasive, pleiotropic, and spatially structured genetic control of white-matter structural connectivity via diverse neurodevelopmental pathways, and support the relevance of this genetic control to healthy brain function.

Admixture mapping implicates 13q33.3 as ancestry-of-origin locus for Alzheimer disease in Hispanic and Latino populations

Alzheimer disease (AD) is the most common form of senile dementia, with high incidence late in life in many populations including Caribbean Hispanic (CH) populations. Such admixed populations, descended from more than one ancestral population, can present challenges for genetic studies, including limited sample sizes and unique analytical constraints. Therefore, CH populations and other admixed populations have not been well represented in studies of AD, and much of the genetic variation contributing to AD risk in these populations remains unknown. Here, we conduct genome-wide analysis of AD in multiplex CH families from the Alzheimer Disease Sequencing Project (ADSP). We developed, validated, and applied an implementation of a logistic mixed model for admixture mapping with binary traits that leverages genetic ancestry to identify ancestry-of-origin loci contributing to AD. We identified three loci on chromosome 13q33.3 associated with reduced risk of AD, where associations were driven by Native American (NAM) ancestry. This AD admixture mapping signal spans the FAM155A, ABHD13, TNFSF13B, LIG4, and MYO16 genes and was supported by evidence for association in an independent sample from the Alzheimer's Genetics in Argentina-Alzheimer Argentina consortium (AGA-ALZAR) study with considerable NAM ancestry. We also provide evidence of NAM haplotypes and key variants within 13q33.3 that segregate with AD in the ADSP whole-genome sequencing data. Interestingly, the widely used genome-wide association study approach failed to identify associations in this region. Our findings underscore the potential of leveraging genetic ancestry diversity in recently admixed populations to improve genetic mapping, in this case for AD-relevant loci.

Neutral sphingomyelinase 2 inhibition attenuates extracellular vesicle release and improves neurobehavioral deficits in murine HIV

People living with HIV (PLH) have significantly higher rates of cognitive impairment (CI) and major depressive disorder (MDD) versus the general population. The enzyme neutral sphingomyelinase 2 (nSMase2) is involved in the biogenesis of ceramide and extracellular vesicles (EVs), both of which are dysregulated in PLH, CI, and MDD. Here we evaluated EcoHIV-infected mice for behavioral abnormalities relevant to depression and cognition deficits, and assessed the behavioral and biochemical effects of nSMase2 inhibition. Mice were infected with EcoHIV and daily treatment with either vehicle or the nSMase2 inhibitor (R)-(1-(3-(3,4-dimethoxyphenyl)-2,6-dimethylimidazo[1,2-b]pyridazin-8-yl)pyrrolidin-3-yl)-carbamate (PDDC) began 3 weeks post-infection. After 2 weeks of treatment, mice were subjected to behavior tests. EcoHIV-infected mice exhibited behavioral abnormalities relevant to MDD and CI that were reversed by PDDC treatment. EcoHIV infection significantly increased cortical brain nSMase2 activity, resulting in trend changes in sphingomyelin and ceramide levels that were normalized by PDDC treatment. EcoHIV-infected mice also exhibited increased levels of brain-derived EVs and altered microRNA cargo, including miR-183-5p, miR-200c-3p, miR-200b-3p, and miR-429-3p, known to be associated with MDD and CI; all were normalized by PDDC. In conclusion, inhibition of nSMase2 represents a possible new therapeutic strategy for the treatment of HIV-associated CI and MDD.

Low peripheral mitochondrial DNA copy number during manic episodes of bipolar disorders is associated with disease severity and inflammation

Mitochondria (Mt) are intra-cellular components essential for cellular energy processes whose dysfunction may induce premature cellular senescence and/or inflammation, both observed in bipolar disorders (BD). We investigated mitochondrial DNA copy number (mtDNAcn) levels in patients with BD being in manic, depressive or euthymic phase and in healthy controls (HC) both characterized for the levels of blood-based inflammatory markers and stigma of pathogens. 312 patients with BD were compared to 180 HC. mtDNAcn were measured using a digital droplet PCR. Serum levels of 14 inflammatory molecules and 3 anti-infectious IgG stigma were respectively evaluated by electro-chemiluminescence, ELISA and dedicated immunoassays. The statistical analyses were performed using Spearman's correlation, Wilcoxon signed-rank and Kruskal-Wallis rank sum tests. P-values were adjusted for multiple testing with Benjamini-Hochberg method. We found low levels of mtDNAcn in BD patients as compared to HC (P = 0.008) especially during manic episodes (P = 0.0002). We also observed that low levels of mtDNAcn are negatively correlated with mood and psychotic scales (PANSS, YMRS and CGI) (adjusted P (Adj P) = 0.02, 0.003 and 0.05 respectively) and positively with the GAF severity scale (Adj P = 0.002). They were also correlated with high levels of both intercellular adhesion molecule (ICAM)-1 and vascular cell adhesion molecule (VCAM)-1 (Adj P = 0.003 and 0.001) along with a trend toward increased IL-2, IL-10 and B2M circulating levels (Adj P = 0.05). Here, we report correlations between marker of mitochondria functioning and both clinical scales and inflammatory markers in BD patients experiencing manic episodes. If replicated, these finding might allow to predict transition between disease phases and to design accurate therapeutic options.

Genetic contributions to bipolar disorder: current status and future directions

Bipolar disorder (BD) is a highly heritable mental disorder and is estimated to affect about 50 million people worldwide. Our understanding of the genetic etiology of BD has greatly increased in recent years with advances in technology and methodology as well as the adoption of international consortiums and large population-based biobanks. It is clear that BD is also highly heterogeneous and polygenic and shows substantial genetic overlap with other psychiatric disorders. Genetic studies of BD suggest that the number of associated loci is expected to substantially increase in larger future studies and with it, improved genetic prediction of the disorder. Still, a number of challenges remain to fully characterize the genetic architecture of BD. First among these is the need to incorporate ancestrally-diverse samples to move research away from a Eurocentric bias that has the potential to exacerbate health disparities already seen in BD. Furthermore, incorporation of population biobanks, registry data, and electronic health records will be required to increase the sample size necessary for continued genetic discovery, while increased deep phenotyping is necessary to elucidate subtypes within BD. Lastly, the role of rare variation in BD remains to be determined. Meeting these challenges will enable improved identification of causal variants for the disorder and also allow for equitable future clinical applications of both genetic risk prediction and therapeutic interventions.

The C-terminal tail extension of myosin 16 acts as a molten globule, including intrinsically disordered regions, and interacts with the N-terminal ankyrin

The lesser-known unconventional myosin 16 protein is essential in proper neuronal functioning and has been implicated in cell cycle regulation. Its longer Myo16b isoform contains a C-terminal tail extension (Myo16Tail), which has been shown to play a role in the neuronal phosphoinositide 3-kinase signaling pathway. Myo16Tail mediates the actin cytoskeleton remodeling, downregulates the actin dynamics at the postsynaptic site of dendritic spines, and is involved in the organization of the presynaptic axon terminals. However, the functional and structural features of this C-terminal tail extension are not well known. Here, we report the purification and biophysical characterization of the Myo16Tail by bioinformatics, fluorescence spectroscopy, and CD. Our results revealed that the Myo16Tail is functionally active and interacts with the N-terminal ankyrin domain of myosin 16, suggesting an intramolecular binding between the C and N termini of Myo16 as an autoregulatory mechanism involving backfolding of the motor domain. In addition, the Myo16Tail possesses high structural flexibility and a solvent-exposed hydrophobic core, indicating the largely unstructured, intrinsically disordered nature of this protein region. Some secondary structure elements were also observed, indicating that the Myo16Tail likely adopts a molten globule-like structure. These structural features imply that the Myo16Tail may function as a flexible display site particularly relevant in post-translational modifications, regulatory functions such as backfolding, and phosphoinositide 3-kinase signaling.

Aging and age-related diseases: from mechanisms to therapeutic strategies

Aging is a physiological process mediated by numerous biological and genetic pathways, which are directly linked to lifespan and are a driving force for all age-related diseases. Human life expectancy has greatly increased in the past few decades, but this has not been accompanied by a similar increase in their healthspan. At present, research on aging biology has focused on elucidating the biochemical and genetic pathways that contribute to aging over time. Several aging mechanisms have been identified, primarily including genomic instability, telomere shortening, and cellular senescence. Aging is a driving factor of various age-related diseases, including neurodegenerative diseases, cardiovascular diseases, cancer, immune system disorders, and musculoskeletal disorders. Efforts to find drugs that improve the healthspan by targeting the pathogenesis of aging have now become a hot topic in this field. In the present review, the status of aging research and the development of potential drugs for aging-related diseases, such as metformin, rapamycin, resveratrol, senolytics, as well as caloric restriction, are summarized. The feasibility, side effects, and future potential of these treatments are also discussed, which will provide a basis to develop novel anti-aging therapeutics for improving the healthspan and preventing aging-related diseases.

Identifying age-specific gene signatures of the human cerebral cortex with joint analysis of transcriptomes and functional connectomes

The human cerebral cortex undergoes profound structural and functional dynamic variations across the lifespan, whereas the underlying molecular mechanisms remain unclear. Here, with a novel method transcriptome-connectome correlation analysis (TCA), which integrates the brain functional magnetic resonance images and region-specific transcriptomes, we identify age-specific cortex (ASC) gene signatures for adolescence, early adulthood and late adulthood. The ASC gene signatures are significantly correlated with the cortical thickness (P-value <2.00e-3) and myelination (P-value <1.00e-3), two key brain structural features that vary in accordance with brain development. In addition to the molecular underpinning of age-related brain functions, the ASC gene signatures allow delineation of the molecular mechanisms of neuropsychiatric disorders, such as the regulation between ARNT2 and its target gene ETF1 involved in Schizophrenia. We further validate the ASC gene signatures with published gene sets associated with the adult cortex, and confirm the robustness of TCA on other brain image datasets. Availability: All scripts are written in R. Scripts for the TCA method and related statistics result can be freely accessed at https://github.com/Soulnature/TCA. Additional data related to this paper may be requested from the authors.

Myosin XVI in the Nervous System

The myosin family is a large inventory of actin-associated motor proteins that participate in a diverse array of cellular functions. Several myosin classes are expressed in neural cells and play important roles in neural functioning. A recently discovered member of the myosin superfamily, the vertebrate-specific myosin XVI (Myo16) class is expressed predominantly in neural tissues and appears to be involved in the development and proper functioning of the nervous system. Accordingly, the alterations of MYO16 has been linked to neurological disorders. Although the role of Myo16 as a generic actin-associated motor is still enigmatic, the N-, and C-terminal extensions that flank the motor domain seem to confer unique structural features and versatile interactions to the protein. Recent biochemical and physiological examinations portray Myo16 as a signal transduction element that integrates cell signaling pathways to actin cytoskeleton reorganization. This review discusses the current knowledge of the structure-function relation of Myo16. In light of its prevalent localization, the emphasis is laid on the neural aspects.

The genetics of bipolar disorder

Bipolar disorder (BD) is one of the most heritable mental illnesses, but the elucidation of its genetic basis has proven to be a very challenging endeavor. Genome-Wide Association Studies (GWAS) have transformed our understanding of BD, providing the first reproducible evidence of specific genetic markers and a highly polygenic architecture that overlaps with that of schizophrenia, major depression, and other disorders. Individual GWAS markers appear to confer little risk, but common variants together account for about 25% of the heritability of BD. A few higher-risk associations have also been identified, such as a rare copy number variant on chromosome 16p11.2. Large scale next-generation sequencing studies are actively searching for other alleles that confer substantial risk. As our understanding of the genetics of BD improves, there is growing optimism that some clear biological pathways will emerge, providing a basis for future studies aimed at molecular diagnosis and novel therapeutics.

Serum miRNA as a possible biomarker in the diagnosis of bipolar II disorder

The diagnosis of Bipolar II disorder (BD-II) is currently based on the patients' description of symptoms and clinical behavioral observations. This study explored the possibility of miRNA in peripheral blood (serum) as a specific biomarker for BD-II. We identified 6 candidate miRNAs to differentiate BD-II patients from controls using next-generation sequencing. We then examined these candidate miRNAs using real-time PCR in the first cohort (as training group) of 79 BD-II and 95 controls. A diagnostic model was built based on these candidate miRNAs and then tested on an individual testing group (BD-II: n = 20, controls: n = 20). We found that serum expression levels of miR-7-5p, miR-23b-3p, miR-142-3p, miR-221-5p, and miR-370-3p significantly increased in BD-II compared with controls in the first cohort, whereas that of miR-145-5p showed no significant difference. The diagnostic power of the identified miRNAs was further analyzed using receiver-operating characteristic (ROC). Support vector machine (SVM) measurements revealed that a combination of the significant miRNAs reached good diagnostic accuracy (AUC: 0.907). We further examined an independent testing group and the diagnostic power reached fair for BD-II (specificity = 90%, sensitivity = 85%). We constructed miRNA panels using SVM model, which may aid in the diagnosis for BD-II.

Myosin XVI

Myosin XVI (Myo16), a vertebrate-specific motor protein, is a recently discovered member of the myosin superfamily. The detailed functionality regarding myosin XVI requires elucidating or clarification; however, it appears to portray an important role in neural development and in the proper functioning of the nervous system. It is expressed in the largest amount in neural tissues in the late embryonic-early postnatal period, specifically the time in which neuronal cell migration and dendritic elaboration coincide. The impaired expression of myosin XVI has been found lurking in the background of several neuropsychiatric disorders including autism, schizophrenia and/or bipolar disorders.Two principal isoforms of class XVI myosins have been thus far described: Myo16a, the tailless cytoplasmic isoform and Myo16b, the full-length molecule featuring both cytoplasmic and nuclear localization. Both isoforms contain a class-specific N-terminal ankyrin repeat domain that binds to the protein phosphatase catalytic subunit. Myo16b, the predominant isoform, exhibits a diverse function. In the cytoplasm, it participates in the reorganization of the actin cytoskeleton through activation of the PI3K pathway and the WAVE-complex, while in the nucleus it may possess a role in cell cycle regulation. Based on the sequence, myosin XVI may have a compromised ATPase activity, implying a potential stationary role.

Myosin XVI Regulates Actin Cytoskeleton Dynamics in Dendritic Spines of Purkinje Cells and Affects Presynaptic Organization

The actin cytoskeleton is crucial for function and morphology of neuronal synapses. Moreover, altered regulation of the neuronal actin cytoskeleton has been implicated in neuropsychiatric diseases such as autism spectrum disorder (ASD). Myosin XVI is a neuronally expressed unconventional myosin known to bind the WAVE regulatory complex (WRC), a regulator of filamentous actin (F-actin) polymerization. Notably, the gene encoding the myosin’s heavy chain (MYO16) shows genetic association with neuropsychiatric disorders including ASD. Here, we investigated whether myosin XVI plays a role for actin cytoskeleton regulation in the dendritic spines of cerebellar Purkinje cells (PCs), a neuronal cell type crucial for motor learning, social cognition and vocalization. We provide evidence that both myosin XVI and the WRC component WAVE1 localize to PC spines. Fluorescence recovery after photobleaching (FRAP) analysis of GFP-actin in cultured PCs shows that Myo16 knockout as well as PC-specific Myo16 knockdown, lead to faster F-actin turnover in the dendritic spines of PCs. We also detect accelerated F-actin turnover upon interference with the WRC, and upon inhibition of Arp2/3 that drives formation of branched F-actin downstream of the WRC. In contrast, inhibition of formins that are responsible for polymerization of linear actin filaments does not cause faster F-actin turnover. Together, our data establish myosin XVI as a regulator of the postsynaptic actin cytoskeleton and suggest that it is an upstream activator of the WRC-Arp2/3 pathway in PC spines. Furthermore, ultra-structural and electrophysiological analyses of Myo16 knockout cerebellum reveals the presence of reduced numbers of synaptic vesicles at presynaptic terminals in the absence of the myosin. Therefore, we here define myosin XVI as an F-actin regulator important for presynaptic organization in the cerebellum.

Multivariate analysis of genome-wide data to identify potential pleiotropic genes for five major psychiatric disorders using MetaCCA

BACKGROUND

Genome-wide association studies have been extensively applied in identifying SNP associated with major psychiatric disorders. However, the SNPs identified by the prevailing univariate approach only explain a small percentage of the genetic variance of traits, and the extensive data have shown the major psychiatric disorders have common biological mechanisms and the overlapping pathophysiological pathways.

METHODS

We applied the genetic pleiotropy-informed metaCCA method on summary statistics data from the Psychiatric Genomics Consortium Cross-Disorder Group to examine the overlapping genetic relations between the five major psychiatric disorders. Furthermore, to refine all genes, we performed gene-based association analyses for the five disorders respectively using VEGAS2. Gene enrichment analysis was applied to explore the potential functional significance of the identified genes.

RESULTS

After metaCCA analysis, 1147 SNPs reached the Bonferroni corrected threshold (p < 1.06 × 10-6) in the univariate SNP-multivariate phenotype analysis, and 246 genes with a significance threshold (p < 3.85 × 10-6) were identified as potentially pleiotropic genes in the multivariate SNP-multivariate phenotype analysis. By screening the results of gene-based p-values, we identified 37 putative pleiotropic genes which achieved significance threshold in metaCCA analyses and were also associated with at least one disorder in the VEGAS2 analyses.

LIMITATIONS

Alternative approaches and experimental studies may be applied to check whether novel genes could still be identified/substantiated with these methods.

CONCLUSIONS

The metaCCA method identified novel variants associated with psychiatric disorders by effectively incorporating information from different GWAS datasets. Our analyses may provide insights for some common therapeutic approaches of these five major psychiatric disorders based on the pleiotropic genes and common mechanisms identified.

Specificity in Etiology of Subtypes of Bipolar Disorder: Evidence From a Swedish Population-Based Family Study

BACKGROUND

Uncertainty remains whether bipolar I disorder (BDI) and bipolar II disorder (BDII) differ etiologically. We used a population-based family sample to examine the etiological boundaries between BDI and BDII by assessing their familial aggregation/coaggregation and by assessing the coaggregation between them and schizophrenia, depression, attention-deficit/hyperactivity disorder, eating disorders, autism spectrum disorder, substance use disorders, anxiety disorders, and personality disorders.

METHODS

By linking Swedish national registers, we established a population-based cohort (N = 15,685,511) and identified relatives with different biological relationships. Odds ratios (ORs) were used to measure the relative risk of BDI and BDII in relatives of individuals diagnosed with BDI (n = 4309) and BDII (n = 4178). The heritability for BDI and BDII and the genetic correlation across psychiatric disorders were estimated by variance decomposition analysis.

RESULTS

Compared with the general population, the OR of BDI was 17.0 (95% confidence interval [CI] 13.1-22.0) in first-degree relatives of BDI patients, higher than that of BDII patients (OR 9.8, 95% CI 7.7-12.5). The ORs of BDII were 13.6 (95% CI 10.2-18.2) in first-degree relatives of BDII patients and 9.8 (95% CI 7.7-12.4) in relatives of BDI patients. The heritabilities for BDI and BDII were estimated at 57% (95% CI 32%-79%) and 46% (95% CI 21%-67%), respectively, with a genetic correlation estimated as 0.78 (95% CI 0.36-1.00). The familial coaggregation of other psychiatric disorders, in particular schizophrenia, showed different patterns for BDI and BDII.

CONCLUSIONS

Our results suggest a distinction between BDI and BDII in etiology, partly due to genetic differences.

Evaluation of the interaction between genetic variants of GAD1 and miRNA in bipolar disorders

BACKGROUND

Glutamic acid dehydrogenase 1 (GAD1) serves as the rate-limiting enzyme for synthesizing GABA, and is reported to be associated with several psychiatric disorders. The present study examined the effects of GAD1 genetic variants on bipolar disorder (BD) and its subtypes. Moreover, we investigated functional interactions between genetic variants and miRNAs via algorithm prediction and experimental validation.

METHODS

A case-control study was conducted with 280 BD patients and 200 healthy controls. Eight tag SNPs in GAD1 were genotyped. For associated markers, we performed in silico prediction for their potential functions through SNP-miRNA interactions by establishing a scoring system to combine information from several miRNA predictive algorithms. We then tested allelic expression differences using Dual-Glo luciferase reporter assays for the selected SNP-miRNA pair. Lastly, we examined the associations of the GAD1 gene and BD in two additional independent datasets with a few thousand samples for replication.

RESULTS

Marker rs3749034 was associated with BD, in particular the BD-II subtype. According to our scoring system, several candidate miRNAs were predicted to interact with rs3749034, and hsa-miR-504 had the highest score. Findings from an in vitro experiment revealed a non-statistically significant trend for lower gene expression level with the A allele of rs3749034 compared with the G allele. The association between rs3749034 and BD was not replicated in either of the independent datasets. Instead, other rarer genetic variants in GAD1 showed suggestive signals (e.g. rs575441409, p-value = 3.8*10^-4, D' = 1 with rs3749034) with BD in the Taiwanese dataset.

LIMITATIONS

The present study considered common genetic variants only. In addition, we only used a 293T cell-line in conducting luciferase reporter assays, as no primary cell-lines from patient samples were available to differentiate the effects between BD subtypes.

CONCLUSIONS

Our results demonstrate a weak effect of the GAD1 gene on the risk of bipolar illness, and the associated marker might represent a proxy for real signals of rare variants.

Genome-wide association study of borderline personality disorder reveals genetic overlap with bipolar disorder, major depression and schizophrenia

Borderline personality disorder (BOR) is determined by environmental and genetic factors, and characterized by affective instability and impulsivity, diagnostic symptoms also observed in manic phases of bipolar disorder (BIP). Up to 20% of BIP patients show comorbidity with BOR. This report describes the first case-control genome-wide association study (GWAS) of BOR, performed in one of the largest BOR patient samples worldwide. The focus of our analysis was (i) to detect genes and gene sets involved in BOR and (ii) to investigate the genetic overlap with BIP. As there is considerable genetic overlap between BIP, major depression (MDD) and schizophrenia (SCZ) and a high comorbidity of BOR and MDD, we also analyzed the genetic overlap of BOR with SCZ and MDD. GWAS, gene-based tests and gene-set analyses were performed in 998 BOR patients and 1545 controls. Linkage disequilibrium score regression was used to detect the genetic overlap between BOR and these disorders. Single marker analysis revealed no significant association after correction for multiple testing. Gene-based analysis yielded two significant genes: DPYD (P=4.42 × 10-7) and PKP4 (P=8.67 × 10-7); and gene-set analysis yielded a significant finding for exocytosis (GO:0006887, PFDR=0.019; FDR, false discovery rate). Prior studies have implicated DPYD, PKP4 and exocytosis in BIP and SCZ. The most notable finding of the present study was the genetic overlap of BOR with BIP (rg=0.28 [P=2.99 × 10-3]), SCZ (rg=0.34 [P=4.37 × 10-5]) and MDD (rg=0.57 [P=1.04 × 10-3]). We believe our study is the first to demonstrate that BOR overlaps with BIP, MDD and SCZ on the genetic level. Whether this is confined to transdiagnostic clinical symptoms should be examined in future studies.

Transcription factor 4 gene rs9960767 polymorphism in bipolar disorder

The transcription factor 4 (TCF4) gene encodes a helix-loop-helix transcription factor protein, which initiates neuronal differentiation and is primarily expressed during nervous system development. The aim of the present study is to investigate the association of the TCF4 rs9960767 polymorphism and bipolar disorder, which is highly heritable. DNA isolation was performed on 95 patients with bipolar disorder and 108 healthy control subjects to examine the TCF4 rs9960767 polymorphism. Genotypic and allelic frequencies were determined using the polymerase chain reaction-restriction fragment length polymorphism method designed in our laboratory. Statistical analysis was performed using χ² test within the 95% confidence interval. Odds ratios were calculated and Hardy-Weinberg equilibrium (HWE) was verified for all control subjects and patients. The A allele frequency was 95.8% in the patients and 94.4% in the control subjects, and 4.2% in the patients and 5.6% in the control subjects for the C allele. The genotype frequencies of the TCF4 gene rs9960767 variant were as follows: AA, 91.6% and AC, 8.4% in patients with bipolar (CC genotype was not observed in cases); AA, 89.8%; AC, 9.3% and CC, 0.9% in the control subjects. No statistically significant difference was identified between the patients and control subjects (χ²=0.937; P=0.626). In addition, gender specific analysis was performed, although no significant association was found according to the gender distrubition. All patients and control subjects were in HWE (P>0.05). Statistical analysis of the data indicates that the TCF4 gene rs9960767 polymorphism is not an independent risk factor for bipolar disorder in the overall population or in terms of gender; however, an increased population size would improve the statistical power. Furthermore, additional gene variants that are specifically involved in neuronal development may be analyzed for revealing the complex genetic architecture of bipolar disorder. An improved approach would be better to evaluate the TCF4 gene in a pathway specific manner due to its role as a transcription factor.