Alzheimer's Disease Blood Biomarkers Associated With Neuroinflammation as Therapeutic Targets for Early Personalized Intervention

Abstract

The definitive diagnosis of Alzheimer's Disease (AD) without the need for neuropathological confirmation remains a challenge in AD research today, despite efforts to uncover the molecular and biological underpinnings of the disease process. Furthermore, the potential for therapeutic intervention is limited upon the onset of symptoms, providing motivation for studying and treating the AD precursor mild cognitive impairment (MCI), the prodromal stage of AD instead. Applying machine learning classification to transcriptomic data of MCI, AD, and cognitively normal (CN) control patients, we identified differentially expressed genes that serve as biomarkers for the characterization and classification of subjects into MCI or AD groups. Predictive models employing these biomarker genes exhibited good classification performances for CN, MCI, and AD, significantly above random chance. The PI3K-Akt, IL-17, JAK-STAT, TNF, and Ras signaling pathways were also enriched in these biomarker genes, indicating their diagnostic potential and pathophysiological roles in MCI and AD. These findings could aid in the recognition of MCI and AD risk in clinical settings, allow for the tracking of disease progression over time in individuals as part of a therapeutic approach, and provide possible personalized drug targets for early intervention of MCI and AD.

Introduction

Alzheimer's Disease (AD) is one of the most widely studied neurodegenerative disorders, and is associated with widespread brain atrophy and cognitive decline. It is clinically characterized by memory deficits, and patients develop progressive neuropsychiatric symptoms such as apathy, delusions, and agitation (1). The World Health Organization (WHO) estimates that AD is the most common form of dementia and that it stands to be the seventh leading cause of death in elderly worldwide. The neuropathology of AD is widely associated with the accumulation of amyloid-beta (Aβ) plaques and neurofibrillary tangles in the brain, and the definitive diagnosis of AD is only possible by post-mortem microscopic examination of brain tissues (2). While neuropathologic changes are well correlated with cognitive decline in AD (3), it remains unfeasible to examine brain tissue for AD diagnosis in clinical settings.

Diagnostic approaches for AD historically involve the interpretation of neuroimaging data, neuropsychological tests and laboratory tests (4). Although much progress has been made in the diagnosis of AD, the methodology of diagnosis relies heavily on clinicians' interpretations of laboratory results and neuropsychological tests for detecting cognitive deficits (5), with cerebrospinal fluid biomarkers such as phosphorylated tau only being recently included as considerations for clinical diagnosis (6). However, the reliability of AD diagnosis remains variable due to several confounders such as human interpretation, age and education. For instance, the Mini-Mental State Exam (MMSE) is widely used by clinicians for the screening of dementia by administering a 30-points questionnaire to evaluate a subject's orientation, recall, attention, language, and comprehension abilities. While the MMSE is a quick way of assessing subjects' cognitive decline, its use in monitoring the progression of AD is limited due to its low sensitivity to intermediate conditions such as mild cognitive impairment (MCI) (7). The complexity due to the different variants of AD also requires clinicians to evaluate a subject's condition by relying on a combination of other non-quantitative factors for diagnosis, which can include, but are not limited to, the medical history of patients, neuroimaging, and interviewing patients' kin (8). This complex combination of diagnostic factors often leads to variability in diagnosis by different clinicians, therefore a diagnostic approach that is based on quantitative measurements of the biological process of AD will be useful to provide more timely and accurate AD screening.

Recent advancements in genome sequencing technology have contributed to improved accessibility of large transcriptomic datasets, leading to greater opportunities for identifying biomarkers associated with complex and rare diseases (9). The quantification of biomarkers from gene expression profiles may be a possible diagnostic approach for AD that reduces reliance on clinicians' experience thereby reducing human error and eliminating subjectivity in diagnosis. Gene expression profiles from several studies have uncovered valuable patterns in AD patients, namely the presence of Aβ and hyperphosphorylated tau in the brain (10). Unfortunately, such approaches still face the limitation of using brain tissue from biopsies, which carries the risk of complications (11), and may not be translatable to all clinical settings. Alternatively, blood gene expression profiles may prove useful in AD screening with significant reduction of risks and greater tissue accessibility. Two large-scale blood gene expression datasets were conducted in recent years that aim to detect biomarkers for early diagnosis of AD: the Alzheimer's Disease Neuroimaging Initiative [ADNI; (12)] and AddNeuroMed (13, 14), which include subjects from North America and Europe respectively Both studies comprise similar protocols and data modalities, namely clinical and cognitive tests, blood transcriptomics, and neuroimaging, although the ADNI includes additional data from positron emission tomography imaging and from genetic studies on cerebrospinal fluid biomarkers of AD (15).

In this study, we identified gene expression profiles associated with inflammation, vascular dementia, MCI, AD, stroke and other cerebrovascular diseases in a Singaporean MCI and AD cohort. The focus on a largely inflammatory set of genes stems from earlier work suggesting that neuroinflammation in AD is closely linked to neurodegeneration, the severity of which can be observed as alterations in white matter hyperintensity through neuroimaging. We therefore aim to identify potential biomarkers that could predict neuroinflammatory changes before severe and irreversible neurodegeneration in AD, or even in MCI, so as to provide early intervention. Our findings were then validated using blood gene expression data from ADNI.

