Surface proteins play pivotal roles in physiological processes, including cell recognition, signal transduction, substance transport, and immune responses. However, challenges persist in characterizing abnormal surface proteins in disease states and identifying therapeutic targets, due to the low abundance of these proteins within the total proteome and the frequent presence of their complex glycosylation modifications. Recent years have witnessed the vigorous development of chemical proteomics, leading to the successful creation of various chemical probes for the labeling and characterization of cell surface proteins. These techniques have subsequently been applied to the detection of disease surface proteins and the discovery of corresponding targets. Surface protein characterization techniques based on chemical proteomics are discussed herein, focusing on the principles of amino acid-targeted labeling, proximity labeling, and glycoprotein capture. The novelty, advantages, and limitations of techniques such as targeted lysine labeling, peroxidase and photocatalytic proximity labeling, and chemical glycan capture and metabolic glycan labeling are elaborated, and their applications across various biological models and disease types are described, aiming to provide some reference for target discovery and drug development targeting surface proteins.

Autism spectrum disorder （ASD）， also known as autism， is a series of neurodevelopmental disorders characterized by social disorders and repetitive stereotyped behaviors/narrow interests. Its pathogenesis is complex， and there is a lack of effective treatment drugs， with some cases having adverse outcomes. Recent studies have consistently revealed that individuals with autism spectrum disorder （ASD） commonly exhibit characteristics such as gut microbiota dysbiosis （abnormal Bacteroidetes/Firmicutes ratio）， impaired intestinal barrier function （elevated serum levels of zonulin and LPS）， and intestinal immune dysregulation （increased pro-inflammatory cytokines including IL-6 and TNF-α）， suggesting that gastrointestinal abnormalities may influence central nervous system development through neuroendocrine， immunoregulatory， and metabolic pathways. Consequently， growing scholarly attention has focused on dietary interventions as potential approaches to alleviate clinical symptoms in children with ASD. This review systematically summarizes the role of gut microbiota and their metabolite alterations in ASD pathogenesis， along with recent advancements in understanding the microbiota-gut-brain axis mechanisms. Additionally， it elaborates on the therapeutic effects and underlying biological basis of restrictive diet therapy， modified diet therapy， and nutritional supplementation therapy in promoting the health of children with ASD. This systematic review reveals that children with ASD exhibit significant gut microbiota dysbiosis （e.g.， increased Clostridium， decreased Faecalibacterium） and abnormal metabolite profiles （e.g.， altered short-chain fatty acid spectra， elevated 4EPS levels）. These alterations exacerbate neuroinflammation and immune dysregulation through the microbiota-gut-brain axis， thereby impacting nervous system development and function. Furthermore， interventions such as ketogenic diets， camel milk， and specific nutritional supplements can alleviate certain ASD symptoms by modulating gut microbiota， restoring intestinal barrier function， and improving metabolic pathways. Future investigations should aim to create multi-omics evaluation systems for pinpointing potential beneficiaries， devise individualized intervention strategies rooted in microbiome characteristics， and verify their therapeutic value and safety in large-scale randomized controlled trials. These efforts are crucial to transitioning ASD treatment from symptomatic control to address disease etiology， thereby paving the way for improving prognoses.

