Email updates

Keep up to date with the latest news and content from BMC Medical Genomics and BioMed Central.

Open Access Highly Accessed Research article

Overlap of expression Quantitative Trait Loci (eQTL) in human brain and blood

Marna McKenzie1, Anjali K Henders2, Anthony Caracella2, Naomi R Wray1 and Joseph E Powell1*

Author Affiliations

1 The University of Queensland, Queensland Brain Institute, Brisbane, Queensland, Australia

2 Queensland Institute of Medical Research, Brisbane, Herston, Australia

For all author emails, please log on.

BMC Medical Genomics 2014, 7:31  doi:10.1186/1755-8794-7-31


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1755-8794/7/31


Received:3 September 2013
Accepted:20 May 2014
Published:3 June 2014

© 2014 Mckenzie et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Abstract

Background

Expression quantitative trait loci (eQTL) are genomic regions regulating RNA transcript expression levels. Genome-wide Association Studies (GWAS) have identified many variants, often in non-coding regions, with unknown functions and eQTL provide a possible mechanism by which these variants may influence observable phenotypes. Limited access and availability of tissues such as brain has led to the use of blood as a substitute for eQTL analyses.

Methods

Here, we evaluate the overlap of eQTL reported in published studies conducted in blood and brain tissues to assess the utility of blood as an alternative to brain tissue in the study of neurological and psychiatric conditions. Expression QTL results from eight published brain studies were compared to blood eQTL identified in from a meta-analysis involving 5,311 individuals. We accounted for differences in SNP platforms and study design by using SNP proxies in high linkage disequilibrium with reported eQTL. The degree of overlap between studies was calculated by ascertaining if an eQTL identified in one study was also identified in the other study.

Results

The percentage of eQTL overlapping for brain and blood expression after adjusting for differences in sample size ranged from 13 - 23% (mean 19.2%). Amongst pairs of brain studies eQTL overlap ranged from 0 - 35%, with higher degrees of overlap found for studies using expression data collected from the same brain region.

Conclusion

Our results suggest that whenever possible tissue specific to the pathophysiology of the disease being studied should be used for transcription analysis.

Background

The combination of expression profiling and genotyping from microarrays has led to estimation of correlations between genetic variants such as Single Nucleotide Polymorphisms (SNPs) and RNA transcript expression levels [1-3]. Expression quantitative trait loci (eQTL) are the genomic loci that influence RNA transcript expression levels. The contribution of eQTL to underlying variation in complex traits such as disease susceptibility has been studied extensively since expression levels were first described as quantitative traits with a genetic basis. Gene expression is the primary detectable phenotype in the development of complex traits; therefore expression levels act as an intermediate phenotype between genetic architecture and observable multifactorial traits such as common diseases. It has been shown that SNPs associated with common diseases identified by genome-wide association studies (GWAS) are enriched for expression-affecting SNPs (eSNPs) [4,5]. Furthermore, most associated variants are not located in protein coding regions and are instead highly enriched for regulatory regions of the genome [6], suggesting that for many variants, the functional mechanism by which they affect disease susceptibility is through gene regulation. The majority of eQTL that have currently been identified are located within the cis-region (typically +/− 1 MB) of the transcription start site (TSS). Expression-QTL located outside of the cis-region are termed trans-eQTL and likely reflect an indirect relationship of the SNP on gene expression. Cis-eQTL tend to have larger effect sizes [6], are more often replicated across multiple studies [7] and more likely to reflect a direct functional relationship between the SNP and the measured expression levels.

The pathophysiology of most common diseases is restricted to a limited number of tissue types or organ systems. Therefore, to understand the mechanisms of disease susceptibility and develop preventative and targeted therapies, we ultimately require knowledge of genetic control of regulatory variation in many different tissues. A common limitation in disease genomics studies is availability of pathologically relevant tissue on which to measure expression. Most studies instead rely on inferences drawn from readily available (‘proxy’) tissues, typically whole blood. However, the genetic control of RNA transcription is known to vary between tissues [5,8-11], and the phenotypic correlations and co-expression can be low [12]. Expression-SNPs that are also associated with disease susceptibility are more likely to affect expression levels in a tissue-specific manner compared to eSNPs with no known disease association [5]. Due to mainly study-design limitations most eQTL analyses focus on gene expression measured in a single tissue type [9,13-18], posing limitations when attempting to study more inaccessible tissues, such as the brain.

Using expression levels measured from brain tissue for eQTL analysis poses several unique problems. First, the brain is a collection of numerous cell types meaning that gene expression levels are not consistent throughout the entire structure [7,19]. This introduces cell type specificity within the brain itself. Secondly, collection of tissue samples in post-mortem requires consideration of several novel factors not relevant when examining other tissue types – such as cause of death, post-mortem RNA integrity and post-mortem interval [20]. This has led to the search for alternate methods for evaluation of eQTL within the brain.

One approach is to use expression levels measured from more accessible tissues, such as blood, as a substitute to gain insight into gene regulation underlying brain-related traits such as psychiatric and neurological disorders [21-25]. When such approaches are used an important question is: To what degree is the genetic control for expression the same in brain and blood tissues? Difficulties in the collection of expression data from both blood and brain tissue for the same individuals mean methods are required to assess the overlap of genetic control. Using results from published studies represents a practical solution. Here we report a comprehensive evaluation of the degree of eQTL overlap between blood and brain and aim to determine the validity of using blood as a surrogate for brain tissue. We have sought to account for differences in expression and genotype platform arrays as well as to correct for differences in sample size and thus differences in statistical power between studies.

Methods

eQTL results

