Network dynamics in the tumor microenvironment

Measuring the composition of the tumor microenvironment with transcriptome analysis: past, present and future

Interactions between tumor cells and immune cells in the tumor microenvironment (TME) play a vital role the mechanisms of immune evasion, by which cancer cells escape immune elimination. Thus, the characterization and quantification of different components in the TME is a hot topic in molecular biology and drug discovery. Since the development of transcriptome sequencing in bulk tissue, single cells and spatial dimensions, there are increasing methods emerging to deconvolute and subtype the TME. This review discusses and compares such computational strategies and downstream subtyping analyses. Integrative analyses of the transcriptome with other data, such as epigenetics and T-cell receptor sequencing, are needed to obtain comprehensive knowledge of the dynamic TME.

Tumor-associated macrophages confer colorectal cancer 5-fluorouracil resistance by promoting MRP1 membrane translocation via an intercellular CXCL17/CXCL22-CCR4-ATF6-GRP78 axis

Chemotherapy represents a major type of clinical treatment against colorectal cancer (CRC). Aberrant drug efflux mediated by transporters acts as a key approach for tumor cells to acquire chemotherapy resistance. Increasing evidence implies that tumor-associated macrophages (TAMs) play a pivotal role in both tumorigenesis and drug resistance. Nevertheless, the specific mechanism through which TAMs regulate drug efflux remains elusive. Here, we discovered that TAMs endow CRC cells with resistance to 5-fluorouracil (5-FU) treatment via a cell-cell interaction-mediated MRP1-dependent drug efflux process. Mechanistically, TAM-secreted C-C motif chemokine ligand 17 (CCL17) and CCL22, via membrane receptor CCR4, activated the PI3K/AKT pathway in CRC tumor cells. Specifically, phosphorylation of AKT inactivated IP3R and induced calcium aggregation in the ER, resulting in the activation of ATF6 and upregulation of GRP78. Accordingly, excessive GRP78 can interact with MRP1 and promote its translocation to the cell membrane, causing TAM-induced 5-FU efflux. Taken together, our results demonstrated that TAMs promote CRC chemotherapy resistance via elevating the expression of GRP78 to promote the membrane translocation of MRP1 and drug efflux, providing direct proof for TAM-induced drug resistance.

Evaluating the relevance of surgical margins. Part one: The problems with current methodology

The goal of cancer surgery is to achieve a "clean" microscopic resection, with no residual tumour remaining in the wound. To achieve that goal, the surgeon typically incorporates a measured buffer of grossly normal tissue about the entire circumference of the tumour. Microscopic analysis of the resection boundaries is then performed to determine if all traces of the tumour have been completely removed. This analysis is thought to provide a surrogate indication as to the likelihood for that tumour to recur after surgery. However, it is recognised that tumour recurrence may not occur even when microscopic evidence of tumour has been identified at the resection margins, and recurrence can also occur when conventional histology has considered the tumour to have been completely removed. The explanations for this dichotomy are numerous and include technical and practical limitations of the processing methodology, and also several surgeon-related and tumour-related reasons. Ultimately, the inability to confidently determine when a tumour has been removed sufficiently to prevent recurrence can impact on the ability to provide owners with confident treatment advice. In this article, the authors describe the challenges with defining the true extent of the tumour margin from the perspective of the surgeon, the pathologist and the tumour. The authors also provide an analysis of why our current efforts to ensure that all traces of the local tumour have been successfully removed may provide an imperfect assessment of the risk of recurrence.

Metabolic rearrangements and intratumoral heterogeneity for immune response in hepatocellular carcinoma

Liver cancer is one of the most common malignant tumors globally. Not only is it difficult to diagnose, but treatments are scarce and the prognosis is generally poor. Hepatocellular carcinoma (HCC) is the most common type of liver cancer. Aggressive cancer cells, such as those found in HCC, undergo extensive metabolic rewiring as tumorigenesis, the unique feature, ultimately causes adaptation to the neoplastic microenvironment. Intratumoral heterogeneity (ITH) is defined as the presence of distinct genetic features and different phenotypes in the same tumoral region. ITH, a property unique to malignant cancers, results in differences in many different features of tumors, including, but not limited to, tumor growth and resistance to chemotherapy, which in turn is partly responsible for metabolic reprogramming. Moreover, the different metabolic phenotypes might also activate the immune response to varying degrees and help tumor cells escape detection by the immune system. In this review, we summarize the reprogramming of glucose metabolism and tumoral heterogeneity and their associations that occur in HCC, to obtain a better understanding of the mechanisms of HCC oncogenesis.

Characterization of pyruvate metabolism and citric acid cycle patterns predicts response to immunotherapeutic and ferroptosis in gastric cancer

Background

Gastric cancer is one of the most common malignancies of the digestive system with a high lethal rate. Studies have shown that inherited and acquired mutations in pyruvate metabolism and citric acid cycle (P-CA) enzymes are involved in tumorigenesis and tumor development. However, it is unclear how different P-CA patterns affect the tumor microenvironment (TME), which is critical for cancer progression.

Methods

This study mainly concentrated on investigating the role of the P-CA patterns in multicellular immune cell infiltration of GC TME. First, the expression levels of P-CA regulators were profiled in GC samples from The Cancer Genome Atlas and Gene Expression Omnibus cohorts to construct a consensus clustering analysis and identify three distinct P-CA clusters. GSVA was conducted to reveal the different biological processes in three P-CA clusters. Subsequently, 1127 cluster-related differentially expressed genes were identified, and prognostic-related genes were screened using univariate Cox regression analysis. A scoring system was then set up to quantify the P-CA gene signature and further evaluate the response of the patients to the immunotherapy.

