ObjectiveTo investigate the effects of targeted silencing of Runt-related Transcription Factor 3 （RUNX3） on the proliferation and migration of Mouse Hepatic Stellate Cells （HSCs）， as well as subsequent collagen deposition. MethodsMouse hepatic stellate cell line （JS-1） was selected and then morphologically observed and identified under a microscope. After the cells had fully adhered， they were treated with 5 ng/mL of transforming growth factor beta 1 （TGF-β₁） for 24 hours to induce hepatic stellate cell activation. Furthermore， a RUNX3 silencing model was established using RUNX3 lentiviral infection. The experiment was divided into four groups： Control group， TGF-β₁ group， TGF-β₁+siRNA-NC group， and TGF-β₁+siRNA-RUNX3 group. Protein expression changes of RUNX3， alpha-smooth muscle actin （α-SMA）， and Alpha 1 type I collagen （Collagen I） were detected using Western blot method. Cellular immunofluorescence assays were employed to investigate the deposition changes of α-SMA and RUNX3 in hepatic stellate cells. RT-qPCR was utilized to examine the mRNA expression changes of RUNX3， α-SMA， and Collagen I. The proliferative capacity of hepatic stellate cells was assessed using Edu staining. The migratory ability of hepatic stellate cells was evaluated through wound healing assays and Transwell migration experiments. ResultsCompared with Control group， a significant elevation in RUNX3 was observed in the TGF-β₁-induced activated HSCs （P<0.01）. Meanwhile， the protein and mRNA levels of fibrosis-related markers and α-SMA and Collagen I were significantly upregulated （P<0.001）. Additionally， the proliferation and migration capabilities of HSCs were significantly enhanced （P<0.001）. In contrast， when compared to TGF-β₁+siRNA-NC group， TGF-β₁+siRNA-RUNX3 group exhibited a notable decrease in RUNX3 and other related indicators， such as the protein and mRNA levels of α-SMA and Collagen I （P<0.05）. Concurrently， the proliferation and migration capabilities of HSCs were significantly inhibited in TGF-β₁+siRNA-RUNX3 group （P<0.01）. ConclusionSilencing RUNX3 can inhibit the deposition of collagen and the proliferation and migration of hepatic stellate cells. Conversely， RUNX3 promotes the proliferation and migration capabilities of HSCs， thereby facilitating the activation of HSC.

ObjectiveCleft palate (CP) is a common congenital deformity often associated with velopharyngeal insufficiency (VPI), which disrupts the physiological coupling between respiration and speech. Conventional clinical assessments, such as nasometry and spirometry, provide limited static data and fail to visualize the dynamic spatiotemporal distribution of lung ventilation during phonation. This study introduces spatiotemporal electrical impedance tomography (ST-EIT) to evaluate speech-respiratory functional features in CP patients compared to normal controls (NC). The aim is to characterize multi-domain respiratory patterns and to validate an interpretable machine learning framework for providing objective, quantitative evidence for clinical assessment. MethodsSeventy-five participants were enrolled in this study, comprising 37 patients with surgically repaired CP and 38 healthy volunteers matched for age, gender, and body mass index (BMI). All subjects performed standardized sustained phonation tasks while undergoing synchronous monitoring with a 16-electrode EIT system and a pneumotachograph. A comprehensive feature engineering pipeline was developed to extract physiological parameters across 3 complementary domains. (1) Temporal domain: including inspiratory/expiratory phase duration (tPhase), time constants (Tau), and inspiratory-to-expiratory time ratios (TI/TE); (2) airflow domain: comprising mean flow, peak flow, and instantaneous flow at 25%, 50%, and 75% of tidal volume; and (3) spatial domain: quantifying global and regional tidal impedance variation (TIV), global inhomogeneity (GI), and center of ventilation (CoV). Extreme Gradient Boosting (XGBoost) classifiers were trained using 5 distinct data sources (Spirometry, Nasometry, Inspiratory-EIT, Expiratory-EIT, and fused ST-EIT). Model performance was rigorously evaluated via stratified 5-fold cross-validation, and Shapley additive explanations (SHAP) were employed to quantify global and local feature contributions. ResultsThe CP group exhibited a distinct respiratory phenotype compared to controls. In the temporal domain, CP patients showed significantly shorter inspiratory (1.60 s vs.1.85 s, P<0.001) and expiratory phase durations (2.45 s vs. 3.95 s, P<0.001), indicating a rapid, shallow breathing rhythm. In the airflow domain, while inspiratory flows were comparable, the CP group demonstrated significantly elevated mean and peak flows during the expiratory phase (P<0.001), reflecting compensatory respiratory effort. Spatially, CP patients presented significant ventilation redistribution, characterized by higher regional TIV in the right-anterior (ROI1) and left-posterior (ROI4) quadrants, but lower TIV in the left-anterior (ROI2) quadrant. In terms of diagnostic accuracy, the multi-modal ST-EIT model achieved the highest performance (AUC: 0.915±0.012, Accuracy: 0.843±0.019, F1-score: 0.872±0.017), substantially outperforming models based on spirometry (AUC: 0.721) or nasometry (AUC: 0.625) alone. Interpretability analysis revealed that spatial domain features were the most critical, contributing 53.4% to the model’s decision-making, followed by temporal (25.0%) and airflow (21.6%) features. ConclusionST-EIT successfully captures the temporal, airflow, and spatial deviations in CP speech respiration that are undetectable by conventional methods—specifically, rapid phase transitions, hyperdynamic expiratory airflow, and regional ventilation heterogeneity. This study validates ST-EIT as a robust, non-invasive, and radiation-free tool for characterizing speech-respiratory dysfunction, offering high clinical value for bedside screening, rehabilitation planning, and longitudinal monitoring of patients with cleft palate.

