Skip to main content

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 01 December 2021
Sec. General Cardiovascular Medicine
This article is part of the Research Topic Cardiovascular Risk Factors: Related vascular injury and new molecular biomarkers View all 23 articles

Identification of CALU and PALLD as Potential Biomarkers Associated With Immune Infiltration in Heart Failure

\nXing Liu&#x;Xing Liu1Shiyue Xu&#x;Shiyue Xu2Ying Li&#x;Ying Li3Qian ChenQian Chen1Yuanyuan Zhang
Yuanyuan Zhang1*Long Peng
Long Peng1*
  • 1Department of Cardiovascular Medicine, The Third Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
  • 2Department of Hypertension and Vascular Disease, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
  • 3Department of Dermatology, Guangzhou Eighth People's Hospital, Guangzhou Medical University, Guangzhou, China

Background: Inflammatory activation and immune infiltration play important roles in the pathologic process of heart failure (HF). The current study is designed to investigate the immune infiltration and identify related biomarkers in heart failure patients due to ischemic cardiomyopathy.

Methods: Expression data of HF due to ischemic cardiomyopathy (CM) samples and non-heart failure (NF) samples were downloaded from gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) between CM and NF samples were identified. Single sample gene set enrichment analysis (ssGSEA) was performed to explore the landscape of immune infiltration. Weighted gene co-expression network analysis (WGCNA) was applied to screen the most relevant module associated with immune infiltration. The diagnostic values of candidate genes were evaluated by receiver operating curves (ROC) curves. The mRNA levels of potential biomarkers in the peripheral blood mononuclear cells (PBMCs) isolated from 10 CM patients and 10 NF patients were analyzed to further assess their diagnostic values.

Results: A total of 224 DEGs were identified between CM and NF samples in GSE5406, which are mainly enriched in the protein processing and extracellular matrix related biological processes and pathways. The result of ssGSEA showed that the abundance of dendritic cells (DC), mast cells, natural killer (NK) CD56dim cells, T cells, T follicular helper cells (Tfh), gammadelta T cells (Tgd) and T helper 2 (Th2) cells were significantly higher, while the infiltration of eosinophils and central memory T cells (Tcm) were lower in CM samples compared to NF ones. Correlation analysis revealed that Calumenin (CALU) and palladin (PALLD) were negatively correlated with the abundance of DC, NK CD56dim cells, T cells, Tfh, Tgd and Th2 cells, but positively correlated with the level of Tcm. More importantly, CALU and PALLD were significantly lower in PBMCs from CM patients compared to NF ones.

Conclusion: Our study revealed that CALU and PALLD are potential biomarkers associated with immune infiltration in heart failure due to ischemic cardiomyopathy.

Introduction

Heart failure refers to the complex clinical syndrome and the end-stage manifestations of cardiovascular disease (1, 2). Large-scale epidemiological analysis shows that the global prevalence of HF is on the rise due to the aging of the population and the progress in the diagnosis and treatment of cardiovascular diseases. In developed countries, the prevalence of HF is 1.5–2.2% (3). The latest report in China shows that the prevalence of HF among residents ≥35 years old is 1.3%, that is, there are ~13.7 million patients (3).Therefore, HF has always been a hot spot in the field of cardiovascular research.

HF is caused by the complicated interaction of myocardial damage, neurohormonal activation, inflammatory response, and renal dysfunction (46). During the process, the cytoskeletal and membrane associated proteins are increased and disorganized, while the contractile myofilaments and sarcomeric proteins are decreased in the heart (7). In addition, cardiomyocytes in the failing heart display impaired excitation-contraction coupling due to the decreased calcium transients, enhanced diastolic sarcoplasmic reticulum (SR) calcium leak and diminished SR calcium sequestration (8). Although the pathogenesis of HF is still perplexing, the persistent inflammation and immune abnormalities are believed to participate in the pathogenesis across the spectrum of HF (9). Elevated and long-lasting leukocyte recruitment mediated by G protein-coupled receptor kinase 5 (GRK5) in the injured heart is reported to be associated with chronic cardiac inflammation and heart failure (10). Moreover, evidence indicates that transcriptome changes in immune cells could affect the prognosis of HF. DNA methyltransferase 3 alpha (DNMT3A) mutations in monocytes significantly increase the expression of inflammatory genes and are correlated with the aggravation of chronic HF (11). Metabolically active genes such as fatty acid binding protein 5 (FABP5) are highly enriched in classical monocytes from heart failure patients, whereas b-catenin expression was significantly higher in another functionally distinct monocyte subset (CD14++CD16+ intermediate monocytes) (12). These studies suggest that further understanding of the inflammatory response and immune cell infiltration in HF is of great significance for optimizing the diagnosis and treatment of heart failure.

