Functional and microstructural plasticity following social and interoceptive mental training

  1. Sofie Louise Valk  Is a corresponding author
  2. Philipp Kanske
  3. Bo-yong Park
  4. Seok-Jun Hong
  5. Anne Böckler
  6. Fynn-Mathis Trautwein
  7. Boris C Bernhardt
  8. Tania Singer
  1. Otto Hahn Group Cognitive Neurogenetics, Max Planck Institute for Human Cognitive and Brain Sciences, Germany
  2. INM-7, FZ Jülich, Germany
  3. Clinical Psychology and Behavioral Neuroscience, Faculty of Psychology, Technische Universität Dresden, Germany
  4. Max Planck Institute for Human Cognitive and Brain Sciences, Germany
  5. Multimodal Imaging and Connectome Analysis Lab, McConnell Brain Imaging Centre, Montreal Neurological Institute and Hospital, McGill University, Canada
  6. Department of Data Science, Inha University, Republic of Korea
  7. Center for Neuroscience Imaging Research, Institute for Basic Science, Republic of Korea
  8. Center for the Developing Brain, Child Mind Institute, United States
  9. Department of Biomedical Engineering, Sungkyunkwan University, Republic of Korea
  10. Department of Psychology, Wurzburg University, Germany
  11. Department of Psychosomatic Medicine and Psychotherapy, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Germany
  12. Social Neuroscience Lab, Max Planck Society, Germany
5 figures, 6 tables and 2 additional files

Figures

Study design.

(A) Training design of the ReSource study; (B) Training modules; (C) Task-based meta-analytical maps, and a legend of the color-coding of the maps; (D) Functional cortical organization: i. functional connectivity matrix, gradient 1–3, eccentricity metric; ii. task-based network embedding; iii. intracortical microstructure, mean qT1 values as a function of task-based meta-analytical maps and cortical depth and relative values (z-scored per depth-compartment).

Figure 2 with 4 supplements
Training-induced changes in cortical functional organization.

(A) upper: T-maps of TM-specific changes in functional eccentricity; lower: TM-specific change in functional eccentricity, p<0.01, FDRq <0.05 outlined in black, below: alterations of eccentricity in the FDRq <0.05 regions, right: mean changes in FDRq <0.05 eccentricity regions as a function of G1-G2-G3; (B) A-priori network functional eccentricity change in networks that showed TM-specific change.

Figure 2—figure supplement 1
TM-specific change in functional eccentricity as a function of training cohort and a-priori network.

X-axes depict days of training, y-axes change in functional eccentricity, colors reflect the TM (yellow = Presence; red = Affect; green = Perspective).

Figure 2—figure supplement 2
TM-specific change in functional Gradient 1.

Trends at p<0.01, FDRq <0.05 outlined in black.

Figure 2—figure supplement 3
TM-specific change in functional Gradient 2.

Trends at p<0.01, FDRq <0.05 outlined in black.

Figure 2—figure supplement 4
TM-specific change in functional Gradient 3.

Trends at p<0.01, FDRq <0.05 outlined in black.

Figure 3 with 2 supplements
Changes in depth-varying microstructure as a function of TM.

(A). Change in cortical microstructure, per TM, red indicates positive change in qT1, blue negative change. FDRq <0.05 findings are outlined in black on top of t-values per parcel; (B) TM specific change in cortical microstructure. Red indicates positive change in qT1, blue negative change. FDRq <0.05 findings are outlined in black in combination with semi-transparent trends (p<0.01).

Figure 3—figure supplement 1
Training-specific change in qT1 from baseline to T3 as a function superficial (1:4) mid (5:8) and deep (8:12) depth compartment.

(A). Left: Change in retest controls, right: Change in training versus retest controls, trends at p<0.01, FDRq <0.05 outlined in black. (B) Network-specific change in qT1 as a function of training cohort.

Figure 3—figure supplement 2
Cohort-specific change in qT1 from baseline to T1 as a function superficial (1:4) mid (5:8) and deep (8:12) depth compartment.

(A). First row: Left: Change in Presence TC1 vs retest controls, middle: change in Presence TC1 versus Affect TC3, right: Change in Affect TC3 versus retest controls; second row: Left: Change in Presence TC2 vs retest controls, right: change in Presence TC2 versus Affect TC3, trends at p<0.01, FDRq <0.05 outlined in black; (B) Network-specific change in qT1 as a function of training cohort in T0-T1.

Figure 4 with 2 supplements
Dissociable microstructural alterations following mental training.

