ObjectiveTo investigate the inhibitory effect of Biejiajian Pills （BJJW） on the growth of liver cancer， as well as its potential mechanism in mediating the AMP-activated protein kinase （AMPK）/mammalian target of rapamycin （mTOR） pathway through mitochondrial energy metabolism. MethodsHuman hepatoma Huh7 cells were used to establish a nude mouse model of subcutaneous xenograft tumor. A total of 18 tumor-bearing nude mice were randomly divided into model group， BJJW group （2.2 g/kg）， and metformin group （250 mg/kg）， and the corresponding drug was given by gavage for 14 consecutive days. Tumor volume and weight were monitored during the experiment； HE staining was used to observe histopathological changes； the levels of reactive oxygen species （ROS） and adenosine triphosphate （ATP） in tumor tissue were measured； immunohistochemistry and Western blotting were used to measure the expression levels of proteins associated with the AMPK/mTOR pathway. A one-way analysis of variance was used for comparison of normally distributed continuous data between multiple groups， and the Tukey’s test was used for further comparison between two groups； the Kruskal-Wallis H test was used for comparison of non-normally distributed continuous data between multiple groups， and the Dunn’s test was used for further comparison between two groups. ResultsCompared with the model group， the BJJW group had a tumor inhibition rate of 45.73%， with significant reductions in both tumor volume and weight （P<0.01）. Pathological examination showed that compared with the model group， the BJJW group had a significant reduction in the number of tumor cells and the presence of extensive necrosis. Mechanistic studies showed that compared with the model group， the BJJW group had a significant increase in ROS level （P<0.001） and a significant reduction in ATP level （P<0.001）， as well as significant increases in p-AMPK/AMPK ratio （0.81±0.20 vs 0.13±0.04， P<0.01） and p-ULK1/ULK1 ratio （0.69±0.17 vs 0.18±0.13， P<0.01） and a significant reduction in p-mTOR/mTOR ratio （1.34±0.16 vs 3.20±0.62， P<0.01）. ConclusionBJJW may inhibit the growth of liver cancer by inducing mitochondrial energy metabolism dysfunction， increasing the level of ROS， reducing the level of ATP， and activating the AMPK/mTOR signaling pathway.

ObjectiveTo establish an ultra-performance liquid fingerprint and multi-components determination method for Naomaili granules. To evaluate the quality of different batches by chemometrics， and the anti-neuroinflammatory effects of water extract and main components of Naomaili granules were tested in vitro. MethodsThe similarity and common peaks of 27 batches of Naomaili granules were evaluated by using Ultra performance liquid chromatography （UPLC） fingerprint detection. Ultra-performance liquid chromatography-tandem mass spectrometry （UPLC-MS/MS） technology was used to determine the content of the index components in Naomaili granules and to evaluate the quality of different batches of Naomaili granules by chemometrics. LPS-induced BV-2 cell inflammation model was used to investigate the anti-neuroinflammatory effects of the water extract and main components of Naomaili granules. ResultsThe similarity of fingerprints of 27 batches of samples was > 0.90. A total of 32 common peaks were calibrated， and 23 of them were identified and assigned. In 27 batches of Naomaili granules， the mass fractions of 14 components that were stachydrine hydrochloride， leonurine hydrochloride， calycosin-7-O-glucoside， calycosin，tanshinoneⅠ， cryptotanshinone， tanshinoneⅡ_A， ginsenoside Rb₁， notoginsenoside R₁， ginsenoside Rg₁， paeoniflorin， albiflorin， lactiflorin， and salvianolic acid B were found to be 2.902-3.498， 0.233-0.343， 0.111-0.301， 0.07-0.152， 0.136-0.228， 0.195-0.390， 0.324-0.482， 1.056-1.435， 0.271-0.397， 1.318-1.649， 3.038-4.059， 2.263-3.455， 0.152-0.232， 2.931-3.991 mg∙g^-1， respectively. Multivariate statistical analysis showed that paeoniflorin， ginsenoside Rg₁， ginsenoside Rb₁ and staphylline hydrochloride were quality difference markers to control the stability of the preparation. The results of bioactive experiment showed that the water extract of Naomaili granules and the eight main components with high content in the prescription had a dose-dependent inhibitory effect on the release of NO in the cell supernatant. Among them， salvianolic acid B and ginsenoside Rb₁ had strong anti-inflammatory activity， with IC₅₀ values of （36.11±0.15） mg∙L^-1 and （27.24±0.54） mg∙L^-1， respectively. ConclusionThe quality evaluation method of Naomaili granules established in this study was accurate and reproducible. Four quality difference markers were screened out， and eight key pharmacodynamic substances of Naomaili granules against neuroinflammation were screened out by in vitro cell experiments.

