Use of Deep-Learning Genomics to Discriminate Healthy Individuals from Those with Alzheimer’s Disease or Mild Cognitive Impairment

Li, Lanlan; Yang, Yeying; Zhang, Qi; Wang, Jiao; Jiang, Jiehui; Neuroimaging Initiative, Alzheimer’s Disease

doi:https://doi.org/10.1155/2021/3359103

Behavioural Neurology

On this page

Abstract Introduction Materials and Methods Results Discussion Limitations Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Applications of Theranostics for Detecting and Targeting CNS Injuries and Diseases 2021

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 3359103 | https://doi.org/10.1155/2021/3359103

Use of Deep-Learning Genomics to Discriminate Healthy Individuals from Those with Alzheimer’s Disease or Mild Cognitive Impairment

Lanlan Li,¹Yeying Yang,²Qi Zhang,¹Jiao Wang,³Jiehui Jiang,¹and Alzheimer’s Disease Neuroimaging Initiative⁴

Academic Editor: Muh-Shi Lin

Received14 May 2021

Accepted11 Jun 2021

Published15 Jul 2021

Abstract

Objectives. Alzheimer’s disease (AD) is the most prevalent neurodegenerative disorder and the most common form of dementia in the elderly. Certain genes have been identified as important clinical risk factors for AD, and technological advances in genomic research, such as genome-wide association studies (GWAS), allow for analysis of polymorphisms and have been widely applied to studies of AD. However, shortcomings of GWAS include sensitivity to sample size and hereditary deletions, which result in low classification and predictive accuracy. Therefore, this paper proposes a novel deep-learning genomics approach and applies it to multitasking classification of AD progression, with the goal of identifying novel genetic biomarkers overlooked by traditional GWAS analysis. Methods. In this study, we selected genotype data from 1461 subjects enrolled in the Alzheimer’s Disease Neuroimaging Initiative, including 622 AD, 473 mild cognitive impairment (MCI), and 366 healthy control (HC) subjects. The proposed deep-learning genomics (DLG) approach consists of three steps: quality control, coding of single-nucleotide polymorphisms, and classification. The ResNet framework was used for the DLG model, and the results were compared with classifications by simple convolutional neural network structure. All data were randomly assigned to one training/validation group and one test group at a ratio of 9 : 1. And fivefold cross-validation was used. Results. We compared classification results from the DLG model to those from traditional GWAS analysis among the three groups. For the AD and HC groups, the accuracy, sensitivity, and specificity of classification were, respectively, , , and using the DLG model, while , , and using traditional GWAS. Similar results were obtained from the other two intergroup classifications. Conclusion. The DLG model can achieve higher accuracy and sensitivity when applied to progression of AD. More importantly, we discovered several novel genetic biomarkers of AD progression, including rs6311 and rs6313 in HTR2A, rs1354269 in NAV2, and rs690705 in RFC3. The roles of these novel loci in AD should be explored in future research.

1. Introduction

Alzheimer’s disease (AD) is the most common type of dementia and is an irreversible, progressive neurological brain disorder typically beginning with mild memory decline; in time, it can seriously impair an individual’s ability to carry out daily activities and lead to loss of autonomy [1, 2]. Mild cognitive impairment (MCI) is a preclinical stage of AD, in which individuals have no obvious cognitive behavioral symptoms but can show subtle prodromal signs of dementia [3, 4]. It is widely recognized that early detection of AD and MCI is essential to slowing progression.

Among factors that influence AD progression, common genetic variants are major risk factors [5]. Currently, the development of cheap comprehensive genetic testing of peripheral blood has brought dramatic changes to studies of the mechanisms of disease development. In recent decades, several genes have been associated with AD risk based on full-genome genotyping arrays using blood samples [6, 7]. For instance, genomics analysis showed APOE to be the most strongly associated AD risk gene [8]. In addition, the CLU, PICALM, SORL1, BIN1, and TOMM40 genes have also been identified as AD risk factors in the literature [7, 9, 10].

Technological advances [11] have allowed analysis of millions of nucleotide polymorphisms from thousands of subjects, including advanced genome-wide association studies (GWAS) and whole genome sequencing [12–16] that have increased our understanding of the genetic complexity of AD susceptibility. For instance, recent GWAS from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) have related known AD risk genes to differences in rates of brain atrophy and biomarkers of AD in the cerebrospinal fluid [17]. Moreover, the International Genomics of Alzheimer’s Project studied 74046 participants, confirming nearly all of the previous genetic risk factors and identifying 12 new susceptibility loci for AD [18]. Therefore, genomics analysis, especially GWAS analysis, has yielded important advances in AD research.

However, there are some limitations of GWAS. Firstly, traditional GWAS intergroup analysis is distorted by differences in sample sizes [19]. Secondly, traditional GWAS analysis is strongly dependent on prior knowledge and hand coding, which requires much time and energy and risks bias or errors in data entry [16] that can result in poor repeatability. Moreover, although traditional GWAS analysis can assure high specificity of disease screening, accuracy, and sensitivity are relatively low. In practice, false positives are preferred over false negatives in order to avoid omissions in disease screening. Therefore, alternative analytical tools would help to drive novel hypotheses and models.