However, blood gene expression studies are limited by two main drawbacks. Firstly, the gene expression profiles are usually difficult to interpret due to the data being highly variable (16). Secondly, the large proportion of genes to a small proportion of subjects creates a statistical hurdle in the search for potential biomarkers. Technological advancements in recent years have provided various supervised and unsupervised models for the applications in features extraction of large expression datasets (17). Supervised machine learning models, such as Random Forests (RF), are particularly useful in such applications due to their flexibility in both classification and regression studies. RF also has the advantage of tackling the non-linear nature of gene expression profiles (18). In this study, we employed a multi-stage machine learning pipeline for the exploration of selecting differentially expressed genes (DEGs) in AD and MCI subjects. The analysis pipeline was built upon Boruta, a statistically driven machine learning method, which differentiated the important genes from noise after rigorous iterations of RF models (19). We showed that machine learning techniques may prove useful in discovering potential biomarkers for AD and MCI detection in the large data pool of blood gene expressions.

Materials and Methods

This section discusses the methodology used for processing of datasets as well as machine learning pipeline used to extract out the genes of interest and building a classification model. The overall framework is illustrated in a simplified diagram as shown in Figure 1.

Figure 1

Study Subjects and Datasets

This study was conducted on blood gene expression datasets consisting of one set of clinical samples recruited from the Memory clinic of the National Neuroscience Institute in Singapore (NNI), and the publicly available ADNI dataset (12). The distribution of subjects from each dataset is shown in Table 1.

In this investigation, the NNI dataset was used as the discovery cohort for the identification of differentially expressed genes as potential biomarkers for MCI and AD, and their respective biological pathways between CN, MCI, and AD subjects. The applicability of these genes to the identification of MCI and AD subjects was then evaluated on the ADNI dataset.

Clinical samples from NNI were collected over a period of 3 years, from 2013 to 2016. Informed consent was obtained for all subjects. Participants underwent clinical evaluation of psychological and cognitive performance, using the MMSE, the Montreal Cognitive Assessment (MoCA), and magnetic resonance imaging (MRI). Healthy control (CN) subjects were further required to have no cognitive complaints, no significant cognitive defects, and a clinical dementia rating (CDR) of 0. Diagnosis of MCI was based on the criteria of the National Institute on Aging–Alzheimer's Association (NIA-AA) Research Framework (20), the clinical presentation of cognitive symptoms and neuropsychological deficits without significant functional impairment, and a CDR score of 0.5. Subjects with AD were diagnosed using criteria from the National Institute of Neurological and Communicative Disorders and Stroke (NINCDS) and the Alzheimer's Disease and Related Disorders Association (ADRDA). According to the NINCDS-ADRDA criteria, patients with mild AD display cognitive symptoms and deficits according to neuropsychological evaluation, with significant functional impairment (4). Patients with mild AD were also defined as having a CDR score of 1. The NNI data therefore constitute 254 subjects (99 CN, 61 MCI, and 94 AD, listed in Table 1).

RNA Isolation and Quantification for NNI Data

Non-fasting venous blood was drawn via antecubital venipuncture. RNA was isolated from 1 ml of whole blood using the QIAamp RNA Blood Mini Kit (QIAGEN no. 52304). Total RNA from subjects was extracted from whole blood on the same day as blood collection. Isolated RNA was stored at −80°C. Complementary DNA (cDNA) was obtained by reverse transcription of 10 ng of total RNA (Fluidigm no. 100-6298). One hundred eighty two target genes were shortlisted from a literature search of genes associated with inflammatory response, vascular dementia, MCI, AD, stroke, and other cerebrovascular diseases (Supplementary Table 1). Primers for the gene targets were designed using Primer3 (21) or with reference to the OriGene database, and validated by real time quantitative polymerase chain reaction (qPCR) using RNA isolated from lipopolysaccharide-induced human lymphoblastoid cells. The 96.96 Dynamic Array Integrated Fluidic Circuits (IFC) for Gene Expression was used for quantifying RNA for each gene. The data was collected using the BioMark HD Image Capture System for further analysis.

Data Preparation and Pre-processing

Blood gene expression data from NNI were normalized against the housekeeping gene RP2 to obtain relative expression values of each gene. Replicates of samples were averaged to obtain a single gene expression value per gene. Genes with more than 25% missing values were removed from the analysis. For the remaining genes, missing value imputation was carried out by the adaptive LSimpute method using the missMethods package (v.0.4.0) in R (v.4.1.1), based on the least squares principle and correlations between genes and between arrays (22–24). This resulted in expression data from 176 genes for analysis in the NNI dataset (Supplementary Table 1).

Likewise, relative gene expression values for the ADNI dataset were obtained by normalizing all gene expression values, which had been obtained by microarray, by that of RP2 for each respective subject, and relative expression values for replicate samples were averaged to obtain a final relative gene expression value for each subject. A subset of the ADNI dataset containing 151 genes in common with NNI was used for analysis. For both datasets, uniform manifold approximation and projection (UMAP) plots were constructed and visualized using the umap package (v.0.2.7.0) and plotly graphing library in R respectively (25, 26).

To verify that the NNI data were not strongly influenced by age and gender of the subjects, preliminary random forest regression models were constructed for regression of age, and for classification of male vs. female, and for CN vs. MCI, CN vs. AD, and MCI vs. AD, using the randomForest package (v.4.7.1) in R (27). The top 30 genes with greatest variable importance, quantified by their influence on the prediction error for each model, were denoted as being either age-dependent, gender-dependent, or diagnosis-dependent variables. The age-dependent and gender-dependent variables were compared against diagnosis-dependent variables to determine the proportion of genes whose expression values are potentially confounded by age and gender.

In preparation for feature selection, the NNI data were split into training and testing datasets. Data splits were conducted using the caTools package in R (v.1.17), such that each split preserves the original relative distributions of CN, MCI, and AD subjects in the training and testing data. For the identification of differentially expressed genes, 80% of the NNI data were used for training, while the remaining 20% were reserved as testing data for internal validation. A 80/20 training/testing split was also conducted on the ADNI data in preparation for external validation.

