Rare variants implicate NMDA receptor signaling and cerebellar gene networks in risk for bipolar disorder

Bipolar disorder is an often-severe mental health condition characterized by alternation between extreme mood states of mania and depression. Despite strong heritability and the recent identification of 64 common variant risk loci of small effect, pathophysiological mechanisms remain unknown. Here, we analyzed genome sequences from 41 multiply-affected pedigrees and identified variants in 741 genes with nominally significant linkage or association with bipolar disorder. These 741 genes overlapped known risk genes for neurodevelopmental disorders and clustered within gene networks enriched for synaptic and nuclear functions. The top variant in this analysis - prioritized by statistical association, predicted deleteriousness, and network centrality - was a missense variant in the gene encoding D-amino acid oxidase (DAOG131V). Heterologous expression of DAOG131V in human cells resulted in decreased DAO protein abundance and enzymatic activity. In a knock-in mouse model of DAOG131, DaoG130V/+, we similarly found decreased DAO protein abundance in hindbrain regions, as well as enhanced stress susceptibility and blunted behavioral responses to pharmacological inhibition of N-methyl-D-aspartate receptors (NMDARs). RNA sequencing of cerebellar tissue revealed that DaoG130V resulted in decreased expression of two gene networks that are enriched for synaptic functions and for genes expressed, respectively, in Purkinje neurons or granule neurons. These gene networks were also down-regulated in the cerebellum of patients with bipolar disorder compared to healthy controls and were enriched for additional rare variants associated with bipolar disorder risk. These findings implicate dysregulation of NMDAR signaling and of gene expression in cerebellar neurons in bipolar disorder pathophysiology and provide insight into its genetic architecture.

Genome-Wide Linkage Study Meta-Analysis of Male Sexual Orientation

Male sexual orientation is a scientifically and socially important trait shown by family and twin studies to be influenced by environmental and complex genetic factors. Individual genome-wide linkage studies (GWLS) have been conducted, but not jointly analyzed. Two main datasets account for > 90% of the published GWLS concordant sibling pairs on the trait and are jointly analyzed here: MGSOSO (Molecular Genetic Study of Sexual Orientation; 409 concordant sibling pairs in 384 families, Sanders et al. (2015)) and Hamer (155 concordant sibling pairs in 145 families, Mustanski et al. (2005)). We conducted multipoint linkage analyses with Merlin on the datasets separately since they were genotyped differently, integrated genetic marker positions, and combined the resultant LOD (logarithm of the odds) scores at each 1 cM grid position. We continue to find the strongest linkage support at pericentromeric chromosome 8 and chromosome Xq28. We also incorporated the remaining published GWLS dataset (on 55 families) by using meta-analytic approaches on published summary statistics. The meta-analysis has maximized the positional information from GWLS of currently available family resources and can help prioritize findings from genome-wide association studies (GWAS) and other approaches. Although increasing evidence highlights genetic contributions to male sexual orientation, our current understanding of contributory loci is still limited, consistent with the complexity of the trait. Further increasing genetic knowledge about male sexual orientation, especially via large GWAS, should help advance our understanding of the biology of this important trait.

NRXN1 is associated with enlargement of the temporal horns of the lateral ventricles in psychosis

Schizophrenia, Schizoaffective, and Bipolar disorders share behavioral and phenomenological traits, intermediate phenotypes, and some associated genetic loci with pleiotropic effects. Volumetric abnormalities in brain structures are among the intermediate phenotypes consistently reported associated with these disorders. In order to examine the genetic underpinnings of these structural brain modifications, we performed genome-wide association analyses (GWAS) on 60 quantitative structural brain MRI phenotypes in a sample of 777 subjects (483 cases and 294 controls pooled together). Genotyping was performed with the Illumina PsychChip microarray, followed by imputation to the 1000 genomes multiethnic reference panel. Enlargement of the Temporal Horns of Lateral Ventricles (THLV) is associated with an intronic SNP of the gene NRXN1 (rs12467877, P = 6.76E-10), which accounts for 4.5% of the variance in size. Enlarged THLV is associated with psychosis in this sample, and with reduction of the hippocampus and enlargement of the choroid plexus and caudate. Eight other suggestively significant associations (P < 5.5E-8) were identified with THLV and 5 other brain structures. Although rare deletions of NRXN1 have been previously associated with psychosis, this is the first report of a common SNP variant of NRXN1 associated with enlargement of the THLV in psychosis.

Genome-wide association study identifies 30 loci associated with bipolar disorder

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Efficient region-based test strategy uncovers genetic risk factors for functional outcome in bipolar disorder

Genome-wide association studies of case-control status have advanced the understanding of the genetic basis of psychiatric disorders. Further progress may be gained by increasing sample size but also by new analysis strategies that advance the exploitation of existing data, especially for clinically important quantitative phenotypes. The functionally-informed efficient region-based test strategy (FIERS) introduced herein uses prior knowledge on biological function and dependence of genotypes within a powerful statistical framework with improved sensitivity and specificity for detecting consistent genetic effects across studies. As proof of concept, FIERS was used for the first genome-wide single nucleotide polymorphism (SNP)-based investigation on bipolar disorder (BD) that focuses on an important aspect of disease course, the functional outcome. FIERS identified a significantly associated locus on chromosome 15 (hg38: chr15:48965004 - 49464789 bp) with consistent effect strength between two independent studies (GAIN/TGen: European Americans, BOMA: Germans; n = 1592 BD patients in total). Protective and risk haplotypes were found on the most strongly associated SNPs. They contain a CTCF binding site (rs586758); CTCF sites are known to regulate sets of genes within a chromatin domain. The rs586758 - rs2086256 - rs1904317 haplotype is located in the promoter flanking region of the COPS2 gene, close to microRNA4716, and the EID1, SHC4, DTWD1 genes as plausible biological candidates. While implication with BD is novel, COPS2, EID1, and SHC4 are known to be relevant for neuronal differentiation and function and DTWD1 for psychopharmacological side effects. The test strategy FIERS that enabled this discovery is equally applicable for tag SNPs and sequence data.

The transcription factor POU3F2 regulates a gene coexpression network in brain tissue from patients with psychiatric disorders

Schizophrenia and bipolar disorder are complex psychiatric diseases with risks contributed by multiple genes. Dysregulation of gene expression has been implicated in these disorders, but little is known about such dysregulation in the human brain. We analyzed three transcriptome datasets from 394 postmortem brain tissue samples from patients with schizophrenia or bipolar disorder or from healthy control individuals without a known history of psychiatric disease. We built genome-wide coexpression networks that included microRNAs (miRNAs). We identified a coexpression network module that was differentially expressed in the brain tissue from patients compared to healthy control individuals. This module contained genes that were principally involved in glial and neural cell genesis and glial cell differentiation, and included schizophrenia risk genes carrying rare variants. This module included five miRNAs and 545 mRNAs, with six transcription factors serving as hub genes in this module. We found that the most connected transcription factor gene POU3F2, also identified on a genome-wide association study for bipolar disorder, could regulate the miRNA hsa-miR-320e and other putative target mRNAs. These regulatory relationships were replicated using PsychENCODE/BrainGVEX datasets and validated by knockdown and overexpression experiments in SH-SY5Y cells and human neural progenitor cells in vitro. Thus, we identified a brain gene expression module that was enriched for rare coding variants in genes associated with schizophrenia and that contained the putative bipolar disorder risk gene POU3F2 The transcription factor POU3F2 may be a key regulator of gene expression in this disease-associated gene coexpression module.

Genome-Wide Association Study of Male Sexual Orientation

Family and twin studies suggest that genes play a role in male sexual orientation. We conducted a genome-wide association study (GWAS) of male sexual orientation on a primarily European ancestry sample of 1,077 homosexual men and 1,231 heterosexual men using Affymetrix single nucleotide polymorphism (SNP) arrays. We identified several SNPs with p < 10-5, including regions of multiple supporting SNPs on chromosomes 13 (minimum p = 7.5 × 10-7) and 14 (p = 4.7 × 10-7). The genes nearest to these peaks have functions plausibly relevant to the development of sexual orientation. On chromosome 13, SLITRK6 is a neurodevelopmental gene mostly expressed in the diencephalon, which contains a region previously reported as differing in size in men by sexual orientation. On chromosome 14, TSHR genetic variants in intron 1 could conceivably help explain past findings relating familial atypical thyroid function and male homosexuality. Furthermore, skewed X chromosome inactivation has been found in the thyroid condition, Graves' disease, as well as in mothers of homosexual men. On pericentromeric chromosome 8 within our previously reported linkage peak, we found support (p = 4.1 × 10-3) for a SNP association previously reported (rs77013977, p = 7.1 × 10-8), with the combined analysis yielding p = 6.7 × 10-9, i.e., a genome-wide significant association.

