SEARCH · Search PubMed

Results for “Integrative omics”

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

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29Linked to original sources

PubMed2006

Evaluation of a non-targeted "omic" approach in the safety assessment of genetically modified plants.

Genetically modified plants must be approved before release in the European Union, and the approval is generally based upon a comparison of various characteristics between the transgenic plant and a conventional counterpart. As a case study, focusing on safety assessment of genetically modified plants, we here report the development and characterisation of six independently transformed ARABIDOPSIS THALIANA lines modified in the flavonoid biosynthesis. Analyses of integration events and comparative analysis for characterisation of the intended effects were performed by PCR, quantitative Real-time PCR, and High Performance Liquid Chromatography. Analysis by cDNA microarray was used as a non-targeted approach for the identification of potential unintended effects caused by the transformation. The results revealed that, although the transgenic lines possessed different types of integration events, no unintended effects were identified. However, we found that the majority of genes showing differential expression were identified as stress-related genes and that environmental conditions had a large impact on the expression of several genes, proteins, and metabolites. We suggest that the microarray approach has the potential to become a useful tool for screening of unintended effects, but state that it is crucial to have substantial information on the natural variation in traditional crops in order to be able to interpret "omics" data correctly within the framework of food safety assessment strategies of novel plant varieties, including genetically modified plant varieties.

Arabidopsis↗

PubMed2005

Towards quantitative biology: integration of biological information to elucidate disease pathways and to guide drug discovery.

Developing a new drug is a tedious and expensive undertaking. The recently developed high-throughput experimental technologies, summarised by the terms genomics, transcriptomics, proteomics and metabolomics provide for the first time ever the means to comprehensively monitor the molecular level of disease processes. The "-omics" technologies facilitate the systematic characterisation of a drug target's physiology, thereby helping to reduce the typically high attrition rates in discovery projects, and improving the overall efficiency of pharmaceutical research processes. Currently, the bottleneck for taking full advantage of the new experimental technologies are the rapidly growing volumes of automatically produced biological data. A lack of scalable database systems and computational tools for target discovery has been recognised as a major hurdle. In this review, an overview will be given on recent progress in computational biology that has an impact on drug discovery applications. The focus will be on novel in silico methods to reconstruct regulatory networks, signalling cascades, and metabolic pathways, with an emphasis on comparative genomics and microarray-based approaches. Promising methods, such as the mathematical simulation of pathway dynamics are discussed in the context of applications in discovery projects. The review concludes by exemplifying concrete data-driven studies in pharmaceutical research that demonstrate the value of integrated computational systems for drug target identification and validation, screening assay development, as well as drug candidate efficacy and toxicity evaluations.

Animals↗

PubMed2026

Genome-wide association and multi-omics functional screens reveal the genetic architecture of foveal development.

Foveal hypoplasia causes visual impairment across congenital eye disorders, yet the genetic programmes governing foveal development remain poorly characterised and no tractable model exists for foveal disease. In the first genome-wide association study of foveal hypoplasia, we identified 42 sentinel variants mapping to 54 effector genes supported by ≥ 2 criteria from a variant-to-gene framework incorporating developmental multi-omics. Disruption of six effector genes using mutant lines and CRISPR knockouts in the zebrafish high acuity zone recapitulates structural, functional, and ultrastructural hallmarks of foveal hypoplasia, establishing the first vertebrate disease model. Integration with human foetal single-cell and spatial transcriptomics reveals two temporal waves of effector gene expression and identifies Müller glia as critical mediators of foveal patterning. Phenome-wide analyses reveal foveal variants are pleiotropic with refractive, lenticular, and metabolic traits, connecting foveal development to anterior segment and systemic disease biology. These findings should inform mechanistic studies of macular disease.

Journal Article↗

PubMed2025

Multi-Omics Analysis of the Potential Mechanisms of Skin Albinism in Edangered Percocypris pingi: Abnormal Ubiquitination and Calcium Signal Inhibition.

Percocypris pingi is an endangered protected fish species in China. Its albino variants exhibit growth retardation and physiological abnormalities. Understanding its albinism mechanism holds significant scientific importance for molecular breeding programs and disease model development. This study integrated transcriptomic and proteomic analyses, combined with histopathological and molecular biological techniques, to systematically compare molecular differences in skin tissues between albino and wild-type P. pingi, with a focus on elucidating the multidimensional regulatory mechanisms underlying skin albinism. Our findings suggest that albinism in P. pingi is synergistically driven by hyperactivation of ubiquitin-mediated proteolysis (which suppressed TYR/TYRP1 enzymatic activity and disrupted the pH homeostasis of melanosomes), and inhibition of calcium signaling (which impeded melanin transport). This discovery provides novel insights into the mechanisms of pigment loss in fish species and offers a valuable reference for molecular breeding of endangered species as well as research on pigmentation-related disorders.

