Strange Tunes—Acoustic Variation and Character Displacement in a Tarsier Hybrid Zone

Natural hybridization provides valuable insights into evolutionary processes, such as speciation and the forces driving or hindering it. Sulawesi tarsiers Tarsius dentatus and T. lariang hybridize within a limited area, suggesting selection against hybrids. Their species- and sex-specific duet songs might serve as a premating barrier in sympatry, especially if differences are strengthened by character displacement. Individuals of mixed origin might face disadvantages if they inherit intermediate song traits. To shed light on the processes shaping this hybrid zone, we analysed 55 duet songs from within and outside the zone. For females and males, we identified temporal and frequency-related parameters that differ between species. We inspected hybrid songs for intermediate characteristics and analysed purebred songs for character displacement in sympatry. Female hybrid songs (N = 2) were intermediate in four to five of six parameters; interpretation of male hybrid songs (N = 2) was inconclusive, because only two parameters were reliably quantifiable. There was no character displacement in female songs in sympatry (N = 11) compared with monospecific areas (N = 17). In male songs, interspecific differences in note rate were significantly larger within the hybrid zone (N = 8) compared with outside (N = 13). Intermediate song traits indicate inheritance and may disadvantage hybrids during mate choice. Character displacement in male songs is consistent with female mate choice, because females should opt for unmistakable signals to avoid costly hybridization. Our findings thus suggest that duet songs of T. lariang and T. dentatus play an important role in limiting this hybrid zone.


