ORIGINAL RESEARCH article

Front. Genet., 17 February 2023

Sec. Applied Genetic Epidemiology

Volume 14 - 2023 | https://doi.org/10.3389/fgene.2023.1124431

Circulatory proteins relate cardiovascular disease to cognitive performance: A mendelian randomisation study

  • 1. Singapore Institute for Clinical Sciences (SICS), Agency for Science, Technology and Research (A*STAR), Singapore, Singapore

  • 2. Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London, London, United Kingdom

  • 3. Department of Brain Sciences, Faculty of Medicine, Imperial College London, London, United Kingdom

  • 4. UK Dementia Research Institute at Imperial College London, London, United Kingdom

  • 5. MRC Centre for Environment and Health, School of Public Health, Imperial College London, London, United Kingdom

  • 6. Department of Hygiene and Epidemiology, University of Ioannina Medical School, Ioannina, Greece

Article metrics

View details

1

Citations

5,2k

Views

2,7k

Downloads

Abstract

Background and objectives: Mechanistic research suggests synergistic effects of cardiovascular disease (CVD) and dementia pathologies on cognitive decline. Interventions targeting proteins relevant to shared mechanisms underlying CVD and dementia could also be used for the prevention of cognitive impairment.

Methods: We applied Mendelian randomisation (MR) and colocalization analysis to investigate the causal relationships of 90 CVD-related proteins measured by the Olink CVD I panel with cognitive traits. Genetic instruments for circulatory protein concentrations were obtained using a meta-analysis of genome-wide association studies (GWAS) from the SCALLOP consortium (N = 17,747) based on three sets of criteria: 1) protein quantitative trait loci (pQTL); 2) cis-pQTL (pQTL within ±500 kb from the coding gene); and 3) brain-specific cis-expression QTL (cis-eQTL) which accounts for coding gene expression based on GTEx8. Genetic associations of cognitive performance were obtained from GWAS for either: 1) general cognitive function constructed using Principal Component Analysis (N = 300,486); or, 2) g Factor constructed using genomic structural equation modelling (N = 11,263–331,679). Findings for candidate causal proteins were replicated using a separate protein GWAS in Icelanders (N = 35,559).

Results: A higher concentration of genetically predicted circulatory myeloperoxidase (MPO) was nominally associated with better cognitive performance (p < 0.05) using different selection criteria for genetic instruments. Particularly, brain-specific cis-eQTL predicted MPO, which accounts for protein-coding gene expression in brain tissues, was associated with general cognitive function (βWald = 0.22, PWald = 2.4 × 10−4). The posterior probability for colocalization (PP.H4) of MPO pQTL with the g Factor was 0.577. Findings for MPO were replicated using the Icelandic GWAS. Although we did not find evidence for colocalization, we found that higher genetically predicted concentrations of cathepsin D and CD40 were associated with better cognitive performance and a higher genetically predicted concentration of CSF-1 was associated with poorer cognitive performance.

Conclusion: We conclude that these proteins are involved in shared pathways between CVD and those for cognitive reserve or affecting cognitive decline, suggesting therapeutic targets able to reduce genetic risks conferred by cardiovascular disease.

Introduction

