The time and place of European admixture in Ashkenazi Jewish history.

The Ashkenazi Jewish (AJ) population is important in genetics due to its high rate of Mendelian disorders. AJ appeared in Europe in the 10th century, and their ancestry is thought to comprise European (EU) and Middle-Eastern (ME) components. However, both the time and place of admixture are subject to debate. Here, we attempt to characterize the AJ admixture history using a careful application of new and existing methods on a large AJ sample. Our main approach was based on local ancestry inference, in which we first classified each AJ genomic segment as EU or ME, and then compared allele frequencies along the EU segments to those of different EU populations. The contribution of each EU source was also estimated using GLOBETROTTER and haplotype sharing. The time of admixture was inferred based on multiple statistics, including ME segment lengths, the total EU ancestry per chromosome, and the correlation of ancestries along the chromosome. The major source of EU ancestry in AJ was found to be Southern Europe (≈60-80% of EU ancestry), with the rest being likely Eastern European. The inferred admixture time was ≈30 generations ago, but multiple lines of evidence suggest that it represents an average over two or more events, pre- and post-dating the founder event experienced by AJ in late medieval times. The time of the pre-bottleneck admixture event, which was likely Southern European, was estimated to ≈25-50 generations ago.

Identification of a missense variant in the WFS1 gene that causes a mild form of Wolfram syndrome and is associated with risk for type 2 diabetes in Ashkenazi Jewish individuals.

AIMS/HYPOTHESIS: Wolfram syndrome is a rare, autosomal recessive syndrome characterised by juvenile-onset diabetes and optic atrophy and is caused by bi-allelic mutations in the WFS1 gene. In a recent sequencing study, an individual with juvenile-onset diabetes was observed to be homozygous for a rare missense variant (c.1672C>T, p.R558C) in the WFS1 gene. The aim of this study was to perform the genetic characterisation of this variant and to determine whether it is causal for young-onset diabetes and Wolfram syndrome. METHODS: We analysed the allele frequency of the missense variant in multiple variant databases. We genotyped the variant in 475 individuals with type 1 diabetes and 2237 control individuals of Ashkenazi Jewish ancestry and analysed the phenotypes of homozygotes. We also investigated the association of this variant with risk for type 2 diabetes using genotype and sequence data for type 2 diabetes cases and controls. RESULTS: The missense variant demonstrated an allele frequency of 1.4% in individuals of Ashkenazi Jewish ancestry, 60-fold higher than in other populations. Genotyping of this variant in 475 individuals diagnosed with type 1 diabetes identified eight homozygotes compared with none in 2237 control individuals (genotype relative risk 135.3, p = 3.4 × 10-15). The age at diagnosis of diabetes for these eight individuals (17.8 ± 8.3 years) was several times greater than for typical Wolfram syndrome (5 ± 4 years). Further, optic atrophy was observed in only one of the eight individuals, while another individual had the Wolfram syndrome-relevant phenotype of neurogenic bladder. Analysis of sequence and genotype data in two case-control cohorts of Ashkenazi ancestry demonstrated that this variant is also associated with an increased risk of type 2 diabetes in heterozygotes (OR 1.81, p = 0.004). CONCLUSIONS/INTERPRETATION: We have identified a low-frequency coding variant in the WFS1 gene that is enriched in Ashkenazi Jewish individuals and causes a mild form of Wolfram syndrome characterised by young-onset diabetes and reduced penetrance for optic atrophy. This variant should be considered for genetic testing in individuals of Ashkenazi ancestry diagnosed with young-onset non-autoimmune diabetes and should be included in Ashkenazi carrier screening panels.

High-depth whole genome sequencing of an Ashkenazi Jewish reference panel: enhancing sensitivity, accuracy, and imputation.

