Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “transcriptomic”

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 1,045 records · Page 58Linked to original sources
PubMed2025

Unsupervised multiscale clustering of single-cell transcriptomes to identify hierarchical structures of cell subtypes.

BACKGROUND: Cell clustering is an essential step in uncovering cellular architectures in single-cell RNA sequencing (scRNA-seq) data. However, the existing cell clustering approaches are not well designed to dissect complex structures of cellular landscapes at a finer resolution. RESULTS: Here, we develop a multiscale clustering (MSC) approach to construct a sparse cell-cell correlation network for unsupervised identification of de novo cell types and subtypes across multiple resolutions. Based upon simulated silver- and gold-standard data as well as real scRNA-seq data in diseases, MSC demonstrates significantly improved performance compared to established benchmark methods and reveals a biologically meaningful cell hierarchy to facilitate the discovery of novel disease-associated cell subtypes and mechanisms. CONCLUSIONS: We present MSC as a new single-cell multiscale clustering framework as a powerful tool for advancing discoveries in disease-associated cell populations using single-cell sequencing data.

Single-Cell Analysis↗
PubMed2025

Integrative Genomic, Transcriptomic and Epigenomic Analysis Reveals cis-regulatory Contributions to High-altitude Adaptation in Tibetan Pigs.

The Qinghai-Tibet Plateau, characterized by its extreme environmental conditions, presents significant challenges to life, making it an ideal region for studying adaptation and evolution. Tibetan pigs, known for their high genetic diversity and exceptional adaptability to high altitudes, serve as excellent models for investigating high-altitude adaptation. While previous studies have extensively identified genetic determinants associated with high-altitude adaptation, the molecular mechanisms, particularly cis-regulatory patterns, remain poorly understood. Here, we conducted a selective sweep analysis using 484 genomes from Chinese and Western pig breeds across various altitudes, revealing 38.56 Mb of genomic regions under selection in Tibetan pigs. Enrichment analysis identified the lung as the primary functional tissue involved in high-altitude adaptation, supported by tissue-specific transcriptional and regulatory patterns observed between Tibetan and Meishan pigs (low altitude). By integrating genomic, RNA-seq, ATAC-seq, and H3K27ac HiChIP data, we constructed comprehensive enhancer-promoter regulatory maps of candidate genes and pinpointed promising genetic determinants associated with high-altitude adaptation, including SNPs in EPAS1, KLF13, SPRED1, and CFD. These loci were predicted to influence chromatin accessibility and the interactions of regulatory elements, with altered binding strength of relevant transcription factors. Further in vitro experiments confirmed that these loci function as allele-specific enhancers, modulating the expression of target genes. Our findings elucidate the regulatory basis of high-altitude adaptation in Tibetan pigs and provide valuable insights for exploring hypoxia-related diseases in livestock and humans.

Animals↗
PubMed2026

Transcriptomic shift in ethanol and amino acid metabolic genes regulated by Med15 during alcoholic fermentation.

Organisms that thrive in extreme environments provide natural experiments in evolution, revealing the genetic regulators that orchestrate complex phenotypic change. Wine yeast (WY) are specialized strains that are adapted to survive in the wine making environment while producing high concentrations of ethanol. In addition to large genomic changes that differentiate WY from yeast used in other industries, SNP and polyglutamine tract polymorphism in the transcriptional regulator Med15 are associated with the fermentation efficiency and stress response phenotypes of WY. In this study, we investigated the transcriptional differences during wine fermentation in transgenic lab strain yeast having integrated WY MED15 alleles. Compared to the unmodified lab strain (MED15 LAB), the same strain in which the MED15 locus was replaced with a MED15 allele from yeast isolated from palm wine, the fermented sap of palm (oil, date, coconut) trees (MED15 WY23), exhibited enhanced expression of amino acid biosynthesis genes as well as stress resistance and metabolic adaptation genes. Our experimental data confirm the role of arginine in efficient fermentation and suggest that certain MED15 alleles alter the expression patterns of arginine pathway genes in some cases improving carbon flux under nitrogen stress. The global benefits conferred by natural polymorphisms in a single transcriptional regulator highlight Med15 as a target for engineering of strains devoted to various types of alcohol production.

