Depression is Associated with Moderate-Intensity Physical Activity Among College Students During the COVID-19 Pandemic: Differs by Activity Level, Gender and Gender Role.

PURPOSE: The novel coronavirus disease (COVID-19) pandemic and associated restrictive measures have implications for depressive symptoms (henceforth depression) of young people and risk may be associated with their reduced physical activity (PA) level. Therefore, we aimed to examine the association between depressive symptoms and PA among college students with different gender and gender role (masculinity traits and femininity traits) during the COVID-19 pandemic. PARTICIPANTS AND METHODS: Cross-sectional study included 628 healthy college students from nineteen different locations. The Center for Epidemiological Studies Depression Scales (CES-D), the International Physical Activity Questionnaire - Short Form (IPAQ-SF), and the 50-item Chinese Sex-Role Inventory (CSRI-50) were used to measure depressive symptoms, PA continuous (weekly metabolic equivalent minutes, MET-minutes/week) and categorical indicators (activity level category) and gender role, respectively. The statistical analyses were used in partial correlation analysis, t-test, one-way ANOVA, moderation model tests, and linear regression model tests.
RESULTS: Total of 34.72% participants had clinically relevant depression (16, CES-D scale). Total of 58.6% participants were classified as a "low" activity level for spending less time on PA. Depression significantly negatively correlated with MET-minutes/week in moderate-intensity PA but not vigorous and walking scores. Of note, the depression-PA association was only moderated by the "low" activity level group in terms of categorical scores across gender groups. Participants with higher masculinity traits were less likely to have depression among all participants. Moreover, more recovered cases and fewer deaths could also predict the lower depression risk in the "high" activity level group.
CONCLUSION: Moderate-intensity PA is beneficial for reducing depression risk among college students at a low activity level. College students with fewer masculinity traits (regardless of gender) are highly vulnerable to depression during the outbreak of COVID-19. Effective control of the COVID-19 pandemic seems critical to alleviating the burden of mental disorders of the public including depression.

Introduction

The Novel Coronavirus disease (COVID-19) has formally been considered as a global pandemic on 11th March 2020.1 Given its high person-to-person transmission rate, multiple appropriate measures (eg, social distancing, quarantine, and isolation) were implemented across the world.2 Such major lifestyle changes during times of enormous uncertainty have been found to associate with mental health problems such as a higher risk of depression3–5 during the COVID-19 pandemic. Depression is a major mental disorder with symptoms including exaggerated and persistent sadness, anhedonia, cognitive impairments, and a negative procession bias as well as somatic symptoms (insomnia and fatigue).6 More than 322 million people worldwide are currently experiencing clinically relevant depression7 and would become one of the leading causes of disability worldwide by 2030.6,8 The risks for increasing levels of depression associated with the COVID-19 pandemic must, therefore, be considered as a major global health issue of increasing relevance.9,10

Compared with the prevalence of depression (3.6%) among young people (18–34 years) in a cross-sectional study conducted by the China Mental Health Survey (CMHS) in 2019,11 the risk of depression in Chinese young adults (22.1%) and college students (23.3%) have steeply risen during the COVID-19 pandemic.12,13 Increased levels of depressive symptomatology up to clinically relevant levels of depression have not only been found among patients and healthcare professionals (eg, surgical staff, medical doctor, and nurse) but also in the general population.14,15 In addition to healthcare professionals, young people were generally more likely to experience strongly increased levels of depression during the pandemic.13 Moreover, students who were affected by prolonged school closure may experience severe psychological stress,16 which may increase the likelihood of depression. Collectively, the currently available data suggest that college students are at an increased vulnerability to developing depression during the COVID-19 pandemic.9,12,17–19

Previous studies demonstrated that lower levels of depression were associated with more time spent on physical activity (PA) including team sports, gym exercises, and outdoor walking.20–22 Also, a recent large-scale meta-analysis of prospective cohort studies suggested that fewer depressive symptoms were associated with more frequent PA regardless of age and geographical region.23 Of note, these studies were not conducted during the COVID-19 pandemic when the public was recommended to isolate themselves and minimize public gatherings. In this context, it is obvious that people have to change their regular exercise patterns (eg, outdoor sports, gym exercises, and team sports). Although exercise scientists have discussed the detrimental impacts of the COVID-19 pandemic on PA,24–27 how those effects reflected in Chinese youth is still unclear.