A Rare Variant in CACNA1D Segregates with 7 Bipolar I Disorder Cases in a Large Pedigree

Whole-genome sequencing was performed on 3 bipolar I disorder (BPI) cases from a multiplex pedigree of European ancestry with 7 BPI cases. Within CACNA1D, a gene implicated by genome-wide association studies, a G to C nucleotide transversion at 53,835,340 base pairs (bps) was found predicting the substitution of proline for alanine at amino acid position 1751 (A1751P). Using Sanger sequencing, the DNA variant was shown to co-segregate with the remaining 4 BPI cases within the pedigree. A high-resolution DNA denaturing curve method was then used to screen for the presence of the A1751P change in 4,150 BPI cases from the NIMH Genetics Initiative. The A1751P variant was found in 4 BPI cases. A second variant within exon 43, a C to T nucleotide transition, was found in 1 case at 53,835,355 bps, predicting the substitution of tryptophan for arginine at amino acid position 1771 (R1771W). In the NHLBI Exome Sequencing Project database, the heterozygous A1751P variant was present in 3 of 4,300 subjects of European ancestry, and the R1771W change was not present in any subject. Given the rarity of these variants, large-scale case/control rare variant sequencing studies will be required for definitive conclusions.

Genome-wide association study of 40,000 individuals identifies two novel loci associated with bipolar disorder

Bipolar disorder (BD) is a genetically complex mental illness characterized by severe oscillations of mood and behaviour. Genome-wide association studies (GWAS) have identified several risk loci that together account for a small portion of the heritability. To identify additional risk loci, we performed a two-stage meta-analysis of >9 million genetic variants in 9,784 bipolar disorder patients and 30,471 controls, the largest GWAS of BD to date. In this study, to increase power we used ∼2,000 lithium-treated cases with a long-term diagnosis of BD from the Consortium on Lithium Genetics, excess controls, and analytic methods optimized for markers on the X-chromosome. In addition to four known loci, results revealed genome-wide significant associations at two novel loci: an intergenic region on 9p21.3 (rs12553324, P =  5.87 × 10 ^{-  9}; odds ratio (OR) = 1.12) and markers within ERBB2 (rs2517959, P =  4.53 × 10 ^{-  9}; OR = 1.13). No significant X-chromosome associations were detected and X-linked markers explained very little BD heritability. The results add to a growing list of common autosomal variants involved in BD and illustrate the power of comparing well-characterized cases to an excess of controls in GWAS.

Pursuit eye movements as an intermediate phenotype across psychotic disorders: Evidence from the B-SNIP study

Smooth pursuit eye tracking deficits are a promising intermediate phenotype for schizophrenia and possibly for psychotic disorders more broadly. The Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP) consortium investigated the severity and familiality of different pursuit parameters across psychotic disorders. Probands with schizophrenia (N=265), schizoaffective disorder (N=178), psychotic bipolar disorder (N=231), their first-degree relatives (N=306, N=217, N=273, respectively) and healthy controls (N=305) performed pursuit tracking tasks designed to evaluate sensorimotor and cognitive/predictive aspects of pursuit. Probands from all diagnostic groups were impaired on all pursuit measures of interest compared to controls (p<0.001). Schizophrenia probands were more impaired than other proband groups on both early pursuit gain and predictive gain. Relatives with and without enhanced psychosis spectrum personality traits were impaired on initial eye acceleration, the most direct sensorimotor pursuit measure, but not on pursuit gain measures. This suggests that alterations in early sensorimotor function may track susceptibility to psychosis even in the absence of psychosis related personality traits. There were no differences in pursuit measures between relatives of the three proband groups. Familiality estimates of pursuit deficits indicate that early pursuit gain was more familial than predictive gain, which has been the most widely used measure in previous family studies of psychotic disorders. Thus, while disease-related factors may induce significant impairments of pursuit gain, especially in schizophrenia, the pattern of deficits in relatives and their familiality estimates suggest that alterations in sensorimotor function at pursuit onset may indicate increased susceptibility across psychotic disorders.

Neuropsychological and dimensional behavioral trait profiles in Costa Rican ADHD sib pairs: Potential intermediate phenotypes for genetic studies

Attention deficit hyperactivity disorder (ADHD) is associated with substantial functional impairment in children and in adults. Many individuals with ADHD have clear neurocognitive deficits, including problems with visual attention, processing speed, and set shifting. ADHD is etiologically complex, and although genetic factors play a role in its development, much of the genetic contribution to ADHD remains unidentified. We conducted clinical and neuropsychological assessments of 294 individuals (269 with ADHD) from 163 families (48 multigenerational families created using genealogical reconstruction, 78 affected sib pair families, and 37 trios) from the Central Valley of Costa Rica (CVCR). We used principal components analysis (PCA) to group neurocognitive and behavioral variables using the subscales of the Child Behavior Checklist (CBCL) and 15 neuropsychological measures, and created quantitative traits for heritability analyses. We identified seven cognitive and two behavioral domains. Individuals with ADHD were significantly more impaired than their unaffected siblings on most behavioral and cognitive domains. The verbal IQ domain had the highest heritability (92%), followed by auditory attention (87%), visual processing speed and problem solving (85%), and externalizing symptoms (81%). The quantitative traits identified here have high heritabilities, similar to the reported heritability of ADHD (70-90%), and may represent appropriate alternative phenotypes for genetic studies. The use of multigenerational families from a genetically isolated population may facilitate the identification of ADHD risk genes in the face of phenotypic and genetic heterogeneity.

Genome-wide scan demonstrates significant linkage for male sexual orientation

BACKGROUND

Findings from family and twin studies support a genetic contribution to the development of sexual orientation in men. However, previous studies have yielded conflicting evidence for linkage to chromosome Xq28.

METHOD

We conducted a genome-wide linkage scan on 409 independent pairs of homosexual brothers (908 analyzed individuals in 384 families), by far the largest study of its kind to date.

RESULTS

We identified two regions of linkage: the pericentromeric region on chromosome 8 (maximum two-point LOD = 4.08, maximum multipoint LOD = 2.59), which overlaps with the second strongest region from a previous separate linkage scan of 155 brother pairs; and Xq28 (maximum two-point LOD = 2.99, maximum multipoint LOD = 2.76), which was also implicated in prior research.

CONCLUSIONS

Results, especially in the context of past studies, support the existence of genes on pericentromeric chromosome 8 and chromosome Xq28 influencing development of male sexual orientation.

Genome-wide association study implicates NDST3 in schizophrenia and bipolar disorder

Schizophrenia and bipolar disorder are major psychiatric disorders with high heritability and overlapping genetic variance. Here we perform a genome-wide association study in an ethnically homogeneous cohort of 904 schizophrenia cases and 1,640 controls drawn from the Ashkenazi Jewish population. We identify a novel genome-wide significant risk locus at chromosome 4q26, demonstrating the potential advantages of this founder population for gene discovery. The top single-nucleotide polymorphism (SNP; rs11098403) demonstrates consistent effects across 11 replication and extension cohorts, totalling 23, 191 samples across multiple ethnicities, regardless of diagnosis (schizophrenia or bipolar disorder), resulting in P_meta=9.49 × 10⁻¹² (odds ratio (OR)=1.13, 95% confidence interval (CI): 1.08–1.17) across both disorders and P_meta=2.67 × 10⁻⁸ (OR=1.15, 95% CI: 1.08–1.21) for schizophrenia alone. In addition, this intergenic SNP significantly predicts postmortem cerebellar gene expression of NDST3, which encodes an enzyme critical to heparan sulphate metabolism. Heparan sulphate binding is critical to neurite outgrowth, axon formation and synaptic processes thought to be aberrant in these disorders.

