Proteomic and phosphoproteomic profiling in heart failure with preserved ejection fraction (HFpEF)

Abstract

Although the prevalence of heart failure with preserved ejection fraction (HFpEF) is increasing, evidence-based therapies for HFpEF remain limited, likely due to an incomplete understanding of this disease. This study sought to identify the cardiac-specific features of protein and phosphoprotein changes in a murine model of HFpEF using mass spectrometry. HFpEF mice demonstrated moderate hypertension, left ventricle (LV) hypertrophy, lung congestion and diastolic dysfunction. Proteomics analysis of the LV tissue showed that 897 proteins were differentially expressed between HFpEF and Sham mice. We observed abundant changes in sarcomeric proteins, mitochondrial-related proteins, and NAD-dependent protein deacetylase sirtuin-3 (SIRT3). Upregulated pathways by GSEA analysis were related to immune modulation and muscle contraction, while downregulated pathways were predominantly related to mitochondrial metabolism. Western blot analysis validated SIRT3 downregulated cardiac expression in HFpEF vs. Sham (0.8 ± 0.0 vs. 1.0 ± 0.0; P < 0.001). Phosphoproteomics analysis showed that 72 phosphosites were differentially regulated between HFpEF and Sham LV. Aberrant phosphorylation patterns mostly occurred in sarcomere proteins and nuclear-localized proteins associated with contractile dysfunction and cardiac hypertrophy. Seven aberrant phosphosites were observed at the z-disk binding region of titin. Additional agarose gel analysis showed that while total titin cardiac expression remained unaltered, its stiffer N2B isoform was significantly increased in HFpEF vs. Sham (0.144 ± 0.01 vs. 0.127 ± 0.01; P < 0.05). In summary, this study demonstrates marked changes in proteins related to mitochondrial metabolism and the cardiac contractile apparatus in HFpEF. We propose that SIRT3 may play a role in perpetuating these changes and may be a target for drug development in HFpEF.

Introduction

Heart failure (HF) is a clinical syndrome caused by abnormalities in the heart that limit its ability to fill or eject blood (1). Heart failure with preserved ejection fraction (HFpEF) is symptomatic clinical HF where left ventricular (LV) ejection fraction (EF) is preserved (LVEF ≥ 50%), and presently accounts for about 50% all HF clinical presentations. However, unlike HF with reduced EF (HFrEF), where LVEF is < 50%, there are limited evidence-based therapies for HFpEF (2–4). In addition to its escalating prevalence, HFpEF morbidity (5) and mortality (6) continues to increase. Central to HFpEF is the involvement of both cardiac and extra-cardiac abnormalities (7, 8). In contrast to HFrEF, HFpEF is highly associated with comorbidities and as such is a heterogenous multisystem disorder involving the heart, pulmonary, renal, adipose tissue, skeletal muscle, immune/inflammatory signaling and the vascular system (9, 10). Patients with HFpEF are generally older, more often female and have a predominance of comorbidities, such as hypertension, obesity, type 2 diabetes, atrial fibrillation, renal dysfunction, etc. (11, 12). However, the specific etiologies by which patients develop HFpEF are variable. Thus, a precision-based approach is needed to identify pathogenic mechanisms in HFpEF (10, 13).

Proteomic studies are powerful tools that allow for large-scale characterization of the entire protein phenotype in a biological system (14). Alterations in proteome patterns, such as global changes in protein expression and post-translational modifications (PTMs), are often indicative of marked changes in functional stages in health and disease (15). Thus, investigating the varying patterns of the proteome may provide insights into pathogenic pathways (16) and these protein signatures may facilitate rapid screening of the efficacy of novel treatments and aid in drug development (17, 18).

Previous proteomic studies have identified protein changes in dilated cardiomyopathy, atherosclerosis, and atrial fibrillation (19–23) and these types of studies likely provided a deeper mechanistic understanding of the molecular pathways in HF. For example, cardiac tissue from patients with HFrEF demonstrated protein modifications associated with cardiac metabolism, cardiac remodeling, and impaired cardiac contractility (24–27). Additionally, differentially regulated pathways by proteomic signatures were observed in HFrEF vs. HFpEF patients, which is consistent with the predominant view that the underlying pathophysiology in these two diseases are largely different, and thus the variable response to therapies. This difference is exemplified by Adamo et al., where blood samples from both HFrEF and HFpEF patients demonstrated increased growth factor signaling and increased angiogenesis markers, while proteomic signatures from only HFpEF patients showed increased humoral immunity and those from HFrEF patients showed increased extracellular matrix remodeling markers, consistent with active cardiac remodeling (28). These findings underscore the potential that high-performance proteomics, in combination with clinical assessment, may identify unique targets in specific groups of HF patients.

Although HFpEF is greatly impacted by the obesity and diabetes pandemic, hypertension remains the most prevalent and modifiable risk factor in HFpEF and is implicated in both its pathogenesis and prognosis (12, 29). Hypertensive HFpEF pathophysiology extends beyond the emphasis on LV hypertrophy development and diastolic dysfunction to impaired myocardial contractility, left atrial myopathy, cardiomyocyte remodeling, macro- and microvascular dysfunction, to systemic inflammation, fibrosis, and collagen deposition. However, despite this knowledge a paucity of therapies exists for HFpEF. Here, we applied a deep quantitative proteomics and phosphoproteomic profiling approach to identify molecular protein signatures that are altered in HFpEF in a well characterized murine model of hypertension-associated HFpEF, the SAUNA model (SAlty drinking water/Unilateral Nephrectomy/Aldosterone), which recapitulates the human HFpEF phenotype (30–37) (Supplementary Figure 1). Using an unbiased and comprehensive analysis, we report a systematic, large-scale study of pathway, metabolic and organelle level changes that occur in the left ventricle of this HFpEF murine model.

Material and methods

All procedures related to the handling and surgery of the mice conformed to the Guide for the Care and Use of Laboratory Animals published by the United States National Institutes of Health and were approved by the Institutional Animal Care and Use Committee at Boston University School of Medicine.

SAUNA model of HFpEF

As previously described (30–33, 35–37), eight-week-old male C57BL/6J mice (Jackson Laboratories) were anesthetized with 80–100 mg/Kg ketamine and 5–10 mg/Kg xylazine intraperitoneally. Mice (20–25 g) then underwent uninephrectomy, received either a continuous infusion of saline (Sham) or d-aldosterone (0.30 μg/h, Sigma-Aldrich, St. Louis, MO, United States; HFpEF) for 4 weeks via osmotic minipumps (Alzet, Durect Corp., Cupertino, CA, United States) and were maintained on 1% sodium chloride drinking water.

Physiological measurements

Blood pressure and echocardiographic measurements were performed at the end of the 4 weeks. Systolic blood pressure was measured using a non-invasive tail-cuff blood pressure analyzer (BP-2000 Blood Pressure Analysis System; Visitech Systems Inc., Apex, NC, United States). Transthoracic echocardiography was performed using a Vevo770 High-Resolution in vivo Micro-Imaging System and a Real-Time Micro Visualization 707B Scanhead (VisualSonic Inc., Toronto, ON, Canada) as previously described (33). Briefly, interventricular septum wall thickness (IVST), left ventricle (LV) posterior wall thickness (LVPWT), LV end-diastolic diameter (LVEDD), LV end-systolic diameter (LVESD), and LV ejection fraction (LVEF) were measured. As a measure of systolic function and cardiac contractility fractional shortening (FS) was calculated as follows (LVEDD-LVESD/LVEDD) × 100. Total wall thickness (TWT) was derived from an average of the IVST and LVPWT. Relative wall thickness (RWT) was calculated as 2× LVPWT/LVEDD. LV mass was calculated using the formula described by Kiatchoosakun et al. (38). As diastolic function is sensitive to heart rate (HR) and loading conditions, HR was maintained at ∼350 bpm during these measurements (39). Pulse wave measurements were then recorded and analyzed blinded to group.

Histopathological analyses

Paraffin-embedded sections (5 μm) of the mid-LV were stained with hematoxylin and eosin (H&E, Sigma-Aldrich) to measure LV cardiac myocyte cross-sectional area. Microscopy images (BZ-9000 BioRevo microscope, Keyence Corp. of America, Itasca, IL, United States) were analyzed blinded to group identity using ImageJ measuring software (National Institutes of Health, Bethesda, MD, United States).

Tissue sample preparation for proteomics and phosphoproteomics

