Abstract:This study aimed to screen differentially expressed genes (DEGs) involved in the influence of antiangiogenic therapy on myeloid-derived suppressor cell (MDSC) infiltration and investigate their mechanisms of action. Data on DEGs after the action of antiangiogenic drugs in a pan-cancer context were obtained from the Gene Expression Omnibus (GEO) database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the clusterProfiler package in R software. Single-sample gene set enrichment analysis was performed using the gene set variation analysis package to evaluate the levels of immune cells and the activity of immune-related pathways. The relationships of DEGs with the infiltration levels of MDSCs and specific immune cell subpopulations were investigated via gene module analysis. The top 10 key genes were subsequently obtained from PPI network analysis using the cytoHubba plugin of the Cytoscape platform. When the DEGs of the four datasets were intersected, a DEG in the intersection of three datasets and 12 DEGs in the intersection of two datasets were upregulated, and 28 DEGs in the intersection of two datasets were downregulated. GO and KEGG pathway enrichment analyses revealed that the DEGs were associated with multiple important signaling pathways closely related to tumor onset and development, including cell differentiation, cell proliferation, the cell cycle, and immune responses. Most downregulated genes in lung adenocarcinoma (LUAD) were positively correlated with MDSC expression. Only MGP was negatively correlated; the correlation between CACNG6 and MDSC expression was statistically insignificant. In lung squamous cell carcinoma (LUSC), the relationships of PMEPA1 , PCDH7 , NEURL1B , and CACNG6 with MDSC expression were statistically insignificant; MGP was negatively correlated with MDSC expression. The top 10 key genes with the highest degree scores obtained using the cytoHubba plugin of Cytoscape were AURKB , RRM2 , BUB1 , NUSAP1 , PRC1 , TOP2A , NCAPH , CENPA , KIF2C , and CCNA2 . Most of these genes were upregulated in LUAD and associated with immune cell infiltration and prognosis in tumors. An analysis of the relationships between DEGs and infiltration by other specific immune cells revealed the presence of consistent patterns in the downregulated genes, which exhibited positive correlations with the levels of Th2 cells, γδ T cells, and CD56dim NK cells, and negative correlations with other infiltrating immune cells. Antiangiogenic therapy may regulate MDSC infiltration through multiple important signaling pathways closely associated with tumor onset and development, such as cell differentiation, cell proliferation, the cell cycle, and immune responses. Antiangiogenic drugs may exert effects by affecting various types of infiltrating cells associated with immune suppression.

Integrated analysis on mRNA microarray and microRNA microarray to screen immune-related biomarkers and pathways in myelodysplastic syndrome

An Integrated Regulatory Network Based on Comprehensive Analysis of mRNA Expression, Gene Methylation and Expression of Long Non-coding RNAs (lncRNAs) in Myelodysplastic Syndromes

Identification of Circular RNAs As Novel Biomarkers and Potentially Functional Competing Endogenous RNA Network for Myelodysplastic Syndrome Patients

Network-based analysis of the molecular mechanisms of multiple myeloma and monoclonal gammopathy of undetermined significance.

Gene Expression Profiling of the Bone Marrow Mononuclear Cells from Patients with Myelodysplastic Syndrome

Soluble PD-L1 as a novel biomarker predicts poor outcomes and disease progression in de novo myelodysplastic syndromes

Bioinformatics-based screening and analysis of the key genes involved in the influence of antiangiogenesis on myeloid-derived suppressor cells and their effects on the immune microenvironment

Biological Analysis of Chronic Lymphocytic Leukemia: Integration of Mrna and Microrna Expression Profiles.

MiR-181a-2-3p as a potential diagnostic and prognostic marker for myelodysplastic syndrome

Increased Expression Of Circulating Mir-22 In The Peripheral Blood Of Patients With Myelodysplastic Syndromes Was Associated With Poorer Prognosis

Identification of potential therapeutic target genes, key miRNAs and mechanisms in acute myeloid leukemia based on bioinformatics analysis

Analyzing Differences in Hematological and Immunological Characteristics Related to Common Gene Mutations in Myelodysplastic Syndromes

Investigating the Microrna-Mrna Regulatory Network in Acute Myeloid Leukemia

The yin-yang of immunity: Immune dysregulation in myelodysplastic syndrome with different risk stratification

Identification of Key Genes and Pathways in Myeloma side population cells by Bioinformatics Analysis

Identification of biomarkers in multiple myeloma: A comprehensive study combining microarray analysis and Mendelian randomization

Identification of a novel monocyte/macrophage-related gene signature for predicting survival and immune response in acute myeloid leukemia

TIM3/CEACAM1 pathway involves in myeloid-derived suppressor cells induced CD8+ T cells exhaustion and bone marrow inflammatory microenvironment in myelodysplastic syndrome

Integrated mRNA and miRNA profiling revealed deregulation of cellular stress response in bone marrow mesenchymal stem cells derived from patients with immune thrombocytopenia

Microarray data analysis on gene and miRNA expression to identify biomarkers in non-small cell lung cancer

Hypomethylation of MIR-378 5'-flanking region predicts poor survival in young patients with myelodysplastic syndrome