Genome-wide association studies in frontotemporal dementia showed limited success in identifying associated loci. This is possibly due to small sample size, allelic heterogeneity, small effect sizes of single genetic variants, and the necessity to statistically correct for testing millions of genetic variants. To overcome these issues, we performed gene-based association studies on 3348 clinically identified frontotemporal dementia cases and 9390 controls (discovery, replication and joint-cohort analyses). We report association of APOE and TOMM40 with behavioural variant frontotemporal dementia, and ARHGAP35 and SERPINA1 with progressive non-fluent aphasia. Further, we found the ɛ2 and ɛ4 alleles of APOE harbouring protective and risk increasing effects, respectively, in clinical subtypes of frontotemporal dementia against neurologically normal controls. The APOE-locus association with behavioural variant frontotemporal dementia indicates its potential risk-increasing role across different neurodegenerative diseases, whereas the novel genetic associations of ARHGAP35 and SERPINA1 with progressive non-fluent aphasia point towards a potential role of the stress-signalling pathway in its pathophysiology.

Introduction

Frontotemporal dementia (FTD) is one of the leading causes of dementia in patients younger than 65 years of age (Rabinovici and Miller, 2010; Seelaar et al., 2011). It is characterized by degeneration of the frontal and anterior temporal lobes leading to a decline in behaviour and language. FTD is a heterogeneous condition clinically, pathologically and genetically (Cairns et al., 2007; Seelaar et al., 2011). Clinically, it is broadly categorized into the behavioural variant (bvFTD) and the language variant or primary progressive aphasia (PPA), which is further categorized into semantic dementia and progressive non-fluent aphasia (PNFA). There is frequent overlap between the presence of FTD and a number of motor diseases such as parkinsonian disorders, corticobasal syndrome, progressive supranuclear palsy and motor neuron disease (FTD-MND) (Rohrer et al., 2009). The underlying pathological spectrum of FTD, termed frontal temporal lobar degeneration (FTLD), is based on neuronal lesions and protein inclusions such as with tau or TAR-DNA binding protein (TDP)-43 pathology. Besides the Mendelian genes MAPT, GRN and C9orf72 that are causal in up to ∼30–50% of familial FTLD cases, rare variability in few other genes has been implicated in <5% of cases (Rohrer et al., 2009; Ferrari et al., 2013; Woollacott and Rohrer, 2016). To date, few large genome wide association studies (GWAS) have been performed for FTD (Van Deerlin et al., 2010; Ferrari et al., 2014, 2015) reporting an association with TMEM106B for FTLD with TDP-43 pathology (Van Deerlin et al., 2010), and with the locus comprising RAB38 and CTSC as well as the HLA-DRA/HLA-DRB5 locus for bvFTD and FTD, respectively (Ferrari et al., 2014).

In a typical GWAS, an association test on a single variant [single nucleotide polymorphisms (SNPs) or Indels] is performed to map genes associated with a phenotype; however, many independent risk alleles for a given phenotype can be localized within a gene (Yang et al., 2012; Tada et al., 2014; Zhang et al., 2014). Hence a classical GWAS approach will be less powered to detect genes containing many independent risk alleles (Hagg et al., 2015). A joint-variant gene-based test that combines independent association signals within a gene while accounting for the linkage disequilibrium (LD) between variants can overcome this limitation. A number of approaches have been reported to perform joint-SNP gene-based analysis: the permutation test—where empirical evidence of association of the combined test statistics is calculated by shuffling the samples while keeping markers intact—is currently considered the golden standard (Liu et al., 2010). However, the requirement of genotype data and computational burden limits its use. Recently, our group developed a new approach called Multi-marker Analysis of GenoMic Annotation (MAGMA) that uses a multiple regression model to perform joint-SNP gene-based analysis using GWAS summary data (de Leeuw et al., 2015).

In this study, we performed a hypothesis free gene-wide association study on FTD subtypes (bvFTD, semantic dementia, PNFA and FTD-MND) using GWAS summary files obtained from the International FTD-Genomics consortium (IFGC) (Ferrari et al., 2014). We used the MAGMA software to perform the gene-based analysis. We report results of discovery, replication and combined cohort analyses for each FTD subtype; we also assessed individual risk variants for associated genes, which can be used for replication in the individual variant genotype setting.

Materials and methods

Samples

The dataset used in the FTD-GWAS was described previously (Ferrari et al., 2014). Briefly, 44 international groups contributed clinical FTD samples. Patients were diagnosed according to the Neary criteria or the revised criteria for bvFTD and language variants of FTD (Neary et al., 1998; Gorno-Tempini et al., 2011; Rascovsky et al., 2011). Approximately 3% of cases were pathologically confirmed. To each individual case two ancestry and age-matched neurologically normal controls were assigned. For the current study we used the GWAS summary datasets of the discovery and replication cohorts of each FTD subtype, bvFTD (discovery: 1377 cases and 2754 controls; replication: 690 cases and 5092 controls), PNFA (discovery: 269 cases and 538 controls; replication: 221 cases and 5092 controls), semantic dementia (discovery: 308 cases and 616 controls; replication: 189 cases and 5092 controls) and FTD-MND (discovery: 200 cases and 400 controls; replication: 94 cases and 5092 controls). Overall bvFTD, PNFA, semantic dementia and FTD-MND constituted 61.74%, 14.64%, 14.85% and 8.78% of the total FTD cases, respectively.

Statistical analysis

We performed the joint-SNP gene-based analysis using MAGMA (de Leeuw et al., 2015). The MAGMA approach is based on a multiple linear principal components regression model. By projecting the multivariate LD matrix of SNPs in a gene it first extracts principal components that explain genetic variation. These principal components are further used as predictors of a phenotype under a linear regression framework. MAGMA then uses Fisher’s test to compute P-values to test association between a gene and the phenotype.

We used 19 418 hg19 annotated protein-coding genes to perform the analysis. As all the samples involved were of European descent, we used the 1000 Genomes phase 1 European reference population to estimate LD between variants (1000 Genomes Project Consortium et al., 2012). We only considered SNPs in the 5’- and 3’-untranslated region (UTR) and the open reading frame for the joint-SNP gene-based tests. This strategy resulted in loss of cis- regulatory variants, but was more stringent and open reading frame-specific.

The schematic representation of the strategy for multi-stage gene-wide association analysis for FTD and subtypes is described in Supplementary Fig. 1. We performed separate gene-based tests using discovery and replication datasets for each FTD subtype reported previously by the IFGC (Ferrari et al., 2014). We performed gene-based tests using only those variants that were either genotyped or imputed with imputation score >0.50. Moreover, we only considered common variants with minor allele frequency >0.01. For individual FTD subtypes, 4 303 460 and 55 375 variants were available for the gene-based analysis in the discovery and replication cohorts, respectively. Further, 16 313 and 10 349 genes that contained at least one variant within the 5’-, 3’-UTR and open reading frame, were tested for association with a given FTD subtype in the discovery and replication cohorts, respectively. To identify additional genes associated with individual FTD subtypes, we meta-analysed the gene-based P-values obtained in the discovery and replication cohorts using the Stouffer’s combination approach for the sample size weighted combination of P-values. For each FTD subtype we tested association of total 16 920 genes either in the discovery or replication cohorts. To correct for multiple association tests performed for 16 920 genes with one of the four subtypes of FTD, we applied the conservative Bonferroni correction method establishing a gene-wide significance threshold at 7.388 × 10−7 [= 0.05 / (4 × 16 920)]. To identify genes associated with any FTD subtype, we combined the gene-based test statistics for either subtype (bvFTD, semantic dementia, PNFA and FTD-MND) using the sample size weighted Stouffer’s combination method.

