ObjectiveTo evaluate the impact of yang-reinforcing and blood-activating therapy on the long-term prognosis for patients with dilated cardiomyopathy （DCM） of yang deficiency and blood stasis syndrome. MethodsA retrospective cohort study was conducted involving 371 DCM patients with yang deficiency and blood stasis syndrome. The yang-reinforcing and blood-activating therapy was defined as the exposure factor. Patients were categorized into exposure group （186 cases） and non-exposure group （185 cases） according to whether they received yang-reinforcing and blood-activating therapy combined with conventional western medicine for 6 months or longer. The follow-up period was set at 48 months， and the Kaplan-Meier survival analysis was used to assess the cumulative incidence of major adverse cardiovascular events （MACE） in both groups. Cox regression analysis was used to explore the impact of yang-reinforcing and blood-activating therapy on the risk of MACE， and subgroup analysis was performed. Changes in traditional Chinese medicine （TCM） syndrome score， left ventricular ejection fraction （LVEF）， left ventricular fractional shortening （LVFS）， left ventricular end-diastolic diameter （LVEDD）， and Minnesota Living with Heart Failure Questionnaire （MLHFQ） score were compared between groups at the time of first combined use of yang-reinforcing and blood-activating therapy （before treatment） and 1 year after receiving the therapy （after treatment）. ResultsMACE occurred in 31 cases （16.67%） in the exposure group and 47 cases （25.41%） in the non-exposure group. The cumulative incidence of MACE in the exposure group was significantly lower than that in the non-exposure group ［HR=0.559， 95%CI（0.361，0.895）， P=0.014］. Cox regression analysis showed that yang-reinforcing and blood-activating therapy was an independent factor for reducing the risk of MACE in DCM patients ［HR=0.623， 95%CI（0.396，0.980）， P=0.041］， and consistent results were observed in different subgroups. Compared with pre-treatment， the exposure group showed decreased TCM syndrome score and MLHFQ score， reduced LVEDD， and increased LVEF and LVFS after treatment （P＜0.05）； in the non-exposure group， TCM syndrome score decreased， LVEF and LVFS increased， and LVEDD reduced after treatment （P＜0.05）. After treatment， the exposure group had higher LVEF and LVFS， smaller LVEDD， and lower TCM syndrome score and MLHFQ score compared with the non-exposure group （P＜0.05）. ConclusionCombining yang-reinforcing and blood-activating therapy with conventional western medicine can reduce the risk of MACE in DCM patients with yang deficiency and blood stasis syndrome， meanwhile improving their clinical symptoms， cardiac function， and quality of life.

Objective To elucidate the molecular mechanism of sirtuin-3 (SIRT3) in regulating mitochondrial biogenesis in human renal tubular epithelial cells. Methods Cells were stimulated with different concentrations of H₂O₂ and divided into four groups: control (NC), 50 μmol/L H₂O₂, 110 μmol/L H₂O₂ and 150 μmol/L H₂O₂. SIRT3 protein expression was then measured. SIRT3 was knocked down with siRNA, and cells were further assigned to five groups: control (NC), negative-control siRNA (NCsi), SIRT3-siRNA (siSIRT3), NCsi+H₂O₂, and siSIRT3+H₂O₂. After 24 h, cellular adenosine triphosphate (ATP) and mitochondrial superoxide anion (O₂•⁻) levels were determined, together with mitochondrial expression of SIRT3, peroxisome proliferator-activated receptor γ coactivator-1α (PGC-1α), nuclear respiratory factor 1 (NRF1), mitochondrial transcription factor A (TFAM), superoxide dismutase 2 (SOD2), acetylated-SOD2 and adenosine monophosphate activated protein kinase α1 (AMPKα1). Results The 110 and 150 μmol/L H₂O₂ decreased SIRT3 protein (both P<0.05). ATP and mitochondrial O₂•⁻ did not differ between NC and NCsi groups (both P>0.05). Compared to the NCsi group, the siSIRT3 group exhibited elevated O₂•⁻ level, decreased SIRT3 protein and increased expression levels of SOD2 and acetylated SOD2 protein (all P<0.05). Compared to the NCsi group, the NCsi+H₂O₂ group exhibited decreased cellular ATP levels, elevated mitochondrial O₂•⁻ levels, and reduced protein expression levels of SIRT3, SOD2, TFAM, AMPKα1, PGC-1α and NRF1 (all P<0.05). Compared with the siSIRT3 group, the siSIRT3+H₂O₂ group showed a decrease in cellular ATP levels, an increase in mitochondrial O₂•⁻ levels, a decrease in SIRT3, SOD2, TFAM, AMPKα1, PGC-1α and NRF1 protein expression levels and a decrease in acetylated SOD2 protein expression levels (all P<0.05). Compared with the NCsi+H₂O₂ group, the siSIRT3+H₂O₂ group showed a decrease in cellular ATP levels, an increase in mitochondrial O₂•⁻ levels, a decrease in SIRT3, AMPKα1, PGC-1α and NRF1, TFAM protein expression levels, and an increase in SOD2 and acetylated SOD2 protein expression levels (all P<0.05). Conclusions SIRT3 promotes mitochondrial biogenesis in tubular epithelial cells via the AMPK/PGC-1α/NRF1/TFAM axis, representing a key mechanism through which SIRT3 ameliorates oxidative stress-induced mitochondrial dysfunction.

