Objective To investigate the heterogeneity and key molecular features of epithelial cells in triple-negative breast cancer (TNBC), identify prognostic biomarkers, and develop a robust survival prediction model. Methods Using TNBC single-cell transcriptomic data, epithelial cells were extracted, normalized, and subclustered to characterize their molecular signatures and functional differences. High-dimensional weighted gene co-expression network analysis (hdWGCNA) was applied to establish co-expression modules in epithelial cells. Multiple machine learning algorithms were integrated to select key prognostic genes and develop a risk-score model, whose performance was evaluated using receiver operating characteristic (ROC) curves and Kaplan-Meier (K-M) survival analysis. In addition, the immune microenvironment features and potential drug-response differences between the high- and low-risk groups were systematically assessed. Finally, PCR was performed to validate the expression differences of the key genes between tumor and normal tissues. Results We characterized the composition and molecular features of TNBC epithelial subpopulations and identified a TNBC-associated epithelial subset. By integrating hdWGCNA with machine learning approaches, 10 key genes were selected to construct a prognostic model, which effectively stratified patients into distinct survival-risk groups and demonstrated favorable predictive performance in ROC and K-M analyses. Immune profiling revealed the differences in the infiltration levels of seven immune cell types and immune function-related features between the high- and low-risk groups. Drug-sensitivity analysis suggested potential differential responses to eight agents across the risk groups. PCR validation further confirmed the differential expression of the ten signature genes between tumor and normal tissues. Conclusion This study reveals epithelial heterogeneity in TNBC at single-cell resolution and establishes a 10-gene prognostic model, which may facilitate the stratification of TNBC risk and the evaluation of immune characteristics and potential therapeutic strategies.

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

Objective This study reports on the obligatory parasitism of water mites Arrenurus sp. on Anopheles messeae at the Manzhouli Port, Inner Mongolia, China. Methods Duing July 2024, a survey on the mosquito diversity was conducted at the Manzhouli Port. Captured mosquitoes and their ectoparasites were identified to species level. Results A total of 1840 adult mosquitoes were collected, representing species from three genera: Culex(Cx. modestus, Cx. pipiens pallens), Aedes(Ae. dorsalis, Ae. flavidorsalis, Ae. flavescens), and Anopheles (An. messeae). Among all the mosqutioes specimens,3 out of 150 captured An. messeae were found to carry ectoparasitic mites, with number of 2,4,27 mites separately. Morphological and molecular identification reached the same result as water mites(Hydrachnidiae, Hydracrina). COI gene sequence showed 94% similarity with the closest species Arrenurus truncatellus. Conlusions Literature review suggests water mites are host-specific parasitism of mosquito species and herein with the first record of Arrenurus sp. parasiting on An. Messeae in the most high-latitude region globally.