Ethanol↗
PubMed2025

Comparative Embryology and Transcriptomics of Asellus infernus, an Isopod Crustacean From Sulfidic Groundwater.

Sulfidic caves are harsh and extreme environments characterized by limited oxygen, low pH, and the presence of hydrogen sulfide. Amazingly, animals can live in sulfidic caves, one such animal being Asellus infernus, a representative of the Asellus aquaticus species complex, originating from Movile Cave and from old wells that represent windows of access to a sulfidic groundwater ecosystem located in southeast Romania. Little previous work has been done on lab-reared populations of A. infernus as they have been historically difficult to raise in the lab. Here, we develop resources for A. infernus, examining questions of timing of morphological differences in cave versus surface individuals, whether the environment (lab-bred vs. wild-caught) influenced size characteristics, and the genes and pathways showing differential expression between cave and surface samples. We found that A. infernus did not develop pigmentation embryonically, and juveniles had increased body length and longer antenna II as compared to surface individuals. Furthermore, we found that some of these measures differed between wild-caught and lab-reared juveniles for a given population, indicating that environmental differences can also influence these size characteristics. In addition, differential expression between cave and surface samples and allele-specific expression studies within F1 hybrids identified multiple genes, including those involved in sulfide metabolism and phototransduction. Strikingly, molecular convergence of genes involved in sulfide detoxification was observed between A. infernus and previous work on a fish that lives in both cave and sulfidic environments, Poecilia mexicana. In sum, we were able to develop embryonic and genomic tools for A. infernus, a model for understanding cave adaptation and adaptation to sulfidic environments.

Animals↗
PubMed2026

Single-Cell Transcriptome-Wide Mendelian Randomization and Colocalization Uncover Potential Immunocytes-Related Therapeutic Targets for Obesity.

Weight-loss treatment is crucial for individuals with obesity to prevent various complications. The role of Immune cells in obesity has been recently recognized, whereas its translation into therapy requires identifying key target genes. We performed Mendelian randomization (MR) analysis to assess causal relationships between expression quantitative trait loci (eQTL) of 14 immune cells and obesity-related traits (obesity, body mass index and body fat percentage), and validated the results in colocalization analysis. For the putative causal genes identified by the MR and colocalization analyses, we conducted pathway enrichment, differential expressed gene (DEG) analysis and search of druggable evidence, and utilized a Tier system to prioritize drug targets for obesity. MR and colocalization evidence was observed for 1630 genes associated with one or more obesity-related traits, mainly expressed in CD4+ naive/central memory T cells and enriched in antigen processing and presentation pathways. Forty-one genes showed causal relationship with all three outcomes, among which 19 genes have not been reported for obesity previously. DEG analysis using single-cell RNA sequencing data of blood or adipose tissue indicated that the differential expression of UBE2Z in monocytes, ZCCHC7 in T cells, and FNBP4 in B cells between lean and obese individuals were consistent with the MR results. By searching drug-gene interaction databases, we found targeted drugs for PYGB and PRUNE1, and PYGB was the top gene ranked in the Tier system. This study provides evidence for the involvement of immune cells in obesity, and the potential cell-specific, immune-related targets for obesity treatment.

Obesity↗
PubMed2025

Transcriptomic and enzymological evidence for plastid peptidoglycan synthesis in the gymnosperm Picea abies.

