Polygenic risk for severe psychopathology among Europeans is associated with major depressive disorder in Han Chinese women.

BACKGROUND: Previous studies have demonstrated that several major psychiatric disorders are influenced by shared genetic factors. This shared liability may influence clinical features of a given disorder (e.g. severity, age at onset). However, findings have largely been limited to European samples; little is known about the consistency of shared genetic liability across ethnicities. METHOD: The relationship between polygenic risk for several major psychiatric diagnoses and major depressive disorder (MDD) was examined in a sample of unrelated Han Chinese women. Polygenic risk scores (PRSs) were generated using European discovery samples and tested in the China, Oxford, and VCU Experimental Research on Genetic Epidemiology [CONVERGE (maximum N = 10 502)], a sample ascertained for recurrent MDD. Genetic correlations between discovery phenotypes and MDD were also assessed. In addition, within-case characteristics were examined. RESULTS: European-based polygenic risk for several major psychiatric disorder phenotypes was significantly associated with the MDD case status in CONVERGE. Risk for clinically significant indicators (neuroticism and subjective well-being) was also associated with case-control status. The variance accounted for by PRS for both psychopathology and for well-being was similar to estimates reported for within-ethnicity comparisons in European samples. However, European-based PRS were largely unassociated with CONVERGE family history, clinical characteristics, or comorbidity. CONCLUSIONS: The shared genetic liability across severe forms of psychopathology is largely consistent across European and Han Chinese ethnicities, with little attenuation of genetic signal relative to within-ethnicity analyses. The overall absence of associations between PRS for other disorders and within-MDD variation suggests that clinical characteristics of MDD may arise due to contributions from ethnicity-specific factors and/or pathoplasticity.

A simple yet accurate correction for winner's curse can predict signals discovered in much larger genome scans.

MOTIVATION: For genetic studies, statistically significant variants explain far less trait variance than 'sub-threshold' association signals. To dimension follow-up studies, researchers need to accurately estimate 'true' effect sizes at each SNP, e.g. the true mean of odds ratios (ORs)/regression coefficients (RRs) or Z-score noncentralities. Naïve estimates of effect sizes incur winner's curse biases, which are reduced only by laborious winner's curse adjustments (WCAs). Given that Z-scores estimates can be theoretically translated on other scales, we propose a simple method to compute WCA for Z-scores, i.e. their true means/noncentralities. RESULTS: WCA of Z-scores shrinks these towards zero while, on P-value scale, multiple testing adjustment (MTA) shrinks P-values toward one, which corresponds to the zero Z-score value. Thus, WCA on Z-scores scale is a proxy for MTA on P-value scale. Therefore, to estimate Z-score noncentralities for all SNPs in genome scans, we propose F: DR I: nverse Q: uantile T: ransformation (FIQT). It (i) performs the simpler MTA of P-values using FDR and (ii) obtains noncentralities by back-transforming MTA P-values on Z-score scale. When compared to competitors, realistic simulations suggest that FIQT is more (i) accurate and (ii) computationally efficient by orders of magnitude. Practical application of FIQT to Psychiatric Genetic Consortium schizophrenia cohort predicts a non-trivial fraction of sub-threshold signals which become significant in much larger supersamples. CONCLUSIONS: FIQT is a simple, yet accurate, WCA method for Z-scores (and ORs/RRs, via simple transformations). AVAILABILITY AND IMPLEMENTATION: A 10 lines R function implementation is available at https://github.com/bacanusa/FIQT CONTACT: sabacanu@vcu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

JEPEGMIX: gene-level joint analysis of functional SNPs in cosmopolitan cohorts.

