The Regulatory Role of Histone Modification on Gene Expression in the Early Stage of Myocardial Infarction

Myocardial infarction (MI) is a fatal heart disease with high morbidity and mortality. Various studies have demonstrated that a series of relatively specific biological events occur within 24 h of MI. However, the roles of histone modifications in this pathological process are still poorly understood. To investigate the regulation of histone modifications on gene expression in early MI, we performed RNA sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) on myocardial tissues 24 h after the onset of MI. The genome-wide profiles of five histone marks (H3K27ac, H3K9ac, H3K4me3, H3K9me3, and H3K27me3) were explored through ChIP-seq. RNA-seq identified 1,032 differentially expressed genes (DEGs) between the MI and sham groups. ChIP-seq analysis found that 195 upregulated DEGs were modified by change of at least one of the three active histone marks (H3K27ac, H3K9ac, and H3K4me3), and the biological processes and pathways analysis showed that these DEGs were significantly enriched in cardiomyocyte differentiation and development, inflammation, angiogenesis, and metabolism. In the transcriptional regulatory network, Ets1, Etv1, and Etv2 were predicted to be involved in gene expression regulation. In addition, by integrating super-enhancers (SEs) with RNA-seq data, 76 DEGs were associated with H3K27ac-enriched SEs in the MI group, and the functions of these SE-associated DEGs were mainly related to angiogenesis. Our results suggest that histone modifications may play important roles in the regulation of gene expression in the early stage of MI, and the early angiogenesis response may be initiated by SEs.

Introduction

Myocardial infarction (MI) is a leading cause of death all over the world, mainly caused by acute coronary artery occlusion, which leads to severe myocardial ischemia, and subsequent myocardium necrosis (1). A considerable number of biological events occur in the early stage of MI, such as inflammation, apoptosis, autophagy, and angiogenesis, leading to ventricular remodeling (2–4). Studies have demonstrated that myocardial pathological changes are mostly reversible, and intervention is more effective during this period (5, 6). Thus, it is of great significance to explore the key molecular events and regulatory mechanisms in the early stage of MI.

Histone modification is a post-translational modification including methylation, acetylation, phosphorylation, ubiquitylation, and SUMOylation, etc. (7), among which the methylation and acetylation of H3K4, H3K9, and H3K27 were common histone marks. Generally, acetylation correlates with transcriptional activation, such as acetylation of histone H3 on lys9 and lys27 (H3K9ac and H3K27ac) (8). Unlike acetylation, the transcriptional activity of methylation modifications on lysine residues varies according to the degree and site of methylation. Tri-methylation of histone H3 on lys9 and lys27 (H3K9me3 and H3K27me3) represents transcriptional suppression, whereas tri-methylation on lys4 (H3K4me3) is associated with transcriptional activation (9, 10).

Super-enhancers (SEs) are characterized as clusters of enhancers, which are enriched with transcription factors, cofactors, and epigenetic modifications, and considered to contain high transcription activity (11). Some enhancer-associated histone modification marks such as H3K27ac can be used to identify SEs (12). Emerging evidences suggest that SEs play important roles in differentiation, development, and tumorigenesis (13–15).

Previous studies had uncovered gene expression changes at the transcription level after MI. The expression of inflammatory response-related genes were reported to be significantly upregulated after MI, which resulted from the activation of NF-κB and TNF signaling pathways (16, 17). However, little is known about the roles of histone modifications and SEs in the early stage of MI. Here, by combining transcriptome and ChIP-seq of key histone modifiers, we set to reveal epigenetic patterns of differentially expressed genes and their roles in the early stage of MI.

Materials and Methods

Mouse MI Model

The animal experiments were approved by the Ethics Committee of Tongji University. The mouse model of MI was established by permanent ligation of the left anterior descending (LAD) coronary artery as described previously (18). Briefly, 8-week male C57BL/6 mice were subjected to LAD coronary artery ligation under anesthesia. Sham-operated mice of the same age underwent the same procedure without LAD occlusion. Twenty-four hours after operation, the infarct border zones were dissected, snap-frozen in liquid nitrogen, and stored at −80°C for subsequent experiments.

RNA Isolation

Myocardial tissues from the infarct border zones of MI and sham group were collected for RNA isolation. Five biological replicates were used for each group. Briefly, total RNA was extracted with 1 ml of RNAiso Plus reagent (Takara, 9108). After adding 200 μl of chloroform, the sample was mixed thoroughly and incubated at room temperature for 3 min, then centrifuged at 12,000 g for 15 min at 4°C. The supernatant was transferred to a new Eppendorf tube on ice, mixed with an equal volume of ice-cooled isopropanol, and then centrifuged at 12,000 g for 15 min at 4°C after placing on ice for 10 min. The supernatant was discarded, and the pellet was washed with 75% ice-cooled ethanol three times, and RNA was dissolved in DEPC water.

Chromatin Immunoprecipitation (ChIP) Assay

ChIP was carried out with the SimpleChIP^® Plus Sonication Chromatin IP kit (Cell Signaling Technology, 56383) according to the manufacturer's protocol. In brief, three biological replicates of myocardial tissues for each group from the infarct border zones were cross-linked with 1% formaldehyde for 10 min, then terminated by adding glycine. After disaggregating the samples into suspensions using Dounce homogenizer, the chromatins were isolated and sonicated to 100–300 bp fragments. Immunoprecipitation (IP) was performed with antibodies against H3K9me3 (Abcam, ab8898), H3K9ac (Active Motif, 39137), H3K4me3 (Abcam, ab8580), H3K27me3 (Cell Signaling Technology, 9733S), and H3K27ac (Active Motif, 39133). Finally, the chromatins were eluted, and then cross-links were reversed. The immunoprecipitated DNA and DNA input were purified using spin columns.