Functional characterization of associated genes

We downloaded the gene expression profiles of the associated genes across 13 human brain tissues (in alphabetical order: amygdala, anterior cingulate cortex, caudate, cerebral hemisphere, cerebellum, cortex, frontal cortex, hippocampus, nucleus accumbens, putamen, spinal cord, substantia nigra) using the GTEx portal (Mele et al., 2015) (GTEx Analysis V6 dbGaP Accession phs000424.v6.p1). We also investigated the functional annotations of variants that are in LD (r2 > 0.8 in the 1000 Genomes phase 1 European panel) with SNPs used in deriving gene-based P-values of the associated genes using software HaploReg (Ward and Kellis, 2012) (version 4.1) and RegulomeDB (Boyle et al., 2012) (version 1.1).

Results

Associations with frontotemporal dementia and its subtypes

Behavioural variant frontotemporal dementia

In the discovery cohort, two genes passed the gene-wide significant P-value threshold 7.39 × 10−7: TOMM40 (P = 5.79 × 10−8) and APOE (P = 1.37 × 10−7). In the replication cohort the P-values were 6.40 × 10−5 for TOMM40 and 1.69 × 10−3 for APOE suggesting consistency of associations across independent bvFTD samples. No other genes passed the significance threshold for the bvFTD subtype (Supplementary Fig. 2A).

Interestingly, in the discovery cohort the SNPs rs7412 and rs429358, which determine three epsilon (ɛ) alleles ɛ2, ɛ3 and ɛ4 of the APOE gene, were among the SNPs driving its association with bvFTD, with P-values 0.023 (rs7412) and 5.04 × 10−6 (rs429358). In the replication cohort rs429358 was not genotyped, whereas information on rs769449, an intronic variant in high LD (r2 = 0.82, 1000 Genomes phase 1 European population) with rs429358, was available; here the P-values were 0.222 for rs7412 and 1.95 × 10−4 for rs769449.

To check whether the association of TOMM40 with bvFTD was independent of the epsilon variants, we reperformed the gene-based test on TOMM40 gene using only those variants in negligible LD (r2 < 0.2) with rs7412 and rs429358 in 1000 Genomes phase 1 European panel. This analysis showed moderate association of TOMM40 with bvFTD (P = 7.51 × 10−6; Table 1) suggesting that the TOMM40 gene harbours signals for the risk of bvFTD that are independent of the epsilon alleles of APOE gene.

Table 1

Top two genes for each FTD subtype observed in combined analyses

SymbolChrStartStopnSNPs discovery cohortnSNPs replicationbvFTDPNFASDFTD-MNDFTD-Any
Discovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueMeta- analysed P-value across all cohorts
BvFTD
TOMM4019453944774540694629135.79 × 10−86.40 × 10−55.43 × 10–11*0.9460.0230.4910.1070.1440.0520.4320.5040.444.98 × 10–11*
APOE194540903945412650521.37 × 10−71.655 × 10−31.69 × 10−9*0.9989.84 × 10−40.6190.0460.0276.48 × 10−30.0610.5840.0962.18 × 10–10*
TOMM40 after removing variants in LD (r2 ≥ 0.2) with the epsilon variants1945394477454069461851.95 × 10−50.0737.51 × 10−6NANANANANANANANANANA
PNFA
ARHGAP3519474219334750833425120.7440.130.5334.26 × 10−67.16 × 10−32.95 × 10−7*0.6230.3050.4860.4160.1180.2430.139
SERPINA11494843084948570295080.9580.4240.9283.71 × 10−51.19 × 10−33.02 × 10−7*0.8730.3550.7650.4430.0290.1750.637
SD
WDR6612122356463122441833149110.4700.2320.3470.6530.3360.5149.89 × 10−32.30 × 10−51.10 10−50.1880.3360.1630.071
TTC2723285308733046118353240.4230.9690.7440.7890.9350.9433.97 10−60.4569.75 10−50.1030.6800.1720.489
FTD-MND
C9orf729275465442757384275250.0930.7690.1960.9890.2770.9190.7840.2520.6031.97 × 10−59.95 × 10−31.23 × 10−6**0.139
IFNK92752431227526496940.0540.9040.1970.9520.4620.8900.9930.3040.9578.80 × 10−60.0391.77 × 10−6**0.220
SymbolChrStartStopnSNPs discovery cohortnSNPs replicationbvFTDPNFASDFTD-MNDFTD-Any
Discovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueMeta- analysed P-value across all cohorts
BvFTD
TOMM4019453944774540694629135.79 × 10−86.40 × 10−55.43 × 10–11*0.9460.0230.4910.1070.1440.0520.4320.5040.444.98 × 10–11*
APOE194540903945412650521.37 × 10−71.655 × 10−31.69 × 10−9*0.9989.84 × 10−40.6190.0460.0276.48 × 10−30.0610.5840.0962.18 × 10–10*
TOMM40 after removing variants in LD (r2 ≥ 0.2) with the epsilon variants1945394477454069461851.95 × 10−50.0737.51 × 10−6NANANANANANANANANANA
PNFA
ARHGAP3519474219334750833425120.7440.130.5334.26 × 10−67.16 × 10−32.95 × 10−7*0.6230.3050.4860.4160.1180.2430.139
SERPINA11494843084948570295080.9580.4240.9283.71 × 10−51.19 × 10−33.02 × 10−7*0.8730.3550.7650.4430.0290.1750.637
SD
WDR6612122356463122441833149110.4700.2320.3470.6530.3360.5149.89 × 10−32.30 × 10−51.10 10−50.1880.3360.1630.071
TTC2723285308733046118353240.4230.9690.7440.7890.9350.9433.97 10−60.4569.75 10−50.1030.6800.1720.489
FTD-MND
C9orf729275465442757384275250.0930.7690.1960.9890.2770.9190.7840.2520.6031.97 × 10−59.95 × 10−31.23 × 10−6**0.139
IFNK92752431227526496940.0540.9040.1970.9520.4620.8900.9930.3040.9578.80 × 10−60.0391.77 × 10−6**0.220

*P-values less than the gene-wide significant threshold 7.388 × 10−7 correcting multiple tests for association of 16 920 genes with one of the four subtypes of FTD.

**P-values less than the Bonferroni corrected P-value threshold 2.955 × 10−6 for association tests of 16 920 genes with a single phenotype.

SD = semantic dementia.

Table 1

Top two genes for each FTD subtype observed in combined analyses