In recent years, microarray technology and integrated bioinformatics analyses have been performed to identify novel genes related to various diseases that might act as diagnostic and prognostic biomarkers (1315). However, the diagnostic value of genes associated with immune infiltration in heart failure remains unclear. Thus, in the current study, we downloaded two microarray datasets of HF from the GEO database and used bioinformatic methods to screen for immune infiltration related biomarkers in heart failure, and to provide a theoretical basis for the diagnosis and treatment of HF patients.

Materials and Methods

Ethics Statement

This study was approved by the institutional review board of the Third Affiliated Hospital, Sun Yat-sen University (IRB: 202102-201-01).

Data Source

In the current study, gene expression profiles of 16 non-failure controls (NF) and 108 heart failure samples caused by ischemic cardiomyopathy (CM) in GSE5406 dataset and 14 NF and 13 CM samples in GSE116250 dataset were downloaded from GEO database.

Identification and Functional Enrichment Analysis of DEGs

Limma R package was used to identify DEGs between NF and CM samples with |log2FC| >0.5 and adjusted p < 0.05 in GSE5406 datasets. ClusterProfiler R package was applied for GO and KEGG pathway enrichment analyses of DEGs. Biological process (BP), molecular function (MF) and cellular component (CC) were included in the GO analysis.

ssGSEA

ssGSEA was performed by Gene Set Variation Analysis (GSVA) R package to analyze the infiltraion of 24 immune cells in NF and CM samples (16). The 24 immune cells were TFH, Th2 cells, B cells, T cells, Tgd, NK CD56dim cells, Tem, macrophages, neutrophils, Th1 cells, mast cells, cytotoxic cells, DC, iDC, eosinophils, T helper cells, aDC, TReg, pDC, NK CD56bright cells, NK cells, Th17 cells, CD8 T cells, and Tcm.

WGCNA Analysis

A sample clustering tree map was first constructed to detect and eliminate outliers. Then, WGCNA was performed based on the gene expression profiles from GSE5406 dataset and sample traits (differentially infiltrated immune cells between NF and CM samples). The pick Soft Threshold function of WGCNA was used to calculate β from 1 to 20 in order to select the best soft threshold. Based on the selected soft threshold, the adjacency matrix was converted to topological overlap matrix to construct the network, and the gene dendrogram and module color were established by using the degree of dissimilarity. We further divided the initial module by dynamic tree cutting and merged similar modules. The Pearson correlation coefficient between the module eigengenes and sample traits were calculated to find out the most relevant module (hub module) associated with sample traits.

Identification of Biomarkers in CM

First, DEGs were intersected with genes from the hub module in WGCNA analysis to obtain immune infiltration related candidate genes. Next, gene signature was selected by least absolute shrinkage and selection operator (LASSO) algorithm using glmnet R package (17) and support vector machine-recursive feature elimination (SVM-RFE) method using e1071 package (18), respectively. Robust gene signature was identified by overlapping gene signature obtained from LASSO and SVM-RFE. The diagnostic values of gene signature were evaluated by receiver operating curves (ROC) curves. Then, the external validation dataset GSE116250 was used to verify the expressions and diagnostic values of gene signature identified in GSE5406. Validated gene signature was identified as robust diagnostic biomarkers in heart failure.

Functional Analysis of Biomarkers in CM

To investigate the potential mechanisms of diagnostic biomarkers in regulating heart failure, 108 patients in GSE5406 were divided into high- and low-expression groups based on the median expression of each diagnostic biomarker. Moreover, to explore the relationship between diagnostic biomarkers and immune infiltration, the correlations between the expressions of diagnostic biomarkers and the abundance of differentially infiltrated immune cells were calculated.

Subject Characteristics and Realtime-PCR

Patients aged 18 and older, diagnosed with CHD by coronary computed tomography angiography or coronary angiography, with ejection fraction of 40% or less were enrolled into CM group. Age-matched CHD patients without heart failure (ejection fraction of 50% or above) were enrolled into NF control group. Patients with a history of malignancy, acute coronary syndrome, pulmonary embolism, renal failure [Glomerular filtration rate <60 ml/(min·1.73 m2)] were excluded. The characteristics of CM and NF patients were shown in Supplementary Table 3.

