Background Meteorological factors are among the key extrinsic triggers for the onset and exacerbation of cardiovascular and cerebrovascular diseases (CVD). Against the backdrop of sustained global warming, elucidating the impact of ambient temperature and atmospheric pressure on CVD, especially on pre-hospital CVD emergent events, has become imperative for evidence-based prevention and emergency preparedness. Objective To quantify the temporal trends of daily mean temperature and atmospheric pressure and their associations with pre-hospital CVD emergent events in Yantai, and to explore effect modification by demographic subgroups and geographic areas, thereby providing an empirical basis for the rational allocation of emergency medical resources. Methods Pre-hospital CVD emergency data from January 1, 2016 to December 31, 2022 were selected from the Yantai 120 Emergency Medical Command System. Synchronous meteorological factors and environmental pollutant data were obtained from the websites of the National Oceanic and Atmospheric Administration and the National Centers for Environmental Information of the United States. Time-series analysis combined with distributed lag non-linear model was used to analyze the association between daily temperature, atmospheric pressure, and pre-hospital CVD emergencies. Average annual percentage changes (AAPC) were calculated using Joinpoint (version 5.2.0.0) to reflect temporal trends. Spearman correlation analysis was employed to screen variables with low collinearity for inclusion in the multi-pollutant adjusted models. Results From 2016 to 2022, a total of 268275 pre-hospital CVD emergency cases were recorded in Yantai City, showing an increasing annual trend (AAPC=7.16%, P=0.002). Emergency volume, daily temperature, atmospheric pressure, and PM₁₀ displayed marked seasonal patterns, with summer bottoms and winter peaks for emergencies, daily temperature higher in summer and lower in winter, and the inverse for atmospheric pressure. Correlation analysis showed that both daily temperature and atmospheric pressure were positively correlated with O₃, with correlation coefficients of 0.745 and 0.723, respectively; negatively correlated with PM_2.5, with correlation coefficients of −0.246 and −0.235, respectively; and negatively correlated with PM₁₀, with correlation coefficients of −0.241 and −0.229, respectively (P < 0.05). A U-shaped correlation was revealed between daily temperature and pre-hospital CVD emergency risk. Using 24 °C as the reference, risk was elevated at both low and high daily temperatures, with the strongest cold effect at –11 ℃ (RR=1.53, 95%CI: 1.35, 1.73) and strongest heat effect at 30 ℃ (RR=1.06, 95%CI: 1.01, 1.11); cold effects were prolonged (lag0–14d) compared to heat effects. Atmospheric pressure exhibited a positive association with emergencies; using a median of 1012.2 hPa as the reference, the maximum risk occurred at 1037 hPa (RR=1.25, 95%CI: 1.03, 1.52), with a 0–5 d lag. In the subgroup analysis, the impact of low temperature was more pronounced among females (RR_{−9 ℃}=1.66), people over 60 years (RR_{−9 ℃}=1.91), and rural residents (RR_{−9 ℃}=1.76). High atmospheric pressure significantly increased the risk of pre-hospital CVD emergency among females (RR_{1035 hPa}=1.54), people under 18 years (RR_{1035 hPa}=3.62), and rural residents (RR_{1035 hPa}=1.46). The sensitivity analysis results indicated that, after excluding the impact of pandemic and the degrees of freedom for time variations, the effects of daily average temperature and atmospheric pressure on pre-hospital CVD emergency volume were consistent with the primary analysis findings. Conclusion Both non-optimal temperature and atmospheric pressure significantly increase the pre-hospital CVD emergency volume; therefore, emergency medical services should develop targeted response protocols, optimize resource allocation, and intensify public health education to mitigate the pre-hospital burden attributable to temperature and pressure fluctuations.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.