It is understood that a cyanobacterium was the progenitor of plastids and that the biosynthesis of cell wall peptidoglycan was lost during chloroplast evolution. However, accumulated data, especially from the moss Physcomitrium patens, suggest that peptidoglycan remains essential for plastid division in some land plants. A fundamental set of peptidoglycan biosynthesis (Mur) genes has been identified in the genomes of these land plants, while many angiosperms no longer encode some core Mur genes, including a bifunctional penicillin-binding protein (PBP). Ten incomplete Mur genes were previously identified in the genome of the gymnosperm Picea abies but these could be pseudogenes or encode proteins that have been repurposed. For instance, mutant albino maize and Arabidopsis seedlings possess a defective UDP-N-acetylmuramoyl-l-alanyl-d-glutamate--2,6-diaminopimelate ligase (MurE), an intact MurE ligase being essential for peptidoglycan synthesis. In this study, we isolated a full set of cDNAs for peptidoglycan biosynthesis from P. abies. GFP fusion proteins with either P. abies (Pa)MurE or PaPBP were detected in chloroplasts. Cross-species complementation assays with PaMurE in Arabidopsis albino MurE mutants and Physcomitrium MurE chloroplast division mutants showed that the gymnosperm MurE completely rescued both mutant phenotypes. Enzymatic assay of recombinant PaMurE proteins revealed they catalyze the same reaction performed by their bacterial MurE homologs. Moreover, the expression of the PaPbp cDNA partially rescued the giant chloroplast phenotype in the moss Pbp knockout line. These results are consistent with the operation of a functional Mur gene set in the Norway spruce genome.

Peptidoglycan↗
PubMed2026

Integrated transcriptomic and immunogenomic analysis unravels the immunological functions and prognostic landscape of WD repeat domain 76.

BackgroundWD Repeat Domain 76 (WDR76) plays a potential role in cellular regulation; however, its comprehensive landscape across human malignancies and its specific biological function in hepatocellular carcinoma (HCC) remain largely unexplored.MethodsWe conducted a systematic pan-cancer analysis utilizing multi-omics data from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Cancer Cell Line Encyclopedia (CCLE) atabases to evaluate WDR76 expression, subcellular localization, and its correlation with clinicopathologic features, genomic instability, and immune infiltration. Diagnostic and prognostic values were assessed via Receiver operating characteristic (ROC) and Kaplan-Meier analyses. Furthermore, the functional role of WDR76 in HCC was validated in vitro using Hep-3B and Huh7 cell lines through siRNA-mediated knockdown, followed by CCK-8, wound-healing, and transwell assays.ResultsWDR76 was significantly upregulated in the majority of tumor types, including LIHC, LUAD, and COAD, while exhibiting nuclear localization. Elevated WDR76 expression correlated with advanced tumor staging, metastasis, and poor clinical outcomes across multiple cohorts, particularly in ACC, KIRP, and LIHC. ROC analysis highlighted its exceptional diagnostic precision in cancers such as GBM and LIHC. Immunologically, WDR76 expression was intricately linked to immune cell infiltration, immune checkpoint markers, and genomic instability parameters, suggesting a role in shaping the tumor microenvironment. Drug sensitivity profiling revealed that high WDR76 levels correlate with resistance to specific chemotherapeutic agents. Experimentally, silencing WDR76 in HCC cells significantly suppressed cell proliferation, migration, and invasion capabilities.ConclusionOur study establishes WDR76 as a robust pan-cancer prognostic biomarker and a potential immunotherapeutic target. Specifically, we provide experimental evidence that WDR76 functions as an oncogenic driver in liver cancer, promoting malignant phenotypes and offering a novel avenue for targeted therapeutic intervention.

Humans↗
PubMed2026

Integrated phenotype, endogenous hormones and transcriptome analysis revealed the mechanism of response of Phoebe bournei seedlings to shade signals.

