Vocal variation of Ferruginous ( ferruginea) is related to geographic distances

Guilherme S Cardoso, Reginaldo J Donatelli  Correspondence Laboratory of Ornithology, School of Sciences/São Paulo State University [email protected] [email protected] Abstract

 Disciplines The vocal variation in that produce innate vocalizations reflects the differences Ecology that occur within the species’ gene pool. These variations can indicate whether the

 Keywords species is undergoing a process of isolation by distance. In this case, the vocalization is Bioacoustics an auxiliary tool for understanding the action of evolutionary processes in birds. Thus, Ornithology we analyzed vocalizations of (a suboscine that has innate vo- Zoology calizations) and correlated the vocal features to the geographical distances among indi-

 Type of Observation viduals. We collected vocalizations from sound archives and analyzed frequency Standalone and time parameters of each sound file. The variables were summarized by a principal

 Type of Link component analysis (PCA) and then we compared the distances among these compo- Standard Data nents with the geographical distances among samples using Mantel’s test. Thus, we have found evidence that vocal differences are related to geographic separation: the larger  Submitted Mar 23, 2016 the distance between individuals, the greater the difference in vocalizations. These dis- ⴑ Published Jun 21, 2016 similarities occur both in time and in frequency of vocalizations, showing that processes of geographic isolation are occurring within the species.

Introduction 3 x Song variation has been mainly studied in the Oscine [1], a group of birds in which learning and auditory feedback exert significant influence on song development [2]. However, there is evidence that suboscine birds, such as the Ferruginous Antbird Triple Blind Peer Review (Drymophila ferruginea), develop their vocalizations normally in the absence of vocal The handling editor, the re- tutors [3] [4]. This fact could suggest that genetic components exert more influence on viewers, and the authors are the vocal development than auditory feedback and extensive learning [5] [6]. all blinded during the review process. If the vocal differences between distant populations of suboscine species were related to their geographical distances, this could imply that these populations are undergoing iso- lation by distance. This can be an indication of allopatric divergence within the species, which is an initial step to the process of speciation [7]. Such vocal variation may be used as an auxiliary tool in taxonomic studies between species and even characterize different subspecies [8].

Full Open Access Supported by the Velux Objective Foundation, the University of We aim to determine if the vocalizations of the Ferruginous have variations Zurich, and the EPFL School that are correlated to the geographical distances among individuals. of Life Sciences.

4.0

Creative Commons 4.0 This observation is dis- tributed under the terms of the Creative Commons Attribution 4.0 International License.

a

DOI: 10.19185/matters.201605000015 Matters (ISSN: 2297-8240) 1

| Vocal variation of Ferruginous Antbird (Drymophila ferruginea) is related to geographic distances

b

DOI: 10.19185/matters.201605000015 Matters (ISSN: 2297-8240) 2

| Vocal variation of Ferruginous Antbird (Drymophila ferruginea) is related to geographic distances

c

DOI: 10.19185/matters.201605000015 Matters (ISSN: 2297-8240) 3

| Vocal variation of Ferruginous Antbird (Drymophila ferruginea) is related to geographic distances

d

Figure Legend Figure 1. (A) Male individuals of Ferruginous Antbird (photo taken by Rafael M. Martins, 2015). (B) Distribuition of Ferruginous Antbird. The red circles represent individuals whose vocalizations were collected. The numbers above and below them represent the number of individuals whose vocalizations were recorded in the same place. (C) Visual representations of the vocalizations. (1) The spectrogram displays the two syllables of the vocalization. (2) The oscillogram shows measurement of the temporal variables. (3) The power spectra of each section (i-iv, as shown in the oscillogram) shows the measurement of spectral variables. (D) Principal components (PC) scores for each variable. The eigenvalue and the per- centage of variance summarized by the PCs are assigned below each component. Scores values greater than|0.5|were highlighted in gray. Study Species The Ferruginous Antbird (Drymophyla ferruginea Temminck, 1822) is a species dis- tributed in areas of dense Evergreen Forest of the Atlantic domain (Fig. 1A). It is found only in the Southeast of Brazil, mostly near coastal areas, but also extending inland from Minas Gerais and Sao Paulo (Fig. 1B). It inhabits the understory and mid-story of thick- ets in primary forests and also in second-growth woodland. It is an active forager that feeds on arthropods. This species is not globally threatened, and it is fairly common throughout its range [18]. File Survey Vocalization recordings were obtained from Laboratory of Ornithology (School of Sci- ences, São Paulo State University), Xeno-Canto (www.xeno-canto.org), and Macauley Library (Cornell Lab of Ornithology, Cornell University). Each sound file represented a different individual. Files with excessive noise, poor recording quality, overlapping vocalizations, or any other factor that could hinder audio and spectrogram inspection were excluded from the sample. We sought to include recordings that cover as much as possible the species’ geographic range that was proposed by [18]. The coordinates and the year of recording were provided by the recordists of each vocalization. A total of 45 sound archives were sampled.

DOI: 10.19185/matters.201605000015 Matters (ISSN: 2297-8240) 4

| Vocal variation of Ferruginous Antbird (Drymophila ferruginea) is related to geographic distances

Vocalization Analysis All the vocalizations were measured blind both to the identity and the location of the vocalization, to prevent the measurements from being biased in any way. The identity of the recording was extracted after the measurements. The digitized sound files were standardized (44.100 Hz, 16 bits of resolution, WAV format) using Sony SoundForge Pro 11 (Sony Creative Software, 2010). We randomly selected 3 vocalizations from each sound file, totaling 135 analyzed vocalizations. The vocalizations were selected accord- ing to the morphology of their syllables, which were visible by the spectrogram (Han- ning Window, FFT = 512 points, Overlap = 93,75%) generated in the Sonic Visualizer software (Queen Mary, University of London, 2013). We employed the oscillogram to measure the time variables and the amplitude spectra (Hanning Window, resolution = 2048 points) to measure the frequencies (Fig. 1C and D). A limit of -42 db below the peak amplitude was taken on the power spectra to avoid the noise frequency bands. They were used to measure these variables because spectrograms have a trade-off be- tween temporal and frequency resolution, which may lead to major problems on mea- surement [19]. Both were generated using Sony SoundForge Pro 11 software and the variables were measured by the same person. The following variables were obtained: T1- duration of the first syllable(s); T2-dura- tion of the interval between the two syllables(s); T3- duration of the first frequency modulation (FM) of the second syllable(s); T4- duration of the second FM of the second syllable(s); T5- duration of the second syllable(s); T6- duration of the entire vocaliza- tion; F1- Maximum frequency of the first syllable (kHz); F2- Minimum frequency of the first syllable (kHz); F3- Frequency bandwidth (MAX-MIN) of the first syllable (kHz); F4- Peak frequency of the first syllable (kHz); F5- Maximum frequency of the second syllable (kHz); F6- Minimum frequency of the second syllable (kHz); F7- Frequency bandwidth (MAX-MIN) of the second syllable (kHz); F8- Peak frequency of the first FM of the sec- ond syllable (kHz); F9- Peak frequency of the second FM of the second syllable (kHz); F10- Difference between the peak frequencies of the second and the first FM of the sec- ond syllable (F9 minus F8) (kHz). Statistical Analysis We calculated the mean of the variables using three vocalizations for each sound file. We performed a principal component analysis (PCA) on the mean of variables to reduce its number, to avoid the colinearity of some variables and to decrease the data redun- dancy. We selected the components that had eigenvalues greater than|1.0|and extracted the scores of the principal component space for each individual. We created a “vocal distance” matrix among the individuals using the Euclidean distance between each pair of individuals in the principal component space. A matrix of “geographical distances” among individuals was created by Geographic Distance Matrix Generator 1.2.3 software (Center for Biodiversity and Conservation, American Museum of Natural History) em- ploying the coordinates of latitude and longitude that were obtained from the sound files. We used the Mantel test to examine if the both matrices are correlated to each other. The factor “year of recording” was analyzed by a second Mantel’s test. Thistest demonstrated that there was no effect of the year of recording on the characteristics of vocalizations. We also tested if the differences among the years of recording affects the vocal distances with another Mantel test. The “temporal distance” matrix were created using the Euclidean distance between the years of recording of each pair of individuals. Results & Discussion

The PCA generated 6 components with eigenvalues greater than|1.0|(Fig. 1D). These components explained 90.19% of the variance among the samples on the 16 acoustic variables. PC1 had high positive loadings for frequency features (F1, F2, F4, F5, F6, F8, and F9). PC2 had positive loadings for three time features (T4, T5 and T6). PC3 had positive loadings for T3 and T5, and negative loadings to F3 and F7. PC4 had negative loadings for time features (T2 and T3). PC5 had a positive loading for F10. PC6 showed no predominance among loading values. The component analysis showed the relevance

