Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “single cell transcriptomics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20Linked to original sources

CoxFormer enables spatial omics inference with multimodal generative modeling.

Gene co-expression maps transcriptome-wide gene-gene relationships, yet high-quality estimates cover less than half the genome. Meanwhile, spatial omics either profiles restricted in situ panels or lacks cellular resolution. Extending co-expression transcriptome-wide could overcome these limitations by inferring unassayed gene expression at subcellular resolution. Here we show that CoxFormer integrates literature-derived gene knowledge with co-expression networks from bulk tissues and large-scale single-cell atlases to learn 512-dimensional representations for 32,016 human genes. These embeddings capture functional gene relationships and serve as a generative prior for spatial inference across platforms and modalities. Without requiring a matched single-cell RNA-sequencing reference, CoxFormer supports four applications beyond measured genes: histology-based expression imputation, gene activity prediction from chromatin accessibility, subcellular super-resolution inference, and pathological region detection. Together, CoxFormer extends gene embedding from gene- and cell-level tasks to whole-transcriptome spatial inference, providing a unified framework for biological analysis beyond the limited gene coverage of current spatial omics technologies.

Humans↗

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma.

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities with adult glioma, it comprises distinct disease entities. In this study, we longitudinally profile a molecularly diverse cohort of 16 pHGG patients through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of neoplastic and microenvironmental features during disease progression and treatment. We define a set of core pHGG neoplastic cell states and observe differential tumor-myeloid interactions between malignant cell phenotypes. We find that essential neuromodulators and the interferon response are upregulated post-therapy, implicating them as malignant cell-intrinsic targets. We observe an increase in oligodendrocytes upon progression and that they coordinate spatial motifs with proneural tumor cells. This multiomic atlas of longitudinal pHGG captures features of therapy response and provides a scalable reference for the study of pediatric brain tumors.

Humans↗

Divergent PTEN-p53 interaction upon DNA damage in a human thyroid organoid model with germline PTEN mutations.

Germline mutations in the tumor suppressor phosphatase and tensin homolog (PTEN) cause PTEN hamartoma tumor syndrome (PHTS). PHTS is characterized by an elevated lifetime risk of differentiated thyroid cancer (DTC), 30 times higher than the general population. However, only 1 in 3 PHTS patients develop DTC, and it remains unknown whether specific PTEN variants are associated with an increased risk of DTC. PTEN antagonizes the phosphatidylinositol 3-kinase (PI3K)-AKT signaling pathway, a frequently affected pathway in sporadic DTC. PTEN also acts as a guardian of the genome by interacting with other tumor suppressors. Here, we report how ionizing radiation, an environmental tumorigenic contributor, modifies the DNA damage response based on the type of germline PTEN variants. We hypothesized that certain PTEN variants associated with DTC create a pro-oncogenic molecular signature upon radiation-induced DNA damage. DTC-associated (PTEN M134R ) or DTC-non-associated (PTEN G132D ) germline PTEN mutant alleles were introduced into a human induced pluripotent cell (hiPSC) line derived from a healthy donor utilizing CRISPR-Cas9 gene editing technology. We determined radiation-induced transcriptomic changes in functional thyroid organoids induced from wild-type and both heterozygous PTEN mutant hiPSCs. Both bulk and single-cell RNA sequencing data indicated that radiation upregulated the p53 network more potently in the thyroid organoids with PTEN WT/G132D than those with PTEN WT/M134R , which could be mediated by AKT-dependent MDM2 inactivation and PTEN-p53 physical interaction. Our data suggest that the lack of p53 pathway activation through PTEN-p53 network interactions explains why PTEN M134R is a DTC-susceptible variant.

Humans↗

Proteomic Characterization of 1000 Human and Murine Neutrophils Freshly Isolated From Blood and Sites of Sterile Inflammation.