Results

We found that GC patients in the high P-CA score group had a higher tumor mutational burden, higher microsatellite instability, and better prognosis. The opposite was observed in the low P-CA score group. Interestingly, we demonstrated P-CA gene cluster could predict the sensitivity to immunotherapy and ferroptosis-induced therapy.

Conclusion

Collectively, the P-CA gene signature in this study exhibits potential roles in the tumor microenvironment and predicts the response to immunotherapeutic. The identification of these P-CA patterns may significantly accelerate the strategic development of immunotherapy for GC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-022-02739-z.

The involvement of homeobox-C 4 in predicting prognosis and unraveling immune landscape across multiple cancers via integrated analysis

Background: There has been growing evidence that the aberrantly expressed Homeobox-C 4 (HOXC4) plays crucial roles in the development of some cancer types. However, it remains unclear as far as its expression patterns and prognostic significance are concerned, as is tumor immunity.

Methods: To investigate the expression levels and prognostic implications of HOXC4, multiple data sources were used in conjunction with quantitative real-time polymerase chain reaction (qRT-PCR) verification. Afterward, diverse immunological-related analyses, along with anti-cancer drug sensitivity, were performed in a number of cancer types. A further exploration of the underlying mechanisms of HOXC4 in tumorigenesis and immunity was carried out using the Gene Set Enrichment Analysis (GSEA) and the Gene Set Variation Analysis (GSVA).

Results: Based on extensive database mining, HOXC4 was ubiquitously expressed across 21 tumor cell lines and significantly higher than that of normal tissues in 21 tumor types. The outcome of survival analysis including overall survival (OS), disease-free interval (DFI), disease-specific survival (DSS) and progression-free interval (PFI) revealed that upregulation of HOXC4 expression in several cancers was associated with worse prognosis. Additionally, HOXC4 was observed to correlate closely with colon adenocarcinoma (COAD), head and neck squamous cell carcinoma (HNSC), lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), rectum adenocarcinoma (READ), and thyroid carcinoma (THCA) in terms of tumor immune cells infiltration. As a result of our comprehensive pan-cancer study, we have identified a significant link between the expression of HOXC4 and the efficacy of immunotherapy-related treatments, together with anti-cancer drug sensitivity. As a final note, HOXC4 was found to modulate multiple signaling pathways involved in tumorigenesis and immunity.

Conclusion: HOXC4 has been implicated in our study for the first time as an oncogene in cancers with a poor prognosis, potentially laying the groundwork for promising clinical biomarkers and immunotherapy approaches.

Recent developments in application of single-cell RNA sequencing in the tumour immune microenvironment and cancer therapy

The advent of single-cell RNA sequencing (scRNA-seq) has provided insight into the tumour immune microenvironment (TIME). This review focuses on the application of scRNA-seq in investigation of the TIME. Over time, scRNA-seq methods have evolved, and components of the TIME have been deciphered with high resolution. In this review, we first introduced the principle of scRNA-seq and compared different sequencing approaches. Novel cell types in the TIME, a continuous transitional state, and mutual intercommunication among TIME components present potential targets for prognosis prediction and treatment in cancer. Thus, we concluded novel cell clusters of cancer-associated fibroblasts (CAFs), T cells, tumour-associated macrophages (TAMs) and dendritic cells (DCs) discovered after the application of scRNA-seq in TIME. We also proposed the development of TAMs and exhausted T cells, as well as the possible targets to interrupt the process. In addition, the therapeutic interventions based on cellular interactions in TIME were also summarized. For decades, quantification of the TIME components has been adopted in clinical practice to predict patient survival and response to therapy and is expected to play an important role in the precise treatment of cancer. Summarizing the current findings, we believe that advances in technology and wide application of single-cell analysis can lead to the discovery of novel perspectives on cancer therapy, which can subsequently be implemented in the clinic. Finally, we propose some future directions in the field of TIME studies that can be aided by scRNA-seq technology.

Gene mining of immune microenvironment in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is a common malignant tumor worldwide with a poor prognosis. Recent studies have shown that the occurrence, development and prognosis of liver cancer are closely related to tumor microenvironment (TME) and tumor immune infiltration.

METHODS

Therefore, important information on various diseases can be obtained from public databases such as The Cancer Gene Atlas (TCGA), and ideas or schemes that may be effective for the treatment of various diseases can be screened and analyzed by screening various conditions. In this study, 424 cases of liver hepatocellular carcinoma (LIHC) in the TCGA database and CIBERSORT algorithm were used to calculate the proportion of tumor-invasive immune cells. Combined with the clinical data from TCGA database, it was concluded that T cells regulatory (Tregs) were correlated with the development and prognosis of HCC. Cox regression analysis was used to screen differentially expressed genes, and survival analysis was performed according to the screened differentially expressed genes to see whether there was a significant association with the prognosis of HCC. Then gene ontology and kyoto encyclopedia of genes and genomes analysis of differentially expressed genes were carried out to explore the possibility of differentially expressed genes becoming potential therapeutic targets of HCC.

RESULTS

Finally, I identified the gene centromere protein o (CENPO), which is associated with immune cells and improve the prognosis of HCC.

CONCLUSION

CENPO may be a potential biological therapeutic target for hepatocellular treatment.

Simulations of tumor growth and response to immunotherapy by coupling a spatial agent-based model with a whole-patient quantitative systems pharmacology model

