- Open Access
Convergent transcriptomic and genomic evidence supporting a dysregulation of CXCL16 and CCL5 in Alzheimer’s disease
Alzheimer's Research & Therapy volume 15, Article number: 17 (2023)
Neuroinflammatory factors, especially chemokines, have been widely reported to be involved in the pathogenesis of Alzheimer’s disease (AD). It is unclear how chemokines are altered in AD, and whether dysregulation of chemokines is the cause, or the consequence, of the disease.
We initially screened the transcriptomic profiles of chemokines from publicly available datasets of brain tissues of AD patients and mouse models. Expression alteration of chemokines in the blood from AD patients was also measured to explore whether any chemokine might be used as a potential biomarker for AD. We further analyzed the association between the coding variants of chemokine genes and genetic susceptibility of AD by targeted sequencing of a Han Chinese case–control cohort. Mendelian randomization (MR) was performed to infer the causal association of chemokine dysregulation with AD development.
Three chemokine genes (CCL5, CXCL1, and CXCL16) were consistently upregulated in brain tissues from AD patients and the mouse models and were positively correlated with Aβ and tau pathology in AD mice. Peripheral blood mRNA expression of CXCL16 was upregulated in mild cognitive impairment (MCI) and AD patients, indicating the potential of CXCL16 as a biomarker for AD development. None of the coding variants within any chemokine gene conferred a genetic risk to AD. MR analysis confirmed a causal role of CCL5 dysregulation in AD mediated by trans-regulatory variants.
In summary, we have provided transcriptomic and genomic evidence supporting an active role of dysregulated CXCL16 and CCL5 during AD development.
Alzheimer’s disease (AD) is a progressive and incurable age-related neurodegenerative disease with cognitive decline caused by neuronal loss and brain atrophy [1, 2]. The neuropathological features of AD include the presence of abundant extracellular amyloid plaques laden with β-amyloid peptide (Aβ), intraneuronal neurofibrillary tangles formed by the hyperphosphorylated tau, neuritic dystrophy, loss of synapses and neuronal somata [1, 2]. Microglia activation and neuroinflammation are also hallmarks of AD [1, 2]. Genetic studies including genome-wide association studies (GWAS) [3,4,5,6,7], whole-exome sequencing (WES) [8, 9], and whole-genome sequencing (WGS)  studies have identified numerous genomic loci associated with AD, with immune genes being highlighted frequently.
The immune factors that are linked to neuroinflammation may accompany and contribute to neurodegenerative pathology . Chemokines are small proteins with the ability to induce targeted chemotaxis of nearby reactive immune cells and have a well-established role in the immune system [12, 13]. Chemokines are a group of chemotactic cytokines and can be classified into two groups: the homeostatic chemokines and the inflammatory chemokines. The homeostatic chemokines are important for lymphoid organ development and immune cell trafficking. The inflammatory chemokines are actively involved in the mobilization of effector cells to the inflammatory sites . Apart from their roles in the immune system, chemokines also take part in the physiological and pathological processes in the central nervous system (CNS). In the human brain, neurons and glial cells are able to express chemokines and also have chemokine receptors [15, 16]. Chemokines participate in the proliferation, differentiation, and migration of neural cells and are important for brain homeostasis . The expression of chemokines and the functioning of their receptors may change in CNS diseases [14,15,16]. Chemokines have been reported to be actively involved in CNS development and neurological diseases . Chemokine signaling affects a variety of cellular activities and functions, including the migration and survival of neuronal precursors , the migration and proliferation of oligodendrocyte progenitors , the maintenance of oligodendrocyte lineage, myelination, and white matter , the central synaptic transmission , glymphatic function and neuroinflammation , and aging-dependent neuronal regenerative decline .
There are several lines of evidence to show expression change, chemokine/chemokine receptor axis signaling, pathological correlation, and genetic regulation of chemokines, each take a part in the disease progression of AD. First, expression changes of chemokines have been implicated in the pathogenesis of AD. Some chemokines were reported to be significantly altered in brain tissues , cerebral microcirculation , cerebrospinal fluid (CSF) [25, 26], plasma  or peripheral blood cells [28, 29] of AD patients, and in a mouse model of AD . Second, chemokine/chemokine receptor axes may affect AD pathologies in multiple ways. The impaired CCL2/CCR2 axis in blood-derived monocytes caused a deficit in cell migration in mild cognitive impairment (MCI) and AD patients . Dysregulation of the CX3CL1/CX3CR1 axis has shown both neuroprotective and neurotoxic effects in different AD mouse models [32, 33]. Third, a growing body of evidence indicates the chemokine/chemokine receptor axes affect Aβ and tau pathologies [32, 34,35,36], glymphatic function, cognition, or neuroinflammation [21, 36] in AD mouse models. Despite the functional involvement, several genetic association studies focusing on a few chemokine genes showed no direct association between common genetic polymorphisms of chemokine genes and the genetic susceptibility to AD [13, 37,38,39]. However, one study has suggested that a haplotype of single nucleotide polymorphism (SNP) within a chemokine gene cluster may modify the age of onset of familial AD . Since genetic variants play a role in the regulation and activation of chemokine/chemokine receptor signaling in AD , a systematic evaluation of the genetic variants and regulation of chemokines in AD is warranted.
In this study, we evaluated the transcriptomic dysregulation of chemokines in publicly available datasets of AD patients and mouse models. A three-stage genetic study of 31 chemokine genes was conducted in a Han Chinese cohort of 1280 AD cases and 5044 cognitively normal control subjects, to evaluate the potential association of genetic variation of chemokines with AD. Mendelian randomization (MR) was further used to evaluate the causally associated chemokine gene(s) in AD (Fig. 1). We found an active involvement of CXCL16 and CCL5 dysregulation in the development of AD.
Materials and methods
We assigned 31 typical chemokine genes, including nine C–C motif chemokine ligands (CCLs) and seven receptors (CCRs), nine C-X-C motif chemokine ligands (CXCLs) and four receptors (CXCRs), C-X3-C motif chemokine receptor (CX3CR1) and X-C motif chemokine receptor (XCR1), as defined by the KEGG “Chemokine signaling pathway” (https://www.gsea-msigdb.org/gsea/msigdb/cards/KEGG_CHEMOKINE_SIGNALING_PATHWAY.html) and the Immport database (https://www.immport.org/shared/genelists) .
mRNA expression profiling of chemokine genes in AD patients
The mRNA expression levels of the chemokine genes were analyzed in four AD-relevant brain regions (the entorhinal cortex, hippocampus, temporal cortex and frontal cortex) . The microarray expression data of the four brain regions from 269 AD patients and 271 controls were retrieved from Gene Expression Omnibus (GEO: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo) and were integrated to generate a normalized expression profile, as described in our previous study (http://www.alzdata.org/) . In brief, for each original microarray data retrieved from the GEO database, we conducted data normalization, log2 transformation, probe filtration, and probe mapping to entrez gene IDs. Expression datasets for the same brain region were then combined and re-normalized to remove batch effects. The normalized expression data was used to detect if there were any differentially expressed chemokine genes between AD patients and controls. More details regarding to data processing of these reported datasets can be retrieved from our previous study ( and references therein).
The expression alterations of chemokines in peripheral blood from individuals with and without AD were explored in two large independent age-matched dementia case–control data sets . The first dataset (GSE63060) includes 329 individuals containing 104 healthy controls, 80 MCI, and 145 AD patients from AddneuroMed Cohort (batch 1, https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE63060). The second dataset (GSE63061) includes 388 individuals, among which 382 were explicitly defined (including 134 healthy controls, 109 MCI, and 139 AD patients) (batch 2, https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE63061) . Patients with MCI and AD had different levels of cognitive impairments. Specifically, the MCI patients had problems with memory, but had normal daily activities. The healthy control subjects had no cognitive impairment. More description of these patients with AD, MCI, and healthy controls can be found in the original study . The expression profiling of the original microarray expression matrix series was loaded using the R package GEOquery. The limma R package was used to process the data and differential expression analysis was conducted under the linear model. More information regarding the related datasets can be found in the original publication .
mRNA expression profiling and pathological correlation in AD mouse models
So as to investigate the dynamic alteration of chemokines before and during the development of AD pathology, we retrieved the spatial–temporal expression data of chemokine genes in AD mouse models from Mouseac (www.mouseac.org) . We retrieved pathological scores of Aβ and tau of the AD mouse models and performed the correlation analysis between the mRNA expression level of chemokine genes and scored AD pathology by using the nonparametric Pearson correlation test with the GraphPad Prism software (GraphPad Software, La Jolla, CA, USA), as described in our previous study . The original study  measured the gene expression changes in the brain tissues (including 113 hippocampus samples, 113 cortex samples, and 111 cerebellum samples) from different AD murine models at the age of 2, 4, 8, and 18 months. Also, the levels of amyloid burden and phosphorylated tau pathology were investigated through immunostaining with antibodies for Aβ40 or phosphorylated tau, and the semi-quantitative scores were based on the pathology severity . In total, the expression data and pathological features of 114 brain tissues from wild-type mice (WILD) and 219 brain tissues from five AD transgenic mice with human APP, PSEN1 or MAPT mutant (TAS10 [with APP K670N/M671L mutant], TPM [with PSEN1 M146V mutant], HO_TASTPM [with homozygous mutant of APP and PSEN1 mentioned above], HET_TASTPM [with heterozygous mutant of APP and PSEN1 mentioned above], and TAU [with MAPT P301L mutant]) were analyzed. More information about these mouse models can be found in the original research .
Targeted sequencing of chemokine genes in Han Chinese
To investigate whether there is a potential genetic association between the chemokine gene variants and AD, we conducted targeted sequencing of 31 chemokine genes in Han Chinese with and without AD. The majority of Han Chinese samples have been described in our previous studies [46,47,48,49]. Briefly, two cohorts containing AD cases (n = 1280) and controls (n = 5044) were enrolled from Southern and Eastern China. All AD patients were confirmed to have no known pathogenic variants in APP, PSEN1, or PSEN2. The Southern cohort contains 635 sporadic AD patients (mean age 79.7 ± 8.2 years, 40.0% male) and 1507 controls (mean age 35.2 ± 15.5 years, 56.2% male) recruited from Sichuan, Hunan, and Yunnan Province of China. The Eastern cohort contains 645 AD patients (mean age 79.2 ± 9.1 years, 41.2% male) recruited from Shanghai and Zhejiang and 3537 controls from the general population with WGS data from the China Metabolic Analytics Project (ChinaMAP) . The patients were diagnosed by at least two clinical psychiatrists following the Diagnostic and Statistical Manual of Mental Disorders (DSM-IV), as have been previously described [46,47,48]. Some of the control samples from the Yunnan Province of China have been reported in our previous study [51, 52]. The Southern and Eastern cohorts were used in stage 1 and stage 2 analyses, respectively. Written informed consents conforming to the tenets of the Declaration of Helsinki were obtained from all participants before the enrollment of this study. The experimental protocols were approved by the institutional review board of Kunming Institute of Zoology, Chinese Academy of Sciences.
Genomic DNA was extracted from peripheral blood by using the AxyPrep Blood Genomic DNA Miniprep Kit (Axygen Scientific). DNA probes were designed using the online NimbleDesign tool (now updated as https://www.hyperdesign.com/#/). Coding regions, exon–intron boundaries, the 5’- and 3’-untranslated regions (UTRs) of the 31 genes were captured with the NimbleGen SeqCap EZ Choice Enrichment Kit (Roche NimbleGen) according to the manufacturer’s protocols, as described in our recent studies [51, 52]. Briefly, all captured DNA libraries were sequenced with the NovaSeq 6000 (150 bp paired-end). Raw reads were trimmed to remove sequencing adapters and low-quality reads by using the Trimmomatic (v0.33) . Clean reads were aligned according to the human reference genome GRCh37/hg19 with the Burrows-Wheeler Aligner . Post alignment quality control and variant calling were performed using the Genome Analysis Toolkit (GATK v4.1) following the best practices pipeline (https://www.broadinstitute.org/gatk/guide/best-practices) . Quality control was performed based on all of the sequencing data. Variants were excluded if the genotype rate were less than 90%, or deviated from Hardy–Weinberg equilibrium (HWE p < 1 × 10−6). Samples with an average genotype rate < 80% were also excluded.
Variants were annotated into different functional categories by using ANNOVAR . A variant was defined as rare if it had a minor allele frequency (MAF) < 0.01, otherwise, it was defined as common. Rare variants were further classified into three categories: loss-of-function (LoF) variants, missense variants, and possibly pathogenic variants. Variants belonging to stop gain/loss, frameshift indels, initiation codon, and splice sites were defined as LoF variants. Rare missense variants with sensitivity score ≤ 0.95 were defined as possibly pathogenic by using the Mendelian Clinically Applicable Pathogenicity (M-CAP) .
Cross-genetic validation in European-ancestry populations
Data from two previous studies [6, 44] were used for cross-validation of the association of the single variants and genes with AD identified in this study. The summary data of the newly published GWAS meta-analysis  contains a total of 111,326 clinically diagnosed/ “proxy” AD cases and 677,663 controls. We used this dataset to validate the potential association of common variants with AD in Han Chinese. The WES data from the discovery case–control association results of the Alzheimer’s Disease Sequencing Project (ADSP)  was used to validate the gene-based burden test discerned in this study. The cumulative minor allele counts from the 5740 AD cases and 5096 cognitively normal controls and the p-values adjusted for different covariates were extracted from the original study. More information regarding patient description and data quality control can be found in the original reports [6, 44].
In order to explore whether there was a causal effect of the observed association between chemokine and AD, we conducted MR analyses by using the R package “TwoSampleMR” . The causal effect of chemokine on AD was measured using chemokine as the exposure and AD as the outcome. Genetic variants (SNPs) associated with chemokine expression (p-value < 1 × 10−5) were set as the instrumental variables and were investigated in GWAS of outcomes. The causal effect of chemokine-related instruments (significant SNPs) on AD outcomes was assessed in large-scale GWAS for AD. In brief, we extracted plasma protein data for chemokines from 3301 healthy blood donors  and 3394 individuals with multiple cardiovascular diseases . GWAS data for AD were extracted from the stage 1 summary results of 21,982 diagnosed AD cases and 41,944 controls in a large GWAS meta-analysis study . Data were retrieved from IEU GWAS database (https://gwas.mrcieu.ac.uk/)  using the gwasglue R package (https://mrcieu.github.io/gwasglue/).
In consideration of the fact that reverse causation is a common confounding factor in observational studies, we also checked the reverse causal effect setting AD as the exposure and chemokine as the outcome in the MR analyses. SNPs associated with AD (p-value < 1 × 10−5) were used as the instrumental variables. The effect of the instruments was estimated in the corresponding GWAS for each outcome with inverse-variance weighted (IVW) linear regression . For the significance of the MR effect, a MR p-value < 0.05 was defined as significant. All MR analyses were conducted in R version 3.6.3 with packages.
We used Quanto software to evaluate the statistical power of our samples under the gene-only hypothesis and log additive model . Fisher’s exact test was applied to test if the allele frequency of a variant was significantly different between AD patients and controls. Meta-analyses combined the Southern and Eastern cohorts were performed using the metafor R package under the fixed effects model . Nominal significance was defined as p < 0.05. The Bonferroni correction for multiple testing was performed based on the corresponding numbers of tested variants or genes. Rare variants within each gene were arranged into a gene-based burden test. Single rare and common variants were subjected to association analysis. The Bonferroni corrected significance requires a p < 0.0016 (0.05/31 genes) for gene-based analysis, p < 5.96 × 10−5 (0.05/839 rare variants identified in this study) for rare variant association analysis, and p < 7.04 × 10−4 (0.05/71 common variants identified in this study) for common variant association analysis. The burden of rare variants in each targeted gene was tested by using the optimized sequence kernel association test (SKAT-O) in the R package SKAT . Rare variants classified into LoF, possibly pathogenic, or missense in each gene, were assessed. As we had no detailed genotype data of each individual from the ChinaMAP , we followed the same strategy as in our recent study  and performed the burden tests using only the 1507 healthy individuals from the Southern cohort as the control sample in the burden test.
Two-tailed Student’s t-test was used to investigate the mRNA expression difference between AD patients and controls with the GraphPad Prism software. Correlations between mRNA levels of the chemokine genes and AD pathology scores from Mouseac  were measured by using the nonparametric Pearson correlation test with the GraphPad Prism software.
Upregulation of mRNA levels of chemokine genes in AD brain tissues
Since expression changes of chemokines are implicated in the pathogenesis of AD , we compared the mRNA expression levels of chemokine genes in four brain tissues (entorhinal cortex, hippocampus, temporal cortex, and frontal cortex) between AD patients and controls, using the normalized microarray data compiled in our previous study . Among these 31 chemokine genes, no expression data of CCL3, CCL14, and CCL15 were available for the analysis. The mRNA expression levels of 15 chemokine genes were nominally upregulated in one or more of the four AD-relevant tissues. Four genes, CCL5, CXCL1, CXCL16, and CXCR4, survived multiple-testing and showed significant differential expression (Table 1 and Fig. 2), suggesting a robust dysregulation of these genes in AD. No significant alteration of mRNA expression was observed for the remaining 13 genes.
Correlation of chemokine mRNA expression levels with Aβ and tau pathology in AD mouse models
We reanalyzed the expression data of the AD mouse models and compared the gene expression levels of chemokines in brain tissues including 114 from WILD mice, 44 from HO_TASTPM mice, and 45 from TAU mice at different ages . The mRNA expression levels of 10 chemokine genes (Ccl3, Ccl5, Ccl8, Ccrl2, Cx3cr1, Cxcl1, Cxcl9, Cxcl10, Cxcl13, and Cxcl16) were upregulated along with age and reached a significant level of differential expression at the late stage of AD. Among these 10 chemokines, six genes (Ccl3, Ccl5, Cxcl1, Cxcl10, Cxcl13, and Cxcl16) were upregulated in both HO_TASTPM and TAU mice (Additional file 1: Fig. S1). We observed a strong positive correlation between the mRNA expression of the above 10 genes and AD pathology. Except for the Ccrl2 gene, whose mRNA expression was positively correlated only with Aβ pathology, the other nine genes showed a positive correlation of mRNA expression with both Aβ and tau pathology (Additional file 1: Fig. S2). These observations support an active involvement of chemokines in AD pathology in mouse models.
We overlapped the genes with significant expression alteration in brain tissues from AD patients and the mouse models, as we speculated that these genes validated in two systems may be more reliable for being causal for AD development. Three genes, including CCL5, CXCL1, and CXCL16 showed a consistent alteration pattern in AD patients and in the mouse models (Figs. 2 and 3).
Upregulation of CXCL16 expression in peripheral blood of patients with MCI and AD
To investigate if any of these dysregulated chemokine genes could be used as a potential biomarker for AD development, we analyzed the mRNA expression alterations of these genes in the peripheral blood from MCI patients, AD patients, and healthy controls. Among the three genes (CCL5, CXCL1, and CXCL16) consistently upregulated in brain tissues of AD patients and the mouse models, only CXCL16 showed an upregulated expression in peripheral blood of MCI and AD patients as compared to controls (Table 2 and Fig. 4). This result suggests the potential of CXCL16 to be a biomarker for AD development.
No association of chemokine genetic variants with AD in Han Chinese
Among these Han Chinese subjected to targeted sequencing, the mean sequencing depth of each gene was higher than 90 × in targeted sequencing (Additional file 1: Table S1), and the sequencing depth for each variant was higher than 25 × (Additional files 2 and 3). A total of 910 genetic variants (including 839 rare and 71 common variants) in 31 chemokine genes passed the quality control and were analyzed subsequently. The SKAT-O analysis  was conducted using rare variants in each gene for each cohort. In the Southern cohort, possibly pathogenic variants in CCR7 (p = 0.004) and missense variants in CCR9 (p = 0.009) were enriched in AD patients. In the Eastern cohort, CCL3, CCR5, and CCR9 showed an association with AD at the gene-burden level. When combined both cohorts together to achieve a better statistic power, we found that five genes (including CCL3, CCR7, CCR9, CCRL2, and XCR1) were associated with AD at the gene-burden level and had an enrichment of rare variants (Additional file 1: Table S2). However, none of these significances survived the multiple testing correction, partially due to the limited sample size.
To confirm the gene-based association between chemokine genes and AD identified in Han Chinese, we used the WES data in the discovery stage of ADSP (5740 AD cases and 5096 cognitively normal controls)  as a validation cohort. None of these chemokine genes associated with AD in our Han Chinese population could be validated in the ADSP dataset , although we found that four genes (CCL2, CCR6, CXCL6, CX3CR1) were associated with AD at the gene-burden level under correction for different covariates, including the principal components, sequencing center, sex, age, and APOE ε4 and ε2 dosages, with the SKAT-O analysis  by using the ADSP dataset  (Additional file 1: Table S3).
Next, we tested the single-variant association of 839 rare and 71 common variants with AD risk in each Han Chinese cohort and the combined sample. We observed nominally significant association with AD of eight rare variants in six of the 31 chemokine genes in meta-analysis, but none of these variants survived the multiple testing correction (Additional file 2). The MAF of 71 common variants ranged from 0.01 to 0.5 in the combined control samples. Under the gene-only hypothesis and log additive model with an average population MAF of 0.1, the statistical power to detect an odds ratio (OR) value of 1.25 for a risk allele using the current sample size (1280 cases and 5044 controls) was above 88.8%. We observed a significant association of rs181868085 (p = 0.006, OR = 1.59) in CXCL1 and rs2304973 (p = 0.045, OR = 0.84) in CXCL16 with AD in the meta-analysis (Additional file 3). The OR values of these two SNPs indicated a consistent direction of genetic effect on disease risk between the Southern and the Eastern cohorts, although the association did not survive the multiple testing correction. Note that we observed a significant association between CCL3 and AD in the Eastern cohort of Han Chinese, even after the Bonferroni correction (p = 4.35 × 10−4), but this association was weakened in the combined Han Chinese sample (p = 0.008; Additional file 1: Table S2), suggesting that population heterogeneity may exist even between our two cohorts under study. A reanalysis of the allele frequency data of 71 common variants in the newly published GWAS study  showed no association with AD in the European population (Additional file 1: Table S4).
Taken together, there seems to be no robust association of rare and common variants of chemokine genes with AD in both Han Chinese and European populations. The weak association between the chemokine gene and AD observed in the gene-based burden test might be caused by population substructure and stratification.
MR analyses prioritized CCL5 as a causal gene for AD
The MR analysis based on large-scale proteomic data of plasma chemokine levels [59, 60] and a genetic study of AD  showed that CCL5 (p = 0.0055, β = − 0.0667, Table 3) was causally linked to AD without the reverse causal effect (p = 0.7864, β = 0.0078). This result, together with the observation of altered expression of CCL5 in AD patients and the mouse models, indicated that CCL5 dysregulations might be actively involved in the development of AD.
Alteration of chemokines has been frequently reported to be involved in AD [28, 31, 32]. However, the exact mechanisms of the upregulation or downregulation of chemokine genes during the development of AD have not been sufficiently clarified, and this has led to a dispute regarding whether the chemokine expression alterations are the drivers, or by-products, of AD pathobiology. In this study, we aimed to determine the role of mRNA expression and genetic variations of 31 chemokine genes in AD, with an intention of defining the involvement of chemokines in AD. By integrating the genetic analyses, mRNA expression alterations, and pathological correlation in both AD patients and mouse models, we found an involvement of CXCL16 and CCL5 in the development of AD. This comprehensive analysis enabled us to provide a systematic view for understanding the roles of chemokines in the development of AD. First, the mRNA expression levels of a small proportion of chemokines under study were upregulated in brain regions both in AD patients (Table 1 and Fig. 2) and AD mice (Fig. 3). This result is consistent with previous reports for a higher level of proinflammatory chemokines in peripheral blood and brain tissues in AD patients [28, 64], and supports the significant role of chemokines in the neuroinflammatory process of AD. Second, we found that the upregulated chemokines were positively correlated with Aβ and tau pathology in AD mice (Fig. 3). This observation suggested an active role of chemokines in AD progression although their function in regulating Aβ or tau pathology remains to be determined . Third, the MR analysis showed that CCL5 was prioritized to be causally linked to AD, indicating that chemokine gene may be involved in AD in different ways.
We did not obtain a firm conclusion regarding the association of genetic variants of 31 chemokine genes with AD. Although we found five genes (CCL3, CCR7, CCR9, CCRL2, and XCR1) were nominally associated with AD in the combined Han Chinese sample in the gene-based burden test (Additional file 1: Table S2), none of the associations can be verified in the populations of European ancestry based on a data-mining of the ADSP dataset  (Additional file 1: Table S3). Instead, four other genes, CCL2, CCR6, CX3CR1, and CXCL6, were suggested to be associated with AD in the populations of European ancestry . Our analysis of the potential association between rare or common variants and AD risk showed the same pattern of mixed signal for positive association. Therefore, we were unable to make a firm conclusion that genetic variants of the chemokine genes had a role in conferring genetic risk to AD, and the current weak association could be real or could be explained by different genetic structures between Asian and European populations. Independent validation analysis with a large sample size in populations of different ancestral origins is needed to clarify this issue.
The dysregulation of CCL5 and CXCL16 in AD patients at the transcriptional level deserved further attention. The CCL5 was reported to be upregulated in peripheral blood mononuclear cells of AD patients [29, 65], and Aβ42 treatment could increase the expression of CCL5 and its receptor CCR5 in peripheral mononuclear cells . The upregulation of CCL5 in the AD brain may play a possible neuroprotective role , as soluble CCL5 activated by Aβ had an ameliorating effect on AD in mice by recruiting bone marrow-induced microglia immune response . In this study, we found that CCL5 increases with AD development in both mouse models and patients, and an increased CCL5 level was causally associated with decreased AD risk in the MR analysis, collectively supporting the active role of CCL5 in AD and the potential utility of CCL5 as a therapeutic target.
Concerning the involvement of CXCL16 in AD, there were only a few reports available until very recently, Piehl et al. highlighted the CXCL16-CXCR6 axis in CSF of aged and AD brain . These researchers found that CXCL16, derived from inflamed microglia and increased in CSF, activated the CD8+ T cell trafficking to the CSF through the CXCL16-CXCR6 pathway . Note that dysregulation of CXCL16 was previously suggested as a possible mechanism of neurodegeneration in AD . The consistent alteration pattern of CXCL16 in serum, CSF, and brain tissues of AD patients suggested that this chemokine might be used as a potential biomarker for monitoring AD development. A clinical observation study is needed to test this possibility.
This study has some limitations. First, the sample size used in genetic analyses was relatively small, and potential population stratification and different population structure might blur the potential association between chemokine genes and AD. Second, although we observed expression alterations of some chemokine genes in brain tissues of AD patients and the mouse models, we have not linked the expression changes to cell types and had no experimental data to discern its potential effect on cellular function of these affected cells. We also did not validate the causal role of the highlighted genes, such as CCL5 and CXCL16, in animal models of AD . Third, we only analyzed a proportion of chemokines in this study, and we could not exclude the possibility that other chemokine genes might have a prominent role in AD pathobiology. Nonetheless, the accumulating knowledge of chemokines’ roles in AD, as exemplified in this study with an intention of comprehensive integrative analysis, is undoubtedly essential for guiding the development of potential novel immunotherapies for AD.
In short, through an extensive analysis of chemokine genes based on expressional, pathological, and genetic analysis data, we provide multiple lines of evidence to support the important role of chemokines CCL5 and CXCL16 in the development of AD. Further genetic studies with larger sample sizes and functional assays are needed to validate our conclusion and to depict the mechanisms of these two chemokines in AD pathogenesis.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
Mild cognitive impairment
Genome-wide association study
Whole-exome sequencing study
Whole-genome sequencing study
Central nervous system
Single nucleotide polymorphism
Homozygous mutant human APP K670N/M671L and PSEN1 M146V mice
Mutant human MAPT P301L mice
China Metabolic Analytics Project
Diagnostic and Statistical Manual of Mental Disorders IV
Minor allele frequency
Mendelian Clinically Applicable Pathogenicity
Inverse-variance weighted linear regression
The Alzheimer’s Disease Sequencing Project
The optimized sequence kernel association test
Querfurth HW, LaFerla FM. Alzheimer’s disease. N Engl J Med. 2010;362:329–44.
Scheltens P, De Strooper B, Kivipelto M, Holstege H, Chetelat G, Teunissen CE, et al. Alzheimer’s disease. Lancet. 2021;397:1577–90.
Kunkle BW, Grenier-Boley B, Sims R, Bis JC, Damotte V, Naj AC, et al. Genetic meta-analysis of diagnosed Alzheimer’s disease identifies new risk loci and implicates Abeta, tau, immunity and lipid processing. Nat Genet. 2019;51:414–30.
Jansen IE, Savage JE, Watanabe K, Bryois J, Williams DM, Steinberg S, et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet. 2019;51:404–13.
Wightman DP, Jansen IE, Savage JE, Shadrin AA, Bahrami S, Holland D, et al. A genome-wide association study with 1,126,563 individuals identifies new risk loci for Alzheimer’s disease. Nat Genet. 2021;53:1276–82.
Bellenguez C, Küçükali F, Jansen IE, Kleineidam L, Moreno-Grau S, Amin N, et al. New insights into the genetic etiology of Alzheimer’s disease and related dementias. Nat Genet. 2022;54:412–36.
Jia L, Li F, Wei C, Zhu M, Qu Q, Qin W, et al. Prediction of Alzheimer’s disease using multi-variants from a Chinese genome-wide association study. Brain. 2021;144:924–37.
Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, Majounie E, et al. TREM2 variants in Alzheimer’s disease. N Engl J Med. 2013;368:117–27.
Zhang DF, Fan Y, Xu M, Wang G, Wang D, Li J, et al. Complement C7 is a novel risk gene for Alzheimer’s disease in Han Chinese. Natl Sci Rev. 2019;6:257–74.
Zhou X, Chen Y, Mok KY, Zhao Q, Chen K, Chen Y, et al. Identification of genetic risk factors in the Chinese population implicates a role of immune system in Alzheimer’s disease pathogenesis. Proc Natl Acad Sci U S A. 2018;115:1697–706.
Heneka MT, Carson MJ, El Khoury J, Landreth GE, Brosseron F, Feinstein DL, et al. Neuroinflammation in Alzheimer’s disease. Lancet Neurol. 2015;14:388–405.
Koelink PJ, Overbeek SA, Braber S, de Kruijf P, Folkerts G, Smit MJ, et al. Targeting chemokine receptors in chronic inflammatory diseases: an extensive review. Pharmacol Ther. 2012;133:1–18.
Su F, Bai F, Zhang Z. Inflammatory cytokines and Alzheimer’s disease: a review from the perspective of genetic polymorphisms. Neurosci Bull. 2016;32:469–80.
Savarin-Vuaillat C, Ransohoff RM. Chemokines and chemokine receptors in neurological disease: raise, retain, or reduce? Neurotherapeutics. 2007;4:590–601.
Bajetto A, Bonavia R, Barbero S, Florio T, Schettini G. Chemokines and their receptors in the central nervous system. Front Neuroendocrinol. 2001;22:147–84.
Réaux-Le Goazigo A, Van Steenwinckel J, Rostène W, Mélik PS. Current status of chemokines in the adult CNS. Prog Neurobiol. 2013;104:67–92.
Dziembowska M, Tham TN, Lau P, Vitry S, Lazarini F, Dubois-Dalcq M. A role for CXCR4 signaling in survival and migration of neural and oligodendrocyte precursors. Glia. 2005;50:258–69.
Tsai HH, Frost E, To V, Robinson S, Ffrench-Constant C, Geertman R, et al. The chemokine receptor CXCR2 controls positioning of oligodendrocyte precursors in developing spinal cord by arresting their migration. Cell. 2002;110:373–83.
Padovani-Claudio DA, Liu L, Ransohoff RM, Miller RH. Alterations in the oligodendrocyte lineage, myelin, and white matter in adult mice lacking the chemokine receptor CXCR2. Glia. 2006;54:471–83.
Bertollini C, Ragozzino D, Gross C, Limatola C, Eusebi F. Fractalkine/CX3CL1 depresses central synaptic transmission in mouse hippocampal slices. Neuropharmacology. 2006;51:816–21.
Da Mesquita S, Herz J, Wall M, Dykstra T, de Lima KA, Norris GT, et al. Aging-associated deficit in CCR7 is linked to worsened glymphatic function, cognition, neuroinflammation, and beta-amyloid pathology. Sci Adv. 2021;7:eabe4601.
Zhou L, Kong G, Palmisano I, Cencioni MT, Danzi M, De Virgiliis F, et al. Reversible CD8 T cell-neuron cross-talk causes aging-dependent neuronal regenerative decline. Science. 2022;376:eabd5926.
Sokolova A, Hill MD, Rahimi F, Warden LA, Halliday GM, Shepherd CE. Monocyte chemoattractant protein-1 plays a dominant role in the chronic inflammation observed in Alzheimer’s disease. Brain Pathol. 2009;19:392–8.
Tripathy D, Thirumangalakudi L, Grammas P. RANTES upregulation in the Alzheimer’s disease brain: a possible neuroprotective role. Neurobiol Aging. 2010;31:8–16.
Kauwe JS, Bailey MH, Ridge PG, Perry R, Wadsworth ME, Hoyt KL, et al. Genome-wide association study of CSF levels of 59 alzheimer’s disease candidate proteins: significant associations with proteins involved in amyloid processing and inflammation. PLoS Genet. 2014;10:e1004758.
Perea JR, Lleo A, Alcolea D, Fortea J, Avila J, Bolos M. Decreased CX3CL1 levels in the cerebrospinal fluid of patients with Alzheimer’s disease. Front Neurosci. 2018;12:609.
Lee WJ, Liao YC, Wang YF, Lin IF, Wang SJ, Fuh JL. Plasma MCP-1 and cognitive decline in patients with Alzheimer’s disease and mild cognitive impairment: a two-year follow-up study. Sci Rep. 2018;8:1280.
Goldeck D, Larbi A, Pellicano M, Alam I, Zerr I, Schmidt C, et al. Enhanced chemokine receptor expression on leukocytes of patients with Alzheimer’s disease. PLoS ONE. 2013;8:e66664.
Iarlori C, Gambi D, Gambi F, Lucci I, Feliciani C, Salvatore M, et al. Expression and production of two selected beta-chemokines in peripheral blood mononuclear cells from patients with Alzheimer’s disease. Exp Gerontol. 2005;40:605–11.
Jorda A, Cauli O, Santonja JM, Aldasoro M, Aldasoro C, Obrador E, et al. Changes in chemokines and chemokine receptors expression in a mouse model of Alzheimer’s disease. Int J Biol Sci. 2019;15:453–63.
Guedes JR, Santana I, Cunha C, Duro D, Almeida MR, Cardoso AM, et al. MicroRNA deregulation and chemotaxis and phagocytosis impairment in Alzheimer’s disease. Alzheimers Dement (Amst). 2016;3:7–17.
Guedes JR, Lao T, Cardoso AL, El Khoury J. Roles of microglial and monocyte chemokines and their receptors in regulating Alzheimer’s disease-associated amyloid-beta and tau pathologies. Front Neurol. 2018;9:549.
Subbarayan MS, Joly-Amado A, Bickford PC, Nash KR. CX3CL1/CX3CR1 signaling targets for the treatment of neurodegenerative diseases. Pharmacol Ther. 2022;231:107989.
Bhaskar K, Konerth M, Kokiko-Cochran ON, Cardona A, Ransohoff RM, Lamb BT. Regulation of tau pathology by the microglial fractalkine receptor. Neuron. 2010;68:19–31.
Liu Z, Condello C, Schain A, Harb R, Grutzendler J. CX3CR1 in microglia regulates brain amyloid deposition through selective protofibrillar amyloid-beta phagocytosis. J Neurosci. 2010;30:17091–101.
Kiyota T, Gendelman HE, Weir RA, Higgins EE, Zhang G, Jain M. CCL2 affects beta-amyloidosis and progressive neurocognitive dysfunction in a mouse model of Alzheimer’s disease. Neurobiol Aging. 2013;34:1060–8.
Vendramini AA, de Labio RW, Rasmussen LT, Minett T, Bertolucci PH, de Arruda Cardoso Smith M, et al. Interleukin-8 gene polymorphism -251T>A and Alzheimer’s disease. J Alzheimers Dis. 2007;12:221–2.
Villa C, Venturelli E, Fenoglio C, Clerici F, Marcone A, Benussi L, et al. CCL8/MCP-2 association analysis in patients with Alzheimer’s disease and frontotemporal lobar degeneration. J Neurol. 2009;256:1379–81.
Balistreri CR, Grimaldi MP, Vasto S, Listi F, Chiappelli M, Licastro F, et al. Association between the polymorphism of CCR5 and Alzheimer’s disease: results of a study performed on male and female patients from Northern Italy. Ann N Y Acad Sci. 2006;1089:454–61.
Lalli MA, Bettcher BM, Arcila ML, Garcia G, Guzman C, Madrigal L, et al. Whole-genome sequencing suggests a chemokine gene cluster that modifies age at onset in familial Alzheimer’s disease. Mol Psychiatry. 2015;20:1294–300.
Xu M, Zhang DF, Luo R, Wu Y, Zhou H, Kong LL, et al. A systematic integrated analysis of brain expression profiles reveals YAP1 and other prioritized hub genes as important upstream regulators in Alzheimer’s disease. Alzheimers Dement. 2018;14:215–29.
Sood S, Gallagher IJ, Lunnon K, Rullman E, Keohane A, Crossland H, et al. A novel multi-tissue RNA diagnostic of healthy ageing relates to cognitive health status. Genome Biol. 2015;16:185.
Matarin M, Salih DA, Yasvoina M, Cummings DM, Guelfi S, Liu W, et al. A genome-wide gene-expression analysis and database in transgenic mice during development of amyloid or tau pathology. Cell Rep. 2015;10:633–44.
Bis JC, Jian X, Kunkle BW, Chen Y, Hamilton-Nelson KL, Bush WS, et al. Whole exome sequencing study identifies novel rare and common Alzheimer’s-associated variants involved in immune response and transcriptional regulation. Mol Psychiatry. 2020;25:1859–75.
Bhattacharya S, Dunn P, Thomas CG, Smith B, Schaefer H, Chen J, et al. ImmPort, toward repurposing of open access immunological assay data for translational and clinical research. Sci Data. 2018;5:180015.
Bi R, Zhang W, Zhang DF, Xu M, Fan Y, Hu QX, et al. Genetic association of the cytochrome c oxidase-related genes with Alzheimer’s disease in Han Chinese. Neuropsychopharmacology. 2018;43:2264–76.
Zhang DF, Li J, Wu H, Cui Y, Bi R, Zhou HJ, et al. CFH variants affect structural and functional brain changes and genetic risk of Alzheimer’s disease. Neuropsychopharmacology. 2016;41:1034–45.
Li GD, Bi R, Zhang DF, Xu M, Luo R, Wang D, et al. Female-specific effect of the BDNF gene on Alzheimer’s disease. Neurobiol Aging. 2017;53(192):e11–9.
Wang HZ, Bi R, Hu QX, Xiang Q, Zhang C, Zhang DF, et al. Validating GWAS-identified risk loci for Alzheimer’s disease in Han Chinese populations. Mol Neurobiol. 2016;53:379–90.
Cao Y, Li L, Xu M, Feng Z, Sun X, Lu J, et al. The ChinaMAP analytics of deep whole genome sequences in 10,588 individuals. Cell Res. 2020;30:717–31.
Wang D, Fan Y, Malhi M, Bi R, Wu Y, Xu M, et al. Missense variants in HIF1A and LACC1 contribute to leprosy risk in Han Chinese. Am J Hum Genet. 2018;102:794–805.
Li X, Xu M, Bi R, Tan LW, Yao YG, Zhang DF. Common and rare variants of EGF increase the genetic risk of Alzheimer’s disease as revealed by targeted sequencing of growth factors in Han Chinese. Neurobiol Aging. 2022;S0197-4580(22)00223-8.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164.
Jagadeesh KA, Wenger AM, Berger MJ, Guturu H, Stenson PD, Cooper DN, et al. M-CAP eliminates a majority of variants of uncertain significance in clinical exomes at high sensitivity. Nat Genet. 2016;48:1581–6.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408.
Sun BB, Maranville JC, Peters JE, Stacey D, Staley JR, Blackshaw J, et al. Genomic atlas of the human plasma proteome. Nature. 2018;558:73–9.
Folkersen L, Fauman E, Sabater-Lleal M, Strawbridge RJ, Franberg M, Sennblad B, et al. Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease. PLoS Genet. 2017;13:e1006706.
Gauderman WJ. Sample size requirements for matched case-control studies of gene-environment interaction. Stat Med. 2002;21:35–50.
Viechtbauer W. Conducting meta-analyses in R with the metafor package. J Stat Softw. 2010;36:1–48.
Lee S, Emond MJ, Bamshad MJ, Barnes KC, Rieder MJ, Nickerson DA, et al. Optimal unified approach for rare-variant association testing with application to small-sample case-control whole-exome sequencing studies. Am J Hum Genet. 2012;91:224–37.
Zuena AR, Casolini P, Lattanzi R, Maftei D. Chemokines in Alzheimer’s disease: new insights into prokineticins, chemokine-like proteins. Front Pharmacol. 2019;10:622.
Pellicano M, Bulati M, Buffa S, Barbagallo M, Di Prima A, Misiano G, et al. Systemic immune responses in Alzheimer’s disease: in vitro mononuclear cell activation and cytokine production. J Alzheimers Dis. 2010;21:181–92.
Reale M, Iarlori C, Feliciani C, Gambi D. Peripheral chemokine receptors, their ligands, cytokines and Alzheimer’s disease. J Alzheimers Dis. 2008;14:147–59.
Lee JK, Schuchman EH, Jin HK, Bae JS. Soluble CCL5 derived from bone marrow-derived mesenchymal stem cells and activated by amyloid beta ameliorates Alzheimer’s disease in mice by recruiting bone marrow-induced microglia immune responses. Stem Cells. 2012;30:1544–55.
Piehl N, van Olst L, Ramakrishnan A, Teregulova V, Simonton B, Zhang Z, et al. Cerebrospinal fluid immune dysregulation during healthy brain aging and cognitive impairment. Cell. 2022. https://0-doi-org.brum.beds.ac.uk/10.1016/j.cell.2022.11.019.
Yeung SS, Ho YS, Chang RC. The role of meningeal populations of type II innate lymphoid cells in modulating neuroinflammation in neurodegenerative diseases. Exp Mol Med. 2021;53:1251–67.
Chen ZY, Zhang Y. Animal models of Alzheimer’s disease: Applications, evaluation, and perspectives. Zool Res. 2022;43:1026–40.
We thank Dr. Ian Logan for language editing and Prof. Yanan Cao for sharing allele frequency data in the ChinaMAP project.
This study was supported by the Ministry of Science and Technology of China (grant number STI2030-Major Projects 2022ZD0213500), the National Natural Science Foundation of China (grant numbers 31730037, 82022017, 31970965); and the Strategic Priority Research Program (B) of the Chinese Academy of Sciences (CAS) (grant number XDB02020003), and the Youth Innovation Promotion Association of CAS.
Ethics approval and consent to participate
This study was approved by the institutional review board of Kunming Institute of Zoology, Chinese Academy of Sciences. Written informed consents were obtained from all participants in targeted sequencing. The study was adhered to the tenets of the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Fig. S1. Up-regulation of mRNA expression of chemokine genes during aging in AD mouse models. Fig. S2. Correlation between the mRNA expression levels of chemokine genes with AD pathology in AD mouse models. Table S1. Sequence coverage of each gene in the targeted sequencing. Table S2. Results of SKAT-O analysis in Han Chinese with and without AD. Table S3. Results of gene-based burden test of chemokine genes from ADSP. Table S4. Association results of common variants in European population.
Additional file 2.
Rare variants of 31 chemokine genes in Han Chinese.
Additional file 3.
Association results of common variants in 31 chemokine genes in Han Chinese of AD.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, X., Zhang, DF., Bi, R. et al. Convergent transcriptomic and genomic evidence supporting a dysregulation of CXCL16 and CCL5 in Alzheimer’s disease. Alz Res Therapy 15, 17 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s13195-022-01159-5
- Alzheimer’s disease
- mRNA expression
- Targeted sequencing
- Mendelian randomization