# Dog-human vocal interactions match dogs' sensory-motor tuning

Eloïse Déaux, Théophane Piette, Florence Gaunet, Thierry Legou, Luc H. Arnal, Anne-Lise Giraud

<figure><img src="image_1.png" /></figure>

To cite this version:

Eloïse Déaux, Théophane Piette, Florence Gaunet, Thierry Legou, Luc H. Arnal, et al. Dog-human vocal interactions match dogs' sensory-motor tuning. PLoS Biology, 2024, 22 (10), pp.e3002789. 10.1371/journal.pbio.3002789. pasteur-04916468

HAL Id: pasteur-04916468

https://pasteur.hal.science/pasteur-04916468v1

Submitted on 28 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Dog-human vocal interactions match dogs' sensory-motor tuning

Eloïse C. Déaux<sup>1*</sup>, Théophane Piette<sup>1</sup>, Florence Gaunet<sup>2</sup>, Thierry Legou<sup>3</sup>, Luc Arnal<sup>4</sup>,
Anne-Lise Giraud<sup>1,4</sup>

1 Department of Basic Neurosciences, Faculty of Medicine, University of Geneva, Geneva, Switzerland,

2 Aix-Marseille University and CNRS, Laboratoire de Psychologie Cognitive (UMR 7290), Marseille, France,

3 Aix Marseille University and CNRS, Laboratoire Parole et Langage (UMR 6057), Aix-en-Provence, France,

4 Université Paris Cité, Institut Pasteur, AP-HP, Inserm, Fondation Pour l'Audition, Institut de l'Audition, IHU reConnect, F-75012 Paris, France

* eloise.deaux@unige.ch

<figure><img src="image_3.png" /></figure>

## Abstract

Within species, vocal and auditory systems presumably coevolved to converge on a critical temporal acoustic structure that can be best produced and perceived. While dogs cannot produce articulated sounds, they respond to speech, raising the question as to whether this heterospecific receptive ability could be shaped by exposure to speech or remains bounded by their own sensorimotor capacity. Using acoustic analyses of dog vocalisations, we show that their main production rhythm is slower than the dominant (syllabic) speech rate, and that human-dog-directed speech falls halfway in between. Comparative exploration of neural (electroencephalography) and behavioural responses to speech reveals that comprehension in dogs relies on a slower speech rhythm tracking (delta) than humans' (theta), even though dogs are equally sensitive to speech content and prosody. Thus, the dog audio-motor tuning differs from humans', and we hypothesise that humans may adjust their speech rate to this shared temporal channel as means to improve communication efficacy.

## Introduction

Acoustic communication dynamically evolves as auditory systems are tuned to vocal signals, while in turn vocal production adapts to exploit the capacity of sensory systems [1-4]. In this fine audio-vocal tuning, temporal acoustic features have a universal ecological relevance, being essential, for example, vocal recognition [5,6], predator avoidance [7], or mate choice [8-10].

Production and perception systems can also evolve as a result of interspecific interactions, e.g., in the prey-predator arm race [11] or during interspecific eavesdropping [12]. One of the most prominent and long-term examples of interspecific interactions is that of dogs, *Canis familiaris*, and humans. It is likely that the dog-human cooperation needed to achieve working tasks, i.e., hunting, hauling, and sentinel work required the development of effective interspecific communication skills [13,14]. In terms of vocal production, quantitative and qualitative hypertrophy of bark vocalisations, for example, is believed to be a reflection of the necessity to develop a human-targeted repertoire [15,16]. Concurrently, dog owners speaking to their pet

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Déaux EC, Piette T, Gaunet F, Legou T, Arnal L, Giraud A-L (2024) Dog-human vocal interactions match dogs' sensory-motor tuning. PLoS Biol 22(10): e3002789. https://doi.org/10.1371/journal.pbio.3002789

**Academic Editor:** Attila Andics, Eotvos Lorand University, HUNGARY

**Received:** November 1, 2023

**Accepted:** August 6, 2024

**Published:** October 1, 2024

**Copyright:** © 2024 Déaux et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting information files.

**Funding:** ALG received funding from the NCCR Evolving Language, Swiss National Science Foundation Agreement #51NF40_180888. https://www.snf.ch/en/2P8KtV2RBhgO0eyv/page/researchinFocus/nccr/evolving-language The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
**Competing interests:** The authors have declared that no competing interests exist.

**Abbreviations:** ADS, adult-directed speech; DDS, dog-directed speech; EEG, electroencephalography; ICA, independent component analysis; IQR, interquartile range; LMM, linear-mixed model; PIC, potential for individual coding.

naturally employ accented speech modulations, known as dog-directed speech [17–19]. Notably, dog-directed speech differs from other forms of animal-directed speech, suggesting that humans adjust their speech based on their target audience, instead of having a universal animal-directed speech register [20]. Such an adaptation is also the basis of infant-directed speech that optimises communication with infants’ immature auditory system whose speech processing relies predominantly on the reception of prosodic rhythms that progressively shape faster rhythms [21,22].

Perceptually, dogs exhibit remarkable receptivity to human speech cues [23–26], can learn extensive word repertoires [27], demonstrate fast mapping [28] and statistical learning abilities [29], and possibly word-elicited mental representations [30]. Moreover, behavioural evidence suggests that humans possess a unique sensitivity to dogs’ acoustic signals, as sightless individuals without prior experience with dogs can accurately assess emotional valence from dog vocalisations [31], and humans generally demonstrate greater proficiency in evaluating the emotional valence of dog vocalisations compared to those of other species [32,33]. Collectively, these findings suggest that humans’ accurate perception of dogs’ acoustic signals, even without ontogenic experience, cannot be solely explained by cross-species rules of emotional information transfer, but may reflect interspecific communicative adaptation.

However, one question that remains open is the extent to which dog-human vocal interactions may have been shaped by either species’ production/perception constraints in the temporal domain. Indeed, exploration of the speech system has provided capital insight into the neural bases of the temporal perception/production tuning required for successful intraspecific vocal communication. Speech rhythms are mechanically constrained by the motor effectors, but also operated within a certain dynamic range to best match perception-action neural rhythms. Thus, the dominant speech rhythm, the syllable rate, is cross-culturally stable [34] because it both arises from the interplay of the different articulators [35,36] and corresponds to neural theta oscillations, involved in active sensing across species [37]. In speech perception, the auditory theta rhythm serves to actively interface the acoustics with endogenous neural processes, and the closer the acoustics to this rhythm the more efficient the information transfer. Crucially, the neural theta rhythm can flexibly adapt to speech quasiperiodicity via a mechanism referred to as “speech tracking” [38], and comprehension critically depends on its precision [39–44]. Thus speech production and reception tuning has led to a common temporal window of analysis centred on the 4 to 8 Hz range [45].

However, dogs lack the vocal/neural system necessary to produce articulated communication [46–48], such that they may not have developed the neural machinery needed to perceive theta-based speech signals. While there may have been no specific dog-human adaptation, it is also possible to hypothesise that either the dog’s neural system has adapted to human speech or conversely that humans have adjusted their vocal production to exploit the dogs’ neural (auditory) capacity.

To address these questions, we first analysed dog vocalisations, as well as adult- (ADS) and dog-directed speech (DDS), to probe whether dogs vocalise at the same or at a different rate than humans, and whether the temporal properties of DDS differ from those of ADS. Second, we compared speech neural processing in dogs and humans using noninvasive electroencephalography (EEG), to investigate (1) how dogs track speech modulations; and (2) if, like in humans, dogs’ speech tracking accuracy predicts comprehension. Unlike previous studies, e.g., [24,30,49–51], we selected command words as speech stimuli, which allowed us to use dogs’ behavioural responses as an index of “intelligibility,” while remaining within the structural definition of the DDS register, i.e., short (3 words on average), mostly one-node, imperative utterances [18].
# Results

## Natural vocal rate in dogs and humans

Using 143 vocal sequences (30 dogs) including all major vocal classes (barks, growls, howls, snarls, and whines [15]), 106 adult-directed (27 individuals, 10 women) and 149 dog-directed speech sequences (22 individuals, 16 women) spanning 5 different languages, we found that dogs vocalise at a slower rate than humans (dogs mean ± SD: 2 ± 1.1 vocalisations/s, ADS: 4 ± 1.9 syllables/s; Tukey-corrected post hoc pairwise comparison: t = 6.8, *p* < 0.001, Fig 1A and 1B). We also found that DDS has a slower rate (3 ± 1.6 Hz) than ADS (t = 3.1, *p* = 0.008), but faster rate than the average dog vocal rate (t = 3.9, *p* = 0.006). For a subset of speakers, we

<figure><img src="image_5.png" /><figcaption>**Fig 1.** Comparison of dog/human vocal production. (A) Oscillograms and, overlaid, envelopes used to compute the vocal rate. (B) Model estimates and their 95% CI of vocal rate in dog and human sequences. Black dots are the original observations. (C) Vocal rate (Hz) and mean F0 (Hz) for matched ADS and DDS speech sentences. (D) Model slope and 95% CI of weight effect on dog VR and DF. (E) Density distribution of vocal rate according to vocal classes for dogs and languages for humans. Overall mean (thick dashed line) and SD (thin dashed lines) statistics are displayed. See S1 Data for the underlying data. ADS, adult-directed speech; DDS, dog-directed speech; DF, dominant acoustic frequency; VR, vocal rate.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3002789.g001
found duration-matched DDS and ADS sentences, allowing us to confirm that pet owners slow their speech rate when talking to their dogs (paired *t* test: *t* = 2.7, *df* = 11, *p* = 0.02, Cohen's *d* = 0.8, Fig 1C). DDS also has higher mean F0 than ADS (*t* = -2.2, *df* = 11, *p* = 0.05, Cohen's *d* = 0.6) confirming previous results [17,19]. Further analyses of vocal sequences returned no significant differences in vocal rate among vocal classes in dogs (*F*<sub>4,16.1</sub> = 1.4, *p* = 0.28) nor among languages in both speech types (*F*<sub>4,31.8</sub> = 2.1, *p* = 0.11, Fig 1E). Thus, the dog's vocal rate is overall slower than human speech and importantly, pet owners modify not only the spectral but also the temporal feature of their output when speaking to their dogs, in a direction that brings them closer to the natural vocal rate of the latter.

