Summary
Cerebrovascular damage coexists with Alzheimer’s disease (AD) pathology and increases AD risk. However, it is unclear whether endothelial progenitor cells reduce AD risk via cerebrovascular repair. By using the Framingham Heart Study (FHS) offspring cohort, which includes data on different progenitor cells, the incidence of AD dementia, peripheral and cerebrovascular pathologies, and genetic data (n = 1,566), we found that elevated numbers of circulating endothelial progenitor cells with CD34+CD133+ co-expressions had a dose-dependent association with decreased AD risk (HR = 0.67, 95% CI: 0.46-0.96, p = 0.03) after adjusting for age, sex, years of education, and APOE ε4. With stratification, this relationship was only significant among those individuals who had vascular pathologies, especially hypertension (HTN) and cerebral microbleeds (CMB), but not among those individuals who had neither peripheral nor central vascular pathologies. We applied a genome-wide association study (GWAS) and found that the number of CD34+CD133+ cells impacted AD risk depending on the homozygous genotypes of two genes: KIRREL3 rs580382 CC carriers (HR = 0.31, 95% CI: 0.17-0.57, p<0.001), KIRREL3 rs4144611 TT carriers (HR = 0.29, 95% CI: 0.15-0.57, p<0.001), and EXOC6B rs61619102 CC carriers (HR = 0.49, 95% CI: 0.31-0.75, p<0.001) after adjusting for confounders. In contrast, the relationship did not exist in their counterpart genotypes, e.g. KIRREL3 TT/CT or GG/GT carriers and EXOC6B GG/GC carriers. Our findings suggest that circulating CD34+CD133+ endothelial progenitor cells can be therapeutic in reducing AD risk in the presence of cerebrovascular pathology, especially in KIRREL3 and EXOC6B genotype carriers.
Introduction
Consistent with the observation that cerebrovascular pathologies co-exist with neuropathological hallmarks of Alzheimer’s disease (AD), both peripheral and central vascular diseases, including hypertension (HTN) and cerebral microbleeds (CMB), have been shown to increase the risk of late-onset AD.1 Microvascular endothelia, which are composed of blood-facing cells of the blood-brain barrier (BBB) in brain capillary beds, may play an important role in the development of AD.2, 3 Recently, we demonstrated that monomeric C-reactive protein (mCRP) and apolipoprotein E (APOE) protein compete to bind with endothelial membrane surface CD31 protein during peripheral inflammation. 4 Although mCRP-CD31 binding causes shortened microvasculature in the brain, thus leading to AD pathology and cognitive impairment, APOE-CD31 binding in the presence of the APOE2 genotype is protective for brain endothelia and against the pathogenesis of AD. Due to the fact that APOE ε4 carriers are vulnerable to peripheral inflammation, as a risk factor, leading to endothelial dysfunction, cerebrovascular pathology and AD pathogenesis,4 we hypothesized that endothelial progenitor cells, as an antagonizing element, can lessen AD risk in persons with particular AD-associated genotypes.
Blood circulating endothelial progenitor cells have different subtypes and express CD34, CD133, or both proteins. When CD133, a member of a novel family of cell surface glycoproteins, is expressed in hematopoietic-derived CD34+ stem cells, these cells can differentiate into endothelia, vascular smooth muscle cells, and neurons and glia in the brain.5 Some studies have suggested that CD34+CD133+ cells are prognostic biomarkers6 and can be therapeutic for ischemic diseases.7–9 Two studies with a small sample size showed that AD patients had a lower number of circulating CD34+CD133+ cells10 or CD34+ cells11 than controls. However, it is unknown whether CD34+CD133+ cells are related to AD risk, especially in the presence of cerebrovascular pathology.
In this study, we examined the relationship between different circulating endothelial progenitor cells, endothelial mature cells and other stem cells, vascular pathologies and AD/dementia risk in participants of the Framingham Heart Study (FHS) Offspring Generation 2 (Gen 2) cohort who have been longitudinally examined for AD dementia incidence. It was found that only CD34+CD133+ endothelial progenitor cells reduced AD risk. We further conducted a genome-wide association study (GWAS) to identify loci that are associated with AD risk through their interaction with CD34+CD133+ endothelial progenitor cells. We used the RNAseq expression data from the Religious Orders Study/Memory and Aging Project (ROSMAP) study to explore the relationship between the expressions of the identified genes and AD.
Results
Circulating CD34+CD133+ progenitor cells and their associated characteristics in the FHS study
This study included 1,566 participants from the FHS Gen2 cohort with an average age (mean±SD) of 65.93±8.91 years at Exam 8, of which 46.36% were females. The median frequency of circulating CD34+CD133+ progenitor cells was 0.032% mononuclear cells measured at baseline (Exam 8). The participants were then divided into four quartiles based on the circulating CD34+CD133+ frequency (Q1 = 0.002-0.020; Q2 = 0.020-0.032; Q3 = 0.032-0.049 and Q4 = 0.049-0.609). We found that younger (p = 0.007) and female (p = 0.007) participants were more likely to be in the higher CD34+CD133+ frequency quartile (Table 1), with no significant differences in education and APOE ε4 carriers across the four quartiles. Vascular disease rates also did not significantly differ across the four CD34+CD133+ quartiles.
The association between increased CD34+CD133+ endothelial progenitor cells and reduced AD/dementia risk
There were 59 cases of incident AD after exam 8. Interestingly, the incident rates of AD (5.85% vs. 3.94% vs. 2.99% vs. 2.08%, p = 0.036) and all-cause dementia (7.66% vs. 5.56% vs. 4.41% vs. 2.86%, p = 0.018) decreased with an increase in the CD34+CD133+ progenitor cell quartiles (Table 1). The average onset age for AD (p = 0.007) and dementia (p = 0.005) slightly decreased as CD34+CD133+ quartiles increased. Cox proportional hazards regression models were used to study the relationship between CD34+CD133+ cell frequency and AD risk or all-cause dementia risk after adjusting for the covariates (Table 2). We found that CD34+CD133+ progenitor cell frequency (log-transformed) was negatively associated with the risk of AD or dementia (Model 1), and this association remained significant after adjusting for covariates, including age, sex, years of education (Model 2), and the addition of APOE ε4 (Model 3), showing that a higher frequency of circulating CD34+CD133+ progenitor cells decreased AD risk (HR = 0.67, 95% CI = 0.46-0.96, p = 0.03) and dementia risk (HR = 0.65, 95% CI = 0.48-0.90, p = 0.008). Furthermore, Kaplan‒Meier survival curves based on four CD34+CD133+ cell quartiles revealed that there were significant differences in AD-free probability among the four CD34+CD133+ progenitor cell quartiles (p = 0.034), thus demonstrating a dose-dependent relationship between higher CD34+CD133+ endothelial progenitor cell counts and reduced AD incidence (Figure 1).
Cox proportional hazards regression models were also used to examine the relationships between other progenitor cells, including CD34+, CD34+CD133-, CD34-CD133+, CD34-CD133-, or CD34+KDR+(VEGFR2) cells, and AD dementia risk. In addition, we studied mature cells including CD31+, CD31-, CD31+DIM, and CD31+lymphoid cells for AD risk using the same analyses. Unlike CD34+CD133+ progenitor cells, none of these cell types impacted the risk of AD dementia (Supplement Table 1).
The relationship between CD34+CD133+ endothelial progenitor cells and reduced AD/dementia risk in the presence of peripheral and central vascular diseases
As expected, compared to those participants without any vascular diseases (0.78%), participants with peripheral vascular diseases, including coronary heart disease (CHD) (8.27%, p = 4.1e-05) and HTN (7.42%, p = 7.5e-05), and participants with cerebrovascular diseases, including clinical diagnoses of stroke (10.83%, p = 3.3e-06), silent infarct (10.97%, p = 1.8e-06), CMB (15.12%, p = 1.8e-08), and high white matter hyperintensities (WMHI) (7.59%, p = 7.0e-05), had significantly higher rates of incident AD (Figure 2A). A similar but stronger trend was observed for the incident rates of all-cause dementia in participants with these vascular diseases, although this can be primarily attributed to the fact that there were more cases of dementia than AD (Figure 2B). When each peripheral vascular or cerebrovascular disease was added as a covariate into the cox proportion regression model, in addition to adjusting for age, sex, years of education and APOE ε4 (Model 3 plus), the negative relationship between CD34+CD133+ progenitor cells and the risk of AD dementia remained (Supplement Table 2).
We further hypothesized that endothelial progenitor cells would play roles in reducing AD risk especially in the presence of vascular pathologies when endothelial cells are damaged. Indeed, after stratification, whereas this trend was not observed in individuals without any vascular diseases, increased CD34+CD133+ endothelial progenitor cell frequency was significantly associated with decreased AD risk among participants who had peripheral vascular diseases (HR = 0.62, 95% CI: 0.40-0.94, p = 0.02) or cerebrovascular disease (HR = 0.58, 95% CI: 0.35-0.96, p = 0.03) (Table 2). Similar trends were observed for all-cause dementia risk among those individuals with peripheral vascular diseases (HR = 0.61, 95% CI: 0.43-0.88, p = 0.008) or with cerebrovascular disease (HR = 0.53, 95% CI: 0.35-0.84, p = 0.003). Consistently, in the absence of any vascular disease, increased CD34+CD133+ frequency quartiles were not associated with a change in AD or dementia incidence (Figure 2C and 2D). In contrast, CD34+CD133+ progenitor cell quartiles significantly reduced AD risk in the presence of vascular diseases in a dose-dependent pattern (Figure 2E-P). Specifically, increased CD34+CD133+ progenitor cell quartiles had a significantly dose-dependent association with decreased AD risk among participants who had HTN (p=3.7e-04) (Figure 2G) and CMB (p=0.02) (Figure 2M). Regarding dementia risk, similar trends for increasing CD34+CD133+progenitor cell quartiles to reduce dementia risk were observed in different vascular diseases in addition to HTN and CMB, specifically among those with CHD (p=0.05) (Figure 2F), HTN (p=1.8e-04) (Figure 2H), stroke (p=0.02) (Figure 2J), silent infarct (p=0.01) (Figure 2L), CMB (p=0.02) (Figure 2N) or in a high WMHI level (p=0.03) (Figure 2P).
Interaction effect between KIRREL3 or EXOC6B homozygous genotype and CD34+CD133+ endothelial progenitor cells on AD risk
Studies have shown that brain endothelial damage increases AD pathogenesis in the brain among APOE ε4 carriers,4 therefore, we hypothesized that peripheral circulating CD34+CD133+ endothelial progenitor cells can further reduce AD risk in persons with particular genotypes. To test this hypothesis, we first conducted a genome-wide association study (GWAS) for AD with the interaction term between CD34+CD133+ cell frequency and single nucleotide polymorphism (SNP) as predictor using a logistic regression model. For SNPs passing a genome-wide significance (p < 5.0e-08), a Cox regression model was applied to validate the interaction effects on AD incidence. As depicted in the Manhattan plot, 8 SNPs showed significant interactions with CD34+CD133+ frequency for AD risk (Figure 3A). Seven of the 8 SNPs with a p value <5.0×10-8 were located in gene KIRREL3 on chromosome 11 (two of the 7 SNPs selected for further analysis: rs580382 and rs4144611) and 1 of the 8 SNPs was in gene EXOC6B (chromosome 2, rs61619102) (Figure 3A, Supplement Figure 2). These associations are supported by evidence with surrounding SNPs (Figures 3B and 3C, Supplement Table 3). Notably, for these 3 SNPs, although their main SNP effects were not associated with AD risk (p>0.40), their interactive effects with CD34+CD133 cells were significantly associated with AD after adjusting for age, sex, years of education, APOE ε4, and principle components (PCs) using both logistic and Cox regression analyses (Table 3).
Consistently, in genotype stratification analysis adjusting for age, sex, years of education, APOE ε4, and PCs, we found that CD34+CD133+ frequency (log-transformed and continuous) was negatively associated with AD incidences among the homozygotes of KIRREL3 rs4144611 TT carriers (HR = 0.29, 95% CI: 0.15-0.57, p<0.001), KIRREL3 rs580382 CC carriers (HR = 0.31, 95% CI: 0.17-0.57, p<0.001) and EXOC6B rs61619102 CC carriers (HR = 0.49, 95% CI: 0.31-0.75, p<0.001) (Supplementary Table 4). Further examination of these findings revealed a dose-dependent effect of increased CD34+CD133+ cell quartiles on the reduction in AD risk among KIRREL3 rs4144611 TT or rs580382 CC homozygotes as well as in EXOC6B rs61619102 CC homozygotes (Table 4, Supplementary Table 5). In contrast, these impacts from CD34+CD133+ cells were not present in the counterpart genotypes of these 3 SNPs. We further tested each genotype e.g. KIRREL3 rs4144611 either TG or GG; rs580382 either CT or TT; as well as EXOC6B rs61619102 either GC or GG carriers, using Cox regression model and found that none of them had interactive effects with CD34+CD133+ progenitor cells on AD risk (Supplementary Table 5).
These relationships were also observed in Kaplan‒Meier survival curves. Among the homozygotes of KIRREL3 rs4144611 TT carriers (p=0.003, Figure 4A), KIRREL3 rs580382 CC carriers (p=1.2e-04, Figure 4B), and EXOC6B rs61619102 CC carriers (p=0.001, Figure 4C), higher CD34+CD133+ frequency quartiles corresponded to a lower incidence of AD development and difference in AD-free probability with statistical significance. Again, among their counterpart genotype carriers (KIRREL3 rs4144611 GG+GT or rs580382 TT+CT carriers and EXOC6B rs61619102 GG+GC carriers), there were no such relationships between CD34+CD133+ quartiles and AD or all-cause dementia risk. Additionally, unlike for AD dementia risk, we did not find associations of these genotypes or the interactive effects of these genotypes and CD34+CD133+ cells for any vascular disease in Table 1 (data not shown).
Exploration of KIRREL3 and EXOC6B genes in Alzheimer’s disease
To address why these two genes could interact with CD34+CD133+ endothelial progenitor cells for AD risk, we explored if the mRNA expressions level of gene KIRREL3 and EXOC6B are associated with AD pathology. We first used the RNAseq dataset generated from the dorsolateral prefrontal cortex (DLPFC) region and monocytes in the ROSMAP study. As shown in Figure 5A-a, brain KIRREL3 expression was associated with AD (p=0.02), but peripheral KIRREL3 expression was not. Compared to controls, brain expression level of KIRREL3 was lower in AD brains specifically in the temporal cortex (TCX) (p = 2.7e-05) regions (Figure 5A-b). To investigate the relationship between KIRREL3 SNPs and lower mRNA levels in the AD brain, we investigated the epigenetics of KIRREL3 SNPs and found that the methylation site cg11751545 was consistently negatively correlated with KIRREL3 gene expression in DLPFC (Supplement Figure 2A). Two KIRREL3 genotypes (rs4144611 GG+GT or rs580382 TT+CT) were negatively associated with DNA methylation level of two CpG sites within KIRREL3 (cg11751545; cg04445570; p<0.05) (Figure 3B, Supplement Figure 2B). Based on these results, we hypothesized that the rescue effect of CD34+CD133+ progenitor cells may reduce AD risk only in vulnerable KIRREL3 rs4144611 TT carriers (Table 4 and Figure 4A) and rs580382 CC carriers (Table 4; Figure 4B and Figure 5A-c) because the low brain expression of KIRREL3 is regulated by the high brain methylation of cg11751545 in these carriers.
While peripheral monocyte expression of EXOC6B was significantly and positively associated with AD (p < 0.001), the brain expression of EXOC6B was not associated with AD after Bonferroni correction (Figure 5B-a). Consistently, EXOC6B rs61619102 G allele was positively associated with EXOC6B gene expression as an eQTL compared to C allele in the following two types of peripheral tissues: visceral and subcutaneous adipose tissues (Figure 5B-b). Thus, we reasoned that the impact of EXOC6B polymorphisms on AD was mainly elicited through the peripheral system and that the EXOC6B rs61619102 homozygous genotype (CC) itself (which was associated with a low level of expression in the peripheral system) interacted with circulating CD34+CD133+ progenitor cells to reduce AD risk (Table 4; Figure 4C; Figure 5B-c).
Discussion
As peripheral and central vascular diseases increase AD risk, it is possible that there are some protective factors in vivo antagonizing or rescuing cerebrovascular pathologies for reducing the risk of AD. However, this kind of in vivo protection is largely unclear. To the best of our knowledge, this study was the first to illustrate that different vascular pathologies increased the risk of consequential AD development which was depended on two factors: 1) the presence of vasculature damage and 2) the ability to repair/reconstitute dysfunctional endothelia by circulating CD34+CD133+ progenitor cells. As endothelial repair for peripheral blood vessels is mainly performed by circulating endothelial progenitor cells with the co-expression of CD34 and CD133,12 circulating CD34+CD133+ progenitor cells have the potential to be therapeutic for AD via repairing damaged endothelia in cerebrovasculature. Furthermore, we discovered that the genotypes in two genes, KIRREL3 and EXOC6B, had interactive effects with the number of circulating CD34+CD133+ cells for AD risk. Interestingly, both KIRREL313, 14 and EXOC6 (the high homolog of EXOC6B)15 have been shown to be risk genes for and are involved in autism spectrum disorder (ASD) and AD13–16.
Progenitor cells have different subtypes with different surface CD expressions. Among different progenitor cells we studied, only CD34+CD133+ endothelial progenitor cells showed dose-dependent protective effects against AD development, especially in the presence of HTN and CMB (Table 2, Figure 1). Unlike the endothelial progenitor cells, circulating mature endothelial cells marked with CD31 expression were not associated with AD risk (Supplement Table 1). Supporting our findings, one clinical study (n = 40) found that AD patients had a lower level of circulating CD34+CD133+ progenitor cells than controls.10 Human cord blood CD133+CD34+ cells have been shown to reduce the apoptosis of endothelial and ganglion cells and increase the number of surviving CD31+ mature endothelial cells in the retina after radiation damage.17 CD34 (also known as sialomucin) is a transmembrane protein expressed on early hematopoietic and vascular-associated tissues.18 CD133 is a stem cell marker, and peripheral blood-derived CD133+ cells can further differentiate into hematopoietic, endothelial, and muscle cells in vivo and in vitro.5 Without counting CD133 expression, our study did not find a relationship between CD34+ cells and AD risk (Supplement Table 1), as another cross-sectional study reported that circulating CD34+ cells were higher, not lower, in AD than in controls.19 It is shown that the leukocyte antigen CD34 is localized to the vascular endothelium throughout the human brain,20 suggesting that CD34+ cells are mixed with endothelial progenitors and leukocytes, but CD34+CD133+ cells are endothelial progenitors only. On the other hand, CD133+ donor cells were detected in several vessels near areas of regeneration, where they expressed human VE-cadherin and CD31, which are biomarkers of mature endothelia.21 Based on our data and others, it is likely that CD34+CD133+ progenitor cells are the specific subtype of progenitors for endothelia that repair damaged brain endothelia and rescue AD risk.
CD34+CD133+ progenitor cells were not associated with different vascular diseases themselves (Supplement Table 2), although they were significantly associated with decreased AD risk in the presence of vascular diseases (Table 2). Cerebrovascular diseases are complex diseases involving multiple cell type pathologies in the brain and endothelia is one of them. The role of cerebral microvasculature pathology in AD is likely through dysfunctional endothelia (reviewed by Dalkara and Alarcon-Martinez)22, thus suggesting that endothelial damage/dysfunction which often occurs in cerebral pathology, but not the other cell types in cerebrovascular pathology, is a key element in AD pathogenesis. In the literature, the relationship between CD34+CD133+ cells and cerebrovascular diseases are not consistent. One study shows that autologous circulating endothelial precursor cells may have the potential for the treatment of ischemic diseases.23 Some studies have shown that the relationship between CD34+CD133+ progenitor cells and cerebral small vessel disease (CSVD) is positive,24, 25 but another one shows a negative relationship.26
Our study (Table 1) and others have shown that increased age may alter the availability of or decrease the number of circulating endothelial progenitor cells27–29. Some studies show that aging diseases like cardiovascular diseases30 and diabetes31 and the absence of healthy aging are associated with decreased CD34+CD133+ progenitor cells. On the other hand, an increased number of circulating CD34+CD133+ endothelial progenitor cells is associated with better outcomes of physical function32 and cardiovascular diseases30. Although we found that females have more CD34+CD133+ progenitor cells in the FHS (Table 1), another study found that females have a lower number of CD34+CD133+ progenitor cells than males,33 this discrepancy could be attributed to different age periods or some other unknown factors in the two cohorts.
By using GWAS analyses, we found that CD34+CD133+ cells differentially affected AD protection in the genotypes of two genes (KIRREL3 and EXOC6B) (Figure 2), thus suggesting that circulating endothelial progenitor cells for AD protection work more efficiently in certain genotypes. CD34+CD133+ cells interact with KIRREL3 SNP to reduce AD risk when the expression of KIRREL3 was low in the brain (Figure 5A). One meta-analysis of GWAS with a large number of participants showed that SNPs of KIRREL3 are a genetic risk factor for AD.13, 14 KIRREL3 has been shown to mediate neuron synapse formation and cell adhesion/guidance34–36, as well as cell differentiation.37, 38 The loss of KIRREL3 leads to changed axon organization, as well as male‒male aggression and cognitive impairment, in mice.39–41 Our data suggest that when brain KIRREL3 gene expression is low and vulnerable due to its genotypes or high methylation status, the impact of circulating CD34+CD133+ progenitor cells on decreasing AD risk is increased; in the counterpart genotypes, the KIRREL3 gene expression in the brain is high, and AD risk is not influenced by circulating endothelial CD34+CD133+ progenitor cells.
The peripheral, but not the brain, expression of EXOC6B was significantly and positively associated with AD risk (Figure 5B). In the EXOC6B rs61619102 CC genotype with low peripheral expression level, increased number of CD34+CD133+ progenitor cells were associated with decreased AD risk (Figure 3 and Table 3). In noncarriers of APOE ε4, the EXOC6 gene was found to be associated with AD risk.15 EXOC6B plays a novel role in promoting stem cell progeny differentiation42 and regulates cellular exocyst formation for protein-membrane trafficking and secretion.43 Decreased EXOC6B expression is associated with low insulin secretion from β cells and an increased risk of type 2 diabetes.44 It is plausible that the EXOC6B rs61619102 CC genotype regulates circulating CD34+CD133+ endothelial progenitor cells to rescue AD pathology in the brain.
Interestingly, both KIRREL3 16, 35, 38, 45–50 and EXOC6B 51–53 genes are related to ASD, another cognitive disease with a young age of onset. Like neurodegenerative diseases, ASD is also associated with cerebral hypoperfusion,54 suggesting the presence of common cerebrovascular pathologies in both diseases. Like AD research studies, post-mortem human studies and some animal ASD models have shown brain neuroinflammation, oxidative stress, and changes in BBB integrity.55 The balance between Wnt-β-catenin and Shh pathways controls angiogenesis, barrier genesis, neurodevelopment, CNS morphogenesis, and neuronal guidance that are related to both AD56, 57 and ASD.55
This study had several limitations. The FHS study was based on a non-Hispanic White population and more large cohorts with different ethnicities are needed to replicate the finding that CD34+CD133+ progenitor cells are associated with reduced risk of AD dementia. Moreover, although we had a one-time point of CD34+CD133+ cell measurement, longitudinal measurements are necessary to characterize the change of CD34+CD133+ progenitor cells during aging and the disease process. Future clinical trials with circulating CD34+CD133+ progenitor cells are key to proving that it can be therapeutic for AD prevention. Nevertheless, our study sheds light for potential cell-based therapy by using circulating CD34+CD133+ endothelial progenitors to reduce AD risk, especially among those individuals with vascular diseases and the vulnerable KIRREL3 and EXOC6B genotypes.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Author contributions
Lindsay A. Farrer and Rhoda Au guided the main mission for the U19 project and this study was part of it.
Yixuan Wang, Jinghan Huang, Ting Fang Alvin Ang, Yibo Zhu and Qiushan Tao had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Yixuan Wang and Jinghan Huang conducted the main data analyses for the manuscript.
Yixuan Wang and Wei Qiao Qiu drafted the initial manuscript; Jinghan Huang and Xiaoling Zhang were actively involved in the manuscript writing.
Ting Fang Alvin Ang, Jesse Mez, Michael Alosco, Gerald V. Denis, Anna Belkina, Ashita Gurnani, Mark Ross, Bin Gong, Jingyan Han, Kathryn L. Lunetta, Thor D. Stein, Rhoda Au and Lindsay A. Farrer were critically involved in the data analyses, the manuscript writing and editing.
Xiaoling Zhang and Wei Qiao Qiu designed the study and supervised the data analyses.
Competing interests
The authors declare no competing interests. The sponsor institutes did not play any role in design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Methods
Study samples and participants
Framingham Heart Study (FHS)
FHS is a single-site, multi-generation, community-based, prospective cohort study of health in Framingham, Massachusetts. The current study included 1,566 participants (Meanage = 65.93±8.91; 46.36% female) from the Offspring cohort (Gen 2) with data on genome wide single nucleotide polymorphisms (SNPs) and circulating CD34+CD133+ cells measured at examination 8 (2005-2007) who were followed subsequently for incident AD dementia. Relevant details about this cohort have been previously described.58 159 of them developed incident AD through the period ending in 2019 (Supplement Figure 1). Informed consent was obtained from all participants and the study protocol was approved by the Institutional Review Board of Boston University Medical Campus.
Religious Orders Study/Memory and Aging Project (ROSMAP)
ROSMAP data were used to investigate the functions of the variants at the DNA methylation levels and gene expression levels and to explore their relationships with AD (https://dss.niagads.org/cohorts/religious-orders-study-memory-and-aging-project-rosmap/).
ROSMAP participants have annual blood draws that provide a repository of stored serum, plasma, and cells. Clinical evaluation, self-report, and medication review are used to document medical conditions. Moreover, subjects are neurologically evaluated every year, and a review of all of the ante-mortem data at the time of death leads to a final clinical diagnosis for each participant; specifically, each individual receives a diagnosis of AD, MCI, or no cognitive impairment (NCI). The details of this methodology have been described in the website.
Circulating endothelial progenitor cells count method
Fasting blood samples were collected from participants in order to count circulating progenitor cells and were stored at 2-8 ℃.59 Peripheral blood mononuclear cells (PBMCs) were used for cell phenotyping.60 Briefly, PBMCs were isolated via Ficoll (GE) gradient centrifugation, counted, incubated with FcR blocking antibodies (Miltenyi Biotech), labeled with anti-human KDR PE antibody (R&D Systems), anti-human CD34 FITC antibody (BD Biosciences), anti-human CD133 APC antibody (BD Biosciences). Samples were washed and fixed with 2% paraformaldehyde. Cells were analyzed on BD FACSCalibur cytometer running CellQuest software (BD Biosciences). Data were analyzed using FlowJo 8 (Flowjo, Inc) and positive populations were identified by comparisons of stained samples with samples incubated with isotype-matching IgG antibodies conjugated to FITC and PE. Compensation was performed using single-stain controls. Percentages of CD34+, CD34+CD133+, CD34+CD133-, CD34-CD133+, CD34-CD133-, CD34+, and CD34+KDR+(VEGFR2) cells were quantified in the data. Similarly, anti-human CD31 PE (BD Bioscences) and anti-human CD45 PerCP antibody (BD Biosciences) were used to analyze CD31+, CD31-, CD31+CD45-, CD31+DIM, and CD31+lymphoid cells.
Dementia/AD diagnoses
All of the Gen 2 participants were followed for incident AD and dementia. The Mini-Mental State Examination (MMSE) was administered beginning at the time of the 5th health exam (i.e., 1991-1995) to monitor changes in cognitive status. A decrease in MMSE performance of 3 or more points from the immediately preceding exam or 5 or more points across all of the exams would indicate a change in cognitive status that warranted review by a dementia diagnostic panel consisting of at least one neurologist and one neuropsychologist. Furthermore, from 1999-2005, all of the surviving Gen 2 participants were invited for an in-depth cognitive examination, which also screened for incident cognitive impairment that warranted review by the dementia diagnostic panel. Consensus diagnostic procedures according to the criteria of Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) and National Institute of Neurological and Communicative Diseases and Stroke/Alzheimer’s Disease and Related Disorders Association (NINCDS–ADRDA) for dementia and AD were conducted and have been previously described.61 Incidences of dementia and AD after the CD34+CD133+ measurement at Exam 8 were used for the analyses.
Peripheral vascular diseases and cerebrovascular diseases
Coronary Heart Disease (CHD)
Participants were diagnosed as having developed coronary heart disease (CHD) if upon review of the case, a panel of three investigators (The Framingham Endpoint Review Committee) agreed on one of the following definite manifestations: myocardial infarction, coronary insufficiency, angina pectoris, sudden death from CHD, or nonsudden death from CHD.
Hypertension (HTN)
We used treatment for HTN as a proxy for having HTN. Treatment for HTN was based on evidence of prescribed hypertension medication use, self-report or FHS M.D. report of treatment for hypertension at the time of the exam. If a discrepancy between these 3 sources of information occurred, a chart review of medical records was conducted to make a final determination of HTN treatment.
Stroke
The diagnosis of stroke was based on the occurrence of a clinically evident stroke documented by clinical records that was reviewed by at least two neurologists. Stroke was defined as the sudden or rapid onset of a focal neurologic deficit persisting for greater than 24 hours.
Silent Brain Infarcts and Cerebral Microbleeds (CMB)
MRI infarcts in persons without a clinical stroke documented by the Framingham Stroke surveillance team, were considered silent brain infarct62 and cerebral microbleeds (CMB) were examined by visual inspection of brain MRI scans as have been previously described63. Exclusion criteria for the brain MRI scan included those who could not undergo MRI due to claustrophobia, presence of a pacemaker or any other source of metal in their body.
White Matter Hyperintensity (WMH)
White matter hyperintensity (WMH) volumetric measures were extracted from a combination of FLAIR and 3D T1 images with a modified Bayesian probability structure based on a previously published method of histogram fitting.64 Volumes are log-transformed to normalize population variance. Details have been described previously.64–68 We dichotomized WMH volume into high and low levels using the median as the cutoff.
Genome-wide association studies (GWAS)
To study the genotypes that were potentially impacted by circulating CD34+CD133+ progenitor cells for AD, genome-wide association studies (GWAS) were performed using the Framingham analytical pipeline. We conducted an interaction GWAS with log transformed CD34+CD133+ cell frequency by using a logistic regression model for AD. The relationship between AD and interaction terms (CD34+CD133+ frequency and genome-wide SNP dosage) was tested using GEEPACK (logistic regression utilizing generalized estimating equations) while adjusting for age, sex, years of education, and the first 10 population principal components (PCs). Additive models were assumed. Only autosomal SNPs on chromosomes 1-22 were considered. A statically filtered set of dosage files for MACH imputation via the 1,000 Genomes reference files from March 2012 of only samples of European ancestry was used. The filtering criterion was r2 ≥0.1. Moreover, only SNPs with minor allele frequency (MAF) ≥10% were chosen for the analysis. In addition, samples with a genotyping rate less than 97% and with an excess number of heterozygote observations or Mendelian errors were removed. Manhattan plot, QQ plot, and genomic control69 were used for visualization and quality control, and LocusZoom70 was used to present the regional information. Only p values <5.0×10-8 were considered to be statistically significant for the SNP-CD34+CD133+ progenitor cell interactive effects for AD.
Gene expression analysis and DNA methylation
We additionally explored the gene expression patterns across different AD brain regions by using Agora (https://agora.adknowledgeportal.org/genes), which included 9 brain regions, e.g. the anterior cingulate cortex (ACC), cerebellum (CBE), dorsolateral prefrontal cortex (DLPFC), frontal pole (FP), inferior frontal gyrus (IFG), posterior cingulate cortex (PCC), parahippocampal gyrus (PHG), superior temporal gyrus (STG), and temporal cortex (TCX). Two AD-related traits, e.g., BRAAK (neurofibrillary tangles) and CERAD (neuritic plaques), were also measured. In addition, GTEx Portal (https://gtexportal.org/home/) was used to study the associations between the selected SNPs and gene expression (eQTLs). The details of these procedures have been previously described.71
We used Brain xQTLServe from ROSMAP (http://mostafavilab.stat.ubc.ca/xQTLServe/) to explore the associations between SNPs and DNA methylation (mQTLs). Genotype data were generated from 2,093 individuals of European descent. Of these individuals, DNA methylation (450K Illumina array) (n=468) was derived from postmortem frozen samples of a single cortical region, DLPFC. Gene expression data were generated by using RNA-sequencing (Illumina HiSeq) from the DLPFC of 494 individuals at an average sequence depth of 90 M reads. The details of this methodology have been previously described.72
Statistical analysis
Data were analyzed by using R version 3.6.2. The Cox models were applied to circulating cells including CD34+CD133+, CD34+, CD34+CD133-, CD34-CD133+, CD34-CD133- or CD34+KDR+(VEGFR2), CD31+, CD31-, CD31+DIM and CD31+lymphoid cells to examine the relationship between CD34+CD133+ or quartiles and dementia/AD diagnosis risk with adjustments for age, sex, years of education, APOE ε4, and different vascular diseases (Supplement Table 1). Only CD34+CD133+ progenitor cells were found to be associated with AD dementia.
Circulating CD34+CD133+ progenitor cells were first divided into four quartiles based on cell numbers. Quartiles of cell numbers were compared on the number of subjects, age at baseline, sex, years of education, APOE ε4, incident AD/dementia status, CHD, HTN and different cerebrovascular diseases using analysis of variance (ANOVA), chi-squared tests and Kaplan‒Meier survival analysis. We further stratified the subjects into the absence and the presence of different vascular pathologies and examined the relationship between CD34+CD133+ quartiles and dementia/AD diagnosis risk by using the chi-squared test in each subgroup. Furthermore, the Cox model was applied to the stratification of 1) no vascular diseases, 2) peripheral vascular diseases, and 3) cerebrovascular diseases with adjustments for age, sex, years of education and APOE ε4. Additionally, the main effects and interaction effects (with CD34+CD133+) of GWAS-selected SNPs were obtained by using logistic regression and Cox proportional hazard regression adjusting for age, sex, years of education, APOE ε4, and PCs. These interaction effects were further assessed by applying stratification analyses of genotypes with Cox proportional hazard regression models and a Kaplan‒Meier survival analysis.
Acknowledgments
This study was supported by National Institute on Aging grants U19-AG068753, RF1-AG057519, and R01-AG048927. The FHS data collection was supported by the National Heart, Lung, and Blood Institute contract (N01-HC-25195). We want to express our thanks to the FHS participants for their decades of dedication and to the FHS staff for their hard work in collecting and preparing the data.