Characterizing the Relation Between Expression QTLs and Complex Traits: Exploring the Role of Tissue Specificity

Measurement of gene expression levels and detection of eQTLs (expression quantitative trait loci) are difficult in tissues with limited sample availability, such as the brain. However, eQTL overlap between tissues might be high, which would allow for inference of eQTL functioning in the brain via eQTLs detected in readily accessible tissues, e.g. whole blood. Applying Stratified Linkage Disequilibrium Score Regression (SLDSR), we quantified the enrichment in polygenic signal of blood and brain eQTLs in genome-wide association studies (GWAS) of 11 complex traits. We looked at eQTLs discovered in 44 tissues by the Genotype-Tissue Expression (GTEx) consortium and two other large representative studies, and found no tissue-specific eQTL effects. Next, we integrated the GTEx eQTLs with regions associated with tissue-specific histone modifiers, and interrogated their effect on rheumatoid arthritis and schizophrenia. We observed substantially enriched effects of eQTLs located inside regions bearing modification H3K4me1 on schizophrenia, but not rheumatoid arthritis, and not tissue-specific. Finally, we extracted eQTLs associated with tissue-specific differentially expressed genes and determined their effects on rheumatoid arthritis and schizophrenia, these analysis revealed limited enrichment of eQTLs associated with gene specifically expressed in specific tissues. Our results pointed to strong enrichment of eQTLs in their effect on complex traits, without evidence for tissue-specific effects. Lack of tissue-specificity can be either due to a lack of statistical power or due to the true absence of tissue-specific effects. We conclude that eQTLs are strongly enriched in GWAS signal and that the enrichment is not specific to the eQTL discovery tissue. Until sample sizes for eQTL discovery grow sufficiently large, working with relatively accessible tissues as proxy for eQTL discovery is sensible and restricting lookups for GWAS hits to a specific tissue for which limited samples are available might not be advisable

Frontotemporal dementia: insights into the biological underpinnings of disease through gene co-expression network analysis

BACKGROUND: In frontotemporal dementia (FTD) there is a critical lack in the understanding of biological and molecular mechanisms involved in disease pathogenesis. The heterogeneous genetic features associated with FTD suggest that multiple disease-mechanisms are likely to contribute to the development of this neurodegenerative condition. We here present a systems biology approach with the scope of i) shedding light on the biological processes potentially implicated in the pathogenesis of FTD and ii) identifying novel potential risk factors for FTD. We performed a gene co-expression network analysis of microarray expression data from 101 individuals without neurodegenerative diseases to explore regional-specific co-expression patterns in the frontal and temporal cortices for 12 genes (MAPT, GRN, CHMP2B, CTSC, HLA-DRA, TMEM106B, C9orf72, VCP, UBQLN2, OPTN, TARDBP and FUS) associated with FTD and we then carried out gene set enrichment and pathway analyses, and investigated known protein-protein interactors (PPIs) of FTD-genes products. RESULTS: Gene co-expression networks revealed that several FTD-genes (such as MAPT and GRN, CTSC and HLA-DRA, TMEM106B, and C9orf72, VCP, UBQLN2 and OPTN) were clustering in modules of relevance in the frontal and temporal cortices. Functional annotation and pathway analyses of such modules indicated enrichment for: i) DNA metabolism, i.e. transcription regulation, DNA protection and chromatin remodelling (MAPT and GRN modules); ii) immune and lysosomal processes (CTSC and HLA-DRA modules), and; iii) protein meta/catabolism (C9orf72, VCP, UBQLN2 and OPTN, and TMEM106B modules). PPI analysis supported the results of the functional annotation and pathway analyses. CONCLUSIONS: This work further characterizes known FTD-genes and elaborates on their biological relevance to disease: not only do we indicate likely impacted regional-specific biological processes driven by FTD-genes containing modules, but also do we suggest novel potential risk factors among the FTD-genes interactors as targets for further mechanistic characterization in hypothesis driven cell biology work

Gene co-expression networks shed light into diseases of brain iron accumulation