Left ventricle samples from 4 mice/group were processed as previously described (22, 40–42). Briefly, freshly thawed samples were homogenized on ice in with a mixer mill MM 400 (Retsch USA Verder Scientific Inc., Newtown, PA, United States) in 10 volumes of 8 M urea, 50 mM ammonium bicarbonate, 2 mM dithiothreitol, and protease and phosphatase inhibitor cocktails (Roche Applied Science, Indianapolis, IN, United States). Tissue homogenate was then sonicated with a probe sonicator (Branson Ultrasonics Corporation, North Billerica, MA, United States) and centrifuged. After centrifugation, supernatant was decanted and total protein in each sample was determined using a modified “microtiter plate” version of the Bradford assay (Sigma-Aldrich). For phosphoproteomics experiments, aliquots containing 300 μg of protein were alkylated with 5 mM iodoacetamide for additional 45 min at room temperature in the dark. Samples were then diluted eight-fold with 50 mM ammonium bicarbonate and digested overnight with sequencing-grade trypsin (#90057, Thermo Fisher Scientific Inc., Waltham, MA, United States). Digestion was stopped by acidification to a final concentration of 1% (v/v) formic acid and the peptide solutions were desalted using disposable C18 Sep-Pak syringes (Waters Corporation, Milford, MA, United States) and lyophilized to dryness following manufacturer’s instructions.

Tandem mass tag (TMT) labeling

Peptide concentrations were determined by a colorimetric peptide assay kit (Thermo Fisher Scientific Inc., Waltham, MA, United States) and an aliquot of 100 μg was placed in 100 μl of 100 mM triethylammonium bicarbonate. Peptides were labeled with 0.4 mg of TMT label (TMT10plex™ Isobaric Label Reagent Set, Thermo Fisher Scientific Inc., Waltham, MA, United States). All samples were labeled in the same TMT-batch, representing reporter tags 126C, 127N, 127C, 128C, 129N, 129C, 130N, and 131N. Labeled samples were pooled, and 95% was set aside for phosphopeptide enrichment. The remaining 5% of labeled peptides and the phosphopeptide enriched samples were analyzed separately by mass spectrometry.

Phosphopeptide enrichment

Phosphopeptides were selectively enriched by binding to titanium dioxide (TiO₂) beads (Titansphere Phos-TiO Bulk 10 μm, GL Sciences, Tokyo, Japan) (43). Briefly, peptides were resuspended in 200 μl 80% acetontirile, 6% trifluoroacetic acid and incubated for 10 min with 10 μl of slurry containing TiO₂ beads. Unbound peptides and supernatant were decanted, and the beads were washed three times with a wash buffer containing 50% acetonitrile and 1% trifluoroacetic acid. After final decanting, the beads were incubated for 10 min with elution solution containing 25% ammonium hydroxide and 50% acetonitrile and the eluate was carefully removed and dried prior to mass spectrometry analysis.

Mass spectrometry analysis

Tryptic peptide mixtures and enriched phosphopeptides were analyzed by nano-scale high-performance liquid chromatography (Proxeon EASY-Nano system, Thermo Fisher Scientific Inc., Waltham, MA, United States) and online nano electrospray ionization tandem mass spectrometry (Q-Exactive HF-X mass spectrometer; Thermo Fisher Scientific Inc., Waltham, MA, United States). Briefly, samples were loaded in aqueous 0.1% (v/v) formic acid via a trap column (75 μm i.d. × 2 cm, Acclaim PepMap100 C18 3 μm, 100 Å, Thermo Fisher Scientific) and peptides were resolved over an Easy-Spray analytical column (50 cm × 75 μm ID, PepMap RSLC C18, Thermo Fisher Scientific) by an increasing mobile phase B. Mobile phase A consisted of 2% acetonitrile and 0.1% formic acid, and organic phase B contained 80% acetonitrile and 0.1% formic acid. Reverse phase separation was performed over 120 min at a flow rate of 300 nl/min. Eluted peptides were ionized directly into the mass spectrometer using a nanospray ion source. The mass spectrometer was operated in positive ion mode with a capillary temperature of 300 C, and with a potential of 2,100 V applied to the frit. Tandem mass spectrometry (MS/MS) was performed using high-energy collision-induced disassociation and 10 MS/MS data-dependent scans (45,000 resolution) were acquired in profile mode alongside each profile mode full-scan mass spectra (120,000 resolution) as reported previously (44). The automatic gain control (AGC) for MS scans was 1 × 10⁶ ions with a maximum fill time of 60 ms. The AGC for MS/MS scans was 3 × 10⁴, with 80 ms maximum injection time, 0.1 ms activation time, and 33% normalized collision energy. To avoid repeated selection of peptides for MS/MS a dynamic exclusion list was enabled to exclude all fragmented ions for 60 s.

Protein identification

Data files (RAW format) were searched using the standard workflow of MaxQuant (version 1.3.0.5)¹ under standard settings using the entire Swiss-Prot mouse database² downloaded January 24, 2019, allowing for two missed trypsin cleavage sites, carbamidomethylation of cysteine (fixed) and variable oxidation of methionine, protein N-terminal acetylation and phosphorylation of STY residues. Precursor ion tolerances were 20 ppm for first search and 4.5 ppm for a second search. The MS/MS peaks were de-isotoped and searched using a 20-ppm mass tolerance. A stringent false discovery rate threshold of 1% was used to filter candidate peptide, protein, and phosphosite identifications. The datasets generated for this study have been deposited and publicly available at the PRIDE Archive, proteomics data repository (European Bioinformatics Institute, European Molecular Biology Laboratory) with the data set identifier PXD033501.

Bioinformatics analysis

The searched intensity data were filtered, normalized, and clustered using Omics Notebook (45). Filtering was performed to remove any proteins or phosphopeptides not quantified in at least 70 percent of samples, with 2,905 and 281 proteins and phosphopeptides passing the filter, respectively. After filtering, both datasets showed low levels of sparsity and no missing value imputation was performed. The LIMMA R package was used for LOESS normalization and differential expression analysis (46). A combined ranked list for both sets was generated where duplicate gene entries were removed to keep the entry with the highest absolute rank value.

GSEA analysis

Gene Set Enrichment analysis (GSEA) software from the fgsea R package was used to compute gene set enrichment after ranking proteins by differential expression in HFpEF vs. Sham (45, 47, 48). Briefly, GSEA was used in rank mode along with gene sets downloaded from the Bader Lab (Mouse_GOBP_AllPathways_no_GO_iea_October_01_2018_ symbol.gmt)³ (49, 50). GSEA results were visualized using the Enrichment Map app (Version 3.1) in Cytoscape (Version 3.6.1) and highly related pathways were grouped into a theme and labeled by AutoAnnotate (version 1.2). For the merged gene set analyses, we applied an enrichment P < 0.01 and FDR ≤ 0.1 cutoffs and calculated overlap between gene set annotations using a combination of Jaccard and overlap coefficients with a cutoff of 0.375.

Titin isoform analysis

Additional studies were performed to investigate changes in titin isoforms in HFpEF. Briefly, LV protein lysates from Sham (N = 7) and HFpEF (N = 11) mice were extracted and electrophoresed in 1% agarose gels using a SE600X vertical gel system (Hoefer Inc., Holliston, MA, United States) as previously described (51). Gels were run at 15 mA constant current, stained with Neuhoff’s Coomassie (52), and then scanned using Epson Perfection V750 PRO scanner (Epson America Inc., Los Alamitos, CA, United States) and analyzed using One-D scan EX analysis software (Scanalytics Inc., Rockville, MD, United States). The integrated optical density of titin and total myosin heavy chain (MHC) was determined as a function of the slope of the linear range between integrated optical density and loaded volume (53). The expression of compliant N2BA titin, stiffer N2B titin and total titin (TT) was normalized to the expression of total MHC. The expression of titin degraded product (T2) was normalized to the TT expression.

SIRT3 immunoblotting analysis

Protein lysates were extracted from LV tissue using in ice-cold RIPA buffer as previously described (54). Equal amounts of protein were then subjected to electrophoresis in SDS-polyacrylamide gel under reducing conditions and blotted to polyvinylidene difluoride (PVDF) membranes using the Bio-Rad Transblot Turbo Transfer System (Hercules, CA, United States). The membranes were blocked in 5% BSA, 0.1% Tween-20 in tris-buffered saline for 1 h at room temperature and then incubated overnight at 4°C with rabbit anti-SIRT3 antibody (Cell Signaling Technology, Inc., Danvers, MA, United States, #5490; 1:1.000). Membranes were then washed with tris-buffered saline and incubated with respective horseradish peroxidase (HRP)-conjugated secondary antibodies for 1 h in room temperature: anti-rabbit antibody (R&D system, HAF008; 1:5,000). Immune complexes were detected with the enhanced chemiluminescence ECL detection system (Bio-Rad, #1705060) in the ImageQuant LAS 4000 biomolecular imaging system (GE Healthcare, Pittsburgh, PA, United States). The intensity of bands for each protein was normalized to the loading control mouse anti-GAPDH (Abcam, Ab8245; 1:10.000).

Statistical analysis

Proteomics and phosphoproteomics differential analysis were based on a moderated t-test and performed using R: A language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria) (45, 55). For histology analysis, titin isoform studies and SIRT3 expression, data are shown as mean ± SEM and statistical significance of differences was assessed using the Student’s t-test (two sided). In those cases when data were not sampled as a normal distribution, non-parametric Mann–Whitney U test was used. P ≤ 0.05 values were considered significant. These statistical tests were performed using GraphPad Prism software (GraphPad Software Inc., La Jolla, CA, United States).

Results

Mouse model of HFpEF

As previously described (30–33, 35–37), salty drinking water, unilateral nephrectomy, and chronic exposure to aldosterone (SAUNA) induced hypertension associated HFpEF in mice. Compared to Sham, HFpEF mice demonstrated a moderate increase in systolic blood pressure (137.8 ± 7.0 mmHg vs. 115.4 ± 6.0 mmHg; P < 0.05), lung congestion (4.5 ± 0.1 vs. 4.0 ± 0.1 P < 0.01), and LV hypertrophy, measured by the LV weight-to-total body weight ratio (3.7 ± 0.1 mg/g vs. 3.3 ± 0.1 mg/g; P < 0.05). Additionally, cardiomyocyte size was increased 1.2-fold in HFpEF mice vs. Sham; P < 0.05 (Supplementary Figure 2).

Echocardiography demonstrated preserved LVEF and increased LV mass (107.5 ± 4.9 mg vs. 78.2 ± 7.9 mg in Sham; P < 0.05; Table 1). Wall thickness was significantly increased in HFpEF and there was evidence of concentric hypertrophy, as demonstrated by the increased relative wall thickness (0.7 ± 0.1 vs. 0.5 ± 0.0 in Sham; P < 0.005). As previously shown (33), LV end-systolic dimensions and end-diastolic dimensions were also decreased in HFpEF (Table 1). HFpEF mice had impaired diastolic function, characterized by an increase in isovolumetric relaxation time (24.3 ± 2.6 ms vs. 14.4 ± 1.6 ms in Sham; P < 0.05).

Comparison to human HFpEF: Recently, two clinical scores (HFA-PEFF and H2FPEF) were developed to standardize the clinical diagnosis of human HFpEF. However, a discrepancy exists between these scores (56). The H2FPEF score largely includes clinical parameters whereas the HFA-PEFF score includes predominantly echocardiographic measures and natriuretic peptides. The HFA-PEFF score can rule in human HFpEF with high specificity (93%) and positive predictive value (98%) when the score is high (5–6 points) (57). As such, the translational utility of the HFpEF SAUNA mouse model was demonstrated in the context of this HFpEF score with a HFA-PEFF score of ≥ 6 as described by Withaar et al. (58), where a score of ≥ 5 is a high probability of clinical HFpEF.

Proteome profile of the left ventricle in HFpEF

To achieve comprehensive evaluation of the cardiac signaling that is seen in HFpEF, a global quantitative proteome and phosphoproteome profile was performed in LV cardiac tissue obtained from HFpEF mice and their respective Shams (N = 4 mice/group; Figure 1).

FIGURE 1

Proteomics analysis found a total of 2,905 identified proteins that were then used for comparative analysis (Supplementary Table 1). Among them, 897 proteins were differentially expressed between HFpEF and Sham LV, with 19% of these being predominantly higher in HFpEF than in Sham (P < 0.05; Figures 2A,B).

FIGURE 2

Systematic evaluation of the datasets revealed abundant changes in sarcomeric proteins, namely skeletal alpha (α)-actin (ACTA1; P = 0.000039), beta (β)-myosin heavy chain (MYH7; P = 0.006963), myosin heavy chain 9 (MYH9; P = 0.000408), tropomyosin alpha (α)-1 chain (TPM1; P = 0.048698); the mitochondria-related proteins mitofusin 1 (MFN1; P = 0.001059), mitochondrial dynamin like GTPase (aka optic atrophy protein 1, OPA1; P = 0.046441) and transcription factor A mitochondrial (TFAM; P = 0.005837); and the NAD-dependent protein deacetylase sirtuin-3 (SIRT3; P = 0.000914), recently implicated in cardiac function and cardiac stress responsiveness in HFpEF (59, 60) (Figure 2B and Supplementary Table 1).

Impaired mitochondrial function and oxidative metabolism of energy substrates in HFpEF

There was an extensive reduction in the abundance of proteins involved in cardiac metabolism in the LV of HFpEF mice, including the oxidation of free fatty acid (FFA), pyruvate, and ketone bodies. Significant changes are summarized in Figure 3. These include:

FIGURE 3

(I) β-oxidation related enzymes, implicated in FFA metabolism to acetyl-CoA, such as acyl-CoA dehydrogenase (ACAD) family member 11 (ACAD11; P = 0.00002), long-chain specific ACAD (ACADL; P = 0.00449), short-chain specific ACAD (ACADS; P = 0.00896), short-branched chain specific ACAD (ACADSB; P = 0.011317), 3-ketoacyl-CoA thiolase (ACAA2, P = 0.00609), enoyl-CoA hydratase (ECHS1, P = 0.02647), hydroxyacyl-CoA dehydrogenase trifunctional multienzyme complex (HADH) beta (β)-subunit (HADHB, P = 0.00729) and HADH alpha (α)-subunit (HADHA, P = 0.03235).

(II) the pyruvate oxidation enzyme pyruvate dehydrogenase X component (PDHX, P = 0.00961), which is part of the pyruvate dehydrogenase complex that catalyzes pyruvate to acetyl-CoA; and

(III) the ketone metabolism enzyme succinyl-CoA:3-keto- acid coenzyme A transferase 1 (OXCT1, P = 0.01237), which catalyzes ketone bodies and produces acetyl-CoA for the tricarboxylic acid (TCA) cycle.

These cumulative results suggest that the energy substrates for mitochondrial oxidative metabolism may be inefficient in HFpEF. Interestingly, although there were no significant alterations in the protein signature of fatty acid and glucose transporters (CD36 and GLUT1 and 4, respectively) in HFpEF, there was an upregulation of the ketone bodies transporter monocarboxylate transporter 1 (SLC16A1, P = 0.00376) in the LV of HFpEF.

Additional analysis revealed that the mitochondrial proteins involved in the TCA cycle were also significantly decreased in the LV of HFpEF mice. These mitochondrial enzymes, namely citrate synthase (CS, P = 0.008068), succinyl-CoA ligase beta subunit (SUCLA2, P = 0.02523), isocitrate dehydrogenase (IDH1, P = 0.00105) and pyruvate carboxylase (PC, P = 0.00206), are required to catalyze acetyl-CoA and produce essential intermediates for the biosynthesis process, and most importantly, high energy molecules such as nicotinamide adenine dinucleotide (NADH) and flavin adenine dinucleotide (FADH₂) for the electron transport chain (ETC). Subsequent analysis then showed that 27 proteins involved in the ETC (namely the respiratory complex I, III, IV, and V) were also differentially expressed between HFpEF and Sham. Of these 27 proteins, 19 proteins were significantly reduced in the HFpEF, suggesting impaired ETC, which was consistent with an additional reduction of the uncoupling protein 3 (UCP3, P = 0.02765).

Lastly, additional proteins involved in mitochondrial biogenesis (transcription factor A, TFAM, P = 0.00584) and fusion (mitofusin-1, MFN1, P = 0.001056 and dynamin-like 120 kDa protein, OPA1, P = 0.04644) were similarly decreased in the LV tissue from HFpEF mice.

These findings (Figure 3) suggest that mitochondrial dysfunction may lead to inefficient metabolism of energy substrates, possibly contributing to an energy deficit and thus affecting cardiac function in HFpEF.

Pathway enrichment analysis

Pathway enrichment analyses of the proteomics and phosphoproteomics combined datasets were performed by means of GSEA, which detects biology-driven gene sets of canonical pathways from databases of molecular signatures (61). These analyses revealed that the most relevant and over-represented (enriched) biological annotations in the LV from HFpEF to be: (I) processes involving immune system modulation, (II) cardiac muscle cell development and differentiation, and (III) muscle contraction (Table 2). These processes included positive regulation of cytokine production (GO:0001819; P = 0.0000), striated muscle contraction (Wikipathway; P = 0.0000), positive regulation of adaptive immune response (GO:0002821; P = 0.00578), cardiac muscle cell development (GO:0055013; P = 0.03158) and cardiac muscle cell differentiation (GO: 0055007; P = 0.03571). In contrast, the downregulated pathways were related to a multitude of GO terms associated with cellular metabolism (Table 3). This is consistent with the earlier data from Figure 3, where pathways and processes involving acetyl-CoA metabolic process (GO:0006084, P = 0.00000), fatty acid metabolic process (GO:0006631, P = 0.00000), acyl-CoA biosynthesis process (GO:0071616, P = 0.0000), fatty acid oxidation (GO:0019395, P = 0.00111), coenzyme metabolic process (GO:0006732, P = 0.01015) were significantly reduced in HFpEF.

Phospho-proteome profile of the left ventricle in HFpEF

We next investigated the phosphoproteomics dataset. Phosphoproteomics analysis profiled 281 mouse reference protein sequences, of which 240 mapped to serine, 37 mapped to threonine and 3 mapped to tyrosine residues, consistent with the expected 90:9:1 cellular distribution ratio (22). The abundance of 72 phosphosites was differentially altered (elevated or reduced) between HFpEF and Sham (P < 0.05; Figures 4A,B).

FIGURE 4

Aberrant phosphorylation patterns occurred on proteins linked to disparate subcellular compartments, ranging from sarcomeric proteins (LIM domain-binding protein 3, LDB3; myozenin 2, MYOZ2; titin, TTN), to nuclear-localized proteins (BAG family molecular chaperone regulator 3, BAG3; high mobility group protein HMG-I/HMG-Y, HMGA1) with established links to cardiac contractile function, cardiac hypertrophy and/or cardiomyopathy (Figure 4B and Supplementary Table 2).

Left ventricular titin expression and phosphorylation in HFpEF

Despite global proteomics not showing a significant change in total titin in the LV between HFpEF and Sham mice, extensive phosphorylation changes across titin were observed in HFpEF vs. Sham. We identified 22 titin phosphosites (including 76% serines, 19% threonines, and 4% tyrosines) in the phospho-dataset, and among them, seven phosphosites were P < 0.05, all considered class 1 (localization probability 0.75–1.00) (Table 4). Interestingly, five of these were located at the z disk binding region (S262, S264, T266, S1411, and S1415) while the remainder were residues in the C-terminal region (S34464, T34467), suggesting changes in the mechano-sensing activity of titin. These regions are known to act as titin “hotspots,” which respond to mechanical stress and regulate specific actions such as activating the hypertrophic gene program or interacting with the protein quality control machinery (62, 63).

Additionally, high-resolution gel electrophoresis was performed to further examine other potential switches in titin isoform expression which may affect titin stiffness, i.e., to quantitatively detect changes in the stiffer N2B or the compliant N2BA isoforms of titin (Figure 5A). As expected, there were no changes in total titin (TT) expression between HFpEF and Sham mice. However, as previously described (64), N2B expression was significantly increased in the LV from HFpEF compared to Sham mice (0.144 ± 0.010 vs. 0.127 ± 0.010; P < 0.05). Neither N2BA expression, N2BA/N2B ratio nor titin degradation were differentially altered between HFpEF and Sham mice (Figure 5B).

FIGURE 5

Left ventricular expression of SIRT3 in HFpEF mice

Accumulating evidence suggests that SIRT3 plays a critical role in the development of HF (65), particularly in HFpEF (60, 66, 67). As the global proteomics dataset showed a decreased in SIRT3 in the LV from HFpEF mice vs. Sham (P = 0.000914), we thus performed additional immunoblot analysis to validate these findings. Indeed, SIRT3 expression was significantly decreased in the LV from HFpEF mice vs. Sham (0.8 ± 0.0 vs. 1.0 ± 0.0; P < 0.001; Figure 6A).

FIGURE 6

Discussion

Heart failure with preserved ejection fraction is a complex disease involving several sub-phenotypes within a heterogeneous HFpEF syndrome (10, 13, 68). Of all the comorbidities in HFpEF, hypertension remains the most common, and is implicated in both the pathogenesis and the prognosis of the disease (12, 29). However, the exact biological mechanisms that underlie hypertension associated HFpEF remain largely unclear. In this study, we investigated the proteomic and phosphoproteomics profile underlying HFpEF in a clinically relevant murine model of hypertension associated HFpEF. The SAUNA model of HFpEF model fulfils the criteria for a “high probability of HFpEF” based on HFA-PEFF diagnostic algorithm for human HFpEF (58, 69).

In the present study, extensive proteomics and phosphoproteomics analysis permitted in-depth screening of the changes in protein expression, post-translational modifications (i.e., phosphorylation), and pathway alterations in HFpEF. These included but were not limited to: (I) changes in cardiac metabolism, where the predominant components were the mitochondrial metabolic processes and mitochondrial dysfunction; (II) alteration in cardiac contractile function-related proteins; (III) overexpression of pathways related to immune modulation; and (IV) a significant decrease in SIRT3 expression, that was validated by immunoblotting.

We found marked changes in signatures of protein expression related to mitochondrial function and oxidative metabolism of energy substrates in HFpEF. There was a significant decrease in targets related to mitochondrial substrate oxidation, suggesting that cardiac mitochondrial metabolic function is impaired in HFpEF. Interestingly, there was an upregulation of the ketone bodies transporter SLC16A1 in the LV of HFpEF, but this was not accompanied by comparable changes in ketone metabolism enzymes. Although not investigated in this study, these findings may contribute to the metabolic impairment seen in HFpEF by increasing the transport of ketone bodies into the mitochondria, but without a compensatory catabolic response. We hypothesize that this mismatch in mitochondrial substrate intake and utilization results in mitochondrial ketone bodies accumulation which may detrimentally affect cardiac function (70). Ketone bodies are thought to be a relevant energy source in both preclinical HFrEF models (71) and advance HFrEF patients (72). Additionally, it has been shown that HFpEF patients have significantly higher circulating ketone levels than HFrEF patients (73) suggesting that some of the beneficial effects of SGLT2 inhibitors in HFpEF may be due to enhanced ketone bodies availability and cardiac utilization (74–76), a process known as “thrifty substrate/fuel hypothesis” (77). We also observed decreased OXCT1 (aka SCOT, succinyl-CoA:3-ketoacid CoA transferase) expression in HFpEF hearts (Figure 3). OXCT1 allows cells to utilize energy stored in ketone bodies thus its decrease in HFpEF hearts supports a role for ketone body cardiac metabolism. Similarly, others have shown worse HF in pre-clinical models with cardio-specific deletion of OXCT1 (78).

Proteomic evaluation of PTMs is essential to understand the function of many proteins in physiological and pathophysiological settings. PTMs are regulators of protein structure and function and, in the heart the predominant PTM is phosphorylation, followed by acetylation (79), and it is also recognized that many proteins are regulated by phosphorylation independently of their expression (80).

Titin is a major cardiac protein regulated by phosphorylation and facilitates myocardial passive tension by conditioning cardiomyocyte-derived stiffness (81). Titin regulates cardiomyocyte stiffness both at the transcriptional and post-transcriptional level. At the transcriptional level, titin shifts from its compliant isoform N2BA toward its stiff isoform N2B, which contributes to the impaired diastolic function that is seen in HFpEF (33, 64, 82). In the present study, translational and PTMs in titin are apparent in the LV of HFpEF hearts. The stiffer N2B isoform was significantly increased in HFpEF mice. However, it is notably that the N2B isoform is also the predominant isoform expressed in the LV of rodents (83). At the post-transcriptional level, despite comparable global proteomics expression between HFpEF mice and Sham, phosphoproteomics analysis showed that titin was one of the proteins with the greatest alterations in phosphorylation in HFpEF mice. Similar to a Dahl salt-sensitive rat study (84), in these SAUNA HFpEF mice most of the significantly hyper-phosphorylated titin residues were located at the Z-disk binding region of the titin protein (Figure 5C). Interestingly, it has been suggested that titin may be part of a Z-disk macromolecular machinery acting as a node for hypertrophic signaling (85). As such, our findings that the myofilament and myofilament-associated proteins viz. ACTA1, MYH7, MHY9, TPM1, and MYOZ2 were differentially expressed in both global and phosphoproteomics dataset, support the premise that alterations in sarcomeric and myofilament regulating proteins play a central role in HFpEF. Of these proteins, MYH7, TPM1 and MYOZ2 are known to be important in hypertrophic cardiomyopathy (86–88), and may play a similar role in HFpEF. However, although their function in muscle contraction is well known (89, 90), their role in cardiac hypertrophy and adverse cardiac remodeling remains elusive. It has been hypothesized that changes in cardiac architecture may be a compensatory response that eventually fails, resulting in a re-induction of fetal genes, fibrosis replacing necrotic and apoptotic cardiac cells, and a shift in metabolic substrates (91). However, additional studies are warranted to identify the precise role these myofilament-associated proteins play in HFpEF.

Because of its size, titin has more phosphorylation sites than other smaller proteins and hundreds of phosphorylation sites have been predicted based on proteomic analysis (83), Z-disk Similarly, multiple kinases are also involved in titin phosphorylation (92), representing more opportunities for the regulation of the cardiomyocyte structure and function. However, the effect that a specific phosphorylation pattern has on the function of titin is largely dependent on the specific structural domain which is modified within the protein (93). For example several studies have focused on the “spring-like” I-domain, including the N2bus and PEVK regions, likely due to the mechanically active nature of this specific domain, where phosphorylation may modulate passive and active tension of the sarcomere (85, 92–94). Conversely, the proline-directed kinases, including extracellular signal-regulated kinase-1/-2 (ERK1/2) and cyclin-dependent protein kinase-2 (Cdc2) were able to regulate the phosphorylation status of non-extensible Z-disk (95, 96) and C-terminal (M-band) (97) regions (98). Although additional studies using site-specific methods are needed (92), it has been suggested that changes in the phosphorylation status of these regions may have an important function, not only during developmental stages, but also regulating the binding of titin to other and M-band proteins, as well as the assembly and turnover of these binding partners (99, 100).

In addition to phosphorylation, HFpEF is also associated with hyperacetylation of mitochondrial proteins in the myocardium (101, 102). In the mitochondria, the acetylation state of key enzymes involved in mitochondrial metabolism, oxidative stress defense and mitochondrial dynamics is regulated by the mitochondrial, NAD-dependent protein deacetylase SIRT3 (103–106). SIRT3 interacts with at least 84 mitochondrial proteins involved in many aspects of mitochondrial biology, such as maintaining mitochondrial integrity and function (67, 106). In the present study, the global proteomics data set showed decreased expression of SIRT3 in the LV of HFpEF mice, which was also confirmed by immunoblotting. Others have shown that reduced SIRT3 expression is related to reduced NAD+ bioavailability in HFpEF, and that cardiomyocyte specific SIRT3 knockout mice developed worse diastolic dysfunction in HFpEF (60). Additional studies using whole-body knockout or transgenic mice similarly showed that SIRT3 is required to maintain cardiac contractile function under pro-hypertrophic or ischemic stress (107–110). A recent study showed that a deficit of cardiac NAD+ exists not only pre-clinical HFpEF models but also in patients with HFpEF, and that increasing NAD+ levels with nicotinamide improved diastolic dysfunction (111). The authors hypothesized the beneficial effects were mediated, partly by increasing deacetylation of proteins that regulate the mechano-elastic properties of cardiac myocytes such as titin and sarco/endoplasmic reticulum Ca²⁺ ATPase 2a (SERCA2a). Although not investigated in the present study, SIRT3 may also play a role in cardiomyocyte stiffness and impaired diastolic function in HFpEF, possibly by titin acetylation (59). SIRT 3 may be a future target since, compared to younger subjects, exercise increases SIRT3 protein expression in muscle, which is decreased in older sedentary individuals (112). Interestingly, of the 7 mammalian sirtuins described, SIRT3 is the only analog whose increased expression associates with longevity in humans (113–115). Since HFpEF is highly associated with aging, and exercise training is effective in improving the quality of life in HFpEF patients (116), future studies are warranted to explore the role of SIRT3 expression in muscle in patients with HFpEF.

In conclusion, untargeted proteomics have demonstrated a key role of protein PTMs in metabolism, cell preservation and sarcomere function in the heart (117). In the present study marked proteomics and phosphoproteomics changes occurred in the heart in HFpEF mice which were related to altered mitochondrial metabolism and sarcomere contractility. It is possible that SIRT3 plays a pivotal role in HFpEF, by regulating mitochondrial metabolism and titin stiffness but this requires further study (Figure 6B).

References

• 1.

  MetraMTeerlinkJR. Heart failure.Lancet. (2017) 390:1981–95. 10.1016/S0140-6736(17)31071-1

  • CrossRef
  • Google Scholar

• 2.

  UdelsonJELewisGDShahSJZileMRRedfieldMMBurnettJJret alEffect of praliciguat on peak rate of oxygen consumption in patients with heart failure with preserved ejection fraction: the CAPACITY HFpEF randomized clinical trial.JAMA. (2020) 324:1522–31. 10.1001/jama.2020.16641

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  ArmstrongPWLamCSPAnstromKJEzekowitzJHernandezAFO’ConnorCMet alEffect of vericiguat vs placebo on quality of life in patients with heart failure and preserved ejection fraction: the VITALITY-HFpEF randomized clinical trial.JAMA. (2020) 324:1512–21. 10.1001/jama.2020.15922

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  AnkerSDButlerJFilippatosGFerreiraJPBocchiEBöhmMet alEmpagliflozin in heart failure with a preserved ejection fraction.N Engl J Med. (2021) 385:1451–61.

  • Google Scholar

• 5.

  SteinbergBAZhaoXHeidenreichPAPetersonEDBhattDLCannonCPet alTrends in patients hospitalized with heart failure and preserved left ventricular ejection fraction: prevalence, therapies, and outcomes.Circulation. (2012) 126:65–75. 10.1161/CIRCULATIONAHA.111.080770

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  ZileMRGottdienerJSHetzelSJMcMurrayJJKomajdaMMcKelvieRet alPrevalence and significance of alterations in cardiac structure and function in patients with heart failure and a preserved ejection fraction.Circulation. (2011) 124:2491–501. 10.1161/CIRCULATIONAHA.110.011031

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  PaulusWJTschöpeC. A novel paradigm for heart failure with preserved ejection fraction: comorbidities drive myocardial dysfunction and remodeling through coronary microvascular endothelial inflammation.J Am Coll Cardiol. (2013) 62:263–71. 10.1016/j.jacc.2013.02.092

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8.

  FerrariRBöhmMClelandJGPaulusWJPieskeBRapezziCet alHeart failure with preserved ejection fraction: uncertainties and dilemmas.Eur J Heart Fail. (2015) 17:665–71. 10.1002/ejhf.304

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  BorlaugBA. Heart failure with preserved and reduced ejection fraction: different risk profiles for different diseases.Eur Heart J. (2013) 34:1393–5. 10.1093/eurheartj/eht117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  ShahSJKitzmanDWBorlaugBAvan HeerebeekLZileMRKassDAet alPhenotype-specific treatment of heart failure with preserved ejection fraction: a multiorgan roadmap.Circulation. (2016) 134:73–90. 10.1161/CIRCULATIONAHA.116.021884

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  van HeerebeekLPaulusWJ. Understanding heart failure with preserved ejection fraction: where are we today?Neth Heart J. (2016) 24:227–36. 10.1007/s12471-016-0810-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  TrompJClaggettBLLiuJJacksonAMJhundPSKøberLet alGlobal differences in heart failure with preserved ejection fraction: the PARAGON-HF trial.Circ Heart Fail. (2021) 14:e007901. 10.1161/CIRCHEARTFAILURE.120.007901

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  ShahSJKatzDHSelvarajSBurkeMAYancyCWGheorghiadeMet alPhenomapping for novel classification of heart failure with preserved ejection fraction.Circulation. (2015) 131:269–79. 10.1161/CIRCULATIONAHA.114.010637

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  LundbyAOlsenJV. Phosphoproteomics taken to heart.Cell Cycle. (2013) 12:2707–8. 10.4161/cc.25883

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  ArrellDKNeverovaIVan EykJE. Cardiovascular proteomics: evolution and potential.Circ Res. (2001) 88:763–73. 10.1161/hh0801.090193

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  FarmakisDPapingiotisGParissisJFilippatosG. Ups and downs in heart failure: the case of proteomics.Eur J Heart Fail. (2018) 20:63–6. 10.1002/ejhf.1065

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  SteinerSWitzmannFA. Proteomics: applications and opportunities in preclinical drug development.Electrophoresis. (2000) 21:2099–104. 10.1002/1522-2683(20000601)21:11<2099::AID-ELPS2099>3.0.CO;2-N

  • CrossRef
  • Google Scholar

• 18.

  SharmaPCosmeJGramoliniAO. Recent advances in cardiovascular proteomics.J Proteomics. (2013) 81:3–14. 10.1016/j.jprot.2012.10.026

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19.

  CorbettJMWhyHJWheelerCHRichardsonPJArchardLCYacoubMHet alCardiac protein abnormalities in dilated cardiomyopathy detected by two-dimensional polyacrylamide gel electrophoresis.Electrophoresis. (1998) 19:2031–42. 10.1002/elps.1150191123

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  DuranMCMasSMartin-VenturaJLMeilhacOMichelJBGallego-DelgadoJet alProteomic analysis of human vessels: application to atherosclerotic plaques.Proteomics. (2003) 3:973–8. 10.1002/pmic.200300389

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  LaiLPLinJLLinCSYehHMTsayYGLeeCFet alFunctional genomic study on atrial fibrillation using cDNA microarray and two-dimensional protein electrophoresis techniques and identification of the myosin regulatory light chain isoform reprogramming in atrial fibrillation.J Cardiovasc Electrophysiol. (2004) 15:214–23. 10.1046/j.1540-8167.2004.03423.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  KuzmanovUGuoHBuchsbaumDCosmeJAbbasiCIsserlinRet alGlobal phosphoproteomic profiling reveals perturbed signaling in a mouse model of dilated cardiomyopathy.Proc Natl Acad Sci U.S.A. (2016) 113:12592–7. 10.1073/pnas.1606444113

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23.

  LuDXiaYChenZChenAWuYJiaJet alCardiac proteome profiling in ischemic and dilated cardiomyopathy mouse models.Front Physiol. (2019) 10:750. 10.3389/fphys.2019.00750

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24.

  BrioschiMPolvaniGFrattoPParolariAAgostoniPTremoliEet alRedox proteomics identification of oxidatively modified myocardial proteins in human heart failure: implications for protein function.PLoS One. (2012) 7:e35841. 10.1371/journal.pone.0035841

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25.

  LiWRongRZhaoSZhuXZhangKXiongXet alProteomic analysis of metabolic, cytoskeletal and stress response proteins in human heart failure.J Cell Mol Med. (2012) 16:59–71. 10.1111/j.1582-4934.2011.01336.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26.

  Roselló-LletíEAlonsoJCortésRAlmenarLMartínez-DolzLSánchez-LázaroIet alCardiac protein changes in ischaemic and dilated cardiomyopathy: a proteomic study of human left ventricular tissue.J Cell Mol Med. (2012) 16:2471–86. 10.1111/j.1582-4934.2012.01565.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27.

  DeAgueroJLMcKownENZhangLKeirseyJFischerEGSamediVGet alAltered protein levels in the isolated extracellular matrix of failing human hearts with dilated cardiomyopathy.Cardiovasc Pathol. (2017) 26:12–20. 10.1016/j.carpath.2016.10.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28.

  AdamoLYuJRocha-ResendeCJavaheriAHeadRDMannDL. Proteomic signatures of heart failure in relation to left ventricular ejection fraction.J Am Coll Cardiol. (2020) 76:1982–94. 10.1016/j.jacc.2020.08.061

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29.

  ShahSJBorlaugBAKitzmanDWMcCullochADBlaxallBCAgarwalRet alResearch priorities for heart failure with preserved ejection fraction: national heart, lung, and blood institute working group summary.Circulation. (2020) 141:1001–26. 10.1161/CIRCULATIONAHA.119.041886

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  Wilson RichardMDe Silva DeepaSSatoKIzumiyaYSamF. Effects of fixed-dose isosorbide dinitrate/hydralazine on diastolic function and exercise capacity in hypertension-induced diastolic heart failure.Hypertension. (2009) 54:583–90. 10.1161/HYPERTENSIONAHA.109.134932

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31.

  GarciaAGWilsonRMHeoJMurthyNRBaidSOuchiNet alInterferon-γ ablation exacerbates myocardial hypertrophy in diastolic heart failure.Am J Physiol Heart Circ Physiol. (2012) 303:H587–96. 10.1152/ajpheart.00298.2012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32.

  TanakaKValero-MuñozMWilsonRMEssickEEFowlerCTNakamuraKet alFollistatin like 1 regulates hypertrophy in heart failure with preserved ejection fraction.JACC Basic Transl Sci. (2016) 1:207–21. 10.1016/j.jacbts.2016.04.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33.

  Valero-MunozMLiSWilsonRMBoldbaatarBIglarzMSamF. Dual endothelin-A/endothelin-B receptor blockade and cardiac remodeling in heart failure with preserved ejection fraction.Circ Heart Fail. (2016) 9:e003381. 10.1161/CIRCHEARTFAILURE.116.003381

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34.

  Valero-MuñozMBackmanWSamF. Murine models of heart failure with preserved ejection fraction: a “fishing expedition”.JACC Basic Transl Sci. (2017) 2:770–89. 10.1016/j.jacbts.2017.07.013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  HulsmansMSagerHBRohJDValero-MuñozMHoustisNEIwamotoYet alCardiac macrophages promote diastolic dysfunction.J Exp Med. (2018) 215:423–40. 10.1084/jem.20171274

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  YangHJKongBShuaiWZhangJJHuangH. Knockout of MD1 contributes to sympathetic hyperactivity and exacerbates ventricular arrhythmias following heart failure with preserved ejection fraction via NLRP3 inflammasome activation.Exp Physiol. (2020) 105:966–78. 10.1113/EP088390

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  YoonSKimMLeeHKangGBediKMarguliesKBet alS-nitrosylation of histone deacetylase 2 by neuronal nitric oxide synthase as a mechanism of diastolic dysfunction.Circulation. (2021) 143:1912–25. 10.1161/CIRCULATIONAHA.119.043578

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  KiatchoosakunSRestivoJKirkpatrickDHoitBD. Assessment of left ventricular mass in mice: comparison between two-dimensional and m-mode echocardiography.Echocardiography. (2002) 19:199–205. 10.1046/j.1540-8175.2002.00199.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39.

  ZhouYQFosterFSParkesRAdamsonSL. Developmental changes in left and right ventricular diastolic filling patterns in mice.Am J Physiol Heart Circ Physiol. (2003) 285:H1563–75. 10.1152/ajpheart.00384.2003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  GuoHIsserlinREmiliABurnistonJG. Exercise-responsive phosphoproteins in the heart.J Mol Cell Cardiol. (2017) 111:61–8. 10.1016/j.yjmcc.2017.08.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41.

  KimSRabhiNBlumBCHekmanRWynneKEmiliAet alTriphenyl phosphate is a selective PPARγ modulator that does not induce brite adipogenesis in vitro and in vivo.Arch Toxicol. (2020) 94:3087–103.

  • Google Scholar

• 42.

  DostAFMMoyeALVedaieMTranLMFungEHeinzeDet alOrganoids model transcriptional hallmarks of oncogenic KRAS activation in lung epithelial progenitor cells.Cell Stem Cell. (2020) 27:663–78.e8. 10.1016/j.stem.2020.07.022

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  CantinGTShockTRParkSKMadhaniHDYatesJRIIIOptimizing TiO2-based phosphopeptide enrichment for automated multidimensional liquid chromatography coupled to tandem mass spectrometry.Anal Chem. (2007) 79:4666–73. 10.1021/ac0618730

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44.

  GuoHIsserlinRChenXWangWPhanseSZandstraPWet alIntegrative network analysis of signaling in human CD34(+) hematopoietic progenitor cells by global phosphoproteomic profiling using TiO2 enrichment combined with 2D LC-MS/MS and pathway mapping.Proteomics. (2013) 13:1325–33. 10.1002/pmic.201200369

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  BlumBCEmiliAOmics notebook: robust, reproducible, and flexible automated multi-omics exploratory analysis and reporting. Bioinform. Adv. (2021) 1:1–4.

  • Google Scholar

• 46.

  RitchieMEPhipsonBWuDHuYLawCWShiWet allimma powers differential expression analyses for RNA-sequencing and microarray studies.Nucleic Acids Res. (2015) 43:e47. 10.1093/nar/gkv007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47.

  SergushichevAA. An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation.bioRxiv [Preprint]. (2016).

  • Google Scholar

• 48.

  KorotkevichGSukhovVBudinNShpakBArtyomovMNSergushichevA. Fast gene set enrichment analysis.bioRxiv [Preprint]. (2021). 10.1101/060012

  • CrossRef
  • Google Scholar

• 49.

  SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAet alGene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proc Natl Acad Sci U.S.A. (2005) 102:15545–50. 10.1073/pnas.0506580102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50.

  MericoDIsserlinRStuekerOEmiliABaderGD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation.PLoS One. (2010) 5:e13984. 10.1371/journal.pone.0013984

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51.

  WarrenCMKrzesinskiPRGreaserML. Vertical agarose gel electrophoresis and electroblotting of high-molecular-weight proteins.Electrophoresis. (2003) 24:1695–702. 10.1002/elps.200305392

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52.

  NeuhoffVAroldNTaubeDEhrhardtW. Improved staining of proteins in polyacrylamide gels including isoelectric focusing gels with clear background at nanogram sensitivity using Coomassie brilliant blue G-250 and R-250.Electrophoresis. (1988) 9:255–62. 10.1002/elps.1150090603

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53.

  van der PijlRJHudsonBGranzier-NakajimaTLiFKnottnerusAMSmithJet alDeleting titin’s C-terminal PEVK exons increases passive stiffness, alters splicing, and induces cross-sectional and longitudinal hypertrophy in skeletal muscle.Front Physiol. (2020) 11:494. 10.3389/fphys.2020.00494

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54.

  SawELPearsonJTSchwenkeDOMunasinghePETsuchimochiHRawalSet alActivation of the cardiac non-neuronal cholinergic system prevents the development of diabetes-associated cardiovascular complications.Cardiovasc Diabetol. (2021) 20:50. 10.1186/s12933-021-01231-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55.

  R Core Team.R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing (2020).

  • Google Scholar

• 56.

  Sanders-Van WijkSBarandiarán AizpuruaABrunner-La RoccaHPHenkensMWeertsJKnackstedtCet alThe HFA-PEFF and H2FPEF scores largely disagree in classifying patients with suspected heart failure with preserved.Eur J Heart Fail. (2021) 23:838–40. 10.1002/ejhf.2019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57.

  Barandiarán AizpuruaASanders-Van WijkSBrunner-La RoccaHPHenkensMHeymansSBeussink-NelsonLet alValidation of the HFA-PEFF score for the diagnosis of heart failure with preserved ejection fraction.Eur J Heart Fail. (2020) 22:413–21. 10.1002/ejhf.1614

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58.

  WithaarCLamCSPSchiattarellaGGDe BoerRAMeemsLMG. Heart failure with preserved ejection fraction in humans and mice: embracing clinical complexity in mouse models.Eur Heart J. (2021) 42:4420–30. 10.1093/eurheartj/ehab389

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59.

  ZengHChenJX. Sirtuin 3, endothelial metabolic reprogramming, and heart failure with preserved ejection fraction.J Cardiovasc Pharmacol. (2019) 74:315–23. 10.1097/FJC.0000000000000719

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60.

  TongDSchiattarellaGGJiangNAltamiranoFSzwedaPAElnwasanyAet alNAD(+) repletion reverses heart failure with preserved ejection fraction.Circ Res. (2021) 128:1629–41. 10.1161/CIRCRESAHA.120.317046

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61.

  ZhangWZhangHYaoWLiLNiuPHuoYet alMorphometric, hemodynamic, and multi-omics analyses in heart failure rats with preserved ejection fraction.Int J Mol Sci. (2020) 21:3362. 10.3390/ijms21093362

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62.

  VoelkelTLinkeWA. Conformation-regulated mechanosensory control via titin domains in cardiac muscle.Pflugers Arch. (2011) 462:143–54. 10.1007/s00424-011-0938-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63.

  KötterSAndresenCKrügerM. Titin: central player of hypertrophic signaling and sarcomeric protein quality control.Biol Chem. (2014) 395:1341–52. 10.1515/hsz-2014-0178

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64.

  ZileMRBaicuCFIkonomidisJSStroudRENietertPJBradshawADet alMyocardial stiffness in patients with heart failure and a preserved ejection fraction: contributions of collagen and titin.Circulation. (2015) 131:1247–59. 10.1161/CIRCULATIONAHA.114.013215

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65.

  ChenJChenSZhangBLiuJ. SIRT3 as a potential therapeutic target for heart failure.Pharmacol Res. (2021) 165:105432. 10.1016/j.phrs.2021.105432

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66.

  LaiYCTabimaDMDubeJJHughanKSVanderpoolRRGoncharovDAet alSIRT3-AMP-activated protein kinase activation by nitrite and metformin improves hyperglycemia and normalizes pulmonary hypertension associated with heart failure with preserved ejection fraction.Circulation. (2016) 133:717–31. 10.1161/CIRCULATIONAHA.115.018935

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67.

  Del CampoAPerezGCastroPFParraVVerdejoHE. Mitochondrial function, dynamics and quality control in the pathophysiology of HFpEF.Biochim Biophys Acta Mol Basis Dis. (2021) 1867:166208. 10.1016/j.bbadis.2021.166208

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68.

  ShahSJ. Precision medicine for heart failure with preserved ejection fraction: an overview.J Cardiovasc Transl Res. (2017) 10:233–44. 10.1007/s12265-017-9756-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69.

  PieskeBTschöpeCDe BoerRAFraserAGAnkerSDDonalEet alHow to diagnose heart failure with preserved ejection fraction: the HFA–PEFF diagnostic algorithm: a consensus recommendation from the Heart Failure Association (HFA) of the European Society of Cardiology (ESC).Eur Heart J. (2019) 40:3297–317. 10.1093/eurheartj/ehz641

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70.

  BerteroESequeiraVMaackC. Hungry hearts.Circ Heart Fail. (2018) 11:e005642. 10.1161/CIRCHEARTFAILURE.118.005642

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71.

  AubertGMartinOJHortonJLLaiLVegaRBLeoneTCet alThe failing heart relies on ketone bodies as a fuel.Circulation. (2016) 133:698–705. 10.1161/CIRCULATIONAHA.115.017355

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72.

  BediKCJrSnyderNWBrandimartoJAzizMMesarosCWorthAJet alEvidence for intramyocardial disruption of lipid metabolism and increased myocardial ketone utilization in advanced human heart failure.Circulation. (2016) 133:706–16. 10.1161/CIRCULATIONAHA.115.017545

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73.

  ZordokyBNSungMMEzekowitzJMandalRHanBBjorndahlTCet alMetabolomic fingerprint of heart failure with preserved ejection fraction.PLoS One. (2015) 10:e0124844. 10.1371/journal.pone.0124844

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74.

  Echouffo-TcheuguiJBLewseySCWeissRG. SGLT2 inhibitors: further evidence for heart failure with preserved ejection fraction as a metabolic disease?J Clin Invest. (2021) 131:e156309. 10.1172/JCI156309

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75.

  MurashigeDJangCNeinastMEdwardsJJCowanAHymanMCet alComprehensive quantification of fuel use by the failing and nonfailing human heart.Science. (2020) 370:364–8. 10.1126/science.abc8861

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76.

  YuristaSRRosenzweigANguyenCT. Ketone bodies: universal cardiac response to stress?J Am Coll Cardiol. (2021) 78:1433–6. 10.1016/j.jacc.2021.08.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77.

  FerranniniEMarkMMayouxE. CV protection in the EMPA-REG OUTCOME trial: a “thrifty substrate” hypothesis.Diabetes Care. (2016) 39:1108–14. 10.2337/dc16-0330

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78.

  SchugarRCMollARAndré d’AvignonDWeinheimerCJKovacsACrawfordPA. Cardiomyocyte-specific deficiency of ketone body metabolism promotes accelerated pathological remodeling.Mol Metab. (2014) 3:754–69. 10.1016/j.molmet.2014.07.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79.

  Van EykJE. Overview: the maturing of proteomics in cardiovascular research.Circ Res. (2011) 108:490–8.

  • Google Scholar

• 80.

  HuttlinELJedrychowskiMPEliasJEGoswamiTRadRBeausoleilSAet alA tissue-specific atlas of mouse protein phosphorylation and expression.Cell. (2010) 143:1174–89.

  • Google Scholar

• 81.

  LamCSPVoorsAAde BoerRASolomonSDvan VeldhuisenDJ. Heart failure with preserved ejection fraction: from mechanisms to therapies.Eur Heart J. (2018) 39:2780–92.

  • Google Scholar

• 82.

  BorlaugBAPaulusWJ. Heart failure with preserved ejection fraction: pathophysiology, diagnosis, and treatment.Eur Heart J. (2011) 32: 670–9.

  • Google Scholar

• 83.

  FukudaNGranzierHIshiwataSMorimotoS. Editorial: recent advances on myocardium physiology.Front Physiol. (2021) 12:697852. 10.3389/fphys.2021.697852

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84.

  SoetkampDGalletRParkerSJHolewinskiRVenkatramanVPeckKet alMyofilament phosphorylation in stem cell treated diastolic heart failure.Circ Res. (2021) 129:1125–40. 10.1161/CIRCRESAHA.119.316311

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85.

  LinkeWAKrugerM. The giant protein titin as an integrator of myocyte signaling pathways.Physiology. (2010) 25:186–98. 10.1152/physiol.00005.2010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86.

  FreyNLueddeMKatusHA. Mechanisms of disease: hypertrophic cardiomyopathy.Nat Rev Cardiol. (2011) 9:91–100. 10.1038/nrcardio.2011.159

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87.

  MaronBJMaronMS. Hypertrophic cardiomyopathy.Lancet. (2013) 381:242–55. 10.1016/S0140-6736(12)60397-3

  • CrossRef
  • Google Scholar

• 88.

  LandstromAPAckermanMJ. Beyond the cardiac myofilament: hypertrophic cardiomyopathy- associated mutations in genes that encode calcium-handling proteins.Curr Mol Med. (2012) 12:507–18. 10.2174/156652412800620020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89.

  MarianAJ. Molecular genetic basis of hypertrophic cardiomyopathy.Circ Res. (2021) 128:1533–53. 10.1161/CIRCRESAHA.121.318346

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90.

  SequeiraVNijenkampLLReganJAvan der VeldenJ. The physiological role of cardiac cytoskeleton and its alterations in heart failure.Biochim Biophys Acta. (2014) 1838:700–22. 10.1016/j.bbamem.2013.07.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91.

  HarveyPALeinwandLA. The cell biology of disease: cellular mechanisms of cardiomyopathy.J Cell Biol. (2011) 194:355–65. 10.1083/jcb.201101100

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92.

  KoserFLoescherCLinkeWA. Posttranslational modifications of titin from cardiac muscle: how, where, and what for?FEBS J. (2019) 286:2240–60. 10.1111/febs.14854

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 93.

  TharpCAHaywoodMESbaizeroOTaylorMRGMestroniL. The giant protein titin’s role in cardiomyopathy: genetic, transcriptional, and post-translational modifications of TTN and their contribution to cardiac disease.Front Physiol. (2019) 10:1436. 10.3389/fphys.2019.01436

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 94.

  LeWinterMMGranzierHL. Cardiac titin and heart disease.J Cardiovasc Pharmacol. (2014) 63:207–12. 10.1097/FJC.0000000000000007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 95.

  GautelMGouldingDBullardBWeberKFurstDO. The central Z-disk region of titin is assembled from a novel repeat in variable copy numbers.J Cell Sci. (1996) 109(Pt 11):2747–54. 10.1242/jcs.109.11.2747

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 96.

  SebestyenMGWolffJAGreaserML. Characterization of a 5.4 kb cDNA fragment from the Z-line region of rabbit cardiac titin reveals phosphorylation sites for proline-directed kinases.J Cell Sci. (1995) 108(Pt 9):3029–37. 10.1242/jcs.108.9.3029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 97.

  GautelMLeonardKLabeitS. Phosphorylation of KSP motifs in the C-terminal region of titin in differentiating myoblasts.EMBO J. (1993) 12:3827–34. 10.1002/j.1460-2075.1993.tb06061.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 98.

  LoescherCMHobbachAJLinkeWA. Titin (TTN): from molecule to modifications, mechanics and medical significance.Cardiovasc Res. (2021) 1–6. 10.1093/cvr/cvab328[Epub ahead of print].

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99.

  HidalgoCGranzierH. Tuning the molecular giant titin through phosphorylation: role in health and disease.Trends Cardiovasc Med. (2013) 23:165–71. 10.1016/j.tcm.2012.10.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 100.

  HamdaniNHerwigMLinkeWA. Tampering with springs: phosphorylation of titin affecting the mechanical function of cardiomyocytes.Biophys Rev. (2017) 9:225–37. 10.1007/s12551-017-0263-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 101.

  HortonJLMartinOJLaiLRileyNMRichardsALVegaRBet alMitochondrial protein hyperacetylation in the failing heart.JCI Insight. (2016) 2:e84897.

  • Google Scholar

• 102.

  LaiLLeoneTCKellerMPMartinOJBromanATNigroJet alEnergy metabolic reprogramming in the hypertrophied and early stage failing heart: a multisystems approach.Circ Heart Fail. (2014) 7:1022–31. 10.1161/CIRCHEARTFAILURE.114.001469

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 103.

  JinLGalonekHIsraelianKChoyWMorrisonMXiaYet alBiochemical characterization, localization, and tissue distribution of the longer form of mouse SIRT3.Protein Sci. (2009) 18:514–25. 10.1002/pro.50

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 104.

  RardinMJNewmanJCHeldJMCusackMPSorensenDJLiBet alLabel-free quantitative proteomics of the lysine acetylome in mitochondria identifies substrates of SIRT3 in metabolic pathways.Proc Natl Acad Sci U.S.A. (2013) 110:6601–6. 10.1073/pnas.1302961110

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 105.

  SamantSAZhangHJHongZPillaiVBSundaresanNRWolfgeherDet alSIRT3 deacetylates and activates OPA1 to regulate mitochondrial dynamics during stress.Mol Cell Biol. (2014) 34:807–19. 10.1128/MCB.01483-13

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 106.

  ZhangJXiangHLiuJChenYHeRRLiuB. Mitochondrial sirtuin 3: new emerging biological function and therapeutic target.Theranostics. (2020) 10:8315–42. 10.7150/thno.45922

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 107.

  SundaresanNRGuptaMKimGRajamohanSBIsbatanAGuptaMP. Sirt3 blocks the cardiac hypertrophic response by augmenting Foxo3a-dependent antioxidant defense mechanisms in mice.J Clin Invest. (2009) 119:2758–71. 10.1172/JCI39162

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 108.

  KoentgesCPfeilKSchnickTWieseSDahlbockRCimolaiMCet alSIRT3 deficiency impairs mitochondrial and contractile function in the heart.Basic Res Cardiol. (2015) 110:36. 10.1007/s00395-015-0493-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 109.

  HeXZengHChenJX. Ablation of SIRT3 causes coronary microvascular dysfunction and impairs cardiac recovery post myocardial ischemia.Int J Cardiol. (2016) 215:349–57. 10.1016/j.ijcard.2016.04.092

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 110.

  HeXZengHChenJX. Emerging role of SIRT3 in endothelial metabolism, angiogenesis, and cardiovascular disease.J Cell Physiol. (2019) 234:2252–65. 10.1002/jcp.27200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 111.

  AbdellatifMTrummer-HerbstVKoserFDurandSAdãoRVasques-NóvoaFet alNicotinamide for the treatment of heart failure with preserved ejection fraction.Sci Transl Med. (2021) 13:eabd7064. 10.1126/scitranslmed.abd7064

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 112.

  LanzaIRShortDKShortKRRaghavakaimalSBasuRJoynerMJet alEndurance exercise as a countermeasure for aging.Diabetes. (2008) 57:2933–42. 10.2337/db08-0349

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 113.

  RoseGDatoSAltomareKBellizziDGarastoSGrecoVet alVariability of the SIRT3 gene, human silent information regulator Sir2 homologue, and survivorship in the elderly.Exp Gerontol. (2003) 38:1065–70. 10.1016/s0531-5565(03)00209-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 114.

  BellizziDRoseGCavalcantePCovelloGDatoSDe RangoFet alA novel VNTR enhancer within the SIRT3 gene, a human homologue of SIR2, is associated with survival at oldest ages.Genomics. (2005) 85:258–63. 10.1016/j.ygeno.2004.11.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 115.

  PillaiVBSundaresanNRJeevanandamVGuptaMP. Mitochondrial SIRT3 and heart disease.Cardiovasc Res. (2010) 88:250–6. 10.1093/cvr/cvq250

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 116.

  PandeyAParasharAKumbhaniDAgarwalSGargJKitzmanDet alExercise training in patients with heart failure and preserved ejection fraction: meta-analysis of randomized control trials.Circ Heart Fail. (2015) 8:33–40. 10.1161/CIRCHEARTFAILURE.114.001615

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 117.

  RaghowR. An ‘omics’ perspective on cardiomyopathies and heart failure.Trends Mol Med. (2016) 22:813–27. 10.1016/j.molmed.2016.07.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar