Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Regulatory network”

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.

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

At least 19 recordsLinked to original sources
PubMed2026

Integrative computational analysis combining network pharmacology, regulatory network modeling, and molecular dynamics reveals the mechanisms of Quanshen compound in ITP.

UNLABELLED: Immune thrombocytopenia (ITP) is a hemorrhagic disorder caused by immune dysfunction. Quanshen Compound (QSC) is an in-house preparation developed by the Uyghur Hospital in Hotan Prefecture. This study primarily investigates and validates the potential pharmacological basis and mechanism of action of QSC in modulating immune thrombopoiesis. Based on the multi-database screening of the QSC and the related targets of ITP, the intersection was obtained to construct a protein-protein interaction (PPI) network and screen the core targets; the intersection targets were analyzed for gene ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using R packages; a component-target-pathway network was constructed to screen the key active components and their mechanisms of action. At the same time, the TF-mRNA-miRNA regulatory network of the core targets was constructed, and chromosome localization and subcellular localization analysis were performed; further, the binding stability of key components and core targets was verified through molecular docking and molecular dynamics simulation. A total of 227 potential target sites were screened out, among which TNF, IL6, AKT1, TP53 and IL1B were the core targets. The enrichment results indicated that these intersecting target sites mainly participated in inflammatory responses, immune regulation and hemostasis-related biological processes, and were significantly enriched in the PI3K-Akt signaling pathway, Toll-like receptor signaling pathway, Th17 cell differentiation and PD-1/PD-L1 signaling pathway. The core target TF-mRNA-miRNA regulatory network contained 184 nodes and 200 edges, suggesting that the core targets were subject to multi-level regulation. Molecular docking results showed that the main active components had good binding activity with the core targets, and molecular dynamics simulation further verified the stability of the complex. QSC may improve ITP through a multi-component, multi-target, and multi-pathway synergistic mechanism involving key targets such as TNF, IL6, AKT1, TP53, and IL1B, as well as the PI3K-Akt signaling pathway. These findings provide new insights into the potential therapeutic mechanisms of QSC against ITP and warrant further experimental validation. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at https://doi.org/10.1007/s40203-026-00718-0.

Immune thrombocytopenia↗
PubMed2025

Deep learning-based cell-specific gene regulatory networks inferred from single-cell multiome data.

Gene regulatory networks (GRNs) provide a global representation of how genetic/genomic information is transferred in living systems and are a key component in understanding genome regulation. Single-cell multiome data provide unprecedented opportunities to reconstruct GRNs at fine-grained resolution. However, the inference of GRNs is hindered by insufficient single omic profiles due to the characteristic high loss rate of single-cell sequencing data. In this study, we developed scMultiomeGRN, a deep learning framework to infer transcription factor (TF) regulatory networks via unique integration of single-cell genomic (single-cell RNA sequencing) and epigenomic (single-cell ATAC sequencing) data. We create scMultiomeGRN to elucidate these networks by conceptualizing TF network graph structures. Specifically, we build modality-specific neighbor aggregators and cross-modal attention modules to learn latent representations of TFs from single-cell multi-omics. We demonstrate that scMultiomeGRN outperforms state-of-the-art models on multiple benchmark datasets involved in diseases and health. Via scMultiomeGRN, we identified Alzheimer's disease-relevant regulatory network of SPI1 and RUNX1 for microglia. In summary, scMultiomeGRN offers a deep learning framework to identify cell type-specific gene regulatory network from single-cell multiome data.

Deep Learning↗
PubMed2026

gMISpy: integration of complex regulatory networks and genome scale metabolic models.

MOTIVATION: Genome-scale metabolic models lack explicit regulatory mechanisms, limiting their predictive accuracy for genetic interventions. Current methods for computing genetic Minimal Cut Sets either ignore regulatory networks entirely or use simplified acyclic representations that cannot capture regulatory feedback loops, ubiquitous features critical in cellular modeling. RESULTS: We developed gMISpy, a Python package that that enables efficient computation of genetic Minimal Intervention Sets (gMISs) in integrated genome-scale metabolic and regulatory networks. gMISpy incorporates cyclic regulatory logic into our previous computational framework using layered Boolean networks and BoNesis framework, resulting in a more accurate modeling of how regulatory interactions affect metabolic genes. Benchmarking across four different regulatory networks with Human-GEM showed consistent improvements in prediction accuracy, with Matthews correlation coefficient gains ranging from 2.50% to 14.42%. Validation against cancer data from DepMap and Project Score confirmed that cyclic integration reduces false positives and better captures biological vulnerabilities compared to acyclic approaches. AVAILABILITY AND IMPLEMENTATION: https://github.com/PlanesLab/cyclic-gMISpy.