Neutrophils are indispensable for defense against pathogens. Injured tissue-infiltrated neutrophils can establish a niche of chronic inflammation and promote degeneration. Studies investigated transcriptome of single-infiltrated neutrophils which could misinterpret molecular states of these post mitotic cells. However, neutrophil proteome characterization has been challenging due to low harvests from affected tissues. Here, we present a workflow to obtain proteome of 1000 murine and human tissue-infiltrated neutrophils. We generated spectral libraries containing ∼6200 mouse and ∼5300 human proteins from circulating neutrophils. 4800 mouse and 3400 human proteins were recovered from 1000 cells with 102-108 copies/cell. Neutrophils from stroke-affected mouse brains adapted to the glucose-deprived environment with increased mitochondrial activity and ROS-production, while cells invading inflamed human oral cavities increased phagocytosis and granule release. We provide an extensive protein repository for resting human and mouse neutrophils, identify proteins lost in low input samples, thus enabling the proteomic characterization of limited tissue-infiltrated neutrophils.

Proteomics↗

Through the lens of bioenergy crops: advances, bottlenecks, and promises of plant engineering.

Advances in engineering of bioenergy crops were driven over the past years by adapting technological breakthroughs and accelerating conventional applications but also exposed intriguing challenges. New tools revealed rich interconnectivity in the exponentially growing and dynamic 'big' omics data' of metabolomes, transcriptomes, and genomes at previously inaccessible magnitude (global, cross-species, meta-) and resolution (single cell). Insights enabled fresh hypotheses and stimulated disciplines such as functional genomics with discovery of broad regulatory networks and their determinants, that is, DNA parts, including promoters, regulatory elements, and transcription factors. Their rational design, assembly into increasingly complex blueprints, and installation into diverse chassis is an existing frontier that may benefit from emerging technologies to address bottlenecks. Interweaving nature-inspired to fully synthetic parts has already allowed building of fine-tuned regulatory circuits, or new-to-nature metabolic routes insulated from the biological context of the chassis species. Similarly, developments and the evolving need for unifying principles in plant transformation and species-agnostic technologies highlight future opportunities for engineering the next generation of bioenergy plants.

Crops, Agricultural↗

2025 Donald Seldin Lecture: Leveraging Diverse Population Genomics and Multiomics Integration for Gene Discovery of Cardiovascular and Kidney Diseases.

This review discusses the implications of frameworks leveraging genetic admixture and multiomics data for advancing gene discovery in cardiovascular and kidney disease research. By broadening gene discovery efforts to additional populations that have a disproportionately high risk of disease and leveraging genetic diversity in admixed populations, studies can identify population-enriched risk variants that traditionally have been missed in genome-wide association studies. The use of multiomics approaches, including the transcriptome, proteome, and metabolome, advances a mechanistic understanding of disease beyond associations. As single-cell omics technologies continue to improve, their integration into gene discovery may help uncover cell-type-specific regulatory pathways and more precise biological contexts. The full potential of these approaches depends on sustained investment in diverse, well-characterized omics data sets, methodological innovation in multiancestry statistical approaches, and interdisciplinary collaboration bridging genomics, epidemiology, and clinical medicine. These efforts will need to be translated into clinically actionable insights, including ancestry-informed risk stratification and targeted therapeutics, to improve outcomes for cardiovascular and kidney diseases.

Humans↗

Complement expression profiles in human glomerular mesangial cells, endothelial cells, podocytes and proximal tubular epithelial cells.

