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.

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.

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.

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

Suggestive linkage of schizophrenia to 5p13 in Costa Rica.

Schizophrenia afflicts roughly 1% of all people worldwide. Remarkably, despite differing cultures and environments, the expression of illness is essentially the same. Family, twin, and adoption studies identify schizophrenia as a genetically influenced disease. Linkage studies suggest many positive regions of interest, but as a complex genetic disorder most of the pathogenic loci have not yet been found. Isolated populations are commonly used to study rare Mendelian inherited diseases due to the more homogenous genetic background of the subjects and are thought to be useful for detecting linkage in complex genetic disorders such as schizophrenia. This study aims to define areas of the genome that exhibit co-inheritance with schizophrenia in one large, Mendelian-like family from the central valley of Costa Rica. The whole genome scan analysis of this pedigree, which included 11 cases of schizophrenia and schizoaffective disorder, identified a number of markers on chromosome 5p that appear to co-segregate with the disease with a maximum lod score of 2.70 at marker D5S426. Current studies include investigating additional Costa Rican pedigrees to replicate these findings and identify additional loci linked to the disease.

Genome-wide multipoint linkage analyses of multiplex schizophrenia pedigrees from the oceanic nation of Palau.

The oceanic nation of Palau has been geographically and culturally isolated over most of its 2000 year history. As part of a study of the genetic basis of schizophrenia in Palau, we genotyped five large, multigenerational schizophrenia pedigrees using markers every 10 cM (CHLC/Weber screening set 6). The number of affected/unaffected individuals genotyped per family ranged from 11/21 to 5/5. Thus the pedigrees varied in their information for linkage, but each was capable of producing a substantial LOD score. We fitted a simple dominant and recessive model to these data using multipoint linkage analysis implemented by Simwalk2. Predictably, the most informative pedigrees produced the best linkage results. After genotyping additional markers in the region, one pedigree produced a LOD = 3.4 (5q distal) under the dominant model. Seven of nine schizophrenics in the pedigree, mostly 3rd-4th degree relatives, share a 15-cM, 7-marker haplotype. For a different pedigree, another promising signal occurred on distal 3q, LOD = 2.6, for the recessive model. For two other pedigrees, the best LODs were modest, slightly better than 2.0 on 5q and 9p, while the fifth pedigree produced no noteworthy linkage signal. Similar to the results for other populations, our results suggest there are multiple genes conferring liability to schizophrenia even in the small population of Palau (roughly 21,000 individuals) in remote Oceania.

Genome-wide distribution of linkage disequilibrium in the population of Palau and its implications for gene flow in Remote Oceania.

Linkage disequilibrium (LD) between alleles on the same human chromosome results from various evolutionary processes and is thus telling about the history of populations. Recently, LD has garnered substantial interest for its value to map and fine-map disease genes. We examine the distribution of LD between short tandem repeat alleles on autosomes and sex chromosomes in the Remote Oceanic population of Palau to evaluate whether the data are consistent with a recent hypothesis about the origins of genetic variation in Palau, specifically that the population experienced extensive male-biased gene flow following initial settlement. Consistent with evolutionary theory based on effective population size, LD between X-linked alleles is stochastically greater than LD between autosomal alleles, however, small but detectable LD occurs for autosomal markers separated by substantial distances. By contrast, while Y-linked alleles experience only one-third the effective population size of X-linked alleles, their mean value for pairwise LD is only slightly larger than X-linked alleles. For a small population known to experience at least two extreme bottlenecks, 56 six-locus Y haplotypes exhibit remarkable diversity (0.96), comparable to Y diversity of Europeans, however, autosomal and X-linked markers display significantly less diversity, as measured by heterozygosity (4.1% less). Palauan Y haplotypes also fall into distinct clusters, again unlike that of Europe. We argue these data are consistent with waves of male-biased gene flow.

Ocular motor delayed-response task performance among patients with schizophrenia and their biological relatives.

Schizophrenia patients and their relatives have saccadic abnormalities characterized by problems inhibiting a response. The dorsolateral prefrontal cortex and its associated circuitry ostensibly mediate inhibition and support correct delayed response performance. In this context, two components of delayed response task performance are of interest: memory saccade metrics and error saccades made during the delay. To evaluate these variables, an ocular motor delayed response task was presented to 23 schizophrenia patients, 25 of their first-degree biological relatives, and 19 normal subjects. The measure that best differentiated groups was an increased frequency of error saccades generated during the delay by schizophrenia subjects and relatives. Decreased memory saccade gain also characterized patients and relatives. The similar pattern of results demonstrated by the patients with schizophrenia and their relatives suggests that performance on ocular motor delayed response tasks, either alone or in combination with other saccadic variables, may provide useful information about neural substrates associated with a liability for developing schizophrenia.

The ascertainment of multiplex schizophrenia pedigrees from Daghestan genetic isolates (Northern Caucasus, Russia).

