ObjectiveHistorically documented Zhejiang Atractylodis Macrocephalae Rhizoma（Baizhu） possesses premium characteristics such as phoenix-like head and crane-like neck， pronounced sweetness， and fragrant aroma. However， its current market circulation is low， and the processed products with Zhejiang-style characteristics are at the risk of being lost. This study aims to preserve the ancient Zhejiang-style processing techniques and evaluate them using modern scientific methods. MethodsMultidimensional intelligent sensory evaluation was used to digitally characterize the "quality-structure" of the external appearance of nine-steamed and nine-processed Baizhu medicinal materials（intermediate processed products） and the "odor-taste" of the internal quality of its decoction pieces（slices）， and the appearance parameters were digitally characterized by colorimeter， texture analyzer， electronic nose and electronic tongue， the chemical composition was analyzed via ultra-performance liquid chromatography-quadrupole-time-of-flight mass spectrometry（UPLC-Q-TOF-MS/MS）. Then， cluster analysis on the differences in odor between the medicinal materials（intermediate processed products） and decoction pieces（slices） of nine-steamed and nine-processed Baizhu was conducted， as well as the differences in taste between water-soluble and alcohol-soluble extracts of the decoction pieces（slices）， and the correlation analysis of chroma value-alcohol-soluble extract content-component response value. ResultsThe nine-steamed and nine-processed Baizhu had a dark brown to black epidermis， a brownish-yellow to brownish-gray cross-section， a slightly tough texture， a faint odor， and a slightly sweet， bitter and pungent taste. Texture analyzer measurements revealed minimal adhesion and maximum recovery in the middle section of the characteristic processed Baizhu， consistent with the processing endpoint of thorough steaming and cooking. The head section showed the highest internal hardness， elasticity and chewiness， indicating a denser texture in this area. The electronic nose sensor could clearly distinguish the difference between the medicinal materials and its decoction pieces， with a more significant clustering effect at 60 ℃ for 30 minutes compared to ambient temperature headspace for 2 hours， highlighting the significant impact of the baking degree before slicing on the quality. The electronic tongue taste signal map clearly distinguished the differences between water-soluble and alcohol-soluble extracts of nine-steamed and nine-processed Baizhu decoction pieces， and the addition of auxiliary materials during processing could enhance its alcohol-soluble extract content. A total of 82 chemical components were identified in the characteristic processed Baizhu. After processing， the contents of 58 components increased， while 24 components decreased. Correlation analysis revealed significant negative correlations（P<0.01） between ethanol-soluble extract content and colorimetric values of brightness（L^*）， yellow-bule value（b^*）， and total color difference（E^*_ab）. E^*_ab showed marked negative correlations（P<0.05） with the response values of isochlorogenic acid A and C. ConclusionThis study establishes a modern intelligent sensory evaluation model for multidimensional holographic characterization of nine-steamed and nine-processed Baizhu， clarifying the correlation between increased isochlorogenic acid content and the visual color appearance after different steaming cycles， as well as its intrinsic alcohol-soluble extracts. This provides a reference for quality evaluation and processing standards of the Zhejiang-style characteristic processed products.

Objective: To investigate the prevalence of undernutrition and its associated factors among left behind and non left behind primary and secondary school students in the Nutrition Improvement Program for Rural Compulsory Education Students (NIPRCES) areas of central and western China, so as to provide evidence for improving the nutritional status of children and adolescents. Methods: A survey was conducted among 123 782 students selected by random cluster sampling method in grades 3-9 from NIPRCES in central (Hebei, Shanxi, Heilongjiang, Jilin, Anhui, Jiangxi, Henan, Hunan, Hubei, and Hainan) and western (Gansu, Guangxi, Inner Mongolia, Ningxia, Tibet, Shaanxi, Guizhou, Sichuan, Xinjiang, the Xinjiang Production and Construction Corps, Yunnan, Qinghai, and Chongqing) China in 2023. Anthropometric measurements and questionnaires were used to assess nutritional and dietary status. The prevalence of undernutrition was compared between left behind and non left behind students by Chi square test, and associated factors were analyzed by three level Logistic mixed effects model. Results: The prevalence of undernutrition was 8.5% (4 326) in left behind students and 8.1% (5 905) in non left behind students. Three level Logistic mixed effect model analysis showed that whether left behind or non left behind, the undernutrition rates of primary and secondary students in western regions were higher than those of students in central regions [ OR (95% CI )=1.72(1.57-1.87),2.25(2.07- 2.43 )]; the undernutrition risk was lower for those whose fathers had a cultural level of high school or above [ OR (95% CI )=0.69(0.62-0.77),0.90(0.82-0.98)] or junior high school [ OR (95% CI )=0.72(0.66-0.79),0.92(0.85-0.99)] compared to those with primary school or below; picky eating or selective eating increased the risk of undernutrition [ OR (95% CI )=2.36(2.07-2.68),2.28(2.04-2.55)], and primary and secondary school students without nutritional content in health education classes had higher rates of undernutrition [ OR (95% CI )=1.12(1.03-1.23),1.09(1.01-1.17)](all P <0.05). Conclusion The prevalence of undernutrition is slightly higher in left behind primary and secondary students than in non left behind primary and secondary students in central and western NIPRCES areas, with variations across different characteristics.

