Alarmin S100A8 imparts chemoresistance of esophageal cancer by reprogramming cancer-associated fibroblasts.

Chemotherapy remains the first-line treatment for advanced esophageal cancer. However, durable benefits are achieved by only a limited subset of individuals due to the elusive chemoresistance. Here, we utilize patient-derived xenografts (PDXs) from esophageal squamous-cell carcinoma to investigate chemoresistance mechanisms in preclinical settings. We observe that activated cancer-associated fibroblasts (CAFs) are enriched in the tumor microenvironment of PDXs resistant to chemotherapy. Mechanistically, we reveal that cancer-cell-derived S100A8 triggers the intracellular RhoA-ROCK-MLC2-MRTF-A pathway by binding to the CD147 receptor of CAFs, inducing CAF polarization and leading to chemoresistance. Therapeutically, we demonstrate that blocking the S100A8-CD147 pathway can improve chemotherapy efficiency. Prognostically, we found the S100A8 levels in peripheral blood can serve as an indicator of chemotherapy responsiveness. Collectively, our study offers a comprehensive understanding of the molecular mechanisms underlying chemoresistance in esophageal cancer and highlights the potential value of S100A8 in the clinical management of esophageal cancer.

Lycopene Promotes Osteogenesis and Reduces Adipogenesis through Regulating FoxO1/PPARγ Signaling in Ovariectomized Rats and Bone Marrow Mesenchymal Stem Cells.

Recent interest in preventing the development of osteoporosis has focused on the regulation of redox homeostasis. However, the action of lycopene (LYC), a strong natural antioxidant compound, on osteoporotic bone loss remains largely unknown. Here, we show that oral administration of LYC to OVX rats for 12 weeks reduced body weight gain, improved lipid metabolism, and preserved bone quality. In addition, LYC treatment inhibited ROS overgeneration in serum and bone marrow in OVX rats, and in BMSCs upon H2O2 stimulation, leading to inhibiting adipogenesis and promoting osteogenesis during bone remodeling. At the molecular level, LYC improved bone quality via an increase in the expressions of FoxO1 and Runx2 and a decrease in the expressions of PPARγ and C/EBPα in OVX rats and BMSCs. Collectively, these findings suggest that LYC attenuates osteoporotic bone loss through promoting osteogenesis and inhibiting adipogenesis via regulation of the FoxO1/PPARγ pathway driven by oxidative stress, presenting a novel strategy for osteoporosis management.

CodLncScape Provides a Self-Enriching Framework for the Systematic Collection and Exploration of Coding LncRNAs.

Recent studies have revealed that numerous lncRNAs can translate proteins under specific conditions, performing diverse biological functions, thus termed coding lncRNAs. Their comprehensive landscape, however, remains elusive due to this field's preliminary and dispersed nature. This study introduces codLncScape, a framework for coding lncRNA exploration consisting of codLncDB, codLncFlow, codLncWeb, and codLncNLP. Specifically, it contains a manually compiled knowledge base, codLncDB, encompassing 353 coding lncRNA entries validated by experiments. Building upon codLncDB, codLncFlow investigates the expression characteristics of these lncRNAs and their diagnostic potential in the pan-cancer context, alongside their association with spermatogenesis. Furthermore, codLncWeb emerges as a platform for storing, browsing, and accessing knowledge concerning coding lncRNAs within various programming environments. Finally, codLncNLP serves as a knowledge-mining tool to enhance the timely content inclusion and updates within codLncDB. In summary, this study offers a well-functioning, content-rich ecosystem for coding lncRNA research, aiming to accelerate systematic studies in this field.

Cnidii Fructus: A traditional Chinese medicine herb and source of antiosteoporotic drugs.

