Reevaluating the neural noise in dyslexia using biomarkers from electroencephalography and high-resolution magnetic resonance spectroscopy

  1. Agnieszka Glica
  2. Katarzyna Wasilewska
  3. Julia Jurkowska
  4. Jarosław Żygierewicz
  5. Bartosz Kossowski
  6. Katarzyna Jednoróg  Is a corresponding author
  1. Laboratory of Language Neurobiology, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Poland
  2. Faculty of Physics, University of Warsaw, Poland
  3. Laboratory of Brain Imaging, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Poland
4 figures, 11 tables and 1 additional file

Figures

Selected results for the electroencephalography (EEG) excitatory and inhibitory (E/I) balance biomarkers.

n = 119 (DYS n = 59, CON n = 60). (A) Power spectral densities averaged across three electrodes (T7, TP7, TP9)corresponding to the left superior temporal sulcus (STS) separately for dyslexic (DYS) and control (CON) groups at rest and (C) during the language task. (B) Plots illustrating results for the exponent, offset, and the beta power from the left STS electrodes at rest and (D) during the language task.

Figure 2 with 1 supplement
Selected results for the magnetic resonance spectroscopy (MRS) excitatory and inhibitory (E/I) balance biomarkers.

(A) Sample MRS voxel localization. (B) Sample MRS spectrum with individual total creatine (tCr), Glu, and GABA+ contributions plotted below. (C) Group results (CON>DYS) from the fMRI localizer task for words compared to the control stimuli (p<0.05 FWE cluster threshold) and overlap of the MRS voxel placement across participants. n = 50 (DYS n = 21, CON n = 29). (D) MRS spectra separately for the dyslexic (DYS) and control (CON) groups. n = 50 (DYS n = 21, CON n = 29). (E) Plots illustrating results for the Glu, n = 50 (DYS n = 21, CON n = 29), GABA+, n = 47 (DYS n = 20, CON n = 27), Glu/GABA+ ratio, n = 47 (DYS n = 20, CON n = 27), and the Glu/GABA+ imbalance, n = 47 (DYS n = 20, CON n = 27).

Figure 2—figure supplement 1
Main effect of group CON>DYS for both visual (words > false fonts) and auditory runs (words > words backward) from the fMRI localizer task.

Note: CON – control group, DYS – dyslexic group, L – left hemisphere, R – right hemisphere. Results reported at p<0.001 height threshold corrected for multiple comparisons using p<0.05 FWE cluster threshold. n = 50 (DYS n = 21, CON n = 29).

Relationships between electroencephalography (EEG) and magnetic resonance spectroscopy (MRS) excitatory and inhibitory (E/I) balance biomarkers.

(A) Semi-partial correlation between offset at rest in the left superior temporal sulcus (STS) electrodes and Glu controlling for age and gray matter volume (GMV), n = 50 (DYS n = 21, CON n = 29). (B) Semi-partial correlation between exponent at rest in the left STS electrodes (controlled for age) and Glu/GABA+ imbalance (controlled for GMV), n = 47 (DYS n = 20, CON n = 27).

Associations between beta power, phonological awareness, and reading.

n = 119 (DYS n = 59, CON n = 60). (A) Semi-partial correlation between phonological awareness controlling for age and beta power (in the left superior temporal sulcus [STS] electrodes) at rest and (B) during the language task. (C) Partial correlation between phonological awareness and reading controlling for age. (D) Mediation analysis results. Unstandardized b regression coefficients are presented. Age was included in the analysis as a covariate. 95% CI – 95% confidence intervals; left STS – values averaged across three electrodes corresponding to the left superior temporal sulcus (T7, TP7, TP9).

Tables

Table 1
Demographic and behavioral characteristics of the entire sample of 120 participants.
DYS28 F, 32 MCON28 F, 32 Mt(df)pCohen’s dBF10
MSDMSD
Demographics
Age19.413.1819.542.960.25
(118)
0.8060.050.20
Mother’s education
(years)
17.203.3616.582.28–1.19
(103.89)
0.235–0.220.37
Father’s education
(years)
16.12a3.10a17.13a3.27a1.71
(114)
0.0910.320.73
IQ103.56b11.83b111.1210.433.70
(117)
<0.0010.6875.31
Nonverbal IQ
(scaled score)
10.402.9411.622.572.42
(118)
0.0170.442.63
ARHQ-PL51.509.7025.478.00–16.04
(113.87)
<0.001–2.93>10,000
Reading and reading-related tasks
Words/min108.3820.93134.5713.298.18
(99.90)
<0.001*1.49>10,000
Pseudowords/min56.7514.1683.4317.049.33
(118)
<0.001*1.70>10,000
RAN objects (s)32.125.1128.704.43–3.92
(118)
<0.001*–0.72149.93
RAN colors (s)35.836.8231.185.73–4.04
(118)
<0.001*–0.74229.96
RAN digits(s)19.324.6116.252.94–4.34
(100.28)
<0.001*–0.79642.86
RAN letters (s)22.704.5319.683.16–4.23
(105.42)
<0.001*–0.77433.23
Reading comprehension (s)64.4720.1343.729.63–7.20
(84.66)
<0.001*–1.32>10,000
Phoneme deletion
(% correct)
76.4124.6891.479.074.44
(74.66)
<0.001*0.81898.25
Spoonerisms phonemes
(% correct)
54.2935.4282.7422.065.28
(98.78)
<0.001*0.96>10,000
Spoonerisms syllables
(% correct)
46.9430.6173.0623.985.20
(111.62)
<0.001*0.95>10,000
Orthographic awareness
(accuracy/time)
0.330.130.530.148.12
(118)
<0.001*1.48>10,000
Perception speed (sten score)3.322.044.501.673.48
(118)
<0.001*0.6438.71
Digits forward5.531.646.981.954.40
(118)
<0.001*0.80792.55
Digits backward5.251.497.332.255.99
(102.59)
<0.001*1.09>10,000
  1. Note: DYS – dyslexic group; CON – control group; F – females, M – males. BF10 – Bayes factor indicating ratio of the likelihood of an alternative hypothesis (H1) to a null hypothesis (H0). ARHQ-PL – Polish version of the Adult Reading History Questionnaire. RAN – rapid automatized naming. Boldface indicates statistical significance at p<0.05 level (uncorrected).

  2. *Significance after Bonferroni correction for 14 planned comparisons for reading and reading-related tasks

  3. an = 58 (two participants did not provide information about the father’s education)

  4. bn = 59 (one participant did not attempt a verbal subtest of the scale, thus we were not able to calculate overall IQ)

Table 2
Semi-partial correlations between magnetic resonance spectroscopy (MRS) and electroencephalography (EEG) biomarkers of excitatory-inhibitory balance.

For Glu the effects of age and gray matter volume (GMV) were regressed, for GABA+, Glu/GABA+ ratio, and Glu/GABA+ imbalance the effect of GMV was regressed, while for exponents and offsets the effect of age was regressed.

Variable1.
r(BF10)
2.3.4.5.6.7.8.
EEG resting state
1. Glu
2. GABA+0.32*a
(1.82)
3. Glu/GABA+ ratio–0.08a
(0.21)
–0.91***a
(>10,000)
4. Glu/GABA+ imbalance0.12a
(0.25)
0.31*a
(1.63)
–0.18a
(0.37)
5. Exponent mean (rest)–0.03b
(0.18)
0.04a
(0.19)
–0.11a
(0.23)
0.29*a
(1.21)
6. Offset mean (rest)–0.21b
(0.49)
0.08a
(0.21)
–0.17a
(0.35)
0.25a
(0.69)
0.68***c
(>10,000)
7. Exponent
left STS (rest)
–0.16b
(0.32)
0.01a
(0.18)
–0.07a
(0.20)
0.35*a
(2.87)
0.68***c
(>10,000)
0.45***c
(>10,000)
8. Offset
left STS (rest)
–0.38**b
(6.28)
–0.10a
(0.23)
0.02a
(0.18)
0.18a
(0.38)
0.18*c
(0.80)
0.47***c
(>10,000)
0.66***c
(>10,000)
9. Beta power
left STS (rest)
–0.12b
(0.25)
0.17a
(0.33)
–0.25a
(0.74)
0.01a
(0.18)
0.18*c
(0.80)
0.21*c
(1.40)
0.51***c
(>10,000)
0.59***c
(>10,000)
EEG language task
5. Exponent mean (task)–0.10b
(0.22)
0.06a
(0.20)
–0.15a
(0.29)
0.21a
(0.49)
6. Offset mean (task)–0.26b
(0.92)
0.09a
(0.22)
–0.20a
(0.43)
0.22a
(0.53)
0.72***c
(>10,000)
7. Exponent
left STS (task)
–0.20b
(0.44)
0.01a
(0.18)
–0.09a
(0.22)
0.22a
(0.51)
0.65***c
(>10000)
0.51***c
(>10,000)
8. Offset
left STS (task)
–0.37**b
(5.05)
–0.08a
(0.21)
0.00a
(0.18)
0.15a
(0.30)
0.29**c
(18.31)
0.55***c
(>10,000)
0.77***c
(>10,000)
9. Beta power
left STS (task)
–0.19b
(0.42)
0.16a
(0.32)
–0.22a
(0.55)
–0.03a
(0.19)
0.05c
(0.13)
0.13c
(0.31)
0.47***c
(>10,000)
0.60***c
(>10,000)
  1. Note: r – Pearson’s correlation coefficient; BF10 – Bayes factor indicating ratio of the likelihood of an alternative hypothesis (H1) to a null hypothesis (H0); mean – values averaged across all electrodes; left STS – values averaged across three electrodes corresponding to the left superior temporal sulcus (T7, TP7, TP9).

  2. ***p < 0.001 (uncorrected); **p < 0.01 (uncorrected); *p < 0.05 (uncorrected)

  3. an = 47; bn = 50; cn = 119

Table 3
Semi-partial correlations between reading, phonological awareness, rapid automatized naming (RAN), multisensory integration, and biomarkers of excitatory-inhibitory balance.

For reading, phonological awareness, and multisensory integration the effect of age was regressed, for Glu the effects of age and gray matter volume (GMV) were regressed, for GABA+, Glu/GABA+ ratio, and Glu/GABA+ imbalance the effect of GMV was regressed, while for exponents and offsets the effect of age was regressed.

Variable1.
r(BF10)
2.3.4.
EEG resting state
1. Reading
2. Phonological awareness0.60***c
(>10,000)
3. RAN0.71***c
(>10,000)
0.48***c
(>10,000)
4. Multisensory integration0.16d
(0.41)
0.25*d
(2.09)
0.02d
(0.14)
5. Glu–0.03b
(0.18)
–0.12b
(0.25)
–0.05b
(0.19)
0.03b
(0.18)
6. GABA+–0.18a
(0.36)
–0.06a
(0.20)
–0.23a
(0.56)
0.31*a
(1.62)
7. Glu/GABA+ ratio0.07a
(0.20)
0.04a
(0.19)
0.09a
(0.22)
–0.32*a
(1.84)
8. Glu/GABA+ imbalance–0.21a
(0.47)
–0.08a
(0.21)
–0.20a
(0.44)
0.17a
(0.33)
9. Exponent mean
(rest)
–0.08c
(0.17)
0.10c
(0.20)
–0.06c
(0.14)
0.02d
(0.14)
10. Offset mean
(rest)
0.06c
(0.14)
0.14c
(0.35)
0.03c
(0.12)
0.16d
(0.38)
11. Exponent left STS
(rest)
–0.08c
(0.16)
0.06c
(0.14)
–0.04c
(0.12)
–0.04d
(0.14)
12. Offset left STS
(rest)
0.12c
(0.25)
0.18*c
(0.77)
0.08c
(0.17)
0.08d
(0.18)
13. Beta power left STS
(rest)
0.04c
(0.13)
0.23*c
(2.73)
–0.04c
(0.12)
0.05d
(0.15)
EEG language task
9. Exponent mean
(task)
–0.07c
(0.16)
0.13c
(0.30)
–0.10c
(0.21)
0.01d
(0.14)
10. Offset mean
(task)
0.05c
(0.13)
0.14c
(0.34)
0.01c
(0.12)
0.18d
(0.50)
11. Exponent left STS
(task)
–0.03c
(0.12)
0.09c
(0.18)
–0.07c
(0.15)
–0.04d
(0.14)
12. Offset left STS
(task)
0.13c
(0.28)
0.18c
(0.71)
0.07c
(0.15)
0.09d
(0.19)
13. Beta power left STS
(task)
0.07c
(0.15)
0.23*c
(2.84)
0.02c
(0.12)
0.15d
(0.33)
  1. Note: r – Pearson’s correlation coefficient; BF10 – Bayes factor indicating ratio of the likelihood of an alternative hypothesis (H1) to a null hypothesis (H0); mean – values averaged across all electrodes; left STS – values averaged across three electrodes corresponding to the left superior temporal sulcus (T7, TP7, TP9).

  2. ***p < 0.001 (uncorrected); **p < 0.01 (uncorrected); *p < 0.05 (uncorrected)

  3. an = 47 bn = 50; cn = 119; dn = 87