While increasingly large reference panels for genome-wide imputation have been recently made available, the degree to which imputation accuracy can be enhanced by population-specific reference panels remains an open question. Here, we sequenced at full-depth (≥ 30×), across two platforms (Illumina X Ten and Complete Genomics, Inc.), a moderately large (n = 738) cohort of samples drawn from the Ashkenazi Jewish population. We developed a series of quality control steps to optimize sensitivity, specificity, and comprehensiveness of variant calls in the reference panel, and then tested the accuracy of imputation against target cohorts drawn from the same population. Quality control (QC) thresholds for the Illumina X Ten platform were identified that permitted highly accurate calling of single nucleotide variants across 94% of the genome. QC procedures also identified numerous regions that are poorly mapped using current reference or alternate assemblies. After stringent QC, the population-specific reference panel produced more accurate and comprehensive imputation results relative to publicly available, large cosmopolitan reference panels, especially in the range of rare variants that may be most critical to further progress in mapping of complex phenotypes. The population-specific reference panel also permitted enhanced filtering of clinically irrelevant variants from personal genomes.

A Frameshift in CSF2RB Predominant Among Ashkenazi Jews Increases Risk for Crohn's Disease and Reduces Monocyte Signaling via GM-CSF.

BACKGROUND & AIMS: Crohn's disease (CD) has the highest prevalence in Ashkenazi Jewish populations. We sought to identify rare, CD-associated frameshift variants of high functional and statistical effects. METHODS: We performed exome sequencing and array-based genotype analyses of 1477 Ashkenazi Jewish individuals with CD and 2614 Ashkenazi Jewish individuals without CD (controls). To validate our findings, we performed genotype analyses of an additional 1515 CD cases and 7052 controls for frameshift mutations in the colony-stimulating factor 2-receptor ß common subunit gene (CSF2RB). Intestinal tissues and blood samples were collected from patients with CD; lamina propria leukocytes were isolated and expression of CSF2RB and granulocyte-macrophage colony-stimulating factor-responsive cells were defined by adenomatous polyposis coli (APC) time-of-flight mass cytometry (CyTOF analysis). Variants of CSF2RB were transfected into HEK293 cells and the expression and functions of gene products were compared. RESULTS: In the discovery cohort, we associated CD with a frameshift mutation in CSF2RB (P = 8.52 × 10(-4)); the finding was validated in the replication cohort (combined P = 3.42 × 10(-6)). Incubation of intestinal lamina propria leukocytes with granulocyte-macrophage colony-stimulating factor resulted in high levels of phosphorylation of signal transducer and activator of transcription (STAT5) and lesser increases in phosphorylation of extracellular signal-regulated kinase and AK straining transforming (AKT). Cells co-transfected with full-length and mutant forms of CSF2RB had reduced pSTAT5 after stimulation with granulocyte-macrophage colony-stimulating factor, compared with cells transfected with control CSF2RB, indicating a dominant-negative effect of the mutant gene. Monocytes from patients with CD who were heterozygous for the frameshift mutation (6% of CD cases analyzed) had reduced responses to granulocyte-macrophage colony-stimulating factor and markedly decreased activity of aldehyde dehydrogenase; activity of this enzyme has been associated with immune tolerance. CONCLUSIONS: In a genetic analysis of Ashkenazi Jewish individuals, we associated CD with a frameshift mutation in CSF2RB. Intestinal monocytes from carriers of this mutation had reduced responses to granulocyte-macrophage colony-stimulating factor, providing an additional mechanism for alterations to the innate immune response in individuals with CD.

Identification of loci associated with schizophrenia by genome-wide association and follow-up.

We carried out a genome-wide association study of schizophrenia (479 cases, 2,937 controls) and tested loci with P < 10(-5) in up to 16,726 additional subjects. Of 12 loci followed up, 3 had strong independent support (P < 5 x 10(-4)), and the overall pattern of replication was unlikely to occur by chance (P = 9 x 10(-8)). Meta-analysis provided strongest evidence for association around ZNF804A (P = 1.61 x 10(-7)) and this strengthened when the affected phenotype included bipolar disorder (P = 9.96 x 10(-9)).

Excess of homozygosity in the major histocompatibility complex in schizophrenia.