Schizophrenia and bipolar disorder are important psychiatric disorders with overlapping genetic components. Here, the authors identify and replicate a genome-wide significant risk locus for the two disorders, and suggest a role for NDST3 in severe psychiatric disease.

Two gene co-expression modules differentiate psychotics and controls

Schizophrenia (SCZ) and bipolar disorder (BD) are highly heritable psychiatric disorders. Associated genetic and gene expression changes have been identified, but many have not been replicated and have unknown functions. We identified groups of genes whose expressions varied together, that is co-expression modules, then tested them for association with SCZ. Using weighted gene co-expression network analysis, we show that two modules were differentially expressed in patients versus controls. One, upregulated in cerebral cortex, was enriched with neuron differentiation and neuron development genes, as well as disease genome-wide association study genetic signals; the second, altered in cerebral cortex and cerebellum, was enriched with genes involved in neuron protection functions. The findings were preserved in five expression data sets, including sets from three brain regions, from a different microarray platform, and from BD patients. From those observations, we propose neuron differentiation and development pathways may be involved in etiologies of both SCZ and BD, and neuron protection function participates in pathological process of the diseases.

Heritability and linkage analysis of personality in bipolar disorder

BACKGROUND

The many attempts that have been made to identify genes for bipolar disorder (BD) have met with limited success, which may reflect an inadequacy of diagnosis as an informative and biologically relevant phenotype for genetic studies. Here we have explored aspects of personality as quantitative phenotypes for bipolar disorder through the use of the Temperament and Character Inventory (TCI), which assesses personality in seven dimensions. Four temperament dimensions are assessed: novelty seeking (NS), harm avoidance (HA), reward dependence (RD), and persistence (PS). Three character dimensions are also included: self-directedness (SD), cooperativeness (CO), and self-transcendence (ST).

METHODS

We compared personality scores between diagnostic groups and assessed heritability in a sample of 101 families collected for genetic studies of BD. A genome-wide SNP linkage analysis was then performed in the subset of 51 families for which genetic data was available.

RESULTS

Significant group differences were observed between BD subjects, their first-degree relatives, and independent controls for all but RD and PS, and all but HA and RD were found to be significantly heritable in this sample. Linkage analysis of the heritable dimensions produced several suggestive linkage peaks for NS (chromosomes 7q21 and 10p15), PS (chromosomes 6q16, 12p13, and 19p13), and SD (chromosomes 4q35, 8q24, and 18q12).

LIMITATIONS

The relatively small size of our linkage sample likely limited our ability to reach genome-wide significance in this study.

CONCLUSIONS

While not genome-wide significant, these results suggest that aspects of personality may prove useful in the identification of genes underlying BD susceptibility.

Heritability and genome-wide SNP linkage analysis of temperament in bipolar disorder

BACKGROUND

The many attempts to identify genes for bipolar disorder (BD) have met with limited success, which has generally been attributed to genetic heterogeneity and small gene effects. However, it is also possible that the categorical phenotypes used in genetic studies of BD are not the most informative or biologically relevant. We have explored aspects of temperament as quantitative phenotypes for BD through the use of the Temperament Evaluation of Memphis, Pisa, Paris, and San Diego Auto-questionnaire (TEMPS-A), which is designed to assess lifelong, milder aspects of bipolar symptomatology and defines five temperaments: hyperthymic, dysthymic, cyclothymic, irritable, and anxious.

METHODS

We compared temperament scores between diagnostic groups and assessed heritability in a sample of 101 families collected for genetic studies of BD. A genome-wide SNP linkage study was then performed in the subset of 51 families for which genetic data was available.

RESULTS

Significant group differences were observed between BD subjects, their first-degree relatives, and independent controls, and all five temperaments were found to be significantly heritable, with heritabilities ranging from 21% for the hyperthymic to 52% for the irritable temperaments. Suggestive evidence for linkage was observed for the hyperthymic (chromosomes 1q44, 2p16, 6q16, and 14q23), dysthymic (chromosomes 3p21 and 13q34), and irritable (chromosome 6q24) temperaments.

LIMITATIONS

The relatively small size of our linkage sample likely limited our ability to reach genome-wide significance in this study.

CONCLUSIONS

While not genome-wide significant, these results suggest that aspects of temperament may prove useful in the identification of genes underlying BD susceptibility.

Enrichment of cis-regulatory gene expression SNPs and methylation quantitative trait loci among bipolar disorder susceptibility variants

We conducted a systematic study of top susceptibility variants from a genome-wide association (GWA) study of bipolar disorder to gain insight into the functional consequences of genetic variation influencing disease risk. We report here the results of experiments to explore the effects of these susceptibility variants on DNA methylation and mRNA expression in human cerebellum samples. Among the top susceptibility variants, we identified an enrichment of cis regulatory loci on mRNA expression (eQTLs), and a significant excess of quantitative trait loci for DNA CpG methylation, hereafter referred to as methylation quantitative trait loci (mQTLs). Bipolar disorder susceptibility variants that cis regulate both cerebellar expression and methylation of the same gene are a very small proportion of bipolar disorder susceptibility variants. This finding suggests that mQTLs and eQTLs provide orthogonal ways of functionally annotating genetic variation within the context of studies of pathophysiology in brain. No lymphocyte mQTL enrichment was found, suggesting that mQTL enrichment was specific to the cerebellum, in contrast to eQTLs. Separately, we found that using mQTL information to restrict the number of single-nucleotide polymorphisms studied enhances our ability to detect a significant association. With this restriction a priori informed by the observed functional enrichment, we identified a significant association (rs12618769, P(bonferroni)<0.05) from two other GWA studies (TGen+GAIN; 2191 cases and 1434 controls) of bipolar disorder, which we replicated in an independent GWA study (WTCCC). Collectively, our findings highlight the importance of integrating functional annotation of genetic variants for gene expression and DNA methylation to advance the biological understanding of bipolar disorder.

Evidence for association of bipolar disorder to haplotypes in the 22q12.3 region near the genes stargazin, IFT27 and parvalbumin

We have previously reported genome-wide significant linkage of bipolar disorder to a region on 22q12.3 near the marker D22S278. Towards identifying the susceptibility gene, we have conducted a fine-mapping association study of the region in two independent family samples, an independent case-control sample and a genome-wide association dataset. Two hundred SNPs were first examined in a 5 Mb region surrounding the D22S278 marker in a sample of 169 families and analyzed using PLINK. The peak of association was a haplotype near the genes stargazin (CACNG2), intraflagellar transport protein homolog 27 (IFT27) and parvalbumin (PVALB; P = 4.69 × 10(-4)). This peak overlapped a significant haplotype in a family based association study of a second independent sample of 294 families (P = 1.42 × 10(-5)). Analysis of the combined family sample yielded statistically significant evidence of association to a rare three SNP haplotype in the gene IFT27 (P = 8.89 × 10(-6)). Twelve SNPs comprising these haplotypes were genotyped in an independent sample of 574 bipolar I cases and 550 controls. Statistically significant association was found for a haplotype window that overlapped the region from the first two family samples (P = 3.43 × 10(-4)). However, analyses of the two family samples using the program LAMP, found no evidence for association in this region, but did yield significant evidence for association to a haplotype 3' of CACNG2 (P = 1.76 × 10(-6)). Furthermore, no evidence for association was found in a large genome-wide association dataset. The replication of association to overlapping haplotypes in three independent datasets suggests the presence of a bipolar disorder susceptibility gene in this region.

Genome-wide linkage analysis of obsessive-compulsive disorder implicates chromosome 1p36

BACKGROUND

Obsessive-compulsive disorder (OCD) has a complex etiology involving both genetic and environmental factors. However, the genetic causes of OCD are largely unknown, despite the identification of several promising candidate genes and linkage regions.

METHODS

Our objective was to conduct genetic linkage studies of the type of OCD thought to have the strongest genetic etiology (i.e., childhood-onset OCD), in 33 Caucasian families with ≥2 childhood-onset OCD-affected individuals from the United States (n = 245 individuals with genotype data). Parametric and nonparametric genome-wide linkage analyses were conducted with Morgan and Merlin in these families using a selected panel of single nucleotide repeat polymorphisms from the Illumina 610-Quad Bead Chip. The initial analyses were followed by fine-mapping analyses in genomic regions with initial heterogeneity logarithm of odds (HLOD) scores of ≥2.0.