Understory tree seedlings are subjected to prolonged shading stress imposed by the canopy foliage, which significantly impedes their growth. A hallmark of shaded environments is a reduced red to far-red light ratio (R: FR). This study elucidates the physiological and molecular responses of the endangered tree species Phoebe bournei to shading signals. Seedlings were exposed to white light (control) and simulated shading environments with R: FR ratios of 1.5, 0.8, and 0.2. The findings reveal that an increase in the proportion of far-red light significantly enhances seedling height, root-collar diameter, internode length, petiole length, leaf surface area, and leaf biomass. Differentially expressed genes (DEGs) in each treatment group predominantly enrich pathways associated with hormone signaling, stress responses, and photosynthesis. Validation experiments demonstrate that shading promotes the activity of Rubisco and RCA enzymes, total chlorophyll (Chl) accumulation, and elevated levels of hormones including indole-3-acetic acid (IAA), gibberellic acid (GA3), salicylic acid (SA)/methyl salicylate (MeSA), cytokinins (CK), abscisic acid (ABA), and jasmonic acid (JA). Weighted Gene Co-expression Network Analysis (WGCNA) identifies seven hub genes linked to photosynthesis and plant hormone regulation: MYB, KSC, SUAR, CESA POD, CESA, and SAUR. Collectively, shading signals induce P. bournei seedlings to elongate their stems and petioles, enhance photosynthetic enzyme activity, and accumulate specific hormones, with pertinent genes actively participating in light signal transduction. This research sheds light on the shading response mechanism of P. bournei, providing a robust theoretical framework for the breeding of shade-tolerant trees and the conservation of endangered species.

Transcriptome↗
PubMed2024

Construction of a molecular diagnostic system for neurogenic rosacea by combining transcriptome sequencing and machine learning.

Patients with neurogenic rosacea (NR) frequently demonstrate pronounced neurological manifestations, often unresponsive to conventional therapeutic approaches. A molecular-level understanding and diagnosis of this patient cohort could significantly guide clinical interventions. In this study, we amalgamated our sequencing data (n = 46) with a publicly accessible database (n = 38) to perform an unsupervised cluster analysis of the integrated dataset. The eighty-four rosacea patients were partitioned into two distinct clusters. Neurovascular biomarkers were found to be elevated in cluster 1 compared to cluster 2. Pathways in cluster 1 were predominantly involved in neurotransmitter synthesis, transmission, and functionality, whereas cluster 2 pathways were centered on inflammation-related processes. Differential gene expression analysis and WGCNA were employed to delineate the characteristic gene sets of the two clusters. Subsequently, a diagnostic model was constructed from the identified gene sets using linear regression methodologies. The model's C index, comprising genes PNPLA3, CUX2, PLIN2, and HMGCR, achieved a remarkable value of 0.9683, with an area under the curve (AUC) for the training cohort's nomogram of 0.9376. Clinical characteristics from our dataset (n = 46) were assessed by three seasoned dermatologists, forming the NR validation cohort (NR, n = 18; non-neurogenic rosacea, n = 28). Upon application of our model to NR diagnosis, the model's AUC value reached 0.9023. Finally, potential therapeutic candidates for both patient groups were predicted via the Connectivity Map. In summation, this study unveiled two clusters with unique molecular phenotypes within rosacea, leading to the development of a precise diagnostic model instrumental in NR diagnosis.

Humans↗
PubMed2024

Transcriptome-wide N6-methyladenosine modification profiling of long non-coding RNAs in patients with recurrent implantation failure.

N6-methyladenosine (m6A) is involved in most biological processes and actively participates in the regulation of reproduction. According to recent research, long non-coding RNAs (lncRNAs) and their m6A modifications are involved in reproductive diseases. In the present study, using m6A-modified RNA immunoprecipitation sequencing (m6A-seq), we established the m6A methylation transcription profiles in patients with recurrent implantation failure (RIF) for the first time. There were 1443 significantly upregulated m6A peaks and 425 significantly downregulated m6A peaks in RIF. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that genes associated with differentially methylated lncRNAs are involved in the p53 signalling pathway and amino acid metabolism. The competing endogenous RNA network revealed a regulatory relationship between lncRNAs, microRNAs and messenger RNAs. We verified the m6A methylation abundances of lncRNAs by using m6A-RNA immunoprecipitation (MeRIP)-real-time polymerase chain reaction. This study lays a foundation for further exploration of the potential role of m6A modification in the pathogenesis of RIF.

Humans↗
PubMed2024

Effects of different temperatures on chondrocyte growth: a transcriptomic analysis.