Animals↗

PubMed2027

Systematic Analysis of Tumor Microenvironment Using IOBR.

The Immuno-Oncology Biological Research (IOBR) package is an R-based analysis tool for exploring the tumor microenvironment (TME) and its influence on anti-tumor immunity. Built for high-throughput data-spanning both transcriptomic and genomic profiles-IOBR integrates six analytical modules, including transcriptomic data preprocessing, TME profiling, TME pattern identification, ligand-receptor interaction analysis, genome-TME interaction assessment, and visualization. In this chapter, we walk through a multi-omics workflow using example datasets, illustrating data preparation, distribution analyses, result interpretation, and graphical output. IOBR is open source and is available at https://github.com/IOBR/IOBR and a detailed GitBook ( https://iobr.github.io/book/ ) offers a complete manual and analysis guide for each function.

Tumor Microenvironment↗

PubMed2026

NMR metabolomics and glycomics for cancer detection in patients with non-specific symptoms: a prospective observational cohort study.

BACKGROUND: Early cancer diagnosis in patients with non-specific symptoms is limited by the lack of discriminatory tests. Within the Oxfordshire Suspected CANcer (SCAN) pathway, exploratory biomarker work showed that serum 1H NMR-based metabolomics can identify cancer with high accuracy. SCAN2 evaluated whether integrating metabolomics with glycomics provides complementary molecular information and improves discrimination in a clinically complex, real-world population. METHODS: Serum from 369 SCAN patients (59 cancers) was analysed using AXINON® System-derived NMR metabolomics and HPLC-MS glycomics. Machine-learning models were trained to predict cancer status, with performance assessed by receiver operating characteristic (ROC) analysis of pooled cross-validated predictions. To place cancer risk in a broader clinical context, a second classifier modelling alternative non-cancer diagnosis was incorporated, and mean predicted probabilities from both models were jointly projected into a two-dimensional space, maintaining strict separation of training and test data. FINDINGS: In the full cohort, integration of glycomics with metabolomics achieved an AUC of 0.814 (95% CI 0.808-0.820). In a refined sub-cohort excluding major comorbidities and selected cancer types (32 cancers, 277 non-cancers), performance improved to an AUC of 0.884 (95% CI 0.879-0.890). Discriminatory features included cancer-associated biantennary fucosylated glycans alongside amino acid metabolites (glutamate, histidine) and lipoprotein-related measures. A classifier distinguishing metastatic from non-metastatic disease (n = 29 vs. 30) achieved an AUC of 0.80. Joint probability analysis in the full cohort preserved cancer-associated signatures across comorbidity burden, with projection-based classification achieving an accuracy of 89.2% (95% CI 85.7-92.6). INTERPRETATION: These findings validate the SCAN1 metabolomic signature in a more clinically complex cohort and indicate that integrating glycomics with metabolomics provides complementary biological information for cancer discrimination. Joint probability analysis provides an interpretable framework for cancer risk stratification within multimorbid diagnostic pathways, supporting the clinical potential of scalable multi-omics blood testing. FUNDING: EPSRC, EU Horizon 2020, Wellcome/MLSTF, Novo Nordisk Foundation.

Humans↗

PubMed2005

VIS-O-BAC: exploratory visualization of functional genome studies from bacteria.

UNLABELLED: The visualization-aided exploration of complex datasets will allow the research community to formulate novel functional hypotheses leading to a better understanding of biological processes at all levels. Therefore, we have developed a web resource termed VIS-O-BAC designed for the functional investigation of expression data for model systems, such as bacterial pathogens based on a graphical display. Genome-scale datasets derived from typical 'omic' approaches can directly be explored with respect to three biologically relevant aspects, the genome structure (operon organization), the organization of genes in pathways (KEGG) and the gene function with Gene Ontology (GO) terms. The integrated viewers can be used in parallel and combine expression data and functional annotations from different external data repositories. The graphical visualizations evidently accelerate both the validation of regulatory information and the detection of affected biological processes. AVAILABILITY: http://leger2.gbf.de/cgi-bin/vis-o-bac.pl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chromosome Mapping↗