BACKGROUND: Osteoporosis (OP) is a prevalent chronic metabolic bone disease for which limited countermeasures are available. Cnidii Fructus (CF), primarily derived from Cnidium monnieri (L.) Cusson., has been tested in clinical trials of traditional Chinese medicine for the management of OP. Accumulating preclinical studies indicate that CF may be used against OP. MATERIALS AND METHODS: Comprehensive documentation and analysis were conducted to retrieve CF studies related to its main phytochemical components as well as its pharmacokinetics, safety and pharmacological properties. We also retrieved information on the mode of action of CF and, in particular, preclinical and clinical studies related to bone remodeling. This search was performed from the inception of databases up to the end of 2022 and included PubMed, China National Knowledge Infrastructure, the National Science and Technology Library, the China Science and Technology Journal Database, Weipu, Wanfang, the Web of Science and the China National Patent Database. RESULTS: CF contains a wide range of natural active compounds, including osthole, bergapten, imperatorin and xanthotoxin, which may underlie its beneficial effects on improving bone metabolism and quality. CF action appears to be mediated via multiple processes, including the osteoprotegerin (OPG)/receptor activator of nuclear factor-κB ligand (RANKL)/receptor activator of nuclear factor-κB (RANK), Wnt/ß-catenin and bone morphogenetic protein (BMP)/Smad signaling pathways. CONCLUSION: CF and its ingredients may provide novel compounds for developing anti-OP drugs.

Genomic alterations driving precancerous to cancerous lesions in esophageal cancer development.

Esophageal squamous cell carcinoma (ESCC) develops through a series of increasingly abnormal precancerous lesions. Previous studies have revealed the striking differences between normal esophageal epithelium and ESCC in copy number alterations (CNAs) and mutations in genes driving clonal expansion. However, due to limited data on early precancerous lesions, the timing of these transitions and which among them are prerequisites for malignant transformation remained unclear. Here, we analyze 1,275 micro-biopsies from normal esophagus, early and late precancerous lesions, and esophageal cancers to decipher the genomic alterations at each stage. We show that the frequency of TP53 biallelic inactivation increases dramatically in early precancerous lesion stage while CNAs and APOBEC mutagenesis substantially increase at late stages. TP53 biallelic loss is the prerequisite for the development of CNAs of genes in cell cycle, DNA repair, and apoptosis pathways, suggesting it might be one of the earliest steps initiating malignant transformation.

Aberrant epithelial cell interaction promotes esophageal squamous-cell carcinoma development and progression.

Epithelial-mesenchymal transition (EMT) and proliferation play important roles in epithelial cancer formation and progression, but what molecules and how they trigger EMT is largely unknown. Here we performed spatial transcriptomic and functional analyses on samples of multistage esophageal squamous-cell carcinoma (ESCC) from mice and humans to decipher these critical issues. By investigating spatiotemporal gene expression patterns and cell-cell interactions, we demonstrated that the aberrant epithelial cell interaction via EFNB1-EPHB4 triggers EMT and cell cycle mediated by downstream SRC/ERK/AKT signaling. The aberrant epithelial cell interaction occurs within the basal layer at early precancerous lesions, which expands to the whole epithelial layer and strengthens along the cancer development and progression. Functional analysis revealed that the aberrant EFNB1-EPHB4 interaction is caused by overexpressed ΔNP63 due to TP53 mutation, the culprit in human ESCC tumorigenesis. Our results shed new light on the role of TP53-TP63/ΔNP63-EFNB1-EPHB4 axis in EMT and cell proliferation in epithelial cancer formation.

Collagen 1-mediated CXCL1 secretion in tumor cells activates fibroblasts to promote radioresistance of esophageal cancer.

