Differentially expressed genes in hypertrophic cardiomyopathy
The difference of gene expression between patients in disease state and healthy people may be closely related to the disease. We analyzed the data of GSE130036 in geo database to explore the genes related to hypertrophic cardiomyopathy. By setting the screening threshold twofold, we identified 920 differentially expressed genes (DEGs) between the hypertrophic cardiomyopathy and the healthy controls (Table S1). There are 636 up-regulated genes with high expression level and 284 down regulated genes with low expression level (Fig. 1A, 1B). It is important to compare these differentially expressed genes with those of HCM and control groups in gse36961 data. We found 18 common differentially expressed genes (Table 1). We suggest that the differentially expressed genes between HCM and control may be dysfunctional genes related to hypertrophic cardiomyopathy.
Table 1
The common differentially expressed genes in the GSE130036 and GSE36961.
| GSE130036 | | GSE36961 | | |
symbol | log2FoldChange | pvalue | padj | logFC | P.Value | adj.P.Val |
CENPA | 2.428995 | 1.71E-06 | 2.99E-05 | 1.699106 | 2.22E-18 | 1.62E-16 |
TUBA3E | -3.20702 | 8.72E-23 | 4.96E-20 | -2.17305 | 3.12E-47 | 1.18E-42 |
TUBA3D | -3.1646 | 6.82E-26 | 6.25E-23 | -2.20486 | 2.59E-44 | 1.63E-40 |
CORIN | -3.2704 | 8.28E-13 | 6.96E-11 | -1.56596 | 5.15E-11 | 1.18E-09 |
SLITRK4 | 3.644176 | 1.91E-19 | 6.24E-17 | 1.260442 | 1.38E-14 | 5.43E-13 |
CA3 | 2.876428 | 5.99E-08 | 1.54E-06 | 1.362298 | 9.87E-13 | 2.93E-11 |
LYVE1 | -2.0283 | 9.92E-28 | 1.06E-24 | -1.87319 | 2.62E-39 | 7.08E-36 |
METTL7B | -2.17555 | 6.86E-20 | 2.45E-17 | -1.35577 | 6E-23 | 9.23E-21 |
SERPINA3 | -3.81907 | 4.52E-09 | 1.53E-07 | -3.56141 | 3.7E-42 | 1.75E-38 |
RASD1 | -2.08726 | 6.49E-12 | 4.42E-10 | -3.51434 | 2.45E-45 | 1.85E-41 |
SLCO4A1 | -2.59057 | 1.24E-24 | 9.11E-22 | -1.00128 | 1.37E-32 | 1.2E-29 |
COMP | 2.861793 | 8.22E-07 | 1.56E-05 | 1.099086 | 7.8E-09 | 1.24E-07 |
C1QB | -2.15322 | 1.30E-30 | 2.13E-27 | -1.46863 | 4.93E-28 | 1.94E-25 |
F13A1 | -2.29595 | 1.12E-38 | 7.01E-35 | -1.5787 | 3.76E-30 | 2.26E-27 |
CD163 | -2.76903 | 4.77E-33 | 1.19E-29 | -2.13651 | 1.05E-37 | 2.34E-34 |
FCN3 | -2.50434 | 2.72E-17 | 6.00E-15 | -2.03805 | 4.98E-40 | 1.57E-36 |
PLA2G2A | -2.93842 | 3.35E-08 | 9.14E-07 | -1.35137 | 2.51E-20 | 2.58E-18 |
CHRDL2 | -3.62053 | 2.89E-10 | 1.33E-08 | -1.17134 | 1.47E-29 | 7.62E-27 |
Mechanism of dysregulation related to differentially expressed genes
In order to explore the biological functions of these dysfunctional genes in HCM, we analyzed the enrichment of DEG with GO and KEGG. From the results of biological process (BP) enrichment, we found that the maladjusted genes are mainly involved in the biological process related to epidermis development (Figure 2A, 2B). In addition, in the enrichment results of cell component (CC), the maladjusted gene is mainly related to hemoglobin complex (Figure 2C, 2D). In the molecular function (MF), the maladjusted gene is mainly enriched in the molecular activity related molecular function (Figure 2E, 2F). On the other hand, we found that the signal pathways involved in the maladjusted genes mainly include arachidonic acid metabolism, aldosterone synthesis and secret, and drug metabolism (Figure 2G, 2H). The results showed that the maladjusted genes were mainly related to metabolism related signaling pathway in the course of HCM.
PPI network identification of key dysregulated molecules
In order to provide the contents and ways of DEG participating in cell biological activities, we constructed a full view of their interacting proteins to elucidate their functional networks. After mapping up-regulated and down-regulated genes into the network, the PPI network of 187 dysregulated genes (Figure 3A) was screened. Currently, we map 18 common differentially expressed genes into PPI network. Six genes were screened, including C1QB, F13A1, CD163, FCN3, PLA2G2A and CHRDL2. Their expression in the two groups of data is similar, with consistent down-regulation behavior (Figure 3B). In addition, ROC analysis showed that the AUC values of the six genes were all over 91%, indicating that they had a certain clinical diagnostic ability (Figure 3C). Therefore, we think these six genes may be potential biological target genes of HCM.
Regulatory mechanism of HCM influenced by mutation
Using gene cohorts in GSE130036 data, we identified mutations in key genes generated using maftools (Table S2, figure 4A). The SNVs in the sample can be divided into six categories (Figure 4B). Thus, we present the mutation types and frequencies of the first 10 mutations (Figure 4C). Including PDE11A, PRMT5−AS1, TSPAN9, MTR−RPL35P1, RBM23, NFKBIZ, PRMT5, RP11−14N7.2, NOTCH2 and RP5−857K21.4. In addition, we found that the total number of samples with PRMT5 mutation was the largest (Figure 4D). It is suggested that the mutation frequency of HCM is the highest, which may have a wider impact. The above results show that the phenomenon of gene mutation often occurs in HCM patients, which means that gene mutation has a significant impact on the disease.