OBJECTIVE To analyze the specific risks and long-term toxicities of four BCR-ABL1 tyrosine kinase inhibitors （TKIs）（imatinib， dasatinib， nilotinib， and bosutinib） in pediatric patients with hematological malignancies. METHODS Adverse drug event （ADE） reports submitted to the the United States FDA Adverse Event Reporting System （FAERS） from January 2012 to December 2024， with imatinib， dasatinib， nilotinib， and bosutinib as the primary suspect drugs， were collected. Data mining was performed using the reporting odds ratio method and proportional reporting ratio method. ADE terms were classified and summarized by system organ class （SOC） and preferred term （PT） according to the Medical Dictionary for Drug Regulatory Activities （MedDRA， version 26.0）. Meanwhile， the ADE reports were divided by age into the adult group （≥18 years） and the pediatric group （＜18 years） to compare the differences in ADE between the two groups. RESULTS A total of 1 512 pediatric ADE reports were included： 993 for imatinib， 391 for dasatinib， 112 for nilotinib， and 16 for bosutinib. Among the reported ADEs， the patients were mainly aged 12-＜18 years； the reports mainly originated from the United States， France， and Japan； and the primary indications were chronic myeloid leukemia and acute lymphoblastic leukemia. A total of 5 256 ADE signals were mined， among which 235 were positive signals， involving 1 103 PT across 27 SOC. The top five PT ranked by the number of positive signals were nausea， febrile neutropenia， abdominal pain， neutropenia， and anemia. The top two SOC were general disorders and administration site conditions， and gastrointestinal disorders. Compared with the adult group， the pediatric group had relatively higher proportions of events related to infections and infestations as well as blood and lymphatic system disorders. Pediatric long-term toxicity signals primarily included growth retardation， accompanied by signals related to endocrine system abnormalities and bone metabolism abnormalities. Specific signals included imatinib-associated septic shock， dasatinib-associated chylothorax， and nilotinib-associated electrocardiographic QT interval prolongation. CONCLUSIONS When pediatric patients use BCR-ABL1 TKIs， priority monitoring of infection risk and hematologic parameters is required， along with long-term follow-up of height， endocrine， and bone metabolism parameters. Targeted screening and management of drug-specific signals should be performed to ensure the long-term safety of pediatric medication.

“Bio-narrative integration” refers to the process in which a life subject with narrative consciousness actively reviews and reintegrates their life stories， or a subject lacking life and health narrative awareness， with the intervention of healthcare professionals， tells their own life stories and integrates them into a coherent and constantly evolving life narrative process. Starting from the keyword of bio-narrative integration， this paper proposed a classification model of narrative integration. From the perspective of life stages， it was divided into “phasic narrative integration” and “holistic narrative integration.” In terms of integrated narrative style， it was categorized as “positive narrative integration style” and “negative narrative integration style.” Regarding subjective initiative， it was classified as “active narrative integration regulation” and “passive narrative integration regulation.” Then it elaborated the significant value of narrative integration for every life subject， especially in pain relief， the improvement of life resilience， the healthy aging of the elderly， and the ultimate peace of the dying. It was advocated that healthcare practitioners should enhance their professional narrative competence， effectively guide patients to engage in bio-narrative integration regulation， and help them overcome narrative closure， thereby improving the quality of medical care.

