- College of Physical Education, Zhengzhou University of Light Industry, Zhengzhou, China
Background: Suicide risk (SR), including suicidal ideation and suicide attempts, is a critical public health issue, particularly among adolescents. Emerging evidence highlights substance use risk (SUR) as a prominent factor elevating SR among youth. Physical activity (PA) has demonstrated protective effects against both SR and SUR, potentially due to its broad benefits for physical and mental well-being. Despite accumulating evidence supporting these associations, the underlying mediating mechanisms remain poorly understood, especially the roles of sleep duration (SD), healthy diet (HD), and mental health (MH). Clarifying the complex interactions and sequential mediating pathways among these variables is essential for developing effective interventions to mitigate SR and SUR among adolescents.
Methods: We analyzed data from the 2023 U.S. Youth Risk Behavior Surveillance System (YRBSS), including 15,734 adolescents after excluding cases with missing key variables. PA, HD, SUR, SR, SD, and MH were measured using validated YRBSS items, with multi-item constructs modeled as latent variables. Structural equation modeling and bootstrapped mediation analyses (5,000 resamples) were conducted in R 4.5, with model fit evaluated using RMSEA (< 0.08) and CFI/TLI (> 0.95).
Results: Higher PA significantly predicted lower SR (β = −0.024, p < 0.001) and SUR (β = −0.021, p < 0.001). Notably, direct pathways from PA to SR and SUR revealed small positive effects, indicating suppression effects and emphasizing the importance of indirect mechanisms. MH emerged as a crucial mediator between PA and both SUR and SR. Additionally, SD significantly mediated the association between PA and SUR. Chain mediation analyses revealed significant sequential pathways: PA → SD → MH → SUR; PA → MH → SUR → SR; and PA → HD → MH → SR, highlighting complex protective mechanisms linking physical activity with adolescent health outcomes. Conversely, HD alone did not significantly mediate the PA–SR relationship, and the PA → SD → SUR → SR pathway was not supported.
Discussion: The beneficial effects of PA on reducing SR and SUR among adolescents primarily operate indirectly through enhanced MH and SD, with minor contributions from dietary habits. These findings underscore the importance of integrated intervention strategies targeting MH and sleep quality to maximize the protective benefits of PA for adolescents.
1 Introduction
Suicide risk (SR), encompassing both suicidal ideation and suicide attempts, represents a major public health concern (1). Although suicide is not classified as a disease, the World Health Organization (WHO) has reported that approximately 727,000 people die by suicide each year. Among individuals aged 15–29 years, suicide was the third leading cause of death overall, ranking second among females and third among males (2). Notably, WHO has indicated that over half of individuals who die by suicide are under 50 years old. In a given year, suicide accounts for approximately 1.1% of deaths worldwide, however, the fatality rate is significantly higher among individuals aged 15–29 years. From a practical perspective, each suicide death represents a personal tragedy and is estimated to indirectly affect numerous individuals, including family members, friends and the broader community (3).
These data highlight the importance of SR intervention, particularly among adolescents. Advances in our understanding of SR assessment and intervention have contributed to reductions in both SR and suicidal behavior. Researchers have identified several predictive factors associated with suicidal behavior and have compared their relative importance. From a psychiatric perspective, SR has been conceptualized as the result of three interrelated factors: perceived burdensomeness, thwarted belongingness, and acquired capability for suicide, these factors may interact to produce a heightened risk of suicide behavior (4). It is worth noting that over 90% of suicide cases involve at least one mental disorder, with major depressive disorder, substance use disorder, and anxiety disorders being the most common (5). From a comprehensive perspective, major contributing factors to SR include psychiatric disorders, childhood trauma, social deprivation, and neurobiological mechanisms (6). For adolescents, several specific risk factors have been identified, including family conflict, academic pressure, peer relationships difficulties, and mental disorders (7). It is worth noting that substance use risk (SUR) is directly associated with SR, and has gradually become a critical focus in to SR research.
SUR among adolescents encompasses behaviors of high-risk substance use, such as misuse of prescription medications, illicit drug experimentation, or injection drug use, which carry elevated potential for adverse outcomes including injury, legal consequences, dropout, or even death (8). Importantly, The Diagnostic and Statistical Manual of Mental Disorders (DSM-5) defines substance use disorder (SUD) as a chronic condition characterized by impaired control over substance use and continued use despite significant harm. Adolescents, whose brains are still in a critical developmental stage, are particularly vulnerable to such risks (9). Individuals with substance use behaviors are approximately three times more likely to experience elevated SR compared to those without such behaviors (10). Alcohol-related disorders significantly increase suicidal ideation, suicide attempt, and completed suicide (11). One investigation reported that approximately 25–33% of adolescents involved in suicide cases had a history of alcohol or other substance use. Moreover, the risk of substance use increases by approximately 30% even after a single episode of alcohol consumption (12). Moreover, a significant relationship has been observed between substance use disorder and SR among youth (13). In addition to alcohol consumption, the use of other illicit substances such as heroin and cocaine is also significantly associated with SR (14). A direct study indicated that individuals with substance us disorders significantly predicted SR, with odds ratios (ORs) ranging from 2.0 to 11.2 (15). According to the DSM-5, both SR and SUR are considered manifestations of mental behavior, or neurodevelopmental disorders (9). These risks are often accompanied by three key characteristics: tolerance, withdrawal, and compulsive use behavior, all of which are indicative of a chronic condition (16). According to the ICD-11, SUR is categorized into three stages of harm: single episodes of harmful use, harmful pattern of use, and substance dependence. These stages contribute to serious and long-term harm across physical, psychological, and social domains (17). We need to attend to these growing longer-term mental and complex consequences of SR and SUR. Identifying effective strategies to intervene in SR and SUR is essential for public health.
Insufficient physical activity (PA) has been identified as a potential risk factor for numerous serious conditions, including cardiovascular disease (18), cancer (19), osteoporosis (20), diabetes, obesity, hypertension (21), anxiety, and depression (22). Regular PA is an essential and effective intervention for a wide range of chronic conditions, and a linear relationship has been observed between PA and mental health (23). Because of PA can effectively promote overall well-being, it provides significant benefits for both physical and mental health (24). Enhancing PA levels has been regarded as a critical strategy for improving overall health (25).
Adolescent substance use is embedded within broader psychosocial systems. Problem Behavior Theory posits that adolescent health-compromising behaviors, such as substance use, often cluster within the same behavioral syndrome, whereas protective factors like sport involvement and PA can buffer these risks (26). The contemporary suicide models such as the Interpersonal Theory of Suicide and the Integrated Motivational–Volitional (IMV) model highlight the role of belongingness, emotional regulation, and self-control—factors that can be enhanced through PA—in reducing SR (4, 27). PA appears to have the potential to reduce SR and SUR due to its beneficial properties. Current research has indicated a significant relationship between PA, SR, and SUR. For example, A comprehensive review demonstrated that higher PA levels can significantly predict lower suicidal ideation (OR = 0.91, 95% CI [0.51, 0.99]) (28). A meta-analysis of 17 randomized controlled trials (RCTs) found that exercise interventions significantly reduced suicide attempts compared to inactive control groups (OR = 0.23, 95% CI [0.09, 0.67]) (29). However, the underlying mechanisms of these complex relationships remain insufficiently understood, particularly among adolescents, who are more likely to contribute key factors to associated with these risks. Therefore, we aim to further clarify the complex relationship between PA, SR, and SUR, with particular attention to identifying the significant mediators involved.
In literature, the term exercise is often used to distinguish it from general PA, as PA includes both structured, scheduled exercise and unstructured daily activities such as household chores (30). Rather than focusing on whether PA is classified as structured exercise or part of daily routines, greater emphasis should be placed on its total volume, frequency, and intensity (31). ACSM (American College of Sports Medicine) defines PA as any bodily movement produced by skeletal muscles that results in energy expenditure, and exercise is defined as a planned, structured, and repetitive form of physical activity intended to improve or maintain one or more components of physical fitness (32). In this study, we do not strictly differentiate PA and exercise, instead, PA is broadly defined as any activity that results in energy expenditure.
1.1 The direct effects of PA on SR and SUR
Previous studies have demonstrated the benefits of PA across various age groups and populations. Among preschool children, even a single bout of moderate-to-vigorous PA has been shown to improve sleep quality, enhance cognitive function, and regulate insulin sensitivity, and regular PA is associated with more pronounced improvements in these domains (33). Moreover, individuals across diverse backgrounds, including pregnant and postpartum women, people of all genders and ages, individuals with disabilities, and those seeking to reduce their risk of disease, can benefit from PA interventions (34). A substantial body of epidemiological research has demonstrated that PA can significantly enhance self-esteem and subjective well-being. Adults who engage in regular PA tend to report fewer symptoms of depression and anxiety. These findings support the view that exercise may play a preventive role in the development of mental disorders (35). The theoretical research of PA suggests that it can reduce body fat, which subsequently enhances cardiopulmonary function, improves muscular strength, and promotes skeletal health, ultimately contributing to better mental health (36). We have discussed the psychological definitions of SR and SUR, both of which can lead to serious consequences. However, accumulating evidence suggests that these risks may be mitigated through PA. Based on these findings, we hypothesize that there may be a significant relationship between PA, SR, and SUR among adolescents. Accordingly, we propose the following research question:
H2a: High levels of PA will significantly predict lower levels of SR.
H2b: High levels of PA will significantly predict lower levels of SUR.
1.2 Mediating roles of sleep, diet, and mental health
The complex mediating roles of sleep duration (SD), healthy diet (HD), and mental health (MH) in the associations between PA, SR, and SUR among adolescents warrant further examination. Although direct significant relationships have been identified among SD, HD, MH, SR and SUR, research specifically focusing on adolescents remains limited. For example, a meta-analysis has identified that insomnia can prospectively predict suicidal thoughts and behaviors (37). Sleep disturbances can also be an important risk factor for SR and suicide attempts among otherwise healthy adolescents (38). A healthy diet, such as the Mediterranean diet, may prospectively reduce depressive symptoms, and depression is one of the key antecedent factors associated with increased suicide risk (39). Moreover, individuals with depression are estimated to have a 20-fold increased risk of suicide compared to the general population. Depression is currently recognized as one of the most well-established and stable predictors of suicide (40). Depression and SUR are highly comorbid and exert bidirectional influences. On one hand, depression may lead to substance misuse through a self-medication mechanism; on the other hand, SUR can exacerbate depressive symptoms and increase the risk of suicide (41). Individuals with SUR exhibit significantly elevated suicide risk, with ORs typically ranging from 2.5 to 7.5. The risk is particularly pronounced among individuals with polysubstance use, as well as those using alcohol and opioids (42). Given the established associations among PA, SR, and SUR, it is worth exploring whether SD, HD, and MH may serve as potential mediators in these relationships. Accordingly, we propose the following research question:
H1a: SD significantly mediates the relationship between PA and SUR.
H1b: MH significantly mediates the relationship among PA, SUR and SR.
H1c: HD significantly mediates the relationship between PA and SR.
1.3 Chain mediation models linking PA to SR and SUR
Chain mediation (also called serial mediation) models posit that an independent variable influences the outcome through a sequence of mediators, implying a causal chain among the mediators themselves (43). This model allows researchers to test more complex mechanisms of indirect effects beyond simple mediation (44). In this study, we aim to elucidate the complex mediating mechanisms underlying the relationship among the key variables in the model. In addition, this study seeks to further identify the sequential mediation pathways that explain how these variables are interrelated. We propose the following research question, and present a hypothesized model (Figure 1).
RQ1a: PA is hypothesized to influence SUR through a sequential mediation pathway involving SD and MH.
RQ1b: PA is hypothesized to influence SR through a sequential mediation pathway involving HD and MH.
RQ1c: PA is hypothesized to influence SR through a sequential mediation pathway involving MH and SUR.
2 Methods
2.1 Data source and participants
The data for this study were derived from the Youth Risk Behavior Surveillance System (YRBSS), a national survey conducted in the United States. The YRBSS was established in 1991 by the Centers for Disease Control and Prevention (CDC) to investigate critical health-risk behaviors among high school students. For the present analysis, we utilized the most recent cross-sectional dataset, collected in 2023.
Institutional review boards at CDC and ICF, the survey contractor, approved the study protocol for YRBSS. All procedures were conducted in accordance with the relevant ethical guidelines and federal regulations. Data collection was carried out in compliance with applicable U.S. federal law and CDC policy. All participants provided voluntary, anonymous responses, and parental consent procedures were followed prior to data collection. The publicly available dataset can be freely accessed from the official CDC website.
Due to the questionnaire targeting to collect disparate information from adolescents, it is not suitable to construct potential variable. However, a recent test–retest study of most of the 2023 survey questions indicated substantial reliability among the questions (45). Although the official guidelines indicate that the high dataset has high reliability, we conducted additional reliability tests across multiple dimensions. The results demonstrated acceptable reliability, as reported in the following sections.
The 2023 TRBSS dataset included 254,675 respondents. Figure 2 illustrates the sample filtering process. Participants with missing values in key demographic covariates were excluded, namely gender (n = 863), age (n = 907), grade (n = 1,003), and race (n = 3,200), resulting in a remaining sample of 248,702. Subsequently, participants with missing values on core study variables were also excluded. These variables included PA, HD, SUR, SR, SD, and MH. After excluding participants with missing data on any of variables, the final analytic sample comprised 15,734 participants.
2.2 Measures
2.2.1 Physical activity
Table 1 presents the variables included in the structural equation model along with their corresponding YRBSS item numbers.
PA was assessed using four items: frequency of being PA for at least 60 min per day in the past week (Q76), frequency of physical education (PE) class attendance in a typical week (Q77), number of sports teams participated in over the past 12 months (Q78), and frequency of muscle-strengthening activities such as push-ups or weight lifting in the past 7 days (Qmusclestrength). All items were coded so that higher scores reflected higher levels of PA. For example, Q76 ranged from 1 (“0 day”) to 8 (“7 days”), with greater values indicating more frequent PA. These items were treated as indicators of the latent PA construct in the SEM. The internal consistency of the PA construct was acceptable in this study (Cronbach’s α = 0.670).
2.2.2 Healthy diet
HD was assessed using three items: frequency of fruit consumption over the past 7 days (Q69), frequency of green salad consumption over the past 7 days (Q70), and frequency of consumption of other vegetables (excluding salad, potatoes, and carrots; Q73). All items were coded such that higher scores indicated higher HD behavior. For instance, Q69 ranged from 1 (“I did not eat fruit during the past 7 days”) to 7 (“4 or more times per day”), with higher values indicating better HD behavior. These items were treated as indicators of the latent HD construct in the SEM. In the present study, the scale demonstrated acceptable internal consistency (Cronbach’s α = 0.667).
2.2.3 Substance use risk
SUR was assessed using three items: number of days using electronic vapor products in the past 30 days (Q36), number of days having at least one alcoholic drink in the past 30 days (Q42), frequency of marijuana use (Q48). All items were coded such that higher scores indicated greater SUR. For instance, Q36 ranged from 1 (“0 days”) to 7 (“all 30 days”), with higher values reflecting more frequent use. These items were treated as indicators of a latent SUR construct in the SEM. In the present study, the internal consistency was acceptable (Cronbach’s α = 0.738).
2.2.4 Suicide risk
SR was assessed using three items: whether the participant had ever seriously considered attempting suicide in the past 12 months (Q27), whether they had made a plan for attempting suicide in the past 12 months (Q28), and number of suicide attempts in the past 12 months (Q29). All items were coded such that higher scores indicated greater SR. For instance, Q27 was ranged from 1 (“NO”) to 2 (“YES”), with higher values reflecting increased risk. These items were also treated as indicators of a latent SR construct in the SEM. In the present study, the internal consistency of the scale was acceptable (Cronbach’s α = 0.931).
2.2.5 Sleep duration
SD was assessed using a single item (Q85): “On an average school night, how many hours of sleep do you get?” Responses were coded such that higher scores reflected longer and healthier SD. The response options ranged from 1 (“4 or less hours”) to 7 (“10 or more hours”), with higher values indicating better SD.
2.2.6 Mental health
MH was assessed using two items. The first item (Q26) asked: “During the past 12 months, did you ever feel so sad or hopeless almost every day for 2 weeks or more in a row that you stopped doing some usual activities?” The second item (Q84) assessed the frequency of poor MH in the past 30 days, defined as experiences of stress, anxiety, or depression: “During the past 30 days, how often was your mental health not good?” Both items were coded such that higher scores indicated poorer MH. For example, Q84 ranged from 1 (“Never”) to 5 (“Always”), with higher values reflecting greater MH problems. These items were treated as indicators of latent MH construct in the SEM.
2.3 Statistical analysis
Data processing and statistical analyses were conducted using R 4. 5. For mediation analysis, the bootstrap method was adopted based on Hayes’ recommendation. Compared to Sobel test, the bootstrap method estimates the sampling distribution directly without assuming normality, providing greater statistical power and a lower Type I error rate (46). Arguably, we employed the bootstrap method with 5,000 resamples to estimate the 95% confidence intervals of the mediation effects. All mediation analyses and model fit test were also conducted using R 4.5.
Model fit was considered acceptable when all of the following criteria were met: a non-significant chi-square (χ2) test result (47). A comparative Fit Index (CFI) and Tucker-Lewis Index (TLI) greater than 0.95. And a Root Mean Square Error of Approximation (RMSEA) being less than 0.08 (48).
3 Results
3.1 Sample characteristics
Table 2 presents the covariate variables included in this study. Among the final sample, the gender distribution was nearly balanced, with 7,672 females and 8,062 males. Age was categorized into seven groups, ranging from 12 or younger to 18 or older. The distribution approximated normality, with the majority of participants aged 15 (24.9%), 16 (26.0%), and 17 (25.2), which aligns with the typical adolescent age range.
Participants were primarily in grades 9 through 12, with a relatively even distribution across grade levels and a mean grade level of 2.49. Regarding race, White participants comprised the largest group (48.7%), followed by Hispanic/Latino, Black or African American, and other racial categories. We reported frequencies, percentages, and means with standard deviations to provide a comprehensive overview of sample characteristics. All variables included in the correlation matrix were either continuous or ordinal and were treated as approximately continuous.
Table 3 presents the means, standard deviations, and intercorrelations among the key study variables, including gender, age, grade level, PA, HD, SUR, SR, SD and MH. PA was significantly correlated with HD (r = 0.26, p < 0.001), SUR (r = 0.02, p < 0.05), SR (r = −0.10, p < 0.001), SD (r = 0.16, p < 0.001) and MH (r = 0.03, p < 0.001). These results indicate that higher levels of PA are associated with higher HD, longer SD, better MH and lower SR. Although a weak positive correlation was observed between PA and SUR, this unexpected association warrants further analysis to clarify the complex relationships among these variables. Race was not included in Table 3, as it is a nominal variable lacking a meaningful numeric scale for correlation analysis. These findings support the hypothesized associations and provide empirical justification for conducting subsequent mediation analysis. While several covariates (e.g., gender, age, grade) showed modest correlations with the main variables, these effects were small and are not considered potential confounders.
3.2 Measurement model evaluation
The YRBSS is designed as a surveillance tool that provides observed behavioral indicators rather than latent constructs, and it does not require factor validation such as AVE or CR (49). However, in order to increase the measuring accuracy, we still tried to construct a latent variable (e.g., mental health, substance use risk) in line with prior SEM research and to enhance construct validity.
The measurement model demonstrated acceptable reliability and convergent validity. Standardized factor loadings ranged from 0.34 to 0.84, with most values above the recommended 0.50 threshold (50). Composite reliability (CR) ranged from 0.675 to 0.795, and AVE ranged from 0.401 to 0.568. Although AVE values for PA and Healthy Diet were slightly below the 0.50 criterion, their CR values remained above 0.60, supporting convergent validity (51). Overall, the results confirmed that the latent constructs were measured reliably (Table 4).
3.3 Model fit
We used maximum likelihood estimation with robust standard errors (MLR) to handle non-normality. All items were retained without parceling to avoid information loss. Model modifications were not applied to maintain theoretical purity. The model demonstrated the overall model acceptable fit: χ2 (93) = 1206.60, p < 0.001, RMSEA = 0.0028, 95% CI [0.026, 0.031], CFI = 0.981, TLI = 0.975, SRMR = 0.023. Although the chi-square statistic was significant, which is common in models with large sample sizes and complex structure, alternative fit indices fell within acceptable ranges. The alternative fit indices fell within conventionally accepted thresholds, indicating satisfactory model fit and construct validity.
3.4 Structural model
Table 5 presents the standardized path coefficients, standard errors, 95% confidence intervals, and significant levels for the direct effects among the key study variables. PA demonstrated significant direct associations with SR, MH, SD, SUR, HD, with all p-values less than 0.001. Specifically, PA significantly predicted better MH (β = −0.241, 95% CI [0.036, 0.077]), longer SD (β = 0.177, 95% CI [0.159, 0.194]) and higher HD (β = 0.395, 95% CI [0.375, 0.415]).
We observed that two the standardized direct paths showed counterintuitive directions: PA significantly predicted SR (β = 0.047, 95% CI [0.036, 0.077]) and SUR (β = 0.085, 95% CI [0.064, 0.106]). Although the direct effect of PA on SR was small and positive, the strong indirect pathway through improved mental health (β = −0.241) and its substantial predictive power on suicide risk (β = 0.739) suggest that the total effect is largely driven by the indirect mechanisms. The pattern is commonly observed in chain mediation models, particularly when suppressor effects or competing pathways exist. A similar explanation may apply to the relationship between PA and SR.
MH was positively associated with both SR (β = 0.395, 95% CI [0.724, 0.754], p < 0.001) and SUR (β = 0.348, 95% CI [0.326, 0.371], p < 0.001). In addition, SD significantly predicted both better MH (β = −0.303, 95% CI [−0.321, −0.286], p < 0.001) and lower SUR (β = −0.064, 95% CI [−0.084, −0.045], p < 0.001). HD was positively associated with MH (β = 0.059, 95% CI [0.033, 0.084], p < 0.001). SUR was a significant positive predictor to SR (β = 0.114, 95% CI [0.091, 0.137], p < 0.001). However, the direct path from HD to SR was not statistically significant (β = 0.016, 95% CI [−0.005, 0.037], p = > 0.05). Nevertheless, this path was retained in the mediation analysis, as on indirect effects can still be significant even when the constituent path is not significant (44).
These results suggest that the effect of PA on SR and SUR may operate primarily through indirect pathways (e.g., via MH, SD, HD), providing initial support for the hypothesized mediation structure to be tested in the subsequent analysis.
Table 6 presents the standardized direct, indirect, and total effects from PA to SUR. We observed that the total effect (β = −0.021, 95% CI [−0.037, −0.006], p < 0.001) shown that higher levels of PA can significantly predict lower SUR, supporting H2a.
The mediation analysis revealed several significant pathways. RQ1a (PA → MH → SUR) was supported, indicating that MH significantly mediated the effect of PA on SUR (β = −0.062, 95% CI [−0.070, −0.055], p < 0.001). More importantly, the pathway (PA → SD → SUR) was also supported, demonstrating that SD significantly mediated the relationship between PA and SUR (β = −0.009, 95% CI [−0.012, −0.006], p < 0.001). In addition, a significant chain mediation pathway (PA → SD → MH → SUR) was observed (β = −0.014, 95% CI [−0.016, −0.012], p < 0.001), highlighting a complex multi-step mechanism through which PA influences SUR. These findings suggest that lower levels of PA may initially lead to deteriorated MH or reduced SD, which in turn contribute to increased SUR. The Chain mediation model illustrates that insufficient PA may negatively impact SD, which subsequently impairs MH, ultimately leading to elevated levels of SUR. Together, these results emphasize that both MH and SD are the critical mediating mechanisms linking PA to SUR among adolescents.
Table 7 presents the standardized direct, indirect, and total effects from PA to SR. The total effect (β = −0.024, 95% CI [−0.028, −0.021], p < 0.001) demonstrates that higher levels of PA significantly predict lower levels of SR, thereby supporting H2b.
The mediation and chain mediation analyses revealed several key pathways linking PA to SR. RQ1b was supported, indicating that MH significantly mediated the effect of PA on SR (β = −0.03, 95% CI [−0.033, −0.027], p < 0.001). This suggesting that higher levels of PA may enhance MH, thereby reducing SR. The pathway from PA to SR visa HD was not supported (β = 0.001, 95% CI [0.000, 0.003], p < 0.001), indicating that HD did not serve as a significant mediator. However, the indirect effect via SUR was significant (β = 0.002, 95% CI [0.001, 0.002], p < 0.001), suggesting that PA may reduce SUR, which in turn decreases SR. A significant chain mediation pathway (PA → MH → SUR → SR) was also identified (β = −0.002, 95% CI [−0.002, −0.001], p < 0.001), indicating that higher levels of PA can improve MH, which subsequently lowers SUR, ultimately leading to minimize SR. Another chain pathway (PA → HD → MH → SR) was supported (β = 0.003, 95% CI [0.002, 0.004], p < 0.001), highlighting the combined mediating role of HD and MH. Similarly, the sequential pathway (PA → SD → MH → SR) was also significant (β = −0.007, 95% CI [−0.007, −0.006], p < 0.001), suggesting that PA promotes SD, which enhances MH and consequently reduces SR. In contrast, the pathway (PA → SD → SUR → SR) although statistically significant (p < 0.001), yielded a near-zero effect size (β ≈ 0.000 95% CI [0.000, 0.000], p < 0.001), indicating a lack of practical significance (Figure 3).
4 Discussion
4.1 Summary of key findings
We examined the direct and indirect relationships between PA, SR, SUR among adolescents, focusing specifically on the mediating roles of SD, HD, and MH. The results generally supported the proposed mediation and chain mediation hypotheses, though with some notable exceptions.
Consistent with hypotheses, higher levels of PA significantly predicted lower overall levels of both SUR and SR. However, direct pathways from PA to SR and SUR were unexpectedly positive, suggesting the presence of suppressor effects and highlighting the critical role of indirect pathways. Indeed, MH emerged as a central mediator, significantly mediating the relationships between PA and both SUR and SR. Additionally, SD was identified as a significant mediator linking PA to SUR, indicating that adequate sleep contributes to reduced substance use risk.
Chain mediation analyses further underscored the complexity of these relationships. PA influenced SUR through a sequential mediation pathway involving SD and MH, highlighting a stepwise protective mechanism (PA → SD → MH → SUR). Similarly, PA influenced SR through multiple significant sequential mediation pathways, including PA → MH → SUR → SR, and PA → SD → MH → SR, reinforcing the central role of MH as a pivotal mediating factor. The sequential mediation pathway via HD and MH (PA → HD → MH → SR) was statistically significant but presented a minor positive effect. Conversely, the sequential mediation via SD and SUR (PA → SD → SUR → SR) were not supported.
These findings illustrate that the beneficial impacts of PA on SR and SUR among adolescents predominantly operate indirectly through improved MH, increased SD, and partially via HD practices. These mechanisms suggest targeted avenues for intervention strategies designed to enhance adolescent risk.
4.2 Interpretation of individual hypotheses
Hypothesis H2a was supported, as higher levels of PA significantly predicted lower SUR (β = −0.021, p < 0.001). This finding aligns with previous studies suggesting physical activity’s protective effects against substance use, possibly through behavioral regulation and healthy lifestyle promotion. For instance, a longitudinal study of adolescents demonstrated that higher engagement in regular PA was associated with reduced likelihood of substance abuse behaviors, supporting our observed protective relationship (52).
Hypothesis H2b was also confirmed; PA significantly predicted reduced SR (β = −0.024, p < 0.001). Although the direct pathway between PA and SR was unexpectedly small and positive (β = 0.047, p < 0.001), the robust negative indirect effects primarily via MH indicate a suppression effect. Previous research supports this indirect protective pathway, emphasizing that PA reduces depression and anxiety, thereby lowering SR (53).
Regarding mediation hypotheses, H1a was supported as SD significantly mediated the relationship between PA and SUR (β = −0.009, p < 0.001). Adequate sleep has been consistently documented as an essential protective factor against adolescent substance use, likely by promoting emotional regulation and reducing impulsivity (54).
Hypothesis H1b was strongly supported, highlighting MH as a significant mediator in the relationships among PA, SUR, and SR (β ranging from −0.03 to −0.062, p < 0.001). Consistent with earlier findings, our results underline MH as a central explanatory factor linking physical activity to adolescent mental health outcomes, including SUR (55) and SR (56).
Hypothesis H1c regarding HD mediating effect of PA on SR was not supported (β = 0.001, p < 0.001). Although HD was positively associated with improved MH, the direct mediating role of HD alone was negligible. This outcome aligns with studies suggesting dietary effects on mental health are subtle and often operate indirectly via complex biopsychosocial pathways rather than simple, direct mediation (57).
Chain mediation hypotheses (RQ1a) received support. Revealing a significant sequential pathway from PA through SD and MH to SUR (β = −0.014, p < 0.001), suggesting a cumulative protective process through improved sleep and subsequent mental health enhancement. This finding is congruent with existing literature emphasizing multi-step health behavior chains among adolescents (58).
RQ1b was also confirmed, demonstrating significant chain mediation through HD and MH to SR (β = −0.007, p < 0.001). This result emphasizes the complex interplay among behavioral health factors in reducing suicide risk, consistent with the previous theoretical model proposed (59).
RQ1c, involving chain mediation via MH and SUR (β = −0.002, p < 0.001), was statistically significant, suggesting that complex interplay among mental health factors in reducing suicide risk. This finding is congruent with existing literature emphasizing multi-step health behavior chains among adolescents (60).
These interpretations underline the central role of indirect pathways—particularly mental health—as the primary mechanisms through which physical activity influences critical adolescent health outcomes.
4.3 Interpretation of mediation effects
The mediation analyses provide robust support for indirect pathways linking PA to SUR and SR, highlighting the critical roles of MH and SD. MH was the most influential mediator, significantly mediating relationships between PA and both SUR and SR. Adolescents with higher PA typically experience improved mental health, consequently reducing SUR and SR, consistent with prior findings emphasizing MH as a key intervention target. SD also significantly mediated the relationship between PA and SUR, reinforcing prior research that adequate sleep contributes to better emotional regulation and reduced substance use among adolescents.
However, the mediating role of HD between PA and SR was negligible, suggesting dietary influences alone may have limited practical significance. Chain mediation analyses further supported multi-step pathways (PA → SD → MH → SUR; PA → MH → SUR → SR), underscoring the interconnectedness of these health behaviors in reducing risk.
4.4 Theoretical contributions
The results indicated several important theoretical contributions to the literature on adolescent health behaviors and risk prevention. It enriches the understanding of the mechanisms linking PA to SR and SUR by highlighting the central mediating roles of MH and SD. The identification of MH as a pivotal mediator underscores its significance within adolescent behavioral health models, reinforcing theories that emphasize emotional and psychological regulation as critical pathways through which lifestyle behaviors impact health outcomes.
By explicitly examining sequential mediation (chain mediation) pathways involving SD, MH, and SUR, this study advances current conceptualizations of health behavior interconnectivity among adolescents. It demonstrates empirically that PA exerts its protective effects through complex, multi-step processes, emphasizing an integrated approach to health promotion rather than isolated behavioral interventions.
In conclusion, the complex positive direct associations observed between PA and both SR and SUR contribute to the theoretical discourse on suppressor effects and competing mechanisms in behavioral health research. These findings highlight the necessity of considering indirect and suppressor effects within health behavior theories and interventions, urging related practitioners to adopt a nuanced understanding of how protective behaviors may sometimes interact with contextual factors in complex ways.
4.5 Practical implications
The findings provide important practical implications for practitioners, educators, and policymakers aiming to promote adolescent well-being. The identification of MH and SD as key mediators between PA and adolescent risk behaviors underscores specific target areas for intervention programs.
Interventions designed to enhance MH among adolescents through regular PA should be prioritized. The related institutions could implement structured PA integrated with emotional and psychological support to foster comprehensive health benefits. Additionally, given SD protective role against SUR, educational programs emphasizing healthy sleep hygiene and routines may effectively complement existing SUR prevention efforts.
Moreover, recognizing the chain mediation effect (PA → SD → MH → SUR) highlights the value of adopting integrated intervention strategies. Practitioners should consider multi-component programs that simultaneously address PA levels, SD, and MH to optimize reductions in SUR and SR among adolescents.
4.6 Limitations
Despite its contributions, the current study has several limitations. First, the cross-sectional design prevents causal inferences from being established, as it does not allow for the determination of temporal sequence among variables. And family environment factors (e.g., parental supervision, economic status) and social support (e.g., partner relationship) may simultaneously influence PA, SD, MH, and SUR. Longitudinal or experimental designs are required to verify the directionality and causal relationships identified in the model. Second, the generalizability of our findings is limited by the sample composition, which included only U.S. adolescents. The lack of data from developing countries and some minority groups may restrict the universality of the conclusions. Future studies should validate these relationships in more diverse, multi-ethnic, and cross-national samples to enhance external validity. Third, the YRBSS does not collect detailed physical activity dose metrics, such as total weekly energy expenditure (e.g., MET-min/week), which limits our ability to examine potential dose–response or threshold effects, including the possibility that excessive activity could have adverse impacts on sleep or mental health. Finally, employing the identified model as a basis for designing and implementing longitudinal intervention studies could yield more definitive insights into the pathways connecting physical activity and adolescent health outcomes.
5 Conclusion
(1) PA significantly reduced SR and SUR, primarily through improved MH and increased SD, highlighting these factors as critical mediators. (2) Although HD alone did not mediate the relationship between PA and SR, it showed modest effects in sequential mediation pathways. (3) Significant chain mediation pathways were identified, demonstrating that PA enhances sleep duration, which subsequently improves mental health, ultimately reducing SUR and SR.
These findings emphasize that MH and SD play crucial indirect roles in linking PA with adolescent health outcomes, suggesting integrated strategies promoting PA, MH support, and sleep quality as effective approaches for reducing SR and SUR among adolescents. Although our study does not directly evaluate intervention programs, these results provide a theoretical basis for developing and testing multi-component approaches in future experimental or longitudinal research.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.cdc.gov/yrbs/about/index.html.
Ethics statement
This study analyzed publicly available, de‐identified data from the Youth Risk Behavior Surveillance System (YRBSS), collected by the U.S. Centers for Disease Control and Prevention (CDC) and its survey contractor (ICF). The original data collection was approved by the CDC and relevant institutional review boards, and parental permission (or parental consent procedures) was obtained at that time as per local requirements. The present secondary analysis was exempt from additional ethical review since the data are anonymous and publicly accessible.
Author contributions
FT: Conceptualization, Data curation, Investigation, Methodology, Project administration, Resources, Software, Writing – original draft, Writing – review & editing. FL: Validation, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported in part by grants from the General project of humanities and social science research in universities of Henan Province (2025-ZZJH-281), Doctor Research Fund from Zhengzhou University of Light Industry (2022BSJJSK12), 2024 Henan Province Philosophy and Social Sciences Planning Annual Project (2024CTY038).
Acknowledgments
The raw data supporting the conclusions of this article is available from YRBSS open website, anybody can check and download the data freely.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Turecki, G, Brent, DA, Gunnell, D, O’Connor, RC, Oquendo, MA, Pirkis, J, et al. Suicide and suicide risk. Nat Rev Dis Primers. (2019) 5:74. doi: 10.1038/s41572-019-0121-0
2. World Health Organization. Suicide worldwide in 2021: Global health estimates. Geneva: World Health Organization (2022).
3. Cerel, J, Brown, MM, Maple, M, Singleton, M, van de Venne, J, Moore, M, et al. How many people are exposed to suicide? Not six. Suicide Life Threat Behav. (2019) 49:529–34. doi: 10.1111/sltb.12450
5. Cavanagh, JTO, Carson, AJ, Sharpe, M, and Lawrie, SM. Psychological autopsy studies of suicide: a systematic review. Psychol Med. (2003) 33:395–405. doi: 10.1017/S0033291702006943
6. Turecki, G, and Brent, DA. Suicide and suicidal behaviour. Lancet. (2016) 387:1227–39. doi: 10.1016/S0140-6736(15)00234-2
7. Bridge, JA, Goldstein, TR, and Brent, DA. Adolescent suicide and suicidal behavior. J Child Psychol Psychiatry. (2006) 47:372–94. doi: 10.1111/j.1469-7610.2006.01615.x
9. American Psychiatric Association. Diagnostic and statistical manual of mental disorders. Arlington, VA: American Psychiatric Publishing (2013).
10. Favril, L, Yu, R, Uyar, A, Sharpe, M, and Fazel, S. Risk factors for suicide in adults: systematic review and meta-analysis of psychological autopsy studies. Evid Based Ment Health. (2022) 25:148–55. doi: 10.1136/ebmental-2022-300549
11. Darvishi, N, Farhadi, M, Haghtalab, T, and Poorolajal, J. Alcohol-related risk of suicidal ideation, suicide attempt, and completed suicide: a meta-analysis. PLoS One. (2015) 10:e0126870. doi: 10.1371/journal.pone.0126870
12. Canadian Centre on Substance Use and Addiction. Forging connections, igniting change: 2023–2024 annual report. Ottawa, ON: CCSA (2024).
13. Rioux, CR, Huet, A-S, Castellanos-Ryan, N, Fortier, L, Le Blanc, M, Hamaoui, S, et al. Substance use disorders and suicidality in youth: a systematic review and meta-analysis with a focus on the direction of the association. PLoS One. (2021) 16:e0255799. doi: 10.1371/journal.pone.0255799
14. Pérez Sosa, MZ, de la Vega-Sánchez, D, Sanz-Gómez, S, Melvin, G, Sánchez-Gómez, A, Guija, JA, et al. Suicide risk among Cannabis, heroin and cocaine users. Int J Ment Health Addict. (2025) 6:1–10. doi: 10.1007/s11469-025-01460-y
15. Lynch, FL, Peterson, EL, Lu, CY, Hu, Y, Rossom, RC, Waitzfelder, BE, et al. Substance use disorders and risk of suicide in a general US population: a case-control study. Addict Sci Clin Pract. (2020) 15:14. doi: 10.1186/s13722-020-0181-1
16. National Institute on Drug Abuse. Drugs, brains, and behavior: the science of addiction. Bethesda, MD: NIDA (2020).
17. World Health Organization. International classification of diseases 11th revision (ICD-11). Geneva: WHO (2019).
18. Blair, SN, Kohl, HW, Gordon, NF, and Paffenbarger, RS. How much physical activity is good for health? Annu Rev Public Health. (1992) 13:99–126. doi: 10.1146/annurev.pu.13.050192.000531
19. Lee, IM. Physical activity and cancer prevention—data from epidemiologic studies. Med Sci Sports Exerc. (1995) 27:975–81. doi: 10.1249/00005768-199507000-00003
20. Vuori, IM. Exercise and physical health: musculoskeletal health and functional capabilities. Res Q Exerc Sport. (1995) 66:276–85. doi: 10.1080/02701367.1995.10607915
21. Bouchard, C, and Deprés, J-P. Physical activity and health: atherosclerotic, metabolic, and hypertensive diseases. Res Q Exerc Sport. (1995) 66:268–75. doi: 10.1080/02701367.1995.10607914
22. Biddle, SJH. Exercise and psychosocial health. Res Q Exerc Sport. (1995) 66:292–7. doi: 10.1080/02701367.1995.10607917
23. Warburton, DER, and Bredin, SSD. Reflections on physical activity and health: what should we recommend? Can J Cardiol. (2016) 32:495–504. doi: 10.1016/j.cjca.2016.01.024
24. Warburton, DER, Nicol, CW, and Bredin, SSD. Health benefits of physical activity: the evidence. CMAJ. (2006) 174:801–9. doi: 10.1503/cmaj.051351
25. Paterson, DH, Jones, GR, and Rice, CL. Ageing and physical activity: evidence to develop exercise recommendations for older adults. Appl Physiol Nutr Metab. (2007) 32:S69–S108. doi: 10.1139/H07-111
26. Jessor, R. The origins and development of problem behavior theory In: SM Reich, K Subrahmanyam, and G Espinoza, editors. Advancing responsible adolescent development. Cham: Springer (2016). 155–80.
27. O’Connor, RC. The integrated motivational–volitional model of suicidal behavior. Crisis. (2011) 32:295–8. doi: 10.1027/0227-5910/a000120
28. Vancampfort, D, Hallgren, M, Firth, J, Rosenbaum, S, Schuch, FB, Mugisha, J, et al. Physical activity and suicidal ideation: a systematic review and meta-analysis. J Affect Disord. (2018) 225:438–48. doi: 10.1016/j.jad.2017.08.070
29. Fabiano, N, Gupta, A, Fiedorowicz, JG, Firth, J, Stubbs, B, Vancampfort, D, et al. The effect of exercise on suicidal ideation and behaviors: A systematic review and meta-analysis of randomized controlled trials. J Affect Disord. (2023) 330:355–66. doi: 10.1016/j.jad.2023.02.071
30. Caspersen, CJ, Powell, KE, and Christenson, GM. Physical activity, exercise, and physical fitness: definitions and distinctions for health-related research. Public Health Rep. (1985) 100:126–31.
31. Taylor, D. Physical activity is medicine for older adults. Postgrad Med J. (2014) 90:26–32. doi: 10.1136/postgradmedj-2012-131366
32. American College of Sports Medicine. ACSM'S guidelines for exercise testing and prescription. Netherlands: Wolters Kluwer (2021).
33. Piercy, KL, Troiano, RP, Ballard, RM, Carlson, SA, Fulton, JE, Galuska, DA, et al. The physical activity guidelines for Americans. JAMA. (2018) 320:2020–8. doi: 10.1001/jama.2018.14854
34. U.S. Department of Health and Human Services. Physical Activity Guidelines for Americans. Washington, DC: USDHHS (2008).
35. van Minnen, A, Hendriks, L, and Olff, M. When do trauma experts choose exposure therapy for PTSD patients? A controlled study of therapist and patient factors. Behav Res Ther. (2010) 48:312–20. doi: 10.1016/j.brat.2009.12.003
36. Zhou, X, Li, J, and Jiang, X. Effects of different types of exercise intensity on improving health-related physical fitness in children and adolescents: a systematic review. Sci Rep. (2024) 14:14301. doi: 10.1038/s41598-024-64830-x
37. Pigeon, WR, Pinquart, M, and Conner, K. Meta-analysis of sleep disturbance and suicidal thoughts and behaviors. J Clin Psychiatry. (2012) 73:e1160–7. doi: 10.4088/JCP.11r07586
38. Harris, LM, Huang, X, Linthicum, KP, Bryen, CP, Ribeiro, JD, et al. Sleep disturbances as risk factors for suicidal thoughts and behaviours: a meta-analysis of longitudinal studies. Sci Rep. (2020) 10:13888. doi: 10.1038/s41598-020-70866-6
39. Lassale, C, Batty, GD, Baghdadli, A, Jacka, F, Sánchez-Villegas, A, Kivimäki, M, et al. Healthy dietary indices and risk of depressive outcomes: a systematic review and meta-analysis of observational studies. Mol Psychiatry. (2019) 24:965–86. doi: 10.1038/s41380-018-0237-8
40. Harris, EC, and Barraclough, B. Suicide as an outcome for mental disorders: a meta-analysis. Br J Psychiatry. (1997) 170:205–28. doi: 10.1192/bjp.170.3.205
41. Swendsen, J, and Merikangas, KR. The comorbidity of depression and substance use disorders. Clin Psychol Rev. (2000) 20:173–89. doi: 10.1016/S0272-7358(99)00026-4
42. Armoon, B, SoleimanvandiAzar, N, Fleury, MJ, Noroozi, A, Bayat, AH, Mohammadi, R, et al. Prevalence, sociodemographic variables, mental health condition, and type of drug use associated with suicide behaviors among people with substance use disorders: a systematic review and meta-analysis. J Addict Dis. (2021) 39:1–29. doi: 10.1080/10550887.2021.1912572
43. Kessler, RC, McGonagle, KA, Zhao, S, Nelson, CB, Hughes, M, Eshleman, S, et al. Lifetime and 12-month prevalence of DSM-III-R psychiatric disorders in the United States: results from the National Comorbidity Survey. Arch Gen Psychiatry. (1994) 51:8–19. doi: 10.1001/archpsyc.1994.03950010008002
44. Hayes, AF. Introduction to mediation, moderation, and conditional process analysis: a regression-based approach. New York, NY: Guilford Press (2017).
45. Jones, SE, Brener, ND, Queen, B, Hershey-Arista, M, Harris, WA, Mpofu, JJ, et al. Reliability of the 2021 national youth risk behavior survey questionnaire. Am J Health Promot. (2024) 38:843–51. doi: 10.1177/08901171241239735
46. Hayes, AF. Beyond baron and Kenny: statistical mediation analysis in the new millennium. Commun Monogr. (2009) 76:408–20. doi: 10.1080/03637750903310360
47. Jöreskog, K, and Sörbom, D. Structural equation modelling: guidelines for determining model fit. New York, NY: University Press of America (1993).
48. Hu, LT, and Bentler, PM. Cutoff criteria for fit indexes in covariance structure analysis: conventional criteria versus new alternatives. Struct Equ Model Multidiscip J. (1999) 6:1–55. doi: 10.1080/10705519909540118
49. Centers for Disease Control and Prevention (CDC). Youth risk behavior survey data summary & trends report: 2011–2021. Department of Health and Human Services: United States (2023).
50. Hair, JF, Black, WC, Babin, BJ, and Anderson, RE. Multivariate data analysis. New York, NY: Pearson (2010).
51. Fornell, C, and Larcker, DF. Evaluating structural equation models with unobservable variables and measurement error. J Mark Res. (1981) 18:39–50.
52. Ashdown-Franks, G, Firth, J, Carney, R, Carvalho, AF, Hallgren, M, Koyanagi, A, et al. Exercise as medicine for mental and substance use disorders: a meta-review of the benefits for neuropsychiatric and cognitive outcomes. Sports Med. (2020) 50:151–70. doi: 10.1007/s40279-019-01187-6
53. Michael, SL, Lowry, R, Merlo, C, Cooper, AC, Hyde, ET, and McKeon, R. Physical activity, sedentary, and dietary behaviors associated with indicators of mental health and suicide risk. Prev Med Rep. (2020) 19:101153. doi: 10.1016/j.pmedr.2020.101153
54. Geoffroy, PA, Tebeka, S, Blanco, C, Dubertret, C, and Le Strat, Y. Shorter and longer durations of sleep are associated with an increased twelve-month prevalence of psychiatric and substance use disorders: findings from a nationally representative survey of US adults (NESARC-III). J Psychiatr Res. (2020) 124:34–41. doi: 10.1016/j.jpsychires.2020.02.009
55. Galaj, E, Barrera, ED, and Ranaldi, R. Therapeutic efficacy of environmental enrichment for substance use disorders. Pharmacol Biochem Behav. (2020) 188:172829. doi: 10.1016/j.pbb.2019.172829
56. Grasdalsmoen, M, Eriksen, HR, Lønning, KJ, and Sivertsen, B. Physical exercise, mental health problems, and suicide attempts in university students. BMC Psychiatry. (2020) 20:175. doi: 10.1186/s12888-020-02529-0
57. Dalwood, P, Marshall, S, Burrows, TL, McIntosh, A, and Collins, CE. Diet quality indices and their associations with health-related outcomes in children and adolescents: an updated systematic review. Nutr J. (2020) 19:118. doi: 10.1186/s12937-020-00637-3
58. Gambardella, J, Morelli, MB, Wang, XJ, and Santulli, G. Pathophysiological mechanisms underlying the beneficial effects of physical activity in hypertension. J Clin Hypertens (Greenwich). (2020) 22:291. doi: 10.1111/jch.13744
59. Firth, J, Solmi, M, Wootton, RE, Vancampfort, D, Schuch, FB, Hoare, E, et al. A meta-review of “lifestyle psychiatry”: the role of exercise, smoking, diet and sleep in the prevention and treatment of mental disorders. World Psychiatry. (2020) 19:360–80. doi: 10.1002/wps.20773
Keywords: physical activity, suicide risk, substance use risk, behavior mechanism, emotional mechanisms
Citation: Tian F and Liu F (2025) Physical activity and risk behaviors among adolescents: the mediating roles of diet, sleep, mental health, and substance use. Front. Public Health. 13:1665023. doi: 10.3389/fpubh.2025.1665023
Edited by:
Francesco Maria Semeraro, University of Studies G. d'Annunzio Chieti and Pescara, ItalyReviewed by:
Mauro Ceccanti, Sapienza University of Rome, ItalyMing Gao, Shanxi Medical University, China
Copyright © 2025 Tian and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fengbo Liu, bGZienp1bGlAMTI2LmNvbQ==