BACKGROUND:In China,the patient population with osteonecrosis is large,and there is an urgent need to find new preventive targets to develop more effective treatment strategies.Metabolomics studies have shown that there is an association between human metabolites and osteonecrosis,but the causal relationship between blood metabolites and osteonecrosis has not yet been clarified.OBJECTIVE:To investigate the causal relationship between blood metabolites and osteonecrosis through two-sample Mendelian randomization analysis.METHODS:The public data of 486 blood metabolites(exposure factors)and osteonecrosis(outcome factors)were collected.Data of 486 blood metabolites were derived from a genome-wide association estimate for blood metabolites published in Nature Genetics in 2014,which covered 7 824 European adults.The single nucleotide polymorphism data for osteonecrosis were obtained from the FinnGen public database R11 dataset,containing information on a total of 431 614 samples and 21 306 430 single nucleotide polymorphism loci,with 1 788 cases of osteonecrosis and 429 826 controls,with all participants being of European descent.Mendelian randomization analysis(inverse variance weighting method,MR-Egger method,and weighted median method)was performed by Rstudio software,and then the heterogeneity test,horizontal pleiotropy test and Steiger directionality test were performed to ensure the robustness and reliability of the results.RESULTS AND CONCLUSION:(1)Sixteen blood metabolites were identified as having a significant causal relationship with osteonecrosis(Pinverse variance weighting＜Pfalse discovery rate＜0.05).(2)Eight blood metabolites increased the risk of osteonecrosis(including four known metabolites and four unknown metabolites),specifically pantothenate,beta-hydroxyisovalerate,hippurate,salicyluric glucuronide,X-08766,X-11452,X-12776 and X-14662.(3)Eight blood metabolites could reduce the risk of osteonecrosis(six known metabolites and two unknown metabolites),including cortisol,1-palmitoylglycerol(1-monopalmitin),pyroglutamyl glycine,2-stearoylglycerophosphocholine,p-cresol sulfate,ergothioneine,X-06307,X-12092.(4)The above results suggest that there is a causal relationship between 16 blood metabolites and osteonecrosis,which is expected to be a potential target for intervention in the occurrence and treatment of osteonecrosis in the future.(5)Despite the lack of relevant data from large-scale Asian populations at present,this study provides important reference value for the field of osteonecrosis in China based on European population data.In the future,domestic medical workers may be able to achieve precise intervention for osteonecrosis by regulating metabolite levels.In addition,based on the results of this study,relevant researchers can further explore the mechanism of action of metabolites in the treatment of osteonecrosis with traditional Chinese medicine,which not only helps to deepen the understanding of traditional Chinese medical therapies but also promotes the progress of integrated traditional Chinese and Western medicine research,driving the development of personalized treatment plans that are more suitable for the characteristics of the Chinese population.

Objective: To propose an improved method for constructing the internal quality control (IQC) framework for ELISA assays and validate its efficacy by statistically analyzing IQC data from nine blood center laboratories. Methods: 1) IQC data was collected from nine blood centers and analyzed using a domestic HBsAg ELISA detection kit as an example. 2) Differences between IQC values across batches within Blood Center 1 were assessed. 3) Statistical analyses were performed on batch usage, number of batches used, days of use, number of QC points, batch-specific means, and coefficients of variation (CV) across all nine centers. 4) Using the improved construction method for IQC framework, provisional and permanent frames were established for batches within Blood Center 1 and Blood Center 9, followed by outlier determination. Results: 1) Statistically significant differences were observed in IQC data between batches within Blood Center 1 (P<0.01). It is recommended that both the control material/reagents and the control chart framework be replaced simultaneously. 2) There were substantial differences among 9 blood centers regarding the control material/reagent lot numbers used, the number of QC runs per batch, and the QC values for identical lots. Therefore, individual laboratories should establish their own IQC chart frameworks. 3) The improved IQC framework construction method for ELISA assays is as follows: provisional frames are established via frame-shifting, using the pre-experimental mean and cumulative coefficient of variation (CV) from the preceding batch. For batches used >20 days with >20 QC points, permanent frames are constructed by aggregating in-control data accumulated over ≥20 days with ≥20 points to calculate cumulative mean and standard deviation. The provisional and permanent frames constructed by this method identified all 26 extreme outliers across Blood Centers 1 and 9 as out-of-control. Among the 218 general outliers, 10 were classified as normal by the provisional frames, while the remainder were designated as warnings or out-of-control. This method effectively monitors assay stability. Conclusion: Based on the statistical analysis of IQC practices across blood centers of varying scales, combined with the inherent characteristics of ELISA assays and the batch-to-batch instability of reagents/QC materials, it is recommended to reconstruct QC charts upon lot changes. The proposed method—utilizing frame-shifting for provisional frames and establishing permanent frames based on cumulative data—is applicable to blood center laboratories of differing sizes and effectively monitors the stability of the ELISA assay process.