RESULTS

We identified five areas of interest (HLOD score ≥2) on chromosomes 1p36, 2p14, 5q13, 6p25, and 10p13. The strongest result was on chromosome 1p36.33-p36.32 (HLOD = 3.77, suggestive evidence for linkage after fine mapping). At this location, several of the families showed haplotypes co-segregating with OCD.

CONCLUSIONS

The results of this study represent the strongest linkage finding for OCD in a primary analysis to date and suggest that chromosome 1p36, and possibly several other genomic regions, may harbor susceptibility loci for OCD. Multiple brain-expressed genes lie under the primary linkage peak (approximately 4 megabases in size). Follow-up studies, including replication in additional samples and targeted sequencing of the areas of interest, are needed to confirm these findings and to identify specific OCD risk variants.

Incorporation of molecular data and redefinition of phenotype: new approaches to genetic epidemiology of bipolar manic depressive illness and schizophrenia

Considerable advances have been made in identifying specific genetic components of bipolar manic depressive illness (BP) and schizophrenia (SZ), despite their complex inheritance. Meta-analysis of all published whole-genome linkage scans reveals overall support for illness genes in several chromosomal regions. In two of these regions, on the lonq arm of chromosome 13 and on the long arm of chromosome 22, the combined studies of BP and SZ are consistent with a common susceptibility locus for the two disorders. This lends some plausibility to the hypothesis of some shared genetic predispositions for BP and SZ. Other linkages are supported by multiple studies of specific chromosomal regions, most notably two regions on chromosome 6 in SZ. The velocardiofacial syndrome is associated with deletions very close to the linkage region on chromosome 22, and with psychiatric manifestations of both BP and SZ. Endophenotypes of SZ, previously demonstrated to be heritable, have been found to have chromosomal linkage in at least one study. These include eye-tracking abnormalities linked to the short arm of chromosome 6, and abnormality of the P50 cortical evoked potential linked to chromosome 15. Variants in specific genes have been associated with susceptibility to illness, and other genes have been associated with susceptibility to side effects of pharmacological treatment. These genetic findings may eventually be part of an integrated genetic, environmental, and interactive-factor epidemiology of the major mental illnesses.

Loci nominally associated with autism from genome-wide analysis show enrichment of brain expression quantitative trait loci but not lymphoblastoid cell line expression quantitative trait loci

Background

Autism spectrum disorder is a severe early onset neurodevelopmental disorder with high heritability but significant heterogeneity. Traditional genome-wide approaches to test for an association of common variants with autism susceptibility risk have met with limited success. However, novel methods to identify moderate risk alleles in attainable sample sizes are now gaining momentum.

Methods

In this study, we utilized publically available genome-wide association study data from the Autism Genome Project and annotated the results (P <0.001) for expression quantitative trait loci present in the parietal lobe (GSE35977), cerebellum (GSE35974) and lymphoblastoid cell lines (GSE7761). We then performed a test of enrichment by comparing these results to simulated data conditioned on minor allele frequency to generate an empirical P-value indicating statistically significant enrichment of expression quantitative trait loci in top results from the autism genome-wide association study.

Results

Our findings show a global enrichment of brain expression quantitative trait loci, but not lymphoblastoid cell line expression quantitative trait loci, among top single nucleotide polymorphisms from an autism genome-wide association study. Additionally, the data implicates individual genes SLC25A12, PANX1 and PANX2 as well as pathways previously implicated in autism.

Conclusions

These findings provide supportive rationale for the use of annotation-based approaches to genome-wide association studies.

A quantitative association study of SLC25A12 and restricted repetitive behavior traits in autism spectrum disorders

Background

SLC25A12 was previously identified by a linkage-directed association analysis in autism. In this study, we investigated the relationship between three SLC25A12 single nucleotide polymorphisms (SNPs) (rs2056202, rs908670 and rs2292813) and restricted repetitive behavior (RRB) traits in autism spectrum disorders (ASDs), based on a positive correlation between the G allele of rs2056202 and an RRB subdomain score on the Autism Diagnostic Interview-Revised (ADI-R).

Methods

We used the Repetitive Behavior Scale-Revised (RBS-R) as a quantitative RRB measure, and conducted linear regression analyses for individual SNPs and a previously identified haplotype (rs2056202-rs2292813). We examined associations in our University of Illinois at Chicago-University of Florida (UIC-UF) sample (179 unrelated individuals with an ASD), and then attempted to replicate our findings in the Simons Simplex Collection (SSC) sample (720 ASD families).

Results

In the UIC-UF sample, three RBS-R scores (ritualistic, sameness, sum) had positive associations with the A allele of rs2292813 (p = 0.006-0.012) and with the rs2056202-rs2292813 haplotype (omnibus test, p = 0.025-0.040). The SSC sample had positive associations between the A allele of rs2056202 and four RBS-R scores (stereotyped, sameness, restricted, sum) (p = 0.006-0.010), between the A allele of rs908670 and three RBS-R scores (stereotyped, self-injurious, sum) (p = 0.003-0.015), and between the rs2056202-rs2292813 haplotype and six RBS-R scores (stereotyped, self-injurious, compulsive, sameness, restricted, sum)(omnibus test, p = 0.002-0.028). Taken together, the A alleles of rs2056202 and rs2292813 were consistently and positively associated with RRB traits in both the UIC-UF and SSC samples, but the most significant SNP with phenotype association varied in each dataset.

Conclusions

This study confirmed an association between SLC25A12 and RRB traits in ASDs, but the direction of the association was different from that in the initial study. This could be due to the examined SLC25A12 SNPs being in linkage disequilibrium with another risk allele, and/or genetic/phenotypic heterogeneity of the ASD samples across studies.

Genetic control of individual differences in gene-specific methylation in human brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island "shores" and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations-the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.

Common and rare variants of DAOA in bipolar disorder

The D-amino acid oxidase activator (DAOA, previously known as G72) gene, mapped on 13q33, has been reported to be genetically associated with bipolar disorder (BP) in several populations. The consistency of associated variants is unclear and rare variants in exons of the DAOA gene have not been investigated in psychiatric diseases. We employed a conditional linkage method-STatistical Explanation for Positional Cloning (STEPC) to evaluate whether any associated single nucleotide polymorphisms (SNPs) account for the evidence of linkage in a pedigree series that previously has been linked to marker D13S779 at 13q33. We also performed an association study in a sample of 376 Caucasian BP parent-proband trios by genotyping 38 common SNPs in the gene region. Besides, we resequenced coding regions and flanking intronic sequences of DAOA in 555 Caucasian unrelated BP patients and 564 mentally healthy controls, to identify putative functional rare variants that may contribute to disease. One SNP rs1935058 could "explain" the linkage signal in the family sample set (P = 0.055) using STEPC analysis. No significant allelic association was detected in an association study by genotyping 38 common SNPs in 376 Caucasian BP trios. Resequencing identified 53 SNPs, of which 46 were novel SNPs. There was no significant excess of rare variants in cases relative to controls. Our results suggest that DAOA does not have a major effect on BP susceptibility. However, DAOA may contribute to bipolar susceptibility in some specific families as evidenced by the STEPC analysis.

Genome-wide association study of bipolar disorder in European American and African American individuals

To identify bipolar disorder (BD) genetic susceptibility factors, we conducted two genome-wide association (GWA) studies: one involving a sample of individuals of European ancestry (EA; n=1001 cases; n=1033 controls), and one involving a sample of individuals of African ancestry (AA; n=345 cases; n=670 controls). For the EA sample, single-nucleotide polymorphisms (SNPs) with the strongest statistical evidence for association included rs5907577 in an intergenic region at Xq27.1 (P=1.6 x 10(-6)) and rs10193871 in NAP5 at 2q21.2 (P=9.8 x 10(-6)). For the AA sample, SNPs with the strongest statistical evidence for association included rs2111504 in DPY19L3 at 19q13.11 (P=1.5 x 10(-6)) and rs2769605 in NTRK2 at 9q21.33 (P=4.5 x 10(-5)). We also investigated whether we could provide support for three regions previously associated with BD, and we showed that the ANK3 region replicates in our sample, along with some support for C15Orf53; other evidence implicates BD candidate genes such as SLITRK2. We also tested the hypothesis that BD susceptibility variants exhibit genetic background-dependent effects. SNPs with the strongest statistical evidence for genetic background effects included rs11208285 in ROR1 at 1p31.3 (P=1.4 x 10(-6)), rs4657247 in RGS5 at 1q23.3 (P=4.1 x 10(-6)), and rs7078071 in BTBD16 at 10q26.13 (P=4.5 x 10(-6)). This study is the first to conduct GWA of BD in individuals of AA and suggests that genetic variations that contribute to BD may vary as a function of ancestry.

