Molecular effects of cardiac contractility modulation in patients with heart failure of ischemic aetiology uncovered by transcriptome analysis
E. Lyasnikova1 
K. Sukhareva1,2 M. Vander1 K. Zaitsev3 M. Firulyova3 A. Sergushichev3 M. Sorokina1 M. Trukshina1 V. Galenko1 T. Lelyavina1 L. Mitrofanova1 K. Simonova1 M. Abramov1 
G. Faggian2 
G. B. Luciani2 D. S. Lebedev1 
E. N. Mikhaylov1 M. Sitnikova1 
A. Kostareva1,4*
- 1Institute of Molecular Biology and Genetics, Almazov National Medical Research Centre, Saint Petersburg, Russia
- 2Graduate School of Life and Health Science, University of Verona, Verona, Italy
- 3Computer Technologies Laboratory, ITMO University, Saint Petersburg
- 4Department of Women’s and Children’s Health and Center for Molecular Medicine, Karolinska Institute, Stockholm, Sweden
Cardiac contractility modulation (CCM) is based on electrical stimulation of the heart without alteration of action potential and mechanical activation, the data on its fundamental molecular mechanisms are limited. Here we demonstrate clinical and physiological effect of 12 months CCM in 29 patients along with transcriptomic molecular data. Based on the CCM effect the patients were divided into two groups: responders (n = 13) and non-responders (n = 16). RNA-seq data were collected for 6 patients before and after CCM including 3 responders and 3 non-responders. The overall effect of CCM on gene expression was mainly provided by samples from the responder group and included the upregulation of the genes involved in the maintenance of proteostasis and mitochondrial structure and function. Using pathway enrichment analysis, we found that baseline myocardial tissue samples from responder group were characterized by upregulation of mitochondrial matrix-related genes, Z disc-protein encoding genes and muscle contraction-related genes. In summary, twelve months of ССM led to changes in signaling pathways associated with cellular respiration, apoptosis, and autophagy. The pattern of myocardial remodeling after CCM is associated with initial expression level of myocardial contractile proteins, adaptation reserves associated with mitochondria and low expression level of inflammatory molecules.
Introduction
Cardiac contractility modulation (CCM) has been actively studied in order to improve the prognosis for patients with heart failure and reduced ejection fraction (HFrEF) (1–3). CCM is relatively new method of electrophysiological therapy for patients with HFrEF which is based on electrical stimulation of the heart in an absolute refractory period (1). CCM does not lead to the emergence of an action potential and does not change the course of electrical and mechanical activation of the heart, but only provides a positive inotropic effect without increasing the myocardial oxygen demand (4–6). The molecular mechanisms underlying the positive inotropic effect of CCM are suggested to be mediated by the increase in intracellular calcium levels and changes in phosphorylation of the key proteins modulating the activity of sarcoplasmic calcium (3–7). Until now, the only study focused on transcriptomic and proteomic changes measured directly in human cardiac tissue under CCM has been published (8). This study included only a target analysis of genes involved in calcium metabolism after tree and six months of 3-month long CCM. The data on the global molecular events in myocardial tissue under long-term exposure to CCM are still missing. The present study aimed to identify differentially expressed genes in myocardium biopsies obtained from patients with HFrEF after twelve months of CCM using whole-transcriptome sequencing approach in order to associate these alterations with a type of CCM response and myocardial remodeling.
Materials and methods
The study protocol was approved by the local ethics committee of the Almazov National Medical Research Centre and complied with the requirements of the Declaration of Helsinki. All patients signed written informed consent. Patient characteristics and therapy are presented in Table 1.
Table 1. Demographic characteristics and clinical course of 6 patients with repeated endomyocardial biopsy samples, mean ± SD or median [Q25; Q75].
Total RNAs were isolated from myocardial biopsy specimens obtained from right ventricular side of interventricular septum and libraries for RNA sequencing were prepared using TruSeq Stranded mRNA kit (Illumina, USA). We used DeSeq2 to perform differential expression analysis and compared patients before CCM implantation against patients after implantation. Raw sequencing data are available at SRA, NCBI, under GSE251971. Details of RNA sequencing and bioinformatic analysis are described in Suppl. File.
Results
The mean follow-up period for all 29 patients with HFrEF of ischemic aetiology after CCM intervention was 11.8 ± 1.5 months, by the end of the first year, 11 patients (38%) had an ICD implanted mainly as a primary prevention of sudden cardiac death. After 6 and 12 months of CCM, a reduction in the NYHA (New York Heart Association) class accompanied by a significant improvement of echocardiographic, clinical, biochemical and quality of life parameters were reported (Supplementary Table S1). Based on the CCM effect on LVESV (Left Ventricular End-Systolic Volume), the study cohort was divided into two groups: responders (n = 13), where LVESV decreased by more than 10% compared to baseline after 6 months of CCM, and non-responders (n = 16), where parameters remained unaltered. These groups did not differ significantly in age, baseline hemodynamic and clinical parameters (Table 2). The degree of left ventricular (LV) reverse remodeling was significantly greater in the responder group compared to the non-responder group (Table 2). Both groups showed a significant reduction in the NYHA class, an increase in the distance 6MWT (6 Minute Walk Test) and quality of life according to EQ-5D data.
Table 2. Comparison of baseline characteristics and clinical course among the responder and non-responder groups (n = 29). Data presented as: n (%), mean ± SD or Me [Q25; Q75].
Repeated RNA-seq analysis was performed to identify common transcriptional changes related to CCM as well as to identify transcriptional changes unique to the responder or non-responder groups of patients, the samples were first collected during CCM implantation, and then in patients with later implantable cardioverter-defibrillator (ICD) insertion during the operation. Two consecutive RNA-seq data were collected for 6 patients, including 3 responders and 3 non-responders (Table 1, Supplementary Table S2). The mean time of ICD implantation with repeated myocardial biopsy sampling was 10.6 ± 2.6 months after CCM implantation and did not differ between the responder and non-responder groups. First, we compared transcriptional changes before and after CCM in all the patients. A total of 242 genes were differentially expressed (Figure 1A), with SYNPO, NRAP, and RORC being the most upregulated after CCM therapy and BTNL9, PTN, OGN, NOTCH4, MCF2l, and ADCY4 as the most downregulated. Many pathways related to Z-disc structure, mitochondria, and macroautophagy were upregulated in patients after CMM (Figure 1B). A further separate analysis of post CCM samples in responder and non-responder groups revealed that most of the genes upregulated in the entire post-CCM group were represented by the samples from the responder group. We revealed 32 differentially expressed genes, with HSPB6, MYO18B, CDKN1A, GALNT17, KCNK6, KCNJ4, and SNAP47 being the most upregulated genes; and BTNL9 and MCF2l being the most downregulated genes in CCM responder group (Figure 1C). These upregulated genes included signal recognition particle (SRP)−dependent co-translational protein targeting to membrane, many metabolic genes as well as genes encoding for actin-binding proteins (Figure 1D). Very few genes were upregulated in CCM non-responder group, while several genes (LUM, PTN, OGN) were significantly downregulated, including genes involved in cell adhesion and integrin-complex, collagen and heparin-binding processes (Figures 1E,F). The overall effect of CCM on global gene expression was mainly provided by samples from the responder group and included the upregulation of the genes involved in the maintenance of proteostasis and mitochondrial structure and function. Since macroautophagy was one of the upregulated pathways, we performed a detailed analysis of the genes involved in a muscle-specific type of autophagy—chaperone-assisted selective autophagy (CASA) and mitophagy. The gene set enrichment analysis revealed the significant upregulation (p < 5e-13) of the genes involved in CASA and mitophagy in patients after CCM (Figure 2).
Figure 1. Gene expression analysis by RNA sequencing of biopsy samples from patient before and after CCM treatment. (A) Volcano plot illustration RNA-Seq differential expression data. Pairwise comparisons is shown between all samples before and after therapy. (B) Gene enrichment analyses for comparison of all before and after CCM therapy samples. (C) Volcano plot illustration RNA-Seq differential expression data for pairwise comparisons samples between responders to therapy samples before and after CCM therapy. (D) Gene enrichment analyses for comparison of responders samples before and after CCM therapy. (E) Volcano plot illustration RNA-Seq differential expression data for pairwise comparisons samples between non-responders to therapy samples before and after CCM therapy. (F) Gene enrichment analyses for comparison of non-responders samples before and after CCM therapy.
Figure 2. Heat map illustrating genes associated with CASA process. P Pairwise comparisons is shown between all samples before and after therapy. Blue, negative log fold-change (log FC) indicates lower expression; red, positive log FC.
To further search the predictive expression markers that could distinguish responders from non-responders based on their baseline characteristics, we compared the baseline transcriptional profile in these groups. The statistical power of this test was lower than in the previous comparison since the number of analyzed samples was less (n = 6), and many factors were included into differential expression design (e.g., sex of the donor, batch), which resulted in the almost complete absence of differentially expressed genes (Supplementary Figure S1). We found that many mitochondrial pathways, mitochondrial matrix-related genes, Z disc-protein encoding genes and muscle contraction-related genes were upregulated in responders. In contrast, immune pathways including leukocyte cell-cell adhesion, and neutrophil degranulation were upregulated in non-responders. To exclude the difference in tissue composition, muscle cell cellularity, and the area of fibrosis, the morphological and voltage electrophysiology analysis were performed. We demonstrated that there were no differences in the total cardiomyocyte area and ECG voltage in myocardial samples from responder and non-responder groups (Table 3, Figure 3). Therefore, the observed differences in myocardial expression profile was not linked to different cellular composition.
Table 3. Histological analyses and intraoperative parameters of electrodes of 6 patients with repeated endomyocardial biopsy samples (pre- and post-fixation electrode values), mean ± SD.
Figure 3. Endomyocardial biopsies. (A–C) responders’ samples. (D) non-responders’ sample. (A) Endocardial fibrosis with a pronounced vacuolar degeneration of cardiomyocytes and fibrosis-braiding of cardiomyocytes. 65 years old, х100. (B) Small-focal fibrosis (or Focal fibrosis) with fibrosis-braiding of cardiomyocytes, endocardial fibrosis. 65 years old, ×50. (C) Fatty infiltration with endocardial fibrosis. 56 years old, ×50. (D) Small-focal fibrosis (Focal fibrosis) with fatty infiltration. 62 years old, х100. Masson's trichrome stain.
Discussion
In our study, we demonstrated that 6–12 months of CCM of ischemic HF patients (II-III NYHA) was associated with improvement of clinical parameters and quality of life. By 6 months of CCM, reduction of LVESV ≥ 10% and LVESV ≥ 15% were registered in 45% and 34% patients, respectively. These data are in line with data on 3 months of CCM in patients with HFrEF (III NYHA) (50% with ischemic aetiology) published by Zhang and co-authors who documented the decrease in LVESV ≥ 15% in 39% of cases (9). Of note, in our study, the long-term CCM did not result in significant dynamics of reverse myocardial remodeling in 55% of patients despite improvement in the clinical symptoms and exercise tolerance.
The reverse remodeling of LV under CCM is determined by molecular pathways and previous studies underlined the cardiomyocyte calcium-handling genes as the primary effectors of CCM impact (8, 10–12). However, we did not observe any effect of CCM on either of calcium-related proteins. One of the possible explanations could be the different time point of biopsy sampling in our and previously published studies (12 and 3 months correspondingly). Additional explanation could be the different technical approach utilizing RNA sequencing in a current study compared to previously reported RT-PCR data. Thus, our study provides the information on the later time points of CCM effect and illuminate new and previously unreported genes and pathways linked to a positive effect of CCM. Most of the gene expression effect on CCM raised from the responder cardiac samples leaving the non-responder group as almost “unreactive”. Among the most upregulated are mitochondrial matrix and mitophagy genes, mitochondrial metabolism-related genes, contractile sarcomeric genes as well as structural cardiomyocyte genes. Allover, the CCM-related improvement of myocardial function is associated with expression increase of cardiomyocyte structural and metabolic genes and downregulation of extracellular matrix and collagen synthesis-related genes. This is well in line with data reported by D'Onofrio et al., who demonstrated the improvement of inflammatory circulating biomarkers and markers of fibrosis such as collagen 3, collagen 4, C-cystatin and IL-6 after 6 and 12 months of CCM therapy in patient with LMNA-associated dilated cardiomyopathy (13). However, the observed modulation of gene expression still leaves the question whether the described molecular effects are specifically linked to CCM-related cardiomyocyte changes or simply reflect the cellular processes under positive cardiac remodeling and HFrEF positive dynamics. For example, upregulation of mitochondrial pathways, structural and Z disc-protein encoding genes and muscle contraction-related genes can represent the molecular signature of increased cardiomyocyte contractility and functional myocyte properties as a result of HF treatment independently on CCM. Further experimental and clinical studies with extended control groups and increased number of patients included will allow to answer these questions.
The involvement of several genes linked to CASA and mitophagy prompted deeper analysis of this gene set under CCM. We detected the relative increase in expression of such genes including BAG3, HSPB8, HSPA8, VSP and SQSTM1 as well as MFN2, VDAC, PINK1 and PRKN (Figure 2G). CASA has already been reported as essential process degrading damaged components of Z-disc and our data further confirms that positive CASA flux can be one of the attributes of restoration of cardiomyocyte contractile function (14–17).
The increase in CASA-related genes, structural, mitochondrial and contractile genes reflect the late CCM-mediated effects which do not include the direct involvement of Ca-operating genes described as the early CCM effects. Importantly, these effects were mostly represented by responder group samples. This underlines the importance of initial patient stratification and identification of those subjects who have the highest probability to benefit from CCM. We speculate the preservation of mitochondrial structure and metabolism as well as contractile and cytoskeletal apparatus determinate greater effect of CCM. In contrast, the increased expression of immune and inflammatory response genes is associated with moderate or no long-term functional response to CCM. We conclude that the baseline cardiomyocyte status—either more pro-contractile or more proinflammatory—is critical for the long term on cite and remote effects of CCM.
Study limitations
Due to the very small number of samples available for RNA sequencing the current study represents the pilot project aiming to underline the possible molecular effects of long-term CCM treatment. Further verification and validation of the described tendencies need to increase the number of samples and, possible, to extend the number of time points analyzed. In addition, the data obtained in frame of multicentre study rather than in single center study will allow more accurate and objective clinical data assessment.
Data availability statement
The data presented in the study are deposited in the SRA, NCBI repository, accession number GSE251971. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE251971
Ethics statement
The studies involving humans were approved by Ethics committee of the Almazov National Medical Research Centre and complied with the requirements of the Declaration of Helsinki. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
EL: Conceptualization, Data curation, Formal Analysis, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing. KSu: Data curation, Formal Analysis, Writing – original draft, Writing – review & editing. MV: Data curation, Methodology, Supervision, Writing – review & editing. KZ: Data curation, Formal Analysis, Investigation, Software, Writing – original draft, Writing – review & editing. MF: Data curation, Formal Analysis, Methodology, Software, Writing – review & editing. AS: Formal Analysis, Project administration, Software, Supervision, Writing – review & editing. MSo: Data curation, Formal Analysis and Writing – review & editing. MT: Investigation, Methodology, Writing – original draft. VG: Investigation, Methodology, Supervision, Validation, Writing – original draft. TL: Data curation, Investigation, Methodology, Supervision, Writing – original draft. LM: Conceptualization, Project administration, Supervision, Writing – review & editing. KSi: Project administration, Supervision, Writing – review & editing. MA: Project administration, Supervision, Writing – review & editing. GF: Project administration, Supervision, Writing – review & editing. GL: Project administration, Supervision, Writing – review & editing. DL: Methodology, Project administration, Supervision, Writing – review & editing. EM: Methodology, Project administration, Supervision, Writing – review & editing. MSi: Project administration, Supervision, Writing – review & editing. AK: Conceptualization, Data curation, Project administration, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article.
The work was supported by Russian Science Foundation grant number 20–15-00271П.
Acknowledgments
This work was supported by Anastasia Knyazeva and Artem Kiselev in the preparation of experiments and the compilation of RNA libraries.
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
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.1321005/full#supplementary-material
Supplementary Figure S1
Gene expression analysis by RNA sequencing of biopsy samples from patient before and after CCM treatment. (A) Volcano plot illustration RNA-seq differential expression data. Pairwise comparisons is shown between all responders and non-responders samples before therapy. (B) Gene enrichment analyses for comparison of all responders and non-responders samples before therapy.
References
1. Rao IV, Burkhoff D. Cardiac contractility modulation for the treatment of moderate to severe HF. Expert Rev Med Devices. (2021) 18(1):15–21. doi: 10.1080/17434440.2020.1853525
2. Fastner C, Yuecel G, Rudic B, Schmiel G, Toepel M, Burkhoff D, et al. Cardiac contractility modulation in patients with ischemic versus non-ischemic cardiomyopathy: results from the MAINTAINED observational study. Int J Cardiol. (2021) 342:49–55. doi: 10.1016/j.ijcard.2021.07.048
3. Lyon AR, Samara MA, Feldman DS. Cardiac contractility modulation therapy in advanced systolic heart failure. Nat Rev Cardiol. (2013) 10(10):584–98. doi: 10.1038/nrcardio.2013.114
4. Butter C, Wellnhofer E, Schlegl M, Winbeck G, Fleck E, Sabbah HN. Enhanced inotropic state of the failing left ventricle by cardiac contractility modulation electrical signals is not associated with increased myocardial oxygen consumption. J Card Fail. (2007) 13(2):137–42. doi: 10.1016/j.cardfail.2006.11.004
5. Goliasch G, Khorsand A, Schutz M, Karanikas G, Khazen C, Sochor H, et al. The effect of device-based cardiac contractility modulation therapy on myocardial efficiency and oxidative metabolism in patients with heart failure. Eur J Nucl Med Mol Imaging. (2012) 39(3):408–15. doi: 10.1007/s00259-011-1977-8
6. Tschope C, Kherad B, Klein O, Lipp A, Blaschke F, Gutterman D, et al. Cardiac contractility modulation: mechanisms of action in heart failure with reduced ejection fraction and beyond. Eur J Heart Fail. (2019) 21(1):14–22. doi: 10.1002/ejhf.1349
7. Yu CM, Chan JY, Zhang Q, Yip GW, Lam YY, Chan A, et al. Impact of cardiac contractility modulation on left ventricular global and regional function and remodeling. JACC Cardiovasc Imaging. (2009) 2(12):1341–9. doi: 10.1016/j.jcmg.2009.07.011
8. Butter C, Rastogi S, Minden HH, Meyhofer J, Burkhoff D, Sabbah HN. Cardiac contractility modulation electrical signals improve myocardial gene expression in patients with heart failure. J Am Coll Cardiol. (2008) 51(18):1784–9. doi: 10.1016/j.jacc.2008.01.036
9. Zhang Q, Chan YS, Liang YJ, Fang F, Lam YY, Chan CP, et al. Comparison of left ventricular reverse remodeling induced by cardiac contractility modulation and cardiac resynchronization therapy in heart failure patients with different QRS durations. Int J Cardiol. (2013) 167(3):889–93. doi: 10.1016/j.ijcard.2012.01.066
10. Ning B, Qi X, Li Y, Liu H, Zhang F, Qin C. Biventricular pacing cardiac contractility modulation improves cardiac contractile function via upregulating SERCA2 and miR-133 in a rabbit model of congestive heart failure. Cell Physiol Biochem. (2014) 33(5):1389–99. doi: 10.1159/000358705
11. Imai M, Rastogi S, Gupta RC, Mishra S, Sharov VG, Stanley WC, et al. Therapy with cardiac contractility modulation electrical signals improves left ventricular function and remodeling in dogs with chronic heart failure. J Am Coll Cardiol. (2007) 49(21):2120–8. doi: 10.1016/j.jacc.2006.10.082
12. Tschope C, Van Linthout S, Spillmann F, Klein O, Biewener S, Remppis A, et al. Cardiac contractility modulation signals improve exercise intolerance and maladaptive regulation of cardiac key proteins for systolic and diastolic function in HFpEF. Int J Cardiol. (2016) 203:1061–6. doi: 10.1016/j.ijcard.2015.10.208
13. D’Onofrio A, Palmiero G, D’Alterio G, De Vivo S, Maione B, Leonardi S. First human implant of the cardiac contractility modulation in patient with dilated cardiomyopathy-related laminopathy. HeartRhythm Case Rep. (2023) 9(6):381–5. doi: 10.1016/j.hrcr.2023.03.011
14. Kanamori H, Yoshida A, Naruse G, Endo S, Minatoguchi S, Watanabe T, et al. Impact of autophagy on prognosis of patients with dilated cardiomyopathy. J Am Coll Cardiol. (2022) 79(8):789–801. doi: 10.1016/j.jacc.2021.11.059
15. Li Y, Meng W, Hou Y, Li D, Wang X, Wu K, et al. Dual role of mitophagy in cardiovascular diseases. J Cardiovasc Pharmacol. (2021) 78(1):e30–e9. doi: 10.1097/FJC.0000000000001046
16. Martin TG, Myers VD, Dubey P, Dubey S, Perez E, Moravec CS, et al. Cardiomyocyte contractile impairment in heart failure results from reduced BAG3-mediated sarcomeric protein turnover. Nat Commun. (2021) 12(1):2942. doi: 10.1038/s41467-021-23272-z
Keywords: CCM, autophagy, optimizer, heart failure, sarcomeric genes, RNA sequencing
Citation: Lyasnikova E, Sukhareva K, Vander M, Zaitsev K, Firulyova M, Sergushichev A, Sorokina M, Trukshina M, Galenko V, Lelyavina T, Mitrofanova L, Simonova K, Abramov M, Faggian G, Luciani GB, Lebedev DS, Mikhaylov EN, Sitnikova M and Kostareva A (2024) Molecular effects of cardiac contractility modulation in patients with heart failure of ischemic aetiology uncovered by transcriptome analysis. Front. Cardiovasc. Med. 11:1321005. doi: 10.3389/fcvm.2024.1321005
Received: 27 October 2023; Accepted: 8 January 2024;
Published: 1 February 2024.
Edited by:
Syamal K. Bhattacharya, University of Tennessee Health Science Center (UTHSC), United StatesReviewed by:
Fuyang Zhang, Air Force Medical University, ChinaAntonio D’Onofrio, AORN dei Colli- Monaldi Hospital, Italy
© 2024 Lyasnikova, Sukhareva, Vander, Zaitsev, Firulyova, Sergushichev, Sorokina, Trukshina, Galenko, Lelyavina, Mitrofanova, Simonova, Abramov, Faggian, Luciani, Lebedev, Mikhaylov, Sitnikova and Kostareva. 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: A. Kostareva anna.kostareva@ki.se