Software↗
PubMed2026

GRNContext: an interactive web platform for contextualized gene regulatory networks visualization across human cancers.

SUMMARY: While current Gene Regulatory Network (GRN) databases provide comprehensive reference maps of potential interactions between transcription factors and target genes, they do not specify which regulatory interactions are active within specific biological contexts. This limitation is particularly critical in cancer, where transcriptional programs are inherently tissue-specific. To address this gap, we developed GRNContext, an interactive web platform designed for the visualization, exploration, and comparative analysis of gene regulatory networks contextualized across 33 cancer types from The Cancer Genome Atlas (TCGA). Our approach uses the TFLink human reference GRN as a starting point and integrates TCGA transcriptomic profiles to infer cancer-specific regulatory activity. Regulatory relevance was assessed using complementary machine learning and statistical methods, which were unified into a consensus score to prioritize and filter the most relevant candidate regulators for each target gene. By providing both curated context-specific GRNs and a user-friendly platform, GRNContext constitutes a comprehensive and accessible resource that supports mechanistic investigations, hypothesis generation, and translational research focused on transcriptional regulation in cancer. AVAILABILITY AND IMPLEMENTATION: GRNContext is supported by all major browsers and freely available on the web at https://apps.cienciavida.org/grncontext. It is implemented as a client-server web application featuring a FastAPI backend and a React frontend utilizing Cytoscape.js for interactive network visualization, all containerized via Docker for cross-platform compatibility.

Humans↗
PubMed2026

Construction of circRNA-miRNA-mRNA regulatory networks in the intestine of turbot (Scophthalmus maximus) following Vibrio anguillarum infection.

Circular RNAs (circRNAs) play pivotal roles in post-transcriptional regulation by acting as molecular sponges for microRNAs (miRNAs) within the competitive endogenous RNA (ceRNA) network. However, the regulatory mechanisms in teleost immune responses remain poorly understood. In this study, circRNA-miRNA-mRNA networks were investigated in turbot (Scophthalmus maximus) following Vibrio anguillarum infection to elucidate host-pathogen interactions. Through high-throughput sequencing of intestinal tissues, a total of 50 differentially expressed circRNAs (DE-circRNAs) (18 at 2 hpi, 16 at 12 hpi, 16 at 48 hpi), 212 DE-miRNAs (11 at 2 hpi, 70 at 12 hpi, 15 at 48 hpi), and 1774 DE-mRNAs were identified. Functional enrichment analyses (GO/KEGG) revealed significant associations with immune pathways, including the MAPK signaling pathway and gap junction. An integrated circRNA-miRNA-mRNA regulatory network was constructed, highlighting key interactions including novel_circ_0002573/DE-miR-27a-3p/FGB and novel_circ_0002423/novel_347/GNE, which may regulate inflammatory and antibacterial responses. The expression patterns of selected circRNAs, miRNAs and mRNAs were validated using qRT-PCR, confirming the reliability of the sequencing results. Importantly, fibrinogen beta chain (FGB) and CXCR4/CXCL12 signaling were identified as critical immune modulators. These findings provide insights of the ceRNA regulatory networks involved in teleost intestinal immunity and provide potential molecular targets for selective breeding of disease resistance in this species.

Animals↗
PubMed2027

Inferring Gene Regulatory Networks in Stem Cells: Methods and Applications.

Gene regulatory networks (GRNs) represent the complex interplay of transcription factors, regulatory elements, and target genes that orchestrate cellular identity and function, playing a crucial role in the differentiation and maintenance of stem cells. This chapter provides an overview of experimental and computational methodologies for inferring GRNs, with particular emphasis on single-cell approaches. We first review key experimental techniques for detecting transcription factor binding sites, chromatin accessibility, and DNA motifs, alongside essential databases that support GRN reconstruction. We then introduce computational inference methods that can be categorized into four principal frameworks: correlation-based approaches, regression and machine learning models, probabilistic and deep learning methods, and integrative or message-passing frameworks. To illustrate practical application, we present a case study applying the pySCENIC workflow to a peripheral blood mononuclear cell single-cell RNA sequencing dataset from mouse, demonstrating how regulon-based analysis can reveal cell-type-specific regulatory programs. This chapter aims to serve as a practical guide for researchers seeking to understand and implement GRN inference methodologies in stem cell biology and related fields.

