Abstract
Very preterm (VPT) young adolescents are at high risk of executive, behavioural and socio-emotional difficulties. Previous research has shown significant evidence of the benefits of mindfulness-based intervention (MBI) on these abilities. This study aims to assess the association between the effects of MBI on neurobehavioral functioning and changes in white-matter microstructure in VPT young adolescents who completed an 8-week MBI program. Neurobehavioural assessments and multi-shell diffusion MRI were performed before and after MBI in 32 VPT young adolescents. Combined diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) measures were extracted on well-defined white matter tracts (TractSeg). A multivariate data-driven approach (partial least squares correlation) was used to explore associations between MBI-related changes on neurobehavioural measures and microstructural changes. Our finding showed an enhancement of global executive functioning after MBI that was associated with a general pattern of significant increase in fractional anisotropy (FA) and decrease in axonal dispersion (ODI) in white-matter tracts involved in executive processes. Young VPT adolescents with lower gestational age at birth showed the greatest gain in white-matter microstructural changes after MBI.
Highlights
Very preterm adolescents (VPT) completed an 8-week mindfulness-based intervention (MBI)
Improvement in overall executive functioning was observed after MBI
Executive gain was associate with white-matter microstructural changes
The increase in microstructural properties was in tracts involved in executive processes
VPT with lower gestational age show bigger gains in microstructural changes
CRediT roles Vanessa Siffredi: Conceptualization; Data curation; Formal analysis; Investigation; Methodology; Project administration; Software; Visualization; Writing - original draft; Writing - review & editing - Maria Chiara Liverani: Conceptualization; Data curation; Investigation; Methodology; Project administration; Writing - review & editing. - Dimitri Van De Ville: Methodology; Resources; Software; Supervision; Writing - review & editing. - Lorena Freitas: Data curation; Investigation; Writing - review & editing. - Cristina Borradori Tolsa: Conceptualization; Funding acquisition; Investigation; Project administration; Resources; Supervision; Validation; Writing - review & editing. - Petra Susan Hüppi: Conceptualization; Funding acquisition; Methodology; Project administration; Resources; Supervision; Validation; Writing - review & editing. -Russia Hà-Vinh Leuchter: Conceptualization; Funding acquisition; Investigation; Methodology; Project administration; Resources; Supervision; Validation; Writing - review & editing.
Introduction
Children and adolescents born very preterm (VPT; <32 completed weeks of gestation) are at increased risk for executive, behavioural and socio-emotional impairments that persist into adolescence and adulthood. Executive functioning (EF) is essential for goal-directed and adaptive problem-solving and behaviour. It has been conceptualised in four distinct subdomains: (i) attentional control, (ii) information processing, (iii) cognitive flexibility, and (iv) goal setting (Anderson, 2002). Behavioural and socio-emotional competences refer to a set of skills related to how individuals identify, express, understand, use and regulate their behaviours as well as their emotions and those of others (Mikolajczak, Quoidbach, Kotsou, & Nelis, 2009). Importantly, these competences are crucial in daily life activities and are closely linked to academic abilities and social behaviours (Best, Miller, & Naglieri, 2011; Carter, Briggs-Gowan, & Davis, 2004; Izard, 2011; Vaughan & Giovanello, 2010). Numerous studies show that white matter alterations are highly common after preterm birth and persist until adolescence and adulthood (Volpe, 2003). Importantly, these white-matter changes have been associated with executive, behavioural and socioemotional deficits in this population (Allin et al., 2011; K. M. Jones, Champion, & Woodward, 2013; Loe, Lee, & Feldman, 2013; Montagna & Nosarti, 2016; Nosarti, Allin, Frangou, Rifkin, & Murray, 2005; Perlman, 1998; Soria-Pastor et al., 2008; Thompson et al., 2014; Woodward, Clark, Pritchard, Anderson, & Inder, 2011).
Mindfulness-based intervention (MBI) - commonly defined as the ongoing monitoring of present-moment experience while attending to it in an open and accepting way and without judgment (Kabat-Zinn, 2003) - has been associated with enhanced executive, behavioural and socioemotional functioning in typically developing children and adolescents (Barnes, Bauza, & Treiber, 2003; Beauchemin, Hutchins, & Patterson, 2008; Biegel, Brown, Shapiro, & Schubert, 2009; Black, Milam, & Sussman, 2009; Broderick & Metz, 2009; Felver, Tipsord, Morris, Racer, & Dishion, 2017; Flook et al., 2010; Geronimi, Arellano, & Woodruff-Borden, 2019; Lee, Semple, Rosa, & Miller, 2008; Napoli, Krech, & Holley, 2004; Schonert-Reichl & Lawlor, 2010; Schonert-Reichl et al., 2015; Semple, Lee, Rosa, & Miller, 2009). In line with behavioural improvements, studies show that MBI may induce structural neuroplastic changes (Marchand, 2014). In this context, diffusion-weighted magnetic resonance imaging (DW-MRI) and the diffusion tensor imaging (DTI) model has been shown to be sensitive to subtle white matter microstructural changes. Fractional anisotropy (FA) has been shown to reflect functionally relevant microstructural properties of white matter, including axonal architecture, extent of myelination and density of axonal fibres comprising axonal bundles (C. Beaulieu, 2002; Boorman, O’Shea, Sebastian, Rushworth, & Johansen-Berg, 2007). Increased FA measures have been associated with MBI in adult populations in different brain regions and tracts, including frontal regions (Kang et al., 2013; Luders, Clark, Narr, & Toga, 2011), callosal regions (Laneri et al., 2016; Luders et al., 2012), anterior cingulate regions (Tang, Lu, Fan, Yang, & Posner, 2012), insular (Sharp et al., 2018), ucinate fasciculus (Luders et al., 2011) and superior longitudinal fasciculus (Hölzel et al., 2016; Luders et al., 2011). In adolescents, a recent study also showed increased FA in the superior longitudinal fasciculus after a 2-weeks mindfulness video game (Patsenko et al., 2019). To our knowledge, the DTI model is the only methods employed so far to explore potential structural neuroplasticity induced by MBI.
We recently reported a randomised controlled trial showing significant benefits of MBI on executive, behavioural and socio-emotional functioning in preterm young adolescents (Siffredi, Liverani, Hüppi, et al., 2021). The current study aims to assess the association between the benefits of MBI on neurobehavioral functioning and changes in microstructural brain properties. To this end, we first compared preterm and full-term young adolescents on a range of executive, behavioural and socio-emotional measures. Secondly, we assessed the benefit of an 8-weeks MBI on the neurobehavioural measures in which preterm young adolescents showed significant difficulties (i.e., significant reduction in neurobehavioural scores compared to the full-term group). Finally, we explored the association between changes on neurobehavioural measures after MBI and changes on white-matter microstructural measures. As described above, previous findings have revealed changes in quantitative DTI metrics after MBI, specifically an increase in FA. One potential confound of the DTI analyses relates to the inherent non-specificity of the measurement, where observed changes in FA may be due to changes to any combination of axonal dispersion, axonal densities or fibre crossing (Christian Beaulieu, 2002; D. K. Jones, Knösche, & Turner, 2013). In an attempt to provide more specific information about white-matter microstructural changes associated with MBI, a biophysical model of diffusion, NODDI (neurite orientation dispersion and density imaging), was used. The NODDI model allows to distinguish two key variables contributing to changes in FA, including neurite density and fibre orientation dispersion (Zhang, Schneider, Wheeler-Kingshott, & Alexander, 2012). Combining the high sensitivity of the DTI model with the high specificity of the NODDI model might allow for a better understanding of white-matter microstructural changes associated with changes on neurobehavioural measures after MBI.
Methods
Participants
This study uses data collected as part of the ‘Mindful preterm teens’ study (Siffredi, Liverani, Magnus-Smith, et al., 2021). One hundred and sixty-five VPT young adolescents, i.e., born before 32 gestational weeks, were invited to participate in the study. They were aged 10 to 14 years born between 01.01.2003 and 31.12.2008 at the Geneva University Hospital, Switzerland, and followed at the Division of Child Development and Growth at the Geneva University Hospital. VPT young adolescents were excluded if they had an intelligence quotient below 70, sensory or physical disabilities (cerebral palsy, blindness, hearing loss), or an insufficient understanding of French. Moreover, some families declined to participate due to lack of time, lack of interest, geographical constraints or unreachability. A total of 63 participants were enrolled in the ‘Mindful preterm teens’ study and 52 of them completed the MBI and neurobehavioural assessment before and after MBI. Of the 52 young adolescents, 39 completed MRI scans before and after intervention, and 32 were included in the current diffusion MRI analyses (diffusion sequences not completed, n=5; high level of motion artefacts, n=2), Figure 1.
Participant flow chart. This figure shows the number of young adolescents who were enrolled into the MBI study, and of these, the number who had usable scans at both pre- and post-MBI time points and were thus included in the current analysis (Siffredi, Liverani, Magnus-Smith, et al., 2021). MBI, Mindfulness-based intervention; MRI, magnetic resonance imaging.
Moreover, 22 term-born young adolescents aged between 10 and 14 years old were recruited through the community and completed a neurobehavioural assessment similar to the VPT group prior to the MBI. This study was approved by the Swiss Ethics Committees on research involving humans, ID: 2015-00175. Written informed consent was obtained from the principal caregiver and from the participant.
Mindfulness-based intervention
The proposed MBI was designed by the authors, adapting well-known MBI programs, including Mindfulness-Based Stress Reduction (MBSR; (Kabat-Zinn, 2003)) and Mindfulness-Based Cognitive Therapy (MBCT; (Segal, Williams, & Teasdale, 2001)), to adolescents’ needs and language. The length of sessions and practices were also adapted to this population. The program consisted of 8 weekly sessions in groups of up to 8 participants, lasting 1h30. Two instructors were present for each group throughout the intervention (for further details on the MBI, see (Siffredi, Liverani, Magnus-Smith, et al., 2021)).
Neurobehavioral measures
Participants’ executive, behavioural and socio-emotional functioning were assessed using parent-report and self-report questionnaires, neuropsychological testing and computerised neurocognitive tasks, see Supplementary Table S1.
(i) Executive competences measures
Behaviours related to executive competences of young adolescents were assessed using the Behaviour Rating Inventory of Executive Function-parent version (BRIEF) (Gioia, Isquith, Guy, & Kenworthy, 2000). The BRIEF comprises 86 items over two standardised subscales, the Behavioural Regulation Index (BRI), which evaluates inhibition, flexibility and emotional control, and the Metacognition Index (MI), which evaluates initiation, working memory, planning/organization, self-regulation as well as a global score called the Global Executive Composite (GEC). Neurocognitive computerised tasks comprised: (i) the Flanker Visual Filtering Task, in which reaction time of the congruent condition was used to assess speed of processing (which belongs to the information processing subdomain), and the inhibition score (reaction time in incongruent conditions – reaction time in congruent conditions) was used as a measure of the subdomain of attentional control (Anderson, 2002; Christ, Kester, Bodner, & Miles, 2011); (ii) the child-adapted version of the Reality Filtering task, in which the temporal context confusion index (TCC) was used as a reality filtering measure, which involves integration of different executive processes (Liverani et al., 2020; Liverani et al., 2017). Neuropsychological testing included the Letter-Number Sequencing subtest from the Wechsler Intelligence Scale for Children, 4th Edition (WISC-IV) assessing working memory, which belongs to the cognitive flexibility subdomain (Anderson, 2002). Given the strong association between executive functions and mathematical abilities in children and adolescents,(Holmes & Adams, 2006; Visu-Petra, Cheie, Benga, & Miclea, 2011) we also used the total score of the Tempo Test Rekenen to assess timed mathematical achievement (De Vos, 1994).
(ii) Behavioural and socio-emotional competences measures
The total score of the Strength and Difficulties Questionnaire – parent version (SDQ) was used to assess behaviour in daily life (Goodman, 1997, 2001). Participants completed three self-reported questionnaires: the KIDSCREEN-27 items questionnaire was used to assess the quality of life of the participants (Robitail et al., 2007); the total score of the Social Goal Scale was used to assess social responsiveness and social relationships (Wentzel, 1994); and the total score of the Self-Compassion Scale – Short form was used to assess the main components of self-compassion (Raes, Pommier, Neff, & Van Gucht, 2011). Neuropsychological testing included the Affect Recognition subtest (A Developmental Neuropsychological Assessment, 2nd Edition - NEPSY-II), giving a total score assessing facial emotional recognition and the Theory of Mind subtest (NEPSY-II), giving a total score measuring the ability to understand mental functions, such as belief, intention or deception (Korkman, Kirk, & Kemp, 2007).
Magnetic resonance imaging
MRI data were acquired at Campus Biotech in Geneva, Switzerland, using a Siemens 3T Magnetom Prisma scanner. All participants completed a simulated “mock” MRI session prior to their first MRI scan. This preparation process was conducted by trained research staff and allowed participants to familiarise themselves with the scanner and the scanning process, eventually raising any concerns they might have had prior to the MRI scan. Furthermore, this process is known to facilitate the acquisition of good quality MRI images in children and adolescents (de Bie et al., 2010; Tamnes, Roalf, Goddings, & Lebel, 2018).
A multi-shell diffusion-weighted (DW) echo planar imaging (EPI) protocol was used and included four shells. The first sequence, referred to as ‘b200’, included 10 gradient directions with b-values of 200 s/mm2; the second one referred to as ‘b1700’, included 30 gradient directions with b-values of 1700 s/mm2; the third one referred to as ‘b4200a’, included 26 gradient directions with b-values of 4200 s/mm2; and the fourth one referred to as ‘b4200b’, included 24 gradient directions with b-values of 4200 s/mm2. Each of four sequences included the acquisition of 4 images with b-value = 0 s/mm2 images, and the parameter set for all sequences were: TR = 7000 ms, TE = 87 ms, FOV = 234 × 243 mm, slice thickness = 1.3 mm, voxel size = 1.3 × 1.3 × 1.3 mm.
Diffusion image preprocessing
Visual inspection of raw data for brain coverage, spike artefacts, severe head motion, and other severe image artefacts was completed and participants were excluded if necessary. The four diffusion shells (b200, b1700, b4200a, b4200b) were preprocessed independently using MRtrix3 (Tournier et al., 2019) and using the following pipeline: a) denoising (Cordero-Grande, Christiaens, Hutter, Price, & Hajnal, 2019; Veraart, Fieremans, & Novikov, 2016; Veraart, Novikov, et al., 2016);, b) Gibbs ringing removal (Kellner, Dhital, Kiselev, & Reisert, 2016), c) correction for movement and eddy current-induced geometric distortions using the eddy tool implemented in FSL (Jenkinson, Beckmann, Behrens, Woolrich, & Smith, 2012). The first b = 0 s/mm2 images of the b1700, b4200a, b4200b sequences were linearly registered to the first b = 0 s/mm2 image of the b200 sequence using FreeSurfer to bring them into b200 space before merging them together. The brain extraction tool (BET) from FSL (Smith, 2002) was then applied to the combined b200, b1700, b4200a, b4200b image to remove non-brain tissue and subsequently intensity normalisation was applied. Following Pines and colleague’s (2020) recommendations (Pines et al., 2020), the resulting multi-shell diffusion weighted image was then used for both DTI and NODDI models fitting and tractography.
Diffusion models fitting
The DTI model was applied to the resulting multi-shell diffusion weighted image and whole-brain maps of FA was calculated for each participant. FA (between 0 and 1) is a measure of the directionality of diffusion that characterise the variance of the three eigenvalues pairs that represent the direction and magnitude of diffusivity along the three orthogonal axes (v1, λ1; v2, λ2; v3, λ3) (Tamnes et al., 2018). In addition, intra-cellular volume fraction and orientation dispersion indices were estimated from the resulting multi-shell diffusion weighted image using the NODDI model (Zhang et al., 2012). The NODDI Matlab Toolbox http://www.nitrc.org/projects/noddi_toolbox was used to extract maps of intra-cellular volume fraction (ICVF) and fibre orientation dispersion (ODI) across the brain for each participant.
Tractography and tractometry measures
Whole-brain fibre orientation distributions (FOD) were estimated using the multi-shell multi-tissue constrained spherical deconvolution (MSMT-CSD) method (Jeurissen, Tournier, Dhollander, Connelly, & Sijbers, 2014), resulting in a condensed representation of diffusion along three principal fibre directions per voxel according to tissue type (grey, white, cortico-spinal fluid). Tractography-based Segmentation (TractSeg) use a supervised-learning approach with a aconvolutional neural network-based that directly segments tracts in the field of fibre orientation distribution function (fODF) peaks without using parcellation (Wasserthal, Neher, & Maier-Hein, 2018). TractSeg has achieved state-of-the-art performance and allows for an accurate reconstruction of fibre tracts in subject space, thus avoiding the problem of inaccurate coregistration of tracts with varying size and shape. Whole-brain fibre orientation distribution function (fODF) peaks map were input into a two stage fully convolutional neural network trained using segmented priors of 72 anatomically well-defined white matter tracts from the Human Connectome Project. Using the tractometry function, along-tract mean FA, ICVF and ODI were calculated for the 50 most consistent white-matter tracts in each subject’s native space (Chandio, Harezlak, & Garyfallidis, 2019), see Figure 2.
Flowchart summarizing the computation and analysis of microstructural properties of white-matter tracts. Briefly, diffusion-weighted images were preprocessed. Whole-brain multi-shell multi-tissue constrained spherical deconvolution (MSMT-CSD) was extracted and fibre orientation distribution (FOD) function peaks map were input into Tractography-based Segmentation (TractSeg) to extract the main white-matter tracts in subject space. Mean fractional anisotropy (FA), intra-cellular volume fraction (ICVF) and orientation dispersion index (ODI) were then extracted along each track by applying the diffusion tensor and NODDI models.
Statistical analyses
Statistical analyses were performed using R version 4.0.3 (R. C. Team, 2019) and RStudio version 1.3.1093 (R. Team, 2020) and employing the following methods:
Neurobehavioural functioning: comparison of the preterm and full-term groups
Group comparisons between the preterm and the full-term groups were conducted for all neurobehavioural measures using Wilcoxon signed rank test, given that assumptions for parametric testing were violated. P values were corrected for multiple comparisons with the false discovery rate (FDR, q values <0.05 (Benjamini & Hochberg, 1995)). Effect size were assessed using Wilcoxon effect size (r).
Impact of MBI on neurobehavioral functioning affected by preterm-birth
Paired-sample t-test were used to evaluate the impact of MBI on the neurobehavioural measures showing a significant deficit in the preterm group as compared to the full-term group. P values were corrected for multiple comparisons with the false discovery rate (FDR, q values <0.05 (Benjamini & Hochberg, 1995)). Effect size were assessed using Wilcoxon effect size (r).
Association between neurobehavioural changes and microstructural changes after MBI
For each neurobehavioural and mean microstructural measures (i.e., FA, ICVF and ODI), score differences between assessment at pre-MBI and assessment at post-MBI were calculated for each participant. Score differences for neurobehavioural measures will be referred to as Δ neurobehavioural measures (neurobehavioural score post-MBI − neurobehavioural score pre-MBI). Score differences for mean FA, ICVF and ODI measures will be referred to as Δ mean FA, ICVF and ODI measures (i.e., mean FA post-MBI − mean FA pre-MBI; mean ICVF post-MBI – mean ICVF pre-MBI; mean ODI post-MBI − mean ODI pre-MBI). Negative Δ indicates a reduction of the scores between two time points, whereas positive Δ indicates an increase in scores between two time points.
Partial least square correlation analyses (PLSC) were performed to evaluate associations between changes in neurobehavioural measures and changes in mean microstructural measures after MBI. PLSC is a data-driven multivariate technique that maximizes the covariance between two matrices by identifying latent components (LCs) which are linear combinations of the two matrices, i.e., Δ neurobehavioural functioning measures and Δ mean microstructural measures (McIntosh & Lobaugh, 2004). A publicly available Matlab PLSC implementation was used: https://github.com/danizoeller/myPLS (Kebets et al., 2019; Zoller et al., 2019).
In a first PLSC, associations between changes in neurobehavioural measures and changes in mean FA measures after MBI were investigated to capture any changes related to white-matter microstructural properties measured by the classical tensor model.
The neurobehavioural functioning data refers to the three Δ neurobehavioural measures showing significant difference between the full-term and preterm group (i.e., BRIEF MI, BRIEF GEC, SDQ total), as well as age at assessment and gestational age at birth. The neurobehavioural functioning data were stored in a 32 × 5 matrix denoted X. Each row of X represents one subject and the matrix’s 5 columns are made up of the three Δ neurobehavioural measures showing significant difference between the full-term and preterm group, as well as age at assessment and gestational age at birth. The Δ mean FA data were gathered in a 32 × 50 matrix denoted Y, with each row matching one subject and each column one Δ mean FA measure for each of the 50 tracts extracted using the TractSeg framework. A cross-covariance matrix (which is effectively a correlation matrix, since the data are z-scores) was then computed between X and Y. Singular value decomposition was then applied to this cross-covariance matrix, resulting in latent components. Each latent component is composed of neurobehavioural functioning saliences and Δ mean FA saliences, and saliences indicate how strongly each neurobehavioural functioning measures and Δ mean FA measures contribute to the multivariate association of neurobehavioral functioning and Δ mean FA. The significance of latent components was determined by permutation testing (1000 permutations). Stability of neurobehavioural functioning saliences and Δ mean FA saliences were estimated using bootstrapping (500 bootstrap samples with replacement). Bootstrap ratio z-scores for each neurobehavioural functioning and Δ mean FA measures were obtained by dividing each neurobehavioural functioning and Δ mean FA salience by its bootstrap-estimated standard deviation, and a p-value was obtained for each bootstrap ratio z-score. Following the PLSC interpretation (Krishnan, Williams, McIntosh, & Abdi, 2011), the contribution of neurobehavioural functioning and Δ mean FA saliences for a given latent component was considered robust at p < 0.01 (i.e., absolute bootstrap ratio z-scores above 2.3 or below −2.3).
In a second PLSC, we investigated further associations between changes in neurobehavioural measures and changes in microstructural measures using the NODDI model, including mean ICVF and ODI measures after MBI. A procedure similar to the first PLSC was employed using a 32 × 5 matrix denoted X containing the neurobehavioural functioning data and a 32 × 100 matrix denoted Y containing the Δ mean ICVF and Δ mean ODI.
Results
Participant characteristics
The final sample included 32 preterm-born and 22 full-term young adolescents between 10 and 14 years of age. Baseline characteristics were similar between preterm and full-term participants for sex, age at the assessment and socio-economic status. The preterm group had significantly lower general ability index compared to the full-term group, see Table 1.
Neonatal and demographic characteristics of the preterm and full-term participants
Neurobehavioural outcomes: Comparison of the preterm and full-term groups
Neurobehavioural measures showing a significant difference and surviving FDR correction between the preterm and the full-term groups are presented in Figure 3. The preterm group showed higher BRIEF GEC (q < 0.001), BRIEF MI (q < 0.001), BRIEF BRI (q < 0.001) and SDQ total scores (q < 0.001) reflecting more executive function and behavioural difficulties in daily life compared to the full-term group. Moreover, the self-compassion score was reduced in the preterm group compared to the full-term group (q < 0.011), see Figure 3 and Supplementary Table S2. The other neurobehavioural scores were comparable between the preterm and full-term groups.
Neurobehavioural scores showing significant group differences between the preterm (prior to the MBI) and the full-term groups using independent sample t-test (FDR, q values <0.05). The mean for each measure and each group is represented by a grey line.
Impact of MBI on neurobehavioral measures affected by preterm-birth
A significant decrease after MBI was observed for the following measures: BRIEF GEC (q < 0.00), BRIEF MI (q < 0.001) and SDQ total scores (q < 0.04), reflecting significantly less executive and behavioural difficulties in daily life. There was no significant difference before and after MBI for the BRIEF BRI and self-compassion scores, see Figure 4 and Supplementary Table S3.
Neurobehavioural scores before and after MBI. Paired-sample t-tests was used to compared scores before and after MBI (FDR, q values <0.05). The mean for each measure and each group is represented by a grey line.
Association between neurobehavioural changes and microstructural changes after MBI
The first PLSC analysis applied on neurobehavioural functioning (i.e., Δ neurobehavioural measures, age at the assessment and gestational age) and Δ mean FA identified one statistically significant latent component: latent component 1 (p = 0.005). Mean saliences as well as their bootstrap-estimated standard deviations for neurobehavioural functioning and Δ mean FA are reported in Supplementary Table S4. Latent component 1 revealed a general pattern of association between neurobehavioural functioning and Δ mean FA, see Figure 5. A decreased in Δ scores of the BRIEF GEC and BRIEF IM were associated with an increase in Δ mean FA values for a range of tracts, including the arcuate fascicle left/right, the anterior thalamic radiation left/right, the superior thalamic radiation right, the genu and isthmus of the corpus callosum, the cingulum left/right, the corticospinal tract left, the fronto-pontine tract left/right, the inferior cerebellar peduncle right, the inferior occipito-frontal fascicle left, the middle cerebellar peduncle, the optic radiation right, the superior cerebellar peduncle left/right, the superior longitudinal fascicle II and III left/right, the uncinate fascicle left, the thalamo-premotor and thalamo-occipital tract left/right, the striato-fronto-orbital left as well as the striato-premotor tract left/right. Moreover, lower gestational age was associated with a robust increased in Δ mean FA values in these same tracts, see Figure 5a. Importantly these microstructural changes were not associated with the age at the assessment of the participants.
Associations between neurobehavioural functioning and Δ mean FA based on the PLSC analysis. A) Neurobehavioural functioning saliences: the diverging graph show mean saliences averaged across bootstrap samples and their bootstrap-estimated standard deviations (x-axis) for each neurobehavioural functioning measure (y-axis); robust saliences are represented in blue. Of note, saliences bellow 0 indicate a decreased in scores after MBI and saliences above 0 indicate an increased in scores after MBI. B) Robust Δ mean FA saliences: the diverging graph show robust mean saliences averaged across bootstrap samples and their bootstrap-estimated standard deviations (x-axis) for each Δ mean FA along a given tract (y-axis). The tracts extracted from TractSeg and showing robust Δ mean FA saliences are shown on the right side and the colour code for principal direction of diffusion is used: red (left – right), green (anterior – posterior), blue (inferior – superior). Based on the results from the second PLSC analyses: pink star indicates a robust negative contribution of Δ mean ODI along the given tract associated with a similar pattern of neurobehavioural functioning saliences; and orange star indicates both a robust negative contribution of Δ mean ODI and a robust positive contribution of Δ mean ICVF along the given tract associated with a similar pattern of neurobehavioural functioning saliences. Of note, saliences below 0 indicate a decreased in scores after MBI and saliences above 0 indicate an increased in scores after MBI; L left, R right.
The second PLSC analysis applied on neurobehavioural functioning (i.e., Δ neurobehavioural measures, age at the assessment and gestational age) and Δ mean ICVF and ODI identified one statistically significant latent component: latent component 1 (p = 0.02). Mean saliences as well as their bootstrap-estimated standard deviations for neurobehavioural functioning and Δ mean ICVF and ODI are reported in Supplementary Table S5 and Supplementary Figure S1. Latent component 1 revealed an association between decreased in Δ scores of the BRIEF GEC and BRIEF IM with a general decrease in Δ mean ODI for all tracts showing a robust increase in Δ mean FA, including: arcuate fascicle left/right, the anterior thalamic radiation left/right, the superior thalamic radiation right, the genu of the corpus callosum, the cingulum left/right, the corticospinal tract left, the fronto-pontine tract left/right, the inferior cerebellar peduncle right, the inferior occipito-frontal fascicle left, the middle cerebellar peduncle, the optic radiation right, the superior cerebellar peduncle left/right, the superior longitudinal fascicle II left/right and III right, the uncinate fascicle left, the thalamo-premotor and thalamo-occipital tract left/right, the striato-fronto-orbital left as well as the striato-premotor tract left/right; with the exception of the isthmus of the corpus callosum and the superior longitudinal fascicle III left, see Figure 5. This decreased in Δ scores of the BRIEF GEC and BRIEF IM were also associated with a combined increased in Δ mean ICVF and decrease in Δ mean ODI for left-lateralised tracts including: the arcuate fascicle, inferior cerebellar peduncle, thalamo-parietal and thalamo-occipital. Moreover, the decreased in Δ scores of the BRIEF GEC and BRIEF IM were associated with a robust reduction of ODI only on the rostrum of the corpus callosum, on the superior longitudinal fascicle I left and on the striato-fronto-occipital right. Finally, the general increase in Δ ICVF and decrease in ODI measures were also associated with higher gestational age in VPT young adolescents. Importantly these microstructural changes were not associated with the age of the participants.
Discussion
In this study, we investigated the benefits of an 8-weeks MBI on neurobehavioral functioning and its association with white-matter microstructural changes in VPT young adolescents using both DTI and NODDI parameters. Our finding showed an enhancement of global executive functioning in daily life after MBI that was associated with a general pattern of significant increase in FA along with a decrease in ODI in a range of white-matter tracts involved in executive processes. This general pattern of increase FA and decrease ODI values was also negatively associated with gestational age, meaning that the microstructural changes in FA and ODI after MBI was particularly marked in young adolescents with lower gestational age.
In the present study, preterm young adolescents showed significant difficulties in executive and behavioural functioning in daily life as measured by the BRIEF and SDQ parent questionnaires, as well as significantly lower self-compassion using a self-reported questionnaire compared to full-term controls. These weaknesses have been previously described in several cohorts of children, adolescents and adults born prematurely (Aarnoudse-Moens, Weisglas-Kuperus, van Goudoever, & Oosterlaan, 2009; Luu, Ment, Allan, Schneider, & Vohr, 2011; Mulder, Pitchford, Hagger, & Marlow, 2009). After an 8-weeks MBI, an enhancement of executive and behavioural functioning in daily life of preterm young adolescents was observed. These behavioural results corroborate the results from our previous study using a gold standard randomised controlled trial design (Siffredi, Liverani, Hüppi, et al., 2021), which indicates that the gain observed after MBI on these executive and behavioural measures is independent of test-retest effect. The benefits of MBI observed on daily executive and behavioural functioning are also consistent with previous studies completed in other populations (Flook, Goldberg, Pinger, & Davidson, 2015; Flook et al., 2010; Geronimi et al., 2019; van de Weijer-Bergsma, Formsma, de Bruin, & Bogels, 2012; van der Oord, Bogels, & Peijnenburg, 2012).
Combining both the classical DTI model and the NODDI model (i.e., a multi-component model that captures neurite morphology), microstructural changes associated with the beneficial effect of MBI on executive and behavioural functioning were explored. Our finding show that the gain observed in executive functioning after an 8-week MBI was associated with a general pattern of increase in FA along with a decrease in ODI. While an increase in FA values reflect changes in the diffusion of water molecules along a given tract, the decrease in axonal dispersion indicate more specifically an enhancement in the coherence of the spatial organisation of the axons (Mah, Geeraert, & Lebel, 2017). This pattern of higher FA coupled with decrease ODI in a specific white-matter tract is commonly interpreted as an indication of increased connectivity (Laneri et al., 2016). Importantly, this pattern of microstructural changes associated with improved executive functioning after MBI was observed in white-matter tracts known to be involved in different executive processes, i.e., attentional control, inhibition, processing speed, updating, fluency and self-regulation (Bettcher et al., 2016; Diao et al., 2015; Frye et al., 2010; Mamiya, Richards, & Kuhl, 2018) or in white-matter tracts implicated in sensorimotor circuits that support executive functioning (Koziol & Lutz, 2013), including the anterior and superior thalamic radiation, the genu of the corpus callosum, the cingulum, the fronto-pontine tract, the inferior occipito-frontal fascicle, the superior longitudinal fascicle II and III, the arcuate fascicle, the thalamo-premotor and thalamo-occipital tract, the striato-fronto-orbital as well as the striato-premotor tract. Moreover, the uncinate fascicle, known to be involved in emotional-regulation processes, also showed this pattern of increase FA and reduced ODI associated with improvement in executive functioning after MBI (Oishi et al., 2015; Sobhani, Baker, Martins, Tuvblad, & Aziz-Zadeh, 2015). These results corroborate previous research studies showing that the practice of MBI or MBI-related approaches such as meditation are associated with increases in FA in these tracts or related-regions of interest (Hölzel et al., 2016; Kang et al., 2013; Laneri et al., 2016; Luders et al., 2011; Luders et al., 2012; Patsenko et al., 2019; Tang et al., 2012). Moreover, there was an association between changes in executive functioning after MBI with changes in cerebellar mean FA including, the inferior cerebellar peduncle, the middle cerebellar peduncle and the superior cerebellar peduncles. Again, these FA changes were more specifically coupled with a decrease in neurite orientation dispersion. Cerebellar regions have been implicated in executive functions and proposed to play a role of “control of behaviour” by linking movement to though (Berquin et al., 1998; Koziol, Budding, & Chidekel, 2012; Schweizer et al., 2008). This is also in line with previous results showing changes in grey matter density of cerebellar regions after an 8-week Mindfulness program (Hölzel et al., 2011). Compared to previous findings, the use of the NODDI model in the current study allows to characterise more precisely this increase in FA explained by a concurrent decrease in neurite orientation dispersion.
Altogether, the enhancement in executive functioning was not only associated with FA changes in white-matter tracts involved in executive functioning per se; but it was also associated with FA changes in white-matter tracts implicated in bottom-up and sensorimotor interaction on which executive functions rely. Importantly, the NODDI model allowed to characterise more precisely this general increase in FA with concurrent decrease in neurite orientation dispersion. As opposed to executive functioning, the changes observed in socio-emotional and behavioural functioning after MBI, as measured by the SDQ total score, were not significantly associated with changes in mean FA.
Additionally, the increase in FA and reduction in ODI values following MBI of these tracts was associated with lower gestational age at birth. In other words, the microstructural gain observed on the tracts mentioned above was bigger in the most preterm young adolescents. It is likely that MBI is most beneficial for the most vulnerable young adolescents as they have more opportunities for further enhancements. This is consistent with the results of our randomised controlled trial showing individual variability in the response to MBI with a greater benefit on overall executive skills in VPT young adolescents born with higher risk (i.e., smaller gestational age and lower weight at birth), in comparison with VPT young adolescents born with lower risk (Siffredi, Liverani, Hüppi, et al., 2021). This pattern of increased benefit of MBI in the most vulnerable individuals has already been observed in children with weaker executive abilities (Flook et al., 2010).
Finally, the reduction in executive difficulties observed after an 8-week MBI was also associated with an increase in neurite density (ICVF) in white-matter tracts implicated in executive functioning, including the arcuate fasciculus as well as thalamic and cerebellar tracts. This pattern reflecting an increase in packing density of axons in white matter in a few tracts was also associated with reduced gestational age at birth. These results go in line with the suggested increased benefit of MBI in the most vulnerable individuals and in individuals with weaker executive abilities (Flook et al., 2010; Siffredi, Liverani, Hüppi, et al., 2021).
Our findings should also be considered in light of the limitations of the current study. The beneficial effect of MBI on executive and socio-emotional competences were observed through parent-questionnaires. These changes observed via parents-reported questionnaires might be questionable as parents were not blind to the intervention. This point might be methodologically improved via the inclusion of an active control condition or placebo condition, as done in previous studies (Britton et al., 2014; Pascoe et al., 2013). Nevertheless, the association found with micro-structural changes on specific tracks suggest that the benefit of MBI observed via parent questionnaire are well-grounded.
Conclusions
This study provides evidence that the enhancement in executive functioning after an MBI in VPT young adolescents is associated with white-matter microstructural changes in tracts involved in executive processes as well as in tracts involved in sensorimotor circuit supporting and interacting with executive processes. Our study also suggests that the most vulnerable VPT young adolescents, i.e., the ones with lowest gestational age at birth, show the biggest gain in terms of white-matter microstructural changes. Finally, MBI appears to be a promising tool for enhancing executive functioning and white-matter brain plasticity in a vulnerable population such as VPT young adolescents.
Data Availability
Ethical restrictions prevent us from making anonymised data available in a public repository. Data may be available for researchers to researchers who meet the criteria for access to confidential data by direct request to the "Mindful preterm teens study" team: Russia.HaVinhLeuchter{at}unige.ch. There are restrictions on data related to identifying participant information and appropriate ethical approval is required prior to release. Only de-identified data will be available.
Funding
This work was supported by the Swiss National Science Foundation, No. 324730_163084 [PI: P.S. Hüppi].
Data availability statement
Ethical restrictions prevent us from making anonymised data available in a public repository. Data may be available for researchers to researchers who meet the criteria for access to confidential data by direct request to the “Mindful preterm teens study” team: Russia.HaVinhLeuchter{at}unige.ch. There are restrictions on data related to identifying participant information and appropriate ethical approval is required prior to release. Only de-identified data will be available.
Acknowledgements
We thank and acknowledge all participating young adolescents and families who made this research possible. We also thank the Fondation Campus Biotech Geneva (FCBG), a foundation of the Swiss Federal Institute of Technology Lausanne (EPFL), the University of Geneva (UniGe), and the University Hospitals of Geneva (HUG); the Research Platform of the University Hospitals of Geneva (HUG) for their practical help; as well as Mariana Magnus Smith and Françoise Stuckelberger-Grobéty for their implication as MBI instructors.
Abbreviations
- NEPSY-II
- A Developmental Neuropsychological Assessment, 2nd Edition
- AD
- Axial Diffusivity
- BRIEF
- Behaviour Rating Inventory of Executive Function-parent version
- BRI
- Behavioural Regulation Index
- DTI
- Diffusion Tensor Imaging
- DW
- Diffusion-Weighted
- FDR
- False Discovery Rate
- FOD
- Fibre Orientation Distributions
- fODF
- Fibre Orientation Distribution Function
- FA
- Fractional Anisotropy
- GAI
- General Ability Index
- GEC
- Global Executive Composite
- ICVF
- Intra-cellular Volume Fraction
- LCs
- Latent Components
- MRI
- Magnetic Resonance Imaging
- MI
- Mean Diffusivity
- MD
- Metacognition Index
- MBCT
- Mindfulness-Based Cognitive Therapy
- MBI
- Mindfulness-Based Intervention
- MBSR
- Mindfulness-Based Stress Reduction
- MSMT-CSD
- Multi-Shell Multi-Tissue Constrained Spherical Deconvolution
- NODDI
- Neurite Orientation Dispersion and Density Imaging
- (ODI)
- Orientation Dispersion Index
- PLSC
- Partial Least Square Correlation Analyses
- RD
- Radial Diffusivity
- SDQ
- Strength and Difficulties Questionnaire
- TCC
- Temporal Context Confusion index
- TractSeg
- Tractography-Based Segmentation
- VPT
- Very Preterm
- WISC-IV
- Wechsler Intelligence Scale for Children, 4th Edition
References
