Sex-related differences of early cardiac functional and proteomic alterations in a rat model of myocardial ischemia
Journal of Translational Medicine volume 19, Article number: 507 (2021)
Reduced cardiovascular risk in premenopausal women has been the focus of research in recent decades. Previous hypothesis-driven experiments have highlighted the role of sex hormones on distinct inflammatory responses, mitochondrial proteins, extracellular remodeling and estrogen-mediated cardioprotective signaling pathways related to post-ischemic recovery, which were associated with better cardiac functional outcomes in females. We aimed to investigate the early, sex-specific functional and proteomic changes following myocardial ischemia in an unbiased approach.
Ischemia was induced in male (M-Isch) and female (F-Isch) rats with sc. injection of isoproterenol (85 mg/kg) daily for 2 days, while controls (M-Co, F-Co) received sc. saline solution. At 48 h after the first injection pressure–volume analysis was carried out to assess left ventricular function. FFPE tissue slides were scanned and analyzed digitally, while myocardial proteins were quantified by liquid chromatography–tandem mass spectrometry (LC–MS/MS) using isobaric labeling. Concentrations of circulating steroid hormones were measured with LC–MS/MS. Feature selection (PLS and PLS-DA) was used to examine associations among functional, proteomic and hormonal datasets.
Induction of ischemia resulted in 38% vs 17% mortality in M-Isch and F-Isch respectively. The extent of ischemic damage to surviving rats was comparable between the sexes. Systolic dysfunction was more pronounced in males, while females developed a more severe impairment of diastolic function. 2224 proteins were quantified, with 520 showing sex-specific differential regulation. Our analysis identified transcriptional, cytoskeletal, contractile, and mitochondrial proteins, molecular chaperones and the extracellular matrix as sources of disparity between the sexes. Bioinformatics highlighted possible associations of estrogens and their metabolites with early functional and proteomic alterations.
Our study has highlighted sex-specific alterations in systolic and diastolic function shortly after ischemia, and provided a comprehensive look at the underlying proteomic changes and the influence of estrogens and their metabolites. According to our bioinformatic analysis, inflammatory, mitochondrial, chaperone, cytoskeletal, extracellular and matricellular proteins are major sources of intersex disparity, and may be promising targets for early sex-specific pharmacologic interventions.
Sex has been shown to have a major impact on the pathophysiology, presentation and outcomes of cardiovascular diseases . This phenomenon is particularly apparent in coronary artery diseases (CAD). The occurrence of myocardial infarction (MI) in men is expected earlier and with higher severity compared to age-matched women . The exact mechanisms underlying this discrepancy are still poorly understood, but a scientific consensus has formed around the possible protective role of estrogen hormones [2, 3]. Although estradiol has been investigated with regard to its protective effects under ischemic conditions [4, 5], the effect of estrogen metabolites remains understudied despite their potential to modulate estrogen hormone signaling.
Several studies have been conducted on rodents to gain morphological, functional and molecular insights into sex-related differences in response to myocardial ischemic insults. A majority of these investigations have found higher survival rates as well as faster recovery of contractility and relaxation in females in the first few days after MI (as reviewed by Regitz-Zagrosek and Kararigas ). Furthermore, myocardial remodeling following MI has been associated with cavity dilation, eccentric hypertrophy and scar thinning in males, while females exhibited concentric hypertrophy with maintained cavity dimensions and stable scar formation .
Most of the data currently available in the literature on sex-dependent cardiac functional outcomes shortly after ischemic periods come from ex vivo experiments performed on Langendorff perfused hearts [7, 8] or echocardiographic measurements . The major drawback of the former approach is that the excised heart might not accurately represent in vivo cardiac function, while the disadvantage of the latter method is that echocardiographic functional parameters are largely dependent on loading conditions (changes in preload or afterload). The gold standard method for assessment of left-ventricular function remains in vivo pressure–volume (P–V) analysis, which offers highly sensitive, load-independent indices of systolic and diastolic function measured under physiologic conditions . Therefore, P–V analysis is uniquely suited to evaluate even subtle sex-specific alterations in cardiac function following ischemia.
Mechanistic insights into the molecular background of observed sex differences have been gained from multiple studies performed on a great variety of animal models and species, as well as on cultured cardiomyocytes. Mitochondrial function, lipid and carbohydrate metabolism, inflammation, regulation of cell death and fibrosis are referred to as crucial sources of sex-dependent distinction [3, 11,12,13]. In each of these biological processes, proteins play a pivotal role. Thus proteomic characterization of sex-related differences in cardiac adaptation after ischemia is indispensable for a complete understanding of underlying pathomechanisms. Traditional 2D gel-based proteomic methods are still the most abundant source of information due to their availability . However, liquid chromatography coupled mass spectrometry (LC–MS/MS) based proteomics is a rapidly advancing field offering substantially higher proteome coverage than traditional methods, while isobaric labeling molecules (iTRAQ, TMT), allow for drastically increased throughput .
Isoproterenol is a non-selective beta-adrenergic agonist that is known to have a dose-dependent effect on the heart. In moderate amounts, it leads to an immediate, but short-term increase in cardiac performance due to its positive chronotropic, dromotropic and inotropic properties, while chronic application induces ventricular hypertrophy and cardiac dysfunction. Grant et al. recently reported a lack of sexual dimorphism in both short-term and long-term changes in cardiac function in response to moderate doses of isoproterenol . In high doses, however, it causes diffuse endomyocardial infarction as a result of overburdening the reserves of the myocardium, upsetting the delicate balance of oxygen delivery and consumption , which makes it ideal for investigation of sex-related differences in post-ischemic adaptation.
In this study, we aimed to provide a comprehensive and unbiased characterization of sex-specific early cardiac functional and proteomic alterations in a rat model of isoproterenol-induced myocardial ischemia. Furthermore, we investigated the steroid hormone profiles of ischemic female rats and conducted a detailed bioinformatic analysis to establish possible associations among functional outcomes, proteomic and hormonal profiles.
Materials and methods
Our investigation was performed as approved by the Ethics Committee of Hungary for Animal Experimentation, Semmelweis University, Budapest (PEI/001/2374-4/2015) conforming to European Directive 2010/63/EU and the Guide for the Care and Use of Laboratory Animals used by the US National Institutes of Health (Publication No. 85-23, revised 1996). Our study has been interpreted according to the Animals in Research: Reporting In Vivo Experiments Guidelines. The animals were kept in standard conditions (22 ± 2 °C with 12-h light/dark cycles) and were allowed access to laboratory rat diet and water ad libitum during the whole experiment. All animals received humane care.
After acclimatization, male and female Wistar rats (180 to 220 g) were assigned to matched control (Co) and ischemic (Isch) groups (Fig. 1A):
Male control (M-Co, n = 8)
Male ischemic (M-Isch, n = 20)
Female control (F-Co, n = 8)
Female ischemic (F-Isch, n = 17)
All animals underwent body weight measurement (BW) and were briefly anaesthetized using 5% isoflurane. Rats from the ischemic groups were injected with isoproterenol (ISO, 85 mg/kg, subcutaneous, Sigma-Aldrich) daily for 2 consecutive days to induce diffuse subendocardial ischemia as previously described , while control rats received equivalent volumes of subcutaneous saline (NaCl 0.9%) injections at the same timepoints. Electrocardiograms were recorded 2 min after injections to assess the acute effect of the drug by a computer-based data acquisition system (PowerLab 16/30 and LabChart Pro software v7; AD Instruments, Colorado Spring, US).
At 48 h after the first injection, invasive left ventricular (LV) pressure–volume (P–V) analysis was performed according to a previously described protocol . In brief, anesthesia was induced with 5% isoflurane and maintained with 1.5% for the duration of the measurements. Animals were tracheotomized, intubated and placed on an automatic heating pad (equipped with a rectal probe) to maintain their body temperature at 37 °C. Fluid administration during the procedure was carried out using a polyethylene catheter inserted into the left external jugular vein. A 2F microtip pressure microcatheter (SPR-838, Millar Instruments, Houston, TX, USA) was advanced into the ascending aorta through the right carotid artery. After 5 min of stabilization, arterial pressure was recorded with a sampling rate of 1000/s. The catheter was then further advanced through the aortic valve to record LV P–V loops. A dedicated P–V analysis program (PVAN, Millar Instruments) was used to calculate the following parameters: systolic arterial blood pressure (SABP), diastolic arterial blood pressure (DABP), heart rate (HR), LV end-systolic pressure (LVESP), LV end-diastolic pressure (LVEDP), LV end-diastolic volume (LVEDV), LV end-systolic volume (LVESV), stroke volume (SV), cardiac output (CO), ejection fraction (EF), stroke work (SW), maximal rate of rise in LV pressure (dP/dtmax), maximal rate of decrease in LV pressure (dP/dtmin), time-constant of LV pressure decay Tau (τ, Glantz). Load-independent parameters of systolic and diastolic function were determined based on P–V loops recorded under conditions of gradually decreasing preload achieved by transient occlusion of the inferior vena cava. As indices of contractility preload recruitable stroke work (PRSW) and maximal slope of dP/dtmax – EDV relationship (dP/dtmax-EDV), while as a measure of LV diastolic stiffness the slope of the LV end-diastolic pressure–volume relationship (EDPVR) were determined. Parallel conductance was measured and volume calibrations were performed as previously described . After the invasive hemodynamic measurements, rats were euthanized by exsanguination. The hearts were perfused with oxygenated ringer solution and excised. Heart weight (HW) and tibia length (TL) were measured.
Histology and immunohistochemistry
Mid-papillary slices of the left ventricle were excised from the explanted hearts and were stored either by snap freezing or by formalin fixation (24 h in 4% buffered paraformaldehyde solution). From the formalin-fixed paraffin-embedded tissue, 5 μm thick slices were sectioned and stained with hematoxylin and eosin (H&E) and picrosirius red. Sections were then digitalized using a whole slide scanner (Ventana DP 200, Roche). The percentage of areas with intact myofibers to the total myocardial area (an indirect measure of the extent of necrosis) was quantified with QuPath 0.2.3  using H&E color deconvolution following the principles set by van Putten et al. . The percentage of picrosirius red stained area was quantified by thresholding in QuPath.
Cryopreserved myocardial samples were homogenized in 0.1% Rapigest and 100 mM HEPES (pH 7.5) solution. Reductive alkylation was performed using 5 mM tris(2-carboxyethyl)phosphine-hydrochloride (TCEP, Sigma-Aldrich) and 20 mM iodoacetamide (Sigma-Aldrich). Proteins were digested overnight by Trypsin (Worthington, Lakewood, NJ, sequencing grade), and the resulting peptides were desalted on PreOmics columns (PreOmics, Bavaria, Germany). Further steps were performed as described previously. Briefly, peptides were labeled with TMT11plex isobaric label reagents (Thermo Fisher Scientific) and combined into four batches including a cohort-wide normalization channel. Reverse-phase prefractionation (pH = 10) was performed on an XBridge C18 column, 150 mm × 1 mm column containing 3.5 µm particles (Waters) inserted into an Agilent 1100 high performance liquid chromatography system (HPLC). Fractions were analyzed on a Q-Exactive Plus (Thermo Scientific, Bremen, Germany) operating in a data-dependent acquisition (DDA) mode. Mass spectra were analyzed using MaxQuant version 18.104.22.168  with the Uniprot rat database downloaded in November 2020. Quantified peptide intensities were then summarized by MSstatsTMT (R package ).
Myocardial LV tissue samples were homogenized in RIPA buffer (Bio-Rad Laboratories, Hercules, CA, USA) containing protease and phosphatase inhibitor cocktail (Roche, Basel, Switzerland), using the Bertin Precellys 24 Tissue Homogenizer with the Bertin Cryolys cooling system (Bertin Technologies). The concentrations of the extracted proteins were measured by BCA assay (Thermo Fisher Scientific). Then, protein homogenates were suspended in sample buffer and heated at 70 °C for 10 min. A total of 20 µg protein for each sample was loaded onto 6–12% acrylamide gels and separated with a sodium dodecyl sulfate polyacrylamide gel electrophoresis system (Bio-Rad Laboratories). Gels were transferred to polyvinylidene fluoride membranes under dry conditions. Membranes were then washed and blocked for 1 h in 5% bovine serum albumin (BSA) in Tris-buffered saline Tween 20 (TBST) at room temperature. Next, membranes were incubated overnight at 4 °C with the following primary antibodies diluted in 2.5% BSA in TBST. The following primary antibodies were used: from Cell Signaling (Cell Signaling Technology, Danvers, MA, USA): VASP (1:1000, ID:#3112), ATP2A2 (1:1000, ID:#4388) and the housekeeping GAPDH (1:5000; ID: #5174); from R&D Systems (R&D Systems Inc., Minneapolis, MN, US): POSTN (1:1000, ID: AF3548); from Abcam (Abcam Inc., Toronto, ON, Canada): OPN (1:1000, ID: ab63856). The blots were washed and incubated with horseradish peroxidase-conjugated secondary antibody (1:5000, 2.5% BSA in TBST) for 2 h at room temperature. The immunoreactive protein bands were developed using Super Signal West Pico Plus (Thermo Fisher Scientific) or SuperSignal West Femto chemiluminescent substrate. The intensity of the immunoblot bands was analyzed with Image J (NIH, Bethesda, MD, USA). The intensity of the bands of the primary targets was normalized to that of the housekeeping GAPDH on the same blot.
The concentrations of circulating endogenous steroid hormones (ALDO aldosterone, PROG progesterone, 2-OHE1 2-hydroxyestrone, 4-OHE1 4-hydroxyestrone, 4-MeE2 4-methoxyestradiol, DHT dihydrotestosterone, 2-MeE1 2-methoxyestrone, 2-MeE2 methoxyestradiol, 16-OHE1 + 16OE2 16-hydroxyestrone + 16-ketoestradiol, E2 estradiol, CORT corticosterone, DOC 11-deoxycorticosterone, E1 estrone) were quantitated using in-house methods developed and validated at the Department of Laboratory Medicine, Semmelweis University, relying on ultra-high performance liquid chromatography-triple quadrupole mass spectrometry with positive electrospray ionization and multiple reaction monitoring (Shimadzu Nexera X2-LCMS-8060, Simkon Kft., Budapest, Hungary). Isotopically labeled internal standards were added to the samples at the beginning of sample preparation to correct for random errors during preparation and analysis. Following deproteinization of 200 µL serum with 600 µL methanol containing the internal standards, ALDO, DHT, DOC, CORT and PROG were extracted by solid phase extraction (Phenomenex Strata-X 60 mg, Gen-Lab Kft, Budapest, Hungary). E1, E2, 2-OHE1, 4-OHE1, 2-MeE1, 2-MeE2, 4-MeE2, 16OHE1 and 16OE2 were first released from their conjugates by incubating another 500-µL fraction of serum with 500 µL β-glucuronidase/aryl sulfatase in pH = 5.0 acetate buffer at 60 °C for 120 min. Subsequently, the mixture was cooled and mixed intensively with 2 × 1 mL ethyl acetate. Following the evaporation of the combined fractions of the organic solvent, 100 µL 1 mg/mL dansyl chloride in acetonitrile and 20 µL 0.5 mol/L sodium carbonate in water were added to the residue. The mixture was incubated at 45 °C for 15 min, followed by the addition of 20 µL 0.5 mol/L hydrochloric acid.
Statistics and bioinformatics
The data are presented as mean ± SEM. Normal distribution was confirmed by the Shapiro-Wilks method. Two-way analysis of variance (ANOVA) on the factors of “sex” and “ischemia”, followed by a Tukey HSD post hoc test to examine intergroup differences. Contingency tables were analyzed with Fisher’s exact test. Correlations were tested using Pearson’s method. Two-sample T-test was used to compare the area of intact myofibers in the F-isch subgroups. Relationships of functional, proteomic and hormonal profiles were established using log- transformed and quantile normalized values to account for differences among datasets obtained with different methodologies in a subgroup analysis of F-Isch animals. A P-value of < 0.05 was deemed significant for all but gene ontology enrichment analyses (executed with clusterProfiler ), where an adjusted P-value of < 0.01 was defined as the threshold.
Differential expression analysis was performed with Linear Models for Microarray Data (Limma R package ) by setting up a 2 × 2 factorial design. Significantly changed proteins were grouped into functional groups based on an extensive literature search. Supervised algorithms from the MixOmics R package  were utilized to establish associations among functional, proteomic and hormonal datasets. Partial Least Squares Discriminant Analysis (PLS-DA) was employed to examine hormonal-functional relationships where a Variable Importance in the Projection (VIP) score of more than one was considered influential in feature selection. sparse PLS-DA (sPLS-DA) was utilized to provide a visual overview of the extent of inter-sex and ischemic-control differences. Projection to Latent Structure (PLS) analysis was performed to screen for hormone-protein associations, where a VIP score of more than 1.3 was considered influential in feature selection. K-fold cross-validation from MixOmics was used to test the classification efficiency of the models. As a measure of model performance, balanced error rate (BER) was calculated. The number of components to be included in the models was based on the lowest BER and lowest SEM of BER that was achievable.
Assessment of effect of isoproterenol treatment
Subcutaneous injection of ISO was followed by an immediate and comparable rise in HR (Fig. 1B) in the Isch groups. The induction of ischemia on two consecutive days has resulted in 38% and 17% mortality in M-Isch and F-Isch rats respectively. The HW/TL ratio was greatly increased, while the area of intact myofibers as quantified on the H&E stained histological slides was markedly decreased in both M-Isch and F-Isch animals (Fig. 1C, D, F). Myocardial collagen content was assessed on picrosirius stained slides and was not found to be different in any of our experimental groups (Fig. 1E, G). These data suggest similar LV myocardial damage in male and female animals that survived the ischemic insult.
LV functional characterization
Although 48 h after the first injection, LVESP was significantly lowered in both Isch groups, MAP, SABP, SADP and SW were found to be markedly diminished only in M-Isch. The decrements of both dP/dtmax and dP/dtmax-EDV were not found to be significant in either sex after ischemia, whereas EF and the gold standard index of contractility (PRSW) have indicated a marked deterioration in systolic function in M-Isch animals. Furthermore, cardiac efficiency was also more severely reduced in males after ischemia (Fig. 2 and Table 1).
Rats in the F-Isch group were burdened by markedly decreased CO. In line with this finding, F-Isch rats experienced a higher, but not significant decrease in SV compared to sex-matched control groups. Aspects of diastolic function: active relaxation as determined by Tau and myocardial stiffness as evidenced by LVEDP and EDPVR have deteriorated to a greater extent in F-Isch (Fig. 2 and Table 1).
2224 proteins were identified and quantified from the myocardial samples. sPLS-DA performed on the entire quantified proteome of the myocardium shows that the four experimental groups can be clearly separated from each other with the inclusion of two components (0.0002 ± 0.0069 BER), thus highlighting that there are consistent protein level differences between our groups (Fig. 3B). According to sPLS-DA analysis, the most important proteins in the successful separation of Isch vs. Co were JPH2, ATP5ME, KNG2, CAVIN2, APOOL, MRRF and COQ8A, while the greatest distinction between males and females could be made based on the abundances of A1BG, LIFR, FKBP7, RCN3.
Fisher’s exact test showed a statistically significant association between the number of sex-specific differentially regulated proteins in response to ischemia and sexes (p < 0.001). Although the number of downregulated proteins was similar (603 vs 607, F-Isch vs. M-Isch), females had 56% more sex-specific upregulated proteins compared to males (455 vs. 291, F-Isch vs. M-Isch) (Fig. 3A, Additional file 2: Table S1). The main biological categories that were found to be differentially regulated were comprised of proteins with a role in transcription, cytoskeleton & contraction, extracellular matrix & wound healing and inflammation. In all these categories, females exhibited a higher number of significantly upregulated proteins. Although under the threshold of significance, more metabolic and mitochondrial proteins were downregulated in ischemic males than in females after ischemia (Fig. 3C). To validate the proteomic measurements in the cases of sex-dependent and sex-independent upregulation, as well as unchanged protein expression after ischemia, four proteins with broadly acknowledged cardiovascular relevance were selected: VASP, POSTN, OPN, ATP2A2. Western blotting has confirmed the female-specific upregulation of VASP and OPN in F-Isch, and the sex-independent upregulation of POSTN in M-Isch, as well as in F-Isch compared to respective controls. Furthermore, the expression of ATP2A2 was consistently unaffected by ischemia according to both western blotting and LC–MS/MS measurements (Additional file 1: Fig. S1).
Association of steroid hormones and cardiac functional outcomes in females after ischemia
A relationship between functional outcomes and circulating steroid hormones was investigated in the rats of the F-Isch group in a detailed subgroup analysis. Hierarchical clustering of all P–V hemodynamic parameters has automatically identified two distinct subgroups of 7–7 F-Isch animals which were labeled as mild and severe functional outcomes (Fig. 4A). We investigated if the functional differences observed between these clusters may be related to different percentages of areas of intact myofibers left after infarction (an indirect measurement of infarct size) in the mild and severe F-Isch groups. Nevertheless, no significant difference was found (79.774 ± 2.963% vs 78.259 ± 2.953%, F-Isch mild vs. F-Isch severe). PLS-DA analysis was performed on the steroid hormone profiles of F-Isch rats in an attempt to predict functional outcomes based on hormonal profiles (0.307 ± 0.063 BER, Fig. 4B). Aldosterone, progesterone, 2-hydroxyestrone, 4-hydroxyestrone and 4-methoxyestradiol were found to be influential in the classification of mild vs. severe outcomes based on VIP scores extracted from the model (Fig. 4C, D). Strong correlations were found between 2-hydroxyestrone and crucial P–V parameters of systolic and diastolic function (PRSW, Tau and dP/dtmin, Fig. 4E–G). The concentration of 2-hydroxyestrone correlated with the HW/TL as well in the F-Isch group (Fig. 4H).
Association of hormonal and proteomic profiles in females after ischemia
PLS was used to establish associations among circulating steroid hormone concentrations and the post-ischemic female myocardial proteome. When VIP scores of proteins per hormone were extracted from the model, a protein with a VIP score of more than 1.3 for a given hormone was considered to show an association with the concentration of the hormone. The top four hormones that were associated with the most proteins are as follows: 4-hydroxyestrone (n = 381), 2-hydroxyestrone (n = 196), estradiol (n = 141) and corticosterone (n = 93). As a second step gene ontology biological process (GO:BP) enrichment analysis was performed on all of the hormone-associated protein groups to find common biological processes, the regulation thereof individual hormones might be responsible for. Only three protein groups have shown significant GO:BP enrichment: 4-hydroxyestrone, 2-hydroxyestrone and estradiol (Fig. 5). Estradiol and 2-hydroxyestrone were associated with proteins responsible for oxidative phosphorylation and aerobic respiration. Furthermore, estradiol was found to associate with proteins from the tricarboxylic acid cycle, while fatty acid beta-oxidation was associated only with 2-hydroxyestrone. Proteins linked to 4-hydroxyestrone were also implicated in cellular respiration, however, elements of supramolecular fibers and the cytoskeleton, in particular proteins related to actin filaments have shown a strong association.
To the best of our knowledge, this is the first study to employ the gold standard methods of pressure–volume analysis, LC–MS/MS based proteomics and steroid hormone measurements to evaluate sex-specific early LV functional alterations following ischemia in a comprehensive manner. Since highly sensitive indices of cardiac performance were measured in vivo, our results reflect LV function independent of loading conditions. We have also provided thorough profiles of sex-dependent proteomic changes by achieving deep proteome coverage on a high number of bioreplicates compared to an earlier cardiac proteomics study performed only on male rats by Lim et al. .
Assessment of sex-independent effects of ischemia
The isoproterenol-induced ischemia model has some unique characteristics compared to surgical alternatives. Instead of permanent or transient ligature of coronary arteries, where the ischemic insult is delimited by the borders of the area that is perfused by the given coronary artery, high dose isoproterenol induces global ischemia that affects the entire myocardium with special emphasis on subendocardial regions. Consequently, this model can be best translated to cases of global ischemia. Subcutaneous administration of isoproterenol led to a sudden and comparable rise in HR in both sexes (Fig. 1). Upon histologic evaluation, small necrotic islands were spread throughout the myocardium, thereby confirming the global nature of ischemia. Using digital whole slide analysis and H&E color deconvolution, the percentage of the area of intact myofibers was calculated, which is an indirect measure of the extent of necrosis. Although previous studies pointed to lower mortality and smaller infarct sizes in females after ischemia (as reviewed by Regitz-Zagrosek and Kararigas ) our investigation did not confirm the latter. Our findings are more in line with the works of Shioura et al. who also found that infarct sizes were not different in male and female mice , thus sex-specific functional outcomes are not a direct consequence of differences in the extent of myocardial damage. HW/TL was increased in both sexes as well, which was most likely due to myocardial edema caused by the initiation of inflammatory processes, as part of the post-ischemic myocardial remodeling. Correspondingly, among the 7 proteins that have shown the most consistent, sex-independent changes in abundance when separating Isch and Co animals as identified by sPLS-DA analysis, we found KNG2 (T-kininogen 2), an important contributor to inflammatory processes. Of the 7, further 4 proteins (Atypical kinase COQ8A, Ribosome-recycling factor MRRF, MICOS complex subunit APOOL, ATP synthase subunit e ATP5ME) have indicated mitochondrial impairment and another 2 membrane-associated proteins (Junctophilin-2 JPH2, Caveolae-associated protein 2 CAVIN2) have pointed to the disintegration of the T-tubular network of cardiomyocytes , that is essential for effective excitation–contraction coupling.
Systolic function and its proteomic background
Female hearts were found to exhibit improved recovery of contractility in experiments performed on Langendorff perfused hearts after isoproterenol treatment [29, 30]. In line with this finding, 48 h after the first isoproterenol injection, a significant decline in LV contractility (PRSW, EF) was observed only in male ischemic rats, while females displayed overall preserved systolic function. Likewise, a marked decrease in arterial pressure (SABP, DABP, MAP) could only be detected in male animals. Although on the proteomic level we were able to observe a comparable number of downregulated proteins in response to ischemia between the sexes, 56% more proteins were upregulated in females compared to males. After grouping sex-specific up- and downregulated proteins according to biological function based on extensive literature and database searches, we have found significantly more cytoskeletal and contractile proteins upregulated in the myocardium of F-Isch than in M-Isch rats. Together with increased transcriptional activity in females, this could indicate faster replenishment of damaged contractile structures and membrane scaffolding proteins and thus could lead to more effective excitation–contraction coupling and better functional recovery. Although not statistically significant, we have found a greater number of DNA repair enzymes and molecular chaperones upregulated exclusively in females. Of note are ORP-150 , HSP90AA1 , HSP90AB1 , CALR , PDIA3  and BIP  which are known to exert protective effects under ischemic conditions. As a further source of intersex disparity, we have found a tendency towards a higher number of mitochondrial proteins which were downregulated only in males. This finding could be associated with impaired myocardial energetics that corresponds to the observed decrease in cardiac efficiency of the M-Isch, group, while females displayed maintained Eff values. Of particular interest were proteins with a pivotal role in the citric acid cycle (ACO2, OGDH, IDH3A, IDH2, FH) and fatty acid oxidation (ECI1, ETFRF1, ETFA, ETFB, ETFDH, CISD1, CPT2, ACOT2), all of which were selectively downregulated in males following ischemia but have not changed significantly in females. When we took a look at glycolytic processes, differing strategies of adaptation became apparent. In males, ischemia led to a downregulation of β-enolase, while in females an upregulation of α-enolase was found. α-enolase has been linked to improved cardiac contractility after ischemia , which is in line with our current findings.
Diastolic function and its proteomic background
Conflicting studies have been published regarding diastolic function. In Langendorff setups, the female heart seems to show greater immediate recovery after ischemia [29, 30]. However, in vivo measurements 10 weeks after ischemia indicated increased myocardial stiffness in females , possibly due to enhanced extracellular remodeling . In our in vivo study of short-term cardiac functional adaptation, we have found significantly increased end-diastolic pressure in the left ventricle of F-Isch rats. This could be explained by a simultaneous rise in myocardial stiffness, as suggested by EDPVR. Furthermore, Tau has indicated a significant deterioration of active relaxation in females, but not in males, following ischemia. Diastolic impairments can result in a filling deficit and thus might limit the volume of blood that can be ejected during systole. Analogously, females were characterized by a numerical decrease in SV and a significant decrease in CO after ischemia, which was not observed in males. On the proteomic level, we have found a significantly higher quantity of upregulated proteins linked to extracellular matrix remodeling, wound healing and inflammation in females after ischemia. Since extracellular remodeling is a long-term process, inflammation and interstitial fluid accumulation might predominantly determine myocardial stiffness in an acute scenario . Indeed, we have found marked increases in HW/TL after isoproterenol treatment in our experiment, which could be attributed to myocardial edema. Thus, upregulation of inflammatory proteins may contribute to the observed short-term deterioration of diastolic function in females. Nevertheless, this may be beneficial overall, as inhibition of inflammatory processes in the first days after MI was found to result in worse long-term outcomes by impeding cardiac remodeling and tissue repair . Moreover, immediate alterations of extracellular proteins in this critical initial phase of myocardial remodeling should also be addressed. Even though we did not find a substantial increase of collagen accumulation on our picrosirius red stained slides, our proteomic results indicated that extracellular remodeling has already begun. Proteins of interest for the structural integrity of the myocardium are COL18A1, inhibition of which was found to deteriorate myocardial remodeling after infarction , and ITGA5, deletion of which was found to inhibit tissue fibrosis in multiple organs , both of which were uniquely upregulated in females. Of particular note were matricellular proteins: we have identified SPARC, SPP1 and CTSB to be sex-specifically upregulated in females, while CST3 was upregulated only in males. A transient early rise in SPARC expression has been linked to both fibroblast migration and effective early extracellular remodeling after ischemia [42, 43], while SPP1 was found to be essential for effective collagen synthesis preventing LV dilation. Furthermore, CTSB deficiency was found to slow myocardial remodeling . Additionally, CST3 was identified as a potent inhibitor of CTSB , which might contribute to the more passive remodeling strategy of males. This is in line with earlier findings of delayed myocardial remodeling in males compared to females after ischemia .
Associations of sex hormone levels with functional outcomes and proteomic alterations in the female ischemic group.
It is established that estrogen hormone signaling greatly influences the outcome of myocardial ischemia . The most studied estrogen in this respect is E2, but it is far from the only sex hormone that could exert potentially beneficial or detrimental effects. The majority of metabolites of both E1 and E2 are known to bind to estrogen receptor subtypes  and may contribute to sex-related differences. To date, gonadectomy and estrogen receptor deletion have been the most commonly used targeted strategies for studying the hormonal influence on post-ischemic recovery [4, 5]. However, both approaches lead to a complete blockade of estrogenic effects, independently of the specific hormones or hormone metabolites. Complete deprivation of estrogen hormones through gonadectomy led to worse outcomes after ischemia, while E2 administration was shown to ameliorate post-ischemic functional impairments and enhance myocardial remodeling in most studies . However, in an investigation by Cavasin et al. E2 supplementation after gonadectomy failed to live up to its protective potential, while testosterone was found to be detrimental to cardiac remodeling and functional recovery . To study the potential influence of sex hormones and their metabolites on post-ischemic recovery, concentrations of a wide range of steroid hormones and metabolites were measured in the F-Isch group. Based on all parameters of left-ventricular function obtained with P–V analysis, automatic hierarchical clustering was performed on rats of the F-Isch group, which identified two distinct clusters of F-Isch rats: one with a mild deterioration of cardiac function and one with a severe impairment of systolic and diastolic performance. To examine the potential role of circulating steroid hormones on functional outcomes, PLS-DA supervised algorithm was utilized for feature selection. Of all the measured hormones, the contributions of PROG, ALDO, 2-OHE1, 4-OHE1 and 4-MeE2 were the most important in the successful classification of mild vs. severe functional outcomes, thus their possible involvement in regulating the pivotal processes of post-ischemic recovery was considered. Out of the five influential hormones, only 2-OHE1 has shown significant and strong correlations with our functional and structural parameters (PRSW, dP/dtmin, Tau and HW/TL). As a follow-up, protein groupings were established using PLS modeling based on the strength of the associations they showed with the concentration of circulating steroid hormones. GO:BP analysis performed on all hormone-related groups of proteins found significant enrichment only in protein groups associated with 4-OHE1, 2-OHE1 and E2. E2 and 2-OHE1 have shown consistent associations with proteins of oxidative phosphorylation and aerobic respiration. As the main sources of disparity in mitochondrial differential regulation after ischemia, the citric acid cycle and fatty acid oxidation were identified. The former showed a strong association with E2, while the latter seemed to be more related to 2-OHE1. Cytoskeletal and contractile proteins were linked principally to the concentration of 4-OHE1 in particular supramolecular fibers and proteins responsible for the build-up of actin filaments.
In our rat model of isoproterenol-induced myocardial ischemia, P–V analysis uncovered sex-specific patterns of LV dysfunction: males were characterized by predominantly systolic, while females principally diastolic dysfunction. Distinct proteomic profiles underlay these functional differences. The female post-ischemic myocardial adaptation involved an overall higher number of upregulated proteins, notably in the domains of inflammation, extracellular remodeling, and molecular chaperones. Males exhibited sex-specific downregulation of several proteins involved in the citric acid cycle and fatty acid oxidation. Detailed subgroup analysis of the ischemic female group highlighted the potential role of E2 and estrogen metabolites in determining post-ischemic functional outcomes and proteomic profiles.
The interpretation of the results from the current study is limited to young rats. All data have been acquired from male and female rats that have survived the two-day induction of intermittent myocardial ischemia. The pharmacologic induction of acute global ischemia utilized in our investigation may not be directly translatable to chronic ischemic heart disease in humans. Blood samples for hormonal profiling were taken after P–V measurements to ensure undisturbed hemodynamic characterization. Steroid hormones can be produced and metabolized locally inside the tissue as well, thus circulating levels of sex hormones might not accurately represent hormone levels inside the myocardium. Further studies are needed to assess the effect of supplementation of the steroid hormones highlighted in this study on gonadectomized ischemic rodent models.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The proteomic raw data files supporting the conclusions of this article are available in the Mass Spectrometry Interactive Visual Repository (MassIVE, MSV000088184 dataset).
Regitz-Zagrosek V, Kararigas G. Mechanistic pathways of sex differences in cardiovascular disease. Physiol Rev. 2017;97:1–37.
Zhai P, Eurell TE, Cotthaus R, Jeffery EH, Bahr JM, Gross DR. Effect of estrogen on global myocardial ischemia-reperfusion injury in female rats. Am J Physiol Heart Circ Physiol. 2000;279:H2766-2775.
Rattanasopa C, Phungphong S, Wattanapermpool J, Bupha-Intr T. Significant role of estrogen in maintaining cardiac mitochondrial functions. J Steroid Biochem Mol Biol. 2015;147:1–9.
Booth EA, Marchesi M, Kilbourne EJ, Lucchesi BR. 17Beta-estradiol as a receptor-mediated cardioprotective agent. J Pharmacol Exp Ther. 2003;307:395–401.
Booth EA, Marchesi M, Knittel AK, Kilbourne EJ, Lucchesi BR. The pathway-selective estrogen receptor ligand WAY-169916 reduces infarct size after myocardial ischemia and reperfusion by an estrogen receptor dependent mechanism. J Cardiovasc Pharmacol. 2007;49:401–7.
Jain M, Liao R, Podesser BK, Ngoy S, Apstein CS, Eberli FR. Influence of gender on the response to hemodynamic overload after myocardial infarction. Am J Physiol Heart Circ Physiol. 2002;283:H2544-2550.
Bae S, Zhang L. Gender differences in cardioprotection against ischemia/reperfusion injury in adult rat hearts: focus on Akt and protein kinase C signaling. J Pharmacol Exp Ther. 2005;315:1125–35.
Brown DA, Lynch JM, Armstrong CJ, Caruso NM, Ehlers LB, Johnson MS, Moore RL. Susceptibility of the heart to ischaemia-reperfusion injury and exercise-induced cardioprotection are sex-dependent in the rat. J Physiol. 2005;564:619–30.
Chen A, Chen Z, Xia Y, Lu D, Jia J, Hu K, Sun A, Zou Y, Qian J, Ge J. Proteomics analysis of myocardial tissues in a mouse model of coronary microembolization. Front Physiol. 2018;9:1318.
Pacher P, Nagayama T, Mukhopadhyay P, Batkai S, Kass DA. Measurement of cardiac function using pressure-volume conductance catheter technique in mice and rats. Nat Protoc. 2008;3:1422–34.
Koricanac G, Tepavcevic S, Romic S, Zivkovic M, Stojiljkovic M, Milosavljevic T, Stankovic A, Petkovic M, Kamceva T, Zakula Z. Estradiol enhances effects of fructose rich diet on cardiac fatty acid transporter CD36 and triglycerides accumulation. Eur J Pharmacol. 2012;694:127–34.
Fang L, Gao XM, Moore XL, Kiriazis H, Su Y, Ming Z, Lim YL, Dart AM, Du XJ. Differences in inflammation, MMP activation and collagen damage account for gender difference in murine cardiac rupture following myocardial infarction. J Mol Cell Cardiol. 2007;43:535–44.
Biondi-Zoccai GG, Abate A, Bussani R, Camilot D, Giorgio FD, Marino MP, Silvestri F, Baldi F, Biasucci LM, Baldi A. Reduced post-infarction myocardial apoptosis in women: a clue to their different clinical course? Heart. 2005;91:99–101.
Lim SH, Lee J, Han MJ. Comprehensive analysis of the cardiac proteome in a rat model of myocardial ischemia-reperfusion using a TMT-based quantitative proteomic strategy. Proteome Sci. 2020;18:2.
Zhang Z, Wu S, Stenoien DL, Pasa-Tolic L. High-throughput proteomics. Annu Rev Anal Chem. 2014;7:427–54.
Grant MKO, Abdelgawad IY, Lewis CA, Seelig D, Zordoky BN. Lack of sexual dimorphism in a mouse model of isoproterenol-induced cardiac dysfunction. PLoS ONE. 2020;15:e0232507.
Korkmaz S, Radovits T, Barnucz E, Hirschberg K, Neugebauer P, Loganathan S, Veres G, Pali S, Seidel B, Zollner S, et al. Pharmacological activation of soluble guanylate cyclase protects the heart against ischemic injury. Circulation. 2009;120:677–86.
Ruppert M, Korkmaz-Icoz S, Loganathan S, Jiang W, Lehmann L, Olah A, Sayour AA, Barta BA, Merkely B, Karck M, et al. Pressure-volume analysis reveals characteristic sex-related differences in cardiac function in a rat model of aortic banding-induced myocardial hypertrophy. Am J Physiol Heart Circ Physiol. 2018;315:H502–11.
Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, McQuaid S, Gray RT, Murray LJ, Coleman HG, et al. QuPath: open source software for digital pathology image analysis. Sci Rep. 2017;7:16878.
van Putten M, de Winter C, van Roon-Mom W, van Ommen GJ. t Hoen PA, Aartsma-Rus A: a 3 months mild functional test regime does not affect disease parameters in young mdx mice. Neuromuscul Disord. 2010;20:273–80.
Baumert HM, Metzger E, Fahrner M, George J, Thomas RK, Schilling O, Schule R. Depletion of histone methyltransferase KMT9 inhibits lung cancer cell proliferation by inducing non-apoptotic cell death. Cancer Cell Int. 2020;20:52.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26:1367–72.
Huang T, Choi M, Tzouros M, Golling S, Pandya NJ, Banfai B, Dunkley T, Vitek O. MSstatsTMT: statistical detection of differentially abundant proteins in experiments with isobaric labeling and multiple mixtures. Mol Cell Proteomics. 2020;19:1706–23.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Phipson B, Lee S, Majewski IJ, Alexander WS, Smyth GK. Robust hyperparameter estimation protects against hypervariable genes and improves power to detect differential expression. Ann Appl Stat. 2016;10:946–63.
Rohart F, Gautier B, Singh A, Le Cao KA. mixOmics: an R package for ’omics feature selection and multiple data integration. PLoS Comput Biol. 2017;13:e1005752.
Shioura KM, Geenen DL, Goldspink PH. Sex-related changes in cardiac function following myocardial infarction in mice. Am J Physiol Regul Integr Comp Physiol. 2008;295:R528-534.
Chan BYH, Roczkowsky A, Cho WJ, Poirier M, Lee TYT, Mahmud Z, Schulz R. Junctophilin-2 is a target of matrix metalloproteinase-2 in myocardial ischemia-reperfusion injury. Basic Res Cardiol. 2019;114:42.
Gabel SA, Walker VR, London RE, Steenbergen C, Korach KS, Murphy E. Estrogen receptor beta mediates gender differences in ischemia/reperfusion injury. J Mol Cell Cardiol. 2005;38:289–97.
Wang M, Tsai BM, Reiger KM, Brown JW, Meldrum DR. 17-beta-Estradiol decreases p38 MAPK-mediated myocardial inflammation and dysfunction following acute ischemia. J Mol Cell Cardiol. 2006;40:205–12.
Aleshin AN, Sawa Y, Kitagawa-Sakakida S, Bando Y, Ono M, Memon IA, Tohyama M, Ogawa S, Matsuda H. 150-kDa oxygen-regulated protein attenuates myocardial ischemia-reperfusion injury in rat heart. J Mol Cell Cardiol. 2005;38:517–25.
Zhu WS, Guo W, Zhu JN, Tang CM, Fu YH, Lin QX, Tan N, Shan ZX. Hsp90aa1: a novel target gene of miR-1 in cardiac ischemia/reperfusion injury. Sci Rep. 2016;6:24498.
Haase M, Fitze G. HSP90AB1: helping the good and the bad. Gene. 2016;575:171–86.
Wang JL, Li YZ, Tao TQ, Wang XR, Wang Y, Song DD, Liu XH. Postconditioning with calreticulin attenuates myocardial ischemia/reperfusion injury and improves autophagic flux. Shock. 2020;53:363–72.
Yoo DY, Cho SB, Jung HY, Kim W, Lee KY, Kim JW, Moon SM, Won MH, Choi JH, Yoon YS, et al. Protein disulfide-isomerase A3 significantly reduces ischemia-induced damage by reducing oxidative and endoplasmic reticulum stress. Neurochem Int. 2019;122:19–30.
Bodalia A, Li H, Jackson MF. Loss of endoplasmic reticulum Ca2+ homeostasis: contribution to neuronal cell death during cerebral ischemia. Acta Pharmacol Sin. 2013;34:49–59.
Diaz-Ramos A, Roig-Borrellas A, Garcia-Melero A, Lopez-Alemany R. alpha-Enolase, a multifunctional protein: its role on pathophysiological situations. J Biomed Biotechnol. 2012;2012:156795.
Jugdutt BI. Ventricular remodeling after infarction and the extracellular collagen matrix: when is enough enough? Circulation. 2003;108:1395–403.
Maleki Dizaji N, Garjani A, Mousavi S, Mohammadi M, Vaez H. Time-dependent influence of infliximab on hemodynamic responses and cardiac injuries of isoproterenol-induced myocardial infarction in rats. Eur J Pharmacol. 2021;903:174122.
Isobe K, Kuba K, Maejima Y, Suzuki J, Kubota S, Isobe M. Inhibition of endostatin/collagen XVIII deteriorates left ventricular remodeling and heart failure in rat myocardial infarction model. Circ J. 2010;74:109–19.
Henderson NC, Arnold TD, Katamura Y, Giacomini MM, Rodriguez JD, McCarty JH, Pellicoro A, Raschperger E, Betsholtz C, Ruminski PG, et al. Targeting of alphav integrin identifies a core molecular pathway that regulates fibrosis in several organs. Nat Med. 2013;19:1617–24.
McCurdy SM, Dai Q, Zhang J, Zamilpa R, Ramirez TA, Dayah T, Nguyen N, Jin YF, Bradshaw AD, Lindsey ML. SPARC mediates early extracellular matrix remodeling following myocardial infarction. Am J Physiol Heart Circ Physiol. 2011;301:H497-505.
Wu RX, Laser M, Han H, Varadarajulu J, Schuh K, Hallhuber M, Hu K, Ertl G, Hauck CR, Ritter O. Fibroblast migration after myocardial infarction is regulated by transient SPARC expression. J Mol Med. 2006;84:241–52.
Wu QQ, Xu M, Yuan Y, Li FF, Yang Z, Liu Y, Zhou MQ, Bian ZY, Deng W, Gao L, et al. Cathepsin B deficiency attenuates cardiac remodeling in response to pressure overload via TNF-alpha/ASK1/JNK pathway. Am J Physiol Heart Circ Physiol. 2015;308:H1143-1154.
Xie L, Terrand J, Xu B, Tsaprailis G, Boyer J, Chen QM. Cystatin C increases in cardiac injury: a role in extracellular matrix protein modulation. Cardiovasc Res. 2010;87:628–35.
Cavasin MA, Tao Z, Menon S, Yang XP. Gender differences in cardiac function during early remodeling after acute myocardial infarction in mice. Life Sci. 2004;75:2181–92.
Nikolic I, Liu D, Bell JA, Collins J, Steenbergen C, Murphy E. Treatment with an estrogen receptor-beta-selective agonist is cardioprotective. J Mol Cell Cardiol. 2007;42:769–80.
Zhu BT, Han GZ, Shim JY, Wen Y, Jiang XR. Quantitative structure-activity relationship of various endogenous estrogen metabolites for human estrogen receptor alpha and beta subtypes: Insights into the structural determinants favoring a differential subtype binding. Endocrinology. 2006;147:4132–50.
Patten RD, Pourati I, Aronovitz MJ, Baur J, Celestin F, Chen X, Michael A, Haq S, Nuedling S, Grohe C, et al. 17beta-estradiol reduces cardiomyocyte apoptosis in vivo and in vitro via activation of phospho-inositide-3 kinase/Akt signaling. Circ Res. 2004;95:692–9.
Cavasin MA, Tao ZY, Yu AL, Yang XP. Testosterone enhances early cardiac remodeling after myocardial infarction, causing rupture and degrading cardiac function. Am J Physiol Heart Circ Physiol. 2006;290:H2043-2050.
Expert technical assistance of Henriett Biró, Gábor Fritz, Urbán Edina, Prokaj Benjamin, Luisa Braun and Bettina Wehrle is greatly acknowledged. We thank Daniella Vörösmarti for her contribution to the graphical abstract.
Open Access funding provided by Semmelweis University. This work was supported by the New National Excellence Program of the Ministry for Innovation and Technology from the source of the National Research and Development Fund [ÚNKP-20-3-I-SE-1 to B.A.B]; the National Heart Program of the National Research, Development and Innovation Fund of Hungary, financed under the NVKP_16 funding scheme [NVKP_16-1–2016-0017]; the Thematic Excellence Programme of the Ministry for Innovation and Technology in Hungary within the framework of the Therapeutic Development and Bioimaging thematic programmes of Semmelweis University [2020-4.1.1.-TKP2020]; the National Research, Development and Innovation Office (NKFIH) of Hungary [K134939 to T.R.]; the Deutsche Forschungsgemeinschaft [SCHI 871/11-1, SCHI 871/15-1, GR 4553/5-1, PA 2807/3-1, INST 39/1244-1 (P12), INST 39/766-3 (Z1), GRK 2606 "ProtPath"]; the ERA PerMed programme [BMBF, 01KU1916, 01KU1915A], the German-Israel Foundation [Grant No. 1444], and the German Consortium for Translational Cancer Research [project Impro-Rec].
Ethics approval and consent to participate
Our investigation was performed as approved by the Ethics Committee of Hungary for Animal Experimentation, Semmelweis University, Budapest (PEI/001/2374-4/2015) conforming to European Directive 2010/63/EU and the Guide for the Care and Use of Laboratory Animals used by the US National Institutes of Health (Publication No. 85-23, revised 1996).
Consent for publication
This study consists of animal data and is devoid of any human data, thus consent for publication is not applicable.
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.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Validation of LC–MS/MS measurements with western blot. A–D: Relative protein expression of VASP, POSTN, OPN and ATP2A2 as measured by LC–MS/MS and western blot. Values were normalized to the mean of the corresponding control group. Statistical significance of post hoc test compared to same-sex control is highlighted as follows: *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. M-Co male control, M-Isch male ischemic, F-Co female control, F-Isch female ischemic.
Results of differential expression analysis performed with Limma.
About this article
Cite this article
Barta, B.A., Ruppert, M., Fröhlich, K.E. et al. Sex-related differences of early cardiac functional and proteomic alterations in a rat model of myocardial ischemia. J Transl Med 19, 507 (2021). https://doi.org/10.1186/s12967-021-03164-y