ObjectiveAge-related macular degeneration （AMD） is one of the leading causes of low vision and blindness in people over 50 years old， and dry AMD （dAMD） is one type for which there is currently no clear treatment. On the basis of the diagnosis and clinical characteristics of dAMD in traditional Chinese and Western medicine， this paper evaluated the fitting degrees of existing animal models of dAMD with clinical characteristics according to the evaluation methods of animal models， and put forward suggestions and prospects. MethodsLiterature on animal models of dAMD was searched against database， and the characteristics of the models were assigned according to the diagnosis criteria of diseases and syndromes of traditional Chinese and Western medicine， and the fitting degrees of the models with clinical characteristics were analyzed and evaluated. ResultsAt present， the animal models of dAMD are mainly established targeting complement factors， chemokines， oxidative damage， lipid/glucose metabolism， and natural strains. Most of the models can simulate the major pathological changes of dAMD， showing the fitting degree of 25%-50% with clinical characteristics in Western medicine. However， the evaluation of traditional Chinese medicine （TCM） syndromes， especially the evaluation of secondary syndromes， is missing， and the models present low fitting degrees with the clinical characteristics in TCM. ConclusionExisting animal models of dAMD are mostly established under the guidance of Western diagnostic standards， which reproduce the main disease characteristics of Western medicine and lack observation of TCM syndromes. Future studies can pay attention to the intervention factors and evaluation systems of spleen deficiency Qi deficiency and liver-kidney Yin deficiency syndrome and build the animal model of dAMD with integration of disease and syndrome based on clinical characteristics of traditional Chinese and Western medicine.

ObjectiveAge-related macular degeneration （AMD） is one of the leading causes of low vision and blindness in people over 50 years old， and dry AMD （dAMD） is one type for which there is currently no clear treatment. On the basis of the diagnosis and clinical characteristics of dAMD in traditional Chinese and Western medicine， this paper evaluated the fitting degrees of existing animal models of dAMD with clinical characteristics according to the evaluation methods of animal models， and put forward suggestions and prospects. MethodsLiterature on animal models of dAMD was searched against database， and the characteristics of the models were assigned according to the diagnosis criteria of diseases and syndromes of traditional Chinese and Western medicine， and the fitting degrees of the models with clinical characteristics were analyzed and evaluated. ResultsAt present， the animal models of dAMD are mainly established targeting complement factors， chemokines， oxidative damage， lipid/glucose metabolism， and natural strains. Most of the models can simulate the major pathological changes of dAMD， showing the fitting degree of 25%-50% with clinical characteristics in Western medicine. However， the evaluation of traditional Chinese medicine （TCM） syndromes， especially the evaluation of secondary syndromes， is missing， and the models present low fitting degrees with the clinical characteristics in TCM. ConclusionExisting animal models of dAMD are mostly established under the guidance of Western diagnostic standards， which reproduce the main disease characteristics of Western medicine and lack observation of TCM syndromes. Future studies can pay attention to the intervention factors and evaluation systems of spleen deficiency Qi deficiency and liver-kidney Yin deficiency syndrome and build the animal model of dAMD with integration of disease and syndrome based on clinical characteristics of traditional Chinese and Western medicine.

Objective: To understand the impact of childhood trauma on psychological distress among upper grade elemetary school students, and to explore the mediating role of leisure activities in the relationship, so as to provide a basis for developing mental health intervention strategies. Methods: From August to November 2024, a combination of convenience sampling and stratified cluster random sampling was employed to recruit 1 373 fourth to sixth grade students from four primary schools in Harbin. The Childhood Trauma Questionnaire(CTQ), a self designed leisure activity scale (including active and passive leisure activities), and the Kessler Psychological Distress Scale (K10) were used to assess childhood trauma experiences, leisure activities, and levels of psychological distress. Spearman correlation analysis and linear regression analysis were conducted to explore the relationships among childhood trauma, leisure types, leisure time, and psychological distress. Based on the mediation analysis framework proposed by Hayes (Model 4), the mediating role of leisure types in the relationship between childhood trauma and psychological distress was examined. Results: Totally 19.1% of the upper elemetary school students exhibited psychological distress, while 30.2% had experienced childhood trauma. During school days, 64.6% of the students were reported of having leisure time concentrated between 1 and 5 hours per day, whereas 67.4% reported leisure time exceeding 5 hours per day on weekends. After controlling for potential demographic confounders such as gender, grade, ethnicity, household registration, being an only child, parents educational level, co residence, and whether parents are first time married,linear regression analysis showed that childhood trauma experience had positive predictive effect on psychological distress in upper primary school students( β =0.20, P <0.01). Leisure time showed no statistically significant association with psychological distress, both on school days ( β =-0.58 to -0.56) and weekends ( β =0.26- 0.98 )(all P >0.05). Active leisure activities were negatively associated with psychological distress ( β =-0.20), while passive leisure activities were positively associated with psychological distress ( β =0.29)(both P <0.01). Leisure type partially mediated the relationship between childhood trauma and psychological distress, accounting for 11.7% of the indirect effect. Conclusion Childhood trauma experiences positively predict psychological distress in upper elementary school students, and affect psychological distress through active leisure and passive leisure.