ObjectiveCleft palate (CP) is a common congenital deformity often associated with velopharyngeal insufficiency (VPI), which disrupts the physiological coupling between respiration and speech. Conventional clinical assessments, such as nasometry and spirometry, provide limited static data and fail to visualize the dynamic spatiotemporal distribution of lung ventilation during phonation. This study introduces spatiotemporal electrical impedance tomography (ST-EIT) to evaluate speech-respiratory functional features in CP patients compared to normal controls (NC). The aim is to characterize multi-domain respiratory patterns and to validate an interpretable machine learning framework for providing objective, quantitative evidence for clinical assessment. MethodsSeventy-five participants were enrolled in this study, comprising 37 patients with surgically repaired CP and 38 healthy volunteers matched for age, gender, and body mass index (BMI). All subjects performed standardized sustained phonation tasks while undergoing synchronous monitoring with a 16-electrode EIT system and a pneumotachograph. A comprehensive feature engineering pipeline was developed to extract physiological parameters across 3 complementary domains. (1) Temporal domain: including inspiratory/expiratory phase duration (tPhase), time constants (Tau), and inspiratory-to-expiratory time ratios (TI/TE); (2) airflow domain: comprising mean flow, peak flow, and instantaneous flow at 25%, 50%, and 75% of tidal volume; and (3) spatial domain: quantifying global and regional tidal impedance variation (TIV), global inhomogeneity (GI), and center of ventilation (CoV). Extreme Gradient Boosting (XGBoost) classifiers were trained using 5 distinct data sources (Spirometry, Nasometry, Inspiratory-EIT, Expiratory-EIT, and fused ST-EIT). Model performance was rigorously evaluated via stratified 5-fold cross-validation, and Shapley additive explanations (SHAP) were employed to quantify global and local feature contributions. ResultsThe CP group exhibited a distinct respiratory phenotype compared to controls. In the temporal domain, CP patients showed significantly shorter inspiratory (1.60 s vs.1.85 s, P<0.001) and expiratory phase durations (2.45 s vs. 3.95 s, P<0.001), indicating a rapid, shallow breathing rhythm. In the airflow domain, while inspiratory flows were comparable, the CP group demonstrated significantly elevated mean and peak flows during the expiratory phase (P<0.001), reflecting compensatory respiratory effort. Spatially, CP patients presented significant ventilation redistribution, characterized by higher regional TIV in the right-anterior (ROI1) and left-posterior (ROI4) quadrants, but lower TIV in the left-anterior (ROI2) quadrant. In terms of diagnostic accuracy, the multi-modal ST-EIT model achieved the highest performance (AUC: 0.915±0.012, Accuracy: 0.843±0.019, F1-score: 0.872±0.017), substantially outperforming models based on spirometry (AUC: 0.721) or nasometry (AUC: 0.625) alone. Interpretability analysis revealed that spatial domain features were the most critical, contributing 53.4% to the model’s decision-making, followed by temporal (25.0%) and airflow (21.6%) features. ConclusionST-EIT successfully captures the temporal, airflow, and spatial deviations in CP speech respiration that are undetectable by conventional methods—specifically, rapid phase transitions, hyperdynamic expiratory airflow, and regional ventilation heterogeneity. This study validates ST-EIT as a robust, non-invasive, and radiation-free tool for characterizing speech-respiratory dysfunction, offering high clinical value for bedside screening, rehabilitation planning, and longitudinal monitoring of patients with cleft palate.