Appendix 1—table 1
Descriptive statistics for electroencephalography (EEG) and magnetic resonance spectroscopy (MRS) results separately for the groups.
DYSCON
MSDMSD
EEG resting statea
Exponent mean (rest)1.540.141.540.18
Exponent left IFG (rest)1.540.161.530.18
Exponent left STS (rest)1.500.181.470.22
Exponent right IFG (rest)1.540.151.540.18
Exponent right STS (rest)1.480.181.450.22
Offset mean (rest)–10.800.19–10.800.24
Offset left IFG (rest)–10.720.34–10.740.33
Offset left STS (rest)–10.970.38–10.980.37
Offset right IFG (rest)–10.790.29–10.810.32
Offset right STS (rest)–10.990.31–11.040.36
Beta power left IFG (rest)0.480.180.480.20
Beta power left STS (rest)0.490.190.480.21
Beta power right IFG (rest)0.490.180.500.19
Beta power right STS (rest)0.510.200.500.21
EEG language taska
Exponent mean (task)1.550.151.560.18
Exponent left IFG (task)1.550.161.550.19
Exponent left STS (task)1.500.201.470.21
Exponent right IFG (task)1.540.171.550.19
Exponent right STS (task)1.470.191.450.22
Offset mean (task)–10.670.25–10.670.28
Offset left IFG (task)–10.580.39–10.600.37
Offset left STS (task)–10.860.44–10.870.42
Offset right IFG (task)–10.650.36–10.680.37
Offset right STS (task)–10.880.36–10.940.41
Beta power left IFG (task)b0.500.230.510.21
Beta power left STS (task)0.540.240.530.23
Beta power right IFG (task)0.510.230.520.21
Beta power right STS (task)0.550.260.550.23
MRS
Gluc1.110.121.070.10
GABA+d0.460.140.440.12
Glu/GABA +ratiod2.670.872.680.75
Glu/GABA +imbalanced0.270.110.240.12
  1. Note: DYS – dyslexic group; CON – control group; mean – values averaged across all electrodes; left IFG – values averaged across three electrodes corresponding to the left inferior frontal gyrus (F7, FT7, FC5); left STS – values averaged across three electrodes corresponding to the left superior temporal sulcus (T7, TP7, TP9); right IFG – values averaged across three electrodes corresponding to the right inferior frontal gyrus (F8, FT8, FC6); right STS – values averaged across three electrodes corresponding to the right superior temporal sulcus (T8, TP8, TP10).

  2. an = 119 (DYS n = 59, CON n = 60); bn = 117 (DYS n = 57, CON n = 60); cn = 50 (DYS n = 21, CON n = 29); dn = 47 (DYS n = 20, CON n = 27).

Appendix 1—table 2
Demographic and behavioral characteristics of the subsample of 47 participants.
DYS(n=20)CON(n=27)TestpEffect sizeBF10
MSDMSD
Demographics
Sex12 F, 8 M11 F, 16 Mχ=1.710.192phi = –0.190.81
Age19.983.9220.333.25U=261.00.846rrb = 0.030.29
Mother’s education
(years)
17.483.6816.632.56U=242.50.551rrb = –0.100.36
Father’s education
(years)
16.74a3.46a16.413.49U=237.00.659rrb = –0.080.31
IQ103.8014.30113.569.23t(30.41)=2.670.012d=0.816.65
Nonverbal IQ
(scaled score)
10.453.1412.372.17U=174.50.038rrb = 0.351.79
ARHQ-PL52.9010.7824.266.71U=6.0<0.001rrb = –0.981060.02
Reading and reading-related tasks
Words/min110.8019.59136.1112.80U=74.0<0.001*rrb = 0.73108.53
Pseudowords/min59.1513.9185.5915.83t(45)=5.96<0.001*d=1.76>10,000
RAN objects (s)32.554.3628.564.46U=122.00.001*rrb = –0.557.92
RAN colors (s)36.404.9130.633.67t(45)=–4.61<0.001*d=–1.36571.37
RAN digits(s)19.354.0016.332.24U=134.00.003*rrb = –0.503.82
RAN letters (s)21.753.2419.522.53t(45)=–2.650.011d=–0.784.52
Reading comprehension (s)63.1018.3043.897.87U=82.0<0.001*rrb = –0.7065.88
Phoneme deletion
(% correct)
74.6229.2094.166.61U=120.00.001*rrb = 0.5613.46
Spoonerisms phonemes
(% correct)
52.1438.2387.838.36U=106.5<0.001*rrb = 0.6117.43
Spoonerisms syllables
(% correct)
41.6729.8677.7820.67U=87.5<0.001*rrb = 0.6827.77
Orthographic awareness
(accuracy/time)
0.350.160.540.13t(45)=4.26<0.001*d=1.26215.19
Perception speed (sten score)2.951.794.331.64U=126.00.002*rrb = 0.534.49
Digits forward5.651.797.042.03U=158.00.014rrb = 0.422.13
Digits backward4.951.737.482.06t(45)=4.45<0.001*d=1.31355.64
  1. Note: DYS – dyslexic group; CON – control group; F – females, M – males. BF10 – Bayes factor indicating ratio of the likelihood of an alternative hypothesis (H1) to a null hypothesis (H0). ARHQ-PL – Polish version of the Adult Reading History Questionnaire. RAN – rapid automatized naming. Non-parametric Mann-Whitney test was performed when assumption of normal distribution was violated. rrb – rank biserial correlation provided as an effect size parameter for Mann-Whitney test. Boldface indicates statistical significance at p<0.05 level (uncorrected).

  2. *Significance after Bonferroni correction for 14 planned comparisons for reading and reading-related tasks;

  3. an = 19 (one participant did not provide information about the father’s education).