Objective To prepare the recombinant Echinococcus granulosus Polo-like kinase 2 (rEgPLK2) protein and evaluate its immunoprotective efficacy against cystic echinococcosis, so as to provide insights into research and development of novel vaccines against echinococcosis. Methods The Polo-like kinase (PLK) protein sequences were retrieved from 12 species in the NCBI protein database, including E. granulosus and E. multilocularis. Multiple sequence alignment was performed using the Clustal Omega program, and structural visualization and homology analysis were conducted using the ESPript 3.2 program. The recombinant plasmid pET-30a-EgPLK2 was transformed into BL21(DE3) competent cells. Protein expression was induced with isopropyl-β-D-thiogalactoside (IPTG), and sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) was performed to characterize the expression and molecular weight of the rEgPLK2 protein. The purified rEgPLK2 protein was thoroughly emulsified with Freund’s complete adjuvant at a 1 : 1 volume ratio. Two New Zealand white rabbits were immunized with multipoint subcutaneous injection on the back at a dose of 300 μg per rabbit for primary immunization. For booster immunizations, the protein was emulsified with Freund’s incomplete adjuvant at a 1 : 1 volume ratio and administered on days 14, 28, and 42 after the primary immunization at a dose of 150 μg per rabbit. Serum was sampled from the rabbit ear vein on day 7 after the final immunization to yield anti-rEgPLK2 polyclonal antibodies. Antibody titer was determined by indirect enzyme-linked immunosorbent assay (ELISA), and antibody specificity was verified by Western blotting. The tissue localization of the EgPLK2 protein was detected in E. granulosus protoscoleces and adult worms using immunofluorescence assay (IFA). Eighteen 6- to 8-week-old female SPF-grade BALB/c mice were randomly divided into three groups, including the blank control group, rEgPLK2-ISA immunization group, and PBS-ISA adjuvant control group, of 6 mice each group. Mice in the rEgPLK2-ISA immunization group and PBSISA group received three primary immunizations via intramuscular injection, and animals in the rEgPLK2-ISA immunization group was inoculated with immunogens prepared by emulsifying rEgPLK2 protein with ISA 201 adjuvant at a 1 : 1 volume ratio (6 μg per mouse), while mice in the PBS-ISA adjuvant control group received an equal volume of PBS emulsified with ISA adjuvant at a 1 : 1 volume ratio. A fourth booster immunization was administered via intraperitoneal injection. Mice in the rEgPLK2-ISA immunization group received a booster immunization with 8 μg of rEgPLK2 protein per mouse, and animals in the PBS-ISA group received an equal volume of PBS, with immunizations given at 2-week intervals. Mice in the blank control group were given no treatment, and housed under standard conditions. Tail vein blood was collected from all mice 7 days after the final immunization, and levels of specific anti-rEgPLK2 IgG antibody and its subclasses (IgG1, IgG2a, IgG2b, IgG3) were measured by indirect ELISA. E. granulosus infection was modelled in mice through injection with 1 000 E. granulosus protoscoleces via intrahepatic portal vein in the rEgPLK2-ISA immunization group and PBS-ISA adjuvant control group 2 weeks after the last immunization. All mice were sacrificed and dissected. The number of cysts was counted in mouse livers, and the cyst reduction rate was calculated. Liver tissues were processed for paraffin sectioning and stained with hematoxylin and eosin (HE), and histopathological changes were examined under a light microscope. Results Sequence analysis revealed that EgPLK2 shared a high amino acid sequence homology with E. multilocularis PLK2 (EmPLK2) and contained the typical domains of the Polo-like kinase family, including the serine/threonine protein kinase catalytic domain (STKc) and Polo-box. The IPTG-induced rEgPLK2 protein was mainly expressed in the form of inclusion bodies, and the purified rEgPLK2 protein showed a relative molecular mass of approximately 70 kDa. The prepared rabbit anti-rEgPLK2 polyclonal antibody had a titer of 1 : 256 000, and Western blotting assay showed that this anti-body specifically recognized the rEgPLK2 protein with a relative molecular mass of approximately 70 kDa. Immunofluorescence assay showed that the EgPLK2 protein was localized in the excretory bladder and rostellum of E. granulosus protoscoleces, as well as the tegument, suckers, and inter-proglottid junctions of adult worms. Immunoprotective assay showed that the serum levels of specific anti-rEgPLK2 IgG, IgG1, IgG2a, and IgG2b antibodies were 2.92 ± 0.49, 0.33 ± 0.10, 0.31 (0.36), and 3.12 (1.73) in mice in the rEgPLK2-ISA immunization group, which were all significantly higher than those in the PBS-ISA adjuvant control group (0.14 ± 0.04, 0.07 ± 0.01, 0.12 ± 0.04, and 0.11 ± 0.04, respectively) (t = 19.28 and 8.46, Z = 3.75 and 4.15; all P values < 0.001); however, there was no significant difference in the serum anti-IgG3 antibody level between the rEgPLK2-ISA immunization group and the PBS-ISA adjuvant control group [0.07 (0.01) vs. 0.073 (0.07); Z = 0.69, P > 0.05)]. In the mouse model of E. granulosus infections, the area of hepatic lesions was reduced and the inflammatory infiltration was alleviated in the rEgPLK2-ISA immunization group than in the PBS-ISA adjuvant control group, and the number of hepatic cysts was higher in the PBS-ISA adjuvant control group than in the rEgPLK2-ISA immunization group [8.00 (2.00) vs. 1.00 (0.75); Z = −2.93, P < 0.01], with a cyst reduction rate of 80.40%. Indirect ELISA assay measured higher serum levels of specific anti-rEgPLK2 IgG (3.28 ± 0.48 vs. 0.11 ± 0.04; t = 15.86, P < 0.01), IgG1 (0.29 ± 0.02 vs. 0.09 ± 0.01; t = 15.67, P < 0.01), IgG2a [3.71 (1.09) vs. 0.08 (0.03); Z = 2.88, P < 0.01], and IgG2b antibodies [3.34 (1.01) vs. 0.08 (0.03); Z = 2.88, P < 0.01] in the rEgPLK2-ISA immunization group than in the PBS-ISA adjuvant control group, and there was no significant difference in the serum level of the specific anti-rEgPLK2 IgG3 antibody between the rEgPLK2-ISA immunization group and the PBS-ISA adjuvant control group (0.07 ± 0.01 vs. 0.07 ± 0.01; t = 1.29, P > 0.05). Conclusions The prokaryotic expression system has been successfully constructed for the EgPLK2 gene and the anti-rEgPLK2 polyclonal antibody has been obtained. The rEgPLK2 protein exhibits a high immunogenicity, and is effective to protect against E. granulosus infection, and inhibits cyst development, which is a promising candidate vaccine target against cystic echinococcosis.

