BACKGROUND:Stress-induced damage to hippocampal neurons may underlie abnormalities in neuronal structure and function,ultimately leading to mood disorders.G protein-coupled receptors in brain tissue play an important role in mood regulation.OBJECTIVE:To analyze the mechanism of depression-like behavior in chronic social defeat stress mice based on high-throughput sequencing and bioinformatics analysis.METHODS:C57BL/6J mice were randomly divided into control group and model group.There was no special treatment in the control group,while a mouse model of chronic social defeat stress was established in the model group.Depression-like behavior was assessed through the sucrose preference test,tail suspension test,and forced swim test.Anxiety behavior was evaluated using the elevated plus-maze,while social behavior was measured through the social interaction test.Cognitive function was assessed with the Y-maze spontaneous alternation test.Immunofluorescence staining was performed to quantify microglia markers in the mouse hippocampus,and Nissl staining was used to examine neuronal damage in mice.High-throughput sequencing was used to identify differentially expressed genes and gene enrichment in the mouse hippocampus,and qPCR was used to measure the expression of G protein-coupled receptors in the mouse hippocampus.RESULTS AND CONCLUSION:(1)Compared with the control group,chronic social defeat stress mice showed significant behavioral impairments,including increased anxiety,depression,and cognitive deficits.(2)Additionally,the Nissl body light density in hippocampal neurons was significantly reduced in chronic social defeat stress mice.(3)Sequencing results revealed synaptic damage in the neurons after chronic social defeat stress.Microglia activation was also markedly increased in the hippocampus of CSDS mice.Furthermore,the expression of G protein-coupled receptors in the hippocampus was significantly higher in chronic social defeat stress mice compared with the control group.These findings suggest that chronic social defeat stress induces anxiety,depression,and cognitive deficits in mice,accompanied by neuropathological changes in the hippocampus,and that altered G protein-coupled receptors expression may play a key role in these behavioral and neuropathological changes.

ObjectiveTo investigate the association between noninvasive liver fibrosis markers and carotid plaque （CP） in patients with lean metabolic dysfunction-associated fatty liver disease （MAFLD）， and to provide a basis for screening high-risk populations. MethodsA total of 957 patients with lean MAFLD who underwent physical examination in Subei People’s Hospital from January 2021 to June 2023 was enrolled as the observation cohort， with the presence or absence of CP as the outcome， and fibrosis-4 （FIB-4） index and nonalcoholic fatty liver disease fibrosis score （NFS） were used to assess liver fibrosis degree. The independent-samples t test was used for comparison of normally distributed continuous data between two groups， and 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. The multivariate logistic regression analysis， the restricted cubic spline analysis， the receiver operating characteristic curve， and the mediation effect analysis were used to investigate the association between liver fibrosis degree and CP. ResultsThe prevalence rate of CP was 36.6% in the lean MAFLD population. Compared with the non-CP group（n=607）， the CP group （n=350） had a significantly higher proportion of male patients， a significantly higher proportion of patients with smoking/diabetes/hypertension， and significantly higher levels of age， creatinine， blood urea nitrogen， triglycerides， fasting blood glucose， aspartate aminotransferase， aspartate aminotransferase/alanine aminotransferase ratio， NFS， and FIB-4 index， as well as significantly lower levels of platelet count and albumin （all P<0.05）. The multivariate logistic regression analysis showed that after adjustment for confounding factors， FIB-4 index （odds ratio［OR］=2.979， 95% confidence interval［CI］：2.141 — 4.219， P<0.001） and NFS （OR=1.747， 95%CI： 1.499 — 2.046， P<0.001） were positively correlated with CP. Both FIB-4 index and NFS had a good value in predicting CP. Hypertension had a significant indirect effect on the prevalence rate of CP through its impact on liver fibrosis markers， and its mediating effect accounted for 39.5% — 40.8% of the total effect （P<0.001）. ConclusionIn patients with lean MAFLD， NFS and FIB-4 index are significantly positively correlated with the prevalence rate of CP， and they can be used as potential epidemiological predictive indicators. Liver fibrosis markers may play a mediating role in the association between hypertension and CP. Interventions targeting hypertension and liver fibrosis markers may help to prevent and delay the progression of CP.

