PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Gene regulatory networks”

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 19 recordsLinked to original sources
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↗
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

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

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

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↗
PubMed2026

Plasma proteins are integral to cross-tissue gene regulatory networks implicated in cardiometabolic disorders and coronary artery disease.

The plasma proteome has demonstrated promise for identifying diagnostic markers for cardiometabolic disorders (CMDs) and coronary artery disease (CAD). However, identifying the organ of origin for these biomarkers is critical for establishing biological relevance. We performed a multi-omic integrative analysis across multiple tissues from the STARNET study by profiling 974 plasma proteins in 532 CAD patients, integrating RNA sequencing (RNA-seq) data from the arterial wall, major metabolic organs, and blood. We identified 144 cis-protein quantitative trait loci in plasma, colocalizing with tissue cis-expression quantitative trait loci. Additionally, by mapping tissue mRNA "seed genes," we traced 262 plasma proteins to their source organs, primarily the liver. Crucially, we found that 851 plasma proteins are associated with the activity of cross-tissue gene regulatory networks (GRNs), including GRNs implicated in CMD and CAD development. Our findings demonstrate that plasma proteins are integral components of GRNs, with potential for developing reliable diagnostics and precise therapeutic targets. A record of this paper's transparent peer review process is included in the supplemental information.

cardiometabolic disorders↗
PubMed2025

Characterization of FLOWERING LOCUS T-related genes and their putative gene regulatory network in semi-winter Brassica napus cultivar Zhongshaung11.

In many species, FLOWERING LOCUS T (FT)-like genes promote the floral transition by integrating environmental signals, in particular photoperiod, and internal cues. Here we show that Brassica napus contains six FT-like genes and two pseudogenes belonging to three orthogroups. All B. napus FT-like genes induce early flowering when expressed at the shoot apical meristems of Arabidopsis thaliana ft mutants; however, BnaFT.C6 and non-orthologous FT-like genes do not encode fully functional mobile florigens. In the case of BnFT.C6, the functional change is associated with a T to C amino acid change that is restricted to semi-winter accessions. Expression of orthologs of FT is photoperiod-dependent, and two distal enhancers are conserved; however, the homeologs BnaFT.A7 and BnaFT.C6 show rearrangements of DNA motifs binding NF-Y/CO and NF-Y transcriptional activator complexes between the promoter and downstream enhancers. Motif rearrangements correlate with differences in tissue-specific expression. Furthermore, homeologs with rearranged motifs could not be transactivated by B. napus CO in transient assays, although they show LD photoperiod-dependent expression. We propose that differential diurnal expression of NF-Y genes contributes to the photoperiod-dependent regulation of B. napus FT genes.

Brassica napus↗
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

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↗
PubMed2024

Multi-organ gene expression analysis and network modeling reveal regulatory control cascades during the development of hypertension in female spontaneously hypertensive rat.

Hypertension is a multifactorial disease with stage-specific gene expression changes occurring in multiple organs over time. The temporal sequence and the extent of gene regulatory network changes occurring across organs during the development of hypertension remain unresolved. In this study, female spontaneously hypertensive (SHR) and normotensive Wistar Kyoto (WKY) rats were used to analyze expression patterns of 96 genes spanning inflammatory, metabolic, sympathetic, fibrotic, and renin-angiotensin (RAS) pathways in five organs, at five time points from the onset to established hypertension. We analyzed this multi-dimensional dataset containing ~15,000 data points and developed a data-driven dynamic network model that accounts for gene regulatory influences within and across visceral organs and multiple brainstem autonomic control regions. We integrated the data from female SHR and WKY with published multiorgan gene expression data from male SHR and WKY. In female SHR, catecholaminergic processes in the adrenal gland showed the earliest gene expression changes prior to inflammation-related gene expression changes in the kidney and liver. Hypertension pathogenesis in male SHR instead manifested early as catecholaminergic gene expression changes in brainstem and kidney, followed by an upregulation of inflammation-related genes in liver. RAS-related gene expression from the kidney-liver-lung axis was downregulated and intra-adrenal RAS was upregulated in female SHR, whereas the opposite pattern of gene regulation was observed in male SHR. We identified disease-specific and sex-specific differences in regulatory interactions within and across organs. The inferred multi-organ network model suggests a diminished influence of central autonomic neural circuits over multi-organ gene expression changes in female SHR. Our results point to the gene regulatory influence of the adrenal gland on spleen in female SHR, as compared to brainstem influence on kidney in male SHR. Our integrated molecular profiling and network modeling identified a stage-specific, sex-dependent, multi-organ cascade of gene regulation during the development of hypertension.

Animals↗
PubMed2025

A mouse organoid platform for modeling cerebral cortex development and cis-regulatory evolution in vitro.

Natural selection has shaped the gene regulatory networks that orchestrate cortical development, leading to structural and functional variation across mammals, but the molecular and cellular mechanisms underpinning these changes have only begun to be characterized. Here, we develop a reproducible protocol for cerebral cortex organoid generation from mouse epiblast stem cells (EpiSCs), which recapitulates the timing and cellular differentiation programs of the embryonic cortex. We generated cortical organoids from F1 hybrid EpiSCs derived from crosses between laboratory mice (C57BL/6J) and four wild-derived inbred strains spanning ∼1 M years of evolutionary divergence to comprehensively map cis-acting transcriptional regulatory variation across developing cortical cell types, using single-cell RNA sequencing (scRNA-seq). We identify hundreds of genes that exhibit dynamic allelic imbalances, providing the first insight into the developmental mechanisms underpinning changes in cortical structure and function between subspecies. These experimental methods and cellular resources represent a powerful platform for investigating gene regulation in the developing cerebral cortex.

Organoids↗
PubMed2026

Insights into glandular trichome biology from analysis of organ-specific gene expression programmes in cannabis, hop and tomato.

Glandular trichomes (GTs) are epidermal outgrowths in which diverse specialised (secondary) metabolites are synthesised and stored. Cannabis (Cannabis sativa L.) and its close relative hop (Humulus lupulus L.) have pharmaceutical and industrial significance due to the presence of these metabolites in their GTs. We examined the conservation or divergence of the specific transcriptional programmes underlying GT biology. To achieve this, we generated transcriptome atlases of trichomes, flower, leaf, stem and root for cannabis, hop and tomato. We found that 12.9, 10.1 and 16.8% of cannabis, hop and tomato genes, respectively, were expressed organ/tissue specifically across all organs/tissues. Transcription factors (TFs) on average accounted for 7.5% of the organ-specific transcriptome and likely regulate organ-specific functions. We also conducted weighted gene co-expression network analysis and gene regulatory network (GRN) analysis to identify key regulators of GT function across the species and validated our predictions by DNA affinity purification sequencing for a subset of the cannabis and tomato GT TFs. The GRNs specific to cannabis or hop GTs were enriched for TFs and target genes associated with specialised metabolism, reflecting their species-specific nature. Conversely, the shared GRN components (identified via orthology analysis) were involved in highly conserved processes, such as flavonoid biosynthesis, solute transport and metabolite storage. Together, these GRNs and the associated transcriptome atlases are valuable resources to improve our knowledge of GT function and organ-specific genome regulation.

Solanum lycopersicum↗