PubMed2026

LLPS-based classification and a novel prognostic signature reveal NRF1 as a therapeutic target in pancreatic cancer.

BACKGROUND: Aberrant liquid-liquid phase separation (LLPS) can alter biomolecular condensate functions and may influence pancreatic tumorigenesis and progression, but the specific role of LLPS regulators in prognosis and the tumor immune microenvironment (TIME) in pancreatic ductal adenocarcinoma (PDAC) remains unclear. METHODS: We integrated transcriptome data of LLPS regulator-related differentially expressed genes (DEGs; n = 298) in a cohort of 176 PDAC patients from TCGA. Three LLPS regulator subtypes (LS1-LS3) were identified through multi-omics analyses, and a prognostic LLPS subtype-related risk model (LRRPC) was developed and validated. Chromatin immunoprecipitation confirmed NRF1 binding to promoters of key risk genes, and in vitro and in vivo experiments assessed the effects of NRF1 targeting on tumor growth. RESULTS: The three LLPS regulator subtypes exhibited significant differences in prognosis, clinical features, genomic alterations, TIME patterns and predicted immunotherapy response. The LRRPC signature predicted prognosis and immunotherapy efficacy across cohorts and was associated with tumor biomarkers and immune infiltration. Nuclear Respiratory Factor 1 (NRF1) directly regulated hub genes such as FAM83A, RHOV and ITGB6, promoting PDAC cell proliferation, while its inhibition induced apoptosis and reduced tumor growth. CONCLUSIONS: This study proposes an LLPS-based stratification framework for PDAC, and the LRRPC model provides an LLPS subtype-related risk score that may assist personalized prognostic assessment and immunotherapy stratification. NRF1 emerges as a promising therapeutic candidate whose targeting can inhibit tumor progression in PDAC experimental models and warrants further evaluation.

Immunotherapy↗

PubMed2026

Multi-omics characterization of flavor profile differences in the Longissimus thoracis between Angus and Hereford cattle.

BACKGROUND: Angus and Hereford cattle are premier breeds widely used in genetic improvement and crossbreeding programs to enhance meat quality, yet the flavor differences between them remain poorly understood. RESULTS: In this study, we performed an integrated analysis using headspace solid-phase microextraction coupled with gas chromatography-mass spectrometry (HS-SPME-GC-MS)-based volatile metabolomics, lipidomics, and untargeted metabolomics to characterize the flavor profiles of the Longissimus thoracis (LT) muscle from both breeds and to identify potential precursor substances underlying flavor formation. In total, we identified 76 differential volatile organic compounds (VOCs) among the 493 candidate VOCs. By combing relative odor activity value (ROAV) and sensory attribute annotation, 2,3-butanedione which may contribute to the creamy aroma was revealed as the core differential VOC between the two breeds. This finding was robustly validated across SHAP (Shapley additive explanations) analysis, KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment, and flavor annotation. Lipidomic analysis revealed 689 differential lipids primarily belonging to classes such as phosphatidylcholine, triglycerides, and phosphatidylethanolamine. Correlation analysis further linked these lipid profiles to flavor, showing that fatty acids (FAs) including FA(19:0), FA(18:2 + O), FA(14:1), FA(16:1), and FA(14:0) were significantly correlated with 2,3-butanedione. Notably, the unsaturated fatty acids (UFAs) in the Longissimus thoracis (LT) of Hereford cattle exhibited higher double bond content compared to Angus cattle, suggesting a greater potential for rich flavor development. Untargeted metabolomics revealed that nine of the 9474 metabolites were significantly correlated with both 2,3-butanedione and FAs, including norepinephrine, l-beta-aspartyl-l-leucine, and artemetin. CONCLUSIONS: Overall, our research has identified differential flavor compounds and potential precursor substances between Angus cattle and Hereford cattle, providing targeted guidance for breed improvement. © 2026 Society of Chemical Industry.

2,3‐butanedione↗

PubMed2026

Multi-omics identification of ZFPM2 and CD44 as key candidates for testicular size in sheep.

