Abstract
Objective: To identify feature immune-related genes (IRGs) in patients with hypertrophic cardiomyopathy (HCM) and verify their ability to diagnose HCM.
Methods: The GSE160997 dataset on cardiac tissue from 18 HCM patients and 5 controls was downloaded from the Gene Expression Omnibus database. A false discovery rate <0.05 and |log2 fold change| >1 were the filters applied to identify the differentially expressed genes (DEGs). The differentially expressed IRGs were the intersection results between the DEGs and an IRG dataset from the IMMPORT database. The protein-protein interaction network of differentially expressed IRGs was constructed, and the top 20 hub genes with the most adjacent nodes in the network were selected. The least absolute shrinkage and selection operator regression algorithm and a random forest algorithm were used to identify the feature IRGs as biomarkers that were then verified against GSE36961.
Results: A total of 1079 DEGs were identified in GSE160997. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses indicated that immune-related mechanisms play an important role in the pathogenesis of HCM. A total of 121 differentially expressed IRGs were identified, and 5 feature IRGs were selected, 4 of which were confirmed as potential biomarkers of HCM by external verification with excellent discrimination ability. A diagnosis model of HCM based on the four feature IRGs was developed and visualized as a nomogram with a C-index of 0.925 (95% confidence interval 0.869–0.981).
Conclusion: Our study identified four feature IRGs as biomarkers for the diagnosis of HCM, offering an innovative perspective of the underlying immune-related pathological molecular mechanisms.
Introduction
Hypertrophic cardiomyopathy (HCM) is a hereditary heart disease that is defined as unexplained, isolated, progressive myocardial hypertrophy, usually accompanied by heart failure or arrhythmia. It is considered the primary underlying condition in sudden cardiac death in young adults and athletes. Mutations in the genes encoding proteins of the cardiomyocytic contractile apparatus are considered the main cause of HCM. However, the relationship between sarcomere mutations and clinical outcomes remains unpredictable in patients with HCM, largely because of genetic and phenotypic heterogeneity (, ). During the past two decades, the continuous improvement of gene sequencing technology has enabled the diagnosis of HCM among populations with a family history and in asymptomatic patients through genetic testing to avoid sudden death. This has become a class I recommendation for supplementary testing according to the European Society of Cardiology and American Heart Association/American College of Cardiology guidelines (, ).
Sarcomere protein gene (such as MYBPC3, MYH7, TNNT2, TNNI3, TPM1, MYL2, MYL3, and ACTC) mutations have become targets in the genetic testing of HCM. Unfortunately, in the genetic testing of different HCM cohorts, up to 50–68% of patients were negative for sarcomere protein gene mutations (). The signaling pathways and regulatory networks underlying the pathogenesis of HCM remain unclear to date. Thus, a systematic analysis of its transcriptional map might provide useful information about the overall underlying molecular mechanisms and pathological processes of HCM and aid in the diagnosis and treatment of HCM.
As the flowchart showed in Figure 1, this study aimed to explore the high-throughput sequencing gene expression data of HCM patients and normal samples, construct the differential gene expression profile of HCM, annotate the gene ontology and potential signal pathways, and identify valuable biomarkers that may further the gene test diagnosis of HCM.
Figure 1
Materials and Methods
Data Source and Differentially Expressed Gene Analyses
There are five comparatively valuable datasets that contain more than 10 samples of HCM in GEO database, namely GSE160997 (), GSE143786 (), GSE148602 (), GSE130036 (), and GSE36961. Among these five datasets, GSE160997, GSE130036, and GSE36961 focus on transcriptome information and a previous study of our team had explored GSE130036 (). The high-throughput sequencing dataset GSE160997 used in this study contained data regarding myocardial tissue of 18 HCM patients and 5 healthy controls. All data were statistically analyzed and visualized using R software (version 4.0.3) and its relevant packages. The original data of all samples were normalized and differentially expressed genes (DEGs) were identified using the following filters: a false discovery rate [FDR ()] <0.05 and |Log2 fold change| >1 analyzed with the limma package (). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs were further explored by the “ClusterProfiler” () and “org.Hs.eg.db” () packages. The results were visualized with the “ggplot2” () package.
Immune-Related Gene Analyses
A dataset with 1,793 unique immune-related genes (IRGs) that had been confirmed in previous studies as being involved in the natural immune process was downloaded from the IMMPORT database [(); https://www.immport.org/shared/]. The differentially expressed IRGs were the intersection results between the DEGs and the 1793 IRGs. The protein-protein interaction (PPI) network of differentially expressed IRGs was constructed using the STRING database [(); https://string-db.org/, Version:11.0], with the highest confidence level defined as an interaction score of more than 0.9. Subsequently, the top 20 hub genes with the most adjacent nodes in the PPI network were selected.
Feature IRGs Identified and Correlation Analysis
The least absolute shrinkage and selection operator (LASSO) regression arithmetic and Boruta arithmetic () were used independently to identify the feature IRGs from the 20 hub genes. The former is a regression analysis algorithm that uses regularization to improve prediction accuracy (). The latter is a random forest method that can provide a numerical estimate of feature importance (). In addition, the R “corrplot” package () was used to analyze a potential correlation among the selected feature IRGs.
Feature IRGs Verification and Diagnosis Model Constructed
The selected feature IRGs were further verified using the independent dataset GSE36961 with 106 HCM patients and 39 healthy controls. The gene expression values were displayed by boxplots, and receiver-operating characteristic (ROC) curves were used to assess the discrimination ability of the feature IRGs. A five-fold cross-validation was conducted 500 times on the GSE36961 cohort to conform the reliability of the area under the ROC curves (AUC) of the selected biomarkers. Moreover, an online webtool NetworkAnalyst () was used to further explore the potential miRNA—feature IRG regulatory network and transcription factor (TF)—feature IRG regulation mechanisms. Finally, a diagnosis model combining the selected biomarkers was developed by logistic regression and visualized as a nomogram (). The C-statistic and calibration curves were used to evaluate the identification ability and calibration of the model.
Results
DEGs Between HCM and Normal Samples
After the normalization of GSE160997 dataset (Figures 2A,B), a total of 1,079 DEGs were identified in the GSE160997 dataset, and filtering resulted in 618 significantly upregulated and 461 significantly downregulated DEGs (Figure 2C). The top 30 DEGs were displayed in a heatmap (Figure 2D). The GO annotation and KEGG pathway enrichment analyses of DEGs are illustrated in Figures 2E,F. It became obvious that inflammation and immune-related signaling pathways and biological regulatory functions played a dominant role in the enrichment results of the DEGs, which prompted us to further explore the IRGs among the DEGs.
Figure 2
The Information of Immune-Related Genes
A total of 121 differentially expressed IRGs between the DEGs and an IRG dataset from the IMMPORT database were identified (Figure 3A). These 121 differentially expressed IRGs formed a close PPI network interaction at a confidence score >0.9 (Figure 3B). The top 20 hub genes with the most adjacent nodes in the PPI network are listed in Figure 3C, and the detailed information is listed in Table 1.
Figure 3
Table 1
| Gene symbol | logFC | FDR | Adjacent nodes | Description |
|---|---|---|---|---|
| C3 | −1.99 | <0.001 | 16 | Complement C3. Also known as: AHUS5, ARMD9, ASP, C3a, C3b, CPAMD1, and HEL-S-62p |
| C3AR1 | −2.09 | <0.001 | 15 | Complement C3a receptor 1. Also known as: AZ3B, C3AR, and HNFAG09 |
| FCER1G | −2.10 | <0.001 | 14 | Fc fragment of IgE receptor Ig. Also known as: FCRG |
| ITGB2 | −1.96 | <0.001 | 14 | Integrin subunit beta 2. Also known as: CD18, LAD, LCAMB, LFA-1, MAC-1, MF17, and MFI7 |
| GPR17 | 2.32 | 0.005 | 14 | G protein-coupled receptor 17 |
| PENK | 2.27 | 0.009 | 14 | Proenkephalin. Also known as: PE and PENK-A |
| SYK | −1.41 | <0.001 | 13 | Spleen associated tyrosine kinase. Also known as: IMD82 and p72-Syk |
| FPR1 | −2.24 | <0.001 | 13 | Formyl peptide receptor 1. Also known as: FMLP and FPR |
| VAV1 | −1.60 | <0.001 | 12 | Vav guanine nucleotide exchange factor 1. Also known as: VAV |
| PTAFR | −1.25 | 0.003 | 12 | Platelet activating factor receptor. Also known as: PAFR |
| AGTR1 | −1.23 | 0.008 | 11 | Angiotensin II receptor type 1. Also known as: AG2S, AGTR1B, AT1, AT1AR, AT1B, AT1BR, AT1R, AT2R1, and HAT1R |
| C5AR1 | −1.59 | <0.001 | 11 | Complement C5a receptor 1. Also known as: C5A, C5AR, C5R1, and CD88 |
| CCR1 | −2.65 | <0.001 | 11 | C-C motif chemokine receptor 1. Also known as: CD191, CKR-1, CKR1, CMKBR1, HM145, MIP1aR, and SCYAR1 |
| SSTR2 | −1.39 | <0.001 | 11 | Somatostatin receptor 2 |
| CYBB | −1.84 | <0.001 | 10 | Cytochrome b-245 beta chain. Also known as: AMCBX2, CGD, CGDX, GP91-1, GP91-PHOX, GP91PHOX, IMD34, NOX2, and p91-PHOX |
| CCL21 | −2.18 | 0.004 | 10 | C-C motif chemokine ligand 21. Also known as: 6Ckine, CKb9, ECL, SCYA21, SLC, and TCA4 |
| EDN1 | −1.31 | 0.003 | 10 | Endothelin 1. Also known as: ARCND3, ET1, HDLCQ7, PPET1, and QME |
| ACKR3 | 1.17 | <0.001 | 10 | Atypical chemokine receptor 3. Also known as: CMKOR1, CXC-R7, CXCR-7, CXCR7, GPR159, RDC-1, and RDC1 |
| JUN | 1.18 | <0.001 | 10 | Jun proto-oncogene, AP-1 transcription factor subunit. Also known as: AP-1, AP1, c-Jun, cJUN, and p39 |
| CXCR6 | 1.73 | 0.002 | 10 | C-X-C motif chemokine receptor 6. Also known as: BONZO, CD186, CDw186, STRL33, and TYMSTR |
The detail information of the top-20 immune-related genes with most adjacent nodes in protein-protein interaction network.
Feature IRGs Identified and Correlation Analysis
The LASSO arithmetic selected 6 feature IRGs (Figures 4A,B) and the Boruta arithmetic 14 feature IRGs (Figure 4C). Five feature IRGs (ACKR3, AGTR1, CCR1, CYBB, and JUN) were identified by both methods (Figure 4D). The correlation analysis result displaying the Pearson correlation values indicated that there was a strong positive correlation between CYBB and CCR1 and a moderate positive correlation between ACKR3 and JUN (Figure 4E).
Figure 4
Feature IRG Verification and Diagnosis Model Constructed
The gene expression values of the five selected IRGs in GSE36961 are illustrated in Figure 5, and the original data are presented in the Supplementary File. It is worth noting that JUN was highly expressed in the normal controls in GSE36961, which was contrary to the previous result. The gene expressions of the remaining four IRGs (ACKR3, AGTR1, CCR1, and CYBB) were consistent with the previous analysis and were statistically significant, suggesting them as potential diagnostic biomarkers of HCM. The AUC of the four biomarkers were all remarkable, indicating their ability to contribute to the diagnosis of HCM (Figure 6). The results of 500 times five-fold cross-validation show that the average AUC scores of CYBB, CCR1, AGTR1 and ACKR3 were 0.883, 0.856, 0.842, and 0.865 respectively.
Figure 5
Figure 6
The miRNA—feature IRG regulatory network and transcription factor (TF)—feature IRG regulatory network were listed as Figure 7. The four IRGs were transformed into binary variables according to the median value of gene expression, based on which a diagnosis model was constructed by logistic regression that was visualized as a nomogram (Figure 8A). The C-index of the diagnosis model was 0.925 (95% confidence interval 0.869–0.981) with appropriate calibration (Figure 8B).
Figure 7
Figure 8
Discussion
HCM is a common inherited cardiovascular disease that occurs in one of 500 people in the general population and leads to heart failure, fatal arrhythmia, and stroke. Mutations in genes encoding proteins of the cardiomyocytic contractile apparatus are considered the main cause of HCM. Genetic tests focusing on sarcomere protein genes have been recommended in current guidelines and experience widespread adoption. However, with the accumulation of experimental data (–), approximately two-thirds of HCM patients were shown not to have sarcomere protein mutations. Therefore, it is necessary to explore alternative molecular mechanisms and identify more valuable biomarkers of HCM by utilizing contemporary methods to analyze biological complexity with an emphasis on network medicine ().
According to bioinformatics analysis, the GO annotation and KEGG enrichment analysis of the DEGs in GSE160997 were dominated by immune and inflammatory response mechanisms, such as neutrophil and macrophage activation, neutrophil degranulation, regulation of leukocyte-mediated immunity, and immune effector processes. Currently, inflammation and immunity attract substantial interest with regard to their role in the pathophysiology of cardiovascular diseases. Some previous studies (–) suggested that immune and inflammatory response mechanisms may play an important role in myocardial hypertrophy and remodeling. In contrast, most research has focused on heart failure and myocardial remodeling after acute myocardial infarction, with extremely limited evidence on HCM.
The preliminary results of bioinformatics analysis enabled further exploration of IRGs as biomarkers of HCM. A dataset with 1,793 unique IRGs was downloaded from the IMMPORT database and used to identify the IRGs among the DEGs in GSE160997. A total of 121 IGRs were selected from the 1,079 DEGs, and 5 feature IRGs were identified as candidate biomarkers by both the LASSO and Boruta arithmetic. Finally, four feature IRGs were confirmed as biomarkers with an excellent AUC score using the GSE36961 dataset of HCM patients and 39 healthy controls. Our novel diagnosis model of HCM constructed based on these four feature IRGs showed high discrimination ability (C-index = 0.925) and favorable calibration. This pioneering research offers a new perspective on the identification of potential biomarkers that might aid in the diagnosis of HCM.
Angiotensin II receptor type 1 (AGTR1) is a significant effector of the renin-angiotensin-aldosterone system, which is associated with cardiovascular diseases, especially cardiac remodeling. The system consists of two regulatory pathways: a pathological pro-inflammatory pathway and a protective anti-inflammatory pathway. AGTR1, involved in the former, can promote vasoconstriction, cell proliferation, inflammation, oxidative stress, hypertrophy, and fibrosis, and angiotensin II–type 1 receptor blockers are widely used to prevent cardiac remodeling after acute myocardial infarction and heart failure (, ).
Atypical chemokine receptor 3 (ACKR3) (), previously known as C-X-C chemokine receptor type 7 (CXCR7), has emerged as a potential hub gene in cardiac development in animal experiments (–). Nevertheless, animal studies have focused on the effects of ACKR3 knockout on cardiac development. In this study, the expression of ACKR3 in HCM patients was significantly higher than that in the normal controls in both GSE160997 and GSE36961. At present, no study has explored the effect of a high expression of ACKR3 on cardiac development or structural remodeling, which may be a promising area of research going forward.
CCR1 gene encodes a protein that acts as a chemokine receptor. This protein and its ligands macrophage inflammatory protein 1 (MIP-1, CCL3) and regulated on activation normal T expressed and secreted protein (RANTES, CCL5) can be expressed in monocytes/macrophages/T lymphocytes, whose signaling is the key to recruit effector immune cells to inflammatory sites. Liehn et al. proved that in model mice with myocardial infarction, the deletion of CCR1 can reduce tissue inflammatory injury and protect the myocardium (). Data from animal and human studies show that CCR1 and its ligands are central mediators of a variety of inflammatory diseases, making the receptor an attractive candidate for therapeutic intervention. At present, CCR1 receptor antagonists are mainly used in clinical trials for the treatment of rheumatoid arthritis (). Batis et al. selectively blocks the CCR1 receptor through met RANTES therapy, confirming that a treatment strategy based on regulating CCR1/CCR5 mediated cell migration and/or effector function may help limit cardiac tissue damage during chronic Chagas disease ().
The protein encoded by cytochrome b-245 beta chain (CYBB) gene, also named NOX-2, is a component of cytochrome b-245, which plays an important role in the development of phagocytic microbicidal oxidase system. A consensus has been reached that the deficiency of CYBB leads to decreased NADPH oxidase activity in phagocytes, and neutrophils phagocytize bacteria but do not effectively kill bacteria, and that this might be the potential pathogenesis of chronic granulomatous disease (). In addition, the gene has also been reported to be associated with ventricular hypertrophy and arrhythmia caused by hypoxia and high-fat diet (–).
The strengths of this study are (a) the innovative perspective that immune and inflammation may underlie the molecular mechanism of HCM and (b) the identification of IRGs as diagnostic biomarkers using bioinformatics analysis of transcriptome datasets of HCM patients and healthy controls. Nonetheless, some limitations of this study need to be recognized. First, the study was a retrospective data analysis and lacked detailed clinical or prognostic information, such as the thickness of the ventricular septum and the occurrence of major adverse cardiovascular events in patients with HCM. This limited the further exploration of the relevant genes determining these clinical features or outcomes. Second, although the four feature biomarkers were verified by an independent dataset of sufficient size and achieved excellent discrimination, a larger dataset is essential for verification in the future. Third, the different transcriptome datasets were detected by different sequencing platforms. The baselines of the different transcriptome datasets were inconsistent, indicating that they may only determine whether the expression of characteristic genes is abnormal without the ability to identify a generally applicable reference range.
Conclusion
The systematic analysis of the transcriptional profiles of HCM and exploration of the potential underlying molecular mechanisms led to the identification of four feature IRGs as biomarkers of HCM. Based on these, a diagnosis model of HCM was developed that brings a novel aspect to the current understanding of the pathological processes and required genetic testing in HCM.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
XZ and GL were involved in the conception and design of the study. XZ was responsible for visualization and article writing. RH provided scientific supervision. All authors reviewed and approved the final manuscript.
Acknowledgments
We would like to thank the Gene Expression Omnibus (GEO) database for the precious data used for free in scientific research.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.752559/full#supplementary-material
- HCM
hypertrophic cardiomyopathy
- GO
gene ontology
- KEGG
Kyoto Encyclopedia of Genes and Genomes
- DEGs
differentially expressed genes
- IRG
immune-related gene
- PPI
protein-protein interaction
- LASSO
least absolute shrinkage and selection operator
- ROC curve
receiver operator characteristics curve.
Abbreviations
References
1.
WolfCM. Hypertrophic cardiomyopathy: genetics and clinical perspectives. Cardiovasc Diagnosis Ther. (2019) 9:S388–415. 10.21037/cdt.2019.02.01
2.
MaronBJMaronMS. Hypertrophic cardiomyopathy. Lancet. (2013) 381:242–5510.1016/S0140-6736(12)60397-3
3.
OmmenSRMitalSBurkeMADaySMDeswalAElliottPet al. 2020 AHA/ACC guideline for the diagnosis and treatment of patients with hypertrophic cardiomyopathy: executive summary: a report of the American college of cardiology/American heart association joint committee on clinical practice guidelines. J Am Coll Cardiol. (2020) 76:3022–55
4.
UKCOM. 2014 esc guidelines on diagnosis and management of hypertrophic cardiomyopathy. Eur Heart J. (2014) 35:2733–79. 10.1093/eurheartj/ehu284
5.
KoCArscottPConcannonMSaberiSDaySMYasharBMet al. Genetic testing impacts the utility of prospective familial screening in hypertrophic cardiomyopathy through identification of a nonfamilial subgroup. Genet Med. (2018) 20:69–75. 10.1038/gim.2017.79
6.
MaronBAWangR-SShevtsovSDrakosSGAronsEWever-PinzonOet al. Individualized interactomes for network-based precision medicine in hypertrophic cardiomyopathy with implications for other clinical pathophenotypes. Nat Commun. (2021) 12:873. 10.1038/s41467-021-21146-y
7.
GuoQWangJSunRGuWHeZChenQet al. Identification of circulating hub long noncoding RNAs associated with hypertrophic cardiomyopathy using weighted correlation network analysis. Mol Med Rep. (2020) 22:4637–44. 10.3892/mmr.2020.11566
8.
GuoQWangJSunRHeZChenQLiuWet al. Comprehensive construction of a circular RNA-associated competing endogenous RNA network identified novel circular RNAs in hypertrophic cardiomyopathy by integrated analysis. Front Genet. (2020) 11:764. 10.3389/fgene.2020.00764
9.
LiuXMaYYinKLiWChenWZhangYet al. Long non-coding and coding RNA profiling using strand-specific RNA-seq in human hypertrophic cardiomyopathy. Sci Data. (2019) 6:90. 10.1038/s41597-019-0094-6
10.
ZhengXYangYFuCHHuangR. Identification and verification of promising diagnostic biomarkers in patients with hypertrophic cardiomyopathy associate with immune cell infiltration characteristics. Life Sci. (2021) 285:119956. 10.1016/j.lfs.2021.119956
11.
KorthauerKKimesPKDuvalletCReyesASubramanianATengMet al. A practical guide to methods controlling false discoveries in computational biology. Genome Biol. (2019) 20:118. 10.1186/s13059-019-1716-1
12.
RitchieMEPhipsonBWuDHuYLawCWShiWet al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. 10.1093/nar/gkv007
13.
YuGWangL-GHanYHeQ-Y. Clusterprofiler: an r package for comparing biological themes among gene clusters. Omics. (2012) 16:284–7. 10.1089/omi.2011.0118
14.
CarlsonMFalconSPagesHLiN. org.Hs.eg.db: Genome Wide Annotation for Human. R package version (2019). p. 3.
15.
GinestetC. Ggplot2: elegant graphics for data analysis. J R Statist Soc Ser A. (2011) 174:245. 10.1111/j.1467-985X.2010.00676_9.x
16.
BhattacharyaSAndorfSGomesLDunnPSchaeferHPontiusJet al. ImmPort: disseminating data to the public for the future of immunology. Immunol Res. (2014) 58:234–9. 10.1007/s12026-014-8516-1
17.
MeringCVHuynenMJaeggiDSchmidtSBorkPSnelB. String: a database of predicted functional associations between proteins. Nucleic Acids Res. (2003) 31:258–61. 10.1093/nar/gkg034
18.
KursaMBRudnickiWR. Feature selection with the Boruta package. J Stat Softw. (2010) 36:1–13. 10.18637/jss.v036.i11
19.
ZhaoEXieHZhangY. Predicting diagnostic gene biomarkers associated with immune infiltration in patients with acute myocardial infarction. Front Cardiovasc Med. (2020) 7:586871. 10.3389/fcvm.2020.586871
20.
ZhengXHuangRLiuGJiaZChenKHeY. Development and verification of a predictive nomogram to evaluate the risk of complicating ventricular tachyarrhythmia after acute myocardial infarction during hospitalization: a retrospective analysis. Am J Emerg Med. (2020) 46:462–8. 10.1016/j.ajem.2020.10.052
21.
WeiTSimkoVLevyMXieYJinYZemlaJ. Corrplot: Visualization of a Correlation Matrix. R package version 0.73. (2013). p. 230.
22.
ZhouGSoufanOEwaldJHancockREBasuNXiaJ. Networkanalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. (2019) 47:W234–41. 10.1093/nar/gkz240
23.
DengMYinYZhangQZhouXHouG. Identification of inflammation-related biomarker Lp-PLA2for patients with COPD by comprehensive analysis. Front Immunol. (2021) 12:670971. 10.3389/fimmu.2021.670971
24.
NeubauerSKolmPHoCYKwongRYDesaiMYDolmanSFet al. Distinct subgroups in hypertrophic cardiomyopathy in the NHLBI HCM registry. J Am Coll Cardiol. (2019) 74:2333–45. 10.1016/j.jacc.2019.08.1057
25.
BosJMWillMLGershBJKruisselbrinkTMOmmenSRAckermanMJ. Characterization of a phenotype-based genetic test prediction score for unrelated patients with hypertrophic cardiomyopathy. Mayo Clinic Proc. (2014) 89:727–37. 10.1016/j.mayocp.2014.01.025
26.
GrunerCIvanovJCareMWilliamsLMoravskyGYangHet al. Toronto hypertrophic cardiomyopathy genotype score for prediction of a positive genotype in hypertrophic cardiomyopathy. Circ Cardiovasc Genet. (2013) 6:19–26. 10.1161/CIRCGENETICS.112.963363
27.
MaronBJMaronMSMaronBALoscalzoJ. Moving beyond the sarcomere to explain heterogeneity in hypertrophic cardiomyopathy: JACC review topic of the week. J Am Coll Cardiol. (2019) 73:1978–86. 10.1016/j.jacc.2019.01.061
28.
ShimadaYJRaitaYLiangLWMaurerMSHasegawaKFiferMAet al. Comprehensive proteomics profiling reveals circulating biomarkers of hypertrophic cardiomyopathy. Circ Heart Fail. (2021) 14:e007849. 10.1161/CIRCHEARTFAILURE.120.007849
29.
ChenRGeTJiangWHuoJChangQGengJet al. Identification of biomarkers correlated with hypertrophic cardiomyopathy with co-expression analysis. J Cell Physiol. (2019) 234:21999–2008. 10.1002/jcp.28762
30.
BeckerRCOwensAPSadayappanS. Tissue-level inflammation and ventricular remodeling in hypertrophic cardiomyopathy. J Thromb Thrombolysis. (2020) 49:177–83. 10.1007/s11239-019-02026-1
31.
KuusistoJKärjäVSipolaPKholováIPeuhkurinenKJääskeläinenPet al. Low-grade inflammation and the phenotypic expression of myocardial fibrosis in hypertrophic cardiomyopathy. Heart. (2012) 98:1007–13. 10.1136/heartjnl-2011-300960
32.
SaraviBLiZLangCNSchmidBLangFKGradSet al. The tissue renin-angiotensin system and its role in the pathogenesis of major human diseases: quo vadis?Cells. (2021) 10:650. 10.3390/cells10030650
33.
De MelloWC. Local renin angiotensin aldosterone systems and cardiovascular diseases. Med Clinics. (2017) 101:117–27. 10.1016/j.mcna.2016.08.017
34.
KoenenJBachelerieFBalabanianKSchlecht-LoufGGallegoC. Atypical chemokine receptor 3 (ackr3): a comprehensive overview of its expression and potential roles in the immune system. Mol Pharmacol. (2019) 96:809–18. 10.1124/mol.118.115329
35.
YuSCrawfordDTsuchihashiTBehrensTWSrivastavaD. The chemokine receptor cxcr7 functions to regulate cardiac valve remodeling. Dev Dyn. (2011) 240:384–93. 10.1002/dvdy.22549
36.
GerritsHvan Ingen SchenauDSBakkerNEvan DisseldorpAJStrikAHermensLSet al. Early postnatal lethality and cardiovascular defects in cxcr7-deficient mice. Genesis. (2008) 46:235–45. 10.1002/dvg.20387
37.
KleinKRKarpinichNOEspenschiedSTWillcocksonHHDunworthWPHoopesSLet al. Decoy receptor cxcr7 modulates adrenomedullin-mediated cardiac and lymphatic vascular development. Dev Cell. (2014) 30:528–40. 10.1016/j.devcel.2014.07.012
38.
LiehnEMerxMPosteaOBecherSDjalali-TalabYShagdarsurenEet al. Ccr1 deficiency reduces inflammatory remodelling and preserves left ventricular function after myocardial infarction. J Cell Mol Med. (2008) 12:496–506. 10.1111/j.1582-4934.2007.00194.x
39.
WhiteGEIqbalAJGreavesDR. Cc chemokine receptors and chronic inflammation—therapeutic opportunities and pharmacological challenges. Pharmacol Rev. (2013) 65:47–89. 10.1124/pr.111.005074
40.
BatistaAMAlvarado-ArnezLEAlvesSMMeloGPereiraIRRuivoLADSet al. Genetic polymorphism at CCL5 is associated with protection in Chagas' heart disease: antagonistic participation of CCR1+ and CCR5+ cells in chronic chagasic cardiomyopathy. Front Immunol. (2018) 9:615. 10.3389/fimmu.2018.00615
41.
AnjaniGVigneshPJoshiVShandilyaJKBhattaraiDSharmaJet al. Recent advances in chronic granulomatous disease. Genes Dis. (2020) 7:84–92. 10.1016/j.gendis.2019.07.010
42.
PenaESiquesPBritoJArribasSMBögerRHannemannJet al. Nox2 upregulation and p38α MAPK activation in right ventricular hypertrophy of rats exposed to long-term chronic intermittent hypobaric hypoxia. Int J Mol Sci. (2020) 21:8576. 10.3390/ijms21228576
43.
BhattiSNLiJ-M. Nox2 dependent redox-regulation of AKT and ERK1/2 to promote left ventricular hypertrophy in dietary obesity of mice. Biochem Biophys Res Commun. (2020) 528:506–13. 10.1016/j.bbrc.2020.05.162
44.
JosephLCAvulaUMRWanEYReyesMVLakkadiKRSubramanyamPet al. Dietary saturated fat promotes arrhythmia by activating NOX2 (NADPH oxidase 2). Circulation Arrhythm Electrophysiol. (2019) 12:e007573. 10.1161/CIRCEP.119.007573
Summary
Keywords
biomarkers, hypertrophic cardiomyopathy, immune-related gene, diagnosis model, random forest—ensemble classifier, LASSO
Citation
Zheng X, Liu G and Huang R (2021) Identification and Verification of Feature Immune-Related Genes in Patients With Hypertrophic Cardiomyopathy Based on Bioinformatics Analyses. Front. Cardiovasc. Med. 8:752559. doi: 10.3389/fcvm.2021.752559
Received
03 August 2021
Accepted
28 September 2021
Published
26 October 2021
Volume
8 - 2021
Edited by
Alexandre Francois Roy Stewart, University of Ottawa, Canada
Reviewed by
Basavaraj Mallikarjunayya Vastrad, Soniya Education Trust's College of Pharmacy, India; Thiago Britto-Borges, Heidelberg University Hospital, Germany
Updates
Copyright
© 2021 Zheng, Liu and Huang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ruina Huang wolidedd@gdmu.edu.cn
This article was submitted to Cardiovascular Genetics and Systems Medicine, a section of the journal Frontiers in Cardiovascular Medicine
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.