SEARCH · Search PubMed

Results for “bulk RNA sequencing (bulk RNA-seq)”

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.

SourceResearch connectionFrom yearThrough year

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

34 records · Page 2Linked to original sources

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

Dysregulation of U12-Type Splicing in Lupus Neutrophils.

OBJECTIVE: Neutrophil dysfunction is a hallmark of systemic lupus erythematosus (SLE), but its molecular basis remains unclear. This study explores transcriptional and posttranscriptional changes in low-density granulocytes (LDGs), a proinflammatory neutrophil subset expanded in SLE, focusing on NADPH oxidase (Nox) function and minor intron splicing. METHODS: LDGs and normal-density granulocytes (NDGs) were isolated from patients with SLE and healthy controls (HCs). CYBA (p22phox) expression was evaluated at transcript and protein levels. Nox activity was measured using luminol assays. Bulk RNA sequencing (RNA-seq) and rMATS software were used to assess alternative splicing, particularly of U12-type intron-containing genes. RESULTS: CYBA expression was reduced in SLE LDGs (n = 11) compared to SLE and HC NDGs (n = 6), with levels resembling those in chronic granulomatous disease neutrophils. SLE LDGs exhibited impaired Nox activity (n = 7 SLE, n = 12 HC). CYBA is a U12 intron-containing gene, and transcriptomic analysis revealed broad down-regulation of this gene class in SLE LDGs, suggesting minor spliceosome dysfunction. rMATS analysis showed increased U12-type intron retention and widespread splicing defects-including exon skipping and mutually exclusive exon use-in genes such as GBP5, MAEA, and STX10. These abnormalities were validated in an independent long-read RNA-seq data set from SLE peripheral blood mononuclear cells. Importantly, splicing disruptions correlated with disease activity and autoantibody profiles. CONCLUSION: Impaired U12-dependent splicing may contribute to neutrophil dysfunction in SLE, potentially via defective oxidative burst and altered immune regulation. These findings highlight the minor spliceosome as a novel player in lupus pathogenesis.

Humans↗

PubMed2026

Integrated Genomic and Proteomic Analysis Reveals T-B Lymphocyte Signatures in the MYCN Driven "Immune Desert" of Specific Neuroblastoma Subtypes.

AIMS: This study aims to systematically dissect how MYCN amplification shapes the immunosuppressive tumor microenvironment (TME) in high-risk neuroblastoma, elucidating key mechanisms underlying immune evasion. METHODS: We performed an integrated multi-omics analysis of bulk RNA-seq (n = 721), single-cell RNA-seq (n = 9), proteomic data (n = 49) and spatial transcriptomics (Visium, with external validation in melanoma). Analyses included unsupervised clustering, cell-cell communication inference, transcriptional regulatory network reconstruction, and spatial proximity assessment to map the immune landscape. RESULTS: A distinct molecular subtype (Class C), defined by MYCN amplification and poor prognosis, exhibited a comprehensive "immune desert" phenotype characterized by low immune scores and minimal leukocyte infiltration. Single-cell analysis confirmed significant depletion of T and B lymphocytes within the Class C TME. Dysregulated transcriptional networks were identified, including upregulation of REL and EOMES in T cells-with EOMES potentially driving exhaustion via regulation of Transient Receptor Potential (TRP) genes, and REL inhibition enhancing cytotoxic function in vitro. A unique immunosuppressive B-cell subset (B7) engaged in enhanced crosstalk with exhausted T cells and harbored a MYC-centered network linked to cell cycle dysregulation and poor survival. Spatial transcriptomics revealed significant proximity between B7-active regions and Treg/exhaustion-enriched areas, externally validated in melanoma. Proteomic data validated elevated REL expression in MYCN-amplified tumors. CONCLUSION: This work delineates the immunosuppressive architecture of MYCN-driven neuroblastoma, revealing novel regulatory nodes within specific lymphocyte compartments. Integrating single-cell, spatial, and proteomic evidence, we propose REL inhibition as a therapeutic candidate, the EOMES/TRP axis as a bioinformatically supported hypothesis, and the B7/MYC hub as a hypothesis supported by transcriptomic and spatial evidence.

