Research Paper Volume 8, Issue 12 pp 3241—3254

Serum protein mediators of dementia and aging proper

Donald R. Royall1,2,3,4, , Safa Al-Rubaye1, , Ram Bishnoi5, , Raymond F. Palmer3, ,

  • 1 Department of Psychiatry, University of Texas Health Science Center at San Antonio, San Antonio, TX 78229, USA
  • 2 Department of Medicine, University of Texas Health Science Center at San Antonio, San Antonio, TX 78229, USA
  • 3 Department of Family and Community Medicine, University of Texas Health Science Center at San Antonio, San Antonio, TX 78229, USA
  • 4 South Texas Veterans’ Health System Audie L. Murphy Division GRECC, San Antonio, TX 78229, USA
  • 5 Department of Psychiatry, The Medical College of Georgia, Augusta, GA 30912, USA

Received: August 31, 2016       Accepted: November 17, 2016       Published: December 3, 2016      

https://doi.org/10.18632/aging.101091
How to Cite

Abstract

The latent variable “δ” (for “dementia”) appears to be uniquely responsible for the dementing aspects of cognitive impairment. Age, depressive symptoms, gender and the apolipoprotein E (APOE) ε4 allele are independently associated with δ. In this analysis, we explore serum proteins as potential mediators of age’s specific association with δ in a large, ethnically diverse longitudinal cohort, the Texas Alzheimer’s Research and Care Consortium (TARCC). 22 serum proteins were recognized as partial mediators of age’s association with δ. These include Insulin-like Growth Factor-Binding Protein 2 (IGF-BP2), which we had previously associated with age-specific cognitive change, and both Pancreatic Polypeptide (PP) and von Willebrand Factor (vWF), previously associated with δ. Nine other δ-related proteins were not confirmed by this ethnicity adjusted analysis. Our findings suggest that age’s association with the disabling fraction of cognitive performance is partially mediated by serum proteins, somatomedins and hormones. Those proteins may offer targets for the specific treatment of age-related effects on dementia severity and conversion risk.

Introduction

Age, depression, and the apolipoprotein E (APOE) ε4 allele are independently associated with the latent dementia phenotype “δ” (for “dementia”) [1]. Their associations with dementia do not necessarily involve neurodegeneration. Depression’s association with cognitive decline in older persons is not mediated by neurodegenerative changes [2], while age’s association with δ has been shown to be fully mediated by a paucity of neurodegenerative changes in pathologically confirmed Alzheimer’s Disease (AD) cases [3]. Brain aging is therefore not AD [4].

On the other hand, clinical “AD” may very well have an aging component. Since δ is essentially the sole cognitive determinant of dementia severity, clinical dementia must arise from the sum of all independent δ-related processes. Age’s small independent effect appears to be linear over the lifespan, and cumulative [5]. Over a 50 year age range, aging might account for up to a standard deviation change in composite “d-scores”. That is not trivial. δ’s intercept and slope are uniquely strong determinants of future dementia status [6,7]. Each quintile in the d-score distribution of non-demented persons increases conversion to clinical "AD" by 50% [three-fold among "Mild Cognitive Impairment (MCI)" cases] [8].

In the “oldest old”, aging alone might sum with comorbid neurodegenerative processes to push d-scores into their demented range. This should effectively reduce the amount of neuropathology required to make the diagnosis of dementia in centenarians, and modulate the apparent associations between various neuro-pathologies and clinical dementia. In fact, dementia at advanced age is associated with lower levels of AD-specific neuropathology [9], and less widely spread pathology [10].

Even in their aggregate, demographic-specific dementia risks explain a minority of δ’s variance [1]. Thus, regardless of whether age’s effect is mediated by neurodegeneration, observed dementia status must be largely determined by age-independent factors. In Age, depression and APOE adjusted models, we have found the majority of δ’s remaining variance to be associated with a large number of pro- and anti-inflammatory serum protein biomarkers [1,1113].

On the other hand, we have reported serum Insulin-like Binding Protein 2 (IGF-BP2) to be a strong correlate of age’s specific cognitive effects [14]. However, age has both direct (δ-independent) and indirect (δ-related) effects on cognition. It has yet to be determined whether IGF-BP2 mediates age’s association with δ, or its δ-independent direct effects instead.

In this analysis, we combine structural equation models (SEM) with longitudinal data from the Texas Alzheimer’s Research and Care Consortium (TARCC) to explore more than 100 serum proteins as potential mediators of age’s specific association with δ. Our models are constructed such that the significant mediators of age’s effect on prospective δ scores can be interpreted causally. The mediators should offer both insights into the pathophysiology of Aging Proper, and potential targets for the remediation of age-specific cognitive impairments.

Results

The demographic characteristics of our sample are presented in Table 1. The ethnicity equivalent unadjusted Visit 2 δ homolog composite score (i.e., “dEQ”) achieved a high Area Under the Receiver Operating Characteristic curve (AUC/ROC) for the discrimination between AD cases and normal controls (NC) (AUC = 0.953; CI: 0.946-0.960). “g’”’s (i.e., δ’s residual in Spearman’s general intelligence factor “g”) AUC for the same discrimination was at a near chance level [AUC = 0.536 (CI: 0.514-0.558)]. This is consistent with past findings, across batteries, in this and other cohorts.

Table 1. Descriptive Statistics.

VariableNMean (SD)
Age (observed)338170.88 (9.48)
APOE e4 alleles (1 = e4+, n = 1223)31540.39 (0.49)
CDR (Sum of Boxes)33062.42 (3.35)
COWA33818.41 (3.49)
DIS33818.89 (3.01)
EDUC (observed)338113.24 (4.25)
Ethnicity (1 = MA, n = 1189)33810.36 (0.47)
GDS30 (observed)30055.60 (5.25)
Gender (♂ = 1, n = 1281)33120.39 (0.49)
IADL (Summed)338110.48 (4.52)
MMSE331125.52 (4.76)
WMS LM II33818.05 (4.30)
WMS VR I33817.88 (3.68)
Complete Cases2861
CDR = Clinical Dementia Rating scale; COWA = Controlled Oral Word Association Test; DIS = Digit Span Test; GDS = Geriatric Depression Scale [66]; IADL = Instrumental Activities of Daily Living; MMSE = Mini-Mental State Exam [67]; SD = standard deviation; WMS LM II = Weschler Memory Scale: Delayed Logical Memory; WMS VR I = Weschler Memory Scale: Immediate Visual Reproduction.