Gene Regulatory Networks↗
PubMed2026

Causal assessment of Bayesian gene regulatory networks from single-cell transcriptomics.

Gene regulatory network (GRN) inference is an essential tool for revealing dysregulated relationships between genes in different cell types from single-cell transcriptomic (SCT) data. GRNs based on Bayesian networks (BNs) learned from SCT data can elucidate directed regulatory relationships representing complex disease mechanisms and their interplay through graphical modeling. However, software for learning BNs from SCT data is not widely available, nor is software for evaluating the BNs' structural accuracy in representing causal relationships between genes. Here, we describe the scstruc R package. This package provides a suite of BN structure learning algorithms specifically designed to handle SCT data, to evaluate the resulting networks based on the causal relationships they represent regardless of the availability of established molecular interaction networks, and to compare regulatory relationships between conditions. We demonstrated that scstruc can identify biologically relevant differential regulatory relationships between groups on a per-cell basis.

Bayesian networks↗
PubMed2024

Whole transcriptome sequencing analyses of islets reveal ncRNA regulatory networks underlying impaired insulin secretion and increased β-cell mass in high fat diet-induced diabetes mellitus.

AIM: Our study aims to identify novel non-coding RNA-mRNA regulatory networks associated with β-cell dysfunction and compensatory responses in obesity-related diabetes. METHODS: Glucose metabolism, islet architecture and secretion, and insulin sensitivity were characterized in C57BL/6J mice fed on a 60% high-fat diet (HFD) or control for 24 weeks. Islets were isolated for whole transcriptome sequencing to identify differentially expressed (DE) mRNAs, miRNAs, IncRNAs, and circRNAs. Regulatory networks involving miRNA-mRNA, lncRNA-mRNA, and lncRNA-miRNA-mRNA were constructed and functions were assessed through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. RESULTS: Despite compensatory hyperinsulinemia and a significant increase in β-cell mass with a slow rate of proliferation, HFD mice exhibited impaired glucose tolerance. In isolated islets, insulin secretion in response to glucose and palmitic acid deteriorated after 24 weeks of HFD. Whole transcriptomic sequencing identified a total of 1324 DE mRNAs, 14 DE miRNAs, 179 DE lncRNAs, and 680 DE circRNAs. Our transcriptomic dataset unveiled several core regulatory axes involved in the impaired insulin secretion in HFD mice, such as miR-6948-5p/Cacna1c, miR-6964-3p/Cacna1b, miR-3572-5p/Hk2, miR-3572-5p/Cckar and miR-677-5p/Camk2d. Additionally, proliferative and apoptotic targets, including miR-216a-3p/FKBP5, miR-670-3p/Foxo3, miR-677-5p/RIPK1, miR-802-3p/Smad2 and ENSMUST00000176781/Caspase9 possibly contribute to the increased β-cell mass in HFD islets. Furthermore, competing endogenous RNAs (ceRNA) regulatory network involving 7 DE miRNAs, 15 DE lncRNAs and 38 DE mRNAs might also participate in the development of HFD-induced diabetes. CONCLUSIONS: The comprehensive whole transcriptomic sequencing revealed novel non-coding RNA-mRNA regulatory networks associated with impaired insulin secretion and increased β-cell mass in obesity-related diabetes.

Mice↗
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↗
PubMed2025

Gene regulatory network structure informs the distribution of perturbation effects.