Furthermore, when exploring other factors known to influence the structure of animal vocal signals [2], we found no evidence of large inter-individual differences in vocal rate unlike for the dominant acoustic frequency (S1 Table) confirming the latter's functional significance in individual discrimination [52,53] and speaking against such selection effects in the former. Concurrently, body weight had no explanatory effect on vocal rate variation (F<sub>1,11.41</sub> = 0.04, *p* = 0.8) while it was inversely related to dominant acoustic frequency (F<sub>1,12.07</sub> = 6.03, *p* = 0.03) confirming the known acoustic allometric relationship between body weight and spectral parameters [54] and speaking to other types of constraints on vocal rate (Fig 1D) [1].

## Neural tracking and speech “intelligibility”

To investigate auditory neural processes in dogs, we adapted typical human protocols, e.g., [39,41,43], where speech intelligibility is altered using spectral and temporal modifications of speech stimuli and neural tracking strength is correlated to behavioural measures of intelligibility (Fig 2). Speech streams were composed of words that the dogs had learned to respond to, i.e., command words (e.g., “sit,” “come”). We used command DDS rather than praising DDS (e.g., “oh that’s a good boy!”) to be able to obtain an objective index of “intelligibility” in the sense of a successful stimulus–action relationship, assessed during the behavioural task. In humans, comprehension was measured by asking participants to rate word streams on an intelligibility scale. We performed EEG and behavioural experiments on 12 dogs (1 to 13 years old, 7 females) and 12 paired human participants (18 to 65 years old, 6 women) with no self-reported hearing deficits. Four dogs and 1 human participant were excluded from analyses due to poor EEG signal quality.

We first confirmed that modifying speech spectral and temporal features altered both species' perceptual performances. Modifying speech rate (main effect: F<sub>2,80</sub> = 46.6, *p* < 0.001) and type (main effect: F<sub>2,80</sub> = 112, *p* < 0.001) affected speech intelligibility in humans, in an interactive way (speech rate by speech type interaction: F<sub>4,80</sub> = 10.5, *p* < 0.001). When content was removed (i.e., the prosody-only condition), participants failed to understand the speech sentence at all speech rates. In the other 2 conditions, increasing speech rate decreased speech intelligibility (see Fig 3A for post hoc significance pairwise tests). In dogs, speech rate (main effect: F<sub>2,64</sub> = 4.9, *p* = 0.01) and type (main effect: F<sub>2,64</sub> = 6.4, *p* = 0.003) also impacted speech intelligibility, again interactively (speech rate by speech type interaction: F<sub>4,64</sub> = 6.9; *p* < 0.001), with intelligibility dropping as speech rate increased, but only in the normal speech type condition (Fig 3A).

We then quantified the 2 species' neural responses, restricting the EEG analyses to the FCz electrode in humans as it showed the strongest response to the acoustic stimulation (S1 Fig) and Cz in dogs, known to capture auditory stimulus-locked activity [55,56]. Both dogs and humans showed increased power activity (relative to the pre-stimulus baseline period) in the low frequency range (<10 Hz, Fig 3B), confirming and characterising the auditory cortex activity reported in fMRI studies of dog speech processing [49,50,57]. However, we noted a
<figure><img src="image_6.png" /><figcaption>**Fig 2. Schematic of the perception study.** Word stream stimuli were first constructed by recording dog-specific command words (mostly disyllabic and monosyllabic, cf. small insert) that were appended into a 5-word stream with approximately 300 ± 50 ms silence intervals. These word streams were altered with regards to (1) speech type: by either removing content (reversed words) or prosodic information (flattened pitch modulation and reversed energy contour); and (2) speech rate: compression by a factor of 2 or 4; amounting to 9 word-stream conditions in total. The behavioural experiment consisted of an intelligibility scoring task for humans who listened to the full word stream, and of a playback task for dogs, who heard each word command separately (45 in total) a maximum of 3 times each, while the experimenter and the owner agreed on a behavioural response score. For the EEG experiment, dogs were first fitted with 1 to 4 electrodes covered by a headband and linked to an amplifier strapped on their back (photo inserts). They were then instructed to lie down and passively listen to an audio track (broadcasted via a speaker) containing 40 repetitions of each word stream condition. For comparability purposes, human EEG recordings were made under the same experimental conditions (see also S2 Data). Photo credit: E. Déaux. EEG, electroencephalography.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3002789.g002

first difference between the 2 species' neural responses in this frequency range. Dogs showed a predominant power increase in the delta band (1 to 3 Hz), as opposed to the theta band (4 to 7 Hz) in humans (Mann-Whitney U test: U = 70, df = 17, *p* = 0.03, *r* = 0.5, Fig 3C), speaking to possible divergent auditory processes.

Given the presence of a stimulus-related and sustained neural response, we then probed whether dogs display evidence of a speech tracking response under normal speech conditions. Cerebro-acoustic coherence, a measure that quantifies the phase-locking of neural signals to speech envelope [41,42], was above the mean random coherence value throughout the 1 to 10 Hz range in humans, but restricted to a 1 to 3 Hz peak in dogs (Fig 4A). Averaged values in the delta band were significantly higher in the real cerebro-acoustic than in the cerebro-randomised acoustic pairings in humans (paired t test: t = -5, df = 10, *p* < 0.001, d = 1.5) and in dogs (t = -3, df = 7, *p* = 0.02, d = 1.1). However, theta cerebro-acoustic coherence was significantly higher than in the cerebro-randomised acoustic pairings in humans (t = -2.87, df = 10, *p* = 0.02,
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Speech stimulus alteration effects on intelligibility and characterisation of neural responses. (A) Mean (± SE) behavioural responses according to speech type and rate in humans and dogs. Tukey-corrected, post hoc pairwise comparisons are shown. ***p < 0.001, **p < 0.01, *p < 0.05. (B) Time-frequency plots averaged across all conditions and individuals within species. Z-score transformed relative power is plotted to ease visual comparison across species. (C) Power spectra (peak normalised and averaged between 0 and 1.3 s) and unpaired *t* test between species on frequency of highest power (range = 1–7 Hz). For the underlying data, see S2 Data.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3002789.g003

d = 0.9) but not in dogs (t = -0.6, df = 7, *p* = 0.5, d = 0.2). In other words, dogs show evidence of auditory tracking capabilities, as do other species [58–60]; however, in the context of speech stimulation, and unlike in humans, such tracking is restricted to the delta band (Fig 4B).