Quantitative systems pharmacology (QSP) models and spatial agent-based models (ABM) are powerful and efficient approaches for the analysis of biological systems and for clinical applications. Although QSP models are becoming essential in discovering predictive biomarkers and developing combination therapies through in silico virtual trials, they are inadequate to capture the spatial heterogeneity and randomness that characterize complex biological systems, and specifically the tumor microenvironment. Here, we extend our recently developed spatial QSP (spQSP) model to analyze tumor growth dynamics and its response to immunotherapy at different spatio-temporal scales. In the model, the tumor spatial dynamics is governed by the ABM, coupled to the QSP model, which includes the following compartments: central (blood system), tumor, tumor-draining lymph node, and peripheral (the rest of the organs and tissues). A dynamic recruitment of T cells and myeloid-derived suppressor cells (MDSC) from the QSP central compartment has been implemented as a function of the spatial distribution of cancer cells. The proposed QSP-ABM coupling methodology enables the spQSP model to perform as a coarse-grained model at the whole-tumor scale and as an agent-based model at the regions of interest (ROIs) scale. Thus, we exploit the spQSP model potential to characterize tumor growth, identify T cell hotspots, and perform qualitative and quantitative descriptions of cell density profiles at the invasive front of the tumor. Additionally, we analyze the effects of immunotherapy at both whole-tumor and ROI scales under different tumor growth and immune response conditions. A digital pathology computational analysis of triple-negative breast cancer specimens is used as a guide for modeling the immuno-architecture of the invasive front.

Identifying Potential Biomarkers of Prognostic Value in Colorectal Cancer via Tumor Microenvironment Data Mining

Colorectal cancer (CRC) is a common cancer that has increased rapidly worldwide in the past decades with a relatively high mortality rate. An increasing body of evidence has highlighted the importance of infiltrating immune and stromal cells in CRC. In this study, based on gene expression data of CRC patients in TCGA database we evaluated immune and stromal scores in tumor microenvironment using ESTIMATE method. Results showed there was potential correlation between these scores and the prognosis, and that patients with higher immune score and lower stromal score had longer survival time. We found that immune score was correlated with clinical characteristics including tumor location, tumor stage, and survival time. Specifically, the right-sided colon cancer had markedly elevated immune score, compared to left-sided colon cancer and rectal cancer. These results might be useful for understanding tumor microenvironment in colorectal cancer. Through the differential analysis we got a list of genes significantly associated with immune and stromal scores. Gene Set Enrichment and protein-protein interaction network analysis were used to further illustrate these differentially expressed genes. Finally, 15 hub genes were identified, and three (CXCL9, CXCL10 and SELL) of them were validated with favorable outcomes in CRC patients. Our result suggested that these tumor microenvironment related genes might be potential biomarkers for the prognosis of CRC.

Disruption of Cell-Cell Communication in Anaplastic Thyroid Cancer as an Immunotherapeutic Opportunity

Thyroid cancer incidence is increasing at an alarming rate, almost tripling every decade. About 44,280 new cases of thyroid cancer (12,150 in men and 32,130 in women) are estimated to be diagnosed in 2021, with an estimated death toll of around 2200. Although most thyroid tumors are treatable and associated with a favorable outcome, anaplastic thyroid cancer (ATC) is extremely aggressive with a grim prognosis of 6-9 months post-diagnosis. A large contributing factor to this aggressive nature is that ATC is completely refractory to mainstream therapies. Analysis of the tumor microenvironment (TME) associated with ATC can relay insight to the pathological realm that encompasses tumors and aids in cancer progression and proliferation. The TME is defined as a complex niche that surrounds a tumor and involves a plethora of cellular components whose secretions can modulate the environment in order to favor tumor progression. The cellular heterogeneity of the TME contributes to its dynamic function due to the presence of both immune and nonimmune resident, infiltrating, and interacting cell types. Associated immune cells discussed in this chapter include macrophages, dendritic cells (DCs), natural killer (NK) cells, and tumor-infiltrating lymphocytes (TILs). Nonimmune cells also play a role in the establishment and proliferation of the TME, including neuroendocrine (NE) cells, adipocytes, endothelial cells (ECs), mesenchymal stem cells (MSCs), and fibroblasts. The dynamic nature of the TME contributes greatly to cancer progression.Recent work has found ATC tissues to be defined by a T cell-inflamed "hot" tumor immune microenvironment (TIME) as evidenced by presence of CD3+ and CD8+ T cells. These tumor types are amenable to immune checkpoint blockade (ICB) therapy. This therapeutic avenue, as of 2021, has remained unexplored in ATC. New studies should seek to explore the therapeutic feasibility of a combination therapy, through the use of a small molecule inhibitor with ICB in ATC. Screening of in vitro model systems representative of papillary, anaplastic, and follicular thyroid cancer explored the expression of 29 immune checkpoint molecules. There are higher expressions of HVEM, BTLA, and CD160 in ATC cell lines when compared to the other TC subtypes. The expression level of HVEM was more than 30-fold higher in ATC compared to the others, on average. HVEM is a member of tumor necrosis factor (TNF) receptor superfamily, which acts as a bidirectional switch through interaction with BTLA, CD160, and LIGHT, in a cis or trans manner. Given the T cell-inflamed hot TIME in ATC, expression of HVEM on tumor cells was suggestive of a possibility for complex crosstalk of HVEM with inflammatory cytokines. Altogether, there is emerging evidence of a T cell-inflamed TIME in ATC along with the expression of immune checkpoint proteins HVEM, BTLA, and CD160 in ATC. This can open doors for combination therapies using small molecule inhibitors targeting downstream effectors of MAPK pathway and antagonistic antibodies targeting the HVEM/BTLA axis as a potentially viable therapeutic avenue for ATC patients. With this being stated, the development of adaptive resistance to targeted therapies is inevitable; therefore, using a combination therapy that targets the TIME can serve as a preemptive tactic against the characteristic therapeutic resistance that is seen in ATC. The dynamic nature of the TME, including the immune cells, nonimmune cells, and acellular components, can serve as viable targets for combination therapy in ATC. Understanding the complex interactions of these associated cells and the paradigm in which their secretions and components can serve as immunomodulators are critical points of understanding when trying to develop therapeutics specifically tailored for the anaplastic thyroid carcinoma microenvironment.

