Page 71 - 南京医科大学自然版
P. 71

第44卷第6期       马晴晴,顾圣玮,王红玉,等. 基于线粒体自噬相关基因的哮喘亚型鉴定及预测模型构建研究[J].
                  2024年6月                     南京医科大学学报(自然科学版),2024,44(6):802-811                        ·807 ·


                ATG5、ATG12、MAP1LC3B、CSNK2A1、TOMM20、               三聚化、亚硝基化、金属离子转运以及对缺氧反应
                UBA52、TOMM22、TOMM5 及 FUNDC1 表达显著高                 的相关功能(图5B)。
                于亚型 2(图 4D),应用预测模型对亚型 1 的预测准                          KEGG 分析结果示,两种亚型均与 IL⁃17(由
                确性更高(图4E);年龄方面两亚型间差异无统计学                          Th17 细胞和 3 型固有淋巴细胞产生的 IL⁃17 可刺激
                意义(图4F),但需指出的是亚型1中女性患者数量                          中性粒细胞性气道炎症)信号通路相关,但每种亚
                略多于男性患者(图4G)。                                     型均有各自独特的通路;亚型 1 特异富集病毒蛋白
                2.5  两种亚型GO及KEGG富集分析                              与细胞因子及受体相互作用、金黄色葡萄球菌感
                    对筛选出的DEG同时进行GO及KEGG富集分                        染、趋化因子信号等通路,亚型2特异富集与卡波西
                析。GO 分析结果示,两亚型间肽酶活性调控的功                           肉瘤相关的疱疹病毒感染、结核病及精氨酸生物合
                能特征差异显著,尽管存在一定共性,但在特定的                            成相关通路(图5C)。
                生物过程中差异明显;亚型 1 显著富集于与体液免                          2.6  小分子药物预测
                疫反应相关功能(图5A),亚型2则富集于与蛋白质                              使用WGCNA识别与哮喘不同亚型相关的关键

                A                        B                                             E
                  Consensus matrix k=2                           Probability  2 Probability
                                                                 Sex                                  ***
                                                                              0.8          1.2
                                                                 Age
                                    1                                     1
                                    2                            Cluster                   0.9
                                                                 SRC          0.2
                                                                 CSNK2B   0
                                                                            Sex            0.6
                                                                 MAP1LC3A     Feale       Probability
                                                                 ULK1     -1  Male
                                                                 MFN1       Age            0.3
                                                                 UBB      -2  70
                                                                 UBC                        0
                                                                 CSNK2A2                        Cluster 1  Cluster 2
                                                                 PINK1        20               (n=157) (n=132)
                                                                 PGAM5      Cluster
                C             ***                                SQSTM1       1 2      F   70
                  GSVA score of MRG geneset  5.0                 MAP1LC3B                  (years)  Age  50
                                                                 TOMM22
                                                                 ATG5
                                                                                           60
                                                                 UBA52
                                                                 FUNDC1
                                                                                           40
                                                                 CSNK2A1
                   4.5
                                                                 TOMM20
                                                                                           30
                                                                 MFN2
                                                                                           20
                   4.0
                                                                                                        Cluster 2
                                                                                                Cluster 1
                                                                 RPS27A
                        Cluster 1
                                Cluster 2
                                                                 TOMM5
                        (n=157) (n=132)                          ATG12                          (n=91)  (n=70)
                D                Cluster 1(n=157) Cluster 2(n=132)    G            P=0.17      P=0.81
                                                                          100
                             ****
                   15.0                                                    90
                                                         ****
                           ****                                            80      43%          51%
                   12.5                    **          ****     ****     (%)  70
                  Expression  10.0  ***  ****  * ****  ***  ***  ****  ****  Percent weight  50            Male
                                                                           60
                                            **
                                                                                                           Female
                                                                           40
                    7.5  *                                                 30      57%          49%
                                                                           20
                                                                           10
                    5.0
                                                                           0
                              CSNK2B
                         RPS27A
                     PGAM5 SRC  UBB MFN2 SQSTM1 ATG5 MAP1LC3B  UBC PINK1 ULK1 MFN1 TOMM20 UBA52  TOMM5  (n=91)  (n=70)
                                                         TOMM22
                                                  CSNK2A1
                                                               FUNDC1
                                                          MAP1LC3A
                                        ATG12
                                    CSNK2A2
                                                                                  Cluster 1
                                                                                                Cluster 2
                                            Gene
                   A:Consensus clustering matrix. B:The heatmap of the expression pattern of 22 MRGs in the two subtypes with the addition of clinical information.
                C:The boxplot of the overall MRGs expression in the two subtypes(cluster 1 and cluster 2). D:The violin plot of 22 MRG expression in the two subtypes
               (cluster 1 and cluster 2). E:The comparison of diagnostic performance of prediction models for the two subtypes(cluster 1 and cluster 2). F:Compari⁃
                son of the age between asthma patients with alternative subtype(cluster 1 and cluster 2). G:The comparison of sex distribution between asthma patients
                                             ***
                                *
                                      **
                with alternative subtype. P < 0.05,P < 0.01, P < 0.001,and  **** P < 0.000 1.
                                               图4   哮喘亚型的鉴定及亚型间的差异分析
                              Figure 4 Identification of the subtypes of asthma and differences between the subtypes
   66   67   68   69   70   71   72   73   74   75   76