ObjectiveHeadspace solid-phase microextraction-gas chromatography-mass spectrometry（HS-SPME-GC-MS） and GC-triple quadrupole MS（GC-QqQ-MS） in combination with non-targeted and targeted metabolomics were employed to systematically analyze the chemical composition differences of Xihuangwan prepared with natural musk and artificial musk， and establish an identification system for them. MethodsThe volatile components of 9 batches of Xihuangwan samples from 8 manufacturers were analyzed by HS-SPME-GC-MS non-targeted metabolomics， and identified by comparing their MS data with the National Institute of Standards and Technology（NIST） spectral library. Orthogonal partial least squares-discriminant analysis（OPLS-DA） was used to identify differential volatile components of Xihuangwan prepared with natural musk and artificial musk. Additionally， GC-QqQ-MS targeted metabolomics was applied to quantify the levels of α-pinene， β-elemene， muscone， dehydroepiandrosterone， bornyl acetate， and octyl acetate in 27 batches of samples from 9 manufacturers. Cluster analysis， principal component analysis（PCA）， and partial least squares-discriminant analysis（PLS-DA） were conducted to further explore the differences in volatile components between Xihuangwan samples prepared with natural musk and artificial musk. ResultsNon-targeted metabolomics identified 291 volatile compounds in Xihuangwan， including alkanes， esters， alkanes， alcohols， ketones， naphthalenes and others. OPLS-DA analysis revealed distinct separation between Xihuangwan samples containing artificial musk（A1， C1， D1， E1， F1， G1， I1） and those containing natural musk（H1， H3）. A total of 30 differential metabolites were identified. The relative contents of these 30 differential metabolites were visualized using a radar chart， revealing significant differences in the levels of octanol， borneol acetate and muscone. Cluster analysis and PCA results from targeted metabolomics indicated that Xihuangwan could be classified into two distinct groups：one composed of natural musk（H1， H3） and the other of artificial musk， sample H2. PLS-DA identified muscone， octyl acetate， and dehydroepiandrosterone as key differential volatile components. Although no significant difference was observed in the content of octyl acetate between the two groups， statistically significant differences were found for muscone and dehydroepiandrosterone（P<0.05）. ConclusionMuscone and dehydroepiandrosterone can be used for the differentiation of Xihuangwan samples containing natural musk from those containing artificial musk. This study systematically and comprehensively analyzed the differences in the types and contents of major volatile components in Xihuangwan prepared with natural musk and artificial musk， providing a scientific basis for quality evaluation and control of Xihuangwan.

ObjectiveHeadspace solid-phase microextraction-gas chromatography-mass spectrometry（HS-SPME-GC-MS） and GC-triple quadrupole MS（GC-QqQ-MS） in combination with non-targeted and targeted metabolomics were employed to systematically analyze the chemical composition differences of Xihuangwan prepared with natural musk and artificial musk， and establish an identification system for them. MethodsThe volatile components of 9 batches of Xihuangwan samples from 8 manufacturers were analyzed by HS-SPME-GC-MS non-targeted metabolomics， and identified by comparing their MS data with the National Institute of Standards and Technology（NIST） spectral library. Orthogonal partial least squares-discriminant analysis（OPLS-DA） was used to identify differential volatile components of Xihuangwan prepared with natural musk and artificial musk. Additionally， GC-QqQ-MS targeted metabolomics was applied to quantify the levels of α-pinene， β-elemene， muscone， dehydroepiandrosterone， bornyl acetate， and octyl acetate in 27 batches of samples from 9 manufacturers. Cluster analysis， principal component analysis（PCA）， and partial least squares-discriminant analysis（PLS-DA） were conducted to further explore the differences in volatile components between Xihuangwan samples prepared with natural musk and artificial musk. ResultsNon-targeted metabolomics identified 291 volatile compounds in Xihuangwan， including alkanes， esters， alkanes， alcohols， ketones， naphthalenes and others. OPLS-DA analysis revealed distinct separation between Xihuangwan samples containing artificial musk（A1， C1， D1， E1， F1， G1， I1） and those containing natural musk（H1， H3）. A total of 30 differential metabolites were identified. The relative contents of these 30 differential metabolites were visualized using a radar chart， revealing significant differences in the levels of octanol， borneol acetate and muscone. Cluster analysis and PCA results from targeted metabolomics indicated that Xihuangwan could be classified into two distinct groups：one composed of natural musk（H1， H3） and the other of artificial musk， sample H2. PLS-DA identified muscone， octyl acetate， and dehydroepiandrosterone as key differential volatile components. Although no significant difference was observed in the content of octyl acetate between the two groups， statistically significant differences were found for muscone and dehydroepiandrosterone（P<0.05）. ConclusionMuscone and dehydroepiandrosterone can be used for the differentiation of Xihuangwan samples containing natural musk from those containing artificial musk. This study systematically and comprehensively analyzed the differences in the types and contents of major volatile components in Xihuangwan prepared with natural musk and artificial musk， providing a scientific basis for quality evaluation and control of Xihuangwan.