The main protein biomarkers identified for dementia include amyloid precursor protein and tau protein (O'Brien and Wong, 2011). Clinical trials of anti-amyloid and anti-tau intervention for AD have been conducted but meaningfully effective drugs are still not available (Congdon and Sigurdsson, 2018; Huang et al., 2020), and their causal role is under question (Kametani and Hasegawa, 2018; Morris et al., 2018; Thomas et al., 2020; Sturchio et al., 2021). Accumulating evidence has suggested a link between cardiovascular disease (CVD) and dementia (Paciaroni and Bogousslavsky, 2013). Classical cardiovascular risk factors, such as smoking, diabetes, hypertension, obesity, and physical inactivity contribute to a higher risk of cognitive impairment and dementia (Juul Rasmussen et al., 2020; Yaffe et al., 2020). Mechanistic research suggests that the impact of CVD risk burden on dementia pathologies could lead to cognitive decline (Attems and Jellinger, 2014; Santos et al., 2017). Plasma proteins for CVD related to several biological pathways that may be responsible for this have been identified. Confident identification of causal relationships between these protein biomarkers and cognition would support the development of pathway-specific treatment (Ho et al., 2018; Wallentin et al., 2021). With strong evidence that such proteins play shared causal roles in the mechanisms underlying CVD and cognition, interventions targeting these proteins and their pathways could be developed to prevent, slow or reverse disease progression and cognitive impairment.

Some recent trials have shifted the focus to the early stages of AD or mild cognitive impairment (Congdon and Sigurdsson, 2018). Given the high failure rate of the translation of preclinical drug candidates in animal models into treatment in humans, a more cost-effective strategy for drug discovery for cognitive impairment is needed. The recent development of genotyping and proteomic technologies has enabled the detection of protein quantitative trait loci (pQTL) in large-scale genome-wide association studies (GWAS) (Folkersen et al., 2020). Albeit most GWAS in human subjects interrogate plasma pQTL (rather than brain or brain cell-specific pQTL), they still can be used as genetic instruments for instrumental variable analyses such as Mendelian randomisation (MR) (Gill et al., 2021). MR analysis uses genetic variants, which are randomised at conception, to mimic random allocation in clinical trials and investigate the causal relationship between a risk factor and an outcome (Harris et al., 2020). It is less susceptible to confounding or reverse causation than conventional observational studies and can be used to identify potential molecular targets for intervention (Harris et al., 2020). In this study, we applied MR analyses to assess the associations of a wide range of plasma proteins, which are known to associate with CVD, with cognitive traits.

Methods

Genetic associations of concentrations of plasma proteins

Genetic associations of 90 proteins in plasma were obtained from a GWAS meta-analysis comprised of 13 cohorts of European ancestry (SCALLOP consortium; average per-protein sample size, 17,747) (Folkersen et al., 2020). Relative protein quantification was measured using the Olink proximity extension assay (PEA) cardiovascular (CVD) I panel (Assarsson et al., 2014). The selected protein biomarkers are associated with cardiovascular risk or prognosis in human observational studies and animal models (Assarsson et al., 2014; Folkersen et al., 2020). The Olink PEA CVD I panel returns normalized protein expression (NPX) values (on the log2 scale so that each one-unit difference in NPX indicates a doubling of protein concentration). NPX values were rank-based inverse-normalised to unit variance for the genetic association analysis. Supplementary Table S1 listed the proteins under investigation in this study.

To replicate the findings, we obtained the summary statistics for the genetic instruments from a separate GWAS of protein concentrations measured by SomaScan assay among 35,559 Icelanders (Ferkingstad et al., 2021). The SomaScan assay provides an aptamer-based measurement of the relative binding of a putative target protein to each aptamer in relative fluorescence units (Ferkingstad et al., 2021). Because of the methodological differences, the effect sizes from analyses based on Olink and Somalogic data cannot be compared directly.

Genetic associations of cognition

We obtained the effect estimates for the association of genetic instruments with two cognitive traits constructed using different plasma protein datasets (Davies et al., 2018; de la Fuente et al., 2021). A GWAS meta-analysis for general cognitive function (N = 300,486, age 16–102 years) was performed based on the first unrotated principal component of multiple cognitive tasks from 57 population-based cohorts of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) and the Cognitive Genomics Consortium (COGENT) consortia and the verbal and numerical reasoning test in UK Biobank (Davies et al., 2018). A GWAS for the general factor of intelligence (g Factor) was performed based on seven different cognitive tests (N = 11,263–331,679, age 40–70  years at first assessment) in UK Biobank (de la Fuente et al., 2021). For this association analysis, a univariate GWAS was performed for each of the seven standardised cognitive test scores and the univariate summary statistics then were used to conduct a multivariate GWAS using genomic structural equation modelling (Genomic SEM) (Grotzinger et al., 2019). In Genomic SEM, the target trait represents the genetic components of the individual GWAS traits (de la Fuente et al., 2021). A higher g Factor score indicates better performance in cognitive tasks. Supplementary Table S2 shows the details of GWAS of cognition.

Bidirectional two-sample univariable mendelian randomisation

We performed a two-sample univariable MR to investigate the associations of each protein with either general cognitive function or the g Factor. For each MR analysis, we selected genetic instruments for each protein using three sets of criteria. Under the first set of criteria, we selected any single-nucleotide polymorphisms (SNPs) throughout the genome associated with the protein of interest at a p-value<5 × 10−8 (both cis- and trans-). For the proteins with less than three SNPs under these criteria, we loosened the p-value threshold to 5 × 10−6. We referred to this set of instruments as pQTL instrumental variables. Under the second set of criteria, we selected cis-pQTL (i.e., pQTL within ±500 kb from the protein coding gene) associated with the protein of interest at p-value<5 × 10−6. To test directly for potentially causal associations with brain protein expression, we, under the third set of criteria, leveraged the information from GTEx8 expression quantitative trait loci (eQTL) by meta-analysing eQTL from 13 brain tissues (Supplementary Table S3; sample sizes ranged from 139 to 255). We selected SNPs located within ±500 kb from the coding gene, associated with expression of the coding gene at p-value<10−4, and associated with circulating protein concentration at p-value<0.05. We defined these as brain-specific cis-eQTL. For all three sets of criteria, we only included SNPs with a minor allele frequency (MAF) greater than 5% and F-statistics greater than 10. Correlated SNPs (r2 ≥ 0.001) were excluded by keeping the one with the smallest p-value.

We also performed a two-sample MR to investigate whether cognitive traits or their genetic liability affect the plasma concentrations of selected marker proteins. We selected independent SNPs (r2 < 0.001) associated with the cognitive traits of interest at p-value<5 × 10−8, MAF greater than 5%, and F-statistics>10 as genetic instruments.

For each MR analysis, genetic associations of the selected genetic instruments with the outcome of interest were obtained from the corresponding GWAS. We estimated the SNP-specific effects using the Wald ratio (Burgess et al., 2017) and pooled the SNP-specific estimates using inverse-variance weighted (IVW) for the MR effect estimates (Lawlor et al., 2008). We used IVW fixed-effects model for the analyses with two or three SNPs and IVW random-effects model for the analyses with more than three SNPs (Bowden et al., 2017). We also performed two sensitivity methods (weighted median (WM) and MR-Egger regression) to assess the robustness of MR estimates and horizontal pleiotropic effects for analyses with at least three SNPs (Bowden et al., 2016; Burgess and Thompson, 2017). Potential outlier SNPs were identified using MR-PRESSO and excluded from the analysis (Verbanck et al., 2018). We accounted for multiple comparisons using Bonferroni correction for 90 plasma proteins and two directions with a p-value threshold of 0.05/180 = 0.0003. We prioritised candidate proteins by accounting for the findings based on all three sets of criteria (i.e., pQTL, cis-pQTL, and brain-specific cis-eQTL)

Colocalization