RNA Sequencing (RNA-Seq) and Analysis

The concentration and quality of RNA were measured by Qubit 4 Fluorometer (Invitrogen™, Q33238). For each group, two biological replicates of isolated RNA above (1 μg total RNA) with RIN > 8 were used for mRNA library construction. RNA sequencing libraries were generated with KAPA mRNA Hyper Prep kits (Roche, KK8581). The libraries were sequenced on an Illumina NovaSeq 6000 platform using a 2^* 150-bp sequencing kit for double-ended sequencing.

The raw data were first evaluated by FastQC (v0.11.3). Trimmomatic (v0.36) (19) was applied to remove adaptors, and information of clean reads in RNA-seq data was listed in Supplementary Table 1. Clean reads were then mapped to the mouse genome (mm10) using Hisat2 (20), and the mapping rates of clean reads are listed at Supplementary Table 2. The sequencing alignment data were evaluated by Qualimap (21), and the coverage of gene bodies is shown in Supplementary Figure 1. The counts of genes in each sample were calculated using HTSeq (v0.9.1) (22) and normalized to fragments per kilobase million (FPKM) for further visualization. The raw counts of genes were used to identify differentially expressed genes (DEGs) between sham and MI groups through DESeq2 R package (v1.36.0) (23) with a false discovery rate (FDR) P_FDR value < 0.05 and |Log₂ (Fold change)| ≥ 1. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted by R package clusterProfiler (v3.12.0) (24).

ChIP Sequencing (ChIP-Seq) and Analysis

VAHTSTM Universal DNA Library Prep Kit for Illumina^® V3 kit (Vazyme, ND607-02) was used to construct the ChIP-DNA library. Three biological replicates for MI group and two biological replicates for sham group were used for ChIP-seq library construction. The libraries were sequenced on an Illumina NovaSeq 6000 platform. Clean ChIP-seq reads were mapped to the mouse genome (mm10) using bowtie2 (25). To identify the H3K4me3, H3K9me3, H3K9ac, H3K27me3, and H3K27ac peaks, peaks of enriched occupancy relative to a background input were called using MACS 1.4 (26). The peaks were then annotated by the R package ChIPseeker (27) (v1.22.1). Peaks of each histone modification were merged by bedtools (v2.29.2) (28); unique peaks for each group were identified by MAnorm (v1.2.0) (29) software and confirmed by Diffbind R package (v2.14.0) (30). Diffbind R package was used to do the clustering. Hypergeometric optimization of motif enrichment (HOMER) (31) was used to discover motifs of specific regions. The bam files were converted into bigwig files using bedtools (v2.29.2) (28), and then bigwig files were visualized through deeptools (v3.3.0) (32) to generate metagene plots.

Quantitative Real-Time Quantitative Polymerase Chain Reaction (qRT-PCR) and ChIP-qPCR

Total RNA was reversely transcribed into cDNA with the PrimeScript RT reagent kit (Takara, RR036Q, Japan). qRT-PCR assays were performed in five biological replicates. The reaction volume of real-time qPCR is 20 μl, containing 10 μl of SYBR Green Master Mix (ABI, 4309155, USA), 4 μl cDNA, and 1 μl forward and reverse primers, the primers are listed in Supplementary Table 3. The relative gene expression was calculated with the 2^−ΔΔCt method. ChIP-qPCR was performed with 20 μl of reaction volume, containing 10 μl of SYBR Green Master Mix, 1 μl of forward and reverse primers, 2 μl of DNA, and 6 μl of H₂O. ChIP DNA enrichment was calculated as % of input, and the primers are listed in Supplementary Table 3. All qPCR assays were performed in three biological replicates.

Identification of Typical and Super-Enhancers

Ranking of super-enhancer (ROSE) (11) was applied to call typical enhancers (TEs) and super-enhancers (SEs). In brief, H3K27ac peaks with the distance shorter than 12.5 kb were merged and ranked by the H3K27ac signal. The inflection point with the tangent line slope equals 1 was the dividing point of TEs and SEs. Enhancers above or below the point were defined as SEs and TEs, respectively. The genes that were closest to the enhancers were considered as enhancer-associated genes. For the expression of enhancer-associated genes, the Wilcoxon rank sum test was performed to test statistical significance between two groups.

Results

Transcriptomic Analysis and Histone Modification Profiling