ObjectiveTo investigate the risk factors for acute kidney injury （AKI） after liver transplantation， and to establish and validate a risk prediction model， and to provide a basis for early identification of high-risk patients and intervention in clinical practice. MethodsA single-center retrospective study was conducted， and clinical data were collected from 162 patients who received liver transplantation in Liver Transplantation Center of The First Hospital of Shanxi Medical University from March 2020 to June 2025. The patients were divided into AKI group with 69 patients and non-AKI group with 93 patients according to the diagnostic criteria for AKI established by the Kidney Disease： Improving Global Outcomes organization and the presence or absence of AKI within 7 days after surgery. The independent-samples t test was used for comparison of normally distributed continuous data between groups， while the Mann-Whitney U test was used for comparison of non-normally distributed continuous data between groups， and the chi-square test or the Fisher’s exact test was used for comparison of categorical data between groups. The univariate differential analysis was used to obtain the factors associated with AKI after liver transplantation， and the multivariate logistic regression analysis was used to identify the independent risk factors and establish a nomogram model； the Bootstrap method with 1 000 repeated samples was used to perform internal validation of the model. The dataset was randomly divided into a training set and a validation set at a ratio of 7∶3， and the receiver operating characteristic （ROC） curve， the calibration curve， and decision curve analysis （DCA） were used to assess the discriminatory ability， calibration， and clinical applicability of the predictive model. ResultsBody mass index （BMI） （odds ratio ［OR］=1.281， 95% confidence interval ［CI］： 1.037　—　1.582， P=0.022）， serum creatinine （OR=1.097， 95%CI： 1.020　—　1.181， P=0.013）， intraoperative blood loss （OR=1.005， 95%CI： 1.002　—　1.009， P=0.004）， and cold ischemia time （OR=0.984， 95%CI： 0.976　—　0.991， P<0.001） were independent risk factors for the development of AKI after liver transplantation. The nomogram prediction model established based on the above factors had an area under the ROC curve （AUC） of 0.964 （95%CI： 0.931　—　0.997）， with an optimal cutoff value of 0.319， a sensitivity of 0.971， and a specificity of 0.903. In the training set （n=113）， the nomogram had an AUC of 0.969 （95% CI： 0.933　—　0.971）， while in the validation set （n=49）， the nomogram had an AUC of 0.941 （95%CI： 0.855　—　0.944）. The calibration curve showed good consistency between the predicted incidence rate and the actual incidence rate， and DCA showed that it had good net clinical benefit. ConclusionBMI， serum creatinine， cold ischemia time， and intraoperative blood loss are independent risk factors for the development of AKI after liver transplantation， and the nomogram prediction model established based on these factors performs well and has a good value in predicting the development of AKI after liver transplantation.

Tuberous sclerosis complex (TSC) is an autosomal dominant genetic disorder primarily caused by pathogenic variants in the TSC1 or TSC2 genes. It is characterized by multisystem hamartomatous lesions and neuropsychiatric manifestations. Here we report a young male patient with TSC involving multiple organs, caused by a heterozygous frameshift mutation in TSC2 (c.2071delC, p.Arg691Alafs^*7). The patient initially presented with classic facial angiofibromas and hypomelanotic macules, and subsequently developed psychiatric and behavioral disturbances with auditory hallucinations. Neuroimaging revealed an intracranial mass lesion, which was confirmed as a subependymal giant cell astrocytoma on postoperative histopathology following surgery performed at an outside institution. After admission, the patient underwent a comprehensive diagnostic workup, and multidisciplinary treatment evaluation revealed extensive multisystem involve-ment, including the nervous system, skin, kidneys, lungs, heart, eyes, pancreas, liver and bones. Combined with relevant literature, this article discusses the molecular mechanisms, clinical phenotypes, and comprehensive management of TSC, in order to provide a reference for the standardized and individualized management of this disease.