Introduction
Speciation is based on the evolution of reproductive barriers between populations. When populations come into secondary contact before reproductive isolation is complete, hybridization-interbreeding-is likely to occur (Abbott et al., 2013;Brandler et al., 2021). Hybridization may thus provide insights into the evolution of species, making hybrid zones "natural laboratories for evolutionary studies" (Hewitt, 1988, p. 158). Therefore, hybridization has been in the focus of evolutionary research ever since Charles Darwin realized the significance of "crossing species" in the evolutionary context (Zinner et al., 2011). The consequences of hybridization depend on various factors, such as habitat properties, hybrid fitness, frequency of assortative mating, and backcrossings (Arnold, 1997;Cortés-Ortiz et al., 2019;den Hartog et al., 2007;Stebbins, 1959;Zinner et al., 2011). Hybridization can push, delay, or even prevent speciation (Abbott et al., 2013). If hybrids exhibit adaptive traits leading to fitness higher than or similar to the parent species (e.g., in the periphery of the latter's ecological range) and if backcrossings into the parental populations are common (den Hartog et al., 2007), the hybrid zone is likely to expand. As a result, both parental species will eventually be superseded or will fuse back into a single species 1 3 Strange Tunes-Acoustic Variation and Character Displacement… ("speciation reversal"; Seehausen, 2006). If hybrids exhibit fitness similar to or even better than the original species but mainly breed assortatively without crossing back, they may establish as a third species ("hybrid speciation") and coexist with both parent taxa (Abbott et al., 2013;Lamichhaney et al., 2018;Mallet, 2007;Zinner et al., 2011). Most frequently, however, hybridization in animals results in no offspring or offspring with lower fitness than the parent species and thus is selected against (Barton & Hewitt, 1985;Hoskin & Higgie, 2013). The latter case results in a stable, narrow hybrid zone, a so-called "tension zone" (Barton, 2001;Barton & Hewitt, 1985;Buggs, 2007;den Hartog et al., 2007;Merker et al., 2009) where backcrosses with the parental populations cause introgression of foreign genes from one species' gene pool into the other.
Selection against disadvantageous heterospecific matings reinforces reproductive barriers in the parental species (Hoskin & Higgie, 2013;Howard et al., 1993). These barriers will preferentially take effect before mating, promoting the recognition of conspecific mates (Abbott et al., 2013;Coyne & Orr, 2004;den Hartog et al., 2007;Ortiz-Barrientos et al., 2009;Servedio & Noor, 2003). An outcome of reinforced premating barriers might be displaced reproductive characters-the phenomenon of diverging phenotypes in sympatric populations of closely related species that experience reduced fitness due to reproductive interaction (Brown & Wilson, 1956;Pfennig & Pfennig, 2009;Wilkins et al., 2013). If premating barriers are not invincible, postmating barriers might contribute to drive reproductive isolation. Complications for individuals of mixed ancestry can be intermediate, possibly unattractive signals (El-Shehaby et al., 2011;Rundle & Nosil, 2005;Segura et al., 2011;Xue et al., 2018). As a result, hybrid ⨯ hybrid matings may be more common than backcrosses (Lamichhaney et al., 2018). However, intermediate signals can as well be more attractive to one or both parental species (Coyne et al., 1994;Vander Meer et al., 1985), which would increase the number of backcrossings.
In many vocalizing species, such as birds and nocturnal mammals, acoustic signals are prevailing traits for mate recognition and mate choice and thus can act as premating reproductive barriers. For this reason, acoustic mate choice signals are prone to reproductive character displacement (Braune et al., 2008;den Hartog et al., 2007;Kenyon et al., 2011): If vocal signals are used for species recognition and if hybridization is disadvantageous in terms of fitness, the calls of closely related species might diverge where populations co-occur. If an individual's vocal repertoire is inherited rather than being learned from its consexual parent or its social neighbourhood, such "strange tunes" of hybrids might be unfavourable in the course of mate choice. While reproductive character displacement in vocal traits has been observed in several animal taxa (insects: Marshall & Cooley, 2000;anurans: Gordon et al., 2017;Micancin & Wiley, 2014;Pfennig & Rice, 2014;birds: Demko et al., 2019;Kirschel et al., 2009;Seddon, 2005;mammals: Campbell et al., 2019), it has not yet been clearly demonstrated in primates. Evidence for ecological or reproductive character displacement occurs in visual signals in guenons (Allen et al., 2014) and ecological character displacement occurs in acoustic signals in tamarins (Sobroza et al., 2021).
One of the most complex and elaborate-and still little understood-forms of vocal communication in animals is duetting, as the necessary coordination is 1 3 cognitively highly demanding (Nieder & Mooney, 2019). Duets are defined as overlapping, more or less precisely coordinated and stereotype song displays of two, usually mated, individuals-a leader/initiator and a follower/responder (Dahlin & Benedict, 2014;Farabaugh, 1982;Hall, 2004;Logue & Krupp, 2016). Duet songs often are loud and emitted mainly at dawn, thus optimizing long distance transmission (Adret et al., 2018;Seibt & Wickler, 1977); they are mostly sex-specific (Todt & Naguib, 2000). Beside the emitters' locations, details transferred in a duet can comprise crucial information for mate choice, such as species and individual identity, sex, reproductive status, mate quality, or pair-bond-status (Hall, 2004). Duetting might have more than one function, such as joint resource (e.g., territory) defence and pair bond strengthening or mate guarding, even within a single species (Hall, 2004).

3
Strange Tunes-Acoustic Variation and Character Displacement… While, despite the recognition of its evolutionary significance, our knowledge about primate hybridization is in general fragmentary , this particularly applies to tarsiers (Zinner et al., 2011). The first evidence of natural hybridization in tarsiers concerns the two Eastern tarsier species Tarsius lariang and T. dentatus . Both species are endemic to Sulawesi, where they are parapatrically distributed. They are similar in morphology, ecology, and behaviour; they sleep in small groups (mostly one adult pair plus their offspring), and their mating system can be described as predominantly monogamous with occasional extra-pair matings Merker et al., 2004;Tremble et al., 1993;Bohr, unpublished data). T. lariang and T. dentatus diverged ca. 1 Mio. years ago (Driller et al., 2015) and came into secondary contact before their reproductive isolation was completed. Genetic evidence points to a narrow hybrid tension zone ) and thus to selection against hybrids. Duet songs might play an important role in shaping the hybrid zone by acting as a reproductive barrierbefore mating (possibly strengthened by character displacement) or after mating in terms of "strange tunes", i.e., unattractive signals, if songs show intermediate characteristics (Rundle & Nosil, 2005).
We focussed on a narrow contact zone of ca. 1 km width between Tarsius lariang and T. dentatus  to further elucidate the role of duet songs in tarsier speciation processes. Subsequent to an assessment of the discriminative ability of parameters measured for male and female tarsier songs, we addressed the following questions: 1) Do individuals of mixed ancestry show purebred or intermediate song traits? As vocal learning in nonhuman primates is highly limited, we hypothesize that also tarsier songs are to a large extent inherited. We thus predict intermediate traits for the songs of hybrids compared with those of purebreds. If tarsiers learn their vocal repertoire from their consexual parent, then songs of hybrid individuals should clearly resemble those of purebreds. 2) Is there evidence for vocal character displacement in male or female songs in areas of sympatry compared with allopatric occurrences, i.e., monospecific areas? Because females have higher reproduction costs, they should opt for clearly discernible signals during mate choice. Therefore, we predict character displacement in male but not or less in female songs.