ObjectiveThe accuracy of Y-chromosome haplogroup assignment is crucial for tracing paternal lineage in male samples. With the advancement of high-throughput sequencing technologies, high-density Y-SNP genotyping from whole-genome or array-based data has become a standard method for determiningY-chromosome haplogroups. This study systematically evaluated the performance of 4 commonly used high-density SNP genotyping systems—namely, the Global Screening Array (GSA), Chinese Genotyping Array (CGA), Affymetrix array, and the 1240K capture panel—for haplogroup assignment. This work provides a reference for data comparison across different systems. MethodsWe extracted genotype data for the 4 Y-SNP panels from 30× whole-genome sequencing (WGS) data of 1 590 male samples from the 1000 Genomes Project. Additionally, GSA array genotype data from 384 relative pairs (spanning 1st- to 12th-degree relationships) from 109 Chinese Han families were collected. Haplogroup assignment was performed using Y-LineageTracker v1.3.0 software. We assessed the concordance and resolution of haplogroup assignments between the four Y-SNP panels and the WGS data. The consistency and resolution of haplogroup assignments were also evaluated for both the 1000 Genomes Project samples and the 109 family samples collected in this study. Furthermore, the impact of varying numbers of Y-SNPs on haplogroup assignment was examined. ResultsThe GSA and CGA panels demonstrated superior resolution and discrimination of haplogroup subclades compared with the other two panels. The haplogroup assignments from the GSA, CGA, and 1240K panels showed high concordance with WGS data, with consistency rates exceeding 88.70%, whereas the Affymetrix platform exhibited a significantly lower consistency rate of 61.89%. Specifically, the GSA and CGA panels consistently demonstrated superior performance compared with the other two panels in the assignment of haplogroups O-M175 and H-L901, achieving complete concordance (100%) for both haplogroups. In contrast, the Affymetrix panel erroneously assigned all individuals belonging to haplogroup O-M175 to haplogroup K2-M526. Furthermore, its accuracy for haplogroup H-L901 was exceedingly low, at merely 1.41%. This poor performance was characterized by the misassignment of 98.59% of H-L901 samples—specifically, 1.41% to J-M304 and a predominant 97.18% to F-M89. For haplogroup R-M207, all four panels exhibited uniformly high levels of consistency, with concordance values exceeding 94.00%. Notably, for haplogroup E-M96, the 1240K and Affymetrix panels outperformed the GSA and CGA panels in terms of concordance, representing the first instance in which these two panels surpassed the latter. Conversely, for haplogroups J-M304, Q-M242, and I-M170, all 4 panels showed relatively elevated misclassification rates, with the Affymetrix array demonstrating the poorest overall performance. None of the four panels showed any discordant haplogroup assignments among the familial relative pairs analyzed. A positive correlation was observed between the number of Y-SNPs (ranging from 1 000 to 10 000) and classification consistency; however, classification consistency plateaued when the number of Y-SNPs exceeded 10 000. Furthermore, a random sampling analysis conducted on the GSA and CGA panels demonstrated that the haplogroup misclassification rate exhibited negligible fluctuation across the Y-SNP range of 500 to 1 000. Conversely, a marked enhancement in classification consistency was observed as the number of markers increased from 1 000 to 5 000, ultimately reaching a plateau within the interval of 5 000 to 8 000 markers. ConclusionThese findings indicate that the GSA and CGA panels provide high resolution and concordance, delivering reliable Y-haplogroup assignment for forensic investigations.

