Tumor immune microenvironment characteristics of papillary thyroid carcinoma are associated with histopathological aggressiveness and BRAF mutation status

CXCL16 promotes tumor metastasis by regulating angiogenesis in the tumor micro-environment of BRAF V600E mutant colorectal cancer

Patients of colorectal cancer (CRC) with BRAF V600E mutation obtain poor prognosis. This study aimed to explore the role and mechanism of BRAF V600E mutation in angiogenesis of tumor micro-environment (TME). It has been reported that CXCL16 expression in TME is closely related to BRAF mutation. Clinicopathological features of CRC with BRAF V600E mutant or wild type were collected in this study. Immunohistochemistry (IHC) assays were conducted to test the expressions of vascular endothelial growth factor (VEGF), CD31 and CXCL16. ROC curve was used to determine the optimal cut off values of CXCL16. A total of 680 patients including 141 BRAF V600E type and 679 wild type were included. BRAF V600E mutant tumors were presented with significant worse clinicopathological features and a shorter overall survival (OS) than wild-type. Besides, chemokines CXCL16 was up-regulated in BRAF V600E mutant tissues and was associated with poorer prognosis. In addition, VEGF levels and vascular endothelial cell density was significantly increased in BRAF mutation. At last, CXCL16 was positively correlated with VEGF expression and vascular endothelial cell density. In conclusion, BRAF V600E mutations may promote metastasis of CRC by regulating CXCL16 expression and promoting angiogenesis in the TME.

Circulating immunophenotypes are potentially prognostic in follicular cell-derived thyroid cancer

Background

Exploring the immune interface of follicular cell-derived thyroid cancer has prognostic and therapeutic potential. The available literature is lacking for comprehensive immunophenotyping in relation to clinical outcomes. In this study, we identify circulating immunophenotypes associated with thyroid cancer prognosis.

Methods

We conducted a pilot observational study of adults with follicular cell-derived thyroid cancer who underwent surgery at our tertiary care referral center and had consented for flow cytometry on peripheral blood collected at the time of thyroidectomy.

Results

Of the 32 included subjects, 20 (62%) had well differentiated, 5 (16%) had poorly differentiated, and 7 (22%) had anaplastic thyroid cancer. The most frequent AJCC stage was 4 (59%) and the ATA risk of recurrence category was high (56%). Patients with AJCC stage 3/4 demonstrated fewer circulating mononuclear cells (CD45+), more monocytes (CD14+), fewer total lymphocytes (CD14-), fewer T cells (CD3+), fewer CD4+ T cells, fewer gamma-delta T cells, fewer natural killer (NK) T-like cells, more myeloid-derived suppressor cells (MDSCs; Lin-CD33+HLADR-), and more effector memory T cells but similar CD8+ T cells compared to stage1/2. Immunophenotype comparisons by ATA risk stratification and course of thyroid cancer were comparable to those observed for stage, except for significant differences in memory T cell subtypes. The median follow-up was 58 months.

Conclusions

Aggressive follicular cell-derived thyroid cancer either at presentation or during follow-up is associated with down-regulation of the T cell populations specifically CD4+ T cells, gamma-delta T cells, and NK T-like cells but up-regulation of MDSCs and altered memory T cells. These immunophenotypes are potential prognostic biomarkers supporting future investigation for developing targeted immunotherapies against advanced thyroid cancer.

Mining Prognostic Biomarkers of Thyroid Cancer Patients Based on the Immune-Related Genes and Development of a Reliable Prognostic Risk Model

PURPOSE

Tumor immunity serves an essential role in the occurrence and development of thyroid cancer (THCA). The aim of this study is to establish an immune-related prognostic model for THCA patients by using immune-related genes (IRGs).

METHODS

Wilcox test was used to screen the differentially expressed immune-related genes (DEIRGs) in THCA and normal tissues, then the DEIRGs related to prognosis were identified using univariate Cox regression analysis. According to The Cancer Genome Atlas (TCGA) cohort, we developed a least absolute shrinkage and selection operator (LASSO) regression prognostic model and performed validation analyses regard to the predictive value of the model in internal (TCGA) and external (International Cancer Genome Consortium) cohorts respectively. Finally, we analyzed the correlation among the prognostic model, clinical variables, and immune cell infiltration.

RESULTS