Study Sites and Sampling
Between 2005 and 2012, we recorded duet songs of Tarsius lariang and T. dentatus at seven different locations in Central Sulawesi, Indonesia ( Fig. 1; Tables S1, S2), at several monospecific sites and within a hybrid zone described by Merker et al. (2009).
To account for a possible movement of the hybrid zone over time, we extensively monitored its soundscape in 2005/2006 and 2012. Yvonne Bohr consistently spent a full year in the hybrid zone, and Stefan Merker spent 7 months there in total; fieldwork that 1 3 included song monitoring took place almost every day. We were thus able to observe in each period where tarsier groups were living in audible range of calls of the other species (considered as hybrid zone) or not (considered as monospecific area). Because we monitored each group several times, we were highly familiar with the local soundscape. We only included vocalizations in the study whose assignment to the hybrid zone or to a monospecific area was unambiguous. We recorded all but one duet between 0500-0630 h when tarsiers rejoined at their sleeping sites. One recording took place at 1735-1745 h when the group dispersed. Recording distance ranged between ca. 5 and 30 m.

Recording Equipment and Acoustic Analysis
Stefan Merker recorded tarsier songs at locations 1-6 ( Fig. 1). Two recordings from Luwuk (location 7) were provided by Christine Driller. Both used a Sony MZ-NH900 Hi-MD Walkman (Sony Corporation, Tokyo, Japan) with a frequency range of 20 Hz -20 kHz, connected to a hand-held Røde NT3 condenser microphone (RØDE For site-specific sample sizes, see Table S2.

3
Strange Tunes-Acoustic Variation and Character Displacement… Microphones, Silverwater, NSW, Australia), covering a range of 20 Hz -20 kHz. After recording, they converted linear PCM files into Waveform Audio File Format (.wav) using Sony's SonicStage 4 software. Yvonne Bohr recorded tarsier songs at Winatu (location 1; Fig. 1). She employed an Olympus LS-3 Linear PCM recorder (Olympus Corporation, Tokyo, Japan), with a frequency range of 40 Hz -21 kHz, connected to the handheld, highly directional gun microphone Sennheiser ME 67 (Sennheiser electronic GmbH & Co. KG, Wedemark, Germany), which covers a frequency range of 40 Hz -20 kHz. Conversion of files was not needed here. We sampled all signals at a rate of 44.1 kHz and digitized them with 16-bit resolution, resulting in a dynamic range of 96 dB. As the dominant frequencies of T. lariang and T. dentatus duets in general range between 3 and 20 kHz (Merker & Groves, 2006;Niemitz et al., 1991;Grow, 2019; see Gursky, 2015 for nonduet ultrasonic vocalizations in Sulawesi tarsiers), our equipment adequately covered the frequency range of the vocalizations examined in this study. We visualized and analyzed the final stereo .wav files in Raven Pro 1.4 and 1.5 (The Cornell Lab of Ornithology, 2011). We only analyzed spectrograms of completely recorded duet songs where all relevant song parameters were clearly discernible, which lowered our sample size considerably. We generated spectrograms in Raven Pro with a fast Fourier transformation (FFT) applying a Hann window shape (sine-squared) with a size of 256 samples, 3 dB bandwidth of 248 Hz, an overlap of 50% and a frequency grid spacing of 172 Hz. The term "duet" refers to the conjoined vocalization of female and male tarsiers. We refer to one individual's part in a duet as "song" or "duet song". A "note" is the basic unit of a song and consists of a single, continuous, up-or downwardly modulated trace in a spectrogram. We analyzed 32 songs from female and 23 songs from male tarsiers. All songs stem from different individuals to ensure independence of the samples. Due to pronounced divocalism (Fig. 2;MacKinnon & MacKinnon, 1980;Merker & Groves, 2006;Niemitz et al., 1991;Nietsch, 1999), it was easy to distinguish female and male songs, both audibly and visually. For female songs, we manually drew selection frames around the emphasized frequency. This energy-richest or dominant frequency of a harmonic series is in general identical to the fundamental or first harmonic, referring to the lowest frequency in a harmonic series (Whitehead, 1995). We then extracted the following measurements: begin time, end time, delta time, lowest frequency, highest frequency, delta frequency. From these measurements, we calculated six parameters applying to the whole duet song or to single notes (Table I). Because we found notes of male songs difficult to visually identify and to frame consistently over a whole song, we restricted framing to the 10 central notes of each male song (the number of notes in male songs ranged between ca. 35 and 240 notes). This limited the number of parameters compared with female songs. We extracted begin time, end time, lowest frequency, and highest frequency of each note, and we then calculated the number of notes per second and the mean frequency height of notes (Table I). In tarsier duet songs, notes of males and females partly overlap. We therefore refrained from an automatic placement of selection frames but framed manually instead. To minimize bias from placing selection frames in Raven Pro by eye (and thereby to reduce bias in the extracted parameters), a single researcher (Yvonne Bohr) checked and, if necessary, adjusted the final placement of frames.