The unequal distribution of CN, MCI, and AD subjects in each dataset presented potential risks of increased bias during feature selection and classification. This was particularly pertinent for the ADNI data used in external validation, where the class imbalance caused by the higher proportion of MCI patients and relatively lower number of AD patients could have led to deceptively high accuracy metrics if the data were employed directly. To overcome this, the Synthetic Minority Oversampling Technique (SMOTE) was applied on the training and testing data from both NNI and ADNI (28), employing functions from the smotefamily package in R (v.1.3.1). SMOTE oversampled underrepresented groups by repeatedly generating a new data point between a randomly selected real data point, and a randomly selected point among its k nearest neighbors. In order to balance our training and testing data, for each underrepresented group, new data points were generated by SMOTE to equal the number of data points of the majority group, such that each class contained the same number of data points at each feature selection or classification step.

Feature Selection for Differentially Expressed Genes

Differentially expressed genes (DEGs) were identified for each pairwise comparison, specifically CN vs. MCI, CN vs. AD, and MCI vs. AD. Boruta was employed to extract DEGs from the balanced training NNI data based on multiple iterations of the RF classifier (19). The RF classifier was selected for its quick performance and unbiasedness in classification through majority voting. The Boruta method selects for DEGs through a statistically rigorous approach by comparing the input features against shadow variables, derived from randomization of input variable values, and collating the variables that have greater feature importance than the best-performing shadow feature for subject classification at every run, denoted as a “hit” for that variable for that run. Over multiple runs, the number of hits for each variable makes up a binomial distribution, from which the list of important genes can be determined. To further eliminate variability, for each pairwise comparison, up to 1,000 runs were performed within each round of Boruta, and Boruta was itself conducted 10 times. Genes that were deemed important more than 50% of the time were defined as being differentially expressed between the two conditions.

Validation and Evaluation of Biomarkers

Validation was conducted by constructing RF classifier models using the DEGs from Boruta for each pairwise comparison, which were trained on SMOTE-balanced data derived from the 80% training data of each dataset, and tested on the reserved 20% testing data, which was also balanced by SMOTE. Apart from internal validation on NNI data, to further validate that the DEGs identified from the NNI data were generalisable as biomarkers for MCI and AD, external validation was carried out by constructing RF classifier models from the NNI DEGs that were also present in ADNI data, and evaluating the performance of these models on ADNI data.

In both internal and external validation, the highest accuracies from 10 iterations of evaluation consisting of 30 RF classifier models each, were recorded and compared against 50% accuracy for each pairwise comparison using one-tailed Student's t-tests.

Protein-protein association networks were constructed based on DEGs from Boruta using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database [v.11.5; (29)].

Gene Set Enrichment Analysis

The roles of the DEGs in a wider biological landscape were explored by annotating these genes with their Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (30), using the enrichR package (v.3.0) in R (31, 32). The KEGG pathways that overlap between the three pairwise comparisons were of particular interest as they could represent biological processes implicated during transitions from CN to MCI to AD.

Results

Heterogeneity in Alzheimer's Disease

Heterogeneity in AD presents difficulties for diagnosis and developing drug treatments (33). This heterogeneity was observed even at the transcriptomic level, as illustrated by UMAP plots of the NNI and ADNI datasets (Figure 2).

Figure 2

There was no clear separation between CN, MCI, and AD subjects within both the NNI and ADNI datasets (Figures 2A,B). This suggests that there is no clear cluster of genes that characterize MCI and AD to differentiate them from the CN subjects. This pattern was retained even after the NNI and ADNI datasets were split into training/testing sets and adjusted for class imbalance using SMOTE, demonstrating that the lack of clear clusters is not an effect of class imbalance and differences between each group of subjects cannot be easily resolved (Supplementary Figure 1).

The lack of clear defining delineations between groups of subjects further complicates diagnosis for these conditions, and demonstrates how there is no “one size fits all” approach for AD treatment, emphasizing the need for personalized treatments based on one's unique disease profile. Therefore, we employed the Boruta algorithm to identify differentially expressed genes and potential biomarkers within this complex dataset.

Differentially Expressed Genes Identified as Potential Biomarkers

According to random forest regression models for age-dependent variables and classification models for gender and pairwise diagnosis conditions, more than 82% of diagnosis-dependent variables were not age-dependent, and 80% of diagnosis-dependent variables were not gender-dependent, suggesting that most differentially expressed genes across diagnosis in this study reflect changes across disease rather than age and gender differences between subjects.

From feature selection on the NNI training data balanced by SMOTE, 17 DEGs were identified between CN and MCI subjects, 10 genes identified between CN and AD subjects, and 16 genes identified between MCI and AD subjects, from a total of 176 genes in the initial dataset (Table 2).

These potential biomarker genes were used to construct predictive RF models for classifying the remaining 20% data from the NNI dataset balanced by SMOTE, as well as subjects from ADNI. Using the list of DEGs as predictive features to classify subjects of the ADNI dataset, the mean highest accuracies obtained were 59.55, 55.96, and 56.65% for the pairwise comparisons CN vs. MCI, CN vs. AD, and MCI vs. AD, respectively (Table 3). The classification accuracies using the list of DEGs as biomarkers were significantly higher than random classification of 50% for all three pairwise comparisons. Similarly, RF models constructed using the list of DEGs evaluated on NNI data classified between subjects from each pairwise comparison, i.e., CN vs. MCI, CN vs. AD, and MCI vs. AD, with max accuracies significantly higher than random classification (Supplementary Table 2).