Having found evidence for a general speech tracking process in both species, albeit at differ-
ing oscillation bands, we then focused on determining how flexible tracking was within these
frequency ranges and whether, similar to humans, dogs' cortical tracking was linked to speech
intelligibility. As word streams differed among individuals, the specific peak rhythm within
each of these frequency ranges (i.e., delta and theta) also differed. As these peak rhythms acous-
tically represent specific granularity levels (i.e., word versus syllable levels [45]), we first
extracted these word and syllable rates and then computed the corresponding cerebro-acoustic
coherence value specific to each participant (hereafter referred as to "word cerebro-acoustic
coherence" and "syllable cerebro-acoustic coherence"). Indeed, if speech comprehension is
linked to tracking strength within these bounded ranges, then as speech rate increases, both cor-
tical tracking and speech comprehension should be negatively impacted. We first confirmed
that increasing speech rate had a negative effect on cerebro-acoustic coherence in both species
and at both granularity levels (S2 Fig). In humans, both syllable and word cerebro-acoustic
coherence decreased as syllable rate (F₁,₉₀.₅ = 9.2, *p* = 0.003, S2A Fig) and word rate (F₁,₈₄ = 4.2,
A

<figure><img src="image_8.png" /></figure>

B

<figure><img src="image_9.png" /><figcaption>**Fig 4.** Speech tracking in both species and both the delta and the theta bands. (A) Mean (± SD) cerebro-acoustic coherence over the 1–20 Hz range, calculated from the normal speech condition. Black dashed line shows mean (± SD) random coherence values for pairings of neural signals with randomised acoustic envelopes. (B) Paired *t* test of coherence in the delta and theta range between the real cerebro-acoustic and cerebro-randomised acoustic pairings (see also S3 Data).</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3002789.g004
<figure><img src="image_10.png" /><figcaption>**Fig 5.** Stronger speech neural tracking (cerebro-acoustic coherence) is linked to increased comprehension. (A) Model slopes and 95% CI for the syllable and word cerebro-acoustic coherence effect on intelligibility in dogs and humans. (B) Mean and SE intercepts for each speech type in humans and dogs, showing that beyond speech tracking, additional processes must be present to explain the differences (see also S2 Data for numerical data).</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3002789.g005

$p = 0.04$, S2B Fig) increased respectively, but speech type had no effect in either model (Syllable model: $F_{2,83} = 0.27$, $p = 0.8$; Word model: $F_{2,83} = 0.5$, $p = 0.6$). The same pattern was found in dogs, with both syllable and word cerebro-acoustic coherence dropping with increasing syllable rate ($F_{1,65.5} = 5$, $p = 0.03$, S2C Fig) and word rate ($F_{1,61} = 9.1$, $p = 0.003$, S2D Fig) respectively, while speech type (Syllable model: $F_{2,59} = 1.5$, $p = 0.2$; Word model: $F_{2,59} = 2.9$, $p = 0.07$) had no effect.

Remarkably however, the 2 species differed with regards to the granularity level at which tracking was most strongly related to behavioural outputs (Fig 5A). Specifically, in humans, word cerebro-acoustic coherence did not explain intelligibility ($F_{1,58.3} = 1.3$, $p = 0.3$) while stronger syllable cerebro-acoustic coherence led to increased intelligibility ($F_{1,89.5} = 5.5$, $p = 0.02$). Conversely in dogs, syllable cerebro-acoustic coherence had no impact on speech intelligibility ($F_{1,62} = 0.21$, $p = 0.6$), while intelligibility increased with stronger word cerebro-acoustic coherence ($F_{1,62} = 4.69$, $p = 0.03$). Interestingly, in both species, the speech type (normal, content-only, prosody-only) main effect remained (humans syllable model: $F_{2,85} = 49.9$, $p < 0.001$; dogs word model: $F_{2,61} = 3.2$, $p = 0.05$, Fig 5B), with significant differences among the intercepts of all speech types in humans (all pairwise comparisons: $p < 0.001$) and higher intelligibility in the normal speech condition compared to the prosody-only condition (Norm. speech–Prosody-only est = 0.1, df = 69, $t = 2.4$, $p = 0.05$) in dogs (all other pairwise
comparisons: *p* > 0.05). In other words, like humans, dogs' comprehension of speech appears
to involve more than stimulus-driven auditory processes [45,56,61].

Discussion

While humans' main speaking rate is centred around 4 ± 1.9 syllables/s (i.e., theta band) [34,62], we show here that dogs' vocal rate is slower centred around 2 ± 1.1 vocalisations/s (i.e., delta range). This rate is conserved across vocalisation types, not influenced by body weight and shows only limited inter-individual differences. To ascertain this negative result was not caused by inadequate data, we used the same sample to probe the well-known individual-related differences in dominant frequency and its allometric relationship to body weight (Fig 1) [1,2,63]. These results suggest that despite variation among calls in length and repetition rate, which may convey context-related information, e.g., [64], dogs exhibit a generic ~2 Hz vocal rhythm. While this rhythm differs from the human dominant speech rhythm, these results are consistent with what is known from the speech production system, namely that despite wide linguistic and spectral variations, vocal rates show remarkable consistency [34,62]. It should be pointed, however, that while we tested several important factors known to influence signal variability other variables, such as subject sex or dog breed could not be explored but may yet be shown to also have an impact on vocal rate.

Interestingly, the theta vocal rhythm is also present in many primate species including both
closely and more distantly related ones [65-67] leading to hypotheses of an exaptation from
masticatory movements [68]. Yet, despite being masticators, dogs do not vocalise in that range
and instead exhibit a lower rate, which suggests that the theta rhythm evolved sometime after
the split between the Laurasiatheria and Euarchontoglires, and begs for a more thorough char-
acterisation of the phylogeny of the theta rhythm and of the possible selective forces behind its
emergence. Nonetheless, although further replication of these results on an independent sample
would be welcome, dogs, like humans, appear to be subject to evolutionary factors that have
kept temporal production constrained within their species-specific ranges (Fig 1). Such factors
could include inherent differences in the function(s) of their respective communication systems
[4] and/or, as the results of our perception experiment suggest, neural constraints [69].

We found that dogs exhibit cortical tracking of acoustic stimuli, confirming its nature as a
basal processing mechanism [58,59,70]. However, they do not syllabify speech as humans do,
but instead primarily track slow amplitude modulations via delta oscillations. Importantly, this
slower, word-level tracking does not mean that dogs only process speech via prosodic cues, as
we also show that eliciting successful responses requires the full integration of prosodic and
content information (Figs 3 and 5). While delta oscillations are also important in human
speech encoding and decoding [71], current models of speech perception place the theta-sylla-
ble rhythm as the centre of this processing machinery [45,72]. Thus, that dogs do not process
speech via theta oscillations but instead primarily rely on delta oscillations, and that this neural
rhythm matches their natural vocal rate, supports the hypothesis that dogs' production/per-
ception tuning is centred around the delta rhythm and has been exapted in the context of
human-dog vocal communication. However, what remains to be demonstrated in order to
test further this hypothesis is whether delta tracking does indeed support conspecific signal
processing. Furthermore, one caveat for these results is that we used concatenated word
streams for analytical purposes, and the unnaturalness of the stimulus could have had an
impact on dogs' neural responses. Yet from human research, we know that unnatural stimuli,
such as backward speech, suppress tracking [73], such that had the "unnaturalness" of the
stimulus strongly affected dogs' neural responses we would have likely failed to detect cortical
tracking.
Importantly, while cortical tracking plays a critical role in speech comprehension in humans [39,41,45,73], this study is the first to establish such a link between neural processes and behavioural outputs in another species. Indeed, while previous research in dog neuroscience has revealed interesting parallels between human and dog neural processes, they typically used speech stimuli that did not allow to objectively measure dogs' behavioural responses and as such had to assume the functionally adaptive significance of their neural data, e.g., [30,49]. By electing to use commanding DDS (e.g., "sit," "come") rather than praising DDS (e.g., "good boy!"), we could objectively quantify stimulus intelligibility and provide conclusive evidence that cortical tracking is a prerequisite of successful acoustic communication not only within but also across species.

Our results additionally lend support to the hypothesis that humans have adjusted their vocal production to exploit the dogs' neural (auditory) capacity. We find that DDS has a slower rhythm than ADS that brings it closer to the dog's production rate. That signal design is tuned to exploit receivers' neural capacities is well established both within species [3] and between species [74] and even professionals across cultures, spontaneously make use of existing sensory biases when training animals [75]. As such, and considering the evolutionary cooperative bond between humans and dogs [13,14], it is plausible that speech temporal tuning to dogs' sensory systems would have led to more efficient communicative interactions and thus been selected for [76,77]. However, what remains to be established is whether this temporal tuning results from the happenstance of an overgeneralisation of infant-directed speech, possibly due to perceived neotonic traits or limited language abilities [20,78] or truly specific to dogs. Future research could aim to investigate interspecific production and perception tuning between humans and other species, to establish whether humans specifically adjust to the neural capacities of the listener or if it is a more general process. If this vocal adjustment is due to the specific human-dog working bond, we may predict that slow speech should be stable cross-culturally and exhibited even more so in vocal interactions that involve directing working dogs.

One possible limitation of this study is that in the production experiment, DDS consists mainly of sentences (e.g., "what a good boy!", "wanna play outside?") while the perceptual experiment used commanding DDS. One could thus argue that they may not have the same acoustic characteristics and thus not be directly comparable. However, when constructing acoustic stimuli, we ensured that the vocal rate would be within the natural DDS range (i.e., 3 ± 1.6 Hz) and complementary analyses revealed no differences in spectral characteristics between the 2 registers (S3 Fig). Furthermore, like previous studies that used praising DDS as stimuli [49,51], we found that eliciting successful responses required the full integration of prosodic and content information (Fig 3). This suggests that commanding and praising DDS may be similarly processed by dogs at least at the auditory level which we explored here. Nonetheless and while it was beyond the scope of this study, it would be greatly beneficial for future comparative work to do a formal acoustic analysis of the 2 types of speech to better understand their similarities and differences and what this may entail for dog perception.

Finally, and for both species, we found that differences in the intercepts among speech types remained even after accounting for the speech tracking effect (Fig 5), suggesting that other processes also influence intelligibility. In humans, cortical tracking is a part of a complex bottom-up and top-down hierarchy known to contribute to comprehension [45,79–81]. Most notably, the hierarchical phase-amplitude coupling of theta and gamma frequency bands allows for phoneme encoding [38,82] while top-down processes involving motor cortex activity [83] are causally related to perception [84]. Such processes have not yet been uncovered in dogs, but have been described in primates [85–88]. Additionally, a recent study in dogs mentioned differential activation in a cortical premotor region when comparing familiar versus unfamiliar language processing, hinting at a possible top-down process [61]. Thus, it will be
greatly interesting for future research to investigate whether and to what extent hierarchical
bottom-up and predictive processes initially linked to within-species acoustic processing can
adapt to or constrain interspecific communication.

Overall, these results reveal that dogs' auditory and vocal systems have aligned on a single temporal processing window that differs from that of humans, and which remains predomi- nant even when dogs process and appropriately respond to human speech. In parallel, we show that humans who speak to their dogs adopt a speech rate that differs from adult-directed speech and more closely aligns with the dog's neural delta oscillatory capacity. These data are consistent with the hypothesis that in the history of the dog-human relationship, the neural constraints of the dogs' reception system may have limited this heterospecific communication to a temporal structure falling midway between the natural speech rate and a slower rate that would perfectly match the dog's analysis capacity. However, future research on different ani- mal-directed speech registers particularly as it pertains to their temporal characterisation is needed before we can fully establish the extent to which humans are sensitive to their target's perceptual constraints and whether the dog-human relationship is as special as it may appear.

**Methods**

Ethics statement

All the dogs used in this study were pet dogs who lived with their caregivers. As the tests took
place in France and involved noninvasive EEG recordings and behavioural tests, no ethical
approval was required under the French law. The human participants all provided informed
written consent prior to the experiments and the procedures were approved by the ethics com-
mission of Geneva University CUREG.202011.18.

Subjects

**Dogs:** Dog owners were recruited by contacting canine clubs located in France. After initial contact with potential participants, dogs were recruited if, based on their owners' reports, they met the following inclusion criteria: being 1 year or older, having no hearing deficits, a good sociability level, high trainability, and a good level of education. This recruitment process resulted in a pool of 12 dogs (7 females) aged 1 to 13 years old being included, all being medium to large dog breeds, the smallest being Shetland sheepdogs and the biggest being the Beauceron. While all dogs had had basic obedience training, their primary canine sport varied. Half of the dogs were doing obedience training, 4 participated in dog dancing classes and the last 2 animals were training in mantrailing.

**Humans:** We also recruited the same number of human participants (6 women) from the clubs who served as paired controls. Inclusion criteria for human participants were: being aged 18 to 65 years old, having no self-reported hearing deficit, no psychiatric or motor disorders, and speaking French fluently.

Procedure

**Perception experiment.** **Preexperiment dog training:** We developed a training protocol using positive reinforcement and behavioural shaping to condition dogs to wear the EEG equipment while remaining still. First, dogs were clicker-trained to lay down while resting their head. Once dogs could maintain this position for at least 15 s, they were habituated to wear a headband (happy hoodie, Zony Pets, United States of America) normally used during toileting, to which we made holes to let the ears out (Fig 2, photo inserts). Then, they were finally conditioned to maintain the position, while wearing the headband and listening to a
variety of noises including music, environmental noises, and voices. Dogs were judged suffi-
ciently experienced once they could maintain this position regardless of noise or other envi-
ronmental disturbances for at least 15 s. Throughout the training, dogs were monitored for
signs of stress, based on the well-known behavioural markers: yawning, excessive panting,
and/or lip licking, as well as body and ear positions. Had a dog exhibited these signs, their
training would have been stopped and they would have been excluded from the study.

**Acoustic stimuli:** Typically, comprehension is assessed in humans by asking participants to rate word sequences on an intelligibility scale. As this was not possible for dogs, we selected words that the dogs had learned to respond to, i.e., command words (e.g., “sit,” “come”), allowing us to use behavioural responses to these words as an index of “comprehension” in the sense of a successful stimulus–action relationship. For each dog, we recorded 5 command words spoken by their owners during a typical training session to obtain original, naturalistic DDS. The words were a mix of mono- and disyllabic words (Fig 2). Each dog listened (EEG task) and responded (behavioural task) to their specific set of command words. Their matched control human participant also listened to the same stimuli. Recordings were made with a Sennheiser ME64 microphone and a K6 module mounted onto a FOSTEX FR-2LE field recorder in 44.1 kHz—16 bit wav format. One exemplar of each command word was selected based on the sound quality and on whether that occurrence resulted in a clear, successful behavioural response. The selected command words were first high-pass filtered at 100 Hz, independently normalised at –2 dB and then concatenated into 1 word stream with 300 ± 50 ms silent intervals in between command words. The decision to concatenate command words was methodological, in the sense that to perform the coherence analysis (see the cerebro-acoustic coherence section below) long stimuli were required (at least 1 s). As a control analysis, we used a Kruskal–Wallis test to compare the mean F0 (H = 9.6, *n* = 8, df = 2, *p* = 0.008) and interquartile range (IQR) of F0 (H = 12.2, *n* = 8, df = 2, *p* = 0.002) of the command streams to French ADS and DDS registers. Post hoc Dunn tests with a Bonferroni correction confirmed that for both measurements, there was no significant difference between the command stream and DDS (*p* > 0.05), while ADS had significantly lower mean and IQR F0 than the other 2 speech registers (*p* < 0.05, S3 Fig).

We then used PRAAT and the VocalToolkit plug-in to construct the acoustic stimuli. In
total, we constructed nine-word stream stimuli using a fully crossed design of the 3 levels of
speech type and the 3 levels of speech rate (Fig 2). In the Content-only condition, we first
changed the pitch median of the original dog-directed word sequence to match that of the
owner's adult-directed speech pitch, and to remove all pitch modulations. In a second step, we
altered the intensity component of prosody, by reversing the natural intensity contour, while
keeping the speech forward. To create the Prosody-only condition, we reversed each individual
word rendering the speech unintelligible, while keeping their order in the word stream.
Because this process also reversed pitch modulation and intensity, we then copied the pitch
and intensity modulation patterns from the original speech word stream, effectively reinstating
the original prosody. Finally, given that these procedures resulted in undesired contingent
effects, such as slightly robotized voice effects, we also created a control Normal speech condi-
tion by first making the pitch monotone and recopying the original pitch contour from the
original recording. This ensured that these contingent effects were also present in the Normal
word stream and thus controlled for.

To accelerate the speech rate, we used the “change tempo” function in Audacity, https://
audacityteam.org/, which accelerates the rate without impacting the pitch. Each stream was
compressed by a factor of 2 (twice as fast) and a factor of 4 (4 times as fast). This resulted in a
total of 9 word streams to which we affixed the dog’s name in its original form, as a way of cap-
turing the dog’s attention throughout the experimental session. We then created experimental
tracks that included all 9 word streams repeated 40 times each, presented in a random order and separated by an inter-stimulus silent interval of 1.5 ± 0.5 s. Experimental tracks lasted on average 23 min.

**Experimental location:** All tests took place at the owner's home whenever possible, or at another place familiar to the dog (such as another participant's house belonging to the same canine club). This avoided having to familiarise animals with new locations and gave us more flexibility during the COVID-19 situation. Typically, this involved using the living room area of the house, with the dog being either positioned on a bedding or a couch, depending on its usual place during the preexperiment training.

