Modeling linkage disequilibrium increases accuracy of polygenic risk scores

Polygenic risk scores have shown great promise in predicting complex disease risk and will become more accurate as training sample sizes increase. The standard approach for calculating risk scores involves linkage disequilibrium (LD)-based marker pruning and applying a p value threshold to association statistics, but this discards information and can reduce predictive accuracy. We introduce LDpred, a method that infers the posterior mean effect size of each marker by using a prior on effect sizes and LD information from an external reference panel. Theory and simulations show that LDpred outperforms the approach of pruning followed by thresholding, particularly at large sample sizes. Accordingly, predicted R(2) increased from 20.1% to 25.3% in a large schizophrenia dataset and from 9.8% to 12.0% in a large multiple sclerosis dataset. A similar relative improvement in accuracy was observed for three additional large disease datasets and for non-European schizophrenia samples. The advantage of LDpred over existing methods will grow as sample sizes increase

Comparative genetic architectures of schizophrenia in East Asian and European populations

Schizophrenia is a debilitating psychiatric disorder with approximately 1% lifetime risk globally. Large-scale schizophrenia genetic studies have reported primarily on European ancestry samples, potentially missing important biological insights. Here, we report the largest study to date of East Asian participants (22,778 schizophrenia cases and 35,362 controls), identifying 21 genome-wide-significant associations in 19 genetic loci. Common genetic variants that confer risk for schizophrenia have highly similar effects between East Asian and European ancestries (genetic correlation = 0.98 ± 0.03), indicating that the genetic basis of schizophrenia and its biology are broadly shared across populations. A fixed-effect meta-analysis including individuals from East Asian and European ancestries identified 208 significant associations in 176 genetic loci (53 novel). Trans-ancestry fine-mapping reduced the sets of candidate causal variants in 44 loci. Polygenic risk scores had reduced performance when transferred across ancestries, highlighting the importance of including sufficient samples of major ancestral groups to ensure their generalizability across populations

Identification of Gene Loci That Overlap Between Schizophrenia and Educational Attainment

There is evidence for genetic overlap between cognitive abilities and schizophrenia (SCZ), and genome-wide association studies (GWAS) demonstrate that both SCZ and general cognitive abilities have a strong polygenic component with many single-nucleotide polymorphisms (SNPs) each with a small effect. Here we investigated the shared genetic architecture between SCZ and educational attainment, which is regarded as a "proxy phenotype" for cognitive abilities, but may also reflect other traits. We applied a conditional false discovery rate (condFDR) method to GWAS of SCZ (n = 82 315), college completion ("College," n = 95 427), and years of education ("EduYears," n = 101 069). Variants associated with College or EduYears showed enrichment of association with SCZ, demonstrating polygenic overlap. This was confirmed by an increased replication rate in SCZ. By applying a condFDR threshold <0.01, we identified 18 genomic loci associated with SCZ after conditioning on College and 15 loci associated with SCZ after conditioning on EduYears. Ten of these loci overlapped. Using conjunctional FDR, we identified 10 loci shared between SCZ and College, and 29 loci shared between SCZ and EduYears. The majority of these loci had effects in opposite directions. Our results provide evidence for polygenic overlap between SCZ and educational attainment, and identify novel pleiotropic loci. Other studies have reported genetic overlap between SCZ and cognition, or SCZ and educational attainment, with negative correlation. Importantly, our methods enable identification of bi-directional effects, which highlight the complex relationship between SCZ and educational attainment, and support polygenic mechanisms underlying both cognitive dysfunction and creativity in SCZ

Transcriptome-wide association study of schizophrenia and chromatin activity yields mechanistic disease insights

Genome-wide association studies (GWAS) have identified over 100 risk loci for schizophrenia, but the causal mechanisms remain largely unknown. We performed a transcriptome-wide association study (TWAS) integrating a schizophrenia GWAS of 79,845 individuals from the Psychiatric Genomics Consortium with expression data from brain, blood, and adipose tissues across 3,693 primarily control individuals. We identified 157 TWAS-significant genes, of which 35 did not overlap a known GWAS locus. Of these 157 genes, 42 were associated with specific chromatin features measured in independent samples, thus highlighting potential regulatory targets for follow-up. Suppression of one identified susceptibility gene, mapk3, in zebrafish showed a significant effect on neurodevelopmental phenotypes. Expression and splicing from the brain captured most of the TWAS effect across all genes. This large-scale connection of associations to target genes, tissues, and regulatory features is an essential step in moving toward a mechanistic understanding of GWAS

A molecule-based genetic association approach implicates a range of voltage-gated calcium channels associated with schizophrenia

