Immune-Related Genes: Potential Regulators and Drug Therapeutic Targets in Hypertrophic Cardiomyopathy

Liu, Wei; Ye, Ju; Cai, Jinqiang; Xie, Feng; Tang, Mengjie; Yang, QinJing

doi:https://doi.org/10.1155/2021/5528347

Journal of Nanomaterials

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Novel Micro- and Nanomaterials for Pharmaceutical Applications

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 5528347 | https://doi.org/10.1155/2021/5528347

Immune-Related Genes: Potential Regulators and Drug Therapeutic Targets in Hypertrophic Cardiomyopathy

Wei Liu,¹Ju Ye,²Jinqiang Cai,³Feng Xie,⁴Mengjie Tang,⁵and QinJing Yang¹

Academic Editor: Tingting Hong

Received29 Jan 2021

Revised05 Feb 2021

Accepted09 Feb 2021

Published03 Mar 2021

Abstract

Background. Accumulating evidence shows that the innate immune system is a key player in cardiovascular repair and regeneration, but little is known about the role of immune-related genes (IRGs) in hypertrophic cardiomyopathy (HCM). Methods. The differential mRNA expression profiles of HCM samples were downloaded from the Gene Expression Omnibus (GEO) dataset (GSE89714), and the IRG expression profile was obtained from the ImmPort database. The regulatory pathways of IRGs in HCM were screened out through discrepantly expressive genes (DEGs) analysis, enrichment of gene function/pathway analysis, and protein-protein interaction (PPI) network. Besides, hub IRGs in the PPI network were selected for drug prediction. Results. A total of 854 genes were differentially expressed in HCM, of which 88 were IRGs. Functional enrichment analysis revealed that 88 IRGs were mainly involved in the biological processes (BP) of SMAD protein pathway, smooth muscle cell proliferation, protein serine/threonine kinase, and mitogen-activated protein kinase (MAPK) cascade. Cytokine-cytokine receptor interaction, TGFβ signaling pathway, PI3K-Akt signaling pathway, and MAPK signaling pathway were enriched in the pathway enrichment analysis of these 88 IRGs. Furthermore, the PPI regulatory network of IRGs was constructed, and 10 hub IRGs were screened out to construct a regulatory network for HCM. 4 transcription factors (TFs) were the major regulator of 10 hub IRGs. Finally, these 10 hub IRGs were entered into the pharmacogenomics database for prediction, and the relevant drugs were obtained. Conclusions. In this study, 10 hub IRGs were coexpressed with 4 TFs to construct a regulatory network for HCM. Drug prediction of these 10 hub IRGs proposed potential therapeutic agents that could be used in HCM. These results indicate that IRGs are potential regulators and drug therapeutic targets in HCM.

1. Introduction

Cardiomyopathy refers to a myocardial disease with cardiac dysfunction, among which hypertrophic cardiomyopathy (HCM) is the most common. HCM is typically characterized by abnormal left and/or right ventricular hypertrophy [1]. HCM is commonly known as hereditary cardiomyopathy [2, 3], which is estimated that more than 16 genes and 900 mutations are responsible for the occurrence of HCM [4–6]. Currently, HCM treatment mainly focuses on gene therapy [7], immune-modulatory therapy [8], and regenerative therapy [9–11].

An abundance of evidence shows that the immune system with the multiple therapeutic targets for HCM plays a pivotal role in cardiac tissue damage and repair [12–15] by blocking the proinflammatory pathway [16, 17] or regulating the circulation of immune cells and the regeneration of resident immune cells [18, 19]. Accumulating evidence indicates that the immune system is also involved in cardiac regeneration and has become a new hotspot in HCM regenerative therapy [20]. To date, treating MYBPC3 mutations is the only feasible gene therapies of HCM [7]. Nevertheless, the effect of IRGs on the process of HCM remains unclear.

In our study, we performed a systematic investigation of the HCM-related IRGs to identify prognostic biomarkers of HCM. The differentially expressed genes were obtained based on HCM, and the prognostic IRGs were further identified. Besides, functional analysis and the PPI network further revealed that these genes were implicated in the regulation of many biological processes, including SMAD protein pathways, smooth muscle cell proliferation, protein serine/threonine kinase, and MAPK cascade. Mechanistically, these regulatory pathways are involved in cytokine-cytokine receptor interaction, TGFβ signaling pathway, PI3K-Akt signaling pathway, and MAPK signaling pathway. Finally, the abundance of 10 hub IRGs was screened, and these ten hub IRGs were used for drug prediction. Meanwhile, these 10 hub IRGs were coexpressed with 4 TFs to construct a regulatory network for HCM. Altogether, we provided new insights into the regulatory mechanism of the IRGs in HCM, which may be a new direction for the treatment of HCM.

2. Materials and Methods

2.1. Data Acquisition