OBJECTIVE: To explore the genetic etiology of 46 Chinese pedigrees affected with Hereditary multiple exostoses (HME) and provide genetic counseling and reproductive intervention. METHODS: Whole-exome sequencing and Sanger sequencing were carried out on 87 patients from the 46 pedigrees to analyze the variants of EXT1 and EXT2 genes. Pathogenicity of the variants was assessed based on the guidelines from the American College of Medical Genetics and Genomics and Association for Molecular Pathology (ACMG/AMP). Prenatal diagnosis and preimplantation genetic testing (PGT) were provided for couples with identified pathogenic mutations. This study was approved by the Medical Ethics Committee of the hospital (Ethics No.: LL-SC-SG-2014-010). RESULTS: In total 17 and 22 pathogenic variants were respectively identified in the EXT1 and EXT2 genes, among which 5 EXT1 and 12 EXT2 variants were unreported previously. Three patients with no family history were found to harbor de novo variants of the EXT1 gene. Twenty nine couples had opted for PGT or underwent prenatal diagnosis following natural conception, and 17 healthy babies were born. CONCLUSION This study has clarified the genetic etiology of 45 HME pedigrees and identified 17 novel variants, which has enriched the mutational spectrum of the EXT1 and EXT2 genes. Reproductive intervention through PGT and prenatal diagnosis have prevented the recurrence of HME in these families.
Humans ; Female ; Male ; Pedigree ; Exostoses, Multiple Hereditary/diagnosis* ; N-Acetylglucosaminyltransferases/genetics* ; Adult ; Exostosin 1 ; Asian People/genetics* ; Genetic Testing ; Exostosin 2 ; Mutation ; China ; Prenatal Diagnosis ; Pregnancy ; Genetic Counseling ; Preimplantation Diagnosis ; Exome Sequencing ; East Asian People

ObjectiveTo clarify the expression of TRIM28 in non-M3 acute myeloid leukemia （AML） and its correlation with clinical indicators and prognosis， and to further explore the effect of TRIM28 expression levels on the proliferation and apoptosis of AML cells using small interfering RNA. MethodsThe GSE34577 dataset was analyzed using R software to compare TRIM28 expression between healthy controls and non-M3 acute myeloid leukemia （AML） patients. Clinical samples from non-M3 AML patients were collected， with TRIM28 expression levels measured using real-time quantitative PCR （qPCR）. The analysis focused on correlations between TRIM28 expression and various clinical indicators， treatment efficacy， and patient prognosis. Furthermore， small interfering RNA （siRNA） technology was employed to downregulate TRIM28 expression in human primary AML cells （HL60 cell line）. The effects on cell proliferation and apoptosis were then assessed through CCK-8 assays and flow cytometry， respectively. ResultsThe results showed that TRIM28 was up-regulated in non-M3 AML of both online database GSE34577 and clinical samples （P<0.000 1）， TRIM28 expression of new diagnosis group and relapsed refractory group was higher than iron deficiency anemia group （P<0.01）， and there was no significance between different French-American-British classification systems subtype. TRIM28 expression was higher in non-M3 AML patients with a poor genetic prognosis stratified as moderate than in the good prognosis group， and TRIM28 expression was associated with NPM1 combined with the FLT3-ITD mutation， positively correlated with age， bone marrow blast， peripheral blood blast and white blood cell， negatively correlated with hemoglobin. In addition， interference TRIM28 greatly inhibited cell proliferation and promoted cell apoptosis. ConclusionThis study reveals that TRIM28 is highly expressed in non-M3 AML and associated with prognosis， and plays a key role in the proliferation and apoptosis of AML cells， suggesting that TRIM28 may serve as a novel therapeutic target for non-M3 AML.

