Objective: To investigate the epidemiological characteristics and spatial-temporal clustering of varicella in Changchun City from 2020 to 2024, so as to provide the evidence for formulating local varicella prevention and control measures. Methods: The individual case data of varicella were collected through the Surveillance and Reporting Management System of the Chinese Disease Prevention and Control Information System in Changchun City from 2020 to 2024. Descriptive epidemiological methods were used to analyze the population ,regional, and temporal distribution. Spatial autocorrelation and spatio-temporal scanning analyses were used to identify the spatial-temporal clustering characteristics. Results: A total of 8 850 varicella cases were reported in Changchun City from 2020 to 2024, with an average annual incidence of 19.64/105. There were 4 929 male cases and 3 921 female cases, with a male-to-female ratio of 1.26∶1. The age was mainly 0-<20 years (6 649 cases, 75.13%), and students were the predominant occupation (6 036 cases, 68.20%). The top three counties (cities, districts) with the highest number of cases were Chaoyang District (1 944 cases), Gongzhuling City (1 054 cases) and Nanguan District (987 cases), accounting for 45.03%. The peak incidence periods were from April to June and from October to December, with 2 166 and 4 226 cases, accounting for 24.47% and 47.75%, respectively. Spatial autocorrelation analysis showed that spatial clustering existed from 2020 to 2024. The high-high clustering areas were mainly some townships (streets) in Chaoyang District, Nanguan District, Changchun New District and Jingyue District. Spatio-temporal scanning analysis identified 6 high-risk clustering areas. The class Ⅰ clustering area was Nanhu Street in Chaoyang District, with the clustering period from September 2020 to February 2022. Conclusions Varicella cases in Changchun City were mainly males and students aged 0-<20 years from 2020 to 2024. The peak incidence was mainly in winter. Chaoyang District was a high-risk area, with obvious spatial-temporal clustering.

ObjectiveHeadspace solid-phase microextraction-gas chromatography-mass spectrometry（HS-SPME-GC-MS） and GC-triple quadrupole MS（GC-QqQ-MS） in combination with non-targeted and targeted metabolomics were employed to systematically analyze the chemical composition differences of Xihuangwan prepared with natural musk and artificial musk， and establish an identification system for them. MethodsThe volatile components of 9 batches of Xihuangwan samples from 8 manufacturers were analyzed by HS-SPME-GC-MS non-targeted metabolomics， and identified by comparing their MS data with the National Institute of Standards and Technology（NIST） spectral library. Orthogonal partial least squares-discriminant analysis（OPLS-DA） was used to identify differential volatile components of Xihuangwan prepared with natural musk and artificial musk. Additionally， GC-QqQ-MS targeted metabolomics was applied to quantify the levels of α-pinene， β-elemene， muscone， dehydroepiandrosterone， bornyl acetate， and octyl acetate in 27 batches of samples from 9 manufacturers. Cluster analysis， principal component analysis（PCA）， and partial least squares-discriminant analysis（PLS-DA） were conducted to further explore the differences in volatile components between Xihuangwan samples prepared with natural musk and artificial musk. ResultsNon-targeted metabolomics identified 291 volatile compounds in Xihuangwan， including alkanes， esters， alkanes， alcohols， ketones， naphthalenes and others. OPLS-DA analysis revealed distinct separation between Xihuangwan samples containing artificial musk（A1， C1， D1， E1， F1， G1， I1） and those containing natural musk（H1， H3）. A total of 30 differential metabolites were identified. The relative contents of these 30 differential metabolites were visualized using a radar chart， revealing significant differences in the levels of octanol， borneol acetate and muscone. Cluster analysis and PCA results from targeted metabolomics indicated that Xihuangwan could be classified into two distinct groups：one composed of natural musk（H1， H3） and the other of artificial musk， sample H2. PLS-DA identified muscone， octyl acetate， and dehydroepiandrosterone as key differential volatile components. Although no significant difference was observed in the content of octyl acetate between the two groups， statistically significant differences were found for muscone and dehydroepiandrosterone（P<0.05）. ConclusionMuscone and dehydroepiandrosterone can be used for the differentiation of Xihuangwan samples containing natural musk from those containing artificial musk. This study systematically and comprehensively analyzed the differences in the types and contents of major volatile components in Xihuangwan prepared with natural musk and artificial musk， providing a scientific basis for quality evaluation and control of Xihuangwan.