ObjectiveTo explore the complex network relationships among pain, kinesiophobia, social participation and knee function in patients after total knee arthroplasty (TKA), and to analyze the moderating effects of different socio-structural factors on the rehabilitation network from an ethical equity perspective. MethodsA convenience sampling method was used to select 291 patients who underwent TKA in Qilu Hospital of Shandong University from May to July, 2023. Pain was assessed using Numerical Rating Scale, kinesiophobia with Chinese short version of the Tampa Scale for Kinesiophobia, social participation with Impact on Participation and Autonomy Questionnaire, and knee function with Hospital for Special Surgery Knee Score. A partial correlation network among pain, kinesiophobia, social participation and knee function was constructed using Graphical Least Absolute Shrinkage and Selection Operator. Key variables were identified through node centrality and bridge centrality analysis. Network Comparison Tests (NCT) were used to analyze network differences among subgroups based on different socio-structural characteristics. ResultsIn the network model, the nodes with the highest strength centrality were indoor participation, activity behavior and activity pain. Bridge centrality analysis indicated that activity pain, knee function, indoor participation and activity cognition were key bridge nodes. NCT revealed no significant differences in overall network structure or global strength among subgroups based on residence, education level or payment method (P > 0.05). However, significant differences in edge weights were found for specific edges such as activity cognition-activity behavior and knee function-indoor participation (P < 0.05). ConclusionThere is a network of interactions among pain, kinesiophobia, social participation and knee function in patients after TKA, with nodes such as indoor participation and activity pain playing key roles in the rehabilitation process. Although the overall rehabilitation network is similar across different socio-structural groups, variations exist in specific relational pathways among patients from rural areas, those with lower education levels, and those with out-of-pocket payment. This suggests that clinical rehabilitation interventions should focus on these core nodes and implement targeted support strategies for socio-structurally disadvantaged groups to promote rehabilitation equity.

ObjectiveTo explore the complex network relationships among pain, kinesiophobia, social participation and knee function in patients after total knee arthroplasty (TKA), and to analyze the moderating effects of different socio-structural factors on the rehabilitation network from an ethical equity perspective. MethodsA convenience sampling method was used to select 291 patients who underwent TKA in Qilu Hospital of Shandong University from May to July, 2023. Pain was assessed using Numerical Rating Scale, kinesiophobia with Chinese short version of the Tampa Scale for Kinesiophobia, social participation with Impact on Participation and Autonomy Questionnaire, and knee function with Hospital for Special Surgery Knee Score. A partial correlation network among pain, kinesiophobia, social participation and knee function was constructed using Graphical Least Absolute Shrinkage and Selection Operator. Key variables were identified through node centrality and bridge centrality analysis. Network Comparison Tests (NCT) were used to analyze network differences among subgroups based on different socio-structural characteristics. ResultsIn the network model, the nodes with the highest strength centrality were indoor participation, activity behavior and activity pain. Bridge centrality analysis indicated that activity pain, knee function, indoor participation and activity cognition were key bridge nodes. NCT revealed no significant differences in overall network structure or global strength among subgroups based on residence, education level or payment method (P > 0.05). However, significant differences in edge weights were found for specific edges such as activity cognition-activity behavior and knee function-indoor participation (P < 0.05). ConclusionThere is a network of interactions among pain, kinesiophobia, social participation and knee function in patients after TKA, with nodes such as indoor participation and activity pain playing key roles in the rehabilitation process. Although the overall rehabilitation network is similar across different socio-structural groups, variations exist in specific relational pathways among patients from rural areas, those with lower education levels, and those with out-of-pocket payment. This suggests that clinical rehabilitation interventions should focus on these core nodes and implement targeted support strategies for socio-structurally disadvantaged groups to promote rehabilitation equity.

