Study Design
This study utilized publicly available data obtained from the National Health and Nutrition Examination Survey (NHANES), a comprehensive cross-sectional survey administered by physicians and highly trained medical personnel. The survey encompasses questionnaires, physical examinations, and laboratory data. NHANES aims to ascertain the prevalence and identify risk factors associated with major diseases in the U.S. population. Released biannually, the survey provides data collected from participants across the United States, selected through a sophisticated multistage, stratified sampling method. Each year, a representative sample of approximately 5000 noninstitutionalized individuals are surveyed across the US; data are released in 2-year cycles. This method has been demonstrated to yield a cohort that accurately represents the diverse demographic composition of the U.S. population[22, 23]. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines for cross-sectional studies.
Study Population
This was a cross-sectional analysis, using six 2-year NHANES cycles (2005–2016). By setting the inclusion criteria as Fig. 1 shown, a total of 16602 adults who had complete examination data on 9-item Patient Health Questionnaire (PHQ-9) questionnaire, sexual identity questionnaire, and other covariables, were included in this study.
Sexual Identity
Participants from 2005–2014 were categorized as heterosexual, gay/lesbian, bisexual, and others based on the following item: “Do you think of yourself as heterosexual or straight, homosexual or lesbian, bisexual, something else, or not sure?”. After the extensive cognitive and field testing, the sexual identity item was updated in 2015–2016 cycle[24]. The new item asked: “Which of the following best represents how you think of yourself?”. We excluded participants who responded “refused”, “don’t know”, or “I don’t know the answer” to the sexual identity item. The others category included participants who identified as something other than gay/lesbian, heterosexual, or bisexual. However, we are unable to determine which identities are represented in this category.
Depression Symptoms
The depression screener was composed of the PHQ-9 depression assessment tool. In our analyses, depression symptoms were defined as the summed score on the PHQ-9 was greater than or equal to 10 points (range, 0–27). The use of a PHQ-9 score cutoff of 10 points has been shown to offer a favorable balance between sensitivity and specificity[25].
Suicidal ideation was assessed by the following item: “Have you had thoughts that you would be better off dead or hurting yourself in some way?”. Then the participants responding “yes” for several days a week or more were subsequently considered as have a suicidal ideation in the further analysis[26].
Covariates
Based on the previous studies[27–29], the study considered the following covariates: age at the interview, sex (male, female), race/ethnicity (Mexican American, Other Hispanic, non-Hispanic white, non-Hispanic black, other race/multiracial), education levels (high school or below, greater than high school), family Poverty Income Ratio (PIR) (< 1.30, 1.31–3.50, ≥ 3.50), marital status (married, never married, living with a partner/ widowed/divorced/separated), body mass index (BMI) (< 25, 25–30, ≥ 30 kg/m2). Serum cotinine concentration was utilized as a proxy for environmental tobacco exposure and categorized into active/secondhand smoker (> 0.011 ng/mL) and nonsmoker (≤ 0.011 ng/mL). Alcohol drinking status was determined by the survey question, “In any year, have you had at least 12 drinks of any type of alcoholic beverage?”, with those responding “yes” subsequently considered as alcohol drinkers. We also considered the covariates related to medical history of hypertension (no or yes), and diabetes (no or yes). The drug use was defined as participants responded “yes” to the item: “Ever used marijuana or hashish?”.
Statistical analysis
Our data were analyzed following analytic guidelines and using the recommended survey weight for NHANES data[30]. All analyses in out study was sex stratified, and heterosexual participants of same sex were the reference group. To describe the characteristic of participants, we used mean (95% CI) for continuous variables and percentage frequency (95% CI) for categorical variables. Continuous data were compared using t-tests, and categorical data were compared using χ2 test. Because the percentage of missing data was small (missing rate range from 0–6.6%) for any variable, no imputation method was used.
Sex stratified multiple logistic regression models were performed to examine sexual identity differences in depression and suicidal ideation. Model 1 was unadjusted (crude model), Model 2 was adjusted for sociodemographic variables including age, sex, and race/ethnicity, education level, marital status, family PIR, Model 3 was fully adjusted model which including age, sex, race/ethnicity, education level, marital status, family PIR, BMI, alcohol drinker, serum cotinine, hypertension, diabetes, and drug use. A Bonferroni correction was used for multiple comparison.
All statistic analyses were performed with R (version 4.1.3, R Project for Statistical Computing, Vienna, Austria) and EmpowerStats (version 4.1, Boston, Massachusetts). In all tests, P < 0.05 (2-sided) was considered to indicate statistical significance.