Considering that more people might spend less time on PA during the COVID-19 pandemic, we further investigated the association of depression with PA across participants with different activity levels (low, moderate, and high) like previous studies did.23 On the other hand, the risk of depression was commonly observed to be higher in females than males,28 but the frequency of PA was generally lower in females than males,29 thus we considered gender as a potential moderator. Besides, accumulate studies suggested that psychosocial gender (gender role: masculinity and femininity traits) could be a sensitive indicator to interpret both mental and physical health among general populations.30,31 A well-known gender-role opinion suggested that everyone has two clusters of independent and orthogonal masculine (eg, dominant, has leadership abilities, and competitive, etc.) and feminine traits (eg, understanding, warm and compassionate, etc.), regardless of gender.32,33 Of note, people with greater orientation in masculinity and femininity were found experienced lower depressive symptoms than those with lower masculinity and femininity level.30 Besides, adolescent girls who increased their physical activity were also higher in masculinity and femininity than physically inactive girls with unhealthier mental functioning.34 In this regard, gender role was also considered as a potential moderator in the association of depression and PA.

Collectively, this study aimed to investigate associations between depression and PA among Chinese college students and to examine the effect of gender and gender role on such relationships during the COVID-19 pandemic. Results of this study would provide references for health professionals and governments to engage in PA-related intervention programs for college students with different risks of depression during the pandemic.

Methods

This study used a cross-sectional design. Data were collected between the 10th to 16th, March in 2020. At that time, the lockdown policies have been implemented in almost all regions of China as a measure to regulate infection rates.35 Of note, all colleges and universities in China closed at that time. All participants took part in this study voluntarily with a written informed online consent form. Each participant was remunerated with CNY10.0 (approximately 1.5 US dollars) upon the completion of their participation. This study was approved by the Human Research Ethics Committee (No: 2,020,005) of Shenzhen University in accordance with the Declaration of Helsinki.

Participants

A total of 628 healthy Chinese college students (64.8% females, mean aged (±SD) 20.17 (±1.87) years) without any symptoms related to the COVID-19 participated in this survey distributed through social media (WeChat) app (), which is one of the most popular online social platforms in China connecting a billion people with free calls and chats. People can access to the survey links from their individual or group chats and advertisement pages of this platform with convenience. All participants (freshmen = 309, sophomore = 143, junior = 75, senior = 86, and post-graduate = 15) lived in nineteen different locations across China. Three participants who missed their location information were excluded, leading a final sample size of 625 participants for statistical analyses.

Measures

Physical Activity

The International Physical Activity Questionnaire-Short Form (IPAQ-SF) is a widely used instrument developed by the International Consensus Group on Physical Activity Measurement to measure physical activity in adults.36,37 The original questionnaire was translated and validated in Chinese.38,39 The IPAQ-SF consists of 6 items about varies dimensions of PA include walking, moderate-intensity (eg, carrying light loads, bicycling at a regular pace) and vigorous-intensity (eg, heavy lifting, digging, aerobics) activity in terms of both the duration (in minutes/day) and frequency (days/week) and one item about the duration (in minutes/day) of sedentary. According to the guidelines for data processing and analysis,40 both continuous (the Metabolic Equivalent of Task-minutes/week (MET-m/w) for walking, moderate-, vigorous-intensity, and sitting activity) and categorical indicators (low, moderate, and high activity level) of PA are possibly obtained from IPAQ-SF after a raw data processing procedure (see ) recommended by the IPAQ group. For the continuous indicators, scores in each activity could be processed into the MET-m/w using MET values for sitting (1.0 MET), walking (3.3 METs), moderate-intensity (4.0 METs), and vigorous-intensity (8.0 METs). For example, if someone spent 45 minutes/day for moderate-intensity PA during last 7 days, his moderate-intensity PA MET-m/w would be 1039.5 via a formula (MET-m/w = MET-value * minutes * days). For the categorical indicators, participants can be categorized into “low” activity level (LAL), “moderate” activity level (MAL), or “high” activity level (HAL) categories using the published criterion (see ) recommended by the IPAQ group. In the current study, both MET-m/w and activity levels were regarded as PA indices.

Depressive Symptoms

