Engineering Multidimensional Evolutionary Forces to Combat Cancer

Pan-cancer analysis and single-cell analysis reveals FAM110B as a potential target for survival and immunotherapy

Background: FAM110B belongs to the family that has a 110 sequence similarity (FAM110) and is located in the centrosome and mitotic spindle. FAM110B has been linked to tumor cell growth in earlier research. Uncertainty exists regarding FAM110B's function within the tumor microenvironment is unclear as well as pan-cancer. Methods: In order to assess the variation in FAM110B expression within normal and pan-cancer tissues, we combined the TCGA and GTEx databases. The cBioPortal database and the GSCALite platform were used to examine the variation in genome and methylation alteration of FAM110B. Cox regression, Kaplan-Meier, and SangerBox were employed to examine the clinical features and prognosis of FAM110B and pan-cancer. The purpose of the correlational research was to investigate the associations within immunerelated genes, tumor mutation burden, microsatellite instability, immune-related genes, and immunological checkpoints and FAM110B expression. ESTIMATE, EPIC, QUANTISEQ, and MCPCOUNTER methods were used to calculate the interaction among FAM110B expression as well as the tumor immune microenvironment. The immunoinfiltration and function of FAM110B were analyzed by single-cell databases (TISCH and CancerSEA). Finally, we evaluated the sensitivity of FAM110B to small-molecule medications through GDSC and CTRP databases. Results: The transcription and protein expression of FAM110B varies significantly throughout cancer types, and this has predictive value for the prognosis of some tumors; including brain lower grade glioma (LGG), stomach adenocarcinoma (STAD), pancreatic adenocarcinoma (PAAD), etc. In the tumor microenvironment, the expression level of FAM110B was associated with immune cell infiltration, immune checkpoint immune regulatory genes, tumor mutational burden, and microsatellite fragility to a certain extent. Conclusion: This work investigates the possibility of utility of FAM110B as a marker to forecast pan-cancer immunotherapy response, providing a theoretical basis for cancer therapy.

Peripheral Blood CD8+ T-Lymphocyte Immune Response in Benign and Subpopulations of Breast Cancer Patients

Peripheral blood CD8+ T lymphocytes play a crucial role in cell-mediated immunity and tumor-related immune responses in breast cancer. In this study, label-free quantification analysis and gene set enrichment analysis (GSEA) of CD8+ T lymphocytes in the peripheral blood of benign patients and patients with different breast cancer (BC) subtypes, i.e., luminal A, luminal B, and triple-negative breast cancer (TNBC), were performed using nano-UHPLC and Orbitrap mass spectrometry. Differential protein expression in CD8+ T lymphocytes revealed significant downregulation (log2 FC ≥ 0.38 or ≤-0.38, adj. p < 0.05), particularly in proteins involved in cytotoxicity, cytolysis, and proteolysis, such as granzymes (GZMs) and perforin 1 (PRF1). This downregulation was observed in the benign group (GZMH, GZMM, and PRF1) and luminal B (GZMA, GZMH) subtypes, whereas granzyme K (GZMK) was upregulated in TNBC in comparison to healthy controls. The RNA degradation pathway was significantly downregulated (p < 0.05, normalized enrichment score (NES) from -1.47 to -1.80) across all BC subtypes, suggesting a potential mechanism for regulating gene expression during T cell activation. Also, the Sm-like proteins (LSM2, LSM3, and LSM5) were significantly downregulated in the RNA degradation pathway. Proteomic analysis of CD8+ T lymphocytes in peripheral blood across different breast cancer subtypes provides a comprehensive view of the molecular mechanisms of the systemic immune response that can significantly contribute to advancements in the diagnosis, treatment, and prognosis of this disease.

Choline metabolism reprogramming mediates an immunosuppressive microenvironment in non-small cell lung cancer (NSCLC) by promoting tumor-associated macrophage functional polarization and endothelial cell proliferation

INTRODUCTION