ObjectiveTo investigate the incidence and influencing factors of posttraumatic stress disorder （PTSD） three months after a traffic accident， and to explore the role of social support and coping strategies. MethodsA total of 117 individuals exposed to trauma following road traffic accidents were recruited. General demographic and clinical information was collected within one week， and the hamilton anxiety rating scale （HAMA）， the hamilton depression rating scale-24 （HAMD-24）， the social support rating scale （SSRS）， and the simplified coping style questionnaire （SCSQ） were administered. A 3-month follow-up was subsequently conducted， during which PTSD symptoms were assessed using the post-traumatic stress disorder checklist for DSM-5 （PCL-5）. Participants were divided into a PTSD group and a non-PTSD group according to whether PTSD occurred. Between-group comparisons were performed using the Mann-Whitney U non-parametric test or the χ² test， as appropriate. Spearman correlation analysis was used to examine the associations between general characteristics and PCL-5 scores. Binary Logistic regression was applied to identify factors influencing PTSD， and receiver operating characteristic （ROC） curve analysis was conducted to evaluate the diagnostic value of the SCSQ and SSRS. ResultsDuring the 3-month follow-up of the 117 trauma-exposed individuals， 17 cases developed PTSD， with a higher proportion of females （70.59%）. Between-group comparisons showed that， compared with the PTSD group， the non-PTSD group had higher scores for positive coping， objective support， and subjective support （P<0.05）， and lower scores for negative coping， HAMA， HAMD， and PCL-5 （P<0.05）. Correlation analysis indicated that female gender， negative coping， and higher HAMA and HAMD scores were associated with greater PTSD severity. Logistic regression analysis demonstrated that educational level （OR=1.715， 95% CI： 1.020-2.883， P=0.042） and negative coping （OR=1.590， 95% CI： 1.003-2.522， P=0.048） were risk factors for PTSD， whereas objective support （OR=0.646， 95% CI： 0.451-0.925， P=0.017） was a protective factor. The ROC analysis showed that the total SCSQ score and its negative and positive coping dimensions， the total SSRS score and its subjective and objective support dimensions， as well as their combined use， all demonstrated good discriminative ability in distinguishing between the PTSD and non-PTSD groups. ConclusionThe results suggest that individuals who are female， with higher HAMA and HAMD scores after a motor vehicle accident， and those with lower social support and negative coping strategies， should be given particular attention. Early interventions for these individuals may reduce the incidence of PTSD.

Objective To evaluate the effectiveness of schistosomiasis surveillance in Jiangsu Province during the stage moving from transmission control to transmission interruption, and to analyze the current risk and challenges, so as to provide the evidence for achieving the target of schistosomiasis elimination. Methods Schistosomiasis surveillance data were collected from Jiangsu Province from 2012 to 2024, and the endemic areas, Schistosoma japonicum infections in humans and livestock, Oncomelania hupensis snail distribution and implementation of integrated interventions were descriptively analyzed. In addition, the trends in areas with snails, seroprevalence of human S. japonicum infections and numbers of advanced schistosomiasis cases were assessed using a Joinpoint regression model. Results The endemic areas of schistosomiasis continued to shrink in Jiangsu Province from 2012 to 2024, with the number of schistosomiasis-eliminated counties (cities, districts) increasing from 53 (75.71%) to 63 (96.92%), and interruption of schistosomiasis transmission was achieved across the province. A total of 4 600 300 person-times were tested for serum antibodies against S. japonicum, with 28 719 person-times positive detected; and 616 500 person-times were tested S. japonicum infections among local residents in Jiangsu Province from 2012 to 2024, with only 3 egg-positives detected, and no egg-positives found since 2017. A total of 187 600 herd-times were tested for schistosomiasis in livestock, and no S. japonicum infections were found. O. hupensis snail survey was performed covering 1 018 408.97 hm², and a total of 35 556.35 hm² was found with snail-infested habitats, including 174.40 hm2 of emerging snail-infested habitats. A total of 1 102 800 O. hupensis snails were identified for S. japonicum infections, and no infections were found. The areas of snail-infested habitats appeared a tendency towards a rise in Jiangsu Province from 2019 to 2023 (APC = 23.67%, P < 0.05), and the actual areas of snail-infested habitats appeared a tendency towards a decline from 2012 to 2015 (APC = −22.77%, P < 0.05), and towards a rise from 2015 to 2023 (APC = 9.76%, P < 0.01). The seroprevalence of anti-S. japonicum antibodies appeared a tendency towards a decline among residents in Jiangsu Province from 2017 to 2023 (APC = −14.92%, P < 0.01). In addition, the number of newly diagnosed advanced schistosomiasis cases appeared a tendency towards a decline from 2012 to 2024 (APC = −12.02%, P < 0.01), and the numbers of advanced schistosomiasis patients requiring treatment showed a tendency towards a decline from 2012 to 2021 (APC = −10.56%, P < 0.01) and from 2021 to 2023 (APC = −20.06%, P < 0.01). Conclusions Great progresses had been achieved in schistosomiasis control in Jiangsu Province following transmission control, and transmission interruption had been achieved; however, there are still snail-infested habitats. High-intensity surveillance and integrated control are required to be maintained to advance the achievement of the target of schistosomiasis elimination in Jiangsu Province.