Association and mutation analyses of 16p11.2 autism candidate genes

Background

Autism is a complex childhood neurodevelopmental disorder with a strong genetic basis. Microdeletion or duplication of a ∼500–700-kb genomic rearrangement on 16p11.2 that contains 24 genes represents the second most frequent chromosomal disorder associated with autism. The role of common and rare 16p11.2 sequence variants in autism etiology is unknown.

Methodology/Principal Findings

To identify common 16p11.2 variants with a potential role in autism, we performed association studies using existing data generated from three microarray platforms: Affymetrix 5.0 (777 families), Illumina 550 K (943 families), and Affymetrix 500 K (60 families). No common variants were identified that were significantly associated with autism. To look for rare variants, we performed resequencing of coding and promoter regions for eight candidate genes selected based on their known expression patterns and functions. In total, we identified 26 novel variants in autism: 13 exonic (nine non-synonymous, three synonymous, and one untranslated region) and 13 promoter variants. We found a significant association between autism and a coding variant in the seizure-related gene SEZ6L2 (12/1106 autism vs. 3/1161 controls; p = 0.018). Sez6l2 expression in mouse embryos was restricted to the spinal cord and brain. SEZ6L2 expression in human fetal brain was highest in post-mitotic cortical layers, hippocampus, amygdala, and thalamus. Association analysis of SEZ6L2 in an independent sample set failed to replicate our initial findings.

Conclusions/Significance

We have identified sequence variation in at least one candidate gene in 16p11.2 that may represent a novel genetic risk factor for autism. However, further studies are required to substantiate these preliminary findings.

Identification of genomic regions contributing to etoposide-induced cytotoxicity

Etoposide is routinely used in combination-based chemotherapy for testicular cancer and small-cell lung cancer; however, myelosuppression, therapy-related leukemia and neurotoxicity limit its utility. To determine the genetic contribution to cellular sensitivity to etoposide, we evaluated cell growth inhibition in Centre d' Etude du Polymorphisme Humain lymphoblastoid cell lines from 24 multi-generational pedigrees (321 samples) following treatment with 0.02-2.5 microM etoposide for 72 h. Heritability analysis showed that genetic variation contributes significantly to the cytotoxic phenotypes (h (2) = 0.17-0.25, P = 4.9 x 10(-5)-7.3 x 10(-3)). Whole genome linkage scans uncovered 8 regions with peak LOD scores ranging from 1.57 to 2.55, with the most significant signals being found on chromosome 5 (LOD = 2.55) and chromosome 6 (LOD = 2.52). Linkage-directed association was performed on a subset of HapMap samples within the pedigrees to find 22 SNPs significantly associated with etoposide cytotoxicity at one or more treatment concentrations. UVRAG, a DNA repair gene, SEMA5A, SLC7A6 and PRMT7 are implicated from these unbiased studies. Our findings suggest that susceptibility to etoposide-induced cytotoxicity is heritable and using an integrated genomics approach we identified both genomic regions and SNPs associated with the cytotoxic phenotypes.

Clock genes may influence bipolar disorder susceptibility and dysfunctional circadian rhythm

Several previous studies suggest that dysfunction of circadian rhythms may increase susceptibility to bipolar disorder (BP). We conducted an association study of five circadian genes (CRY2, PER1-3, and TIMELESS) in a family collection of 36 trios and 79 quads (Sample I), and 10 circadian genes (ARNTL, ARNTL2, BHLHB2, BHLHB3, CLOCK, CRY1, CSNK1D, CSNK1E, DBP, and NR1D1) in an extended family collection of 70 trios and 237 quads (Sample II), which includes the same 114 families but not necessarily the same individuals as Sample I. In Sample II, the Sibling-Transmission Disequilibrium Test (sib-tdt) analysis showed nominally significant association of BP with three SNPs within or near the CLOCK gene (rs534654, P = 0.0097; rs6850524, P = 0.012; rs4340844, P = 0.015). In addition, SNPs in the ARNTL2, CLOCK, DBP, and TIMELESS genes and haplotypes in the ARNTL, CLOCK, CSNK1E, and TIMELESS genes showed suggestive evidence of association with several circadian phenotypes identified in BP patients. However, none of these associations reached gene-wide or experiment-wide significance after correction for multiple-testing. A multi-locus interaction between rs6442925 in the 5' upstream of BHLHB2, rs1534891 in CSNK1E, and rs534654 near the 3' end of the CLOCK gene, however, is significantly associated with BP (P = 0.00000172). It remains significant after correcting for multiple testing using the False Discovery Rate method. Our results indicate an interaction between three circadian genes in susceptibility to bipolar disorder.

Neurotransmission and bipolar disorder: a systematic family-based association study

Neurotransmission pathways/systems have been proposed to be involved in the pathophysiology and treatment of bipolar disorder for over 40 years. In order to test the hypothesis that common variants of genes in one or more of five neurotransmission systems confer risk for bipolar disorder, we analyzed 1,005 tag single nucleotide polymorphisms in 90 genes from dopaminergic, serotonergic, noradrenergic, GABAergic, and glutamatergic neurotransmitter systems in 101 trios and 203 quads from Caucasian bipolar families. Our sample has 80% power to detect ORs >or= 1.82 and >or=1.57 for minor allele frequencies of 0.1 and 0.5, respectively. Nominally significant allelic and haplotypic associations were found for genes from each neurotransmission system, with several reaching gene-wide significance (allelic: GRIA1, GRIN2D, and QDPR; haplotypic: GRIN2C, QDPR, and SLC6A3). However, none of these associations survived correction for multiple testing in an individual system, or in all systems considered together. Significant single nucleotide polymorphism associations were not found with sub-phenotypes (alcoholism, psychosis, substance abuse, and suicide attempts) or significant gene-gene interactions. These results suggest that, within the detectable odds ratios of this study, common variants of the selected genes in the five neurotransmission systems do not play major roles in influencing the risk for bipolar disorder or comorbid sub-phenotypes.

Association study of Wnt signaling pathway genes in bipolar disorder

CONTEXT

The Wnt signaling pathways promote cell growth and are best known for their role in embryogenesis and cancer. Several lines of evidence suggest that these pathways might also be involved in bipolar disorder.

OBJECTIVE

To test for an association between candidate genes in the Wnt signaling pathways and disease susceptibility in a family-based bipolar disorder study.

DESIGN

Two hundred twenty-seven tagging single- nucleotide polymorphisms (SNPs) from 34 genes were successfully genotyped. Initial results led us to focus on the gene PPARD, in which we genotyped an additional 13 SNPs for follow-up.

SETTING

Nine academic medical centers in the United States.

PARTICIPANTS

Five hundred fifty-four offspring with bipolar disorder and their parents from 317 families.

MAIN OUTCOME MEASURES