Testicular size is a key determinant of ram fertility, yet its genetic architecture in sheep remains poorly understood. Crossbred sheep exhibit substantial testicular developmental variation and serve as ideal models for screening fertility-related genes. Here, we performed whole-genome sequencing (WGS)-based genome-wide association study (GWAS) on 115 rams, including seven crossbred populations (each derived from a distinct sire breed crossed with Hu sheep) and a purebred Hu sheep population. We identified two variants within ZFPM2, including an intronic single nucleotide polymorphism (SNP) rs421073404 and a missense SNP rs1086332841, both significantly correlated with testis weight, length, and width. Importantly, these variants exerted testis-specific effects, consistent with their weak correlation with body weight (r < 0.2). Genomic selection scans (FST and π-ratio) between rams with large (>150 g/side) and small (<75 g/side) testes revealed multiple divergent genomic regions, with prominent haplotype differentiation at the CD44 locus on chromosome 15. Notably, two linked missense variants (rs160734053 and rs414318703) in CD44 exerted opposing effects on testicular size, suggesting allelic heterogeneity at this locus. Integrative RNA-seq and ATAC-seq across 0-12 months further uncovered prepubertal stage-specific expression and chromatin accessibility patterns of ZFPM2 and CD44, providing mechanistic insights into their regulatory roles. Collectively, this study provides candidate SNPs for sheep marker-assisted selection and multi-omic evidence for dissecting the genetic basis of testicular development in ovine breeding.

Animals↗

PubMed2026

engGNN: a dual-graph neural network for omics-based disease classification and feature selection.

Omics data, such as transcriptomics, proteomics, and metabolomics, provide critical insights into disease mechanisms and clinical outcomes. However, their high dimensionality, small sample sizes, and intricate biological networks pose major challenges for reliable prediction and meaningful interpretation. Graph neural networks offer a promising way to integrate prior knowledge by encoding feature relationships as graphs. Yet, existing methods typically rely solely on either an externally curated feature graph or a data-driven generated graph, which limits their ability to capture complementary information. To address this, we propose the external and generated Graph Neural Network (engGNN), a dual-graph framework that jointly leverages both external biological networks and data-driven generated graphs. Specifically, engGNN constructs a biologically informed undirected feature graph from established network databases and complements it with a directed feature graph derived from tree-ensemble models. This dual-graph design produces more comprehensive representations, thereby improving predictive performance and interpretability. Through extensive simulation studies and real-world applications to three independent gene expression datasets, engGNN consistently demonstrates strong classification performance compared with competitive baselines. Beyond classification, engGNN provides feature- and source-level interpretability, enabling biologically meaningful analyses such as pathway enrichment analysis. Taken together, these results highlight engGNN as a robust, flexible, and interpretable framework for disease classification and biomarker discovery in high-dimensional omics contexts.

Graph Neural Networks↗

PubMed2026

Integrated histone and proteome analyses reveal convergent and distinct hepatotoxic mechanisms of tenuazonic acid and deoxynivalenol.

Mycotoxins are widespread dietary contaminants whose health impacts are expected to intensify under climate change. Although their mechanisms of toxicity remain incompletely understood, epigenetic dysregulation has been increasingly implicated. Here, mass spectrometry-based multi-omics was used to profile histone post-translational modifications and proteome dynamics in HepG2 cells exposed to seven mycotoxin conditions. Time-resolved analyses identified tenuazonic acid as the dominant cellular disruptor, inducing alterations in H3K27 and H1 variants, and revealing a previously unrecognized oxidative modification of the H1.0 N-terminal methionine (H1.0N-term0AcM0Ox) that retains the protein's N-terminal acetylation. An Alternaria toxin mixture induced similar H1 responses, largely driven by tenuazonic acid, while deoxynivalenol produced convergent chromatin and proteomic alterations. Proteomic remodeling was characterized by increased protein translation, reduced mitochondrial complex IV expression, and impaired cholesterol biosynthesis, whereas sterigmatocystin activated DNA replication and repair pathways. Together, these findings demonstrate that mycotoxins disrupt chromatin organization, protein synthesis, and lipid metabolism, providing toxicological insight into hepatocellular dysfunction. These findings warrant further validation and mechanistic investigation in future hypothesis-driven studies of mycotoxin exposure.

Trichothecenes↗

PubMed2026

Identification and genetic validation of potential therapeutic targets for pulmonary hypertension through multi-omics causal inference.