MOTIVATION: To increase detection power, gene level analysis methods are used to aggregate weak signals. To greatly increase computational efficiency, most methods use as input summary statistics from genome-wide association studies (GWAS). Subsequently, gene statistics are constructed using linkage disequilibrium (LD) patterns from a relevant reference panel. However, all methods, including our own Joint Effect on Phenotype of eQTL/functional single nucleotide polymorphisms (SNPs) associated with a Gene (JEPEG), assume homogeneous panels, e.g. European. However, this renders these tools unsuitable for the analysis of large cosmopolitan cohorts. RESULTS: We propose a JEPEG extension, JEPEGMIX, which similar to one of our software tools, Direct Imputation of summary STatistics of unmeasured SNPs from MIXed ethnicity cohorts, is capable of estimating accurate LD patterns for cosmopolitan cohorts. JEPEGMIX uses this accurate LD estimates to (i) impute the summary statistics at unmeasured functional variants and (ii) test for the joint effect of all measured and imputed functional variants which are associated with a gene. We illustrate the performance of our tool by analyzing the GWAS meta-analysis summary statistics from the multi-ethnic Psychiatric Genomics Consortium Schizophrenia stage 2 cohort. This practical application supports the immune system being one of the main drivers of the process leading to schizophrenia. AVAILABILITY AND IMPLEMENTATION: Software, annotation database and examples are available at http://dleelab.github.io/jepegmix/. CONTACT: donghyung.lee@vcuhealth.org SUPPLEMENTARY INFORMATION: Supplementary material is available at Bioinformatics online.

Meta-analysis of genome-wide association studies of anxiety disorders.

Anxiety disorders (ADs), namely generalized AD, panic disorder and phobias, are common, etiologically complex conditions with a partially genetic basis. Despite differing on diagnostic definitions based on clinical presentation, ADs likely represent various expressions of an underlying common diathesis of abnormal regulation of basic threat-response systems. We conducted genome-wide association analyses in nine samples of European ancestry from seven large, independent studies. To identify genetic variants contributing to genetic susceptibility shared across interview-generated DSM-based ADs, we applied two phenotypic approaches: (1) comparisons between categorical AD cases and supernormal controls, and (2) quantitative phenotypic factor scores (FS) derived from a multivariate analysis combining information across the clinical phenotypes. We used logistic and linear regression, respectively, to analyze the association between these phenotypes and genome-wide single nucleotide polymorphisms. Meta-analysis for each phenotype combined results across the nine samples for over 18 000 unrelated individuals. Each meta-analysis identified a different genome-wide significant region, with the following markers showing the strongest association: for case-control contrasts, rs1709393 located in an uncharacterized non-coding RNA locus on chromosomal band 3q12.3 (P=1.65 × 10(-8)); for FS, rs1067327 within CAMKMT encoding the calmodulin-lysine N-methyltransferase on chromosomal band 2p21 (P=2.86 × 10(-9)). Independent replication and further exploration of these findings are needed to more fully understand the role of these variants in risk and expression of ADs.

JEPEG: a summary statistics based tool for gene-level joint testing of functional variants.

MOTIVATION: Gene expression is influenced by variants commonly known as expression quantitative trait loci (eQTL). On the basis of this fact, researchers proposed to use eQTL/functional information univariately for prioritizing single nucleotide polymorphisms (SNPs) signals from genome-wide association studies (GWAS). However, most genes are influenced by multiple eQTLs which, thus, jointly affect any downstream phenotype. Therefore, when compared with the univariate prioritization approach, a joint modeling of eQTL action on phenotypes has the potential to substantially increase signal detection power. Nonetheless, a joint eQTL analysis is impeded by (i) not measuring all eQTLs in a gene and/or (ii) lack of access to individual genotypes. RESULTS: We propose joint effect on phenotype of eQTL/functional SNPs associated with a gene (JEPEG), a novel software tool which uses only GWAS summary statistics to (i) impute the summary statistics at unmeasured eQTLs and (ii) test for the joint effect of all measured and imputed eQTLs in a gene. We illustrate the behavior/performance of the developed tool by analysing the GWAS meta-analysis summary statistics from the Psychiatric Genomics Consortium Stage 1 and the Genetic Consortium for Anorexia Nervosa. CONCLUSIONS: Applied analyses results suggest that JEPEG complements commonly used univariate GWAS tools by: (i) increasing signal detection power via uncovering (a) novel genes or (b) known associated genes in smaller cohorts and (ii) assisting in fine-mapping of challenging regions, e.g. major histocompatibility complex for schizophrenia. AVAILABILITY AND IMPLEMENTATION: JEPEG, its associated database of eQTL SNPs and usage examples are publicly available at http://code.google.com/p/jepeg/. CONTACT: dlee4@vcu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DISTMIX: direct imputation of summary statistics for unmeasured SNPs from mixed ethnicity cohorts.