Appendix 1—table 3
Zero-order correlations between magnetic resonance spectroscopy (MRS) and electroencephalography (EEG) biomarkers of excitatory-inhibitory balance.
Variable1.
r(BF10)
2.3.4.5.6.7.8.
EEG resting state
1. Glu
2. GABA+0.40**a
(8.08)
3. Glu/GABA +ratio–0.11
(0.24)a
–0.90***a
(>10,000)
4. Glu/GABA +imbalance0.16
(0.32)a
0.33*a
(2.14)
–0.19
(0.42)a
5. Exponent mean (rest)0.13
(0.25)b
0.10
(0.23)a
–0.11
(0.24)a
0.31*
(1.58)a
6. Offset mean (rest)0.10
(0.23)b
0.17
(0.35)a
–0.15
(0.30)a
0.26
(0.79)a
0.70***
(>10,000)c
7. Exponent
left STS (rest)
0.04
(0.18)b
0.09
(0.22)a
–0.08
(0.21)a
0.37*
(4.06)a
0.70***
(>10,000)c
0.52***
(>10,000)c
8. Offset
left STS (rest)
–0.09
(0.21)b
0.01
(0.18)a
0.00
(0.18)a
0.21
(0.48)a
0.24**
(3.92)c
0.54***
(>10,000)c
0.70***
(>10,000)c
9. Beta power
left STS (rest)
–0.06
(0.19)b
0.22
(0.53)a
–0.28
(1.04)a
0.03
(0.19)a
0.19*
(0.85)c
0.20*
(1.29)c
0.50***
(>10,000)c
0.56***
(>10,000)c
EEG language task
5. Exponent mean (task)0.14
(0.27)b
0.15
(0.29)a
–0.14
(0.28)a
0.24
(0.64)a
6. Offset mean (task)0.10
(0.23)b
0.19
(0.41)a
–0.16
(0.31)a
0.23
(0.59)a
0.75***
(>10,000)c
7. Exponent
left STS (task)
0.05
(0.19)b
0.09
(0.22)a
–0.09
(0.21)a
0.24
(0.64)a
0.69***
(>10000)c
0.58***
(>10,000)c
8. Offset
left STS (task)
–0.07
(0.20)b
0.03
(0.19)a
–0.01
(0.18)a
0.18
(0.38)a
0.36***
(375.45)c
0.62***
(>10,000)c
0.80***
(>10,000)c
9. Beta power
left STS (task)
–0.07
(0.20)b
0.22
(0.52)a
–0.25
(0.77)a
–0.01
(0.18)a
0.07
(0.15)c
0.15
(0.43)c
0.47***
(>10,000)c
0.59***
(>10,000)c
  1. Note: r – Pearson’s correlation coefficient; BF10 – Bayes factor indicating ratio of the likelihood of an alternative hypothesis (H1) to a null hypothesis (H0); mean – values averaged across all electrodes; left STS – values averaged across three electrodes corresponding to the left superior temporal sulcus (T7, TP7, TP9).

  2. ***p < 0.001 (uncorrected); **p < 0.01 (uncorrected); *p < 0.05 (uncorrected).

  3. an = 47; bn = 50; cn = 119.