AIM:To assess the efficacy of intravitreal conbercept for treating neovascular age-related macular degeneration(nAMD)under a treat-and-extend(T & E)regimen.METHODS: A retrospective analysis was conducted on nAMD patients followed over a 2-year period(May 2020 to May 2022). All eyes received three monthly loading intravitreal injections of conbercept, followed by a T& E regimen in which the injection interval was extended by 2 or 4 wk according to disease activity, up to a maximum of 16 wk. When disease activity recurred, the interval was shortened. Patients were divided into initial and non-initial treatment groups based on treatment history. Best-corrected visual acuity(BCVA), central macular thickness(CMT), injection frequency, and intervals between injections over the 24-month follow-up were compared.RESULTS:Totally 27 patients(15 males and 12 females, 33 eyes)were enrolled. In the initial treatment group(18 eyes, mean age 65.72±12.32 y), BCVA significantly improved at 1, 3, and 6 mo(P<0.05), and CMT significantly improved at 1 and 3 mo(P<0.05). In the non-initial treatment group(15 eyes, mean age 69.00±9.21 y), BCVA improved significantly at 3 mo(P<0.05), whereas CMT remained stable(P >0.05). Baseline CMT was similar between the groups(P>0.05). However, significant differences were observed at multiple post-injection time points(P<0.05). The total number of injections did not differ between the groups(P>0.05). Intervals between injections varied, with the majority at 4 and 3-4 mo in the initial and non-initial treatment groups, respectively.CONCLUSION:Initiating intravitreal conbercept therapy under a T & E regimen results in superior visual and anatomical outcomes compared with non-initial treatment.

ObjectiveTo investigate the therapeutic mechanism of Xuefu Zhuyutang （XFZYT） for metabolic-associated fatty liver disease （MAFLD） through integrated network pharmacology and animal experiments. MethodsNetwork pharmacology was utilized to predict the core components， key therapeutic targets， and signaling pathways of XFZYT in the treatment of MAFLD. For animal experiments， a rat model of MAFLD was established by feeding a high-cholesterol diet for 4 weeks. Intervention was then administered with low-dose （2 g·kg^-1） and high-dose （4 g·kg^-1） XFZYT for 2 weeks. Biochemical assays were performed to measure the serum levels of aspartate aminotransferase （AST）， alanine aminotransferase （ALT）， total cholesterol （TC）， triglycerides （TG）， high-density lipoprotein （HDL）， and low-density lipoprotein （LDL）. In addition， the activities of superoxide dismutase （SOD） and catalase （CAT） and levels of malondialdehyde （MDA） and glutathione （GSH） in the serum were measured. The same way was adopted to measure the levels of TC and TG in the liver tissue. Enzyme-linked immunosorbent assay （ELISA） was employed to quantify the serum levels of interleukin （IL）-6， IL-1β， and tumor necrosis factor-alpha （TNF-α）. Histopathological evaluations included hematoxylin and eosin （HE） staining for liver tissue morphology， Oil Red O staining for lipid deposition， and dihydroethidium （DHE） probe staining for reactive oxygen species （ROS） levels. Western blot analysis was conducted to assess the protein levels of AMP-activated protein kinase （AMPK）， phosphorylated （p）-AMPK， nuclear factor erythroid 2-related factor 2 （Nrf2）， heme oxygenase-1 （HO-1）， nuclear factor-kappa B （NF-κB）， and p-NF-κB in the liver tissue. Untargeted metabolomics analysis of the serum was performed by liquid chromatography-tandem mass spectrometry （LC-MS/MS）. ResultsNetwork pharmacology analysis predicted 155 potential targets of XFZYT for MAFLD treatment， with core targets including signal transducer and activator of transcription 3 （STAT3）， protein kinase B1 （Akt1）， TNF， and IL-6. Kyoto Encyclopedia of Genes and Genomes （KEGG） pathway enrichment primarily implicated the AMPK signaling pathway. Animal experiments demonstrated that compared with the normal group， the model group exhibited dyslipidemia， hepatic function impairment， pronounced hepatic lipid deposition， and inflammatory manifestations， with elevated serum levels of AST， ALT， TC， TG， LDL， and MDA （P<0.05）， reduced HDL and GSH levels plus decreased SOD and CAT activities （P<0.05）， downregulated protein levels of Nrf2， HO-1， and p-AMPK （P<0.05）， and upregulated protein level of p-NF-κB （P<0.05） in the liver tissue. Compared with the model group， XFZYT intervention groups showed significant amelioration of dyslipidemia and hepatic function impairment， markedly reduced hepatic lipid deposition and inflammatory cell infiltration， decreased serum levels of AST， ALT， TC， TG， LDL， and MDA （P<0.05）， increased HDL and GSH levels plus enhanced SOD and CAT activities （P<0.05）， upregulated protein levels of Nrf2， HO-1， and p-AMPK （P<0.05）， and downregulated protein level of p-NF-κB （P<0.05）. Serum metabolomics revealed 511 differentially expressed metabolites （231 upregulated and 280 downregulated） between normal and model groups， while XFZYT groups versus model group showed 94 differential metabolites （51 upregulated and 43 downregulated）. Among them， 11 metabolites displayed the most significant alterations， with enriched pathways including glycerolipid metabolism， cholesterol metabolism， and insulin resistance， multiple of which demonstrated AMPK association. ConclusionXFZYT alleviates MAFLD by regulating the AMPK signaling pathway and associated metabolic networks.

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.

