ABSTRACT
The behavioural ecology of ultrasonic-singing katydids is not well understood, and the general bioacoustics, barely known for a few Neotropical Meconematinae, tends to be overlooked for species from Southeast Asia. These include Asiatic species of Phlugidini, commonly known as crystal predatory katydids. One of its genera, Asiophlugis consists of 16 species for which acoustic signals and stridulum anatomy are broadly unknown. These characters can be used to understand species boundaries. Here, we sampled Asiophlugis from five sites in Malay Peninsula and Borneo Island, recorded the acoustic signals of five species plus one subspecies using ultrasound sensitive equipment, and examined their stridulum anatomy. The calling songs of the taxa involved were documented for the first time. We found that the stridulum anatomy (e.g., tooth distributions, tooth length and tooth density) is distinct between species but less so between subspecies. In contrary, songs of different taxa are different based on acoustic parameters (e.g., pulse duration, peak frequency) and descriptive patterns, even between the subspecies. We also did not observe that song signals are more different in sympatry than in allopatry. Whether this can be generalised requires further sampling, highlighting the need for more research on the ultrasonic acoustic communication in Asiatic katydids.
Introduction
In katydids (Insecta: Orthoptera: Tettigonioidea), ~70% of the species investigated so far communicate using calls with ultrasonic components, and these extend to the extreme ultrasonic frequency of 150 kHz (Montealegre-Z et al. Citation2006, Citation2017). Sound production in katydids is achieved by tegmino-tegminal stridulation, where a scraper in the right forewing or tegmina is swept over a serrated vein (the stridulatory file) on the left forewing (Montealegre-Z Citation2009). Since males produce sound to attract conspecific females, calling song parameters can be used to understand species boundaries. The use of acoustics for examining species boundaries is more prominent for species from the Neotropics than that of relatives from the similarly diverse Southeast Asia. This is probably an artefact of the limited sampling and a lack of acoustic experts working in this region, leading to many species which are known only from their original description (Tan et al. Citation2017). Considering that many katydid clades are more congruent with biogeography (Mugleston et al. Citation2018), lacking knowledge on the acoustics of Southeast Asian species impedes our understanding on the evolution (and biomechanics) of acoustic communication.
One group of poorly known katydids from Southeast Asia belongs to the Phlugidini (Meconematinae). These crystal predatory katydids, which comprise about 100 species globally (Cigliano et al. Citation2019), are characterised by large protruding eyes and movable spines on anterior legs. They have diversified more in South America, with as many as 59 species known to science (Cigliano et al. Citation2019). Except for the acoustic behaviour of two species of Phlugis (see Chamorro-Rengifo et al. Citation2014; Chamorro-Rengifo and Braun Citation2016; Sarria-S et al. Citation2017), very little is known about the acoustics and biomechanics of sound production in this group.
This is also the case for their relatives from Southeast Asia, which include species from the genus Asiophlugis Gorochov Citation1998 (). Asiophlugis currently consists of 16 described species (Cigliano et al. Citation2019) but calling song was known for only one species (Helfert and Sänger Citation1998). Species are delimited using morphological characters on the abdominal apex, pronotum and tegmina length (Kevan and Jin Citation1993; Gorochov Citation1998, Citation2012, Citation2013, Citation2019; Gorochov and Tan Citation2011; Tan Citation2011). Unlike relatives from the Neotropics (see Nickle Citation2003, Citation2005; Chamorro-Rengifo and Olivier Citation2017), their stridulatory apparatus are rarely examined and used in taxonomy (but see Helfert and Sänger Citation1998; Gorochov and Tan Citation2011).
To address the lack of information on the bioacoustics of Asiophlugis, we collected specimens from the Malay Peninsula (which includes Peninsular Malaysia and Singapore) and Borneo (i.e., Brunei Darussalam and Sabah), recorded the songs, analysed and described the acoustic signatures, and examined and compared their stridulum anatomy. Then, we attempt to answer the following research question: Can the call structure and stridulum anatomy assist with inferring species identity? We predict that the acoustic parameters in the domains of frequency (e.g., peak frequency) and time (e.g., pulse repetition rates) and/or organisation of stridulatory file (i.e., inter-tooth distance, tooth length and tooth density) can reflect the species and subspecies delimited by morphology since males should produce species-specific acoustic signals and that stridulum anatomy is responsible for sound production. Alternatively, songs are also sexual signals, which tend to be more different in sympatry than in allopatry (Fitzpatrick and Gray Citation2001; Izzo and Gray Citation2004). If this is true, we can expect syntopic species from Singapore to have very different songs than allopatric species from Peninsular Malaysia, Brunei Darussalam and Sabah.
Materials and methods
Collection and husbandry of katydids
Surveys were conducted by MKT in the lowland forests of (1) Singapore; (2) Pulau Tioman, Pahang, Peninsular Malaysia; (3) Ulu Gombak Field Studies Centre, Selangor, Peninsular Malaysia; (4) Andulau Forest Reserve, Belait, Brunei Darussalam; and (5) Rainforest Discovery Centre, Sandakan, Sabah, East Malaysia () between July 2018 and March 2019. Opportunistic collection was conducted at night-time (1930 to 2300 h). The katydids were then housed in standardised insect cages (25 cm in diameter and 33 cm tall) with nylon cover. They were fed with living fruit flies and regularly sprayed to avoid dehydration. Fresh plants, typically Dillenia suffruticosa, were also provided. The katydids typically survive in captivity between two to four weeks.
Acoustic recordings and analysis
All recordings were obtained in laboratory conditions or in biological stations. Calling songs of males in the cages were recorded using a portable ultrasound recorder Echo Meter Touch (Wildlife Acoustics, Inc. USA, Massachusetts), attached to an iPhone, and using a sampling frequency of 256 kHz-samples/s. The Echo Meter Touch 1 (based on Knowles FG sensor) was placed at about 1 m horizontally away from the cage. We observed that the katydid usually (but not always) sing while clinging onto the side of the cage (thus positioning with dorsal surface facing the recording device). Triggered recording was used with the Trigger Minimum Frequency set at 20 kHz. However, this is only a trigger and will not affect lower frequencies to be recorded. As temperature can influence the frequency of calls, a temperature-humidity metre (Smartsensor AR867) was used to record the ambient temperature (°C) and relative humidity (%) at the time of recording. All recordings were done at an ambient temperature between 27.0°C and 29.8°C.
The recorded signals were saved in 12-bit WAV format. All acoustic analyses were done using the open source R package WarbleR version 1.1.14 (Araya‐Salas and Smith‐Vidaurre Citation2017) in the R software version 3.5.1 (R Development Core Team Citation2018). To filter out low or high background noise before performing measurements, the lower and upper limits of a frequency bandpass filter (in kHz) were set at 20 kHz and 120 kHz, respectively, so measurements of the peak frequency and mean frequency were obtained only between the frequency range. The limits are beyond the frequency range of the Asiophlugis songs and will not affect the measurements of the frequencies. The peak frequency (frequency with highest energy from the mean spectrum, in kHz) and the mean frequency (weighted average of frequency by amplitude) were quantified using the function ‘specan’. Power spectra were generated using custom-designed Matlab code in Matlab (Version 9.5.0.1,049,112, R2018b; The MathWorks, Inc., Natick, MA, USA). The following parameters associated with the time domain of the songs were also measured: (1) pulse duration, which refers to duration of the rapid-decay pulse, (2) pulse repetition rate (number of discrete pulses per second), and (3) down-time, which refers to the silent interval between calls and/or pulses. All sound files were uploaded to the Orthoptera Species File Online Version 5.0/5.0 (Cigliano et al. Citation2019).
Specimen identification and curation
Identification was done using descriptions and keys: Hebard (Citation1922); Gorochov (Citation1998, Citation2013, Citation2019); Gorochov and Tan (Citation2011); Tan (Citation2011). The specimens were first preserved in absolute analytic-grade ethanol and later pinned and dry-preserved (after dissection of the tegmen). A single hind leg was preserved in absolute analytic-grade ethanol for future molecular work.
Material examined
Asiophlugis longiuncus Gorochov Citation2013
Material with acoustic data (n = 3): One male (BRU.19.17), Brunei Darussalam, Belait, Andulau Forest Reserve, mixed dipterocarp forest, N4.62,488, E114.51,228, 76.0 ± 6.6 m, 25 February 2019, coll. M.K. Tan and H. Yeo; two males (BRU.19.29, 30), same locality, N4.62,502, E114.51,218, 90.1 ± 8.0 m, 26 February 2019, coll. M.K. Tan and H. Yeo.
Additional material examined (n = 3): One male (BRU.19.36), same locality, N4.63,355, E114.51,096, 75.6 ± 6.6 m, 27 February 2019, coll. M.K. Tan and H. Yeo; two males (BRU.19.49, 50), same locality, N4.28,437, E114.41,873, 39.1 ± 6.5 m, 1 March 2019, coll. M.K. Tan and H. Yeo.
Asiophlugis rete Gorochov Citation1998
Material with acoustic data (n = 1): One male, Singapore, Dairy Farm Nature Park, along Wallace Trail, abandoned plantation, 17 August 2018, coll. M.K. Tan and H. Yeo.
Additional material examined (n = 1): One male (UG.17.14), Peninsular Malaysia, Selangor, Ulu Gombak Field Studies Centre, secondary lowland forest, N3.32,344, E101.74,861, 267.3 ± 6.5 m, 24 February 2017, M.K. Tan, D.M Belabut, Sofwan Badr.
Asiophlugis temasek temasek Gorochov and Tan Citation2011
Material with acoustic data (n = 3): One male, Singapore, Rifle Range Link, secondary lowland forest, 13 July 2018, coll. M.K. Tan; one male, Singapore, Dairy Farm Nature Park, along Wallace Trail, 16 November 2018, coll. M.K. Tan and H. Yeo; one male, Singapore, Upper Seletar Trail, secondary lowland forest, 16 January 2019, coll. M.K. Tan, H. Yeo and X.B. Jin.
Additional material examined (n = 4): One male, Singapore, Rifle Range Link, secondary lowland forest, 13 July 2018, coll. M.K. Tan; one male, Singapore, Dairy Farm Nature Park, near Dairy Farm Quarry, 16 September 2018, coll. M.K. Tan; one male, Singapore, Dairy Farm Nature Park, near Dairy Farm Quarry, 6 November 2018, coll. M.K. Tan and H. Yeo; one male, Singapore, Dairy Farm Nature Park, near Dairy Farm Quarry, 16 September 2018, coll. M.K. Tan
Asiophlugis temasek tioman Gorochov, Citation2011
Material with acoustic data (n = 1): One male, Peninsular Malaysia, Pahang, Pulau Tioman, Paya Village, N2.785, E104.1236, coastal secondary forest, 8 August 2018, coll. M.K. Tan and H. Yeo.
Additional material examined (n = 3): Three males (PT.17.01, 06, 14), same locality, coastal secondary forest, N2.78,524, E104.12,331, 30.3 ± 6.2 m, N2.78,481, E104.12,357, 26.4 ± 6.0 m and N2.78,463, E104.12,400, 28.2 ± 6.1 m, 5 April 2017, coll. M.K. Tan and S.T. Toh
Asiophlugis thaumasia (Hebard 1922)
Material with acoustic data (n = 2): Two males, Singapore, Dairy Farm Nature Park, near Car Park A, abandoned plantation, 9 December 2018 and 31 January 2019, coll. M.K. Tan.
Additional material examined (n = 2): Two males, Singapore, Dairy Farm Nature Park, near Car Park A, abandoned plantation, 21 August 2018, coll. M.K. Tan.
Asiophlugis trusmadi Gorochov, 2011
Material with acoustic data (n = 2): Two males (SDK.19.64, 65), East Malaysia, Sabah, Sandakan, Rainforest Discovery Centre, secondary/primary lowland forest, N5.87,469, E117.94,057, 53.9 ± 6.9 m and N5.87,529, E117.94,171, 50.7 ± 5.7 m, 12 January 2019, coll. M.K. Tan and S.T. Toh.
Measurement of morphological traits
The left tegmen was removed using micro-scissors for imaging and measurements of the tooth distribution. Photomicrograph of the stridulatory apparatus was done using a Canon EOS 50D digital SLR camera with an Olympus OM system extension tube and an Olympus NFK 6.7 × LD (125) lens attached to Olympus BH2 research compound light microscope. Canon Macro Ring Lite MR-14EX were used for lighting and flash. Image-editing and measurements were accomplished using Adobe Photoshop CC 2014 (Adobe Systems Incorporated, San Jose, CA, USA) and/or ImageJ 1.51j8 (Wayne Rasband, Research Services Branch, National Institute of Mental Health, Bethesda, MD, USA) respectively:
The number of teeth on the file was counted.
The inter-tooth distance (mm) was measured from the edge of the cusp of one tooth to the cusp of the next one (Montealegre-Z and Mason Citation2005), from the anal end to the last tooth at the basal end of the left stridulatory file.
The tooth length (mm) as an indication of tooth size was also measured from the anterior to the posterior ends of the cusp of each tooth ().
The tooth density was calculated by dividing the summation of inter-tooth distance by the number of teeth.
Data analysis
To describe the songs for each Asiophlugis, we obtained the arithmetic means and standard deviations of pulse duration, pulse repetition rate, down-time, mean peak frequency and mean frequency. To compare these acoustics parameters, we fitted the linear mixed effects models (LMMs) for pulse duration, pulse repetition rate, down-time and mean peak and peak frequencies using the ‘lmer’ function from the R package ‘lme4ʹ (Bates et al. Citation2014). The frequency and time parameters were natural log-transformed. In each model, we used species as a fixed effect and individual katydid as a random intercept since more than one recording were obtained for each individual. Temperature was also used as a random slope since it is known to affect the acoustic parameters. Estimated marginal means (or least-squares means) were computed as the sampling was imbalanced. These were done using the ‘emmeans’ function from the R package ‘emmeans’ (Lenth Citation2018). To compare the stridulum anatomy of Asiophlugis, we compared the inter-tooth distance and tooth length between the different Asiophlugis by fitting trend lines using generalised additive models to visualise the tooth distributions on the stridulatory file based on inter-tooth distances and tooth length variation for each species. We also compared the tooth density between taxa using Kruskal–Wallis rank sum test.
Results
Part I: calling songs
We analysed the songs of Asiophlugis from a total of 163 sound files: A. longiuncus (n = 42), A. rete (n = 2), A. temasek temasek (n = 24), A. temasek tioman (n = 24), A. thaumasia (n = 42) and A. trusmadi (n = 29). The songs for the six Asiophlugis are similar in that they consist of rapid-decay pulses with broad-based spectrum (). Such rapid-decay pulse consists of a train of complex waves showing rapid decay after only one or a few cycles at maximum amplitude. The songs are also completely inaudible to the human ear (). We excluded A. rete for the statistical modelling of the calling songs since there were only two recordings, much fewer than other species.
Figure 3. Oscillograms of the Asiophlugis songs: A. longiuncus (a), A. rete (b), A. temasek temasek (c), A. temasek tioman (d), A. thaumasia (e), and A. trusmadi (f). Scale bar: 0.1 s