BACKGROUND: Our previous study demonstrated that temperature-related microwave ablation (MWA) can safely modulate growth plates of piglets' vertebrae. Therefore, this study is designed to investigate the effects of different temperatures on chondrocyte viability and the underlying molecular mechanisms in vitro. METHODS: Following a 10-minute treatment at different temperatures (37 °C, 40 °C, 42 °C, 44 °C, 46 °C, 48 °C, and 50 °C), CCK-8 assay was used to examine the viability of ATDC5 cells at 12 h. Differentially expressed genes (DEGs) and the hub genes in ATDC5 cells treated at 37 °C, 40 °C and 44 °C were identified using RNA-seq. The expression of hub genes in ATDC5 cells was validated using RT-qPCR. RESULTS: Compared with 37 °C, exposure to 40 °C significantly increased the viability of ATDC5 cells, while 42 °C had no significant effect. Additionally, exposure to 44 °C, 46 °C, 48 °C, and 50 °C exhibited the opposite pattern, with ATDC5 cells being particularly less than 50% active after treatment at 46 °C, 48 °C, and 50 °C. Differential expression analysis identified 179, 374 and 221 DEGs in the comparisons of 40 °C vs. 37 °C, 44 °C vs. 37 °C, and 44 °C vs. 40 °C, respectively. These DEGs predominantly regulated proliferation, differentiation, necrosis, inflammatory and immune responses, and ECM synthesis/degradation. Furthermore, they were associated with the Ras, PI3K/AKT, mTOR, cAMP, and MAPK pathways. Agt, Hspa1a, Hspb1, and Nlrc4 were identified as hub genes in DEGs, and RT-qPCR confirmed that the mRNA expression patterns of these hub genes in ATDC5 cells were largely consistent with the RNA-seq results. CONCLUSION: The regulation of chondrocyte viability by temperature is associated with Ras, PI3K/AKT, mTOR, cAMP, and MAPK pathways. Additionally, Agt, Hspa1a, Hspb1, and Nlrc4 may be the key regulatory genes in this process.

Chondrocytes↗
PubMed2025

Transcriptomics-based exploration of ubiquitination-related biomarkers and potential molecular mechanisms in laryngeal squamous cell carcinoma.

BACKGROUND: One of the most common and prevalent cancers is laryngeal squamous cell carcinoma (LSCC), which poses a great threat to the life and health of the patient. Nonetheless, it has been demonstrated that ubiquitination is crucial for the development and course of LSCC. Therefore, it is particularly important to identify biomarkers for ubiquitination-related genes (UbRGs) in LSCC. METHODS: Differentially expressed genes (DEGs) in the LSCC versus controls were obtained by differential expression analysis. Also, key modular genes associated with LSCC were obtained using weighted gene co-expression network analysis (WGCNA). Next, DEGs, key module genes, and UbRGs were taken to intersect to obtain candidate genes. And then machine algorithms were to screen potential biomarkers, further their diagnostic value were analyzed and validated. Then, therapeutic agents for biomarkers were predict. In addition, the regulatory networks of the biomarkers were mapped. The expression levels of biomarkers were detected in clinical samples using reverse transcription-quantitative PCR (RT-qPCR). RESULTS: A total of eight candidate genes were acquired by the overlap 1,911 DEGs, the key modular genes of WGCNA, and 1,393 UbRGs. A sum of four biomarkers (WDR54, KAT2B, NBEAL2 and LNX1) were identified by two machine learning, then these four biomarkers were validated in GSE127165 and the expression trend was consistent with TCGA-LSCC, they were recorded as biomarkers. Moreover, the accuracy of the biomarkers in predicting clinical aspects of LSCC was confirmed by the receiver operating characteristic (ROC) curves. Subsequently, cancers such as malignant neoplasms, colorectal cancers, tumors, and primary malignant neoplasms were significantly associated with the biomarkers, which further suggests that these four biomarkers were strongly associated with cancer. Meanwhile, the drugs garcinol, cocaine, and triazolam, among others, used for LSCC treatment were predicted. Finally, transcription factors (TFs) (BRD4, MYC, AR, and CTCF) were predicted to regulate the biomarkers. RT-qPCR assays illustrated that the expression trends of KAT2B, LNX1 and NBEAL2 remained consistent with the dataset. CONCLUSION: The identification of four biomarkers (WDR54, KAT2B, NBEAL2 and LNX1) associated with UbRGs could ultimately serve as a predictive clinical diagnosis of LSCC and provide insight into the molecular mechanisms of LSCC.