Objective To establish an oxidative stress mouse model of atopic dermatitis (AD) by applying 2,4-dinitrofluorobenzene (DNFB) to the back and post-auricular skin of KM mice, and to evaluate the regulatory role of the RAGE-NLRP3 axis (receptor for advanced glycation end products-NOD-like receptor family, pyrin domain containing 3 axis) in AD-related oxidative stress, thereby providing a potential therapeutic target for AD treatment. Methods Twenty SPF-grade female KM mice were randomly divided into a control group (Control group) and an experimental group (DNFB group), with 10 mice in each group. Mice in the Control group were treated with an acetone-olive oil vehicle (acetone: olive oil = 3:1) on their back and post-auricular skin. Mice in the DNFB group were treated with 0.5% DNFB (prepared by adding 0.5 g DNFB per 100 mL of acetone-olive oil vehicle) on the same areas, once daily for 14 consecutive days. The severity of skin lesions was scored on days 2, 4, 6, 9, 12, and 14 of treatment. On day 14, scratching behavior and ear thickness were evaluated. Ear swelling was evaluated on the final day by measuring bilateral ear thickness three times with a vernier caliper; the three measurements were averaged. HE staining was used to observe morphological and structural changes of cells in the back skin tissues. The mRNA and protein expression levels of RAGE (receptor for advanced glycation end products) in skin tissues were detected by quantitative real-time PCR, Western blot, and immunohistochemical staining. The mRNA expression levels of oxidative stress-related molecules, including NLRP3 (NOD-like receptor family, pyrin domain containing 3), caspase-1 (cysteine-dependent aspartate-specific protease 1), and IL-1β (Interleukin-1β), were detected by quantitative real-time PCR. Results On day 14, the back skin lesion scores of the Control group and DNFB group were (0.20±0.42) and (9.93±1.30) (P<0.000 1), respectively. Scratching behavior scores were (5.00±2.05) and (49.26±8.49) episodes, respectively (P<0.000 1), and ear thicknesses were (213.00±11.87) μm and (765.93±140.47) μm (P<0.000 1), respectively. The DNFB group exhibited marked skin dryness, desquamation, and thickening. HE staining results showed that skin inflammation was obvious in the DNFB group, consistent with the pathological features of AD. Quantitative real-time PCR and Western blot results showed that compared with the Control group, the mRNA expression level of RAGE in skin tissues of the DNFB group was significantly increased (P<0.05), and the protein expression level of RAGE was also significantly increased (P<0.01). Immunohistochemical staining results showed that compared with the Control group, skin tissue sections of the DNFB group exhibited thickened stratum corneum and fibrotic proliferation of fibroblasts in the interstitium under microscopic observation, with a significant increase in RAGE protein expression in the skin tissues (P<0.01). Quantitative real-time PCR results showed that the mRNA expression levels of NLRP3, caspase-1, and IL-1β in skin tissues of the DNFB group were all significantly increased (P<0.01). Conclusion The AD mouse oxidative stress model has been successfully established by topical DNFB application. RAGE may promote the development of AD by regulating the NLRP3 inflammasome and IL-1β release, forming an oxidative-inflammatory cascade, suggesting that it could be a potential therapeutic target for AD.

