PubMed HealthSearch

SEARCH · PubMed Health

Results for “Single cell transcriptomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 181 records · Page 10Linked to original sources

Deciphering the Genetic Underpinnings of Liver Cirrhosis-Heart Failure Comorbidity Through Multi-Omics: CRIM1 as a Key Endothelial Mediator.

The co-occurrence of liver cirrhosis (LC) and heart failure (HF) poses considerable clinical challenges, yet the cellular and molecular determinants of this comorbidity remain poorly characterized. To address this, we developed an integrative multi-omics pipeline encompassing GWAS meta-analysis, gsMap-based spatial transcriptomic projection, GeneEnrich functional annotation, single-cell atlas construction, seismicGWAS and ECLIPSER cell-type scoring, eCAVIAR and fastenloc colocalization, hdWGCNA network inference, scTenifoldKnk in silico gene perturbation, and GCTA-COJO fine-mapping. Quality-controlled meta-analysis yielded 12,347,758 and 9,256,862 variant-level associations for LC and HF, respectively. Spatial projection confirmed preferential enrichment of disease signals within embryonic hepatic and cardiac compartments. Pathway analyses disclosed that LC-linked loci were concentrated in lipid metabolic programs, whereas HF-linked loci implicated mitochondrial bioenergetics and lysosomal degradation. At the cellular level, endothelial cells emerged as the dominant HF-associated population. Convergent evidence from five orthogonal algorithms pinpointed CRIM1 as the sole robustly supported shared gene, selectively enriched in HF endothelial cells; virtual perturbation further identified LCP1 and PTPRC as downstream regulatory nodes. Fine-mapping of the chromosome 2 locus harboring rs12476437 revealed multiple statistically independent signals in the vicinity of CRIM1. Collectively, these findings computationally prioritize the endothelial-CRIM1 axis as a previously unappreciated candidate mechanistic bridge between LC and HF requiring experimental validation.

Humans

GZMK+CD8+ T cells target a specific acinar cell type in Sjögren's disease.

OBJECTIVES: Sjögren's disease (SjD) is a systemic autoimmune disorder characterized by dysfunction of exocrine glands, particularly the salivary and lacrimal glands, with no clear etiology or effective therapy. This study explores the complex interplay of varied cell types in the salivary glands and their role in the pathology of Sjögren's disease. METHODS: Utilizing single-cell and spatial transcriptomics alongside spatial immunophenotyping to analyze human minor salivary glands, we developed a comprehensive understanding of the cellular landscape of non-SjD salivary glands and how that landscape changes in SjD patients. In vitro cellular assays and novel patient-derived primary epithelial cells were co-cultured with autologous T cells to confirm effector states and the delivery and effect of disease-associated granzymes. RESULTS: We identified previously unrecognized heterogeneity among acinar cells, including a PRR4⁺CST3⁺WFDC2⁻ seromucous acinar population that is selectively lost in Sjögren's disease. Expression and organizational changes were linked to clinical features: (i) T cells in the glands of SSA⁺, high-focus score patients showed increased transcriptional signatures of activation, antigen presentation, and apoptosis resistance compared with patients with mild or moderate disease, and (ii) patients with low immune infiltration exhibited distinct epithelial organization. Notably, GZMK⁺CD8⁺ T cells, which accumulate with disease severity, displayed a cytotoxic transcriptional program, degranulated upon stimulation ex vivo, and localized spatially with immune-engaged epithelial cells. Functional assays demonstrated that GZMK activates interferon signaling in vitro, and autologous co-cultures of patient-derived T cells and epithelial cells validated these findings. CONCLUSIONS: Using single-cell and spatial transcriptomics and proteomics, this study identifies a selective loss of PRR4⁺CST3⁺WFDC2⁻ seromucous acinar cells and a rise in GZMK⁺CD8⁺ T cells in Sjögren's disease, revealing distinct immune-mediated epithelial remodeling and interferon-driven dysfunction across diverse clinical presentations. These findings uncover a novel sub-cytolytic effector mechanism by which GZMK⁺CD8⁺ T cells impair mitochondrial integrity and activate innate immune signaling, linking epithelial injury to type I interferon responses and offering new therapeutic targets.

Humans

Combined somatic mutation and transcriptome analysis reveals region-specific differences in clonal architecture in human cortex.