Gene regulatory networks (GRNs) govern many core developmental and biological processes underlying human complex traits. Even with broad-scale efforts to characterize the effects of molecular perturbations and interpret gene coexpression, it remains challenging to infer the architecture of gene regulation in a precise and efficient manner. Key properties of GRNs, like hierarchical structure, modular organization, and sparsity, provide both challenges and opportunities for this objective. Here, we seek to better understand properties of GRNs using a new approach to simulate their structure and model their function. We produce realistic network structures with a novel generating algorithm based on insights from small-world network theory, and we model gene expression regulation using stochastic differential equations formulated to accommodate modeling molecular perturbations. With these tools, we systematically describe the effects of gene knockouts within and across GRNs, finding a subset of networks that recapitulate features of a recent genome-scale perturbation study. With deeper analysis of these exemplar networks, we consider future avenues to map the architecture of gene expression regulation using data from cells in perturbed and unperturbed states, finding that while perturbation data are critical to discover specific regulatory interactions, data from unperturbed cells may be sufficient to reveal regulatory programs.

Gene Regulatory Networks↗
PubMed2025

Selective targeting of TBXT with DARPins identifies regulatory networks and therapeutic vulnerabilities in chordoma.

The embryonic transcription factor TBXT (brachyury) drives chordoma, a spinal neoplasm without effective drug therapies. TBXT's regulatory network is poorly understood, and strategies to disrupt its activity for therapeutic purposes are lacking. We developed designed ankyrin repeat proteins that block TBXT-DNA binding (T-DARPins). In chordoma cells, T-DARPins reduced cell cycle progression, spheroid formation, and tumor growth in mice and induced signs of senescence and differentiation. Transcriptomic and proteomic analyses identified gene networks involved in cell cycle regulation, embryonic cell identity, and interferon response and revealed features of regulome components, such as susceptibility to pharmacologic inhibition and the fine-tuning of TBXT downstream effectors through IGFBP3. Finally, we found high interferon signaling in chordoma cell lines and patient tumors, which was promoted by TBXT and associated with sensitivity to JAK2 inhibitors. These findings demonstrate the potential of DARPins for probing nuclear proteins to understand the regulatory networks of transcription factor-driven cancers, including entry points for therapies that warrant testing in patients.

Humans↗
PubMed2025

Genome-Wide Analysis of DtxR and HrrA Regulons Reveals Novel Targets and a High Level of Interconnectivity Between Iron and Heme Regulatory Networks in Corynebacterium glutamicum.

Iron is vital for most organisms, serving as a cofactor in enzymes, regulatory proteins, and respiratory cytochromes. In Corynebacterium glutamicum , iron and heme homeostasis are tightly interconnected and controlled by the global regulators DtxR and HrrA. While DtxR senses intracellular Fe2+, HrrSA is activated by heme. This study provides the first genome-wide analysis of DtxR and HrrA binding dynamics under varying iron and heme conditions using chromatin affinity purification and sequencing (ChAP-Seq). We revealed 25 novel DtxR targets and 210 previously unrecognized HrrA targets. Among these, metH, encoding homocysteine methyltransferase, and xerC, encoding a tyrosine recombinase, were bound by DtxR exclusively under heme conditions, underscoring condition-dependent variation. Activation of metH by DtxR links iron metabolism to methionine synthesis, potentially relevant for the mitigation of oxidative stress. Beyond novel targets, 16 shared targets between DtxR and HrrA, some with overlapping operator sequences, highlight their interconnected regulons. Strikingly, we demonstrate the significance of weak ChAP-Seq peaks that are often disregarded in global approaches, but feature an impact of the regulator on differential gene expression. These findings emphasize the importance of genome-wide profiling under different conditions to uncover novel targets and shed light on the complexity and dynamic nature of bacterial regulatory networks.

Corynebacterium glutamicum↗
PubMed2025

Decoding nitrogen uptake efficiency in maize and sorghum: insights from comparative gene regulatory networks.

