SEARCH · PubMed Health

Results for “Gene co-expression network”

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

Construction of molecular signatures based on the co-expression network of NECSO-related gene TRPM4 and its prognostic value in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) demonstrates significant prognostic variability that is not entirely accounted for by traditional staging systems. Necrosis by sodium overload (NECSO) is an emerging programmed cell death pathway, but its clinical relevance in HCC remains undefined. Therefore, this study aimed to identify TRPM4-associated core genes, develop and validate a prognostic signature, and investigate its relationship with the tumor immune microenvironment, tumor mutational burden, and single-cell expression patterns in HCC. METHODS: We integrated transcriptomic, clinical, and mutational datasets from The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) (n=421) and Gene Expression Omnibus (GEO) cohorts (n=115) to identify genes co-expressed with TRPM4-a key NECSO mediator-and those differentially expressed in HCC. A prognostic signature was developed using least absolute shrinkage and selection operator (LASSO)-Cox regression and validated through survival analysis, time-dependent receiver operating characteristic (ROC) curves, and multivariate Cox regression analysis. The immune landscape was characterized using CIBERSORT, somatic mutation data were used to calculate tumor mutational burden (TMB) and assess its correlation with the risk score, and single-cell RNA sequencing (scRNA-seq) resolved cell-type-specific expression patterns. RESULTS: From 294 TRPM4-associated core genes, we identified an 11-gene signature (BRSK1, MMP1, GRIN2D, GP6, MYOM2, N4BP3, CCDC112, TSEN54, MAP3K9, SPP1, B3GNT4) that independently predicted overall survival (OS) (hazard ratio =5.419, P<0.001) with areas under the curve (AUCs) of 0.779, 0.693, and 0.701 at 1, 3, and 5 years. These values were superior or comparable to conventional clinicopathologic variables after direct comparison. High-risk patients exhibited an immunosuppressive microenvironment, characterized by enrichment of M0 macrophage, a higher M2/M1 ratio (P<0.001) and distinct immune checkpoint profiles. When integrated with TMB, the prognostic stratification was further refined: high-TMB/high-risk patients had poorest outcomes (median OS, 15.3 months), while low-TMB/low-risk patients had the most favorable survival (median OS, 68.7 months). Single-cell analysis revealed that MMP1 was induced in cancer-associated fibroblasts (CAFs) and SPP1 was downregulated in macrophages, single-cell risk scores confirmed TAFs and macrophages as the main contributors to the prognostic model. CONCLUSIONS: The TRPM4-centered 11-gene signature provides robust and independent prognostic stratification in HCC by integrating immune, mutational, and single-cell features. This signature serves as a potential tool for prognostic evaluation and may help inform immunotherapeutic strategies for HCC.

Hepatocellular carcinoma (HCC)↗

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↗

PubMed2025

Multi-Omics insights into OsZFP252-OsGA20ox5 mediated drought tolerance in rice through stomatal and vascular regulation.

Rice growth is highly dependent on water availability, and drought stress significantly impacts its entire life cycle. However, previous studies lack systematic investigations into drought-responsive candidate genes across the full life cycle of rice. This study integrates transcriptomic and phenotypic data from two rice lines, IR64 (drought-sensitive) and DK151 (drought-tolerant), under varied environmental conditions at distinct growth stages. Using k-means clustering, 13 369 genes were categorized into 17 distinct expression patterns, revealing drought-responsive genes specifically upregulated or downregulated under drought stress. Weighted co-expression network analysis (WGCNA) further identified four gene modules strongly correlated with drought-related phenotypes, co-localizing 2859 drought-responsive genes through both approaches. Proteomics and metabolomics were supplemented at the booting stage, where phenotypic and transcriptomic differences under drought were most pronounced. Integrated omics results demonstrate gibberellin (GA) and abscisic acid (ABA) pathways play a key role during drought tolerance in rice, and 79 high-confidence drought-resistant candidate genes were prioritized from the 2859 drought-responsive genes. Among these, Gibberellin 20-oxidase 5 (OsGA20ox5) was identified as a key negative regulator of drought tolerance. Furthermore, the transcription factor zinc finger protein 252 (OsZFP252) directly binds to the OsGA20ox5 promoter, repressing its expression and enhancing ABA biosynthesis, thereby improving drought tolerance by increasing stomatal closure and expanding vascular bundle water transport capacity. Notably, the drought-tolerant haplotype 2-4 (Hap2-4) of OsGA20ox5 provides valuable insights for drought-resistant breeding.

Oryza↗

PubMed2026

Evidence for dual pathways of Tc1/mariner domestication in Drosophila.