SymbolChrStartStopnSNPs discovery cohortnSNPs replicationbvFTDPNFASDFTD-MNDFTD-Any
Discovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueMeta- analysed P-value across all cohorts
BvFTD
TOMM4019453944774540694629135.79 × 10−86.40 × 10−55.43 × 10–11*0.9460.0230.4910.1070.1440.0520.4320.5040.444.98 × 10–11*
APOE194540903945412650521.37 × 10−71.655 × 10−31.69 × 10−9*0.9989.84 × 10−40.6190.0460.0276.48 × 10−30.0610.5840.0962.18 × 10–10*
TOMM40 after removing variants in LD (r2 ≥ 0.2) with the epsilon variants1945394477454069461851.95 × 10−50.0737.51 × 10−6NANANANANANANANANANA
PNFA
ARHGAP3519474219334750833425120.7440.130.5334.26 × 10−67.16 × 10−32.95 × 10−7*0.6230.3050.4860.4160.1180.2430.139
SERPINA11494843084948570295080.9580.4240.9283.71 × 10−51.19 × 10−33.02 × 10−7*0.8730.3550.7650.4430.0290.1750.637
SD
WDR6612122356463122441833149110.4700.2320.3470.6530.3360.5149.89 × 10−32.30 × 10−51.10 10−50.1880.3360.1630.071
TTC2723285308733046118353240.4230.9690.7440.7890.9350.9433.97 10−60.4569.75 10−50.1030.6800.1720.489
FTD-MND
C9orf729275465442757384275250.0930.7690.1960.9890.2770.9190.7840.2520.6031.97 × 10−59.95 × 10−31.23 × 10−6**0.139
IFNK92752431227526496940.0540.9040.1970.9520.4620.8900.9930.3040.9578.80 × 10−60.0391.77 × 10−6**0.220
SymbolChrStartStopnSNPs discovery cohortnSNPs replicationbvFTDPNFASDFTD-MNDFTD-Any
Discovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueDiscovery cohort P-valueReplication cohort P-valueCombined P-valueMeta- analysed P-value across all cohorts
BvFTD
TOMM4019453944774540694629135.79 × 10−86.40 × 10−55.43 × 10–11*0.9460.0230.4910.1070.1440.0520.4320.5040.444.98 × 10–11*
APOE194540903945412650521.37 × 10−71.655 × 10−31.69 × 10−9*0.9989.84 × 10−40.6190.0460.0276.48 × 10−30.0610.5840.0962.18 × 10–10*
TOMM40 after removing variants in LD (r2 ≥ 0.2) with the epsilon variants1945394477454069461851.95 × 10−50.0737.51 × 10−6NANANANANANANANANANA
PNFA
ARHGAP3519474219334750833425120.7440.130.5334.26 × 10−67.16 × 10−32.95 × 10−7*0.6230.3050.4860.4160.1180.2430.139
SERPINA11494843084948570295080.9580.4240.9283.71 × 10−51.19 × 10−33.02 × 10−7*0.8730.3550.7650.4430.0290.1750.637
SD
WDR6612122356463122441833149110.4700.2320.3470.6530.3360.5149.89 × 10−32.30 × 10−51.10 10−50.1880.3360.1630.071
TTC2723285308733046118353240.4230.9690.7440.7890.9350.9433.97 10−60.4569.75 10−50.1030.6800.1720.489
FTD-MND
C9orf729275465442757384275250.0930.7690.1960.9890.2770.9190.7840.2520.6031.97 × 10−59.95 × 10−31.23 × 10−6**0.139
IFNK92752431227526496940.0540.9040.1970.9520.4620.8900.9930.3040.9578.80 × 10−60.0391.77 × 10−6**0.220

*P-values less than the gene-wide significant threshold 7.388 × 10−7 correcting multiple tests for association of 16 920 genes with one of the four subtypes of FTD.

**P-values less than the Bonferroni corrected P-value threshold 2.955 × 10−6 for association tests of 16 920 genes with a single phenotype.

SD = semantic dementia.

Table 2

The odds of clinical subtypes of FTD and epsilon alleles, and homozygous ɛ4/ɛ4 genotype

ɛ2 alleleɛ3 alleleɛ4 alleleɛ4/ɛ4 genotype
Case cohortORCIP-valueORCIP-valueORCIP-valueORCIP-value
BvFTD discovery0.8340.696–0.9940.0410.9080.817–1.0110.0761.3411.182–1.5175.00 × 10−6*1.7371.090–2.6820.018
BvFTD replication0.7710.598–0.9800.0310.9430.815–1.0940.4311.1540.961–1.3800.1101.4070.679–2.6370.270
BvFTD combined0.7720.664–0.8943.88 × 10−4*0.9200.841–1.0060.0661.2781.147–1.4218.14 × 10−6*1.6271.089–2.3840.012
PNFA discovery0.7920.526–1.1500.2610.9500.755–1.2050.6391.2460.938–1.6310.1202.3200.902–5.0150.040
PNFA replication0.6020.358–0.9560.0280.9710.754–1.2660.7961.3631.006–1.8160.0412.4230.860 5.5370.046
PNFA combined0.7060.515–0.9470.0190.960.808–1.1450.6281.2981.056–1.5840.0112.3671.211–4.2618.52 × 10−3**
SD discovery0.8080.554–1.1430.2580.7640.625–0.9400.0101.6491.300–2.0734.33 × 10−5*2.6141.152–5.2180.011
SD replication0.4020.198–0.7309.47 × 10−4*1.2170.906–1.6640.2081.1120.777–1.5520.5371.8780.497–5.0310.176
SD combined0.6510.470–0.8803.64 × 10−3*0.8980.759–1.0660.2051.4381.181–1.7412.93 × 10−4*2.3331.194–4.1999.08 × 10−3**
FTD–MND discovery0.5220.289–0.8769.30 × 10−31.0420.793–1.3880.8381.3110.948–1.780.0871.3240.266–4.0280.500
FTD–MND replication0.6740.303–1.3130.3111.1340.754–1.7640.6220.9340.53–1.5440.9010.0000.000–3.5360.629
FTD–MND combined0.5710.361–0.8626.01 × 10−31.0700.852–1.3570.6121.1880.901–1.5440.2020.8950.181–2.7141
ɛ2 alleleɛ3 alleleɛ4 alleleɛ4/ɛ4 genotype
Case cohortORCIP-valueORCIP-valueORCIP-valueORCIP-value
BvFTD discovery0.8340.696–0.9940.0410.9080.817–1.0110.0761.3411.182–1.5175.00 × 10−6*1.7371.090–2.6820.018
BvFTD replication0.7710.598–0.9800.0310.9430.815–1.0940.4311.1540.961–1.3800.1101.4070.679–2.6370.270
BvFTD combined0.7720.664–0.8943.88 × 10−4*0.9200.841–1.0060.0661.2781.147–1.4218.14 × 10−6*1.6271.089–2.3840.012
PNFA discovery0.7920.526–1.1500.2610.9500.755–1.2050.6391.2460.938–1.6310.1202.3200.902–5.0150.040
PNFA replication0.6020.358–0.9560.0280.9710.754–1.2660.7961.3631.006–1.8160.0412.4230.860 5.5370.046
PNFA combined0.7060.515–0.9470.0190.960.808–1.1450.6281.2981.056–1.5840.0112.3671.211–4.2618.52 × 10−3**
SD discovery0.8080.554–1.1430.2580.7640.625–0.9400.0101.6491.300–2.0734.33 × 10−5*2.6141.152–5.2180.011
SD replication0.4020.198–0.7309.47 × 10−4*1.2170.906–1.6640.2081.1120.777–1.5520.5371.8780.497–5.0310.176
SD combined0.6510.470–0.8803.64 × 10−3*0.8980.759–1.0660.2051.4381.181–1.7412.93 × 10−4*2.3331.194–4.1999.08 × 10−3**
FTD–MND discovery0.5220.289–0.8769.30 × 10−31.0420.793–1.3880.8381.3110.948–1.780.0871.3240.266–4.0280.500
FTD–MND replication0.6740.303–1.3130.3111.1340.754–1.7640.6220.9340.53–1.5440.9010.0000.000–3.5360.629
FTD–MND combined0.5710.361–0.8626.01 × 10−31.0700.852–1.3570.6121.1880.901–1.5440.2020.8950.181–2.7141

