Construction of a Macrophage Infiltration Regulatory Network and Related Prognostic Model of High-Grade Serous Ovarian Cancer

Background High-grade serous ovarian cancer (HGSOC) carries the highest mortality in the gynecological cancers; however, therapeutic outcomes have not significantly improved in recent decades. Macrophages play an essential role in the occurrence and development of ovarian cancer, so the mechanisms of macrophage infiltration should be elucidated. Method We downloaded transcriptome data of ovarian cancers from the Gene Expression Omnibus and The Cancer Genome Atlas. After rigorous screening, 1566 HGSOC were used for data analysis. CIBERSORT was used to estimate the level of macrophage infiltration and WGCNA was used to identify macrophage-related modules. We constructed a macrophage-related prognostic model using machine learning LASSO algorithm and verified it using multiple HGSOC cohorts. Results In the GPL570-OV cohort, high infiltration level of M1 macrophages was associated with a good outcome, while high infiltration level of M2 macrophages was associated with poor outcomes. We used WGCNA to select genes correlated with macrophage infiltration. These genes were used to construct protein-protein interaction maps of macrophage infiltration. IFL44L, RSAD2, IFIT3, MX1, IFIH1, IFI44, and ISG15 were the hub genes in the network. We then constructed a macrophage-related prognostic model composed of CD38, ACE2, BATF2, HLA-DOB, and WARS. The model had the ability to predict the overall survival rate of HGSOC patients in GPL570-OV, GPL6480-OV, TCGA-OV, GSE50088, and GSE26712. In exploring the immune microenvironment, we found that CD4 memory T cells and activated mast cells showed that the degree of infiltration was higher in the high-risk group, while M1 macrophages were the opposite, and HLA molecules were overexpressed in the high-risk group. Conclusion We constructed a macrophage infiltration-related protein interaction network that provides a basis for studying macrophages in HGSOC. Our macrophage-related prognostic model is robust and widely applicable. It predicts overall survival in HGSOC patients and may improve HGSOC treatment.


Introduction
Ovarian cancer is a highly malignant gynecological tumor usually found in the advanced stage. Global cancer data statistics in 2020 showed that there were 313959 new cases of ovarian cancer and 207252 new ovarian cancer deaths [1]. High-grade serous ovarian cancer (HGSOC) is the most common subtype of ovarian cancer, accounting for 75% of ovarian cancer and 70% of deaths [2]. HGSOC has a specific genetic susceptibility, and about 15-20% of HGSOC patients show BRCA1 or BRCA2 germline mutations [3]. HGSOC frequently acquires and loses DNA, making chromosomes in these cancers unstable and increasing the risk of acquired chemotherapy resistance [4]. Platinum chemotherapeutic drugs are the first-line treatment of HGSOC, and immunerelated therapy is a new treatment modality [5]. However, the long-term survival from ovarian cancer has not significantly increased in the past three decades. Currently, prognostic factors of HGSOC include the FIGO stage, residual disease, BRCA1/2 germline mutations, and tumorinfiltrating lymphocyte score. However, these prognostic factors have great limitations, and the prediction effect is not ideal. erefore, the construction of a molecular prognostic model of ovarian cancer is the focus of future research [5].
Macrophages play a complex and important role in tumor microenvironment. Resting macrophage are polarized into proinflammatory or anti-inflammatory forms [6]. M1 macrophages have properties of cytotoxicity, tumor inhibition, and immunostimulation functions, while M2 macrophages promote tumor growth and invasion [7]. Tumor-associated macrophages (TAM) are macrophages infiltrating around tumor cells, which are closely related to the malignant progression and clinical prognosis of ovarian cancer. In recent years, research on the targeting strategy of TAM has achieved great success. TAM targeting strategies include inhibition of macrophage recruitment, increase in M1 macrophages, and inhibition of the tumor-promoting activity of M2 macrophage [8]. Studies have also illustrated a relationship between the outcome of ovarian cancer and macrophages. e ratio of M1/M2 to M2/TAM was positively associated with overall survival and disease-free survival [9,10]. e density of M2 macrophage in tumor samples was associated with decreased recurrence-free survival [11].
ese studies suggest that macrophages or related markers are potential prognostic factors for ovarian cancer.
In this study, we calculated the level of macrophage infiltration and evaluated the potential of macrophages as prognostic markers. We constructed a gene coexpression network and identified the macrophage-related gene modules. We used the module genes to build a regulatory network related to macrophage infiltration. en, a prognostic model related to macrophage infiltration was constructed using machine learning and verified several datasets. Finally, we explored the relationship between the model and immune factors and carried out enrichment analysis to determine the differences in signal pathways under different groups.