Assignment of Songs
Before acoustic and statistical analyses, we assigned songs to sex, species (or hybrid origin), and area. We could easily distinguish male and female duet songs based on spectrograms. To assign female and male songs to either Tarsius lariang, T. dentatus, or hybrid origin, we used geographic information (calls close to or far from the common species border, i.e., hybridization possible or unlikely) and-in the hybrid zone-genetic data (microsatellite-based assignment scores) and acoustic monitoring findings. We only labelled songs to be of hybrid origin if the emitting animal was unambiguously identified and genotyped as a hybrid. Both female and both male hybrids were adult at the time of recording. We labelled songs that we could not assign to T. lariang, T. dentatus, or hybrid origin as "unassigned". If not stated otherwise, samples used in statistical analyses comprise assigned (not unassigned) songs. Concerning the area, we assigned recordings as stemming from the hybrid zone if T. lariang and T. dentatus both occurred in that area, i.e., when we had heard duets of both species. We denoted recordings as stemming from a monospecific area if we had heard vocalizations of only one of both species. For more details of the assignment procedure, see Supplement.

Statistical Analyses
We performed all statistical analyses using IBM SPSS Statistics, Version 26.0. If not stated otherwise, tests were two-tailed with α set at 0.05. We report means with standard errors ( x ± SE). We tested data for normal distribution using 1 3

Strange Tunes-Acoustic Variation and Character Displacement…
Shapiro-Wilk test and by visually examining Q-Q plots. To test homogeneity of variances, we used Levene's test based on medians and visually inspected scatterplots of standardized predicted values versus standardized residuals (zpred vs. zresid).
We used independent t-tests to test for interspecific differences in song parameters in both sexes. As appropriate for small sample sizes, we report effect size for t-tests as Hedge's g* (Hedges & Olkin, 1985). For a qualitative interpretation of g*, we use Cohen's (1992) and Rosenthal's (1996) classification (0.2 being a small effect, 0.5 being a medium effect, and 0.8 being a large effect).
To identify the most significant song parameters discriminating between Tarsius lariang and T. dentatus females, to classify yet unassigned female songs, and to evaluate female hybrid songs, we applied a linear discriminant function analysis (DFA), embedding the six variables (Table I) in a stepwise manner. We opted for the stepwise method in order to exclude variables that do not enhance the discriminant function's power. Variables entered the stepwise DFA according to their impact on Wilks' lambda with "F to enter" = 3.84 and "F to remove" = 2.71. Classification results are based on equal prior probabilities. To examine the classification's robustness, we applied a leave-one-out cross-validation. For female songs, we used the described assignment as "actual group membership". We interpreted squared canonical correlation coefficients (R c 2 ) as effect sizes for DFA. Because we used only two parameters for male tarsiers (Table I), we evaluated male hybrid songs by examination of scatterplots and raw data of male song characteristics.
To assess potential reproductive character displacement in the hybrid zone, we tested whether the difference between purebred T. dentatus and T. lariang songs was larger for individuals in the hybrid zone (HZ) compared with monospecific areas (MSA). We therefore created a new variable ("displacement"), calculated as the absolute difference between an individual's parameter value and the other species' mean for this parameter in the respective area. We compared "displacement" between the monospecific areas and the hybrid zone by means of a one-tailed, independent samples t-test where H 0 : displacement HZ ≤ displacement MSA, and H a : displacement HZ > displacement MSA.
In addition, to examine the details of intraspecific character displacement, we used Mann-Whitney U tests to compare MSA with HZ for each song parameter. We report r as effect size estimate with 0.1 being a small effect, 0.3 a medium effect, and 0.5 a large effect (Fritz et al., 2012). For correct interpretation of effect sizes, we provide formulas for g* and r in the Supplement.

Ethical Note
This study complied to the legal requirements for foreign researchers in Indonesia In addition, our research followed the IPS Code of Best Practices for Field Primatology; all methods for data acquisition were non-invasive. The authors affirm that they have no conflict of interest with any entities described herein.  Fig. S1. Detailed descriptions of the species' spectrograms can be found in Niemitz et al. (1991), Nietsch (1999), and Nietsch and Kopp (1998) for T. dentatus and for T. lariang in Merker and Groves (2006).

3
Strange Tunes-Acoustic Variation and Character Displacement…

