SEARCH · PubMed Health

Results for “transcriptomic profiling”

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

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

At least 487 records · Page 27Linked to original sources

PubMed2026

Comparative transcriptomic analysis of the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress.

To investigate the differences in molecular responses between the gills and hepatopancreas of freshwater-cultured Litopenaeus vannamei under chronic nitrite stress, a 30-day chronic stress experiment was conducted with a control group and a stress group. Transcriptomic analysis of the gills and hepatopancreas was performed using Illumina sequencing; differentially expressed genes (DEGs) were identified, and GO, KEGG, GSEA, PPI, and RT-qPCR validation were carried out. The results showed that 196 DEGs (161 up-regulated and 35 down-regulated) were identified in the gills, and 287 DEGs (199 up-regulated and 88 down-regulated) in the hepatopancreas, with only 18 DEGs shared between the two tissues. DEGs in the gills were enriched in oxidoreductase activity, glycerophospholipid metabolism, and tyrosine metabolism; DEGs in the hepatopancreas were enriched in lipid transporter activity, phagosome, ECM-receptor interaction, and riboflavin metabolism. GSEA revealed significant suppression of the mTOR pathway in the gills and the Polycomb complex pathway in the hepatopancreas. PPI network analysis identified hub genes P5CS and eEF2 in the gills, and PER, TUBB1, SHMT, and TUBB4B in the hepatopancreas. RT-qPCR validation was consistent with the RNA-seq results (R2 = 0.764). This study indicates that, under chronic nitrite stress, the gill response is centered on redox regulation and inhibition of growth metabolism, whereas the hepatopancreas response primarily involves lipid transport, cytoskeletal remodeling, and phagosome activation. The two tissues synergistically adapt through fundamental biosynthetic and motor protein pathways. This research provides molecular evidence for deciphering the nitrite tolerance mechanisms in freshwater-cultured shrimp.

Animals↗

PubMed2026

Identification of Critical Genes for Recurrent Aphthous Ulcer by Transcriptome Data Analysis and Mendelian Randomization.

PURPOSE: Recurrent aphthous ulcer (RAU) is a common oral mucosal disorder with a poorly understood etiology, significantly affecting patients' quality of life. This study aims to investigate critical genes linked to RAU and explore their biological mechanisms using transcriptomic data and Mendelian randomization (MR) analysis. MATERIALS AND METHODS: RAU-related gene expression data from the GEO database (GSE37265) were analyzed to identify differentially expressed genes (DEGs). A two-sample MR approach was used to assess the causal impact of expression quantitative trait loci (eQTL) on RAU. Critical genes were identified by intersecting DEGs with significant MR findings. GO and KEGG pathway enrichment analyses were performed, along with GSEA and immune cell infiltration analysis, to investigate the functions and mechanisms of these genes in RAU. RESULTS: A total of 184 differentially expressed genes (DEGs) were identified, while 339 RAU-associated genes were screened through MR analysis. Cross-validation further identified 7 critical genes. Among these, CCR1, ERP27, HCK, MICB, and SLC2A3 showed protective associations with RAU risk, whereas CD177 and IFITM1 were positively associated with increased risk. Enrichment analysis revealed that these genes are involved in specific biological processes, including cell migration, immune response, and metabolic regulation, which are closely linked to RAU pathogenesis. CONCLUSION: This systematic study comprehensively investigates the critical causative genes underlying RAU, emphasizing the intricate relationships between immune regulation and metabolic disturbances in its pathology. These findings lay a solid foundation for the development of novel biomarkers and may inform future research on targeted therapeutic strategies for RAU.

Stomatitis, Aphthous↗

PubMed2025

Inflammatory pathways and immune dysregulation in pediatric postoperative septic shock: A study integrating transcriptomics, machine learning and molecular docking.

