Functional effects of schizophrenia-linked genetic variants on intrinsic single-neuron excitability: A modeling study

Background: Recent genome-wide association studies (GWAS) have identified a large number of genetic risk factors for schizophrenia (SCZ) featuring ion channels and calcium transporters. For some of these risk factors, independent prior investigations have examined the effects of genetic alterations on the cellular electrical excitability and calcium homeostasis. In the present proof-of-concept study, we harnessed these experimental results for modeling of computational properties on layer V cortical pyramidal cell and identify possible common alterations in behavior across SCZ-related genes. Methods: We applied a biophysically detailed multi-compartmental model to study the excitability of a layer V pyramidal cell. We reviewed the literature on functional genomics for variants of genes associated with SCZ, and used changes in neuron model parameters to represent the effects of these variants. Results: We present and apply a framework for examining the effects of subtle single nucleotide polymorphisms in ion channel and Ca2+ transporter-encoding genes on neuron excitability. Our analysis indicates that most of the considered SCZ- related genetic variants affect the spiking behavior and intracellular calcium dynamics resulting from summation of inputs across the dendritic tree. Conclusions: Our results suggest that alteration in the ability of a single neuron to integrate the inputs and scale its excitability may constitute a fundamental mechanistic contributor to mental disease, alongside with the previously proposed deficits in synaptic communication and network behavior

Cross-tissue eQTL enrichment of associations in schizophrenia.

The genome-wide association study of the Psychiatric Genomics Consortium identified over one hundred schizophrenia susceptibility loci. The number of non-coding variants discovered suggests that gene regulation could mediate the effect of these variants on disease. Expression quantitative trait loci (eQTLs) contribute to variation in levels of mRNA. Given the co-occurrence of schizophrenia and several traits not involving the central nervous system (CNS), we investigated the enrichment of schizophrenia associations among eQTLs for four non-CNS tissues: adipose tissue, epidermal tissue, lymphoblastoid cells and blood. Significant enrichment was seen in eQTLs of all tissues: adipose (β = 0.18, p = 8.8 × 10−06), epidermal (β = 0.12, p = 3.1 × 10−04), lymphoblastoid (β = 0.19, p = 6.2 × 10−08) and blood (β = 0.19, p = 6.4 × 10−06). For comparison, we looked for enrichment of association with traits of known relevance to one or more of these tissues (body mass index, height, rheumatoid arthritis, systolic blood pressure and type-II diabetes) and found that schizophrenia enrichment was of similar scale to that observed when studying diseases in the context of a more likely causal tissue. To further investigate tissue specificity, we looked for differential enrichment of eQTLs with relevant Roadmap affiliation (enhancers and promoters) and varying distance from the transcription start site. Neither factor significantly contributed to the enrichment, suggesting that this is equally distributed in tissue-specific and cross-tissue regulatory elements. Our analyses suggest that functional correlates of schizophrenia risk are prevalent in non-CNS tissues. This could be because of pleiotropy or the effectiveness of variants affecting expression in different contexts. This suggests the utility of large, single-tissue eQTL experiments to increase eQTL discovery power in the study of schizophrenia, in addition to smaller, multiple-tissue approaches. Our results conform to the notion that schizophrenia is a systemic disorder involving many tissues

Bromodomain protein 4 discriminates tissue-specific super-enhancers containing disease-specific susceptibility loci in prostate and breast cancer.

Background Epigenetic information can be used to identify clinically relevant genomic variants single nucleotide polymorphisms (SNPs) of functional importance in cancer development. Super-enhancers are cell-specific DNA elements, acting to determine tissue or cell identity and driving tumor progression. Although previous approaches have been tried to explain risk associated with SNPs in regulatory DNA elements, so far epigenetic readers such as bromodomain containing protein 4 (BRD4) and super-enhancers have not been used to annotate SNPs. In prostate cancer (PC), androgen receptor (AR) binding sites to chromatin have been used to inform functional annotations of SNPs.Results Here we establish criteria for enhancer mapping which are applicable to other diseases and traits to achieve the optimal tissue-specific enrichment of PC risk SNPs. We used stratified Q-Q plots and Fisher test to assess the differential enrichment of SNPs mapping to specific categories of enhancers. We find that BRD4 is the key discriminant of tissue-specific enhancers, showing that it is more powerful than AR binding information to capture PC specific risk loci, and can be used with similar effect in breast cancer (BC) and applied to other diseases such as schizophrenia.Conclusions This is the first study to evaluate the enrichment of epigenetic readers in genome-wide associations studies for SNPs within enhancers, and provides a powerful tool for enriching and prioritizing PC and BC genetic risk loci. Our study represents a proof of principle applicable to other diseases and traits that can be used to redefine molecular mechanisms of human phenotypic variation

Estimating effect sizes and expected replication probabilities from GWAS summary statistics