RNA of PBMCs from CM (n = 10, 7 male and 3 female) and NF (n = 10, 8 male and 2 female) patients were extracted using Nuclezol LS RNA Isolation Reagent (ABP Biosciences Inc.) according to manufacturer's instructions. Collected RNA was diluted using nuclease-free water and electrophoresed on a denaturing formaldehyde agarose gel to visualize rRNA and ensure overall sample quality. RNA concentrations and purity were detected on ultraviolet spectrophotometer (Jinghua, Shanghai, China). Reverse transcription was performed on 1 μg total RNA from each sample using the SureScript-First-strand-cDNA-synthesis-kit (GeneCopoeia) according to manufacturer's instructions, and a CFX96 Real-time PCR System (Bio-Rad) was utilized to conduct the real time quantitative PCR (qPCR) reactions. BlazeTaq™ SYBR® Green qPCR Mix 2.0 (GeneCopoeia) was used for qPCR reactions, using 4 μL cDNA and appropriate volumes of specific primers in a final 10 μL volume. Triplicate reactions were performed to ensure accuracy. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as the reference gene, and the relative gene expression was quantified by the 2−ΔΔCT method (19). The primer sequences were given in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Primers sequence.

Statistical Analysis

All data were analyzed by R software (version 4.0.0). Wilcoxon test was used to compare the data between two groups, and significant difference was considered as p < 0.05 unless specified.

Results

Transcriptome Profile Analyses of NF and CM Samples

A total of 224 DEGs were identified in GSE5406, including 93 up-regulated and 131 down-regulated genes in CM group compared to NF group (Figure 1A). The expression profile of top 50 up-regulated DEGs and top 50 down-regulated DEGs were shown in the heatmap (Figure 1B). To investigate the biological function of DEGs, we performed GO and KEGG pathway analysis. A total of 217 BP, 42 CC, 37 MF, and 16 KEGG pathways were significantly enriched (Supplementary Tables 1, 2). As shown in Figure 1C, DEGs were mainly enriched into protein processing and extracellular matrix (ECM) related BPs, including response to topologically incorrect protein, response to unfolded protein (UPR), “de novo” protein folding, protein folding, chaperone-mediated protein folding, extracellular matrix organization, extracellular structure organization, “de novo” posttranslational protein folding, chaperone cofactor-dependent protein refolding, response to mechanical stimulus. Consistent with the results of GO analysis, protein processing and ECM related pathways were significantly enriched, including protein processing in endoplasmic reticulum, ECM-receptor interaction, and focal adhesion. In addition, estrogen signaling and MAPK signaling pathways showed to have close relationship with heart failure (Figure 1D).

FIGURE 1
www.frontiersin.org

Figure 1. Comprehensive analyses of the transcriptome profiles of NF and CM samples. (A) Volcano plot of significant DEGs between NF and CM samples. (B) A heatmap of the top 50 significantly upregulated or downregulated DEGs. (C) Bar plot of top 10 enriched GO terms of DEGs in each category. BP, biological process; CC, cellular components; MF, molecular functions. (D) Bubble plot of significantly enriched KEGG pathways of DEGs.

Identification of Immune Infiltration Pattern in CM

Mounting evidence suggest that immune cells play important roles in heart failure (1416). Thus, we explored the profile of immune cell infiltration in CM and NF samples by ssGSEA. Twenty-four subpopulations of infiltrated immune cells in CM and NF samples were identified and presented in the heatmap (Figure 2A). Interestingly, we found that the abundance of DC, mast cells, NK CD56dim cells, T cells, Tfh, Tgd, and Th2 cells were significantly higher, while the infiltration of eosinophils and Tcm were significantly lower in CM samples compared to NF ones (Figure 2B). These results indicate that the inflammatory response of these immune cells may be critical for the etiology of heart failure.

FIGURE 2
www.frontiersin.org

Figure 2. Landscape of the immune infiltration in NF and CM samples. (A) Heatmap of the immune infiltration profiles of NF and CM samples analyzed by ssGSEA score-based method. (B) Comparison of immune cell infiltration between NF and CM samples. *p < 0.05, **p < 0.01, ****p < 0.001.

Screening for Gene Signature of Immune Infiltration in CM

To further explore the genes mostly correlated with the inflammatory response in CM, we performed WGCNA to screen for the hub module associated with above infiltrated immune cells. After eliminating the outlier samples (Supplementary Figure 1), we built the sample dendrogram and trait heatmap (Supplementary Figure 2). By using the pick Soft Threshold function of WGCNA, we found the optimal soft threshold power was 9, in which R2 was 0.85 (Supplementary Figure 3A). After merging similar modules, eight modules from the co-expression network were identified (Supplementary Figure 3B). According to the module-trait relationships in Figure 3, we found that the MEbrown module was the most relevant module associated with DC (Cor = 0.3, p < 0.01), NK CD56dim cells (Cor = 0.62, p < 0.01), T cells (Cor = 0.81, p < 0.01), Tcm (Cor = −0.5, p < 0.01), TFH (Cor = 0.52, p < 0.01), Tgd (Cor = 0.65, p < 0.01) and Th2 cells (Cor =0.49, p < 0.01). Thus, MEbrown module was selected for downstream analysis.