This study elucidates the molecular and immune regulatory mechanisms of pediatric postoperative septic shock. Transcriptomic data were obtained from the Gene Expression Omnibus database. Differentially expressed genes were identified using the limma package, and gene co-expression modules were constructed using Weighted Gene Co-expression Network Analysis. Functional enrichment was performed via gene set enrichment analysis, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes analyses. Immune cell infiltration was assessed using ESTIMATE and CIBERSORT. Mendelian randomization was applied to explore causal relationships between gene expression and septic shock. Feature genes were selected using machine learning algorithms, and a diagnostic nomogram model was constructed. Finally, molecular docking analysis was performed to screen and evaluate the binding affinity of traditional Chinese medicine monomers to core target proteins. A total of 1331 differentially expressed genes were identified, and the turquoise module was strongly correlated with septic shock. Enrichment analysis revealed significant activation of IL-6/JAK/STAT3, TNF-α/NF-κB, and PI3K/Akt/mTOR pathways. Immune infiltration analysis indicated suppressed immune scores and imbalances in neutrophils, macrophages, T cells, and B cells. Mendelian randomization confirmed causal associations for 6 genes, including PIM3. The predictive model based on feature genes demonstrated high diagnostic performance. Molecular docking suggested that quercetin and astramembrannin I could stably bind PIM3. This study systematically identified core genes, dysregulated immune pathways, and candidate small-molecule interventions in pediatric septic shock, providing novel insights for early diagnosis and targeted therapy.

Humans↗

PubMed2026

Integrative machine learning and transcriptomic analysis reveals molecular mechanisms underlying low survival rate in larval Chinese Bahaba (Bahaba taipingensis).

Chinese Bahaba (Bahaba taipingensis) is a Class I protected marine fish endemic to China. Low larvae survival during artificial breeding severely hinder population recovery. To investigate the molecular mechanism of high mortality in larval fish, this study performed RNA-seq on liver from naturally deceased (ND) and mass-dead (MD) individuals, combined with least absolute shrinkage and selection operator (LASSO) regression and random forest (RF) algorithms to screen for core signature genes. A total of 873 differentially expressed genes (DEGs) were identified, including 112 upregulated and 761 downregulated genes. GO and KEGG enrichment analyses revealed significant enrichment in amino acid metabolism disorders, one‑carbon folate pool impairment, PPAR signaling abnormalities, ECM-receptor interaction, focal adhesion pathway, indicating widespread metabolic suppression accompanied by extracellular matrix remodeling and signaling disturbances in the livers of MD fish. MAD pre-filtering combined with dual machine learning algorithms yielded 18 robust core signature genes, among which SLC38A4, MMP1, FADD, FKBP5, and APOB were consistently identified as high-frequency core genes by both algorithms. SLC38A4 exhibited the highest importance score in the RF model and was significantly downregulated, making it the primary molecule distinguishing ND from MD phenotypes. ROC curve analysis showed that both models achieved an AUC of 1.000 (95% CI lower bound: 0.610), confirming the precise discriminatory ability of the core genes. GSEA further demonstrated significant enrichment of this core gene set in ND samples. This study provides the first systematic elucidation of the molecular mechanisms underlying liver dysfunction in low survival rate B. taipingensis, characterized by amino acid transport impairment, metabolic reprogramming, and structural remodeling, offering theoretical foundations for health assessment, early mortality risk warning, and artificial breeding conservation of this species.

Animals↗

PubMed2025

Transcriptome-wide analysis reveals potential roles of CFD and ANGPTL4 in fibroblasts regulating B cell lineage for extracellular matrix-driven clustering and novel avenues for immunotherapy in breast cancer.