Several previously published studies aiming to identify blood biomarkers for AD focus on distinguishing between healthy and disease states, and demonstrate variations in DEGs between studies depending on the analysis methods employed (34, 35). In contrast, in addition to identifying DEGs associated only with AD, this investigation aims to identify genes that are differentially expressed in MCI, compared to CN and AD, specifically in the context of inflammatory genes. Nevertheless, we evaluated the performance of DEGs from Li et al. (34, 35) in pairwise classifications of CN vs. MCI, CN vs. AD, and MCI vs. AD, compared with the three sets of DEGs from this study. The DEGs identified in this study performed similarly to those from previous studies in classifying between CN and AD subjects (Supplementary Table 3). Our DEGs identified performed better in the classification of CN vs. MCI and MCI vs. AD subjects. Overall, pairwise classification using DEGs identified in this study is likely more representative of differential expression across CN, MCI, and AD, compared to genes in past studies that focus on AD classification.

Protein-protein association networks obtained by STRING demonstrated that most of the DEGs across the three pairwise comparisons have been shown or hypothesized to share protein functions in past studies, particularly for genes from CN vs. AD and MCI vs. AD (Figure 3). This indicates that DEGs identified by Boruta are functionally related, thus supporting the possibility that the biological regulation of protein networks involving these genes is affected during MCI and AD development.

Figure 3

Differentially Expressed Pathways in MCI and AD

In our analysis of the gene annotations of each set of genes obtained from pairwise Boruta analyses, we were interested in signaling pathways that were found in the intersection between the three pairwise comparisons. These overlapping regions indicate gene annotations that could be implicated during transition across CN, MCI, and AD development, and could present potential therapeutic targets. From the KEGG pathway annotations, we identified five such pathways, namely the PI3K-Akt, IL-17, JAK-STAT, TNF, and Ras signaling pathways (Figure 4).

Figure 4

The location of each pathway in the Venn diagram in Figure 4 indicates the stage at which the pathway is implicated in the transition between conditions. Pathways at the center of the Venn diagram, where all three pairwise comparisons overlap, likely undergo expression changes at the transcriptomic level throughout CN, MCI, and AD development. Pathways that overlap between CN vs. MCI and CN vs. AD are likely to undergo changes during MCI, which are retained into AD, while pathways that overlap between CN vs. MCI and MCI vs. AD likely undergo changes that are more detectable during MCI but not necessarily AD. Pathways that overlap between CN vs. AD and MCI vs. AD indicate changes that are likely to be more representative of AD, rather than MCI.

Discussion

Visualization of the blood gene expression profile using UMAP did not show distinctive clusters of AD, MCI and CN, illustrating the heterogeneous nature of AD. The visualization under UMAP also demonstrated that the underlying structure of detecting biomarkers in AD could not be resolved easily even when using non-linear approaches. Therefore, conventional classification approaches will likely be ineffective and a robust machine learning algorithm will be required for a more accurate biomarker identification. In this regard, RF was chosen for its versatile applications in classification and regression, and Boruta was employed to further increase robustness. In this study, we employed a machine learning pipeline for the identification of AD-related genes, which in turn highlighted five signaling pathways that could potentially be used for discriminating between stages of AD progression.

We recognize that the ambiguity of a patient's true condition, under the MCI classification, may influence the feature selection as not all MCI patients evolve to AD, and some individuals can even revert back to CN. This suggests that the disease progression from MCI to AD is non-linear, as previously suggested in literature (34). The clinical diagnosis of MCI, upon which our feature selection algorithm depends, is also challenging due to the lack of clear distinguishing symptoms. In addition, the worsening of MCI may develop into other forms of dementia besides AD. Nevertheless, the association between MCI and AD development provides motivation for studying both of them in conjunction. By considering MCI and AD as different conditions, we aim to identify not only potential blood biomarkers for AD diagnosis, but also biomarkers for MCI diagnosis, with the goal to provide early intervention for AD by detection of MCI subjects before they progress into AD. Early detection of MCI subjects and treatment may alleviate their symptoms in the early stages and reduce the probability of MCI patients progressing into AD.

The genes identified in this study are able to classify subjects as belonging to either CN, MCI, or AD groups, with greater than random accuracy, thus supporting their potential application in biomarker panels for diagnostic purposes in clinical settings. However, the classification accuracies observed in this study also suggest that MCI and AD are highly complex and heterogeneous. This could have been a result of the gene expression landscape between healthy and diseased groups being relatively similar overall (35, 36). To reduce variability resulting from class imbalance in our data during feature selection and evaluation, we performed SMOTE on both our training and testing data, and employed several iterations of Boruta, with majority voting to reduce the likelihood of false positives being included in our final set of DEGs. Further validation and investigation of these DEGs as clinically relevant biomarkers for MCI and AD will ideally be conducted in the future using a much larger dataset (than the current NNI data used in this study), that is also supported by common AD biomarkers such as amyloid and tau.

Despite the discovery cohort being relatively small, our classification models using DEGs identified from NNI data also classified subjects with greater than random accuracy when validated on external ADNI data, with comparable prediction accuracies across all three pairwise comparisons of CN vs. MCI, CN vs. AD, and MCI vs. AD. Therefore, our DEGs can be generalized to different cohorts of subjects from varying geographical regions and data collection methods.

The five pathways (PI3K-Akt, IL-17, JAK-STAT, TNF, and Ras signaling) identified from gene enrichment analysis of the DEGs represent potential changes in biochemical signaling over different stages of MCI and AD onset and development. It is likely that all of the five pathways undergo some dysregulation during the course of MCI development and AD progression. The finding of inflammatory pathways in the differential gene expression analysis is not particularly surprising, given that the initial set of NNI genes were also associated with several inflammatory processes. Despite this, the localization of the different pathways to different regions of the Venn diagram in Figure 4 indicates that at different stages of disease progression, different pathways could experience more detectable transcriptomic changes. This has applications in diagnostic practice and personalized medicine–although MCI is associated with increased systemic inflammation in general (37), greater dysregulation of specific pathways could provide indications of disease stages in each patient, and provide potential targets for that disease stage.