ObjectiveHeadspace solid-phase microextraction-gas chromatography-mass spectrometry（HS-SPME-GC-MS） and GC-triple quadrupole MS（GC-QqQ-MS） in combination with non-targeted and targeted metabolomics were employed to systematically analyze the chemical composition differences of Xihuangwan prepared with natural musk and artificial musk， and establish an identification system for them. MethodsThe volatile components of 9 batches of Xihuangwan samples from 8 manufacturers were analyzed by HS-SPME-GC-MS non-targeted metabolomics， and identified by comparing their MS data with the National Institute of Standards and Technology（NIST） spectral library. Orthogonal partial least squares-discriminant analysis（OPLS-DA） was used to identify differential volatile components of Xihuangwan prepared with natural musk and artificial musk. Additionally， GC-QqQ-MS targeted metabolomics was applied to quantify the levels of α-pinene， β-elemene， muscone， dehydroepiandrosterone， bornyl acetate， and octyl acetate in 27 batches of samples from 9 manufacturers. Cluster analysis， principal component analysis（PCA）， and partial least squares-discriminant analysis（PLS-DA） were conducted to further explore the differences in volatile components between Xihuangwan samples prepared with natural musk and artificial musk. ResultsNon-targeted metabolomics identified 291 volatile compounds in Xihuangwan， including alkanes， esters， alkanes， alcohols， ketones， naphthalenes and others. OPLS-DA analysis revealed distinct separation between Xihuangwan samples containing artificial musk（A1， C1， D1， E1， F1， G1， I1） and those containing natural musk（H1， H3）. A total of 30 differential metabolites were identified. The relative contents of these 30 differential metabolites were visualized using a radar chart， revealing significant differences in the levels of octanol， borneol acetate and muscone. Cluster analysis and PCA results from targeted metabolomics indicated that Xihuangwan could be classified into two distinct groups：one composed of natural musk（H1， H3） and the other of artificial musk， sample H2. PLS-DA identified muscone， octyl acetate， and dehydroepiandrosterone as key differential volatile components. Although no significant difference was observed in the content of octyl acetate between the two groups， statistically significant differences were found for muscone and dehydroepiandrosterone（P<0.05）. ConclusionMuscone and dehydroepiandrosterone can be used for the differentiation of Xihuangwan samples containing natural musk from those containing artificial musk. This study systematically and comprehensively analyzed the differences in the types and contents of major volatile components in Xihuangwan prepared with natural musk and artificial musk， providing a scientific basis for quality evaluation and control of Xihuangwan.

Objective: To characterize longitudinal trajectories of testicular development in boys and breast development in girls, so as to provide reference data for understanding patterns of pubertal sexual maturation. Methods: Based on the Shanghai Pudong New Area Cohort Study on Growth, Development and Health in Children and Adolescents, a baseline survey was conducted in 2020 using a mult stage cluster random sampling method. A total of 2 184 children who completed all follow ups during the primary school period from 13 elementary schools in Pudong New Area,Shanghai,with annual follow ups during 2021-2025. Testicular volume and Tanner stage of breast development were assessed by professional physicians using standardized visual inspection and palpation. The age distribution of testicular volume and breast development was fitted by using cumulative link mixed models and Turnbull s nonparametric maximum likelihood estimation method. Results: Median ages for testicular volumes of 2, 3, 4 and 5 mL in boys were 7.07, 9.24, 10.29, and 11.57 years old, respectively. Median ages for Tanner breast stages Ⅱ, Ⅲ, Ⅳ, and Ⅴ in girls were 8.55 , 10.17, 11.18, and 13.78 years old, respectively. Based on overweight and obesity, stratified analysis showed that earlier pubertal onset among overweight/obesity children, and the key milestones for pubertal initiation were testicular volume reaching 4 mL in boys and breast Tanner II in girls for 10.29, 10.83; 8.18, 9.00 years. Conclusion Overweight and obesity are associated with earlier pubertal initiation,but there are certain gender and developmental stage specific patterns.