Bioinformatics prediction and experimental verification identify MAD2L1 and CCNB2 as diagnostic biomarkers of rhabdomyosarcoma

Background

Rhabdomyosarcoma (RMS) is a malignant soft-tissue tumour. In recent years, the tumour microenvironment (TME) has been reported to be associated with the development of tumours. However, the relationship between the occurrence and development of RMS and TME is unclear. The purpose of this study is to identify potential tumor microenvironment-related biomarkers in rhabdomyosarcoma and analyze their molecular mechanisms, diagnostic and prognostic significance.

Methods

We first applied bioinformatics method to analyse the tumour samples of 125 patients with rhabdomyosarcoma (RMS) from the Gene Expression Omnibus database (GEO). Differential genes (DEGs) that significantly correlate with TME and the clinical staging of tumors were extracted. Immunohistochemistry (IHC) was applied to validate the expression of mitotic arrest deficient 2 like 1 (MAD2L1) and cyclin B2 (CCNB2) in RMS tissue. Then, we used cell function and molecular biology techniques to study the influence of MAD2L1 and CCNB2 expression levels on the progression of RMS.

Results

Bioinformatics results show that the RMS TME key genes were screened, and a TME-related tumour clinical staging model was constructed. The top 10 hub genes were screened through the establishment of a protein–protein interaction (PPI) network, and then Gene Expression Profiling Interactive Analysis (GEPIA) was conducted to measure the overall survival (OS) of the 10 hub genes in the sarcoma cases in The Cancer Genome Atlas (TCGA). Six DEGs of statistical significance were acquired. The relationship between these six differential genes and the clinical stage of RMS was analysed. Further analysis revealed that the OS of RMS patients with high expression of MAD2L1 and CCNB2 was worse and the expression of MAD2L1 and CCNB2 was related to the clinical stage of RMS patients. Gene set enrichment analysis (GSEA) revealed that the genes in MAD2L1 and CCNB2 groups with high expression were mainly related to the mechanism of tumour metastasis and recurrence. In the low-expression MAD2L1 and CCNB2 groups, the genes were enriched in the metabolic and immune pathways. Immunohistochemical results also confirmed that the expression levels of MAD2L1 (30/33, 87.5%) and CCNB2 (33/33, 100%) were remarkably higher in RMS group than in normal control group (0/11, 0%). Moreover, the expression of CCNB2 was related to tumour size. Downregulation of MAD2L1 and CCNB2 suppressed the growth, invasion, migration, and cell cycling of RMS cells and promoted their apoptosis. The CIBERSORT immune cell fraction analysis indicated that the expression levels of MAD2L1 and CCNB2 affected the immune status in the TME.

Conclusions

The expression levels of MAD2L1 and CCNB2 are potential indicators of TME status changes in RMS, which may help guide the prognosis of patients with RMS and the clinical staging of tumours.

Tumor-associated macrophages (TAMs) depend on MMP1 for their cancer-promoting role

The complex interaction between tumor-associated macrophages (TAMs) and tumor cells through several soluble factors and signaling is essential for colorectal cancer (CRC) progression. However, the molecular mechanism involved remains elusive. In this study, we demonstrated that MMP1 derived from TAMs markedly facilitated colon cancer cell proliferation via accelerating cell cycle transition from G0/G1 to S and G2/M phase. Moreover, exogenous MMP1 activated cdc25a/CDK4-cyclin D1 and p21/cdc2-cyclin B1 complexes through altering c-Myc and ETV4. Mechanistic studies indicated that inhibition of PAR1 or blockage of MAPK/Erk signaling eliminated the proliferation induced by exogenous MMP1 in vitro and in vivo. In addition, ETV4 could bind to the promoter of MMP1 and activate MMP1 transcription, which confirmed the MMP1/ETV4/MMP1 positive feedback. Altogether, our study identified a cytokine paracrine manner between colon cancer cells and TAMs. MMP1/PAR1/Erk1/2/ETV4 positive feedback loop may represent to be a therapeutic target and prognostic marker in CRC.

Pan-cancer characterization of lncRNA modifiers of immune microenvironment reveals clinically distinct de novo tumor subtypes

The emerging field of long noncoding RNA (lncRNA)-immunity has provided a new perspective on cancer immunity and immunotherapies. The lncRNA modifiers of infiltrating immune cells in the tumor immune microenvironment (TIME) and their impact on tumor behavior and disease prognosis remain largely uncharacterized. In the present study, a systems immunology framework integrating the noncoding transcriptome and immunogenomics profiles of 9549 tumor samples across 30 solid cancer types was used, and 36 lncRNAs were identified as modifier candidates underlying immune cell infiltration in the TIME at the pan-cancer level. These TIME lncRNA modifiers (TIL-lncRNAs) were able to subclassify various tumors into three de novo pan-cancer subtypes characterized by distinct immunological features, biological behaviors, and disease prognoses. Finally, a TIL-lncRNA-derived immune state index (TISI) that was reflective of immunological and oncogenic states but also predictive of patients' prognosis was proposed. Furthermore, the TISI provided additional prognostic value for existing tumor immunological and molecular subtypes. By applying the TISI to tumors from different clinical immunotherapy cohorts, the TISI was found to be significantly negatively correlated with immune-checkpoint genes and to have the ability to predict the effectiveness of immunotherapy. In conclusion, the present study provided comprehensive resources and insights for future functional and mechanistic studies on lncRNA-mediated cancer immunity and highlighted the potential of the clinical application of lncRNA-based immunotherapeutic strategies in precision immunotherapy.