Additional challenges remain to be addressed regarding blood transcriptomes as biomarkers for disease diagnosis. Relating the dynamics of gene expression in the brain to blood can be confounded by changes in RNA availability rather than specific pathway dysregulation, such as the possibility that extensive neuronal death in highly inflammatory environments could lead to less detectable RNA changes at later AD stages (38). Therefore, it would be prudent to consider if the observed transcriptional changes detected were indeed associated with disease progression.

Each of the five pathways identified in this study has individually been implicated in cognitive impairment, particularly when they are altered in the central nervous system, even in neurological conditions not directly linked to AD (39–43). Should their expression levels in blood aid in identifying individuals with cognitive impairment accurately, dysregulation of these pathways would be possible biomarkers for MCI. An association between brain and blood levels of their expression, with corroboration from behavioral studies, would further support the utility of these pathways as tools for assessing patients' risks of developing MCI and AD onset, providing opportunities for therapeutic intervention.

The five inflammatory pathways observed from gene enrichment have each been associated with AD, albeit usually in the context of pathway dysregulation in the brain. Of note, IL-17A demonstrates a strong correlation with the pathogenesis of AD, and evidence of elevated IL-17A was observed in AD patients (44). Studies performed on rodents with Aβ-induced neurodegeneration demonstrated improved memory function after treatment with IL-17 antibodies (45). This provides strong evidence that the inhibition of IL-17 reduces the degenerative effects of Aβ in the glial cells, preventing further complications in AD patients.

Another pathway that holds potential in the treatment of AD is the TNF signaling pathway. Numerous studies have described the elevation of TNF-α in AD patients (46, 47) and the intervention of this pathway has been shown to alleviate brain pathology in rodent models (48). Additionally, the role of the TNF signaling pathway is supported by a recent genome-wide association study, which suggests that regulation of this pathway and its interactions with other signaling pathways could be implicated in AD development and neuropathology (49). Several drugs have been proposed to target different parts of the TNF signaling pathway and there is some reported clinical evidence suggesting that intervention of this pathway lowers AD neuropathology (50). Currently there are a few FDA approved drugs that target the IL-17A signaling and TNF signaling pathways. A few drug examples include Secukinumab, an antibody that selectively binds and neutralizes IL-17A, and Etanercept, a biological antagonist to TNF-α, that has been shown in a pilot study to improve cognitive function (51). This information may be applicable in the treatment of MCI patients, if early detection and diagnosis is conducted through routine blood tests, to prevent the potential development of AD.

Additionally, the JAK-STAT signaling and PI3K-Akt signaling pathways were found to be differentially expressed between AD, MCI and CN, serving as potential grounds for future research or drug targets. Although the mechanism behind how these pathways lead to development of AD is unclear, past investigations have found correlations between altered activity of these pathways and development of AD. The JAK-STAT signaling pathway was found to be activated in reactive astrocytes present in rodent models (52), although the mechanism of how JAK-STAT is involved in AD has not been studied extensively. The PI3K-Akt signaling pathways were found to be inhibited by Aβ, leading to neuronal death (53), but the mechanism of how this occurs is still not well understood. The involvement of these two pathways in AD pathophysiology could be explored in further research, with the possibility of them being therapeutic targets.

Ras/ERK signaling and its associated MAPK signaling pathway have also been investigated widely in contexts of disease due to their diverse regulatory roles in processes such as cell survival, migration, proliferation, and differentiation. In AD, ERK proteins are hypothesized to play a role in mediating Tau hyperphosphorylation and β-secretase expression which influences Aβ aggregation (54). Ras/ERK signaling has also been shown to be activated by Aβ, with the resultant aberrant signaling leading to neurodegeneration in AD (41). Therefore, early detection and targeted treatment for Ras signaling dysregulation could provide a means for preventing prolific Aβ aggregation and neurodegeneration in early AD.

Differential expression in these five pathways indicates that these pathways are dysregulated at different stages of disease development over MCI and AD, although it should be noted that the complex relationship between blood and brain expression makes interpretation of blood gene expression difficult, specifically in understanding how differential expression in blood is indicative of transcriptional patterns in the brain. The transcriptional signature of some gene modules in the brain, particularly those with roles in basic cellular processes such as gene expression regulation and infection, have been observed to be preserved in the blood of healthy subjects (55). In AD patients, especially at advanced stages of disease, strong blood-brain correlation for transcription has been observed, including brain-specific genes and inflammation-associated genes (56). However, blood-brain transcriptional correlation in MCI subjects remains to be established. While further investigation is needed before it can be concluded if these patterns of dysregulation definitively reflect neurophysiological changes in the brain, the findings of past human and animal studies suggest that the blood and brain transcriptome are correlated for some pathways, including inflammatory genes such as those studied here, especially in neurodegeneration (57).

Another challenge in identifying blood biomarkers of disease relates to population heterogeneity. Highly generalisable biomarkers are rarely available, due to heterogeneity between individuals, even within the healthy control population–the blood transcriptome is highly dynamic and inflammatory processes may vary in response to age, underlying disease, viral infection, or even seasonal changes (16).

