SEARCH · Search PubMed

Results for “Co-expression analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 19 recordsLinked to original sources

PubMed2026

Integrated bioinformatics analysis reveals cross-talking hub genes and therapeutic agents between sepsis and acute myocardial infarction.

BACKGROUND: Sepsis and acute myocardial infarction (AMI) are two significant diseases that may share overlapping etiological mechanisms. This study aims to systematically identify core genes common to both conditions and to explore their potential as therapeutic targets and drug candidates through an integrative analysis of clinical data and bioinformatics. METHODS: The AMI dataset was obtained from the GEO database, and RNA sequencing data were collected from blood samples of patients with sepsis at our hospital. Common genes were identified using differential expression gene analysis (DEG) and weighted gene co-expression network analysis (WGCNA). Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were performed. A protein-protein interaction (PPI) network was constructed, and hub genes were identified using the MCC/Degree algorithm. Diagnostic value was assessed via receiver operating characteristic curve analysis. Immune infiltration patterns, single-cell sequencing data, and molecular docking simulations were employed to evaluate immune relevance and identify potential therapeutic compounds. RESULTS: A total of 417 genes were identified between sepsis and AMI, with enrichment analysis revealing significant involvement in inflammatory responses. Three hub genes-JAK2, MYD88, and TIMP1-were selected for further investigation. ROC curves confirmed their strong diagnostic performance for both diseases. Immune infiltration analysis showed that these core genes were significantly correlated with the infiltration levels of various immune cell types. Molecular docking indicated that quercetin exhibited stable binding affinity with the proteins encoded by these genes. qPCR validation further confirmed the upregulation of these three genes, supporting the anti-inflammatory effects of quercetin as a potential targeted therapy. CONCLUSION: JAK2, MYD88, and TIMP1 were identified as shared core genes in sepsis and AMI. These genes not only serve as potential diagnostic biomarkers but also offer novel targets for developing common therapeutic strategies for both conditions. Furthermore, quercetin emerges as a promising candidate for targeted treatment.

Humans↗

PubMed2026

New insights into diagnostic values and mechanisms of ferroptosis associated with immune infiltration in diabetic kidney disease.

The pathogenesis of diabetic kidney disease (DKD) is complex and closely related to ferroptosis and immune dysregulation, but the relevance is unclear. The present study investigates the potential mechanisms of ferroptosis-related genes (FRGs) in DKD and their relationship with the immune-inflammatory response. It searches for new diagnostic biomarkers to help diagnose and treat DKD. Four Gene Expression Omnibus (GEO) datasets, GSE30528, GSE30529 and GSE30122 as the test set, and GSE96804 for validation, were analyzed. FRGs were obtained from GeneCards, and 47 ferroptosis-related differentially expressed genes (FRDEGs) were identified by intersecting with DKD-related differentially expressed genes. Functional enrichment analyses, including Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Gene Set Enrichment Analysis and Gene Set Variation Analysis, revealed that these FRDEGs are primarily associated with ferroptosis, hypoxia response and immune inflammation. Subsequently, the weighted gene co-expression network analysis (WGCNA) was employed to expand the ferroptosis-related gene network, and intersection of the 47 FRDEGs with key WGCNA module genes yielded 10 key genes. Based on the 10 key genes, the least absolute shrinkage and selection operator and support vector machine algorithms identified three hub genes [chemokine ligand 5 (CCL5), forkhead box C1 (FOXC1) and lactotransferrin (LTF)] for DKD diagnosis. Receiver operating characteristic curves confirmed their diagnostic value, with FOXC1 and LTF validated in the independent dataset. Immune infiltration analysis via CIBERSORT revealed eight immune cell types with significantly different infiltration levels between the DKD and control group in the integrated GEO datasets. Notably, both LTF and CCL5 showed a significant positive correlation with gamma delta T cells (γδT). Quantitative PCR results confirmed differential expression of the three hub genes in the DKD group, with elevated expression observed in DKD mice following intervention with rosiglitazone and hyperoside.