The Base Model had excellent fit [χ2 = 5.59 (13), p = 0.960; CFI = 1.00; RMSEA = 0.00]. Independently of the covariates [i.e., APOE ε4 allelic burden, depressive symptoms, education, ethnicity, gender, homocysteine (HCY), and hemoglobin A1c (HgbA1c)], baseline age was significantly directly associated with Visit 2 dEQ (r = -0.25, p<0.001), and weakly with the Visit 2 g’ composite (r = -0.11, p ≤ 0.001). Age’s significant association with Visit 2 dEQ scores was in a negative direction suggesting an adverse effect on observed cognitive performance.

The mediation models all had acceptable fit [e.g., IGF-BP2: χ2 = 387.90 (17), p < 0.001; CFI = 0.927; RMSEA = 0.044 (Figure 1)]. 22 proteins achieved statistically significant mediation effects after Bonferroni correction for multiple comparisons (Table 2). IGF-BP2 had previously been recognized as an age-specific serum protein biomarker [14]. Pancreatic Polypeptide (PP) and von Willebrand Factor (vWF) had previously been recognized as δ-related serum protein biomarkers [13]. Table 3 presents the mediation effects. All the identified proteins were partial mediators, but several had relatively large effects (range 9.9 – 45.2%). We did not test multivariate mediations or interactions.

IGF-BP2 Mediates Age’s Direct Association with Future Dementia Severity, as measured by dEQ. APOE = apolipoprotein e4 status; CFI = Comparative Fit Index; GDS = Geriatric Depression Scale; HCY = serum homocysteine; HgbA1c = serum hemoglobin A1c; IGF-BP2 = Insulin-like Growth Factor Binding Protein 2; RMSEA = Root Mean Square Error of Association.*All observed variables except AGE are adjusted for APOE, education, ethnicity, gender, GDS, HCY, and HgbA1c (paths not shown for clarity). The covariates are densely intercorrelated.

Figure 1. IGF-BP2 Mediates Age’s Direct Association with Future Dementia Severity, as measured by dEQ. APOE = apolipoprotein e4 status; CFI = Comparative Fit Index; GDS = Geriatric Depression Scale; HCY = serum homocysteine; HgbA1c = serum hemoglobin A1c; IGF-BP2 = Insulin-like Growth Factor Binding Protein 2; RMSEA = Root Mean Square Error of Association.
*All observed variables except AGE are adjusted for APOE, education, ethnicity, gender, GDS, HCY, and HgbA1c (paths not shown for clarity). The covariates are densely intercorrelated.

Table 2. Potential Mediators of Age’s-Specific Dementing Effect.

Adiponectin (APN)
Angiopoetin-2N (ANG-2)
Compliment 3 (C3)
Creatinine Kinase-MB (CK-MB)
Epidermal Growth Factor Receptor 1 (EGFR)
FAS
Follicle stimulating hormone (FSH)
Glutathione S-Transferase (GST)
granulocyte colony stimulating factor (G-CSF)
Insulin-like Growth Factor-1 (IGF-I)
Insulin-like Growth Factor-Binding Protein 2 (IGF-BP2)*
Interleukin 5 (IL-5)
Myoglobin (MyG)
Pancreatic Polypeptide (PP)
Plasminogen Activator Inhibitor type 1(PAI-1)
Platelet-Derived Growth Factor (PDGF)
Progesterone
Resistin
S100b
Serum Amyloid P (SAP)
Thyroxine Binding Globulin (TBG)
von Willebrand Factor (vWF)
*Previously recognized aging biomarker [14]. Previously recognized biomarker of δ [13].

Table 3. Mediation Effects (Class 1).

Mediating BiomarkersAdjusted Path a
(Figure 1)
Z (p)Effect (%)
Adiponectin (APN)-0.22, p < 0.001-4.26 (<0.001)11.7
Angiopoetin-2 (ANG-2)-0.20, p < 0.001-4.72 (<0.001)13.8
Compliment 3 (C3)-0.27, p < 0.0014.46 (<0.001)11.9
Creatinine Kinase-MB (CK-MB)-0.21, p < 0.001-5.50 (<0.001)13.1
Epidermal Growth Factor Receptor 1 (EGFR)-0.19, p < 0.001-6.15 (<0.001)22.5
Fatty Acid Synthase (FAS)-0.22, p < 0.001-2.65 (0.004)14.1
Follicle stimulating hormone (FSH)-0.20, p < 0.001-5.27 (<0.001)14.0
Glutathione S-Transferase (GST)-0.28, p < 0.0013.80 (<0.001)17.2
granulocyte colony stimulating factor (G-CSF)-0.21, p < 0.001-5.22 (<0.001)12.6
Insulin-like Growth Factor-1 (IGF-I)-0.25, p < 0.0011.83 (0.03)9.9
Insulin-like Growth Factor-Binding Protein 2 (IGF-BP2)-0.13, p < 0.001-8.85 (<0.001)45.2
Interleukin 5 (IL-5)-0.27, p < 0.0012.57 (0.005)17.6
Myoglobin (MyG)-0.30, p < 0.0014.87 (<0.001)21.6
Pancreatic Polypeptide (PP)-0.21, p < 0.001-4.71 (<0.001)14.1
Plasminogen Activator Inhibitor type 1(PAI-1)-0.21, p < 0.001-5.50 (<0.001)15.0
Platelet-Derived Growth Factor (PDGF)-0.29, p < 0.0014.22 (<0.001)17.8
Progesterone-0.27, p < 0.0012.56 (0.005)12.2
Resistin-0.20, p < 0.001-4.39 (<0.001)13.8
S100b-0.28, p < 0.0014.97 (<0.001)18.1
Serum Amyloid P (SAP)-0.21, p < 0.001-5.66 (<0.001)13.5
Thyroxine Binding Globulin (TBG)-0.22, p < 0.001-4.60 (<0.001)9.9
von Willebrand Factor (vWF)-0.22, p < 0.001-3.92 (<0.001)11.0

Table 4 presents other age-related proteins, unrelated to δ by path b (Figure 1). Alpha2-macroglobulin (α2M), Interferon-gamma (IFN-γ), Interleukin 10 (IL-10), Interleukin 12-p40 (IL-12p40), Interleukin 15 (IL-15), Prolactin (PRL), Stem Cell Factor (SCF), Thrombopoietin (THPO), and Tumor Necrosis Factor alpha (TNF-α) had previously been associated with δ in non-Hispanic White (NHW) TARCC participants [1,13]. None were associated with δ in these ethnicity adjusted models.

Table 4. Other Age-Related Biomarkers (unrelated to the dEQ by Path b).