ObjectiveTo explore the clinical and phenomics characteristics of patients with phlegm-stasis binding syndrome and its accompanied patterns in stable angina pectoris （SAP） of coronary heart disease （CHD）. MethodsA multicenter cross-sectional study design was adopted. A total of 300 patients with SAP of CHD were enrolled and classified into 120 cases of phlegm-stasis binding syndrome， 125 cases of qi deficiency-accompanied syndrome， 38 cases of qi stagnation-accompanied syndrome， and 17 cases of toxin accumulation-accompanied syndrome according to traditional Chinese medicine （TCM） patterns. Data of patients with different TCM patterns were collected， including general condition， TCM symptoms， blood lipids， coagulation function， immune indicators， and serum metabolomics. Metabolic pathway enrichment analysis was used to compare differences in phenomics characteristics among groups. Metabolites with variable importance in projection （VIP） ≥1 and fold change （FC） ≥2 were considered as representative differential metabolites. ResultsPatients in each TCM pattern type were most commonly in the 60-75 years age group， with a relatively high proportion of males. Regarding TCM symptoms， patients with phlegm-stasis binding syndrome most commonly presented with wiry-choppy or wiry-slippery pulse， chest pain， and chest tightness； in patients with qi deficiency-accompanied syndrome， fatigue， chest pain， and weak pulse were most common； in patients with qi stagnation-accompanied sydnrome， wiry-choppy or wiry-slippery pulse， chest pain， and symptoms that increase or decrease with emotional changes， belching， or flatulence were most common； in patients with toxin accumulation-accompanied syndrome， bitter taste in the mouth， chest tightness， irritability， restlessness or manic delirium， and dry and hard stools or foul-smelling diarrhea were most common. Comparisons among the different TCM patterns showed statistically significant differences in coagulation parameters （P＜0.05）， whereas no statistically significant difference was found in blood lipid levels and immune indicators （P＞0.05）.Metabolomics analysis suggested that disorders of glycerophosphate metabolism and valine， leucine， and isoleucine metabolism are characteristic features of phlegm-stasis binding syndrome and its accompanied patterns. The representative differential metabolites between phlegm-stasis binding syndrome and qi deficiency-accompanied syndrome were 7,8-dihydrobiopterin （FC = 3.58） and oxypurinol （FC = 124.50）. Those between phlegm-stasis binding syndrome and qi stagnation-accompanied syndrome were cytidine 5′-diphosphocholine （FC = 2.62） and D-mannosamine （FC = 2.99）. Those between phlegm-stasis binding syndrome and toxin accumulation-accompanied syndrome were anserine （FC = 7.83） and canrenone （FC = 8.94）. ConclusionThere are certain differences in the clinical characteristics and phenomics characteristics among patients with SAP due to CHD exhibiting phlegm-stasis binding syndrome and its accompanied patterns. The phenomics characteristics mainly involve biological alterations such as lipid metabolism disorders， amino acid metabolism abnormalities， and multiple immune indicators activation， which may provide a reference for precise differentiation and treatment of SAP of CHD in TCM clinical practice.