Humans↗

PubMed2026

Integrated single-cell and bulk transcriptomic analysis identifies a novel senescent fibroblast subtype associated with poor prognosis in acral melanoma.

BACKGROUND: Acral melanoma (AM) exhibits significant intratumoral heterogeneity, but its tumor microenvironment (TME) and immune regulation remain unclear. This study aims to dissect TME heterogeneity and establish a prognostic model based on key cell subpopulations. METHODS: We collected AM single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). Unsupervised clustering, CellChat, and Scissor analysis were performed to characterize cellular heterogeneity, cell-cell communication, and prognosis-related cell subpopulations. Kaplan-Meier analysis was used to assess the prognostic value of key genes, which were further validated by multiplex immunohistochemistry (mIHC). RESULTS: In AM, Mel_C2, C7, and C9 with high SEMA6A and KIT expression were strongly linked to poor prognosis. We further identified a senescent fibroblast subpopulation (sCAF_CDKN2A) characterized by high fibroblast senescence signature (FSS) scores. Integrating Scissor analysis of fibroblast subtypes with bulk prognostic data, we identified COL3A1, VCAN, and KIT as prognosis-associated genes upregulated in poor-outcome-related fibroblast subsets. Cell-cell communication analysis revealed that sCAF_CDKN2A engages in an immunosuppressive network, interacting with regulatory T cells (Tregs) via MIF signaling and receiving signals from exhausted CD8+ T cells through PPIA-BSG interactions. Using transcription factor expression patterns from these fibroblast subtypes, we constructed a prognostic model that effectively stratified patients into distinct risk groups with significant differences in overall survival (OS). mIHC confirmed significantly higher protein levels of SEMA6A and COL3A1 in tumor tissues compared to matched normal tissues. CONCLUSIONS: We established a novel prognostic model for AM and identified sCAF_CDKN2A as an immunosuppressive senescent fibroblast subpopulation driving poor prognosis.

Acral melanoma↗

PubMed2026

Single-cell multiomics reveals exosome-mediated reprogramming and clonotypic remodeling of T cells in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive and immunogenic subtype lacking targeted therapies. While tumor-derived exosomes are known to modulate immune function, their direct impact on human T cell plasticity and antigen specificity remains poorly defined. Here, we conducted a comprehensive single-cell multiomic analysis of primary human T cells exposed to exosomes derived from 17 genomically diverse TNBC cell lines and 35 patient samples. Integrating single-cell RNA-seq, V(D)J sequencing, non-coding RNA profiling, bulk and single-cell cytokine analyses, we uncovered conserved and subtype-specific immunomodulatory programs induced by TNBC exosomes. Exosome-treated T cells displayed skewing toward regulatory and dysfunctional phenotypes, including Th17-like, Treg, and PD-1⁺/PD-L1⁺ Tfh cells. Functional profiling revealed suppression of early activation markers and cytokine responses, alongside selective preservation of cytotoxic features in γδ T and NKT subsets. Transcriptomic and miRNA network analyses demonstrated widespread downregulation of immune effector genes (e.g., HBEGF and TNFSF9) mediated by exosome-delivered regulatory miRNAs (has-miR-98-5p). Notably, exosome-stimulated T cells displayed distinct clonotypic expansions, characterized by the emergence of five tumor-specific γδ TCR clonotypes and 30 unique αβ TCR CDR3 sequences that were absent in mock-treated controls, underscoring the role of exosomes in shaping TCR repertoire dynamics.

Humans↗

PubMed2026

Translation of scRNA-seq to a clinical blood test for infection diagnostics.