Glaucoma， characterized by pathological elevation of intraocular pressure， progressive optic nerve damage， and apoptosis of retinal ganglion cells， is one of the leading causes of irreversible blindness. It is strongly associated with significant vision loss and a decline in quality of life. Although pharmacological therapy remains the primary approach to managing glaucoma， clinical outcomes are often suboptimal， highlighting the urgent need for safe and effective alternative treatments. In traditional Chinese medicine （TCM）， glaucoma is categorized as part of the "five wind internal obstruction" syndrome， and TCM has amassed substantial experience in the prevention and treatment of this condition. Therefore， this article provided a comprehensive review of recent findings on the relationship between glaucoma and relevant signaling pathways， as well as the regulatory effects of TCM on these pathways in the treatment of glaucoma. TCM can exert therapeutic effects by modulating key signaling pathways， including the Toll-like receptor 4/myeloid differentiation factor 88/nuclear factor-κB （TLR4/MyD88/NF-κB） signaling pathway， Janus kinase （JAK）/signal transducer and activator of transcription （JAK/STAT） pathway， phosphoinositide 3-kinase/serine-threonine protein kinase/mammalian target of rapamycin （PI3K/Akt/mTOR） pathway， nucleotide-binding oligomerization domain-like receptor protein 3 （NLRP3） inflammasome pathway， and nuclear factor erythroid 2-related factor 2/heme oxygenase-1 （Nrf2/HO-1） pathway. These pathways are involved in reducing inflammation， inhibiting apoptosis and ferroptosis， and ameliorating oxidative stress. By synthesizing current research， this article offers theoretical insights and practical references for advancing the understanding of the pathological mechanisms underlying glaucoma， innovating strategies for optic nerve protection， and promoting integrative TCM and Western medical approaches in glaucoma management.

Objective:To investigate the association between small vulnerable newborn (SVN) phenotypes and the risk of neurodevelopmental delay at the age of 1 year.Methods:A prospective cohort study was conducted. A total of 25 860 singleton infants from "The Born in Guangzhou Cohort Study" who completed the Gesell developmental scale assessment at 1 year of age between January 2013 and June 2025 were included. Maternal sociodemographic characteristics, and other information were collected using a self-administered questionnaire, and maternal pregnancy-related information and neonatal birth data were extracted from medical records. Global developmental delay (GDD) was defined as a developmental quotient below 86 in ≥3 domains of the Gesell developmental scale, which assesses the adaptive, gross motor, fine motor, language, and personal-social domains. The random forest algorithm was employed for missing data imputation. Based on prematurity, small for gestational age (SGA), and low birth weight (LBW), newborns were categorized into 6 phenotypes: preterm-SGA-LBW, preterm-appropriate for gestational age (AGA)-LBW, preterm-AGA-nonLBW, term-SGA-LBW, term-LBW-only or term-SGA-only, and term-AGA-nonLBW phenotype. Among these, the first 5 were classified as SVN phenotypes, and the last one served as the reference group. Inter-group comparisons were performed using analysis of variance (ANOVA), χ2 tests, or Kruskal-Wallis test, as appropriate.?? Multivariable robust Poisson regression models were applied to analyze the association of different SVN phenotypes with the risks of GDD and developmental delays in specific domains, with stratified analyses by sex. Results:Among the 25 860 infants, 13 719 (53.1%) were male and 12 141 (46.9%) were female. The gestational age at birth was 39.4 (38.6, 40.0) weeks. The overall detection rate of GDD at 1 year of age was 3.7% (962/25 860). The rates of delay across developmental domains, in descending order, language in 8 134 cases (31.5%), gross motor in 4 488 cases (17.4%), personal-social in 1 271 cases (4.9%), adaptive in 1 262 cases (4.9%), and fine motor in 621 cases (2.4%). Compared with the reference group, preterm-AGA-LBW, preterm-SGA-LBW, preterm-AGA-noneLBW, and term-SGA-LBW phenotypes were all associated with an increased risk of GDD, with the adjusted RR (95% CI) of 6.07(5.01-7.35), 4.81(3.11-7.46), 2.10(1.54-2.88) and 1.89(1.29-2.76) respectively.The preterm-AGA-noneLBW phenotype was all associated with an increased risk of delay in gross motor, language and personal-social functional domains (all P<0.05). The term-SGA-LBW phenotype was associated with an increased risk of delay in gross motor, fine motor and personal-social functional domains (all P<0.01). Whereas the term-LBW-only or term-SGA-only phenotype showed no statistically association with developmental delay in any functional domain (all P≥0.05). Conclusion:The combined classification based on gestational age and birth weight helps identify infants at high risk for neurodevelopmental delay at 1 year of age, suggesting that it may offer a reference for the rational allocation of clinical resources.

