- Open Access
Gene expression endophenotypes: a novel approach for gene discovery in Alzheimer's disease
Molecular Neurodegeneration volume 6, Article number: 31 (2011)
Uncovering the underlying genetic component of any disease is key to the understanding of its pathophysiology and may open new avenues for development of therapeutic strategies and biomarkers. In the past several years, there has been an explosion of genome-wide association studies (GWAS) resulting in the discovery of novel candidate genes conferring risk for complex diseases, including neurodegenerative diseases. Despite this success, there still remains a substantial genetic component for many complex traits and conditions that is unexplained by the GWAS findings. Additionally, in many cases, the mechanism of action of the newly discovered disease risk variants is not inherently obvious. Furthermore, a genetic region with multiple genes may be identified via GWAS, making it difficult to discern the true disease risk gene. Several alternative approaches are proposed to overcome these potential shortcomings of GWAS, including the use of quantitative, biologically relevant phenotypes. Gene expression levels represent an important class of endophenotypes. Genetic linkage and association studies that utilize gene expression levels as endophenotypes determined that the expression levels of many genes are under genetic influence. This led to the postulate that there may exist many genetic variants that confer disease risk via modifying gene expression levels. Results from the handful of genetic studies which assess gene expression level endophenotypes in conjunction with disease risk suggest that this combined phenotype approach may both increase the power for gene discovery and lead to an enhanced understanding of their mode of action. This review summarizes the evidence in support of gene expression levels as promising endophenotypes in the discovery and characterization of novel candidate genes for complex diseases, which may also represent a novel approach in the genetic studies of Alzheimer's and other neurodegenerative diseases.
It is well established that the risk for many neurodegenerative diseases such as Alzheimer's disease (AD)[1, 2], Parkinson's disease (PD)[2, 3], frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS)[4, 5] is under substantial genetic control. Identification of deterministic mutations with a Mendelian pattern of inheritance in familial forms of these diseases has led to significant progress in our understanding of their underlying pathophysiology. Nonetheless, such monogenic forms constitute only a fraction of these conditions (e.g. < 1% of all AD, 5-10% of all PD); and risk for most neurodegenerative diseases, like other common diseases/traits, is thought to arise from multiple genetic variants and their interaction with each other, as well as environmental factors. Whereas investigation of large families with Mendelian inheritance patterns, using linkage followed by sequence analyses have generally proven useful in gene discoveries for monogenic forms of neurodegenerative diseases[1, 5], association studies in unrelated case-control series have emerged as a more viable strategy in the common, multigenic and complex forms of these conditions[1–3]. The relative lack of success from candidate gene association studies coupled with technologic advances led to the emergence of genome-wide association studies (GWAS) as a promising approach in gene discovery for common, complex diseases/traits with heterogeneous and multigenic underpinnings.
The past several years have witnessed an explosion of GWAS that survey hundreds of thousands of genetic variants across the whole genome for association with disease risk in a hypothesis-independent manner. The recent widespread use of this approach became possible with the generation of a linkage disequilibrium (LD) map of common single nucleotide polymorphisms (SNPs) across the whole genome as part of the International HapMap Project[6, 7]. This effort combined with the technological advances in high-throughput genotyping allowed for the production of chips that contain up to 1 million SNPs which provide information about most (but not all) of the common genetic variation (usually defined as having a minor allele frequency > 0.05) in the genome by tagging (marking) the un-genotyped SNPs through LD. GWAS emerged based on the underlying "common disease-common variant" hypothesis which posits that the risk for many common diseases such as coronary artery disease, diabetes mellitus and AD is conferred by multiple common variants . This approach proved to be successful in the discovery of candidate risk variants or regions for many common diseases, including neurodegenerative diseases[1–3].
It became evident from these studies that although GWAS have identified numerous genetic loci for common diseases, they fell short of accounting for all of the genetic component of these conditions[8, 9]. The reasons for this have been discussed in recent reviews[8, 9] and include a) modest power of GWAS that utilize disease phenotype given relatively small effect sizes of risk variants, genetic heterogeneity of the disease (different disease risk variants at play in different subjects), and heterogeneity of the populations (different extent of LD); b) presence of rare disease risk variants or structural variants (such as insertions and deletions) that are not captured by the GWAS SNP platforms, c) presence of gene-gene and/or gene-environment interactions that are yet unexplored.
The SNPs identified through GWAS may reside near genes with no known or disease-relevant function. Some disease-associating variants may reside in LD regions harboring multiple genes, making it impossible to discern the true disease-risk gene. These represent relative weaknesses of GWAS resulting in part from the hypothesis-independent nature of this approach. On the contrary, this may also be perceived as a relative strength, since it may allow for the discovery of unexpected genes that would have been missed with a hypothesis-based approach. Nevertheless, the plethora of loci identified through GWAS need both functional confirmation and characterization, since their roles in disease risk are usually not inherently obvious.
Multiple different approaches are proposed to overcome these shortcomings of GWAS that are beyond the scope of this review. One of the proposed approaches involves the use of biologically relevant, quantitative phenotypes (endophenotypes) for discovery and characterization of disease risk genes[8, 9]. This review focuses on gene expression levels as a potentially powerful group of endophenotypes and discusses the endophenotype concept, the evidence in support of a substantial genetic component for human gene expression, the GWAS examples that combine disease phenotype and expression endophenotypes, the use of gene expression endophenotype in the Alzheimer's and other neurodegenerative disease literature to-date and future directions.
The Endophenotype Concept
The term "endophenotype" was first introduced in 1966 in the context of Drosophila genetics to contrast the phenotypes that are "obvious and external" (i.e. exophenotypes) with those that are "microscopic and internal" (i.e. endophenotypes). The term was applied in psychiatric genetics in the 1970s, specifically in schizophrenia genetics to describe molecular outcomes of gene activity, which lead to disease. It was not until the 21st century that this term reached more widespread usage, mainly in psychiatric genetics initially, but subsequently also in other complex conditions such as neurodegenerative[12, 13], cardiovascular and atopic  diseases. That the endophenotype approach was first advocated in psychiatric disorders partly stems from the need to have objectively quantifiable phenotypes associated with disease to avoid the imprecise and therefore heterogeneous nature of psychiatric diagnostic criteria, which is thought to contribute to the failure of genetic studies of psychiatric disorders. In a hypothetical construct that defines genotype-disease relationship, endophenotypes were proposed as quantitative phenotypes that are intermediate between genes and the disease outcome, that are directly influenced by a smaller number of genes than the disease phenotype, and that represent one of many facets of a disease. Thus, another premise for this approach is the assumption that the underlying genetic architecture of the endophenotype will be simpler than that of the disease phenotype, though this may clearly not be the case. Nevertheless, their closer proximity to gene action could lead to greater "genetic signal-to-noise" ratios[17, 18], which may translate into larger effect sizes for the genetic variants under investigation, and thereby increased power to detect genetic loci of interest. Furthermore, given that endophenotypes are measurable traits regardless of disease state, their use can allow the inclusion of unaffected as well as affected subjects in genetic studies, which can also enhance power, especially in family-based studies[18, 19].
Genes identified via the endophenotype approach may be more amenable to the study and discovery of genetic pathways underlying disease. This is partly due to the quantifiable nature of the endophenotype, which makes possible objective testing of the downstream effects of genetic perturbations, including development of animal models. A potential caveat of the endophenotype approach is the fact that it represents only one of many aspects of the disease pathology. This may be a limitation for animal models that study the endophenotype rather than disease symptomatology as the outcome. That said, given the improbability of recapitulating the complete spectrum of a complex disease in an animal model and the potential difficulty in drawing parallels between animal and human symptomatology, a focused study of a "good" endophenotype could enhance our understanding of complex disease pathophysiology.
What constitutes a "good" endophenotype? A variety of criteria have previously been proposed to define "valid and useful" endophenotypes in a review, where a core set of "necessary and sufficient" criteria was identified. Although a detailed discussion of these criteria is beyond the scope of this review, we will briefly mention those that are commonly proposed. Endophenotypes should be measured reliably and reproducibly upon multiple measurements of the same subject/sample. Given the proposition that endophenotypes are traits that can be used to map genetic variants underlying disease risk, heritability or having a measurable genetic component is a sine qua non of a useful endophenotype. Neuroimaging and cognitive[17, 21] endophenotypes are examples of heritable traits that are proposed for genetic studies of neuropsychiatric diseases. The underlying assumption is that genes influence the endophenotype in a way that is detectable prior to the clinical onset of disease. These changes in the endophenotype in turn leads to increased disease risk. Consequently, endophenotypes should be associated with the disease in the general population and also co-segregate with it within families. They should manifest changes that are detectable in the clinically-unaffected but at-risk subjects, such as family members of patients. We determined that plasma Aβ levels, which show variation in the general population, show significant elevations in cognitively-normal first-degree relatives of patients with late-onset AD and are highly heritable traits. Thus, we used plasma Aβ as an excellent endophenotype in a linkage study of late-onset AD pedigrees and mapped an AD risk locus on chromosome 10. Endophenotypes need not, however, be disease-specific. For example, both neuroimaging and cognitive endophenotypes are quantitative traits that are variable and detectable in the general population. These endophenotypes may show changes that are associated with more than one psychiatric or neurologic disorder. Ideally, endophenotypes should be state-independent, and not change with the disease state or environmental factors. For example, although potentially useful endophenotypes for AD and possibly also other neurodegenerative diseases, both neuroimaging and cognitive endophenotypes are influenced by the disease state, as well as other variables such as age and gender. It is therefore important to recognize and either statistically or technically control for such variables in genetic analysis of endophenotypes.
The evidence supporting a strong genetic component that influences human gene expression is discussed in the next section.
Genetics of Human Gene Expression
In this section, 15 key studies using expression levels as an endophenotype to identify genetic loci or variants that influence human gene expression are reviewed (Tables 1 and 2). This approach is known as expression quantitative trait locus (eQTL) mapping via genetic linkage in families or association studies in unrelated populations. This section describes these studies highlighting their results which provide support for the existence of a substantial genetic component for gene expression in humans and how this component could be utilized to study the genetics of complex diseases and traits.
Genetic linkage analysis of whole transcriptome expression levels in yeast determined that 1,528 of the 6,215 genes tested had natural variation in their expression levels and 570 of these gene expression phenotypes showed linkage to ≥ 1 loci in the yeast genome. This study demonstrated that genetic factors account for a substantial proportion of variation in gene expression levels, even in a single cell organism. The first study in humans to evaluate genetics of gene expression (also known as the genetical genomics approach) utilized lymphoblastoid cell lines (LCL) from the CEPH (Centre d'Etude du polymorphisme humain) repository of Caucasian, multi-generational families. Yan et al. screened 96 subjects and identified 17-37 individuals who were heterozygous for a SNP in 13 target genes. Measuring the mRNA of these 13 genes in the same cellular sample, they assessed relative expression of the two alleles of the SNPs and identified evidence of allele-specific expression for 6 of 13 genes in a subset of their samples. Examining the gene expression levels as a phenotype, they identified evidence of Mendelian inheritance for two expression phenotypes that co-segregated with nearby genetic markers in two families. These results provided the foundation to utilize the "genetical genomics" approach in a more high-throughput and systematic fashion.
Schadt et al. provided a comparative analysis of gene expression genetics in mouse, Zea mays and human, in a pioneering study which showed the generalizability of this approach to different organisms. In this study, an eQTL linkage analysis was performed using 23,574 transcript expression levels measured in livers of 111 mice from the F2 generation of two standard inbred strains, C57BL/6J and DBA/2J and >100 microsatellite markers. Differential gene expression was detected for >33% of the genes in the parental and ≥10% genes in the F2 strain. 9-16% genes have eQTLs with LOD scores > 4.3. To identify transcripts that influence the complex obesity trait in these mice, the authors compared the gene expression profiles of mice at the lower vs. upper 25th percentile of an obesity trait and identified three distinct expression patterns for 280 genes. Importantly, these 280 genes were enriched for eQTLs and linkage analysis of the obesity trait performed on the subgroups with distinct expression patterns improved the linkage signal. These results established a paradigm for the combined use of gene expression traits and another complex trait of clinical relevance for improved mapping of genetic factors that influence the complex clinical trait by affecting gene expression levels.
Analysis of corn (Z. mays) ear leaf tissue using the eQTL approach in the same study identified that 77% of genes had differential expression in this organism, 26% had ≥ 1 eQTL with LOD score > 3.0 and there appeared to be genetic interactions between some of the eQTLs. Schadt et al. also studied a small number of human LCLs of 56 subjects from four CEPH families by variance components analysis and identified differential expression for 11% of the genes assayed of which about a third had detectable heritability. Overall, these findings demonstrate the complexity of gene expression genetics ranging from a "simpler" plant organism to humans, but also the applicability of genetical genomics approach in tracking eQTLs in different organisms.
One of the first studies which established the natural variation of gene expression levels in humans, evaluated LCLs of 35 unrelated subjects from CEPH families vs. 1 reference pool of 10 subjects utilizing 3-4 replicate measurements per person. They determined that for many genes (n = 813) between-person gene expression variations were higher than within-person variations, which are due to technical variability. Evaluation of 5 genes, which revealed highest variance among unrelateds, then sibs then monozygotic twins, provided proof of principle for a genetic component underlying at least some of the variability in human gene expression.
Seven of the studies[31–37] following these initial reports utilized human LCL for eQTL mapping. Three of these studies[31–33] assessed cell lines from the CEPH repository, similar to the prior reports[28–30], three utilized LCL from the HapMap consortium[6, 34–36] and one studied samples from families with an asthma proband. Morley et al. measured levels of 8,500 transcripts in 94 unrelated subjects from CEPH families and identified 3,554 (42%) transcripts with greater between-subject (biological) variation than within-subject (technical). Genome-wide linkage analysis of these 3,554 expression phenotypes in 14 CEPH families (n~140) detected evidence of significant linkage for 984 (28%) transcripts with genome-wide p < 0.05 and for 142 (~4%) with genome-wide p < 0.001, which far exceeds the numbers expected by chance alone. When these eQTLs were distinguished as cis- or trans-regulators, defined in this study as linkage regions within or outside 5 Megabases (5 Mb) of the target gene, respectively, most of the top 142 eQTLs were found to be trans-regulators. Some expression traits had multiple significant eQTLs. There were linkage regions with multiple expressions linking to it, termed as expression "hotspots" or "master regulators". Importantly, genes that map to one hotspot had expression levels with higher than expected correlations and some of them had close physical locations.
Monks et al., also assessed LCLs from CEPH families (n = 167), though they measured a larger number of transcripts than Morley et al., with 23,499 expression phenotypes that were evaluated by variance components analysis for their heritability, as well as eQTL analysis. They determined that 10% of the genes were differentially expressed in ≥50% of children, and a third of these were heritable at a false discovery rate (FDR) of 0.05 and median heritability of 0.34. The heritable transcripts were enriched for immunity pathways. Twenty-two genes had significant eQTLs at genome-wide level, eight of which were within 5 Mb of the gene (i.e. cis-regulating). In contrast to Morley et al., Monks et al. did not find an enrichment for transcription "hotspots" over what would be expected by chance alone, based on simulation studies. The authors studied a subset of 574 transcripts with heritability for pairwise genetic correlations (GC) and Pearson's correlations (PC), which does not account for a genetic component, followed by pathway analysis using the Gene Ontology (GO; http://www.geneontology.org/) and Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.genome.jp/kegg/) databases. This analysis revealed that genes clustered by both GC and PC had similar pathway coherence for GO, but that GC gene clusters had better pathway coherence for KEGG pathways. This suggests that analysis of pairwise genetic correlations between gene expression phenotypes may identify novel biological pathways that may not be possible by approaches which neglect the genetic component of correlations.
In a follow up to their prior linkage study, Cheung et al. assessed 57 unrelated CEPH subjects in an association study for 374 expression phenotypes with prior evidence of linkage at p < 0.02 and cis-SNPs near linkage peaks (± 50 kb of target gene). Additionally, they evaluated the top 27 cis-linkage expression phenotypes (p < 3.7 × 10-5), in a GWAS using >770,000 SNPs. They determined that 65 of 374 expression levels have ≥1 SNP that associates at nominal p < 0.001, 12 with p < 10-10 and 133 with p < 0.01. Same proportions of associations were found for the 5', 3' and genic regions. Fourteen out of the top 27 cis-linkage regions showed associations significant at the genome-wide level after Bonferroni corrections for the number of SNPs tested. Twelve of those top 14 associations were cis only, one was cis+trans and one was trans only. The authors also performed functional analysis for one of the top genes and confirmed the presence of a functional variant that influenced gene expression by modifying strength of RNA polymerase II binding. This study demonstrated that strong linkage also predicts strong association for expression levels and expression GWAS (eGWAS) may be a feasible and powerful approach to identify genetic determinants of expression phenotypes.
In the first of three studies[34–36] performed using HapMap LCL, Stranger et al. analyzed 60 samples for levels of 1,433 transcripts (630 genes), then performed eGWAS for 688 transcripts (374 genes) with highly variable expression signals above the background. They compared and identified good concordance between three methods for multiple corrections: Bonferroni, FDR and permutations. There were 10-40 genes which had cis-SNPs (defined as 1 Mb from genomic midpoint of gene) significant at genome-wide level by ≥1 statistical method, whereas trans signals were detected only for 3 genes.
In a subsequent larger study, Stranger et al., investigated LCL from four HapMap populations: 30 Caucasian trios (CEU), 30 Yoruba trios (YRI), 45 unrelated Chinese (CHB) and 45 unrelated Japanese (JPT). Analysis of 14,456 transcripts (13,643 genes) revealed a genetic component with heritabilities > 0.2 in 10% (4,829) and 13% (6,482) of the transcript probes analyzed in the CEU and YRI trios, respectively, with 958 overlapping genes. There were 154 CEU and 217 YRI genes with substantial heritability >0.5, with overlap of 9 genes. There were 831 genes (6% genes tested) with a cis-association significant in at least 1 population at p < 0.001 after 10,000 permutations; 310 in at least 2 and 62 in all four populations. There were many less trans-associations, with 108 significant genes in ≥1 population, 16 genes in ≥2 and 5 in all 4 populations. These findings demonstrate the presence of both overlap and diversity in significant eSNPs across different ethnic groups. More than 50% of the significant eQTLs had heritability estimates > 0.2, though pooling populations to increase sample size captured some additional genes with smaller effect sizes. Most trans-SNPs were also found to have cis-effects and most cis-associations were in genic and immediate intergenic regions, suggesting that cis-variants may be more abundant and/or may have stronger effects than trans-variants. This study also evaluated the influence of technical variability on gene expression measurements by testing 60 cell lines on two different arrays, that led to results with high correlations, suggesting that transcript measurements are stable across different experiments, measurement times and platforms. They also identified that CEU population, which had the most aged cell lines, also had the most divergent expression profile from other populations, drawing attention to the potential technical concerns with LCL.
To investigate the contribution of a different type of genetic variation, namely copy number variants (CNVs) to gene expression levels, Stranger et al., performed association of 14,925 transcripts (14,072 genes) with CNVs in the same four HapMap populations and determined that there are significant CNV associations that replicate across ethnic groups as well as those that are unique to one ethnic group. Mapping of CNVs vis à vis the genes suggested that CNVs exert their effects by disrupting both regulatory regions as well as the genic regions and altering gene dosage. Notably, most CNV associations cannot be detected by SNPs (87%), indicating the importance of the survey of structural variants in addition to SNPs in eGWAS.
In a study of 400 LCL from affected and unaffected children of families with an asthma proband, Dixon et al., studied 54,675 transcripts (20,599 genes). First, there were no significant differences in the expression profiles between asthmatics and non-asthmatics (unchallenged cells). Twenty-eight percent of the transcripts (15,084) corresponding to 6,660 genes have heritabilities > 0.3, which is in good agreement with prior reports[29, 32]. 1,989 transcripts had significant associations at genome-wide level, where many of the strongest associations were in cis. This study did not identify master regulators with strong effects. Combined investigation of these eGWAS results with those of published GWAS on some diseases, identified potential candidate genes for asthma, Crohn's disease, diabetes and thalassemia implying that eGWAS is a useful approach to detect disease SNPs with a functional role.
In one of the largest eQTL studies to date, Goring et al. analyzed 1240 lymphocyte (not LCL) samples collected from multigenerational Mexican-American families from the San Antonio Family Heart Study (SAFHS) for 19,648 transcripts (18,519 genes). Nearly 85% of the transcripts (16,678) were significantly heritable with a median heritability estimate of 22.5%. This estimate is similar to those from some LCL studies[29, 32, 37]. There were 1,345 significant cis-regulated transcripts with median effect size of 24.6%. There were many less trans-regulators, again reinforcing the idea that trans effects likely have smaller effect sizes than those in cis. This study was able to replicate the cis-findings of Morley et al., but not those in trans, which suggests higher stability of cis-regulatory effects across studies, populations and cell types, contrary to trans-effects, which may either be more cell/tissue-specific, may have smaller effect sizes or may be false positives. There was also no evidence of "master regulators" in this study. The authors demonstrated the utility of the eQTL approach in complex trait mapping, by identifying eSNPs which associate with HDL-cholesterol levels and expression levels of VNN1. They ultimately identified variants in this gene with functional consequences on transcription binding.
The first eQTL study to systematically compare two different tissues and pursue combined genetic analysis of expression levels and complex traits -in this case obesity-related phenotypes- assessed 20,877 transcripts in the blood (IFB cohort; n = 1,002) and adipose tissue (IFA cohort; n = 673) of Icelandic families, genotyped for 1,732 microsatellite markers for linkage analysis and also for 317,503 SNPs for association analysis in a subset of 150 unrelated subjects. In this study, adipose tissue expression levels (63-72%) correlated better with obesity traits than do blood expression levels (3-9%), and this effect was more pronounced when the analysis was confined to the subset of subjects that overlapped in the blood and adipose tissue cohorts. Fifty-five percent of blood and 75% of adipose tissue transcripts were significantly heritable, with average heritability estimates of 30%, similar to the estimates from prior studies[29, 32, 37, 38]. Variance components linkage analysis revealed 2,529 (12%) significant cis-eQTLs in blood, and 1,489 (7%) in adipose tissue, where cis-eQTL region corresponded to the microsatellite located nearest the expression probe of interest. Greater than 50% of the significant adipose tissue cis-eQTLs were also significant in blood. Expression traits with higher heritability in both tissues had greater reproducibility of eQTL signals. Thus, although there is significant overlap in genetic factors underlying gene expression in two different tissue types, expression levels from the clinically-relevant tissue appears to correlate better with clinical-phenotypes. There were many less significant trans-eQTLs and no evidence of "master regulators" above what would be expected by chance alone. There were 2,714 (12.9%) significant cis-SNPs in blood and 3,364 (16%) in adipose tissue, where cis-SNPs were defined as those residing within a 2 Mb window centered at the location of the probe corresponding to the transcript of interest.
This study, by Emilsson et al., also characterized the transcriptional network by evaluating all pair-wise correlations in the most differentially expressed genes and generated a "connectivity map", which defined "connectivity" of a given gene as "the sum of its connection strengths with all other genes in the network". This approach identified a group of highly correlated gene-expressions (termed "network module") in human adipose tissue, which significantly overlapped with a "network module" in mouse adipose tissue. This module was found to be significantly enriched for genes in macrophage activation pathways. These genes also had cis-eSNPs which significantly influenced expression levels as well as obesity-related traits. These results collectively demonstrate that expression correlation networks combined with cis-eSNPs could potentially identify genes/pathways underlying complex clinical phenotypes.
Using similar approaches with Emilsson et al., but human liver tissue from 427 unrelated subjects, Schadt et al., performed an eGWAS using levels of 39,280 transcripts (34,266 genes) and 782,476 SNPs. At Bonferroni adjusted p < 0.05, there were 1,350 expression traits (1,273 genes, 3.7%) and at FDR <10%, 3,210 traits (3,043 genes; 8.8%) which had ≥1 significant cis-eSNP, explaining 2-90% of variation in expression levels. In contrast, significant trans-eQTLs were far fewer with only 242 traits (236 genes) or 491 traits (474 genes) showing significance after Bonferroni or FDR corrections, respectively. Comparison of their former blood and adipose tissue eQTLs to the liver eQTLs determined 30% overlap of cis-eQTLs from each tissue with significant liver cis-eQTLs. By combining their liver eGWAS with publicly available disease GWAS data and transcriptional network approach, Schadt et al. identified SORT1 and CELSR2 as candidate genes for coronary artery disease and LDL cholesterol levels, and RPS26 for Type 1 diabetes. These findings have implications for mapping those loci that influence disease risk by affecting gene expression levels. They demonstrate how the combined eQTL and disease mapping approaches can help distinguish the actual disease susceptibility gene in regions of LD, where SNPs in more than one gene may associate with disease risk, but only the disease-related transcript levels will associate with the eSNPs.
There are two published eGWAS utilizing transcript levels from human brain tissue[41, 42]. The first study by Myers et al. assessed neuropathologically "normal" cerebral cortical tissue from 193 subjects for 14,078 transcripts and 366,140 SNPs. They determined that 58% of the transcriptome had expression in ≥5% of control brains. Of these, 21% had significant associations with a cis- or trans-eSNP. Contrary to prior studies[34, 35, 37–40] performed in human LCL, lymphocyte, adipose and liver tissue, human brain eGWAS identified less significant cis-eSNPs (433 eSNPs for 99 transcripts) than trans-eSNPs (16,701 eSNPs for 2,876 transcripts). There was, nevertheless, enrichment of significant cis- vs. trans-associations over chance expectations at distances close to the transcripts and maximizing at distances <70 kb. Myers et al. determined that MAPT cis-eSNP alleles that are on the major H1 haplotype are associated with higher MAPT transcript levels, in agreement with their prior study. Compared to Cheung et al., and Stranger et al., there were only two common results (1 cis and 1 trans-association), in contrast to 30% overlap of cis-eQTLs between liver and blood or adipose tissue transcriptomes. These differences could arise from differences in platform and/or may suggest a more distinct genetic control for the brain transcriptome. Technical challenges arising from measuring levels of potentially degraded transcripts in post-mortem brain tissue, in comparison to other tissues that may have more well-preserved RNA, may also underlie some of these discrepancies.
In the follow-up brain eGWAS, Webster et al., assessed 176 brain samples with AD neuropathology, and performed a joint evaluation with 188 control brains from their prior study. The analyses were restricted to 8,650 transcripts out of 24,357 (35.5%), which were detected in >90% of the cases and controls. Similar to control brains, 58% of the transcriptome was detectable in ≥5% of AD brains. Assessment of technical and biological covariates revealed that hybridization date and APOE had the strongest influence and post-mortem interval had the least influence on brain expression levels. There were 1,829 significant cis-eSNPs (within gene ± 1 Mb of 3' or 5' end) in their combined sample and 656 significant trans-eSNPs. Twenty-seven percent of all eQTLs were found to have significant interaction term with diagnosis. Of the cis-eSNPs with significant diagnosis interaction, 37% were also found in just the control brains. Thus, a subset of cis-eSNPs that influence transcript levels differentially for AD subjects, also influence gene expression in control brains, suggesting that a portion of disease-relevant eSNPs may be captured in disease-free tissue and also that presence of other factors besides the disease-relevant eSNPs are likely necessary to predispose to AD. The authors did network analysis on transcripts with a significant eQTL and those that were differentially expressed between ADs and controls, but did not have a significant eQTL. This led to the identification of some transcript clusters that were significantly enriched for gene ontology groups; and "hub" genes with expression levels that correlate with many other transcripts. The authors highlighted their findings of significant cis-eSNPs that influence expression of GSTO2, but not GSTO1, which were previously implicated in AD risk[46, 47], thereby implicating GSTO2 as the likely AD risk gene in this region.
Collective evidence from the eQTL studies discussed in this section lead to the following conclusions:
Transcriptome expression levels can be repliably and replicably detected in human cell lines and multiple types of tissues.
Genetic factors account for a substantial proportion (3-85% depending on study size, tissue source, expression platform) of the variation in human gene expression, with median heritability estimates of 20-30% in most studies. This is similar to the genetics of gene expression for other organisms.
Expression QTLs can be mapped by linkage or eGWAS approaches.
Most eQTLs appear to be cis-regulating, suggesting that trans-eQTLs may have smaller effect sizes.
Both SNPs and structural variations, namely CNVs, appear to underlie eQTL effects.
Many eQTLs are common across different ethnic groups, though there appear to be eQTLs that are unique to one ethnic group.
Many eQTLs are common across different tissue types, though there appear to be eQTLs that are unique to one tissue type.
"Gene expression endophenotype" approach is powerful and can detect hundreds to thousands of significant eQTLs with hundreds of subjects, unlike "complex disease mapping" approaches. Increasing sample size has the most impact on power for eQTL studies (vs. increasing SNP markers).
"Gene expression endophenotype" approach can be utilized in conjunction with mapping for complex diseases or other disease-related phenotypes to identify or confirm novel genes with functional eSNPs that confer disease risk. "Network" analysis of human eQTLs alone or in conjunction with those from other organisms can identify novel biological pathways that may be "disease-relevant".
"Disease-relevant" eQTLs are more likely to be captured in "disease-relevant" tissue (e.g. obesity-related eQTLs in adipose tissue), although many such eQTLs may also be identified in other tissues. There appears to be substantial overlap in eQTLs detected from subjects with and without disease, although detection of some eQTLs may require assessment of "disease-relevant" tissue from subjects with disease.
Gene Expression Endophenotypes in Neurodegenerative Diseases-Current Status
Gene expression endophenotypes have thus far been utilized in two types of studies in neurodegenerative diseases: "Candidate gene" studies, where influence of variants on expression levels of one or a few candidate genes are assessed; and "Expression profiling" studies, which compare levels of transcriptome expression in tissues from patients vs. controls in an effort to identify novel disease genes and pathways. This section provides examples of both types of studies to demonstrate the most common use of gene expression endophenotypes in neurodegenerative diseases to date. While not a comprehensive review on these studies, it will provide knowledge on these current approaches, including their uses and limitations.
Gene Expression Endophenotypes in Neurodegenerative Diseases-"Candidate gene" approach
The initial studies utilizing the gene expression endophenotype in neurodegenerative diseases, have been conducted for known disease risk genes. A prime example of this is Apolipoprotein E gene (APOE), which is a well-established risk factor for AD, where the coding ApoEε4 polymorphism is associated with increased risk, but is neither sufficient nor necessary for its development (reviewed in). Polymorphisms in the promoter region of APOE, which influence its expression, have received some attention as potential modifiers of AD risk that may be independent of ApoEε4 (reviewed in ). Although gene expression levels have not been directly used as endophenotypes in most of these studies, these polymorphisms were investigated because of their predicted influence on ApoE expression based on their location and in-vitro functional assays. One of the most well-studied of these polymorphisms is -491A/T, which was initially found to confer risk for AD in its -491AA homozygote form, even in subjects who lacked the risky ApoEε4 allele. Functional transcriptional studies identified a stronger promoter activity for APOE -491A vs. -491T containing constructs. ApoEε4 with the -491AA genotype had greater AD risk compared to those with one or no copies of the -491A allele, suggesting that both the isoform and level of expression of ApoE may be important in conferring disease risk and that ApoEε4 subjects with higher expression of this protein may be at highest risk. While this notion is in agreement with findings from animal studies, the influence of ApoE promoter polymorphisms on brain ApoE levels have not been conclusive.
Another gene implicated in multiple neurodegenerative diseases, which was investigated with the expression endophenotype approach, is microtubule associated protein tau (MAPT). Missense and exon 10 splicing mutations in MAPT lead to frontotemporal dementia with parkinsonism linked to chromosome 17 (FTDP-17), whereas a region of LD within MAPT, known as the H1 haplotype is associated with increased risk of taupathies, namely, corticobasal degeneration (CBD) and progressive supranuclear palsy (PSP) compared with the (reviewed in ). Given the increase in MAPT exon 10-containing transcripts (also known as 4 repeat or 4R tau) in affected brain regions in PSP and CBD, Caffrey et al.[55, 56] performed allele-specific gene expression studies in both human neuronal cell lines and brain tissue of MAPT H1/H2 heterozygous subjects. They determined that the risky MAPT H1 haplotype is associated with significantly higher MAPT exon 10-containing transcript expression without significant total MAPT levels compared to MAPT H2 haplotype. Myers et al. identified MAPT H1c as the sub-haplotype associated with an increase in both total and 4R tau levels in the human brain, with higher in-vitro transcriptional activity and with increased AD risk.
While the above examples demonstrate how the gene expression endophenotypes can be exploited to uncover the underlying biology of well-established genes in neurodegenerative diseases, our work utilizing the cerebellar expression levels of 12 late-onset AD (LOAD) candidate genes in 200 AD brains illustrate the use of this approach in identification of novel functional disease risk polymorphisms. In this study, Zou et al. investigated association of 619 cis-SNPs with cerebellar expression levels of 12 LOAD candidate genes and identified three significant cis-SNPs in insulin-degrading enzyme (IDE). The top cis-SNP (rs7910977) reached genome-wide significance where the minor allele led to ~twofold increase in cerebellar IDE mRNA levels, reduced AD risk, and reduced plasma Aβ levels, which is biologically congruent. This IDE cis-SNP was in complete LD with an IDE SNP in a conserved region (rs6583817) and increased reporter gene expression in an in-vitro assay, providing additional evidence for a functional effect of this polymorphism on gene expression.
Gene Expression Endophenotypes in Neurodegenerative Diseases-"Expression Profiling" approach
Large-scale comparisons of gene expression levels in subjects with disease vs. controls, known as "Expression Profiling" is one of the most commonly used approaches implementing the gene expression endophenotype. While potentially illuminating, this approach is vulnerable to technical confounders that may influence gene expression differentially in subjects with disease vs. controls, including tissue and RNA integrity. To minimize these confounds, disease vs. control tissues need to be carefully matched and/or these variables need to be accounted for in the statistical analyses. Additional technical confounders to consider and control are microarray platform and quality and batch effects for experiments conducted on different dates. There are important statistical considerations including variability in expression levels, which may lead to false positive findings especially for low expressing genes with high variability, as well as false negative results for those with small, but reproducible changes. Perhaps the most important biological caveat is that, expression profiling design, especially if conducted in the disease-relevant tissue, does not distinguish expression changes that are a result of the disease process from those that are underlying causes of it. In this aspect, it is inferior to the eQTL approach, which can be designed to uncover the genetic factors underlying expression changes and disease risk. Despite these pitfalls, the potential utility of the expression profiling approach will be discussed in this section, highlighting the results from several expression profiling studies in neurodegenerative diseases.
In a hippocampal gene expression profiling study of 9 control and 22 AD subjects of various severity, determined by the bedside cognitive test, Mini Mental State Examination (MMSE), Blalock et al. correlated expression levels with both cognitive (MMSE) and neuropathology (neurofibrillary tangle) scores, in all subjects as well as the subset of nine controls or mild ADs with MMSE scores of 20-26 (collectively termed as the "Incipient ADs"). They identified upregulation of genes that pertain to transcription factor and tumor suppressor pathways among others. The small sample size and the inability to discriminate expression changes "due to" vs. "underlying" disease are the main concerns with this otherwise novel approach, which utilized three endophenotypes.
Bossers et al. utilized a similar approach by generating expression profiles of 49 prefrontal cortex samples from subjects with different severity of underlying AD neuropathology detected by Braak staging for neurofibrillary tangles. There were 1,071 transcripts (922) genes which showed significant changes in their expression levels with Braak stage. Some of these findings were validated with quantitative PCR (qPCR). There were clusters of genes which appeared to show concerted changes with advancing Braak stage, such as increasing early and then decreasing (UPDOWN clusters) or vice versa (DOWNUP clusters), where the biggest changes coincide with the appearance of amyloid plaques at Braak stage III. Functional annotation and pathway analysis of these clusters, revealed an enrichment for synaptic genes in the UPDOWN cluster and those involved in proliferation, differentiation and inflammation in the DOWNUP clusters. Importantly, the synaptic gene expression changes correlated with Aβ levels, which led to the conclusion that synaptic activity and Aβ production may be part of a feed-back loop that ultimately leads to AD. Given the appearance of changes even before significant neuropathology or clinical decline, the authors suggest that these findings are not a result of neuronal loss due to disease process.
While the previous examples are aimed at understanding the role of gene expression changes in the pathophysiology of a single neurodegenerative disease[60, 61], Bronner et al. performed expression profiling in the medial temporal cortex of 5 patients each from four distinct disease categories, namely, PSP, FTD, AD and PiD (Pick's disease) as well as a control group. Comparison of gene expression profiles between each disease group against controls identified a set of 166 transcript probes, which could discriminate PSP, FTD/PiD from controls and each other. The FTD and PiD groups have similar gene expression profiles. AD had the most similar profile to control group, in this small study. Although given its very limited sample size, this study should be considered as a pilot, it also demonstrates another potential application of gene expression profiling as a means to discriminate between neurodegenerative disorders at a molecular level.
Combining genetic association studies with gene expression profiling represents another paradigm utilizing this approach[64, 65]. Taguchi et al. identified 35 genes that were significantly up- or down-regulated in the hippocampus of AD vs. control subjects, which they tested for genetic association with AD risk in 376 AD patients vs. 376 controls. This study identified nine nominally significant AD risk associations, a higher significance rate than most such studies, with the POU2F1 association also reaching study-wide significance. Similarly, Chapuis et al. performed case-control genetic association analysis for 82 genes that were found to be differentially expressed between 12 control and 9 AD brains, and determined nominally significant associations with AD risk for 17 genes. Of these, the association for IL-33 achieved study-wide significance, which was replicated in three additional series and was found to interact with ApoEε4. The rare allelic variants in IL-33 were associated with decreased risk of AD risk and reduced levels of cerebral amyloid angiopathy. Furthermore, overexpression of IL-33 was associated with decreased Aβ40 secretion in-vitro. Though replication studies are necessary, these studies illustrates the utility of gene expression profiling studies in identifying candidate disease genes, some of which may subsequently be shown to associate with both disease risk[64, 65] and other disease-related phenotypes.
Gene Expression Endophenotypes in Neurodegenerative Diseases-The Future
Available eQTL studies reviewed above illustrate the power of utilizing gene expression endophenotypes in conjunction with the disease phenotype, to identify novel genes, variants and pathways implicated in complex diseases. This approach, already applied to HDL-cholesterol levels and obesity, can provide an alternative to the current genetic mapping approaches in neurodegenerative diseases. Distinct from candidate gene and expression profiling studies, eQTL investigations in neurodegenerative diseases will aim to identify genetic loci that influence both gene expression and disease risk at the genome-wide level.
The motivation to exploit the gene expression endophenotype as an alternative in the genetic mapping of neurodegenerative disease loci stems in part from the relative shortcomings of GWAS of complex diseases. First, though multiple genetic loci have been identified for complex diseases through disease GWAS, these genetic variants fail to account for a substantial proportion of the underlying genetic risk[1, 9]. To give AD as an example, there are 12 published LOAD GWAS[1, 66], which led to the identification of a handful of genes that achieved genome-wide significance. These novel LOAD candidate loci collectively explain only a modest proportion of AD risk (reviewed in), despite evidence that genetics account for ~80% of the risk for AD. Gene expression endophenotypes appear to be more powerful than the disease phenotype, given the identification of hundreds to thousands of eQTLs in studies of only hundreds of subjects and may therefore identify genetic variants that may be missed by the classical disease mapping approach. Furthermore, they may allow for selecting a subset of the suggestive disease risk associations for further follow-up, by providing additional evidence for their functionality.
A relative concern with the eQTL approach is that since genetic studies may well be powered to detect eQTLs but underpowered to detect disease loci, complementary studies may be necessary to demonstrate that a particular eQTL also confers disease risk. Such studies may include concomitant mapping of other disease QTLs (e.g. serum/CSF biomarkers such as Aβ, neuropathology scores, cognitive measures); identification of eQTLs that are enriched for pathways previously implicated in disease; accepting less stringent significance for disease risk, but seeking replication in additional disease cohorts; functional studies with the identified genes to show their role in a disease-related in-vitro paradigm. It should also be noted that the eQTL approach can only capture those variants that confer disease risk via affecting gene expression and will miss coding changes that do not change transcript levels.
Another advantage of the eQTL approach is that the mechanism of action of the identified variants is already evident, which may allow for immediate downstream validation experiments such as measurement of protein levels in the brain or in-vitro transcriptional activity studies. Finally, combined eQTL and disease mapping approaches may identify the true disease-risk gene in regions of high LD spanning multiple genes, since the gene with the transcript that is influenced by the risk variant is most likely to be the disease-risk gene.
To achieve maximal benefit from the expression endophenotype approach in mapping novel neurodegenerative disease loci, attention should be given to technical aspects of the study design to minimize experimental confounds. These include choice and quality of tissue and transcriptome measuring platform, RNA integrity, statistical analysis to control for confounders and to detect significance[45, 58], to name a few.
It should be mentioned that in addition to eSNPs and coding variants, epigenetic mechanisms are also recognized as an additional source of influence for risk of complex diseases such as cancer and AD[68, 69]. Epigenetics comprise those reversible and dynamic mechanisms that influence gene expression, usually independent of DNA sequence, and include processes such as DNA methylation, histone acetylation and microRNAs[68, 69]. Seminal work in DNA methylation profiles from monozygotic and dizygotic twins provided evidence for an important role of epigenetics in heritability that is also environmentally regulated. Genome-wide epigenomic approach of investigating DNA-methylation differences in brains of subjects with schizophrenia and bipolar disorder identified epigenetic differences at numerous loci associated with psychosis, thereby providing evidence for the potential utility of this approach in CNS diseases. A comprehensive review of epigenomic approaches in gene discovery for neurodegenerative diseases is beyond the scope of this review. Suffice it to say that the combined eQTL and disease GWAS approaches advocated in this review can be further enhanced by additional incorporation of the epigenomic approach as reviewed elsewhere[68, 69].
Understanding the underlying genetic component of complex diseases, including Alzheimer's and other neurodegenerative diseases, has proven to be a challenge, despite the advances made mainly via GWAS of the dichotomous disease traits. Powerful approaches that constitute an alternative and are complementary to the current disease mapping algorithms are needed to overcome this challenge. Gene expression endophenotypes, which have a substantial genetic component, have already been used in mapping and functional validation of a few complex diseases and traits. While not all disease variants are expected to operate by changing transcript levels, it is expected that there will be many that confer disease risk by influencing gene expression. Such variants for neurodegenerative diseases may be captured by combining gene expression endophenotypes with existing disease GWAS to a) identify novel disease genes/pathways; b) validate suggestive findings from disease GWAS; c) elucidate the mechanism of action of newly discovered disease genes. Given their potential, gene expression endophenotypes are expected to be utilized in gene discovery for neurodegenerative diseases in the years to come.
Funding: R01032990 and P50AG016574 (PI: Ronald Petersen) to NET
Amyotrophic lateral sclerosis
Centre d'Etude du polymorphisme humain
Copy number variant
Expression genome-wide association study
Expression quantitative trait locus
False discovery rate
Frontotemporal dementia with parkinsonism linked to chromosome 17
Genome-wide association study(ies)
Kyoto Encyclopedia of Genes and Genomes
Lymphoblastoid cell line
Late-onset Alzheimer's disease
Mini Mental State Examination
Microtubule associated protein tau
Progressive supranuclear palsy
Quantitative polymerase chain reaction
Single nucleotide polymorphism
Ertekin-Taner N: Genetics of Alzheimer disease in the pre- and post-GWAS era. Alzheimers Res Ther. 2010, 2: 3-10.1186/alzrt26.
Gandhi S, Wood NW: Genome-wide association studies: the key to unlocking neurodegeneration?. Nat Neurosci. 2010, 13: 789-794. 10.1038/nn.2584.
Simon-Sanchez J, Singleton A: Genome-wide association studies in neurological disorders. Lancet Neurol. 2008, 7: 1067-1072. 10.1016/S1474-4422(08)70241-2.
Mackenzie IR, Rademakers R, Neumann M: TDP-43 and FUS in amyotrophic lateral sclerosis and frontotemporal dementia. Lancet Neurol. 2010, 9: 995-1007. 10.1016/S1474-4422(10)70195-2.
Sleegers K, Cruts M, Van Broeckhoven C: Molecular pathways of frontotemporal lobar degeneration. Annu Rev Neurosci. 2010, 71-88. 33:
A haplotype map of the human genome. Nature. 2005, 437: 1299-1320. 10.1038/nature04226.
Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, et al: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449: 851-861. 10.1038/nature06258.
Altshuler D, Daly MJ, Lander ES: Genetic mapping in human disease. Science. 2008, 322: 881-888. 10.1126/science.1156409.
Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al: Finding the missing heritability of complex diseases. Nature. 2009, 461: 747-753. 10.1038/nature08494.
John B, Lewis KR: Chromosome variability and geographic distribution in insects. Science. 1966, 152: 711-721. 10.1126/science.152.3723.711.
Shields J: Genetic theorizing and schizophrenia. Br J Psychiatry. 1973, 122: 15-30. 10.1192/bjp.122.1.15.
Zou F, Carrasquillo MM, Pankratz VS, Belbin O, Morgan K, Allen M, Wilcox SL, Ma L, Walker LP, Kouri N, et al: Gene expression levels as endophenotypes in genome-wide association studies of Alzheimer disease. Neurology. 2010, 74: 480-486. 10.1212/WNL.0b013e3181d07654.
Seshadri S, DeStefano AL, Au R, Massaro JM, Beiser AS, Kelly-Hayes M, Kase CS, D'Agostino RB, Decarli C, Atwood LD, Wolf PA: Genetic correlates of brain aging on MRI and cognitive test measures: a genome-wide association and linkage analysis in the Framingham Study. BMC Med Genet. 2007, 8 (Suppl 1): S15-10.1186/1471-2350-8-S1-S15.
Ikram MK, Xueling S, Jensen RA, Cotch MF, Hewitt AW, Ikram MA, Wang JJ, Klein R, Klein BE, Breteler MM, et al: Four novel Loci (19q13, 6q24, 12q24, and 5q14) influence the microcirculation in vivo. PLoS Genet. 2010, 6: e1001184-10.1371/journal.pgen.1001184.
Weidinger S, Baurecht H, Naumann A, Novak N: Genome-wide association studies on IgE regulation: are genetics of IgE also genetics of atopic disease?. Curr Opin Allergy Clin Immunol. 2010, 10: 408-417. 10.1097/ACI.0b013e32833d7d2d.
Bearden CE, Freimer NB: Endophenotypes for psychiatric disorders: ready for primetime?. Trends Genet. 2006, 22: 306-313. 10.1016/j.tig.2006.04.004.
Gould TD: The endophenotype concept in psychiatry: etymology and strategic intentions. Am J Psychiatry. 2003, 160: 636-645. 10.1176/appi.ajp.160.4.636.
Glahn DC, Thompson PM, Blangero J: Neuroimaging endophenotypes: strategies for finding genes influencing brain structure and function. Hum Brain Mapp. 2007, 28: 488-501. 10.1002/hbm.20401.
Ertekin-Taner N, Graff-Radford N, Younkin LH, Eckman C, Baker M, Adamson J, Ronald J, Blangero J, Hutton M, Younkin SG: Linkage of plasma Abeta42 to a quantitative locus on chromosome 10 in late-onset Alzheimer's disease pedigrees. Science. 2000, 290: 2303-2304. 10.1126/science.290.5500.2303.
Gould TD, Gottesman II: Psychiatric endophenotypes and the development of valid animal models. Genes Brain Behav. 2006, 5: 113-119. 10.1111/j.1601-183X.2005.00186.x.
Glahn DC, Almasy L, Blangero J, Burk GM, Estrada J, Peralta JM, Meyenberg N, Castro MP, Barrett J, Nicolini H, et al: Adjudicating neurocognitive endophenotypes for schizophrenia. Am J Med Genet B Neuropsychiatr Genet. 2007, 144: 242-249.
Ertekin-Taner N, Younkin LH, Yager DM, Parfitt F, Baker MC, Asthana S, Hutton ML, Younkin SG, Graff-Radford NR: Plasma amyloid beta protein is elevated in late-onset Alzheimer disease families. Neurology. 2008, 70: 596-606. 10.1212/01.wnl.0000278386.00035.21.
Ertekin-Taner N, Graff-Radford N, Younkin L, Eckman C, Adamson J, Schaid D, Blangero J, Hutton M, Younkin S: Heritability of plasma Aβ levels in typical late onset Alzheimer's disease pedigrees. Genetic Epidemiology. 2001, 21: 19-30. 10.1002/gepi.1015.
Jack CR, Knopman DS, Jagust WJ, Shaw LM, Aisen PS, Weiner MW, Petersen RC, Trojanowski JQ: Hypothetical model of dynamic biomarkers of the Alzheimer's pathological cascade. Lancet Neurol. 2010, 9: 119-128. 10.1016/S1474-4422(09)70299-6.
Brem RB, Yvert G, Clinton R, Kruglyak L: Genetic dissection of transcriptional regulation in budding yeast. Science. 2002, 296: 752-755. 10.1126/science.1069516.
Jansen RC, Nap JP: Genetical genomics: the added value from segregation. Trends Genet. 2001, 17: 388-391. 10.1016/S0168-9525(01)02310-1.
Dausset J, Cann H, Cohen D, Lathrop M, Lalouel JM, White R: Centre d'etude du polymorphisme humain (CEPH): collaborative genetic mapping of the human genome. Genomics. 1990, 6: 575-577. 10.1016/0888-7543(90)90491-C.
Yan H, Yuan W, Velculescu VE, Vogelstein B, Kinzler KW: Allelic variation in human gene expression. Science. 2002, 297: 1143-10.1126/science.1072545.
Schadt EE, Monks SA, Drake TA, Lusis AJ, Che N, Colinayo V, Ruff TG, Milligan SB, Lamb JR, Cavet G, et al: Genetics of gene expression surveyed in maize, mouse and man. Nature. 2003, 422: 297-302. 10.1038/nature01434.
Cheung VG, Conlin LK, Weber TM, Arcaro M, Jen KY, Morley M, Spielman RS: Natural variation in human gene expression assessed in lymphoblastoid cells. Nat Genet. 2003, 33: 422-425. 10.1038/ng1094.
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: 743-747. 10.1038/nature02797.
Monks SA, Leonardson A, Zhu H, Cundiff P, Pietrusiak P, Edwards S, Phillips JW, Sachs A, Schadt EE: Genetic inheritance of gene expression in human cell lines. Am J Hum Genet. 2004, 75: 1094-1105. 10.1086/426461.
Cheung VG, Spielman RS, Ewens KG, Weber TM, Morley M, Burdick JT: Mapping determinants of human gene expression by regional and genome-wide association. Nature. 2005, 437: 1365-1369. 10.1038/nature04244.
Stranger BE, Forrest MS, Clark AG, Minichiello MJ, Deutsch S, Lyle R, Hunt S, Kahl B, Antonarakis SE, Tavare S, et al: Genome-wide associations of gene expression variation in humans. PLoS Genet. 2005, 1: e78-10.1371/journal.pgen.0010078.
Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, Ingle CE, Dunning M, Flicek P, Koller D, et al: Population genomics of human gene expression. Nat Genet. 2007, 39: 1217-1224. 10.1038/ng2142.
Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird CP, de Grassi A, Lee C, et al: Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007, 315: 848-853. 10.1126/science.1136678.
Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, Wong KC, Taylor J, Burnett E, Gut I, Farrall M, et al: A genome-wide association study of global gene expression. Nat Genet. 2007, 39: 1202-1207. 10.1038/ng2109.
Goring HH, Curran JE, Johnson MP, Dyer TD, Charlesworth J, Cole SA, Jowett JB, Abraham LJ, Rainwater DL, Comuzzie AG, et al: Discovery of expression QTLs using large-scale transcriptional profiling in human lymphocytes. Nat Genet. 2007, 39: 1208-1216. 10.1038/ng2119.
Emilsson V, Thorleifsson G, Zhang B, Leonardson AS, Zink F, Zhu J, Carlson S, Helgason A, Walters GB, Gunnarsdottir S, et al: Genetics of gene expression and its effect on disease. Nature. 2008, 452: 423-428. 10.1038/nature06758.
Schadt EE, Molony C, Chudin E, Hao K, Yang X, Lum PY, Kasarskis A, Zhang B, Wang S, Suver C, et al: Mapping the genetic architecture of gene expression in human liver. PLoS Biol. 2008, 6: e107-10.1371/journal.pbio.0060107.
Myers AJ, Gibbs JR, Webster JA, Rohrer K, Zhao A, Marlowe L, Kaleem M, Leung D, Bryden L, Nath P, et al: A survey of genetic human cortical gene expression. Nat Genet. 2007, 39: 1494-1499. 10.1038/ng.2007.16.
Webster JA, Gibbs JR, Clarke J, Ray M, Zhang W, Holmans P, Rohrer K, Zhao A, Marlowe L, Kaleem M, et al: Genetic control of human brain transcript expression in Alzheimer disease. Am J Hum Genet. 2009, 84: 445-458. 10.1016/j.ajhg.2009.03.011.
Myers AJ, Pittman AM, Zhao AS, Rohrer K, Kaleem M, Marlowe L, Lees A, Leung D, McKeith IG, Perry RH, et al: The MAPT H1c risk haplotype is associated with increased expression of tau and especially of 4 repeat containing transcripts. Neurobiol Dis. 2007, 25: 561-570. 10.1016/j.nbd.2006.10.018.
Stranger BE, Forrest MS, Clark AG, Minichiello MJ, Deutsch S, Lyle R, Hunt S, Kahl B, Antonarakis SE, Tavare S, et al: Genome-wide associations of gene expression variation in humans. PLoS Genet. 2005, 1: e78-10.1371/journal.pgen.0010078.
Atz M, Walsh D, Cartagena P, Li J, Evans S, Choudary P, Overman K, Stein R, Tomita H, Potkin S, et al: Methodological considerations for gene expression profiling of human brain. J Neurosci Methods. 2007, 163: 295-309. 10.1016/j.jneumeth.2007.03.022.
Li YJ, Oliveira SA, Xu P, Martin ER, Stenger JE, Scherzer CR, Hauser MA, Scott WK, Small GW, Nance MA, et al: Glutathione S-transferase omega-1 modifiesage-at-onset of Alzheimer disease and Parkinson disease. Hum Mol Genet. 2003, 12: 3259-3267. 10.1093/hmg/ddg357.
Li YJ, Scott WK, Zhang L, Lin PI, Oliveira SA, Skelly T, Doraiswamy MP, Welsh-Bohmer KA, Martin ER, Haines JL, et al: Revealing the role of glutathione S-transferase omega in age-at-onset of Alzheimer and Parkinson diseases. Neurobiol Aging. 2006, 27: 1087-1093. 10.1016/j.neurobiolaging.2005.05.013.
Corder EH, Saunders AM, Strittmatter WJ, Schmechel DE, Gaskell PC, Small GW, Roses AD, Haines JL, Pericak-Vance MA: Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. Science. 1993, 261: 921-923. 10.1126/science.8346443.
Ertekin-Taner N: Genetics of Alzheimer's disease: a centennial review. Neurol Clin. 2007, 25: 611-667. 10.1016/j.ncl.2007.03.009.
Laws SM, Hone E, Gandy S, Martins RN: Expanding the association between the APOE gene and the risk of Alzheimer's disease: possible roles for APOE promoter polymorphisms and alterations in APOE transcription. J Neurochem. 2003, 84: 1215-1236. 10.1046/j.1471-4159.2003.01615.x.
Bullido MJ, Artiga MJ, Recuero M, Sastre I, Garcia MA, Aldudo J, Lendon C, Han SW, Morris JC, Frank A, et al: A polymorphism in the regulatory region of APOE associated with risk for Alzheimer's dementia. Nat Genet. 1998, 18: 69-71. 10.1038/ng0198-69.
Wang JC, Kwon JM, Shah P, Morris JC, Goate A: Effect of APOE genotype and promoter polymorphism on risk of Alzheimer's disease. Neurology. 2000, 55: 1644-1649.
Holtzman DM, Bales KR, Tenkova T, Fagan AM, Parsadanian M, Sartorius LJ, Mackey B, Olney J, McKeel D, Wozniak D, Paul SM: Apolipoprotein E isoform-dependent amyloid deposition and neuritic degeneration in a mouse model of Alzheimer's disease. Proc. Natl Acad Sci USA. 2000, 97: 2892-2897. 10.1073/pnas.050004797.
Pittman AM, Fung HC, de Silva R: Untangling the tau gene association with neurodegenerative disorders. Hum Mol Genet. 2006, 15 (Spec No 2): R188-195.
Caffrey TM, Joachim C, Paracchini S, Esiri MM, Wade-Martins R: Haplotype-specific expression of exon 10 at the human MAPT locus. Hum Mol Genet. 2006, 15: 3529-3537. 10.1093/hmg/ddl429.
Caffrey TM, Joachim C, Wade-Martins R: Haplotype-specific expression of the N-terminal exons 2 and 3 at the human MAPT locus. Neurobiol Aging. 2008, 29: 1923-1929. 10.1016/j.neurobiolaging.2007.05.002.
Carrasquillo MM, Belbin O, Zou F, Allen M, Ertekin-Taner N, Ansari M, Wilcox SL, Kashino MR, Ma L, Younkin LH, et al: Concordant association of insulin degrading enzyme gene (IDE) variants with IDE mRNA, Abeta, and Alzheimer's disease. PLoS One. 2010, 5: e8764-10.1371/journal.pone.0008764.
Verducci JS, Melfi VF, Lin S, Wang Z, Roy S, Sen CK: Microarray analysis of gene expression: considerations in data mining and statistical treatment. Physiol Genomics. 2006, 25: 355-363. 10.1152/physiolgenomics.00314.2004.
Folstein MF, Folstein SE, McHugh PR: "Mini-mental state". A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res. 1975, 12: 189-198. 10.1016/0022-3956(75)90026-6.
Blalock EM, Geddes JW, Chen KC, Porter NM, Markesbery WR, Landfield PW: Incipient Alzheimer's disease: microarray correlation analyses reveal major transcriptional and tumor suppressor responses. Proc Natl Acad Sci USA. 2004, 101: 2173-2178. 10.1073/pnas.0308512100.
Bossers K, Wirz KT, Meerhoff GF, Essing AH, van Dongen JW, Houba P, Kruse CG, Verhaagen J, Swaab DF: Concerted changes in transcripts in the prefrontal cortex precede neuropathology in Alzheimer's disease. Brain. 2010, 133: 3699-3723. 10.1093/brain/awq258.
Braak H, Braak E: Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol (Berl). 1991, 82: 239-259. 10.1007/BF00308809.
Bronner IF, Bochdanovits Z, Rizzu P, Kamphorst W, Ravid R, van Swieten JC, Heutink P: Comprehensive mRNA expression profiling distinguishes tauopathies and identifies shared molecular pathways. PLoS One. 2009, 4: e6826-10.1371/journal.pone.0006826.
Taguchi K, Yamagata HD, Zhong W, Kamino K, Akatsu H, Hata R, Yamamoto T, Kosaka K, Takeda M, Kondo I, Miki T: Identification of hippocampus-related candidate genes for Alzheimer's disease. Ann Neurol. 2005, 57: 585-588. 10.1002/ana.20433.
Chapuis J, Hot D, Hansmannel F, Kerdraon O, Ferreira S, Hubans C, Maurage CA, Huot L, Bensemain F, Laumet G, et al: Transcriptomic and genetic studies identify IL-33 as a candidate gene for Alzheimer's disease. Mol Psychiatry. 2009, 14: 1004-1016. 10.1038/mp.2009.10.
Seshadri S, Fitzpatrick AL, Ikram MA, DeStefano AL, Gudnason V, Boada M, Bis JC, Smith AV, Carassquillo MM, Lambert JC, et al: Genome-wide analysis of genetic loci associated with Alzheimer disease. JAMA. 2010, 303: 1832-1840. 10.1001/jama.2010.574.
Bergem AL, Engedal K, Kringlen E: The role of heredity in late-onset Alzheimer disease and vascular dementia. A twin study. Arch Gen Psychiatry. 1997, 54: 264-270.
Mastroeni D, Grover A, Delvaux E, Whiteside C, Coleman PD, Rogers J: Epigenetic mechanisms in Alzheimer's disease. Neurobiol Aging. 2011.
Mill J: Toward an integrated genetic and epigenetic approach to Alzheimer's disease. Neurobiol Aging. 2011.
Kaminsky ZA, Tang T, Wang SC, Ptak C, Oh GH, Wong AH, Feldcamp LA, Virtanen C, Halfvarson J, Tysk C, et al: DNA methylation profiles in monozygotic and dizygotic twins. Nat Genet. 2009, 41: 240-245. 10.1038/ng.286.
Mill J, Tang T, Kaminsky Z, Khare T, Yazdanpanah S, Bouchard L, Jia P, Assadzadeh A, Flanagan J, Schumacher A, et al: Epigenomic profiling reveals DNA-methylation changes associated with major psychosis. Am J Hum Genet. 2008, 82: 696-711. 10.1016/j.ajhg.2008.01.008.
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to 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 cited.
About this article
Cite this article
Ertekin-Taner, N. Gene expression endophenotypes: a novel approach for gene discovery in Alzheimer's disease. Mol Neurodegeneration 6, 31 (2011). https://doi.org/10.1186/1750-1326-6-31
- Progressive Supranuclear Palsy
- Lymphoblastoid Cell Line
- Human Gene Expression
- eQTL Study
- CEPH Family