Bulk and Single-Cell Next-Generation Sequencing: Individualizing Treatment for Colorectal Cancer

Molecular subtypes of colorectal cancer in the era of precision oncotherapy: Current inspirations and future challenges

BACKGROUND

Colorectal cancer (CRC) is among the most hackneyed malignancies. Even patients with identical clinical symptoms and the same TNM stage still exhibit radically different clinical outcomes after receiving equivalent treatment regimens, indicating extensive heterogeneity of CRC. Myriad molecular subtypes of CRC have been exploited for decades, including the most compelling consensus molecular subtype (CMS) classification that has been broadly applied for patient stratification and biomarker-drug combination formulation. Encountering barriers to clinical translation, however, CMS classification fails to fully reflect inter- or intra-tumor heterogeneity of CRC. As a consequence, addressing heterogeneity and precisely managing CRC patients with unique characteristics remain arduous tasks for clinicians.

REVIEW

In this review, we systematically summarize molecular subtypes of CRC and further elaborate on their clinical applications, limitations, and future orientations.

CONCLUSION

In recent years, exploration of subtypes through cell lines, animal models, patient-derived xenografts (PDXs), organoids, and clinical trials contributes to refining biological insights and unraveling subtype-specific therapies in CRC. Therapeutic interventions including nanotechnology, clustered regulatory interspaced short palindromic repeat/CRISPR-associated nuclease 9 (CRISPR/Cas9), gut microbiome, and liquid biopsy are powerful tools with the possibility to shift the immunologic landscape and outlook for CRC precise medicine.

Unveiling the best predictive models for early‑onset metastatic cancer: Insights and innovations (Review)

Cancer metastasis is the primary cause of cancer deaths. Metastasis involves the spread of cancer cells from the primary tumors to other body parts, commonly through lymphatic and vascular pathways. Key aspects include the high mutation rate and the capability of metastatic cells to form invasive tumors even without a large initial tumor mass. Particular emphasis is given to early metastasis, occurring in initial cancer stages and often leading to misdiagnosis, which adversely affects survival and prognosis. The present review highlighted the need for improved understanding and detection methods for early metastasis, which has not been effectively identified clinically. The present review demonstrated the clinicopathological and molecular characteristics of early‑onset metastatic types of cancer, noting factors such as age, race, tumor size and location as well as the histological and pathological grade as significant predictors. In conclusion, the present review underscored the importance of early detection and management of metastatic types of cancer and called for improved predictive models, including advanced techniques such as nomograms and machine learning, so as to enhance patient outcomes, acknowledging the challenges and limitations of the current research as well as the necessity for further studies.

Clinical implications of single cell sequencing for bladder cancer

Bladder cancer (BC) is the 10th most common cancer worldwide, with about 0.5 million reported new cases and about 0.2 million deaths per year. In this scoping review, we summarize the current evidence regarding the clinical implications of single-cell sequencing for bladder cancer based on PRISMA guidelines. We searched PubMed, CENTRAL, Embase, and supplemented with manual searches through the Scopus, and Web of Science for published studies until February 2023. We included original studies that used at least one single-cell technology to study bladder cancer. Forty-one publications were included in the review. Twenty-nine studies showed that this technology can identify cell subtypes in the tumor microenvironment that may predict prognosis or response to immune checkpoint inhibition therapy. Two studies were able to diagnose BC by identifying neoplastic cells through single-cell sequencing urine samples. The remaining studies were mainly a preclinical exploration of tumor microenvironment at single cell level. Single-cell sequencing technology can discriminate heterogeneity in bladder tumor cells and determine the key molecular properties that can lead to the discovery of novel perspectives on cancer management. This nascent tool can advance the early diagnosis, prognosis judgment, and targeted therapy of bladder cancer.

Single cell metabolism: current and future trends

Single cell metabolomics is an emerging and rapidly developing field that complements developments in single cell analysis by genomics and proteomics. Major goals include mapping and quantifying the metabolome in sufficient detail to provide useful information about cellular function in highly heterogeneous systems such as tissue, ultimately with spatial resolution at the individual cell level. The chemical diversity and dynamic range of metabolites poses particular challenges for detection, identification and quantification. In this review we discuss both significant technical issues of measurement and interpretation, and progress toward addressing them, with recent examples from diverse biological systems. We provide a framework for further directions aimed at improving workflow and robustness so that such analyses may become commonly applied, especially in combination with metabolic imaging and single cell transcriptomics and proteomics.

Capturing tumour heterogeneity in pre- and post-chemotherapy colorectal cancer ascites-derived cells using single-cell RNA-sequencing

Malignant ascites is an abnormal accumulation of fluid within the peritoneal cavity, caused by metastasis of several types of cancers, including colorectal cancer (CRC). Cancer cells in ascites reflect poor prognosis and serve as a good specimen to study tumour heterogeneity, as they represent a collection of multiple metastatic sites in the peritoneum. In the present study, we have employed single-cell RNA-sequencing (scRNA-seq) to explore and characterise ascites-derived cells from a CRC patient. The samples were prepared using mechanical and enzymatic dissociations, and obtained before and after a chemotherapy treatment. Unbiased clustering of 19,653 cells from four samples reveals 14 subclusters with unique transcriptomic patterns in four major cell types: epithelial cells, myeloid cells, fibroblasts, and lymphocytes. Interestingly, the percentages of cells recovered from different cell types appeared to be influenced by the preparation protocols, with more than 90% reduction in the number of myeloid cells recovered by enzymatic preparation. Analysis of epithelial cell subpopulations unveiled only three out of eleven subpopulations with clear contraction after the treatment, suggesting that the majority of the heterogeneous ascites-derived cells were resistant to the treatment, potentially reflecting the poor treatment outcome observed in the patient. Overall, our study showcases highly heterogeneous cancer subpopulations at single-cell resolution, which respond differently to a particular chemotherapy treatment. All in all, this work highlights the potential benefit of single-cell analyses in planning appropriate treatments and real-time monitoring of therapeutic response in cancer patients through routinely discarded ascites samples.

Single-cell RNA sequencing and bioinformatics as tools to decipher cancer heterogenicity and mechanisms of drug resistance

It is well known that cancer is an aggressive disease, often associated with relapse, in many cases due to drug resistance. Cancer stem cell and clonal evolution are frequently causes of innate or acquired drug resistance. Current RNA sequencing technologies do not distinguish gene expression of different cell lineages because they are based on bulk cell studies. Single-cell RNA sequencing technologies and related bioinformatics clustering and differential expression analysis represent a turning point in cancer research. They are emerging as essential tools for dissecting tumors at single-cell resolution and represent novel tools to understand carcinogenesis and drug response. In this review, we will outline the role of these new technologies in addressing cancer heterogeneity and cell lineage-dependent drug resistance.

Touch Imprint Intraoperative Flow Cytometry as a Complementary Tool for Detailed Assessment of Resection Margins and Tumor Biology in Liver Surgery for Primary and Metastatic Liver Neoplasms

Liver resection is the main treatment for primary and metastatic liver tumors in order to achieve long-term survival with good quality of life. The ultimate goal of surgical oncology is to achieve complete tumor removal with adequate clear surgical margins. Flow cytometry is a powerful analytical technique with applications such as phenotypic analysis and quantification of DNA content. Intraoperative flow cytometry (iFC) is the application of flow cytometry for DNA content/ploidy and cell cycle distribution analysis during surgery for tumor cell analysis and margin evaluation. It has been used for cell analysis of intracranial tumors and recently of head and neck carcinomas and breast carcinomas, as well as for tumor margin evaluation. Herein, we present a novel touch imprint iFC protocol for the detailed assessment of tumor margins during excision of malignant hepatic lesions. The protocol aims to offer information on surgical margins after removal of malignant liver tumors based on DNA content of cancer cells and to corroborate the results of iFC with that of histopathological analysis. Based on the established role of iFC in other types of malignancies, our specialized protocol has the potential, through characterization of cells in liver transection surface post hepatectomy, to offer significant information on the type of resection and tumor biology. This information can be used to effectively guide intra- and postoperative patient management.

A BAYESIAN NONPARAMETRIC MODEL FOR INFERRING SUBCLONAL POPULATIONS FROM STRUCTURED DNA SEQUENCING DATA

There are distinguishing features or "hallmarks" of cancer that are found across tumors, individuals, and types of cancer, and these hallmarks can be driven by specific genetic mutations. Yet, within a single tumor there is often extensive genetic heterogeneity as evidenced by single-cell and bulk DNA sequencing data. The goal of this work is to jointly infer the underlying genotypes of tumor subpopulations and the distribution of those subpopulations in individual tumors by integrating single-cell and bulk sequencing data. Understanding the genetic composition of the tumor at the time of treatment is important in the personalized design of targeted therapeutic combinations and monitoring for possible recurrence after treatment. We propose a hierarchical Dirichlet process mixture model that incorporates the correlation structure induced by a structured sampling arrangement and we show that this model improves the quality of inference. We develop a representation of the hierarchical Dirichlet process prior as a Gamma-Poisson hierarchy and we use this representation to derive a fast Gibbs sampling inference algorithm using the augment-and-marginalize method. Experiments with simulation data show that our model outperforms standard numerical and statistical methods for decomposing admixed count data. Analyses of real acute lymphoblastic leukemia cancer sequencing dataset shows that our model improves upon state-of-the-art bioinformatic methods. An interpretation of the results of our model on this real dataset reveals co-mutated loci across samples.

The roles and prognostic significance of ABI1-TSV-11 expression in patients with left-sided colorectal cancer

Abnormally expressed and/or phosphorylated Abelson interactor 1 (ABI1) participates in the metastasis and progression of colorectal cancer (CRC). ABI1 presents as at least 12 transcript variants (TSVs) by mRNA alternative splicing, but it is unknown which of them is involved in CRC metastasis and prognosis. Here, we firstly identified ABI1-TSV-11 as a key TSV affecting the metastasis and prognosis of left-sided colorectal cancer (LsCC) and its elevated expression is related to lymph node metastasis and shorter overall survival (OS) in LsCC by analyzing data from The Cancer Genome Atlas and TSVdb. Secondly, ABI1-TSV-11 overexpression promoted LoVo and SW480 cells adhesion and migration in vitro, and accelerated LoVo and SW480 cells lung metastasis in vivo. Finally, mechanism investigations revealed that ABI1-isoform-11 interacted with epidermal growth factor receptor pathway substrate 8 (ESP8) and regulated actin dynamics to affect LoVo and SW480 cells biological behaviors. Taken together, our data demonstrated that ABI1-TSV-11 plays an oncogenic role in LsCC, it is an independent risk factor of prognosis and may be a potential molecular marker and therapeutic target in LsCC.

Pathologic evolution-related Gene Analysis based on both single-cell and bulk transcriptomics in Colorectal Cancer

Purpose: The patients diagnosed with colorectal cancer (CRC) are likely to undergo differential outcomes in clinical survival owing to different pathologic stages. However, signatures in association with pathologic evolution and CRC prognosis are not clearly defined. This study aimed to identify pathologic evolution-related genes in CRC based on both single-cell and bulk transcriptomics.

Patients and methods: The CRC single-cell transcriptomic dataset (GSE81861, n=590) with clinical information and tumor microenvironmental tissues was collected to identify the pathologic evolution-related genes. The colonic adenocarcinoma and rectum adenocarcinoma transcriptomics from The Cancer Genome Atlas were obtained as the training dataset (n=363) and 5 other CRC transcriptomics cohorts from Gene Expression Omnibus (n=1031) were acquired as validation data. Graph-based clustering analysis algorithm was applied to identify pathologic evolution-related cell populations. Pseudotime analysis was performed to construct the trajectory plot of pathologic evolution and to define hub genes in the evolution process. Cell-type identification by estimating relative subsets of RNA transcripts was then executed to build a novel cell infiltration classifier. The prediction efficacy of this classifier was validated in bulk transcriptomic datasets.

Results: Epithelial and T cells were elucidated to be related to the pathologic stages in CRC tissues. Pseudotime analysis and survival analysis indicated that HOXC5, HOXC8 and BMP5 were the marker genes in pathologic evolution process. Our cell infiltration classifier exhibited excellent forecast efficacy in predicting pathologic stages and prognosis of CRC patients.

Conclusion: We identified pathologic evolution-related genes in single-cell transcriptomic and proposed a novel specific cell infiltration classifier to forecast the prognosis of CRC patients based on pathologic stage-related hub genes HOXC6, HOXC8 and BMP5.

Analysis of prognostic genes in the tumor microenvironment of lung adenocarcinoma

Background

Prognostic genes in the tumor microenvironment play an important role in immune biological processes and the response of cancer to immunotherapy. Thus, we aimed to assess new biomarkers that are associated with immune/stromal cells in lung adenocarcinomas (LUAD) using the ESTIMATE algorithm, which also significantly affects the prognosis of cancer.

Methods

The RNA sequencing (RNA-Seq) and clinical data of LUAD were downloaded from the the Cancer Genome Atlas (TCGA ). The immune and stromal scores were calculated for each sample using the ESTIMATE algorithm. The LUAD gene chip expression profile data and the clinical data (GSE37745, GSE11969, and GSE50081) were downloaded from the Gene Expression Omnibus (GEO) for subsequent validation analysis. Differentially expressed genes were calculated between high and low score groups. Univariate Cox regression analysis was performed on differentially expressed genes (DEGs) between the two groups to obtain initial prognosis genes. These were verified by three independent LUAD cohorts from the GEO database. Multivariate Cox regression was used to identify overall survival-related DEGs. UALCAN and the Human Protein Atlas were used to analyze the mRNA /protein expression levels of the target genes. Immune cell infiltration was evaluated using the Tumor Immune Estimation Resource (TIMER) and CIBERSORT methods, and stromal cell infiltration was assessed using xCell.

Results

In this study, immune scores and stromal scores are significantly associated with the clinical characteristics of LUAD, including T stage, M stage, pathological stage, and overall survival time. 530 DEGs (18 upregulated and 512 downregulated) were found to coexist in the difference analysis with the immune scores and stromal scores subgroup. Univariate Cox regression analysis showed that 286 of the 530 DEGs were survival-related genes (p < 0.05). Of the 286 genes initially identified, nine prognosis-related genes (CSF2RB, ITK, FLT3, CD79A, CCR4, CCR6, DOK2, AMPD1, and IGJ) were validated from three separate LUAD cohorts. In addition, functional analysis of DEGs also showed that various immunoregulatory molecular pathways, including regulation of immune response and the chemokine signaling pathways, were involved. Five genes (CCR6, ITK, CCR4, DOK2, and AMPD1) were identified as independent prognostic indicators of LUAD in specific data sets. The relationship between the expression levels of these genes and immune genes was assessed. We found that CCR6 mRNA and protein expression levels of LUAD were greater than in normal tissues. We evaluated the infiltration of immune cells and stromal cells in groups with high and low levels of expression of CCR6 in the TCGA LUAD cohort. In summary, we found a series of prognosis-related genes that were associated with the LUAD tumor microenvironment.

Resolving Metabolic Heterogeneity in Experimental Models of the Tumor Microenvironment from a Stable Isotope Resolved Metabolomics Perspective

The tumor microenvironment (TME) comprises complex interactions of multiple cell types that determines cell behavior and metabolism such as nutrient competition and immune suppression. We discuss the various types of heterogeneity that exist in solid tumors, and the complications this invokes for studies of TME. As human subjects and in vivo model systems are complex and difficult to manipulate, simpler 3D model systems that are compatible with flexible experimental control are necessary for studying metabolic regulation in TME. Stable Isotope Resolved Metabolomics (SIRM) is a valuable tool for tracing metabolic networks in complex systems, but at present does not directly address heterogeneous metabolism at the individual cell level. We compare the advantages and disadvantages of different model systems for SIRM experiments, with a focus on lung cancer cells, their interactions with macrophages and T cells, and their response to modulators in the immune microenvironment. We describe the experimental set up, illustrate results from 3D cultures and co-cultures of lung cancer cells with human macrophages, and outline strategies to address the heterogeneous TME.