DOI: 10.19185/matters.201605000015 Matters (ISSN: 2297-8240) 5

| Vocal variation of Ferruginous Antbird (Drymophila ferruginea) is related to geographic distances

of each variable. PC1 was strongly related to spectral features while PC2 was related to temporal ones. The two components together summarized almost 50% of the variation, showing that both time and frequency features varied among samples and that they are equally important to define the vocal variations within the species. The Mantel test showed a significant and positive relationship between geographic and vocal distances (N = 10000 permutations; r = 0.4502; p = 0.0003), but it did not find any relationship between the vocal distances and differences among the years of recording (N = 10000 permutations; r = -0.0658; p = 0.6949). Several studies demonstrate vocal variation across geographic distances. However, they are focused on variation in the oscine passerines [9] [10] [11] [12] [13]. The papers that have looked for variations in vocalizations on non-oscine birds are mainly restricted to comparing between shorter distances (less than 50 km) [14] [15] [16]. Our work is one of the few studies demonstrating geographical variation on suboscine vocalizations over large distances (average distance between pairs of individuals = 483 km) (e.g. [6] on Alder Flycatcher (Empidonax alnorum) and [17] on Common Cuckoo (Cuculus canorus)) and the first one to explore such effect in the Thamnophilidae family. These results suggest that vocal differences and geographic distances are closely cor- related. This could be explained by the “isolation by distance” model. Differences in vocalizations of Thamnophilidae species are mainly affected by genetic factors, andex- tensive learning plays no significant role in vocal ontogeny [4]. Therefore, the corre- lation between distance and vocal differences may be a reflection of variations on the genetic structure originated by isolation of Ferruginous Antbird populations.

Conclusions The vocal variation on Ferruginous Antbird is related to the geographic distance be- tween individuals. This fact may indicate that there is a certain degree of genetic isola- tion between populations of this species, which may in turn be a reflection of geographic distance.

Limitations Our study does not cover molecular analysis, despite making inferences about this ge- netic aspects based on references. Such studies in this species need to be done to inves- tigate if the vocal variation is indeed a reflection of genetic variation.

Additional Information

Methods Study Species The Ferruginous Antbird (Drymophyla ferruginea Temminck, 1822) is a species dis- tributed in areas of dense Evergreen Forest of the Atlantic domain (Fig. 1A). It is found only in the Southeast of Brazil, mostly near coastal areas, but also extending inland from Minas Gerais and Sao Paulo (Fig. 1B). It inhabits the understory and mid-story of thick- ets in primary forests and also in second-growth woodland. It is an active forager that feeds on arthropods. This species is not globally threatened, and it is fairly common throughout its range [18]. File Survey Vocalization recordings were obtained from Laboratory of Ornithology (School of Sci- ences, São Paulo State University), Xeno-Canto (www.xeno-canto.org), and Macauley Library (Cornell Lab of Ornithology, Cornell University). Each sound file represented a different individual. Files with excessive noise, poor recording quality, overlapping vocalizations, or any other factor that could hinder audio and spectrogram inspection were excluded from the sample. We sought to include recordings that cover as much as possible the species’ geographic range that was proposed by [18]. The coordinates and the year of recording were provided by the recordists of each vocalization. A total of 45 sound archives were sampled. Vocalization Analysis

DOI: 10.19185/matters.201605000015 Matters (ISSN: 2297-8240) 6

| Vocal variation of Ferruginous Antbird (Drymophila ferruginea) is related to geographic distances

