SEARCH · Search PubMed

Results for “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.

At least 55 records · Page 3Linked to original sources

PubMed2024

Meningioma transcriptomic landscape demonstrates novel subtypes with regional associated biology and patient outcome.

Meningiomas, although mostly benign, can be recurrent and fatal. World Health Organization (WHO) grading of the tumor does not always identify high-risk meningioma, and better characterizations of their aggressive biology are needed. To approach this problem, we combined 13 bulk RNA sequencing (RNA-seq) datasets to create a dimension-reduced reference landscape of 1,298 meningiomas. The clinical and genomic metadata effectively correlated with landscape regions, which led to the identification of meningioma subtypes with specific biological signatures. The time to recurrence also correlated with the map location. Further, we developed an algorithm that maps new patients onto this landscape, where the nearest neighbors predict outcome. This study highlights the utility of combining bulk transcriptomic datasets to visualize the complexity of tumor populations. Further, we provide an interactive tool for understanding the disease and predicting patient outcomes. This resource is accessible via the online tool Oncoscape, where the scientific community can explore the meningioma landscape.

Meningioma↗

PubMed2026

A novel deep learning-driven framework for improving lncRNA comprehensive annotation with LncADeep 2.0.

MOTIVATION: Long non-coding RNAs (lncRNAs) have emerged as crucial players in diverse physiological and pathological processes, yet the biological mechanisms of the vast majority of lncRNAs remain elusive. To fill this gap, it is necessary to improve the accuracy of lncRNA identification and functional annotation. RESULTS: Here, we introduce LncADeep 2.0, an integrated deep learning framework designed to meet these needs. In the identification module, LncADeep 2.0 incorporated novel peptide features along with sequence and structural information, demonstrating superior performance over our previous LncADeep and other existing tools on both annotated transcripts from GENCODE and RNA-seq data. For functional annotation, LncADeep 2.0 leveraged lncRNA-centric interaction networks and gene ontology terms through the transfer learning strategy to achieve robust annotation performance with limited functional data. Compared to LncADeep, LncADeep 2.0 could accurately elucidate the general functions of given lncRNA sequences, predict tissue- or cell-type-specific functions from bulk and single-cell RNA-seq data, and establish connections between tumor-associated lncRNAs and genomic markers. Overall, LncADeep 2.0 stands out as an efficient and reliable tool for lncRNA identification and functional annotation across a wide spectrum of biological processes. AVAILABILITY AND IMPLEMENTATION: LncADeep 2.0 is available for use at https://github.com/Jefferson-Chou/LncADeep2 and https://doi.org/10.5281/zenodo.17164767.

RNA, Long Noncoding↗

PubMed2025

The offonome reveals on and off states of gene expression near the detection limit of RNA-seq.

RNA-seq, widely used for gene expression profiling, provides nucleotide level genome coverage and summary gene expression values. Generally, low-expressed genes are ignored due to their unfavorable signal-to-noise ratio, however, these genes may offer crucial information, such as detecting rare cells in bulk tissues. In this study, we applied an approach that transforms the expression levels of low-expressed genes into a robust dichotomized on/off state by leveraging similarities in transcript coverage shape. Applied to three human cancer cohorts from the Cancer Genome Atlas (TCGA), chosen based on tissue morphology and anatomic site, we identified genes, the "offonome" near the detection limit, consistently or occasionally off across samples. Genes in the offonome spectrum proved useful for supervised and unsupervised applications, including characterizing oncogenic pathways, and identifying rare populations of cells in bulk tissue. Interrogating the offonome is relevant to bulk tumor analyses like TCGA, potentially expediting gene investigation in low-input situations like single cell RNA-seq.

Humans↗

PubMed2026

A transcription factor regulatory atlas for activity inference and perturbation prediction.

Inferring transcription factor (TF) activity from transcriptomes and predicting transcriptome-wide responses to TF perturbations remain challenging, in part because available TF-mRNA resources often face a trade-off between precision and coverage and typically lack signed regulatory information. Here, we present TFActProfiler, a TF-mRNA resource and computational framework that learns signed, quantitative TF-mRNA regulatory coefficients by integrating heterogeneous prior evidence (ChIP-based, motif-based, and curated TF-mRNA annotations) with large-scale bulk and single-cell RNA-seq atlases. TFActProfiler contains 2 606 176 signed TF-mRNA interactions and improves TF activity inference in TF knockdown benchmarks relative to widely used regulon resources while retaining broad TF and target coverage. In addition, because the same learned regulatory coefficients can be used to model downstream transcriptional effects, TFActProfiler enables prediction of transcriptome-wide gene expression responses to TF knockdown without training on task-matched perturbation data. When perturbation datasets are available, TFActProfiler can be further refined to achieve performance comparable to state-of-the-art machine-learning baselines. By providing a direction-aware representation of TF-mRNA regulation for both activity inference and perturbation-response modeling, TFActProfiler supports systematic dissection of gene regulatory programs across diverse cellular contexts.

Transcription Factors↗

PubMed2026

Activation of transposable elements is linked to a region- and cell type-specific interferon response in Parkinson's disease.

Parkinson's disease (PD) is a neurodegenerative disorder involving a neuroinflammatory response, the cause of which remains unclear. Transposable elements (TEs) have been linked to inflammation, but their potential role in PD remains unexplored. Using bulk- and single-nuclei RNA-seq of postmortem brain tissue from four brain regions, we studied TE transcription and its correlation with PD neuroinflammation. Over a thousand TEs, including LINE-1 s and ERVs, were expressed in a cell type- and region-specific manner in the human brain. Increased TE expression was found in microglia and neurons in the substantia nigra and putamen of PD brains, but not amygdala or prefrontal cortex, compared to controls. This TE activation correlated with an innate immune response in the same brain regions. The link between an interferon response and TE activation was mechanistically confirmed using human pluripotent stem cell-derived microglia and neurons. Our findings provide insights into TE transcription in the PD brain and suggest that TEs may contribute to neuroinflammation and pathological progression in PD.

Humans↗

PubMed2025

Comprehensive Transcriptome Annotation of Thousands of HIV-1 Genomes.

