Abstract
Background:
Bone morphogenetic protein (BMP) plays important roles in the pathology of Alzheimer’s disease (AD).
Objective:
We sought blood BMP6 involved in the processes underlying cognitive decline and detected them in association with AD.
Methods:
A total of 309 participants in Shanghai Mental Health Center (SMHC) and 547 participants in Alzheimer’s disease Neuroimaging Initiative (ADNI) cohort were included. Blood BMP6 and cognitive functions were measured in all subjects of both cohorts at baseline, and in 482 subjects of ADNI cohort after one year. A total of 300 subjects in ADNI cohort were detected cerebrospinal fluid (CSF) tau biomarker, and 244 received 1-year follow-up.
Results:
AD patients had lower levels of blood BMP6 compared to normal controls, and BMP6 was positively associated with cognitive functions. Longitudinal BMP6 combing with APOE genotype could distinguish probable AD from normal controls. The influence of blood BMP6 on cognition was modulated by tau pathology.
Conclusion:
Blood BMP6 was associated with cognitive performance and identified as a potential predictor for probable AD.
INTRODUCTION
Alzheimer’s disease (AD) is the most common form of dementia worldwide and is characterized by extracellular deposits of amyloid-β (Aβ) plaques and intracellular neurofibrillary tangles of tau protein in brain [1]. Easily accessible biomarkers of AD pathology are essential for detecting individuals with the greatest risk of developing mild cognitive impairment (MCI) or AD. The ratios of cerebrospinal fluid (CSF) total tau/Aβ42 and phosphorylated tau/Aβ42 can predict cognitive decline [2], but require invasive procedures and patient acceptance. Blood-based cytokines associated AD pathology have significant advantages of time efficiency and reduced invasiveness. It has been demonstrated that neither plasma nor serum Aβ1 - 40 and Aβ1 - 42 consistently differ between AD patients and normal controls [3]. Elevated levels of phosphorylated tau is thought to be more specific to AD than total tau [4]. Besides above biomarkers, there are still lack of other valuable and meaningful markers.
Bone morphogenetic protein (BMP), a member of the β-transforming growth factor (TGF-β) subfamily, has important effects on neuronal differentiation and axonal growth [5]. The BMPs, more than 20 at the last count, are discovered in bone tissue, which activate the canonical small mother against decapentaplegic (Smad) pathway in the brain via their type I and type II Serine/Threonine kinase receptors. BMP6, is demonstrated increased in AD pathology, accompanied by impaired neurogenesis [6]. This information implicates that BMP6 has an important role in AD pathology, but there is still lack of research on characteristics of peripheral BMP6 in AD patients. Therefore, we aimed to identify blood BMP6 expression in AD patients and relationship with global cognition. We analyzed the data for blood BMP6 of 309 participants in Shanghai Mental Health Center (SMHC) cohort and 547 participants in Alzheimer’s disease Neuroimaging Initiative (ADNI) cohort. To additionally test the prediction effect of blood BMP6 for AD and whether the influences of blood BMP6 on cognition were modulated by AD core pathology in ADNI cohort.
MATERIALS AND METHODS
309 participants were derived from SMHC cohort, which included 103 probable AD, 92 mild cognitive impairment (MCI), and 114 cognitively normal controls (NC) between March 2011 and April 2018. All participants were Mandarin-speaking Han Chinese and were older than 55 years old. Subjects underwent a screening process that included medical history, physical and neurological examination, and cognitive assessment by a face-to-face interview at baseline. Some were assessed at a 1-year follow-up. All participants had scores on the Hachinski ischemia scale of < 4, and no history of significant systemic or psychiatric conditions or traumatic brain injuries that could compromise brain function. The Beijing version of the Montreal Cognitive Assessment (MoCA) [7] and Petersen Mini-Mental States Examination (MMSE) [8] were used to assess cognitive function. Probable AD was diagnosed according to the criteria from the Diagnostic and Statistical Manual for Mental Disorders IV (DSM IV), and from the National Institute of Neurologic and Communicative Disorders and Stroke/Alzheimer’s Disease and Related Disorders Association (NINCDS/ADRDA) [9]. MCI was diagnosed according to Petersen’s criteria [10]. Based on cognitive test scores, subjects displaying memory deficit or additional deficit in another domain were included in the MCI group [11]. Cognitively normal controls had no history of cognitive decline, neurologic disorders, or uncontrolled systemic medical disorders.
Data collection in SMHC cohort was carried out in accordance with the recommendations of the Shanghai Mental Health Center ethical standards committee on human experimentation.
To validate the findings from our cohort and to explore the causal relationships, data of 482 elderly including 54 NC, 544 MCI, and 84 definite AD was downloaded from ADNI database (https://adni.loni.usc.edu). ADNI is a multi-site dataset designed to test clinical, imaging, genetic, and biochemical biomarkers of AD, which was launched in 2003. Data collection and sharing in ADNI were approved by institutional review boards of all participating institutions, and written informed consent was obtained from all participants or their guardians according to the Declaration of Helsinki. The participants are older adults aged 55–90 years. Each participant underwent an in-person interview for health and neuropsychological assessments at baseline and at annual follow-up. MMSE and Alzheimer’s Disease Assessment Scale-Cognitive section (ADAS-cog) were used to assess cognitive function. The inclusion criteria were as follows: NC: MMSE scores between 24–30, Clinical Dementia Rating (CDR) of 0; MCI subjects: MMSE scores between 24–30, with memory complaint or other cognitive impairments, a CDR of 0.5, and preserved activities of daily living; Mild AD subjects: MMSE scores between 20–26, CDR of 0.5–1.0, meeting NINCDS/ADRDA criteria for probable AD [9], and CSF Aβ < 976.6 pg/ml. In this study, we chose subjects from ADNI database who had blood BMP6 protein measurement at baseline and 1-year follow-up.
Measurement of blood BMP6 and CSF AD biomarkers
In SMHC cohort, peripheral blood samples were collected in the fasting state from all subjects by venipuncture into a coagulation promoting tube. Samples were collected and centrifuged at 3000 g for 20 min at 4°C, and serum was stored at –80°C. Enzyme linked immunosorbent assay kits for BMP6 (HEB037) and noggin (HEN036) (Bogoo Biological, Shanghai, China) were used to determine the serum levels of BMP6. Concentrations were expressed as ng/ml for BMP6.
In ADNI cohort, procedure of plasma protein data collection and measurement was explained in detail elsewhere (https://adni.loni.ucla.edu/wp-content/uploads/2010/11/BC_Plasma_Proteomics_Data_Primer.pdf). In brief, plasma proteins including BMP6 protein were measured in a subset of ethylene diamine tetraacetic acid (EDTA) plasma samples (obtained in the morning following an overnight fast), using a 190-analyte multiplex immunoassay panel. The panel, referred to as the human discovery map, was developed on the Luminex xMAP platform by Rules-Based Medicine (RBM) to contain multiple proteins [12]. Furthermore, CSF AD biomarkers including Aβ42, phosphorylated tau (p-tau), and total tau, were detected. The ADNI used the fully automated and highly standardized Roche Elecsys immunoassay to assess AD biomarkers. Aβ+ individuals had a CSF Aβ1 - 42 < 976.6 pg/ml [13].
APOE genotype
DNA was extracted with the QIAamp®DNA Blood Mini Kit and amplified by the polymerase chain reaction (PCR) with forward primers 14 5’-ACGGCTGTCCAAGGAGCTG-3’ (rs429358) and 5’-CTCCGCGATGCCGATGAC-3’ 15 (rs7412). APOE genotype was performed through Restriction Fragment Length Polymorphism (RFLP) technology.
Statistical analyses
The data was not normally distributed; therefore, statistical significance was assessed using nonparametric tests. Demographics, physical disease, cognitive scores, and blood cytokine were analyzed using Mann-Whitney U test for continuous variables and a χ2 test for categorical variables. BMP6 change rate was defined as the ratio of (follow-up BMP6 - Baseline BMP6)/Baseline BMP6, which reflected the annual change in BMP6 in ADNI cohort. Longitudinal change of BMP6 was plotted against age to see the annual difference between NC and AD groups. The Spearman correlation coefficient was used to explore the associations of blood BMP6 with cognitive scores.
Modulation analysis was used to test whether the association between BMP6 and follow-up cognitive function was modulated by AD core pathology. This analysis was performed using the PROCESS macro for SPSS [14]. Age, education years, sex, and APOE4 genotype were included as covariates. Mean center was used for construction of continuous variables. Significance was determined through 95% bias corrected confidence intervals from bootstrapping of 1000 iterations. Significant interactions were probed using simple slopes analyses, which provided information about the deferential associations between variables, depending on the levels (i.e., high or low) of predictor variables. ‘High’ levels of each variable were computed as one standard deviation above the mean, and ‘low’ levels were computed as zero [15].
The least absolute shrinkage and selection operator (LASSO) regression through R software was applied for selection of BMP6 and demographics data to distinguish NC and AD subjects in ADNI cohort. This machine-learning algorithm identified the variables which could predict a given dependent variable and allowed optimal variable weights for this prediction. Randomly assigned 65% of all subjects, were involved in the train set, and the remaining 35% were involved in the test set. A 5-fold cross-validation procedure served to optimize the penalization and weighted parameters of this LASSO model. We estimated the discriminatory efficacy of this model in the train set by quantifying the area under curve (AUC) of the receiver-operator characteristic (ROC), and verified the generalizability of this model in the test set [16].
The statistical significance of all tests was set at a two-sided p value < 0.05. All analyses were performed using SPSS 17.0 or R version 4.0.3.
RESULTS
Demographic and clinical variables
As for theSMHC cohort, demographics and cognitive scores for a total of 309 subjects including probable AD (n = 103), MCI (n = 92), and NC (n = 114) subjects at baseline were shown Table 1. The mean age of this study sample was 71.02 years, the mean educational years was 10.77, and males accounted for 40.78%. After 1 year, 163 subjects were assessed by MoCA test and 127 were assessed by MMSE test, and some subjects showed changes of cognitive functions including the conversion of NC to MCI, the conversion of MCI to probable AD, and the reversion of MCI to NC (Table 1). As expected, differences in cognitive scores were observed.
Characteristics of participants in SMHC cohort and ADNI cohort
p < 0.05 The significance of difference among groups was examined by Mann-Whitney test (for continuous variable) and Pearson’s Chi-squared test (for categorical variable). BMP6 change rate: (1-year BMP6-Baseline BMP6)/Baseline BMP6.
As for ADNI, a total of 482 subjects including definite AD (n = 84), MCI (n = 344), and NC (n = 54) were involved at baseline (Table 1). Compared with SMHC cohort, the study sample was older (mean = 75.17 years), had more educational years (mean = 15.64 years), and had more males (61.41%). Individuals in definite AD group tended to be APOE4 positive. All subjects received 1-year follow-up, and 425 subjects received 2-year follow-up which included 411 subjects were assessed by MMSE and 409 were assessed by ADAS-cog (Table 1). Some subjects showed changes of cognitive functions including the conversion of the conversion of MCI to AD and the reversion of MCI to NC (Table 1). Similarly, differences in cognitive scores were observed.
Blood levels of BMP6
As for SMHC cohort, blood levels of BMP6 protein were decreased significantly in probable AD than NC subjects (Table 1, Fig. 1). However, in ADNI database, there were no difference of baseline BMP6 levels between NC and definite AD groups. After a year, significant lower BMP6 levels were found in definite AD and MCI groups than NC group (Table 1, Fig. 1). Slopes representing the change of blood BMP6 between 2 years plotted against age. Individual longitudinal changes showed decreased rate of BMP6 in definite AD subjects when compared to NC subjects (Fig. 1E).

There was significant difference of blood BMP6 between NC and probable AD groups in SMHC cohort (A). There was not significant difference of blood BMP6 between NC and definite AD groups at baseline (B1), but significant differences of 1-year BMP6 (B2) and BMP6 change rate (B3) between NC and definite AD groups in ADNI cohort. Correlation plots showed positive correlations between blood BMP6 and cognitive scores including baseline MoCA (C1), baseline MMSE (C2), 1-year MoCA (C3), and 1-year MMSE (C4) in SMHC cohort. Correlation plots showed positive correlation between blood BMP6 change rate and 1-year MMSE (D1), and negative correlation between blood BMP6 change rate and 1-year ADAS-cog (D2) in ADNI cohort. Longitudinal BMP6 plotted against age in definite AD and NC groups between baseline and 1-year follow-up (E) in ADNI cohort.
In ANDI, a feature set including baseline BMP6, 1-year BMP6, APOE4, demographics, diabetes, and hypertension, was used to distinguish definite AD and NC subjects at 1-year and 2-year follow-up. The total sample was randomly divided into the train set (65%) and test set (35%), and feature selection was conducted by LASSO logistic regression. First, BMP6, 1-year BMP6, and the change rate of BMP6 were involved, and 1-year BMP6 was selected by LASSO model, which showed the prediction efficacy with AUC of 56–67% at 1-year follow-up (Fig. 2A1, A2) and 2-year follow-up (Fig. 2F1, F2). Second, APOE4 and demographics as features were involved, and APOE4 was selected by LASSO model which showed better prediction efficacy with AUC of 76–82% at 1-year follow-up (Fig. 2B1, B2) and 2-year follow-up (Fig. 2F1, F2). Last, all features were involved, and an overview of the weighted coefficients of the feature set (Fig. 2C2, G2) and the features selection process (Fig. 2C1, G1) were showed. 1-year BMP6 levels combining with APOE4 displayed favorable prediction efficacy. ROC analysis revealed an AUC of 0.854 (95% confidence interval (CI), 0.789–0.919) in the train set and 0.892 (95% CI, 0.818–0.965) in the test set (Fig. 2D1, D2) after 1 year, and an AUC was 0.832 (95% CI, 0.764–0.899) in the train set and 0.803 (95% CI, 0.711–0.895) in the test set (Fig. 2H1, H2) after two years.

In ADNI cohort, ROC curve analysis revealed the AUC for the diagnostic accuracy of 1-year blood BMP6 (A1, A2) or APOE4 (B1, B2) to categorize definite AD patients from NC controls in train and test sets respectively at 1-year follow-up. ROC curve analysis revealed the AUC for the diagnostic accuracy of 1-year blood BMP6 (E1, E2) or APOE4 (F1, F2) to categorize AD patients from NC controls in train and test sets respectively at 2-year follow-up. The vertical dotted line points to the optimal λ value and the number of optimal predictors, and the pathway of coefficients among all parameters at 1-year follow-up (C1, C2) and at 2-year follow-up (G1, G2) respectively. The parameters including demographics, APOE4, baseline BMP6, 1-year BMP6, and BMP6 change rate were selected by LASSO model. In the train and test sets, ROC curve analysis revealed the AUC for the diagnostic accuracy of 1-year blood BMP6 combing APOE4 to categorize AD patients from NC controls at 1-year follow-up (D1, D2) and 2-year follow-up (H1, H2) respectively.
Association of blood BMP6 with cognitive functions
As SMHC cohort, we observed a significant correlation between baseline BMP6 and cognitive functions including MoCA and MMSE scores at baseline and 1 year follow-up (Fig. 1C1-C4). Similarly, as for ADNI cohort, a significant correlation between BMP6 change rate and cognitive functions including MMSE and ADAS-cog scores at 1-year follow-up was observed (Fig. 1D1, D2).
Moderator analyses
Blood BMP6 was a significant risk factor for cognitive impairment, and we next investigated whether BMP6 contributed to cognitive impairments via modulating tau pathology. As for ADNI cohort, a total of 300 subjects with both blood BMP6 and CSF tau tests were selected, which included 53 NC, 167 MCI, and 80 definite AD. The characteristics were showed in Table 2. We found the relationship between BMP6 change rate and 1-year global cognition was modulated by tau pathology including baseline tau (Fig. 3A, C), baseline p-tau (Fig. 3B, D), 1-year tau (Fig. 3E, G), and 1-year p-tau (Fig. 3F, H). First, the subjects with lower BMP6 change rate tended to have lower follow-up MMSE and ADAS-cog scores. Second, subjects with lower BMP6 change rate, higher baseline, or follow-up tau burden, tended to report lower follow-up MMSE and ADAS-cog scores.
Characteristics of participants with blood BMP6 and CSF AD biomarkers in ADNIcohort

In ADNI cohort, the relationship between blood BMP6 change rate and cognitive measures including 1-year global cognition measured by MMSE (A, B, E, F) as well as ADAS-cog (C, D, G, H), was modulated by tau pathology including baseline tau (A, C), baseline p-tau (B, D), 1-year tau (E, F), and 1-year p-tau (G, H). Subjects who experienced lower BMP6 change rate, higher baseline, or follow-up tau burden, reported lower follow-up cognitive function.
DISCUSSION
The major findings of this study were that 1) AD patients had lower levels of blood BMP6 compared to normal controls; 2) BMP6 was positively associated with baseline or follow-up cognitive functions; 3) Longitudinal BMP6 combining with APOE4 could distinguish AD from NC elderly; 4) the influence of BMP6 change rate on cognition was modulated by tau pathology. These findings consolidated the close relationships of BMP6 with cognitive function and tau pathology, supporting the hypothesis that BMP6 was identified as a potential predictor of AD.
The focus on blood-based AD biomarkers had grown exponentially during the past decade. A series of studies had applied plasma profiling to detect AD from control subjects, which were related to the inflammatory response, lipid metabolism, and the microcirculation. These measurements showed substantial promise, but it remained uncertain whether these protein profiles had a definitive correlation with brain function [17]. In the present study, BMP6, as a member of TGF-β subfamily, influenced brain development and neurogenesis, and had a great impact on brain degenerative diseases, particularly AD [18]. Previous studies found BMP6 significantly increased in the hippocampus of AD patients and transgenic mouse model of familial AD [6]. The increase in BMP6 had been reported during healthy aging in wild-type mice, and the process might be exacerbated in pathological conditions such as AD that led to Aβ accumulation, microglia activation, and unfavorable inflammatory microenvironment [19, 20]. At present, there is lack of research on the expression of peripheral BMP6 in AD patients. We found a decrease of blood BMP6 both in probable AD of SMHC cohort and definite AD of ADNI cohort, compared to normal controls. As for ADNI cohort, BMP6 levels in definite AD did not differ from NC at baseline, but reduced significantly at 1-year follow-up, while in SMHC cohort, BMP6 levels in probable AD decreased significantly at baseline. The possible reason was that AD patients involved in ADNI cohort were mild, but severe in SMHC cohort. The changing pattern of BMP6 in peripheral blood and brain tissue was opposite, which was similar to Aβ in AD pathology. The possible reason was that BMP6 existed as hydrophobic homodimer/heterodimer in brain tissue and as hydrophilic monomer in body fluids [21]. Furthermore, the baseline levels of BMP6 in SMHC cohort and the changing rates of BMP6 in ADNI cohort, were both associated with global cognition. Above findings suggested that blood BMP6 had lower expression in AD patients and had a definitive correlation with cognitive function.
Previous research found APOE4 participants had 3.34 times the odds of developing AD within 17 years than APOE ɛ4- participants [22]. In ADNI cohort, we found APOE4 alone could distinguish definite AD from normal elderly with the AUC > 75%, and blood 1-year BMP6 alone with the AUC of 56–67%. However, blood 1-year BMP6 combining with APOE4 had diagnostic and prediction for definite AD with the AUC > 85% at 1-year follow-up and AUC > 80% at 2-year follow-up, which suggested longitudinal BMP6 combining APOE4 as a potential diagnostic and predictive biomarker of cognitive decline in AD. These findings were thus less driven by the complex status of dementia. Although the dementia patients involved at baseline were diagnosed as AD through CSF AD biomarkers, dementia patients transforming from MCI were not assessed by CSF.
In ADNI cohort, we revealed that increased blood BMP6 contributed to cognitive impairment under a modulation of tau pathology. The elderly who experienced decreased BMP6 rate and high tau burden showed the lowest levels of cognitive function. As one of AD characteristics, tau protein was the most abundant microtubule-associated protein in the brain, and BMP/TGF-β pathway was not commonly discussed in relation to tau pathology [23]. Previous publications found the link between BMP and tau protein. Wu et al. revealed that BMP2 might affect injured facial nerve regeneration and levels of tau protein [24]. Lauzon et al. developed a small peptide derived from BMP9, which could inactivate GSK3β, a tau kinase [25]. The above research supported our finding that tau pathology could modulate the influence of blood BMP6 on cognitive impairment.
There were limitations in this study. First, the study sample of definite AD and NC was limited in ADNI cohort, which might introduce population heterogeneity bias. Second, CSF AD biomarkers and APOE genotype were adapted in ADNI cohort, but not in SMHC cohort, which made the loss of the validation about the modulation and ROC analysis. Third, dementia patients transforming from MCI were not assessed by CSF, which made the complex status of dementia at follow-up.
Conclusions
The present study indicates that AD patients have lower levels of blood BMP6 compared to normal controls, and BMP6 was positively associated with global cognition. Longitudinal BMP6 combing with APOE4 could distinguish AD from NC elderly. The influence of BMP6 change rate on cognition was modulated by tau pathology. These findings demonstrate the relationship of blood BMP6 with cognitive performance and AD diagnosis.
Footnotes
ACKNOWLEDGMENTS
Data collection and sharing for ADNI data section was funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie, Accelerate Cure/Treatments for Alzheimer’s Disease (ACT-AD), Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; Alector, Inc.; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Denali, Inc.; Diamir, Inc.; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; Fujifilm Toyama chemical Co., Ltd.; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Magqu Co., Ltd.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; PeopleBio; Pfizer Inc.; Piramal Imaging; Servier; Saladax Biomedical; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (
). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.
This research was also supported by NIH grants P30 AG010129 and K01 AG030514. The principal investigators of ADNI includes Michael W. Weiner (UCSF, NCIRE, VA), Paul Aisen (University of Southern California), Laurel Beckett (University of California), Robert Green (Brigham and Women’s Hospital), Clifford Jack (Mayo Clinic, Rochester, Minnesota), William Jagust (University of California, Berkeley), John C. Morris (Washington University), Ronald Petersen (Mayo Clinic), Andrew J. Saykin (Indiana University), Leslie Shaw (University of Pennsylvania), Arthur W. Toga (University of Southern California), John Q. Trojanowski (University of Pennsylvania).
This study was supported by grants of National Natural Science Foundation of China (81301139), Science Foundation of the Chinese Academy of Sciences (XDA1204101), Clinical research center project of Shanghai Mental Health Center (CRC2017ZD02), and National Science and Technology Support Program of China (2009BAI77B03).