We performed colocalization to estimate the posterior probability of the hypothesis that a cognitive trait shares the same causal variant with a candidate protein (i.e., PP.H4) using both pQTL and brain-specific eQTL (Zuber et al., 2022). Specifically, we focused on the genomic region within ±50 kb from the protein coding gene of interest. Given that we performed colocalization for the candidate MR findings, here we tested only whether the common causal variant is more likely than other hypotheses (see Supplementary Table S8). Evidence of colocalization provides complementary information on causal relationships since distinct causal variants for the exposure and outcome of interest could lead to a violation of the exchangeability assumption in MR analysis (Zuber et al., 2022). We considered a PP.H4>0.5 as evidence for colocalization (Giambartolomei et al., 2014).

Replication analysis

To validate our findings, we repeated the analysis of findings in our main MR analyses using summary statistics from a GWAS of putative protein concentrations measured using the SomaScan assay (Ferkingstad et al., 2021). Specifically, we conducted MR and colocalisation analysis for macrophage colony-stimulating factor 1 (CSF-1), cathepsin D (CTSD), and myeloperoxidase (MPO). We were not able to extend the analysis to tumour necrosis factor receptor superfamily member 5 (TNFRSF5, commonly known as CD40) as this protein was not measured in the Icelandic GWAS. We also performed additional analyses for IL-34 and CSF-1 receptor (CSF1R), since IL-34 and CSF-1 are both ligands of CSF-1 receptor.

Statistical software

All analyses were performed in R 3.6.1. Bidirectional MR analyses were performed using the TwoSampleMR and MR-PRESSO packages (Verbanck et al., 2018). Colocalization was performed using the coloc package (Giambartolomei et al., 2014).

Results

Figure 1; Supplementary Tables S4, S5 show the relationships between plasma proteins and cognitive traits defined by the bidirectional univariable MR analyses. Among the 180 univariable MR analyses for the associations of pQTL-predicted plasma protein concentrations with cognitive traits, 115 had at least three SNPs with p-value<5 × 10−8 and r2 < 0.001 (F-statistics ranged from 30 to 7,577). For the remaining 65 univariable MR analyses (33 proteins), we selected SNPs with p-value<5 × 10−6 and r2 < 0.001 (F-statistics ranged from 21 to 1,811). We successfully identified cis-pQTL instruments for 75 circulating proteins and brain-specific cis-eQTL instruments for 50 circulating proteins (Figures 2, 3; Supplementary Tables S6, S7). Table 1 shows that genetically MPO, CSF-1, CTSD, and CD40 were associated with cognitive performance under different instrument selection criteria.

FIGURE 1

FIGURE 2

FIGURE 3

TABLE 1

ExposureInstrument typeNumber of SNPsMethodBetaSEp-valuea95% CI (lower)95% CI (upper)Heterogeneity p-valuePleiotropy p-value
Protein → General cognitive function
 CD40pQTL3IVW-FE0.0190.0060.001*0.0070.0300.696-
WM0.0190.0060.001*0.0080.030--
MR-Egger0.0240.0080.2160.0070.0400.9230.553
cis-pQTL3IVW-FE0.0170.0060.003*0.0060.0290.021-
WM0.0190.0060.001*0.0080.030--
MR-Egger0.0310.0110.2220.0090.0530.1660.332
cis-eQTL2IVW-FE0.0210.0064.4E-04*0.0090.032--
CSF-1pQTL8IVW-RE−0.0330.0084.5E-05**−0.048−0.0170.814-
WM−0.0360.0130.004*−0.061−0.011--
MR-Egger−0.0420.0250.141−0.0910.0070.7410.690
cis-pQTL1Wald ratio−0.0370.0150.013*−0.065−0.008--
CTSDpQTL7IVW-RE−0.0080.0100.409−0.0280.0110.081-
WM0.0010.0090.935−0.0170.019--
MR-Egger0.0100.0160.579−0.0220.0420.1430.234
cis-pQTL4IVW-RE0.0010.0070.828−0.0110.0140.614-
WM0.0040.0090.622−0.0130.022--
MR-Egger0.0080.0150.650−0.0220.0380.4650.653
cis-eQTL3IVW-FE0.0050.0100.608−0.0140.0240.945-
WM0.0050.0100.603−0.0140.024--
MR-Egger0.0030.0140.852−0.0240.0300.7740.888
MPOpQTL10IVW-RE0.0200.0140.157−0.0080.0470.051-
WM0.0380.0140.006*0.0110.065--
MR-Egger0.0520.0280.103−0.0030.1070.0820.233
cis-pQTL1Wald ratio0.0450.0180.010*0.0110.080--
cis-eQTL1Wald ratio0.2190.0592.0E-04**0.1040.335--
Protein → g Factor
 CD40pQTL3IVW-FE0.0160.0060.012*0.0040.0280.332-
WM0.0160.0060.014*0.0030.028--
MR-Egger0.0150.0140.470−0.0120.0420.1390.945
cis-pQTL3IVW-FE0.0170.0060.009*0.0040.0290.469-
WM0.0160.0060.011*0.0040.029--
MR-Egger0.0140.0100.399−0.0060.0350.2410.804
cis-eQTL2IVW-FE0.0160.0070.012*0.0040.029--
CSF-1pQTL10IVW-RE−0.0230.0130.082−0.0490.0030.215-
WM−0.0350.0150.017*−0.064−0.006--
MR-Egger−0.0480.0310.161−0.1090.0130.2080.398
cis-pQTL1Wald ratio−0.0350.0160.032*−0.067−0.003--
CTSDpQTL6IVW-RE0.0180.0120.145−0.0060.0420.104-
WM0.0280.0100.008*0.0070.048--
MR-Egger0.0410.0160.0680.0090.0730.2790.149
cis-pQTL4IVW-RE0.0260.0058.4E-08**0.0170.0360.853-
WM0.0280.0100.005*0.0080.048--
MR-Egger0.0300.0170.213−0.0030.0630.7070.790
cis-eQTL3IVW-FE0.0320.0110.004*0.0100.0530.189-
WM0.0310.0110.004*0.0100.053--
MR-Egger0.0210.0250.546−0.0270.0690.1180.655
MPOpQTL9IVW-RE0.0310.0140.026*0.0040.0590.139-
WM0.0520.0160.001*0.0200.084--
MR-Egger0.0790.0230.011*0.0340.1240.4850.047
cis-pQTL1Wald ratio0.0640.0190.001*0.0280.101--
cis-eQTL1Wald ratio0.1380.0640.030*0.0130.264--

