2.1. Cell Cycle Regulation of Cardiomyocytes in during the Early Neonatal Period
We began our analysis with the assessment of heart weight/body weight (HBW) ratios in neonatal CM in the first 7d after birth. We found that HBW did not significantly change over this period (
Figure 1A). However, with the observation that mRNAs involved in DNA replication and mitosis showed a biphasic periodicity, we investigated DNA synthesis in CM in the first 7d after birth. Neonatal C57BL/6J wild-type (wt) mice were injected subcutaneously with a single dose of the 5-ethynyl-2′-deoxyuridine (EDU), a thymidine analogue[
37]. Cardiac specimens were prepared for indirect immunofluorescence microscopy employing an EDU assay, and anti-phospho Histone Pi.H3-Ser28 (H3.Pi-S28) antibodies (
Figure 1B)[
38] to identify CM in S- and M-phase, respectively, in conjunction with anti-cardiac actinin and 4′,6-diamidino-2-phenylindole (Dapi). We also co-stained neonatal cardiac sections with an antibody to pericentriolar material 1 (Pcm1), a CM-specific peri-nuclear marker[
39,
40], and EDU to specifically identify CM in S-phase (
Figure 1C). Using this approach, the CM labeling index was monitored from embryonic day 21 (0d), neonatal 12h and 1d through 7d (
Figure 1C). Exceedingly high labeling indexes were detected at 0d and 12h (
Figure 1D). In contrast, levels of EDU/Pcm-1 double-positive CM nuclei dropped dramatically at 1d (
Figure 1D,E). This proliferative block was corroborated by immunofluorescence analysis of Pi.H3-Ser28 positive CM in M-phase (
Figure 1E), and aurora kinase B- (AurkB) positive CM in cytokinesis (
Figure 1F,G)[
41]. Re-initiation of DNA synthesis occurred at 4d, with a peak labeling index occurring at 5d (49±4.3%;
P<0.001 vs. 2d). The number of EDU-positive CM nuclei decreased again markedly by 6d (5.1±2.1%;
P<0.01 vs. 2d) (
Figure 1D). there were no sex-dependent differences noted.
To ascertain that this second wave of DNA synthesis represented CMBN, we analyzed this morphological developmental change in 2d- versus 5d-old CM by 3D reconstruction of confocal immunofluorescence micrographs of Dapi-stained cardiac sections[
42]. CMBN was essentially absent in 1d old CM (
Figure 1H). In contrast, 579.1% of CM were binucleated at the end of 5d (
P<0.001 vs. 1d). Combined, all these findings demonstrate that CM proliferation ceases in a sharply defined time window after birth and is followed by a distinct second phase of CMBN.
2.2. Differential mRNA Expression of Cell Cycle Factors in the Early Postnatal Period
To elucidate the molecular mechanisms responsible for cessation of proliferation, DNA binucleation and cell cycle exit, we performed genome-wide mRNA microarray profiling[
43]. Unsupervised hierarchical cluster analysis at a high confidence threshold revealed that approximately 1,900 individual transcripts changed significantly in murine hearts at 1d through 15d, relative to 0d (
Figure 2A). Three-dimensional principal component analysis was performed to compare gene expression profiles across the postnatal days and revealed a clustering between 1d, 3d, 7d, 10d samples, and 0d and 5d samples, consistent with our microscopic results showing cell cycle re-entry and peaks of DNA synthesis (
Figure 2B). The highest ranked ‘Gene Ontology’ (GO) terms comprised transcripts involved in the regulation of nuclear division, negative regulation of cell division and regulation of mitotic cell cycle (
Figure 2C,D,E).
Next, we validated the transcriptomic results of well-characterized cell cycle regulatory genes by qRT-qPCR analysis (
Figure 2F)[
42]. We found that the mRNA levels of a set of factors involved in the transition from S-phase into mitosis (AurkB, Skp2, AurkB), G1-Cdks (Cdk1, Cdk2, Cdk4), G1/G2-cyclins (Cyclins A, B, E, F), and key cell cycle activators (e.g. Cdc25A) and inhibitors (Cdkn2b, Cdnkn3) were significantly downregulated by postnatal day 7-10 coinciding with terminal cell cycle exit (
Figure 2F). These analyses revealed that distinct cell cycle factors are differentially regulated in the immediate postnatal period. Our transcriptomic findings were confirmed by Western blotting of CM cell cycle inhibitors, cell cycle regulators, and Mapk signaling pathways, and cardiac-specific genes (
Figure 2G)[
42]. Among all these factors, Cdk1 was amongst the top ten enriched cell cycle genes that were selectively upregulated at 1d and 5d post-birth (
Figure 2F). Thus, we decided to examine the potential role of Cdk1 in the regulation of neonatal CM proliferation
in vivo.
2.4. Loss of Cdk1 Prolongs the Proliferative Window of Cardiomyocytes during the Postnatal Stage
Initially, we investigated whether genetic silencing of Cdk1 affects DNA synthesis in
Cdk1KOc hearts in the early postnatal period. For analysis by confocal immunofluorescence microscopy, mice were injected subcutaneously with EDU, and hearts were excised, fixed, and sectioned. Subsequently, specimens were co-stained with antibodies to α-actinin, EDU, Pi.H3-Ser28, wheat germ agglutinin (WGA), and Dapi to visualize nuclear DNA. We detected a very high EDU-labeling index in CM from
control mice at 1d that dropped dramatically at 2d (
P<0.001 vs. 0d) (
Figure 4A and Figure S1,S2). This was followed by a re-induction of DNA synthesis that peaked at 5d. A progressive decrease in the number of CM in S-phase was observed over 7d and 10d, until DNA synthesis in this cell type ceased completely by 15d. Intriguingly, the numbers of EDU-positive CM in
Cdk1KOc mice remained markedly higher on 1d, 2d and 5d when compared to
controls (P<0.001). Subsequently, DNA synthesis in Cdk1-deficient CM progressively decreased over 7d-10d and was undetectable at 15d (
Figure 4A).
To analyse the distribution of CM in mitosis, cardiac sections were co-stained with antibodies to H3.Pi-S28, an M-phase-specific nuclear marker, in combination with α-actinin immunostaining to identify CM. Microscopic inspection showed dramatically elevated numbers of H3.Pi-S28-positive CM in
control animals at 0d and 5d (
Figure 4B and Figure S3). This effect was very much lower in
control CM at 7d and undetectable at 10d and 15d (P<0.01 vs 0d). Again, the number of CM in mitosis remained significantly higher in
Cdk1KOc mice at between 1d-2d and at 7d post-birth in comparison to
control animals (P<0.01). In this strain, H3.Pi-S28-positive CM were never observed at 15d (
Figure 4B). Additionally,
Cdk1KOc hearts at 4w and 6w of age failed to show any signs of cell cycle activity (Figure S4). Inspection of non-cardiomyocytes (NCM) during the first week after birth revealed no significant differences in the number of cycling NCM in
Cdk1KOc compared to
control hearts (Figure S5). Importantly, all these cell cycle events in CM are sex-independent processes, since no difference was noted in the number of H3.Pi-S28-positive CM positive CM in female mice in comparison to male mice (Figure S6).
Reportedly, CM undergo a binucleation process at 5d postnatally implying that in these cells M-phase can occur independently of cytokinesis[
8]. Therefore, we analyzed whether Cdk1 mutant CM completed cytokinesis by forming two daughter cells. Cytokinesis and polyploidization were distinguished by immunostaining with anti-AurkB antibodies to detect midbody structures between daughter CM in fixed cardiac specimen by confocal microscopy. In the presence of Cdk1, CM cytokinesis was detected only at 0d and 1d post-birth based on the analysis of AurkB positive midbody structures between CM (P<0.01 vs 2d) (
Figure 4C). The second phase of DNA synthesis at 5d was associated only with CMBN in
control hearts as indicated by the absence of AurkB-positive CM (P<0.01 vs 1d). In contrast, ablation of Cdk1 induced ongoing CM cytokinesis from 1d to 7d (P<0.001 vs 1d-7d
controls) (
Figure 4C). In particular, there was a significant increase in the percentage of binucleated CM and a decrease in mononucleated CM in
control hearts at 7d when compared to hearts from
Cdk1KOc mice (P<0.001) (
Figure 4D). Thus, the loss of Cdk1 not only extended the normal window of postnatal CM proliferation by 6d but also markedly impaired CMBN. This effect was dependent on the absence of Cdk1 and never observed in
controls. Intriguingly, the high numbers of dividing CM in hearts
Cdk1KOc mice at 0d are compatible with the concept, that Cdk1 is not absolutely required for prenatal CM proliferation.
2.5. Loss of Cdk1 Triggers Cell Cycle Re-Entry and Proliferation of Adult Cardiomyocytes Post-MI
To investigate the potential role of Cdk1 loss in CM proliferation and cardiac remodeling in a clinically relevant model of ischemic stress, 10-12 weeks old male and female
Cdk1KOc and
control mice were subjected to MI by permanent ligation of the Left Anterior Descending Artery. Next, RNA-seq analysis on total left ventricular (LV) tissues from
Cdk1KOc and
control mice was performed at 4 days post-MI for genome-wide mRNA profiling. Our previous study identified the 4-day post-MI time point as a critical ‘window’ that drives the proliferative cardiac phenotype
in vivo [
42]. Three-dimensional principal component analysis (3D-PCA) revealed distinct clustering of
Cdk1KOc and
control samples, both with and without MI, as analyzed using Strand-NGS software (
Figure 5A). Unsupervised hierarchical clustering at a high-confidence threshold identified 1,433 transcripts that significantly differed between
Cdk1KOc, and
control mice post-MI compared to sham-operated
controls (
Figure 5B,C). This finding was further supported by Volcano plot analysis, which showed that Cdk1 ablation led to the destabilization of mRNAs in
control versus
Cdk1KOc mice (
Figure 5D).
To explore gene clusters with similar biological functions in mutant and wild type transcriptomes, we performed Gene Set Enrichment Analysis (GSEA) at 4 days post-MI treatment. We found that transcripts related to ventricular remodeling, mitochondrial function, cardiac contraction, and mitosis were among the most enriched GSEA terms (
Figure 5E). Notably, transcripts for cell cycle-promoting factors, such as Mcm2, Mcm3, and Tert, were selectively enriched at 4 days post-MI (
Figure 5F). At 4d post-MI, the area at risk (AR) in
Cdk1KOc and
control groups was of equal size and involved greater than 80% of the cross-sectional area of the LV, measured at the level of the papillary muscles (Figure S7). We also noted a decreased sarcomeric α-actinin expression and lower extracellular matrix content, as evaluated by wheat germ agglutinin (WGA) staining. The border zone (BZ) showed α-actinin expression comparable to that of the area at risk with a markedly higher level of WGA reactivity. This is in stark contrast to the remote area (RA) where both α-actinin and WGA are observed at levels comparable to sham
Cdk1KOc and sham
controls (Figure S7, S8).
To further define the role of Cdk1 loss in cardiac remodeling after MI, DNA synthesis was assessed using EdU labeling in adult mice, followed by confocal immunofluorescence microscopy. Intriguingly, we observed the following numbers of S-phase CM in
Cdk1KOc mice at 4d post-MI: (1) AR 61±10.9 nuclei/mm
2, (2) BZ 33±7.6 CM nuclei/mm
2, (3) RA 1.2±3.4 CM nuclei/mm
2, and (4) total number (TN) 98±11.3 CM nuclei/mm
2 (
P<0.01) (
Figure 5G,H and Figure S7,S8). These results clearly demonstrate that induction of DNA synthesis in CM lacking Cdk1 requires an ischemic stress.
Assessment of numbers of CM in
Cdk1KOc mice in M-phase by immunofluorescence microscopy employing anti-phospho-histone 3 phosphorylated at H3.Pi-S28 showed: (1) AR 19.8±7.6 CM nuclei/mm
2, (2) BZ 9.7±3.2 CM nuclei/mm
2, (3) RA zero CM nuclei/mm
2, and (4) TN 30±8.6 CM nuclei/mm
2 (
P<0.01) (
Figure 5I,J,K and Figure S9). In all these experiments, EDU-positive and H3.Pi-S28-positive CM were absent from
control mice. In addition, the induction of CM in M-phase post-MI was restricted to CM residing in the BZ and was absent from the right ventricle.
Next, we analyzed whether CM from
Cdk1KOc mice completed cytokinesis by dividing into 2 daughter cells by detecting midbody structures in the final phase of daughter cell separation employing anti-AurkB antibodies. MI in
Cdk1KOc animals induced cytokinesis in 15.2±4.6 CM/mm
2 in the AR and minimally in the BZ (
P<0.001) (
Figure 5L,M,N and Figure S10) that was not observed in
controls. This supports the notion, that Cdk1 deletion restricts proliferation to a specific subpopulation located in an ischemic milieu. Of note, CM proliferation in hearts from sham
Cdk1KOc and sham
control mice was never observed.
Infarct ‘wound healing’ is characterized by initiation of non-CM (NCM) proliferation that mainly comprises endothelial cells and fibroblasts[
48]. We noticed a significant increase in the numbers of NCM in S- and M-phase in
Cdk1KOc and
control mice at 4d following MI (
P<0.001 vs. sham
Cdk1KOc and sham
controls) (
Figure 5G,J and Figure S7). Numbers of proliferating EDU-positive NCM was markedly lower in hearts of
Cdk1KOc mice compared to
controls (1506.1 vs 3519.6 NCM nuclei/mm
2). All these findings support the view that Cdk1 ablation induces adult CM proliferation after ischemic injury. We and others have previously shown that cell cycle re-entry in adult CM is tightly regulated by Cdk2[
42]. At 4d post MI, sham control hearts lacked an induction of Cdk2 protein relative to
Cdk1KOc post-MI (
Figure 5O). Elevated Cdk2 protein levels were accompanied by sustained induction of proliferative Cyclins A, B, D2, and E (
Figure 5O). Intriguingly, ablation of Cdk1 also led to extremely low protein levels of Cdk2-inhibitory p21 and p19, a cdk4/6 inhibitor.
We have established that adult
Cdk1KOc mice have developed a hypertrophic heart (
Figure 3I) that mainly consists of mononucleated CM (78±7.6%;
P<0.001 vs.
control) (
Figure 4D). Thus, we examined whether Cdk1-deficiency can further influence the ploidy of CM post-MI. We observed a small but significant increase (11%) of mononucleated CM, and a decrease (−75.4%) in binucleated CM in the AR/BZ areas in
Cdk1KOc post-MI (
P<0.001) (
Figure 5P). This data provide additional evidence that genetic Cdk1 ablation leads to CM proliferation rather than polyploidy.
2.6. Activation of EGFR1 Signaling Promotes Cardiomyocyte Proliferation in Cdk1KOc Mice Post-MI
RNA-seq analysis at 4d post-MI identified differential pathway activation between
Cdk1KOc and
control mice, as analyzed employing the pathway analysis module of Strand-NGS (
Figure 5Q). The EGFR1 signaling pathway was significantly upregulated in the
Cdk1KOc MI group, as indicated by the enrichment of 88 out of 159 genes within this pathway (
P< 0.001). This activation is associated with enhanced CM proliferation in the
Cdk1KOc mice post-MI. In contrast, the p53- and p38-MAPK signaling pathways were predominantly activated in the
control group. The activation of these pathways is known to inhibit CM division, suggesting that
control mice respond to MI with signals that block cell cycle re-entry[
49,
50]. In contrast, the absence of Cdk1 allows for the activation of proliferative pathways, such as EGFR1, potentially facilitating CM proliferation (
Figure 5Q)[
51]. The heatmap in
Figure 5R further corroborates the pathway analysis by displaying the expression profiles of key genes involved in the EGFR1, p38, and p53 signaling pathways across the four different experimental groups. Notably, key genes involved in the EGFR1 pathway, such as ErbB4, Erk1, and Sos2, show marked upregulation in the
Cdk1KOc MI group when compared to
controls, aligning with the observed pathway activation. Conversely, key genes associated with the p38-MAPK and p53 pathways, including p38α, Cdkn1, and Gadd45g, are upregulated in the control group post-MI, which is consistent with the inhibition of CM proliferation in this cohort (
Figure 5R). Western blot analysis confirms the activation of these pathways at the protein level (
Figure 5S). In
Cdk1KOc mice post-MI, there is an increase in phosphorylated ERK1/2 (Pi-ERK1/2), a key downstream effector of EGFR1 signaling, supporting the role of this pathway in promoting CM proliferation. In contrast, control mice show increased levels of p53 and phosphorylated p38α (Pi-38α), further validating the RNA-seq findings that these pathways are involved in blocking CM division following MI. The presence of these protein markers in the
Cdk1KOc and
control groups provides strong evidence for the differential activation of proliferative vs. inhibitory pathways depending on the genetic background of the experimental strains.
2.7. Cdk1-Deficiency Mitigates Hypertrophy, Boosts Survival and Improves Heart Function Post-MI
Next, we determined the physiological consequences of cardiac-specific Cdk1 knockout after ischemic injury. Cdk1 ablation significantly improve cardiac function in
Cdk1KOc mice post-MI versus sham
controls as assessed by echocardiographic measurement of FS at 21d post-MI (
Figure 6A). Measurements of LVEDD and LVESD dimensions revealed that
Cdk1KOc animals were protected from developing significant LV dilation post-MI that was not seen in hearts of sham
control mice (
Figure 6A).
Western blot analysis in
Figure 6B shows several key cardiac proteins, including Serca2a and RyR2, along with their phosphorylated activated forms. The
Cdk1KOc cohort demonstrates altered protein expression post-MI that correlates with improved calcium handling and contractility when compared to
controls.
Figure 6C shows levels of sarcomeric genes, which are crucial for cardiac muscle contraction, across different experimental groups. Notably, key genes such as Ttn, Myh7b and Tnnt1 are significantly upregulated in the
Cdk1KOc MI group relative to
controls. This elevated expression of sarcomeric genes in
Cdk1KOc mice post-MI suggests an enhanced contractile machinery, contributing to improved cardiac function. The heatmap in
Figure 6D shows the differential expression of genes involved in the cardiac cycle. Key genes such as Kcnj3, Scn4b, and Prkca, which are involved in cardiac rhythm and contraction regulation, are specifically upregulated in
Cdk1KOc mice post-MI. This upregulation implies that the
Cdk1 knockout promotes a more efficient cardiac cycle post-MI, likely contributing to the observed preservation of heart function in these mice.
Figure 6E illustrates the significant upregulation of transcript levels of various cardiac transcription factors, highlighting Gata4, Srf, and Nkx2-5 in
Cdk1KOc mice post-MI compared to
controls. These transcription factors are known to regulate genes critical for cardiac function and remodeling, and their increased expression in the knockout mice may play a pivotal role in the observed enhancement of cardiac function post-MI. Collectively, these findings indicate that the upregulation of sarcomeric and cardiac cycle genes, along with key transcription factors, contributes to the superior cardiac function observed in
Cdk1KOc mice post-MI.
Increases in HBW ratios in
Cdk1KOc mice were significantly less post-MI versus MI/
controls (
Figure 6F). Moreover, lung weight/body weight (LBW) ratios, an index of pulmonary congestion, was also less elevated in
Cdk1KOc mice as opposed to
controls post-MI (
Figure 6G). Importantly, Mantel-Cox test showed significantly improved survival in
Cdk1KOc at 21d post-MI relative to
controls (
P<0.001) (
Figure 6H). RT-qPCR analysis of canonical hypertrophic marker genes revealed decreased transcript levels ANP, BNP, and -MHC in
Cdk1KOc mice post-MI against
controls (
P<0.01) at 3d and 21d(
Figure 6I,J)). Moreover, CM cross-sectional area (
Figure 6K) and infarct size (
Figure 6L,M) as measured at the peri-infarct zone, were significantly lower in
Cdk1KOc mice in contrast to
controls post-MI (
P<0.01).
Inflammation typically precedes angiogenesis in the sequence of biological responses following MI[
3,
4,
52]. The heatmap in
Figure 6N shows the expression levels of inflammation-related genes across different experimental groups. In
Cdk1KOc mice post-MI, several key pro-inflammatory genes (Stat1, Il6ra, Nfkb1, Tgfbr1) are significantly upregulated against the
control cohort. This indicates a heightened inflammatory response in
Cdk1KOc mice following MI. Notably, genes involved in promoting angiogenesis (Vegfa, Flt1, Pecam1, Angpt1) are upregulated in the
Cdk1KOc mice relative to
controls post-MI (
Figure 6O). This suggests that
Cdk1KOc mice may have enhanced angiogenic responses post-MI, potentially contributing to better tissue repair and remodeling. Additionally, immunofluorescence images of fixed
Cdk1KOc-derived cardiac tissue samples exhibit significantly higher vWF expression post-MI (
P < 0.01) versus
controls, indicating increased endothelial cell activity and angiogenesis (
Figure 6P). All these findings indicate a higher degree of coronary angiogenesis in
Cdk1KOc mice post-MI as opposed to their wt counterparts.
Our heatmap analysis reveals that key apoptotic genes (Bax, Casp3, Fas, and Bcl2) are upregulated in the
controls in comparison to
Cdk1KOc animals post-MI (
Figure 6Q). Quantification of apoptotic CM through TUNEL staining and fluorescence microscopic analysis, demonstrates a significantly higher number of apoptotic CMs in
controls versus.
Cdk1KOc mice post-MI (
P<0.01) (
Figure 6R) suggesting that Cdk1 deletion protects against CM apoptosis post-MI.
2.8. Cdk1 Loss Preserves Mitochondrial Energetics by Protecting Against Ischemic Oxidative Stress
We have previously shown that MI-induced oxidative stress decreases in ATP levels negatively impact the contractile function and survival of CM[
42]. Thus, we analyzed gene expression of mt-DNA replication and mt-dynamics in
Cdk1KOc mice post-MI. The heatmap in
Figure 7A shows that transcripts levels of key genes involved in mt-DNA replication such as Tfam, PrimPol, and Tfb2m are upregulated in
Cdk1KOc mice post-MI. Notably, the expression of genes related to mt-fission and fusion like Fis1, Mfn2, and Mff are also significantly higher in the
Cdk1KOc mice post-MI (
Figure 7B). Next, mt-DNA copy number was quantified by qPCR using mt-Cytb that were normalized to the nuclear-encoded gene B2m. In
Cdk1KOc mice post-MI mice post-MI, we found a significant restoration of mt DNA copy number (P<0.01) (
Figure 7C) in contrast to
controls. This restoration was accompanied by an increase in mitochondrial respiratory capacity (
P<0.01) (
Figure 7D) and notably higher ATP/ADP ratios in the myocardium of
Cdk1KOc mice vs.
controls (
P<0.01) (
Figure 7E,F,G).
Moreover, we also observed significantly decreased levels of cytotoxic 4-HAE, an indicator of ROS-dependent lipid peroxidation[
27] (
Figure 7H,I) with upregulated glutathione/oxidized glutathione (GSH/GSSG) ratios in
Cdk1KOc mice post-MI (
Figure 7J,K,L). The heatmap in
Figure 7M illustrates transcript levels of genes involved in the pentose phosphate pathway (PPP), which is critical for generating glutathione, a key molecule for producing NADPH to combat oxidative stress[
53]. In
Cdk1KOc mice post-MI, there is an upregulation of genes such as Tkt and G6pdx
, indicating increased PPP activity compared to
controls. This upregulation correlates with higher NADPH/NADP+ ratios in
Cdk1KOc hearts subjected to MI (
Figure 7N) and implies an enhanced capacity to neutralize oxidative stress (
P< 0.01). This view agrees with our Western blot analysis showing elevated levels of antioxidative enzymes such as Sod2, Nqo1, and Pink1 in
Cdk1KOc mice post-MI, underscoring the improved oxidative stress response associated with Cdk1 knockout (
Figure 7O).