Humans↗
PubMed2026

TNFα-induced endothelial extracellular vesicles regulate astrocyte function: an integrated transcriptomic and proteomic study.

Endothelial cells and astrocytes are critical structural and functional components of the blood-brain barrier. In many neuroinflammatory diseases, endothelial cells are among the first to respond to inflammatory stimuli and release extracellular vesicles (EVs). However, whether inflammatory stimulation alters EV RNA cargo and subsequently regulates astrocyte function remains unclear. In this study, we performed integrated RNA sequencing and proteomic analyses to investigate the effects of TNFα-stimulated endothelial EVs on astrocytes. RNA profiling revealed significant alterations in EV cargo after TNFα stimulation, including 867 upregulated and 577 downregulated mRNAs, 317 upregulated and 15 downregulated lncRNAs, and 88 upregulated and 62 downregulated miRNAs. The results of functional enrichment analysis suggested that altered EV RNAs may primarily promote inflammatory responses, cell migration, and RNA splicing in astrocytes while reducing their regulatory effects on neuronal projection and calcium homeostasis. Further integrative analysis of EV RNAs and astrocytic proteomics revealed key overlapping targets, including upregulated expression of ICAM1, SOD2, TFPI2, and TNFAIP8, whereas NFKBIA expression was consistently decreased. Network analysis revealed NF-κB as the central regulatory node. Reduced levels of EV-derived NFKBIA mRNA were associated with decreased IκBα protein levels in astrocytes, which promoted NF-κB activation and inflammatory cytokine release. Finally, overexpression of IκBα in astrocytes significantly attenuated TNFα EV-induced IL-1β and IL-6 secretion. Collectively, these findings demonstrate that TNFα-stimulated endothelial EVs coordinately regulate astrocyte function through mRNA, lncRNA, and miRNA cargo and that the IκBα/NF-κB axis may be a key mechanism underlying endothelial EV-mediated inflammatory disruption of the blood-brain barrier.

Astrocytes↗
PubMed2026

Transcriptome-based epigenetic screening identifies DNA hypermethylation signatures as prognostic biomarkers in oral squamous cell carcinoma.

Promoter DNA hypermethylation is a key epigenetic mechanism of gene silencing in cancer, yet the DNA hypermethylome of oral squamous cell carcinoma (OSCC) and its prognostic relevance remain poorly characterized. Here, we systematically identified and validated novel hypermethylated genes with prognostic significance in OSCC using a genome-wide discovery and multi-platform validation strategy. Candidate genes were first identified by pharmacologic demethylation combined with RNA sequencing across OSCC cell lines, then validated by quantitative RT-PCR, methylation-specific PCR, and bisulfite sequencing in OSCC cell lines, normal oral mucosa, and primary OSCC tumors, with independent confirmation in the TCGA-HNSC dataset. Immunohistochemistry confirmed protein-level silencing, and Kaplan-Meier survival analysis assessed prognostic significance across both cohorts. This pipeline identified five candidate genes, GPX3, ANG, CTGF, GPRC5B, and BAMBI, exhibiting cancer-specific promoter hypermethylation associated with transcriptional and protein silencing in OSCC. Validation in oral cavity tumor samples extracted from the TCGA-HNSC dataset confirmed tumor-specific hypermethylation and revealed significant inverse correlations between methylation and expression for GPX3, GPRC5B, and CTGF. Notably, CTGF hypermethylation was independently associated with poor overall survival in both cohorts (institutional cohort, p=0.03; oral tumor subset from TCGA-HNSC, p=0.01), and a combined ANG+CTGF methylation signature showed superior and reproducible prognostic performance across both platforms. Pathway analysis linked these genes to epithelial-mesenchymal transition and interferon response signaling. This study establishes the first validated DNA methylation biomarker panel for OSCC prognosis, identifying CTGF hypermethylation as a robust prognostic driver with translational potential for clinical risk stratification.

