TY - JOUR AU - Iorgu, Ionuţ, Ştefan AB - Abstract Divergence of the acoustic signals used in mate calling may be an important driver of speciation. Male song may vary accordingly to, in mismatch with or independently of female preferences. Therefore, to estimate the importance of male signal variation between subspecies, female preferences against subspecies-specific signal variants must be tested. We examined the female response probability in the two subspecies of the bush-cricket Isophya kraussii for consubspecific and heterosubspecific male signals, and also for the song of a closely related species (I. camptoxypha), which is sympatric with one, but allopatric to the other subspecies. Performing no-choice playback experiments, we found that females of both subspecies responded to the male song of their own subspecies with significantly higher probability than to heterosubspecific and heterospecific songs. Response specificity for consubspecific vs. heterosubspecific signals was not significantly asymmetric comparing the two subspecies. A significant difference was found, however, in the discrimination of heterospecific calls: females showed stronger discrimination against the song of I. camptoxypha in the subspecies sympatric with that species. Our results are best explained by a tightly covarying, stepwise coevolution of male signals and female preferences. acoustic communication, mate recognition, signal evolution, speciation INTRODUCTION Interpopulation divergence in sexual acoustic communication can be a crucial element of speciation (Wilkins, Seddon & Safran, 2013). It may be either an initial step when sexual acoustic communication evolves along different trajectories in allopatric populations due to sexual selection (e.g. Grace & Shaw, 2011; de Oliveira Gordinho et al., 2015), ecological constraints (Zuk, Simmons & Cupp, 1993; Ballentine, 2006; Cocroft, Rodriguez & Hunt, 2010), genetic drift or mutation (Campbell et al., 2010; Goodman et al., 2015) and cultural evolution (Lachlan & Servedio, 2004; Linossier et al., 2016) or it may arise during the secondary contact of previously isolated populations diverged in traits unrelated to sexual communication during allopatry. In those cases, divergence in sexual acoustic communication may take part in the reinforcement of the genetic isolation between them (Gerhardt, 2013). The resulting interpopulation divergence of male signals is often the first detected sign of the process of divergence. Thus, the examination of acoustic signals in morphologically similar populations led to the discovery of ‘cryptic’ species in several taxa of animals (Henry, 1994; Jones, 1997; Ng et al., 2016; Heller et al., 2017). However, to estimate the significance of signal divergence regarding reproductive isolation and to get closer to the understanding of the evolutionary process causing the divergence, we also need to examine the receiver side of the communication system: female preferences against male signal variants. If females prefer the male song variant of their own population, then it is reasonable to presume that female preferences are the dominant drivers of male signal divergence (Gray & Cade, 2000; Grace & Shaw, 2011; Barbosa, Rebar & Greenfield, 2016) and the differences in acoustic communication are expected to be the significant components of prezygotic isolation between the populations. However, if females do not show definite preference for the male signal variants of their own population, then strongly directional sexual selection, sensory exploitation (Ryan & Rand, 1993) or processes other than intersexual selection are likely to have significant roles in the male signal divergence (Simmons, Zuk & Rotenberry, 2001; Velásquez et al., 2015). In those cases, the divergence in acoustic communication would not bar interpopulation crossings in case of their secondary contact. Acoustic signalling is a characteristic behaviour of many species in the insect order Orthoptera (Robinson & Hall, 2002). Males emit their acoustic signals spontaneously to attract distant females or to stimulate their acoustic response. The majority of studies in which female preferences for male song patterns were experimentally investigated found definite preferences for conspecific male signals (e.g. Perdeck, 1958; Stumpner & von Helversen, 1994; Guerra & Morris, 2002). However, in some species, the behavioural response of females against different male song variants has proved to be non-discriminative (Heller, von Helversen & Sergejeva, 1997; Bush & Schul, 2010). Thus, the high level of interspecies acoustic diversity is paired with a hidden variation in female response specificity (Schul, von Helversen & Weber, 1998; Kowalski & Lakes-Harlan, 2011). Female response specificity may be even more variable during speciation. In this study, we examined female preferences for subspecies-specific male songs and heterospecific calling songs in the two subspecies of the bush-cricket Isophya kraussii Brunner von Wattenwyl. Isophya is the second most species-rich genus within the subfamily Phaneropterinae in Europe (Heller et al., 1998). Most species in the genus show subtle morphological differences, but differ conspicuously in their male calling songs (Heller et al., 2004; Chobanov et al., 2013; Zhantiev, Korsunovskaya & Benediktov, 2017), suggesting that acoustic diversification is a dominant component of the evolution of species richness in this genus. In spite of the conspicuous acoustic diversity of Isophya, experimental studies testing the effectiveness of male song differences as premating barriers in the genus are rare (Zhantiev & Dubrovin, 1977; Zhantiev & Korsunovskaya, 1990; Orci, 2007). As it is usual in the subfamily of Phaneropterinae, Isophya males and females perform an acoustic duet during their mate-finding behaviour: females emit short acoustic signals (a single click or group of a few impulses) in response to the male calling songs (Heller, 1990). The impulse repetition pattern of those short response signals is variable, but females emit their response during a species-specific time window in relation to the male signal. Among the European species of the genus, I. kraussii has the largest distribution range (Heller et al., 2004; Cigliano et al., 2017), and recently, Iorgu & Heller (2013) described a new subspecies, I. kraussii moldavica, near the eastern border of the area of the nominotypic subspecies. The two taxa differ mainly in the rhythmic features of the male calling songs: I. kraussii moldavica males produce significantly shorter syllables (see Material and Methods for the bioacoustic terminology applied in this study) containing a smaller number of impulses in comparison to the male song of the nominotypic subspecies (Iorgu & Heller, 2013). In addition to that difference, the male calling songs are very similar in oscillographic structure in the two subspecies: calling males produce long sequences of evenly repeated syllables (Fig. 1A, C), where each syllable is composed of a main impulse series and a few after clicks (Fig. 1B, D). Male songs of the two subspecies have similar, wide-band power spectra containing dominant frequency components between 20 and 40 kHz (Iorgu & Heller, 2013). The main aim of our study was to examine how much the attractiveness of subspecies-specific male calling songs differs for females in both subspecies. The examination of that question may help us to estimate whether the divergence of acoustic communication could cause a non-random mating pattern in case of a secondary contact of the two taxa. Furthermore, our results provide some important pieces of information about the signaller–receiver coevolution in this insect. Another question examined in this study was whether females in the two subspecies of I. kraussii discriminate against the heterospecific calls of Isophya camptoxypha (Fieber), a species sympatric with one subspecies, but allopatric with the other (Kenyeres & Bauer, 2005; Bauer & Kenyeres, 2006; Iorgu & Heller, 2013). The male calling song in that species is composed of the same oscillographic elements as in I. kraussii (Fig. 1E, F) and has a similar, wide-band frequency spectrum (15–45 kHz), but differs in syllable duration, syllable repetition rate and in the number of impulses per syllable (Orci, 2007). The attractiveness of that heterospecific signal was examined to test acoustic heterospecific discrimination in I. kraussii and to see whether the two subspecies differ in that. Figure 1. Open in new tabDownload slide Oscillograms showing the rhythmic patterns of male calling songs in the two subspecies of Isophya kraussii (A–D) and in I. camptoxypha (E, F) at two time scale resolutions. A, C, E, series of syllables; B, D, F, single syllables. Oscillograms on the right are time axis magnified portions (as indicated by dashed line rectangles) of the left ones. Time scale bars at the bottom are for all oscillograms in the given column of oscillograms. Figure 1. Open in new tabDownload slide Oscillograms showing the rhythmic patterns of male calling songs in the two subspecies of Isophya kraussii (A–D) and in I. camptoxypha (E, F) at two time scale resolutions. A, C, E, series of syllables; B, D, F, single syllables. Oscillograms on the right are time axis magnified portions (as indicated by dashed line rectangles) of the left ones. Time scale bars at the bottom are for all oscillograms in the given column of oscillograms. MATERIAL AND METHODS Study animals In both subspecies of I. kraussii, females were collected as large nymphs from their natural habitats. Isophya k. kraussii females originated from meadows and forest clearings in the Pilis Mountains (Hungary), and I. k. moldavica females were collected in the hillside meadows of Suceava county (Romania). All of them were reared to adults without any contact with males, nor could they hear any male signal. That way we tested virgin females with high motivation level concerning mate finding. Playback tests were done successfully with 16 females in I. k. kraussii and with 16 females in I. k. moldavica. During the whole period of the study, females were fed with dicotyledonous plants ad libitum. Bioacoustic terminology The bioacoustic terminology used in this study follows Ragge & Reynolds (1998) with some additional terms that we found useful when describing the acoustic signals of phaneropterine bush-crickets (Heller et al., 2004). Calling song: spontaneous song produced by an isolated male or the song produced by a male singing a duet with a female. Generally, this term is used in the literature only for the spontaneous song of an isolated male, but in the Isophya species examined in this paper, males use the same song type when calling alone and when duetting with a female. Syllable: the song produced by one opening–closing movement cycle of the tegmina. Syllable repetition period (SRP): one syllable and the following silent intersyllable interval. Duration of syllable (DS): the time elapsed from the first to the last impulse of the main impulse series of a syllable. Impulse: a simple, undivided, transient train of sound waves (here: the highly damped sound impulse arising as the impact of one tooth of the stridulatory file). Click: an isolated, distinct impulse. Playback stimuli Playback stimuli were made from original male calling song recordings. Male song recordings originated from the same population from which females were collected in the two subspecies of I. kraussii. Song samples of I. camptoxypha males were recorded from populations living in the Bieszczady Mountains (Poland) and in the Mecsek Mountains (Hungary). Male song stimuli (96 kHz sample rate, 16 bit) were modified with minimal acoustic processing: all recordings were high pass filtered (cut-off frequency: 3 kHz) and normalized to 90% of the available 16-bit amplitude range. Sound recordings from five males in both subspecies of I. kraussii and in I. camptoxypha were chosen (altogether, samples of 15 males were included). Each stimulus contained a series of 30 syllables from a single male. Those 15 playback stimuli were arranged in three stimulus series with a randomized stimulus order. In each stimulus series, every stimulus was present in one copy and stimuli were separated from each other by a silent interval of 45 s. Playback protocol Each female was tested separately. Stimulus series composed of original male song samples were played back to the tested female and her acoustic response songs were recorded together with the playback stimuli. To measure the attractiveness of the tested male song samples, the number of syllables answered vs. not answered by the tested female was counted by examining the sound recording of the playback trial (Fig. 2). We applied a no-choice (single-stimulus) playback design. In comparison to two-choice playback, this design gives less-sensitive results regarding the difference of attractiveness of the tested stimuli (Tauber et al., 2001), but provides a more rigorous test of the ability of females to choose a consubspecific/conspecific reproductive partner on the basis of its acoustic signal (Gerhardt & Huber, 2002). Each examined female was placed in a cylinder-shaped mesh cage (diameter 8 cm, height 12 cm). The cage was enclosed in an acoustic isolator box made of 18-mm-thick sheets of laminated wood. The interior of the box was lined with a 5-cm-thick acoustic foam. The box contained a speaker (Alpine SPS-130A) emitting the playback stimuli and a microphone (Brüel & Kjaer, type: 4191) for recording the acoustic response of the female. The microphone was placed between the speaker and the cage. Both the microphone and speaker were directed towards the cage of the tested female from a distance of 5 and 20 cm, respectively (Supporting Information, Fig. S1). Playback stimuli were broadcasted to the experimental animal from a Zoom H4n digital recorder through a Mac Audio (MPX 2000) amplifier driving the speaker [Alpine SPS-130A, frequency response flat (±3 dB) between 150 and 30 000 Hz] within the isolator box. The amplitude of playback stimuli was 69 (±2) dB SPL at the centre of the cage of the experimental animal, based on the root mean square amplitude of the main impulse series of the syllables. The acoustic response of the tested female was recorded using the software Avisoft Recorder 4.2.15 (Avisoft Bioacoustics, Germany) running on a computer connected to a National Instruments data acquisition card (NI USB 6122), which received its input from a 1/2-inch Brüel & Kjaer microphone (type: 4191) through a signal conditioning amplifier (Brüel & Kjaer, Nexus). All three stimulus series were presented to each tested female in a random order so that we could examine whether or not a female response was elicited in the case of 30*15*3 syllables (1350 observations per female) contained by the three stimulus series. We succeeded in performing the playback experiment with 16 females in I. k. kraussii and 16 females in I. k. moldavica. The software Adobe Audition 1.5 (San Jose, CA, USA) was used to count the number of female responses elicited by the playback stimuli. Oscillograms presented in this paper were produced using the R-package ‘seewave’ (Sueur, Aubin & Simonis, 2008). Figure 2. Open in new tabDownload slide Two-traced oscillograms of male–female pair-forming duets: A, Isophya kraussii kraussii; B, I. kraussii moldavica (in order to make the duet structure easier to observe, male song was removed from the female tracks). Figure 2. Open in new tabDownload slide Two-traced oscillograms of male–female pair-forming duets: A, Isophya kraussii kraussii; B, I. kraussii moldavica (in order to make the duet structure easier to observe, male song was removed from the female tracks). Statistical analysis Our main question was whether male song type (consubspecific, heterosubspecific, heterospecific) has an effect on female response probability in the two subspecies of I. kraussii. To examine that question, a mixed-effect logistic regression model was calculated using the lme4 package (Bates et al., 2015) running in the R statistical environment (R Core Team, 2014). A logistic regression model describes the behaviour of a categorical, dependent variable in response to one or more explanatory (independent) variables. In its basic version, the examined dependent variable is binary and can be given as a series of values 0 and 1 or the proportion of those two outcomes (Hilbe, 2009). The dependent variable of our model was the ratio of answered vs. not answered syllables in the playback stimuli. Our observations (whether a female response song occurred) were not independent of each other because of the repeated observations of females, the repeated presentation of the same playback stimuli and the repeated use of the same three randomized stimulus series. Therefore, a mixed-effect model was calculated where random factors can account for the grouped pattern of observations (Faraway, 2016). The calculated mixed-effect logistic regression model included two fixed-effect variables and their interaction: (1) subspecies identity of female and (2) male song type as a factor with three levels such as consubspecific, heterosubspecific and heterospecific. Furthermore, the following random factors were included in the model: specimen identity of the tested females, specimen identity of the males from which the playback stimuli were recorded and the identity of the stimulus series. RESULTS A mixed-effect logistic regression model calculated on the basis of our playback results (Supporting Information, Table S1) showed that females of both subspecies responded to the male song of their own subspecies with significantly higher probability than to heterosubspecific and heterospecific male songs (Table 1; Fig. 3). Response proportion (syllables responded/all presented syllables) for consubspecific songs was similar in the two subspecies (Fig. 3), and the female’s subspecies identity had no significant effect on the probability of response to own-subspecies male song (Table 1). No significant interaction between female subspecific identity and the effect of heterosubspecific male song on response probability was detected (Table 1). In contrast, interaction between female’s subspecies identity and the effect of heterosubspecific male song on response probability was significant: I. k. kraussii females discriminated against heterospecific (I. camptoxypha) songs more strongly than I. k. moldavica females did. Figure 3. Open in new tabDownload slide Variation of female response ratio measured during single-stimulus playback experiments in the females of the two subspecies of Isophya kraussii in response to consubspecific, heterosubspecific and heterospecific (I. camptoxypha) male song samples. Statistical significance of differences (NS, not significant, ***P < 0.001) is indicated based on the results of a mixed-effect logistic regression model calculated on the basis of the playback results (Table 1). The figure is based on the playback trials of 32 females (16 females in each subspecies). Figure 3. Open in new tabDownload slide Variation of female response ratio measured during single-stimulus playback experiments in the females of the two subspecies of Isophya kraussii in response to consubspecific, heterosubspecific and heterospecific (I. camptoxypha) male song samples. Statistical significance of differences (NS, not significant, ***P < 0.001) is indicated based on the results of a mixed-effect logistic regression model calculated on the basis of the playback results (Table 1). The figure is based on the playback trials of 32 females (16 females in each subspecies). Table 1. Fixed-effect estimates of a mixed-effect logistic regression model testing the effect of male song taxonomic identity on the probability of acoustic response in females of the two subspecies of Isophya kraussii Fixed effect . Estimate . SE . z . P . Intercept 0.217 0.295 0.733 0.464 Male song is heterosubspecific (reference level: male song is consubspecific) −2.118 0.274 −7.741 <0.0001 Male song is heterospecific (ref. consubspecific male song) −4.576 0.286 −16.027 <0.0001 Female’s subspecies identity: I. kraussii moldavica (ref. I. kraussii kraussii) 0.075 0.417 0.18 0.857 Interaction: male song heterosubspecific * female’s subspecies identity is I. k. moldavica 0.241 0.543 0.443 0.657 Interaction: male song heterospecific * female’s subspecies identity is I. k. moldavica 3.456 0.288 12.001 <0.0001 Fixed effect . Estimate . SE . z . P . Intercept 0.217 0.295 0.733 0.464 Male song is heterosubspecific (reference level: male song is consubspecific) −2.118 0.274 −7.741 <0.0001 Male song is heterospecific (ref. consubspecific male song) −4.576 0.286 −16.027 <0.0001 Female’s subspecies identity: I. kraussii moldavica (ref. I. kraussii kraussii) 0.075 0.417 0.18 0.857 Interaction: male song heterosubspecific * female’s subspecies identity is I. k. moldavica 0.241 0.543 0.443 0.657 Interaction: male song heterospecific * female’s subspecies identity is I. k. moldavica 3.456 0.288 12.001 <0.0001 Model estimates are given in logits. Model formula: probability of female acoustic response ~ male song type * female subspecific identity + (1|female_id) + (1|male song sample_id) + (1|male song sample series_id). Nobservations = 1440, Nfemales = 32, Nmale song samples = 15, Nmale song sample series = 3; male song type is a factor with three levels: consubspecific, heterosubspecific and heterospecific; one observation means a data pair: the number of answered vs. not answered syllables in case of a stimulus containing 30 syllables. Open in new tab Table 1. Fixed-effect estimates of a mixed-effect logistic regression model testing the effect of male song taxonomic identity on the probability of acoustic response in females of the two subspecies of Isophya kraussii Fixed effect . Estimate . SE . z . P . Intercept 0.217 0.295 0.733 0.464 Male song is heterosubspecific (reference level: male song is consubspecific) −2.118 0.274 −7.741 <0.0001 Male song is heterospecific (ref. consubspecific male song) −4.576 0.286 −16.027 <0.0001 Female’s subspecies identity: I. kraussii moldavica (ref. I. kraussii kraussii) 0.075 0.417 0.18 0.857 Interaction: male song heterosubspecific * female’s subspecies identity is I. k. moldavica 0.241 0.543 0.443 0.657 Interaction: male song heterospecific * female’s subspecies identity is I. k. moldavica 3.456 0.288 12.001 <0.0001 Fixed effect . Estimate . SE . z . P . Intercept 0.217 0.295 0.733 0.464 Male song is heterosubspecific (reference level: male song is consubspecific) −2.118 0.274 −7.741 <0.0001 Male song is heterospecific (ref. consubspecific male song) −4.576 0.286 −16.027 <0.0001 Female’s subspecies identity: I. kraussii moldavica (ref. I. kraussii kraussii) 0.075 0.417 0.18 0.857 Interaction: male song heterosubspecific * female’s subspecies identity is I. k. moldavica 0.241 0.543 0.443 0.657 Interaction: male song heterospecific * female’s subspecies identity is I. k. moldavica 3.456 0.288 12.001 <0.0001 Model estimates are given in logits. Model formula: probability of female acoustic response ~ male song type * female subspecific identity + (1|female_id) + (1|male song sample_id) + (1|male song sample series_id). Nobservations = 1440, Nfemales = 32, Nmale song samples = 15, Nmale song sample series = 3; male song type is a factor with three levels: consubspecific, heterosubspecific and heterospecific; one observation means a data pair: the number of answered vs. not answered syllables in case of a stimulus containing 30 syllables. Open in new tab DISCUSSION Females in both subspecies of I. kraussii responded with significantly higher probability to consubspecific male calling songs than to heterosubspecific signals. However, with lower probability, heterosubspecific male songs were also effective in eliciting female response (Fig. 3; Table 1). These findings suggest that the divergence of male–female acoustic communication would cause an assortative mating pattern in the case of secondary contact of the two taxa, but acoustic divergence alone could not maintain a complete reproductive isolation between them. Therefore, our results are in accordance with the present subspecific treatment of the two taxa. Sometimes, male signals vary significantly between populations, but their conspicuous differences are not important to females. For example, in the case of the frog Pleurodema thaul, a remarkably high, interpopulation call divergence was detected (Velásquez et al., 2013). That divergence seems to be unimportant for females and is likely to be the outcome of inter-male competition (Velásquez et al., 2015). New features of male signals may evolve whether or not those are important to their females, as in the bush-crickets Neoconocephalus retusus and N. maxillosus, where females show no preferences for the derived, double-pulsed calls of the conspecific males (Bush & Schul, 2010). Another possibility is that the variation in male signals significantly affects their attractiveness, but female preferences mismatch male signal variation. For example, in the frog species group Physalaemus pustulosus, females show preferences for signal traits not produced by their males, but those are present in some of the closely related species in the group, suggesting that the evolution of male signals is driven by pre-existing female sensory biases (Ryan & Rand, 1993). Our study shows the most frequently reported pattern, where variation in male signals matches the variation of female preferences when closely related taxa are compared (e.g. Henry, 1985; Barth & Schmitt, 1991; Rodríguez, Ramaswamy & Cocroft, 2006; Boumans & Johnsen, 2014). Our results provide only a snapshot of the evolution of male–female acoustic communication in the two subspecies of I. kraussii, and therefore, we cannot make well-supported inferences about the evolutionary history of signaller–receiver divergence in this taxon pair; nevertheless, the supported concordance between female preferences and male signal pattern is most likely to result from tight, gradual coevolution of male song and female preferences. Reproductive isolation between closely related allopatric or parapatric populations is often asymmetric, which may help us to make inferences about the sequence of ancestry (Kaneshiro, 1980), the progress of the speciation process (Arnold, Verrell & Tilley, 1996) and gene flow during hybridization (Hochkirch & Lemke, 2011). Our results suggest that reproductive isolation associated with acoustic communication is not strongly asymmetric in the case of this pair of subspecies (Tab. 1; Fig. 3). A quantitative genetic model (Arnold et al., 1996) predicts that high levels of asymmetry in reproductive isolation are most likely to occur during the intermediate phase of speciation, suggesting that the two subspecies examined here are at an initial or at a highly progressed phase of speciation. Morphology of the stridulatory files suggests that the latter may be the case, since the stridulatory files show clear differences in the two subspecies (Iorgu & Heller, 2013), indicating that behavioural divergence [which often appears to be faster than morphological differentiation (Henry, 1994; Jones, 1997)] is accompanied by the divergence of stridulatory organs. Location-specific interaction with sympatric species is a potential driving force of speciation. Concerning the evolution of mate recognition, the presence or absence of closely related species is of special importance (Hamao, 2016). We found that females of the two subspecies of I. kraussii showed significantly different levels of discrimination against the song of I. camptoxypha, a closely related species (Fig. 3). Isophya camptoxypha is sympatric with I. k. kraussii, but allopatric with I. k. moldavica. As expected, I. k. kraussii females discriminated more strongly against the male song of I. camptoxypha than females of I. k. moldavica did. A possible explanation for this is that females are not exposed to a selection pressure favouring discrimination against a signal type of an allopatric relative, since they have no chance to meet and hybridize. Similarly, weaker discrimination against allopatric relatives than against sympatric ones has been reported in a wide range of animals (Barth & Schmitt, 1991; Honda-Sumi, 2005; Braune, Schmidt & Zimmermann, 2008; Bewick & Dyer, 2014). Another possible explanation, not excluding the previous one, for the weak discrimination of I. k. moldavica females against I. camptoxypha song is that the male signals of the two taxa are similar to each other (see the oscillograms in Fig. 1). The male songs of I. k. kraussii and I. camptoxypha differ conspicuously in the DS, and I. k. moldavica is intermediate regarding this signal parameter. Therefore, even with the same accuracy of signal recognition, I. k. moldavica females may be more responsive to I. camptoxypha songs than I. k. kraussii females are. Our study is one step towards a better understanding of the evolution and behaviour of these insects. A number of interesting questions arise and await for further examination. Which signal traits are important and which ones are unimportant regarding mate recognition or sexual selection in the two subspecies? The chemical communication channel may also be involved in mate recognition, as in other orthopterans (e.g. Ritchie, 1990; Tyler et al., 2015; Finck & Ronacher, 2017); furthermore, postzygotic isolation may develop simultaneously with prezygotic isolation in allopatric closely related taxa (Gray et al., 2016). Geographic patterns in acoustic signalling and molecular genetic divergence (Pecsenye, Vadkerti & Varga, 2003) could provide us information needed to answer further interesting questions regarding the evolution of these subspecies. For example, is the larger male song difference between I. k. kraussii and I. camptoxypha a result of character displacement? We hope that our study draws attention to this species, which may be a suitable subject for further research on speciation, sexual selection and the relationship between those two processes. SUPPORTING INFORMATION Additional Supporting Information may be found in the online version of this article at the publisher's web-site: Figure S1. A schematic drawing showing the configuration of playback devices within the isolator box from above. S, speaker; M, microphone; A, cage of the tested bush-cricket female; arrows show the direction of speaker and microphone. Table S1. The results of playback experiments examining the response rate of females to consubspecific, heterosubspecific and heterospecific male songs in the two subspecies of Isophya kraussii. ACKNOWLEDGEMENTS We thank three anonymous reviewers for their helpful comments. We are grateful to L. M. E. Sutcliffe for checking the English of an earlier version of the manuscript. Our sincere thanks are due to C. S. Henry for his valuable comments on the manuscript and for linguistic corrections of the text. This work was supported by a grant from the Hungarian National Research Fund (OTKA/NKFI K81929) and by a grant from the Romanian National Authority for Scientific Research and Innovation, CNCS-UEFISCDI, project number PN-II-RU-TE-2014-4-2093. REFERENCES Arnold SJ , Verrell PA , Tilley SG . 1996 . The evolution of asymmetry in sexual isolation: a model and a test case . Evolution 50 : 1024 – 1033 . Google Scholar Crossref Search ADS WorldCat Ballentine B . 2006 . Morphological adaptation influences the evolution of a mating signal . Evolution 60 : 1936 – 1944 . Google Scholar Crossref Search ADS WorldCat Barbosa F , Rebar D , Greenfield MD . 2016 . Female preference functions drive interpopulation divergence in male signalling: call diversity in the bushcricket Ephippiger diurnus . Journal of Evolutionary Biology 29 : 2219 – 2228 . Google Scholar Crossref Search ADS WorldCat Barth FG , Schmitt A . 1991 . Species recognition and species isolation in wandering spiders (Cupiennius spp.; Ctenidae) . Behavioral Ecology and Sociobiology 29 : 333 – 339 . Google Scholar Crossref Search ADS WorldCat Bates D , Mächler M , Bolker B , Walker S . 2015 . Fitting linear mixed-effects models using lme4 . Journal of Statistical Software 67 : 1 – 48 . Bauer N , Kenyeres Z . 2006 . Habitat preference studies of some species of the genus Isophya Brunner von Wattenwyl, 1878 (Orthoptera: Phaneropteridae) in the western part of the Carpathian Basin . Journal of Orthoptera Research 15 : 175 – 185 . Google Scholar Crossref Search ADS WorldCat Bewick ER , Dyer KA . 2014 . Reinforcement shapes clines in female mate discrimination in Drosophila subquinaria . Evolution 68 : 3082 – 3094 . Google Scholar Crossref Search ADS WorldCat Boumans L , Johnsen A . 2014 . Species-specific communication bars interspecific mating between syntopic species of Zwicknia stoneflies (Plecoptera: Capniidae) . Biological Journal of the Linnean Society 113 : 969 – 980 . Google Scholar Crossref Search ADS WorldCat Braune P , Schmidt S , Zimmermann E . 2008 . Acoustic divergence in the communication of cryptic species of nocturnal primates (Microcebus ssp.) . BMC Biology 6 : 19 . Google Scholar Crossref Search ADS WorldCat Bush SL , Schul J . 2010 . Evolution of novel signal traits in the absence of female preferences in Neoconocephalus katydids (Orthoptera, Tettigoniidae) . PLoS ONE 5 : e12457 . Google Scholar Crossref Search ADS WorldCat Campbell P , Pasch B , Pino JL , Crino OL , Phillips M , Phelps SM . 2010 . Geographic variation in the songs of Neotropical singing mice: testing the relative importance of drift and local adaptation . Evolution 64 : 1955 – 1972 . OpenURL Placeholder Text WorldCat Chobanov DP , Grzywacz B , Iorgu IŞ , Ciplak B , Ilieva MB , Warchałowska-Śliwa E . 2013 . Review of the Balkan Isophya (Orthoptera: Phaneropteridae) with particular emphasis on the Isophya modesta group and remarks on the systematics of the genus based on morphological and acoustic data . Zootaxa 3658 : 1 – 81 . Google Scholar Crossref Search ADS WorldCat Cigliano MM , Braun H , Eades DC , Otte D . 2017 . Orthoptera species file, Version 5.0/5.0 . Available at: http://Orthoptera.SpeciesFile.org Cocroft RB , Rodriguez RL , Hunt RE . 2010 . Host shifts and signal divergence: mating signals covary with host use in a complex of specialized plant-feeding insects . Biological Journal of the Linnean Society 99 : 60 – 72 . Google Scholar Crossref Search ADS WorldCat de Oliveira Gordinho L , Matheu E , Hasselquist D , Neto JM . 2015 . Song divergence between subspecies of reed bunting is more pronounced in singing styles under sexual selection . Animal Behaviour 107 : 221 – 231 . Google Scholar Crossref Search ADS WorldCat Faraway JJ . 2016 . Extending the linear model with R: generalized linear, mixed effects and nonparametric regression models, 2nd edn . Boca Raton : CRC Press . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Finck J , Ronacher B . 2017 . Components of reproductive isolation between the closely related grasshopper species Chorthippus biguttulus and C. mollis . Behavioral Ecology and Sociobiology 71 : 70 . Google Scholar Crossref Search ADS WorldCat Gerhardt HC . 2013 . Geographic variation in acoustic communication: reproductive character displacement and speciation . Evolutionary Ecology Research 15 : 605 – 632 . OpenURL Placeholder Text WorldCat Gerhardt HC , Huber F . 2002 . Acoustic communication in insects and anurans: common problems and diverse solutions . Chicago : University of Chicago Press . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Goodman KR , Kelley JP , Welter SC , Roderick GK , Elias DO . 2015 . Rapid diversification of sexual signals in Hawaiian Nesosydne planthoppers (Hemiptera: Delphacidae): the relative role of neutral and selective forces . Journal of Evolutionary Biology 28 : 415 – 427 . Google Scholar Crossref Search ADS WorldCat Grace JL , Shaw KL . 2011 . Coevolution of male mating signal and female preference during early lineage divergence of the Hawaiian cricket, Laupala cerasina . Evolution 65 : 2184 – 2196 . Google Scholar Crossref Search ADS WorldCat Gray DA , Cade WH . 2000 . Sexual selection and speciation in field crickets . Proceedings of the National Academy of Sciences 97 : 14449 – 14454 . Google Scholar Crossref Search ADS WorldCat Gray DA , Gutierrez NJ , Chen TL , Gonzalez C , Weissman DB , Cole JA . 2016 . Species divergence in field crickets: genetics, song, ecomorphology, and pre- and postzygotic isolation . Biological Journal of the Linnean Society 117 : 192 – 205 . Google Scholar Crossref Search ADS WorldCat Guerra PA , Morris GK . 2002 . Calling communication in meadow katydids (Orthoptera, Tettigoniidae): female preferences for species-specific wingstroke rates . Behaviour 139 : 23 – 43 . Google Scholar Crossref Search ADS WorldCat Hamao S . 2016 . Asymmetric response to song dialects among bird populations: the effect of sympatric related species . Animal Behaviour 119 : 143 – 150 . Google Scholar Crossref Search ADS WorldCat Heller K-G . 1990 . Evolution of song pattern in East Mediterranean Phaneropterinae . In: Bailey WJ , Rentz DCF , eds. The Tettigoniidae: biology, systematics, evolution . Bathurst : Crawford House Press , 130 – 151 . Google Scholar Crossref Search ADS Google Scholar Google Preview WorldCat COPAC Heller K-G , Ingrisch S , Liu CX , Shi FM , Hemp C , Warchałowska-Śliwa E , Rentz DC . 2017 . Complex songs and cryptic ethospecies: the case of the Ducetia japonica group (Orthoptera: Tettigonioidea: Phaneropteridae: Phaneropterinae) . Zoological Journal of the Linnean Society .doi: 10.1093/zoolinnean/zlw019 . OpenURL Placeholder Text WorldCat Heller K-G , Korsunovskaya O , Ragge DR , Vedenina V , Willemse F , Zhantiev RD , Frantsevich L . 1998 . Check-list of European Orthoptera . Articulata 7 : 1 – 61 . OpenURL Placeholder Text WorldCat Heller K-G , Orci KM , Grein G , Ingrisch S . 2004 . The Isophya species of Central and Western Europe (Orthoptera: Tettigonioidea: Phaneropteridae) . Tijdschrift voor Entomologie 147 : 237 – 258 . Google Scholar Crossref Search ADS WorldCat Heller K-G , von Helversen O , Sergejeva M . 1997 . Indiscriminate response behaviour in a female buschcricket: sex role reversal in selectivity of acoustic mate recognition? Naturwissenschaften 84 : 252 – 255 . Google Scholar Crossref Search ADS WorldCat Henry CS . 1985 . Sibling species, call differences, and speciation in green lacewings (Neuroptera: Chrysopidae: Chrysoperla) . Evolution 39 : 965 – 984 . OpenURL Placeholder Text WorldCat Henry CS . 1994 . Singing and cryptic speciation insects . Trends in Ecology & Evolution 9 : 388 – 392 . Google Scholar Crossref Search ADS WorldCat Hilbe JM . 2009 . Logistic regression models . Boca Raton : CRC Press . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Hochkirch A , Lemke I . 2011 . Asymmetric mate choice, hybridization, and hybrid fitness in two sympatric grasshopper species . Behavioral Ecology and Sociobiology 65 : 1637 – 1645 . Google Scholar Crossref Search ADS WorldCat Honda-Sumi E . 2005 . Difference in calling song of three field crickets of the genus Teleogryllus: the role in premating isolation . Animal Behaviour 69 : 881 – 889 . Google Scholar Crossref Search ADS WorldCat Iorgu IŞ , Heller KG . 2013 . The bush-cricket Isophya kraussii (Orthoptera: Phaneropteridae): bioacoustics, distribution and description of a new subspecies from Romania . Zootaxa 3640 : 258 – 269 . Google Scholar Crossref Search ADS WorldCat Jones G . 1997 . Acoustic signals and speciation: the roles of natural and sexual selection in the evolution of cryptic species . Advances in the Study of Behaviour 26 : 317 – 354 . Google Scholar Crossref Search ADS WorldCat Kaneshiro KY . 1980 . Sexual isolation, speciation and the direction of evolution . Evolution 34 : 437 – 444 . Google Scholar Crossref Search ADS WorldCat Kenyeres Z , Bauer N . 2005 . Untersuchung des Lebensraumes von Isophya camptoxypha (Fieber, 1853) im Kőszeg Gebirge (Westungarn) . Articulata 20 : 1 – 15 . OpenURL Placeholder Text WorldCat Kowalski K , Lakes-Harlan R . 2011 . Temporal patterns of intra-and interspecific acoustic signals differ in two closely related species of Acanthoplus (Orthoptera: Tettigoniidae: Hetrodinae) . Zoology 114 : 29 – 35 . Google Scholar Crossref Search ADS WorldCat Lachlan RF , Servedio MR . 2004 . Song learning accelerates allopatric speciation . Evolution 58 : 2049 – 2063 . Google Scholar Crossref Search ADS WorldCat Linossier J , Zsebők S , Baudry E , Aubin T , Courvoisier H . 2016 . Acoustic but no genetic divergence in migratory and sedentary populations of blackcaps, Sylvia atricapilla . Biological Journal of the Linnean Society 119 : 68 – 79 . Google Scholar Crossref Search ADS WorldCat Ng EY , Eaton JA , Verbelen P , Hutchinson RO , Rheindt FE . 2016 . Using bioacoustic data to test species limits in an Indo-Pacific island radiation of Macropygia cuckoo doves . Biological Journal of the Linnean Society 118 : 786 – 812 . Google Scholar Crossref Search ADS WorldCat Orci KM . 2007 . Female preferences for male song characters in the bush-cricket Isophya camptoxypha (Orthoptera, Tettigonioidea) . Journal of Insect Behavior 20 : 503 – 513 . Google Scholar Crossref Search ADS WorldCat Pecsenye K , Vadkerti E , Varga Z . 2003 . Temporal and spatial pattern of genetic differentiation in Isophya kraussii (Orthoptera: Tettigonoidea) in NE Hungary . Acta Zoologica Academiae Scientiarum Hungaricae 49 : 167 – 178 . OpenURL Placeholder Text WorldCat Perdeck AC . 1958 . The isolating value of specific song patterns in two sibling species of grasshoppers (Chorthippus brunneus Thunb. and C. biguttulus L.) . Behaviour 12 : 1 – 75 . Google Scholar Crossref Search ADS WorldCat Ragge DR , Reynolds WJ . 1998 . The songs of the grasshoppers and crickets of Western Europe . Colchester : Harley Books . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC R Core Team . 2014 . R: a language and environment for statistical computing . Vienna : R Foundation for Statistical Computing . Available at: http://www.R-project.org/ Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Ritchie MG . 1990 . Are differences in song responsible for assortative mating between subspecies of the grasshopper Chorthippus parallelus (Orthoptera: Acrididae)? Animal Behaviour 39 : 685 – 691 . Google Scholar Crossref Search ADS WorldCat Robinson DJ , Hall MJ . 2002 . Sound signalling in Orthoptera . Advances in Insect Physiology 29 : 151 – 278 . Google Scholar Crossref Search ADS WorldCat Rodríguez RL , Ramaswamy K , Cocroft RB . 2006 . Evidence that female preferences have shaped male signal evolution in a clade of specialized plant-feeding insects . Proceedings of the Royal Society of London B: Biological Sciences 273 : 2585 – 2593 . Google Scholar Crossref Search ADS WorldCat Ryan MJ , Rand AS . 1993 . Sexual selection and signal evolution: the ghost of biases past . Philosophical Transactions of the Royal Society of London B: Biological Sciences 340 : 187 – 195 . Google Scholar Crossref Search ADS WorldCat Schul J , von Helversen D , Weber T . 1998 . Selective phonotaxis in Tettigonia cantans and T. viridissima in song recognition and discrimination . Journal of Comparative Physiology A 182 : 687 – 694 . Google Scholar Crossref Search ADS WorldCat Simmons LW , Zuk M , Rotenberry JT . 2001 . Geographic variation in female preference functions and male songs of the field cricket Teleogryllus oceanicus . Evolution 55 : 1386 – 1394 . Google Scholar Crossref Search ADS WorldCat Stumpner A , von Helversen O . 1994 . Song production and song recognition in a group of sibling grasshopper species (Chorthippus dorsatus, Ch. dichrous and Ch. loratus: Orthoptera, Acrididae) . Bioacoustics 6 : 1 – 23 . Google Scholar Crossref Search ADS WorldCat Sueur J , Aubin T , Simonis C . 2008 . Seewave, a free modular tool for sound analysis and synthesis . Bioacoustics 18 : 213 – 226 . Google Scholar Crossref Search ADS WorldCat Tauber E , Cohen D , Greenfield MD , Pener MP . 2001 . Duet singing and female choice in the bushcricket Phaneroptera nana . Behaviour 138 : 411 – 430 . Google Scholar Crossref Search ADS WorldCat Tyler F , Fisher D , d’Ettorre P , Rodríguez-Muñoz R , Tregenza T . 2015 . Chemical cues mediate species recognition in field crickets . Frontiers in Ecology and Evolution 3 : 48 . Google Scholar Crossref Search ADS WorldCat Velásquez NA , Marambio J , Brunetti E , Méndez MA , Vásquez RA , Penna M . 2013 . Bioacoustic and genetic divergence in a frog with a wide geographical distribution . Biological Journal of the Linnean Society 110 : 142 – 155 . Google Scholar Crossref Search ADS WorldCat Velásquez NA , Valdés JL , Vásquez RA , Penna M . 2015 . Lack of phonotactic preferences of female frogs and its consequences for signal evolution . Behavioural Processes 118 : 76 – 84 . Google Scholar Crossref Search ADS WorldCat Wilkins MR , Seddon N , Safran RJ . 2013 . Evolutionary divergence in acoustic signals: causes and consequences . Trends in Ecology & Evolution 28 : 156 – 166 . Google Scholar Crossref Search ADS WorldCat Zhantiev RD , Dubrovin NN . 1977 . Sound communication in the grasshopper genus Isophya (Orthoptera, Tettigoniidae) . Zoologicheskii Zhurnal 56 : 38 – 51 (in Russian, with English summary). OpenURL Placeholder Text WorldCat Zhantiev RD , Korsunovskaya OS . 1990 . Sound communication of Phaneropteridae (Orthoptera) . In: Gribakin FG , Wiese K , Popov AV , eds. Sensory systems and communication in arthropods . Basel : Birkhäuser , 402 – 406 . Google Scholar Crossref Search ADS Google Scholar Google Preview WorldCat COPAC Zhantiev RD , Korsunovskaya OS , Benediktov A . 2017 . Acoustic signals of the bush-crickets Isophya (Orthoptera: Phaneropteridae) from Eastern Europe, Caucasus and adjacent territories . European Journal of Entomology 114 : 301 – 3011 . Google Scholar Crossref Search ADS WorldCat Zuk M , Simmons LW , Cupp L . 1993 . Calling characteristics of parasitized and unparasitized populations of the field cricket Teleogryllus oceanicus . Behavioral Ecology and Sociobiology 33 : 339 – 343 . OpenURL Placeholder Text WorldCat © 2017 The Linnean Society of London, Biological Journal of the Linnean Society TI - Subspecies-specific song preferences and diverged heterospecific discrimination in females of the bush-cricket Isophya kraussii (Orthoptera: Phaneropterinae) JO - Biological Journal of the Linnean Society DO - 10.1093/biolinnean/blx094 DA - 2017-10-25 UR - https://www.deepdyve.com/lp/oxford-university-press/subspecies-specific-song-preferences-and-diverged-heterospecific-SxfD7XAVdf SP - 505 VL - 122 IS - 3 DP - DeepDyve ER -