Nitrogen (N) is an essential macronutrient for plant growth and yield, yet optimizing nitrogen use efficiency remains a challenge in agriculture. To better understand the regulatory basis of plant responses to N availability, we constructed a maize-specific nitrogen uptake efficiency gene regulatory network (mNUEGRN) comprising 1625 protein-DNA interactions (PDI) between 70 promoters and 301 transcription factors using enhanced yeast one-hybrid assays. We also projected a sorghum NUE GRN (spNUEGRN) based on maize orthologs and analyzed N-responsive subnetworks in both species using transcriptome profiling under N stress of early deprivation and recovery. Cross-species comparison with an existing Arabidopsis GRN revealed about 18% conserved interaction, corresponding to 11% of the mNUEGRN, particularly within the nitrate assimilation pathways. Notably, bZIP18 and bZIP30 emerged as central regulators in mNUEGRN, forming highly connected feed-forward loops (FFLs). From our time series data, we identified 19 236 and 23 864 differentially expressed genes in maize and sorghum, respectively. Gini correlation analysis uncovered 764 and 638 FFLs in mNUEGRN and spNUEGRN, respectively, of which 22 FFLs in maize and 35 in sorghum were identified in both leaf and root for each species. These FFLs may represent candidate regulatory motifs that contribute to modulating transcriptional responses under fluctuating N conditions, but their potential roles require further investigation. Together, our findings reveal evolutionarily conserved and species-specific regulatory strategies that mediate early N responsiveness, offering a foundation for engineering crops with improved NUE.

Sorghum↗
PubMed2026

Crosstalk mediators implicated in the Stevens-Johnson Syndrome through gene regulatory network analysis.

Stevens-Johnson syndrome (SJS) is a rare and severe mucocutaneous disorder often triggered by medications or infections. Our previous research identified that four key genes, Ikzf1, Ptger3, Mavs, and Tlr3 are involved in SJS susceptibility and the conjunctival epithelial innate immune response, demonstrating their role in regulating interferon-stimulated genes. However, the interplay among these regulatory factors remains unclear. This study aimed to elucidate the crosstalk mechanisms between the pathways regulated by these four genes in conjunctival epithelial cells. We constructed a comprehensive gene regulatory network using transcriptomic data from murine conjunctival epithelial cells under 16 distinct conditions, including polyI:C stimulation across wild-type, knockout, and transgenic backgrounds for the key genes. A targeted network analysis systematically identified numerous candidate genes mediating the crosstalk between the regulatory pathways initiated by Ikzf1, Ptger3, Mavs, and Tlr3. The identified candidates suggest the involvement of diverse signaling pathways previously unlinked to SJS pathology. Our findings suggest that the pathogenesis of SJS may arise not from the dysfunction of isolated genes but from the disruption of a balance maintained by intricate pathway crosstalk.

Animals↗
PubMed2024

Organ-delimited gene regulatory networks provide high accuracy in candidate transcription factor selection across diverse processes.

Organ-specific gene expression datasets that include hundreds to thousands of experiments allow the reconstruction of organ-level gene regulatory networks (GRNs). However, creating such datasets is greatly hampered by the requirements of extensive and tedious manual curation. Here, we trained a supervised classification model that can accurately classify the organ-of-origin for a plant transcriptome. This K-Nearest Neighbor-based multiclass classifier was used to create organ-specific gene expression datasets for the leaf, root, shoot, flower, and seed in Arabidopsis thaliana. A GRN inference approach was used to determine the: i. influential transcription factors (TFs) in each organ and, ii. most influential TFs for specific biological processes in that organ. These genome-wide, organ-delimited GRNs (OD-GRNs), recalled many known regulators of organ development and processes operating in those organs. Importantly, many previously unknown TF regulators were uncovered as potential regulators of these processes. As a proof-of-concept, we focused on experimentally validating the predicted TF regulators of lipid biosynthesis in seeds, an important food and biofuel trait. Of the top 20 predicted TFs, eight are known regulators of seed oil content, e.g., WRI1, LEC1, FUS3. Importantly, we validated our prediction of MybS2, TGA4, SPL12, AGL18, and DiV2 as regulators of seed lipid biosynthesis. We elucidated the molecular mechanism of MybS2 and show that it induces purple acid phosphatase family genes and lipid synthesis genes to enhance seed lipid content. This general approach has the potential to be extended to any species with sufficiently large gene expression datasets to find unique regulators of any trait-of-interest.

Arabidopsis↗
PubMed2026

Disagreement-informed arbitration for gene regulatory network inference: A score-level meta-classifier and a diagnostic typology of inter-method conflict.