Aberrant brain iron deposition is observed in both common and rare neurodegenerative disorders, including those categorized as Neurodegeneration with Brain Iron Accumulation (NBIA), which are characterized by focal iron accumulation in the basal ganglia. Two NBIA genes are directly involved in iron metabolism, but whether other NBIA-related genes also regulate iron homeostasis in the human brain, and whether aberrant iron deposition contributes to neurodegenerative processes remains largely unknown. This study aims to expand our understanding of these iron overload diseases and identify relationships between known NBIA genes and their main interacting partners by using a systems biology approach. We used whole-transcriptome gene expression data from human brain samples originating from 101 neuropathologically normal individuals (10 brain regions) to generate weighted gene co-expression networks and cluster the 10 known NBIA genes in an unsupervised manner. We investigated NBIA-enriched networks for relevant cell types and pathways, and whether they are disrupted by iron loading in NBIA diseased tissue and in an in vivo mouse model. We identified two basal ganglia gene co-expression modules significantly enriched for NBIA genes, which resemble neuronal and oligodendrocytic signatures. These NBIA gene networks are enriched for iron-related genes, and implicate synapse and lipid metabolism related pathways. Our data also indicates that these networks are disrupted by excessive brain iron loading. We identified multiple cell types in the origin of NBIA disorders. We also found unforeseen links between NBIA networks and iron-related processes, and demonstrate convergent pathways connecting NBIAs and phenotypically overlapping diseases. Our results are of further relevance for these diseases by providing candidates for new causative genes and possible points for therapeutic intervention

Hypermethylation in the ZBTB20 gene is associated with major depressive disorder.

This is the final version of the article. Available from BioMed Central via the DOI in this record.BACKGROUND: Although genetic variation is believed to contribute to an individual's susceptibility to major depressive disorder, genome-wide association studies have not yet identified associations that could explain the full etiology of the disease. Epigenetics is increasingly believed to play a major role in the development of common clinical phenotypes, including major depressive disorder. RESULTS: Genome-wide MeDIP-Sequencing was carried out on a total of 50 monozygotic twin pairs from the UK and Australia that are discordant for depression. We show that major depressive disorder is associated with significant hypermethylation within the coding region of ZBTB20, and is replicated in an independent cohort of 356 unrelated case-control individuals. The twins with major depressive disorder also show increased global variation in methylation in comparison with their unaffected co-twins. ZBTB20 plays an essential role in the specification of the Cornu Ammonis-1 field identity in the developing hippocampus, a region previously implicated in the development of major depressive disorder. CONCLUSIONS: Our results suggest that aberrant methylation profiles affecting the hippocampus are associated with major depressive disorder and show the potential of the epigenetic twin model in neuro-psychiatric disease.The study was funded by the Wellcome Trust; European Community’s Seventh Framework Programme (FP7/2007-2013). The study also receives support from the National Institute for Health Research (NIHR) Clinical Research Facility at Guy’s & St Thomas’ Davies et al. Genome Biology 2014, 15:R56 Page 9 of 12 http://genomebiology.com/2014/15/4/R56 NHS Foundation Trust and NIHR Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust and King's College London. Matthew Davies is supported by the EU FP7 grant EuroBATS (No. 259749). Tim Spector is an NIHR senior Investigator and is holder of an ERC Advanced Principal Investigator award. Further funding support for this project was obtained from the European Research Council (project number 250157). The members of the UK Brain Expression Consortium (UKBEC) are: (1) Department of Molecular Neuroscience, UCL Institute of Neurology, London, UK: John A Hardy, Mina Ryten, and Daniah Trabzuni; (2) Department of Medical and Molecular Genetics, King's College London, UK: Michael E Weale, Adaikalavan Ramasamy and Paola Forabosco; (3) Department of Pathology, The University of Edinburgh, Wilkie Building, Teviot Place, Edinburgh, UK: Colin Smith and Robert Walker. Australia: funding for phenotype and blood collection was from NHMRC grants to Nick Martin and NIH grants to Andrew Heath and Pamela Madden. We thank David Smyth for database management, Lisa Bowdler for sample preparation, and the twins for their cooperation

Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption.

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91 462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10-8).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee

Epilepsy, hippocampal sclerosis and febrile seizures linked by common genetic variation around SCN1A

