Materials and methods
78
Participants 79
Thirty-two volunteers participated in this study (17 female, 19-28 years old; M = 22.2, SD 80
= 2.50). Participants declared never having received a diagnosis of neurological/psychological 81
disorders or dyslexia. All participants were right-handed (Oldfield, 1971). Self-reports of hearing 82
ability (Five-minute Hearing test, revised version) (Koike et al., 1994) indicated that 31 83
participants had no hearing impairments, while one participant was recommended a hearing test 84
(score of 23/60, with 20 being the cutoff) . This study was approved by the School of Social 85
Sciences Research Ethics Committee at the University of Dundee (approval number: UoD-SoSS-86
PSY-UG-2021-263) and adhered to the guidelines for the treatment of human participants in the 87
Declaration of Helsinki . Volunteers received monetary compensation of £10/h. Research 88
questions, hypotheses, measured variables and analyses were pre -registered on the OSF 89
website (https://osf.io/xrq36). Deviations from the pre-registration are detailed where necessary. 90
Procedure 91
Online procedure 92
Prior to taking part in the EEG experiment, p articipants were asked to complete the first 93
part of the study online, using the experiment builder Gorilla (Anwyl-Irvine et al., 2020) . In this 94
online session, participants completed a demographics questionnaire, a musicality assessment 95
(Müllensiefen et al., 2013), a handedness questionnaire (Oldfield, 1971) and the lexical decision 96
task. The musicality questionnaire included 14 items from the Goldsmiths Musical 97
Sophistication Index (GMSI) covering perceptual abilities, active engagement, singing abilit ies 98
and musical training. Each item included a statement such as “I spend a lot of my free time doing 99
music-related activities ”, which participants rated on a scale of 1 to 7 . The online task and 100
questionnaires had to be completed on a PC and in a single, uninterrupted session. 101
In-person procedure 102
After completing the online portion of the study, participants were invited to take part in 103
the in-person EEG session. Participants performed the EEG experiment in a soundproof booth 104
(160 cm x 110 cm) where they were seated approximately 65cm away from a 24 -inch monitor. 105
Participants were equipped with high-quality, wired headphones (Sennheiser, HD 25, 75Ω). The 106
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
experiment was run using Psychtoolbox (version ‘3.0.17)(Brainard & Vision, 1997; Pelli & Vision, 107
1997) through MATLAB (MATLAB, 2021). During the presentation of speech and music stimuli, 108
participants were instructed to relax, limit motor movement to avoid movement artefacts and 109
fixate on a small green circle on the screen. Blocks were self-paced, and the order of the auditory 110
stimuli was counterbalanced. 111
Stimuli and tasks 112
Lexical Decision Task 113
The lexical decision task was adapted from a pre -existing task (no reference was found 114
for the task or the study in which it was used) on the experiment-building website Gorilla (Anwyl-115
Irvine et al., 2020). The words and non-words used in this task were taken from a previous study 116
(Yao et al., 2018) and are part of the corpus of Glasgow Norms (Scott et al., 2019) . Both words 117
and nonwords ranged from 3 to 11 letters in length (M = 6.11, SD = 1.75). All words were concrete 118
nouns with neutral valence, as words with differing emotional connotations can affect response 119
times (Yao et al., 2018) . Word frequency (occurrences per million as per the British National 120
Corpus) was on average M = 23.95 (SD = 36.49). In the task, participants were asked to use the 121
left and right arrow keys to indicate whether the string of letters on the screen was a word (e.g., 122
‘statue’) or a non-word (e.g., ‘depane’, which does not have any meaning in the English language). 123
Participants were given 12 practice trials with feedback. At the beginning of each block, there 124
were two brief screens, appearing for 1000 ms each, saying “Ready?” and “GO!” to prepare the 125
participant for the trials. A fixation cross was also shown for 500 milliseconds before each word 126
and nonword in the task. The experimental task consisted of 90 trials, given in 2 blocks of 45 trials 127
each, and feedback was not provided. All tasks and questionnaires used Open Sans font. 128
Halfway through the lexical decision task, participants were given the option to take a break and 129
recommence when they preferred. 130
Passive Listening to speech and music 131
All auditory stimuli were presented at a sampling rate of 44,100 Hz. The short story 132
selected for this study was “The Elves and the Shoemaker”, originally written by the brothers 133
Grimm, read by a female speaker with a pleasant voice (https://librivox.org/). The length of the 134
story was 300 s ( 5 minutes ). The articulation rate of the stimulus was calculated using Praat 135
(Boersma, 2001), through the automatic detection of syllable nuclei based on intensity peaks in 136
the speech signal (De Jong & Wempe, 2009) . The stimulus had an articulation rate of 3.61 137
syllables per second. Similarly, the modulation spectrum of the speech stimulus showed a peak 138
at 3.75 Hz (Figure 1A). 139
The music piece selected was “Fluid” by Lin Rountree, a pleasant jazz piece at 95 BPM. 140
This piece featured bass, keyboard and trumpet; there were no vocals included. The duration was 141
255 s (4 minutes and 25 seconds), and the modulation spectrum showed a peak at 3.15Hz (Figure 142
1A). After each piece, participants rated how familiar and how pleasant they had found the 143
stimuli by using a Visual Analog Scale: this involved placing a vertical marker between “did not 144
enjoy at all” and “enjoyed a lot” for enjoyment and “not at all familiar” and “extremely familiar” 145
for familiarity (each analysed in arbitrary units between 0 and 100). Participants rated the story 146
(M = 45.67, SD = 35.17) as significantly more familiar than the music piece (M = 27.27, SD = 25.17) 147
(t(31) = 2.55, p = .015, Figure 1C). However, enjoyment ratings did not differ significantly between 148
the story (M = 58.72, SD = 24.38) and the music piece ( M = 64.80, SD = 18.98) (t(31) = -1.08, p = 149
.287, Figure 1B). All stimuli are available on the OSF server (https://osf.io/xrq36). 150
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
151
Figure 1 . (A) Modulation spectrum of speech (green) and music (purple) stimuli. Thick lines indicate 152
average values across 6 -s segments, and shaded areas represent the standard error of the mean. Peak 153
frequency is shown with dotted lines. (B) Enjoyment ratings for speech and music excerpts . There was no 154
significant difference between stimulus ratings (p = .287 ). (C) Familiarity ratings for speech and music 155
excerpts. The story was rated as more familiar than the music piece (t(31) = 2.55, p = .015). Note: Dots in 156
(B) and (C) show individual data points, violin plots show kernel density estimates , and boxplots show 157
median interquartile ranges and minimum/maximum. 158
Analysis of behavioural data 159
For analysis of the lexical decision task , incorrect trials (M = 5.69, SD = 3.25) were 160
excluded from analysis. For correct trials, responses faster than 250 ms or exceeding 1500 ms 161
were excluded (M = 6.58, SD = 10.45), as these could reflect inadvertent movements or lapses in 162
attention (Yao et al., 2018). Three participants were excluded from further analyses involving the 163
lexical decision task as their total number of rejected trials exceeded 2 SDs from the group mean. 164
For the remaining participants, median reaction times and accuracy rates were calculated 165
separately for word and nonword trials. These values were then averaged to yield total accuracy 166
and reaction time measures. 167
Acoustic envelope pre-processing 168
To analyse the neural tracking of speech and music signals, the wideband envelope of 169
the stimuli was extracted. The acoustic waveforms were filtered into eight frequency bands 170
(between 100 and 8000 Hz, 3rd order Butterworth filter, forward and reverse) equally distant on 171
the cochlear frequency map (Smith et al., 2002). The signal in each of these frequency bands was 172
then Hilbert -transformed and the magnitude extracted before being averaged to obtain the 173
wideband music and speech envelopes used in further analyses. Lastly, envelopes were down-174
sampled to a sampling rate of 150 Hz (Keitel et al., 2018). 175
EEG acquisition and pre-processing 176
EEG was recorded from 64 scalp electrodes and digitally sampled at 512 Hz , using a 177
BioSemi ActiveTwo system. Scalp electrodes were positioned according to the international 10-178
20 system. Electrodes with an offset of greater/less than ±20 mV were adjusted prior to starting 179
the recording. Ultimately, electrode offset for all electrodes was below an absolute value of 30 180
mV before the experiment began. Lateral eye movements were monitored by two electro -181
oculographic electrodes placed at the outer canthus of each eye. Vertical eye movements and 182
blinks were monitored by two electro -oculographic electrodes positioned below and above the 183
left eye. 184
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
Data pre-processing was conducted using the Fieldtrip Toolbox (Oostenveld et al., 2011) 185
and custom -made scripts in MATLAB (R2025a) (The Mathworks, 2025) . The data were cut 186
according to the length of the stimul i, with an additional 2 -s leading and trailing window . Data 187
was first re-referenced to Cz, then, fourth-order Butterworth low-pass (60 Hz) and high-pass (0.2 188
Hz) filters were applied. Noisy (e.g., with higher variance than most) channels were visually 189
identified and interpolated through triangulation. A maximum of five channels was interpolated 190
per participant (M = 1.25, SD = 1.34). To remove eye artifacts and blinks, an independent 191
component analysis (ICA) was then carried out for 30 principal components . Components, 192
including eye movements or blinks, were selected and removed (M = 1.40, SD = 0.49). The EEG 193
signal was then down-sampled to 150 Hz to match the envelope signals (Keitel et al., 2018). 194
MI analysis 195
The correspondence between the continuous EEG signal and t he acoustic envelope 196
signals was analysed using a Gaussian copula mutual information framework (Ince et al., 2017). 197
The Mutual Information (MI) between the continuous L1-normalised EEG signal s and the 198
envelopes of the auditory stimuli (as well as their derivat ives, see Fig. 2 ) was calculated in the 199
frequency domain by applying a continuous wavelet transform (Chalas et al., 2022) , for 63 200
logarithmically spaced frequencies between .25 and 20 Hz. We used a participant-specific 201
optimal brain-stimulus lag computed by identifying, for each participant and condition, the lag at 202
which the MI value was highest (identified at electrode Cz for slow frequencies, averaged 203
between 0.5 – 4 Hz) . Each MI value was computed per participant, condition, frequency and 204
channel. To normalise MI, we first created random (surrogate) MI distributions. For this, we 205
segmented the continuous envelope signals into 5 -s segments and shuffled the segments 206
randomly. This kept the statistical properties of the signal but destroyed the temporal 207
relationship between the acoustic stimuli and brain signals. MI was then computed between the 208
brain signal and the shuffled envelope signals. Normalised MI was computed by z -scoring the 209
observed MI values against a surrogate distribution, subtracting the mean and dividing by the 210
standard deviation of randomised MI estimates. The normalised MI values were used for all 211
further analyses. 212
213
Figure 2. Brief excerpt (10 s) of acoustic envelope , waveform, and first derivative for the speech (green) 214
and music stimuli (purple). 215
216
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
Statistical analysis 217
For each participant, normalised MI values at each channel and frequency were tested 218
against a null baseline of zero using dependent -samples t -tests. Cluster -based multiple -219
comparison correction (implemented in FieldTrip (Oostenveld et al., 2011) ) was applied using 220
5000 Monte Carlo permutations, in which the condition labels ( normalised MI vs. zero) were 221
randomly exchanged within subjects to form the null distribution. Clusters were defined as 222
spanning more than 3 frequency bins , and the cluster -level statistic was the sum of t-values. 223
Observed clusters were considered significant if their cluster statistic exceeded the 95th 224
percentile of the permutation distribution . This controlled for the family -wise error rate at the 225
cluster level. To avoid reporting spurious results, we only report clusters that cover a range of 226
frequencies of more than 0.5 Hz. 227
For the comparison between the two tracking conditions, dependent samples t-tests 228
were computed across channels and frequencies with MI values of one condition compared 229
against MI values of the other , with the same cluster -permutation approach used for multiple 230
comparison as described above. 231
To test the relationship between cortical tracking and behavioural measures (reaction 232
times from the lexical decision task) , Pearson’s correlations were computed between the 233
behavioural measures and the normalised MI values , across all electrodes and frequencies . 234
Before comparing the r values with the permutation distribution using cluster-based permutation 235
(using a minimum cluster size of 3 channels and tested against the 95 th percentile of the 236
permutation distribution), Pearson’s r values were transformed to be normally distributed using 237
Fisher’s z-transformation (Gorsuch & Lehmann, 2010). As an indicator of effect sizes, we report 238
Cohen’s d for peak electrodes. 239
To further compare the contribution of cortical tracking in both speech and music 240
conditions on performance in the lexical decision task , a robust multiple linear regression was 241
computed (R 4.5.1). The model included: (i) MI values, (ii) stimulus type (speech or music), (iii) 242
frequency band (delta or alpha), (iv) musical sophistication scores, (v) age , and (vi) reading 243
enjoyment as predictors , along with the three -way interaction between MI × stimulus type × 244
musical sophistication and the two -way interaction between MI × frequency band. Reaction 245
times in the lexical decision task served as the outcome variable . All continuous variables were 246
z-scored. 247
To explicitly test for hemispheric lateralisation in our significant correlation clusters, we 248
computed a lateralisation index (LI) following the procedure outlined by Haegens et al. (2011) 249
which is comparable to the method used by Thut et al. (2006). The LI was calculated as: 250
LI = (ipsilateral ROI – contralateral ROI) / (ipsilateral ROI + contralateral ROI) 251
where the region of interest (ROI) comprised the electrodes within each significant lateral 252
cluster. For each participant, correlation values ( r) from electrodes in one hemisphere were 253
paired with their hemispheric counterparts. The resulting LI values were then tested against zero 254
using a two -sided one -sample t-test to determine whether correlations differed significantly 255
between hemispheres. 256
257
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
References
555
Abrams, D. A., Nicol, T., Zecker, S., & Kraus, N. (2008). Right-hemisphere auditory cortex is 556
dominant for coding syllable patterns in speech. Journal of Neuroscience, 28(15), 3958-557
3965. 558
Ahissar, E., Nagarajan, S., Ahissar, M., Protopapas, A., Mahncke, H., & Merzenich, M. M. (2001). 559
Speech comprehension is correlated with temporal response patterns recorded from 560
auditory cortex. Proceedings of the National Academy of Sciences, 98(23), 13367-561
13372. 562
Anwyl-Irvine, A. L., Massonnié, J., Flitton, A., Kirkham, N., & Evershed, J. K. (2020). Gorilla in our 563
midst: An online behavioral experiment builder. Behavior research methods, 52(1), 388-564
407. 565
Araújo, J., Simons, B. D., Peter, V., Mandke, K., Kalashnikova, M., Macfarlane, A., Gabrielczyk, 566
F., Wilson, A., Di Liberto, G. M., & Burnham, D. (2024). Atypical low-frequency cortical 567
encoding of speech identifies children with developmental dyslexia. Frontiers in Human 568
Neuroscience, 18, 1403677. 569
Asano, R. (2022). The evolution of hierarchical structure building capacity for language and 570
music: a bottom-up perspective. Primates, 63(5), 417-428. 571
Atanasova, T., Gross, J., Rimmele, J. M., & Keitel, A. (2026). The involvement of endogenous 572
brain rhythms in speech processing. Neuroscience & Biobehavioral Reviews, 106568. 573
Baudouin, A., Vanneste, S., & Isingrini, M. (2004). Age-related cognitive slowing: The role of 574
spontaneous tempo and processing speed. Experimental aging research, 30(3), 225-575
239. 576
Bhide, A., Power, A., & Goswami, U. (2013). A rhythmic musical intervention for poor readers: A 577
comparison of efficacy with a letter‐based intervention. Mind, Brain, and Education, 578
7(2), 113-123. 579
Boersma, P. (2001). Praat, a system for doing phonetics by computer. Glot. Int., 5(9), 341-345. 580
Brainard, D. H., & Vision, S. (1997). The psychophysics toolbox. Spatial vision, 10(4), 433-436. 581
Cason, N., Astésano, C., & Schön, D. (2015). Bridging music and speech rhythm: Rhythmic 582
priming and audio–motor training affect speech perception. Acta psychologica, 155, 43-583
50. 584
Chalas, N., Daube, C., Kluger, D. S., Abbasi, O., Nitsch, R., & Gross, J. (2022). Multivariate 585
analysis of speech envelope tracking reveals coupling beyond auditory cortex. 586
NeuroImage, 258, 119395. 587
De Jong, N. H., & Wempe, T. (2009). Praat script to detect syllable nuclei and measure speech 588
rate automatically. Behavior research methods, 41(2), 385-390. 589
der Nederlanden, C. M. V. B., Joanisse, M. F., & Grahn, J. A. (2020). Music as a scaffold for 590
listening to speech: Better neural phase-locking to song than speech. NeuroImage, 214, 591
116767. 592
Di Liberto, G. M., Peter, V., Kalashnikova, M., Goswami, U., Burnham, D., & Lalor, E. C. (2018). 593
Atypical cortical entrainment to speech in the right hemisphere underpins phonemic 594
deficits in dyslexia. NeuroImage, 175, 70-79. 595
Ding, N., Patel, A. D., Chen, L., Butler, H., Luo, C., & Poeppel, D. (2017). Temporal modulations 596
in speech and music. Neuroscience & Biobehavioral Reviews, 81, 181-187. 597
Doelling, K. B., & Poeppel, D. (2015). Cortical entrainment to music and its modulation by 598
expertise. Proceedings of the National Academy of Sciences, 112(45), E6233-E6242. 599
Drakoulaki, K., Anagnostopoulou, C., Guasti, M. T., Tillmann, B., & Varlokosta, S. (2024). 600
Situating language and music research in a domain‐specific versus domain‐general 601
framework: A review of theoretical and empirical data. Language and Linguistics 602
Compass, 18(2), e12509. 603
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
Flaugnacco, E., Lopez, L., Terribili, C., Zoia, S., Buda, S., Tilli, S., Monasta, L., Montico, M., Sila, 604
A., & Ronfani, L. (2014). Rhythm perception and production predict reading abilities in 605
developmental dyslexia. Frontiers in Human Neuroscience, 8, 392. 606
Garcia-de-Soria, M. C., Mathias, B., Keitel, A., & Klimovich-Gray, A. (2025). Benefits of music 607
training for learning to read: Evidence from cortical tracking of speech in children. 608
bioRxiv, 2025.2009. 2005.674218. 609
Giraud, A.-L., & Poeppel, D. (2012). Cortical oscillations and speech processing: emerging 610
computational principles and operations. Nature neuroscience, 15(4), 511-517. 611
Gorsuch, R. L., & Lehmann, C. S. (2010). Correlation coefficients: Mean bias and confidence 612
interval distortions. Journal of Methods and Measurement in the Social Sciences, 1(2), 613
52-65. 614
Goswami, U. (2007). Phonological representations for reading acquisition across languages. In 615
Single-Word Reading (pp. 78-97). Psychology Press. 616
Goswami, U. (2011). A temporal sampling framework for developmental dyslexia. Trends in 617
cognitive sciences, 15(1), 3-10. 618
Goswami, U., Huss, M., Mead, N., Fosker, T., & Verney, J. P. (2013). Perception of patterns of 619
musical beat distribution in phonological developmental dyslexia: Significant 620
longitudinal relations with word reading and reading comprehension. Cortex, 49(5), 621
1363-1376. 622
Gross, J., Hoogenboom, N., Thut, G., Schyns, P., Panzeri, S., Belin, P., & Garrod, S. (2013). 623
Speech rhythms and multiplexed oscillatory sensory coding in the human brain. PLoS 624
biology, 11(12), e1001752. 625
Haegens, S., & Golumbic, E. Z. (2018). Rhythmic facilitation of sensory processing: A critical 626
review. Neuroscience & Biobehavioral Reviews, 86, 150-165. 627
Haegens, S., Händel, B. F., & Jensen, O. (2011). Top-down controlled alpha band activity in 628
somatosensory areas determines behavioral performance in a discrimination task. 629
Journal of neuroscience, 31(14), 5197-5204. 630
Harding, E. E., Sammler, D., Henry, M. J., Large, E. W., & Kotz, S. A. (2019). Cortical tracking of 631
rhythm in music and speech. NeuroImage, 185, 96-101. 632
Ince, R. A., Giordano, B. L., Kayser, C., Rousselet, G. A., Gross, J., & Schyns, P. G. (2017). A 633
statistical framework for neuroimaging data analysis based on mutual information 634
estimated via a gaussian copula. Human brain mapping, 38(3), 1541-1573. 635
Jerrim, J., Lopez-Agudo, L. A., & Marcenaro-Gutierrez, O. D. (2020). Does it matter what children 636
read? New evidence using longitudinal census data from Spain. Oxford Review of 637
Education, 46(5), 515-533. 638
Keitel, A., Gross, J., & Kayser, C. (2018). Perceptually relevant speech tracking in auditory and 639
motor cortex reflects distinct linguistic features. PLoS biology, 16(3), e2004473. 640
Keitel, A., Pelofi, C., Guan, X., Watson, E., Wight, L., Allen, S., Mencke, I., Keitel, C., & Rimmele, 641
J. (2025). Cortical and behavioral tracking of rhythm in music: Effects of pitch 642
predictability, enjoyment, and expertise. Annals of the New York Academy of Sciences, 643
1546(1), 120-135. 644
Keshavarzi, M., Moore, B. C., & Goswami, U. (2025). Atypical low-frequency and high-frequency 645
neural entrainment to rhythmic audiovisual speech in adults with dyslexia. bioRxiv, 646
2025.2009. 2014.675562. 647
Koike, K. J., Hurst, M. K., & Wetmore, S. J. (1994). Correlation between the American Academy 648
of Otolaryngology—Head and Neck Surgery Five-Minute Hearing Test and Standard 649
Audiologic Data. Otolaryngology—Head and Neck Surgery, 111(5), 625-632. 650
Kotz, S. A., Ravignani, A., & Fitch, W. T. (2018). The evolution of rhythm processing. Trends in 651
cognitive sciences, 22(10), 896-910. 652
Kumagai, Y., Arvaneh, M., & Tanaka, T. (2017). Familiarity affects entrainment of EEG in music 653
listening. Frontiers in Human Neuroscience, 11, 384. 654
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
Lehongre, K., Morillon, B., Giraud, A.-L., & Ramus, F. (2013). Impaired auditory sampling in 655
dyslexia: further evidence from combined fMRI and EEG. Frontiers in Human 656
Neuroscience, 7, 454. 657
Lehongre, K., Ramus, F., Villiermet, N., Schwartz, D., & Giraud, A.-L. (2011). Altered low-gamma 658
sampling in auditory cortex accounts for the three main facets of dyslexia. Neuron, 659
72(6), 1080-1090. 660
Lerdahl, F., & Jackendoff, R. (1983). An overview of hierarchical structure in music. Music 661
Perception, 229-252. 662
Lubineau, M., Watkins, C. P., Glasel, H., & Dehaene, S. (2024). Examining the Impact of Reading 663
Fluency on Lexical Decision Results in French 6th Graders. Open Mind, 8, 535-557. 664
Miyagawa, S., Berwick, R. C., & Okanoya, K. (2013). The emergence of hierarchical structure in 665
human language. Frontiers in Psychology, 4, 71. 666
Molinaro, N., Lizarazu, M., Lallier, M., Bourguignon, M., & Carreiras, M. (2016). Out‐of‐667
synchrony speech entrainment in developmental dyslexia. Human brain mapping, 668
37(8), 2767-2783. 669
Müllensiefen, D., Gingras, B., Stewart, L., & Musil, J. J. (2013). Goldsmiths Musical 670
Sophistication Index (Gold-MSI) v1. 0: Technical Report and Documentation Revision 671
0.3. In. 672
Nora, A., Rinkinen, O., Renvall, H., Arkkila, E., Smolander, S., Laasonen, M., & Salmelin, R. 673
(2024). Impaired cortical tracking of speech in children with developmental language 674
disorder. Journal of Neuroscience, 44(22). 675
Obleser, J., Eisner, F., & Kotz, S. A. (2008). Bilateral speech comprehension reflects differential 676
sensitivity to spectral and temporal features. Journal of Neuroscience, 28(32), 8116-677
8123. 678
Obleser, J., & Kayser, C. (2019). Neural entrainment and attentional selection in the listening 679
brain. Trends in cognitive sciences, 23(11), 913-926. 680
Oldfield, R. C. (1971). The assessment and analysis of handedness: the Edinburgh inventory. 681
Neuropsychologia, 9(1), 97-113. 682
Olszewska, A. M., Gaca, M., Herman, A. M., Jednoróg, K., & Marchewka, A. (2021). How musical 683
training shapes the adult brain: Predispositions and neuroplasticity. Frontiers in 684
neuroscience, 15, 630829. 685
Oostenveld, R., Fries, P., Maris, E., & Schoffelen, J.-M. (2011). FieldTrip: open source software 686
for advanced analysis of MEG, EEG, and invasive electrophysiological data. 687
Computational intelligence and neuroscience, 2011. 688
Osorio, S., & Assaneo, M. F. (2025). Anatomically distinct cortical tracking of music and speech 689
by slow (1–8Hz) and fast (70–120Hz) oscillatory activity. PLoS One, 20(5), e0320519. 690
Pantaleo, M. M., Arcuri, G., Manfredi, M., & Proverbio, A. M. (2024). Music literacy improves 691
reading skills via bilateral orthographic development. Scientific Reports, 14(1), 3506. 692
Patel, A. D. (2010). Music, language, and the brain. Oxford university press. 693
Patel, A. D. (2011). Why would musical training benefit the neural encoding of speech? The 694
OPERA hypothesis. Frontiers in psychology, 2, 142. 695
Peelle, J. E., & Davis, M. H. (2012). Neural oscillations carry speech rhythm through to 696
comprehension. Frontiers in Psychology, 3, 320. 697
Pelli, D. G., & Vision, S. (1997). The VideoToolbox software for visual psychophysics: 698
Transforming numbers into movies. Spatial vision, 10, 437-442. 699
Poeppel, D. (2003). The analysis of speech in different temporal integration windows: cerebral 700
lateralization as ‘asymmetric sampling in time’. Speech communication, 41(1), 245-701
255. 702
Power, A. J., Colling, L. J., Mead, N., Barnes, L., & Goswami, U. (2016). Neural encoding of the 703
speech envelope by children with developmental dyslexia. Brain and Language, 160, 1-704
10. 705
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
Retelsdorf, J., Köller, O., & Möller, J. (2011). On the effects of motivation on reading 706
performance growth in secondary school. Learning and instruction, 21(4), 550-559. 707
Schellenberg, E. G., & Lima, C. F. (2024). Music Training and Nonmusical Abilities. Annu Rev 708
Psychol, 75, 87-128. https://doi.org/10.1146/annurev-psych-032323-051354 709
Scott, G. G., Keitel, A., Becirspahic, M., Yao, B., & Sereno, S. C. (2019). The Glasgow Norms: 710
Ratings of 5,500 words on nine scales. Behavior research methods, 51(3), 1258-1270. 711
Smith, Z. M., Delgutte, B., & Oxenham, A. J. (2002). Chimaeric sounds reveal dichotomies in 712
auditory perception. Nature, 416(6876), 87-90. 713
Sohoglu, E., Peelle, J. E., Carlyon, R. P., & Davis, M. H. (2012). Predictive top-down integration of 714
prior knowledge during speech perception. Journal of Neuroscience, 32(25), 8443-8453. 715
Swaminathan, S., Schellenberg, E. G., & Venkatesan, K. (2018). Explaining the association 716
between music training and reading in adults. Journal of Experimental Psychology: 717
Learning, Memory, and Cognition, 44(6), 992. 718
Symons, A. E., Dick, F., & Tierney, A. T. (2021). Dimension-selective attention and dimensional 719
salience modulate cortical tracking of acoustic dimensions. NeuroImage, 244, 118544. 720
Te Rietmolen, N., Mercier, M. R., Trébuchon, A., Morillon, B., & Schön, D. (2024). Speech and 721
music recruit frequency-specific distributed and overlapping cortical networks. Elife, 722
13, RP94509. 723
The Mathworks, I. (2025). MATLAB version: 25.1.0.2973910 (R2025a) Update 1. In The 724
MathWorks, Inc. https://www.mathworks.com 725
Thut, G., Nietzel, A., Brandt, S. A., & Pascual-Leone, A. (2006). α-Band electroencephalographic 726
activity over occipital cortex indexes visuospatial attention bias and predicts visual 727
target detection. Journal of neuroscience, 26(37), 9494-9502. 728
Tierney, A., & Kraus, N. (2013). Music training for the development of reading skills. Progress in 729
brain research, 207, 209-241. 730
Weems, S. A., & Zaidel, E. (2004). The relationship between reading ability and lateralized 731
lexical decision. Brain and Cognition, 55(3), 507-515. 732
Yao, B., Keitel, A., Bruce, G., Scott, G. G., O'Donnell, P. J., & Sereno, S. C. (2018). Differential 733
emotional processing in concrete and abstract words. Journal of Experimental 734
Psychology: Learning, Memory, and Cognition, 44(7), 1064. 735
736
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
Supplemental material 737
Correlation between cortical tracking and word trials in the lexical decision task 738
To test whether envelope tracking during listening to the story and music predicted participants’ 739
behavioural performance in the lexical decision task, we correlated the MI values per electrode 740
and frequency with participants’ median reaction time across word trials. 741
For the story condition, we found one negative cluster, indicating that participants who showed 742
stronger envelope tracking to speech at low frequencies had faster reaction times in the lexical 743
decision task (Figure S1A(i), highlighted in red). This negative cluster (Figure S1A(i)) at low 744
frequencies ( 0.83 Hz - 1.66 Hz), significantly predicted reaction times to words in the lexical 745
decision task (Cohen’s dpeak = 1.79, pcluster < .001, 10 electrodes). We also found a positive cluster 746
at higher frequencies (highlighted in black in Figure S1A(i) spanning 10.05 Hz – 13.26 Hz, (Cohen’s 747
dpeak = 3.11, pcluster < .001, 3 electrodes), indicating that participants who showed stronger 748
envelope tracking to speech had slower reaction times to words in the lexical decision task. 749
For the music condition, we found one negative cluster, similarly indicating that participants who 750
showed stronger envelope tracking to music at low frequencies had faster reaction times in the 751
lexical decision task (Figure S1A(ii). The negative cluster spanned 0.89 Hz - 1.44 Hz, also 752
predicting reaction times (Cohen’s dpeak = 2.32, pcluster < .001, 12 electrodes). We also found one 753
positive cluster spanning 3.55 – 4.37 Hz, indicating that participants who showed stronger 754
envelope tracking to music had slower reaction times in the lexical decision task (Cohen’s dpeak 755
= 1.24, pcluster < .001, 3 electrodes). 756
Correlation between cortical tracking and nonword trials in the lexical decision task 757
To test whether envelope tracking during listening to the story and music predicted participants’ 758
behavioural performance in the lexical decision task, we correlated the MI values per electrode 759
and frequency with participants’ median reaction time across nonword trials. 760
For the story condition we found three clusters, two negative and one positive. The first negative 761
cluster spanned 1.02 Hz - 1.66 Hz (Figure S1B(i) outlined in black), that predicted reaction times 762
in the lexical decision task (Cohen’s dpeak = 1.07, pcluster < .001, 3 electrodes). The second negative 763
cluster at 7.61 Hz – 8.75 Hz also predicted reaction times in the lexical decision task (Cohen’s 764
dpeak = 1.30, pcluster < .001, 3 electrodes, Figure S1B(i) outlined in yellow), indicating that 765
participants who showed stronger envelope tracking to speech at low frequencies had faster 766
reaction times in the lexical decision task. We also found one positive cluster (Figure S1B(i) 767
outlined in green), indicating that participants who showed stronger envelope tracking to speech 768
at low frequencies had slower reaction times in the lexical decision task. The positive cluster 769
spanned 8.75 Hz – 13.26 Hz (Cohen’s dpeak = 1.93, pcluster < .001, 6 electrodes). 770
For the music condition, we found one negative cluster in frontal electrodes at low frequencies 771
(0.89 Hz - 1.35 Hz, Figure S1B(ii)) that predicted reaction times in the lexical decision task 772
(Cohen’s dpeak = 1.84, pcluster < .001, 11 electrodes). This indicates that participants who showed 773
stronger envelope tracking to music at low frequencies had faster reaction times in the lexical 774
decision task. 775
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint
Allen et al. Tracking of speech and music predicts reading ability
776
Figure S1. Relationship between cortical tracking and performance on the lexical decision task, quantified 777
through reaction times for word and nonword trials separately. (A) Channel -by-frequency heatmap of the 778
correlation between cortical tracking of speech and participants’ reaction times to word trials. Two 779
clusters were found for speech, one negative cluster at low frequencies, and one positive cluster at higher 780
frequencies. Two clusters were found for music: one negative cluster indicating that stronger tracking of 781
music correlated with faster reaction times to words in the lexical decision task, and one positive cluster, 782
indicating that stronger tracking of music at higher frequencies correlated with slower reaction times to 783
words. (B) Channel -by-frequency heatmap of the correlation between cortical tracking of music and 784
participants' reaction times to nonword trials. Three clusters were found for speech, two negative clusters 785
indicating that stronger tracking of speech correlated with faster reaction times to nonwords in the lexical 786
decision task, and one positive cluster at higher frequencies (outlined in green), indicating that stronger 787
tracking of speech correlated with slower reaction times to nonwords. One negative cluster was found for 788
music, indicating that stronger tracking of music correlated with faster reaction times to nonwords in the 789
lexical decision task. 790
791
.CC-BY-NC 4.0 International licenseperpetuity. It is made available under a
preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in
The copyright holder for thisthis version posted February 18, 2026. ; https://doi.org/10.64898/2026.02.18.706526doi: bioRxiv preprint