Genome-wide Association Studies (GWAS) result in millions of summary statistics (“z-scores”) for single nucleotide polymorphism (SNP) associations with phenotypes. These rich datasets afford deep insights into the nature and extent of genetic contributions to complex phenotypes such as psychiatric disorders, which are understood to have substantial genetic components that arise from very large numbers of SNPs. The complexity of the datasets, however, poses a significant challenge to maximizing their utility. This is reflected in a need for better understanding the landscape of z-scores, as such knowledge would enhance causal SNP and gene discovery, help elucidate mechanistic pathways, and inform future study design. Here we present a parsimonious methodology for modeling effect sizes and replication probabilities, relying only on summary statistics from GWAS substudies, and a scheme allowing for direct empirical validation. We show that modeling z-scores as a mixture of Gaussians is conceptually appropriate, in particular taking into account ubiquitous non-null effects that are likely in the datasets due to weak linkage disequilibrium with causal SNPs. The four-parameter model allows for estimating the degree of polygenicity of the phenotype and predicting the proportion of chip heritability explainable by genome-wide significant SNPs in future studies with larger sample sizes. We apply the model to recent GWAS of schizophrenia (N = 82,315) and putamen volume (N = 12,596), with approximately 9.3 million SNP z-scores in both cases. We show that, over a broad range of z-scores and sample sizes, the model accurately predicts expectation estimates of true effect sizes and replication probabilities in multistage GWAS designs. We assess the degree to which effect sizes are over-estimated when based on linear-regression association coefficients. We estimate the polygenicity of schizophrenia to be 0.037 and the putamen to be 0.001, while the respective sample sizes required to approach fully explaining the chip heritability are 106 and 105. The model can be extended to incorporate prior knowledge such as pleiotropy and SNP annotation. The current findings suggest that the model is applicable to a broad array of complex phenotypes and will enhance understanding of their genetic architectures

Identification of shared genetic variants between schizophrenia and lung cancer.

Epidemiology studies suggest associations between schizophrenia and cancer. However, the underlying genetic mechanisms are not well understood, and difficult to identify from epidemiological data. We investigated if there is a shared genetic architecture between schizophrenia and cancer, with the aim to identify specific overlapping genetic loci. First, we performed genome-wide enrichment analysis and second, we analyzed specific loci jointly associated with schizophrenia and cancer by the conjunction false discovery rate. We analyzed the largest genome-wide association studies of schizophrenia and lung, breast, prostate, ovary, and colon-rectum cancer including more than 220,000 subjects, and included genetic association with smoking behavior. Polygenic enrichment of associations with lung cancer was observed in schizophrenia, and weak enrichment for the remaining cancer sites. After excluding the major histocompatibility complex region, we identified three independent loci jointly associated with schizophrenia and lung cancer. The strongest association included nicotinic acetylcholine receptors and is an established pleiotropic locus shared between lung cancer and smoking. The two other loci were independent of genetic association with smoking. Functional analysis identified downstream pleiotropic effects on epigenetics and gene-expression in lung and brain tissue. These findings suggest that genetic factors may explain partly the observed epidemiological association of lung cancer and schizophrenia

Leveraging genomic annotations and pleiotropic enrichment for improved replication rates in schizophrenia GWAS

Most of the genetic architecture of schizophrenia (SCZ) has not yet been identified. Here, we apply a novel statistical algorithm called Covariate-Modulated Mixture Modeling (CM3), which incorporates auxiliary information (heterozygosity, total linkage disequilibrium, genomic annotations, pleiotropy) for each single nucleotide polymorphism (SNP) to enable more accurate estimation of replication probabilities, conditional on the observed test statistic (“z-score”) of the SNP. We use a multiple logistic regression on z-scores to combine information from auxiliary information to derive a “relative enrichment score” for each SNP. For each stratum of these relative enrichment scores, we obtain nonparametric estimates of posterior expected test statistics and replication probabilities as a function of discovery z-scores, using a resampling-based approach that repeatedly and randomly partitions meta-analysis sub-studies into training and replication samples. We fit a scale mixture of two Gaussians model to each stratum, obtaining parameter estimates that minimize the sum of squared differences of the scale-mixture model with the stratified nonparametric estimates. We apply this approach to the recent genome-wide association study (GWAS) of SCZ (n = 82,315), obtaining a good fit between the model-based and observed effect sizes and replication probabilities. We observed that SNPs with low enrichment scores replicate with a lower probability than SNPs with high enrichment scores even when both they are genome-wide significant (p < 5x10-8). There were 693 and 219 independent loci with model-based replication rates ≥80% and ≥90%, respectively. Compared to analyses not incorporating relative enrichment scores, CM3 increased out-of-sample yield for SNPs that replicate at a given rate. This demonstrates that replication probabilities can be more accurately estimated using prior enrichment information with CM3

GBA and APOE ε4 associate with sporadic dementia with Lewy bodies in European genome wide association study

Dementia with Lewy Bodies (DLB) is a common neurodegenerative disorder with poor prognosis and mainly unknown pathophysiology. Heritability estimates exceed 30% but few genetic risk variants have been identified. Here we investigated common genetic variants associated with DLB in a large European multisite sample. We performed a genome wide association study in Norwegian and European cohorts of 720 DLB cases and 6490 controls and included 19 top-associated single-nucleotide polymorphisms in an additional cohort of 108 DLB cases and 75545 controls from Iceland. Overall the study included 828 DLB cases and 82035 controls. Variants in the ASH1L/GBA (Chr1q22) and APOE ε4 (Chr19) loci were associated with DLB surpassing the genome-wide significance threshold (p < 5 × 10). One additional genetic locus previously linked to psychosis in Alzheimer's disease, ZFPM1 (Chr16q24.2), showed suggestive association with DLB at p-value < 1 × 10. We report two susceptibility loci for DLB at genome-wide significance, providing insight into etiological factors. These findings highlight the complex relationship between the genetic architecture of DLB and other neurodegenerative disorders