OBJECTIVE To explore the ameliorative effect and potential mechanism of vitexin on inflammation in ulcerative colitis （UC） mice. METHODS The UC mice model was established by continuous administration of 3% dextran sulfate sodium solution for 5 days. Mice with successful modeling were randomly divided into UC group， vitexin low- and high-dose groups （vitexin-L and vitexin-H groups， 40， 80 mg/kg）， mesalazine group （400 mg/kg）， and vitexin-H+recombinant Jagged canonical Notch ligand 1 （rJagged-1） group （vitexin-H+rJagged-1 group， 80 mg/kg vitexin+1 mg/kg rJagged-1）， with 12 mice in each group. Another 12 normal mice were used as the control （CK） group. Mice in each group were administered the corresponding drugs or the corresponding drugs and normal saline by gavage and intraperitoneal injection once daily for 7 consecutive days. General conditions were observed during the experiment. At 24 h after the last administration， the disease activity index （DAI） score was evaluated. Colonic histopathological morphology was observed and scored. Macrophage polarization levels in the spleen and colon tissues were measured. The protein expressions of interleukin-6 （IL-6）， IL-10， tumor necrosis factor-α （TNF-α）， transforming growth factor-β 1 （TGF-β 1 ）， Jagged-1， Notch1 and Notch intracellular domain （NICD） in colonic tissues were determined. RESULTS Compared with the UC group， the symptoms （reduced food and water intake， dull fur， etc.） and pathological changes （epithelial cell shedding， inflammatory cell infiltration， etc.） were significantly improved in the vitexin-L， vitexin-H and mesalazine groups. DAI scores， colonic histopathological scores， M1 macrophage contents in spleen tissue， M1/M2 macrophage ratios， M1 macrophage proportions in colon tissue， and protein expressions of IL-6， TNF-α， Jagged-1， Notch1 and NICD in colon tissue were significantly decreased （ P ＜0.05）. Meanwhile， the M2 macrophage contents in spleen tissue， M2 macrophage proportions in colon tissue， and protein expressions of IL-10 and TGF-β 1 in colon tissue were significantly increased （ P ＜0.05）. Moreover， the improvement effects in the vitexin-H and mesalazine groups were significantly superior to those in the vitexin-L group （ P ＜0.05）. Compared with the vitexin-H group， the above symptoms and pathological changes were aggravated， and all quantitative indicators were significantly reversed in the vitexin-H+rJagged-1 group （ P ＜0.05）. CONCLUSIONS Vitexin can ameliorate the inflammation of UC mice， which is associated with its inhibition of the Jagged-1/Notch1 pathway and regulation of macrophage polarization （inhibition of M1-type polarization and promotion of M2-type polarization）.