The Type III Secretion System (T3SS) serves as a pivotal virulence apparatus for numerous Gram-negative bacterial pathogens, enabling them to infect both animal and plant hosts. Functioning as a molecular syringe, the T3SS directly translocates bacterial effector proteins from the bacterial cytoplasm into the interior of eukaryotic host cells. These effectors are central weapons that precisely manipulate a wide spectrum of host cellular physiological processes, ranging from cytoskeletal dynamics to immune signaling, to establish a favorable niche for bacterial survival and proliferation. Among the diverse arsenal of T3SS effectors, the YopJ family constitutes a critical group of virulence factors. Members of this family are characterized by a conserved catalytic triad structure—a hallmark of the CE clan of cysteine proteases that has been evolutionarily repurposed to confer acetyltransferase activity. A defining and intriguing feature of these enzymes is their stringent dependence on a host-derived eukaryotic cofactor, inositol hexakisphosphate (IP₆), for allosteric activation. This requirement acts as a sophisticated molecular safeguard, ensuring enzymatic activity only within the appropriate host environment, thereby preventing detrimental effects on the bacterium itself. While seminal studies on individual members such as Yersinia’s YopJ and Salmonella’s AvrA have provided deep mechanistic insights, a systematic and integrative understanding of the structure-function relationships across the entire family remains fragmented. Key questions persist regarding how a conserved catalytic core has diverged to recognize distinct host substrates in different kingdoms of life. To address this gap, this article provides a systematic review of the YopJ family, focusing on three interconnected aspects: their structural features, their catalytic mechanism, and their divergent immunosuppressive strategies in animal versus plant hosts. By conducting a comparative analysis of the sequences and resolved three-dimensional structures of three representative members (e.g., HopZ1a, PopP2, AvrA), we elucidate regions of significant variation embedded within the conserved core catalytic architecture. These variable regions, often involving surface loops and substrate-binding interfaces, are crucial determinants of target specificity and functional specialization. The functional divergence of this effector family is most apparent when comparing their modes of action in different hosts. In animal hosts, YopJ-family effectors primarily sabotage innate immune signaling pathways. They achieve this by acetylating key serine and threonine residues within the activation loops of critical kinases in the MAPK and NF‑κB pathways. This post-translational modification blocks the phosphorylation and subsequent activation of these kinases, leading to potent suppression of inflammatory cytokine production. Conversely, in plant hosts, the strategy broadens to dismantle the two-tiered plant immune system. YopJ homologs target a more diverse set of substrates, including immune-associated receptor-like cytoplasmic kinases (RLCKs), microtubule networks via tubulin acetylation (which disrupts cellular trafficking and signaling), and transcription factors central to defense gene regulation. This multi-target approach effectively suppresses both Pattern-Triggered Immunity (PTI) and Effector-Triggered Immunity (ETI). In conclusion, this synthesis aims to deepen the mechanistic understanding of YopJ family-mediated pathogenesis by integrating structural biology with cellular function across host kingdoms. Elucidating the precise molecular basis for substrate selection—how conserved platforms achieve target diversity—is a major frontier. Furthermore, this knowledge provides a vital theoretical foundation for developing novel anti-virulence strategies. Targeting the conserved IP₆-binding pocket or the catalytic acetyltransferase activity itself represents a promising avenue for designing broad-spectrum inhibitors that could disarm this critical family of bacterial effectors, potentially offering new therapeutic approaches against a range of pathogenic bacteria.