BACKGROUND: Local expression of complement components in the kidney has been reported sporadically in both diseased and normal kidneys. This study aimed to comprehensively characterize the expression of complement components in human glomerular mesangial cells (GMCs), glomerular endothelial cells (GECs), podocytes, and proximal tubular epithelial cells (PTECs) in non-diseased renal tissue. METHODS: Complement expression in cultured human renal intrinsic cells was initially evaluated using reverse transcription polymerase chain reaction and immunofluorescence staining. These findings were further examined using publicly available single-cell RNA-sequencing datasets and 10×Genomics single-cell RNA sequencing of non-diseased human kidney tissue. The analyses focused on complement components involved in the initiation of the classical, lectin, and alternative pathways, as well as components shared among these activation pathways, terminal pathway components, complement regulators, and complement receptors. RESULTS: Complement components unique to the initial phase for classical pathway (C1S, C1R, C2, C4), lectin pathway (MBL2, FCN1, MASP1), alternative pathway (CFB, CFD), and the C3 component shared by the three activation pathways were detected in these cells. The components shared by the terminal pathways including C5, C6, C7, C8 and C9 exhibited lower expression, while complement regulators (CFH, CFI, CD55/DAF, CD46/MCP, CD59, C4BPB, PROS1/Protein S) or receptors (CD93/C1QR1, CR1), particularly membrane-bound proteins, such as DAF, MCP and CD59, which inhibit complement activation and the formation of the membrane attack complex, showed relatively high expression. CONCLUSION: These results showed that all four types of intrinsic renal cells expressed multiple complement components associated with the classical, lectin, and alternative pathways. In non-diseased kidney tissue, complement regulatory molecules involved in the control of complement activation showed relatively higher expression, whereas components of the terminal complement pathway were expressed at relatively lower levels, suggesting that renal intrinsic cells maintain a locally poised but tightly regulated complement system.

Humans↗

Paired analysis of primary adenoid cystic carcinoma and derived cell lines reveals a mesenchymal and stem-like shift associated with therapy resistance.

Adenoid cystic carcinoma (ACC) is a salivary gland malignancy characterized by slow but persistent growth, frequent local recurrence, and late metastatic progression. Patients with unresectable, recurrent, or metastatic disease have limited therapeutic options. Efforts to identify effective therapeutic targets have been hindered by the limited availability of well-characterized ACC models. In this study, we established 11 ACC cell lines and performed RNA sequencing of nine cell lines and their matched primary tumors to evaluate the preservation and evolution of molecular and lineage-associated characteristics during cell line establishment. Comparative transcriptomic analysis revealed reduced epithelial and luminal differentiation programs in the cell lines, accompanied by enrichment of myoepithelial, EMT-, and cancer stem cell-associated transcriptional programs. Digital deconvolution and single-sample gene set enrichment analysis supported enrichment of hybrid EMT/stem-like states during in vitro propagation, while comparison with publicly available primary-recurrent ACC data demonstrated partial preservation of recurrence-associated plasticity and invasion programs. Protein-level validation of representative epithelial, myoepithelial, EMT, and stemness markers supported the major transcriptomic changes. In addition, a cell line with a higher stemness signature showed reduced sensitivity to cisplatin. Together, these findings indicate that ACC cell line establishment is associated with transcriptional reprogramming and enrichment of plastic, EMT/stem-like states while retaining selected ACC lineage characteristics. These models provide experimentally tractable platforms for investigating ACC progression, therapeutic response, and mechanisms of treatment resistance.

Adenoid cystic carcinoma↗

A High-Resolution Stereo-Seq Spatial Transcriptomic Resource for Adult Holstein Cattle Liver.

The bovine liver is a highly compartmentalized organ that plays essential roles in continuous gluconeogenesis and nitrogen recycling; however, its spatial molecular architecture has remained largely uncharacterized due to the limitations of traditional bulk and single-cell approaches. To address this gap, Spatial Enhanced Resolution Omics-sequencing (Stereo-seq) was utilized to generate a subcellular-resolution (500 nm) transcriptomic map of an adult Holstein cattle liver, and a refined reference-guided workflow was implemented to overcome standard annotation limitations in livestock. Raw sequencing data were processed using the Stereo-seq Analysis Workflow and analyzed with Stereopy, Seurat, SingleR, and reference-guided workflows. Spatial aggregation was evaluated at Bin20, Bin50, Bin100, Bin150, and Bin200. Increasing bin size increased molecular identifier counts and detected-gene complexity while progressively reducing spatial granularity. Bin50, corresponding to 50 × 50 DNA nanoballs and an approximate nominal footprint of 25 × 25 µm, was therefore selected as a practical intermediate aggregation level for the primary analyses. Quality-control assessment, Leiden clustering, UMAP visualization, reference-based cell-type annotation, cluster-marker analysis, and spatial mapping of canonical hepatic genes demonstrated preservation of biologically interpretable liver transcriptional organization. Raw sequencing data processed spatial matrices, annotated objects, and analysis code are publicly available to support reanalysis and computational benchmarking. In summary, we present a Stereo-seq spatial transcriptomic resource generated from liver tissue of an adult Holstein cow. This initial resource provides a valuable foundation for future studies of bovine liver biology, comparative genomics, and the spatial basis of livestock health and production traits.