A. longiuncus song is made up of a rapid-decay pulse, within which impulses typically paired in time ()). Each pulse consisting of two rapid-decay impulses is about 16.7 ± 2.8 ms in duration ()). The second impulse is of much lower amplitude than the first one (), )). The average pulse repetition rate is 5.1 ± 1.5 s−1. The average down-time between consecutive pulses is 175.5 ± 84.9 ms. The mean peak frequency is 55.5 ± 1.4 kHz and mean frequency is 58.1 ± 1.3 kHz ()).
Figure 4. Pulse structure of the Asiophlugis: A. longiuncus (a), A. rete (b), A. temasek temasek (c), A. temasek tioman (d), A. thaumasia (e), and A. trusmadi (f). Scale bar: 0.2 ms

Figure 5. Power spectra of the Asiophlugis songs: A. longiuncus (a), A. rete (b), A. temasek temasek (c), A. temasek tioman (d), A. thaumasia (e), and A. trusmadi (f)

A. rete song consists of a rapid-decay pulse which repeats over time ()). Each pulse is 2.7–4.4 ms in duration ()). The average pulse repetition rate is 4–9 s−1. The average down-time between consecutive pulses is 22.8 ms. The mean peak frequency is 33.3–35.3 kHz and mean frequency is 47.9–52.0 kHz. ())
A. temasek temasek song is similar to A. rete song. It is made up of a rapid-decay pulse and repeats in a series ()). Each pulse is about 7.7 ± 3.5 ms in duration ()). The song has two modes, even within the same individual ()). The first mode has a high average pulse repetition rate of 52.2 ± 28.1 s−1, and low average down-time between pulses of 13.5 ± 8.0 ms. The second mode has a low average pulse repetition rate of 12.3 ± 4.0 s−1 and high average down-time between pulses of 67.6 ± 22.7 ms. The mean peak frequency is 57.8 ± 2.9 kHz and mean frequency is 57.2 ± 1.2 kHz ()).
Figure 6. Histograms of down time and pulse repetition rate of A. temasek temasek (a, b) and A. temasek tioman (c, d), Arrow indicates a gap between the two modes of song