BACKGROUND:In China,the patient population with osteonecrosis is large,and there is an urgent need to find new preventive targets to develop more effective treatment strategies.Metabolomics studies have shown that there is an association between human metabolites and osteonecrosis,but the causal relationship between blood metabolites and osteonecrosis has not yet been clarified.OBJECTIVE:To investigate the causal relationship between blood metabolites and osteonecrosis through two-sample Mendelian randomization analysis.METHODS:The public data of 486 blood metabolites(exposure factors)and osteonecrosis(outcome factors)were collected.Data of 486 blood metabolites were derived from a genome-wide association estimate for blood metabolites published in Nature Genetics in 2014,which covered 7 824 European adults.The single nucleotide polymorphism data for osteonecrosis were obtained from the FinnGen public database R11 dataset,containing information on a total of 431 614 samples and 21 306 430 single nucleotide polymorphism loci,with 1 788 cases of osteonecrosis and 429 826 controls,with all participants being of European descent.Mendelian randomization analysis(inverse variance weighting method,MR-Egger method,and weighted median method)was performed by Rstudio software,and then the heterogeneity test,horizontal pleiotropy test and Steiger directionality test were performed to ensure the robustness and reliability of the results.RESULTS AND CONCLUSION:(1)Sixteen blood metabolites were identified as having a significant causal relationship with osteonecrosis(Pinverse variance weighting＜Pfalse discovery rate＜0.05).(2)Eight blood metabolites increased the risk of osteonecrosis(including four known metabolites and four unknown metabolites),specifically pantothenate,beta-hydroxyisovalerate,hippurate,salicyluric glucuronide,X-08766,X-11452,X-12776 and X-14662.(3)Eight blood metabolites could reduce the risk of osteonecrosis(six known metabolites and two unknown metabolites),including cortisol,1-palmitoylglycerol(1-monopalmitin),pyroglutamyl glycine,2-stearoylglycerophosphocholine,p-cresol sulfate,ergothioneine,X-06307,X-12092.(4)The above results suggest that there is a causal relationship between 16 blood metabolites and osteonecrosis,which is expected to be a potential target for intervention in the occurrence and treatment of osteonecrosis in the future.(5)Despite the lack of relevant data from large-scale Asian populations at present,this study provides important reference value for the field of osteonecrosis in China based on European population data.In the future,domestic medical workers may be able to achieve precise intervention for osteonecrosis by regulating metabolite levels.In addition,based on the results of this study,relevant researchers can further explore the mechanism of action of metabolites in the treatment of osteonecrosis with traditional Chinese medicine,which not only helps to deepen the understanding of traditional Chinese medical therapies but also promotes the progress of integrated traditional Chinese and Western medicine research,driving the development of personalized treatment plans that are more suitable for the characteristics of the Chinese population.

Objective:To investigate the association between small vulnerable newborn (SVN) phenotypes and the risk of neurodevelopmental delay at the age of 1 year.Methods:A prospective cohort study was conducted. A total of 25 860 singleton infants from "The Born in Guangzhou Cohort Study" who completed the Gesell developmental scale assessment at 1 year of age between January 2013 and June 2025 were included. Maternal sociodemographic characteristics, and other information were collected using a self-administered questionnaire, and maternal pregnancy-related information and neonatal birth data were extracted from medical records. Global developmental delay (GDD) was defined as a developmental quotient below 86 in ≥3 domains of the Gesell developmental scale, which assesses the adaptive, gross motor, fine motor, language, and personal-social domains. The random forest algorithm was employed for missing data imputation. Based on prematurity, small for gestational age (SGA), and low birth weight (LBW), newborns were categorized into 6 phenotypes: preterm-SGA-LBW, preterm-appropriate for gestational age (AGA)-LBW, preterm-AGA-nonLBW, term-SGA-LBW, term-LBW-only or term-SGA-only, and term-AGA-nonLBW phenotype. Among these, the first 5 were classified as SVN phenotypes, and the last one served as the reference group. Inter-group comparisons were performed using analysis of variance (ANOVA), χ2 tests, or Kruskal-Wallis test, as appropriate.?? Multivariable robust Poisson regression models were applied to analyze the association of different SVN phenotypes with the risks of GDD and developmental delays in specific domains, with stratified analyses by sex. Results:Among the 25 860 infants, 13 719 (53.1%) were male and 12 141 (46.9%) were female. The gestational age at birth was 39.4 (38.6, 40.0) weeks. The overall detection rate of GDD at 1 year of age was 3.7% (962/25 860). The rates of delay across developmental domains, in descending order, language in 8 134 cases (31.5%), gross motor in 4 488 cases (17.4%), personal-social in 1 271 cases (4.9%), adaptive in 1 262 cases (4.9%), and fine motor in 621 cases (2.4%). Compared with the reference group, preterm-AGA-LBW, preterm-SGA-LBW, preterm-AGA-noneLBW, and term-SGA-LBW phenotypes were all associated with an increased risk of GDD, with the adjusted RR (95% CI) of 6.07(5.01-7.35), 4.81(3.11-7.46), 2.10(1.54-2.88) and 1.89(1.29-2.76) respectively.The preterm-AGA-noneLBW phenotype was all associated with an increased risk of delay in gross motor, language and personal-social functional domains (all P<0.05). The term-SGA-LBW phenotype was associated with an increased risk of delay in gross motor, fine motor and personal-social functional domains (all P<0.01). Whereas the term-LBW-only or term-SGA-only phenotype showed no statistically association with developmental delay in any functional domain (all P≥0.05). Conclusion:The combined classification based on gestational age and birth weight helps identify infants at high risk for neurodevelopmental delay at 1 year of age, suggesting that it may offer a reference for the rational allocation of clinical resources.