All the vocalizations were measured blind both to the identity and the location of the vocalization, to prevent the measurements from being biased in any way. The identity of the recording was extracted after the measurements. The digitized sound files were standardized (44.100 Hz, 16 bits of resolution, WAV format) using Sony SoundForge Pro 11 (Sony Creative Software, 2010). We randomly selected 3 vocalizations from each sound file, totaling 135 analyzed vocalizations. The vocalizations were selected accord- ing to the morphology of their syllables, which were visible by the spectrogram (Han- ning Window, FFT = 512 points, Overlap = 93,75%) generated in the Sonic Visualizer software (Queen Mary, University of London, 2013). We employed the oscillogram to measure the time variables and the amplitude spectra (Hanning Window, resolution = 2048 points) to measure the frequencies (Fig. 1C and D). A limit of -42 db below the peak amplitude was taken on the power spectra to avoid the noise frequency bands. They were used to measure these variables because spectrograms have a trade-off be- tween temporal and frequency resolution, which may lead to major problems on mea- surement [19]. Both were generated using Sony SoundForge Pro 11 software and the variables were measured by the same person. The following variables were obtained: T1- duration of the first syllable(s); T2-dura- tion of the interval between the two syllables(s); T3- duration of the first frequency modulation (FM) of the second syllable(s); T4- duration of the second FM of the second syllable(s); T5- duration of the second syllable(s); T6- duration of the entire vocaliza- tion; F1- Maximum frequency of the first syllable (kHz); F2- Minimum frequency of the first syllable (kHz); F3- Frequency bandwidth (MAX-MIN) of the first syllable (kHz); F4- Peak frequency of the first syllable (kHz); F5- Maximum frequency of the second syllable (kHz); F6- Minimum frequency of the second syllable (kHz); F7- Frequency bandwidth (MAX-MIN) of the second syllable (kHz); F8- Peak frequency of the first FM of the sec- ond syllable (kHz); F9- Peak frequency of the second FM of the second syllable (kHz); F10- Difference between the peak frequencies of the second and the first FM of the sec- ond syllable (F9 minus F8) (kHz). Statistical Analysis We calculated the mean of the variables using three vocalizations for each sound file. We performed a principal component analysis (PCA) on the mean of variables to reduce its number, to avoid the colinearity of some variables and to decrease the data redun- dancy. We selected the components that had eigenvalues greater than|1.0|and extracted the scores of the principal component space for each individual. We created a “vocal distance” matrix among the individuals using the Euclidean distance between each pair of individuals in the principal component space. A matrix of “geographical distances” among individuals was created by Geographic Distance Matrix Generator 1.2.3 software (Center for Biodiversity and Conservation, American Museum of Natural History) em- ploying the coordinates of latitude and longitude that were obtained from the sound files. We used the Mantel test to examine if the both matrices are correlated to each other. The factor “year of recording” was analyzed by a second Mantel’s test. Thistest demonstrated that there was no effect of the year of recording on the characteristics of vocalizations. We also tested if the differences among the years of recording affects the vocal distances with another Mantel test. The “temporal distance” matrix were created using the Euclidean distance between the years of recording of each pair of individuals.

Supplementary Material Please see https://sciencematters.io/articles/201605000015.

Funding Statement This work was supported by National Council of Technological and Scientific Develop- ment (CNPq).

Acknowledgements We are very thankful to Rafael M. Martins for the photo, and to all the recordists that registered the vocalizations used in this work.

DOI: 10.19185/matters.201605000015 Matters (ISSN: 2297-8240) 7

| Ethics Statement Not applicable.

Citations