Esophageal squamous-cell carcinoma (ESCC) is commonly treated with radiotherapy; however, radioresistance hinders its clinical effectiveness, and the underlying mechanism remains elusive. Here, we develop patient-derived xenografts (PDXs) from 19 patients with ESCC to investigate the mechanisms driving radioresistance. Using RNA sequencing, cytokine arrays, and single-cell RNA sequencing, we reveal an enrichment of cancer-associated fibroblast (CAF)-derived collagen type 1 (Col1) and tumor-cell-derived CXCL1 in non-responsive PDXs. Col1 not only promotes radioresistance by augmenting DNA repair capacity but also induces CXCL1 secretion in tumor cells. Additionally, CXCL1 further activates CAFs via the CXCR2-STAT3 pathway, establishing a positive feedback loop. Directly interfering with tumor-cell-derived CXCL1 or inhibiting the CXCL1-CXCR2 pathway effectively restores the radiosensitivity of radioresistant xenografts in vivo. Collectively, our study provides a comprehensive understanding of the molecular mechanisms underlying radioresistance and identifies potential targets to improve the efficacy of radiotherapy for ESCC.

Body color plasticity of Diaphorina citri reflects a response to environmental stress.

Body color polyphenism is common in Diaphorina citri. Previous studies compared physiological characteristics in D. citri, but the ecological and biological significance of its body color polyphenism remains poorly understood. We studied the ecological and molecular effects of stressors related to body color in D. citri. Crowding or low temperature induced a high proportion of gray morphs, which had smaller bodies, lower body weight, and greater susceptibility to the insecticide dinotefuran. We performed transcriptomic and metabolomics analysiis of 2 color morphs in D. citri. Gene expression dynamics revealed that the differentially expressed genes were predominantly involved in energy metabolism, including fatty acid metabolism, amino acid metabolism, and carbohydrate metabolism. Among these genes, plexin, glycosidase, phospholipase, take out, trypsin, and triacylglycerol lipase were differentially expressed in 2 color morphs, and 6 hsps (3 hsp70, hsp83, hsp90, hsp68) were upregulated in gray morphs. The metabolome data showed that blue morphs exhibited a higher abundance of fatty acid and amino acid, whereas the content of carbohydrates was elevated in gray morphs. This study partly explains the body color polyphenism of D. citri and provides insights into the molecular changes of stress response of D. citri.

Comprehensive Analyses Reveal Effects on Tumor Immune Infiltration and Immunotherapy Response of APOBEC Mutagenesis and Its Molecular Mechanisms in Esophageal Squamous Cell Carcinoma.

The apolipoprotein B mRNA editing enzyme catalytic polypeptide (APOBEC) mutagenesis is prevalent in esophageal squamous cell carcinoma (ESCC). However, the functional role of APOBEC mutagenesis has yet to be fully delineated. To address this, we collect matched multi-omics data of 169 ESCC patients and evaluate characteristics of immune infiltration using multiple bioinformatic approaches based on bulk and single-cell RNA sequencing (scRNA-seq) data and verified by functional assays. We find that APOBEC mutagenesis prolongs overall survival (OS) of ESCC patients. The reason for this outcome is probably due to high anti-tumor immune infiltration, immune checkpoints expression and immune related pathway enrichment, such as interferon (IFN) signaling, innate and adaptive immune system. The elevated AOBEC3A (A3A) activity paramountly contributes to the footprints of APOBEC mutagenesis and is first discovered to be transactivated by FOSL1. Mechanistically, upregulated A3A exacerbates cytosolic double-stranded DNA (dsDNA) accumulation, thus stimulating cGAS-STING pathway. Simultaneously, A3A is associated with immunotherapy response which is predicted by TIDE algorithm, validated in a clinical cohort and further confirmed in mouse models. These findings systematically elucidate the clinical relevance, immunological characteristics, prognostic value for immunotherapy and underlying mechanisms of APOBEC mutagenesis in ESCC, which demonstrate great potential in clinical utility to facilitate clinical decisions.

Simiao San alleviates hyperuricemia and kidney inflammation by inhibiting NLRP3 inflammasome and JAK2/STAT3 signaling in hyperuricemia mice.

