Corrigendum to "Reticulate evolution, introgression, and recent diversification in Epimedium sect. Macroceras" [Molecul. Phylog. Evol. 223 (2026) 108646].
Explore the source record for details and available documents.
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.
Paramecium bursaria maintains a stable endosymbiosis with green algae, yet the evolutionary consequences of this association remain unclear. Here, we screened the host genome for algal-derived horizontally transferred genes (HTGs) using a lineage-aware workflow designed to detect horizontal gene transfer (HGT) between two defined lineages. We identified 16 candidate HTGs, including four putative newly transferred genes and 12 homologous transferred genes, most of which were functionally associated with redox homeostasis and metabolism. Five HTGs showed symbiosis-dependent expression. RNAi knockdown of GH32s and SATs reduced host proliferation, total cell area, and motility, while GH32s knockdown also reduced endosymbiont load. Duplication patterns suggest that most transfers may have occurred after the P. bursaria lineage diverged from the sampled Paramecium species but before its lineage-specific whole-genome duplication (WGD). The HTGs also showed host-associated shifts in GC content and gene length, while representative HTGs retained conserved domains and functional motifs. Together, our results support algae-to-host HGT in P. bursaria and suggest that some transferred genes may contribute to metabolic integration during endosymbiosis.
The razor clam Sinonovacula constricta is an economically important cultured mollusk in China, but the molecular mechanism of its gonadal development and sexual differentiation remains unclear. This study integrated gonadal transcriptomic, proteomic, and metabolomic analysis to identify key sex-related molecules. Transcriptome analysis identified 2795 DELs and 6497 DEGs between sexes, including the sex-related genes Fem-1b, Fem-1c, GUCY1B2 and FAT4, as well as a regulatory network of 39 lncRNA-mRNA pairs involving Tektin-4, Ropporin-1, Histone H1, and FoxN4. Proteomic analysis revealed 3217 DEPs: Tektin family members, Ropporin-1 and Tssk proteins were upregulated in the testis, while histone H1 and FAT4 were upregulated in the ovary. Metabolomic analysis detected 409 DEMs, with uridine identified as a potential sex differential marker (upregulated in the ovary), and 23 gonadal development-related DEMs showed sex-specific upregulation. Integrative transcriptome-proteome analysis identified 1543 co-expressed DEGs/DEPs enriched in nucleosome assembly, oxidative phosphorylation, and carbon metabolism, including key sex-related genes AKAP14, Tektin/Tssk families, Histone H1, and FAT4. Transcriptome-metabolome integration identified 32 shared KEGG pathways (e.g., biosynthesis of unsaturated fatty acids, pyrimidine metabolism), while proteome-metabolome integration revealed 5 (positive ion) and 6 (negative ion) co-enriched pathways, with alanine, aspartate and glutamate metabolism and oxidative phosphorylation being functionally relevant to gonadal development. Collectively, these results reveal the molecular basis of gonadal development, highlight critical sex-related genes and steroid metabolic pathways, and provide valuable resources for future reproduction and breeding in S. constricta.
Erinacine A, a cyathane-type diterpenoid mainly from Hericium erinaceus mycelia, exhibits prominent neurotrophic and neuroprotective activities, making it a promising candidate for managing neurodegenerative diseases. However, its low abundance and unclear genetic regulatory mechanisms hinder its application as a nutraceutical. This study aimed to decipher its regulatory mechanisms and enhance production. Four exogenous inducers were screened, with salicylic acid (SA) and ergosterol (ERG) significantly increasing erinacine A content by 62.21% and 146.70% at 20 days, respectively. Transcriptome and WGCNA of inducer-treated sample identified darkorange and magenta modules associated with erinacine A biosynthesis, with the eri gene cluster enriched in the darkorange module and eriG and eriF as hub genes. Forward genetic analysis via QTL mapping of the HeD127 dikaryon population revealed significant phenotypic variation in erinacine A content (0.341-13.085 mg/g) and identified two loci (erA-1 and erA-2) explaining 18.63% of phenotypic variation. Integrating these forward and reverse genetic analyses revealed that salicylic acid and ergosterol synergistically regulate core carbon metabolic pathways to augment acetyl-CoA supply for the mevalonate pathway, suppressed competitive metabolism, enhanced diterpene skeleton construction and structural modification. These results deepen our understanding of the genetic and molecular basis governing accumulation of erinacine A, and facilitate its application in neuroprotective pharmaceuticals.
Acute cold stress can elicit distinct molecular responses even when bay scallop populations show similar phenotypic outcomes. We compared a seventh-generation fast-growing bay scallop line (BS) with a commercial control population (CC) during a 72-h acute cold exposure at -1 ± 0.3 °C. RNA-seq was used as the discovery layer, representative BS cold-responsive genes were evaluated by qRT-PCR, and paired LC-MS profiles provided a comparative metabolic layer. At baseline, 138 genes differed between BS and CC; after cold exposure, 134 of these baseline differences disappeared and 61 of 65 cold-state differences newly emerged. BS showed a larger transcriptomic response magnitude than CC, with 1129 cold-responsive genes compared with 28 genes in CC, and this ordering remained robust across multiple sensitivity analyses. Survival after 72 h was identical in BS and CC (83/90, 92.2% in each population). Biochemical responses were time-dependent and marker-specific: CAT, LZM, T-SOD and T-AOC showed population-by-time interactions, whereas GSH-Px and MDA did not, and the 72-h differences were not consistently favourable to BS. Metabolomic cold effects were strongly concordant between populations, and no feature showed a significant population-by-cold interaction. Features putatively assigned to arachidonic acid metabolism were enriched, but this provider-annotated pathway signal remains exploratory because authentic-standard confirmation was not performed. These findings indicate population-specific differences in molecular responsiveness but do not establish superior cold tolerance in BS.
In aquatic environments, natural and anthropogenic factors commonly reduce dissolved oxygen (DO) and trigger hypoxia, which threatens the health and survival of aquatic organisms. As an important economic fish species in China, pikeperch (Sander lucioperca) is extremely sensitive to hypoxia. However, there are relatively few reports on how hypoxia and reoxygenation affect its intestinal physiology and microbial community. Three treatment groups were set for pikeperch: normoxia (DO = 8.5 ± 0.5 mg/L), 48 h hypoxia (DO = 2.5 ± 0.1 mg/L), and reoxygenation (48 h hypoxia followed by 6 h reoxygenation at normal DO), to evaluate alterations in intestinal histopathology, tight junction gene expression, oxidative stress, hypoxia signaling molecules and intestinal microbiota composition. The results showed that hypoxia significantly decreased muscularis thickness by approximately 32.5% and reduced the expression of tight junction genes (Occludin, Claudin2, and ZO-2). Moreover, hypoxia significantly increased oxidative stress index levels (GSH-Px, CAT, and MDA), markedly upregulated the expression of Bax, Caspase3, and HIF-1α, while significantly downregulating the expression of Bcl-2, Egln1, and Egln2. Notably, reoxygenation elicited partial compensatory effects against these hypoxia-induced changes. 16S rRNA sequencing analysis revealed that hypoxic stress altered the intestinal microbial community composition of pikeperch and increased its diversity. In the hypoxia group, the abundance of the phylum Bacillota, along with the genera Halomonas and Acinetobacter, was significantly elevated, whereas in the reoxygenation group, the genus Lactobacillus increased approximately 180-fold. The results indicated that hypoxia caused intestinal oxidative damage, cell apoptosis, and intestinal microbiota dysbiosis in pikeperch, while short-term reoxygenation achieved partial recovery from these hypoxia-triggered intestinal injuries. The present research provides valuable references for in-depth exploration of the molecular mechanisms behind the response of pikeperch to acute hypoxia and reoxygenation stress, while also offering a novel perspective to understand the mechanism by which hypoxia impacts intestinal health in fish.
Mactra veneriformis is an economically important bivalve mollusc in China, but its aquaculture is frequently threatened by Vibrio infections, particularly Vibrio alginolyticus. To investigate the molecular immune response of M. veneriformis to V. alginolyticus, we performed RNA-seq analysis of hepatopancreatic tissues collected at 48 h post-infection, the peak mortality time point, with PBS-injected individuals used as controls. Infection with V. alginolyticus caused severe histopathological damage in the hepatopancreas and resulted in a cumulative mortality of 53.3% over 14 d, compared with 3.3% in the control group. Transcriptomic analysis identified 2623 differentially expressed genes (DEGs), including 1585 significantly up-regulated genes and 1038 down-regulated genes. KEGG enrichment analysis demonstrated that DEGs were significantly enriched in immune related and metabolism pathways, including the JAK-STAT signaling pathway, RIG-I-like receptor (RLR) signaling pathway, and cytochrome P450 (CYP450) signaling pathway. Collectively, these findings revealed candidate immune related genes (tlr3, tlr5, myd88, nfkb1, il-17d, and ifi44l), a putative TLR-MyD88-NF-κB signaling axis, and KEGG signaling pathways, including JAK-STAT, RLR and CYP450, that may be involved in the innate immune response of M. veneriformis to V. alginolyticus infection. These results provide a transcriptomic basis for understanding host-pathogen interactions in this species and highlight candidate genes and pathways for future functional validation and potential application in disease-resistance breeding.
Using a metal-organic framework (MOF) - supported Ru(bpy)32+ (MIL-88@Ru) composite luminescent material, this study innovatively adopted electropolymerization to fabricate a molecularly imprinted polymer-based electrochemiluminescent (MIP-ECL) sensor for enrofloxacin (ENR) detection in animal-derived foods. Systematic investigation of the ECL luminescence and ENR's quenching mechanism confirmed that the sensor integrates ECL's high sensitivity and MIP's high specificity, enabling rapid and accurate recognition of ENR. Experimental results show a good linear response in the range of 1 nmol/L-20 μmol/L (R2 = 0.99), a limit of detection (LOD) as low as 0.28 nmol/L, as well as excellent selectivity and stability. Recoveries of ENR in all investigated matrices ranged from 97.7% to 106.4%, confirming the reliability of the established method. This ECL-MIP coupling strategy provides a new technical approach and application references for the efficient detection of trace pollutants in food safety and environmental monitoring fields.
Parasitic mites Varroa destructor and Tropilaelaps mercedesae represent major threats to global honey bee (Apis mellifera) health and productivity, yet comparative molecular insights into host responses remain limited. To address this, we systematically compiled published studies (2015-2025) reporting genes associated with honey bee interactions with V. destructor (11 studies, 87 genes), T. mercedesae (4 studies, 35 genes), and hygienic behavior (6 studies, 44 genes). Gene identifiers were harmonized to the Amel_HAv3.1 genome assembly, yielding three non-redundant sets: 64 Varroa-associated, 34 Tropilaelaps-associated, and 44 hygienic behavior-associated genes. Venn analysis identified 10 overlapping genes (including A0A088A8D5, A0A088ADL8, ABAE_APIME, Def1, Def2, Gapdh, HYTA_APIME, Imd, LOC726783, and Vg), suggesting conserved defense mechanisms, while 41 and 24 genes were uniquely associated with Varroa and Tropilaelaps, respectively. Enrichment analyses revealed Varroa-responsive genes were enriched in immune processes, chitin catabolism, and signaling pathways (Toll/Imd, MAPK, Wnt). Tropilaelaps-associated genes were enriched for antibacterial defense and stress response, with Toll/Imd signaling as the sole significantly enriched pathway. Overlapping genes reinforced core innate immunity activation. Protein-protein interaction network centrality analysis identified key hub genes: Def1, HYTA_APIME, ABAE_APIME, PPO, Imd, PGRP-LC, Vg for Varroa; and ACPH1_APIME, MRJP1, Vg, LOC726783 for Tropilaelaps. Results demonstrate that, despite differences in mite biology, honey bees show a conserved immune response against both parasites, centered on antibacterial defense, humoral immunity, and activation of the Toll/Imd pathway. Although limited by the in-silico nature and research asymmetries reflecting Tropilaelaps' emergence, this curated resource establishes a comprehensive framework for elucidating shared and distinct molecular defense mechanisms. Ultimately, this approach prioritizes diagnostic markers and candidate genes for functional validation and breeding strategies to enhance colony resilience against mite‑driven disease globally.
Chinese Bahaba (Bahaba taipingensis) is a Class I protected marine fish endemic to China. Low larvae survival during artificial breeding severely hinder population recovery. To investigate the molecular mechanism of high mortality in larval fish, this study performed RNA-seq on liver from naturally deceased (ND) and mass-dead (MD) individuals, combined with least absolute shrinkage and selection operator (LASSO) regression and random forest (RF) algorithms to screen for core signature genes. A total of 873 differentially expressed genes (DEGs) were identified, including 112 upregulated and 761 downregulated genes. GO and KEGG enrichment analyses revealed significant enrichment in amino acid metabolism disorders, one‑carbon folate pool impairment, PPAR signaling abnormalities, ECM-receptor interaction, focal adhesion pathway, indicating widespread metabolic suppression accompanied by extracellular matrix remodeling and signaling disturbances in the livers of MD fish. MAD pre-filtering combined with dual machine learning algorithms yielded 18 robust core signature genes, among which SLC38A4, MMP1, FADD, FKBP5, and APOB were consistently identified as high-frequency core genes by both algorithms. SLC38A4 exhibited the highest importance score in the RF model and was significantly downregulated, making it the primary molecule distinguishing ND from MD phenotypes. ROC curve analysis showed that both models achieved an AUC of 1.000 (95% CI lower bound: 0.610), confirming the precise discriminatory ability of the core genes. GSEA further demonstrated significant enrichment of this core gene set in ND samples. This study provides the first systematic elucidation of the molecular mechanisms underlying liver dysfunction in low survival rate B. taipingensis, characterized by amino acid transport impairment, metabolic reprogramming, and structural remodeling, offering theoretical foundations for health assessment, early mortality risk warning, and artificial breeding conservation of this species.
Triple-negative breast cancer (TNBC) remains a highly aggressive and therapeutically challenging subtype, defined by the absence of oestrogen, progesterone, and HER2 expression. Tumour Protein 53 (TP53) mutations represent the most frequent genetic alteration, occurring in over 80% of cases and driving tumour initiation, progression, and therapeutic resistance. Mutant p53 proteins not only lose canonical tumour-suppressive functions but also often acquire gain-of-function (GOF) oncogenic properties that promote metastasis, genomic instability, and resistance to mechanisms like ferroptosis. This review examines the biological role of TP53 in TNBC pathogenesis and evaluates emerging pharmacological strategies aimed at targeting these vulnerabilities. Key approaches include the pharmacological reactivation of mutant p53 using small molecules such as APR-246, COTI-2, and the mutation-specific reactivator rezatapopt (PC14586), which has shown significant clinical tumour reduction in Y220C-mutant patients. Other strategies involve targeted protein degradation, the exploitation of synthetic lethal interactions (e.g., Chk1 or Aurora kinase B inhibition), and the use of natural products like cryptolepine or piperine derivatives. Recent clinical evidence further highlights the potential of combining epigenetic agents like decitabine with chemotherapy in TP53-mutant populations. Integrating TP53 mutation status into biomarker-driven treatment paradigms is a pivotal step toward achieving precision oncology and improving clinical outcomes for patients with TNBC.
Staphylococcus aureus (S. aureus) represents a major foodborne pathogen whose enterotoxin production poses significant challenges to food safety due to its high environmental resistance and limited efficacy of conventional sterilization. Since the expression of enterotoxins is predominantly governed by the agr quorum sensing system, targeting this regulatory pathway has become a strategic choice for virulence control. This study elucidated the mechanism by which phloretin, a potential quorum sensing inhibitor, interferes with the agr system to attenuate virulence. To achieve this, the recombinant AgrA LytTR domain was expressed and purified, and its interaction with phloretin was characterized using thermal shift assays (TSA), electrophoretic mobility shift assays (EMSA), and molecular dynamics (MD) simulations. The results showed that phloretin specifically binds to the AgrA LytTR domain, enhancing its thermal stability and disrupting AgrA LytTR-agr operon binding by reducing the free energy of interaction between them, without causing significant structural rearrangement. Mechanistic analysis indicated that phloretin sterically hinders key β-sheet turn residues (HIS169, ASN201, ARG233), thereby impairing DNA recognition, downregulating RNAIII transcription, and inhibiting agr signaling. In cooked beef, phloretin significantly inhibited the secretion of enterotoxins and α-hemolysin, while delaying lipid oxidation and protein degradation, and maintaining the meat texture. These findings suggested that phloretin is a multifunctional substance with anti-virulence, antioxidant, and preservative properties, demonstrating its potential as a natural food preservative.
Recumbent cycling provides an adjustable configuration in which backrest angle and crank-axis height can alter lower-limb loading. Backrest angle changes the orientation of the active musculature relative to gravity, whereas crank-axis height modifies the foot-pedal interface and the ankle lever arm. Whether such positioning is associated with differences in lower-limb electromyographic (EMG) activation across both factors together has not been examined. Using a randomized crossover design, fifty-nine healthy young men cycled at 100 W and 50 rpm under four conditions combining two backrest angles (30°, 75°) and two crank-axis heights (Rh, high, 77 cm; Rl, low, 45 cm). Surface EMG from the rectus femoris (RF), biceps femoris (BF), and gastrocnemius medialis (GAS) was analyzed for root mean square (RMS) amplitude, absolute onset sequence, and relative recruitment timing (%Sequence, derived using Telescan software and referenced to RF). RMS amplitude did not differ across conditions for any muscle (all p ≥ 0.356), and absolute onset sequence showed no condition effect (all p ≥ 0.197). For GAS %Sequence, the Greenhouse-Geisser-corrected omnibus test was significant (F(1.32, 76.56) = 28.45, p < 0.001), driven by a within-30° contrast in which %Sequence was lower under Rl than Rh (paired t(58) = 10.39, Holm-adjusted p < 0.0001; mean difference 22.0 percentage points, 95 % CI [17.8, 26.3]; dz = 1.35). Because %Sequence is referenced to RF, this exploratory shift should be interpreted cautiously. Lower-limb activation may be constrained by the closed-chain nature of pedaling, whereas distal recruitment timing may be sensitive to crank-axis height and the foot-pedal interface.
Rapid fluctuations in water temperature driven by global warming have become a major abiotic stressor affecting muscle function in teleost fish. This study examined the effects of long-term thermal conditions on heat tolerance in Clarias fuscus. Fish were maintained for 90 days at either a normal temperature group (NT, 26 °C) or a high-temperature group (HT, 34 °C). Subsequently, muscle histology, and transcriptomic profiles were observed following acute high-temperature exposure (34 °C) and after temperature recovery (26 °C). Histological analysis showed that fish from the NT under acute high-temperature stress exhibited severe muscle damage (atrophy, myofilament disruption, and myolysis), whereas fish from the HT displayed markedly reduced lesions. RNA-seq profiling revealed 5769 differentially expressed genes (DEGs) in the NT and 3292 DEGs in the HT following acute temperature challenges. Functional enrichment indicated that, in the HT, modulation of key cell cycle regulators (e.g., ccna, ccnb, cdk1, cdk2) contributed to alleviating muscle damage caused by temperature fluctuations. In the NT, genes associated with ribosome biogenesis (e.g., nop56, riok2, riok1) were up-regulated and then down-regulated during temperature fluctuation, whereas p53 in the cell cycle pathway showed the opposite expression pattern. These findings demonstrate long-term heat exposure reshapes molecular expression and regulatory mechanisms in the muscle of C. fuscus, thereby enhancing thermal tolerance and adaptability, and providing a theoretical basis for breeding heat-resistant, high-quality aquaculture strains.
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.
Underwater noise is a widespread environmental pollutant in marine ecosystems, yet the effects of low- to mid-frequency acoustic exposure on immune physiology and molecular responses in echinoderms remain unclear. In this study, the sea urchin Strongylocentrotus intermedius was exposed to continuous pure-tone acoustic stimulation at 80, 125, 250, 500, 750, and 1000 Hz for 3 h. Results showed that acoustic exposure significantly affected redox homeostasis, energy metabolism, and immune function in S. intermedius coelomocytes: the antioxidant system and glutathione redox balance were altered, as indicated by increased superoxide dismutase (SOD) and catalase (CAT) activities, elevated reduced glutathione (GSH) content, and a higher GSH/GSSG ratio; glycolysis-related enzyme activities were enhanced, with increased pyruvate kinase (PK) activity under 125-500 Hz exposure and elevated hexokinase (HK) activity at 250 Hz; and immune function was impaired, as shown by increased coelomocyte mortality, reduced phagocytic activity, and inhibited acid phosphatase (ACP) and alkaline phosphatase (AKP) activities, whereas respiratory burst activity showed no significant change. Among all treatments, 250 Hz induced the most pronounced physiological responses. Transcriptomic analysis of coelomocytes from the 250 Hz group identified 663 differentially expressed genes, including 537 upregulated and 126 downregulated genes, mainly enriched in pathways related to apoptosis, phagosome, lysosome, glutathione metabolism, arachidonic acid metabolism, and carbohydrate metabolism. These findings indicate that low- to mid-frequency acoustic exposure can act as a physiological and molecular stressor to S. intermedius coelomocytes by affecting redox homeostasis, enhancing energy metabolism, and suppressing immune effector processes, with 250 Hz showing the strongest effect under the present exposure conditions. This study provides experimental evidence for evaluating the potential biological effects of low- to mid-frequency acoustic exposure on benthic echinoderms.
Lynch syndrome (LS) carriers are at elevated cancer risk. Emerging evidence suggests that exercise may serve as a non-pharmacologic preventive strategy, yet the epigenetic and immunological mechanisms underlying its protective effects in this population remain unclear. Here, we perform integrative multi-omics profiling of DNA methylation, gene expression, and the T cell receptor (TCR) repertoire in LS carriers undergoing a 52-week aerobic cycling intervention. We identify compartment-specific DNA methylation changes, including innate immune activation in cfDNA and oncogenic pathway repression in tissue. Integrative transcriptomic analysis highlights ISL1 as a key exercise-repressed, epigenetically regulated gene, and identifies FLCN as a colorectal cancer (CRC)-associated methylation target. TCR analysis reveals an exercise-associated increase in systemic repertoire diversity and tissue-specific clonal convergence, thus suggesting antigen-driven recruitment. Collectively, these findings uncover epigenetic and immune remodeling as potential mechanisms of exercise-mediated protection in LS.
Yaks (Bos grunniens) are crucial for the livelihoods of pastoral communities in cold regions, where feed scarcity during the cold season poses challenges to their health and productivity, underscoring the necessity of understanding how dietary management influences the gut microbiota. In this study, 24 yak steers matched for body weight and health status were randomly allocated to four groups: natural grazing or indoor feeding with roughage-to-concentrate ratios of 50:50, 70:30, or 90:10. Fecal samples were collected for 16S rDNA sequencing and subsequent functional prediction of the microbiota. The results showed that Firmicutes and Bacteroidetes were the dominant phyla across all groups, and UCG-005 and Rikenellaceae_RC9_gut_group were the predominant genera. Concentrate supplementation during the cold season significantly enhanced microbial richness and diversity, with the 70:30 ratio exerting the most pronounced beneficial effects on microbiota structure and key taxa enrichment. These findings highlight the critical role of dietary management in shaping the yak gut microbiota during cold seasons and suggest that the 70:30 ratio optimally improves microbial community structure, thereby promoting yak health and productivity under harsh climatic conditions. Future research should explore the long-term implications of such dietary strategies.