Background Meteorological factors are among the key extrinsic triggers for the onset and exacerbation of cardiovascular and cerebrovascular diseases (CVD). Against the backdrop of sustained global warming, elucidating the impact of ambient temperature and atmospheric pressure on CVD, especially on pre-hospital CVD emergent events, has become imperative for evidence-based prevention and emergency preparedness. Objective To quantify the temporal trends of daily mean temperature and atmospheric pressure and their associations with pre-hospital CVD emergent events in Yantai, and to explore effect modification by demographic subgroups and geographic areas, thereby providing an empirical basis for the rational allocation of emergency medical resources. Methods Pre-hospital CVD emergency data from January 1, 2016 to December 31, 2022 were selected from the Yantai 120 Emergency Medical Command System. Synchronous meteorological factors and environmental pollutant data were obtained from the websites of the National Oceanic and Atmospheric Administration and the National Centers for Environmental Information of the United States. Time-series analysis combined with distributed lag non-linear model was used to analyze the association between daily temperature, atmospheric pressure, and pre-hospital CVD emergencies. Average annual percentage changes (AAPC) were calculated using Joinpoint (version 5.2.0.0) to reflect temporal trends. Spearman correlation analysis was employed to screen variables with low collinearity for inclusion in the multi-pollutant adjusted models. Results From 2016 to 2022, a total of 268275 pre-hospital CVD emergency cases were recorded in Yantai City, showing an increasing annual trend (AAPC=7.16%, P=0.002). Emergency volume, daily temperature, atmospheric pressure, and PM₁₀ displayed marked seasonal patterns, with summer bottoms and winter peaks for emergencies, daily temperature higher in summer and lower in winter, and the inverse for atmospheric pressure. Correlation analysis showed that both daily temperature and atmospheric pressure were positively correlated with O₃, with correlation coefficients of 0.745 and 0.723, respectively; negatively correlated with PM_2.5, with correlation coefficients of −0.246 and −0.235, respectively; and negatively correlated with PM₁₀, with correlation coefficients of −0.241 and −0.229, respectively (P < 0.05). A U-shaped correlation was revealed between daily temperature and pre-hospital CVD emergency risk. Using 24 °C as the reference, risk was elevated at both low and high daily temperatures, with the strongest cold effect at –11 ℃ (RR=1.53, 95%CI: 1.35, 1.73) and strongest heat effect at 30 ℃ (RR=1.06, 95%CI: 1.01, 1.11); cold effects were prolonged (lag0–14d) compared to heat effects. Atmospheric pressure exhibited a positive association with emergencies; using a median of 1012.2 hPa as the reference, the maximum risk occurred at 1037 hPa (RR=1.25, 95%CI: 1.03, 1.52), with a 0–5 d lag. In the subgroup analysis, the impact of low temperature was more pronounced among females (RR_{−9 ℃}=1.66), people over 60 years (RR_{−9 ℃}=1.91), and rural residents (RR_{−9 ℃}=1.76). High atmospheric pressure significantly increased the risk of pre-hospital CVD emergency among females (RR_{1035 hPa}=1.54), people under 18 years (RR_{1035 hPa}=3.62), and rural residents (RR_{1035 hPa}=1.46). The sensitivity analysis results indicated that, after excluding the impact of pandemic and the degrees of freedom for time variations, the effects of daily average temperature and atmospheric pressure on pre-hospital CVD emergency volume were consistent with the primary analysis findings. Conclusion Both non-optimal temperature and atmospheric pressure significantly increase the pre-hospital CVD emergency volume; therefore, emergency medical services should develop targeted response protocols, optimize resource allocation, and intensify public health education to mitigate the pre-hospital burden attributable to temperature and pressure fluctuations.

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.

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.