CIBERSORT.
CIBERSORT is an algorithm based on the principle of linear support vector regression, which uses immune cell subtype expression matrix for deconvolution [20]. LM22 in CIBERSORT is a signature gene expression matrix used to estimate the proportion of leukocytes in bulk RNA. Under the operation of the R package "cibersort" and LM22 matrix, we can calculate the composition proportion of 22 kinds of leukocytes of new samples.

Survival Analysis.
Kaplan-Meier analysis is a univariate survival analysis. rough using the R packages "survival" and "survminer," we identified the best cut-off value by the function "res.cut." Receiver operating characteristic curves (ROC) were drawn according to a series of binary classifications. e area under the curve is the AUC value. When the AUC value is greater than 0.5, the result is statistically significant. Univariate Cox regressions were used to identify features related to survival, and multivariate Cox regressions were used to identify multiple features related to survival. Nomograms are multivariate regression analyses that integrate several predictive indicators and use graduated line segments to perform predictive analysis. Calibration curves are scatter plots of actual incidence and predicted incidence to evaluate logistic regression models. ese analyses were performed by the R package "survival."

Weighted Gene Correlation Network Analysis (WGCNA).
WGCNA is a method for analyzing gene expression patterns in multiple samples [21]. e calculation of WGCNA is based on the R package "WGCNA." WGCNA was used to identify highly related gene sets and to construct coexpression networks of gene sets. First, the Pearson correlation between every two genes is calculated. Second, the most suitable β value is calculated to make the network satisfy the scale-free distribution; then the weighted correlation adjacency matrix is constructed by β value. ird, the topological overlap matrix is calculated to add some indirect interactions. Finally, the dynamic cut tree method is used to cluster genes to form gene modules. Genes clustered into the same module indicate that they may have similar functions. Pearson correlation analysis was performed between gene module and clinical trait data. en we can explore the most relevant gene modules for specific clinical traits.

