SEARCH · PubMed Health

Results for “multi-omics data”

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 145 records · Page 8Linked to original sources

PubMed2026

A Multi-omics Regulated Cell Death Framework Defines Immune Phenotypes and Guides Precision Therapy in Colorectal Cancer.

Colorectal cancer (CRC) is molecularly and immunologically heterogeneous, contributing to variable treatment response. Because regulated cell death (RCD) intersects with tumor metabolism, immune regulation, and therapeutic susceptibility, we built an RCD-centered framework for CRC stratification. Multi-cohort transcriptomic data were used to infer RCD subtypes with non-negative matrix factorization (NMF) and non-negative least squares (NNLS). Genomic, bulk RNA-seq, single-cell RNA-seq, and spatial transcriptomic datasets were integrated to characterize subtype-associated biology. Machine-learning models were developed for immunotherapy response and survival-risk estimation. Candidate compounds were screened by GDSC2-based drug-sensitivity modeling and molecular docking, and FSTL3 was functionally assessed in vitro. The framework separated CRC samples into two RCD-related phenotypes resembling immune-hot and immune-cold states. RCD1 showed immune activation and higher mutational burden, whereas RCD2 showed immune-suppressed features, intratumoral heterogeneity, and aggressive biology. RCD-associated signatures showed potential for predicting immunotherapy response and survival risk. Dasatinib was prioritized for immune-cold, high-risk tumors, with preliminary evidence supporting its activity in CRC cells, while functional assays suggested a role for FSTL3 in growth, invasion, epithelial-mesenchymal transition, and apoptosis regulation. These findings suggest that RCD-based multi-omics analysis may refine CRC stratification and help generate therapeutic hypotheses.

Colorectal cancer↗

PubMed2026

Clinicopathological response and survival outcomes of HER2-low versus HER2-zero early breast Cancer: A systematic review and Meta-analysis.

BACKGROUND: Breast cancer is the most common malignant tumor in women. Human epidermal growth factor receptor 2 (HER2) is a key biomarker for classification and treatment. A subgroup with HER2-low expression has been identified, but existing evidence is heterogeneous. This systematic review and meta-analysis compared pathological response and survival outcomes between HER2-low and HER2-zero early-stage breast cancer to clarify prognostic features. METHODS: This study followed PRISMA guidelines and was registered in PROSPERO (CRD420251120506). PubMed, Embase, Web of Science, ClinicalTrials.gov, and major oncology conferences were searched through September 2025. Cohort studies of early-stage breast cancer comparing HER2-low (IHC 1+/2+ and ISH-negative) vs. HER2-zero with extractable pCR, DFS, or OS data were included. Studies involving HER2-positive patients or inconsistent definitions were excluded. Meta-analyses were performed using RevMan 5.3. RESULTS: Twenty-eight studies involving 115,182 patients were included. HER2-low patients showed significantly lower pCR rates (OR = 0.58, 95% CI: 0.52-0.65). DFS favored HER2-low (multivariate HR = 0.75, 95% CI: 0.69-0.83), especially in HR+ tumors, with a weaker effect in HR- cases. OS also favored HER2-low (HR = 0.80, 95% CI: 0.72-0.89), mainly driven by the HR- subgroup; no OS difference was seen in HR+ tumors. Sensitivity analyses and funnel plots indicated robust results with no apparent publication bias. Overall study quality was high (17 high-quality, 11 moderate-quality). CONCLUSION: HER2-low early breast cancer shows lower pCR after neoadjuvant therapy but better long-term survival. These findings support the clinical relevance of HER2-low as a biologically meaningful subgroup within HER2-negative disease, while its status as a stable and independent subtype still requires further validation through prospective studies, standardized testing, and multi-omics investigation.

Humans↗

PubMed2026

Reduced legacy precipitation decreases microbial community growth efficiency and alters soil organic carbon in a California grassland.

BACKGROUND: Changes in global patterns can leave a lasting legacy in semiarid grasslands by reshaping microbial growth dynamics and carbon cycling during the first wet-up in the autumn-a period known for intense microbial activity and significant carbon emissions. To study the lasting impacts of decreased winter rain, we implemented two precipitation regimes (100% vs. 50% mean annual precipitation) in California Mediterranean-climate grassland field plots. After the dry season, soils were rewetted in the laboratory with H218O and sampled at 0 h, 3 h, 24 h, 48 h, 72 h, and 168 h post rewet. We quantified CO2 efflux, measured microbial growth and mortality via quantitative 18O stable isotope probing and 16S rRNA gene amplicon sequencing, and characterized the soil organic carbon chemical composition, metagenomes, and metatranscriptomes. RESULTS: We found that reduced winter precipitation imposed a strong legacy effect on microbial turnover; despite maintaining similar respiration rates, microbial growth declined by ~1 order of magnitude, yielding decreased community growth efficiency (CGE = new biomass growth/respiration), and microbial mortality declined by ~2 orders of magnitude. Soil organic carbon also shifted from lipid-like, amino-sugar-like, and protein-like compounds (indicative of microbial necromass) to more oxidized lignin-like and tannin-like compounds (indicative of decomposing plant-derived compounds). Meta-omics revealed distinct metabolic strategies linked to CGE. At high-CGE, microbes appeared to consume more energetically favorable N-rich necromass (released via high microbial turnover); this allowed for increased amino acids and peptidoglycan biosynthesis and greater aromatic compound degradation, fueling further energy production and growth efficiency. At low CGE, communities had elevated carbohydrate metabolism and lipid turnover, consistent with increased investment in plant detritus degradation and membrane repair and maintenance rather than growth. CONCLUSIONS: Together, our findings demonstrate that reduced winter rainfall decreases microbial turnover following rewetting without a concurrent reduction in CO2 emissions. This shift results in persistently lower CGE, which has the potential to increase soil carbon loss as CO2. If such conditions are maintained over multiple years, these changes could reshape soil organic carbon stocks and alter the balance of grassland ecosystems under future climate scenarios. While our data suggest that sustained reductions in CGE may drive SOC decline, the magnitude and persistence of these effects depend on long-term environmental dynamics and warrant further investigation. Video Abstract.

Soil Microbiology↗

PubMed2026

Decoding age-stratified clinical and molecular heterogeneity in male breast cancer through multiomic profiling.

OBJECTIVE: Age-associated molecular heterogeneity is well described in female breast cancer but remains insufficiently characterized in male breast cancer (MBC). We profiled age-stratified clinical and molecular differences between younger (≤55 years) male breast cancer (YMBC) and older (>55 years) male breast cancer (OMBC). METHODS: We retrospectively analyzed 347 patients with MBC diagnosed at Fudan University Shanghai Cancer Center by integrating clinicopathological data, RNA sequencing, and whole-exome sequencing (WES). Survival, differential expression, and mutational signature analyses were performed. Tumor microenvironment features were inferred using xCell and ESTIMATE, and weighted gene co-expression network analysis (WGCNA) was conducted to identify age-associated co-expression modules. Candidate therapeutics were prioritized using the Genomics of Drug Sensitivity in Cancer (GDSC) resource and evaluated using patient-derived organoids (PDOs). RESULTS: Compared with OMBC, YMBC more frequently had human epidermal growth factor receptor 2 (HER2)-positive status (14.91% vs. 4.02%) and triple-negative tumors (4.92% vs. 1.78%), and had worse 5-year recurrence-free survival (hazard ratio=2.19, P=0.018). Transcriptomic analyses indicated enrichment of neural-related programs and reduced immune-related signaling in YMBC, and xCell/ESTIMATE supported lower immune infiltration. Consistently, WGCNA identified age-associated modules linking neural-related programs with reduced immune infiltration. Immunohistochemistry supported increased perineural invasion and lower CD8+ T cell infiltration in YMBC. GDSC-guided prioritization with PDO testing nominated sepantronium bromide (YM155) as a candidate vulnerability in YMBC. WES showed a higher NBPF10 mutation frequency in YMBC (54.5% vs. 14.3%, P<0.05). CONCLUSIONS: Integrated multi-omics profiling revealed age-stratified clinical and molecular heterogeneity in MBC. YMBC patients demonstrated inferior recurrence-free survival, neural signaling enrichment, an immune-cold microenvironment, and enriched NBPF10 mutations. These findings support age as a meaningful stratification variable in MBC risk assessment and treatment planning, and highlight the need for caution when considering treatment de-escalation in younger patients, while nominating YM155 as a candidate agent for prospective evaluation.

Male breast cancer↗

PubMed2026

MET Exon 14 Skipping Mutation in NSCLC: From Genomic Discovery to Biomarker-Guided Therapeutic Innovation.

INTRODUCTION: Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, and the MET exon 14 skipping mutation is a key oncogenic driver, which promotes tumor progression and provides a new direction for precision therapy. METHODS: A systematic search of English-language literature and clinical trial data related to the MET exon 14 skipping mutation from 2020-2025 was performed to summarize the role of the mutation and therapeutic advances. RESULTS: DNA-based next-generation sequencing (NGS), RNA-based NGS, and RT-qPCR were employed as the main detection methods. Preclinical models confirmed that mutations promote tumor progression by activating the RAS/MAPK pathway. Clinical trials have reported objective remission rates (ORR) of 46-68% for first-line treatment with MET inhibitors in NSCLC patients harboring MET exon 14 skipping mutations. DISCUSSION: MET exon 14 skipping mutation as a therapeutic target for NSCLC has made significant progress, and MET inhibitors are more advantageous than chemotherapy and immunotherapy, and have been recommended by national and international guidelines as a first-line treatment option. Additionally, NGS technology has the potential to dynamically monitor tumor evolution and drugresistant mutations, thereby helping to realize precision medicine. CONCLUSION: The MET exon 14 skipping mutation is an important target for the precision treatment of NSCLC, and MET-TKIs have remarkable efficacy but a prominent problem with drug resistance. The construction of a precision medicine system encompassing diagnosis, treatment, and drug resistance management through multi-omics research, technological innovation, and international collaboration is a key direction for improving prognosis.

Humans↗

PubMed2026

Multi-omics integrative analysis provides insight into potential molecular responses to sustained high water flow in common carp (Cyprinus carpio) cultured in recirculating aquaculture.

To investigate the potential molecular responses by which water flow intensity affects the growth of common carp (Cyprinus carpio) in a recirculating aquaculture system (RAS), a control group (CG, actual water velocity 0.3 cm/s) and three sustained flow treatment groups were established, including a low-flow group (LF, 1 body length per second, bl/s), a medium-flow group (MF, 2 bl/s), and a high-flow group (HF, 3 bl/s). After 12 weeks of culture in the RAS, growth performance was compared among groups under different flow intensities. The best-performing group and the control group were then selected for the determination of intestinal digestive enzyme activities, as well as transcriptomic and whole-genome bisulfite sequencing analyses of muscle tissue. The results showed that the specific growth rate and feed intake of the HF group were significantly higher than those of the other groups (P < 0.05), whereas no significant difference in feed conversion ratio was observed among groups. Compared with the CG group, lipase activity was significantly higher in the HF group (P < 0.05), while α-amylase and trypsin activities showed increasing trends without significant differences. RNA-seq identified a total of 273 differentially expressed genes, including 72 upregulated genes and 201 downregulated genes in the HF group relative to the CG group. These genes were mainly enriched in glycolysis, pyruvate metabolism, ATP metabolism, the pentose phosphate pathway, the insulin signaling pathway, the PPAR signaling pathway, and the adipocytokine signaling pathway, indicating that sustained high water flow induced a muscle transcriptional response characterized by remodeling of energy metabolism and substrate utilization. Whole-genome bisulfite sequencing analysis showed that DNA methylation in common carp muscle occurred predominantly in the CpG context. Differentially methylated regions between the HF and CG groups were mainly distributed in transcription-related regulatory regions, including promoters, CpG islands, and CpG island shores. In promoter regions, the number of hypermethylated regions in the HF group relative to the CG group was markedly higher than that of hypomethylated regions. Integrated analysis further identified two candidate genes showing both promoter differential methylation and differential expression, namely LOC109094644 and bcorl1, suggesting that adaptation to high water flow may involve IGF-related growth regulation and remodeling of upstream transcriptional programs. The qPCR results were consistent with the transcriptomic data. Taken together, within the tested range, a sustained water flow of 3 bl/s was more conducive to the growth of common carp in the RAS, which may be associated with enhanced lipid digestion and utilization, remodeling of the muscle energy metabolic network, changes in promoter methylation, and the coordinated regulation of key candidate genes. This study provides a theoretical basis for clarifying the exercise adaptation mechanism of common carp in recirculating aquaculture and for optimizing flow velocity parameters.

Animals↗

PubMed2025

Potential biomarkers for human Ascending aortic aneurysm identified through metagenomic and metabolomic analyses: A case-control study.

INTRODUCTION: Ascending aortic aneurysm (AsAA) is a high-risk cardiovascular condition; recent research indicates a possible association between gut microbiota, plasma metabolites, and the pathogenesis of AsAA. OBJECTIVE: This study aims to investigate the compositional and metabolic alterations in the gut microbiota of AsAA patients to identify potential biomarkers for AsAA. METHODS: This study enlisted 72 participants, comprising 44 individuals with AsAA and 28 healthy controls. All participants underwent examination for clinical features, and fecal and plasma samples were obtained for metagenomic and metabolomic studies. RESULTS: Metagenomic analysis revealed a significant reduction of 23 bacterial species in AsAA patients, including Bifidobacterium adolescentis, Bifidobacterium longum, Lactiplantibacillus plantarum, Enterococcus faecalis, and Streptococcus thermophilus, while 52 bacterial species, such as Prevotella copri, Phascolarctobacterium faecium, and Eubacterium ventriosum, were found to be enriched. Furthermore, we identified seven microbial co-abundance groups (CAGs), of which three (predominantly comprising Roseburia, Agathobacter, and Prevotella) were significantly elevated in AsAA patients, whereas one (predominantly comprising Escherichia) was substantially diminished. KEGG pathway enrichment analysis indicated that the biosynthesis of unsaturated fatty acids pathway displayed the most pronounced differences between groups. Metabolomics data revealed that 22 metabolites, including ceramides, were significantly elevated, while 8 metabolites, such as threonine, were notably downregulated. Moreover, clinical indicators like C-reactive protein (CRP) and complement components C3 and C4 have shown strong correlations with specific gut microbiota (Streptococcus, Prevotella) and plasma metabolites (threonine, ceramides). These findings indicate that inflammatory responses, metabolic dysregulation, and gut microbiota imbalance are pivotal in the etiology of AsAA. CONCLUSION: This study demonstrates substantial alterations in gut microbiota composition and plasma metabolites in patients with AsAA. Prevotella and ceramides exhibit potential as biomarkers for AsAA diagnosis. Furthermore, a synergy of Prevotella and ceramides may function as a potent disease prediction classifier, offering novel perspectives on the early diagnosis and targeted treatment of AsAA.

Humans↗

PubMed2026

Integrated single-cell and spatial transcriptomic analyses reveal malignant epithelial glycolytic heterogeneity and spatial niche remodeling during colorectal cancer progression.

Colorectal cancer (CRC) progression is shaped by metabolic reprogramming and complex interactions within the tumor microenvironment. However, the cellular heterogeneity, spatial organization, and clinical relevance of glycolytic activity in CRC remain incompletely understood. In this study, we integrated single-cell RNA sequencing, bulk transcriptomics, and spatial transcriptomics data to systematically characterize glycolytic heterogeneity in CRC. Glycolytic activity was quantified using five independent scoring methods, consistently showing that epithelial cells exhibited the highest glycolytic activity across the two single-cell cohorts. Stratification of CopyKAT-verified aneuploid malignant epithelial cells into high-glycolysis (HG) and low-glycolysis (LG) subgroups by glycolysis scores revealed that HG cells exhibited higher stemness scores and chromosomal copy number variations. Cell-cell communication analysis revealed that, compared with LG cells, HG cells exhibited increased interaction frequency and strength with immune and stromal populations, indicating enhanced malignant epithelial-microenvironment crosstalk. Spatial transcriptomics analyses further revealed that glycolytic activity varied across normal colorectal tissue, primary CRC, and colorectal liver metastases, accompanied by progressive remodeling of epithelial-associated spatial niches and MIF-mediated intercellular communication. Bulk transcriptomic analysis identified a glycolysis-related prognostic signature with robust predictive performance, which served as an independent prognostic factor for overall survival in CRC cohorts. Collectively, these findings indicate that glycolytic heterogeneity is a key feature of CRC malignant epithelial cells and is closely associated with tumor progression, microenvironmental remodeling, and clinical outcomes.

Humans↗

PubMed2026

A regulatory network underlying idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a progressive interstitial lung disease in which genetic susceptibility interacts with epithelial, immune, and mesenchymal remodeling. Although the chromosome 11p15.5 locus contains established IPF susceptibility signals near MUC5B and TOLLIP, the broader regulatory architecture of this region remains incompletely resolved. METHODS: We integrated IPF genome-wide association study summary statistics with methylation, expression, and protein quantitative trait loci using summary-data-based Mendelian randomization (SMR). SMR-prioritized candidates were evaluated in independent transcriptomic and methylation cohorts and further contextualized using microRNA, transcription-factor, protein-interaction, machine-learning, single-cell, and spatial transcriptomic analyses. Fibrosis-associated expression patterns were assessed in a bleomycin-induced pulmonary fibrosis rat model. RESULTS: The analyses recovered the established MUC5B and TOLLIP signals and prioritized BRSK2 as a comparatively underexplored candidate supported by eQTL-based SMR and independent molecular evidence. The BRSK2 pQTL association did not pass the HEIDI test and was therefore not interpreted as convergent protein-level genetic evidence. Network analyses linked BRSK2 to cell-cycle, metabolic-stress, and senescence-related programs, while cross-cohort machine learning prioritized FOXA2, CDC25B, and NFE2 as informative network features. Single-cell and spatial analyses localized BRSK2 preferentially to fibroblast and myofibroblast compartments and to regions with greater histological fibrosis severity. In fibrotic rat lungs, BRSK2 expression increased, whereas FOXA2 and CDC25B decreased at the transcript and protein levels. CONCLUSIONS: These findings refine the molecular landscape of the chromosome 11p15.5 IPF susceptibility locus and prioritize BRSK2 as a candidate component of an IPF-associated profibrotic fibroblast state. Its causal contribution, direct regulatory relationships, and therapeutic tractability require targeted mechanistic validation.

Idiopathic Pulmonary Fibrosis↗

PubMed2026

Multi-Omics Integration Identifies a Five-Gene Metabolic Signature With Experimental Validation in Clear Cell Renal Cell Carcinoma.

BACKGROUND: Clear cell renal cell carcinoma (ccRCC) is hallmarked by profound metabolic reprogramming; however, its intricate crosstalk with the tumor immune microenvironment (TIME) and its clinical ramifications remain inadequately elucidated. This study aims to systematically decipher the metabolic-immune interplay in ccRCC through multi-omics integration, with the goal of identifying robust prognostic biomarkers and actionable therapeutic vulnerabilities. AIMS: This study aims to systematically decipher the metabolic-immune interplay in clear cell renal cell carcinoma (ccRCC) through multi‑omics integration, and to identify robust prognostic biomarkers and actionable therapeutic vulnerabilities that can inform precision risk stratification and individualized treatment strategies. METHODS: We integrated bulk transcriptomic, genomic, and clinical data from multiple ccRCC cohorts. Differential expression and functional enrichment analyses were performed to characterize metabolic pathway alterations. Mendelian randomization (MR) was employed to infer causal relationships between metabolic disorders and ccRCC risk. A machine learning-based prognostic framework, incorporating SHAP (SHapley Additive exPlanations) for feature interpretability, was constructed and rigorously validated. TIME heterogeneity was dissected using deconvolution algorithms, while drug sensitivity, tumor mutation burden (TMB), and TIDE scores were utilized to assess therapeutic responses and immune evasion. Candidate gene function was evaluated through in vitro gain- and loss-of-function assays, with expression validated via TCGA, HPA, western blot, and qRT-PCR. RESULTS: Enrichment analysis identified coordinated dysregulation in lipid metabolism, energy homeostasis, and hypoxia response pathways. MR analysis confirmed lipid metabolism disorders as a causal risk factor for ccRCC. Our machine-learning model, centered on five core SHAP-identified features (SUCLA2, ACAT1, PC, SUCLG1, and HMGCS2), demonstrated superior predictive accuracy over conventional clinical staging. Immune profiling unveiled dichotomous TIME states: the low-risk group retained active immune surveillance, whereas the high-risk group was enriched with immunosuppressive subsets. Drug sensitivity screening pinpointed LY2109761 and carmustine as high-risk-specific candidate agents. Furthermore, TMB and TIDE analyses stratified high-risk patients displaying genomic instability and immune evasion phenotypes. Functionally, SUCLA2 knockdown significantly enhanced ccRCC cell proliferation and invasion, while its overexpression suppressed these malignant phenotypes, corroborating its tumor-suppressive role. Expression patterns of the hub genes were consistently validated across multi-level datasets and experimental assays. CONCLUSION: This study establishes a precision oncology framework for ccRCC by functionally linking metabolic biomarkers, immunophenotypes, and stratified therapeutic strategies. Importantly, we identify SUCLA2 as a potential functional tumor suppressor and a promising target for further mechanistic and translational investigation.

Humans↗

PubMed2026

Multi-omic biomarkers in cardiovascular disease: Discovery to clinical translation.

Cardiovascular disease (CVD) remains the leading cause of mortality worldwide, necessitating improved risk stratification and early detection strategies. Multiomics approaches that integrate genomics, transcriptomics, proteomics, metabolomics, and epigenomics offer unprecedented opportunities for biomarker discovery and precision medicine in cardiovascular care. This narrative review examines the current landscape of multiomics biomarkers for CVD, tracing their evolution from discovery to clinical translation. We synthesize evidence from recent studies evaluating the clinical utility of integrated omics approaches across diverse cardiovascular conditions, including atherosclerotic cardiovascular disease, heart failure, and atrial fibrillation. High-throughput proteomics has identified novel protein signatures that enhance cardiovascular risk prediction beyond traditional risk factors. Metabolomics has revealed pathway-specific biomarkers, including trimethylamine N-oxide and lipid species, associated with atherogenesis. Polygenic risk scores derived from genomic data demonstrate incremental value when combined with clinical risk scores. Multiomics biomarkers represent a transformative approach to cardiovascular risk assessment and disease management.

Humans↗

PubMed2026

Multi-omics panorama of glaucoma: Pathogenesis, biomarkers, and novel therapeutic strategies.

Glaucoma is a group of irreversible, blinding eye diseases characterized by progressive loss of retinal ganglion cells, leading to gradual visual field defects that severely impact patients' quality of life. Its complex pathophysiological mechanisms remain incompletely understood, limiting the development of early diagnostic and effective therapeutic strategies. Advances in omics technologies have provided new insights into elucidating the pathophysiology of glaucoma. We summarize specific alterations in genomics, transcriptomics, proteomics, metabolomics, epigenomics, and microbiomics associated with glaucoma. We emphasize the systematic analysis of disease mechanisms, identification of clinically applicable biomarkers, and discovery of novel therapeutic targets through the integration of these data. This approach paves new pathways for glaucoma subtype diagnosis and personalized treatment, while also outlining future research directions and challenges.

Humans↗

PubMed2026

Decoding the molecular basis of blue grain color codominance in Qingke: Integrative analysis of RNA-seq, DNA methylation, and miRNA-seq.

The grains on single spike of the F1 generation from the cross between blue- and white-grained Qingke (Hordeum vulgare L. var. nudum Hook. f.) are randomly distributed in blue and white colors. This study integrated data from RNA-seq, DNA methylation, and miRNA-seq to analyze this trait. The results showed that the HvF3'5'H gene is likely central to the development of this codominant phenotype. Through cross-validation of three omics approaches, it was found that the HvMYB gene targeted by miR858-z, as well as the WRKY24 and At3g44326 genes targeted by novel-m0152-5p, novel-m0153-5p, and novel-m0154-5p, are correlated with DNA methylation. qRT-PCR analysis confirmed that the four aforementioned genes exhibited variety-specific and developmental stage-specific expression patterns. This study dissects the regulatory network underlying the codominant blue and white grain color divergence on a single Qingke spike from a multi-omics perspective.

DNA Methylation↗

PubMed2025

PLSKO: a robust knockoff generator to control false discovery rate in omics variable selection.

MOTIVATION: Integrating the knockoff framework with any variable-selection method delivers stringent false discovery rate (FDR) control without recourse to p-values, offering a powerful alternative for differential expression analysis of high-throughput omics datasets. However, existing knockoff generators rely on restrictive modelling assumptions or coarse approximations that often inflate the FDR when applied to real-world data. RESULTS: We introduce Partial Least Squares Knockoff (PLSKO), an efficient, assumption-free generator that remains robust across diverse omics platforms. Our extensive simulations show that PLSKO is the only method to maintain FDR control with sufficient power in complex non-linear settings. Our semi-simulation studies drawn from RNA-seq, proteomics, metabolomics, and microbiome experiments confirm PLSKO generates valid knockoff variables. In pre-eclampsia multi-omics case studies, we combine PLSKO with Aggregation Knockoff to address the randomness of knockoffs and improve power, and demonstrate the method's ability to recover biologically meaningful features. AVAILABILITY AND IMPLEMENTATION: Our proposed algorithm is available on Github (https://github.com/guannan-yang/PLSKO) and Zenodo (https://doi.org/10.5281/zenodo.16879594).

Algorithms↗

PubMed2025

scBSP: a fast and accurate tool for identifying spatially variable features from high-resolution spatial omics data.

MOTIVATION: Emerging spatial omics technologies empower comprehensive exploration of biological systems from multi-omics perspectives in their native tissue location in 2D and 3D space. However, the limited sequencing depth, increasing spatial resolution, and growing spatial spots in spatial omics technologies present significant computational challenges in identifying biologically meaningful molecules with variable spatial distributions across various omics modalities. RESULTS: We introduce scBSP, an open-source, versatile, and user-friendly package for identifying spatially variable features in large-scale spatial omics data. scBSP demonstrates significantly enhanced computational efficiency, processing high-resolution spatial omics data within seconds, and exhibits robust cross-platform performance by consistently identifying spatially variable features with high reproducibility across various sequencing platforms. AVAILABILITY AND IMPLEMENTATION: scBSP is available for download from R CRAN at https://cran.r-project.org/web/packages/scBSP/index.html and PyPI at https://pypi.org/project/scbsp/.

Software↗

PubMed2025

HoloFoodR: a statistical programming framework for holo-omics data integration workflows.

SUMMARY: Holo-omics is an emerging research area that integrates multi-omic datasets from the host organism and its microbiome to study their interactions. Recently, curated and openly accessible holo-omic databases have been developed. The HoloFood database, for instance, provides nearly 10 000 holo-omic profiles for salmon and chicken under controlled treatments. However, bridging the gap between holo-omic data resources and algorithmic frameworks remains a challenge. Combining the latest advances in statistical programming with curated holo-omic data sets can facilitate the design of open and reproducible research workflows in the emerging field of holo-omics. AVAILABILITY AND IMPLEMENTATION: HoloFoodR R/Bioconductor package and the source code are available under the open-source Artistic License 2.0 at the package homepage https://doi.org/10.18129/B9.bioc.HoloFoodR.

Software↗

PubMed2025

Integrative metabolomic and proteomic analysis of diabetic kidney disease progression with younger-onset type 2 diabetes.

AIM: Younger-onset type 2 diabetes (YT2D) confers a disproportionately high risk of diabetic kidney disease (DKD), yet early biomarkers and underlying mechanisms remain poorly defined. We aimed to identify metabolites associated with DKD progression and integrate metabolomic and proteomic data to elucidate pathways involved in a multi-ethnic Asian cohort. MATERIALS AND METHODS: In this prospective study, 787 YT2D patients (diagnosed at ≤ age 40) were followed for a median of 5.7 years. DKD progression was defined as an annual decline in estimated glomerular filtration rate (eGFR) of ≥3 mL/min/1.73 m2 or ≥ 40% reduction in eGFR from baseline. Plasma metabolites were measured by nuclear magnetic resonance spectroscopy. Multivariable regression analysis was performed in a discovery (N = 550) and internal validation cohort (N = 237). Integrative metabolomic-proteomic analysis (N = 428) was performed using sparse partial least squares discriminant analysis (sPLS-DA). RESULTS: Ninety-eight metabolites were differentially expressed between DKD progressors and non-progressors, of which total branched-chain amino acids (BCAAs) (OR = 0.60, 95% CI 0.46-0.79), valine (OR = 0.62, 95% CI 0.48-0.81), and leucine (OR = 0.56, 95% CI 0.43-0.74) associated with DKD progression, independent of metabolic risk factors. Integrative analysis identified three components comprising 23 proteins and 30 metabolites, involved in the citrate cycle and apoptosis, which improved prediction of DKD progression beyond clinical risk factors (AUC 0.69-0.83). CONCLUSION: Lower plasma BCAA levels are independently associated with DKD progression in YT2D. Integrative multi-omics analysis highlights disruptions in metabolic and apoptotic pathways, providing insights into DKD pathophysiology and potential biomarkers for early risk stratification.

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↗