Genome-wide association studies (GWAS) in schizophrenia have focused on additive allelic effects to identify disease risk loci. In order to examine potential recessive effects, we applied a novel approach to identify regions of excess homozygosity in an ethnically homogenous cohort: 904 schizophrenia cases and 1640 controls drawn from the Ashkenazi Jewish (AJ) population. Genome-wide examination of runs of homozygosity identified an excess in cases localized to the major histocompatibility complex (MHC). To refine this signal, we used the recently developed GERMLINE algorithm to identify chromosomal segments shared identical-by-descent (IBD) and compared homozygosity at such segments in cases and controls. We found a significant excess of homozygosity in schizophrenia cases compared with controls in the MHC (P-value = 0.003). An independent replication cohort of 548 schizophrenia cases from Japan and 542 matched healthy controls demonstrated similar effects. The strongest case-control recessive effects (P = 8.81 × 10(-8)) were localized to a 53-kb region near HLA-A, in a segment encompassing three poorly annotated genes, TRIM10, TRIM15 and TRIM40. At the same time, an adjacent segment in the Class I MHC demonstrated clear additive effects on schizophrenia risk, demonstrating the complexity of association in the MHC and the ability of our IBD approach to refine localization of broad signals derived from conventional GWAS. In sum, homozygosity in the classical MHC region appears to convey significant risk for schizophrenia, consistent with the ecological literature suggesting that homozygosity at the MHC locus may be associated with vulnerability to disease.

Genome-wide mapping of IBD segments in an Ashkenazi PD cohort identifies associated haplotypes.

The recent series of large genome-wide association studies in European and Japanese cohorts established that Parkinson disease (PD) has a substantial genetic component. To further investigate the genetic landscape of PD, we performed a genome-wide scan in the largest to date Ashkenazi Jewish cohort of 1130 Parkinson patients and 2611 pooled controls. Motivated by the reduced disease allele heterogeneity and a high degree of identical-by-descent (IBD) haplotype sharing in this founder population, we conducted a haplotype association study based on mapping of shared IBD segments. We observed significant haplotype association signals at three previously implicated Parkinson loci: LRRK2 (OR = 12.05, P = 1.23 × 10(-56)), MAPT (OR = 0.62, P = 1.78 × 10(-11)) and GBA (multiple distinct haplotypes, OR > 8.28, P = 1.13 × 10(-11) and OR = 2.50, P = 1.22 × 10(-9)). In addition, we identified a novel association signal on chr2q14.3 coming from a rare haplotype (OR = 22.58, P = 1.21 × 10(-10)) and replicated it in a secondary cohort of 306 Ashkenazi PD cases and 2583 controls. Our results highlight the power of our haplotype association method, particularly useful in studies of founder populations, and reaffirm the benefits of studying complex diseases in Ashkenazi Jewish cohorts.

Expanded genetic screening panel for the Ashkenazi Jewish population.

PURPOSE: Carrier screening programs that identify the presence of known mutations have been effective for reducing the incidence of autosomal recessive conditions in the Ashkenazi Jewish (AJ) population and other populations. Yet, these programs have not realized their full potential. Furthermore, many known autosomal recessive and dominant conditions are not screened for and the molecular basis of other conditions for which screening might be offered is unknown. METHODS: Through literature review and annotation of full sequenced genomes from healthy individuals, we expanded the list of mutations. Mutations were identified in a sample of 128 fully sequenced AJ genomes that were filtered through clinical databases and curated manually for clinical validity and utility using the American College of Medical Genetics and Genomics scoring (ACMG) system. Other known mutations were identified through literature review. RESULTS: A panel of 163 mutations was identified for 76 autosomal recessive, 24 autosomal dominant, and 3 X-linked disorders. CONCLUSION: Screening for a broader range of disorders not only could further reduce the incidence of autosomal recessive disorders but also could offer the benefits of early or presymptomatic diagnosis.Genet Med 18 5, 522-528.

Length distributions of identity by descent reveal fine-scale demographic history.

Data-driven studies of identity by descent (IBD) were recently enabled by high-resolution genomic data from large cohorts and scalable algorithms for IBD detection. Yet, haplotype sharing currently represents an underutilized source of information for population-genetics research. We present analytical results on the relationship between haplotype sharing across purportedly unrelated individuals and a population's demographic history. We express the distribution of IBD sharing across pairs of individuals for segments of arbitrary length as a function of the population's demography, and we derive an inference procedure to reconstruct such demographic history. The accuracy of the proposed reconstruction methodology was extensively tested on simulated data. We applied this methodology to two densely typed data sets: 500 Ashkenazi Jewish (AJ) individuals and 56 Kenyan Maasai (MKK) individuals (HapMap 3 data set). Reconstructing the demographic history of the AJ cohort, we recovered two subsequent population expansions, separated by a severe founder event, consistent with previous analysis of lower-throughput genetic data and historical accounts of AJ history. In the MKK cohort, high levels of cryptic relatedness were detected. The spectrum of IBD sharing is consistent with a demographic model in which several small-sized demes intermix through high migration rates and result in enrichment of shared long-range haplotypes. This scenario of historically structured demographies might explain the unexpected abundance of runs of homozygosity within several populations.

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci.