Animals↗

An immune exhaustion signature predicts prognosis and identifies patients with diffuse large B-cell lymphoma (DLBCL) who derive preferential benefit from chimeric antigen receptor (CAR)-T cell therapy.

BACKGROUND: The tumor microenvironment (TME) is a key determinant of prognosis in diffuse large B-cell lymphoma (DLBCL). While T-cell exhaustion is implicated in therapeutic failure, its precise molecular hallmarks and utility for predicting response to modern immunotherapies, such as chimeric antigen receptor (CAR)-T cell therapy, remain unclear. METHODS: We performed an integrative analysis of transcriptomic and clinical data from multiple DLBCL cohorts (The Cancer Genome Atlas [TCGA], GSE181063, GSE10846, GSE248835, GSE182434). We used unsupervised clustering, exploratory analysis of single-cell RNA sequencing data, and the least absolute shrinkage and selection operator for variable selection (LASSO-Cox) regression to characterize the exhausted TME, construct a prognostic model, and evaluate its predictive value for CAR-T cell therapy. The model's dynamic behavior was assessed in a proof-of-concept longitudinal cohort of patients treated with the T-cell-engaging bispecific antibody glofitamab. RESULTS: We identified a "high-exhaustion" subtype associated with significantly poorer overall survival (OS; log-rank P = 0.016). Based on this, we developed a five-gene immune exhaustion-Related Prognostic Score (IERPS) that served as a robust independent predictor of poor OS across multiple cohorts. Critically, in a cohort of 256 relapsed/refractory patients, the IERPS was strongly prognostic for event-free survival (EFS) in the standard-of-care (SOC) arm (HR = 2.02, 95% confidence interval [95% CI]: 1.07-3.81, P = 0.029) but lost prognostic significance in the CAR-T arm (HR = 0.70, 95 % CI: 0.35-1.40, P = 0.314). This significant interaction suggests that CAR-T cell therapy may abrogate the poor prognosis associated with a high IERPS. Biologically, exploratory single-cell analysis (n = 4 samples) defined the high-IERPS state by hallmarks of classical T-cell exhaustion, and a descriptive case study showed the score dynamically tracked clinical response to glofitamab. CONCLUSIONS: A state of active T-cell exhaustion and a suppressive TME drive the adverse immune phenotype in DLBCL. Our IERPS model captures this dysfunctional state, acting as a powerful prognostic tool and, more importantly, as a potential predictive biomarker to identify high-risk patients who appear to overcome their inherently poor prognosis through CAR-T cell therapy.

Biomarkers↗

Integrated bulk and single-cell RNA sequencing reveals a prognostic neuro-mimicry signature in papillary thyroid carcinoma.

BACKGROUND: Cancer cells can acquire neuron-like characteristics ("neural mimicry") to promote progression. However, the role of specific ion channel genes in Papillary Thyroid Carcinoma (PTC) and their clinical significance remains unclear. METHODS: We included transcriptomic data from 521 PTC patients in the TCGA cohort. A neuron-specific gene set was used to screen for potential targets. We constructed a prognostic model using LASSO logistic regression. To verify the cellular origin of the signature, we performed single-cell RNA sequencing (scRNA-seq) analysis on the GSE184362 dataset. RESULTS: We established an 8-gene signature involving KCNN4, KCNN1, KCNT2, SNAP25, KCNK16, GABRG1, GABRG2, and GABRB2. The model demonstrated good predictive performance for lymph node metastasis, with an AUC of 0.721 (95% CI 0.677-0.765). Single-cell analysis of seven integrated tumor samples (N = 65,744 cells) confirmed that GABRB2 was specifically enriched in malignant thyrocytes (EPCAM+/KRT18+) at 200-fold higher detection rates than immune cells (20.0% vs. 0.1%, P ≈ 0), supporting tumor-intrinsic neural mimicry. High-risk patients showed immunosuppressive features with altered immune cell infiltration patterns. CONCLUSION: This study identifies a malignant cell-intrinsic signature for predicting PTC prognosis. Validated by single-cell data, our findings suggest that targeting ion channels may represent a potential therapeutic strategy for modulating neuro-immune interactions in thyroid cancer, pending experimental validation.