Adult C57/BL6 mice (8 weeks old) were subjected to ligation of left anterior descending (LAD) or sham surgery randomly. The mice were sacrificed 24 h after the surgery, and the left ventricular tissues in the infarct border zone and normal tissues in the corresponding position were dissected to examine the expression profile of MI and sham groups by RNA-seq. The quality of the clean reads and mapping rates of RNA-seq are shown in Supplementary Tables 1, 2, and Supplementary Figure 1. Differentially expressed genes (DEGs), 1,032, between the MI and sham groups were identified by RNA-seq with an absolute log-fold change over or equal to 1 (adjusted P-value < 0.05) (Figure 1B, Supplementary Table 4). Among them, 423 (40.99%) and 609 (59.01%) genes were upregulated and downregulated, respectively, in the MI group. DEGs, 428 (41.47%), in our RNA-seq data, were overlapped with the DEGs from cardiac organoid model of human myocardial infarction (GSE113871, Supplementary Figure 2). Moreover, the expression of representative cardiac remodeling-related genes (Nppa and Nppb) was significantly increased in the MI group (Supplementary Figure 3A), which demonstrated the credibility of our MI models. We further performed qRT-PCR to validate the mRNA expression of Nppa, Nppb, and Myh7 in samples both 24 and 48 h after MI. The expression of Nppa and Nppb increased in 24 h after MI, and then decreased in 48 h after MI (Supplementary Figure 3B). The upregulated genes were enriched in the function of muscle structure development, angiogenesis, positive regulation of cell death, muscle cell development, and ribosomal large subunit biogenesis. In contrast, downregulated genes were enriched in the establishment of organelle localization, oocyte meiosis, positive regulation of potassium ion transmembrane transport, flavonoid glucuronidation, and mitotic cell cycle process (Figure 1C).

FIGURE 1

Figure 1. Transcriptome analysis and histone modification profiling in the early stage of myocardial infraction (MI), (A) Schematic illustration of the present study. (B) Volcano plots showing the differentially expressed genes (DEGs) between sham and MI groups. Blue dots represent downregulated genes, and red dots represent upregulated genes. (C) Gene ontology (GO) analysis of the differentially expressed genes (DEGs) in RNA-seq. Orange bras represent the enrichment of upregulated genes; blue bars represent the enrichment of downregulated genes. (D) Hierarchical clustering of H3K9ac, H3K9me3, H3K27ac, H3K27me3, and H3K4me3 occupancy at genome wide in sham and MI group. Spearman correlation coefficient was used to calculate sample correlations. (E) The genomic distribution of different histone modification regions in sham and MI groups.

To profile epigenetic changes in the early stage of MI, we selected H3K27ac, H3K9ac, and H3K4me3 (three marks associated with active regulatory regions) and H3K9me3 and H3K27me3 (two marks associated with repressive regulatory regions) antibodies to perform ChIP-Seq (Figure 1A). Hierarchical clustering demonstrated that these histone marks were mainly divided into three clusters (Figure 1D). While H3K9me3 and H3K27me3 (two marks associated with repressive regulatory regions) formed two separate clusters, interestingly, all the marks associated with active regulatory regions (H3K27ac, H3K9ac, and H3K4me3) formed a big cluster, suggesting that the active regions might be modified by multiple active histone marks. We then quantified peaks of different histone marks in genomic regions, the number of which was similar between MI and sham groups (Figure 1D). The numbers of peaks of different histone marks and the statistical information of sequencing reads are shown in Table 1 and Supplementary Figure 4. The distribution of genomic regions modified by H3K27ac, H3K27me3, H3K4me3, H3K9ac, and H3K9me3 was classified into six regions (promoter, exon, 3′ UTR, 5′ UTR, intron, and intergenic). The H3K4me3 and H3K9ac marks were mainly mapped to promoters; H3K27ac and H3K27me3 marks were mapped to promoters, intron, and intergenic regions, and H3K9me3 marks were mostly mapped to intergenic regions (Figure 1E), which were consistent with the reported distribution modes of these histone marks (33).

TABLE 1

Table 1. Information of the peaks called in different samples.

Epigenetic Profile Changes

To reveal the difference in histone modifications between MI and sham, we drew the average signals of histone marks relative to ±3 kb of the transcription start site (TSS) (Figure 2A). The average signals of active histone marks (H3K27ac, H3K4me3, and H3K9ac) were all higher in MI than those in the sham group. The signals of H3K9me3 in MI group were relatively lower than those in the sham group, while no apparent changes in H3K27me3 were observed between sham and MI. We identified 3,986 MI unique peaks for H3K27ac mark, 5,417 peaks for H3K27me3 mark, 1,910 peaks for H3K4me3 mark, 12,084 peaks for H3K9ac mark, and 937 peaks for H3K9me3 mark in the TSS ± 3 kb (Figure 2B and Supplementary Tables 5–9). We then performed gene ontology (GO) analysis to investigate the biological function of these MI unique peaks of each histone mark and discovered that the MI unique peaks of H3K27ac (2,778 mapped genes) were highly enriched in small GTPase-mediated signal transduction, H3K27me3 (3,325 mapped genes) in cell-substrate adhesion, H3K4me3 (1,500 mapped genes) in the process of inflammatory response, and H3K9ac (6,547 mapped genes) in small GTPase mediated signal transduction (Figure 2C).

FIGURE 2

Figure 2. Epigenetic profile change in the early stage of MI. (A) The average signals of each histone mark relative to ±3 kb of the transcription start site (TSS) between the sham and MI group. (B) The number of MI unique peaks in promoters of H3K9ac, H3K9me3, H3K27ac, H3K27me3, and H3K4me3 marks. (C) GO analysis of the MI unique peaks in promoters of H3K9ac, H3K27ac, H3K27me3, and, H3K4me3 marks.

Active Chromatin Regions Associated With the Upregulated Genes