BACKGROUND: The domestication of transposable elements is a key source of evolutionary innovation, yet the pathways by which their functional modules are repurposed by the host remain poorly understood. The Tc1/mariner superfamily is a widespread group of DNA transposons, but the prevalence and patterns of their domestication are underexplored. RESULTS: We performed a systematic genomic screen across 43 drosophilid species using stringent criteria for molecular domestication. This analysis identified five high-confidence, evolutionarily conserved genes derived from Tc1/mariner transposases. Phylogenetic and structural analyses suggest domestication via two distinct molecular pathways: co-option of the DNA-binding module and co-option of the catalytic domain. The DNA-binding module pathway includes CG4570, the previously known genes cag and toy (the latter fused with a homeodomain), and a lineage-restricted gene in the Drosophila obscura group that exhibits signatures of recent domestication. In contrast, the catalytic domain pathway is represented solely by CG14478. Structural modeling reveals that CG14478 protein preserves a canonical DDE endonuclease fold. Co-expression network analysis suggests potential cellular roles of these genes: CG14478 is linked to RNA/chromatin-related processes, CG4570 to cell cycle/chromosome functions, cag to ciliary and nuclear functions, and toy to neuronal development. CONCLUSIONS: This study establishes a stringent framework for identifying domesticated TEs, demonstrating that Tc1/mariner elements are co-opted via two distinct pathways: retention of either catalytic or DNA-binding modules. Our findings suggest that domestication is a dynamic continuum, ranging from recent, lineage-specific events to ancient, conserved genes, and underscore how genomic conflict with TEs can drive eukaryotic evolution and regulatory complexity.

Animals↗

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↗

PubMed2025

Genome imbalance modulates the expression of long non-coding RNAs in maize.

Genome imbalance, resulting from varying the dosage of individual chromosomes (aneuploidy), has a more detrimental effect than changes in complete sets of chromosomes (haploidy/polyploidy). This imbalance is likely due to disruptions in stoichiometry and interactions among macromolecular assemblies. Previous research has shown that aneuploidy causes global modulation of protein-coding genes (PCGs), microRNAs, and transposable elements (TEs), affecting both the varied chromosome (cis-located) and unvaried genome regions (trans-located) across various taxa. While long non-coding RNAs (lncRNAs) are important gene expression regulators, their roles in the context of genomic imbalance remain largely unexplored. In this study, we analyzed and compared the impact of aneuploidy and haploidy/polyploidy on lncRNA expression using RNA-seq data from maize mature leaf tissue. Our results indicate that cis-located lncRNAs are modulated from dosage compensation to a gene dosage effect, while trans-located lncRNAs exhibit trends ranging from an inverse effect to a positive correlation with chromosomal dosage. Remarkably, the ploidy series showed a lesser degree of lncRNA modulation. LncRNAs and TEs display a similar trend of inverse modulation but exhibit greater sensitivity to dosage changes compared to PCGs. The construction of cis-acting and trans-acting lncRNA co-expression networks indicates that lncRNAs likely function as dosage-sensitive regulators of gene expression under conditions of genomic imbalance. Overall, this study not only elucidates the dosage effect of plant lncRNAs but also serves as a valuable resource for exploring potential regulators of PCGs that play significant biological functions.

Zea mays↗

PubMed2026

Temporal DIA-MS proteomics reveals coordinated metabolic reprogramming associated with oil accumulation in oil palm mesocarp.

Oil palm (Elaeis guineensis Jacq.) is the most productive oil-bearing crop globally, yet the molecular basis of mesocarp development and lipid accumulation remains poorly understood. Ultra-deep data-independent acquisition mass spectrometry (DIA-MS) was applied to characterize proteome dynamics in two contrasting genotypes, seedless (KS) and thin-shelled (TS), across five developmental stages (P1-P5) spanning fruit development to mature oil accumulation. Phenotypic analysis revealed higher mesocarp proportion and oil content in KS during late maturation. A total of 137,615 peptides corresponding to 12,163 protein groups were identified, providing a temporal proteomic landscape of mesocarp development. Multivariate analysis indicated that developmental progression was the primary contributor to proteomic variation, whereas genotype-associated differences increased during lipid accumulation. Differentially abundant proteins were mainly associated with carbohydrate metabolism, photosynthesis, proteolysis, antioxidant responses, and lipid biosynthesis. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and KOG analyses suggested extensive remodeling of metabolic networks, including developmental changes in photosynthesis-associated proteins and increased representation of lipid-associated pathways during maturation. Weighted protein co-expression network analysis identified 17 modules associated with developmental progression and lipid accumulation, highlighting candidate proteins involved in carbon metabolism, energy production, and cellular protection. Genes encoding selected hub protein candidates were further examined by RT-qPCR. Biochemical analyses supported these proteomic patterns, showing increased acetyl-CoA availability, enhanced antioxidant enzyme activities (SOD, CAT, APX, and GR), improved GSH/GSSG balance, and reduced oxidative damage in KS. Together, these findings provide a temporal proteomic and biochemical framework for understanding genotype-associated differences in oil accumulation and identify candidate metabolic networks for functional studies.