Lung cancer is a prevalent malignancy globally, and immunotherapy has revolutionized its treatment. However, resistance to immunotherapy remains a challenge. Abnormal cholinesterase (ChE) activity and choline metabolism are associated with tumor oncogenesis, progression, and poor prognosis in multiple cancers. Yet, the precise mechanism underlying the relationship between ChE, choline metabolism and tumor immune microenvironment in lung cancer, and the response and resistance of immunotherapy still unclear.

METHODS

Firstly, 277 advanced non-small cell lung cancer (NSCLC) patients receiving first-line immunotherapy in Sun Yat-sen University Cancer Center were enrolled in the study. Pretreatment and the alteration of ChE after 2 courses of immunotherapy and survival outcomes were collected. Kaplan-Meier survival and cox regression analysis were performed, and nomogram was conducted to identify the prognostic and predicted values. Secondly, choline metabolism-related genes were screened using Cox regression, and a prognostic model was constructed. Functional enrichment analysis and immune microenvironment analysis were also conducted. Lastly, to gain further insights into potential mechanisms, single-cell analysis was performed.

RESULTS

Firstly, baseline high level ChE and the elevation of ChE after immunotherapy were significantly associated with better survival outcomes for advanced NSCLC. Constructed nomogram based on the significant variables from the multivariate Cox analysis performed well in discrimination and calibration. Secondly, 4 choline metabolism-related genes (MTHFD1, PDGFB, PIK3R3, CHKB) were screened and developed a risk signature that was found to be related to a poorer prognosis. Further analysis revealed that the choline metabolism-related genes signature was associated with immunosuppressive tumor microenvironment, immune escape and metabolic reprogramming. scRNA-seq showed that MTHFD1 was specifically distributed in tumor-associated macrophages (TAMs), mediating the differentiation and immunosuppressive functions of macrophages, which may potentially impact endothelial cell proliferation and tumor angiogenesis.

CONCLUSION

Our study highlights the discovery of ChE as a prognostic marker in advanced NSCLC, suggesting its potential for identifying patients who may benefit from immunotherapy. Additionally, we developed a prognostic signature based on choline metabolism-related genes, revealing the correlation with the immunosuppressive microenvironment and uncovering the role of MTHFD1 in macrophage differentiation and endothelial cell proliferation, providing insights into the intricate workings of choline metabolism in NSCLC pathogenesis.

FAM family gene prediction model reveals heterogeneity, stemness and immune microenvironment of UCEC

Background: Endometrial cancer (UCEC) is a highly heterogeneous gynecologic malignancy that exhibits variable prognostic outcomes and responses to immunotherapy. The Familial sequence similarity (FAM) gene family is known to contribute to the pathogenesis of various malignancies, but the extent of their involvement in UCEC has not been systematically studied. This investigation aimed to develop a robust risk profile based on FAM family genes (FFGs) to predict the prognosis and suitability for immunotherapy in UCEC patients. Methods: Using the TCGA-UCEC cohort from The Cancer Genome Atlas (TCGA) database, we obtained expression profiles of FFGs from 552 UCEC and 35 normal samples, and analyzed the expression patterns and prognostic relevance of 363 FAM family genes. The UCEC samples were randomly divided into training and test sets (1:1), and univariate Cox regression analysis and Lasso Cox regression analysis were conducted to identify the differentially expressed genes (FAM13C, FAM110B, and FAM72A) that were significantly associated with prognosis. A prognostic risk scoring system was constructed based on these three gene characteristics using multivariate Cox proportional risk regression. The clinical potential and immune status of FFGs were analyzed using CiberSort, SSGSEA, and tumor immune dysfunction and rejection (TIDE) algorithms. qRT-PCR and IHC for detecting the expression levels of 3-FFGs. Results: Three FFGs, namely, FAM13C, FAM110B, and FAM72A, were identified as strongly associated with the prognosis of UCEC and effective predictors of UCEC prognosis. Multivariate analysis demonstrated that the developed model was an independent predictor of UCEC, and that patients in the low-risk group had better overall survival than those in the high-risk group. The nomogram constructed from clinical characteristics and risk scores exhibited good prognostic power. Patients in the low-risk group exhibited a higher tumor mutational load (TMB) and were more likely to benefit from immunotherapy. Conclusion: This study successfully developed and validated novel biomarkers based on FFGs for predicting the prognosis and immune status of UCEC patients. The identified FFGs can accurately assess the prognosis of UCEC patients and facilitate the identification of specific subgroups of patients who may benefit from personalized treatment with immunotherapy and chemotherapy.

The ferroptosis signature predicts the prognosis and immune microenvironment of nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is a cancer with a high metastatic rate and poor prognosis. Growing studies suggest that ferroptosis take part in the development of tumours. At the same time, the connection between ferroptosis-related genes (FRGs) and the prognosis of NPC remains unclear. In this study, we explored the dysregulated FRGs between normal control and tumour samples of NPC. Firstly, 14 of 36 differentially expressed FRGs were identified in NPC tissues compared to normal tissues, among which ABCC1, GLS2, CS and HMGCR were associated with poor prognosis for patients. The four ferroptosis genes were used for consensus cluster analysis and two risk-related FRGs (ABCC1 and GLS2) were used in a risk model. The ROC curve revealed the good predictive performance of this risk signature. Multivariate analysis revealed that risk score and intratumoral TILs were independent risk factors linked to prognosis. Additionally, our results suggested that the risk signature was attached to the immune microenvironment. Moreover, the NPC patients with high risk were sensitive to chemotherapeutic drugs including axitinib, docetaxel, embelin, epothilone.B, parthenolide, thapsigargin, tipifarnib, vinorelbine. Finally, the expression of ABCC1 and GLS2 was validated in NPC tissues using immunohistochemistry. Together, these results revealed ferroptosis may be a potential biomarker in NPC and representing a promising future direction in prognosis and therapeutic strategy for the treatment of NPC.

Towards precision oncology with patient-derived xenografts

Under the selective pressure of therapy, tumours dynamically evolve multiple adaptive mechanisms that make static interrogation of genomic alterations insufficient to guide treatment decisions. Clinical research does not enable the assessment of how various regulatory circuits in tumours are affected by therapeutic insults over time and space. Likewise, testing different precision oncology approaches informed by composite and ever-changing molecular information is hard to achieve in patients. Therefore, preclinical models that incorporate the biology and genetics of human cancers, facilitate analyses of complex variables and enable adequate population throughput are needed to pinpoint randomly distributed response predictors. Patient-derived xenograft (PDX) models are dynamic entities in which cancer evolution can be monitored through serial propagation in mice. PDX models can also recapitulate interpatient diversity, thus enabling the identification of response biomarkers and therapeutic targets for molecularly defined tumour subgroups. In this Review, we discuss examples from the past decade of the use of PDX models for precision oncology, from translational research to drug discovery. We elaborate on how and to what extent preclinical observations in PDX models have confirmed and/or anticipated findings in patients. Finally, we illustrate emerging methodological efforts that could broaden the application of PDX models by honing their predictive accuracy or improving their versatility.

Small-molecule targeted therapies induce dependence on DNA double-strand break repair in residual tumor cells

Residual cancer cells that survive drug treatments with targeted therapies act as a reservoir from which eventual resistant disease emerges. Although there is great interest in therapeutically targeting residual cells, efforts are hampered by our limited knowledge of the vulnerabilities existing in this cell state. Here, we report that diverse oncogene-targeted therapies, including inhibitors of epidermal growth factor receptor (EGFR), anaplastic lymphoma kinase (ALK), KRAS, and BRAF, induce DNA double-strand breaks and, consequently, ataxia-telangiectasia mutated (ATM)-dependent DNA repair in oncogene-matched residual tumor cells. This DNA damage response, observed in cell lines, mouse xenograft models, and human patients, is driven by a pathway involving the activation of caspases 3 and 7 and the downstream caspase-activated deoxyribonuclease (CAD). CAD is, in turn, activated through caspase-mediated degradation of its endogenous inhibitor, ICAD. In models of EGFR mutant non-small cell lung cancer (NSCLC), tumor cells that survive treatment with small-molecule EGFR-targeted therapies are thus synthetically dependent on ATM, and combined treatment with an ATM kinase inhibitor eradicates these cells in vivo. This led to more penetrant and durable responses in EGFR mutant NSCLC mouse xenograft models, including those derived from both established cell lines and patient tumors. Last, we found that rare patients with EGFR mutant NSCLC harboring co-occurring, loss-of-function mutations in ATM exhibit extended progression-free survival on first generation EGFR inhibitor therapy relative to patients with EGFR mutant NSCLC lacking deleterious ATM mutations. Together, these findings establish a rationale for the mechanism-based integration of ATM inhibitors alongside existing targeted therapies.

Experiments from unfinished Registered Reports in the Reproducibility Project: Cancer Biology

As part of the Reproducibility Project: Cancer Biology, we published Registered Reports that described how we intended to replicate selected experiments from 29 high-impact preclinical cancer biology papers published between 2010 and 2012. Replication experiments were completed and Replication Studies reporting the results were submitted for 18 papers, of which 17 were accepted and published by eLife with the rejected paper posted as a preprint. Here, we report the status and outcomes obtained for the remaining 11 papers. Four papers initiated experimental work but were stopped without any experimental outcomes. Two papers resulted in incomplete outcomes due to unanticipated challenges when conducting the experiments. For the remaining five papers only some of the experiments were completed with the other experiments incomplete due to mundane technical or unanticipated methodological challenges. The experiments from these papers, along with the other experiments attempted as part of the Reproducibility Project: Cancer Biology, provides evidence about the challenges of repeating preclinical cancer biology experiments and the replicability of the completed experiments.

Massively parallel quantification of phenotypic heterogeneity in single-cell drug responses

Single-cell analysis tools have made substantial advances in characterizing genomic heterogeneity; however, tools for measuring phenotypic heterogeneity have lagged due to the increased difficulty of handling live biology. Here, we report a single-cell phenotyping tool capable of measuring image-based clonal properties at scales approaching 100,000 clones per experiment. These advances are achieved by exploiting a previously unidentified flow regime in ladder microfluidic networks that, under appropriate conditions, yield a mathematically perfect cell trap. Machine learning and computer vision tools are used to control the imaging hardware and analyze the cellular phenotypic parameters within these images. Using this platform, we quantified the responses of tens of thousands of single cell–derived acute myeloid leukemia (AML) clones to targeted therapy, identifying rare resistance and morphological phenotypes at frequencies down to 0.05%. This approach can be extended to higher-level cellular architectures such as cell pairs and organoids and on-chip live-cell fluorescence assays.

A mathematical study on the drug resistant virus emergence with HIV/AIDS treatment cases

HIV/AIDS drug treatments, one of which is highly active anti-retroviral ther-apy (HAART), often fail by the emergence of drug resistant virus. In this paper we study a quantitative method to evaluate the chance of resistant virus gen-eration. To this end we develop a mathematical description of the possibility of the emergence of resistant virus species against drug treatments, depend-ing on the trajectories of the state variables of HIV infection dynamic model. By simulation studies of mathematical models we apply the proposed analy-sis method to HIV/AIDS drug therapies, improved gradual dosage reduction (iGDR) and structured treatment interruption (STI). Based on the analysis it can be explained the reason why STI therapy often fails. Moreover it is con-cluded that iGDR is desirable particularly by decreasing the threat of resistant virus emergence.

Therapy-Induced Evolution of Human Lung Cancer Revealed by Single-Cell RNA Sequencing

Lung cancer, the leading cause of cancer mortality, exhibits heterogeneity that enables adaptability, limits therapeutic success, and remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) of metastatic lung cancer was performed using 49 clinical biopsies obtained from 30 patients before and during targeted therapy. Over 20,000 cancer and tumor microenvironment (TME) single-cell profiles exposed a rich and dynamic tumor ecosystem. scRNA-seq of cancer cells illuminated targetable oncogenes beyond those detected clinically. Cancer cells surviving therapy as residual disease (RD) expressed an alveolar-regenerative cell signature suggesting a therapy-induced primitive cell-state transition, whereas those present at on-therapy progressive disease (PD) upregulated kynurenine, plasminogen, and gap-junction pathways. Active T-lymphocytes and decreased macrophages were present at RD and immunosuppressive cell states characterized PD. Biological features revealed by scRNA-seq were biomarkers of clinical outcomes in independent cohorts. This study highlights how therapy-induced adaptation of the multi-cellular ecosystem of metastatic cancer shapes clinical outcomes.