Objective Objective To analyze the potential spatial factors affecting the genetic differentiation of Oncomelania hupensis robertsoni in Yunnan Province. Methods A total of 13 administrative villages were selected from schistosomiasis-endemic areas of Yunnan Province as O. hupensis snail sampling sites. At least 200 snails were collected in each site, and the spatial variable data of each site were recorded, including longitude, latitude and altitude. Thirty active and Schistosoma japonicum uninfected O. hupensis snails were selected from each sampling site by means of the crawling method and the cercarial shedding method. Genomic DNA was extracted from O. hupensis snails. Following PCR amplification, purification of PCR amplification products and sequencing, the gene sequences of O. hupensis snail samples were spliced and edited using the DNAstar software and the NCBI database to yield the complete mitochondrial sequences of O. hupensis snails at each sampling site, and the mitochondrial genetic distance matrix of O. hupensis robertsoni was calculated at each sampling site. The geographical coordinates of each sampling site were marked using the software ArcGIS 10.2, and the straight-line geographical distance between each sampling site was calculated. The altitude difference, longitude difference and latitude difference between each sampling site were calculated using the Excel software, and the correlation between the mitochondrial genetic distance matrix of O. hupensis robertsoni and each spatial variable matrix was examined by using the Mantel test at 13 sampling sites in Yunnan Province. Results Among the 13 O. hupensis snail sampling sites in Yunnan Province, the largest mitochondrial genetic distance of O. hupensis robertsoni snail populations was seen between Anding Village, Nanjian Yi Autonomous County and Caizhuang Village, Midu County (26.244 2), and the largest geographical distance was seen between Dongyuan Village, Gucheng District and Cangling Village, Chuxiong County (272.64 km). The highest altitude difference was seen between Anding Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1 086.10 m), and the largest longitude difference was found between Qiandian Village, Eryuan County and Cangling Village, Chuxiong County (1.86°), while the largest latitude difference was measured between Leqiu Village, Nanjian Yi Autonomous County and Dongyuan Village, Gucheng District (1.81°). In addition, the mitochondrial genetic distance of O. hupensis robertsoni snail populations was positively correlated with altitude at 13 snail sampling sites in Yunnan Province (r = 0.542 8, P < 0.001), and showed no significant correlations with geographical distance (r = 0.093 4, P > 0.05), longitude (r = −0.199 5, P > 0.05) or latitude (r = 0.205 7, P > 0.05). Conclusion Altitude may be a potential spatial factor affecting the genetic differentiation of O. hupensis robertsoni in Yunnan Province.

ObjectiveTo analyze the characteristics of 150 patients with spinal cord injury complicated with spasticity. MethodsA cross-sectional survey was conducted on 150 patients with spinal cord injury accompanied by spasticity from September, 2019 to December, 2024. Their age, gender, cause of injury, injury site, severity of injury, spasticity severity and other indicators were recorded. The relationships between different characteristics were analyzed, and a correlation analysis of disease duration, spasticity grade, injury level, injury severity and age were conducted. ResultsThere was no significant difference in age distribution between patients with tetraplegia and paraplegia (Z = 0.806, P = 0.420). The proportions of trauma (χ² = 3.982, P = 0.046) and tetraplegia (χ² = 10.559, P = 0.010) were higher in males than in females. Trauma was the main cause of injury in both tetraplegia and paraplegia patients; the proportion of tetraplegia was higher than paraplegia in trauma patients, while paraplegia was higher than tetraplegia in non-trauma patients (χ² = 11.885, P < 0.001). Patients with tetraplegia was dominated by incomplete injury, whereas patients with paraplegia was dominated by complete injury (χ² = 10.885, P = 0.012). Grade A injury was predominant in trauma patients (P = 0.003). Spasticity grade showed a very weak positive correlation with disease duration (r = 0.175, P = 0.032) and age (r = 0.168, P = 0.040). Injury severity showed a very weak positive correlation with age (r = 0.183, P = 0.025). ConclusionCharacteristics of patients with spinal cord injury complicated with spasticity is different with gender, cause of injury, injury level, injury severity.

