Exome Sequencing of Extended Families with Alzheimer's Disease Identifies Novel Genes Implicated in Cell Immunity and Neuronal Function

¹1John P. Hussman Institute for Human Genomics, University of Miami Miller School of Medicine, Miami, FL, USA

²2Department of Neurology, University of Miami Miller School of Medicine, Miami, FL, USA

³3John T. Macdonald Foundation, Department of Human Genetics, University of Miami Miller School of Medicine, Miami, FL, USA

⁴4Mental Health and Behavioral Sciences Service, Miami Veterans Affairs, Miami, FL, USA

⁵Departments of Medicine, Neurology, Ophthalmology, Genetics and Genomics, Epidemiology and Biostatistics, Boston University, Boston, MA, USA

⁶Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, OH, USA

^†Authors contributed equally

*Corresponding Author:

Pericak VMA
John T. Macdonald Foundation
Department of Human Genetics
University of Miami Miller School of Medicine
Miami, FL, USA
Tel: 3052432308
Fax: 3052432396
E-mail: mpericak@med.miami.edu

Received date: June 16, 2017; Accepted date: July 24, 2017; Published date: July 31, 2017

Citation: Holly NC, Brian WK, Hamilton NKL, Sophie R, Martin AK, et al. (2017) Exome Sequencing of Extended Families with Alzheimer’s Disease Identifies Novel Genes Implicated in Cell Immunity and Neuronal Function. J Alzheimers Dis Parkinsonism 7:355. doi:10.4172/2161-0460.1000355

Copyright: © 2017 Holly NC, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Visit for more related articles at Journal of Alzheimers Disease & Parkinsonism

Abstract

Objective: Alzheimer’s disease (AD) is a neurodegenerative disorder for which more than 20 genetic loci have been implicated to date. However, studies demonstrate not all genetic factors have been identified. Therefore, in this study we seek to identify additional rare variants and novel genes potentially contributing to AD. Methods: Whole exome sequencing was performed on 23 multi-generational families with an average of eight affected subjects. Exome sequencing was filtered for rare, nonsynonymous and loss-of-function variants. Alterations predicted to have a functional consequence and located within either a previously reported AD gene, a linkage peak (LOD>2), or clustering in the same gene across multiple families, were prioritized. Results: Rare variants were found in known AD risk genes including AKAP9, CD33, CR1, EPHA1, INPP5D, NME8, PSEN1, SORL1, TREM2 and UNC5C. Three families had five variants of interest in linkage regions with LOD>2. Genes with segregating alterations in these peaks include CD163L1 and CLECL1, two genes that have both been implicated in immunity, CTNNA1, which encodes a catenin in the cerebral cortex and MIEF1, a gene that may induce mitochondrial dysfunction and has the potential to damage neurons. Four genes were identified with alterations in more than one family include PLEKHG5, a gene that causes Charcot-Marie-Tooth disease and THBS2, which promotes synaptogenesis. Conclusion: Utilizing large families with a heavy burden of disease allowed for the identification of rare variants co-segregating with disease. Variants were identified in both known AD risk genes and in novel genes.

Keywords

Alzheimer’s disease; Dominant inheritance; Linkage; Multiplex; Whole exome sequencing

Introduction

Alzheimer’s disease (AD) is the leading cause of dementia in the elderly [1]. The majority of individuals present with late-onset AD (≥ 65 years), but early-onset (<65 years) has also been reported in ~5% of cases. Both common genetic variants, such as the APOE ε4 allele, and rare variants, have been found to impact risk for both early- and late- onset AD [2-5]. While more than 20 genetic loci have been connected with late-onset AD to date, the underlying genetic architecture is complex and new risk genes are still being identified [6].

While genome-wide association studies (GWAS) have been key in identifying a majority of the novel regions of genetic risk in the past ten years, by design, GWAS are unlikely to recognize risk variants with rare frequencies in the population and necessitate the use of large cohorts of hundreds or even thousands of individuals to reach statistically significant conclusions [6]. In contrast, whole exome sequencing (WES) provides an alternative and complementary method to locate rare alterations in genes which may have medium to large effects on disease risk and require far fewer participants [6-8]. WES studies have identified new mutations in both known AD genes and novel risk genes, including AKAP9, PLD3, TREM2 and UNC5C, as well as protective variants, such as those in TREML2 [7-17]. Moreover, studying families with a heavy burden of AD and searching for genetic changes that segregate with disease can provide a unique opportunity to locate rare variants in novel risk genes such as NOTCH3, PLD3 and TTC3 [9,13,18].

