Figures & data
Figure 1 Workflow chart of the study.
![Figure 1 Workflow chart of the study.](/cms/asset/49576978-8f0c-488b-94d6-183cb0583206/dijg_a_12168273_f0001_b.jpg)
Table 1 Homologous Between CNT Exposure and Controls of Human and Mouse
Figure 2 Identification of DEGs between CNT-exposed and control human pleural mesothelial cells (MeT5A) in the GSE48855 dataset. (A), the heatmap of DEGs identified by comparing SWCNTs-exposed with dispersant-treated MeT5A cells; (B), the heatmap of DEGs identified by comparing MWCNTs-exposed with dispersant-treated MeT5A cells; (C), Venn diagram to screen the common upregulated genes between SWCNTs and MWCNTs; (D), Venn diagram to the common downregulated genes between SWCNTs and MWCNTs. The heatmap was generated based on the Z-score.
![Figure 2 Identification of DEGs between CNT-exposed and control human pleural mesothelial cells (MeT5A) in the GSE48855 dataset. (A), the heatmap of DEGs identified by comparing SWCNTs-exposed with dispersant-treated MeT5A cells; (B), the heatmap of DEGs identified by comparing MWCNTs-exposed with dispersant-treated MeT5A cells; (C), Venn diagram to screen the common upregulated genes between SWCNTs and MWCNTs; (D), Venn diagram to the common downregulated genes between SWCNTs and MWCNTs. The heatmap was generated based on the Z-score.](/cms/asset/52c0ccc9-b47e-4b1c-8326-559d3d8fcb04/dijg_a_12168273_f0002_c.jpg)
Table 2 Modules Based on WGCNA of Genes in the GSE48855 Dataset
Figure 3 Identification of key modules associated with CNT-induced malignant transformation of MeT5A cells by weighted gene co-expression network analysis. (A), calculation of the scale-free topology fit index according to different soft-threshold power values; (B), calculation of mean connectivity according to different soft-threshold power values; (C), a clustering dendrogram of samples; (D), the correlations between modules and disease status (r, correlation coefficient); (E), scatter plot of module eigengenes in the blue module; (F), scatter plot of module eigengenes in the brown module.
![Figure 3 Identification of key modules associated with CNT-induced malignant transformation of MeT5A cells by weighted gene co-expression network analysis. (A), calculation of the scale-free topology fit index according to different soft-threshold power values; (B), calculation of mean connectivity according to different soft-threshold power values; (C), a clustering dendrogram of samples; (D), the correlations between modules and disease status (r, correlation coefficient); (E), scatter plot of module eigengenes in the blue module; (F), scatter plot of module eigengenes in the brown module.](/cms/asset/c2822b45-cdd0-4fcd-ad8e-337afc0a5ff9/dijg_a_12168273_f0003_c.jpg)
Figure 4 The protein–protein interaction network established by interaction pairs of four hub genes. Genes in circular represent hub genes. Gene color is corresponding to module color; regular triangle, upregulated genes; inverted triangle, downregulated genes.
![Figure 4 The protein–protein interaction network established by interaction pairs of four hub genes. Genes in circular represent hub genes. Gene color is corresponding to module color; regular triangle, upregulated genes; inverted triangle, downregulated genes.](/cms/asset/a913d1c0-4698-4cd5-a0cb-19e82541d6bf/dijg_a_12168273_f0004_c.jpg)
Figure 5 Kaplan-Meier survival curves for patients with a high and low expression of five hub module genes in mesothelioma samples. (A–E), GEPIA analysis results; F-J, TIMER analysis results. (A and F) for FN1; (B and G) for ALCAM; (C and H) for THBS2; (D and I) for UGCG; (E and J) for CHPF2.
![Figure 5 Kaplan-Meier survival curves for patients with a high and low expression of five hub module genes in mesothelioma samples. (A–E), GEPIA analysis results; F-J, TIMER analysis results. (A and F) for FN1; (B and G) for ALCAM; (C and H) for THBS2; (D and I) for UGCG; (E and J) for CHPF2.](/cms/asset/7abdf84f-f109-4357-9713-c56fcde0683a/dijg_a_12168273_f0005_c.jpg)
Table 3 Validation of the Expression of Hub Genes
Table 4 The Predictive Power of Hub Genes for MPM
Figure 6 ROC curves for evaluating the diagnostic values of four hub genes for distinguishing mesothelioma from normal patients. (A and B), ROC curves were plotted based on the mRNA expression levels of four genes in the GSE2549 dataset. (A), single gene; (B), gene-combination; (C and D), ROC curves were plotted based on the mRNA expression levels of four genes in the GSE51024 dataset. (C), single gene; (D), gene-combination; (E), ROC curve was plotted based on the mRNA expression levels of two genes in all mesothelioma samples of the GSE42977 dataset; (F), ROC curve was plotted based on the mRNA expression levels of three genes in the sarcomatoid and biphasic histological subtypes of the GSE42977 dataset.
![Figure 6 ROC curves for evaluating the diagnostic values of four hub genes for distinguishing mesothelioma from normal patients. (A and B), ROC curves were plotted based on the mRNA expression levels of four genes in the GSE2549 dataset. (A), single gene; (B), gene-combination; (C and D), ROC curves were plotted based on the mRNA expression levels of four genes in the GSE51024 dataset. (C), single gene; (D), gene-combination; (E), ROC curve was plotted based on the mRNA expression levels of two genes in all mesothelioma samples of the GSE42977 dataset; (F), ROC curve was plotted based on the mRNA expression levels of three genes in the sarcomatoid and biphasic histological subtypes of the GSE42977 dataset.](/cms/asset/d368c6ac-aed6-49d6-97c4-3ee56226c5d6/dijg_a_12168273_f0006_c.jpg)
Table 5 Correlation Analysis Between Hub Genes and Markers of Immune Cells
Figure 7 Correlations between the expression of hub genes and the abundance of six immune infiltration cells of mesothelioma. (A), FN1; (B), UGCG; (C), CHPF2; (D), THBS2.
![Figure 7 Correlations between the expression of hub genes and the abundance of six immune infiltration cells of mesothelioma. (A), FN1; (B), UGCG; (C), CHPF2; (D), THBS2.](/cms/asset/b612cb34-1828-4772-ae82-6f1de475c673/dijg_a_12168273_f0007_c.jpg)