Prognostic value of CCR2 as an immune indicator in lung adenocarcinoma: A study based on tumor-infiltrating immune cell analysis

Background

Prognostic indicators in lung adenocarcinoma (LUAD) have been seeking under database analysis, and remarkable advance is on the way.

Methods

This study calculated the scores of stromal and immune components of the tumor microenvironment (TME) in 551 LUAD samples using the ESTIMATE algorithm on The Cancer Genome Atlas (TCGA) database. R package ''limma'' was used to selected differentially expressed genes (DEG). We have analyzed the DEGs by means of Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichments. The protein‐protein network, univariate Cox analysis, and Lasso regression analysis were performed to selected survival‐related genes. Gene Set Enrichment Analysis (GSEA) represented the enriched pathway of CC chemokine receptor 2 (CCR2). The ratios of immune cells in the TME of each LUAD sample were obtained using the R package "limma" and CIBERSORT algorithm in R 4.0.2.

Results

The ImmuneScore was positively correlated with prognosis regarding survival rate, T classification of TNM stages, and clinicopathological staging characteristics. GO and KEGG enrichments showed DEGs were associated with immune‐related activities. Three genes of LUAD were selected from the PPI network and Cox proportional hazards regression analysis. CCR2 was the most survival correlated gene by Lasso regression analysis. GSEA results showed that C2 kegg gene sets in the CCR2 high‐expression group were mainly enriched in the B cell or T cell receptor signaling pathway and natural killer cell‐mediated cytotoxicity. Correlation of CCR2 expression with prognosis was conducted, implicating a positive correlation with the prognosis of survival rate and M classification, negative correlation with the prognosis of T and N classifications. The correlation between CCR2 and tumor‐infiltrating immune cells (TICs) was analyzed, and 14 kinds of TICs were found closely correlated with CCR2 expression through difference analysis.

Conclusion

Therefore, CCR2 has prognostic value as an immune indicator in LUAD.

LINC02257, an Enhancer RNA of Prognostic Value in Colon Adenocarcinoma, Correlates With Multi-Omics Immunotherapy-Related Analysis in 33 Cancers

Accumulated evidence supports that long non-coding RNAs (lncRNAs) are involved significantly in the development of human cancers. Enhancer RNAs (eRNAs), a subtype of lncRNAs, have recently attracted much attention about their roles in carcinogenesis. Colon adenocarcinoma is one of the most commonly diagnosed tumors with unfavorable prognosis. It highlights the great significance of screening and identifying novel biomarkers. More importantly, it remains to be elucidated with respect to the function of eRNAs in colon adenocarcinoma, as is in pan-cancers. The expression of LINC02257 was determined based on the data obtained from The Cancer Genome Atlas (TCGA). Further evaluation was performed on the basis of the following analyses: clinicopathology and survival analysis, gene ontology (GO) terms, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, as well as multi-omics immunotherapy-related analysis and co-expression analysis. The statistical analysis was conducted in R software, and immune cell infiltration of LINC02257 expression in cancers was investigated by using the CIBERSORT algorithm. By large-scale data mining, our study highlighted that a total of 39 eRNA genes were associated with colon adenocarcinoma prognosis, among which 25 eRNAs showed significant associations with their predicted target genes. LINC02257 was identified as the most significant survival-associated eRNA, with DUSP10 as its target gene. Besides, the high expression of LINC02257 in colon adenocarcinoma was more vulnerable to unfavorable prognosis and correlated with various clinical characteristics. GO and KEGG analyses revealed that LINC02257 was closely correlated with extracellular matrix organization via the PI3K-Akt signaling pathway. Besides, LINC02257 expression correlated with a multi-omics analysis of 33 cancer types, such as survival analysis [overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI)] and immunotherapy-related analysis [tumor microenvironment (TME), tumor mutational burden (TMB), and microsatellite instability (MSI)]. Finally, we investigated the co-expression genes of LINC02257 and its potential signaling pathways across different cancer types. LINC02257 is screened and can function as an independent prognostic biomarker through the PI3K-Akt signaling pathway for colon adenocarcinoma. Simultaneously, LINC02257 may be a multifaceted and significant immunotherapy-related eRNA in different cancers.

A model study for the manufacture and validation of clinical-grade deciduous dental pulp stem cells for chronic liver fibrosis treatment

BACKGROUND

Human deciduous pulp stem cells (hDPSCs) have remarkable stem cell potency associated with cell proliferation, mesenchymal multipotency, and immunosuppressive function and have shown beneficial effects in a variety of animal disease models. Recent studies demonstrated that hDPSCs exhibited in vivo anti-fibrotic and anti-inflammatory action and in vivo hepatogenic-associated liver regeneration, suggesting that hDPSCs may offer a promising source with great clinical demand for treating liver diseases. However, how to manufacture ex vivo large-scale clinical-grade hDPSCs with the appropriate quality, safety, and preclinical efficacy assurances remains unclear.

METHODS