SD = semantic dementia.

*P-values less than significant threshold 4.17 × 10−3 correcting multiple tests for association of three epsilon alleles with one of the four subtypes of FTD.

**P-values less than significant threshold 1.25 × 10−2 correcting multiple tests for association of homozygous ε4/ε4 genotypes with one of the four subtypes of FTD.

Table 2

The odds of clinical subtypes of FTD and epsilon alleles, and homozygous ɛ4/ɛ4 genotype

ɛ2 alleleɛ3 alleleɛ4 alleleɛ4/ɛ4 genotype
Case cohortORCIP-valueORCIP-valueORCIP-valueORCIP-value
BvFTD discovery0.8340.696–0.9940.0410.9080.817–1.0110.0761.3411.182–1.5175.00 × 10−6*1.7371.090–2.6820.018
BvFTD replication0.7710.598–0.9800.0310.9430.815–1.0940.4311.1540.961–1.3800.1101.4070.679–2.6370.270
BvFTD combined0.7720.664–0.8943.88 × 10−4*0.9200.841–1.0060.0661.2781.147–1.4218.14 × 10−6*1.6271.089–2.3840.012
PNFA discovery0.7920.526–1.1500.2610.9500.755–1.2050.6391.2460.938–1.6310.1202.3200.902–5.0150.040
PNFA replication0.6020.358–0.9560.0280.9710.754–1.2660.7961.3631.006–1.8160.0412.4230.860 5.5370.046
PNFA combined0.7060.515–0.9470.0190.960.808–1.1450.6281.2981.056–1.5840.0112.3671.211–4.2618.52 × 10−3**
SD discovery0.8080.554–1.1430.2580.7640.625–0.9400.0101.6491.300–2.0734.33 × 10−5*2.6141.152–5.2180.011
SD replication0.4020.198–0.7309.47 × 10−4*1.2170.906–1.6640.2081.1120.777–1.5520.5371.8780.497–5.0310.176
SD combined0.6510.470–0.8803.64 × 10−3*0.8980.759–1.0660.2051.4381.181–1.7412.93 × 10−4*2.3331.194–4.1999.08 × 10−3**
FTD–MND discovery0.5220.289–0.8769.30 × 10−31.0420.793–1.3880.8381.3110.948–1.780.0871.3240.266–4.0280.500
FTD–MND replication0.6740.303–1.3130.3111.1340.754–1.7640.6220.9340.53–1.5440.9010.0000.000–3.5360.629
FTD–MND combined0.5710.361–0.8626.01 × 10−31.0700.852–1.3570.6121.1880.901–1.5440.2020.8950.181–2.7141
ɛ2 alleleɛ3 alleleɛ4 alleleɛ4/ɛ4 genotype
Case cohortORCIP-valueORCIP-valueORCIP-valueORCIP-value
BvFTD discovery0.8340.696–0.9940.0410.9080.817–1.0110.0761.3411.182–1.5175.00 × 10−6*1.7371.090–2.6820.018
BvFTD replication0.7710.598–0.9800.0310.9430.815–1.0940.4311.1540.961–1.3800.1101.4070.679–2.6370.270
BvFTD combined0.7720.664–0.8943.88 × 10−4*0.9200.841–1.0060.0661.2781.147–1.4218.14 × 10−6*1.6271.089–2.3840.012
PNFA discovery0.7920.526–1.1500.2610.9500.755–1.2050.6391.2460.938–1.6310.1202.3200.902–5.0150.040
PNFA replication0.6020.358–0.9560.0280.9710.754–1.2660.7961.3631.006–1.8160.0412.4230.860 5.5370.046
PNFA combined0.7060.515–0.9470.0190.960.808–1.1450.6281.2981.056–1.5840.0112.3671.211–4.2618.52 × 10−3**
SD discovery0.8080.554–1.1430.2580.7640.625–0.9400.0101.6491.300–2.0734.33 × 10−5*2.6141.152–5.2180.011
SD replication0.4020.198–0.7309.47 × 10−4*1.2170.906–1.6640.2081.1120.777–1.5520.5371.8780.497–5.0310.176
SD combined0.6510.470–0.8803.64 × 10−3*0.8980.759–1.0660.2051.4381.181–1.7412.93 × 10−4*2.3331.194–4.1999.08 × 10−3**
FTD–MND discovery0.5220.289–0.8769.30 × 10−31.0420.793–1.3880.8381.3110.948–1.780.0871.3240.266–4.0280.500
FTD–MND replication0.6740.303–1.3130.3111.1340.754–1.7640.6220.9340.53–1.5440.9010.0000.000–3.5360.629
FTD–MND combined0.5710.361–0.8626.01 × 10−31.0700.852–1.3570.6121.1880.901–1.5440.2020.8950.181–2.7141

SD = semantic dementia.

*P-values less than significant threshold 4.17 × 10−3 correcting multiple tests for association of three epsilon alleles with one of the four subtypes of FTD.

**P-values less than significant threshold 1.25 × 10−2 correcting multiple tests for association of homozygous ε4/ε4 genotypes with one of the four subtypes of FTD.

The summary statistics of variants used for deriving gene-based P-values for TOMM40 and APOE are given in Supplementary Tables 1A and B, respectively, and the regional plots are shown in Fig. 1A and B. The regional plots show many variants in TOMM40 with P-values < 0.05 that are in negligible LD (r2 < 0.2) with rs769449, a proxy of epsilon variant rs429358.

Regional plots at the TOMM40 /APOE locus in bvFTD cohorts. (A) Discovery; and (B) replication.
Figure 1

Regional plots at the TOMM40 /APOE locus in bvFTD cohorts. (A) Discovery; and (B) replication.

Furthermore we did not find significant gene-based association of RAB38 gene, which was identified as associated with bvFTD using the SNP-based GWAS. The association P-value of RAB38 gene in our joint-SNP gene-based analysis might be diluted due to inclusion of many non-risk variants, refer to Supplementary Fig. 3A and B for regional association plots at the RAB38 gene showing many non-risk variants within the gene’s transcription site.

Progressive non-fluent aphasia

The joint-cohort (discovery and replication) analysis revealed association for ARHGAP35 (P = 2.95 × 10−7) and SERPINA1 (P = 3.02 × 10−7) with PNFA (Table 1 and Supplementary Table 2). The regional plots for ARHGAP35 and SERPINA1 (Fig. 2A and B, respectively) show a robust LD block only for ARHGAP35 for which all variants show association P-values < 0.05 with PNFA (also refer to Supplementary Table 2A). In SERPINA1 many LD independent variants with PNFA association P-values < 0.05 can be observed.

Regional plots at (A) ARHGAP35 and (B) SERPINA1 loci in PNFA cohorts. (1) and (2) are discovery and replication cohorts, respectively.
Figure 2

Regional plots at (A) ARHGAP35 and (B) SERPINA1 loci in PNFA cohorts. (1) and (2) are discovery and replication cohorts, respectively.

Semantic dementia

No gene exceeded the gene-wide significance threshold 7.39 × 10−7, possibly because of the smaller sample size for this subtype, thus reduced power. The top gene identified in the combined analysis was WDR66 (P = 9.50 × 10−6, Table 1).