ObjectiveThis Study was conducted to investigate the interaction mechemisms between gutmicrobiota dysregulation and ischemic stroke by establishing a rat model of ischemic stroke and employing fecal microbiota transplantation (FMT). MethodsA preliminary experiment was conducted to establish an antibiotic-induced pseudo-sterile (ABX) rat model through antibiotic treatment, and a cerebral ischemia model was prepared using the middle cerebral artery occlusion (MCAO) method. Fecal microbiota from stroke patients and healthy individuals were transplanted via FMT, followed by behavioral testing. 16S rRNA sequencing was used to analyze the microbial community, hematoxylin and eosin (HE) staining to observe histopathological status, transmission electron microscopy (TEM) to examine the tight junction structure of the small intestine, and enzyme-linked immunosorbent assay (ELISA) to detect levels of inflammatory factors and intestinal barrier-related markers. Results16S rRNA sequencing of fecal samples showed that compared with the normal control group and the metronidazole group, the abundance and diversity of fecal microorganisms in the quadruple antibiotic group were significantly reduced, indicating successful establishment of the ABX model. After transplanting fecal microbiota from stroke patients into ABX rats, significant changes in gut microbiota composition were observed. Behavioral tests revealed that the MCAO model group showed significant decreases in both horizontal movement and vertical exploration abilities. ELISA results indicated that IL-17 concentration in the ABX+mFMT (antibiotic-treated+model fecal microbiota transplantation) group was lower than in the ABX+cFMT (antibiotic-treated+control fecal microbiota transplantation) group, suggesting that IL-17 may serve as a key inflammatory indicator for evaluating the impact of stroke intervention on gut microbiota. Triphenyltetrazolium chloricle staining (TTC) staining suggested that gut microbiota intervention may increase the risk of stroke. HE staining showed that, except for the control group, all groups exhibited ischemic changes and inflammatory infiltration in brain tissues. TEM revealed that microvilli of small intestinal epithelial cells in the ABX+mFMT group were sparser than those in the ABX+cFMT group, indicating that microbial intervention affects intestinal barrier function. ConclusionThe ABX model established using broad-spectrum antibiotics showed no significant differences in physiological characteristics compared to normal rats, and the findings were consistent with those from germ-free rat models. Stroke prognosis appears to be influenced by intestinal dysbiosis, accompanied by significantly elevated levels of the pro-inflammatory cytokine IL-17, which may exacerbate neural injury via the gut-brain axis. Behavioral experiments indicated that transplantation of gut microbiota from stroke rats impaired cognitive function. Furthermore, IL-17 demonstrated sensitivity to alterations in the gut microbiota, suggesting its potential as a key therapeutic target for stroke intervention.

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.