INTRODUCTION: Early and accurate triage of patients with febrile illness is crucial for appropriate treatment. While standard inflammatory biomarkers are often nonspecific, transcriptome analysis of peripheral blood has diagnostic potential. However, bulk gene expression data is often confounded by changes in cell count proportions, a more robust quantification of gene expression in specific single-cell types, such as monocytes, is required to serve as a reliable clinical biomarker. AREAS COVERED: Various methods to obtain single-cell-type gene expression results, including the gold standard of gene expression analysis after cell sorting and single-cell RNA sequencing, which are difficult to implement in the routine settings are discussed. Other method to interrogate gene expression of a single cell-type is needed. Finally, monocyte cell-type specific ratio-based biomarker (RBB, called Direct Leukocyte Single cell-type Transcript Abundance, or DIRECT LS-TA) which can estimate single cell-type (monocyte) specific gene expression without cell sorting is introduced. EXPERT OPINION: Traditional diagnostic test for differentiating infection has several limitations requiring breakthrough including turn-around time and cost. DIRECT LS-TA provides a reliable way to quantify monocyte-specific gene expression that strongly correlates with gold-standard methods. It is more affordable than single-cell RNA sequencing and can be readily implemented in clinical laboratories using widely available quantitative PCR or digital PCR machines.

Humans↗

PubMed2026

Nested co-expression network analysis identifies compact gene clusters in a black box.

MOTIVATION: Digital analysis of biological systems requires methods capable of identifying both broad and nested gene modules reflecting complex biological processes. Existing transcriptomic methods often miss compact gene sets corresponding to subprocesses in specialized cell types, limiting insights into functional heterogeneity. RESULTS: We present Nested-WGCNA, a two-stage unsupervised network analysis algorithm designed to identify coarse-grained and fine-grained gene modules. Applied to bulk RNA-Seq data, Nested-WGCNA reveals stable modules reproducible across datasets. When validated against scRNA-Seq data, these modules correspond to both major and minor immune cell subtypes. Application to immunotherapy response datasets uncovers predictive and prognostic biomarkers, highlighting its utility in treatment stratification and biomarker discovery. AVAILABILITY: The NestedWGCNA source code and analysis pipeline are available on GitHub (https://github.com/ilyada/NestedWGCNA) and archived on Zenodo (https://doi.org/10.5281/zenodo.18959244).

Algorithms↗

PubMed2026

Multiomics approaches reveal direct NF-κB p65 target genes in pancreatic islets during cytokine exposure and in type 1 diabetes.

Autoimmune diseases, including Type 1 diabetes (T1D), are often characterized by overactive inflammatory signaling pathways. The proinflammatory cytokine interleukin-1β (IL-1β) elicits global gene expression changes in islet β-cells which overlap with islets obtained from human donors with T1D. The direct transcriptional link between NF-κB subunit p65 and target genes involved with autoimmune events was investigated. We used a multiomics approach including bulk RNA-sequencing (RNA-Seq), single-cell RNA-sequencing (scRNA-Seq), and chromatin immunoprecipitation coupled to deep sequencing (ChIP-Seq), alongside molecular docking simulations, and transcriptional assays. Through the various experimental modalities, we identified early response genes driven by IL-1β that were differentially expressed in pancreatic islets from human T1D donors and also conserved across mouse, rat, and human tissues. ChIP-Seq revealed genes that are direct genomic targets of the NF-κB p65 transcription factor. Moreover, regions that gained RNA polymerase II binding following cellular exposure to IL-1β were identified, complementing the early response gene profile induced by β-cell exposure to IL-1β. Molecular docking simulations predicted that mutations reducing p65 transcriptional capacity do not alter DNA binding ability. These findings clearly show that IL-1β signaling in pancreatic β-cells directs p65 to specific genomic regions congruent with increased gene expression relevant to T1D in β-cell lines as well as mouse and human islets exposed to cytokines. Islets from human donors with T1D express genes identified as direct p65 targets using unbiased approaches, implicating heightened NF-κB activity as a critical component of autoimmune disease etiology.NEW & NOTEWORTHY Using multiple Seq-based approaches, this study identified genes expressed in human pancreatic tissue from donors with Type 1 diabetes that are regulated acutely by exposure to the cytokine interleukin-1beta. The NF-kB transcription factor p65 (RelA) was determined via ChIP-Seq to be a major control node regulating this immediate early response. These collective datasets are consistent with a paradigm of overactive NF-kB signaling as a critical component of autoimmunity in both rodents and humans.

Humans↗

PubMed2026

FZD5 drives macrophage-mediated immunomodulation and predicts prognosis in glioma: evidence from single-cell sequencing.

BACKGROUND: Gliomas are highly malignant brain tumors characterized by an immunosuppressive microenvironment, which limits therapeutic efficacy and contributes to poor clinical outcomes. The WNT/β-catenin signaling pathway is critically involved in tumor progression, and FZD5, a key receptor within this pathway, may participate in immune regulation. However, its specific role and underlying mechanisms in glioma remain unclear. METHODS: RNA-seq and microarray datasets from the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA), together with single-cell RNA sequencing (scRNA-seq) datasets from GEO, were comprehensively analyzed. The Seurat package was used to identify macrophage-related clusters and mitophagy-associated pathways. Cox and LASSO regression analyses, along with a prognostic nomogram, were applied to evaluate the prognostic significance of FZD5. Immune infiltration, functional enrichment, and immunotherapy response analyses were conducted, followed by validation using spatial transcriptomics, immunohistochemistry, and in vitro assays. RESULTS: In bulk glioma transcriptomes, FZD5 emerged as an independent predictor of poor prognosis. Crucially, single-cell and spatial analyses revealed that the biologically significant FZD5 signal originated predominantly within tumor-associated macrophages (TAMs), where it colocalized with the M2 marker CD163. Consistently, elevated FZD5 levels correlated with increased myeloid infiltration and an immunosuppressive tumor microenvironment. Functionally, macrophage-expressed FZD5 was associated with mitophagy-related programs and promoted an M2-skewed phenotype, thereby enhancing glioma cell proliferation, migration, and invasion via macrophage-glioma crosstalk. CONCLUSION: FZD5 is a TAM-enriched marker in glioma tissues and a potential regulator of macrophage-associated immunosuppressive programs, supporting its utility as a prognostic biomarker and a candidate target for microenvironment-oriented interventions in glioma.