We isolated hDPSCs from human deciduous dental pulp tissues formed by the colony-forming unit-fibroblast (CFU-F) method and expanded them under a xenogeneic-free and serum-free (XF/SF) condition; hDPSC products were subsequently stored by two-step banking including a master cell bank (MCB) and a working cell bank (WCB). The final products were directly thawed hDPSCs from the WCB. We tested the safety and quality check, stem cell properties, and preclinical potentials of final hDPSC products and hDPSC products in the MCB and WCB.

RESULTS

We optimized manufacturing procedures to isolate and expand hDPSC products under a XF/SF culture condition and established the MCB and the WCB. The final hDPSC products and hDPSC products in the MCB and WCB were validated the safety and quality including population doubling ability, chromosome stability, microorganism safety, and stem cell properties including morphology, cell surface marker expression, and multipotency. We also evaluated the in vivo immunogenicity and tumorigenicity and validated in vivo therapeutic efficacy for liver regeneration in a CCl₄-induced chronic liver fibrosis mouse model in the final hDPSC products and hDPSC products in the WCB.

CONCLUSION

The manufacture and quality control results indicated that the present procedure could produce sufficient numbers of clinical-grade hDPSC products from a tiny deciduous dental pulp tissue to enhance clinical application of hDPSC products in chronic liver fibrosis.

Mining TCGA Database for Tumor Microenvironment-Related Genes of Prognostic Value in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most common and lethal malignancies. Recent studies reveal that tumor microenvironment (TME) components significantly affect HCC growth and progression, particularly the infiltrating stromal and immune cells. Thus, mining of TME-related biomarkers is crucial to improve the survival of patients with HCC. Public access of The Cancer Genome Atlas (TCGA) database allows convenient performance of gene expression-based analysis of big data, which contributes to the exploration of potential association between genes and prognosis of a variety of malignancies, including HCC. The "Estimation of STromal and Immune cells in MAlignant Tumors using Expression data" algorithm renders the quantification of the stromal and immune components in TME possible by calculating the stromal and immune scores. Differentially expressed genes (DEGs) were screened by dividing the HCC cohort of TCGA database into high- and low-score groups according to stromal and immune scores. Further analyses of functional enrichment and protein-protein interaction networks show that the DEGs are mainly involved in immune response, cell adhesion, and extracellular matrix. Finally, seven DEGs have significant association with HCC poor outcomes. These genes contain FABP3, GALNT5, GPR84, ITGB6, MYEOV, PLEKHS1, and STRA6 and may be candidate biomarkers for HCC prognosis.

Tumor Microenvironment and Cell Fusion

Cell fusion is a highly regulated biological process that occurs under both physiological and pathological conditions. The cellular and extracellular environment is critical for the induction of the cell-cell fusion. Aberrant cell fusion is initiated during tumor progression. Tumor microenvironment is a complex dynamic system formed by the interaction between tumor cells and their surrounding cells. Cell-cell fusion mediates direct interaction between tumor cells and their surrounding cells and is associated with tumor initiation and progression. Various microenvironmental factors affect cell fusion in tumor microenvironment and generate hybrids that acquire genomes of both parental cells and exhibit novel characteristics, such as tumor stem cell-like properties, radioresistance, drug resistance, immune evasion, and enhanced migration and invasion abilities, which are closely related to the initiation, invasion, and metastasis of tumor. The phenotypic characteristics of hybrids are based on the phenotypes of parental cells, and the fusion of tumor cells with diverse types of microenvironmental fusogenic cells is concomitant with phenotypic heterogeneity. This review highlights the types of fusogenic cells in tumor microenvironment that can fuse with tumor cells and their specific significance and summarizes the various microenvironmental factors affecting tumor cell fusion. This review may be used as a reference to develop strategies for future research on tumor cell fusion and the exploration of cell fusion-based antitumor therapies.

Breast tumors educate the proteome of stromal tissue in an individualized but coordinated manner

Cancer forms specialized microenvironmental niches that promote local invasion and colonization. Engrafted patient-derived xenografts (PDXs) locally invade and colonize naïve stroma in mice while enabling unambiguous molecular discrimination of human proteins in the tumor from mouse proteins in the microenvironment. To characterize how patient breast tumors form a niche and educate naïve stroma, subcutaneous breast cancer PDXs were globally profiled by species-specific quantitative proteomics. Regulation of PDX stromal proteins by breast tumors was extensive, with 35% of the stromal proteome altered by tumors consistently across different animals and passages. Differentially regulated proteins in the stroma clustered into six signatures, which included both known and previously unappreciated contributors to tumor invasion and colonization. Stromal proteomes were coordinately regulated; however, the sets of proteins altered by each tumor were highly distinct. Integrated analysis of tumor and stromal proteins, a comparison made possible in these xenograft models, indicated that the known hallmarks of cancer contribute pleiotropically to establishing and maintaining the microenvironmental niche of the tumor. Education of the stroma by the tumor is therefore an intrinsic property of breast tumors that is highly individualized, yet proceeds by consistent, nonrandom, and defined tumor-promoting molecular alterations.

A novel individual-cell-based mathematical model based on multicellular tumour spheroids for evaluating doxorubicin-related delivery in avascular regions

BACKGROUND AND PURPOSE

Effective drug delivery in the avascular regions of tumours, which is crucial for the promising antitumour activity of doxorubicin-related therapy, is governed by two inseparable processes: intercellular diffusion and intracellular retention. To accurately evaluate doxorubicin-related delivery in the avascular regions, these two processes should be assessed together. Here we describe a new approach to such an assessment.

EXPERIMENTAL APPROACH