FIGURE 3
www.frontiersin.org

Figure 3. Screening of the hub module associated with immune infiltration in CM. Heatmap of the correlation between module eigengene and differentially infiltrated immune cells. Each row represents a color-coded module eigengene, each column represents a type of infiltrated immune cells. The number in each cell means the correlation coefficient and p-value.

Next, we overlapped DEGs with genes in the MEbrown module and obtained 10 candidate genes (Figure 4A). Ten candidate genes were input into LASSO and SVM-RFE to identify gene signature, respectively. LASSO identified nine gene signatures under lambda.min = 0.0021, including PALLD, DexD-box helicase 39A (DDX39A), stress induced phosphoprotein 1 (STIP1), solute carrier family 38 member 2 (SLC38A2), CALU, CD164 molecule (CD164), selenoprotein T (SELT), four and a half LIM domain 1 (FHL1) and claudin domain containing 1 (CLDND1) (Figures 4B,C). The accuracy of LASSO was evaluated by ROC curve that the area under the ROC curve (AUC) was 0.902 (Figure 4D). Meanwhile, we identified 7 gene signatures by SVM-RFE, including PALLD, DDX39A, CD164, CLDND1, SLC38A2, CALU, and heat shock protein family D member 1 (HSPD1) with the accuracy of 0.962 (Figures 5A,B). To get the robust gene signature in heart failure, we overlapped genes from LASSO and SVM-RFE and got six gene signatures (Figure 5C), including PALLD, DDX39A, CD164, CLDND1, SLC38A2, and CALU. The expression levels of PALLD, DDX39A, CD164, CLDND1, SLC38A2, and CALU were all significantly higher in NF samples compared to CM ones in GSE5406 (Figure 5D).

FIGURE 4
www.frontiersin.org

Figure 4. Identification of candidate genes by LASSO regression model. (A) Venn diagram of 10 overlapped candidate genes shared by DEGs and MEbrown module. (B) LASSO coefficient profiles of candidate genes. (C) Cross-validation to select the optimal tuning parameter log (Lambda) in LASSO regression analysis. (D) ROC curve evaluation of LASSO regression analysis.

FIGURE 5
www.frontiersin.org

Figure 5. Identification of candidate genes by SVM-RFE. (A) 7 gene signatures are identified by SVM-RFE analysis with the accuracy of 0.962 and (B) error of 0.038. (C) Venn diagram of six overlapped candidate genes shared by the LASSO and SVM-RFE algorithms. (D) The expressions of candidate diagnostic biomarkers in the NF and CM samples from GSE5406 dataset.

Verification of the Potential Biomarkers for CM

We further evaluated the diagnostic values of PALLD, DDX39A, CD164, CLDND1, SLC38A2, and CALU in GSE5406 by ROC curves. We found that they all had high accuracy with AUC >0.7 (Figure 6A). To verify the diagnostic values of the candidate biomarkers, we examined their expressions in an external validation dataset GSE116250 to get robust diagnostic biomarkers. We found that the expression trends of CALU and PALLD in GSE116250 were consistent with those in GSE5406 (Figure 6B), and that CALU and PALLD also had high accuracy in classifying heart failure samples, as evidenced by AUC >0.7 (Figure 6C). Moreover, we found that the expressions of CALU and PALLD were significantly negative correlated with DC, NK CD56dim cells, T cells, Tfh, Tgd, and Th2 cells, and positive correlated with Tcm (Figure 6D).

FIGURE 6
www.frontiersin.org

Figure 6. Verification of biomarkers for CM. (A) ROC curve evaluation of the diagnostic effectiveness of candidate biomarkers using GSE5406 dataset. (B) The expressions of candidate diagnostic biomarkers in the GSE116250 dataset. (C) ROC curve evaluation of the diagnostic effectiveness of candidate biomarkers using GSE116250 dataset. (D) Heatmap of correlations between PALLD, CALU and differentially infiltrated immune cells. (E) Real-time PCR analyses of the expression levels of PALLD and CALU in PBMCs isolated from CM and NF patients *p < 0.05, **p < 0.01, ***p < 0.001

To further evaluate the value of CALU and PALLD as biomarkers, the levels of CALU and PALLD were assessed in PBMCs isolated from CM and NF patients. In agreement with the results in GSE116250 and GSE5406, the levels of CALU and PALLD were also significantly lower in PBMCs collected from CM patients compared to NF ones (Figure 6E) Collectively, these results suggest that CALU and PALLD may be potential diagnostic biomarkers for heart failure due to ischemic cardiomyopathy.