ObjectiveTo explore the clinical and phenomics characteristics of patients with phlegm-stasis binding syndrome and its accompanied patterns in stable angina pectoris （SAP） of coronary heart disease （CHD）. MethodsA multicenter cross-sectional study design was adopted. A total of 300 patients with SAP of CHD were enrolled and classified into 120 cases of phlegm-stasis binding syndrome， 125 cases of qi deficiency-accompanied syndrome， 38 cases of qi stagnation-accompanied syndrome， and 17 cases of toxin accumulation-accompanied syndrome according to traditional Chinese medicine （TCM） patterns. Data of patients with different TCM patterns were collected， including general condition， TCM symptoms， blood lipids， coagulation function， immune indicators， and serum metabolomics. Metabolic pathway enrichment analysis was used to compare differences in phenomics characteristics among groups. Metabolites with variable importance in projection （VIP） ≥1 and fold change （FC） ≥2 were considered as representative differential metabolites. ResultsPatients in each TCM pattern type were most commonly in the 60-75 years age group， with a relatively high proportion of males. Regarding TCM symptoms， patients with phlegm-stasis binding syndrome most commonly presented with wiry-choppy or wiry-slippery pulse， chest pain， and chest tightness； in patients with qi deficiency-accompanied syndrome， fatigue， chest pain， and weak pulse were most common； in patients with qi stagnation-accompanied sydnrome， wiry-choppy or wiry-slippery pulse， chest pain， and symptoms that increase or decrease with emotional changes， belching， or flatulence were most common； in patients with toxin accumulation-accompanied syndrome， bitter taste in the mouth， chest tightness， irritability， restlessness or manic delirium， and dry and hard stools or foul-smelling diarrhea were most common. Comparisons among the different TCM patterns showed statistically significant differences in coagulation parameters （P＜0.05）， whereas no statistically significant difference was found in blood lipid levels and immune indicators （P＞0.05）.Metabolomics analysis suggested that disorders of glycerophosphate metabolism and valine， leucine， and isoleucine metabolism are characteristic features of phlegm-stasis binding syndrome and its accompanied patterns. The representative differential metabolites between phlegm-stasis binding syndrome and qi deficiency-accompanied syndrome were 7,8-dihydrobiopterin （FC = 3.58） and oxypurinol （FC = 124.50）. Those between phlegm-stasis binding syndrome and qi stagnation-accompanied syndrome were cytidine 5′-diphosphocholine （FC = 2.62） and D-mannosamine （FC = 2.99）. Those between phlegm-stasis binding syndrome and toxin accumulation-accompanied syndrome were anserine （FC = 7.83） and canrenone （FC = 8.94）. ConclusionThere are certain differences in the clinical characteristics and phenomics characteristics among patients with SAP due to CHD exhibiting phlegm-stasis binding syndrome and its accompanied patterns. The phenomics characteristics mainly involve biological alterations such as lipid metabolism disorders， amino acid metabolism abnormalities， and multiple immune indicators activation， which may provide a reference for precise differentiation and treatment of SAP of CHD in TCM clinical practice.