MOTIVATION: To increase the signal resolution for large-scale meta-analyses of genome-wide association studies, genotypes at unmeasured single nucleotide polymorphisms (SNPs) are commonly imputed using large multi-ethnic reference panels. However, the ever increasing size and ethnic diversity of both reference panels and cohorts makes genotype imputation computationally challenging for moderately sized computer clusters. Moreover, genotype imputation requires subject-level genetic data, which unlike summary statistics provided by virtually all studies, is not publicly available. While there are much less demanding methods which avoid the genotype imputation step by directly imputing SNP statistics, e.g. Directly Imputing summary STatistics (DIST) proposed by our group, their implicit assumptions make them applicable only to ethnically homogeneous cohorts. RESULTS: To decrease computational and access requirements for the analysis of cosmopolitan cohorts, we propose DISTMIX, which extends DIST capabilities to the analysis of mixed ethnicity cohorts. The method uses a relevant reference panel to directly impute unmeasured SNP statistics based only on statistics at measured SNPs and estimated/user-specified ethnic proportions. Simulations show that the proposed method adequately controls the Type I error rates. The 1000 Genomes panel imputation of summary statistics from the ethnically diverse Psychiatric Genetic Consortium Schizophrenia Phase 2 suggests that, when compared to genotype imputation methods, DISTMIX offers comparable imputation accuracy for only a fraction of computational resources. AVAILABILITY AND IMPLEMENTATION: DISTMIX software, its reference population data, and usage examples are publicly available at http://code.google.com/p/distmix. CONTACT: dlee4@vcu.edu SUPPLEMENTARY INFORMATION: Supplementary Data are available at Bioinformatics online.

DIST: direct imputation of summary statistics for unmeasured SNPs.

MOTIVATION: Genotype imputation methods are used to enhance the resolution of genome-wide association studies, and thus increase the detection rate for genetic signals. Although most studies report all univariate summary statistics, many of them limit the access to subject-level genotypes. Because such an access is required by all genotype imputation methods, it is helpful to develop methods that impute summary statistics without going through the interim step of imputing genotypes. Even when subject-level genotypes are available, due to the substantial computational cost of the typical genotype imputation, there is a need for faster imputation methods. RESULTS: Direct Imputation of summary STatistics (DIST) imputes the summary statistics of untyped variants without first imputing their subject-level genotypes. This is achieved by (i) using the conditional expectation formula for multivariate normal variates and (ii) using the correlation structure from a relevant reference population. When compared with genotype imputation methods, DIST (i) requires only a fraction of their computational resources, (ii) has comparable imputation accuracy for independent subjects and (iii) is readily applicable to the imputation of association statistics coming from large pedigree data. Thus, the proposed application is useful for a fast imputation of summary results for (i) studies of unrelated subjects, which (a) do not provide subject-level genotypes or (b) have a large size and (ii) family association studies. AVAILABILITY AND IMPLEMENTATION: Pre-compiled executables built under commonly used operating systems are publicly available at http://code.google.com/p/dist/. CONTACT: dlee4@vcu.edu .

Statistical properties of single-marker tests for rare variants.