bioinformatics analysis↗

PubMed2026

Pan-cancer single-cell atlas of immunotherapy response identifies ZNF385A as a regulator of immune evasion in small cell lung cancer.

Although immune checkpoint inhibitors (ICIs) have revolutionized the treatment landscape of solid tumors, response rates in patients with small cell lung cancer (SCLC) remain limited, and acquired resistance is highly prevalent. The underlying mechanisms of this immunotherapy resistance remain to be fully elucidated. Clinically, SCLC typically manifests as an "immune-cold" tumor, characterized by a low abundance of CD8+ T cell infiltration and the rare formation of tertiary lymphoid structures (TLS). While DNA damage repair (DDR) is closely linked to innate immune responses, how DDR networks orchestrate the SCLC immune microenvironment remains obscure. In this study, we integrated single-cell transcriptomic data (comprising 344,447 high-quality cells) from six cancer types (BCC, CRC, HCC, HNSCC, iCCA, and SCLC). Our comparative analysis revealed a fundamental depletion of TLS-associated cellular subpopulations (e.g., CXCL13+ CD8+ T cells, HLA-DRB5+ B cells, and CXCL9+ dendritic cells) in SCLC, which was significantly correlated with aberrant DDR activity. Through high-dimensional weighted gene co-expression network analysis (hdWGCNA), we identified ZNF385A as the core hub gene within the DDR-associated module. ZNF385A is highly expressed in SCLC and is associated with poorer prognosis. In vitro, ZNF385A depletion suppressed SCLC cell proliferation and induced apoptosis, accompanied by R-loop accumulation and activation of cGAS-STING signaling, indicating a potential link between ZNF385A, genomic stability and tumor-intrinsic innate immune signaling. Collectively, these findings identify ZNF385A as a potential regulator associated with TLS deficiency and immune evasion in SCLC.

Immunotherapy resistance↗

PubMed2026

Identifying Co-Expressed lncRNAs Correlated With Traits of Interest in an Animal Model for Metabolic Diseases in Humans.

Nutrigenomics investigates how nutrients modulate gene expression. Among them, fatty acids (FA) play important roles in regulating gene transcription, while long non-coding RNAs (lncRNAs) may be associated with gene regulation and metabolic diseases. This study aimed to analyze the hepatic transcriptome of pigs, a species frequently used as a model for nutrigenomic studies, to identify novel lncRNAs and their potential target genes in response to diets containing different sources of FA. Seventy-two pigs were fed four diets supplemented with 1.5% soybean oil (control), 3% canola oil, 3% fish oil, and 3% soybean oil. RNA sequencing of liver samples was performed to identify novel lncRNAs. Weighted Gene Co-expression Network Analysis (WGCNA) was used to identify modules associated with phenotypic traits related to lipid metabolism and inflammation. Functional enrichment analyses were then conducted to annotate genes within these modules using Gene Ontology (GO) terms and to assess overlap with Quantitative Trait Loci (QTL). The results revealed 106 novel lncRNAs potentially regulating genes associated with lipid metabolism and immune responses in pigs fed diets with different FA sources. These findings enhance understanding of the regulatory role of lncRNAs in pigs and reinforce their relevance as models for human metabolic diseases.

Animals↗

PubMed2026

ACE2 and Parkinsonism‑related bone metabolic alterations: signaling pathways and hub gene analysis.