Published eQTL results and supporting information were obtained from eight eQTL studies carried out on expression levels measured in the brain using high-throughput arrays [7,10,26-31]. Full details of each study are given in Additional file 1: Table S1. Each of these studies evaluated the association between transcript probe expression levels and SNP genotype using linear regression models. Significant eQTL were determined at a p-value threshold that was specific to the study. Studies used different metrics and protocols to identify eQTL and so we sought to standardize results between studies by producing lists of eQTL and eSNPs that reached the study-wide significance level of α  =  0.05. The number of significant eQTL reported ranged from 52 to 2,975, which represents between <1 - 19% of the probes analyzed in each study. This range most likely reflects differences in sample size, protocols and statistical analyses between the studies. Three studies [10,26,30] reported cis associations only, while the remaining five reported both cis and trans associations (trans eQTL are those located outside of the cis-region definition of the study), including SNPs distal to the associated transcript [7,27-29,31]. Details of the significant eQTL and eSNPs reported in each study are given in Table 1.

Additional file 1: Table S1. Complete summary of each study included in comparison.

Format: DOCX Size: 87KB Download fileOpen Data

Table 1. Summary of each study included in the comparison

All brain tissue samples were collected post-mortem. Some studies reported eQTL for expression collected in multiple brain regions or from samples with different neuropathologies. Where possible, eQTL overlap was evaluated for separate tissues and neuropathologies, however, most studies pooled results from different tissues or from individuals with different neurological/psychiatric conditions. Therefore, separate analysis was not always possible. For example, two of the studies included samples from individuals diagnosed with schizophrenia, bipolar disorder, or major depressive disorder [27,28], whilst others pooled samples with pathologies such as Alzheimer’s disease and progressive supranuclear palsy [29,30]. Only four studies used expression data collected from a single region of the brain [10,26,28,30], whilst the remainder, with the exception of Gibbs et al. [7], combined data acquired from several brain regions.

The eQTL results from the eight brain studies were compared to results from a large meta-analysis of eQTL in for expression levels measured in whole blood [32]. The Westra et al. analysis is the largest eQTL study published to date and comprises of a meta-analysis of 5,311 individuals in a discovery phase and 2,775 individuals in a replication phase. All expression levels were measured in peripheral blood and corrected for batch effects and cell counts. SNP effects were estimated using a weighted Z-score of the beta values calculated using a linear regression model in each study cohort. Significance was determined using permutations and a study-wide threshold of false discovery rate (FDR) of 0.05. Cis-eQTL were identified for a total of 4,909 genes (33% of all tested) and trans-eQTL for 430 genes (3%).

This study used previously published data. The research was approved by the University of Queensland Human Ethics Review Board and the QIMR Berghofer Medical Research Institute Institutional Review Board for Research on Human Subjects.

eQTL overlap

The eQTL studies included in this analysis used a variety of high-throughput expression arrays (Additional file 1: Table S1). We sought to overcome the degree of probe overlap between by comparing eQTL for probes that tag the same gene exon. We evaluated the overlap of eQTL identified in blood [32] and each of the brain studies as well as between each pair of brain studies. The latter allows us to draw conclusions on eQTL overlap for different brain regions or tissues as well as potential impact of psychiatric or neurological disorders. The following procedures were followed to ascertain the overlap of eQTL between peripheral blood and each brain study and between each pair-wise comparison of brain studies.