Family-based association using FBAT and HBAT (http://www.biostat.harvard.edu/~fbat/default.html; Harvard School of Public Health, Boston, Massachusetts). Exploratory analyses testing for interactions of PPARD SNPs with clinical covariates and with other Wnt genes were conducted with GENASSOC (Stata Corp, College Station, Texas).

RESULTS

In the initial analysis, the most significantly associated SNP was rs2267665 in PPARD (nominal P < .001). This remained significant at P = .05 by permutation after accounting for all SNPs tested. Additional genotyping in PPARD yielded 4 SNPs in 1 haplotype block that were significantly associated with bipolar disorder (P < .01), the most significant being rs9462082 (P < .001). Exploratory analyses revealed significant evidence (P < .01) for interactions of rs9462082 with poor functioning on the Global Assessment Scale (odds ratio [OR], 3.36; 95% confidence interval [CI], 1.85-6.08) and with SNPs in WNT2B (rs3790606: OR, 2.56; 95% CI, 1.67-4.00) and WNT7A (rs4685048: OR, 1.79; 95% CI, 1.23-2.63).

CONCLUSIONS

We found evidence for association of bipolar disorder with PPARD, a gene in the Wnt signaling pathway. The consistency of this result with one from the Wellcome Trust Case-Control Consortium encourages further study. If the finding can be confirmed in additional samples, it may illuminate a new avenue for understanding the pathogenesis of severe bipolar disorder and developing more effective treatments.

Novel submicroscopic chromosomal abnormalities detected in autism spectrum disorder

BACKGROUND

One genetic mechanism known to be associated with autism spectrum disorders (ASD) is chromosomal abnormalities. The identification of copy number variants (CNV), i.e., microdeletions and microduplications that are undetectable at the level of traditional cytogenetic analysis, allows the potential association of submicroscopic chromosomal imbalances and human disease.

METHODS

We performed array comparative genomic hybridization (aCGH) utilizing a 19K whole genome tiling path bacterial artificial chromosome (BAC) microarray on 397 unrelated subjects with autism spectrum disorder. Common CNV were excluded using a control group comprised of 372 individuals from the National Institute of Mental Health (NIMH) Genetics Initiative Control samples. Confirmation studies were performed on all remaining CNV using fluorescence in situ hybridization (FISH), microsatellite analysis, and/or quantitative polymerase chain reaction (PCR) analysis.

RESULTS

A total of 51 CNV were confirmed in 46 ASD subjects. Three maternal interstitial duplications of 15q11-q13 known to be associated with ASD were identified. The other 48 CNV ranged in size from 189 kilobase (kb) to 5.5 megabase (Mb) and contained from 0 to approximately 40 National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) genes. Seven CNV were de novo and 44 were inherited.

CONCLUSIONS

Fifty-one autism-specific CNV were identified in 46 of 397 ASD patients using a 19K BAC microarray for an overall rate of 11.6%. These microdeletions and microduplications cause gene dosage imbalance in 272 genes, many of which could be considered as candidate genes for autism.

Recurrent 16p11.2 microdeletions in autism

Autism is a childhood neurodevelopmental disorder with a strong genetic component, yet the identification of autism susceptibility loci remains elusive. We investigated 180 autism probands and 372 control subjects by array comparative genomic hybridization (aCGH) using a 19K whole-genome tiling path bacterial artificial chromosome microarray to identify submicroscopic chromosomal rearrangements specific to autism. We discovered a recurrent 16p11.2 microdeletion in two probands with autism and none in controls. The deletion spans approximately 500-kb and is flanked by approximately 147-kb segmental duplications (SDs) that are >99% identical, a common characteristic of genomic disorders. We assessed the frequency of this new autism genomic disorder by screening an additional 532 probands and 465 controls by quantitative PCR and identified two more patients but no controls with the microdeletion, indicating a combined frequency of 0.6% (4/712 autism versus 0/837 controls; Fisher exact test P = 0.044). We confirmed all 16p11.2 deletions using fluorescence in situ hybridization, microsatellite analyses and aCGH, and mapped the approximate deletion breakpoints to the edges of the flanking SDs using a custom-designed high-density oligonucleotide microarray. Bioinformatic analysis localized 12 of the 25 genes within the microdeletion to nodes in one interaction network. We performed phenotype analyses and found no striking features that distinguish patients with the 16p11.2 microdeletion as a distinct autism subtype. Our work reports the first frequency, breakpoint, bioinformatic and phenotypic analyses of a de novo 16p11.2 microdeletion that represents one of the most common recurrent genomic disorders associated with autism to date.

Sequence variation in DOCK9 and heterogeneity in bipolar disorder

BACKGROUND

Linkage of bipolar disorder to a broad region on chromosome 13q has been supported in several studies including a meta-analysis on genome scans. Subsequent reports have shown that variations in the DAOA (G72) locus on 13q33 display association with bipolar disorder but these may not account for all of the linkage evidence in the region.

OBJECTIVE

To identify additional susceptibility loci on 13q32-q33 by linkage disequilibrium mapping and explore the impact of phenotypic heterogeneity on association.

METHODS

In the initial phase, 98 single nucleotide polymorphism (SNPs) located on 13q32-q33 were genotyped on 285 probands with bipolar disorder and their parents were drawn from families in the NIMH Genetics Initiative consortium for bipolar disorder (NIMH1-4) and two other series. Fine scale mapping using one family series (NIMH1-2) as the test sample was targeted on a gene that displayed the highest evidence of association. A secondary analysis of familial component phenotypes of bipolar disorder was conducted.

RESULTS

Three of seven SNPs in DOCK9, a gene that encodes an activator of the Rho-GTPase Cdc42, showed significant excess allelic transmission (P=0.0477-0.00067). Fine scale mapping on DOCK9 yielded evidence of association at nine SNPs in the gene (P=0.02-0.006). Follow-up tests detected excess transmission of the same allele of rs1340 in two out of three other sets of families. The association signals were largely attributable to maternally transmitted alleles (rs1927568: P=0.000083; odds ratio=3.778). A secondary analysis of familial component phenotypes of bipolar disorder detected significant association across multiple DOCK9 markers for racing thoughts, psychosis, delusion during mania and course of illness indicators.

CONCLUSION

These results suggest that DOCK9 contributes to both risk and increased illness severity in bipolar disorder. We found evidence for the effect of phenotypic heterogeneity on association. To our knowledge this is the first report to implicate DOCK9 or the Rho-GTPase pathway in the etiology of bipolar disorder.

No evidence for association between 19 cholinergic genes and bipolar disorder

Cholinergic dysfunction has been proposed for the pathogenesis of bipolar disorder (BD), and we have therefore performed a systematic association study of cholinergic system genes in BD (including schizoaffective disorder bipolar type). We genotyped 93 single nucleotide polymorphisms (SNPs) in 19 genes (CHAT, CHRM1-5, CHRNA1-7, CHRNA9, CHRNA10, and CHRNB1-4) in two series of samples: the National Institute of Mental Health (NIMH) Genetics Initiative pedigrees with 474 samples from 152 families, and the Clinical Neurogenetics (CNG) pedigrees with 83 samples from 22 multiplex families. Sib-transmission/disequilibrium test (sib_TDT) analysis showed nominally significant transmission bias for four SNPs (CHRNA2: rs7017417, P = 0.024; CHRNA5: rs514743, P = 0.031; CHRNB1: rs2302762, P = 0.049; CHRNB4: rs1948, P = 0.031). Haploview analyses showed nominally significant transmission bias of several haplotypes in CHRNA2, CHRNA7, CHRNB1, and CHRNB4, respectively. However, none of these associations reached gene-wide significance after correction by permutation. Alcohol dependence (including alcohol abuse) was not a significant covariate in the present genetic association analysis. Thus, it is unlikely that these 19 cholinergic genes play a major role in the pre-disposition to BD in these pedigrees.

Mapping genes that contribute to daunorubicin-induced cytotoxicity

Daunorubicin is an anthracycline antibiotic agent used in the treatment of hematopoietic malignancies. Toxicities associated with this agent include myelosuppression and cardiotoxicity; however, the genes or genetic determinants that contribute to these toxicities are unknown. We present an unbiased genome-wide approach that incorporates heritability, whole-genome linkage analysis, and linkage-directed association to uncover genetic variants contributing to the sensitivity to daunorubicin-induced cytotoxicity. Cell growth inhibition in 324 Centre d' Etude du Polymorphisme Humain lymphoblastoid cell lines (24 pedigrees) was evaluated following treatment with daunorubicin for 72 h. Heritability analysis showed a significant genetic component contributing to the cytotoxic phenotypes (h2 = 0.18-0.63 at 0.0125, 0.025, 0.05, 0.1, 0.2, and 1.0 mumol/L daunorubicin and at the IC50, the dose required to inhibit 50% cell growth). Whole-genome linkage scans at all drug concentrations and IC50 uncovered 11 regions with moderate peak LOD scores (> 1.5), including 4q28.2 to 4q32.3 with a maximum LOD score of 3.18. The quantitative transmission disequilibrium tests were done using 31,312 high-frequency single-nucleotide polymorphisms (SNP) located in the 1 LOD confidence interval of these 11 regions. Thirty genes were identified as significantly associated with daunorubicin-induced cytotoxicity (P < or = 2.0 x 10(-4), false discovery rate < or = 0.1). Pathway and functional gene ontology analysis showed that these genes were overrepresented in the phosphatidylinositol signaling system, axon guidance pathway, and GPI-anchored proteins family. Our findings suggest that a proportion of susceptibility to daunorubicin-induced cytotoxicity may be controlled by genetic determinants and that analysis using linkage-directed association studies with dense SNP markers can be used to identify the genetic variants contributing to cytotoxicity.