Clinical co-occurrence of Parkinson's disease (PD) and age-related bone loss in elderly patients has garnered increasing attention, yet its molecular mechanisms remain incompletely elucidated. This study used an 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-induced PD model in Ace2-/y mice to investigate the regulatory mechanisms of bone-brain axis-related genes and signaling pathways. Behavioral tests assessed motor and non-motor symptoms. Immunohistochemistry, Western blot, and histopathological staining analyzed dopaminergic neuron activity, microglial activation, and bone metabolic abnormalities. GEO dataset transcriptomics and weighted gene co-expression network analysis (WGCNA) identified key hub genes, with receiver operating characteristic (ROC) curves evaluating their diagnostic value in public single-disease transcriptome data. MPTP significantly exacerbated motor dysfunction and depression-like behaviors; Ace2 deletion lowered total Wnt, β-catenin, BMP and IGF-1 protein abundance alongside reduced phosphorylation ratios of their downstream kinases in brain and bone, while upregulating RANKL/RANK/OPG-associated inflammatory mediators, accompanied by elevated total α-synuclein, Casp3 and Bax protein levels. The parallel reduction of these signaling proteins only suggests potential perturbation of related cascades; WGCNA identified 10 hub genes (e.g., DNM1, OCRL, OPA1), whose dysregulation was linked to synaptic dysfunction and inflammation. ROC analysis based on single-disease datasets showed high diagnostic accuracy for PD and `osteoporosis (OP) (AUC: 0.683-0.981), with core genes influencing synaptic, MAPK, Rap1, and Ras pathways. These preclinical findings indicate that Ace2 deficiency is associated with concurrent pathological abnormalities in the brain and transient bone metabolic disturbance under short-term MPTP treatment in growing young male mice; coordinated dysregulation of shared signaling pathways was observed in the two tissues, consistent with a potential bone-brain axis pathological phenotype, though causal bidirectional tissue cross-talk cannot be confirmed in the current experimental design, providing candidate targets that warrant further validation.

Animals↗

PubMed2026

Uncovering molecular regulatory networks of low-temperature stress response in Trachinotus ovatus via integrated transcriptome and metabolome analyses.

Golden pompano (Trachinotus ovatus) is one of the most economically important marine fish species in China. It is susceptible to low-temperature stress, which significantly challenges its production and supply. Nevertheless, study on the regulatory mechanisms underlying low-temperature stress responses in golden pompano remains limited. Here, we firstly performed a time-series transcriptome analysis to reconstruct dynamic response patterns under low-temperature stress in golden pompano. Transcriptome profiling identified common differentially expressed genes (DEGs), including fos, hlf, and hmgb1, as well as condition-specific DEGs across distinct low-temperature stress groups. Based on cluster analysis, all DEGs were classified into five distinct expression patterns, reflecting diversified regulation of expression in golden pompano during low-temperature stress. Furthermore, condition-specific regulatory modules were explored via weighted gene co-expression network analysis (WGCNA), highlighting that the two module hub genes, serbf2 and lipc, might respond to low-temperature stress by regulating the lipid catabolic process. Subsequently, untargeted metabolomic analysis revealed that glycerophospholipid metabolism was a significantly enriched common pathway, highlighting its crucial role in mediating the response to low-temperature stress. Finally, by integrating transcriptomic and metabolomic analyses, a gene-metabolite interaction network associated with glycerophospholipid metabolism under low-temperature stress was established. These findings underscore the significance of multiple candidate genes and glycerophospholipid metabolism in golden pompano's response to low-temperature stress, thereby laying a solid molecular foundation for the development of low-temperature-tolerant fish strains.

Animals↗

PubMed2026

Transcriptomics reveals species-specific adaptive strategies to calorie restriction in two Argopecten scallops with distinct lifespans.

