Abstract
Short-term memory dysfunction is a key early feature of Alzheimer’s disease (AD). Psychiatric patients may be at higher risk for memory dysfunction and subsequent AD due to the negative effects of stress and depression on the brain. We carried out longitudinal within-subject studies in male and female psychiatric patients to discover blood gene expression biomarkers that track short term memory as measured by the retention measure in the Hopkins Verbal Learning Test. These biomarkers were subsequently prioritized with a convergent functional genomics approach using previous evidence in the field implicating them in AD. The top candidate biomarkers were then tested in an independent cohort for ability to predict state short-term memory, and trait future positive neuropsychological testing for cognitive impairment. The best overall evidence was for a series of new, as well as some previously known genes, which are now newly shown to have functional evidence in humans as blood biomarkers: RAB7A, NPC2, TGFB1, GAP43, ARSB, PER1, GUSB, and MAPT. Additional top blood biomarkers include GSK3B, PTGS2, APOE, BACE1, PSEN1, and TREM2, well known genes implicated in AD by previous brain and genetic studies, in humans and animal models, which serve as reassuring de facto positive controls for our whole-genome gene expression discovery approach. Biological pathway analyses implicate LXR/RXR activation, neuroinflammation, atherosclerosis signaling, and amyloid processing. Co-directionality of expression data provide new mechanistic insights that are consistent with a compensatory/scarring scenario for brain pathological changes. A majority of top biomarkers also have evidence for involvement in other psychiatric disorders, particularly stress, providing a molecular basis for clinical co-morbidity and for stress as an early precipitant/risk factor. Some of them are modulated by existing drugs, such as antidepressants, lithium and omega-3 fatty acids. Other drug and nutraceutical leads were identified through bioinformatic drug repurposing analyses (such as pioglitazone, levonorgestrel, salsolidine, ginkgolide A, and icariin). Our work contributes to the overall pathophysiological understanding of memory disorders and AD. It also opens new avenues for precision medicine- diagnostics (assement of risk) as well as early treatment (pharmacogenomically informed, personalized, and preventive).
Similar content being viewed by others
Introduction
“If you find ways to repair the memory damaged by Alzheimer’s disease or dementia and so forth, it is very likely that the same methods could be used to upgrade the memory of completely healthy people”
- Yuval Noah Harari
Alzheimer disease (AD) is a clear and present danger to older adults, and has a profound socio-economic impact. Existing therapies are limited in efficacy. Early identification of individuals at risk may open the door to preventive approaches. Short-term memory dysfunction is a key early feature of AD. We proposed to identify blood biomarkers that track a relevant related quantitative phenotype for short-term memory, the retention measure of recall in Hopkins Verbal Learning Test (HVLT) (Fig. S1). Blood biomarkers are easily accessible clinically, as opposed to brain imaging or even CSF changes. The nervous and immune system have some common developmental roots, and there are bi-directional brain-immune interactions [1]. Similar gene expression changes in brain and blood cells may also be due to similar gene promoter regulation in response to shared internal milieu and external environment changes [2,3,4,5,6]. Previous work by our group has identified blood gene expression biomarkers that track suicidal ideation using longitudinal within-subject designs, validated them in suicide completers, and tested them in independent cohorts for ability to assess state (suicidal ideation), and ability to predict trait (future hospitalizations for suicidality) [7,8,9,10]. Those biomarkers were also useful in pharmacogenomic and drug repurposing analyses. We recently successfully used a similar approach for pain [11], and for stress [12]. In the current study, we endeavored to use this comprehensive step-wise approach to identify biomarkers for short-term memory that may have relevance to AD. We conducted our longitudinal studies in psychiatric disorder subjects, a population enriched in memory retention abnormalities, and which may be at increased risk of AD and other aging-related disorders, due at least in part to the effects of stress and depression [13, 14]. The subjects had blood gene expression data at multiple testing visits, and were phenotyped at each visit, including with HVLT. They also had electronic medical records for long term follow-up of subsequent outcomes, including future neuropsychological testing as part of standard clinical care.
Materials and methods
Cohorts
We used two independent cohorts of psychiatric disorders patients, one for Discovery of candidate biomarkers, and one for validation/testing of top biomarkers (for predicting memory state, and predicting future positive neuropsychological testing for cognitive impairment) (Fig. 1, Table 1, and Table S1).
Similar to our previous studies [7,8,9], the psychiatric subjects are part of a larger longitudinal cohort of adults that we are continuously collecting. Subjects were recruited from the patient population at the Indianapolis VA Medical Center. All subjects understood and signed informed consent forms detailing the research goals, procedure, caveats and safeguards, per IRB approved protocol. Subjects completed diagnostic assessments by an extensive structured clinical interview—Diagnostic Interview for Genetic Studies, and up to six testing visits, 3–6 months apart or whenever a new psychiatric hospitalization occurred. At each testing visit, they received a series of rating scales, including a Hopkins Verbal Learning Test (HVLT-R), and the blood was drawn. We collected whole blood (10 ml) in two RNA-stabilizing PAXgene tubes, labeled with an anonymized ID number, and stored at −80 °C in a locked freezer until the time of future processing. Whole-blood RNA was extracted for microarray gene expression studies from the PAXgene tubes, as detailed below.
For this study, our within-subject longitudinal discovery cohort, from which the biomarker data were derived, consisted of 159 subjects (131 males, 28 females) with multiple testing visits (a total of 496), who each had at least one 20% change in the Retention measure of HVLT (Fig. S1) from one consecutive testing visit to another.
Our independent test cohort for predicting state (low memory retention) consisted of 127 subjects (97 males, 30 females), demographically matched with the discovery cohort, with one or more testing visits in our lab (for a total of 238 visits). Low Memory Retention was defined as a score of ≤ 40 (Fig. 1, Table 1 and Table S1). Subjects were matched for age between low memory group and the rest group.
Our test cohort for predicting trait (future positive neuropsychological testing for cognitive impairment) is a subset of the independent test cohort, that consisted of 56 subjects (47 males, 9 females), demographically matched with the discovery cohort, with one or more testing visits in our lab (for a total of 111 visits). Positive neuropsychological testing was defined as a diagnosis of MCI, ADRD (Alzheimer Disorder Related Dementia), or other dementia upon neuropsychological testing done in a clinical setting, triggered by clinical concerns as part of regular clinical care (Fig. 1, Table 1, and Table S1). Subjects were matched for age at time of lab visit between the positive neuropsychological testing group and the rest group.
Medications
The subjects in the discovery cohort were all diagnosed with various psychiatric disorders (Table 1), and had various medical co-morbidities. Their medications were listed in their electronic medical records, and documented by us at the time of each testing visit. Medications can have a strong influence on gene expression. However, our discovery of differentially expressed genes was based on within-subject analyses, which factor out not only genetic background effects but also minimizes medication effects, as the subjects rarely had major medication changes between visits. Moreover, there was no consistent pattern of any particular type of medication, as our subjects were on a wide variety of different medications, psychiatric and nonpsychiatric. Furthermore, the independent validation/testing cohorts’ gene expression data was Z-scored by gender and diagnosis before being combined, to normalize for any such effects. Some subjects may be noncompliant with their treatment and may thus have changes in medications or drug of abuse not reflected in their medical records. That being said, our goal is to find biomarkers that track memory retention, regardless if the reason for it is endogenous biology or driven by substance abuse or medication noncompliance. In fact, one would expect some of these biomarkers to be direct or indirect targets of medications, as we show in this paper. Overall, the discovery, prioritization, and validation/ replication by testing in independent cohorts of the biomarkers, with our design, occurs despite the subjects having different genders, diagnoses, being on various different medications, and other lifestyle variables.
Blood gene expression experiments
RNA extraction
Whole blood (2.5–5 ml) was collected into each PaxGene tube by routine venipuncture. PaxGene tubes contain proprietary reagents for the stabilization of RNA. RNA was extracted and processed as previously described [7,8,9].
Microarrays
Microarray work was carried out using previously described methodology [7,8,9,10].
The dataset is available at GEO (GSE127711).
Biomarkers
Step 1: Discovery
We have used the subject’s score from the HVLT- DR Retention measure, assessed at the time of blood collection (Fig. 1). Using a 20% change threshold in Retention, we analyzed gene expression differences between visits, using a powerful within-subject design, then an across-subjects summation (Fig. 1).
We analyzed the data in two ways: an absent–present (AP) approach, and a differential expression (DE) approach, as in previous work by us on suicide biomarkers [7,8,9]. The AP approach may capture turning on and off of genes, and the DE approach may capture gradual changes in expression. We used a within-subject design, then an across-subjects summation score for probesets. Analyses were performed as previously described [8,9,10]. In brief, we imported all Affymetrix microarray data as CEL. files into Partek Genomic Suites 6.6 software package (Partek Incorporated, St Louis, MI, USA). Using only the perfect match values, we ran a robust multi-array analysis (RMA) by gender and diagnosis, background corrected with quantile normalization and a median polish probeset summarization of all chips, to obtain the normalized expression levels of all probesets for each chip. Then, to establish a list of differentially expressed probesets we conducted a within- subject analysis, using a fold change in expression of at least 1.2 between high- and low-memory visits within each subject. For each comparison, probesets that have a 1.2-fold change are then assigned either a 1 (increased in high memory) or a -1 (decreased in high memory; 0.5 or −0.5 if the change is between 1.1 and 1.2 fold; and 0 if the change is less than 1.1 fold. These values were then summed for each probeset across all the comparisons and subjects, yielding a range of raw scores. The raw scores were converted into internal scores for the next step, CFG prioritization: probesets above the 33.3% of raw scores received an internal score of 2 points, those above 50% 4 points, and those above 80% 6 points [8,9,10]. This was done in order to normalize and be able to easily combine internal (discovery) and external evidence (prioiritzation) scores, and to avoid over-fitting to the discovery cohort. We have developed in our lab R scripts to automate and conduct all these large dataset analyses in bulk, checked against human manual scoring [10].
Gene Symbol for the probesets were identified using NetAffyx (Affymetrix) for the Affymetrix HG-U133 Plus 2.0 Arrays which were used in this study, followed by GeneCards to confirm the primary gene symbol. In addition, for those probesets that were not assigned a gene symbol by NetAffyx, we used GeneAnnot or UCSC to obtain gene symbol for these uncharacterized probesets, followed by GeneCard. Genes were then scored using our manually curated CFG databases as described below (Fig. 1).
Step 2: Prioritization using Convergent Functional Genomics (CFG)
Databases
We have established in our laboratory (Laboratory of Neurophenomics, www.neurophenomics.info) manually curated databases of the human gene expression/protein expression studies (postmortem brain, peripheral tissue/fluids: CSF, blood and cell cultures), human genetic studies (association, copy number variations and linkage), and animal model gene expression and genetic studies, published to date on psychiatric disorders. Only findings deemed significant in the primary publication, by the study authors, using their particular experimental design and thresholds, are included in our databases. Our databases include only primary literature data and do not include review papers or other secondary data integration analyses to avoid redundancy and circularity. These large and constantly updated databases have been used in our CFG cross validation and prioritization platform (Fig. 1). For this study, data from 213 papers on AD were present in the databases at the time of the CFG analyses (August 2018) (human genetic studies-62, human brain tissue studies-49, human peripheral tissue/fluids- 83, nonhuman genetic studies-4, nonhuman brain tissue studies-13, nonhuman peripheral tissue/fluids- 2). Analyses were performed as previously described, using these databases [8, 9].
CFG scoring
For CFG analyses [15, 16], external cross-validating lines of evidence were weighted such that findings in human studies are prioritized over animal model studies, and brain studies are prioritized over peripheral fluid or genetic findings, by giving them twice as many points. Each line of evidence is capped in such a way that any positive findings within that line of evidence result in maximum points, regardless of how many different studies support that single line of evidence, to avoid potential popularity biases. Thus, genetic data had a maximum of two points, human brain expression data 4 points, human peripheral expression data 2 points, nonhuman genetic data 1 point, nonhuman brain expression data 2 points, and nonhuman peripheral expression data had 1 point, for a total of 12 points (Table S2). In addition to our external score, we also prioritized genes based upon the initial differential expression analyses used to identify them. Probesets identified by within-subject differential expression analyses can receive a maximum of 6 internal points. Thus, the maximum possible total CFG score for each gene is 18 points (6 points for the internal score and 12 points for the external score), with the external evidence weighted twice as much as the internal evidence to prevent a fit to cohort effect. It has not escaped our attention that other ways of scoring the lines of evidence may give slightly different results in terms of prioritization, if not in terms of the list of genes per se. Nevertheless, this simple scoring system provides a good separation of genes based on differential expression and on independent cross-validating evidence in the field.
Choice of biomarkers to be carried forward
We carried forward into testing the candidate biomarkers after the prioritization step, using as threshold a CFG score ≥ 10 (n = 138 probesets, 112 genes). A CFG score of 10 or above reflects an empirical cutoff of over half of the maximum possible total CFG score of 18, and permits the inclusion of potentially novel genes (with a maximal internal score of 6) and with some solid external prior literature evidence for involvement in AD. Of these, the top candidate biomarkers had a CFG score ≥ 12 (n = 23 probesets, 18 genes). In Step 3, testing, we then predict in independent cohorts state (low memory retention), and trait (future positive neuropsychological testing for cognitive impairment).
Diagnostics
In Step 3, testing, the test cohort for predicting low memory retention (state), and the test cohort for predicting Future Positive Neuropsychological Testing (trait), were completely independent from the discovery and validation cohorts; there was no subject overlap with them. They were assembled out of data that was RMA normalized by gender and diagnosis. The expression values of markers used for predictions were furthermore Z scored by gender and diagnosis in order to avoid potential artefacts due to different ranges of expression in different gender and diagnoses, and to be able to combine different markers into panels. Markers in panels were combined by simple summation of the increased risk markers minus the decreased risk markers. Predictions were performed using R-studio. For cross-sectional analyses, we used marker expression levels, z-scored by gender and diagnosis. For longitudinal analyses, we combined four measures: marker expression levels, slope (defined as ratio of levels at current testing visit vs. previous visit, divided by time between visits), maximum levels (at any of the current or past visits), and maximum slope (between any adjacent current or past visits). For decreased markers, we used the minimum rather than the maximum for level calculations. All four measures were each Z-scored by gender and diagnosis, then combined in an additive fashion into a single measure. The longitudinal analysis was carried out in a sub-cohort of the testing cohort consisting of subjects that had at least two test visits.
Predicting state—low memory
Receiver-operating characteristic (ROC) analyses between marker levels and memory state were performed by assigning subjects visits with a HVLT Retention score of ≤40 into the Low Memory category. We used the pROC package of R [17] (Fig. 2).
Predicting trait-future positive neuropsychological testing for cognitive impairment
We conducted analyses for predicting future positive neuropsychological testing performed as part of routine clinical care in subjects that had follow-up in the VA system, in which we have access to complete electronic medical records, using electronic medical records follow-up data of our study subjects (up to 12.81 years from initial visit at the time of the analyses). Analyses between marker measures (cross-sectional, longitudinal) at a specific testing visit and future positive neuropsychological test were performed a as described below, based on assigning if subjects had a future positive neuropsychological test for cognitive impairment or not. A Cox regression was performed using the time in days from the lab testing visit date to the positive neuropsychological testing date. The hazard ratio was calculated such that a value greater than 1 always indicates increased risk for positive neuropsychological testing, regardless if the biomarker is increased or decreased in expression.
Biological understanding
Clock gene database
We compiled a database of genes associated with circadian function, by using a combination of review papers and searches of existing databases CircaDB, GeneCards, and GenAtlas. Using the data we compiled from these sources we identified a total of 1468 genes that show circadian functioning. We further classified genes into “core” clock genes, i.e. those genes that are the main engine driving circadian function (n = 18), “immediate” clock genes, i.e., the genes that directly input or output to the core clock (n = 331), and “distant” clock genes, i.e., genes that directly input or output to the immediate clock genes (n = 1119).
Pathway analyses
IPA (Ingenuity Pathway Analysis, Qiagen), David Functional Annotation Bioinformatics Microarray Analysis (National Institute of Allergy and Infectious Diseases), and Kyoto Encyclopedia of Genes and Genomes (KEGG) (through DAVID) were used to analyze the biological roles, including top canonical pathways and diseases (Table 2), of the candidate genes resulting from our work. We conducted analyses for the 112 unique genes (from 138 probesets) that came out of the prioritization step, and for the top 18 unique genes (from 23 probesets). For Network analysis of the 112 unique genes we performed STRING Interaction Network by inputting the genes into the search window and performed Multiple Proteins Homo sapiens analysis.
CFG beyond Alzheimer’s: evidence for involvement in other psychiatric and related disorders
We also used a CFG approach to examine evidence from other psychiatric and related disorders, for the list of top biomarkers after Step 3 testing (n = 36 genes, from 42 probesets) (Table S3).
Therapeutic
Pharmacogenomics
We analyzed which of our top biomarkers from Table 3 (n = 36 genes, from 42 probesets) are known to be modulated by existing drugs using our CFG databases, and using Ingenuity Drugs analyses (Table 3 and Table S4).
New drug discovery/repurposing
We also analyzed which drugs and natural compounds are a match for the gene expression profiles of panels of our top biomarkers, using the Connectivity Map (https://portals.broadinstitute.org, Broad Institute, MIT), and the NIH LINCS L1000 database (Table 4).
Convergent functional evidence (CFE)
For the top biomarkers (Table 3 and Fig. 3), we tabulated into a convergent functional evidence (CFE) score all the evidence from discovery (up to 6 points), prioritization (up to 12 points), testing (state, trait - up to 6 points each if significantly predicts in all subjects, 4 points if predicts by gender, 2 points if predicts in gender/diagnosis). The total score can be up to 30 points: 18 from our data and 12 from literature data. We weigh our data more than the literature data, as ours is direct functional evidence for the involvement of the gene/biomarker in the phenotype. The goal is to highlight, based on the totality of our data and of the evidence in the field to date, biomarkers that have all around evidence: track memory, are reflective of AD and related pathology, and predict it. Such biomarkers merit priority evaluation in future clinical trials.
Results
In Step 1 Discovery, we identified blood gene expression biomarkers that track memory using the Retention measure from HVLT. At a phenotypic level, the Retention measure normalizes memory retention measurements in each subject during a testing visit, comparing the latter delayed recall trials with the earlier ones in that particular testing visit. We then used a powerful within –subject design in a cohort of subjects who displayed at least a 20% change in the Retention measure between different testing visits (n = 159 subjects, with 496 visits), to identify differentially expressed genes that track memory retention. Using our 33% of maximum raw score threshold (internal score of 2 pts, we had 10,941 unique probesets (Fig. 1). These were carried forward to the prioritization step. This represents approximately a 5-fold enrichment of the 54,625 probesets on the Affymetrix array.
In Step 2 Prioritization, we used a Convergent Functional Genomics (CFG) approach to prioritize the the candidate biomarkers identified in the discovery step (33% cutoff, internal score of ≥ 2 pt.) by using published literature evidence (genetic, gene expression and proteomic), from human and animal model studies, for involvement in AD (Fig. 1 and Table S2). There were 138 probesets that had a CFG score (combined internal and external score) of 10 and above. These were carried forward to the testing in independent cohorts step (see Supplementary Information- Complete Data and Analyses). This represents approximately a 400-fold enrichment of the probesets on the Affymetrix array.
In Step 3 Testing, we examined in an independent cohort (n = 127) from the one used for discovery whether the longer list of biomarkers prioritized in Step 2 (n = 138, CFG ≥ 10) are predictive of memory retention measure (state), and of future positive neuropsychological testing for MCI, AD or other dementia (trait), using electronic medical records follow-up data of our study subjects (up to 12.81 years from initial visit at the time of the analyses- December 2018) (Fig. 1, Table 1, and Table S1). The gene expression data in the test cohorts was normalized (Z-scored) across genders and various psychiatric diagnoses, before those different demographic groups were combined. We used biomarker levels information cross-sectionally, as well as expanded longitudinal information about biomarker levels at multiple visits, as predictors. We tested the biomarkers in all subjects in the independent test cohort, as well as in a more personalized fashion by gender and psychiatric diagnosis. We were successful in identifying gene expression biomarkers that were predictive in the independent cohorts for memory state, and for future neuropsychological testing positive for cognitive decline. We show increased accuracy with the personalized approach. In general, the longitudinal information was more predictive than the cross-sectional information.
Top predictive biomarkers for state were ITPKB, NKTR, RGS10, and PER1 (Fig. 2a and Table 3). The AUC ROCs ranged from over 0.6 for all subjects tested, to over 0.8 personalized by gender, and over 0.9 personalized by gender and diagnosis. NKTR is illustrative in this regard. It has a very modest AUC across all (0.59/p = 1.40E−02), higher in males (0.62/p = 5.55E−03), even more so in males with schizophrenia using cross-sectional levels of the biomarker (0.72/p = 1.38E-02), and the highest in male with schizophrenia using longitudinal levels of the biomarker (1/p = 1.35E−03). NKTR (Natural Killer Cell Triggering Receptor) is increased in expression in blood in low memory states in our study. There is some previous evidence showing it is increased in expression in the hippocampus in Alzheimer [18], but not as a top finding, and there is no previous functional evidence as a blood biomarker for memory or AD. Interestingly, it is a top biomarker for stress, increased in expression in high stress states stress states in a recent study from our group [12]. NKTR was also reported increased in expression in blood in studies of social isolation in humans [19], and in brain in studies of chronic variable stress in mice by Nestler and colleagues [20]. NKTR is also increased in expression in our previous blood biomarker studies of suicide, in both males [8, 21], and females [22], as well as increased in expression in postmortem brain studies in depression [23] and in schizophrenia [24], possibly underlying the effect of stress on those disorders and, based on our new data, on decreasing memory retention and promoting AD.
Top predictive biomarkers for trait were FCGR1A, GAP43, and MAPT (Fig. 2b and Table 3). Their Cox Regression Odds Ratios are at least 2-fold or higher, and significant. Caution should be used in interpreting these results because of the small Ns. We illustrate the virtues of personalization by gender and diagnosis with the example of RHEB in male schizophrenia, as our only future ADRD conversion to date was a male with schizophrenia (Figure S2A, B).
Based on our studies and analyses, the biomarkers with the best overall convergent functional evidence (CFE) for relevance to memory and AD were some new genes such as RAB7A, NPC2, TGFB1, GAP43, ARSB, PER1, GUSB, and MAPT (tau), as well as the well-known GSK3B, PTGS2, APOE, BACE1, PSEN1, and TREM2 (Table 3 and Fig. 3), with previous genetic and/or animal model evidence, that now have functional evidence in humans from our studies. The fact that key genes for AD brain pathology came out of our unbiased whole-genome discovery step tracking memory is reassuring, and they serve as de facto positive controls for our approach. Some of our new genes have previous supportive convergent evidence, buried in previous datasets; they were not considered strong candidate genes to be involved in AD before, and had no functional evidence as blood biomarkers for memory. An example is ARSB, a top biomarker from discovery, and our strongest predictor for low memory state in subjects with psychosis (SZ, SZA), with an AUC of 0.88/p-value 1.04E-03. ARSB (arylsulfatase B) is decreased in expression in blood low memory states in our studies. It has previous evidence for being decreased in expression in hippocampus and other brain regions in MCI [25], and genetic evidence for association with hippocampal volume in AD [26]. Additional new genes/biomarkers, with strong evidence for tracking memory state in the discovery step, are NKTR (natural killer cell triggering receptor), RHEB (Ras homolog enriched in brain), UBE2L3 (ubiquitin conjugating enzyme E2L 3), FOXO3 (forkhead box O3), SERTAD3 (SERTA domain containing 3), and UBE2I (ubiquitin conjugating enzyme E2I).
Some of the biomarkers are targets of existing drugs, such as lithium, antidepressants, and omega-3 fatty acids (Fig. 4, Table S4), of potential utility in patient stratification and pharmacogenomics approaches. Moreover, the top biomarkers gene expression signature, upon bioinformatics drug repurposing analyses (Table 4), yielded new drug candidates (such as pioglitazone and levonorgestrel), and natural compounds (such as salsolidine, ginkgolide A and icariin). The signature could be used for targeted enrollment of patients in clinical trials for these compounds, which would increase the odds of success, and for objectively measuring response to treatment.
Discussion
We describe a novel approach for discovering blood biomarkers of relevance to AD, and testing them in independent cohorts. The biomarkers also open a window into understanding the biology of memory disorders in general, and of AD and related disorders in particular.
Brain–blood
Blood biomarkers offer real-world clinical practice advantages. As the brain cannot be readily biopsied in live individuals, and CSF is less easily accessible than blood, we have endeavored over the years to identify blood biomarkers for neuropsychiatric disorders. A whole–blood approach facilities field deployment of sample collection. The assessment of gene expression changes focuses our approach on immune cells. The ability to identify peripheral gene expression changes that reflect brain activities is likely due to the fact that the brain and immune system have developmental commonalities, marked by shared reactivity and ensuing gene expression patterns. There is also a bi-directional interaction between the brain and immune system. Not all changes in expression in peripheral cells are reflective of or germane to brain activity. By carefully tracking a phenotype with our within subject design, and using convergent functional genomics prioritization, we are able to extract the peripheral changes that do track and are relevant to the brain activity studied, in this case memory and AD.
Power considerations
Based on our work for the last two decades in genetics and gene expression, along with the results of others in the field, we estimate that the within-subject longitudinal design, by factoring out all genetic and some environmental variability, is up to 3 orders of magnitude more powerful than an inter-subject case-control cross-sectional design. Moreover, gene expression, by integrating the effects of many SNPs and environment, is up to 3 orders of magnitude more powerful than a genetic study. Combined, our approach may be up to 6 orders of magnitude more powerful than a GWAS study, even prior to the CFG literature-based prioritization step.
As an illustration of this, after the completion of our analyses, two major genetic studies came out. The first was a large scale genetic meta-analysis of clinically diagnosed AD and AD-by-proxy (71,880 cases, 383,378 controls) (Jansen et al. Nature Genetics 2019) [27]. 14 of their 29 genome-wide significant loci/ genes were present in the candidate gene expression biomarkers for memory list that had survived our initial whole-genome unbiased within-subject Discovery step, before any CFG literature prioritization (APOE, TREM2, ABCA7, APH1B, BIN1, MS4A6A, PICALM, SLC24A4, HLA-DRB1, CNTNAP2, CASS4, CD33, CD2AP, and ADAM10). Two more of the 29 had anti-sense transcripts in our candidate biomarker list from Discovery (EPHA1-AS1 and ABI3-AS1) (see Supplementary Information- Complete Data and Analyses). The second was a large scale genetic meta-analysis of clinically diagnosed late-onset Alzheimer’s Disease (94,437 individuals) [28]. 15 of their 24 genome-wide significant loci/ genes were present in the candidate biomarker for memory list that had survived our initial whole-genome unbiased within-subject Discovery step, before any CFG literature prioritization (TREM2, ABCA7,INPP5D, BIN1, MS4A6A, PICALM, FERMT2, SLC24A4, HLA-DRB1, CASS4, CD2AP, IQCK, WWOX, ACE and ADAM10). One more of the 24 had anti-sense transcripts in our candidate biomarker list from Discovery (EPHA1-AS1) (see Supplementary Information- Complete Data and Analyses). This independent reproducibility of findings between our studies and the genetic studies, which are done in independent cohorts from ours, with independent methodologies, is reassuring, and provides strong convergent evidence for the validity and relevance of our approach and of the GWAS approach. Our work also provides functional evidence for some of the GWAS top hits. For example, APOE, a strong finding from genetic studies of AD, seems to be involved in pathology by levels of expression also, not just mutant forms. In fact, there is prior evidence that APOE promoter region SNPs may be involved in the pathology of AD [29]. APOE may be involved in particular in memory impairment [30].
Pathophysiological insights
A number of top biomarkers identified by us have biological roles that are related to the circadian clock (Table S3). To be able to ascertain all the genes in our dataset that were circadian and do estimates for enrichment, we compiled from the literature a database of all the known circadian genes, numbering a total of 1468 genes. Using an estimate of about 21 000 genes in the human genome, that gives about 7% of genes having some circadian pattern. Out of our 112 candidate biomarker genes, 19 had circadian evidence (16.96%), including the clock core gene PER1, suggesting a 2.4-fold enrichment for circadian genes. This provides a molecular underpinning for the epidemiological data between disrupted sleep and risk of AD [31, 32], and for the clinical phenomenology of “sundowning”.
The top candidate biomarkers also had prior evidence of involvement in other psychiatric and related disorders (Table S3), providing a molecular underpinning for the possible precursor effects of these disorders in AD. In particular, a majority of them have an overlap with stress (84%), aging (49%) and suicide (75%), consistent with them being part of the effects of stress on aging and the “life switch”, as described in a previous study by us [14]. The primary medications used to treat stress disorders are serotonin reuptake inhibitors (SSRIs). A recent study has shown that long-term SSRI treatment may delay progress from MCI to AD [33].
The top biological pathways where the candidate biomarkers from discovery and prioritization map were related to LXR/RXR activation, neuroinflammation signaling and atherosclerosis signaling (Table 2B). Agonists of nuclear receptors LXR:RXR and PPAR:RXR are known to ameliorate AD-related cognitive impairment and amyloid accumulation in murine models of AD. Interestingly, the ω-3 fatty acid docosahexaenoic acid (DHA), in combination with RXR agonist bexarotene, enhances LXR:RXR target gene expression of Abca1 and ApoE, reduces soluble forms of Aβ, and abrogates release of pro-inflammatory cytokines and mediators both in vitro and in a mouse model of AD [34]. Neuroinflammation is postulated to be a central mechanism in AD [35]. Atherosclerosis, and more general cardiovascular disorders, are recognized clinical risk factors for AD. The STRING gene interaction analysis (Fig. S3) revealed at least 3 networks. Network 1 (red) includes TREM2, along with ITPKB and FCGR1A; it may be involved in reactivity and inflammatory responses. Network 2 (green) includes GSK3B, along with MAPT (tau) and BACE1; it may be involved in activity and cellular trophicity. Network 3 (blue) includes PSEN1, along with FOXO3 and RHEB; it may be involved in connectivity and synaptic integrity. APOE is at the overlap of Networks 1 and 2. It is conceivable mechanistically that stress can lead to inflammation, followed by neuronal death (apoptosis) and synaptic connectivity loss, resulting in decreased memory and ultimately dementia.
The majority of top blood biomarkers we have identified have prior evidence in postmortem brain datasets from AD, which indicates their relevance to the pathophysiology of AD (Table S2). The co-directionality of blood changes in our work and brain changes reported in the literature needs to be interpreted with caution, as it may depend on brain region and on disease stage. Nevertheless, it is intriguing that three well known genes for AD, TREM2, MAPT (tau), and BACE1 are changed in the same direction (increased in expression) in blood in high memory states in our cohorts as they are in AD samples (Table S2). Consistent with the above discussion, RHEB is also changed in same direction (decreased) in high memory states in our data, and in AD brains in previous datasets, which increases BACE1 levels and amyloid β generation [36]. APP itself does not make the threshold as a candidate biomarker for memory in our discovery approach, but BACE 1 does, and it is increased in expression in high memory states. So is the APP-related gene APLP2 (data not shown). This opens the possibility that these genes are part of normal memory function, and their increase in AD is compensatory and/or excessive, like a scar formation. If so, treatments that target them may backfire in early stages of the disease, or need to be very carefully dosed in later stages of the disease. It is notable that various efforts to target amyloid accumulation have not been successful in multiple clinical trials to date.
The fact that GSK3B is a top candidate biomarker, has interesting neurobiological and therapeutic implications. GSK3B is decreased in expression in high memory states in our current work, and is mostly increased in expression in psychiatric disorders, including depression and stress related ones (Table S3), suggesting an avenue to their impact on memory and on later life AD. It is a known target of lithium [37], a medication that could be used for memory disorders improvement and prevention of AD, particularly in individuals with mood disorders earlier in life. Lithium is known to have anti-apoptotic roles [38, 39], so it is possible that it may prevent the neuronal loss of viability that is a precursor and part of AD, and which in turn triggers and is accentuated by the scar-like deposition of MAPT and APP.
Pharmacogenomics and drug repurposing
Besides GSK3B, lithium also modulates expression of MAPT, PER1, PTGS2 and IGF1 (Fig. 4). GSK3B and MAPT are also modulated by omega-3 fatty acids. Such information could be useful for personalized treatment selection, and for monitoring response to treatment. Moreover, the drug repurposing analyses using the biomarker signatures from our studies (Table 4) yielded drug candidates such as pioglitazone (a diabetes medication, previously studied for AD prevention in clinical trials that may have been too heterogeneous in enrollment), levonorgestrel (a progesterone derivative), and the natural compounds salsolidine (a natural compound with acetylcholinesterase inhibitory properties), ginkgolide A (a natural compound with neuroprotective effects) and icariin (a natural compound reported to improve memory impairment in Alzheimer’s disease model mice). We now have human functional data-based evidence for the potentially utility of these drugs. Various abbreviated signatures present in individual patients could be used in the future to choose the best repurposed drug fit for each patient.
Diagnostics
Our data for assessment and predictions needs to be interpreted with caution, particularly for the trait predictions, due to the small Ns. Nevertheless, the results are obtained in an independent cohort from the one used in discovery, and suggest interesting leads, as well as the fact that gender differences will be important for clinical applications. We would propose that a panel of our top genes, such as those in Table 3, be tested for predictions in other independent cohorts, and potentially be used in clinical settings. A comination of phenomic data (clinical information, imaging) and and biomarker data may work best. Of note, however, the HVLT was not a significant trait predictor in our cohort (data not shown).
Conclusions
Our work has provided evidence for novel possible precision medicine approaches, diagnostic and therapeutic. In particular, it may lead to improved early objective assessment of state, of future risk, and to targeted preventive treatments (in essence, a risk evaluation and mitigation system) for memory disorders in general, and AD in particular, that result in decreased quality and quantity of life, at a massive cost to individuals, families and society. It also opens a novel window into disease pathophysiology, and may lead to new targets for drug development. Given the growing world-wide burden of AD, and the unsuccessful approaches to date, such new avenues should be pursued with vigor and alacrity.
References
Bilbo SD, Schwarz JM. The immune system and developmental programming of brain and behavior. Front Neuroendocrinol. 2012;33:267–86.
Le-Niculescu H, McFarland MJ, Ogden CA, Balaraman Y, Patel S, Tan J, et al. Phenomic, convergent functional genomic, and biomarker studies in a stress-reactive genetic animal model of bipolar disorder and co-morbid alcoholism. Am J Med Genet B Neuropsychiatr Genet. 2008;147B:134–66.
Le-Niculescu H, Kurian SM, Yehyawi N, Dike C, Patel SD, Edenberg HJ, et al. Identifying blood biomarkers for mood disorders using convergent functional genomics. Mol Psychiatry. 2009;14:156–74.
Kurian SM, Le-Niculescu H, Patel SD, Bertram D, Davis J, Dike C, et al. Identification of blood biomarkers for psychosis using convergent functional genomics. Mol Psychiatry. 2011;16:37–58.
Le-Niculescu H, Balaraman Y, Patel SD, Ayalew M, Gupta J, Kuczenski R, et al. Convergent functional genomics of anxiety disorders: translational identification of genes, biomarkers, pathways and mechanisms. Transl Psychiatry. 2011;1:e9.
Witt SH, Sommer WH, Hansson AC, Sticht C, Rietschel M, Witt CC. Comparison of gene expression profiles in the blood, hippocampus and prefrontal cortex of rats. In Silico Pharmacol. 2013;1:15.
Le-Niculescu H, Levey DF, Ayalew M, Palmer L, Gavrin LM, Jain N, et al. Discovery and validation of blood biomarkers for suicidality. Mol Psychiatry. 2013;18:1249–64.
Niculescu AB, Levey DF, Phalen PL, Le-Niculescu H, Dainton HD, Jain N, et al. Understanding and predicting suicidality using a combined genomic and clinical risk assessment approach. Mol Psychiatry. 2015;20:1266–85.
Levey DF, Niculescu EM, Le-Niculescu H, Dainton HL, Phalen PL, Ladd TB, et al. Towards understanding and predicting suicidality in women: biomarkers and clinical risk assessment. Mol Psychiatry. 2016;21:768–85.
Niculescu AB, Le-Niculescu H, Levey DF, Phalen PL, Dainton HL, Roseberry K, et al. Precision medicine for suicidality: from universality to subtypes and personalization. Mol Psychiatry. 2017;22:1250–73.
Niculescu AB, Le-Niculescu H, Levey DF, Roseberry K, Soe KC, Rogers J, et al. Towards precision medicine for pain: diagnostic biomarkers and repurposed drugs. Mol Psychiatry. 2019;24:501–22.
Le-Niculescu H, Roseberry K, Levey DF, Rogers J, Kosary K, Prabha S et al. Towards precision medicine for stress disorders: diagnostic biomarkers and targeted drugs. Mol Psychiatry. 2019. [Epub ahead of print].
Nho K, Ramanan VK, Horgusluoglu E, Kim S, Inlow MH, Risacher SL, et al. Comprehensive gene- and pathway-based analysis of depressive symptoms in older adults. J Alzheimer’s Dis. 2015;45:1197–206.
Rangaraju S, Levey DF, Nho K, Jain N, Andrews KD, Le-Niculescu H, et al. Mood, stress and longevity: convergence on ANK3. Mol Psychiatry. 2016;21:1037–49.
Bertsch B, Ogden CA, Sidhu K, Le-Niculescu H, Kuczenski R, Niculescu AB. Convergent functional genomics: a Bayesian candidate gene identification approach for complex disorders. Methods. 2005;37:274–9.
Le-Niculescu H, Patel SD, Niculescu AB. Convergent integration of animal model and human studies of bipolar disorder (manic-depressive illness). Curr Opin Pharmacol. 2010;10:594–600.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinf. 2011;12:77.
Blalock EM, Geddes JW, Chen KC, Porter NM, Markesbery WR, Landfield PW. Incipient Alzheimer’s disease: microarray correlation analyses reveal major transcriptional and tumor suppressor responses. Proc Natl Acad Sci USA. 2004;101:2173–8.
Cole SW, Hawkley LC, Arevalo JM, Sung CY, Rose RM, Cacioppo JT. Social regulation of gene expression in human leukocytes. Genome Biol. 2007;8:R189.
Labonte B, Engmann O, Purushothaman I, Menard C, Wang J, Tan C, et al. Sex-specific transcriptional signatures in human depression. Nat Med. 2017;23:1102–11.
Karege F, Perroud N, Burkhardt S, Fernandez R, Ballmann E, La Harpe R, et al. Protein levels of beta-catenin and activation state of glycogen synthase kinase-3beta in major depression. A study with postmortem prefrontal cortex. J Affect Disord. 2012;136:185–8.
Canli T, Wen R, Wang X, Mikhailik A, Yu L, Fleischman D, et al. Differential transcriptome expression in human nucleus accumbens as a function of loneliness. Mol Psychiatry. 2017;22:1069–78.
Blasi G, Napolitano F, Ursini G, Di Giorgio A, Caforio G, Taurisano P, et al. Association of GSK-3beta genetic variation with GSK-3beta expression, prefrontal cortical thickness, prefrontal physiology, and schizophrenia. Am J Psychiatry. 2013;170:868–76.
Glatt SJ, Everall IP, Kremen WS, Corbeil J, Sasik R, Khanlou N, et al. Comparative gene expression analysis of blood and brain provides concurrent validation of SELENBP1 up-regulation in schizophrenia. Proc Natl Acad Sci USA. 2005;102:15533–8.
Berchtold NC, Sabbagh MN, Beach TG, Kim RC, Cribbs DH, Cotman CW. Brain gene expression patterns differentiate mild cognitive impairment from normal aged and Alzheimer’s disease. Neurobiol Aging. 2014;35:1961–72.
Potkin SG, Guffanti G, Lakatos A, Turner JA, Kruggel F, Fallon JH, et al. Hippocampal atrophy as a quantitative trait in a genome-wide association study identifying novel susceptibility genes for Alzheimer’s disease. PLoS ONE. 2009;4:e6501.
Jansen IE, Savage JE, Watanabe K, Bryois J, Williams DM, Steinberg S et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nature Genet. 2019;51:404–13.
Kunkle BW, Grenier-Boley B, Sims R, Bis JC, Damotte V, Naj AC, et al. Genetic meta-analysis of diagnosed Alzheimer’s disease identifies new risk loci and implicates Abeta, tau, immunity and lipid processing. Nat Genet. 2019;51:414–30.
Maloney B, Ge YW, Petersen RC, Hardy J, Rogers JT, Perez-Tur J, et al. Functional characterization of three single-nucleotide polymorphisms present in the human APOE promoter sequence: Differential effects in neuronal cells and on DNA-protein interactions. Am J Med Genet Part B, Neuropsychiatr Genet. 2010;153B:185–201.
Mukherjee S, Mez J, Trittschuh EH, Saykin AJ, Gibbons LE, Fardo DW et al. Genetic data and cognitively defined late-onset Alzheimer’s disease subgroups. Mol Psychiatry. 2018. [Epub ahead of print].
Noble W, Spires-Jones TL. Sleep well to slow Alzheimer’s progression? Science. 2019;363:813–4.
Irwin MR, Vitiello MV. Implications of sleep disturbance and inflammation for Alzheimer’s disease dementia. Lancet Neurol. 2019;18:296–306.
Bartels C, Wagner M, Wolfsgruber S, Ehrenreich H, Schneider A, Alzheimer’s Disease, et al. Impact of SSRI therapy on risk of conversion from mild cognitive impairment to Alzheimer’s dementia in individuals with previous depression. Am J Psychiatry. 2018;175:232–41.
Casali BT, Corona AW, Mariani MM, Karlo JC, Ghosal K, Landreth GE. Omega-3 fatty acids augment the actions of nuclear receptor agonists in a mouse model of Alzheimer’s disease. J Neurosci. 2015;35:9173–81.
Kinney JW, Bemiller SM, Murtishaw AS, Leisgang AM, Salazar AM, Lamb BT. Inflammation as a central mechanism in Alzheimer’s disease. Alzheimers Dement. 2018;4:575–90.
Shahani N, Pryor W, Swarnkar S, Kholodilov N, Thinakaran G, Burke RE, et al. Rheb GTPase regulates beta-secretase levels and amyloid beta generation. J Biol Chem. 2014;289:5799–808.
Mendes CT, Mury FB, de Sa Moreira E, Alberto FL, Forlenza OV, Dias-Neto E, et al. Lithium reduces Gsk3b mRNA levels: implications for Alzheimer disease. Eur Arch psychiatry Clin Neurosci. 2009;259:16–22.
Fan M, Jin W, Zhao H, Xiao Y, Jia Y, Yin Y, et al. Lithium chloride administration prevents spatial learning and memory impairment in repeated cerebral ischemia-reperfusion mice by depressing apoptosis and increasing BDNF expression in hippocampus. Behavioural Brain Res. 2015;291:399–406.
Liechti FD, Studle N, Theurillat R, Grandgirard D, Thormann W, Leib SL. The mood-stabilizer lithium prevents hippocampal apoptosis and improves spatial memory in experimental meningitis. PLoS ONE. 2014;9:e113607.
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, et al. The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006;313:1929–35.
Lamb J. The Connectivity Map: a new tool for biomedical research. Nat Rev Cancer. 2007;7:54–60.
Acknowledgements
This work is, de facto, a field-wide collaboration. We would like to acknowledge the many other groups, cited in our paper, who have conducted and published studies in Alzheimer and related disorders. With their arduous and careful work, a convergent approach such as ours is possible. We would like to thank Andrew Burkhart for literature searches and database work. We would particularly like to thank the subjects in our studies, their families and their caregivers. Without their contribution, such work to advance the understanding of memory disorders and Alzheimer’s would not be possible. This work was supported by an NIH Directors’ New Innovator Award (1DP2OD007363) and a VA Merit Award (2I01CX000139) to ABN.
Author information
Authors and Affiliations
Contributions
ABN conceived and designed the study, and wrote the manuscript. HLN, KR, AK and HR analyzed the data. TJ, AS, AW, JH and SW organized, scored, and conducted testing in psychiatric subjects. SMK conducted microarray experiments. BL, AS, DKL and AJS provided input on methods and results interpretation. All authors reviewed the final manuscript and agreed with it.
Corresponding author
Ethics declarations
Conflict of interest
ABN is listed as inventor on a patent application being filed by Indiana University, and is a co-founder of MindX Sciences. The other authors declare no conflict of interest related to this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Niculescu, A.B., Le-Niculescu, H., Roseberry, K. et al. Blood biomarkers for memory: toward early detection of risk for Alzheimer disease, pharmacogenomics, and repurposed drugs. Mol Psychiatry 25, 1651–1672 (2020). https://doi.org/10.1038/s41380-019-0602-2
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41380-019-0602-2
This article is cited by
-
Precision medicine for psychotic disorders: objective assessment, risk prediction, and pharmacogenomics
Molecular Psychiatry (2024)
-
Leveraging electronic health records and knowledge networks for Alzheimer’s disease prediction and sex-specific biological insights
Nature Aging (2024)
-
Aging gene signature of memory CD8+ T cells is associated with neurocognitive functioning in Alzheimer’s disease
Immunity & Ageing (2023)
-
Lithium rescues dendritic abnormalities in Ank3 deficiency models through the synergic effects of GSK3β and cyclic AMP signaling pathways
Neuropsychopharmacology (2023)
-
Towards precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs
Molecular Psychiatry (2023)