
==== Front
Cereb Cortex
Cereb Cortex
cercor
Cerebral Cortex (New York, NY)
1047-3211
1460-2199
Oxford University Press

10.1093/cercor/bhae330
bhae330
Original Article
AcademicSubjects/MED00310
AcademicSubjects/MED00385
AcademicSubjects/SCI01870
Using mindfulness-based intervention to promote executive function in young children: a multivariable and multiscale sample entropy study
https://orcid.org/0000-0001-5160-6799
Xie Sha Faculty of Education, Shenzhen University, Shenzhen 518055, China

https://orcid.org/0009-0002-9985-9968
Lu Shuqi School of Biomedical Engineering, Shenzhen University Medical School, Shenzhen University, Lihu Campus, Shenzhen 518055, China

Lu Jiahao School of Biomedical Engineering, Shenzhen University Medical School, Shenzhen University, Lihu Campus, Shenzhen 518055, China

Gong Chaohui School of Biomedical Engineering, Shenzhen University Medical School, Shenzhen University, Lihu Campus, Shenzhen 518055, China

https://orcid.org/0000-0003-1172-6491
Chang Chunqi School of Biomedical Engineering, Shenzhen University Medical School, Shenzhen University, Lihu Campus, Shenzhen 518055, China
Peng Cheng Laboratory, Shenzhen 518055, China

Corresponding author: Chunqi Chang, School of Biomedical Engineering, Health Science Center, Shenzhen University, Lihu Campus, Shenzhen 518055, China. Email: cqchang@szu.edu.cn
Sha Xie and Shuqi Lu are the co-first authors.

9 2024
05 9 2024
05 9 2024
34 9 bhae33013 3 2024
09 7 2024
10 7 2024
© The Author(s) 2024. Published by Oxford University Press.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Early childhood marks a pivotal period in the maturation of executive function, the cognitive ability to consciously regulate actions and thoughts. Mindfulness-based interventions have shown promise in bolstering executive function in children. This study used the functional near-infrared spectroscopy technique to explore the impact of mindfulness-based training on young children. Brain imaging data were collected from 68 children (41 boys, aged 61.8 ± 10.7 months) who were randomly assigned to either an intervention group (N = 37, aged 60.03 ± 11.14 months) or a control group (N = 31, aged 59.99 ± 10.89 months). Multivariate and multiscale sample entropy analyses were used. The results showed that: (1) brain complexity was reduced in the intervention group after receiving the mindfulness-based intervention in all three executive function tasks (ps < 0.05), indicating a more efficient neural processing mechanism after the intervention; (2) difference comparisons between the intervention and control groups showed significant differences in relevant brain regions during cognitive shifting (left dorsolateral prefrontal cortex and medial prefrontal cortex) and working memory tasks (left dorsolateral prefrontal cortex), which corroborates with improved behavioral results in the intervention group (Z = −3.674, P < 0.001 for cognitive shifting; Z = 2.594, P < 0.01 for working memory). These findings improve our understanding of early brain development in young children and highlight the neural mechanisms by which mindfulness-based interventions affect executive function. Implications for early intervention to promote young children’s brain development are also addressed.

brain complexity
multivariable and multiscale sample entropy
executive function
mindfulness-based intervention
National Natural Science Foundation of China 10.13039/501100001809 61971289 Guangdong Planning Office of Philosophy and Social Science, China GD21YJY09 Shenzhen Education Science Planning Project, China zdzz21002 Shenzhen-Hong Kong Institute of Brain Science-Shenzhen Fundamental Research Institutions 2023SHIBS0003 SZU Top Ranking Project 86000000210
==== Body
pmcIntroduction

Executive function (EF) encapsulates the cognitive prowess essential for the deliberate regulation of one’s actions and thoughts, representing a sophisticated mental process (Zelazo and Müller 2002; Griffin et al. 2016). Recognized as integral to children’s preparatory cognitive abilities for formal education (Blair 2002; Guedes et al. 2023), early childhood EF undergoes pivotal developmental stages, playing a foundational role in behavior regulation (Moriguchi and Hiraki 2011). The cultivation of EF during this critical period lays the groundwork for subsequent advancements in children’s social, emotional, and cognitive competencies (Griffin et al. 2016).

Empirical evidence underscores the potential of mindfulness training in significantly enhancing EF among young children and simultaneously mitigating problematic behavioral manifestations (Flook et al. 2015; Razza et al. 2015). Despite this promising evidence, existing studies often provide a limited examination of the tripartite components of EF—working memory (WM), cognitive flexibility (CF), and inhibitory control (IC). Moreover, there is a paucity of research that substantiates the efficacy of these interventions through neuroimaging evidence, leaving a critical gap in understanding the neural mechanisms underlying these improvements.

To address this gap, our study aims to provide a comprehensive exploration of the effects of a mindfulness-based intervention on EF in young children. By employing multivariable and multiscale sample entropy analysis, we seek to elucidate the changes in brain complexity associated with enhanced EF. This approach not only advances the current literature by integrating detailed neuroimaging evidence but also offers novel insights into how mindfulness training can be effectively utilized to promote cognitive development in early childhood.

EF and its three components

EF is commonly divided into three cores: IC (covering self-control, such as behavioral inhibition, and interference control, such as selective attention and cognitive inhibition), WM (ability to temporarily store and manipulate information), and CF (the ability to flexibly shift the focus of one’s mental frame) (Blair 2016; Moriguchi 2017). Numerous studies have noted significant activation of the prefrontal cortex of the brain when performing EF-related cognitive tasks (Huttenlocher and Dabholkar 1997; Scheibel and Levin 1997; Anderson 2008). During infancy and preschool, these core components of executive functioning develop, forming a critical foundation that sets the stage for the development of higher cognitive processes in adulthood. A series of studies in structural equation modeling have shown that each subdomain of EF plays a different role in performance on a range of executive outcome metrics, underscoring the importance of recognizing the diversity of these sub-processes (Garon et al. 2008). Functional near-infrared spectroscopy (fNIRS) measures task-related changes in cerebral hemodynamics. Because of its portability, fNIRS is widely used in the research of EF. Therefore, the present study used the fNIRS technique to measure preschoolers’ task performance in three different EF cores: IC (Cragg and Nation 2008; Inoue et al. 2012), CF (Davidson et al. 2006), and WM (Baddeley and Hitch 1994; Smith and Jonides 1999).

Advances in studying the EF’s three components

CF, one of the primary EFs, typically begins to develop at approximately 3–4 years of age. During the first 2 years of life, early CF skills are primarily supported by attentional control and language skills, both of which are critical to the early development of goal-directed behavior (Hendry et al. 2016). Research suggests that IC and WM underlie CF development, beginning at approximately 3 years of age, and that they each contribute to successful switching behaviors in different ways (Blakey et al. 2016). In particular, IC promotes CF by allowing children to ignore task-irrelevant information, while WM functions by allowing children to maintain and update new sequencing rules. Recent research has noted that cognitive performance is highly dependent on the task and context employed (Perry et al. 2014; Yee and Thompson-Schill 2016).

The Dimension-Converting Card Sorting Test (DCCS) is a classic task for assessing CF. Findings show that when three-year-olds are asked to sort complex stimuli, they typically sort according to an initial rule, whether that rule is based on shape or color. However, when asked to switch to a new rule, most three-year-olds failed to adapt and continued to sort according to the initial rule (Zelazo 2006). In other words, they tended to stick to their initial response pattern when confronted with a request for a rule change. In contrast, children around the age of four usually can sort stimuli according to the new rules. Thus, research on the DCCS task tends to suggest that the critical period of CF development is around three to four years of age, with the most significant developmental advances manifested in the ability to overcome perseveration (Munakata et al. 2012). Significant neural activity occurred in the bilateral inferior prefrontal region during the conduct of the experiment. Thus, developmental psychologists believe that maturation of the lateral prefrontal cortex (LPFC) may play a key role in the development of cognitive switching abilities (Davidson et al. 2006; Gupta et al. 2009).

During preschool years, children’s ability to control potentially disruptive thought processes and behaviors, namely IC ability, also develops rapidly. Three-year-olds typically have difficulty performing tasks that require IC of attentional and motor responses, such as inhibiting dominant responses to follow rules. However, by the age of five, their performance on these tasks improves significantly (Carlson 2016). Indeed, IC is thought to make an important contribution to individual differences and developmental changes in various cognitive abilities, such as attention, memory, reading comprehension, and theory of mind (Carlson et al. 2004). The Go/No-Go task is widely recognized as an important measure of IC, and the research of Moffitt et al. (2011) has shown that children from the ages of 3–11 years perform better in terms of IC. However, with the normal aging process, people’s IC declines significantly (Hasher and Zacks 1988; Hasher et al. 1991). In the field of neuroimaging, studies on response inhibition tasks have shown that such tasks persistently and specifically activate the right inferior frontal cortex. In addition, analyses of functional connectivity have found that short-range connections between the right frontal and right parietal cortex show stronger partial coherence in 4- to 6-year-old children than in adults (Mehnert et al. 2013). These findings further emphasize the uniqueness and developmental dynamics of children’s brain networks during IC.

