Abstract
The complement system is readily triggered by the presence of damage-associated molecular patterns on the surface of tumor cells. The complement alternative pathway provides rapid amplification of the molecular stress signal, leading to complement cascade activation to deal with pathogens or malignant cells. Properdin is the only known positive regulator of the alternative pathway. In addition, properdin promotes the phagocytic uptake of apoptotic T cells by macrophages and dendritic cells without activating the complement system, thus, establishing its ability to recognize “altered-self”. Dysregulation of properdin has been implicated in substantial tissue damage in the host, and in some cases, chronic unresolved inflammation. A corollary of this may be the development of cancer. Hence, to establish a correlation between properdin presence/levels in normal and cancer tissues, we performed bioinformatics analysis, using Oncomine and UALCAN. Survival analyses were performed using UALCAN and PROGgeneV2 to assess if properdin can serve as a potential prognostic marker for human lung adenocarcinoma (LUAD), liver hepatocellular carcinoma (LIHC), cervical squamous cell carcinoma (CESC), and pancreatic adenocarcinoma (PAAD). We also analyzed levels of tumor-infiltrating immune cells using TIMER, a tool for characterizing immune cell composition in cancers. We found that in LUAD and LIHC, there was a lower expression of properdin in the tumors compared to normal tissues, while no significant difference was observed in CESC and PAAD. Survival analysis demonstrated a positive association between properdin mRNA expression and overall survival in all 4 types of cancers. TIMER analysis revealed that properdin expression correlated negatively with tumor purity and positively with levels of infiltrating B cells, cytotoxic CD8+ T cells, CD4+ helper T cells, macrophages, neutrophils and dendritic cells in LUAD, CESC and PAAD, and with levels of B cells, CD8+ T cells and dendritic cells in LIHC. Immunohistochemical analysis revealed that infiltrating immune cells were the most likely source of properdin in the tumor microenvironment. Thus, complement protein properdin shows promise as a prognostic marker in cancer and warrants further study.
Introduction
Properdin, a ~50 kDa glycoprotein found at 4–25 μg/ml level in plasma, is the only known positive regulator of the complement system (1, 2). The complement system is an integral part of the innate immune system. It consists of three pathways, based on the recognition subcomponents as well as pathogen-associated molecular patterns (PAMPs) or danger-associated molecular patterns (DAMPs) on activators such as pathogens and IgG/IgM (3). The complement classical pathway is triggered by C1q interaction with IgG/IgM bound to target ligands; the lectin pathway, a homolog of the classical pathway, is triggered by the binding of mannan-binding lectin (MBL) to the carbohydrate patterns on the pathogen surface (4, 5). The complement alternative pathway is triggered by the spontaneous hydrolysis of complement protein C3 into C3(H2O) (6, 7), which in the presence of surface PAMPs, binds to the zymogen form of Factor B to form the fluid-phase C3 proconvertase, C3(H2O)B (8). Factor B is then cleaved by the serine protease, Factor D, into fragments Ba and Bb; Bb remains bound to C3(H2O) non-covalently producing C3(H2O)Bb, which is stabilized by binding to properdin yielding C3(H2O)BbP, the alternative pathway fluid-phase C3 convertase (2). C3(H2O)Bb further cleaves additional molecules of C3 into C3a and C3b (2); C3b gets deposited on the pathogen cell surface, where it associates with Factor B to generate C3bB, the membrane-bound C3 proconvertase (4). Factor D then cleaves Factor B in the C3bB complex and binds to properdin to generate the membrane-bound C3 convertase, C3bBbP (4), which cleaves several C3 molecules, leading to C3b deposition on the pathogen surface (9). C3b deposition is essential for opsonisation while the chemoattractant, C3a, helps mount an inflammatory response (9). Additionally, C3b also associates with C3bBb and properdin, forming the C5 convertase, C3bBbP3b (10) that cleaves C5 into C5a (a potent chemoattractant) and C5b. C5b sequentially binds C6, C7, C8 and multiple molecules of C9 to form the membrane attack complex (MAC) (10). The MAC disrupts the cell membrane of the pathogen causing cell lysis (11–13).
The alternative pathway is tightly regulated by the interactions of Factor H, C3b and properdin. In the presence of polyanionic markers that are present on “Self” cells, Factor H binds to their surface and dissociates Bb from the proconvertase (C3bBb and C3(H2O)Bb), thereby accelerating the decay of the convertases (2, 14). Structural analysis of properdin has revealed that it can interfere with this decay activity of Factor H (15). Factor H also prevents the formation of new convertases by playing a role as a cofactor for Factor I-mediated cleavage of C3b to “inactive” C3b (iC3b) and C3(H2O) to iC3(H2O) (4, 16). This cofactor activity has also been shown to be limited by properdin (17, 18). Properdin has also been shown to trigger the alternative pathway by directly recruiting C3b or C3(H2O) and Factor B (19). Properdin is found in circulation as cyclic dimers, trimers and tetramers in a 26:54:20 ratio via head-to-tail association of its monomers; it can interact with surface-bound iC3b, C3b and C3bBb with increasing efficiency (1, 17, 20).
Properdin is locally produced in secondary granules of stimulated peripheral blood neutrophils, dendritic cells and T cells (21–25). In addition to the C3 convertase-stabilizing function, properdin can also act as a humoral pattern recognition receptor (PRR) for ligands such as lipopolysaccharide, acetylated low-density lipoprotein, and zymosan (2). Surface bound properdin has also been shown to activate the alternative pathway by recruiting C3b and factor B to form C3bBbP (19). Properdin deficiency, an X-chromosomal recessive trait, is associated with ~250-fold more susceptibility to Neisseria meningitidis infection (26). In the case of non-infectious diseases such as rheumatoid arthritis, patients show decreased levels of properdin, C3 and C5, but high levels of Ba and C3d in the synovial fluid (27). Recent demonstration of properdin’s complement-independent functions as a PRR and modulator of pro-inflammatory immune response have suggested its important roles beyond the upregulation of the complement alternative pathway (28, 29).
Macrophages derived from properdin knockout mice showed a shift to M2 phenotype, with increased type II cellular immune responses and decreased MHC class II expression (30). The lack of properdin has also been shown to protect the host from an aggravated inflammatory response, increasing the survival of properdin gene knock-out mice compared to wild type mice (31–34). In colitis models, properdin deficiency intensified colonic injury, and more severe renal injury was reported in the presence of a dysfunctional Factor H and low properdin levels (35, 36). Properdin has also been reported to promote the uptake of apoptotic cells by phagocytes (37, 38).
Since properdin is a modulator of pro-inflammatory response and ~80%–90% of the terminal pathway is activated by the alternative pathway, the role of properdin in malignancy has been recently examined. Analysis of 13,023 genes for mutations between primary breast tumors and normal tissue revealed mutations at a significant frequency in the properdin (CFP) gene (39), hinting at mutants contribution to tumor progression (40). The tumorigenic consequences of dysregulated alternative pathway has been reaffirmed by a recent study that showed that aging factor H-deficient mice spontaneously develop hepatic tumor (~50% of male mice) (41); an increased expression of Factor H mRNA seemed associated with improved survival in hepatocellular carcinoma patients (41).
Here, we analyzed properdin expression and its correlation with prognosis using databases (Oncomine, UALCAN and PROGgeneV2) in lung adenocarcinoma (LUAD), liver hepatocellular carcinoma (LIHC), cervical squamous cell carcinoma (CESC), and pancreatic adenocarcinoma (PAAD), and interrogated if properdin could be an independent marker for diagnosis and a surrogate marker of severity and overall survival. Moreover, we investigated the correlation of properdin with tumor-infiltrating immune cells (TIICs) in all 4 types of cancers via Tumor IMmune Estimation Resource (TIMER) and performed immunohistochemical analysis to identify the most likely sources of properdin in the tumor microenvironment.
Materials and Methods
Gene Expression and Survival Analysis
The expression levels of CFP gene in different carcinomas were analyzed using Oncomine (www.oncomine.org) and UALCAN (http://ualcan.path.uab.edu). Oncomine is a cancer microarray database and web-based data mining platform from the genome-wide expression analyses. UALCAN, in addition to being a web resource for analyzing cancer transcriptome data, estimates the effects of gene expression levels on the patient survival (42–44). We compared the differences in mRNA levels between carcinomas and normal tissues. The mRNA expression levels in the neoplastic tissues, as compared to healthy tissues, were obtained as parameters of p-value < 0.05, fold change > 2, and gene ranking in the top 10%. Further information about the datasets can be found in Table 1.
Table 1
| Datasets | Study Description | Experiment Type |
|---|---|---|
| Hou Lung | 91 non-small cell lung carcinoma and 65 adjacent normal lung samples were analyzed. Sample data includes age, sex, cancer sample site, and survival. | mRNA |
| Okayama Lung | 226 lung adenocarcinoma samples and 20 normal lung tissues were analyzed. Sample data includes EGFR mutation, KRAS mutation, EML4-ALK gene fusion, stage, recurrence, survival status, and others. | mRNA |
| Bhattacharjee Lung | 139 lung adenocarcinoma, 21 squamous cell lung carcinoma, 20 lung carcinoid tumor, six small cell lung carcinoma, and 17 normal lung samples were analyzed on Affymetrix U95A microarrays. Sample data includes type, age, M stage, max tumor percentage, N stage, primary/metastatic, recurrence, sex, site of metastasis, smoking rate (packs per year), stage, survival, and T stage. | mRNA |
| Su Lung | 91 non-small cell lung carcinoma and 65 adjacent normal lung samples were analyzed. Sample data includes age, sex, cancer sample site, and survival. | mRNA |
| TCGA Lung | 503 squamous cell lung carcinoma samples and 515 lung adenocarcinoma samples were analyzed on a custom Agilent microarray. Sample data includes age, TNM stage, sex, survival, smoking status, and others. This dataset is a combination of Lung Adenocarcinoma [LUAD] and Lung Squamous Cell Carcinoma [LUSC] data from the TCGA data portal and consists of Level 2 (processed) data. | mRNA |
| Roessler Liver | 225 hepatocellular carcinoma and 220 normal liver samples were analyzed. | mRNA |
| Roessler 2 Liver | 22 hepatocellular carcinoma, 19 normal liver, and two pooled normal liver tissue samples were analyzed. | mRNA |
| Wurmbach Liver | 75 liver tissue samples were analyzed on Affymetrix U133 Plus 2.0 microarrays. This set of samples represents the stepwise carcinogenic process from preneoplastic lesions to hepatocellular carcinoma, in the following order: Normal Liver (n=10), Cirrhotic Liver (n=13), Dysplastic Liver (n=17), and Hepatocellular Carcinoma (n=35). All non-normal samples were hepatitis virus C (HCV) positive, but negative for other known causal agents of hepatocellular cancer. Sample data includes HCC stage, tumor size, differentiation, vascular invasion, and satellites. | mRNA |
| TCGA Liver | 371 hepatocellular carcinoma, 56 paired normal blood, and 59 paired normal liver tissue samples were analyzed. Sample data includes sex, TNM stage, race/ethnicity, and others. This dataset consists of Level 3 data (segmented using CBS) from the TCGA data portal. The resulting segments were mapped to RefSeq gene coordinates as provided by UCSC. The samples were originally run on the Affymetrix SNP 6.0 platform. | mRNA |
| TCGA Cervix | 305 cervical squamous cell carcinoma, 95 paired blood derived normal, and three paired normal cervix tissue samples were analyzed. Sample data includes age, TNM stage, grade, FIGO stage, survival and others. This dataset consists of Level 3 data (segmented using CBS) from the TCGA data portal. The resulting segments were mapped to RefSeq gene coordinates as provided by UCSC. The samples were originally run on the Affymetrix SNP 6.0 platform. | mRNA |
| TCGA Pancreas | 178 pancreatic adenocarcinoma, 17 pancreatic ductal adenocarcinoma, one colloid carcinoma of the pancreas, 42 normal blood, and eight normal pancreas tissue samples were analyzed. This dataset consists of Level 3 data (segmented using CBS) from the TCGA data portal. The resulting segments were mapped to RefSeq gene coordinates as provided by UCSC. The samples were originally run on the Affymetrix SNP 6.0 platform. | mRNA |
Characteristics of the datasets used in bioinformatics analysis.
The prognostic significance of CFP mRNA expression with respect to survival in carcinomas was analyzed by UALCAN, using genomic data from “The Cancer Genome Atlas” (TCGA) to generate survival probability plots in the given carcinomas via Kaplan–Meier survival analysis (44). Additional survival analyses were performed with PROGgeneV2 (www.genomics.jefferson.edu/proggene/index.php) using data from TCGA (45, 46). Hazard ratio with a 95% confidence interval and log-rank p-value were also computed. In order to evaluate the prognostic effect of properdin, the four carcinomas we considered were LUAD, LIHC, CESC and PAAD.
Protein Expression Analysis
The differential protein expression of properdin between carcinoma and normal tissues was analyzed via UALCAN, which allowed protein expression analysis using the Clinical Proteomic Tumor Analysis Consortium (CPTAC) Confirmatory/Discovery dataset (47).
TIMER Database Analysis
TIMER is a comprehensive resource for systematic analysis of immune infiltrates across diverse cancer types (www.cistrome.shinyapps.io/timer) (48). TIMER applies a statistical method to infer the abundance of TIICs from the gene expression profiles using TCGA dataset (49). We analyzed, in LUAD, LIHC, CESC, PAAD, the correlation between CFP expression and the abundance of immune infiltrates, such as B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells, via gene modules (50–52). The correlation module generated the expression scatter plots between the CFP gene and defined genes of TIICs in selected carcinomas, together with the Spearman’s correlation and the estimated statistical significance. CFP was plotted on the x-axis and the related marker genes were represented on the y-axis as genes of TIICs. The gene expression level was displayed with log2 RNA-Seq by Expectation-Maximization (RSEM).
Immunohistochemical Analysis
Surgical normal tissue samples of lung, liver, uterine cervix and pancreas along with their malignant counterparts were selected for immunohistochemical analysis for properdin expression. Invasive malignant neoplasma specimens included the histotypes- LUAD, LIHC, CESC, and PAAD. The study was approved by the Institutional review board of the University of Palermo (09/2018). A specific informed consent was not required at the time of tissue sample collection for the immunohistochemical analysis of archival tissue sections since the patients were not identified and genetic analysis was not carried out. Immunohistochemistry (IHC) was performed using a polymer detection method. Briefly, tissue samples were fixed in 10% (v/v) neutral buffered formalin and then paraffin embedded. Tissue sections (4 μm thick) were deparaffinized and rehydrated. The antigen unmasking technique was carried out using Novocastra Epitope Retrieval Solutions, pH 9 (Leica Biosystems) in thermostatic bath at 98°C for 30 min. Sections were then brought to room temperature and washed in PBS. After neutralization of the endogenous peroxidase with 3% (v/v) H2O2 and Fc blocking by a specific protein block (Novocastra, Leica Biosystems), samples were incubated for 1 h at room temperature with mouse anti-human properdin (Clone C-4, 1:50, Santa Cruz) monoclonal antibody. Staining was revealed using Novolink Polymer Detection System (Novocastra, Leica Biosystems) and Romulin AEC Chromogen Kit (BioCare) as substrate chromogen, following the manufacturer’s instructions. Slides were counterstained with Harris Hematoxylin (BioOptica). Slides were analyzed under the Axio Scope A1 optical microscope (Zeiss) and microphotographs were collected through the Axiocam 503 color digital camera (Zeiss) using the Zen2 software.
Statistical Analysis
Survival curves were generated by UALCAN and/or PROGgeneV2. All results were displayed with p-values from a log-rank test. p-values < 0.05 were considered significant (53, 54). Similarly, in the case of Oncomine, the program provided the statistical significance of the data (p-values). In TIMER, the correlation of gene expression was evaluated by Spearman’s correlation and statistical significance, and the strength of the correlation was determined using the following guide for the absolute value: 0.00-0.19: very weak; 0.20-0.39: weak; 0.40-0.59: moderate; 0.60-0.79: strong; and 0.80-1.0: very strong.
Results
Lower Properdin Expression Correlates With Poor Prognosis in Lung Adenocarcinoma and Vice Versa
We initially evaluated the differences in the mRNA level of properdin in LUAD as compared to normal lung tissue, using the Oncomine platform via Hou’s, Okayama’s, Su’s and Bhattacharjee’s datasets (Table 1). We noted a significantly lower properdin mRNA expression in LUAD in comparison with healthy lung tissue counterpart (p < 0.05) (Figures 1A–D). Similar result was revealed by UALCAN tool analysis using the TCGA dataset (p < 0.05) (Figure 1E). UALCAN proteomics analysis also validated lower expression of properdin in LUAD, as compared to the corresponding healthy tissue (p < 0.05) (Figure 1F). As shown in Figure 2A, properdin mRNA expression positively correlated with the overall survival (OS) rate of LUAD patients (p < 0.05). Based on the TIMER analysis, CFP expression showed a significantly weak negative correlation with tumor purity (p < 0.05). It revealed a partial but significant and positive correlation with infiltrating levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells (p < 0.05) (Figure 2B). The strength of the correlation was very weak for CD8+ T cells, whereas it was relatively weaker for B cells, CD4+ T cells, macrophages, neutrophils and dendritic cells.
Figure 1
Figure 2
In lung IHC staining, with membrane as well as cytoplasmic labeling, properdin was found to be expressed mainly by the alveolar macrophages in a diffuse pattern; endothelial cells labeled mildly (Figures 3A, B). In LUAD, we observed properdin expression within the tumor-associated immune infiltrates (Figures 3C, D). Properdin positive cells were mainly represented by macrophages and somewhat minimally by granulocytes, present within the tumoral nests as well as in the intratumoral stromal microenvironment. Moreover, the density of properdin positive cells showed a variable distribution in the neoplastic tissue; in fact, it was higher in the inter-tumoral stroma compared to the one detected in the neoplastic proliferation. Thus, the immunohistochemical assessment confirmed that the immune cells were the main source of properdin.
Figure 3
Properdin Expression in Liver Hepatocellular Carcinoma Correlates Positively With Overall Survival and Immune Infiltration in the Tumor Microenvironment
Bioinformatics analysis of properdin mRNA expression was performed in the context of the LIHC using the Roessler, Roessler 2, Wurmbach and TCGA dataset. Similar to LUAD, a lower expression level of properdin was detected as compared to normal liver tissue (p < 0.05) (Figure 4). According to UALCAN tool, properdin mRNA expression positively associated with an OS rate in LIHC patients (p < 0.05) (Figure 5A). A significant and moderate negative correlation between properdin expression and tumor purity was also observed (p < 0.05) (Figure 5B). It showed partial but significant and positive correlations with infiltrating levels of B cells, CD8+ T cells and dendritic cells (p < 0.05). The strength of the correlation was very weak for B cells and dendritic cells, whereas it was weak for CD8+ T cells. No significant correlation between the infiltrating levels of CD4+ T cells, macrophages and neutrophils and CFP expression was noted in LIHC.
Figure 4
Figure 5
The IHC staining of liver parenchyma for properdin highlighted a variable pattern: ~30%–40% of hepatocytes exhibited cytoplasmic expression of properdin. A few immune infiltrates normally present within the portal space and in the portal stroma also showed properdin positivity (Figures 6A, B). In LIHC, immunohistochemical staining revealed properdin expression in the TIICs within the stromal microenvironment as well as a few neoplastic components (~2% of total tumor cells) (Figures 6C, D).
Figure 6
Properdin Is Expressed in Cervical Squamous Cell Carcinoma and Normal Cervical Tissues at Similar Levels
No significant difference regarding properdin expression level in CECS was detected as compared to normal cervical tissues using different datasets (Data Not Shown). However, among CECS patients, a higher expression of properdin positively correlated with an OS (p < 0.05) (Figure 7A). According to the TIMER tool, CFP gene expression moderately, but significantly and negatively, correlated with the tumor purity (p < 0.05). It revealed a partial but significant and positive correlation with the levels of infiltrated B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells (p < 0.05) (Figure 7B). The strength of the correlation was very weak for CD8+ T cells, whereas it was weak for B cells, CD4+ T cells, macrophages, neutrophils and dendritic cells. In the uterine cervical tissue, immunohistochemical assessment for the presence of properdin protein revealed a high expression in both the normal as well as the tumoral tissues (Figure 8). Properdin appeared to be diffusely expressed not only within the infiltrating cells (macrophages and few cells of lymphoid morphology), but also within the endothelial cells of small vessels located in the subepithelial connective tissue.
Figure 7
Figure 8
Properdin in Pancreatic Adenocarcinoma
No significant difference regarding properdin expression level in PAAD was detected as compared to normal pancreas using various datasets (Data Not Shown). According to UALCAN tool, no correlation between properdin mRNA expression and OS rate was observed in patients with PAAD (Figure 9A). However, when we used PROGgeneV2 tool, properdin mRNA expression was found to be positively associated with an OS rate in PAAD patients (p < 0.05) (Figure 9B). With respect to immune infiltration, CFP gene expression showed a partial but significant and negative relationship with the tumor purity (p < 0.05). This means that properdin is unlikely to be expressed by tumor cells but by the immune infiltrates. It suggested significant positive partial correlations with infiltrating levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells (p < 0.05) (Figure 9C). The strength of the correlation was weak for B cells, CD8+ T cells and macrophages, whereas it was found to be moderate for CD4+ T cells, neutrophils and dendritic cells. IHC staining of the healthy pancreatic tissue revealed that ~10% of acinar cells of the exocrine part of pancreas had a high expression of properdin (strong cytoplasmic labeling) (Figure 10A). The immune infiltrate, consisting of scant macrophages present in the connective tissue contexture also stained for properdin (Figure 10B). A marked expression of properdin by the infiltrating immune cells, mostly consisting of macrophages and granulocytes, was evident; however, we observed no properdin expression in the neoplastic subclones (Figures 10C, D).
Figure 9
Figure 10
Discussion
The innate immunity plays a critical role in the protection against cancer. The complement system is an important arsenal of the innate immune system against pathogens as well as cancer. The complement system is initiated in response to the tumor-associated antigens, and leads to increased deposition of complement activation fragments on the surface of the tumor. C1q, for instance, has been shown to bind to phospholipids of lung tumor cell lines and evidence of activation of the classical pathway in various cancer such as lung cancer, thyroid cancer, astrocytoma, lymphoma, leukemia, and oropharyngeal cancer have been reported (55, 56). C1q has recently been shown to be involved in a range of pathophysiological functions that are not dependent on complement activation. In fact, it is expressed in the microenvironment of various types of human tumors, where it can exert a protective or a harmful effect on cancer progression (57–60). Similarly, in cases of acute lymphoblastic leukemia, Burkitt lymphoma or multiple myeloma, dysregulation of the complement alternative pathway has been reported (61). Exaggerated or aberrant complement activation in the tumor microenvironment has been shown to increase inflammation/infiltration via C5a. Anti-tumor immune responses are stifled by the C5a-driven recruitment of myeloid-derived suppressor cells, which support tumorigenesis by upregulating anti-inflammatory (inhibitory) molecules such as Programmed Death-Ligand 1 (PD-L1), transforming growth factor-β (TGF-β) and IL-10 (62–64). C3a has also been implicated in tumorigenesis; C3aR−/− mice show improved disease outcome in melanoma (65). Factor H gene-deficient mice appear to develop spontaneous hepatocellular carcinoma with aging (41). An in silico analysis appeared to suggest a protective role of properdin in breast cancer by interacting with glycosaminoglycan (GAG) structures on the tumor cell surface (39, 41, 66).
In the current study, we analyzed the prognostic value of properdin, the positive regulator of the complement alternative pathway, in four types of cancer, namely LUAD, LIHC, CESC, and PAAD. The choice of the four tumors was made based on the prognostic effect of properdin, which was assessed via UALCAN using TCGA RNA-sequencing and patients’ clinical data from 33 different cancer types, including several metastatic tumors. Of the 33 cancers, properdin had a significant prognostic effect in only 4 neoplasms. Subsequently, CFP gene expression in various neoplastic tissues was analyzed using Oncomine and UALCAN, and compared with their healthy counterparts. Lower levels of CFP expression were observed in both LUAD and LIHC tissues compared to their respective normal tissues. No significant difference was observed in the expression levels in the case of CESC and PAAD. Proteomic analysis of properdin levels in LUAD was also analyzed by UALCAN and found to be lower in the neoplastic tissue compared to the healthy tissues. However, currently, proteomic analysis of the other three cancers are not available on UALCAN. Kaplan–Meier survival analysis of LUAD, LIHC and CESC by UALCAN revealed that properdin mRNA expression positively correlated with an OS rate of the patients. In the case of PAAC, UALCAN survival analysis did not find any correlation between the mRNA expression levels and OS; however, analysis by PROGgeneV2 revealed a positive correlation between mRNA expression levels and OS. This variation can be attributed to the difference in the way the tools define patient cohorts, as both tools use the same TCGA Pancreas dataset. UALCAN divides patients into a high expression of properdin versus low/medium expression. PROGgeneV2 divides patients at the median of gene expression between the ones with low expression and those with high expression. This is a limitation of the study involving bioinformatics tools that can be used to evaluate the prognostic implications of genes in various cancers. Thus, using the same cohort of patients (for example, TCGA Pancreas dataset), the results on the prognostic effect can differ, depending on whether the court is divided by the median of expression or not.
It is highly plausible that the expression of the CFP gene during tumorigenesis may lead to properdin export to the cell surface (66), which may then promote phagocytosis by macrophages and dendritic cells (37, 38), in addition to recruiting immune cells. Thus, the neoplastic cells have evolved to express low levels of properdin, and hence, higher expression of the protein leads to increased overall survival. Since TIICs are independent predictors of sentinel lymph node status and survival in cancer (67, 68), we investigated whether CFP expression correlated with immune infiltration levels in various cancer types. Tumor purity is defined as the proportion of cancer cells in the tumor tissue or the percentage of cancer cells in a solid tumor sample. CFP expression negatively correlated with the tumor purity in all cancers. This suggests that the expression of properdin likely occurs by the immune infiltrates and not by the actual neoplastic cells. Infiltrating levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells positively correlated in the case of LUAD, CESC and PAAD. In the case of LIHC, a positive correlation with infiltrating levels of B cells, CD8+ T cells and dendritic cells was observed, while no significant correlation was observed between the expression of CFP and the levels of infiltrating CD4+ T cells, macrophages and neutrophils. These findings suggest that CFP plays a specific role in immune infiltration in our carcinomas, in particular as regard CD4+ T cells, neutrophils and dendritic cells. Therefore, CFP could exert an important function in recruitment and regulation of immune infiltrating cells in several carcinomas.
Interestingly, preliminary studies conducted in properdin knock-out mice showed a difference in chemokine levels and cell mobility in a tumor model (CM Stover, unpublished Data). The role of immune cell infiltration in the tumor microenvironment and its consequences on tumor development is an intense area of research: immune infiltration and cell phenotype polarization can be a double-edged sword. It can promote or impede tumor growth, and hence, overall survival. Given a positive correlation with immune infiltration, the properdin expression levels reflect the presence of immune cells that are less immunosuppressive but can engage in the anti-tumor response. Another possibility is that properdin acts as a balancing factor between macrophage polarization between pro-inflammatory M1 and anti-inflammatory M2 phenotypes (30). NK cells, a key player in the anti-tumor immunity, have recently been shown to bind properdin through NKp46 receptor (69). However, how properdin can impact on NK cell-mediated properties in the tumor microenvironment is unclear.
Anti-tumor immunity has been recently categorized into mainly three phenotypes: immune-desert (immune ignorance), immune-excluded (inhibitory chemokine/soluble factor-mediated), and inflamed (where infiltration is immunosuppressive and/or dysregulated) (70). The immune cell infiltration here includes regulatory T cells, myeloid-derived suppressor cells, suppressor B cells and fibroblasts. Tumor-infiltrating CD8+ T lymphocytes can also be aberrantly active and immunologically exhausted. In a recent study, Troiano et al. identified a specific subgroup of squamous cell carcinoma of the oral tongue with poor prognosis based on the density of tumor‐infiltrating lymphocytes and localization (71). Thus, the tumor immune phenotyping can be a better approach to stratify patients and decide on precision medicine.
Understanding how the level of properdin expressed in the tumor-infiltrating immune cells alters the tumor microenvironment and manipulates immune cells versus tumor cells should be of great relevance for immunotherapeutic development. Hence, further study to examine the differential expression and distribution of properdin along with immune cells in the neoplastic tissues is urgently required. As properdin has also been shown to function in a complement-independent manner (27, 28, 37), experiments in the absence of the complement system and with an individual pathway blocked must also be undertaken. In conclusion, our in-silico analysis highlights a possible role of properdin as a novel marker for tumor prognosis and increased survival in a range of cancers.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary materials; further inquiries can be directed to the corresponding authors.
Author contributions
AlM, PV, BB, AnM, CA, and GR generated and analyzed the data. SA, HK, and CS interpreted and contributed to the draft. RB and UK led the work. AnM and PV did the first draft. RB and UK finalized the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors acknowledge the funding from the Italian Ministry of Health (RC 20/2016- Institute for Maternal and Child Health IRCCS Burlo Garofolo, Trieste, Italy) and International Scientific Partnership Programme (ISPP) at the King Saud University, Riyadh for funding via ISPP-145.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
PangburnMK. Analysis of the natural polymeric forms of human properdin and their functions in complement activation. J Immunol (1989) 142:202–7.
2
ChenJYCortesCFerreiraVP. Properdin: A multifaceted molecule involved in inflammation and diseases. Mol Immunol (2018) 102:58–72. doi: 10.1016/j.molimm.2018.05.018
3
PedersenDVGadebergTAFThomasCWangYJoramNJensenRKet al. Structural basis for properdin oligomerization and convertase stimulation in the human complement system. Front Immunol (2019) 10:2007. doi: 10.3389/fimmu.2019.02007
4
KouserLAbdul-AzizMNayakAStoverCMSimRBKishoreU. Properdin and Factor H: Opposing Players on the Alternative Complement Pathway “See-Saw.” Front Immunol (2013) 4:93. doi: 10.3389/fimmu.2013.00093
5
KishoreUReidKBM. Modular organization of proteins containing C1q-like globular domain. Immunopharmacology (1999) 42:15–21. doi: 10.1016/S0162-3109(99)00011-9
6
AbbasALichtmanAPillaiS. Cellular and Molecular Immunology. 6th ed.Philadelphia, PA, USA: Saunders/Elsevier (2010). p. 267–301.
7
Muller-EberhardH. Molecular Organization And Function Of The Complement System. Annu Rev Biochem (1988) 57:321–47. doi: 10.1146/annurev.biochem.57.1.321
8
NorisMRemuzziG. Overview of complement activation and regulation. Semin Nephrol (2013) 33:479–92. doi: 10.1016/j.semnephrol.2013.08.001
9
BarnumSScheinT. The Complement FactsBook.London, United Kingdom:Elsevier (2018) pp. 480. doi: 10.1016/C2015-0-06595-9
10
ZwarthoffSABerendsETMMolSRuykenMAertsPCJózsiMet al. Functional Characterization of Alternative and Classical Pathway C3/C5 Convertase Activity and Inhibition Using Purified Models. Front Immunol (2018) 9:1691. doi: 10.3389/fimmu.2018.01691
11
MorganBP. The membrane attack complex as an inflammatory trigger. Immunobiology (2016) 221:747–51. doi: 10.1016/j.imbio.2015.04.006
12
TeglaCACudriciCPatelSTrippeRRusVNiculescuFet al. Membrane attack by complement: the assembly and biology of terminal complement complexes. Immunol Res (2011) 51:45. doi: 10.1007/s12026-011-8239-5
13
NautaAJDahaMRTijsmaOVan De WaterBTedescoFRoosA. The membrane attack complex of complement induces caspase activation and apoptosis. Eur J Immunol (2002) 32:783–92. doi: 10.1002/1521-4141(200203)32:3<783
14
FerreiraVPPangburnMKCortésC. Complement control protein factor H: The good, the bad, and the inadequate. Mol Immunol (2010) 47:2187–97. doi: 10.1016/j.molimm.2010.05.007
15
AlcorloMTortajadaARodriguez de CordobaSLlorcaO. Structural basis for the stabilization of the complement alternative pathway C3 convertase by properdin. Proc Natl Acad Sci (2013) 110:13504–9. doi: 10.1073/pnas.1309618110
16
SimRBDayAJMoffattBEFontaineM. Complement factor I and cofactors in control of complement system convertase enzymes. Methods Enzymol (1993) 223:13–35. doi: 10.1016/0076-6879(93)23035-l
17
FarriesTCLachmannPJHarrisonRA. Analysis of the interactions between properdin, the third component of complement (C3), and its physiological activation products. Biochem J (1988) 252:47–54. doi: 10.1042/bj2520047
18
MedicusRGGötzeOMüller-EberhardHJ. Alternative pathway of complement: Recruitment of precursor properdin by the labile C3/C5 convertase and the potentiation of the pathway. J Exp Med (1976) 144:1076–93. doi: 10.1084/jem.144.4.1076
19
HourcadeDE. The role of properdin in the assembly of the alternative pathway C3 convertases of complement. J Biol Chem (2006) 281:2128–32. doi: 10.1074/jbc.M508928200
20
SmithCAPangburnMKVogelCWMuller-EberhardHJ. Molecular architecture of human properdin, a positive regulator of the alternative pathway of complement. J Biol Chem (1984) 259:4582–8.
21
DixonKOO’FlynnJKlar-MohamadNDahaMRvan KootenC. Properdin and factor H production by human dendritic cells modulates their T-cell stimulatory capacity and is regulated by IFN-γ. Eur J Immunol (2017) 47:470–80. doi: 10.1002/eji.201646703
22
O’FlynnJDixonKKlar MohamadNDahaMRvan KootenC. Dendritic cells express alternative pathway components properdin and factor H: Differential regulation by Interferons and IL-27. Mol Immunol (2013) 56:263. doi: 10.1016/j.molimm.2013.05.069
23
WirthmuellerUDewaldBThelenMSchäferMKStoverCWhaleyKet al. Properdin, a positive regulator of complement activation, is released from secondary granules of stimulated peripheral blood neutrophils. J Immunol (1997) 158(9):4444–51.
24
SchwaebleWDippoldWGSchäferMKPohlaHJonasDLuttigBet al. Properdin, a positive regulator of complement activation, is expressed in human T cell lines and peripheral blood T cells. J Immunol (1993) 151(5):2521–8.
25
VuagnatBBMachJ-PLe DoussalJ-M. Activation of the alternative pathway of human complement by autologous cells expressing transmembrane recombinant properdin. Mol Immunol (2000) 37:467–78. doi: 10.1016/S0161-5890(00)00069-9
26
FijenCAPVan Den BogaardRSchipperMMannensMSchlesingerMNordinFGet al. Properdin deficiency: Molecular basis and disease association. Mol Immunol (1999) 36:863–7. doi: 10.1016/S0161-5890(99)00107-8
27
El-GhobareyAFWhaleyK. Alternative pathway complement activation in rheumatoid arthritis. J Rheumatol (1980) 7(4):453–60.
28
Al-MozainiMATsolakiAGAbdul-AzizMAbozaidSMAl-AhdalMNPathanAAet al. Human properdin modulates macrophage: Mycobacterium bovis BCG Interaction via Thrombospondin Repeats 4 and 5. Front Immunol (2018) 9:533. doi: 10.3389/fimmu.2018.00533
29
KouserLPaudyalBKaurAStenbeckGJonesLAAbozaidSMet al. Human properdin opsonizes nanoparticles and triggers a potent pro-inflammatory response by macrophages without involving complement activation. Front Immunol (2018) 9:131. doi: 10.3389/fimmu.2018.00131
30
Al-RayahiIAMBrowningMJStoverC. Tumour cell conditioned medium reveals greater M2 skewing of macrophages in the absence of properdin. Immunity Inflammation Dis (2017) 5:68–77. doi: 10.1002/iid3.142
31
WangYMiwaTDucka-KokalariBRedaiIGSatoSGullipalliDet al. Properdin Contributes to Allergic Airway Inflammation through Local C3a Generation. J Immunol (2015) 195:1171–81. doi: 10.4049/jimmunol.1401819
32
KimuraYZhouLMiwaTSongWC. Genetic and therapeutic targeting of properdin in mice prevents complement-mediated tissue injury. J Clin Invest (2010) 120:3545–54. doi: 10.1172/JCI41782
33
DimitrovaPIvanovskaNBelenskaLMilanovaVSchwaebleWStoverC. Abrogated RANKL expression in properdin-deficient mice is associated with better outcome from collagen-antibody-induced arthritis. Arthritis Res Ther (2012) 14:R173. doi: 10.1186/ar3926
34
ZhouHFYanHStoverCMFernandezTMRodriguez De CordobaSSongWCet al. Antibody directs properdin-dependent activation of the complement alternative pathway in a mouse model of abdominal aortic aneurysm. Proc Natl Acad Sci USA (2012) 109:E415–22. doi: 10.1073/pnas.1119000109
35
JainUMidgenCASchwaebleWJStoverCMStadnykAW. Properdin Regulation of Complement Activation Affects Colitis in Interleukin 10 Gene–Deficient Mice. Inflammation Bowel Dis (2015) 21:1519–28. doi: 10.1097/MIB.0000000000000398
36
RusevaMMVernonKALesherAMSchwaebleWJAliYMBottoMet al. Loss of Properdin Exacerbates C3 Glomerulopathy Resulting from Factor H Deficiency. J Am Soc Nephrol (2013) 24:43–52. doi: 10.1681/ASN.2012060571
37
KemperCMitchellLMZhangLHourcadeDE. The complement protein properdin binds apoptotic T cells and promotes complement activation and phagocytosis. Proc Natl Acad Sci USA (2008) 105:9023–8. doi: 10.1073/pnas.0801015105
38
XuWBergerSPTrouwLAde BoerHCSchlagweinNMutsaersCet al. Properdin Binds to Late Apoptotic and Necrotic Cells Independently of C3b and Regulates Alternative Pathway Complement Activation. J Immunol (2008) 180:7613–21. doi: 10.4049/jimmunol.180.11.7613
39
SjöblomTJonesSWoodLDParsonsDWLinJBarberTDet al. The consensus coding sequences of human breast and colorectal cancers. Sci (80- ) (2006) 314:268–74. doi: 10.1126/science.1133427
40
FusterMMEskoJD. The sweet and sour of cancer: glycans as novel therapeutic targets. Nat Rev Cancer (2005) 5:526–42. doi: 10.1038/nrc1649
41
LaskowskiJRennerBPickeringMCSerkovaNJSmith-JonesPMClambeyETet al. Complement factor H–deficient mice develop spontaneous hepatic tumors. J Clin Invest (2020) 130(8):4039–54. doi: 10.1172/JCI135105
42
RhodesDRKalyana-SundaramSMahavisnoVVaramballyRYuJBriggsBBet al. Oncomine 3.0: Genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia (2007) 9:166–80. doi: 10.1593/neo.07112
43
RhodesDRYuJShankerKDeshpandeNVaramballyRGhoshDet al. ONCOMINE: A Cancer Microarray Database and Integrated Data-Mining Platform. Neoplasia (2004) 6:1–6. doi: 10.1016/s1476-5586(04)80047-2
44
ChandrashekarDSBashelBBalasubramanyaSAHCreightonCJPonce-RodriguezIChakravarthiBVSKet al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia (2017) 19:649–58. doi: 10.1016/j.neo.2017.05.002
45
GoswamiCNakshatriH. PROGgene: gene expression based survival analysis web application for multiple cancers. J Clin Bioinform (2013) 3:22. doi: 10.1186/2043-9113-3-22
46
GoswamiCPNakshatriH. PROGgeneV2: Enhancements on the existing database. BMC Cancer (2014) 14:970. doi: 10.1186/1471-2407-14-970
47
ChenFChandrashekarDSVaramballySCreightonCJ. Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers. Nat Commun (2019) 10:5679. doi: 10.1038/s41467-019-13528-0
48
LiTFanJWangBTraughNChenQLiuJSet al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res (2017) 77:e108–10. doi: 10.1158/0008-5472.CAN-17-0307
49
LiBSeversonEPignonJ-CZhaoHLiTNovakJet al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol (2016) 17:174. doi: 10.1186/s13059-016-1028-7
50
SiemersNOHollowayJLChangHChasalowSDRoss-MacDonaldPBVolivaCFet al. Genome-wide association analysis identifies genetic correlates of immune infiltrates in solid tumors. PloS One (2017) 12:e0179726. doi: 10.1371/journal.pone.0179726
51
DanaherPWarrenSDennisLD’AmicoLWhiteADisisMLet al. Gene expression markers of Tumor Infiltrating Leukocytes. J Immunother Cancer (2017) 5:18. doi: 10.1186/s40425-017-0215-8
52
AranDSirotaMButteAJ. Systematic pan-cancer analysis of tumour purity. Nat Commun (2015) 6:8971. doi: 10.1038/ncomms9971
53
MangognaABelmonteBAgostinisCRicciGGulinoAFerraraIet al. Pathological Significance and Prognostic Value of Surfactant Protein D in Cancer. Front Immunol (2018) 9:1748. doi: 10.3389/fimmu.2018.01748
54
KumarJMurugaiahVSotiriadisGKaurAJeyaneethiJSturnioloIet al. Surfactant Protein D as a Potential Biomarker and Therapeutic Target in Ovarian Cancer. Front Oncol (2019) 9:542. doi: 10.3389/fonc.2019.00542
55
BerraondoPMinuteLAjonaDCorralesLMeleroIPioR. Innate immune mediators in cancer: between defense and resistance. Immunol Rev (2016) 274:290–306. doi: 10.1111/imr.12464
56
AjonaDPajaresMJCorralesLPerez-GraciaJLAgorretaJLozanoMDet al. Investigation of Complement Activation Product C4d as a Diagnostic and Prognostic Biomarker for Lung Cancer. JNCI J Natl Cancer Inst (2013) 105:1385–93. doi: 10.1093/jnci/djt205
57
MangognaAAgostinisCBonazzaDBelmonteBZacchiPZitoGet al. Is the Complement Protein C1q a Pro- or Anti-tumorigenic Factor? Bioinformatics Analysis Involving Human Carcinomas. Front Immunol (2019) 10:865. doi: 10.3389/fimmu.2019.00865
58
MangognaABelmonteBAgostinisCZacchiPIacopinoDGMartoranaAet al. Prognostic Implications of the Complement Protein C1q in Gliomas. Front Immunol (2019) 10:2366. doi: 10.3389/fimmu.2019.02366
59
AgostinisCVidergarRBelmonteBMangognaAAmadioLGeriPet al. Complement Protein C1q Binds to Hyaluronic Acid in the Malignant Pleural Mesothelioma Microenvironment and Promotes Tumor Growth. Front Immunol (2017) 8:1559. doi: 10.3389/fimmu.2017.01559
60
BullaRTripodoCRamiDLingGSAgostinisCGuarnottaCet al. C1q acts in the tumour microenvironment as a cancer-promoting factor independently of complement activation. Nat Commun (2016) 7:1–11. doi: 10.1038/ncomms10346
61
ReisESMastellosDCRicklinDMantovaniALambrisJD. Complement in cancer: Untangling an intricate relationship. Nat Rev Immunol (2018) 18:5–18. doi: 10.1038/nri.2017.97
62
WarrenerPBonnellJMustelinTFungMKolbeckR. Complement C5a induces PD-L1 expression and acts in synergy with LPS through Erk1/2 and JNK signaling pathways. Sci Rep (2016) 6:33346. doi: 10.1038/srep33346
63
PiaoCCaiLQiuSJiaLSongWDuJ. Complement 5a enhances hepatic metastases of colon cancer via monocyte chemoattractant protein-1-mediated inflammatory cell infiltration. J Biol Chem (2015) 290:10667–76. doi: 10.1074/jbc.M114.612622
64
MarkiewskiMMDeAngelisRABenenciaFRicklin-LichtsteinerSKKoutoulakiAGerardCet al. Modulation of the antitumor immune response by complement. Nat Immunol (2008) 9:1225–35. doi: 10.1038/ni.1655
65
NabizadehJAMantheyHDSteynFJChenWWidiapradjaAMd AkhirFNet al. T Cell Responses + Neutrophil and CD4 to Melanoma Tumorigenesis by Inhibiting The Complement C3a Receptor Contributes. Am Assoc Immnol (2020) 196(11):4783–92. doi: 10.4049/jimmunol.1600210
66
KemperCHourcadeDE. Properdin: New roles in pattern recognition and target clearance. Mol Immunol (2008) 45:4048–56. doi: 10.1016/j.molimm.2008.06.034
67
OhtaniH. Focus on TILs: prognostic significance of tumor infiltrating lymphocytes in human colorectal cancer. Cancer Immun (2007) 7:4.
68
AzimiFScolyerRARumchevaPMoncrieffMMuraliRMcCarthySWet al. Tumor-Infiltrating Lymphocyte Grade Is an Independent Predictor of Sentinel Lymph Node Status and Survival in Patients With Cutaneous Melanoma. J Clin Oncol (2012) 30:2678–83. doi: 10.1200/JCO.2011.37.8539
69
Narni-MancinelliEGauthierLBaratinMGuiaSFenisADeghmaneAEet al. Complement factor P is a ligand for the natural killer cell–activating receptor NKp46. Sci Immunol (2017) 2(10):eaam9628. doi: 10.1126/sciimmunol.aam9628
70
ChenDSMellmanI. Elements of cancer immunity and the cancer-immune set point. Nature (2017) 541(7637):321–30. doi: 10.1038/nature21349
71
TroianoGRubiniCTogniLCaponioVCAZhurakivskaKSantarelliAet al. The immune phenotype of tongue squamous cell carcinoma predicts early relapse and poor prognosis. Cancer Med (2020) 9(22):8333–44. doi: 10.1002/cam4.3440
Summary
Keywords
properdin, innate immunity, bioinformatics, cancer, complement, prognosis
Citation
Mangogna A, Varghese PM, Agostinis C, Alrokayan SH, Khan HA, Stover CM, Belmonte B, Martorana A, Ricci G, Bulla R and Kishore U (2021) Prognostic Value of Complement Properdin in Cancer. Front. Immunol. 11:614980. doi: 10.3389/fimmu.2020.614980
Received
07 October 2020
Accepted
08 December 2020
Published
19 January 2021
Volume
11 - 2020
Edited by
Lorenzo Lo Muzio, University of Foggia, Italy
Reviewed by
Marco Mascitti, Marche Polytechnic University, Italy; Kushagra Bansal, Harvard Medical School, United States; Maria Eleonora Bizzoca, University of Foggia, Italy
Updates
Copyright
© 2021 Mangogna, Varghese, Agostinis, Alrokayan, Khan, Stover, Belmonte, Martorana, Ricci, Bulla and Kishore.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Roberta Bulla, rbulla@units.it; Uday Kishore, uday.kishore@brunel.ac.uk; ukishore@hotmail.com
This article was submitted to Molecular Innate Immunity, a section of the journal Frontiers in Immunology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.