Deep-learning algorithms implemented via deep neural networks can automatically embed computational features to yield end-to-end models that facilitate discovery of relevant highly complex features [20]. Seminal studies in 2015 demonstrated the applicability of deep neural networks to DNA sequence data [21, 22]. Deep convolutional neural networks (CNNs) have been used in recent studies to predict various molecular phenotypes on the basis of DNA sequence alone. Applications include classifying transcription factor binding sites, predicting molecular phenotypes such as DNA methylation, microRNA targets, and gene expression [23–27]. In addition, CNNs have been utilized to call genetic variants [28] and classify genetic mutations in tumors [29]. Multitask and multimodal models and transfer learning have also been developed in genomics [30, 31]. In this work, we hypothesize that deep-learning genomics (DLG) can be applied to AD and outperform traditional GWAS analysis. We propose a DLG method to replace traditional GWAS analysis for multitasking classification of AD progression and use this approach to seek novel genetic biomarkers of AD susceptibility.

2. Materials and Methods

The experimental workflow of this study consisted of three steps as shown in Figure 1. First, we conducted quality control and SNP genotype coding for SNP genotype data. Second, we used the deep residual network ResNet for DLG. The goal of the deep residual network was to obtain a model by supervised learning for prediction and extraction of DLG features. The details of this process are described in detail in the following sections. Finally, we investigated interpretability of the DLG model by applying Gradient-weighted Class Activation Mapping (Grad-CAM).

2.1. Subjects