BACKGROUND: The remodeling of the extracellular matrix (ECM) plays a pivotal role in tumor progression and drug resistance. However, the compositional patterns of ECM in breast cancer and their underlying biological functions remain elusive. METHODS: Transcriptome and genome data of breast cancer patients from TCGA database was downloaded. Patients were classified into different clusters by using non-negative matrix factorization (NMF) based on signatures of ECM components and regulators. Weighted Gene Co-expression Network Analysis (WGCNA) was used to identify core genes related to ECM clusters. Additional 10 independent public cohorts including Metabric, SCAN_B, GSE12276, GSE16446, GSE19615, GSE20685, GSE21653, GSE58644, GSE58812, and GSE88770 were collected to construct Training or Testing cohort, following machine learning calculating ECM correlated index (ECI) for survival analysis. Pathway enrichment and correlation analysis were used to explore the relationship among ECM clusters, ECI and TME. Single-cell transcriptome data from GSE161529 was processed for uncovering the differences among ECM clusters. RESULTS: Using NMF, we identified three ECM clusters in the TCGA database: C1 (Neuron), C2 (ECM), and C3 (Immune). Subsequently, WGCNA was employed to pinpoint cluster-specific genes and develop a prognostic model. This model demonstrated robust predictive power for breast cancer patient survival in both the Training cohort (n = 5,392, AUC = 0.861) and the Testing cohort (n = 1,344, AUC = 0.711). Upon analyzing the tumor microenvironment (TME), we discovered that fibroblasts and B cell lineage were the core cell types associated with the ECM cluster phenotypes. Single-cell RNA sequencing data further revealed that angiopoietin like 4 (ANGPTL4)+ fibroblasts were specifically linked to the C2 phenotype, while complement factor D (CFD)+ fibroblasts characterized the other ECM clusters. CellChat analysis indicated that ANGPTL4+ and CFD+ fibroblasts regulate B cell lineage via distinct signaling pathways. Additionally, analysis using the Kaplan-Meier Plotter website showed that CFD was favorable for immunotherapy response, whereas ANGPTL4 negatively impacted the outcomes of cancer patients receiving immunotherapy. CONCLUSION: We identified distinct ECM clusters in breast cancer patients, irrespective of molecular subtypes. Additionally, we constructed an effective prognostic model based on these ECM clusters and recognized ANGPTL4+ and CFD+ fibroblasts as potential biomarkers for immunotherapy in breast cancer.

Humans↗

PubMed2025

Colorectal Liver Metastasis Pathomics Model: Integrating Single-Cell and Spatial Transcriptome Analysis With Pathomics for Predicting Liver Metastasis in Colorectal Cancer.

The liver is the primary target organ for hematologic metastasis of colorectal cancer (CRC), and CRC liver metastasis (CRLM) often precludes radical resection, making it the leading cause of death in patients with CRC. To improve the identification and prediction of liver metastasis risk, we identified a cell type of liver metastasis--triggering malignant cells (LMTMCs) through integrating single-cell RNA sequencing and spatial transcriptome analysis. Multiomics cell communication analysis indicated that the interaction between fibroblasts and LMTMCs through the COL1A1-CD44/SDC4 and LAMA4-CD44 signaling axes could promote CRLM. By applying the one-class logistic regression algorithm, we developed a CRLM scoring system in the bulk RNA-sequencing data according to the abundance of LMTMCs in each individual. Using the grouping labels derived from the CRLM scoring system in the bulk data and the corresponding whole-slide images without any manual annotations at the region or pixel level, processed via slide-level weakly supervised learning, a deep-learning model based on the ResNet18 architecture, called Colorectal Liver Metastasis Pathomics Model, was developed to predict the risk of liver metastasis in patients with CRC. The Colorectal Liver Metastasis Pathomics Model achieved an area under the curve of 0.84 at the internal test set of The Cancer Genome Atlas-CRC histology images. In the external independent validation sets, namely the Affiliated Hospital of Southwest Medical University and the Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University cohorts, the areas under the curve were 0.89 and 0.72, respectively, indicating effective classification performances. This study provided new insights and tools for the early identification of CRLM and demonstrated the potential of combining multiomics with deep learning-based pathomics in cancer research.

Humans↗

PubMed2026

Spatial transcriptomics of primary and metastatic ALK-rearranged NSCLC reveals site-specific adaptations.

INTRODUCTION: Genetic alterations and the tumor microenvironment (TME) influence treatment response in anaplastic lymphoma kinase-rearranged non-small cell lung cancer (ALK+ NSCLC). This study maps site-specific TME adaptations and exploratory risk-associated signatures in lymph node metastases (LNT) to investigate metastatic evolution. METHOD: We applied spatial transcriptomics to profile tumor (PanCK+) and stromal (PanCK-) compartments in a pilot cohort of 16 cases: primary lung tumors (LT, n = 3), LNT (n = 10), and brain metastases (BT, n = 3), with three site-matched non-tumor controls. LNT-derived prognostic signatures were evaluated using The Cancer Genome Atlas-Lung Adenocarcinoma (TCGA LUAD) cohorts. RESULTS: Distinct, site-specific TME features were observed. LNT stroma was enriched in fibroblasts and macrophages, while tumor segments showed increased neutrophils. BT exhibited a macrophage-associated immunosuppressive TME. Tumor cells evolved divergently: LT retained pulmonary identity and showed trend towards translation-associated programs, LNT cells shifted toward senescence and epigenetic remodeling, and BT cells showed activation of Class A/1 (Rhodopsin-like) receptor, GPCR and drug metabolism pathways. In LNT, exploratory risk-associated differences were observed. Low-risk cases (n = 6) showed adaptive immune signatures, whereas high-risk cases (n = 4) showed enrichment for stromal MET signaling and stress-response pathways. Because treatment exposure differed markedly between the risk groups, these observations should be interpreted as hypothesis-generating. TCGA LUAD analysis suggested the broader biological relevance of immune-associated markers, but reflected general LUAD rather than ALK+ specific biology. Discordant associations for GCLC and TIMP1 underscored the importance of spatial context. CONCLUSION: Site-specific microenvironments may influence tumor adaptation across metastatic niches in ALK+ NSCLC. The exploratory risk-associated findings require validation in larger, uniformly treated cohorts.

