Identification of Alzheimer’s disease– associated rare coding variants in the ECE2 gene

Accumulation of amyloid b protein (Ab) due to increased generation and/or impaired degradation plays an important role in Alzheimer’s disease (AD) pathogenesis. In this report, we describe the identification of rare coding mutations in the endothelin-converting enzyme 2 (ECE2) gene in 1 late-onset AD family, and additional case-control cohort analysis indicates ECE2 variants associated with the risk of developing AD. The 2 mutations (R186C and F751S) located in the peptidase domain in the ECE2 protein were found to severely impair the enzymatic activity of ECE2 in Ab degradation. We further evaluated the effect of the R186C mutation in mutant APP–knockin mice. Overexpression of wild-type ECE2 in the hippocampus reduced amyloid load and plaque formation, and improved learning and memory deficits in the AD model mice. However, the effect was abolished by the R186C mutation in ECE2. Taken together, the results demonstrated that ECE2 peptidase mutations contribute to AD pathogenesis by impairing Ab degradation, and overexpression of ECE2 alleviates AD phenotypes. This study indicates that ECE2 is a risk gene for AD development and pharmacological activation of ECE2 could be a promising strategy for AD treatment. Research Article Neuroscience

Here, we performed WES in a large LOAD pedigree without any known AD causative gene mutations and found that the R186C mutation encoded in the ECE2 gene was segregated with disease status in 1 branch of the pedigree. Furthermore, an excess burden of ECE2 deleterious coding mutations that contributed to AD pathogenesis was found through screening whole ECE2 gene exomes in a small cohort Accumulation of amyloid β protein (Aβ) due to increased generation and/or impaired degradation plays an important role in Alzheimer's disease (AD) pathogenesis. In this report, we describe the identification of rare coding mutations in the endothelin-converting enzyme 2 (ECE2) gene in 1 late-onset AD family, and additional case-control cohort analysis indicates ECE2 variants associated with the risk of developing AD. The 2 mutations (R186C and F751S) located in the peptidase domain in the ECE2 protein were found to severely impair the enzymatic activity of ECE2 in Aβ degradation. We further evaluated the effect of the R186C mutation in mutant APP-knockin mice. Overexpression of wild-type ECE2 in the hippocampus reduced amyloid load and plaque formation, and improved learning and memory deficits in the AD model mice. However, the effect was abolished by the R186C mutation in ECE2. Taken together, the results demonstrated that ECE2 peptidase mutations contribute to AD pathogenesis by impairing Aβ degradation, and overexpression of ECE2 alleviates AD phenotypes. This study indicates that ECE2 is a risk gene for AD development and pharmacological activation of ECE2 could be a promising strategy for AD treatment.
of Han Chinese ancestry. Subsequently, functional studies demonstrated that ECE2 mutations impaired its enzymatic activity in degrading Aβ in HEK293 cells. We further demonstrated that the total Aβ load in AD model mice was reduced, and impaired learning and memory functions were rescued by overexpression of wild-type but not mutant ECE2 protein. In summary, our study demonstrates that ECE2 rare variants may increase AD risk through interfering with the degradation of Aβ degradation, which suggests that ECE2 is a potential target for AD treatment.

Results
The clinical phenotypes of the family with LOAD. We identified a large family with 4 AD patients over 2 generations (family A) ( Figure 1A). All affected individuals and healthy members of the third generation were assessed in at least 1 medical history review and thorough neurological examination by 2 experienced neurologists, and the clinical features of the 3 patients are summarized in Supplemental Table 1 (supplemental material available online with this article; https://doi.org/10.1172/jci.insight.135119DS1). The average age at onset of the patients was 66 ± 7.9 years, and the average disease duration was 11 ± 1.7 years. The proband (III:3) began to develop episodic memory loss at the age of 69. He often repeated words, misplaced things, and had difficulty in remembering the names of friends. Over the next 7 years, his symptoms worsened with a gradual loss of calculation abilities, orientation in time and space, and perceptivity. He also presented with personality changes characterized by irritability, depression, and barely talked to others. He relied entirely on others for care after 8 years of onset, including bathing, dressing, eating and other daily tasks. He died 10 years after the onset of the disease. The clinical features of patients III:1 and III:5 were similar to those of the proband, with progressive cognitive decline and personality changes. Neurological examination of the 3 patients revealed no myoclonic jerks, seizures, extrapyramidal or upper motor neuron signs. Structural magnetic resonance imaging (MRI) showed widespread cortical and hippocampal atrophy in patient III:3 ( Figure 1, B and C), and lacunar stroke and medial temporal lobe atrophy in patient III:5 (Supplemental Figure 1). The proband had increased 11 C-Pittsburgh compound B (PiB) accumulation in the frontal lobe, temporal lobe, and occipital lobe by PET imaging ( Figure 1D). 18 F-fluorodeoxyglucose ( 18 F-FDG) PET examination showed hypometabolic activity in temporoparietal, frontal, and occipital cortices of patient III:3 ( Figure 1E). Senile plaques and neurofibrillary tangles were clearly detected in the proband by autopsy (Figure 1, F and G). There was no accumulation of α-synuclein in these brain tissues (Supplemental Figure 2) Table 2).
Exome sequencing identified ECE2 as a candidate risk gene for AD. Initial genotyping results showed that the proband AD patient carried no pathogenic mutations in the APP, PSEN1, PSEN2, GRN, MAPT, NOTCH3, PPNR, or C9ORF72 gene (Supplemental Figure 3). The APOE genotypes of the 3 patients (III:1, III:3, and III:5) and 1 control (III:2) were all ε3/ε3, which means that the APOE genotype cannot account for the phenotype of this pedigree (Figure 2A). To further determine the genetic contributions, exome sequencing was performed on the 3 affected individuals (III:1, III:3, and III:5) and 1 control (III:2) in family A (Supplemental Table 3). After SNP and indel calling, we identified over 66 million reads/person, of which 98% mapped to the exome. Variants identified with a quality score lower than 20 and a depth of coverage lower than 5 were removed. No previously reported AD-associated pathogenic mutation was identified in this family by screening the Human Gene Mutation Database (HGMD) and ClinVar database. Because family A showed possible autosomal dominant inheritance, the dementia-associated mutation in this family should not be common in the general population. We then removed the nonsynonymous/ splice acceptor and donor site/insertions or deletions (NS/SS/indel) variants with minor allele frequency (MAF) greater than 0.5% in the following databases: Genome Aggregation Database (gnomAD), ExAc database, the 1000 Genomes Project, the 8 previously exome-sequenced HapMap samples (HapMap 8), the Beijing Genomics Institute (BGI) in-house exome variant data set, and the PVFD database (variant frequency in normal population database). More than 2200 rare variants were identified per individual. Twenty variants were shared by the 3 affected individuals and absent in the unaffected individual within the family. Four variants were further removed due to their known associated clinical symptoms not matching the phenotype of family A (Supplemental Table 4). Genotyping the remaining 16 variants in 10 controls (III:4, III:6, and III:7-III:14) of family A using Sanger sequencing excluded an additional 13 variants because they did not cosegregate with disease status. Finally, 3 variants in 3 genes segregated with disease were selected as candidate mutations (Supplemental Table 5). The potential pathogenicity of the identified missense variants was evaluated by ReVe Score (23), PolyPhen-2 (Polymorphism Phenotyping v2; http://genetics.bwh.harvard.edu/pph/) (24), Sorting Intolerant from Tolerant (SIFT) (http:// sift.jcvi.org/) (25), and Mutation Taster (26). The variants, D378N in YEATS2 and H462Y in APEH, were evaluated to be "tolerable" by ReVe and "Polymorphism" by Mutation Taster. The variant rs368866385 (Arg186Cys, ECE2, NM_001100120.1) with the highest ReVe score and evaluated to be damaging by 3 in silico prediction software packages was picked up as a promising candidate AD risk variant ( Figure 2B). Arginine 186 is located in the M13 peptidase domain of ECE2 and a highly conserved residue across species (Figure 2, C and D).
Gene-based association test confirmed ECE2 as a promising candidate risk gene for AD. To determine if there are additional AD-associated variants in the ECE2 gene, we sequenced the coding region of ECE2 in 741 AD cases and 545 controls. Among the 741 cases, 72 were diagnosed with AD through cerebrospinal fluid testing (n = 33), PET-CT (n = 38), or autopsy (n = 1). We identified 29 high-quality SNPs that met predefined quality criteria. None of the polymorphic ECE2 coding variants was in linkage disequilibrium in  (Table 1). When using more strict criteria to extract likely damaging variants (ReVe value > 0.7), the SKAT-O result remained significant (P = 0.01; OR = 9.71, 95%CI = 1.27-74.48) (Supplemental Table 6). In contrast, we did not observe a difference in synonymous variants for AD cases in our data set compared with controls (1.12-fold, P = 0.51).
Our results indicate that the association between nonsynonymous variant burden in ECE2 and AD risk is positive, which further confirmed that ECE2 is a promising candidate risk gene for AD.
ECE2 R186C and ECE2 F751S mutants reduced Aβ degradation in vitro. ECE2 was reported to degrade Aβ (28), and ECE2 homozygous knockout mice display significantly elevated endogenous Aβ40 and Aβ42 (29). As shown in Supplemental Table 6, F751S is the only variant predicted to be damaging by bioinformatics analysis, and is carried by both AD patients and normal control; therefore, further biological experiments were warranted. To investigate the effect of ECE2 R186C and F751S mutations on Aβ degradation, the Swedish-mutant-APP stable cell line 20E2 was transfected with ECE2 variants. Compared with empty vector control, overexpression of wild-type ECE2 (ECE2 WT ) significantly decreased the levels of secreted Aβ40 and Aβ42 in conditioned media to 64.93% ± 9.77% and 45.10% ± 6.38%, respectively (P < 0.01) ( Figure 3, A and B), whereas ECE2 R186C and ECE2 F751S had no effect on Aβ40 and Aβ42 levels (P > 0.05). In cell extracts, ECE2 WT also showed markedly reduced intracellular Aβ40 levels, to 58.07% ± 3.48% (P < 0.05), while the intracellular Aβ40 levels in the ECE2 R186C -and ECE2 F751S -transfected cells were similar to that in the control group (P > 0.05) ( Figure 3C). The data suggested that unlike ECE2 WT , overexpression of ECE2 R186C and ECE2 F751S mutants had no effect on Aβ degradation in 20E2 cells.
R186C and F751S mutations are located in the M13 peptidase domain of the ECE2 protein. Immunohistochemical analyses showed that these 2 mutations did not alter distribution and subcellular location of the mutant ECE2 proteins (Supplemental Figures 4 and 5). To examine whether the mutations impair ECE2's metalloprotease activity, leading to the deficiency in Aβ degradation, the activity of ECE2 variants was measured using an McaBk2 assay (30). The results showed that the relative fluorescence units (RFU) value of ECE2 WT was significantly higher than that of ECE2 R186C and ECE2 F751S ( Figure 3D, P < 0.0001), indicating ECE2 WT degraded the substrate much faster. There was no significant difference in RFU value between ECE2 R186C , ECE2 F751S , and empty vector control, indicating that R186C and F751S diminished ECE2's enzymatic activity (Figure 3, D and E). To investigate whether the ECE2 R186C and ECE2 F751S mutants had effects on Aβ production, the levels of APP and its C-terminal fragments were assayed by Western blot with C20 antibody. No significant difference was detected in C83 and C99 level between groups, suggesting that overexpressed ECE2 WT , ECE2 R186C , and ECE2 F751S did not alter APP processing to generate Aβ (Figure 3, F-H). These data indicated that both mutations R186C and F751S caused loss of ECE2's ability to degrade Aβ, resulting in Aβ accumulation.
ECE2 WT , not ECE2 R186C , reduced plaque formation in AD model mice. To investigate the effects of ECE-2 WT and the ECE2 R186C mutant on amyloid degradation and plaque formation, adeno-associated viruses (AAVs) carrying human ECE2 WT or ECE2 R186C were injected bilaterally into the ventricles of neonatal APP NL-G-F/NL-G-F -knockin mice on postnatal day 0 (P0). At 5 months after injection, the exogenous expression of Myc-tagged ECE2 proteins was detected throughout the brain, with dense labeling apparent in cerebral cortex, hippocampus, and subcortical nuclei (Figure 4, A and C). The average immunofluorescence intensity of ECE2 WT was similar to ECE2 R186C in whole brain sections (P > 0.05, Figure 4B). Overexpression of ECE2 WT significantly reduced plaque formation in rAAV9-ECE2 WT -injected mice compared with rAAV9-GFP-injected mice (2.73% ± 0.08% vs. 3.16% ± 0.09%, P < 0.01) (Figure 4, A and D). However, APP NL-G-F/NL-G-F -knockin mice injected with rAAV9-ECE2 R186C exhibited plaque loads similar to those in mice that received rAAV9-GFP (P > 0.05), and an increased plaque formation when compared with rAAV9-ECE2 WT -injected mice (3.36% ± 0.10% vs. 2.73% ± 0.08%, P < 0.001) (Figure 4, A and D). ELISAs confirmed that ECE2 WT overexpression resulted in a significant decrease in Aβ42 levels of 63% ± 1.9% in the hippocampal tissues from rAAV9-ECE2 WT -injected mice compared with GFP control (P < 0.001) ( Figure 4E). No differences in Aβ42 level were observed between mice that received ECE2 R186C and GFP control virus injection (P > 0.05). The expression level of ECE2 WT and ECE2 R186C in bilateral hippocampus was verified to be similar by Western blot (P > 0.05) ( Figure 4F). These data demonstrated that the expression of ECE2 WT reduced the amyloid plaque formation but ECE2 R186C abolished ECE2's Aβ-hydrolyzing function and its inhibitory effect on plaque formation in the brains of APP NL-G-F/NL-G-F -knockin mice in vivo. ECE2 WT , not ECE2 R186C , rescued the cognitive deficits in AD model mice. To examine the effects of ECE2 WT and the ECE2 R186C mutant on learning and memory in vivo, a Morris water maze test for spatial learning and memory and a Y-maze test for working memory were performed on the virus-injected APP NL-G-F/NL-G-Fknockin mice at 5 months of age. In the visible platform test of the Morris water maze test, no significant differences were observed in escape latency (P > 0.05) ( Figure 5A) and path length (P > 0.05) ( Figure 5B) between the GFP, ECE2 WT , and ECE2 R186C groups, indicating that all mice had similar motility and vision function. In the hidden platform tests, the rAAV9-ECE2 WT -injected mice had relatively shorter escape latency and path length than the rAAV9-ECE2 R186C -injected mice and the GFP control, but no significant differences were observed ( Figure 5, C and D). However, in the probe trial, the rAAV9-ECE2 WT -injected mice traveled significantly more times into the third quadrant where the hidden platform was previously placed than control mice (2.33 ± 0.41 vs. 0.80 ± 0.20, P < 0.05) ( Figure 5E). In contrast, the rAAV9-ECE2 R186C -injected mice had similar passing times compared with rAAV9-GFP-injected control mice (P > 0.05). We also performed a Y-maze test on the mice, and the results showed that the rAAV9-ECE-2 WT -injected mice exhibited a significantly higher percentage of spontaneous alternation than the rAAV9-GFP-injected mice (74.00% ± 3.47% vs. 57.53% ± 1.78%, P < 0.001), whereas rAAV9-ECE2 R186C -injected mice and rAAV9-GFP-injected mice performed similarly in the Y-maze test (59.57% ± 0.95% vs. 57.53% ± 1.78%, P > 0.05) ( Figure 5F). No difference between total arm entries was observed between these groups (P > 0.05) ( Figure 5G). Taken together, these data demonstrated that overexpression of ECE-2 WT protein, but not ECE2 R186C , alleviated learning and memory deficits in APP NL-G-F/NL-G-F -knockin mice. . All results are expressed as mean ± SEM. Two-tailed Student's t test was used to analyze the difference between 2 groups, and multiple comparisons were analyzed by ANOVA followed by Bonferroni's multiple-comparisons test.

Discussion
Recent studies showed that low-frequency coding variants not detected by GWAS may be a source of functional variants, with significant effects on LOAD risk. Next-generation sequencing technology provides powerful tools for identifying such rare variations. However, the frequency of rare coding variants is relatively low, and the impact of a single rare variation is often moderate in disease progression, which poses challenges to sample size and statistical method. One solution is to use highly selected pedigrees with increased genetic risk for disease to identify candidate risk genes, followed by a targeted resequencing of the gene in a case-control cohort. Here, we used WES to identify the ECE2 R186C mutation segregating with AD phenotypes in one branch of a LOAD pedigree. Although we cannot exclude the possibility that other variants may also contribute to the AD feature of this family, ECE2 R186C was the only variant we identified to be an AD-associated gene variant. Mutation analysis by using 7 in silico missense

Figure 5. Overexpressed ECE2 WT , but not ECE2 R186C mutant, alleviated cognitive deficits in APP NL-G-F/NL-G-F -knockin mice. (A-E)
A Morris water maze test consisted of 1 day of visible platform trials and 4 days of hidden platform trials, plus a probe trial 24 hours after the last hidden platform trial. Animal movement was traced and recorded by ANY-maze tracking software. In the visible platform tests, mice in each group exhibited a similar (A) escape latency and (B) swimming path length to escape onto the visible platform. In hidden platform tests, mice were trained with 5 trials per day for 4 days. ECE2 WT group mice showed shorter latency (C) and path length (D) to escape onto the visible platform on the third and fourth days, but without statistically significant differences. (E) In the probe trial on the sixth day, ECE2 WT group mice traveled into the third quadrant, where the hidden platform was previously placed, significantly more times than GFP and ECE2 R186C groups. (F and G) The Y-maze test was performed. ECE2 WT -injected mice exhibited significantly higher percentage alternation than GFP-injected mice, whereas ECE2 R186C -injected mice and GFP-injected mice performed similarly (F). No difference between arm entries was observed between groups (G). Data represented as mean ± SEM (for behavior tests, n = 10 GFP, 9 ECE2 WT , and 10 ECE2 R186C ). Significance was assessed by 1-way ANOVA with Newman-Keuls post hoc test.
prediction programs (SIFT, Polyphen-2_HDIV, Polyphen-2_HVAR LRT, MutationTaster, Provean, and ReVe) suggests that this mutation would potentially result in altered ECE2 protein function. Additionally, this mutation was absent in 545 Chinese controls with normal cognition, and showed an ultrarare allele frequency in a public database (0.0003 and 0 in gnomAD East Asian and African American groups, respectively). Subsequently, our study showed that the R186C mutation impaired ECE2's ability to degrade Aβ. By analyzing all ECE2 rare coding mutations identified in an AD case-control cohort of Han Chinese ancestry, the SKAT-O test revealed that nonsynonymous rare coding mutations in the ECE2 gene were enriched in AD patients by 5-fold compared with control, whereas synonymous mutations in the ECE2 gene were comparable between AD cases and the control cohort (P = 0.51). This observation provided robust evidence that ECE2 is an AD risk gene, and nonsynonymous mutations may play an important role in AD pathogenesis. Although a targeted exome and genome sequencing analysis for APP-Aβ metabolism genes in a cohort from North America and the United Kingdom did not identify an obvious correlation between ECE2 and AD pathogenesis (22), this difference might be due to the different ethnicities and small sample size of the 2 projects. Indeed, investigations on a larger sample size from multiple populations will be required in a future study to clarify the role of ECE2 in AD pathogenesis.
ECE family member proteins include ECE1 and -2, and encode membrane-bound metalloproteases mediating big endothelin conversion into vasoactive endothelins (31,32). ECE2 is a type II transmembrane protein with a short N-terminal cytoplasmic domain, a transmembrane domain, and a long M13 peptidase domain within the C-terminal extracellular region (33). The mutations R186C and F751S identified in the AD patients were located in the M13 domain that contains an active peptidase domain. Consistent with previous reports (18,34), we showed that overexpression of ECE2 WT possessed sufficient Aβ-degrading ability to reduce Aβ40 and Aβ42 levels. However, overexpression of the R186C or F751S mutant could not reduce Aβ40 and Aβ42 levels. In vitro enzymatic activity analysis revealed that the R186C or F751S mutation diminished the enzymatic activity of ECE2, contributing to the impaired Aβ degradation. Furthermore, our in vivo studies demonstrated that the R186C mutation caused loss of ECE2 function in the mouse model, which indicated that ECE2 R186C protein might lead to Aβ accumulation to form neuritic plaques and cognitive deficits.
It has been demonstrated that both Aβ40 and Aβ42 could enhance vascular tone by inducing endothelin 1 (ET-1) expression in rat aorta rings (35), which indicated that Aβ might reduce blood flow in cerebral vessels, leading to Aβ accumulation, cerebral vessel destruction, and finally neurodegeneration (36)(37)(38)(39)(40). Given that ECE2 inactivation could lead to Aβ accumulation in the central nervous system, it is possible that mutant ECE2 reduces cerebral blood flow via an Aβ-mediated ET-1 pathway and finally impairs cognitive function. Additionally, it has been reported that ECE2 is mainly expressed in somatostatin-expressing (SOM + ) interneurons (41), which are essential for cognitive function. In addition, Aβ aggregates may be toxic to SOM + interneurons (42)(43)(44) and SOM innervation is significantly impaired during AD (45). These results suggested that Aβ accumulation induced by ECE2 mutations may lead to the loss of SOM + interneurons and cognitive deficits in AD. Therefore, it is important to explore the effect of ECE2 mutations on cerebral blood flow and the functions of SOM + interneurons in the mutation carriers, which might provide valuable insights into the ECE2-associated AD pathogenesis.
We have identified an AD risk gene, ECE2, and provided evidence that low-frequency variants exert an intermediate to high effect on facilitating AD pathogenesis. Our study suggested that ECE2 is a potential therapeutic target for AD treatment.

Methods
A LOAD family for candidate gene discovery. Family A was a 3-generation Chinese pedigree with 4 members presenting as possible autosomal dominant inheritance ( Figure 1A). All family members of generation III completed a clinical assessment, neurological examination, and neuropsychological tests (including MoCA and MMSE). Patients III-3 and III-5 had blood tests to measure thyroid hormones, vitamin B12, and treponema pallidum hemagglutination test (TPHA). Patients III-3 and III-5 also underwent MRI to evaluate brain atrophy and other conditions, while only patient III-3 received 18 F-FDG/PiB-PET scans. Brain autopsy was performed on patient III-3 after death within 3 hours. Immunohistochemistry was performed using the following antibodies: tau (AT8, Invitrogen, catalog MN1020); β-amyloid (6E10, Sigma-Aldrich, catalog A1474); and α-synuclein (LB509, Abcam, catalog ab27766).
Exome sequencing. Exome-sequencing was performed as previously reported (46,47). Genomic DNA from patients was sheared by the Covaris system, then hybridized to NimbleGen's SeqCap EZ exome library v2.0 for enrichment according to the manufacturer's instructions. The enriched library targeting the exome was sequenced on the Illumina HiSeq 2000 platform for paired-end reads with read length of 90 bp. The reads were aligned to the reference coding sequence using the BWA (Burrows Wheeler Aligner) software. After alignment, the output files were used to perform sequencing coverage and depth. Genome Analysis Toolkit (GATK) software was used to detect single-nucleotide variants (SNVs) and indels. Variants identified with a quality score lower than 20 and a depth of coverage lower than 5 were removed. The HGMD (http:// www.hgmd.cf.ac.uk/ac) and ClinVar database (48) were used to screen reported mutations in previously published studies.
Participants for targeted resequencing study. The study recruited 741 AD patients (445 female and 296 male) with a mean age at onset of 66.20 ± 10.65 years from the Memory Clinic of Xiangya Hospital. Three hundred twelve (42.10%) AD patients carried the ApoE-ε4 allele. All patients were assessed in a thorough neurologic examination and neuropsychological assessment by 2 experienced neurologists. A diagnosis of probable or definite AD was determined following the NINCDS-ADRDA criteria (49). Ninety-six of 741 patients (12.96%) had a positive family history for AD, but only the probands were included in the SKAT-O analysis. AD patients carrying known mutations in APP, PSEN1, PSEN2, GRN, MAPT, NOTCH3, C9ORF72, or PPNR were excluded. Five hundred forty-five (261 female and 284 male) controls were recruited from the partners of patients visiting our Memory Clinic or geographically matched community-recruited volunteers. The mean age at inclusion was 68 ± 7.36 years, and 73 (13.39%) of the control subjects were ApoE-ε4 positive. Subjective memory complaints, neurological or psychiatric antecedents, and familial history of neurodegeneration were ruled out in the controls by a systematic medical interview. Cognitive screening was done in the controls with the MMSE (50,51).
Targeted resequencing of the ECE2 gene. Resequencing of ECE2 exons was performed on 741 AD cases and 545 controls using targeted capture and MPS. The baits, a pool of 18,167 individually synthesized 5′-biotinylated 120-bp RNA oligonucleotides, covered the untranslated regions (UTRs), coding regions, and exon-intron boundaries of 93 dementia-related genes, including APP, PSEN1, PSEN2, GRN, MAPT, NOTCH3, PPNR, ECE2, etc. Target enrichment was done with HaloPlex technology (Agilent Technologies). After quality control, captured library sequencing was carried out on Illumina HiSeq × Ten Analyzers following the manufacturer's standard sequencing protocols for 150 cycles per read to generate paired-end reads. Image analysis, error estimation, and base calling were performed using Illumina Pipeline software to generate raw data. We used AfterQC (52) to generate clean reads for further analysis. The clean reads (with a length of 150 bp) derived from targeted sequencing and filtering were then aligned to the human genome reference (hg19) using BWA software (53). GATK software was used to detect SNVs and indels. Both variant and reference genotypes were called using the same parameters. When reaching a quality score greater than 20 and coverage greater than 5×, a variant or reference gene call was made. Numbering of variations at the gDNA level was based on the GenBank accession number NC_000003.11, transcript level on NM_001100120.1.
Gene-based analysis with SKAT-O. The SKAT-O test (27,54) was implemented in R using SKAT v1.0.9 to determine the difference in the aggregate burden of ECE2 gene variants between AD cases and controls. An empirical P value was derived from the distribution of null results based on 10,000 permutation trials in which case/control assignment was randomized. SKAT-O statistics were performed on 2 subsets of variants that were created based on impact of the variants. One is synonymous variants, the other one is nonsynonymous variants, including frameshift, missense, in-frame deletion, and nonsense variants. Also, any bias introduced in the test due to an unbalanced case-control set was compared to observations from the synonymous variants class that was used to set the background expectation.
Plasmids and virus. Human ECE2 variant 5, preferentially expressed in brain (33), was subcloned into the pcDNA4/myc-His mammalian expression vector (Invitrogen) to generate pcDNA4/Myc-HisA-ECE2-WT. Plasmids carrying the ECE2 R186C or F751S mutation were generated by PCR-based site-directed mutagenesis using pcDNA4/Myc-HisA-ECE2-WT as a backbone according to the manufacturer's instructions. All constructs were verified by Sanger sequencing. The WT and R186C mutant ECE2 with Myc tag were introduced into the pFB-AAV9-CAG-WPRE-SV40pA vector under the control of the chicken β-actin promoter for AAV virus production. The recombinant viruses were generated and purified by Virovek. Infectious rAAV particles were expressed as vector genomes (vg/mL).
Western blot and immunohistochemical analysis. Cells were lysed in cold RIPA-DOC buffer (0.1% SDS, 1% sodium deoxycholate, 1% Triton X-100, 0.15 M NaCl, 0.05 M Tris-HCl, pH 7.2) plus protease inhibitor cocktail (Roche). Cell lysates were analyzed by SDS polyacrylamide gel electrophoresis (SDS-PAGE) with 16% Tris-tricine or 7.5% Tris-glycine gels, and then transferred to Immobilon-FL PVDF membranes (Millipore). Antibody 9E10 (Abcam, catalog ab32) was used to detect Myc-tagged proteins. APP and its C-terminal fragment products were detected with C20 antibody (57). β-Actin was detected with monoclonal antibody AC-15 (Sigma-Aldrich). The membranes were incubated with IRDye 800CW-labeled goat anti-mouse (Li-COR Biosciences, 926-32212) or anti-rabbit secondary antibodies (Li-COR Biosciences, 925-32213) and scanned on the Odyssey system (LI-COR Biosciences). For amyloid plaque detection, the mouse brains were fixed in 4% paraformaldehyde and sectioned at 30 μm in a Leica cryostat. Monoclonal antibody 6E10 (Sigma-Aldrich, catalog A1474) against amino acids 1-16 of human Aβ was used for plaque detection. Monoclonal antibody 9E10 was used to detect ECE2 expression. 4′,6-Diamidino-2-phenylindole (DAPI) was used for visualization of nuclei. Images were acquired using an inverted fluorescence microscope (Zeiss Axio Observer). A total of 24 coronal brain slices (180 μm apart; every sixth section) per mouse and a total of 10 animals per group were used for staining. Plaques were quantified by percentage of staining area per slice for each mouse. For quantification of WT and R186C mutant ECE2 in brain sections, the average signal intensity of Myc-tagged protein was measured by ImageJ (NIH) with subtraction of mean background intensity.
Aβ ELISA. Conditioned medium from cell culture was collected, and protease inhibitors and 4-(2-aminoethyl)-benzenesulfonyl fluoride (AEBSF) (Roche) were added to prevent Aβ degradation. For intracellular Aβ detection, cells were lysed in cold cell extraction buffer (1× PBS with 0.1% Triton X-100 and protease inhibitors), briefly sonicated, and cleared by centrifugation at 12000 g at 4°C for 15 minutes. For mouse brain tissue Aβ42 levels, bilateral hippocampus was homogenized in 5 M guanidine-HCl buffer with 50 mM Tris, pH 8.0. Samples were diluted 10-fold with cold PBS with 1× protease inhibitor cocktail and centrifuged at 16,000 g for 20 minutes at 4°C before application to the ELISA plate. Quantification of Aβ40 and Aβ42 in these fractions was performed using a Human Aβ ELISA kit (Invitrogen) according to the manufacturer's instructions.
Morris water maze and Y-maze tests. The mice were assessed in behavioral tests at 5 months of age. The Morris water maze test was performed as previously described (12). In the visible platform test performed on the first day, the mice were tested for 5 continuous trials, with an intertrial interval of 60 minutes. In the hidden platform tests on the following 4 days, a 10-cm-diameter platform was placed in the southwestern (SW) quadrant of the pool. Mice were trained for 5 trials, with an intertrial interval of 60 minutes. Each mouse was allowed 60 seconds to search for the platform. On the last day of the test, the platform was removed and the mice were given 60 seconds to locate where the platform was originally placed. Mouse movement was recorded by automated video tracking (ANY-maze, Stoelting). Y-maze testing occurred in a Y-shaped maze (Stoelting Co.) with 3 gray plastic arms at a 120° angle from each other, each 5 cm wide and 35 cm in length, and shielded with 10 cm high walls. The apparatus was positioned 60 cm above the floor. Each arm was differentiated using visual cues attached to the inner top portion of the wall: white (white squares), red (red triangles), or blue (blue dots). In this test, each mouse was placed in the center of the maze facing toward one of the arms and then allowed to explore freely for 8 minutes. Tracking of animals was achieved automatically using ANY-maze Video Tracking Software (version 4.63). The maze was cleaned before commencing with the next mouse. The spontaneous alternation was analyzed as a measure of spatial working memory. Percentage spontaneous alternation was calculated as the percentage of the number of alternations performed out of the total number of alternations possible, based on the number of arm entries per mouse.
Statistics. The number of times each experiment was repeated was at least 3 times as indicated in the figure legends. All results were expressed as mean ± SEM. Two-tailed Student's t test was used to analyze the difference between 2 groups, and multiple comparisons were analyzed by ANOVA followed by Bonferroni's multiple-comparisons test unless otherwise specified in figure legends. A significance of P < 0.05 was set for acceptance: *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001; NS, not significant.
Study approval. Human subject study was approved by the institutional ethics board of Xiangya Hospital. All subjects provided written informed consent for both study enrollment and sample collection. All animal experiments were conducted in accordance with institutional guidelines, with ethical approval from the UBC ACC (The University of British Columbia, Animal Care Center).