The list of significant (study-wide α = 0.05 eQTL for studies i and j were first assessed for probes containing an eQTL in both studies. For each eQTL present in studies i and j, overlapping eQTL were determined using two approaches. Firstly, SNPs with the smallest p-value (eSNPs) from eQTL were compared based on their recorded presence in both studies. This provides a basic comparison of all overlapping eQTL that can be identified from standard reported results. However, the mean overlap between i and j is expected to be biased downwards due to differences in SNP array platforms and inconsistencies in linkage disequilibrium (LD) between study populations. We accounted for this by applying a second approach that used known SNP proxies [33] from the 1000 Genomes study (CEU panel). SNP proxies were defined as SNPs in high LD (r2 > 0.8) and within 100 kb of the original genotyped eSNPs. A list of SNP proxies for each eSNP was generated for all studies. The SNP proxies for eQTL in studies i and j were then compared to find overlapping proxy eQTL (hereon termed eQTLproxy). This approach provides a realistic estimate of the true eQTL overlap that would be expected if all studies had used the same genotyping platform.

The power to detect eQTL is partly a function of sample size. Thus, comparison of eQTL between two studies of different sample sizes is expected to downwardly bias the estimate of the true proportion overlapping due to the lower power of the smaller study. We addressed this problem by using a method presented by Ding et al. [32] that adjusts for the difference in samples sizes, and thus statistical power, to re-estimate eQTL overlap expected under the scenario of equal sample size [34]. Following the notation of Ding et al., the estimate of power-adjusted overlap, we start with the lists of eQTL identified in studies i and j. EQTL are determined based on significance threshold, α, which corresponds to a false discovery rate (FDR) of 0.05. From these lists we determined the observed percentage of eQTL overlapping between studies (πraw). Because both eQTL lists are inevitably incomplete, πraw will be an underestimate of the true level of overlap, π. If study j is of smaller sample size than i, then the power to detect ‘true’ overlapping eQTL will be lower. The power-adjusted expected overlap in significant eQTL (<a onClick="popup('http://www.biomedcentral.com/1755-8794/7/31/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1755-8794/7/31/mathml/M1">View MathML</a>) is calculated as [32]

<a onClick="popup('http://www.biomedcentral.com/1755-8794/7/31/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1755-8794/7/31/mathml/M2">View MathML</a>

(1)

where powerj is determined from the effect sizes of the overlapping eQTL, the sample size and aj,

<a onClick="popup('http://www.biomedcentral.com/1755-8794/7/31/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1755-8794/7/31/mathml/M3">View MathML</a>

(2)

powerj_raw is defined as the statistical power of study i with a sample size matching study j to detect all identified eQTL from study i. Since the blood-based eQTL analysis was the largest study by far, we used this sample to estimate powerj_raw by:

1. From the Westra et al. list of eQTL we sampled two mutually exclusive and independent datasets (i1 and i2), with the sample size of both i1 and i2 equal to that of study j.

2. EQTL were identified in datasets i1 and i2, assuming significance threshold of study-wide FDR of 0.05.

3. powerj_raw was calculated as the proportion of eQTL identified in i1 that are also significant in i2.

4. <a onClick="popup('http://www.biomedcentral.com/1755-8794/7/31/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1755-8794/7/31/mathml/M4">View MathML</a> was calculated (equations 1 and 2) between each brain study and the Westra et al. blood eQTL study.

Functional characterization of overlapping eQTL

A KEGG pathway enrichment analysis was used to assess shared biological functionality amongst genes with either overlapping or independent eQTL. From each pairwise comparison of blood and brain eQTL, genes listed as having either overlapping or non-overlapping eQTL were analyzed for shared functionality using DAVID [35]. Significance of pathway enrichment was determined from a modified Fisher’s test, which represents the probability that a set of genes of related terms are presented at a given percentage in the list. Multiple testing was accounted for using a Benjamini-Hochberg FDR of 0.05 [36].

Results

Blood and brain eQTL overlap

For cis-acting eQTL only, the overlap of eQTL between blood and brain studies ranged 0.1–14.2% (mean 4.3%) (Table 2). In total 4,909 cis-eQTL were found for blood gene expression, of which 156 were also identified for cerebellum and temporal cortex expression [29]. Overlapping trans-eQTL were identified from 5 studies, but the number of genes was low (2–35) (Table 2). Trans-eQTL typically have smaller effect sizes than cis-eQTL and suffer from a greater multiple testing burden, resulting in lower power for detection compared to cis-effects. These estimates of overlap are likely to be an under-estimate of the true level of overlap as the ‘replication’ studies in brain tissue have smaller sample sizes and thus lower power than Westra et al. [32]. Indeed, if we use the brain studies as the ‘discovery’ and Westra et al. as ‘replication’ the proportion of cis-eQTL overlap ranges 5.7-70% (mean 22%). These results are from analyses where the probes are matched to gene exons and using SNP proxies of eSNPs. Verification of the robustness of our results to deviation in the SNP proxy LD threshold is shown in Additional file 2: Table S2. We believe our results provide an estimate of the proportion of overlap that is as accurate as possible given the limitations of the published data. Using the method published by Ding et al. [34] we attempted to account for differences in the sample sizes between the blood and brain studies (see Methods). Accounting for the smaller sample sizes of brain studies, the expected proportion of blood eQTL that would be identified had the brain studies be of equal size (n = 862), <a onClick="popup('http://www.biomedcentral.com/1755-8794/7/31/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1755-8794/7/31/mathml/M5">View MathML</a>, ranges 13% -23% (mean = 19.2%) (Table 3), suggesting that observed estimates (πraw) are under-estimates due to lower sample sizes in brain studies. It is important to note that the method used to estimate the adjusted overlap assumes that the distributions of eQTL effect sizes are similar for overlapping and non-overlapping eQTL [34]. There is evidence to suggest this assumption may be incorrect and the implications of this are discussed below.

Table 2. Summary of blood eQTL overlap with each brain eQTL study

Additional file 2: Table S2. Overlap of eQTL was evaluated using the SNP proxy lists generated under two linkage disequilibrium threshold (r2 = 0.8 and r2 = 0.5).

Format: DOCX Size: 35KB Download fileOpen Data

Table 3. Estimation of the expected degree of overlap between blood eQTL and each of the brain studies should the sample sizes be equal

Functional characterization

To investigate if both overlapping and non-overlapping eQTL have a shared functionality we performed a pathway analysis using DAVID [35]. The lists of overlapping and non-overlapping eQTL genes are given in Additional file 3: Table S4. For both sets of genes with overlapping and non-overlapping eQTL no statistically significant enrichment was found for functional annotation or Gene Ontology (GO) terms. This suggests that the genes with overlapping eQTL have similar functional roles in brain and blood cell types and tissues. Amongst the eight comparisons, 268 independent genes were found to have overlapping cis-eQTL with blood. Of these 55 (21%) had an eQTL observed in more than two brain studies. The genes that most consistently reported an overlapping eQTL were NSFL1C and PEX6 (reported in six of the brain eQTL studies), and CDK5RAP2, CDS2, CHURC1, CRIPT, HMBOX1, MRPL43, NAPRT1, NSUN2, RABEP1, ZNF266, ABHD12, and PILRB (reported in five of the brain eQTL studies) (Additional file 4: Table S3).

Additional file 3: Table S4. Gene lists with overlapping eQTL found for whole blood. Overlapping eQTL genes are given separately for each brain study.

Format: XLSX Size: 22KB Download fileOpen Data

Additional file 4: Table S3. Genes with an overlapping eQTL reported in Westra et al. and brain studies.

Format: DOCX Size: 39KB Download fileOpen Data

Brain region eQTL overlap

The brain studies examined here report eQTL for expression levels measured in cells collected from several different tissues and brain regions. Given specific differences in the aetiology of many neurological disorders, one important consideration is whether eQTL overlap between different brain regions. To investigate this we analyzed overlap using our observed and SNP proxy methods for each pairwise combination of brain studies. Cis-eQTL overlap is shown in Table 4, while the only observable trans-eQTL overlap was in Myers et al. [31] and Webster et al. [29] at 0.2% and 1.5% respectively (including SNP proxies). The variability in overlap between studies remains high even when comparing the eQTL data generated from the same brain region. For example, in the two studies examining the prefrontal cortex, the proportion of overlapping genes with a significant eQTL was only 1% to 7% (Table 4) [26,28]. However, the proportion of overlapping genes with a significant cis-eQTL from cerebellar tissues in Gibbs et al. [7] and Zou et al. [30] was much greater. Of the genes from Gibbs et al. [7] with a significant cis-eQTL in the cerebellum, 39% were also observed in Zou et al. [30]. Conversely, 8% of the genes with a significant cis-QTL from cerebellar tissue in Zou et al. [30] were reported in Gibbs et al. [7].

Table 4. Overlap of eQTL from each of the pairwise comparisons of brain studies

To evaluate eQTL overlap between different brain regions we calculated the average overlap for studies that used expression levels measured from the same region compared to the average for studies using expression measured in different regions. For each comparison we chose the dataset with the largest sample size as the ‘discovery’ sample (Table 2). The mean overlap between studies with the same tissue is 20% and for studies using different tissues is 12%. Following this, we chose pairs of studies that both collected samples from healthy or normal neuropathology and compared those to studies using healthy versus diseased patients. The mean overlap between healthy-healthy studies was 19% and between healthy-disease was 13% (Figure 1).

thumbnailFigure 1. Venn diagram showing the overlap in genes with common eQTL between brain studies. In figure (A) all studies used samples with healthy or normal neuropathology, and tissue samples were collected from different cortical brain regions; Cortex (Myers), Prefrontal Cortex (Colantuoni) and Temporal Cortex (Gibbs). In (B) samples are collected from the Temporal Cortex and Cerebellum in individuals with normal or healthy neuropathology (Zou) and Cortex in individuals with late onset Alzheimer’s disease (Webster).

Of the 765 genes found to have a cis-eQTL in one of the two studies involving cerebellar tissue [7,30], 137 (18%) were also reported to have a cis-eQTL in blood. This proportion increased considerably when evaluating eQTL reported by both cerebellar studies – 25 of the 68 (37%) genes reported in both cerebellar studies also displayed a significant cis-eQTL association in blood.

Through contrasting distinctive neuropathologies, Liu et al. [28] showed that different psychiatric conditions had little effect on eQTL mapping compared to healthy controls. Consequently, most studies included here had pooled eQTL results from cohorts where samples have differing psychiatric conditions [27,28]. However, Webster et al. [29], reported eQTL that had a significant effect only amongst individuals with late-onset Alzheimer’s. The overlap of these eQTL with those found in blood is summarized in Table 5. Notably, blood eQTL overlap with Alzheimer’s disease interacting eQTL was three times lower than eQTL without a significant interaction with disease status.

Table 5. Summary of Westra et al.[32]blood eQTL overlap with eQTL from Webster et al.[29]which were found to have an interaction with late-onset Alzheimer’s disease (AD) status as well as those independent of disease status

Discussion

We have quantified the overlap in eQTL identified from eight brain studies with those found in blood. A key motivation of this work was to provide information on the utility of using blood as a surrogate for brain tissue when trying to understand the role of gene expression in neurological and psychiatric conditions. Previous work had shown that the transcriptome organization is poorly preserved between brain and blood, although the conservation was stronger for certain hub genes [12].

If we solely consider the overlap amongst the eight brain studies, there is relatively low, but highly significant, levels of replication (p = 1.4e−8). The greatest degree of overlap (35%) of cis-eQTL is between two studies that focused on cortex regions [28,30]. If trans-eQTL are included in our analysis then the percentage overlap falls to 15%. The low level of overlap is likely due to a number of contributing factors. Firstly, we are comparing individuals from different population demographics whose expression levels have been measured using different laboratory procedures and protocols. Difference in expression profiles as a result of these variables is well documented [26,37,38]. Secondly, individuals from different disease cohorts were included together in the analyses. The effect of disease status on eQTL is still unclear; Liu et al. [28] found that psychiatric disease status had minimal effects on eQTL mapping, whilst Webster et al. [29] showed that most of the Alzheimer’s disease associated eQTL effects were not found in healthy controls. Thirdly, the studies sampled and often combined data from different regions of the brain. Yet it has been shown that expression profiles generated from different brain regions are distinct [7,19]. Additional and larger studies are needed to clarify and quantify the true overlap.

Where possible we used approaches to account for differences in microarray and genotype platforms that could contribute to a low concordance in eQTL overlap between studies. EQTL and eSNP were matched for expression probes located within the same exon and SNP proxies were generated and analyzed for all eSNPs. The proxy SNPs included all known common SNPs within 100 kb region that are linked to the eQTL by LD. Therefore, if the eQTL was present in each study, it should be detected using this method.

Brain eQTL studies, with smaller sample sizes, lack statistical power compared to the blood meta-analysis [32]. To address this, we employed an approach [32] which re-estimates the expected degree of eQTL overlap, should the sample size of the brain study match that of Westra et al. Although it has been estimated that a sample size of 100 individuals is sufficient for 80% power in eQTL studies [39], the power to detect an effect across multiple studies is reduced when an eSNP has a small effect size or when multiple loci are controlling transcript expression levels.

The highest proportion of cis-eQTL overlap between blood and any brain study was 14.2% [7]. Under a null hypothesis of no eQTL overlap, these observed values are much greater than would be expected by chance (p = 2.7e−11). We can calculate the expected amount of overlap under the null hypothesis that there is no true overlap using the proportion of probes with identified eQTL in each study. For example, the amount of overlap expected between Westra et al. [32] and Gibbs et al. [7] based on chance would be equal to the percentage of probes analyzed in Westra et al. that were found to have an eQTL (33%), multiplied by the percentage of probes examined in Gibbs et al. [7] that were found to have an eQTL (5%), which is 1.6%.

Consistent with findings from other tissues, the genetic regulation of transcript expression levels within the brain and blood appears to be largely tissue specific [8,10]. Our results presented here use data collected among different individuals and thus will underestimate the true overlap because of inconsistencies in study design, environmental effects and allele frequencies of SNPs. Studies such as GTEX [40], which is collecting samples from multiple tissues in the same individuals, can be used to provide a more complete understanding of identifiable eQTL overlap between these tissues. Importantly, it would also allow an accurate quantification of the direction of allelic effects between tissues. Here, due to unavailability of beta estimates, we have assumed the direction of eQTL effects between tissues is the same. As is expected, the overlap between brain and blood was consistently lower than the amount of overlap among brain eQTL studies. This provides support for the methods used for eQTL comparison in this study. The high level of overlap identified between studies using samples with disparate medical history provides support for the hypothesis that neurological conditions have little impact on eQTL effect estimates [28].

Table 2 provides information on the genes with an eQTL in both brain and blood, presented as a proportion of eQTL identified in the studies. The proportion of eQTL overlap was on average greater when the brain eQTL study was used as the ‘discovery’ cohort. This is likely due to the larger sample size of the blood eQTL study.

The mean genomic distance of overlapping blood–brain cis-eSNP from the probe TSS was 38.5 kb (not including proxy SNPs) – substantially less than the 121 kb average distance for brain eQTL reported by Gibbs et al. [7]. This is in agreement with previous findings, which show that eQTL found in multiple tissues tended to localize closer to the TSS than tissue-specific eQTL [9]. It is also known that eQTL found in multiple tissues have larger effect sizes than average, which could lead to a upwards bias in the reported percentage of overlapping eQTL between brain and blood [5].

The pathophysiology of many neurological and psychiatric conditions is often localized to specific brain regions. For example, prefrontal and temporal cortex abnormalities have been repeatedly associated with schizophrenia [41-43] while affective disorders such as bipolar disorder and major depressive disorder have been linked with dysfunction in several brain regions – cingulate cortex, amygdala, thalamus, hippocampus and the frontal lobe [44-48]. Similarly, neurological disorders tend to act within particular brain structures – degeneration in the substantia nigra is a hallmark of Parkinson’s disease [49] and one of the primary characteristics of Alzheimer’s disease is the prevalence of neurofibrillary tangles and amyloid plaques, particularly in medial temporal lobe structures [50,51]. Thus, analysis of regions specifically affected by the disorder of interest is likely to be more relevant when attempting to understand the contribution of eQTL to disease susceptibility. Evaluating brain-region specific overlap among the brain studies demonstrated highly variable results. Amongst the brain study comparisons there was a higher concordance of eQTL overlap when comparing results from expression levels measured in the same brain region, notably cerebellar tissue.

To consider the effects of neuropathologies on eQTL overlap between brain and blood, individuals with several neurological and psychiatric conditions were included in our comparisons. As most studies pooled data from healthy and diseased individuals, a separate analysis of blood eQTL overlap was not always possible. The exception is Webster et al. [29], who reported eQTL with a significant effect only amongst individuals with Alzheimer’s disease diagnosis. The overlap of blood eQTL with Alzheimer’s disease associated eQTL was lower than eQTL that showed no disease specific effects, suggesting the possibility of increased tissue specificity of genetic regulation of expression levels in individuals with Alzheimer’s disease.

Conclusion

There are several recognizable limitations in this study such as demographic differences, disease status, brain cell heterogeneity, sample size and potential differences in protocols and array platforms. Where possible we have used methods to address these limitations, leading to a picture of eQTL overlap that represents the best-case scenario should studies have used the same array platforms and equal sample sizes. Ideally, measuring expression levels for multiple brain regions and blood in the same individuals would provide the best approach to more fully evaluate the overlap. In summary, although the genetic regulation of expression levels appears to act in a primarily tissue-dependent manner, overlap is still observed although there appears to be no functional differences in the genes with overlapping eQTL. Our results suggest that whenever possible tissue specific to the pathophysiology of the disease being studies should be used for transcription analysis. However, given the availability of blood, and the likely increases in sample size, analysis should not be deemed worthless for informing on brain eQTL associated with neurological and psychiatric conditions.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JEP and NRW designed the study. MM and JEP performed analyses and drafted the manuscript. AKH and AC performed experiments to generate data. All authors read and approved the final manuscript.

Acknowledgements

This research was directly supported by the Australian Research Council (ARC) (FT0991360, DE130100691) and the Australian National Health and Medical Research Council (1050218, 496667). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We thank Jake Gratten for providing critical comment on the manuscript.

References

  1. Schadt EE, Monks SA, Drake TA, Lusis AJ, Che N, Colinayo V, Ruff TG, Milligan SB, Lamb JR, Cavet G, Linsley PS, Mao M, Stoughton RB, Friend SH: Genetics of gene expression surveyed in maize, mouse and man.

    Nature 2003, 422(6929):297-302. PubMed Abstract | Publisher Full Text OpenURL

  2. Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, Spielman RS, Cheung VG: Genetic analysis of genome-wide variation in human gene expression.

    Nature 2004, 430(7001):743-747. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Stranger BE, Forrest MS, Clark AG, Minichiello MJ, Deutsch S, Lyle R, Hunt S, Kahl B, Antonarakis SE, Tavare S, Deloukas P, Dermitzakis ET: Genome-wide associations of gene expression variation in humans.

    PLoS Genet 2005, 1(6):e78. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Nicolae DL, Gamazon E, Zhang W, Duan S, Dolan ME, Cox NJ: Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS.

    PLoS Genet 2010, 6(4):e1000888. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Fu J, Wolfs MG, Deelen P, Westra HJ, Fehrmann RS, Te Meerman GJ, Buurman WA, Rensen SS, Groen HJ, Weersma RK, van den Berg LH, Veldink J, Ophoff RA, Snieder H, van Heel D, Jansen RC, Hofker MH, Wijmenga C, Franke L: Unraveling the regulatory mechanisms underlying tissue-dependent genetic variation of gene expression.

    PLoS Genet 2012, 8(1):e1002431. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Powell JE, Henders AK, McRae AF, Kim J, Hemani G, Martin NG, Dermitzakis ET: Congruence of additive and non-additive effects on gene expression estimated from pedigree and SNP data.

    PLoS Genet 2013, 9:1003502. Publisher Full Text OpenURL

  7. Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, Lai SL, Arepalli S, Dillman A, Rafferty IP, Troncoso J, Johnson R, Zielke HR, Ferrucci L, Longo DL, Cookson MR, Singleton AB: Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain.

    PLoS Genet 2010, 6(5):e1000952. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Hernandez DG, Nalls MA, Moore M, Chong S, Dillman A, Trabzuni D, Gibbs JR, Ryten M, Arepalli S, Weale ME, Zonderman AB, Troncoso J, O'Brien R, Walker R, Smith C, Bandinelli S, Traynor BJ, Hardy J, Singleton AB, Cookson MR: Integration of GWAS SNPs and tissue specific expression profiling reveal discrete eQTLs for human traits in blood and brain.

    Neurobiol Dis 2012, 47(1):20-28. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Dimas AS, Deutsch S, Stranger BE, Montgomery SB, Borel C, Attar-Cohen H, Ingle C, Beazley C, Gutierrez Arcelus M, Sekowska M, Gagnebin M, Nisbett J, Deloukas P, Dermitzakis ET, Antonarakis SE: Common regulatory variation impacts gene expression in a cell type-dependent manner.

    Science 2009, 325(5945):1246-1250. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Heinzen EL, Ge D, Cronin KD, Maia JM, Shianna KV, Gabriel WN, Welsh-Bohmer KA, Hulette CM, Denny TN, Goldstein DB: Tissue-specific genetic control of splicing: implications for the study of complex traits.

    PLoS Biol 2008, 6(12):e1. Publisher Full Text OpenURL

  11. Powell JE, Henders AK, McRae AF, Wright MJ, Martin NG, Dermitzakis ET, Montgomery GW, Visscher PM: Genetic control of gene expression in whole blood and lymphoblastoid cell lines is largely independent.

    Genome Res 2012, 22(3):456-466. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Cai CC, Langfelder P, Fuller TF, Oldham MC, Luo R, van den Berg LH, Ophoff RA, Horvath S: Is human blood a good surrogate for brain tissue in transcriptional studies?

    BMC Genomics 2010, 11:589. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  13. Zeller T, Wild P, Szymczak S, Rotival M, Schillert A, Castagne R, Maouche S, Germain M, Lackner K, Rossmann H, Eleftheriadis M, Sinning CR, Schnabel RB, Lubos E, Mennerich D, Rust W, Perret C, Proust C, Nicaud V, Loscalzo J, Hubner N, Tregouet D, Munzel T, Ziegler A, Tiret L, Blankenberg S, Cambien F: Genetics and beyond–the transcriptome of human monocytes and disease susceptibility.

    PLoS One 2010, 5(5):e10693. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Gregersen PK: Cell type-specific eQTLs in the human immune system.

    Nat Genet 2012, 44(5):478-480. PubMed Abstract | Publisher Full Text OpenURL

  15. Innocenti F, Cooper GM, Stanaway IB, Gamazon ER, Smith JD, Mirkov S, Ramirez J, Liu W, Lin YS, Moloney C, Aldred SF, Trinklein ND, Schuetz E, Nickerson DA, Thummel KE, Rieder MJ, Rettie AE, Ratain MJ, Cox NJ, Brown CD: Identification, replication, and functional fine-mapping of expression quantitative trait loci in primary human liver tissue.

    PLoS Genet 2011, 7(5):e1002078. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Powell JE, Henders AK, McRae AF, Caracella A, Smith S, Wright MJ, Whitfield JB, Dermitzakis ET, Martin NG, Visscher PM, Montgomery GW: The Brisbane Systems Genetics Study: genetical genomics meets complex trait genetics.

    PLoS One 2012, 7(4):e35430. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Schadt EE, Molony C, Chudin E, Hao K, Yang X, Lum PY, Kasarskis A, Zhang B, Wang S, Suver C, Zhu J, Millstein J, Sieberts S, Lamb J, GuhaThakurta D, Derry J, Storey JD, Avila-Campillo I, Kruger MJ, Johnson JM, Rohl CA, van Nas A, Mehrabian M, Drake TA, Lusis AJ, Smith RC, Guengerich FP, Strom SC, Schuetz E, Rushmore TH, et al.: Mapping the genetic architecture of gene expression in human liver.

    PLoS Biol 2008, 6(5):e107. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Wheeler HE, Metter EJ, Tanaka T, Absher D, Higgins J, Zahn JM, Wilhelmy J, Davis RW, Singleton A, Myers RM, Ferrucci L, Kim SK: Sequential use of transcriptional profiling, expression quantitative trait mapping, and gene association implicates MMP20 in human kidney aging.

    PLoS Genet 2009, 5:1000685. Publisher Full Text OpenURL

  19. Khaitovich P, Muetzel B, She X, Lachmann M, Hellmann I, Dietzsch J, Steigele S, Do HH, Weiss G, Enard W, Heissig F, Arendt T, Nieselt-Struwe K, Eichler EE, Paabo S: Regional patterns of gene expression in human and chimpanzee brains.

    Genome Res 2004, 14(8):1462-1473. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Liu C: Brain expression quantitative trait locus mapping informs genetic studies of psychiatric diseases.

    Neurosci Bull 2011, 27(2):123-133. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Borovecki F, Lovrecic L, Zhou J, Jeong H, Then F, Rosas HD, Hersch SM, Hogarth P, Bouzou B, Jensen RV, Krainc D: Genome-wide expression profiling of human blood reveals biomarkers for Huntington’s disease.

    Proc Natl Acad Sci U S A 2005, 102(31):11023-11028. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Martin MV, Rollins B, Sequeira PA, Mesen A, Byerley W, Stein R, Moon EA, Akil H, Jones EG, Watson SJ, Barchas J, DeLisi LE, Myers RM, Schatzberg A, Bunney WE, Vawter MP: Exon expression in lymphoblastoid cell lines from subjects with schizophrenia before and after glucose deprivation.

    BMC Med Genomics 2009, 2:62. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  23. Vawter MP, Atz ME, Rollins BL, Cooper-Casey KM, Shao L, Byerley WF: Genome scans and gene expression microarrays converge to identify gene regulatory loci relevant in schizophrenia.

    Hum Genet 2006, 119(5):558-570. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. de Jong S, van Eijk KR, Zeegers DW, Strengman E, Janson E, Veldink JH, van den Berg LH, Cahn W, Kahn RS, Boks MP, Ophoff RA: Expression QTL analysis of top loci from GWAS meta-analysis highlights additional schizophrenia candidate genes.

    Eur J Hum Genet 2012, 20(9):1004-1008. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Diekstra FP, Saris CG, van Rheenen W, Franke L, Jansen RC, van Es MA, van Vught PW, Blauw HM, Groen EJ, Horvath S, Estrada K, Rivadeneira F, Hofman A, Uitterlinden AG, Robberecht W, Andersen PM, Melki J, Meininger V, Hardiman O, Landers JE, Brown RH Jr, Shatunov A, Shaw CE, Leigh PN, Al-Chalabi A, Ophoff RA, van den Berg LH, Veldink JH: Mapping of gene expression reveals CYP27A1 as a susceptibility gene for sporadic ALS.

    PLoS One 2012, 7(4):e35333. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Colantuoni C, Lipska BK, Ye T, Hyde TM, Tao R, Leek JT, Colantuoni EA, Elkahloun AG, Herman MM, Weinberger DR, Kleinman JE: Temporal dynamics and genetic control of transcription in the human prefrontal cortex.

    Nature 2011, 478(7370):519-523. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Kim S, Cho H, Lee D, Webster MJ: Association between SNPs and gene expression in multiple regions of the human brain.

    Transl Psychiatry 2012, 2(Journal Article):e113. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Liu C, Cheng L, Badner JA, Zhang D, Craig DW, Redman M, Gershon ES: Whole-genome association mapping of gene expression in the human prefrontal cortex.

    Mol Psychiatry 2010, 15(8):779-784. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Webster JA, Gibbs JR, Clarke J, Ray M, Zhang W, Holmans P, Rohrer K, Zhao A, Marlowe L, Kaleem M, McCorquodale DS 3rd, Cuello C, Leung D, Bryden L, Nath P, Zismann VL, Joshipura K, Huentelman MJ, Hu-Lince D, Coon KD, Craig DW, Pearson JV, Heward CB, Reiman EM, Stephan D, Hardy J, Myers AJ: Genetic control of human brain transcript expression in Alzheimer disease.

    Am J Hum Genet 2009, 84(4):445-458. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Zou F, Chai HS, Younkin CS, Allen M, Crook J, Pankratz VS, Carrasquillo MM, Rowley CN, Nair AA, Middha S, Maharjan S, Nguyen T, Ma L, Malphrus KG, Palusak R, Lincoln S, Bisceglio G, Georgescu C, Kouri N, Kolbert CP, Jen J, Haines JL, Mayeux R, Pericak-Vance MA, Farrer LA, Schellenberg GD, Petersen RC, Graff-Radford NR, Dickson DW, Younkin SG, et al.: Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants.

    PLoS Genet 2012, 8(6):e1002707. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Myers AJ, Gibbs JR, Webster JA, Rohrer K, Zhao A, Marlowe L, Kaleem M, Leung D, Bryden L, Nath P, Zismann VL, Joshipura K, Huentelman MJ, Hu-Lince D, Coon KD, Craig DW, Pearson JV, Holmans P, Heward CB, Reiman EM, Stephan D, Hardy J: A survey of genetic human cortical gene expression.

    Nat Genet 2007, 39(12):1494-1499. PubMed Abstract | Publisher Full Text OpenURL

  32. Westra HJ, Peters MJ, Esko T, Yaghootkar H, Schurmann C, Kettunen J, Christiansen MW, Fairfax BP, Schramm K, Powell JE, Zhernakova A, Zhernakova DV, Veldink JH, Van den Berg LH, Karjalainen J, Withoff S, Uitterlinden AG, Hofman A, Rivadeneira F, ‘tHoen PAC, Reinmaa E, Fischer K, Nelis M, Milani L, Melzer D, Ferrucci L, Singleton AB, Hernandez DG, Nalls MA, Homuth G, et al.: Systematic identification of trans eQTLs as putative drivers of known disease associations.

    Nat Genet 2013, 45(10):1238-1243. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O’Donnell CJ, de Bakker PIW: SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap.

    Bioinformatics 2008, 24(24):2938-2939. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Ding J, Gudjonsson JE, Liang L, Stuart PE, Li Y, Chen W, Weichenthal M, Ellinghaus E, Franke A, Cookson W, Nair RP, Elder JT, Abecasis GR: Gene expression in skin and lymphoblastoid cells: refined statistical method reveals extensive overlap in cis-eQTL signals.

    Am J Hum Genet 2010, 87(6):779-789. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. da Huang W, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists.

    Nucleic Acids Res 2007, 35:169-175. OpenURL

  36. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing.

    J R Stat Soc Ser B 1995, 57(1):289-300. OpenURL

  37. Dimas A, Nica A, Montgomery S, Stranger B, Raj T, Buil A, Giger T, Lappalainen T, Gutierrez-Arcelus M, McCarthy M, Dermitzakis E: Sex-biased genetic effects on gene regulation in humans.

    Genome Res 2012, 23:1947-1956. OpenURL

  38. Spielman RS, Bastone LA, Burdick JT, Morley M, Ewens WJ, Cheung VG: Common genetic variants account for differences in gene expression among ethnic groups.

    Nat Genet 2007, 39(2):226-231. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Schliekelman P: Statistical power of expression quantitative trait loci for mapping of complex trait loci in natural populations.

    Genetics 2008, 178(4):2201-2216. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, Hasz R, Walters G, Garcia F, Young N, Foster B, Moser M, Karasik E, Gillard B, Ramsey K, Sullivan S, Bridge J, Magazine H, Syron J, Fleming J, Siminoff L, Traino H, Mosavel M, Barker L, Jewell S, Rohrer D, Maxim D, Filkins D, Harbach P, Cortadillo E, et al.: The Genotype-Tissue Expression (GTEx) project.

    Nat Genet 2013, 45(6):580-585. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Hill K, Mann L, Laws KR, Stephenson CM, Nimmo-Smith I, McKenna PJ: Hypofrontality in schizophrenia: a meta-analysis of functional imaging studies.

    Acta Psychiatr Scand 2004, 110(4):243-256. PubMed Abstract | Publisher Full Text OpenURL

  42. Glatt SJ, Everall IP, Kremen WS, Corbeil J, Sasik R, Khanlou N, Han M, Liew CC, Tsuang MT: Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia.

    Proc Natl Acad Sci U S A 2005, 102(43):15533-15538. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Honea R, Crow TJ, Passingham D, Mackay CE: Regional deficits in brain volume in schizophrenia: a meta-analysis of voxel-based morphometry studies.

    Am J Psychiatry 2005, 162(12):2233-2245. PubMed Abstract | Publisher Full Text OpenURL

  44. Kempton MJ, Salvador Z, Munafo MR, Geddes JR, Simmons A, Frangou S, Williams SC: Structural neuroimaging studies in major depressive disorder. Meta-analysis and comparison with bipolar disorder.

    Arch Gen Psychiatry 2011, 68(7):675-690. PubMed Abstract | Publisher Full Text OpenURL

  45. Meyer-Lindenberg A, Weinberger DR: Intermediate phenotypes and genetic mechanisms of psychiatric disorders.

    Nat Rev Neurosci 2006, 7(10):818-827. PubMed Abstract | Publisher Full Text OpenURL

  46. Savitz J, Drevets WC: Bipolar and major depressive disorder: neuroimaging the developmental-degenerative divide.

    Neurosci Biobehav Rev 2009, 33(5):699-771. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Greicius MD, Flores BH, Menon V, Glover GH, Solvason HB, Kenna H, Reiss AL, Schatzberg AF: Resting-state functional connectivity in major depression: abnormally increased contributions from subgenual cingulate cortex and thalamus.

    Biol Psychiatry 2007, 62(5):429-437. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Fountoulakis KN, Giannakopoulos P, Kovari E, Bouras C: Assessing the role of cingulate cortex in bipolar disorder: neuropathological, structural and functional imaging data.

    Brain Res Rev 2008, 59(1):9-21. PubMed Abstract | Publisher Full Text OpenURL

  49. Lehericy S, Sharman MA, Dos Santos CL, Paquin R, Gallea C: Magnetic resonance imaging of the substantia nigra in Parkinson’s disease.

    Mov Disord 2012, 27(7):822-830. PubMed Abstract | Publisher Full Text OpenURL

  50. Querfurth HW, LaFerla FM: Alzheimer’s disease.

    N Engl J Med 2010, 362(4):329-344. PubMed Abstract | Publisher Full Text OpenURL

  51. Guillozet AL, Weintraub S, Mash DC, Mesulam MM: Neurofibrillary tangles, amyloid, and memory in aging and mild cognitive impairment.

    Arch Neurol 2003, 60(5):729-736. PubMed Abstract | Publisher Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1755-8794/7/31/prepub