Gene regulatory network inference methods routinely disagree about individual edges, and practitioners resolve those conflicts by choosing one method or averaging them all. We ask whether the conflict can instead be arbitrated per edge. A gradient-boosted classifier is trained on the raw scores that ten inference methods-correlation-based, information-theoretic, sparse-regression and tree-ensemble, including GENIE3, GRNBoost2, CLR and ARACNe-assign to each candidate regulator-target pair, so that the weight given to each method varies from edge to edge. Across six single-cell perturbation screens spanning four cell types, arbitration improves on mean ensembling by +0.056 AUROC on Adamson and +0.083 on Shifrut under target-grouped cross-validation. The evaluation protocol turns out to matter more than the model. Edge-level cross-validation, standard in this literature, inflates apparent gains by 0.060 AUROC through target-gene leakage-comparable to the entire honest improvement. The effect is far larger for methods that represent genes implicitly: a supervised graph-attention link predictor trained on identical folds scores AUROC 0.930 under edge-level cross-validation, better than anything else we evaluate, and 0.533 once target genes are held out. Any method that parameterises genes is exposed, which covers most graph- and embedding-based approaches. A five-category typology of inter-method conflict localises where arbitration pays off, with the largest gains on edges where the methods disagree and the smallest where they already agree, while adding nothing as model input; we therefore report it as a diagnostic instrument rather than a modelling contribution. We also characterise what the ground truth measures: most perturbed genes in widely used screens are not transcription factors, and a mediation screen bounds how much of the perturbation response can be direct.

Ensemble methods↗
PubMed2026

Comparative analysis of conserved non-coding elements identifies gene regulatory networks rewired during the water-to-land transition in vertebrates.

The conquest of land by vertebrates has been a pivotal moment in evolutionary history. Adapting to the new habitats necessitated numerous changes in vertebrate anatomy and physiology, creating an enduring imprint on the developmental gene regulatory networks (GRNs) of tetrapods. The increase of high-quality genomic resources over the past decade has made it possible to study the genomic legacy of the water-to-land transition. While much attention has been given to the highly conserved non-coding elements (CNEs) of the genome that share high levels of similarity across evolutionarily diverged clades, recent evidence suggests that perhaps comparable attention should be given to "missing" CNE-s, conserved sequence patches present in extant stem gnathostomes and actinopterygian fishes that have become undetectable in tetrapods during the adaptation to terrestrial life, whether through true sequence loss or divergence beyond alignability. These sequences could help us reveal the relaxation of certain developmental constraints, related to the aquatic lifestyle, that made reaching new adaptive peaks in the developmental landscape possible. In this paper, we search for such CNEs and characterize them in comparison with pan-Gnathostome CNEs, using the zebrafish (Danio rerio) genome as a reference. Our results suggest that the rewiring of developmental networks related to pigmentation and muscle structure formation has left the largest genomic imprint. We also find that components of canonical Wnt and Hedgehog signalling, are enriched among CNEs retained in fish.

cis-regulatory evolution↗
PubMed2025

Integrated multi-omics analyses provide new insights into genomic variation landscape and regulatory network candidate genes associated with walnut endocarp.

Persian walnut (Juglans regia) is an economically important nut oil tree; the fruit has a hard endocarp/shell to protect seeds, thus playing a key role in its evolution, and the shell thickness is an important trait for walnut breeding. However, the genomic landscape and the gene regulatory networks associated with walnut shell development remain to be systematically elucidated. Here, we report a high-quality genome assembly of the walnut cultivar 'Xiangling' and construct a graphic structure pan-genome of eight Juglans species to reveal the genetic variations at the genome level. We re-sequence 285 accessions to characterize the genomic variation landscape. Through genome-wide association studies (GWAS), we identified 19 loci associated with more than 268 loci that underwent selection during walnut domestication and improvement. Multi-omics analyses, including transcriptomics, metabolomics, DNA methylation, and spatial transcriptomics across eleven developmental stages, revealed several candidate genes related to secondary cell biosynthesis and lignin accumulation. This integrated multi-omics approach revealed several candidate genes associated with secondary cell biosynthesis and lignin accumulation, such as UGP, MYB308, MYB83, NAC043, NAC073, CCoAOMT1, CCoAOMT7, CHS2, CESA7, LAC7, COBL4, and IRX12. Overexpression of JrUGP and JrMYB308 in Arabidopsis thaliana confirmed their roles in lignin biosynthesis and cell wall thickening. Consequently, our comprehensive multi-omics findings offer novel insights into walnut genetic variation and network regulation of endocarp development and shell thickness, which enable further genome-informed breeding strategies for walnut cultivar improvement.

Juglans↗