Associations of circulating proteins with cognitive performance based on Mendelian randomisation analyses.

a

“**” indicates a p-value smaller than 0.05/180, and “*” indicates a p-value smaller than 0.05.

Abbreviations: tumour necrosis factor receptor superfamily member 5 (TNFRSF5, commonly known as CD40); confidence interval (CI); macrophage colony-stimulating factor 1 (CSF-1); cathepsin D (CTSD); expression quantitative trait loci (eQTL); inverse-variance weighted fixed-effect (IVW-FE); inverse-variance weighted random-effect (IVW-RE); myeloperoxidase (MPO); protein quantitative trait loci (pQTL); standard error (SE); single-nucleotide polymorphisms (SNPs); weighted median (WM).

Supplementary Figures S1, S2 show that nominal significant associations for pQTL predicted plasma protein concentration with cognitive traits (p-value for IVW <0.05) were consistent in the direction using sensitivity methods (WM and MR-Egger). After accounting for multiple comparisons (PIVW<0.05/180), a higher concentration of pQTL-predicted CSF-1 was associated with poorer general cognitive function (βIVW = −0.03; 95% confidence interval (CI) −0.05, −0.02; PIVW = 4.5 × 10−5). A higher concentration of cis-pQTL-predicted CSF-1 was nominally associated with both poorer general cognitive function (βWald = −0.04; 95% CI -0.065, −0.008; PWald = 0.01) and lower g Factor score (βWald = −0.04; 95% CI -0.067, −0.003; PWald = 0.03). The only two brain-specific cis-eQTL for CSF-1 were both rare variants (MAF = 2%), thus MR analysis using brain-specific cis-eQTL was not performed. However, there was weak evidence for colocalization for CSF-1 with cognitive traits (Supplementary Table S8; PP.H4<0.1 for both pQTL and eQTL).

A higher concentration of pQTL-predicted IL-8 was associated with better general cognitive function (βIVW = 0.03; 95% CI 0.02, 0.05; PIVW = 2.4 × 10−5). Cis-pQTL and brain-specific cis-eQTL instruments were not identified for IL-8. We did not find supporting evidence for colocalization for IL-8 with cognitive traits (Supplementary Table S8; PP.H4<0.1 for both pQTL and eQTL).

Although pQTL-predicted CTSD plasma concentration was not associated with cognitive traits, a higher concentration of cis-pQTL-predicted CTSD was associated with a higher g Factor score after accounting for multiple comparisons (βIVW = 0.03; 95% CI 0.02, 0.04; PIVW = 8.4 × 10−8), and a higher concentration of brain-specific cis-eQTL-predicted CTSD was nominally associated with higher g Factor score (βIVW = 0.03; 95% CI 0.01, 0.05; PIVW = 0.004). MR estimates using sensitivity methods showed a consistent direction and similar effect size for these associations and MR-Egger did not suggest horizontal pleiotropy. However, colocalization for CTSD with g Factor was not found (Supplementary Table S8; PP.H4<0.1 for both pQTL and eQTL).

The pQTL, cis-pQTL, and brain-specific cis-eQTL predicted CD40 and MPO were consistently nominally associated with better cognitive performance in the MR analyses (Figure 4; Supplementary Tables S4–S6). Particularly, brain-specific cis-eQTL-predicted MPO was associated with general cognitive function after accounting for multiple comparisons (βWald = 0.22; 95% CI 0.10, 0.33; PWald = 2.4 × 10−4). PP.H4 was 0.577 for MPO pQTL with g Factor (Supplementary Table S8).

FIGURE 4

Genetically predicted cognition was nominally associated with concentrations of plasma proteins (Figure 1; Supplementary Table S5). However, we did not find evidence supporting an influence of genetically predicted cognition on plasma concentrations of CSF-1, IL-8, CTSD, CD40, and MPO.

Using genetic associations from the Icelandic GWAS, a higher concentration of cis-pQTL-predicted MPO was nominally associated with better cognitive performance (Supplementary Table S9) measured as the general cognitive function (βIVW = 0.02; 95% CI 0.003, 0.04; PIVW = 0.02) or g Factor (βIVW = 0.02; 95% CI 0.002, 0.04; PIVW = 0.03). Consistent with our main analysis based on the SCALLOP GWAS, the PP.H4 was 0.577 for MPO pQTL (Icelandic GWAS) with the g Factor. In addition, a higher concentration of pQTL-predicted IL-34 (a ligand of CSF-1 receptor) was nominally associated with better cognitive performance but there was no evidence of colocalization (PP.H4<0.1).

Discussion

By aggregating evidence from MR (with three sets of instrument selection criteria) and colocalization, our study found that a higher concentration of genetically predicted MPO was associated with better cognitive performance. Genetically predicted CSF-1, CTSD, and CD40 were associated with cognitive performance, but causal relationships were not supported by colocalization. Conversely, we found no evidence suggesting cognitive traits or their genetic liability to affect circulating concentrations of these proteins.

