ObjectiveTo evaluate the antitumor activity of Periplaneta americana extract（PAE） against human-derived lung adenocarcinoma organoids（LUAD-PDOs） and to elucidate its potential mechanism based on transcriptomics. MethodsFresh tumor and adjacent normal tissues from patients with LUAD were collected to construct LUAD-PDOs and normal lung organoid（Nor-PDOs） models using 3D organoid culture technology. The effective intervention concentration of PAE was determined using the cell counting kit-8（CCK-8） assay. Experimental groups included the model group（LUAD-PDOs）， normal group， model administration group（LUAD-PDOs+PAE）， and normal administration group（Nor-PDOs+PAE）. Hematoxylin-eosin（HE） staining was used to observe the pathological structures of PDOs， immunohistochemistry（IHC） was performed to detect the expressions of the proliferation marker Ki-67 and lung adenocarcinoma differentiation markers cytokeratin-7（CK-7） and Napsin A， TUNEL staining was applied to detect cell apoptosis. RNA sequencing（RNA-Seq） was conducted to identify differentially expressed genes（DEGs）， followed by Gene Ontology（GO）， Kyoto Encyclopedia of Genes and Genomes（KEGG）， and Gene Set Enrichment Analysis（GSEA）， alongside protein-protein interaction（PPI） network analysis to screen core mechanisms. Finally， key targets were validated by integrating external database analysis with immunofluorescence（IF）. ResultsNor-PDOs and LUAD-PDOs that highly recapitulated the pathological characteristics of the primary tissues were successfully established. The CCK-8 assay determined that the effective intervention concentration of PAE was 16 g·L^-1. Morphological observation showed that Nor-PDOs exhibited lumen-forming structures， whereas LUAD-PDOs displayed dense， solid structures. CCK-8 and TUNEL assays revealed that， compared with the model group， PAE intervention inhibited the proliferation of LUAD-PDOs and promoted apoptosis in LUAD cells， while showing no significant effect on the viability of Nor-PDOs. Transcriptomic analysis identified 719 DEGs that were significantly reversed after PAE intervention（347 up-regulated and 372 down-regulated）（P<0.05）. GO enrichment analysis indicated that DEGs in the model administration group were significantly enriched in biological processes related to cell cycle regulation compared to the model group. KEGG pathway analysis revealed that PAE affected pathways related to proliferation and metabolism， including pathways in cancer and the p53 signaling pathway. GSEA further confirmed that PAE significantly enhanced the activity of the p53 signaling pathway（P<0.05）. PPI network analysis indicated that breast cancer type 1 susceptibility protein（BRCA1） and checkpoint kinase 1（CHEK1） were the core down-regulated targets in the p53 pathway. IF verified the high expression of BRCA1 and CHEK1 in LUAD-PDOs and their significant downregulation after PAE intervention（P<0.05）. Furthermore， survival analysis based on The Cancer Genome Atlas（TCGA） database indicated that low expression of BRCA1 and CHEK1 was significantly associated with prolonged overall survival in patients with LUAD（P<0.05）. ConclusionPAE effectively inhibits proliferation of LUAD-PDOs and promotes their apoptosis， its anti-tumor mechanism is potentially associated with the activation of the p53 signaling pathway， with BRCA1 and CHEK1 genes likely serving as key downstream targets for the effects of PAE.