Pulmonary hypertension (PH) underscores the urgent need for novel therapeutic targets. This study aimed to employ a proteome-wide Mendelian randomization (MR) approach to systematically identify circulating proteins causally associated with PH, thereby providing genetically validated candidate targets for drug development. We adopted a 2-sample MR design, integrating large-scale plasma proteomic quantitative trait loci (pQTL) data (encompassing 4148 proteins) and summary statistics from a large-scale PH genome-wide association study (2047 cases, 8301 controls). Candidate targets were screened through a multilayered analytical pipeline comprising proteomic MR, transcriptomic MR, and summary-data-based Mendelian randomization. The ultimately identified MR-Identified Causal Candidate Targets (MR-ICTs) underwent rigorous Bayesian colocalization analysis, followed by biological characterization through functional enrichment analysis, single-cell transcriptomics, and phenome-wide association studies. Through robust genetic causal inference, this study provides that circulating proteins such as LYZ, GREM2, NID1, and PF4V1 play causal roles in PH pathogenesis. These findings offer a set of rigorously genetically validated, high-priority therapeutic targets for developing novel PH treatments, specifically addressing key pathological mechanisms such as innate immunity, BMP signaling pathway dysregulation, and platelet activation. Our multi-dimensional analysis ultimately identified 6 MR-ICTs causally associated with PH. Notably, the causal associations for lysozyme C (LYZ), gremlin-2 (GREM2), nidogen-1 (NID1), and platelet factor 4 variant 1 (PF4V1) were stringently validated by Bayesian colocalization analysis (posterior probability for hypothesis 4 [PPH4], indicating a shared causal variant, > 0.99). Functional enrichment analysis revealed significant involvement of these targets in immune response and TGF-β signaling pathways. Single-cell analysis further elucidated their cell-type-specific expression, with LYZ predominantly expressed in monocytes and PF4V1 almost exclusively in platelets.

Hypertension, Pulmonary↗

PubMed2026

Epigenetics and In Silico Transcriptome Analysis of Pediatric Acute Myeloid Leukemia.

Pediatric acute myeloid leukemia (AML) is a heterogeneous hematologic malignancy that accounts for about 15%-20% of childhood leukemias. Despite therapeutic advances, relapses remain common, and survival for high-risk patients is below 60%. Unlike adult AML, pediatric AML displays distinct genetic mutations, including FLT3-ITD, NPM1, KMT2A rearrangements, and core-binding factors (CBF) fusions, as well as extensive epigenetic dysregulation. Aberrant DNA methylation, histone modifications, and altered non-coding RNA expressions disrupt hematopoietic differentiation and activate oncogenic transcriptional networks. Recent advances in silico transcriptomic analysis have transformed the study of pediatric AML by integrating gene expression and epigenetic data to identify molecular drivers and regulatory networks. Computational RNA-seq pipelines and pathway analyses have highlighted key epigenetic regulators, including DNMT3A, TET2, and HDACs, as potential therapeutic targets. Multi-omics approaches combining transcriptomic, methylomic, and chromatin accessibility data are increasingly used to define biomarkers for diagnosis, prognosis, and therapeutic response. This review provides a comprehensive overview of the molecular and epigenetic landscape of pediatric AML, emphasizing the power of in silico transcriptome analysis to uncover disease mechanisms, refine patient stratification, and guide the development of precision-based epigenetic therapies aimed at improving long-term outcomes in children with AML.

Humans↗

PubMed2026

A Landscape of Drosophila melanogaster Disease Models: From Genetic Platforms to Cross-Disease Mechanisms and Translational Research.

Modeling human diseases using the fruit fly (Drosophila melanogaster) has established itself as a cornerstone of functional genomics and preclinical medicine. Despite its anatomical simplicity, the Drosophila genome shares remarkable functional conservation with human disease-related genes, enabling the study of complex physiological traits through accessible tissue models. Furthermore, beyond individual disease models, we propose a framework demonstrating how these diseases converge at common molecular centers, such as the breakdown of protein homeostasis, mitochondrial dysfunction, chronic inflammation, and organ-to-organ communication. Finally, we discuss strategies for integrating the Drosophila platform into drug development pipelines and establishing standards to enhance inter-laboratory reproducibility. Overall, this review highlights the enduring value of fruit flies as a model system, particularly when combined with AI-omics approaches to transform complex biological datasets into actionable therapeutic strategies.

Drosophila↗

PubMed2026

A nucleolar stress gene signature enables quantitative scoring across multi-omics contexts.