ObjectiveTo investigate the contamination status of ochratoxin A (OTA) in commercially available food products in Shanghai, and to assess OTA exposure levels and the associated non-carcinogenic and carcinogenic risks among pregnant women by integrating dietary consumption data of this population. MethodsThe levels of OTA contamination in 1 520 food samples collected in Shanghai from 2022 to 2023 were determined using liquid chromatography-tandem mass spectrometry. An exposure assessment model was developed based on the dietary consumption levels of pregnant women from the 2016‒2017 Shanghai Pregnant Women Dietary Monitoring Survey to calculate the estimated daily intake (EDI) of OTA, the margin of exposure for non-carcinogenic toxicity (MOE₁), and the margin of exposure for carcinogenic toxicity (MOE₂). An MOE₁ greater than 200 and an MOE₂ greater than 10 000 indicate that the non-carcinogenic toxicity and carcinogenic toxicity resulting from exposure are negligible, respectively. For samples with OTA contamination levels below the limit of detection (LOD), which accounted for more than 80% of the samples, the OTA levels were assigned values of 0 and LOD, respectively, for subsequent calculations. ResultsThe detection rates of OTA in cereals, nuts, dried fruits, and alcohol samples collected in 2022 were 2.03%, 0, 0, and 0, respectively. The OTA detection rates in cereals, nuts, dried fruits, beans, and alcohol samples collected in 2023 were 2.50%, 0.39%, 2.47%, 1.67%, and 13.33%, respectively. For pregnant women in Shanghai in 2022, simulation results indicated that when assigning a value of 0 and the LOD, theP₅₀ values of EDI for dietary OTA exposure were 0.05 and 0.72 ng·(kg·d)^-1, respectively, and the P₉₅ values of EDI for dietary OTA exposure were 0.25 and 2.40 ng·(kg·d)^-1, respectively. For pregnant women in Shanghai in 2023, the P₅₀ values of EDI for dietary OTA exposure were 0.04 and 1.00 ng·(kg·d)^-1, respectively, and the P₉₅ values of EDI for dietary OTA exposure were 0.23 and 2.67 ng·(kg·d)^-1, respectively, both substantially below the tolerable daily intake (TDI) for OTA [17 ng·(kg·d)^-1]. The EDI for dietary OTA exposure in 100.0% of Shanghai pregnant women was lower than the TDI, indicating an overall low level of dietary OTA exposure among this population. For 100.0% of pregnant women, the MOE₁ for dietary OTA exposure exceeded 200. When assigned a value of 0, the MOE₂ for 100.0% of pregnant women in both 2022 and 2023 exceeded10 000. When assigned the LOD value, 72.3% and 81.8% of pregnant women in 2022 and 2023, respectively, had an MOE₂ exceeding 10 000. ConclusionFrom 2022 to 2023, samples of cereals, nuts, dried fruits, beans, and alcohol sold in Shanghai exhibited varying degrees of OTA contamination. The overall EDI of OTA exposure among pregnant women in Shanghai remained at a low level. The non-carcinogenic and carcinogenic risks associated with OTA exposure were generally low and at controllable levels.

