Shuqi Yao
Guochen Wang2
Longhai Zhang
Puyan Chi- 1School of Physical Education, Shaanxi Normal University, Xi’an, China
- 2SMS Longhua School, Shenzhen, China
- 3Department of Physical Education, Shanghai Maritime University, Shanghai, China
- 4School of Physical Education, Qinghai Nationalities University, Xining, Qinghai, China
Introduction: To compare the immediate effects of a single bout of virtual reality (VR) aerobic exercise versus traditional aerobic exercise on depressive mood in middle-school students and to explore the underlying Electroencephalogram (EEG) microstate and power-spectrum mechanisms.
Methods: Forty middle-school students were classified into depressed and healthy groups based on the PHQ-9 and completed 15 min of moderate-intensity conventional cycling and VR cycling in a crossover design. Mood was assessed with the Brief Mood Scale (BFS) before and after each intervention, and resting-state EEG was recorded. EEG signals were processed using power spectrum analysis and microstate analysis, and correlation analysis was conducted between BFS questionnaire scores and microstate parameters.
Results: Both interventions significantly increased vigor and pleasure while reducing depression and lethargy (p < 0.001); VR was superior to traditional cycling in improving vigor, depression, and lethargy (p < 0.01). At baseline, the depressed group showed elevated occurrence and contribution of microstate C with abnormal centroids. After exercise, microstate distributions normalized in both groups; VR specifically reduced microstate C occurrence (p < 0.05), prolonged microstate D duration (p < 0.05), and enhanced B → C and D → B transitions (p < 0.05). Post-VR B → C transition rate correlated negatively with depression (r = −0.462), whereas microstate D occurrence correlated positively with pleasure (r = 0.450). Relative theta, alpha, and beta power decreased after exercise (p < 0.01), without additional VR-mediated power-spectrum gains.
Conclusion: A single bout of VR-assisted aerobic exercise more effectively enhances immediate mood in middle-school students with depressive symptoms than traditional cycling, likely owing to stronger modulation of microstate dynamics linked to emotion and attention networks and of oscillatory activity in specific frequency bands.
1 Introduction
Adolescence is a crucial stage for physical, psychological, and social development of individuals (Blakemore, 2008; Crone and Dahl, 2012; Dumith et al., n.d.). However, modern teenagers are experiencing unprecedented pressures from academics, social situations, and family dynamics. This has led to a rise in mental health issues, with depression being a significant concern (Huang X. et al., 2023; Farmakopoulou et al., 2024; Muthelo et al., 2024; Samsudin et al., 2024; Urbańska-Grosz et al., 2024). Global research shows that depression is now a major cause of disability in adolescents, with its prevalence increasing over the past decade (Mojtabai et al., 2016). Adolescent depression is marked by persistent feelings of sadness, loss of interest, cognitive difficulties, sleep problems, and challenges in relationships (Whitehouse et al., 2009; Asarnow et al., 2011; Lovato and Gradisar, 2014). In severe cases, it can even result in self-harm and suicidal tendencies. It is therefore imperative to find safe, effective, and easily accessible interventions to address symptoms of adolescent depression.
The adolescent developmental stage, particularly among early high school students aged 15–16, warrants focused investigation (Nelson et al., 2008). This period extends from and builds upon the critical developmental window for prefrontal cortex maturation–a core neural substrate for emotion regulation (Wang et al., 2024) and cognitive control that remains incompletely mature yet retains considerable neuroplasticity into later adolescence (Koechlin et al., 2003; MacDonald et al., 2000). This sustained plasticity provides a continued physiological basis for modulating brain function through interventions such as physical exercise (Thapar et al., 2012). Furthermore, this phase immediately follows the first peak incidence of depressive symptoms, making the study of this age group crucial for understanding the persistence and progression of early-onset vulnerabilities (Desai Boström et al., 2024; Kessler et al., 2001). Consequently, implementing and evaluating interventions within this population holds significant practical importance.
Physical exercise is commonly recommended as a non-pharmacological intervention for depression because of its affordability, minimal side effects, and easy accessibility (El-khalek et al., 2025). Many studies have shown that regular aerobic exercise can enhance mood, decrease anxiety, and benefit physiological functions such as the cardiovascular system (Ruiz et al., 2011). However, traditional exercise regimens may be seen as dull, posing a challenge for individuals with depression who struggle with motivation and enjoyment.
Virtual Reality (VR) technology has emerged as a promising solution to the challenge of engaging users in physical exercise (Jones and Wheat, 2023). By immersing users in interactive 3D virtual environments, VR enhances user engagement and enjoyment (Wagler and Hanus, 2018; Kim and Ko, 2019; Keller et al., 2022). Combining VR with exercise, known as “exergames,” not only maintains the physical benefits of sports training but also makes exercise more appealing and fun through gamification. Research has demonstrated that the audiovisual stimulation provided by VR has been shown to positively impact emotions by diverting attention and enhancing positive emotional experiences (Lavoie et al., 2021; Li et al., 2021; Tao et al., 2024). The application of VR technology in the intervention of depression is a novel and promising auxiliary treatment method (Liu et al., 2025). It provides patients with a safe and engaging intervention experience by creating an immersive and controllable digital environment (Wang et al., 2025a). VR immersive experience modulates large-scale brain network functions through multiple neural mechanisms. The core “sense of presence” originates from the significant suppression of the default mode network (DMN)–when users fully engage their attention in the virtual environment, novel and abundant external sensory inputs continuously occupy cognitive resources, thereby reducing the allocation of resources for internal processes such as self-referential thinking and memory retrieval. This downregulation of DMN activity aligns with neural manifestations observed during high cognitive load tasks (Sridharan et al., 2008), and the depth of immersion is positively correlated with the degree of DMN suppression (Chang et al., 2025). Simultaneously, the continuously changing salient stimuli in the virtual environment enhance the activation of the salience network (SN), enabling it to persistently detect and coordinate the brain’s responses to external events (Hermans et al., 2014). The heightened activity of the SN further facilitates the switching of cognitive resources to the executive control network (ECN), which is consistently activated through goal-directed interactive tasks such as virtual navigation and object manipulation (Menon, 2011; Petersen et al., 2022). Moreover, the VR experience reshapes the dynamic interactions between these networks: the typical antagonistic relationship between the DMN and ECN is reinforced, resulting in stronger anti-correlations (Seeburger et al., 2024), while the SN acts as a dynamic hub, flexibly regulating the brain’s transition from a resting state to a task state in response to real-time fluctuations in event salience and task demands within the virtual environment (An et al., 2024). This synergistic modulation across multiple networks collectively constitutes the neural basis of VR immersion.
Although VR exercise shows promising prospects in improving mental health, its advantages over traditional exercise and the underlying neural mechanisms still need further exploration. Electroencephalogram (EEG) technology provides a high temporal resolution (Coronel-Oliveros et al., 2024; Liu et al., 2024), non-invasive tool (Hodnik et al., 2024; Xie et al., 2024) to objectively assess changes in emotional states (Gong et al., 2024) and explore their neurophysiological basis (Babiloni et al., 2020; Hack et al., 2024). Previous studies have found significant abnormalities in EEG microstates in patients with depression, such as changes in parameters of microstate C related to the default mode network (DMN) and microstate D related to attention networks (Zhao et al., 2022), which are considered as electrophysiological manifestations of neural circuit dysfunction in depression. Additionally, depression is also associated with abnormalities in power spectra of specific frequency bands (e.g., Theta waves) (Bress et al., 2013). However, whether VR exercise can more effectively regulate the brain functional state of depressed adolescents compared to traditional exercise is a scientific question that urgently needs to be answered.
Therefore, this study aims to compare the immediate emotional effects of single-session VR aerobic exercise and traditional aerobic exercise on high school students with depressive symptoms through a randomized crossover controlled trial. Additionally, we will use resting-state EEG technology, specifically microstate and power spectrum analysis, to explore the differences in the effects of these two interventions on brain neural activity, with the goal of uncovering the potential neural mechanisms underlying the improvement of mood through VR exercise.
2 Materials and methods
2.1 Participants
This study distributed health questionnaires (PHQ-9) (Ettman et al., 2020) to a total of 1608 students in the first and second grades of a high school in Xi’an. The inclusion criteria were: (1) students aged 15–16 years old; (2) right-handed; (3) normal vision or corrected vision. Additional inclusion criteria for the depression group were: PHQ-9 total score ≥ 10 points, confirmed by structured interviews by qualified school counselors, indicating that significant depressive emotional distress exists. The inclusion criteria for the normal control group were: PHQ-9 total score ≤ 4 points. Participants with intermediate PHQ-9 scores (5–9) were excluded from the study. The exclusion criteria were: (1) a history of mental illness (excluding depressive symptoms), neurological disorders, or a family history; (2) taking medications that may affect the central nervous system; (3) the presence of cardiovascular or other physical illnesses unsuitable for moderate-intensity exercise; (4) contraindications related to VR use, such as severe motion sickness.
Using G*Power3.1.9.7 software for sample size estimation, with f = 0.25, α = 0.05, Power = 0.90 (Ezegbe et al., 2019), the calculated minimum sample size for a single group is 18 cases. Considering sample loss, ultimately, 40 high school students were involved in this research. The depression group and the healthy control group consisted of 20 individuals each, with 12 females and 8 males in each group. There were no statistically significant differences in age, body mass index (BMI), resting heart rate, and PHQ-9 scores between the two groups (p > 0.05) (Table 1). After fully comprehending the research procedures and potential risks. The Ethics Review Committee of Shaanxi Normal University (Approval No. 202516035).
Table 1. Comparison of baseline characteristics of participants.
2.2 Experimental procedure
In this study, two groups of participants were required to complete two different types of exercise interventions as follows: (1) “Pure aerobic exercise”: participants performed moderate-intensity (50% ∼ 80% HRmax) cycling exercise (15 min) on a power bike (Sweat Horse JTB616); (2) “VR aerobic exercise”: participants wore a VR head-mounted display (PICO 4, resolution 4320 × 2160, refresh rate 90 Hz) while cycling on the same model of stationary bike. The intervention employed a unified cycling control mode, set in a virtual seaside city environment of the “independent exploration” type, and incorporated the “Ride Bar 2” game program. The VR setup featured the following standardized characteristics: cycling speed was real-time mapped to the forward movement speed in the virtual environment, with the system providing real-time distance prompts and zone-triggered visual feedback (Yin et al., 2016); exercise intensity was strictly controlled within the same range as the pure aerobic exercise condition via heart rate monitoring (Gillen et al., 2012); the visual component presented a seaside landscape under daytime lighting conditions, accompanied by auditory stimuli that included cycling rhythm sounds, waves, and occasional seagull calls to create a multi-channel sensory experience (Moullec et al., 2022). The two protocols were fully consistent in terms of exercise duration, equipment configuration, and intensity control, with the uniqueness of the VR condition lying solely in the enhanced experiential dimension achieved through multi-sensory immersion and gamified interaction.
Participants were randomly and evenly assigned to the two intervention modes, and after a washout period (Mills et al., 2009; Bougrine et al., 2024) (6 h between the two exercise interventions), they underwent the other mode of exercise intervention. EEG data were collected at rest before the first intervention, after “Traditional aerobic exercise,” and after “VR aerobic exercise.” A physical activity monitor (GT3-X+, USA) monitored heart rate changes during exercise. Additionally, after both exercise interventions, participants must fill out the Brief Mood Scale (BFS) questionnaire (Altinsoy and Dikmen, 2025).
2.3 Data acquisition and processing
Electroencephalogram data were collected using a 32-channel EEG signal acquisition system from Neuroscan (Brain Vision Recorder, Neuroscan, USA). The electrode cap followed the international 10–20 system for electrode placement (Annaka et al., 2024; Desdentado and Pollatos, 2025). The experiment occurred in a quiet, soundproof room with reduced indoor lighting to minimize distractions. Participants washed their scalp before the experiment, and electrode impedances were adjusted to below 5 kΩ. Participants were instructed to avoid movements such as swallowing, eye movements, frowning, leg shaking, and to maintain a quiet, upright sitting position with their eyes closed and remain awake (Ahsan Awais et al., 2024). Resting-state EEG data were collected before the participant’s first movement and immediately after two subsequent movements, each lasting 6 min (Ng et al., 2024). The EEG signal was sampled at a rate of 1024 Hz, with the mastoid electrodes (M1/M2) used as a reference.
Data processing steps include: (1) importing continuous EEG data; (2) localizing channels; (3) downsampling to 256 Hz; (4) bandpass filtering (0.5–100 Hz) and notch filtering (48–52 Hz) (Alexander et al., 2024); (5) checking waveform plots to identify bad electrodes and interpolating abnormal electrodes; (6) segmenting continuous EEG data into 2-s epochs; (7) using Independent Component Analysis (ICA) to correct data artifacts caused by eye movements, blinks, muscle activity (Electromyogram, EMG), heart activity (Electrocardiogram, ECG), or other non-physiological sources. The selection of independent components was performed manually based on topographic and spectral features; (8) removing electrode data with amplitudes exceeding ±100 μV (Kappenman and Luck, 2010); (9) manually removing segments with large drift; (10) saving the processed EEG data.
2.4 EEG microstate analysis
Electroencephalogram microstate analysis is based on the Temporal Atomized Agglomerative Hierarchical Clustering (T-AAHC) algorithm (Chen et al., 2021). (1) Calculate the Global Field Power (GFP), which is the standard deviation of the average potential across all electrodes with a common reference. This is defined as:
Among them, ui is the voltage value on the i electrode, is the average voltage of all electrodes, and N is the number of electrodes in the topography map. GFP is a single quantification index of the topography map for each time point, reflecting the overall intensity of brain electrical activity. Studies have shown that the topographical structure of EEG remains stable when GFP is high, while it changes rapidly near the minimum GFP. Therefore, conducting microstate analysis at GFP peak values can provide the optimal signal-to-noise ratio. In the analysis process, ignoring the polarity of the topography map, each result cluster represents a microstate category. (2) Utilize the microstate analysis software Cartool 3.70 to submit each subject’s EEG topographic map data to the T-AAHC algorithm. This algorithm can identify clusters of topographic maps with similar configurations. The arrangement of microstate categories is determined by the spatial correlation between group-level microstate topographic maps (Chu et al., 2020). (3) Based on the maximum spatial correlation coefficient between each original topographic map and the group-level microstate topographic map, the extracted global microstate topographic maps are “matched” to each subject’s original EEG data stream. The microstate category that matches most at each time point is determined based on spatial correlation (Li et al., 2023). Calculate the following four core microstate parameters (He et al., 2021): (1) Duration: the average duration (ms) of a single microstate appearing continuously. (2) Occurrence: the average number of occurrences of a certain microstate per second. (3) Contribution: the percentage of total time a certain microstate occupies. (4) Transition: the probability of transitioning from one microstate to another.
2.5 EEG power spectrum analysis
Electroencephalogram spectral analysis is conducted using MATLAB software for each participant and segment. The analysis method involves using a batch processing approach, with the following steps: For each participant and each segment, Fast Fourier Transform (FFT) is used to obtain the EEG spectrum, transforming the signal from the time domain to the frequency domain, generating a power spectrum (μV2) ranging from 1 ∼ 100 Hz. FFT relies on the Discrete Fourier Transform (DFT), where for an N-item complex sequence x(s), its DFT is:
X(K) represents the data after DFT transformation, x(s) is the sampled analog signal, where x (s) can be imaginary. In reality, x (s) is all real signals, with the imaginary part being 0. Then the formula can be expanded as:
In the time domain, n represents the s-th sample; x represents the signal time series (s = 0,1,2,…,S−1); X stands for the h-th frequency domain representation of the signal x; S represents the signal x corresponding to the k-th frequency component; k corresponds to the k-th frequency component (k = 0, 1, 2…, S−1).
Calculate the absolute power of each electrode in different frequency bands (Delta: 1–4 Hz, Theta: 4–8 Hz, Alpha: 8–13 Hz, Beta: 13–30 Hz) in the MATLAB command window. Then, calculate the relative power (percentage of power in a specific frequency band relative to the total power). Increased power in the slow wave δ frequency band during the resting state is associated with enhanced brain plasticity (Castellanos et al., 2010). Frontal θ wave activity may reflect the level of attention (Gevins et al., 1997; Ishii et al., 1999). α waves, as the fundamental rhythm of brain electrical activity, are closely related to the excitatory processes of cortical neurons, alertness (Silber et al., 2007), and working memory (Hsueh et al., 2016). β waves are involved in oscillatory synchronization between motor and somatosensory cortices (Brovelli et al., 2004), and sensorimotor interactions (Tan et al., 2016; Spitzer and Haegens, 2017). Depression symptoms are closely linked to executive and emotional regulation functions. Hence, the experiment focuses on frontal and central regions as regions of interest, and C3, C4, F3, F4, FP1, and FP2 as electrodes of interest.
2.6 Statistical analysis
Statistical analyses were performed using SPSS 26.0 software. The data were organized in Excel 2010 and then imported into the statistical software. Normality was assessed using the Shapiro-Wilk test. Data that followed a normal distribution were presented as mean ± standard deviation (M ± SD). A 2 (Group: normal group, depression group) × 3 (Condition: pre-exercise, Post-Traditional Aerobic Exercise, Post-VR-Aerobic Exercise) mixed-design analysis of variance was used to analyze the scores of BFS dimensions and EEG data. The group was the between-subjects variable, while the condition was the within-subjects variable. Simple effects analysis was conducted in case of a significant interaction (p < 0.05). For significant main effects, independent samples t-tests were used for post hoc comparisons. Effect sizes were indicated by partial η2 for Analysis of Variance (ANOVA) and Cohen’s d for t-tests. The Pearson correlation coefficient was utilized to examine the relationship between changes in the mood scale depression index and changes in EEG microstate temporal parameters and the transition rate index. The significance level for all analyses was set at α = 0.05. The Greenhouse-Geisser correction was applied for variables that did not meet the sphericity assumption. The Bonferroni correction was employed for post-hoc tests.
3 Results
3.1 BFS mood scale results
Results of the two-way mixed design ANOVA (Table 2) revealed that BFS showed significant differences in the dimensions of activity, pleasure, depression, and vigor, with significant main effects of condition, Group, and the interaction between group and condition (p < 0.05). Further simple effects analysis revealed: in terms of within-subject effects analysis focusing on condition effects, compared to before exercise, the depression group showed significant increases in both arousal and pleasure ratings after aerobic exercise (p < 0.001), while depression and lethargy ratings significantly decreased (p < 0.001). Additionally, compared to aerobic exercise alone, the depression group showed significant increases in arousal ratings after VR aerobic exercise and significant decreases in depression and lethargy ratings (p < 0.01). There were no significant differences in the normal group (p > 0.05). Analysis of inter-subject effects with group as the main factor: the depression group exhibited significantly lower arousal and pleasure scores compared to the normal group at pre-exercise (p < 0.001), after aerobic exercise alone (p < 0.01), and after VR aerobic exercise (p = 0.004, p = 0.001). Conversely, the depression group showed significantly higher scores for depression and lethargy than the normal group at pre-exercise (p < 0.001), after aerobic exercise alone (p < 0.001), and after VR aerobic exercise (p < 0.01). No significant differences were observed between the two groups in the remaining components of the scale (p > 0.05).
Table 2. Results of the BFS mood scale.
3.2 EEG microstate results
3.2.1 Microstate topographic maps
As shown in Figure 1, before the exercise intervention, both groups exhibited a similar distribution of centroids in microstates A and D. The centroids of microstate A were located in the right frontal lobe and left occipital lobe, while those of microstate D were in the central frontal lobe and occipital lobe for both groups. However, there were significant differences in the centroid distribution of microstates B and C. In the control group, the centroid of microstate B was in the left frontal lobe and right occipital lobe, while in the depression group, it was in the parietal lobe and left frontal lobe. Similarly, the centroid of microstate C in the control group was in the anterior frontal lobe and the occipital lobe. In contrast, in the depression group, it was in the central region and the left occipital lobe.
Figure 1. Microstate topographic maps.
After pure aerobic exercise, the centroids of microstates A, B, C, and D in both groups tend to converge toward specific regions in the brain. The centroids of microstate A are located in the right frontal lobe and the left occipital lobe, the centroids of microstate B are located in the left frontal lobe and the right occipital lobe, the centroids of microstate C are located in the frontal lobe and the occipital lobe, and the centroids of microstate D are evenly distributed in the central frontal area and the occipital lobe. Following the VR aerobic exercise intervention, the centroids of microstate topography maps A, B, C, and D in both groups also show similar distributions to those observed after aerobic exercise. This indicates that aerobic exercise intervention can significantly improve the brain’s functional state of depressed high school students, normalizing the distribution of microstate centroids.
3.2.2 Microstate temporal parameters
For Duration A, the main effect of condition is significant [F (1.699, 64.552) = 8.623, p = 0.001, η2 = 0.185], with a highly significant difference between pre-exercise and aerobic exercise alone (p < 0.01). The condition × group interaction is not significant (p = 0.365), but simple effect analysis shows a significant difference for the normal group between pre and post aerobic exercise (p < 0.05). For Duration B, neither the main effect of condition nor the condition × group interaction is significant. For Duration C, the main effect of condition is significant [F (2, 76) = 3.835, p = 0.026, η2 = 0.092], with a significant difference between aerobic exercise and pre-exercise (p = 0.048). The condition × group interaction is not significant (p = 0.495). For Duration D, the main effect of condition is significant [F (2, 76) = 4.019, p = 0.022, η2 = 0.096]. The condition × group interaction is not significant (p = 0.317), but simple effect analysis reveals a significant difference between pre-exercise for the normal group and the depression group (p = 0.032), and a significant difference between pre - and post - aerobic exercise for the depression group (p = 0.047) (Figure 2).
Figure 2. Duration parameters for microstates A–D in the two groups under three. *p < 0.05.
Occurrence A and Occurrence B: neither the main effect of condition nor the interaction effect of group × condition was significant. Occurrence C: the main effect of condition was not significant [F (2, 76) = 2.041, p = 0.137, η2 = 0.051]; however, there was a significant interaction between group and condition [F (2, 76) = 4,430, p = 0.015, η2 = 0.104]. Simple effect analysis revealed that the depression group scored significantly higher than the normal group before exercise (p = 0.016); the depression group showed a significant decrease in VR aerobic exercise post-exercise compared to pre-exercise (p = 0.027). Occurrence D: the main effect of condition was statistically significant [F (1.629, 61.898) = 10.019, p < 0.01, η2 = 0.209]. Post hoc comparisons indicated a significant improvement in VR aerobic exercise post-exercise compared to pre-exercise (p = 0.039); aerobic exercise also showed a significant improvement post-exercise compared to pre-exercise (p < 0.01). The interaction between group and condition was not significant (p = 0.281). Further simple effect analysis showed that the depression group scored significantly lower than the normal group before exercise (p = 0.014); the normal group showed a significant improvement in aerobic exercise post-exercise compared to pre-exercise (p = 0.027); the depression group showed a very significant improvement in aerobic exercise post-exercise compared to pre-exercise (p < 0.01); and the depression group also showed a significant improvement in VR aerobic exercise post-exercise compared to pre-exercise (p = 0.033) (Figure 3).
Figure 3. Occurrence parameters for microstates A–D in the two groups under three. *p < 0.05; **p < 0.01.
Contribution A: the main effect of condition and the group × condition interaction were both not significant. Contribution B: the main effect of condition and the group × condition interaction were both not significant (p = 0.922, p = 0.053). Still, simple effect analysis showed a significant difference between the pre-exercise measurements of the depression group and the normal group (p = 0.037). Contribution C: the main effect of condition was significant [F (2, 76) = 4.917, p = 0.010, η2 = 0.115], with a significant difference between VR aerobic exercise and pre-exercise measurements (p = 0.003). The group × condition interaction was not significant (p = 1.542). Still, simple effect analysis revealed a significant difference between the pre-exercise measurements of the depression group and the normal group (p = 0.048), as well as significant differences between post-aerobic exercise and post-VR aerobic exercise compared to pre-exercise measurements in the depression group (p = 0.042, p = 0.002) (Figure 4). Contribution D: the main effect of condition was significant [F (1.386, 52.686) = 4.337, p = 0.030, η2 = 0.102], with a significant difference between aerobic exercise and pre-exercise measurements (p = 0.018). The main effect of group was significant [F (1, 38) = 6.460, p = 0.015, η2 = 0.145], with the depression group scoring higher than the normal group (p = 0.015). The group × condition interaction was not significant (p = 0.955).
Figure 4. Contribution parameters for microstates A–D in the two groups under three. *p < 0.05; **p < 0.01.
3.2.3 Microstate transition probabilities
According to Figure 5 there are no significant findings for the conversion rates from A → B and A → D: neither the main effect of condition nor the group × condition interaction (p > 0.05). As for the conversion rate from A → C, both the main effect of condition and the group × condition interaction are not significant (p = 0.612, p = 0.135), but further analysis reveals a substantial difference between the depression group and the normal group in pre-exercise (p = 0.027).
Figure 5. Microstates A–D transition rates: two groups across three conditions. *p < 0.05; **p < 0.01.
Conversion rates B → A and B → D: neither the main effect of condition nor the group × condition interaction was significant (p > 0.05). Conversion rate B → C: the main effect of condition was significant [F (2, 76) = 6.460, p = 0.003, η2 = 0.145]. Further simple effect analysis revealed that there were significant differences between the VR aerobic exercise post-exercise and pre-exercise, as well as between the VR aerobic exercise post-exercise and the pure aerobic exercise post-exercise (p < 0.05); The group × condition interaction was significant [F (2, 76) = 4.865, p = 0.010, η2 = 0.113]. The pre-exercise differences between the depression group and the normal group were significant (p = 0.044), and the differences between the post-exercise of VR aerobic exercise and the pre-exercise, as well as between the post-exercise of VR aerobic exercise and the traditional aerobic exercise post-exercise in the depression group were both highly significant (p = 0.001, p = 0.002).
Conversion rate C → A: the main effect of condition showed significance [F (2, 76) = 3.357, p = 0.040, η2 = 0.081], while the group × condition interaction was not significant (p = 0.367). Further analysis of simple effects revealed a substantial difference in VR aerobic exercise post-exercise compared to pre-exercise for the depression group (p < 0.05). Conversion rates C → B, C → D: neither the main effect of condition nor the group × condition interaction was significant (p > 0.05). However, simple effect analysis showed significant differences between the depression and normal Groups in pre-exercise measurements for these two conversion rates (p < 0.05).
Conversion rates from D → A and from D → C: the main effect of condition is not significant, and the group × condition interaction is also not significant (p > 0.05). Conversion rate from D → B: the main effect of condition is substantial [F (2, 76) = 4.082, p = 0.021, η2 = 0.097]. Post hoc comparisons showed a significant difference in VR aerobic exercise compared to the pre-exercise (p < 0.05). The group × condition interaction is not essential, but simple effect analysis revealed a significant difference in the depression group after VR aerobic exercise compared to the pre-exercise (p < 0.05).
3.2.4 Correlations between microstates and mood scale scores
In comparison with the pre-exercise assessment, significant changes were observed in the post-exercise BFS mood scale in terms of liveliness, pleasantness, depression, and lack of vitality. After aerobic exercise alone, Occurrence D and Contribution C showed significant changes. Following VR aerobic exercise, Duration D, Occurrence C, Occurrence D, Contribution C, transition rate B → C, and transition rate D → B all exhibited significant changes.
During pure aerobic exercise, as shown in Figure 6, there was a moderate positive correlation between pre-exercise pleasure and Occurrence D (r = 0.464, p = 0.039), but this correlation disappeared after exercise (r = 0.129, p = 0.587). This suggests that pure aerobic exercise may weaken the linear relationship between the two by restructuring the brain’s functional network, supporting the “network flexibility hypothesis” of exercise improving mood – that mood enhancement may stem from optimizing brain information processing patterns from static network activation to dynamic network interaction, rather than changes in a single microstate parameter. The other BFS mood component scale indicators and microstate time parameters showed no statistical significance.
Figure 6. Correlations between mood scale scores and EEG microstate changes.
Before the VR aerobic exercise intervention, there was no significant correlation between the B → C transition rate and depression (r = 0.008, p = 0.973), indicating no clear relationship between the severity of depressive symptoms at rest and the dynamic switching ability of the brain microstates. However, after the VR aerobic exercise intervention, there was a significant moderate negative correlation between them (r = −0.462, p = 0.04), suggesting that individuals with higher levels of depression had lower B → C transition rates after exercise, indicating that VR aerobic exercise may improve emotional regulation in individuals with depression by enhancing the flexibility of brain functional networks. Before the exercise, there was a significant moderate negative correlation between the level of inactivity and Contribution C (r = −0.466, p = 0.038), indicating that individuals with more severe symptoms of inactivity had lower Contribution C rates at rest. However, this correlation disappeared after the exercise (r = 0.271, p = 0.248), suggesting that VR aerobic exercise weakened the linear relationship between inactivity and Contribution C, improving individuals’ energy levels and subjective vitality. Before the exercise, there was no significant correlation between pleasure and Occurrence D (r = −0.118, p = 0.622). However, after the exercise, there was a significant moderate positive correlation between them (r = 0.450, p = 0.046), indicating that individuals with higher levels of pleasure had higher occurrence rates of Occurrence D after the exercise.
3.3 EEG relative power results
Theta frequency band: the main effect of the condition was highly significant for channels C3, C4, F3, F4, FP1, and FP2 (p < 0.01). Post hoc comparisons showed that the EEG values after pure aerobic and VR aerobic exercise were significantly lower than pre-exercise values (p < 0.01) for all channels. There was no significant interaction between group and condition for any channel. Regarding group differences, the depression group had significantly higher EEG values than the normal group on channels F4 (p = 0.004) and FP1 (p = 0.013). Simple effect analysis revealed that the depression group had significantly higher pre-exercise EEG values than the normal group on channels C3, C4, F3, F4, and FP1 (p < 0.05). The normal and depression groups showed significantly lower EEG values after pure aerobic exercise and VR aerobic exercise compared to their respective pre-exercise values (p < 0.01) on these channels. On channel FP2, both groups had significantly lower EEG values after pure aerobic exercise and VR aerobic exercise compared to their pre-exercise values (p < 0.01), with no significant difference between the groups in the pre-exercise values (p > 0.05) (Table 3).
Table 3. Relative power values of normal and depressed groups across three states.
Alpha frequency band: the main effects of C3, C4, F3, F4, FP1, and FP2 channels on the condition are all highly significant (p < 0.01). The EEG values after pure aerobic exercise and VR aerobic exercise are both significantly lower than the pre-exercise (p < 0.01). The group effect is significant for FP1 and FP2 channels (p = 0.022, p = 0.038), with the depression group showing significantly higher EEG values than the normal group (p < 0.05). There are no significant group × condition interactions for all channels. Simple effect analysis shows that in C3, C4, F3, F4, FP1, and FP2 channels, the EEG values after pure aerobic exercise and VR aerobic exercise for both normal and depression groups are significantly lower than their respective pre-exercise (p < 0.01), and both groups show consistent trends in changes in each condition.
The beta frequency band showed highly significant main effects on the state for C3, C4, F3, F4, FP1, and FP2 channels (p < 0.01). The EEG values after traditional aerobic exercise and VR aerobic exercise were both significantly lower than those before the exercise (p < 0.01). The main effect of group was highly significant for the FP1 channel (p = 0.001), with the EEG values of the depression group significantly lower than those of the normal group (p < 0.01). However, the main effect of group was not significant for the other channels. There were no significant interactions between group and state for any of the channels. Simple effect analysis showed that in the C3, C4, F3, F4, and FP2 channels, the EEG values of both the normal and depression groups after simple/VR aerobic exercise were significantly lower than those before the exercise (p ≤ 0.025). In the FP1 channel, the EEG values of the depression group after simple and VR aerobic exercise were not only significantly lower than their respective pre-exercise values (p < 0.01), but also significantly lower than the corresponding post-exercise values of the normal group (p ≤ 0.048).
4 Discussion
The study investigated the effects of a single session of VR aerobic exercise on middle school students with depressive symptoms. The results showed that VR aerobic exercise was more effective than traditional aerobic exercise in improving the emotional state of students with depression, particularly in increasing activity and reducing depression. EEG analysis revealed that VR aerobic exercise could normalize abnormal resting-state brain network dynamics in depressed adolescents and modulate the brain’s electrical power spectrum related to emotions and arousal.
4.1 Analysis of mood scale results
A single session of moderate-intensity aerobic exercise can effectively improve negative emotions and enhance positive emotions in adolescents with depression, as shown in this study. This aligns with the many past studies on the antidepressant benefits of exercise (Jacka et al., 2017; Firth et al., 2019; Thapa et al., 2023). Students in the depression group had higher activity scores after completing VR aerobic exercise compared to traditional aerobic exercise, with lower scores of depression and lethargy. Combining immersive virtual reality with exercise can significantly improve the emotional state of the depression group compared to exercise alone.
Virtual reality’s immersive and engaging experience may contribute to its superiority over traditional exercise environments (Bird et al., 2021). The virtual world created by VR can distract users from negative emotions and physical discomfort, enhancing pleasure during exercise (Ekkekakis, 2003). Gamified forms of exercise in VR can increase fun and motivation, which is beneficial for individuals with motivation deficits or symptoms of depression (Muñoz et al., 2022). The dynamic visual and auditory stimuli in VR environments provide richer sensory inputs, potentially positively regulating brain function and emotions (Tao et al., 2024). A recent study with college students found that an 8-weeks VR exercise intervention was more effective than traditional aerobic exercise in improving mood disturbances, tension, anger, depression, and vitality (Wang et al., 2025b). This evidence suggests that VR can amplify the psychological benefits of exercise.
4.2 Analysis of microstate findings
The research findings revealed that the depressed group displayed notable abnormalities in the distribution of microstate B and C centroids at baseline, supporting the theory of disrupted brain network functional connectivity in depression (Zhang et al., 2011; Guo et al., 2025). Following the exercise intervention, the centroid distribution of the four microstates in both groups appeared to normalize, suggesting that exercise plays a crucial role in regulating brain function states.
The study found that after VR exercise, the occurrence rate of microstate C decreased significantly in the depression group. Microstate C is linked to the salience network (Nishida et al., 2025), which includes the anterior cingulate cortex and insula (Menon and Uddin, 2010), important for emotion regulation (Fang et al., 2024; Roberts and Mulvihill, 2024). The decrease in the occurrence of Microstate C reflects how the VR experience redistributes attentional resources, reducing the excessive monitoring of internal negative states–a process potentially mediated by the normalization of anterior cingulate cortex and insula activity–thereby helping to alleviate rumination in individuals with depression. Additionally, the occurrence rate and contribution of microstate D increased after exercise, indicating improved function of the dorsal attention network. This improvement may be due to the sensory stimulation and task demands in the VR environment, which help shift individuals from internal contemplation to external task engagement. The depression group showed a significant increase in B → C transition rate after VR exercise, which was negatively correlated with the improvement of depressive symptoms. This suggests enhanced efficiency of information transmission between visual networks and salience networks. The increased flexibility in network collaboration reflected by the increase in transition rate may be a key mechanism in breaking the characteristic cognitive rigidity of depression. The occurrence rate of microstate D after VR exercise was positively correlated with feelings of pleasure (Bai et al., 2024), while the B → C transition rate was negatively correlated with the severity of depression (Luo et al., 2025). This links specific changes in network dynamics directly to emotional improvement (Weiss et al., 2024). In contrast, single aerobic exercise could also change some microstate parameters, but did not show similar systematic neural-behavioral correlations, highlighting the unique advantage of VR intervention in regulating brain network dynamics.
These findings collectively support the applicability of the “network flexibility hypothesis” in explaining the mechanism of how exercise improves mood. VR exercise not only regulates the activity levels of specific brain networks but also promotes flexible transitions of brain functional states by enhancing dynamic interactions between networks. The synergistic effect of normalizing network activity and enhancing attention network function drives individuals to shift from internal reflective modes to external engaged modes, providing a new theoretical perspective for understanding the neural mechanisms of how VR-enhanced exercise improves depressive mood.
4.3 Analysis of power spectrum results
The depression group showed significantly higher Theta relative power in the F4 and FP1 channels in the frontal lobe, and higher Alpha relative power in the FP1 and FP2 channels in the anterior frontal lobe compared to the control group. The increase in the Theta band may be related to emotional regulation dysfunction and weakened cognitive control in depression. While some studies suggest a decrease in theta activity among individuals with depression (Zandbagleh et al., 2024), other research has indicated an increase in theta power (Huang Y. et al., 2023). This discrepancy could be attributed to sample heterogeneity or variations in analysis techniques. The higher Alpha power in the anterior frontal lobe in the depression group may reflect an enhanced inhibitory or deactivating state in that region, consistent with reduced executive function and motivation in depression patients (Gotlib, 1998; Huang Y. et al., 2023). In addition, the depression group had significantly lower Beta power in the left frontal lobe (FP1 channel) compared to the healthy group, which may indicate an insufficient cortical arousal level or decreased cognitive processing ability, consistent with symptoms of cognitive slowing and lack of concentration in depression.
Both exercise interventions significantly decreased relative power in the Theta, Alpha, and Beta frequency bands for all participants. A decrease in low-frequency power is a common finding after acute aerobic exercise and is generally interpreted as an indicator of increased cortical arousal and activation. The baseline differences we observed, such as higher theta and alpha power in the depressed group, are also consistent with literature suggesting altered cortical arousal and information processing in depression. However, we found no significant additional effect of VR on these power spectrum measures. This negative finding is critical because it indicates that the superior mood-enhancing properties of VR are not simply due to a greater change in overall cortical arousal. Instead, the advantage of VR appears to lie in its ability to more specifically reorganize the temporal dynamics of large-scale brain networks, a phenomenon captured by microstate analysis but overlooked by traditional spectral analysis. This highlights the methodological value of employing microstate analysis to uncover subtle yet functionally significant changes in brain activity.
4.4 Study limitations and future directions
Although this study systematically explored the immediate effects of VR-enhanced aerobic exercise on mood and brain neural activity in depressed adolescents and yielded some positive findings, several limitations exist, which also point the way for future research.
(a) This study is an acute-effect investigation of a single bout of exercise intervention; thus, its results cannot address the long-term benefits of VR exercise. The satisfaction derived from VR might be transient, and the short-term design of this study cannot verify the persistence of its effects, nor can it assess whether long-term use might lead to dependency or potential negative mood impacts due to increased screen time. Future research needs to test the cumulative effects and long-term stability of VR exercise through trials involving multiple interventions and long-term follow-ups. (b) The intervention protocol in this study (e.g., VR scenario, exercise intensity, and duration) was standardized. While this design ensured internal validity, it did not account for individual differences in depressive symptom profiles, exercise preferences, and baseline neural activity levels, thus failing to achieve truly personalized intervention. Future research could explore adaptive intervention protocols based on individual FMS or EEG baseline characteristics. (c) Caution is warranted in interpreting the neural mechanisms. The observed changes in EEG microstates and power spectra were associated with mood improvement, but this correlation does not equate to causation. It remains difficult to fully disentangle to what extent these changes in neural indicators are due to the physiological effects of exercise itself versus the psychological effects induced by VR immersion. Future studies could incorporate neuroimaging techniques with higher spatial resolution, such as fMRI, to more precisely localize activity changes in relevant brain networks. (d) Regarding the washout period implemented in this study, although the 6-h interval was determined based on physiological recovery time and precedents in relevant research–and substantially exceeds the time required for the normalization of key physiological parameters–we acknowledge that residual effects on subjective psychological states such as mood may persist beyond this duration. This represents a limitation of the present study. Future research employing longer washout periods or between-subjects designs would help to clarify this issue. (e) The spectral analysis of this study focused on comparing the relative power of ICA-derived source components (reflecting large-scale brain networks) across groups and conditions. Therefore, we did not perform cluster-based or permutation-based statistical mapping of scalp topography. This approach was chosen to prioritize physiologically interpretable source-level signals and to account for the limited spatial resolution and volume conduction effects inherent in 32-channel EEG data. Future studies may employ high-density EEG or combined EEG-fMRI to better resolve the spatial topography and neural origins of these spectral changes.
5 Conclusion
This study demonstrates that a single session of VR-enhanced aerobic exercise is superior to traditional aerobic exercise in immediately improving mood states in adolescents with depressive symptoms. While both interventions enhanced activation and pleasure while reducing depression and lethargy, VR elicited significantly greater emotional benefits. At the neurophysiological level, VR exercise specifically modulated key microstate parameters–reducing the occurrence of salience-network-related microstate C, prolonging attention-network-related microstate D duration, and enhancing transitions between visual-salience networks. These changes correlated with mood improvement and suggest enhanced dynamic coupling between large-scale brain networks. In contrast, spectral power changes indicated general cortical arousal without modality-specific differences. The findings support VR-enhanced exercise as an effective, engaging intervention for acutely alleviating depressive mood in adolescents. Future studies should investigate its long-term efficacy, optimal dosing parameters, and personalized application through longitudinal designs and mechanistic exploration.
Data availability statement
The original contributions presented in this study are included in this article/supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Academic Committee of Shaanxi Normal University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
SY: Writing – original draft. GW: Conceptualization, Funding acquisition, Writing – original draft. TP: Formal analysis, Writing – original draft. LZ: Investigation, Writing – original draft. FM: Writing – review & editing. PC: Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by Qinghai Provincial Key R&D and Transformation Program (2023-SF-116, “Scientific Selection, Training Theory and Methods for High-Altitude “Big-Heart” Athletes”).
Acknowledgments
We thank all the funders for their support and the valuable contributions of all investigators and participants.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Ahsan Awais, M., Ward, T., Redmond, P., and Healy, G. (2024). From lab to life: Assessing the impact of real-world interactions on the operation of rapid serial visual presentation-based brain-computer interfaces. J. Neural Eng. 21:046011. doi: 10.1088/1741-2552/ad5d17
Alexander, K. E., Estepp, J. R., and Elbasiouny, S. M. (2024). Adaptive filtering with fitted noise estimate (AFFiNE): Blink artifact correction in simulated and real P300 data. Bioengineering 11:707. doi: 10.3390/bioengineering11070707
Altinsoy, C., and Dikmen, D. (2025). How are brain fog symptoms related to diet, sleep, mood and gastrointestinal health? A cross-sectional study. Medicina 61:344. doi: 10.3390/medicina61020344
An, Z., Tang, K., Xie, Y., Tong, C., Liu, J., Tao, Q., et al. (2024). Aberrant resting-state co-activation network dynamics in major depressive disorder. Transl. Psychiatry 14:1. doi: 10.1038/s41398-023-02722-w
Annaka, H., Hiraoka, T., Nomura, T., Annaka, H., Hiraoka, T., and Nomura, T. (2024). Effects of a virtual reality game with leap motion controller on brain activity related to attentional function in healthy adults - A pilot EEG study. Cureus 16:e71838. doi: 10.7759/cureus.71838
Asarnow, J. R., Porta, G., Spirito, A., Emslie, G., Clarke, G., Wagner, K. D., et al. (2011). Suicide attempts and nonsuicidal self-injury in the treatment of resistant depression in adolescents: Findings from the TORDIA study. J. Am. Acad. Child Adolesc. Psychiatry 50, 772–781. doi: 10.1016/j.jaac.2011.04.003
Babiloni, C., Barry, R. J., Başar, E., Blinowska, K. J., Cichocki, A., Drinkenburg, W. H. I. M., et al. (2020). International Federation of Clinical Neurophysiology (IFCN) – EEG research workgroup: Recommendations on frequency and topographic analysis of resting state EEG rhythms. Part 1: Applications in clinical research studies. Clin. Neurophysiol. 131, 285–307. doi: 10.1016/j.clinph.2019.06.234
Bai, Y., Yu, M., Li, Y., Bai, Y., Yu, M., and Li, Y. (2024). Dynamic neural patterns of human emotions in virtual reality: Insights from EEG microstate analysis. Brain Sci. 14:113. doi: 10.3390/brainsci14020113
Bird, J. M., Karageorghis, C. I., Baker, S. J., Brookes, D. A., and Nowicky, A. V. (2021). Ready exerciser one: Effects of music and virtual reality on cycle ergometer exercise. Br. J. Health Psychol. 26, 15–32. doi: 10.1111/bjhp.12445
Blakemore, S.-J. (2008). The social brain in adolescence. Nat. Rev. Neurosci. 9, 267–277. doi: 10.1038/nrn2353
Bougrine, H., Ammar, A., Salem, A., Trabelsi, K., Jahrami, H., Chtourou, H., et al. (2024). Optimizing short-term maximal exercise performance: The superior efficacy of a 6 mg/kg caffeine dose over 3 or 9 mg/kg in young female team-sports athletes. Nutrients 16:640. doi: 10.3390/nu16050640
Bress, J. N., Foti, D., Kotov, R., Klein, D. N., and Hajcak, G. (2013). Blunted neural response to rewards prospectively predicts depression in adolescent girls. Psychophysiology 50, 74–81. doi: 10.1111/j.1469-8986.2012.01485.x
Brovelli, A., Ding, M., Ledberg, A., Chen, Y., Nakamura, R., and Bressler, S. L. (2004). Beta oscillations in a large-scale sensorimotor cortical network: Directional influences revealed by Granger causality. Proc. Natl. Acad. Sci. U S A. 101, 9849–9854. doi: 10.1073/pnas.0308538101
Castellanos, N. P., Paúl, N., Ordóñez, V. E., Demuynck, O., Bajo, R., Campo, P., et al. (2010). Reorganization of functional connectivity as a correlate of cognitive recovery in acquired brain injury. Brain 133, 2365–2381. doi: 10.1093/brain/awq174
Chang, S., Suh, J., Chang, S., and Suh, J. (2025). The impact of digital storytelling on presence, immersion, enjoyment, and continued usage intention in VR-based museum exhibitions. Sensors 25:2914. doi: 10.3390/s25092914
Chen, J., Li, H., Ma, L., Bo, H., Soong, F., and Shi, Y. (2021). Dual-threshold-based microstate analysis on characterizing temporal dynamics of affective process and emotion recognition from EEG signals. Front. Neurosci. 15:689791. doi: 10.3389/fnins.2021.689791
Chu, C., Wang, X., Cai, L., Zhang, L., Wang, J., Liu, C., et al. (2020). Spatiotemporal EEG microstate analysis in drug-free patients with Parkinson’s disease. NeuroImage Clin. 25:102132. doi: 10.1016/j.nicl.2019.102132
Coronel-Oliveros, C., Gómez, R. G., Ranasinghe, K., Sainz-Ballesteros, A., Legaz, A., Fittipaldi, S., et al. (2024). Viscous dynamics associated with hypoexcitation and structural disintegration in neurodegeneration via generative whole-brain modeling. Alzheimer’s Dement. 20, 3228–3250. doi: 10.1002/alz.13788
Crone, E. A., and Dahl, R. E. (2012). Understanding adolescence as a period of social–affective engagement and goal flexibility. Nat. Rev. Neurosci. 13, 636–650. doi: 10.1038/nrn3313
Desai Boström, A. E., Cars, T., Hellner, C., and Lundberg, J. (2024). Adolescent and adult transitions from major depressive disorder to bipolar disorder. JAMA Psychiatry 81, 863–869. doi: 10.1001/jamapsychiatry.2024.1133
Desdentado, L., and Pollatos, O. (2025). Expressive suppression of emotions in bulimia nervosa: An electroencephalography study. J. Clin. Psychol. 81, 158–170. doi: 10.1002/jclp.23761
Dumith, S. C., Gigante, D. P., Domingues, M. R., and Kohl, H. W. (n.d.). Physical activity change during adolescence: A systematic review and a pooled analysis. Available online at: https://dx.doi.org/10.1093/ije/dyq272 (accessed August 14, 2025).
Ekkekakis, P. (2003). Pleasure and displeasure from the body: Perspectives from exercise. Cogn. Emot. 17, 213–239. doi: 10.1080/02699930302292
El-khalek, W. O. A. A., Abdulrahman, R. S., Mahmoud, H., Saeed, D., Elias, M. A. G., Ebid, A. A., et al. (2025). Effects of dietary advice, aerobic exercise, and virtual reality games on the quality of life of obese adolescent females. SPORT TK-EuroAmerican J. Sport Sci. 14:42. doi: 10.6018/sportk.660601
Ettman, C. K., Abdalla, S. M., Cohen, G. H., Sampson, L., Vivier, P. M., and Galea, S. (2020). Prevalence of depression symptoms in US adults before and during the COVID-19 pandemic. JAMA Netw. Open 3:e2019686. doi: 10.1001/jamanetworkopen.2020.19686
Ezegbe, B. N., Eseadi, C., Ede, M. O., Igbo, J. N., Anyanwu, J. I., Ede, K. R., et al. (2019). Impacts of cognitive-behavioral intervention on anxiety and depression among social science education students: A randomized controlled trial. Medicine 98:e14935. doi: 10.1097/MD.0000000000014935
Fang, F., Teixeira, A. L., Li, R., Zou, L., and Zhang, Y. (2024). The control patterns of affective processing and cognitive reappraisal: Insights from brain controllability analysis. Cereb. Cortex 34:bhad500. doi: 10.1093/cercor/bhad500
Farmakopoulou, I., Lekka, M., and Gkintoni, E. (2024). Clinical symptomatology of anxiety and family function in adolescents—The self-esteem mediator. Children 11:338. doi: 10.3390/children11030338
Firth, J., Marx, W., Dash, S., Carney, R., Teasdale, S. B., Solmi, M., et al. (2019). The effects of dietary improvement on symptoms of depression and anxiety: A meta-analysis of randomized controlled trials. Biopsychosoc. Sci. Med. 81:265. doi: 10.1097/PSY.0000000000000673
Gevins, A., Smith, M. E., McEvoy, L., and Yu, D. (1997). High-resolution EEG mapping of cortical activation related to working memory: Effects of task difficulty, type of processing, and practice. Cereb. Cortex 7, 374–385. doi: 10.1093/cercor/7.4.374
Gillen, J. B., Little, J. P., Punthakee, Z., Tarnopolsky, M. A., Riddell, M. C., and Gibala, M. J. (2012). Acute high-intensity interval exercise reduces the postprandial glucose response and prevalence of hyperglycaemia in patients with type 2 diabetes. Diabetes Obesity Metabolism 14, 575–577. doi: 10.1111/j.1463-1326.2012.01564.x
Gong, P., Wang, P., Zhou, Y., Wen, X., and Zhang, D. (2024). TFAC-Net: A temporal-frequential attentional convolutional network for driver drowsiness recognition with single-channel EEG. IEEE Trans. Intelligent Transportation Syst. 25, 7004–7016. doi: 10.1109/TITS.2023.3347075
Gotlib, I. H. (1998). EEG alpha asymmetry, depression, and cognitive functioning. Cogn. Emot. 12, 449–478. doi: 10.1080/026999398379673
Guo, Y., Chu, T., Li, Q., Gai, Q., Ma, H., Shi, Y., et al. (2025). Diagnosis of major depressive disorder based on individualized brain functional and structural connectivity. J. Magnetic Resonance Imaging 61, 1712–1725. doi: 10.1002/jmri.29617
Hack, R. L., Aigner, M., Musalek, M., Crevenna, R., and Konicar, L. (2024). Brain regulation training improves emotional competences in patients with alcohol use disorder. Soc. Cogn. Affect. Neurosci. 19:nsae048. doi: 10.1093/scan/nsae048
He, Y., Yu, Q., Yang, T., Zhang, Y., Zhang, K., Jin, X., et al. (2021). Abnormalities in electroencephalographic microstates among adolescents with first episode major depressive disorder. Front. Psychiatry 12:775156. doi: 10.3389/fpsyt.2021.775156
Hermans, E. J., Henckens, M. J. A. G., Joëls, M., and Fernández, G. (2014). Dynamic adaptation of large-scale brain networks in response to acute stressors. Trends Neurosci. 37, 304–314. doi: 10.1016/j.tins.2014.03.006
Hodnik, T., Roytman, S., Bohnen, N. I., and Marusic, U. (2024). Beta–gamma phase-amplitude coupling as a non-invasive biomarker for Parkinson’s disease: Insights from electroencephalography studies. Life 14:391. doi: 10.3390/life14030391
Hsueh, J.-J., Chen, T.-S., Chen, J.-J., and Shaw, F.-Z. (2016). Neurofeedback training of EEG alpha rhythm enhances episodic and working memory. Hum. Brain Mapp. 37, 2662–2675. doi: 10.1002/hbm.23201
Huang, X., Zhang, Y., Wu, X., Jiang, Y., Cai, H., Deng, Y., et al. (2023). A cross-sectional study: Family communication, anxiety, and depression in adolescents: The mediating role of family violence and problematic internet use. BMC Public Health 23:1747. doi: 10.1186/s12889-023-16637-0
Huang, Y., Yi, Y., Chen, Q., Li, H., Feng, S., Zhou, S., et al. (2023). Analysis of EEG features and study of automatic classification in first-episode and drug-naïve patients with major depressive disorder. BMC Psychiatry 23:832. doi: 10.1186/s12888-023-05349-9
Ishii, R., Shinosaki, K., Ukai, S., Inouye, T., Ishihara, T., Yoshimine, T., et al. (1999). Medial prefrontal cortex generates frontal midline theta rhythm. NeuroReport 10, 675–679. doi: 10.1097/00001756-199903170-00003
Jacka, F. N., O’Neil, A., Opie, R., Itsiopoulos, C., Cotton, S., Mohebbi, M., et al. (2017). A randomised controlled trial of dietary improvement for adults with major depression (the ‘SMILES’ trial). BMC Med. 15:23. doi: 10.1186/s12916-017-0791-y
Jones, L., and Wheat, J. (2023). Green and pleasant lands: The affective and cerebral hemodynamic effects of presence in virtual environments during exercise. Percept. Mot. Skills 130, 826–843. doi: 10.1177/00315125221146614
Kappenman, E. S., and Luck, S. J. (2010). The effects of electrode impedance on data quality and statistical significance in ERP recordings. Psychophysiology 47, 888–904. doi: 10.1111/j.1469-8986.2010.01009.x
Keller, N., Whittle, R. S., McHenry, N., Johnston, A., Duncan, C., Ploutz-Snyder, L., et al. (2022). Virtual Reality “exergames”: A promising countermeasure to improve motivation and restorative effects during long duration spaceflight missions. Front. Physiol. 13:932425. doi: 10.3389/fphys.2022.932425
Kessler, R. C., Avenevoli, S., and Ries Merikangas, K. (2001). Mood disorders in children and adolescents: An epidemiologic perspective. Biol. Psychiatry 49, 1002–1014. doi: 10.1016/S0006-3223(01)01129-5
Kim, D., and Ko, Y. J. (2019). The impact of virtual reality (VR) technology on sport spectators’ flow experience and satisfaction. Comput. Hum. Behav. 93, 346–356. doi: 10.1016/j.chb.2018.12.040
Koechlin, E., Ody, C., and Kouneiher, F. (2003). The architecture of cognitive control in the human prefrontal cortex. Science 302, 1181–1185. doi: 10.1126/science.1088545
Lavoie, R., Main, K., King, C., and King, D. (2021). Virtual experience, real consequences: The potential negative emotional consequences of virtual reality gameplay. Virtual Reality 25, 69–81. doi: 10.1007/s10055-020-00440-y
Li, H., Dong, W., Wang, Z., Chen, N., Wu, J., Wang, G., et al. (2021). Effect of a virtual reality-based restorative environment on the emotional and cognitive recovery of individuals with mild-to-moderate anxiety and depression. Int. J. Environ. Res. Public Health 18:9053. doi: 10.3390/ijerph18179053
Li, Y., Wang, L., Han, Q., Han, Q., Jiang, L., Wu, Y., et al. (2023). Preoperative resting-state microstate as a marker for chronic pain after breast cancer surgery. Brain Behav. 13:e3196. doi: 10.1002/brb3.3196
Liu, H., Lou, T., Zhang, Y., Wu, Y., Xiao, Y., Jensen, C. S., et al. (2024). EEG-based multimodal emotion recognition: A machine learning perspective. IEEE Trans. Instrumentation Meas. 73, 1–29. doi: 10.1109/TIM.2024.3369130
Liu, N., Deng, J., Lu, F., and Xiao, J. (2025). Virtual reality enhanced mindfulness and yoga intervention for postpartum depression and anxiety in the post COVID era. Sci. Rep. 15:11766. doi: 10.1038/s41598-025-96165-6
Lovato, N., and Gradisar, M. (2014). A meta-analysis and model of the relationship between sleep and depression in adolescents: Recommendations for future research and clinical practice. Sleep Med. Rev. 18, 521–529. doi: 10.1016/j.smrv.2014.03.006
Luo, Y., Shen, Y., and Fan, X. (2025). EEG microstates in adolescent depression: Effects of depression severity and overall symptoms. J. Affect. Disord. 390:119819. doi: 10.1016/j.jad.2025.119819
MacDonald, A. W., Cohen, J. D., Stenger, V. A., and Carter, C. S. (2000). Dissociating the role of the dorsolateral prefrontal and anterior cingulate cortex in cognitive control. Science 288, 1835–1838. doi: 10.1126/science.288.5472.1835
Menon, V. (2011). Large-scale brain networks and psychopathology: A unifying triple network model. Trends Cogn. Sci. 15, 483–506. doi: 10.1016/j.tics.2011.08.003
Menon, V., and Uddin, L. Q. (2010). Saliency, switching, attention and control: A network model of insula function. Brain Struct. Funct. 214, 655–667. doi: 10.1007/s00429-010-0262-0
Mills, E. J., Chan, A.-W., Wu, P., Vail, A., Guyatt, G. H., and Altman, D. G. (2009). Design, analysis, and presentation of crossover trials. Trials 10:27. doi: 10.1186/1745-6215-10-27
Mojtabai, R., Olfson, M., and Han, B. (2016). National trends in the prevalence and treatment of depression in adolescents and young adults. Pediatrics 138:e20161878. doi: 10.1542/peds.2016-1878
Moullec, Y., Saint-Aubert, J., Manson, J., Cogné, M., and Lécuyer, A. (2022). Multi-sensory display of self-avatar’s physiological state: Virtual breathing and heart beating can increase sensation of effort in VR. IEEE Trans. Visualization Comput. Graphics 28, 3596–3606. doi: 10.1109/TVCG.2022.3203120
Muñoz, J., Mehrabi, S., Li, Y., Basharat, A., Middleton, L. E., Cao, S., et al. (2022). Immersive virtual reality exergames for persons living with dementia: User-centered design study as a multistakeholder team during the COVID-19 pandemic. JMIR Serious Games 10:e29987. doi: 10.2196/29987
Muthelo, L., Mbombi, M. O., Mphekgwana, P., Mabila, L. N., Dhau, I., Tlouyamma, J., et al. (2024). Exploring mental health problems and support needs among pregnant and parenting teenagers in rural areas Of Limpopo, South Africa. BMC Women’s Health 24:236. doi: 10.1186/s12905-024-03040-z
Nelson, M. C., Story, M., Larson, N. I., Neumark-Sztainer, D., and Lytle, L. A. (2008). Emerging adulthood and college-aged youth: An overlooked age for weight-related behavior change. Obesity 16, 2205–2211. doi: 10.1038/oby.2008.365
Ng, H.-Y. H., Wu, C. W., Hsu, H.-C., Huang, C.-M., Hsu, A.-L., Chao, Y.-P., et al. (2024). Neurological evidence of diverse self-help breathing training with virtual reality and biofeedback assistance: Extensive exploration study of electroencephalography markers. JMIR Formative Res. 8:e55478. doi: 10.2196/55478
Nishida, K., Minami, S., Yamane, T., Ueda, S., Tsukuda, B., Ikeda, S., et al. (2025). A single session of tDCS stimulation can modulate an EEG microstate associated with anxiety in patients with depression. Brain Behav. 15:e70580. doi: 10.1002/brb3.70580
Petersen, G. B., Petkakis, G., and Makransky, G. (2022). A study of how immersion and interactivity drive VR learning. Comput. Educ. 179:104429. doi: 10.1016/j.compedu.2021.104429
Roberts, A. C., and Mulvihill, K. G. (2024). Multiple faces of anxiety: A frontal lobe perspective. Trends Neurosci. 47, 708–721. doi: 10.1016/j.tins.2024.07.001
Ruiz, R. J., Simão, R., Saccomani, M. G., Casonatto, J., Alexander, J. L., Rhea, M., et al. (2011). Isolated and combined effects of aerobic and strength exercise on post-exercise blood pressure and cardiac vagal reactivation in normotensive men. J. Strength Conditioning Res. 25:640. doi: 10.1519/JSC.0b013e3181c1fcc7
Samsudin, N., Bailey, R. P., Ries, F., Hashim, S. N. A. B., and Fernandez, J. A. (2024). Assessing the impact of physical activity on reducing depressive symptoms: A rapid review. BMC Sports Sci. Med. Rehabil. 16:107. doi: 10.1186/s13102-024-00895-5
Seeburger, D. T., Xu, N., Ma, M., Larson, S., Godwin, C., Keilholz, S. D., et al. (2024). Time-varying functional connectivity predicts fluctuations in sustained attention in a serial tapping task. Cogn. Affect. Behav. Neurosci. 24, 111–125. doi: 10.3758/s13415-024-01156-1
Silber, M. H., Ancoli-Israel, S., Bonnet, M. H., Chokroverty, S., Grigg-Damberger Madeleine, M., Hirshkowitz, M., et al. (2007). The visual scoring of sleep in adults. J. Clin. Sleep Med. 3, 121–131. doi: 10.5664/jcsm.26814
Spitzer, B., and Haegens, S. (2017). Beyond the status Quo: A role for beta oscillations in endogenous content (Re)activation. eNeuro 4:ENEURO.0170-17.2017. doi: 10.1523/ENEURO.0170-17.2017
Sridharan, D., Levitin, D. J., and Menon, V. (2008). A critical role for the right fronto-insular cortex in switching between central-executive and default-mode networks. Proc. Natl. Acad. Sci. 105, 12569–12574. doi: 10.1073/pnas.0800005105
Tan, H., Wade, C., and Brown, P. (2016). Post-movement beta activity in sensorimotor cortex indexes confidence in the estimations from internal models. J. Neurosci. 36, 1516–1528. doi: 10.1523/JNEUROSCI.3204-15.2016
Tao, M., Huang, H., Gao, J., Cao, Y., and Zhuang, J. (2024). The effects of combining visual-auditory stimuli with exercise on short-term affect improvement: A randomized controlled trial. Sci. Rep. 14:18995. doi: 10.1038/s41598-024-69578-y
Thapa, N., Park, H.-J., Yang, J.-G., Park, K. W., and Park, H. (2023). The effect of virtual reality and physical exercise training on cognitive function in older adults with mild cognitive impairment. Alzheimer’s Dement. 19:e066305. doi: 10.1002/alz.066305
Thapar, A., Collishaw, S., Pine, D. S., and Thapar, A. K. (2012). Depression in adolescence. Lancet 379, 1056–1067. doi: 10.1016/S0140-6736(11)60871-4
Urbańska-Grosz, J., Sitek, E. J., Pakalska, A., Pietraszczyk-Kędziora, B., Skwarska, K., and Walkiewicz, M. (2024). Family functioning, maternal depression, and adolescent cognitive flexibility and its associations with adolescent depression: A cross-sectional study. Children 11:131. doi: 10.3390/children11010131
Wagler, A., and Hanus, M. D. (2018). Comparing virtual reality tourism to real-life experience: Effects of presence and engagement on attitude and enjoyment. Commun. Res. Rep. 35, 456–464. doi: 10.1080/08824096.2018.1525350
Wang, J., Yang, Z., Klugah-Brown, B., Zhang, T., Yang, J., Yuan, J., et al. (2024). The critical mediating roles of the middle temporal gyrus and ventrolateral prefrontal cortex in the dynamic processing of interpersonal emotion regulation. NeuroImage 300:120789. doi: 10.1016/j.neuroimage.2024.120789
Wang, Y., Duan, B., Chen, X., Song, Y., and Liu, X. (2025a). The application of metaverse in mental health. Front. Public Health 13:1463494. doi: 10.3389/fpubh.2025.1463494
Wang, Y., Zhang, F., Gao, Z., Zhang, Z., and Liu, W. (2025b). Effects of virtual reality-based exercise intervention on college students’ mood states: An 8-week randomized controlled trial. Games Health J. doi: 10.1177/2161783X251360803 Online ahead of print.
Weiss, B., Ginige, I., Shannon, L., Giribaldi, B., Murphy-Beiner, A., Murphy, R., et al. (2024). Personality change in a trial of psilocybin therapy v. escitalopram treatment for depression. Psychol. Med. 54, 178–192. doi: 10.1017/S0033291723001514
Whitehouse, A. J. O., Durkin, K., Jaquet, E., and Ziatas, K. (2009). Friendship, loneliness and depression in adolescents with Asperger’s Syndrome. J. Adolesc. 32, 309–322. doi: 10.1016/j.adolescence.2008.03.004
Xie, X.-M., Sha, S., Cai, H., Liu, X., Jiang, I., Zhang, L., et al. (2024). Resting-state alpha activity in the frontal and occipital lobes and assessment of cognitive impairment in depression patients. PRBM 17, 2995–3003. doi: 10.2147/PRBM.S459954
Yin, C., Hsueh, Y.-H., Yeh, C.-Y., Lo, H.-C., and Lan, Y.-T. (2016). A virtual reality-cycling training system for lower limb balance improvement. BioMed Res. Int. 2016:9276508. doi: 10.1155/2016/9276508
Zandbagleh, A., Sanei, S., and Azami, H. (2024). Implications of aperiodic and periodic EEG components in classification of major depressive disorder from source and electrode perspectives. Sensors 24:6103. doi: 10.3390/s24186103
Zhang, J., Wang, J., Wu, Q., Kuang, W., Huang, X., He, Y., et al. (2011). Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder. Biol. Psychiatry 70, 334–342. doi: 10.1016/j.biopsych.2011.05.018
Keywords: adolescents, aerobic exercise, depressive symptoms, electroencephalography, virtual reality
Citation: Yao S, Wang G, Peng T, Zhang L, Ma F and Chi P (2026) Brain neural mechanisms underlying VR enhanced aerobic exercise for mood enhancement in depressed adolescents. Front. Hum. Neurosci. 19:1706793. doi: 10.3389/fnhum.2025.1706793
Received: 16 September 2025; Revised: 14 December 2025; Accepted: 18 December 2025;
Published: 12 January 2026.
Edited by:
Keying Zhang, Southeast University, ChinaCopyright © 2026 Yao, Wang, Peng, Zhang, Ma and Chi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Puyan Chi, Y2hpcDJAY292ZW50cnkuYWMudWs=