Discussion

The prevalence of HF is increasing worldwide and has reached epidemic proportions. Although significant progress has been made in the medications and interventions for HF, the mortality and hospitalization rates remain high (20). Recently, immune cell infiltration has been confirmed to play a vital role in the occurrence and development of cardiovascular diseases (10, 13). Cardiac inflammation and subsequent tissue damage are orchestrated by the infiltration and activation of various immune cells in the myocardium which result in heart failure eventually (20). Thus, it is essential to comprehensively investigate the contributions of infiltrated immune cells in HF. In this study, we analyzed the immune infiltration profiles of CM and NF samples, and identified CALU and PALLD as potential diagnostic biomarkers for heart failure due to ischemic cardiomyopathy by using integrated bioinformatics analyses.

Consistent with a previous report (21), we found that DEGs between NF and CM were mainly enriched into biological processes and pathways related to protein processing and ECM. Structural changes occur in the level of ER and UPR components in cardiomyocytes of patients with HF (22). When maladaptive UPR fails to restore ER homeostasis, it might induce risk factors for HF including increased reactive oxygen species (ROS) production, inflammation and apoptosis which further aggravate HF (22). In addition, numerous studies have explored the relationship between ECM and heart diseases (23, 24). Alterations in the architecture, composition, and distribution of interstitial ECM play a major role in pathological myocardial structural remodeling and left ventricular diastolic dysfunction (25, 26). Myocardial fibrosis is characterized by accumulation of collagen-rich ECM, such as collagen type I and III fibers, results from the predominance of fiber formation and deposition over its degradation and removal (27). ECM dyshomeostasis is also postulated to occur during the development and progression of HF, including changes in the synthesis, processing, degradation, and turnover of proteins such as fibrillar collagen (28). Therefore, those DEGs we found may regulate HF through protein processing and ECM.

Increasing evidence show that immune cell infiltration in the myocardium has adverse effect on heart function (2931). Single-cell sequencing analyses reveal that the immune cell profiles are remarkably different in healthy and diseased hearts (15, 32, 33). In this study, we found that the abundance of DC, mast cells, NK CD56dim cells, T cells, Tfh, Tgd, and Th2 cells were higher, while the infiltration of eosinophils and Tcm were lower in HF samples, indicating their important roles in the etiology of HF. Consistent with our findings, Patella V et al. and Abdolmaleki F et al. found an increase in numbers of mast cells and T cells in HF, respectively (34, 35). Mast cells initiate adverse myocardial remodeling by activating matrix metalloproteinase (MMP) and fibrosis in the heart (36). Moreover, the profibrotic and antiangiogenic functions of Th17, Th2, and dysfunctional Treg cells are indispensable for the progression to ischemic heart failure (37, 38). On the contrary, inducible depletion of eosinophils exacerbates cardiac dysfunction, cell death, and fibrosis, by producing IL-4 and cationic protein mEar1, fibroblast activation, and neutrophil adhesion (39). Studies from our and other groups suggest a critical role of immune infiltration in the development of heart failure, and further understandings of the effect of each type of infiltrated immune cells may provide us clues for developing novel therapeutic strategies for heart failure.

PALLD encodes a cytoskeleton protein involved in actin reorganization (40), which plays an important role in heart development (41). PALLD was reported to be related to vein graft stenosis after coronary artery bypass grafting. A single-nucleotide polymorphism (SNP) in the PALLD has been reported to be associated with coronary heart disease (CHD) (42). CALU produces a Ca2+-binding protein that is localized in the endoplasmic reticulum (ER) (43). During the excitation-contraction coupling process, CALU regulates Ca2+ uptake and plays an important role in maintaining normal heart function (44). Our analyses show that the expressions of CALU and PALLD were significantly negative correlated with the abundances of DC, NK CD56dim cells, T cells, Tfh, Tgd and Th2 cells, and Tcm, suggesting that CALU and PALLD may regulate HF via immune-related pathways mediated by the infiltration of these immune cells. More importantly, we confirmed that the expression levels of CALU and PALLD were markedly lower in PBMCs from CM patients compared with NF ones. Collectively, these findings indicate that CALU and PALLD are potential biomarkers for the diagnosis of HF.