ETHNOPHARMACOLOGICAL RELEVANCE: Simiao San (SmS), a famous traditional Chinese formula, is clinically used to treat patients with hyperuricemia (HUA). However, its mechanism of action on lowering uric acid (UA) and inhibiting inflammation still deserves further investigation. AIM OF THE STUDY: To examine the effect and its possible underlying mechanism of SmS on UA metabolism and kidney injury in HUA mouse. MATERIALS AND METHODS: The HUA mouse model was constructed with the combined administration of both potassium oxalate and hypoxanthine. The effects of SmS on UA, xanthine oxidase (XOD), creatinine (CRE), blood urea nitrogen (BUN), interleukin-10 (IL-10), interleukin-1ß (IL-1ß), interleukin-6 (IL-6), and tumor necrosis factor-α (TNF-α) were determined by ELISA or biochemical assays. Hematoxylin and eosin (H&E) was used to observe pathological alterations in the kidneys of HUA mice. The expression levels of organic anion transporter 1 (OAT1), recombinant urate transporter 1 (URAT1), glucose transporter 9 (GLUT9), nucleotide binding domain and leucine rich repeat pyrin domain containing 3 (NLRP3), Cleaved-Caspase 1, apoptosis-associated speck like protein (ASC), nuclear factor kappa-B (NF-κB), IL-6, janus kinase 2 (JAK2), phosphor (P)-JAK2, signal transducers and activators of transcription 3 (STAT3), P-STAT3, suppressor of cytokine signaling 3 (SOCS3) were examined by Western blot and/or immunohistochemical (IHC) staining. The major ingredients in SmS were identified by a HPLC-MS assay. RESULTS: HUA mouse exhibited an elevation in serum levels of UA, BUN, CRE, XOD, and the ratio of urinary albumin to creatinine (UACR), and a decline in urine levels of UA and CRE. In addition, HUA induces pro-inflammatory microenvironment in mouse, including an increase in serum levels of IL-1ß, IL-6, and TNF-α, and renal expressions of URAT1, GULT9, NLRP3, ASC, Cleaved-Caspase1, P-JAK2/JAK2, P-STAT3/STAT3, and SOCS3, and a decrease in serum IL-10 level and renal OAT1 expression as well as a disorganization of kidney pathological microstructure. In contrast, SmS intervention reversed these alterations in HUA mouse. CONCLUSION: SmS could alleviate hyperuricemia and renal inflammation in HUA mouse. The action mechanisms behind these alterations may be associated with a limitation of the NLRP3 inflammasome and JAK2/STAT3 signaling pathways.

Epithelial cells activate fibroblasts to promote esophageal cancer development.

Esophageal squamous-cell carcinoma (ESCC) develops through multistage epithelial cancer formation, i.e., from normal epithelium, low- and high-grade intraepithelial neoplasia to invasive carcinoma. However, how the precancerous lesions progress to carcinoma remains elusive. Here, we report a comprehensive single-cell RNA sequencing and spatial transcriptomic study of 79 multistage esophageal lesions from 29 patients with ESCC. We reveal a gradual and significant loss of ANXA1 expression in epithelial cells due to its transcription factor KLF4 suppression along the lesion progression. We demonstrate that ANXA1 is a ligand to formyl peptide receptor type 2 (FPR2) on fibroblasts that maintain fibroblast homeostasis. Loss of ANXA1 leads to uncontrolled transformation of normal fibroblasts into cancer-associated fibroblasts (CAFs), which can be enhanced by secreted TGF-ß from malignant epithelial cells. Given the role of CAFs in cancer, our study underscores ANXA1/FPR2 signaling as an important crosstalk mechanism between epithelial cells and fibroblasts in promoting ESCC.

LncReader: identification of dual functional long noncoding RNAs using a multi-head self-attention mechanism.

