Suicide prevention is a priority of the Veterans Health Administration (VHA), which serves a patient population with elevated suicide rates compared with the general adult population of the United States (
1–
3). A screening tool that is quick to administer and provides valid information could help VHA clinicians to identify patients at greater risk of suicide and direct them toward appropriate treatment. Many suicide assessment measures have been developed, and their reliability and concurrent validity compared with previous suicide attempts or other psychiatric or mood scales have been documented (
4). However, the literature lacks prospective evidence regarding whether a screening instrument can predict suicide mortality across a range of settings (
5,
6).
The nine-item Patient Health Questionnaire (PHQ-9) is a reliable and validated measure of depression severity (
7). The final item of the PHQ-9 is commonly used in research studies to indicate the presence of suicidal ideation (
8–
10). It asks, “Over the past two weeks, how often have you been bothered by thoughts that you would be better off dead or of hurting yourself in some way?” Response options are “not at all,” “several days,” “more than half the days,” and “nearly every day.”
Suicidal ideation has been found to be a predictor of both immediate and long-term suicide risk (
11,
12). Studies have assessed the validity of item 9 of the PHQ-9 by comparing the results with those from a more detailed clinical interview (
13–
15). However, comparing item 9 to another self-reported assessment has limitations, given that both rely on a respondent’s willingness to report suicidal thoughts. To guide clinical care and enhance VHA suicide prevention, it is important to assess how suicidal thoughts expressed through responses to item 9 of the PHQ-9 may relate to risk of suicide mortality.
To our knowledge, only one study, by Simon and colleagues (
16), has evaluated whether item 9 of the PHQ-9 is associated with subsequent suicide mortality. The investigators analyzed a cohort of 84,418 outpatients ages 13 and older who were enrolled in an integrated health care system providing primary and mental health care. Providers in the system, with members in two states, had been instructed to conduct PHQ-9 assessments at each depression treatment encounter. The study documented positive associations between item 9 and risk of suicide death (hazard ratio [HR]=1.92, 95% confidence interval [CI]=1.53–2.41).
Replication studies are important in the scientific literature in order to establish an evidence base to guide clinical practice. As the largest integrated health system in the United States, the VHA maintains national data to evaluate and improve patient care. U.S. Department of Veterans Affairs (VA)/Department of Defense (DOD) clinical practice guidelines (CPGs) recommend annual assessment for depression in primary care and use of the PHQ-9 if patients screen positive for possible depression (
17). Among PHQ-9 recipients, information from item 9 may help VHA providers to identify those at elevated risk of suicide, enhance care, and prevent suicide deaths.
Methods
The VA’s Corporate Data Warehouse (CDW) was used to identify unique PHQ-9 assessments in fiscal year (FY) 2010 (October 1, 2009–September 30, 2010). Descriptive analyses were conducted at the patient level, and cohort construction and survival analyses were conducted at the assessment level. Assessment records were excluded if they were missing PHQ-9 score information or had improperly computed scores (N=2,609); were for patients whose records indicated age less than 18 years (N=864); had an assessment date after the patient’s date of death (N=7); were administered in the Philippines, American Samoa, or Guam (N=220) (the National Death Index [NDI] does not include deaths from these areas); or were administered to individuals who did not have recorded VHA outpatient or inpatient encounters in FY 2010 (N=79). Together, these exclusions resulted in removal of 3,779 (.8%) of the unique PHQ-9 assessments identified in FY 2010, with 447,245 assessments remaining for 391,492 individual patients. Patients excluded did not differ significantly from those who were retained with regard to suicide mortality. This study was conducted as part of ongoing VHA program evaluation and did not require informed consent and institutional review board review.
Suicide mortality was assessed from the time of each PHQ-9 assessment through September 30, 2011. Using previously reported procedures (
2), we identified vital status and cause of death data by NDI search records. Suicide death was defined by
ICD-10 codes X60–X84 and Y87.0 (
18).
Covariates were extracted from the CDW and VA national encounter data sets. At the patient level, risk factors collected were age, sex, and psychiatric diagnoses. Data on race-ethnicity were not collected because this information was not consistently available. Age was measured at the initial PHQ assessment date in FY 2010. The presence of psychiatric diagnoses within 365 days prior to this index assessment was assessed.
At the assessment level, six indicator variables were created to indicate the type of encounter for each PHQ-9 assessment: standard primary care (PC) encounters; PC–mental health integration (PC-MHI) encounters; specialty, non–PC-MHI mental health outpatient encounters; other outpatient encounters; mental health inpatient encounters; and other inpatient encounters not associated with mental health.
The sum of scores from the first eight items of the PHQ-9 was categorized into groups, consistent with procedures of Simon and colleagues (
16). Sums ranging from 0 to 4 were categorized as minimal depression; 5 to 9, mild depression; 10 to 14, moderate depression; and ≥15, severe depression. These groupings are similar to those used in other studies (
7).
Analyses were conducted by using SAS 9.3 and SAS Enterprise Guide 6.1. For descriptive analyses, frequencies were calculated for all variables at the patient level (N=391,492), overall and by suicide status. Chi square tests were used to test for differences in proportions of patients with various characteristics by suicide status. Suicide rates were calculated at the patient level by dividing the number of suicides by risk time, assessed as days from the initial assessment date until the end of the observation period or death (whichever came first for each individual), and dividing the quotient by 365. The result was multiplied by 100,000 to express rates per 100,000 person-years.
In proportional-hazards regression analyses, each individual PHQ-9 assessment was used as the unit of analysis. Risk time started on the day of the PHQ-9 assessment and ended at death or the end of the observation period (September 30, 2011), whichever came first. Consistent with methods of Simon and colleagues (
16), we used partly conditional proportional-hazards regression models (
19), controlling for nonsuicide mortality as a competing risk. By treating each assessment as the unit of analysis, the models were conditioned on the baseline covariates (age, sex, and psychiatric diagnoses) but not on the time-varying covariates (PHQ-9 score and type of encounter). Mirroring the models used in the analyses by Simon and colleagues (
16), the first four models assessed the association between item 9 of the PHQ-9 and suicide mortality while controlling for several covariates (age, sex, psychiatric diagnoses, and items 1–8 from the PHQ-9). A fifth model also adjusted for inpatient and outpatient encounters on each assessment date.
Covariance sandwich estimators were used to adjust for the nested nature of the data. The proportional-hazards assumption was evaluated by examining the significance of interaction terms between risk time and all covariates in regression models.
Two sensitivity analyses were conducted. First, we restricted the analyses to patients’ initial PHQ-9 in FY 2010 to examine whether eliminating repeated assessments would alter the observed associations between item 9 of the PHQ-9 and suicide mortality. Also, because we had assessed mental health diagnoses at an index date, and not the date of each assessment, we wanted to assess whether results differed when restricted to initial assessments. Second, replicating Simon and others (
16), to eliminate overlapping risk time for VHA users with more than one assessment, we ran all models with risk time censored at each user’s subsequent PHQ-9 assessment. These models produced results similar to those reported for the main analysis.
To assess potential differences in suicide risk in the study population, supplemental analyses compared suicide rates among FY 2010 VHA users who did or did not have a PHQ-9 assessment in that year.
Results
Most patients were male (92.7%) and age 50 or older (75.9%) (
Table 1). Sixteen percent had a diagnosis of depression within 365 days prior to their first PHQ-9 assessment, and 29.6% had moderate or severe depressive symptoms on the basis of items 1–8 from their first PHQ-9 assessment. On the date of their first PHQ-9 assessment, most patients (73.0%) had a standard primary care encounter, and 19.7% had a non-PC-MHI outpatient mental health or inpatient mental health encounter. A majority (87.6%) responded “not at all” to item 9 from their first PHQ-9 assessment.
Most patients (90.8%) had only a single documented PHQ-9 assessment in FY 2010. The maximum number of recorded assessments by a VHA user was 30. The number of risk days ranged from 1 to 730 (mean±SD=530.9±118.9). Among the 391,492 unique VHA users in the study, there were 310 suicide deaths, with a crude suicide mortality rate of 53.9 per 100,000 person-years (
Table 1). According to patients’ first PHQ-9 assessments, 7.4% of individuals who died by suicide in the follow-up period had reported that they were bothered nearly every day by thoughts that they would be better off dead or thoughts of hurting themselves (item 9), compared with 2.3% of patients who did not die by suicide (p<.001).
A total of 31 individuals died by suicide within 30 days of a PHQ-9 assessment (
Table 2); 18.2% of patients who responded “nearly every day” to item 9 from their most recent PHQ-9 assessment died by suicide within 30 days of the assessment, compared with 8.6% who responded “not at all.” In total, 71.6% of all suicide decedents responded “not at all” to item 9 from their most recent PHQ-9 assessment.
Without considering other covariates, analyses based only on item 9 of the PHQ-9 found that the relative hazard of suicide mortality increased with increased frequency of suicidal thoughts (
Table 3, model 1). A post hoc test for trend revealed a 64% increased risk of suicide for every unit increase in response to item 9 (HR=1.64, CI=1.47–1.82). [A figure showing the unadjusted suicide-specific cumulative incidence function by item-9 response is available as an
online supplement to this article.]
Although HRs were attenuated, a similar trend persisted after the analyses controlled for demographic characteristics, previous psychiatric diagnoses, and scores on PHQ-9 items 1–8 (HR=1.47, CI=1.29–1.68). Results were similar after the analyses also controlled for presence of inpatient and outpatient encounters on each assessment date.
Supplemental analyses indicated that the suicide rate for all FY 2010 VHA users was 53.3 per 100,000 person-years for users who completed PHQ-9 assessments, compared with 37.0 per 100,000 person-years for users who did not complete PHQ-9 assessments.
Discussion
Among VHA users who received PHQ-9 assessments in FY 2010, higher scores on item 9 were associated with increased risk of suicide mortality after analyses adjusted for demographic characteristics, psychiatric diagnoses, PHQ-9 items 1–8, and type of VHA encounter on the assessment date. These results were similar to analyses reported by Simon and colleagues (
16), who also documented a significantly elevated risk of suicide associated with higher scores on item 9 of the PHQ-9 while controlling for other risk factors. However, Simon and colleagues (
16) did not report a significant HR for responses of “nearly every day” compared with “not at all.” In this study, which had greater statistical power because of the large sample of assessments among VHA users, we observed a significant increased hazard for all item-9 responses compared with a response of “not at all,” after adjustment for other predictors.
The crude HRs representing risk of suicidality associated with item 9 of the PHQ-9 reported in this study were lower compared with HRs reported by Simon and colleagues (
16); however, their CIs overlapped. We noted differences in the population characteristics of the two samples. It is likely that the VHA patient sample was older and had more males compared with the sample analyzed by Simon and others (
16). Only 33.8% of the patients in this study had at least one psychiatric diagnosis, compared with 91% of patients in Simon and colleagues’ analyses who were receiving or had recently received mental health treatment. Also, only 18.3% of the total assessments in this study were categorized as indicating severe depression on the basis of items 1–8 of the PHQ-9 compared with 31% of assessments analyzed by Simon’s group (
16). Patient populations outside the VHA may use the PHQ-9 to primarily assess depression treatment, whereas VHA users who are not receiving depression-specific treatment may be asked to complete the PHQ-9 because of VA/DOD CPGs recommending use of this assessment in primary care (
17). Item 9 of the PHQ-9 assessment may better predict suicide among persons who already have mental health symptoms. After the analyses controlled for demographic characteristics, prior mental health care, and PHQ-9 items 1–8, the study results presented in this article were more similar to those reported by Simon and colleagues (
16).
Depression diagnosis was not found to be a significant predictor of suicide mortality after the analyses controlled for self-reported depressive severity (PHQ-9 items 1–8) (
Table 3, model 3). This suggests that patient reports of depressive-related symptoms were better predictors of suicide mortality compared with diagnosis. Similar to findings by Simon and colleagues (
16), self-reported depressive symptoms from items 1–8 appeared less significant after the analyses controlled for item 9 (
Table 3, models 4 and 5). Most important, even after the analyses controlled for psychiatric diagnoses, which are known risk factors for suicide (
20), as well as self-reported depressive severity (items 1–8), item 9 of the PHQ-9 was associated with suicide mortality in this population.
Although higher scores for item 9 were associated with increased suicide risk, the majority of suicides (71.6%) were by patients who responded “not at all” to item 9 on their most recent PHQ-9 assessment (
Table 2). These results were consistent with previous studies, suggesting that relying solely on item 9 of the PHQ-9 to assess suicide risk may miss a large portion of patients at risk of suicide mortality (
14).
Patients who died by suicide and did not endorse item 9 of the PHQ-9 may have chosen not to report suicidal ideation. A study designed to explore veterans’ perceptions about suicide screening found that veterans denied suicidal ideation for several reasons, including thoughts that suicidal ideation was shameful, belief that their thoughts were too private to share with strangers, and concerns regarding the consequences of divulging their suicidal thoughts (
21). Addressing these barriers could improve detection of suicidal ideation in clinic settings, enhancing the predictive utility of assessments like the PHQ-9.
We note that the study had limitations. First, although analyses examined all PHQ-9 assessments available in the CDW, these may not encompass all PHQ-9 assessments given to VHA users in FY 2010. Second, although our aim was to replicate the design used by Simon and colleagues (
16), we did not use the same covariates. Specifically, instead of controlling for prior psychiatric care, we controlled for prior psychiatric diagnoses. This may compromise the ability to compare results; however, psychiatric diagnoses are established suicide risk factors, and they indicate provider recognition of mental health treatment needs. Finally, the overall suicide rate was higher among VHA users who received PHQ-9 assessments compared with those who did not; therefore, findings may not be generalizable to the larger population of VHA users.