- 1The First Clinical College of Chinese Medicine, Hunan University of Chinese Medicine, Changsha, China
- 2International Medical Department, The First Affiliated Hospital of Hunan University of Chinese Medicine, Changsha, China
- 3College of Integrative Medicine, Hunan University of Chinese Medicine, Changsha, China
- 4Health Management Department, The First Affiliated Hospital of Hunan University of Chinese Medicine, Changsha, China
- 5Department of Cardiovascular Medicine, The First Affiliated Hospital of Hunan University of Chinese Medicine, Changsha, China
- 6Preventive Treatment Center, The First Affiliated Hospital of Hunan University of Chinese Medicine, Changsha, China
- 7Intensive Care Unit, The First Affiliated Hospital of Hunan University of Chinese Medicine, Changsha, China
Background: Chronic heart failure (CHF) patients exhibit alterations in cerebral cortical structure and cognitive function. However, the mechanisms by which CHF affects cortical structure and functional regions remain unknown. This study aims to investigate potential causal relationship between CHF and cerebral cortical structure through Mendelian randomization (MR).
Methods: The research utilized genome-wide association studies (GWAS) to explore the causal association between CHF and cerebral cortical structure. The results were primarily analyzed using the inverse-variance weighted (IVW). The reliability of the data was verified through horizontal pleiotropy and heterogeneity analysis by MR-Egger intercept test and Cochran's Q-test, respectively. Replication analysis was conducted in the Integrative Epidemiology Unit (IEU) OpenGWAS project for further validation. In addition, we collected mediator genes that mediate causality to reveal potential mechanisms. Integrated bioinformatics analysis was conducted using the Open Target Genetics platform, the STRING database, and Cytoscape software.
Results: The IVW results did not reveal any significant causal association between genetically predicted CHF and the overall structure of the cerebral cortex or the surface area (SA) of the 34 functional regions of the cerebral cortex (P > 0.05). However, the results revealed that CHF increased the thickness (TH) of pars opercularis (IVW: β = 0.015, 95% CI: 0.005–0.025, P = 3.16E-03). Replication analysis supported the causal association between CHF and pars opercularis TH (IVW: β = 0.02, 95% CI: 0.010–0.033, P = 1.84E-04). We examined the degree centrality values of the top 10 mediator genes, namely CDKN1A, CELSR2, NME5, SURF4, PSMA5, TSC1, RPL7A, SURF6, PRDX3, and FTO.
Conclusion: Genetic evidence indicates a positive correlation between CHF and pars opercularis TH.
1 Introduction
Chronic heart failure (CHF) is the severe and final stage in the development of most cardiovascular diseases (CVDs), characterized by dyspnea, decreased exercise tolerance, and edema in patients. According to epidemiological surveys, there are currently over 64 million patients with CHF worldwide (1). Within 5 years, the readmission and mortality rates of patients can be as high as 80% and 50%, respectively (2). With the aging trend of society, the current prevalence of CHF among Chinese residents aged ≥35 years is 1.3%, and the number of patients is as high as 8.9 million (3, 4). By 2030, the prevalence rate in the elderly population is expected to reach 8.5% (5). Statistics indicate that CHF patients in China are hospitalized an average of 3.3 times per year, with a per capita cost of $4,982 for inpatient and outpatient treatment (6). Despite the progress that has been made in treating CHF, plenty of patients still experience unsatisfactory outcomes and poor prognosis.
The cerebral cortex is a layer of gray matter that covers the surface of the cerebral hemispheres and is composed primarily of neuronal cells (7). It is widely acknowledged that the structure of the cerebral cortex is responsible for a number of higher cognitive functions in humans, including the processes of thinking, learning, memory, and language (8, 9). The human cerebral cortex is characterized by two main parameters: surface area (SA) and thickness (TH), which are regulated by multiple genes (10). In-depth studies related to the cerebral cortex can enhance our comprehension of alterations in cortical structure during disease progression and throughout the lifespan. The “heart-brain axis” refers to the bidirectional communication network between the heart and the central nervous system (CNS), involving hemodynamic changes and neuronal signaling. This concept involves a tight heart-brain feedback interaction and is highly complex (11, 12). Researches have demonstrated that more than 50% of patients with CHF experience secondary brain damage, which can manifest as autonomic damage, cognitive dysfunction, and neuropsychological deficits (13). In addition to their relation to the grey matter structure of the subcortex and brainstem, these varying degrees of brain damage also involve the integrity of cerebral cortical structure (14). Several scientific studies have examined the structure of the cerebral cortex, and the cortical regions of patients with CHF have been roughly examined using magnetic resonance imaging (MRI) (15–17). However, the assessment of altered cerebral cortical structure in specific regions has not been effective. The causal associations between CHF and cerebral cortex remain unknown.
Mendelian randomization (MR) analysis, as a naturalistic randomized controlled trial, is an effective approach to evaluate potential causal associations between specific exposure factors and outcomes. In contrast to randomized controlled trials (RCTs), MR analysis employs single nucleotide polymorphism (SNP) as an instrumental variable (IV) based on Mendel's law of independent assortment. This approach can effectively aid in establishing a causal association between phenotype and disease (18, 19). Additionally, MR analysis can address the limitations of traditional epidemiological studies that are susceptible to potential confounding factors. As a result, MR study allows for obtaining more plausible causal associations and a stronger ability to argue for etiological inferences (20). Currently, the joint application of genome-wide association studies (GWAS) and biological big data has become a new trend in scientific research. This approach provides an opportunity to explore etiological associations from a genetic perspective and lays the foundation for the wide application of MR analysis (21). In recent years, researchers from both domestic and international settings have investigated the impact of diseases, circulating biomarkers, and behaviors on the structure of the cerebral cortex through the use of MR studies. These studies have confirmed that obesity (body mass index (BMI) and waist-to-hip ratio (WHR)), elevated blood lipid levels, and sleep disorders (insomnia and shorter sleep duration) are associated with alterations in the structure of the cerebral cortex (22–24). In this study, we employed publicly available GWAS summary statistics to preliminarily investigate the effects of CHF on the structure of the cerebral cortex and to explore its correlation with cerebral cortical SA and TH. Subgroup analyses were conducted according to different functional regions of the cerebral cortex, providing new insights into the effects of CHF on cerebral cortical structure and its underlying mechanisms.
2 Materials and methods
2.1 Study design
This study was conducted in accordance with the established standards for reporting on two-sample MR analyses (25). The genetic IVs for phenotypes should meet the three crucial assumptions (26): (a) There is a robust and strong correlation between IVs and CHF. (b) IVs are not affected by confounding factors that may influence the relationship between CHF and cerebral cortical structure. (c) IVs affect cerebral cortical structure only through the onset of CHF, but not through other pathways. The study investigated CHF as the exposure factor and measured three outcomes in the following order: overall cortical structure, SA of functional regions of the cerebral cortex, and TH of functional regions of the cerebral cortex. We initially conducted MR studies with the objective of comprehensively analyzing the causal associations between CHF and the cerebral cortex. Subsequently, we proceeded to perform replication analysis with the intention of validating the identified causal association. Figure 1 shows the study design flowchart.
Figure 1 Research design flow chart. CHF, chronic heart failure; SA, surface area; TH, thickness; SNP, single nucleotide polymorphism; MR, Mendelian randomization; IVW, inverse-variance weighted.
2.2 Data sources
2.2.1 Exposure: CHF
The GWAS summary statistics for CHF were obtained from the Cardiovascular Disease Knowledge Portal (CVDKP) database, which included 26 cohort studies (27). The study consisted of 17 population-based cohorts that participated in the Heart Failure Molecular Epidemiology for Therapeutic Targets Consortium (HERMES), comprising 38,780 cases and 893,657 controls, as well as nine case-control studies with 8,529 cases and 36,357 controls. The study participants in this GWAS study were of European ancestry, totaling 47,309 cases and 93,014 healthy controls. In this meta-analysis study, Shah et al. employed logistic regression (LR) to analyze the association between autosomal SNPs and heart failure (HF). The overall effect model was evaluated using inverse-variance weighted (IVW). SNPs exhibiting linkage disequilibrium (LD) were excluded based on the principle of “P < 5 × 10−8, r 2 < 0.1”. The researchers ultimately identified 12 independent genetic variants associated with HF risk by analyzing 11 genomic loci, providing new insights into the understanding of HF etiology.
2.2.2 Outcome: cerebral cortex structure
The outcome dataset utilized in this study was derived from the GWAS summary statistics published in March 2020 by the Enhancing Neuro Imaging Genetics Through Meta-Analysis (ENIGMA) consortium (http://enigma.ini.usc.edu). The study performed a meta-analysis of MRI data from 51,665 individuals from 60 cohorts worldwide (28). The GWAS summary statistics included the total cerebral cortex SA and average TH, as well as 34 cortical regions with known function defined according to the Desikan-Killiany cortical atlas. The subject population was approximately 94% of European ancestry. Grasby et al. identified a total of 70 categories of cortical phenotypes and 306 statistically significant loci through GWAS analyses and multiple statistical comparisons (P < 5 × 10−8). After multiple statistical correction tests, the study confirmed that only 187 genomic loci were associated with cortical SA and 12 genomic loci were associated with cortical TH. The proportion of SA phenotypic variation that could be explained by common genetic variants was greater than that of TH phenotypic variation, both in terms of the overall structure of the cerebral cortex and in terms of the 34 specific functional regions.
2.2.3 CHF replication samples
The CHF replication samples were derived from GWAS summary statistics (ebi-a-GCST90018806) included in September 2021 by the Integrative Epidemiology Unit (IEU) OpenGWAS project (https://gwas.mrcieu.ac.uk) (29). A total of 14,262 cases and 47,189 controls were included in the study, with data from BioBank Japan (BBJ), the UK Biobank (UKB) and FinnGen. The CHF replication samples comprised 10,540 cases and 168,186 controls derived from BioBank Japan, and 6,526 cases and 350,289 controls derived from the UK Biobank and FinnGen databases. The subject population consisted of individuals of Japanese and European ancestry. In this study, Sakaue et al. conducted a genome-wide association study of 220 human phenotypes from BioBank Japan, successfully replicating 94.2% of the variant loci in the European population. A cross-population meta-analysis was conducted with 196 human phenotypes from the UK Biobank and 128 human phenotypes from FinnGen. The meta-analysis identified 1,730 disease-associated loci, 12,066 biomarker-associated loci, and 1,018 drug-associated loci, respectively, which extended the genetic association map in non-European populations.
2.3 Selection of instrumental variables
Instrumental variables significantly associated with CHF were screened from the GWAS summary statistics with a threshold of “P < 5 × 10−8”. SNPs that were not affected by linkage disequilibrium and were independent of each other were screened using “r² < 0.001 and Clump distance >10,000 kb” as thresholds (30). Additionally, the formula was employed to calculate the F statistic in order to detect any potential bias due to weak IVs. Only IVs with F statistics greater than 10 were included in this study. In accordance with the exclusivity assumption of MR analysis, SNPs associated with cerebral cortical structure should be excluded (P < 5 × 10−8). Regarding the independence assumption, secondary phenotypes were identified for each SNP using the PhenoScanner V2 website (http://www.phenoscanner.medschl.cam.ac.uk/) (31), and SNPs associated with confounding factors (type 2 diabetes, body mass index, blood pressure, and coronary heart disease) (32) were excluded.
2.4 Functional exploration of mediator genes
In order to identify the core mediator genes regulating the causal association between CHF and cerebral cortex structure, and to reveal the potential biological pathways between CHF and brain dysfunction, we conducted an integrated bioinformatics analysis. Genes that are functionally associated with genetic variation are referred to as mediator genes. The Open Target Genetics platform (33) is a free and open-source tool that highlights statistical evidence centered on genetic variation. It utilizes human genetics and genomics data for systematic drug target identification and prioritization, exploring the intrinsic associations between traits, variants, and genes, and identifying potential drug targets (34). The Open Target Genetics platform (https://genetics.opentargets.org/) was employed to genetically annotate the genetic variants utilized as IVs, with the subsequent construction of a protein-protein interaction (PPI) network based on these mediator genes in the STRING database (version 12.0) (https://cn.string-db.org/). The Cytoscape software (version 3.9.1) was used to visualize and analyze the complex network based on the node degree and integrated score of target proteins in the PPI network. The topology of the network graph was analyzed using the CytoHubba plug-in. Each mediator gene was assigned a value by the topological network algorithm, and the hub genes and the subnetwork were identified after sorting. We conducted Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis using the Bioinformatics Online Cloud Platform (https://www.bioinformatics.com.cn/) in order to present the results of the functional annotation of mediator genes in a more intuitive manner.
2.5 Statistical analysis
R studio software (version 4.3.2) and the TwoSampleMR package (version 0.5.6) were used for statistical analysis. We utilized selected SNPs as IVs to evaluate the causal relationship between CHF and cerebral cortex structure. Five methods were employed: MR-Egger, inverse-variance weighted, weighted median, simple mode, and weighted mode. The results of the MR study were mainly dominated by the IVW analysis, while the MR-Egger method and the weighted median approach were able to enhance the estimation accuracy of the IVW analysis. The cluster-heatmap format of IVW analysis results was generated by the Bioinformatics Online Cloud Platform. The replication MR analysis procedure was identical to that employed in the discovery stage.
Instrumental variables from different analytical platforms, experiments, populations, etc., may be heterogeneous and thus affect the results of MR analyses. In addition, if IVs influence the occurrence of an outcome through factors other than the exposure factor, this indicates that the IVs are pleiotropic. Pleiotropy can lead to a violation of the assumptions of independence and exclusivity. Therefore, if there is a causal relationship between CHF and cerebral cortex structure, it is necessary to confirm the reliability of the findings through the test of heterogeneity and horizontal pleiotropy. In this study, Cochran's Q-test was used to test for heterogeneity among the included studies (P-value < 0.05 indicates heterogeneity). The MR-Egger regression test was used to determine the presence of horizontal pleiotropy (P-value < 0.05 indicates horizontal pleiotropy). A leave-one-out (LOO) analysis was employed in the sensitivity analysis to evaluate whether a single SNP had an effect on the association between CHF and cerebral cortex structure. The leave-one-out method entailed the exclusion of a single SNP one at a time, with the remaining SNPs subsequently subjected to an aggregate effect calculation.
Furthermore, the Bonferroni correction was employed to account for the multiplicity of test results. For each multiple test, the new significance level was set at P < 0.05/n, where n represents the number of results (35). Consequently, a P-value of less than 2.5 × 10−2 (0.05/2) was deemed statistically significant in the estimation of the overall structure of the cerebral cortex, whereas a P-value of less than 7.35 × 10−4 (0.05/68) was considered statistically significant in the estimation of the 34 specific functional regions of the SA and TH. A P-value of less than 0.05 was deemed suggestive of an association and a P-value of less than 0.01 was considered suggestive of a stronger association.
3 Results
3.1 Details of instrumental variables
This study screened for SNPs that differed significantly on a genome-wide basis using the completed GWAS database for CHF. A total of nine SNPs were finally included (Supplementary Table S1). The median value of the F-statistics for the IVs was 34.87, with a range of 30.89–83.10. It is noteworthy that all F-statistics were greater than 10. This suggests a significant correlation between the IVs and the exposed factors represented by CHF. The results of the MR analyses are not significantly affected by the inclusion of weak IVs.
3.2 Results of MR analysis
The study conducted comprehensive MR analyses to assess the causal associations between genetically predicted CHF and the overall structure of the cerebral cortex, SA and TH of the 34 functional regions of the cerebral cortex. The results showed no statistically significant association between CHF and the overall structure of the cerebral cortex (IVW: βSA = 247.489, 95% CISA: −1,799.191–2,294.168, PSA = 0.813; IVW: βTH = −0.009, 95% CITH: −0.019–0.001, PTH = 0.066) (Supplementary Table S2). Figure 2 presents the changes in IVW analysis results between CHF and the SA/TH of the 34 functional regions of the cerebral cortex using a cluster-heatmap format. The IVW analysis did not find a significant causal relationship between genetically predicted CHF and the SA of the 34 functional regions (P > 0.05) (Supplementary Table S3). Nevertheless, in specific functional region analyses, a potential causal relationship between CHF and pars opercularis TH was identified (IVW: β = 0.015, 95% CI: 0.005–0.025, P = 3.16E-03) (Supplementary Table S4). The IVW results indicated a positive correlation between CHF and alterations of pars opercularis TH. Specifically, the results of the weighted median method and the IVW analysis were similar and both were statistically significant (weighted median: β = 0.017, 95% CI: 0.004–0.030, P = 9.50E-03). Figures 3A,B illustrate the scatter plot and the funnel plot of the causal relationship between CHF and pars opercularis TH, respectively. The asymmetry in the funnel plot may be due to factors such as clinical or methodological heterogeneity between studies that do not introduce a significant degree of bias into the results. The results of the meta-analysis of the causal relationship between CHF and pars opercularis TH are presented in Figure 4.
Figure 2 Clustering heat maps of IVW results of causality between CHF and SA and TH of cerebral cortex. SA, surface area; TH, thickness; IVW, inverse-variance weighted.
Figure 3 Causal relationship of CHF in predicting pars opercularis TH. (A) Scatter plot; (B) funnel diagram; (C) individual SNP causal effect diagram; (D) leave-one-out sensitivity analysis. MR, Mendelian randomization; SNP, single nucleotide polymorphism; CHF, chronic heart failure.
Figure 4 Forest diagram of causal association between CHF and pars opercularis TH. MR, Mendelian randomization; CI, confidence interval.
The causal relationship between CHF and pars opercularis TH was successfully replicated in the MR analysis of the IEU OpenGWAS project. The results of the IVW analysis demonstrated a positive correlation between CHF and pars opercularis TH (β = 0.02, 95% CI: 0.010–0.033, P = 1.84E-04) (Supplementary Table S5). The replication analysis yielded results consistent with those of previous MR studies, thereby providing support for the reliability of the present findings.
3.3 Sensitivity analysis
We assessed the magnitude of heterogeneity and horizontal pleiotropy by using Cochran's Q-test and MR-Egger intercept test. The Cochran's Q-test statistic was not statistically significant (P > 0.05), indicating no heterogeneity among the SNPs associated with CHF (Supplementary Tables S6–S8). The MR-Egger intercept results indicated no horizontal pleiotropy between CHF and pars opercularis TH (P > 0.05) (Supplementary Tables S9–S11), providing strong evidence for a direct causal relationship. The leave-one-out analysis demonstrated that no SNP could potentially drive the null causal effect of CHF on pars opercularis TH (Figures 3C,D).
3.4 Results of integrated bioinformatics analysis
This study identified the dominant instrumental loci for the causal associations between CHF and cerebral cortex structure and analyzed their genetic profiles. After eliminating duplicated sequences, we matched 9 independent SNPs with genetic variations of the mediator genes, and identified 107 mediator genes. The correspondence between SNPs and mediator genes is shown in Supplementary Table S12. The PPI network was constructed using the STRING database with a screening criterion of “Minimum Required Interaction Score >0.150” (Figure 5A). The network illustrates the correlation between mutation-associated mediator genes (Supplementary Table S13). The PPI network consisted of 107 nodes and 483 edges, with an average node degree of 9.03, an average local clustering coefficient of 0.427, and an enriched P-value < 1.0E-16. The PPI network data were detailed in Supplementary Table S14. The top 10 mediator genes with the highest degree centrality values were screened, which included CDKN1A, CELSR2, NME5, SURF4, PSMA5, TSC1, RPL7A, SURF6, PRDX3, and FTO (Figure 5B).
Figure 5 Results of integrated bioinformatics analysis. (A) Protein-protein interaction network of mediator genes; (B) diagram of the top 20 mediator genes interactions; (C) histogram of GO biological function analysis; (D) bubble diagram for KEGG pathway enrichment analysis. GO, gene ontology; KEGG, Kyoto encyclopedia of genes and genomes; BP, biological process; CC, cellular component; MF, molecular function.
A total of 210 entries were obtained from GO terms analysis under the screening condition of “P < 0.05”, including 149 Biological Process (BP), 16 Cellular Component (CC), and 45 Molecular Function (MF) entries (Supplementary Table S15). Furthermore, 16 entries were obtained from KEGG pathway enrichment analysis (Supplementary Table S16). Figures 5C,D show the results of the GO terms histogram and the KEGG pathway enrichment analysis bubble plot. The mediator genes in the BP analysis were primarily concentrated in “glutathione derivative metabolic process”, “glutathione derivative biosynthetic process” and “xenobiotic catabolic process”. The CC analysis revealed that mediator genes were mainly enriched in “intercellular bridge”, “lysosomal membrane” and “lytic vacuole membrane”. Regarding MF analysis, the mediator genes were principally involved in “glutathione transferase activity”, “glutathione binding” and “oligopeptide binding”. The KEGG pathway enrichment analysis indicated that the mediator genes were primarily enriched in pathways related to “platinum drug resistance”, “glutathione metabolism” and “epatocellular carcinoma”.
4 Discussion
Patients with CHF frequently exhibit concomitant brain damage, which typically involves the integrity of cortical structure. In this study, a causal relationship between CHF and pars opercularis TH was determined through the use of two-sample MR analysis. The results of this study indicate that patients with CHF exhibit increased pars opercularis TH. The integrated bioinformatics analysis revealed that CHF may affect pars opercularis TH through 10 key mediator genes. Further KEGG enrichment analysis indicates that glutathione (GSH) metabolism may be a potential biological pathway by which CHF affects the structure of the cerebral cortex. This establishes a foundation for further elucidation of the relationship between CHF and neuropsychiatric disorders, including cognitive dysfunction, dementia, and depression.
The pars opercularis is a key component of the motor-linguistic area, situated posterior to the inferior frontal gyrus of the frontal lobe, as designated by Brodmann as area 44 (BA44) (36). The pars opercularis is essential for lexical extraction, language synthesis, memory, coordination of the oral muscles, and various cognitive functions (37–39). The normal function of the pars opercularis depends on the integrity of the cortex, while its volume is determined by both total SA and average TH. A review of the literature indicates that several factors regulate the SA and TH of the human cerebral cortex. Based on the radial unit hypothesis (40), cortical SA is determined by the number of neural progenitor cell proliferation units, while cortical TH is dependent on the number of cell divisions within each proliferation unit. In functional regions of the cerebral cortex, changes in cortical TH are primarily influenced by genetic factors and tend to decrease with age (41–43). Most neurodegenerative diseases (NDDs) and cerebrovascular diseases (CVDs) are characterized by cortical TH thinning (44). The study findings indicate that there is no significant causal relationship between CHF and the overall structure of the cerebral cortex. This may be attributed to the fact that cerebral glucose metabolism in patients with early CHF maintains a relative balance of compensatory increases and decreases, thus preventing a significant alteration in the overall structure of the cerebral cortex (45). In recent years, researchers have devoted considerable attention to the potential causal association between HF and the structure of functional cortical regions of the brain (46, 47). These findings have confirmed that frontal brain activity is reduced in HF, and that there is a significant positive correlation between left ventricular ejection fraction (LVEF) and overall density of the frontal cortex (48, 49). Notably, the conclusions drawn from subgroup analyses were opposite to the results of previous findings. Kumar et al. conducted a study that revealed a decrease in regional cortical TH in areas regulating automatic memory, cognitive, emotional, linguistic, and visual functions in patients with HF (50). Nevertheless, the present study found that CHF is associated with an increase in pars opercularis TH. This indicates that neuropsychiatric disorders, such as cognitive dysfunction, dementia, or depression, may not be the primary cause of thickening in the pars opercularis in individuals with CHF. Most patients with CHF often experience low cardiac output and sleep apnea, which can cause hypoxia/ischemia-induced brain damage. This damage can affect the TH of cortical structural regions of the brain (51). Baril et al. discovered obstructive sleep apnea generally alters gray matter structure and that high levels of hypoxemia are associated with increased TH of the left prefrontal cortex (52). Therefore, we propose the hypothesis that CHF may lead to hypoxemia, resulting in neuronal cell dysfunction, increased capillary permeability, accelerated division of neural progenitor cells, and ultimately compensatory hypertrophy of cerebral cortical structure. In addition, prolonged hypoxia in the brain stimulates the secretion of adrenaline (A) and norepinephrine (NE), leading to vasoconstriction. This exacerbates cerebral hypoxia, resulting in cerebral oedema and abnormal thickening of the cerebral cortex. In summary, further in-depth exploration is required to investigate the potential molecular mechanisms behind the complex changes that may occur in pars opercularis TH or other specific functional regions in the CHF population.
Integrated bioinformatics analysis identified 10 key genes, namely CDKN1A, CELSR2, NME5, SURF4, PSMA5, TSC1, RPL7A, SURF6, PRDX3, and FTO, which could help elucidate the potential mechanisms underlying the association between CHF and structural alterations in the cerebral cortex. The KEGG pathway enrichment analysis revealed that five mediator genes, which were closely related to the genetic variations, were enriched in the “glutathione metabolism” signaling pathway (P = 9.36E-06). GSH is a biologically active polypeptide composed of glutamic acid, cysteine, and glycine, containing a γ-amide bond and sulfhydryl groups (53). GSH is widely distributed in animals and plants and plays a crucial role in various biological processes, including DeoxyriboNucleic Acid (DNA) and protein synthesis, amino acid transport, gene expression, and cell proliferation and apoptosis (54, 55). In the context of CVDs, GSH serves as a marker for prophylactic antioxidant therapy and the risk of adverse cardiovascular events (56). Researches have demonstrated that GSH can not only scavenge reactive oxygen species via the glutathione peroxidase (GSH-Px) and glutathione S-transferase (GST) pathways, but also bind to free radicals to exert detoxifying and antioxidant effects, thereby attenuating myocardial cells apoptosis (57, 58). Several CVDs, including myocardial ischemia/reperfusion injury (MI/RI), coronary atherosclerotic heart disease (CHD), hypertension (HT), and myocardial infarction (MI), have been found to have a negative correlation with GSH levels (59–61). Meanwhile, GSH plays a vital role in the antioxidant defense system and in maintaining redox homeostasis in neuronal cells. A reduction in the concentration of GSH in neuronal cells has been demonstrated to be a significant contributing factor in the development of Alzheimer's disease (AD), Huntington's disease (HD), and Parkinson's disease (PD) (62–64). Based on this hypothesis, if GSH metabolism is abnormal and GSH levels are reduced, it may affect both the cardiovascular and central nervous systems, increasing the susceptibility of myocardial cells and neuronal cells to oxidative stress. This can trigger cellular degeneration, necrosis, and apoptosis, ultimately leading to alterations in cortical structure.
To the best of our knowledge, a similar MR study has been conducted previously. The results of the study by Hu et al. demonstrated a causal effect of HF on the SA of the caudal middle frontal lobule, insula lobule, precuneus lobule and superior parietal lobule (65). In contrast to the present study, the investigators employed a screening criterion of “P < 5 × 10−6” to identify SNPs that predicted HF as the final IVs included in the study. The discrepancy in screening criteria may have contributed to the disparate outcomes observed in the two studies. In addition to employing MR analysis, our study also investigated the potential mechanism of CHF on the structure of the cerebral cortex through integrated bioinformatics techniques. The utilization of MR analysis in this study enabled the avoidance of residual confounding factors and reverse causality, thereby compensating for the inherent limitations of observational studies. Genetic variants were identified at the time of conception and followed the principle of random assignment. In this study, genetic variants associated with CHF were employed as the exposure and cerebral cortical structure as the outcome in order to ascertain whether there was a causal effect between CHF and cerebral cortex structure. This approach effectively avoided reverse causality in this MR study. Besides, the replication analysis was conducted within the IEU OpenGWAS project with the objective of validating the MR results.
Inevitably, the current study still has certain limitations. The GWAS sample data utilized in this study were predominantly derived from populations of European ancestry. Therefore, the results cannot be used to establish causal associations between CHF and cerebral cortex structure in ethnically diverse populations. Moreover, the statistics on exposure factors selected during the study were pooled and did not categorize specific forms of CHF. The causal associations between cerebral cortical structure and CHF type and grade remains unclear. What's more, the present study did not address the severity of altered SA or TH of functional regions of the cerebral cortex. The present study employed integrated bioinformatics analysis to clarify the critical targets and potential molecular mechanisms involved in regulating the TH alteration of pars opercularis. Nevertheless, further experimental studies are required to substantiate these findings.
5 Conclusion
In summary, this study provides preliminary genetic evidence supporting a positive correlation between CHF and pars opercularis TH through MR analysis. This study contributes to the in-depth examination of the “heart-brain axis” theory. Further investigation is necessary to comprehend the association and potential mechanism of action between CHF, cerebral cortical structure, and altered brain function.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
LP: Conceptualization, Resources, Methodology, Writing – original draft. HC: Conceptualization, Software, Writing – original draft. YT: Conceptualization, Software, Writing – original draft. FZ: Data curation, Formal Analysis, Writing – original draft. YL: Data curation, Formal Analysis, Writing – original draft. ZX: Visualization, Validation, Writing – original draft. QC: Writing – review & editing, Supervision. XC: Funding acquisition, Writing – review & editing, Project administration.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article.
This study is supported by the National Famous Elderly Chinese Medicine Expert Chen Xinyu Inheritance Workshop Construction Project [National Chinese Medicine Human Education Letter (2022) No. 75]; National Natural Science Foundation of China (81704061, 81173213); Hunan Provincial Postgraduate Research and Innovation Project (CX20230806); Hunan Innovative Province Building Special Project against Novel Coronavirus Pneumonia (2020SK3009); Hunan Provincial Development and Reform Commission Innovation Guidance Special Project (Hunan Development and Reform Investment No. 2019-412); Hunan Provincial Engineering Research Centre for the “Four Seasons of Yang Regulation” for Preventive Treatment of Disease [Hunan Development and Reform High-tech (2020) No. 1006]; Hunan University of Chinese Medicine First-class Discipline Open Fund for Traditional Chinese Medicine (2022ZYX01, 2022ZYX11, 2021ZYX41); Hunan University of Chinese Medicine “Popularisation of Experimental Animal Science” Special Project (2023DWKP11).
Acknowledgments
We thank the CVDKP database, the HERMES consortium, the ENIGMA consortium, the IEU OpenGWAS project, and the Open Targets Genetics platform for sharing the GWAS summary statistics, human genetics, and genomics data. We are equally grateful to Shanghai NewCore Biotechnology Co., Ltd. (https://www.bioinformatics.com.cn) for providing data analysis and visualisation support for this study. Finally, we thank all researchers involved for their contributions.
Conflict of interest
The 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/fcvm.2024.1396311/full#supplementary-material
References
1. Bragazzi NL, Zhong W, Shu J, Abu Much A, Lotan D, Grupper A, et al. Burden of heart failure and underlying causes in 195 countries and territories from 1990 to 2017. Eur J Prev Cardiol. (2021) 28(15):1682–90. doi: 10.1093/eurjpc/zwaa147
2. Savarese G, Becher PM, Lund LH, Seferovic P, Rosano GMC, Coats AJS. Global burden of heart failure:a comprehensive and updated review of epidemiology. Cardiovasc Res. (2023) 118(17):3272–87. doi: 10.1093/cvr/cvac013
3. Huang J. Heart failure 2019–progress and reflections. Chin J Heart Fail Cardiomyopathy. (2020) 04(1):52–9. doi: 10.3760/cma.j.cn101460-20200221-00016
4. In China TWCOTROCHAD, Hu SS. Report on cardiovascular health and diseases in China 2021:an updated summary. J Geriatr Cardiol. (2023) 20(6):399–430. doi: 10.26599/1671-5411.2023.06.001
5. Roger VL. Epidemiology of heart failure:a contemporary perspective. Circ Res. (2021) 128(10):1421–34. doi: 10.1161/circresaha.121.318172
6. Wang H, Chai K, Du M, Wang S, Cai JP, Li Y, et al. Prevalence and incidence of heart failure among urban patients in China:a national population-based analysis. Circ Heart Fail. (2021) 14(10):e008406. doi: 10.1161/circheartfailure.121.008406
7. Molnár Z, Clowry G. Human cerebral cortex development. J Anat. (2019) 235(3):431. doi: 10.1111/joa.13000
8. Pang JC, Aquino KM, Oldehinkel M, Robinson PA, Fulcher BD, Breakspear M, et al. Geometric constraints on human brain function. Nature. (2023) 618(7965):566–74. doi: 10.1038/s41586-023-06098-1
9. Nakai T, Nishimoto S. Representations and decodability of diverse cognitive functions are preserved across the human cortex, cerebellum, and subcortex. Commun Biol. (2022) 5(1):1245. doi: 10.1038/s42003-022-04221-y
10. Strike LT, Hansell NK, Couvy-Duchesne B, Thompson PM, de Zubicaray GI, McMahon KL, et al. Genetic complexity of cortical structure:differences in genetic and environmental factors influencing cortical surface area and thickness. Cereb Cortex. (2019) 29(3):952–62. doi: 10.1093/cercor/bhy002
11. Doehner W, Čelutkienė J, Yilmaz MB, Coats AJS. Heart failure and the heart-brain axis. QJM. (2023) 116(11):897–902. doi: 10.1093/qjmed/hcad179
12. Wang M, Xu B, Hou X, Shi Q, Zhao H, Gui Q, et al. Altered brain networks and connections in chronic heart failure patients complicated with cognitive impairment. Front Aging Neurosci. (2023) 15:1153496. doi: 10.3389/fnagi.2023.1153496
13. Ovsenik A, Podbregar M, Fabjan A. Cerebral blood flow impairment and cognitive decline in heart failure. Brain Behav. (2021) 11(6):e02176. doi: 10.1002/brb3.2176
14. Göpfert D, Traub J, Sell R, Homola GA, Vogt M, Pham M, et al. Profiles of cognitive impairment in chronic heart failure-a cluster analytic approach. Front Hum Neurosci. (2023) 17:1126553. doi: 10.3389/fnhum.2023.1126553
15. Frey A, Homola GA, Henneges C, Mühlbauer L, Sell R, Kraft P, et al. Temporal changes in total and hippocampal brain volume and cognitive function in patients with chronic heart failure-the COGNITION.MATTERS-HF cohort study. Eur Heart J. (2021) 42(16):1569–78. doi: 10.1093/eurheartj/ehab003
16. Ferro D, van den Brink H, Amier R, van Buchem M, de Bresser J, Bron E, et al. Cerebral cortical microinfarcts:a novel MRI marker of vascular brain injury in patients with heart failure. Int J Cardiol. (2020) 310:96–102. doi: 10.1016/j.ijcard.2020.04.032
17. Doehner W, Čelutkienė J, Haeusler KG. Central command in heart failure:was there effect of hemispheric lateralization in insular cortex activation? Reply. Eur J Heart Fail. (2018) 20(9):1370–1. doi: 10.1002/ejhf.1281
18. Birney E. Mendelian randomization. Cold Spring Harb Perspect Med. (2022) 12(4):a041302. doi: 10.1101/cshperspect.a041302
19. Ference BA, Holmes MV, Smith GD. Using Mendelian randomization to improve the design of randomized trials. Cold Spring Harb Perspect Med. (2021) 11(7):a040980. doi: 10.1101/cshperspect.a040980
20. Sekula P, Del Greco MF, Pattaro C, Köttgen A. Mendelian randomization as an approach to assess causality using observational data. J Am Soc Nephrol. (2016) 27(11):3253–65. doi: 10.1681/asn.2016010098
21. Zuber V, Grinberg NF, Gill D, Manipur I, Slob EAW, Patel A, et al. Combining evidence from Mendelian randomization and colocalization:review and comparison of approaches. Am J Hum Genet. (2022) 109(5):767–82. doi: 10.1016/j.ajhg.2022.04.001
22. Chen W, Feng J, Guo J, Dong S, Li R, Ngo JCK, et al. Obesity causally influencing brain cortical structure:a Mendelian randomization study. Cereb Cortex. (2023) 33(15):9409–16. doi: 10.1093/cercor/bhad214
23. Zeng Y, Guo R, Cao S, Yang H. Causal associations between blood lipids and brain structures:a Mendelian randomization study. Cereb Cortex. (2023) 33(21):10901–8. doi: 10.1093/cercor/bhad334
24. Wang Q, Hu S, Qi L, Wang X, Jin G, Wu D, et al. Causal associations between sleep traits and brain structure:a bidirectional Mendelian randomization study. Behav Brain Funct. (2023) 19(1):17. doi: 10.1186/s12993-023-00220-z
25. Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, et al. Strengthening the reporting of observational studies in epidemiology using Mendelian randomization:the STROBE-MR statement. JAMA. (2021) 326(16):1614–21. doi: 10.1001/jama.2021.18236
26. Skrivankova VW, Richmond RC, Woolf BAR, Davies NM, Swanson SA, VanderWeele TJ, et al. Strengthening the reporting of observational studies in epidemiology using Mendelian randomisation (STROBE-MR):explanation and elaboration. Br Med J. (2021) 375:n2233. doi: 10.1136/bmj.n2233
27. Shah S, Henry A, Roselli C, Lin H, Sveinbjörnsson G, Fatemifar G, et al. Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure. Nat Commun. (2020) 11(1):163. doi: 10.1038/s41467-019-13690-5
28. Grasby KL, Jahanshad N, Painter JN, Colodro-Conde L, Bralten J, Hibar DP, et al. The genetic architecture of the human cerebral cortex. Science. (2020) 367(6484):eaay6690. doi: 10.1126/science.aay6690
29. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. (2021) 53(10):1415–24. doi: 10.1038/s41588-021-00931-x
30. Zou M, Liang Q, Zhang W, Liang J, Zhu Y, Xu Y. Diet-derived circulating antioxidants and risk of inflammatory bowel disease:a Mendelian randomization study and meta-analysis. Front Immunol. (2024) 15:1334395. doi: 10.3389/fimmu.2024.1334395
31. Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. Phenoscanner V2:an expanded tool for searching human genotype-phenotype associations. Bioinformatics. (2019) 35(22):4851–3. doi: 10.1093/bioinformatics/btz469
32. Triposkiadis F, Xanthopoulos A, Parissis J, Butler J, Farmakis D. Pathogenesis of chronic heart failure:cardiovascular aging, risk factors, comorbidities, and disease modifiers. Heart Fail Rev. (2022) 27(1):337–44. doi: 10.1007/s10741-020-09987-z
33. Carvalho-Silva D, Pierleoni A, Pignatelli M, Ong C, Fumis L, Karamanis N, et al. Open targets platform:new developments and updates two years on. Nucleic Acids Res. (2019) 47(D1):D1056–65. doi: 10.1093/nar/gky1133
34. Ochoa D, Hercules A, Carmona M, Suveges D, Gonzalez-Uriarte A, Malangone C, et al. Open targets platform:supporting systematic drug-target identification and prioritisation. Nucleic Acids Res. (2021) 49(D1):D1302–10. doi: 10.1093/nar/gkaa1027
35. Tomaiuolo F, MacDonald JD, Caramanos Z, Posner G, Chiavaras M, Evans AC, et al. Morphology, morphometry and probability mapping of the pars opercularis of the inferior frontal gyrus:an in vivo MRI analysis. Eur J Neurosci. (1999) 11(9):3033–46. doi: 10.1046/j.1460-9568.1999.00718.x
36. Zhang H, Ma L, Peng W, Wang B, Sun Y. Association between gut microbiota and onset of type 2 diabetes mellitus:a two-sample Mendelian randomization study. Front Cell Infect Microbiol. (2024) 14:1327032. doi: 10.3389/fcimb.2024.1327032
37. Foundas AL, Eure KF, Luevano LF, Weinberger DR. MRI asymmetries of broca’s area:the pars triangularis and pars opercularis. Brain Lang. (1998) 64(3):282–96. doi: 10.1006/brln.1998.1974
38. Kühn S, Brass M, Gallinat J. Imitation and speech:commonalities within broca’s area. Brain Struct Funct. (2013) 218(6):1419–27. doi: 10.1007/s00429-012-0467-5
39. Rakic P. Specification of cerebral cortical areas. Science. (1988) 241(4862):170–6. doi: 10.1126/science.3291116
40. Ozzoude M, Varriano B, Beaton D, Ramirez J, Adamo S, Holmes MF, et al. White matter hyperintensities and smaller cortical thickness are associated with neuropsychiatric symptoms in neurodegenerative and cerebrovascular diseases. Alzheimers Res Ther. (2023) 15(1):114. doi: 10.1186/s13195-023-01257-y
41. Fjell AM, Grydeland H, Krogsrud SK, Amlien I, Rohani DA, Ferschmann L, et al. Development and aging of cortical thickness correspond to genetic organization patterns. Proc Natl Acad Sci U S A. (2015) 112(50):15462–7. doi: 10.1073/pnas.1508831112
42. Holland MA, Budday S, Li G, Shen D, Goriely A, Kuhl E. Folding drives cortical thickness variations. Eur Phys J Spec Top. (2020) 229(17–18):2757–78. doi: 10.1140/epjst/e2020-000001-6
43. Frangou S, Modabbernia A, Williams SCR, Papachristou E, Doucet GE, Agartz I, et al. Cortical thickness across the lifespan:data from 17,075 healthy individuals aged 3–90 years. Hum Brain Mapp. (2022) 43(1):431–51. doi: 10.1002/hbm.25364
44. Sampedro F, Puig-Davi A, Martinez-Horta S, Pagonabarraga J, Horta-Barba A, Aracil-Bolaños I, et al. Cortical macro and microstructural correlates of cognitive and neuropsychiatric symptoms in Parkinson’s disease. Clin Neurol Neurosurg. (2023) 224:107531. doi: 10.1016/j.clineuro.2022.107531
45. Yang T, Lu Z, Wang L, Zhao Y, Nie B, Xu Q, et al. Dynamic changes in brain glucose metabolism and neuronal structure in rats with heart failure. Neuroscience. (2020) 424:34–44. doi: 10.1016/j.neuroscience.2019.10.008
46. Roy B, Woo MA, Wang DJJ, Fonarow GC, Harper RM, Kumar R. Reduced regional cerebral blood flow in patients with heart failure. Eur J Heart Fail. (2017) 19(10):1294–302. doi: 10.1002/ejhf.874
47. Moazzami K, Wittbrodt MT, Lima BB, Nye JA, Mehta PK, Pearce BD, et al. Higher activation of the rostromedial prefrontal cortex during mental stress predicts major cardiovascular disease events in individuals with coronary artery disease. Circulation. (2020) 142(5):455–65. doi: 10.1161/CIRCULATIONAHA.119.044442
48. Lu Z, Teng Y, Wang L, Jiang Y, Li T, Chen S, et al. Abnormalities of hippocampus and frontal lobes in heart failure patients and animal models with cognitive impairment or depression:a systematic review. PLoS One. (2022) 17(12):e0278398. doi: 10.1371/journal.pone.0278398
49. Ichijo Y, Kono S, Yoshihisa A, Misaka T, Kaneshiro T, Oikawa M, et al. Impaired frontal brain activity in patients with heart failure assessed by near-infrared spectroscopy. J Am Heart Assoc. (2020) 9(3):e014564. doi: 10.1161/JAHA.119.014564
50. Kumar R, Yadav SK, Palomares JA, Park B, Joshi SH, Ogren JA, et al. Reduced regional brain cortical thickness in patients with heart failure. PLoS One. (2015) 10(5):e0126595. doi: 10.1371/journal.pone.0126595
51. Mueller K, Thiel F, Beutner F, Teren A, Frisch S, Ballarini T, et al. Brain damage with heart failure:cardiac biomarker alterations and gray matter decline. Circ Res. (2020) 126(6):750–64. doi: 10.1161/circresaha.119.315813
52. Baril AA, Gagnon K, Brayet P, Montplaisir J, De Beaumont L, Carrier J, et al. Gray matter hypertrophy and thickening with obstructive sleep apnea in middle-aged and older adults. Am J Respir Crit Care Med. (2017) 195(11):1509–18. doi: 10.1164/rccm.201606-1271OC
53. Sies H. Glutathione and its role in cellular functions. Free Radic Biol Med. (1999) 27(9–10):916–21. doi: 10.1016/s0891-5849(99)00177-x
54. Wu G, Fang YZ, Yang S, Lupton JR, Turner ND. Glutathione metabolism and its implications for health. J Nutr. (2004) 134(3):489–92. doi: 10.1093/jn/134.3.489
55. Lapenna D. Glutathione and glutathione-dependent enzymes:from biochemistry to gerontology and successful aging. Ageing Res Rev. (2023) 92:102066. doi: 10.1016/j.arr.2023.102066
56. Tan M, Yin Y, Ma X, Zhang J, Pan W, Tan M, et al. Glutathione system enhancement for cardiac protection:pharmacological options against oxidative stress and ferroptosis. Cell Death Dis. (2023) 14(2):131. doi: 10.1038/s41419-023-05645-y
57. Averill-Bates DA. The antioxidant glutathione. Vitam Horm. (2023) 121:109–41. doi: 10.1016/bs.vh.2022.09.002
58. Georgiou-Siafis SK, Tsiftsoglou AS. The key role of GSH in keeping the redox balance in mammalian cells:mechanisms and significance of GSH in detoxification via formation of conjugates. Antioxidants (Basel). (2023) 12(11):1953. doi: 10.3390/antiox12111953
59. Rashdan NA, Shrestha B, Pattillo CB. S-glutathionylation, friend or foe in cardiovascular health and disease. Redox Biol. (2020) 37:101693. doi: 10.1016/j.redox.2020.101693
60. Bertero E, Maack C. Ins and outs of glutathione in cardiac ischemia/reperfusion injury. Circ Res. (2023) 133(10):877–9. doi: 10.1161/circresaha.123.323715
61. Burns M, Rizvi SHM, Tsukahara Y, Pimentel DR, Luptak I, Hamburg NM, et al. Role of glutaredoxin-1 and glutathionylation in cardiovascular diseases. Int J Mol Sci. (2020) 21(18):6803. doi: 10.3390/ijms21186803
63. Asanuma M, Miyazaki I. Glutathione and related molecules in parkinsonism. Int J Mol Sci. (2021) 22(16):8689. doi: 10.3390/ijms22168689
64. Haddad M, Hervé V, Ben Khedher MR, Rabanel JM, Ramassamy C. Glutathione:an old and small molecule with great functions and new applications in the brain and in Alzheimer’s disease. Antioxid Redox Signal. (2021) 35(4):270–92. doi: 10.1089/ars.2020.8129
Keywords: chronic heart failure, cerebral cortical structure, Mendelian randomization, causal association, integrated bioinformatics
Citation: Peng L, Cai H, Tang Y, Zhou F, Liu Y, Xu Z, Chen Q and Chen X (2024) Causal associations between chronic heart failure and the cerebral cortex: results from Mendelian randomization study and integrated bioinformatics analysis. Front. Cardiovasc. Med. 11:1396311. doi: 10.3389/fcvm.2024.1396311
Received: 5 March 2024; Accepted: 11 June 2024;
Published: 4 July 2024.
Edited by:
Aifeng Zhang, Boston Medical Center, United StatesReviewed by:
Marisol Herrera Rivero, University of Münster, GermanyYanguo Xin, Capital Medical University, China
Si Cao, Reproductive and Genetic Hospital of CITIC-Xiangya, China
© 2024 Peng, Cai, Tang, Zhou, Liu, Xu, Chen and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingyang Chen, NjIzNDMzNzQ0QHFxLmNvbQ==; Xinyu Chen, Y2hlbnhpbnl1Y2hlbkAxNjMuY29t