Epilepsy comprises several syndromes, amongst the most common being mesial temporal lobe epilepsy with hippocampal sclerosis. Seizures in mesial temporal lobe epilepsy with hippocampal sclerosis are typically drug-resistant, and mesial temporal lobe epilepsy with hippocampal sclerosis is frequently associated with important co-morbidities, mandating the search for better understanding and treatment. The cause of mesial temporal lobe epilepsy with hippocampal sclerosis is unknown, but there is an association with childhood febrile seizures. Several rarer epilepsies featuring febrile seizures are caused by mutations in SCN1A, which encodes a brain-expressed sodium channel subunit targeted by many anti-epileptic drugs. We undertook a genome-wide association study in 1018 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 7552 control subjects, with validation in an independent sample set comprising 959 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 3591 control subjects. To dissect out variants related to a history of febrile seizures, we tested cases with mesial temporal lobe epilepsy with hippocampal sclerosis with (overall n = 757) and without (overall n = 803) a history of febrile seizures. Meta-analysis revealed a genome-wide significant association for mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures at the sodium channel gene cluster on chromosome 2q24.3 [rs7587026, within an intron of the SCN1A gene, P = 3.36 × 10(-9), odds ratio (A) = 1.42, 95% confidence interval: 1.26-1.59]. In a cohort of 172 individuals with febrile seizures, who did not develop epilepsy during prospective follow-up to age 13 years, and 6456 controls, no association was found for rs7587026 and febrile seizures. These findings suggest SCN1A involvement in a common epilepsy syndrome, give new direction to biological understanding of mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures, and open avenues for investigation of prognostic factors and possible prevention of epilepsy in some children with febrile seizures

Identification of expression quantitative trait loci associated with schizophrenia and affective disorders in normal brain tissue.

Schizophrenia and the affective disorders, here comprising bipolar disorder and major depressive disorder, are psychiatric illnesses that lead to significant morbidity and mortality worldwide. Whilst understanding of their pathobiology remains limited, large case-control studies have recently identified single nucleotide polymorphisms (SNPs) associated with these disorders. However, discerning the functional effects of these SNPs has been difficult as the associated causal genes are unknown. Here we evaluated whether schizophrenia and affective disorder associated-SNPs are correlated with gene expression within human brain tissue. Specifically, to identify expression quantitative trait loci (eQTLs), we leveraged disorder-associated SNPs identified from 11 genome-wide association studies with gene expression levels in post-mortem, neurologically-normal tissue from two independent human brain tissue expression datasets (UK Brain Expression Consortium (UKBEC) and Genotype-Tissue Expression (GTEx)). Utilizing stringent multi-region meta-analyses, we identified 2,224 cis-eQTLs associated with expression of 40 genes, including 11 non-coding RNAs. One cis-eQTL, rs16969968, results in a functionally disruptive missense mutation in CHRNA5, a schizophrenia-implicated gene. Importantly, comparing across tissues, we find that blood eQTLs capture 30% of brain-associated eQTLs are significant in tibial nerve. This study identifies putatively causal genes whose expression in region-specific tissue may contribute to the risk of schizophrenia and affective disorders

Dystonia genes functionally converge in specific neurons and share neurobiology with psychiatric disorders

Dystonia is a neurological disorder characterized by sustained or intermittent muscle contractions causing abnormal movements and postures, often occurring in absence of any structural brain abnormality. Psychiatric comorbidities, including anxiety, depression, obsessive-compulsive disorder and schizophrenia, are frequent in patients with dystonia. While mutations in a fast-growing number of genes have been linked to Mendelian forms of dystonia, the cellular, anatomical, and molecular basis remains unknown for most genetic forms of dystonia, as does its genetic and biological relationship to neuropsychiatric disorders. Here we applied an unbiased systems-biology approach to explore the cellular specificity of all currently known dystonia-associated genes, predict their functional relationships, and test whether dystonia and neuropsychiatric disorders share a genetic relationship. To determine the cellular specificity of dystonia-associated genes in the brain, single-nuclear transcriptomic data derived from mouse brain was used together with expression-weighted cell-type enrichment. To identify functional relationships among dystonia-associated genes, we determined the enrichment of these genes in co-expression networks constructed from 10 human brain regions. Stratified linkage-disequilibrium score regression was used to test whether co-expression modules enriched for dystonia-associated genes significantly contribute to the heritability of anxiety, major depressive disorder, obsessive-compulsive disorder, schizophrenia, and Parkinson’s disease. Dystonia-associated genes were significantly enriched in adult nigral dopaminergic neurons and striatal medium spiny neurons. Furthermore, 4 of 220 gene co-expression modules tested were significantly enriched for the dystonia-associated genes. The identified modules were derived from the substantia nigra, putamen, frontal cortex, and white matter, and were all significantly enriched for genes associated with synaptic function. Finally, we demonstrate significant enrichments of the heritability of major depressive disorder, obsessive-compulsive disorder and schizophrenia within the putamen and white matter modules, and a significant enrichment of the heritability of Parkinson’s disease within the substantia nigra module. In conclusion, multiple dystonia-associated genes interact and contribute to pathogenesis likely through dysregulation of synaptic signalling in striatal medium spiny neurons, adult nigral dopaminergic neurons and frontal cortical neurons. Furthermore, the enrichment of the heritability of psychiatric disorders in the co-expression modules enriched for dystonia-associated genes indicates that psychiatric symptoms associated with dystonia are likely to be intrinsic to its pathophysiology