Depressive symptoms were measured using the self-reported Center for Epidemiological Studies Depression Scales (CES-D) in Chinese with good reliability and validity of Cronbach’s alpha coefficients of 0.90 for full scale.41,42 It consists of 20 items within four dimensions (depressed affect, positive affect, somatic complaints, and interpersonal problems).43 Each item can be scored on a 4-point Likert scale (0 = rarely or not at all to 3 = most of the time). In this regard, individuals can obtain a maximum of 60 points, with higher scores indicating more severe depression. In the present study, the Cronbach’s alpha coefficients were 0.95, 0.88, 0.84, 0.83, and 0.79 for the full scale of CES-D, depressed affect, positive affect, somatic complaints and interpersonal problems subscales, respectively. The continuous scores in CES-D scale were adopted as depression indices.

Gender Role

Gender-role traits were measured using the 50-item Chinese Sex-Role Inventory (CSRI-50).44 It included three dimensions related to masculinity (eg, leadership, rationality, and generosity, etc.), femininity (eg, empathy, thrift, and careful, etc.), and neutral traits. Each item is scored from 1 (never or almost never true) to 7 (almost always or always true), with higher scores in masculinity or femininity indicating more masculine or feminine traits. Additionally, continuous gender-role scores can be classified into four gender-role identities (masculine = high in masculinity and low in femininity, feminine = high in femininity and low in masculinity, androgynous = both high in masculinity and femininity, and undifferentiated = both low in masculinity and femininity) based on the median of distributions of masculinity and femininity scales. In the present study, Cronbach’s alpha coefficients were 0.91 and 0.88 for items on masculinity and femininity subscales, respectively. The recommended median score of masculinity (4.8) and femininity (5.0) subscales were adopted in this study.44 Also, both gender-role continuous scores and categories were regarded as gender-role indices.

Demographic Variables and COVID-19 Data

Demographic data were also collected, including gender, age, grade, and location of current residence. In addition, from the official website of the National Health Commission in China,45,46 we had extracted data related to the COVID-19, such as the number of new cases, confirmed cases, recovered patients, and death toll, which was matched with location of participants and the date when they filled out the questionnaire.

Statistical Analysis

Step 1: Pearson correlation analyses were performed to investigate associations of depression (sum scores of the full scale) with MET-m/w in different levels (walking, moderate-intensity, and vigorous-intensity) of PA and sedentary activity. Significant associations of depression and PA were further investigated in terms of four dimensions of the CES-D when controlling for age, grade, sedentary MET-m/w, and the COVID-19-related data using partial correlation analyses.

Step 2: Between-group differences on PA MET-m/w (walking, moderate-intensity, vigorous-intensity, and sedentary activity) and depression were tested across gender (female and male), and gender-role identities (undifferentiated, androgynous, masculine, and feminine identity). Likewise, we further tested between-group differences in depression across activity levels (LAL, MAL, and HAL). This is followed by moderation analyses where the CES-D score, MET-m/w (walking, moderate-intensity, vigorous-intensity, and sedentary activity), and potential confounders (activity level, gender, and gender-role continuous scores) were considered as the dependent variable, independent variables, and mediating variables, respectively. Gender differences in PA MET-m/w and depression were conducted using independent sample t-test, respectively. Between-group comparisons regarding activity level and gender-role identities categories in PA MET-m/w and/or depression were conducted using one-way ANOVA. The moderation analyses were conducted using the PROCESS by model 1.

Step 3: To better understand whether and how depression (sum scores of full scale) were impacted by PA MET-m/w (walking, moderate-intensity, vigorous-intensity, and sedentary activity), age, gender, gender-role continuous scores, grade, the COVID-19-related data among all participants and across different activity levels (LAL, MAL, and HAL), we, therefore, used linear regression to test these models, by the stepwise method. If the significant results were found in One-way ANOVA, the post-hoc comparison would be conducted using Bonferroni-correction. All statistical analyses were conducted using SPSS v22.0 with PROCESS v3.3.47 A significance threshold of p < 0.05 was set when a two-tailed test was used.

Results

Characteristics of Participants

Of 625 participants (mean aged 20.18 years, 64.8% female), 217 (34.72%) reported depression scores over the cut-off points (≥16, CES-D). The median (percentiles 25 to 75) MET-m/w in vigorous-, moderate-intensity, walking and sedentary were 0 (0 to 240), 0 (0 to 450), 99 (0 to 462) and 2940 (2100 to 3360), respectively. Fifty-eight (9.28%) participants were included in the high activity level (HAL) group and 201 (32.16%) were included in the moderate activity level (MAL) group, while the remaining 368 (58.6%) participants were included in low activity level (LAL) group. For gender-role identities, 211 (33.76%), 92 (14.72%), 80 (12.8%) and 242 (38.72%) participants were included in undifferentiated, masculine, feminine and androgynous identities, respectively. The characteristics of the participants are reported in Table 1.

