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.

Genome-wide linkage analysis of 972 bipolar pedigrees using single-nucleotide polymorphisms.

Because of the high costs associated with ascertainment of families, most linkage studies of Bipolar I disorder (BPI) have used relatively small samples. Moreover, the genetic information content reported in most studies has been less than 0.6. Although microsatellite markers spaced every 10 cM typically extract most of the genetic information content for larger multiplex families, they can be less informative for smaller pedigrees especially for affected sib pair kindreds. For these reasons we collaborated to pool family resources and carried out higher density genotyping. Approximately 1100 pedigrees of European ancestry were initially selected for study and were genotyped by the Center for Inherited Disease Research using the Illumina Linkage Panel 12 set of 6090 single-nucleotide polymorphisms. Of the ~1100 families, 972 were informative for further analyses, and mean information content was 0.86 after pruning for linkage disequilibrium. The 972 kindreds include 2284 cases of BPI disorder, 498 individuals with bipolar II disorder (BPII) and 702 subjects with recurrent major depression. Three affection status models (ASMs) were considered: ASM1 (BPI and schizoaffective disorder, BP cases (SABP) only), ASM2 (ASM1 cases plus BPII) and ASM3 (ASM2 cases plus recurrent major depression). Both parametric and non-parametric linkage methods were carried out. The strongest findings occurred at 6q21 (non-parametric pairs LOD 3.4 for rs1046943 at 119 cM) and 9q21 (non-parametric pairs logarithm of odds (LOD) 3.4 for rs722642 at 78 cM) using only BPI and schizoaffective (SA), BP cases. Both results met genome-wide significant criteria, although neither was significant after correction for multiple analyses. We also inspected parametric scores for the larger multiplex families to identify possible rare susceptibility loci. In this analysis, we observed 59 parametric LODs of 2 or greater, many of which are likely to be close to maximum possible scores. Although some linkage findings may be false positives, the results could help prioritize the search for rare variants using whole exome or genome sequencing.

An international two-stage genome-wide search for schizophrenia susceptibility genes.

Schizophrenia is thought to be a multifactorial disease with complex mode of inheritance. Using a two-stage strategy for another complex disorder, a number of putative IDDM-susceptibility genes have recently been mapped. We now report the results of a two-stage genome-wide search for genes conferring susceptibility to schizophrenia. In stage I, model-free linkage analyses of large pedigrees from Iceland, a geographical isolate, revealed 26 loci suggestive of linkage. In stage II, ten of these were followed-up in a second international collaborative study comprising families from Austria, Canada, Germany, Italy, Scotland, Sweden, Taiwan and the United States. Potential linkage findings of stage I on chromosomes 6p, 9 and 20 were observed again in the second sample. Furthermore, in a third sample from China, fine mapping of the 6p region by association studies also showed evidence for linkage or linkage disequilibrium. Combining our results with other recent findings revealed significant evidence for linkage to an area distal of the HLA region on chromosome 6p. However, in a fourth sample from Europe, the 6p fine mapping finding observed in the Chinese sample could not be replicated. Finally, evidence suggestive of locus heterogeneity and oligogenic transmission in schizophrenia was obtained.

Singleton deletions throughout the genome increase risk of bipolar disorder.

An overall burden of rare structural genomic variants has not been reported in bipolar disorder (BD), although there have been reports of cases with microduplication and microdeletion. Here, we present a genome-wide copy number variant (CNV) survey of 1001 cases and 1034 controls using the Affymetrix single nucleotide polymorphism (SNP) 6.0 SNP and CNV platform. Singleton deletions (deletions that appear only once in the dataset) more than 100 kb in length are present in 16.2% of BD cases in contrast to 12.3% of controls (permutation P=0.007). This effect was more pronounced for age at onset of mania

Meta-analysis of 32 genome-wide linkage studies of schizophrenia.

A genome scan meta-analysis (GSMA) was carried out on 32 independent genome-wide linkage scan analyses that included 3255 pedigrees with 7413 genotyped cases affected with schizophrenia (SCZ) or related disorders. The primary GSMA divided the autosomes into 120 bins, rank-ordered the bins within each study according to the most positive linkage result in each bin, summed these ranks (weighted for study size) for each bin across studies and determined the empirical probability of a given summed rank (P(SR)) by simulation. Suggestive evidence for linkage was observed in two single bins, on chromosomes 5q (142-168 Mb) and 2q (103-134 Mb). Genome-wide evidence for linkage was detected on chromosome 2q (119-152 Mb) when bin boundaries were shifted to the middle of the previous bins. The primary analysis met empirical criteria for 'aggregate' genome-wide significance, indicating that some or all of 10 bins are likely to contain loci linked to SCZ, including regions of chromosomes 1, 2q, 3q, 4q, 5q, 8p and 10q. In a secondary analysis of 22 studies of European-ancestry samples, suggestive evidence for linkage was observed on chromosome 8p (16-33 Mb). Although the newer genome-wide association methodology has greater power to detect weak associations to single common DNA sequence variants, linkage analysis can detect diverse genetic effects that segregate in families, including multiple rare variants within one locus or several weakly associated loci in the same region. Therefore, the regions supported by this meta-analysis deserve close attention in future studies.

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.

