Genetic influences on gene expression in the human fetal brain plausibly impact upon a variety of postnatal brain-related traits, including susceptibility to neuropsychiatric disorders. However, to ...date, there have been no studies that have mapped genome-wide expression quantitative trait loci (eQTL) specifically in the human prenatal brain.
We performed deep RNA sequencing and genome-wide genotyping on a unique collection of 120 human brains from the second trimester of gestation to provide the first eQTL dataset derived exclusively from the human fetal brain. We identify high confidence cis-acting eQTL at the individual transcript as well as whole gene level, including many mapping to a common inversion polymorphism on chromosome 17q21. Fetal brain eQTL are enriched among risk variants for postnatal conditions including attention deficit hyperactivity disorder, schizophrenia, and bipolar disorder. We further identify changes in gene expression within the prenatal brain that potentially mediate risk for neuropsychiatric traits, including increased expression of C4A in association with genetic risk for schizophrenia, increased expression of LRRC57 in association with genetic risk for bipolar disorder, and altered expression of multiple genes within the chromosome 17q21 inversion in association with variants influencing the personality trait of neuroticism.
We have mapped eQTL operating in the human fetal brain, providing evidence that these confer risk to certain neuropsychiatric disorders, and identifying gene expression changes that potentially mediate susceptibility to these conditions.
Schizophrenia is a highly polygenic disorder with important contributions from both common and rare risk alleles. We analyzed exome sequencing data for de novo variants (DNVs) in a new sample of 613 ...schizophrenia trios and combined this with published data to give a total of 3,444 trios. In this new data, loss-of-function (LoF) DNVs were significantly enriched among 3,471 LoF-intolerant genes, which supports previous findings. In the full dataset, genes associated with neurodevelopmental disorders (n = 159) were significantly enriched for LoF DNVs. Within these neurodevelopmental disorder genes, SLC6A1, which encodes a γ-aminobutyric acid transporter, was associated with missense-damaging DNVs. In 1,122 trios for which genome-wide common variant data were available, schizophrenia and bipolar disorder polygenic risk were significantly overtransmitted to probands. Probands carrying LoF or deletion DNVs in LoF-intolerant or neurodevelopmental disorder genes had significantly less overtransmission of schizophrenia polygenic risk than did non-carriers, which provides a second robust line of evidence that these DNVs increase liability to schizophrenia.
A recent genome-wide association study (GWAS) of autism spectrum disorder (ASD) (ncases = 18,381, ncontrols = 27,969) has provided novel opportunities for investigating the etiology of ASD. Here, we ...integrate the ASD GWAS summary statistics with summary-level gene expression data to infer differential gene expression in ASD, an approach called transcriptome-wide association study (TWAS).
Using FUSION software, ASD GWAS summary statistics were integrated with predictors of gene expression from 16 human datasets, including adult and fetal brains. A novel adaptation of established statistical methods was then used to test for enrichment within candidate pathways and specific tissues and at different stages of brain development. The proportion of ASD heritability explained by predicted expression of genes in the TWAS was estimated using stratified linkage disequilibrium score regression.
This study identified 14 genes as significantly differentially expressed in ASD, 13 of which were outside of known genome-wide significant loci (±500 kb). XRN2, a gene proximal to an ASD GWAS locus, was inferred to be significantly upregulated in ASD, providing insight into the functional consequence of this associated locus. One novel transcriptome-wide significant association from this study is the downregulation of PDIA6, which showed minimal evidence of association in the GWAS, and in gene-based analysis using MAGMA. Predicted gene expression in this study accounted for 13.0% of the total ASD single nucleotide polymorphism heritability.
This study has implicated several genes as significantly up/downregulated in ASD, providing novel and useful information for subsequent functional studies. This study also explores the utility of TWAS-based enrichment analysis and compares TWAS results with a functionally agnostic approach.
Abstract
Schizophrenia is a complex highly heritable disorder. Genome-wide association studies (GWAS) have identified multiple loci that influence the risk of developing schizophrenia, although the ...causal variants driving these associations and their impacts on specific genes are largely unknown. We identify a significant correlation between schizophrenia risk and expression at 89 genes in the dorsolateral prefrontal cortex (P ≤ 9.43 × 10−6), including 20 novel genes. Genes whose expression correlate with schizophrenia were enriched for those involved in abnormal CNS synaptic transmission (PFDR = 0.02) and antigen processing and presentation of peptide antigen via MHC class I (PFDR = 0.02). Within the CNS synaptic transmission set, we identify individual significant candidate genes to which we assign direction of expression changes in schizophrenia. The findings provide strong candidates for experimentally probing the molecular basis of synaptic pathology in schizophrenia.
The majority of common risk alleles identified for neuropsychiatric disorders reside in noncoding regions of the genome and are therefore likely to impact gene regulation. However, the genes that are ...primarily affected and the nature and developmental timing of these effects remain unclear. Given the hypothesized role for early neurodevelopmental processes in these conditions, we here define genetic predictors of gene expression in the human fetal brain with which we perform transcriptome-wide association studies (TWASs) of attention deficit hyperactivity disorder (ADHD), autism spectrum disorder, bipolar disorder, major depressive disorder, and schizophrenia. We identify prenatal cis-regulatory effects on 63 genes and 166 individual transcripts associated with genetic risk for these conditions. We observe pleiotropic effects of expression predictors for a number of genes and transcripts, including those of decreased DDHD2 expression in association with risk for schizophrenia and bipolar disorder, increased expression of a ST3GAL3 transcript with risk for schizophrenia and ADHD, and increased expression of an XPNPEP3 transcript with risk for schizophrenia, bipolar disorder, and major depression. For the protocadherin alpha cluster genes PCDHA7 and PCDHA8, we find that predictors of low expression are associated with risk for major depressive disorder while those of higher expression are associated with risk for schizophrenia. Our findings support a role for altered gene regulation in the prenatal brain in susceptibility to various neuropsychiatric disorders and prioritize potential risk genes for further neurobiological investigation.
Chronic pain is highly prevalent and a significant source of disability, yet its genetic and environmental risk factors are poorly understood. Its relationship with major depressive disorder (MDD) is ...of particular importance. We sought to test the contribution of genetic factors and shared and unique environment to risk of chronic pain and its correlation with MDD in Generation Scotland: Scottish Family Health Study (GS:SFHS). We then sought to replicate any significant findings in the United Kingdom Biobank study.
Using family-based mixed-model analyses, we examined the contribution of genetics and shared family environment to chronic pain by spouse, sibling, and household relationships. These analyses were conducted in GS:SFHS (n = 23,960), a family- and population-based study of individuals recruited from the Scottish population through their general practitioners. We then examined and partitioned the correlation between chronic pain and MDD and estimated the contribution of genetic factors and shared environment in GS:SFHS. Finally, we used data from two independent genome-wide association studies to test whether chronic pain has a polygenic architecture and examine whether genomic risk of psychiatric disorder predicted chronic pain and whether genomic risk of chronic pain predicted MDD. These analyses were conducted in GS:SFHS and repeated in UK Biobank, a study of 500,000 from the UK population, of whom 112,151 had genotyping and phenotypic data. Chronic pain is a moderately heritable trait (heritability = 38.4%, 95% CI 33.6% to 43.9%) that is significantly concordant in spouses (variance explained 18.7%, 95% CI 9.5% to 25.1%). Chronic pain is positively correlated with depression (ρ = 0.13, 95% CI 0.11 to 0.15, p = 2.72x10-68) and shows a tendency to cluster within families for genetic reasons (genetic correlation = 0.51, 95%CI 0.40 to 0.62, p = 8.24x10-19). Polygenic risk profiles for pain, generated using independent GWAS data, were associated with chronic pain in both GS:SFHS (maximum β = 6.18x10-2, 95% CI 2.84 x10-2 to 9.35 x10-2, p = 4.3x10-4) and UK Biobank (maximum β = 5.68 x 10-2, 95% CI 4.70x10-2 to 6.65x10-2, p < 3x10-4). Genomic risk of MDD is also significantly associated with chronic pain in both GS:SFHS (maximum β = 6.62x10-2, 95% CI 2.82 x10-2 to 9.76 x10-2, p = 4.3x10-4) and UK Biobank (maximum β = 2.56x10-2, 95% CI 1.62x10-2 to 3.63x10-2, p < 3x10-4). Limitations of the current study include the possibility that spouse effects may be due to assortative mating and the relatively small polygenic risk score effect sizes.
Genetic factors, as well as chronic pain in a partner or spouse, contribute substantially to the risk of chronic pain for an individual. Chronic pain is genetically correlated with MDD, has a polygenic architecture, and is associated with polygenic risk of MDD.
Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease characterised by progressive destruction of intrahepatic bile ducts. The strongest genetic association is with HLA-DQA1*04:01, ...but at least three additional independent HLA haplotypes contribute to susceptibility. We used dense single nucleotide polymorphism (SNP) data in 2861 PBC cases and 8514 controls to impute classical HLA alleles and amino acid polymorphisms using state-of-the-art methodologies. We then demonstrated through stepwise regression that association in the HLA region can be largely explained by variation at five separate amino acid positions. Three-dimensional modelling of protein structures and calculation of electrostatic potentials for the implicated HLA alleles/amino acid substitutions demonstrated a correlation between the electrostatic potential of pocket P6 in HLA-DP molecules and the HLA-DPB1 alleles/amino acid substitutions conferring PBC susceptibility/protection, highlighting potential new avenues for future functional investigation.
Accumulating evidence suggests that genetic factors have a role in major depressive disorder (MDD). However, only limited MDD risk loci have been identified so far. Here we perform a meta-analysis (a ...total of 90,150 MDD cases and 246,603 controls) through combing three genome-wide association studies of MDD, including 23andMe (cases were self-reported with a clinical diagnosis or treatment of depression), CONVERGE (cases were diagnosed using the Composite International Diagnostic Interview) and PGC (cases were diagnosed using direct structured diagnostic interview (by trained interviewers) or clinician-administered DSM-IV checklists). Genetic variants from two previously unreported loci (rs10457592 on 6q16.2 and rs2004910 on 12q24.31) showed significant associations with MDD (P < 5 × 10
) in a total of 336,753 subjects. SNPs (a total of 171) with a P < 1 × 10
in the meta-analysis were further replicated in an independent sample (GS:SFHS, 2,659 MDD cases (diagnosed with DSM-IV) and 17,237 controls) and one additional risk locus (rs3785234 on 16p13.3, P = 1.57 × 10
) was identified in the combined samples (a total of 92,809 cases and 263,840 controls). Risk variants on the identified risk loci were associated with gene expression in human brain tissues and mRNA expression analysis showed that FBXL4 and RSRC1 were significantly upregulated in brains of MDD cases compared with controls, suggesting that genetic variants may confer risk of MDD through regulating the expression of these two genes. Our study identified three novel risk loci (6q16.2, 12q24.31, and 16p13.3) for MDD and suggested that FBXL4 and RSRC1 may play a role in MDD. Further functional characterization of the identified risk genes may provide new insights for MDD pathogenesis.
Antidepressants demonstrate modest response rates in the treatment of major depressive disorder (MDD). Despite previous genome-wide association studies (GWAS) of antidepressant treatment response, ...the underlying genetic factors are unknown. Using prescription data in a population and family-based cohort (Generation Scotland: Scottish Family Health Study; GS:SFHS), we sought to define a measure of (a) antidepressant treatment resistance and (b) stages of antidepressant resistance by inferring antidepressant switching as non-response to treatment. GWAS were conducted separately for antidepressant treatment resistance in GS:SFHS and the Genome-based Therapeutic Drugs for Depression (GENDEP) study and then meta-analysed (meta-analysis n = 4213, cases = 358). For stages of antidepressant resistance, a GWAS on GS:SFHS only was performed (n = 3452). Additionally, we conducted gene-set enrichment, polygenic risk scoring (PRS) and genetic correlation analysis. We did not identify any significant loci, genes or gene sets associated with antidepressant treatment resistance or stages of resistance. Significant positive genetic correlations of antidepressant treatment resistance and stages of resistance with neuroticism, psychological distress, schizotypy and mood disorder traits were identified. These findings suggest that larger sample sizes are needed to identify the genetic architecture of antidepressant treatment response, and that population-based observational studies may provide a tractable approach to achieving the necessary statistical power.
Few replicable genetic associations for Major Depressive Disorder (MDD) have been identified. Recent studies of MDD have identified common risk variants by using a broader phenotype definition in ...very large samples, or by reducing phenotypic and ancestral heterogeneity. We sought to ascertain whether it is more informative to maximize the sample size using data from all available cases and controls, or to use a sex or recurrent stratified subset of affected individuals. To test this, we compared heritability estimates, genetic correlation with other traits, variance explained by MDD polygenic score, and variants identified by genome-wide meta-analysis for broad and narrow MDD classifications in two large British cohorts - Generation Scotland and UK Biobank. Genome-wide meta-analysis of MDD in males yielded one genome-wide significant locus on 3p22.3, with three genes in this region (CRTAP, GLB1, and TMPPE) demonstrating a significant association in gene-based tests. Meta-analyzed MDD, recurrent MDD and female MDD yielded equivalent heritability estimates, showed no detectable difference in association with polygenic scores, and were each genetically correlated with six health-correlated traits (neuroticism, depressive symptoms, subjective well-being, MDD, a cross-disorder phenotype and Bipolar Disorder). Whilst stratified GWAS analysis revealed a genome-wide significant locus for male MDD, the lack of independent replication, and the consistent pattern of results in other MDD classifications suggests that phenotypic stratification using recurrence or sex in currently available sample sizes is currently weakly justified. Based upon existing studies and our findings, the strategy of maximizing sample sizes is likely to provide the greater gain.