Humans↗

PubMed2026

SCMO: a deep learning model integrating the single-cell resolution TME ecosystem and multi-omics for survival prediction in CRC patients.

BACKGROUND: Colorectal cancer (CRC) remains a leading cause of global cancer mortality, highlighting the need for precise survival prediction to guide clinical decisions. Although tissue-level multi-omics is widely utilized for survival prediction, its limited resolution cannot capture tumor heterogeneity. Single-cell RNA sequencing (scRNA-seq) enables dissection of the tumor microenvironment (TME) at cellular resolution, supporting personalized prognostic assessment. METHODS: We collected 213 CRC scRNA-seq samples and established a CRC-specific TME atlas comprising 339,060 cells. Using this atlas as a reference, we deconvolved bulk RNA-seq data from TCGA-CRC cohort with the EcoTyper algorithm to reconstruct TME features. Clinical, genomic, and transcriptomic data were obtained from the Xena platform; microbial data were sourced from the BIC database. We integrated TME and multi-omics features through a self-normalizing neural network to construct a deep learning model (single-cell resolution TME ecosystem with multi-omics data [SCMO]) for survival prediction. To enhance interpretability, we utilized the Integrated Gradients algorithm and spatial transcriptomic data to analyze multi-omics and TME features. We performed anticancer drug screening with tumor necrosis factor receptor-associated protein 1 (TRAP1), a critical feature according to the Integrated Gradients algorithm, as a potential target. RESULTS: We identified 13 survival-related TME features from the CRC-specific atlas: 12 cell states and one multi-cellular ecosystem. SCMO, which combined TME and multi-omics features, improved survival prediction and outperformed existing methods, achieving a concordance index of 0.762. The SCMO demonstrated robust performance for long-term predictions, achieving areas under the curve (AUCs) of 0.752, 0.772, and 0.869 for 1-, 3-, and 5-year predictions in the training set, with corresponding test set AUCs of 0.639, 0.756, and 0.772. TME features from the SCMO model revealed that ecosystem density increased with CRC malignancy. Multi-omics features included TRAP1 as a potential drug target. Drug screening identified saikosaponin A as a novel TRAP1 inhibitor, and its anticancer activity was validated in vitro. We developed SCMO-Lite, a simplified model incorporating 12 high-attribution-weight multi-omics features, which demonstrated robust risk stratification. CONCLUSIONS: SCMO combines analytical precision with biological interpretability, offering novel insights for oncology survival prediction.