Frontotemporal dementia-motor neuron disease

No gene reached gene-wide significant association in the FTD-MND subtype. However, the top genes for FTD-MND were C9orf72 and IFNK with gene-based association P-value in joint-cohort analysis 1.23 × 10−6 and 1.77 × 10−6, respectively. Neither gene showed associations with any other subtypes of FTD (Table 1).

Frontotemporal dementia meta-analysis

The meta-analysis across all subtypes (bvFTD, semantic dementia, PNFA and FTD-MND) identified association of TOMM40 and APOE. It is worth noting that the bvFTD samples make nearly two-thirds of the total sample; hence, P-values for association with bvFTD dominated the meta-analysis of FTD subtypes.

Risk of APOE alleles on frontotemporal dementia subtypes

Based on the gene-based association results with TOMM40 and APOE we extended our analysis to the epsilon alleles and genotypes. We compared each FTD case cohort (discovery, replication and combined) against a total of 9390 ancestry-matched controls using Fisher’s exact test. For replication cohorts, we used rs769449 as a proxy for rs429358. The distribution of epsilon alleles and genotypes in our cohort is given in Supplementary Tables 3 and 4, respectively. We established the significance threshold for allele associations as 4.17 × 10−3 (0.05/12) correcting for three epsilon alleles and four FTD subtypes. We identified that the ɛ2 allele significantly reduces the risk of bvFTD [odds ratio (OR) = 0.772, P = 3.88 × 10−4] and semantic dementia (OR = 0.651, P = 3.64 × 10−3). We observed marginal association (P < 0.05) of ɛ2 allele with PNFA (OR = 0.706, P = 0.019) and moderate with FTD-MND (OR = 0.571, P = 6.01 × 10−3). The ɛ4 allele significantly increased risk of bvFTD (OR = 1.278, P = 8.14 × 10−6) and semantic dementia (OR = 1.438, P = 2.93 × 10−4). The association for the disease increasing effect of ɛ4 allele was marginal (P < 0.05) for PNFA (OR = 1.298, P = 0.011), and the result was inconclusive for FTD-MND (OR = 1.188, P = 0.202) possibly due to underpowered sample size.

We also quantified the risk of homozygous ɛ4/ɛ4 genotype on FTD subtypes. We used 1.25 × 10−2 as a significance threshold for association testing of four subtypes with homozygous ɛ4/ɛ4 genotype. The homozygous ɛ4/ɛ4 genotype showed significant association with increased risk for bvFTD (OR = 1.627, P = 0.012), PNFA (OR = 2.367, P = 8.52 × 10−3) and semantic dementia (OR = 2.333, P = 9.08 × 10−3) with notable P- and OR-values for PNFA and semantic dementia compared to the effect size of a single copy of ɛ4 allele for respective FTD subtypes. We did not perform association between homozygous ɛ2/ɛ2 genotypes with FTD subtypes due to its low frequency in our cohort.

Functional characterization of associated genes

We extracted the gene expression profiles of APOE, TOMM40, ARHGAP35 and SERPINA1 across different human brain tissues from the GTeX database (Mele et al., 2015) (see Supplementary Fig. 4A–D for respective genes). The APOE, TOMM40 and ARHGAP35 genes are strongly expressed in different brain tissues. Notably the anterior cingulate cortex (Brodmann area 24) and the frontal cortex (Brodmann area 9) are the top tissues for ARHGAP35 gene expression. The anterior cingulate cortex is one of the early affected regions in FTD patients (Seeley et al., 2006, 2008); this area is reported to be involved in language control and resolving non-verbal conflict (Abutalebi et al., 2012). The SERPINA1 gene did not show strong expression in the brain tissues.

We used the HaploReg (version 4.1) software (Ward and Kellis, 2012) to investigate functionally annotated variants linked with variants used in deriving gene-based P-values of TOMM40 (those in negligible LD r2 < 0.2 with epsilon variants), ARHGAP35, and SERPINA1, respectively. We found that all SNPs used in deriving gene-based P-values of TOMM40, ARHGAP35 and SERPINA1 are in strong LD (r2 > 0.8, in 1000 Genomes phase 1 European panel) with at least one variant residing in the regulatory regions such as chromatin marks or DNase hypersensitive sites, suggesting a possible regulatory role (see Supplementary Tables 5A–C for HaploReg results for variants in TOMM40, ARHGAP35 and SERPINA1, respectively). Overall we identified 21, 56 and 93 regulatory variants in LD with SNPs driving gene-based P-values of TOMM40, ARHGAP35 and SERPINA1 genes, respectively. We further ranked these regulatory variants based on their functional relevance using the RegulomeDB (version 1.1) software (Boyle et al., 2012) (see Supplementary Tables 6A–C for detailed RegulomeDB results of these variants mapped to TOMM40, ARHGAP35 and SERPINA1 genes, respectively).

Discussion

Here we report novel genetic insight into FTD and its clinical subtypes using a joint-SNP gene-based approach. We identified association of the TOMM40 and APOE genes with bvFTD, and the ARHGAP35 and SERPINA1 genes with PNFA.

Our study suggested TOMM40 as the top gene in bvFTD. The TOMM40 gene encodes a channel forming subunit of the translocase of the mitochondrial outer membrane (TOM complex), which facilitates translocation of unfolded proteins from the cytosol into the mitochondrial intermembrane space for use in oxidative phosphorylation (Mager et al., 2011). Recently, Bannwarth et al. (2014) reported mitochondrial origin in pathogenesis of FTD-ALS diseases through association of variants in CHCHD10. There is growing evidence suggesting a role of mitochondria in neurodegenerative disorders, also including Parkinson’s disease (Parker et al., 1989; Schapira et al., 1990), Huntington’s disease (Moran et al., 2012) and Alzheimer’s disease (Petrozzi et al., 2007; Hroudova et al., 2014). Our conditional analysis suggests the TOMM40 association with bvFTD being independent from the epsilon alleles; however, interdependence with the epsilon alleles cannot be fully excluded at this stage. Future, more sophisticated, sequencing studies might further confirm the independence we highlight here.

The association of the APOE gene with bvFTD was primarily driven by SNPs rs7412 and rs429358 (or variants in strong LD such as rs769449) suggesting that, if on one hand the TOMM40 association is independent from the epsilon alleles, on the other the APOE association (which is, however, lower comparatively to that of TOMM40) is dependent on the epsilon alleles lending support for the inference to be made that two different haplotypes at this locus might confer risk differently. The SNPs rs7412 and rs429358 determine the APOE epsilon alleles ɛ2, ɛ3 and ɛ4. We quantified the risk of epsilon alleles across clinical subtypes of FTD diagnosed using the Neary’s criteria against neurologically normal controls, and saw that the ɛ2 and ɛ4 alleles showed protective and increased disease risk effects, respectively, for FTD subtypes (strong associations for bvFTD and semantic dementia, and marginal associations for PNFA and FTD-MND). Interestingly, individuals carrying homozygous copies of the ɛ4 allele revealed higher risk for PNFA and semantic dementia (OR > 2.3) suggesting dose-dependent effect for each copy of a gene. The pattern of association of epsilon alleles with FTD subtypes might reflect the potential overlap between patients diagnosed with clinical FTD and Alzheimer’s disease (van der Zee et al., 2008; Bang et al., 2015) or a genuine association with FTD and its subtypes given the increasing number of studies arguing in favour of the latter hypothesis (Engelborghs et al., 2006; Agosta et al., 2009; Seripa et al., 2011; Rubino et al., 2013; Ferrari et al., 2015).