Objective To assess the survival vulnerability of Oncomelania hupensis in Jiangxi Province under future climate scenarios, and to identify low-vulnerability areas for its survival in this province. Methods Village-level O. hupensis snail survey and O. hupensis snail control with chemical treatments in Jiangxi Province from 2016 to 2024 were captured from the Parasitic Disease Prevention and Control Information Management System of China Disease Prevention and Control Information System. Climatic data were primarily sourced from the Resource and Environmental Science Data Platform, Chinese Academy of Sciences (http://www.resdc.cn/), including annual average temperature, annual average precipitation, annual accumulated temperature above 10 °C, annual accumulated temperature above 0 °C, annual maximum temperature, annual minimum temperature, and annual average relative humidity, and nineteen bioclimatic variables were downloaded from the WorldClim website (https://www.worldclim.org/), including mean diurnal range, isothermality, temperature seasonality, and so on. Elevation and normalized difference vegetation index were catprued from the Resource and Environmental Science Data Platform, Chinese Academy of Sciences (http://www.resdc.cn/), and distance to rivers was downloaded from the WorldPop website (http://www.worldpop.org), and land use and land cover (LULC) data were downloaded from the Big Earth Data Center, Chinese Academy of Sciences (https://data.casearth.cn/), and nature reserve data were obtained from the China Nature Reserve Specimen Resource Sharing Platform (http://www.papc.cn/). Three Shared Socioeconomic Pathways (SSPs) from the Beijing Climate Center-Climate System Model version 2-Medium Resolution (BCC-CSM2-MR) global climate model were employed as future climate scenarios, including SSP126, SSP245, SSP585, and the biomod2 ensemble model in R package was used to simulate suitable habitats for O. hupensis snails in Jiangxi Province in 2050 and 2070 under these scenarios. A snail survival vulnerability index was constructed based on the area of suitable snail habitats, area covered by snail control through chemical treatment, area covered by nature reserves, and changes in snail habitat fragmentation, and a map of snail survival vulnerability distribution was plotted. Results The real area of snail habitats ranged from 78 486.76 to 85 309.47 hm², and the area of snail control with chemical treatment ranged from 10 138.98 to 13 240.16 hm² in Jiangxi Province from 2016 to 2024. There were 429 to 531 villages detected with snails during the nine-year period, and the number of actually snail-infested villages ranged from 645 to 686. A total of 818 snail-present points and 1 996 snail-absent points were obtained from snail survey records. The best performance of the biomod2 ensemble model was achieved if a weighted mean approach was used as the ensemble strategy, with a true skill statistic value of 0.799 and an area under the receiver operating characteristic curve of 0.957, and modeling identified annual average relative humidity and annual average precipitation as two most influencing climatic variables for snail distribution. Relative to the current areas of suitable snail habitats under present climate conditions, the area of suitable snail habitats was projected to expand by 24.49% to 46.28% in Jiangxi Province under future climate scenarios, and the proportion of nature reserves areas in the areas of suitable snail habitats was projected to decrease slightly from the current 2.77% to approximately 2.52%, while the proportion of areas of snail control through chemical treatment in areas of suitable snail habitats varied from 0.64% to 19.57%, and the percentage of changes in snail habitat fragmentation ranged from 3.86% to 12.23%. Based on these four indicators, the snail survival vulnerability index was estimated to range from –1.96 to 0.62 in Jiangxi Province. The arithmetic mean of the snail survival vulnerability index differed under three SSP scenarios (SSP126, SSP245 and SSP585), with the highest mean value (–0.69) in 2070 under SSP126, and the lowest mean value (–0.78) in 2070 under SSP585. Conclusions The snail survival vulnerability index ranges from –1.96 to 0.62 in Jiangxi Province under future climate scenarios, and the suitable habitats for O. hupensis snails appear an overall tendency towards expansion. Low-vulnerability snail habitats are mainly distributed along the shores of Poyang Lake and the Yangtze River in Jiangxi Province, partially overlapping with nature reserves. Intensified surveillance of O. hupensis snails is recommended in these areas in the future.

Objective To reveal the scientific output and trends in pulmonary nodules/early-stage lung cancer prediction models. Methods Publications on predictive models of pulmonary nodules/early lung cancer between January 1, 2002 and June 3, 2023 were retrieved and extracted from CNKI, Wanfang, VIP and Web of Science database. CiteSpace 6.1.R3 and VOSviewer 1.6.18 were used to analyze the hotspots and theme trends. Results A marked increase in the number of publications related to pulmonary nodules/early-stage lung cancer prediction models was observed. A total of 12581 authors from 2711 institutions in 64 countries/regions published 2139 documents in 566 academic journals in English. A total of 282 articles from 1256 authors were published in 176 journals in Chinese. The Chinese and English journals which published the most pulmonary nodules/early-stage lung cancer prediction model-related papers were Journal of Clinical Radiology and Frontiers in Oncology, respectively. Chest was the most frequently cited journal. China and the United States were the leading countries in the field of pulmonary nodules/early-stage lung cancer prediction models. The institutions represented by Fudan University had significant academic influence in the field. Analysis of keywords revealed that multi-omics, nomogram, machine learning and artificial intelligence were the current focus of research. Conclusion Over the last two decades, research on risk-prediction models for pulmonary nodules/early-stage lung cancer has attracted increasing attention. Prognosis, machine learning, artificial intelligence, nomogram, and multi-omics technologies are both current hotspots and future trends in this field. In the future, in-depth explorations using different omics should increase the sensitivity and accuracy of pulmonary nodules/early-stage lung cancer prediction models. More high-quality future studies should be conducted to validate the efficacy and safety of pulmonary nodules/early-stage lung cancer prediction models further and reduce the global burden of lung cancer.

Objective:To investigate the treatment and prognosis of esophageal cancer patients with pulmonary resection history.Methods:The retrospective and descriptive study was conducted. The clinicopathological data of 58 esophageal cancer patients with pulmonary resection history who were admitted to Chest Hospital Affiliated to Shanghai Jiaotong University School of Medicine and Jiangxi Provincial People's Hospital from May 2019 to April 2024 were collected. There were 52 males and 6 females, aged (69±3)years. Observation indicators: (1) surgical and postopera-tive conditions; (2) postoperative pathological examination results; (3) follow-up; (4) stratified analysis. Comparison of measurement data with normal distribution between groups was conducted using the independent sample t test. Comparison of count data between groups was conducted using the chi-square test or Fisher exact probability. Comparison of ordinal data between groups was conducted using the non-parametric rank sum test. The Kaplan-Meier method was used to plot survival curve and calculate survival rate, and the Log-rank test was used for survival analysis. Results:(1) Surgical and postoperative conditions. Of the 58 esophageal cancer patients, 49 patients underwent transthoracic approach (26 cases of ipsilateral approach and 23 cases of contralateral approach of pulmonary resection history), and 9 patients underwent mediastinoscopic-laparoscopic approach. There were 57 cases with R 0 resection and 1 case with R 2 resection because of tumor invading carina. The total operation time of 58 patients was (246±27)minutes, and the volume of intraoperative blood loss was (114±29)mL. There was no unplanned reoperation or perioperative death for all patients. The duration of postoperative hospital stay of 58 patients was (10.4±4.6)days, and time for intensive care unit stay was (1.4±0.5)days, and no patient readmitted to intensive care unit due to changes in conditions. The postoperative total incidence of complications of 58 patients was 41.4%(24/58). The Clavien-Dindo grading of complications for all patients was 1-2 grade. (2) Postoperative pathological examination results. Results of postoperative pathological examination showed there were 51 cases of squamous cell carcinoma, 6 cases of adenocarcinoma, and 1 case of melanoma. Number of lymph node dissected of 58 patients was 27±6. The ratio of patient with positive lymph node was 37.9%(22/58). One patient may experience more than 1 region of positive lymph node metastasis. Results of postoperative pathological staging showed 5 cases of ⅠA stage, 2 cases of ⅠB stage, 13 cases of ⅡA stage, 15 cases of ⅡB stage, 4 cases of ⅢA stage, 16 cases of ⅢB stage, and 3 cases of ⅣA stage. Thirteen of the 58 patients underwent neoadjuvant therapy, with the pathological staging as 6 cases of Ⅰ stage, 4 cases of Ⅱ stage, 3 cases of ⅢB stage after therapy. Results of postoperative tumor regression grade for the 13 patients with neoadjuvant therapy showed 4 cases of grad 0, 3 cases of grade 1, 6 cases of grade 2. (3) Follow-up. All 58 patients were followed for 24 (4, 50)months, and no patient died within 90 days after surgery. During the follow-up period, 19 patients experienced tumor recurrence and metastasis and 17 patients died. Twenty-one patients underwent postoperative adjuvant therapy, including 7 cases with chemoradiotherapy, 7 cases with chemotherapy, 3 cases with chemotherapy and immunotherapy, 2 cases with immuno-therapy, 2 cases with radiotherapy. The postoperative 1-, 2-year overall survival rates of the 58 patients were 91.3%, 78.7%, respectively, of whom undergoing McKeown surgery and mediastinoscopic-laparoscopic surgery with postoperative 1-, 2-year overall survival rates as 89.2%, 83.1% and 85.7%, 53.6%, respectively. The postoperative 1-, 2-year esophageal cancer specific survival rates for patients undergoing McKeown surgery and mediastinoscopic-laparoscopic surgery were 94.4%, 87.9% and 85.7%, 71.4%, respectively. There was no significant difference in postoperative 1-, 2-year overall survival rates and postoperative 1-, 2-year esophageal cancer specific survival rates between patients undergoing McKeown surgery and mediastinoscopic-laparoscopic surgery ( P>0.05). (4) Stratified analysis. Of the 49 patients underwent transthoracic approach for esophageal cancer, there were significant differences in surgical method, surgical type, time of chest surgery, cases with upper mediastinal lymph node dissection, and duration of postoperative hospital stay between patients with pulmonary resection history as ipsilateral approach and contralateral approach ( χ2=11.74, 11.68, t=-2.25, χ2=8.45, t=-2.17, P<0.05), and there was no significant difference in total operation time, volume of intraoperative blood loss, the number of lymph node dissected, post-operative total complications, and postoperative pathological TNM staging ( P>0.05). For patients with pulmonary resection history as ipsilateral approach and contralateral approach, the postopera-tive 1-, 2-year esophageal cancer specific survival rates were 95.5%, 95.5% and 81.4%, 71.1%, showing a significant difference between them ( χ2=5.63, P<0.05). Conclusions:The transthoracic approach and mediastinoscopic-laparoscopic approach are safe and feasible for esophageal cancer patients with pulmonary resection history. Compared with patients with pulmonary resection history as contralateral approach, patients with pulmonary resection history as ipsilateral approach have a higher ratio of McKeown surgery, minimally invasive surgery and upper mediastinal lymph node dissection, shorter time of chest surgery and duration of postoperative hospital stay, better esophageal cancer specific survival rate. And there is no increase in perioperative risk.

Objective:To explore the RNA expression and alterations in brain structure in individuals who have experienced stressful life events (SLE), as well as the correlation patterns between them and their association with the occurrence of depression.Methods:Prospectively, a total of 80 SLE subjects were recruited from the psychiatry and psychology clinic of the Jiangsu University Affiliated Yixing Hospital between January 2021 and December 2022, with 16 normal controls (NC) enrolled concurrently. The 17 items Hamilton depression scale (HAMD-17) and social readjustment rating scale (SRRS) were used to assess depressive symptoms and stress levels. RNA sequencing information of peripheral blood and imaging data at baseline were collected. Based on whether depression occurred during the 2-year follow-up period, SLE subjects were divided into the SLE-depression group ( n=15) and the SLE-non-depression group ( n=65). Differentially expressed genes (DEGs) were screened using differential analysis and protein-protein interaction (PPI) networks. Fractional anisotropy (FA) of white matter tracts and gray matter volume (GMV) were extracted using tract-based spatial statistics and voxel-based morphometry.Using analysis of variance compared inter-group differences in gene expression, GMV and white matter FA values. Partial correlation analysis was used to explore correlations between DEGs, altered GMV and white matter microstructure. Gene set enrichment analysis (GSEA) was performed on key genes to identify potential biological pathways. Propensity score matching constructed sensitivity subgroups to verify result robustness. Results:The SLE-depression group showed significantly higher SRRS and HAMD-17 scores at baseline and at the end of follow-up compared to the SLE-non-depression group and the NC group ( H=47.773, 35.427, 41.114, all P<0.05). Expression levels of IL-10 (2.12±0.28, 2.43±0.44), EZH2 (2.11±0.43, 2.45±0.51), NCAM1 (3.60±0.30, 3.03±0.39), CD3E (4.95±0.37, 4.57±0.48), CCK (3.29±0.28, 3.02±0.42), and CX3CR1 (5.55±0.40, 5.91±0.34) were significantly different between the SLE-depression group and SLE-non-depression group( F=5.549~28.371, all P<0.05). Compared with the SLE-non-depression group, the SLE-depression group exhibited significantly lower FA values in the genu of the corpus callosum (0.29±0.04, 0.31±0.04) and the left uncinate fasciculus (0.31±0.02, 0.33±0.02), as well as significantly smaller GMV in the right hippocampus (0.29±0.07, 0.33±0.06), bilateral middle frontal gyrus (left: 0.27±0.05, 0.31±0.05; right: 0.28±0.06, 0.32±0.06), right insula (0.36±0.03, 0.38±0.04), and left precentral gyrus (0.19±0.04, 0.24±0.05) ( F=4.593-12.064, all P<0.05, FDR correction). GMV in the right anterior cingulate and paracingulate gyri was significantly larger than that in the SLE-non-depression group (0.34±0.05, 0.29±0.06) ( F=6.704, P=0.034, FDR correction). Partial correlation analysis revealed significantly stronger correlations between hub DEGs and altered brain regions in the SLE-depression group ( r=0.017-0.801) compared to the SLE-non-depression group ( r=0.002-0.382), with a statistically significant difference ( U=629, P<0.001; Cliff's Delta=0.454). GSEA indicated that the aforementioned genes were primarily involved in pathways including the ribosome, spliceosome, ribosome biogenesis in eukaryotes, and neuroactive ligand-receptor interaction. Sensitivity analysis confirmed that the above results remained statistically significant after balancing sample sizes (all P<0.05). Conclusion:The SLE-depression group showed specific RNA expression and brain structure alterations compared to the SLE-non-depression group, and the correlation between RNA and brain structure was significantly enhanced in the SLE-depression group. This suggests that the correlation between genes and brain structure in the SLE population may be related to their susceptibility to depression.

In recent years, with the exploration and development of machine learning algorithms in the medical scenario, their potential for application in clinical diagnosis and treatment has been gradually explored. Image-based machine learning algorithms have made a lot of progress in urological cancer-related research, and have been proven to have a positive effect on the diagnosis, treatment and prognosis of urological cancer. Despite the great impact of machine learning algorithms clinically, they are still in the exploratory stage, with many deficiencies and limitations.

Objective To compare the difference,agreement,and axial length measurement success rate between biometers ZW-30 and IOLMaser 700 based on swept-source optical coherence tomography for the measurement of ocular bi-ological parameters in patients with cataracts.Methods A total of 126 cataract patients(233 eyes)who were advised to undergo cataract surgery at the Department of Ophthalmology at Beijing Tongren Hospital,Capital Medical University from January to February 2024 were included in this study.Two biometers were used to measure the axial length(AL),mean keratometry(Km),anterior chamber depth(ACD),lens thickness(LT),central corneal thickness(CCT),and horizontal corneal diameter[namely,the white-to-white(WTW)distance].The axial measurement success rate of the two biometers and the difference and agreement between the parameters were calculated.Results The mean difference between ZW-30 and IOLMaster 700 was(-0.006±0.042)mm for AL,(-0.074±0.204)D for Km,(0.031±0.051)mm for ACD,(0.001±0.005)mm for CCT,and(-0.286±0.337)mm for WTW,and the differences were statistically significant(all P＜0.05).The mean difference between ZW-30 and IOLMaster 700 was(0.008±0.215)mm for LT,and the difference was not statis-tically significant(t=0.579,P=0.563).The 95％limits of agreement range was between-0.011 mm and 0.000 mm for AL,between-0.474 D and 0.326 D for Km,between-0.010 mm and 0.012 mm for CCT,between-0.068 mm and 0.131 mm for ACD,between-0.116 mm and 0.159 mm for LT,and between-0.947 mm and 0.376 mm for WTW.The intra-class correlation coefficient of all measurements ranged from 0.790 to 1.000.The AL measurement success rate of IOLMaster 700 and ZW-30 was 95.3％and 95.7％,respectively.The latter had an AL measurement success rate of 98.7％after manually marking the position of the retinal identification line.Conclusion There were statistically significant differences between ZW-30 and IOLMaster 700 in the measurement of the AL,Km,ACD,and CCT,which,however,were not clinically significant.The agreement between both was good.ZW-30 had a higher AL measurement success rate,espe-cially for the manual identification function of eyes with opacified refractive media,which can further improve the AL meas-urement success rate and provide reference for clinical work.

BACKGROUND:Myocardial cell hypoxic injury is closely associated with oxidative stress and ferroptosis.Previous studies have shown that aloin has various effects such as antioxidant,anti-inflammatory,and anti-tumor activities.OBJECTIVE:To investigate the effects of aloin on oxidative stress and ferroptosis in hypoxia-induced H9C2 cells.METHODS:A hypoxia model was established using H9C2 myocardial cells.Firstly,cell viability was determined to confirm the lack of cytotoxicity of aloin and to determine its optimal therapeutic concentration.Subsequently,the effects of aloin on hypoxia-induced lactate dehydrogenase release,reactive oxygen species production,and mitochondrial oxidative stress in H9C2 cells were evaluated using assay kits,dihydroethidium fluorescent probes,and MitoSOX?Red fluorescent probes,respectively.To verify the effect of aloin on ferroptosis,intracellular Fe2+content and lipid peroxidation level were detected using fluorescence staining and flow cytometry,respectively.Then,the expression levels of ferroptosis regulatory factors glutathione peroxidase 4,acyl-CoA synthetase long-chain family member 4,and nuclear factor E2-related factor 2 were detected using western blot assay and real-time fluorescence quantitative PCR techniques.Finally,the role of ferroptosis in aloin-mediated myocardial protection was further confirmed by using the ferroptosis inducer Erastin.RESULTS AND CONCLUSION:(1)Compared with the control group,the viability of H9C2 cells in the hypoxia group was significantly decreased,lactate dehydrogenase release,reactive oxygen species level,mitochondrial oxidative stress degree,Fe2+content,and lipid peroxidation degree were significantly increased,while glutathione peroxidase 4 and nuclear factor E2-related factor 2 mRNA and protein expression levels were significantly decreased,and acyl-CoA synthetase long-chain family member 4 mRNA and protein expression were significantly increased(all P＜0.05).(2)Compared with the hypoxia group,both low and high doses of aloin reversed the changes in above indicators(all P＜0.05).(3)Compared with the hypoxia+aloin group,the hypoxia+aloin+Erastin group showed a significant decrease in H9C2 cell viability and a significant increase in lactate dehydrogenase release(both P＜0.01).The results indicate that aloin has a protective effect on hypoxia-treated H9C2 cells in a dose-dependent manner,mainly achieved by inhibiting oxidative stress and ferroptosis.

Porcine reproductive and respiratory syndrome virus(PRRSV)mainly causes sow abor-tion,stillbirth,mummified fetus and respiratory symptoms in piglets.Since first reported in China in 1996,the virus complexity has increased significantly in more than 20 years of genetic evolution,bringing huge economic losses to the pig industry.In recent years,with the emergence of various PRRSV recombinant virus strains,preventing and controlling this epidemic became increasingly difficult.The purpose of this article is to comprehensively review the genome structure and func-tion of PRRSV,RNA virus recombination mechanism,main types of recombination,and the epi-demic status and recombination for the dominant epidemic PRRSV strains,in order to provide clues for in-depth research on gene recombination of PRRSV,thus providing the theoretical sup-port for formulating scientific prevention and control strategies.

Objective:To explore the RNA expression and alterations in brain structure in individuals who have experienced stressful life events (SLE), as well as the correlation patterns between them and their association with the occurrence of depression.Methods:Prospectively, a total of 80 SLE subjects were recruited from the psychiatry and psychology clinic of the Jiangsu University Affiliated Yixing Hospital between January 2021 and December 2022, with 16 normal controls (NC) enrolled concurrently. The 17 items Hamilton depression scale (HAMD-17) and social readjustment rating scale (SRRS) were used to assess depressive symptoms and stress levels. RNA sequencing information of peripheral blood and imaging data at baseline were collected. Based on whether depression occurred during the 2-year follow-up period, SLE subjects were divided into the SLE-depression group ( n=15) and the SLE-non-depression group ( n=65). Differentially expressed genes (DEGs) were screened using differential analysis and protein-protein interaction (PPI) networks. Fractional anisotropy (FA) of white matter tracts and gray matter volume (GMV) were extracted using tract-based spatial statistics and voxel-based morphometry.Using analysis of variance compared inter-group differences in gene expression, GMV and white matter FA values. Partial correlation analysis was used to explore correlations between DEGs, altered GMV and white matter microstructure. Gene set enrichment analysis (GSEA) was performed on key genes to identify potential biological pathways. Propensity score matching constructed sensitivity subgroups to verify result robustness. Results:The SLE-depression group showed significantly higher SRRS and HAMD-17 scores at baseline and at the end of follow-up compared to the SLE-non-depression group and the NC group ( H=47.773, 35.427, 41.114, all P<0.05). Expression levels of IL-10 (2.12±0.28, 2.43±0.44), EZH2 (2.11±0.43, 2.45±0.51), NCAM1 (3.60±0.30, 3.03±0.39), CD3E (4.95±0.37, 4.57±0.48), CCK (3.29±0.28, 3.02±0.42), and CX3CR1 (5.55±0.40, 5.91±0.34) were significantly different between the SLE-depression group and SLE-non-depression group( F=5.549~28.371, all P<0.05). Compared with the SLE-non-depression group, the SLE-depression group exhibited significantly lower FA values in the genu of the corpus callosum (0.29±0.04, 0.31±0.04) and the left uncinate fasciculus (0.31±0.02, 0.33±0.02), as well as significantly smaller GMV in the right hippocampus (0.29±0.07, 0.33±0.06), bilateral middle frontal gyrus (left: 0.27±0.05, 0.31±0.05; right: 0.28±0.06, 0.32±0.06), right insula (0.36±0.03, 0.38±0.04), and left precentral gyrus (0.19±0.04, 0.24±0.05) ( F=4.593-12.064, all P<0.05, FDR correction). GMV in the right anterior cingulate and paracingulate gyri was significantly larger than that in the SLE-non-depression group (0.34±0.05, 0.29±0.06) ( F=6.704, P=0.034, FDR correction). Partial correlation analysis revealed significantly stronger correlations between hub DEGs and altered brain regions in the SLE-depression group ( r=0.017-0.801) compared to the SLE-non-depression group ( r=0.002-0.382), with a statistically significant difference ( U=629, P<0.001; Cliff's Delta=0.454). GSEA indicated that the aforementioned genes were primarily involved in pathways including the ribosome, spliceosome, ribosome biogenesis in eukaryotes, and neuroactive ligand-receptor interaction. Sensitivity analysis confirmed that the above results remained statistically significant after balancing sample sizes (all P<0.05). Conclusion:The SLE-depression group showed specific RNA expression and brain structure alterations compared to the SLE-non-depression group, and the correlation between RNA and brain structure was significantly enhanced in the SLE-depression group. This suggests that the correlation between genes and brain structure in the SLE population may be related to their susceptibility to depression.