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.

[1]  Thomas W. Mühleisen,et al.  Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease , 2011, Nature Genetics.

[3]  Jakob Grove,et al.  Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection , 2018, Nature Genetics.

[4]  Ming D. Li,et al.  Genome-wide meta-analyses identify multiple loci associated with smoking behavior , 2010, Nature Genetics.

[5]  K. Karalis,et al.  Mechanisms of obesity and related pathology: linking immune responses to metabolic stress , 2009, The FEBS journal.

[6]  Simon C. Potter,et al.  The Architecture of Gene Regulatory Variation across Multiple Human Tissues: The MuTHER Study , 2011, PLoS genetics.

[7]  Benjamin D. Greenberg,et al.  Partitioning the Heritability of Tourette Syndrome and Obsessive Compulsive Disorder Reveals Differences in Genetic Architecture , 2013, PLoS genetics.

[8]  Andrew D. Johnson,et al.  Parent-of-origin specific allelic associations among 106 genomic loci for age at menarche , 2014, Nature.

[9]  P. Sullivan,et al.  Heritability and Genomics of Gene Expression in Peripheral Blood , 2014, Nature Genetics.

[10]  M. Stephens,et al.  A Statistical Framework for Joint eQTL Analysis in Multiple Tissues , 2012, PLoS genetics.

[11]  C. Molony,et al.  Genetic analysis of genome-wide variation in human gene expression , 2004, Nature.

[12]  Ross M. Fraser,et al.  Defining the role of common variation in the genomic and biological architecture of adult human height , 2014, Nature Genetics.

[13]  Todd L Edwards,et al.  Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics , 2018, Nature Communications.

[14]  A. Singleton,et al.  Genetic variability in the regulation of gene expression in ten regions of the human brain , 2014, Nature Neuroscience.

[15]  B. Penninx,et al.  Differential gene expression patterns between smokers and non‐smokers: cause or consequence? , 2015, Addiction biology.

[16]  T. Lehtimäki,et al.  Integrative approaches for large-scale transcriptome-wide association studies , 2015, Nature Genetics.

[17]  Chapel Hill.  Genome-wide meta-analyses identify multiple loci associated with smoking behavior. , 2010 .

[18]  Mary Goldman,et al.  Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics , 2016, Nature Communications.

[19]  P. Cuijpers,et al.  The Netherlands Study of Depression and Anxiety (NESDA): rationale, objectives and methods , 2008, International journal of methods in psychiatric research.

[20]  David C. Wilson,et al.  Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease , 2012, Nature.

[21]  Ivan Rusyn,et al.  An empirical Bayes approach for multiple tissue eQTL analysis , 2013, Biostatistics.

[22]  Alexander Gusev,et al.  Functional Architectures of Local and Distal Regulation of Gene Expression in Multiple Human Tissues. , 2017, American journal of human genetics.

[23]  Jun S. Liu,et al.  The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans , 2015, Science.

[24]  Wei Chen,et al.  Gene Expression in Skin and Lymphoblastoid Cells: Refined Statistical Method Reveals Extensive Overlap in Cis-eqtl Signals , 2022 .

[25]  P. Sullivan,et al.  Genome-wide association of major depression: description of samples for the GAIN Major Depressive Disorder Study: NTR and NESDA biobank projects , 2008, European Journal of Human Genetics.

[26]  Henrik,et al.  Association analyses of 249,796 individuals reveal eighteen new loci associated with body mass index , 2012 .

[27]  Jonathan P. Beauchamp,et al.  Genome-wide association study identifies 74 loci associated with educational attainment , 2016, Nature.

[28]  Simon C. Potter,et al.  Mapping cis- and trans-regulatory effects across multiple tissues in twins , 2012, Nature Genetics.

[29]  B. Bernstein,et al.  Charting histone modifications and the functional organization of mammalian genomes , 2011, Nature Reviews Genetics.

[30]  Jun S. Liu,et al.  Genetics of rheumatoid arthritis contributes to biology and drug discovery , 2013 .

[31]  Fred A. Wright,et al.  Conditional eQTL analysis reveals allelic heterogeneity of gene expression , 2017, Human molecular genetics.

[32]  Y Wang,et al.  Genetic pleiotropy between multiple sclerosis and schizophrenia but not bipolar disorder: differential involvement of immune-related gene loci , 2014, Molecular Psychiatry.

[33]  Tanya M. Teslovich,et al.  Biological, Clinical, and Population Relevance of 95 Loci for Blood Lipids , 2010, Nature.

[34]  P. Visscher,et al.  10 Years of GWAS Discovery: Biology, Function, and Translation. , 2017, American journal of human genetics.

[35]  C. Allis,et al.  The molecular hallmarks of epigenetic control , 2016, Nature Reviews Genetics.

[36]  J. Hirschhorn,et al.  Candidate genes for obesity-susceptibility show enriched association within a large genome-wide association study for BMI. , 2012, Human molecular genetics.

[37]  M. Daly,et al.  Genome-wide association studies for common diseases and complex traits , 2005, Nature Reviews Genetics.

[38]  Evan Z. Macosko,et al.  Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types , 2017, Nature Genetics.

[39]  N. Cox,et al.  Trait-Associated SNPs Are More Likely to Be eQTLs: Annotation to Enhance Discovery from GWAS , 2010, PLoS genetics.

[40]  Robert M. Maier,et al.  Causal associations between risk factors and common diseases inferred from GWAS summary data , 2017, Nature Communications.

[41]  Ellen T. Gelfand,et al.  The Genotype-Tissue Expression (GTEx) project , 2013, Nature Genetics.

[42]  Kaanan P. Shah,et al.  A gene-based association method for mapping traits using reference transcriptome data , 2015, Nature Genetics.

[43]  Nicola J. Rinaldi,et al.  Genetic effects on gene expression across human tissues , 2017, Nature.

[44]  M. Daly,et al.  LD Score regression distinguishes confounding from polygenicity in genome-wide association studies , 2014, Nature Genetics.

[45]  Adan Valladares-Salgado,et al.  Cross-tissue and tissue-specific eQTLs: partitioning the heritability of a complex trait. , 2014, American journal of human genetics.

[46]  Timothy E. Reddy,et al.  Genomic approaches for understanding the genetics of complex disease , 2015, Genome research.

[47]  Yakir A Reshef,et al.  Partitioning heritability by functional annotation using genome-wide association summary statistics , 2015, Nature Genetics.