Osteoarthritis (OA), a highly prevalent degenerative joint disease worldwide, is defined by articular cartilage degradation, abnormal bone remodeling, and persistent chronic inflammation. It severely compromises patients’ quality of life, and currently, there is no radical cure. Abnormal mechanical stress is widely regarded as a core driver of OA pathogenesis, and the exploration of mechanical signal perception and transduction mechanisms has become crucial for deciphering OA’s pathophysiological processes. Piezo1, a key mechanosensitive cation channel belonging to the Piezo protein family, has recently gained significant attention due to its pivotal role in mediating cellular responses to mechanical stimuli in joint tissues. This review systematically examines Piezo1’s expression patterns, regulatory mechanisms, and pathological functions in OA, with a particular focus on its dual roles in modulating chondrocyte homeostasis and bone metabolism disorders, while also delving into the underlying molecular signaling pathways and potential therapeutic implications. Piezo1, consisting of approximately 2 500 amino acids and forming a unique trimeric propeller-like structure, is widely expressed in chondrocytes, osteocytes, mesenchymal stem cells, and synovial cells. It exhibits permeability to cations such as Ca²⁺, K⁺, and Na⁺, and directly responds to membrane tension changes induced by mechanical stimuli like fluid shear stress and mechanical overload. In OA patients and animal models, Piezo1 expression is significantly upregulated, especially in cartilage regions subjected to abnormal mechanical stress (e.g., human temporomandibular joint cartilage). This overexpression is closely associated with aggravated cartilage degeneration, increased chondrocyte apoptosis, accelerated cellular senescence, and intensified inflammatory responses. Mechanical overload and pro-inflammatory cytokines (e.g., IL-1β) are key inducers of Piezo1 upregulation: IL-1β activates the PI3K/AKT/mTOR signaling pathway to enhance Piezo1 expression, forming a pathogenic positive feedback loop that inhibits chondrocyte autophagy, promotes apoptosis, and further accelerates joint degeneration. Mechanistically, Piezo1 mediates OA progression through multiple interconnected pathways. When activated by mechanical stress, Piezo1 triggers excessive Ca²⁺ influx, leading to endoplasmic reticulum stress (ERS) and mitochondrial dysfunction, which directly induce chondrocyte apoptosis. This process involves the activation of downstream signaling cascades such as cGAS-STING and YAP-MMP13/ADAMTS5. YAP, a transcriptional regulator, upregulates the expression of matrix metalloproteinase 13 (MMP13) and aggrecanase (ADAMTS5), thereby accelerating cartilage matrix degradation. Additionally, Piezo1-driven Ca²⁺ overload promotes the accumulation of reactive oxygen species (ROS) and upregulates senescence markers (p16 and p21), accelerating chondrocyte senescence via the p38MAPK and NF-κB pathways. Senescent chondrocytes secrete senescence-associated secretory phenotype (SASP) factors (e.g., IL-6, IL-1β), further amplifying joint inflammation. In terms of bone metabolism, Piezo1 maintains joint homeostasis by promoting the differentiation of fibrocartilage stem cells into chondrocytes and balancing bone formation and resorption through regulating the FoxC1/YAP axis and RANKL/OPG ratio. Therapeutically, targeting Piezo1 shows promising potential. Preclinical studies have demonstrated that Piezo1 inhibitors (e.g., GsMTx4) can reduce joint damage and alleviate pain in OA mice. Simultaneously, siRNA-mediated co-silencing of Piezo1 and TRPV4 (another mechanosensitive channel) decreases intracellular Ca²⁺ concentration, inhibits chondrocyte apoptosis, and promotes cartilage repair. Conditional knockout of Piezo1 using Gdf5-Cre transgenic mice alleviates cartilage degeneration in post-traumatic OA models by downregulating MMP13 and ADAMTS5 expression. Despite existing challenges, such as off-target effects of inhibitors, inefficient local drug delivery, and interindividual genetic variability, strategies like developing selective Piezo1 antagonists, optimizing targeted nanocarriers, and combining Piezo1-targeted therapy with physical therapy provide viable avenues for clinical translation. The authors propose that Piezo1 serves as a critical therapeutic target for OA, and future research should focus on deciphering its context-dependent regulatory networks, developing tissue-specific intervention strategies, and validating their efficacy and safety in clinical trials to address the unmet medical needs of OA patients.

Hepatic alveolar echinococcosis is a highly invasive zoonotic parasitic disease with poor prognosis. Surgical intervention serves as the pivotal approach to achieve radical cure and improve the prognosis of hepatic alveolar echinococcosis patients. In recent years, with the popularization of the concept of precision surgery and the development of the multidisciplinary diagnosis and treatment model, the surgical treatment strategies for hepatic alveolar echinococcosis have been continuously enriched, and the selection of surgical procedures has become increasingly diversified. Although key surgical techniques such as radical hepatectomy, autologous liver transplantation and allogeneic liver transplantation have achieved remarkable progress in clinical application, many insurmountable challenges still remain. Therefore, by sorting out the latest evidence-based advances in the field of surgical treatment for hepatic alveolar echinococcosis, this article focuses on discussing the application status and bottlenecks of radical hepatectomy, autologous liver transplantation and allogeneic liver transplantation in hepatic alveolar echinococcosis, aiming to provide a reference for the clinical treatment of hepatic alveolar echinococcosis.