A. temasek tioman song is also made up of a rapid-decay pulse and repeats in a series ()). We only observed one mode for this subspecies ()), unlike subspecies temasek. Each pulse is about 5.0 ± 0.9 ms in duration ()), shorter than that of subspecies temasek. The average pulse repetition rate is 43.1 ± 21.6 s−1. The average down-time between consecutive pulses is 9.3 ± 5.1 ms. The mean peak frequency is 46.2 ± 6.2 kHz and mean frequency is 50.9 ± 2.1 kHz, lower than that of subspecies temasek ()).
A. thaumasia song is quite distinct from that of the above taxa in that it consists of a series of four rapid-decay pulses of very similar amplitude ()). The duration of each series is 166.9 ± 49.5 ms and consecutive series are separated by a long and highly variable down-time of average 744.8 ± 378.4 ms. Each pulse is about 30.0 ± 4.3 ms in duration ()), distinctively higher than the other taxa. The average pulse repetition rate is 4.3 ± 1.1 s−1. The average down-time between consecutive rapid-delay pulses within a series is about 10.9 ± 3.2 ms. The mean peak frequency is 44.8 ± 2.1 kHz and mean frequency is 48.3 ± 0.9 kHz ()).
A. trusmadi song, similar to A. longiuncus, is made up of a rapid-decay pulse, within which impulses typically occur in triplet in time ()). Each pulse consisting of three rapid-decay impulses is about 18.5 ± 3.0 ms in duration ()). The average pulse repetition rate is 1.5 ± 0.5 s−1. The average down-time between consecutive pulses is 687.5 ± 303.2 ms. The mean peak frequency is 59.5 ± 5.2 kHz and mean frequency is 61.7 ± 0.7 kHz ()).
The time domain of the songs was not different between the two subspecies of A. temasek and between A. longiuncus, A. thaumasia and A. trusmadi. The pulse duration of A. longiuncus, A. thaumasia and A. trusmadi were longer than that of A. temasek (), ). Pulse repetition rate was highly variable for A. temasek, which is drastically different from the other species, and greater than the other species (), ). This can be attributed to the presence of two song modes for A. temasek temasek. Pulse repetition rate also did not appear to be different between A. longiuncus, A. thaumasia and A. trusmadi (), ). We observed similar patterns for down-time since longer down-time is associated with higher pulse repetition rate (), ). Down-time for A. longiuncus and A. thaumasia was also highly variable compared to the other taxa (), ), which can be explained by the irregular calling behaviour in these two species.
Table 1. Pairwise comparison of the song signatures between the five taxa of Asiophlugis. The values indicate the ratio of the parameter between two taxa. A. rete was excluded from the LMMs since there were only two recordings, much fewer than other species Significance is denoted as follows: * P < 0.05; ** P < 0.01; *** P < 0.001; ns not significant
Figure 7. Mean plots on the comparison of the time and frequency domains of the Asiophlugis song. The dot represents the least-square means and the horizontal line represents the 95% confidence interval (CI). The mean of each parameter was included for A. rete but the 95% CIs were not included since it was excluded from the statistical analysis