Eighty-two of 2,498 IRGs were differentially expressed between THCA and normal tissues, and 18 of them were related to prognosis. LASSO Cox regression analysis identified seven DEIRGs with the greatest prognostic value to construct the prognostic model. The risk model showed high predictive value for the survival of THCA in two independent cohorts. The risk score according to the risk model was positively associated with poor survival and the infiltration levels of immune cells, it can evaluate the prognosis of THCA patients independent of any other clinicopathologic feature. The prognostic value and genetic alternations of seven risk genes were evaluated separately.

CONCLUSION

Our study established and verified a dependable prognostic model associated with immune for THCA, both the identified IRGs and immune-related risk model were clinically significant, which is conducive to promoting individualized immunotherapy against THCA.

Receptor–ligand pair typing and prognostic risk model for papillary thyroid carcinoma based on single-cell sequencing

The papillary thyroid carcinoma (PTC) microenvironment consists of various cancer and surrounding cells, and the communication between them is mainly performed through ligand–receptor (LR) interactions. Single-cell RNA sequencing (scRNA-seq) has been performed to investigate the role of intercellular communication networks in tumor progression. In addition, scRNA-seq can accurately identify the characteristics of immune cell subsets, which is of great significance for predicting the efficacy of immunotherapy. In this study, the cell–cell communication network was analyzed through LR pairs, and a new PTC molecular phenotype was developed based on LR pairs. Furthermore, a risk model was established to predict patient response to PD-1 blockade immunotherapy. The scRNA-seq dataset was obtained from GSE184362, and the bulk tumor RNA-seq dataset was obtained from The Cancer Genome Atlas. CellPhoneDB was used for cellular communication analysis. LR pair correlations were calculated and used to identify molecular subtypes, and the least absolute shrinkage and selection operator (Lasso) Cox regression was used to develop a risk model based on LR pairs. The IMvigor210 and GSE78220 cohorts were used as external validations for the LR.score to predict responses to PD-L1 blockade therapy. A total of 149 LR pairs with significant expression and prognostic correlation were included, and three PTC molecular subtypes were obtained from those with significant prognostic differences. Then, five LR pairs were selected to construct the risk scoring model, a reliable and independent prognostic factor in the training set, test set, and whole dataset. Furthermore, two external validation sets confirmed the predictive efficacy of the LR.score for response to PD-1 blockade therapy.

Spatially resolved immune microenvironmental profiling for follicular thyroid carcinoma with minimal capsular invasion

Spatial profiles of the tumor-immune microenvironment are associated with disease progression and clinicopathological factors in various cancers. Follicular thyroid carcinoma (FTC) is the second most common thyroid cancer, where the presence of capsular invasion or angioinvasion determines the pathological diagnosis; however, little is known about the immune microenvironment profiles associated with the acquisition of invasive potential of FTC. In this study, we focused on FTC with minimal capsular invasion, and the spatially resolved immune microenvironment of FTC was studied in the discovery (n = 13) and validation cohorts (n = 40). CD8⁺ T cells, helper T cells, regulatory T cells, B cells, natural killer cells, tumor-associated macrophages, CD66⁺ granulocytes, mature dendritic cells, and mast cells were quantitatively evaluated in single tissue sections, via a 12-marker multiplex immunohistochemistry and image cytometry. Cell densities and compositions of immune cells were spatially stratified by six tissue regions including tumor center, subcapsular region, capsular invasion, adjacent stroma of capsular invasion, peritumoral stroma, and adjacent normal. Lymphoid cell lineages in the tumor center and subcapsular regions were significantly lower than those in adjacent normal and peritumoral stroma, potentially related to the lymphoid lineage exclusion from the intratumoral regions of FTC. Interestingly, immune cell composition profiles in the capsular invasive front were distinct from those of intratumoral region. The ratios of T cells to CD66b⁺ granulocytes with capsular invasion were significantly higher than those without capsular invasion, suggesting the presence of a unique immune microenvironment at the invasive front between tumor foci and stroma. In addition, tumor cells at the capsular invasive front showed significantly higher expression of tumor programmed cell death ligand 1 (PD-L1) than those at the tumor center. This study revealed spatial immune profiles associated with capsular invasion of FTC, providing new insights into the mechanisms underlying its development and initial invasion.