(A).TM-specific changes in cortical microstructure; qT1 in regions showing eccentricity change (y-axis: depth, x-axis: qT1 change); (B) Network-specific change in cortical microstructure as a function of depth, mean change per TM, pFDR <0.05 have black outline (y-axis: depth, x-axis: qT1 change). The boxes on the right of each plot display the statistics (t-values) of the respective difference between TM, with the contrast color coded as upper minus lower TM (defined by color); (C) Correspondence of functional versus microstructural change; i. Spatial correlation of mean alterations in each TM, black outline indicates pspin <0.05, as a function of cortical depth.

Figure 4—figure supplement 1
TM-specific change in superficial (1:4) mid (5:8) and deep (8:12) depth compartment microstructure as a function of training cohort and a-priori network.

Yellow = Presence, Red = Affect, Green = Perspective, Blue = Retest Control. Changes are displayed as a function of days past baseline (MRI measurement points).

Figure 4—figure supplement 2
TM-specific change in qT1 as a function superficial (1:4) mid (5:8) and deep (8:12) depth compartment.

Network specific changes in microstructure as a function of depth per training cohort/retest control cohort, and time point.

Behavioral change prediction.

(A).Attention change; from left to right: attention task, behavioral change, predicting weights, nMAE and holdout r distribution; (B). Compassion change; from left to right: compassion task, behavioral change, predicting weights, nMAE and holdout r distribution; (C). Perspective-taking change; from left to right: perspective-taking task, behavioral change, predicting weights, nMAE and holdout r distribution (left side; within TM, right side; other TM (yellow: Presence; red: Affect; Green: Perspective) or RCC (grey)).

Tables

Table 1
Participant inclusion in resting-state analysis and quantitative T1 analysis.
Recruited (N, mean age, % female)T0 (N)T1 (N)T2 (N)T3 (N)
Total (N=332)
TC1 (N=80; 41.3; 58.8)
TC2 (N=81; 41.2; 59.3)
RCC (N=90; 40.0; 58.9)
TC3 (N=81; 40.4; 60.5)
268
69
67
65
67
259
65
59
68
67
182
57
61
64
184
55
64
65
Table 2
Reason for missing data across the study duration.

MR incidental findings are based on T0 radiological evaluations; participants who did not survive MRI quality control refers to movement and/or artefacts in the T1-weighted MRI; dropout details can be found in Singer et al., 2016; no MRT: due to illness / scheduling issues / discomfort in scanner; other: non-disclosed; functional MRI missing: no complete functional MRI; functional MRI quality:>0.3 mm movement (low quality in volume +surface).

Reason for dropout(TC1, TC2, RCC)T0T1T2T3
Structural MR incidental finding
Structural MRI quality control
Dropout
Medical reasons
Other
Functional MRI missing/QC qT1 missing
5
7
2
1
4
18
13
(5 based on T0)
6
7 (2 based on T0)
7 (1 based on T0)
10
14
7
(5 based on T0)
4
7 (9 based on T01)
7 (8 based on T01)
7
16
6
(5 based on T0)
2
7 (16 based on T012)
(15 based on T012)
7
8
7
Table 3
Reason for missing data across the study duration.

MR incidental findings are based on T0 radiological evaluations; participants who did not survive MRI quality control refers to movement and/or artefacts in the T1-weighted MRI; dropout details can be found in Singer et al., 2016; no MRT: due to illness / scheduling issues / discomfort in scanner; other: non-disclosed.

Reason for dropout (TC3)T0T1
MR incidental finding
MRI quality control
Dropout
Medical reasons
Other
Functional MRI missing/QC qT1 missing
3
0
0
1
5
1
4
(3 based on T0)
0
3
2
3
3
0
Table 4
TM-specific changes in eccentricity per a priori networks.
Presence (n=109) vs Perspective (n=96)Presence (n=109) vs Affect (n=104)Perspective (n=96)vs Affect (n=104)
Attentiont=2.842, p=0.005*, d=0.247t=1.458, p>0.05, d=0.127t=−1.692, p>0.05, d=−0.147
Interoceptiont=2.765, p=0.006*, d=0.240t=1.043, p>0.05, d=0.091t=−2.008, p=0.045, d=−0.174
Emotiont=0.387, p>0.05, d=0.035t=−0.135, p>0.05, d=−0.011t=−0.552, p>0.05, d=−0.048
Empathyt=2.218, p=0.027, d=193t=0.879, p>0.05, d=0.076t=−1.569, p>0.05, d=−0.136
Theory of Mindt=1.721, p>0.05, d=0.149t=1.324, p>0.05, d=0.115t=−0.601, p>0.05, d=−0.052
  1. *

    signifies FDR corrected differences.