The human cerebral cortex is specialized into regions, but little is known about how human cellular lineages shape cortical regional variation and neuronal cell-type distribution during development. Here, we map single-cell lineages of human cortical regions and neuronal subtypes using >1,000 somatic single-nucleotide variants (sSNVs) identified from deep bulk whole-genome sequencing and analyzed over 25 regions and >72,000 single cells. In the fronto-parietal cortex, sSNVs are rarely restricted, marking neuron-generating clones that disperse into neighboring regions. In contrast, the primary visual cortex harbors 30%-70% more sSNVs than the neighboring secondary visual cortex. Clones at this border exhibit more restricted dispersion, suggesting late developmental lineage segregation. Single-nucleus sSNV and whole-transcriptome analysis reveal glutamatergic neuron clones with modest regional restrictions that share low-mosaic sSNVs with some GABAergic neurons, suggesting a recent dorsal cortical progenitor. Our analysis reveals human-specific cortical lineage patterns, regional differences in clonal patterns, and late divergence of some glutamatergic/GABAergic lineages.

Humans

Integrated multi-omics analyses identify an RAS-SLC11A2-associated molecular framework linking iron metabolism with PCOS-related cardiometabolic risk.

INTRODUCTION: PCOS is a common endocrine disorder with elevated cardiometabolic risk, yet the role of the renin-angiotensin system (RAS)-iron metabolism axis in this comorbidity remains unclear. We explored its underlying mechanisms and evaluated the therapeutic potential of gentiopicroside. METHODS: Integrated multi-omics analyses combining transcriptomics, single-cell RNA sequencing, Mendelian randomization, machine learning, molecular docking, and in vitro functional assays were performed to identify shared molecular pathways and therapeutic targets across PCOS, hypertension, NAFLD, and T2DM. RESULTS: SLC11A2 was consistently dysregulated in PCOS transcriptomic datasets, and associated with iron metabolism, inflammatory response and oxidative stress pathways. Genetic analyses validated RAS-related regulation in hypertension susceptibility and revealed shared genetic architecture between PCOS and cardiometabolic traits. Network and single-cell analyses characterized SLC11A2-associated molecular patterns in disease-relevant cell types; machine learning identified disease-classifying molecular signatures. Gentiopicroside alleviated inflammatory and oxidative stress phenotypes, including reduced IL-6 expression and reactive oxygen species accumulation. CONCLUSION: This study defines an RAS-SLC11A2 molecular framework linking iron metabolism dysregulation to PCOS-related cardiometabolic risk, elucidating the mechanisms connecting ovarian dysfunction, inflammation, oxidative stress and hypertension, and supports gentiopicroside as a promising therapeutic candidate.

Humans

Spatiotemporal mapping of tertiary lymphoid structure heterogeneity shapes immune niches and clinical outcomes in intrahepatic cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (iCCA) is a highly lethal malignancy with limited therapeutic options. The spatial architecture and functional diversity of tertiary lymphoid structures (TLSs) in iCCA remain unclear. Here, we present a multimodal spatial atlas of TLSs and identified intratumoral TLSs (iTLSs) as independent prognostic markers. Bulk proteomic profiling of 214 discovery and 155 validation cases identified a four-tier TLS-based tumor microenvironment classification system and supported development of a TLS-predictive random forest classifier. Imaging mass cytometry revealed that iTLS+ tumors harbor structured immune architectures, where M1-like tissue-resident macrophages (RTMs), dendritic cells, and CXCL13+ CD4+ T cells colocalize to form antigen-presenting neighborhoods (apc-CNs) spatially coupled to TLS core regions (TLScore-CNs). Single-cell spatial transcriptomics further resolved 61 TLSs into 14 spatial niches and defined a pseudotemporal maturation continuum: aggregated, activated, and postactivated. Intraniche communication, primarily mediated by ifnCAFs, iCAFs, and CXCL12+ macrophages, evolved dynamically with maturation. Single-nucleus RNA sequencing combined with Tangram-based spatial mapping revealed CXCL12+ macrophages and iCAFs forming a peripheral band in aggregated TLSs, whereas ifnCAFs infiltrated TLS interiors during activation. These findings define TLS heterogeneity and provide insights for stroma-directed immunotherapy.

Cholangiocarcinoma

Cell-type-specific genetic associations in Lewy body dementia identified using single-cell eQTL-based Mendelian randomization.