ObjectiveThis Study was conducted to investigate the interaction mechemisms between gutmicrobiota dysregulation and ischemic stroke by establishing a rat model of ischemic stroke and employing fecal microbiota transplantation (FMT). MethodsA preliminary experiment was conducted to establish an antibiotic-induced pseudo-sterile (ABX) rat model through antibiotic treatment, and a cerebral ischemia model was prepared using the middle cerebral artery occlusion (MCAO) method. Fecal microbiota from stroke patients and healthy individuals were transplanted via FMT, followed by behavioral testing. 16S rRNA sequencing was used to analyze the microbial community, hematoxylin and eosin (HE) staining to observe histopathological status, transmission electron microscopy (TEM) to examine the tight junction structure of the small intestine, and enzyme-linked immunosorbent assay (ELISA) to detect levels of inflammatory factors and intestinal barrier-related markers. Results16S rRNA sequencing of fecal samples showed that compared with the normal control group and the metronidazole group, the abundance and diversity of fecal microorganisms in the quadruple antibiotic group were significantly reduced, indicating successful establishment of the ABX model. After transplanting fecal microbiota from stroke patients into ABX rats, significant changes in gut microbiota composition were observed. Behavioral tests revealed that the MCAO model group showed significant decreases in both horizontal movement and vertical exploration abilities. ELISA results indicated that IL-17 concentration in the ABX+mFMT (antibiotic-treated+model fecal microbiota transplantation) group was lower than in the ABX+cFMT (antibiotic-treated+control fecal microbiota transplantation) group, suggesting that IL-17 may serve as a key inflammatory indicator for evaluating the impact of stroke intervention on gut microbiota. Triphenyltetrazolium chloricle staining (TTC) staining suggested that gut microbiota intervention may increase the risk of stroke. HE staining showed that, except for the control group, all groups exhibited ischemic changes and inflammatory infiltration in brain tissues. TEM revealed that microvilli of small intestinal epithelial cells in the ABX+mFMT group were sparser than those in the ABX+cFMT group, indicating that microbial intervention affects intestinal barrier function. ConclusionThe ABX model established using broad-spectrum antibiotics showed no significant differences in physiological characteristics compared to normal rats, and the findings were consistent with those from germ-free rat models. Stroke prognosis appears to be influenced by intestinal dysbiosis, accompanied by significantly elevated levels of the pro-inflammatory cytokine IL-17, which may exacerbate neural injury via the gut-brain axis. Behavioral experiments indicated that transplantation of gut microbiota from stroke rats impaired cognitive function. Furthermore, IL-17 demonstrated sensitivity to alterations in the gut microbiota, suggesting its potential as a key therapeutic target for stroke intervention.

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 crystalline lens of the eye is recognized as one of the most radiosensitive tissues in the human body. While the International Commission on Radiological Protection (ICRP) has classified ionizing radiation (IR)-induced cataracts as a tissue reaction (deterministic effect) and subsequently reduced the occupational equivalent dose limit for the lens, significant uncertainties remain regarding the precise dose threshold and the complex biological pathways driving lens opacification. This review provides a comprehensive synthesis of current knowledge concerning radiation-induced lens damage, integrating epidemiological exposure characteristics with dose-response modeling and mechanistic molecular insights. First, we analyze exposure characteristics through four epidemiological dimensions: dose, time, space, and population. Clinical evidence suggests that radiation cataracts—particularly posterior subcapsular opacities—exhibit a distinct latency period that is inversely correlated with dose. We highlight that risk is not confined to acute high-dose scenarios (such as in atomic bomb survivors) but is increasingly relevant in chronic low-dose occupational settings (e.g., interventional radiology) and medical diagnostics (e.g., CT scans). Crucially, individual susceptibility is modified by genetic background, age, and environmental co-factors, complicating risk assessment. Second, we critically examine the dose-effect relationship. Although the ICRP suggests a threshold of 0.5 Gy, emerging data challenge the traditional threshold model, with some studies advocating for a linear non-threshold (LNT) relationship. We further discuss the critical roles of radiation quality and dose rate. High linear energy transfer (LET) radiation demonstrates a significantly higher relative biological effectiveness (RBE) for cataractogenesis compared to low-LET radiation. Paradoxically, and unlike many other tissues, the lens may exhibit an “inverse dose-rate effect,” where fractionated or protracted exposures potentially enhance biological damage—a finding that challenges classical radiobiological paradigms. Third, drawing upon the “cataractogenic load” hypothesis and the unique physiological constraints of the lens, this review elucidates the multidimensional molecular mechanisms driving radiation-induced opacification. Key mechanisms include four aspects. (1) DNA damage and repair: IR induces DNA double-strand breaks (DSBs) that, due to the lens’ limited repair capacity (modulated by genes such as ATM, Ptch1, and Ercc2), lead to the accumulation of damage. (2) Antioxidant defense system: dysfunction of the Nrf2/HO-1 antioxidant axis results in redox imbalances, triggering NF-κB-mediated inflammation and protein aggregation. (3) Cell proliferation and senescence: IR disrupts cell cycle regulation, causing a dichotomy of effects—driving premature senescence in some cell populations (evidenced by ATM nuclear foci) while inducing aberrant proliferation via growth factor upregulation (FGF2, TGFβ) in others. (4) Cell migration and adhesion: activation of the Wnt/β‑catenin pathway and alterations in the E-cadherin complex promote the abnormal migration of epithelial cells to the posterior capsule, a hallmark of radiation-induced cataracts. In conclusion, radiation-induced cataractogenesis is a multifactorial process in which genetic susceptibility and environmental stressors converge to overwhelm the lens’ homeostatic thresholds. Future research must prioritize longitudinal cohort studies to refine dose thresholds and employ multi-omics approaches to map the crosstalk between DNA damage responses and matrix remodeling. Establishing a robust mechanistic model is essential for developing targeted radioprotective strategies and optimizing radiation protection standards for occupational and medical safety.

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.