Humans↗
PubMed2025

Revealing cancer driver genes through integrative transcriptomic and epigenomic analyses with Moonlight.

Cancer involves dynamic changes caused by (epi)genetic alterations such as mutations or abnormal DNA methylation patterns which occur in cancer driver genes. These driver genes are divided into oncogenes and tumor suppressors depending on their function and mechanism of action. Discovering driver genes in different cancer (sub)types is important not only for increasing current understanding of carcinogenesis but also from prognostic and therapeutic perspectives. We have previously developed a framework called Moonlight which uses a systems biology multi-omics approach for prediction of driver genes. Here, we present an important development in Moonlight2 by incorporating a DNA methylation layer which provides epigenetic evidence for deregulated expression profiles of driver genes. To this end, we present a novel functionality called Gene Methylation Analysis (GMA) which investigates abnormal DNA methylation patterns to predict driver genes. This is achieved by integrating the tool EpiMix which is designed to detect such aberrant DNA methylation patterns in a cohort of patients and further couples these patterns with gene expression changes. To showcase GMA, we applied it to three cancer (sub)types (basal-like breast cancer, lung adenocarcinoma, and thyroid carcinoma) where we discovered 33, 190, and 263 epigenetically driven genes, respectively. A subset of these driver genes had prognostic effects with expression levels significantly affecting survival of the patients. Moreover, a subset of the driver genes demonstrated therapeutic potential as drug targets. This study provides a framework for exploring the driving forces behind cancer and provides novel insights into the landscape of three cancer sub(types) by integrating gene expression and methylation data.

Humans↗
PubMed2026

Integrated single-cell transcriptomics, Mendelian randomization, and machine learning identify CEBPZ as an immune-related biomarker in oral lichen planus.

BACKGROUND: Oral lichen planus (OLP) is a chronic, immune-mediated oral mucosal disease with complex pathophysiology and potential for malignant transformation. Understanding its molecular basis is critical for the development of precise diagnostic and therapeutic strategies. OBJECTIVES: We aimed to identify key immune-related biomarkers and characterize cellular dynamics in OLP, with a particular focus on the role of CEBPZ in disease pathogenesis. MATERIAL AND METHODS: We analyzed single-cell RNA sequencing (scRNA-seq) data from OLP lamina propria samples (GSE211630) to identify disease-specific T-cell subpopulations using high-dimensional weighted gene co-expression network analysis (hdWGCNA) for oxidative stress-related gene modules.-data-based Mendelian randomization (SMR) integrated FinnGen genome-wide association study (GWAS; 342,499 Europeans) data with Genotype-Tissue Expression (GTEx) expression quantitative trait loci (eQTL) data to identify causal genes. Machine learning (ML) models (least absolute shrinkage and selection operator (LASSO) and convolutional neural network (CNN)) were developed using bulk RNA-seq datasets (GSE52130 and GSE38616) for diagnostic purposes. RESULTS: We identified OLP-specific T-cell populations (clusters 0, 3, 5, 7, 13, and 15) with enhanced migration inhibition factor (MIF) pathway signaling toward B cells and monocytes. Two oxidative stress-associated modules contained hub genes, including CEBPZ. Summary-data-based Mendelian randomization analysis identified 231 OLP-associated genes, with CEBPZ uniquely intersecting LASSO-selected markers (odds ratio (OR) = 1.057, 95% confidence interval (95% CI) = 1.013-1.102, p = 0.010). Machine learning models achieved area under the curve (AUC) values ranging from 0.653 to 0.745, with the CNN model reaching a validation accuracy of 0.735. CEBPZ showed elevated expression in OLP T cells and correlated with enhanced MIF-(CD74+CXCR4) signaling. CONCLUSIONS: This integrative approach identifies CEBPZ as a pivotal biomarker linking genetic susceptibility, oxidative stress, and immune dysregulation in OLP. Our diagnostic models offer promising tools for OLP management.