ObjectiveTo explore the effect of oral sodium butyrate on skeletal muscle atrophy in CT26 tumor mice through the gut microbiota-skeletal muscle axis and its potential mechanism. MethodsSixty SPF BALB/c male mice aged 8 weeks were randomly divided into a normal control group (NC, n=18) and a ABX-depleted group (ABX, n=42). The ABX mice were pretreated with a quadruple antibiotic cocktail via oral gavage (0.2 ml per administration, once daily, 6 d per week, for 2 weeks), whereas NC received an equal volume of sterile water. The quadruple antibiotic cocktail consisted of metronidazole (1 g/L), vancomycin (0.5 g/L), ampicillin (1 g/L), and gentamicin (1 g/L). Following successful pretreatment, six mice from each group were randomly selected for gut microbiota sequencing analysis and designated as the Abx group and the NC0 group, respectively. Theremaining mice in ABX were subcutaneously inoculated in the dorsum with 0.2 ml of CT26 cell suspension (at a cell density of 1×10⁷/ml). Then these mice were randomly allocated into three subgroups: a control tumor bearing model group (0_NaB, n=12), a tumor-bearing model group receiving low-dose oral sodium butyrate (L_NaB, n=12), a tumor-bearing model group receiving high-dose oral sodium butyrate (H_NaB, n=12). And mice in NC were inoculated at the same site with 0.2 ml of normal saline. The administration dose for L_NaB was 0.3 g/(kg·d), that for H_NaB was 0.5 g/(kg·d), while NC and 0_NaB were given the same volume of normal saline (0.2ml per time, once daily, 6 d per week, for 4 weeks). The general condition of mice was monitored, and forelimb grip strength gastrocnemius muscle mass and its muscle fiber cross-sectional area were measured for each group. The structural changes in gut microbiota were assessed by 16S rRNA sequencing of cecal contents. Pathological alterations in the intestinal wall were examined via HE staining. Serum and gastrocnemius muscle levels of TNF‑α, IL-6, IL-1β, and LPS were quantified using ELISA. The protein expression of ZO-1 and occludin in the small intestine, as well as proteins associated with the TLR4/MyD88/NF-κB signaling pathway in the gastrocnemius muscle, were detected by Western blot analysis. Results(1) The alpha-diversity in Abx was significantly lower than that in NC0 (P<0.01), a significant decrease of the mass and muscle fiber cross-sectional area of the gastrocnemius (P<0.01), with the majority of gut microbiota being effectively depleted. (2) Compared with NC, the subcutaneous tumors of mice in 0_NaB were prominent, a significant increase of the mass and muscle fiber cross-sectional area of the gastrocnemius, accompanied by a significant decrease in body weight at the end of the 3th and 4th week (P<0.05), and a significant weakening of the forelimb grasping strength at the 5th and 6th week (P<0.01). Compared with 0_NaB, the tumor mass of mice in L_NaB and H_NaB showed a significant decreasing trend, and the grip strength of the forelimbs significantly increased at the 5th and 6th week (P<0.05, P<0.01). (3) Compared with 0_NaB, the Shannon and Observed species indices in α diversity of L_NaB and H_NaB were significantly increased (P<0.05). At the genus level, compared with 0_NaB, L_NaB exhibited a significant decrease in the relative abundance of Parasutterella (P< 0.01), while H_NaB showed significant reductions in the relative abundances of both Escherichia-Shigella and Parasutterella (P < 0.01). (4) Compared with 0_NaB, the small intestinal tissue structure in L_NaB and H_NaB was more intact, the infiltration of inflammatory cells was significantly reduced, and the capillaries were slightly dilated. The expression levels of ZO-1 and occludin proteins in L_NaB were significantly increased (P<0.01). (5) The LPS concentration in the gastrocnemius muscle and the protein expression levels of TLR4, MyD88, p-IκBα, and p-NF‑κB p65 in L_NaB and H_NaB were significantly lower than those in 0_NaB (P<0.05). The serum TNF‑α concentration in H_NaB and TNF-α concentration in the gastrocnemius muscle of the L_NaB and H_NaB were significantly lower than those in 0_NaB (P<0.05, P<0.01, P<0.01). ConclusionOral administration of NaB can improve gut microbiota α diversity, adjusting its composition, improving intestinal mucosal barrier function, reducing the LPS-induced pro-inflammatory response, and delaying skeletal muscle atrophy. The underlying mechanism may involve down regulation of TLR4/MyD88/NF-κB signaling in skeletal muscle.

