- 1Liaison Psychiatry, “Pius Brinzeu” County Emergency Hospital, Timisoara, Romania
- 2Department of Functional Sciences, Center for Modeling Biological Systems and Data Analysis, “Victor Babes” University of Medicine and Pharmacy, Timisoara, Romania
- 3Discipline of Psychiatry, Department of Neuroscience, NEUROPSY-COG Center for Cognitive Research in Neuropsychiatric Pathology, “Victor Babes” University of Medicine and Pharmacy, Timisoara, Romania
- 4Psychiatry Compartment, “Dr. Victor Popescu” Emergency Military Clinical Hospital, Timisoara, Romania
- 5“Eduard Pamfil” Psychiatry Clinic, “Pius Brinzeu” County Emergency Hospital, Timisoara, Romania
- 6International Centre for Education and Research in Neuropsychiatry (ICERN), Samara State Medical University, Samara, Russia
- 7Department of Psychiatry, Narcology, Psychotherapy and Clinical Psychology, Samara State Medical University, Samara, Russia
- 83rd Department of Psychiatry, School of Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
- 9Mental Health Section, Research Institute, Panhellenic Medical Association, Thessaloniki, Greece
Background: Suicidality is a serious public health concern at a global scale. Suicide itself is considered to be preventable death; worldwide, suicide rates and their trends are under constant scrutiny. As part of the international COMET-G cross-sectional study, we conducted a national level investigation to examine the individual disturbances (such as anxiety, depression, or history of life-threatening attempts) and contextual factors (such as adherence to conspiracy theories or Internet use) associated with suicidality related to the COVID-19 lockdown in a lot of Romanian adults.
Participants and Methods: One thousand four hundred and forty-six adults responded to an anonymous on-line questionnaire, with mean age ± standard deviation of 47.03 ± 14.21 years (1,142 females, 292 males, 12 identified themselves as non-binary). Data were analyzed using descriptive statistics and structural equation modeling (SEM).
Results: Univariate analysis showed strong significant correlation between anxiety and depression scorings among the respondents (Spearman R = 0.776, p < 0.001). Both the suicidality scorings and the Internet use correlated fairly with anxiety and depression, with two-by-two Spearman coefficients between R = 0.334 and R = 0.370 (p < 0.001 for each). SEM analysis substantiated the emotional disturbances, previous life-threatening attempts, and younger age as significant predictors for suicidality. The patterns of reality reading (including religious inquiries, Internet use, and beliefs in conspiracy theories) did not reach the statistical significance as influential factors in the suicidality of these respondents. There was no covariance between the Internet use and belief in conspiracy theories.
Conclusion: The study confirmed the suicidality risk initially hypothesized as being associated with the history of life-threatening attempts, increased depression within the younger population, and higher anxiety during the first year of the COVID-19 pandemic and its related lockdown. National strategies for effective interventions at various levels of the healthcare system should be developed.
Introduction
Suicidality is a serious public health concern at the global scale, affecting millions of people, their families, and society itself (1). The term “suicidality” includes suicidal ideation (SI, such as serious thoughts about taking one's own life), suicide plans, and suicide attempts (2). Significant resources and efforts have been focused on a better understanding of its underlying etiology, assessing the risks, and designing effective solutions at different levels of interventions (3). Suicide itself has been linked to the well-documented psychopathological risk predictors (such as suicidal behavior, history of self-harm and suicidal attempts), but there also are wide variations in suicidality indicators and suicide rates across countries and cultural environments (2, 4–13). Compared to other parts of the world, Europe is characterized by relatively high suicide rates, namely 10.5 (8.3–13.6) per 100,000 people per year (11, 12, 14). Depression has been acknowledged as a major risk factor for suicidality (15) and several studies have pointed toward anxiety as a major risk factor as well (16, 17). Demographic factors (e.g., young age, male gender, or ethnicity), social status (e.g., low income, income inequality, unemployment, low education, and low social support), social changes, neighborhood (e.g., inadequate housing, overcrowding, or violence), and adverse environmental events (e.g., climate change, natural catastrophe, war, conflict, and migration) were also linked to suicidality. Surprisingly, reported global trends for suicide rates and suicidal behavior demonstrated a stability not only before, but also in the early months of the COVID-19 pandemic (18).
On the other hand, the COVID-19 pandemic onset had a disruptive impact on societies, with global devastating consequences (19). As more and more countries instituted total lockdown, various reports pointed out that such measures exacerbated mental health issues, although the interventions were acknowledged as necessary and effective in stopping the spread of the virus (20–22). Some researchers focused on segments of the population at a higher risk, such as the youth or the frontline healthcare workers (23–26). Furthermore, the general population experienced exacerbated anxiety, with additional symptoms of depression, psychosis, panic attacks, trauma and suicidal ideation that seemed to exceed the experience in the previous SARS and MERS outbreaks (27). There were case reports of unusual neuropsychiatric manifestations like catatonia (28), but results regarding the rates of suicide behavior, attempts, ideation, and self-harm during the COVID-19 pandemic have varied and have been inconclusive (18, 29). The dramatic societal changes, serious environmental incidents, and a rise in family violence have been registered among the most influential factors highly correlated with the suicide risks during the period and after the COVID-19 total lockdown (30–32). Consequently, there were several warnings issued regarding the mental health in general (33) and suicidality in particular (34). In extraordinary times, such as the COVID-19 pandemic and the associated lockdown, suicide was rather unpredictable, with questionable dynamic rates. The suitable timeframe for assessing causal psychological changes and factors' inter-relationships arouse controversy over the gauging limitations, although memory-based retrospective assessment on behavioral and complex emotionality would offer the means to circumvent the distorting irrelevant momentary details and grant a respite for the emotions to settle and restructure (35).
Suicidality in Romania
In 2019, Romania reported an age-standardized suicide rate of 7.3 per 100,000 people per year, thus falling under the global age-standardized suicide rate of 10.5 per 100,000 people per year (14). Reported trends for suicide rates had been constantly decreasing since 2012, though with a consistent difference between sexes (i.e., females had a much lower rate than men) (36). To our knowledge, data on suicide risk factors in Romanian adult population has been scarce and of suboptimal quality.
On 16 March 2020, a state of national emergency was declared in Romania and total lockdown was instituted for 60 days, which brought a considerable burden of mental health consequences. A large community of migrant workforce in the Western Europe (over three million citizens), who massively returned home when the pandemic began, made Romania unique among the countries in European Union. Additional hurdles challenged the implementation of the protective measures: intrinsic weaknesses of the national healthcare system (e.g., aging infrastructure, low national health expenditure, and reported corruption), and one of the most religious populations in Europe (37, 38). Most Romanians identify themselves as Orthodox Christian, a highly conservative denomination, which was slow to react during this crisis (39). Notwithstanding these characteristics, psychological investigations in this period have reported the general population as being stable (40, 41), although actual information on suicide and suicidality is still too little.
Objective of the Study
In this paper we report the results of a national sub-set analysis comprised in the international COMET-G study (COVID-19 MEntal health inTernational for the General population) and based on the data from the Romanian population. In the pandemic context, the COMET-G study (22) aimed at investigating levels of depression, changes in anxiety, distress, suicidal ideation, and spreading of conspiracy theories in relation with a number of personal and interpersonal variables. Some national level findings have already been reported (21, 23, 42–45) along with the comprehensive report of the international study (22).
The specific target of this national level investigation was to examine the individual and contextual factors associated with suicidality in the Romanian adult population in the context of the COVID-19 pandemic related lockdown, which provoked major societal turmoil.
The main objective was to investigate the association of suicidality with individual proximal disturbances (such as anxiety and depression) and a history of life-threatening events. We also hypothesized the following secondary aims to be scrutinized: (a) contextual factors such as adherence to conspiracy theories propagated through the classical media and the Internet would play a significant role in suicidality; in addition, traditional cultural factors such as religiosity would also influence the individual pattern of reality reading and subsequent suicidal ideation; (b) socio-demographic factors (such as age and level of education) would play a role in suicidality.
Figure 1 illustrates the main objective and the secondary aims of this analysis. The conceptual framework included: suicidality, emotional disturbances, life threatening attempts and reality reading patterns. Suicidality refers to the “risk of suicide, usually indicated by suicidal ideation or intent, especially as evident in the presence of a well-elaborated suicidal plan” (46). Emotional disturbances comprise of three theoretical dimensions: emotional disturbances, emotion intensity/regulation disturbances, and emotion disconnections (47). Emotion intensity/regulation disturbances were mostly captured in the COMET-G study. Life threatening attempts encompassed the suicide attempts and the history of self-harm. Reality reading patterns would arise from the philosophical debate over the nature of conscious experience (48). The notion of indirect realism was extended to the reality perception in regard to the arising conspiracy theories, Internet use, and change in religiosity during the unfolding pandemic.
Figure 1. Main objective and secondary aims of this national level analysis of data collected on Romanian adult participants in the cross-sectional international COMET-G study.
Materials and Methods
Participants and Procedures
The study followed the cross-sectional COMET-G study protocol (22). The anonymous questionnaires (available in the Supplementary Material 1) gathered demographic data, general health data, previous psychiatric history, current symptoms of anxiety, depression and suicidality, and data regarding the changes caused by the lockdown in sleep patterns, sexual life, family relationships, finances, eating behavior, physical exercising, and religiousness/spirituality. Beliefs regarding the COVID-19 outbreak, perceived efficacy of the lockdown measures, and conspiracy theories were also investigated.
The international questionnaire was translated into Romanian according to established standards (49). Independent translation and back translation were conducted by two Romanian- English speaking authors. Following the Delphi technique, a panel of professionals agreed upon the final version that was deployed.
Retrospective data were collected online from 1 June to 23 December, 2020 (total lockdown had been instituted in Romania from 16 March until 15 May, 2020). Participants were instructed to give answers referring to their state and mindset during the total lockdown. No identification information was collected. Participants were able to access the survey and complete their responses only after reading and acknowledging the information regarding the study (i.e., the cover story): aim of the research, organizations involved and their contact information, and planned use of collected data. This acknowledgment served as the on-line form of informed consent. Announcement and advertisements were placed on social media, and distributed via e-mail and other instant messaging Apps.
Ethical approval (no. 194/ 4 June, 2020) was issued by the Ethics Committee of the “Pius Brinzeu” County Emergency Clinical Hospital, in Timisoara, Romania.
Instruments for Data Collection, Measures
Symptoms of anxiety were evaluated with State-Trait Anxiety Inventory (STAI), the S-Anxiety scale (STAI-Y1) (50). The STAI consists of 20 items that evaluate the respondent's current feelings on a 4-point Likert type scale. It is often employed for general and clinical populations (51) and had been used in Romania (52, 53).
Depression was evaluated with the Center for Epidemiological Studies-Depression Scale (CES-D), a popular and widely used instrument, based on self-reporting (54–56). It consists of 20 items that cover affective, psychological, and somatic symptoms (57), and had also been applied in Romanian population (56, 58–60).
Suicidality was evaluated with the Risk Assessment Suicidality Scale (RASS) (61), a self-assessment instrument. The last two RASS items were separately analyzed: RASS_11, “Have you ever hurt yourself in any way deliberately, during your whole life so far?”; RASS_12, “Have you ever attempted suicide, during your whole life so far?”. Each statement employed a 4-point Likert-type scale: 1 = not at all, 2 = a little bit, 3 = much, 4 = very much. An additional RASS-related item was included: “SI (Suicidal ideation) change,” “How much has your tendency to think about death and/or suicide changed, compared to before the outbreak of COVID-19?”. This instrument used a 5-point range: 2 = Very much increased, 1 = Increased a bit, 0 = Neither increased, nor decreased, −1 = Decreased a bit, −2 = Very much decreased. RASS had not been previously adapted for the Romanian population. Therefore, a confirmatory factor analysis (CFA) was conducted to verify its validity, based on the originally reported factors: “fear,” “intention,” and “life” (61).
Three additional Likert-type scales (designed for this study and included in the Supplementary Material 1) measured the extent of Internet use, belief in conspiracy theories, and the individual's religiosity. These scales underwent only face analysis prior to their deployment. CFA was conducted for the variables of Internet use and beliefs in conspiracy theories, which have been taken together as contributors to the patterns of reality reading.
Definition of Latent Variables Based on the Manifest Exogenous Variables
Suicidality (S) was inferred from the total score of the first 10 items of RASS (RASS tot) and from the change in suicidal ideation (SI change). Observable emotional disturbances (ED) were measured by the total score for STAI-Y1 (STAI tot) and total score for CES-D (CES tot). Life threatening attempts (LTA) was a latent variable based on the items of RASS_11 (RASS 11 self-harm) and RASS_12 (RASS 12 suicide). Reality reading patterns (RRP) were inferred from the change in individual's religiosity or spirituality inquiries (Relig increase), belief in conspiracy theories (Consp theories), and Internet use.
Data Analysis
Descriptive and Exploratory Statistics
Scale scores were treated as rank variables and described by the median (Inter Quartile Range). Descriptive statistics included the observed frequency counts (percent) for categorical variables or particular scales' selected items of interest. Normality of numerical variables was tested with the Shapiro-Wilk statistical test; these variables were described by the sample's mean and standard deviation (SD) when normally distributed, or by the sample's median (Inter Quartile Range–IQR) accordingly. The actual reliability of scale measurements was assessed based on the Cronbach's alpha: values >0.8 were considered to indicate good internal consistency, but scales with very few items were not discarded solely based on this coefficient. Cronbach's alpha actual values were reported for each scale with more than one item. The Harman's single factor method was applied to examine the amount of common method variance affecting the multi-item scales which had not undergone previous validation, other than face validation during the development stage. Harman's single factor method indicates possibly problematic common method bias (62). Separate application of the CFA marker technique to quantify the actual common method variance was unsuitable for scales taken in isolation, with possible additional issues related to the post-hoc choice of the marker (63). Non-parametric Spearman correlation approach was used to explore the covariances between various scales' scores employed in this study.
All reported probability values were two-tailed. A 0.05 level of significance was set, and highly significant values were also marked. Data were analyzed with the statistical software IBM SPSS v. 20.0 (Armonk, New York, USA) and the software packages R v. 4.0.5 (https://cran.r-project.org/).
Structural Equation Modeling and Confirmatory Factor Analysis
Based on the study specific target, structural equation modeling (SEM) was employed to investigate the structural connections between latent variables underlying the actual scores measured in the collected data. SEM was the method of choice for this analysis for its mathematical and statistical characteristics, i.e., a combination of model's structural features defined by equations, followed by their estimation across the available data based on the matrix algebra and generalized linear models. SEM is commonly used in the fields of social and psychological sciences for identifying hypothesized latent variables, which cannot be directly observed and measured. It also allows a simultaneous statistical estimation procedure, rather than separately estimating each part of a model, an approach which is believed to increase the overall accuracy (64).
We started with a nucleus model based on the main objective and its associated research hypotheses regarding the individual proximal disturbances (i.e., anxiety and depression) combined with a history of life-threatening attempts which would increase the suicidality related to the lockdown, thus including the endogenous latent variables of ED, LTA, and S. This model comprised previously validated scales as exogenous variables. In the following step, based on the secondary aims, we added the additional latent variable RRP in the model, which included the one-item change in individual's religiosity or spirituality inquiries, and the two multi-item scales for belief in conspiracy theories and Internet use (all three with only face validation). Furthermore, to this extended model we added two additional variables describing socio-demographic individual characteristics as potential independent predictors in the regression with S as an outcome. This approach yielded three SEM models, reflecting the results with reference to the main objective and the two secondary aims, respectively.
For all observed variables included in the models, the min-max rescaling was applied in order to preserve the shape of the original distributions and to retain the importance of outliers. The features would range as [0, 1] for all observed variables except for the change in suicidal thoughts, which was rescaled in the range [−1, 1] such that “no change” would correspond to a nil score. For model fitting, the maximum likelihood (ML) with robust estimators was used, with adjustments for non-normality of some variables (64, 65). The non-linear box-constrained optimization using PORT routines (NLMIB) was employed as the optimization method. When defining the SEM models, we placed the focus on the theoretical basis and meaningfulness of the variables' inter-relations. Nevertheless, the models were compared regarding their fit statistics and the Akaike information criterion (AIC). The Vuong's closeness test based on likelihood ratio was applied for determining the statistical significance of the change in AIC values.
The CFA and SEM models' goodness of fit indices and their corresponding [cut-off] values were: model Chi-square test and the resulting p-value, [<0.05]; Comparative Fit Index (CFI), [> 0.90]; Root Mean Square Error of Approximation (RMSEA), [<0.08 for a good fit, and up to 0.1 for marginal fit]; Standardized Root Mean Square Residual (SRMR), [<0.08].
The levels of statistical confidence and significance were 0.95 and 0.05, respectively, except for the RMSEA fit index, for which the confidence was explicitly specified to be 0.90. All reported probability values were two-tailed. We conducted the analysis with the statistical packages R v. 4.0.5 (including “lavaan” v. 0.6-9, “semPlot” v. 1.1.2, and “nonnest2” v. 2020-07-05).
Results
Confirmatory Factor Analysis and Harman's Test for the Scales at Their First Deployment in Romanian Population
CFA was conducted for the first 10 items of the RASS scale based on the three factors originally identified: fear, intention, and life. Results are presented in Table 1. All indices proved a good fit, except for the RMSEA, which was marginal. Figure 2 illustrates the CFA path diagram and factors' loadings, confirming the balanced contribution of all items to the overall score, with reversed effect for items #3 and #9. The actual RASS scale measurements were confirmed as consistent with scale's hypothesized construct.
Table 1. Confirmatory factor analysis (CFA) for the first 10 items of the RASS scale based on the three factors originally identified (fear, intention, and life) and the newly developed scales for belief in conspiracy theories (7 items) and Internet use (3 items).
Figure 2. The path diagram for the confirmatory factor analysis for the first 10 items of the RASS scale, based on the three factors originally identified: fear, intention, and life. Latent variables are drawn in circles and manifest variables are drawn in squares. The edge labels indicate the parameter estimates.
For the belief in conspiracy theories and Internet use, the Harman's single factor method resulted in 47.84 and 48.21%, respectively, of variance explained by one factor in exploratory factor analysis. These results on forced one factor model (namely less than 50% each) supported the further inclusion of the two scales in a SEM model. Table 1 also includes the CFA results for these two scales. Similarly to the 10-item RASS scale, the fit indices were good, except for the RMSEA, which was marginal.
Descriptive Analysis of Socio-Demographic Characteristics, Self-Reported Health Data and Mental Disturbances of the Respondents
One thousand, four hundred and forty-six (N = 1,446) adults responded to the anonymous questionnaire: 1,142 were females (aged 46.83 ± 14.16 years), 292 were males (aged 47.64 ± 14.36 years), and 12 self-identified as “non-binary” (aged 51.58 ± 15.45 years). Details of the respondents' socio-demographic information, and data regarding education and employment are presented in Tables 2A,B, respectively. Additional self-reported health related data were included in Supplementary Material 2.
Table 3 presents the descriptive statistics for the scales' scoring totals, and the corresponding values of Cronbach's alpha for individual scales or RASS sub-scales (i.e., as they resulted from the CFA). Lower values of internal consistency can be noted for the three-item “Internet use” (alpha = 0.456). For the three individual RASS items (O11, O12, and O13), the median (IQR) statistics were all nil. Table 4 shows the distribution of the scores for these individual items and also includes the scoring distribution for the change in religious/spiritual inquiries, where the actual spread over the whole range is apparent.
Table 3. Descriptive statistics and internal consistency of scale scorings for STAI, CES, RASS, belief in conspiracy theories, and internet use.
Table 4. The scorings' distributions for the individual items on suicidality change, personal history of self-harm, and increase in religious/spiritual inquiries.
Correlation between the scale scorings is presented in Table 5. The strong correlation between STAI and CES-D is noteworthy, although all scorings (except for the belief in conspiracy theories) were significantly two-by-two correlated. The belief in conspiracy theories showed a very weak or no relation to anxiety, depression, suicidality and Internet use.
Table 5. Associations between the scales total scorings on anxiety, depression, suicidality, conspiracy beliefs and Internet use.
Structural Equation Models
Table 6 shows the SEM models along with their statistical fit indices. We started with the nucleus Model 1, with reference to the main objective, which included three latent variables (ED, LTA, and S) and a regression (ED and LTA as predictors for S). In Model 2, we added an additional latent variable (RRP), which was included in the regression, as well. Model 3 kept the same latent variables, and also incorporated age and education as independent predictors in the regression. For all three models, we also investigated meaningful covariance. According to the Vuong's statistical test and the AIC, each model gave successively better description of the variables inter-relations, when compared to the previous one. For all three models, the fits indices reflected good reliability.
Table 6. The structural equation modeling of the multivariable relationships between the mental health indicators, beliefs and life changes.
Table 7 presents the parameters for the SEM model 3 in detail. ED and previous LTA were significant predictors for S, while the RRP were not. In addition, the participants' age was a significant predictor (with negative regression coefficient), but the level of education was not. It is important to note the significant covariance between each of the three latent variables considered as predictors in the regression, namely ED, LTA, and RRP; there was a significant negative covariance between the previous LTA (RASS_11 and RASS_12 items) and the reported change in SI during the COVID-19 pandemic lockdown. There was no covariance between Internet use and beliefs in conspiracy theories.
Table 7. The parameters of SEM Model 3 examining the relationships between anxiety, depression, life-threatening attempts, suicidality, religion/spirituality, conspiracy theories, Internet use scorings, age, and education.
Figure 3 shows the path diagrams for the SEM model 3. The latent variables are drawn in circles; the manifest variables are drawn in squares.
Figure 3. The path diagram for the SEM Model 3. Latent variables are drawn in circles and manifest variables are drawn in squares. The edge labels indicate the parameter estimates. CES tot, Center for Epidemiologic Studies Depression Scale (CES-D), 20-item total; ED, Emotional Disturbances; LTA, Life Threatening Attempts; RASS 11 self-harm, 4-point score of RASS 11; RASS 12 suicide, 4-point score of RASS 12; RASS tot, Risk Assessment Suicidality Scale (RASS), 10-item total; S, Suicidality; SI change, Suicidal Ideation change; STAI tot, State-Trait Anxiety Inventory (STAI), 20-item total.
The parameters of the SEM model 1 and model 2, and their corresponding path diagrams are presented in Supplementary Material 3.
Discussion
The present study on a lot of 1,446 Romanian adult participants in the international cross-sectional COMET-G study included persons aged between 19 and 84 years, with a mean of 47.03 years. More than 50% of the 1,446 respondents self-declared an increased level of religiosity and spiritual inquiries during the COVID-19 lockdown in the pandemic outbreak. Eighty-one percent self-reported no suicidality change, but more than 11% reported increased suicidal ideation during the lockdown. More than 10% of the 1,446 respondents admitted having a history of self-harm and more than 7% reported previous suicide attempts. In the structural models of suicidality, emotional disturbances and previous life-threatening attempts acted as significant predictors, while the patterns of reality reading were not. Young age was also a significant predictor for suicidality. The construct of suicidality was based on the RASS total scoring and the change in suicidal ideation. Of the two, the change in suicidal ideation played a more consistent role. For emotional disturbances, both STAI-Y1 total scoring (anxiety) and CES-D total scoring (depression) contributed in a similar way. Previous life-threatening attempts were observed in terms of two items of the RASS scale regarding the self-harm and previous suicide attempts, both contributing to life-threatening attempts in almost equal terms. The patterns of reality reading encompassed the adherence to conspiracy theories, Internet use, and change in spirituality inquiries. The Internet use had the highest estimate and the conspiracy beliefs the lowest, although both had high statistical significance for reality reading patterns. Although the SEM models 2 and 3 which included them were significantly better compared to the nucleus model, their contribution to the suicidality proved insignificant. They might only indirectly contribute through their significant covariance with the emotional disturbances.
We compared the socio-demographic characteristics for our responders with the officially reported data on the general population of Romania (66–69): median age of 43.2 years, rural residence of 43.6% in the general population (compared to 18.1% among the respondents), 61.1% married (59.5% in our data set), 4.84% unemployed (1.1% in our data set), 51.4 % females (78.9% in our data set). Summing up, compared to the general population of Romania, the respondents in the present study were of similar age, higher urban representation, similar marital status, higher employment status, and higher female representation. In particular, the dissimilarities in females' proportion and unemployment rates could have an impact on the models' validity, due to their previously reported effect on suicidality. Despite these concerns, the rate of Internet users in Romania is high and 12 million people use social media in Romania, a country with 19.18 million citizens (66, 70), therefore we confidently chose the on-line means to promote the COMET-G study. Although the respondents' sample was not totally representative for the general population, the number of respondents was high, compared to other countries cited in the COMET-G project (22).
Suicidality and self-harm history are widely acknowledged as substantial predictors for suicidal risk (71, 72), but the evidence is largely based on data from high-income countries (18). Data regarding the Romanian population is particularly scarce. In our models, previous life-threatening attempts proved to be significant predictors for suicidality. An intriguing finding was that both factors – suicide attempts and self-harm history – were negatively correlated with the change in suicidal ideation, albeit the correlation was weak (but statistically significant).
There are published reports of decreased suicidal ideation in association with the pandemic outbreak in Europe and the United States (18, 73, 74). Suicidal ideation might decrease when people are confronted with immediate potentially existential dangers, such as the risk of illness and the sense of incertitude during the COVID-19 pandemic. There are established observations of this phenomenon in other situations of immediate threats, like the First World War or terrorist attacks (75, 76). Depression was long seen as a suicide-related factor, but the effect of anxiety has not been separately investigated until recently (77, 78). The context of the COVID-19 pandemic might have also mediated a more direct connection between the increased anxiety and suicidality, as it unmasked and developed multiple anxiety-generating factors such as the fear of contamination, general insecurity, fear for the loved ones' health, and subject's overexposure on the media. What seems especially intriguing in this specific context is that the suicide rates were stable and suicidality was reported as decreasing; in a context when the general rates of the risk factors for suicidality (such as depression, anxiety, contextual and social vulnerabilities) increased in most of the reports on the COVID-19 pandemic and the contribution of these factors is well-established in the literature, alerts for constant vigilance regarding the suicidal dynamic were issued (18).
Age is usually inversely correlated with suicidal risk (4), and our SEM model 3 also put younger people at a higher risk. Studies on Romanian population showed a significant rise in suicide for young people and the elderly, even before the pandemic crisis (79). Lower education levels is typically seen as a general risk factor for suicidality, but it loses influence when adding other dominant factors, such as preexisting mental health issues, ancestry information, and demographic factors (7, 80, 81). In our SEM model 3, the level of education was not a significant predictor for suicidality. On the other hand, this lack of education significance in our model might be due to the disequilibrium in the level of education among the respondents: fewer than 30% of them did not graduate a form of post-high- school education. As these data were somewhat incongruous and unforeseen in the pandemic context, further investigations for long-term consequences in stationary societal circumstances are necessary, accompanied by national policies aimed at this public health issue.
There is a consensus that most people are quite resilient in face of negative changes or potentially traumatic events (82–84). Onset of societal or economic instability (for example, a recession) may have unstructured effects on suicide rates (85, 86). Nevertheless, vulnerability factors (such as previous mental health issues, suicide attempts, a history of self-harm, male sex, age, unemployment and belonging to disadvantaged social groups) may influence the life-long mental health risks, and indeed play an important role in the suicidality dynamics (77, 87–89). Religion generally plays a protective role regarding suicide (90), while religious turmoil is associated with a greater suicidal risk (91), albeit moderated by specific cultural differences (92). Responses to the COMET-G questionnaire showed that self-reported change in spiritual inquiries may have acted as a signal that previously successful coping mechanisms might have been exhausted, and the increased religiosity could thus be viewed as an attempt to regain emotional balance. Conspiracist ideation is also grounded on psychological mechanisms (93) and tends to increase during times of crises (94). Moreover, current media misinformation seems to generate a specific dynamic that exacerbates and promotes conspiracy thinking (95). These mechanisms were initially hypothesized to also work in the Romanian adult population, but in our structural models the conspiracy beliefs did not correlate with the degree of Internet use and did not demonstrate a significant influence on suicidality. This might be explained by the methodology we used that raises issues of consistency and common method bias for the variables related to the reality reading patterns. Both findings need further investigation and additional channels for proliferation of conspiracy theories should be considered. Specific scales, thoroughly validated, are also needed for clarification.
The pandemic context calls for consideration of new factors related to suicide. This implies that measures already established as being protective might require reconsideration and adjustment in the near future. For example, anxiety disorder and anxiety related distress emerged as a significant suicidality factor in the present study, thus needing deeper scrutiny in further research. We put forward a particular need for consolidation of the presently proposed structural models of suicidality.
Worldwide, several different studies have proposed vulnerability models for mental health issues (21, 42, 96) while suicide-related studies of the Romanian population have found it as atypical and therefore faced difficulties in applying models from other Eastern European Countries (97). Romania presents with a set of challenges regarding the medical system and with several cultural and socio-economic particularities, some of which are widely acknowledged as associating with higher suicidal risk. However, the national suicide rates have slightly declined over the past years, recently falling below the annual global age-standardized suicide rate (14). Precaution was recommended in regard to the pandemic consequences (98), but recent results showed a degree of psychological stability during the lockdown in Romanian population, and studies have indicated no change in suicide rates for some regions of Romania (99). However, little overall data is available, so the present results may bring valuable contribution toward moving forward with the novel understanding of suicidality.
Limitations
The main limitation of this investigation is that the proposed SEM models of suicidality were generated based on anonymous self-reported data, which were retrospectively collected in regard to the COVID-19 lockdown, within a limited time window and based on natural self-selection of respondents. Its cross-sectional design with no previous baseline and no follow-up prevented us from obtaining the risk estimates. Nevertheless, the SEM procedure in data analysis allowed the combination of the structural features with a general linear model for regression, and increased the overall accuracy and subsequent reliability of the findings.
The common method bias (CMB) implied by the cross-sectional design and the one time single-administration questionnaire (with its associated actual effect of the common method variance, CMV) is a major concern that cannot be overlooked. On the other hand, appropriate procedural measures were taken and carefully observed to limit the shared variance and control the method biases: different scales (such as those corresponding to predictors and criterion constructs) were included in non-adjacent sections, separated by questions collecting factual data (e.g., about diet or physical exercising); the scales included both positively and negatively (i.e., reverse) worded items; the wording was kept clear, concise and accurate; at the beginning of the questionnaire, respondents were provided comprehensive information on the COMET-G study and were assured of the anonymity; different scale formats were alternated, such as 4-point and 5-point Likert-type scales, or even dichotomy items.
In addition to these preventative measures, we explored the CMB possible impact on the performance of the measuring instruments and subsequent results by the post-hoc statistical techniques. Moreover, the approach with three SEM models (the nucleus including only previously validated scales and widely acknowledged constructs) and the stability of these nuclei regression coefficients' estimates (i.e., proximal predictors) across the three models proved the robustness of the results: significant and balanced interrelationship between the nuclei constructs (i.e., emotional disturbances and previous life-threatening attempts on the one hand, and suicidality on the other hand). We acknowledge that CMV, as a systematic error variance, could have a confounding influence on empirical results and produce potentially misleading conclusions, but this issue was improbable in our case. In our data set, there was a weak and insignificant relationship between the beliefs in conspiracy theories and the use of Internet – an actual CMB issue should have resulted in a stronger relationship.
Additional concerns might arise from employing scales not previously applied in the Romanian population, such as the RASS scale. Not only was the translation endorsed following a Delphi technique, but the CFA did support the original structure of the scale. The scales comprised in the construct of reality reading patterns proved to be less consistent and this issue should be addressed more carefully in the future. In addition, certain recall bias was possibly included in the answers.
Furthermore, this cross-sectional survey of self-reported perceived changes selectively recruited respondents who habitually navigate on the Internet, so the response rate was difficult to estimate, and acceptable rate was also problematic to anticipate or gauge. Moreover, this approach in questionnaire distribution led to a certain bias toward the population favorable toward on-line instruments, and this might have affected their appreciation toward the information and communication technology, thus the inconsistency on the scale of Internet use (it might have been too simplistic for many respondents).
The lack of follow-up imposes limits on the proposed models' external validity. An additional caveat regarding the validity originates in the pronounced gender disequilibrium among the respondents, which presumably reflect the degree of Internet engagement, but would affect the models' cross-gender validity for the rate of suicide completion is greater among males.
Conclusion
Suicidality has specific particularities for each country, region, or cultural context and environment, and the results we report bring evidence toward improving the insights into the Romanian population. Suicidality also has a context related inner dynamic, but affective disturbances, history of suicide attempts and self-harm remain the main factors related to suicide risk even in the special context of the COVID-19 pandemic in Romania. The final suicidality construct we developed also related to the COVID-19 lockdown, for this was specifically mentioned in the cover story of the questionnaire.
Our results confirmed anxiety and depression as significant proximal predictors in suicidality. In spite of every effort we made to answer the secondary aims of this study, the issue of quantifying the reality reading patterns' influence on suicidality remains open and must yet be further investigated.
Because suicide has a disastrous impact on the immediate family, it brings trans-generational mental health vulnerability. This investigation contributes to a better understanding of suicidality in a specific context, and may thus serve as a guide for assessing risks and identifying effective interventions.
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 human participants were reviewed and approved by Pius Brinzeu County Emergency Clinical Hospital, in Timisoara, Romania. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
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.
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.
Acknowledgments
We gratefully acknowledge the collective contribution of all Romanian participants who kindly agreed to respond to the COMET-G questionnaire. Authors express their gratitude to Prof. Paul Cumming of the Institute of Nuclear Medicine, Inselspital, Bern University, Bern, Switzerland, School of Psychology, Counseling and IHBI, Queensland University of Technology, Brisbane, Australia, and the International Centre for Education and Research in Neuropsychiatry, Samara State Medical University, Samara, Russia, for his valuable recommendations and review of the paper.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2022.818712/full#supplementary-material
References
1. WHO. Suicide: Key Facts. World Health Organization. (2021). Available online at: https://www.who.int/news-room/fact-sheets/detail/suicide/ (accessed January 19, 2021).
2. House A, Kapur N, Knipe D. Thinking about suicidal thinking. Lancet Psychiatry. (2020) 7:997–1000. doi: 10.1016/S2215-0366(20)30263-7
3. WHO. LIVE LIFE: An Implementation Guide for Suicide Prevention in Countries. Geneva: World Health Organization (2021). Available online at: https://www.who.int/publications/i/item/9789240026629 (accessed January 19, 2021).
4. Bachmann S. Epidemiology of suicide and the psychiatric perspective. Int J Environ Res Public Health. (2018) 15:1425. doi: 10.3390/ijerph15071425
5. Bauer M, Glenn T, Alda M, Andreassen OA, Angelopoulos E, Ardau R, et al. Association between solar insolation and a history of suicide attempts in bipolar I disorder. J Psychiatr Res. (2019) 113:1–9. doi: 10.1016/j.jpsychires.2019.03.001
6. Fountoulakis KN. Suicides in Greece before and during the period of austerity by sex and age group: relationship to unemployment and economic variables. J Affect Disord. (2020) 260:174–82. doi: 10.1016/j.jad.2019.09.001
7. Fountoulakis KN, Gonda X. Ancestry and different rates of suicide and homicide in European countries: a study with population-level data. J Affect Disord. (2018) 232:152–62. doi: 10.1016/j.jad.2018.02.030
8. Fountoulakis KN, Chatzikosta I, Pastiadis K, Zanis P, Kawohl W, Kerkhof AJ, et al. Relationship of suicide rates with climate and economic variables in Europe during 2000-2012. Ann Gen Psychiatry. (2016) 15:19. doi: 10.1186/s12991-016-0106-2
9. Fountoulakis KN, Savopoulos C, Zannis P, Apostolopoulou M, Fountoukidis I, Kakaletsis N, et al. Climate change but not unemployment explains the changing suicidality in Thessaloniki Greece (2000-2012). J Affect Disord. (2016) 193:331–8. doi: 10.1016/j.jad.2016.01.008
10. Fountoulakis KN. Suicide and the economic situation in Europe: are we experiencing the development of a 'reverse stigma'? Br J Psychiatry. (2016) 209:273–4. doi: 10.1192/bjp.bp.115.171785
11. Fountoulakis KN, Gonda X, Dome P, Theodorakis PN, Rihmer Z. Possible delayed effect of unemployment on suicidal rates: the case of Hungary. Ann Gen Psychiatry. (2014) 13:12. doi: 10.1186/1744-859X-13-12
12. Fountoulakis KN, Kawohl W, Theodorakis PN, Kerkhof AJ, Navickas A, Hoschl C, et al. Relationship of suicide rates to economic variables in Europe: 2000-2011. Br J Psychiatry. (2014) 205:486–96. doi: 10.1192/bjp.bp.114.147454
13. Turecki G, Brent DA. Suicide and suicidal behaviour. Lancet. (2016) 387:1227–39. doi: 10.1016/S0140-6736(15)00234-2
14. WHO. Age-Standardized Suicide Rates (per 100 000 Population), Both Sexes. (2019). Available online at: https://www.who.int/teams/mental-health-and-substance-use/suicide-data (accessed August 27, 2021).
15. Yoshimasu K, Kiyohara C, Miyashita K, The The Stress Research Group of the Japanese Society for Hygiene. Suicidal risk factors and completed suicide: meta-analyses based on psychological autopsy studies. Environ Health Prev Med. (2008) 13:243–56. doi: 10.1007/s12199-008-0037-x
16. Kanwar A, Malik S, Prokop LJ, Sim LA, Feldstein D, Wang Z, et al. The association between anxiety disorders and suicidal behaviors: a systematic review and meta-analysis. Depress Anxiety. (2013) 30:917–29. doi: 10.1002/da.22074
17. Khan A, Leventhal RM, Khan S, Brown WA. Suicide risk in patients with anxiety disorders: a meta-analysis of the FDA database. J Affect Disord. (2002) 68:183–90. doi: 10.1016/S0165-0327(01)00354-8
18. Pirkis J, John A, Shin S, DelPozo-Banos M, Arya V, Analuisa-Aguilar P, et al. Suicide trends in the early months of the COVID-19 pandemic: an interrupted time-series analysis of preliminary data from 21 countries. Lancet Psychiatry. (2021) 8:579–88. doi: 10.1016/S2215-0366(21)00091-2
19. WHO. Director-General's Opening Remarks at the Media Briefing on COVID19 -March 2020 (2020). Available online at: https://www.who.int/director-general/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19-−11-march-2020 (accessed August 27, 2021).
20. Wu T, Jia X, Shi H, Niu J, Yin X, Xie J, et al. Prevalence of mental health problems during the COVID-19 pandemic: a systematic review and meta-analysis. J Affect Disord. (2021) 281:91–8. doi: 10.1016/j.jad.2020.11.117
21. Fountoulakis KN, Apostolidou MK, Atsiova MB, Filippidou AK, Florou AK, Gousiou DS, et al. Self-reported changes in anxiety, depression and suicidality during the COVID-19 lockdown in Greece. J Affect Disord. (2021) 279:624–9. doi: 10.1016/j.jad.2020.10.061
22. Fountoulakis KN, Karakatsoulis G, Abraham S, Adorjan K, Ahmed HU, Alarcón RD, et al. Results of the COVID-19 mental health international for the general population (COMET-G) study. Eur Neuropsychopharmacol. (2022) 54:21–40. doi: 10.1016/j.euroneuro.2021.10.004
23. Patsali ME, Mousa DV, Papadopoulou EVK, Papadopoulou KKK, Kaparounaki CK, Diakogiannis I, et al. University students' changes in mental health status and determinants of behavior during the COVID-19 lockdown in Greece. Psychiatry Res. (2020) 292:113298. doi: 10.1016/j.psychres.2020.113298
24. Luo M, Guo L, Yu M, Jiang W, Wang H. The psychological and mental impact of coronavirus disease 2019 (COVID-19) on medical staff and general public – a systematic review and meta-analysis. Psychiatry Res. (2020) 291:113190. doi: 10.1016/j.psychres.2020.113190
25. Bredicean C, Tamasan SC, Lungeanu D, Giurgi-Oncu C, Stoica IP, Panfil AL, et al. Burnout toll on empathy would mediate the missing professional support in the COVID-19 outbreak. Risk Manag Healthc Policy. (2021) 14:2231–44. doi: 10.2147/RMHP.S300578
26. Kaparounaki CK, Patsali ME, Mousa DV, Papadopoulou EVK, Papadopoulou KKK, Fountoulakis KN. University students' mental health amidst the COVID-19 quarantine in Greece. Psychiatry Res. (2020) 290:113111. doi: 10.1016/j.psychres.2020.113111
27. Salari N, Hosseinian-Far A, Jalali R, Vaisi-Raygani A, Rasoulpoor S, Mohammadi M, et al. Prevalence of stress, anxiety, depression among the general population during the COVID-19 pandemic: a systematic review and meta-analysis. Global Health. (2020) 16:57. doi: 10.1186/s12992-020-00589-w
28. Kopishinskaia S, Cumming P, Karpukhina S, Velichko I, Raskulova G, Zheksembaeva N, et al. Association between COVID-19 and catatonia manifestation in two adolescents in Central Asia: incidental findings or cause for alarm? Asian J Psychiatr. (2021) 63:102761. doi: 10.1016/j.ajp.2021.102761
29. Dubé JP, Smith MM, Sherry SB, Hewitt PL, Stewart SH. Suicide behaviors during the COVID-19 pandemic: a meta-analysis of 54 studies. Psychiatry Res. (2021) 301:113998. doi: 10.1016/j.psychres.2021.113998
30. Seddighi H. COVID-19 as a natural disaster: focusing on exposure and vulnerability for response. Disaster Med Public Health Prep. (2020) 14:E42–3. doi: 10.1017/dmp.2020.279
31. Kelman I. COVID-19: what is the disaster? Soc Anthropol. (2020) 28:296–7. doi: 10.1111/1469-8676.12890
32. Prati G, Mancini A. The psychological impact of COVID-19 pandemic lockdowns: a review and meta-analysis of longitudinal studies and natural experiments. Psychol Med. (2021) 51:201–11. doi: 10.1017/S0033291721000015
33. Rajkumar RP. COVID-19 and mental health: a review of the existing literature. Asian J Psychiatr. (2020) 52:102066. doi: 10.1016/j.ajp.2020.102066
34. Banerjee D, Kosagisharaf JR, Sathyanarayana Rao TS. 'The dual pandemic' of suicide and COVID-19: a biopsychosocial narrative of risks and prevention. Psychiatry Res. (2021) 295:113577. doi: 10.1016/j.psychres.2020.113577
35. Murphy KR, Balzer WK. Systematic distortions in memory-based behavior ratings and performance evaluations: consequences for rating accuracy. J Appl Psychol. (1986) 71:39–44. doi: 10.1037/0021-9010.71.1.39
36. World Bank. Suicide Mortality Rate (per 100,000 Population) – Romania. (2019). Available online at: https://data.worldbank.org/indicator/SH.STA.SUIC.P5?locations=RO (accessed March 28, 2022).
37. Bratosin S, Tudor MA. Church, religion and belief: paradigms for understanding the political phenomenon in post-communist Romania. J Study Relig Ideol. (2009) 8:3–17. Available online at: http://www.jsri.ro/ojs/index.php/jsri/article/viewFile/262/261 (accessed April 10, 2022).
38. Tudor MA, Filimon Benea A, Bratosin S. COVID-19 pandemic lockdown and religious mediatization of social sustainability a case study of Romania. Int J Environ Res Public Health. (2021) 18:2287. doi: 10.3390/ijerph18052287
39. Dascalu S. The successes and failures of the initial covid-19 pandemic response in Romania, Front. Public Health. (2020) 8:344. doi: 10.3389/fpubh.2020.00344
40. Vancea F, Apostol MS. Changes in mental health during the COVID-19 crisis in Romania: A repeated cross-section study based on the measurement of subjective perceptions and experiences. Sci Prog. (2021) 104:368504211025873. doi: 10.1177/00368504211025873
41. MKOR Studiul impactul epidemiei de Coronavirus asupra vietii Românilor (in Romanian),. [Study impact of the Coronavirus epidemic on Romanian life 2020] (2020). Available online at: https://mkor.ro/studii/studiu-coronavirus-romania-2020/ (accessed June 10, 2020).
42. Šurina S, Martinsone K, Perepjolkina V, Kolesnikova J, Vainik U, RuŽa A, et al. Factors related to COVID-19 preventive behaviors: a Structural equation model. Front Psychol. (2021) 12:1821. doi: 10.3389/fpsyg.2021.676521
43. Vrublevska J, Sibalova A, Aleskere I, Rezgale B, Smirnova D, Fountoulakis KN, et al. Factors related to depression, distress, and self-reported changes in anxiety, depression, and suicidal thoughts during the COVID-19 state of emergency in Latvia. Nord J Psychiatry. (2021) 75:614–23. doi: 10.1080/08039488.2021.1919200
44. Kulig B, Erdelyi-Hamza B, Elek LP, Kovacs I, Smirnova D, Fountoulakis KN, et al. Effects of COVID-19 on psychological well-being, lifestyle and attitudes towards the origins of the pandemic in psychiatric patients and mentally healthy subjects: first Hungarian descriptive results from a large international online study. Neuropsychopharmacol Hung. (2020) 22:154-65.
45. Smirnova D, Syunyakov T, Pavlichenko A, Bragin D, Fedotov I, Filatova V, et al. Interactions between anxiety levels and life habits changes in general population during the pandemic lockdown: Decreased physical activity, falling asleep late and internet browsing about COVID-19 are risk factors for anxiety, whereas social media is not. Psychiatr Danub. (2021)33:119–29.
46. APA. American Psychological Association. Suicidality (in APA Dictionary of Psychology) (2022). Available online at: https://dictionary.apa.org/suicidality (accessed March 28, 2022).
47. Berenbaum H, Raghavan C, Le H-N, Vernon LL, Gomez JJ. A taxonomy of emotional disturbances. Clin Psychol. (2003) 10:206–26. doi: 10.1093/clipsy.bpg011
49. Flaherty JA, Gaviria FM, Pathak D, Mitchell T, Wintrob R, Richman JA, et al. Developing instruments for cross-cultural psychiatric research. J Nerv Ment Dis. (1988) 176:257–63. doi: 10.1097/00005053-198805000-00001
50. Spielberger CD, Gorsuch RL, Lushene RE. Manual for the State-Trait Anxiety Inventory (Self-Evaluation Questionnaire). Palo Alto, CA; Consulting Psychologists press (1970).
51. Fountoulakis KN, Papadopoulou M, Kleanthous S, Papadopoulou A, Bizeli V, Nimatoudis I, et al. Reliability and psychometric properties of the Greek translation of the State-Trait Anxiety Inventory form Y: preliminary data. Ann Gen Psychiatry. (2006) 5:2. doi: 10.1186/1744-859X-5-2
52. Miclea S, Albu M, Ciuca A. The romanian adaptation of endler multidimensional anxiety scale (EMAS). Cogn Brain Behav: Interdiscip J. (2009) 59–77. Available online at: https://www.academia.edu/ 16286637/The_Romanian_adaptation_of_Endler_Multidimensional_Anxiety _Scale_EMAS (accessed April 10, 2022).
53. Pitariu H, Peleaşă C. Manual for the Romanian Version of State-Trait Anxiety Inventory form Y. Synapsis. Cluj–Napoca, Romania. (2007).
54. Fountoulakis K, Iacovides A, Kleanthous S, Samolis S, Kaprinis SG, Sitzoglou K, et al. Reliability, validity and psychometric properties of the greek translation of the center for epidemiological studies-depression (CES-D) scale. BMC Psychiatry. (2001) 1:3. doi: 10.1186/1471-244X-1-3
55. Radloff LS. The CES-D scale: a self-report depression scale for research in the general population. Appl Psychol Meas. (1977) 1:385–401. doi: 10.1177/014662167700100306
56. Stevens MJ, Uscatescu L, Butucescu A, Constantinescu PM, Ion BC, Sandu CG. The Romanian CES-D scale: a promising instrument for clinical and non-clinical use. Int J Educ Psychol Commun. (2013) 3:7–21. Available online at: http://www.marianjournals.com/files/IJEPC_articles/Vol_3_no_1_2013/Stevens_and_al_IJEPC_3_1_2013.pdf (accessed April 10, 2022).
57. Carroll BJ, Fielding JM, Blashki TG. Depression rating scales: a critical review. Arch Gen Psychiatry. (1973) 28:361–6. doi: 10.1001/archpsyc.1973.01750330049009
58. Contantinescu PM, Glont ICN, Stevens M. The factorial structure of the romanian center for epidemiological studies depression scale: a single measurement of well-being and depression on a continuum. In; Milcu M, Stevens M, Neves de Jesus S, editors. Mind, Body and Behaviour: New Trends and Prospects in Health, Education and Social Sciences: 85-90. Editura Universitara. (2018).
59. Robila M. Effects of economic pressure on marital conflict in Romania. J Fam Psychol. (2005) 19:246–51. doi: 10.1037/0893-3200.19.2.246
60. Wallis AB, Fernandez R, Oprescu F, Chereches R, Zlati A, Dungy CI. Validation of a romanian scale to detect antenatal depression. Cent Eur J Med. (2012) 7:216–23. doi: 10.2478/s11536-011-0130-1
61. Fountoulakis KN, Pantoula E, Siamouli M, Moutou K, Gonda X, Rihmer Z, et al. Development of the risk assessment suicidality scale (RASS): a population-based study. J Affect Disord. (2012) 138:449–57. doi: 10.1016/j.jad.2011.12.045
62. Podsakoff PM, MacKenzie SB, Lee JY, Podsakoff NP. Common method biases in behavioral research: a critical review of the literature and recommended remedies. J Appl Psychol. (2003) 88:879–903. doi: 10.1037/0021-9010.88.5.879
63. Williams LJ, Hartman N, Cavazotte F. Method variance and marker variables: a review and comprehensive CFA marker technique. Organ Res Methods. (2010) 13:477–514. doi: 10.1177/1094428110366036
64. Kline RB. Principles and Practice of Structural Equation Modeling, Fourth Edition. New York, NY; Guilford Press ISBN. (2016).
65. Rosseel Y. Lavaan: an R package for structural equation modeling. J Stat Softw. (2012) 48:1–36. doi: 10.18637/jss.v048.i02
66. Worldometer. Romania Population. (2020). Available online at: https://www.worldometers.info/world-population/romania-population/ (accessed March 28, 2022).
67. Baetica R,. Nuptialitatea şi natalitatea în România (in Romanian). Eurostat Statistics Explained. (2015). Available online at: https://ec.europa.eu/eurostat/statistics-explained/index.php?title=Archive:Marriages_and_births_in_Romania/ro&oldid=255157 (accessed March 28, 2022).
68. O'Neill A,. Unemployment Rate in Romania 2020. (2020). Available online at: https://www.statista.com/statistics/373188/unemployment-rate-in-romania/ (accessed March 28, 2022).
69. World Bank. Population Female (% of Total Population) – Romania. (2019). Available online at: https://data.worldbank.org/indicator/SP.POP.TOTL.FE.ZS?locations=RO (accessed March 28, 2022).
70. DataReportal. Digital 2021: Romania. (2021). Available online at: https://datareportal.com/reports/digital-2021-romania (accessed March 28, 2022).
71. Victor SE, Klonsky ED. Correlates of suicide attempts among self-injurers: a meta-analysis. Clin Psychol Rev. (2014) 34:282–97. doi: 10.1016/j.cpr.2014.03.005
72. Guzmán EM, Cha CB, Ribeiro JD, Franklin JC. Suicide risk around the world: a meta-analysis of longitudinal studies. Soc Psychiatry Psychiatr Epidemiol. (2019) 54:1459–70. doi: 10.1007/s00127-019-01759-x
73. Cobo A, Porras-Segovia A, Pérez-Rodríguez M, Artés-Rodríguez A, Barrigón M, Courtet P, et al. Patients at high risk of suicide before and during a COVID-19 lockdown: Ecological momentary assessment study. B J Psych Open. (2021) 7:E82. doi: 10.1192/bjo.2021.43
74. Jacobson NC, Lekkas D, Price G, Heinz MV, Song M, O'Malley AJ, et al. Flattening the mental health curve: COVID-19 stay-at-home orders are associated with alterations in mental health search behavior in the United States. JMIR Ment Health. (2020) 7:e19347. doi: 10.2196/19347
75. Osman M, Parnell AC. Effect of the first world war on suicide rates in Ireland: an investigation of the 1864-1921 suicide trends. BJ Psych Open. (2015) 1:164–5. doi: 10.1192/bjpo.bp.115.000539
76. Claassen CA, Carmody T, Stewart SM, Bossarte RM, Larkin GL, Woodward WA, et al. Effect of 11 September 2001 terrorist attacks in the USA on suicide in areas surrounding the crash sites. Br J Psychiatry. (2010) 196:359–64. doi: 10.1192/bjp.bp.109.071928
77. Li Y, Li Y, Cao J. Factors associated with suicidal behaviors in mainland China: a meta-analysis. BMC Public Health. (2012) 12:524. doi: 10.1186/1471-2458-12-524
78. Bentley KH, Franklin JC, Ribeiro JD, Kleiman EM, Fox KR, Nock MK. Anxiety and its disorders as risk factors for suicidal thoughts and behaviors: a meta-analytic review. Clin Psychol Rev. (2016) 43:30–46. doi: 10.1016/j.cpr.2015.11.008
79. Dumitru M, Papari A, Seceleanu A, Sunda I. Suicide in Romania compared to the EU-28 countries. Soc Sci Edu Res Rev. (2019) 6:131–48. Available online at: https://econpapers.repec.org/article/edtjsserr/v_3a6_3ay_3a2019_3ai_3a2_3ap_3a131-148.htm (accessed April 10, 2022).
80. Huang X, Ribeiro JD, Musacchio KM, Franklin JC. Demographics as predictors of suicidal thoughts and behaviors: a meta-analysis. PLoS ONE. (2017) 12:e0180793. doi: 10.1371/journal.pone.0180793
81. Fountoulakis KN, Gonda X. Differential correlation of suicide and homicide rates according to geographical areas: a study with population-level data. Psychiatry Res. (2017) 249:167–71. doi: 10.1016/j.psychres.2016.12.059
82. Perkonigg A, Kessler RC, Storz S, Wittchen H.-U. Traumatic events and post-traumatic stress disorder in the community: prevalence,risk factors and comorbidity. Acta Psychiatr Scand. (2000) 101:46–59. doi: 10.1034/j.1600-0447.2000.101001046.x
83. Steel Z, Chey T, Silove D, Marnane C, Bryant RA, van Ommeren M. Association of torture and other potentially traumatic events with mental health outcomes among populations exposed to mass conflict and displacement: a systematic review and meta-analysis. JAMA. (2009) 302:537–49. doi: 10.1001/jama.2009.1132
84. Wu X, Kaminga AC, Dai W, Deng J, Wang Z, Pan X, et al. The prevalence of moderate-to-high posttraumatic growth: a systematic review and meta-analysis. J Affect Disord. (2019) 243:408–15. doi: 10.1016/j.jad.2018.09.023
85. Fountoulakis KN, Koupidis SA, Siamouli M, Grammatikopoulos IA, Theodorakis PN. Suicide, recession, and unemployment. Lancet. (2013) 381:721–2. doi: 10.1016/S0140-6736(13)60573-5
86. Ntountoulaki E, Paika V, Papaioannou D, Guthrie E, Kotsis K, Fountoulakis KN, et al. The relationship of the perceived impact of the current Greek recession with increased suicide risk is moderated by mental illness in patients with long-term conditions. J Psychosom Res. (2017) 96:98–105. doi: 10.1016/j.jpsychores.2017.03.008
87. Cao XL, Zhong BL, Xiang YT, Ungvari GS, Lai KYC, Chiu HFK, et al. Prevalence of suicidal ideation and suicide attempts in the general population of China: a meta-analysis. Int J Psychiatry Med. (2015) 49:296–308. doi: 10.1177/0091217415589306
88. Fountoulakis KN, Savopoulos C, Apostolopoulou M, Dampali R, Zaggelidou E, Karlafti E, et al. Rate of suicide and suicide attempts and their relationship to unemployment in Thessaloniki Greece (2000-2012). J Affect Disord. (2015) 174:131–6. doi: 10.1016/j.jad.2014.11.047
89. Xiang J, Kaminga AC, Wu XY, Lai Z, Yang J, Lian Y, et al. Lifetime prevalence of suicidal attempt among homeless individuals in North America: a meta-analysis. J Affect Disord. (2021) 287:341–9. doi: 10.1016/j.jad.2021.03.052
90. Wu A, Wang JY, Jia CX. Religion and completed suicide: a meta-analysis. PLoS ONE. (2015) 10:e0131715. doi: 10.1371/journal.pone.0131715
91. Paika V, Andreoulakis E, Ntountoulaki E, Papaioannou D, Kotsis K, Siafaka V, et al. The greek-orthodox version of the brief religious coping (B-RCOPE) instrument: psychometric properties in three samples and associations with mental disorders, suicidality, illness perceptions, and quality of life. Ann Gen Psychiatry. (2017) 16:13. doi: 10.1186/s12991-017-0136-4
92. Sisask M, Varnik A, Kolves K, Bertolote JM, Bolhari J, Botega NJ, et al. Is religiosity a protective factor against attempted suicide: a cross-cultural case-control study. Arch Suicide Res. (2010) 14:44–55. doi: 10.1080/13811110903479052
93. van Prooijen JW, van Vugt M. Conspiracy theories: evolved functions and psychological mechanisms. Perspect Psychol Sci. (2018) 13:770–88. doi: 10.4324/9781315525419
94. van Prooijen JW, Douglas KM. Conspiracy theories as part of history: the role of societal crisis situations. Mem Stud. (2017) 10:323–33. doi: 10.1177/1750698017701615
95. Stein RA, Ometa O, Pachtman Shetty S, Katz A, Popitiu MI, Brotherton R. Conspiracy theories in the era of COVID-19: a tale of two pandemics. Int J Clin Pract. (2021) 75:e13778. doi: 10.1111/ijcp.13778
96. Franklin JC, Ribeiro JD, Fox KR, Bentley KH, Kleiman EM, Huang X, et al. Risk factors for suicidal thoughts and behaviors: a meta-analysis of 50 years of research. Psychol Bull. (2017) 143:187–232. doi: 10.1037/bul0000084
97. Mäkinen IH. Eastern European transition and suicide mortality. Soc Sci Med. (2000) 51:1405–20. doi: 10.1016/S0277-9536(00)00105-2
99. Vuscan ME, Buciuta A, Vica ML, Balici S, Rusu SI, Siserman CV, et al. Impact of the COVID-19 pandemic on the suicidal behavior in Romania, Arch. Suicide Res. (2022). doi: 10.1080/13811118.2021.2022052. [Epub ahead of print].
Abbreviations and Nodes in Sem Models
• Age scale, age on a 6-group scale of adult age;
• CES tot, Center for Epidemiologic Studies Depression Scale (CES-D), 20-item total;
• Consp theories, 7-item total of belief in conspiracy theories;
• ED, Emotional Disturbances, latent variable describing self-reported fear, anxiety, depression;
• Edu scale, level of education on a 5-item scale;
• Internet use, 3-item total of Internet use;
• LTA, Life Threatening Attempts, latent variable describing suicide-related personal history;
• RASS 11 self-harm, 4-point score of RASS 11 (“Have you ever harmed yourself in any way deliberately, during your whole life so far?”);
• RASS 12 suicide, 4-point score of RASS 12 (“Have you ever attempted suicide, during your whole life so far?”);
• RASS tot, Risk Assessment Suicidality Scale (RASS), 10-item total;
• Relig increase, 4-point score of increase of religious/spiritual inquiries (“Over the last two-three weeks, my religious/spiritual inquiries have been increased.”);
• RRP, Reality Reading Patterns, latent variable describing personal vulnerabilities related to social and spiritual issues;
• S, Suicidality, latent variable describing suicidality;
• SI change, Suicidal Ideation change, 5-point score of change in suicidal ideation;
• STAI tot, State-Trait Anxiety Inventory (STAI), 20-item total.
Keywords: suicide, suicidal ideation, SEM, anxiety, depression, self-harm behavior
Citation: Panfil A-L, Lungeanu D, Tamasan S, Bredicean C, Papava I, Smirnova D and Fountoulakis KN (2022) Suicidality Related to the COVID-19 Lockdown in Romania: Structural Equation Modeling. Front. Psychiatry 13:818712. doi: 10.3389/fpsyt.2022.818712
Received: 19 November 2021; Accepted: 19 April 2022;
Published: 17 May 2022.
Edited by:
Umberto Volpe, Marche Polytechnic University, ItalyReviewed by:
Delia Virga, West University of Timişoara, RomaniaGratiela Noja, West University of Timişoara, Romania
Mikhail Sorokin, St. Petersburg V. M. Bekhterev Psychoneurological Research Institute, Russia
Copyright © 2022 Panfil, Lungeanu, Tamasan, Bredicean, Papava, Smirnova and Fountoulakis. 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: Diana Lungeanu, ZGx1bmdlYW51JiN4MDAwNDA7dW1mdC5ybw==
†These authors share senior authorship