Crohn's disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD-susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2-4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻6). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻8; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻9; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻8; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻8; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻9; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.

Genome-wide association study of schizophrenia in Ashkenazi Jews.

Schizophrenia is a common, clinically heterogeneous disorder associated with lifelong morbidity and early mortality. Several genetic variants associated with schizophrenia have been identified, but the majority of the heritability remains unknown. In this study, we report on a case-control sample of Ashkenazi Jews (AJ), a founder population that may provide additional insights into genetic etiology of schizophrenia. We performed a genome-wide association analysis (GWAS) of 592 cases and 505 controls of AJ ancestry ascertained in the US. Subsequently, we performed a meta-analysis with an Israeli AJ sample of 913 cases and 1640 controls, followed by a meta-analysis and polygenic risk scoring using summary results from Psychiatric GWAS Consortium 2 schizophrenia study. The U.S. AJ sample showed strong evidence of polygenic inheritance (pseudo-R(2) ∼9.7%) and a SNP-heritability estimate of 0.39 (P = 0.00046). We found no genome-wide significant associations in the U.S. sample or in the combined US/Israeli AJ meta-analysis of 1505 cases and 2145 controls. The strongest AJ specific associations (P-values in 10(-6) -10(-7) range) were in the 22q 11.2 deletion region and included the genes TBX1, GLN1, and COMT. Supportive evidence (meta P < 1 × 10(-4) ) was also found for several previously identified genome-wide significant findings, including the HLA region, CNTN4, IMMP2L, and GRIN2A. The meta-analysis of the U.S. sample with the PGC2 results provided initial genome-wide significant evidence for six new loci. Among the novel potential susceptibility genes is PEPD, a gene involved in proline metabolism, which is associated with a Mendelian disorder characterized by developmental delay and cognitive deficits.

Independent evidence for an association between general cognitive ability and a genetic locus for educational attainment.

Cognitive deficits and reduced educational achievement are common in psychiatric illness; understanding the genetic basis of cognitive and educational deficits may be informative about the etiology of psychiatric disorders. A recent, large genome-wide association study (GWAS) reported a genome-wide significant locus for years of education, which subsequently demonstrated association to general cognitive ability ("g") in overlapping cohorts. The current study was designed to test whether GWAS hits for educational attainment are involved in general cognitive ability in an independent, large-scale collection of cohorts. Using cohorts in the Cognitive Genomics Consortium (COGENT; up to 20,495 healthy individuals), we examined the relationship between g and variants associated with educational attainment. We next conducted meta-analyses with 24,189 individuals with neurocognitive data from the educational attainment studies, and then with 53,188 largely independent individuals from a recent GWAS of cognition. A SNP (rs1906252) located at chromosome 6q16.1, previously associated with years of schooling, was significantly associated with g (P = 1.47 × 10(-4) ) in COGENT. The first joint analysis of 43,381 non-overlapping individuals for this a priori-designated locus was strongly significant (P = 4.94 × 10(-7) ), and the second joint analysis of 68,159 non-overlapping individuals was even more robust (P = 1.65 × 10(-9) ). These results provide independent replication, in a large-scale dataset, of a genetic locus associated with cognitive function and education. As sample sizes grow, cognitive GWAS will identify increasing numbers of associated loci, as has been accomplished in other polygenic quantitative traits, which may be relevant to psychiatric illness.

The architecture of long-range haplotypes shared within and across populations.