AIM:To investigate the prevalence of visual impairment and its correction status among junior high school students in Xi'an, so as to provide evidence for the development of targeted myopia prevention and control strategies.METHODS: A stratified cluster sampling design was adopted. From March to May 2025, students in grades 7-9 were recruited from three schools in Xi'an, Shaanxi Province, China: Dongfang Middle School, the Middle School Attached to Xi'an University of Technology, and the Xingqing Campus of the High School Affiliated to Xi'an Jiaotong University. In total, 3 974 students were invited, including 1 726 in grade 7, 1 206 in grade 8, and 1 042 in grade 9. The visual acuity was measured monocularly using a 5 m standard logarithmic visual acuity chart, with the fellow eye occluded; the line corresponding to the smallest optotype that could be correctly identified was recorded as the visual acuity value. Non-cycloplegic autorefraction was performed with a desktop autorefractor to obtain spherical equivalent(SE)values for refractive error screening.RESULTS: This study initially included 3 974 students, of whom 32 did not participate in the vision test, resulting in 3 942 students being included in the final analysis. Among them, 3 067(77.80%)were identified with poor vision. The prevalence of myopia was 81.47%(1 746)in males and 87.55%(1 575)in females(P<0.01). A stratified analysis by grade showed myopia rates of 81.72%(1 386)in junior grade one, 84.47%(1 017)in junior grade two, and 88.10%(918)in junior grade three, demonstrating a significant upward trend with increasing grade level(χ2=19.8484, P<0.01). Among the 3 321 myopic students, 2 287 adopted corrective measures. The rates of full correction, under-correction, and non-correction among all myopic students were 48.15%(1 599), 20.71%(688), and 31.14%(1 034), respectively. The rate of non-correction was significantly higher in male students than in females(32.70% vs 29.40%, χ2=4.2222, P<0.05).CONCLUSION: The findings indicate a high prevalence of visual impairment among junior high school students in Xi'an, coupled with suboptimal spectacle-wearing and full-correction rates. There is an urgent need for collaborative efforts across society, schools, and families to implement effective interventions to slow the onset and progression of myopia in this population.

ObjectiveTo investigate the effects and underlying mechanisms of Shenqi Wenfei prescription （SQWF） on chronic obstructive pulmonary disease （COPD）. MethodsA rat model of COPD with lung Qi deficiency was established using lipopolysaccharide （LPS） combined with cigarette smoke. Forty-eight SD rats were randomly divided into a blank group， a model group， low-， medium-， and high-dose SQWF groups （2.835， 5.67， 11.34 g·kg^-1）， and a Yupingfeng group （1.35 g·kg^-1）. Drug administration began on day 29 after modeling and continued for 2 weeks. The general condition of the rats was observed， and the lung function in each group was assessed. Hematoxylin-eosin （HE） staining was used to observe pathological changes in lung tissue. The proportion of inflammatory cells in bronchoalveolar lavage fluid （BALF） was measured. Apoptosis in lung tissue was examined by terminal deoxynucleotidyl transferase-mediated dUTP nick end-labeling （TUNEL） staining. The release level of lactate dehydrogenase （LDH） in BALF was detected by a microplate assay. Reactive oxygen species （ROS） levels in lung tissue were detected using fluorescent probes. The levels of malondialdehyde （MDA）， total superoxide dismutase （SOD）， and reduced glutathione （GSH） in BALF were measured by biochemical methods. Ultrastructural changes in lung cells were observed via transmission electron microscopy. Double immunofluorescence staining was performed to detect the expression of thioredoxin-interacting protein （TXNIP） and nucleotide-binding oligomerization domain-like receptor protein 3 （NLRP3） in lung tissue. Western blot analysis was used to detect the protein expression of TXNIP， NLRP3， apoptosis-associated speck-like protein containing a CARD （ASC）， cysteinyl aspartate-specific protease-1 （Caspase-1）， Caspase-1 p20， gasdermin D （GSDMD）， GSDMD N-terminal active fragment （GSDMD-N）， interleukin-1β （IL-1β）， and IL-18 in lung tissue. Serum IL-1β and IL-18 levels were measured by ELISA. ResultsCompared with the blank group， the model group showed lassitude， fatigue， tachypnea， and audible phlegm sounds， and lung function significantly declined （P0.01）. Pulmonary emphysema and inflammatory cell infiltration were obvious. The level of inflammatory cells in BALF increased significantly （P0.05）. The number of TUNEL-positive cells increased （P0.01）. Levels of LDH， ROS， and MDA in BALF increased significantly （P0.01）， while GSH and SOD activities decreased significantly （P0.01）. Lung tissue cells showed irregular morphology， swollen mitochondria， disrupted cell membranes， and abundant vesicles， i.e.， pyroptotic bodies. Protein levels of TXNIP， NLRP3， ASC， Caspase-1， Caspase-1 p20， GSDMD， GSDMD-N， IL-1β， and IL-18 in lung tissue were significantly elevated （P0.01）， and serum IL-1β and IL-18 levels also increased significantly （P0.01）. Compared with the model group， each medication group showed alleviation of qi deficiency symptoms and improved lung function （P0.01）. Pulmonary emphysema and inflammatory cell infiltration were reduced. Inflammatory cell levels decreased （P0.05）. The number of TUNEL-positive cells decreased significantly （P0.01）. Levels of LDH， ROS， and MDA decreased significantly （P0.05）， while GSH and SOD activities significantly increased （P0.01）. Morphological and structural damage in lung tissue was improved to varying degrees. Protein levels of TXNIP， NLRP3， ASC， Caspase-1， Caspase-1 p20， GSDMD， GSDMD-N， IL-1β， and IL-18 in lung tissue significantly decreased （P0.01）， and serum IL-1β and IL-18 levels also decreased significantly （P0.05）. ConclusionSQWF can improve lung function and alleviate inflammatory responses in COPD rats. Its mechanism may be related to regulating the ROS/TXNIP/NLRP3 pathway and inhibiting pyroptosis.

BACKGROUND:One of the advantages of clear aligner treatment is molar distalization.However,tooth tilting movement and loss of anterior anchorage may occur during treatment.There are few studies on whether these problems can be improved by selecting clear aligners with different thicknesses and edges to improve the clinical treatment effect.OBJECTIVE:To analyze the control ability of clear aligners with different thickness and edges on the central incisor,lateral incisor,and second molar when pushing the maxillary second molar distally by three-dimensional finite element analysis.METHODS:Three-dimensional finite element analysis models of bilateral maxillary second molar distalization with clear aligner,maxillary dentition,periodontal ligament,and alveolar bone with different thicknesses and margins were established by Mimics,Geomagic Wrap,3-matic and SolidWorks software,respectively.There were 16 combinations of four thicknesses(0.4,0.5,0.625,and 0.75 mm)and four margins(scallop,straight,straight extension 2 mm and straight extension 4 mm).The data were imported into Ansys Workbench software for design and solution.The mean value,peak value and distribution of the periodontal ligament equivalent stress of the second molar,the equivalent stress and the maximum initial displacement of the second molar,and the control ability of each appliance on the second molar,central incisor,and lateral incisor were analyzed.RESULTS AND CONCLUSION:(1)The mean equivalent stress of periodontal ligament of the second molar,the equivalent stress of the second molar and the maximum initial displacement of the second molar increased with the extension of the appliance edge and the increase of the thickness of the appliance in the 16 groups of models.(2)When the thickness of appliances was the same,the maximum equivalent stress of the second molar in the linear appliance group was the highest,and the maximum equivalent stress of the second molar in the linear extended appliance group was greater than that in the scallop appliance group.When the edge of the appliance was the same,the periodontal ligament equivalent stress peak of the second molar increased with the increase of the thickness of the appliance.The equivalent stress distribution in the periodontal ligament of the second molar in the linear extendable appliance group was more uniform than that in the scallop appliance group and the linear appliance group.(3)When the thickness of the appliance was the same,the scallop-shaped appliance had the worst control on the second molar.When the edge of the appliance was the same,with the increase of the thickness of the appliance,the control of the second molar by the linear extender appliance was gradually stronger than that by the linear appliance.The control of the central incisor was stronger and more stable with the linear extended 2 mm appliance,while the control of the lateral incisor was stronger and more stable with the linear appliance.(4)The results showed that when using clear aligners to push molars distally,extending the edge of the appliance could improve the control of the molars and reduce the tilting movement of the teeth.The design of a straight extension margin of 2 mm for the central incisor and a straight edge for the lateral incisor can enhance the control of the anchorage incisor and reduce the labial inclination of the anterior teeth.

BACKGROUND:Previous studies have found that neuronal damage caused by continuous excessive fluoride exposure is related to Ca2+overload,but the mechanism of Ca2+flow conversion between intracellular calcium stores and cell apoptosis damage is still unclear.OBJECTIVE:To investigate the effect of fluoride exposure on Ca2+transport channel proteins and apoptosis levels in the mitochondria-associated endoplasmic reticulum membrane of primary cultured neural cells.METHODS:Primary nerve cells of neonatal SD rats were cultured in vitro and identified by immunofluorescence staining with neuronal nucleus-specific antibody up to day 7.The nerve cells were divided into control group(containing 0 mmol/L sodium fluoride),low fluoride group(containing 0.5 mmol/L sodium fluoride),and high fluoride group(containing 1 mmol/L sodium fluoride).The cell morphological changes were observed by light microscope 24 hours after fluorine exposure.The expression levels of apoptosis-related protein BAX/BCL-2 and calcium transfer-related pathways VDAC1,GRP 75,and IP3R were detected using western blot assay.The expression levels of VDAC1,GRP 75,and IP3R mRNA were detected by RT-PCR.Ca2+levels were detected by Rhood-2AM Ca2+probe.Mitochondrial membrane potential detection kit was used to detect the change in mitochondrial membrane potential.The level of apoptosis was determined by flow cytometry and TUNEL staining.RESULTS AND CONCLUSION:(1)The purity of neurons cultured on day 7 had been determined to be over 90％,with few impurities,good growth status,and tight cell network connections,meeting the requirements of subsequent experiments.(2)Compared with the control group,growth of neural cell clusters in the low-fluoride group and the high-fluoride group increased;the processes were broken;the cell body was rounded,and the connection network between cells was destroyed.Compared with the low-fluoride group,the cell damage changes in the high-fluoride group were more obvious.(3)Compared with the control group,the protein expressions of VDAC1,GRP75,and IP3R were increased in the low-fluoride group and the high-fluoride group(P＜0.05),and the ratio of apoptosis-related protein BAX/BCL-2 was increased(P＜0.05).Compared with the control group,the expression of VDAC1 and GRP75 mRNA in the low-fluoride group was significantly increased(P＜0.05);the expression levels of VDAC1,GRP75,and IP3R mRNA in the high-fluoride group were significantly increased(P＜0.01).(4)The level of cell apoptosis increased significantly after fluoride exposure,and the high-fluoride group was significantly higher than the control and low-fluoride groups(P＜0.01).(5)After fluoride exposure,the concentration of mitochondrial Ca2+in nerve cells increased significantly(P＜0.05),the mitochondrial membrane potential decreased(P＜0.01),and the degree of damage in the high-fluoride group was more obvious(P＜0.05).The results show that fluoride exposure impairs the morphological structure of primary neural cells,resulting in upregulation of Ca2+transfer pathway protein expression between the endoplasmic reticulum and mitochondria,mitochondrial Ca2+overload,mitochondrial damage,and increased levels of apoptosis.

BACKGROUND:Studies have confirmed that Pterocarya hupehensis skan total flavonoids(PHSTF)can improve the level of collagen-induced arthritis in rats,but there is still a lack of research on the regulation of Wnt/β-catenin signaling pathway in fibroblast-like synoviocytes and its effect on related cell functions.OBJECTIVE:To investigate the effect and mechanism of PHSTF on lipopolysaccharide-induced proliferation,migration and apoptosis of fibroblast-like synoviocytes based on the Wnt/β-catenin signaling pathwayMETHODS:Fibroblast-like synoviocytes were divided into control group,lipopolysaccharide group,lipopolysaccharide+low-,medium-,and high-dose PHSTF groups(10,20,and 40 μg/mL),lipopolysaccharide+Wnt pathway inhibitor DKK1 group,and lipopolysaccharide+Wnt pathway inhibitor DKK1+high-dose PHSTF group(40 μg/mL).The cell counting kit-8 method was used to detect the effect of PHSTF on the viability of fibroblast-like synoviocytes,and the final drug concentration and time were screened.Flow cytometry was used to detect the apoptosis of fibroblast-like synoviocytes.Cell scratch assay,EDU staining and cell cloning assay were used to detect the migration and proliferation of fibroblast-like synoviocytes.Western blot assay was used to detect the protein expression levels of Wnt3a,β-catenin,tumorigenic genes,matrix metalloproteinase 2,matrix metalloproteinase 9,Bax and Bcl-2 in fibroblast-like synoviocytes.RESULTS AND CONCLUSION:(1)Compared with the control group,the cell viability decreased significantly when the concentration of PHSTF was＞40 μg/mL(P＜0.01).Therefore,the drug concentration of≤40 μg/mL was selected for subsequent experiments.(2)Compared with the lipopolysaccharide group,the wound healing rate,cell clone formation rate and the number of EDU-positive cells in the low-,medium-and high-dose PHSTF groups were significantly reduced,while the apoptosis rate was significantly increased(P＜0.05-0.01).(3)Western blot results showed that compared with the lipopolysaccharide group,low-,medium-and high-dose PHSTF significantly inhibited cellular Wnt3a,β-catenin,cellular tumorigenic genes,matrix metalloproteinase 2,matrix metalloproteinase 9,and Bcl-2 protein expression,and promoted the expression of Bax protein(P＜0.01).(4)Compared with the DKK1 group,the combination of DKK1 and high-dose PHSTF significantly inhibited the protein expression of Wnt3a,β-catenin,matrix metalloproteinase 2,matrix metalloproteinase 9 and Bcl-2 protein expression and promoted the protein expression of Bax(P＜0.01).To conclude,PHSTF may inhibit the proliferation and migration of fibroblast-like synoviocytes and promote apoptosis by inhibiting the Wnt/β-catenin signaling pathway.

BACKGROUND:Arthroscopic anchor repair has become the main treatment method for rotator cuff tears at present.Among them,the insertion status of the anchor is a key factor in the success or failure of the operation.However,currently,the impact of the insertion depth of the anchor on the stress of the bone tunnel and the anchor under different bone density conditions remains unclear.OBJECTIVE:To explore the stress distribution of the bone tunnel and the anchor when the insertion depth of the anchor varies under different bone density conditions by using three-dimensional finite element analysis technology.METHODS:The CT image data of the humerus of volunteers were collected,and the models of the humerus and the anchor were constructed by using Mimics,3-Matic,and Solidworks software.In 3-Matic,holes with distances of 0,2,4,6,and 8 mm from the surface of the humerus were respectively created at the same position of the humerus and assembled with the anchor.In Mimics,values were assigned based on the CT gray value to obtain a model with normal bone mass(T value ≥-1.0).The parameters were changed to construct models with reduced bone mass(-2.5＜T value＜-1.0)and osteoporosis(T value＜-2.5).In each model,a 70 N pulling force was applied to the anchor along the direction tangent to the inner edge of the bone tunnel.The stress distribution and magnitude of the bone tunnel and the anchor when inserted at different depths under different bone density conditions were observed.RESULTS AND CONCLUSION:(1)When the insertion depth was the same,as the bone density decreased,the maximum equivalent stress of the anchor increased,while the maximum equivalent stress of the bone tunnel decreased.(2)When the bone density was the same,as the insertion depth of the anchor increased,the maximum equivalent stress of the anchor decreased.When the insertion depth was 4 mm,the stress of the bone tunnel was the smallest and the distribution was relatively uniform.The stress of the anchor was mainly distributed around the lower anchor hole and the proximal thread,and the stress of the bone tunnel was mainly at the part in contact with the proximal thread.The increase in the insertion depth would change the uniformity and pattern of the stress distribution,while the bone density had a relatively small impact on the stress distribution pattern.(3)It is concluded that the bone density of the humerus is crucial for the anchor repair of rotator cuff tears.It is recommended that clinicians measure the bone density of the greater tuberosity of the humerus before the operation.Excessive insertion depth of the anchor does not significantly increase its stability.Clinicians can conduct personalized preoperative assessments by using the finite element analysis method in combination with the actual situation of patients to achieve the best surgical results.

ObjectiveProtein-protein interactions (PPIs) are fundamental to the execution of biological functions within living cells. However, traditional biochemical methods, such as co-immunoprecipitation (Co-IP), often fail to capture transient, weak, or membrane-associated interactions due to the stringent detergent requirements for cell lysis. Proximity labeling (PL) has emerged in recent years as a transformative technology for mapping the proteomes of specific subcellular compartments and identifying dynamic interactomes in situ. Golgi protein 73 (GP73, also known as GOLPH2), a resident type II Golgi transmembrane protein, is a well-recognized clinical biomarker for liver diseases, including hepatocellular carcinoma (HCC). Despite its clinical significance, the comprehensive physiological and pathological functions of GP73 remain partially understood. This study aims to establish an APEX2-mediated proximity labeling system specifically targeting GP73 to map its interactome in a living cellular environment, thereby providing new insights into its molecular roles and regulatory mechanisms. MethodsTo achieve spatial specificity, we first constructed a stable cell line expressing a fusion protein consisting of GP73 and the engineered soybean peroxidase APEX2. The localization of the GP73-APEX2 fusion protein was validated to ensure it correctly targeted the Golgi apparatus. The proximity labeling reaction was initiated by incubating the cells with biotin-phenol (BP) for 30 min, followed by a brief (1 min) treatment with1 mmol/L hydrogen peroxide (H₂O₂). This catalytic reaction converts BP into highly reactive, short-lived biotin-phenoxyl radicals that covalently attach to endogenous proteins within a small labeling radius of the GP73-APEX2 enzyme. Subsequently, the cells were quenched, and biotinylated proteins were enriched using high-affinity streptavidin-coated magnetic beads. The captured “neighbor” proteins were subjected to on-bead digestion and analyzed via liquid chromatography-tandem mass spectrometry (LC-MS/MS) for high-throughput identification. Rigorous bioinformatics analysis, including Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and protein-protein interaction network mapping, was performed to interpret the biological significance of the identified candidates. ResultsOur results demonstrate the successful establishment of a robust and sensitive APEX2-based proximity labeling system for GP73. We identified a total of 95 high-confidence interacting proteins that were significantly enriched in the GP73 proximity proteome compared to control groups. Bioinformatics analysis revealed that these interactors were predominantly associated with biological processes such as vesicular transport, protein localization, and, most notably, molecular functions related to “ribosome binding” and “translation regulation”. This suggested an unexpected role for the Golgi-resident GP73 in the cellular translation machinery. To validate these findings, we performed targeted biochemical assays which confirmed a direct interaction between GP73 and the subunits of the eukaryotic translation initiation factor 3 (eIF3) complex, specifically EIF3G and EIF3I. Furthermore, functional validation using the surface sensing of translation (SUnSET) assay—a non-radioactive method to monitor protein synthesis—revealed that the overexpression of GP73 significantly promoted global protein translation levels in the cell, whereas its depletion or inhibition resulted in reduced translation efficiency. ConclusionThis study successfully utilized APEX2-mediated proximity labeling to provide the first systematic map of GP73 interactome in living cells. Our findings uncover a novel, unconventional function of GP73 as a regulator of cellular protein translation, likely mediated through its interaction with the eIF3 complex. This discovery significantly broadens our understanding of the biological roles of GP73 beyond its traditional function in the Golgi apparatus and suggests that it may act as a bridge between Golgi-related trafficking and the protein synthesis machinery. Furthermore, the technical framework established in this study provides a valuable template for investigating other complex organelle-associated protein networks and resolving transient macromolecular interactions in various physiological and pathological contexts.