•

scRNA-seq is feasible in metastatic human NSCLCs and reveals a rich tumor ecosystem

•

Individual tumors and cancer cells exhibit substantial molecular diversity

•

Cancer and tumor microenvironment cells exhibit marked therapy-induced plasticity

•

scRNA-seq of metastatic NSCLCs unveils new opportunities to improve clinical outcomes

Acquired Resistance in Lung Cancer

The last decade has witnessed a transformation in the treatment of advanced-stage lung cancer from a largely palliative approach to one where long-term durable remissions and even cures might be within reach. In this review, we discuss the current state of oncogene-directed precision medicine therapies in lung cancer and focus on the major cause of mortality for lung cancer patients: acquired resistance. We consider the multifaceted resistance mechanisms tumors utilize, often simultaneously. We then present areas for future scientific and clinical investigation with an emphasis on population dynamics, early detection, combinatorial therapies targeting resistance mechanisms, and understanding the drug-tolerant persister state.

Comprehensive analysis of the association between tumor glycolysis and immune/inflammation function in breast cancer

Background

Metabolic reprogramming, immune evasion and tumor-promoting inflammation are three hallmarks of cancer that provide new perspectives for understanding the biology of cancer. We aimed to figure out the relationship of tumor glycolysis and immune/inflammation function in the context of breast cancer, which is significant for deeper understanding of the biology, treatment and prognosis of breast cancer.

Methods

Using mRNA transcriptome data, tumor-infiltrating lymphocytes (TILs) maps based on digitized H&E-stained images and clinical information of breast cancer from The Cancer Genome Atlas projects (TCGA), we explored the expression and prognostic implications of glycolysis-related genes, as well as the enrichment scores and dual role of different immune/inflammation cells in the tumor microenvironment. The relationship between glycolysis activity and immune/inflammation function was studied by using the differential genes expression analysis, gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, gene set enrichment analyses (GSEA) and correlation analysis.

Results

Most glycolysis-related genes had higher expression in breast cancer compared to normal tissue. Higher phosphoglycerate kinase 1 (PGK1) expression was associated with poor prognosis. High glycolysis group had upregulated immune/inflammation-related genes expression, upregulated immune/inflammation pathways especially IL-17 signaling pathway, higher enrichment of multiple immune/inflammation cells such as Th2 cells and macrophages. However, high glycolysis group was associated with lower infiltration of tumor-killing immune cells such as NKT cells and higher immune checkpoints expression such as PD-L1, CTLA4, FOXP3 and IDO1.

Conclusions

In conclusion, the enhanced glycolysis activity of breast cancer was associated with pro-tumor immunity. The interaction between tumor glycolysis and immune/inflammation function may be mediated through IL-17 signaling pathway.

Treating Cancer as an Invasive Species

To cure a patient's cancer is to eradicate invasive cells from the ecosystem of the body. However, the ecologic complexity of this challenge is not well understood. Here we show how results from eradications of invasive mammalian species from islands-one of the few contexts in which invasive species have been regularly cleared-inform new research directions for treating cancer. We first summarize the epidemiologic characteristics of island invader eradications and cancer treatments by analyzing recent datasets from the Database of Invasive Island Species Eradications and The Cancer Genome Atlas, detailing the superior successes of island eradication projects. Next, we compare how genetic and environmental factors impact success in each system. These comparisons illuminate a number of promising cancer research and treatment directions, such as heterogeneity engineering as motivated by gene drives and adaptive therapy; multiscale analyses of how population heterogeneity potentiates treatment resistance; and application of ecological data mining techniques to high-throughput cancer data. We anticipate that interdisciplinary comparisons between tumor progression and invasive species would inspire development of novel paradigms to cure cancer.