Work-related musculoskeletal disorders (WMSDs) represent a significant occupational health challenge among healthcare professionals globally, posing substantial threats to physical and mental well-being as well as work sustainability. Adopting preventive behaviors—including ergonomic postural adjustments, optimized work-rest scheduling, proper use of protective and assistive equipment, and regular physical activity—is essential for mitigating the risk of WMSDs. Guided by the social ecological model, the review synthesized current evidence on the determinants of WMSDs preventive behaviors across four levels: intrapersonal characteristics, work environment conditions, interpersonal support, and policy/institutional factors. The findings suggest that higher educational attainment, favorable health-related behavioral patterns, optimized ergonomic work environments, adoption of supportive collaborative systems, strong organizational support, as well as policy safeguards facilitate preventive behavior adoption. Conversely, limited prevention-related knowledge, low risk perception, insufficient physical activity, excessive workload, lack of appropriate protective equipment, inadequate ergonomic training, a prevailing culture of presenteeism, and inadequate policy implementation constitute significant barriers. Multi-dimensional intervention strategies targeting these determinants are warranted to enhance preventive behaviors, reduce the risk of WMSDs, and strengthen occupational health protection for healthcare professionals.

Exercise-induced metabolic remodeling is a fundamental adaptive process whereby the body reorganizes systemic and cellular metabolism to meet the dynamic energy demands posed by physical activity. Emerging evidence reveals that such remodeling not only enhances energy homeostasis but also profoundly influences immune function through complex molecular interactions involving glucose, lipid, and protein metabolism. This review presents an in-depth synthesis of recent advances, elucidating how exercise modulates immune regulation via metabolic reprogramming, highlighting key molecular mechanisms, immune-metabolic signaling axes, and the authors’ academic perspective on the integrated “exercise-metabolism-immunity” network. In the domain of glucose metabolism, regular exercise improves insulin sensitivity and reduces hyperglycemia, thereby attenuating glucose toxicity-induced immune dysfunction. It suppresses the formation of advanced glycation end-products (AGEs) and interrupts the AGEs-RAGE-inflammation positive feedback loop in innate and adaptive immune cells. Importantly, exercise-induced lactate, traditionally viewed as a metabolic byproduct, is now recognized as an active immunomodulatory molecule. At high concentrations, lactate can suppress immune function through pH-mediated effects and GPR81 receptor activation. At physiological levels, it supports regulatory T cell survival, promotes macrophage M2 polarization, and modulates gene expression via histone lactylation. Additionally, key metabolic regulators such as AMPK and mTOR coordinate immune cell energy balance and phenotype; exercise activates the AMPK-mTOR axis to favor anti-inflammatory immune cell profiles. Simultaneously, hypoxia-inducible factor-1α (HIF-1α) is transiently activated during exercise, driving glycolytic reprogramming in T cells and macrophages, and shaping the immune landscape. In lipid metabolism, exercise alleviates adipose tissue inflammation by reducing fat mass and reshaping the immune microenvironment. It promotes the polarization of adipose tissue macrophages from a pro-inflammatory M1 phenotype to an anti-inflammatory M2 phenotype. Moreover, exercise alters the secretion profile of adipokines—raising adiponectin levels while reducing leptin and resistin—thereby influencing systemic immune balance. At the circulatory level, exercise improves lipid profiles by lowering pro-inflammatory free fatty acids (particularly saturated fatty acids) and triglycerides, while enhancing high-density lipoprotein (HDL) function, which has immunoregulatory properties such as endotoxin neutralization and macrophage cholesterol efflux. Regarding protein metabolism, exercise triggers the expression of heat shock proteins (HSPs) that act as intracellular chaperones and extracellular immune signals. Exercise also promotes the secretion of myokines (e.g., IL-6, IL-15, irisin, FGF21) from skeletal muscle, which modulate immune responses, facilitate T cell and macrophage function, and support immunological memory. Furthermore, exercise reshapes amino acid metabolism, particularly of glutamine, arginine, and branched-chain amino acids (BCAAs), thereby influencing immune cell proliferation, biosynthesis, and signaling. Leucine-mTORC1 signaling plays a key role in T cell fate, while arginine metabolism governs macrophage polarization and T cell activation. In summary, this review underscores the complex, bidirectional relationship between exercise and immune function, orchestrated through metabolic remodeling. Future research should focus on causative links among specific metabolites, signaling pathways, and immune phenotypes, as well as explore the epigenetic consequences of exercise-induced metabolic shifts. This integrated perspective advances understanding of exercise as a non-pharmacological intervention for immune regulation and offers theoretical foundations for individualized exercise prescriptions in health and disease contexts.