ObjectiveThis paper aims to find the identified and validated clinical biomarker data building upon a clinical study of early-phase phase Ⅱ and investigate the correlation analysis of Huanglian Jiedu Wan on syndrome improvement and clinical biomarkers in the treatment of "excess heat-toxicity" based on a machine learning model. Additionally， the effective prediction of clinical biomarker values for the main symptoms of the "excess heat-toxicity" syndrome was assessed. MethodsA total of 229 patients meeting the inclusion criteria for "excess heat-toxicity" syndrome were randomly divided into the Huanglian Jiedu Wan group and the placebo group. Syndrome score transition matrices were constructed for the Huanglian Jiedu Wan group and the placebo group based on three main symptoms of "excess heat-toxicity" syndrome， such as oral ulcers， sore throat， and gum swelling and pain. Data from the patients with these three syndromes were also integrated for an overall analysis. The corresponding syndrome score transition matrices were further constructed to visualize symptom change trends of the patients in the two groups via heatmaps. Based on the identified and validated clinical biomarkers related to inflammation， oxidative stress， and energy metabolism in the early phase， Spearman correlation analysis was employed to analyze and evaluate the associations between clinical biomarkers and syndrome improvement. Key clinical biomarkers reflecting the effect of Huanglian Jiedu Wan were screened through the comparison of differences between groups. An extreme gradient boosting （XGBoost） algorithm was used to develop a prediction model for main symptom classification， with classification performance evaluated through 10-fold cross-validation. Feature importance analysis was applied to identify variables with the greatest contribution to the prediction result. ResultsThe syndrome transition matrix results indicated that the Huanglian Jiedu Wan group showed a superior effect to the placebo group in improving oral ulcers， sore throat， and overall symptoms， with significant effects observed especially in sore throat and overall symptom analyses （P<0.01）. Spearman correlation analysis revealed that several clinical biomarkers positively correlated with "excess heat-toxicity" syndrome and its main symptom improvement， were also called "heat-related biomarkers"， including succinic acid， α-ketoglutaric acid， glycine， lactic acid， adenosine monophosphate （AMP）， tumor necrosis factor-α （TNF-α）， interferon-γ （IFN-γ）， interleukin-1β （IL-1β）， interleukin-4 （IL-4）， interleukin-6 （IL-6）， interleukin-8 （IL-8）， interleukin-10 （IL-10）， and so on. Conversely， clinical biomarkers negatively correlated with symptom severity， were also called "heat-clearing related biomarkers" after administration of Huanglian Jiedu Wan， including malic acid， fumaric acid， cis-aconitic acid， adrenocorticotropic hormone （ACTH）， IL-1β， IL-4， IL-8， succinic acid， and citric acid. The XGBoost classification model using all 52 biomarkers as variables achieved an average test accuracy of 0.754 and an average F1 score of 0.777. Feature importance analysis identified the scores of glutamic acid in saliva and IL-6 were the highest in all the variables， with importance scores of 0.081 and 0.080， respectively. After screening out 14 key variables and optimizing the parameters， model performance improved to an average accuracy of 0.758 and an F1 score of 0.798. Feature importance analysis further determined that the glutamic acid in saliva and IL-6 showed obvious changes after screening the variables， confirming the good syndrome prediction ability of the model constructed by these key clinical biomarkers. ConclusionThis study systematically elucidates the correlation between syndrome improvement and clinical biomarkers of Huanglian Jiedu Wan in the treatment of "excess heat-toxicity" syndrome. An XGBoost classification model based on key clinical biomarkers is successfully established， achieving effective prediction of the symptoms related to the "excess heat-toxicity" syndrome such as oral ulcers and sore throat and providing a new insight for objective identification of traditional Chinese medicine syndromes.