[1] Jeffrey Podos and Paige S. Warren and. “The Evolution of [10] Todd T. Tracy and Myron C. Baker. “Geographic Variation in Geographic Variation in Birdsong”. In: Advances in the Study of Syllables of House Finch Songs”. In: The Auk 116.3 (July 1999), Behavior 37 (2007), pp. 403–458. doi: pp. 666–676. doi: 10.2307/4089328. url: 10.1016/s0065-3454(07)37009-5. url: http://dx.doi.org/10.2307/4089328. http://dx.doi.org/10.1016/S0065- [11] Rachel T. Bolus. “Geographic variation in songs of the Common 3454(07)37009-5. Yellowthroat”. In: The Auk 131.2 (Apr. 2014), pp. 175–185. doi: [2] Kroodsma and D. E. “Learning and the Ontogeny of Sound 10.1642/auk-12-187.1. url: Signals in Birds”. In: Acoustic Communication in Birds and http://dx.doi.org/10.1642/AUK-12-187.1. Academic Press and New York 2 (1982), pp. 1–23. [12] Kim G Mortega, Heiner Flinks, and Barbara Helm. “Behavioural [3] Donald E. Kroodsma. “Songs of the Alder Flycatcher (Empidonax response of a migratory songbird to geographic variation in song alnorum) and Willow Flycatcher (Empidonax traillii) Are Innate”. and morphology”. In: Frontiers in Zoology 11.1 (Nov. 2014), p. 85. In: The Auk 101 (1984), pp. 13–24. doi: 10.1186/s12983-014-0085-6. url: http: [4] Janeene M. Touchton, Nathalie Seddon, and Joseph A. Tobias. //dx.doi.org/10.1186/s12983-014-0085-6. “Captive Rearing Experiments Confirm Song Development [13] Scott M. Ramsay and Ken A. Otter. “Geographic variation in without Learning in a Tracheophone Suboscine Bird”. In: PLOS White-throated Sparrow song may arise through cultural drift”. ONE 9.4 (Apr. 2014), e95746. doi: In: Journal of Ornithology 156.3 (Mar. 2015), pp. 763–773. doi: 10.1371/journal.pone.0095746. url: http: 10.1007/s10336-015-1183-8. url: http: //dx.doi.org/10.1371/journal.pone.0095746. //dx.doi.org/10.1007/s10336-015-1183-8. [5] Luis F. Baptista and Sandra L. L. Gaunt. “Advances in Studies of Avian Sound Communication”. In: The Condor 96.3 (Aug. 1994), [14] Robert B. Payne and Paul Budde. “Song Differences and Map pp. 817–830. doi: 10.2307/1369489. url: Distances in a Population of Acadian Flycatchers”. In: The Wilson http://dx.doi.org/10.2307/1369489. Bulletin 91 (1979), pp. 29–41. [6] Scott F. Lovell and M. Ross Lein. “Geographical Variation In [15] Scott F. Lovell and M. Ross Lein. “Song variation in a population Songs of A Suboscine , the Alder Flycatcher ( of Alder Flycatchers”. In: Journal of Field Ornithology 75.2 (Apr. Empidonax alnorum )”. In: The Wilson Journal of Ornithology 2004), pp. 146–151. doi: 125.1 (Mar. 2013), pp. 15–23. doi: 10.1676/12-087.1. url: 10.1648/0273-8570-75.2.146. url: http: http://dx.doi.org/10.1676/12-087.1. //dx.doi.org/10.1648/0273-8570-75.2.146. [7] NATHALIE SEDDON and JOSEPH A. TOBIAS. “Song divergence [16] M. Ross Lein. “Song Variation In Buff-Breasted Flycatchers at the edge of Amazonia: an empirical test of the peripatric (Empidonax fulvifrons)”. In: The Wilson Journal of Ornithology speciation model”. In: Biological Journal of the Linnean Society 120.2 (June 2008), pp. 256–267. doi: 10.1676/07-067.1. 90.1 (Jan. 2007), pp. 173–188. doi: url: http://dx.doi.org/10.1676/07-067.1. 10.1111/j.1095-8312.2007.00753.x. url: http://dx.doi.org/10.1111/j.1095- [17] Chentao Wei et al. “Geographic variation in the calls of the 8312.2007.00753.x. Common Cuckoo (Cuculus canorus): isolation by distance and divergence among subspecies”. In: Journal of Ornithology 156.2 [8] James A. Sedgwick. “GEOGRAPHIC VARIATION IN THE SONG (Jan. 2015), pp. 533–542. doi: OF WILLOW FLYCATCHERS: DIFFERENTIATION BETWEEN 10.1007/s10336-014-1153-6. url: http: EMPIDONAX TRAILLII ADASTUS AND E. T. EXTIMUS”. In: The //dx.doi.org/10.1007/s10336-014-1153-6. Auk 118.2 (2001), pp. 366–379. doi: 10.1642/0004- 8038(2001)118[0366:gvitso]2.0.co;2. url: [18] K. J. Zimmer and M. L. Isler. “Family Thamnophilidae (Typical http://dx.doi.org/10.1642/0004- Antbirds)”. In: Handbook of the Birds of the World: Broadbills to 8038(2001)118[0366:GVITSO]2.0.CO;2. Tapaculos and Lynx Editions 8 (2003), pp. 448–681. [9] Michael A. Cunningham, Myron C. Baker, and [19] Sue Anne Zollinger et al. “On the relationship between, and Thomas J. Boardman. “Microgeographic Song Variation inthe measurement of, amplitude and frequency in birdsong”. In: Nuttall’s White-Crowned Sparrow”. In: The Condor 89.2 (May Animal Behaviour 84.4 (Oct. 2012), e1–e9. doi: 1987), pp. 261–275. doi: 10.2307/1368479. url: 10.1016/j.anbehav.2012.04.026. url: http:// http://dx.doi.org/10.2307/1368479. dx.doi.org/10.1016/j.anbehav.2012.04.026.

8