In the CNS, APOE is synthesized in response to neuronal injury or stress to initiate the neuronal repair mechanisms. The ɛ4 carriers are hypothesized to have reduced neuronal repair capacity compared to the other alleles (Rubino et al., 2013). The protein products of APOE were also reported to modulate neuroinflammation (Tai et al., 2015). The hypothesis of enhanced inflammatory response in FTD patients is supported by both neuroimaging and genetic studies (Cagnin et al., 2004; Rainero et al., 2009). Tau pathology is found in up to ∼50% of FTLD cases (Halliday et al., 2012); interestingly the knock-in study in mice showed association between epsilon alleles and the concentration of hyper-phosphorylated tau in neurons: ɛ4 knock in mice showed higher concentration of hyper-phosphorylated tau than ɛ3 knock in mice (Inbar et al., 2010). It is worth noting that in this scenario the APOE locus and the epsilon allelic variability might impact processes such as modulation of neuronal repair mechanisms, neuroinflammation, broad lipid metabolism, synaptic plasticity, neuronal toxicity and tau phosphorylation (Verghese et al., 2011). It is likely that variability in the genes or isoforms turnover or larger haplotype blocks at this locus, coupled with ageing, might influence negative outcomes in brain and thus support our findings from a biological and functional perspective. More work should be directed towards testing these possibilities in the future.

Our study is the first to report association of the ARHGAP35 and SERPINA1 genes with PNFA. The ARHGAP35 gene encodes the glucocorticoid receptor DNA-binding factor 1, which is a repressor of glucocorticoid receptor (hGR) transcription. At the cellular level, the glucocorticoid receptor mediates the maintenance of basal and stress-related homeostasis. The second gene we found associated with PNFA was SERPINA1, which was previously reported to be associated with cortisol level (Bolton et al., 2014) (top variants: rs11621961, rs12589136, rs2749527) and serum lipid profile (Inouye et al., 2012) (top variant: rs1303). The non-synonymous variant rs1303 (Glu400Asp) in SERPINA1, which is also in moderate LD with morning plasma cortisol level associated variant rs12589136 (Bolton et al., 2014), showed PNFA association P-values < 0.05 in both discovery and replication cohorts (Supplementary Table 2B). The top SERPINA1 variant in the PNFA discovery cohort rs11628917 (this variant in moderate LD with rs1303) is an established blood eQTL (Westra et al., 2013), with the C allele increasing SERPINA1 expression in blood (Inouye et al., 2012; Westra et al., 2013). We observed that the C allele at rs11628917 increased the risk of PNFA in both discovery (OR = 2.893, P = 2.58 × 10−6) and replication (OR = 1.385, P = 0.077) cohorts. The SERPINA1 gene encodes protease inhibitor α1-antitrypsin enzyme, which inhibits cleavage of the reactive centre loop of the corticosteroid binding globulin (CBG) by ceasing neutrophil elastase activity (Lewis and Elder, 2014). The reactive centre loop cleavage by neutrophil elastase reduces the CBG binding affinity to cortisol. During stress CBG activity is positively correlated with the glucocorticoid access to the brain (Moisan et al., 2014). Increased glucocorticoid level in brain can activate glucocorticoid signalling through binding to the low affinity glucocorticoid receptors and result in reduced neurogenesis and impaired neuroplasticity (Anacker et al., 2011). This hypothesis suggesting the role of enhanced glucocorticoid signalling leading to neurodegeneration in PNFA patients is based on our current preliminary report on association of ARHGAP35 and SERPINA1 with PNFA; this finding will need to be further explored and replicated in an independent cohort.

In conclusion, we report genetic associations for FTD and its subtypes—notably of the TOMM40 and APOE genes with bvFTD and the ARHGAP35 and SERPINA1 genes with PNFA—using the joint-SNP gene-based approach. This approach improves power of the association test by combining signals across variants in a functional unit such as a gene. Replication and functional characterization of these findings will further establish their role in pathology of FTD and help towards a better management of the disease.

Abbreviations

    Abbreviations
     
  • bvFTD

    behavioural variant of frontotemporal dementia

  •  
  • GWAS

    genome-wide association study

  •  
  • FTD

    frontotemporal dementia

  •  
  • MND

    motor neuron disease

  •  
  • LD

    linkage disequilibrium

  •  
  • MAGMA

    Multi-marker Analysis of GenoMic Annotation

  •  
  • PNFA

    progressive non-fluent aphasia

  •  
  • SNP

    single nucleotide polymorphism

Acknowledgements

We acknowledge the investigators of the original GWAS by IFGC, (Appendix 1 and Supplementary material). We also acknowledge contributors to the Genotype-Tissue Expression (GTEx) Project (Supplementary material).

Funding

The Netherlands Organization for Scientific Research (NWO VICI 453-14-005) funded this work. R.F. is supported by the Alzheimer’s Society (grant number 284). P.H. and J.H. are funded by the JPND project RiMod-FTD.

Supplementary material

Supplementary material is available at Brain online.

Appendix 1

R. Ferrari; D. G. Hernandez; M. A. Nalls; J. D. Rohrer; A. Ramasamy; J. B. J. Kwok; C. Dobson-Stone; P. R. Schofield; G. M. Halliday; J. R. Hodges; O. Piguet; L. Bartley; E. Thompson; E. Haan; I. Hernández; A. Ruiz; M. Boada; B. Borroni; A. Padovani; C. Cruchaga; N. J. Cairns; L. Benussi; G. Binetti; R. Ghidoni; G. Forloni; D. Albani; D. Galimberti; C. Fenoglio; M. Serpente; E. Scarpini; J. Clarimón; A. Lleó; R. Blesa; M. Landqvist Waldö; K. Nilsson; C. Nilsson; I. R. A. Mackenzie; G-Y. R. Hsiung; D. M. A. Mann; J. Grafman; C. M. Morris; J. Attems; T. D. Griffiths; I. G. McKeith; A. J. Thomas; P. Pietrini; E. D. Huey; E. M. Wassermann; A. Baborie; E. Jaros; M. C. Tierney; P. Pastor; C. Razquin; S. Ortega-Cubero; E. Alonso; R. Perneczky; J. Diehl-Schmid; P. Alexopoulos; A. Kurz; I. Rainero; E. Rubino; L. Pinessi; E. Rogaeva; P. St George-Hyslop; G. Rossi; F. Tagliavini; G. Giaccone; J. B. Rowe; J. C. M. Schlachetzki; J. Uphill; J. Collinge; S. Mead; A. Danek; V. M. Van Deerlin; M. Grossman; J. Q .Trojanowski; J. van der Zee; M. Cruts; C. Van Broeckhoven; S. F. Cappa; I. Leber; D. Hannequin; V. Golfier; M. Vercelletto; A. Brice; B. Nacmias; S. Sorbi; S. Bagnoli; I. Piaceri; J. E. Nielsen; L. E. Hjermind; M. Riemenschneider; M. Mayhaus; B. Ibach; G .Gasparoni; S. Pichler; W. Gu; M. N. Rossor; N. C. Fox; J. D. Warren; M. G. Spillantini; H. R. Morris; P. Rizzu; P. Heutink; J. S. Snowden; S. Rollinson; A. Richardson; A. Gerhard; A. C. Bruni; R. Maletta; F. Frangipane; C. Cupidi; L. Bernardi; M. Anfossi; M. Gallo; M. E. Conidi; N. Smirne; R. Rademakers; M. Baker; D. W. Dickson; N. R. Graff-Radford; R. C. Petersen; D. Knopman; K. A. Josephs; B. F. Boeve; J. E. Parisi; W. W. Seeley; B. L. Miller; A. M. Karydas; H. Rosen; J. C. van Swieten; E. G. P. Dopper; H. Seelaar; Y. A. L. Pijnenburg; P. Scheltens; G. Logroscino; R. Capozzo; V. Novelli; A. A. Puca; M. Franceschi; A. Postiglione; G. Milan; P. Sorrentino; M. Kristiansen; H-H. Chiang; C. Graff; F. Pasquier; A. Rollin; V. Deramecourt; T. Lebouvier; D. Kapogiannis; L. Ferrucci; S. Pickering-Brown; A. B. Singleton; J. Hardy; P. Momeni.