The concept of WM originates from an in-depth understanding of short-term memory, which is primarily responsible for the temporary storage of a limited amount of information over a short period. WM is closely linked to IC, and the two are often interdependent and work together. WM declines progressively with age and gradually improves during early development. There is a close correlation between this change in WM and the decline in processing speed with age and its enhancement during early development. (Fiore et al. 2012; Fournet et al. 2012). When performing WM tasks, the frontal and parietal cortical regions show significant activation (Tsujimoto et al. 2004; Buss et al. 2014). Furthermore, there is a positive correlation between children’s age and activity, response accuracy, and speed in the LPFC (Perlman et al. 2016). These findings highlight the developmental dynamics of children’s WM capacity and its neurophysiological basis in brain networks.

Advances in early intervention studies

EF can be improved with a targeted training (Klingberg 2010; Diamond and Lee 2011). Research has demonstrated positive outcomes from various interventions, including aerobic exercise (Bryant 2019), mindfulness (Farb et al. 2007), yoga (Gothe et al. 2014), and the Tools of the Mind early childhood program (Manjunath and Telles 2001). Notably, mindfulness-based training has been widely implemented in early childhood classrooms (Zelazo et al. 2018; Li et al. 2019; Xie et al. 2022), emphasizing the potential of mindfulness to enhance EF by reducing mind-wandering and improving focus (Hölzel et al. 2011).

Mindfulness interventions incorporate key components such as physical relaxation, breathing exercises, mental imagery, and mind–body awareness to foster an integrated mind–body connection (Thera 1972; Kabat-Zinn and Hanh, 2009; Tang et al. 2015). These interventions have shown positive effects on various components of EF, including attentional control, emotion regulation, and stress response (Tang, 2005, 2007, Tang et al. 2012; Tang and Posner 2009). Importantly, these effects are attributed specifically to the mindfulness components rather than mere contact with a trained interventionist.

Neuroimaging studies have demonstrated increased activation in the right dorsolateral and ventral lateral prefrontal cortex (PFC) following mindfulness training (Farb et al. 2007; Tang et al. 2015). Additionally, higher levels of mindfulness correlate with greater PFC activation (Xie et al. 2022; Lorenzetti et al. 2023; Yue et al. 2023), and experienced meditators exhibit enhanced PFC activity during various tasks (Hölzel et al. 2007). Functional magnetic resonance imaging (fMRI) studies further reveal that mindfulness training enhances performance on attentional control (Tang et al. 2015), IC (Jack et al. 2013), and WM tasks (Mrazek et al. 2013).

Despite these promising findings, a significant gap remains in the literature. Many studies have focused on the behavioral outcomes of mindfulness interventions without thoroughly investigating the underlying neural mechanisms, especially in young children. Moreover, while fMRI studies provide valuable insights, they are often limited by their constraints and less feasible with preschoolers. fNIRS offers a more portable and child-friendly alternative for neuroimaging, yet its application in mindfulness intervention studies remains underexplored.

Our study aims to address this gap by using fNIRS to examine the neural correlates of improved EF following mindfulness-based intervention in young children. By integrating multivariable and multiscale sample entropy analysis, we seek to provide a comprehensive understanding of how mindfulness training influences brain complexity and EF, offering novel insights that extend beyond current literature.

Multivariate multiscale sample entropy

Entropy is a powerful tool for characterizing the regularity or complexity of a system (Elbert et al. 1994) and can be used to analyze systems at a single time scale (Costa et al. 2005). In the context of neuroimaging, Sample Entropy (SampEn) has been recognized as a biomarker for Alzheimer’s disease diagnosis when applied to fNIRS signaling in the frontal cortex (Li et al. 2018; Perpetuini et al. 2018, 2019). Additionally, SampEn has been used to study the effects of WM tasks on fNIRS signals (Angsuwatanakul et al. 2015).

Building on this, Ahmed and Mandic (2011) extended the concept of Multivariate Sample Entropy (MSE) by integrating it with the principles of multiscale entropy, resulting in Multivariate Multiscale Sample Entropy (MMSE). This innovative approach allows for the assessment of structural complexity across multiple time scales and variables, making it particularly suitable for analyzing the nonlinear dynamic behavior of the human brain (Zhu et al. 2020). MMSE has been successfully applied to various multichannel biomedical signals, such as electroencephalograms and fNIRS, providing deeper insights into the complex dynamics of these data (Ding et al. 2023).

Despite the demonstrated potential of MMSE in analyzing fNIRS signals, its application to the study of EF in children remains underdeveloped. Most existing research has focused on its use in diagnosing neurological conditions or assessing cognitive tasks in adults. The lack of studies employing MMSE to investigate the neural mechanisms underlying EF improvements in young children following mindfulness-based interventions represents a significant gap in the literature.

Our study aims to fill this gap by utilizing MMSE to analyze fNIRS data collected from young children undergoing mindfulness-based training. This approach will enable us to explore the complex neural dynamics associated with EF and assess the impact of mindfulness interventions on brain complexity. By doing so, we aim to provide novel insights into the neural underpinnings of EF development in early childhood, thereby advancing the field and addressing a critical gap in the current research landscape.

The objective of this study

In the present study, we aimed to explore the effects of mindfulness-based training on children’s executive functioning and the changes in its neural mechanisms. Specifically, we administered a 15-session, twice per week mindfulness-based intervention to a group of 3- to 6-year-olds and used fNIRS to measure brain activity before and after the intervention, alongside MMSE analysis. This approach addresses the identified literature gap regarding the application of MMSE in studying EF in young children and the need for neuroimaging evidence to support the benefits of mindfulness interventions.

Based on this, this study proposed and examined the following two hypotheses:

Hypothesis 1: Children participating in the mindfulness-based intervention will demonstrate significant differences in executive functioning, as measured by the level of brain area activation, compared to pre-intervention baseline levels.

Hypothesis 2: After the mindfulness-based intervention, there will be a significant difference between the control and intervention groups in terms of accuracy and response time in the behavioral tasks and the level of activation of brain regions related to executive functioning, suggesting greater enhancement in the intervention group.

Materials and methods

Participants

This study recruited 112 typically developed preschool children. The parents of all the participating children provided written consent and were informed of the purpose of the study and the safety of the fNIRS experiment. The experiments were approved by the university ethics committee (PN-2021-038). Of these, 44 preschool children were excluded from the final analysis because of front-end or back-end fNIRS data tuning, poor data collection quality, or refusal to wear the cap. Ultimately, 68 children (Mage = 60.00 months, SD = 10.94) were included in the MMSE analysis. There was no statistically significant difference in the age of the participants between the intervention (n = 37, Mage = 60.03 months, SD = 11.14) and control (n = 31, Mage = 60.00 months, SD = 10.89) groups (t = 0.22, P = 0.98). The intervention and data collection occurred from October 2022 to June 2023.

Measures

This study employed three tasks to assess EFs in children. The first task, Dimensional Change Card Sort (DCCS), is commonly used to study cognitive shifting in children. The second task, Go/No-Go, was used to measure IC in children. The third behavioral task was used to evaluate WM in children.

Stimuli were displayed on a computer screen, and the participants used a keyboard to respond. Before starting each task, the participants underwent training sessions to ensure that they understood the rules of the task. During the testing phase, the experimenter recorded each participant’s response time and actions, with participants clicking on the keyboard to make the corresponding responses. In the rest phases, participants were asked to look at a “+” on the computer screen and remain still.

DCCS task

The DCCS task is commonly used to assess children’s CF, also known as task-switching or set-shifting. This experiment, designed by Zelazo’s (2006) team, has been widely applied in studying the development of EF during early childhood. In the traditional DCCS task, two target cards are displayed on a computer screen, such as a blue rabbit and a red boat, and participants are required to categorize a series of cards that have two attributes. These stimulus cards possess two dimensions: color and shape. The target cards (like a red boat and a blue rabbit) match the test cards (such as a blue boat or a red rabbit) on one dimension but not on the other.

Based on Moriguchi et al’.s (Moriguchi et al. 2007; Moriguchi and Itakura 2008) modification and our previous study (Xie et al. 2022), the present study further modified the DCCS task by introducing a mixed sorting card task to increase the task difficulty and a blank card task to exclude the effects from keystrokes (Fig. 1). In the mixed sorting card task phase, participants were asked to sort test cards based on verbal instructions (shape or color). In the blank card task phase, participants were asked to randomize key presses to eliminate the effects of key press actions on the brain. The experiment consisted of four consecutive test blocks interspersed with four rest phases. During the rest phases, a “+” symbol appeared in the center of the screen and participants were instructed to focus on the symbol and remain still. The first three test blocks were mixed blocks containing eight randomized trials. The fourth test block was a blank task lasting 50 s. Participants responded by pressing keys on the keyboard, and the computer recorded the answers and reaction times for the tests. The total number of correct answers in each block was calculated as an indicator of overall performance for subsequent analysis.