These large AD families can reveal how multiple genetic variants may act in concert to influence risk [19-21]. For example, the APOE ε2 allele was found to delay the age of onset by ~12 years in carriers of the E280A mutation in the PSEN1 gene in the early-onset ‘Paisa’ pedigree [19]. In addition, genetic linkage can assist in narrowing genomic regions of interest potentially related to disease in large families [22]. In an effort to discover novel genes that may contribute to late-onset AD risk, we performed WES in 23 multiplex families that present with dominant inheritance patterns and prioritized variants that were inherited from common ancestors.

Materials and Methods

Patient ascertainment of extended AD families

240 individuals (77 AD subjects, 4 individuals with mild cognitive impairment (MCI) and 159 unaffected relatives) from 23 families of European ancestry heavily affected with late-onset AD were utilized in this study (Supplementary Table 1). All family members were recruited after providing informed consent and with approval by the relevant institutional review boards. Affected individuals meet the standard NINCDS-ADRDA criteria for AD and MCI [23-25]. In addition, cognitive and neuropsychiatric data were collected on all affected indivduals using the NCRAC LOAD battery, the Geriatric Depression Scale (GDS15), the Cornell Scale for Depression in Dementia (CSDD) and the Neuropsychiatric Inventory Questionnaire (NPIQ).

Whole exome sequencing and variant detection

99 individuals (77 AD patients, 4 individuals with MCI, and 18 unaffected relatives) from 23 AD extended families underwent WES (Supplementary Table 1). Three micrograms of DNA from each sample were prepared using the SureSelect Human All Exon 50Mb Kit (Agilent Technologies) and the Paired-End Multiplexed Sequencing library kit (Illumina). Exome capture and sequence library construction was performed on a Sciclone G3 NGS Workstation (Caliper Life Sciences) and DNA was tested for uniform enrichment of targets with qPCR following established protocols provided by Agilent. Two exome sample libraries were sequenced per lane on a HiSeq 2000 Sequencing System (Illumina) in paired-end 2 × 100 bp runs. Sequencing data was processed using the Illumina RTA base calling pipeline v1.8. Reads were aligned to the human reference genome (hg19) with the Burrows- Wheeler Aligner (BWA) and variant calling performed with the Genome Analysis Toolkit (GATK) version 2.8 [26,27]. GATK parameters for variant quality control included duplicate sequence read removal, minimum read depth of 5, genotype quality (GQ) ≥ 20, variant quality score recalibration (VQSR, VQSLOD>0) and Genome Mappability Scores equal to 1 for the 35 base pair (bp) track and greater than or equal to 0.5 for the 20 bp track from the Duke Uniqueness Track [28]. The Duke uniqueness scores, generated for the ENCODE project and available as tracks in the University of California, Santa Cruz (UCSC) Genome Browser, report how unique a sequence is, where scores of 1 represent a completely unique sequence, a score of 0.5 indicates the sequence occurs exactly twice, and 0 represents the sequence occurs >4 times in the genome [29,30]. Small insertions and deletions were recognized by aligning the data with Bowtie2 and analyzing with the Pindel program [31,32].

Genotyping and variant filtering

234 individuals, including all 99 samples that had WES , were evaluated by genome-wide SNP (single nucleotide polymorphism) arrays including the Human 1Mv1 BeadChip, the 1M-DuoV3 BeadChip, the HumanOmniExpress-12 v1.0 BeadChip, and the HumanOmni2.5- 4v1 BeadChip. All chips were processed using the Tecan EVO-1 robot and BeadChips were scanned with either the Illumina BeadArray Reader or iScan. Data was extracted by the Genome Studio software and a GenCall cutoff score of 0.15 was used. Samples were required to have a genotyping call rate of 98% or higher, and SNPs a call rate of 95% or greater, to pass quality control. SNPs were only included in the analysis if they were present in at least 60% of samples across all platforms. Checks for relatedness, Mendelian inconsistencies, gender based on X-chromosome heterozygosity, and concordance between the genotypes of the variants identified through exome sequencing and genotyping were evaluated with PLINK version 1.07 [33]. All samples passed the quality control metrics.