Traditional genome-wide association studies (GWAS) have successfully detected genetic variants associated with schizophrenia. However, only a small fraction of heritability can be explained. Gene-set/pathway-based methods can overcome limitations arising from single nucleotide polymorphism (SNP)-based analysis, but most of them place constraints on size which may exclude highly specific and functional sets, like macromolecules. Voltage-gated calcium (Cav ) channels, belonging to macromolecules, are composed of several subunits whose encoding genes are located far away or even on different chromosomes. We combined information about such molecules with GWAS data to investigate how functional channels associated with schizophrenia. We defined a biologically meaningful SNP-set based on channel structure and performed an association study by using a validated method: SNP-set (sequence) kernel association test. We identified eight subtypes of Cav channels significantly associated with schizophrenia from a subsample of published data (N = 56,605), including the L-type channels (Cav 1.1, Cav 1.2, Cav 1.3), P-/Q-type Cav 2.1, N-type Cav 2.2, R-type Cav 2.3, T-type Cav 3.1, and Cav 3.3. Only genes from Cav 1.2 and Cav 3.3 have been implicated by the largest GWAS (N = 82,315). Each subtype of Cav channels showed relatively high chip heritability, proportional to the size of its constituent gene regions. The results suggest that abnormalities of Cav channels may play an important role in the pathophysiology of schizophrenia and these channels may represent appropriate drug targets for therapeutics. Analyzing subunit-encoding genes of a macromolecule in aggregate is a complementary way to identify more genetic variants of polygenic diseases. This study offers the potential of power for discovery the biological mechanisms of schizophrenia

A polygenic resilience score moderates the genetic risk for schizophrenia.

Based on the discovery by the Resilience Project (Chen R. et al. Nat Biotechnol 34:531–538, 2016) of rare variants that confer resistance to Mendelian disease, and protective alleles for some complex diseases, we posited the existence of genetic variants that promote resilience to highly heritable polygenic disorders1,0 such as schizophrenia. Resilience has been traditionally viewed as a psychological construct, although our use of the term resilience refers to a different construct that directly relates to the Resilience Project, namely: heritable variation that promotes resistance to disease by reducing the penetrance of risk loci, wherein resilience and risk loci operate orthogonal to one another. In this study, we established a procedure to identify unaffected individuals with relatively high polygenic risk for schizophrenia, and contrasted them with risk-matched schizophrenia cases to generate the first known “polygenic resilience score” that represents the additive contributions to SZ resistance by variants that are distinct from risk loci. The resilience score was derived from data compiled by the Psychiatric Genomics Consortium, and replicated in three independent samples. This work establishes a generalizable framework for finding resilience variants for any complex, heritable disorder

The role of polygenic risk score gene-set analysis in the context of the omnigenic model of schizophrenia

A recent development in the genetic architecture of schizophrenia suggested that an omnigenic model may underlie the risk for this disorder. The aim of our study was to use polygenic profile scoring to quantitatively assess whether a number of experimentally derived sets would contribute to the disorder above and beyond the omnigenic effect. Using the PGC2 secondary analysis schizophrenia case-control cohort (N = 29,125 cases and 34,836 controls), a robust polygenic signal was observed from gene sets based on TCF4, FMR1, upregulation from MIR137 and downregulation from CHD8. Additional analyses revealed a constant floor effect in the amount of variance explained, consistent with the omnigenic model. Thus, we report that putative core gene sets showed a significant effect above and beyond the floor effect that might be linked with the underlying omnigenic background. In addition, we demonstrate a method to quantify the contribution of specific gene sets within the omnigenic context

Complement genes contribute sex-biased vulnerability in diverse disorders

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren’s syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren’s syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3,4,5,6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren’s syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren’s syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren’s syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women’s greater risk of SLE and Sjögren’s syndrome and men’s greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses

Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs (hg(2)) across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of hg(2) from imputed SNPs (5.1× enrichment; p = 3.7 × 10(-17)) and 38% (SE = 4%) of hg(2) from genotyped SNPs (1.6× enrichment, p = 1.0 × 10(-4)). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of hg(2) despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease

Genetic Markers of Human Evolution Are Enriched in Schizophrenia

BACKGROUND: Why schizophrenia has accompanied humans throughout our history despite its negative effect on fitness remains an evolutionary enigma. It is proposed that schizophrenia is a by-product of the complex evolution of the human brain and a compromise for humans' language, creative thinking, and cognitive abilities. METHODS: We analyzed recent large genome-wide association studies of schizophrenia and a range of other human phenotypes (anthropometric measures, cardiovascular disease risk factors, immune-mediated diseases) using a statistical framework that draws on polygenic architecture and ancillary information on genetic variants. We used information from the evolutionary proxy measure called the Neanderthal selective sweep (NSS) score. RESULTS: Gene loci associated with schizophrenia are significantly (p = 7.30 × 10(-9)) more prevalent in genomic regions that are likely to have undergone recent positive selection in humans (i.e., with a low NSS score). Variants in brain-related genes with a low NSS score confer significantly higher susceptibility than variants in other brain-related genes. The enrichment is strongest for schizophrenia, but we cannot rule out enrichment for other phenotypes. The false discovery rate conditional on the evolutionary proxy points to 27 candidate schizophrenia susceptibility loci, 12 of which are associated with schizophrenia and other psychiatric disorders or linked to brain development. CONCLUSIONS: Our results suggest that there is a polygenic overlap between schizophrenia and NSS score, a marker of human evolution, which is in line with the hypothesis that the persistence of schizophrenia is related to the evolutionary process of becoming human