Long noncoding ribonucleic acids (RNAs; LncRNAs) endowed with both protein-coding and noncoding functions are referred to as 'dual functional lncRNAs'. Recently, dual functional lncRNAs have been intensively studied and identified as involved in various fundamental cellular processes. However, apart from time-consuming and cell-type-specific experiments, there is virtually no in silico method for predicting the identity of dual functional lncRNAs. Here, we developed a deep-learning model with a multi-head self-attention mechanism, LncReader, to identify dual functional lncRNAs. Our data demonstrated that LncReader showed multiple advantages compared to various classical machine learning methods using benchmark datasets from our previously reported cncRNAdb project. Moreover, to obtain independent in-house datasets for robust testing, mass spectrometry proteomics combined with RNA-seq and Ribo-seq were applied in four leukaemia cell lines, which further confirmed that LncReader achieved the best performance compared to other tools. Therefore, LncReader provides an accurate and practical tool that enables fast dual functional lncRNA identification.

Enhanced degradation of poly(ethylene terephthalate) by the addition of lactic acid / glycolic acid: composting degradation, seawater degradation behavior and comparison of degradation mechanism.

The degradability improvement of poly(ethylene terephthalate) (PET), one of the most widely used but non-degradable disposable packaging material, is of great significance. However, the balance between degradability and mechanical properties remains a huge challenge. Herein, simple hydroxy acids, lactic acid (LA) and glycolic acid (GA) as easy hydrolysis sites were introduced into non-degradable PET via melt polycondensation. A series of high molecular weight poly(ethylene terephthalate-co-L­lactide) (PETL) and poly(ethylene terephthalate-co-glycolate) (PETG) copolyesters were synthesized with an excellent tensile strength greater than 50 MPa, much higher than that of most commercially available degradable polymers. The introduction of hydroxy acid endows PET with significantly improved composting and seawater degradation performance. Furtherly, the degradation rate of PETG with hydrophilic GA unit was faster than that of PETL, and the mineralization rate of PETG80 reaches 22.0%. The density of functional theory (DFT) calculation revealed that adding hydroxy acid to the PET molecular chain reduced the energy barrier of the hydrolysis reaction. The molecular polarity index (MPI) analysis furtherly confirmed that the higher affinity between the GA unit and water may be the primary reason for the faster degradation of PETG.

Apigenin Ameliorates Hyperuricemia and Renal Injury through Regulation of Uric Acid Metabolism and JAK2/STAT3 Signaling Pathway.

Hyperuricemia (HUA) is a kind of metabolic disease with high incidence that still needs new countermeasures. Apigenin has uric-lowering and kidney-protective activities, but how apigenin attenuates HUA and renal injury remains largely unexploited. To this end, an acute HUA mouse model was established by intraperitoneal injection of potassium oxazinate and oral administration with hypoxanthine for 7 consecutive days. Apigenin intervention decreased serum uric acid (UA), creatinine (CRE), blood urea nitrogen (BUN), interleukin-1ß (IL-1ß), interleukin-6 (IL-6), tumor necrosis factor (TNF-α), interleukin-18 (IL-18), liver xanthine oxidase (XOD), and urine protein levels, and increased serum interleukin-10 (IL-10) and urine UA and CRE levels in HUA mice. Moreover, administration of apigenin to HUA mice prevented renal injury, decreased renal glucose transporter 9 (GLUT9) and urate anion transporter 1 (URAT1) levels, and increased renal organic anion transporter 1 (OAT1). These alterations were associated with an inhibition of IL-6, phospho-janus kinase 2 (P-JAK2), phospho-signal transducer, and activator of transcription 3 (P-STAT3), and suppression of cytokine signaling 3 (SOCS3) expression in the kidneys. Additionally, the molecular docking results showed that apigenin had strong binding capacity with UA transporters and JAK2 proteins. In summary, apigenin could improve UA metabolism and attenuate renal injury through inhibiting UA production, promoting excretion, and suppressing the JAK2/STAT3 signaling pathway in HUA mice. The results suggest that apigenin may be a suitable drug candidate for management of HUA and its associated renal injury.

Characterization of genomic alterations and neoantigens and analysis of immune infiltration identified therapeutic and prognostic biomarkers in adenocarcinoma at the gastroesophageal junction.