BACKGROUND: Lewy body dementia (LBD) is a complex neurodegenerative disorder marked by α-synuclein aggregation and dual impairment of cognitive and motor function.While genome-wide association studies have identified risk loci, the cellular mechanisms linking genetic variation to disease susceptibility remain largely unexplored. METHODS: We performed single-cell transcriptome-wide Mendelian randomization using brain cell-type-specific eQTLs across eight major cell types. Genetic associations were evaluated using inverse-variance weighted models, followed by Bayesian colocalization analysis. Replication was performed in independent stratified LBD cohorts based on APOE ε4 carrier status. Phenome-wide association analysis was included as a supplementary, descriptive assessment of cross-trait associations. RESULTS: Expression of ANKRD65 in excitatory neurons was significantly associated with reduced LBD risk (odds ratio = 0.65, 95 % CI: 0.52-0.81, p = 0.00013). This association passed a false discovery rate of 0.1 and showed strong evidence of colocalization (posterior probability = 0.93). Effect direction was consistent across APOE ε4+ and ε4- LBD subgroups in independent cohorts. No genome-wide significant associations were observed with non-neurological traits in the phenome-wide analysis. CONCLUSIONS: Our findings identify a genetically supported, cell-type-resolved association between ANKRD65 expression in excitatory neurons and LBD risk. This study demonstrates the value of integrating cell-resolved transcriptomic regulation with genetic inference to pinpoint functionally relevant targets in neurodegenerative diseases.

Humans

Multi-omics analyses reveal DjTcf4 critical for proper timing of differentiation in planarian regeneration.

The blastema is key to forming complete tissues in regenerating Dugesia japonica (D. japonica). However, the dynamic changes in cellular compositions and transcription landscapes in blastema during regeneration are understudied. Here, through genome reannotation, 3D spatial transcriptome construction, single-cell RNA sequencing (scRNA-seq), and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) analyses of changes in gene expression and chromatin structures, we delineate key transcription factors regulating the developmental trajectories of major cell clusters in the regenerating head. Importantly, we find that the T cell factor 4 (DjTcf4)-positive cells highly accumulate at wound areas, and its gene network is critical for the proper timing of development during regeneration in multiple progenitor cells. Depletion of DjTcf4 and its target genes leads to singular eye and/or dull tail phenotypes and delays regeneration. Taken together, we build multi-omics atlases in D. japonica and reveal the noncanonical function of the DjTcf4 network in developmental pattern formation, laying a foundation for studies of regeneration in D. japonica.

Animals

Spatial biology reveals altered macrophage states in immunosuppressed non-melanoma skin cancer.

Immunosuppressed patients with non-melanoma skin cancer experience worse clinical outcomes, yet the tumor immune microenvironment associated with systemic immunosuppression remains incompletely defined. Using integrated single-cell, spatial transcriptomic, multiplex immunofluorescence, and spatial epigenomic profiling across immunocompetent and immunosuppressed tumors, we found that overall immune-cell composition was largely preserved despite differences in immune-cell distribution, spatial organization, and T cell clonality. Immunosuppressed tumors demonstrated reduced intratumoral macrophage densities, decreased T cell clonal diversity, altered antigen-presenting cell and T cell spatial interactions, and distinct fibroblast- and macrophage-associated spatial niches. Multi-cohort validation across complementary spatial and single-cell platforms identified consistent alterations in innate-adaptive immune organization in immunosuppressed tumors. Together, these findings define spatial and functional remodeling of the tumor immune microenvironment under systemic immunosuppression and provide a framework for future therapeutic investigation in high-risk patients.

Humans

Hypoxia Response Is Associated with Reduced HPV Activity and Tumor Microenvironment Remodeling in Cervical Cancer.

Human papillomavirus (HPV) significantly influences cervical cancer progression and treatment, yet its interactions with the tumor microenvironment remain incompletely understood. We performed single-cell and spatial transcriptomic sequencing on cervical cancer samples to explore these interactions. By aligning sequencing reads to a merged HPV16-human reference genome, we characterized HPV16 heterogeneity and its association with host states at the single-cell and single-gene levels. E5 transcriptional activity was negatively associated with the host interferon response, indicating a role in immune evasion. A hypoxic environment was correlated with the downregulation of E5 activity and elevated MHC-I expression, which may contribute to stronger interactions between hypoxic cancer cells and cytotoxic CD8⁺ T cells. Additionally, HPV16 integration in host cells was associated with increased fatty acid metabolism. These findings suggest that combining anti-angiogenic drugs and fatty acid metabolism inhibitors has the potential to improve cervical cancer treatment.

Cervical cancer

The protective role of γδ T cells in endometrial cancer.