GABRB2↗

Spatiotemporal single-cell roadmap of human skin wound healing.

Wound healing is vital for human health, yet the details of cellular dynamics and coordination in human wound repair remain largely unexplored. To address this, we conducted single-cell multi-omics analyses on human skin wound tissues through inflammation, proliferation, and remodeling phases of wound repair from the same individuals, monitoring the cellular and molecular dynamics of human skin wound healing at an unprecedented spatiotemporal resolution. This singular roadmap reveals the cellular architecture of the wound margin and identifies FOSL1 as a critical driver of re-epithelialization. It shows that pro-inflammatory macrophages and fibroblasts sequentially support keratinocyte migration like a relay race across different healing stages. Comparison with single-cell data from venous and diabetic foot ulcers uncovers a link between failed keratinocyte migration and impaired inflammatory response in chronic wounds. Additionally, comparing human and mouse acute wound transcriptomes underscores the indispensable value of this roadmap in bridging basic research with clinical innovations.

Humans↗

Spatial Transcriptomics Identifies Characteristic Immunological Niches in Atopic Dermatitis.

BACKGROUND: Atopic dermatitis (AD) is primarily driven by a Type 2 immune response, with T helper (TH2) cells producing IL-4 and IL-13, thereby promoting inflammation, itch, and a compromised skin barrier. Yet, the spatial organization of pathogenic immune cells and their interactions with stromal and epithelial compartments in human AD skin remain incompletely understood. METHODS: We performed 10× Genomics Visium spatial transcriptomics on FFPE skin biopsies from patients with AD (n = 6), psoriasis (n = 2), and healthy controls (n = 5). Data were integrated with AD single-cell RNA sequencing (scRNA-seq) datasets and complemented by imaging mass cytometry (IMC) and multiplex immunofluorescence (IF) to validate the spatial localization of immune cells. Cell-cell communication analysis revealed putative signaling interactions within immune niches. RESULTS: Spatial clustering resolved tissue compartments and demonstrated transcriptional dysregulation in keratinocytes in AD and psoriasis. AD lesions showed a conserved spatial organization of immune aggregates within the superficial dermis. Integration of scRNA-seq signatures revealed spatially organized co-localization of T cells and mature migratory dendritic cells (mmDCs). We developed a ring-based neighborhood analysis to characterize the cellular organization of the immune-stromal niches, revealing T cell-enriched regions surrounded by inflammatory fibroblasts and activated keratinocytes. Intercellular communication analysis further identified putative signaling within mmDC-T cell niches that may promote pathogenic T cell recruitment and activation. Application of tertiary lymphoid structure (TLS) signatures indicated the presence of TLS-like regions. IMC and IF validated the close spatial proximity between activated TH2 cells and mmDCs. CONCLUSION: AD lesions contain spatially organized TLS-like immune niches at the dermal-epidermal junction, characterized by the close association of T cells and mmDCs and coordinated interactions with surrounding stromal and epithelial compartments. These mmDC-T cell niches may represent potential targets for future therapeutic strategies aimed at disrupting persistent local inflammatory pathways and improving long-term disease control.

atopic dermatitis↗

Single-cell profiling of trabecular meshwork identifies mitochondrial dysfunction in a glaucoma model that is protected by vitamin B3 treatment.