Metabolic associated fatty liver disease (MAFLD) is fundamentally driven by an imbalance in hepatic fatty-acid flux: the influx of fatty acids exceeds the liver’s capacity for disposal, resulting in excessive hepatic lipid accumulation, predominantly in the form of triglycerides (TGs). The occurrence and progression of MAFLD depend on disordered regulation across multiple metabolic steps, including fatty-acid uptake, de novo lipogenesis (DNL), fatty-acid oxidation (FAO), and very low-density lipoprotein (VLDL) export. Forkhead box protein O1 (FOXO1) is a key transcriptional regulator within the hepatic network coordinating glucose and lipid metabolism. Under metabolic stress and insulin resistance (IR), FOXO1 expression is frequently increased, whereas its inhibitory phosphorylation is reduced. These changes enhance FOXO1 nuclear localization and transcriptional activity, thereby reprogramming the expression of genes related to metabolism in the liver. Because hepatic lipid deposition is the central pathological feature of MAFLD, the functional status of FOXO1 directly influences hepatic lipid homeostasis. Growing evidence suggests that FOXO1 can exert bidirectional, environment-dependent effects on hepatic lipid accumulation; however, the molecular basis for this functional switch remains incompletely understood. This review systematically summarizes the biological functions and regulatory mechanisms of FOXO1 and its roles in hepatic lipid metabolism, with a particular focus on its crosstalk with insulin signaling. FOXO1 expression is shaped by RNA modifications and epigenetic regulation mediated by non-coding RNAs. Its transcriptional output is precisely governed by post-translational modifications—such as phosphorylation and acetylation—as well as by coordinated nucleocytoplasmic shuttling. Notably, these regulatory patterns vary markedly across nutritional states, degrees of insulin resistance, and stages of disease. In the fed state, insulin/IGF-1 signaling activates the PI3K-AKT pathway, promoting the inhibitory phosphorylation of FOXO1 and facilitating additional modifications, including acetylation, methylation, and ubiquitination. Together, these events drive FOXO1 export from the nucleus and dampen its transcriptional activity, suppressing gluconeogenesis and constraining lipogenic programs. Conversely, during fasting or when insulin signaling is weakened, FOXO1 inhibition is relieved. FOXO1 accumulates in the nucleus, binds to DNA, and regulates the transcription of downstream target genes. Mechanistically, FOXO1 can aggravate hepatic lipid accumulation by activating genes involved in TG synthesis while repressing FAO-related pathways, thereby favoring storage over oxidation. However, under specific conditions, FOXO1 may also alleviate the hepatic lipid burden by promoting TG hydrolysis and enhancing VLDL secretion, thereby reducing the net hepatic lipid load. In addition, lipotoxic signals mediated by ceramides and diacylglycerols (Cer/DAG) activate atypical protein kinase C (aPKC), further exacerbating the disruption of the AKT-FOXO1 axis. This vicious cycle ultimately produces a metabolic paradox in which increased hepatic glucose output coexists with persistent, insulin-independent lipogenesis, accelerating MAFLD progression. Importantly, FOXO1 regulation is not uniform: during early metabolic overload, insulin-mediated suppression may remain effective, whereas in advanced insulin resistance, the loss of AKT control permits sustained FOXO1 activity. Such stage-dependent dynamics may help explain why FOXO1 can either promote steatosis or, in certain contexts, support programs that facilitate lipid turnover. Accordingly, interventions should be liver-specific and tuned to the disease stage, aiming to curb maladaptive FOXO1 signaling while preserving its capacity to promote triglyceride hydrolysis and VLDL secretion when advantageous. Overall, this review offers an important perspective on MAFLD pathogenesis, emphasizing FOXO1 as a potential therapeutic target and providing a theoretical basis for developing liver-specific, disease-course-dependent precision interventions.

Objective To explore the effectiveness of the generative large language model MedGo in nursing decision-making for elderly patients with multimorbidity. Methods A quasi-randomized controlled trial study was conducted involving 6 junior nurses, 6 senior nurses and the MedGo model from January 1, 2025 to March 31, 2025 at the Emergency Internal Medicine Ward of Shanghai East Hospital Affiliated to Tongji University. Clinical data of 120 elderly patients with multimorbidity were analyzed to compare the performance of the three groups in four tasks (nursing diagnosis assessment, nursing intervention formulation, complication identification, and complication prevention) from three evaluation dimensions: decision-making time consumption, decision accuracy, and decision-making quality. Results In terms of decision-making time, the senior nurse group completed all four tasks faster than the junior nurse group (P＜0.01), and the MedGo group completed all four tasks faster than the junior nurse group (P＜0.001) and the senior nurse group (P＜0.001). In terms of decision-making accuracy, senior nurse group scored higher than junior nurse group in all four tasks (P＜0.001), while the MedGo group outperformed the senior nurse group only in complication identification (P＜0.001). In terms of decision-making quality, the MedGo group scored higher than junior nurse group (P＜0.001) and senior nurse group (P＜0.001) in all four tasks. Conclusions The MedGo model demonstrates advantages of high efficiency, accuracy, and quality in nursing decision-making for elderly patients with multimorbidity; senior nurses outperform junior nurses in decision-making, providing diverse references for clinical nursing decision-making.