References

Abutalebi
J
,
Della Rosa
PA
,
Green
DW
,
Hernandez
M
,
Scifo
P
,
Keim
R
et al.
Bilingualism tunes the anterior cingulate cortex for conflict monitoring
.
Cereb Cortex
2012
;
22
:
2076
86
.

Agosta
F
,
Vossel
KA
,
Miller
BL
,
Migliaccio
R
,
Bonasera
SJ
,
Filippi
M
et al.
Apolipoprotein E epsilon4 is associated with disease-specific effects on brain atrophy in Alzheimer’s disease and frontotemporal dementia
.
Proc Natl Acad Sci USA
2009
;
106
:
2018
22
.

Anacker
C
,
Zunszain
PA
,
Carvalho
LA
,
Pariante
CM
.
The glucocorticoid receptor: pivot of depression and of antidepressant treatment?
.
Psychoneuroendocrinology
2011
;
36
:
415
25
.

Bang
J
,
Spina
S
,
Miller
BL
.
Frontotemporal dementia
.
Lancet
2015
;
386
:
1672
82
.

Bannwarth
S
,
Ait-El-Mkadem
S
,
Chaussenot
A
,
Genin
EC
,
Lacas-Gervais
S
,
Fragaki
K
et al.
A mitochondrial origin for frontotemporal dementia and amyotrophic lateral sclerosis through CHCHD10 involvement
.
Brain
2014
;
137
(
Pt 8
):
2329
45
.

Bolton
JL
,
Hayward
C
,
Direk
N
,
Lewis
JG
,
Hammond
GL
,
Hill
LA
et al.
Genome wide association identifies common variants at the SERPINA6/SERPINA1 locus influencing plasma cortisol and corticosteroid binding globulin
.
PLoS Genet
2014
;
10
:
e1004474
.

Boyle
AP
,
Hong
EL
,
Hariharan
M
,
Cheng
Y
,
Schaub
MA
,
Kasowski
M
et al.
Annotation of functional variation in personal genomes using RegulomeDB
.
Genome Res
2012
;
22
:
1790
7
.

Cagnin
A
,
Rossor
M
,
Sampson
EL
,
Mackinnon
T
,
Banati
RB
.
In vivo detection of microglial activation in frontotemporal dementia
.
Ann Neurol
2004
;
56
:
894
7
.

Cairns
NJ
,
Bigio
EH
,
Mackenzie
IR
,
Neumann
M
,
Lee
VM
,
Hatanpaa
KJ
et al.
Neuropathologic diagnostic and nosologic criteria for frontotemporal lobar degeneration: consensus of the consortium for frontotemporal lobar degeneration
.
Acta Neuropathol
2007
;
114
:
5
22
.

de Leeuw
CA
,
Mooij
JM
,
Heskes
T
,
Posthuma
D
.
MAGMA: generalized gene-set analysis of GWAS data
.
PLoS Comput Biol
2015
;
11
:
e1004219
.

Engelborghs
S
,
Dermaut
B
,
Marien
P
,
Symons
A
,
Vloeberghs
E
,
Maertens
K
et al.
Dose dependent effect of APOE epsilon4 on behavioral symptoms in frontal lobe dementia
.
Neurobiol Aging
2006
;
27
:
285
92
.

Ferrari
R
,
Grassi
M
,
Salvi
E
,
Borroni
B
,
Palluzzi
F
,
Pepe
D
et al.
A genome-wide screening and SNPs-to-genes approach to identify novel genetic risk factors associated with frontotemporal dementia
.
Neurobiol Aging
2015
;
36
:
2904.e13
26
.

Ferrari
R
,
Hernandez
DG
,
Nalls
MA
,
Rohrer
JD
,
Ramasamy
A
,
Kwok
JB
et al.
Frontotemporal dementia and its subtypes: a genome-wide association study
.
Lancet Neurol
2014
;
13
:
686
99
.

Ferrari
R
,
Thumma
A
,
Momeni
P
.
Molecular genetics of frontotemporal dementia
.
eLS
:
John Wiley & Sons Ltd
,
2013
. doi: 10.1002/9780470015902.a0024457

Genomes Project Consortium
,
Abecasis
GR
,
Auton
A
,
Brooks
LD
,
DePristo
MA
,
Durbin
RM
et al.
An integrated map of genetic variation from 1,092 human genomes
.
Nature
2012
;
491
:
56
65
.

Gorno-Tempini
ML
,
Hillis
AE
,
Weintraub
S
,
Kertesz
A
,
Mendez
M
,
Cappa
SF
et al.
Classification of primary progressive aphasia and its variants
.
Neurology
2011
;
76
:
1006
14
.

Hagg
S
,
Ganna
A
,
Van Der Laan
SW
,
Esko
T
,
Pers
TH
,
Locke
AE
et al.
Gene-based meta-analysis of genome-wide association studies implicates new loci involved in obesity
.
Hum Mol Genet
2015
;
24
:
6849
60
.

Halliday
G
,
Bigio
EH
,
Cairns
NJ
,
Neumann
M
,
Mackenzie
IR
,
Mann
DM
.
Mechanisms of disease in frontotemporal lobar degeneration: gain of function versus loss of function effects
.
Acta Neuropathol
2012
;
124
:
373
82
.

Hroudova
J
,
Singh
N
,
Fisar
Z
.
Mitochondrial dysfunctions in neurodegenerative diseases: relevance to Alzheimer’s disease
.
Biomed Res Int
2014
;
2014
:
175062
.

Inbar
D
,
Belinson
H
,
Rosenman
H
,
Michaelson
DM
.
Possible role of tau in mediating pathological effects of apoE4 in vivo prior to and following activation of the amyloid cascade
.
Neurodegener Dis
2010
;
7
:
16
23
.

Inouye
M
,
Ripatti
S
,
Kettunen
J
,
Lyytikainen
LP
,
Oksala
N
,
Laurila
PP
et al.
Novel loci for metabolic networks and multi-tissue expression studies reveal genes for atherosclerosis
.
PLoS Genet
2012
;
8
:
e1002907
.

Lewis
JG
,
Elder
PA
.
The reactive centre loop of corticosteroid-binding globulin (CBG) is a protease target for cortisol release
.
Mol Cell Endocrinol
2014
;
384
:
96
101
.