With the dramatic technological developments of genome-wide association single-nucleotide polymorphism (SNP) chips and next generation sequencing, human geneticists now have the ability to assay genetic variation at ever-rarer allele frequencies. To fully understand the impact of these rare variants on common, complex diseases, we must be able to accurately assess their statistical significance. However, it is well established that classical association tests are not appropriate for the analysis of low-frequency variation, giving spurious findings when observed counts are too few. To further our understanding of the asymptotic properties of traditional association tests, we conducted a range of simulations of a typical rare variant (~1%) under the null hypothesis and tested the allelic χ2, Cochran-Armitage trend, Wald, and Fisher's exact tests. We demonstrate that rare variation shows marked deviation from the expected distributional behavior for each test, with fewer minor alleles corresponding to a greater degree of test statistics deflation. The effect becomes more pronounced at progressively smaller α levels. We also show that the Wald test is particularly deflated at α levels consistent with genome-wide association significance, much more so than the other association tests considered. In general, these classical association tests are inappropriate for the analysis of variants for which the minor allele is observed fewer than 80 times, largely irrespective of sample size.

Association study of 83 candidate genes for bipolar disorder in chromosome 6q selected using an evidence-based prioritization algorithm.

BACKGROUND: Prior genome-scans of bipolar disorder have revealed chromosome 6q22 as a promising candidate region. However, linkage disequilibrium (LD) mapping studies have yet to identify replicated susceptibility loci. METHODS: We analyzed 1,422 LD-tagging single nucleotide polymorphisms (SNPs) in 83 genes to test single-marker and locus-wide evidence of association with bipolar disorder in the NIMH Genetics Initiative bipolar pedigrees and the Portuguese Island Collection (PIC) (N = 1,093 in 528 informative pairs). Both studies previously demonstrated significant evidence of linkage to 6q. SNPs were genotyped using an Illumina iSelect genotyping array which employs the Infinium assay. Evidence of single-marker association was assessed using the generalized disequilibrium test (GDT). Empirical estimates of gene-wide significance were obtained by permutation (via 100,000 gene-dropping simulations) of Fisher's combined test of P-values for each locus. RESULTS: No single variant yielded significant experiment-wide evidence of association, for either the combined sample or in each subsample. Our gene-dropping simulations identified nominally significant gene-wide associations with multiple loci, of which NT5DC1 in the NIMH subsample and CCNC in the PIC were the strongest candidates. However, no one gene consistently exceeded empirical significance criteria in both independent samples or survived Bonferroni correction for the number of genes tested. CONCLUSIONS: Using a gene-based approach to family-based association, we identified gene-wide associations with several genes, though no single locus was significantly associated with bipolar disorder in both cohorts. This suggests that chromosome 6q may harbor multiple susceptibility loci or that complex patterns of LD in this region may confound approaches based on common SNPs. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.

Replication of association between schizophrenia and ZNF804A in the Irish Case-Control Study of Schizophrenia sample.

A recent genome-wide association study reported association between schizophrenia and the ZNF804A gene on chromosome 2q32.1. We attempted to replicate these findings in our Irish Case-Control Study of Schizophrenia (ICCSS) sample (N=1021 cases, 626 controls). Following consultation with the original investigators, we genotyped three of the most promising single-nucleotide polymorphisms (SNPs) from the Cardiff study. We replicate association with rs1344706 (trend test one-tailed P=0.0113 with the previously associated A allele) in ZNF804A. We detect no evidence of association with rs6490121 in NOS1 (one-tailed P=0.21), and only a trend with rs9922369 in RGRIP1L (one-tailed P=0.0515). On the basis of these results, we completed genotyping of 11 additional linkage disequilibrium-tagging SNPs in ZNF804A. Of 12 SNPs genotyped, 11 pass quality control criteria and 4 are nominally associated, with our most significant evidence of association at rs7597593 (P=0.0013) followed by rs1344706. We observe no evidence of differential association in ZNF804A on the basis of family history or sex of case. The associated SNP rs1344706 lies in approximately 30 bp of conserved mammalian sequence, and the associated A allele is predicted to maintain binding sites for the brain-expressed transcription factors MYT1l and POU3F1/OCT-6. In controls, expression is significantly increased from the A allele of rs1344706 compared with the C allele. Expression is increased in schizophrenic cases compared with controls, but this difference does not achieve statistical significance. This study replicates the original reported association of ZNF804A with schizophrenia and suggests that there is a consistent link between the A allele of rs1344706, increased expression of ZNF804A and risk for schizophrenia.