Genetic mapping of adrenergic receptor genes in humans.

We have genetically mapped the genes encoding four human adrenergic receptors (ARs) of subtypes alpha 1C, alpha 2A, alpha 2B, and beta 1, which are prototypic G protein coupled receptors that mediate the physiological effects of neurotransmitters, hormones, and drugs. We placed these genes onto the Cooperative Human Linkage Center (CHLC) and Genethon framework maps, within confidence intervals with greater than 1000:1 odds. With multipoint analysis the alpha 1C gene (locus ADRA1C) mapped to the interval between NEFL and D8S283; alpha 2-C4, the gene encoding the alpha 2C AR (locus ADRA2C), mapped to the interval between D4S126 and D4S62; and the alpha 2-C10 (alpha 2A AR)/beta 1 haplotype (loci ADRA2A/ADRB1) mapped to the interval between D10S259 and D10S187. A fifth AR gene, beta 2, yielded significant LOD scores with markers on the long arm of chromosome 5; however, this locus (ADRB2) could not be mapped to any specific interval with odds of greater than 1000:1. The two AR genes that are completely linked, alpha 2-C10 and beta 1, were oriented on their shared 225-kb genomic fragment relative to the direction of transcription, with beta 1 being 5' to alpha 2-C10. The positioning of these genes on high-density framework maps allows them to be tested as candidates in a spectrum of diseases that might involve AR dysfunction.

Developmental and genetic influences on the P50 sensory gating phenotype.

Evoked potentials to pairs of click stimuli were recorded from 127 subjects ranging in age from 10 to 39 years to examine the developmental course of auditory sensory gating. The ratio of the amplitude of the second response to that of the first provides a quantitative measure of auditory sensory gating. Contrary to earlier results, the distribution of P50 ratios was unchanged between children and younger adolescents (10-14 years), older adolescents (15-19 years), and adults (20-29 and 30-39 years). Included in the sample were 39 adolescent twins, allowing assessment for possible genetic effects underlying the P50 sensory gating phenotype, by comparison of the similarity of the measure in monozygotic and same-sex dizygotic twin pairs. The monozygotic twins had significantly higher similarity for the P50 ratio within each twin pair than the dizygotic twins. These results are consistent with the presence of genetic influences on the P50 sensory gating phenotype.

Analysis of chromosome 18 DNA markers in multiplex pedigrees with manic depression.

Six pedigrees segregating manic-depressive illness (MDI) were analyzed for linkage to 21 highly polymorphic microsatellite DNA markers on chromosome 18. These markers span almost the entire length of the chromosome, and gaps between markers are less than 20 cM. In particular, we analyzed several markers localizing to the pericentromeric region of chromosome 18 which generated lod scores suggestive of linkage in an independent study. Lod score analysis was performed and results were examined by family. One region produced positive lod scores, though at 18q23 and not in the pericentromeric region. We additionally used two nonparametric methods because the true mode of transmission of MDI is unknown; results were again somewhat suggestive for markers in the region of 18q23 but not in the pericentromeric region.

Use of a neurophysiological trait in linkage analysis of schizophrenia.

Traditional diagnostic techniques may not provide all the information necessary to reveal the genetic causes of schizophrenia through linkage analysis. Use of neurophysiological indicator variables that are associated with the disease may increase the probability of detecting linkage. Such variables not only produce simpler phenotypes for analysis, but they also may be more proximal to the gene products involved in neurological dysfunctions underlying schizophrenia. We have used a previously characterized neurophysiological variable, the P50 evoked-auditory response, to search for chromosomal regions that may be of interest in the study of schizophrenia. Although our scan of over 300 markers did not show strong evidence for linkage to P50 in nine families, this exploratory analysis has revealed several chromosomal regions that may deserve further study.

Search for a gene predisposing to manic-depression on chromosome 21.

Six kindreds containing multiple cases of Manic-depressive illness (MDI) were genotyped with seven highly polymorphic microsatellite loci used in the construction of an index map for chromosome 21. The kindreds were also genotyped with a microsatellite polymorphism for PFKL, a chromosome 21 locus that has shown suggestive linkage to MDI in one pedigree [Straub et al., 1993: The American Society of Human Genetics]. Evidence of linkage was not found assuming either autosomal dominant or recessive inheritance. The nonparametric affected sib pair test did not yield significant evidence of linkage.

