Metabolic Profiling of Cognitive Aging in Midlife

Introduction

Human brain aging and Alzheimer’s dementia (AD) impose a huge health and economic burden on modern society. Recent research suggests that cognitive decline begins in middle age (45–65 years old) (Singh-Manoux et al., 2012; Karlamangla et al., 2017). However, there is a substantial heterogeneity in the trajectories of cognitive aging, with some individuals showing minimal or no decline while others of the same age experiencing rapid decline or even dementia. Discovery of sensitive, specific and non-invasive biomarkers is a prerequisite for identifying high-risk individuals, and developing effective strategies that can monitor, delay or prevent the onset of dementia.

A growing body of evidence suggests that brain and cognitive aging are accompanied by extensive metabolic perturbations (Clarke et al., 2018). The altered metabolic profiles can be quantified by metabolomics, a high-throughput biochemical technology that can identify hundreds to thousands of small molecules (metabolites) in biofluids or tissues (Rochfort, 2005). These metabolites represent the endpoints of metabolic processes encompassing the interaction between the internal genome and the external environment (Harrigan and Goodacre, 2012), and thus are closer to disease phenotypes compared to genomic, transcriptomic and proteomic profiles.

While the majority of the existing metabolomics studies in human (Han et al., 2011; Orešič et al., 2011; Björkqvist et al., 2012; Mapstone et al., 2014; Inoue et al., 2015; Casanova et al., 2016; Varma et al., 2018; Huo et al., 2020) have focused on older individuals (aged 65+), several recent studies reported metabolic changes associated with cognitive aging (Bressler et al., 2017; Chouraki et al., 2017; Darst et al., 2018; van der Lee et al., 2018; Proitsi et al., 2018; Tynkkynen et al., 2018) among individuals in midlife (45–65 years). However, most of these studies have examined a small number of metabolites (in general less than 300), leaving the full spectrum of blood metabolome largely unexplored. To date, no consensus metabolites that can predict brain and cognitive aging in midlife have been reported.

Leveraging the comprehensive clinical and cognitive phenotypes in a well-characterized longitudinal twin cohort, the current study aimed to identify metabolites predictive of incident mild cognitive impairment (MCI) or cognitive change using untargeted metabolomics in 30 middle-aged (mean age 57.2 at baseline), male-male twin pairs [15 monozygotic (MZ) pairs, 15 dizygotic (DZ) pairs]. All twins were cognitively normal at baseline (Wave 1, 2003–2007), but one twin in each of these pairs developed incident MCI, while his co-twin remained cognitively normal at the end of 5.5-year follow-up (Wave 2, 2008–2013). Our initial study was followed by replication in an independent cohort of middle age adults (mean age 59.1 at blood draw). Although preliminary, findings of this research are likely to generate novel hypotheses that may lead to the identification of novel metabolic markers predictive of cognitive aging in midlife and provide insights into our understanding of the role of metabolic disturbance in early cognitive aging.

Materials and Methods

Study Populations

Discovery Cohort – Vietnam Era Twin Study of Aging (VETSA)

Our discovery sample included 30 middle-aged, male-male twin pairs (15 MZ pairs and 15 DZ pairs) participating in the VETSA, a longitudinal observational study examining the role of genetic and environmental factors for cognitive and brain aging beginning in midlife using a twin design. Detailed methods for twin recruitment, longitudinal follow-up and phenotypes for the twins had been described previously (Kremen et al., 2006, 2013a). Briefly, 1,237 twins (349 MZ pairs, 265 DZ pairs, 9 singletons, age range 51–60, mean age 56) attended the initial examination at Wave 1 (2003–2007), and 1,016 twins were re-examined at Wave 2 (2008–2013, mean follow-up 5.5 years). The only two recruitment criteria for the VETSA were: (1) the ages of participants at Wave 1 were between 51 and 59 years at the time of initial recruitment (2003–2007), and (2) both twins in a pair were willing to participate in the baseline assessment. Twins enrolled in VETSA were predominately Caucasians (86%). Although all twins served in the United States military between 1965 and 1975, nearly 80% of them did not serve in combat or Vietnam (Kremen et al., 2006, 2011). Since 25% of the men nationwide within this age range served in the military, the VETSA participants were generally representative of middle-aged men living across the United States, with respect to health and sociodemographic characteristics (Schoenborn and Heyman, 2009). Participants were administered identical protocols at the University of California, San Diego, or Boston University. Individual twins chose their site, although brothers most often chose the same site. The complete protocols for twin enrollment and exams had been described previously (Kremen et al., 2013a). The protocols were approved by the University of California and Boston University Institutional Review Boards.

The current analysis included 30 complete pairs (mean age 57.2, age range 52.6–59.6, all Caucasians) attending the clinical exam at Wave 1 (baseline, 2003–2007) followed through Wave 2 (2008–2013). All twins were cognitively normal at baseline, but one twin per pair developed incident MCI by Wave 2 (mean follow-up period 5.5 years) while their co-twin did not. Both twins of the same pair were examined at the same site. The baseline characteristics of these twins are shown in Table 1. There was no significant difference in the listed characteristics at baseline between twins who converted to MCI (converters) versus their co-twins who did not (non-converters).

TABLE 1

Table 1. Baseline characteristics of twins in the VETSA (N = 60).

Replication Cohort – Wisconsin Registry for Alzheimer’s Prevention (WRAP)

We replicated our initial findings in the WRAP (La Rue et al., 2008), which is an ongoing prospective cohort that enrolls middle-aged adults and examines the transition from middle age to early older adults (between 40 and 65). Participants had no dementia at baseline. The WRAP allows for the enrollment of siblings and is enriched for a parental history of Alzheimer’s disease. Detailed information for the study design and methods has been described elsewhere (Sager et al., 2005; Johnson et al., 2018). The current analysis included 232 WRAP participants (mean age 59.1, age range 40–72, 32% males, all Caucasians) with complete clinical and metabolomics data at baseline (2011–2015) and follow-up (average follow-up period = 4.1 years). Supplementary Table 1 displays the baseline characteristics of the WRAP participants.

Cognitive Phenotypes and Diagnostics

Vietnam Era Twin Study of Aging (VETSA)