**EEG listening task:** On the day of the EEG test, dogs were fitted with 1 to 4 golden cup electrodes (at least Cz and if possible, C3, C4, and POz) using gel and a conductive paste, and connected to a g.Nautilus amplifier (g.tec medical engineering GmbH, Austria) secured on the back of the dog, which wirelessly transmitted data to a receiver connected to a recording DELL laptop. For 8 dogs (out of the 12), only 1 electrode could be positioned, such that we selected Cz as it is the most reliably located [89] and is known to show acoustic stimulus-locked responses [55], as our results confirmed. The reference electrode was placed at the nap of the neck (Fig 2, photo inserts). Electrodes were then secured by the headband to prevent any movement during the experiment. Electrode impedance was kept under 30 kΩ and data were recorded at a 500 Hz sampling rate. Dogs laid down facing a PREMIO 8 speaker (T.A.G Montpellier, France) placed 2 m away. The experimental track was then broadcast at 60 ± 5 dBC. The experiment was paused regularly to reward the dog for maintaining the position or when the dog became restless. On average the dog EEG listening task lasted 39.6 ± 15.8 min.

Human recordings were made as similar as possible, using the same recording device and the same set-up. The only differences being that we used 7 to 8 gel-based g.SCARABEO (g.tec medical engineering GmbH, Austria) active electrodes (FCz, AFz, CP3, CP4, CPz, FC3, FC4, and POz) inserted in a cap and ear-referenced, and that the participants were asked to sit in a chair and instructed to avoid movements and blinking during the stimulus presentation. No breaks were given during the presentation.

## Behavioural task

**Humans.** Participants were asked before the EEG listening task, to score the linguistic material. For that, they were equipped with headphones and listened to each stimulus and were prompted to score on a scale of 0 to 5 how many words they understood. The word streams were randomly ordered but only presented once to avoid learning effects.

**Dogs.** To obtain a comparable index for dogs, we used a playback experiment where dogs were made to listen to each command word separately (45 words in total) and scored on how well they responded to the command. To do so, we installed the speaker at the mouth level of the dog owner, who stood quietly next to it while wearing sunglasses and a face mask, holding their arms along their body or behind their back. This procedure ensured that the experiment was as realistic as possible while preventing dogs from using visual cues to answer the command. Prior to each command word, the dog was positioned in front of the speaker 1 to 2 m away in a position that allowed it to display the appropriate behavioural response (e.g., standing up if the next command was a "sit"). Each command word was played a maximum of 3 times with 10 s of silent interval in between. The first time the word command was played, it was preceded by the name of the dog, to grab her attention and replicate typical training settings. After the command word was played, we scored on a scale of 1 to 5 how accurately the dog responded to the command (Table 1). If the dog obtained a score of 4 or 5 (i.e., perfect response within the 10 s scoring interval), she was rewarded with her usual treat, the playback
<table><caption>Table 1. Scoring scale for dogs' behavioural responses to command words.</caption><thead><tr><td>Score</td><td>Behavioural response</td></tr></thead><tbody><tr><td>5</td><td>Complete response within 5 s</td></tr><tr><td>4</td><td>Complete response after 5 s and before 10 s</td></tr><tr><td>3</td><td>Incomplete response</td></tr><tr><td>2</td><td>Nonspecific or wrong response</td></tr><tr><td>1</td><td>No response within 10 s</td></tr></tbody></table>

https://doi.org/10.1371/journal.pbio.3002789.t001

series for that command was stopped and we moved on to the next command word series, again mimicking a typical training session. Scoring was performed by the experimenter and the dog owner. If the two disagreed, the lower response score was given. If the dog became restless and/or inattentive, the experiment was interrupted by a play and/or walk session and then resumed. On average, the task lasted 53.5 +/- 20.3 min.

## Production experiment

**Dogs.** We collected vocal sequences from YouTube videos using the freely available Audio Set database [90]. A total of 143 sequences (30 individuals) lasting >1.5 s were extracted, spanning the range of basic vocalisations in canids: barks (n = 54, 38%), growls (n = 18, 13%), whines (n = 21, 15%), snarls (n = 17, 12%), and howls (n = 33, 23%) [15]), with vocalisation classification being done by ear. If the sequences contained more than 1 vocal type, i.e., were “mixed sounds” (which concerns approximately a third of the recordings), we classified these according to the vocalisation that was most present in the sequence. Inter- and intra-observer reliability in vocal classification was assessed using the Kappa measure, based on a random sample of 50 recordings (approximately a third of all recordings). The values obtained, i.e., inter-observer (between 2 researchers, ED and TP) Kappa = 0.8 and intra-observer (with a 6-month interval between the 2 classifications) Kappa = 0.9 showed very high reliability in both cases.