Sensory gating deficits in parents of schizophrenics.

Although schizophrenia clusters in families, it is not inherited in Mendelian fashion. This suggests that there may be alternative phenotypic expressions of genes that convey risk for schizophrenia, such as more elementary physiological or biochemical defects. One proposed phenotype is impaired inhibitory gating of the auditory evoked potential to repeated stimuli. Normally, the amplitude of the P50 response to the second stimulus is significantly less than the response to the first, but this gating of response is generally impaired in schizophrenia. Clinically unaffected individuals within a pedigree who have both an ancestral and descendant history of schizophrenia may be useful for studying whether this physiological defect is a possible alternative phenotype. We have studied inhibitory gating of the auditory P50 response to pairs of auditory stimuli in 17 nuclear families. In 11, there was one parent who had another relative with a chronic psychotic illness, in addition to the schizophrenic proband. All of the parents with family histories of schizophrenia had gating of the P50 response similar to their schizophrenia offspring, whereas only 7% of the parents without family history had gating of the P50 response in the abnormal range. These results support loss of gating of the auditory P50 wave as an inherited deficit related to schizophrenia and suggest that studies of parents may help elucidate the neurobiological expression of genes that convey risk for schizophrenia.

Genotype-to-phenotype analysis: search for clinical characteristics of a missense change in the GABAA-beta 1 receptor gene.

Genotype-to-phenotype analysis reverses the classical approach to genetic disease in which an unknown genotype is sought for a known phenotype. This paper provides an example of genotype-to-phenotype analysis for the possible psychiatric effects of a missense mutation (H396Q) at a highly conserved residue of the beta 1 subunit gene of the gamma aminobutyric acid type A receptor. DNA samples from 1,507 Caucasians of Western European descent were screened, and 10 heterozygotes for H396Q were identified. These individuals were matched to homozygous normal individuals by age, gender, and length of available medical records. The complete medical records of these 20 individuals were reviewed blindly by two psychiatrists (D.C.S., L.L.H.) to assess psychiatric symptomatology, with an emphasis on anxiety and related disorders. However, no association was found between this missense change at a conserved amino acid and a dominant neuropsychiatric disease phenotype. Thus, this missense change may be neutral or only mildly deleterious, may only cause recessive disease in rare individuals, or may interact epistatically with some other gene(s).

Mutation scan of the D1 dopamine receptor gene in 22 cases of bipolar I disorder.

In a previous study [Coon et al., Am J Hum Genet 52:1234-1249, 1993], we found suggestive evidence of linkage between manic-depressive illness (MDI) in eight multiplex pedigrees and D5S62, a DNA marker mapping to the telomeric region of 5q. As the D1 dopamine receptor gene (DRD1) maps to this region and as alterations in dopaminergic neurotransmission have been indirectly implicated in the pathogenesis of MDI, we directly searched for mutations in the coding region of the DRD1 gene in 22 unrelated cases of bipolar I (BPI) disorder derived from multiplex families. Using single strand conformation polymorphism (SSCP) analysis, we did not observe any abnormal SSCP variants in the BPI cases that differed from controls.

Search for a schizophrenia susceptibility locus on human chromosome 22.

We used 10 highly informative DNA polymorphic markers and genetic linkage analysis to examine whether a gene locus predisposing to schizophrenia is located on chromosome 22, in 105 families with schizophrenia and schizoaffective disorder. The LOD score method, including analysis for heterogeneity, provided no conclusive evidence of linkage under a dominant, recessive, or penetrance free model of inheritance. Affected sib-pair analysis was inconclusive. Affected pedigree member analysis gave only suggestive evidence for linkage. Multipoint APM analysis, using 4 adjacent loci including D22S281 and IL2RB, a region of interest from the APM analysis, gave non-significant results for the three different weighting functions.

Analysis of chromosome 22 markers in nine schizophrenia pedigrees.

Previous results of a genome-wide survey for schizophrenia susceptibility genes in nine multiplex families indicated a possible region of linkage on chromosome 22. We therefore tested for linkage using ten highly polymorphic chromosome 22 DNA markers. Lod score analyses were suggestive of linkage for several markers on the distal end of the chromosome; however, no lod score exceeded 3 assuming either autosomal dominant or autosomal recessive transmission. The highest lod score was 2.09 (theta = 0.10) for marker D22S276 under autosomal recessive inheritance. Based on simulation analyses, this result is unlikely to represent a false positive. Analyses using information from affected individuals only resulted in reduced lod scores, with a maximum of 1.40 (theta = 0.05) for D22S276 assuming autosomal recessive inheritance. Two nonparametric methods, sib pair analysis and the Affected-Pedigree-Member method, also yielded suggestive but inconclusive findings; results were positive, but strict thresholds of significance were not met. Additionally, we tested one candidate gene, the Arylsulfatase A gene, located in the region of 22q13.31-qter. Results were again inconclusive, though the DNA marker available for this gene was a 2-allele RFLP with heterozygosity of 0.5, and therefore not maximally informative. Further investigation of this chromosomal region and this and other candidate genes may be warranted.

