Fedorova - Shabanov - 2022 - Differences in Release Calls
Fedorova - Shabanov - 2022 - Differences in Release Calls
Fedorova - Shabanov - 2022 - Differences in Release Calls
https://doi.org/10.1007/s11756-022-01256-8
ORIGINAL ARTICLE
Abstract
Vocalization in anurans is mainly studied in the context of male-specific advertisement calls, while release calls – charac-
teristic of both sexes – usually attract less attention. Here we compare release calls of both sexes of hybrid Pelophylax escu-
lentus (considering triploid hybrids) and its parental species P. ridibundus in Ukrainian population systems. The principal
component analysis resulted in 83.3% of total variance explained by the first three components, while discriminant analysis
correctly differentiated 69.6% of frogs only. Triploid male hybrids had intermediate positions between diploid hybrids
and parental species. Triploid females had intermediate positions only in terms of temporal parameters because frequency
parameters mainly depended on their body size. A significant correlation between peak frequency and body size was shown
for females of both taxonomical groups and P. ridibundus males. Although we identified differences between the studied
groups, these differences were not unambiguous, with overlapping frequency parameters in males. Including more triploids
and the other parental species is important for further analysis.
13
Vol.:(0123456789)
Biologia
Here we compare the release calls of the hybrid P. escu- Ploidy identification
lentus with its parental species P. ridibundus of both sexes
in a Ukrainian population to find specific traits that can help To identify ploidy, we measured erythrocyte lengths for
distinguish between hybrids and their parental species. It was all studied frogs. We sampled blood either from the facial
shown that advertisement call characteristics of Pelophylax vein (Forzán et al. 2012) or cut fingertips while marking
frogs could depend on geography and types of population individuals for annual monitoring and studying population
system (Wycherley et al. 2002; Hoffmann and Reyer 2013). size. We prepared an air-dried blood smear for each frog
Since our research was exclusively carried out in the Siver- and measured lengths of at least 30 erythrocytes under × 40
skyi Donets River basin (Eastern Ukraine, Kharkiv region), magnification, either with a Carl Zeiss Amplival or Leica
where the other parental species P. lessonae is absent (Biriuk DM 2000 microscope. Bondarieva et al. (2012) identified
et al. 2016), we did not include P. lessonae in our analysis. the margin between diploid and triploid erythrocytes of
hybrids in the Siversky Donets River as 26 µm. However,
our further studies showed that this margin differs between
Materials and methods population systems of this river basin, and attention should
be paid to the “gap” in cell lengths measured for all the
Sample collection individuals in a sample (Drohvalenko et al. 2019).
Some frogs had erythrocyte lengths close to the mar-
Frogs were collected at four localities in the National Park gin or “gap”. For those individuals, we performed kar-
“Homilshanski lisy” in Ukraine (Table 1). Frogs were col- yoanalysis. Frogs were injected with 0.05% colchicine
lected from 2014 to 2021 during the breeding season (spring solution for 12–18 h and then sacrificed with vapors
and summer). Taxonomical and sex identification was car- of ETAC or injection of 2% lidocaine. Further, we dis-
ried out by examining external morphology (Atemasova sected intestine tissue and bone marrow, put it into a
et al. 2019). We examined the shape of the metatarsal tuber- hypotonic solution of 0.07 M KCl for 30 min, centri-
cule, coloration of inguinal region, underside of thighs and fuged at 3000 rpm, and replaced the hypotonic solution
vocal sacs, and spotting pattern. While taxonomical identifi- with Carnoy’s fixative (3 parts of ethanol or methanol
cation is relatively easy for diploid frogs, it could be errone- and 1 part of glacial acetic acid). Small pieces (roughly
ous for triploid hybrids. Therefore, we also identified ploidy 5 m m 2) of intestine tissue or bone marrow were homog-
for all the frogs, even those identified as parental species. enized in 100 µl of 70% acetic acid. The resulting cell
For sex identification, we examined the presence of vocal suspension was put on slides heated to 60° C and left to
sacs and nuptial pads. dry. Slides were stained with A gNO 3 following Birstein
Frogs were fed and kept in 20–50 L plastic tanks for (1984) with modifications and analyzed under × 100
1–2 days. Recordings were performed under laboratory con- magnification, either with a Carl Zeiss Amplival or
ditions at ~ 25 °C. Calls were recorded with Nikon Coolpix Leica DM 2000 microscope.
P500 and Nikon D3100 cameras placed at a distance of Pelophylax species used in this study are not listed in
20 cm from frogs. Release calls were stimulated by gently the IUCN Red List of Threatened Species or Red Data
pressing the frog flanks to simulate amplexus. Since it was Book of Ukraine. Tissue collection was performed to mini-
shown that call traits could depend on the body size (e.g., mize animal harm according to the Guidelines for Use of
Davies and Halliday 1978; Gingras et al. 2013), ninety-one Live Amphibians and Reptiles in Field and Laboratory
frogs were measured for snout-vent length (SVL) with a Research of HACC and The Ukrainian Law on the Protec-
vernier caliper. tion of Animals from Cruelty.
13
Biologia
Fig. 1 Oscillograms and spectrograms of release calls of Pelophy- tus triploid females, Em 3n – P. esculentus triploid males, Rf – P.
lax males (a) and females (b): Ef 2n – P. esculentus diploid ridibundus females, Rm – P. ridibundus males
females, Em 2n – P. esculentus diploid males, Ef 3n – P. esculen-
13
Biologia
comparisons between studied groups of frogs. Pearson cor- Fig. 2 Results of principal component analysis: Ef 2n – P. esculentus
diploid females, Em 2n – P. esculentus diploid males, Ef 3n – P. escu-
relation was applied to see the connection between frequency
lentus triploid females, Em 3n – P. esculentus triploid males, Rf – P.
parameters and body size. ridibundus females, Rm – P. ridibundus males
Results from the main part of the call by a pause. Such low amplitude
peaks were not observed in P. esculentus female calls.
We analyzed 732 release calls from 183 frogs, among which Principal component (PCA) analysis resulted in three
110 were P. esculentus (79 males and 31 females), and 73 principal components with eigenvalues > 1.0 and 83.3% of
were P. ridibundus (44 males and 29 females). Among the total variance explained (Table 2).
hybrids, eight frogs were identified as triploids. The general The discriminant analysis using all seven temporal and fre-
structure of the release calls differed between the two taxo- quency parameters resulted in two discriminant functions that
nomical groups. Pelophylax esculentus calls of both sexes were significant (peak frequency and peak rates included into
consisted of single or double peaks with approximately the the model). Only 69.6% of frogs were classified correctly, and
same period of silence between them (Fig. 1a, b). Such sepa- all triploid hybrids were wrongly classified (Table 3), which
rate peaks were not observed in P. ridibundus release calls. corresponds to the results of the PCA (Fig. 2).
In most P. ridibundus male calls, groups of peaks could be Higher values of peak numbers and peak rates were typi-
distinguished, as well as in female calls. The number of peaks cal for P. ridibundus of both sexes, while triploid hybrids
was higher for P. esculentus females and lower for males had an intermediate position between diploid P. esculentus
(Fig. 1). Similar observations were made for P. ridibun- and P. ridibundus (Fig. 3). Frequency parameters differed
dus calls: the number of peaks was significantly higher for only between females, with P. ridibundus having higher
females than males. Most P. esculentus male release calls values (Fig. 3). Post hoc analysis with the Kruskal–Wallis
(71%) start with several peaks of low amplitude separated test showed that call duration was the least variable between
Table 3 Results of discriminant analysis based on three temporal and four frequency parameters
Ef 2n 15 0 2 6 0 5 53.57
Ef 3n 1 0 0 0 0 2 0.00
Rf 2 0 23 0 0 5 76.66
Em 2n 7 0 0 67 0 0 91.89
Em 3n 0 0 0 4 0 1 0.00
Rm 2 0 7 13 0 22 50.00
Bold face – correctly classified frogs, Ef2n – P. esculentus diploid females, Em2n – P. esculentus diploid males, Ef3n – P. esculentus triploid
females, Em3n – P. esculentus triploid males, Rf – P. ridibundus females, Rm – P. ridibundus males
13
Biologia
Fig. 3 Means of the three temporal and four frequency parameters: Ef2n – P. esculentus diploid females, Em2n – P. esculentus diploid males,
Ef3n – P. esculentus triploid females, Em3n – P. esculentus triploid males, Rf – P. ridibundus females, Rm – P. ridibundus males
studied groups (Fig. 3). The most variable parameters were A significant correlation between peak frequency and body
peak numbers and peak rates (Fig. 3). size was shown for females (Fig. 4a, b) of both taxonomical
Call durations in males were lower in our study than groups and P. ridibundus males (Fig. 4d), but not for P. escu-
those reported by Zamfirescu (2002): 0.42 ± 0.071 vs. lentus males (Fig. 4c). However, the significance level for P.
0.48 ± 0.119 in P. esculentus (p = 0.0154) and 0.34 ± 0.079 ridibundus was relatively low, which may be the consequence
vs. 0.45 ± 0.051 in P. ridibundus (p < 0.0001). of the low number of individuals included in the analysis. A high
All seven measured parameters were dynamic (evalu- correlation was also shown between body size and minimum fre-
ated according to Gerhardt 1991) with CV values above quency (p < 0.001) in P. esculentus of both sexes and P. ridibun-
15% (Table 4). Only for triploid hybrids, peak frequencies dus females, but not P. ridibundus males (p = 0.3854). No cor-
and minimum frequencies were evaluated as stable param- relation was found between body size and all other parameters.
eters, which might be caused by the low number of sam- Such correlation explains why triploid females differ drastically
ples. In general, females tended to have higher coefficients from other forms by frequency parameters (Fig. 3): all three trip-
of variation, but the difference was not significant (KW- loid females had relatively small body sizes and, consequently,
H(5;42) = 1.7675; p = 0.8803). rather high “peak frequency” and “min frequency” parameters.
13
Biologia
Table 4 Coefficients of variation (CV) of seven temporal and frequency parameters calculated for six groups of Pelophylax frogs
Call duration Peak frequency Frequency range Min frequency Max frequency Peaks number Peaks rate
Bold – CV ≤ 5%. CV = SD*100/mean. Ef 2n – P. esculentus diploid females, Em 2n – P. esculentus diploid males, Ef 3n – P. esculentus triploid
females, Em 3n – P. esculentus triploid males, Rf – P. ridibundus females, Rm – P. ridibundus males
Fig. 4 Correlation between peak frequency and snout-vent length (SVL). a P. esculentus females (dot – diploid, circle – triploid), b P. ridibun-
dus females, c P. esculentus males (dot – diploid, circle – triploid), d P. ridibundus males
13
Biologia
analysis of triploid P. esculentus shows that they have (Fig. 4a). More triploid females with different body sizes are
an intermediate position between P. ridibundus and dip- needed for a comprehensive analysis. Our analysis also did
loid P. esculentus. Still, they were incorrectly differenti- not include the other parental species – P. lessonae – because
ated by the discriminant analysis. The best classification it is absent in the studied territory (Siversky Donets River
worked for P. esculentus males (90.5%) and P. ridibundus basin), but it should be included in further studies.
females (75.9%). Unambiguous differences between the Although this analysis did not reveal any characteristics
release calls of P. ridibundus and P. esculentus should of release calls that can exclusively differentiate between
not be expected since P. esculentus is a hybrid, and its hybrid P. esculentus (both diploid and triploid) and one of
release call may overlap with the parental ones. However, their parental species P. ridibundus, we showed general sex-
for mating calls, a strict difference between hybrids and and species-specific trends in frogs’ calls that may become
parental species was shown (Hoffmann and Reyer 2013). a basis for further analysis.
For the release calls, we see similar tendency, but compar-
Supplementary Information The online version contains supplemen-
isons with the other parental species P. lessonae and more tary material available at https://d oi.o rg/1 0.1 007/s 11756-0 22-0 1256-8.
triploid hybrids are needed. Some temporal parameters
for both parental species were measured by Chernyshov Acknowledgements The authors thank the students at the School of
(2005). Even though Chernyshov did not compare the two Biology (V. N. Karazin Kharkiv National University) who helped catch
and record frogs during the fieldwork. Anna Fedorova is very grate-
species directly, disparities were observed, e.g., a differ- ful to Eleonora Pustovalova, whose support and interest in this work
ence in call durations of P. ridibundus and P. lessonae, helped to complete this manuscript. The authors are also sincerely
which also depended on frog size and water temperature. grateful to the Armed Forces of Ukraine, whose bravery and profes-
Following Chernyshov, we expect to see a more evident sionalism provided a chance to finish the work on this manuscript.
difference between the release calls of P. lessonae and P.
ridibundus in our further analysis. It is possible that the
Declarations
temporal and frequency data we extracted are not sufficient Conflict of interest The authors declare that they have no conflict of
to distinguish between P. esculentus and P. ridibundus of interest.
both sexes. We suggest that more sophisticated compari-
son methods, which will not be based on individual dimen-
sions but the comparison of patterns as a whole, should
be more effective. References
There is also a possibility that release calls are less spe-
cies-specific than mating calls (Hoffmann and Reyer 2013). Atemasova TA, Vlashchenko AS, Goncharov GL, Zinenko OI, Kors-
hunov OV, Tokarsky VA, Shabanov DA, Shandikov GA (2019)
Also, the function of the release call could be different for Educational and field practice in vertebrate zoology. VN Karazin
different sexes. Males use the call as a direct release signal KhNU, Kharkiv
to indicate their sex, so the difference between species might Audacity Team (2015) Audacity®: free audio editor and recorder
not be very important. Such differences should arise from [Computer program] version 211. http://audacity.sourceforge.
net/. Accessed 20 Aug 2015
selection against interspecific mating. But such selection was Biriuk OV, Shabanov DA, Korshunov AV, Borkin LJ, Lada GA,
ineffective in this group of frogs. Pasynkova RA, Rosanov JM, Litvinchuk SN (2016) Gamete
According to Gerhardt (1991) static parameters are more production patterns and mating systems in water frogs of the
taxonomically relevant than dynamic ones. Call duration and hybridogenetic Pelophylax esculentus complex in north-eastern
Ukraine. J Zool Syst Evol Res 54:215–225. https://doi.org/10.
peak rate were identified as static parameters for many dif- 1111/jzs.12132
ferent species (Gerhardt 1991; Castellano and Rosso 2006). Birstein VJ (1984) Localization of NORs in karyotypes of four Rana
However, in our analysis, these parameters seemed to be species. Genetica 64(3):149–154
dynamic, i.e., intragroup variation was higher than inter- Bondarieva AA, Bibik YuS, Samilo SM, Shabanov DA (2012) Eryth-
rocytes cytogenetic characteristics of green frogs from Siversky
group, which was also confirmed by a Kruskal–Wallis test Donets centre of Pelophylax esculentus complex diversity (in Rus-
explaining the absence of differences between the studied sian). The Journal of V.N.Karazin Kharkiv National University.
groups. Ser Biol 15:116–123
The significant correlation between body size and fre- Castellano S, Rosso A (2006) Variation in call temporal properties and
female preferences in Hyla intermedia. Behaviour 143:405–424
quency parameters corresponds with previous data for anu- Center for Conservation Bioacoustics (2019) Raven pro: interactive
rans (e.g., Davies and Halliday 1978; Gingras et al. 2013). sound analysis software (Version 1.6.1) [Computer software].
Such a correlation explains the results of our post hoc analy- Ithaca, NY: The Cornell Lab of Ornithology. http://www.birds.
sis (Fig. 3). Moreover, it explains why triploid females dif- cornell.edu/raven. Accessed 21 March 2020
Chernyshov KI (2005) Variability of advertisement calls and release
fered so extremely from diploid P. esculentus by frequency calls of green frogs in the Moscow oblast’, Russia, In: Ananjeva N,
parameters. All three triploid females had rather small body Tsinenko O (eds) Proceedings of the12th Ordinary General Meeting
sizes and, therefore, relatively high frequency parameters of the Societas Europaea Herpetologica, Saint-Petersburg, pp 20–26
13
Biologia
Davies NB, Halliday TR (1978) Deep croaks and fighting assessment in Mable BK, Bogart JP (1991) Call analysis of triploid hybrids result-
toads Bufo. Nature 274:683–685. https://d oi.o rg/1 0.1 038/274683a0 ing from diploid-tetraploid species crosses of hylid tree frogs.
Drohvalenko M, Pustovalova E, Fedorova A (2019) Pelophylax escu- Bioacoustics 3:111–119. https://doi.org/10.1080/09524622.1991.
lentus complex from Iskiv pond: One more step of long-term 9753167
monitoring. In: Proceedings of the international zoological Roberts JD (2010) Natural hybrid between the frogs Crinia pseud-
conference "Ukraine fauna at the turn of XX-XXI centuries". insignifera and Crinia subinsignifera (Myobatrachidae) from
CПOЛOM, Shatsk, Lviv, pp 8–10 Southwestern Australia defined by allozyme phenotype and call.
Forzán MJ, Vanderstichel RV, Ogbuah CT, Barta JR, Smith TG (2012) J Herpetol 44:654–657
Blood collection from the facial (maxillary) / musculo-cutaneous Wahl M (1969) Untersuchungen zur Bio-Akustik des Wasserfrosches
vein in true frogs (family Ranidae). J Wildl Dis 48(1):176–80. Rana esculenta (L.). Oecologia 3:14–55. https://doi.org/10.1007/
https://doi.org/10.7589/0090-3558-48.1.176 BF00344632
Gerhardt HC (1974) The vocalizations of some hybrid treefrogs: acous- Weber E (1976) Alterations in the release calls of six European Anura
tic and behavioral analyses. Behaviour 49:130–151 (Amphibia) after partial or total extirpation of the vocal cords. Behav
Gerhardt HC (1991) Female mate choice in treefrogs: static and Process 1:197–216. https://doi.org/10.1016/0376-6357(76)90022-X
dynamic acoustic criteria. Anim Behav 42:615–635. https://doi. Wycherley J, Doran S, Beebee TJC (2002) Male advertisement call
org/10.1016/S0003-3472(05)80245-3 characters as phylogeographical indicators in European water
Gerhardt HC, Guttman SI, Karlin AA (1980) Natural hybrids between frogs. Biol J Linn Soc 77:355–365. https://doi.org/10.1046/j.
Hyla cinerea and Hyla gratiosa: Morphology, vocalization and 1095-8312.2002.00104.x
electrophoretic analysis. Copeia 1980:577–584. https://doi.org/ Zamfirescu Ş (2000) Bioacoustic aspects of the release calls of some
10.2307/1444432 anuran species. Analele Ştiinţifice ale Universităţii “Alexandru
Gingras B, Boeckle M, Herbst ST, Fitch WT (2013) Call acoustics Ioan Cuza” din Iaşi, s. Biologie Animală 46:101–106
reflect body size across four clades of anurans. J Zool 289:143– Zamfirescu Ş (2001) The bioacoustic analysis of some water frogs’
150. https://doi.org/10.1111/j.1469-7998.2012.00973.x release calls (Rana ridibunda, Rana kl. esculenta). Analele
Graf JD, Polls Pelaz M (1989) Evolutionary genetics of the Rana escu- Ştiinţifice ale Universităţii “Alexandru Ioan Cuza” din Iaşi, s.
lenta complex. In: Dawley RM, Bogart JP (eds) Evolution and Biologie Animală 47:109–112
Ecology of Unisexual Vertebrates. New York State Museum, New Zamfirescu Ş (2002) Comparison between water frogs (Rana escu-
York, pp 289–301 lenta complex) release calls’. Analele Ştiinţifice ale Universității
Grenat PR, Martino L (2013) The release call as a diagnostic character “Alexandru Ioan Cuza” din Iaşi, s. Biologie Animală 48:182–194
between cryptic related species Odontophrynus cordobae and O.
americanus (Anura: Cycloramphidae). Zootaxa 3635:583–586. Publisher's note Springer Nature remains neutral with regard to
https://doi.org/10.11646/zootaxa.3635.5.8 jurisdictional claims in published maps and institutional affiliations.
Hoffmann A, Reyer H-U (2013) Genomic effects on advertisement
call structure in diploid and triploid hybrid waterfrogs (Anura, Springer Nature or its licensor (e.g. a society or other partner) holds
Pelophylax esculentus). BMC Ecol 13:47. https://d oi.o rg/1 0.1 186/ exclusive rights to this article under a publishing agreement with the
1472-6785-13-47 author(s) or other rightsholder(s); author self-archiving of the accepted
Köhler J, Jansen M, Rodriguez A, Kok PJR, Toledo LF (2017) The use manuscript version of this article is solely governed by the terms of
of bioacoustics in anuran taxonomy: theory, terminology, meth- such publishing agreement and applicable law.
ods and recommendations for best practice. Zootaxa 4251:1–124.
https://doi.org/10.11646/zootaxa.4251.1.1
13