Notably, there were few overlapped genes between downregulated genes in RNA-seq and unique peaks in repressive histone marks (Supplementary Figure 5B). Thus, we mainly study the regulation of histone modification on upregulated genes; we intersected the 423 activated genes in RNA-seq [Log2 (fold change) ≥ 1, adjusted P-value < 0.05] with 7,226 open regions (the gain of active marks H3K9ac, H3K27ac, and H3K4me3 in the ±2 kb from TSS). A total of 195 activated genes were enriched for at least one of H3K9ac (145/195), H3K27ac (79/195), or H3K4me3 (59/195) marks (Figure 3A, Supplementary Figure 6, and Supplementary Table 10). According to our pipeline, we further validated the 195 activated genes revealed by another method (Diffbind), which also identifies differential peaks. Of the genes, 42.05% (82/195) were shared by the two methods (Supplementary Figure 5A and Supplementary Table 11). We also checked the other 113 genes, which were not identified by Diffbind in genome browser. Most of the promoter regions of the 113 genes did have more active modifications in MI group than that in sham group. We then performed GO analysis to categorize the biological function of the 195 upregulated DEGs and found that they were mainly enriched in muscle cell-related processes (striated muscle cell development, muscle cell differentiation, muscle cell development), inflammatory response processes (leukocyte migration, negative regulation of response to external stimulus, myeloid leukocyte migration, leukocyte chemotaxis) and angiogenesis-related processes (regulation of vasculature development, positive regulation of vasculature development, and regulation of angiogenesis) (Figure 3B). KEGG analysis indicated that inflammatory-related signaling pathways, including TNF signaling pathway, IL-17 signaling pathway, and cytokine–cytokine receptor interaction; cell signal-related pathways (PI3K-Akt signaling pathway and MAPK signal pathway); ECM–receptor interaction and HIF-1 signaling pathways were activated in the process of MI (Figure 3C). The inflammatory response processes and pathways were revealed by enrichment analysis, and inflammatory response-related genes are shown in Figure 3D. Two typical genes (Cxcl2 and Cxcl3) are shown in Figures 3E,F, with increasing H3K4me3, H3K9ac, and H3K27ac signals at the promoter or increasing H3K27ac signals at the enhancer region of Cxcl2 and Cxcl3 loci in MI group. ChIP-qPCR and RT-qPCR verified elevated H3K27ac modification in promoter regions and elevated expression of Cxcl2 and Cxcl3, respectively (Figures 3G,H).

FIGURE 3

Figure 3. The upregulated genes associated with gain of active histone marks. (A) Venn diagram showing upregulated DEGs associated with the open regions in the promoters (gain of H3K9ac, H3K27ac, and H3K4me3 marks in the ±2 kb from TSS). (B) GO analysis of the upregulated DEGs enriched in promoters (±2 kb from TSS) of H3K9ac, H3K27ac, and H3K4me3 marks. (C) Kyoto Encyclopedia of Genes and Genomes (KEGGs) pathway enrichment analysis of the upregulated DEGs enriched in promoters (±2 kb from TSS) of H3K9ac, H3K27ac, and H3K4me3 marks. (D) Heatmap showing the relative expression fragments per kilobase million (FPKM) of genes related to inflammatory response processes. (E,F) The increased active histone marks at promoters and enhancers of Cxcl2 and Cxcl3 in MI compared with the sham group. (G) ChIP-qPCR validation of the distribution of H3K27ac changes in Cxcl2 and Cxcl3. (H) qPCR validation of the expression of Cxcl2 and Cxcl3. (I) The transcription factors and binding motifs identified at promoter regions of upregulated DEGs enriched for H3K9ac, H3K27ac, and H3K4me3. n = 3; *P < 0.05, ***P < 0.001.

Moreover, to identify the transcription factors, which may regulate the upregulated genes in MI, we performed motif analysis at the promoter regions of unique peaks of three active marks (H3K27ac, H3K9ac, and H3K4me3) in the MI group. Hypergeometric optimization of motif enrichment (HOMER) was employed to identify top enriched motifs, which showed that motifs for transcription factors including ETS1 (P-value = 1e-15), ETV1 (P-value = 1e-14), and ETV2 (P-value = 1e-11) were significantly enriched in promoters of upregulated genes associated with active histone marks (Figure 3I and Supplementary Table 12).

Identification of MI-Specific SEs

The above data showed that the signal of H3K27ac was notably higher in the early stage of MI (Figure 2A). H3K27ac is also an active TE and SE marker; however, little is known about the function of SEs in the early stage of MI. We identified 20,337 TEs and 1,672 SEs in the sham group, and 21,444 TEs and 1,723 SEs in the MI group (Figure 4A), showing more enhancers activated in the MI group. Many SE-associated genes, such as Per2, Adamts1, Ccl2, Pgf , Thbs1, and Tgfb2, were known to be related to angiogenesis. Expression levels of SE-associated transcripts in MI were significantly higher than those of TE-associated transcripts (Figure 4B). SEs, 1,723, in the MI group were mapped to 1,565 SE-associated genes, and 76 SEs in the MI group were in accordance with DEGs, which can be further categorized into 32 MI unique SEs and 44 common SEs (Figure 4C). We then performed GO analysis for the 76 SEs associated DEGs and found that about half of the differentially expressed SE-associated genes were enriched in the function of angiogenesis (Figure 4D), and the specific genes in the angiogenesis-related GO terms are listed in Figure 4E. In addition, SEs that were not connected with DEGs were mainly enriched in actin filament-based process (actin filament organization, muscle organ development, regulation of actin filament-based process, lamellipodium organization, and assembly) and cell migration (Supplementary Figure 7). Two representative MI SE-associated genes Itgb2 and Ccl2 are shown in Figure 4F, whose H3K27ac occupancy of some component enhancers in the MI group was higher than that in the sham group. Itgb2- and Ccl2-associated SEs can also be validated on the human umbilical vein endothelial cells (HUVEC) and human skeletal muscle myoblasts (HSMM) from ENCODE database (Figures 4G,H). Our analysis identified a set of SEs, which might be essential for the early response of MI and suggested that MI-induced global changes in SEs might coordinate with angiogenesis.