This not only presents an opportunity for personalized medicine, but also highlights its importance in MCI and AD diagnosis. Since the expression of these inflammatory pathways is dynamic and highly variable even in healthy subjects, there is a risk of false positives in MCI and AD diagnosis if one only considers differences in expression levels of genes in the five pathways between individuals. Consequently, in evaluating MCI and AD risk of an individual, it may be beneficial to track longitudinal changes in expression of the five pathways for each subject, rather than making a diagnosis based on data from a single time point.

There are several previously published studies that aim to identify blood biomarkers of AD using blood gene expression data (14, 34–36, 58–60). A number of these studies employed publicly available data from AddNeuroMed for biomarker discovery and evaluation, although the different methods for feature selection and classification employed for each study resulted in varying differentially expressed features between AD and CN subjects (14, 34, 35, 59, 60). To overcome variations in individual differentially expressed genes, researchers may employ pathway-based classification models instead (14, 36), which we also include in our investigation. Through pathway analysis, there is more consensus between studies–for instance, the JAK-STAT pathway we identify in this investigation was also mentioned as a pathway of interest by Li et al. (35). Here, our study differs from previously published work in that we focus on biomarkers of MCI as indicators of AD risk, and on the pathways that are implicated across the development of MCI and AD, rather than just differential expression at each stage compared to CN.

The findings of this study provide potential blood biomarker genes and pathways for assessing MCI and AD risk in individuals. However, the complex molecular landscape of MCI and AD presents opportunities for further research, with the aim of making definitive diagnoses based on RNA extracted from blood. For instance, the relationship between blood RNA levels and RNA and protein levels in the brain could be further explored to determine how the differences observed in blood between conditions are linked to physiological changes in the brain during AD progression. Another aspect that could be considered is how the biomarkers identified in this study compare against blood biomarkers for other diseases, especially inflammatory conditions.

Conclusions

A machine learning pipeline employing Boruta effectively identifies DEGs in AD, which can classify subjects as either MCI or AD patients. These genes are enriched in five pathways–the PI3K-Akt, IL-17, JAK-STAT, TNF, and Ras signaling pathways–which are possibly dysregulated in MCI and AD. These pathways are also potential biomarkers that classify MCI and AD patients with reasonable accuracy. Our study demonstrates the potential to make use of these discovered biomarkers for early diagnosis of MCI and AD patients through routine blood testing, thereby providing biological insights toward early intervention through targeted drug treatment development for preventing deterioration into AD. These biomarker pathways may also provide targets for disease monitoring over time as an approach for personalized medicine as their transcriptional patterns are altered as MCI develops into AD.

Funding

This study was supported by the Ministry of Education (MOE) Tier 3 grant (MOE2017-T3-1-002) to EG.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1.

  Bruen PD McGeown WJ Shanks MF Venneri A . Neuroanatomical correlates of neuropsychiatric symptoms in Alzheimer's disease. Brain. (2008) 131:2455–63. 10.1093/brain/awn151

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2.

  DeTure MA Dickson DW . The neuropathological diagnosis of Alzheimer's disease. Mol Neurodegener. (2019) 14:32. 10.1186/s13024-019-0333-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  Hyman BT Phelps CH Beach TG Bigio EH Cairns NJ Carrillo MC et al . National Institute on Aging–Alzheimer's Association guidelines for the neuropathologic assessment of Alzheimer's disease. Alzheimers Dement J Alzheimers Assoc. (2012) 8:1–13. 10.1016/j.jalz.2011.10.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  McKhann G Drachman D Folstein M Katzman R Price D Stadlan EM . Clinical diagnosis of Alzheimer's disease: report of the NINCDS-ADRDA work group under the auspices of department of health and human services task force on Alzheimer's disease. Neurology. (1984) 34:939–44. 10.1212/WNL.34.7.939

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  Weintraub S Besser L Dodge HH Teylan M Ferris S Goldstein FC et al . Version 3 of the Alzheimer Disease Centers' Neuropsychological Test Battery in the Uniform Data Set (UDS). Alzheimer Dis Assoc Disord. (2018) 32:10–7. 10.1097/WAD.0000000000000223

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  Porsteinsson AP Isaacson RS Knox S Sabbagh MN Rubino I . Diagnosis of early Alzheimer's disease: clinical practice in 2021. J Prev Alzheimers Dis. (2021) 8:371–86. 10.14283/jpad.2021.23

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  Arevalo-Rodriguez I Smailagic N Roqué i Figuls M Ciapponi A Sanchez-Perez E Giannakou A et al . Mini-Mental State Examination (MMSE) for the detection of Alzheimer's disease and other dementias in people with mild cognitive impairment (MCI). Cochrane Database Syst Rev. (2015) 2015:CD010783. 10.1002/14651858.CD010783.pub2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8.

  Dubois B Feldman HH Jacova C DeKosky ST Barberger-Gateau P Cummings J et al . Research criteria for the diagnosis of Alzheimer's disease: revising the NINCDS–ADRDA criteria. Lancet Neurol. (2007) 6:734–46. 10.1016/S1474-4422(07)70178-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  Fertig EJ Slebos R Chung CH . Application of genomic and proteomic technologies in biomarker discovery. Am Soc Clin Oncol Educ Book. (2012) 2012:377–82. 10.14694/EdBook_AM.2012.32.156

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  Yates D . Frontal cortex biopsy samples can predict Alzheimer disease. Nat Rev Neurol. (2011) 7:5–5. 10.1038/nrneurol.2010.190

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  Schott JM Reiniger L Thom M Holton JL Grieve J Brandner S et al . Brain biopsy in dementia: clinical indications and diagnostic approach. Acta Neuropathol. (2010) 120:327–41. 10.1007/s00401-010-0721-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  Petersen RC Aisen PS Beckett LA Donohue MC Gamst AC Harvey DJ et al . Alzheimer's Disease Neuroimaging Initiative (ADNI). Neurology. (2010) 74:201–9. 10.1212/WNL.0b013e3181cb3e25

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  Lovestone S Francis P Kloszewska I Mecocci P Simmons A Soininen H et al . AddNeuroMed-the European collaboration for the discovery of novel biomarkers for Alzheimer's disease. Ann N Y Acad Sci. (2009) 1180:36–46. 10.1111/j.1749-6632.2009.05064.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  Voyle N Keohane A Newhouse S Lunnon K Johnston C Soininen H et al . a pathway based classification method for analyzing gene expression for Alzheimer's disease diagnosis. J Alzheimers Dis. (2016) 49:659–69. 10.3233/JAD-150440

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  Kim S Swaminathan S Shen L Risacher SL Nho K Foroud T et al . Genome-wide association study of CSF biomarkers Aβ1-42, t-tau, and p-tau181p in the ADNI cohort. Neurology. (2011) 76:69–79. 10.1212/WNL.0b013e318204a397

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  Schmidt M Hopp L Arakelyan A Kirsten H Engel C Wirkner K et al . The human blood transcriptome in a large population cohort and its relation to aging and health. Front Big Data. (2020) 3:36. 10.3389/fdata.2020.548873

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  McDermott JE Wang J Mitchell H Webb-Robertson B-J Hafen R Ramey J et al . Challenges in biomarker discovery: combining expert insights with statistical analysis of complex omics data. Expert Opin Med Diagn. (2013) 7:37–51. 10.1517/17530059.2012.718329

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  Qi Y . Random forest for bioinformatics. Ensemble Mach Learn. (2012) 7:307–23. 10.1007/978-1-4419-9326-7_11

  • CrossRef
  • Google Scholar

