- Research
- Open access
- Published:
APOC1 is a prognostic biomarker associated with M2 macrophages in ovarian cancer
BMC Cancer volume 24, Article number: 364 (2024)
Abstract
Background
Recent studies have demonstrated that APOC1 is associated with cancer progression, exerting cancer-promoting and immune infiltration-promoting effects. Nevertheless, there is currently no report on the presence of APOC1 in ovarian cancer (OV).
Method
In this study, we conducted data analysis using the GEO and TCGA databases. We conducted a thorough bioinformatics analysis to investigate the function of APOC1 in OV, utilizing various platforms including cBioPortal, STRING, GeneMANIA, LinkedOmics, GSCALite, TIMER, and CellMarker. Additionally, we performed immunohistochemical staining on tissue microarrays and conducted in vitro cellular assays to validate our findings.
Result
Our findings reveal that APOC1 expression is significantly upregulated in OV compared to normal tissues. Importantly, patients with high APOC1 levels show a significantly poorer prognosis. Furthermore, our study demonstrated that APOC1 exerted a crucial function in promoting the capacity of ovarian cancer cells to proliferate, migrate, and invade. Additionally, we have identified that genes co-expressed with APOC1 are primarily associated with adaptive immune responses. Notably, the levels of APOC1 in OV exhibit a correlation with the presence of M2 Tumor-associated Macrophages (TAMs).
Conclusion
APOC1 emerges as a promising prognostic biomarker for OV and exhibits a significant association with M2 TAMs in OV.
Introduction
Ovarian cancer (OV) is the second most common malignant tumor and has the highest mortality rate among reproductive system tumors in women [1]. The most prevalent subtype of ovarian cancer is High-Grade Serous Ovarian Cancer (HGSOC), with the majority of cases originating from the fallopian tube epithelium. HGSOC has a low 5-year survival rate of only 30% primarily because most patients are diagnosed at stage III (51%) or stage IV (29%) [2, 3]. As a result, the therapy for HGSOC has changed little over the decades and remains surgical treatment combined with radiation and chemotherapy. Patients with advanced ovarian cancer and mutations in the tumor suppressor genes BRCA1/2 are often treated with maintenance therapy involving Poly (ADP-ribose) polymerase (PARP) inhibitors, such as Olaparib [4, 5]. Nevertheless, not all patients will respond favorably to immunotherapy, and those with OV have limited responses. Therefore, it is crucial to deeply understand the mechanisms of ovarian cancer development and progression to develop improved early diagnostic and therapeutic strategies. Currently, immunotherapy for OV mainly utilizes programmed death-ligand 1 (PD-L1) immune checkpoint inhibitor, but the effect is not obvious. The reason for this lack of clarity is that the targets and mechanisms of effect of OV drugs are poorly understood. Recent research has indicated that the effectiveness of immunotherapy in OV could be linked to the abundance of CD8 + T cells [6, 7]. Even based on the level of T cells, OV has been categorized into three distinct types: T-cell inflamed or “hot” tumors, characterized by infiltration of T cells within deposits (islets) of malignant cells and peripheral stroma; excluded tumors, where T cells are confined to the stroma and absence in the deposits of malignant cells; and non-inflammatory or “immune deserts” or “cold” tumors. These different types exhibit varying degrees of responsiveness to immune checkpoint blockade. Among ovarian cancer subtypes, HGSOC exhibits the highest infiltration of CD8 + T cells, which is associated with a more favorable prognosis [8]. However, the current clinical application of immunotherapy does not take into consideration the immunophenotype. During a phase III clinical trial, the first-line treatment involving the combination of the anti-PDL1 antibody atezolizumab, carboplatin, paclitaxel, and bevacizumab revealed that only approximately 20% of patients who derived benefits exhibited high PD-L1 positivity [9]. Nevertheless, the immune environment is a complex and synergistic system that cannot be accurately explained by a single factor. For instance, ovarian cancer recruits immature myeloid cells and M2 macrophages, which downregulate the expression of CD80, CD86, and IL-12 in effector T cells. Consequently, the identification of additional biomarkers will aid in the selection of appropriate candidates for immunotherapy and facilitate the investigation of mechanisms underlying tumor immunity.
Apolipoprotein C (APOC) is the most abundant apolipoprotein in very low-density lipoprotein (VLDL) cholesterol and is also present in small amounts in low-density lipoprotein (LDL) and high-density lipoprotein (HDL) cholesterol. In humans, there are four different isoforms. Apolipoprotein C1 (APOC1) is a single-chain protein consisting of 83 amino acids initially synthesized in the endoplasmic reticulum. Its molecular weight is 9.3 kDa. The peptide chain undergoes shearing to produce mature APOC1, which contains 57 amino acids and has a molecular weight of 6.6 kDa, making it the smallest apolipoprotein [10]. APOC1 is engaged in normal biological processes and disease development. It is found in chylomicrons (CM), VLDL, and HDL, and serves as an exchange apolipoprotein between various lipoproteins, playing a vital role in lipid homeostasis [11]. Additionally, APOC1 plays an essential role in diabetes, atherosclerosis, and Alzheimer’s disease [12, 13]. Recently, studies on APOC1 have increasingly emphasized its involvement in cancer. In the study of breast cancer, APOC1 has been identified as a potential diagnostic and classification biomarker [14]. Furthermore, APOC1 has been found to promote the development of gastric cancer [15]. In a study on glioblastoma, APOC1 was discovered to enhance tumorigenesis by reducing KEAP1/NRF2 and CBS-regulated ferroptosis [16]. Recently, it has been revealed that APOC1 can be used as an immunological biomarker, influencing macrophage polarization and facilitating renal cell cancer [17]. In hepatocellular carcinoma, a single-cell RNA sequencing study demonstrated that inhibiting APOC1 promoted the transformation of M2 macrophages into M1 macrophages via the ferroptosis pathway, thereby enhancing anti-PD1 immunotherapy [18]. However, the function of APOC1 in OV hasn’t been reported.
The aim of this research was to survey the function of APOC1 in OV as well as to evaluate its potential as an immune biomarker.
Materials and methods
Data resources
The fragments per kilobases per million mapped fragments (FPKM) files containing clinical information were derived from The Cancer Genome Atlas (TCGA). Expression matrix data in SOFT format, including GSE12470 [19], GSE66957 (Cheng et al., unpublished data, 2015), GSE10971 [20] and GSE52037 [21], were obtained from the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) (GEO). Somatic mutation and copy number variation (CNV) analyses of APOC1 in OV were carried out using the cBioPortal website (www.cbioportal.org) [22] (Date of data acquisition and analysis: 2023-06-11). Additionally, Kaplan-Meier overall survival analysis of APOC1 expression in ovarian cancer patients was performed using the Kaplan-Meier Plotter tool (https://kmplot.com/analysis/) [23, 24] (Date of data acquisition and analysis: 2023-06-11).
Analysis of differentially expressed genes
The differential analysis of genes was conducted using the limma R package (version 3.52.4) [25]. DEGs were identified for genes with p < 0.05 &|log2FC|≥1.
Protein-protein interaction analysis
The protein-protein interaction (PPI) network analysis of APOC1 was performed using the online tool STRING 11.0 b (https://string-db.org/) [26] (Date of data and analysis and analysis: 2023-06-11). In order to further investigate these relationships between interacting genes further, we utilized the online tool GeneMANIA (http://www.genemania.org) [27] (Date of data acquisition and analysis: 2023-06-13).
LinkedOmics and GSCALite
The genes significantly associated with APOC1 were identified through screening TCGA data, and these genes were enriched using the online tool LinkedOmics (www.linkedomics.org) [28] (Date of data acquisition and analysis: 2023-06-13). Subsequently, pathway activity analysis was conducted using the online tool GSCALite (www.bioinfo.life.hust.edu.cn/web/GSCALite/) [29] (Date of data acquisition and analysis: 2023-06-13).
Immune infiltration analysis
The immune correlation assessment was performed using the ssGSEA algorithm provided in the GSVA R package (version 1.44.5) [30]. In order to calculate the immune infiltration corresponding to the TCGA data, markers for 24 immune cells were utilized [31]. The analysis of the correlation between APOC1 expression and immunity cells was performed utilizing the TIMER 2.0 platform (http://timer.cistrome.org/) [32] (Date of data acquisition and analysis: 2023-08-18). Specifically, CD206, CD163, IL10, and ARG1 were selected as markers for macrophage M2.
Single-cell RNAseq data analysis
The single-cell sequencing data (GSE147082 [33] and GSE165404 [34]) were obtained from the GEO. Cell marker analysis of single cell sequencing data utilizing the online tool CellMarker 2.0 (http://biocc.hrbmu.edu.cn/CellMarker/) [35] (Date of data acquisition and analysis: 2023-08-25).
Clinical sample
Between 2017 and 2022, a total of 61 HGSOC tissues and 13 fallopian tubes (FT) tissues were obtained from Shanghai General Hospital. Prior to the operations, patients were informed about the procedures and provided their consent. Ethical approval for the study was obtained from the Medical Ethics Committee of Shanghai First People’s Hospital. Informed consent was taken from participants to participate in the study.
Production of issue microarray
The tissue samples were fixed in tissue fixative for more than 24 h. Subsequently, the tissues were trimmed using a scalpel in a fume hood and placed in embedding cassettes with appropriate labels. To remove water from the tissues, they were sequentially immersed in different concentrations of ethanol and xylene for dehydration. The dehydrated tissues were then immersed in preheated liquefied wax and left to soak overnight. Using the HistoCore Arcadia system (HistoCore Arcadia, Leica, Germany), tissue wax blocks were created. Finally, the prepared tissue wax blocks were sent to Servicebio (Wuhan, China) for tissue microarray production. The tissue samples were initially stained with HE and then localized by two pathologists before being sampled for tissue microarrays. Subsequently, the scoring of the IHC staining was performed by ImageJ software, as described in the Immunohistochemistry section in methods.
Immunohistochemistry (IHC)
To melt the surface sealing wax, the tissue microarray was placed in a Thermostatic Incubator (DNP-9052, JINGHONG, Shanghai) at 60 °C overnight. Deparaffinization and rehydration of the tissue microarrays were carried out by immersing them in xylene and ethanol concentration gradients in a fume hood. For antigen retrieval Citrate Antigen Extraction Solution (P0081; Beyotime, China) was used in boiling water for 7 min. To eliminate endogenous peroxidase activity, a 3% hydrogen peroxide solution was applied, and the tissues were then sealed with 10% goat serum (C0265; Beyotime, China) for 1 h. After rinsing with PBS, the tissue microarrays were incubated with the anti-APOC1 antibody (EPR16813; Abcam, UK) diluted to 1:800, the anti-CD163 antibody (16646-1-AP, Proteintech, China) diluted to 1:1000 and the anti-CD206 antibody (60143-1-Ig, Proteintech, China) diluted to 1:10000 at 4 °C overnight. The Immunohistochemistry Kit (GK500705; Genentech, China) was used for secondary antibody incubation and staining of the tissue microarrays. Hematoxylin (C0107; Beyotime, China) was applied for staining the tissue microarrays for 3 min, followed by termination in distilled water. Sequential dehydration of the tissue microarrays was achieved by immersing them in a gradient solution of xylene and ethanol. Finally, the slices were closed with neutral resin (GT100519; GeneTech, China). The figures were captured utilizing a microscope (Leica, London, UK). The obtained results were analyzed using ImageJ software (version 1.52a; National Institutes of Health) [36] and the IHC Profiler plugin. Cytoplasmic staining scores were calculated using Plugins -> IHC Profiler->Cytoplasmic Stained Image->H DAB, while nuclear staining scores were obtained by Plugins -> IHC Profiler->Nuclear Stained Image->H DAB->Set Threshold, followed by Plugins -> Macros -> IHC Profiler. Finally, Immunohistochemical score = cytoplasmic staining score * nuclear staining score. The four IHC Profiler scores were defined as follows: 4: High positive (gray scale: 0–60); 3: Positive (gray scale value: 61–120); 2: Low Positive (gray scale value: 121–180); 1: Negative (gray scale value: 181–236).
Given that our tissue microarrays primarily target regions with a high concentration of tumor cells, which could potentially influence the IHC results, we conducted an analysis of the correlation between APOC1 expression and the expression of CD163 and CD206 using the average optical density (AOD) values obtained through Image-Pro Plus (version 6.0; Media Cybernetics, Rockville, MD). The procedure is as follows:
-
1.
Opening the picture, click measure -> intensity, click new -> std. optical density -> options -> image in the intensity box, and then selecting the blank place in the figure->ok.
-
2.
Changing the incidental level to the value of blank place and then clicking measure->count/size->selecte colours->click the icon of the pen.
-
3.
Selecting the immunohistochemistry picture of the positive protein expression region, and clicking close after the selection.
-
4.
Finally getting the mean of area and mean of Integrated Optical Density (IOD) by clicking Measure->Select Measurements, selecting iod->ok->count inside and clicking view->statistic.
AOD = mean of IOD/ mean of area.
Cell growth and cell culture
Hey, Caov3, 293T, and tubal epithelial cells OE E6/E7 were acquired from the National Collection of Authenticated Cell Cultures (Shanghai, China). THP-1 cells were provided by Servicebio (Wuhan, China). Hey, Caov3, 293T, and OE E6/E7 were cultured in high-sugar DMEM medium (319-005-CL; MUTICELL, China) with 10% fetal bovine serum (086–150; MUTICELL, China). THP-1 cells were cultured in 1640 DMEM medium (350-000-CL; MUTICELL, China) with 10% fetal bovine serum (086–150; MUTICELL, China). The cells were cultured in a growth environment of 37 °C and 5% CO2 in a CO2-Incubator (51,023,126, Thermo Fisher Scientific, USA).
Plasmid construction and transfection
The plasmids designed for APOC1 knockdown were synthesized by Genomeditech (Shanghai, China) and were validated through DNA sequencing. To generate the shAPOC1 plasmid, the shRNA primer pair specific for APOC1 was derived from the shRNAlibrary (TRC) and inserted into the pLKO.1-Puro vector. Cells transfected with the APOC1 knockdown plasmids were labeled as shAPOC1-1 and shAPOC1-2, while cells transfected with the control plasmid were labeled as shNC. 293T cells were cultured in 6 cm dishes until reaching a cell density of 50%. Transfection was carried out using Lipofectamine 3000 reagent (Invitrogen Life Technologies, USA) based on the protocol provided by the manufacturer. The 293T cells were then incubated at 37 °C with 5% CO2 for 48 h. The supervisory fluid was gathered and centrifuged to obtain the lentiviral solution. This lentiviral solution was added to Hey and Caov3 cells in six-well plates at a 40% cell density with varying concentration gradients. The cells were subsequently incubated at 37 °C with 5% CO2 for 48 h, followed by replacement with 4 ml of DMEM containing puromycin (2 mg/ml) for screening. After 24–48 h, the cells with the highest viability in the six-well plate were selected for expansion. Subsequent validation was conducted using the Western Blot assay. The sequence of the short hairpin ribonucleic acid (shRNA) used is as follows:
shAPOC1-1: 5′- GACATTTCAGAAAGTGAAGGA − 3′.
shAPOC1-2: 5′- GCTGAAGGAGTTTGGAAACAC − 3′.
Western blotting
Protein extraction was performed using RIPA lysis buffer (P0013C; Beyotime, China) supplemented with 1% phenylmethanesulfonyl fluoride (ST507-10 ml; Beyotime, China). The protein lysates were separated with a 15% PAGE Gel Fast Preparation Kit (Epizyme, China) and subsequently transferred to methanol-activated polyvinylidene fluoride membranes (PVDF; Millipore, USA). After membrane transfer, the membranes were placed in 5% skimmed milk (E504BA0014, BBI Life Sciences, China) at room temperature for 1 h. Subsequently, the membranes were rinsed three times with 0.1% Tris-HCl plus Tween-20 (TBST) for 5 min each. After specific detection of the target antigen confirmation, we cut between 25 kDa − 35 kDa and incubated the upper half with anti-GAPDH polyclonal antibody and the lower half with anti-APOC1 polyclonal antibody. In the article show images we marked with red boxes. Images of specific detection of the target antigen and original blots can be seen in the Supplementary file. The membranes were then placed with primary and corresponding secondary antibodies, followed by three additional rinses. Detection was carried out using the ECL luminescence kit (BioVision, USA) on a chemiluminescence imaging system (Tanon 5200; Tanon, China). The antibodies and their respective dilutions used in the experiments were as follows:
Anti-APOC1 antibody (ab205718; Abcam, UK) diluted to 1:1000.
Anti-GAPDH antibody (60004-1-lg; Proteintech, China) diluted to 1:50000.
Anti-rabbit (SA00001-2; Proteintech, China) diluted to 1:2000.
Anti-mouse (SA00001-1; Proteintech, China) diluted to 1:2000.
Quantitative comparison of the relative protein levels was performed using ImageJ software (Version 1.52a; National Institutes of Health).
RNA extraction and qPCR
Total RNA extraction from THP-1 cells was performed using TRIeasy™ LS Total RNA Extraction Reagent (19201ES60; Yeasen, China) following the manufacturer’s protocol. Subsequently, cDNA synthesis was carried out using a reverse transcription reagent (R202-02; EnzyArtisan, China). The mRNA expression levels were normalized to GAPDH and calculated using the 2-ΔΔCt method [37]. The primer sequences for qPCR and the cycling conditions can be found in Table S1.
Cell proliferation assay
The sh-APOC1 cells were seeded into 96-well plates with a density of 2000 cells per well and cultured for 0, 24, 48, and 72 h. Cell proliferation capacity was evaluated with the Cell Counting Kit-8 (C0038; Beyotime, China) according to the manufacturer’s instructions. After 1 h of incubation, the optical density (OD) values were measured at 450 nm using a multifunctional enzyme marker (VLBLATGD2; Thermo Fisher Scientific, USA). A blank control consisting of DMEM with 10% FBS was used.
Colony-formation assay
The sh-APOC1 cells were cultured in six-well plates with a seeding density of 800 cells/well for Hey cells and 1500 cells/well for Caov3 cells. The plates were placed in a 37 °C incubator with 5% CO2 for 14 days. After the cultured period, the cells were fixed with 4% paraformaldehyde (BL539A; Biosharp, China) for 30 min. Subsequently, the cells were stained with crystal violet staining solution (C0121; Beyotime, China) for 20 min.
Transwell cell migration and invasion assays
For the Transwell cell migration assay, 200 µl of serum-free medium containing 1 × 105 Hey cells or 200 µl of serum-free medium containing 1 × 105 Caov3 cells were added to the upper chamber (14,421,030; Corning Incorporated, USA). 700 µl of DMEM medium with 10% FBS was placed in the lower chamber and cultured for 24 h. (Caov3 cells: 48 h). After the incubation period, the cells were fixed with 4% paraformaldehyde for 30 min and subsequently stained with crystal violet staining solution for 20 min.
For the Transwell cell invasion assay, 60 µl of Matrigel substrate (356,234; Corning Incorporated, USA), diluted 1:6 in serum-free medium, was added to the upper chamber of the transwell and allowed to solidify for 1 h at 37 °C before use. Cells and medium were added to the transwell chamber as in the cell migration assay. Finally, after culturing for 48 h (Caov3 cells: 72 h), the cells were fixed and stained.
Generation and differentiation of macrophages
THP-1 cells were pre-treated with 100 ng/ml PMA (S1819-1 mg; Beyotime, China) for 24 h to induce the generation of M0 macrophages. These M0 macrophages were subsequently stimulated with 20 ng/ml IL-13 (P5178-10 µg; Beyotime, China) and 20 ng/ml IL-4 (P5129-5 µg; Beyotime, China) for 48 h to differentiate into M2 macrophages. The M0 macrophages were then co-cultured with HGSOC cells (Hey or Caov3) in a 6-well transwell chamber with a pore size of 0.4 μm (3470; Corning Incorporated, USA). After co-culturing for 48 h, the macrophages were collected to obtain Tumor-associated Macrophages (TAMs) for qPCR analysis.
Flow cytometry
Induced macrophages from 6-well plates were collected, centrifuged, washed once with PBS, and then centrifuged again to obtain a single-cell suspension. The cells were subsequently treated with Permeabilization Buffer (00-8333-56, Thermo Fisher Scientific, USA) following the manufacturer’s protocol. Following this, they were incubated with FITC Mouse Anti-human CD163 (333,617, BioLegend, USA), APC Mouse anti-Human CD206 (321,109, BioLegend, USA), and PE Mouse anti-Human ARG1 (369,703, BioLegend, USA) at 4 °C for 30 min. Finally, the stained cells were washed once, resuspended in 200 µl PBS solution, and analyzed on a Fortessa flow cytometer (BD, Biosciences, USA) using FlowJo software (version 10.0, FlowJo, USA).
Statistical analysis
All experiments in this study were repeated three times. All of our data were initially assessed for normal distribution using the Shapiro-Wilk test. For normally distributed data, unpaired and two-tailed student’s t-tests were conducted. For non-normally distributed data, the Mann-Whitney test was used. All statistical analyses were performed using SPSS software (version 19.0; Chicago, IL, USA). R (version 4.0.3; https://www.r-project.org/) was utilized to analyze data from TCGA and GEO. All graphs were plotted using GraphPad Prism software (version 8.0; San Diego, CA). All experiments were repeated three times. p < 0.05 was deemed statistically significant (*p < 0.05, **p < 0.01, ***p < 0.001).
Results
The expression of APOC1 was upregulated in human OV compared to FT
First, the mutation and amplification of APOC1 rank as the sixth highest in OV compared to pan-cancer (Fig. 1A). By analyzing the GEO datasets GSE12470, GSE66957, GSE10971, and GSE52037, we observed elevated expression of APOC1 in OV (Fig. 1B). Furthermore, through Western Blot assay, APOC1 expression was found to be elevated in the OV cell lines Hey and Caov3 compared to the fallopian tube epithelial cell line OE E6/E7 (Fig. 1C). Consequently, the Hey and Caov3 cell lines were selected for vitro assays. Additionally, immunohistochemical staining of the HGSOC tissue microarray demonstrated increased APOC1 expression in HGSOC compared to FT (Fig. 1D-E, S1A). According to the immunohistochemical scores acquired by IHC, we classified these patients in two groups: the APOC1 hyper-expression group (n = 35) and the lower-expression group n = 26). Furthermore, there was no significant correlation between APOC1 expression and age, FIGO stage, CA125, HE4, lymph node metastasis, large omental metastasis, ascites, BRCA1/2 mutation and residual tumor size (P > 0.05) (Table 1). The Kaplan-Meier overall survival analysis revealed that ovarian cancer patients with high expression of APOC1 exhibited poorer overall survival compared to patients with low expression of APOC1 (Fig. 1F). In conclusion, these findings indicate that APOC1 is high-expressed in OV and is related to poor patient prognosis.
APOC1 enhances the proliferation, migration and invasion of OV cells in vitro
First, a Western Blot experiment was carried out to validate the effective knockdown of APOC1 in Hey and Caov3 cells (Fig. 2A). Subsequently, a CCK-8 assay was conducted to evaluate cell proliferation, demonstrating that decreased APOC1 expression hindered cell proliferation (Fig. 2B). Furthermore, a colony formation assay indicated that reduced APOC1 expression impaired the cells’ ability to form colonies (Fig. 2C). Additionally, the cell migration assay revealed a significant inhibition of cell migration with reduced APOC1 expression. Moreover, the cell invasion assay demonstrated that decreased APOC1 expression suppressed cell invasion (Fig. 2D-E). In conclusion, these findings indicate that APOC1 plays a critical function in regulating the proliferative, migratory, and invasive capacities of OV cells in vitro.
Potential mechanism of APOC1 in OV
We utilized STRING and GeneMANIA to conduct PPI network analysis of APOC1, aiming to explore potential interactions between proteins associated with APOC1. The PPI network yielded 11 nodes and 55 edges (Figure S2A). The identified APOC1-related genes were functionally linked to several processes, including cholesterol metabolism, vitamin digestion and absorption, fat digestion and absorption, and lipid particle composition. Additionally, GeneMANIA analysis revealed that the functions of these APOC1-related genes were primarily related to the remodeling of protein-lipid complexes, plasma lipoprotein particle remodeling, remodeling of protein-containing complexes, organization of plasma lipoprotein particles, organization of protein-lipid complex subunits, and lipoprotein particles (Figure S2B). To further investigate the potential mechanism of APOC1 in OV, we employed LinkedOmics and GSCALite online platforms to discover co-expressed genes and associated pathways in the TCGA data. The LinkedOmics platform results revealed that APOC1 positively correlated with 8,299 genes (depicted as red dots) and negatively correlated with 11,732 genes (depicted as green dots) in OV (Figure S2C). We have provided the top 10 genes that positively and negatively associate with APOC1 in OV in Tables 2 and 3 respectively. Notably, TYROBP (cor = 0.849397, p = 1.53e-85), FCER1G (cor = 0.838395, p = 2.56e-81), CD48 (cor = 0.812535, p = 1.58e-72), LST1 (cor = 0.812054, p = 1.00e-85), and APOE (cor = 0.808613, p = 2.58e-71) displayed the strongest correlation with APOC1 in OV (Figure S2D). We explored them utilizing the pathway activity module of the GSCALite platform to determine whether these six genes (TYROBP, FCER1G, CD48, LST1, APOE, and APOC1) act through specific cancer pathways. The results indicated that APOC1 could promote processes such as apoptosis, cell cycle, estrogen receptor (ER), and signaling epithelial-mesenchymal transition (EMT), while it inhibited androgen receptor (AR), PI3K/AKT, receptor tyrosine kinase (RTK), and TSC/mTOR signaling pathways (Fig. 3A). Moreover, we conducted enrichment analysis on the top 500 genes most relevant to APOC1, according to KEGG pathway analysis. The upregulated pathways included, natural killer cell-mediated cytotoxicity, osteoclast differentiation, cell adhesion molecules (CAMs), phagosome, cytokine-cytokine receptor interaction, chemokine signaling pathway, and complement and coagulation cascades. Conversely, the downregulated pathways included the AMP-activated protein kinase (AMPK) signaling pathway, signaling pathways regulating pluripotency of stem cells, and RNA transport (Fig. 3B). Additionally, GO terms revealed that the co-expressed genes of APOC1 primarily activated activities related to adaptive immune response, neutrophil-mediated immunity, response to molecule of bacterial origin, leukocyte differentiation, and positive regulation of defense response. Activities such as appendage development, protein alkylation, regulation of mRNA metabolic process, cilium organization, and microtubule-based movement were inhibited (Fig. 3C).
Association of APOC1 expression with immune cells in OV
The analysis conducted on our LinkedOmics platform suggests that the function of APOC1 in OV primarily revolves around the adaptive immune response. Therefore, it is possible that APOC1 is linked to the immune response in OV. Through immune infiltration analysis of TCGA data (Fig. 4A), we observed that T cells, cytotoxic cells, and macrophages exhibited the highest correlation with APOC1 expression. To explore further the association between the expression of APOC1 and immune cells, we utilized the TIMER 2.0 platform. Using the EPIC algorithm, we discovered a strong positive correlation between APOC1 expression and macrophages (Rho = 0.705, p = 1.03e-38) (Fig. 4B). Moreover, employing the XCELL algorithm, we determined that APOC1 expression positively correlated with B cells (Rho = 0.394, p = 1.10e-10), myeloid dendritic cells (Rho = 0.452, p = 5.74e-14), macrophages (Rho = 0.755, p = 2.97e-47), M1 macrophages (Rho = 0.737, p = 7.33e-44), and M2 macrophages (Rho = 0.613, p = 4.35e-27) (Fig. 4C). These findings suggest a specific association between APOC1 expression and macrophage cells.
APOC1 is associated with M2 phenotypic TAMs in OV
Based on these findings, we conducted further investigations into the relationship between APOC1 and macrophages in OV. Through the analysis of single-cell sequencing data from GSE147082 and GSE150864 (Fig. 5A), we discovered that APOC1 was predominantly enriched in macrophages in OV. To gain a deeper understanding, we established an in vitro model of TAMs. Initially, THP-1 cells were treated with PMA to convert them into M0 macrophages. Subsequently, a portion of these cells were further treated with IL-13 and IL-4 for 48 h to induce the conversion into M2 macrophages. Another portion of the M0 macrophages were co-cultured with ovarian cancer cells (Hey and Caov3 cells) for 48 h, leading to their transformation into TAMs. Using a qPCR assay, we assessed the expression levels of common markers of M2 macrophages, namely CD163, CD206, IL10, and ARG1. The results indicated a high expression of CD163, CD206, IL10, and ARG1 in TAMs (Fig. 5B). Additionally, we employed an APOC1 knockdown cell line derived from OV to transform TAMs. Notably, the qPCR results showed a decrease in the expression of CD163, CD206, IL10, and ARG1 in these TAMs (Fig. 5C). Collectively, these findings suggest an association between APOC1 and M2 TAMs in OV.
Subsequently, flow cytometry assays were performed, and the results validated the conclusions of our previous assays (Fig. 6A). Following this, IHC staining of tissue microarrays for CD163 and CD206 was conducted, and the AOD values of APOC1, CD163, and CD206 were calculated (Fig. 6B-C, S1B-C). A correlation analysis was then conducted, revealing a significant correlation between the expression of APOC1 and CD163 (p = 0.0083), while no significant correlation was observed with the expression of CD206 (p = 0.4184) (Table 4).
Discussion
Tumor promotion inflammation represents a significant characteristic of tumors, whereby the inflammatory response provides the tumor microenvironment with various bioactivating molecules. These molecules encompass growth factors that sustain cancer cell proliferation and angiogenic factors that facilitate vascular growth, tumor metastasis, and invasion [38]. Among these factors, TAMs are a type of macrophages that infiltrate tumor tissues or accumulate in the microenvironment of solid tumors. TAMs possess the ability to influence tumor growth, angiogenesis, immunomodulation, and drug resistance [39]. While several current articles have established the association between APOC1 and tumor immunomodulation, as well as the phenotypic transition of macrophages, there is a lack of exploration regarding the function of APOC1 in OV. In this study, combining biological information analysis with in vitro cell assays, we aim to investigate the potential association between APOC1 and tumor immunity, as well as cancer promotion.
To begin with, we conducted an analysis of somatic mutations and CNV of APOC1 across multiple cancer types. Among these, ovarian cancer ranked sixth in terms of APOC1 alterations. Furthermore, by analyzing data from GEO, including GSE12470, GSE52037, GSE66957, and GSE10971, we observed higher expression levels of APOC1 in OV compared to normal tissues. To validate this finding, we performed immunohistochemical staining on a tissue microarray, which confirmed the elevated expression of APOC1. Previous studies have suggested that APOC1 may serve as a potential biomarker for various cancers [17, 40,41,42], and our analysis indicates its potential as a biomarker for OV as well. Additionally, there have been reports indicating that APOC1 promotes the progression of gastric cancer and the proliferation of prostate cancer [15, 43]. In our in vitro cellular assays, we discovered that APOC1 promotes the proliferation, migration, and invasive abilities of ovarian cancer cells.
Following that, we performed an analysis of functional enrichment on the top 500 APOC1-related genes. This analysis revealed that these APOC1-related genes are primarily involved in adaptive immune response, neutrophil mediated immunity, leukocyte differentiation, and the response to molecules of bacterial origin. These findings led us to suspect that APOC1 may be associated with the immune response in OV. In our analysis of the association between APOC1 and immunocytes, we discovered a strong association between APOC1 and macrophages in OV. Additionally, using single-cell sequencing data (GSE147082 and GSE150864), we found that APOC1 was predominantly enriched in macrophages within the OV microenvironment. Subsequently, we constructed an in vitro model of TAMs in OV and analyzed the expression of common markers in M2 macrophages through qPCR and flow cytometry. Remarkably, we observed a correlation between APOC1 and markers such as CD163, CD206, IL10, and ARG1 in M2 macrophages. Based on these findings, we hypothesized that APOC1 may be associated with M2 macrophages in OV. We subsequently conducted IHC staining of tissue microarrays for CD163 and CD206, followed by a correlation analysis with APOC1. The analysis revealed a correlation between the expression of APOC1 and CD163, but not with the expression of CD206.
In studies on inflammation, M1 macrophages are typically believed to promote inflammation, whereas M2 macrophages are generally believed to inhibit inflammation. However, this is different in TAMs. Due to the crucial role of immune escape in tumors, M2 macrophages could suppress inflammation and enhance tumor proliferation. Consequently, M2 macrophages have garnered significant attention in tumor research. Activation of M2 macrophages can give rise to four distinct types: M2a, M2b, M2c, and M2d. M2a macrophages enhance endocytosis, M2b macrophages produce pro-inflammatory and anti-inflammatory cytokines, M2c macrophages phagocytose apoptotic cells, and M2d macrophages promote angiogenesis and tumor progression [44, 45]. This highlights the significance of M2 macrophage polarization in tumors. It remains to be investigated in future studies whether there is an association between APOC1 and the polarization of M2 macrophages in OV. However, our study demonstrated an association between APOC1 and M2 macrophages in OV. APCO1 holds potential as a prognostic biomarker.
Data availability
The FPKM files of RNA-seq transcriptome data used to support the findings of this study were obtained from the TCGA database. The Matrix files in SOFT format for four GEO datasets (GSE12470, GSE66957, GSE10971 and GSE52037) used to support the findings of this study were obtained from in the Gene Expression Omnibus. The two single-cell sequencing datasets (GSE147082 and GSE165404) used to support the results of this study were obtained from the Gene Expression Omnibus. The data of somatic mutation and CNV of APOC1 in OV were carried out from the TCGA database.
Abbreviations
- HGSOC :
-
High grade serous ovarian cancer
- OV:
-
Ovarian Cancer
- TAMs:
-
Tumor-associated Macrophages
- APOC:
-
Apolipoprotein C
- VLDL:
-
low-density lipoprotein
- HDL:
-
high-density
- LDL:
-
low-density lipoprotein
- APOC1:
-
Apolipoprotein C1
- PPI:
-
Protein and protein interaction
- shRNA:
-
short hairpin ribonucleic acid
- FPKM:
-
Fragments per kilobase of exon model per million mapped fragments
- TCGA:
-
The Cancer Genome Atlas
- GEO:
-
Gene Expression Omnibus
- CNV:
-
Copy number variation
- GO:
-
Gene Ontology
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- GSEA:
-
Gene set enrichment analysis
- RNA-seq :
-
RNA sequencing
- FT:
-
Fallopian tubes
- FTE:
-
Fallopian tube epithelium
- IOD:
-
Integrated Optical Density
- AOD:
-
Average Optical Density
References
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. CA Cancer J Clin. 2022;72(1):7–33.
Matulonis UA, Sood AK, Fallowfield L, Howitt BE, Sehouli J, Karlan BY. Ovarian cancer. Nat Rev Dis Primers. 2016;2:16061.
Armstrong DK, Alvarez RD, Backes FJ, Bakkum-Gamez JN, Barroilhet L, Behbakht K, Berchuck A, Chen LM, Chitiyo VC, Cristea M, et al. NCCN Guidelines® insights: ovarian Cancer, Version 3.2022. J Natl Compr Canc Netw. 2022;20(9):972–80.
Tattersall A, Ryan N, Wiggans AJ, Rogozińska E, Morrison J. Poly(ADP-ribose) polymerase (PARP) inhibitors for the treatment of ovarian cancer. Cochrane Database Syst Rev. 2022;2(2):Cd007929.
DiSilvestro P, Banerjee S, Colombo N, Scambia G, Kim BG, Oaknin A, Friedlander M, Lisyanskaya A, Floquet A, Leary A, et al. Overall survival with maintenance olaparib at a 7-Year Follow-Up in patients with newly diagnosed Advanced Ovarian Cancer and a BRCA mutation: the SOLO1/GOG 3004 Trial. J Clin Oncol. 2023;41(3):609–17.
Kandalaft LE, Dangaj Laniti D, Coukos G. Immunobiology of high-grade serous ovarian cancer: lessons for clinical translation. Nat Rev Cancer. 2022;22(11):640–56.
Wan C, Keany MP, Dong H, Al-Alem LF, Pandya UM, Lazo S, Boehnke K, Lynch KN, Xu R, Zarrella DT, et al. Enhanced efficacy of simultaneous PD-1 and PD-L1 Immune Checkpoint Blockade in High-Grade Serous Ovarian Cancer. Cancer Res. 2021;81(1):158–73.
Chen DS, Mellman I. Elements of cancer immunity and the cancer-immune set point. Nature. 2017;541(7637):321–30.
Moore KN, Bookman M, Sehouli J, Miller A, Anderson C, Scambia G, Myers T, Taskiran C, Robison K, Mäenpää J, et al. Atezolizumab, Bevacizumab, and Chemotherapy for newly diagnosed stage III or IV ovarian Cancer: placebo-controlled Randomized Phase III Trial (IMagyn050/GOG 3015/ENGOT-OV39). J Clin Oncol. 2021;39(17):1842–55.
Trenchevska O, Schaab MR, Nelson RW, Nedelkov D. Development of multiplex mass spectrometric immunoassay for detection and quantification of apolipoproteins C-I, C-II, C-III and their proteoforms. Methods. 2015;81:86–92.
Swertfeger DK, Li H, Rebholz S, Zhu X, Shah AS, Davidson WS, Lu LJ. Mapping atheroprotective functions and related Proteins/Lipoproteins in size fractionated human plasma. Mol Cell Proteom. 2017;16(4):680–93.
Rouland A, Masson D, Lagrost L, Vergès B, Gautier T, Bouillet B. Role of apolipoprotein C1 in lipoprotein metabolism, atherosclerosis and diabetes: a systematic review. Cardiovasc Diabetol. 2022;21(1):272.
Zhou Q, Peng D, Yuan X, Lv Z, Pang S, Jiang W, Yang C, Shi X, Pang G, Yang Y, et al. APOE and APOC1 gene polymorphisms are associated with cognitive impairment progression in Chinese patients with late-onset Alzheimer’s disease. Neural Regen Res. 2014;9(6):653–60.
Beňačka R, Szabóová D, Guľašová Z, Hertelyová Z, Radoňák J. Classic and New Markers in Diagnostics and Classification of Breast Cancer. Cancers (Basel) 2022, 14(21).
Gu Q, Zhan T, Guan X, Lai C, Lu NA, Wang G, Xu L, Gao X, Zhang J. Apolipoprotein C1 promotes tumor progression in gastric cancer. Oncol Res. 2023;31(3):287–97.
Zheng XJ, Chen WL, Yi J, Li W, Liu JY, Fu WQ, Ren LW, Li S, Ge BB, Yang YH, et al. Apolipoprotein C1 promotes glioblastoma tumorigenesis by reducing KEAP1/NRF2 and CBS-regulated ferroptosis. Acta Pharmacol Sin. 2022;43(11):2977–92.
Ren L, Yi J, Yang Y, Li W, Zheng X, Liu J, Li S, Yang H, Zhang Y, Ge B, et al. Systematic pan-cancer analysis identifies APOC1 as an immunological biomarker which regulates macrophage polarization and promotes tumor metastasis. Pharmacol Res. 2022;183:106376.
Hao X, Zheng Z, Liu H, Zhang Y, Kang J, Kong X, Rong D, Sun G, Sun G, Liu L, et al. Inhibition of APOC1 promotes the transformation of M2 into M1 macrophages via the ferroptosis pathway and enhances anti-PD1 immunotherapy in hepatocellular carcinoma based on single-cell RNA sequencing. Redox Biol. 2022;56:102463.
Yoshihara K, Tajima A, Komata D, Yamamoto T, Kodama S, Fujiwara H, Suzuki M, Onishi Y, Hatae M, Sueyoshi K, et al. Gene expression profiling of advanced-stage serous ovarian cancers distinguishes novel subclasses and implicates ZEB2 in tumor progression and prognosis. Cancer Sci. 2009;100(8):1421–8.
Szabova L, Yin C, Bupp S, Guerin TM, Schlomer JJ, Householder DB, Baran ML, Yi M, Song Y, Sun W, et al. Perturbation of rb, p53, and Brca1 or Brca2 cooperate in inducing metastatic serous epithelial ovarian cancer. Cancer Res. 2012;72(16):4141–53.
Hill CG, Matyunina LV, Walker D, Benigno BB, McDonald JF. Transcriptional override: a regulatory network model of indirect responses to modulations in microRNA expression. BMC Syst Biol. 2014;8:36.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.
Lánczky A, Győrffy B. Web-based Survival Analysis Tool tailored for Medical Research (KMplot): development and implementation. J Med Internet Res. 2021;23(7):e27633.
Győrffy B. Discovery and ranking of the most robust prognostic biomarkers in serous ovarian cancer. Geroscience. 2023;45(3):1889–98.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res 2010, 38(Web Server issue):W214–220.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63.
Liu CJ, Hu FF, Xia MX, Han L, Zhang Q, Guo AY. GSCALite: a web server for gene set cancer analysis. Bioinformatics. 2018;34(21):3771–2.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7.
Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, Angell H, Fredriksen T, Lafontaine L, Berger A, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.
Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, Li B, Liu XS. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48(W1):W509–14.
Olalekan S, Xie B, Back R, Eckart H, Basu A. Characterizing the tumor microenvironment of metastatic ovarian cancer by single-cell transcriptomics. Cell Rep. 2021;35(8):109165.
Binnewies M, Pollack JL, Rudolph J, Dash S, Abushawish M, Lee T, Jahchan NS, Canaday P, Lu E, Norng M, et al. Targeting TREM2 on tumor-associated macrophages enhances immunotherapy. Cell Rep. 2021;37(3):109844.
Zhang X, Lan Y, Xu J, Quan F, Zhao E, Deng C, Luo T, Xu L, Liao G, Yan M, et al. CellMarker: a manually curated resource of cell markers in human and mouse. Nucleic Acids Res. 2019;47(D1):D721–8.
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
Greten FR, Grivennikov SI. Inflammation and Cancer: triggers, mechanisms, and consequences. Immunity. 2019;51(1):27–41.
Pan Y, Yu Y, Wang X, Zhang T. Tumor-Associated macrophages in Tumor Immunity. Front Immunol. 2020;11:583084.
Shi X, Wang J, Dai S, Qin L, Zhou J, Chen Y. Apolipoprotein C1 (APOC1): a Novel Diagnostic and Prognostic Biomarker for Cervical Cancer. Onco Targets Ther. 2020;13:12881–91.
Cui Y, Miao C, Hou C, Wang Z, Liu B. Apolipoprotein C1 (APOC1): a Novel Diagnostic and Prognostic Biomarker for Clear Cell Renal Cell Carcinoma. Front Oncol. 2020;10:1436.
Hilbert M, Kuzman P, Mueller WC, Meixensberger J, Nestler U. Feasibility of ApoC1 serum levels as tumor biomarker in glioblastoma patients: a pilot study. Sci Rep. 2022;12(1):16981.
Su WP, Sun LN, Yang SL, Zhao H, Zeng TY, Wu WZ, Wang D. Apolipoprotein C1 promotes prostate cancer cell proliferation in vitro. J Biochem Mol Toxicol. 2018;32(7):e22158.
Bosco MC. Macrophage polarization: reaching across the aisle? J Allergy Clin Immunol. 2019;143(4):1348–50.
Kerneur C, Cano CE, Olive D. Major pathways involved in macrophage polarization in cancer. Front Immunol. 2022;13:1026954.
Funding
This work was sponsored by the National Natural Science Foundation of China (81772767).
Author information
Authors and Affiliations
Contributions
DZ and XX proposed the idea of the article and supervised its completion. SY drafted and reviewed the article. WW performed the bioinformatics analysis. JD completed the experiments related to the article. SY and JD helped to interpret the data and prepare the figures. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethical approval and consent to participate
Prior to the operations, patients were informed about the procedures and provided their consent. Ethical approval for the study was obtained from the Medical Ethics Committee of Shanghai First People’s Hospital. Informed consent was taken from participants to participate in the study.
Competing interests
The authors declare no competing interests.
Consent for publication statement
Not applicable.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Yang, S., Du, J., Wang, W. et al. APOC1 is a prognostic biomarker associated with M2 macrophages in ovarian cancer. BMC Cancer 24, 364 (2024). https://doi.org/10.1186/s12885-024-12105-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12885-024-12105-z