Enhancing Psychosis-Spectrum Nosology Through an International Data Sharing Initiative.

The latent structure of schizotypy and psychosis-spectrum symptoms remains poorly understood. Furthermore, molecular genetic substrates are poorly defined, largely due to the substantial resources required to collect rich phenotypic data across diverse populations. Sample sizes of phenotypic studies are often insufficient for advanced structural equation modeling approaches. In the last 50 years, efforts in both psychiatry and psychological science have moved toward (1) a dimensional model of psychopathology (eg, the current Hierarchical Taxonomy of Psychopathology [HiTOP] initiative), (2) an integration of methods and measures across traits and units of analysis (eg, the RDoC initiative), and (3) powerful, impactful study designs maximizing sample size to detect subtle genomic variation relating to complex traits (the Psychiatric Genomics Consortium [PGC]). These movements are important to the future study of the psychosis spectrum, and to resolving heterogeneity with respect to instrument and population. The International Consortium of Schizotypy Research is composed of over 40 laboratories in 12 countries, and to date, members have compiled a body of schizotypy- and psychosis-related phenotype data from more than 30000 individuals. It has become apparent that compiling data into a protected, relational database and crowdsourcing analytic and data science expertise will result in significant enhancement of current research on the structure and biological substrates of the psychosis spectrum. The authors present a data-sharing infrastructure similar to that of the PGC, and a resource-sharing infrastructure similar to that of HiTOP. This report details the rationale and benefits of the phenotypic data collective and presents an open invitation for participation.

Genome-wide gene pathway analysis of psychotic illness symptom dimensions based on a new schizophrenia-specific model of the OPCRIT.

Empirically derived phenotypic measurements have the potential to enhance gene-finding efforts in schizophrenia. Previous research based on factor analyses of symptoms has typically included schizoaffective cases. Deriving factor loadings from analysis of only narrowly defined schizophrenia cases could yield more sensitive factor scores for gene pathway and gene ontology analyses. Using an Irish family sample, this study 1) factor analyzed clinician-rated Operational Criteria Checklist items in cases with schizophrenia only, 2) scored the full sample based on these factor loadings, and 3) implemented genome-wide association, gene-based, and gene-pathway analysis of these SCZ-based symptom factors (final N=507). Three factors emerged from the analysis of the schizophrenia cases: a manic, a depressive, and a positive symptom factor. In gene-based analyses of these factors, multiple genes had q<0.01. Of particular interest are findings for PTPRG and WBP1L, both of which were previously implicated by the Psychiatric Genomics Consortium study of SCZ; results from this study suggest that variants in these genes might also act as modifiers of SCZ symptoms. Gene pathway analyses of the first factor indicated over-representation of glutamatergic transmission, GABA-A receptor, and cyclic GMP pathways. Results suggest that these pathways may have differential influence on affective symptom presentation in schizophrenia.

Molecular validation of the schizophrenia spectrum.

BACKGROUND: Early descriptive work and controlled family and adoption studies support the hypothesis that a range of personality and nonschizophrenic psychotic disorders aggregate in families of schizophrenic probands. Can we validate, using molecular polygene scores from genome-wide association studies (GWAS), this schizophrenia spectrum? METHODS: The predictive value of polygenic findings reported by the Psychiatric GWAS Consortium (PGC) was applied to 4 groups of relatives from the Irish Study of High-Density Schizophrenia Families (ISHDSF; N = 836) differing on their assignment within the schizophrenia spectrum. Genome-wide single nucleotide polymorphism data for affected and unaffected relatives were used to construct per-individual polygene risk scores based on the PGC stage-I results. We compared mean polygene scores in the ISHDSF with mean scores in ethnically matched population controls (N = 929). RESULTS: The schizophrenia polygene score differed significantly across diagnostic categories and was highest in those with narrow schizophrenia spectrum, lowest in those with no psychiatric illness, and in-between in those classified in the intermediate, broad, and very broad schizophrenia spectrum. Relatives of all of these groups of affected subjects, including those with no diagnosis, had schizophrenia polygene scores significantly higher than the control sample. CONCLUSIONS: In the relatives of high-density families, the observed pattern of enrichment of molecular indices of schizophrenia risk suggests an underlying, continuous liability distribution and validates, using aggregate common risk alleles, a genetic basis for the schizophrenia spectrum disorders. In addition, as predicted by genetic theory, nonpsychotic members of multiply-affected schizophrenia families are significantly enriched for replicated, polygenic risk variants compared with the general population.