Carbon metabolism↗

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

MetaCCI: meta cell-cell interaction inference and its application to CCIs characteristics of MDS.

MOTIVATION: Cell-cell interactions (CCIs) are fundamental to multicellular organisms and play crucial roles in diverse biological processes and disease mechanisms. Understanding CCIs is vital for deciphering disease pathogenesis and developing therapeutic strategies. Although numerous computational methods have been developed to infer CCIs from complex biological data, most existing approaches rely primarily on single-gene expression levels and ligand-receptor databases, often failing to capture the nuanced network-wide changes characteristic of disease states. RESULT: We propose MetaCCI, a novel computational strategy that integrates meta-information into CCI inference by extending the traditional gene expression-based analysis to a gene regulatory network framework. MetaCCI meticulously combines established ligand-receptor pairs with quantitative insights into gene behavior within complex gene networks, enabling the precise extraction of relevant targets for CCI inference. Subsequently, CCI inference was performed using an eigen cell co-expression network, providing a more holistic view of cell-cell communication. Monte Carlo simulations demonstrated that MetaCCI consistently outperforms existing methods in CCI inference. We applied MetaCCI to characterize cell-cell communication in Myelodysplastic Syndromes (MDS). Our results identified distinct interaction patterns in MDS compared with normal cell populations, specifically highlighting the loss of CCIs between "Dendritic cells and Hematopoietic precursor cells" and between "Dendritic cells and Hematopoietic multipotent progenitor cells" as characteristic features of MDS. Furthermore, FABP5, CD63, and HMGB1 were identified as MDS-specific markers. These findings suggest that diminished CCIs involving dendritic cells, hematopoietic precursor cells, and multipotent progenitor cells are pivotal to MDS pathogenesis. AVAILABILITY AND IMPLEMENTATION: The MetaCCI software is freely available at https://github.com/HeewonGitHub/MetaCCI. An archived version of the software and example datasets used in this study is available at Zenodo: https://doi.org/10.5281/zenodo.20101527.

Myelodysplastic Syndromes↗

PubMed2026

Quantum computing-assisted validation of a conserved macrophage suppression module shared by ASFV and PEDV.

BACKGROUND: African swine fever virus (ASFV) and porcine epidemic diarrhea virus (PEDV) differ in viral biology and cellular tropism, yet both pathogens suppress macrophage-mediated immune responses in pigs. OBJECTIVE: To identify a conserved macrophage suppression module shared by ASFV and PEDV and evaluate quantum computing as an independent framework for biological network validation. METHODS: Integrated analysis of publicly available GEO datasets (GSE231435 for ASFV and GSE306895) identified 471 shared downregulated genes. A network- and multi-omics-informed 20-gene core was selected and encoded as a 20-qubit modularity-based Quadratic Unconstrained Binary Optimization (QUBO) problem. Community detection was benchmarked using the Quantum Approximate Optimization Algorithm (QAOA) on both the IBM Quantum Aer simulator and the 156-qubit IBM Fez (Heron r2) quantum processor and compared with brute-force enumeration and simulated annealing. RESULTS: A conserved macrophage suppression module shared by ASFV and PEDV was identified. For the STRING protein-protein interaction network, QAOA at circuit depth p = 3 reproduced the brute-force optimum with an approximation ratio of 1.000. In contrast, performance progressively declined in the denser co-expression network with increasing circuit depth, consistent with noise accumulation under current Noisy Intermediate-Scale Quantum (NISQ) conditions. Multi-run consensus analysis identified stable hub genes, including MMP9 and SLA-DOA, as well as genes exhibiting variable community assignments. CONCLUSION: These findings reveal a conserved macrophage suppression module shared between ASFV and PEDV and demonstrate that quantum computing can serve as an independent validation framework for biologically meaningful host-response networks. Network topology emerged as a key determinant of QAOA performance on real NISQ hardware.

Animals↗

PubMed2026

Genes associated with translation and oxidative phosphorylation as components of the translational response in nodulated and water-restricted soybean.