ObjectiveTo investigate the association of liver fibrosis scores and inflammation markers with gallstones in patients with metabolic dysfunction-associated fatty liver disease （MAFLD）， as well as the mediating role of liver fibrosis scores in the relationship between inflammation markers and gallstones. MethodsA total of 14 567 patients who received physical examination and were diagnosed with MAFLD in Subei People’s Hospital from January 2014 to June 2023 were enrolled in this study， and according to the results of abdominal color Doppler ultrasound， they were divided into gallstone group with 1 724 patients and non-gallstone group with 12 843 patients. Related clinical data were collected from all patients， including demographic data， medical history， family history， physical examination， Color Doppler ultrasound， and biochemical parameters. The biomarkers associated with metabolic disorders and insulin resistance included triglyceride-glucose index （TyG）， TyG-body mass index （BMI） index， atherogenic index of plasma （AIP）， and non-high-density lipoprotein cholesterol-to-high-density lipoprotein cholesterol ratio （NHHR）； the biomarkers associated with inflammation and nutritional status included neutrophil-to-lymphocyte ratio （NLR）， neutrophil percentage-to-albumin ratio （NPAR）， and monocyte-to-lymphocyte ratio （MLR）； the biomarkers for assessing liver fibrosis degree and liver function included albumin-bilirubin （ALBI） score， NAFLD fibrosis score （NFS）， fibrosis-4 （FIB-4） index， and aspartate aminotransferase-to-platelet ratio index （APRI）. The independent-samples t test was used for comparison of normally distributed continuous data between two groups， while the Mann-Whitney U test was used for comparison of non-normally distributed continuous data between two groups； the chi-square test was used for comparison of categorical data between two groups. Multivariate Logistic regression analysis， restricted cubic spline analysis， and mediating effect analysis were used to assess the association of liver fibrosis markers and inflammation markers with the risk of gallstones. ResultsThe prevalence rate of gallstones was 11.8% among the MAFLD patients. There were significant differences between the gallstone group and the non-gallstone group in sex， age， smoking history， diabetes， hypertension， lymphocytes， platelets， glucose， albumin， serum uric acid， alanine aminotransferase， aspartate aminotransferase， red blood cell， NLR， NPAR， MLR， NFS， FIB-4 index， and ALBI score （all P<0.05）. The multivariate Logistic regression analysis showed that NLR （odds ratio ［OR］=1.091， 95% confidence interval ［CI］： 1.028　—　1.160， P<0.05）， NPAR （OR=1.073， 95%CI： 1.042　—　1.105， P<0.05）， MLR （OR=1.142， 95%CI： 1.057　—　1.232， P<0.05）， NFS （OR=1.239， 95%CI： 1.190　—　1.291， P<0.05）， and FIB-4 index （OR=1.326， 95%CI： 1.241　—　1.417， P<0.05） were influencing factors for the prevalence rate of gallstones. The restricted cubic spline analysis showed a significant non-linear association between NFS/FIB-4 index and the risk of gallstone （non-linear P<0.05）. The mediating effect analysis further showed that the association of NLR， MLR， and NPAR with gallstones was partially mediated by NFS or FIB-4 index， with a mediating effect accounting for 36.79%、28.09%、29.67% and 18.31%、17.70、11.57%， respectively. ConclusionNFS and FIB-4 index have a non-linear association with the prevalence rate of gallstones in MAFLD patients， and they also mediate the association of NLR， NPAR， and MLR with the risk of gallstone.