Adrenocorticotropic Hormone (ACTH)
alpha1-antitrypsin (A1AT)
alpha2-macroglobulin (α2M)*
alpha-Fetoprotein (α-FP)
Amphiregulin (AREG)
Angiotensinogen
AXL
Betacellulin
Bone Morphogenic Protein 6 (BMP6)
Cortisol
Eotaxin-3
Epiregulin (EREG)
FAS-Ligand (FAS-L)
Heparin-binding EGF-like growth factor (HB-EGF)
Hepatocyte Growth Factor (HGF)
Interferon-gamma (IFN-γ)*
Interleukin 1 receptor (IL-1r)
Interleukin 3 (IL-3)
Interleukin 7 (IL-7)
Interleukin 10 (IL-10)*
Interleukin 12-p40 (IL-12p40)*
Interleukin 13 (IL-13)
Interleukin 15 (IL-15)††
Interleukin 16 (IL-16)
Lipoprotein a
Luteinizing Hormone (LH)
Matrix Metalloproteinase type 3 (MMP-3)
Prolactin (PRL)*
Prostatic Acid Phosphatase (PAP)
Pulmonary and Activation-Regulated Chemokine (PARC)
Serum Glutamic Oxaloacetic Transaminase (SGOT)
Stem Cell Factor (SCF) *
Thrombopoietin (THPO)*
Thrombospondin-1 (THBS1)
Thymus-Expressed Chemokine (TECK)
Tissue Factor (TF)
Tissue Growth Factor alpha (TGF-α)
Tissue Inhibitor of Metalloproteinase type 1 (TIMP-1)
Tumor Necrosis Factor alpha (TNF-α)*
Tumor Necrosis Factor-Related Apoptosis-Inducing Ligand Receptor 3 (TRAIL-R3)
Vascular Endothelial Growth Factor (VEGF)‡
*Previously recognized δ biomarkers in non-Hispanic Whites (NHW) only [13], [1] (i.e., unconfirmed as a biomarker of dEQ in this ethnicity adjusted analysis). ††Previously recognized ethnicity adjusted δ biomarker [11]. ‡δ-related trend, p = 0.002.

Table 5 presents δ-related proteins unrelated to age by path c (Figure 1). Beta2-Microglobulin (β2M) was the only previously recognized δ biomarker [13]. The remainders are newly recognized as such. Table 6 lists biomarkers that were related neither to age nor to δ.

Table 5. Age-independent dEQ Biomarkers (unrelated to Age by Path c).

beta2-Microglobulin (β2M)*
Brain-Derived Neurotrophic Factor (BDNF)
Carcinoembryonic antigen (CEA)
CD40
Chromogranin A
Fatty Acid Binding Protein (FABP)
Growth Hormone
Immunoglobulin M (IgM)
Insulin
Interleukin 8 (IL-8)
Interleukin 18 (IL-18)
Macrophage Inflammatory Protein type 1 alpha (MIP-1α)
RANTES
Sex Hormone Binding Globulin (SHBG)
Tenascin C
Testosterone
Tumor Necrosis Factor receptor type II (TNF-RII)
Vascular Cell Adhesion Molecule type 1 (VCAM-1)
*Previously recognized δ biomarker [13].

Table 6. Unrelated Biomarkers.

Angiotensin Converting Enzyme (ACE)
Agouti-Related Protein (AgRP)
Apolipoprotein A1 (APOA1)
Apolipoprotein CIII (APOCIII)
Apolipoprotein H (apoH)
B Lymphocyte Chemoattractant (BLC)
Cancer Antigen 125 (CA 125)
Cancer Antigen 19-9 (CA 19-9)
CD40 Ligand
Connective Tissue Growth Factor (CTGF)
C Reactive Protein (CRP)
ENA-78 (ENA-78)
EN-RAGE (EN-RAGE)
Epidermal Growth Factor (EGF)
Eotaxin
Factor VII
Ferritin
Fibrinogen
GRO alpha (GROa)
Haptoglobin
Human CC Cytokine (HCC-4)
I-309
Immunoglobulin A (IgA)
Immunoglobulin E (IgE)
Intercellular Adhesion Molecule, type 1 (ICAM-1)
Interleukin 1 receptor antagonist (IL-1ra)
Leptin
Macrophage Inflammatory Protein type 1 beta (MIP-1b)
Macrophage Derived Chemokine (MDC)
Macrophage Migration Inhibitory Factor (MMIF)
Monocyte Chemotactic Protein type 1 (MCP-1)
Prostate Specific Antigen (PSA)
Serum Amyloid P (SAP)
Soluble Advanced Glycosylation End Product-Specific Receptor) (sRAGE)
Sortilin
Thyroid Stimulating Hormone (TSH)
Tumor Necrosis Factor beta (TNFb)
Vitamin D Binding Protein (VDBP)
Previously recognized δ biomarker (Bishnoi, Palmer & Royall, 2015).

Discussion

We have surveyed more than 100 potential mediators of age’s specific and significant association with the latent dementia phenotype, δ. Our sample size was large, and we were powered to detect even statistically weak effects. All our findings have been replicated in random subsets of TARCC’s data. We have replicated our previously reported association between age and IGF-BP2, and three of our previously observed age-independent associations with δ, even though 1) TARCC’s sample size has increased over time, 2) we are using a new δ homolog, 3) the biomarkers are being used here to predict future cognitive performance, and 4) the prior associations were obtained using raw biomarker data, while these employ normalized variables.

We have identified four classes of proteins: 1) potential mediators of age’s significant direct effect on δ, 2) δ-independent age-related proteins, 3) age-independent predictors of δ, and 4) proteins related neither to age nor to δ.

While many proteins were related to age, only a subset was also associated with δ (Class 1, Table 2). δ in turn has been associated with atrophy in the Default Mode Network (DMN) [15]. This suggests that the mediators in Table 2 may effect aging-specific changes to the structure or function of the DMN.

The DMN is a network of interconnected brain regions that are particularly active in the resting state [16]. Functional connectivity studies in older subjects have shown decreased DMN connectivity [1720] and less deactivation during task performance [21,22]. The impact of aging-related serum biomarkers on the integrity and functioning of the DMN is not well-studied. However, Thompson et al. [23] found that elevated Serum protein S100B levels significantly correlated with DMN activity. S100B has been confirmed by this analysis to mediate age’s specific-effect on a DMN-related cognitive construct (i.e., δ).

Our observations may help further clarify age’s specific effects on cognitive function. First, although age has both direct and indirect effects on observed cognitive performance [5], only its indirect effects, mediated by δ, are functionally salient, and thus “dementing.” This constrains “senility” and its biology to an effect on intelligence.

Second, δ has been shown to be “agnostic” to dementia’s etiology [6]. Age’s association with δ suggests that it too may have a role in determining all cause dementia risk, not just AD risk. This risk may not be conveyed through neurodegeneration. Age’s specific association with δ is characterized by lesser levels of AD-specific lesions [3].

Age accounts for only 5% of δ’s variance in this sample (per the base model). Regardless, correcting any δ-related pathology might improve dementia status, including age’s small effect. The mediators identified in Table 2 then, may offer targets for the remediation of age’s specific contribution.