Humans↗

PubMed2026

Radiogenomics predicts immune microenvironment heterogeneity and response to combination immunotherapy in hepatocellular carcinoma.

BACKGROUND: The combination of immune checkpoint inhibitors (ICIs) with anti-angiogenic agents is the preferred first-line therapy option for patients with advanced hepatocellular carcinoma (HCC), yet only a subset of patients responds, urging the quest for prediction biomarkers. We aimed to integrate genomics with radiology to propose an immune-derived radiogenomics biomarker of response to such combination immunotherapy and evaluate its added value in clinical context. METHODS: We integrated bulk RNA sequencing (RNA-seq) and proteomics data of 994 HCC patients with single-cell RNA-seq data of 11 samples across multiple datasets to identify an immune-related signature (IRS) that may influence sensitivity or resistance to such combined immunotherapy strategy, followed by verification of selected marker genes using immunohistochemistry and cytological experiments. We then trained/validated a cross-modality radiogenomics biomarker using machine learning based on TCIA database that was further tested in multi-scale independent cohorts covering 754 HCC patients. RESULTS: Integrative multi-omics analysis identifed a parsimonious 2-gene prognostic signature including KPNA2 and SMG5 that was significantly associated with immune heterogeneity and response to combination immunotherapy. Machine-learning pipeline exported the optimal 4-feature radiogenomics biomarker using support vector machine that significantly discriminated prognosis (hazard ratio 1.415–1.890; p < 0.05 for all) and modestly predicted response to ICI plus anti-angiogenic therapy (area under the curve 0.720–0.829) in independent retrospective series across major imaging modalities (computed tomography/magnetic resonance imaging). In a prospective neoadjuvant cohort, this biomarker also showed favorable performance for predicting pathological response and tumor recurrence, accompanied by biological validation through single-cell RNA-seq analysis of pre-treatment biopsies. CONCLUSIONS: Our study provides a cross-device-cross-modal radiogenomics biomarker that can improve patient selection for emerging ICI plus anti-angiogenic therapy with novel potential therapeutic targets in HCC.

Humans↗

PubMed2026

Transcriptome changes in circulating immune cells of critical COVID-19 patients predict a specific metabolic and epigenetic imprint.

BACKGROUND: The progression to critical COVID-19 arises predominantly from a dysregulated host immune response although the underlying regulatory mechanisms still remain partially elusive. This limits a prompt prediction of the disease progression, reduces the therapeutic options and restrains our understanding of “long COVID”. METHODS: Here, we analyzed the transcriptome of peripheral blood mononuclear cells (PBMCs) collected from COVID-19 patients experiencing different degrees of the disease (mild and critical), and control patients enrolled in the clinical trial COntAGIouS as well as independent bulk RNA-seq, single-cell RNA-seq and proteomic datasets. RESULTS: In critical COVID-19 patients, the integrative analysis of transcriptomic data revealed an altered regulatory network involving microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and coding genes that control mRNA translation-related genes, epigenetics, and metabolism. In parallel, we observed an upregulation of tRNA aminoacylation genes in critical COVID-19 patients by the analysis of either bulk or single-cell RNA-seq data from publicly available independent cohorts. Additionally, we found increased expression of coding genes enriched for the cognate amino acids (glycine, alanine, isoleucine and tyrosine), all related to protein localization, post-translational modifications, and cell metabolism in our cohort. Similar alterations in amino acid frequency were found in an independent proteomic dataset. CONCLUSIONS: Collectively, our findings indicate a broad perturbation of the gene expression landscape that characterizes the aberrant host immune response in critical COVID-19 patients and is potentially coordinated by miRNA and tRNA metabolism alterations. TRIAL REGISTRATION: COntAGIouS, NCT04327570. Registered 26 March 2020, https://clinicaltrials.gov/ct2/show/NCT04327570 .

Female↗

PubMed2026

