Expression pattern of genes associated with tumor microenvironment in prostate cancer
Summary. Aim: To assess relative expression (RE) levels of CAF-, TAM-specific, immune defense-associated genes in prostate tumors and to show correlation of RE with clinical, pathological and molecular characteristics, with the aim to define clinically significant specific alterations in a gene expression pattern. Methods: RE of 23 genes was analyzed by a quantitative polymerase chain reaction in 37 freshly frozen samples of prostate cancer tissues of a different Gleason score (GS) and at various tumor stages, compared with RE in 37 paired conventionally normal prostate tissue (CNT) samples and 20 samples of prostate adenomas. Results: Differences in RE were shown for 11 genes out of 23 studied, when tumor samples were compared with corresponding CNTs. 7 genes, namely ACTA2, CXCL14, CTGF, THY1, FAP, CD163, CCL17 were upregulated in tumors. 4 genes, namely CCR4, NOS2A, MSMB, IL1R1 were downregulated in tumors. 14 genes demonstrated different RE in TNA at different stages: CXCL12, CXCL14, CTGF, FAP, HIF1A, THY1, CCL17, CCL22, CCR4, CD68, CD163, NOS2A, CTLA4, IL1R1. RE changes of 9 genes — CXCL12, CXCL14, HIF1A, CCR4, CCL17, NOS2A, CTLA4, IL1R1, IL2RA — were found in tumors with different GS. Moreover, 9 genes showed differences in RE in TNA, dependently on the presence or absence of the TMPRSS2/ERG fusion and 7 genes showed differences in RE of groups with differential PTEN expression. Significant correlations were calculated between RE of 9 genes in adenocarcinomas and the stage, and GS; also, between RE of 2 genes and the fusion presence; and between RE of 4 genes and PTEN expression. Conclusions: Several gene expression patterns were identified that correlated with the GS, stage and molecular characteristics of tumors, i.e. presence of the TMPRSS2/ERG fusion and alterations in PTEN expression. These expression patterns can be used for molecular profiling of prostate tumors, with the aim to develop personalized medicine approaches. However, the proposed profiling requires a more detailed analysis and a larger cohort of patients with prostate tumor.
Submitted: July 05, 2018.
*Correspondence: E-mail: email@example.com
Abbreviations used: A — sample of prostate adenomas; CAF — cancer-associated fibroblasts; CNT/N — conventionally normal prostate tissue; CPC — clinical and pathological characteristics; FDR — false discovery rate; GS — Gleason score, IAG — immune-associated genes; PCR — polymerase chain reaction; RE — relative gene expression; st — stage; T – sample of prostate cancer, adenocarcinoma; TAM — tumor-associated macrophages; TNM — International System of Classiﬁcation of Tumors, based on tumor-node-metastasis.
It is well known that tumor initiation and progression, including prostate cancer, depend on genetic and epigenetic changes, leading to transcriptomic and proteomic disorders in cells upon malignant transformation [1–3]. Also, many parameters are altered, that is hierarchically regulated at different levels (organism, organs and systems, tissues and cells) [4–7]. All these processes constitute a common phenomenon named as a tumor-host interaction. One of the disorders in the organism that contributes to the appearance of tumors is chronic inflammation [8, 9]. It could be of different etiologies, such as viral, bacterial, parasitic, chemical or radiation-induced and that of unidentified etiology [8, 10]. However, regardless of its origin, both, the epithelial cells of the organs and the cells of the stroma (fibroblasts, macrophages, endothelial cells, etc.) are involved. Later, all these cells become the tumor microenvironment , changing their properties from tumor-suppressive to tumor-stimulating and supporting [12, 13]. For example, fibroblasts are transformed into cancer-associated fibroblasts (CAF) [14, 15] and macrophages change their phenotype to tumor-associated macrophages (TAM) [16–19]. Tumors are filled with tumor-infiltrating leukocytes which suppress immune defense and reactions [20, 21].
We have to mention, that the tumor growth is affected by antimicrobial proteins, secreted by epithelial cells of the prostate. One of such proteins is a product of the MSMB gene . Furthermore, elements of a major histocompatibility complex have as immunostimulating as immunosuppressive properties. One of the immunosuppressive glycoprotein is HLA-G .
Immunohistochemistry, the new generation sequencing and a real-time polymerase chain reaction (PCR) have been used to identify markers, typical for different types of cells in a tumor microenvironment [15, 16, 20]. We selected a number of genes specific for fibroblasts and CAF, macrophages and TAM and various subpopulations of T-lymphocytes, NK cells. Also, markers of inflammation and nonspecific immune reactions were chosen. The levels of relative gene expression (RE) were assessed and putative correlation with clinical and pathological characteristics (CPC) and molecular features was analyzed. Further, we wanted to monitor the specific alterations in an expression pattern of selected genes that are clinically significant, for the possible application in the molecular profiling of prostate cancers.
MATERIALS AND METHODS
Prostate tissue collection. Tissues of prostate cancer (T) and the paired conventionally normal prostate tissues (CNT, or N from a site opposite to cancer) were frozen in liquid nitrogen directly after surgery. All samples were collected at National Cancer Institute (Kyiv, Ukraine). Benign prostate tumors (prostate adenoma samples — A) were collected with the same procedure at the Institute of Urology (Kyiv, Ukraine) after radical prostatectomy. The samples were collected in accordance with the Declaration of Helsinki and the guidelines, issued by the Ethic Committee of the Institute of Urology of the National Academy of Medical Sciences of Ukraine and the National Cancer Institute and the Ethic Committee of the Institute of Molecular Biology and Genetics of the National Academy of Sciences of Ukraine. 37 prostate adenocarcinomas of different Gleason score (GS) and stages, 37 paired CNT and 20 samples of benign prostate tumors (A, adenomas) were studied. Tumors were characterized, according to the International System of Classiﬁcation of Tumors based on the tumor-node-metastasis (TNM) and the World Health Organization (WHO) criteria classiﬁcation. CPC of prostate cancer samples were described earlier .
Total RNA isolation and cDNA synthesis. 50–70 mg of frozen prostate tissues were mashed to powder in liquid nitrogen. Total RNA was extracted by TRI-reagent (Sigma, USA), according to the manufacturer’s protocol. The total RNA concentration was analyzed by a spectrophotometer (NanoDrop Technologies Inc., USA). The quality of the total RNA was determined in a 1% agarose gel by band intensity of 28S and 18S rRNA (28S/18S ratio). cDNA was synthesized from 1 µg of the total RNA treated with RNase free DNase I (Thermo Fisher Scientific, USA) using RevertAid H Minus M-MuLV Reverse Transcriptase (Thermo Fisher Scientific, USA) according to the manufacturer’s protocol.
Quantitative PCR. Levels of a RE of 23 genes were assessed by quantitative PCR using a Maxima SYBR Green Master mix (Thermo Fisher Scientific, USA) on Bio-Rad CFX96 Real-Time PCR Detection System (USA) under the following conditions: 95 °C — 10 min, following 40 cycles of 95 °C — 15 s, 60 °C — 30 s, elongation 72 °C — 30 s. Primers for all genes were selected from a qPrimerDepot (https://primerdepot.nci. nih.gov/) database and confirmed, using an https:// http://www.ncbi.nlm.nih.gov/tools/primer-blast/algorithm.
Four reference genes TBP, HPRT, ALAS1 and TUBA1B were used for gene expression normalization  in comparison with TBP . Two main models (2-ΔCt and 2-ΔΔCt methods), described earlier [25, 27] were used for calculation and an analysis of RE levels.
Statistical analysis. The Kolmogorov — Smirnov test was applied to assess the normality of distribution. The Wilcoxon Matched Pairs test was performed for comparison of RE in prostate adenocarcinomas and paired normal tissues samples. RE fold differences in 2-ΔΔCt model were considered increased in cases of changes more than 2 fold (RE fold changes > 2.01) and decreased when changes were more than 2 fold (RE fold changes < 0.49), when gene RE were compared in adenocarcinoma versus CNT. The Fisher exact test was calculated to find out differences between these sample groups . The Benjamini — Hochberg procedure with false discovery rate (FDR) 0.10–0.25 was used under multiple comparisons detection . The Kruskal — Wallis test was used to determine differences between experimental groups. The Dunn — Bonferroni post hoc test for multiple comparisons was performed to determine RE differences between pairs of prostate samples. The Spearman’s rank correlation test was used to find possible correlations between gene RE and CPC of prostate tumors .
We have analyzed RE of 23 genes in prostate T, CNT and A. Eight genes are associated with fibroblasts and CAFs, we called them a CAF group. Six genes are associated with pro-inflammatory and immuno-suppressive macrophages (the TAM group). Nine genes RE is associated with immune defense and reactions (the immune-associated genes (IAG) group).
From the beginning, we have found differences in RE in T, compared with CNT (N) samples. Statistically significant differences in RE in the paired T/N samples are shown in Table 1 (FDR = 0.2).
Table 1. Genes with significant RE differences between prostate adenocarcinoma samples and paired CNT (2-ΔΔCt model)
Differences in RE were calculated for 11 genes out of 23 in the paired T/N samples. Increased RE (> 2.01) in T in comparison with the paired N (p < 0.05) was observed for 7 genes: ACTA2, CXCL14, CTGF, THY1, FAP, CD163, CCL17. Decreased RE (< 0.49) in T in comparison with the paired N (p < 0.05) was observed for 4 genes: CCR4, NOS2A, MSMB, IL1R1.
Next, RE was analyzed in 3 sample groups, namely T, the paired CNT (N) and A. Furthermore, the GS, stage, presence of the TMPRSS2/ERG fusion , levels of PTEN expression (low (l) or high (h)) [25, 29] was taken into consideration for grouping of the samples. Summarized data are presented in Table 2.
Table 2. Significant RE changes according to Dunn — Bonferroni post hoc test for multiple comparisons in investigated genes
Note: #Presence of RE differences only in TBP normalized genes.
These data show that the greatest number of changes in RE, for 14 out of 23 genes could be found in TNA groups, when the stage of tumor was analyzed as well. The expression pattern of these genes and also another 2 genes (ACTA2 and MSMB), for which RE was changed only in the paired T/N samples is shown on Figure. It should be noted, that for two genes (FAP and CD68) the use of four reference genes and one reference gene for normalization of RE resulted in various fold change differences.
Fig. 1. RE profiles of genes with differences in different stages of disease in prostate adenocarcinomas (carcinomas), CNT and adenomas (a–e), f — genes with RE differences in paired T/N
Descriptive statistics of the data on RE in T, N and A samples at the various stages and also significant differences between groups are shown in Table 3. Among genes showing the significant differences in RE (14 genes) eight genes showed the lowest RE in the A group, namely CTGF, CXCL14, FAP, THY1, CD163, CCL17, CCL22, and CTLA4. Increased RE in T at stages 1–2 (p < 0.05) was observed for CTGF, FAP, THY1, CCL22, CTLA4, whereas increased RE in T stages 3–4 (p < 0.05) was found for CD163 and CCL17. CXCL14 was upregulated in all tumors.
Table 3. Descriptive statistics of the RE data in T, N, A sample groups at the various of disease stages and significant differences between pairs of groups
Table 3 (breakover)
Note: *Dunn — Bonferroni post hoc method for multiple comparisons with FDR = 0.2 (p < 0.05); #presence of RE differences only in TBP normalized genes; st — stage; 25.000th — 25th percentile, 75.000th — 75th percentile.
CXCL12, HIF1A, CD68, CCR4, NOS2A, and IL1R showed the highest RE in the A group 1. The significantly decreased RE (more than two-fold) in T stages 3–4 (p < 0.05) was demonstrated for five genes, namely HIF1A, CD68, CCR4, NOS2A and IL1R1. Only CXCL12 was downregulated in all tumors. Significant changes in RE (up- and downregulation) in the CNT group of patients with tumors of the various stages was also found for several of these genes (Table 4). For example, decreased RE of HIF1A and IL1R1 was detected in N compared with A in patients with T stages 3–4. The CXCL12 was downregulated in all normal tissues. Increased RE in the CNT group was found for four genes. The CXCL14, CD163 and CCL17 genes were upregulated in T stages 3–4. One gene (CTLA4) was upregulated in all tumors.
Table 4. RE differences between pairs in TNA groups with TMPRSS2/ERG fusion status (4A) and different PTEN RE (4B)
Note: *Dunn — Bonferroni post hoc method for multiple comparisons with FDR = 0.2, F+ and F– groups with presence (+) and absence (–) of TMPRSS2/ERG fusion, h/l — groups with high (h) and low (l) PTEN RE.
Nine genes in groups with different GS showed similar RE changes as genes in groups with different stages (Table 2). IL2RA was altered only in groups with different GS. Pairs with RE variability in the TNA group and depending on the presence/absence of the TMPRSS2/ERG fusion and also of PTEN RE are shown in Table 4. We have found 10 genes in TNA groups with the different TMPRSS2/ERG fusion status, showing significant differences in RE (p < 0.05). The most interesting are genes, which were differently expressed in T with and without the TMPRSS2/ERG fusion. Three genes — THY1, FAP and CCR4 — showed alterations only in T without the fusion, compared with the A group. Four genes (CTLA4, CCL17, CCL22) were expressed differently in the CNT group without the fusion, compared with the A group. One gene was (CD163) expressed at various levels in the fusion negative compared with the fusion-positive CNT groups.
Eight genes expressed at the differential levels in TNA groups (p < 0.05), when PTEN mRNA was analyzed. The most interesting among them are THY1 and FAP that were expressed at significantly different RE levels in adenocarcinomas with high PTEN (T h) and in adenomas with high PTEN RE (A h). CCR4 and IL1R1 showed differences in adenocarcinomas with low PTEN RE (T l) and adenomas with high PTEN RE (A h).
Correlation analysis of RE related to CPC and molecular characteristics of tumors. Using the Spearman’s rank correlation (rs) test, a number of correlations was suggested between the gene expression pattern and CPC, the TMPRSS2/ERG fusion status and PTEN levels in prostate adenocarcinomas (Table 5). RE of the nine genes, namely CXCL14, HIF1A, S100A4, CCL17, CCR4, NOS2A, CTLA4, IL2RA and MSMB correlated with GS. RE of the nine genes (HIF1A, S100A4, CD68, CD163, CCL17, CCL22, CCR4, NOS2A, IL1R1) showed significant correlations with the stage. Six genes out of nine belong to the TAM group. CD163 RE demonstrated the maximal positive correlation with the stage (rs = 0.615). RE of the THY1 gene correlated significantly with an age. RE of CXCL14, IL1RL1 and IRF1 correlated with PSA levels. RE of ACTA2 and CXCL12 negatively correlated with the TMPRSS2/ERG fusion status. These both genes belong to the CAF group. RE of other genes from the CAF group, namely ACTA2, CXCL12, FAP and S100A4 positively correlated with PTEN levels.
Table 5. Spearman rank order correlations (rs) of CPC and molecular characteristics with genes RE
Note: p < 0.05 (red); p < 0.01 (red bold italic); p < 0.001 (red bold).
We have to mention, that several genes show correlation in their expression levels. For example, in the T group, expression of CD163 positively correlated with S100A4 levels, as did FAP and THY1, CXCL12 and IL2RA, and also CTLA4 and IL2RA. The correlation indexes were high — rs = |0.633–0.712|, and correlation was p < 0.0001. The largest number of correlations out of 23 investigated genes (p < 0.01) demonstrated CD68, CCR4 and CCL22, i.e. each of these genes showed 8 correlations. For NOS2A, seven dependences were found, for HIF1A — six and for CTLA4, CIAS, IL2RA — five.
A surgical material of prostate cancer is very heterogeneous, as a rule . Tumor samples usually show a different structure with various sets of cells (i.e. cancerous cells, normal epithelial cells, stromal and immune cells, etc.) [31, 32]. Therefore, a study on specific tumor-host interactions is a very difficult task. To solve this problem, at least, partially, specific markers of different cell types were selected and gene RE was assessed in benign and malignant prostate tumors. We aimed to find the clinically significant alterations in the gene expression pattern. We found significant changes in RE only in a proportion of the investigated genes. Thus, 6 out of 8 genes of the CAF group showed differences in RE in the TNA group, containing tumors at various stages. Of note, the known CAF marker ACTA2  showed RE changes only in the group of the paired T/N samples.
The CXCL14, CTGF, FAP and THY1 also differentially expressed in tumors and CNT (N), and they showed alterations in RE in the TNA group, including tumors of the various stages. These four genes that are upregulated in tumors encode proteins secreted by CAFs. These proteins play an important role in paracrine interactions in the tumor microenvironment, enhancing tumor progression and metastasizing . Moreover, the FAP protein exhibits immune-suppressive properties . RE of HIF1A and S100A4 (also the CAF genes) correlated with the stage of tumors. Only CAF genes demonstrated correlation of RE with molecular characteristics of prostate tumors, namely the TMPRSS2/ERG fusion status and PTEN levels. Probably, this indicates the closest and selective type of interaction of the tumor and fibroblasts, as matrix-regulating elements of the prostate stroma [13, 15, 33].
Genes of the TAM group showed altered RE in the TNA group, including tumors at different stages. Some of genes demonstrated RE changes in the TNA group, when GS and the TMPRSS2/ERG fusion status varied. All of the 6 TAM genes correlated with tumor stage. We may speculate on the clinical significance of the selected markers of the TAM group [10, 16, 19]. Genes CD68 and NOS2A are the markers of normal, non-cancerous macrophages , whereas CD163 is the TAM marker . CCR4, CCL17 and CCL22 are associated with transition of the macrophage into TAM, i.e. these genes play a role in tumor progression .
The smallest number of genes with altered RE is found in the IAG group. Only two genes (CTLA4 and IL1R1) showed changes in RE in the TNA group, including the tumors of different stages. RE of IL1R1and MSMB was changed in the paired T/N group. It is known that inhibition of CTLA4 expression is very important in immune checkpoint therapy of prostate cancer, due to the fact that CTLA4 transmits the immune inhibitory signals to T cells . An IL1R1 signaling suppresses proliferation of mammary tumor cells, early in tumorigenesis . We found the significant decrease of IL1R1 RE in prostate cancer at the advanced stages of tumor development.
As mentioned earlier, the histological analysis indicated focal manifestations of leukocyte infiltration in the prostate tumors . Therefore, expression of leukocyte genes could vary quite much within even one sample [21, 32]. Based on this, it was important for us to assess the gene expression in the various immune cells. Also, the immune score indexes could have their own independent diagnostic and prognostic values . More suitable markers of prostate carcinogenesis should be studied, to reveal clinically significant alterations of immune response.
The observed high dispersion of the RE values did not make it possible to find the real difference between the groups of samples in all genes. In addition, the lack of the adequate control (a healthy donor’s prostate with the relevant age and nationality) also affected experimental results. Moreover, the changes in expression of the studied genes in adenomas and CNT might have a specific character, showing the effect of hypertrophy and hyperplasia in adenomas and the influence of malignant tumors on prostate CNT, especially at the late stages of tumor development [31, 32]. This assumption is confirmed by the results of the present study. We got the different number of genes with RE changes in the group of paired T/N (11 genes) and in adenomas, taking into account the tumor stage (14 genes). Of note, ACTA2 and MSMB showed RE changes only in the paired samples. Unexpectedly, ACTA2 was highly upregulated in adenomas. Probably, this is due to the specific expression of various ACTA2 isoforms upon tumor development and the specific expression of ACTA2 in different cell types, like smooth muscle cells [33, 34].
Earlier, we have shown that, based on RE of prostate-specific genes and EMT-associated genes, samples with the same clinical parameters could be grouped into specific RE subtypes [24, 25]. To confirm this, more experiments should be done on the larger patient cohort.
Several gene expression patterns were identified that correlated with the GS, stage and molecular characteristics of tumors, i.e. presence of the TMPRSS2/ERG fusion and alterations in PTEN expression. These expression patterns can be used for molecular profiling of prostate tumors, with the aim to develop personalized medicine approaches. However, the proposed profiling requires a more detailed analysis and a larger cohort of patients with prostate tumor.
1. Saraon P, Jarvi K, Diamandis EP. Molecular alterations during progression of prostate cancer to androgen independence. Clin Chem 2011; 57: 1366–75.
No Comments » Add comments