ObjectiveTo explore the complex network relationships among pain, kinesiophobia, social participation and knee function in patients after total knee arthroplasty (TKA), and to analyze the moderating effects of different socio-structural factors on the rehabilitation network from an ethical equity perspective. MethodsA convenience sampling method was used to select 291 patients who underwent TKA in Qilu Hospital of Shandong University from May to July, 2023. Pain was assessed using Numerical Rating Scale, kinesiophobia with Chinese short version of the Tampa Scale for Kinesiophobia, social participation with Impact on Participation and Autonomy Questionnaire, and knee function with Hospital for Special Surgery Knee Score. A partial correlation network among pain, kinesiophobia, social participation and knee function was constructed using Graphical Least Absolute Shrinkage and Selection Operator. Key variables were identified through node centrality and bridge centrality analysis. Network Comparison Tests (NCT) were used to analyze network differences among subgroups based on different socio-structural characteristics. ResultsIn the network model, the nodes with the highest strength centrality were indoor participation, activity behavior and activity pain. Bridge centrality analysis indicated that activity pain, knee function, indoor participation and activity cognition were key bridge nodes. NCT revealed no significant differences in overall network structure or global strength among subgroups based on residence, education level or payment method (P > 0.05). However, significant differences in edge weights were found for specific edges such as activity cognition-activity behavior and knee function-indoor participation (P < 0.05). ConclusionThere is a network of interactions among pain, kinesiophobia, social participation and knee function in patients after TKA, with nodes such as indoor participation and activity pain playing key roles in the rehabilitation process. Although the overall rehabilitation network is similar across different socio-structural groups, variations exist in specific relational pathways among patients from rural areas, those with lower education levels, and those with out-of-pocket payment. This suggests that clinical rehabilitation interventions should focus on these core nodes and implement targeted support strategies for socio-structurally disadvantaged groups to promote rehabilitation equity.

Objective To construct an index system for assessment of schistosomiasis transmission risk following natural disasters such as rainstorms, floods, earthquakes, mudslides, and landslides, so as to provide insights into rapid identification of schistosomiasis transmission risk post-disasters and formulation of targeted schistosomiasis control strategies. Methods An initial framework for the index system for assessment of schistosomiasis transmission risk following natural disasters was drafted through literature review, brainstorming, and focus group discussions. Two rounds of expert correspondence consultations were conducted using the Delphi method to refine and finalize the system, and the degrees of expert activeness, authority and endorse ment, and consensus were evaluated. In addition, the weights of each index were calculated using the analytic hierarchy process. Results A total of 18 experts participated in the consultation. The expert positive coefficients were 100.00% and 94.44% for two rounds of consultations, with authority coefficients of 0.92 and 0.94, respectively. The coefficients of coordination on the index importance, rationality and operability were 0.209, 0.185, 0.222 and 0.407, 0.214, 0.257 for two rounds of consultations, respectively, and all consistency tests were statistically significant (χ² = 246.771 to 505.278, all P values < 0.001). Following two rounds of expert consultations, an index system consisting of 6 first-level indicators, 15 second-level indicators, and 49 third-level indicators was ultimately constructed. In terms of first-level indicators, “disaster situation”, “previous epidemics”, “healthcare guarantee”, “response capacity” and “emergency recovery” had the highest weights, each at 18.18%. Regarding second-level indicators, “Schistosoma japonicum infections in animals”, “S. japonicum infections in snails” and “medical treatment” had the highest weights, each at 7.35%. In terms of third-level indicators, ten items had the highest weights, including “identification of schistosomiasis cases”, “detection of S. japonicum infections in wild feces”, “detection of S. japonicum infections in snails”, “reserves of schistosomiasis diagnostic/testing reagents and consumables”, “reserves of chemotherapy agents for human and animal schistosomiasis”, “reserves of cercariacides”, “periodical surveillance on schistosomiasis”, “identification of schistosomiasis transmission risk and timely response”, “normal provision of diagnosis and treatment services” and “post-disaster schistosomiasis surveillance”, each at 2.40%. Conclusion A scientific, systematic, and practical index system has been constructed for assessment of schistosomiasis transmission risk following natural disasters, which may provide insights into rapid post-disaster identification of schistosomiasis transmission risk, formulation of targeted schistosomiasis control strategies and optimization of resource allocation.