Several limitations of the present study should be noted. Firstly, the study was retrospective, further prospective studies are needed to assess the diagnostic and prognostic value of CALU and PALLD in HF. Secondly, the relationships between the biomarkers and immune regulation in heart failure were only verified by assessing their levels in PBMCs from CM and NF patients. Further in vitro and in vivo experiments are required to explore the detailed mechanisms through which CALU and PALLD regulate inflammatory responses in HF. Taken together, our analyses delineate the potential etiology of HF due to CM and identified immune infiltration related biomarkers in HF, which may provide guidance for diagnosis and treatment of patients with HF.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Review Board of the Third Affiliated Hospital, Sun Yat-sen University (IRB: 202102-201-01). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

YZ and LP contributed to the study concepts and study design and helped in revising the manuscript. XL, SX, and YL drafted the manuscript, performed data management, and bioinformatics analysis. QC were responsible for clinical sample collection. All authors were involved in reporting the results of this study and approved the final version of the submitted manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (82000466), the Medical Science and Technology Research Project of Guangzhou (Grant No. 202002020030), and the cultivation project of National Natural Science Foundation of the Third Affiliated Hospital, Sun Yat-sen University (Grant No. 2021G2RPYQN11).

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.2021.774755/full#supplementary-material

References

1. Ponikowski P, Voors AA, Anker SD, Bueno H, Cleland JGF, Coats AJS, et al. 2016 ESC Guidelines for the diagnosis treatment of acute chronic heart failure: The Task Force for the diagnosis treatment of acute chronic heart failure of the European Society of Cardiology (ESC)Developed with the special contribution of the Heart Failure Association (HFA) of the ESC. Eur Heart J. (2016) 37:2129–200. doi: 10.1093/eurheartj/ehw128

CrossRef Full Text | Google Scholar

2. Yancy CW, Jessup M, Bozkurt B, Butler J, Casey DE, Colvin MM, et al. 2017 ACC/AHA/HFSA Focused Update of the 2013 ACCF/AHA guideline for the management of heart failure: a report of the American College of Cardiology/American heart association task force on clinical practice guidelines and the heart failure society of America. Circulation. (2017) 136:e137–61. doi: 10.1161/CIR.0000000000000509

CrossRef Full Text | Google Scholar

3. Cunningham JW, Claggett BL, O'Meara E, Prescott MF, Pfeffer MA, Shah SJ, et al. Effect of sacubitril/valsartan on biomarkers of extracellular matrix regulation in patients with HFpEF. J Am Coll Cardiol. (2020) 76:503–14. doi: 10.1016/j.jacc.2020.05.072

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Mann DL. Inflammatory mediators and the failing heart: past, present, and the foreseeable future. Circ Res. (2002) 91:988–98. doi: 10.1161/01.RES.0000043825.01705.1B

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zimmet JM, Hare JM. Nitroso-redox interactions in the cardiovascular system. Circulation. (2006) 114:1531–44. doi: 10.1161/CIRCULATIONAHA.105.605519

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Braunwald E. Biomarkers in heart failure. N Engl J Med. (2008) 358:2148–59. doi: 10.1056/NEJMra0800239

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Hein S, Kostin S, Heling A, Maeno Y, Schaper J. The role of the cytoskeleton in heart failure. Cardiovasc Res. (2000) 45:273–8. doi: 10.1016/S0008-6363(99)00268-0

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Luo M, Anderson ME. Mechanisms of altered Ca(2)(+) handling in heart failure. Circ Res. (2013) 113:690–708. doi: 10.1161/CIRCRESAHA.113.301651

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Dick SA, Epelman S. Chronic heart failure and inflammation what do we really know? Circul Res. (2016) 119:159–76. doi: 10.1161/CIRCRESAHA.116.308030

PubMed Abstract | CrossRef Full Text | Google Scholar