All twins enrolled in the VETSA underwent extensive neurocognitive testing including 13 neuropsychological tests (23 scores) covering 7 cognitive domains (Kremen et al., 2019) that were designed to avoid ceiling effects in middle-age adults. Specifically, general cognitive ability was measured by the Armed Forces Qualification Test (AFQT), a 50-min paper-and-pencil test with 100 multiple-choice items on vocabulary, arithmetic, spatial processing, and mechanical ability (Uhlaner and Bolanovich, 1952; Orme et al., 2001). The AFQT is highly correlated (∼0.85) with other measures of general cognitive ability, including the Wechsler scales of intelligence (Lyons et al., 2009, 2017). Executive function was assessed by a factor comprising 7 measures of prepotent response inhibition, task-set shifting, and working memory: Stroop interference; AX-Continuous Performance Test; Delis-Kaplan Executive System (D-KEFS) Trails switching; D-KEFS category switching; Wechsler Memory Scale-III (WMS-III) Letter-number sequencing; Reading span; and WMS-III Digit span (Gustavson et al., 2018). Episodic memory was measured by summarizing verbal and visual-spatial episodic memory measures, compromising 6 individual episodic memory measures (Kremen et al., 2014), including California Verbal Learning Test-II short/long delayed free recall; WMS-III logical memory immediate/delayed free recall, and visual reproductions immediate/delayed recall. Short-term memory was derived based on standardized and averaged scores from WMS-III digit span (forward condition) and spatial span (forward condition) (Franz et al., 2011). MCI status was determined via the Jak-Bondi approach, requiring at least 2 tests within a cognitive domain to each be more than 1.5 SDs below normative expectations in order to define impairment in that domain (Jak et al., 2009; Kremen et al., 2013b; Granholm et al., 2017).

WRAP

As previously described (Dowling et al., 2010; Koscik et al., 2014), the WRAP has also collected detailed neurocognitive phenotypes for all participants. Since the WRAP did not have the AFQT, we used the average of verbal and visuospatial ability from the Wechsler Abbreviated Scale of Intelligence (Koscik et al., 2014) as a surrogate measure for the general cognitive ability. Executive function in the WRAP was calculated by averaging the z-score of the following four components (Clark et al., 2016), including trail making test part B (TMT B) total time to completion; Stroop neuropsychological screening test color-word interference total items completed in 120 s; Wechsler abbreviated intelligence scale-revised (WAIS-R); and digit symbol coding total items completed in 90 s. We used a composite of delayed recall scores (Clark et al., 2016) as a surrogate measure for episodic memory: Rey Auditory Verbal Learning Test long-delay free recall, WMS-R logical memory delayed recall, Brief Visuospatial Memory Test-Revised delayed recall. Since WRAP did not collect short-term memory, we used immediate memory, assessed by a factor compromising Trials 1 and 2 of Rey Auditory Verbal Learning Test (Schmidt, 1996), in our replication analysis.

Metabolomics Data Acquisition

Fasting blood samples were drawn in the morning of the same day when cognitive tests were performed. Plasma samples were collected (EDTA tube), aliquoted and then stored at −80°C. Untargeted metabolomics data in fasting plasma samples were assayed in both VETSA (n = 60) and WRAP (n = 232) using the same protocols on the same platform at Metabolon Inc. Detailed procedures for metabolomics analysis have been described elsewhere (Evans et al., 2009, 2014). Briefly, automated MicroLab STAR^® system from Hamilton Company was used to remove protein and recover chemically diverse metabolites. The consequent extract was equally aliquoted into five parts for different analyses, including two separate reverse phase (RP)/UPLC-MS/MS methods for positive ion mode electrospray ionization (ESI), one RP/UPLC-MS/MS method for negative ion mode ESI, one HILIC/UPLC-MS/MS method for negative ion mode ESI, and the remaining one for backup purposes. Quality control samples were also included to ensure the Metabolon procedures were conducted within the specified protocol. Internal standards (endogenous compounds) were spiked into each analyzed sample to monitor instrument performance and aid chromatographic alignment. Positive and negative samples were randomized and QC samples were evenly spaced among the injections. After drying, the samples were reconstituted according to the pre-mentioned methods, followed by utilization of Ultrahigh Performance Liquid Chromatography-Tandem Mass Spectroscopy (UPLC-MS/MS). Metabolon’s software was used to extract the raw data and perform peak-identification, QC processing, and biochemical identification. A series of curation procedures were performed to guarantee the quality of the data. The relative intensity of each metabolite was quantified using area-under-the curve (AUC) of its corresponding peak.

Metabolomics Data Pre-processing and Quality Control

For the 30 twin pairs in VETSA, we obtained a total of 1,228 metabolites, including 960 known and 268 unknown compounds. The AUCs of all samples were first normalized to their extracted volume, followed by median normalization (average of the median levels in all samples) as previously described (van den Berg et al., 2006). The normalized data were further log2 transformed to improve normality. Metabolites with missing values in >50% of the samples were removed. Metabolites with missing value ≤ 50% of the samples were replaced by limit of detection (LOD). A total of 1,131 metabolites (879 knowns and 252 unknowns) was included in the final analysis in the VETSA. Similar procedures were used to pre-process and QC the metabolomics data in the WRAP (Darst et al., 2019), resulting in a total of 1,097 metabolites in our replication analysis.

Statistical Analyses

Identifying Metabolites Predictive of MCI Onset and Cognitive Change in VETSA

Discordant twin analyses

To identify plasma metabolites that can predict incident MCI, we first performed paired t-tests by comparing the relative abundance of metabolites between MCI converters and non-converters. These discordant twin analyses can be particularly powerful because twins within pairs are nearly perfectly matched for demographic, developmental, and environmental factors, resulting in considerably reduced error variance compared with non-twin studies. In addition, MZ pairs are perfectly matched for genetic background, and DZ pairs share, on average, 50% of their genetic background.

Non-twin analyses

While matched pair analysis is powerful, it does not allow for adjustments of other covariates. Thus, we also conducted non-twin analyses, meaning that the unit of analysis was the individual rather than the twin pair. This was done by constructing generalized linear mixed models via the R package lme4 (Bates et al., 2014), in which MCI status (y/n) was the dependent variable, and baseline metabolite level was the independent variable, adjusting for twin age, BMI, smoking, and education at baseline and zygosity (MZ vs. DZ). To identify metabolites that can predict the longitudinal change in cognitive function over time, we fitted linear mixed models, in which change in a cognitive measure (difference in a cognitive measure between Wave 1 and Wave 2) was the dependent variable and baseline metabolite level was the independent variable, adjusting for twin age, zygosity (MZ vs. DZ), education, BMI, smoking, and the cognitive measure at baseline. The effect size and P-value were estimated for each metabolite for each cognitive measure. In the above described models, twin pair was included as a random effect to adjust for the correlated observations. Multiple testing was controlled by false discovery rate (FDR) (Benjamini and Hochberg, 1995) (adjusting for a total of 1,228 metabolites × 5 cognitive phenotypes including MCI, general cognitive ability, executive function, episodic memory, and short-term memory), and FDR-adjusted P-value was used to determine statistical significance.

Replication in the WRAP