Calorie restriction (CR) is a well-established non-genetic intervention for lifespan extension in multiple model organisms. Seasonal food shortage in cold and temperate seas may mimic CR, inducing in bivalves a response similar to that in vertebrates and thereby prolonging life expectancy. However, the relationship and the mechanism underlying the food availability and lifespan in bivalves remain largely unexplored. Two closely related scallop species the short-lived warm-water Argopecten irradians (lifespan <2 years) and the longer-lived cold-water Argopecten purpuratus (7-10 years) provide an ideal comparative system to investigate species-specific adaptive strategies. In this study, we subjected both species to CR for 30 and 56 days and performed comparative transcriptomic profiling, weighted gene co-expression network analysis (WGCNA), and physiological assays to elucidate their distinct molecular responses. Transcriptomic analysis revealed that A. purpuratus exhibited substantially more DEGs than A. irradians at both time points under CR, with both species showing downregulation of metabolic pathways but to different extents. A. irradians mounted an early nutrient-sensing response at 30 days (IGF1R, PIK3R3, INSR suppression), indicating acute sensitivity to limitation; by contrast, A. purpuratus displayed delayed FoxO activation at 56 days, along with its downstream effectors NFKBIA, CREB3L4, and SMAD4, suggesting a gradual adaptive program may link to its extended lifespan. WGCNA identified three negatively correlated modules in each species, with coral2 being the most prominent in A. irradians and darkolivegreen in A. purpuratus. The former was dominated by ciliary motility genes, whereas the latter featured coordinated repression of oxidative phosphorylation. Additionally, both species exhibited conserved suppression of mTOR/S6K growth signaling and activation of cellular maintenance programs. Collectively, these findings expand the understanding of CR-mediated longevity regulation in bivalves and provide candidate gene resources for future functional studies and breeding programs.

Pectinidae↗

PubMed2026

Multi-omics analysis reveals coordinated epigenetic dysregulation in atrazine-induced dopaminergic neurotoxicity.

Atrazine (ATR), a widely used triazine herbicide, has been linked to neurotoxicity, yet the epigenetic mechanisms underlying its dopaminergic effects remain unclear. This study investigated whether coordinated miRNA dysregulation and DNA methylation alterations contribute to ATR-induced Parkinson's disease (PD)-like neurotoxicity. Male Sprague-Dawley rats were administered ATR (50 mg/kg/day) for 90 days, resulting in motor and cognitive deficits with dopaminergic dysfunction, including increased α-synuclein and reduced tyrosine hydroxylase expression. Small RNA sequencing identified 72 differentially expressed miRNAs in the substantia nigra, enriched in PI3K-Akt, MAPK, and Ras signaling pathways. In a cohort of six PD patients and six matched controls, genome-wide DNA methylation profiling revealed 4694 differentially methylated positions, predominantly hypomethylated, with overlapping enrichment in neuronal signaling pathways. Weighted gene co-expression network analysis identified a PD-associated module strongly correlated with disease status (r = -0.95, P < 0.001). Multi-omics integration identified CASP3 as a central hub gene. External validation supported CASP3 relevance in PD (AUC = 0.833), and molecular docking suggested potential ATR-CASP3 interaction. Further analysis predicted upregulated miR-3552 as a potential upstream regulator of CASP3. These findings indicate that ATR-induced neurotoxicity may be mediated through the miR-3552/CASP3 signaling axis, ultimately regulating apoptosis and contributing to neurodegeneration.

Animals↗

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

Dynamic evolution of chaperone-mediated autophagy is associated with tumor microenvironment remodeling and prognostic stratification in lung adenocarcinoma: insights from single-cell transcriptomics, ensemble machine learning, and experimental validation.

BACKGROUND: Lung adenocarcinoma (LUAD) shows prognostic heterogeneity, and tumor-node-metastasis (TNM) staging is limited for individualized management. Chaperone-mediated autophagy (CMA) maintains proteostasis, but its role during adenocarcinoma in situ (AIS)-minimally invasive adenocarcinoma (MIA)-invasive adenocarcinoma (IAC) progression remains unclear. METHODS: Single-cell RNA sequencing (scRNA-seq) data from GSE189357 and bulk transcriptomes from The Cancer Genome Atlas (TCGA)-LUAD and Gene Expression Omnibus (GEO) cohorts were integrated. CMA activity, cell-cell communication, weighted gene co-expression network analysis (WGCNA), tumor-normal differential expression, machine-learning survival modeling, tumor microenvironment (TME) features, drug sensitivity, and EPC1 function were analyzed. RESULTS: CMA-high tumor epithelial cells increased from AIS (58.1%) to MIA (65.7%) but declined in IAC (44.4%; p < 0.001). CMA-low cells preferentially received fibroblast-derived extracellular matrix cues. A CMA-negatively correlated module identified 69 core genes. Random survival forest (RSF) performed best among 117 machine-learning combinations (mean concordance index > 0.873). High-risk patients had worse survival across cohorts, and the risk score was independently associated with overall survival (hazard ratio = 16.013, 95% confidence interval: 9.579-26.768, p < 0.001). High-risk tumors showed proliferative activation and M0 macrophage enrichment, whereas low-risk tumors showed stronger immune-related signaling. EPC1 overexpression suppressed malignant phenotypes in A549 cells. CONCLUSION: CMA dynamics are associated with stromal and immune remodeling during LUAD progression. A CMA-based model provides robust prognostic stratification and may offer a basis for future TME-guided studies.

