Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “omics analysis”

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 271 records · Page 15Linked to original sources

Systematic discovery of retina-enriched Rik genes identifies 1190005I06Rik as a novel modulator of visual signalling.

BACKGROUND: High‑throughput transcriptome projects have revealed thousands of mammalian genes with little or no functional annotation. Among these are hundreds of loci assigned provisional “Rik” identifiers following discovery in the RIKEN cDNA annotation effort. Although often dismissed as genomic dark matter, such genes may encode tissue‑restricted proteins that modulate physiologic functions and influence disease. The retina is a highly specialised neural tissue and a common site of inherited disorders; understanding its molecular repertoire could illuminate novel therapeutic avenues. METHODS: We integrated bulk RNA‑seq from ten adult mouse tissues, evolutionary and domain analysis, single‑cell RNA‑seq, and CRISPR/Cas9 gene disruption to systematically catalogue protein‑coding Rik genes enriched in the retina and test the function of a representative gene. RESULTS: A rigorous differential expression analysis identified 44 Rik genes with robust retina‑specific expression compared with nine non‑retinal tissues. Many of these genes lack orthologues beyond rodents, while others show broad conservation, illustrating a continuum from lineage‑restricted to conserved retinopathy candidates. Single‑cell transcriptomics revealed that these genes are expressed across retinal cell types, with the highest aggregate expression in cone photoreceptors and inner interneurons. To evaluate physiological significance, we generated a 1190005I06Rik knockout mouse. Although retinal architecture appeared normal, loss of 1190005I06Rik enhanced electroretinogram b‑wave amplitudes and altered light‑avoidance behaviour, indicating that this previously uncharacterised gene acts as a negative modulator of visual signalling. CONCLUSIONS: We present a curated atlas of retina‑enriched Rik genes and demonstrate that 1190005I06RIK modulates retinal circuit function. This resource expands the molecular landscape of the retina and provides new candidates for the genetic basis of inherited retinal disease. Our findings underscore that unannotated genes may exert measurable effects on sensory processing and warrant systematic exploration in the context of human ocular disorders.

Animals↗

Integromic analysis of the NCI-60 cancer cell lines.

Microarray-based transcript profiling has become exceedingly popular, particularly for breast cancer. However, other 'omic' profiling technologies at the DNA, RNA, protein, functional, and pharmacological levels are also becoming increasingly practical. We define 'integromics' as the melding of such diverse types of data from different experimental platforms. The whole can sometimes be more than the sum of its parts. We describe here a set of integromic studies in which we have profiled the 60 human cancer cell lines (the NCI-60) used by the National Cancer Institute to screen >100,000 chemical compounds over the last 13 years. Patterns of potency in the screen can be mapped into molecular structures of the compounds or into molecular characteristics of the cells. Here we discuss conceptual and experimental aspects of the profiling, as well as a number of bioinformatic computer programs (CIMminer, MedMiner, MatchMiner, and GoMiner) that we have developed for biological interpretation of the profiles. As briefly reviewed here, we have used the combination of NCI-60 data types to identify markers for distinguishing tumor types and to obtain pharmacogenomic clues for possible individualization of a cancer therapy.

Journal Article↗

Multi-omics uncovers the pleiotropic genetic mechanisms linking MASLD and cardiometabolic syndromes.

BACKGROUND: Metabolic dysfunction-associated steatotic liver disease (MASLD) and cardiovascular-kidney-metabolic (CKM) syndrome are interrelated conditions with shared pathophysiological features; however, the genetic architecture underlying their relationship has not been fully elucidated. Deciphering this shared genetic basis holds promise for advancing mechanistic insights and therapeutic discovery. METHODS: We performed an integrated genome-wide cross-trait analysis using GWAS summary statistics for MASLD and 38 CKM traits. Our analysis estimated genetic correlations, inferred causal relationships, and identified pleiotropic variants. Candidate causal genes and druggable targets were subsequently prioritized through integrating multi-omics data. RESULTS: MASLD exhibited significant genetic correlations with 16 CKM traits, especially metabolic and cardiovascular conditions. Bidirectional causal relationships were observed between MASLD and T2D, adiposity, and lipid traits. We discovered 116 pleiotropic loci, including 65 shared causal variants such as rs429358 near APOE, which exerted influence across multiple traits. Gene-based analyses prioritized 152 unique candidate pleiotropic genes, enriched in lipid and cholesterol metabolism, and highly expressed in the liver, adipose, and immune-related cell types, such as macrophages and endothelial cells. Multi-omics integration validated 131 genes using eQTL and pQTL data from multiple tissues and cohorts. Notably, FTO and APOE emerged as central pleiotropic hubs, and druggability evaluation highlighted APOE, LPL, PPARG, and GPBAR1 as established therapeutic targets for metabolic diseases. CONCLUSION: This study provides a comprehensive map of the shared genetic architecture between MASLD and CKM syndrome, reveals novel causal genes and repurposable drug targets, and offers insights into precision medicine approaches for cardiometabolic and liver diseases.