Fig. 1 Task paradigm of the DCCS task.

Go/No-Go Task

The Go/No-Go task, modified from the paradigm of Lahat et al. (2010) and our previous study (Xie et al. 2022) by replacing the fixation picture with the same picture as the screen background to reduce interference, aims to measure children’s IC. This task was chosen for its proven effectiveness and solid neural basis mapping (Wiebe et al. 2012). In each trial (Fig. 2), images of animals (such as cows, horses, tigers, or dogs) appear on the screen. The participants were instructed to press the space bar on the keyboard immediately upon seeing any animal except for a dog (Go stimulus) and not to press when a dog image appeared (No-Go stimulus). During the practice phase, there were four Go trials and four No-Go trials, with participants receiving reminders when they made incorrect responses. The entire task was divided into three blocks, each containing 16 Go trials and 4 No-Go trials, which were randomly assigned. Participants had resting periods after completing each block of trials. For subsequent analysis, accuracy was calculated by dividing the number of correct responses by the total number of trials.

Fig. 2 Task paradigm of the Go/No-Go task.

WM task

The third behavioral task is an adaptation from a previous study (Perlman et al. 2016) designed to measure the WM capacity of preschool children aged 3–6 years. In this experiment (Fig. 3), each trial displayed six trees on the screen, with a monkey randomly appearing on one of the six trees. Participants were required to remember the location where the monkey appeared. After 2 s, the monkey disappeared, leaving the image of the six trees on display for another 2 s. Subsequently, a question mark appeared at the location of each tree on the screen, and the participants had 4 s to indicate the monkey’s previous position. The experimenter recorded the answer by pressing the corresponding key on the basis of the participant’s indication. During the practice phase, participants received feedback that indicated whether they responded correctly or not. The entire task was divided into three blocks, each containing eight trials. Between each trial, an unrelated image appeared for 1 s to remind participants to prepare for the next trial. There was a 20-s rest period for participants between the three blocks of the task. For subsequent analysis, accuracy was calculated by dividing the number of times the correct location was indicated by the total number of trials.

Fig. 3 Task paradigm of the WM task.

fNIRS data recording and processing

During cognitive tasks, brain regions that are more active require more oxygen, leading to an increase in oxygenated hemoglobin and a decrease in deoxygenated hemoglobin in these areas. fNIRS (Ferrari and Quaresima 2012) is a non-invasive imaging technique used to monitor brain activity. It indicates neural activity by measuring changes in cerebral blood flow, specifically the levels of oxygenated and deoxygenated hemoglobin. In this study, a multi-channel fNIRS system (OxyMon MK III, Artinis, The Netherlands) was used to measure the concentration changes of oxygenated hemoglobin (HbO) and deoxygenated hemoglobin (HbR) in participants during the experiment. Based on EF research in children (Diamond 2013) and the experimental design of the current study, we employed a 3 × 10 light pole template located at the forehead, covering the frontal area with 44 channels. To ensure consistent light pole array placement for all participants, the central-lower position of the array was located at the Fpz position, consistent with the 10-20 measurement system. Accordingly, the regions of interest (ROIs) were the left ventrolateral prefrontal cortex (VLPFC), right VLPFC, left dorsolateral prefrontal cortex (DLPFC), right DLPFC, left posterior superior frontal cortex (PSFC), right PSFC, left temporal cortex (TC), right TC, and medial prefrontal cortex (MPFC) (Fig. 4). Previous studies have shown (Menon and D’Esposito 2022; Buckner and DiNicola 2019) that the frontal regions are actively involved in the maintenance and regulation of EF. The sampling rate was set at 50 Hz for data collection. The differential pathlength factor (DPF) value for each participant was calculated based on their age using the formula (DPF = 4.99 + 0.0678 × Age^0.814).

Fig. 4 Localization of regions of interest.

Procedure

First, the principal of the target preschool was contacted and explained of the study purpose. Then, print questionnaires and consent forms were enveloped and carried home by the children to their parents, who gave consent and filled the questionnaires. The children then brought them back to the preschool, and class teachers collected these forms and passed them to the research team. Next, the children were invited to complete the three EF tasks in a quiet classroom in the preschool. Before the tasks, an experienced NIRS technician placed the child-sized NIRS cap and installed the optodes. At the same time, a research assistant who majored in early childhood education engaged in story-book reading with the child. All three tasks were computerized using Psychophysics Toolbox extensions and displayed on a 52.71 cm × 29.65 cm Huawei monitor. The children were trained to perform the tasks before each experiment began. The order of the three tasks was randomly assigned.

After the pre-intervention data collection, the participants were randomly assigned to the intervention and control groups. This mindfulness intervention was adapted from the MindUP training package to be both playful and mindful for preschoolers. A certified MindUP trainer and a group of researchers in early childhood education adapted the training course. The class teachers of the intervention groups received training from the certified teacher and delivered mindfulness training during the school day. The interventions were structured in three parts: my mindful brain, my five senses, and my mindful world. For the intervention group, the children received 20 min per session, twice a week, all together 15 sessions of mindfulness intervention. After the intervention, all participants from both the intervention and control groups were assessed for their EF again.

Data analysis plan

Behavioral results analysis

In Matlab, we analyzed the behavioral results of the participants in the three experimental tasks. Specifically, for the DCCS and GO/NO-GO tasks, we calculated the response time and accuracy for each participant, while for the WM task, only accuracy was calculated. Then, because the behavioral data did not fit a normal distribution, we used Wilcoxon Signed-Rank Test to compare the differences before and after the mindfulness-based training as well as the differences in the intervention group before and after the intervention in order to explore whether the changes in the brain manifested in the behavioral data. All of our statistics use the false discovery rate correction.

fNIRS data preprocessing

For preprocessing of fNIRS data, we used the Homer2 (Huppert et al. 2009) software package to process the multi-channel near-infrared data. First, we manually removed channels with poor quality or bad contact to ensure data quality, converted raw signals into optical density signals, and applied a wavelet transform algorithm (Molavi and Dumont 2012) to remove motion artifacts. These steps helped reduce signal interference caused by participant movement. Next, we filtered the data using a third-order Butterworth filter to remove frequencies outside the range of 0.01–0.08 Hz. This aided in reducing high-frequency noise and low-frequency drift. We then converted optical density signals into concentration changes of oxygenated and deoxygenated hemoglobin and conducted correlation-based signal improvement to further correct for motion artifacts (Engelmann et al. 2012).

Statistical analysis

We performed MMSE calculations on the preprocessed fNIRS signals to measure brain complexity. MMSE is an extension based on the MSE. MSE is used to calculate the complexity of a single time series at different temporal scales.

The calculation of the MMSE is divided into two steps:

Coarse-Graining of Multivariate Time Series: the multivariate time series \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${\left\{{x}_k,i\right\}}_{i=1}^N$\end{document}, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $k=1,2,\dots, p$\end{document}, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $p$\end{document} represents the number of variables (channels), and\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $N$\end{document} is the number of samples per variable, is coarse-grained to define increasing time scales. For a scale factor\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\tau$\end{document}, the coarse-graining process comes in two forms: non-overlapping, where every \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\tau$\end{document} data points are jumped and averaged to produce a new data point; and overlapping, where the jump is 1 − \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\tau$\end{document} data points, averaging \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\tau$\end{document} data points each time (Busa and Van Emmerik 2016). The mathematical definition of a non-overlapping multivariate coarse-grained time series is as follows:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $$ {y}_{j,k}^{\left(\tau \right)}=\frac{1}{\tau}\sum_{i=\left(j-1\right)\tau +1}^{j\tau}{x}_{i,k,} $$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $1\le j\le \frac{N}{\tau }$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $k=1,\dots, p$\end{document}.

Coarse-grain the original multivariate time series and calculate the MSE for each coarse-grained sequence. Plot the relationship between MSampEn and the scale factor \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\tau$\end{document}.

MSE is the calculation of multiscale entropy simultaneously across multiple channels. For a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $p$\end{document}-variate time series \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${\left\{{x}_k,i\right\}}_{\mathrm{i}=1}^{\mathrm{N}}$\end{document}, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $k=1,2,\dots, p$\end{document}, the mathematical definition of MSE is:

Construct \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\left(N-n\right)$\end{document} composite delay vectors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${X}_m(i)\in{R}^{\mathrm{m}}$\end{document}, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $i=1,2,\dots, N-n$\end{document}, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $n=\mathit{\max}\left\{M\right\}\times \mathit{\max}\left\{\tau \right\}$\end{document}.