ObjectiveTo establish the mouse models of psoriasis-like lesions induced by continuous cold exposure or intermittent cold exposure combined with imiquimod （IMQ）， and to evaluate the interventional effects of Kaixuan Jiedu core prescription （KXJD） on the two models. MethodsMale C57BL/6J mice were selected and classified into two experimental batches. The first batch of 36 mice was randomized into a room temperature group， a continuous cold exposure （10 ℃/24 h） group， and an intermittent cold exposure （10 ℃/6 h） group. Each group was further divided into a normal subgroup and a model subgroup （topical application of IMQ to induce skin lesions）， with 6 mice in each subgroup， for modeling and evaluation. The second batch of 54 mice， with 6 in each group， were subjected to the same temperature grouping with an additional KXJD （30.42 g·kg^-1， continuous gavage for 5 days） group. Comprehensive evaluation of model characteristics and KXJD efficacy was conducted through Psoriasis Area and Severity Index （PASI） scoring， skin temperature measurement by infrared thermography， histopathological observation by hematoxylin-eosin （HE） staining， detection of vascular endothelial growth factor （VEGF） and platelet endothelial cell adhesion molecule 1 （CD31） by immunohistochemistry， detection of Claudin-1 and Occludin by immunofluorescence assay， determination of serum levels of tumor necrosis factor-α （TNF-α） and interleukin （IL）-10 by enzyme-linked immunosorbent assay （ELISA）， and quantification of mRNA levels of IL-17A， IL-23， IL-6， and chemokine ligand 20 （CCL20） in skin lesions by quantitative Real-time polymerase chain reaction （Real-time PCR）. ResultsModel mice in all temperature groups exhibited typical psoriasis-like skin lesions. Compared with the normal groups， the model groups showed increased PASI scores， decreased skin temperatures （P<0.05）， obvious epidermal thickening， parakeratosis， and dermal inflammatory cell infiltration， as well as elevated mRNA levels of IL-17A， IL-23， IL-6， and CCL20 （P<0.05）. Cold exposure further aggravated psoriasis. The total PASI score of the intermittent cold exposure model group was higher than that of the room temperature model group （P<0.05）. The serum IL-10 did not show a compensatory elevation， and the blood vessels presented a characteristic of elevated CD31 expression （P<0.05） without a synchronous increase in VEGF. The continuous cold exposure model group exhibited more significant dermal capillary tortuosity and dilation， with the highest mRNA levels of IL-17A， IL-23， IL-6， and CCL20 among all groups. Compared with the respective model groups， KXJD intervention alleviated skin lesions， reduced epidermal thickness and inflammatory cell infiltration， and increased skin temperature， with the temperature increase being particularly significant in the intermittent cold exposure+KXJD group （P<0.05）. Furthermore， KXJD down-regulated the expression of VEGF and CD31， restored the expression of Claudin-1 and Occludin， decreased the mRNA levels of IL-17A and IL-23 （P<0.05）， and reduced the serum TNF-α level. ConclusionThis study successfully established compound psoriasis-like mouse models induced by cold exposure combined with IMQ. It confirms that cold aggravates the severity of psoriasis by exacerbating the closure of Xuanfu （sweat pores）， microcirculation disorders， and immune imbalance. Moreover， different cold exposure patterns have distinct mechanism differences. Continuous cold exposure focuses on enhancing the inflammatory response via the IL-23/IL-17 axis and angiogenesis， simulating chronic aggravation under a long-term cold environment. Intermittent cold exposure tends to impair immune regulation and induce microvascular endothelial stress， corresponding to acute exacerbations caused by sudden temperature drops. KXJD can effectively alleviate psoriasis-like skin lesions under cold conditions by unblocking Xuanfu， regulating vasomotor function， and correcting abnormal immune-inflammatory responses.

In recent years, remarkable advances have been made in lung transplantation regarding donor evaluation, perioperative support, and immune regulation, laying a foundation for improving the long-term survival of recipients. The rational application of extended criteria donor lungs and the popularization of ex vivo lung perfusion technology have enhanced the utilization efficiency and safety of marginal donor lungs. Refined assessment of risk factors such as donors’ smoking history and infectious status has reduced the incidence of primary graft dysfunction and postoperative infection. In terms of recipient management, individualized assessment systems have been gradually optimized. The standardized application of intraoperative extracorporeal membrane oxygenation and the development of minimally invasive surgical techniques have effectively alleviated perioperative injuries. The optimization of postoperative immunosuppressive regimens and advances in rejection monitoring technologies have further improved the long-term survival of grafts. Nevertheless, chronic graft dysfunction and recurrent or refractory infections remain the major bottlenecks restricting long-term prognosis. Future research should focus on expanding donor sources, improving the repair quality of marginal donor lungs, promoting precise immunosuppression and anti-infection strategies, and comprehensively enhancing the long-term survival rate and quality of life of lung transplant recipients.