Homologous long segments along the genomes of close or remote relatives that are identical by descent (IBD) from a common ancestor provide clues for recent events in human genetics. We set out to extensively map such IBD segments in large cohorts and investigate their distribution within and across different populations. We report analysis of several data sets, demonstrating that IBD is more common than expected by naïve models of population genetics. We show that the frequency of IBD pairs is population dependent and can be used to cluster individuals into populations, detect a homogeneous subpopulation within a larger cohort, and infer bottleneck events in such a subpopulation. Specifically, we show that Ashkenazi Jewish individuals are all connected through transitive remote family ties evident by sharing of 50 cM IBD to a publicly available data set of less than 400 individuals. We further expose regions where long-range haplotypes are shared significantly more often than elsewhere in the genome, observed across multiple populations, and enriched for common long structural variation. These are inconsistent with recent relatedness and suggest ancient common ancestry, with limited recombination between haplotypes.

Susceptibility to chronic pain following nerve injury is genetically affected by CACNG2.

Chronic neuropathic pain is affected by specifics of the precipitating neural pathology, psychosocial factors, and by genetic predisposition. Little is known about the identity of predisposing genes. Using an integrative approach, we discovered that CACNG2 significantly affects susceptibility to chronic pain following nerve injury. CACNG2 encodes for stargazin, a protein intimately involved in the trafficking of glutamatergic AMPA receptors. The protein might also be a Ca(2+) channel subunit. CACNG2 has previously been implicated in epilepsy. Initially, using two fine-mapping strategies in a mouse model (recombinant progeny testing [RPT] and recombinant inbred segregation test [RIST]), we mapped a pain-related quantitative trait locus (QTL) (Pain1) into a 4.2-Mb interval on chromosome 15. This interval includes 155 genes. Subsequently, bioinformatics and whole-genome microarray expression analysis were used to narrow the list of candidates and ultimately to pinpoint Cacng2 as a likely candidate. Analysis of stargazer mice, a Cacng2 hypomorphic mutant, provided electrophysiological and behavioral evidence for the gene's functional role in pain processing. Finally, we showed that human CACNG2 polymorphisms are associated with chronic pain in a cohort of cancer patients who underwent breast surgery. Our findings provide novel information on the genetic basis of neuropathic pain and new insights into pain physiology that may ultimately enable better treatments.

The Collaborative Cross, a community resource for the genetic analysis of complex traits.

The goal of the Complex Trait Consortium is to promote the development of resources that can be used to understand, treat and ultimately prevent pervasive human diseases. Existing and proposed mouse resources that are optimized to study the actions of isolated genetic loci on a fixed background are less effective for studying intact polygenic networks and interactions among genes, environments, pathogens and other factors. The Collaborative Cross will provide a common reference panel specifically designed for the integrative analysis of complex systems and will change the way we approach human health and disease.

Genome-wide association identifies a common variant in the reelin gene that increases the risk of schizophrenia only in women.

Sex differences in schizophrenia are well known, but their genetic basis has not been identified. We performed a genome-wide association scan for schizophrenia in an Ashkenazi Jewish population using DNA pooling. We found a female-specific association with rs7341475, a SNP in the fourth intron of the reelin (RELN) gene (p = 2.9 x 10(-5) in women), with a significant gene-sex effect (p = 1.8 x 10(-4)). We studied rs7341475 in four additional populations, totaling 2,274 cases and 4,401 controls. A significant effect was observed only in women, replicating the initial result (p = 2.1 x 10(-3) in women; p = 4.2 x 10(-3) for gene-sex interaction). Based on all populations the estimated relative risk of women carrying the common genotype is 1.58 (p = 8.8 x 10(-7); p = 1.6 x 10(-5) for gene-sex interaction). The female-specific association between RELN and schizophrenia is one of the few examples of a replicated sex-specific genetic association in any disease.

Novel ultra-rare exonic variants identified in a founder population implicate cadherins in schizophrenia.

The identification of rare variants associated with schizophrenia has proven challenging due to genetic heterogeneity, which is reduced in founder populations. In samples from the Ashkenazi Jewish population, we report that schizophrenia cases had a greater frequency of novel missense or loss of function (MisLoF) ultra-rare variants (URVs) compared to controls, and the MisLoF URV burden was inversely correlated with polygenic risk scores in cases. Characterizing 141 "case-only" genes (MisLoF URVs in ≥3 cases with none in controls), the cadherin gene set was associated with schizophrenia. We report a recurrent case mutation in PCDHA3 that results in the formation of cytoplasmic aggregates and failure to engage in homophilic interactions on the plasma membrane in cultured cells. Modeling purifying selection, we demonstrate that deleterious URVs are greatly overrepresented in the Ashkenazi population, yielding enhanced power for association studies. Identification of the cadherin/protocadherin family as risk genes helps specify the synaptic abnormalities central to schizophrenia.