To replicate the metabolites identified in the VETSA, we fitted linear mixed models in the WRAP, in which cognitive change (difference in a cognitive measure between baseline and follow-up) was the dependent variable and baseline plasma level was the independent variable. The model adjusted for age, sex, education, BMI, and the cognitive measure under investigation at baseline. Family ID was included as a random effect in the model to account for relatedness between family members.

Network analysis

To identify metabolic networks (modules) associated with cognitive functions, we conducted the weighted gene correlation network analysis (WGCNA) (Langfelder and Horvath, 2008) using all 1,131 metabolites (both known and unknown) in the VETSA. The structure of each module was constructed using data from MCI non-converters (n = 30) via the ARACNE algorithm (Margolin et al., 2006) using R package minet (Meyer et al., 2008). To examine whether module structure predicts incident MCI, we performed modular differential connectivity (MDC) analysis (McKenzie et al., 2016) by comparing the network connectivity between MCI converters and non-converters. The statistical significance of the MDC analysis was assessed by 1,000 permutations. The network structure of the modules containing significant MDC was visualized by CytoScape (Shannon et al., 2003).

Functional annotation

To functionally annotate the metabolites predictive of MCI conversion or cognitive change over time, we conducted pathway enrichment analysis using the software MetaboAnalyst 4.0 (Chong et al., 2019). Metabolites with raw P < 0.05 were used as input, and hypergeometric test was employed to obtain P-values for metabolic pathways.

Results

Metabolites Predictive of Incident MCI or Cognitive Change in VETSA

At raw P < 0.001, we identified four plasma metabolites that can predict MCI conversion (Figure 1). Specifically, higher levels of three known metabolites [gulonate, β = 0.36, P = 5.5 × 10^–4; 1-palmitoyl-GPE (16:0), β = 0.35, P = 7.5 × 10^–4; and 1-palmitoyl-2-arachidonoyl-GPI (16:0/20:4) (16:0/20:4), β = 0.58, P = 9.7 × 10^–4] and a lower level of one unknown metabolite (X − 21959, β = −0.47, P = 7.2 × 10^–4) predicted MCI conversion. These P-values remained similar after adjustments for age, cognition, zygosity, education, BMI, and smoking status at baseline (Supplementary Table 2). However, none of them passed multiple comparison at FDR-adjusted P < 0.05.

FIGURE 1

Figure 1. Plasma metabolites predictive of incident MCI in twins (P < 0.001). P-values were obtained by paired t-test.

At FDR-adjusted P < 0.05, we found that baseline levels of four plasma metabolites significantly predicted longitudinal change in cognitive function over time (Table 2). Specifically, higher baseline levels of hydroxy-CMPF, linoleoyl-arachidonoyl-glycerol [DAG] (18:2/20:4), sphingomyelin [SM] (d18:1/20:1 and d18:2/20:0), and SM (d18:1/22:1, d18:2/22:0, and d16:1/24:1) significantly predicted a slower decrease in general cognitive ability (β = 0.16, q = 0.04), short-term memory (β = 0.16, q = 0.03), and executive function (β = 0.04, q = 0.03;β = 0.04, q = 0.04), respectively. Figure 2 depicts the associations of these metabolites with cognitive functions. No metabolites significantly predicted changes in episodic memory at FDR-adjusted P < 0.05.

TABLE 2

Table 2. Plasma metabolites whose baseline level predict cognitive changes over time in twins (q-value < 0.05).

FIGURE 2

Figure 2. Significant metabolites predictive of cognitive change over time in twins (FDR-adjusted P < 0.05). X-axis indicates the relative abundance in plasma metabolite (log2 transformed) at baseline. Y-axis indicates cognition change over time (Wave 1 to Wave 2). P-values were obtained by linear mixed models, in which change in a cognitive measure (difference in a cognitive measure between Wave 1 and Wave 2) was the dependent variable and metabolite level was the independent variable, adjusting for twin age, zygosity (MZ vs. DZ), education, BMI, smoking, and the cognitive measure at baseline. Twin pair was included as a random effect in the model.

Replication in the WRAP

Of the four metabolites predictive of cognitive change in the VETSA, we were able to replicate the association of SM (d18:1/20:1 and d18:2/20:0) (P = 0.003) with longitudinal change in executive function in the WRAP. The other 3 metabolites showed consistent effect size direction, but statistically non-significant.

Metabolic Networks

We identified 16 metabolite modules (Supplementary Figure 1 and Supplementary Table 3). Of these, the network connectivity (a measure quantifying the connections among metabolites) for the pink module (Figure 3) was significantly higher in MCI non-converters compared to that in converters (117 vs. 75, MDC = 0.64, p = 0.01), suggesting that reduced connectivity between metabolites may be associated with MCI onset. In addition, we identified several hub metabolites that differ significantly between MCI converters and non-converters. For example, 2,2′-methylenebis(6-tert-butyl-p-cresol), ceramide (d16:1/24:1, d18:1/22:1) and gamma-glutamylglutamate were the hub metabolites only in non-converters but not in converters.

FIGURE 3

Figure 3. Differential network analysis showing significant structural changes of the pink module by comparing twins who developed incident MCI versus their co-twin brothers who did not. Only hub metabolites were labeled. The hub metabolites for NCI network were 2,2′-Methylenebis(6-tert-butyl-p-cresol), ceramide (d16:1/24:1, d18:1/22:1), and gamma-glutamylglutamate. The hub metabolites for the MCI network were 2,2′-Methylenebis(6-tert-butyl-p-cresol), iminodiacetate (IDA), and sebacate (C10-DC).

Functional Annotation

The aminoacyl-tRNA biosynthesis pathway showed the most significant association with short-term memory (P = 4.59 × 10^–6), followed by the sphingolipid metabolism pathway associated with MCI (P = 1.56 × 10^–4), and the valine, leucine, and isoleucine biosynthesis pathway in relation to short-term memory (P = 5.02 × 10^–4). The top five metabolic pathways are shown in Table 3.

TABLE 3

Table 3. Top metabolic pathways derived using plasma metabolites predictive of MCI onset or cognitive change over time.

Discussion

In an untargeted metabolomic analysis including 30 middle-aged twin pairs discordant for progression to MCI, we identified four fast plasma metabolites whose baselines significantly predict cognitive change over time after adjusting for multiple testing and known covariates. Of these, the association between a sphingomyelin (SM [d18:1/20:1, d18:2/20:0]) and the longitudinal change in executive function was replicated in an independent cohort of middle-aged adults. Metabolic network analysis revealed that a reduced metabolic connection could be associated with progression to MCI. The identified metabolites were enriched in several metabolic pathways. Together, our results provide initial evidence that metabolic perturbation affects cognitive aging in midlife and suggest that plasma metabolites may serve as potential biomarkers for early cognitive aging and AD prediction in midlife.