Each Class 1 protein is a partial mediator of age’s contribution, ranging from Thyroxine Binding Globulin (TBG) (9.9%) to IGF-BP2 (45.2%)(Table 3). In their aggregate, they may have interacting effects. For example, S100b is elevated after cardiac surgery and correlated with post-operative cognitive impairments [24]. It binds to the receptor for advanced glycation end products (RAGE), which induces nuclear factor kappa-B (NF-kappaB)-regulated cytokines, including Compliment 3 (C3) [25]. However, we did not test multivariate interactions.

We note that not all of the mediators attenuate age’s direct effect. C3, Glutathione S-Transferase (GST), Interleukin 5 (IL-5), Myoglobin (MyG), Platelet-Derived Growth Factor (PDGF), Progesterone, and S100b accentuated age’s adverse effect on δ.

We had previously identified IGF-BP2 as a strong predictor of a δ ortholog targeting age itself instead of IADL [14]. That ortholog was significantly associated with δ, suggesting overlap between δ and Aging Proper. It has since been shown that δ mediates the majority of age’s effect on cognition, but not all [5]. IGF-BP2’s appearance in Class 1 confirms its contribution to age’s dementing aspect (i.e., “Senility”).

Class 1 also contains Insulin-like Growth Factor-1 (IGF-1). The appearance of both IGF-1 and IGF-BP2 among the Class 1 mediators strongly implicates the insulin-like growth factor (IGF) system in Aging Proper. The IGF system is comprised of two growth factors (IGF-I and 2), six high affinity binding proteins (IGF-BP 1 to 6) and four receptors [26,27]. Most of these are not available in TARCC.

Serum levels of IGFs I and II appear to mediate growth hormone (GH)-related somatotrophic changes in humans. These “somatomedins” circulate in non-covalent associations with IGF-BP2. It has been suggested that decreased function of the GH-somatomedin axis is responsible for age-specific anabolic changes (e.g., the “somatopause”) [28]. Interestingly, GH itself appears to be an age-independent δ-related protein (Class 3) (Table 5).

Serum IGF-BP2 increases with age, and high serum levels have been associated with greater disability, poorer physical performance, reduced muscle strength and lower mineral bone density [29]. Serum IGF-I declines with age [30,31]. Consistent with those findings, IGF-BP2’s association with age was positive (Figure 1) while IGF-I’s association with age was inverse. Insulin itself is related to δ, but not to age (in this HgbA1c adjusted analysis) (Table 5).

MyG and Creatinine Kinase-MB (CK-MB)’s appearance among the Class 1 mediators, recent associations between simple motor tasks and dementia risk [32], and the age-related somatomedins in Classes 1 and 2 lend credence to the hypothesis that there is a cognitive ortholog of somatic “frailty” [33].

IGF-I had an inverse adverse effect on dEQ. Other adverse mediators were Adiponectin (APN), Angiopoetin-2N (ANG-2), C3, CK-MB, Fatty Acid Synthase (FAS), Follicle stimulating hormone (FSH), Glutathione S-Transferase (GST), IGF-BP2, PP, PDGF, Progesterone, Resistin, and vWF. Like IGF-1, C3, GST, PDGF, and Progesterone declined with age [while also accentuating age’s effect on dEQ (see above)]. The others increased significantly with age. Thus all might contribute to age’s adverse effect on δ.

Epidermal Growth Factor Receptor 1 (EGFR), granulocyte colony stimulating factor (G-CSF), IL-5, MyG, Plasminogen Activator Inhibitor type 1(PAI-1), S100b, Serum Amyloid P (SAP) and Thyroxine Binding Globulin (TBG) had positive associations with δ and might offer some protection from age’s otherwise adverse effects.

The mechanism(s) by which the other Class 1 proteins affect δ remain to be elucidated. However, aging’s pathophysiology will be necessarily constrained, by δ’s mediation of its dementing effects, to the physiological processes that mediate intelligence. Two candidate processes might be synaptogenesis and network connectivity. C3, IGF-I, Progresterone, PAI-1, and S100b, all Class 1 mediators, are modulators of synaptic structure and function [3438].

Class 2 (Table 4) comprises proteins that although age-related, never the less fail to contribute to dementia via δ scores. They may mediate non-dementing age-related cognitive changes via g’s “domain-specific” residuals (e.g., memory, etc.). Alternatively, they may contribute to Aging Proper’s manifestation in other tissues or organs.

Notable among these are multiple EGFR agonist ligands, including Amphiregulin (AREG), Betacellulin, Epiregulin (EREG), Heparin-binding EGF-like growth factor (HB-EGF), and Tissue Growth Factor alpha (TGF-α) [39]. Epidermal Growth Factor (EGF), another EGFR agonist ligand, showed statistically insignificant trends as a potential mediator. The EGFR itself is a Class 1 Mediator (Table 2). These findings suggest the EGFR family of agonist ligands may have potential roles as therapeutic agents for age-specific cognitive and /or somatic decline.

On the other hand, several EGFR antagonists are approved by the Federal Drug Administration (FDA) for the treatment of certain cancers. These might be expected to have adverse effects, according to our findings. Chemotherapy has been noted to adversely impact connectivity in the DMN [40]. Such effects might explain reports of disability due to “chemobrain” in the literature [41]. They also illustrate the potential for reciprocal relationships between cognitive performance and cancer risks. It has been suggested both that chemotherapy is a risk factor for cognitive decline in late life [42], and that AD cases are relatively protected from cancer [43].

Class 2 also includes almost all of TARCC’s interleukin panel. The interleukins’ appearance in Table 4 suggests that inflammatory mechanisms may mediate age-specific changes outside the brain (and /or non-dementing aspects of cognition). IL-6 has been reported to protect cognition in centenarians [44], but is not in TARCC’s biomarker panel. The interleukin 1 receptor antagonist is related neither to age nor to δ (Table 6).

Eight of the eleven proteins we previously associated with δ in TARCC [13] are also in Class 2 (Table 5), including IFN-γ, Interleukins 10, 15, 12p40, and the Interleukin 1 receptor (IL-1r). Five of those eight [i.e., alpha2-macroglobulin (α2M), IFN-γ, IL-10, IL-12-p40 (IL-12p40), and SCF], exhibited otherwise significant trends in their associations with δ, which could not survive Bonferroni correction. Their previously reported associations were specific to NHW, while the current models were ethnicity adjusted. It remains to be seen whether ethnicity-specific effects on δ can be confirmed for any of the Class 2 proteins in Table 4.

Table 5 identifies many newly recognized age-independent determinants of dEQ (Class 3). Their relationships with δ are beyond the scope of this manuscript. However, GH’s appearance on this list is of interest given the prominence of other somatomedins among the Class 1 and 2 proteins (Tables 2 and 4).