Γδ T cells are non-conventional T cells that are not MHC restricted and have T cell receptors (TCRs) that are stimulated by phosphoantigens, stress-induced proteins, lipids, and other antigens. These cells are prognostic across cancer types in The Cancer Genome Atlas (TCGA) but have not been well studied in endometrial cancer, which has a rising incidence and mortality rate. Endometrial cancer patients have variable responses to checkpoint inhibitors which are related to the molecular subtype of their cancer. As such, there is a pressing need to understand the immune microenvironment in endometrial cancer. This study addresses this gap in knowledge by investigating γδ T cell repertoires and transcriptomes in this disease site. γδ T cell repertoires were obtained for 543 endometrial cancer patients within the TCGA and from 5 endometrial cancer patients in the single cell dataset SRP349751 using TRUST4. GLIPH2 was used to identify TCRs predicted to bind the same antigen. Transcriptomes were investigated in the single cell dataset. DNA Polymerase Epsilon Exonuclease (POLE) and Microsatellite Instability High (MSI-H) endometrial cancer subtypes had the most γδ T cell infiltration. Vδ1 and Vδ3 γδ T cell infiltration was prognostic independent of stage and molecular subtype. GLIPH2 analysis revealed TCRδ motifs for TDK, YTD, and GEL were public across all four molecular subtypes and were present in the single cell data set. Vδ1 γδ T cell transcriptomes were associated with cytotoxicity and recent TCR stimulation. These data support further investigation of immunotherapies targeting γδ T cells in endometrial cancer.

Humans

Human and mouse adrenal glands are characterized by species-specific steroidogenic states and tissue turnover.

The adult adrenal cortex undergoes constant renewal, yet underlying human-specific mechanisms remain poorly understood. Here we generated single-cell and spatial transcriptomic atlases of adult human and mouse adrenal glands, leveraging single-cell-resolution spatial data and a rare clonal mosaic case for lineage inference. In humans, we identified age-associated zona glomerulosa (ZG) cell states with direct cortisol synthesis capacity and sex-specific differences in inferred cholesterol balance. Cross-species comparison revealed conserved aldosterone-producing ZG but notable divergence in zona fasciculata markers, absence of zona reticularis homologs in mice and differential SHH-WNT4 signaling in proliferating cells. We uncovered human WT1- capsule-to-ZG transition and vascular smooth muscle cell-to-steroidogenic transitions supported by mosaic lineage evidence. We revealed dispersed proliferating cortical SF1+EZH2+ cells throughout the human cortex in contrast with ZG restriction in mice. Taken together, our data expand the centripetal renewal model and establish a comparative framework for human adrenocortical biology.

Animals

Mapping the Immune cell-specific gene regulatory network in bipolar disorder: A framework from scTWMR to exploratory drug-target annotation.

BACKGROUND: Although the involvement of the immune system in the genetic susceptibility of bipolar disorder (BD) is widely acknowledged, the causal relationship between gene expression in specific immune cell subtypes and BD requires systematic elucidation. METHODS: We implemented an analytical framework integrating single-cell transcriptome-wide Mendelian randomization (scTWMR) with colocalization analysis. This approach utilized cis-expression quantitative trait loci (cis-eQTLs) derived from 14 distinct immune cell types as instrumental variables to interrogate BD genome-wide association study (GWAS) summary statistics (comprising 41,917 cases and 371,549 controls). Subsequent investigations encompassed functional enrichment analysis, protein-protein interaction (PPI) network construction, phenome-wide association study (PheWAS), and performed an exploratory drug-target annotation. RESULTS: Our analysis identified 33 gene-immune cell associations. Colocalization analysis provided robust evidence (PPH4 > 90%) for shared causal variants implicating the MAD1L1, APOM, and NFKBIL1 loci. Significantly enriched biological pathways included cell cycle regulation, circadian rhythm entrainment, and neuroinflammation. The PPI network revealed a core regulatory module centered on histone-encoding and immune-related genes. Exploratory drug-target annotation nominated compounds for further investigation for compounds targeting APOM, TMEM258, and NFKBIL1. CONCLUSION: This study systematically delineates a genetically supported regulatory network of immune cell-specific gene expression in BD, predominantly implicating CD8⁺ effector T cells, plasma cells, and B cells. The findings corroborate established pathological pathways while uncovering novel cell type-specific therapeutic targets, thereby providing a genetic framework for prioritizing candidate targets for future investigation.

Bipolar disorder

Pan-genomics and multi-omics for deciphering genetic variation and accelerating genetic improvement in ruminant livestock.