ObjectiveThis paper aims to find the identified and validated clinical biomarker data building upon a clinical study of early-phase phase Ⅱ and investigate the correlation analysis of Huanglian Jiedu Wan on syndrome improvement and clinical biomarkers in the treatment of "excess heat-toxicity" based on a machine learning model. Additionally， the effective prediction of clinical biomarker values for the main symptoms of the "excess heat-toxicity" syndrome was assessed. MethodsA total of 229 patients meeting the inclusion criteria for "excess heat-toxicity" syndrome were randomly divided into the Huanglian Jiedu Wan group and the placebo group. Syndrome score transition matrices were constructed for the Huanglian Jiedu Wan group and the placebo group based on three main symptoms of "excess heat-toxicity" syndrome， such as oral ulcers， sore throat， and gum swelling and pain. Data from the patients with these three syndromes were also integrated for an overall analysis. The corresponding syndrome score transition matrices were further constructed to visualize symptom change trends of the patients in the two groups via heatmaps. Based on the identified and validated clinical biomarkers related to inflammation， oxidative stress， and energy metabolism in the early phase， Spearman correlation analysis was employed to analyze and evaluate the associations between clinical biomarkers and syndrome improvement. Key clinical biomarkers reflecting the effect of Huanglian Jiedu Wan were screened through the comparison of differences between groups. An extreme gradient boosting （XGBoost） algorithm was used to develop a prediction model for main symptom classification， with classification performance evaluated through 10-fold cross-validation. Feature importance analysis was applied to identify variables with the greatest contribution to the prediction result. ResultsThe syndrome transition matrix results indicated that the Huanglian Jiedu Wan group showed a superior effect to the placebo group in improving oral ulcers， sore throat， and overall symptoms， with significant effects observed especially in sore throat and overall symptom analyses （P<0.01）. Spearman correlation analysis revealed that several clinical biomarkers positively correlated with "excess heat-toxicity" syndrome and its main symptom improvement， were also called "heat-related biomarkers"， including succinic acid， α-ketoglutaric acid， glycine， lactic acid， adenosine monophosphate （AMP）， tumor necrosis factor-α （TNF-α）， interferon-γ （IFN-γ）， interleukin-1β （IL-1β）， interleukin-4 （IL-4）， interleukin-6 （IL-6）， interleukin-8 （IL-8）， interleukin-10 （IL-10）， and so on. Conversely， clinical biomarkers negatively correlated with symptom severity， were also called "heat-clearing related biomarkers" after administration of Huanglian Jiedu Wan， including malic acid， fumaric acid， cis-aconitic acid， adrenocorticotropic hormone （ACTH）， IL-1β， IL-4， IL-8， succinic acid， and citric acid. The XGBoost classification model using all 52 biomarkers as variables achieved an average test accuracy of 0.754 and an average F1 score of 0.777. Feature importance analysis identified the scores of glutamic acid in saliva and IL-6 were the highest in all the variables， with importance scores of 0.081 and 0.080， respectively. After screening out 14 key variables and optimizing the parameters， model performance improved to an average accuracy of 0.758 and an F1 score of 0.798. Feature importance analysis further determined that the glutamic acid in saliva and IL-6 showed obvious changes after screening the variables， confirming the good syndrome prediction ability of the model constructed by these key clinical biomarkers. ConclusionThis study systematically elucidates the correlation between syndrome improvement and clinical biomarkers of Huanglian Jiedu Wan in the treatment of "excess heat-toxicity" syndrome. An XGBoost classification model based on key clinical biomarkers is successfully established， achieving effective prediction of the symptoms related to the "excess heat-toxicity" syndrome such as oral ulcers and sore throat and providing a new insight for objective identification of traditional Chinese medicine syndromes.

ObjectiveTo observe the short-term and long-term efficacy of traditional Chinese medicine （TCM） surgical operations in treating mixed hemorrhoids. MethodsA multi-center retrospective cohort study was conducted， collecting clinical data from 17，831 mixed hemorrhoid surgery patients in 8 top-tier TCM hospitals in Jiangsu Province. Standardized and structured datasets were obtained through artificial intelligence models. Patients who underwent western surgical treatment were categorized into the western surgery group （11,646 cases）， and those receiving TCM surgical operations were categorized into the TCM surgery group （6185 cases）. Propensity score matching （1∶1 matching） was used to balance baseline data between groups. The primary outcome was the one-year recurrence rate， and secondary outcomes included the main symptoms （rectal bleeding， degree of prolapse） and secondary symptoms （anal distension， anal edema， wound secretion and exudation， anal stenosis， residual skin tags， perianal itching， and anal pain） measured on days 7， 28， and 60 after discharge. ResultsAfter matching， 2194 patients were included in each group. Symptom scores showed that at 28 days after discharge， the TCM surgical group had superior improvement in rectal bleeding ［OR=5.786， 95%CI （3.092，10.827）］， degree of prolapse ［OR=4.510， 95%CI （1.649，12.333）］， and anal edema ［OR=3.188， 95%CI （1.295，7.845）］ compared to the western surgical group. At 60 days post-discharge， the TCM group still showed advantages in improving rectal bleeding ［OR=5.237， 95%CI （1.077，25.464）］ and anal pain ［OR=11.697， 95%CI （1.186，115.336）］ （P＜0.05）. Long-term follow-up showed that the one-year recurrence rate in the TCM surgery group was 1.1% （8/726）， while that in the western surgery group was 2.3% （10/444）， with no statistically significant difference between the two groups （P＞0.05）. ConclusionBased on real-world data， TCM surgical treatment for mixed hemorrhoids shows significant short-term symptom improvement， particularly in terms of hemostasis， reducing swelling， and alleviating prolapse of anal masses.