Of the four metabolites predictive of cognitive change identified in twins, higher baseline levels of two sphingomyelins predict a slower decline in executive function, and one of these associations was replicated in WRAP. This finding is in agreement with findings from previous epidemiologic study among older adults (age 65+) (Mielke et al., 2010, 2011; Han et al., 2011; Orešič et al., 2011). The relationship between sphingomyelins and cognitive function was also consistent with previous research demonstrating that a disrupted sphingomyelin pathway is involved in Aβ pathology (Jana and Pahan, 2004), and a reduced level of sphingomyelin in human AD brain than that in cognitively healthy individuals (He et al., 2010). We were unable to replicate the other three metabolites in the WRAP, probably due to the small sample size, the heterogeneous clinical phenotypes across different studies, or other differences such as WRAP included women as well and participants who were enriched with a parental history of AD. The lack of replication could also be attributed to false positives as a result of the small sample size in our discovery cohort.

We also found a higher baseline level of plasma DAG (18:2/20:4) was associated with a slower decrease in short-term memory. While the precise mechanism underlying this association is unclear and awaits further investigation, the association appears to be consistent with previous studies showing that neuronal glycerolipid biosynthesis promotes axon regeneration after optic nerve injuries (Yang et al., 2020), and that diffuse axonal injury causes cognitive impairment including short-term memory (Wallesch et al., 2001). In addition, a higher plasma level of hydroxy-CMPF was associated with a slower decrease in general cognitive ability. The CMPF is a furan fatty acid that inhibits insulin secretion in both human and mouse models (Nagy et al., 2017). A higher level of CMPF has been associated with type 2 diabetes (Prentice et al., 2014) and biological aging (Menni et al., 2013) in human. While these findings are seemingly contradictory to the protective effect of hydroxy-CMPF on cognitive function identified in the current study, it is possible that hydroxylation of CMPF may reduce the plasma concentration of CMPF, which could likely improve neuronal insulin resistance, thereby exerting beneficial effects on cognition.

Our network analysis revealed that loss of metabolites connectivity (the pink module) may be related to cognitive change, and we identified several hub metabolites that differ significantly between MCI converters and non-converters. Although the precise mechanisms are unclear, it is likely that the differences in network connectivity and hub metabolites may explain, at least in part, why some twins progressed to MCI while his co-twin brother did not. For example, we found that 2,2′-methylenebis(6-tert-butyl-p-cresol), ceramide (d16:1/24:1, d18:1/22:1), and gamma-glutamylglutamate were the hub metabolites only in non-converters but not in converters. Ceramides are sphingolipids that play important roles in cellular processes including differentiation, proliferation, and apoptosis (Hannun and Obeid, 2008). Dysregulation of ceramides in the brain has been implied in Alzheimer’s and other neurodegenerative disease (Filippov et al., 2012; Jazvinšćak Jembrek et al., 2015). Gamma-glutamylglutamate acid is a dipeptide consisting of gamma-glutamate and glutamic acid. Glutamic acid is the main excitatory neurotransmitter in the brain, which can be released to neurons and form glutamate; the latter can convert to glutamic acid via astrocytes, which completes the glutamate-glutamine cycle (Daikhin and Yudkoff, 2000). Dysregulation of glutamate signaling has been implicated in neurodegeneration and cognitive behaviors (Vogel et al., 1966; Rahn et al., 2012).

Among the metabolic pathways we identified to be associated with cognitive changes, the aminoacyl-tRNA biosynthesis pathway showed the most significant association with changes in short-term memory. The process of aminoacyl tRNA biosynthesis catalyzes the ligation of amino acids to their cognate tRNAs, and dysregulation of this pathway has been associated with neurodegeneration (Park et al., 2008; Kapur et al., 2017). The observed association of the amino acid (valine, leucine, and isoleucine) pathway with short-term memory is consistent with previous research demonstrating that branched-chain amino acids (BCAAs) were associated with dementia (Tynkkynen et al., 2018) and metabolic conditions such as insulin resistance and diabetes (Roberts et al., 2014), both of which have been linked to cognitive impairment (Ma et al., 2015) and AD (Arnold et al., 2018).

In this study, we identified associations of circulating plasma metabolites with cognitive change over time, independent of known risk factors. The mechanisms underlying the observed associations are unclear but believed to be highly complex. While genome-wide association studies (GWAS) have identified many genetic variants associated with cognitive functions, cognitive impairment and AD in various human populations under different clinical settings (Trampush et al., 2017; Jansen et al., 2019), the identified genetic variants explain only a small proportion of the disease risk. A larger contributor to the phenotypic variation may be due to the complex interactions between environmental and inherited factors. Metabolomics provides a snapshot of hundreds to thousands of small molecules (metabolites) in a biological sample (e.g., plasma) at a given time. The abundance of these metabolites is influenced by both genetic and environmental factors, such as diet, lifestyle, behavior, drug, microbiome, hormones, etc., probably through various epigenetic mechanisms. Metabolomic study of circulating metabolites could facilitate the identification of environmental factors contributing to cognition changes during aging, and provides mechanistic insight into how environmental exposures contribute to cognitive impairment in AD progression. Moreover, we identified that different plasma lipids were associated with different cognitive domains, including general cognitive ability, short term memory, executive memory, but not episodic memory. While the mechanism underlying this differential associations remain to be determined, the lack of association with episodic memory might suggest that the underlying pathological mechanisms are dementia-related but not necessarily specific to AD. However, it is worth noting that the VETSA and WRAP samples are both relatively young. Although prospective studies of AD usually emphasize episodic memory, there is also evidence suggesting that executive function can predict progression to MCI or AD as early or possibly earlier than memory deficits (Clark et al., 2012; Ewers et al., 2012; Gibbons et al., 2012). Executive function deficits, which were associated with metabolites in the present study, are also common in MCI and AD, and executive control functions do affect functioning in other cognitive domains (Buckner, 2004; Baudic et al., 2006). Indeed, there is also evidence for an executive-prominent AD subtype (Mukherjee et al., 2012). Thus, although it remains uncertain, it is possible that the associated deficits that were observed will contribute to AD-related deficits.

Limitations of our study include the small sample size and the inclusion of only Caucasian men in the twin study, the sample differences across studies such as the selection of WRAP biased toward parental history of AD, and the different cognitive measures in our discovery and replication cohorts. The uncertainty of the unknown metabolites represents another limitation. However, our study has several strengths, including a focus on midlife, the well-matched MCI discordant twin pairs (which eliminates confounding by age and sex and minimized confounding by multiple background factors), the inclusion of a replication cohort comprised of middle-aged adults in a similar age range, the prospective design of both discovery and replication cohorts, the comprehensive cognitive phenotypes, and the comprehensive coverage of plasma metabolome by using an untargeted metabolomics approach on the same platform in both cohorts.

Conclusion

In summary, in an untargeted metabolomics analysis of middle-aged adults in two well characterized prospective cohorts, we identified four plasma metabolites that are significantly predictive of cognitive change over time, after adjusting for multiple testing and known clinical covariates. The observed association of a sphingomyelin (SM [d18:1/20:1, d18:2/20:0]) with executive function was replicated in an independent cohort including participants in a similar age range. We also identified metabolic networks and hub metabolites that may be involved in the relationship between metabolic dysregulation and cognitive impairment in midlife. Our results support that metabolic perturbation occurs many years before cognitive impairment, and suggest the possibility of using blood metabolites as non-invasive biomarkers in early prediction of cognitive aging and AD in midlife.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: The authors will accept reasonable requests for data access, which will be referred to the Veteran Era Twin Registry where the data are maintained and overseen. Requests to access these datasets should be directed to http://www.vetsatwins.org.

Ethics Statement

The studies involving human participants were reviewed and approved by twins participating in the VETSA were administered identical protocols at the University of California, San Diego, or Boston University. Participants in the WRAP were administered at the University of Wisconsin. Written consents were obtained from all participants enrolled in both studies. The protocols were approved by Institutional Review Boards at the participating universities. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

ZH, BR, and JE performed the statistical analysis. RD, CE, and SJ performed the replication analysis in WRAP. ML, CF, WK, and JZ designed and coordinated the study. ZH and JZ wrote the manuscript. All authors reviewed the manuscript.

Funding

The Cooperative Studies Program of the United States Department of Veterans Affairs provided financial support for development and maintenance of the Vietnam Era Twin Registry (VETSA). The VETSA was supported by National Institute on Aging R01s AG022381 and AG050595. WK received additional support from the VA San Diego Center of Excellence for Stress and Mental Health. The WRAP was supported by National Institute on Aging R01s AG027161 and AG054047 and computational resources were supported by a core grant to the Center for Demography and Ecology at the University of Wisconsin–Madison (P2CHD047873). JZ and ZH were supported by National Institute on Aging grants RF1AG052476 and R01DK107532. The authors would also like to thank all participants for their participation to both studies. The content is the responsibility of the authors and does not necessarily represent official views of the NIA, NIH, or VA.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer LY declared a past co-authorship with several of the authors ZH and JZ to the handling editor.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2020.555850/full#supplementary-material

References

Arnold, S. E., Arvanitakis, Z., Macauley-Rambach, S. L., Koenig, A. M., Wang, H.-Y., Ahima, R. S., et al. (2018). Brain insulin resistance in type 2 diabetes and Alzheimer disease: concepts and conundrums. Nat. Rev. Neurol. 14:168. doi: 10.1038/nrneurol.2017.185

PubMed Abstract | CrossRef Full Text | Google Scholar

Bates, D., Maechler, M., Bolker, B., and Walker, S. (2014). lme4: Linear Mixed-Effects Models using Eigen and S4. R Package Version 1, 1–23.

Google Scholar

Baudic, S., Barba, G. D., Thibaudet, M. C., Smagghe, A., Remy, P., and Traykov, L. (2006). Executive function deficits in early Alzheimer’s disease and their relations with episodic memory. Arch. Clin. Neuropsychol. 21, 15–21. doi: 10.1016/j.acn.2005.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Björkqvist, M., Ohlsson, M., Minthon, L., and Hansson, O. (2012). Evaluation of a previously suggested plasma biomarker panel to identify Alzheimer’s disease. PLoS One 7:e29868. doi: 10.1371/journal.pone.0029868

PubMed Abstract | CrossRef Full Text | Google Scholar

Bressler, J., Yu, B., Mosley, T. H., Knopman, D. S., Gottesman, R. F., Alonso, A., et al. (2017). Metabolomics and cognition in African American adults in midlife?: the atherosclerosis risk in communities study. Transl. Psychiatry 7:e1173. doi: 10.1038/tp.2017.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L. (2004). Memory and executive function in aging and AD: multiple factors that cause decline and reserve factors that compensate. Neuron 44, 195–208. doi: 10.1016/j.neuron.2004.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Casanova, R., Varma, S., Simpson, B., Kim, M., An, Y., Saldana, S., et al. (2016). Blood metabolite markers of preclinical Alzheimer’s disease in two longitudinally followed cohorts of older individuals. Alzheimers Dement. 12, 815–822. doi: 10.1016/j.jalz.2015.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Chong, J., Wishart, D. S., and Xia, J. (2019). Using metaboanalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr. Protoc. Bioinformatics 68:e86.

Google Scholar

Chouraki, V., Preis, S. R., Yang, Q., Beiser, A., Li, S., Vasan, R. S., et al. (2017). Association of amine biomarkers with incident dementia and Alzheimer’s disease in the Framingham Study. Alzheimers Dement. 13, 1327–1336. doi: 10.1016/j.jalz.2017.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Clark, L. R., Racine, A. M., Koscik, R. L., Okonkwo, O. C., Engelman, C. D., Carlsson, C. M., et al. (2016). Beta-amyloid and cognitive decline in late middle age: findings from the Wisconsin Registry for Alzheimer’s prevention study. Alzheimers Dement. 12, 805–814. doi: 10.1016/j.jalz.2015.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Clark, L. R., Schiehser, D. M., Weissberger, G. H., Salmon, D. P., Delis, D. C., and Bondi, M. W. (2012). Specific measures of executive function predict cognitive decline in older adults. J. Int. Neuropsychol. Soc. 18, 118–127. doi: 10.1017/s1355617711001524

PubMed Abstract | CrossRef Full Text | Google Scholar

Clarke, J. R., Ribeiro, F. C., Frozza, R. L., De Felice, F. G., and Lourenco, M. V. (2018). Metabolic dysfunction in Alzheimer’s disease: from basic neurobiology to clinical approaches. J. Alzheimers Dis. 64, S405–S426.

Google Scholar

Daikhin, Y., and Yudkoff, M. (2000). Compartmentation of brain glutamate metabolism in neurons and glia. J. Nutr. 130, 1026S–1031S.

Google Scholar

Darst, B. F., Jonaitis, E. M., Koscik, R. L., Clark, L. R., Lu, Q., Hogan, K. J., et al. (2018). Metabolites associated with early cognitive changes implicated in Alzheimer’s disease. bioRxiv [Preprint] doi: 10.1101/436667

CrossRef Full Text | Google Scholar

Darst, B. F., Lu, Q., Johnson, S. C., and Engelman, C. D. (2019). Integrated analysis of genomics, longitudinal metabolomics, and Alzheimer’s risk factors among 1,111 cohort participants. Genet. Epidemiol. 43, 657–674.

Google Scholar

Dowling, N. M., Hermann, B., La Rue, A., and Sager, M. A. (2010). Latent structure and factorial invariance of a neuropsychological test battery for the study of preclinical Alzheimer’s disease. Neuropsychology 24, 742–756. doi: 10.1037/a0020176

PubMed Abstract | CrossRef Full Text | Google Scholar

Evans, A. M., Bridgewater, B. R., Liu, Q., Mitchell, M. W., Robinson, R. J., Dai, H., et al. (2014). High resolution mass spectrometry improves data quantity and quality as compared to unit mass resolution mass spectrometry in high-throughput profiling metabolomics. Metabolomics 4:1. doi: 10.4172/2153-0769.1000132

CrossRef Full Text | Google Scholar

Evans, A. M., DeHaven, C. D., Barrett, T., Mitchell, M., and Milgram, E. (2009). Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal. Chem. 81, 6656–6667. doi: 10.1021/ac901536h

PubMed Abstract | CrossRef Full Text | Google Scholar

Ewers, M., Walsh, C., Trojanowski, J. Q., Shaw, L. M., Petersen, R. C., and Jack, C. R. Jr., et al. (2012). Prediction of conversion from mild cognitive impairment to Alzheimer’s disease dementia based upon biomarkers and neuropsychological test performance. Neurobiol. Aging 33, 1203–1214. doi: 10.1016/j.neurobiolaging.2010.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Filippov, V., Song, M. A., Zhang, K., Vinters, H. V., Tung, S., Kirsch, W. M., et al. (2012). Increased ceramide in brains with Alzheimer’s and other neurodegenerative diseases. J. Alzheimers Dis. 29, 537–547. doi: 10.3233/jad-2011-111202

PubMed Abstract | CrossRef Full Text | Google Scholar

Franz, C. E., Lyons, M. J., O’Brien, R., Panizzon, M. S., Kim, K., Bhat, R., et al. (2011). A 35-year longitudinal assessment of cognition and midlife depression symptoms: the vietnam era twin study of aging. Am. J. Geriatr. Psychiatry 19, 559–570. doi: 10.1097/jgp.0b013e3181ef79f1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibbons, L. E., Carle, A. C., Mackin, R. S., Harvey, D., Mukherjee, S., Insel, P., et al. (2012). A composite score for executive functioning, validated in Alzheimer’s Disease Neuroimaging Initiative (ADNI) participants with baseline mild cognitive impairment. Brain Imaging Behav. 6, 517–527. doi: 10.1007/s11682-012-9176-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Granholm, E. L., Panizzon, M. S., Elman, J. A., Jak, A. J., Hauger, R. L., Bondi, M. W., et al. (2017). Pupillary responses as a biomarker of early risk for Alzheimer’s disease. J. Alzheimers Dis. 56, 1419–1428. doi: 10.3233/jad-161078

PubMed Abstract | CrossRef Full Text | Google Scholar

Gustavson, D. E., Panizzon, M. S., Franz, C. E., Friedman, N. P., Reynolds, C. A., Jacobson, K. C., et al. (2018). Genetic and environmental architecture of executive functions in midlife. Neuropsychology 32, 18–30. doi: 10.1037/neu0000389

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, X., Rozen, S., Boyle, S. H., Hellegers, C., Cheng, H., Burke, J. R., et al. (2011). Metabolomics in early Alzheimer’s disease: identification of altered plasma sphingolipidome using shotgun lipidomics. PLoS One 6:e21643. doi: 10.1371/journal.pone.0021643

PubMed Abstract | CrossRef Full Text | Google Scholar

Hannun, Y. A., and Obeid, L. M. (2008). Principles of bioactive lipid signalling: lessons from sphingolipids. Nat. Rev. Mol. Cell Biol. 9, 139–150. doi: 10.1038/nrm2329

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrigan, G. G., and Goodacre, R. (2012). Metabolic Profiling: Its Role in Biomarker Discovery and Gene Function Analysis. Berlin: Springer Science & Business Media.

Google Scholar

He, X., Huang, Y., Li, B., Gong, C.-X., and Schuchman, E. H. (2010). Deregulation of sphingolipid metabolism in Alzheimer’s disease. Neurobiol. Aging 31, 398–408. doi: 10.1016/j.neurobiolaging.2008.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Huo, Z., Yu, L., Yang, J., Zhu, Y., Bennett, D. A., and Zhao, J. (2020). Brain and blood metabolome for Alzheimer’s dementia: findings from a targeted metabolomics analysis. Neurobiol. Aging 86, 123–133. doi: 10.1016/j.neurobiolaging.2019.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Inoue, K., Tsuchiya, H., Takayama, T., Akatsu, H., Hashizume, Y., Yamamoto, T., et al. (2015). Blood-based diagnosis of Alzheimer’s disease using fingerprinting metabolomics based on hydrophilic interaction liquid chromatography with mass spectrometry and multivariate statistical analysis. J. Chromatogr. B 974, 24–34. doi: 10.1016/j.jchromb.2014.10.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Jak, A. J., Bondi, M. W., Delano-Wood, L., Wierenga, C., Corey-Bloom, J., Salmon, D. P., et al. (2009). Quantification of five neuropsychological approaches to defining mild cognitive impairment. Am. J. Geriatr. Psychiatry 17, 368–375. doi: 10.1097/jgp.0b013e31819431d5

PubMed Abstract | CrossRef Full Text | Google Scholar

Jana, A., and Pahan, K. (2004). Fibrillar amyloid-$β$ peptides kill human primary neurons via NADPH oxidase-mediated activation of neutral sphingomyelinase implications For Alzheimer’s disease. J. Biol. Chem. 279, 51451–51459. doi: 10.1074/jbc.m404635200

PubMed Abstract | CrossRef Full Text | Google Scholar

Jansen, I. E., Savage, J. E., Watanabe, K., Bryois, J., Williams, D. M., Steinberg, S., et al. (2019). Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat. Genet. 51, 404–413.

Google Scholar

Jazvinšćak Jembrek, M., Hof, P. R., and Šimić, G. (2015). Ceramides in Alzheimer’s disease: key mediators of neuronal apoptosis induced by oxidative stress and A$β$ accumulation. Oxid. Med. Cell. Longev. 2015:346783.

Google Scholar