The frequency domain showed different patterns with the time domain. The frequency domain between the two subspecies of A. temasek were drastically different (), ). Yet, there was no evidence of such difference between A. temasek temasek and A. longiuncus and A. trusmadi; and there was no evidence of such difference observed between A. temasek tioman and A. thaumasia (), ). Both peak and mean frequencies exhibited similar patterns between the Asiophlugis (), ).
Syntopic species from Singapore had different song signatures but there was weak evidence that song signatures of sympatric species are more different than allopatric species. The time and frequency domains of the Bornean species (i.e., A. logniuncus and A. trusmadi) were very similar () although they should exhibit greater differences than with allopatric species from the Malay Peninsula. Moreover, frequency domain of A. temasek temasek was just as different as that of sympatric and allopatric species from Malay Peninsula (,)) and the time domain of A. temasek temasek was more different in sympatric species than allopatric A. temasek tioman ()).
Part II: stridulum anatomy
Stridulatory files on the left tegmen are visibly different between the six Asiophlugis species (). The stridulatory files are generally straight and slender, although that of A. longiuncus ()) and A. trusmadi ()) are slightly sigmoidal.
Figure 8. Stridulatory file on the left tegmen of Asiophlugis: A. longiuncus (a), A. rete (b), A. temasek temasek (c), A. temasek tioman (d), A. thaumasia (e), and A. trusmadi (f). TL refers to tooth length. Scale bars: 0.1 mm