Since the trabecular meshwork (TM) is central to intraocular pressure (IOP) regulation and glaucoma, a deeper understanding of its genomic landscape is needed. We present a multimodal, single-cell resolution analysis of mouse limbal cells (includes TM). In total, we sequenced 9,394 wild-type TM cell transcriptomes. We discovered three TM cell subtypes with characteristic signature genes validated by immunofluorescence on tissue sections and whole-mounts. The subtypes are robust, being detected in datasets for two diverse mouse strains and in independent data from two institutions. Results show compartmentalized enrichment of critical pathways in specific TM cell subtypes. Distinctive signatures include increased expression of genes responsible for 1) extracellular matrix structure and metabolism (TM1 subtype), 2) secreted ligand signaling to support Schlemm's canal cells (TM2), and 3) contractile and mitochondrial/metabolic activity (TM3). ATAC-sequencing data identified active transcription factors in TM cells, including LMX1B. Mutations in LMX1B cause high IOP and glaucoma. LMX1B is emerging as a key transcription factor for normal mitochondrial function and its expression is much higher in TM3 cells than other limbal cells. To understand the role of LMX1B in TM function and glaucoma, we single-cell sequenced limbal cells from Lmx1b V265D/+ mutant mice (2,491 TM cells). In V265D/+ mice, TM3 cells were uniquely affected by pronounced mitochondrial pathway changes. Mitochondria in TM cells of V265D/+ mice are swollen with a reduced cristae area, further supporting a role for mitochondrial dysfunction in the initiation of IOP elevation in these mice. Importantly, treatment with vitamin B3 (nicotinamide), to enhance mitochondrial function and metabolic resilience, significantly protected Lmx1b mutant mice from IOP elevation.

Journal Article↗

Lithium deficiency and the onset of Alzheimer's disease.

The earliest molecular changes in Alzheimer's disease (AD) are poorly understood1-5. Here we show that endogenous lithium (Li) is dynamically regulated in the brain and contributes to cognitive preservation during ageing. Of the metals we analysed, Li was the only one that was significantly reduced in the brain in individuals with mild cognitive impairment (MCI), a precursor to AD. Li bioavailability was further reduced in AD by amyloid sequestration. We explored the role of endogenous Li in the brain by depleting it from the diet of wild-type and AD mouse models. Reducing endogenous cortical Li by approximately 50% markedly increased the deposition of amyloid-β and the accumulation of phospho-tau, and led to pro-inflammatory microglial activation, the loss of synapses, axons and myelin, and accelerated cognitive decline. These effects were mediated, at least in part, through activation of the kinase GSK3β. Single-nucleus RNA-seq showed that Li deficiency gives rise to transcriptome changes in multiple brain cell types that overlap with transcriptome changes in AD. Replacement therapy with lithium orotate, which is a Li salt with reduced amyloid binding, prevents pathological changes and memory loss in AD mouse models and ageing wild-type mice. These findings reveal physiological effects of endogenous Li in the brain and indicate that disruption of Li homeostasis may be an early event in the pathogenesis of AD. Li replacement with amyloid-evading salts is a potential approach to the prevention and treatment of AD.

Alzheimer Disease↗

Ovarian development is driven by early spatiotemporal priming of the coelomic epithelium.

Ovarian organogenesis requires the coordinated specification of supporting and steroidogenic cell lineages from multipotent coelomic epithelium (CE) progenitors. A longstanding question is whether the CE contains transcriptionally distinct, spatially organized progenitor subpopulations with predetermined lineage biases, or whether specification into supporting and steroidogenic lineages occurs only after delamination and integration into the bipotential gonad. The developmental origins of granulosa cells and the emergence of ovarian steroidogenic/stromal progenitors (SPs) also remain poorly defined. Here, we show that CE cells covering the fetal mouse ovary are transcriptionally heterogeneous and spatially organized into subdomains already primed toward supporting or steroidogenic fates. CE priming is dynamic, with transient coexistence of supporting- and steroidogenic-biased CE progenitors before resolving into a predominantly supporting-biased CE. Local delamination of these primed cells seeds intragonadal niches where pre-granulosa cells and SPs mirror the spatio-temporal arrangements of CE-primed progenitors. We further demonstrate a dual origin for the supporting lineage, with granulosa cells deriving from both the CE and supporting-like cells (SLCs). In parallel, we show that SPs arise from steroidogenic-primed CE cells, expand to represent 52% of ovarian somatic cells at birth, persist into adulthood and contribute to both theca and steroidogenic stromal cells. Together, these findings reveal transcriptionally and spatially distinct CE subpopulations that shape somatic lineage emergence with important implications for ovarian pathophysiology.