The pathological changes of myocardial infarction （MI） are mainly characterized by progressive myocardial ischemic necrosis， decline in cardiac diastolic function， thinning of the ventricular wall， and enlargement of the ventricles. The clinical manifestations include myocardial ischemia， heart failure， arrhythmia， shock， and even sudden cardiac death， rendering MI one of the most perilous cardiovascular diseases. Currently， the clinical treatment for MI primarily involves interventional procedures and drug therapy. However， due to their significant side effects and high complication rates associated with these treatments， they fail to ensure a satisfactory quality of life and long-term prognosis for patients. On the other hand， traditional Chinese medicine has demonstrated remarkable potential in improving patient prognosis while reducing side effects. Research has elucidated that various signaling pathways such as nuclear transcription factor-κB （NF-κB）， adenosine 5̒-monophosphate-activated protein kinase （AMPK）， transforming growth factor-β （TGF-β）/Smads， mitogen-activated protein kinase （MAPK）， Wnt/β-catenin （β-catenin）， and phosphatidylinositol 3-kinase （PI3K）/protein kinase B（Akt） play crucial roles in regulating the occurrence and development of MI. Effectively modulating these signaling pathways through its therapeutic interventions， traditional Chinese medicine can enhance MI management by inhibiting apoptosis， providing anti-inflammatory properties， alleviating oxidative stress levels， and resisting myocardial ischemia. Due to its notable efficacy and favorable safety， it has become an area of focus in clinical practice.

BackgroundInvasive vagus nerve stimulation therapy has been approved for the adjunctive treatment of treatment-resistant depression， which may contribute to the anti-inflammatory properties of vagus nerve stimulation （VNS）， whereas the efficacy of non-invasive transcutaneous cervical vagus nerve stimulation （tcVNS） in treating major depressive disorder （MDD） and its impact on plasma inflammatory factors remain unclear. ObjectiveTo observe the effect of escitaloprom combined with tcVNS on the status of depression， anxiety and sleep quality as well as the plasma levels of interleukin-6 （IL-6） and interleukin-10 （IL-10） in MDD patients， in order to provide references for the recovery and treatment of MDD patients. MethodsFrom August 21， 2019 to April 17， 2024， 45 patients who met the diagnostic criteria for MDD in the Diagnostic and Statistical Manual of Mental Disorders， fifth edition （DSM-5） were recruited from the psychosomatic outpatient clinic of the First Affiliated Hospital of Air Force Military Medical University. Subjects were divided into study group （n=23） and control group （n=22） using random number table method. All patients were treated with escitalopram. On this basis， study group added a 30-minute tcVNS therapy once a day for 4 weeks. While control group was given corresponding sham stimulation， and the duration of each stimulation lasted 30 seconds. Before and after 4 weeks of treatment， Hamilton Depression Scale-17 item （HAMD-17） was used to assess depressive symptoms， and HAMD-17 anxiety/somatization subfactor and insomnia subfactor were used to assess patients' anxiety/somatization symptoms and sleep quality. Levels of plasma IL-6 and IL-10 were measured by enzyme-linked immunosorbent assay （ELISA）. ResultsThe generalized estimating equation model yielded a significant time effect for HAMD-17 total score， anxiety/somatization subfactor score and insomnia subfactor score in both groups （Wald χ²=315.226， 495.481， 82.420， P<0.01）. After 4 weeks of treatment， HAMD-17 total score and anxiety/somatization subfactor score of study group were lower than those of control group， with statistically significant differences （Wald χ²=4.967， 32.543， P<0.05 or 0.01）， while no statistically significant difference was found in the insomnia subfactor score between two groups （Wald χ²=0.819， P=0.366）. Significant time effects were reported on plasma IL-6 and IL-10 levels in both groups （Wald χ²=21.792， 5.242， P<0.05 or 0.01）. Compared with baseline data， a reduction in plasma IL-6 levels was detected in both groups （Wald χ²=22.015， 6.803， P<0.01）， and an increase in plasma IL-10 levels was reported in study group （Wald χ²=5.118， P=0.024） after 4 weeks of treatment. ConclusionEscitalopram combined with tcVNS therapy is effective in improving depressive symptoms， anxiety/somatization symptoms and sleep quality in patients with MDD. Additionally， it helps reduce plasma IL-6 levels and increase IL-10 levels. ［Funded by Shaanxi Provincial Key Research and Development Program-General Project （number， 2023-YBSF-185）， www.clinicaltrials.gov number， NCT04037111］