Chaperone-mediated autophagy↗

PubMed2026

Analysis of end-stage renal disease mediated by cuproptosis-related genes.

OBJECTIVE: The complex pathophysiological mechanism of end-stage renal disease (ESRD) has not been fully understood. Cuproptosis is a newly discovered type of programmed cell death. Therefore, this study attempts to clarify the relationship between cuproptosis-related genes (CRGs) and the phenotype of ESRD. MATERIALS AND METHODS: The National Center for Biological Information Gene Expression Omnibus database was applied to obtain the GSE37171 dataset comprising whole-genome microarray analysis of peripheral blood samples. A 3 : 1 case-control design was employed with 75 ESRD patients and 20 healthy controls who were frequency-matched for age, sex, and ethnicity. Based on differentially expressed genes (DEGs) and genes related to cuproptosis, CRGs were identified. Thereafter, we explored two different subpopulations based on the cuproptosis gene and analyzed their expression and immune infiltration. Genes specific to the CRG cluster were identified through the weighted gene co-expression network analysis algorithm, and the best prediction model was determined and verified by four machine learning methods. RESULTS: The study identified 14 differentially expressed CRGs, among which ATP7B, SLC31A1, LIAS, LIPT1, DLD, MTF1, CDKN2A, DBT, and DLST had relatively high expression levels in the ESRD samples. Compared with the control group, expression levels of FDX1, DLAT, PDHA1, PDHB, and GLS were significantly lower in the ESRD group, and CRGs played a key role in the regulation of immune infiltration in ESRD. Two cuproptosis-related molecular clusters were identified in the ESRD samples. Cluster2 was more correlated with the immune infiltration of ESRD. By analyzing the intersection points between CRG cluster and key genes of ESRD, a total of 888 specific DEGs were identified. Functional differences related to specific DEGs were further explored using gene set variation analysis. Five significant genes (SMC5, USP47, USP53, AGA, and DMXL1) were identified by the support vector machine model as key predictors for ESRD disease risk, achieving an area under the curve (AUC) of 1.00 in internal validation. However, external validation in independent cohorts is required prior to clinical application. Individual gene analysis showed an AUC > 0.81 in discriminating ESRD patients from healthy controls, and the expression of all 5 genes in ESRD patients was significantly lower than in the control group. CONCLUSION: This study clarified the relationship between CRGs and the phenotype of ESRD, analyzed their specific roles in the immune microenvironment, and obtained a predictive model, providing new insights for the study of its potential therapeutic targets.

Humans↗

PubMed2026

Network based approach identifies miR-145-3p as a central regulatory hub associated to the progression from localized to metastatic medullary thyroid carcinoma.