Define the distance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $d$\end{document}[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${\mathrm{X}}_m\left(\mathrm{i}\right),{\mathrm{X}}_m\left(\mathrm{j}\right)$\end{document} between any two composite delay vectors \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${\mathrm{X}}_m\left(\mathrm{i}\right)$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${\mathrm{X}}_m\left(\mathrm{j}\right)$\end{document} as the absolute value of the largest difference between their corresponding elements. That is:

(1) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*} \mathrm{d}\left[{X}_m(i),{X}_m(j)\right]={\mathit{\max}}_{k=0,\dots, m-1}\left(\left|x\left(i+k\right)-x\left(j+k\right)\right|\right). \end{equation*}\end{document}

For a given \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${X}_m(i)$\end{document} and threshold \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document}, count the number of instances when the distance between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${X}_m(i)$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${X}_m(j)$\end{document} is less than or equal to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document} (where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\left(1\le j\le N-n,j\ne i\right)$\end{document}), denoted as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${P}_i$\end{document}. Then, calculate the frequency of occurrence:

(2) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*} {B}_i^m(r)=\frac{1}{N-n-1}{P}_i, \end{equation*}\end{document}

And define a global quantity:

(3) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*} {B}^m(r)=\frac{1}{N-n}\sum_{i=1}^{N-n}{B}_i^m(r). \end{equation*}\end{document}

Increase the dimension to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $m+1$\end{document} and calculate the number of instances where the distance between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${X}_{m+1}(i)$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${X}_{m+1}(j)$\end{document} (for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $1\le j\le N-n,j\ne i$\end{document}) is less than or equal to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document}, denoted as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} ${A}_i$\end{document}. Then, compute the frequency of occurrence:

(4) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*} {B}_i^m(r)=\frac{1}{p\left(N-n\right)-1}{Q}_i, \end{equation*}\end{document}

And define a global quantity:

(5) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*} {B}^{m+1}(r)=\frac{1}{p\left(N-n\right)}\sum_{i=1}^{p\left(N-n\right)}{B}_i^{m+1}(r). \end{equation*}\end{document}

Thus, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $ {B}^m(r) $\end{document} is the probability of any two composite delay sequences matching at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $ m $\end{document} points under the similarity tolerance \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $ r $\end{document}, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $ {B}^{m+1}(r) $\end{document} is the probability of any two composite delay sequences matching at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $ m+1 $\end{document} points.

For a tolerance level \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document}, Multivariate Sample Entropy (MSamEn) is calculated as the negative natural logarithm of the conditional probability that two composite delay vectors, which are close to each other in an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $m$\end{document}-dimensional space, remain close when the dimension is increased by 1. The calculation formula is as follows:

(6) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*} M\ SampEn\left(M,\tau, r,N\right)=-\mathit{\ln}\left[\frac{B^{m+1}(r)}{B^m(r)}\right]. \end{equation*}\end{document}

Multivariate Sample Entropy (MSampEn) adopts the algorithm by Cao (Ahmed and Mandic 2011). Regarding the choice of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document}, it is important to remember that to avoid significant contributions from noise to the sample estimation, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document} must be larger than the noise level present in most of the signal. Another criterion for selecting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document} is based on the dynamic characteristics of the signal (signal dynamics) (Ahmed and Mandic 2012). Typically, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $r$\end{document} is generally taken as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $0.15\ast std\left({x}_k\right)$\end{document}. The scale factor \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $\tau$\end{document} is calculated using the autocorrelation function (Cao 1997).

Next, we divided the 44 fNIRS channels into nine ROIs, including the left VLPFC, right VLPFC, left DLPFC, right DLPFC, left PSFC, right PSFC, left TC, right TC, and MPFC. Our aim was to investigate whether there were significant differences in the complexity of these ROIs across different tasks.

All calculations were analyzed by t-test using SPSS software. In the preliminary analysis, we compared the complexity of each ROI before and after mindfulness training using paired t-tests to assess the impact of mindfulness training on brain regions in the intervention group (Hypothesis 1). In processing the fNIRS data from the DCCS task, we implemented an innovative methodology by calculating the difference between the fNIRS data from each task phase and that from a null task (random key presses). This approach effectively eliminated neural activity unrelated to the task, thus more accurately reflecting brain activity changes induced by specific tasks. This finding is significant because it validates the precision of the adopted method in capturing neurophysiological changes triggered by task engagement. We also compared the complexity of each ROI before and after mindfulness training using paired t-tests to assess the impact of mindfulness training on brain regions in the control group. In a more in-depth analysis, we calculated the difference in the MMSE results and behavioral performance of the participants before and after training. We then compared the differences between the intervention and control groups using t-tests to explore the specific effects of positive thinking training on the brain (Hypothesis 2). All of our statistics use the FDR correction。.

Results

The t-test analysis was conducted on the MMSE results of the intervention and control groups before and after the tests. The MMSE comparison results for the DCCS task are shown in Fig. 5. In the DCCS task, significant differences were observed in the intervention group in the left DLPFC (t = 4.713, P < 0.001) (Fig. 5a), right DLPFC (t = 4.964, P < 0.001) (Fig. 5b), right PSFC (t =3.741, P < 0.001) (Fig. 5c), and MPFC (t =3.262, P < 0.01) (Fig. 5d), but no significant changes in the control group.

Fig. 5 Comparison of MMSE before and after the intervention for the DCCS task (full sample).

The MMSE comparison results for the GO/NO-GO task are shown in Fig. 6. In the intervention group during the GO/NO-GO task, significant differences were observed in the left DLPFC (t = 4.348, P < 0.001) (Fig. 6a), right DLPFC (t = 44.477, P < 0.001) (Fig. 6b), and right PSFC (t = 3.079, P < 0.01) (Fig. 6c). Significant differences were also observed in the control group in the right PSFC (t = 2.147, P < 0.05) (Fig. 6c) and right TC (t = 2.132, P < 0.05) (Fig. 6d).

Fig. 6 Comparison of MMSE before and after the intervention for the Go/No-Go task (full sample).

The MMSE comparison results for the WM task are shown in Fig. 7. In the intervention group during the WM task, significant differences were observed in the left DLPFC (t = 4.453, P < 0.001) (Fig. 7a), right DLPFC (t = 2.662, P < 0.05) (Fig. 7b), left PSFC (t = 2.906, P < 0.05) (Fig. 7d), right PSFC (t = 3.125, P < 0.01) (Fig. 7e), and MPFC (t = 2.697, P < 0.01) (Fig. 7f). Significant differences were also observed in the control group in the right DLPFC (t = 2.301, P < 0.05) (Fig. 7b) and right TC (t = 2.140, P < 0.05) (Fig. 7c).

Fig. 7 Comparison of MMSE before and after the intervention for the WM task (full sample).

Notably, in all brain regions where significant differences were observed, there was a trend of decreased complexity (i.e. post-test MMSE results were lower than pre-test), indicating that the complexity of relevant brain regions decreased after mindfulness-based training. Putting together, the findings indicated that mindfulness-based intervention demonstrated significant differences in executive functioning, as measured by the level of brain area activation, compared to pre-intervention baseline levels, supporting H1.

More direct evidence of the positive effect of mindfulness training came from difference comparison in MMSE and behavioral results between the intervention and control groups. The difference comparison results for the DCCS task are presented in Fig. 8. For the DCCS task, significant differences were observed between the intervention and control groups in the left DLPFC (t = − 3.178, P < 0.01) and MPFC (t = 2.232, P < 0.05). These results suggest that mindfulness training had a significant impact on the activity of the left DLPFC and MPFC in the DCCS task. Appendix 1 displayed the results of difference comparison for all ROIs for the DCCS task.

Fig. 8 Pre-test and post-test difference comparison between intervention and control groups for the DCCS task.

The difference comparison results for the GO/NO-GO task showed no significant differences were observed during the GO/NO-GO task (P < 0.05), suggesting that mindfulness training might not have differentially affected brain activity between the intervention and control groups. Appendix 2 displayed the results of difference comparison for all ROIs in the GO/NO-GO task.

The difference comparison results for the WM task are also displayed in Fig. 9. During the task, a significant difference was observed in the left DLPFC (t = − 2.426, P < 0.05), indicating that mindfulness training significantly impacted the activity of the left DLPFC in the WM task. Appendix 3 displayed the results of difference comparison for all ROIs in the WM task.

Fig. 9 Pre-test and post-test difference comparison between intervention and control groups for the WM task.