The association of polygenic risk for schizophrenia, bipolar disorder, and depression with neural connectivity in adolescents and young adults: examining developmental and sex differences.

Neurodevelopmental abnormalities in neural connectivity have been long implicated in the etiology of schizophrenia (SCZ); however, it remains unclear whether these neural connectivity patterns are associated with genetic risk for SCZ in unaffected individuals (i.e., an absence of clinical features of SCZ or a family history of SCZ). We examine whether polygenic risk scores (PRS) for SCZ are associated with functional neural connectivity in adolescents and young adults without SCZ, whether this association is moderated by sex and age, and if similar associations are observed for genetically related neuropsychiatric PRS. One-thousand four-hundred twenty-six offspring from 913 families, unaffected with SCZ, were drawn from the Collaborative Study of the Genetics of Alcoholism (COGA) prospective cohort (median age at first interview = 15.6 (12-26), 51.6% female, 98.1% European American, 41% with a family history of alcohol dependence). Participants were followed longitudinally with resting-state EEG connectivity (i.e., coherence) assessed every two years. Higher SCZ PRS were associated with elevated theta (3-7 Hz) and alpha (7-12 Hz) EEG coherence. Associations differed by sex and age; the most robust associations were observed between PRS and parietal-occipital, central-parietal, and frontal-parietal alpha coherence among males between ages 15-19 (B: 0.15-0.21, p < 10-4). Significant associations among EEG coherence and Bipolar and Depression PRS were observed, but differed from SCZ PRS in terms of sex, age, and topography. Findings reveal that polygenic risk for SCZ is robustly associated with increased functional neural connectivity among young adults without a SCZ diagnosis. Striking differences were observed between men and women throughout development, mapping onto key periods of risk for the onset of psychotic illness and underlining the critical importance of examining sex differences in associations with neuropsychiatric PRS across development.

Association analysis of the PIP4K2A gene on chromosome 10p12 and schizophrenia in the Irish study of high density schizophrenia families (ISHDSF) and the Irish case-control study of schizophrenia (ICCSS).

Molecular studies support pharmacological evidence that phosphoinositide signaling is perturbed in schizophrenia and bipolar disorder. The phosphatidylinositol-4-phosphate-5-kinase type-II alpha (PIP4K2A) gene is located on chromosome 10p12. This region has been implicated in both diseases by linkage, and PIP4K2A directly by association. Given linkage evidence in the Irish Study of High Density Schizophrenia Families (ISHDSF) to a region including 10p12, we performed an association study between genetic variants at PIP4K2A and disease. No association was detected through single-marker or haplotype analysis of the whole sample. However, stratification into families positive and negative for the ISHDSF schizophrenia high-risk haplotype (HRH) in the DTNBP1 gene and re-analysis for linkage showed reduced amplitude of the 10p12 linkage peak in the DTNBP1 HRH positive families. Association analysis of the stratified sample showed a trend toward association of PIP4K2A SNPs rs1417374 and rs1409395 with schizophrenia in the DTNBP1 HRH positive families. Despite this apparent paradox, our data may therefore suggest involvement of PIP4K2A in schizophrenia in those families for whom genetic variation in DTNBP1 appears also to be a risk factor. This trend appears to arise from under-transmission of common alleles to female cases. Follow-up association analysis in a large Irish schizophrenia case-control sample (ICCSS) showed significant association with disease of a haplotype comprising these same SNPs rs1417374-rs1409395, again more so in affected females, and in cases with negative family history of the disease. This study supports a minor role for PIP4K2A in schizophrenia etiology in the Irish population.

Association study of SNAP25 and schizophrenia in Irish family and case-control samples.