Medullary thyroid carcinoma (MTC) is a neuroendocrine tumor originating from calcitonin producing C-cells and accounts for 1-5% of thyroid cancers. Total thyroidectomy is curative in localized disease (N0), whereas lymph node metastases (N1) are associated with poorer prognosis. However, the molecular mechanisms driving the metastatic shift remain poorly understood. This study aimed to identify miRNA features linked to metastatic spread in MTC, focusing on the transition from N0 to N1. Co-expression networks were constructed for N0 and N1 tumors, and differential connectivity analysis was used to identify key miRNAs acting as regulatory hubs. Functional annotation of their target genes was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Reactome pathway analyses. Validation experiments were carried out in MTC cells to evaluate the effects of selected miRNAs on cell proliferation, survival, and MAPK pathway activation. Network analysis revealed distinct miRNA co-expression patterns between N0 and N1 tumors. Differential network analysis highlighted miR-145-3p as a central regulatory hub, exhibiting 29 altered co-expression changes and a marked loss of connectivity in N1. Target enrichment identified 59 validated genes, including key oncogenic drivers such as MYC, PTEN, BCL2, PIK3CA, AKT1, and MAPK7. In MTC cells, simultaneous inhibition of miR-145-3p together with its top co-expressed miRNAs increased proliferation and survival, and enhanced ERK phosphorylation, indicating MAPK pathway activation and a shift toward a more aggressive phenotype. In conclusion, this study identifies a miRNA regulatory hub centered on miR-145-3p that is associated with metastatic progression and highlights the value of network-based approaches in uncovering mechanisms of cancer dissemination. © 2026 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

MAPK signaling↗

PubMed2026

Unveiling the molecular basis of gonadal development: Multi-omics uncovers sex-related genes and steroid pathways in Sinonovacula constricta.

The razor clam Sinonovacula constricta is an economically important cultured mollusk in China, but the molecular mechanism of its gonadal development and sexual differentiation remains unclear. This study integrated gonadal transcriptomic, proteomic, and metabolomic analysis to identify key sex-related molecules. Transcriptome analysis identified 2795 DELs and 6497 DEGs between sexes, including the sex-related genes Fem-1b, Fem-1c, GUCY1B2 and FAT4, as well as a regulatory network of 39 lncRNA-mRNA pairs involving Tektin-4, Ropporin-1, Histone H1, and FoxN4. Proteomic analysis revealed 3217 DEPs: Tektin family members, Ropporin-1 and Tssk proteins were upregulated in the testis, while histone H1 and FAT4 were upregulated in the ovary. Metabolomic analysis detected 409 DEMs, with uridine identified as a potential sex differential marker (upregulated in the ovary), and 23 gonadal development-related DEMs showed sex-specific upregulation. Integrative transcriptome-proteome analysis identified 1543 co-expressed DEGs/DEPs enriched in nucleosome assembly, oxidative phosphorylation, and carbon metabolism, including key sex-related genes AKAP14, Tektin/Tssk families, Histone H1, and FAT4. Transcriptome-metabolome integration identified 32 shared KEGG pathways (e.g., biosynthesis of unsaturated fatty acids, pyrimidine metabolism), while proteome-metabolome integration revealed 5 (positive ion) and 6 (negative ion) co-enriched pathways, with alanine, aspartate and glutamate metabolism and oxidative phosphorylation being functionally relevant to gonadal development. Collectively, these results reveal the molecular basis of gonadal development, highlight critical sex-related genes and steroid metabolic pathways, and provide valuable resources for future reproduction and breeding in S. constricta.

Animals↗

PubMed2026

Context matters: coordinated transcriptional regulation and root plasticity under multinutrient conditions.

Plants often encounter simultaneous imbalances in multiple nutrients, but the regulatory logic coordinating their responses remains poorly understood. We aimed to uncover shared transcriptional programs and regulatory nodes underpinning multinutrient adaptation in Arabidopsis thaliana roots. We analyzed publicly available RNA-seq datasets spanning 15 nutrient and beneficial element conditions using differential expression, co-expression network (WGCNA), and gene regulatory network analysis. Selected transcription factors (TFs) were validated via root phenotyping, suberin staining, and ionomic profiling under two-nutrient stress conditions. We identified a core set of 2050 genes responsive to multiple nutrient treatments, enriched for suberin biosynthesis, and structured into modular co-expression clusters. Eight prioritized candidate TFs (ARR10, GBF3, HHO5, NAC32, NF-YA3, NF-YB2, SARD1, and WRKY33) were shown to modulate root system architecture under specific nutrient combinations. WRKY33 and NF-YB2, in particular, regulated nutrient-responsive suberin deposition and ionomic plasticity. These findings reveal suberin remodeling as a shared downstream process in multinutrient responses and suggest that plasticity is not a fixed trait but a modular, polygenic, and context-dependent outcome. Repurposed TFs with pleiotropic functions coordinate structural and physiological traits, providing regulatory entry points for improving nutrient resilience.