Correlational analysis for identifying genes whose regulation contributes to chronic neuropathic pain.

BACKGROUND: Nerve injury-triggered hyperexcitability in primary sensory neurons is considered a major source of chronic neuropathic pain. The hyperexcitability, in turn, is thought to be related to transcriptional switching in afferent cell somata. Analysis using expression microarrays has revealed that many genes are regulated in the dorsal root ganglion (DRG) following axotomy. But which contribute to pain phenotype versus other nerve injury-evoked processes such as nerve regeneration? Using the L5 spinal nerve ligation model of neuropathy we examined differential changes in gene expression in the L5 (and L4) DRGs in five mouse strains with contrasting susceptibility to neuropathic pain. We sought genes for which the degree of regulation correlates with strain-specific pain phenotype. RESULTS: In an initial experiment six candidate genes previously identified as important in pain physiology were selected for in situ hybridization to DRG sections. Among these, regulation of the Na+ channel alpha subunit Scn11a correlated with levels of spontaneous pain behavior, and regulation of the cool receptor Trpm8 correlated with heat hypersensibility. In a larger scale experiment, mRNA extracted from individual mouse DRGs was processed on Affymetrix whole-genome expression microarrays. Overall, 2552 +/- 477 transcripts were significantly regulated in the axotomized L5DRG 3 days postoperatively. However, in only a small fraction of these was the degree of regulation correlated with pain behavior across strains. Very few genes in the "uninjured" L4DRG showed altered expression (24 +/- 28). CONCLUSION: Correlational analysis based on in situ hybridization provided evidence that differential regulation of Scn11a and Trpm8 contributes to across-strain variability in pain phenotype. This does not, of course, constitute evidence that the others are unrelated to pain. Correlational analysis based on microarray data yielded a larger "look-up table" of genes whose regulation likely contributes to pain variability. While this list is enriched in genes of potential importance for pain physiology, and is relatively free of the bias inherent in the candidate gene approach, additional steps are required to clarify which transcripts on the list are in fact of functional importance.

Association between golli-MBP and schizophrenia in the Jewish Ashkenazi population: are regulatory regions involved?

Multiple studies have reported oligodendrocyte and myelin abnormalities, as well as dysregulation of their related genes, in brains of schizophrenia patients. One of these genes is the myelin-basic-protein (MBP) gene, which encodes two families of proteins: classic-MBPs and golli-MBPs. While the classic-MBPs are predominantly located in the myelin sheaths of the nervous system, the golli proteins are more widely expressed and are found in both the immune and the nervous systems. In the present study we performed a case-control association analysis of golli-MBP in two separate Jewish Ashkenazi cohorts (cohort I: 120 patients, 236 controls; cohort II: 379 patients, 380 controls). In addition we performed an expression analysis of golli-MBP mRNA in post-mortem dorsolateral prefrontal cortex samples of schizophrenia patients, and matched controls. In the first cohort we observed association between six (out of 26 genotyped) single nucleotide polymorphisms (SNPs) and the disease (p<0.05). Of these, three are from one linkage disequilibrium (LD) block which contains a CTCF binding region. Haplotype analysis revealed significant 'risk'/'protective' haplotypes (strongest p=0.005, each) for schizophrenia. The three SNPs (rs12458282, rs2008323, rs721286) were then genotyped in the second cohort. The combined results showed strong effects, both in the single marker and in haplotype analyses (strongest OR 1.77, p=0.0005; OR 1.61, p=0.00001, respectively). Sequencing the CTCF binding region revealed three SNPs in complete LD with the associated haplotypes, located in close proximity to the CTCF binding site. Expression analysis found no significant differences in golli-MBP mRNA levels. These findings suggest that golli-MBP is a possible susceptibility gene for schizophrenia.