Our associations of circulating proteins with general cognitive function and g Factor score were not always consistent. This may be attributed to the different methods utilised to construct the two cognitive traits. While g Factor score was constructed based on seven cognitive tests using Genomic SEM, it may have captured greater contributions from reaction time and memory, given the larger sample sizes for these tests (e.g., N = 331,679 for the memory pairs-matching test, N = 330,024 for reaction time, and N = 11,356 for the matrix pattern recognition test) (de la Fuente et al., 2021). General cognitive function was constructed based on a wider range of cognitive tests using principal component analysis, which showed a stronger genetic correlation with educational attainment (Davies et al., 2018; de la Fuente et al., 2021). Nevertheless, we observed consistent associations for genetically predicted MPO with both these two cognitive traits in our MR analyses.

In this study, we observed a direct association between circulatory concentrations of MPO and cognitive function. In agreement with this finding, studies have reported an association of a functional polymorphism in the promoter region of the MPO gene (rs2333227, G-463A) with cognitive function and risk of Alzheimer’s disease (Crawford et al., 2001; Combarros et al., 2002; Pope et al., 2006; Talarowska et al., 2015). The MPO AA genotype, which decreases the production of myeloperoxidase, was associated with cognitive decline among older adults (aged 70–79 years) (Pope et al., 2006). However, both MPO expression and protein levels in middle-aged adults (aged 20–67 years, mean = 40) were found to be associated with worse cognitive function as assessed using Trail Making Test, Stroop Test, Verbal Fluency Test, and Auditory-Verbal Learning Test (Talarowska et al., 2015). Inconsistent findings have been reported for MPO and the risk of Alzheimer’s disease. The MPO GG genotype (increasing production of myeloperoxidase) was associated with a higher risk for Alzheimer’s disease among Caucasians but not Hispanics (Crawford et al., 2001). Sex-specific risk for Alzheimer’s disease among individuals with the MPO GG genotype also has been reported (Reynolds et al., 1999). In addition, circulating MPO was higher among individuals with mild cognitive impairment and Alzheimer’s patients compared to healthy controls (Folkersen et al., 2020). This may suggest the role of MPO in cognitive performance varies by age and sex. We were not able to perform sex-specific MR analysis given that sex-specific GWAS for Alzheimer’s disease with a large sample size was not available. A higher concentration of circulating MPO was associated with an increased incidence and poorer prognosis of CVD (Ramachandra et al., 2020). Further investigation is needed to elucidate the mechanisms underlying the effects of MPO on cognition and CVD.

We also found a direct association between circulatory concentrations of CD40 and cognitive performance. A published study by Ye et al. (2019) showed that the CD40 concentrations in the cerebrospinal fluid were lower among patients with mild Alzheimer’s disease compared with healthy controls and MCI patients. Increased expression of CD40 is associated with microglial activation, which has been found to contribute to cognitive impairments and the pathogenesis of Alzheimer’s disease (Hamelin et al., 2018; Zhang et al., 2021). Inhibition of the CSF-1 receptor reduces microglial activation (Olmos-Alonso et al., 2016). A previous study reported a higher concentration of plasma CTSD among AD patients with less severe cognitive impairment (Kim et al., 2021). Increased expression of CTSD, as lysosomal protease that degrades both amyloid-beta and tau proteins, is part of an adaptive response to AD-related neurodegenerative pathology (Cataldo et al., 1995; Suire et al., 2020). Our findings of higher genetically predicted CD40 and CTSD associated with better cognitive performance and higher genetically predicted CSF-1 with poorer general cognitive function are consistent with this. CD40, CSF-1, and CTSD also have been reported to play a role in CVD, possibly through their roles in mediating inflammation (Gorelick, 2010; Ozawa et al., 2017; Sjaarda et al., 2018; Daub et al., 2020; Hoes et al., 2020). This supports the shared mechanisms underlying CVD, cognitive impairment, and Alzheimer’s disease. Further research on repurposing CVD drugs for cognitive impairment and Alzheimer’s disease is warranted.

In brief, our findings highlight a number of CVD-related proteins that may play a role in building the cognitive reserve or cognitive decline based on MR analyses with rigorously selected genetic instruments. We improved the validity of genetic instruments by using stringent p-value thresholds and performed sensitivity analyses using cis-pQTL located near the protein coding gene of interest and brain-specific cis-eQTL leveraging gene expression in brain tissues. Sensitivity MR methods showed consistent results and bidirectional MR did not suggest reverse causation. However, we recognise the limitations of our analyses. First, our analysis was based on the GWAS of circulatory protein concentrations measured using Olink CVD I panel. Large GWAS of other CVD-related proteins measured using CVD II or III panel was not available at the time of analysis. Second, the effect sizes from analyses based on the Olink protein GWAS and the Icelandic GWAS (SomaScan assay) cannot be compared directly. However, MPO measured by the two technologies are expected to be highly correlated (Pietzner et al., 2021). Third, colocalization provides moderate evidence for a common causal variant between MPO or CD40 and cognitive traits (PP.H4: 0.3–0.6). However, this is likely to be hindered by the limited power of colocalization. Our analysis showed that it was more likely for the candidate causal proteins and cognitive traits to share a single causal variant than have different causal variants (PP.H4>PP.H3). Fourth, we observed opposite effects of the two ligands of CSF-1 receptor (CSF-1 and IL-34) on general cognitive function, but we were not able to disentangle which interventions would be most effective in the clinical setting. Animal experiments and intervention trials should be considered to bring light to this. In addition, further research is warranted to decipher the role of protein-protein interaction networks in cognitive function. Fifth, this study was based on individuals of European ancestry, thus our findings may not be generalized to other ethnic groups.