OBJECTIVE To evaluate the anticoagulant efficacy and safety of nafamostat mesylate （NM） in the treatment of uremic patients at high risk of bleeding undergoing continuous veno-venous hemofiltration （CVVH） with different methods （pre- dilution and post-dilution）. METHODS A total of 130 uremic patients at high risk of bleeding who underwent CVVH treatment in the nephrology department of Chongqing University Three Gorges Hospital from July 2023 to September 2024 were selected. They were divided into pre-dilution group and post-dilution group according to the random number table method， with 65 cases in each group. Both groups of patients received CVVH treatment under NM anticoagulation. The pre-dilution group adopted the pre-dilution replacement method， while the post-dilution group adopted the post-dilution replacement method. The coagulation， pressure， and usage duration of the filter and dialysis circuit venous reservoirs were compared between the two groups. The changes in prothrombin time （PT）， prothrombin time-international normalized ratio （PT-INR）， activated partial thromboplastin time （APTT）， and fibrinogen （FIB） in the peripheral venous blood before the heparin pump and after the filter at 1， 4 and 7 h of CVVH treatment， as well as 20 min after the end of treatment， were compared between the two groups. The single-compartment urea clearance rate （spKt/V）， β2-microglobulin （β2-MG） clearance rate and the incidence of adverse reactions were duni2007@foxmail.com compared between the two groups. RESULTS Both the pre-dilution and post-dilution groups had 60 patients who completed the study. The incidence of grade Ⅱ-Ⅲ coagulation of the filter and venous reservoirs， as well as the number of patients with transmembrane and venous pressure alarm intervention in the post- dilution group were significantly higher or more than those in the pre-dilution group （P＜0.05）， while usage time of the filter and the pipeline in the post-dilution group was significantly shorter than that in the pre-dilution group （P＜0.05）. The APTT values before the heparin pump as well as PT and APTT values after the filter at 1 h， 4 h， and 7 h of CVVH treatment in the post-dilution group were significantly higher than those in the pre-dilution group （P＜0.001）. There were no significant differences in PT， PT- INR， APTT and FIB between the two groups of patients 20 min after the end of treatment （P＞0.05）. The spKt/v and β2-MG clearance rates in the post-dilution group were significantly higher than those in the pre-dilution group （P＜0.001）. There was no significant difference in the incidence of adverse reactions between the two groups （P＞0.05）. CONCLUSIONS When NM is used as an anticoagulant in the CVVH treatment of uremic patients at high risk of bleeding， compared with the pre-dilution treatment method， the post-dilution treatment method has a higher incidence of filter and dialysis tubing venous reservoir， a shorter usage time of the filter and pipeline， and a greater impact on extracorporeal coagulation， but has a higher solute clearance rate. Clinically， different dilution methods can be selected according to the different treatment needs of patients.

The theory of constitution in traditional Chinese medicine （TCM） has emerged as a new discipline in recent years. Constitution plays a vital role in the onset，progression，transformation，and prognosis of diseases. At present，some clinical scholars have adopted a novel diagnostic and treatment model of "constitution differentiation-disease identification-syndrome differentiation"，in which constitution is regarded as a core element throughout the diagnostic and therapeutic process. Constitution is closely associated with etiology，onset，pathogenesis，syndrome differentiation，and treatment. Against this background，the construction of animal models based on constitution holds far-reaching significance for advancing clinical research. This paper focuses on the construction and evaluation of rodent models with Qi-deficiency constitution，aiming to explore how to further induce Qi-deficiency syndromes and related disease states on the basis of Qi-deficiency constitution models，thereby developing an integrated animal model that embodies the trinity of "constitution-disease-syndrome". The establishment of this model not only provides a solid experimental foundation for the development of new therapies and drugs in TCM targeting specific constitutions，diseases，and syndromes，but also greatly promotes the modernization and scientific advancement of TCM theory. By comprehensively applying multidisciplinary technologies and methods，the study evaluates the model's validity，reliability，and practicality，with the aim of opening new avenues for future research in TCM and promoting the development of the field.