In summary, we have surveyed over 100 serum proteins for their possible roles as mediators of age’s specific association with a latent dementia phenotype. 22 potential mediators were identified. These may offer targets for the disabling aspects of Aging Proper. An additional 41 age-related proteins were identified. These may mediate age’s effects on other organs. Notable among them are the EGFR and many of its ligands. Some EGFR ligands may protect the brain and other organs from age-related changes. However, this may occur at a risk of incurring cancer. Conversely, the use of EGFR antagonists in cancer treatment may accelerate the effects of Aging Proper in the brain and other organs.

Methods

Subjects

Subjects included n = 3385 TARCC participants, including 1240 cases of AD, 688 MCI cases, and 1384 NC. Each underwent serial annual standardized clinical examinations, culminating in a consensus clinical diagnosis of NC, MCI or AD. Institutional Review Board approval was obtained at each site and written informed consent was obtained from all participants.

δ’s Indicators included Logical Memory II (LMII) [45], Visual Reproduction I (VRI) [45], the Controlled Oral Word Association (COWA) [46], Digit Span Test (DST) [45] and Instrumental Activities of Daily Living (IADL) [47]. All tests were available in Spanish translation. The indicators were not adjusted for this analysis. The resulting unadjusted δ homolog was validated by its association with dementia severity, as measured by the Clinical Dementia Rating Scale sum of boxes (CDR) [48] and by ROC analysis.

TARCC’s methodology has been described elsewhere [49]. Serum samples were sent frozen to Rules-Based Medicine (RBM) in Austin, TX. There they were assayed without additional freeze-thaw cycles. RBM conducted multiplexed immunoassay via their human multi-analyte profile (human MAP). A complete listing of the biomarker panel we employed is available at http://www.rulesbasedmedicine.com.

Raw biomarker data were inspected to ascertain their normality. Data points beyond 3.0 standard deviations (SD) about the mean were labeled as “outliers” and deleted. Logarithmic transformation was used to normalize highly skewed distributions. The data were then standardized to a mean of zero and unit variance.

Covariates

All observed measures in the structural models were adjusted for APOE ε4 burden, education, ethnicity, gender, HCY, and HgbA1c. Measurements of HCY, HgbA1c and APOE ε4 genotyping were performed in the Ballantyne laboratory at the Baylor College of Medicine. HgbA1c was measured in whole blood by the turbidimetric inhibition immunoassay (TINIA). HCY was measured in serum using the recombinant enzymatic cycling assay (i.e., Roche Hitachi 911).

APOE genotyping

APOE genotyping was conducted using standard polymerase chain reaction (PCR) methods [50]. APOEε4 status was coded dichotmously based on the presence of one or more ε4 alleles. TARCC’s RBM biomarkers exhibit significant batch effects. Therefore, each biomarker was additionally adjusted for dichotomous dummy variables coding batch.

Statistical analysis

Analysis Sequence

This analysis was performed using Analysis of Moment Structures (AMOS) software [51]. The maximum likelihood estimator was chosen. All observed indicators were adjusted for age, education, ethnicity and gender. Co-variances between the residuals were estimated if they were significant and improved fit.

We used the ethnicity equivalent δ homolog (“dEQ”) as previously described [1]. That homolog has been reported to 1) have excellent fit (i.e., χ2/df = 181/24, p < 0.001; CFI = 0.97; RMSEA = 0.05), 2) have acceptable factor determinacy by Grice’s Method [52], 3) exhibit factor equivalence across ethnicity, 4) to be strongly correlated with dementia severity as measured by the CDR (r = 0.99, p <0.001) and 5) to exhibit an AUC of 0.97 (CI: 0.97-0.98) for the discrimination between AD cases and controls (in Visit 2 TARCC data). For the purposes of this analysis, dEQ was again constructed in Visit 2 data, but without any covariates, specifically age, ethnicity, GDS, gender, HCY, HGbA1c and APOE ε4 status.

dEQ and g’ factor weights were applied to Visit 2 observed data to generate Visit 2 dEQ and g’ composite scores (i.e., dEQ v2 and g’ v2, respectively). g’ is dEQ’s residual in Spearman’s g [53]. The composite scores were used as observed dependent variables in multiple regression models of age’s direct association with covariate adjusted Visit 2 g’ and dEQ.

Next, we constructed a longitudinal mediation model in SEM (Figure 1). Such models can arguably be interpreted causally [54]. Path “a” represents age’s direct association with Visit 2 dEQ scores. Path “b” represents the Visit 1 biomarker’s independent effect on dEQ. Bonferroni correction to p <0.001 was used to offset the potential for Type 2 error after multiple comparisons. When both paths were significant, we considered path “c”. The biomarker’s mediation effect on age’s direct association can then be calculated by MaKinnon’s method [55].

The mediation models were constructed in a randomly selected subset of TARCC participants, comprising approximately 50% of the subjects (i.e., Class 1: n = 1691). As a test of each model’s generalizability to the remainder (n = 1694), each mediation path’s significant direct association was constrained across the two groups, and model fit compared across constrained and unconstrained conditions [56,57]. Mediation effects were calculated in the constrained models.

Missing data

We used the newest instance of TARCC’s dataset (circa 2016). The entire dataset was employed. Clinical diagnoses were available on 3385 subjects, 2861 of whom had complete data for δ’s cognitive indicators and covariates. Modern Missing Data Methods were automatically applied by the AMOS software [58]. AMOS employs Full information Maximum Likelihood (FIML) [59,60]. Only the ROC analyses, performed in Statistical Package for the Social Sciences (SPSS) [61], were limited to complete cases.

Fit indices

Fit was assessed using four common test statistics: chi-square, the ratio of the chi-square to the degrees of freedom in the model (CMIN /DF), the comparative fit index (CFI), and the root mean square error of approximation (RMSEA). A non-significant chi-square signifies that the data are consistent with the model [62]. However, in large samples, this metric conflicts with other fit indices (insensitive to sample size) show that the model fits the data very well. A CMIN/DF ratio < 5.0 suggests an adequate fit to the data [63]. The CFI statistic compares the specified model with a null model [64]. CFI values range from 0 to 1.0. Values below 0.90 suggest model misspecification. Values approaching 1.0 indicate adequate to excellent fit. An RMSEA of 0.05 or less indicates a close fit to the data, with models below 0.05 considered “good” fit, and up to 0.08 as “acceptable“ [65]. All fit statistics should be simultaneously considered when assessing the adequacy of the models to the data.

Acknowledgements