Alternative splicing in HIV-1 has been a central focus of decades of research, uncovering key mechanisms of viral gene regulation, immune evasion, and therapeutic response - yet, no reference resource has existed to support transcriptome-wide analysis, limiting adoption of modern computational methods. We present HIV Atlas (https://ccb.jhu.edu/HIV_Atlas), the first reference-quality annotation of HIV-1 and SIV transcriptional diversity. We manually curated transcriptomes for HIV-1HXB2 and SIVmac239 and developed Vira, an automated annotation-transfer method specifically designed to address unique challenges of viral genome biology, to generate high-quality annotations for 2,077 complete HIV-1 genomes. Using the resources presented in our work, we evaluated conservation of splice sites, revealing near-perfect preservation of major donors and acceptors. Furthermore, using several public datasets, we demonstrate how HIV Atlas enhances methodology, improves the quality and novelty of results, and opens novel avenues for research, supporting more accurate and comprehensive analyses of bulk, single-cell, and spatial RNA-seq in HIV-1 studies.

Journal Article↗

PubMed2026

Integrative analysis of single-cell sequencing identifies CD8+ TIM3+ CD101+ T cell-associated genes as prognostic biomarkers in breast cancer.

BACKGROUND: Breast cancer is a prevalent and deadly malignancy that significantly impacts women's quality of life and imposes financial burdens. Despite therapeutic advancements, tumour heterogeneity and frequent relapses remain major challenges. Accordingly, this study aimed to characterize immune features associated with CD8+ TIM3+ CD101+ T cells and develop a prognostic signature for breast cancer. METHODS: This study integrated single-cell and bulk transcriptomic datasets to characterize CD8+ TIM3+ CD101+ T cell (CCT)-related immune features and construct a prognostic signature in breast cancer. Single-cell RNA-seq data were sourced from the Gene Expression Omnibus (GEO) repository, and bulk transcriptomic data were from The Cancer Genome Atlas (TCGA) and GEO databases. Analytical methods included pseudo-time trajectory reconstruction (Monocle2), intercellular signalling analysis (CellChat), functional enrichment (ClusterProfiler), and immune profiling (ssGSEA). Prognostic modeling was conducted using least absolute shrinkage and selection operator (LASSO) Cox regression, with validation via Kaplan-Meier and time-dependent receiver operating characteristic (ROC) analyses. RESULTS: Single-cell analysis identified 17 clusters spanning seven cell types, including T cells, myeloid cells, and epithelial cells. T-cell sub-clustering revealed four subtypes. Pseudotime analysis suggested a potential state-transition relationship between CD8+ CD101- TIM3+ and CD8+ CD101+ TIM3+ T-cell states. A total of 121 differentially expressed genes were enriched in vital biological processes. An 11-gene prognostic model showed strong predictive power across cohorts. Single-cell T-cell reclustering identified a CD8+ CD101+ TIM3+ T-cell subpopulation, which was primarily characterized by the expression of markers such as CD101 and HAVCR2/TIM3. CONCLUSIONS: This study maps cellular heterogeneity and molecular networks in breast cancer, offering insights for targeted therapy and improved prognosis.

Breast invasive carcinoma↗

PubMed2026

Integration of Genome-Wide Association Studies With Single-Cell and Bulk Expression Quantitative Trait Locus to Identify Stroke Susceptibility Genes.

BACKGROUND: Previous studies have integrated genome-wide association studies with expression quantitative trait locus (eQTL) data from bulk tissues to identify stroke susceptibility genes. However, eQTL data exhibit high cell-type specificity, and genetic variants may have distinct effects across stroke subtypes. METHODS: We applied the summary-data-based Mendelian randomization (MR) method to integrate eQTL data from 7 brain cell types with genome-wide association studies data for 5 stroke phenotypes (stroke, ischemic stroke, cardioembolic stroke, large artery stroke, and small vessel stroke). Results were compared with summary-data-based MR using eQTL data from 49 tissues in the Genotype-Tissue Expression project. Robustness of significant single-cell summary-data-based MR associations was assessed via MR and colocalization analyses. Further evaluations included single-cell RNA-seq differential expression, protein-protein interaction, druggability, and phenome-wide association studies. RESULTS: Single-cell summary-data-based MR identified many novel significant genes not detected using bulk tissue eQTL data. Validated associations revealed 2 stroke risk genes (LRCH1, ICA1L), 3 stroke protective genes (AHI1, LYRM9, CENPQ), 2 large artery stroke risk genes (LIPA, ELL), and 1 ischemic stroke protective gene (CENPQ). Single-cell RNA-seq showed significantly increased LIPA expression in mouse stroke samples compared with controls. Protein-protein interaction and druggability analyses, along with phenome-wide association studies, prioritized LIPA and LRCH1 as potential therapeutic targets for stroke while indicating possible adverse effects. CONCLUSIONS: Integrating single-cell eQTL with stroke-subtype genome-wide association studies uncovers novel cell-type-specific causal genes and highlights promising therapeutic targets, advancing understanding of stroke pathogenesis.

Genome-Wide Association Study↗

PubMed2025

Twisted Sister1: an agravitropic mutant of bread wheat (Triticum aestivum) with altered root and shoot architectures.

We identified a mutant of hexaploid wheat (Triticum aestivum) with impaired responses to gravity. The mutant, named Twisted Sister1 (TS1), had agravitropic roots that were often twisted along with altered shoot phenotypes. Roots of TS1 were insensitive to externally applied auxin, with the genetics and physiology suggestive of a mutated AUX/IAA transcription factor gene. Hexaploid wheat possesses over 80 AUX/IAA genes, and sequence information did not identify an obvious candidate. Bulked segregant analysis of an F2 population mapped the mutation to chromosome 5A, and subsequent mapping located the mutation to a 41 Mbp region. RNA-seq identified the TraesCS5A03G0149800 gene encoding a TaAUX/IAA protein to be mutated in the highly conserved domain II motif. We confirmed TraesCS5A03G0149800 as underlying the mutant phenotype by generating transgenic Arabidopsis thaliana. Analysis of RNA-seq data suggested broad similarities between Arabidopsis and wheat for the role of AUX/IAA genes in gravity responses, although there were marked differences. Here we show that the sequenced wheat genome, along with previous knowledge of the physiology of gravity responses from other plant species, gene mapping, RNA-seq, and expression in Arabidopsis have enabled the cloning of a key wheat gene that defines plant architecture.

Triticum↗

PubMed2026

Single-cell transcriptional profiling identifies the swimming crab Portunus trituberculatus in response to bacterial infection.

Crustaceans rely entirely on innate immunity, yet the cellular composition, functional specialization, and pathogen-induced remodeling of their immune system remain poorly resolved. Here, we generated a high-resolution single-cell transcriptomic atlas of hemocytes from the swimming crab Portunus trituberculatus following Vibrio parahaemolyticus infection using 10× Genomics scRNA-seq. Seven putatively distinct hemocyte clusters were identified, including granulocytes, semigranular hemocytes, prohemocytes, unresolved hemocytes, hyalinocyte-like hemocytes, biosynthetically active secretory hemocytes, and regulatory hemocytes. Although the overall cellular composition remained relatively stable after infection, hemocytes exhibited pronounced cluster-specific transcriptional reprogramming involving Toll/NF-κB signaling, antimicrobial peptide synthesis and metabolic rewiring. By integrating single-cell and bulk transcriptomes, we identified multiple anti-lipopolysaccharide factors (ALFs) as key secretory effectors and experimentally validated their antibacterial activities. FITC-based bacterial engulfment assays and RNA-seq of sorted phagocytes demonstrated that phagocytic capability was shared across multiple hemocyte clusters. Notably, the immunoglobulin superfamily receptor DSCAM displayed extensive alternative splicing and strong infection-induced activation in unresolved hemocytes. Immune-training experiments showed that prior bacterial exposure was associated with altered DSCAM expression and reduced early cumulative mortality upon secondary challenge, suggesting a memory-like immune phenotype. These findings provide a foundational framework for understanding crustacean immunity and advancing disease-resistant breeding in aquaculture.

Antimicrobial peptides↗

PubMed2026

FANCI promotes esophageal squamous cell carcinoma progression and cell cycle regulation and interacts with FANCD2.

BACKGROUND: Esophageal squamous cell carcinoma (ESCC) is an aggressive malignancy with poor clinical outcomes, and reliable molecular biomarkers and therapeutic targets remain limited. Fanconi anemia group I protein (FANCI) is a core component of the Fanconi anemia (FA) pathway, but its expression pattern, clinical significance, and functional role in ESCC have not been comprehensively defined. This study aimed to investigate FANCI expression and prognostic value in ESCC, assess its effects on malignant cellular phenotypes and tumor growth, and explore its potential mechanistic relationship with Fanconi anemia group D2 protein (FANCD2) and cell-cycle regulation. METHODS: Multi-cohort analyses were performed using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets, together with ESCC single-cell RNA sequencing (RNA-seq) data. FANCI functions were assessed by bidirectional gain- and loss-of-function experiments in vitro (proliferation, colony formation, migration, invasion, apoptosis, and cell-cycle assays) and by xenograft models in vivo. Mechanistic studies included protein-protein interaction (PPI) analyses, co-immunoprecipitation (Co-IP), and immunofluorescence (IF) colocalization. RESULTS: FANCI was consistently upregulated in ESCC across bulk transcriptomic datasets and was further supported by quantitative polymerase chain reaction (qPCR), Western blotting, and immunohistochemistry (IHC). FANCI discriminated ESCC from normal tissues in TCGA-ESCC and was independently validated in GSE53624 [area under the curve (AUC) =0.940 and 0.975, respectively]. FANCI was associated with poorer overall survival (OS) and shorter disease-free interval (DFI), and these findings were validated in an independent GEO cohort. Functionally, FANCI promoted ESCC cell proliferation, migration, and invasion, while inhibiting apoptosis; FANCI knockdown suppressed tumor growth in vivo and induced G2/M cell-cycle arrest. Mechanistically, FANCI physically interacted with FANCD2, colocalized with FANCD2 in the nucleus, and was associated with altered FANCD2 protein abundance, consistent with cell-cycle and DNA repair-related programs. Single-cell analysis indicated that FANCI was enriched in epithelial cells and associated with higher activity of malignant functional programs. In TCGA-ESCC, FANCI-high tumors showed distinct mutation profiles, a trend toward increased tumor mutation burden (TMB), and altered immune-associated signatures. CONCLUSIONS: FANCI is upregulated in ESCC and is associated with diagnostic and prognostic value. It promotes malignant phenotypes and tumor growth, potentially through a FANCI-FANCD2-linked cell-cycle/DNA repair program, supporting FANCI as a candidate biomarker and therapeutic target in ESCC.

Esophageal squamous cell carcinoma (ESCC)↗

PubMed2026

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase↗

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

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

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

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↗