Johnson, S. C., Koscik, R. L., Jonaitis, E. M., Clark, L. R., Mueller, K. D., Berman, S. E., et al. (2018). The Wisconsin Registry for Alzheimer’s Prevention: a review of findings and current directions. Alzheimers Dement. 10, 130–142. doi: 10.1016/j.dadm.2017.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Kapur, M., Monaghan, C. E., and Ackerman, S. L. (2017). Regulation of mRNA translation in neurons—a matter of life and death. Neuron 96, 616–637. doi: 10.1016/j.neuron.2017.09.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Karlamangla, A. S., Lachman, M. E., Han, W., Huang, M., and Greendale, G. A. (2017). Evidence for cognitive aging in midlife women: study of women’s health across the nation. PLoS One 12:e0169008. doi: 10.1371/journal.pone.0169008

PubMed Abstract | CrossRef Full Text | Google Scholar

Koscik, R. L., La Rue, A., Jonaitis, E. M., Okonkwo, O. C., Johnson, S. C., Bendlin, B. B., et al. (2014). Emergence of mild cognitive impairment in late middle-aged adults in the wisconsin registry for Alzheimer’s prevention. Dement. Geriatr. Cogn. Dis. 38, 16–30. doi: 10.1159/000355682

PubMed Abstract | CrossRef Full Text | Google Scholar

Kremen, W. S., Beck, A., Elman, J. A., Gustavson, D. E., Reynolds, C. A., Tu, X. M., et al. (2019). Influence of young adult cognitive ability and additional education on later-life cognition. Proc. Natl. Acad. Sci. U.S.A. 116, 2021–2026. doi: 10.1073/pnas.1811537116

PubMed Abstract | CrossRef Full Text | Google Scholar

Kremen, W. S., Franz, C. E., and Lyons, M. J. (2013a). VETSA: the vietnam era twin study of aging. Twin Res. Hum. Genet. 16, 399–402. doi: 10.1017/thg.2012.86

PubMed Abstract | CrossRef Full Text | Google Scholar

Kremen, W. S., Jak, A. J., Panizzon, M. S., Spoon, K. M., Franz, C. E., Thompson, W. K., et al. (2013b). Early identification and heritability of mild cognitive impairment. Int. J. Epidemiol. 43, 600–610.

Google Scholar

Kremen, W. S., Panizzon, M. S., Franz, C. E., Spoon, K. M., Jacobson, K. C., Vasilopoulos, T., et al. (2014). Genetic complexity of episodic memory?: a twin approach to studies of aging. Psychol. Aging 29, 404–417. doi: 10.1037/a0035962

PubMed Abstract | CrossRef Full Text | Google Scholar

Kremen, W. S., Panizzon, M. S., Xian, H., Barch, D. M., Franz, C. E., Grant, M. D., et al. (2011). Genetic architecture of context processing in late middle age: more than one underlying mechanism. Psychol. Aging 26, 852–863. doi: 10.1037/a0025098

PubMed Abstract | CrossRef Full Text | Google Scholar

Kremen, W. S., Thompson-Brenner, H., Leung, Y.-M. J., Grant, M. D., Franz, C. E., Eisen, S. A., et al. (2006). Genes, environment, and time: the Vietnam era twin study of aging (VETSA). Twin Res. Hum. Genet. 9, 1009–1022. doi: 10.1375/twin.9.6.1009

PubMed Abstract | CrossRef Full Text | Google Scholar

La Rue, A., Hermann, B., Jones, J. E., Johnson, S., Asthana, S., and Sager, M. A. (2008). Effect of parental family history of Alzheimer’s disease on serial position profiles. Alzheimers Dement. 4, 285–290. doi: 10.1016/j.jalz.2008.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9:559. doi: 10.1186/1471-2105-9-559

PubMed Abstract | CrossRef Full Text | Google Scholar

Lyons, M. J., Panizzon, M. S., Liu, W., McKenzie, R., Bluestone, N. J., Grant, M. D., et al. (2017). A longitudinal twin study of general cognitive ability over four decades. Dev. Psychol. 53, 1170–1177. doi: 10.1037/dev0000303

PubMed Abstract | CrossRef Full Text | Google Scholar

Lyons, M. J., York, T. P., Franz, C. E., Grant, M. D., Eaves, L. J., Jacobson, K. C., et al. (2009). Genes determine stability and the environment determines change in cognitive ability during 35 years of adulthood. Psychol. Sci. 20, 1146–1152. doi: 10.1111/j.1467-9280.2009.02425.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, L., Wang, J., and Li, Y. (2015). Insulin resistance and cognitive dysfunction. Clin. Chim. Acta 444, 18–23. doi: 10.1016/j.cca.2015.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Mapstone, M., Cheema, A. K., Fiandaca, M. S., Zhong, X., Mhyre, T. R., MacArthur, L. H., et al. (2014). Plasma phospholipids identify antecedent memory impairment in older adults. Nat. Med. 20, 415–418. doi: 10.1038/nm.3466

PubMed Abstract | CrossRef Full Text | Google Scholar

Margolin, A. A., Nemenman, I., Basso, K., Wiggins, C., Stolovitzky, G., Dalla Favera, R., et al. (2006). ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics 7(Suppl. 1):S7. doi: 10.1186/1471-2105-7-S1-S7

PubMed Abstract | CrossRef Full Text | Google Scholar

McKenzie, A. T., Katsyv, I., Song, W.-M., Wang, M., and Zhang, B. (2016). DGCA: a comprehensive R package for differential gene correlation analysis. BMC Systems Biol. 10:106. doi: 10.1186/s12918-016-0349-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Menni, C., Kastenmüller, G., Petersen, A. K., Bell, J. T., Psatha, M., Tsai, P.-C., et al. (2013). Metabolomic markers reveal novel pathways of ageing and early development in human populations. Int. J. Epidemiol. 42, 1111–1119. doi: 10.1093/ije/dyt094

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, P. E., Lafitte, F., and Bontempi, G. (2008). Minet: an open source R/Bioconductor package for mutual information based network inference. BMC Bioinformatics 9:461.

Google Scholar

Mielke, M. M., Bandaru, V. V. R., Haughey, N. J., Rabins, P. V., Lyketsos, C. G., and Carlson, M. C. (2010). Serum sphingomyelins and ceramides are early predictors of memory impairment. Neurobiol. Aging 31, 17–24. doi: 10.1016/j.neurobiolaging.2008.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Mielke, M. M., Haughey, N. J., Bandaru, V. V. R., Weinberg, D. D., Darby, E., Zaidi, N., et al. (2011). Plasma sphingomyelins are associated with cognitive progression in Alzheimer’s disease. J. Alzheimers Dis. 27, 259–269.

Google Scholar