Characterization of the tumor immune microenvironment of sinonasal squamous-cell carcinoma

BACKGROUND

Treatment and prognosis of sinonasal squamous-cell carcinoma (SNSCC) have not significantly improved despite improvements in radical therapy. Characterization of the tumor immune microenvironment (TiME) may identify patient subgroups associated with disease recurrence, and provide new biomarkers for improved patient stratification and treatment.

METHODS

The TiME was quantitatively evaluated by multiplex immunohistochemistry (mIHC) in archived tissue sections from 38 patients with SNSCC, and were assessed for differences between recurrent (n = 20) and nonrecurrent (n = 18) groups. Hierarchical clustering analyses were performed to identify phenotypic TiME subgroups within the cohort and were used to compare survival outcomes.

RESULTS

Our mIHC analysis revealed increased T-cell populations and decreased myeloid-cell populations in SNSCC patients without recurrent disease, as compared with patients with recurrent disease. Within T-cell subsets, there was a significantly higher percentage of granzyme B⁺ , T-bet⁺ , Eomes⁺ T cells, as well as higher proliferation of CD8⁺ T cells within the nonrecurrent group relative to the recurrent group. Furthermore, immune-cell complexity profiles of SNSCC revealed hyper- and hypo-T-cell-inflamed, myeloid-inflamed, B-cell-inflamed, and broadly hypoinflamed subtypes not previously identified by gene expression analyses. Our study revealed that presence of either hyper- or hypo-T-cell-inflamed TiME subtypes were associated with increased survival outcomes as compared with broadly hypoinflamed TiME subtypes (p = 0.035 and 0.0376, respectively).

CONCLUSIONS

The TiME of SNSCC reveals distinct subtypes, which may correlate with recurrence and survival outcomes.

BRAFV600E mutant colorectal cancer cells mediate local immunosuppressive microenvironment through exosomal long noncoding RNAs

BACKGROUND

BRAF^V600E mutated colorectal cancer (CRC) is prone to peritoneal and distant lymph node metastasis and this correlates with a poor prognosis. The BRAF^V600E mutation is closely related to the formation of an immunosuppressive microenvironment. However, the correlation between BRAF^V600E mutation and changes in local immune microenvironment of CRC is not clear.

AIM

To explore the effect and mechanism of BRAF^V600E mutant on the immune microenvironment of CRC.

METHODS

Thirty patients with CRC were included in this study: 20 in a control group and 10 in a treatment group. The density of microvessels and microlymphatic vessels, and M2 subtype macrophages in tumor tissues were detected by immunohistochemistry. Screening and functional analysis of exosomal long noncoding RNAs (lncRNAs) were performed by transcriptomics. The proliferation and migration of human umbilical vein endothelial cells (HUVECs) and human lymphatic endothelial cells (HLECs) were detected by CCK-8 assay and scratch test, respectively. The tube-forming ability of endothelial cells was detected by tube formation assay. The macrophage subtypes were obtained by flow cytometry. The expression of vascular endothelial growth factor (VEGF)-A, basic fibroblast growth factor (bFGF), transforming growth factor (TGF)-β1, VEGF-C, claudin-5, occludin, zonula occludens (ZO)-1, fibroblast activation protein, and α-smooth muscle actin was assessed by western blot analysis. The levels of cytokines interleukin (IL)-6, TGF-β1, and VEGF were assessed by enzyme-linked immunosorbent assay.

RESULTS