Multimodal Analysis Reveals Aberrant Expression of SUMO2 and Its Significant Association With Key Mechanisms of Metabolic Pathways in Hepatocellular Carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related deaths worldwide. However, the role of small ubiquitin-like modifier 2 (SUMO2), a core member of the small ubiquitin-like modifier (SUMO) family, regarding its expression patterns and metabolism-related functions in HCC remains inadequately understood. METHODS: A multidimensional analytical framework was applied, integrating immunohistochemistry (153 HCC vs. 21 non-HCC samples), proteomics (159 paired samples), bulk transcriptomics (3240 HCC vs. 2267 non-HCC samples), single-cell RNA sequencing (RNA-seq) (10 HCC vs. 8 non-HCC samples), spatial transcriptomics, and external CRISPR/Cas9 functional genomics data. Systematic analyses included standardized mean difference (SMD), pathway enrichment, pseudotime trajectory inference, in silico knockout, cell-cell communication, metabolic flux scoring, immune infiltration, clinical correlation, drug sensitivity prediction, and molecular docking. RESULTS: At the protein level, immunohistochemistry (nuclear positivity) and external proteomic data collectively demonstrated consistent SUMO2 overexpression in HCC. Consistent upregulation was also observed at the mRNA level across large-scale cohorts. Single-cell RNA-seq and spatial transcriptomics localized SUMO2 enrichment to malignant hepatocytes and tumor-dominant regions. CRISPR-mediated SUMO2 knockout suppressed proliferation in multiple HCC cell lines. Mechanistically, high SUMO2 expression was significantly associated with metabolic reprogramming involving glycolysis/gluconeogenesis, pyruvate metabolism, and the tricarboxylic acid cycle. SUMO2-high malignant hepatocyte subpopulations exhibited enhanced activity of the macrophage migration inhibitory factor signaling axis and enhanced iron-sensor interactions. Further, the immune infiltration analysis revealed a negative correlation between SUMO2 expression and M1 macrophages and a positive correlation between follicular helper T cells and regulatory T cells. Clinically, elevated SUMO2 levels were found to be associated with adverse prognostic features. Furthermore, high SUMO2 expression was associated with increased sensitivity to dasatinib, and molecular docking simulations predicted potential binding between SUMO2 and dasatinib, with a Vina score of -8.5 kcal/mol. CONCLUSIONS: SUMO2 is aberrantly expressed at the protein, mRNA, single-cell, and spatial transcriptomic levels in HCC and is significantly associated with metabolic reprogramming and altered migration inhibitory factor (MIF)-mediated intercellular communication, suggesting its potential as a novel biomarker for diagnosis and treatment.

Humans↗

PubMed2026

Machine learning-integrated multi-omics risk prediction for pulmonary fungal infection in COPD and lung cancer: a transcriptomic and immune profiling study.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) and lung cancer are major risk factors for invasive pulmonary fungal infection (IPFI), carrying an attributable mortality of 30%-80%. Their coexistence further amplifies immunosuppression, while current diagnostic criteria remain inadequate for early risk identification. METHODS: Transcriptomic data from the GEO dataset GSE296912 (scRNA-seq; 12,078 cells from normal and COPD lung tissue) and The Cancer Genome Atlas (TCGA)-lung adenocarcinoma (LUAD) bulk RNA-seq cohort (539 tumor and 59 normal samples) underwent differential expression and cross-omics integration analysis. Five machine learning models were constructed: logistic regression, SVM, random forest, XGBoost, and LASSO. Candidate genes were validated by qRT-PCR in A549 cells and THP-1-derived macrophages stimulated with heat-inactivated Aspergillus fumigatus conidia, a protocol selected to ensure BSL-2 biosafety compliance and isolate PAMP-mediated innate immune signaling. Model performance was evaluated using 5-fold stratified cross-validation with AUC, calibration curves, and decision curve analysis. RESULTS: Single-cell transcriptomic analysis of 12,078 cells identified 14 distinct cell populations, with marked myeloid expansion and immune dysregulation in COPD lung tissue. Cross-omics integration with TCGA-LUAD data identified 1,145 shared genes (79 immune-related), converging on NF-κB, TLR4, and cytokine receptor signaling. The random forest model achieved excellent discriminative performance (5-fold CV AUC = 0.988), with Treg infiltration, TLR4, and MMP9 as the top predictors. qRT-PCR confirmed significant upregulation of all five candidate genes (DEFB4A, S100A8, IL-8, MMP9, and TLR4) in both A549 and THP-1 cells following fungal stimulation. CONCLUSION: This multi-omics machine learning model integrating scRNA-seq and TCGA transcriptomic data demonstrates excellent discriminative performance (AUC = 0.988), with mechanistic convergence of NF-κB, TLR4, and oncogenic signaling pathways identified across shared immune gene signatures. In vitro qRT-PCR validation confirms the biological relevance of five key antifungal immune genes, providing a transcriptomic foundation for future prospective IPFI risk stratification in patients with COPD and lung cancer.