BACKGROUND: Soybean primarily acquires nitrogen through symbiosis with nitrogen-fixing bacteria. Water deficit (WD) is a major stress limiting crop yield. Nodulation may enhance drought tolerance in legumes by modulating nitrogen and hormone metabolism, osmotic adjustment, and antioxidant defenses; however, the molecular basis underlying the differential WD responses between N-fix and N-fed plants remain unclear. Translational control of gene expression is a key regulatory mechanism during stress. RESULTS: We compared the transcriptome and translatome of soybean roots from N-fix and N-fed plants exposed to WD across four combined treatments. N-fix plants under WD exhibited more complex responses in terms of total differentially expressed genes (DEGs) compared to N-fed plants. This increased complexity was also evident among translationally regulated DEGs and differentially expressed transcription factors, whose involvement in WD responses of N-fix plants is novel. Co-expression network analysis identified modules associated with core biological processes encompassing nodulation, WD, and notably, their interplay was particularly prominent in Module 1, which was enriched in genes related to ribosomal protein synthesis and oxidative phosphorylation (OXPHOS). Guilt-by-Association analysis enabled the prediction of novel functions for differentially expressed, uncharacterized hub genes related to stress and/or nodulation responses. CONCLUSIONS: Translational regulation of genes involved in OXPHOS and translation initiation emerged as a central response in N-fix plants under WD. These findings reveal distinct molecular adaptations in N-fix soybean roots facing WD and highlight translational control as a key regulatory layer. We also identified promising candidate genes-including transcription factors and uncharacterized hub genes under translational regulation-that represent potential targets for improving drought tolerance in legumes once validated functionally.

Glycine max↗

PubMed2026

Modelling the effects of biological intervention in a dynamical gene network.

Cellular response to environmental and internal signals can be modeled by dynamical gene regulatory networks (GRN). In the literature, three main classes of gene network models can be distinguished: (1) non-quantitative (or data-based) models which do not describe the probability distribution of gene expressions; (2) quantitative models which fully describe the probability distribution of all genes co-expression; and (3) mechanistic models which allow for a causal interpretation of gene interactions. We propose two rigorous frameworks to model gene alteration in a dynamical GRN, depending on whether the network model is quantitative or mechanistic. We explain how these models can be used for design of experiment, or, if additional alteration data are available, for validation purposes or to improve the parameter estimation of the original model. We apply these methods to the Gaussian graphical model, which is quantitative but non-mechanistic, and to mechanistic models of Bayesian networks and penalized linear regression.

Gene Regulatory Networks↗

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

Integrated transcriptomic and functional characterization of Claudin-1 reveals its oncogenic and immunomodulatory roles in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) remains among the deadliest malignancies, driven by its invasive nature and lack of effective biomarkers. Disruption of the epithelial barrier, mediated by tight junction components, is a critical yet underexplored contributor to PDAC progression. Claudins, integral regulators of tight junction integrity, display altered expression across cancers, but their prognostic and immunomodulatory roles in PDAC remain unclear. We performed an integrative analysis of 177 RNA-Seq datasets from TCGA and GTEx to characterize Claudin family alterations in PDAC. Differential expression, copy number variation, methylation, and co-expression networks were analyzed alongside clinical and survival data. Prognostic significance was assessed using Kaplan - Meier and Cox regression analyses, while immune cell infiltration was examined using deconvolution algorithms. Functional validation of Claudin-1 was conducted in Capan-1 cells using CRISPR/Cas9 knockout, followed by proliferation, wound-healing, and Western blot assays. Ten Claudin genes were significantly dysregulated, with Claudin-1 and Claudin-4 frequently amplified and associated with advanced stage and poor survival. High Claudin-1 expression correlated with reduced immune infiltration, indicating an immune-excluded phenotype characterized by immune cells retained in the tumor stroma but largely absent from the tumor parenchyma. Claudin-1 knockout markedly inhibited proliferation, migration, and EMT, evidenced by downregulation of Snail and Slug and restoration of E-cadherin expression. This integrative transcriptomic and functional study identifies Claudin-1 as a key driver of PDAC aggressiveness and immune modulation. These findings establish Claudin-1 as a promising prognostic biomarker and therapeutic target for restoring epithelial integrity and counteracting immune evasion in pancreatic cancer.

Humans↗

PubMed2026

Natural leaf shape variation reveals diverse transcriptional targets of GmJAG1 during soybean leaf development.