FIGURE 4

Figure 4. Identification of MI specific super-enhancers (SEs). (A) Ranked plots of typical enhancers (TEs) and super enhancers (SEs) in the sham and MI group by increasing H3K27ac signal (normalized by the ratio to the highest H3K27ac reads of the enhancer). The axes values represent the ranks of enhancers by H3K27ac signal. Angiogenesis-related genes are marked out in the figures. (B) The expression (Log₂) of TE-associated genes vs. SE-associated genes in the MI group. P-value was determined by the Wilcoxon test. (C) Overlapping genes of DEGs and SE-associated genes in the MI group. The MI unique SE-associated genes and common SE-associated genes are listed. (D) GO analysis of the 76 overlapping genes between DEGs and SE-associated genes in the MI group. (E) The specific genes in the angiogenesis-related GO terms. (F) The H3K27ac occupancy of two representative angiogenesis-related genes Itgb2 and Ccl2 is shown. The red and gray peaks, respectively, mean MI and sham group. (G,H) Itgb2- and Ccl2-associated SEs on human umbilical vein endothelial cells (HUVECs) and human skeletal muscle myoblasts (HSMM) from ENCODE project.

Discussion

In this study, we correlated histone modifications to gene expression in the early stage of MI via integrative analysis of the transcriptome and ChIP-seq of histone marks. The biological processes of the upregulated DEGs modified by active marks (H3K9ac, H3K27ac, and H3K4me3) were enriched in cardiomyocyte differentiation and development, inflammation, and angiogenesis. Pathways implicated in the pathological processes of MI were activated, such as inflammation and metabolism-related pathways. Transcription factors, including ETS1, ETV1, and ETV2, may be involved in the transcriptional regulation of molecular events in the early stage of MI. Moreover, a set of SE-associated genes related to angiogenesis were also uncovered.

It is worth noting that inflammatory response was enriched in both function (GO) and pathway (KEGG) analysis, and H3K9ac, H3K27ac, and H3K4me3 regulated the expression of critical chemokines in the early stage of MI. After the occurrence of MI, chemokines were generated in injured myocardium, bound to the endothelial surface and extracellular matrix through glycosaminoglycans, and interacted with receptors on the surface of leukocytes to promote trafficking (34). CC and CXC, two main subtypes of chemokines, were significantly increased in infarcted tissue, which mediated leukocyte migration into the injured myocardium after MI (35, 36). Chemokines played a vital role in the inflammatory response and were associated with cardiac injury, repair, and remodeling following MI (36, 37). Our study showed that the expression of Ccl2, Ccl6, Cxcl2, Cxcl3, was upregulated and modified by H3K9ac, H3K27ac, and H3K4me3 marks. Ccl2 and Cxcl2 were two widely studied chemokines, which belong to the CC and CXC subtypes. Ccl2 played an important role in recruiting and activating mononuclear cells, and inhibiting Ccl2 reduced the infarct size and ameliorated the cardiac remolding after MI (38, 39). Meanwhile, Cxcl2 mediated the recruitment of neutrophils in the inflammatory response and is involved in the occurrence and development of many cardiovascular diseases, including atherosclerosis, ischemic stroke, and myocardial infarction (40). Studies have demonstrated that highly expressed Cxcl2 may exacerbate myocardial injury and inhibition of Cxcl2 reduced neutrophil-mediated tissue injury and infarct size after MI (41–43). Our results revealed that these active histone marks increased at the promoter regions and may promote the expression of these key chemokines after MI. Moreover, biological processes related to cardiomyocyte differentiation and development were also enriched. Among them, Csrp3, Pdlim5, Sorbs2, Rcan1, and Acta1, which were related to cardiac hypertrophy and remodeling, were upregulated and modified by at least one active histone mark. Altogether, these results indicated that H3K9ac, H3K27ac, and H3K4me3 modifications might participate in cardiac remodeling by regulating the inflammatory response.

In order to get the transcription factors that coordinated the transcriptional regulatory network in the early stage of MI, we performed motif analysis and identified ETS1, ETV1, and ETV2. All three transcription factors were members of the ETS family. ETS1 was associated with fibrotic remodeling of the myocardium (44). ETV1 was involved in atrial remodeling and arrhythmia (45). ETV2 was reported to promote angiogenesis and improve myocardial repair after MI (46). These transcription factors may work with active histone modifications to regulate gene expression in the early stage of MI.