Appendix 1—table 4
Zero-order correlations between reading, phonological awareness, rapid automatized naming, multisensory integration, and biomarkers of excitatory-inhibitory balance.
Variable1.
r(BF10)
2.3.4.
EEG resting state
1. Reading
2. Phonological awareness0.62***
(>10,000)c
3. RAN0.73***
(>10,000)c
0.50***
(>10,000)c
4. Multisensory integration0.24*
(1.44)d
0.33**
(16.94)d
0.08
(0.17)d
5. Glu–0.10
(0.22)b
–0.16
(0.32)b
–0.08
(0.20)b
–0.08
(0.21)b
6. GABA+–0.17
(0.35)a
–0.02
(0.18)a
–0.19
(0.41)a
0.24
(0.62)a
7. Glu/GABA+ ratio0.03
(0.19)a
–0.02
(0.18)a
0.08
(0.21)a
–0.31*
(1.62)a
8. Glu/GABA+ imbalance–0.21
(0.47)a
–0.06
(0.20)a
–0.19
(0.40)a
0.13
(0.27)a
9. Exponent mean
(rest)
–0.13
(0.30)c
0.05
(0.13)c
–0.08
(0.17)c
–0.05
(0.15)d
10. Offset mean
(rest)
–0.03
(0.12)c
0.04
(0.13)c
–0.02
(0.12)c
–0.01
(0.14)d
11. Exponent left STS
(rest)
–0.14
(0.37)c
–0.01
(0.12)c
–0.07
(0.16)c
–0.15
(0.33)d
12. Offset left STS
(rest)
0.03
(0.12)c
0.09
(0.18)c
0.03
(0.12)c
–0.07
(0.17)d
13. Beta power left STS`
(rest)
0.03
(0.12)c
0.22*
(1.96)c
–0.04
(0.12)c
0.04
(0.14)d
EEG language task
9. Exponent mean
(task)
–0.13
(0.32)c
0.06
(0.14)c
–0.14
(0.34)c
–0.09
(0.19)d
10. Offset mean
(task)
–0.05
(0.13)c
0.04
(0.12)c
–0.05
(0.13)c
–0.01
(0.13)d
11. Exponent left STS
(task)
–0.11
(0.23)c
0.01
(0.12)c
–0.11
(0.24)c
–0.17
(0.48)d
12. Offset left STS
(task)
0.04
(0.12)c
0.09
(0.18)c
0.01
(0.12)c
–0.07
(0.16)d
13. Beta power left STS
(task)
0.05
(0.13)c
0.21*
(1.61)c
0.02
(0.12)c
0.11
(0.22)d
  1. Note: r – Pearson’s correlation coefficient; BF10 – Bayes factor indicating ratio of the likelihood of an alternative hypothesis (H1) to a null hypothesis (H0); mean – values averaged across all electrodes; left STS – values averaged across three electrodes corresponding to the left superior temporal sulcus (T7, TP7, TP9).

  2. ***p < 0.001 (uncorrected); **p < 0.01 (uncorrected); *p < 0.05 (uncorrected).

  3. an = 47; bn = 50; cn = 119; dn = 87.

Appendix 1—table 5
One-sample t-tests separately for visual (words>false fonts) and auditory runs (words>words backward) within control (CON) and dyslexic groups (DYS).
Brain regionsHemispherePeak of cluster coordinatest-ValueNumber of voxels
xyz
CON (n=29) visual runs (FWEc = 127)
Middle temporal gyrus, superior temporal gyrus, inferior frontal gyrus (pars triangularis, orbitalis, opercularis), supramarginal gyrus, temporal pole (superior and middle temporal gyri), postcentral gyrus, hippocampus, posterior orbital gyrus, amygdala, parahippocampal gyrus, inferior parietal gyrus, rolandic operculum, lateral orbital gyrus, pallidum, anterior orbital gyrus, inferior frontal gyrus, putamen, middle frontal gyrus, insula, angular gyrusL–62–54109.785360
Precuneus, posterior cingulate gyrus, middle cingulate and paracingulate gyri, calcarine fissure, cuneusL/R-4–52207.17957
Middle temporal gyrus, superior temporal gyrus, temporal pole (superior and middle temporal gyri), inferior temporal gyrusR60–3407.101117
Cuneus, superior occipital gyrus, calcarine fissureL/R14–92265.98334
Superior frontal gyrus (dorsolateral and medial)L-854345.95678
Superior frontal gyrus (medial orbital), anterior cingulate cortex (pregenual), superior frontal gyrus (medial)L/R-850-85.33194
Insula, putamen, rolandic operculum, inferior frontal gyrus (pars opercularis)L–36465.09127
Angular gyrus, middle temporal gyrus, middle occipital gyrusL–40–58264.56167
DYS (n=21) visual runs (FWEc = 665)
Inferior frontal gyrus (pars triangularis, orbitalis, opercularis), posterior orbital gyrusL–482406.54665
Middle temporal gyrus, superior temporal gyrus, supramarginal gyrusL–62–3226.40677
CON (n=29) auditory run (FWEc = 124)
Middle temporal gyrus, postcentral gyrus, cingulate gyrus (mid part), superior parietal gyrus, precentral gyrus, precuneus, middle occipital gyrus, supramarginal gyrus, inferior parietal gyrus, inferior temporal gyrus, calcarine fissure, fusiform gyrus, superior temporal gyrus, angular gyrus, temporal pole (superior and temporal gyri), cuneus, supplementary motor area, lingual gyrus, inferior frontal gyrus (pars triangularis, orbitalis), posterior cingulate gyrus, superior frontal gyrus (dorsolateral), inferior occipital gyrus, posterior orbital gyrus, paracentral lobule, rolandic operculum, parahippocampal gyrus, lateral orbital gyrus, middle frontal gyrus, superior occipital gyrus, anterior orbital gyrusL/R–52–12–109.7513,444
Middle temporal gyrus, superior temporal gyrus, temporal pole (superior temporal gyrus), inferior temporal gyrusR50–4–207.981038
Putamen, insula, rolandic operculum, pallidum, Heschl’s gyrus, amygdala, hippocampus, thalamus (lateral geniculate)L–34–2086.15463
Thalamus (mediodorsal medial magnocellular, intralaminar, pulvinar medial, ventral posterolateral, mediodorsal lateral parvocellular, ventral lateral, pulvinar anterior, lateral posterior, pulvinar inferior, medial geniculate), hippocampus, lingual gyrus, parahippocampal gyrusL–10–2005.82259
Superior frontal gyrus (dorsolateral and medial), middle frontal gyrusL/R–6,60345.81438
Thalamus (ventral posterolateral, pulvinar medial, pulvinar anterior, lateral geniculate, intralaminar, pulvinar inferior, pulvinar lateral, ventral lateral, medial geniculate), hippocampusR14–2045.68149
Middle occipital gyrus, middle temporal gyrus, angular gyrus, superior occipital gyrus, superior temporal gyrus, superior parietal gyrus, cuneus, supramarginal gyrus, inferior parietal gyrusR36–76405.46978
Cerebellar hemispheres (lobules IV, V, VI, VIII, IX, crus I, crus II), vermis (lobules VI, VII, VIII, IX)L/R12,–50–305.35630
Inferior frontal gyrus (pars triangularis)L–5022204.54124
DYS (n=21) auditory run (FWEc = 192)
Middle temporal gyrus, superior temporal gyrus, temporal pole (superior temporal gyrus)L–52–12–127.00403
Middle temporal gyrusL–56–3805.33192
  1. Note: L – left, R – right. All results are reported at p<0.001 height threshold corrected for multiple comparisons using p<0.05 FWE cluster threshold.

Appendix 1—table 6
Logical conjunction results from paired t-tests for both visual (words>false fonts) and auditory runs (words>words backward) within control (CON) and dyslexic groups (DYS).
Brain regionsHemispherePeak of cluster coordinatest-ValueNumber of voxels
xyz
CON (n=29) visual and auditory runs conjunction (FWEc = 143)
 Middle temporal gyrus, inferior frontal gyrus (pars triangularis, orbitalis, opercularis), superior temporal gyrus, temporal pole (superior and middle temporal gyri), supramarginal gyrus, posterior orbital gyrus, inferior parietal gyrus, angular gyrus, lateral orbital gyrus, anterior orbital gyrus, middle occipital gyrus L–54–58129.583164
 Middle temporal gyrus, superior temporal gyrus R50–34–27.98186
 Superior frontal gyrus (dorsolateral and medial) L-656365.94430
 Superior temporal gyrus, temporal pole (superior and middle temporal gyri), middle temporal gyrus R5210–165.76175
 Precuneus, posterior cingulate gyrus, middle cingulate and paracingulate gyri L/R–2–56305.48448
 Supramarginal gyrus, postcentral gyrus, rolandic operculum, superior temporal gyrus L–50–24225.35217
 Middle cingulate and paracingulate gyri L/R–2–10405.17143
DYS (n=21) visual and auditory runs conjunction (FWEc = 545)
Middle temporal gyrus, superior temporal gyrus L–52–12–66.73545
  1. Note: L – left, R – right. All results are reported at p<0.001 height threshold corrected for multiple comparisons using p<0.05 FWE cluster threshold.

Appendix 1—table 7
Results from the flexible factorial model for both visual (words>false fonts) and auditory runs (words>words backward) between control (CON) and dyslexic groups (DYS).
Brain regionsHemispherePeak of cluster coordinatest-ValueNumber of voxels
xyz
CON>DYS main effect of group for visual and auditory runs (FWEc = 121)
 Middle temporal gyrus, supramarginal gyrus, superior temporal gyrus, inferior temporal gyrus, inferior parietal gyrus, middle occipital gyrus, fusiform gyrus, angular gyrus, postcentral gyrus, rolandic operculum L–58–30307.192181
 Rolandic operculum, superior temporal gyrus, middle temporal gyrus, supramarginal gyrus, insula, postcentral gyrus, precentral gyrus, inferior frontal gyrus (pars opercularis), putamen, inferior temporal gyrus R46–4147.071512
 Superior frontal gyrus (dorsolateral and medial), middle frontal gyrus R2234546.56464
 Insula, rolandic operculum, superior temporal gyrus, Heschl’s gyrus L–40–6186.21266
 Middle cingulate and paracingulate gyri, paracentral lobule, supplementary motor area, precuneus L/R–8–26425.961108
 Precentral gyrus, superior frontal gyrus (dorsolateral) L–32–8685.70191
 Angular gyrus, supramarginal gyrus, inferior parietal gyrus, middle occipital gyrus R52–62405.37516
 Superior frontal gyrus (dorsolateral and medial), middle frontal gyrus L–1838345.13241
 Precentral gyrus, postcentral gyrus, middle frontal gyrus R58–12465.09140
 Lingual gyrus, cerebellar hemisphere (lobule VI), fusiform gyrus R14–68–105.04209
 Temporal pole (superior temporal gyrus), middle temporal gyrus, superior temporal gyrus L–5410–165.00156
 Supplementary motor area, middle cingulate and paracingulate gyri L/R46484.93174
 Postcentral gyrus, superior parietal gyrus R24–48584.85121
 Postcentral gyrus, supramarginal gyrus R38–32364.78137
DYS>CON main effect of group for visual and auditory runs
 No suprathreshold clusters
  1. Note: L – left, R – right. All results are reported at p<0.001 height threshold corrected for multiple comparisons using p<0.05 FWE cluster threshold.

Appendix 1—table 8
Magnetic resonance spectroscopy (MRS) checklist.
Hardware
Field strength [T]7T
ManufacturerGE Healthcare
Model (software version if available)Discovery MR 950
RF coils: nuclei (transmit/ receive), number of channels, type, body partH1 32 channel head coil
Additional hardwareMR Safe Response PAD
Acquisition
Pulse sequenceSemi-Laser
Volume of interest (VOI) locationsLeft superior temporal sulcus (STS)
Nominal VOI size [cm3, mm3]15×15×15 mm3
Repetition time (TR), echo time (TE) [ms, s]TR = 4000 ms TE = 28 ms
Total number of excitations or acquisitions per spectrum320 averages
Additional sequence parameters
(spectral width in Hz, number of spectral points, frequency offsets)
5000 Hz, 4096 points
Water suppression methodVAPOR
Shimming method, reference peak, and thresholds for ‘acceptance of shim’ chosenAutomated linear shims adjustment, 0 and first shim order only, water peak, <20 Hz
Triggering or motion correction method
(respiratory, peripheral, cardiac triggering, incl. device used and delays)
N/A
Data analysis methods and outputs
Analysis softwarefsl-mrs (version 2.0.7)
Processing steps deviating from quoted reference or productfsl_mrs default pipeline+simulated basis set
Output measure
(e.g. absolute concentration, institutional units, ratio)
Processing steps deviating from quoted reference or product
Ratio to total creatine
Quantification references and assumptions, fitting model assumptionsThe customized basis set includes: Ala, Asc, Asp, Cit, Cr, EtOH, GABA, GPC, GSH, Glc, Gln, Glu, Gly, Ins, Lac, NAA, NAAG, PCh, PCr, PE, Phenyl, Scyllo, Ser, Tau, Tyros, bHB, bHG. Macromolecules MM09, MM12, MM14, MM17, MM21 were added.
Data quality
Reported variables
(SNR, linewidth
(with reference peaks))
Linewidth (for metabolite group reported by fsl_mrs): mean 11.11 Hz, SD = 2.79 Hz SNR (for NAA reported by fsl_mrs): mean 7.05, SD = 23.95
Data exclusion criteriaLinewidth >20 Hz, CRLB >20% and Visual inspection (baseline, residuals)
Glu – 4 subjects excluded GABA – 7 subjects excluded
Quality measures of postprocessing
Model fitting (e.g. CRLB, goodness of fit, SD of residual)
%CRLB of Glu: mean 2.96, SD = 0.79
%CRLB of GABA: mean 10.59, SD = 2.76
%CRLB of NAA: 1.76 SD = 0.46
Sample spectrum

Additional files

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. Agnieszka Glica
  2. Katarzyna Wasilewska
  3. Julia Jurkowska
  4. Jarosław Żygierewicz
  5. Bartosz Kossowski
  6. Katarzyna Jednoróg
(2025)
Reevaluating the neural noise in dyslexia using biomarkers from electroencephalography and high-resolution magnetic resonance spectroscopy
eLife 13:RP99920.
https://doi.org/10.7554/eLife.99920.4