Genome-wide linkage scan of 98 bipolar pedigrees and analysis of clinical covariates

Despite compelling evidence that genetic factors contribute to bipolar disorder (BP), attempts to identify susceptibility genes have met with limited success. This may be due to the genetic heterogeneity of the disorder. We sought to identify susceptibility loci for BP in a genome-wide linkage scan with and without clinical covariates that might reflect the underlying heterogeneity of the disorder. We genotyped 428 subjects in 98 BP families at the Center for Inherited Disease Research with 402 microsatellite markers. We first carried out a non-parametric linkage analysis with MERLIN, and then reanalyzed the data with LODPAL to incorporate clinical covariates for age at onset (AAO), psychosis and comorbid anxiety. We sought to further examine the top findings in the covariate analysis in an independent sample of 64 previously collected BP families. In the non-parametric linkage analysis, three loci were nominally significant under a narrow diagnostic model and seven other loci were nominally significant under a broader model. The top findings were on chromosomes 2q24 and 3q28. The covariate analyses yielded additional evidence for linkage on 3q28 with AAO in the primary and independent samples. Although none of the linked loci were genome-wide significant, their congruence with prior results and, for the covariate analyses, their identification in two separate samples increases the likelihood that they are true positives and deserve further investigation. These findings further demonstrate the value of considering clinical features that may reflect the underlying heterogeneity of disease in order to facilitate gene mapping.

Attempted suicide in bipolar disorder pedigrees: evidence for linkage to 2p12

BACKGROUND

We are interested in identifying susceptibility genes that predispose subjects to attempted suicide.

METHODS

We conducted a secondary analysis of genome-wide linkage data from 162 bipolar pedigrees that incorporated attempted suicide as a clinical covariate.

RESULTS

The strongest covariate-based linkage signal was seen on 2p12 at marker D2S1777. The logarithm of odds (LOD) score at marker D2S1777 rose from 1.56 to 3.82 after inclusion of the suicide covariate, resulting in significant chromosome-wide empirically derived p-values for the overall linkage finding (p = .01) and for the change in LOD score after the inclusion of the covariate (p = .02).

CONCLUSIONS

The finding on chromosome 2 replicates results from two previous studies of attempted suicide in pedigrees with alcohol dependence and in pedigrees with recurrent early-onset depression. Combined, these three studies provide compelling evidence for a locus influencing attempted suicide on 2p12.

Are the arginine vasopressin V1a receptor microsatellites related to hypersexuality in children with a prepubertal and early adolescent bipolar disorder phenotype?

OBJECTIVE

To examine family-based transmission of the number of 5' flanking arginine vasopressin V1a receptor (AVPR1A) microsatellites, which include [(GATA)(14)] and complex [(CT)(4)-TT-(CT)(8)-(GT)(24)] repeats, in probands with a prepubertal and early adolescent bipolar disorder phenotype (PEA-BP). Preferential transmission of the number of AVPR1A microsatellite repeats to hypersexual and uninhibited people-seeking probands was hypothesized, based on reports from preclinical work in the literature.

METHODS

Probands were 83 participants in an ongoing controlled study of PEA-BP. The PEA-BP phenotype was defined by DSM-IV mania with at least one of the cardinal symptoms of mania (elation and/or grandiosity) to avoid diagnosing mania only by symptoms that overlapped with those for attention-deficit hyperactivity disorder (ADHD). Comprehensive assessment of the probands included separate Washington University in St. Louis Kiddie Schedule for Affective Disorders and Schizophrenia (WASH-U-KSADS) interviews of parents about their children and of children about themselves. Hypersexuality and uninhibited people-seeking were assessed from the corresponding WASH-U-KSADS items. Microsatellite genotyping of the AVPR1A repeats was conducted using fluorescently labeled primers and detected by laser-induced fluorescence. Alleles were determined with the assistance of semi-automated allele-calling software. There were 32 complete, biological trios (28 informative families) for the GATA repeat and 34 complete, biological trios (30 informative families) for the complex repeat. Data were analyzed using case-control and family-based association methods.

RESULTS

Preferential transmission of AVPR1A GATA or complex repeats was not significant for hypersexuality or uninhibited people-seeking, using the transmission disequilibrium test. Similarly, case-control analyses found no significant associations between hypersexuality or uninhibited people-seeking and the number of AVPR1A GATA or complex repeats. For p < 0.05, there was about 80% power to detect odds ratios of 5.0 and 4.0 (in the family-based analyses) and 3.5 and 2.6 (in the case-control analyses), for allele frequencies of 0.1 and 0.5, respectively.

CONCLUSION

Preferential transmission of AVPR1A to hypersexual or uninhibited people-seeking probands was not supported.

Familial variation in episode frequency in bipolar affective disorder

OBJECTIVE

Bipolar affective disorder is a familial illness characterized by recurrent episodes of mania and depression, but little is known about the familial nature of episode recurrence or its associated clinical features. The authors analyzed the recurrence frequency of affective episodes (episode frequency), along with associated clinical and demographic variables, in families with at least three members with a major affective disorder.

METHOD

Members of 86 families ascertained through probands with bipolar affective disorder who had two or more first-degree relatives with a major affective disorder were interviewed by psychiatrists and assigned an all-sources diagnosis. Data for 407 subjects with a major affective disorder were analyzed. Episode frequency was estimated as the number of episodes of major depression, mania, and hypomania per year of illness.

RESULTS

Episode frequency was smoothly distributed over the range of 0.02-20.2 episodes/year. Episode frequency was significantly correlated among relatives (r=0.56, p<0.004). Earlier age at onset, bipolar II disorder, hallucinations or delusions, alcoholism, and suicidal behavior were all more prevalent in the highest than in the lowest quartiles of episode frequency. Female gender and recurrent major depression were more prevalent in the lowest quartile. Panic disorder, substance abuse, and thyroid disease were all unrelated to episode frequency. Subjects with DSM-IV rapid cycling did not differ from other affected subjects for most of the variables tested.

CONCLUSIONS

Episode frequency is a highly familial trait in bipolar affective disorder, associated with several indicators of severity, and may be useful in defining clinical subtypes of bipolar affective disorder with greater genetic liability. DSM-IV rapid cycling was not supported by these data as the best predictor of familiality or severity.

Linkage disequilibrium analysis in the LOC93081-KDELC1-BIVM region on 13q in bipolar disorder

Genome-wide scans in bipolar disorder and a meta analysis on published data have provided evidence for linkage to chromosome 13q, although the reported peaks from various studies have not converged in a narrow region. Recently, single nucleotide polymorphisms (SNPs) at the G72/G30 locus have been shown to be associated with bipolar disorder suggesting its potential role in increasing disease risk. The proposed linkage region on 13q extends over a wide span, and could provide a clue to the existence of other susceptibility variants. In the present study, SNPs in the LOC93081-KDELC1-BIVM, a region proximal to G72, were interrogated in two bipolar family series. KDELC1 has a predicted filamin domain and BIVM contains an immunoglobulin-like motif. The small pedigree series yielded a nominally significant global P-value due to under-transmission of a rare haplotype but this finding was not supported by results from the larger series and in the case-control study that compared 278 cases and 277 controls.

Linkage disequilibrium of the brain-derived neurotrophic factor Val66Met polymorphism in children with a prepubertal and early adolescent bipolar disorder phenotype

OBJECTIVE

Transmission of the brain-derived neurotrophic factor (BDNF) Val66 allele in children with a prepubertal and early adolescent bipolar disorder phenotype was examined.

METHOD