The theory of constitution in traditional Chinese medicine （TCM） has emerged as a new discipline in recent years. Constitution plays a vital role in the onset，progression，transformation，and prognosis of diseases. At present，some clinical scholars have adopted a novel diagnostic and treatment model of "constitution differentiation-disease identification-syndrome differentiation"，in which constitution is regarded as a core element throughout the diagnostic and therapeutic process. Constitution is closely associated with etiology，onset，pathogenesis，syndrome differentiation，and treatment. Against this background，the construction of animal models based on constitution holds far-reaching significance for advancing clinical research. This paper focuses on the construction and evaluation of rodent models with Qi-deficiency constitution，aiming to explore how to further induce Qi-deficiency syndromes and related disease states on the basis of Qi-deficiency constitution models，thereby developing an integrated animal model that embodies the trinity of "constitution-disease-syndrome". The establishment of this model not only provides a solid experimental foundation for the development of new therapies and drugs in TCM targeting specific constitutions，diseases，and syndromes，but also greatly promotes the modernization and scientific advancement of TCM theory. By comprehensively applying multidisciplinary technologies and methods，the study evaluates the model's validity，reliability，and practicality，with the aim of opening new avenues for future research in TCM and promoting the development of the field.

Objective: To investigate the status of antibiotic exposure in third grade primary school students in Shenzhen,so as to provide evidence for the scientific management of antibiotic use and reduction of population health risks. Methods: From 1 September to 30 October 2021, 200 third grade students from 8 primary schools in Luohu District of Shenzhen were selected by cluster random sampling as research subjects. The body composition was measured, urine samples were collected, and the contents of 35 antibiotics in the samples were detected by mass spectrometry. Relevant dietary habit information of the subjects was collected via questionnaires. The Chi square test was used to compare the detection rate of antibiotics among different genders and weight grades. The Logistic regression model was adopted to evaluate the correlation between the target antibiotic detection rate and dietary habits. Results: At least one type of antibiotic was detected in 198 of the subjects with an overall detection rate of 99.0% . Among the 35 target antibiotics, 23 were detected with detection rates ranging from 0.5%-69.5%. Quinolones had the highest detection rate of 86.5% , followed by macrolides and sulfonamides with detection rates of 77.5% and 76.5%, respectively. The detection rate of antibiotics was 98.3% in boys and 100.0% in girls with no statistically significant difference ( χ 2=1.35, P >0.05). The detection rates of quinolones, macrolides, and sulfonamides varied significantly among children with different BMI categories ( χ 2=38.18, 12.45, 9.76 , all P <0.05). The multivariate Logistic regression model analysis showed that the macrolide detection rate was affected by genders( OR =0.42) and the sulfonamide detection rate was significantly correlated with the frequency of dairy product consumption and being overweight( OR =2.01)(both P <0.05). Enrofloxacin was associated with the weekly consumption frequency of livestock meat such as pork, beef and mutton, as well as the weekly consumption frequency of poultry meat such as chicken, duck and goose ( OR = 2.81,2.17,both P <0.05). Trimethoprim was associated with the weekly frequency of drinking pure milk ( OR =5.49, P < 0.05 ). Conclusions Third grade primary school students in Shenzhen are generally exposed to low dose antibiotics. Macrolides, quinolones, and sulfonamides may be associated with the risk of obesity in primary school students.