Objective: To analyze the trends of the frequency of meat, egg, and milk consumption among rural primary and junior high school students in central and western China covered by the Nutrition Improvement Program for Rural Compulsory Education Students (NIPRCES) from 2015 to 2023, so as to provide basis for formulating more targeted nutrition intervention policies and health education strategies. Methods: Using data from six rounds of monitoring and evaluation (2015-2021 and 2023), the study included 323 870 students from grade 3 to 9 across 22 provinces (autonomous regions and municipalities) in central and western China. The consumption frequencies of meat, egg, and milk over the past week were collected via questionnaires. The Cochran-Armitage trend test was used to analyze temporal trends, and multivariable Logistic regression models were employed to analyze factors associated with the frequency of meat, egg and milk consumption and to test for interaction effects between the year and gender, region, and grade level. Results: From 2015 to 2023, the proportion of students consuming meat, egg, and milk ≥1 time/day increased from 23.20 %, 10.71%, and 0.74% to 35.53%, 22.09%, and 26.63%, respectively. Trend tests indicated a significant upward trend for the daily intake of all three food categories for meat, egg and milk over the years ( Z =67.18, 64.90, 93.14, all P <0.01). Multivariable Logistic regression analysis showed that the daily meat intake was lower in the central region than in the western region ( OR=0.77, 95%CI =0.76-0.78), whereas the daily intake of eggs ( OR=1.19, 95%CI =1.17-1.22) and milk ( OR= 1.27 , 95%CI =1.24-1.29) was higher in the central region (all P <0.05). Compared with grade 3-4 students, junior high school students had lower daily intake of meat, eggs, and milk≥1 time/day ( OR =0.95, 0.77, 0.77, all P <0.05), with a declining trend as grade increased. Girls also had lower daily intake of meat, eggs, and milk ≥1 time/day than boys ( OR =0.95，0.93，0.91, all P < 0.05). Significant interactions were observed between year and region, as well as between year and grade (all P <0.05). Conclusion From 2015 to 2023, the NIPRCES improved the intake level of among rural students, but the situation of relatively insufficient intake of egg and milk among females, junior high school students and those in the western region still exists.

ObjectiveTo investigate the impact of aberrant SKI expression and its mutations on the biological characteristics of cholangiocarcinoma cell lines QBC939 and RBE， and to explore the underlying molecular mechanisms. MethodsThe Gene Expression Profiling Interactive Analysis 2 （GEPIA2） database was utilized to analyze SKI expression and its clinical relevance in cholangiocarcinoma patients. Lentiviral transduction was employed to establish QBC939 and RBE cell lines with stable SKI overexpression， mutation， or knockdown. Cell proliferation was assessed using CCK-8， colony formation， and EdU assays. Apoptosis and cell cycle distribution were analyzed by flow cytometry. Cell migration was evaluated using Transwell and wound healing assays. The effect of SKI over-expression， mutation， or knockdown on key proteins （SMAD2， SMAD3， SMAD4） in the transforming growth factor-β （TGF-β）/Small mothers against decapentaplegic （SMAD） signaling pathway was examined by Western blot. ResultsCompared to SKI overexpression alone， the introduction of SKI mutations significantly promoted S-phase progression， enhanced proliferation and migration， and inhibited apoptosis. Mechanistically， SKI mutations suppressed the phosphorylation of SMAD2 and SMAD3 proteins， thereby inhibiting the transcriptional activity of the TGF-β signaling pathway. Conversely， SKI knockedown produced the opposite effects. ConclusionSKI gene mutation acts as a gain-of-function genetic alteration， exerting an oncogenic role in cholangiocarcinoma cells. The primary mechanism involves the inhibition of the TGF-β/SMAD signaling pathway， which in turn promotes proliferation and cell cycle progression， and suppresses apoptosis in QBC939 and RBE cells， ultimately driving tumor progression.

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.