Linkage of a composite inhibitory phenotype to a chromosome 22q locus in eight Utah families.

Eight Utah multigenerational families, each with three to six cases of schizophrenia, were phenotyped with two specific measures of inhibitory neurophysiological functioning, P50 auditory sensory gating (P50), and antisaccade ocular motor performance (AS). A genomewide linkage analysis was performed to screen for loci underlying a qualitative phenotype combining the P50 and AS measures. For this composite inhibitory phenotype, the strongest evidence for linkage was to the D22s315 marker on chromosome 22q (lod score = 3.55, theta = 0) under an autosomal dominant model. Simulation analyses indicate that this 3.55 lod score is unlikely to represent a false positive result. Lod scores were 2.0 or greater for markers flanking D22s315. A nonparametric linkage (NPL) analysis of the chromosome 22 data showed evidence for allele sharing over the broad region surrounding D22s315 with a maximum NPL score of 3.83 (p = .002) for all pedigrees combined.

Genetic epidemiological study of schizophrenia in Palau, Micronesia: prevalence and familiality.

We are studying the genetic etiology of schizophrenia in the Republic of Palau, a remote island nation in Micronesia that has been geographically and ethnically isolated for approximately 2,000 years. The first epidemiological phase sought to estimate the lifetime prevalence of schizophrenia and evaluate the familiality of the illness based on complete ascertainment of cases and families segregating schizophrenia. A total of 160 strictly defined cases of schizophrenia were ascertained in a population of 13,750 adults who were 15 years of age and older. The lifetime prevalence of strictly defined schizophrenia in Palau was 1.99% overall and 2.77% in males vs. 1.24% in females. This greater than 2:1 male-to-female risk ratio for schizophrenia was accompanied by an earlier mean age of onset for males (23.3 years) than for females (27.5 years). These 160 cases of strict schizophrenia represent 59 separate families each identified by a single common founder. Eleven of these families have 5 to 14 cases and represent nearly half of the strict schizophrenia cases in Palau. Although schizophrenia is clearly aggregating in these 11 families, cases are distributed sparsely throughout the large sibships. In the entire sample of 160 cases of strict schizophrenia, there were only 11 sib-pairs and 2 sib-trios. When a family was defined to include third-degree relatives, only 11 cases (6.9%) were nonfamilial. The majority of the ascertained cases can be linked together into extended pedigrees with complex multilineal inheritance patterns. These intricately interconnected families may pose challenges for traditional linkage techniques. However, these Palauan families represent a valuable resource for studying the genetic etiology of schizophrenia because there may be fewer susceptibility genes for schizophrenia in this genetic isolate than in the heterogeneous populations that are common throughout the world today.

A combined analysis of D22S278 marker alleles in affected sib-pairs: support for a susceptibility locus for schizophrenia at chromosome 22q12. Schizophrenia Collaborative Linkage Group (Chromosome 22).

Several groups have reported weak evidence for linkage between schizophrenia and genetic markers located on chromosome 22q using the lod score method of analysis. However these findings involved different genetic markers and methods of analysis, and so were not directly comparable. To resolve this issue we have performed a combined analysis of genotypic data from the marker D22S278 in multiply affected schizophrenic families derived from 11 independent research groups worldwide. This marker was chosen because it showed maximum evidence for linkage in three independent datasets (Vallada et al., Am J Med Genet 60:139-146, 1995; Polymeropoulos et al., Neuropsychiatr Genet 54:93-99, 1994; Lasseter et al., Am J Med Genet, 60:172-173, 1995. Using the affected sib-pair method as implemented by the program ESPA, the combined dataset showed 252 alleles shared compared with 188 alleles not share (chi-square 9.31, 1df, P = 0.001) where parental genotype data was completely known. When sib-pairs for whom parental data was assigned according to probability were included the number of alleles shared was 514.1 compared with 437.8 not shared (chi-square 6.12, 1df, P = 0.006). Similar results were obtained when a likelihood ratio method for sib-pair analysis was used. These results indicate that may be a susceptibility locus for schizophrenia at 22q12.

Treatment of seasonal affective disorder with morning light.

Three subjects afflicted with seasonal affective disorder (winter depression) were treated with 2-hour morning light exposures. Within 2 to 5 days, all responded, and remission of their symptoms was sustained during the 2-month treatment period.