Liu
JZ
,
McRae
AF
,
Nyholt
DR
,
Medland
SE
,
Wray
NR
,
Brown
KM
et al.
A versatile gene-based test for genome-wide association studies
.
Am J Hum Genet
2010
;
87
:
139
45
.

Mager
F
,
Gessmann
D
,
Nussberger
S
,
Zeth
K
.
Functional refolding and characterization of two Tom40 isoforms from human mitochondria
.
J Membr Biol
2011
;
242
:
11
21
.

Mele
M
,
Ferreira
PG
,
Reverter
F
,
DeLuca
DS
,
Monlong
J
,
Sammeth
M
et al.
Human genomics
.
The human transcriptome across tissues and individuals. Science
2015
;
348
:
660
5
.

Moisan
MP
,
Minni
AM
,
Dominguez
G
,
Helbling
JC
,
Foury
A
,
Henkous
N
et al.
Role of corticosteroid binding globulin in the fast actions of glucocorticoids on the brain
.
Steroids
2014
;
81
:
109
15
.

Moran
M
,
Moreno-Lastres
D
,
Marin-Buera
L
,
Arenas
J
,
Martin
MA
,
Ugalde
C
.
Mitochondrial respiratory chain dysfunction: implications in neurodegeneration
.
Free Radic Biol Med
2012
;
53
:
595
609
.

Neary
D
,
Snowden
JS
,
Gustafson
L
,
Passant
U
,
Stuss
D
,
Black
S
et al.
Frontotemporal lobar degeneration: a consensus on clinical diagnostic criteria
.
Neurology
1998
;
51
:
1546
54
.

Parker
WD
Jr,
Boyson
SJ
,
Parks
JK
.
Abnormalities of the electron transport chain in idiopathic Parkinson’s disease
.
Ann Neurol
1989
;
26
:
719
23
.

Petrozzi
L
,
Ricci
G
,
Giglioli
NJ
,
Siciliano
G
,
Mancuso
M
.
Mitochondria and neurodegeneration
.
Biosci Rep
2007
;
27
:
87
104
.

Rabinovici
GD
,
Miller
BL
.
Frontotemporal lobar degeneration: epidemiology, pathophysiology, diagnosis and management
.
CNS Drugs
2010
;
24
:
375
98
.

Rainero
I
,
Rubino
E
,
Cappa
G
,
Rota
E
,
Valfre
W
,
Ferrero
P
et al.
Pro-inflammatory cytokine genes influence the clinical features of frontotemporal lobar degeneration
.
Dement Geriatr Cogn Disord
2009
;
27
:
543
7
.

Rascovsky
K
,
Hodges
JR
,
Knopman
D
,
Mendez
MF
,
Kramer
JH
,
Neuhaus
J
et al.
Sensitivity of revised diagnostic criteria for the behavioural variant of frontotemporal dementia
.
Brain
2011
;
134
(
Pt 9
):
2456
77
.

Rohrer
JD
,
Guerreiro
R
,
Vandrovcova
J
,
Uphill
J
,
Reiman
D
,
Beck
J
et al.
The heritability and genetics of frontotemporal lobar degeneration
.
Neurology
2009
;
73
:
1451
6
.

Rubino
E
,
Vacca
A
,
Govone
F
,
De Martino
P
,
Pinessi
L
,
Rainero
I
.
Apolipoprotein E polymorphisms in frontotemporal lobar degeneration: a meta-analysis
.
Alzheimers Dement
2013
;
9
:
706
13
.

Schapira
AH
,
Cooper
JM
,
Dexter
D
,
Clark
JB
,
Jenner
P
,
Marsden
CD
.
Mitochondrial complex I deficiency in Parkinson’s disease
.
J Neurochem
1990
;
54
:
823
7
.

Seelaar
H
,
Rohrer
JD
,
Pijnenburg
YA
,
Fox
NC
,
van Swieten
JC
.
Clinical, genetic and pathological heterogeneity of frontotemporal dementia: a review
.
J Neurol Neurosurg Psychiatry
2011
;
82
:
476
86
.

Seeley
WW
,
Carlin
DA
,
Allman
JM
,
Macedo
MN
,
Bush
C
,
Miller
BL
et al.
Early frontotemporal dementia targets neurons unique to apes and humans
.
Ann Neurol
2006
;
60
:
660
7
.

Seeley
WW
,
Crawford
R
,
Rascovsky
K
,
Kramer
JH
,
Weiner
M
,
Miller
BL
et al.
Frontal paralimbic network atrophy in very mild behavioral variant frontotemporal dementia
.
Arch Neurol
2008
;
65
:
249
55
.

Seripa
D
,
Bizzarro
A
,
Panza
F
,
Acciarri
A
,
Pellegrini
F
,
Pilotto
A
et al.
The APOE gene locus in frontotemporal dementia and primary progressive aphasia
.
Arch Neurol
2011
;
68
:
622
8
.

Tada
H
,
Won
HH
,
Melander
O
,
Yang
J
,
Peloso
GM
,
Kathiresan
S
.
Multiple associated variants increase the heritability explained for plasma lipids and coronary artery disease
.
Circ Cardiovasc Genet
2014
;
7
:
583
7
.

Tai
LM
,
Ghura
S
,
Koster
KP
,
Liakaite
V
,
Maienschein-Cline
M
,
Kanabar
P
et al.
APOE-modulated Abeta-induced neuroinflammation in Alzheimer’s disease: current landscape, novel data, and future perspective
.
J Neurochem
2015
;
133
:
465
88
.

Van Deerlin
VM
,
Sleiman
PM
,
Martinez-Lage
M
,
Chen-Plotkin
A
,
Wang
LS
,
Graff-Radford
NR
et al.
Common variants at 7p21 are associated with frontotemporal lobar degeneration with TDP-43 inclusions
.
Nat Genet
2010
;
42
:
234
9
.

van der Zee
J
,
Sleegers
K
,
Van Broeckhoven
C
.
Invited article: the Alzheimer disease-frontotemporal lobar degeneration spectrum
.
Neurology
2008
;
71
:
1191
7
.

Verghese
PB
,
Castellano
JM
,
Holtzman
DM
.
Apolipoprotein E in Alzheimer’s disease and other neurological disorders
.
Lancet Neurol
2011
;
10
:
241
52
.

Ward
LD
,
Kellis
M
.
HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants
.
Nucleic Acids Res
2012
;
40
(
Database issue
):
D930
4
.

Westra
HJ
,
Peters
MJ
,
Esko
T
,
Yaghootkar
H
,
Schurmann
C
,
Kettunen
J
et al.
Systematic identification of trans eQTLs as putative drivers of known disease associations
.
Nat Genet
2013
;
45
:
1238
43
.

Woollacott
IO
,
Rohrer
JD
.
The clinical spectrum of sporadic and familial forms of frontotemporal dementia
.
J Neurochem
2016
;
138
(Suppl 1):
6
31
.

Yang
J
,
Ferreira
T
,
Morris
AP
,
Medland
SE
Genetic Investigation of ATC, Replication DIG
et al.
Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits
.
Nat Genet
2012
;
44
:
369
75, S1–3
.

Zhang
J
,
Zhang
Y
,
Yang
J
,
Zhang
L
,
Sun
L
,
Pan
HF
et al.
Three SNPs in chromosome 11q23.3 are independently associated with systemic lupus erythematosus in Asians
.
Hum Mol Genet
2014
;
23
:
524
33
.

Supplementary data