Female Songs
In total (monospecific areas and hybrid zone combined), we assigned seven female songs to Tarsius dentatus, 14 songs to T. lariang, and two songs to females of mixed origin (Hybrid 1 and 2). Songs of nine individuals remained "unassigned" (Table S2).
Songs of female T. dentatus (N = 7) and T. lariang (N = 14) differed significantly in all six song parameters examined, with very strong effect sizes (Table II; Fig. S2a-f). T. dentatus females uttered shorter notes at a higher rate; their songs started at a higher frequency, they reached higher frequencies over the entire song, and both their songs and single notes had a larger frequency range (Table II) (Fig. 3) retained three predictor variables that discriminate best between populations: Max. frq. of 1st note, Mean note length, and Mean frq. span of notes. Goodness of fit statistics show that predictor variables significantly discriminate between groups (Table III). Discriminant function (DF) 1 can be regarded as the "frequency function", DF2 as the "temporal function". Effect sizes (R c 2 ) point towards strong effects, especially for the "frequency function" DF1.
Based on original cases, DFA classified all individuals correctly. The nine unassigned individuals were classified as follows: two as Tarsius dentatus; five as T. lariang; and two as hybrids (Table IV; Fig. S3a-f). Cross-validation classified 96% of all cases correctly. It assigned all T. dentatus and all hybrids to the right group, but classified one T. lariang as a hybrid (Table IV).

Character Displacement in Female Songs
We recorded 17 female songs in monospecific areas. Of these, we assigned five to Tarsius dentatus and 12 to T. lariang. In the hybrid zone, we originally assigned only two female songs to T. dentatus and two songs to T. lariang. In order to increase sample size for the hybrid zone, we included the initially unassigned recordings according to their DFA classification. This resulted in 11 songs for the hybrid zone, including four T. dentatus and seven T. lariang songs. The one-tailed t-test to compare "displacement" between the monospecific areas and the hybrid zone did not reveal reproductive character displacement in any of the six parameters (i.e., "displacement" was not larger in the hybrid zone; Table V; Fig. S5a-f).

Male Songs
We assigned (monospecific areas and the hybrid zone combined) 13 male songs to Tarsius dentatus, eight songs to T. lariang, and two songs to males of mixed origin (Hybrid 3 and 4); see also Table S2. Male duet songs differed significantly between species: Songs of T. lariang males were faster and lower in frequency than songs of T. dentatus males ( Fig. 4; Table VI; Fig. S6a-b). The effect of species identity was very strong for both parameters, but it was larger for frequency height than for note rate. As sample size for male hybrids was too small (N = 2) for statistical analyses, we give only descriptive values. While Table IV Classification results of stepwise discriminant function analysis on female T. dentatus, T. lariang, and hybrid song parameters. Prior probabilities calculated as "all groups equal". We recorded songs 2005-2012 in Central Sulawesi, Indonesia.

Group
Predicted group membership
Most of the shift in Notes/s can be attributed to Tarsius lariang (Fig. 6). In songs of male T. dentatus, a Mann-Whitney U test showed no significant difference in Notes/s between monospecific areas (N = 9, M rank = 7.33) and the hybrid zone (N = 4, M rank = 6.25), U = 15.00, Z = − 0.46, P = 0.71 (using the exact sampling distribution of U, Dinneen & Blakesley, 1973) and a small effect of r = 0.13. For T. lariang, although there was no statistically significant difference between monospecific areas (N = 4, M rank = 3.00) and the hybrid zone (N = 4, M rank = 6.00), U = 2.00, Z = − 1.73, P = 0.114, there was a large effect of r = 0.61.

