Research Paper Volume 16, Issue 6 pp 5207—5223

Figure 1. The dataset GSE171821 related to MIR was downloaded from the GEO database, and the data were subjected to quantile normalization (A, B). The volcano plot of visual grouping of DEGs in the dataset GSE171821 was constructed using the ggplot2 package of R software (C), and the cluster analysis heat map of DEGs was plotted using the pheatmap package of R software (D). Similarly, the dataset GSE108940 was subjected to quantile normalization (E, F). The volcano plot of visual grouping of DEGs in the dataset GSE108940 was constructed using the ggplot2 package of R software (G), and the cluster analysis heat map of DEGs was plotted using the pheatmap package of R software (H). The diagrams of up-regulated GO pathways (I, J) and down-regulated pathways (K, L) of DEGs were plotted using the R language. Then the DEGs were subjected to KEGG enrichment analysis, and the KEGG pathway diagram was plotted (M). The DEGs were imported into the STRING database to obtain the PPI network (N). Then the PPI network was imported into Cytoscape software, and the target genes with a score <10 were obtained using the plugin cytoHubba (O). GSEA was conducted on all genes using the clusterProfiler package of R software, and the GSEA pathway enrichment analysis diagram was also plotted (P).