• 19.

  Kursa MB Rudnicki WR . Feature selection with the boruta package. J Stat Softw. (2010) 36. 10.18637/jss.v036.i11

  • CrossRef
  • Google Scholar

• 20.

  Petersen RC Wiste HJ Weigand SD Fields JA Geda YE Graff-Radford J et al . NIA-AA Alzheimer's disease framework: clinical characterization of stages. Ann Neurol. (2021) 89:1145–56. 10.1002/ana.26071

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  Untergasser A Cutcutache I Koressaar T Ye J Faircloth BC Remm M et al . Primer3—new capabilities and interfaces. Nucleic Acids Res. (2012) 40:e115. 10.1093/nar/gks596

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  Bø TH Dysvik B Jonassen I . LSimpute: accurate estimation of missing values in microarray data with least squares methods. Nucleic Acids Res. (2004) 32:e34. 10.1093/nar/gnh026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23.

  R Core Team . R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing (2021).

  • Google Scholar

• 24.

  Rockel T . Missmethods: Methods for MissingData (2020).

  • Google Scholar

• 25.

  Konopka T . umap: Uniform Manifold Approximation andProjection (2020).

  • Google Scholar

• 26.

  Sievert C . Interactive Web-Based Data Visualization with R, plotly, and shiny. Chapman and Hall/CRC (2020). 10.1201/9780429447273

  • CrossRef
  • Google Scholar

• 27.

  Liaw A Wiener M . Classification and Regression by Randomforest2, 5 (2002).

  • Google Scholar

• 28.

  Chawla NV Bowyer KW Hall LO Kegelmeyer WP . SMOTE: synthetic minority over-sampling technique. J Artif Intell Res. (2002) 16:321–57. 10.1613/jair.953

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29.

  Szklarczyk D Gable AL Lyon D Junge A Wyder S Huerta-Cepas J et al . STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. (2019) 47:D607–13. 10.1093/nar/gky1131

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  Kanehisa M Goto S . KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. 10.1093/nar/28.1.27

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31.

  Chen EY Tan CM Kou Y Duan Q Wang Z Meirelles GV et al . Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. (2013) 14:128. 10.1186/1471-2105-14-128

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32.

  Jawaid W . enrichR: Provides an R Interface to “Enrichr.” (2021).

  • Google Scholar