Humans↗

PubMed2026

An oxidative stress - and immunotherapy-related six-gene signature defines immune subtypes and predicts prognosis and immunotherapy response in hepatocellular carcinoma.

BACKGROUND: Oxidative stress and the tumor immune microenvironment jointly shape hepatocellular carcinoma (HCC) progression and response to immunotherapy, yet integrated biomarkers linking these processes are lacking. METHODS: Transcriptomic and clinical data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were used to identify oxidative stress- and immunotherapyrelated differentially expressed genes (OSIRDEGs). Functional enrichment, weighted gene co-expression network analysis (WGCNA) and LASSO-Cox regression were used to construct a prognostic signature. Consensus clustering, TIDE, CIBERSORT and ssGSEA characterized immune phenotypes. Somatic mutation, copy-number and drug-response data were integrated to assess genomic alterations and drug sensitivity. Expression of model genes was validated by qRT-PCR and western blotting in HCC cell lines. RESULTS: We identified 24 OSIRDEGs enriched in cell-cycle and mitotic pathways. WGCNA intersection yielded 18 module genes, from which a six-gene signature (BUB1B, CDKN2A, CENPE, HMMR, PTTG1, SPP1) was derived. The signature robustly stratified patients into high- and low-risk groups with significantly different progression-free and disease-free survival in both TCGA-LIHC and GSE14520. Based on signature expression, two molecular subtypes were defined, exhibiting distinct survival, immune landscapes and predicted immunotherapy responsiveness. Model genes harbored recurrent alterations and showed significant correlations with anticancer agents. All six genes were upregulated at mRNA and protein levels in metastatic HCC cell lines versus normal hepatocytes. CONCLUSIONS: We systematically explored the landscape of OSIRDEGs in HCC, and proposed a validated six-gene signature that refines prognostic stratification, delineates immunerelevant HCC subtypes and highlights candidate biomarkers for therapeutic selection and mechanistic investigation.

Humans↗

PubMed2026

Mineral-driven molecular signatures of energy metabolism underpin sperm motility in buffalo.

The success of spermatogenesis depends on the interplay of various biomolecules that ultimately determine sperm quality. In this study, RNA-seq analysis of frozen-thawed buffalo sperm (n = 8) revealed the presence of 263 mineral-associated genes (>1 FPKM) in high (n = 4) and 181 in low motile (n = 4) sperm groups. Among these, 177 mineral-associated genes were commonly expressed between them, and the majority were upregulated (>1 fold), LOC102391588 (ncRNA; 37-fold), ZNF699 (26.5-fold), MYZAP (13-fold), etc., in the high motile group. The expression of selected mineral-associated genes was validated. The top enriched functions in commonly expressed genes were regulation of transcription by RNA polymerase II (FDR: 4.7 × 10⁻2; ZNF331, ZNF692, ZNF180, etc.), followed by spermatogenesis (FDR: 2.9 × 10⁻2; CALR3, ADAM18, ADAM29, etc.), proton transmembrane transport (FDR: 4.0 × 10⁻2; ATP6V0E1, ATP1A4, ATP6V1B2, etc.) and flagellated sperm motility (FDR: 2.3 × 10⁻1; CATSPERD, EFCAB6, CABS1, etc.). Additionally, the chromatin remodeling pathway (FDR: 3 × 10⁻2; PTP4A1, PPM1A, DUSP1, etc.) emerged as the most significant and may suggest that these minerals influence genome packaging and sperm functionality. Mineral-associated genes were predominantly associated with zinc (49%), followed by calcium (20%), phosphorus (10%), iron (5%), sodium (2%), potassium (2%), copper (1%) and other trace elements (11%). Although the current study uses frozen-thawed sperm, the findings indicate that mineral-associated genes are crucial for promoting membrane stability, energy production, motility and chromatin integrity, which may contribute to the superior fertilizing ability of sperm.