BRAF^V600E mutation was positively correlated with the increase of preoperative serum carbohydrate antigen 19-9 (P < 0.05), and with poor tumor tissue differentiation in CRC (P < 0.01). Microvascular density and microlymphatic vessel density in BRAF^V600E mutant CRC tissues were higher than those in BRAF wild-type CRC (P < 0.05). The number of CD163⁺ M2 macrophages in BRAF^V600E mutant CRC tumor tissue was markedly increased (P < 0.05). Compared with exosomes from CRC cells with BRAF gene silencing, the expression of 13 lncRNAs and 192 mRNAs in the exosomes from BRAF^V600E mutant CRC cells was upregulated, and the expression of 22 lncRNAs and 236 mRNAs was downregulated (P < 0.05). The biological functions and signaling pathways predicted by differential lncRNA target genes and differential mRNAs were closely related to angiogenesis, tumor cell proliferation, differentiation, metabolism, and changes in the microenvironment. The proliferation, migration, and tube formation ability of HUVECs and HLECs induced by exosomes in the 1627 cell group (HT29 cells with BRAF gene silencing) was greatly reduced compared with the HT29 cell group (P < 0.05). Compared with the HT29 cell group, the expression levels of VEGF-A, bFGF, TGF-β1, and VEGF-C in the exosomes derived from 1627 cells were reduced. The expression of ZO-1 in HUVECs, and claudin-5, occludin, and ZO-1 in HLECs of the 1627 cell group was higher. Compared with the 1627 cell group, the exosomes of the HT29 cell group promoted the expression of CD163 in macrophages (P < 0.05). IL-6 secretion by macrophages in the HT29 cell group was markedly elevated (P < 0.05), whereas TGF-β1 was decreased (P < 0.05). The levels of IL-6, TGF-β1, and VEGF secreted by fibroblasts in the 1627 cell group decreased, compared with the HT29 cell group (P < 0.05).

CONCLUSION

BRAF^V600E mutant CRC cells can reach the tumor microenvironment by releasing exosomal lncRNAs, and induce the formation of an immunosuppressive microenvironment.

m5C Regulator-Mediated Methylation Modification Patterns and Tumor Microenvironment Infiltration Characterization in Papillary Thyroid Carcinoma

Recently, immune response modulation at the epigenetic level is illustrated in studies, but the possible function of RNA 5-methylcytosine (m⁵C) modification in cell infiltration within the tumor microenvironment (TME) is still unclear. Three different m⁵C modification patterns were identified, and high differentiation degree was observed in the cell infiltration features within TME under the above three identified patterns. A low m⁵C-score, which was reflected in the activated immunity, predicted the relatively favorable prognostic outcome. A small amount of effective immune infiltration was seen in the high m⁵C-score subtype, indicating the dismal patient survival. Our study constructed a diagnostic model using the 10 signature genes highly related to the m⁵C-score, discovered that the model exhibited high diagnostic accuracy for PTC, and screened out five potential drugs for PTC based on this m5C-score model. m⁵C modification exerts an important part in forming the TME complexity and diversity. It is valuable to evaluate the m⁵C modification patterns in single tumors, so as to enhance our understanding towards the infiltration characterization in TME.

Genomic Instability-Related LncRNA Signature Predicts the Prognosis and Highlights LINC01614 Is a Tumor Microenvironment-Related Oncogenic lncRNA of Papillary Thyroid Carcinoma

Background

Genomic instability (GI) is among the top ten characteristics of malignancy. Long non-coding RNAs (lncRNAs) are promising cancer biomarkers that are reportedly involved in GI. So far, the clinical value of GI-related lncRNAs (GIlncs) in papillary thyroid cancer (PTC) has not been clarified.

Methods

Integrative analysis of lncRNA expression and somatic mutation profiles was performed to identify GIlncs. Analysis of differentially expressed lncRNAs in the group with high- and low- cumulative number of somatic mutations revealed significant GIlncs in PTC. Univariate and multivariate Cox proportional hazard regression analyses were performed to identify hub-GIlncs.

Results

A computational model based on four lncRNAs (FOXD2-AS1, LINC01614, AC073257.2, and AC005082.1) was identified as a quantitative index using an in-silicon discovery cohort. GILS score was significantly associated with poor prognosis, as validated in the TCGA dataset and further tested in our local RNA-Seq cohort. Moreover, a combination of clinical characteristics and the composite GILS-clinical prognostic nomogram demonstrates satisfactory discrimination and calibration. Furthermore, the GILS score and FOXD2-AS1, LINC01614, AC073257.2, and AC005082.1 were also associated with driver mutations and multiple clinical-pathological variables, respectively. Moreover, RNA-Seq confirmed the expression patterns of FOXD2-AS1, LINC01614, AC073257.2, and AC005082.1 in PTC and normal thyroid tissues. Biological experiments demonstrated that downregulated or overexpressed LINC01614 affect PTC cell proliferation, migration, and invasion in vitro. Activation of the stromal and immune cell infiltration was also observed in the high LINC01614 group in the PTC microenvironment.

Conclusion