Data used in the preparation of this study was obtained from the ADNI database (http://adni.loni.usc.edu/). ADNI was launched in 2004 by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, the Food and Drug Administration, private pharmaceutical companies, and nonprofit organizations, as a $60 million, 5-year public-private partnership. In this study, 1461 individuals (622 AD, 473 MCI, and 366 healthy controls (HCs)) from the ADNI 1, ADNI 2, and ADNI GO cohorts of the ADNI database were included. Meanwhile, the following data from the 1461 ADNI participants was downloaded: Illumina SNP genotyping data, demographic information, and diagnosis information. Written informed consent was obtained from all participants, and the study was conducted with prior institutional review board approval. Clinical characteristics, including age, sex, education, and Montreal Cognitive Assessment (MoCA) results, were collected and are listed in Table 1.

The subjects were of age 55-90 (inclusive) years. The detailed ADNI eligibility criteria are available from http://adni.loni.usc.edu/methods/documents/. In brief, eligibility criteria for these participants were as follows: (1) normal subjects: a Clinical Dementia Rating (CDR) of 0, nondepressed, non-MCI, and nondemented; (2) MCI subjects: a memory complaint, objective memory loss measured by education adjusted scores on Wechsler Memory Scale 7/Logical Memory II, a CDR of 0.5, absence of significant levels of impairment in other cognitive domains, essentially preserved activities of daily living, and an absence of dementia; (3) AD: CDR of 0.5 or 1.0 and met the National Institute of Neurological and Communicative Disorders and Stroke and Alzheimer’s disease and Related Disorders Association criteria for probable AD [32]. Specific psychoactive medications were excluded.

We investigated two groups of subjects using SNP genotype data collected from the ADNI databases. Our training and validation group contained 560 subjects with AD, 426 subjects with MCI, and 330 HC subjects. We used the SNP genotype data from this group to establish and test the validity of our predictive models. Our test group consisted of 62 AD subjects, 47 MCI subjects, and 36 HC controls, and we used the SNP genotype data to evaluate the diagnostic value of the predictive models.

2.2. DNA Isolation and SNP Genotyping

SNP genotyping for more than 620,000 target SNPs was completed on all ADNI participants using the following protocol. First, a total of 7 mL of blood was taken from each participant and stored in EDTA-containing Vacutainer tubes, and genomic DNA was extracted using the QIAamp DNA Blood Maxi Kit following the manufacturer’s protocol. Second, lymphoblastoid cell lines were established by transforming B lymphocytes with Epstein-Barr virus [33]. Fourteen genomic DNA samples were analyzed using the Human 610-Quad BeadChip according to the manufacturer’s protocols. Before starting the assay, a 50 ng sample of genomic DNA from each participant was examined qualitatively on a 1% Tris-acetate-EDTA agarose gel to check for degradation. Degraded DNA samples were excluded from further analysis. Third, samples were quantitated in triplicate with PicoGreen® reagent and diluted to 50 ng/L in TrisEDTA buffer (10 mM Tris, 1 mM EDTA, pH 8.0). A total of 200 ng of DNA was denatured, neutralized, and amplified for 22 hours at 37°C, and then fragmented with FMS reagent (Illumina) at 37°C for 1 hour, precipitated with 2-propanol, and incubated at 4°C for 30 minutes. Fourth, the resulting blue precipitate was resuspended in RA1 reagent (Illumina) at 48°C for 1 hour. Samples were then denatured (95°C for 20 minutes) and immediately hybridized onto BeadChips at 48°C for 20 hours. The BeadChips were washed and subjected to single base extension and staining. Finally, the BeadChips were coated with XC4 reagent (Illumina), desiccated, and imaged on a BeadArray Reader (Illumina). Illumina BeadStudio 3.2 software was used to generate SNP genotypes from bead intensity data.

2.3. Quality Control and APOE Genotype

The following quality control (QC) steps were performed on the 1461 samples using PLINK v1.07 software. QC processes were conducted separately between the AD and HC groups, the HC and MCI groups, and the AD and MCI groups. SNPs and participants were excluded from the analysis if they failed to meet any of the following criteria [34]: call rate per ; ; gender check; ; Hardy–Weinberg equilibrium test of ; . After the QC procedure, the numbers of features considered for future analysis of each subject in the paired groups were as follows: 301,388 in the HC and MCI groups, 301,853 in the HC and MCI groups, and 301,138 in the MCI and AD groups. The overall genotyping rate for the remaining dataset was over 99.5%.

In addition, although the APOE gene is an important target gene in AD research, it was not available for all identified APOE SNPs on the Illumina array. Therefore, based on the reported APOE ε2/ε3/ε4 status, the genotypes of the unavailable APOE SNPs were added manually to ADNI genotype data before assessing sample quality.

2.4. SNP Genotype Coding

A single-nucleotide polymorphism is a DNA sequence variation which occurs when a single nucleotide (A, T, C, or G) in the genome differs among members of a biological species or across paired chromosomes. Based on the satisfactory ADNI GWAS SNP data of this study, we encoded SNPs using the following coding scheme: 1 refers to A, 2 refers to T, 3 refers to C, and 4 refers to G.

2.5. GWAS Analysis

In the multitasking classification of this study, GWAS analysis, which has emerged as a popular tool for identifying genetic variants associated with disease risk, was designed to be compared with deep-learning models. Standard analysis of a case-control GWAS involves assessing the association between each individual genotyped SNP and disease risk. Manhattan and quantile–quantile (–) plots were used to visualize the GWAS results. All association results surviving the significance threshold of were saved and prepared for subsequent pattern analysis.

2.6. Deep-Learning Genomics Model Based on ResNet

The DLG model acted as a feature encoder, which had a significant impact on classification. In this study, we applied ResNet, a deep residual network, to the classification between AD and HC groups, AD and MCI groups, and HC and MCI groups. Residual units were added to the deep residual network on the basis of CNNs.

A CNN, the most effective type deep-learning model, is generally composed of three types of layers: convolutional, pooling, and fully connected. The following describes the operation of a CNN. The first step is to convolve the input sequences with a set of filter kernels; all the features active at different positions after convolution constitute the feature map [35]. A nonlinear activation function, typically a rectified linear unit (ReLU), is applied on each layer and on the sum of the feature maps. The operation of the convolutional layer and ReLU can be expressed as follows: where is the output of the convolutional layer, represents the number of filters in the layer, and are, respectively, the weight and bias of the filter of the layer, is the output of previous layer , and denotes the convolutional operation.

Next, the resulting feature map is processed through the pooling layer by taking either the mean or maximum activation over disjoint regions for each channel [20, 35]. By sequential combination of convolutional and pooling layers, a multilayer structure is built for feature description. Lastly, the fully connected layers are employed for classification. In total, when given a training set , the learning process of a CNN with convolutional layers, whose filter parameters are , the bias values are , and refers to classification layers, can be represented as an optimization learning task:

where is the loss function that represents the cost difference between the true label and the predictive label from the CNN model .

Based on the CNN model, the greatest advantage of the ResNet framework lies in adding identity mapping that is performed by the shortcut connections, the outputs of which are added to the outputs of the stacked layers [36]. Therefore, the ResNet addressed the degradation problem and added neither extra parameters nor computational complexity. The formula for residual learning was designed as follows: the desired underlying mapping is denoted as , and the stacked nonlinear layers were allowed to fit a separate mapping of . The original mapping was recast into . Thus, the overall representation of the residual block was as follows:

The formulation of can be realized by feedforward neural networks using “shortcut connections.” A deep residual network can be established by stacking a series of residual blocks. Specifically, there were two steps in the process: forward computation and backward propagation. When residual blocks are chosen to stack, the forward propagation of such a structure can be expressed by where and are the input and the output of the residual network, respectively, and is the weight related to the residual block, being the number of layers within the block.

Likewise, the back propagation of the overall loss of the neural network to can be denoted as where is the whole loss function of the neural network.

Before modeling using the above procedures, each subject’s SNP genotype data was cropped after quality control and mapped to . The pathology type was encoded to one-hot, which was the label. Thereafter, in the training stage, SNP genotype data was fed into the network to update model parameters via backward propagation with the Adam algorithm, a first-order gradient-based optimization algorithm which has been proven to be computationally efficient and appropriate for training deep neural networks. The outputs of the network were used as the classification results, and the crossentropy of the outputs was calculated as the loss function. More specifically, the output of the network for each individual SNP could be a binary value. 1 represented the highest probability of being AD subjects, while 0 represented highest probability of being HC subjects.

We adopted ResNet18 and ResNet34 frameworks in this study. Meanwhile, we also utilized a traditional CNN model for the comparative experiments of classification. In the ResNet models, we set learning rate into 1⁻³ and applied the Adam optimizer to update the model parameters with the batch size of 8. The maximum number of iterations was set into 20. Note that we used L2 regularization in this step to prevent the overfit of our model. For adjusting the CNN model parameters, we set learning rate into 1⁻² and applied the Adam optimizer to update the model parameters with the batch size of 8. The maximum number of iterations was set to 30. Above deep-learning models were processed on a GPU (graphics processing unit, GTX 1080 Ti acceleration of PyCharm 3.5).

For investigating the interpretability of the DLG model, the last convolutional layer of the last res-block was made transparent in order to extract DLG features by applying Grad-CAM and two-sample -tests. For the first step, the last convolutional layer of the last res-block was chosen to extract normalized DLG features. Subsequently, using a two-sample -test with a false discovery rate [37, 38], we compared the -coefficients of the AD and HC groups, the HC and MCI groups, and the MCI and AD groups.

2.7. Classification

In this study, the subjects of multitasking classification were randomly divided into one training group and one independent test group at a ratio of 9 : 1 as shown in Table 1. The training group was then used to optimize the model parameters. We also randomly chose 25% of training group to form a validation group to guide the choice of hyperparameters.

On the one hand, we conducted training of several deep-learning models, including ResNet18, ResNet34, and a traditional CNN, and compared classification performance in order to screen for the optimum DLG. On the other hand, in order to verify the diagnostic capabilities of the DLG model compared with traditional GWAS analysis, we also designed comparative trials. Among all the gene indicators, theta proved to be the most directly related to SNP changes. APOE ε4 status and the normalized theta-value of the significant SNP loci found in this study were seen to be genetic predictors, and we used the support vector machine (SVM) with the linear kernel 500 times for classification of traditional GWAS.

To evaluate classification performance, we repeatedly conducted 5-fold crossvalidation in the training group. Accuracy, sensitivity, and specificity were used to evaluate the results. The mathematical expression of the three parameters was as follows: where , , , and denote, respectively, true negatives, true positives, false negatives, and false positives.

A receiver-operating characteristic (ROC) curve was produced to intuitively compare the results of the different approaches, and the area under the curve (AUC) of the ROC was computed to quantitatively evaluate classification performance.

2.8. Statistical Analysis

Demographic characteristics were compared between groups using a two-sample -test or the chi-square test. In addition, a two-sample -test of the extracted features was applied as a criterion to estimate the differences in DLG features between AD patients and HCs, AD patients and MCIs, and HCs and MCIs. All statistical analyses were performed using SPSS Version 22.0 software (SPSS Inc., Chicago, IL) and Matlab2016b (Mathworks Inc., Sherborn, MA, United States). All values < 0.05 were considered significant.

3. Results

3.1. Outcomes of GWAS Analysis

We carried out case-control GWAS analysis between the AD and HC groups and observed two genome-wide significant loci on chromosome 19, including rs429358 (APOE, the epsilon 4 marker) and rs2075650 (TOMM40). Figures 2 and 3 show the resulting Manhattan and – plots, and Table 2 summaries the SNPs that achieved genome-wide significance. The value used to assess significant differences was calculated as , where indicates the number of satisfied SNPs.

3.2. Classification Performance

Table 3 lists the performance of the different multitasking classification methods, including classification accuracy, sensitivity, specificity, and AUC. Taking the result of classification between the AD and HC test group subjects as an example, accuracy, sensitivity, specificity, and AUC were, respectively, , , , and 0.744, for the GWAS analysis, , , , and 0.915 for the CNN model, , , , and 0.980 for ResNet18, and 98, , , and 0.981 for ResNet34. We found that the deep-learning model exhibited high accuracy, sensitivity, and specificity, whereas accuracy and sensitivity were low for the GWAS analysis. Therefore, we concluded that deep-learning models were superior to traditional GWAS analysis for classification. And compared with the CNN model, the results using ResNet were more robust and stable. These results were the same using the other two group-level classifications. Based on these results, ResNet34 was chosen for the DLG model because the observed classification performance was optimal among the several deep-learning models. A more intuitive comparison is provided by the ROC curves of the multitasking classification shown in Figure 4.

(a)

(b)

(c)

3.3. Interpretability of the DLG Model

Setting a threshold of , more than ten thousand SNP loci showed differences between the groups, and even the significance of the most frequently identified loci was below 0.001.

Firstly, we compared the significant SNPs with those previously identified by GWAS as genetic susceptibility factors. Almost one hundred SNP loci between AD patients and HCs were consistent with findings from previous studies. Likewise, more than one hundred associated SNP loci were also found between the AD and MCI groups and between the HC and MCI groups.

Secondly, we sought significant SNP loci among three classification tasks. The gene regions of sixty-six SNP loci were shared in different stages of AD progression. Table 4 summarizes the sixty-six shared significant SNP loci among the three classifications, including, e.g., the well-known CLU, PICALM, and SORL1 gene regions. For rs11136000 (CLU) in chromosome 8, its values were 6.63⁻⁴ between the AD and HC groups, 8.37⁻⁶ between the MCI and HC groups, and 1.49⁻⁷ between MCI and AD groups. In addition, three SNP loci, rs543293, rs10501602, and rs3851179, were found in the PICALM gene region of chromosome 11, and the values of rs3851179 for the comparisons of the three groups were 6.00⁻³, 1.06⁻⁶, and 1.51⁻²⁰, while the values of rs543293 were 4.65⁻³, 3.43⁻⁸, and 2.27⁻¹³, and the values of rs3851179 were also much less than 0.001. These results are well supported by previous studies. Other significant results are detailed in Table 4, and the heatmaps of significant SNPs in chromosomes 8, 11, and 13 are shown in Figure 5. The horizontal axis represents major and minor alleles, and the vertical axis represents the value of SNP loci in the chromosomes. We observed some distinct differences, for example, between rs11136000, rs3851179, and surrounding loci.

(a)

(b)

(c)

(d)

In addition, except for those in Table 4, there were also several SNP loci showing an association with AD progression in their respective classifications. Several also have been reported and confirmed in previous large-scale GWAS studies, including APOE, BIN1, CHRM1, and TOMM40 with values much less than 0.001. Furthermore, it is notable that rs6311 and rs6313 in the HTR2A gene region, rs1354269 in the NAV2 gene, and rs690705 in the RFC3 gene all exhibited significant differences among the three classifications. For instance, the values of rs6311 were 1.96⁻⁵, 2.52⁻³, and 1.48⁻¹¹ between the respective groups, and the values of rs6313 were 3.21⁻⁵, 4.55⁻³, and 2.05⁻¹². An understanding of the roles of these novel loci in AD requires future study.

All of the information above was deposited in the DisGeNET database (http://www.disgenet.org/home/), a discovery platform containing one of the largest publicly available collections of genes and variants associated with human disease.

4. Discussion

This study used a comparison of the performance of several different deep-learning models as a basis for proposing a deep-learning genomics method based on ResNet34. The classification results indicate that the DLG model offers a higher diagnostic value than traditional GWAS analysis.

4.1. Outcomes of GWAS Analysis

In GWAS analyses, two SNPs have been identified at the significance level: APOE SNP rs429358 was determined to be the most significant genetic risk factor for AD. And the second most significant factor, TOMM40 SNP rs2075650, was found to be adjacent to the APOE SNP [10]. These results are consistent with previous studies. Although these SNP loci were identified by GWAS, traditional GWAS analysis suffers from being influenced by small sample size. Because other common genetic risk factors may have a much smaller impact on risk than the APOE gene, novel risk factors present in small samples may go undetected by GWAS analysis. Several previous studies have also demonstrated an explicit relationship between sample size and the number of significant differences in traits identified by genome-wide association studies [18, 19].

4.2. Classification Performance

In this study, in order to construct a deep-learning genomics model, we compared the performance of several deep-learning classification methods, including a simple CNN model, ResNet18, and ResNet34. As shown in Table 3, we observed that the results of the deep residual network were superior to those of a simple CNN, and in the process of training the model, the ResNet models exhibited robustness and stability superior to those of CNNs, and furthermore, ResNet34 was superior to RseNet18. Therefore, we chose ResNet34 as the final DLG model. More importantly, we compared the performance of the DLG model and traditional GWAS analysis under the same conditions and found the classification results of the DLG model to be superior. These results suggest that the deep-learning algorithm is effective in genome applications and that development of deep-learning genomics is worthy of further exploration.

4.3. Interpretability of the DLG Model

When we interpreted the DLG model, we found more than one thousand SNP loci with significant differences between AD patients and HCs, between the MCI and AD groups, and between the HC and MCI subjects. As is well known, rs11136000 (CLU), rs3851179 (PICALM), rs2070045 (SORL1), and rs1699102 (SORL1) have previously been identified as risk factors for AD [7, 9, 39]. Notably, they were all included among the sixty-six significant SNP loci shared in the three classification tasks in this study (as shown in Table 4 and Figure 5). For example, previous studies have shown that CLU modulates Aβ metabolism and is involved in Aβ clearance or acts as a chaperon for protein degradation [40]. PICALM, as an adaptor protein involved in clathrin-mediated endocytosis, regulates amyloid precursor protein (APP) internalization and subsequent Aβ generation, contributing to brain amyloid plaque load via its effect on Aβ metabolism [41, 42]. In addition to the analysis of the above identical SNP loci found among the three classification tasks, several differential loci were identified among one or two classification tasks, which are also consistent with previous research. Rs10194375 (BIN1), a protein that may be associated with tau-mediated pathology was identified as being significant between the AD and HC groups and the AD and MCI groups. In addition, rs2075650 (TOMM40), rs405509 (APOE), and rs429358 (APOE) were identified as significant between the HC and MCI groups and the MCI and AD groups. In summary, the DLG model is able to identify differential genomics in multitasking classification.

Most importantly, in addition to those shown to associate with AD in the past, we found several new SNP loci, including rs6311 (HTR2A), rs6313 (HTR2A), rs1354269 (NAV2), rs1946518 (IL18), rs1799986 (LRP1), rs690705 (RFC3), and rs7943454 (LUZP2), whose values were highly significant(as shown in Table 4). Rs6311 and rs6313 are in the HTR2A gene region. The HTR2A gene in humans is located on chromosome 13 and consists of exons separated by only two introns and encodes one of the receptors for serotonin. According to previous publications, HTR2A has received much attention in many psychiatric disorders such as mood disorders, attention deficit hyperactivity disorder, anxiety disorders, and schizophrenia. On the one hand, some studies have shown that medications for mood disorders and related conditions work by blocking 5-HT2A and altering the function of certain brain circuits. And blocking 5-HTR2A also seems to improve the effects of some antidepressants [43]. On the other hand, the numbers of the postsynaptic receptor HTR2A are reduced in the neocortex, and it seems to be involved in memory via its role in cortical pyramidal cells. For example, in AD research, HTR2A receptor densities in the brains of AD subjects were found to be reduced compared with age-matched controls, and the researchers also found this reduction correlated with the rate of decline of cognitive scores [44]. Hence, since subjects with AD or mild cognitive impairment exhibit depression and anxiety to various degrees, it is worth exploring whether rs6311 and rs6313 of the HTR2A gene contribute to AD susceptibility. Another significant locus identified here was rs1354269 located in the NAV2 gene region. The NAV2 gene, which encodes a member of the neuron navigator gene family, is highly expressed in brain and is involved in the development of the nervous system. Hence, the role of the NAV2 gene in AD is also worthy of future investigation. In addition, rs690705 of the RFC3 gene region also exhibited a significant difference in group-level classifications, and its impact on AD should be examined in the future.

5. Limitations

It is worth noting some limitations of this study. Firstly, only gene sequences were used as inputs to the DLG classification. In the future work, we plan to combine gene sequences with clinical data and brain imaging [45] to facilitate investigation of the mechanisms of AD progression by deep-learning genomics and deep-learning radiomics approaches. Secondly, we only classified information from the ADNI dataset in this study, so the results could be strengthened by including other datasets such as the Chinese populations. Thirdly, the number of subjects represented in this study may be limiting. Lastly, although this study has demonstrated the feasibility of DLG approach, it will be important to further explore the interpretability of deep-learning genomics.

6. Conclusions

In conclusion, the current study suggests that the deep-learning genomics approach is effective for multitasking classification research on AD progression and outperforms traditional GWAS analysis. Moreover, the several novel SNP loci identified in the DLG approach including rs6311 and rs6313 in HTR2A, rs1354269 in NAV2, and rs690705 in RFC3 are worthy of further exploration to better understand the mechanisms of AD.

Data Availability

The datasets presented in this study were obtained from ADNI (http://adni.loni.usc.edu/).

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Lanlan Li and Yeying Yang contributed equally to this work.

Acknowledgments

Data collection and dissemination for this project were funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI): the National Institutes of Health (grant number U01 AG024904) and the Department of Defense (award number W81XWH-12-2-0012). ADNI is funded by the National Institute of Aging and the National Institute of Biomedical Imaging and Bioengineering as well as through generous contributions from the following organizations: AbbVie, Alzheimer’s Association, Alzheimer’s Drug Discovery Foundation, Araclon Biotech, BioClinica Inc., Biogen, Bristol-Myers Squibb Company, CereSpir Inc., Eisai Inc., Elan Pharmaceuticals Inc., Eli Lilly and Company, EuroImmun, F. Hoffmann-La Roche Ltd. and its affiliated company Genentech Inc., Fujirebio, , IXICO Ltd., Janssen Alzheimer Immunotherapy Research & Development LLC., Johnson & Johnson Pharmaceutical Research & Development LLC., Lumosity, Lundbeck, Merck & Co. Inc., Meso Scale Diagnostics LLC., NeuroRx Research, Neurotrack Technologies, Novartis Pharmaceuticals Corporation, Pfizer Inc., Piramal Imaging, Servier, Takeda Pharmaceutical Company, and Transition Therapeutics. The Canadian Institutes of Health Research are providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (http://www.fnih.org/). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Disease Cooperative Study at the University of California, San Diego, CA, USA. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California, CA, USA. This study was supported by grants received from the National Natural Science Foundation of China (grant numbers 82020108013); the 111 Project (grant number D20031); the Shanghai Municipal Science and Technology Major Project (grant number 2017SHZDZX01).

References

“2019 Alzheimer's disease facts and figures,” Alzheimer's & Dementia, vol. 15, no. 3, pp. 321–387, 2019.
View at: Publisher Site | Google Scholar
“2020 Alzheimer's disease facts and figures,” Alzheimer's & Dementia, vol. 16, no. 3, pp. 391–460, 2020.
View at: Publisher Site | Google Scholar
B. Dubois, H. Hampel, H. H. Feldman et al., “Preclinical Alzheimer's disease: definition, natural history, and diagnostic criteria,” Alzheimer's & Dementia, vol. 12, no. 3, pp. 292–323, 2016.
View at: Publisher Site | Google Scholar
Q. Y. Dong, T. R. Li, X. Y. Jiang, X. N. Wang, Y. Han, and J. H. Jiang, “Glucose metabolism in the right middle temporal gyrus could be a potential biomarker for subjective cognitive decline: a study of a Han population,” Alzheimer's Research & Therapy, vol. 13, no. 1, p. 74, 2021.
View at: Publisher Site | Google Scholar
Y. Huang and L. Mucke, “Alzheimer mechanisms and therapeutic strategies,” Cell, vol. 148, no. 6, pp. 1204–1222, 2012.
View at: Publisher Site | Google Scholar
C. M. Karch, C. Cruchaga, and A. M. Goate, “Alzheimer's disease genetics: from the bench to the clinic,” Neuron, vol. 83, no. 1, pp. 11–26, 2014.
View at: Publisher Site | Google Scholar
D. Harold, R. Abraham, P. Hollingworth et al., “Erratum: Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's disease,” Nature Genetics, vol. 41, no. 10, p. 1156, 2009.
View at: Publisher Site | Google Scholar
A. Montagne, D. A. Nation, A. P. Sagare et al., “APOE4 leads to blood -brain barrier dysfunction predicting cognitive decline,” Nature, vol. 581, no. 7806, pp. 71–76, 2020.
View at: Publisher Site | Google Scholar
B. N. Vardarajan, Y. Zhang, J. H. Lee et al., “Coding mutations inSORL1and Alzheimer disease,” Annals of Neurology, vol. 77, no. 2, pp. 215–227, 2015.
View at: Publisher Site | Google Scholar
H. Huang, J. Zhao, B. Xu et al., “The TOMM40 gene rs2075650 polymorphism contributes to Alzheimer's disease in Caucasian, and Asian populations,” Neuroscience Letters, vol. 628, pp. 142–146, 2016.
View at: Publisher Site | Google Scholar
N. R. Wray, M. E. Goddard, and P. M. Visscher, “Prediction of individual genetic risk of complex disease,” Current Opinion in Genetics & Development, vol. 18, no. 3, pp. 257–263, 2008.
View at: Publisher Site | Google Scholar
V. Escott-Price, C. Bellenguez, L. S. Wang et al., “Gene-wide analysis detects two new susceptibility genes for Alzheimer's disease,” PLoS One, vol. 9, no. 6, article e94661, 2014.
View at: Publisher Site | Google Scholar
V. Escott-Price, R. Sims, C. Bannister et al., “Common polygenic variation enhances risk prediction for Alzheimer's disease,” Brain, vol. 138, no. 12, pp. 3673–3684, 2015.
View at: Publisher Site | Google Scholar
the Alzheimer's Disease Neuroimaging Initiative, CHARGE consortium, EADI1 consortium et al., “Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease,” Nature Genetics, vol. 43, no. 5, pp. 429–435, 2011.
View at: Publisher Site | Google Scholar
Alzheimer’s Disease Sequencing Project, J. C. Bis, X. Jian et al., “Whole exome sequencing study identifies novel rare and common Alzheimer's-associated variants involved in immune response and transcriptional regulation,” Molecular Psychiatry, vol. 25, no. 8, pp. 1859–1875, 2020.
View at: Publisher Site | Google Scholar
Karolinska Schizophrenia Project (KaSP) consortium, T. Elvsåshagen, S. Bahrami et al., “The genetic architecture of human brainstem structures and their involvement in common brain disorders,” Nature Communications, vol. 11, no. 1, article 4016, 2020.
View at: Publisher Site | Google Scholar
S. Kim, S. Swaminathan, L. Shen et al., “Genome-wide association study of CSF biomarkers Abeta1-42, t-tau, and p-tau181p in the ADNI cohort,” Neurology, vol. 76, no. 1, pp. 69–79, 2011.
View at: Publisher Site | Google Scholar
European Alzheimer's Disease Initiative (EADI), Genetic and Environmental Risk in Alzheimer's Disease (GERAD), Alzheimer's Disease Genetic Consortium (ADGC) et al., “Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer's disease,” Nature Genetics, vol. 45, no. 12, pp. 1452–1458, 2013.
View at: Publisher Site | Google Scholar
R. D. Ball, “Designing a GWAS: power, sample size, and data structure,” Methods in Molecular Biology, vol. 1019, pp. 37–98, 2013.
View at: Publisher Site | Google Scholar
G. Eraslan, Ž. Avsec, J. Gagneur, and F. J. Theis, “Deep learning: new computational modelling techniques for genomics,” Nature Reviews Genetics, vol. 20, no. 7, pp. 389–403, 2019.
View at: Publisher Site | Google Scholar
B. Alipanahi, A. Delong, M. T. Weirauch, and B. J. Frey, “Predicting the sequence specificities of DNA- and RNA-binding proteins by deep learning,” Nature Biotechnology, vol. 33, no. 8, pp. 831–838, 2015.
View at: Publisher Site | Google Scholar
J. Zhou and O. G. Troyanskaya, “Predicting effects of noncoding variants with deep learning-based sequence model,” Nature Methods, vol. 12, no. 10, pp. 931–934, 2015.
View at: Publisher Site | Google Scholar
M. Wang, C. Tai, W. E, and L. Wei, “DeFine: deep convolutional neural networks accurately quantify intensities of transcription factor-DNA binding and facilitate evaluation of functional non-coding variants,” Nucleic Acids Research, vol. 46, no. 11, article e69, 2018.
View at: Publisher Site | Google Scholar
H. Zeng and D. K. Gifford, “Predicting the impact of non-coding variants on DNA methylation,” Nucleic Acids Research, vol. 45, no. 11, article e99, 2017.
View at: Publisher Site | Google Scholar
C. Angermueller, H. J. Lee, W. Reik, and O. Stegle, “DeepCpG: accurate prediction of single-cell DNA methylation states using deep learning,” Genome Biology, vol. 18, no. 1, article 67, 2017.
View at: Publisher Site | Google Scholar
S. Cheng, M. Guo, C. Wang, X. Liu, Y. Liu, and X. Wu, “MiRTDL: a deep learning approach for miRNA target prediction,” IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 13, no. 6, pp. 1161–1169, 2016.
View at: Publisher Site | Google Scholar
J. Zhou, C. L. Theesfeld, K. Yao, K. M. Chen, A. K. Wong, and O. G. Troyanskaya, “Deep learning sequence-based ab initio prediction of variant effects on expression and disease risk,” Nature Genetics, vol. 50, no. 8, pp. 1171–1179, 2018.
View at: Publisher Site | Google Scholar
R. Luo, F. J. Sedlazeck, T. W. Lam, and M. C. Schatz, “A multi-task convolutional deep neural network for variant calling in single molecule sequencing,” Nature Communications, vol. 10, no. 1, article 998, 2019.
View at: Publisher Site | Google Scholar
P. Chang, J. Grinband, B. D. Weinberg et al., “Deep-learning convolutional neural networks accurately classify genetic mutations in gliomas,” AJNR. American Journal of Neuroradiology, vol. 39, no. 7, pp. 1201–1207, 2018.
View at: Publisher Site | Google Scholar
D. R. Kelley, Y. A. Reshef, M. Bileschi, D. Belanger, C. Y. McLean, and J. Snoek, “Sequential regulatory activity prediction across chromosomes with convolutional neural networks,” Genome Research, vol. 28, no. 5, pp. 739–750, 2018.
View at: Publisher Site | Google Scholar
D. Quang and X. Xie, “FactorNet: a deep learning framework for predicting cell type specific transcription factor binding from nucleotide-resolution sequential data,” Methods, vol. 166, pp. 40–47, 2019.
View at: Publisher Site | Google Scholar
G. M. McKhann, D. S. Knopman, H. Chertkow et al., “The diagnosis of dementia due to Alzheimer's disease: recommendations from the National Institute on Aging-Alzheimer's Association workgroups on diagnostic guidelines for Alzheimer's disease,” Alzheimer's & Dementia, vol. 7, no. 3, pp. 263–269, 2011.
View at: Publisher Site | Google Scholar
H. Neitzel, “A routine method for the establishment of permanent growing lymphoblastoid cell lines,” Human Genetics, vol. 73, no. 4, pp. 320–326, 1986.
View at: Publisher Site | Google Scholar
L. Shen, S. Kim, S. L. Risacher et al., “Whole genome association study of brain-wide imaging phenotypes for identifying quantitative trait loci in MCI and AD: a study of the ADNI cohort,” NeuroImage, vol. 53, no. 3, pp. 1051–1063, 2010.
View at: Publisher Site | Google Scholar
T. Han, C. Liu, W. Yang, and D. Jiang, “Deep transfer network with joint distribution adaptation: a new intelligent fault diagnosis framework for industry application,” ISA Transactions, vol. 97, pp. 269–281, 2020.
View at: Publisher Site | Google Scholar
K. Chen, K. Chen, Q. Wang, Z. He, J. Hu, and J. He, “Short-term load forecasting with deep residual networks,” IEEE Transactions on Smart Grid, vol. 10, no. 4, pp. 3943–3952, 2019.
View at: Publisher Site | Google Scholar
D. Titov, J. Diehl-Schmid, K. Shi et al., “Metabolic connectivity for differential diagnosis of dementing disorders,” Journal of Cerebral Blood Flow & Metabolism, vol. 37, no. 1, pp. 252–262, 2017.
View at: Publisher Site | Google Scholar
Alzheimer’s Disease Neuroimaging Initiative, M. Wang, J. Jiang et al., “Individual brain metabolic connectome indicator based on Kullback-Leibler Divergence Similarity Estimation predicts progression from mild cognitive impairment to Alzheimer's dementia,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 47, no. 12, pp. 2753–2764, 2020.
View at: Publisher Site | Google Scholar
the European Alzheimer's Disease Initiative Investigators, J. C. Lambert, S. Heath et al., “Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease,” Nature Genetics, vol. 41, no. 10, pp. 1094–1099, 2009.
View at: Publisher Site | Google Scholar
R. Zhu, X. Liu, and Z. He, “Association between CLU gene rs11136000 polymorphism and Alzheimer's disease: an updated meta-analysis,” Neurological Sciences, vol. 39, no. 4, pp. 679–689, 2018.
View at: Publisher Site | Google Scholar
W. Xu, L. Tan, and J. T. Yu, “The role of PICALM in Alzheimer's disease,” Molecular Neurobiology, vol. 52, no. 1, pp. 399–413, 2015.
View at: Publisher Site | Google Scholar
Z. Wang, H. Lei, M. Zheng, Y. Li, Y. Cui, and F. Hao, “Meta-analysis of the association between Alzheimer disease and variants in GAB2, PICALM, and SORL1,” Molecular Neurobiology, vol. 53, no. 9, pp. 6501–6510, 2016.
View at: Publisher Site | Google Scholar
I. M. Badamasi, M. S. Lye, N. Ibrahim, and J. Stanslas, “Genetic endophenotypes for insomnia of major depressive disorder and treatment-induced insomnia,” Journal of Neural Transmission, vol. 126, no. 6, pp. 711–722, 2019.
View at: Publisher Site | Google Scholar
A. Drago and D. de Ronchi, “HTR2A gene variants and psychiatric disorders: a review of current literature and selection of SNPs for future studies,” Current Medicinal Chemistry, vol. 14, no. 19, pp. 2053–2069, 2007.
View at: Publisher Site | Google Scholar
X. Q. Wang, W. J. Huang, L. Su et al., “Neuroimaging advances regarding subjective cognitive decline in preclinical Alzheimer's disease,” Molecular Neurodegeneration, vol. 15, no. 1, p. 55, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Lanlan Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1592

Downloads

1493

Citations