Super enhancers (SEs) were first identified as a large cluster of enhancers with the strong binding of transcription factors that drive the expression of cell identity genes (11). It has been reported that SEs mediate angiotensin II-induced vascular smooth muscle cell dysfunction (47), whereas their roles in MI were less known. In our study, we investigated the function of SEs in MI. We intersected SEs with DEGs and found that 76 DEGs may be regulated by SEs. These SEs can be further categorized into 32 MI unique SEs and 44 common SEs. We then performed GO analysis for the 76 DEGs and found that these genes were significantly enriched in angiogenesis. Angiogenesis after MI is considered to be an important factor leading to ventricular remodeling (48, 49), and neovascularization induced by angiogenic growth factors reduces the infarct size and improves ventricular remodeling (50, 51). Our results showed that SEs might mediate the transcription of angiogenesis-related genes at the border zone in the early stage of MI, consistent with the previous report that angiogenesis initiated in the border zone at the early stage post-MI (52). Moreover, some of these genes have been reported to be associated with angiogenesis in MI, such as Pgf and Ccl2. Pgf (placental growth factor) was reported to promote angiogenesis and myocardial repair post-MI (53, 54). Ccl2 may be implicated in angiogenesis in ischemic myocardial tissue by attracting macrophages capable of producing angiogenic factors and by exerting direct pro-angiogenesis effects through endothelial cells (36, 55). Apart from SEs, which were associated with DEGs, SEs that were not connected with DEGs were significantly enriched in actin filament-based process and cell migration. These SE-associated genes may remain in a poised status in early MI and participate in the process of cardiac remodeling in the later stage of MI. A representative gene with broad and high levels of H3K27ac signal is Itgb2, which mediated the migration and recruitment of endothelial progenitor cells and angiogenesis in the infarcted myocardium (56). Taken together, the transcriptional regulation of angiogenesis-related genes by SEs in the early stage of MI may initiate vascular regeneration and promote repair of damaged myocardium post-MI, which may prevent pathological cardiac remodeling.

Our study also has limitations. We mainly focused on promoters and SEs, while other regulatory elements of gene expression in the early stage of MI were not investigated, which might lead to underestimation of the roles of histone modifications in the early stage of MI. Besides, we only selected five histone marks to study the regulation of histone modification on gene expression in MI; the roles of other histone marks in MI were not investigated. Moreover, the results were based on a mouse model of MI, so further observation and validation in human samples are needed.

In summary, our study revealed the histone modification profile in the early stage of MI in mice and demonstrated that the modulation of histone modifications (H3K9ac, H3K27ac, and H3K4me3) might be involved in inflammation and angiogenesis by regulating promoters and SEs, and participating in the pathological processes of cardiac remodeling.

Data Availability Statement

All the sequencing results have been uploaded and deposited in the NCBI Sequence Read Archive database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA657342) under accession No. SRP277594.

Ethics Statement

The animal study was reviewed and approved by Ethics Committee of Tongji University.

Author Contributions

Y-HC and DS conceived and coordinated the study. YZ and JW constructed the MI models of mice and collected the samples. JW, YZ, LN, and LH constructed the library for RNA-seq and ChIP-seq and performed some experiments. BL made the bioinformatic analysis of RNA-seq and ChIP-seq and drew the figures. JW drafted the manuscript. JY and LX conducted and gave advice on the study. Y-HC, JY, DS, and BL reviewed and revised the manuscript. All authors have approved the final version of the manuscript.

Funding

This work was funded by the Grants from the National Key Research and Development Plan (2019YFA0801501 to Y-HC, 2018YFC1312504 to LX), the Programs of National Natural Science Foundation of China (81530017, 81930013, and 81770397 to Y-HC), the Key Disciplines Group Construction Project of Pudong Health Bureau of Shanghai (PWZxq2017-05), the Top-level Clinical Discipline Project of Shanghai Pudong District (PWYgf2018-02), the Research Unit of Origin and Regulation of Heart Rhythm, and the Chinese Academy of Medical Sciences (2019RU045). Y-HC was a Fellow at the Collaborative Innovation Center for Cardiovascular Disease Translational Medicine, Nanjing Medical University.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Feizhen Wu for comments on the data analyses and insightful scientific discussions.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2020.594325/full#supplementary-material

References

1. Anderson JL, Morrow DA. Acute myocardial infarction. N Engl J Med. (2017) 376:2053–64. doi: 10.1056/NEJMra1606915

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Kanamori H, Takemura G, Goto K, Maruyama R, Tsujimoto A, Ogino A, et al. The role of autophagy emerging in postinfarction cardiac remodelling. Cardiovasc Res. (2011) 91:330–9. doi: 10.1093/cvr/cvr073

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Wang X, Guo Z, Ding Z, Mehta JL. Inflammation, autophagy, and apoptosis after myocardial infarction. J Am Heart Assoc. (2018) 7:e008024. doi: 10.1161/JAHA.117.008024

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Xiang F, Shi Z, Guo X, Qiu Z, Chen X, Huang F, et al. Proteomic analysis of myocardial tissue from the border zone during early stage post-infarct remodelling in rats. Eur J Heart Fail. (2011) 13:254–63. doi: 10.1093/eurjhf/hfq196

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Opie LH, Commerford PJ, Gersh BJ, Pfeffer MA. Controversies in ventricular remodelling. Lancet. (2006) 367:356–67. doi: 10.1016/S0140-6736(06)68074-4

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Cohn JN, Ferrari R, Sharpe N. Cardiac remodeling–concepts and clinical implications: a consensus paper from an international forum on cardiac remodeling. Behalf of an international forum on cardiac remodeling. J Am Coll Cardiol. (2000) 35:569–82. doi: 10.1016/S0735-1097(99)00630-0

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Bannister AJ, Kouzarides T. Regulation of chromatin by histone modifications. Cell Res. (2011) 21:381–95. doi: 10.1038/cr.2011.22

CrossRef Full Text | Google Scholar

