2.1. Design and participants
This study is a multicenter cross-sectional study conducted during the COVID-19. A total of 1236 questionnaires were retrieved from seven sperm banks in China. After screening according to the inclusion and exclusion criteria, 896 qualified questionnaires were finally obtained. The exclusion criteria were as follows: patients with genital tract disease, genetic disease diagnosis or family genetic disease, and infectious diseases; patients who are pathogen carriers; patients with hypertension, diabetes, and mental illness; and patients who lack semen quality reports. The inclusion criteria were as follows: age 18 to 42 years old, no history of mental illness, and no history of drug and organic disease in the last 3 months. The details of the study design and analysis of related processes are shown in Fig. 1. The study was approved by the Association of Ethical Review of Tongji Medical College, Huazhong University of Science and Technology. All participants completed questionnaires, including demographic characteristics, living habits, and psychological depression, and signed an informed consent. The semen concentration, semen volume, total sperm count, and forward motility were used as outcome indicators.
2.2. Depressive Symptoms
Patient Health Questionnaire
In this study, the Patient Health Questionnaire was used to evaluate the psychological depression status of subjects [19]. The scale consists of nine items, which belong to the criteria of major depression in the Diagnostic and Statistical Manual of Mental Disorders. This scale has high sensitivity and has been widely used in the adult population. Each item was scored on a scale of 0–3. Participants completed a mental scale assessment based on their psychological feelings according to their mental state in the last two weeks while completing the semen collection. On the basis of previous studies, a score of 5–9 was defined as mild depression, and a score of 10 or greater was defined as severe depressive symptoms [20, 21]. The reliability and validity of the simplified Chinese version of PHQ-9 have been demonstrated [22].
2.3. Semen Sample Collection And Analysis
The semen of each eligible subject was collected for routine semen analysis according to the inclusion criteria of the questionnaire survey, and the participants provided semen samples via masturbation in the sperm collection room. At the same time, the participant’s abstinence time and time of ejaculation were recorded. According to the latest guidelines of the World Health Organization, the semen concentration, semen volume, total number of sperm, and forward motility of sperm were detected and recorded for subsequent statistical analysis [23]. Immediately after ejaculation, the sample was collected and placed in an incubator at 37°C for incubation. The liquefaction time (< 1 h) was recorded. This step was completed within 1 h. The color of semen was observed and recorded, and the semen volume was estimated by weighing method (1 g = 1 mL). About 10 µL of semen was taken for smear to observe the morphology of sperm. The procedure was carried out in strict accordance with the latest WHO guidelines [24]. A computer-aided sperm analysis system (SCA CASA System, MicroChip S.L., Barcelona, Spain) was used to assist in the analysis of semen concentration and sperm motility indicators. In order to reduce errors caused by human manipulation, a skilled technician was selected for each sperm bank to perform the process in strict accordance with the implementation specifications.
The normal ranges for semen volume, total sperm count, semen concentration, and progressive motility were assessed strictly in accordance with the WHO fifth edition manual [25].
2.4. Statistical analysis
Classification variables in the data were presented in the form of frequency and percentage, while continuous variables were represented as mean ± standard deviation. The respondents were divided into depressed group and non-depressed group according to the scoring principle of PHQ-9. The scoring rules are described in the methods section. Table 1 describes the basic characteristics of the studying population, including weight, BMI, age, educational level, smoking, drinking, and other factors. Table 2 compares the differences in semen parameters (including semen concentration, total sperm count, semen volume, and progressive motility of sperm) between the two groups.
To verify the stability of the results, multivariate linear regression analysis and logistic regression analysis were used to examine the association between depression and sperm concentration, semen volume, total sperm count, and progressive motility. The odds ratio (OR) and 95% confidence interval (95% CI) were calculated, and the distribution of depression scale scores was significantly skewed to the right. Logarithmic transformation was performed before multiple linear regression analysis, and the estimates were reconstructed using exponential formulas [26–29].
To analyze the relationship between depression grade and semen quality, logistic regression model was used after converting the sperm concentration, semen volume, total sperm count, and progressive motility from continuous variable data to binary variables. The “WHO Human Semen Inspection and Processing Laboratory Manual” (fifth edition) was used to determine whether the indicators are qualified or not [7]. The values of semen quality indicators were directly incorporated into the linear regression model as continuous variable data, and depression symptoms were incorporated into the model as dichotomous variable data [30]. The model were adjusted for latent variables known to have an impact on semen quality (smoking[26], drinking[31], BMI[32] ) .When analyzing the association between depression and four outcome variables of semen quality, several indexes influencing semen quality reported by previous studies, including smoking, drinking and BMI, were set as fixed adjustment variables. Other covariables were filtered according to the screening principle of published articles. We selected these confounders on the basis of their associations with the outcomes of interest or a change of more than 10% in effect estimate [33, 34]. Supplementary Tables S1–S4 show the associations of each confounder with the outcomes of interest. The regression analysis results are shown in Tables 3 and 4.
The depression scale scores of all participants were divided into three categories according to the grading standard. Participants with scores between 0 and 4 were considered to have no depressive symptoms, those with scores between 5 and 9 were considered to be mildly depressed, and those with scores over 10 were considered to be severely depressed. The Cochran–Armitage trend test was used to determine whether there was a certain trend between depressive symptom grade and semen quality, as shown in Table 5. Curve fitting was used to test whether there was a linear relationship between the trend of depressive symptom score and various parameters of semen quality. The dose–response curve of each index is shown in Figs. 2–5. Meanwhile, the variables of each parameter index were adjusted according to the covariate screening principle. A two-stage linear regression model was used to test the threshold effect of depression score on semen quality by smoothing function.
For all analyses, the percentage of covariate’s missing values was less than 20%. We used multiple estimations to calculate the missing data of covariates, and used multiple regression algorithm to process the covariant data [35]. P-values < 0.05 were interpreted as statistically significant. All statistical analyses were performed using EmpowerStats (http://www.empowerstats .com) and the statistical package R(4.03 version).