- 1Department of Hematology, The Affiliated Hospital of North Sichuan Medical College, Nanchong, Sichuan, China
- 2School of Psychiatry, North Sichuan Medical College, Nanchong, Sichuan, China
Introduction: Parental marital status is a pivotal determinant of adolescent maladaptive behaviors. However, the specific correlation between parental marital status and Non-Suicidal Self-Injury (NSSI) among adolescents in western China has not been thoroughly investigated.
Methods: A cross-sectional study was conducted involving 7,274 students from non-divorced families and 1,332 from divorced families. Participants were assessed using the Ottawa Self-Injury Inventory (OSI) and the Depression, Anxiety, and Stress Scale-21 Items (DASS-21). Propensity score matching (PSM) was applied to control for confounders, and logistic regression was used to analyze associations.
Results: The prevalence of non-suicidal self-injury (NSSI) was higher in adolescents from non-divorced families (4.1%) than in those from divorced families (1.6%). However, propensity score matching analysis revealed that parental divorce was associated with a significantly increased odds of NSSI after controlling for confounders. Furthermore, depression and anxiety were identified as independent risk factors for NSSI. Mediation analysis indicated that the emotional distress (depression, anxiety, and stress) did not significantly mediate the relationship between parental divorce and NSSI. Additional significant predictors of NSSI included girls and age.
Discussion: Adolescents from divorced families in western China show a higher incidence of NSSI, though this relationship is not mediated by emotional distress. These findings highlight the need for families and society to enhance attention to the psychological well-being of this vulnerable adolescent population.
Introduction
NSSI is defined as the intentional and socially unacceptable infliction of harm to one’s own body without the intent to commit suicide (1). Numerous studies have demonstrated that engagement in NSSI significantly increases the risk of suicidal ideation, suicide attempts, and a range of mental disorders (2, 3). The detrimental consequences of NSSI extend beyond the individual, profoundly impacting familial dynamics, interpersonal relationships, and in extreme cases, resulting in fatal outcomes (4). Adolescents are particularly susceptible to NSSI (5), rendering self-injury a critical global public health issue that demands immediate attention and intervention.
Globally,a previous meta-analysis indicated that globally, 8.3% of adolescents engage in occasional non-suicidal self-injury, while 20.3% exhibit repeated non-suicidal self-injury (6). In China, a meta-analysis of middle school populations reveals an aggregate NSSI prevalence of 22.4% (7), though marked regional disparities are evident. Provincial-level studies illustrate this heterogeneity: the incidence among rural primary and secondary school students in Jiangxi Province is 14.8% (8), while the lifetime prevalence of NSSI among adolescents in Yunnan Province is 47.1% (9). Notably, despite these documented variations, the epidemiological profile of NSSI in Western China’s adolescent population remains poorly characterized.
Previous studies have identified a range of factors associated with NSSI, including personal and family factors, social environment, and psychological conditions (10, 11), all of which significantly affect adolescents’ physical and psychological well-being. Research has revealed a strong link between NSSI and depression, with adolescents exhibiting depressive symptoms showing significantly higher vulnerability to NSSI compared to their non-depressed counterparts (12, 13). Furthermore, anxiety has been identified as a crucial mediating factor in the relationship between adverse life events and NSSI (14), and elevated anxiety levels often act as significant triggers for NSSI episodes (15, 16). When individuals encounter stressful life events, the likelihood of engaging in NSSI increases, particularly if they are in a negative mental health state, as NSSI can serve as a coping strategy to alleviate the negative emotions stemming from these events (17). Beyond these established factors, researchers have also explored other potential risk factors for NSSI, such as borderline personality disorder, externalizing symptoms (e.g., aggression and behavioral problems), impulsivity, self-prediction (i.e., self-reported likelihood of engaging in NSSI in the future), and gender. However, studies have reported significant variations in these factors, and intervention programs targeting known risk factors often yield only moderate results. These limitations underscore the necessity for further research to identify novel risk factors and develop more targeted prevention strategies, thereby advancing our understanding and improving clinical outcomes in NSSI management.
Divorce has emerged as a widespread and increasingly prevalent social phenomenon, with substantial variations in divorce rates observed across different countries and regions (18). In China, the divorce rate has experienced a rapid escalation in recent years, a trend that can be attributed to evolving societal attitudes, enhanced educational attainment among women, and improvements in female economic independence (19, 20). Research has consistently demonstrated that parental separation or divorce constitutes a significant risk factor for children’s psychological well-being, predisposing them to various mental health challenges (21).These include increased susceptibility to depression (21, 22), elevated risk of suicidal ideation (21), and greater vulnerability to suicide attempts. Such psychological distress may subsequently contribute to a higher propensity for NSSI. This pattern is particularly evident in South Korea, where adolescents from divorced or remarried families demonstrate significantly elevated risks of engaging in NSSI behaviors (23, 24). In Lebanon, adolescents from separated families demonstrate significantly higher levels of social anxiety, avoidance behaviors, depression, and suicidal ideation (25).A Taiwanese study further emphasized that dysfunctional family dynamics substantially increase the risk of NSSI among high school students (26). Similarly, research conducted in Hong Kong has established parental marital discord as a critical risk factor for the chronic persistence of emotional and behavioral issues in adolescents (27). These findings collectively underscore the significant influence of parental marital status on adolescent NSSI behaviors. Middle school students, who are navigating a critical phase of emotional and psychological development, are particularly vulnerable to external stressors such as academic pressure and emotional instability. During this sensitive period, the marital status of their parents can exert an even more pronounced impact on their mental health, especially in relation to the development and maintenance of NSSI behaviors.
The western region of China faces slower economic and social development compared to other areas, coupled with a substantial population of children and adolescents. However, research on NSSI among adolescents in Sichuan remains limited, particularly studies investigating the influence of parental marital status on adolescent NSSI behaviors. To address this gap, and after controlling for potential confounding variables using Propensity Score Matching (PSM), we propose the following two hypotheses regarding NSSI among adolescents in rural western China: Hypothesis 1: Adolescents from divorced families show a significantly higher rate of non-suicidal self-injury compared to those from intact families. Hypothesis 2: Emotional distress (including depression, anxiety, and stress) mediates the association between parental divorce and non-suicidal self-injury among adolescents.
Methods
Survey and inclusion criteria
This study employs a cross-sectional design to investigate the prevalence and characteristics of NSSI among adolescents. The statistical analysis utilizes a two-tailed test with a significance level (α) of 0.05 and an allowable margin of error (δ) of 0.02. Based on the most recent epidemiological data from China, which reports an adolescent NSSI prevalence rate of 22.4% (7). We calculated the minimum required sample size via PASS15 software. The initial calculation yielded a minimum sample size of 6,678 participants. To account for potential data attrition (estimated at 20%) due to incomplete responses, participant withdrawal, and other methodological factors, we adjusted the target sample size to 8,343 participants.
The study data were collected during March and April 2020 through a comprehensive survey of middle school students across 20 towns in the target region. To ensure representativeness, the study design incorporated stratification based on key demographic variables, including population size, economic development levels, and educational resource distribution. A stratified random sampling approach was implemented, with schools serving as the primary sampling unit. Specifically, one school was randomly selected from each of the 17 smaller towns, while two schools were chosen from each of the three larger towns, yielding a total of 23 participating institutions.
The research team utilized the Wenjuanxing (WJX) platform (website: https://www.wjx.cn) for data collection. After uploading the questionnaire, a unique QR code was generated for survey access. School administrators distributed this QR code to parents or legal guardians through established communication channels. Participants were required to provide informed consent before participation, with the consent form displayed as the initial page of the digital survey interface. Only after electronically signing the consent form could participants proceed to complete the questionnaire. To minimize reporting bias, an anonymous online survey was conducted. Participants were assured of the strict confidentiality of their data and were informed that the information would be used solely for scientific research purposes. This study was in accordance with the ethical principles, reviewed by the Ethics Committee of Nanchong Psychosomatic Hospital.
Data collection and measurement
The data were primarily collected from middle schools in Yingshan County, a remote region in western China. This location was selected due to its representativeness of rural areas in western China, where the majority of families live in poverty, engage in agricultural labor, and reside in either villages or nearby towns.
Data collection was conducted in two stages. In the first stage, researchers underwent comprehensive training to familiarize themselves with the study procedures, questionnaire content, and form completion protocols. They also collaborated with class teachers of the participating classes to discuss the research objectives and questionnaire content, aiming to improve response rates and ensure data quality. Following unanimous consent from all stakeholders, the questionnaires were distributed. In the second stage, two researchers collected the completed questionnaires and conducted a preliminary review to verify data accuracy. The results were exported from the WJX platform into an Excel spreadsheet for further processing. Prior to analysis, missing values, obvious errors, and inconsistent responses were systematically removed, retaining only data that met the predefined inclusion criteria for analysis.
The research team comprised practicing physicians from the Departments of Psychiatry and Psychology of Affiliated Hospital of North Sichuan Medical College, as well as graduate students from the psychiatry program at the same institution. Before initiating the survey, all team members received systematic training to ensure familiarity with the procedures, content, and subsequent steps. Throughout the study, researchers strictly adhered to Standard Operating Procedures (SOP) to maintain consistency and reliability.
The study collected 8,785 questionnaires, retaining 8,606 (98.0% validity rate) after excluding 176 based on predefined criteria: lack of consent, missing data, or inconsistent responses. Analysis focused on adolescents aged 12-18 (inclusive of 12 and 18 years old) from divorced families, with inclusion requiring: (1) age 12-18, (2) parental divorce/separation, and (3) ≥1 year since separation to ensure stable family structure and measurable psychological impact.
Demographic characteristics
A self-designed questionnaire was utilized to collect socio-demographic information from middle school students. Including age, gender, only-child status, home location (rural or urban), boarding school(in China, middle schools provide dormitories, and students may choose to board based on their family circumstances), left-behind status (defined as individuals under 18 years old who reside in their household registration location for six months or longer due to one or both parents working away and being unable to live with both parents), education level ((junior high school, senior high school)), father’s occupation, paternal education, maternal education, mother’s occupation, and parental marital status (individuals not living with both parents due to separation or divorce were categorized as having experienced parental divorce).
NSSI assessment tool
The Ottawa Self-Injury Inventory (OSI) (28), originally developed by Nxion and Paula Cloutier, is a comprehensive psychometric instrument designed to assess multiple dimensions of NSSI, including behavioral characteristics, underlying motivations, frequency patterns, and addictive components. Extensive international validation studies have established its reliability and validity as an effective assessment tool for NSSI (28). For the current study, we employed the Chinese version of the Ottawa Self-Injury Inventory (OSI Chinese version), which has demonstrated robust psychometric properties in previous research. Specifically, the Chinese adaptation has shown excellent test-retest reliability (coefficient > 0.400) and strong construct validity, with subscale internal consistency reaching a Cronbach’s α coefficient of 0.952 (29). In our sample, the questionnaire demonstrated high reliability, with a Cronbach’s α coefficient of 0.864. NSSI identification was operationalized according to the DSM-5 diagnostic criteria for NSSI, supplemented by relevant items from the OSI inventory. Participants were classified as exhibiting NSSI behaviors if they reported engaging in one or more of the following behaviors (e.g., cutting, burning, scratching, hitting, biting) within the preceding 12 months.
Mental health issues
The Depression, Anxiety, and Stress Scale (DASS-21) was employed to evaluate negative emotional states among left-behind children. The Likert scale scores range from 0 to 3 ((0 = never to 3 = almost always), This validated instrument comprises 21 items across three subscales: depression (items 3, 5, 10, 13, 16, 17, and 21), anxiety (items 2, 4, 7, 9, 15, 19, and 20), and stress (items 1, 6, 8, 11, 12, 14, and 18). The maximum and minimum scores for each dimension are 21 and 0, respectively, and they are positively correlated with the corresponding emotional levels. The internal consistency of this scale, measured by Cronbach’s alpha, was 0.956, with the subscale coefficients exceeding 0.800 (30). The Chinese version of the scale effectively reflects the emotional states of residents, university students, and adolescents in mainland China (31, 32), In this study, the alpha coefficients for the anxiety, depression, and stress dimensions of the DASS-21 were all greater than 0.830, further supporting its cross-cultural validity (33).making it a valuable instrument for assessing negative emotional states.
Statistical analyses
Propensity score matching
In this study, significant differences in demographic characteristics were observed between students from divorced and non-divorced parental groups. To address these differences and mitigate potential confounding effects on the study outcomes, we employed Propensity Score Matching (PSM) to balance the baseline characteristics between the two groups (34). PSM is a robust statistical method that minimizes discrepancies in demographic characteristics between the experimental and control groups, thereby reducing bias introduced by confounding variables (35).Specifically, the PSM procedure was implemented in three steps. First, a set of observed covariates (i.e., confounders) was identified based on theoretical and empirical considerations. Second, a Logit model was used to estimate the probability of participants being assigned to the divorced parental group, known as the propensity score. Finally, a matching algorithm was applied to pair each participant in the divorced parental group with a participant in the non-divorced parental group who had a similar propensity score. This matching process ensured that the two groups were balanced with respect to potential confounding factors, enhancing the validity of the comparative analysis.
PSM was performed using R 4.2.2 with the MatchIt, tableone, and cobalt packages. The Logit model was selected for the matching function, with a caliper set at 0.2*SD (36, 37). A greedy nearest neighbor matching approach without replacement was used for 1:1 matching. The standardized mean difference (SMD) was employed to assess the balance of demographic characteristics after matching (34, 38). To achieve a satisfactory matching balance, the absolute value of the SMD was set to be less than 0.1 (36, 38, 39). The balance of confounding factors before and after matching is illustrated using an SMD plot (Figure 1), while the distribution of propensity scores before and after matching for both the divorced and non-divorced parental groups is shown in Figures 2, 3.
Data analysis
Univariate analyses were conducted using SPSS 26.0 both before and after matching. The Kolmogorov-Smirnov test (K-S test) revealed that the data were not normally distributed. Consequently, non-parametric tests were utilized for between-group comparisons to evaluate median differences, the Wilcoxon rank-sum test was applied to compare non-normally distributed continuous and ordinal variables across groups. Additionally, binary logistic regression analysis was performed to identify potential risk factors, with model adequacy verified through the Hosmer-Lemeshow test (P > 0.05 indicating satisfactory fit), Mediation analyses were conducted using SPSS PROCESS Macro (v4.1). All statistical tests were two-tailed, with a significance level set at α = 0.05.
Results
Analysis of NSSI
In this study, 7,274 middle school students with non-divorced parents and 1,332 middle school students with divorced parents participated in the survey and completed the assessment. The prevalence of NSSI was found to be 4.1% in the non-divorced parental group, and 1.6% in the divorced parental group. There were statistically significant differences between the two groups (divorced vs. non-divorced parents) in terms of age, grade, left-behind experience, Boarding at school, only child status, residence place, Paternal educational qualification, Maternal educational qualification, father’s occupation, mother’s occupation (P<0.01, Table 1).
Potential confounder selection for propensity score matching
Based on the relationship between the grouping variables and outcome variables, potential confounding factors such as grade, left-behind experience, residence place, Paternal educational qualification, Maternal educational qualification, and age were selected as covariates and matched in the propensity score model.
Propensity score matching
The PSM procedure resulted in a final matched sample of 2,664 students, comprising 1,332 from divorced parental families and 1,332 from non-divorced parental families. After PSM, the incidence rate of NSSI was 3.4% in the non-divorced parental group and 5.0% in the divorced parental group. The standard deviation of the propensity scores for both groups was 0.0001, indicating a high level of balance achieved through the matching process.
Potential confounding factors, such as grade, left-behind experience, residence place, Paternal educational qualification, Maternal educational qualification, age, and the variables not included as potential confounders, demonstrated excellent balance (all SMD < 0.1; see Table 2). This suggests that the matching process effectively balanced all potential confounders between the two groups.
The SMD values for each potential confounding factor before and after matching are illustrated in Figure 1. After matching, all confounders exhibited excellent balance (all SMD < 0.1). Additionally, the SMD values for all confounders and their subgroups indicated strong matching both before and after PSM (all SMD < 0.1). Furthermore, the symmetry of propensity scores between the divorced and non-divorced groups improved significantly after matching, as visually demonstrated in Figures 2, 3.
Demographic and clinical characteristics of the matched sample
Table 2 shows the demographic and clinical characteristics of the matched samples after propensity score matching. Univariate analysis indicates that after matching, the divorced and non-divorced groups are well-balanced in grade, left-behind experience, residence place, Paternal educational qualification, Maternal educational qualification, and age (P > 0.05, all SMD < 0.1). Regarding the DASS-21 dimensions, namely depression (Z = 39.283, P = 0.011), anxiety (Z = 42.824, P = 0.048), and stress (Z = 30.797, P = 0.013), there were statistically significant differences between the divorced and non-divorced groups. Specifically, the divorced group exhibited higher scores in depression, anxiety, and stress compared to the non-divorced group.
Multivariable analyses
We employed binary logistic regression to examine the relationship between parental marital status and NSSI, adjusting for depression, anxiety, and stress as covariates. The results of the analysis revealed that several factors were significantly associated with NSSI among middle school students. Parental divorce (OR = 4.052, 95% CI: 2.629~6.245, P < 0.001), depression (OR = 1.061, 95% CI: 1.029~1.095, P < 0.001), anxiety (OR = 1.074, 95% CI: 1.038~1.111, P < 0.001), girls (OR = 2.463, 95% CI: 1.704~3.559, P < 0.001), and age (OR = 1.127, 95% CI: 1.006~1.262, P = 0.039) were all significantly linked to an increased likelihood of engaging in NSSI (Table 3). The Hosmer-Lemeshow test (P>0.05) confirmed that the model had a good fit.
Mediation analyses
Using the SPSS PROCESS Macro (v4.1), separate mediation models were constructed to examine the potential mediating roles of depression, anxiety, and stress. Parental divorce (YES vs. No) was included as the independent variable, and NSSI (YES vs. No) as the dependent variable, with depression, anxiety, and stress tested as individual mediators. All analyses were conducted using a propensity score-matched sample to ensure comparability between groups.
The results indicated that parental marital status did not significantly predict levels of depression (β = 0.173, P = 0.584), anxiety (β = 0.234, P = 0.423), or stress (β = 0.302, P = 0.348). However, depression (β = 0.108, P < 0.001), anxiety (β = 0.112, P < 0.001), and stress (β = 0.101, P < 0.001) each significantly positively predicted NSSI. After controlling for these emotional variables, the direct effect of parental marital status on NSSI remained significant (all P < 0.05). Furthermore, all indirect effect Bootstrap confidence intervals included zero, indicating that depression, anxiety, and stress did not show significant mediating roles in the relationship between parental marital status and NSSI.
Discussion
This study employs PSM methodology to examine the relationship between parental marital status and NSSI among Chinese adolescents in western regions, effectively addressing selection bias concerns that have limited previous observational research. Through rigorous analysis of matched samples, we found parental divorce is independently associated with NSSI risk, even after controlling for key psychological covariates. Initial comparisons in unmatched samples revealed statistically significant disparities (P < 0.05) between adolescents from divorced and non-divorced families across multiple mental health dimensions, including NSSI frequency, depressive symptoms, anxiety levels, and stress. Following rigorous PSM adjustment for potential confounders, robust between-group differences persisted, with adolescents from divorced families demonstrating a 4-fold increased risk of NSSI engagement compared to their counterparts from intact families. Interestingly, although a significant association was observed between parental divorce and increased risk of NSSI, unlike previous studies (40), the present study did not find that this relationship was mediated by elevated levels of emotional distress in adolescents. This suggests that the link between parental divorce and NSSI is more likely attributable to direct pathways or moderated by other contextual factors. Future research should further explore these potential moderating mechanisms to more comprehensively understand how family structure risks translate into individual behavioral issues. These methodologically rigorous findings not only address a critical gap in regional mental health research but also advance our understanding of family system dynamics in adolescent NSSI within China’s unique sociocultural context.
Compared to students from non-divorced families, those from divorced families exhibit a higher prevalence of NSSI behaviors (41). Divorce is the outcome of a separation process that often begins even before parents leave the household, with long-lasting effects (42). The quality of the emotional bond between parents and children plays a significant role in the onset of adolescent NSSI (43).According to Family Process Theory, factors such as post-divorce family conflict (44),chaotic parenting practices, and diminished parental monitoring (45), exert a more critical influence on adolescents’ psychological adjustment than divorce itself. After divorce, decreased parental supervision may increase the risk of NSSI (46). Middle school students from divorced families are more likely to experience both internalizing and externalizing behavioral problems than their peers from non-divorced families. Additionally, societal attitudes toward divorce remain largely negative, and adolescents from divorced families may experience higher levels of distress (21). The divorce process typically exposes both parents and children to distinct stressors at different stages, often resulting in emotional detachment and progressive marital dissatisfaction (42). This familial disruption can lead to impaired family functioning and adverse psychological consequences, potentially compromising children’s mental health outcomes (47). From an economic perspective (48), divorce often leads to reduced household financial resources, and the resulting financial strain may further contribute to risk behaviors among adolescents. Evidence suggests that adolescents from divorced families are particularly vulnerable to various psychological difficulties, including heightened anxiety (22), depression (40), aggression, stress, behavioral disturbances, substance abuse, and suicidal ideation (22). Meanwhile, social stigma and support mechanisms warrant equal attention: adolescents may experience a sense of shame associated with changes in family structure, and the breakdown of social support networks can amplify psychological distress and engagement in NSSI (41).These findings highlight the significant impact of parental marital status on NSSI behaviors among middle school students.
Our findings substantiate a significant association between NSSI and depressive symptoms, which aligns with findings from previous research (49, 50). According to the emotion regulation disorder theory, individuals who lack effective strategies for regulating negative emotions are more likely to engage in NSSI behaviors as a way to alleviate emotional pain and distress (51). Over time, persistent difficulties in emotion regulation can lead to the development of depressive symptoms, creating a vicious cycle. Klonsky (52) found that individuals who engage in NSSI typically experience higher levels of depressive symptoms, and there is a strong, bidirectional relationship between the two, with mutual reinforcement. Family environmental factors also play a crucial role in this process. Research suggests that family dysfunction can make individuals more vulnerable to emotion regulation difficulties, which in turn increases the risk of both NSSI behaviors and depressive symptoms (53). Living in a single-parent household often results in limited economic and social resources (such as support and parental supervision), which can be linked to poorer mental health outcomes in children, including depression (54) and NSSI (55). Previous studies have shown that 40.4% of adolescents from separated families in Sweden suffer from long-term depression (54), while a study in Galicia, Spain, revealed that parental separation leads to an average increase of 20% in depressive symptoms, generalized anxiety, hostility, paranoid thoughts, and interpersonal alienation (56).
Our findings corroborate previous research demonstrating a strong association between anxiety symptoms and NSSI in adolescents (22). From a cognitive-affective perspective, adolescence represents a critical developmental period characterized by rapid physical and emotional changes, during which self-identity and emotion regulation capacities are still evolving. When confronted with stressors or negative emotions, adolescents often lack effective coping mechanisms, leading to heightened anxiety (57). This emotional dysregulation, coupled with low self-directedness in managing challenges, may manifest in self-harming behaviors (55). Adolescents from divorced families appear particularly vulnerable to anxiety-related difficulties. Research conducted in Lebanon revealed that parental separation is associated with elevated social anxiety and generalized anxiety symptoms, often stemming from feelings of insecurity and environmental instability (25). These anxiety-provoking experiences may contribute to social avoidance behaviors and, ultimately, increase the risk of NSSI episodes (15, 16). Consistent with prior studies, parental divorce has been shown to foster chronic anxiety in children, characterized by persistent uncertainty and future-oriented fears (47).
In addition, we analyzed several significant risk factors associated with NSSI, with gender emerging as a particularly noteworthy correlate. Compared to boys, girls are more likely to engage in NSSI behaviors. However, the nature and magnitude of gender differences in NSSI prevalence remain a subject of ongoing academic debate. While some studies have consistently reported higher NSSI prevalence rates among female adolescents (58, 59), other investigations have found minimal or nonsignificant gender differences (60). Notably, a comprehensive meta-analysis focusing on Chinese middle school populations revealed an unexpected pattern, with boys exhibiting higher NSSI detection rates than girls (61). These inconsistent findings across studies may be attributable to methodological variations in assessment tools, sampling strategies, and cultural factors.
Moreover, our findings indicate that the average age of NSSI onset was 15 years, consistent with previous reports (62). This study further confirms the relationship between age and NSSI: the incidence of NSSI tends to increase with age, which is consistent with previous studies (24). However, some researchers suggest that the prevalence of NSSI gradually declines from adolescence into early adulthood (63). The varying results may be due to differences in age groups or the lack of longitudinal data.
Therefore, our study contributes to the literature by demonstrating that in the context of Rural Western China, the association between parental divorce and NSSI remains significant even when accounting for emotional symptoms, pointing to the necessity of exploring these alternative, more direct pathways in future research. Clinically, this suggests that interventions for adolescents from divorced families should not only focus on emotional symptom reduction but also directly address issues like family conflict and improve parenting competencies.
Strengths and limitations
This study exhibits several notable strengths. The primary methodological strength lies in the application of the PSM approach to control for potential confounding factors, including gender, age, living environment, and parental educational qualification. This approach allows for the effective matching of individuals with similar demographic characteristics. The study’s design, which incorporates class-level analysis complemented by PSM, represents a robust methodological framework that significantly enhances the accuracy and reliability of the research outcomes. Furthermore, the application of PSM in investigating NSSI among adolescents constitutes a novel methodological contribution, as this approach remains relatively underutilized in this specific research domain.
However, the study also has some limitations. First, the cross-sectional nature of the study design inherently limits the ability to establish causal relationships between variables. Second, the reliance on self-reported measures for assessing self-harm behaviors may introduce validity concerns, particularly regarding potential social desirability bias that could affect participants’ willingness to disclose sensitive information. Nevertheless, the school-based survey setting likely mitigates this bias to some extent. Third, the study’s sample size limitations precluded the ability to differentiate between non-divorced parents who have lost a spouse, although this limitation is unlikely to substantially impact the overall findings.
Additionally, the reliability of PSM is dependent on the observed covariates, and observational studies may suffer from omitted variables, leading to potential selection bias. Therefore, it is crucial to include as many relevant factors as possible in the demographic data to reduce selection bias in PSM. Despite these limitations, the study’s large sample size is a significant strength, which enhances the generalizability of the findings.
In conclusion, our findings demonstrate a significant association between parental divorce and increased risk of NSSI among Chinese adolescents, with depression, anxiety, and stress serving as important factors. These results have several important implications for both clinical practice and future research. First, they highlight the need for routine screening for psychological distress and NSSI behaviors among adolescents from divorced families in clinical settings. Second, they suggest that interventions targeting emotion regulation and coping skills may be particularly beneficial for this vulnerable population. Finally, our study underscores the importance of developing family-centered support programs that address the unique challenges faced by divorced families in China. Future longitudinal studies are needed to confirm these associations and explore potential cultural factors influencing the relationship between family structure and adolescent mental health outcomes.
Conclusion
In summary, our findings indicate that Divorced family, depression, anxiety, and stress each serve as significant and independent risk factors for NSSI among Chinese adolescents. Although no mediating effect through emotional pathways was established, depression, anxiety, and stress remain strong independent predictors of NSSI, underscoring their crucial role in prevention and intervention efforts. Importantly, adolescents from divorced families demonstrated a substantially higher risk of engaging in NSSI compared to their peers from intact families. While this association warrants further validation through longitudinal studies, our results robustly support the importance of familial structural factors in adolescent mental health. These findings highlight the necessity of incorporating targeted support systems—such as school-based psychological services and community-based programs—specifically designed for adolescents experiencing family disruptions.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Nanchong Psychosomatic Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
XR: Funding acquisition, Writing – original draft, Visualization, Conceptualization. LP: Investigation, Formal Analysis, Writing – original draft. YS: Supervision, Investigation, Writing – original draft. PH: Software, Writing – original draft, Formal Analysis, Investigation. DW: Project administration, Validation, Writing – original draft, Investigation. JL: Conceptualization, Methodology, Supervision, Visualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was financially supported by the Primary Health Development Research Center of Sichuan Province Program (Investigator Grant: SWFZ24-Y-60, SWFZ24-Y-69), Nanchong City Social Science Research ‘14th Five-Year Plan’ Project (NC24B241) and Project of Affiliated Hospital of North Sichuan Medical College (2024MPRW005).
Acknowledgments
We thank School of Psychiatry, North Sichuan Medical College for their strong support of this study, and the psychiatric students for their great contribution to this study during data collection and collation. We thank all adolescent students who participated in this study.
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 author(s) declare that no Generative 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.
Abbreviations
LBC, Left-Behind Children; NLBC, Non-Left-Behind Children; NSSI, Non-Suicidal Self-Injury; PSM, Propensity Score Matching; SMD, Standardized Mean Difference.
References
1. Brown RC and Plener PL. Non-suicidal self-injury in adolescence. Curr Psychiatry Rep. (2017) 19:20. doi: 10.1007/s11920-017-0767-9
2. Xia Q, Zhang Y, and Huang X. Psychological well-being and associated factors among parents of adolescents with non-suicidal self-injury: a cross-sectional study. Front Psychiatry. (2023) 14:1253321. doi: 10.3389/fpsyt.2023.1253321
3. Porras-Segovia A, Pascual-Sanchez A, Greenfield G, Creese HM, Saxena S, Hargreaves D, et al. Early risk factors for self-injurious thoughts and behaviours: A UK population-based study of 219,581 people. Behav Sci (Basel Switzerland). (2023) 14:16. doi: 10.3390/bs14010016
4. Mars B, Heron J, Klonsky ED, Moran P, O’Connor RC, Tilling K, et al. Predictors of future suicide attempt among adolescents with suicidal thoughts or non-suicidal self-harm: a population-based birth cohort study. Lancet Psychiatry. (2019) 6:327–37. doi: 10.1016/S2215-0366(19)30030-6
5. Kaligis F, Ismail RI, Wiguna T, Prasetyo S, Indriatmi W, Gunardi H, et al. Mental health problems and needs among transitional-age youth in Indonesia. Int J Environ Res Public Health. (2021) 18:4046. doi: 10.3390/ijerph18084046
6. Xiao Q, Song X, Huang L, Hou D, and Huang X. Global prevalence and characteristics of non-suicidal self-injury between 2010 and 2021 among a non-clinical sample of adolescents: A meta-analysis. Front Psychiatry. (2022) 13:912441. doi: 10.3389/fpsyt.2022.912441
7. Lang J and Yao Y. Prevalence of nonsuicidal self-injury in chinese middle school and high school students: A meta-analysis. Medicine. (2018) 97:e12916. doi: 10.1097/MD.0000000000012916
8. Tang J, Li G, Chen B, Huang Z, Zhang Y, Chang H, et al. Prevalence of and risk factors for non-suicidal self-injury in rural China: Results from a nationwide survey in China. J Affect Disord. (2020) 226:188–95. doi: 10.1016/j.jad.2017.09.051
9. Liu Y, Xiao Y, Ran H, He X, Jiang L, Wang T, et al. Association between parenting and non-suicidal self-injury among adolescents in Yunnan, China: a cross-sectional survey. PeerJ. (2020) 8:e10493. doi: 10.7717/peerj.10493
10. Calvo N, Lugo-Marín J, Oriol M, Pérez-Galbarro C, Restoy D, Ramos-Quiroga JA, et al. Childhood maltreatment and non-suicidal self-injury in adolescent population: A systematic review and meta-analysis. Child Abuse Negl. (2024) 157:107048. doi: 10.1016/j.chiabu.2024.107048
11. Tang J, Li G, Chen B, Huang Z, Zhang Y, Chang H, et al. Prevalence of and risk factors for non-suicidal self-injury in rural China: Results from a nationwide survey in China. J Affect Disord. (2018) 226:188–95. doi: 10.1016/j.jad.2017.09.051
12. Lei H, Yang Y, Zhu T, Zhang X, and Dang J. Network analysis of the relationship between non-suicidal self-injury, depression, and childhood trauma in adolescents. BMC Psychol. (2024) 12:234. doi: 10.1186/s40359-024-01729-2
13. Lu YL, Ge Y, Li M, Liang S, Zhang X, Sui Y, et al. Cognitive Phenotype Shifts in Risk Taking: Interplay of Nonsuicidal Self-Injury Behaviors and Intensified Depression. Biol Psychiatry Cogn Neurosci Neuroimaging. (2025) 10:504–12. doi: 10.1016/j.bpsc.2024.05.011
14. Zhang Y, Suo X, Zhang Y, Zhang S, Yang M, Qian L, et al. The relationship between negative life events and nonsuicidal self-injury among Chinese adolescents: A moderated-mediation model. Neuropsychiatr Dis Treat. (2022) 18:2881–90. doi: 10.2147/NDT.S386179
15. Robinson K, Brocklesby M, Garisch JA, O’Connell A, Langlands R, Russell L, et al. Socioeconomic deprivation and non-suicidal selfinjury in New Zealand adolescents: The mediating role of depression and anxiety. New Z J Psychol (Online). (2017) 46:126–36.
16. Tatnell R, Hasking P, Newman L, Taffe J, and Martin G. Attachment, emotion regulation, childhood abuse and assault: examining predictors of NSSI among adolescents. Arch suicide Res. (2017) 21:610–20. doi: 10.1080/13811118.2016.1246267
17. Haliczer LA and Dixon-Gordon KL. Social stressors, emotional responses, and NSSI urges and behaviors in daily life. J Affect Disord. (2023) 338:601–9. doi: 10.1016/j.jad.2023.06.043
18. Wang C-TL and Schofer E. Coming out of the penumbras: world culture and cross-national variation in divorce rates. Soc Forces. (2018) 97:675–704. doi: 10.1093/sf/soy070
19. Chen M, Rizzi EL, and Yip PS. Divorce trends in China across time and space: an update. Asian Population Stud. (2021) 17:121–47. doi: 10.1080/17441730.2020.1858571
20. Lan X, Marci T, and Moscardino U. Parental autonomy support, grit, and psychological adjustment in Chinese adolescents from divorced families. J Fam Psychol. (2019) 33:511–20. doi: 10.1037/fam0000514
21. Auersperg F, Vlasak T, Ponocny I, and Barth A. Long-term effects of parental divorce on mental health - A meta-analysis. J Psychiatr Res. (2019) 119:107–15. doi: 10.1016/j.jpsychires.2019.09.011
22. Sands A, Thompson EJ, and Gaysina D. Long-term influences of parental divorce on offspring affective disorders: A systematic review and meta-analysis. J Affect Disord. (2017) 218:105–14. doi: 10.1016/j.jad.2017.04.015
23. Park YS and Park EC. Association between parental marital status and types of suicidal behavior among Korean adolescents: A cross-sectional study. J Prev Med Public Health = Yebang Uihakhoe chi. (2020) 53:419–28. doi: 10.3961/jpmph.20.004
24. Kim H, Ryu JM, and Kim HW. Characteristics and trends of suicide attempt or non-suicidal self-injury in children and adolescents visiting emergency department. J Korean Med Sci. (2020) 35:e276. doi: 10.3346/jkms.2020.35.e276
25. Obeid S, Al Karaki G, Haddad C, Sacre H, Soufia M, Hallit R, et al. Association between parental divorce and mental health outcomes among Lebanese adolescents: results of a national study. BMC Pediatr. (2021) 21:455. doi: 10.1186/s12887-021-02926-3
26. Tang W-C, Tseng H-Y, Lin M-P, Lee Y-T, Wu JY-W, Cheng L-H, et al. Depression fully mediates the effects of problematic internet use on nonsuicidal self-injury among adolescents during the COVID-19 outbreak. J Psychiatr Res. (2024) 178:236–42. doi: 10.1016/j.jpsychires.2024.08.005
27. Shek DT and Yu L. Self-harm and suicidal behaviors in Hong Kong adolescents: prevalence and psychosocial correlates. TheScientificWorldJournal. (2012) 2012:932540. doi: 10.1100/2012/932540
28. Nixon MK, Levesque C, Preyde M, Vanderkooy J, and Cloutier PF. The Ottawa Self-Injury Inventory: Evaluation of an assessment measure of nonsuicidal self-injury in an inpatient sample of adolescents. Child Adolesc Psychiatry Ment Health. (2015) 9:26. doi: 10.1186/s13034-015-0056-5
29. Chen H, Pan B, Zhang C, Guo Y, Zhou J, and Wang X. Revision of the non-suicidal self-injury behavior scale for adolescents with mental disorder. Zhong nan da xue xue bao Yi xue ban = J Cent South Univ Med Sci. (2022) 47:301–8. doi: 10.11817/j.issn.1672-7347.2022.210549
30. Henry JD and Crawford JR. The short-form version of the Depression Anxiety Stress Scales (DASS-21): construct validity and normative data in a large non-clinical sample. Br J Clin Psychol. (2005) 44:227–39. doi: 10.1348/014466505X29657
31. Luo J, Xian G, Yao C, Jin J, Zhao J, Song X, et al. Further validating of the Depression Anxiety Stress Scales-21 (Dass-21) in Chinese adolescents: factor structure and measurement invariance. BMC Public Health. (2025) 25:2073. doi: 10.1186/s12889-025-23287-x
32. Cao CH, Dang CY, Zheng X, Chen WG, Chen IH, and Gamble JH. The psychometric properties of the DASS-21 and its association with problematic internet use among Chinese college freshmen. Healthcare (Basel Switzerland). (2023) 11:700. doi: 10.3390/healthcare11050700
33. Cao CH, Liao XL, Gamble JH, Li LL, Jiang XY, Li XD, et al. Evaluating the psychometric properties of the Chinese Depression Anxiety Stress Scale for Youth (DASS-Y) and DASS-21. Child Adolesc Psychiatry Ment Health. (2023) 17:106. doi: 10.1186/s13034-023-00655-2
34. Austin PC. A critical appraisal of propensity-score matching in the medical literature between 1996 and 2003. Stat Med. (2008) 27:2037–49. doi: 10.1002/sim.3150
35. VanderWeele TJ and Shpitser I. On the definition of a confounder. Ann Stat. (2013) 41:196–220. doi: 10.1214/12-AOS1058
36. Austin PC. Optimal caliper widths for propensity-score matching when estimating differences in means and differences in proportions in observational studies. Pharm Stat. (2011) 10:150–61. doi: 10.1002/pst.433
37. Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivariate Behav Res. (2011) 46:399–424. doi: 10.1080/00273171.2011.568786
38. Yao XI, Wang X, Speicher PJ, Hwang ES, Cheng P, Harpole DH, et al. Reporting and guidelines in propensity score analysis: A systematic review of cancer and cancer surgical studies. J Natl Cancer Inst. (2017) 109:djw232. doi: 10.1093/jnci/djw323
39. Woo A, Lee SW, Koh HY, Kim MA, Han MY, and Yon DK. Incidence of cancer after asthma development: 2 independent population-based cohort studies. J Allergy Clin Immunol. (2021) 147:135–43. doi: 10.1016/j.jaci.2020.04.041
40. Ai M, Xu XM, Wang W, Chen J, He J, Zhang Q, et al. Impact of parental marital status on self-harm in Chinese primary school students: the mediating role of depression and the moderating effect of classmate relationships. PeerJ. (2025) 13:e19307. doi: 10.7717/peerj.19307
41. Wang F, Lu J, Lin L, Cai J, Xu J, and Zhou X. Impact of parental divorce versus separation due to migration on mental health and self-injury of Chinese children: a cross sectional survey. Child Adolesc Psychiatry Ment Health. (2021) 15:71. doi: 10.1186/s13034-021-00424-z
42. Amato PR. The impact of family formation change on the cognitive, social, and emotional well-being of the next generation. Future Child. (2005) 15:75–96. doi: 10.1353/foc.2005.0012
43. Eren Sarıkaya BD, Ermiş Ç, Baykara HB, and Demirgören BS. Non-suicidal self injury: relationship with attachment, childhood trauma and dissociation. Turk psikiyatri dergisi = Turkish J Psychiatry. (2024) 36:1–7. doi: 10.5080/u27163
44. Cao H, Fine MA, and Zhou N. The divorce process and child adaptation trajectory typology (DPCATT) model: the shaping role of predivorce and postdivorce interparental conflict. Clin Child Family Psychol Rev. (2022) 25:500–28. doi: 10.1007/s10567-022-00379-3
45. Cano T and Gracia P. The gendered effects of divorce on mothers’ and fathers’ Time with children and children’s developmental activities: A longitudinal study. Eur J population = Rev europeenne demographie. (2022) 38:1277–313. doi: 10.1007/s10680-022-09643-2
46. Victor SE, Hipwell AE, Stepp SD, and Scott LN. Parent and peer relationships as longitudinal predictors of adolescent non-suicidal self-injury onset. Child Adolesc Psychiatry Ment Health. (2019) 13:1. doi: 10.1186/s13034-018-0261-0
47. D’Onofrio B and Emery R. Parental divorce or separation and children’s mental health. World Psychiatry. (2019) 18:100–1. doi: 10.1002/wps.20590
48. Smith-Etxeberria K, Corres-Medrano I, and Fernandez-Villanueva I. Parental divorce process and post-divorce parental behaviors and strategies: examining emerging adult children’s attachment-related anxiety and avoidance. Int J Environ Res Public Health. (2022) 19:10383. doi: 10.3390/ijerph191610383
49. Jiao T, Guo S, Zhang Y, Li Y, Xie X, Ma Y, et al. Associations of depressive and anxiety symptoms with non-suicidal self-injury and suicidal attempt among Chinese adolescents: The mediation role of sleep quality. Front Psychiatry. (2022) 13. doi: 10.3389/fpsyt.2022.1018525
50. Lee JY, Kim H, Kim SY, Kim JM, Shin IS, and Kim SW. Non-suicidal self-injury is associated with psychotic like experiences, depression, and bullying in Korean adolescents. Early intervention Psychiatry. (2021) 15:1696–704. doi: 10.1111/eip.13115
51. Wolff JC, Thompson E, Thomas SA, Nesi J, Bettis AH, Ransford B, et al. Emotion dysregulation and non-suicidal self-injury: A systematic review and meta-analysis. Eur Psychiatry. (2019) 59:25–36. doi: 10.1016/j.eurpsy.2019.03.004
52. Klonsky ED. The functions of deliberate self-injury: a review of the evidence. Clin Psychol Rev. (2007) 27:226–39. doi: 10.1016/j.cpr.2006.08.002
53. Bhattarai A, King N, Adhikari K, Dimitropoulos G, Devoe D, Byun J, et al. Childhood Adversity and Mental Health Outcomes Among University Students: A Longitudinal Study. Can J Psychiatry Rev Can Psychiatr. (2023) 68:510–20. doi: 10.1177/07067437221111368
54. Bohman H, Låftman SB, Päären A, and Jonsson U. Parental separation in childhood as a risk factor for depression in adulthood: a community-based study of adolescents screened for depression and followed up after 15 years. BMC Psychiatry. (2017) 17:117. doi: 10.1186/s12888-017-1252-z
55. Koposov R, Stickley A, and Ruchkin V. Non-suicidal self-injury among incarcerated adolescents: prevalence, personality, and psychiatric comorbidity. Front Psychiatry. (2021) 12:652004. doi: 10.3389/fpsyt.2021.652004
56. Seijo D, Fariña F, Corras T, Novo M, and Arce R. Estimating the epidemiology and quantifying the damages of parental separation in children and adolescents. Front Psychol. (2016) 7:1611. doi: 10.3389/fpsyg.2016.01611
57. Jiang Z, Wang Z, Diao Q, Chen J, Tian G, Cheng X, et al. The relationship between negative life events and non-suicidal self-injury (NSSI) among Chinese junior high school students: the mediating role of emotions. Ann Gen Psychiatry. (2022) 21:45. doi: 10.1186/s12991-022-00423-0
58. Hu R, Peng LL, Hong BX, Jiang LH, Deng CF, Song GY, et al. Effect of family environment on non-suicidal self-injury among middle school students during the COVID-19 epidemic: the mediating role of depression. Sichuan Da Xue Xue Bao Yi Xue Ban. (2023) 54:374–9. doi: 10.12182/20230160106
59. Jeong JY and Kim DH. Gender differences in the prevalence of and factors related to non-suicidal self-injury among middle and high school students in South Korea. Int J Environ Res Public Health. (2021) 18:5965. doi: 10.3390/ijerph18115965
60. Xin M, Yang X, Liu K, Naz Boke B, and Bastien L. Impact of negative life events and social support on nonsuicidal self-injury among Chinese middle school students. Am J men’s Health. (2020) 14:1557988320937124. doi: 10.1177/1557988320937124
61. Yang X and Feldman MW. A reversed gender pattern? A meta-analysis of gender differences in the prevalence of non-suicidal self-injurious behaviour among Chinese adolescents. BMC Public Health. (2017) 18:66. doi: 10.1186/s12889-017-4614-z
62. Gandhi A, Luyckx K, Adhikari A, Parmar D, Desousa A, Shah N, et al. Non-suicidal self-injury and its association with identity formation in India and Belgium: A cross-cultural case-control study. Transcultural Psychiatry. (2021) 58:52–62. doi: 10.1177/1363461520933759
Keywords: parental divorce, non-suicidal self-injury, mental health, middle school students, propensity score
Citation: Ren X, Pan L, Song Y, He P, Wu D and Luo J (2025) The effect of parental marital status on non-suicidal self-injurious behavior among middle school students in rural western China: a multicenter cross-sectional study based on propensity score matching method. Front. Psychiatry 16:1644937. doi: 10.3389/fpsyt.2025.1644937
Received: 11 June 2025; Accepted: 22 September 2025;
Published: 07 October 2025.
Edited by:
Yi Yesen Zhong, Peking University Health Science Center, ChinaReviewed by:
Zhenpeng Huang, Guangxi University of Chinese Medicine, ChinaDwi Indah Iswanti, Airlangga University, Indonesia
Copyright © 2025 Ren, Pan, Song, He, Wu and Luo. 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: Jiaming Luo, amlhbWluZ2x1b0Buc21jLmVkdS5jbg==
†These authors have contributed equally to this work
‡ORCID: Jiaming Luo, orcid.org/0000-0003-0025-0379