Genotyping information was further used to delineate identical by descent (IBD) regions within each multiplex AD family. IBD filtering was implemented through the extended haplotype procedure in MERLIN version 1.1.2 [34]. Regions shared across all available AD individuals within a family were used to determine the IBD sharing segments and were, therefore, unique within each family. To determine the start and stop positions of IBD sharing regions within each family, the MERLIN output was evaluated in a sliding window of ten SNPs, defining IBD as sharing at each location with a threshold >50%.

Linkage analysis

Nonparametric and parametric two-point and multipoint linkage analyses were performed using MERLIN [32]. A disease allele frequency of 0.0001 was used in an affecteds-only model for parametric analysis. PLINK was employed for LD pruning in the multipoint analysis, with CEU HapMap data as the reference population and the following settings: the indep-pairwise option with a window size of 50, a step of 5 and an r² threshold of 0.5 [33,34].

Variant annotation and prioritization

Alterations passing quality measurements were annotated with the KGGSeq and ANNOVAR programs [35,36]. Variants were normalized prior to annotation [37]. Ensembl, RefSeq, and Gencode transcripts were all annotated, and the top consequence per gene was used for prioritization. CADD v1.3 scores were downloaded from the CADD server (http://cadd.gs.washington.edu/home) [38]. Figure 1 is an overview of the filtering and prioritization strategies used in this study. Brief descriptions of our three prioritization strategies are described below.

Figure 1: Study design. Strategy for processing the samples and prioritizing the variants that were resulting from whole exome sequencing.

Table 1: Known AD genes and loci with rare, potentially damaging variants.

Variants in reported AD genes or loci: For all of the families, we evaluated whether variants were located in known AD risk genes; this includes genes identified in both early (APP, PSEN1, PSEN2, GRN and TREM2) and late-onset AD (Supplemental Table 2) [3,4]. Variants of interest were restricted to those with a minor allele frequency (MAF) ≤ 2% in the Kaviar Genomic Variant Database (version 160204-Public, 77,781 individuals) since these genes are known loci for AD [39]. The top variants of interest were validated by traditional Sanger sequencing.

Table 2: Families with segregating, rare, potentially damaging variants in high LOD regions.

Families with LOD scores >2: For each of the families, variants that segregated in all sequenced, affected individuals within areas LOD>2 were evaluated. Variants with a global MAF ≤ 1% in the Kaviar Database were prioritized. A MAF cutoff of ≤ 1% was implemented because variants with a MAF>1% in any ethnic population are unlikely to be a highly penetrant risk variant for AD [39-41]. This stricter MAF criteria was utilized to attempt to identify novel risk genes as opposed to variants in known AD genes. Variants were also prioritized based on their potential pathogenicity with the Combined Annotation- Dependent Depletion (CADD) score; scores ≥ 15 are predicted to be more likely to contribute to a disease risk as this score represents “the median value for all possible canonical splice site changes and nonsynonymous variants” [38].

Variants and genes shared across families: Analysis across all 23 families was performed to identify if there were any genes with rare, nonsynonymous or loss-of-function (LOF) variants in more than one family. Variants were selected that had a MAF ≤ 1% in the Kaviar Database and CADD scores ≥ 15 to try to identify potentially damaging alterations [39].

Association testing of top candidates

All top variants and genes from the three separate analyses described above were evaluated as potential risk variants using genomewide association statistics for two family study cohorts (NIA-LOAD and MIRAGE) in the Alzheimer Disease Genetics Consortium [42,43]. Both gene and SNP-based tests were adjusted for age, sex and principal components (PCs). SNP-based logistic regression tests in each study were performed in the SNPTest program, and meta-analysis of these results was conducted using METAL [21,44]. Gene-based tests were conducted on meta-analysis summary statistics using VEGAS [44]. Variants tested in the gene-based analysis included all variants with a MAF<5%.

Results

Variants identified in known AD genes

Each sequenced family contained between 4 and 16 individuals diagnosed with AD. Mean age-at-onset across all families was 74.3 years. We identified 14 potentially damaging variants in 10 known AD genes and GWAS implicated loci (Table 1). Seven of the variants were observed in multiple affected individuals in the same family, while the remaining variants were observed only once. All alterations were single nucleotide changes with the exception of a four base pair deletion in CD33. This deletion is potentially the most deleterious as it is predicted to causes a frameshift that encodes two incorrect amino acids before terminating prematurely, thus failing to generate over 40% of the protein. In addition, multiple variants were observed in four genes: AKAP9, INPP5D, SORL1 and UNC5C. Each gene had at least two families with a variant identified in it, while family 191 have a single affected individual with two alterations in UNC5. One of the variants in UNC5C, Ala860Thr, was identified in two different families; this alteration has a CADD score of 33, the highest score in this category.

Segregating variants in linkage regions

Linkage scans aggregating all families identified one primary linkage region, a parametric multipoint peak on chromosome 1q23 (161.9-165.6 MB). Two families had strong linkage in this region (family specific LOD>2). However, no variants met our filtering criteria for these two families, suggesting the causal variant(s) may be noncoding changes either removed from by our filtering criteria or not present in our WES. Three of the 23 families also have family-specific parametric LOD scores >2; rare, potentially damaging alterations in five genes occurred within these regions and may potentially be the strongest novel AD candidate genes (Table 2). The five alterations were all missense changes in CD163L1, CLECL1, CTNNA1, GALR3 and MIEF1.

Genes with variants in more than one family

We identified four genes that had rare (MAF ≤ 1%), segregating, and potentially deleterious variants in at least two families (Table 3). Three of these genes had the same missense alteration identified in distinct families: MKL2, PLEKHG5 and THBS2.

Table 3: Genes with rare, potentially damaging variants in more than one family.

Association testing of variants and genes

From our prioritized sets, a total of 9 SNPs and 14 genes were available for testing in the ADGC family-based meta-analysis datasets (NIA-LOAD and Mirage). None of the variants tested were significantly associated with disease (Supplemental Table 3). The gene MIEF1, identified as a candidate gene in a family 1201 with rare, potentially damaging segregating variants in a region with a LOD score of 2.22, reached nominal significance (p=0.049, Supplemental Table 4).

Discussion

Through WES of large families with a heavy burden of AD, variants in both known and novel loci were identified that could contribute to risk. Filtering for rare, segregating, and potentially damaging variants identified five novel candidate genes (Table 2). These genes encompass a variety of functions that are suggestive of a link to AD. For example, two of these genes are involved in regulating immunity: CD163L1 and CLECL1 [45-47]. CD163L1 is expressed in macrophages, upregulated in response to IL-10 and acts as an endocytic receptor [48]. CLECL1 is highly expressed in B cells and dendritic cells and may enhance the immune response through upregulation of IL-4 [46]. Neuroinflammation has been shown to occur in AD patients, possibly through misregulation of microglia and triggered by amyloid beta plaques [49]. Additionally, established AD risk genes, such as ABCA7, CD33 and TREM2, have also been linked to the immune system [4]. Another gene identified through this study, CTNNA1, encodes a catenin expressed at elevated levels in the nervous system [50]. GALR3 is a receptor for the neuropeptide galanin, which has been shown to modulate a variety of processes, including cognition and memory, functions disrupted in AD [51,52]. MIEF1 was nominally associated with late onset AD in a meta-analysis of two family datasets from the ADGC, thereby suggesting that it may play a wider role in AD that extends beyond a single multiplex AD family. MIEF1 may play a role in dysfunctional mitochondria and their potential to damage neurons [53,54]. Thus, each of the genes in the peak linkage regions are connected to known AD functions or neuronal pathways.

Four genes had rare, potentially damaging variants in more than one family (Table 3). When evaluating known functions of these genes, two are of particular interest, PLEKHG5 and THBS2. PLEKHG5 has been previously implicated in both Charcot-Marie-Tooth disease and spinal muscular atrophy [55,56]. PLEKHG5 is ubiquitously expressed throughout the nervous system and murine studies demonstrated lowered expression can alter the velocity of nerve conduction [55,56]. In addition, THBS2 is an intriguing novel AD candidate gene due to its involvement in synaptogenesis in immature astrocytes [57]. Further investigation into each of these novel AD candidates and the variants identified in this study is required.

After evaluating our families for rare alterations in known AD genes and loci, variants were discovered in genes previously connected to both early and late-onset AD (Table 1). Four genes had multiple alterations: AKAP9, INPP5D, SORL1 and UNC5C. Some of these alterations have the potential to interfere with a protein’s function due to their location within specific domains. For example, a rare alteration in UNC5C identified in two distinct families, Ala860Thr, falls within the highly conserved DEATH domain, a region composed of alphahelices and involved in apoptotic functions. Another study identified a different alteration in the same region in AD families and proposed that the alteration may increase the susceptibility of neurons to death [17]. A single affected individual in family 2349 was found to carry a frameshift deletion in CD33 predicted to remove over 40% of the protein. There is evidence that higher expression of CD33 in brains is associated with cognitive decline [58]. However, it may be that dysregulation of the protein, either through over or under expression, could contribute to AD risk. In addition, a rare alteration in SORL1, Thr588Ile, was identified within the vacuolar protein sorting 10 (VPS10) domain and may influence the processing of amyloid beta fragments, as has been shown for other AD associated variants in this gene [59]. Moreover, a potentially pathogenic alteration was identified in PSEN1 in a single individual, Glu318Gly; this variant was previously reported to result in higher tau and phosphorylated tau levels in cerebrospinal fluids [60]. Three affected individuals from family 1893 were found to share the Arg336His alteration in the NME8 gene, a change that fell within the first NDK domain of the protein. This gene has been associated with clinical features of AD including atrophy of the hippocampus and occipital gyrus [6,61]. These alterations, while not segregating within all affected individuals in the families, may play a contributing role in AD risk.

Conclusion

In conclusion, this study demonstrates how using large, extended families to evaluate exome data identifies segregating risk variants in potentially novel AD candidate genes. In contrast to GWAS studies that have grown from hundreds to thousands and tens of thousands of participants, this study design requires far fewer participants. Indeed, a single extended family may be sufficient to identify a novel AD candidate gene [18]. Moreover, WES has the sensitivity to directly detect both common and rare variants that may confer a risk to AD, while GWAS findings are limited to pinpointing a region of interest, but not necessarily the causative alterations. In the study presented here, rare changes potentially contributing to AD risk were found in genes implicated in the immune response, CD163L1 and CLECL1, and neuronal function, CTNNA1, GALR3, MIEF1, PLEKHG5 and THBS2. Variants were also identified in genes previously connected to both early and late-onset AD including AKAP9, INPP5D, SORL1 and UNC5C. Further investigation will be required to fully assess the cellular and molecular consequences of the alterations identified here as well as determine whether the novel genes found are involved in AD risk across larger datasets.

Acknowledgement

This work was supported by the National Institutes of Health (R01 AG027944, R01 AG028786 to MAPV, R01 AG019085 to JLH, P20 MD000546); a joint grant from the Alzheimer’s Association (SG-14-312644) and the Fidelity Biosciences Research Initiative to MAPV; the BrightFocus Foundation (A2011048 to MAPV). NIA-LOAD Family-Based Study supported the collection of samples used in this study through NIH grants U24 AG026395 and R01 AG041797 and the MIRAGE cohort was supported through the NIH grants R01 AG025259 and R01 AG048927. We thank contributors, including the Alzheimer’s disease Centers who collected samples used in this study, as well as patients and their families, whose help and participation made this work possible. Study design: HNC, BWK, JLH, MAPV; Sample collection: MLC, JMV, RMC, LAF, JLH, MAPV; Whole exome sequencing and Sanger sequencing: SR, PLW; Sequencing data analysis: HNC, BWK, KLHN, SR, MAK, JRG, ERM, GWB, MAPV; Statistical analysis: BWK, KLHN, JMJ, MAPV; Preparation of manuscript: HNC, BWK. The authors jointly discussed the experimental results throughout the duration of the study. All authors read and approved the final manuscript.

References

1. Qiu C, Kivipelto M, Von Strauss E (2009) Epidemiology of Alzheimer's disease: Occurrence, determinants and strategies toward intervention. Dialogues Clin Neurosci 11: 111-128.
2. Gatz M, Reynolds CA, Fratiglioni L, Johansson B, Mortimer JA et al. (2006) Role of genes and environments for explaining Alzheimer disease. Arch Gen Psychiatry 63: 168-174.
3. Van Cauwenberghe C, Van Broeckhoven C, Sleegers K (2015) The genetic landscape of Alzheimer disease: Clinical implications and perspectives. Genet Med 18: 421-430.
4. Corder EH, Saunders AM, Strittmatter WJ, Schmechel DE, Gaskell PC, et al. (1993) Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. Science 261: 921-923.
5. Strittmatter WJ, Saunders AM, Schmechel D, Pericak-Vance M, Enghild J, et al. (1993) ApolipoproteinE: High-avidity binding to beta-amyloid and increased frequency of type 4 allele in late-onset familial Alzheimer disease. Proc Natl Acad Sci U S A 90: 1977-1981.
6. Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, et al. (2013) Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer's disease. Nat Genet 45: 1452-1458.
7. Pottier C, Hannequin D, Coutant S, Rovelet-Lecrux A, Wallon D, et al. (2012) High frequency of potentially pathogenic SORL1 mutations in autosomal dominant early-onset Alzheimer disease. Mol Psychiatry 17: 875-879.
8. Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, et al. (2013) TREM2 variants in Alzheimer's disease. N Engl J Med 368: 117-127.
9. Guerreiro RJ, Lohmann E, Kinsella E, Bras JM, Luu N, et al. (2012) Exome sequencing reveals an unexpected genetic cause of disease: NOTCH3 mutation in a Turkish family with Alzheimer's disease. Neurobiol Aging 33: 1008.e17-1008.e23.
10. Jonsson T, Stefansson H, Steinberg S, Jonsdottir I., Jonsson, P.V., et al. (2013) Variant of TREM2 associated with the risk of Alzheimer's disease. N Engl J Med 368: 107-116.
11. Benitez BA, Jin SC, Guerreiro R, Graham R, Lord J, et al. (2014) Missense variant in TREML2 protects against Alzheimer's disease. Neurobiol Aging 35: 1510.e19-1510.e26.
12. Carney RM, Kohli MA, Kunkle BW, Naj AC, Gilbert JR, et al. (2014) Parkinsonism and distinct dementia patterns in a family with the MAPT R406W mutation. Alzheimers Dement 10: 360-365.
13. Cruchaga C, Karch CM, Jin SC, Benitez BA, Cai Y, et al. (2014) Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease. Nature 505: 550-554.
14. Dolzhanskaya, N, Gonzalez MA, Speriani F, Stefl S, Messing J, et al. (2014) A novel p.Leu(381)Phe mutation in presenilin 1 is associated with very early onset and unusually fast progressing dementia as well as lysosomal inclusions typically seen in Kufs disease. J Alzheimers Dis 39: 23-27.
15. Logue MW, Schu M, Vardarajan BN, Farrell J, Bennett DA, et al. (2014) Two rare AKAP9 variants are associated with Alzheimer's disease in African Americans. Alzheimers Dement 10: 609-618.e11.
16. Sassi, C, Guerreiro R, Gibbs R, Ding J, Lupton MK, et al. (2014) Investigating the role of rare coding variability in Mendelian dementia genes (APP, PSEN1, PSEN2, GRN, MAPT and PRNP) in late-onset Alzheimer's disease. Neurobiol Aging 35: 2881.e1-2881.e6.
17. Wetzel-Smith MK, Hunkapiller J, Bhangale TR, Srinivasan K, Maloney JA, et al. (2014) A rare mutation in UNC5C predisposes to late-onset Alzheimer's disease and increases neuronal cell death. Nat Med 20: 1452-1457.
18. Kohli MA, Cukier HN, Hamilton-Nelson KL, Rolati S, Kunkle, BW, et al. (2016) Segregation of a rare TTC3 variant in an extended family with late-onset Alzheimer disease. Neurol Genet 2: e41.
19. Velez JI, Lopera F, Sepulveda-FAlla D, Patel HR, Jphar AS, et al. (2016) APOE*E2 allele delays age of onset in PSEN1 E280A Alzheimer�s disease. Mol Psychiatry 21: 916-924.
20. Ehling R, Noskov�¡ L, Str�¡neck�½ V, Hartmannov�¡ H, Pristoupilov�¡ A, et al. (2013) Cerebellar dysfunction in a family harboring the PSEN1 mutation co-segregating with a cathepsin D variant p.A58V. J Neurol Sci 326: 75-82.
21. Louwersheimer E, Cohn-Hokke PE, Pijnenburg YA, Weiss MM, et al. (2017) Rare genetic variant in SORL1 may increase penetrance of Alzheimer's disease in a family with several generations of APOE-e4 Homozygosity. J Alzheimers Dis 56: 63-74.
22. Kunkle BW, Jaworski J, Barral S, Vardarajan B, Beecham GW, et al. (2016) Genome-wide linkage analyses of non-Hispanic white families identify novel loci for familial late-onset Alzheimer's disease. Alzheimers Dement 12: 2-10.
23. McKhann G, Drachman D, Folstein M, Katzman R, Price D, et al. (1984) Clinical diagnosis of Alzheimer's disease: Report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer's Disease. Neurology 34: 939-944.
24. McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack CR, et al. (2011) The diagnosis of dementia due to Alzheimer�s disease: Recommendations from the National Institutes on Aging-Alzheimer�s Association workgroups on diagnostic guidelines for Alzheimer�s disease. Alzheimers Dement 7: 263-269.
25. Albert MS, DeKosky ST, Dickson D, Dubois B, Feldman HH, et al. (2011) The diagnosis of mild cognitive impairment due to Alzheimer�s disease: Recommendations from the National Institute on Aging-Alzheimer�s Association workgroups on diagnostic guidelines for Alzheimer�s disease. Alzheimers Dement 7: 270-279.
26. Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754-1760.
27. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, et al. (2010) The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20: 1297-1303.
28. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, et al. (2011) A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 43: 491-498.
29. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, et al. (2002) The human genome browser at UCSC. Genome Res 12: 996-1006.
30. Raney BJ, Dreszer TR, Barber GP, Clawson H, Fujita PA, et al. (2014) Track data hubs enable visualization of user-defined genome-wide annotations on the UCSC genome browser. Bioinformatics 30: 1003-1005
31. Ye K, Schulz MH, Long Q, Apweiler R, Ning Z (2009) Pindel: A pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics 25: 2865-2871
32. Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nat Methods 9: 357-359
33. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559-575
34. Abecasis GR, Cherny SS, Cookson WO, Cardon LR (2002) Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 30: 97-101
35. Wang K, Li M, Hakonarson H (2010) ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 38: e164
36. Li MX, Gui HS, Kwan JS, Bao SY, Sham PC (2012) A comprehensive framework for prioritizing variants in exome sequencing studies of Mendelian diseases. Nucleic Acids Res 40: e53
37. Tan A, Abecasis GR, Kang HM (2015) Unified representation of genetic variants. Bioinformatics 31: 2202-2204
38. Kircher M, Witten DM, Jain P, O'Roak BJ, Cooper GM, et al. (2014) A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet 46: 310-315
39. Glusman G, Caballero J, Mauldin DE, Hood L, Roach JC (2011) Kaviar: An accessible system for testing SNV novelty. Bioinformatics 27: 3216-3217
40. Shearer AE, Eppsteiner RW, Booth KT, Ephraim SS, Gurrola J, et al. (2014) Utilizing ethnic-specific differences in minor allele frequency to recategorize reported pathogenic deafness variants. Am J Hum Genet 95: 445-453
41. Richards S, Aziz N, Bale S, Bick D, Das S, et al. (2015) Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med 17: 405-424
42. Green RC, Cupples LA, Go R, Benke KS, Edeki T, et al. (2002) Risk of dementia among white and African American relatives of patients with Alzheimer disease. JAMA 287: 329-336
43. Lee JH, Cheng R, Graff-Radford N, Foroud T, Mayeux R, et al. (2008) Analysis of the National Institute on Aging Late-Onset Alzheimer�s Disease Family Study: Implication of additional loci. Arch Neurol 65: 1518-1526
44. Willer CJ, Li Y, Abecasis GR (2010) METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26: 2190-2191
45. Liu JZ, McRae AF, Nyholt DR, Medland SE, Wray NR, et al. (2010) A versatile gene-based test for genome-wide association studies. Am J Hum Genet 87: 139-145
46. Ryan EJ, Marshall AJ, Magaletti D, Floyd H, Draves KE, et al. (2002) Dendritic cell-associated lectin-1: A novel dendritic cell-associated, C-type lectin-like molecule enhances T cell secretion of IL-4. J Immunol 169: 5638-5648
47. Gonzalez-Dominguez E, Samaniego R, Flores-Sevilla JL, Campos-Campos SF, Gomez-Campos G, et al. (2015) CD163L1 and CLEC5A discriminate subsets of human resident and inflammatory macrophages in vivo. J Leukoc Biol 98: 453-466
48. Moeller JB, Nielson MJ, Reichhardt MP, Schlosser A, Sorensen GL, et al. (2012) CD163-L1 is an endocytic macrophage protein strongly regulated by mediators in the inflammatory response. J Immunol 188: 2399-2409
49. Minter MR, Taylor JM, Crack PJ (2016) The contribution of neuroflammation to amyloid toxicity in Alzheimer'd disease. J Neurochem 136: 457-474
50. Kask M, Pruunsild P, Timmusk T (2011) Bidirectional transcription from human LRRTM2/CTNNA1 and LRRTm1/CTNNA2 gene loci leads to expression of N-terminally truncated CTNNA1 and CTNNA2 isoforms. Biochem Biophys Res Commun 411: 56-61
51. Kolakowski LF, O'Neill GP, Howard AD, Broussard SR, Sullivan KA, et al. (1998) Molecular characterization and expression of cloned human galanin receptors GALR2 and GALR3. J Neurochem 71: 2239-2251
52. Steiner RA, Hohmann JG, Holmes A, Wrenn CC, Cadd G, et al. (2001) Galanin transgenic mice display cognitive and neurochemical deficits characteristic of Alzheimer's disease. Proc Natl Acad Sci U S A 98: 4184-4189
53. Wu F, Yao PJ (2009) Clathrin-mediated endocytosis and Alzheimer's disease: An update. Ageing Res Rev 8: 147-149
54. Garcia-Escudero V, Martin-Maestro P, Perry G, Avila J (2013) Deconstructing mitochondrial dysfunction in Alzheimer disease. Oxid Med Cell Longev 2013: 162152
55. Maystadt I, Rezsohazy R, Barkats M, Duque S, Vannuffel P, et al. (2007) The nuclear factor kappaB-activator gene PLEKHG5 is mutated in a form of autosomal recessive lower motor neuron disease with childhood onset. Am J Hum Genet 81: 67-76
56. Azzedine H, Zavadakova P, Plante-Bordeneuve V, Vaz Pato M, Pinto N, et al. (2013) PLEKHG5 deficiency leads to an intermediate form of autosomal-recessive Charcot-Marie-Tooth disease. Hum Mol Genet 22: 4224-4232
57. Christopherson KS, Ullian EM, Stokes CC, Mullowney CE, Hell JW, et al. (2005) Thrombospondins are astrocyte-secreted proteins that promote CNS synaptogenesis. Cell 120: 421-433
58. Karch CM, Jeng AT, Nowotny P, Cady J, Cruchaga C et al. (2012) Expression of novel Alzheimer's disease risk genes in control and Alzheimer's disease brains. PLoS ONE 7: e50976
59. Vardarajan BN, Zhang Y, Lee JH, Cheng R, Bohm C, et al. (2015) Coding mutations in SORL1 and Alzheimer disease. Ann Neurol 77: 215-227
60. Benitez BA, Karch CM, Cai Y, Jin SC, Cooper B, et al. (2013) The PSEN1, p.E318G variant increases the risk of Alzheimer's disease in APOE-epsilon4 carriers. PLoS Genet 9: e1003685
61. Liu Y, Yu JT, Wang HF, Hao XK, Yang YF, et al. (2014) Association between NME8 locus polymorphism and cognitive decline, cerebrospinal fluid and neuroimaging biomarkers in Alzheimer's disease. PLoS ONE 9: e114777