Plant Roots↗

PubMed2026

Integrative multi-omics reveals a fibroblast-centered, ZFHX3-prioritized regulatory framework linking sick sinus syndrome and atrial fibrillation.

OBJECTIVE: To define shared genetic and multi-scale mechanisms underlying comorbidity between sick sinus syndrome (SSS) and atrial fibrillation (AF). METHODS: We integrated genome-wide association study (GWAS) summary statistics for SSS and AF with Genotype-Tissue Expression (GTEx) expression and splicing quantitative trait loci (eQTL/sQTL), atrial single-cell and spatial transcriptomics, and epigenomics. We identified trait-relevant tissues and pathways, prioritized shared cell types, quantified genome-wide and local genetic sharing, detected joint loci by cross-trait meta-analysis, and linked loci to regulatory programs via colocalization and cell-prioritized co-expression networks. RESULTS: Both traits showed strongest enrichment in cardiac tissue, especially Heart Atrial Appendage. Fibroblasts from the left atrial appendage were consistently prioritized as the key shared cell population. SSS and AF displayed significant positive genome-wide genetic correlation, with multiple locally shared regions, including six major loci. Cross-trait meta-analysis identified eight joint-phenotype SNPs implicating four susceptibility genes. ZFHX3 was the leading tissue-cell-gene candidate, acting as a hub in fibroblast co-expression modules and colocalizing with cardiac regulatory signals. CONCLUSION: Shared liability for SSS and AF is highly tissue- and cell-specific, converging on regulatory networks in atrial appendage fibroblasts, with ZFHX3 serving as a central mechanistic and biomarker node.

Humans↗

PubMed2026

Brachyury expression levels predict lineage potential and axis-forming ability of in vitro-derived neuromesodermal progenitors.

Neuromesodermal progenitors (NMPs) produce the spinal cord and musculoskeleton in the elongating anterior-posterior axis. In vivo, NMPs possess dual potency, coinciding with regions co-expressing SOX2 and Brachyury (TBXT). In vitro, SOX2/TBXT co-expressing cells can be produced from pluripotent cells and, like their in vivo counterparts, can produce neural tube and somitic mesoderm. However, the functional characteristics of in vitro SOX2/TBXT co-expressing cells remain unclear, confounding comparisons with in vivo data. To address this, we developed a dual Sox2/Tbxt reporter mouse ESC line. SOX2/TBXT reporter-positive cells emerge in vitro from pluripotent populations with dynamics that mirror their appearance in the embryo. Purified SOX2/TBXT co-expressing populations can differentiate towards neurectoderm or mesoderm, including lateral mesoderm upon BMP stimulation. In gastruloids, quantitative live imaging shows that WNT or NOTCH inhibition rapidly leads to downregulation of TBXT expression and diminished axial extension. We show that clonally plated SOX2/TBXT co-expressing cells are bipotent NMPs that can also self-propagate. By combining clonal analysis with mathematical inference, we identify two thresholds of TBXT and/or SOX2 expression, switching clonal output from neural- to mesoderm-biased, and from mesoderm-biased to mesoderm-specified. Image analysis of embryonic NMPs supports a model whereby SOX2 and TBXT independently influence neuromesodermal differentiation. Thus, this Sox2/Tbxt double reporter cell line highlights unsuspected heterogeneity in NMPs, and together with image analysis of embryonic SOX2/TBXT levels, challenges the assumption that neuromesodermal fate choice is primarily governed by mutual antagonism between SOX2/TBXT.