Investigators from the Texas Alzheimer’s Research and Care Consortium: Baylor College of Medicine: Rachelle Doody MD, PhD, Mimi M. Dang MD, Valory Pavlik PhD, Wen Chan PhD, Paul Massman PhD, Eveleen Darby, Monica Rodriguear RN, Aisha Khaleeq; Texas Tech University Health Sciences Center: Chuang-Kuo Wu MD, PhD, Matthew Lambert PhD, Victoria Perez, Michelle Hernandez; University of North Texas Health Science Center: Thomas Fairchild PhD, Janice Knebl DO, Sid E. O’Bryant PhD, James R. Hall PhD, Leigh Johnson PhD, Robert C. Barber PhD, Douglas Mains, Lisa Alvarez, Rosemary McCallum; University of Texas Southwestern Medical Center: Perrie Adams PhD, Munro Cullum PhD, Roger Rosenberg MD, Benjamin Williams MD, PhD, Mary Quiceno MD, Joan Reisch PhD, Ryan Huebinger PhD, Natalie Martinez, Janet Smith; University of Texas Health Science Center – San Antonio: Donald Royall MD, Raymond Palmer PhD, Marsha Polk; Texas A&M University Health Science Center: Farida Sohrabji PhD, Steve Balsis PhD, Rajesh Miranda, PhD; Essentia Institute of Rural Health: Stephen C. Waring DVM, PhD; University of North Carolina: Kirk C. Wilhelmsen MD, PhD, Jeffrey L. Tilson PhD, Scott Chasse PhD.

Conflicts of Interest

DRR and RFP have disclosed the results of these analyses to the University of Texas Health Science Center at San Antonio (UTHSCSA), which has filed patent application 2012.039.US1.HSCS and provisional patents 61/603,226 and 61/671,858 relating to the latent variable δ’s construction and biomarkers.

Funding

This study was made possible by the Texas Alzheimer’s Research and Care Consortium (TARCC) funded by the state of Texas through the Texas Council on Alzheimer’s Disease and Related Disorders.