• 33.

  Devi G Scheltens P . Heterogeneity of Alzheimer's disease: consequence for drug trials?Alzheimers Res. Ther. (2018) 10:122. 10.1186/s13195-018-0455-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  Li H Hong G Lin M Shi Y Wang L Jiang F et al . Identification of molecular alterations in leukocytes from gene expression profiles of peripheral whole blood of Alzheimer's disease. Sci Rep. (2017) 7:14027. 10.1038/s41598-017-13700-w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  Li X Wang H Long J Pan G He T Anichtchik O et al . Systematic analysis and biomarker study for Alzheimer's disease. Sci Rep. (2018) 8:17394. 10.1038/s41598-018-35789-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  Lee T Lee H . Prediction of Alzheimer's disease using blood gene expression data. Sci Rep. (2020) 10:3485. 10.1038/s41598-020-60595-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  Trollor JN Smith E Baune BT Kochan NA Campbell L Samaras K et al . Systemic inflammation is associated with MCI and its subtypes: the sydney memory and aging study. Dement Geriatr Cogn Disord. (2010) 30:569–78. 10.1159/000322092

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  Christov A Ottman JT Grammas P . Vascular inflammatory, oxidative and protease-based processes: implications for neuronal cell death in Alzheimer's disease. Neurol Res. (2004) 26:540–6. 10.1179/016164104225016218

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39.

  Cipollini V Anrather J Orzi F Iadecola C . Th17 and cognitive impairment: possible mechanisms of action. Front Neuroanat. (2019) 13:95. 10.3389/fnana.2019.00095

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  Ding S Hu J Yang J Liu L Huang W Gu X et al . The inactivation of JAK2/STAT3 signaling and desensitization of M1 mAChR in minimal hepatic encephalopathy (MHE) and the protection of naringin against MHE. Cell Physiol Biochem. (2014) 34:1933–50. 10.1159/000366391

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41.

  Kirouac L Rajic AJ Cribbs DH Padmanabhan J . Activation of Ras-ERK signaling and GSK-3 by amyloid precursor protein and amyloid beta facilitates neurodegeneration in Alzheimer's disease. eNeuro. (2017) 4:e0149-16.2017. 10.1523/ENEURO.0149-16.2017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42.

  Pillai JA Maxwell S Bena J Bekris LM Rao SM Chance M et al . Key inflammatory pathway activations in the MCI stage of Alzheimer's disease. Ann Clin Transl Neurol. (2019) 6:1248–62. 10.1002/acn3.50827

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  Shu Y Zhang H Kang T Zhang J Yang Y Liu H et al . PI3K/Akt signal pathway involved in the cognitive impairment caused by chronic cerebral hypoperfusion in rats. PLoS ONE. (2013) 8:e81901. 10.1371/journal.pone.0081901

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44.

  Chen J Liu X Zhong Y . Interleukin-17A: the key cytokine in neurodegenerative diseases. Front Aging Neurosci. (2020) 12:307. 10.3389/fnagi.2020.566922

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  Cristiano C Volpicelli F Lippiello P Buono B Raucci F Piccolo M et al . Neutralization of IL-17 rescues amyloid-β-induced neuroinflammation and memory impairment. Br J Pharmacol. (2019) 176:3544–57. 10.1111/bph.14586

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46.

  Brosseron F Krauthausen M Kummer M Heneka MT . Body fluid cytokine levels in mild cognitive impairment and Alzheimer's disease: a comparative overview. Mol Neurobiol. (2014) 50:534–44. 10.1007/s12035-014-8657-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47.

  Tarkowski E Liljeroth A-M Minthon L Tarkowski A Wallin A Blennow K . Cerebral pattern of pro- and anti-inflammatory cytokines in dementias. Brain Res Bull. (2003) 61:255–60. 10.1016/S0361-9230(03)00088-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48.

  He P Zhong Z Lindholm K Berning L Lee W Lemere C et al . Deletion of tumor necrosis factor death receptor inhibits amyloid β generation and prevents learning and memory deficits in Alzheimer's mice. J Cell Biol. (2007) 178:829–41. 10.1083/jcb.200705042

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49.

  Bellenguez C Küçükali F Jansen IE Kleineidam L Moreno-Grau S Amin N et al . New insights into the genetic etiology of Alzheimer's disease and related dementias. Nat Genet. (2022) 54:412–36. 10.1038/s41588-022-01024-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50.

  Decourt B Lahiri DK Sabbagh MN . Targeting tumor necrosis factor alpha for Alzheimer's disease. Curr Alzheimer Res. (2017) 14:412–25. 10.2174/1567205013666160930110551

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51.

  Tobinick EL Gross H . Rapid cognitive improvement in Alzheimer's disease following perispinal etanercept administration. J Neuroinflammation. (2008) 5:2. 10.1186/1742-2094-5-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52.

  Haim LB Ceyzériat K Sauvage MAC Aubry F Auregan G Guillermier M et al . The JAK/STAT3 pathway is a common inducer of astrocyte reactivity in alzheimer's and huntington's diseases. J Neurosci. (2015) 35:2817–29. 10.1523/JNEUROSCI.3516-14.2015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53.

  Gabbouj S Ryhänen S Marttinen M Wittrahm R Takalo M Kemppainen S et al . Altered insulin signaling in Alzheimer's disease brain – special emphasis on PI3K-Akt pathway. Front Neurosci. (2019) 13:629. 10.3389/fnins.2019.00629

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54.

  Kim EK Choi E-J . Pathological roles of MAPK signaling pathways in human diseases. Biochim Biophys Acta. (2010) 1802:396–405. 10.1016/j.bbadis.2009.12.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55.

  Cai C Langfelder P Fuller TF Oldham MC Luo R van den Berg LH et al . Is human blood a good surrogate for brain tissue in transcriptional studies?BMC Genomics. (2010) 11:589. 10.1186/1471-2164-11-589

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56.

  Naughton BJ Duncan FJ Murrey DA Meadows AS Newsom DE Stoicea N et al . Blood genome-wide transcriptional profiles reflect broad molecular impairments and strong blood-brain links in Alzheimer's disease. J Alzheimers Dis JAD. (2015) 43:93–108. 10.3233/JAD-140606

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57.

  Tylee DS Kawaguchi DM Glatt SJ . On the outside. looking in: a review and evaluation of the comparability of blood and brain “-omes”. Am J Med Genet B Neuropsychiatr Genet. (2013) 162:595–603. 10.1002/ajmg.b.32150

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58.

  Booij BB Lindahl T Wetterberg P Skaane NV Sæbø S Feten G et al . A Gene Expression pattern in blood for the early detection of Alzheimer's disease. J Alzheimers Dis. (2011) 23:109–19. 10.3233/JAD-2010-101518

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59.

  Lunnon K Sattlecker M Furney SJ Coppola G Simmons A Proitsi P et al . A Blood gene expression marker of early Alzheimer's disease. J Alzheimers Dis. (2013) 33:737–53. 10.3233/JAD-2012-121363

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60.

  Sood S Gallagher IJ Lunnon K Rullman E Keohane A Crossland H et al . A novel multi-tissue RNA diagnostic of healthy ageing relates to cognitive health status. Genome Biol. (2015) 16:185. 10.1186/s13059-015-0750-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar