ScTrio-seq identifies colon cancer lineages: To better design treatments for cancer, it is important to understand the heterogeneity in tumors and how this contributes to metastasis. To examine this process, Bian et al. used a single-cell triple omics sequencing (scTrio-seq) technique to examine the mutations, transcriptome, and methylome within colorectal cancer tumors and metastases from 10 individual patients. The analysis provided insights into tumor evolution, linked DNA methylation to genetic lineages, and showed that DNA methylation levels are consistent within lineages but can differ substantially among clones. Science , this issue p. 1060

Colorectal cancer (CRC) is a highly heterogeneous cancer and exploring novel therapeutic options is a pressing issue that needs to be addressed. Here, we established human CRC tumor-derived organoids that well represent both morphological and molecular heterogeneities of original tumors. To efficiently identify repurposed drugs for CRC, we developed a robust organoid-based drug screening system. By combining the repurposed drug library and computation-based drug prediction, 335 drugs were tested and 34 drugs with anti-CRC effects were identified. More importantly, we conducted a detailed transcriptome analysis of drug responses and divided the drug response signatures into five representative patterns: differentiation induction, growth inhibition, metabolism inhibition, immune response promotion, and cell cycle inhibition. The anticancer activities of drug candidates were further validated in the established patient-derived organoids-based xenograft (PDOX) system in vivo. We found that fedratinib, trametinib, and bortezomib exhibited effective anticancer effects. Furthermore, the concordance and discordance of drug response signatures between organoids in vitro and pairwise PDOX in vivo were evaluated. Our study offers an innovative approach for drug discovery, and the representative transcriptome features of drug responses provide valuable resources for developing novel clinical treatments for CRC.

Human gastric cancer is a highly lethal disease, but the underlying multiomic molecular signatures remain largely unclear. Here, we performed multi-regional sampling, parallel single-cell multiomics sequencing and integrated analyses of human gastric cancer. We identified common transcriptomic alterations of gastric cancer cells, such as aberrant down-regulation of genes associated with normal stomach function and up-regulation of KRT7, PI3, S100A4, etc. Surprisingly, aberrant and prevalent up-regulation of genes highly-expressed in normal colorectal epithelial cells were also identified in cancer cells, which may be partially regulated by promoter chromatin accessibility and DNA methylation levels. We revealed the single-cell DNA methylome landscape of gastric cancer, and identified candidate DNA methylation biomarkers, such as hypermethylated promoters of TMEM240 and HAGLROS, and hypomethylated promoters of TRPM2-AS and HRH1. Additionally, the relationships between genetic lineages, DNA methylation and transcriptomic clusters were systematically revealed at single-cell level. We showed that DNA methylation heterogeneities were mainly among different genetic lineages of cancer cells. Moreover, we found that the DNA methylation levels of cancer cells with poorer differentiation states tend to be higher than those of cancer cells with better differentiation states in the primary tumor within the same patient, although still lower than in normal gastric epithelial cells. Cancer cells with poorer differentiation states also prevalently down-regulated MUC1 expression and immune-related pathways, and had poor infiltration of CD8⁺ T cells. Our study dissected molecular signatures of intratumoral heterogeneities and differentiation states of human gastric cancer using integrative single-cell multiomics analyses.

Colorectal cancer is a highly heterogeneous disease, with well-characterized subtypes based on genome, DNA methylome, and transcriptome signatures. To chart the epigenetic landscape of colorectal cancers, we generated a high-quality single-cell chromatin accessibility atlas of epithelial cells for 29 patients. Abnormal chromatin states acquired in adenomas were largely retained in colorectal cancers, which were tightly accompanied by opposite changes of DNA methylation. Unsupervised analysis on malignant cells revealed two epigenetic subtypes, exactly matching the iCMS classification, and key iCMS-specific transcription factors (TFs) were identified, including HNF4A and PPARA for iCMS2 tumors and FOXA3 and MAFK for iCMS3 tumors. Notably, subtype-specific TFs bind to distinct target gene sets and contribute to both interpatient similarities and diversities for both chromatin accessibilities and RNA expressions. Moreover, we identified CpG-island methylator phenotypes and pinpointed chromatin state signatures and TF regulators for the CIMP-high subtype. Our work systematically revealed the epigenetic basis of the well-known iCMS and CIMP classifications of colorectal cancers. Significance: Our work revealed the epigenetic basis of the well-known iCMS and CIMP classifications of colorectal cancers. Moreover, interpatient minor similarities and major diversities of chromatin accessibility signatures of TF target genes can faithfully explain the corresponding interpatient minor similarities and major diversities of RNA expression signatures of colorectal cancers, respectively. This article is featured in Selected Articles from This Issue, p. 897

Background: Heterogeneity in colorectal cancer (CRC) patients provides novel strategies in clinical decision-making. Identifying distinctive subgroups in patients can improve the screening of CRC and reduce the cost of tests. Metabolism-related long non-coding RNA (lncRNA) can help detection of tumorigenesis and development for CRC patients.Methods: RNA sequencing and clinical data of CRC patients which extracted and integrated from public databases including The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were set as training cohort and validation cohort. Metabolism-related genes were acquired from Kyoto Encyclopedia of Genes and Genomes (KEGG) and the metabolism-related lncRNAs were filtered using correlation analysis. The risk score was calculated based on lncRNAs with prognostic value and verified through survival curve, receiver operating characteristic (ROC) curve and risk curve. Prognostic factors of CRC patients were also analyzed. Nomogram was constructed based on the results of cox regression analyses. The different immune status was observed in the single sample Gene Set Enrichment Analysis (ssGSEA).Results: The training cohort and the validation cohort enrolled 432 and 547 CRC patients respectively. A total of 23 metabolism-related lncRNAs with prognostic value were screened out and 10 of which were significantly differentially expressed between tumour and normal tissues. Finally, 8 lncRNAs were used to establish a risk score (DICER1-AS1, PCAT6, GAS5, PRR7-AS1, MCM3AP-AS1, GAS6-AS1, LINC01082 and ADIRF-AS1). Patients were divided into high-risk and low-risk groups according to the median of risk scores in training cohort and the survival curves indicated that the survival prognosis was significantly different. The area under curve (AUC) of the ROC curve in two cohorts were both greater than 0.6. The age, tumour stage and risk score were selected as independent factors and used to construct a nomogram to predict CRC patients' survival rate with the c-index of 0.806. The ssGSEA indicated that the risk score was associated with immune cells and functions.Conclusions: Our systematic study established a metabolism-related lncRNA signature to predict outcomes of CRC patients which may contribute to individual prevention and treatment.