In summary, we identified a signature for clinical outcome prediction in PTC comprising four lncRNAs associated with GI. A better understanding of the GI providing an alternative evaluation of the progression risk of PTC. Our study also demonstrated LINC01614 as a novel oncogenic lncRNA and verified its phenotype in PTC.

A Novel Immune-Related Prognostic Signature for Thyroid Carcinoma

Background:

Recent studies have confirmed that immune-associated genes perform a crucial function in recurrence and metastasis of thyroid carcinoma. A reliable immune-related prognostic signature for patients with thyroid cancer is needed. This study constructed a novel immune-related prognostic signature for thyroid cancer and evaluated its prognostic value by bioinformatics analysis.

Methods:

In this study, we anatomized differentially expressed immune-associated genes from The Cancer Genome Atlas database. The samples from The Cancer Genome Atlas database were randomly divided into training set and test set. A novel immune-related prognostic signature for thyroid cancer was developed by least absolute shrinkage and selection operator and Cox regression analysis: Risk score = (0.6846 × expression value of C-X-C motif chemokine ligand 5 [CXCL5]) + (1.1556 × expression value of Azurocidin 1 [AZU1]) + (−0.3156 × expression value of nucleotide binding oligomerization domain containing 1 [NOD1] + (0.0542 × expression value of TNF Receptor Superfamily Member 11b [TNFRSF11B]) + (0.0952 × expression value of VGF nerve growth factor inducible [VGF]). The established prognostic signature was evaluated based on training set and test set by survival curves, receiver–operator characteristic curves, risk score, survival status, heatmap, and independent prognostic analysis. Meanwhile, we appraised the correlation between target immune-associated genes and clinical stage, tumor-infiltrating immune cells respectively.

Results:

Five immune-associated genes were used for constructing an immune-related prognostic signature by least absolute shrinkage and selection operator, univariate, and multivariate analysis. Survival curves, receiver–operator characteristic curves, and independent prognostic analysis showed the signature had significant prediction value. Clinical and immune cell correlation analyses indicated that target immune-associated genes may participate in tumor immune infiltration and tumor progression.

Conclusions:

We constructed a novel 5 immune-associated genes signature for predicting the prognosis of patients with thyroid cancer, which may help clinical workers evaluate individualized therapy and prognosis.

Prognostic significance of spatial immune profiles in human solid cancers

Immune-based tumor characteristics in the context of tumor heterogeneity are associated with suppression as well as promotion of cancer progression in various tumor types. As immunity typically functions based on intercellular contacts and short-distance cytokine communications, the location and spatial relationships of the tumor immune microenvironment can provide a framework to understand the biology and potential predictive biomarkers related to disease outcomes. Immune spatial analysis is a newly emerging form of cancer research based on recent methodological advances in in situ single-cell analysis, where cell-cell interaction and the tissue architecture can be analyzed in relation to phenotyping the tumor immune heterogeneity. Spatial characteristics of tumors can be stratified into the tissue architecture level and the single-cell level. At the tissue architecture level, the prognostic significance of the density of immune cell lineages, particularly T cells, is leveraged by understanding longitudinal changes in cell distribution in the tissue architecture such as intra-tumoral and peri-tumoral regions, and invasive margins. At the single-cell level, the proximity of the tumor to the immune cells correlates with disease aggressiveness and therapeutic resistance, providing evidence to understand biological interactions and characteristics of the tumor immune microenvironment. In this review, we summarize recent findings regarding spatial information of the tumor immune microenvironment and review advances and challenges in spatial single-cell analysis toward developing tissue-based biomarkers rooted in the immune spatial landscape.

Mechanisms of the Impact of Hashimoto Thyroiditis on Papillary Thyroid Carcinoma Progression: Relationship with the Tumor Immune Microenvironment

BACKGROUND

The relationship between Hashimoto thyroiditis (HT) and papillary thyroid carcinoma (PTC) remains uncertain. We assessed the impact of HT on the tumor immune microenvironment (TIME) in PTC.

METHODS

Thirty patients with PTC (group 1) and 30 patients with PTC and HT (group 2) were enrolled in this pilot study. The distribution and number of CD8+ lymphocytes, plasma cells (CD138+), regulatory T cells (forkhead box P3 [FOXP3+)], mast cell tryptase (MCT+), and M2 macrophages (CD163+) were evaluated. To test the hypothesis that HT impacts PTC development via signal transducer and activator of transcription 6 (STAT6) activation and M2 macrophage polarization, we investigated STAT6 expression in tumor and stromal cells. We also evaluated vascular endothelial growth factor (VEGF) expression by lymph node metastasis (LNM) status.

RESULTS

TIME showed significant between-group differences. Group 1 patients demonstrated immune desert or immune-excluded immunophenotypes, while an inflamed phenotype with more CD8+ cells (P&lt;0.001) predominated in group 2. Immune-excluded TIME was associated with the highest LNM rate. In PTC, LNM was associated with more numerous CD163+ cells. Moreover, LNM in group 1 was associated with increased numbers of mast cells peritumorally and FOXP3+ cells intratumorally and peritumorally. Group 2 demonstrated higher STAT6 but not higher VEGF expression in tumor cells. High VEGF expression was associated with LNM regardless of HT status.

CONCLUSION

Concomitant HT impacted PTC signaling via STAT6 and TIME by increasing the number of CD8+ cells. LNM is associated with increases in CD163+ cells and VEGF expression in PTC, whereas HT affected LNM through different mechanisms.

Hashimoto's thyroiditis attenuates progression of papillary thyroid carcinoma: deciphering immunological links

Although some studies have investigated the clinicopathologic relationships between papillary thyroid carcinoma (PTC) and Hashimoto's thyroiditis (HT), there is still no clear understanding of differences in tumor immune microenvironment for PTC with coexisting HT and HT effect on PTC progression. The aim of this study was to clarify immune-mediated mechanisms of coexisting HT, which might influence PTC progression. 30 patients with histologically confirmed conventional-type PTC and 30 patients with PTC and coexisting HT were enrolled in the study. To analyze the role of immune-mediated links between PTC and HT, immunohistochemical investigation was conducted to count the number of different immune cells including T-cytotoxic cells (CD8), plasma cells (CD138), Treg cells (FOXP3), mast cells (MCT), and M2 macrophages (CD163). It was shown that despite the high number of immune cells in the intact thyroid tissues of PTC patients with coexisting HT there were no significant differences in M2 macrophages, mast cells and Treg counts inside PTC with or without HT. PTC with HT was associated with a higher number of CD8+ cells (P < 0.001) reflecting the ability of immune system to generate and recruit T-cytotoxic cells in tumor area, which can explain the protective effect of HT on PTC progression. Lymph node metastases development was associated with an increased number of mast cells, M2 macrophages and Treg along with a decreased plasma cells count regardless of coexisting HT. However, we did not find significant differences in T-cytotoxic cells quantity in node-positive and node-negative patients with or without HT, which encourages further investigation of immune escape mechanisms in PTC.

BCL2 and hsa-miR-181a-5p are potential biomarkers associated with papillary thyroid cancer based on bioinformatics analysis

Background

The morbidity of thyroid carcinoma has been rising worldwide and increasing faster than any other cancer type. The most common subtype with the best prognosis is papillary thyroid cancer (PTC); however, the exact molecular pathogenesis of PTC is still not completely understood.

Methods

In the current study, 3 gene expression datasets (GSE3678, GSE3467, and GSE33630) and 2 miRNA expression datasets (GSE113629 and GSE73182) of PTC were selected from the Gene Expression Omnibus (GEO) database and were further used to identify differentially expressed genes (DEGs) and deregulated miRNAs between normal thyroid tissue samples and PTC samples. Then, Gene Ontology (GO) and pathway enrichment analyses were conducted, and a protein-protein interaction (PPI) network was constructed to explore the potential mechanism of PTC carcinogenesis. The hub gene detection was performed using the CentiScaPe v2.0 plugin, and significant modules were discovered using the MCODE plugin for Cytoscape. In addition, a miRNA-gene regulatory network in PTC was constructed using common deregulated miRNAs and DEGs.

Results

A total of 263 common DEGs and 12 common deregulated miRNAs were identified. Then, 6 significant KEGG pathways (P < 0.05) and 82 significant GO terms were found to be enriched, indicating that PTC was closely related to amino acid metabolism, development, immune system, and endocrine system. In addition, by constructing a PPI network and miRNA-gene regulatory network, we found that hsa-miR-181a-5p regulated the most DEGs, while BCL2 was targeted by the most miRNAs.

Conclusions

The results of this study suggested that hsa-miR-181a-5p and BCL2 and their regulatory networks may play important roles in the pathogenesis of PTC.

PSMA expression level predicts differentiated thyroid cancer aggressiveness and patient outcome

Background

Prostate-specific membrane antigen (PSMA) is overexpressed on the endothelial cells of tumor neo-vessels of several solid malignancies, including differentiated thyroid cancer (DTC). We aimed to test the potential role of PSMA as a biomarker for DTC aggressiveness and outcome prediction.

We retrospectively screened all patients who underwent thyroidectomy between 1 January 2010 and 31 December 2017 in our institution. Applying the inclusion (histological diagnosis of thyroid cancer and tissue availability) and exclusion criteria (no clinical or follow-up data or diagnosis of medullary thyroid cancer), a cohort of 59 patients was selected. The monoclonal mouse anti-human PSMA antibody was used to stain tissue sections. A 3-point scale was used to score PSMA positivity: 0–5% expression was considered as negative (score 0), 6–50% as moderately positive (score 1), and 51–100% as highly positive (score 2). A cumulative score (0–10%, 11–79%, and 80–100%) was also explored. Univariate and multivariate logistic regression analyses were performed to predict the presence of distant metastases, chosen as endpoint of aggressiveness. The area under the curve (AUC) was calculated. Cox models were built to predict patient outcome in terms of recurrence, iodine refractoriness, and status at last follow-up, which were calculated using the Kaplan-Meier failure function.

Results

At immunostaining, 12, 25, and 22 patients had scores of 0, 1, and 2, respectively. According to the cumulative score, PSMA expression was ≤ 10% in 17 cases, 11–79% in 31 cases, and ≥ 80% in 11 cases. At multivariate analysis, age, sex, histotype, vascular invasion, T and N parameters, and PSMA positivity were significant predictors of distant metastases. The AUC was 0.92. Recurrence or progression occurred in 19/59 patients. Twelve patients developed radioiodine (RAI) refractoriness, after a median time of 17 months (range 2–32). One patient died of DTC; 46 of the 58 patients alive at last follow-up were disease free. Median DFS was 23 months (range 3–82). The final multivariate model to predict RAI refractoriness included as covariates the stage, high PSMA expression (≥ 80%), and the interaction between moderate PSMA expression (11–79%) and stage.

Conclusions

PSMA, a marker of neovasculature formation expressed by DTC, contributes in the prediction of tumor aggressiveness and patient outcome.

Electronic supplementary material

The online version of this article (10.1186/s13550-019-0559-9) contains supplementary material, which is available to authorized users.

High-dimensional multiplexed immunohistochemical characterization of immune contexture in human cancers

Biomarker assessments of tumor specimens is widely used in cancer research to audit tumor cell intrinsic as well as tumor cell extrinsic features including the diversity of immune, stromal, and mesenchymal cells. To comprehensively and quantitatively audit the tumor-immune microenvironment (TiME), we developed a novel multiplex immunohistochemistry (mIHC) platform and computational image processing workflow using a single formalin-fixed paraffin-embedded (FFPE) tissue section. Herein, we validated this platform using nine matched primary newly diagnosed and recurrent head and neck squamous cell carcinoma (HNSCC) sections sequentially subjected to immunodetection with a panel of 29 antibodies identifying malignant tumor cells, and 17 distinct leukocyte lineages and their functional states. Image cytometric analysis was applied to interpret chromogenic signals from digitally scanned and coregistered light microscopy-based images enabling identification and quantification of individual tumor cells, structural features, immune cell phenotypes and their functional state. In agreement with our previous study via a 12-plex imaging mIHC platform, myeloid-inflamed status in newly diagnosed primary tumors associated with significantly short progression free survival, independent of lymphoid-inflamed status. Spatial distribution of tumor and immune cell lineages in TiME was also examined and revealed statistically significant CD8⁺ T cell exclusion from tumor nests, whereas regulatory T cells and myeloid cells, when present in close proximity to tumor cells, highly associated with rapid cancer recurrence. These findings indicate presence of differential immune-spatial profiles in newly diagnosed and recurrent HNSCC, and establish the robustness of the 29-plex mIHC platform and associated analytics for quantitative analysis of single tissue sections revealing longitudinal TiME changes.