Tooth distributions based on inter-tooth distances () and tooth lengths () between subspecies of A. temasek are similar. The inter-tooth distances of A. temasek temasek and A. temasek tioman are fairly even throughout the file, in contrast to other species in which the inter-tooth distance tends to increase to a maximum before decreasing drastically (). Likewise, A. temasek temasek and A. temasek tioman have evenly-sized teeth, except at both ends of the files (). This was also observed in A. rete, whereas A. longiuncus, A. thaumasia and A. trusmadi have a few very large teeth in the middle and much smaller teeth at both ends of the files (). Tooth density was significantly different between all species (χ = 13.3, p-value = 0.021), but the differences are slightly smaller between the two subspecies of A. temasek ().
Figure 9. Tooth distribution on the stridulatory file of Asiophlugis based on inter-tooth distances. Trend lines are generalised additive model fits. A. longiuncus, r2 = 0.91; A. rete, r2 = 0.89; A. temasek temasek, r2 = 0.28; A. temasek tioman, r2 = 0.69; A. thaumasia, r2 = 0.87; and A. trusmadi, r2 = 0.63

Figure 10. Tooth distribution on the stridulatory file of Asiophlugis based on tooth length variation. Trend lines are generalised additive model fits. A. longiuncus, r2 = 0.88; A. rete, r2 = 0.94; A. temasek temasek, r2 = 0.57; A. temasek tioman, r2 = 0.70; A. thaumasia, r2 = 0.95; and A. trusmadi, r2 = 0.95