Construction of Protein-Protein Interaction Network.
STRING data is a database for searching protein interaction networks (https://string-db.org/). Cytoscape is the network map visualization software, and we import the STRING results into the Cytoscape to draw protein-protein interaction network [22].

Enrichment Analysis.
Metascape is a powerful gene annotation software package (https://metascape.org/) [23] used to recognize protein or gene function cognition. e database is updated monthly. Gene set enrichment analysis (GSEA) is an enrichment method used to study whether there is a statistically significant difference in the expression level of 2 Journal of Oncology a given group of genes between the two biological states. Gene set variation analysis (GSVA) is a nonparametric analysis method used to calculate pathway enrichment.

Least Absolute Shrinkage and Selection Operator (LASSO).
LASSO is a linear regression method using L1 regularization. rough the operation of a penalty function for compressing partial regression coefficients, a more refined model is obtained.
e process of building a model using LASSO is based on the R language, mainly using R packages "glmnet" and "survival." First, the "glmnet" function was randomly simulated 1000 times to construct the model. en the relationship between penalty coefficient lambda and gene coefficient was established. With the increase of lambda, some gene coefficients become zero, indicating that the gene is an invalid gene of the model. en the 10x cross-validation is repeated 1000 times using the random simulation function "cv.glmnet." When the deviation is minimum, the constructed model was the best, and then the corresponding lambda value was used to calculate the gene coefficient. Finally, we obtained the following prognostic model: risk score � n i (exp i ·coef i ) (where n is the number of genes, exp i is the expression of the ith gene, and coef i is the regression coefficient of the ith gene). We use the R package "survminer" to obtain the best cut-off value of the risk score and separate the cohort into high-risk and low-risk groups [22].
e Human Protein Atlas is a database (https://www.proteinatlas.org/) that provides protein expression in human tissues. Cancer Cell Line Encyclopedia (CCLE, https://portals.broadinstitute. org/ccle/) is used to query gene expression, mutation, copy number, and methylation of thousands of cell lines from tissue origin [25]. We used those tools to query genes expression characteristics.

Statistical Analysis.
e statistical analysis used in this study is based on the environment of R language software (Rx64 3.5.1). Venn diagrams were drawn using online tools (http://bioinformatics.psb.ugent.be/webtools/Venn/). e risk score's distribution map was generated by the R package "pheatmap." e box and correlation scatter charts were drawn by R packages "limma" and "ggplot2."

Calculating the Infiltration Level of 22 Kinds of Leukocytes
in GPL570-OV. GPL570-OV had a total of 351 HGSOC samples, and the number of genes was 20486. We used CIBERSORT to calculate 22 leukocytes infiltration ratios in each sample (Figure 2(a)). We then extracted the infiltration levels of M0, M1, and M2 macrophages and used Kaplan-Meier analysis to determine whether the infiltration of the three kinds of macrophages affected the outcome. We found that outcomes in HGSOC with high infiltration level of M1 macrophages were better (Figure 2(c)), while that of HGSOC with high infiltration level of M2 macrophages was poor ( Figure 2(d)). Although the result of M0 macrophages was not statistically significant, the outcome was better in the case of high infiltration (Figure 2(b)).

Construction of Coexpression Network
Related to Macrophage Infiltration. In the expression network, genes with small gene expression variation represent noise; therefore, we calculated the coefficient of variation of each gene in the GPL570-OV cohort and took the coefficient of variation of the first 8000 genes according to their size to enter the subsequent analysis. First, we clustered the samples, eliminated the samples with high discrete degrees (Figure 3(a)), and then calculated the soft threshold β. We found that when β � 3, R 2 was greater than 0.8, and the mean connectivity was as large as possible (Figure 3(b)). We use soft threshold to build scale-free network and then use dynamic tree cutting method to divide genes into different modules. We obtained 14 gene modules in this manner (Figure 3(c)). We calculated the Pearson correlation between the modules and macrophages ( Figure 3(d)) and found that the black module had the highest correlation with M1 macrophage (R � 0.67). In contrast, the green-yellow module negatively correlated with M0 macrophage (R � -0.26) and positively correlated with M2 macrophage (R � 0.26). Figures 3(e)-3(g) are scatter diagrams of module membership and gene significance.

Construction of Macrophage Infiltration-Associated
Protein Interaction Network. We took the genes in the green-yellow and black modules as macrophage-related genes to enter the following analysis. ere were 89 genes in the green-yellow module and 234 genes in the black module. We input these 323 genes into the STRING database and built a gene interaction network. We found that IFL44L, RSAD2, IFIT3, MX1, IFIH1, IFI44, and ISG15 showed high connectivity and were at the center of the network (Figure 4(a)). After 323 genes were inputted into Metascape for enrichment analysis (Figures 4(b)-4(c)), we found that these genes were significantly enriched in immune-related pathways, such as "defense response to virus" and "regulation of response to biotic stimulus."

Construction of a Macrophage-Related Prognostic Model.
We performed univariate Cox analysis of 323 macrophage infiltration-related genes in GPL570-OV, GPL6480-OV, and TCGA-OV cohorts and selected the statistically significant results (Figures 5(a)-5(b)). Results showed 15 protective genes (HR < 1 and P < 0.05) and no risk genes (HR > 1 and P < 0.05) were obtained. We selected these 15 genes for LASSO analysis in the GPL570-OV cohort as the training set. With the increase of lambda value, the coefficient of some genes decreased to 0, suggesting that the contribution of these genes to the model is small and should be abandoned ( Figure 5(c)). en 10x cross-validation was carried out; when the number of genes was 5, the model reached the optimal solution (Figure 5(d) and Table 1). e formula of the model was as follows: risk score � CD38 e ROC curve showed that the AUC values of 3/5/7 years were all greater than 0.5, and the predictive ability of 7 years was the best ( Figure 5(f )). We drew the expression heat map of five genes ( Figure 5(g)) and the distribution map of risk scores ( Figure 5(h)). Univariate and multivariate analysis showed that the model was predictive and independent of clinical stage and pathological grade ( Figure 5(i)). Finally, we measured the ability of the risk score to predict progressionfree survival (PFS) in GPL570-OV. e model maintains a good ability of hazard stratification ( Figure 5(j), P < 0.05).
e AUC values of three years and five years were greater than 0.5 ( Figure 5(k)).

Analysis of the Expression Characteristics of Five Genes.
We used the TIMER database to analyze the differential expressions of CD38, ACE2, BATF2, HLA-DOB, and WARS between various tumor tissues and normal tissues. e results show that the five genes had different expression characteristics in different tumors. For example, CD38, BATF2, HLA-DOB, and WARS were overexpressed in Head and Neck   Journal of Oncology squamous cell cancer, while they were lower in Kidney Chromophobe (Supplementary Figure 1). We examined the differential expression value of five genes in 426 ovarian cancer tissues and 88 normal tissues in GEPIA. It was found that only the differential expression of WARS was statistically significant, and it was overexpressed in ovarian cancer. Although CD38, BATF2, and HLA-DOB have no statistical significance, they show an increasing trend of expression in ovarian cancer ( Supplementary Figures 2A-2E). e mRNA level of CD38, ACE2, BATF2, and HLA-DOB is relatively low, with an average of 0-3, while the mRNA level of WARS is relatively high, with an average of more than 6. e results were consistent with the immunohistochemistry results (Supplementary Figures 2F-2J), and the protein level of WARS was relatively higher than the other four genes. WARS is mainly expressed in nuclear, while the other four genes are expressed in cytoplasmic/membranous. We also investigated the expression values of five genes in ovarian cancer cell lines (Supplementary Figure 2K). It was found that the expression of mRNA in ovarian cancer cell lines showed the same result,     Journal of Oncology and the mRNA expression level of WARS was significantly higher than those of the other four genes.

Analysis of Model Gene Expression at the Single-Cell Level.
In the CancerSEA database, we queried the related station of ACE2, BATF2, and WARS (Supplementary Figure 3A). We found a positive correlation between ACE2 and angiogenesis/hypoxia/metastasis/quiescence/inflammation pathway (cor > 0.3, P < 0.05) and a negative correlation between ACE2 and DNA-damage/DNA-repair/Invasion (cor < −0.3, P < 0.05). ere was a negative correlation between BATF2 and the invasion/stemness pathway (cor < −0.3, P < 0.05) and there was a positive correlation between BATF2 and the quiescence pathway (cor > 0.3, P < 0.05). We extracted the GSE115007 dataset from the TISCH database and measured expression levels of CD38, BATF2, HLA-DOB, and WARS in immune cells. We found that expression levels of CD38 and BATF2 were low in immune cells ( Supplementary  Figures 3B-3C), while WARS was expressed in M2 macrophages, monocyte, plasma, cDC1, and cDC2 (Supplementary Figure 3D). HLA-DOB was significantly overexpressed in cDC1 and slightly expressed in cDC2 (Supplementary Figure 3E).

Construction of a Nomogram for Clinical Application.
We built a nomogram diagram in the GPL570-OV cohort (Figure 7(a)). In this chart, there were three variables: clinical stage, pathological grade, and risk score used to calculate the 3-/5-/7-year survival rates. ROC and calibration curves were used to appraise the prediction ability of the line chart. e AUC values were all greater than 0.5, indicating that the nomogram had strong prognostic capacity (Figure 7(b)). e calibration curve shows that there was no obvious difference between the predicted value and the measured value (Figure 7(c)).

e Relationship between Macrophage-Related Model and
Immunity. We measured differences of infiltration in 22 leucocytes between high-and low-risk groups. We found that "CD4+ memory resting T cells," "CD4+ memory activated T cells," and "activated mast cells" showed higher infiltration levels in the high-risk group. In contrast, "M1 macrophages" showed higher infiltration levels in the lowrisk group (Figure 8(a)). We also calculated the differences in the expression of 19 HLA (human lymphocyte antigen) molecules in high-and low-risk groups and found that these molecules were overexpressed in the high-risk group (Figure 8(b)). We then calculated the relationship between CD274/PDCD1 and risk scores, and CD274 showed higher expression in the low-risk group (Figure 8(c)    statistical correlation between PDCD1 and risk score (Figure 8(d)). CD14 and CD163 are markers of macrophages. We explored their relationship with risk scores and found that the expression of these two genes was higher in the high-risk group and negatively correlated with the risk score (Figures 8(e)-8(f )).

Enrichment Analysis.
We divided the GPL570-OV cohort into two groups for GSEA, and the gene set was "c2.cp.kegg.v7.0.symbols." e pathways enriched in the high-risk group were "BASAL CELL CARCINOMA" and "RIBOSOME." e low-risk group was enriched with "ANTIGEN PROCESSING AND PRESENTATION," "AUTOIMMUNE THYROID DISEASE," "CYTOSOLIC DNA SENSING PATHWAY," "NATURAL KILLER CELL MEDIATED CYTOTOXICITY," "PRIMARY IMMUNO-DEFICIENCY," "PROTEASOME," "RIG I LIKE RECEP-TOR SIGNALING PATHWAY," "SYSTEMIC LUPUS ERYTHEMATOSUS," "TOLL LIKE RECEPTOR SIG-NALING PATHWAY," and "VIRAL MYOCARDITIS" (Figure 9(a)). Finally, we performed GSVA analysis of the high-risk and low-risk groups and found that "NOTCH SIGNALINGNOTCH SIGNALING" and "WNT BETA CATENIN SIGNALING" were enriched in the high-risk group, while "INTERFERON GAMMA RESPONSE" and "INTERFERON ALPHA RESPONSE" were enriched in the low-risk group (Figure 9(b)).

Discussion
HGSOC has the highest mortality among all gynecological cancers, and most cases show platinum-resistant recurrence [28]. Tumor microenvironment cannot be ignored in the pathogenesis and treatment of ovarian cancer [5]. In this study, we used the GPL570-OV cohort to explore the predictive ability of macrophages and found that both M1 and M2 macrophages predicted the OS of HGSOC. We used the WGCNA to identify the macrophage-related infiltration module and used the genes in the module to construct a prognostic model. e model can carry out risk stratification for 1566 HGSOC samples (n � 1566). We constructed a macrophage-associated infiltration network. We calculated the degree of infiltration of M0/M1/ M2 macrophages using CIBERSORT in the GPL570-OV cohort and measured the prognostic ability of these cells. We found that outcomes in patients with high levels of M1 macrophage infiltration were better, while those of patients with high M2 macrophage infiltration levels were worse. * * * *** *** ** ** *** * * * ***    is is consistent with the fact that M1 macrophages play a role in tumor inhibition, while M2 macrophages help tumor immune escape. en we used WGCNA to construct a gene coexpression network and found that black (M1) and greenyellow (M0/2) associated with macrophage infiltration. We selected the genes to construct the macrophage infiltration network, in which seven genes with high connectivity played a leading role in this network. RSAD2 is an antiviral protein that is significantly upregulated in M1 macrophages [29]. IFIT3 is a marker of M1 macrophage polarization and is highly upregulated in atherosclerosis and other inflammatory diseases [30]. In pneumonia, the microRNA network controls the replication of human macrophages through LGALS8 and MX1 [31]. IFIH1 contributes to the polarization of M1 macrophages in acute respiratory distress syndrome [32]. IFI44 is related to the migration and activation of macrophages [33]. ISG15 secreted by tumor cells increases tumor cell migration and immune escape by inducing M2 macrophage polarization [34]. In the case of ISG15 depletion, M1 macrophages show a robust proinflammatory cytokine expression pattern. ese studies indicated that ovarian cancer cells may inhibit the polarization of M1 macrophages by inhibiting RSAD2, IFIT3, MX1, IFIH1, and IFI44 and induce the polarization of M2 macrophages by secreting ISG15, thereby promoting tumor progression. ese genes are potential therapeutic targets for ovarian cancer.
We constructed a macrophage-related prognostic model of ovarian cancer. e model was made up of CD38, ACE2, BATF2, HLA-DOB, and WARS. Among these genes, WARS is highly expressed in the mRNA and protein levels of patients, as well as the mRNA level of cell lines, and the gene coefficient of WARS is the largest, indicating that WARS plays a core role in the model, while other genes play an auxiliary role. Studies showed that CD38 expression correlated with favorable outcomes by enhancing immune infiltration in the microenvironment of epithelial ovarian cancer [35]. ACE2, also known as ACEH and as angiotensin-converting enzyme 2, is a novel coronavirus cell surface receptor. Studies showed that ACE2 expression positively correlated with immunotherapy response and is a potential protective factor for ovarian cancer [36]. e ACE2/MAS1 axis is involved in the complex regulation of ovarian cancer function [37]. BATF2 has antitumor effects in many tumors: BATF2 combined with p53 to enhance protein stability in gastric cancer, thereby inhibiting ERK phosphorylation [38]. Upregulation of BATF2 inhibited human colorectal cancer cells' growth and epithelialmesenchymal transformation [39]. BATF2 induced an antitumor effect in TAM by upregulating IL-12 expression [40]. e expression of HLA-DOB in multiple myeloma is significantly higher than that in normal plasma cells, suggesting that it is a potential target for immunotherapy [41]. Research shows that WARS compensates for the depletion of IFN-c, thereby inhibiting tumor growth [42]. Single-cell data analysis of ovarian cancer showed that WARS was expressed in various immune cells, which may promote the infiltration of immune cells and inhibiting tumor. WARS may be an important tumor suppressor in OV, and further research on the mechanism is needed in the future.
Our macrophage-related model can carry out risk stratification of GPL570, GPL6480-OV, TCGA-OV, GSE50088, and GSE26712, with a sample size of 1566, indicating that the model is robust and widely applicable. In most previous studies, the scale of ovarian research was less than 1000 [43,44], suggesting that these models are insufficiently accurate. e data from our study were derived from ten independent ovarian cancer cohorts worldwide, suggesting that our model is generalizable. ese cohorts were strictly screened, leaving only HGSOC samples undergoing platinum therapy. Finally, we constructed a nomogram diagram in GPL570-OV which intuitively calculates the OS of HGSOC patients. e AUC values of 3/ 5/7 years were greater than 0.65, indicating that the nomogram diagram has a good predictive ability.
is model can help to explore the immune infiltration mechanism of HGSOC. CD4 memory T cells and activated mast cells showed higher infiltration levels in the high-risk group. In contrast, the degree of M1 macrophage infiltration increased significantly in low-risk group, indicating that CD4+ memory T cells and activated mast cells may be markers of poor outcome and may help tumor immune escape in the immune microenvironment. In contrast, M1 macrophage may be inhibited in high-risk patients. We then found that HLA molecules are overexpressed in high-risk groups, consistent with other studies on HLA molecules in ovarian cancer; HLA-G is a potential biomarker of advanced and complex ovarian cancer [45,46]. HLA-G and HLA-E in ovarian cancer are potentially associated with the mechanism of disease progression [47]. Immune checkpoint inhibitors are new treatments for ovarian cancer. We found that the expression of CD247 increased in the low-risk group. CD247 is mainly expressed by tumor-infiltrating macrophages, not by malignant cells [48], which explains why CD247 is significantly overexpressed in low-risk patients.
Low-risk patients are enriched in immune-signaling pathways, suggesting that when immune-signaling pathways are active, patients are in a low-risk state. e results of GSEA analysis showed that the high-risk group was enriched in "BASAL_CELL_CARCINOMA" and "RIBOSOME," suggesting that ovarian cancer and basal cell carcinoma have similar mechanisms. Recent studies showed that ribosomal ADP-ribosylation inhibits translation and maintains protein homeostasis in ovarian cancer; it is suggested that the ribosome may be an accomplice in the deterioration of ovarian cancer [49]. In the low-risk group of GSEA, many immune-signaling pathways were enriched, including NK cell killing, Toll-like receptor signals, antigen processing, and presentation. It is indicated that the immune system is highly activated in low-risk patients. At present, many studies have shown that the cancer vaccine triggers the immune response of ovarian cancer [5]. According to our research, we can try to use cancer vaccines in high-risk patients to activate the immune pathway. In GSVA analysis, NOTCH and WNT signals were enriched in the high-risk group, while interferon signal was activated in the low-risk group. PARP inhibitors have a specific therapeutic effect on ovarian cancer and show antitumor immunity, which occurs in a manner dependent on interferon gene stimuli and is enhanced by immune checkpoint blockade [50]. erefore, PARP inhibitors combined with immune checkpoint blocking may be an effective treatment for ovarian cancer.
Although our model demonstrated excellent predictive power in 1566 patients with HGSOC, limitations remain. ese ovarian cancer cohorts come from an online database and do not include our cohort, and samples should be expanded for verification before clinical application. e expression of some genes in the cohort is relatively low; therefore, highly sensitive detection techniques are needed in the future. e genes we screened did not further explore their mechanism in ovarian cancer, and further in vivo and in vitro experiments are needed in the future.

Conclusions
We constructed a relationship network related to macrophage infiltration that is helpful to explore the mechanisms of macrophages in HGSOC. We created a prognostic LASSO model based on macrophage-related genes. e model successfully predicted the OS of HGSOC in 1566 samples, which might be useful for assessing the condition and proper treatment of ovarian cancer.

Abbreviations
HGSOC: High-grade serous ovarian cancer TCGA: e Cancer Genome Atlas GEO: Gene Expression Omnibus WGCNA: Weighted Gene Correlation Network Analysis AUC: Area under the curve ROC: Receiver operating characteristic LASSO: Least absolute shrinkage and selection operator coef: Regression coefficient PFS: Progression-free survival OS: Overall survival.
Data Availability e datasets used in the present study are available from the Cancer Genome Atlas database (http://cancergenome.nih. gov/) and GEO database (https://www.ncbi.nlm.nih.gov/ geo/).

Conflicts of Interest
e authors declare that there are no conflicts of interest.