CEBPZ↗
PubMed2026

Meta-Merging the Transcriptomes of Gastric Tumors Redefines the Connections among Molecular and Clinical Subtypes.

INTRODUCTION: The availability of a large number of cancer expression profiles presents an excellent opportunity to re-investigate various biological and clinical questions. While several expression profiles have been established for different cancers, merging them may provide a more powerful platform for extensively extrapolating molecular and clinical features across multiple cohorts. MATERIALS AND METHODS: In this study, five gastric tumor expression profiles from the Gene Expression Omnibus [GEO] and one in-house cohort comprising a total of 1,060 samples were merged. The batch effect was removed using non-parametric ComBat analysis, and the seamless merging of datasets was confirmed through various parameters. RESULTS: Extrapolation of ACRG [Asian Cancer Research Group] and TCGA [The Cancer Genome Atlas] molecular subtypes in the merged cohort of 1,060 gastric tumors revealed nine distinct clusters. Notably, the following patterns were observed: [i] mutual exclusivity between Epithelial to Mesenchymal Transition [EMT] and Microsatellite Instability [MSI] subtypes in 90% of tumors; [ii] overlapping occurrence of EMT and MSI subtypes in the remaining tumors; [iii] overlap between MSI and Epstein-Barr Virus [EBV] subtype tumors; [iv] both commonalities and differences between EMT and Genomically Stable [GS] subtypes; and [v] an association between EBV positivity and PI3K mutation. CONCLUSION: The current study demonstrates that compiling a larger expression profile is valuable for revisiting the molecular features and epidemiology associated with molecular subtypes, thereby aiding in the development of novel diagnostics and targeted therapeutics.

Humans↗
PubMed2026

Single-cell transcriptomics reveals distinct microglial state remodeling associated with the (R)-nicotine/diosmetin combination and galantamine in LPS-challenged BV2 cells.

BACKGROUND: Microglial neuroinflammation contributes to the progression of neurodegenerative diseases, yet it remains challenging to attenuate inflammatory responses while preserving cellular function. The effects of (R)-nicotine, diosmetin, their combined administration, and galantamine on heterogeneous BV2 transcriptional states have not been compared at single-cell resolution. METHODS: LPS-stimulated BV2 microglial-like cells were treated with (R)-nicotine, diosmetin, their combination (DR), or galantamine. Single-cell RNA sequencing was performed with three biological replicates per group and integrated with RNA velocity and SCENIC regulon inference to characterize treatment-associated state redistribution, inferred local transcriptional directionality and regulon-activity patterns. Functional validation included CCK-8 metabolic activity assays, multiplex cytokine ELISA, BDNF/GDNF quantification, qPCR, and high-content immunofluorescence analysis of iNOS and Arg1 at single-cell resolution. RESULTS: LPS decreased the relative abundance of the Itgae+/Plk4+ cluster while increasing the Nmur1+/Limk2+ cluster and inflammatory effector programs. DR treatment suppressed pro-inflammatory cytokine release without significantly reducing CCK-8-assessed metabolic activity, increased the Itgae+/Plk4+ cluster proportion, and increased BDNF/GDNF relative to LPS. RNA velocity and SCENIC analyses suggested that DR and galantamine showed distinct transcriptional and regulatory patterns: DR attenuated Batf-associated inflammatory regulons and was associated with increased Atf3-linked stress-response activity, whereas galantamine preferentially engaged DNA repair and genome-maintenance programs. CONCLUSION: These findings indicate that the DR condition was associated with remodeling of LPS-challenged BV2 microglial-like states, attenuation of inflammatory programs, and increased neurotrophic outputs relative to LPS, without significantly reducing CCK-8-assessed metabolic activity. This study provides a single-cell characterization of distinct treatment-associated responses to (R)-nicotine, diosmetin, their combined administration, and galantamine.

Microglia↗