Discussion
Songs of Asiophlugis
We recorded and analysed the calling songs of six taxa of Asiophlugis (among 16 known species), where previously the song of only one species from the Asiatic genera of Phlugidini was known (see Helfert and Sänger Citation1998). The calling frequencies used by all species recorded here are ultrasonic (>40 kHz) and have a broad distribution of carrier frequencies (as shown by the spectrograms). We also examined whether ultrasonic songs and stridulum anatomy of Asiophlugis and Phlugidini are useful evidence for species delimitation, which has been challenging so far.
The songs of Phlugidini are similar in that they are entirely ultrasonic, supporting the fact that Meconematinae species usually employ ultrasonic calls (Helfert and Sänger Citation1998; Montealegre-Z et al. Citation2006, Citation2017; Chamorro-Rengifo et al. Citation2014; Chamorro-Rengifo and Braun Citation2016; Sarria-S et al. Citation2017). A. rete (A. thai junior synonym) sings at frequencies between 30 and 50 kHz (Helfert and Sänger Citation1998) and neotropical species such as Phlugis ocraceovittata Piza, 1960 also sings at frequencies between 40 and 60 kHz (Chamorro-Rengifo and Braun Citation2016). The call of species of the tribes Phlugidini and Phisidini (e.g., Phugis spp, Arachnoscelis, Arachnoides, Supersonus spp.) incorporates a sequence of short impulses or very discrete pulses. These impulses likely result from scraper pauses along the file teeth, its deformation, storage of elastic energy, and slipping free and forward drive over a few file teeth at a higher rate that likely exceeds the normal scraper velocity being driven by muscles in a sustained closing wing phase. The elastic energy that was stored then released as kinetic energy leads to increase in scraper velocity, thus enhancing extremely fast tooth strikes, which end in ultrasonic frequency generated beyond what the wing muscle alone can generate (Montealegre-Z et al. Citation2006). Different from other singing Ensifera species, like crickets (Koch et al. Citation1988) or other katydids (Montealegre-Z and Mason Citation2005; Montealegre-Z and Postles Citation2010) where a systematic distribution of teeth is necessary for maintaining a constant tooth strike rate and a coherent tonal pulse, such systematic organisation is not required in the extreme ultrasonic Phlugidini. Instead, a few large and separated teeth seem to be appropriate for scraper distortion and recoiling at high speeds.
Can stridulum anatomy assist with inferring species identity?
Since the tooth distributions are similar between the two subspecies of A. temasek, but different between species of Asiophlugis (especially for inter-tooth distance), the morphology of stridulatory files seems to be useful characters in differentiating species and subspecies. Prior to this study, the stridulatory files of only A. rete (as A. thai junior synonym) and A. temasek temasek were described and illustrated using scanning electron microscopy and hand-drawing, respectively (Helfert and Sänger Citation1998; Gorochov and Tan Citation2011). However, they were never used for species identification. Here, we support previous work on other katydids (e.g., Naskrecki and Rentz Citation2010; Sarria-S et al. Citation2014; Heller and Hemp Citation2018) that propose that the morphology of the stridulatory file is congruent with general morphology in delimiting species in Phlugidini.
Can the call structure assist with inferring species identity?
Overall, both the descriptive comparison and statistical analysis of the song parameters indicate that song signatures of the six Asiophlugis are all different, even though the time and frequency domains do not necessarily differ between the taxa in a similar manner. Furthermore, songs of A. longiunus and A. trusmadi differ by the paired and tripled rapid-decay impulses, respectively, although they are not different in the other acoustic parameters, including the frequency domain and pulse duration and repetition rate. The songs between the two subspecies of A. temasek can also be considered as different as congeneric differences, as corroborated by the differences in the frequency domain and the presence of two modes in only the subspecies temasek. These underscore that one acoustic parameter cannot give an accurate and reliable indication for species identity and it is crucial to examine as many parameters as possible to avoid biases when interpreting the songs.
The difference in the frequency domain of the A. temasek subspecies songs signifies that song signatures do not entirely reflect the greater similarities between subspecies and greater differences between species we initially predicted and have observed for the morphology and stridulum anatomy. The incongruency between song signatures and stridulum anatomy is not unexpected since the differences in stridulum anatomy do not necessarily produce songs of different frequencies. This has been observed in the Aerotegmina katydids (e.g., Heller and Hemp Citation2018). Carrier frequency is mainly dictated by the resonator anatomy (mirror and harp but see; Montealegre-Z and Postles Citation2010; Montealegre-Z Citation2012; Montealegre-Z et al. Citation2017).
Subspecies status of A. temasek based on songs and stridulum anatomy
This led us to question whether the two subspecies of A. temasek should more rightly be considered as two different species. Although differences in their morphology (see Gorochov and Tan Citation2011) and stridulum anatomy are probably not as clear-cut as species differences in Asiophlugis, one can also argue that the differences in the calling songs may be considered species specific. Our hypothesis is that A. temasek temasek and A. temasek tioman probably represent two recently diverged sister species. But there is currently a lack of phylogenetic analysis that may inform about their relationships. We also only know that the subspecies temasek occurs in Singapore and subspecies tioman occurs in Pulau Tioman. The differences in their songs may be an artefact of sympatry but there is currently no way to confirm this and infer species boundaries until sampling in Malay Peninsula between these two islands and congeneric data from other parts of the genus’ distribution (e.g., east of Wallace Line) have been done. Thus, to avoid confusion, we tentatively keep the status of the two taxa as status quo until we have more evidence to suggest otherwise.
Synotypic species of Asiophlugis from Singapore
The song signatures of syntopic species from Singapore (they can be found in the same forest patch of Dairy Farm Nature Park) are also drastically different, even though there was no clear evidence that songs are more different in sympatry than allopatry. Such difference in the songs of A. temasek temasek, A. thaumasia and A. rete can be a crucial pre-mating isolating mechanism to avoid hybridisation between these syntopic species. We have also observed the occurrence of syntopic species within Brunei Darussalam, Pulau Tioman and Sandakan but could not find the males and/or record their songs. This implies that further sampling is needed to validate whether the song differences between syntopic species observed in Singapore can be generalised.
Conclusions
The present study provides the most in-depth investigation into the bioacoustics of Asiatic Phlugidini to date. This is despite our study spanning only across Asiophlugis species from Malay Peninsula and Borneo. There are other representatives of the Asiatic Phlugidini, including Papuaphlugis Gorochov Citation2012 and Stenophlugis, Citation2012, but of which we do not know anything about their acoustics. These signify that the study of acoustic communication in Asiatic Phlugidini is still in its early stage. We demonstrate that the stridulum morphology and acoustics can provide additional meaningful evidence in species delimitation. As such, we recommend that species description and species delimitation exercises in the future should include acoustic information, since it is expected that more species (see Gorochov Citation2019) and perhaps genera of Phlugidini to be discovered in the understudied regions of Southeast Asia (Tan et al. Citation2017). Our study aims to provide a baseline for more research on ultrasound communication in Asiophlugis, Phlugidini and other understudied ultrasonic-singing katydids.
Author contributions
MKT and FMZ contributed equally to this study.
Acknowledgements
The authors are thankful to Huiqing Yeo and Xingbao Jin (in Singapore and Brunei Darussalam), Siew Tin Toh (in Pulau Tioman and Sandakan), Amira Aqilah Muhammad and Sofwan Badr (in Ulu Gombak), Momin Binti, John Lee Yukang and Saudi Bintang (in Sandakan) for field assistance. The permissions for collecting material were granted by the Forestry Department, Ministry of Primary Resources and Tourism, Brunei Darussalam (JPH/PDK/01 Pt 2); the Sabah Biodiversity Centre (JKM/MBS.1000–2/3 JLD.3 (99)) (for Sandakan); the National Parks Board (NP/RP18–064), Singapore; and the Research Promotion and Co-Ordination Committee, Economic Planning Unit, Prime Minister’s Department (UPE: 40/200/19/3395), Malaysia and the University Malaya (for Pulau Tioman and Ulu Gombak).
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- Araya‐Salas M, Smith‐Vidaurre G. 2017. warbleR: an R package to streamline analysis of animal acoustic signals. Methods Ecol Evol. 8(2):184–191.
- Bates D, Maechler M, Bolker B, Walker S, Christensen RHB, Singmann H, Dai B 2014. lme4: linear mixed-effects models using Eigen and S4 (Version 1.1–7).
- Chamorro-Rengifo J, Braun H. 2016. Phlugis ocraceovittata and its ultrasonic calling song (Orthoptera, Tettigoniidae, Phlugidini). Zootaxa. 4107(3):439–443.
- Chamorro-Rengifo J, Braun H, Lopes-Andrade C. 2014. The secret stridulatory file under the right tegmen in katydids (Orthoptera, Ensifera, Tettigonioidea). Zootaxa. 3821(5):590–596.
- Chamorro-Rengifo J, Olivier RDS. 2017. A new genus of Phlugidini (Orthoptera: Tettigoniidae: Meconematinae) with asymmetrical mandibles. Zootaxa. 4286(3):391–400.
- Cigliano MM, Braun H, Eades DC, Otte D 2019. Orthoptera species file online. version 5 (5.0), (accessed 2018 Mar 14). http://orthoptera.speciesfile.org/HomePage/Orthoptera/HomePage.aspx(open in a new window)
- Core Team R. 2018. R: A language and environment for statistical computing. Vienna (Austria): R Foundation for Statistical Computing.
- Fitzpatrick MJ, Gray DA. 2001. Divergence between the courtship songs of the field crickets Gryllus texensis and Gryllus rubens (Orthoptera, Gryllidae). Ethology. 107(12):1075–1085.
- Gorochov AV. 1998. New and little known Meconematinae of the tribes Meconematini and Phlugidini (Orthoptera: Tettigoniidae). Zoosyst Ross. 7(1):101–131.
- Gorochov AV. 2012. Systematics of the American katydids (Orthoptera: Tettigoniidae). Communication 2. Zoosyst Ross. 316(4):285–306.
- Gorochov AV. 2013. Taxonomy of the katydids (Orthoptera: Tettigoniidae) from East Asia and adjacent islands. Communication 6. Far East Entomol. 259:1–12.
- Gorochov AV. 2019. Taxonomy of the katydids (Orthoptera: Tettigoniidae) from East Asia and adjacent islands. Communication 12. Far East Entomol. 379:1–24.
- Gorochov AV, Tan MK. 2011. New katydids of the genus Asiophlugis Gorochov (Orthoptera: Tettigoniidae: Meconematinae) from Singapore and Malaysia. Russ Entomol J. 20(2):129–133.
- Hebard M. 1922. Studies in Malayan, Melanesian and Australian Tettigoniidae (Orthoptera). Proceedings of the Academy of Natural Sciences of Philadelphia. 74:121–299.
- Helfert B, Sänger K. 1998. Phlugis thai sp. n., a new Phlugidini (Insecta: Ensifera, Tettigoniidae, Meconematinae) from Thailand. Annalen des Naturhistorischen Museums in Wien. Serie B Für Botanik Und Zoologie. 100:61–68.
- Heller KG, Hemp C. 2018. Extremely divergent song types in the genus Aerotegmina Hemp (Orthoptera: Tettigoniidae: Hexacentrinae) and the description of a new species from the Eastern Arc Mountains of Tanzania (East Africa). Bioacoustics. 28(3):269–285.
- Izzo AS, Gray DA. 2004. Cricket song in sympatry: species specificity of song without reproductive character displacement in Gryllus rubens. Ann Entomol Soc Am. 97(4):831–837.
- Kevan DKM, Jin XB. 1993. Remarks in the tribe Phlugidini Eichler and recognition of the new taxa from the Indo-Malayan region and east Africa (Grylloptera: Tettigonioidea: Meconematidae). Invertebr Syst. 7(6):1589–1610.
- Koch UT, Elliott CJH, Schäffner KH, Kleindienst HU. 1988. The mechanics of stridulation of the cricket Gryllus campestris. J Comp Physiol A. 162:213–223.
- Lenth R. 2018. Emmeans: estimated marginal means, aka least-squares means. R Package Version. 1:2.
- Montealegre-Z F. 2009. Scale effects and constraints for sound production in katydids (Orthoptera: Tettigoniidae): generator size constrains signal parameters. J Evol Biol. 22(2):355–366.
- Montealegre-Z F. 2012. Reverse stridulatory wing motion produces highly resonant calls in a neotropical katydid (Orthoptera: Tettigoniidae: Pseudophyllinae). J Insect Physiol. 58:116–124.
- Montealegre-Z F, Mason AC. 2005. The mechanics of sound production in Panacanthus pallicornis (Orthoptera: Tettigoniidae: Conocephalinae): the stridulatory motor patterns. J Exp Biol. 208(7):1219–1237.
- Montealegre-Z F, Morris GK, Mason AC. 2006. Generation of extreme ultrasonics in rainforest katydids. J Exp Biol. 209:4923–4937.
- Montealegre-Z F, Ogden J, Jonsson T, Soulsbury CD. 2017. Morphological determinants of carrier frequency signal in katydids (Orthoptera): a comparative analysis using biophysical evidence of wing vibration. J Evol Biol. 30(11):2068–2078.
- Montealegre-Z F, Postles M. 2010. Resonant stridulation in Copiphora gorgonensis (Tettigoniidae: Copiphorini) an endemic species from the natural National Park Gorgona, Colombia. J Orthoptera Res. 19(2):345–353.
- Mugleston JD, Naegle M, Song H, Whiting MF. 2018. A comprehensive phylogeny of Tettigoniidae (Orthoptera: Ensifera) reveals extensive ecomorph convergence and widespread taxonomic incongruence. Insect Syst Divers. 2(4):5.
- Naskrecki P, Rentz DC. 2010. Studies in the orthopteran fauna of Melanesia: new katydids of the tribe Agraeciini from Papua New Guinea (Orthoptera: Tettigoniidae: Conocephalinae). Zootaxa. 2664:1–35.
- Nickle DA. 2003. New neotropical species of the genus Phlugis (Orthoptera: Tettigoniidae: Meconematinae). J Orthoptera Res. 21(1):37–56.
- Nickle DA. 2005. Additional notes on the genus Phlugis (Orthoptera: Tettigoniidae: Meconematinae) with the descriptions of two new arboreal species from Costa Rica. J Orthoptera Res. 14(1):57–62.
- Sarria-S FA, Chivers BD, Soulsbury CD, Montealegre-Z F. 2017. Non-invasive biophysical measurement of travelling waves in the insect inner ear. R Soc Open Sci. 4(5):170171.
- Sarria-S FA, Morris GK, Jackson J, Windmill JFC, Montealegre-Z F. 2014. Shrinking wings for ultrasonic pitch production: hyperintense ultra-short-wavelength calls in a new genus of neotropical katydids (Orthoptera: Tettigoniidae). PLoS One. 6:e98708.
- Tan MK. 2011. The species of Asiophlugis Gorochov, 1998 in Singapore (Orthoptera: Tettigoniidae: Meconematinae). Nat Singapore. 4:233–239.
- Tan MK, Choi J, Shankar N. 2017. Trends in new species discovery of Orthoptera (Insecta) from Southeast Asia. Zootaxa. 4238(1):127–134.