8. Shahbazian MD, Grunstein M. Functions of site-specific histone acetylation and deacetylation. Annu Rev Biochem. (2007) 76:75–100. doi: 10.1146/annurev.biochem.76.052705.162114

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Black JC, Van Rechem C, Whetstine JR. Histone lysine methylation dynamics: establishment, regulation, biological impact. Mol Cell. (2012) 48:491–507. doi: 10.1016/j.molcel.2012.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Kouzarides T. Chromatin modifications and their function. Cell. (2007) 128:693–705. doi: 10.1016/j.cell.2007.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Whyte WA, Orlando DA, Hnisz D, Abraham BJ, Lin CY, Kagey MH, et al. Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell. (2013) 153:307–19. doi: 10.1016/j.cell.2013.03.035

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Hnisz D, Abraham BJ, Lee TI, Lau A, Saint-Andre V, Sigova AA, et al. Super-enhancers in the control of cell identity and disease. Cell. (2013) 155:934–47. doi: 10.1016/j.cell.2013.09.053

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sun X, Ren Z, Cun Y, Zhao C, Huang X, Zhou J, et al. Hippo-YAP signaling controls lineage differentiation of mouse embryonic stem cells through modulating the formation of super-enhancers. Nucleic Acids Res. (2020) 48:7182–96. doi: 10.1093/nar/gkaa482

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Wang X, Cairns MJ, Yan J. Super-enhancers in transcriptional regulation and genome organization. Nucleic Acids Res. (2019) 47:11481–96. doi: 10.1093/nar/gkz1038

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Thandapani P. Super-enhancers in cancer. Pharmacol Ther. (2019) 199:129–38. doi: 10.1016/j.pharmthera.2019.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Li Y, Wang C, Li T, Ma L, Fan F, Jin Y, et al. The whole transcriptome and proteome changes in the early stage of myocardial infarction. Cell Death Discov. (2019) 5:73. doi: 10.1038/s41420-019-0152-z

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zhao Q, Wu K, Li N, Li Z, Jin F. Identification of potentially relevant genes for myocardial infarction using RNA sequencing data analysis. Exp Ther Med. (2018) 15:1456–64. doi: 10.3892/etm.2017.5580

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Li L, Weng Z, Yao C, Song Y, Ma T. Aquaporin-1 deficiency protects against myocardial infarction by reducing both edema and apoptosis in mice. Sci Rep. (2015) 5:13807. doi: 10.1038/srep13807

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics. (2014) 30:2114–20. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. (2015) 12:357–60. doi: 10.1038/nmeth.3317

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Garcia-Alcalde F, Okonechnikov K, Carbonell J, Cruz LM, Gotz S, Tarazona S, et al. Qualimap: evaluating next-generation sequencing alignment data. Bioinformatics. (2012) 28:2678–9. doi: 10.1093/bioinformatics/bts503

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Anders S, Pyl PT, Huber W. HTSeq–a python framework to work with high-throughput sequencing data. Bioinformatics. (2015) 31:166–9. doi: 10.1093/bioinformatics/btu638

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. (2012) 16:284–7. doi: 10.1089/omi.2011.0118

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. (2012) 9:357–9. doi: 10.1038/nmeth.1923

CrossRef Full Text | Google Scholar

26. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. (2008) 9:R137. doi: 10.1186/gb-2008-9-9-r137

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Yu G, Wang LG, He QY. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics. (2015) 31:2382–3. doi: 10.1093/bioinformatics/btv145

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. (2010) 26:841–2. doi: 10.1093/bioinformatics/btq033

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Shao Z, Zhang Y, Yuan GC, Orkin SH, Waxman DJ. MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets. Genome Biol. (2012) 13:R16. doi: 10.1186/gb-2012-13-3-r16

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ross-Innes CS, Stark R, Teschendorff AE, Holmes KA, Ali HR, Dunning MJ, et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature. (2012) 481:389–93. doi: 10.1038/nature10730

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. (2010) 38:576–89. doi: 10.1016/j.molcel.2010.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Ramirez F, Ryan DP, Gruning B, Bhardwaj V, Kilpert F, Richter AS, et al. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. (2016) 44:W160–5. doi: 10.1093/nar/gkw257

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Papait R, Cattaneo P, Kunderfranco P, Greco C, Carullo P, Guffanti A, et al. Genome-wide analysis of histone marks identifying an epigenetic signature of promoters and enhancers underlying cardiac hypertrophy. Proc Natl Acad Sci USA. (2013) 110:20164–9. doi: 10.1073/pnas.1315155110

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Jung S, Littman DR. Chemokine receptors in lymphoid organ homeostasis. Curr Opin Immunol. (1999) 11:319–25. doi: 10.1016/S0952-7915(99)80051-X

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Frangogiannis NG Chemokines in ischemia and reperfusion. Thromb Haemost. (2007) 97:738–47. doi: 10.1160/TH07-01-0022

CrossRef Full Text | Google Scholar

36. Chen B, Frangogiannis NG. Chemokines in myocardial infarction. J Cardiovasc Transl Res. (2020). doi: 10.1007/s12265-020-10006-7. [Epub ahead of print].

CrossRef Full Text | Google Scholar