Conclusion

By applying MR analysis and colocalization with a rigorous selection of genetic instrument accounting for both pQTL and brain-specific eQTL, our study investigated the relationship of a wide range of plasma proteins that were known to be related to cardiovascular risk with cognitive traits and found supporting evidence for MPO, CSF-1, CTSD, and CD40. These proteins are involved in shared pathways between CVD and those for cognitive reserve or affecting cognitive decline, suggesting they might be used as therapeutic targets for cognitive impairment. Functional work and further investigations in more relevant tissues are needed.

Statements

Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

JH designed the study, curated the data, performed the formal analysis, and drafted the manuscript. IT and AD contributed to the design of the study and methodology. All authors reviewed and edited the manuscript.

Acknowledgments

PM acknowledges generous personal support from the Edmond J Safra Foundation and Lily Safra and an NIHR Senior Investigator Award. PM and PE are supported by the UK Dementia Research Institute, which receives its funding from UK DRI Ltd., funded by the UK Medical Research Council, Alzheimer’s Society and Alzheimer’s Research UK. All of the authors are grateful for support from the Imperial College Healthcare Trust (ICHT) NIHR Biomedical Research Centre. JH received funding from the National Medical Research Council Singapore under the Open Fund—Young Individual Research Grant

Conflict of interest

DG is employed part-time by Novo Nordisk. PM has received consultancy fees from Novartis and Biogen. PM also has received honoraria or speakers’ fees from Novartis and Biogen and has received research or educational funds from Biogen, Novartis, Merck and Bristol Myers Squibb.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1124431/full#supplementary-material

Abbreviation

CVD, Cardiovascular disease; CTSD, Cathepsin D; eQTL, Expression quantitative trait loci; g Factor, General factor of intelligence; GWAS, Genome-wide association studies; Genomic SEM, Genomic structural equation modelling; IVW, Inverse-variance weighted; CSF-1, Macrophage colony-stimulating factor 1; MR, Mendelian randomisation; MAF, Minor allele frequency; MPO, Myeloperoxidase; NPX, Normalized protein expression; PP, Posterior probability; pQTL, Protein quantitative trait loci; PEA, Proximity extension assay; SNPs, Single-nucleotide polymorphisms; TNFRSF5, commonly known as CD40, Tumour necrosis factor receptor superfamily member 5; WM, Weighted median.