ObjectiveTo investigate the present situation of input, output, outcome and impact of all registered community-based rehabilitation stations in Inner Mongolia in China, and analyze how the input predict the output, outcome and impact. MethodsFrom March 1st to April 30th, 2025, a questionnaire survey was conducted on all registered community-based rehabilitation stations in Inner Mongolia, covering four dimensions: input, output, outcome and impact. A total of 1 365 questionnaires were distributed. The input included four items: laws and policies, human resources, equipment and facilities, and rehabilitation information management. The output included two items: technical paths and benefits/effectiveness. The outcome included three items: coverage rates, rehabilitation interventions and functional results. The impact included two items: health and sustainability. Each item contained several questions, all of which were described in a positive way. Each question was scored from one to five. A lower score indicated that the situation of the community-based rehabilitation station was more in line with the content described in the question. Regression analysis was performed using the total score of each item of input dimension as independent variables, and the total scores of the output, outcome and impact dimensions as dependent variables. ResultsA total of 1 262 valid questionnaires were collected. The mean values of input, output, outcome and impact of community-based rehabilitation stations were 1.827 to 1.904, with coefficient of variation of 45.892% to 49.239%. The regression analysis showed that, rehabilitation information management, human resources, and laws and policies significantly predicted the output dimension (R² = 0.910, P < 0.001). Meanwhile, all four items in the input dimension predicted both the outcome (R² = 0.850, P < 0.001) and impact dimensions (R² = 0.833, P < 0.001). ConclusionInput, output, outcome and impact of the community-based rehabilitation stations in Inner Mongolia were generally in line with the content of the questions, although some imbalances were observed. Additionally, the input of community-based rehabilitation stations could significantly predict their output, outcome and impact.

Oral squamous cell carcinoma (OSCC) is a common head and neck malignancy. Approximately 50% to 60% of patients with OSCC are diagnosed at a locally advanced stage (clinical staging III-IVa). Even with comprehensive and sequential treatment primarily based on surgery, the 5-year overall survival rate remains below 50%, and patients often suffer from postoperative functional impairments such as difficulties with speaking and swallowing. Programmed death receptor-1 (PD-1) inhibitors are increasingly used in the neoadjuvant treatment of locally advanced OSCC and have shown encouraging efficacy. However, clinical practice still faces key challenges, including the definition of indications, optimization of combination regimens, and standards for efficacy evaluation. Based on the latest research advances worldwide and the clinical experience of the expert group, this expert consensus systematically evaluates the application of PD-1 inhibitors in the neoadjuvant treatment of locally advanced OSCC, covering combination strategies, treatment cycles and surgical timing, efficacy assessment, use of biomarkers, management of special populations and immune related adverse events, principles for immunotherapy rechallenge, and function preservation strategies. After multiple rounds of panel discussion and through anonymous voting using the Delphi method, the following consensus statements have been formulated: 1) Neoadjuvant therapy with PD-1 inhibitors can be used preoperatively in patients with locally advanced OSCC. The preferred regimen is a PD-1 inhibitor combined with platinum based chemotherapy, administered for 2-3 cycles. 2) During the efficacy evaluation of neoadjuvant therapy, radiographic assessment should follow the dual criteria of Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1 and immune RECIST (iRECIST). After surgery, systematic pathological evaluation of both the primary lesion and regional lymph nodes is required. For combination chemotherapy regimens, PD-L1 expression and combined positive score need not be used as mandatory inclusion or exclusion criteria. 3) For special populations such as the elderly (≥ 70 years), individuals with stable HIV viral load, and carriers of chronic HBV/HCV, PD-1 inhibitors may be used cautiously under the guidance of a multidisciplinary team (MDT), with close monitoring for adverse events. 4) For patients with a poor response to neoadjuvant therapy, continuation of the original treatment regimen is not recommended; the subsequent treatment plan should be adjusted promptly after MDT assessment. Organ transplant recipients and patients with active autoimmune diseases are not recommended to receive neoadjuvant PD-1 inhibitor therapy due to the high risk of immune related activation. Rechallenge is generally not advised for patients who have experienced high risk immune related adverse events such as immune mediated myocarditis, neurotoxicity, or pneumonitis. 5) For patients with a good pathological response, individualized de escalation surgery and function preservation strategies can be explored. This consensus aims to promote the standardized, safe, and precise application of neoadjuvant PD-1 inhibitor strategies in the management of locally advanced OSCC patients.

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.