Pathological relationships involving iron and myelin may constitute a shared mechanism linking various rare and common brain diseases.

We previously demonstrated elevated brain iron levels in myelinated structures and associated cells in a hemochromatosis Hfe (-/-) xTfr2 (mut) mouse model. This was accompanied by altered expression of a group of myelin-related genes, including a suite of genes causatively linked to the rare disease family 'neurodegeneration with brain iron accumulation' (NBIA). Expanded data mining and ontological analyses have now identified additional myelin-related transcriptome changes in response to brain iron loading. Concordance between the mouse transcriptome changes and human myelin-related gene expression networks in normal and NBIA basal ganglia testifies to potential clinical relevance. These analyses implicate, among others, genes linked to various rare central hypomyelinating leukodystrophies and peripheral neuropathies including Pelizaeus-Merzbacher-like disease and Charcot-Marie-Tooth disease as well as genes linked to other rare neurological diseases such as Niemann-Pick disease. The findings may help understand interrelationships of iron and myelin in more common conditions such as hemochromatosis, multiple sclerosis and various psychiatric disorders

Dystonia genes functionally converge in specific neurons and share neurobiology with psychiatric disorders

Dystonia is a neurological disorder characterized by sustained or intermittent muscle contractions causing abnormal movements and postures, often occurring in absence of any structural brain abnormality. Psychiatric comorbidities, including anxiety, depression, obsessive-compulsive disorder and schizophrenia, are frequent in dystonia patients. While mutations in a fast-growing number of genes have been linked to Mendelian forms of dystonia, the cellular, anatomical, and molecular basis remains unknown for most genetic forms of dystonia, as does its genetic and biological relationship to neuropsychiatric disorders. Here we applied an unbiased systems-biology approach to explore the cellular specificity of all currently known dystonia-associated genes, predict their functional relationships, and test whether dystonia and neuropsychiatric disorders share a genetic relationship. To determine the cellular specificity of dystonia-associated genes in the brain, single-nuclear transcriptomic data derived from mouse brain was used together with expression-weighted cell-type enrichment. To identify functional relationships amongst dystonia-associated genes, we determined the enrichment of these genes in co-expression networks constructed from ten human brain regions. Stratified linkage-disequilibrium score regression was used to test whether co-expression modules enriched for dystonia-associated genes significantly contribute to the heritability of anxiety, major depressive disorder, obsessive-compulsive disorder, schizophrenia, and Parkinson’s disease. Dystonia-associated genes were significantly enriched in adult nigral dopaminergic neurons and striatal medium spiny neurons. Furthermore, four of the 220 gene co-expression modules tested were significantly enriched for the dystonia-associated genes. The identified modules were derived from the substantia nigra, putamen, frontal cortex, and white matter, and were all significantly enriched for genes associated with synaptic function. Finally, we demonstrated significant enrichments of the heritability of major depressive disorder, obsessive-compulsive disorder and schizophrenia within the putamen and white matter modules, and a significant enrichment of the heritability of Parkinson’s disease within the substantia nigra module. In conclusion, multiple dystonia-associated genes interact and contribute to pathogenesis likely through dysregulation of synaptic signalling in striatal medium spiny neurons, adult nigral dopaminergic neurons and frontal cortical neurons. Furthermore, the enrichment of the heritability of psychiatric disorders in the co-expression modules enriched for dystonia associated genes indicates that psychiatric symptoms associated with dystonia are likely to be intrinsic to its pathophysiology