Discussion
In this study, we determined specific vocal traits for males and females that distinguish the duet songs of Tarsius lariang and T. dentatus. We also provide an initial assessment of how songs of tarsier hybrids relate to songs of parental species. Songs of female hybrids were intermediate in most parameters. A clear interpretation of male hybrid songs was not possible. While Hybrid 4 was intermediate in frequency (lying between confidence intervals), Hybrid 3 grouped with T. lariang. For note rate, Hybrid 4 grouped with T. lariang, while Hybrid 3 was intermediate. Furthermore, we provide evidence for asymmetric reproductive character displacement in a hybrid zone of nonhuman primates: T. lariang males seem to sing faster in the mixed zone compared to outside, whereas the song characteristics of females and T. dentatus males are not reinforced in the hybrid zone.
Our analyses of female and male duet songs illustrate the discriminative power of selected temporal and frequency related song traits. They further corroborate the significance of duet songs in discriminating tarsier species. The notion that duet songs of Sulawesi tarsiers are species-specific has long been held (MacKinnon & MacKinnon, 1980;Niemitz, 1984). Commonly, identification of new tarsier species begins with the detection of unfamiliar vocalizations in the field and their comparison with spectrograms of known species. Morphological and anatomical studies follow (Merker & Groves, 2006;Shekelle et al., 2008Shekelle et al., , 2017Shekelle et al., , 2019, eventually combined with or followed by genetic analyses (Driller et al., 2015;Merker et al., 2010;Shekelle et al., 2010). As was already known from the descriptions of Tarsius dentatus (formerly T. dianae, Niemitz et al., 1991) and T. lariang (Merker & Groves, 2006) and as visual evaluation of spectrograms suggest, our analyses of single parameters prove songs of female T. dentatus and T. lariang to be highly distinct. There are striking differences in temporal as well as frequency-related traits. Discrimination of male songs, both in the field and visually, is less obvious. Our study provides two easily measurable song parameters for male T. dentatus and T. lariang-one temporal and one frequency-related trait-to discern male songs clearly.
T. dentatus females and males sing at higher frequencies compared with T. lariang. The use of different frequency levels might relate to different habitat preferences of species, as vocalizations at lower frequencies reach farther in dense vegetation than higher frequencies (Brown & Waser, 2017;Brown et al., 1995;Masali et al., 1992;Peters & Peters, 2010). So far, however, there is no evidence of striking differences in natural habitats used by T. lariang and T. dentatus Yvonne Bohr and Stefan Merker, personal observation). Vocalization frequency might be also related to body size or body mass of the different taxa; in mammals, body size often is negatively correlated with frequency Martin et al., 2017;Peters & Peters, 2010). Correspondingly, species with a smaller head size are better at locating signals of higher frequencies (Heffner, 2004;Masterton et al., 1969;Ramsier & Rauschecker, 2017). Although T. lariang apparently has a slightly larger mean skull length compared with T. dentatus (Merker & Groves, 2006), there were no significant differences in head length, nor in head-body-size in our sampling of both species in Winatu (Yvonne Bohr, unpublished data). While T. dentatus males had a higher body mass than T. lariang males, females did not differ. Therefore, the correlation effect of body mass and frequency does not seem to be applicable here. The effect of body size on the fundamental frequency is based on the correlation between body size and larynx size: a larger body holds a longer larynx, which produces sounds of lower fundamental frequency . However, the size of the larynx can be decoupled from body size . Therefore, a shorter larynx of T. dentatus compared with that of T. lariang might explain the higher frequencies in former species' songs. The true rationale behind the use of different frequencies in T. dentatus and T. lariang remains unknown for now. The clear interspecific differences in the vocal output of females and males ( Fig. 2; Burton & Nietsch, 2010), however, might facilitate the tarsiers' mate choice or avoid interspecific conflicts in areas of contact, such as in the studied hybrid zone. Playback experiments could further unveil the role of tarsier songs in these and other contexts.
Genetic analyses done after the song recordings showed that we had good quality records from only two female and two male hybrids. This sample size is too low for sound statistical analyses. Therefore, the performed statistical analysis (DFA) and our interpretation of hybrid songs have to be considered preliminary. While we tentatively labelled one of the genetically identified female hybrids (Hybrid 1) as T. lariang during monitoring in the field, the other female hybrid (Hybrid 2) had already attracted attention during monitoring, because her song appeared to start like T. dentatus and to end like T. lariang. Likewise, we already noted both previously "unassigned" songs of non-genotyped females that DFA later characterized as hybrid songs to be peculiar in the field. Acoustic and visual inspection of spectrograms as well as statistical analyses point towards a song structure of female hybrid calls that differs from their parental species. A comparison of their parameter means with confidence intervals of T. lariang and T. dentatus and the DFA point towards an intermediate vocal pattern. While songs of male hybrids partly showed intermediate values with respect to the parent species, the two parameters analyzed are not sufficient to characterize male hybrids' songs clearly. Especially because male tarsier notes are in general more similar between species than are female notes.
There are several examples of intermediate song traits in avian hybrids (de Kort et al., 2002;Moore & Coulson, 2020;Shipilina et al., 2017). For mammals, intermediate vocal traits have been shown for mice (Hahn et al., 1998), fur seals (Page et al., 2001), deer (Long et al., 1998), and primates. In nonhuman primates, intermediate song traits are so far known only from gibbons (Brockelman & Schilling, 1984;Geissmann, 1984;Maples & Haraway, 1982;Marler & Tenaza, 1977;Tenaza, 1985) and howler monkeys (Kitchen et al., 2019). As in ours, the sample sizes in those studies were low, each including one or two hybrid individuals. A female Hylobates muelleri (father) ⨯ H. agilis (mother) hybrid's song resembled the females' song of her father's species, which she had never heard before (Maples & Haraway, 1982). While some song traits of a male H. pileatus (father) ⨯ H. lar (mother) hybrid were very similar to purebred songs of either parental species, others were intermediate or even showed new structures. In the female hybrid's song, some traits equaled those of her father's species, but others were intermediate (Geissmann, 1984). The song of a female H. lar (father) ⨯ H. muelleri (mother) hybrid differed from both parental songs, whereas the male hybrid's song was lar-like except for one trait that resembled H. muelleri (Tenaza, 1985). A genetically intermediate male hybrid of Alouatta palliata and A. pigra vocalized like A. palliata in temporal traits, while frequency-related song characteristics were intermediate (Kitchen et al., 2019). For tarsiers, our results show a similar "mixed" pattern: the two female hybrids were, compared with confidence intervals of purebred T. dentatus and T. lariang, intermediate for the majority of the parameters. For the non-intermediate song characters, there was individual variation in the similarity to purebreds. One male hybrid was intermediate in the temporal parameter and grouped with T. lariang concerning the frequency trait. The other male hybrid was intermediate in the frequency parameter and grouped with T. lariang concerning the temporal trait.
Intermediate song traits indicate genetic inheritance, rather than learned vocal patterns (Page et al., 2001). Across the animal kingdom, there is far less evidence for learned vocalizations, i.e., by imitation of other individuals, compared with genetically determined sound, even in birds (Boves et al., 2010;Janik & Slater, 1997;Nieder & Mooney, 2019;Päckert, 2018;ten Cate, 2021;Woolley & Sakata, 2019). In mammals, the clearest evidence of vocal production learning, i.e., "the ability to modify the structure of vocalizations as a result of hearing those of others" (Janik & Knörnschild, 2021), which includes modifications of frequency properties (Boughman & Moss, 2003;Janik & Slater, 1997), derives from cetaceans (Janik, 2014), pinnipeds (Reichmuth & Casey, 2014), bats (Knörnschild, 2014;Vernes & Wilkinson, 2020), elephants (Stoeger & Manger, 2014), and the human genus (Janik & Slater, 1997;ten Cate, 2021;Tyack, 2020). In nonhuman primates, this ability is debated and seems to be, at best, very limited (Fischer & Hammerschmidt, 2019;Janik & Knörnschild, 2021;Lameira et al., 2013). Several researchers concluded that vocalizations of nonhuman primates are to a large extent innate and have a highly conserved structure within species (Cheney & Seyfarth, 2018;Fischer & Hammerschmidt, 2019;Fischer & Price, 2017;Janik & Knörnschild, 2021;Jürgens, 2002;Owren et al., 1992). Neurobiological studies suggest that nonhuman primates are not able to produce new sounds (in terms of voiced calls) as they lack a direct connection between the primary motor cortex and the laryngeal motoneurons (Cheney & Seyfarth, 2018;Jürgens, 2009). Our analysis gives a first glance at vocalization patterns in hybrid tarsiers and points towards intermediate song traits. This finding suggests that tarsier duet songs are largely innate and that tarsier hybrids may not be capable of learning purebred structured vocalizations from their consexual parent. This is in accordance with current knowledge on vocal learning in nonhuman primates and on the acquisition of their vocal repertoire.
In the hybrid zone of T. lariang and T. dentatus, purebreds may perceive intermediate songs, or more precisely songs with intermediate characteristics, as "strange tunes". Not being perceived as attractive mates will result in fitness disadvantages for hybrids ("extrinsic behavioral hybrid dysfunction"; Irwin, 2020;Servedio & Noor, 2003). In various taxa, purebred females respond less to hybrid vocalizations compared with those of conspecific males (e.g., anurans: Doherty & Gerhardt, 1984;Höbel & Gerhardt, 2003;birds: Derégnaucourt & Guyomarc'h, 2003). Nevertheless, responses to mixed individuals' signals can still be stronger than to heterospecific individuals due to higher similarity or even single "sexy traits" in the intermediate vocalization (Wyman et al., 2016;e.g., anurans: Gerhardt, 1974;Höbel & Gerhardt, 2003;birds: Collins & Goldsmith, 1998;Derégnaucourt & Guyomarc'h, 2003). The resulting reduced discrimination against hybrids enables backcrossings and can lead to intensive introgression (Randler, 2002;Wyman et al., 2016). In the case studied here, the limited width of the hybrid zone between T. lariang and T. dentatus  suggests selection and thus discrimination against hybrids. It remains unclear whether the intermediate/mixed character of hybrid songs goes along with an increased attractiveness to parental species as compared to heterospecific songs, possibly promoting backcrossing.
A larger number of good quality recordings from mixed tarsier individuals, at best in combination with detailed genetic information, would provide more insights into how and to what extent hybrid vocalizations differ from purebreds. As species have, over time, acquired characteristic vocal repertoires and mate recognition systems, purebred song traits are well-predictable over time and space. In contrast, the multitude of ways hybridization can take provides for a great potential of variability in hybrid song traits. An increased sample size of hybrid songs is therefore important. We also encourage future studies to investigate how hybrid tarsiers perform with respect to other characteristics of their duetting performance, such as timing or coordination of song length (Clink et al., 2020). Playback experiments to explore perception of hybrid, con-and heterospecific signals of the same or the opposite sex by purebred and mixed individuals could shed more light on premating reproductive isolation between T. lariang and T. dentatus in the common hybrid zone.
The accentuated interspecific contrast in male song characteristics in the hybrid zone, compared with monospecific sites, is consistent with a key role of tarsier song in mate recognition. Post-mating barriers to gene flow, e.g., reduced fitness of hybrids, can lead to selection of premating barriers, like the reinforcement of mate choice relevant signals. Our finding of character displacement in males but not in females might be interpreted as a result of female mate choice; due to their higher investment in reproduction (Gursky, 2002), female Sulawesi tarsiers should opt for high quality males, i.e., for good genes. If mating with hybrids is disadvantageous, purebred females should choose males with unmistakable signals. Their preference for vocalizations that can be clearly assigned to a conspecific male would then drive character displacement. For males, in contrast, the large interspecific differences in evolved female songs may simply make further differentiation unnecessary. In addition, males generally encounter fewer disadvantages from unfavorable matings. In view of observed extra-pair offspring Yvonne Bohr, unpublished data), this might hold true for T. lariang and T. dentatus. In addition, contrary to the predictions of reproductive character displacement, female songs seem to be more similar within the hybrid zone than in monospecific areas. This may be due to a lack of need for reinforcement, on the one hand, and introgression in the course of hybridization, on the other Moran & Fuller, 2018).
Interestingly, T. lariang males seem to be responsible for the notable shift in Notes/s (Fig. 6), which is underlined by the effect size estimate of r = 0.61 compared with 0.13 in T. dentatus. T. lariang males in general "sing faster" than T. dentatus males, i.e., they show higher values for Notes/s. In the hybrid zone, they seem to sing even faster compared to monospecific areas, whereas T. dentatus males show no difference. This might indicate that T. lariang females have a greater interest in avoiding hybridization than T. dentatus females and choose conspecific males that differ most clearly from T. dentatus. The rationale behind this asymmetric pattern might include differences in selection pressure against hybridization between the two species, possibly because of "unequal hybridization costs, biases in likelihood of hybridization, asymmetrical effects of reproductive interference, evolutionary constraints, or historical accidents" (Cooley, 2007). As Smadja and Ganem (2005) reason for house mice, the potential competitive asymmetry in favour of the probably more opportunistic T. dentatus compared with T. lariang ) might promote migration of T. dentatus into T. lariang range, which would favour stronger selectivity in female T. lariang. In view of previous observations (sampling period 2001-2006) that T. dentatus males successfully mate (or have mated) with T. lariang females but no indication that T. lariang males reproduce with T. dentatus females , our findings of stronger (but probably not yet sufficient) discrimination of T. lariang females against T. dentatus males may thus be interpreted as a reaction to the invasion of the T. lariang range first and foremost by T. dentatus males (and only subsequently by T. dentatus females). This tentative movement of the hybrid zone further into T. lariang range might explain the asymmetric pattern of mitochondrial introgression between the two species . Further genetic analyses and playback experiments may shed more light on this matter.
The observed displacement in male song traits in the hybrid zone concerned the temporal but not the frequency parameter. This matches the prediction that temporal traits can be easier modified than frequency-related traits (Boughman & Moss, 2003). The former can be modified by behavioral changes and are thus under control of the singing individual, whereas the latter require control of the vocal apparatus, which is restricted in nonhuman primates (Janik & Slater, 1997ten Cate, 2021).
The results of this study are in accordance with previous observations that tarsier vocalizations play an important role in mate recognition (Burton & Nietsch, 2010;Merker et al., 2009;Nietsch, 1999;Shekelle et al., 1997). In summary, we identified temporal and frequency-related parameters to discriminate between duet songs of Tarsius lariang and T. dentatus and likely between other Sulawesi tarsier species as well. We report intermediate patterns of hybrid tarsier vocalizations in terms of quantitative acoustic characters and provide evidence of asymmetric reproductive character displacement in male songs. These findings suggest that duets songs play a key role in directing hybridization between T. lariang and T. dentatus.
Strange Tunes-Acoustic Variation and Character Displacement… are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.