ObjectiveTo study the regulatory effect of the partial sequence within the 3’ untranslated region （3’UTR） of slit guidance ligand 1 （Slit1） （Slit1-3’UTR） on the fibrotic phenotypes of cardiac fibroblasts （CFs） and its potential mechanism. MethodsThe adenovirus vector was used to overexpress the 1526nt sequence of Slit1-3’UTR in ICR neonatal mouse CFs （mCFs）. The expression of fibrosis-related genes in mCFs， such as collagen type 1 alpha1（COL1A1）， collagen type 3 alpha3 （COL3A1） and alpha smooth muscle actin （α-SMA） were detected by Western blot assay. The effect of Slit1-3’UTR 1526nt on the proliferation and migration of mCFs was assessed by EdU staining and Trans-well assays. Angiotensin Ⅱ （Ang Ⅱ） was used to treat mCFs， and the impact of Slit1-3’UTR 1526nt on the fibrotic phenotypes of Ang Ⅱ-induced mCFs was evaluated. After overexpression of Slit1-3’UTR 1526nt， miR-34a-5p mimic was transfected into mCFs， followed by actinomycin D treatment to detect the mRNA stability of Slit1-3’UTR 1526nt， and the levels of miR-34a-5p and its target gene SIRT1（si-SIRT1） in mCFs were determined. The effects of miR-34a-5p and small interfering RNA targeting SIRT1 on the Slit1-3’UTR 1526nt-mediated regulation of fibrotic phenotypes were also determined. ResultsAdenovirus-mediated overexpression of Slit 1-3’UTR 1526nt was achieved in mCFs. Overexpression of Slit 1-3’UTR 1526nt markedly inhibited the expression of the fibrosis-related genes， proliferation and migration of mCFs and fibrotic phenotypes of Ang Ⅱ. The results of actinomycin D assay showed that miR-34a-5p inhibited the stability of Slit1-3’UTR 1526nt in mCFs， while the level of miR-34a-5p was reduced in mCFs with overexpression of Slit1-3’UTR 1526nt. Transfection of miR-34a-5p promoted the fibrotic phenotypes， and reversed the inhibitory effect of Slit1-3’UTR 1526nt on the fibrotic phenotypes of mCFs. Overexpression of Slit1-3’UTR 1526nt significantly increased the level of miR-34a-5p target gene SIRT1 in mCFs. Transfection of miR-34a-5p and si-SIRT1 consistently reversed the inhibitory effects of Slit1-3’UTR 1526nt on the fibrotic phenotypes of mCFs. ConclusionSlit1-3’UTR1526nt inhibits the fibrotic phenotypes of mCFs by binding to miR-34a-5p and increasing the expression of its target gene of SIRT1.