Objectives: Adenocarcinoma at the gastroesophageal junction (ACGEJ) refers to a malignant tumor that occurs at the esophagogastric junction. Despite some progress in targeted therapies for HER2, FGFR2, EGFR, MET, Claudin 18.2 and immune checkpoints in ACGEJ tumors, the 5-year survival rate of patients remains poor. Thus, it is urgent to explore genomic alterations and neoantigen characteristics of tumors and identify CD8+ T-cell infiltration-associated genes to find potential therapeutic targets and develop a risk model to predict ACGEJ patients' overall survival (OS). Methods: Whole-exome sequencing (WES) was performed on 55 paired samples from Chinese ACGEJ patients. Somatic mutations and copy number variations were detected by Strelka2 and FACETS, respectively. SigProfiler and SciClone were employed to decipher the mutation signature and clonal structure of each sample, respectively. Neoantigens were predicted using the MuPeXI pipeline. RNA sequencing (RNA-seq) data of ACGEJ samples from our previous studies and The Cancer Genome Atlas (TCGA) were used to identify genes significantly associated with CD8+ T-cell infiltration by weighted gene coexpression network analysis (WGCNA). To construct a risk model, we conducted LASSO and univariate and multivariate Cox regression analyses. Results: Recurrent MAP2K7, RNF43 and RHOA mutations were found in ACGEJ tumors. The COSMIC signature SBS17 was associated with ACGEJ progression. CCNE1 and VEGFA were identified as putative CNV driver genes. PI3KCA and TP53 mutations conferred selective advantages to cancer cells. The Chinese ACGEJ patient neoantigen landscape was revealed for the first time, and 58 potential neoantigens common to TSNAdb and IEDB were identified. Compared with Siewert type II samples, Siewert type III samples had significant enrichment of the SBS17 signature, a lower TNFRSF14 copy number, a higher proportion of samples with complex clonal architecture and a higher neoantigen load. We identified 10 important CD8+ T-cell infiltration-related Hub genes (CCL5, CD2, CST7, GVINP1, GZMK, IL2RB, IKZF3, PLA2G2D, P2RY10 and ZAP70) as potential therapeutic targets from the RNA-seq data. Seven CD8+ T-cell infiltration-related genes (ADAM28, ASPH, CAMK2N1, F2R, STAP1, TP53INP2, ZC3H3) were selected to construct a prognostic model. Patients classified as high risk based on this model had significantly worse OS than low-risk patients, which was replicated in the TCGA-ACGEJ cohort. Conclusions: This study provides new neoantigen-based immunotherapeutic targets for ACGEJ treatment and effective disease prognosis biomarkers.

Resolving mutational signatures in cancer development.

In a recent Science study, an advanced mutational signature analysis was performed on an unprecedentedly large set of cancer whole-genome sequences across multiple organs. This work introduces a computational framework to identify common and rare mutational signatures and provides insights into their biological interpretation.

Computational Identification of Preneoplastic Cells Displaying High Stemness and Risk of Cancer Progression.

Evidence points toward the differentiation state of cells as a marker of cancer risk and progression. Measuring the differentiation state of single cells in a preneoplastic population could thus enable novel strategies for early detection and risk prediction. Recent maps of somatic mutagenesis in normal tissues from young healthy individuals have revealed cancer driver mutations, indicating that these do not correlate well with differentiation state and that other molecular events also contribute to cancer development. We hypothesized that the differentiation state of single cells can be measured by estimating the regulatory activity of the transcription factors (TF) that control differentiation within that cell lineage. To this end, we present a novel computational method called CancerStemID that estimates a stemness index of cells from single-cell RNA sequencing data. CancerStemID is validated in two human esophageal squamous cell carcinoma (ESCC) cohorts, demonstrating how it can identify undifferentiated preneoplastic cells whose transcriptomic state is overrepresented in invasive cancer. Spatial transcriptomics and whole-genome bisulfite sequencing demonstrated that differentiation activity of tissue-specific TFs was decreased in cancer cells compared with the basal cell-of-origin layer and established that differentiation state correlated with differential DNA methylation at the promoters of these TFs, independently of underlying NOTCH1 and TP53 mutations. The findings were replicated in a mouse model of ESCC development, and the broad applicability of CancerStemID to other cancer-types was demonstrated. In summary, these data support an epigenetic stem-cell model of oncogenesis and highlight a novel computational strategy to identify stem-like preneoplastic cells that undergo positive selection. SIGNIFICANCE: This study develops a computational strategy to dissect the heterogeneity of differentiation states within a preneoplastic cell population, allowing identification of stem-like cells that may drive cancer progression.