ObjectiveTo clarify whether METTL14 mediates the core role of acupuncture at Neiguan (PC6) in promoting myelination and improving behavior in young autistic rats through gene intervention technology. MethodsThe ASD model was established by intraperitoneal injection of valproic acid (VPA) in pregnant rats. Male offspring were intracerebroventricularly injected with adenovirus-packaged METTL14 shRNA (sh-METTL14) or its control (sh-NC) on postnatal day 1, with a model group set as well. Subsequently, the juvenile rats were divided into model group, acupuncture group, acupuncture+sh-NC group, and acupuncture+sh-METTL14 group. The acupuncture group received acupuncture at Neiguan (PC6) from postnatal day 7, once daily for 21 consecutive days. Neurobehavioral changes were evaluated by behavioral tests; METTL14 knockdown efficiency and the expression of METTL14, METTL3, and PTEN were detected by quantitative real-time PCR (qRT-PCR) and Western blot (WB); PTEN m⁶A levels were measured by RNA immunoprecipitation-qPCR (RIP-qPCR); myelin ultrastructure, expression of myelin basic protein (MBP) and neurofascin 155 (NF155), and dendritic spine density were observed using transmission electron microscopy (TEM), enzyme-linked immunosorbent assay (ELISA), immunofluorescence, qRT-PCR, and primary neuron culture. ResultsBehaviorally, knockdown of METTL14 significantly counteracted the beneficial effects of acupuncture in improving self-grooming, open field exploration, three-chamber social interaction, and Morris water maze learning and memory (P<0.05, P<0.01). Compared with the acupuncture+sh-NC group, the acupuncture+sh-METTL14 group showed significantly decreased mRNA and protein expression of hippocampal METTL14 (P<0.01), and the upregulating effects of acupuncture on METTL3 and PTEN expression were reversed (P<0.01). Meanwhile, knockdown of METTL14 significantly inhibited the acupuncture-induced increase in PTEN m⁶A levels (P<0.01). Morphologically, knockdown of METTL14 attenuated the improvement of myelin structure by acupuncture, reversed the downregulation of MBP and upregulation of NF155 induced by acupuncture, and blocked the increase in dendritic spine density (P<0.05, P<0.01). ConclusionMETTL14 is a key molecule mediating the therapeutic effect of acupuncture at Neiguan. Acupuncture at Neiguan upregulates METTL14, thereby enhancing m⁶A methylation modification of PTEN mRNA to stabilize its expression, ultimately promoting myelin development and improving behavioral symptoms in ASD juvenile rats. This preliminarily reveals the modern biological connotation of “opening Xuanfu and dredging myelin”.