Animals↗

PubMed2026

Hepatic metabolic adaptation to endurance exercise: temporal and sex differences by multiomics integration and validation.

BACKGROUND: Although endurance exercise benefits liver health, sex-specific adaptive trajectories remain unclear. This study mapped dynamic liver adaptation in males and females during prolonged training and identified underlying molecular programs. METHODS: Using publicly available time-resolved liver multi-omics data generated by the Molecular Transducers of Physical Activity Consortium (MoTrPAC), we established a computational pipeline for differential analysis of transcriptomic, proteomic, phosphoproteomic, and metabolomic data with FDR correction, followed by FGSEA pathway enrichment. Kinase activities were inferred through ortholog mapping and PhosphoSitePlus. Cross-omics co-expression networks were constructed using WGCNA and topological overlap to link omics features with physiological phenotypes. For experimental validation, liver tissues were collected from endurance-trained Sprague-Dawley rats, and key nodes were confirmed by Western blotting, qRT-PCR, and immunofluorescence/immunohistochemical staining. Public scRNA-seq data were further integrated to map multi-omics signals to single-cell resolution and assess functional changes in specific cell types. RESULTS: The hepatic response to exercise stress was stage-specific, shifting from early transcriptional activation to later proteomic and metabolic remodeling. Multi-omics integration revealed distinct sex-associated adaptive trajectories: males were more strongly associated with energy metabolism, redox-related programs, and amino acid/organic acid catabolism, whereas females showed prominent membrane lipid remodeling, proteostasis -related programs, and mitochondrial/ribosomal translational features. Single-cell analysis showed that tissue remodeling occurred without major lineage turnover, instead involving altered communication among pre-existing cell communities. Validation of PPP1R3G identified a protein-dominant exercise-responsive marker, supporting the contribution of post-transcriptional or protein-level regulation. CONCLUSIONS: Hepatic adaptation to endurance stress follows a cross-omics evolutionary pattern with sex-specific reprogramming of energy supply and homeostatic maintenance. This time-resolved framework clarifies how exercise improves liver function and supports sex-oriented metabolic interventions and therapeutic target discovery.

Animals↗

PubMed2026

Dynamic and non-additive gene regulation shapes maize responses to simultaneous salt and cold stress.

Salt and cold stresses often occur together in nature and severely impact crop productivity, yet their transcriptional regulation remains poorly understood. Here, we conducted a time-series transcriptomic analysis of maize under salt, cold, and their combination at 0, 6, 12, and 24 h. Differential expression analysis revealed dynamic, condition-specific gene responses grouped into eight distinct temporal patterns. Promoter motif analysis of genes within each pattern identified 5-39 significantly enriched motifs, with over 40% lacking known counterparts, suggesting the involvement of previously uncharacterized cis-regulatory elements in stress-responsive transcriptional regulation. By comparing combined stress responses to the sum of single-stress effects, we found that about 74% of DEGs showed non-additive patterns, suggesting that combined stress triggers a distinct transcriptional program. Evolutionary analysis showed that additive DEGs tend to be more recently evolved, subject to weaker purifying selection, and enriched in transposed duplications, contrasting with the stronger constraint observed in non-additive DEGs. WGCNA identified 24 co-expression modules, among which 65 hub DEGs were detected in modules significantly correlated with specific stress conditions. Furthermore, we reconstructed 228, 20, and 200 sequential transcription factor cascades spanning 6 h, 12 h, and 24 h under cold, salt, and combined stress, respectively, with no cascade shared across all three conditions. Together, these results reveal that maize responses to combined salt and cold stress are largely non-additive and temporally dynamic, with distinct evolutionary patterns underlying different response types, offering insights and candidate regulators for enhancing crop stress resilience.

Zea mays↗