The prepubertal and early adolescent bipolar disorder phenotype was defined as current DSM-IV bipolar I disorder (manic or mixed phase) with at least one cardinal mania criterion (i.e., euphoria and/or grandiosity) to ensure differentiation from attention deficit hyperactivity disorder. Probands (mean age=10.7 years, SD=2.7) were obtained by consecutive new case ascertainment from designated pediatric and psychiatric venues. Parents and probands were interviewed separately by research nurses who were blind to the probands' diagnoses. Genotyping was done with TaqMan Assay-on-Demand. Analysis was done with the Family Based Association Test program.

RESULTS

There were 53 complete, independent trios. The BDNF Val66 allele was preferentially transmitted (Family Based Association Test: chi(2)=6.0, df=1, p=0.014).

CONCLUSIONS

This finding in child bipolar disorder is consistent with data for adults with bipolar disorder that show preferential transmission of the Val66 allele.

Loci on chromosomes 6q and 6p interact to increase susceptibility to bipolar affective disorder in the national institute of mental health genetics initiative pedigrees

BACKGROUND

We have reported genetic linkage between bipolar disorder and markers on chromosome 6q16.3-22.1 in the National Institute of Mental Health Genetics Initiative wave 3 pedigrees. Here we test for: 1) robustness of the linkage to differing analysis methods, genotyping error, and gender-specific maps; 2) parent-of-origin effects; and 3) interaction with markers within the schizophrenia linkage region on chromosome 6p.

METHODS

Members of 245 families ascertained through a sibling pair affected with bipolar I or schizoaffective-bipolar disorder were genotyped with 18 markers spanning chromosome 6. Nonparametric linkage analysis was performed.

RESULTS

Linkage to 6q is robust to analysis method, gender-specific map differences, and genotyping error. The locus confers a 1.4-fold increased risk. Affected siblings share the maternal more often than the paternal chromosome (p =.006), which could reflect a maternal parent-of-origin effect. There is a positive correlation between family-specific linkage scores on 6q and those on 6p22.2 (r =.26; p <.0001). Linkage analysis for each locus conditioned on evidence of linkage to the other increases the evidence for linkage at both loci (p <.0005). Logarithm of the odds (LOD) scores increased from 2.26 to 5.42 on 6q and from.35 to 2.26 on 6p22.2.

CONCLUSIONS

These results support linkage of bipolar disorder to 6q, uncover a maternal parent-of-origin effect, and demonstrate an interaction of this locus with one on chromosome 6p22.2, previously linked only to schizophrenia.

Heritability and Linkage Analysis of Sensitivity to Cisplatin-Induced Cytotoxicity

Little is known about the genetic determinants explaining variation in sensitivity to chemotherapeutic cytotoxicity. We characterized the degree of cisplatin sensitivity, using lymphoblastoid cell lines derived from 10 Centre d'Etude du Polymorphisme Humain pedigrees. We estimated the heritability for susceptibility to cisplatin-induced cytotoxicity to be approximately 0.47; therefore, sensitivity to the cytotoxic effects of cisplatin is under appreciable genetic influence. Linkage analysis was performed, and the strongest signal (lod score, 2.16; empirical P = 0.0005) was found on chromosome 1 at 44 cM. Susceptibility to cisplatin-induced cytotoxicity is likely due to multiple loci, with low locus-specific heritability contributing to the trait. These data show the power of using large pedigrees that have been extensively genotyped for evaluating the genetic contribution to sensitivity to cell growth inhibition by anticancer agents.

Genome-wide scan and conditional analysis in bipolar disorder: evidence for genomic interaction in the National Institute of Mental Health genetics initiative bipolar pedigrees

BACKGROUND

In 1989 the National Institute of Mental Health began a collaborative effort to identify genes for bipolar disorder. The first 97 pedigrees showed evidence of linkage to chromosomes 1, 6, 7, 10, 16, and 22 (Nurnberger et al 1997). An additional 56 bipolar families have been genotyped, and the combined sample of 153 pedigrees studied.

METHODS

Three hierarchical affection status models were analyzed with 513 simple sequence repeat markers; 298 were common across all pedigrees. The primary analysis was a nonparametric genome-wide scan. We performed conditional analyses based on epistasis or heterogeneity for five regions.

RESULTS

One region, on 16p13, was significant at the genome-wide p <.05 level. Four additional chromosomal regions (20p12, 11p15, 6q24, and 10p12) showed nominally significant linkage findings (p </=.01). Conditional analysis assuming epistasis identified a significant increase in linkage at four regions. Families linked to 6q24 showed a significant increase in nonparametric logarithms of the odds (NPL) scores at 5q11 and 7q21. Epistasis also was observed between 20p12 and 13q21, and 16p13 and 9q21.

CONCLUSIONS

The findings are presented in rank order of nominal significance. Several of these regions have been previously implicated in independent studies of either bipolar disorder or schizophrenia. The strongest finding is at 16p13 at D16S748 with an NPL of 3.3, there is evidence of epistasis between this locus and 9q21. Application of conditional analyses is potentially useful in larger sample collections to identify susceptibility genes of modest influence that may not be identified in a genome-wide scan aimed to identify single gene effects.

Genome scan meta-analysis of schizophrenia and bipolar disorder, part III: Bipolar disorder

Genome scans of bipolar disorder (BPD) have not produced consistent evidence for linkage. The rank-based genome scan meta-analysis (GSMA) method was applied to 18 BPD genome scan data sets in an effort to identify regions with significant support for linkage in the combined data. The two primary analyses considered available linkage data for "very narrow" (i.e., BP-I and schizoaffective disorder-BP) and "narrow" (i.e., adding BP-II disorder) disease models, with the ranks weighted for sample size. A "broad" model (i.e., adding recurrent major depression) and unweighted analyses were also performed. No region achieved genomewide statistical significance by several simulation-based criteria. The most significant P values (<.01) were observed on chromosomes 9p22.3-21.1 (very narrow), 10q11.21-22.1 (very narrow), and 14q24.1-32.12 (narrow). Nominally significant P values were observed in adjacent bins on chromosomes 9p and 18p-q, across all three disease models on chromosomes 14q and 18p-q, and across two models on chromosome 8q. Relatively few BPD pedigrees have been studied under narrow disease models relative to the schizophrenia GSMA data set, which produced more significant results. There was no overlap of the highest-ranked regions for the two disorders. The present results for the very narrow model are promising but suggest that more and larger data sets are needed. Alternatively, linkage might be detected in certain populations or subsets of pedigrees. The narrow and broad data sets had considerable power, according to simulation studies, but did not produce more highly significant evidence for linkage. We note that meta-analysis can sometimes provide support for linkage but cannot disprove linkage in any candidate region.

Polymorphisms at the G72/G30 gene locus, on 13q33, are associated with bipolar disorder in two independent pedigree series

Linkage evidence suggests that chromosome 13 (13q32-33) contains susceptibility genes for both bipolar disorder and schizophrenia. Recently, genes called "G72" and "G30" were identified, and polymorphisms of these overlapping genes were reported to be associated with schizophrenia. We studied two series of pedigrees with bipolar disorder: the Clinical Neurogenetics (CNG) pedigrees (in which linkage to illness had been previously reported at 13q32-33), with 83 samples from 22 multiplex families, and the National Institute of Mental Health (NIMH) Genetics Initiative pedigrees, with 474 samples from 152 families. Sixteen single-nucleotide polymorphisms (SNPs) were genotyped at and around the G72/G30 locus, which covered a 157-kb region encompassing the entire complementary DNA sequences of G72 and G30. We performed transmission/disequilibrium testing (TDT) and haplotype analysis, since a linkage-disequilibrium block was present at this gene locus. In the CNG and NIMH data sets, the results of global TDT of the entire haplotype set were significant and consistent (P=.0004 and P=.008, respectively). In the CNG series, the associated genotypes divided the families into those with linkage and those without linkage (partitioned by the linkage evidence). Analysis of the decay of haplotype sharing gave a location estimate that included G72/G30 in its 95% confidence interval. Although statistically significant association was not detected for individual SNPs in the NIMH data set, the same haplotype was consistently overtransmitted in both series. These data suggest that a susceptibility variant for bipolar illness exists in the vicinity of the G72/G30 genes. Taken together with the earlier report, this is the first demonstration of a novel gene(s), discovered through a positional approach, independently associated with both bipolar illness and schizophrenia.