OBJECTIVE To explore the ameliorative effect and potential mechanism of vitexin on inflammation in ulcerative colitis （UC） mice. METHODS The UC mice model was established by continuous administration of 3% dextran sulfate sodium solution for 5 days. Mice with successful modeling were randomly divided into UC group， vitexin low- and high-dose groups （vitexin-L and vitexin-H groups， 40， 80 mg/kg）， mesalazine group （400 mg/kg）， and vitexin-H+recombinant Jagged canonical Notch ligand 1 （rJagged-1） group （vitexin-H+rJagged-1 group， 80 mg/kg vitexin+1 mg/kg rJagged-1）， with 12 mice in each group. Another 12 normal mice were used as the control （CK） group. Mice in each group were administered the corresponding drugs or the corresponding drugs and normal saline by gavage and intraperitoneal injection once daily for 7 consecutive days. General conditions were observed during the experiment. At 24 h after the last administration， the disease activity index （DAI） score was evaluated. Colonic histopathological morphology was observed and scored. Macrophage polarization levels in the spleen and colon tissues were measured. The protein expressions of interleukin-6 （IL-6）， IL-10， tumor necrosis factor-α （TNF-α）， transforming growth factor-β 1 （TGF-β 1 ）， Jagged-1， Notch1 and Notch intracellular domain （NICD） in colonic tissues were determined. RESULTS Compared with the UC group， the symptoms （reduced food and water intake， dull fur， etc.） and pathological changes （epithelial cell shedding， inflammatory cell infiltration， etc.） were significantly improved in the vitexin-L， vitexin-H and mesalazine groups. DAI scores， colonic histopathological scores， M1 macrophage contents in spleen tissue， M1/M2 macrophage ratios， M1 macrophage proportions in colon tissue， and protein expressions of IL-6， TNF-α， Jagged-1， Notch1 and NICD in colon tissue were significantly decreased （ P ＜0.05）. Meanwhile， the M2 macrophage contents in spleen tissue， M2 macrophage proportions in colon tissue， and protein expressions of IL-10 and TGF-β 1 in colon tissue were significantly increased （ P ＜0.05）. Moreover， the improvement effects in the vitexin-H and mesalazine groups were significantly superior to those in the vitexin-L group （ P ＜0.05）. Compared with the vitexin-H group， the above symptoms and pathological changes were aggravated， and all quantitative indicators were significantly reversed in the vitexin-H+rJagged-1 group （ P ＜0.05）. CONCLUSIONS Vitexin can ameliorate the inflammation of UC mice， which is associated with its inhibition of the Jagged-1/Notch1 pathway and regulation of macrophage polarization （inhibition of M1-type polarization and promotion of M2-type polarization）.