Animals↗

PubMed2026

Correlative analysis of endogenous miRNA expression profiles underlying brown planthopper adaptation to resistant rice.

The brown planthopper (Nilaparvata lugens Stål, BPH) is a major insect pest threatening global rice production. However, the molecular mechanisms underlying the adaptation of BPH populations with different virulence levels to resistant rice cultivars remain poorly understood. MicroRNAs (miRNAs), as key post-transcriptional regulators, play critical roles in host adaptation in herbivorous insects. In this study, we analyzed the miRNA expression profiles of a high-virulent population (IR56p) and a low-virulence population (TN1p) after feeding on susceptible (TN1) and resistant (IR56) rice cultivars. Our findings reveal distinct miRNA-mediated regulatory strategies employed by the two populations. The IR56p population showed downregulation of miRNAs including miR-10, miR-124, and miR-316, showing an inverse correlation with increased expression of predicted target genes involved in detoxification (carboxylesterase, UDP-glycosyltransferase) and effector function (calmodulin). In contrast, several miRNAs highly expressed in IR56p, including miR-307, miR-317, and miR-275, were predicted to target rice genes associated with hormone signaling, cell wall biosynthesis, and oxidative homeostasis, suggesting a possible but unproven inter-species regulatory role that requires functional validation. Collectively, these descriptive and correlative findings provide hypothesis generating insights into insect-plant coevolution and identifies candidate molecular targets for future functional validation and RNA interference-based pest management strategies.

Animals↗

PubMed2026

Meta-PseU: A meta-classifier for robust prediction of RNA pseudouridine modification sites from long sequences.

BACKGROUND AND OBJECTIVES: Pseudouridine (Ψ) represents one of the most abundant and conserved RNA modifications. Ψ provides an additional hydrogen-bond donor that enhances RNA structural stability and modulates translation. It participates in diverse biological processes, including RNA-protein interactions, splicing, translational control, and stress responses. Aberrant pseudouridylation is implicated in cancer, neurodegenerative disorders, and autoimmune diseases. Despite its biological importance, experimental identification of Ψ sites remains time-consuming and costly, limiting the feasibility of transcriptome-wide profiling. Computational approaches have therefore become essential complements to experimental techniques. However, state-of-the-art machine-learning and deep-learning predictors often suffer from limited generalizability due to small training datasets. To overcome these issues, we aim at constructing new long-sequence datasets and developing a novel Ψ site predictor. METHODS: New long-sequence datasets were constructed as benchmarks for RNA Ψ-site prediction. The Ψ modification sites in RMBase 3.0 were mapped to the reference genomes across three species of human, mouse, and yeast, and the RNA sequences with a length of 201 were generated by extending the upstream and downstream from the mapped, central sites. To eliminate sequence redundancy, the sequences were clustered using CD-HIT with a 70% sequence identity threshold. We developed Meta-PseU, a logistic regression-based meta-classifier that considered 118 machine learning and deep learning classifiers. The datasets and programs are freely accessible at https://github.com/kuratahiroyuki/MetaPseU. RESULTS: By optimizing model configuration, we proposed the Meta-PseU model stacking 32 machine learning and deep learning classifiers out of 118 classifiers. Meta-PseU substantially improved model generalizability, overcoming a key limitation of existing approaches. It greatly outperformed state-of-the-art predictors and achieved increasing accuracy with increasing sequence length. CONCLUSIONS: Long-sequence datasets were newly constructed as benchmarks for RNA Ψ-site prediction. Meta-PseU offers a new framework for robust Ψ-site identification by using long sequences.

Pseudouridine↗

PubMed2004

Frontiers in nutrigenomics, proteomics, metabolomics and cancer prevention.

While dietary habits continue to surface as a significant factor that may influence cancer incidence and tumor behavior, there is considerable scientific uncertainty about who will benefit most. Adequate [corrected] knowledge about how the responses depend on an individual's genetic background (nutrigenetic effects), the cumulative effects of food components on genetic expression profiles (nutritional transcriptomics and nutritional epigenomics effects), the occurrence and activity of proteins (proteomic effects) and/or the dose and temporal changes in cellular small molecular weight compounds (metabolomics effects) will [corrected] assist in identifying responders and non-responders. Expanding the information about similarities and differences in the "omic" responses across tissues will not only provide clues about specificity in response to bioactive food components but assist in the identification of surrogate tissues and biomarkers that can be used for predicting a response. Deciphering the importance of each of these potential sites of regulation will be particularly challenging but does hold promise in explaining many of the inconsistencies in the literature.

Anticarcinogenic Agents↗

PubMed2005

Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters.

cis-Acting regulatory elements are important molecular switches involved in the transcriptional regulation of a dynamic network of gene activities controlling various biological processes, including abiotic stress responses, hormone responses and developmental processes. In particular, understanding regulatory gene networks in stress response cascades depends on successful functional analyses of cis-acting elements. The ever-improving accuracy of transcriptome expression profiling has led to the identification of various combinations of cis-acting elements in the promoter regions of stress-inducible genes involved in stress and hormone responses. Here we discuss major cis-acting elements, such as the ABA-responsive element (ABRE) and the dehydration-responsive element/C-repeat (DRE/CRT), that are a vital part of ABA-dependent and ABA-independent gene expression in osmotic and cold stress responses.

Base Sequence↗

PubMed2026

ModiCal: A Targeted Calibration Workflow for Site-Specific m5C Validation by Nanopore Direct RNA Sequencing.

Accurate identification of RNA 5-methylcytidine (m5C) at the single-nucleotide resolution remains a central challenge in nanopore direct RNA sequencing (DRS). Current global scanning and modification-aware basecalling methods enable transcriptome-wide profiling but often yield high false-positive rates and lack site-specific accuracy. To address this, we repurposed ModiDeC, originally a de novo multimodification classifier, into a targeted, high-precision validation tool for RNA modification sites with prior biochemical knowledge. This was implemented through a three-step calibration workflow that alternates between biochemical and computational modules using the well-characterized m5C2278 site in 25S rRNA as a starting point. Baseline training uses short synthetic RNAs carrying either a methylated or unmodified C2278 as ground truth, followed by IVT-derived calibration and validation in methyltransferase knockout yeast. The baseline model accurately detected the bona fide m5C2278 site but initially produced off-target predictions. Iterative retraining with unmodified IVT signals progressively reduced and ultimately eliminated false positives while maintaining a strong signal at the bona fide site. The final model retained enzyme-dependent detection in wild-type versus knockout yeast and, when explicitly targeted, was also able to detect the second rRNA site, C2870, which remained invisible in the initial analysis. Application to native human prerRNA processing intermediates further resolved two distinct m5C deposition regimes on 28S rRNA, while generalization to dengue virus genomic RNA confirmed that the same calibration logic transfers across diverse RNA contexts. Together, this study establishes a reproducible and transferable framework that integrates biochemical validation with iterative neural network refinement, providing a route toward reliable site-specific m5C confirmation by nanopore direct RNA sequencing.

RNA Methylation↗

PubMed2026

Scalable, generalizable and uncertainty-aware integration of spatial multiomics across diverse modalities and platforms with SCIGMA.

Recent advances in spatial omics technologies have enabled simultaneous profiling of transcriptomic, proteomic, epigenomic, metabolomic and imaging data at high spatial resolution, offering unprecedented opportunities to dissect tissue complexity. However, integrating these diverse and large-scale spatial multimodal datasets remains a major computational challenge. We present SCIGMA, a scalable and generalizable deep learning framework for spatial multiomics integration. SCIGMA introduces an uncertainty-aware contrastive learning objective and multiview graph neural networks to preserve modality-specific signals while learning biologically meaningful joint representations. Unlike previous methods, SCIGMA provides spatially resolved uncertainty estimates, interpretably identifying regions of biological or technical heterogeneity. SCIGMA supports integration of up to five modalities, and its modular framework is extensible to future technologies with even more modalities. It also scales to more than 1 million spatial locations, enabling analysis of high-resolution datasets such as Visium HD and Xenium Prime. We evaluated SCIGMA across 19 datasets spanning 8 modalities, 10 tissues and 9 platforms. On benchmarkable datasets, SCIGMA outperformed other methods in spatial domain detection, modality preservation, feature reconstruction and reproducibility. SCIGMA identifies biologically meaningful structures, refined spatial domains and modality-specific regulatory programs, providing a robust, flexible and future-ready solution for scalable spatial multimodal integration.

Multiomics↗

PubMed2026

Dissecting spatial patterning and signaling with directional diffusion in spatial multi-omics.

Spatial multi-omics sequencing enables the simultaneous profiling of transcriptomics, proteomics, and epigenomics at a spatial resolution, offering insights into complex tissue organization and molecular regulation. However, the effective integration of multiple omics modalities in a spatial context remains a major challenge. Here, we present SpaDDM, a spatial multi-omics integration framework based on directional diffusion models (DDMs), which supports spatial pattern identification, cross-omics alignment, and inter-and intracellular signaling flow analysis. SpaDDM employs DDM-based graph networks to learn omics-specific representations by jointly incorporating spatial coordinates and molecular measurements within each modality, followed by an attention mechanism to align features across modalities. We benchmarked SpaDDM on diverse spatial multi-omics datasets, including transcriptomics-epigenomics and transcriptomics-proteomics combinations across multiple tissues and species. SpaDDM consistently outperformed existing methods by more accurately deciphering spatial tissue patterns and effectively reducing the boundary noise between spatial regions. Moreover, the learned low-dimensional coembedded representations of individual cells serve as integral mediators for inferring the signaling flows that underlie spatial patterning. Finally, we demonstrated that SpaDDM alignment of complementary information across multi-omics layers facilitates cross-omics translation and significantly improves the prediction of cell state alignments.

Multiomics↗

PubMed2026

Penicillium melinii promotes root growth through subtle host reprogramming across model and crop species.

Root development is highly responsive to microbial interactions, yet the mechanisms by which beneficial fungi promote root growth remain incompletely understood. Here, we identified Penicillium melinii 'isolate 2' through a screen of endophytic fungi isolated from Arabidopsis and characterized it as a promoter of root development in both Arabidopsis and crop species. We combined phenotyping in vitro, rhizotron, greenhouse and field assays with reporter and mutant analyses, transcriptomics, phytohormone profiling and sequencing and annotation of the fungal genome to investigate the basis of this interaction. P. melinii consistently stimulated root growth and modified root architecture across experimental systems and host species. These effects were associated with subtle but reproducible host transcriptional reprogramming, supporting a model in which the fungus fine-tunes endogenous developmental programmes rather than broadly perturbing stress or growth pathways. Genetic and reporter analyses further suggested that this interaction modulates root branching through localized developmental reprogramming. Genomic analysis provided a framework for understanding the fungal traits associated with this beneficial interaction. The conservation of the response across model and crop species supports the relevance of P. melinii as both a useful experimental system to study beneficial plant-fungus interactions and a promising candidate for improving root traits and crop performance.

Penicillium melinii↗

PubMed2004

Transcriptome analysis and gene expression profiles of early apoptosis-related genes in Streptococcus pyogenes-infected epithelial cells.

Epithelial cells are the initial sites of host invasion by group A Streptococcus pyogenes (GAS), and GAS infection of epithelial cells has been suggested to induce apoptosis. We previously reported that the induction of apoptosis is strongly associated with the protein F1-mediated invasion. We present here the gene expression profiles of the human epithelial HEp-2 cells during GAS-induced apoptosis, using serial gene analysis of expression (SAGE) analysis and macroarray analysis of apoptosis-related genes. Serial gene analysis of expression revealed the downregulation of voltage-dependent anion channels 1 and 2 genes and the upregulation of the cytochrome c oxidase and calcium binding protein genes (calpactin, calgizzarin and programmed cell death 6). Macroarray analysis and quantitative RT-PCR analysis also revealed that the genes for IL-1beta, IL-12 p35, IL12 p40, and GM-CSF are also markedly induced by GAS invasion. Furthermore, caspase-1, -9, and -14 genes are significantly upregulated during GAS invasion. These observations indicated that apoptosis associated with GAS invasion is mainly induced by mitochondrial dysfunction and calcium regulation as well as by stress, and that these transcriptional controls may regulate the cellular response to GAS invasion.

Apoptosis↗