An individual-cell-based mathematical model based on multicellular tumour spheroids was developed that describes the different intercellular diffusion and intracellular retention kinetics of doxorubicin in each cell layer. The different effects of a P-glycoprotein inhibitor (LY335979) and a hypoxia inhibitor (YC-1) were quantitatively evaluated and compared, in vitro (tumour spheroids) and in vivo (HepG2 tumours in mice). This approach was further tested by evaluating in these models, an experimental doxorubicin derivative, INNO 206, which is in Phase II clinical trials.

KEY RESULTS

Inhomogeneous, hypoxia-induced, P-glycoprotein expression compromised active transport of doxorubicin in the central area, that is, far from the vasculature. LY335979 inhibited efflux due to P-glycoprotein but limited levels of doxorubicin outside the inner cells, whereas YC-1 co-administration specifically increased doxorubicin accumulation in the inner cells without affecting the extracellular levels. INNO 206 exhibited a more effective distribution profile than doxorubicin.

CONCLUSIONS AND IMPLICATIONS

The individual-cell-based mathematical model accurately evaluated and predicted doxorubicin-related delivery and regulation in the avascular regions of tumours. The described framework provides a mechanistic basis for the proper development of doxorubicin-related drug co-administration profiles and nanoparticle development and could avoid unnecessary clinical trials.

Inhibition of Lysyl Oxidases Impairs Migration and Angiogenic Properties of Tumor-Associated Pericytes

Pericytes are important cellular components of the tumor microenviroment with established roles in angiogenesis and metastasis. These two cancer hallmarks are modulated by enzymes of the LOX family, but thus far, information about LOX relevance in tumor-associated pericytes is lacking. Here, we performed a comparative characterization of normal and tumoral pericytes and report for the first time the modulatory effects of LOX enzymes on activated pericyte properties. Tumoral pericytes isolated from childhood ependymoma and neuroblastoma specimens displayed angiogenic properties in vitro and expressed typical markers, including CD146, NG2, and PDGFRβ. Expression of all LOX family members could be detected in both normal and tumor-associated pericytes. In most pericyte samples, LOXL3 was the family member displaying the highest transcript levels. Inhibition of LOX/LOXL activity with the inhibitor β-aminopropionitrile (βAPN) significantly reduced migration of pericytes, while proliferation rates were kept unaltered. Formation of tube-like structures in vitro by pericytes was also significantly impaired upon inhibition of LOX/LOXL activity with βAPN, which induced more prominent effects in tumor-associated pericytes. These findings reveal a novel involvement of the LOX family of enzymes in migration and angiogenic properties of pericytes, with implications in tumor development and in therapeutic targeting tumor microenvironment constituents.

Weakly coupled map lattice models for multicellular patterning and collective normalization of abnormal single-cell states

We present a weakly coupled map lattice model for patterning that explores the effects exerted by weakening the local dynamic rules on model biological and artificial networks composed of two-state building blocks (cells). To this end, we use two cellular automata models based on (i) a smooth majority rule (model I) and (ii) a set of rules similar to those of Conway's Game of Life (model II). The normal and abnormal cell states evolve according to local rules that are modulated by a parameter κ. This parameter quantifies the effective weakening of the prescribed rules due to the limited coupling of each cell to its neighborhood and can be experimentally controlled by appropriate external agents. The emergent spatiotemporal maps of single-cell states should be of significance for positional information processes as well as for intercellular communication in tumorigenesis, where the collective normalization of abnormal single-cell states by a predominantly normal neighborhood may be crucial.

Whole tumor RNA-sequencing and deconvolution reveal a clinically-prognostic PTEN/PI3K-regulated glioma transcriptional signature

The concept that solid tumors are maintained by a productive interplay between neoplastic and non-neoplastic elements has gained traction with the demonstration that stromal fibroblasts and immune system cells dictate cancer development and progression. While less studied, brain tumor (glioma) biology is likewise influenced by non-neoplastic immune system cells (macrophages and microglia) which interact with neoplastic glioma cells to create a unique physiological state (glioma ecosystem) distinct from that found in the normal tissue. To explore this neoplastic ground state, we leveraged several preclinical mouse models of neurofibromatosis type 1 (NF1) optic glioma, a low-grade astrocytoma whose formation and maintenance requires productive interactions between non-neoplastic and neoplastic cells, and employed whole tumor RNA-sequencing and mathematical deconvolution strategies to characterize this low-grade glioma ecosystem as an aggregate of cellular and acellular elements. Using this approach, we demonstrate that optic gliomas generated by altering the germline Nf1 gene mutation, the glioma cell of origin, or the presence of co-existing genetic alterations represent molecularly-distinct tumors. However, these optic glioma tumors share a 25-gene core signature, not found in normal optic nerve, that is normalized by microglia inhibition (minocycline), but not conventional (carboplatin) or molecularly-targeted (rapamycin) chemotherapy. Lastly, we identify a genetic signature conferred by Pten reduction and corrected by PI3K inhibition. This signature predicts progression-free survival in patients with either low-grade or high-grade glioma. Collectively, these findings support the concept that gliomas are composite ecological systems whose biology and response to therapy may be best defined by examining the tumor as a whole.

Tumorigenicity analysis of heterogeneous dental stem cells and its self-modification for chromosome instability