Objective To clone and express the heat shock cognate protein 20 (SjHsc20) of Schistosoma japonicum, and to preliminarily investigate its biological characteristics. Methods The target fragment of the SjHsc20 gene was amplified using PCR assay and cloned into the pET-28a(+) expression plasmid to generate the recombinant expression vector pET-28a(+)-SjH-sc20, which was then transformed into Escherichia coli BL21 (DE3) competent cells. The recombinant SjHsc20 (rSjHsc20) protein was induced with isopropyl β-D-thiogalactopyranoside (IPTG) and purified, and the expression of the rSjHsc20 protein was checked with sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). The immunogenicity of the rSjHsc20 protein was detected using Western blotting, and the transcriptional levels of SjHsc20 were quantified in S. japonicum worms at different developmental stages and in male and female adult worms using real-time quantitative PCR (RT-qPCR) assay. Thirty female BALB/c mice at ages 6 to 8 weeks were divided into three groups, including the rSjHsc20 immunization group, the PBS control group, and the ISA 206 adjuvant group, of 10 mice in each group. Mice in the rSjHsc20 immunization group were subcutaneously immunized with 20 μg rSjHsc20 on days 1, 15 and 31, and animals in the PBS control group were subcutaneously injected with the same volume of PBS on days 1, 15 and 31, while mice in the ISA 206 adjuvant group were subcutaneously immunized with the same volume of ISA 206 adjuvant on days 1, 15 and 31, respectively. All mice in each group were infected with (40 ± 2) S. japonicum cercariae via the abdomen 14 day following the last immunization. Levels of serum specific IgG and its subtypes IgG1 and IgG2 antibodies against rSjHsc20, and the serum titers of anti-rSjHsc20 antibody were detected in mice using indirect enzyme-linked immunosorbent assay (ELISA). All mice were sacrifice 42 days post-infection, and S. japonicum worms were collected from the hepatic portal vein and counted. The eggs per gram (EPG), worm burden reductions and egg burden reductions were estimated to evaluate the protective efficacy of the rSjHsc20 protein. Results The SjHsc20 gene had an open reading frame (ORF) with 756 bp in length and encoded 252 amino acids, and the rSjHsc20 protein had a relative molecular mass of approximately 29 kDa. The rSjHsc20 protein was recognized by the serum of mice infected with S. japonicum and the serum of mice immunized with the rSjHsc20 protein, indicating that rSjHsc20 had a good immunogenicity. There was a significant difference in the transcriptional levels of the SjHsc20 gene among the 7-day (1.001 4 ± 0.065 7), 12-day (2.268 3 ± 0.129 2), 21-day (1.378 5 ± 0.160 4), 28-day (1.196 4 ± 0.244 0), 35-day (1.646 3 ± 0.226 1), 42-day worms of S. japonicum (1.758 0 ± 0.611 1) (F = 38.45, P < 0.000 1), and the transcriptional level of the SjHsc20 gene was higher in the 12-day worms than in worms at other developmental stages (all P values < 0.000 1). The serum levels of anti-rSjHsc20 IgG antibody were 0.106 6 ± 0.010 7, 0.108 3 ± 0.010 4, and 0.553 2 ± 0.069 1 in the PBS control group, ISA 206 adjuvant group, and rSjHsc20 immunization group following the last immunization, respectively, and the serum levels of IgG1 antibody were 0.137 3 ± 0.054 0, 0.181 1 ± 0.096 8, and 1.765 8 ± 0.221 1, while the levels of IgG2a antibody were 0.280 3 ± 0.197 6, 0.274 0 ± 0.146 3, and 1.560 4 ± 0.106 0, respectively. There were significant differences in the serum levels of anti-rSjHsc20 IgG (F = 397.70, P < 0.000 1), IgG1 (F = 401.00, P < 0.000 1) and IgG2a antibodies (F = 229.70, P < 0.000 1) among the three groups, and the serum levels of anti-rSjHsc20 IgG, IgG1 and IgG2a antibodies were higher in the rSjHsc20 immunization group than in the PBS control group and the ISA 206 adjuvant group (all P values < 0.000 1). There was a significant difference in the IgG1/IgG2a ratio among the rSjHsc20 immunization group (1.177 2 ± 0.143 6), the PBS control group (0.428 4 ± 0.199 8) and the ISA 206 adjuvant group (0.559 9 ± 0.181 1) (F = 43.97, P < 0.000 1), and the IgG1/IgG2a ratio was > 1 in the rSjHsc20 immunization group, which was higher than in the PBS control group and the ISA 206 adjuvant group (both P values < 0.000 1). The titers of serum anti-rSjHsc20 antibody were all above 1∶16 384 in the rSjHsc20 immunization group following immunizations on days 1, 15 and 31, indicating that the rSjHsc20 protein had a strong immunogenicity. The mean worm burdens were (16.60±5.75), (15.80±5.58) worms per mouse and (14.40±5.75) worms per mouse in the PBS control group, the ISA 206 adjuvant group and the rSjHsc20 immunization group 42 days post-infection with S. japonicum cercariae (F = 0.50, P > 0.05), and the EPG were 68 370 ± 22 690, 67 972 ± 19 502, and 41 075 ± 13 251 in the PBS control group, the ISA 206 adjuvant group and the rSjHsc20 immunization group (F = 4.55, P < 0.05), with lower EPG in the PBS control group and the ISA 206 adjuvant group than in the rSjHsc20 immunization group (both P values < 0.05). Immunization with the rSjHsc20 protein resulted in a worm burden reduction of 13.25% and an egg burden reduction of 39.92% relative to the PBS control group. Conclusions SjHsc20 is successfully cloned and expressed, and the rSjHsc20 protein induces partial immunoprotective effects in mice, which provides a basis for deciphering the biological functions of SjHsc20 and assessing the potential of SjH-sc20 as a vaccine candidate.