No association of dysbindin with symptom factors of schizophrenia in an Irish case-control sample.

Robust associations between the dysbindin gene (DTNBP1) and schizophrenia have been demonstrated in many but not all samples, and evidence that this gene particularly predisposes to negative symptoms in this illness has been presented. The current study sought to replicate the previously reported negative symptom associations in an Irish case-control sample. Association between dysbindin and schizophrenia has been established in this cohort, and a factor analysis of the assessed symptoms yielded three factors, Positive, Negative, and Schneiderian. The sequential addition method was applied using UNPHASED to assess the relationship between these symptom factors and the high-risk haplotype. No associations were detected for any of the symptom factors indicating that the dysbindin risk haplotype does not predispose to a particular group of symptoms in this sample. Several possibilities, such as differing risk haplotypes, may explain this finding.

The dystrobrevin binding protein 1 (DTNBP1) gene is associated with schizophrenia in the Irish Case Control Study of Schizophrenia (ICCSS) sample.

BACKGROUND: DTNBP1 is associated with schizophrenia in many studies, but the associated alleles and haplotypes vary between samples. METHOD: We assessed nine single nucleotide polymorphisms (SNPs) in this gene for association with schizophrenia in a new sample of 1021 cases and 626 controls from Ireland. RESULTS: Four SNPs give evidence of association (0.000018

Apoptotic engulfment pathway and schizophrenia.

BACKGROUND: Apoptosis has been speculated to be involved in schizophrenia. In a previously study, we reported the association of the MEGF10 gene with the disease. In this study, we followed the apoptotic engulfment pathway involving the MEGF10, GULP1, ABCA1 and ABCA7 genes and tested their association with the disease. METHODOLOGY/PRINCIPAL FINDINGS: Ten, eleven and five SNPs were genotyped in the GULP1, ABCA1 and ABCA7 genes respectively for the ISHDSF and ICCSS samples. In all 3 genes, we observed nominally significant associations. Rs2004888 at GULP1 was significant in both ISHDSF and ICCSS samples (p = 0.0083 and 0.0437 respectively). We sought replication in independent samples for this marker and found highly significant association (p = 0.0003) in 3 Caucasian replication samples. But it was not significant in the 2 Chinese replication samples. In addition, we found a significant 2-marker (rs2242436 * rs3858075) interaction between the ABCA1 and ABCA7 genes in the ISHDSF sample (p = 0.0022) and a 3-marker interaction (rs246896 * rs4522565 * rs3858075) amongst the MEGF10, GULP1 and ABCA1 genes in the ICCSS sample (p = 0.0120). Rs3858075 in the ABCA1 gene was involved in both 2- and 3-marker interactions in the two samples. CONCLUSIONS/SIGNIFICANCE: From these data, we concluded that the GULP1 gene and the apoptotic engulfment pathway are involved in schizophrenia in subjects of European ancestry and multiple genes in the pathway may interactively increase the risks to the disease.

Polymorphisms in SLC6A4, PAH, GABRB3, and MAOB and modification of psychotic disorder features.

We tested four genes [phenylalanine hydroxylase (PAH), the serotonin transporter (SLC6A4), monoamine oxidase B (MAOB), and the gamma-aminobutyric acid A receptor beta-3 subunit (GABRB3)] for their impact on five schizophrenia symptom factors: delusions, hallucinations, mania, depression, and negative symptoms. In a 90 family subset of the Irish Study of High Density Schizophrenia Families, the PAH 232 bp microsatellite allele demonstrated significant association with the delusions factor using both QTDT (F=8.0, p=.031) and QPDTPHASE (chi-square=12.54, p=.028). Also, a significant association between the GABRB3 191 bp allele and the hallucinations factor was detected using QPDTPHASE (chi-square=15.51, p=.030), but not QTDT (chi-square=2.07, p=.560).

The trace amine associated receptor (TAAR6) gene is not associated with schizophrenia in the Irish Case-Control Study of Schizophrenia (ICCSS) sample.

To replicate previous association between TAAR6 and schizophrenia, including our own finding in the Irish Study of High Density Schizophrenia Families (ISHDSF) sample, we genotyped 12 single nucleotide polymorphisms (SNPs) in the Irish Case-Control Study of Schizophrenia (ICCSS) sample. Only rs9389020 provided nominal evidence for association (p<0.0228), which did not withstand the permutation testing (p<0.2196). The combined odds ratio from ISHDSF and ICCSS samples [OR (95%CI)=1.0564 (1.0078-1.1074); p=0.02], while nominally significant, did not survive correction for multiple testing. Here we demonstrate that TAAR6 is not associated with schizophrenia in the ICCSS sample.

Genetic variation in the serotonin 2A receptor and suicidal ideation in a sample of 270 Irish high-density schizophrenia families.

Genetic variation in the serotonin 2A receptor (HTR2A) has been associated with both schizophrenia and suicidal behavior. Our sample comprised 270 Irish high-density schizophrenia families (n = 1,408 subjects, including 755 with psychotic illness). Diagnoses were generated using a modified SCID. All patients who had at least one episode of psychosis were rated on the Operation Criteria Checklist for Psychotic Illness (OPCRIT). Lifetime history of suicidal ideation was determined from medical records and psychiatric interviews and was scored in the OPCRIT. Twelve SNPs were selected for study. Ten of these were tagSNPs derived from HapMap data, along with His452Tyr and T102C. We tested for association with psychotic illness as a whole, as well as stratified by the presence of suicidal ideation, using FBAT and PDTPHASE. Single-marker as well as haplotype-based tests using a "sliding window" approach were performed. We observed several 2, 3, and 4 marker haplotypes near the 3' end of the gene that were over-transmitted to psychotic subjects (0.02

Interaction between interleukin 3 and dystrobrevin-binding protein 1 in schizophrenia.

Schizophrenia is a common psychotic mental disorder that is believed to result from the effects of multiple genetic and environmental factors. In this study, we explored gene-gene interactions and main effects in both case-control (657 cases and 411 controls) and family-based (273 families, 1,350 subjects) datasets of English or Irish ancestry. Fifty three markers in 8 genes were genotyped in the family sample and 44 markers in 7 genes were genotyped in the case-control sample. The Multifactor Dimensionality Reduction Pedigree Disequilibrium Test (MDR-PDT) was used to examine epistasis in the family dataset and a 3-locus model was identified (permuted p=0.003). The 3-locus model involved the IL3 (rs2069803), RGS4 (rs2661319), and DTNBP1 (rs2619539) genes. We used MDR to analyze the case-control dataset containing the same markers typed in the RGS4, IL3 and DTNBP1 genes and found evidence of a joint effect between IL3 (rs31400) and DTNBP1 (rs760761) (cross-validation consistency 4/5, balanced prediction accuracy=56.84%, p=0.019). While this is not a direct replication, the results obtained from both the family and case-control samples collectively suggest that IL3 and DTNBP1 are likely to interact and jointly contribute to increase risk for schizophrenia. We also observed a significant main effect in DTNBP1, which survived correction for multiple comparisons, and numerous nominally significant effects in several genes.

FBXL21 association with schizophrenia in Irish family and case-control samples.

FBXL21 gene encodes an F-box containing protein functioning in the SCF ubiquitin ligase complex. The role of the F-box protein is to recruit proteins designated for degradation to the ligase complex so they would be ubiquitinated. Using both family and case-control samples, we found consistent associations in and around FBXL21 gene. In the family sample (Irish study of high density schizophrenia families, ISHDSF, 1,350 subjects from 273 families), a minimal PDT P-value of 0.0011 was observed at rs31555. In the case-control sample (Irish case-control study of schizophrenia, ICCSS, 814 cases and 625 controls), significant associations were observed at two markers (rs1859427 P = 0.0197, and rs6861170 P = 0.0197). In haplotype analyses, haplotype 1-1 (C-T) of rs1859427-rs6861170 was overtransmitted in the ISHDSF (P = 0.0437) and was overrepresented in the ICCSS (P = 0.0177). For both samples, the associated alleles and haplotypes were identical. These data suggested that FBXL21 may be associated with schizophrenia in the Irish samples.

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.

MEGF10 association with schizophrenia.

BACKGROUND: The 5q21-31 region has been implicated as harboring risk genes for schizophrenia in many linkage studies. In our previous report of stepwise fine mapping of this region, the MEGF10 gene was one of the genes showing consistent associations in our screening subsample. In this study, we carried out independent replication and expression studies of the MEGF10 gene. METHODS: Association studies with 8 SNPs in the MEGF10 gene were performed in the Irish case-control study of schizophrenia (ICCSS) sample (652 case and 617 control subjects). The expression of MEGF10 was also compared between healthy control subjects and schizophrenia patients using postmortem brain cDNA libraries. RESULTS: In the ICCSS sample, associations with the disease were found in the same risk alleles and haplotypes as that observed in our fine-mapping studies. The major allele (A) of rs27388 was overrepresented in affected individuals (p = .0169), which remained significant after correction for multiple testing. In expression studies, MEGF10 had higher expression levels in the affected than the unaffected (p = .015). Schizophrenia patients with a 1/1 genotype at rs27388 had higher expressions than those patients with 1/2 and 2/2 genotypes (p = .0008). CONCLUSIONS: Evidence from both association and expression studies suggests that MEGF10 is likely associated with schizophrenia. The major allele and 1/1 genotype at rs27388 impose higher risks for the disease.

A genome-wide scan for modifier loci in schizophrenia.

The purpose of this study was to detect genetic loci that influence clinical features of, but not necessarily susceptibility to, psychotic illness. In the Irish Study of High-Density Schizophrenia Families (n = 270 families, n = 1,408 individuals), subjects with non-affective psychosis were rated using the Operational Criteria Checklist for Psychotic Illness. Factor analysis identified hallucinations, delusions, and negative, manic, and depressive symptom factors. We performed autosomal genome-wide multipoint non-parametric quantitative trait locus linkage analysis, in affected individuals only, using these five factors, as well as age at onset, and course of illness. Determination of empirical significance and correction for multiple testing was implemented using 200 simulated genome scans. We also tested for pleiotropic loci by examining the sums of -log(10)'s of the empirical P values of multiple traits in selected regions. LODs of 2.42 and 2.35 were obtained near D9S934 (9q33.1) and D14S587 (14q24.2), respectively, for course of illness, and of 2.26 between D6S1040-D6S2420 (6q23.1-25.1) and age at onset. No other regions met criteria for suggestive linkage to any one trait. No loci were significant after correction for multiple testing. On 6q, however, the joint linkage of age of onset, course, delusions, and depressive symptoms resulted in a genome-wide P = 0.06. We conclude that genes located near 9q33.1 and 14q24.2 may modify the clinical course and severity of schizophrenia. A gene in 6q may affect several clinical features of illness.

Association between the 5q31.1 gene neurogenin1 and schizophrenia.

Multiple lines of evidence suggest that schizophrenia results from aberrant neurodevelopment. The neurogenin1 gene (neurog1) consists of a single 1,666 bp exon that encodes a basic helix-loop-helix (bHLH) transcription factor that causes neuronal differentiation and induces cortical and glutamatergic differentiation programs. Because of its function and its location in 5q31.1, which has been linked to schizophrenia in multiple samples, we tested it for association with the disorder. We sequenced neurog1 in 25 affected subjects from the Irish Study of High-Density Schizophrenia Families. We observed a 5'-UTR SNP at position -60, already present in databases as rs8192558, and tested it along with rs2344485, rs8192559, and rs2344484. Narrow, intermediate, and broad diagnostic definitions were used. The major alleles of rs8192558 and rs2344484 were over-transmitted to affected subjects using both Pedigree Disequilibrium Test (PDT) (0.01 < or = P < or = 0.06) and FBAT (0.02 < or = P < or = 0.07). A haplotype consisting of the major alleles of all four SNPs was significantly over-transmitted in FBAT to the broad definition (P = 0.049), with trend significance to the narrow and intermediate definitions, and with trend significance in PDT. In confirmatory tests using 657 cases and 411 controls, this haplotype was slightly but not significantly over-represented in cases (81% vs. 77%, P = 0.21). These results, along with a priori evidence for the involvement of neurog1 in neurodevelopment, suggest that variants in neurog1 might have a small effect on susceptibility to schizophrenia. This gene should be tested in additional and larger samples.

Haplotypes spanning SPEC2, PDZ-GEF2 and ACSL6 genes are associated with schizophrenia.

Chromosome 5q22-33 is a region where studies have repeatedly found evidence for linkage to schizophrenia. In this report, we took a stepwise approach to systematically map this region in the Irish Study of High Density Schizophrenia Families (ISHDSF, 267 families, 1337 subjects) sample. We typed 289 SNPs in the critical interval of 8 million basepairs and found a 758 kb interval coding for the SPEC2/PDZ-GEF2/ACSL6 genes to be associated with the disease. Using sex and genotype-conditioned transmission disequilibrium test analyses, we found that 19 of the 24 typed markers were associated with the disease and the associations were sex-specific. We replicated these findings with an Irish case-control sample (657 cases and 414 controls), an Irish parent-proband trio sample (187 families, 564 subjects), a German nuclear family sample (211 families, 751 subjects) and a Pittsburgh nuclear family sample (247 families, 729 subjects). In all four samples, we replicated the sex-specific associations at the levels of both individual markers and haplotypes using sex- and genotype-conditioned analyses. Three risk haplotypes were identified in the five samples, and each haplotype was found in at least two samples. Consistent with the discovery of multiple estrogen-response elements in this region, our data showed that the impact of these haplotypes on risk for schizophrenia differed in males and females. From these data, we concluded that haplotypes underlying the SPEC2/PDZ-GEF2/ACSL6 region are associated with schizophrenia. However, due to the extended high LD in this region, we were unable to distinguish whether the association signals came from one or more of these genes.

Relationship between a high-risk haplotype in the DTNBP1 (dysbindin) gene and clinical features of schizophrenia.

OBJECTIVE: The purpose of this study was to determine whether a haplotype in the dystrobrevin binding protein 1 (DTNBP1) gene previously associated with schizophrenia not only increases the susceptibility to psychotic illness but also to a more or less clinically specific form of psychotic illness. METHOD: In the Irish Study of High-Density Schizophrenia Families, subjects with psychotic illness (N=755) were given lifetime ratings of clinical features according to the Operational Criteria Checklist for Psychotic Illness. Exploratory and confirmatory factor analyses were used to extract five factors-hallucinations, delusions, negative, manic, and depressive symptoms-and to create factor-derived scores. The family-based transmission disequilibrium test operationalized in the program TRANSMIT was used to determine whether a high-risk haplotype in the DTNBP1 gene was overtransmitted to subjects in the upper 20th and 40th percentiles for each factor score. These results were compared to baseline overtransmission by examining the empirical distribution of chi-square statistics in groups of 5,000 replicates in which 20% and 40% of ill subjects were randomly selected. This analysis was done for both narrow and broad definitions of psychotic illness. RESULTS: Subjects in the upper 40th percentile for the negative symptom factor--in both the narrowly (p=0.004) and broadly (p=0.01) defined illness groups--were more likely to inherit the high-risk haplotype than would be expected by chance. No other significant relationships between clinical features and high-risk haplotype transmission were observed. CONCLUSIONS: The etiologically relevant variation in DTNBP1, which is in presumptive linkage disequilibrium with the high-risk haplotype, may predispose individuals to a form of psychotic illness associated with high levels of negative symptoms. This finding supports previous evidence suggesting that genetic factors influence the clinical heterogeneity of schizophrenia.

Evaluation of genetic substructure in the Irish Study of High-Density Schizophrenia Families.

The presence of genetic substructure has the potential to diminish the chances of detecting a linkage signal. Using a Markov chain Monte Carlo procedure developed by Pritchard and colleagues and implemented in the program STRUCTURE, we evaluated the evidence for genetic substructure using genotypes from 37 microsatellite markers in affected individuals selected at random from 263 multiplex families in the Irish Study of High-Density Schizophrenia Families. We found no evidence for the presence of genetic substructure in this sample.

Regulator of G-protein signaling 4 (RGS4) gene is associated with schizophrenia in Irish high density families.

The regulator of the G-protein signaling 4 (RGS4) gene was shown to have a different expression pattern in schizophrenia patients in a microarray study. A family-based study subsequently implicated the association of this gene with schizophrenia. We replicated the study with our sample from the Irish Study of High Density Schizophrenia Families (ISHDSF). Single marker transmission disequilibrium tests (TDT) for the four core SNPs showed modest association for SNP 18 (using a narrow diagnostic approach with FBAT P = 0.044; with PDT P = 0.0073) and a trend for SNP 4 (with FBAT P = 0.1098; with PDT P = 0.0249). For SNP 1 and 7, alleles overtransmitted to affected subjects were the same as previously reported. Haplotype analyses suggested that haplotype G-G-G for SNP1-4-18, which is the most abundant haplotype (42.3%) in the Irish families, was associated with the disease (narrow diagnosis, FBAT P = 0.0061, PDT P = 0.0498). This was the same haplotype implicated in the original study. While P values were not corrected for multiple testing because of the clear prior hypothesis, these results could be interpreted as supporting evidence for the association between RGS4 and schizophrenia.

Genome scan meta-analysis of schizophrenia and bipolar disorder, part II: Schizophrenia.

Schizophrenia is a common disorder with high heritability and a 10-fold increase in risk to siblings of probands. Replication has been inconsistent for reports of significant genetic linkage. To assess evidence for linkage across studies, rank-based genome scan meta-analysis (GSMA) was applied to data from 20 schizophrenia genome scans. Each marker for each scan was assigned to 1 of 120 30-cM bins, with the bins ranked by linkage scores (1 = most significant) and the ranks averaged across studies (R(avg)) and then weighted for sample size (N(sqrt)[affected casess]). A permutation test was used to compute the probability of observing, by chance, each bin's average rank (P(AvgRnk)) or of observing it for a bin with the same place (first, second, etc.) in the order of average ranks in each permutation (P(ord)). The GSMA produced significant genomewide evidence for linkage on chromosome 2q (PAvgRnk<.000417). Two aggregate criteria for linkage were also met (clusters of nominally significant P values that did not occur in 1,000 replicates of the entire data set with no linkage present): 12 consecutive bins with both P(AvgRnk) and P(ord)<.05, including regions of chromosomes 5q, 3p, 11q, 6p, 1q, 22q, 8p, 20q, and 14p, and 19 consecutive bins with P(ord)<.05, additionally including regions of chromosomes 16q, 18q, 10p, 15q, 6q, and 17q. There is greater consistency of linkage results across studies than has been previously recognized. The results suggest that some or all of these regions contain loci that increase susceptibility to schizophrenia in diverse populations.