Humans↗

Urinary multi-omics reveal non-invasive diagnostic biomarkers in clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most common kidney malignancy. Yet, no rapid, non-invasive biomarkers are available for diagnosis or screening. Urine represents an ideal analyte matrix due to its accessibility, low invasiveness, longitudinal sampling, and the kidney's central role in filtration. Here, we integrated proteomic, lipidomic, and metabolomic analyses of urine from ccRCC patients and controls to identify diagnostic biomarkers. Multi-omics profiling revealed urogenital metabolic dysregulation in ccRCC, including increased lipid metabolism, altered mitochondrial respiration signatures, and elevated urinary lipid content. We identified three urinary protein biomarkers: serum amyloid A1 (SAA1), haptoglobin (HP), and lipocalin 15 (LCN15). Using a parallel reaction monitoring mass spectrometry workflow, we developed a rapid and sensitive assay and combined these markers into a diagnostic UrineScore. The UrineScore achieved 0.96 in an area under the receiver operating characteristic curve analysis in the discovery cohort, and 0.95 in an independent validation cohort. Together, these results support the feasibility of multi-omics-guided urinary biomarker discovery and represent a step toward accessible diagnostic platforms for ccRCC.

Humans↗

Multi-omics integration uncovers epigenetic control of metabolic reprogramming in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive subtype characterized by the absence of estrogen, progesterone, and HER2 receptors, limiting effective targeted therapies. Increasing evidence suggests that metabolic reprogramming, a hallmark of TNBC progression, is driven by underlying epigenetic mechanisms such as DNA methylation. The represented study performed an integrative analysis of transcriptomic (RNA-seq) and methylome data to uncover the metabolic-epigenetic interplay in TNBC. Differential gene expression analysis using DESeq2 revealed significant dysregulation of key metabolic genes, including upregulation of genes encoding glycolytic and serine biosynthesis enzymes and downregulation of metabolic tumor suppressors. Genome-wide methylation profiling identified extensive cytosine-phosphate-guanine (CpG) hypermethylation events associated with transcriptional repression, particularly in promoter regions. Integrative analysis pinpointed a subset of metabolism-related genes exhibiting both differential expression and methylation, such as FBP1, RASSF1A, and PHGDH. Pathway enrichment analysis highlighted aberrations in glycolysis/gluconeogenesis, fatty acid metabolism, and one-carbon pathways (adjusted p&#x2009;<&#x2009;0.01). Importantly, TNBC patients with hypermethylated metabolic gene signatures displayed significantly shorter overall survival (log-rank p&#x2009;<&#x2009;0.05). These findings reveal that DNA methylation-driven metabolic dysregulation contributes to TNBC aggressiveness and may provide novel biomarkers and therapeutic targets at the metabolic-epigenetic interface.

Humans↗

Biological master games: using biologists' reasoning to guide algorithm development for integrated functional genomics.

We review some powerful new algorithms that build on the intuitive biological interpretation techniques for statistical analysis of functional genomics experiments. Although they were originally designed for transcriptomics, we argue that these algorithms are applicable to any type of -omics study (transcriptomics, proteomics, metabolomics). Rank Products (RP), a strictly non-parametric test statistic to detect differentially regulated elements (genes, proteins, metabolites) in genome-wide screens. RP is particularly powerful for noisy data and low numbers of replicates and makes full use of the availability of a large number of parallel measurements that is typical of modern large-scale experiments. Iterative Group Analysis (iGA), a statistical method that makes the transition from regulated single elements to significant classes of elements, and thus provides an automatic functional annotation of an experiment. Graph-based iGA (GiGA), an extension of iGA that combines experimental data with a broad variety of biological annotations to highlight physiologically relevant regions in a given "evidence graph" (e.g., metabolic networks, signaling pathway diagrams, protein interaction maps). The sequential application of these techniques yields an increasingly abstract interpretation of experimental data that is at the same time quantitative, statistically rigorous, and biologically significant. The results can be used either as helpful tools to guide data visualization and exploration, or as the input for downstream computational applications in a systems biology framework.

Algorithms↗

Early-Onset Atrial Fibrillation and the Prevalence of Rare Variants in Cardiomyopathy and Arrhythmia Genes.

IMPORTANCE: Early-onset atrial fibrillation (AF) can be the initial manifestation of a more serious underlying inherited cardiomyopathy or arrhythmia syndrome. OBJECTIVE: To examine the results of genetic testing for early-onset AF. DESIGN, SETTING, AND PARTICIPANTS: This prospective, observational cohort study enrolled participants from an academic medical center who had AF diagnosed before 66 years of age and underwent whole genome sequencing through the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine program. Participants were enrolled from November 23, 1999, to June 2, 2015. Data analysis was performed from October 24, 2020, to March 11, 2021. EXPOSURES: Rare variants identified in a panel of 145 genes that are included on cardiomyopathy and arrhythmia panels used by commercial clinical genetic testing laboratories. MAIN OUTCOMES AND MEASURES: Sequencing data were analyzed using an automated process followed by manual review by a panel of independent, blinded reviewers. The primary outcome was classification of rare variants using American College of Medical Genetics and Genomics criteria: benign, likely benign, variant of undetermined significance, likely pathogenic, or pathogenic. Disease-associated variants were defined as pathogenic/likely pathogenic variants in genes associated with autosomal dominant or X-linked dominant disorders. RESULTS: Among 1293 participants (934 [72.2%] male; median [interquartile range] age at enrollment, 56 [48-61] years; median [interquartile range] age at AF diagnosis, 50 [41-56] years), genetic testing identified 131 participants (10.1%) with a disease-associated variant, 812 (62.8%) with a variant of undetermined significance, 92 (7.1%) as heterozygous carriers for an autosomal recessive disorder, and 258 (20.0%) with no suspicious variant. The likelihood of a disease-associated variant was highest in participants with AF diagnosed before the age of 30 years (20 of 119 [16.8%; 95% CI, 10.0%-23.6%]) and lowest after the age of 60 years (8 of 112 [7.1%; 95% CI, 2.4%-11.9%]). Disease-associated variants were more often associated with inherited cardiomyopathy syndromes compared with inherited arrhythmias. The most common genes were TTN (n&#x2009;=&#x2009;38), MYH7 (n&#x2009;=&#x2009;18), MYH6 (n&#x2009;=&#x2009;10), LMNA (n&#x2009;=&#x2009;9), and KCNQ1 (n&#x2009;=&#x2009;8). CONCLUSIONS AND RELEVANCE: In this cohort study, genetic testing identified a disease-associated variant in 10% of patients with early-onset AF (the percentage was higher if diagnosed before the age of 30 years and lower if diagnosed after the age of 60 years). Most pathogenic/likely pathogenic variants are in genes associated with cardiomyopathy. These results support the use of genetic testing in early-onset AF.

Adult↗

Unintended effects and their detection in genetically modified crops.

The commercialisation of GM crops in Europe is practically non-existent at the present time. The European Commission has instigated changes to the regulatory process to address the concerns of consumers and member states and to pave the way for removing the current moratorium. With regard to the safety of GM crops and products, the current risk assessment process pays particular attention to potential adverse effects on human and animal health and the environment. This document deals with the concept of unintended effects in GM crops and products, i.e. effects that go beyond that of the original modification and that might impact primarily on health. The document first deals with the potential for unintended effects caused by the processes of transgene insertion (DNA rearrangements) and makes comparisons with genetic recombination events and DNA rearrangements in traditional breeding. The document then focuses on the potential value of evolving "profiling" or "omics" technologies as non-targeted, unbiased approaches, to detect unintended effects. These technologies include metabolomics (parallel analysis of a range of primary and secondary metabolites), proteomics (analysis of polypeptide complement) and transcriptomics (parallel analysis of gene expression). The technologies are described, together with their current limitations. Importantly, the significance of unintended effects on consumer health are discussed and conclusions and recommendations presented on the various approaches outlined.

Animals↗

Markers of microvascular instability predict severity and survival in idiopathic pulmonary fibrosis.

INTRODUCTION: Most research on idiopathic pulmonary fibrosis (IPF) has focused on the interplay among fibroblasts, the immune system and epithelial cells. There is growing evidence that microvascular dysfunction also plays a role in disease progression, but large human translational studies are lacking. In this research, we aim to identify a proteomic signature of microvascular instability and assess the impact of current therapeutics on the microvasculature. METHODS: Olink proteomic data from patients with IPF were obtained from the Pulmonary Fibrosis Foundation Patient Registry (PFF-PR) (n=914) and an independent validation cohort (n=366). Among the PFF-PR, 640 patients also have whole-blood RNA sequencing data available. A subset of 79 microvascular-associated proteins was curated, and their associations with disease severity and transplant-free survival were examined. An adaptive least absolute shrinkage and selection operator was used to generate a novel microvascular risk score. RESULTS: Higher plasma levels of five microvascular-associated proteins (SDC1, MMP10, THBS2, HGF and SERPINA5) were associated with lung function and survival in both cohorts. Whole-blood RNA sequencing of patients with microvascular risk revealed enrichment of immune-mediated processes. Patients with higher microvascular risk who were subsequently put on nintedanib in the following year had significantly better 3-year transplant-free survival compared with patients who did not receive antifibrotic intervention (HR 0.56, 95%&#x2009;CI 0.35 to 0.89, p=0.0142). DISCUSSION: Integrative multi-omics analyses suggest that perturbations to microvascular remodelling contribute to disease severity and progression in IPF. This analysis offers a framework for a precision medicine approach for IPF.

Idiopathic pulmonary fibrosis↗

Genomic analysis of Dasiphora on the Qinghai-Tibet Plateau provides insights into genetic divergence and flower color variation.

The Qinghai-Tibet Plateau (QTP) harbors diverse alpine flora, including the ecologically significant shrubs Dasiphora fruticosa and D. glabra, for which taxonomic uncertainties remain and adaptive mechanisms are still poorly understood. Based on high-quality genome assembly, population resequencing, and multi-omics integration, we elucidated their evolutionary divergence and flower color genetics. Chromosome-level haplotype-resolved genomes were assembled: autotetraploid D. fruticosa (929.99&#x2009;Mb) and diploid D. glabra (450.89&#x2009;Mb). Phylogenetic analysis showed that the tetraploid D. fruticosa and D. glabra in this study clustered together, while the diploid D. fruticosa sequenced by previous research formed a distinct lineage clustered outside. Consistently, population structure analysis of 55 samples revealed three major clades, with D. fruticosa further subdivided into two divergent branches. Additionally, hybridization events detected by Admixture, coupled with ploidy complexity identified via flow cytometry highlight the intricate genetic relationships within this genus. Adaptive gene families expanded in antioxidant (flavonoid synthesis) and secondary metabolism pathways, adapting to ultraviolet radiation and cold stress. Natural selection analysis identified 193 candidate genes (e.g., TFB5 in the DNA repair pathway), predominantly localized to chromosome 5, which are potential candidates for high-altitude adaptation. Transcriptome and metabolome analyses showed D. fruticosa's yellow petals derive from flavonol (quercetin) accumulation, while D. glabra's white petals result from proanthocyanidin biosynthesis via high LAR/ANR expression. This study provides insights into the taxonomic revision and adaptive genetic divergence of alpine plants, and offers a foundation for horticultural improvement of Dasiphora.

Flowers↗

Metabolic engineering in the -omics era: elucidating and modulating regulatory networks.

The importance of regulatory control in metabolic processes is widely acknowledged, and several enquiries (both local and global) are being made in understanding regulation at various levels of the metabolic hierarchy. The wealth of biological information has enabled identifying the individual components (genes, proteins, and metabolites) of a biological system, and we are now in a position to understand the interactions between these components. Since phenotype is the net result of these interactions, it is immensely important to elucidate them not only for an integrated understanding of physiology, but also for practical applications of using biological systems as cell factories. We present some of the recent "-omics" approaches that have expanded our understanding of regulation at the gene, protein, and metabolite level, followed by analysis of the impact of this progress on the advancement of metabolic engineering. Although this review is by no means exhaustive, we attempt to convey our ideology that combining global information from various levels of metabolic hierarchy is absolutely essential in understanding and subsequently predicting the relationship between changes in gene expression and the resulting phenotype. The ultimate aim of this review is to provide metabolic engineers with an overview of recent advances in complementary aspects of regulation at the gene, protein, and metabolite level and those involved in fundamental research with potential hurdles in the path to implementing their discoveries in practical applications.

Bacteria↗

Giotto Suite: a multiscale and technology-agnostic spatial multiomics analysis ecosystem.

Emerging spatial multiomics technologies provide an increasingly large amount of information content at multiple scales. However, it remains challenging to efficiently represent and harmonize diverse spatial datasets. Here we present Giotto Suite, a suite of modular packages that provides scalable and extensible end-to-end solutions for multiscale and multiomic data analysis, integration and visualization. At its core, Giotto Suite is centered around an innovative data framework, allowing the representation and integration of spatial omics data in a technology-agnostic manner. Giotto Suite integrates molecular, morphology, spatial and annotated feature information to create a responsive and flexible workflow, as demonstrated by applications to several state-of-the-art spatial technologies. Furthermore, Giotto Suite builds upon interoperable interfaces and data structures that bridge the established fields of genomics and spatial data science in R, thereby enabling independent developers to create custom-engineered pipelines. As such, Giotto Suite creates an immersive and multiscale ecosystem for spatial multiomic data analysis.

Genomics↗

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

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

Animals↗

NeuroOmics-Net: An interpretable multimodal deep learning framework for Alzheimer's disease diagnosis and progression prediction using neuroimaging, EEG, and genomic data.

Accurate diagnosis and progression prediction of Alzheimer's disease (AD) remain challenging due to the heterogeneous nature of the disease, which involves structural brain degeneration, electrophysiological dysfunction, and molecular dysregulation. Most existing deep learning approaches rely on a single modality or limited multimodal combinations, thereby failing to capture the complex cross-domain interactions underlying AD progression. Furthermore, the scarcity of large-scale datasets containing synchronized neuroimaging, electrophysiological, and genomic measurements restricts the development of comprehensive multimodal diagnostic systems. To address these challenges, this study proposes NeuroOmics-Net, a multimodal deep learning framework for Alzheimer's disease analysis that integrates structural magnetic resonance imaging (sMRI), electroencephalography (EEG), and gene expression data. The proposed framework combines a Hierarchical Multi-View Encoder (HME) for modality-specific feature extraction, a Cross-Omics Attention Fusion (CAF) module for adaptive integration of complementary biomarkers, and a Disease Progression Graph Learning (DPGL) module for modeling progression-related relationships across biological domains. To facilitate cross-modal integration from independent cohorts, Regularized Canonical Correlation Analysis (RCCA) is employed to align heterogeneous feature representations within a shared latent space. Experiments were conducted using publicly available datasets from ADNI, PhysioNet, and GEO repositories comprising 1120 diagnosis-aligned samples. The proposed framework achieved 94.3% classification accuracy and an AUC of 0.975 for distinguishing normal controls (NC), mild cognitive impairment (MCI), and Alzheimer's disease subjects, while attaining 93.7% accuracy for predicting conversion from stable mild cognitive impairment (sMCI) to progressive mild cognitive impairment (pMCI). However, a fairness sensitivity analysis using stratified demographic reweighting revealed accuracy ranging from 90.8% (low-education, high-comorbidity proxy subgroup) to 96.1% (low-risk, high-reserve proxy subgroup), a demographic parity gap of 5.3 percentage points, indicating that overall accuracy reflects a performance ceiling in a relatively homogeneous research cohort rather than a realistic estimate for demographically diverse clinical populations. Comparative evaluations demonstrated consistent improvements over state-of-the-art unimodal and multimodal deep learning models. Interpretability analysis further identified clinically relevant biomarkers, including hippocampal and entorhinal atrophy, theta-alpha EEG alterations, and APOE-associated molecular pathways. Because sMRI, EEG, and gene expression data were sourced from separate, unpaired cohorts with no subjects possessing all three synchronized measurements, all reported cross-modal associations reflect population-level statistical correspondence across diagnosis-matched groups rather than within-subject physiological coupling; no claim of intra-individual causal cross-modal interaction is made. These findings demonstrate that NeuroOmics-Net provides an effective computer-aided framework for multimodal biomedical data processing and Alzheimer's disease analysis. By integrating neuroimaging, electrophysiological, and genomic information, the proposed approach enables accurate diagnosis, progression prediction, and biologically interpretable decision support for clinical and translational applications.

Humans↗

Genome-wide DNA methylation regulation analysis provides novel insights on post-radiation breast cancer.

Breast cancer (BC) is the most common malignancy with a poor prognosis. Radiotherapy is one of the leading traditional treatments for BC. However, radiotherapy-associated secondary diseases are severe issues for the treatment of BC. The present study integrated multi-omics data to investigate the molecular and epigenetic mechanisms involved in post-radiation BC. The differences in the expression of radiation-associated genes between post-radiation and pre-radiation BC samples were determined. Enrichment analysis revealed that these radiation-associated genes involved diverse biological functions and pathways in BC. Combining epigenetic data, we identified radiation-associated genes whose transcriptional changes might be associated with aberrant methylation. Then, we identified potential therapeutic targets and chemical drugs for post-radiation BC patient treatment by constructing a drug-target association network. Specifically, four radiation-associated genes (CD248, CCDC80, GADD45B, and MMP2) whose increased expression might be regulated by hypomethylation of the corresponding enhancer region were found to have excellent diagnostic effects and clinical prognostic value. Finally, we further used independent samples to verify CD248 expression and established a simple epigenetic regulatory model. In summary, this study provides novel insights for understanding the regulation of target genes mediated by DNA methylation and developing potential biomarkers for radiation-associated secondary diseases in BC.

Humans↗

Association between the gut microbiome and plasma metabolites linked to vocalization-based temperament in Merino sheep.

BACKGROUND: Temperament, as a determinant of behavioural and emotional responses, has a substantial adaptive value in different environments. This study aims to investigate the association between the gut microbiota and temperament plasticity, and clarify the potential metabolic mechanism that underpins that association by running a multi-omics study in sheep. METHODS: The TrackSheep research cohort was generated using 200 healthy juvenile Merino ewes, and the rumen microbiota, plasma metabolome, and temperament phenotype was measured. RESULTS: Rumen metagenomic analysis identified 25 microbial species and 16 MetaCyc pathways that explained 37.5% and 11.1%, respectively, of the variation in temperament as estimated using the vocal reactivity to stress. Among these, the &#x3b3;-aminobutyric acid (GABA) shunt and allantoin degradation pathways showed the strongest associations with vocal behaviour. Multi-omic integration linked these microbial pathways to plasma metabolites that are involved in neurotransmission, antioxidant defense, and energy metabolism, including acetyl-L-carnitine (ALCAR) and urocortisone, which partially mediated the effects of microbial pathways on vocalisations. Notably, functional genomic and mediation analyses indicated that the abundance of Cryptobacteroides sp902761655 was associated with the activity of GABA shunt pathway, where GABA co-occurred with succinate production, in turn correlating with reduced inhibitory effects of ALCAR on stress-susceptible temperament. Although plasma metabolite shifts observed immediately after behavioural tests reflected stress exposure, their associations with rumen microbiota highlight microbiome-metabolite interplay that could underly behavioural variation. CONCLUSIONS: Our study provides the first large-scale multi-omics evidence linking the rumen microbiome to a dimension of emotional reactivity in livestock, while underscoring the need for longitudinal and experimental validation to establish causal mechanisms. Video Abstract.

Animals↗