This article describes the preliminary ascertainment of multiplex schizophrenia pedigrees from the isolated mountain region of Daghestan (Northern Caucasus, Russia). Daghestan has a population of two million people and contains 26 aboriginal ethnic groups. Many of the ethnic groups reside in remote mountain villages that can be classified as 'primary isolates'. Prolonged reproductive isolation and severe environmental conditions in the highlands have created diverse, genetically isolated ethnic populations in Daghestan. A number of the isolates in this region contain large extended multiplex schizophrenia pedigrees that are ideal for genetic analyses. During summer expeditions of 1996 and 1997, 14 separate large multiplex schizophrenia pedigrees were ascertained from 14 different mountain villages. Of the 14 kindreds, one had 50 schizophrenic cases available for ascertainment, one had 32, and another had 24. Seven of the remaining pedigrees had between 11 and 23 living cases. Within the kindreds, the number of males with chronic schizophrenia was at least twice that of females. The average age of onset of schizophrenia is 21.2 years for offspring of consanguineous marriages and 17.4 years for offspring of nonconsanguineous marriages (P = 0.033). Although the pedigrees ascertained from the remote mountain villages may not be representative of the general population, they are unique kindreds for mapping schizophrenia susceptibility genes.

Mutation screening of a neutral amino acid transporter, ASCT1, and its potential role in schizophrenia.

In a previous study, a genome scan of a subset of schizophrenia families from Palau, Micronesia gave evidence suggestive of linkage to microsatellite markers at 2p13-14. In addition, in a large extended multiplex pedigree (K1583), an 11 cM 2p13-14 haplotype segregated with the illness in eight distantly related schizophrenics. The haplotype region includes a neutral amino acid transporter, ASCT1. We mutation-screened the coding region, flanking intronic sequence and 5'-untranslated region of this transporter in affected members of K1583, two Palauan controls and one Caucasian control. Most polymorphisms were found to be silent or common to all samples scanned. A G/A heterozygote within intron 3 was found in one schizophrenic member of K1583, but was not found in any of the other affected members of K1583. A G/A heterozygote within intron 6 was found in two of six schizophrenics tested in K1583, and in one control. As none of the sequence polymorphisms segregated with illness in the eight schizophrenics, it is unlikely that changes in the 5'-untranslated region, coding sequence or flanking intronic sequence of the ASCT gene predispose to schizophrenia in these families.

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.

Measuring liability for schizophrenia using optimized antisaccade stimulus parameters.

The ability to identify unaffected gene carriers within families may be crucial to the success of schizophrenia genetics studies. Data collected from three family samples (N = 365) demonstrated that poor antisaccade performance is an exceptionally promising indicator of liability for schizophrenia. A particular antisaccade task version provides large separations (5-6 sigma) between proband and normal groups. Poor antisaccade performance alone correctly identified 70% of patients in California, Utah, and Micronesia schizophrenia samples. Twenty-five to 50% of these patients' nonpsychotic first-degree relatives also had poor antisaccade performance, yielding risk ratios around 20:1 for simplex and 50:1 for multiplex schizophrenia families. Poor antisaccade performance is associated with dorsolateral prefrontal cortex pathology, suggesting that dysfunction of this circuitry also may predispose individuals to developing this disease.

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.

Evidence for a chromosome 2p13-14 schizophrenia susceptibility locus in families from Palau, Micronesia.

A large multiplex schizophrenia pedigree ascertained from the Micronesian nation of Palau was genotyped with 406 microsatellite DNA markers evenly distributed throughout the genome. Assuming autosomal dominant inheritance, the highest genome-wide lod scores were found for DNA loci mapping to 2p13-14; the maximum lod score was 2.17 (theta = 0.05) at D2S441. A nonparametric APM analysis was also suggestive at D2S441 (APM score = 2.96, P = 0.011). Of the 14 affected cases in this extended family, eight share a large haplotype in this region spanning approximately 11 cM. When 16 other families containing 65 schizophrenic cases were typed in a follow-up study of this region, the maximum lod score remained positive (maximum at D2S441 1.69, theta = 0.20). APM results also remained positive at D2S441 for all 17 families (APM score = 4.87, P = 0.0006). The linkage and haplotype sharing results provide suggestive evidence for a 2p locus predisposing to schizophrenia in a subset of families in the Palauan population.

Analysis of chromosome 6 markers in eight Utah schizophrenia pedigrees.

Eight multiplex Utah schizophrenia pedigrees were screened for linkage by applying a non-parametric linkage program, GENEHUNTER, using 30 chromosome 6 DNA markers. The overall maximum NPL score of the combined pedigrees was 1.50 (P = 0.079) at marker D6s281, located near the q terminus. The highest overall maximum NPL score for an individual pedigree was 2.81 (P = 0.043). In the chromosome 6p region, where numerous positive findings have been reported, we obtained no positive results.

The sensitivity of the Spontaneous Selective Attention Test (SSAT): a study of schizophrenic inpatients and outpatients versus normal controls.

The Spontaneous Selective Attention Task (SSAT) is a visual word-identification task that measures the type of selective attention that occurs spontaneously when there are multiple stimuli, all potentially relevant, and insufficient time to process each of them fully. The present study was designed to examine the sensitivity of the SSAT by comparing the performance of 40 schizophrenic inpatients and 30 schizophrenic outpatients to that of 70 normal controls. The pattern of results reported previously for schizophrenic inpatients versus normals was replicated, and these findings were extended to include schizophrenic outpatients in partial symptom remission. Schizophrenic inpatients and outpatients were just as accurate in identifying words as normals, but spontaneous selective attention under conditions of predictability was abnormal in both patient groups. Furthermore, the ability of the SSAT to discriminate between schizophrenic patients and controls was confirmed. A ratio measure of spontaneous selective attention had a sensitivity of 77% and a base rate of 9% in a normal population (when a cutoff value was set to minimize false positives and false negatives). Thus, the SSAT is a sensitive measure of an attentional phenotype that may be useful in genetic studies of schizophrenia.