References

  • 1. Royall DR, Palmer RF. Thrombopoietin is associated with δ’s intercept, and only in Non-Hispanic Whites. Alzheimers Dement (Amst). 2016; 3:35–42. https://doi.org/10.1016/j.dadm.2016.02.003 [PubMed]
  • 2. Royall DR, Palmer RF. Alzheimer’s disease pathology does not mediate the association between depressive symptoms and subsequent cognitive decline. Alzheimers Dement. 2013; 9:318–25. https://doi.org/10.1016/j.jalz.2011.11.009 [PubMed]
  • 3. Gavett BE, John SE, Gurnani AS, Bussell CA, Saurman JL. The role of Alzheimer’s and cerebrovascular pathology in mediating the effects of age, race, and apolipoprotein E genotype on dementia severity in pathologically confirmed Alzheimer’s disease. J Alzheimers Dis. 2016; 49:531–45. https://doi.org/10.3233/JAD-150252 [PubMed]
  • 4. Nelson PT, Head E, Schmitt FA, Davis PR, Neltner JH, Jicha GA, Abner EL, Smith CD, Van Eldik LJ, Kryscio RJ, Scheff SW. Alzheimer’s disease is not “brain aging”: neuropathological, genetic, and epidemiological human studies. Acta Neuropathol. 2011; 121:571–87. https://doi.org/10.1007/s00401-011-0826-y [PubMed]
  • 5. Royall DR, Palmer RF. Aging is a weak but relentless determinant of dementia severity. Oncotarget. 2016; 7:13307–18. [PubMed]
  • 6. Gavett BE, Vudy V, Jeffrey M, John SE, Gurnani AS, Adams JW. The δ latent dementia phenotype in the uniform data set: cross-validation and extension. Neuropsychology. 2015; 29:344–52. https://doi.org/10.1037/neu0000128 [PubMed]
  • 7. Palmer RF, Royall DR. Future dementia status is almost entirely explained by the latent variable δ’s intercept and slope. J Alzheimers Dis. 2016; 49:521–29. https://doi.org/10.3233/JAD-150254 [PubMed]
  • 8. Royall DR, Palmer RF. δ scores predict MCI and AD conversions from non-demented states. J Prev Alzheimers Dis. 2015a; 2:337–38. Abstract.
  • 9. von Guten A, Ebbing K, Imhof A, Giannakopoulos P, Kövari E. Brain aging in the oldest-old. Curr Gerontol Geriatr Res. 2010. pii: 358531. doi: 10.1155/2010/358531.
  • 10. Crary JF, Trojanowski JQ, Schneider JA, Abisambra JF, Abner EL, Alafuzoff I, Arnold SE, Attems J, Beach TG, Bigio EH, Cairns NJ, Dickson DW, Gearing M, et al. Primary age-related tauopathy (PART): a common pathology associated with human aging. Acta Neuropathol. 2014; 128:755–66. https://doi.org/10.1007/s00401-014-1349-0 [PubMed]
  • 11. Bishnoi R, Palmer RF, Royall DR. Serum interleukin (IL)-15 as a biomarker of Alzheimer’s disease. PLoS One. 2015; 10:e0117282. https://doi.org/10.1371/journal.pone.0117282 [PubMed]
  • 12. Bishnoi RJ, Palmer RF, Royall DR. Vitamin D binding protein as a serum biomarker of Alzheimer’s disease. J Alzheimers Dis. 2015; 43:37–45. https://doi.org/10.3233/JAD-140042 [PubMed]
  • 13. Royall DR, Palmer RF. Ethnicity moderates dementia’s biomarkers. J Alzheimers Dis. 2015; 43:275–87. https://doi.org/10.3233/JAD-140264 [PubMed]
  • 14. Royall DR, Bishnoi RJ, Palmer RF. Serum IGF-BP2 strongly moderates age’s effect on cognition: a MIMIC analysis. Neurobiol Aging. 2015; 36:2232–40. https://doi.org/10.1016/j.neurobiolaging.2015.04.003 [PubMed]
  • 15. Royall DR, Palmer RF, Vidoni ED, Honea RA, Burns JM. The default mode network and related right hemisphere structures may be the key substrates of dementia. J Alzheimers Dis. 2012; 32:467–78. https://doi.org/10.3233/JAD-2012-120424 [PubMed]
  • 16. Raichle ME, MacLeod AM, Snyder AZ, Powers WJ, Gusnard DA, Shulman GL. A default mode of brain function. Proc Natl Acad Sci USA. 2001; 98:676–82. https://doi.org/10.1073/pnas.98.2.676 [PubMed]
  • 17. Andrews-Hanna JR, Snyder AZ, Vincent JL, Lustig C, Head D, Raichle ME, Buckner RL. Disruption of large-scale brain systems in advanced aging. Neuron. 2007; 56:924–35. https://doi.org/10.1016/j.neuron.2007.10.038 [PubMed]
  • 18. Damoiseaux JS, Beckmann CF, Arigita EJ, Barkhof F, Scheltens P, Stam CJ, Smith SM, Rombouts SA. Reduced resting-state brain activity in the “default network” in normal aging. Cereb Cortex. 2008; 18:1856–64. https://doi.org/10.1093/cercor/bhm207 [PubMed]
  • 19. Biswal BB, Mennes M, Zuo XN, Gohel S, Kelly C, Smith SM, Beckmann CF, Adelstein JS, Buckner RL, Colcombe S, Dogonowski AM, Ernst M, Fair D, et al. Toward discovery science of human brain function. Proc Natl Acad Sci USA. 2010; 107:4734–39. https://doi.org/10.1073/pnas.0911855107 [PubMed]
  • 20. Koch W, Teipel S, Mueller S, Buerger K, Bokde AL, Hampel H, Coates U, Reiser M, Meindl T. Effects of aging on default mode network activity in resting state fMRI: does the method of analysis matter? Neuroimage. 2010; 51:280–87. https://doi.org/10.1016/j.neuroimage.2009.12.008 [PubMed]
  • 21. Lustig C, Snyder AZ, Bhakta M, O’Brien KC, McAvoy M, Raichle ME, Morris JC, Buckner RL. Functional deactivations: change with age and dementia of the Alzheimer type. Proc Natl Acad Sci USA. 2003; 100:14504–09. https://doi.org/10.1073/pnas.2235925100 [PubMed]
  • 22. Grady CL, Springer MV, Hongwanishkul D, McIntosh AR, Winocur G. Age-related changes in brain activity across the adult lifespan. J Cogn Neurosci. 2006; 18:227–41. https://doi.org/10.1162/jocn.2006.18.2.227 [PubMed]
  • 23. Thompson WH, Thelin EP, Lilja A, Bellander BM, Fransson P. Functional resting-state fMRI connectivity correlates with serum levels of the S100B protein in the acute phase of traumatic brain injury. Neuroimage Clin. Epub ahead of print.
  • 24. Ali MS, Harmer M, Vaughan R. Serum S100 protein as a marker of cerebral damage during cardiac surgery. Br J Anaesth. 2000; 85:287–98. https://doi.org/10.1093/bja/85.2.287 [PubMed]
  • 25. Lian H, Yang L, Cole A, Sun L, Chiang AC, Fowler SW, Shim DJ, Rodriguez-Rivera J, Taglialatela G, Jankowsky JL, Lu HC, Zheng H. NFκB-activated astroglial release of complement C3 compromises neuronal morphology and function associated with Alzheimer’s disease. Neuron. 2015; 85:101–15. https://doi.org/10.1016/j.neuron.2014.11.018 [PubMed]
  • 26. Le Roith D. The insulin-like growth factor system. Exp Diabesity Res. 2003; 4:205–12. https://doi.org/10.1155/EDR.2003.205 [PubMed]
  • 27. Berryman DE, Christiansen JS, Johannsson G, Thorner MO, Kopchick JJ. Role of the GH/IGF-1 axis in lifespan and healthspan: lessons from animal models. Growth Horm IGF Res. 2008; 18:455–71. https://doi.org/10.1016/j.ghir.2008.05.005 [PubMed]
  • 28. Hoffman AR, Pyka G, Lieberman SA, Ceda GP, Marcus R. The somatopause. In: Muller EE, Cocchi D, Locatelli V, eds. Growth hormone and somatomedins during lifespan. Berlin: Springer Verlag, 1993.
  • 29. Hu D, Pawlikowska L, Kanaya A, Hsueh WC, Colbert L, Newman AB, Satterfield S, Rosen C, Cummings SR, Harris TB, Ziv E, and Health, Aging, and Body Composition Study. Serum insulin-like growth factor-1 binding proteins 1 and 2 and mortality in older adults: the Health, Aging, and Body Composition Study. J Am Geriatr Soc. 2009; 57:1213–18. https://doi.org/10.1111/j.1532-5415.2009.02318.x [PubMed]
  • 30. Aleman A, Verhaar HJ, De Haan EH, De Vries WR, Samson MM, Drent ML, Van der Veen EA, Koppeschaar HP. Insulin-like growth factor-I and cognitive function in healthy older men. J Clin Endocrinol Metab. 1999; 84:471–75. https://doi.org/10.1210/jcem.84.2.5455 [PubMed]
  • 31. Aleman A, Torres-Alemán I. Circulating insulin-like growth factor I and cognitive function: neuromodulation throughout the lifespan. Prog Neurobiol. 2009; 89:256–65. https://doi.org/10.1016/j.pneurobio.2009.07.008 [PubMed]
  • 32. Camargo EC, Weinstein G, Beiser AS, Tan ZS, DeCarli C, Kelly-Hayes M, Kase C, Murabito JM, Seshadri S. Association of physical function with clinical and subclinical brain disease: The Framingham Offspring Study. J Alzheimers Dis. 2016; 53:1597–608. https://doi.org/10.3233/JAD-160229 [PubMed]
  • 33. Panza F, Solfrizzi V, Barulli MR, Santamato A, Seripa D, Pilotto A, Logroscino G. Cognitive Frailty: A systematic review of epidemiological and neurobiological evidence of an age-related clinical condition. Rejuvenation Res. 2015; 18:389–412. https://doi.org/10.1089/rej.2014.1637 [PubMed]
  • 34. McAdory BS, Van Eldik LJ, Norden JJ. S100B, a neurotropic protein that modulates neuronal protein phosphorylation, is upregulated during lesion-induced collateral sprouting and reactive synaptogenesis. Brain Res. 1998; 813:211–17. https://doi.org/10.1016/S0006-8993(98)01014-2 [PubMed]
  • 35. Tournell CE, Bergstrom RA, Ferreira A. Progesterone-induced agrin expression in astrocytes modulates glia-neuron interactions leading to synapse formation. Neuroscience. 2006; 141:1327–38. https://doi.org/10.1016/j.neuroscience.2006.05.004 [PubMed]
  • 36. Deak F, Sonntag WE. Aging, synaptic dysfunction, and insulin-like growth factor (IGF)-1. J Gerontol A Biol Sci Med Sci. 2012; 67:611–25. https://doi.org/10.1093/gerona/gls118 [PubMed]
  • 37. Cho H, Joo Y, Kim S, Woo RS, Lee SH, Kim HS. Plasminogen activator inhibitor-1 promotes synaptogenesis and protects against aβ(1-42)-induced neurotoxicity in primary cultured hippocampal neurons. Int J Neurosci. 2013; 123:42–49. https://doi.org/10.3109/00207454.2012.724127 [PubMed]
  • 38. Mastellos DC. Complement emerges as a masterful regulator of CNS homeostasis, neural synaptic plasticity and cognitive function. Exp Neurol. 2014; 261:469–74. https://doi.org/10.1016/j.expneurol.2014.06.019 [PubMed]
  • 39. Ronan T, Macdonald-Obermann JL, Huelsmann L, Bessman NJ, Naegle KM, Pike LJ. Different Epidermal Growth Factor Receptor (EGFR) agonists produce unique signatures for the recruitment of downstream signaling proteins. J Biol Chem. 2016; 291:5528–40. https://doi.org/10.1074/jbc.M115.710087 [PubMed]
  • 40. Kesler SR. Default mode network as a potential biomarker of chemotherapy-related brain injury. Neurobiol Aging. 2014 (Suppl 2); 35:S11–19. https://doi.org/10.1016/j.neurobiolaging.2014.03.036 [PubMed]
  • 41. Wang X-M, Walitt B, Saligan L, Tiwari AF, Cheung CW, Zhang Z-J. Chemobrain: a critical review and causal hypothesis of link between cytokines and epigenetic reprogramming associated with chemotherapy. Cytokine. 2015; 72:86–96. https://doi.org/10.1016/j.cyto.2014.12.006 [PubMed]
  • 42. Small BJ, Scott SB, Jim HS, Jacobsen PB. Is cancer a risk factor of cognitive decline in late life? Gerontology. 2015; 61:561–66. https://doi.org/10.1159/000381022 [PubMed]
  • 43. Zhang Q, Guo S, Zhang X, Tang S, Shao W, Han X, Wang L, Du Y. Inverse relationship between cancer and Alzheimer’s disease: a systemic review meta-analysis. Neurol Sci. 2015; 36:1987–94. https://doi.org/10.1007/s10072-015-2282-2 [PubMed]
  • 44. Adriaensen W, Matheï C, Vaes B, van Pottelbergh G, Wallemacq P, Degryse JM. Interleukin-6 predicts short-term global functional decline in the oldest old: results from the BELFRAIL study. Age (Dordr). 2014; 36:9723. https://doi.org/10.1007/s11357-014-9723-3 [PubMed]
  • 45. Wechsler D. Wechsler Memory Scale – Third Edition. San Antonio, TX: The Psychological Corporation, 1997.
  • 46. Benton A, Hamsher K. Multilingual Aphasia Examination. AJA Associates, Iowa City, Iowa, 1989.
  • 47. Lawton MP, Brody EM. Assessment of older people: self-maintaining and instrumental activities of daily living. Gerontologist. 1969; 9:179–86. https://doi.org/10.1093/geront/9.3_Part_1.179 [PubMed]
  • 48. Hughes CP, Berg L, Danziger WL, Coben LA, Martin RL. A new clinical scale for the staging of dementia. Br J Psychiatry. 1982; 140:566–72. https://doi.org/10.1192/bjp.140.6.566 [PubMed]
  • 49. Waring S, O’Bryant SE, Reisch JS, Diaz-Arrastia R, Knebl J, Doody R; for the Texas Alzheimer’s Research Consortium. The Texas Alzheimer’s Research Consortium longitudinal research cohort: study design and baseline characteristics. Texas Public Health Journal. 2008; 60:9–13.
  • 50. Koch W, Ehrenhaft A, Griesser K, Pfeufer A, Müller J, Schömig A, Kastrati A. TaqMan systems for genotyping of disease-related polymorphisms present in the gene encoding apolipoprotein E. Clin Chem Lab Med. 2002; 40:1123–31. https://doi.org/10.1515/cclm.2002.197 [PubMed]
  • 51. Arbuckle JL. Analysis of Moment Structures-AMOS (Version 7.0) [Computer Program], SPSS, Chicago, 2006.
  • 52. Grice JW. Computing and evaluating factor scores. Psychol Methods. 2001; 6:430–50. https://doi.org/10.1037/1082-989X.6.4.430 [PubMed]
  • 53. Spearman C, Wynn Jones LL. Human Ability, Macmillan & Co., London, 1951.
  • 54. Kraemer HC, Stice E, Kazdin A, Offord D, Kupfer D. How do risk factors work together? Mediators, moderators, and independent, overlapping, and proxy risk factors. Am J Psychiatry. 2001; 158:848–56. https://doi.org/10.1176/appi.ajp.158.6.848 [PubMed]
  • 55. MacKinnon D. Analysis of mediating variables in prevention and intervention research. In: Czarees A, Beatty L (editors). Scientific methods for prevention intervention research. NIDA Research Monograph. 1994;139:137-153.
  • 56. Metz CE. Basic principles of ROC analysis. Semin Nucl Med. 1978; 8:283–98. https://doi.org/10.1016/S0001-2998(78)80014-2 [PubMed]
  • 57. Zweig MH, Campbell G. Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine. Clin Chem. 1993; 39:561–77. [PubMed]
  • 58. Palmer RF, Royall DR. Missing data? Plan on it!. J Am Geriatr Soc. 2010 (Suppl 2); 58:S343–48. https://doi.org/10.1111/j.1532-5415.2010.03053.x [PubMed]
  • 59. Schafer JL, Graham JW. Missing data: our view of the state of the art. Psychol Methods. 2002; 7:147–77. https://doi.org/10.1037/1082-989X.7.2.147 [PubMed]
  • 60. Graham JW. Missing data analysis: making it work in the real world. Annu Rev Psychol. 2009; 60:549–76. https://doi.org/10.1146/annurev.psych.58.110405.085530 [PubMed]
  • 61. Statistics PA. 18, Release Version 18.0.0, SPSS, Inc., 2009, Chicago, IL.
  • 62. Bollen KA, Long JS. Testing Structural Equation Models. Sage Publications, Thousand Oaks, CA, 1993.
  • 63. Wheaton B, Muthén B, Alwin DF, Summer GF. Assessing reliability and stability in panel models. In D.R.Heise (Ed.) Sociology Methodology San Francisco, CA: Jossey-Bass, 1977.
  • 64. Bentler PM. Comparative fit indexes in structural models. Psychol Bull. 1990; 107:238–46. https://doi.org/10.1037/0033-2909.107.2.238 [PubMed]
  • 65. Browne M, Cudeck R. Alternative ways of assessing model fit. In Testing structural equation models, Bollen KA, Long JS, eds. Sage Publications, Thousand Oaks, CA, pp. 136–162, 1993.
  • 66. Yesavage JA, Brink TL, Rose TL, Lum O, Huang V, Adey M, Leirer VO. Development and validation of a geriatric depression screening scale: A preliminary report. J Am Geriatr Soc. 1982; 29:164–71. https://doi.org/10.1111/j.1532-5415.1981.tb01759.x [PubMed]
  • 67. Folstein MF, Folstein SE, McHugh PR. “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res. 1975; 12:189–98. https://doi.org/10.1016/0022-3956(75)90026-6 [PubMed]