Author response table 1
Eccentricity and non-overlapping networks (Figure 2).
Presence-PerspectivePresence-AffectPerspective-Affect
attention2,776p=0,0061,451p=0,147-1,624p=0,105
interoception2,637p=0,0091,279p=0,201-1,641p=0,101
emotion0,034p=0,973-0,353p=0,724-0,377p=0,706
empathy1,700p=0,0900,510p=0,610-1,360p=0,174
ToM1,901p=0,0581,172p=0,242-0,943p=0,346
Author response table 2
Microstructure and non-overlapping networks as a function of cortical depth in qT1 (Figure 4).
Presence-Perspective
attentioninteroceptionemotionempathyToM
t-valuep-valuet-value p-valuet-valuep-valuet-valuep-valuet-valuep-value
5,6800,0003,6730,000-0,7630,4463,7370,0004,9950,000
5,5930,0004,4540,000-0,8240,4103,8920,0004,7530,000
5,1850,0004,7430,000-0,9850,3253,8260,0004,0860,000
4,5650,0004,5820,000-1,2110,2263,6010,0003,1710,002
3,8450,0004,1350,000-1,4440,1493,3060,0012,3310,020
3,1190,0023,5090,000-1,6050,1093,0200,0031,7410,082
2,4360,0152,7700,006-1,6560,0982,6900,0071,4400,150
1,8520,0651,9200,055-1,6580,0982,3200,0211,3210,187
1,4340,1521,0380,300-1,6960,0901,9520,0511,2680,205
1,2150,2250,2240,823-1,7960,0731,6830,0931,2450,214
1,1610,246-0,4460,656-1,8990,0581,5600,1191,2340,218
1,2160,225-0,9490,343-2,0180,0441,5460,1231,2140,225
Presence-Affect
3,8680,0002,6850,0070,9130,3623,0550,0023,1330,002
4,1030,0003,1790,0020,9600,3383,2570,0013,0360,003
4,0910,0003,4200,0010,9150,3603,3610,0012,8200,005
3,9110,0003,4470,0010,8280,4083,4030,0012,5600,011
3,6560,0003,3130,0010,8090,4193,4320,0012,3950,017
3,3860,0013,0900,0020,8920,3733,5060,0002,3720,018
3,1180,0022,8470,0051,0410,2983,5880,0002,4180,016
2,8670,0042,6180,0091,2010,2303,5520,0002,4340,015
2,6630,0082,4040,0171,2680,2053,3630,0012,3810,018
2,5520,0112,1810,0301,2510,2113,1310,0022,2950,022
2,5260,0122,0170,0441,1690,2432,9570,0032,2120,027
2,5550,0111,8870,0600,9920,3222,8760,0042,1760,030
Perspective-Affect
-2,4470,015-1,4060,1601,7090,088-1,1240,262-2,4180,016
-2,1270,034-1,7780,0761,8210,069-1,1010,271-2,2490,025
-1,6970,090-1,8580,0641,9530,051-0,9320,352-1,7330,084
-1,1960,232-1,6580,0982,1150,035-0,6480,517-0,9890,323
-0,6600,509-1,3010,1942,3500,019-0,2990,765-0,2330,816
-0,1320,895-0,8370,4032,6060,0090,0840,9330,3860,699
0,3520,725-0,2700,7872,8050,0050,5240,6000,7590,448
0,7440,4570,4320,6662,9570,0030,8940,3720,9040,367
1,0020,3171,1820,2383,0620,0021,1110,2670,9100,363
1,1330,2581,8540,0643,1530,0021,1820,2380,8520,394
1,1660,2442,4260,0163,1860,0021,1480,2520,7840,433
1,1340,2572,8490,0053,1460,0021,0850,2780,7720,441

Additional files

Supplementary file 1

Complete summary of supplementary tables and analysis.