ObjectiveTo investigate the epidemiological characteristics and influencing factors of depressive symptoms among middle school (junior, senior, and vocational high school) students in Yunnan Province, China, and to inform evidence-based intervention strategies for adolescent mental health. MethodsA cross-sectional survey was conducted between October and November 2023 using stratified random cluster sampling. Students from eight counties (districts) across four prefectures (cities) in Yunnan Province were included. Depressive symptoms were assessed using the Center for Epidemiologic Studies Depression Scale (CES-D). Multivariate logistic regression analysis was performed to examine factors associated with depressive symptoms, with stratified analyses conducted by gender, educational stage, and ethnicity. ResultsA total of 4 672 questionnaires were distributed, with 4 670 valid questionnaires retrieved, yielding a valid response rate of 99.96%. The surveyed participants were predominantly female students (50.81%), junior high school students (49.13%), ethnic minorities (52.78%), and urban residents (79.29%). The mean CES-D score for middle school students in Yunnan Province was (15.31±10.83). Female students had a significantly higher mean score (16.63±11.41) than male students (13.95±10.02) (P<0.001). Senior high school students had a significantly higher mean score (16.61±10.61) compared to both junior high school students (14.74±11.45) and vocational high school students (13.10±7.71) (all pairwise comparisons P<0.001). The prevalence of depressive symptoms among middle school students in Yunnan Province was 28.18%. The prevalence was significantly higher in females (34.09%) than in males (22.07%). By school type, the detection rate was highest among senior high school students (33.39%), followed by junior high school students (26.29%) and vocational high school students (17.27%) (P<0.05). Multivariate logistic regression analysis showed that female gender (OR=2.16, 95%CI: 1.86‒2.50), being in junior high school (OR=2.43, 95%CI: 1.84‒3.20) or senior high school (OR=2.27, 95%CI: 1.73‒2.98), not living with parents (OR=1.24, 95%CI: 1.07‒1.44), irregular breakfast consumption (OR=1.52, 95%CI: 1.33‒1.75), lack of moderate-to-vigorous physical activity (MVPA) (OR=1.69, 95%CI: 1.37‒2.09), sleep duration ≤5 h per night (OR=2.52, 95%CI: 2.02‒3.14) or 6‒7 h per night (OR=1.47, 95%CI: 1.25‒1.73), smoking (OR=1.86, 95%CI: 1.56‒2.23), and alcohol consumption (OR=1.81, 95%CI: 1.54‒2.13) were positively associated with depressive symptoms. In contrast, screen time ≤1 h (OR=0.71, 95%CI: 0.59‒0.86) was negatively associated with depressive symptoms. Stratified analyses showed that female students not living with parents (OR=1.29, 95%CI: 1.06‒1.58), senior high school students (OR=1.51, 95%CI: 1.21‒1.88), and Han Chinese students (OR=1.37, 95%CI: 1.11‒1.69) were more likely to experience depressive symptoms. Han Chinese students who smoked were also more likely to have depressive symptoms (OR=1.72, 95%CI: 1.34‒2.21). In contrast, male students with screen time ≤1 h (OR=0.71, 95%CI: 0.53‒0.95) and ethnic minority students (OR=0.74, 95%CI: 0.58‒0.95) were less likely to experience depressive symptoms. Regardless of gender, irregular breakfast consumption, lack of MVPA, sleep duration less than 8 h per night, smoking, and alcohol consumption were all positively associated with depressive symptoms (P<0.05). Among both junior and senior high school students, irregular breakfast consumption, lack of MVPA, smoking, and alcohol consumption were positively associated with depressive symptoms (all P<0.05), while screen time ≤1 h was negatively associated with depressive symptoms (all P<0.05). For junior high school students, engaging in 1‒2 days of MVPA per week, screen time more than 2 h per day, and sleep duration 6‒7 h per night were all positively associated with depressive symptoms (all P<0.05). Among junior high, senior high, and vocational high school students, sleep duration ≤5 h per night was positively associated with depressive symptoms (P<0.05). For both Han Chinese and ethnic minority students, irregular breakfast consumption, lack of MVPA, sleep duration less than 8 h per night, and alcohol consumption were all positively associated with depressive symptoms (all P<0.05). ConclusionThe prevalence of depressive symptoms among middle school students in Yunnan Province is comparable to that in central China and higher than that in northern regions. Prevention and control efforts should prioritize female students and those in junior and senior high school stages. Universal improvements in lifestyle behaviors among middle school students, such as regular breakfast consumption, MVPA, sufficient sleep (≥8 h), and abstinence from smoking and alcohol. Particular attention should be given to limiting excessive screen time among junior high school students and addressing the mental health needs of females not living with their parents, senior high school students, and Han Chinese students.