Transcriptome RNA-sequencing data of HCM (GSE89714, analyzed on Illumina HiSeq 2000) was downloaded from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo), which contained data from 5 hypertrophic heart tissue samples and 4 normal samples. The IRGs were obtained from the ImmPort database (http://www.immport.org), including 17 immune categories according to their different molecular function, such as T-cell receptor signaling pathway, B-cell receptor signaling pathway, Natural Killer Cell Cytotoxicity, antimicrobials, cytokine, and TNF family receptors. The list of TFs was getting from the Cistrome database. The function of the IRGs was acquired from Genecards (http://www.genecards.org/).

2.2. Identification of Differentially Expressed Genes (DEGs)

The R package “limma” and “voom” functions are used to normalize the data. DEGs between HCM and normal samples were screened using the limma package under the condition of Log fold and adj. value < 0.05.

2.3. KEGG Pathway Enrichment Analysis

Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource for understanding advanced functions and biological systems from large-scale molecular data generated by high-throughput experimental techniques. KEGG database was used to identify upregulation pathways by DAVID (version 6.8, https://david.ncifcrf.gov/) in the background of Homo sapiens. The pathway with a corrected value < 0.05 was considered a significant difference.

2.4. GO Enrichment Analysis

Gene ontology (GO) is a major bioinformatics tool for annotating genes and analyzing the biological processes (BPs), cellular components (CCs), and molecular functions (MFs), respectively. GO analysis was performed by DAVID (version 6.8, https://david.ncifcrf.gov/), and value < 0.05 was considered statistically significant.

2.5. PPI Network Analysis

The PPI regulatory network of IRGs was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING; http://string-db.org) database (version 10.0), and interactions with a combined were considered statistically significant. The PPI network result was displayed using the Cytoscape software version 3.4.0, and the most important modules were identified by Cytoscape’s plugin molecular complex detection (MCODE) (version 1.4.2). The selection criteria of differentially expressed IRGs in HCM were as follows: MCODE , degree , node score , maximum , and . value < 0.01 was considered statistically significant. The top 10 differentially expressed IRGs act as the hub genes by ranking the degree number of genes in the PPI network.

2.6. Drug Prediction

We used three pharmacogenomic databases for drug prediction of 10 hub genes: Drugbank online website (https://go.drugbank.com/), Drug-Gene interaction database (https://dgidb.genome.wustl.edu/), and Connectivity map (https://portals.broadinstitute.org/cmap/).

2.6.1. Drugbank Online Website

The proteins encoded by the 10 hub genes mentioned above were used in Drugbank online website for drug prediction (https://go.drugbank.com/), and the hub gene-protein-drug network was constructed by the Cytoscape software.

2.6.2. Drug-Gene Interaction Database

These 10 hub genes directly predicted drugs through the Drug-Gene interaction database (https://dgidb.genome.wustl.edu/) and visualized with the Cytoscape software.

2.6.3. Connectivity Map

These 10 hub genes were used for drug prediction with Connectivity map (https://portals.broadinstitute.org/cmap/); drugs with connectivity score below -0.879 (i.e., drugs most likely to reverse gene differential expression therapy for DCM) were selected. The connectivity score of each drug for each hub gene in the specific cell line was summarized; a heat map was made with the R software.

3. Results

3.1. Screening Results of Differential Genes and the Enrichment of Gene Function Analysis in HCM

GSE89714 was downloaded from the GEO database, which contained data from 5 hypertrophic heart tissue samples and 4 normal myocardial samples. 854 differentially expressed genes (DEGs) were screened, including 675 upregulated genes and 179 downregulated genes (Figure 1).

(a)

(b)

After functional enrichment analysis of 854 DEGs in HCM, we found that these 854 DEGs mainly participated in 10 biological processes (BP), including extracellular matrix organization, extracellular structure organization, skeletal system development, ossification, collagen fibril organization, urogenital system development, axon development, bone development, connective tissue development, and renal system development (Figure 2(a)). Additionally, these 854 DEGs were mainly expressed in the extracellular matrix, platelet alpha granule, basement membrane, endoplasmic reticulum lumen, collagen trimer, the complex of collagen trimers, and the cell-substrate junction (Figure 2(b)). Specifically, the molecular function (MF) of these 854 DEGs was mainly related to extracellular matrix structural constituent, collagen binding, glycosaminoglycan binding, extracellular matrix structural constituent conferring tensile strength, extracellular matrix structural constituent conferring compression resistance, heparin-binding, growth factor binding, growth factor activity, sulfur compound binding, and integrin binding (Figure 2(c)).

(a)

(b)

(c)

(d)

Subsequently, the results of the KEGG pathway show that there are 11 significantly upregulated pathways, including ECM-receptor interaction, focal adhesion, TGFβ signaling pathway, protein digestion and absorption, proteoglycans in cancer, PI3K-Akt signaling pathway, dilated cardiomyopathy (DCM), AGE-RAGE signaling pathway in diabetic complications, hypertrophic cardiomyopathy (HCM), amoebiasis, and cytokine-cytokine receptor interaction. The upregulated expression of the pathway is shown in Figure 2(d).

3.2. IRG Expression Profile of HCM and the Enrichment of Gene Function Analysis in HCM

The list of IRGs was downloaded through the ImmPort database. Meanwhile, the differentially expressed IRGs were extracted from GSE89714. 88 IRGs were differentially expressed in HCM. The expression levels of these 88 genes are listed in Supplementary Table S1, of which 16 are low expression and 72 are high expression (Figure 3(a)).

(a)

(b)

(c)

(d)

(e)

Figure 3

Differentially expressed IRGs in HCM and gene functional enrichment. (a) The differentially expressed IRGs between HCM and normal tissues. Green dots represent lower expressed genes, and red dots represent highly expressed genes. (b) Gene ontology analysis: biological process of differentially expressed IRGs. (c) Gene ontology analysis: cellular component of differentially expressed IRGs. (d) Gene ontology analysis: molecular function of differentially expressed IRGs. (e) Encyclopedia of Genes and Genomes analysis: the top 10 significant pathways of differentially expressed IRGs.

Through GO function enrichment analysis of 88 differentially expressed IRGs in HCM, we found that these 88 IRGs are mainly associated with the regulation of 10 biological processes (BP), including 4 processes of SMAD protein, 3 processes of smooth muscle cell proliferation, 2 processes of protein serine/threonine kinase, and MAPK cascade process (Figure 3(b)). Moreover, these 88 IRGs are mainly enriched in 10 Cell Components (CC), especially in platelet, extracellular matrix, vesicle lumen, granule lumen, mast cell, and so on (Figure 3(c)). Furthermore, the MF of these 88 IRGs is mainly correlated with receptor-ligand/regulator activity, growth factor activity/binding, cytokine activity/binding, hormone activity, G protein-coupled peptide receptor activity, and transmembrane receptor protein serine/threonine kinase binding (Figure 3(d)).

Thereafter, KEGG pathway enrichment analysis indicated that there are 30 significantly upregulated pathways in the microarray datasets, including cytokine-cytokine receptor interaction, TGFβ signaling pathway, viral protein interaction with cytokine and cytokine receptor, PI3K-Akt signaling pathway, pluripotency of stem cells regulation, MAPK signaling pathway, Hippo signaling pathway, proteoglycans in cancer, Ras signaling pathway, malaria, chemokine signaling pathway, calcium signaling pathway, hematopoietic cell lineage, amoebiasis, Rap1 signaling pathway, neuroactive ligand-receptor interaction, toxoplasmosis, cGMP-PKG signaling pathway, renal cell carcinoma, melanoma, axon guidance, Leishmaniasis, EGFR tyrosine kinase inhibitor resistance, bladder cancer, focal adhesion, gastric cancer, dilated cardiomyopathy (DCM), cAMP signaling pathway, AGE-RAGE signaling pathway in diabetic complications, and JAK-STAT signaling pathway. The results are shown in Figure 3(e). Importantly, cytokine-cytokine receptor interaction, TGFβ signaling pathway, PI3K-Akt signaling pathway, and MAPK signaling pathway are the most obvious upregulation pathways (Figure 3(e)).

3.3. Screening of Differentially Expressed TFs and Coexpression Analysis of TFs and IRGs

The PPI regulatory network of IRGs was based on data collected from the STRING online database (https://string-db.org/). We found that IRGs interacted directly with each other through the analysis of the PPI network (Figure 4). The list of TFs was downloaded through the Cistrome database and screened out 11 differentially expressed TFs in HCM, including POLR3G, EHF, SOX2, TP63, SOX4, KLF5, ESR1, RUNX1, PRDM1, MAFK, and MYH11(Figure 5(a)). Next, coexpression analysis results of TFs and IRGs showed that there were 5 TFs coexpressed with 84 IRGs, including RUNX1, SOX4, KLF5, MAFK, and MYH11(Figure 5(b)). Notably, the 5 TFs are highly expressed in HCM. The expression levels are shown in Table 1. Finally, the regulatory network between 5 TFs and 84 IRGs was constructed by using the software of Cytoscape.

(a)

(b)

3.4. Identification of Hub IRGs and TF Regulatory Network of IRGs

The top 10 IRGs (EGF, IGF1, HGF, BMP4, CXCR4, CTGF, THBS1, TGFB2, JAK2, TGFB3) were selected as the hub genes. Obviously, they also interacted with each other (Figure 6(a)), especially for IGF1 and HGF, THBS1 and IGF1, THBS1 and HGF, THBS1 and BMP4, CTGF and IGF1, CTGF and HGF, TGFB3 and THBS1, and TGFB3 and CTGF (Figure 6(c)). Except for EGF, the other 9 genes were highly expressed in HCM (Table 1). The functional roles of 10 hub genes were downloaded from the online database (https://www.genecards.org/) and are shown in Table 2. Besides, coexpression analysis results showed that there were 4 TFs with a crucial regulatory role in 10 hub IRGs network, including SOX4, KLF5, MAFK, and MYH11 (Figure 6(b)).

(a)

(b)

(c)

3.5. Drug Prediction of Hub Genes

We used multiple pharmacogenomic databases for drug prediction of 10 hub genes and found that different prediction methods have different results. In Drugbank online website (https://go.drugbank.com/), there are proteins encoded by only 7 hub genes yielded related drugs; proteins encoded by other 3 hub genes have not related drugs (Figure 7(a)). In the Drug-Gene interaction database (https://dgidb.genome.wustl.edu/), all hub genes resulted in the corresponding drugs; only FRESOLIMUMAB and ABT-510 were mainly correspond to two genes, while other drugs are all only correspond to one gene (Figure 7(b)). In the Connectivity map (https://portals.broadinstitute.org/cmap/), 13 drugs with connectivity score below -0.879 were predicted, summarizing the connectivity scores of these 13 drugs with 10 hub genes and making a heat map with the R software (Figure 7(c)).

(a)

(b)

(c)

Figure 7

Drug prediction of 10 hub IRGs. (a) Hub gene-protein-drug network. Red ovals represent hub genes, green triangles represent proteins encoded by hub genes, and blue rectangles are drugs. (b) Hub gene-drug network. Red indicates hub genes, and blue indicates drug. Among them, the 78 drugs corresponding to JAK2 showed only the top 30, and all 78 drug names can be obtained in Supplementary Table S2. (c) Heat map of the effect of drugs on hub genes expression. Green indicates that the drug can downregulate the hub gene expression, red marks that the drug can upregulate the hub gene expression, and black means that the drug and the hub gene expression have little correlation.

4. Discussion

HCM is an inherited myocardial disease with cardiac dysfunction and genetic disorder [21–23]. Accumulating evidence revealed that the immune system plays a key role in heart injury responses [24–26]. It has been reported that myocardial aging is largely associated with changes in gene expression patterns of immune response [27]. Therefore, it is plausible that the development of HCM is attributed to abnormal IRG expression. In our results, a total of 854 mRNAs were differentially expressed in HCM patients, of which 88 were IRGs, including 16 low-expressed IRGs and 72 high-expressed IRGs. 10 hub IRGs (IGF1, CTGF, TGFB2, TGFB3, HGF, BMP4, CXCR4, THBS1, EGF, JAK2) were selected to construct a regulatory network in HCM, suggesting that these 10 hub IRGs might be the new therapeutic targets for the intervention of HCM.

Insulin-like growth factors 1 (IGF1) is a hormone with pleiotropic effects, which consists of 70 amino acids and is regulated by IGF-binding proteins (IGFBPs) [28–31]. It is generally accepted that IGF1 exhibits a protective effect on myocardium by regulating autophagy, precursor cell differentiation, extracellular matrix activity, and myocardial fibrosis [32, 33]. In our results, insulin-like growth factor 1 (IGF1) and IGF-binding protein (CTGF) are two hub IRGs in HCM. Moreover, CTGF and IGF1 were highly expressed in HCM and significantly interacted with each other(). Accordingly, IGF1 may provide a novel therapeutic direction for HCM. A plethora of evidence showed that IGF1 could regulate the metabolism of cardiomyocytes through MAPK [34–36] and PI3K-Akt pathway [37]. In the present study, GO function enrichment analysis showed that IRGs participate in the process of MAPK cascade, and the KEGG pathways are enriched in the PI3K-Akt signaling pathway and MAPK signaling pathway. Therefore, PI3K-Akt and MAPK signaling pathway may be the regulatory mechanism of IGF1 on the pathogenesis of HCM.

TGFβs are multifunctional cytokines with three isoforms (TGFβ1, 2, and 3) encoded by three distinct genes, which are essential mediators of cardiac repair [38–42]. Extensive evidence suggests that TGFβ is produced by various cells in response to tissue injury, especially in injured myocardial regions. In general, TGFβ isoforms show different expression patterns in the injured myocardium: TGFβ1 and TGFβ2 are induced in the early stage, whereas TGFβ3 shows a delayed and prolonged upregulation. Importantly, our results reveal that transforming growth factor β2 proprotein (TGFB2) and transforming growth factor β3 (TGFB3) are two hub IRGs in HCM. Moreover, TGFB2 and TGFB3 were highly expressed in HCM, indicating a promising therapeutic target of TGFβ for HCM. TGFβ can be activated by SMAD [43, 44]. Meanwhile, TGFβ-SMAD3 signaling pathway is an important mediator of cardiac fibroblasts and is expected to be an attractive therapeutic target for HCM [45–48]. In our results, these 88 IRGs are involved in the regulation of 4 biological processes of SMAD protein, including positive regulation of pathway-restricted SMAD protein phosphorylation, pathway-restricted SMAD protein phosphorylation, regulation of pathway−restricted SMAD protein phosphorylation, and SMAD protein signal transduction. Furthermore, the results of KEGG pathway enrichment analysis showed that the TGFβ signaling pathway is upregulated. Collectively, the above results suggest that two hub IRGs (TGFB2 and TGFB3) regulated the TGFβ-SMAD signaling pathway and participate in the pathological process of HCM.

In the regulatory networks of the 10 hub IRGs, 4 TFs (SOX4, KLF5, MAFK, and MYH11) were coexpressed with these 10 hub IRGs, which suggests that SOX4, KLF5, MAFK, and MYH11 are the major regulators of IRGs during HCM. In detail, SOX4 is a member of the Sox (for Sry-box) family, which is essential for the regulation of embryonic development and cell fate. Also, SOX4 is required for the development of cardiac outflow tract. SOX4 deficiency even causes congenital heart defects in embryonic mice [49–51]. As a member of the Krüppel-like factors (KLFs) family, KLF5 is intimately correlated with cardiovascular remodeling, cardiac energetic, and cardiac fibroblasts [52–55]. Notably, KLF5 has recently been proved as a novel disease gene of familial dilated cardiomyopathy (DCM) and negatively regulates HCM [56, 57]. MAFK, a member of the small Maf family, is an important contributor to hematopoiesis [58]. MAFK is well known for its role in modulating NF-κB activity and controlling immune-inflammatory responses [59, 60]. MYH11 is a smooth muscle myosin belonging to the myosin heavy chain family, which plays a vital role in human disease familial hypertrophic cardiomyopathy (FHC) [61, 62]. In our results, these 4 TFs were not only coexpressed with 10 hub genes but also correlated with each other to form a regulatory network. The network may be a potential therapeutic target for HCM.

Recent studies revealed that the immune system plays a pivotal but complex role in myocardial injury and repair [63, 64]. During the cardiac injury, the immune mechanism is activated, which induces an inflammatory response, and subsequently triggers reparative pathways to promote wound healing [65]. Surprisingly, the overall mobilization of the immune system follows a consistent pattern in response to tissue injury: when an injury occurs, the dying cardiomyocytes trigger an acute inflammatory response and are recognized by resident immune cells, such as macrophages, innate lymphoid cells (ILC), and mast cells (MC) [66]. Subsequently, these resident immune cells also participate in the removal of dying cardiomyocytes, the reconstruction of tissue matrix, the restriction of the spread of inflammatory molecules, and the maintenance of local homeostasis [67, 68]. Besides, macrophages/monocytes secrete various chemokines and cytokines that recruit immune cells to promote wound healing [69]. These processes are highly consistent with the cell constituents (CC) of 88 IRGs, which are differentially expressed in our results. In our results, 88 IRGs are primarily enriched in the collagen-containing extracellular matrix, secretory granule lumen, cytoplasmic vesicle lumen, mast cell granule, and cell-substrate adherens junction.

In this study, we used three pharmacogenomic databases for drug prediction of 10 hub genes and found that different prediction methods have different results. In Drugbank online website (Figure 7(a)), proteins encoded by only 7 hub genes yielded related drugs. In the hub gene-drug network (Figure 7(b)), there are two drugs FRESOLIMUMAB and ABT-510, corresponding to two hub genes at the same time, indicating that these two drugs may be vital therapeutic molecules. Connectivity map prediction results (Figure 7(c)) clearly showed that all 13 drugs could upregulate the differentially underexpressed hub genes, that is, EGF. All these 13 drugs could downregulate most of the differentially highly expressed 9 hub genes. These drugs may be potentially effective therapeutic molecules in the future and may inspire the design of similarly structured active drug molecules.

Generally, cardiac injury activates immune mechanisms, which depend on the close interaction of cells, the regeneration of cardiac resident immune cells, and immune cells [70, 71]. In the process of myocardial injury and repair, immune cells, like CD4⁺ T cells, CD8⁺ T cells, B cells, natural killer cells (NK), and macrophages, are the main regulator in wound healing, which includes resident immune cells and circulating immune cells from the blood [72–74]. As mentioned previously, the resident immune cells coordinate the removal of dying cardiomyocytes, the reconstruction of tissue matrix, the restriction of the spread of inflammatory molecules, and the maintenance of local homeostasis. In addition, excessive inflammatory signals can also trigger extramedullary hematopoiesis of the spleen and produce new immune cells that are recruited to the injured site to alleviate tissue injury [75]. Therefore, the balance between self-maintenance of the cardiac resident immune cells and the recruitment of circulating immune cells is crucial for the potential of cardiac regeneration [76]. At present, the replacement of lost cardiomyocytes is the most important strategy for the treatment of cardiomyopathy. In the past, the regenerative therapy of HCM is mainly focused on the “stem cell,” especially for pluripotent stem cell-derived cardiomyocytes [77, 78]. However, it seems that “stromal cell types” have got more attention, particularly including resident immune cells and innate immune system. Therefore, the cardiac resident immune cells become a hotspot in cardiac regeneration of HCM [20]. Moreover, macrophages are also implicated in heart regeneration [79, 80]. In our results, the KEGG pathway of 88 differentially expressed IRGs is also mostly enriched in signaling pathways that regulate pluripotency of stem cells and hematopoietic cell lineage, suggesting that IRGs might be another hotspot in cardiac regeneration of HCM by regulating pluripotency of stem cells and hematopoietic cell lineage.

5. Conclusions

IRGs are a key player in the progression of hypertrophic cardiomyopathy, which is mainly involved in immune regulation, cardiac regeneration, and gene therapy. Meanwhile, 10 hub IRGs were coexpressed with 4 TFs to construct a regulatory network for HCM. Mechanistically, PI3K-Akt and MAPK signaling pathway are the main regulatory mechanism of IGF1, and TGFβ-SMAD signaling pathway is another regulatory mechanism in the pathological process of HCM. Drug prediction of these 10 hub IRGs proposed potential therapeutic agents that could be used in HCM. All these data suggest that IRGs may be a new direction for the treatment of HCM.

Data Availability

The data used to support the findings of this study are available in the Gene Expression Omnibus (GEO) repository (https://www.ncbi.nlm.nih.gov/gds).

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

This study was supported by the Scientific research project of Hunan Health Commission (No. 202102041763, No. 20200985), the Natural Science Foundation of Changsha (No. 202045504), and the Hunan Cancer Hospital Climb Plan.

Supplementary Materials

Supplementary 1. Table S1: the expression levels of 88 IRGs.

Supplementary 2. Table S2: 78 drug names corresponding to JAK2.

References

G. Thiene, D. Corrado, and C. Basso, “Revisiting definition and classification of cardiomyopathies in the era of molecular medicine,” European Heart Journal, vol. 29, no. 2, pp. 144–146, 2008.
View at: Publisher Site | Google Scholar
Y. Hata, K. Hirono, Y. Yamaguchi, F. Ichida, Y. Oku, and N. Nishida, “Minimal inflammatory foci of unknown etiology may be a tentative sign of early stage inherited cardiomyopathy,” Modern Pathology, vol. 32, no. 9, pp. 1281–1290, 2019.
View at: Publisher Site | Google Scholar
S. Abdullah, M. L. Lynn, M. T. McConnell et al., “FRET-based analysis of the cardiac troponin T linker region reveals the structural basis of the hypertrophic cardiomyopathy-causing Δ160E mutation,” The Journal of Biological Chemistry, vol. 294, no. 40, pp. 14634–14647, 2019.
View at: Publisher Site | Google Scholar
C. Lu, W. Wu, F. Liu et al., “Molecular analysis of inherited cardiomyopathy using next generation semiconductor sequencing technologies,” Journal of translational medicine, vol. 16, no. 1, p. 241, 2018.
View at: Google Scholar
N. Adadi, F. Z. Radi, N. Lahrouchi et al., “Inherited dilated cardiomyopathy in a large Moroccan family caused by LMNA mutation,” Anatolian Journal of Cardiology, vol. 20, no. 1, pp. 65–68, 2018.
View at: Publisher Site | Google Scholar
C. D. Vera, C. A. Johnson, J. Walklate et al., “Myosin motor domains carrying mutations implicated in early or late onset hypertrophic cardiomyopathy have similar properties,” The Journal of Biological Chemistry, vol. 294, no. 46, pp. 17451–17462, 2019.
View at: Publisher Site | Google Scholar
M. Prondzynski, G. Mearini, and L. Carrier, “Gene therapy strategies in the treatment of hypertrophic cardiomyopathy,” Pflügers Archiv, vol. 471, no. 5, pp. 807–815, 2019.
View at: Publisher Site | Google Scholar
B. Kardaszewicz, E. Rogala, M. Tendera, P. Kardaszewicz, and J. Jarzab, “Circulating immune complexes in hypertrophic cardiomyopathy and ischemic heart disease,” Kardiologia Polska, vol. 34, no. 1, pp. 21–24, 1991.
View at: Google Scholar
Y. J. Woo, C. M. Panlilio, R. K. Cheng et al., “Myocardial regeneration therapy for ischemic cardiomyopathy with cyclin A2,” The Journal of Thoracic and Cardiovascular Surgery, vol. 133, no. 4, pp. 927–933, 2007.
View at: Publisher Site | Google Scholar
G. D. Richardson, S. Laval, and W. A. Owens, “Cardiomyocyte regeneration in the mdx mouse model of nonischemic cardiomyopathy,” Stem Cells and Development, vol. 24, no. 14, pp. 1672–1679, 2015.
View at: Publisher Site | Google Scholar
F. Pagano, V. Picchio, I. Chimenti et al., “On the road to regeneration: "tools" and "routes" towards efficient cardiac cell therapy for ischemic cardiomyopathy,” Current cardiology reports, vol. 21, no. 11, p. 133, 2019.
View at: Google Scholar
A. Saparov, V. Ogay, T. Nurgozhin et al., “Role of the immune system in cardiac tissue damage and repair following myocardial infarction,” Inflammation Research, vol. 66, no. 9, pp. 739–751, 2017.
View at: Publisher Site | Google Scholar
P. van den Hoogen, S. C. A. de Jager, M. M. H. Huibers et al., “Increased circulating IgG levels, myocardial immune cells and IgG deposits support a role for an immune response in pre- and end-stage heart failure,” Journal of Cellular and Molecular Medicine, vol. 23, no. 11, pp. 7505–7516, 2019.
View at: Publisher Site | Google Scholar
S. Sattler and N. Rosenthal, “The neonate versus adult mammalian immune system in cardiac repair and regeneration,” Biochimica et Biophysica Acta (BBA)-Molecular Cell Research, vol. 1863, no. 7 Part B, pp. 1813–1821, 2016.
View at: Google Scholar
N. G. Frangogiannis, “The immune system and cardiac repair,” Pharmacological Research, vol. 58, no. 2, pp. 88–111, 2008.
View at: Publisher Site | Google Scholar
S. C. Latet, V. Y. Hoymans, P. L. Van Herck, and C. J. Vrints, “The cellular immune system in the post-myocardial infarction repair process,” International Journal of Cardiology, vol. 179, pp. 240–247, 2015.
View at: Publisher Site | Google Scholar
A. Zimmer, A. K. Bagchi, K. Vinayak, A. Bello-Klein, and P. K. Singal, “Innate immune response in the pathogenesis of heart failure in survivors of myocardial infarction,” American Journal of Physiology. Heart and Circulatory Physiology, vol. 316, no. 3, pp. H435–H445, 2019.
View at: Publisher Site | Google Scholar
S. Sattler, P. Fairchild, F. M. Watt, N. Rosenthal, and S. E. Harding, “The adaptive immune response to cardiac injury-the true roadblock to effective regenerative therapies?” NPJ Regenerative medicine, vol. 2, no. 1, p. 19, 2017.
View at: Publisher Site | Google Scholar
R. Gentek and G. Hoeffel, “The innate immune response in myocardial infarction, repair, and regeneration,” Advances in Experimental Medicine and Biology, vol. 1003, pp. 251–272, 2017.
View at: Publisher Site | Google Scholar
S. J. Forbes and N. Rosenthal, “Preparing the ground for tissue regeneration: from mechanism to therapy,” Nature Medicine, vol. 20, no. 8, pp. 857–869, 2014.
View at: Publisher Site | Google Scholar
L. Fang, A. H. Ellims, A. L. Beale, A. J. Taylor, A. Murphy, and A. M. Dart, “Systemic inflammation is associated with myocardial fibrosis, diastolic dysfunction, and cardiac hypertrophy in patients with hypertrophic cardiomyopathy,” American Journal of Translational Research, vol. 9, no. 11, pp. 5063–5073, 2017.
View at: Google Scholar
W. E. Moody, M. Schmitt, and P. Arumugam, “Coronary microvascular dysfunction in hypertrophic cardiomyopathy detected by Rubidium-82 positron emission tomography and cardiac magnetic resonance imaging,” Journal of Nuclear Cardiology, vol. 26, no. 2, pp. 666–670, 2019.
View at: Publisher Site | Google Scholar
A. Doolan, L. Nguyen, and C. Semsarian, “Hypertrophic cardiomyopathy: from “heart tumour” to a complex molecular genetic disorder,” Heart, Lung and Circulation, vol. 13, no. 1, pp. 15–25, 2004.
View at: Publisher Site | Google Scholar
S. Epelman and D. L. Mann, “Communication in the heart: the role of the innate immune system in coordinating cellular responses to ischemic injury,” Journal of Cardiovascular Translational Research, vol. 5, no. 6, pp. 827–836, 2012.
View at: Publisher Site | Google Scholar
G. Baumgarten, S. C. Kim, H. Stapel et al., “Myocardial injury modulates the innate immune system and changes myocardial sensitivity,” Basic Research in Cardiology, vol. 101, no. 5, pp. 427–435, 2006.
View at: Publisher Site | Google Scholar
A. A. Knowlton, “Paying for the tolls: the high cost of the innate immune system for the cardiac myocyte,” Advances in Experimental Medicine and Biology, vol. 1003, pp. 17–34, 2017.
View at: Publisher Site | Google Scholar
B. Bartling, K. Niemann, R. U. Pliquett, H. Treede, and A. Simm, “Altered gene expression pattern indicates the differential regulation of the immune response system as an important factor in cardiac aging,” Experimental Gerontology, vol. 117, pp. 13–20, 2019.
View at: Publisher Site | Google Scholar
T. Halmos and I. Suba, “A növekedési hormon és az inzulinszerű növekedési faktorok élettani szerepe,” Orvosi Hetilap, vol. 160, no. 45, pp. 1774–1783, 2019.
View at: Publisher Site | Google Scholar
S. Pouriamehr, H. Barmaki, M. Rastegary, F. Lotfi, and M. Nabi Afjadi, “Investigation of insulin-like growth factors/insulin-like growth factor binding proteins regulation in metabolic syndrome patients,” BMC research notes, vol. 12, no. 1, p. 653, 2019.
View at: Google Scholar
N. R. Bhakta, A. M. Garcia, E. H. Frank, A. J. Grodzinsky, and T. I. Morales, “The insulin-like growth factors (IGFs) I and II bind to articular cartilage via the IGF-binding proteins,” The Journal of Biological Chemistry, vol. 275, no. 8, pp. 5860–5866, 2000.
View at: Publisher Site | Google Scholar
J. L. Martin and R. C. Baxter, “Signalling pathways of insulin-like growth factors (IGFs) and IGF binding protein-3,” Growth Factors, vol. 29, no. 6, pp. 235–244, 2011.
View at: Publisher Site | Google Scholar
C. Iosef Husted and M. Valencik, “Insulin-like growth factors and their potential role in cardiac epigenetics,” Journal of Cellular and Molecular Medicine, vol. 20, no. 8, pp. 1589–1602, 2016.
View at: Publisher Site | Google Scholar
J. Balcells, A. Moreno, L. Audi, J. Roqueta, J. Iglesias, and A. Carrascosa, “Growth hormone/insulin-like growth factors axis in children undergoing cardiac surgery,” Critical Care Medicine, vol. 29, no. 6, pp. 1234–1238, 2001.
View at: Publisher Site | Google Scholar
K. J. Pollard and J. M. Daniel, “Nuclear estrogen receptor activation by insulin-like growth factor-1 in Neuro-2A neuroblastoma cells requires endogenous estrogen synthesis and is mediated by mutually repressive MAPK and PI3K cascades,” Molecular and Cellular Endocrinology, vol. 490, pp. 68–79, 2019.
View at: Publisher Site | Google Scholar
W. Wang, D. Wen, W. Duan et al., “Systemic administration of scAAV9-IGF1 extends survival in SOD1G93A ALS mice via inhibiting p38 MAPK and the JNK-mediated apoptosis pathway,” Brain Research Bulletin, vol. 139, pp. 203–210, 2018.
View at: Publisher Site | Google Scholar
T. L. C. Wolters, M. G. Netea, A. Hermus, J. W. A. Smit, and R. T. Netea-Maier, “IGF1 potentiates the pro-inflammatory response in human peripheral blood mononuclear cells via MAPK,” Journal of Molecular Endocrinology, vol. 59, no. 2, pp. 129–139, 2017.
View at: Publisher Site | Google Scholar
P. Massoner, M. Ladurner-Rennau, I. E. Eder, and H. Klocker, “Insulin-like growth factors and insulin control a multifunctional signalling network of significant importance in cancer,” British Journal of Cancer, vol. 103, no. 10, pp. 1479–1484, 2010.
View at: Publisher Site | Google Scholar
G. Nicolini, F. Forini, C. Kusmic, L. Pitto, L. Mariani, and G. Iervasi, “Early and short-term triiodothyronine supplementation prevents adverse postischemic cardiac remodeling: role of transforming growth factor-β1 and antifibrotic miRNA signaling,” Molecular Medicine, vol. 21, no. 1, pp. 900–911, 2016.
View at: Publisher Site | Google Scholar
M. Dobaczewski, W. Chen, and N. G. Frangogiannis, “Transforming growth factor (TGF)-β signaling in cardiac remodeling,” Journal of Molecular and Cellular Cardiology, vol. 51, no. 4, pp. 600–606, 2011.
View at: Publisher Site | Google Scholar
H. E. Mewhort, B. D. Lipon, D. A. Svystonyuk et al., “Monocytes increase human cardiac myofibroblast-mediated extracellular matrix remodeling through TGF-β1,” American Journal of Physiology. Heart and Circulatory Physiology, vol. 310, no. 6, pp. H716–H724, 2016.
View at: Publisher Site | Google Scholar
A. V. Shinde and N. G. Frangogiannis, “Fibroblasts in myocardial infarction: a role in inflammation and repair,” Journal of Molecular and Cellular Cardiology, vol. 70, pp. 74–82, 2014.
View at: Publisher Site | Google Scholar
C. Gerarduzzi and J. A. Di Battista, “Myofibroblast repair mechanisms post-inflammatory response: a fibrotic perspective,” Inflammation Research, vol. 66, no. 6, pp. 451–465, 2017.
View at: Publisher Site | Google Scholar
C. M. Tinoco-Veras, A. Santos, J. Stipursky et al., “Transforming growth factor beta1/SMAD signaling pathway activation protects the intestinal epithelium from Clostridium difficile toxin A-induced damage,” Infection and immunity, vol. 85, no. 10, 2017.
View at: Google Scholar
H. L. Ma, X. F. Zhao, G. Z. Chen, R. H. Fang, and F. R. Zhang, “Silencing NLRC5 inhibits extracellular matrix expression in keloid fibroblasts via inhibition of transforming growth factor-β1/Smad signaling pathway,” Biomedicine & Pharmacotherapy, vol. 83, pp. 1016–1021, 2016.
View at: Publisher Site | Google Scholar
M. Liu, J. Ai, J. Feng et al., “Effect of paeoniflorin on cardiac remodeling in chronic heart failure rats through the transforming growth factor β1/Smad signaling pathway,” Cardiovascular diagnosis and therapy, vol. 9, no. 3, pp. 272–280, 2019.
View at: Publisher Site | Google Scholar
C. Y. Zhan, J. H. Tang, D. X. Zhou, and Z. H. Li, “Effects of tanshinone IIA on the transforming growth factor β1/Smad signaling pathway in rat cardiac fibroblasts,” Indian journal of pharmacology, vol. 46, no. 6, pp. 633–638, 2014.
View at: Publisher Site | Google Scholar
D. Seo and J. M. Hare, “The transforming growth factor-beta/Smad3 pathway: coming of age as a key participant in cardiac remodeling,” Circulation, vol. 116, no. 19, pp. 2096–2098, 2007.
View at: Publisher Site | Google Scholar
F. Yang, A. C. Chung, X. R. Huang, and H. Y. Lan, “Angiotensin II induces connective tissue growth factor and collagen I expression via transforming growth factor-beta-dependent and -independent Smad pathways: the role of Smad3,” Hypertension, vol. 54, no. 4, pp. 877–884, 2009.
View at: Publisher Site | Google Scholar
M. H. Paul, R. P. Harvey, M. Wegner, and E. Sock, “Cardiac outflow tract development relies on the complex function of Sox4 and Sox11 in multiple cell types,” Cellular and Molecular Life Sciences, vol. 71, no. 15, pp. 2931–2945, 2014.
View at: Publisher Site | Google Scholar
J. Ya, M. W. Schilham, P. A. de Boer, A. F. Moorman, H. Clevers, and W. H. Lamers, “Sox4-deficiency syndrome in mice is an animal model for common trunk,” Circulation Research, vol. 83, no. 10, pp. 986–994, 1998.
View at: Publisher Site | Google Scholar
A. Restivo, G. Piacentini, S. Placidi, C. Saffirio, and B. Marino, “Cardiac outflow tract: a review of some embryogenetic aspects of the conotruncal region of the heart,” The Anatomical Record. Part A, Discoveries in Molecular, Cellular, and Evolutionary Biology, vol. 288, no. 9, pp. 936–943, 2006.
View at: Google Scholar
K. Drosatos, N. M. Pollak, C. J. Pol et al., “Cardiac myocyte KLF5 regulates Ppara expression and cardiac function,” Circulation Research, vol. 118, no. 2, pp. 241–253, 2016.
View at: Publisher Site | Google Scholar
N. D. Roe, S. W. Standage, and R. Tian, “The relationship between KLF5 and PPARα in the heart,” Circulation Research, vol. 118, no. 2, pp. 193–195, 2016.
View at: Publisher Site | Google Scholar
N. Takeda, I. Manabe, Y. Uchino et al., “Cardiac fibroblasts are essential for the adaptive response of the murine heart to pressure overload,” The Journal of Clinical Investigation, vol. 120, no. 1, pp. 254–265, 2010.
View at: Publisher Site | Google Scholar
R. Nagai, T. Suzuki, K. Aizawa, T. Shindo, and I. Manabe, “Significance of the transcription factor KLF5 in cardiovascular remodeling,” Journal of Thrombosis and Haemostasis, vol. 3, no. 8, pp. 1569–1576, 2005.
View at: Publisher Site | Google Scholar
R. Nagai, T. Shindo, I. Manabe, T. Suzuki, and M. Kurabayashi, “KLF5/BTEB2, a Kruppel-like zinc-finger type transcription factor, mediates both smooth muscle cell activation and cardiac hypertrophy,” Advances in Experimental Medicine and Biology, vol. 538, pp. 57–65, 2003.
View at: Google Scholar
R. M. Di, C. X. Yang, C. M. Zhao et al., “Identification and functional characterization of KLF5 as a novel disease gene responsible for familial dilated cardiomyopathy,” European journal of medical genetics, vol. 63, article 103827, 2019.
View at: Google Scholar
K. Igarashi, K. Itoh, H. Motohashi et al., “Activity and expression of murine small Maf family protein MafK,” The Journal of Biological Chemistry, vol. 270, no. 13, pp. 7615–7624, 1995.
View at: Publisher Site | Google Scholar
F. Katsuoka, H. Motohashi, K. Onodera, N. Suwabe, J. D. Engel, and M. Yamamoto, “One enhancer mediates mafK transcriptional activation in both hematopoietic and cardiac muscle cells,” The EMBO Journal, vol. 19, no. 12, pp. 2980–2991, 2000.
View at: Publisher Site | Google Scholar
Y. J. Hwang, E. W. Lee, J. Song, H. R. Kim, Y. C. Jun, and K. A. Hwang, “MafK positively regulates NF-κB activity by enhancing CBP-mediated p65 acetylation,” Scientific Reports, vol. 3, no. 1, p. 3242, 2013.
View at: Publisher Site | Google Scholar
X. Feng, T. He, J. G. Wang, and P. Zhao, “Asn391Thr mutation of β-myosin heavy chain in a hypertrophic cardiomyopathy family,” International Heart Journal, vol. 59, no. 3, pp. 596–600, 2018.
View at: Publisher Site | Google Scholar
D. I. Keller, J. Schwitter, E. R. Valsangiacomo, P. Landolt, and C. H. Attenhofer Jost, “Hypertrophic cardiomyopathy due to beta-myosin heavy chain mutation with extreme phenotypic variability within a family,” International Journal of Cardiology, vol. 134, no. 3, pp. e87–e93, 2009.
View at: Publisher Site | Google Scholar
S. Epelman, P. P. Liu, and D. L. Mann, “Role of innate and adaptive immune mechanisms in cardiac injury and repair,” Nature Reviews. Immunology, vol. 15, no. 2, pp. 117–129, 2015.
View at: Publisher Site | Google Scholar
T. Kennedy-Lydon, “Immune functions and properties of resident cells in the heart and cardiovascular system: pericytes,” Advances in Experimental Medicine and Biology, vol. 1003, pp. 93–103, 2017.
View at: Publisher Site | Google Scholar
E. Turillazzi, C. Pomara, S. Bello, M. Neri, I. Riezzo, and V. Fineschi, “The meaning of different forms of structural myocardial injury, immune response and timing of infarct necrosis and cardiac repair,” Current Vascular Pharmacology, vol. 13, no. 1, pp. 6–19, 2015.
View at: Publisher Site | Google Scholar
M. Hulsmans, F. Sam, and M. Nahrendorf, “Monocyte and macrophage contributions to cardiac remodeling,” Journal of Molecular and Cellular Cardiology, vol. 93, pp. 149–155, 2016.
View at: Publisher Site | Google Scholar
D. Hashimoto, A. Chow, C. Noizat et al., “Tissue-resident macrophages self-maintain locally throughout adult life with minimal contribution from circulating monocytes,” Immunity, vol. 38, no. 4, pp. 792–804, 2013.
View at: Publisher Site | Google Scholar
S. C. Macri, C. C. Bailey, N. M. de Oca et al., “Immunophenotypic alterations in resident immune cells and myocardial fibrosis in the aging rhesus macaque (Macaca mulatta) heart,” Toxicologic Pathology, vol. 40, no. 4, pp. 637–646, 2012.
View at: Publisher Site | Google Scholar
F. K. Swirski and M. Nahrendorf, “Cardioimmunology: the immune system in cardiac homeostasis and disease,” Nature Reviews. Immunology, vol. 18, no. 12, pp. 733–744, 2018.
View at: Publisher Site | Google Scholar
N. V. Krylova, T. P. Smolina, and G. N. Leonova, “Molecular mechanisms of interaction between human immune cells and far eastern tick-borne encephalitis virus strains,” Viral Immunology, vol. 28, no. 5, pp. 272–281, 2015.
View at: Publisher Site | Google Scholar
F. Bonner, N. Borg, S. Burghoff, and J. Schrader, “Resident cardiac immune cells and expression of the ectonucleotidase enzymes CD39 and CD73 after ischemic injury,” PLoS One, vol. 7, no. 4, article e34730, 2012.
View at: Google Scholar
X. Xue and D. M. Falcon, “The role of immune cells and cytokines in intestinal wound healing,” International journal of molecular sciences, vol. 20, no. 23, 2019.
View at: Google Scholar
D. M. Zaiss, C. M. Minutti, and J. A. Knipper, “Immune- and non-immune-mediated roles of regulatory T-cells during wound healing,” Immunology, vol. 157, no. 3, pp. 190–197, 2019.
View at: Publisher Site | Google Scholar
M. V. Krasnoselsky, L. I. Simonova, V. Z. Gertman, E. S. Pushkar, and T. S. Zavadskaya, “Tissue immune cells and their role in the healing process of infected radiation ulcers under the impact of photodynamic therapy (experimental study),” Problemy radiatsiinoi medytsyny ta radiobiolohii, vol. 24, pp. 250–260, 2019.
View at: Publisher Site | Google Scholar
F. C. Simoes, T. J. Cahill, A. Kenyon et al., “Macrophages directly contribute collagen to scar formation during zebrafish heart regeneration and mouse heart repair,” Nature communications, vol. 11, no. 1, p. 600, 2020.
View at: Google Scholar
G. Caputa, A. Castoldi, and E. J. Pearce, “Metabolic adaptations of tissue-resident immune cells,” Nature Immunology, vol. 20, no. 7, pp. 793–801, 2019.
View at: Publisher Site | Google Scholar
X. Li, W. J. Lu, Y. Li et al., “MLP-deficient human pluripotent stem cell derived cardiomyocytes develop hypertrophic cardiomyopathy and heart failure phenotypes due to abnormal calcium handling,” Cell death & disease, vol. 10, no. 8, p. 610, 2019.
View at: Google Scholar
N. Aoyama, J. Liu, S. Lor, T. K. Feaster, S. Hilcove, and E. Jones, “Assessment of hypertrophic cardiomyopathy using human induced pluripotent stem cell-derived cardiomyocytes reveals abnormal excitation contraction coupling,” Journal of Pharmacological and Toxicological Methods, vol. 99, p. 106595, 2019.
View at: Publisher Site | Google Scholar
J. Leor, D. Palevski, U. Amit, and T. Konfino, “Macrophages and regeneration: lessons from the heart,” Seminars in Cell & Developmental Biology, vol. 58, pp. 26–33, 2016.
View at: Publisher Site | Google Scholar
A. B. Aurora, E. R. Porrello, W. Tan et al., “Macrophages are required for neonatal heart regeneration,” The Journal of Clinical Investigation, vol. 124, no. 3, pp. 1382–1392, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Wei Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

698

Downloads

878

Citations