BACKGROUND:Studies have shown that rheumatoid arthritis and osteoporosis are positively correlated,but the causal relationship and related mechanisms have not yet been confirmed.With the cross-fertilization of computer science and life sciences,Mendelian randomization and bioinformatics analyses based on genome-wide association study(GWAS)and transcriptome sequencing data can assess the causal relationship between two diseases,explore the related mechanisms,and mine the therapeutic targets,which will be beneficial to the precision treatment of rheumatoid arthritis combined with osteoporosis.OBJECTIVE:To explore the causal relationship between rheumatoid arthritis and osteoporosis using two-sample Mendelian randomization and to mine potential co-morbid targets and potential targeted drugs through summary-data-based Mendelian randomization and bioinformatics analyses,aiming to provide theoretical basis for mechanism exploration and precision treatment in the field of rheumatoid arthritis combined with osteoporosis.METHODS:(1)Firstly,GWAS data of rheumatoid arthritis,osteoporosis,and cis-expression quantitative trait locus(cis-eQTL)in Asian and European populations were downloaded from the GWAS Catalog,IEU Open GWAS,FinnGen,and eQTLGen databases,and were used for two-sample Mendelian randomization analysis and summary-data-based Mendelian randomization analysis.(2)Transcriptome sequencing data of rheumatoid arthritis(GSE93272 and GSE15573)were downloaded from the GEO database for bioinformatics analysis.(3)Subsequently,forward and inverse Mendelian randomization analyses between rheumatoid arthritis and osteoporosis were performed,and inverse variance weighted was used as the main metric for the analyses,and the results were corroborated with MR Egger,simple mode,weighted median and weighted mode.(4)Then,the genes closely related to rheumatoid arthritis and osteoporosis were identified based on the summary-data-based Mendelian randomization analysis,and the co-disease targets of rheumatoid arthritis and osteoporosis were mined based on cross-analysis.Meanwhile,the biological functions of the co-morbid targets were verified based on bioinformatics analysis and cellular experiments.(5)In addition,a rheumatoid arthritis risk prediction nomogram was constructed based on DYRK2,and its prediction performance was verified by receiver operating characteristic curve,correction curve and decision curve.Finally,the target potential drugs were mined based on Enrichr database and molecular docking was performed.RESULTS AND CONCLUSION:(1)Forward Mendelian randomization analysis of rheumatoid arthritis and osteoporosis showed statistically significant results except for GCST90044540 and GCST90086118,and all other results indicated a significant causal relationship and positive correlation between rheumatoid arthritis and osteoporosis.(2)Inverse Mendelian randomization analysis suggested that no significant causal relationship was seen between osteoporosis and rheumatoid arthritis.(3)Summary-data-based Mendelian randomization analysis identified a total of 412 and 344 genes positively associated with rheumatoid arthritis and osteoporosis,and 421 and 347 genes negatively associated.Based on the cross-analysis,26 co-morbid genes were subsequently obtained.Among them,DYRK2 was a potential therapeutic target,and subsequent bioinformatics analysis and cellular experiments confirmed its important role in the progression of rheumatoid arthritis and osteoporosis.(4)Furthermore,the constructed nomogram has excellent predictive performance.Finally,four potential DYRK2-targeting drugs(undecanoic acid,metyrapone,JNJ-38877605,and ACA)were discovered and molecular docking also demonstrated reliable targeting ability.(5)In conclusion,based on GWAS data from Asian and European populations,we successfully demonstrated that rheumatoid arthritis and osteoporosis are causally related at the genetic level,DYRK2 is a potential therapeutic target,and four small molecules are potential target drugs.