Background Prolonged vibration exposure can lead to vascular endothelial cell dysfunction and cellular injury. However, research on the association between vibration and ferroptosis in vascular endothelial cells remains insufficient. Objective To explore whether occupational vibration exposure is associated with alterations in serum markers related to ferroptosis in patients with hand-arm vibration disease (HAVD), and to further investigate, through in vitro cell experiments, whether vibration exposure may induce ferroptosis in vascular endothelial cells. Methods ①A judgmental sampling method was employed to select 50 workers with HAVD (the HAVD group), 50 vibration-exposed workers without HAVD (the vibration exposure group), and 50 non–hand-transmitted vibration-exposed workers (the control group). Serum iron levels, malondialdehyde (MDA) content, and superoxide dismutase (SOD) levels were measured using serum iron assay kits, MDA detection kits, and SOD detection kits, respectively. One-way analysis of variance and binary logistic regression analysis were performed to examine the relationships between these indicators and HAVD. ②Human umbilical vein endothelial cells (HUVEC) were divided into a vibration group and a control group. The vibration group was subjected to vibration at 120 Hz with an acceleration of 6.5 m·s⁻² and further subdivided into four subgroups: 1 d 2 h, 1 d 4 h, 2 d 2 h, and 2 d 4 h. The control group was treated identically except for vibration exposure. Cellular iron (Fe²⁺) content and reduced glutathione (GSH) levels in HUVEC were measured using ferrous iron colorimetric assay kits and GSH colorimetric assay kits, respectively, to assess the effects of different vibration exposure schedules. Real-time quantitative polymerase chain reaction (RT-qPCR) was performed to detect the mRNA expression levels of ferroptosis-related genes, including acyl-CoA synthetase long-chain family member 4 (ACSL4), tumor suppressor protein P53 (P53), ferritin heavy chain 1 (FTH1), and glutathione peroxidase 4 (GPX4). Western blot analysis was conducted to determine the protein expression levels of ferroptosis-related markers in HUVEC. Results ①Compared with the control group, the patients in the HAVD group showed increased serum iron and MDA levels, along with decreased SOD levels (P<0.05). The logistic regression analysis indicated that elevated serum iron levels were significantly associated with an increased risk of HAVD (OR=4.034; 95%CI: 2.063, 7.887), and elevated MDA levels were also associated with an increased risk of HAVD (OR=1.523; 95%CI: 1.026, 1.936). ②Compared with the control group, increased intracellular Fe²⁺ content and decreased GSH content were observed in HUVECs in the 1 d 4 h and 2 d 4 h vibration subgroups (P<0.05). The RT-qPCR results showed that, compared with the control group, vibration exposures of 1 d 4 h and 2 d 4 h significantly upregulated the expression of ACSL4 and P53 (P<0.05), whereas the mRNA expression levels of GPX4 and FTH1 were downregulated in all vibration-exposed endothelial cells (P<0.05). The Western blot results revealed that, compared with the control group, the vibration exposure schedules of 1 d 2 h and 1 d 4 h significantly upregulated the protein expression levels of ACSL4 and P53 (P<0.05), while the vibration exposure schedules of 1 d 4 h, 2 d 2 h, and 2 d 4 h significantly downregulated the protein expression levels of FTH1 and GPX4 (P<0.05). Conclusion Occupational vibration exposure is associated with alterations in iron metabolism and oxidative stress status in workers with HAVD. The in vitro experiments further demonstrates that vibration stimulation induces intracellular iron accumulation and reduces antioxidant capacity in vascular endothelial cells, accompanied by dysregulated expression of ferroptosis-related molecules. These findings suggest that ferroptosis may play a role in vibration-induced vascular injury and the pathogenesis of HAVD.

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.