37. Cavalera M, Frangogiannis NG. Targeting the chemokines in cardiac repair. Curr Pharm Des. (2014) 20:1971–9. doi: 10.2174/13816128113199990449

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Liehn EA, Piccinini AM, Koenen RR, Soehnlein O, Adage T, Fatu R, et al. A new monocyte chemotactic protein-1/chemokine CC motif ligand-2 competitor limiting neointima formation and myocardial ischemia/reperfusion injury in mice. J Am Coll Cardiol. (2010) 56:1847–57. doi: 10.1016/j.jacc.2010.04.066

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Dewald O, Zymek P, Winkelmann K, Koerting A, Ren G, Abou-Khamis T, et al. CCL2/Monocyte chemoattractant protein-1 regulates inflammatory responses critical to healing myocardial infarcts. Circ Res. (2005) 96:881–9. doi: 10.1161/01.RES.0000163017.13772.3a

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Guo LY, Yang F, Peng LJ, Li YB, Wang AP. CXCL2, a new critical factor and therapeutic target for cardiovascular diseases. Clin Exp Hypertens. (2020) 42:428–37. doi: 10.1080/10641963.2019.1693585

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Liehn EA, Postea O, Curaj A, Marx N. Repair after myocardial infarction, between fantasy and reality: the role of chemokines. J Am Coll Cardiol. (2011) 58:2357–62. doi: 10.1016/j.jacc.2011.08.034

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Montecucco F, Bauer I, Braunersreuther V, Bruzzone S, Akhmedov A, Luscher TF, et al. Inhibition of nicotinamide phosphoribosyltransferase reduces neutrophil-mediated injury in myocardial infarction. Antioxid Redox Signal. (2013) 18:630–41. doi: 10.1089/ars.2011.4487

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Mylonas KJ, Turner NA, Bageghni SA, Kenyon CJ, White CI, McGregor K, et al. 11beta-HSD1 suppresses cardiac fibroblast CXCL2, CXCL5 neutrophil recruitment to the heart post MI. J Endocrinol. (2017) 233:315–27. doi: 10.1530/JOE-16-0501

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Hao G, Han Z, Meng Z, Wei J, Gao D, Zhang H, et al. Ets-1 upregulation mediates angiotensin II-related cardiac fibrosis. Int J Clin Exp Pathol. (2015) 8:10216–27.

PubMed Abstract | Google Scholar

45. Rommel C, Rosner S, Lother A, Barg M, Schwaderer M, Gilsbach R, et al. The transcription factor ETV1 induces atrial remodeling and arrhythmia. Circ Res. (2018) 123:550–63. doi: 10.1161/CIRCRESAHA.118.313036

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Lee S, Lee DH, Park BW, Kim R, Hoang AD, Woo SK, et al. In vivo transduction of ETV2 improves cardiac function and induces vascular regeneration following myocardial infarction. Exp Mol Med. (2019) 51:1–14. doi: 10.1038/s12276-019-0206-6

CrossRef Full Text | Google Scholar

47. Das S, Senapati P, Chen Z, Reddy MA, Ganguly R, Lanting L, et al. Regulation of angiotensin II actions by enhancers and super-enhancers in vascular smooth muscle cells. Nat Commun. (2017) 8:1467. doi: 10.1038/s41467-017-01629-7

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Saraste A, Koskenvuo JW, Saraste M, Parkka J, Toikka J, Naum A, et al. Coronary artery flow velocity profile measured by transthoracic doppler echocardiography predicts myocardial viability after acute myocardial infarction. Heart. (2007) 93:456–7. doi: 10.1136/hrt.2006.094995

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Oostendorp M, Douma K, Wagenaar A, Slenter JM, Hackeng TM, van Zandvoort MA, et al. Molecular magnetic resonance imaging of myocardial angiogenesis after acute myocardial infarction. Circulation. (2010) 121:775–83. doi: 10.1161/CIRCULATIONAHA.109.889451

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Yanagisawa-Miwa A, Uchida Y, Nakamura F, Tomaru T, Kido H, Kamijo T, et al. Salvage of infarcted myocardium by angiogenic action of basic fibroblast growth factor. Science. (1992) 257:1401–3. doi: 10.1126/science.1382313

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Garbern JC, Minami E, Stayton PS, Murry CE. Delivery of basic fibroblast growth factor with a pH-responsive, injectable hydrogel to improve angiogenesis in infarcted myocardium. Biomaterials. (2011) 32:2407–16. doi: 10.1016/j.biomaterials.2010.11.075

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Zhao W, Zhao T, Chen Y, Ahokas RA, Sun Y. Reactive oxygen species promote angiogenesis in the infarcted rat heart. Int J Exp Pathol. (2009) 90:621–9. doi: 10.1111/j.1365-2613.2009.00682.x

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Iwama H, Uemura S, Naya N, Imagawa K, Takemoto Y, Asai O, et al. Cardiac expression of placental growth factor predicts the improvement of chronic phase left ventricular function in patients with acute myocardial infarction. J Am Coll Cardiol. (2006) 47:1559–67. doi: 10.1016/j.jacc.2005.11.064

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Luttun A, Tjwa M, Moons L, Wu Y, Angelillo-Scherrer A, Liao F, et al. Revascularization of ischemic tissues by PlGF treatment, and inhibition of tumor angiogenesis, arthritis and atherosclerosis by anti-Flt1. Nat Med. (2002) 8:831–40. doi: 10.1038/nm731

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Salcedo R, Ponce ML, Young HA, Wasserman K, Ward JM, Kleinman HK, et al. Human endothelial cells express CCR2 and respond to MCP-1: direct role of MCP-1 in angiogenesis and tumor progression. Blood. (2000) 96:34–40. doi: 10.1182/blood.V96.1.34.013a49_34_40

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Wu Y, Ip JE, Huang J, Zhang L, Matsushita K, Liew CC, et al. Essential role of ICAM-1/CD18 in mediating EPC recruitment, angiogenesis, and repair to the infarcted myocardium. Circ Res. (2006) 99:315–22. doi: 10.1161/01.RES.0000235986.35957.a3

PubMed Abstract | CrossRef Full Text | Google Scholar