SNAP25 occurs on chromosome 20p12.2, which has been linked to schizophrenia in some samples, and recently linked to latent classes of psychotic illness in our sample. SNAP25 is crucial to synaptic functioning, may be involved in axonal growth and dendritic sprouting, and its expression may be decreased in schizophrenia. We genotyped 18 haplotype-tagging SNPs in SNAP25 in a sample of 270 Irish high-density families. Single marker and haplotype analyses were performed in FBAT and PDT. We adjusted for multiple testing by computing q values. Association was followed up in an independent sample of 657 cases and 411 controls. We tested for allelic effects on the clinical phenotype by using the method of sequential addition and 5 factor-derived scores of the OPCRIT. Nine of 18 SNPs had P values <0.05 in either FBAT or PDT for one or more definitions of illness. Several two-marker haplotypes were also associated. Subjects inheriting the risk alleles of the most significantly associated two-marker haplotype were likely to have higher levels of hallucinations and delusions. The most significantly associated marker, rs6039820, was observed to perturb 12 transcription-factor binding sites in in silico analyses. An attempt to replicate association findings in the case-control sample resulted in no SNPs being significantly associated. We observed robust association in both single marker and haplotype-based analyses between SNAP25 and schizophrenia in an Irish family sample. Although we failed to replicate this in an independent sample, this gene should be further tested in other samples.

Genetic effects influencing risk for major depressive disorder in China and Europe.

Major depressive disorder (MDD) is a common, complex psychiatric disorder and a leading cause of disability worldwide. Despite twin studies indicating its modest heritability (~30-40%), extensive heterogeneity and a complex genetic architecture have complicated efforts to detect associated genetic risk variants. We combined single-nucleotide polymorphism (SNP) summary statistics from the CONVERGE and PGC studies of MDD, representing 10 502 Chinese (5282 cases and 5220 controls) and 18 663 European (9447 cases and 9215 controls) subjects. We determined the fraction of SNPs displaying consistent directions of effect, assessed the significance of polygenic risk scores and estimated the genetic correlation of MDD across ancestries. Subsequent trans-ancestry meta-analyses combined SNP-level evidence of association. Sign tests and polygenic score profiling weakly support an overlap of SNP effects between East Asian and European populations. We estimated the trans-ancestry genetic correlation of lifetime MDD as 0.33; female-only and recurrent MDD yielded estimates of 0.40 and 0.41, respectively. Common variants downstream of GPHN achieved genome-wide significance by Bayesian trans-ancestry meta-analysis (rs9323497; log10 Bayes Factor=8.08) but failed to replicate in an independent European sample (P=0.911). Gene-set enrichment analyses indicate enrichment of genes involved in neuronal development and axonal trafficking. We successfully demonstrate a partially shared polygenic basis of MDD in East Asian and European populations. Taken together, these findings support a complex etiology for MDD and possible population differences in predisposing genetic factors, with important implications for future genetic studies.

SNP-based heritability estimates of the personality dimensions and polygenic prediction of both neuroticism and major depression: findings from CONVERGE.

Biometrical genetic studies suggest that the personality dimensions, including neuroticism, are moderately heritable (~0.4 to 0.6). Quantitative analyses that aggregate the effects of many common variants have recently further informed genetic research on European samples. However, there has been limited research to date on non-European populations. This study examined the personality dimensions in a large sample of Han Chinese descent (N=10 064) from the China, Oxford, and VCU Experimental Research on Genetic Epidemiology study, aimed at identifying genetic risk factors for recurrent major depression among a rigorously ascertained cohort. Heritability of neuroticism as measured by the Eysenck Personality Questionnaire (EPQ) was estimated to be low but statistically significant at 10% (s.e.=0.03, P=0.0001). In addition to EPQ, neuroticism based on a three-factor model, data for the Big Five (BF) personality dimensions (neuroticism, openness, conscientiousness, extraversion and agreeableness) measured by the Big Five Inventory were available for controls (n=5596). Heritability estimates of the BF were not statistically significant despite high power (>0.85) to detect heritabilities of 0.10. Polygenic risk scores constructed by best linear unbiased prediction weights applied to split-half samples failed to significantly predict any of the personality traits, but polygenic risk for neuroticism, calculated with LDpred and based on predictive variants previously identified from European populations (N=171 911), significantly predicted major depressive disorder case-control status (P=0.0004) after false discovery rate correction. The scores also significantly predicted EPQ neuroticism (P=6.3 × 10-6). Factor analytic results of the measures indicated that any differences in heritabilities across samples may be due to genetic variation or variation in haplotype structure between samples, rather than measurement non-invariance. Findings demonstrate that neuroticism can be significantly predicted across ancestry, and highlight the importance of studying polygenic contributions to personality in non-European populations.