Objective To investigate the current distribution and predict the future suitable habitats of Haemaphysalis longicornis in Nanjing City, so as to provide insights into control and early warning of ticks and management of tick-borne diseases in Nanjing City. Methods The electronic map of Nanjing City was obtained from the National Platform for Common GeoSpatial Information Services. The distribution of H. longicornis and the longitude and latitude of distribution points from 2022 to 2024 were obtained from centers for disease control and prevention across each district in Nanjing City. Climatic and environmental variable data in Nanjing City were captured from the Worldclim database. Initially, 19 bioclimatic variables in this database were selected, including annual mean temperature, mean diurnal range, isothermality, temperature seasonality, maximum temperature of the warmest month, minimum temperature of the warmest month, temperature annual range, mean temperature of the wettest quarter, mean temperature of the driest quarter, mean temperature of the warmest quarter, mean temperature of the coldest quarter, annual precipitation, precipitation of the wettest month, precipitation of the driest month, precipitation seasonality, precipitation of the wettest quarter, precipitation of the driest quarter, precipitation of the warmest quarter, and precipitation of the coldest quarter. The elevation and normalized difference vegetation index were obtained from Data Sharing Platform of the Center for Resources and Environmental Sciences, Chinese Academy of Sciences. Then, the distribution points of H. longicornis, elevation, vegetation index and 19 bioclimatic variables were loaded into the software MaxEnt 3.4.4 to evaluate and screen out the variables with a contribution rate of 1% and higher. ArcGIS 10.8.1 software was used to extract the elevation, vegetation index and 19 bioclimatic variables of the distribution points of H. longicornis for a correlation analysis. If the absolute value of the correlation coefficient was 0.8 and higher, the variable with the higher contribution was retained. The 2050 dataset of the BCCCSM2-MR atmospheric circulation model in the coupled model intercomparison project phase 6 (CMIP6) were obtained from the Worldclim database as climate data for 2050. Screened H. longicornis species data and environmental and climate data were loaded into the maximum entropy (MaxEnt) model with the software MaxEnt 3.4.4 for training and validation, and then, all data generated from the model were imported into the software ArcGIS 10.8.1 to generate raster data and yield the map pertaining to the distribution of H. longicornis risk in Nanjing City. The accuracy of the model was evaluated with a receiver operating characteristic (ROC) curve, and the predictive effect of the model was assessed with area under the ROC curve (AUC). The suitable habitats of H. longicornis were classified in Nanjing City with the software ArcGIS 10.8.1, and the areas of distribution of suitable habitats in various categories were recorded to create the map of current H. longicornis suitable habitats classification in Nanjing City. The climatic and geographic information data in 2050 were employed as future environmental and climatic factors, and current environmental and climatic factors and current H. longicornis distribution data were additionally used to predict the future suitable habitats of H. longicornis in Nanjing City. In addition, the contributions of environmental and climatic factors to distribution of suitable habitats of H. longicornis was evaluated with the Jackknife method in Nanjing City. Results A total of 10 environmental and climatic variables were screened for analysis of the suitability of H. longicornis in Nanjing City based on correlation analyses and contributions of the MaxEnt model, including annual mean temperature, precipitation of the warmest quarter, vegetation index, precipitation of the wettest month, temperature annual range, annual precipitation, mean temperature of the warmest quarter, elevation, mean temperature of the wettest quarter, and maximum temperature of the warmest month, and annual mean temperature (34.8%), precipitation of the warmest quarter (17.3%), vegetation index (13.1%), and precipitation of the wettest month (10.8%) contributed relatively highly to the distribution of suitable habitats of H. longicornis in Nanjing City. The mean AUC of the ROC curve was 0.810 ± 0.055 for 10 repeated modeling results of the MaxEnt model, indicating high predictive performance of the model. The potential distribution areas of H. longicornis were predicted to be mainly located in Luhe District, Pukou District, Jiangning District, Lishui District, and Gaochun District in Nanjing City with the MaxEnt model. Under current climatic conditions, the area of potential suitable habitats of H. longicornis was 4 182.42 km² in Nanjing City, including 1 252.94 km² highly suitable habitats, which accounted for 19.00% of the total area of Nanjing City. Under the climate scenario in 2050, the area of potential suitable habitats of H. longicornis was projected to increase to 5 467.58 km² in Nanjing City, accounting for 82.95% of the total area of the city, and these habitats were mainly concentrated in Luhe District, Pukou District, Jiangning District, Lishui District, and Gaochun District. The areas of suitable habitats of H. longicornis at various categories were predicted to vary greatly in 2050, and the area of highly suitable habitats of H. longicornis was projected to increase to 2 378.82 km², accounting for 36.08% of the total area of Nanjing City. Based on jackknife tests and contributions of environmental and climatic variables, 6 dominant environmental and climatic factors were screened, including annual mean temperature (34.8% contribution), precipitation of the warmest quarter (17.3% contribution), vegetation index (13.1% contribution), precipitation of the wettest month (10.8% contribution), temperature annual range (5.4% contribution), and mean temperature of the warmest quarter (5.0% contribution), with cumulative contributions of 86.4%. Conclusion The distribution of H. longicornis is strongly associated with vegetation, temperature and precipitation in Nanjing City. Future climate change may lead to an expansion of the distribution area of H. longicornis in Nanjing City.

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.