Plasma neurofilament light and its association with all-cause mortality risk among urban middle-aged men and women
BMC Medicine volume 20, Article number: 218 (2022)
Neurofilament light chain (NfL) is released into the blood during neuronal damage. NfL is linked to mortality in neurological disorders, remaining unexplored in population studies. We investigated whether initial (v1) and annualized change (δ) in plasma NfL can predict all-cause mortality in middle-aged dementia-free urban adults.
Longitudinal data were from 694 participants in the Healthy Aging in Neighborhoods of Diversity Across the Life Span study (HANDLS, mean agev1: 47.8 years, 42% male, 55.8% African American). Plasma NfL was measured prospectively at three visits. Analyses included Cox proportional hazards models for all-cause mortality risk and 4-way decomposition testing for interaction and mediation.
Unlike men, women exhibited a direct association between δNfL (above vs. below median) and all-cause mortality risk in both the minimally (HR = 3.91, 95% CI 1.10–13.9, p = 0.036) and fully adjusted models (HR = 4.92, 95% CI 1.26–19.2, p = 0.022), and for δNfL (per unit increase) in the full model (HR = 1.65, 95% CI 1.04–2.61, p = 0.034). In both models, and among women, 1 standard deviation of NfLv1 was associated with an increased all-cause mortality risk (reduced model: HR = 2.01, 95% CI 1.24–3.25, p = 0.005; full model: HR = 1.75, 95% CI 1.02–2.98, p = 0.041). Only few interactions were detected for cardio-metabolic risk factors. Notably, NfLv1 was shown to be a better prognostic indicator at normal hsCRP values among women, while HbA1c and δNfL interacted synergistically to determine mortality risk, overall.
These findings indicate that plasma NfL levels at baseline and over time can predict all-cause mortality in women and interacts with hsCRP and HbA1c to predict that risk.
Neurofilament light chain (NfL) is a cytoskeletal protein component exclusively expressed in neurons that is released into the extracellular fluids, including blood, during neuroaxonal damage. Emerging data indicates the clinical utility of blood-based measurements of NfL as a novel biomarker for neurodegenerative diseases. This methodological development for assaying plasma NfL has stimulated opportunities for large-scale applications in clinical practice and in randomized clinical trials as a method for identifying patients at risk for dementias, including Alzheimer’s disease (AD) . Thus far, NfL reflects sub-cortical large-caliber axonal degeneration [2, 3]. Plasma NfL levels correlate strongly with cerebrospinal fluid (CSF) NfL levels [1, 4], adding to its clinical utility in differential diagnoses for dementias [5,6,7,8] and other neurodegenerative diseases [9,10,11,12]. Therefore, plasma NfL measurements are advantageous given the invasiveness of CSF assessments and feasibility for long-term monitoring.
Based on a meta-analysis of > 60 studies, cognitive impairment, including overt dementia, has been linked to an increased risk of all-cause mortality . Despite this evidence, the relationship between mortality and blood assay predictors for dementia, including plasma NfL, remains largely unexplored, particularly among dementia-free middle-aged adults. A few studies have examined plasma NfL and found associations with mortality in patients with stroke [14, 15], sporadic Creutzfeldt-Jakob disease , and spontaneous subarachnoid and intracerebral hemorrhages [17, 18] and thus focused on prognostic outcomes for specific clinical diagnoses. More recently, attention was drawn as to the possible prognostic value of plasma NfL in the general population. For instance, a study of elderly adults in the Memory and Morbidity in Augsburg Elderly (MEMO) study found that serum NfL levels were associated with all-cause mortality as well as neuropsychological test and brain atrophy scores . Moreover, Kaeser et al. reported increasing plasma NfL levels with age in humans (n = 122; 21–107 years of age) and that plasma NfL correlated with plasma proteins involved in neural pathways. Importantly, they detected a positive association between plasma NfL levels and mortality among centenarians (n = 135) with a predictive value that exceeded that of cognitive and physical functioning measures, and this finding was replicated among nonagenarians in an independent cohort (n = 180) . These data indicate that plasma NfL may be a novel blood-based biomarker of mortality in various neurological diseases but may also be an indicator of neuronal damage in non-demented adults.
In general, men consistently experience a lower survival probability at each age compared to women . Thus, potential associations between plasma NfL and mortality should be tested both overall and separately among men and women. Moreover, the cardio-metabolic risk may also play a mediating and/or interactive role in this association, as they have been recently linked with elevated plasma NfL over time [22, 23, 24]. We therefore hypothesize that NfL and change in NfL over time interact with baseline cardio-metabolic conditions to affect mortality risk, while cardio-metabolic conditions’ effect on mortality risk may be mediated through NfL or change in NfL over time. A common framework examining both mediation and interaction can help test those associations simultaneously.
The present study examines the association between plasma NfL and all-cause mortality in a socio-economically diverse sample of middle-aged urban White and African American adults, overall and by sex. As a secondary objective, the study also tests the potential interactive and mediating effects of the body mass index (BMI), the allostatic load (AL) index, and other measures of cardio-metabolic risk on this relationship.
Participants in this study were chosen from the Healthy Neighborhoods of Diversity Across the Life Span study . HANDLS is an ongoing, prospective study of African American and White adults living in Baltimore, MD. Data from participants were collected from visit 1 (v1) from 2004 to 2009, visit 2 (v2) from 2009 to 2013, and visit 3 (v3) from 2013 to 2017. Each visit consisted of physical examinations, cognitive testing, and collection of fasting blood samples. Participants provided written informed consent.
Up to three repeats of plasma NfL concentrations were available from v1, v2, and/or v3. Among the 3720 initially recruited HANDLS participants, N = 731 had complete v1, v2, and/or v3 data of plasma NfL, of whom 694 had data on v1 plasma NfL (Additional file 1: Fig. S1). Given that all other covariates were either complete or imputed, the final analytic sample was N = 694 participants. The sample selection methodology for NfL is detailed in Additional file 1: Method S1 [26,26,27,29]. Compared to the unselected group (N = 3026 of 3720), the final eligible sample (N = 694) had a significantly higher proportion of White adults (44.2 vs. 40.1, p = 0.048) and of individuals with household incomes above poverty (70.9% vs. 56.0, p < 0.001), with no significant age or sex differences.
The HANDLS study is approved by the Institutional Review Board of the National Institutes of Health, National Institute of Environmental Health Sciences.
Plasma neurofilament light measurement
Fasting blood samples were collected between 9:30 am and 11:30 am into EDTA blood collection tubes. The tubes were centrifuged at 600g for 15 min followed by the removal of the buffy coat. The steps were repeated twice and visually inspected for hemolysis. Plasma samples were aliquoted and stored at − 80 °C. Plasma NfL levels were quantified using the Simoa® NF-light Advantage Kit on a Simoa® HD-X analyzer by Quanterix (Billerica, MA, USA) following the kit instructions. Samples from the different visits were run on the same plate for each individual, and plates were balanced for individuals within each demographic group (race/sex/poverty). Plasma samples were diluted fourfold, and concentrations were adjusted for this dilution correction. Pooled plasma samples from two individuals were run in duplicate on all plates, and the average intra-assay and inter-assay coefficients of variations were 4.5% and 7%, respectively. Information on the limits of detection and quantification for this assay has previously been described . Plasma NfL was measured with ≤ 3 repeats/participant, at v1, v2, and/or v3. The NfLv1 and δNfL exposures are detailed in Additional file 1: Method S2 . NfLv1 is the baseline NfL measured at v1 and was Loge-transformed to approximate a normal distribution. Previously, we reported the δNfL as the annualized rate of change between NfLv1 and NfLv2 measurements . Here, we included an additional measurement of NfL at v3. Therefore, in this study, the δNfL is the annualized rate of change between NfLv1, NfLv2, and NfLv3, on average, when these measurements were Loge-transformed.
Mortality status in the HANDLS cohort was obtained through linkages to the National Death Index (NDI), National Center for Health Statistics. Information about the underlying cause of death was obtained from death certificates and classified in accordance with the International Statistical Classification of Diseases, version 10 (ICD-10). Deaths attributed to CVD included CVD-related diagnosis codes (ICD-10 codes I00–99.9) listed as the underlying or contributing cause of death, respectively defined by the US Department of Health and Human Services and the Centers for Disease Control and Prevention as “the disease or injury that initiated the chain of morbid events that led directly and inevitably to death” and “all other significant diseases, conditions, or injuries that contributed to death but which did not result in the underlying cause of death” . Vital status information for all participants is available from enrollment (2004–2009) to December 31, 2018 (last date of death available). However, in order to exclude participants who did not survive between visits 1 and 3, all death events occurred after the completion of the v3 exam (2013–2018) (Additional file 1: Method S1).
BMI, AL, and co-morbidity
Visit 1 BMI was calculated as weight (kg) divided by height (m2). We relied on a previously reported method to compute the total AL score , with components also measured at visit 1. This method sums cardiovascular (systolic and diastolic blood pressure, pulse rate), metabolic (total cholesterol, high-density lipoprotein (HDL) cholesterol, glycosylated Hb (HbA1C), sex-specific waist-to-hip ratio), and inflammatory (serum albumin and high-sensitivity C-reactive protein (hsCRP)) risk indicators, as summarized in Additional file 1: Table S1 [33,33,34,35,36,37,38,39,41]. The total AL score (ALtotal) is equally weighted and may range from 0 to 9. The higher the ALtotal, the more the overall cardio-metabolic risk. Total cholesterol (mg/dl), HDL cholesterol (mg/dl), hsCRP (mg/dl), albumin (g/dl), and glycosylated hemoglobin (%) were quantified by contract laboratories (Quest Diagnostics, Chantilly, VA), using reference analytical methods. Using standard protocols, trained examiners measured waist-to-hip ratio, radial pulse (beats/min), and systolic and diastolic blood pressure (mmHg). Specifically, blood pressure was measured using a mercury sphygmomanometer, and the arithmetic mean of left and right systolic and diastolic pressures was used in this analysis (Additional file 1: Table S2).
Cardio-metabolic and co-morbidity were assessed using self-reported, and measured components were assessed at visit 1. These included hypertension (0 = no, 1 = yes), diabetes (0 = non-diabetic, 1 = pre-diabetic, 2 = diabetic), dyslipidemia (or statin use) (0 = no, 1 = yes), and self-reported history of any of several cardiovascular disease conditions (0 = no, 1 = yes). The latter component screened for the occurrence of several conditions, namely atrial fibrillation, angina, coronary artery disease, congestive heart failure, and myocardial infarction.
BMI, AL (total score and parameters from which components were derived), and cardio-metabolic co-morbidity binary indices were considered both a potential mediator and an effect modifier in the association between NfL exposures and mortality.
We assessed multiple other visit 1 covariates as potential confounders, given previous significant associations with plasma NfL and are considered antecedent risk factors to AL and cardio-metabolic risk. These included v1 age (continuous, years), sex (male, female), race (White, African American), poverty status (below vs. above 125% of the federal poverty line), and educational attainment (less than high school, high school, more than high school). We operationalized poverty status using the 2004 US Census Bureau poverty thresholds  based on household income and total family size (including children < 18 years). Few lifestyle and health-related factors were considered as potential confounders, namely current smoking status (0 = no vs. 1 = yes), illicit drug use (0 = no vs. 1 = yes, using any of marijuana, opiates, and cocaine), and the Healthy Eating Index 2010 (HEI-2010) , whereby overall diet quality was measured based on food and macronutrient-related dietary guidelines for Americans, total energy intake based on the average of two 24-h recalls (kcal/day), and the 20-item CES-D total score for depressive symptoms. Sex was the main effect modifier in our analyses.
We used Stata release 16  to conduct all analyses. Data aside from outcome and exposures was imputed, using chained equations (5 imputations, 10 iterations) [45, 46], with most covariates having < 10% missing data compared to the final eligible sample (i.e., N = 694). We characterized the overall analytic sample at baseline using means and proportions; Student’s t-tests were used to examine the sex differences in baseline characteristics. We used a series of bivariate and multivariable regression models to evaluate whether baseline characteristics varied by sex. To examine the association between plasma NfL exposures, we estimated a series of Cox proportional hazard regression models with sequential covariate adjustment. Age (years) at the visit was used as the underlying time scale, while v1 age was the time of entry (i.e., delayed entry). Sex-specific Kaplan-Meier survival curves were presented across binary NfL exposures (> vs. ≤ median) based on the distributions in the final selected sample with a cutoff of 1.966279 for NfLv1 and 0.0466214 for δNfL on the Loge-transformed scale, examining the time on study as the analytic time variable. In Cox proportional hazards models, heterogeneity by sex of the association between NfL exposures and mortality was tested by adding a two-interaction term (NfLv1 × sex; δNfL × sex) in separate models. These models were also stratified by sex. The general modeling strategy consisted of a basic model, adjusted for age, sex, race, and poverty status (model 1), to which other lifestyle and health-related covariates (listed in the “Covariates” section) were subsequently added (model 2).
BMI, AL (total score and components), and cardio-metabolic co-morbidity were separately assessed as mediating/interactive factors in the total effect of NfL exposures on all-cause mortality. Continuous potential mediators (e.g., BMI, AL total score) were transformed into a standardized z-score, while binary cardio-metabolic co-morbidity indices were coded as 0 = no and 1 = yes/any, for ease of interpretation. Diabetes was recoded as 0 = no and 1 = pre-diabetic or diabetic. Continuous AL components were also tested individually as potential mediators/effect modifiers. Specifically, the overall effect of each of our main exposures on all-cause mortality, in the presence of a mediator with which the exposure may interact, was decomposed into four distinctive components: (i) neither mediation nor interaction, (ii) interaction alone (and not mediation), (iii) both mediation and interaction, and (iv) only mediation (but not interaction). This four-way decomposition unifies the methods to attribute effects to interactions and methods that assess mediation, and this method has recently been introduced in Stata, allowing to estimate four-way decomposition using parametric or semi-parametric regression models. Importantly, the Med4way command  [https://github.com/anddis/med4way] was used to test the mediation and interaction of the total effects of NfL exposure on mortality with several mediators/effect modifiers, using Cox PH models for the outcome and linear or logistic regression models for each mediator/effect modifier (BMI, AL total score, continuous AL components, and binary cardio-metabolic co-morbidity indices). Four-way decomposition was applied to the total sample, and among men and women, separately, combining the findings from 5 imputations using Rubin’s rule . A logit link was specified for the mediating variable equation when mediators were binary. In a sensitivity analysis, NfL exposures (NfLv1 and δNfL) were entered as potential mediators while BMI, AL, AL continuous parameters, and co-morbidity indices were considered as the main exposure to assess the bidirectionality of the associations. In this analysis, all mediator equations were linear regression models. Total effects were interpreted as hazard ratios on the Loge scale based on Cox PH models, per SD of exposures if the exposure was continuous and for “exposed” vs. “unexposed” if the exposure was binary, which were then decomposed into four components. An effect size that would result in a hazard ratio > 1.5 was considered as moderate to strong.
In all models, sample selectivity due to missing exposure and outcome data, relative to the initially recruited sample, was adjusted for using a two-stage Heckman selection strategy . Thus, we first predicted an indicator of selection with socio-demographic factors, namely, v1 age, race, sex, and poverty status using a probit regression model, which yielded an inverse Mills ratio (IMR), a function of the probability of being selected given those socio-demographic factors. At a second stage, we estimated our Cox PH hazards regression models adjusted for the IMR in addition to the aforementioned covariates [30, 49]. A sensitivity analysis was also conducted to examine the sample selectivity by co-morbid conditions, while adjusting for socio-demographic factors, comparing the final analytic sample to those excluded from the original sample of N = 3720. We also compared the sex differences in co-morbidity in the final sample and the sample that was excluded.
Overall, 694 HANDLS participants were included in this study (Additional file 1: Fig. S1) and 42% were men and 55.8% were African American. The mean age at v1 was 47.8 years (30–65 years). After a mean follow-up of 11.2 years (range 3.86–14.31), 43 all-cause deaths occurred in the 694 participants (20 deaths among women, 23 deaths among men). The mean follow-up between v1 and v3 for δNfL was 7.77 years (range 4.9–12.5), when NfL was available for all 3 waves (N = 682). Overall, the incidence rate for all-cause death in the final eligible sample (N = 694) was 555 per 100,000 person-years (P-Y) (447 per 100,000 P-Y in women; 702 per 100,000 P-Y in men) and a median survival age of 78.1 years. Table 1 displays the sample characteristic distribution by sex, including the main exposures, covariates, potential mediators/effect modifiers, and the main outcome of interest. Most notably, men had significantly higher plasma NfL (Loge-transformed) compared to women at both visits 1 and 3 (p < 0.05). These differences remained comparable after further adjustment for age, race, and poverty status. Moreover, men were more likely than women to be pre-diabetic, with the reverse being observed in the case of diabetes. Self-reported history of CVD was more prevalent in women (15.9% vs. 8.7%, p = 0.007). While there was no sex difference detected in the mean ALtotal, several continuous components of AL were higher among women, namely hsCRP, total cholesterol, and HDL cholesterol, while the reverse was observed for albumin and DBP. Similar patterns were observed for binary AL components. Moreover, women had higher mean BMI compared to men, while having lower total caloric intakes and lower prevalence of current illicit drug use (p < 0.001). In a sensitivity analysis comparing those participants who were included in the analysis to those who were excluded, it was shown that inclusion was associated with a lower likelihood of diabetes, hypertension, and cardiovascular disease (p < 0.05). Moreover, among those excluded from the study, men were more likely than women to be pre-diabetic, as observed among those who were included. However, no sex differences were detected with respect to cardiovascular disease (p = 0.78), a finding that differed from the sex difference detected in the selected sample whereby men had a lower risk of CVD compared with women.
Figure 1 displays the findings from the Kaplan-Meier survival curve analyses across levels of δNfL (> vs. ≤ median: see Fig. 1 footnote for definitions), stratified by sex. The results show that only among women there was a marked difference in survivorship probability over time across those two levels, with log-rank test (df of 1) = 12.07, p = 0.0005, suggesting higher all-cause mortality risk in “> median” δNfL group vs. “≤ median.” Cox proportional hazards model and minimally and fully adjusted results are shown in Table 2, for both NfLv1 and δNfL, stratified by sex and expressed both as z-scores and binary exposures (above vs. below median). Similar to the findings from Fig. 1, binary δNfL among women was associated with increased risk for all-cause mortality in both the minimally adjusted model 1 (HR = 3.91, 95% CI 1.10–13.9, p = 0.036) and the fully adjusted model 2 (HR = 4.92, 95% CI 1.26–19.2, p = 0.022). When expressed as a standardized z-score, 1 SD increase in δNfL was associated with an increased risk of all-cause mortality among women, particularly in the fully adjusted model (HR = 1.65, 95% CI 1.04–2.61, p = 0.034). In both models 1 and 2, and among women, 1 SD of NfLv1 was associated with an increased risk of all-cause mortality (model 1: HR = 2.01, 95% CI 1.24–3.25, p = 0.005; model 2: HR = 1.75, 95% CI 1.02–2.98, p = 0.041), a relationship not detected for the binary NfLv1 exposure. Men did not exhibit a relationship between NfL exposure and mortality, and sex differences were detected for z-scored NfLv1 (model 1) and binary δNfL (both models).
Tables 3 and 4 show the findings from 4-way decomposition models (secondary analyses), which decompose the total effect (TE) of the specified NfL exposure on all-cause mortality into components attributed to mediation alone (pure indirect effect [PIE]), interaction alone (reference interaction [IR]), to both mediation and interaction (mediated Interaction [IM]), and neither mediation nor interaction (controlled direct effect [CDE]). Many of the total effects observed among women were either non-significant, and most significant associations (p < 0.05) were direct ones that were not explained by or interacted with cardio-metabolic risk factors. Nevertheless, there was some evidence of an antagonistic interaction between hsCRP and NfL, whereby the effect of NfL at v1 was stronger at lower levels of hsCRP among women and while the controlled direct effect was stronger than the total effect (HRcde,crp = 3.0, 95% CI 1.05–8.78, p = 0.04 vs. HRte,crp = 2.3, 95% CI 0.73–7.21, p = 0.11).
Moreover, we also examined the potential effects of each cardio-metabolic risk and co-morbidity factors in relation to all-cause mortality by sex, exploring the mediating and modifying effects with NfLv1 and δNfL as the main exposures (Additional file 1: Tables S1 and S2). Based on the findings presented in Additional file 1: Table S2, among men, when NfLv1 is set at its mean, BMI becomes potentially protective against all-cause mortality (β ± SE − 0.441 ± 0.196, p = 0.024). A potential protective effect of CVD history on mortality would be observed among women, if NfLv1 is set at this mean (β ± SE − 0.803 ± 0.249, p = 0.001). When δNfL was tested as the mediating and interactive factor for the association between cardio-metabolic risk and co-morbidity vs. all-cause mortality, several notable findings were observed. First, among men, when δNfL is set at its mean, BMI has an inverse association with all-cause mortality. The marginally significant inverse relationship between HbA1c and all-cause mortality observed when δNfL is set at its mean (CDE < 0, p < 0.10) becomes non-significant when δNfL is not fixed in the total population with this difference in effect caused by pure interaction as well as mediated interaction (IR > 0, IM > 0, p < 0.05), indicative of synergism between HbA1c and δNfL in determining mortality risk. Finally, among men, a CDE of − 0.684 (p = 0.037) for hypercholesterolemia was found when δNfL was set at its mean, indicating that at an average rate of change in NfL over time, hypercholesterolemia is potentially protective against all-cause mortality.
Few studies have examined the association of plasma NfL and mortality in population-based studies, and these studies have been conducted in elderly White adults . In this study, we examined the association between plasma NfL and all-cause mortality in a socio-economically diverse sample of middle-aged urban White and African American adults. Furthermore, we analyzed these associations both overall and by sex. In addition, it is also the first to test the potential interactive and mediating effects of BMI, AL index, and other measures of cardio-metabolic risk on this relationship. Here, we report that in women, δNfL was associated with an increased risk of all-cause mortality. In women, 1 SD of NfLv1 was associated with an increased risk of all-cause mortality, a relationship not detected for the binary NfLv1 exposure. Men did not exhibit a relationship between NfL exposure and mortality. We further explored these associations and found that in women, most of these associations were direct and were not explained by or did not interact with cardio-metabolic risk factors. We found one exception with a possible antagonistic interaction for hsCRP and NfLv1, indicating that NfLv1 is a better prognostic indicator at normal hsCRP values. Moreover, there was some evidence of synergistic interaction between HbA1c and δNfL in determining mortality risk, overall. These novel results add informative insight that plasma NfL can potentially be used as a biomarker to predict all-cause mortality in middle-aged women across different races.
Plasma NfL has recently gained significant attention due to its association with neurological diseases including sporadic and familial AD [5, 7, 50], frontotemporal degeneration , multiple sclerosis , traumatic brain injury , Parkinson’s disease , and other neurological disorders . In addition, levels of plasma NfL predict the future onset of dementia [8, 51]. In non-demented adults, plasma NfL may also have clinical utility as it has been shown to be associated with a faster decline in normalized mental status scores in White adults and in older adults . In older adults (median age 75 years) in the Multidomain Alzheimer’s Preventive Trial (MAPT), plasma NfL was associated with cognitive scores and executive function in adults with MCI but not in adults without cognitive impairment . Blood-based biomarkers can be powerful tools for assessing the risk and diagnosing of disease and have the advantage of being easier to obtain, less expensive, and not requiring highly specialized clinical staff as well as equipment that is typically needed for the acquisition of CSF and neuroimaging modalities.
As more and more studies test the predictive value of plasma NfL, we develop a clearer understanding of what blood levels of this biomarker may indicate. Yet, we still do not fully understand what a rise in plasma NfL levels may signify in the absence of a disease diagnosis. To gain further insight, recent studies have examined the plasma NfL levels and their association with mortality. Higher plasma NfL levels are associated with mortality in patients with stroke [14, 15], sporadic Creutzfeldt-Jakob disease , and spontaneous subarachnoid and intracerebral hemorrhages [17, 18]. In these patient cohorts, brain injury initiates neuroaxonal damage which would lead to the release of neuronal NfL into the peripheral blood . Higher NfL in these cases of acute neuroaxonal damage predicts short-term mortality as most of these studies had ~ 30-day follow-up periods, with the exception of sporadic Creutzfeldt-Jakob disease which had a longer follow-up period (mean 14.8 months) . In these instances of acute brain injury, plasma NfL has the utility to predict short-term mortality. However, little evidence exists as to whether NfL is associated with future risk of all-cause mortality in the general population. A recent report found that plasma NfL was associated with mortality in centenarians and nonagenarians, with no reported sex differences in these associations or in plasma NfL levels . Data from the MEMO study of elderly adults (N = 386) found that serum NfL levels were associated with all-cause mortality . In this study, the adjusted hazard ratio obtained from Cox proportional hazards analysis was stronger among men compared with that among women. This is in contrast to our study where significant associations between plasma NfL and mortality were only observed in women. These differences may be multi-factorial in that samples from the participants in the MEMO study were from ~ 1997 and the participants were White and leaner (mean BMI = 27.7) with a mean age of 73 years living in Augsburg, Germany. Our study consisted of African American and White adults with a mean age of 47.75 and a higher BMI (mean = 30.3). Nevertheless, in these population-based studies, blood-based NfL was associated with mortality in women. These studies shed light that a blood assay predictor of dementia is also associated with all-cause mortality in the absence of neurological disease. This suggests that plasma/serum NfL may be used as a clinical biomarker to identify individuals at high mortality risk.
Moreover, a recent analysis of the same cohort as this current study, examined the associations of BMI, allostatic load index (ALindex and its continuous components) with change in plasma NfL over time. This analysis indicated that HbA1c plays an important role in mediating the associations of BMI and the ALindex with δNfL . This partly explains our finding that HbA1c may interact with change in NfL over time to determine mortality risk. Another finding with continuous AL components indicated that NfL at v1 may be a prognostic indicator in women at normal levels of hsCRP close to the mean and that its value diminishes at higher hsCRP levels. This is a novel finding that was not found in previous studies and requires further replication in larger samples. Our findings among men indicate that when the annual rate of change in NfL is set at its mean, both BMI and hypercholesterolemia may consistently become protective against all-cause mortality. This finding did not apply to women and needs to be replicated in larger studies of the comparable population but with longer follow-up periods.
Here, we have analyzed the longitudinal levels of plasma NfL. Very few studies have examined the longitudinal changes in the NfL levels in relation to health outcomes. Moreover, our study consisted of community-dwelling non-demented participants who are racially and socio-economically diverse. This scope broadens the current knowledge about the usefulness of plasma NfL in the population and also in the absence of neurological disease. Both a strength and a limitation of our study are the age of our study participants. In general, this middle-aged cohort is younger with a lower mean NfL at baseline compared to other studies, which typically are in older adults. However, given that we do observe significant associations with mortality in this younger population leads credence to the weathering and accelerated aging phenotype and premature mortality experienced by urban-dwelling African American and White adults living in Baltimore, MD [53,54,55]. Therefore, this cohort represents opportunities to identify early biomarkers of mortality that may identify middle-aged individuals that are at risk for premature mortality. This is especially relevant given the current mortality statistics showing a rise in midlife mortality . Because of the baseline age of our cohort and the short follow-up, we are limited in the number of deaths, which may reduce the statistical power to identify additional associations with mortality. Moreover, the small sample size coupled with a limited number of deaths at follow-up reduced statistical power to detect interactive effects by binary effect modifiers with a relatively low prevalence, such as self-reported cardiovascular disease, in the four-way decomposition models, particularly when stratified by sex. Thus, the extrapolation of our findings to the larger target population may be limited, pending further studies with longer follow-up time after the last measure of plasma NfL. Moreover, selection bias due to unequal distribution by factors that were not accounted for compared to the initial sample is a possibility, explaining in part the sex-specific finding for NfL and all-cause mortality. We used 2-stage Heckman selection models to adjust for selection by key socio-demographic factors, including age, sex, race, and poverty status. However, further analysis suggested that women who reported cardio-vascular disease were selected to a greater extent compared to men, which explains the association between NfL change and mortality being restricted to women. Moreover, the role of chance cannot be ruled with multiplicity in hypotheses and mediators as well as effect modifiers. Finally, residual confounding cannot be discounted in explaining some of our key findings.
In conclusion, we report that plasma NfL levels measured both at baseline and over time can predict all-cause mortality in women. These findings merit further investigation in larger comparable samples of middle-aged adults, which will add to its usefulness as a potential prognostic marker at varying degrees of cardio-metabolic risk, particularly in terms of HbA1c and hsCRP levels. Nevertheless, it is important to identify biomarkers of mortality in middle-aged adults given the alarming recent increase in midlife all-cause mortality . The sex-specific predictive value needs to be examined further both by attempting to replicate these findings in other observational studies and through uncovering the mechanisms behind these differences.
Availability of data and materials
Upon request, data can be made available to researchers with approved proposals, after they have agreed to confidentiality as required by our IRB. Policies are publicized on https://handls.nih.gov. Data access request can be sent to principal investigators (PIs) or the study manager, Jennifer Norbeck at firstname.lastname@example.org. These data are owned by the National Institute on Aging at the NIH. The PIs have made those data restricted to the public for two main reasons: “(1) the study collects medical, psychological, cognitive, and psychosocial information on racial and poverty differences that could be misconstrued or willfully manipulated to promote racial discrimination, and (2) although the sample is fairly large, there are sufficient identifiers that the PIs cannot guarantee absolute confidentiality for every participant as we have stated in acquiring our confidentiality certificate.” Codebook and statistical analysis script can be readily obtained from the corresponding author, upon request, by e-mail contact at email@example.com.
Apolipoprotein E genotype
Body mass index
Center for Epidemiologic Studies-Depression
Healthy Aging in Neighborhoods of Diversity Across the Life Span
Healthy Eating Index 2010 version
Inverse Mills ratio
Institutional Review Board
Medical research vehicles
Plasma neurofilament light
Ordinary least squares
Raket LL, Kuhnel L, Schmidt E, Blennow K, Zetterberg H, Mattsson-Carlgren N. Utility of plasma neurofilament light and total tau for clinical trials in Alzheimer’s disease. Alzheimers Dement (Amst). 2020;12(1):e12099.
Hoffman PN, Cleveland DW, Griffin JW, Landes PW, Cowan NJ, Price DL. Neurofilament gene expression: a major determinant of axonal caliber. Proc Natl Acad Sci U S A. 1987;84(10):3472–6.
Norgren N, Rosengren L, Stigbrand T. Elevated neurofilament levels in neurological diseases. Brain Res. 2003;987(1):25–31.
Hansson O, Janelidze S, Hall S, Magdalinou N, Lees AJ, Andreasson U, et al. Blood-based NfL: a biomarker for differential diagnosis of parkinsonian disorder. Neurology. 2017;88(10):930–7.
Preische O, Schultz SA, Apel A, Kuhle J, Kaeser SA, Barro C, et al. Serum neurofilament dynamics predicts neurodegeneration and clinical progression in presymptomatic Alzheimer’s disease. Nat Med. 2019;25(2):277–83.
de Wolf F, Ghanbari M, Licher S, McRae-McKee K, Gras L, Weverling GJ, et al. Plasma tau, neurofilament light chain and amyloid-beta levels and risk of dementia; a population-based cohort study. Brain. 2020;143(4):1220–32.
Mattsson N, Cullen NC, Andreasson U, Zetterberg H, Blennow K. Association between longitudinal plasma neurofilament light and neurodegeneration in patients with Alzheimer disease. JAMA Neurol. 2019;76(7):791–9.
Weston PSJ, Poole T, Ryan NS, Nair A, Liang Y, Macpherson K, et al. Serum neurofilament light in familial Alzheimer disease: a marker of early neurodegeneration. Neurology. 2017;89(21):2167–75.
Khalil M, Pirpamer L, Hofer E, Voortman MM, Barro C, Leppert D, et al. Serum neurofilament light levels in normal aging and their association with morphologic brain changes. Nat Commun. 2020;11(1):812.
Scherling CS, Hall T, Berisha F, Klepac K, Karydas A, Coppola G, et al. Cerebrospinal fluid neurofilament concentration reflects disease severity in frontotemporal degeneration. Ann Neurol. 2014;75(1):116–26.
Shahim P, Gren M, Liman V, Andreasson U, Norgren N, Tegner Y, et al. Serum neurofilament light protein predicts clinical outcome in traumatic brain injury. Sci Rep. 2016;6:36791.
Teunissen CE, Dijkstra C, Polman C. Biological markers in CSF and blood for axonal degeneration in multiple sclerosis. Lancet Neurol. 2005;4(1):32–41.
Dewey ME, Saz P. Dementia, cognitive impairment and mortality in persons aged 65 and over living in the community: a systematic review of the literature. Int J Geriatr Psychiatry. 2001;16(8):751–61.
Gendron TF, Badi MK, Heckman MG, Jansen-West KR, Vilanilam GK, Johnson PW, et al. Plasma neurofilament light predicts mortality in patients with stroke. Sci Transl Med. 2020;12(569):eaay1913.
Uphaus T, Bittner S, Groschel S, Steffen F, Muthuraman M, Wasser K, et al. NfL (neurofilament light chain) levels as a predictive marker for long-term outcome after ischemic stroke. Stroke. 2019;50(11):3077–84.
Staffaroni AM, Kramer AO, Casey M, Kang H, Rojas JC, Orru CD, et al. Association of blood and cerebrospinal fluid tau level and other biomarkers with survival time in sporadic Creutzfeldt-Jakob disease. JAMA Neurol. 2019;76(8):969–77.
Hviid CVB, Gyldenholm T, Lauridsen SV, Hjort N, Hvas AM, Parkner T. Plasma neurofilament light chain is associated with mortality after spontaneous intracerebral hemorrhage. Clin Chem Lab Med. 2020;58(2):261–7.
Hviid CVB, Lauridsen SV, Gyldenholm T, Sunde N, Parkner T, Hvas AM. Plasma neurofilament light chain is associated with poor functional outcome and mortality rate after spontaneous subarachnoid hemorrhage. Transl Stroke Res. 2020;11(4):671–7.
Rubsamen N, Maceski A, Leppert D, Benkert P, Kuhle J, Wiendl H, et al. Serum neurofilament light and tau as prognostic markers for all-cause mortality in the elderly general population-an analysis from the MEMO study. BMC Med. 2021;19(1):38.
Kaeser SA, Lehallier B, Thinggaard M, Häsler LM, Apel A, Bergmann C, et al. A neuronal blood marker is associated with mortality in old age. Nat Aging. 2021;1:218–25.
Lemaitre JF, Ronget V, Tidiere M, Allaine D, Berger V, Cohas A, et al. Sex differences in adult lifespan and aging rates of mortality across wild mammals. Proc Natl Acad Sci U S A. 2020;117(15):8546–53.
Manouchehrinia A, Piehl F, Hillert J, Kuhle J, Alfredsson L, Olsson T, et al. Confounding effect of blood volume and body mass index on blood neurofilament light chain levels. Ann Clin Transl Neurol. 2020;7(1):139–43.
Akamine S, Marutani N, Kanayama D, Gotoh S, Maruyama R, Yanagida K, et al. Renal function is associated with blood neurofilament light chain level in older adults. Sci Rep. 2020;10(1):20350.
Beydoun MA, Noren Hooten N, Maldonado AI, Beydoun HA, Weiss J, Evans MK, et al. Body mass index and allostatic load are directly associated with longitudinal increase in plasma neurofilament light among urban middle-aged adults. J Nutr. 2021; In press.
Evans MK, Lepkowski JM, Powe NR, LaVeist T, Kuczmarski MF, Zonderman AB. Healthy Aging in Neighborhoods of Diversity Across the Life Span (HANDLS): overcoming barriers to implementing a longitudinal, epidemiologic, urban study of health, race, and socioeconomic status. Ethn Dis. 2010;20(3):267–75.
Waldstein SR, Dore GA, Davatzikos C, Katzel LI, Gullapalli R, Seliger SL, et al. Differential associations of socioeconomic status with global brain volumes and white matter lesions in African American and White adults: the HANDLS SCAN Study. Psychosom Med. 2017;79(3):327–35.
Beydoun MA, Hossain S, Chitrala KN, Tajuddin SM, Beydoun HA, Evans MK, et al. Association between epigenetic age acceleration and depressive symptoms in a prospective cohort study of urban-dwelling adults. J Affect Disord. 2019;257:64–73.
Beydoun MA, Shaked D, Tajuddin SM, Weiss J, Evans MK, Zonderman AB. Accelerated epigenetic age and cognitive decline among urban-dwelling adults. Neurology. 2020;94(6):e613–25.
Tajuddin SM, Hernandez DG, Chen BH, Noren Hooten N, Mode NA, Nalls MA, et al. Novel age-associated DNA methylation changes and epigenetic age acceleration in middle-aged African Americans and whites. Clin Epigenetics. 2019;11(1):119.
Beydoun MA, Noren Hooten N, Beydoun HA, Maldonado AI, Weiss J, Evans MK, et al. Plasma neurofilament light as a potential biomarker for cognitive decline in a longitudinal study of middle-aged urban adults. Transl Psychiatry. 2021;11(1):436.
Blackwell E, de Leon CF, Miller GE. Applying mixed regression models to the analysis of repeated-measures data in psychosomatic medicine. Psychosom Med. 2006;68(6):870–8.
Centers for Disease Control and Prevention. Instructions for completing the cause-of-death section of the death certificate. Hyattsville: National Center for Health Statistics; 2004.
Seeman T, Merkin SS, Crimmins E, Koretz B, Charette S, Karlamangla A. Education, income and ethnic differences in cumulative biological risk profiles in a national sample of US adults: NHANES III (1988-1994). Soc Sci Med. 2008;66(1):72–87.
Alberti KG, Zimmet PZ. Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med. 1998;15(7):539–53.
Visser M, Kritchevsky SB, Newman AB, Goodpaster BH, Tylavsky FA, Nevitt MC, et al. Lower serum albumin concentration and change in muscle mass: the Health, Aging and Body Composition Study. Am J Clin Nutr. 2005;82(3):531–7.
Ridker PM. Cardiology patient page. C-reactive protein: a simple test to help predict risk of heart attack and stroke. Circulation. 2003;108(12):e81–5.
Expert Panel on Detection E, Treatment of High Blood Cholesterol in A. Executive summary of the third report of the National Cholesterol Education Program (NCEP) Expert Panel on detection, evaluation, and treatment of high blood cholesterol in adults (Adult Treatment Panel III). JAMA. 2001;285(19):2486–97.
Golden S, Boulware LE, Berkenblit G, Brancati F, Chander G, Marinopoulos S, et al. Use of glycated hemoglobin and microalbuminuria in the monitoring of diabetes mellitus. Evid Rep Technol Assess. 2003;(84):1–6. https://pubmed.ncbi.nlm.nih.gov/15040142/.
Osei K, Rhinesmith S, Gaillard T, Schuster D. Is glycosylated hemoglobin A1c a surrogate for metabolic syndrome in nondiabetic, first-degree relatives of African-American patients with type 2 diabetes? J Clin Endocrinol Metab. 2003;88(10):4596–601.
Seccareccia F, Pannozzo F, Dima F, Minoprio A, Menditto A, Lo Noce C, et al. Heart rate as a predictor of mortality: the MATISS project. Am J Public Health. 2001;91(8):1258–63.
Lenfant C, Chobanian AV, Jones DW, Roccella EJ, Joint National Committee on the Prevention DE, Treatment of High Blood P. Seventh report of the Joint National Committee on the Prevention, Detection, Evaluation, and Treatment of High Blood Pressure (JNC 7): resetting the hypertension sails. Hypertension. 2003;41(6):1178–9.
Bureau UC. US Census Bureau, Social, Economic, and Housing Statistics Division. Poverty thresholds, 2004. 2004. Accessed Jan 2016.
Beydoun MA, Hossain S, Beydoun HA, Shaked D, Weiss J, Evans MK, et al. Red cell distribution width is directly associated with poor cognitive performance among nonanemic, middle-aged, urban adults. J Nutr. 2020;150(1):128–39.
STATA. Statistics/data analysis: release 16.0. Texas: Stata Corporation; 2019.
Beydoun HA, Huang S, Beydoun MA, Hossain S, Zonderman AB. Mediating-moderating effect of allostatic load on the association between dietary approaches to stop hypertension diet and all-cause and cause-specific mortality: 2001-2010 National Health and Nutrition Examination Surveys. Nutrients. 2019;11(10):2311.
Beydoun MA, Beydoun HA, Mode N, Dore GA, Canas JA, Eid SM, et al. Racial disparities in adult all-cause and cause-specific mortality among us adults: mediating and moderating factors. BMC Public Health. 2016;16(1):1113.
Discacciati A, Bellavia A, Lee JJ, Mazumdar M, Valeri L. Med4way: a Stata command to investigate mediating and interactive mechanisms using the four-way effect decomposition. Int J Epidemiol. 2018. https://doi.org/10.1093/ije/dyy236.
Marshall A, Altman DG, Holder RL, Royston P. Combining estimates of interest in prognostic modelling studies after multiple imputation: current practice and guidelines. BMC Med Res Methodol. 2009;9:57.
Heckman JJ. Sample selection bias as a specification error. Econometrica. 1979;47:153–61.
Weston PSJ, Poole T, O’Connor A, Heslegrave A, Ryan NS, Liang Y, et al. Longitudinal measurement of serum neurofilament light in presymptomatic familial Alzheimer’s disease. Alzheimers Res Ther. 2019;11(1):19.
Sanchez-Valle R, Heslegrave A, Foiani MS, Bosch B, Antonell A, Balasa M, et al. Serum neurofilament light levels correlate with severity measures and neurodegeneration markers in autosomal dominant Alzheimer’s disease. Alzheimers Res Ther. 2018;10(1):113.
He L, Morley JE, Aggarwal G, Nguyen AD, Vellas B, de Souto Barreto P, et al. Plasma neurofilament light chain is associated with cognitive decline in non-dementia older adults. Sci Rep. 2021;11(1):13394.
Zonderman AB, Costa PT Jr, McCrae RR. Depression as a risk for cancer morbidity and mortality in a nationally representative sample. JAMA. 1989;262(9):1191–5.
Zonderman AB, Mode NA, Ejiogu N, Evans MK. Race and poverty status as a risk for overall mortality in community-dwelling middle-aged adults. JAMA Intern Med. 2016;176(9):1394–5.
Noren Hooten N, Pacheco NL, Smith JT, Evans MK. The accelerated aging phenotype: the role of race and social determinants of health on aging. Ageing Res Rev. 2022;73:101536.
Woolf SH, Chapman DA, Buchanich JM, Bobby KJ, Zimmerman EB, Blackburn SM. Changes in midlife death rates across racial and ethnic groups in the United States: systematic analysis of vital statistics. BMJ. 2018;362:k3096.
The authors would like to thank Ms. Nicolle Mode for her contribution in selecting participants for plasma NfL analyses and related data management. The authors would also like to thank all HANDLS participants, staff, and investigators, as well as internal reviewers of the manuscript at NIA/NIH/IRP.
Open Access funding provided by the National Institutes of Health (NIH). This work was supported in part by the Intramural Research Program of the NIH, National Institute on Aging, and National Institutes of Health project number AG000513.
Ethics approval and consent to participate
The HANDLS study is approved by the Institutional Review Board of the National Institutes of Health, National Institute of Environmental Health Sciences. Participants provided written informed consent.
Consent for publication
The views expressed in this article are those of the authors and do not necessarily reflect the official policy or position of the Fort Belvoir Community Hospital, the Defense Health Agency, the Department of Defense, or the US Government. Reference to any commercial products within this publication does not create or imply any endorsement by the Fort Belvoir Community Hospital, the Defense Health Agency, the Department of Defense, or the US government.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
MAB had full access to the data used in this manuscript and completed all the statistical analyses.
Participant flowchart illustrating plasma NfL measurements and mortality for this study. Abbreviations: δ = Annualized change; HANDLS = Healthy Aging Neighborhoods of Diversity Across the Life Span; NfL = Neurofilament Light Chain. Method S1. NfL sample selection. Table S1. Allostatic load indicator criteria. Method S2. Mixed-effects regression models. Fig. S2. Observed vs. empirical bayes estimator for annualized rate of change in Loge transformed NfL, Pearson’s r = 0.83, p < 0.001. Table S2. BMI, AL, AL continuous parameters and cardio-metabolic co-morbidity indices and their relation to all-cause mortality by sex: mediating and moderating effects of NfLv1 using 4-way decomposition a,b,c. Abbreviations: AL = Allostatic Load; ALB = Albumin; BMI = Body Mass Index; CDE = Controlled Direct Effect; CES-D = Center for Epidemiological Studies-Depression; CHOL = Total cholesterol; CRP = C-reactive protein (high sensitivity), Loge transformed; CVD = Cardiovascular Disease; DBP = Diastolic Blood Pressure; DIAB = Diabetes; HBA1C = Glycated Hemoglobin; HDL = High Density Lipoprotein-Cholesterol; HEI-2010 = Healthy Eating Index-2010 version; HYPERT = Hypertension; HYPERCHOL = Hypercholesterolemia; IM = Interaction, mediated; IR = Interaction, Reference; M = Mediators/Moderator; NfL = Plasma Neurofilament Light Chain, Loge transformed; PM-Pure Mediation; RHR = Resting Heart Rate; SBP = Systolic Blood Pressure; WHR = Waist-Hip Ratio; X = Exposure. a See Methods and Table 1 for definition of each NfL exposure (i.e., NfLv1 and δNfL). All exposures (X) and potential mediators/moderators (M) were z-scored for ease of interpretation, with the exception of binary M (coded as 0/1), namely DIAB, HYPERT, HYPERCHOL, CVD. Control variables were set at their means. b Cox models for which 4-way decomposition was conducted is equivalent to Model 2, Table 2, for continuous exposures, to which M was added and considered as a potential mediator/moderator. Control variables included age at v1, sex, race, poverty status, education, HEI-2010 total score, mean energy intake (kcal/d), current tobacco use, current illicit drug use, the CES-D total score and the inverse mills ratio. c Total effects are beta = Loge(HR)±SE with associated p-values from Cox PH hazards models associated with each exposure of interest. Hazard Ratios (HR) point estimates exponent of beta. 95% CI for HR can be calculated as follows: Lower confidence limit, LCL: exp[LogeHR-1.96SE(LogeHR)], upper confidence limit, UCL: exp[LogeHR+1.96SE(LogeHR)]. Table S3. BMI, AL, AL continuous parameters and cardio-metabolic co-morbidity indices and their relation to all-cause mortality by sex: mediating and moderating effects of δNfL using 4-way decomposition a,b,c. Abbreviations: AL = Allostatic Load; ALB = Albumin; BMI = Body Mass Index; CDE = Controlled Direct Effect; CES-D = Center for Epidemiological Studies-Depression; CHOL = Total cholesterol; CRP = C-reactive protein (high sensitivity), Loge transformed; CVD = Cardiovascular Disease; DBP = Diastolic Blood Pressure; DIAB = Diabetes; HBA1C = Glycated Hemoglobin; HDL = High Density Lipoprotein-Cholesterol; HEI-2010 = Healthy Eating Index-2010 version; HYPERT = Hypertension; HYPERCHOL = Hypercholesterolemia; IM = Interaction, mediated; IR = Interaction, Reference; M = Mediators/Moderator; NfL = Plasma Neurofilament Light Chain, Loge transformed; PM-Pure Mediation; RHR = Resting Heart Rate; SBP = Systolic Blood Pressure; WHR = Waist-Hip Ratio; X = Exposure. a See Methods and Table 1 for definition of each NfL exposure (i.e., NfLv1 and δNfL). All exposures (X) and potential mediators/moderators (M) were z-scored for ease of interpretation, with the exception of binary M (coded as 0/1), namely DIAB, HYPERT, HYPERCHOL, CVD. Control variables were set at their means. b Cox models for which 4-way decomposition was conducted is equivalent to Model 2, Table 2, for continuous exposures, to which M was added and considered as a potential mediator/moderator. Control variables included age at v1, sex, race, poverty status, education, HEI-2010 total score, mean energy intake (kcal/d), current tobacco use, current illicit drug use, the CES-D total score and the inverse mills ratio. c Total effects are beta = Loge(HR)±SE with associated p-values from Cox PH hazards models associated with each exposure of interest. Hazard Ratios (HR) point estimates exponent of beta. 95% CI for HR can be calculated as follows: Lower confidence limit, LCL: exp[LogeHR-1.96SE(LogeHR)], upper confidence limit, UCL: exp[LogeHR+1.96SE(LogeHR)].
About this article
Cite this article
Beydoun, M.A., Noren Hooten, N., Weiss, J. et al. Plasma neurofilament light and its association with all-cause mortality risk among urban middle-aged men and women. BMC Med 20, 218 (2022). https://doi.org/10.1186/s12916-022-02425-x