References

  • 1

    AssarssonE.LundbergM.HolmquistG.BjorkestenJ.ThorsenS. B.EkmanD.et al (2014). Homogenous 96-plex PEA immunoassay exhibiting high sensitivity, specificity, and excellent scalability. PLoS One9, e95192. 10.1371/journal.pone.0095192

  • 2

    AttemsJ.JellingerK. A. (2014). The overlap between vascular disease and Alzheimer's disease--lessons from pathology. BMC Med.12, 206. 10.1186/s12916-014-0206-2

  • 3

    BowdenJ.Davey SmithG.HaycockP. C.BurgessS. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol.40, 304314. 10.1002/gepi.21965

  • 4

    BowdenJ.Del GrecoM. F.MinelliC.Davey SmithG.SheehanN.ThompsonJ. (2017). A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat. Med.36, 17831802. 10.1002/sim.7221

  • 5

    BurgessS.SmallD. S.ThompsonS. G. (2017). A review of instrumental variable estimators for Mendelian randomization. Stat. Methods Med. Res.26, 23332355. 10.1177/0962280215597579

  • 6

    BurgessS.ThompsonS. G. (2017). Interpreting findings from Mendelian randomization using the MR-Egger method. Eur. J. Epidemiol.32, 377389. 10.1007/s10654-017-0255-x

  • 7

    CataldoA. M.BarnettJ. L.BermanS. A.LiJ.QuarlessS.BursztajnS.et al (1995). Gene expression and cellular content of cathepsin D in alzheimer's disease brain: Evidence for early up-regulation of the endosomal-lysosomal system. Neuron14, 671680. 10.1016/0896-6273(95)90324-0

  • 8

    CombarrosO.InfanteJ.LlorcaJ.PenaN.Fernandez-ViaderoC.BercianoJ. (2002). The myeloperoxidase gene in alzheimer's disease: A case-control study and meta-analysis. Neurosci. Lett.326, 3336. 10.1016/s0304-3940(02)00303-8

  • 9

    CongdonE. E.SigurdssonE. M. (2018). Tau-targeting therapies for Alzheimer disease. Nat. Rev. Neurol.14, 399415. 10.1038/s41582-018-0013-z

  • 10

    CrawfordF. C.FreemanM. J.SchinkaJ. A.MorrisM. D.AbdullahL. I.RichardsD.et al (2001). Association between Alzheimer's disease and a functional polymorphism in the Myeloperoxidase gene. Exp. Neurol.167, 456459. 10.1006/exnr.2000.7560

  • 11

    DaubS.LutgensE.MunzelT.DaiberA. (2020). CD40/CD40L and related signaling pathways in cardiovascular health and disease-the pros and cons for cardioprotection. Int. J. Mol. Sci.21, 8533. 10.3390/ijms21228533

  • 12

    DaviesG.LamM.HarrisS. E.TrampushJ. W.LucianoM.HillW. D.et al (2018). Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function. Nat. Commun.9, 2098. 10.1038/s41467-018-04362-x

  • 13

    De La FuenteJ.DaviesG.GrotzingerA. D.Tucker-DrobE. M.DearyI. J. (2021). A general dimension of genetic sharing across diverse cognitive traits inferred from molecular data. Nat. Hum. Behav.5, 4958. 10.1038/s41562-020-00936-2

  • 14

    FerkingstadE.SulemP.AtlasonB. A.SveinbjornssonG.MagnussonM. I.StyrmisdottirE. L.et al (2021). Large-scale integration of the plasma proteome with genetics and disease. Nat. Genet.53, 17121721. 10.1038/s41588-021-00978-w

  • 15

    FolkersenL.GustafssonS.WangQ.HansenD. H.HedmanÅ. K.SchorkA.et al (2020). Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals. Nat. Metab.2, 11351148. 10.1038/s42255-020-00287-2

  • 16

    GiambartolomeiC.VukcevicD.SchadtE. E.FrankeL.HingoraniA. D.WallaceC.et al (2014). Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet.10, e1004383. 10.1371/journal.pgen.1004383

  • 17

    GillD.GeorgakisM. K.WalkerV. M.SchmidtA. F.GkatzionisA.FreitagD. F.et al (2021). Mendelian randomization for studying the effects of perturbing drug targets. Wellcome Open Res.6, 16. 10.12688/wellcomeopenres.16544.2

  • 18

    GorelickP. B. (2010). Role of inflammation in cognitive impairment: Results of observational epidemiological studies and clinical trials. Ann. N. Y. Acad. Sci.1207, 155162. 10.1111/j.1749-6632.2010.05726.x

  • 19

    GrotzingerA. D.RhemtullaM.De VlamingR.RitchieS. J.MallardT. T.HillW. D.et al (2019). Genomic structural equation modelling provides insights into the multivariate genetic architecture of complex traits. Nat. Hum. Behav.3, 513525. 10.1038/s41562-019-0566-x

  • 20

    HamelinL.LagardeJ.DorotheeG.PotierM. C.CorlierF.KuhnastB.et al (2018). Distinct dynamic profiles of microglial activation are associated with progression of Alzheimer's disease. Brain141, 18551870. 10.1093/brain/awy079

  • 21

    HarrisS. E.CoxS. R.BellS.MarioniR. E.PrinsB. P.PattieA.et al (2020). Neurology-related protein biomarkers are associated with cognitive ability and brain volume in older age. Nat. Commun.11, 800. 10.1038/s41467-019-14161-7

  • 22

    HoJ. E.LyassA.CourchesneP.ChenG.LiuC.YinX.et al (2018). Protein biomarkers of cardiovascular disease and mortality in the community. J. Am. Heart Assoc.7, e008108. 10.1161/JAHA.117.008108

  • 23

    HoesM. F.TrompJ.OuwerkerkW.BomerN.Oberdorf-MaassS. U.SamaniN. J.et al (2020). The role of cathepsin D in the pathophysiology of heart failure and its potentially beneficial properties: A translational approach. Eur. J. Heart Fail22, 21022111. 10.1002/ejhf.1674

  • 24

    HuangL.-K.ChaoS.-P.HuC.-J. (2020). Clinical trials of new drugs for Alzheimer disease. J. Biomed. Sci.27, 18. 10.1186/s12929-019-0609-7

  • 25

    Juul RasmussenI.RasmussenK. L.NordestgaardB. G.Tybjærg-HansenA.Frikke-SchmidtR. (2020). Impact of cardiovascular risk factors and genetics on 10-year absolute risk of dementia: Risk charts for targeted prevention. Eur. Heart J.41, 40244033. 10.1093/eurheartj/ehaa695

  • 26

    KametaniF.HasegawaM. (2018). Reconsideration of amyloid hypothesis and tau hypothesis in alzheimer's disease. Front. Neurosci.12, 25. 10.3389/fnins.2018.00025

  • 27

    KimJ. W.JungS. Y.KimY.HeoH.HongC. H.SeoS. W.et al (2021). Identification of cathepsin D as a plasma biomarker for alzheimer's disease. Cells10, 138. 10.3390/cells10010138

  • 28

    LawlorD. A.HarbordR. M.SterneJ. A.TimpsonN.Davey SmithG. (2008). Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Stat. Med.27, 11331163. 10.1002/sim.3034

  • 29

    MorrisG. P.ClarkI. A.VisselB. (2018). Questions concerning the role of amyloid-β in the definition, aetiology and diagnosis of Alzheimer’s disease. Acta Neuropathol.136, 663689. 10.1007/s00401-018-1918-8

  • 30

    O'brienR. J.WongP. C. (2011). Amyloid precursor protein processing and Alzheimer's disease. Annu. Rev. Neurosci.34, 185204. 10.1146/annurev-neuro-061010-113613

  • 31

    Olmos-AlonsoA.SchettersS. T.SriS.AskewK.MancusoR.Vargas-CaballeroM.et al (2016). Pharmacological targeting of CSF1R inhibits microglial proliferation and prevents the progression of Alzheimer's-like pathology. Brain139, 891907. 10.1093/brain/awv379

  • 32

    OzawaM.ShipleyM.KivimakiM.Singh-ManouxA.BrunnerE. J. (2017). Dietary pattern, inflammation and cognitive decline: The Whitehall II prospective cohort study. Clin. Nutr.36, 506512. 10.1016/j.clnu.2016.01.013

  • 33

    PaciaroniM.BogousslavskyJ. (2013). Connecting cardiovascular disease and dementia: Further evidence. J. Am. Heart Assoc.2, e000656. 10.1161/JAHA.113.000656

  • 34

    PietznerM.WheelerE.Carrasco-ZaniniJ.KerrisonN. D.OertonE.KopruluM.et al (2021). Synergistic insights into human health from aptamer- and antibody-based proteomic profiling. Nat. Commun.12, 6822. 10.1038/s41467-021-27164-0

  • 35

    PopeS. K.KritchevskyS. B.AmbrosoneC.YaffeK.TylavskyF.SimonsickE. M.et al (2006). Myeloperoxidase polymorphism and cognitive decline in older adults in the health, aging, and body composition study. Am. J. Epidemiol.163, 10841090. 10.1093/aje/kwj146

  • 36

    RamachandraC. J. A.JaK.ChuaJ.CongS.ShimW.HausenloyD. J. (2020). Myeloperoxidase as a multifaceted target for cardiovascular protection. Antioxid. Redox Signal32, 11351149. 10.1089/ars.2019.7971

  • 37

    ReynoldsW. F.RheesJ.MaciejewskiD.PaladinoT.SieburgH.MakiR. A.et al (1999). Myeloperoxidase polymorphism is associated with gender specific risk for Alzheimer's disease. Exp. Neurol.155, 3141. 10.1006/exnr.1998.6977

  • 38

    SantosC. Y.SnyderP. J.WuW.-C.ZhangM.EcheverriaA.AlberJ. (2017). Pathophysiologic relationship between alzheimer's disease, cerebrovascular disease, and cardiovascular risk: A review and synthesis. Alzheimer's Dementia Diagnosis, Assess. Dis. Monit.7, 6987. 10.1016/j.dadm.2017.01.005

  • 39

    SjaardaJ.GersteinH.ChongM.YusufS.MeyreD.AnandS. S.et al (2018). Blood CSF1 and CXCL12 as causal mediators of coronary artery disease. J. Am. Coll. Cardiol.72, 300310. 10.1016/j.jacc.2018.04.067

  • 40

    SturchioA.DwivediA. K.YoungC. B.MalmT.MarsiliL.SharmaJ. S.et al (2021). High cerebrospinal amyloid-beta 42 is associated with normal cognition in individuals with brain amyloidosis. EClinicalMedicine38, 100988. 10.1016/j.eclinm.2021.100988

  • 41

    SuireC. N.Abdul-HayS. O.SaharaT.KangD.BrizuelaM. K.SaftigP.et al (2020). Cathepsin D regulates cerebral Aβ42/40 ratios via differential degradation of Aβ42 and Aβ40. Alzheimer's Res. Ther.12, 80. 10.1186/s13195-020-00649-8

  • 42

    TalarowskaM.SzemrajJ.GaleckiP. (2015). Myeloperoxidase gene expression and cognitive functions in depression. Adv. Med. Sci.60, 15. 10.1016/j.advms.2014.06.001

  • 43

    ThomasK. R.BangenK. J.WeigandA. J.EdmondsE. C.WongC. G.CooperS.et al (2020). Objective subtle cognitive difficulties predict future amyloid accumulation and neurodegeneration. Neurology94, e397e406. 10.1212/WNL.0000000000008838

  • 44

    VerbanckM.ChenC. Y.NealeB.DoR. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet.50, 693698. 10.1038/s41588-018-0099-7

  • 45

    WallentinL.ErikssonN.OlszowkaM.GrammerT. B.HagströmE.HeldC.et al (2021). Plasma proteins associated with cardiovascular death in patients with chronic coronary heart disease: A retrospective study. PLOS Med.18, e1003513. 10.1371/journal.pmed.1003513

  • 46

    YaffeK.BahorikA. L.HoangT. D.ForresterS.JacobsD. R.LewisC. E.et al (2020). Cardiovascular risk factors and accelerated cognitive decline in midlife: The CARDIA Study. CARDIA Study95, e839e846. 10.1212/WNL.0000000000010078

  • 47

    YeX.ZhouW.ZhangJ.for Alzheimer’s Disease Neuroimaging Initiative (2019). Association of CSF CD40 levels and synaptic degeneration across the Alzheimer’s disease spectrum. Neurosci. Lett.694, 4145. 10.1016/j.neulet.2018.11.019

  • 48

    ZhangD.LiS.HouL.JingL.RuanZ.PengB.et al (2021). Microglial activation contributes to cognitive impairments in rotenone-induced mouse Parkinson’s disease model. J. Neuroinflammation18, 4. 10.1186/s12974-020-02065-z

  • 49

    ZuberV.GrinbergN. F.GillD.ManipurI.SlobE. a. W.PatelA.et al (2022). Combining evidence from Mendelian randomization and colocalization: Review and comparison of approaches. Am. J. Hum. Genet.109, 767782. 10.1016/j.ajhg.2022.04.001

Summary

Keywords

proteins, cardiovascular disease, cognition, mendelian randomisation, myeloperoxidase (MPO)

Citation

Huang J, Gill D, Zuber V, Matthews PM, Elliott P, Tzoulaki I and Dehghan A (2023) Circulatory proteins relate cardiovascular disease to cognitive performance: A mendelian randomisation study. Front. Genet. 14:1124431. doi: 10.3389/fgene.2023.1124431

Received

18 January 2023

Accepted

08 February 2023

Published

17 February 2023

Volume

14 - 2023

Edited by

Marijana Vujkovic, University of Pennsylvania, United States

Reviewed by

Haoxiang Cheng, Icahn School of Medicine at Mount Sinai, United States

Manoj Kumar Gupta, Helmholtz Association of German Research Centers (HZ), Germany

Updates

Copyright

*Correspondence: Jian Huang,

This article was submitted to Applied Genetic Epidemiology, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics