SEARCH · PubMed Health

Results for “Functional enrichment analysis”

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 55 records · Page 3Linked to original sources

PubMed2026

Integrative Network Analysis of Bioactive Compounds from Punica granatum L. Peel: Multi-Target Mechanisms in Wound Healing.

BACKGROUND: Wound-healing agents often have limited efficacy and require prolonged recovery times, prompting growing interest in developing herbal-based formulations. Among these, Punica granatum L. has attracted considerable attention because of its high polyphenolic content. Despite its widespread use, the precise pharmacological targets underlying its wound-healing effects remain poorly understood and require systematic investigation. OBJECTIVES: This study aimed to elucidate the underlying pharmacological mechanisms of the topical wound-healing properties of P. granatum L. using a network pharmacology approach. METHODS: Bioactive compounds of P. granatum L. and their potential target genes were identified using the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP), Similarity Ensemble Approach (SEA), and SwissTargetPrediction databases. Wound healing-related genes were retrieved from the GeneCards database. Genes intersecting P. granatum L. targets and wound healing-associated genes were subjected to functional enrichment analyses, including protein-protein interaction (PPI), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. The PPI network was further analyzed using Cytoscape, and the phytoconstituent-target interaction network was visualized using Gephi. These findings were validated using molecular docking. RESULTS: A total of 40 intersecting genes were identified as potential P. granatum L. targets involved in wound healing. Among these, EGFR, PTPN11, HRAS, IGF1R, and ESR1 were identified as key hub genes. Functional enrichment analysis indicated that the most significantly enriched signaling pathways included the MAPK, PI3K-Akt, EGFR tyrosine kinase inhibitor resistance, focal adhesion, and FoxO signaling pathways. Molecular docking analysis confirmed favorable binding of quercetin and ellagic acid to the hub targets EGFR, IGF1R, and ESR1. CONCLUSIONS: These findings elucidate the pharmacological pathways underlying P. granatum-mediated wound healing and suggest that P. granatum L. acts as a multi-target modulator in the wound-healing process.

Focal Adhesion↗

PubMed2025

Construction of a prognostic model for gastric cancer based on immune infiltration and microenvironment, and exploration of MEF2C gene function.

BACKGROUND: Advanced gastric cancer (GC) exhibits a high recurrence rate and a dismal prognosis. Myocyte enhancer factor 2c (MEF2C) was found to contribute to the development of various types of cancer. Therefore, our aim is to develop a prognostic model that predicts the prognosis of GC patients and initially explore the role of MEF2C in immunotherapy for GC. METHODS: Transcriptome sequence data of GC was obtained from The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO) and PRJEB25780 cohort for subsequent immune infiltration analysis, immune microenvironment analysis, consensus clustering analysis and feature selection for definition and classification of gene M and N. Principal component analysis (PCA) modeling was performed based on gene M and N for the calculation of immune checkpoint inhibitor (ICI) Score. Then, a Nomogram was constructed and evaluated for predicting the prognosis of GC patients, based on univariate and multivariate Cox regression. Functional enrichment analysis was performed to initially investigate the potential biological mechanisms. Through Genomics of Drug Sensitivity in Cancer (GDSC) dataset, the estimated IC50 values of several chemotherapeutic drugs were calculated. Tumor-related transcription factors (TFs) were retrieved from the Cistrome Cancer database and utilized our model to screen these TFs, and weighted correlation network analysis (WGCNA) was performed to identify transcription factors strongly associated with immunotherapy in GC. Finally, 10 patients with advanced GC were enrolled from Sun Yat-sen University Cancer Center, including paired tumor tissues, paracancerous tissues and peritoneal metastases, for preparing sequencing library, in order to perform external validation. RESULTS: Lower ICI Score was correlated with improved prognosis in both the training and validation cohorts. First, lower mutant-allele tumor heterogeneity (MATH) was associated with lower ICI Score, and those GC patients with lower MATH and lower ICI Score had the best prognosis. Second, regardless of the T or N staging, the low ICI Score group had significantly higher overall survival (OS) compared to the high ICI Score group. For its mechanisms, consistently, for Camptothecin, Doxorubicin, Mitomycin, Docetaxel, Cisplatin, Vinblastine, Sorafenib and Paclitaxel, all of the IC50 values were significantly lower in the low ICI Score group compared to the high ICI Score group. As a result, based on univariate and multivariate Cox regression, ICI Score was considered to be an independent prognostic factor for GC. And our Nomogram showed good agreement between predicted and actual probabilities. Based on CIBERSORT deconvolution analysis, there was difference of immune cell composition found between high and low ICI Score groups, probably affecting the efficacy of immunotherapy. Then, MEF2C, a tumor-related transcription factor, was screened out by WGCNA analysis. Higher MEF2C expression is significantly correlated with a worse OS. Moreover, its higher expression is also negatively correlated with tumor mutation burden (TMB) and microsatellite instability (MSI), but positively correlated with several immunosuppressive molecules, indicating MEF2C may exert its influence on tumor development by upregulating immunosuppressive molecules. Finally, based on transcriptome sequencing data on 10 paired tumor tissues from Sun Yat-sen University Cancer Center, MEF2C expression was significantly lower in paracancerous tissues compared to tumor tissues and peritoneal metastases, and it was also lower in tumor tissues compared to peritoneal metastases, indicating a potential positive association between MEF2C expression and tumor invasiveness. CONCLUSIONS: Our prognostic model can effectively predict outcomes and facilitate stratification GC patients, offering valuable insights for clinical decision-making. The identified transcription factor MEF2C can serve as a biomarker for assessing the efficacy of immunotherapy for GC.

Humans↗

PubMed2025

Identification of potential biomarkers and mechanisms for keloid disorder based on comprehensive bioinformatics analysis and machine learning algorithms.

BACKGROUND: Keloid disorder (KD) encompasses a spectrum of fibroproliferative dermal conditions, the pathogenesis remains complex and incompletely understood. This study sought to identify biomarkers and potential therapeutic targets for KD through an integrative bioinformatics approach and machine learning analysis of RNA sequencing data. METHODS: RNA sequencing was performed on skin tissue samples from 13 patients with KD and 14 healthy controls. Using weighted gene co-expression network analysis and differential expression analysis revealed differentially expressed key module genes, and the CytoHubba plugin identified candidate genes. Subsequently analyzed using least absolute shrinkage and selection operator (LASSO) and support vector machine recursive feature elimination (SVM-RFE) methods to pinpoint feature genes associated with KD. Following this, biomarkers were determined through expression level validation, enrichment analysis, and immune infiltration analysis. RESULTS: A total of 420 differentially expressed key module genes were identified, and the top 10 genes with DMNC values were selected as candidate genes. Five feature genes were selected through LASSO and SVM-RFE, with NID2, MFAP2, COL8A1, and P4HA3 showing significant expression differences between KD and control samples, along with consistent expression patterns across datasets, identified as potential biomarkers. These four biomarkers were proved to possess high diagnostic potential, and they were found to exhibit significant positive correlations with one another. Functional enrichment analysis indicated that the primary KEGG pathways associated with these biomarkers included "steroid hormone biosynthesis" and "cytokine-cytokine receptor interaction." Moreover, immune infiltration analysis revealed that the four biomarkers were negatively correlated with type 17 T helper cells and positively correlated with 15 immune cell types, including activated B cells and central memory CD4 T cells. CONCLUSION: In conclusion, NID2, MFAP2, COL8A1, and P4HA3 were identified as key biomarkers for KD, offering new avenues for more targeted and effective diagnostic and therapeutic strategies for managing this condition.

Humans↗

PubMed2026

Semaglutide treatment in MOSH is associated with altered DNA methylation patterns of genes related to glycolipid metabolism.

Male obesity-associated secondary hypogonadism(MOSH) is a common disease among severely obese male patients. Although surgical interventions have demonstrated clinical benefits, a subset of patients continue to experience MOSH following surgery. Therefore, this study aims to investigate epigenetic changes associated with the use of the weight-loss drug Semaglutide in MOSH, focusing on DNA methylation and miRNA expression. In this exploratory study, samples were classified into three groups: a control group (n = 2), a MOSH group (n = 7), and a follow-up group (n = 4). DNA methylation analysis was performed on all samples, while miRNA sequencing was conducted on a subset of the samples: 2 from the control group, 7 from the MOSH group, and 2 from the follow-up group. Differentially expressed miRNAs (DEMs) were analyzed through the R package "limma", and the methylation level of CpG sites was analyzed based on the methylation β value, obtaining differentially methylated genes (DMGs). The functional enrichment analysis of miRNA target genes and methylation change genes was conducted using the R package "clusterProfiler". Finally, the regulatory networks of miRNA and methylation genes as well as the protein-protein interaction (PPI) network were analyzed. A total of 6 DEMs were screened out. The target genes of these DEMs were mainly enriched in pathways such as ATP binding, phosphorylation, cell adhesion, and Glycosphingolipid biosynthesis. Eighty DMGs were identified, and the largest number of DMGs were found in the X chromosome. In the regulatory network of DMGs and DEMs, hsa-miR-423-5p regulates most of these DMGs. Moreover, the PPI network shows that DPP6, DPP10, CACNA1C, and CNTNAP2 are the proteins with the strongest connectivity. Notably, differential CpG methylation changes were observed on chromosome 7, indicating a potential region of epigenetic alteration in MOSH; however, the biological and functional relevance of these changes remains unclear. Collectively, these findings suggest that Semaglutide treatment in MOSH may be associated with concurrent alterations in DNA methylation and miRNA expression, implicating genes related to energy and glycolipid metabolism, including DPP6, DPP10, CACNA1C, and CNTNAP2. These results are exploratory and hypothesis-generating, providing preliminary observations to inform future validation studies.

Semaglutide↗

PubMed2026

Mapping the Immune cell-specific gene regulatory network in bipolar disorder: A framework from scTWMR to exploratory drug-target annotation.

BACKGROUND: Although the involvement of the immune system in the genetic susceptibility of bipolar disorder (BD) is widely acknowledged, the causal relationship between gene expression in specific immune cell subtypes and BD requires systematic elucidation. METHODS: We implemented an analytical framework integrating single-cell transcriptome-wide Mendelian randomization (scTWMR) with colocalization analysis. This approach utilized cis-expression quantitative trait loci (cis-eQTLs) derived from 14 distinct immune cell types as instrumental variables to interrogate BD genome-wide association study (GWAS) summary statistics (comprising 41,917 cases and 371,549 controls). Subsequent investigations encompassed functional enrichment analysis, protein-protein interaction (PPI) network construction, phenome-wide association study (PheWAS), and performed an exploratory drug-target annotation. RESULTS: Our analysis identified 33 gene-immune cell associations. Colocalization analysis provided robust evidence (PPH4 > 90%) for shared causal variants implicating the MAD1L1, APOM, and NFKBIL1 loci. Significantly enriched biological pathways included cell cycle regulation, circadian rhythm entrainment, and neuroinflammation. The PPI network revealed a core regulatory module centered on histone-encoding and immune-related genes. Exploratory drug-target annotation nominated compounds for further investigation for compounds targeting APOM, TMEM258, and NFKBIL1. CONCLUSION: This study systematically delineates a genetically supported regulatory network of immune cell-specific gene expression in BD, predominantly implicating CD8⁺ effector T cells, plasma cells, and B cells. The findings corroborate established pathological pathways while uncovering novel cell type-specific therapeutic targets, thereby providing a genetic framework for prioritizing candidate targets for future investigation.

Bipolar disorder↗

PubMed2026

Integrative analysis of the roles and prognostic value of RNA-binding proteins in papillary renal cell carcinoma.

RNA-binding proteins (RBPs) serve essential roles in various cancer types, but their functions in papillary renal cell carcinoma (pRCC) have not been elucidated to date. In our work, differentially expressed RBPs in pRCC were identified after acquisition of RNA-sequencing and clinical data related to pRCC from The Cancer Genome Atlas database(TCGA). Functional enrichment analysis and protein interaction network analysis, along with univariate and multivariate Cox regression analyses, were performed to uncover potential biological effects of the identified RBPs and screen the hub RBPs for pRCC prognosis. We identified 251 up-regulated and 129 down-regulated RBPs, and filtered out seven hub RBPs, namely, SRSF8, CD3EAP, HBS1L, ELAC2, MRPL34, NOP2 and IGF2BP2, for their prognostic relevance. A prognostic risk score model for overall survival of pRCC patients was constructed based on the seven hub RBPs. Further analysis showed that the low-risk group had higher survival rate than the high-risk group in both training and validation cohorts. The predictive accuracy was verified in the Human Protein Atlas database.In addition, we introduced the GSE15641 dataset from the Gene Expression Omnibus (GEO) database for independent external validation, and confirmed the expression levels of HBS1L, MRPL34 and IGF2BP2 through real-time quantitative PCR (RT-qPCR) and Western blotting (WB) using human renal tubular epithelial cell line HK-2 and human papillary renal cell carcinoma cell line Caki-2. In pRCC, CD3EAP was significantly elevated, while ELAC2, IGF2BP2, MRPL34, SRSF8 and HBS1L were significantly reduced. There was no significant difference between tumor and normal tissues in NOP2 expression. Risk score and tumor grade were independent prognostic factors associated with overall survival. In addition, we established a nomogram based on the seven prognostic RBPs to help predict overall survival at 1-3 years. In conclusion, seven differentially expressed hub RBPs were identified as potential prognostic biomarkers for pRCC. Our prognostic model might serve as a support for better treatment decision-making. Our work could provide potential new ideas for diagnosis and research on targeted drugs for pRCC.

Bioinformatics↗

PubMed2026

Integrative analysis of transcriptome and chromatin accessibility reveals promoter-proximal regulation and identifies candidate ABC transporters associated with cold stress responses in maize.

BACKGROUND: Low-temperature stress is a formidable environmental constraint that severely limits the growth and productivity of maize (Zea mays L.), particularly during the highly vulnerable early seedling stage. While cold tolerance is a critical agronomic objective, the integrated transcriptional and epigenetic regulatory mechanisms that govern this trait remain largely elusive. Characterizing these coordinated molecular networks is fundamental to the genetic enhancement of cold resilience in maize. METHODS: Using two maize inbred lines contrasting in chilling response (ZHB12 tolerant, B73 sensitive), we performed integrative time‑course RNA‑seq and ATAC‑seq to thoroughly and systematically characterize the precise dynamic interplay between gene expression and chromatin accessibility under cold stress conditions at the seedling stage. RESULTS: Physiological assessments confirmed that ZHB12 possesses superior cold tolerance, manifested by significantly attenuated electrolyte leakage and reduced foliar damage compared to B73. Transcriptomic profiling revealed a massive, time-dependent divergence in gene expression between the two genotypes, with a major regulatory transition identified at 24 h of cold exposure. Functional enrichment analysis demonstrated that ZHB12 preferentially activates a robust defense repertoire, including Photosystem II electron transport, diterpenoid biosynthesis, and ATP biosynthetic pathways. Notably, multiple ATP-binding cassette (ABC) transporter genes were coordinately upregulated under chilling, suggesting their potential involvement in cellular homeostasis. ATAC-seq analysis indicated that cold stress is associated with chromatin remodeling in ZHB12, with increased accessibility observed in proximal promoter regions. Integrative analysis identified a core set of dual-responsive genes, in which increased promoter accessibility coincided with transcriptional upregulation. These genes were predominantly enriched in transporter activity and transcriptional regulation, suggesting potential epigenetic link to the superior stress response of ZHB12. CONCLUSION: Our findings reveal extensive transcriptional and chromatin accessibility changes in ZHB12 under cold stress. The observed associations between promoter accessibility and gene activation, particularly in genes involved in transport processes, highlight candidate regulators potentially contributing to cold tolerance. This study provides a molecular framework and identifies high-value candidate genes that may inform future efforts in breeding cold-tolerant maize, pending functional validation.

Zea mays↗

PubMed2025

Identification and validation of the important role of KIF11 in the development and progression of endometrial cancer.

BACKGROUND: Human kinesin family member 11 (KIF11) plays a vital role in regulating the cell cycle and is implicated in the tumorigenesis and progression of various cancers, but its role in endometrial cancer (EC) is still unclear. Our current research explored the prognostic value, biological function and targeting strategy of KIF11 in EC through approaches including bioinformatics, machine learning and experimental studies. METHODS: The GSE17025 dataset from the GEO database was analyzed via the limma package to identify differentially expressed genes (DEGs) in EC. Functional enrichment analysis of the DEGs was conducted using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. DEGs were further screened for hub genes through protein-protein interaction (PPI) network analysis and machine learning. The role of the hub gene KIF11 in EC was analyzed using clinical data from the TCGA database. The expression of KIF11 in EC was subsequently validated in clinical samples. In vitro experiments were utilized to evaluate the effects of KIF11 on biological functions such as proliferation, migration, apoptosis, and the cell cycle in endometrial cancer cells. RESULTS: A total of 877 DEGs, which are widely involved in important biological processes such as cell division, tubulin binding, and the cell cycle, were identified. Through PPI network analysis and machine learning, KIF11 was selected as the hub gene for subsequent analysis and experimental validation. An analysis of TCGA data revealed that KIF11 is highly expressed in EC and is associated with tumor grade, stage, and a low survival rate. The overexpression of KIF11 in tumor tissues was further confirmed in EC patient samples. KIF11 knockdown had inhibitory effects on cell proliferation, migration and invasion. Flow cytometry analysis revealed that KIF11 knockdown induced G2/M phase arrest and promoted apoptosis in EC cells. CONCLUSION: Our study demonstrated that KIF11 was upregulated in EC and was strongly associated with a poor prognosis. Notably, we found that reduced KIF11 expression inhibited EC cell proliferation, migration and invasion. KIF11 knockdown caused more EC cells to arrest in the G2/M phase and undergo apoptosis. The findings of our study emphasized that KIF11 may be a promising prognostic biomarker and therapeutic target for EC patients.

Humans↗

PubMed2026

Phosphoproteomic Profiling of Early-Stage Non-Small Cell Lung Cancer Provides Preliminary Evidence of Phosphorylation-Regulated Rho GTPase Signaling Driving Cytoskeletal Remodeling, Angiogenesis, and Cell Cycle Progression.

Non-small cell lung cancer (NSCLC) is the primary cause of cancer-related deaths worldwide. This can be attributed to the difficulty in early detection and the limited efficacy of available treatments, partly due to an incomplete understanding of the disease biology. Identification of key proteins involved in early-stage progression and understanding the underlying mechanisms can greatly contribute to the development of diagnostic and treatment strategies for NSCLC. Quantitative phosphoproteomic analysis was done on paired tumor tissues and adjacent normal lung tissues from early-stage NSCLC adenocarcinoma (LUAD) patients to allow for the identification of proteins with differential phosphorylation and their associated pathways. A total of 6483 phosphoproteins were identified, with 1229 proteins having significantly higher phosphorylation and 701 proteins having significantly lower phosphorylation in the tumor tissues. All MS data were deposited in ProteomeXchange with the identifier PXD071583. Function enrichment analysis showed that the differentially phosphorylated proteins and phosphosites were primarily involved in Rho GTPase signaling and cytoskeleton remodeling. Analysis of protein interaction networks suggests that the predicted kinase activity likely drives malignant transformation in NSCLC LUAD, presumably through Rho GTPase-mediated angiogenesis and cell cycle progression. More importantly, this study identified several protein phosphosites with differential phosphorylation and inferred kinase-phosphosite activities that have not previously been reported in NSCLC LUAD.

Humans↗

PubMed2026

Genomic and proteogenomic insights into Spontaneous Coronary Artery Dissection (SCAD): A systematic review of emerging multi-omic evidence.

BACKGROUND: Spontaneous coronary artery dissection (SCAD) is a major cause of myocardial infarction in young women without traditional cardiovascular risk factors (Hayes et al., 2018; Adlam et al., 2018 [1, 2]). Despite growing awareness, its biological underpinnings remain incompletely understood, and clinical management is largely based on observational evidence rather than mechanistic insight (Saw et al., 2014; Lettieri et al., 2015; Steg et al., 2024 [3-5]). OBJECTIVES: To systematically integrate genomic, epitranscriptomic, proteomic, and metabolomic data in order to characterize the multi-omic architecture of SCAD and identify potential biomarkers and therapeutic targets. METHODS: A systematic review was conducted in accordance with the PRISMA 2020 statement (Arbelo et al., 2023 [6]). PubMed/MEDLINE was searched for original studies investigating genomic and multi-omic features of SCAD. Data were extracted on study design, patient characteristics, identified variants, circulating biomarkers, and implicated biological pathways. Functional enrichment analysis was performed using the DAVID bioinformatics resource (Page et al., 2021 [7]). RESULTS: A total of 16 studies were included. Genome-wide association studies consistently identified susceptibility loci related to arterial structure and extracellular matrix integrity, including ADAMTSL4, PHACTR1/EDN1, LRP1, and FBN1 (Huang et al., 2009; Saw et al., 2020; Turley et al., 2020 [8-10]). Rare variant analyses further supported the role of genes involved in extracellular matrix remodeling and vascular smooth muscle cell function, including COL3A1, COL4A1/2, SMAD3, and TLN1 (Adlam et al., 2023; Turley et al., 2021, 2019; Carss et al., 2020; Zekavat et al., 2022; Wang et al., 2022 [11-16]), while ancestry-specific signals such as TSR1 variants were observed in distinct populations (Turley et al., 2023 [17]). Proteogenomic approaches linked genetic susceptibility loci to circulating proteins involved in matrix remodeling and inflammation, including cathepsin B and ECM1 (Maioli et al., 2010 [18]). Epitranscriptomic analyses identified differential microRNA expression profiles associated with vascular injury and repair pathways (Sun et al., 2019 [19]). CONCLUSIONS: SCAD is characterized by a complex, multi-layered biological architecture involving genetic susceptibility, extracellular matrix dysregulation, and vascular signaling pathways. Integration of multi-omic data provides novel insights into disease mechanisms and highlights potential biomarkers and targets for precision medicine approaches in SCAD.

Animals↗

PubMed2026

Comprehensive Analysis of miRNAs and Predicted Protein Interaction Networks in Skeletal Muscle Development of Myostatin-Deficient Rabbits.

Myostatin (MSTN), encoded by the MSTN gene, is a critical negative regulator of skeletal muscle mass. This study aims to identify and characterize the miRNAs involved in the development of the double-muscling phenotype in MSTN-deficient rabbits. We performed high-throughput sequencing to analyze the miRNA expression profiles in gluteus maximus tissue from wild type (MSTN+/+) and MSTN-KO (MSTN+/- and MSTN-/- inclusive) rabbits. Differentially expressed miRNAs (DEmiRNAs) were identified, and their potential target genes were predicted. Functional enrichment analysis of these target mRNAs was conducted using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database to elucidate the involved biological pathways and regulatory networks. A total of 25 DEmiRNAs (13 downregulated and 12 upregulated, |log2FC| ≥ 1.0, adjusted p < 0.05) and 1178 differentially expressed mRNAs (408 upregulated and 770 downregulated, |log2FC| ≥ 2.0, adjusted p < 0.05) were identified in MSTN-KO compared to MSTN+/+ rabbits. Bioinformatics analysis revealed that the target genes of these DEmiRNAs were significantly enriched in key pathways governing muscle growth and metabolism, including the PI3K-Akt signaling pathway, MAPK signaling pathway, and pathways related to ECM-receptor interaction and insulin signaling. Notably, many predicted target mRNAs are expressed by genes that encode key inhibitors of myogenesis (e.g., HDAC4) and major extracellular matrix components (e.g., COL4A3, POSTN). Our results demonstrate that MSTN deficiency induces a distinct and widespread change in the miRNA expression landscape of skeletal muscle.

Animals↗

PubMed2026

Unveiling novel transcriptomic prognostic biomarkers for specific breast cancer subtypes and treatment regimens.

BACKGROUND: Breast cancer (BRCA) is the most common cancer in women worldwide, yet current gene expression panels offer limited insight into treatment responses across different subtypes and therapies. This study aimed to identify reliable biomarkers for predicting treatment outcomes in specific BRCA subtypes and treatment regimens. METHODS: This study analyzed transcriptomic data from The Cancer Genome Atlas to identify differentially expressed genes (DEGs) in patient groups treated with different combinations of hormone therapy (H), chemotherapy (C), radiotherapy (R), and targeted therapy (T). Non-negative matrix factorization clustering was performed to stratify patients into clusters representing different BRCA subtypes. Functional enrichment analysis was performed, and survival assessments were conducted using the METABRIC dataset. RESULTS: A total of 1,148 DEGs were identified across treatment regimens, with 75 common DEGs shared across multiple regimens. Among these, 12 candidate biomarkers were associated with luminal subtypes treated with H, including LRP1B, of which high expression predicted cancer recurrence. In triple-negative breast cancer (TNBC) treated with C, 76 candidate biomarkers were identified, including TTYH1 for recurrence and ANXA8L1 and MPZ for non-recurrence. Functional analyses identified intermediate filament organization and keratinization as pathways associated with specific candidate biomarkers of TNBC following C. Survival analysis using METABRIC strengthened the prognostic ability of LRP1B and TTYH1 to predict worse survival and ANXA8L1 and MPZ to predict prolonged survival, with four additional prognostic biomarkers. CONCLUSION: This study identified gene expression prognostic biomarkers for luminal and TNBC subtypes, thereby supporting personalized therapies. Further experimental validation is required to confirm these findings for clinical application. CLINICAL TRIAL REGISTRY: No.

Breast cancer↗

PubMed2025

Chromosome-Level Genome Assembly of Eden's Whale Clarifies the Taxonomy and Speciation of Bryde's Whale Complex.

Eden's whale (Balaenoptera edeni), a poorly understood baleen cetacean, has long been shrouded in taxonomic ambiguity due to limited genomic resources, obscuring its distinction from closely related species and its position within the cetacean Tree of Life. In this paper, we present a high-quality chromosomal-level genome of B. edeni and conduct comparative genomic analyses to address long-standing taxonomic confusion and elucidate speciation of balaenopterids. Our phylogenomic analysis and demographic reconstruction reveal that B. edeni is a distinct sister to Bryde's whale (Balaenoptera brydei), sharing a common ancestor that diverged approximately 7.84 million years ago during the late Miocene. Their genetic divergence exceeds typical intraspecific variation in whales, supporting the reinstatement of B. brydei as a valid species. Chromosomal syntenic analyses suggest that macro-fragment inversions contributed to speciation in balaenopterid whales and uncover unexpected large-scale complex genome rearrangements in Bryde's whale, offering novel insights into cetacean genome evolution. Functional enrichment analysis of inverted regions between B. edeni and Balaenoptera musculus indicates their predominant association with metabolism and biosynthesis, as well as responses to various substances, stress, and stimuli. These genomic resources for B. edeni not only lay a critical foundation for comparative genetic and evolutionary research of cetaceans but also advance our understanding of the taxonomy and evolutionary dynamics of the Bryde's whale complex, with broader implications for baleen whale conservation and biodiversity.

Animals↗

PubMed2026

Zhiling Jiangya decoction treats hypertension in rats: An integrative study of network pharmacology, immune infiltration, molecular simulation, and 16S rDNA sequencing.

OBJECTIVE: This study integrated network pharmacology, immune infiltration analysis, molecular docking, molecular dynamics simulation, ADMET prediction, 16S rDNA sequencing, and rat experiments to elucidate the potential mechanisms underlying the antihypertensive effects of Zhiling Jiangya Decoction (ZLJYD). METHODS: Active compounds and their potential targets were screened from the PubChem, TCMSP, NovoPro, and SwissTargetPrediction databases. Hypertension-related targets were retrieved from the OMIM and GeneCards databases, and overlapping targets were identified. The STRING database and Cytoscape 3.10.1 software were used to construct a protein-protein interaction network and a herb-component-target-disease network. Gene Ontology functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were performed to identify the key biological processes and signaling pathways involved. Using the CIBERSORT algorithm combined with correlation analysis, we investigated the association between key targets and immune cell infiltration. Molecular docking, molecular dynamics simulations, and ADMET predictions were performed to assess the binding stability and pharmacokinetic properties of the main compounds with their corresponding targets. Finally, the antihypertensive efficacy of ZLJYD was validated using a spontaneously hypertensive rat model, and alterations in gut microbiota were analyzed using 16S rDNA sequencing. RESULTS: A total of 123 active compounds and 267 hypertension-related targets of ZLJYD were identified. Enrichment analysis revealed that these targets were primarily associated with the PI3K-Akt signaling pathway and lipid and atherosclerosis pathways. Immune infiltration analysis suggested that the therapeutic effects of ZLJYD may involve the regulation of follicular helper T cells, naïve B cells, and naïve CD4⁺ T cells. Molecular docking and dynamics simulations supported the stable binding of key compounds to their target proteins, while ADMET predictions indicated favorable pharmacokinetic properties and safety profiles. Rat experiments demonstrated that ZLJYD significantly reduced blood pressure in spontaneously hypertensive rats, partially alleviated gut microbiota dysbiosis, and altered microbial community structure and phylogenetic diversity. CONCLUSION: This study systematically elucidates the potential mechanisms underlying the antihypertensive effects of ZLJYD through multiple components, targets, and pathways, particularly immune regulation and gut microbiota remodeling. These findings provide mechanistic insights into its potential therapeutic application.

16S rDNA sequencing↗

PubMed2026

Identification of candidate genes for reproductive traits in Chinese Holstein cattle using single-step genome-wide association study.

In dairy farming, reproductive efficiency is vital to both profitability and sustainability. However, years of selective breeding for increased milk yield have adversely affected reproductive potential. This study aimed to pinpoint genomic regions and identify potential candidate genes associated with reproductive traits in Chinese Holstein cattle. In this study, a single-step genome-wide association study (ssGWAS) was conducted using 33,202 phenotypic records from 16,379 animals, 55,244 pedigree records, and genomic data from 1,698 cows. These data were integrated into the ssGWAS analysis, resulting in a total pedigree structure of 21,635 animals. A total of 12 significant markers were identified for calving interval (IC), days open (DO), number of services per conception (NS), and conception rate (CR). Among these significant SNPs, three SNPs were for IC, two SNPs were for DO, three SNPs were for NS, and four SNPs were for CR. Several promising candidate genes located near these SNPs have been identified, including SFXN4, B3GAT2, GRK5, PRDX3, and MTHFD1L, highlighting their potential involvement in fertility-related biological processes. Furthermore, functional enrichment analysis identified significant enrichment of pathways associated with cell adhesion and embryonic development, suggesting a potential mechanistic role for DSG family members (DSG1, DSG2, DSG3, and DSG4) in fertility regulation. Collectively, our findings enhance understanding of the complex genetic basis of reproductive traits in dairy cattle and may offer a valuable set of genomic targets for precision breeding of Chinese Holsteins. Integrating these markers into genomic selection programs may contribute to genetic improvements in reproductive efficiency and support the long-term sustainability of dairy production.

Animals↗

PubMed2024

Comprehensive analysis of circRNA-miRNA-mRNA network related to angiogenesis in recurrent implantation failure.

BACKGROUND: Abnormal endometrial blood flow causes a decrease in endometrial receptivity and is considered a relatively independent risk factor for recurrent implantation failure (RIF). This study aimed to explore the potentially functional circRNA-miRNA-mRNA network in RIF, and further explore its mechanism. METHODS: Datasets were downloaded from the GEO database to identify differentially expressed circRNAs, miRNAs and mRNAs. The circRNA-miRNA-mRNA and PPI networks were constructed using Cytoscape 3.6.0 and the STRING database, the hub genes were identified with the cytoHubba plug-in, and a circRNA-miRNA-hub mRNA regulatory sub-network was constructed. Then, GO and KEGG pathway enrichment analyses of the hub genes were performed to comprehensively analyze the mechanism of hub mRNAs in RIF. Due to the results of circRNAs-miRNAs-hub mRNAs regulatory network, we verified the expression of circRNA_0001721, circRNA_0000714, miR-17-5p, miR-29b-3p, HIF1A and VEGFA in the RIF mouse model by qRT‒PCR and western blotting. RESULTS: We initially identified 175 DEmRNAs, 48 DEmiRNAs and 56 DEcircRNAs in RIF associated with angiogenesis and constructed a circRNA-miRNA‒mRNA network and PPI network. We further identified six hub genes in the acquired network. Based on these genes, functional enrichment analysis revealed that the HIF-1 signaling pathway plays a vital role in endometrial angiogenesis in RIF. In addition, the interaction networks of circRNA_0001721/miR-17-5p/HIF1A and the circRNA_0000714/miR-29b-3p/VEGFA axis were predicted. In the RIF mouse model, circRNA_0001721, circRNA_0000714, HIF1A and VEGFA were down-regulated, whereas miR-17-5p and miR-29b-3p were up-regulated according to qRT‒PCR and western blotting. CONCLUSION: This study revealed that the HIF-1 signaling pathway plays a vital role in endometrial angiogenesis in RIF. The circRNA_0001721/miR-17-5p/HIF1A and circRNA_0000714/miR-29b-3p/VEGFA axes might play a role in the pathogenesis of endometrial angiogenesis in RIF.

MicroRNAs↗

PubMed2026

Bidirectional causal relationships between plasma proteins, neuroimaging metrics and risk of Alzheimer's disease.

BACKGROUND: Changes in neuroimaging metrics are among the first detectable pathophysiological alterations in Alzheimer's disease (AD). Proteins are closely linked to fluctuations in neuroimaging metrics. Therefore, the analysis of the proteomic signature associated with neuroimaging metrics holds significant promise for uncovering therapeutic targets that contribute to AD. METHODS: GWAS data concerning the Brain Imaging Data Structure (BIDs). The AD cohort comprised a total of 401,661 individuals diagnosed with AD, alongside 10,520 control participants. For a bidirectional MR analysis involving neuroimaging metrics, proteomics, and AD, the methods utilized included inverse variance weighted (IVW), MR Egger, weighted median, weighted mode, and the Wald ratio approaches. RESULTS: We identified 12 neuroimaging metrics that demonstrate significant relevance to AD (thickness of the left total hemisphere, volume of the right thalamus, and et al.). These metrics are structural magnetic resonance imaging (MRI) biomarkers that remain stable throughout the entire course of AD, from the preclinical stage through mild cognitive impairment (MCI) to dementia. Additionally, we found a substantial number of 1633 proteins that also show a noteworthy causal relationship with AD. Functional enrichment analysis indicated that these proteins were predominantly focused within various pathways linked to AD, encompassing those involved in the synaptic vesicle cycle, synaptic membranes, neurotransmitter release, and the activity of GABA receptors. In addition, our research indicates that the significant relationships observed between the identified proteins and AD are influenced by neuroimaging metrics. Notably, we found that these neuroimaging metrics play a crucial role in mediating a substantial 67% of the inverse relationship that exists between PTPRC and the phenotypic characteristics associated with AD. CONCLUSIONS: This study successfully establishes a connection between proteomic and neuroimaging metrics, as well as the AD that influence them. By creating this relationship, the research offers important information that aids in comprehending the intricate mechanisms involved in AD.

Alzheimer Disease↗

PubMed2025

Integrative proteomic analysis provides novel therapeutic insights for etiological subtypes of diabetes.

AIMS: Type 2 diabetes (T2D) is a highly heterogeneous disease characterised by subtypes with variations in aetiology, disease progression, and risk of complications. However, potential drug targets for these subtypes have not been explored. This study aims to investigate potential drug targets by integrating proteomics. MATERIALS AND METHODS: Summary-level data of circulating proteins were extracted from the UK Biobank and the deCODE Health Study. Genetic associations with five diabetes subtypes were obtained from Swedish All New Diabetics in Scania and Malmö Diet and Cancer cohort, including severe autoimmune diabetes (SAID), severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD), mild obesity-related diabetes (MOD), and mild age-related diabetes (MARD). The associations between circulating proteins and diabetes subtypes were assessed through Mendelian randomisation, followed by multiple sensitivity and colocalization analyses. Additionally, tissue-specific, pathway and functional enrichment analysis, assessment of protein druggability, and the protein-protein interaction (PPI) networks were used to further explore biological mechanisms and therapeutic potential. RESULTS: Genetically predicted levels of 2, 2, 9, 3, and 5 circulating proteins were associated with SIRD, SIDD, MARD, MOD, and SAID, respectively. Colocalization analyses further revealed links between GRN with MARD/SIRD, LILRB5 with SIDD/MARD, CR1 with MARD, TNFSF12 with MOD, and DAPK2 with SAID. Enrichment analysis suggested that these proteins were mainly enriched in blood and adipose tissues and involved in immune and inflammatory related pathways. PPI analysis revealed GRN, TNFSF12, and DAPK2 are associated with known T2D targets. CONCLUSIONS: Our study identified several potential drug targets for different subtypes of diabetes using an integrated genetic approach, yielding new insights for precision medicine of diabetes.

Humans↗