Livestock reference genomes have transformed the discovery of variants associated with production, reproduction, health, and environmental adaptation. Nevertheless, a single linear reference represents only one mosaic haplotype and incompletely captures sequence diversity within a species, particularly structural variants, copy-number changes, repeat-rich regions, and breed-specific sequences. Pangenomes address this limitation by integrating multiple high-quality assemblies or population-scale variants into a unified sequence or graph representation. Concurrently, multi-omics approaches connect genomic variation with transcriptomic, epigenomic, manuscriptproteomic, metabolomic, and microbiome responses, thereby improving biological interpretation of genotype-phenotype relationships. This review synthesizes recent progress in livestock pangenomics and multi-omics, with emphasis on cattle, goats, sheep, water buffalo, and chickens. It describes advances in long-read and haplotype-resolved sequencing, graph construction, structural-variant discovery and genotyping, functional annotation, and integrative analysis. Recent pangenome studies have uncovered substantial non-reference sequence, reduced reference bias, identified breed- and population-specific structural variants, and resolved candidate variants underlying pigmentation, body size, tail morphology, cashmere production, altitude adaptation, and other economically relevant traits. However, translation into routine breeding remains constrained by uneven population representation, inconsistent structural-variant definitions, limited functional annotation, computational demands, and insufficient validation across environments. Future progress will depend on diverse near-complete assemblies, graph-aware imputation and genomic prediction, long-read transcriptomics, single-cell and spatial omics, rigorous causal validation, and open, interoperable resources. Together, these developments can support more accurate, resilient, and biologically informed livestock improvement. Importantly, current dairy-cattle evidence indicates that pangenome-derived structural variants can substantially improve variant discovery and functional interpretation while yielding only marginal average gains in routine genomic prediction, favoring targeted augmentation rather than wholesale replacement of established SNP-based evaluations.

Animals

In vivo differentiation of embryonic cells devoid of key reprogramming factors.

Embryonic cell differentiation depends on reprogramming of the oocyte and sperm nucleus into a transient totipotent state. In zebrafish, this coincides with genome activation, which is regulated by the pioneer factors Nanog, Pou5f3, and Sox19b (NPS). Here, we investigate the role of NPS in developmental reprogramming and differentiation by analyzing the fate of NPS mutant cells in a wild-type embryo using single-cell RNA-seq. We find that many cells fail to activate transcription or undergo cell death, while others acquire gene expression profiles that resemble germ cells, neural progenitors, and motoneuron states. These cells achieve intermediate transcriptional states, revealing the essential role of NPS in coordinating nuclear and cytoplasmic reprogramming and preventing the premature activation of lineage-specific differentiation programs. These results demonstrate that most developmental programs require developmental reprogramming by NPS, yet some cells can bypass transient totipotency to achieve intermediate developmental states resembling wild-type states in vivo.

Animals

Spatially defined microenvironmental niches are associated with clinical outcome and tumor ecosystem diversity in head and neck cancer.

BACKGROUND: Head and neck squamous cell carcinoma (HNSCC) exhibits substantial biological heterogeneity that is not fully explained by human papillomavirus (HPV) status. The spatial organization of tumor, immune, and stromal cell populations and its relationship to clinical outcome remain incompletely understood. METHODS: We performed single-cell spatial transcriptomic and proteomic profiling of 44 primary HNSCC tumors, generating a spatial atlas of 19,471,501 cells across whole-slide tissue sections. Spatial niches and ecosystem states were identified through integrated computational analyses and evaluated for associations with tumor programs, clinicopathologic features, and patient outcomes. FINDINGS: HPV-negative tumors were enriched for fibroblast-rich, immune-poor niches associated with epithelial-mesenchymal transition and hypometabolic tumor programs, whereas HPV-positive tumors displayed more diverse immune, stromal, and vascular niche combinations and were enriched for immunogenic ecosystem states. Approximately 20% of HPV-positive tumors exhibited fibroblast-rich ecosystem architectures resembling HPV-negative disease and were associated with less favorable outcomes than other HPV-positive tumors of similar stage. In patient-derived co-culture models, extracellular matrix-associated fibroblasts were associated with epithelial-mesenchymal transition (EMT)-like tumor states, CD8+ T cell dysfunction, and chemotherapy resistance-associated phenotypes. CONCLUSIONS: Spatial ecosystem architecture is associated with clinically relevant heterogeneity beyond conventional HPV-based classification. Fibroblast-rich, immune-poor ecosystem states characterize a high-risk subset of HPV-positive tumors and may provide a framework for improved biological classification and risk stratification in HNSCC. FUNDING: This work was supported by the National Institutes of Health (R01CA291607 and R21CA267527-01) and the Feldstein Medical Foundation.