ObjectiveThe malaria parasites remodel the host erythrocyte structure by exporting parasite proteins that interact with the membrane skeleton proteins of red blood cells (RBCs), facilitating their intracellular survival and pathogenicity. Skeleton-binding protein 1 (SBP1) is a conserved exported protein across Plasmodium species. In Plasmodium falciparum, SBP1 has been reported to interact with erythrocyte membrane skeleton proteins 4.1R and spectrin, while its contribution to erythrocyte remodeling and parasite virulence in Plasmodium berghei (Pb) remains unclear. This study aims to determine whether PbSBP1 associates with the host cytoskeletal protein 4.1R and to investigate its role in the remodeling of host RBCs and the pathogenicity of Plasmodium berghei. MethodsIn Plasmodium berghei, the relationship between PbSBP1 and the erythrocyte cytoskeletal protein 4.1R was examined using co-immunoprecipitation. A Pbsbp1 gene knockout mutant of Plasmodium berghei (Pbsbp1∆) was generated based on the principle of double crossover homologous recombination. The deformability of erythrocytes infected with Pbsbp1∆ parasites was assessed using microfluidic methods. Microchannels with an array of cylindrical pillars were used to detect modifications in infected RBC deformability. The infected RBCs were squashed between the rows and recovered between the columns and the transit velocity (μm/s) of infected RBCs travelling through the microchannel was recorded. The component of the erythrocyte membrane skeleton junctional complex, tropomodulin (TMOD), was fluorescently labeled, and the cytoskeletal network of infected erythrocytes was imaged using super-resolution stochastic optical reconstruction microscopy (STORM) to analyze ultrastructural changes in the cytoskeleton of wild-type (WT) and Pbsbp1∆-infected erythrocytes. Actin-based junctional complexes were displayed as individual clusters by the labeled TMOD in the STORM images, and the cluster densities and distances between adjacent clusters of infected RBCs were calculated. Additionally, rodent malaria models (BALB/c mice) and experimental cerebral malaria models (C57BL/6 mice) were employed to monitor the growth of Pbsbp1∆ and WT parasites during the intraerythrocytic stage and their capacity to induce cerebral malaria in mice. ResultsPbSBP1 may participate in the remodeling of infected erythrocytes through direct or indirect interaction with the erythrocyte cytoskeletal protein 4.1R. Microfluidic assays revealed that the deformability of erythrocytes infected with Pbsbp1∆ parasites was significantly enhanced compared to those infected with WT parasites. STORM imaging further demonstrated that the ultrastructure of the erythrocyte cytoskeleton in Pbsbp1∆-infected cells was altered relative to that in WT-infected erythrocytes. The distances between nearest neighbors of clusters had a tendency to increase while the cluster densities were decreased in Pbsbp1∆-infected RBCs compared to WT-infected RBCs. Subsequent phenotypic analysis indicated that the growth rate of Pbsbp1∆ parasites during the intraerythrocytic stage was significantly slower than that of WT parasites, and their ability to induce cerebral malaria in mice was also attenuated. These findings suggest that PbSBP1 is involved in the remodeling of the erythrocyte membrane skeleton, likely through its direct or indirect interaction with protein 4.1R, thereby regulating the deformability of infected erythrocytes and influencing the pathogenicity of the blood-stage parasites. ConclusionThis study establishes a role for PbSBP1 in host erythrocyte remodeling and parasite virulence, providing new research strategies for the prevention and treatment of malaria.