The JAGGED transcription factor family regulates lateral organ development across angiosperms. In soybean (Glycine max Merr.), a D9H mutation in the EAR repression motif of GmJAG1 causes a narrow-leaflet phenotype and explains over 70% of phenotypic variance in leaf shape. Because this mutation does not affect the zinc finger DNA-binding domain, both alleles bind identical targets but differ in repressor recruitment. Previous studies mapped GmJAG1 binding sites, but the functional targets controlling leaf morphology are uncharacterized. Here, we used comparative transcriptomics across four soybean genotypes with contrasting leaf shape, spanning a developmental time series from shoot apex to mature leaf, and identified 1567 putative candidate target genes. GmJAG1 expression was confined to the shoot apex, yet 99.1% of candidate targets maintained differential expression throughout development. We found that neither Kip-Related Protein (KRP) cell cycle inhibitors nor Cyclin-Dependent Kinases (CDKs) showed differential expression despite binding evidence in Arabidopsis. However, D-type cyclins were upregulated in narrow-leaf genotypes suggesting that soybean GmJAG1 acts through cyclin-mediated rather than KRP-mediated cell cycle regulation described in Arabidopsis- a divergence in regulatory logic between the two species. Pathway analysis revealed enrichment of auxin (1.8-fold, P = 0.02) and salicylic acid (fourfold, P = 0.016) genes among JAG1D9H targets. Filtering by differential expression, binding data, phenotype correlation, and co-expression network membership identified 79 high-confidence targets, including orthologs of NPH3 (phototropin-mediated leaf flattening), MIK2 (cell wall integrity sensing), RD22 (ABA-responsive stress signaling), and SCL23 (GRAS transcription factor in bundle sheath development). These candidates provide targets for functional validation and breeding in legumes.

Glycine max↗

PubMed2025

Multi-Omics and Integrative Analytics in Natural Products Discovery.

Natural products (NPs) have long been an essential source of new bioactive compounds for drug discovery; however, traditional methods for screening and isolating these compounds can be slow and often yield diminishing returns. Fortunately, advanced multi-omics and computational approaches present powerful solutions to these challenges. This review highlights innovative methodologies that integrate metabolomics, genomics, transcriptomics, and proteomics with bioinformatics and analytical chemistry to accelerate NP discovery. For instance, untargeted metabolomics platforms like high-resolution liquid chromatography-tandem mass spectrometry (LC-MS/MS) and Global Natural Products Social (GNPS) molecular networking allow for comprehensive profiling of new compounds, while targeted isotope-labeling strategies enhance this process. Additionally, genome and metagenome mining tools such as antibiotics and secondary metabolite analysis shell (antiSMASH), Deep Biosynthetic Gene Cluster (DeepBGC), and Pipeline for Reconstructing Integrated Syntheses of Metabolites (PRISM) quickly identify biosynthetic gene clusters (BGCs) in both cultured and uncultured organisms, often using heterologous expression to validate products. Transcriptomic analyses, including RNA sequencing (RNA-seq), co-expression networks, and fluxomics, help clarify how pathways are regulated, while quantitative proteomics techniques like tandem mass tags/isobaric tags for relative and absolute quantitation (TMT/iTRAQ) and label-free methods, along with chemoproteomics approaches such as cellular thermal shift assay and thermal proteome profiling (TPP), uncover molecular targets and their mechanisms of action. This review also places significant emphasis on the role of artificial intelligence (AI) and machine learning (ML) in integrating multi-omics data, spanning activities from constructing gene-metabolite correlation networks to leveraging knowledge graphs and graph neural networks for data fusion and functional prediction. Finally, this review concludes by discussing the synergistic benefits of multi-omics for natural-product discovery, addressing current technical challenges, and exploring future directions toward high-throughput, intelligent data integration for next-generation NP research.

Biological Products↗

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↗

PubMed2025

A powerful framework for differential co-expression analysis of general risk factors.

MOTIVATION: Differential co-expression analysis (DCA) aims to identify genes in a pathway whose shared expression depends on a risk factor. While DCA provides insights into the biological activity of diseases, existing methods are limited to categorical risk factors and/or suffer from bias due to batch and variance-specific effects. We propose a new framework, Kernel-based DCA (KDCA), that harnesses correlation patterns between genes in a pathway to detect differential co-expression arising from general (i.e. continuous, discrete, or categorical) risk factors. RESULTS: Using various simulated pathway architectures, we find that KDCA accounts for common sources of bias to control the type I error rate while substantially increasing the power compared to the standard eigengene approach. We then applied KDCA to The Cancer Genome Atlas thyroid data set and found several differentially co-expressed pathways by age of diagnosis and BRAF mutation status that were undetected by the eigengene method. Collectively, our results demonstrate that KDCA is a powerful testing framework that expands DCA applications in expression studies. AVAILABILITY AND IMPLEMENTATION: KDCA is publicly available in the R package kdca. The package can be downloaded at https://github.com/ajbass/kdca.

Humans↗