ObjectiveThis study aims to explore the mechanism by which Yishen Huoxue Tongqiao Formula improves metabolism-neuroplasticity and treats unilateral vestibular labyrinth destruction by regulating the metabolic balance of glutamate （Glu）/γ-aminobutyric acid （GABA）. Methods48 Sprague-Dawley （SD） adult rats were randomly divided into the sham operation group， model group， Yishen Huoxue Tongqiao Formula groups with low， medium， and high doses （9.20， 18.39， 36.78 g·kg^-1）， and betahistine group （1.62 mg·kg^-1）. A unilateral vestibular labyrinth destruction （vestibular dysfunction） model was established by intratympanic injection of chloroform into the right ear， while the control group received intratympanic injection of normal saline. Drugs were administered once daily for seven consecutive days. During the period， behavioral tests were performed to evaluate the behaviors of rats after unilateral vestibular labyrinth destruction. Hematoxylin-eosin （HE） staining and Nissl staining were used to observe the neuronal morphology in the medial vestibular nucleus. Golgi staining was employed to assess the number of dendritic spines of neurons in the medial vestibular nucleus. Ultra-performance liquid chromatography-tandem mass spectrometry （LC-ESI-MS/MS） was utilized to detect Glu/GABA. Immunofluorescence and immunohistochemistry were used to detect the expressions of neuronal nuclei （NeuN）， growth-associated protein 43 （GAP-43）， and glial fibrillary acidic protein （GFAP）. Western blot and real-time fluorescent quantitative polymerase chain reaction （Real-time PCR） were applied to determine the expressions of glutamate-immunoreactive （Glu-IR）， GABA， GFAP， postsynaptic density protein 95 （PSD-95）， and GAP-43. ResultsCompared with the sham operation group， the model group presented with head deviation， balance disorder， increased tail suspension score， nuclear consolidation of medial vestibular nerve neurons， and decreased Nissl bodies （P<0.01）. The number of dendritic spines in neurons and NeuN-positive cells decreased. The content of Glu decreased. The content of GABA increased （Glu/GABA decreased）. The expression of GAP-43 was down-regulated， and GFAP was up-regulated （P<0.05， P<0.01）. The expressions of Glu-IR， PSD-95， and GAP-43 proteins， as well as Glu-IR mRNA decreased， while the expressions of GABA and GFAP proteins and mRNA increased （P<0.05， P<0.01）. Compared with those in the model group， the head deviation， imbalanced behavior， and tail suspension scores in each treatment group decreased， with alleviated neuronal injury and recovered Nissl bodies （P<0.01）. The number of dendritic spines of neurons increased， and the number of NeuN-positive cells rebounded. The content of Glu increased， and the content of GABA decreased （Glu/GABA increased）. GFAP was down-regulated， and GAP-43 was up-regulated （P<0.05， P<0.01）. The expressions of Glu-IR， PMD-95， and GAP-43 proteins， as well as Glu-IR mRNA increased， while the expressions of GABA and GFAP proteins and mRNA decreased. The effect was more significant in the high-dose group （P<0.01）. ConclusionThe Yishen Huoxue Tongqiao Formula can alleviate vestibular dysfunction， and its mechanism may be associated with regulating the metabolic balance of Glu/GABA， mitigating neural damage， improving synaptic plasticity （promoting GAP-43 expression and inhibiting GFAP expression）， and facilitating vestibular compensation.

ObjectiveTo investigate the present situation of input, output, outcome and impact of all registered community-based rehabilitation stations in Inner Mongolia in China, and analyze how the input predict the output, outcome and impact. MethodsFrom March 1st to April 30th, 2025, a questionnaire survey was conducted on all registered community-based rehabilitation stations in Inner Mongolia, covering four dimensions: input, output, outcome and impact. A total of 1 365 questionnaires were distributed. The input included four items: laws and policies, human resources, equipment and facilities, and rehabilitation information management. The output included two items: technical paths and benefits/effectiveness. The outcome included three items: coverage rates, rehabilitation interventions and functional results. The impact included two items: health and sustainability. Each item contained several questions, all of which were described in a positive way. Each question was scored from one to five. A lower score indicated that the situation of the community-based rehabilitation station was more in line with the content described in the question. Regression analysis was performed using the total score of each item of input dimension as independent variables, and the total scores of the output, outcome and impact dimensions as dependent variables. ResultsA total of 1 262 valid questionnaires were collected. The mean values of input, output, outcome and impact of community-based rehabilitation stations were 1.827 to 1.904, with coefficient of variation of 45.892% to 49.239%. The regression analysis showed that, rehabilitation information management, human resources, and laws and policies significantly predicted the output dimension (R² = 0.910, P < 0.001). Meanwhile, all four items in the input dimension predicted both the outcome (R² = 0.850, P < 0.001) and impact dimensions (R² = 0.833, P < 0.001). ConclusionInput, output, outcome and impact of the community-based rehabilitation stations in Inner Mongolia were generally in line with the content of the questions, although some imbalances were observed. Additionally, the input of community-based rehabilitation stations could significantly predict their output, outcome and impact.

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.