ObjectiveThe exacerbating trend of global population aging poses profound socioeconomic and public health challenges, making the comprehensive elucidation of biological aging mechanisms and the discovery of effective anti-aging interventions an urgent priority in the life sciences. Based on our previous serum metabolomics findings that dimethylglycine, an intermediate metabolite of amino acid metabolism naturally present in the human body, was significantly enriched in the serum of longevity families, this study aimed to systematically investigate the anti-aging effects of dimethylglycine both in living organisms and in controlled laboratory environments, and to preliminarily elucidate its underlying molecular mechanisms. While existing literature indicates that dimethylglycine possesses antioxidant and immunomodulatory properties, its direct anti-aging efficacy and the specific molecular pathways through which it operates remain largely unexplored. MethodsTo comprehensively evaluate the anti-aging properties of dimethylglycine, we utilized replicative senescent human embryonic lung fibroblasts, specifically the WI-38 cell line, as an experimental model in a controlled laboratory environment. Cell viability and safety were thoroughly assessed using Cell Counting Kit-8 and lactate dehydrogenase release assays across various concentrations of dimethylglycine. The impact of dimethylglycine on cellular senescence phenotypes, oxidative stress, and proliferative capacity was evaluated via senescence-associated beta-galactosidase staining, reactive oxygen species fluorescence detection, and 5-ethynyl-2'-deoxyuridine incorporation assays. Furthermore, the molecular alterations of senescence-associated secretory phenotype factors and core senescence signaling pathways were quantified using quantitative reverse transcription polymerase chain reaction for the messenger RNA levels of interleukin-6, interleukin-8, p21, and matrix metalloproteinase-1, and enzyme-linked immunosorbent assay for the measurement of p16 and p21 protein expression levels. For the living organism model, the wild-type nematode Caenorhabditis elegans was used to evaluate systemic physiological effects. We conducted a comprehensive lifespan analysis at 20°C, heat stress resistance survival assays at 35℃, senescence-associated beta-galactosidase staining, lipofuscin accumulation tracking, intracellular reactive oxygen species measurement, and Oil Red O staining to ascertain systemic lipid accumulation. Additionally, network pharmacology bioinformatics tools, including PharmMapper and STRING databases, and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were utilized to predict target pathways, alongside highly detailed molecular docking simulations utilizing SwissDock and Protein-Ligand Interaction Profiler to examine interactions with the cytochrome P450 family 2 subfamily C member 9 protein. ResultsThe experimental outcomes robustly demonstrate the potent anti-aging capabilities of dimethylglycine. At the cellular level, toxicity analyses firmly confirmed that dimethylglycine is highly safe; continuous treatment with 50 mol/L and 70 mol/L of dimethylglycine for 5 d did not induce any cellular membrane damage or cytotoxicity, but rather actively promoted cellular proliferation. Utilizing the optimal standardized concentration of 50 mol/L, dimethylglycine treatment significantly ameliorated senescent phenotypic markers in human embryonic lung fibroblasts, which was evidenced by a drastic and highly significant reduction in the senescence-associated beta-galactosidase positive cell percentage (P<0.000 1) and intracellular reactive oxygen species levels (P<0.000 1), alongside a marked increase in the 5-ethynyl-2'-deoxyuridine-positive proliferation rate (P=0.003 5). On a molecular expression scale, dimethylglycine significantly downregulated the messenger RNA expression of multiple core senescence-associated secretory phenotype inflammatory factors, including interleukin-6, interleukin-8, p21, and matrix metalloproteinase-1. Concurrently, it effectively suppressed the protein expression of critical cell cycle arrest markers, diminishing p16 protein levels by 57.3% (P=0.000 4) and p21 protein levels by 27.2% (P=0.000 7). In the nematode Caenorhabditis elegans animal model, dimethylglycine significantly extended the mean lifespan from 20.402 d to an impressive 23.066 d (P<0.000 1) and notably enhanced overall survival rates under severe heat stress environmental conditions (P=0.017). Furthermore, systemic dimethylglycine intervention significantly mitigated age-related physiological decline by decreasing bodily lipofuscin accumulation (P<0.000 1), significantly reducing senescence-associated beta-galactosidase activity, lowering systemic reactive oxygen species fluorescence (P=0.008), and effectively alleviating overall fat accumulation (P<0.000 1). Mechanistically, extensive network pharmacology and Kyoto Encyclopedia of Genes and Genomes analyses strongly revealed that the potential targets of dimethylglycine are significantly enriched in fundamental drug metabolism and oxidative stress response pathways. Precision molecular docking simulations conclusively demonstrated that dimethylglycine forms highly stable structural interactions with the cytochrome P450 family 2 subfamily C member 9 protein, specifically highlighting the definitive formation of 5 stable hydrogen bonds involving serine 365, leucine 366, and serine 429 residues, as well as two critical salt bridge formations with arginine 97 and histidine 368 residues. It is additionally predicted to interact favorably with glutathione S-transferase family proteins. ConclusionDimethylglycine exhibits a profoundly significant and multifaceted anti-aging activity at both the cellular and entire living animal levels. By powerfully alleviating oxidative stress, heavily suppressing the core p16 and p21-dependent cellular senescence signaling pathways, and substantially mitigating the detrimental senescence-associated secretory phenotype, dimethylglycine effectively delays fundamental cellular senescence processes and drastically extends whole-organism lifespan. The biological mechanisms driving these robust protective effects are highly likely closely associated with its direct stable interactions with crucial metabolic and detoxifying enzyme systems, such as cytochrome P450 family 2 subfamily C member 9 and glutathione S-transferase family proteins, thereby systemically improving metabolic dysregulation and restoring critical redox homeostasis. This comprehensive study provides highly solid experimental evidence supporting dimethylglycine as a highly potent and safe potential anti-aging intervention agent, while simultaneously offering a clear molecular mechanistic explanation for the previously documented high abundance of dimethylglycine observed within exceptionally long-lived human populations.

