SEARCH · Search PubMed
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.
Explore the source record for details and available documents.
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.
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.
Although microRNA-22-3p (miR-22-3p) is abundantly expressed in the avian liver, its epigenetic role in lipid homeostasis remains largely uncharacterized. To elucidate its in vivo function, 14-day-old female Qingyuan Partridge chickens were intravenously injected with lentiviral vectors to establish miR-22-3p overexpression and knockdown models. Phenotypic analysis demonstrated that miR-22-3p knockdown significantly elevated hepatic triglyceride (TG) levels (p < 0.05) and drove marked steatosis, whereas its overexpression reduced TG content. Transcriptome sequencing (RNA-Seq) revealed profound metabolic remodeling, identifying 23 core lipid-associated genes (e.g., ELOVL6, FADS2, ACSBG2, and PTGIS) heavily enriched in steroid biosynthesis, fatty acid metabolism, and elongation pathways. In conclusion, miR-22-3p functions as a bidirectional epigenetic rheostat that negatively regulates hepatic lipid deposition by orchestrating a multilayered polygenic network, providing novel molecular targets for mitigating avian metabolic disorders and optimizing production traits in indigenous poultry breeds.
Golden pompano (Trachinotus ovatus) is one of the most economically important marine fish species in China. It is susceptible to low-temperature stress, which significantly challenges its production and supply. Nevertheless, study on the regulatory mechanisms underlying low-temperature stress responses in golden pompano remains limited. Here, we firstly performed a time-series transcriptome analysis to reconstruct dynamic response patterns under low-temperature stress in golden pompano. Transcriptome profiling identified common differentially expressed genes (DEGs), including fos, hlf, and hmgb1, as well as condition-specific DEGs across distinct low-temperature stress groups. Based on cluster analysis, all DEGs were classified into five distinct expression patterns, reflecting diversified regulation of expression in golden pompano during low-temperature stress. Furthermore, condition-specific regulatory modules were explored via weighted gene co-expression network analysis (WGCNA), highlighting that the two module hub genes, serbf2 and lipc, might respond to low-temperature stress by regulating the lipid catabolic process. Subsequently, untargeted metabolomic analysis revealed that glycerophospholipid metabolism was a significantly enriched common pathway, highlighting its crucial role in mediating the response to low-temperature stress. Finally, by integrating transcriptomic and metabolomic analyses, a gene-metabolite interaction network associated with glycerophospholipid metabolism under low-temperature stress was established. These findings underscore the significance of multiple candidate genes and glycerophospholipid metabolism in golden pompano's response to low-temperature stress, thereby laying a solid molecular foundation for the development of low-temperature-tolerant fish strains.
Non-coding RNAs (ncRNAs) have been demonstrated to be involved in muscle development and to function as key regulators. However, the molecular mechanism underlying muscle crispiness in grass carp (GC) remains poorly understood, and whether these ncRNAs are involved in its regulation is still unknown. In the current investigation, differentially expressed (DE) RNAs (including lncRNAs, circRNAs, miRNAs, and mRNAs) were identified; concomitantly, target genes prediction was conducted, and functional and signaling pathway enrichment analyses were performed. Pathways related to muscle crispiness were identified, and the competitive endogenous RNA (ceRNA) (lncRNA/circRNA-miRNA-mRNA) regulatory network was further constructed. The results showed that a total of 126 DE-lncRNAs, 17 DE-circRNAs, 329 DE-miRNAs, and 442 DE-mRNAs were identified in muscle tissues of both the GC and crisp grass carp (CGC). GO and KEGG enrichment analyses revealed that target genes of DE-ncRNAs were significantly enriched in signaling pathways, including structural constituents of muscle, apoptosis, oxidative phosphorylation, and regulation of actin cytoskeleton, suggesting that these pathways may be involved in muscle texture remodeling. Subsequently, DE-RNAs enriched in related pathways were identified, and a core ceRNA regulation network comprising 3 lncRNAs, 4 circRNAs, 3 miRNAs, and 17 mRNAs was constructed. Additionally, 10 DE-RNAs from randomly selected groups were validated by qRT-PCR. Our findings not only provide scientific evidence elucidating the molecular mechanisms underlying muscle crispiness in GC but also establish a foundation for studying changes in muscle textural qualities across other fish species.
Soil salinization severely threatens global food security, necessitating systematic investigations of halophytes like Portulaca oleracea to decode the molecular mechanisms of environmental resilience. Utilizing an integrated framework of deep learning-based genome annotation (58,817 predicted genes; 96.5% BUSCO completeness), multi-tissue RNA-Seq, phylogenomics, and gene regulatory network (GRN) inference, the synergistic orchestration of 78 aquaporins (AQPs), 525 heat shock proteins (HSPs), and 119 late embryogenesis abundant (LEA) proteins was elucidated. The active transcriptome, encompassing 39,065 expressed loci, revealed a systemic growth-defense trade-off. Tissues displayed distinct adaptive mechanisms: leaves modulated intracellular water balance via specialized AQPs, whereas adult roots maintained proteostasis through robust HSP20/HSP70 induction. Phylogenomic clustering across 154 species demonstrated that salinity tolerance constitutes an evolutionary mosaic, identifying 81 halophyte-exclusive orthogroups and 1129 species-specific clusters. Comparative topology across six independent GRNs (4.2M-5.3 M edges) unmasked a highly modular transcriptional reprogramming strategy governed by a core apparatus of 22 stress-exclusive regulators, with functional enrichment heavily prioritizing protein dimerization and chromatin remodeling. Theoretically, the distinct convergence of Trihelix transcription factors with guard cell differentiation pathways offers a candidate transcriptomic framework to explain the plant's characteristic C4-CAM photosynthetic plasticity under severe osmotic pressure. Practically, these evolutionary blueprints and specific master switches transcend single-gene transgenic limitations. Utilizing these root-sustained and stress-inducible targets under localized promoters provides a naturally optimized, network-level precision engineering roadmap to transfer robust, compartmentalized halotolerance to sensitive glycophytic crops.
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.
Long non-coding RNAs (lncRNAs) stand as newly-arisen molecular types that exert regulatory effects, able to operate as competitive endogenous RNAs (ceRNAs) to engage microRNAs (miRNAs) in interaction, resulting in the recovery of target mRNA expression and activity. Increasing evidences indicate that the ceRNA network affects various biological processes in mammals, including development, cellular differentiation, metabolism, immune response, and disease pathogenesis. In teleost fish, the lncRNA-miRNA-mRNA regulatory networks have been reported occasionally. However, up to now, the roles of lncRNAs in the big-belly seahorse (Hippocampus abdominalis) remains unclear. In this study, we reported for the first time, via whole-transcriptome RNA sequencing, the lncRNA mediated ceRNA regulatory network in Vibrio harveyi-infected H. abdominalis. A total of 4197 differentially expressed mRNAs (DE-mRNAs), 1317 DE-lncRNAs, and 183 DE-miRNAs were identified. Furthermore, the crosstalk between miRNAs and lncRNAs as well as between miRNAs and mRNAs was inferred based on the negative correlations between miRNAs and their target lncRNAs/mRNAs. A core immune associated lncRNA-miRNA-mRNA putative regulatory network was thus constructed, comprising 211 lncRNA-miRNA and 224 mRNA-miRNA pairs. In conclusion, our findings provide an integrative overview of the ceRNA regulatory networks on the underlying immune responses to V. harveyi infection in the big-belly seahorse, and offer a solid theoretical foundation for the comparative immunological research of teleost fish.
Immune function across development, tissue repair, aging, and disease depends not only on signaling pathways but also on epigenetic architectures that determine whether coordinated transcriptional programs can be accessed and resolved. Increasing evidence indicates that epigenetic gene networks regulate the accessibility and reversibility of semi-stable immune states, shaping plastic, homeostatic, reparative, and degenerative configurations. We propose the concept of epigenetic transition windows, defined as temporally and contextually restricted intervals during which epigenetic constraints are relaxed, permitting coordinated and reversible transitions between immune states. During development, these windows are broad and support immune tolerance and adaptive plasticity. In adulthood they become spatially and temporally restricted, preserving stability while enabling conditional adaptation. With aging, they progressively narrow, contributing to chronic inflammation, impaired repair, and increased vulnerability to neurodegeneration. Conversely, pathological persistence of regulatory permissiveness may underlie immune evasion and sustained plasticity in cancer. We outline operational genomic readouts for quantifying transition windows, including chromatin accessibility variance, enhancer switching dynamics, reversibility metrics, and cross-cell coordination indices, and derive experimentally testable predictions that distinguish this model from pathway-centric or damage-centric explanations. By reframing immune dysfunction as a failure of regulated state transition rather than excessive signaling alone, this framework integrates inflammaging, trained immunity, immune resolution failure, and tumor immune escape within a unified regulatory architecture and provides a systems-level perspective on immune adaptability across the lifespan.
Glaucomas are heterogeneous optic neuropathies associated with extracellular matrix dysregulation, abnormal ocular morphogenesis, and inflammatory signaling. Targeted deep sequencing of 586 primary congenital glaucoma (PCG) cases and 1,757 controls identified rare pathogenic variants in multiple genes, including WT1 and MMP9. Notably, WT1 variants clustered within the nuclear export sequence. Further, functional analyses showed that combined wt1-pax6 suppression in zebrafish disrupted ocular morphogenesis, highlighting developmental interdependence. In human trabecular meshwork cells, WT1 acted as a transcriptional repressor of MMP9, while TNF-α signaling triggered nitric oxide-dependent nuclear export of WT1, resulting in delayed MMP9 upregulation. This effect was reversible by inhibiting nuclear export or nitric oxide synthase. A patient-derived mutation in the nuclear-export region of WT1, disrupted this regulatory switch, causing abnormal MMP9 expression. These findings position WT1 as an important regulator linking developmental and inflammatory mechanisms in glaucoma pathogenesis.
High-fat diets are increasingly used in aquaculture due to their protein-sparing effects; however, the post-transcriptional regulatory mechanisms of fish muscle in response to high-fat diets (HFD) remain unclear. In this study, juvenile black seabream were fed either a normal-fat diet (NFD) or a HFD to investigate the miRNA-mRNA regulatory network associated with diet-induced muscle lipid deposition. Oil Red O staining and biochemical analysis showed that high-fat diet feeding markedly increased lipid droplet accumulation and crude lipid content in muscle, indicating significant induction of muscle lipid deposition. Integrated mRNA and miRNA expression profiling revealed substantial transcriptomic and post-transcriptional responses to high-fat diet challenge. A total of 271 differentially expressed genes were identified, including 120 upregulated and 151 downregulated genes. Through combined target prediction and expression correlation analysis, thirteen candidate inverse miRNA-mRNA relationships were subsequently identified, and RT-qPCR supported the expression patterns of selected miRNAs and mRNAs. These pairs included miR-499-x-dmgdh, miR-499-y-gatm, miR-727-y-ass1, miR-4649-x-foxo4, miR-9129-z-myl7, and several novel miRNA-mediated interactions involving adk, chst11, lypla2, frem2, kcnc4, wars1, bag2, and capn2. Functional analysis suggested that these regulatory pairs were mainly associated with metabolic adaptation, structural remodeling, and cellular stress responses. In particular, gatm, dmgdh, ass1, and adk were associated with energy metabolism-related processes, including pathways previously linked to Ampk regulation, whereas myl7, frem2, and kcnc4 may contribute to muscle structural maintenance and excitability regulation. Overall, this study provides candidate miRNA-mRNA regulatory relationships potentially involved in high-fat diet-induced muscle lipid deposition and adaptive remodeling in black seabream, offering a basis for future functional studies on muscle metabolism and quality regulation in marine fish.
In hybrid plants, phenotypic outcomes are governed by interactions between the two parental genomes. However, the mechanisms underlying the interplay of divergent regulatory networks from these genomes remain poorly understood. In this study, we compared gene-level and allele-specific expression patterns, as well as differentially enriched pathways between F₁ and complex backcross (CBC) lines derived from a natural interspecific hybrid population of Populus fremontii (Pf) and P. angustifolia (Pa). Metabolic differences between Pf and Pa which exhibit low and high levels respectively of phenylpropanoid-derived condensed tannins were leveraged. Using individualized transcriptome references, differential expression and clustering analyses revealed CBC-biased and F₁-biased expression for genes involved in phenylpropanoid metabolism and photosynthesis, respectively. Biased expression of these genes at the allele level was also observed in F1. At the whole-transcriptome level, Pa-biased genes predominated in F₁ hybrids, and Pa alleles displayed more conserved expression patterns than Pf alleles across examined samples. Further analyses indicated that allelic expression bias was significantly associated with parental origin, which could be driven by sequence variations in cis-regulatory elements and differences in CpG island length. Our findings demonstrate strong parent-of-origin effects on divergent regulatory networks governing gene expression in poplar hybrids and provide clues for strategic parental selection tailored to specific metabolic pathways of interest.
BACKGROUND: Endocrine-disrupting chemicals (EDCs) are linked to coronary heart disease (CHD), but underlying mechanisms remain unclear. We aimed to identify EDC-related genes and evaluate their causal roles in CHD. METHODS: We curated EDC-related genes from a compound-gene interaction database and integrated them with CHD genome-wide association study (GWAS) summary statistics and tissue-specific expression quantitative trait loci (eQTL) data. Two-sample Mendelian randomization (MR) and Bayesian colocalization were applied to infer causality. Functional enrichment, single-cell RNA sequencing of human coronary arteries, and EDC-gene networks were further analyzed. RESULTS: After FDR correction, 39 genes were significantly associated with CHD risk via MR. Four genes-ZNF827, FCHO1, IPO9 (protective), and RPL13 (risk-increasing)-showed strong colocalization (PPH4 > 0.9). Pathway and single-cell analyses of coronary artery tissue indicated that vascular and immune pathways mediate these effects. An interaction network highlighted associations between specific EDCs and candidate genes implicated in CHD susceptibility. CONCLUSION: This integrative genomic study provides evidence that EDCs influence CHD susceptibility through distinct gene networks, revealing potential mechanisms and molecular targets for prevention and therapy.
The irreversibility of complex trait loss has long been a tenet of evolutionary biology. However, this idea is increasingly at odds with the numerous documented exceptions across the Tree of Life. We synthesise this growing body of evidence across a diverse array of taxa and traits, exploring the evolutionary conditions that enable evolutionary reversal. By integrating macroevolutionary, genetic, and developmental information, we argue that trait reversal is commonly fostered by some form of persistence in the generative developmental pathway of the lost trait. We identify three overarching modes of trait reversal and support them with multiple case studies: by pleiotropy (the involvement of the same generative components in other traits and/or functions), by plasticity (environment-dependent expression of the trait) and by hemiplasy (persistence in another lineage, followed by reticulate evolution). We also examine important affinities between trait reversal and evolutionary novelties, undermining a neat distinction between what is old and what is new in evolution. This survey may provide a useful framework for future explorations of the developmental mechanisms underlying these still overlooked macroevolutionary dynamics.
BACKGROUND: This study aimed to profile the molecular signatures of post-traumatic elbow heterotopic ossification (HO) to identify key regulators and potential therapeutic targets. METHODS: Total RNA from post-traumatic elbow HO tissues (n=4) and normal bone tissues (n=6) was subjected to high-throughput sequencing to identify differentially expressed mRNAs (DEGs), microRNAs (DEMs), and lncRNAs (DELs). Bioinformatics analyses included Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, protein-protein interaction network construction, and transcription factor (TF)-microRNA-mRNA network analysis. The expression trends of four most upregulated and four most downregulated DEGs were validated by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR). RESULTS: We identified 2,138 DEGs, 40 DEMs, and 905 DELs. DEGs were significantly enriched in biological process "bone mineralization," cellular component "plasma membrane," molecular function "integrin binding," and pathways including PI3K-Akt, NF-κB, JAK-STAT, and TNF signaling pathways. Hub genes with high connectivity included MMP9, IL6, MMP3, CTSK, and BGLAP. Integrated network analysis highlighted the transcription factor JUN and key microRNAs (hsa-miR-124-3p, hsa-miR-548c-3p, and hsa-miR-135b). The qRT-PCR results confirmed the expression trends of selected DEGs. CONCLUSIONS: This study, for the first time, profiled the differentially expressed mRNAs, microRNAs, and lncRNAs in post-traumatic elbow HO using high-throughput RNA sequencing. These findings provide valuable insights into the molecular mechanisms of HO following elbow trauma. The identified hub genes (MMP9, IL6, MMP3, CTSK, and BGLAP), key TF (JUN), and key microRNAs (hsa-miR-124-3p, hsa-miR-548c-3p, and hsa-miR-135b) may serve as potential therapeutic targets for preventing and treating post-traumatic elbow HO.
Papillary thyroid cancer (PTC) is the most prevalent endocrine malignancy, accounting for over 90% of thyroid cancers. While differentiated thyroid cancers (DTCs) typically have favorable outcomes, a significant subset progresses to radioactive iodine-refractory (RAIR) disease, characterized by impaired iodine uptake and a 10-year survival rate below 10%. Genetic alterations and dysregulated signaling pathways underlie this transition. Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs), play critical regulatory roles in tumor biology and may be transported via exosomes, facilitating intercellular communication and contributing to RAIR-PTC. This systematic review, conducted according to PRISMA 2020 guidelines, evaluated the role of exosomal ncRNAs in RAIR-PTC. A comprehensive search of PubMed, PubMed Central, and Google Scholar identified studies published within the past 15 years in English. Following stringent quality appraisal, studies with a non-bias score above 40% were included. Of 961 identified publications, 96 high-quality studies met inclusion criteria. Evidence indicates that therapy resistance in RAIR-PTC is driven by convergent ncRNA regulatory networks that suppress sodium-iodide symporter (NIS) expression and activate oncogenic pathways, most notably MAPK, PI3K/AKT/mTOR, and Wnt/β-catenin signaling. Multiple ncRNAs converge on key regulatory nodes, forming redundant circuits that sustain dedifferentiation, metabolic adaptation, and impaired iodide transport. Several consistently dysregulated ncRNAs directly or indirectly regulate NIS expression and trafficking, highlighting actionable targets. Exosomes emerge as biologically compatible, programmable delivery vehicles capable of transporting therapeutic ncRNA payloads independent of endogenous packaging mechanisms. These findings support a precision therapeutic paradigm in which engineered exosomes reprogram ncRNA networks to restore iodine-handling pathways and overcome therapy resistance in RAIR-PTC.
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.