Methods
Study population
This study used data from the 2007–2018 National Health and Nutrition Examination Surveys (NHANES), a cross-sectional survey conducted by the National Center for Health Statistics, part of the Centers for Disease Control and Prevention (CDC). NHANES uses a highly stratified multistage probability sampling design to assess the health and nutrition status of non-institutionalised US civilians. Descriptions of the survey protocol and sampling procedures can be found on the CDC website (https://wwwn.cdc.gov/nchs/nhanes/analyticguidelines.aspx%23).
Male and non-pregnant female participants ≥20 years of age with a BMI ≥18.5 kg/m2 were included in this study. Participants were required to complete the Depression Screener (DPQ) and Prescription Medications Questionnaire. Participants who reported using combinations of prescription psychotropic medications (ie, bupropion with naltrexone, fluoxetine with olanzapine, amitriptyline with perphenazine or amitriptyline with chlordiazepoxide) or more than one class of antidepressant were excluded from all analyses.
Exposure variables
Participants self-reported the frequency of their depressive symptoms over the past 2 weeks in response to items DPQ010 through DPQ090, the nine questions of the Patient Health Questionnaire-9 (PHQ-9).10 The PHQ-9 is a reliable and valid measure of depression, aligned with the diagnostic criteria for major depressive disorder in the Diagnostic and Statistical Manual of Mental Disorders-Fourth Edition.10 Each symptom was scored on a 4-point scale ranging from 0 (not experienced) to 3 (experienced nearly every day). The responses to each of the items were summed, with a score ≥10 indicating the presence of depression.10
Participants were asked if they had taken prescription medications within the past 30 days. Those who answered ‘yes’ were asked to either present their medication containers or, if unavailable, verbally report the name of the medication. The full list of generic names of antidepressants that were permitted for inclusion in this study, along with their classes, is provided in online supplemental table 1. Participants using ‘unspecified’ antidepressants were categorised as ‘yes’ for antidepressant use but were not included in analyses subset by class. Two participants who were using MAOIs were not included in antidepressant class analyses but were categorised as ‘yes’ for antidepressant use. Binary variables were constructed for each antidepressant class, where ‘yes’ indicated use of the specific antidepressant class, and ‘no’ indicated no use of the specific antidepressant class or no use of any antidepressant. Additionally, participants were asked to report their duration of antidepressant use, which was treated as a continuous variable measured in months.
Outcome variable
The outcome of this study was BMI (kg/m2). Weight and height were measured by trained personnel in the Mobile Examination Center (MEC) and were used to calculate BMI. Values were categorised into healthy (≥18.5–<25 kg/m2), overweight (≥25–<30 kg/m2) and obese (≥30 kg/m2).
Statistical analyses
All statistical analyses were performed using R V.4.2.2 and the package ‘survey’ to account for MEC survey weights. Survey weights were divided by six to account for the merging of six survey cycles. Continuous variables were presented as weighted mean (standard deviation, SD), while categorical variables were presented as weighted per cent. A χ2 test of independence was used to assess statistically significant (p<0.05) differences in categorical demographic characteristics across participants with and without depression, and a t-test was used to test for differences in continuous variables. Multivariable ordinal logistic regression was used to examine the relationship between BMI (healthy vs overweight or obese) and depression, antidepressant use (yes/no, antidepressant class and duration of use) or antidepressant use (yes/no) among those with depression along with covariates. Proportional odds assumptions were met for all models, which were run based on an available case analysis. P values were adjusted using Holm’s method in the antidepressant class analysis to account for multiple comparisons.
Potential covariates (ie, age, sex, race, poverty–income ratio, education, cigarette smoking status, alcohol use and minutes of sedentary activity per day) were selected based on prior work.9 Age was definitively included in the model.9 Univariate feature selection was then used to determine which of the remaining potential covariates were to be included in the final model if p<0.10. Final covariates for all models included age (continuous), sex (male, female), race (Mexican American, other Hispanic, non-Hispanic white, non-Hispanic black, other race (including multiracial)), education (≤high school, >high school), cigarette smoking status (not at all, some days or every day), alcohol use (never in lifetime or in past 12 months, at least once in past 12 months) and minutes of sedentary activity per day (continuous). In addition to selected covariates, inclusion of depression status, total PHQ-9 score or antidepressant use as a control was implemented when appropriate.