ObjectiveDonkey hide is the sole legally designated raw material for the preparation of the traditional Chinese medicine Ejiao. The quality stability of donkey hide during preservation directly determines the efficacy and safety of Ejiao. This study focuses on the dynamic succession of microbial communities during the preservation of donkey hides from different origins, aiming to clarify the correlation between microbial biodiversity difference and the degradation profiles of hide collagen and critical biochemical components, thereby providing a theoretical foundation for developing targeted preservation strategies based on microbial regulation. MethodsDonkey hides originating from four different regions were subjected to an accelerated microbial aging assay to simulate the spoilage process. The microbial community succession was analyzed using high-throughput sequencing. Microstructure changes and pore structure characteristics were assessed by scanning electron microscopy and mercury intrusion porosimetry, respectively. Additionally, the content of major components, including lipids, proteins, and sugars were determined by biochemical methods. ResultsAfter 96 h of aging, the collagen fiber structure in Africa donkey hides (ADH) exhibited significant degradation and collapse, followed by Xinjiang donkey hides (XDH). Instead, the microstructure of Dong’e black donkey hides (DDH) and Peru donkey hides (PDH) remained relatively intact. The porosities of DDH, XDH, PDH, and ADH increased from 27.9%, 15.7%, 30.3%, and 46.2% to 36.5%, 52.6%, 42.8%, and 57.7%, respectively, during the aging process, which suggested that the originally compact fiber structure was disrupted by microbial aging. Fourier transform infrared spectrometer analysis revealed the amide bands in XDH exhibited relatively weak intensity, and no collagen amide I band was observed in ADH. Meanwhile, the lipid and protein contents decreased in all four types of donkey hides, indicating that these components served as the primary nutrient sources for the growth of microorganism. Notably, the most severe collagen degradation was observed in XDH and ADH. A substantial increase was detected in the total soluble sugar in PDH aging solution and hydroxyproline in the ADH aging solution, respectively. These results indicated that donkey hides exhibit distinct patterns of structural degradation and nutrient utilization. Furthermore, the viable cells number of donkey hides increased sharply after 48 h of aging. Metagenomic analysis revealed that the relative abundance of Euryarchaeota in ADH, PDH and XDH declining from initial 93.19%, 97.73% and 30.08% to 0.79%, 1.43% and 0.02% after 96 h, respectively. Conversely, a significantly increase was observed in the abundance of Bacillota, with a marked increase in ADH, peaking at 92.75%. Additionally, the abundance of Pseudomonadota in PDH increased from 0.10% to 87.84%, suggesting that Bacillota and Pseudomonadota may be key factors exacerbating donkey hide spoilage. Unlike the other three types of donkey hides, the dominant bacterial phylum in DDH shifted from Pseudomonadota to Bacteroidota, characterized by a substantial abundance increase of Bacteroidota from 0.13% to 44.22%. ConclusionRegional variation in origin significantly influence the microbial aging of donkey hides, leading to distinct patterns of structural deterioration and differential nutrient utilization. Therefore, implementing origin-specific preservation strategies, through the precisely controlling environmental factors to suppress harmful phyla such as Bacillota and Pseudomonadota, is crucial for enhancing the storage quality of donkey hides.