(a). Descriptive statistics Presence (b). Descriptive statistics Affect (c). Descriptive statistics Affect, excluding active controls (TC3) (d). Descriptive statistics Perspective (e). Descriptive statistics Retest controls. (f). Functional eccentricity changes GSR controlled per a-priori network. T-values and p-values below P<0.05, * indicates FDRp <0.05. (g). Functional eccentricity changes in training cohort 1 per a-priori network. T-values and p-values below P<0.05, * indicates FDRp <0.05. (h). Functional eccentricity changes in training cohort 2 per a-priori network. T-values and p-values below P<0.05, * indicates FDRp <0.05. (i). Functional eccentricity changes per a-priori network baseline to T1. T-values and p-values below P<0.05, * indicates FDRp <0.05. (j). Functional eccentricity changes per a-priori network T1 to T3. T-values and p-values below P<0.05, * indicates FDRp <0.05. (k). G1-G3 change per a-priori network Presence vs Perspective. T-values and p-values below P<0.05, * indicates FDRp <0.05. (l). G1-G3 change per a-priori network Presence vs Affect. T-values and p-values below P<0.05, * indicates FDRp <0.05. (m). G1-G3 change per a-priori network Perspective vs Affect. T-values and p-values below P<0.05, * indicates FDRp <0.05. (n). Functional eccentricity changes per a-priori network controlling for cortical thickness change. T-values and p-values below P<0.05, * indicates FDRp <0.05. (o). Functional eccentricity changes per a-priori network from baseline to T3. T-values and p-values below P<0.05, * indicates FDRp <0.05. (p). Depth-dependent qT1 change per a-priori network Training vs Retest Control. T-values and p-values below P<0.05. * indicates FDRp <0.05.(q). Descriptives of retest-control change (mean change over T0-T1; T1-T2; T2-T3) as a function of depth-dependent qT1. (r). Descriptive of Presence change (mean change over T0-T1) as a function of depth-dependent qT1 (s). Descriptives of Affect change (mean change over T0-T1, T1-T2 and T2-T3) as a function of depth-dependent qT1. (t). Descriptives of Affect change (mean change over T1-T2 and T2-T3) as a function of depth-dependent qT1. (u). Descriptives of Perspective change (mean change over T1-T2 and T2-T3) as a function of depth-dependent qT1. (v). Depth-dependent qT1 change per a-priori network Presence vs Perspective. T-values and p-values below P<0.05, and Cohen’s D effect size, * indicates FDRp <0.05. (w). Depth-dependent qT1 change per a-priori network Presence vs Affect. T-values and p-values below P<0.05, and Cohen’s D effect size, * indicates FDRp <0.05. (x). Depth-dependent qT1 change per a-priori network Perspective vs Affect. T-values and p-values below P<0.05, and Cohen’s D effect size, * indicates FDRp <0.05. (y). Depth-dependent qT1 change per a-priori network Presence vs Perspective in TC1. T-values and p-values below P<0.05. * indicates FDRp <0.05. (z). Depth-dependent qT1 change per a-priori network Presence vs Affect in TC1. T-values and p-values below P<0.05. * indicates FDRp <0.05. (za). Depth-dependent qT1 change per a-priori network Perspective vs Affect in TC1. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zb). Depth-dependent qT1 change per a-priori network Presence vs Perspective in TC2. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zc). Depth-dependent qT1 change per a-priori network Presence vs Affect in TC2. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zd). Depth-dependent qT1 change per a-priori network Perspective vs Affect in TC2. T-values and p-values below P<0.05. * indicates FDRp <0.05. (ze). Depth-dependent qT1 change per a-priori network baseline – T1: TC1 (Presence) versus Retest Control. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zf). Depth-dependent qT1 change per a-priori network baseline – T1: TC2 (Presence) versus Retest Control. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zg). Depth-dependent qT1 change per a-priori network baseline – T1: Affect TC3 vs Retest Control. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zh). Depth-dependent qT1 change per a-priori network baseline – T1: Presence TC1 vs Affect TC3. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zi). Depth-dependent qT1 change per a-priori network baseline – T1: Presence TC2 vs Affect TC3. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zj). Depth-dependent qT1 change per a-priori network T1 – T3: Perspective versus Affect (TC1 +TC2). T-values and p-values below P<0.05. * indicates FDRp <0.05. (zk). Depth-dependent qT1 change per a-priori network T1 – T3: Affect vs Retest Control. * T-values and p-values below P<0.05. * indicates FDRp <0.05. (zl). Depth-dependent qT1 change per a-priori network T1 – T3: Perspective vs Retest Control. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zm). Depth-dependent qT1 change per a-priori network controlling for CTX: Presence versus Perspective. T-values and p-values below P<0.05. * indicates FDRp <0.05.(zn). Depth-dependent qT1 change per a-priori network controlling for CTX: Presence versus Affect. T-values and p-values below P<0.05. * indicates FDRp <0.05. (zo). Depth-dependent qT1 change per a-priori network controlling for CTX: Perspective versus Affect. T-values and p-values below P<0.05. * indicates FDRp <0.05.

https://cdn.elifesciences.org/articles/85188/elife-85188-supp1-v2.docx
MDAR checklist
https://cdn.elifesciences.org/articles/85188/elife-85188-mdarchecklist1-v2.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Sofie Louise Valk
  2. Philipp Kanske
  3. Bo-yong Park
  4. Seok-Jun Hong
  5. Anne Böckler
  6. Fynn-Mathis Trautwein
  7. Boris C Bernhardt
  8. Tania Singer
(2023)
Functional and microstructural plasticity following social and interoceptive mental training
eLife 12:e85188.
https://doi.org/10.7554/eLife.85188