Transcriptome Analysis of Hypertrophic Heart Tissues from Murine Transverse Aortic Constriction and Human Aortic Stenosis Reveals Key Genes and Transcription Factors Involved in Cardiac Remodeling Induced by Mechanical Stress
1. Introduction
Heart failure, the end stage for most cardiac diseases, is a clinical syndrome in which the heart is unable to provide sufficient blood flow to meet physiologic requirements of the body. Prior to clinical symptoms or signs of heart failure, the body tries to maintain adequate tissue perfusion using several mechanisms, including the Frank–Starling mechanism and neurohormonal activation, which lead to cardiac remodeling [1].
Cardiac remodeling is a process in which genomic changes occur. Physiologically, signaling and transcriptional control involve precise programs of gene activation and suppression [2]. Transcriptional changes in response to pathological stress might promote deterioration of cardiac remodeling. It has been shown that preventing the genomic changes may be a promising therapeutic approach [2–4].
The murine transverse aortic constriction (TAC) is a commonly used experimental model for mechanical stress-induced cardiac remodeling, which clinically mimics the aortic stenosis (AS). TAC initially leads to compensated hypertrophy of the heart and is often associated with a temporary enhancement of cardiac contractility. In the end stage, the response to chronic hemodynamic overload becomes maladaptive, leading to cardiac dilatation and heart failure. The murine TAC model has since been extensively used as a valuable tool to mimic human cardiovascular diseases and elucidate fundamental signaling processes involved in the cardiac hypertrophic response and heart failure development. It provides a more reproducible model of cardiac hypertrophy and a more gradual time course for the development of heart failure [5].
DNA microarrays facilitate measurement of the expression levels of large numbers of genes simultaneously. Recent data underscored the significance of genomic mechanisms in regulating gene expression programs in cardiac pathology [2]. A number of studies have investigated the genomic changes of the heart in the process of cardiac remodeling.
To investigate the genomic changes in the process of cardiac remodeling induced by mechanical stress systematically and without bias, we collected and analyzed DNA microarray data for cardiac remodeling induced by TAC and AS from the Gene Expression Omnibus (GEO) database and European Bioinformatics Institute (EBI). As a result, we found a set of gene expression changes in the cardiac pathologic remodeling induced by mechanical stress that shared some common transcription factors (TFs) with each other.
2. Methods
2.1. Microarray Data Collection and Preprocessing
The gene expression profiles were screened and downloaded from the National Center for Biotechnology Information GEO database and the EMBL-EBI. To explore cardiac remodeling under mechanical stress, the murine TAC datasets and the human AS datasets were included. The TAC datasets in which hypertrophic genes NPPA, NPPB, ACTA1, and MYH7/MYH6 remained unchanged were excluded from analysis. Datasets with the number of samples in each group of <3 were also excluded.
2.2. DEG Analysis
GEO series were analyzed separately using the online GEO2R tool with default parameters (https://www.ncbi.nlm.nih.gov/geo/geo2r/), in which the empirical Bayes algorithm (function “eBayes”) in the limma package was used to detect differentially expressed genes between the TAC model or AS patients and controls. In the murine model analysis, the genes with a P value (Bayes test) of <0.05 were considered as DEGs for the multiple intersection of different datasets.
Since datasets were from different research centers, group variation was present. It was not possible to conduct the data analysis on interdatasets. Considering these limitations, we obtained only the average values of logFC from each dataset to represent the expression levels [6]. In the analysis of AS patients, significantly changed genes were defined by a logarithmic-transformed fold-change absolute value (log2(FC)) ≥ 1 and a P value of ≤0.05.
2.3. Venn Analysis
Comparative analysis was carried out with the InteractiVenn tool (http://www.interactivenn.net/) [7] and Bioinformatics and Evolutionary Genomics tool (http://bioinformatics.psb.ugent.be/webtools/Venn/).
2.4. GO Analysis
DAVID was employed to perform the GO analysis for biological processes and pathway enrichment. To plot the BPs of the DEGs involved, we used the clusterProfiler package [8].
2.5. PPI Network Construction Analysis
STRING online tool (string-db.org) [9] was used to establish a PPI network for the murine TAC model. Cytoscape software [10] was used to establish a PPI network for DEGs of AS patients, with the cutoff of a combined score of >0.4. The network analyzer plug-in for the Cytoscape software was used to analyze the topological property of the networks [6]. Genes with the edge degree of ≥7 were defined as hub genes in this article.
2.6. TF Analysis
The Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining version 2 database (https://www.grnpedia.org/trrust/) [11] was used to predict regulation of TFs based on the lists of upregulated and downregulated genes generated across the microarray datasets. Significant TFs and potentially regulated genes were identified based on a multiple parameters, P < 0.05 [12]. We used the “igraph” package in R to visualize the output results.
2.1. Microarray Data Collection and Preprocessing
The gene expression profiles were screened and downloaded from the National Center for Biotechnology Information GEO database and the EMBL-EBI. To explore cardiac remodeling under mechanical stress, the murine TAC datasets and the human AS datasets were included. The TAC datasets in which hypertrophic genes NPPA, NPPB, ACTA1, and MYH7/MYH6 remained unchanged were excluded from analysis. Datasets with the number of samples in each group of <3 were also excluded.
2.2. DEG Analysis
GEO series were analyzed separately using the online GEO2R tool with default parameters (https://www.ncbi.nlm.nih.gov/geo/geo2r/), in which the empirical Bayes algorithm (function “eBayes”) in the limma package was used to detect differentially expressed genes between the TAC model or AS patients and controls. In the murine model analysis, the genes with a P value (Bayes test) of <0.05 were considered as DEGs for the multiple intersection of different datasets.
Since datasets were from different research centers, group variation was present. It was not possible to conduct the data analysis on interdatasets. Considering these limitations, we obtained only the average values of logFC from each dataset to represent the expression levels [6]. In the analysis of AS patients, significantly changed genes were defined by a logarithmic-transformed fold-change absolute value (log2(FC)) ≥ 1 and a P value of ≤0.05.
2.3. Venn Analysis
Comparative analysis was carried out with the InteractiVenn tool (http://www.interactivenn.net/) [7] and Bioinformatics and Evolutionary Genomics tool (http://bioinformatics.psb.ugent.be/webtools/Venn/).
2.4. GO Analysis
DAVID was employed to perform the GO analysis for biological processes and pathway enrichment. To plot the BPs of the DEGs involved, we used the clusterProfiler package [8].
2.5. PPI Network Construction Analysis
STRING online tool (string-db.org) [9] was used to establish a PPI network for the murine TAC model. Cytoscape software [10] was used to establish a PPI network for DEGs of AS patients, with the cutoff of a combined score of >0.4. The network analyzer plug-in for the Cytoscape software was used to analyze the topological property of the networks [6]. Genes with the edge degree of ≥7 were defined as hub genes in this article.
2.6. TF Analysis
The Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining version 2 database (https://www.grnpedia.org/trrust/) [11] was used to predict regulation of TFs based on the lists of upregulated and downregulated genes generated across the microarray datasets. Significant TFs and potentially regulated genes were identified based on a multiple parameters, P < 0.05 [12]. We used the “igraph” package in R to visualize the output results.
3. Results
3.1. Datasets Involved in This Study
We searched a total of 14 datasets, which included a model of murine cardiac remodeling induced by TAC and utilized a microarray to detect differential expression genes (DEGs) in an unbiased manner. Four datasets for the early cardiac remodeling and seven datasets for the late cardiac remodeling were used (Table 1).
Table 1
3.2. Genomic Changes in the Early Hypertrophic Response Stage
The period within two weeks after the TAC operation was defined as the early stage of cardiac remodeling, characterized by compensated hypertrophic remodeling. Four datasets included the microarray data from analysis of this period.
The four datasets shared 251 significant DEGs (Figure 1(a)), among which only 74 exhibited similar trends. The heatmap showed DEGs with the same trends across the four datasets (Figure 1(b)). The 74 DEGs were analyzed using the STRING online tools (Figure 1(c)). To show the main biological processes involving DEGs, we performed Gene Ontology (GO) analysis using the Database for Annotation, Visualization and Integrated Discovery (DAVID) (), with the results represented in Figure 1(d).
Intersections among the four datasets comprised only a small portion of each dataset. However, among the intersections, the number of DEGs with same trends was even smaller. STRING analysis showed that the 74 DEGs were mainly concentrated in metabolic changes, cell differentiation and growth, cell cycling, and collagen fibril organization (Figure 1(c)). The BP enrichment of the DEGs mainly occurred during the lipid metabolic change ().
3.3. Genomic Changes in the Late Stage of Cardiac Remodeling
We then analyzed the datasets detecting DEGs more than four weeks post-TAC, which represented gene changes in the late stage of cardiac remodeling. The analysis involved seven datasets. The Venn analysis of DEGs is shown in . In the seven datasets, only 16 DEGs exhibited the same trends, which is shown as a heatmap in Figure 2(a). We consequently performed the protein–protein interaction (PPI) analysis of the 16 DEGs using the STRING tools (Figure 2(b)). Biological processes involving these 16 DEGs were concentrated mainly in collagen biosynthesis and hypertrophic marker molecules, such as NPPA, NPPB, and ACTA1 (Table 2).
Table 2
Term | Count | % | P value | Genes |
---|---|---|---|---|
GO:0071260~cellular response to mechanical stimulus | 4 | 25 | 2.72E − 05 | NPPB, COL1A1, ANKRD1, NPPA |
GO:0071560~cellular response to transforming growth factor beta stimulus | 3 | 18.75 | 9.92E − 04 | POSTN, COL1A1, ANKRD1 |
GO:0035582~sequestering of BMP in extracellular matrix | 2 | 12.5 | 3.09E − 03 | NBL1, FBN1 |
GO:0071356~cellular response to tumor necrosis factor | 3 | 18.75 | 3.18E − 03 | POSTN, COL1A1, ANKRD1 |
GO:0030308~negative regulation of cell growth | 3 | 18.75 | 4.09E − 03 | NPPB, FLCN, NPPA |
GO:0007168~receptor guanylyl cyclase signaling pathway | 2 | 12.5 | 6.95E − 03 | NPPB, NPPA |
GO:0061049~cell growth involved in cardiac muscle cell development | 2 | 12.5 | 8.49E − 03 | NPPB, NPPA |
GO:0001666~response to hypoxia | 3 | 18.75 | 9.39E − 03 | NPPB, POSTN, NPPA |
GO:0003085~negative regulation of systemic arterial blood pressure | 2 | 12.5 | 1.00E − 02 | NPPB, NPPA |
GO:0006182~cGMP biosynthetic process | 2 | 12.5 | 1.31E − 02 | NPPB, NPPA |
Although intersection of the seven datasets credits the genes involved in the TAC-induced cardiac remodeling, its comprehensiveness may be attenuated for the multiple intersections. Additionally, we performed PPI and BP analyses for the DEGs in the intersection of at least six sets ( and ). The results showed collagen biosynthesis process, innate immune response, metabolic changes, and ion transmembrane transport to be the main changes involved in the late remodeling stage.
3.4. Microarray Data Analysis of the Human Heart Tissue from AS Patients
TAC is a common model used to investigate cardiac remodeling and heart failure. Clinically, heart failure is a syndrome with multiple heterogeneous etiologies. Hypertension and AS are the main heart failure types induced by mechanical stress, a model of TAC.
To investigate the DEGs involved in human heart failure induced by mechanical stress, we analyzed {"type":"entrez-geo","attrs":{"text":"GSE1145","term_id":"1145"}}GSE1145, in which datasets for the heart tissues from AS patients were utilized. The total DEG count was 252. Some of these genes are represented by a heatmap in Figure 3(a). BP analysis showed that the genes mainly enriched the inflammatory process, in addition to playing a role in muscular hypertrophic changes (Figure 3(c)). In the PPI analysis, four genes were identified as hub genes with the edge degree ≥ 7. According to the edge degree rank, the four hub genes were IL-8, JAK2, AGTR1, and BCR. IL-8, in particular, might play an important role in the development of mechanical stress induced by AS. However, these four genes were not involved in the analysis of ischemic cardiomyopathy [6], implying a distinct pathogenesis between these two cardiomyopathies.
Compared with the mechanical stress-induced cardiac remodeling in mice, few DEGs or BPs overlapped between the murine TAC model and AS patients, in which the effects of clinical medication had to be excluded.
3.5. TF Analysis
We also predicted the TFs regulating DEGs using the data from the TAC model and AS patients. Although the DEGs shared little overlap between human and murine mechanical stress-induced hypertrophic heart tissue, there were four TFs (SP1, CEBPB, PPARG, and NFKB1) in common between early cardiac remodeling (Figure 4(a)) and AS patients (Figure 4(c)). The TFs predicted in the late cardiac remodeling were few for a little set of DEGs (Figure 4(b)).
The most prominent TF was SP1, which is involved in many cellular processes, including cell differentiation, cell growth, apoptosis, immune response, response to DNA damage, and chromatin remodeling. Activity of CEPPB and NFKB1 is important in the regulation of genes involved in immune and inflammatory responses. PPARG is a regulator of metabolic changes.
3.1. Datasets Involved in This Study
We searched a total of 14 datasets, which included a model of murine cardiac remodeling induced by TAC and utilized a microarray to detect differential expression genes (DEGs) in an unbiased manner. Four datasets for the early cardiac remodeling and seven datasets for the late cardiac remodeling were used (Table 1).
Table 1
3.2. Genomic Changes in the Early Hypertrophic Response Stage
The period within two weeks after the TAC operation was defined as the early stage of cardiac remodeling, characterized by compensated hypertrophic remodeling. Four datasets included the microarray data from analysis of this period.
The four datasets shared 251 significant DEGs (Figure 1(a)), among which only 74 exhibited similar trends. The heatmap showed DEGs with the same trends across the four datasets (Figure 1(b)). The 74 DEGs were analyzed using the STRING online tools (Figure 1(c)). To show the main biological processes involving DEGs, we performed Gene Ontology (GO) analysis using the Database for Annotation, Visualization and Integrated Discovery (DAVID) (), with the results represented in Figure 1(d).
Intersections among the four datasets comprised only a small portion of each dataset. However, among the intersections, the number of DEGs with same trends was even smaller. STRING analysis showed that the 74 DEGs were mainly concentrated in metabolic changes, cell differentiation and growth, cell cycling, and collagen fibril organization (Figure 1(c)). The BP enrichment of the DEGs mainly occurred during the lipid metabolic change ().
3.3. Genomic Changes in the Late Stage of Cardiac Remodeling
We then analyzed the datasets detecting DEGs more than four weeks post-TAC, which represented gene changes in the late stage of cardiac remodeling. The analysis involved seven datasets. The Venn analysis of DEGs is shown in . In the seven datasets, only 16 DEGs exhibited the same trends, which is shown as a heatmap in Figure 2(a). We consequently performed the protein–protein interaction (PPI) analysis of the 16 DEGs using the STRING tools (Figure 2(b)). Biological processes involving these 16 DEGs were concentrated mainly in collagen biosynthesis and hypertrophic marker molecules, such as NPPA, NPPB, and ACTA1 (Table 2).
Table 2
Term | Count | % | P value | Genes |
---|---|---|---|---|
GO:0071260~cellular response to mechanical stimulus | 4 | 25 | 2.72E − 05 | NPPB, COL1A1, ANKRD1, NPPA |
GO:0071560~cellular response to transforming growth factor beta stimulus | 3 | 18.75 | 9.92E − 04 | POSTN, COL1A1, ANKRD1 |
GO:0035582~sequestering of BMP in extracellular matrix | 2 | 12.5 | 3.09E − 03 | NBL1, FBN1 |
GO:0071356~cellular response to tumor necrosis factor | 3 | 18.75 | 3.18E − 03 | POSTN, COL1A1, ANKRD1 |
GO:0030308~negative regulation of cell growth | 3 | 18.75 | 4.09E − 03 | NPPB, FLCN, NPPA |
GO:0007168~receptor guanylyl cyclase signaling pathway | 2 | 12.5 | 6.95E − 03 | NPPB, NPPA |
GO:0061049~cell growth involved in cardiac muscle cell development | 2 | 12.5 | 8.49E − 03 | NPPB, NPPA |
GO:0001666~response to hypoxia | 3 | 18.75 | 9.39E − 03 | NPPB, POSTN, NPPA |
GO:0003085~negative regulation of systemic arterial blood pressure | 2 | 12.5 | 1.00E − 02 | NPPB, NPPA |
GO:0006182~cGMP biosynthetic process | 2 | 12.5 | 1.31E − 02 | NPPB, NPPA |
Although intersection of the seven datasets credits the genes involved in the TAC-induced cardiac remodeling, its comprehensiveness may be attenuated for the multiple intersections. Additionally, we performed PPI and BP analyses for the DEGs in the intersection of at least six sets ( and ). The results showed collagen biosynthesis process, innate immune response, metabolic changes, and ion transmembrane transport to be the main changes involved in the late remodeling stage.
3.4. Microarray Data Analysis of the Human Heart Tissue from AS Patients
TAC is a common model used to investigate cardiac remodeling and heart failure. Clinically, heart failure is a syndrome with multiple heterogeneous etiologies. Hypertension and AS are the main heart failure types induced by mechanical stress, a model of TAC.
To investigate the DEGs involved in human heart failure induced by mechanical stress, we analyzed {"type":"entrez-geo","attrs":{"text":"GSE1145","term_id":"1145"}}GSE1145, in which datasets for the heart tissues from AS patients were utilized. The total DEG count was 252. Some of these genes are represented by a heatmap in Figure 3(a). BP analysis showed that the genes mainly enriched the inflammatory process, in addition to playing a role in muscular hypertrophic changes (Figure 3(c)). In the PPI analysis, four genes were identified as hub genes with the edge degree ≥ 7. According to the edge degree rank, the four hub genes were IL-8, JAK2, AGTR1, and BCR. IL-8, in particular, might play an important role in the development of mechanical stress induced by AS. However, these four genes were not involved in the analysis of ischemic cardiomyopathy [6], implying a distinct pathogenesis between these two cardiomyopathies.
Compared with the mechanical stress-induced cardiac remodeling in mice, few DEGs or BPs overlapped between the murine TAC model and AS patients, in which the effects of clinical medication had to be excluded.
3.5. TF Analysis
We also predicted the TFs regulating DEGs using the data from the TAC model and AS patients. Although the DEGs shared little overlap between human and murine mechanical stress-induced hypertrophic heart tissue, there were four TFs (SP1, CEBPB, PPARG, and NFKB1) in common between early cardiac remodeling (Figure 4(a)) and AS patients (Figure 4(c)). The TFs predicted in the late cardiac remodeling were few for a little set of DEGs (Figure 4(b)).
The most prominent TF was SP1, which is involved in many cellular processes, including cell differentiation, cell growth, apoptosis, immune response, response to DNA damage, and chromatin remodeling. Activity of CEPPB and NFKB1 is important in the regulation of genes involved in immune and inflammatory responses. PPARG is a regulator of metabolic changes.
4. Discussion
Using the data from the high throughput DNA microarray analysis, we were able to systemically reveal genomic changes in a disease, so that potential therapeutic targets could be identified in the future.
In this study, we investigated the datasets for TAC, a typical model to explore cardiac remodeling. We divided the datasets into early and late phases of cardiac remodeling, according to data from mice that succumbed days after the TAC operation. After analysis of the data, we found common gene changes within different datasets, which mainly converged on matrix remodeling, metabolic changes, and mechanical response.
Genomic changes in cardiac remodeling have recently gained attention from researchers and their modulation has been widely investigated. The methylation of DNA [23] and histones [24], acetylation of chromatin and facilitation of transcriptional activation [3, 25], and chromatin structural remodeling [26] all result in genomic changes and lead to heart failure. Suppression of genomic changes could ameliorate cardiac remodeling. Thus, it is important to determine the genomic changes taking place during heart failure. Our study represents the first attempt to systematically elucidate these changes.
In a murine model, the genome is altered in the early stages of metabolic changes, cell differentiation and growth, cell cycling, and collagen fibril organization. A recent study has revealed that cyclins and TGF-β that have terminally exited the cell cycle can unlock the proliferative potential in the myocardium. Moreover, their overexpression could improve the cardiac function [27]. The PPI analysis showed that SLC2A4 and TOP2A are the two centers of genomic change. SLC2A4, also known as GLUT4, is a glucose transporter that facilitates the metabolic switch to glucose in cardiac remodeling. TOP2A (DNA topoisomerase II-alpha) controlled the topological states of DNA by transient breakage and subsequent rejoining of DNA strands that facilitated cellular mitosis, chromatin remodeling, and gene transcription [28–30].
There were fewer gene changes in the late stages of cardiac remodeling, where only 16 genes exhibited similar trends in all of the datasets. The upregulated genes ACTA1, NPPA, NPPB, POSTN, COLIA1, and COL8A1 were regarded as molecular markers in the pathologic process of cardiac remodeling. Ces1d was downregulated in all datasets. It has been shown to be involved in lipolysis, the process whereby the adipocyte hydrolyzes stored triglycerides into fatty acids to be used as fuel in times of need [31, 32], in correlation with the opinion that the switch from fat to glucose is an approach that could be taken to improve cardiac remodeling [33]. FLCN, the inactivation of which could potentially lead to cardiac remodeling [34], was also downregulated. Research reports involving other genes from the set of 16 identified in this study, including P3H4, ANKRD1, CPXM2, FBN1, FXYD5, MFAP5, NBL1, PFKP, and SLMAP, were rare. These genes are therefore worth exploring further.
To delineate the correlation between the murine model and clinical patients, we analyzed datasets from AS patients mimicked by TAC [5]. The results showed that the genomic changes in AS biological processes were mainly in inflammation. Results from the PPI analysis identified IL-8, JAK2, AGTR1, and BCR to be the centers of genomic changes, in which AGTR1 was the target of hypertension, the common cause of cardiac remodeling triggered by mechanical stress. Sartans, antagonists of AGTR1, are the cornerstone of medication for hypertension. Accordingly, IL-8, JAK2, and BCR might be therapeutic targets for hypertension, which were not detected in the microarray data from the TAC model.
The DEG intersections between the murine TAC model and AS patients were few. However, predicted TFs from DEGs SP1, CEBPB, PPARG, and NFKB1 were the common TFs between the two species. Unsurprisingly, they were either regulators in metabolic changes or pivotal hubs in inflammatory response. SP1 has been reported to contribute to the regulation of critical molecules involved in cardiac remodeling [35, 36]. The mice downregulated of CEBPB has been reported to display substantial resistance to cardiac failure upon pressure overload, indicating its repression of cardiomyocyte growth and proliferation in the adult mammalian heart [37]. The PPAR gene pathway coordinately act to regulate cellular processes central to glucose and lipid metabolism [38]. NFKB1 signaling also is critical for both cardiac remodeling and hypertrophy [39].
5. Conclusion
In conclusion, we offer a novel and comprehensive analysis of gene expression profiles using microarray DNA datasets in cardiac remodeling induced by mechanical stress, mimicking hypertension. Genes involved in metabolic changes, extracellular matrix remodeling, and cell differentiation and growth were significantly changed in the heart tissue from the murine TAC model. Compared to results from the TAC model, the significantly changed genes in patients suffering from AS were mostly enriched during the inflammatory biological processes. The analysis will provide valuable information for future research on the molecular mechanisms of cardiac remodeling and offer clues for the discovery of novel therapeutic strategies.
6. Limitations
Although our analysis was comprehensive, with high throughput and a large sample size, some limitations were still present.
Classification of the cardiac remodeling stages in the murine TAC model was performed on the basis of time passed postoperation. No echo or histological standard was used. The TAC model was performed using a standard operating procedure [5]. The period within two weeks after the operation was regarded as hypertrophic or compensatory stage, while the period past four weeks was considered to be the dilated or decompensated stage.
Furthermore, the data for patients suffering from AS were extracted from {"type":"entrez-geo","attrs":{"text":"GSE1145","term_id":"1145"}}GSE1145, which lacks detailed clinical information as no research was published using this dataset. However, the DEGs mainly involved in the inflammatory biological processes were similar to a previous study of cardiac remodeling [6, 40].
Despite these limitations, the comprehensive analysis of microarray data in this study makes the results compelling.
Abstract
Background
Mechanical stress-induced cardiac remodeling that results in heart failure is characterized by transcriptional reprogramming of gene expression. However, a systematic study of genomic changes involved in this process has not been performed to date. To investigate the genomic changes and underlying mechanism of cardiac remodeling, we collected and analyzed DNA microarray data for murine transverse aortic constriction (TAC) and human aortic stenosis (AS) from the Gene Expression Omnibus database and the European Bioinformatics Institute.
Methods and Results
The differential expression genes (DEGs) across the datasets were merged. The Venn diagrams showed that the number of intersections for early and late cardiac remodeling was 74 and 16, respectively. Gene ontology and protein–protein interaction network analysis showed that metabolic changes, cell differentiation and growth, cell cycling, and collagen fibril organization accounted for a great portion of the DEGs in the TAC model, while in AS patients' immune system signaling and cytokine signaling displayed the most significant changes. The intersections between the TAC model and AS patients were few. Nevertheless, the DEGs of the two species shared some common regulatory transcription factors (TFs), including SP1, CEBPB, PPARG, and NFKB1, when the heart was challenged by applied mechanical stress.
Conclusions
This study unravels the complex transcriptome profiles of the heart tissues and highlighting the candidate genes involved in cardiac remodeling induced by mechanical stress may usher in a new era of precision diagnostics and treatment in patients with cardiac remodeling.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (grant numbers 81670355, 81600294, and 81521001).
Abbreviations
TAC: | Transverse aortic constriction |
AS: | Aortic stenosis |
DEG: | Differential expression gene |
PPI: | Protein–protein interaction |
GO: | Gene Ontology |
BP: | Biological process |
DAVID: | The Database for Annotation, Visualization, and Integrated Discovery. |
References
- 1. Kemp C. D., Conte J. V. The pathophysiology of heart failure. Cardiovascular Pathology. 2012;21(5):365–371. doi: 10.1016/j.carpath.2011.11.007.] [[PubMed]
- 2. Mathiyalagan P., Keating S. T., Du X. J., El-Osta A. Chromatin modifications remodel cardiac gene expression. Cardiovascular Research. 2014;103(1):7–16. doi: 10.1093/cvr/cvu122.] [[PubMed]
- 3. Anand P., Brown J. D., Lin C. Y., et al. BET bromodomains mediate transcriptional pause release in heart failure. Cell. 2013;154(3):569–582. doi: 10.1016/j.cell.2013.07.013.] [
- 4. Hannenhalli STranscriptional genomics associates FOX transcription factors with human heart failure. Circulation. 2006;114(12):1269–1276. doi: 10.1161/CIRCULATIONAHA.106.632430.] [[PubMed]
- 5. De Almeida A. C., van Oort R. J., Wehrens X. H. T. Transverse aortic constriction in mice. Journal of Visualized Experiments. 2010;38(38, article e1729) doi: 10.3791/1729.] [
- 6. Li Y., Jiang Q., Ding Z., et al Identification of a common different gene expression signature in ischemic cardiomyopathy. Genes. 2018;9(1):p. 56. doi: 10.3390/genes9010056.] [
- 7. Heberle H., Meirelles G. V., Da Silva F. R., Telles G. P., Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. 2015;16(1):p. 169. doi: 10.1186/s12859-015-0611-3.] [
- 8. Yu G., Wang L. G., Yan G. R., He Q. Y. DOSE: an R/Bioconductor package for disease ontology semantic and enrichment analysis. Bioinformatics. 2015;31(4):608–609. doi: 10.1093/bioinformatics/btu684.] [[PubMed]
- 9. Szklarczyk D., Morris J. H., Cook H., et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Research. 2017;45(D1):D362–D368. doi: 10.1093/nar/gkw937.] [
- 10. Doncheva N. T., Assenov Y., Domingues F. S., Albrecht M. Topological analysis and interactive visualization of biological networks and protein structures. Nature Protocols. 2012;7(4):670–685. doi: 10.1038/nprot.2012.004.] [[PubMed]
- 11. Han H., Cho J., Lee S., et al TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Research. 2018;46(D1):D380–D386. doi: 10.1093/nar/gkx1013.] [
- 12. Coleman M. A., Sasi S. P., Onufrak J., et al. Low-dose radiation affects cardiac physiology: gene networks and molecular signaling in cardiomyocytes. American Journal of Physiology-Heart and Circulatory Physiology. 2015;309(11):H1947–H1963. doi: 10.1152/ajpheart.00050.2015.] [
- 13. Liao X., Zhang R., Lu Y., et al Kruppel-like factor 4 is critical for transcriptional control of cardiac mitochondrial homeostasis. The Journal of Clinical Investigation. 2015;125(9):3461–3476. doi: 10.1172/JCI79964.] [
- 14. Zhao M., Chow A., Powers J., Fajardo G., Bernstein DMicroarray analysis of gene expression after transverse aortic constriction in mice. Physiological Genomics. 2004;19(1):93–105. doi: 10.1152/physiolgenomics.00040.2004.] [[PubMed]
- 15. Bisping E., Ikeda S., Kong S. W., et al. Gata4 is required for maintenance of postnatal cardiac function and protection from pressure overload-induced heart failure. Proceedings of the National Academy of Sciences. 2006;103(39):14471–14476. doi: 10.1073/pnas.0602543103.] [
- 16. Mirotsou M., Watanabe C. M. H., Schultz P. G., Pratt R. E., Dzau V. J. Elucidating the molecular mechanism of cardiac remodeling using a comparative genomic approach. Physiological Genomics. 2003;15(2):115–126. doi: 10.1152/physiolgenomics.00071.2003.] [[PubMed]
- 17. Colston J. T., Boylston W. H., Feldman M. D., et al. Interleukin-18 knockout mice display maladaptive cardiac hypertrophy in response to pressure overload. Biochemical and Biophysical Research Communications. 2007;354(2):552–558. doi: 10.1016/j.bbrc.2007.01.030.] [
- 18. Smeets P. J. H., de Vogel-van D. B. H., Willemsen P. H. M., et al. Transcriptomic analysis of PPARalpha-dependent alterations during cardiac hypertrophy. Physiological Genomics. 2008;36(1):15–23. doi: 10.1152/physiolgenomics.90296.2008.] [[PubMed]
- 19. Mirotsou M., Dzau V. J., Pratt R. E., Weinberg E. O. Physiological genomics of cardiac disease: quantitative relationships between gene expression and left ventricular hypertrophy. Physiological Genomics. 2006;27(1):86–94. doi: 10.1152/physiolgenomics.00028.2006.] [[PubMed]
- 20. Mayer S. C., Gilsbach R., Preissl S., et al. Adrenergic repression of the epigenetic reader MeCP2 facilitates cardiac adaptation in chronic heart failure. Circulation Research. 2015;117(7):622–633. doi: 10.1161/CIRCRESAHA.115.306721.] [
- 21. Nuhrenberg T. G., Hammann N., Schnick T., et al. Cardiac myocyte de novo DNA methyltransferases 3a/3b are dispensable for cardiac function and remodeling after chronic pressure overload in mice. PLoS One. 2015;10(6, article e131019) doi: 10.1371/journal.pone.0131019.] [
- 22. Lai L., Leone T. C., Keller M. P., et al. Energy metabolic reprogramming in the hypertrophied and early stage failing heart: a multisystems approach. Circulation Heart Failure. 2014;7(6):1022–1031. doi: 10.1161/CIRCHEARTFAILURE.114.001469.] [
- 23. Greco C. M., Kunderfranco P., Rubino M., et al. DNA hydroxymethylation controls cardiomyocyte gene expression in development and hypertrophy. Nature Communications. 2016;7(1):p. 12418. doi: 10.1038/ncomms12418.] [
- 24. Papait R., Serio S., Pagiatakis C., et al Histone methyltransferase G9a is required for cardiomyocyte homeostasis and hypertrophy. Circulation. 2017;136(13):1233–1246. doi: 10.1161/circulationaha.117.028561.] [[PubMed]
- 25. Haldar S. M., McKinsey T. A. BET-ting on chromatin-based therapeutics for heart failure. Journal of Molecular and Cellular Cardiology. 2014;74:98–102. doi: 10.1016/j.yjmcc.2014.05.002.] [
- 26. Rosa-Garrido M., Chapski D. J., Schmitt A. D., et al. High-resolution mapping of chromatin conformation in cardiac myocytes reveals structural remodeling of the epigenome in heart failure. Circulation. 2017;136(17):1613–1625. doi: 10.1161/CIRCULATIONAHA.117.029430.] [
- 27. Mohamed T. M. A., Ang Y., Radzinsky E., et al. Regulation of cell cycle to stimulate adult cardiomyocyte proliferation and cardiac regeneration. Cell. 2018;173(1):104–116.e12. doi: 10.1016/j.cell.2018.02.014.] [
- 28. King I. F., Yandava C. N., Mabb A. M., et al. Topoisomerases facilitate transcription of long genes linked to autism. Nature. 2013;501(7465):58–62. doi: 10.1038/nature12504.] [
- 29. Dykhuizen E. C., Hargreaves D. C., Miller E. L., et al. BAF complexes facilitate decatenation of DNA by topoisomerase IIαNature. 2013;497(7451):624–627. doi: 10.1038/nature12146.] [
- 30. Dong K. C., Berger J. M. Structural basis for gate-DNA recognition and bending by type IIA topoisomerases. Nature. 2007;450(7173):1201–1205. doi: 10.1038/nature06396.] [[PubMed]
- 31. Dominguez E., Galmozzi A., Chang J. W., et al. Integrated phenotypic and activity-based profiling links Ces3 to obesity and diabetes. Nature Chemical Biology. 2014;10(2):113–121. doi: 10.1038/nchembio.1429.] [
- 32. Wei E., Gao W., Lehner RAttenuation of adipocyte triacylglycerol hydrolase activity decreases basal fatty acid efflux. The Journal of Biological Chemistry. 2007;282(11):8027–8035. doi: 10.1074/jbc.M605789200.] [[PubMed]
- 33. Zhang L., Lu Y., Jiang H., et al Additional use of trimetazidine in patients with chronic heart failure: a meta-analysis. Journal of the American College of Cardiology. 2012;59(10):913–922. doi: 10.1016/j.jacc.2011.11.027.] [[PubMed]
- 34. Hasumi Y., Baba M., Hasumi H., et al Folliculin (Flcn) inactivation leads to murine cardiac hypertrophy through mTORC1 deregulation. Human Molecular Genetics. 2014;23(21):5706–5719. doi: 10.1093/hmg/ddu286.] [
- 35. Li T., Chen Y. H., Liu T. J., et al. Using DNA microarray to identify Sp1 as a transcriptional regulatory element of insulin-like growth factor 1 in cardiac muscle cells. Circulation Research. 2003;93(12):1202–1209. doi: 10.1161/01.RES.0000104085.76261.02.] [[PubMed]
- 36. Reddy V. S., Prabhu S. D., Mummidi S., et al. Interleukin-18 induces EMMPRIN expression in primary cardiomyocytes via JNK/Sp1 signaling and MMP-9 in part via EMMPRIN and through AP-1 and NF-kappaB activation. American Journal of Physiology Heart and Circulatory Physiology. 2010;299(4):H1242–H1254. doi: 10.1152/ajpheart.00451.2010.] [
- 37. Bostrom P., Mann N., Wu J., et al C/EBPβ Controls Exercise-Induced Cardiac Growth and Protects against Pathological Cardiac Remodeling. Cell. 2010;143(7):1072–1083. doi: 10.1016/j.cell.2010.11.036.] [
- 38. Mistry N. F., Cresci S. PPAR transcriptional activator complex polymorphisms and the promise of individualized therapy for heart failure. Heart Failure Reviews. 2010;15(3, article 9114):197–207. doi: 10.1007/s10741-008-9114-x.] [
- 39. Gaspar-Pereira S., Fullard N., Townsend P. A., et al. The NF-κB subunit c-Rel stimulates cardiac hypertrophy and fibrosis. The American Journal of Pathology. 2012;180(3):929–939. doi: 10.1016/j.ajpath.2011.11.007.] [
- 40. Zhang L., Liu M., Jiang H., et al Extracellular high‐mobility group box 1 mediates pressure overload‐induced cardiac hypertrophy and heart failure. Journal of Cellular and Molecular Medicine. 2016;20(3):459–470. doi: 10.1111/jcmm.12743.] [