10. de Lucia C, Grisanti LA, Borghetti G, Piedepalumbo M, Ibetti J, Maria Lucchese A, et al. GRK5 contributes to impaired cardiac function and immune cell recruitment in post-ischemic heart failure. Cardiovasc Res. (2021) 9:cvab044. doi: 10.1093/cvr/cvab044

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Abplanalp WT, Cremer S, John D, Hoffmann J, Schuhmacher B, Merten M, et al. Clonal hematopoiesis-driver DNMT3A mutations alter immune cells in heart failure. Circ Res. (2021) 128:216–28. doi: 10.1161/CIRCRESAHA.120.317104

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Abplanalp WT, John D, Cremer S, Assmus B, Dorsheimer L, Hoffmann J, et al. Single-cell RNA-sequencing reveals profound changes in circulating immune cells in patients with heart failure. Cardiovasc Res. (2021) 117:484–94. doi: 10.1093/cvr/cvaa101

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zhao E, Xie H, Zhang Y. Predicting diagnostic gene biomarkers associated with immune infiltration in patients with acute myocardial infarction. Front Cardiovasc Med. (2020) 7:586871. doi: 10.3389/fcvm.2020.586871

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Cao Y, Tang W, Tang W. Immune cell infiltration characteristics and related core genes in lupus nephritis: results from bioinformatic analysis. BMC Immunol. (2019) 20:37. doi: 10.1186/s12865-019-0316-x

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Martini E, Kunderfranco P, Peano C, Carullo P, Cremonesi M, Schorn T, et al. Single-cell sequencing of mouse heart immune infiltrate in pressure overload-driven heart failure reveals extent of immune activation. Circulation. (2019) 140:2089–107. doi: 10.1161/CIRCULATIONAHA.119.041694

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. (2013) 14:7. doi: 10.1186/1471-2105-14-7

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Vasquez MM, Hu CC, Roe DJ, Chen Z, Halonen M, Guerra S. Least absolute shrinkage and selection operator type methods for the identification of serum biomarkers of overweight and obesity: simulation and application. Bmc Med Res Methodol. (2016) 16:254. doi: 10.1186/s12874-016-0254-8

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Noble WS. What is a support vector machine? Nat Biotechnol. (2006) 24:1565–7. doi: 10.1038/nbt1206-1565

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C-T method. Nat Protoc. (2008) 3:1101–8. doi: 10.1038/nprot.2008.73

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Strassheim D, Dempsey EC, Gerasimovskaya E, Stenmark K, Karoor V. Role of inflammatory cell subtypes in heart failure. J Immunol Res. (2019) 2019:2164017. doi: 10.1155/2019/2164017

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ren J, Bi YG, Sowers JR, Hetz C, Zhang YM. Endoplasmic reticulum stress and unfolded protein response in cardiovascular diseases. Nat Rev Cardiol. (2021) 18:499–521. doi: 10.1038/s41569-021-00511-w

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Glembotski CC. Endoplasmic reticulum stress in the heart. Circ Res. (2007) 101:975–84. doi: 10.1161/CIRCRESAHA.107.161273

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Frangogiannis NG. The extracellular matrix in ischemic and nonischemic heart failure. Circul Res. (2019) 125:117–46. doi: 10.1161/CIRCRESAHA.119.311148

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Guo Y, Gupte M, Umbarkar P, Singh AP, Sui JY, Force T, et al. Entanglement of GSK-3beta, beta-catenin and TGF-beta1 signaling network to regulate myocardial fibrosis. J Mol Cell Cardiol. (2017) 110:109–20. doi: 10.1016/j.yjmcc.2017.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kasner M, Westermann D, Lopez B, Gaub R, Escher F, Kuhl U, et al. Diastolic tissue Doppler indexes correlate with the degree of collagen expression and cross-linking in heart failure and normal ejection fraction. J Am Coll Cardiol. (2011) 57:977–85. doi: 10.1016/j.jacc.2010.10.024

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zile MR, Baicu CF, Ikonomidis JS, Stroud RE, Nietert PJ, Bradshaw AD, et al. Myocardial stiffness in patients with heart failure and a preserved ejection fraction: contributions of collagen and titin. Circulation. (2015) 131:1247–59. doi: 10.1161/CIRCULATIONAHA.114.013215

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lopez B, Ravassa S, Gonzalez A, Zubillaga E, Bonavila C, Berges M, et al. Myocardial collagen cross-linking is associated with heart failure hospitalization in patients with hypertensive heart failure. J Am College Cardiol. (2016) 67:251–60. doi: 10.1016/j.jacc.2015.10.063

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Gonzalez A, Schelbert EB, Diez J, Butler J. Myocardial interstitial fibrosis in heart failure: biological and translational perspectives. J Am Coll Cardiol. (2018) 71:1696–706. doi: 10.1016/j.jacc.2018.02.021

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Carrillo-Salinas FJ, Ngwenyama N, Anastasiou M, Kaur K, Alcaide P. Heart inflammation immune cell roles and roads to the heart. Am J Pathol. (2019) 189:1482–94. doi: 10.1016/j.ajpath.2019.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Lafuse WP, Wozniak DJ, Rajaram MVS. Role of cardiac macrophages on cardiac inflammation, fibrosis and tissue repair. Cells. (2020) 10:51. doi: 10.3390/cells10010051

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Kawada J, Takeuchi S, Imai H, Okumura T, Horiba K, Suzuki T, et al. Immune cell infiltration landscapes in pediatric acute myocarditis analyzed by CIBERSORT. J Cardiol. (2021) 77:174–8. doi: 10.1016/j.jjcc.2020.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Farbehi N, Patrick R, Dorison A, Xaymardan M, Janbandhu V, Wystub-Lis K, et al. Single-cell expression profiling reveals dynamic flux of cardiac stromal, vascular and immune cells in health and injury. Elife. (2019) 8:43882. doi: 10.7554/eLife.43882

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Suryawanshi H, Clancy R, Morozov P, Halushka MK, Buyon JP, Tuschl T. Cell atlas of the foetal human heart and implications for autoimmune-mediated congenital heart block. Cardiovasc Res. (2020) 116:1446–57. doi: 10.1093/cvr/cvz257

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Patella V, Marino I, Arbustini E, Lamparter-Schummert B, Verga L, Adt M, et al. Stem cell factor in mast cells and increased mast cell density in idiopathic and ischemic cardiomyopathy. Circulation. (1998) 97:971–8. doi: 10.1161/01.CIR.97.10.971

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Abdolmaleki F, Gheibi Hayat SM, Bianconi V, Johnston TP, Sahebkar A. Atherosclerosis and immunity: A perspective. Trends Cardiovasc Med. (2019) 29:363–71. doi: 10.1016/j.tcm.2018.09.017

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Ngkelo A, Richart A, Kirk JA, Bonnin P, Vilar J, Lemitre M, et al. Mast cells regulate myofilament calcium sensitization and heart function after myocardial infarction. J Exp Med. (2016) 213:1353–74. doi: 10.1084/jem.20160081

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Bansal SS, Ismahil MA, Goel M, Patel B, Hamid T, Rokosh G, et al. Activated T lymphocytes are essential drivers of pathological remodeling in ischemic heart failure. Circ-Heart Fail. (2017) 10:3688. doi: 10.1161/CIRCHEARTFAILURE.116.003688

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Bansal SS, Ismahil MA, Goel M, Zhou GH, Rokosh G, Hamid T, et al. Dysfunctional and proinflammatory regulatory T-lymphocytes are essential for adverse cardiac remodeling in ischemic cardiomyopathy. Circulation. (2019) 139:206–21. doi: 10.1161/CIRCULATIONAHA.118.036065

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Liu J, Yang CZ, Liu TX, Deng ZY, Fang WQ, Zhang X, et al. Eosinophils improve cardiac function after myocardial infarction. Nat Commun. (2020) 11:5. doi: 10.1038/s41467-020-19297-5

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Mykkanen OM, Gronholm M, Ronty M, Lalowski M, Salmikangas P, Suila H, et al. Characterization of human palladin, a microfilament-associated protein. Mol Biol Cell. (2001) 12:3060–73. doi: 10.1091/mbc.12.10.3060

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Mori M, Nakagami H, Koibuchi N, Miura K, Takami Y, Koriyama H, et al. Zyxin mediates actin fiber reorganization in epithelial-mesenchymal transition and contributes to endocardial morphogenesis. Mol Biol Cell. (2009) 20:3115–24. doi: 10.1091/mbc.e09-01-0046

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Bare LA, Morrison AC, Rowland CM, Shiffman D, Luke MM, Iakoubova OA, et al. Five common gene variants identify elevated genetic risk for coronary heart disease. Genet Med. (2007) 9:682–9. doi: 10.1097/GIM.0b013e318156fb62

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Sahoo SK, Kim DH. Characterization of calumenin in mouse heart. BMB Rep. (2010) 43:158–63. doi: 10.5483/BMBRep.2010.43.3.158

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Sahoo SK, Kim DH. Calumenin interacts with SERCA2 in rat cardiac sarcoplasmic reticulum. Mol Cells. (2008) 26:265–9. doi: 10.1186/1476-4598-7-74

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: heart failure, ischemic cardiomyopathy, inflammatory activation, biomarker, diagnosis

Citation: Liu X, Xu S, Li Y, Chen Q, Zhang Y and Peng L (2021) Identification of CALU and PALLD as Potential Biomarkers Associated With Immune Infiltration in Heart Failure. Front. Cardiovasc. Med. 8:774755. doi: 10.3389/fcvm.2021.774755

Received: 13 September 2021; Accepted: 08 November 2021;
Published: 01 December 2021.

Edited by:

Ji Bihl, Marshall University, United States

Reviewed by:

Jyotsna Joshi, Mayo Clinic Arizona, United States
Owais Bhat, Virginia Commonwealth University, United States
Zijian Xiao, University of South China, China

Copyright © 2021 Liu, Xu, Li, Chen, Zhang and Peng. 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: Yuanyuan Zhang, zhangyy67@mail.sysu.edu.cn; Long Peng, pengl5@mail.sysu.edu.cn

These authors have contributed equally to this work

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.