ObjectiveTo clarify whether METTL14 mediates the core role of acupuncture at Neiguan (PC6) in promoting myelination and improving behavior in young autistic rats through gene intervention technology. MethodsThe ASD model was established by intraperitoneal injection of valproic acid (VPA) in pregnant rats. Male offspring were intracerebroventricularly injected with adenovirus-packaged METTL14 shRNA (sh-METTL14) or its control (sh-NC) on postnatal day 1, with a model group set as well. Subsequently, the juvenile rats were divided into model group, acupuncture group, acupuncture+sh-NC group, and acupuncture+sh-METTL14 group. The acupuncture group received acupuncture at Neiguan (PC6) from postnatal day 7, once daily for 21 consecutive days. Neurobehavioral changes were evaluated by behavioral tests; METTL14 knockdown efficiency and the expression of METTL14, METTL3, and PTEN were detected by quantitative real-time PCR (qRT-PCR) and Western blot (WB); PTEN m⁶A levels were measured by RNA immunoprecipitation-qPCR (RIP-qPCR); myelin ultrastructure, expression of myelin basic protein (MBP) and neurofascin 155 (NF155), and dendritic spine density were observed using transmission electron microscopy (TEM), enzyme-linked immunosorbent assay (ELISA), immunofluorescence, qRT-PCR, and primary neuron culture. ResultsBehaviorally, knockdown of METTL14 significantly counteracted the beneficial effects of acupuncture in improving self-grooming, open field exploration, three-chamber social interaction, and Morris water maze learning and memory (P<0.05, P<0.01). Compared with the acupuncture+sh-NC group, the acupuncture+sh-METTL14 group showed significantly decreased mRNA and protein expression of hippocampal METTL14 (P<0.01), and the upregulating effects of acupuncture on METTL3 and PTEN expression were reversed (P<0.01). Meanwhile, knockdown of METTL14 significantly inhibited the acupuncture-induced increase in PTEN m⁶A levels (P<0.01). Morphologically, knockdown of METTL14 attenuated the improvement of myelin structure by acupuncture, reversed the downregulation of MBP and upregulation of NF155 induced by acupuncture, and blocked the increase in dendritic spine density (P<0.05, P<0.01). ConclusionMETTL14 is a key molecule mediating the therapeutic effect of acupuncture at Neiguan. Acupuncture at Neiguan upregulates METTL14, thereby enhancing m⁶A methylation modification of PTEN mRNA to stabilize its expression, ultimately promoting myelin development and improving behavioral symptoms in ASD juvenile rats. This preliminarily reveals the modern biological connotation of “opening Xuanfu and dredging myelin”.

GPR126, also known as ADGRG6, is one of the most deeply studied aGPCRs. Initially, GPR126 was thought to be a receptor associated with muscle development and was primarily expressed in the muscular and skeletal systems. With the deepening of research, it was found that GPR126 is expressed in multiple mammalian tissues and organs, and is involved in many biological processes such as embryonic development, nervous system development, and extracellular matrix interactions. Compared with other aGPCRs proteins, GPR126 has a longer N-terminal domain, which can bind to ligands one-to-one and one-to-many. Its N-terminus contains five domains, a CUB (complement C1r/C1s, Uegf, Bmp1) domain, a PTX (Pentraxin) domain, a SEA (Sperm protein, Enterokinase, and Agrin) domain, a hormone binding (HormR) domain, and a conserved GAIN domain. The GAIN domain has a self-shearing function, which is essential for the maturation, stability, transport and function of aGPCRs. Different SEA domains constitute different GPR126 isomers, which can regulate the activation and closure of downstream signaling pathways through conformational changes. GPR126 has a typical aGPCRs seven-transmembrane helical structure, which can be coupled to Gs and Gi, causing cAMP to up- or down-regulation, mediating transmembrane signaling and participating in the regulation of cell proliferation, differentiation and migration. GPR126 is activated in a tethered-stalk peptide agonism or orthosteric agonism, which is mainly manifested by self-proteolysis or conformational changes in the GAIN domain, which mediates the rapid activation or closure of downstream pathways by tethered agonists. In addition to the tethered short stem peptide activation mode, GPR126 also has another allosteric agonism or tunable agonism mode, which is specifically expressed as the GAIN domain does not have self-shearing function in the physiological state, NTF and CTF always maintain the binding state, and the NTF binds to the ligand to cause conformational changes of the receptor, which somehow transmits signals to the GAIN domain in a spatial structure. The GAIN domain can cause the 7TM domain to produce an activated or inhibited signal for signal transduction, For example, type IV collagen interacts with the CUB and PTX domains of GPR126 to activate GPR126 downstream signal transduction. GPR126 has homology of 51.6%-86.9% among different species, with 10 conserved regions between different species, which can be traced back to the oldest metazoans as well as unicellular animals.In terms of diseases, GPR126 dysfunction involves the pathological process of bone, myelin, embryo and other related diseases, and is also closely related to the occurrence and development of malignant tumors such as breast cancer and colon cancer. However, the biological function of GPR126 in various diseases and its potential as a therapeutic target still needs further research. This paper focuses on the structure, interspecies differences and conservatism, signal transduction and biological functions of GPR126, which provides ideas and references for future research on GPR126.

