Oncogenomic studies indicate that copy number variation (CNV) alters genes involved in tumor progression; however, identification of specific driver genes affected by CNV has been difficult, as these rearrangements are often contained in large chromosomal intervals among several bystander genes. Here, we addressed this problem and identified a CNV-targeted oncogene by performing comparative oncogenomics of human and zebrafish melanomas. We determined that the gene encoding growth differentiation factor 6 (GDF6), which is the ligand for the BMP family, is recurrently amplified and transcriptionally upregulated in melanoma. GDF6-induced BMP signaling maintained a trunk neural crest gene signature in melanomas. Additionally, GDF6 repressed the melanocyte differentiation gene MITF and the proapoptotic factor SOX9, thereby preventing differentiation, inhibiting cell death, and promoting tumor growth. GDF6 was specifically expressed in melanomas but not melanocytes. Moreover, GDF6 expression levels in melanomas were inversely correlated with patient survival. Our study has identified a fundamental role for GDF6 and BMP signaling in governing an embryonic cell gene signature to promote melanoma progression, thus providing potential opportunities for targeted therapy to treat GDF6-positive cancers.

We generated cortical interneurons (cINs) from induced pluripotent stem cells derived from 14 healthy controls and 14 subjects with schizophrenia. Both healthy control cINs and schizophrenia cINs were authentic, fired spontaneously, received functional excitatory inputs from host neurons, and induced GABA-mediated inhibition in host neurons in vivo. However, schizophrenia cINs had dysregulated expression of protocadherin genes, which lie within documented schizophrenia loci. Mice lacking protocadherin-α showed defective arborization and synaptic density of prefrontal cortex cINs and behavioral abnormalities. Schizophrenia cINs similarly showed defects in synaptic density and arborization that were reversed by inhibitors of protein kinase C, a downstream kinase in the protocadherin pathway. These findings reveal an intrinsic abnormality in schizophrenia cINs in the absence of any circuit-driven pathology. They also demonstrate the utility of homogenous and functional populations of a relevant neuronal subtype for probing pathogenesis mechanisms during development.

We generated cortical interneurons (cINs) from induced pluripotent stem cells derived from 14 healthy controls and 14 subjects with schizophrenia. Both healthy control cINs and schizophrenia cINs were authentic, fired spontaneously, received functional excitatory inputs from host neurons, and induced GABA-mediated inhibition in host neurons in vivo. However, schizophrenia cINs had dysregulated expression of protocadherin genes, which lie within documented schizophrenia loci. Mice lacking protocadherin-α showed defective arborization and synaptic density of prefrontal cortex cINs and behavioral abnormalities. Schizophrenia cINs similarly showed defects in synaptic density and arborization that were reversed by inhibitors of protein kinase C, a downstream kinase in the protocadherin pathway. These findings reveal an intrinsic abnormality in schizophrenia cINs in the absence of any circuit-driven pathology. They also demonstrate the utility of homogenous and functional populations of a relevant neuronal subtype for probing pathogenesis mechanisms during development.

Genome-wide association studies (GWASs) are highly effective at identifying common risk variants for schizophrenia. Rare risk variants are also important contributors to schizophrenia etiology but, with the exception of large copy number variants, are difficult to detect with GWAS. Exome and genome sequencing, which have accelerated the study of rare variants, are expensive so alternative methods are needed to aid detection of rare variants. Here we re-analyze an Irish schizophrenia GWAS dataset (n = 3,473) by performing identity-by-descent (IBD) mapping followed by exome sequencing of individuals identified as sharing risk haplotypes to search for rare risk variants in coding regions. We identified 45 rare haplotypes (>1 cM) that were significantly more common in cases than controls. By exome sequencing 105 haplotype carriers, we investigated these haplotypes for functional coding variants that could be tested for association in independent GWAS samples. We identified one rare missense variant in PCNT but did not find statistical support for an association with schizophrenia in a replication analysis. However, IBD mapping can prioritize both individual samples and genomic regions for follow-up analysis but genome rather than exome sequencing may be more effective at detecting risk variants on rare haplotypes.

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.

Schizophrenia is a highly heritable disorder. Genetic risk is conferred by a large number of alleles, including common alleles of small effect that might be detected by genome-wide association studies. Here we report a multi-stage schizophrenia genome-wide association study of up to 36,989 cases and 113,075 controls. We identify 128 independent associations spanning 108 conservatively defined loci that meet genome-wide significance, 83 of which have not been previously reported. Associations were enriched among genes expressed in brain, providing biological plausibility for the findings. Many findings have the potential to provide entirely new insights into aetiology, but associations at DRD2 and several genes involved in glutamatergic neurotransmission highlight molecules of known and potential therapeutic relevance to schizophrenia, and are consistent with leading pathophysiological hypotheses. Independent of genes expressed in brain, associations were enriched among genes expressed in tissues that have important roles in immunity, providing support for the speculated link between the immune system and schizophrenia.

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.