Mukherjee, S., Trittschuh, E., Gibbons, L. E., Mackin, R. S., Saykin, A., Crane, P. K., et al. (2012). Dysexecutive and amnesic AD subtypes defined by single indicator and modern psychometric approaches: relationships with SNPs in ADNI. Brain Imaging Behav. 6, 649–660. doi: 10.1007/s11682-012-9207-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagy, E., Liu, Y., Prentice, K. J., Sloop, K. W., Sanders, P. E., Batchuluun, B., et al. (2017). Synthesis and characterization of urofuranoic acids: in vivo metabolism of 2-(2-carboxyethyl)-4-methyl-5-propylfuran-3-carboxylic acid (cmpf) and effects on in vitro insulin secretion. J. Med. Chem. 60, 1860–1875. doi: 10.1021/acs.jmedchem.6b01668

PubMed Abstract | CrossRef Full Text | Google Scholar

Orešič, M., Hyötyläinen, T., Herukka, S. K., Sysi-Aho, M., Mattila, I., Seppänan-Laakso, T., et al. (2011). Metabolome in progression to Alzheimer’s disease. Transl. Psychiatry 1:e57.

Google Scholar

Orme, D. R., Brehm, W., and Ree, M. J. (2001). Armed forces qualification test as a measure of premorbid intelligence. Mil. Psychol. 13, 187–197. doi: 10.1207/s15327876mp1304_1

CrossRef Full Text | Google Scholar

Park, S. G., Schimmel, P., and Kim, S. (2008). Aminoacyl tRNA synthetases and their connections to disease. Proc. Natl. Acad. Sci. U.S.A. 105, 11043–11049.

Google Scholar

Prentice, K. J., Luu, L., Allister, E. M., Liu, Y., Jun, L. S., Sloop, K. W., et al. (2014). The furan fatty acid metabolite CMPF is elevated in diabetes and induces $β$ cell dysfunction. Cell Metabol. 19, 653–666. doi: 10.1016/j.cmet.2014.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Proitsi, P., Kuh, D., Wong, A., Maddock, J., Bendayan, R., Wulaningsih, W., et al. (2018). Lifetime cognition and late midlife blood metabolites: findings from a British birth cohort. Transl. Psychiatry 8:203.

Google Scholar

Rahn, K. A., Slusher, B. S., and Kaplin, A. I. (2012). Glutamate in CNS neurodegeneration and cognition and its regulation by GCPII inhibition. Curr. Med. Chem. 19, 1335–1345. doi: 10.2174/092986712799462649

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberts, L. D., Koulman, A., and Griffin, J. L. (2014). Towards metabolic biomarkers of insulin resistance and type 2 diabetes: progress from the metabolome. Lancet Diabetes Endocrinol. 2, 65–75. doi: 10.1016/s2213-8587(13)70143-8

CrossRef Full Text | Google Scholar

Rochfort, S. (2005). Metabolomics reviewed: a new “omics” platform technology for systems biology and implications for natural products research. J. Nat. Prod. 68, 1813–1820. doi: 10.1021/np050255w

PubMed Abstract | CrossRef Full Text | Google Scholar

Sager, M. A., Hermann, B., and La Rue, A. (2005). Middle-aged children of persons with Alzheimer’s disease: APOE genotypes and cognitive function in the Wisconsin Registry for Alzheimer’s Prevention. J. Geriatr. Psychiatry Neurol. 18, 245–249. doi: 10.1177/0891988705281882

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, M. (1996). Rey Auditory Verbal Learning Test: A Handbook. Los Angeles, CA: Western Psychological Services.

Google Scholar

Schoenborn, C. A., and Heyman, K. M. (2009). Health characteristics of adults aged 55 years and over: United States, 2004–2007. Natl. Health Stat. Report 16, 1–31.

Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh-Manoux, A., Kivimaki, M., Glymour, M. M., Elbaz, A., Berr, C., Ebmeier, K. P., et al. (2012). Timing of onset of cognitive decline: results from Whitehall II prospective cohort study. BMJ 344, d7622. doi: 10.1136/bmj.d7622

PubMed Abstract | CrossRef Full Text | Google Scholar

Trampush, J. W., Yang, M. L. Z., Yu, J., Knowles, E., Davies, G., Liewald, D. C., et al. (2017). GWAS meta-analysis reveals novel loci and genetic correlates for general cognitive function: a report from the COGENT consortium. Mol. Psychiatry 22, 336–345.

Google Scholar

Tynkkynen, J., Chouraki, V., van der Lee, S. J., Hernesniemi, J., Yang, Q., Li, S., et al. (2018). Association of branched-chain amino acids and other circulating metabolites with risk of incident dementia and Alzheimer’s disease: a prospective study in eight cohorts. Alzheimers Dement. 14, 723–733. doi: 10.1016/j.jalz.2018.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Uhlaner, J. E., and Bolanovich, D. J. (1952). Development of Armed Forces Qualification Test and Predecessor Army Screening Tests. 1946–1950. Washington, DC: Personnel Research Section, Department of the Army.

Google Scholar

van den Berg, R. A., Hoefsloot, H. C. J., Westerhuis, J. A., Smilde, A. K., and van der Werf, M. J. (2006). Centering, scaling, and transformations: improving the biological information content of metabolomics data. BMC Genomics 7:142. doi: 10.1186/1471-2164-7-142

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Lee, S. J., Teunissen, C. E., Shipley, M. J., Teumer, A., Chouraki, V., Melo, D., et al. (2018). Circulating metabolites and general cognitive ability and dementia?: evidence from 11 cohort studies. Alzheimers Dement. 14, 707–722. doi: 10.1016/j.jalz.2017.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Varma, V. R., Oommen, A. M., Varma, S., Casanova, R., An, Y., Andrews, R. M., et al. (2018). Brain and blood metabolite signatures of pathology and progression in Alzheimer disease: a targeted metabolomics study. PLoS Med. 15:e1002482. doi: 10.1371/journal.pmed.1002482

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogel, W., Broverman, D. M., Draguns, J. G., and Klaiber, E. L. (1966). The role of glutamic acid in cognitive behaviors. Psychol. Bull. 65, 367–382. doi: 10.1037/h0023351

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallesch, C.-W., Curio, N., Kutz, S., Jost, S., Bartels, C., and Synowitz, H. (2001). Outcome after mild-to-moderate blunt head injury: effects of focal lesions and diffuse axonal injury. Brain Inj. 15, 401–412. doi: 10.1080/02699050116886

CrossRef Full Text | Google Scholar

Yang, C., Wang, X., Wang, J., Siniossoglou, S., Yao, Z., Liu, K., et al. (2020). Rewiring neuronal glycerolipid metabolism determines the extent of axon regeneration. Neuron 105, 276.e5–292.e5. doi: 10.1016/j.neuron.2019.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar