ABSTRACT
Introduction AD is the most complex disorder leading to dementia worldwide. Despite the disease burden among Indians the mutation spectrum in our subcontinent is not well examined.
Methods To identify probable causal variants for AD, in a total of 29 clinically diagnosed, majorly familial and young onset AD cases from Eastern India, whole exome analysis in Illumina NovoSeq. 6000 sequencing platform was performed as per standard methodology. Next, the coding variants were prioritised based on the literature and our bioinformatic analyses for genotype-to-phenotype correlation.
Results A total of 25 missense variants and 4 nonsense variants in 17 genes among 23 AD cases were identified as probable damaging ones in our study cohort. Amongst all, ABCA7 account for the maximum number of pathogenic variants (5/29). A lowering in the age of onset was observed for mutation carriers and cases belonging to the posterior cortical atrophy (PCA) subgroup. Our further comparative analyses suggested that variants in APP metabolism pathway genes are more common in PCA, frontal AD, and young-onset multi-domain amnestic phenotypes than dysexecutive and typical AD.
Conclusion Our study suggests that whole exome sequencing among Indian AD patients holds the potential to identify the most common gene as well as clinically relevant new causal variants which may highlight new insight into disease mechanism through future research.
INTRODUCTION
Alzheimer’s disease (AD) is the most common progressive neurodegenerative disorder, with varying degrees of manifestation between individuals of the same ethnicity. Among Indians, the prevalence rate of dementia is 7.4% (around 8.8 million individuals) over the age of 60 years (Lee et al., 2023) of which AD alone or in combination with vascular pathology constitutes the majority. However, in contrast to European Ancestry, little is known about the Indian genetic architecture conferring AD susceptibility. The maximum number of reports are on the genetic association of AD with ApoE4 (Bharath et al., 2010, Shankarappa et al., 2017). However, an earlier Whole exome sequencing (WES) has been successful in the identification of pathogenic AD-associated novel variants among Indians (Syama et al., 2018). Therefore, considering complex disease mechanisms, little knowledge of AD genetics, and highly heritable polygenic AD susceptibility (h2=0.58-0.79) (Karlsson et al., 2022), the present study aimed to adopt a WES strategy to highlight the genes probably associated with Familial Alzheimer’s Diseases (FAD) from eastern India.
MATERIALS AND METHODS
Study subjects
Based on radiological findings and standard diagnostic criteria for AD by National Institute on Aging and Alzheimer’s Association (NIA-AA) Workgroup criteria (Jack et al., 2024) here a total of 29 AD subjects with either positive family history or younger age of onset were recruited by Cognitive Neurologists and Neuropsychologists from the Bangur Institute of Neurosciences and National Neurosciences Centre Calcutta, Kolkata, India. Among them, 22 had a younger age of onset (< 65 years), and 24 had a positive family history of AD or dementia.
Whole exome sequencing
To perform WES, genomic DNA was extracted using the QIAamp DNA Blood kit (QIAGEN, Germany) followed by targeted gene capture using a custom capture kit. Paired-end sequencing was performed with 2×100/2×150 chemistry on Illumina NovoSeq 6000 sequencing platform. Reads were assembled and aligned to reference sequences based on NCBI RefSeq transcripts and human genome build GRCh37/UCSC hg19.
In-silico analyses
The variation calling, annotation, and prioritization were performed to identify the variants of interest and interpreted in the context of a single clinically relevant transcript. All variants identified in known 36 genes for AD and other dementia-related genes/loci were manually evaluated based on documented literature, online database and bioinformatics prediction. The nonsynonymous variants that showed damaging effects by more than 60% of tools used (13/21) or more would be considered probable damaging and reported here.
Results
Demographic details of the study cohort
To understand the genetic basis of AD in the Eastern Indian AD population, we performed a pilot WES in a total of 29 AD patients with cognitive impairment for an average of more than 1.5 years at the time of presentation. The mean age at onset is 55.03±11.59 years and age ranges from 35 to 74 years. The male-to-female ratio is 0.6111. An except for 5 cases, all of them have a positive family history.
Identification of genetic variants
Among the plethora of variants identified, here we report a total of 25 missense variants and 4 nonsense variants with high pathogenicity in 17 genes among 23 AD cases representing ∼79.31% (23/29) cases as carriers.
It includes majorly five missense mutations in ABCA7 (distributed throughout the protein), three in SORL1 (a novel and two reported), two novel variants like Phe383Val and Leu258Ser in PSEN1 (in TM6, TM7) and the pathogenic Val717Ile in APP gene. Apart from these, dementia and PD-related genes also appeared to carry pathogenic variants among AD patients (Table 1a). The Ser15ProfsTer39 of EIF4G1, Trp218Ter of GBA, Phe7LeufsTer80 of UNC13B, and Lys484Ter of VPS13C variants account for the nonsense variant identified here.
In contrast, missense variants mentioned in Table 1b (Arg566Cys in EIF4G1, Pro455Leu & Ala2075Glu in ABCA7 and Arg202Gln in GBA genes) showed less damaging effects upon analyses.
Phenotypic variability and relation to genotype
To explain the genotype-to-phenotype relation, it is noteworthy to mention that sixteen patients had a monogenic mutation and 7 cases were digenic/double mutants for these genes.
Typical AD, Posterior Cortical Atrophy (PCA), and Frontal variant of AD appeared as common clinical features followed by young-onset multi-domain amnestic and dysexecutive syndrome (Table 1b) in our study cohort.
Our preliminary comparisons of the age of onset between the groups showed that it is a) lowest for PCA and highest for typical AD (P-value=0.0423 and 0.0001297 respectively) and b) lower among the mutation carriers than non-carriers (53±11.21 vs 62.83±10.40; P=0.08438) (Table 2a).
Our further inter-group analysis showed that the genetic variants in the APP metabolism (for PCA 78%, for frontal variant 75%, for young-onset multi-domain amnestic 67%, and 36% for typical AD) and non-APP metabolism pathway (100% for the dysexecutive syndrome, 64% for typical AD, 33% for young-onset multi-domain amnestic, 25% for frontal variant, and 22% for PCA) are present across phenotypes in different frequencies (Table 2b). Therefore, the categorization of genes linked with AD according to metabolic pathway suggests that alterations in genes of the APP metabolism pathway are likely to result in PCA, frontal AD and young-onset multi-domain amnestic, while, non-APP metabolism pathway involvement may show dysexecutive, and typical AD.
The APOE4 load is quite high here; 18 individuals are found with the APOE4 allele in heterozygous condition and two were homozygous for the same. The E4 allele is present irrespective of the different AD variants.
DISCUSSION
Therefore, in the present study, after screening 29 AD patients from Eastern India we have identified several probable pathogenic variants in a total of 17 genes majorly in ABCA7, SORL1, APP and PSEN1.
The London variant in the APP gene that we identified here too, is known as the mutational hotspot for autosomal dominant AD. An increased Aβ1–42/Aβ1–40 leading to more pathological protein accumulation is the underlying mechanism for it. On the other hand, PSEN1 variants located in its transmembrane domains, TM6 and TM7 respectively are novel and uncharacterised. However, since these variants share the same domains with the two major catalytic aspartic acids (Asp275 and Asp385) which are important for their endoproteolytic activity and γ-secretase activity, here we speculate for a deficiency in protein’s function as their patho-mechanism towards early-onset AD. Amongst all the genes, ABCA7 was identified to carry the maximum number of variants which are clustered in Extracellular loop1 and NBD motif followed by 2nd transmembrane helix. A variability in the age of onset for ABCA7 carriers was observed in our study. Several factors like the location of variants, APOE genotype, and the presence of additional variants in other AD or PD-related genes could be reasons for such variable clinical manifestations. Unlike a previous study examining ABCA7 variants in AD cases with memory and non-memory cognitive impairment, here we found severe memory impairment in all carriers (Campbell et al., 2022). In line with growing evidence from candidate and GWAS studies, our study has also identified SORL1 as an important gene conferring risk to AD in familial cases. Our intergroup analysis also highlighted the contribution of non-APP pathway genes across different AD varieties even in a very small sample size.
An earlier study from India on AD-exome analysis already identified variants in genes involved in APP metabolism like ours, but the profile is different concerning location in protein and frequency (Syama et al., 2018). However, considering these two studies, we can say APP-metabolism genes contribute significantly to familial AD among Indians. Moreover, our Indian data regarding the influence of APOE genotype on phenotypic heterogeneity is in accordance with global data.
Therefore, in conclusion, genes related to Aβ formation and clearance can be termed as major players in Indian familial AD cases and future exome analyses in larger cohorts have the potential to identify new variants in common genes which can provide new disease insight through model study.
Funding
Supported by grants from the Department of Science & Technology, Govt. of India, under the Cognitive Science Research Initiative Programme (DST/CSRI-P/2017/22), and Department of Biotechnology, Ministry of Science & Technology, Govt. of India (BT/NIDAN/01/05/2018).
Credit authorship contribution statement
DS and AB were responsible for the concept, study design, experimental work, data analysis and manuscript preparation.
AM, AB, GD, US, and TKB were responsible for the clinical diagnosis, concept, study design and manuscript preparation.
BB and SM were responsible for data acquisition and manuscript preparation. SPH and SG were responsible for manuscript preparation.
All authors read the draft, provided their inputs and agreed on the final version of the manuscript.
Conflict of Interest
There is no conflict of interest.
Ethics statement
All procedures performed in studies involving human participants were in accordance with the ethical standards of the IPGME&R, Kolkata, and National Neurosciences Centre Calcutta, Kolkata, India. The Ethics Committees (Institutional Ethics Committee, IPGME&R, Kolkata and National Neurosciences Centre Calcutta, Kolkata) of the abovementioned Institutes approved the study protocol. Informed consent was taken as per guidelines of the Indian Council of Medical Research, National Ethical Guidelines for Biomedical and Health Research involving human participants, India.
Informed consent
Informed consent from all the participants was received before clinical data and sample collection.
Data Availability Statement
The data described in this study are available from the corresponding author upon reasonable request.
Acknowledgement
The authors thank the patients who participated in the study.
Footnotes
↵* (Joint first author)