OBJECTIVE To optimize the salt-processing technology for Anemarrhena asphodeloides. METHODS Taking soaking time， stir-frying temperature， and stir-frying time as factors， Box-Behnken response surface methodology was employed to optimize the salt-processing technology of A. asphodeloides using the contents of mangiferin， neomangiferin， isomangiferin， timosaponin BⅡ， timosaponin AⅢ， timosaponin BⅢ， total flavonoids， and total saponins as evaluation indicators. The entropy weight method was applied to determine the weight of each indicator and calculate the comprehensive score. Based on the 17 sets of Box-Behnken response surface methodology results， a genetic algorithm （GA）-back propagation （BP） neural network was used to further optimize the salt-processing technology， with soaking time， stir-frying temperature， and stir-frying time as input layers and the comprehensive score as the output layer. The salt-processing parameters obtained from the two methods were validated and compared to determine the optimal salt-processing technology for A. asphodeloides. RESULTS The optimal salt-processing conditions obtained via the Box-Behnken response surface methodology were as follows： soaking time of 23 min， stir-frying temperature of 160 ℃ ， and stir-frying time of 12 min， yielding a comprehensive score of 63.370 2. The GA-BP neural network optimization resulted in the following conditions： soaking time of 24 min， stir-frying temperature of 163 ℃， and stir-frying time of 12 min， with a comprehensive score of 65.163 8. The GA-BP neural network optimization outperformed the results obtained by Box-Behnken response surface methodology. CONCLUSIONS This study successfully optimized the salt-processing technology for A. asphodeloides. Specifically， the technology involves adding 15 mL of 0.1 g/mL saline solution to 50 g of the herbal slices， allowing them to moisten for 24 minutes， and then stir-frying at 163 ℃ for 12 minutes.

Objective: To investigate the detection of depressive symptoms and its influencing factors among middle school students in Huzhou City, so as to provide insights for improving the mental health levels among middle school students. Methods: From September to November 2024, a total of 4 729 middle school students from five counties (districts) in Huzhou City were selected through the stratified cluster random sampling method. Demographic information, lifestyle, and school bullying were collected through questionnaire surveys. Depressive symptoms were assessed using the Center for Epidemiological Studies Depression Scale (CES-D). Factors affecting depressive symptoms among middle school students were analyzed using a multivariable logistic regression model. Results: A total of 4 729 middle school students were surveyed, including 2 200 boys (46.52%) and 2 529 girls (53.48%). Depressive symptoms were detected in 1 026 students, with a detection rate of 21.70%. Multivariable logistic regression analysis showed that girl (OR=1.960, 95%CI: 1.659-2.317), high school (ordinary high school, OR=1.789, 95%CI: 1.465-2.186; vocational high school, OR=1.581, 95%CI: 1.105-2.263), consumption of sugar-sweetened beverages >0 time/day (<1 time/day, OR=1.363, 95%CI: 1.009-1.841; ≥1 time/day, OR=1.568, 95%CI: 1.098-2.239), fried food intake ≥1 time/day (OR=1.890, 95%CI: 1.291-2.769), skipping breakfast daily (OR=2.178, 95%CI: 1.825-2.599), TV viewing time ≥2 hours/day (OR=1.457, 95%CI: 1.154-1.838), insufficient sleep duration (OR=1.761, 95%CI: 1.422-2.181), smoking (OR=2.798, 95%CI: 1.834-4.269), alcohol consumption (OR=2.282, 95%CI: 1.861-2.798), experiencing school bullying (OR=5.440, 95%CI: 3.148-9.402) and parental physical/verbal abuse (OR=3.954, 95%CI: 3.189-4.902) were associated with a higher risk of depressive symptoms among middle school students. Conversely, the middle school students who engaged in moderate-to-vigorous physical activity ≥3 times/week (OR=0.784, 95%CI: 0.668-0.921) and attended physical education classes ≥3 sessions/week (OR=0.736, 95%CI: 0.613-0.884) were associated with a lower risk of depressive symptoms. Conclusion The prevalence of depressive symptoms among middle school students in Huzhou City was lower than national average, and was influenced by dietary habits, physical exercise, sleep duration, smoking, alcohol consumption, and experiencing school bullying.