Genome-wide association study of lifetime cannabis use based on a large meta-analytic sample of 32 330 subjects from the International Cannabis Consortium.

Cannabis is the most widely produced and consumed illicit psychoactive substance worldwide. Occasional cannabis use can progress to frequent use, abuse and dependence with all known adverse physical, psychological and social consequences. Individual differences in cannabis initiation are heritable (40-48%). The International Cannabis Consortium was established with the aim to identify genetic risk variants of cannabis use. We conducted a meta-analysis of genome-wide association data of 13 cohorts (N=32 330) and four replication samples (N=5627). In addition, we performed a gene-based test of association, estimated single-nucleotide polymorphism (SNP)-based heritability and explored the genetic correlation between lifetime cannabis use and cigarette use using LD score regression. No individual SNPs reached genome-wide significance. Nonetheless, gene-based tests identified four genes significantly associated with lifetime cannabis use: NCAM1, CADM2, SCOC and KCNT2. Previous studies reported associations of NCAM1 with cigarette smoking and other substance use, and those of CADM2 with body mass index, processing speed and autism disorders, which are phenotypes previously reported to be associated with cannabis use. Furthermore, we showed that, combined across the genome, all common SNPs explained 13-20% (P<0.001) of the liability of lifetime cannabis use. Finally, there was a strong genetic correlation (rg=0.83; P=1.85 × 10(-8)) between lifetime cannabis use and lifetime cigarette smoking implying that the SNP effect sizes of the two traits are highly correlated. This is the largest meta-analysis of cannabis GWA studies to date, revealing important new insights into the genetic pathways of lifetime cannabis use. Future functional studies should explore the impact of the identified genes on the biological mechanisms of cannabis use.

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.

Genome-wide association data suggest ABCB1 and immune-related gene sets may be involved in adult antisocial behavior.

Adult antisocial behavior (AAB) is moderately heritable, relatively common and has adverse consequences for individuals and society. We examined the molecular genetic basis of AAB in 1379 participants from a case-control study in which the cases met criteria for alcohol dependence. We also examined whether genes of interest were expressed in human brain. AAB was measured using a count of the number of Antisocial Personality Disorder criteria endorsed under criterion A from the Diagnostic and Statistical Manual of Mental Disorders, 4th Edition (DSM-IV). Participants were genotyped on the Illumina Human 1M BeadChip. In total, all single-nucleotide polymorphisms (SNPs) accounted for 25% of the variance in AAB, although this estimate was not significant (P=0.09). Enrichment tests indicated that more significantly associated genes were over-represented in seven gene sets, and most were immune related. Our most highly associated SNP (rs4728702, P=5.77 × 10(-7)) was located in the protein-coding adenosine triphosphate-binding cassette, sub-family B (MDR/TAP), member 1 (ABCB1). In a gene-based test, ABCB1 was genome-wide significant (q=0.03). Expression analyses indicated that ABCB1 was robustly expressed in the brain. ABCB1 has been implicated in substance use, and in post hoc tests we found that variation in ABCB1 was associated with DSM-IV alcohol and cocaine dependence criterion counts. These results suggest that ABCB1 may confer risk across externalizing behaviors, and are consistent with previous suggestions that immune pathways are associated with externalizing behaviors. The results should be tempered by the fact that we did not replicate the associations for ABCB1 or the gene sets in a less-affected independent sample.

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.