Ovarian development↗

Overexpression of TCF7L2 promotes the viability and migration of MHCC-97H human hepatocellular carcinoma cells by upregulating MT-ND4L.

BACKGROUND: Hepatocellular carcinoma (HCC) is a highly aggressive cancer with high metabolic adaptability. TCF7L2, a transcription factor implicated in type 2 diabetes and cancer, is overexpressed in HCC. However, its specific role in HCC metabolic reprogramming is not well defined. We aimed to elucidate the previously unrecognized molecular mechanisms through which TCF7L2 impacts HCC progression. METHODS: To investigate the function of TCF7L2, a stable MHCC-97H cell line with TCF7L2 overexpression was established via lentiviral transduction. Cell viability and migration were assessed by Cell Counting Kit-8 (CCK-8) and Transwell assays. Transcriptomic profiling [RNA sequencing (RNA-seq)] was performed to identify differentially expressed genes (DEGs). Functional enrichment analysis [Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA)] and bioinformatics promoter analysis (the JASPAR CORE database) were conducted. Clinical correlations, survival analysis, and tumor microenvironment (TME) interrogation were performed using The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort and single-cell datasets [Human Protein Atlas (HPA), CellChat]. Drug sensitivity was predicted via the Genomics of Drug Sensitivity in Cancer (GDSC) database. RESULTS: TCF7L2 overexpression significantly promoted HCC cell proliferation and migration. Transcriptomic analysis revealed that TCF7L2 drives a profound metabolic shift, with key enrichments in lipid homeostasis, fatty acid β-oxidation, and the PI3K/Akt pathway. Mechanistically, TCF7L2 directly binds to the promoter of CPT1A, the rate-limiting enzyme of fatty acid oxidation, and indirectly upregulates the mitochondrial gene MT-ND4Lvia a strong positive correlation with the mitochondrial transcription factor TFAM. In clinical cohorts, TCF7L2 was overexpressed in HCC and its expression correlated positively with MT-ND4L, MKI67, and SNAI1, and served as a predictor of poor overall survival (OS). Furthermore, TCF7L2-high tumors were enriched in hepatic progenitor cell (HPC)-like niches, mediated by enhanced ANGPTL4 signaling. High TCF7L2 expression predicted increased sensitivity to PI3K/mTOR pathway inhibitors. CONCLUSIONS: TCF7L2 acts as a master metabolic regulator in HCC, coordinating lipid catabolism and mitochondrial biogenesis to drive aggressive tumor behavior. It further remodels the TME towards an HPC-like state and predicts sensitivity to metabolic-targeted therapies. These findings identify TCF7L2 as a key prognostic biomarker and a promising therapeutic target.

MHCC-97H hepatocellular carcinoma cells (MHCC-97H ↗

BIWT: a bioinformatics walkthrough for embedding spatial multiomics in agent-based models for virtual cells.

SUMMARY: Whereas transcriptomic and spatial profiling offer static snapshots of tissue structure, mechanistic models use biological rules to predict how tissues evolve. We present the BioInformatics WalkThrough (BIWT) software to directly initialize spatial agent-based models from single-cell and spatial molecular data. We demonstrate how initialization strategies affect tumor-immune dynamics and spatial clustering, positioning BIWT as a software suite to generate data-driven virtual cells representing both experimental and clinical contexts. AVAILABILITY AND IMPLEMENTATION: The BIWT software is available at https://github.com/PhysiCell-Tools/PhysiCell-Studio. The sample dataset for running the BIWT is available at https://zenodo.org/records/16365625. The code and instructions for reproducing the use case example is available at https://github.com/drbergman/BIWT-Paper.

Software↗