Table 1

Characteristics of Participants

Continuous Variables		n	Mean	SD	Median	Range(Percentiles 25 to 75)
Age
		625	20.18	1.88	/	/
IPAQ-SF MET-minutes/week
	Vigorous PA	625	399.65	979.84	0	0 to 240
	Moderate PA	625	318.56	527.22	0	0 to 450
	Walking	625	354.39	605.16	99	0 to 462
	Sedentary	625	2929.14	1286.61	2940	2100 to 3360
	Moderate-to-Vigorous PA	625	718.21	1321.03	/	/
IPAQ-SF frequency (days/week, last 7 days)
	Vigorous PA	625	1.03	1.92	/	/
	Moderate PA	625	1.89	2.35	/	/
	Walking	625	2.35	2.65	/	/
IPAQ-SF duration (minutes/day)
	Vigorous PA	625	12.95	31.09	/	/
	Moderate PA	625	20.87	30.27	/	/
	Walking	625	24.09	33.94	/	/
	Sedentary	625	418.45	183.80	/	/
	Moderate-to-Vigorous PA	625	34.82	50.82	/	/
Depressive symptoms
	CES-D score	625	13.48	12.26	/	/
	Depressed affect	625	5.56	5.20	/	/
	Positive affect	625	2.63	2.89	/	/
	Somatic complaints	625	4.48	3.97	/	/
	Interpersonal problems	625	0.8	1.316	/	/
Gender-role scores
	Masculinity scores	625	4.80	0.85	/	/
	Femininity scores	625	4.96	0.78	/	/
COVID-19 data (province level)
	New cases	625	0.51	1.24	/	/
	Accumulated confirmed cases	625	1292.63	243.92	/	/
	Accumulated recovery cases	625	1227.01	233.50	/	/
	Accumulated death cases	625	7.70	1.29	/	/
Categorical Variables		n	%
PA level categories
	Low level	366	58.56
	Moderate level	201	32.16
	High level	58	9.28
Gender
	Male	220	35.20
	Female	405	64.80
Gender-Role identities
	Undifferentiated	211	33.76
	Masculine	92	14.72
	Feminine	80	12.80
	Androgynous	242	38.72
Grade
	Grade1	308	49.28
	Grade2	142	22.72
	Grade3	74	11.84
	Grade4	86	13.76
	Graduate students	15	2.40
Depression
	Undepressed (<16)	408	65.28
	Depressed (≥16)	217	34.72

Risk of Depression Associated with PA

As shown in Table 2, depressive symptoms were significantly negatively related to MET-m/w in both vigorous- and moderate-intensity PA, whereas no significant associations of depression with MET-m/w in walking and sedentary were observed.

Table 2

Correlations Analyses of CES-D Depressive Symptoms Scores and IPAQ MET-Minutes/Week

		Vigorous PA MET-Minutes/Week	Moderate PA MET-Minutes/Week	Walking MET-Minutes/Week	Sedentary MET-Minutes/Week
Pearson Correlations (No Covariate)
All participants
	CES-D Score	−.08*	−.11**	−.060	0.059
Partial Correlations (Age, Educational Attainment, Sedentary MET-m/w, and COVID-19 Data as Covariates)
All participants
	CES-D Score	−.066	−.10*	−.048
	Depressed Affect Score	−.059	−.09*	−.054
	Positive Affect Score	−.067	−.10*	−.067
	Somatic Score	−.062	−.09*	−.037
	Interpersonal Score	−.047	−.09*	0.022
Partial Correlations(Sex, Gender-Role Scores, Age, Educational Attainment, Sedentary MET-m/w, and COVID-19 Data as Covariates)
Low activity level
	CES-D Score	−0.04	−0.15*	−0.04
	Depressed Affect Score	−0.04	−0.15*	−0.04
	Positive Affect Score	−0.05	−0.14*	−0.06
	Somatic Score	−0.01	−0.11*	−0.02
	Interpersonal Score	−0.03	−0.11*	−0.01