ObjectiveFor prepubertal and urgently treated malignant tumor patients, ovarian tissue cryopreservation and transplantation represent more appropriate fertility preservation methods. Current clinical practices often involve freezing ovarian tissue with high concentrations of cryoprotectants (CPAs) and thawing with water baths. These processes lead to varying degrees of toxicity and devitrification damage to ovarian tissue. Therefore, this paper proposes optimized methods for vitrification of ovarian tissues based on sodium alginate hydrogel encapsulation and magnetic induction nanowarming technology. MethodsFirstly, the study investigated the effects of sodium alginate concentration, the sequence of hydrogel encapsulation and CPAs loading on vitrification efficiency of encapsulated ovarian tissue. Additionally, the capability of sodium alginate hydrogel encapsulation to reduce the required concentration of CPAs was validated. Secondly, a platform combining water bath and magnetic induction nanowarming was established to rewarm ovarian tissue under various concentrations of magnetic nanoparticles and magnetic field strengths. The post-warming follicle survival rate, antioxidant capacity, and ovarian tissue integrity were evaluated to assess the efficacy of the method. ResultsThe study found that ovarian tissue encapsulated with 2% sodium alginate hydrogel exhibited the highest follicle survival rate after vitrification. The method of loading CPAs prior to encapsulation proved more suitable for ovarian tissue cryopreservation, effectively reducing the required concentration of CPAs by 50%. A combination of 8 g/L Fe₃O₄ nanoparticles and an alternating magnetic field of 300 Gs showed optimal warming effectiveness for ovarian tissue. Combining water bath rewarming with magnetic induction nanowarming yielded the highest follicle survival rate, enhanced antioxidant capacity, and preserved tissue morphology. ConclusionSodium alginate hydrogel encapsulation of ovarian tissue reduces the concentration of CPAs required during the freezing process. The combination of magnetic induction nanowarming with water bath provides an efficient method ovarian tissue rewarming. This study offers novel approaches to optimize ovarian tissues vitrification.

ObjectiveTo characterize the changes of metabolites of Blumea balsamifera in the process of sweating by non-targeted metabolomics， and to investigate the influence of sweating processing on the constituents of B. balsamifera. MethodsUltra performance liquid chromatography-quadrupole/electrostatic field orbitrap high resolution mass spectrometry（UPLC-Q-Exactive Orbitrap-MS） metabolomics was used to identify the metabolites in no sweating group（F1）， sweating 2 d group（F2） and sweating 4 d group（F3）， the differences of metabolites between the groups were compared by principal component analysis（PCA） and orthogonal partial least squares-discriminant analysis（OPLS-DA）， and differential metabolites were screened according to the variable importance in the projection（VIP） value>1 and P<0.05， and the pathway enrichment of the differential metabolites was analyzed by Kyoto Encyclopedia of Genes and Genomes（KEGG）. ResultsThe results of PCA and OPLS-DA showed a clear distinction between the three groups of samples， indicating significant differences in the compositions of the three groups of samples. A total of 433 differential metabolites were screened between the F1 and F2， with 154 up-regulated and 279 down-regulated， the significant up-regulated metabolites were tangeritin， 5-O-demethylnobiletin and so on， while the metabolites with significant down-regulation included alternariol， fortunellin， etc. A total of 379 differential metabolites were screened between the F2 and F3， with 150 up-regulated and 229 down-regulated， the significant up-regulated metabolites were isoimperatorin， helianyl octanoate and so on， and the significant down-regulated metabolites were hovenoside I， goyasaponin Ⅲ， etc. KEGG pathway enrichment analysis showed that tyrosine metabolism， isoquinoline alkaloid biosynthesis， phenylalanine， tyrosine and tryptophan biosynthesis， tryptophan metabolism， valine， leucine and isoleucine biosynthesis， pantothenate and coenzyme A biosynthesis may be the key pathways affecting metabolite differences of B. balsamifera after sweating treatment. ConclusionSweating can reduce the content of endophytic mycotoxins in B. balsamifera and has a great impact on the synthesis and metabolic pathways of total flavonoids and auxin. This study can provide a reference for the process research on the sweating conditions of B. balsamifera.