The results of the different comparisons of the pre- and post-tests of the behavioral data are shown in Fig. 10. In the DCCS task (Fig. 10a), there was a significant difference in correct response time for both the intervention (Z = − 3.674, P < 0.001) and control groups (Z = − 2.841, P < 0.01), but the change was more pronounced in the intervention group. In the GO/NO-GO task (Fig. 10b), there were no differences in either group. This may indicate that although some changes occurred at the neurophysiological level, these changes may not yet be reflected at the behavioral level, or the task design may not have been effective in capturing changes at the behavioral level. In the WM task (Fig. 10(c), the intervention group (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} $Z$\end{document}=2.594, P < 0.01) showed a significant increase in correctness. Putting together, significant differences between the control and intervention groups in terms of the level of activation of brain regions related to executive functioning and behavioral tasks were observed, with greater enhancement in the intervention group, thus supporting H2.

Fig. 10 Comparison of behavioral data between pre- and post-tests.

Discussion

In this study, we found MMSE changes before and after mindfulness training in all three EF tasks (H1) and confirmed the positive effects of mindfulness training on children’s brain complexity and on the promotion of EF in specific brain regions (H2). We will discuss these findings in details in the following section.

Decreased brain complexity after intervention

The current study found significant decreases in MMSE in the LDLPFC, RDLPFC, RPSFC, and MPFC during the DCCS task, decreases in the LDLPFC, RDLPFC, RPSFC, and RTC during the GO/NO-GO task, as well as decreases in the LDLPFC, RDLPFC, LPSFC, RPSFC, and MPFC during the WM task, in the intervention group. From an informational perspective, fNIRS is a non-smooth, non-linear time series, so its entropy can be used to distinguish between different healthy brain states (Ghouse et al. 2020). Based on this, MMSE is proposed based on sample entropy, as a measure of time series complexity, with smaller the entropy, smaller the time series complexity and higher the self-similarity (Wang et al. 2021). In children, the decrease in complexity of brain regions may be related to different stages of brain development. During early childhood, the brain undergoes a phase of rapid development and shaping, which involves the formation of a large number of neurons and synapses, as well as the establishment and remodeling of neural networks (Gilmore et al. 2018; McLaughlin et al. 2019). During this process, brain regions may exhibit relatively complex structures and functions in response to various learning and cognitive tasks. Thus, the high complexity of the brain implies that more cognitive functions are involved (Ding et al. (2023).

However, as children grow, the brain progressively optimizes and consolidates these connections and may have fewer redundant connections or functions in specific brain regions (McLaughlin et al. 2014; Cao et al. 2017). This reduction in brain region complexity may reflect the brain’s streamlining and optimisation during learning and development, allowing it to perform tasks more efficiently. This phenomenon may play an important role in the development of a number of cognitive abilities, such as improved concentration, memory and problem-solving skills. A decrease in brain complexity, as measured by entropy, might reflect this increased efficiency and reduced neural variability. Previous research suggests that higher entropy in neural signals can indicate a more disorganized and less efficient brain state (Li et al. 2018). Therefore, the observed reduction in complexity could signify a shift towards more streamlined and effective neural processing in the DLPFC and PSFC, regions crucial for EF.

Furthermore, MMSE was able to discover specific brain region activity in oxygenated hemoglobin that could not be resolved using mean estimates of any of the three hemoglobin concentrations (HbR, HbO, HbT). This enhanced sensitivity may stem from the additional quantification of nonlinear and complex kinetic properties provided by entropy analysis. Even when linear effects are weak or insignificant, nonlinear and complex behavioral patterns may remain and can be detected (Perpetuini et al. 2020; Zhang et al. 2021). Thus, the MMSE method is more sensitive in identifying changes in brain regions and more effective in revealing the effects of interventions compared with traditional methods. When the entropy of brain activity decreases (complexity decreases), this could mean that the activity in the brain network becomes more sequential or that the way the brain processes information becomes more efficient and focused. Thus, the reduction in complexity might also reflect the optimization of neural resources. As children become more adept at mindfulness practices, their brains may require fewer resources to maintain focused attention and executive control, leading to decreased variability in neural signals.

In contrast, although the control group also showed reduced complexity of brain regions in the WM and GO/NO-GO tasks, the degree and extent of change was not as significant as in the intervention group. This situation may be due to the fact that even though the children were not trained in mindfulness, they experienced the development of brain functioning over a two-month period, resulting in differential changes in the pre- and post-tests.

Positive effects of mindfulness training on cognitive functioning

To further explore the positive effects of mindfulness training on children’s cognitive functioning, brain complexity differences between the intervention and control groups before and after mindfulness training was examined. It was found that mindfulness training showed a positive effect on the performance of the DCCS and WM tasks, particularly significantly in the left DLPFC. Previous research has shown that the left DLPFC develops more rapidly during early childhood, when language skills develop rapidly. This development is critical for tasks involving verbal WM and logical reasoning. On the other hand, the right DLPFC continues to develop throughout childhood, supporting the maturation of spatial WM, attentional control, and emotion regulation (Goto et al. 2015). Our study provides evidence that mindfulness-based training has a positive effect on the development of CF and WM.

According to the IAA model (Shapiro et al. 2006) and the view of attentional modulation as a core component of the mindfulness mechanism (Tang et al. 2015), mindfulness training works by enhancing sustained attention, better monitoring, and effective switching between task sets (i.e. CF; Vieth and von Stockhausen 2022). The results of the present study are consistent with previous research showing that mindfulness training stabilizes attention and improves CF (Lee and Orsillo 2014; Zou et al. 2020). Individuals typically activate the DLPFC and parietal cortex during the early stages of mindfulness training, a phenomenon that supports the findings of studies conducted with preschoolers who are new to mindfulness training (Posner et al. 2015; Tang et al. 2015).

Furthermore, the current study showed that children in the mindfulness training group showed improved number of consecutive correct responses compared to the control group, suggesting that mindfulness interventions can promote sustained accuracy in WM tasks. This trend may imply that mindfulness training helps to reduce the stress-induced brain load in WM tasks, which in turn reduces the intensity of brain activity. This finding corroborates with a study by Jha et al. (2010), which found that mindfulness training was effective in reducing stress levels and facilitating better performance on a WM task. Similarly, Zeidan et al. (2010) explored the effects of short-term mindfulness training relative to an active control group (i.e. listening to an audiobook) on performance on an N-back task. It is worth noting that Mrazek et al. (2013) also found that mindfulness training was able to reduce distraction during the performance of a WM task. These about findings suggest that mindfulness training may reduce brain activity during WM tasks by reducing the burden of distracting thoughts.

However, similar behavioral and neural results were not observed in the Go/No-Go task, suggesting that mindfulness training may not have significantly improved children’s IC, although a large body of research has shown that the practice of mindfulness meditation can significantly enhance individuals’ IC on both behavioral and neural levels (Allen et al. 2012; Heeren et al. 2009; Moore and Malinowski 2009; Sahdra et al. 2011; Teper et al. 2013). This non-significant change may stem from design issues with the Go to No-Go ratio in the current study and a previous study (Xie et al. 2022), limiting children’s opportunities to perform on No-Go trials, which are typically associated with higher levels of mindfulness (Logemann-Molnár et al. 2022). Another reason might be that the spatiotemporal dynamics of IC change with age (Knežević and Marinković, 2017; Lin et al. 2018). Thus, IC is a long-term developmental function that requires longer interventions. Future research should conduct further studies using time-sensitive methods across the lifespan, exploring in more detail the methods of task design and behavioral outcome assessment to gain a more comprehensive understanding of the combined effects of cognitive training on the brain and behavior. Future research should examine in greater depth the role of these changes in actual behavioral performance, as well as their potential benefits for improving specific cognitive dysfunctions such as attention deficits.

Finally, the use of fNIRS in this study provided valuable insights into the task-related changes in cerebral hemodynamics associated with mindfulness-based training in young children. However, it is important to acknowledge the limitations inherent in interpreting fNIRS data as direct indicators of underlying neurophysiological or functional changes in cortical organization. While the findings suggest alterations in brain activation patterns following the intervention, caution must be exercised when extrapolating these hemodynamic responses to infer specific changes in cortical function or organization. fNIRS primarily measures changes in cerebral blood flow and oxygenation, which are influenced by various factors, including vascular dynamics, systemic physiological changes, and motion artifacts (Scholkmann et al. 2014). As such, observed hemodynamic responses may not always reflect direct neural activity changes and should be interpreted with caution. Furthermore, the spatial resolution of fNIRS is limited compared to other neuroimaging techniques such as fMRI (Boas et al. 2014). This limitation may restrict our ability to precisely localize brain activation changes and infer their neuroanatomical correlates. Therefore, future studies could benefit from integrating fNIRS data with complementary neuroimaging modalities to provide a more comprehensive understanding of the neural mechanisms underlying mindfulness-based interventions (Eggebrecht et al. 2014).

Conclusion

Although our study demonstrated that the intervention had an effect on executive functioning, there were many factors involved (e.g. age, gender, IQ, etc.). In future studies, we should categorize children by age and gender and investigate whether mindfulness training has different improvement effects at different ages and genders. Further, we can explore whether mindfulness training is more effective in children with high IQ. In studying the development of executive functioning in children, we need to realize that children respond in a variety of ways; therefore, tasks to assess executive functioning should reflect this diversity. In particular, in the go/No-Go task, the present study failed to observe a significant effect of mindfulness training, which may be related to problems with the task design and children’s tendency to become tired and have difficulty sustaining their concentration during the task. Future studies should appropriately increase the proportion of No-Go trials and consider toing children’s interest and concentration. Follow-up studies should also determine whether these interventions are effective in increasing inhibitory efficiency and improving cognitive domains mediated by IC. Stronger evidence for the improvement of brain function by mindfulness-based training in the other two tasks is needed. Further exploration is needed to determine what factors may influence the effectiveness of mindfulness interventions. Most of the current studies vary widely in the duration of the intervention; therefore, determining the optimal amount of mindfulness practice is important for studies that promote executive functioning (or other outcomes of interest).

Another limitation of this study is the reliance on fNIRS as the primary neuroimaging modality. While fNIRS offers several advantages, including portability and fewer participant constraints compared to fMRI, it is essential to acknowledge its limitations in accurately capturing the complex dynamics of cortical function. The study’s interpretation of observed hemodynamic changes as indicative of real neurophysiological or functional changes in cortical organization should be approached with caution. Future research should aim to validate fNIRS findings using complementary techniques and consider alternative methods for investigating cortical organization changes more directly.

The current study used a new method to reveal changes in brain regions by calculating the entropy of fNIRS data during the experiment. It was found that entropy reveals areas that produce neuronal correlations and is consistent with the results of traditional methods of analyzing neuronal correlations, as well as providing information not found in mean estimates. The current study found that 15-session of mindfulness training enhanced CF and WM in young children. These brain region changes in EF were evidenced by significant decreases in brain complexity during CF and WM tasks, suggesting that preschoolers can benefit from mindfulness training. Research on EF interventions in schools is in its infancy, and future studies should include generalization and optimal developmental period factors to maximize intervention effects.

Supplementary Material

Supplementary_bhae330

Acknowledgments

The authors thank Principal Tan Jingjie and the participating teachers and children.

Author contributions

Sha Xie (Conceptulization, Data curation, Funding acquisition, Writing—review & editing), Shuqi Lu (Formal analysis, Investigation, Methodology, Writing—original draft), Jiahao Lu (Data curation, Validation), Chaohui Gong (Project administration), Chunqi Chang (Supervision, Funding acquisition).

Funding

This research was funded by the National Natural Science Foundation of China (No. 61971289), the Guangdong Planning Office of Philosophy and Social Science, China (No. GD21YJY09), the Shenzhen Education Science Planning Project, China (No. zdzz21002), the Shenzhen-Hong Kong Institute of Brain Science-Shenzhen Fundamental Research Institutions (No. 2023SHIBS0003), and the SZU Top Ranking Project (No. 86000000210).

 

Conflict of interest statement: The authors have no conflict of interests to declare.
==== Refs
References

Ahmed  MU, Mandic  DP. Multivariate multiscale entropy: a tool for complexity analysis of multichannel data. Phys Rev E. 2011:84 (6 ):061918. 10.1103/PhysRevE.84.061918.
Ahmed  MU, Mandic  DP. Multivariate multiscale entropy analysis. IEEE Signal Process Lett. 2012:19 (2 ):91–94. 10.1109/LSP.2011.2180713.
Allen  M, Dietz  M, Blair  KS, van Beek  M, Rees  G, Vestergaard-Poulsen  P, et al.  Cognitive-affective neural plasticity following active-controlled mindfulness intervention. J Neurosci. 2012:32 (44 ):15601–15610.23115195
Anderson  PJ . Towards a developmental model of executive function. In: Anderson V, Jacobs R, Anderson PJ, editors. Executive functions and the frontal lobes: a lifespan perspective. New York: Taylor & Francis; 2008, pp. 3–21.
Angsuwatanakul  T, Iramina  K, Kaewkamnerdpong  B. Brain complexity analysis of functional near infrared spectroscopy for working memory study. In: 2015 8th Biomedical Engineering International Conference (BMEiCON), Pattaya, Thailand: IEEE; 2015, November, pp. 1–5.
Baddeley  AD, Hitch  GJ. Developments in the concept of working memory. Neuropsychology. 1994:8 (4 ):485.
Blair  C . School readiness: integrating cognition and emotion in a neurobiological conceptualization of children's functioning at school entry. Am Psychol. 2002:57 (2 ):111.11899554
Blair  C . Developmental science and executive function. Curr Dir Psychol Sci. 2016:25 (1 ):3–7.26985139
Blakey  E, Visser  I, Carroll  DJ. Different executive functions support different kinds of cognitive flexibility: evidence from 2-, 3-, and 4-year-olds. Child Dev. 2016:87 (2 ):513–526.26659697
Boas  DA, Elwell  CE, Ferrari  M, Taga  G. Twenty years of functional near-infrared spectroscopy: introduction for the special issue. NeuroImage. 2014:85 :1–5.24321364
Bryant  L . Physical activity, structured sport participation, executive function in preschools. Purdue University, Doctoral dissertation; 2019.
Buckner  RL, DiNicola  LM. The brain’s default network: updated anatomy, physiology and evolving insights. Nat Rev Neurosci. 2019:20 (10 ):593–608.31492945
Busa  MA, Van Emmerik  REA. Multiscale entropy: a tool for understanding the complexity of postural control. J Sport Health Sci. 2016:5 (1 ):44–51. 10.1016/j.jshs.2016.01.018.30356502
Buss  AT, Fox  N, Boas  DA, Spencer  JP. Probing the early development of visual working memory capacity with functional near-infrared spectroscopy. NeuroImage. 2014:85 :314–325.23707803
Cao  L . Practical method for determining the minimum embedding dimension of a scalar time series. Phys D: Nonlinear Phenomena. 1997:110 (1–2 ):43–50. 10.1016/S0167-2789(97)00118-8.
Cao  M, Huang  H, He  Y. Developmental connectomics from infancy through early childhood. Trends Neurosci. 2017:40 (8 ):494–506.28684174
Carlson SM, Faja S, Beck DM.  Incorporating early development into the measurement of executive function: The need for a continuum of measures across development. In: Griffin JA, McCardle P, Freund LS, editors. Executive function in preschool-age children: Integrating measurement, neurodevelopment, and translational research (pp. 45–64). American Psychological Association: Washington, D.C.; 2016.
Carlson  SM, Mandell  DJ, Williams  L. Executive function and theory of mind: stability and prediction from ages 2 to 3. Dev Psychol. 2004:40 (6 ):1105.15535760
Costa  M, Goldberger  AL, Peng  CK. Multiscale entropy analysis of biological signals. Phys Rev E. 2005:71 (2 ):021906.
Cragg  L, Nation  K. Go or no-go? Developmental improvements in the efficiency of response inhibition in mid-childhood. Dev Sci. 2008:11 (6 ):819–827.19046150
Davidson  MC, Amso  D, Anderson  LC, Diamond  A. Development of cognitive control and executive functions from 4 to 13 years: evidence from manipulations of memory, inhibition, and task switching. Neuropsychologia. 2006:44 (11 ):2037–2078.16580701
Diamond  A . Executive functions. Annu Rev Psychol. 2013:64 (1 ):135–168. 10.1146/annurev-psych-113011-143750.23020641
Diamond  A, Lee  K. Interventions shown to aid executive function development in children 4 to 12 years old. Science. 2011:333 (6045 ):959–964.21852486
Ding  K, Wang  H, Li  C, Li  H. Decreased frontal lobe complexity in left-behind children during joint attention: a fNIRS study with multivariable and multiscale sample entropy analysis. Cereb Cortex. 2023:33 (21 ):10949–10958. 10.1093/cercor/bhad341.37727984
Eggebrecht  AT, Ferradal  SL, Robichaux-Viehoever  A, Hassanpour  MS, Dehghani  H, Snyder  AZ, et al.  Mapping distributed brain function and networks with diffuse optical tomography. Nat Photonics. 2014:8 (6 ):448–454.25083161
Elbert  T, Ray  WJ, Kowalik  ZJ, Skinner  JE, Graf  KE, Birbaumer  N. Chaos and physiology: deterministic chaos in excitable cell assemblies. Physiol Rev. 1994:74 (1 ):1–47.8295931
Engelmann  JM, Versace  F, Robinson  JD, Minnix  JA, Lam  CY, Cui  Y, Brown  VL, Cinciripini  PM. Neural substrates of smoking cue reactivity: a meta-analysis of fMRI studies. NeuroImage. 2012:60 (1 ):252–262. 10.1016/j.neuroimage.2011.12.024.22206965
Farb  NA, Segal  ZV, Mayberg  H, Bean  J, McKeon  D, Fatima  Z, Anderson  AK. Attending to the present: mindfulness meditation reveals distinct neural modes of self-reference. Soc Cogn Affect Neurosci. 2007:2 (4 ):313–322.18985137
Ferrari  M, Quaresima  V. A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application. NeuroImage. 2012:63 (2 ):921–935. 10.1016/j.neuroimage.2012.03.049.22510258
Fiore  F, Borella  E, Mammarella  IC, De Beni  R. Age differences in verbal and visuo-spatial working memory updating: evidence from analysis of serial position curves. Memory. 2012:20 (1 ):14–27.22133192
Flook  L, Goldberg  SB, Pinger  L, Davidson  RJ. Promoting prosocial behavior and self-regulatory skills in preschool children through a mindfulness-based Kindness Curriculum. Dev Psychol. 2015:51 (1 ):44.25383689
Fournet  N, Roulin  JL, Vallet  F, Beaudoin  M, Agrigoroaei  S, Paignon  A, et al.  Evaluating short-term and working memory in older adults: French normative data. Aging Ment Health. 2012:16 (7 ):922–930.22533476
Garon  N, Bryson  SE, Smith  IM. Executive function in preschoolers: a review using an integrative framework. Psychol Bull. 2008:134 (1 ):31.18193994
Ghouse  A, Nardelli  M, Valenza  G. fNIRS complexity analysis for the assessment of motor imagery and mental arithmetic tasks. Entropy. 2020:22 (7 ):761.33286533
Gilmore  JH, Knickmeyer  RC, Gao  W. Imaging structural and functional brain development in early childhood. Nat Rev Neurosci. 2018:19 (3 ):123–137.29449712
Gothe  NP, Kramer  AF, McAuley  E. The effects of an 8-week Hatha yoga intervention on executive function in older adults. J Gerontol A Biol Sci Med Sci. 2014:69 (9 ):1109–1116.25024234
Goto  T, Kita  Y, Suzuki  K, Koike  T, Inagaki  M. Lateralized frontal activity for Japanese phonological processing during child development. Front Hum Neurosci. 2015:9 :417.26236223
Griffin  JA, Freund  LS, McCARDLE  P, DelCarmen-Wiggins  R, Haydon  A. Introduction to executive function in preschool-age children. Washington, D.C.: American Psychological Association; 2016.
Guedes  C, Ferreira  T, Leal  T, Cadima  J. Unique and joint contributions of behavioral and emotional self-regulation to school readiness. Appl Dev Sci. 2023:27 (2 ):136–155.
Gupta  R, Kar  BR, Srinivasan  N. Development of task switching and post-error-slowing in children. Behav Brain Funct. 2009:5 (1 ):1–13.19121215
Hasher  L, Zacks  RT. Working memory, comprehension, and aging: a review and a new view. Psychol Learn Motiv. 1988:22 :193–225.
Hasher  L, Stoltzfus  ER, Zacks  RT, Rypma  B. Age and inhibition. J Exp Psychol Learn Mem Cogn. 1991:17 (1 ):163.1826730
Heeren  A, Van Broeck  N, Philippot  P. The effects of mindfulness on executive processes and autobiographical memory specificity. Behav Res Ther. 2009:47 (5 ):403–409.19232573
Hendry  A, Jones  EJ, Charman  T. Executive function in the first three years of life: precursors, predictors and patterns. Dev Rev. 2016:42 :1–33.
Hölzel  BK, Ott  U, Hempel  H, Hackl  A, Wolf  K, Stark  R, Vaitl  D. Differential engagement of anterior cingulate and adjacent medial frontal cortex in adept meditators and non-meditators. Neurosci Lett. 2007:421 (1 ):16–21.17548160
Hölzel  BK, Lazar  SW, Gard  T, Schuman-Olivier  Z, Vago  DR, Ott  U. How does mindfulness meditation work? Proposing mechanisms of action from a conceptual and neural perspective. Perspect Psychol Sci. 2011:6 (6 ):537–559.26168376
Huppert  TJ, Diamond  SG, Franceschini  MA, Boas  DA. HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain. Appl Opt. 2009:48 (10 ):D280. 10.1364/AO.48.00D280.19340120
Huttenlocher  PR, Dabholkar  AS. Developmental anatomy of prefrontal cortex. In: Krasnegor NA, Reid Lyon G, Goldman-Rakic PS, editors. Development of the prefrontal cortex: evolution, neurobiology, and behavior. Baltimore: Paul H Brookes Pub Co., 1997, pp. 69–83.
Inoue  Y, Sakihara  K, Gunji  A, Ozawa  H, Kimiya  S, Shinoda  H, et al.  Reduced prefrontal hemodynamic response in children with ADHD during the Go/NoGo task: a NIRS study. Neuroreport. 2012:23 (2 ):55–60.22146580
Jack  AI, Dawson  AJ, Begany  KL, Leckie  RL, Barry  KP, Ciccia  AH, Snyder  AZ. fMRI reveals reciprocal inhibition between social and physical cognitive domains. NeuroImage. 2013:66 :385–401.23110882
Jha  AP, Stanley  EA, Kiyonaga  A, Wong  L, Gelfand  L. Examining the protective effects of mindfulness training on working memory capacity and affective experience. Emotion. 2010:10 (1 ):54.20141302
Kabat-Zinn  J, Hanh  TN. Full catastrophe living: using the wisdom of your body and mind to face stress, pain, and illness. New York, NY: Dell Pub., a division of Bantam Doubleday Dell Pub. Group, 1991.
Klingberg  T . Training and plasticity of working memory. Trends Cogn Sci. 2010:14 (7 ):317–324.20630350
Knežević  M, Marinković  K. Neurodynamic correlates of response inhibition from emerging to mid adulthood. Cogn Dev. 2017:43 :106–118.29081593
Lahat  A, Todd R, Mahy CE, Lau K, Zelazo PD. Neurophysiological correlates of executive function: a comparison of European-Canadian and Chinese-Canadian 5-year-olds. Front Hum Neurosci. 2010:3 :956.
Lee  JK, Orsillo  SM. Investigating cognitive flexibility as a potential mechanism of mindfulness in generalized anxiety disorder. J Behav Ther Exp Psychiatry. 2014:45 (1 ):208–216.24239587
Li  X, Zhu  Z, Zhao  W, Sun  Y, Wen  D, Xie  Y, et al.  Decreased resting-state brain signal complexity in patients with mild cognitive impairment and Alzheimer’s disease: a multi-scale entropy analysis. Biomed Opt Express. 2018:9 (4 ):1916–1929.29675329
Li  Q, Song  Y, Lian  B, Feng  T. Mindfulness training can improve 3-and 4-year-old children’s attention and executive function. Acta Psychol Sin. 2019:51 (3 ):324.
Lin  MY, Tseng  YJ, Cheng  CH. Age effects on spatiotemporal dynamics of response inhibition: an MEG study. Front Aging Neurosci. 2018:10 :386.30515093
Logemann-Molnár  Z, Akıl  AM, Cserjési  R, Nagy  T, Veres-Székely  A, Demetrovics  Z, Logemann  A. Dispositional mindfulness and inhibitory control after perceiving neutral, food and money related stimuli. Int J Environ Res Public Health. 2022:19 (4 ):2201.35206388
Lorenzetti  V, Gaillard  A, Beyer  E, Kowalczyk  M, Kamboj  SK, Manning  V, Gleeson  J. Do mindfulness-based interventions change brain function in people with substance dependence? A systematic review of the fMRI evidence. BMC Psychiatry. 2023:23 (1 ):407.37286936
Manjunath  NK, Telles  S. Improved performance in the Tower of London test following yoga. Indian J Physiol Pharmacol. 2001:45 (3 ):351–354.11881575
McLaughlin  KA, Sheridan  MA, Lambert  HK. Childhood adversity and neural development: deprivation and threat as distinct dimensions of early experience. Neurosci Biobehav Rev. 2014:47 :578–591.25454359
McLaughlin  KA, Weissman  D, Bitrán  D. Childhood adversity and neural development: a systematic review. Annu Rev Dev Psychol. 2019:1 :277–312.32455344
Mehnert  J, Akhrif  A, Telkemeyer  S, Rossi  S, Schmitz  CH, Steinbrink  J, et al.  Developmental changes in brain activation and functional connectivity during response inhibition in the early childhood brain. Brain and Development. 2013:35 (10 ):894–904.23265620
Menon  V, D’Esposito  M. The role of PFC networks in cognitive control and executive function. Neuropsychopharmacology. 2022:47 (1 ):90–103.34408276
Moffitt  TE, Arseneault  L, Belsky  D, Dickson  N, Hancox  RJ, Harrington  H, et al.  A gradient of childhood self-control predicts health, wealth, and public safety. Proc Natl Acad Sci. 2011:108 (7 ):2693–2698.21262822
Molavi  B, Dumont  GA. Wavelet-based motion artifact removal for functional near-infrared spectroscopy. Physiol Meas. 2012:33 (2 ):259–270. 10.1088/0967-3334/33/2/259.22273765
Moore  A, Malinowski  P. Meditation, mindfulness and cognitive flexibility. Conscious Cogn. 2009:18 (1 ):176–186.19181542
Moriguchi  Y . Neural mechanisms of executive function development during early childhood. In: Executive function. New York: Routledge; 2017, pp. 75–90.
Moriguchi  Y, Hiraki  K. Longitudinal development of prefrontal function during early childhood. Dev Cogn Neurosci. 2011:1 (2 ):153–162.22436437
Moriguchi  Y, Itakura  S. Young children’s difficulty with inhibitory control in a social context 1. Jpn Psychol Res. 2008:50 (2 ):87–92. 10.1111/j.1468-5884.2008.00364.x.
Moriguchi  Y, Lee  K, Itakura  S. Social transmission of disinhibition in young children. Dev Sci. 2007:10 (4 ):481–491. 10.1111/j.1467-7687.2007.00601.x.17552937
Mrazek  MD, Franklin  MS, Phillips  DT, Baird  B, Schooler  JW. Mindfulness training improves working memory capacity and GRE performance while reducing mind wandering. Psychol Sci. 2013:24 (5 ):776–781.23538911
Munakata  Y, Snyder  HR, Chatham  CH. Developing cognitive control: three key transitions. Curr Dir Psychol Sci. 2012:21 (2 ):71–77.22711982
Perlman  SB, Huppert  TJ, Luna  B. Functional near-infrared spectroscopy evidence for development of prefrontal engagement in working memory in early through middle childhood. Cereb Cortex. 2016:26 (6 ):2790–2799.26115660
Perpetuini  D, Bucco  R, Zito  M, Merla  A. Study of memory deficit in Alzheimer’s disease by means of complexity analysis of fNIRS signal. Neurophotonics. 2018:5 (1 ):011010–011010.28983489
Perpetuini  D, Chiarelli  AM, Cardone  D, Filippini  C, Bucco  R, Zito  M, Merla  A. Complexity of frontal cortex fNIRS can support Alzheimer disease diagnosis in memory and visuo-spatial tests. Entropy. 2019:21 (1 ):26.33266742
Perpetuini  D, Chiarelli  AM, Filippini  C, Cardone  D, Croce  P, Rotunno  L, et al.  Working memory decline in Alzheimer’s disease is detected by complexity analysis of multimodal EEG-fNIRS. Entropy. 2020:22 (12 ):1380.33279924
Perry  LK, Samuelson  LK, Burdinie  JB. Highchair philosophers: the impact of seating context-dependent exploration on children's naming biases. Dev Sci. 2014:17 (5 ):757–765.24289734
Posner  MI, Rothbart  MK, Tang  YY. Enhancing attention through training. Curr Opin Behav Sci. 2015:4 :1–5.
Razza  RA, Bergen-Cico  D, Raymond  K. Enhancing preschoolers’ self-regulation via mindful yoga. J Child Fam Stud. 2015:24 :372–385.
Sahdra  BK, MacLean  KA, Ferrer  E, Shaver  PR, Rosenberg  EL, Jacobs  TL, et al.  Enhanced response inhibition during intensive meditation training predicts improvements in self-reported adaptive socioemotional functioning. Emotion. 2011:11 (2 ):299.21500899
Scheibel  RS, Levin  HS. Frontal lobe dysfunction following closed head injury in children: Findings from neuropsychology and brain imaging. In Krasnegor NA, Lyon GR, Goldman-Rakic PS, editors. Development of the prefrontal cortex: Evolution, neurobiology, and behavior (pp. 241–263). Paul H. Brookes Publishing Co.: Baltimore, 1997.
Scholkmann  F, Kleiser  S, Metz  AJ, Zimmermann  R, Pavia  JM, Wolf  U, Wolf  M. A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology. NeuroImage. 2014:85 :6–27.23684868
Shapiro  SL, Carlson  LE, Astin  JA, Freedman  B. Mechanisms of mindfulness. J Clin Psychol. 2006:62 (3 ):373–386.16385481
Smith  EE, Jonides  J. Storage and executive processes in the frontal lobes. Science. 1999:283 (5408 ):1657–1661.10073923
Tang  Y . Health from brain, wisdom from brain. Dalian: Dalian University of Technology Press; 2005  (including training and practice VCDs) (in Chinese).
Tang  YY . Multi-intelligence and unfolding the full potentials of brain. Dalian: Dalian University of Technology Press; 2007  (in Chinese).
Tang  YY, Posner  MI. Attention training and attention state training. Trends Cogn Sci. 2009:13 (5 ):222–227.19375975
Tang  YY, Yang  L, Leve  LD, Harold  GT. Improving executive function and its neurobiological mechanisms through a mindfulness-based intervention: advances within the field of developmental neuroscience. Child Dev Perspect. 2012:6 (4 ):361–366.25419230
Tang  YY, Hölzel  BK, Posner  MI. The neuroscience of mindfulness meditation. Nat Rev Neurosci. 2015:16 (4 ):213–225.25783612
Teper  R, Segal  ZV, Inzlicht  M. Inside the mindful mind: how mindfulness enhances emotion regulation through improvements in executive control. Curr Dir Psychol Sci. 2013:22 (6 ):449–454.
Thera  N . The power of mindfulness. Sri Lanka: Buddhist Publication Society; 1972.
Tsujimoto  S, Yamamoto  T, Kawaguchi  H, Koizumi  H, Sawaguchi  T. Prefrontal cortical activation associated with working memory in adults and preschool children: an event-related optical topography study. Cereb Cortex. 2004:14 (7 ):703–712.15084489
Vieth  E, von Stockhausen  L. Mechanisms underlying cognitive effects of inducing a mindful state. J Cogn. 2022:5 (1 ):13.
Wang  S, Tang  H, Wang  B, Mo  J. Analysis of fatigue in the biceps brachii by using rapid refined composite multiscale sample entropy. Biomed Signal Process Control. 2021:67 :102510.
Wiebe  SA, Sheffield  TD, Espy  KA. Separating the fish from the sharks: a longitudinal study of preschool response inhibition. Child Dev. 2012:83 (4 ):1245–1261. 10.1111/j.1467-8624.2012.01765.x.22533453
Xie  S, Gong  C, Lu  J, Li  H, Wu  D, Chi  X, Chang  C. Enhancing Chinese preschoolers’ executive function via mindfulness training: an fNIRS study. Front Behav Neurosci. 2022:16 :961797.36090651
Yee  E, Thompson-Schill  SL. Putting concepts into context. Psychon Bull Rev. 2016:23 :1015–1027.27282993
Yue  WL, Ng  KK, Koh  AJ, Perini  F, Doshi  K, Zhou  JH, Lim  J. Mindfulness-based therapy improves brain functional network reconfiguration efficiency. Transl Psychiatry. 2023:13 (1 ):345.37951943
Zeidan  F, Johnson  SK, Diamond  BJ, David  Z, Goolkasian  P. Mindfulness meditation improves cognition: evidence of brief mental training. Conscious Cogn. 2010:19 (2 ):597–605.20363650
Zelazo  PD . The Dimensional Change Card Sort (DCCS): a method of assessing executive function in children. Nat Protoc. 2006:1 (1 ):297–301. 10.1038/nprot.2006.46.17406248
Zelazo  PD, Müller  U. Executive function in typical and atypical development. In: Goswami  U, editors. Blackwell handbook of childhood cognitive development. Chichester: Blackwell Publishing; 2002. pp. 445–469.
Zelazo  PD, Forston  JL, Masten  AS, Carlson  SM. Mindfulness plus reflection training: effects on executive function in early childhood. Front Psychol. 2018:9 :208.29535661
Zhang  T, Huang  W, Wu  X, Sun  W, Lin  F, Sun  H, Li  J. Altered complexity in resting-state fNIRS signal in autism: a multiscale entropy approach. Physiol Meas. 2021:42 (8 ):085004. 10.1088/1361-6579/ac184d.
Zhu  L, Haghani  S, Najafizadeh  L. On fractality of functional near-infrared spectroscopy signals: analysis and applications. Neurophotonics. 2020:7 (2 ):025001–025001.32377544
Zou  Y, Li  P, Hofmann  SG, Liu  X. The mediating role of non-reactivity to mindfulness training and cognitive flexibility: a randomized controlled trial. Front Psychol. 2020:11 :1053.32670135