ObjectivePhotoacoustic tomography (PAT) holds significant potential for high-resolution deep-tissue imaging. In preclinical research, custom-designed concave arc-shaped ultrasound transducer arrays are often used to maximize the detection aperture. However, manufacturing limitations and assembly tolerances frequently cause the actual physical positions of array elements to deviate from their theoretical design. Additionally, concave arrays are typically covered with an acoustic lens, which introduces a mismatch in the speed of sound between the coupling medium and the lens material. The combination of these geometric and acoustic-phase errors leads to severe image artifacts, reduced contrast, and degraded resolution. This study proposes a systematic two-step calibration strategy to address these issues and substantially improve image quality. MethodsFirst, a high-intensity isotropic photoacoustic point source was constructed using a multi-mode optical fiber coated with carbon nanotubes (CNTs) to acquire high signal-to-noise ratio calibration data. The Akaike information criterion (AIC) was employed to accurately determine the time of arrival (ToA) of photoacoustic signals. Subsequently, a geometric calibration algorithm based on nonlinear least-squares (NLS) estimation was developed. This algorithm iteratively solves for the true spatial coordinates of each array element by minimizing the residual between theoretical and measured acoustic path lengths. To further address sound-speed inhomogeneity caused by the acoustic lens, a phase compensation algorithm based on bilinear interpolation was proposed. This algorithm computes a pixel-specific phase delay map across the imaging region and performs point-by-point signal correction during delay-and-sum (DAS) reconstruction. The proposed methods were validated using a custom 96-channel concave arc-shaped array (center frequency: 12  MHz) through both phantom imaging and in vivo mouse tumor models. ResultsPhantom experiments showed that at an imaging depth of14 mm, the reconstruction position deviation of the point source in the uncalibrated system reached up to 1  mm. After applying the combined calibration, the lateral resolution (full width at half maximum, FWHM) at the focal point of the arc array reached 95  μm—representing a 85% reduction compared to the uncalibrated state and a 79% reduction compared to geometric calibration alone without phase compensation. In vivo experiments demonstrated that the calibrated system clearly resolved the microvascular network of subcutaneous tumors in mice. Photoacoustic signals were strictly confined within tumor boundaries delineated by ultrasound imaging (USI), eliminating the vascular spillover artifacts commonly observed in uncalibrated images. Furthermore, after intravenous injection of indocyanine green (ICG), the system successfully detected weak photoacoustic signals at a depth of 5 mm, performing significantly better than the uncalibrated system. ConclusionThe proposed calibration method, which integrates nonlinear least-squares estimation with phase compensation, significantly improves image fidelity and spatial resolution consistency across a wide field of view by correcting systemic geometric errors and acoustic phase aberrations. This approach demonstrates high robustness and provides a reliable technical foundation for the clinical translation of photoacoustic probes with non-standard geometries.

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.