Moderate activity level
	CES-D Score	0.08	−0.03	0.05
	Depressed Affect Score	0.06	0.00	0.02
	Positive Affect Score	0.15	−0.04	0.05
	Somatic Score	0.03	−0.03	0.06
	Interpersonal Score	0.11	−0.06	0.14
High activity level
	CES-D Score	−0.25	−0.01	0.10
	Depressed Affect Score	−0.17	0.01	0.05
	Positive Affect Score	−0.24	0.06	0.03
	Somatic Score	−0.28	−0.03	0.14
	Interpersonal Score	−0.08	−0.13	0.15
Partial Correlations(Gender-Role Scores, Age, Educational Attainment, Sedentary MET-m/w, and COVID-19 Data as Covariates)
Female
	CES-D Score	−0.01	−0.05	0.02
	Depressed Affect Score	0.01	−0.02	0.01
	Positive Affect Score	−0.04	−0.06	−0.02
	Somatic Score	−0.02	−0.05	0.05
	Interpersonal Score	−0.01	−0.09	0.03
Male
	CES-D Score	−0.09	−0.14*	−0.05
	Depressed Affect Score	−0.12	−0.17*	−0.07
	Positive Affect Score	−0.05	−0.11	−0.04
	Somatic Score	−0.07	−0.12	−0.06
	Interpersonal Score	−0.07	−0.09	0.07

Notes: **Correlation is significant at the 0.01 level (2-tailed). *Correlation is significant at the 0.05 level (2-tailed).

After controlling for demographic variables (age and grade), MET-m/w in sedentary activity, and data related to the COVID-19 data, results indicated that more time spent on moderate-intensity PA MET-m/w was significantly associated with a lower risk of developing depressive symptoms.

Differences by Activity Level

Considering that most of the participants reported a lower level of PA in this study, we further compared between-group differences of activity levels on depressive symptoms and a significant difference was found (F(2, 622) = 4.32, p < 0 0.05, ηp2 = 0.014), indicated that LAL group reported a significantly higher level of depression than the HAL group (Cohen’s d = 0.39, 95%C.I. 0.11 to 0.67). As showed in Table 2, significant associations of depression (sum scores and subscale scores related to four dimensions) and MET-m/w of moderate-intensity PA were only observed in the LAL group, but not in the MAL and HAL groups, after controlling for age, gender, gender role, grade, sedentary MET-m/w and COVID-19 data. Furthermore, participants in the LAL group who had a greater amount of MET-m/w in moderate-intensity PA were associated with less depressive symptoms. The moderation effects analysis showed that when LAL group was coded by one and others by zero, the moderation model was significant (F (5, 619) = 3.69, p < 0.01, model-R = 0.17, model-R2 = 0.03).

Differs by Gender and Gender Role

No significant gender differences with respect to depression, MET-m/w of walking and sedentary activity was found, but significant differences according to gender were observed for MET-m/w of moderate-intensity (t (390.97) = −3.3, p < 0.01) and vigorous-intensity PA (t (398.89) = −2.74, p < 0.01). Specifically, females had spent less time on moderate- and vigorous-intensity PA than males. As shown in Table 2, there were significant negative associations between depression and MET-m/w of moderate-intensity PA in males but not in the females. Males with higher frequency in moderate-intensity PA reported a lower risk of developing depressive mood symptoms (depressed and positive affect). No significant moderation model was found when using gender as the moderator.

For the gender role, significant between-group differences of four identities (undifferentiated, masculine, feminine, and androgynous identities) on depression (F (3, 621) = 18.08, p < 0.001, ηp2 = 0.08) and MET-m/w of walking (F (3, 621) = 5.46, p < 0.01, ηp2 =0 0.026) were found. For depression, individuals with undifferentiated (17.28±12.98) and feminine (16.9±12.6) identities reported a significantly higher risk of depression than those endorsed with the masculine (10.95±10.57) and androgynous (9.98±10.82) identities. For MET-m/w of activities, the individuals with androgynous identity spent more time walking than those endorsed with undifferentiated identity (Cohen’s d = 0.35, 95%C.I. 0.17 to 0.54). However, no significant moderation model was found when using gender role as the moderator.

Linear Regression Models for Depression

Considering that activity level could directly moderate the associations of depression and MET-m/w, the linear regression model using stepwise method across activity levels were conducted after a similar analysis regarding all participants (Table 3).

Table 3

Linear Regression Models for Depression

PA Level Categories			Unstandardized Coefficients	95.0% C.I. for B		Standardized Coefficients	p-value	Collinearity Statistics		Adjusted R2	Durbin-Watson	Model p-value
			B	Lower	Upper	β		Tolerance	VIF
All participants
	Model 1	Masculinity score	−4.07	−5.15	−2.98	−.28	<0.001	1.00	1.00	0.08	2.00	<0.001
	Model 2	Masculinity score	−4.30	−5.40	−3.2	−0.30	<0.001	0.97	1.03	0.09		<0.001
		Sex	2.40	0.44	4.36	0.93	0.016	0.97	1.03
	Model 3	Masculinity score	−4.21	−5.30	−3.11	−0.29	<0.001	0.97	1.04	0.10		<0.001
		Sex	2.71	0.74	4.68	0.11	<0.01	0.96	1.05
		Moderate PA MET-m/w	−0.00	−0.00	0.00	−0.10	0.013	0.98	1.02
Low activity level
	Model 1	Masculinity score	−4.64	−6.19	−3.09	−0.30	<0.001	1.00	1.00	0.08	2.04	<0.001
	Model 2	Masculinity score	−4.93	−6.47	−3.39	−0.31	<0.001	0.98	1.02	0.10		<0.001
		Sex	4.19	1.43	6.95	0.15	<0.001	0.98	1.02
	Model 3	Masculinity score	−4.84	−6.37	−3.31	−0.31	<0.001	0.98	1.02	0.12		<0.001
		Sex	4.86	2.10	7.63	0.17	<0.001	0.96	1.04
		Moderate PA MET-m/w	−0.01	−0.02	0.00	−0.15	<0.001	0.97	1.03
Moderate activity level
	Model 1	Masculinity score	−2.06	−4.01	−0.11	−0.15	0.04	1.00	1.00	0.02	2.15	0.04
High activity level
	Model 1	Masculinity score	−4.89	−7.11	−2.67	−0.51	<0.001	1.00	1.00	0.24	2.67	<0.001
	Model 2	Masculinity score	−5.13	−7.24	−3.01	−0.53	<0.001	0.99	1.01	0.32		<0.001
		Accumulated recovery cases	−0.01	−0.02	0.00	−0.29	0.01	0.99	1.01
	Model 3	Masculinity score	−4.91	−6.84	−2.98	−0.51	<0.001	0.99	1.01	0.44		<0.001
		Accumulated recovery cases	−0.07	−0.10	−0.04	−1.63	<0.001	0.06	15.45
		Accumulated death cases	10.58	4.59	16.58	1.38	<0.001	0.06	15.49

Notes: For gender, females were coded with zero, and males were coded with 1.

For all participants, model 3 that included masculinity score, gender, and moderate-intensity PA as predictors reflects the highest degree of interpretation (Adjusted R2 = 0.10). In this model, masculinity score (β = –0.29) and MET-m/w of moderate-intensity PA (β = –0.10) were negatively predictive of depression and gender (females were coded as zero, males as 1) (β = 0.11) were positively predictive of depression.

For the LAL group, model 3 that included masculinity score, gender, and MET-m/w of moderate-intensity PA as predictors reflects the highest degree of interpretation (Adjusted R2 = 0.12). In this model, masculinity score (β = –0.31) and MET-m/w of moderate-intensity PA (β = –0.15) were negatively predictive of depression and gender (females were coded as zero, males as 1) (β = 0.17) were positively predictive of depression.

For the MAL group, there was only one significant model (F (1, 199) = 4.35, p < 0.05, Adjusted R2 = 0.02), in which masculine score was included as predictor variable negatively predicted (β = –0.15) depression.

For the HAL group, model 3 that included masculinity score, the number of recovered patients, and death toll as predictors showed the highest degree of interpretation (Adjusted R2 = 0.44). In this model, the masculinity score (β = –0.51) and the number of recovered patients (β = –1.63) negatively predicted depression, and death toll (β = 1.38) could positively predict depression.

Discussion

Results from the present study indicate that more than one-third of college students scored over the cut-off points for clinically relevant depression (≥16, CES-D scale), which was consistent with recent publications examining mental health status during the COVID-19 pandemic among young people and college students.12,13 Such figures of depression are much higher than the prevalence rate (3.6%) at the same age population before the COVID-19 pandemic.11 Moreover, results indicate that less time spent on physical activities, which were far below weekly 150–300 minutes of moderate-to-vigorous physical activities during the COVID-19 pandemic recently recommended by the American College of Sports Medicine (ACSM, 2020). Accordingly, extremely increased depression risk and closure of almost all public places (parks, gym, and schools, etc.) could have significantly decreased time in PA during the COVID-19 pandemic.

More than half of participants were categorized as “low” activity level group in terms of self-reported activity duration and frequency who reported to have a higher level of depression as compared to those in the “high” activity level group. These results were consistent with previous studies investigating the beneficial effects of PA,23,48,49 indicating that even under the impact of the COVID-19, PA can also confer protection against increasing levels of depressive symptoms. Furthermore, for the LAL group, moderate MET-m/w showed a robust and significant negative correlation with depressive symptoms regardless of their sedentary time, age, grade, and COVID-19 situations. Such findings suggest that individuals should spend more moderate-intensity PA times per week for reducing potential risks of developing depressive symptoms. In contrast, the association between sedentary and depression was not significant. The social withdrawal hypothesis suggests that sedentary behaviors might be associated with an increased risk of depression, which may be attributed to reduced social interactions.50 Perhaps this was not the case among the participants in this study, for whom social withdrawal and distancing acted as healthier behavior that reduced the risks of being infected with the COVID-19 and associated psychosocial distress, so the negative effect of higher sedentary might be offset.

In this study, the gender role directly predicted the onset of depression among college students. It is worth noting that this age group who had higher masculinity traits was associated with a lower risk of being depressed. Such results are supported by previous studies,51,52 suggesting that masculine traits that include physical and emotional self-control may help them to effectively manage and regulate negative emotions during episodes of increased stress exposure such as the pandemic. Multidimensional masculine norms influence individual psychological responses and depressive symptomology over time,53 which may inform how pre-existing masculine norms may have contributed to the psychosocial resilience among the study participants. Recently, an imaging study indicated that individuals with higher scores on masculine traits were less sensitive to negative emotional stimuli,54 which may in turn, reduce the susceptibility to depression. Most importantly, the predication patterns varied across activity levels.

For the LAL group, less moderate MET-m/w together with low scores on masculinity traits might increase the risk of depression. Moreover, for the HAL, the number of recovered patients and deaths is also predictive of depression in the negative and positive direction, respectively. As indicated by a recent study, specific local outbreak situation (like the number of confirmed cases and recovered cases) was associated with depression in general populations.55 Of note, effective control (increased number of recovered cases and decreased number of deaths) of the outbreak could help reduce the risk of developing depression among physically active college students. However, further investigations are still needed to keep a close eye on the dynamic changes for monitoring the onset of depression during life-changing events. Accordingly, college students like LAL who were susceptible to depression are more dependent on their internal factors including masculine traits and moderate-intensity PA frequency. In contrast, college students who showed a higher level of physical activity are also susceptible to external factors such as pandemic situations. Collectively, results indicate that individuals with low masculinity traits who showed a “low” level of physical activity may have the greatest vulnerability regardless of gender, which suggests that timely and targeted psychosocial interventions should be implemented for this group during the COVID-19 pandemic.

Practical Implications

Public health, education departments and school administrators should jointly pay more attention to college students with a lower level of masculinity traits (not specifically females) because they are observed to associate with increased risk of developing depression. Regular moderate-intensity exercises (eg, walking or jogging on a treadmill, using an elliptical trainer, cleaning heavy like washing windows, vacuuming, mopping) should be highly encouraged for depression prevention. Furthermore, findings exercise type that fits those college students with lower masculinity trait should be prioritized as well.

Effective control of the COVID-19 pandemic must be the most critical and ultimate goal of the current policies, however, at the same time measure that may allow attenuating the detrimental effects on mental health such as increasing levels of depression need to flank the measure that aims at reducing the transmission rates by physical distancing. Therefore, international cooperation is essential to stop the spread of the COVID-19 and to prevent long-term negative effects on mental health.

Conclusion

An increased risk of developing depression is associated with a decreased level of PA during the COVID-19 pandemic. Results suggest that more moderate-intensity MET-m/w are conducive to achieving mental health benefits college students during this life-changing event. Furthermore, college students endorse fewer masculine traits (regardless of gender) had the greatest susceptibility to depression during this outbreak. The decreased number of deaths and increased recovered patients of COVID-19 would directly lead to a reduced risk of being depressed among general populations.