We categorised the dogs according to their body size as either small (a terrier-like dog or below) or large and to their age class (adult versus juvenile). Whenever available, we recorded the breed of the dog and obtained the corresponding mean breed weight using the American Kennel Club website (https://www.akc.org/). For the Cane corso, data were unavailable on the AKC website, so we used the French equivalent, the Société Centrale Canine website (https://www.centrale-canine.fr). Finally, for those 3 individuals whose breed was known and who were pups, we first estimated the age (in months) of the pup from the video and then used the weight curve of the corresponding weight category provided in [91] to obtain the mean weight (50% centile) at that age.

**Humans.** To keep data sets as comparable as possible, we extracted ADS and DDS sequences from YouTube videos (ADS: 106 sequences, 27 individuals, 10 women; DDS: 149 sequences, 22 individuals, 16 women). We selected speech sequences from 5 different languages: English, French, Italian, Japanese, and Vietnamese to cover the range of stress-, syllable-, and mora-timed speech patterns. DDS sentences included both typical praising and command utterances. For 12 individuals (9 women) that produced DDS sequences, we were able to match 1 DDS and 1 ADS exemplar (matched for duration), either extracted from the same video or by looking at other videos published by that user. For this analysis, we were not able to find matching ADS and DDS sequences in Vietnamese.
## Measurements

**Perception experiment.** **Intelligibility index:** For humans, the intelligibility score corresponded to the proportion of correctly comprehended words. To obtain a comparable index for dogs, we calculated the mean response score from the maximum behavioural score obtained in response to each command word of a given condition (Table 1), and then scaled this variable between 0 and 1.

**Audio signal:** We computed the speech envelope (from the onset of the first command word) using the Hilbert transform, low-pass filtered below 30 Hz using an eighth order Butterworth filter, to extract, for each participant, the word and syllable rate for each of the 9 word streams from the power spectrum of the envelope. These word and syllable rate variables were then z-scored and subsequently used as regressors in the statistical analyses investigating the relationships between neural, acoustic, and behavioural data.

**EEG data:** All EEG preprocessing steps were done in MATLAB using the fieldtrip toolbox [92] and custom-written scripts. EEG data were bandpass filtered between 1 and 70 Hz and a DFT filter was applied at 50, 100, and 150 Hz. Signals were then epoched from 1 s pre-stimulus onset to the end of the word sequence. Human data were re-referenced to average and an independent component analysis (ICA) was used to remove eye blink data. ICA was not used on dog data, as for most subjects (8 out of 12), we only had the Cz recording electrode. Artefact rejection (eye blinks, muscle, and jumps) was automatically done using fieldtrip functions, with species-specific cut-off z-values (more stringent for humans). A final visual inspection of all trials was used to remove any other trial that failed the rejection procedure. During these initial procedures, we had to exclude 4 dogs and 1 human participant due to poor signal quality, leaving 8 dogs and 11 humans for the analyses. On average 27.35 ± 5.48 trials were kept per condition in dogs and 35.64 ± 2.84 trials in humans (S2 Table and S4 and S5 Figs).

**Electrode selection.** For dogs, the EEG analyses had to be restricted to the Cz electrode, as it was the only one available for most of them. Thus, we decided to similarly restrict further data analyses to one electrode for humans. To select which electrode to keep, we used a decoding approach, using the mTRF model [93]. Briefly, mTRF models use regularised linear regression to find the latent relationships between the stimulus features (in our case the speech envelope) and the neural response. We ran mTRF models for each subject and each electrode separately, restricting the shifting lag from 100 ms pre-stimulus onset to 500 ms post-stimulus onset. We then calculated the correlation between the reconstructed and the actual stimulus and saved the mTRF *r* value obtained as our measure of how well each electrode responded to the task. A linear-mixed model with electrode as a fixed effect and subject ID as a random term, followed by post hoc analyses showed that FCz had a significantly higher correlation value compared to the other electrodes, and was thus selected for further analyses (S1 Fig).

**Cerebro-acoustic coherence.** To assess the extent of cortical phase-locking to the speech temporal structure, we used the cerebro-acoustic coherence index. Focusing on the control normal speech condition, we first obtained the cross-spectral density between neural signal and the speech envelope using a wavelet method between 1 and 20 Hz in 0.1 Hz frequency steps and 0.01 ms time steps, from 0.6 s post-stimulus onset to 1.3 s. This time window was selected to exclude ERP components resulting from the first word of the sentence, which was always the dog's name, and to allow keeping trial length equal across subjects. We then used the coherence function in fieldtrip to compute the phase coherence between the speech envelope and neural signal. To evaluate how well subjects tracked the speech signal, we compared the actual coherence to random coherence values obtained from the pairings of neural data with randomised acoustic envelopes averaged over 100 runs. To further characterise neural
tracking in the 2 most relevant auditory frequency bands, i.e., delta and theta bands, we extracted mean coherence values (delta: 1 to 3 Hz; theta: 4 to 7 Hz) in both real and random data sets and compared them using paired *t* tests. Then, to explore how tracking was influenced by speech type and rate and how it related to behavioural data, we calculated, for each subject in each condition, the mean word and syllable cerebro-acoustic coherence (time window: 0.6 to 1.3 s post-stimulus onset, time steps: 0.01 s, frequency steps: 0.5 Hz) value centred around the subject-specific stimulus word and syllable rate (+/− 0.5 Hz).

**Production experiment.** **Vocal rate and dominant acoustic frequency.** Acoustic analyses were performed using the seewave package in R [94]. To extract the peak vocal rate, i.e., the predominant rhythm at which vocalisations in a sequence are produced, we first bandpass filtered the sequence between 0.1 and 10 kHz and then computed the signal's envelope using the Hilbert transform. This envelope was further low-pass filtered below 20 Hz using a fourth order Butterworth filter and a wavelet method was used to obtain the frequency decomposition of the signal and extract the frequency of the highest peak. This method was used rather than the more traditional approaches that are based on counting call units and/or on inter-call intervals, because it makes no assumption with regards to the underlying process that produces changes in the amplitude envelope and is thus more adequate when looking at acoustic signals that vary greatly in length and structure.

As a control analysis, we also extracted the dominant acoustic frequency of one vocalisation per sequence (selected based on its signal-to-noise ratio) for the dogs and the sentence's mean fundamental frequency (F0) for humans. For the dogs, the vocal unit was first bandpass filtered between 50 Hz and 2 kHz, and the averaged frequency spectrum was then computed to extract the frequency of the peak amplitude. We focused on the dominant acoustic frequency rather than the fundamental frequency, because the latter is not always quantifiable, particularly in noisy and chaotic vocalisations such as barks. For humans, we used PraAT (with standard settings) to extract the mean F0 in each sequence, as in human speech, F0 is both easy to compute and a better characterisation of pitch than dominant frequency. Pitch values were visually inspected on the spectrograms before extraction to ensure accurate measures.

**Potential for individual coding (PIC).** Among the numerous selective pressures that can impact signal structure, one is the need for increased individual recognition [63]. To assess whether such a process could explain dog vocal rate variation, we calculated the potential for identity coding (PIC) index of each vocalisation for both the vocal rate and dominant frequency parameters. We calculated the within- and between-individual coefficients of variation (CVw and CVb, respectively) using the formula for small samples and obtained the feature's PIC value, which is the CVb/meanCVw ratio, where meanCVw is the mean value of the CVw for all individuals [95].

**Statistical analyses**

**Behavioural and EEG data:** To investigate how the experimental conditions influenced neural and behavioural responses, we used linear-mixed models (LMMs). Models always included participant ID as a random term. Fixed effects varied depending on the question being addressed and were always first specified as the full model, then interaction terms were dropped if they did not reach significance using the likelihood-ratio test which is suitable in the hypothesis testing framework [96]. For final models, statistical significance of fixed effects was assessed using F-tests and the Kenward–Roger method of degrees-of-freedom approximation, as it has been shown to be a reliable method when LMMs are balanced [97]. Post hoc pairwise comparisons were Tukey corrected. Visual inspection of plots showed that the normality and homoscedasticity of residuals and random effects assumptions were met in all
cases. For full reporting of these models and all other statistical analyses, refer to the S6 Rmd and S7 models tables documents. Unpaired or paired *t* tests, with unequal variance (or their nonparametric equivalent when the normality assumption was not met) were used in comparisons when they were warranted.

**Vocal production data:** We tested for differences in vocal rates between dogs and humans using an LMM with species as a fixed effect and subjects within vocalisation/language types as random effects. We used LMMs to assess whether the acoustic measurements varied with vocalisation/language type adding weight class and subject ID as random terms in dogs, while in humans, the random effects were subject ID and sex. Finally, to assess whether acoustic parameters were allometrically related to body weight in dogs, we first log-transformed the variables and then used an LMM with vocal class and subject ID as random terms. A paired *t* test was used to compare DDS and ADS speech of the same speakers, with the normal distribution and homogeneity of variance assumptions having been met.

All statistical analyses were done in R version 4.2.3 [98] as well as most graphic outputs (except Figs 2, 3D and 4A) and involved the packages: scales [99], lme4 [100], lmerTest [101], Matrix [102], effectsize [103], emmeans [104], MuMln [105], parameters [106], rstatix [107], sjPlot [108], cowplot [109], ggpubr [110], ggsignif [111], gridExtra [112], ggplot2 [113], ggrides [114], and tuneR [115].

**Supporting information**

**S1 Fig. Mean + SE mTRF r values for each electrode, from the human neural data.** The mTRF r values were obtained by first using a decoding model to reconstruct acoustic stimuli from neural data then by correlating this reconstructed acoustic data to the actual stimulus envelope. Thus, higher r values indicate that reconstructed data from that electrode better match the original stimulus. A linear-mixed model using electrodes as fixed effects and human ID as a random term, revealed significant differences among electrodes ($F_{7,67.1} = 3.17, p = 0.006$). Post hoc tests (FDR corrected) were done to compare the mean value of one electrode to the average value of all other electrodes. FCz was the only electrode that showed significantly higher mTRF r values compared to all others (see S4 Data for the corresponding data).

(TIFF)

**S2 Fig. Effect of speech rate on speech neural tracking (cerebro-acoustic coherence).** (A) Slope estimate and 95% CI of syllable rate effect on syllabic coherence in humans for each speech type. (B) Slope estimate and 95% CI of word rate effect on word coherence in humans for each speech type. (C) Slope estimate and 95% CI of syllable rate effect on syllabic coherence in dogs for each speech type. (D) Slope estimate and 95% CI of word rate effect on word coherence in dogs for each speech type. The underlying data can be found in S2 Data.

(TIFF)

**S3 Fig. Comparison of acoustic characteristics of command, DDS and ADS registers.** (A) mean F0 and (B) interquartile range of F0 across the 3 speech types. These analyses are based on 8 participants in each group (2 men and 6 women in all cases). Following a Kruskal-Wallis test, pairwise post hoc significance testing was done using a Dunn test and Bonferroni correction (see S5 Data for data).

(PNG)

**S4 Fig. Time-frequency plots of evoked power for each dog participant.** Data are averaged across all conditions and baselined between -1 s to stimulus onset. Averaged (and SD) evoked
power across all 8 dogs is shown in the last panel (see S8 Data).
(TIFF)

**S5 Fig. Time-frequency plots of evoked power for each human participant.** Data are averaged across all conditions and baselined between -1 s to stimulus onset. Averaged (and SD) evoked power across all 11 humans is shown in the last panel (see S9 Data).
(TIFF)

**S1 Table. Summary statistics of vocal rate (VR) and dominant acoustic frequency (DF) in dog vocalisations.** The potential for individual coding (PIC) is a measure that quantifies the ratio between the inter- and the intra-individual coefficients of variation, with values >1 indicating high individual distinctiveness. Typical vocal contexts are provided, although in some cases (e.g., barks, howls) vocalisations can be used in a range of situations spanning the affiliative-agonistic continuum.
(DOCX)

**S2 Table. Summary statistics of the number of trials remaining after preprocessing EEG data.** Data are shown according to subjects and conditions.
(DOCX)

**S1 Data. Dominant acoustic frequency and speech rate of paired ADS and DDS sentences.**
(XLSX)

**S2 Data. Perception experiment data, including speech rates, cerebro-acoustic coherence, and behavioural variables for humans and dogs.**
(XLSX)

**S3 Data. Delta and theta cerebro-acoustic data for real and permutated acoustic-neural pairings for both humans and dogs.**
(XLSX)

**S4 Data. Results of the mtrf analysis for each electrode and each human subject used to establish which electrode best responded to the auditory stimulation.**
(XLSX)

**S5 Data. Mean and interquartile range pitch data for French speakers of adult-directed, dog-directed, and command speech type.**
(XLSX)

**S6 Data. Rmd file reporting the statistical analyses.**
(PDF)

**S7 Data. Full model results and post hoc analyses for the analyses reported in-text.**
(PDF)

**S8 Data. EEG data (time by frequency) for each dog subject.**
(XLSX)

**S9 Data. EEG data (time by frequency) for each human subject.**
(XLSX)

## Acknowledgments

We are thankful to Silvia Marchesotti and Johanna Nicolle for their help in the initial piloting stage of this work.
## Author Contributions

**Conceptualization:** Eloïse C. Déaux, Florence Gaunet, Thierry Legou, Luc Arnal, Anne-Lise Giraud.

**Data curation:** Eloïse C. Déaux.

**Formal analysis:** Eloïse C. Déaux.

**Funding acquisition:** Anne-Lise Giraud.

**Investigation:** Eloïse C. Déaux, Théophane Piette.

**Methodology:** Eloïse C. Déaux, Luc Arnal.

**Project administration:** Anne-Lise Giraud.

**Resources:** Théophane Piette, Thierry Legou.

**Software:** Eloïse C. Déaux.

**Visualization:** Eloïse C. Déaux.

**Writing – original draft:** Eloïse C. Déaux.

**Writing – review & editing:** Eloïse C. Déaux, Florence Gaunet, Luc Arnal, Anne-Lise Giraud.

## References

1. Charlton BD, Reby D. The evolution of acoustic size exaggeration in terrestrial mammals. Nat Commun. 2016; 7:12739. https://doi.org/10.1038/ncomms12739 PMID: 27598835

2. Taylor AM, Reby D. The contribution of source-filter theory to mammal vocal communication research. J Zool. 2010; 280(3):221-236. https://doi.org/10.1111/j.1469-7998.2009.00661.x

3. Ryan MJ, Fox JH, Wilczynski W, Rand AS. Sexual selection for sensory exploitation in the frog *Physalæmus pustulosus*. Nature. 1990; 343(6253):66-67. https://doi.org/10.1038/343066a0 PMID: 2296291

4. Arnal LH, Flinker A, Kleinschmidt A, Giraud AL, Poeppel D. Human screams occupy a privileged niche in the communication soundscape. Curr Biol. 2015; 25(15):2051-2056. Epub 2015/07/21. https://doi.org/10.1016/j.cub.2015.06.043 PMID: 26190070.

5. Hoy RR, Pollack GS, Moiseff A. Species-Recognition in the Field Cricket, *Teleogryllus oceanicus*: Behavioral and Neural Mechanisms. Am Zool. 2015; 22(3):597-607. https://doi.org/10.1093/icb/22.3.597

6. Ghazanfar AA, Smith-Rohrberg D, Hauser MD. The role of temporal cues in rhesus monkey vocal recognition: Orienting asymmetries to reversed calls. Brain Behav Evol. 2001; 58(3):163-172. https://doi.org/10.1159/000047270 PMID: 11910173

7. Blumstein DT, Armitage KB. Alarm calling in yellow-bellied marmots: I. The meaning of situationally variable alarm calls. Anim Behav. 1997; 53(1):143-171.

8. Drăgănoiu TI, Nagle L, Kreutzer M. Directional female preference for an exaggerated male trait in canary (*Serinus canaria*) song. Proc Biol Sci. 2002; 269(1509):2525-2531. https://doi.org/10.1098/rspb.2002.2192 PMID: 12573066

9. Galeotti P, Sacchi R, Rosa DP, Fasola M. Female preference for fast-rate, high-pitched calls in Hermann's tortoises *Testudo hermanni*. Behav Ecol. 2005; 16(1):301-308. https://doi.org/10.1093/beheco/arh165

10. Pauly GB, Bernal XE, Rand AS, Ryan MJ. The vocal sac increases call rate in the Túngara frog *Physalæmus pustulosus*. Physiol Biochem Zool. 2006; 79(4):708-719. https://doi.org/10.1086/504613 PMID: 16826497

11. Page RA, Bernal XE. The challenge of detecting prey: Private and social information use in predatory bats. Funct Ecol. 2020; 34(2):344-363. https://doi.org/10.1111/1365-2435.13439

12. Magrath RD, Haff TM, Igcic B. Interspecific communication: gaining information from heterospecific alarm calls. Coding strategies in vertebrate acoustic communication. 2020:287-314.
13. Lupo KD. Hounds follow those who feed them: What can the ethnographic record of hunter-gatherers reveal about early human-canid partnerships? J Anthropol Archaeol. 2019; 55:101081. https://doi.org/10.1016/j.jaa.2019.101081

14. Serpell JA. Commensalism or cross-species adoption? A critical review of theories of wolf domestication. Front Vet Sci. 2021; 8:662370. https://doi.org/10.3389/fvets.2021.662370 PMID: 33937382

15. Cohen J, Fox M. Vocalizations in wild canids and possible effects of domestication. Behav Processes. 1976; 1(1):77-92. https://doi.org/10.1016/0376-6357(76)90008-5 PMID: 24923546

16. Pongrácz P. Modeling evolutionary changes in information transfer. European Psychologist. 2017.

17. Burnham D, Kitamura C, Vollmer-Conna U. What's new, pussycat? On talking to babies and animals. Science. 2002; 296(5572):1435. https://doi.org/10.1126/science.1069587 PMID: 12029126

18. Hirsh-Pasek K, Treiman R. Doggerel: Motherese in a new context. J Child Lang. 1982; 9(1):229-237. Epub 2008/09/26. https://doi.org/10.1017/s0305000900003731 PMID: 7061632

19. Ben-Aderet T, Gallego-Abenza M, Reby D, Mathevon N. Dog-directed speech: why do we use it and do dogs pay attention to it? Proc Biol Sci. 1846; 2017(284):20162429. https://doi.org/10.1098/rspb.2016.2429 PMID: 28077769

20. Xu N, Burnham D, Kitamura C, Vollmer-Conna U. Vowel Hyperarticulation in Parrot-, Dog- and Infant-Directed Speech. Anthrozoös. 2013; 26(3):373-380. https://doi.org/10.2752/175303713X13697429463592

21. Attaheri A, Choisdealbha ÁN, Di Liberto GM, Rocha S, Brusini P, Mead N, et al. Delta- and theta-band cortical tracking and phase-amplitude coupling to sung speech by infants. Neuroimage. 2022; 247:118698. https://doi.org/10.1016/j.neuroimage.2021.118698 PMID: 34798233

22. Nallet C, Gervain J. Neurodevelopmental Preparedness for Language in the Neonatal Brain. Annu Rev Dev Psychol. 2021; 3:41-58. https://doi.org/10.1146/annurev-devpsych-050620-025732

23. Hare B, Tomasello M. Human-like social skills in dogs? Trends Cogn Sci. 2005; 9(9):439-444. https://doi.org/10.1016/j.tics.2005.07.003 PMID: 16061417

24. Andics A, Gacsi M, Farago T, Kis A, Miklosi A. Voice-sensitive regions in the dog and human brain are revealed by comparative fMRI. Curr Biol. 2014; 24(5):574-578. Epub 2014/02/25. https://doi.org/10.1016/j.cub.2014.01.058 PMID: 24560578.

25. Ratcliffe VF, Reby D. Orienting asymmetries in dogs' responses to different communicatory components of human speech. Curr Biol. 2014; 24(24):2908-2912. Epub 2014/12/03. https://doi.org/10.1016/j.cub.2014.10.030 PMID: 25454584.

26. Fukuzawa M, Mills DS, Cooper JJ. The effect of human command phonetic characteristics on auditory cognition in dogs (Canis familiaris). J Comp Psychol. 2005; 119(1):117-120. Epub 2005/03/03. https://doi.org/10.1037/0735-7036.119.1.117 PMID: 15740436.

27. Pilley JW, Reid AK. Border collie comprehends object names as verbal referents. Behav Processes. 2011; 86(2):184-195. Epub 2010/12/15. https://doi.org/10.1016/j.beproc.2010.11.007 PMID: 21145379.

28. Kaminski J, Call J, Fischer J. Word learning in a domestic dog: evidence for "fast mapping". Science. 2004; 304(5677):1682-1683.

29. Boros M, Magyari L, Török D, Bozsik A, Deme A, Andics A. Neural processes underlying statistical learning for speech segmentation in dogs. Curr Biol. 2021; 31(24):5512-5521.e5. https://doi.org/10.1016/j.cub.2021.10.017 PMID: 34717832

30. Boros M, Magyari L, Morvai B, Hernández-Pérez R, Dror S, Andics A. Neural evidence for referential understanding of object words in dogs. Curr Biol. https://doi.org/10.1016/j.cub.2024.02.029 PMID: 38521063

31. Molnár C, Pongrácz P, Miklósi Á. Seeing with ears: Sightless humans' perception of dog bark provides a test for structural rules in vocal communication. Q J Exp Psychol. 2010; 63(5):1004-1013. https://doi.org/10.1080/17470210903168243 PMID: 19760535.

32. Scheumann M, Hasting AS, Kotz SA, Zimmermann E. The Voice of Emotion across Species: How Do Human Listeners Recognize Animals' Affective States? PLoS ONE. 2014; 9(3):e91192. https://doi.org/10.1371/journal.pone.0091192 PMID: 24621604

33. Scheumann M, Hasting AS, Zimmermann E, Kotz SA. Human Novelty Response to Emotional Animal Vocalizations: Effects of Phylogeny and Familiarity. Front Behav Neurosci. 2017: 11. https://doi.org/10.3389/fbehav.2017.00204 PMID: 29114210

34. Coupé C, Oh YM, Dediu D, Pellegrino F. Different languages, similar encoding efficiency: Comparable information rates across the human communicative niche. Sci Adv. 2019; 5(9):eaaw2594. https://doi.org/10.1126/sciadv.aaw2594 PMID: 32047854

35. Poeppel D, Assaneo MF. Speech rhythms and their neural foundations. Nat Rev Neurosci. 2020; 21(6):322-334. https://doi.org/10.1038/s41583-020-0304-4 PMID: 32376899
36. MacNeilage PF, Davis BL, Kinney A, Matyear CL. The motor core of speech: A comparison of serial organization patterns in infants and languages. Child Dev. 2000; 71(1):153-163. https://doi.org/10.1111/1467-8624.00129 PMID: 10836569

37. Schroeder CE, Wilson DA, Radman T, Scharfman H, Lakatos P. Dynamics of Active Sensing and perceptual selection. Curr Opin Neurobiol. 2010; 20(2):172-176. https://doi.org/10.1016/j.conb.2010.02.010 PMID: 20307966

38. Hyafil A, Fontolan L, Kabdebon C, Gutkin B, Giraud A-L. Speech encoding by coupled cortical theta and gamma oscillations. Elife. 2015; 4:e06213. https://doi.org/10.7554/eLife.06213 PMID: 26023831

39. Ahissar E, Nagarajan S, Ahissar M, Protopapas A, Mahncke H, Merzenich MM. Speech comprehension is correlated with temporal response patterns recorded from auditory cortex. Proc Natl Acad Sci U S A. 2001; 98(23):13367-13372. https://doi.org/10.1073/pnas.201400998 PMID: 11698688

40. Luo H, Poeppel D. Phase patterns of neuronal responses reliably discriminate speech in human auditory cortex. Neuron. 2007; 54(6):1001-1010. Epub 2007/06/22. https://doi.org/10.1016/j.neuron.2007.06.004 PMID: 17582338.

41. Peelle JE, Gross J, Davis MH. Phase-locked responses to speech in human auditory cortex are enhanced during comprehension. Cereb Cortex. 2013; 23(6):1378-1387. Epub 2012/05/23. https://doi.org/10.1093/cercor/bhs118 PMID: 22610394.

42. Doelling KB, Arnal LH, Ghitza O, Poeppel D. Acoustic landmarks drive delta-theta oscillations to enable speech comprehension by facilitating perceptual parsing. Neuroimage. 2014; 85(Pt 2):761-768. Epub 2013/06/25. https://doi.org/10.1016/j.neuroimage.2013.06.035 PMID: 23791839.

43. Pefkou M, Arnal LH, Fontolan L, Giraud AL. theta-band and beta-band neural activity reflects independent syllable tracking and comprehension of time-compressed speech. J Neurosci. 2017; 37 (33):7930-7938. Epub 2017/07/22. https://doi.org/10.1523/JNEUROSCI.2882-16.2017 PMID: 28729443.

44. Doelling KB, Arnal LH, Assaneo MF. Adaptive oscillators provide a hard-coded Bayesian mechanism for rhythmic inference. bioRxiv. 2022.02:2022.06.18.496664. https://doi.org/10.1101/2022.06.18.496664

45. Giraud AL, Poeppel D. Cortical oscillations and speech processing: emerging computational principles and operations. Nat Neurosci. 2012; 15(4):511-517. Epub 2012/03/20. https://doi.org/10.1038/nn.3063 PMID: 22426255.

46. Fitch T. The phonetic potential of nonhuman vocal tracts: Comparative cineradiographic observations of vocalizing animals. Phonetica. 2000; 57(2-4):205-218. https://doi.org/10.1159/000028474 PMID: 10992141

47. Lieberman P. Vocal tract anatomy and the neural bases of talking. J Phon. 2012; 40(4):608-622.

48. Boë L-J, Sawallis TR, Fagot J, Badin P, Barbier G, Captier G, et al. Which way to the dawn of speech?: Reanalyzing half a century of debates and data in light of speech science. Sci Adv. 2019; 5(12): eaaw3916. https://doi.org/10.1126/sciadv.aaw3916 PMID: 32076631

49. Andics A, Gabor A, Gacsi M, Farago T, Szabo D, Miklosi A. Neural mechanisms for lexical processing in dogs. Science. 2016; 353:1030-1032. https://doi.org/10.1126/science.aaf3777 PMID: 27576923

50. Bálint A, Szabó Á, Andics A, Gácsi M. Dog and human neural sensitivity to voicelikeness: A comparative fMRI study. Neuroimage. 2023; 265:119791. https://doi.org/10.1016/j.neuroimage.2022.119791 PMID: 36476565

51. Benjamin A, Slocombe K. 'Who's a good boy?!" Dogs prefer naturalistic dog-directed speech. Anim Cogn. 2018; 21(3):353-364. https://doi.org/10.1007/s10071-018-1172-4 PMID: 29500713

52. Yin S, McCowan B. Barking in domestic dogs: context specificity and individual identification. Anim Behav. 2004; 68(2):343-355.

53. Molnár C, Pongrácz P, Faragó T, Dóka A, Miklósi Á. Dogs discriminate between barks: The effect of context and identity of the caller. Behav Processes. 2009; 82(2):198-201. https://doi.org/10.1016/j.beproc.2009.06.011 PMID: 19596426

54. Bowling D, Garcia M, Dunn J, Ruprecht R, Stewart A, Frommolt K-H, et al. Body size and vocalization in primates and carnivores. Sci Rep. 2017; 7:41070. https://doi.org/10.1038/srep41070 PMID: 28117380

55. Howell TJ, Conduit R, Toukhsati S, Bennett P. Auditory stimulus discrimination recorded in dogs, as indicated by mismatch negativity (MMN). Behav Processes. 2012; 89(1):8-13. Epub 2011/10/18. https://doi.org/10.1016/j.beproc.2011.09.009 PMID: 22001730.

56. Magyari L, Huszár Z, Turzó A, Andics A. Event-related potentials reveal limited readiness to access phonetic details during word processing in dogs. R Soc Open Sci. 2020; 7(12):200851. https://doi.org/10.1098/rsos.200851 PMID: 33489257
57. Gábor A, Gácsi M, Szabó D, Miklósi Á, Kubinyi E, Andics A. Multilevel fMRI adaptation for spoken word processing in the awake dog brain. Sci Rep. 2020; 10(1):11968. https://doi.org/10.1038/s41598-020-68821-6 PMID: 32747731

58. Lakatos P, Gross J, Thut G. A new unifying account of the roles of neuronal entrainment. Curr Biol. 2019; 29(18):R890–R905. https://doi.org/10.1016/j.cub.2019.07.075 PMID: 31550478

59. Boari S, Mindlin GB, Amador A. Neural oscillations are locked to birdsong rhythms in canaries. Eur J Neurosci. 2022; 55(2):549–565. https://doi.org/10.1111/ejn.15552 PMID: 34852183

60. Theunissen FE, Shaevitz SS. Auditory processing of vocal sounds in birds. Curr Opin Neurobiol. 2006; 16(4):400–407. https://doi.org/10.1016/j.conb.2006.07.003 PMID: 16842993

61. Cuaya LV, Hernández-Pérez R, Boros M, Deme A, Andics A. Speech naturalness detection and language representation in the dog brain. Neuroimage. 2022; 248:118811. https://doi.org/10.1016/j.neuroimage.2021.118811 PMID: 34906714

62. Pellegrino F, Coupé C, Marsico E. A cross-language perspective on speech information rate. Language. 2011:539–558.

63. Jouventin P, Aubin T. Acoustic systems are adapted to breeding ecologies: individual recognition in nesting penguins. Anim Behav. 2002; 64(5):747–757. https://doi.org/10.1006/anbe.2002.4002

64. Taylor AM, Reby D, McComb K. Context-related variation in the vocal growling behaviour of the domestic dog (Canis familiaris). Ethology. 2009; 115(10):905–915. https://doi.org/10.1111/j.1439-0310.2009.01681.x

65. Morrill RJ, Paukner A, Ferrari PF, Ghazanfar AA. Monkey lipsmacking develops like the human speech rhythm. Dev Sci. 2012; 15(4):557–568. https://doi.org/10.1111/j.1467-7687.2012.01149.x PMID: 22709404

66. Risueno-Segovia C, Hage SR. Theta Synchronization of Phonatory and Articulatory Systems in Mar-moset Monkey Vocal Production. Curr Biol. 2020. https://doi.org/10.1016/j.cub.2020.08.019 PMID: 32888481

67. Pereira AS, Kavanagh E, Hobaiter C, Slocombe KE, Lameira AR. Chimpanzee lip-smacks confirm pri-mate continuity for speech-rhythm evolution. Biol Lett. 2020; 16(5):20200232. https://doi.org/10.1098/rsbl.2020.0232 PMID: 32453963

68. MacNeilage PF. The frame/content theory of evolution of speech production. Behav Brain Sci. 1998; 21(4):499–511 Epub 1998/08/01. https://doi.org/10.1017/s0140525x98001265 PMID: 10097020

69. Buzsáki G, Logothetis N, Singer W. Scaling Brain Size, Keeping Timing: Evolutionary Preservation of Brain Rhythms. Neuron. 2013; 80(3):751–764. https://doi.org/10.1016/j.neuron.2013.10.002 PMID: 24183025

70. Doelling KB, Poeppel D. Cortical entrainment to music and its modulation by expertise. Proc Natl Acad Sci U S A. 2015; 112(45):E6233–E6242. Epub 2015/10/28. https://doi.org/10.1073/pnas.1508431112 PMID: 26504238.

71. Meyer L. The neural oscillations of speech processing and language comprehension: state of the art and emerging mechanisms. Eur J Neurosci. 2018; 48(7):2609–2621. https://doi.org/10.1111/ejn. 13748 PMID: 29055058

72. Ghitza O. The theta-syllable: a unit of speech information defined by cortical function. Front Psychol. 2013; 4. https://doi.org/10.3389/fpsyg.2013.00138 PMID: 23519170

73. Ding N, Melloni L, Zhang H, Tian X, Poeppel D. Cortical tracking of hierarchical linguistic structures in connected speech. Nat Neurosci. 2015; 19:158. Available from: https://www.nature.com/articles/nn. 4186#supplementary-information. PMID: 26642090

74. De Serrano AR, Weadick CJ, Price AC, Rodd FH. Seeing orange: prawns tap into a pre-existing sen-sory bias of the Trinidadian guppy. Proc Biol Sci. 2012; 279(1741):3321–3328. https://doi.org/10.1098/rspb.2012.0633 PMID: 22593111

75. Owren MJ, Rendall D. Sound on the rebound: Bringing form and function back to the forefront in under-standing nonhuman primate vocal signaling. Evol Anthropol. 2001; 10(2):58–71. https://doi.org/10. 1002/evan.1014

76. Guilford T, Dawkins MS. Receiver psychology and the evolution of animal signals. Anim Behav. 1991; 42(1):1–14. https://doi.org/10.1016/S0003-3472(05)80600-1

77. Kalashnikova M, Peter V, Di Liberto GM, Lalor EC, Burnham D. Infant-directed speech facilitates seven-month-old infants’ cortical tracking of speech. Sci Rep. 2018; 8(1):13745. https://doi.org/10. 1038/s41598-018-32150-6 PMID: 30214000

78. Leong V, Kalashnikova M, Burnham D, Goswami U. The temporal modulation structure of infant-directed speech. Open Mind. 2017; 1(2):78–90. https://doi.org/10.1162/OPMI_a_00008
79. Ghitza O. Linking Speech Perception and Neurophysiology: Speech Decoding Guided by Cascaded Oscillators Locked to the Input Rhythm. Front Psychol. 2011; 2(130). https://doi.org/10.3389/fpsyg.2011.00130 PMID: 21743809

80. Hovsepyan S, Olasagasti I, Giraud A-L. Combining predictive coding and neural oscillations enables online syllable recognition in natural speech. Nat Commun. 2020; 11(1):3117. https://doi.org/10.1038/s41467-020-16956-5 PMID: 32561726

81. Su Y, MacGregor LJ, Olasagasti I, Giraud A-L. A deep hierarchy of predictions enables online meaning extraction in a computational model of human speech comprehension. PLoS Biol. 2023; 21(3): e3002046. https://doi.org/10.1371/journal.pbio.3002046 PMID: 36947552

82. Lizarazu M, Carreiras M, Molinaro N. Theta-gamma phase-amplitude coupling in auditory cortex is modulated by language proficiency. Hum Brain Mapp. 2023; 44(7):2862-2872. https://doi.org/10.1002/hbm.26250 PMID: 36852454

83. Cheung C, Hamilton LS, Johnson K, Chang EF. The auditory representation of speech sounds in human motor cortex. Elife. 2016; 5:e12577. https://doi.org/10.7554/eLife.12577 PMID: 26943778

84. Meister IG, Wilson SM, Deblieck C, Wu AD, Iacoboni M. The Essential Role of Premotor Cortex in Speech Perception. Curr Biol. 2007; 17(19):1692-1696. https://doi.org/10.1016/j.cub.2007.08.064 PMID: 17900904

85. Williams H, Nottebohm F. Auditory responses in avian vocal motor neurons: a motor theory for song perception in birds. Science. 1985; 229(4710):279-282. https://doi.org/10.1126/science.4012321 PMID: 4012321

86. Márton CD, Fukushima M, Camalier CR, Schultz SR, Averbeck BB. Signature Patterns for Top-Down and Bottom-Up Information Processing via Cross-Frequency Coupling in Macaque Auditory Cortex. eneuro. 2019; 6(2):ENEURO.0467-18.2019. https://doi.org/10.1523/ENEURO.0467-18.2019 PMID: 31088914

87. Archakov D, DeWitt I, Kuśmierek P, Ortiz-Rios M, Cameron D, Cui D, et al. Auditory representation of learned sound sequences in motor regions of the macaque brain. Proc Natl Acad Sci. 2020; 117 (26):15242-15252. https://doi.org/10.1073/pnas.1915610117 PMID: 32541016

88. Kikuchi Y, Attaheri A, Wilson B, Rhone AE, Nourski KV, Gander PE, et al. Sequence learning modulates neural responses and oscillatory coupling in human and monkey auditory cortex. PLoS Biol. 2017; 15(4):e2000219. https://doi.org/10.1371/journal.pbio.2000219 PMID: 28441393

89. Pellegrino F. Canine electroencephalographic recording technique: findings in normal and epileptic dogs. Clin Neurophysiol. 2004; 115(2):477-487. https://doi.org/10.1016/S1388-2457(03)00347-X PMID: 14744591

90. Gemmeke JF, Ellis DP, Freedman D, Jansen A, Lawrence W, Moore RC, et al., editors. Audio set: An ontology and human-labeled dataset for audio events. 2017 IEEE international conference on acoustics, speech and signal processing (ICASSP). 2017: IEEE.

91. Salt C, Morris PJ, German AJ, Wilson D, Lund EM, Cole TJ, et al. Growth standard charts for monitoring bodyweight in dogs of different sizes. PLoS ONE. 2017; 12(9):e0182064. https://doi.org/10.1371/journal.pone.0182064 PMID: 28873413

92. Oostenveld R, Fries P, Maris E, Schoffelen J-M. FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput Intell Neurosci. 2011; 2011:1-9.

93. Crosse MJ, Di Liberto GM, Bednar A, Lalor EC. The Multivariate Temporal Response Function (mTRF) Toolbox: A MATLAB Toolbox for Relating Neural Signals to Continuous Stimuli. Front Hum Neurosci. 2016; 10. https://doi.org/10.3389/fnhum.2016.00604 PMID: 27965557

94. Sueur J, Aubin T, Simonis C. Seewave: a free modular tool for sound analysis and synthesis. Bioacoustics. 2008; 18:213-226.

95. Robisson P, Aubin T, Bremond J-C. Individuality in the voice of the emperor penguin Aptenodytes forsteri: adaptation to a noisy environment. Ethology. 1993; 94(4):279-290.

96. Lewis F, Butler A, Gilbert L. A unified approach to model selection using the likelihood ratio test. Methods Ecol Evol. 2011; 2(2):155-162. https://doi.org/10.1111/j.2041-210X.2010.00063.x

97. Schaalje GB, McBride JB, Fellingham GW. Adequacy of approximations to distributions of test statistics in complex mixed linear models. J Agric Biol Environ Stat. 2002; 7(4):512-524. https://doi.org/10.1198/108571102726

98. R Core T. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/. 2014.

99. Wickham H, Wickham MH. RColorBrewer I. Package 'scales'. R package version. 2016;1(0).

100. Bates D, Maechler M, Bolker B, Walker S, Christensen RHB, Singmann H, et al. Package 'lme4'. Convergence. 2015;12(1):2.
101. Kuznetsova A, Brockhoff PB, Christensen RH. ImerTest package: tests in linear mixed effects models. J Stat Softw. 2017; 82(13):1–26.

102. Bates D, Maechler M. Matrix: Sparse and dense matrix classes and methods. R package version 0999375–43. http://cran.r-project.org/package=Matrix. 2010.

103. Ben-Shachar MS, Lüdecke D, Makowski D. effectsize: Estimation of effect size indices and standardized parameters. J Open Source Softw. 2020; 5(56):2815.

104. Lenth R, Singmann H, Love J, Buerkner P, Herve M. emmeans: estimated marginal means. R package version. 2020;1(4):4.

105. Barton K. MuMln: Multi-model inference. R package version 1.9.5. R Project for Statistical Computing, Vienna, Austria. 2013.

106. Lüdecke D, Ben-Shachar MS, Patil I, Makowski D. Extracting, computing and exploring the parameters of statistical models using R. J Open Source Softw. 2020; 5(53):2445.

107. Kassambara A. Pipe-friendly framework for basic statistical tests [R Package Rstatix Version 0.7.0]. Free Software Foundation Inc: Boston, MA, USA. 2021.

108. Lüdecke D. sjPlot: Data visualization for statistics in social science. R package version. 2018;2(1).

109. Wilke CO. cowplot: streamlined plot theme and plot annotations for 'ggplot2'. R package version. 2019;1(0).

110. Kassambara A. ggpubr: 'ggplot2' based publication ready plots. R package version. 2018:2.

111. Ahlmann-Eltze C. ggsignif: Significance brackets for 'ggplot2'. R package version 0.5.0. 2019.

112. Auguie B, Antonov A. gridExtra: miscellaneous functions for "grid" graphics. R package version. 2017;2(601):602.

113. Wickham H. ggplot2: elegant graphics for data analysis: Springer; 2016.

114. Wilke CO. Ggrides: Ridgeline plots in 'ggplot2'. R package version 05. 2018;1:483.

115. Ligges U, Preusser A, Thieler A, Mielke J, Weihs C. Package 'tuneR'. Version; 2018.