Multi-omic characterization of genome-wide abnormal DNA methylation reveals diagnostic and prognostic markers for esophageal squamous-cell carcinoma.

This study investigates aberrant DNA methylations as potential diagnosis and prognosis markers for esophageal squamous-cell carcinoma (ESCC), which if diagnosed at advanced stages has <30% five-year survival rate. Comparing genome-wide methylation sites of 91 ESCC and matched adjacent normal tissues, we identified 35,577 differentially methylated CpG sites (DMCs) and characterized their distribution patterns. Integrating whole-genome DNA and RNA-sequencing data of the same samples, we found multiple dysregulated transcription factors and ESCC-specific genomic correlates of identified DMCs. Using featured DMCs, we developed a 12-marker diagnostic panel with high accuracy in our dataset and the TCGA ESCC dataset, and a 4-marker prognostic panel distinguishing high-risk patients. In-vitro experiments validated the functions of 4 marker host genes. Together these results provide additional evidence for the important roles of aberrant DNA methylations in ESCC development and progression. Our DMC-based diagnostic and prognostic panels have potential values for clinical care of ESCC, laying foundations for developing targeted methylation assays for future non-invasive cancer detection methods.

CellCall: integrating paired ligand-receptor and transcription factor activities for cell-cell communication.

With the dramatic development of single-cell RNA sequencing (scRNA-seq) technologies, the systematic decoding of cell-cell communication has received great research interest. To date, several in-silico methods have been developed, but most of them lack the ability to predict the communication pathways connecting the insides and outsides of cells. Here, we developed CellCall, a toolkit to infer inter- and intracellular communication pathways by integrating paired ligand-receptor and transcription factor (TF) activity. Moreover, CellCall uses an embedded pathway activity analysis method to identify the significantly activated pathways involved in intercellular crosstalk between certain cell types. Additionally, CellCall offers a rich suite of visualization options (Circos plot, Sankey plot, bubble plot, ridge plot, etc.) to present the analysis results. Case studies on scRNA-seq datasets of human testicular cells and the tumor immune microenvironment demonstrated the reliable and unique functionality of CellCall in intercellular communication analysis and internal TF activity exploration, which were further validated experimentally. Comparative analysis of CellCall and other tools indicated that CellCall was more accurate and offered more functions. In summary, CellCall provides a sophisticated and practical tool allowing researchers to decipher intercellular communication and related internal regulatory signals based on scRNA-seq data. CellCall is freely available at https://github.com/ShellyCoder/cellcall.

Single-cell transcriptomic analysis in a mouse model deciphers cell transition states in the multistep development of esophageal cancer.

Esophageal squamous cell carcinoma (ESCC) is prevalent in some geographical regions of the world. ESCC development presents a multistep pathogenic process from inflammation to invasive cancer; however, what is critical in these processes and how they evolve is largely unknown, obstructing early diagnosis and effective treatment. Here, we create a mouse model mimicking human ESCC development and construct a single-cell ESCC developmental atlas. We identify a set of key transitional signatures associated with oncogenic evolution of epithelial cells and depict the landmark dynamic tumorigenic trajectories. An early downregulation of CD8+ response against the initial tissue damage accompanied by the transition of immune response from type 1 to type 3 results in accumulation and activation of macrophages and neutrophils, which may create a chronic inflammatory environment that promotes carcinogen-transformed epithelial cell survival and proliferation. These findings shed light on how ESCC is initiated and developed.