Humans

Cell-type-specific dysregulation of gene expression due to Chd8 haploinsufficiency during mouse cortical development.

Disruptive variants in the chromodomain helicase CHD8 are associated with risk for autism spectrum disorder (ASD). CHD8 haploinsufficiency is hypothesized to contribute to ASD by perturbing neurodevelopmental gene expression. However, insight into cell-type-specific transcriptional effects of CHD8 haploinsufficiency remains limited. We used single-cell and single-nucleus RNA sequencing to identify dysregulated genes in the embryonic and juvenile Chd8+/- mouse cortex. Chd8 and other ASD risk-associated genes showed a convergent expression trajectory conserved between mouse and human developing cortex, increasing from progenitor zones to the cortical plate. Genes associated with neurodevelopmental disorders or involved in chromatin remodeling and neuron projection development were dysregulated in Chd8+/- embryonic radial glia. Genes implicated in synaptic activity and organization were dysregulated in Chd8+/- postnatal excitatory cortical neurons, suggesting impaired synaptogenesis. Our findings reveal complex patterns of transcriptional dysregulation due to Chd8 haploinsufficiency, potentially with distinct impacts on progenitors and maturing neurons in the excitatory neuronal lineage.

Animals

Clinical translation of senescence-related pan-cancer multi-omics: tools for assessment and immunotherapy prediction.

Cellular senescence (CS) exerts dual roles in tumorigenesis, yet its pan-cancer molecular characteristics and clinical value remain unclear, hindering its translation to oncology and personalized therapy. To address the lack of specific and universal tools for senescence assessment and immunotherapy response prediction, this study systematically analyzed 1259 CS-related genes from the CellAge database across 31 cancer types by integrating multi-omics data, including bulk RNA-seq, single-cell/spatial transcriptomics, and CRISPR screening. We developed a rank-based algorithm SenScoreR (publicly available at https://gxhub.shinyapps.io/SenScoreR/ ) for senescence quantification, validated with 10 independent datasets, and constructed a machine learning-based predictive model CS.Sig for immunotherapy response. Results showed that tumors had significantly lower Rank-based Senescence Score (RSS) than normal tissues across 31 cancers (average diagnostic AUC = 0.895), with low RSS linked to poor survival; high RSS correlated with reduced genomic instability, enriched CD8⁺ T/NK cell/macrophage infiltration, upregulated PD-L1 expression, and elevated immune cytolytic activity. CS.Sig demonstrated robust performance in predicting ICI response (AUC = 0.716 across 10 cohorts), outperforming 13 existing signatures, while CRISPR screening identified 17 senescence-related targets (e.g., CEP55, PPP1CC) whose knockout enhanced anti-tumor immunity. Our findings clarify CS's role in maintaining tumor genomic stability and shaping immune microenvironments, and the developed SenScoreR, CS.Sig, and identified targets bridge basic CS research with clinical oncology, providing a translational resource and hypothesis basis for future experimental and clinical validation.

Journal Article

YAP1 induces hepatocellular carcinoma via DNA demethylation rather than by canonical driver gene mutations.

Large-scale genome sequencing analyses have identified driver gene mutations (DGMs) in most cancers as well as their associated tumorigenic mechanisms. However, a small fraction of cancers are not positive for these canonical DGMs, leaving the mechanisms underpinning their formation a mystery. We hypothesized that canonical DGM-negative cancers might be driven by activation of the transcriptional coactivator YAP1 that led to the induction of epigenetic changes. To test this theory, we established a mouse mosaic model of hepatocellular carcinoma (HCC) in which we induced YAP1-TEAD activation in a few hepatocytes. Whole-exome sequencing did not identify canonical DGMs in HCCs, but bisulfite sequencing revealed widespread DNA demethylation leading to the transcriptional activation of multiple oncogenes. Knockdown of the DNA demethylation-promoting gene, Tet1, attenuated HCC formation in these mice. Single-cell spatial transcriptomics identified a Tet1-high subpopulation of HCC cells that interacted with other hepatic cell types. Our mechanistic mouse data align with the observation that YAP1-TEAD-TET1-associated signatures were also elevated in hepatocytes from patients with Fontan-associated liver disease (FALD), a condition associated with the development of HCCs with lower frequencies of canonical DGMs. Our study suggests that the YAP1-TEAD-TET1 axis promotes canonical DGM-negative HCC development, and provides new insights into the molecular processes involved.

Animals