TLR4↗

PubMed2026

Identifying Single-Cell Expression Quantitative Trait Loci Using a Bootstrap Penalized Hurdle Model.

BACKGROUND: Expression quantitative trait loci (eQTL) analysis links genetic variants to gene expression levels, helping to uncover how genetic variation contributes to gene regulation. While traditional eQTL analyses rely on bulk RNA-seq data, recent advances in single-cell RNA sequencing (scRNA-seq) have made it possible to detect cell-type-specific eQTLs. However, the inherent sparsity and heterogeneity of scRNA-seq data present major challenges for standard modeling approaches. METHODS: In this paper, we propose a novel statistical framework, Bootstrap Penalized Hurdle regression model (BPHurdle), designed specifically for scRNA-seq data. BPHurdle employs a hurdle modeling framework, where a logistic component accounts for the excess zeros in single-cell expression data, and a Poisson component jointly evaluates the effects of multiple SNPs on positive gene expression levels. RESULTS: Through simulation studies, we show that BPHurdle achieves high accuracy and robustness in identifying regulatory variants. We further demonstrate its utility on a real dataset through a case study focusing on a subset of differentially expressed genes, where it successfully identifies reliable cell-type-specific eQTLs. CONCLUSIONS: Overall, BPHurdle offers an advanced and flexible approach for single-cell eQTL mapping, providing deeper insight into the genetic regulation of gene expression at cellular resolution.

Quantitative Trait Loci↗

PubMed2026

Monocarboxylate Transporter 2 (MCT2) Reduction Is Associated with Increased Lung Tumor Growth and Alterations in the Immune Microenvironment in a Subcutaneous Tumor Model.

Monocarboxylate transporter 2 (MCT2; SLC16A7) is a high-affinity pyruvate transporter implicated in cancer metabolism. However, its role in lung cancer progression and the tumor microenvironment remains unclear. This study examined the effects of MCT2 reduction on tumor growth and cell-type-specific transcriptional changes within the tumor microenvironment. MCT2 loxP/loxP mice were crossed with mCre-Tg mice, and MCT2 deletion was induced by tamoxifen. Control (CO) mice received vehicle treatment. TC1 cells (100,000 cells/mouse) were injected subcutaneously, and tumors were harvested after 24 days. Single-nucleus RNA sequencing (snRNA-seq) was performed on isolated tumor nuclei (4000 nuclei/sample; n = 3 per group) using the 10x Genomics Chromium platform. Data were processed with Cell Ranger v3.0.2 and Seurat v5.2.1, followed by differential expression and pathway enrichment analyses integrated with macrophage bulk RNA-seq data. Tumors in mice with systemic MCT2 reduction grew significantly faster than those in control mice, demonstrating an association between host MCT2 reduction and increased tumor growth. Transcriptomic analysis generated high-quality profiles from 6864 CO and 10,055 KO nuclei. Clustering identified 12 cellular populations and cell types. MCT2 reduction altered pathways involved in glycolysis, the tricarboxylic acid cycle, oxidative phosphorylation, and fatty acid metabolism across multiple populations. Macrophages showed prominent transcriptional changes, including enrichment of MAPK, PI3K-Akt, IgSF-CAM, ECM, and cytokine-cytokine signaling pathways. These findings were supported by macrophage bulk RNA-seq data. Systemic MCT2 reduction was associated with increased tumor growth and broad transcriptional alterations within the tumor micro-environment. Differences in metabolic and immune-related transcriptional programs, particularly in macrophages, identify potential mechanisms associated with tumor progression that warrant further functional investigation.

Animals↗