ObjectiveThis paper aims to observe the syndrome improvement and negative antigen conversion rate of Qingxuan Daozhi formula in the treatment of influenza in children （heat accumulation in the lung and stomach syndrome）. MethodsThrough a multi-center randomized controlled methodology design，confirmed influenza cases were collected from October 2022 to April 2023 in the pediatrics department of eight hospitals，such as Dongfang Hospital of Beijing University of Chinese Medicine. A total of 180 children with influenza and heat accumulation in the lung and stomach syndrome conforming to the standard were recruited through the clinic. The sick children meeting the inclusion criteria were randomly divided into groups by a block-randomized method. The children in the experimental group were treated with Qingxuan Daozhi formula for five days，and those in the control group were treated with Oseltamivir Phosphate Granules for five days. The primary efficacy indicator was the negative conversion rate of influenza antigen detection. Secondary efficacy indicators were the Canadian acute respiratory illness and flu scale （CARIFS） and the incidence of complications，severe cases， and critical cases. Follow-up observation was conducted on the day of enrollment，48 hours after medication，72 hours after medication， and （6+1） d after medication. ResultsOne hundred and eighty participants were randomly assigned to the experimental group （90 cases） or the control group （90 cases）. All participants were followed up during the study. Comparison of influenza antigen detection results in the primary efficacy indicators showed that the average time of negative influenza antigen conversion in the experimental group was （5.29±1.25） d，and that in the control group was （5.40±1.68） d，without a statistically significant difference. After five days of intervention，52 cases in the experimental group and 51 cases in the control group converted to negative，without a statistically significant difference. CARIFS score results in the secondary efficacy indicators showed that during 72 hours after intervention，there were statistically significant differences between the experimental group and the control group in three dimensions， including headache，muscle soreness， and the need for extra care （P<0.05）. On the （6+1） days after the intervention，the differences in both the experimental group and the control group were statistically significant in 10 dimensions， including sore throat，bad sleep，uncomfortable feeling，poor spirit and fatigue，crying more than usual，the need for extra care，symptom，function，influence on parents，and total score （P<0.05）. The comparison results within the group in the dimensional scores of symptom， function， and influence on parents，as well as the CARIFS total score showed that with the delay of follow-up time，scores of both groups decreased significantly，with a statistically significant difference （P<0.01）. Inter-group comparison results showed that the mean score of the experimental group was higher than that of the control group at the time of enrollment. With the progress of intervention，the score of the experimental group was significantly decreased compared with that of the control group. At the end of follow-up，the mean score of the experimental group was lower than that of the control group，with no statistically significant difference. In terms of the incidence of complications，severe cases， and critical cases， there were no complications，severe cases， and critical cases in the two groups，without a statistically significant difference. ConclusionThe symptom improvement effect and negative antigen conversion rate of Qingxuan Daozhi formula in the treatment of influenza in children （heat accumulation in the lung and stomach syndrome） are not inferior to Oseltamivir Phosphate granules， and children's acceptance is better. It can be more widely used in clinical treatment of influenza in children （heat accumulation in the lung and stomach syndrome）.

Background/Aims: Lenvatinib resistance remains a critical barrier in advanced hepatocellular carcinoma (HCC) therapy. However, the underlying mechanisms and strategies for reversing resistance remain incompletely understood. Methods: Integrated transcriptomics of lenvatinib-resistant patient tumors and an acquired-resistance murine model identified a novel macrophage subpopulation. Functional validation employed CRISPR-SAM screening, conditioned medium (CM) assays, subcutaneous/orthotopic xenografts, patient-derived organoids (PDOs), and patient-derived xenografts (PDXs). Mechanistic studies included ChIP-qPCR, co-immunoprecipitation, and pharmacologic targeting. Clinical relevance was assessed in a retrospective cohort. Results: Resistant HCC exhibited significant enrichment of a COLEC12high TAM subset , which correlated with poor survival and treatment response. These TAMs secreted neuregulin-1 (NRG1) , activating HER2/HER3-AKT signaling in tumor cells to drive cancer stemness and lenvatinib resistance. Mechanistically, in TAMs COLEC12 sequestered STAT1 in the cytoplasm, preventing its phosphorylation, and thereby derepressing STAT3-mediated NRG1 transcription. Depletion of NRG1 reversed the stemness phenotypes and resensitized tumors to lenvatinib both in vitro and in vivo. Clinically, high NRG1 expression predicted an inferior lenvatinib response and shorter survival. Crucially, the bispecific anti-HER2/HER3 antibody zenocutuzumab restored lenvatinib efficacy in PDOs, PDXs, and murine models. Conclusions Our work establishes the COLEC12high TAM/NRG1 axis as a master regulator of therapeutic resistance and identifies NRG1 as a predictive biomarker, providing a clinically actionable strategy to overcome lenvatinib resistance in HCC.