Heterogeneity demonstrates that stem cells are constituted by several sub-clones in various differentiation states. The heterogeneous state is maintained by cross-talk among sub-clones, thereby ensuring stem cell adaption. In this study, we investigated the roles of heterogeneity on genetic stability. Three sub-clones (DF2, DF8 and DF18) were isolated from heterogeneous dental stem cells (DSCs), and were proved to be chromosome instability (CIN) after long term expansion. Cell apoptosis were not detected in sub-clones, which exhibited strong tumorigenesis tendency, coupled with weak expression of p53 and aberrant ultra-structure. However, 3 sub-clones did not overexpress tumor related markers or induce tumorigenesis in vivo. The mixed-culture study suggested that 3-clone-mixed culturing cells (DF1) presented apparent decrease in the ratio of aneuploidy. The screening experiment further proved that 3 sub-clones functioned separately in this modification procedure but only mixed culturing all 3 sub-clones, simulated heterogeneous microenvironment, could achieve complete modification. Additionally, osteogenesis capability of 3 sub-clones was partially influenced by CIN while DSCs still kept stronger osteogenesis than sub-clones. These results suggested aberrant sub-clones isolated from heterogeneous DSCs were not tumorigenesis and could modify CIN by cross-talk among themselves, indicating that the heterogeneity played a key role in maintaining genetic stability and differentiation capability in dental stem cells.

The Significance of an Enhanced Concept of the Organism for Medicine

Recent developments in evolutionary biology, comparative embryology, and systems biology suggest the necessity of a conceptual shift in the way we think about organisms. It is becoming increasingly evident that molecular and genetic processes are subject to extremely refined regulation and control by the cell and the organism, so that it becomes hard to define single molecular functions or certain genes as primary causes of specific processes. Rather, the molecular level is integrated into highly regulated networks within the respective systems. This has consequences for medical research in general, especially for the basic concept of personalized medicine or precision medicine. Here an integrative systems concept is proposed that describes the organism as a multilevel, highly flexible, adaptable, and, in this sense, autonomous basis for a human individual. The hypothesis is developed that these properties of the organism, gained from scientific observation, will gradually make it necessary to rethink the conceptual framework of physiology and pathophysiology in medicine.

Macrophages play an essential role in the long effects of low-dose photodynamic therapy on vessel permeability

Low-dose photodynamic therapy (L-PDT) has been used to transiently increase the permeability of tumor vessels to improve the delivery of chemotherapeutic drugs to lung tumors. However, the effects of L-PDT in a long-term on delivery of chemotherapeutic drugs are unknown. In this study, we studied this question as well as the underlying mechanisms. We found that the effects of L-PDT on tumor vessel permeability appeared to be prolonged. Moreover, L-PDT significantly increased the number of tumor associated macrophages, and appeared to induce macrophage polarization to a M1 phenotype. Further analyses showed that L-PDT upregulated stromal cell-derived factor 1 (SDF-1) in tumor to recruit macrophages through a SDF-1/Chemokine (C-X-C Motif) Receptor 4 (CXCR4) axis, which accounted for the prolonged effects of L-PDT on vessel permeability. Application of a specific CXCR4 inhibitor significantly suppressed the L-PDT-induced macrophage recruitment, resulting in abolishment of the prolonged effects of L-PDT on vessel permeability. Furthermore, the inhibitory effects of Liporubicin™ on the growth of the implanted tumor in L-PDT-treated mice were significantly attenuated by CXCR4 inhibition. Thus, our data demonstrate a previously unappreciated long-lasting effect of L-PDT on vessel permeability, and suggest that this long-lasting effects of L-PDT treatment on vessel permeability may result from modulation of macrophage recruitment and polarization. Hence, L-PDT may be a promising method to assist chemotherapeutic approaches.

Effect and Molecular Mechanisms of Traditional Chinese Medicine on Regulating Tumor Immunosuppressive Microenvironment

Traditional Chinese medicine (TCM) is an important complementary strategy for treating cancer in China. The mechanism is related to regulating the internal environment and remodeling the tumor immunosuppressive microenvironment (TIM). Herein we illustrate how TIM is reformed and its protumor activity on promoting tumor cell proliferation, angiogenesis and lymphangiogenesis, tumor invasion, and the oncogenicity of cancer stem cells. Furthermore we summarize the effects and mechanism of TCM on regulating TIM via enhancing antitumor immune responses (e.g., regulating the expression of MHC molecules and Fas/FasL, attenuating cancerigenic ability of cancer stem cells) and remolding immunosuppressive cells (e.g., reversing immune phenotypes of T lymphocytes and tumor associated macrophages, promoting dendritic cells mature, restraining myeloid derived suppressor cells function, and regulating Th1/Th2 factors). We also reveal the bidirectional and multitargeting functions of TCM on regulating TIM. Hopefully, it provides new theoretical basis for TCM clinical practice in cancer treatment and prevention.

Electrical coupling in ensembles of nonexcitable cells: modeling the spatial map of single cell potentials

We analyze the coupling of model nonexcitable (non-neural) cells assuming that the cell membrane potential is the basic individual property. We obtain this potential on the basis of the inward and outward rectifying voltage-gated channels characteristic of cell membranes. We concentrate on the electrical coupling of a cell ensemble rather than on the biochemical and mechanical characteristics of the individual cells, obtain the map of single cell potentials using simple assumptions, and suggest procedures to collectively modify this spatial map. The response of the cell ensemble to an external perturbation and the consequences of cell isolation, heterogeneity, and ensemble size are also analyzed. The results suggest that simple coupling mechanisms can be significant for the biophysical chemistry of model biomolecular ensembles. In particular, the spatiotemporal map of single cell potentials should be relevant for the uptake and distribution of charged nanoparticles over model cell ensembles and the collective properties of droplet networks incorporating protein ion channels inserted in lipid bilayers.