The nucleolus is essential for ribosome biogenesis and cellular homeostasis, and its dysfunction can induce nucleolar stress, a process implicated in cancer and other diseases. However, nucleolar stress is commonly inferred from morphological changes or a limited set of functional assays, and quantitative approaches based on gene expression profiles remain lacking. Here, we integrate literature curation with multi-dataset screening to define a nucleolar stress gene signature and develop a nucleolar stress score (NuS) applicable to bulk transcriptomics, single-cell transcriptomics, proteomics, and spatial transcriptomics. Using this framework, we show in colorectal cancer models that oxaliplatin induces nucleolar stress, suppresses nascent rRNA synthesis, and activates p53 signaling, whereas these responses are attenuated in oxaliplatin-resistant cells. Combined with a ribosome biogenesis activity score (RiboSis), NuS captures related but distinct dimensions of nucleolar function and stratifies tumors into functional states associated with clinical outcomes. NuS-based analysis of perturbational transcriptomes further prioritizes compounds with putative nucleolar stress-inducing activity. Collectively, this study provides a quantitative framework for evaluating nucleolar stress and illustrates its applications in disease stratification and drug mechanism discovery.

Cell Nucleolus↗

PubMed2025

Roles of microbial interactions in determining the establishment and function of synthetic consortium inoculants for soil applications.

Synthetic microbial consortium inoculants are emerging nature-based solutions for promoting sustainable agriculture and mitigating environmental challenges. However, despite promising results in simpler lab-scale trials, many inoculants fail to establish or perform satisfactorily in field conditions. One most critical yet least understood factor influencing inoculant effectiveness is the complex microbial interactions, both within consortium inoculants ("within-community" interactions) and between consortium inoculants and native soil communities ("cross-community" interactions). Here, we first discuss major negative and positive "within-community" interactions and highlight the importance to design consortium inoculants with positive interactions for improved stability and functionality. We then examine the bidirectional "cross-community" interactions once introducing consortium inoculants to soils. Soil native communities often create strong resistance to the invasion of inoculants. We discuss major drivers controlling the invasibility of native communities and various strategies increasing the invasiveness of consortium inoculants. We then discuss how consortium inoculants can reshape native communities, with implications for long-term ecosystem resilience and functioning. We propose future research efforts including advancing strategies for harnessing natural species from relatively untapped soil reservoirs and using high-throughput interaction profiling with multi-omics and computational tools to build compatible synthetic consortia with desirable functions; leveraging positive interactions and prebiotics to facilitate inoculant establishment; and assessing fully soil functional resilience over longer terms, including recognizing the importance of rare keystone taxa. By integrating with ecological theory, this review provides a comprehensive insight into microbial interactions to advance the design, application, and monitoring of synthetic consortium inoculants for enhancing soil health and ecosystem sustainability.

establishment↗

PubMed2026

Integrated methylome and transcriptome analysis provides insight into DNA methylation-mediated networks in sexual dimorphism of Vernicia montana.

BACKGROUND: Sexual dimorphism is fundamental to reproduction in dioecious plants and is regulated by both genetic and epigenetic mechanisms. DNA methylation is a central epigenetic mark known to influence phenotypic variation in plants. However, its specific role in shaping sexual dimorphism in dioecious trees remains poorly understood. To address this question, we performed integrated genome-wide DNA methylome and transcriptome analyses of four tissue types in the dioecious tung tree (Vernicia montana), including male and female flower buds and their corresponding leaves. RESULTS: Our analysis revealed distinct DNA methylation patterns between male and female tissues. Notably, the coordination between DNA methylation reprogramming and transcriptional regulation appeared to be more strongly associated with reproductive development than with vegetative growth in V. montana. We identified a set of sex-biased genes that may reflect different reproductive strategies between the sexes. Further analysis identified several key transcription factors (TFs) potentially associated with promoter differentially methylated regions (DMRs), including flowering-time regulators (e.g., FRS5, REM16, and VRN1) and TFs involved in hormone signaling pathways such as jasmonic acid, auxin, and salicylic acid signaling. Cis-regulatory element analysis showed that some promoter DMRs overlapped with hormone response elements related to abscisic acid, auxin, and gibberellin. Co-expression network analysis further revealed potential regulatory correlations among promoter DMR-mediated TFs, hormone-responsive pathways, and key floral development regulators. CONCLUSIONS: Collectively, our results suggest that interactions among DNA methylation, transcriptional regulation, and hormone-responsive pathways may contribute to the establishment of sexual dimorphism in V. montana. This study provides the first integrated view of these regulatory layers in V. montana and supports a species-specific regulatory framework for understanding the epigenetic basis of sexual dimorphism in this economically important dioecious tree. The proposed framework is based on multi-omics analyses and warrants further validation through targeted functional studies.

DNA Methylation↗