Introduction

According to sexual selection theory (Andersson, 1994), many behavioural and morphological traits specifically associated with mating evolved as a result of intra- and intersexual selection. In the latter process one sex, generally the female, mates preferentially with certain members of the other sex. Such preferences may be based on material benefits that a female directly obtains by mating with certain males: resources defended on a territory or transferred at mating, parental care, avoidance of parasite infection and ease of localization. Nonetheless, in many animal species females discriminate among potential mates even when they do not receive any apparent material benefit directly influencing their survival or fecundity. In such cases, the origin and maintenance of female mate preference remain a major problem in evolutionary biology.

Where no apparent direct, material benefits are conferred at mating, it is proposed that discriminating females receive indirect, ‘genetic’ benefits from the chosen male’s gametes. In the ‘arbitrary’ (or ‘Fisherian’) model of indirect sexual selection, the genetic benefits are ‘attractiveness genes’ that are expressed in the next generation. By choosing attractive males, a female is likely to produce sons that display their father’s aesthetically appealing features and later succeed in the ‘sexual selection arena’ (Lande, 1981). Alternatively, in ‘good-genes’ (or ‘viability indicator’) models the genetic benefits also include ‘viability genes’ that enhance offspring survival and growth (Pomiankowski, 1988). Here, females who mate with attractive males produce offspring that not only display their father’s attractive features (if they are sons) but have superior survivorship and maturation characteristics as well.

Both the arbitrary and good-genes models are coevolutionary in that the male character and the female preference coevolve in a reciprocally reinforcing manner (Houde, 1997). Coevolutionary models of indirect sexual selection assume three fundamental genetic conditions: significant additive genetic variance (heritability) in the male character, heritable female preference traits, and genetic covariance between male and female traits generated by female preference for particular male phenotypes. Such covariance reinforces the sexual selection process and may lead to exaggeration of both the male character and the female preference. Various recent studies employing quantitative genetic approaches have confirmed heritability for male characters (e.g. see references in Andersson, 1994). Nonetheless, the factors that maintain additive genetic variance in natural populations remain largely unknown. This problem, often termed ‘the paradox of the lek’ (Ritchie, 1996), is receiving renewed theoretical treatment by various workers (e.g. Pomiankowski & Møller, 1995; Rowe & Houle, 1996).

Here, we report a quantitative genetic study of male signal characters in the lesser waxmoth Achroia grisella (Lepidoptera: Pyralidae). We conducted two parallel breeding experiments and calculated the heritability of and correlations among those signal characters subject to directional sexual selection by female choice. Because previous work has suggested that a good-genes mechanism of sexual selection may operate in A. grisella under certain environmental conditions (Jia & Greenfield, 1997), we also determined the correlations between signal characters and development time. Heritabilities were measured under two different environmental conditions, and variation in these values between conditions suggests ways in which the observed additive genetic variance is maintained.

Natural history and acoustic signalling in Achroia grisella

Achroia grisella is a cosmopolitan symbiont of honeybee (Apis mellifera) colonies where the larvae feed mostly on residual stored food and organic detritus. Adults, having atrophied proboscises, neither feed nor drink but remain in the vicinity of the honeybee colony during their brief (7–14 days) life spans. Both males and females exit the hive immediately after eclosion, which normally occurs close to dusk, expand their wings, and commence sexual activities shortly thereafter on or near the hive (Greenfield & Coffelt, 1983).

Unlike most moths, in which pair formation is typified by female olfactory signalling and male searching, A. grisella males produce ultrasonic advertisement signals that attract receptive females within a radius of several metres (Spangler et al., 1984). Both males and females hear sounds extending over a frequency range from 30 to 120 kHz with a pair of abdominal tympanal organs. Females usually run on the substrate, rather than fly, during phonotaxis towards a signalling male. The males remain perched and stationary while signalling and do not move towards females until the pair contact. An individual male may continue acoustic signalling without interruption for 6–10 h each night. Males also emit an odour (a combination of two aldehydes) from forewing glands, but its function is unclear. Experiments using simulated signals show that the ultrasonic signal alone is as attractive as a live male, i.e. the acoustic signal is both necessary and sufficient for pair formation (Spangler et al., 1984).

Achroia grisella males produce ultrasonic signals by wing fanning, which causes a pair of tymbals on the tegulae (axillary sclerites on the thorax) to buckle in and out. Females, who are silent, lack such tymbals. Males remain perched and stationary while wing fanning and produce ≈50 wingstrokes s–1 at 25°C. A male’s tymbals buckle twice during each wingstroke (complete cycle of wing movement), once on the upstroke and once on the downstroke. Each tymbal buckling results in emission of a highly damped 100-μs pulse of sound. Owing to slightly asynchronous wing movement, a male’s two tymbals produce a pair of pulses during each upstroke and downstroke that are usually separated by a brief (100–400 μs) gap. The sound pulses are loud (amplitude ≈95 dB sound pressure level at 1 cm; 0 dB re 20 μPa), broadcast rather omnidirectionally, and include frequencies from 70 to 130 kHz (peak frequency ≈100 kHz) (Spangler et al., 1984; Snedden et al., 1994). Females prefer male signals that are louder (measured as peak pulse amplitude), delivered at a faster rate (greater wingstroke frequency), and include longer asynchrony intervals (the temporal gap between pulses produced by right and left tymbals) (Jang & Greenfield, 1996, 1998; Jang, 1997). These three signal characters vary greatly between males but are highly repeatable within individuals (r≥ 0.66; Jang et al., 1997).

Respirometric studies have shown that a relatively high expenditure of metabolic energy is necessary, but not sufficient, to produce attractive signals in A. grisella (Reinhold et al., 1998). This requirement reflects the positive relationship between pulse rate and energy expenditure. Under certain favourable diet and temperature conditions the male offspring of attractive males develop faster than the offspring of unattractive males (Jia & Greenfield, 1997). However, this relationship is context dependent, being absent under other environmental conditions.

Females normally mate only once and are then unreceptive (Greenfield & Coffelt, 1983). Males can mate approximately once every 24 h, and female fecundity is unaffected by male mating history. Thus, it is unlikely that material benefits in the form of critical spermatophore materials, which represent less than 1% of the male’s body mass, are transferred to the female at mating. Because paternal care and other material benefits appear to be inconsequential, female choice for male signals in A. grisella may be based primarily on genetic benefits and maintained by an indirect sexual selection mechanism.

Materials and methods

Population studied and general conditions

A laboratory population of A. grisella was founded from several hundred adults collected at honeybee colonies near Auburn, AL. We employed a breeding protocol that maximized effective population size in this stock population. Larvae were reared on a standardized diet of cereals, Brewer’s yeast, glycerol, honey, beeswax and water. Temperature (≈25.5°C), photoperiod (12:12 L:D), and density (≈1 larva per 1 g diet) were maintained such that generations were approximately 45 days long. The experiments reported below were conducted on insects from the 16th to the 18th generation. To insure that all adults were unmated prior to testing, we placed late-instar larvae individually in 30-mL clear plastic containers with ≈1 g of diet to complete their development.

We recorded male signals in one of two similarly constructed semianechoic rooms (3 × 3.5 × 2.5 m) illuminated by 25-W red incandescent bulbs ( <1 lux). Tests were conducted at 25–26°C during the first 4 h of the scotophase, the period of peak signalling and sexual activity (Greenfield & Coffelt, 1983; Jang et al., 1997).

Experimental design

Quantitative genetic parameters for signalling and development were measured for insects taken from the stock population described above. Heritabilities and genetic correlations were estimated using half-sib analysis (Falconer, 1981). Both of the experiments described below were conducted during the same 4-month period (October 1996–January 1997).

Experiment 1: Collective rearing A North Carolina Experiment I design (NCI) (Kearsey & Pooni, 1996) was used to estimate genetic variance in signal characters and development time. Virgin males and females were collected randomly from the stock population and mated when they were 1–3 days old. Each male was paired with a different female on three consecutive days. Approximately 40 eggs from each pairing were chosen randomly and placed on 25–30 g of diet in a clear plastic box (16 × 12 × 6.5 cm deep). If one or more of the three females mated to a given male failed to produce eggs, none of the offspring of that male was tested.

The offspring of each male–female pair were reared collectively until they were pupae or last-instar larvae (approximately 35 days after oviposition). Pupae (or larvae) were then placed individually in 30-mL clear plastic cups with ≈1 g diet. These cups were checked daily until adult eclosion. Development times (oviposition to adult eclosion) for all male and female offspring were recorded for 47 sires paired with 3 dams each. Mean (SD) offspring per sire and per dam were 32.8 (7.0) and 10.9 (3.4), respectively.

The ultrasonic signals of all male offspring were recorded on the 1st or 2nd day post emergence when they were 24–48 or 48–72 h old. Although male-emitted signal characters change with age, these changes occur primarily when males are 3 or more days old; signal characters of 1- and 2-day-old males do not differ significantly (Snedden et al., 1994;).

Experiment 2: Individual rearing Males were paired with females as above, but offspring were transferred as 2nd instar larvae to individual 30-mL plastic cups with 1 g of diet. Thus, most larval development occurred in the absence of competition. We noted the development times of all eclosing offspring and recorded the ultrasonic signals of every third eclosing male. Thirty-three males were paired with an average of 1.9 females, and the mean (SD) number of offspring per mating was 7.1 (1.2). The basic NCI design was also used in this experiment, but the offspring of males that mated with fewer than three females were reared and included in the analysis.

In both experiments larvae were reared at the same mean temperature (≈25.5°C). However, the daily temperature range in expt 1 (24–27°C) exceeded that in expt 2 (25–26°C). All signal recordings in a given experiment were made in the same semianechoic room. Temperatures were ≈25.5°C (range 25–26°C) at the time of recording in both experiments.

Signal recording and analysis

We recorded ultrasonic signals of males held individually in cylindrical screen cages (2.5 × 2.5 cm diameter). The cage screen had no effect on transmission of the ultrasonic signals (Jang et al., 1997). Cages were separated by at least 25 cm to minimize any potential influence of neighbours on a male’s signalling.

We recorded each male with a 0.64-cm condenser microphone (ACO Pacifica, model 7016; frequency response ±6 dB from 10 to 160 Hz) fitted with a preamplifier (ACO model 4012). The microphone was positioned outside the male’s cage and orientated towards its dorsum at a standard distance. Output from the microphone–preamplifier unit was amplified with an ACO Pacifica model PS9200 amplifier and sent to the input port of a sound card installed in a personal computer (80486 processor; 66 MHz clock speed). The sound card (SuperSound Engineering Version, SiliconSoft, San Jose, CA) included 8-bit analogue–digital and digital–analogue converters, each with a maximum sampling rate of 1 MHz. In operation, signals were sampled at 297 or 398 kHz and saved to a computer file using custom software developed for ultrasound recording.

We recorded the signal of each male for 0.164 s (7–10 complete cycles of wing movement). Three consecutive cycles were analysed. Peak amplitude was calculated as the average maximum amplitude (measured on a linear scale) of all pulses in the three cycles. Signal rate was computed as the reciprocal of mean cycle duration. Asynchrony interval was measured as the average delay between pulses produced by the right and left wings during the up and down wing strokes.

Genetic parameter estimation

Genetic parameters were estimated separately for each experiment. Heritability estimates were computed using nested analysis of variance (Model II) for unequal sample sizes to allocate observed variation among sire and dam components (Sokal & Rohlf, 1981, pp. 294–299). Mean squares for the sire and dam components were pooled when the dam component was not significant (Sokal & Rohlf, 1981; pp. 285–286). Heritability estimates were then made using the method described by Becker (1975), pp. 50–52. Variances of the variance components were estimated with the procedures of Searle 1971, (pp. 447–448) and used to estimate standard errors for the heritabilities (Becker, 1975; p. 45).

We also calculated the coefficients of additive genetic variance for signal characters and development time. These values provide additional information about the ability of a population to respond to selection (Houle, 1992).

Phenotypic, genetic and environmental correlations were estimated for each experiment using methods described by Becker 1975; (pp. 105–110). Standard errors for genetic correlations were estimated using a jackknife method. Separate estimates were made by omitting the offspring of one male and calculating genetic correlations based on the offspring of the remaining males. Estimates were made for all possible subsets thus constructed. Variation in this set was then used to calculate variance in the genetic correlation for the entire data set. Standard errors for phenotypic and environmental correlations were estimated using the methods outlined by Becker 1975; (pp. 105–110).

Genetic correlations are reported for all pairs of characters regardless of their heritabilities. However, genetic correlations between traits that are not significantly heritable are difficult to interpret and could reflect Type I errors (rejection of a true null hypothesis) or nonadditive genetic variance.

Distributions for each variable within each experiment were compared to the normal distribution using the Kolmogorov–Smirnov test (Sokal & Rohlf, 1981; p. 440). Where necessary and possible, data were transformed to approximate normality. Transformations used for particular variables are listed in

Analysis of variance was used to assess the significance of differences in mean signal characters between the two experiments. Homogeneity of variance was tested using the Levene test (Norusis, 1990). Transformed data were used when appropriate for all comparisons.

Results

Signal characters and development

Although the insects used in the two experiments were derived from the same stock population at approximately the same time, we observed significant differences in signal character means and variances (Tables 1). Males reared collectively produced signals with lower peak amplitudes and longer asynchrony intervals than those reared individually. Collectively reared males also developed more slowly. Because available food was more than adequate for the number of larvae reared, the observed difference in development may have resulted from extended larval interactions under collective rearing. This result is consistent with that reported by Jia & Greenfield (1997).

Table 1 Phenotypic means† (±SE) of ultrasonic signal characters and development time for male offspring of Achroia grisella reared collectively (expt 1) and individually (expt 2) using NCI half-sib designs

Signal characters and development times were more variable among males reared individually. Variances for peak amplitude, signal rate and development time in the two experiments were not homogeneous (Levene test; P< 0.001). Possibly, the reduced variance among males reared collectively resulted from the common environments established by this rearing condition.

Heritability of signal characters

All of the measured signal characters were significantly heritable for males reared collectively (P< 0.05; (Table 2). Heritabilities for peak amplitude and signal rate both exceeded 0.40. Although heritability for asynchrony interval was only 0.08, it was statistically significant.

Table 2 Heritabilities (h2 ± SE) and coefficients of additive genetic variation (CVA) of ultrasonic signal characters and development time for male offspring of Achroia grisella reared collectively (expt 1) and individually (expt 2) using NCI half-sib designs

Among males reared individually, heritability based on the sire variance component (paternal heritability, hs2) was significant for peak amplitude only (hs2=0.69). For signal rate, hs2 was low and nonsignificant, but heritability based on the dam component (maternal heritability) was both high (0.45) and significant. This discrepancy suggests a significant nonadditive genetic component or common environmental variation resulting from the collective rearing of each full-sib family (see Falconer, 1981). Heritability for asynchrony interval was low and nonsignificant, although marginally so (P=0.07).

Estimated coefficients of additive genetic variation (CVA) for peak amplitude were high (>19%) and consistent with estimated heritabilities for this character (Table 2). This indicates the potential for significant response to selection by female preferences for louder signals.

Although the heritability (hs2) of signal rate was 0.40 for males reared collectively, CVA estimates were low (<6.2) in both rearing conditions. Mean signal rate did not differ among environments, but total phenotypic variance was greater for males reared individually. Therefore, additive genetic variance was comparable in the two environments, but environmental variance was greater for males reared individually.

Heritability for asynchrony interval was low and not statistically significant. However, CVA estimates were 18.5 and 33.7 for males reared collectively and individually, respectively. Therefore, female preference for longer asynchrony intervals could yield a significant selection response even though heritability estimates were low.

Heritability of development time

Heritability of development time was significant (hs2 = 0.66 ± 0.36) only among males reared individually (Table 2). A previous study using the same stock population, however, found significant heritability of development time among males reared either individually (hs2 = 0.43 ± 0.20) or collectively (hs2 = 0.60 ± 0.22) (Jia & Greenfield, 1997). Causes for the low heritability among males reared collectively are unknown. Taken together, these data suggest that variance in development time may contain a significant genetic component. However, relatively low CVA estimates (Table 2) indicate that observed heritabilities may overstate the responsiveness of development time to selection.

Estimated paternal and maternal heritabilities of development time differed greatly (Table 2). Again, this discrepancy suggests the presence of nonadditive genetic variance or common environmental variation.

Correlations

Phenotypic, genetic, and environmental correlations were estimated for all pairs of signal characters and development time (Table 3). Correlations between signal characters were generally low and nonsignificant, except between signal rate and asynchrony interval. The significant positive phenotypic correlation between peak amplitude and signal rate reported by (Jang et al., 1997) for another A. grisella population was not found in the current study.

Table 3 Phenotypic (P), genetic (G), and environmental (E) correlations (r± SE) among ultrasonic signal characters (PA = peak amplitude†, SR = signal rate, AI = asynchrony interval†) and development time (DT) for male offspring of Achroia grisella reared collectively (expt 1) and individually (expt 2) using NCI half-sib designs

Phenotypic and environmental correlations between signal rate and asynchrony interval were negative for males in both experiments. This relationship is consistent with that reported by Jang et al. (1997) and may result from physiological constraints. The asynchrony interval reflects a delay between movement of the right and left wings on the up- or down-strokes, and it may become difficult to maintain lengthy delays as wing beat frequency increases.

Large differences between the genetic and environmental correlations between signal rate and asynchrony interval were observed. These differences suggest that the effects of genetic and environmental variation may affect these signal characters via different mechanisms (Falconer, 1981). The identity of these mechanisms remains unknown.

The only significant genetic correlation was found between signal rate and asynchrony interval among males reared collectively. This genetic correlation could reflect a physiological linkage between the signal characters, suggesting that signal rate and asynchrony interval may not be independent characters in either a behavioural or a statistical sense.

Low but significant negative phenotypic and environmental correlations between development time and both peak amplitude and signal rate were found for males reared collectively. Males that take longer to develop tended to produce louder signals delivered at faster signal rates, features that are more attractive to females (Jang & Greenfield, 1996). Development time is also correlated with adult weight (r=0.52; P< 0.001; K. Reinhold, unpublished data). Thus, males that take longer to develop are larger and produce more attractive signals.

Discussion

Our finding that the heritabilities of two characters, peak amplitude and signal rate, of ultrasonic advertisement signalling in male A. grisella are high (Table 2) is consistent with earlier work showing that these signal characters are highly repeatable (Jang et al., 1997). These previous measurements, however, showed that the asynchrony interval, also critical to a male’s attractiveness, is highly repeatable (r=0.84). The low heritability of asynchrony interval may reflect a large environmental variance component. Similarly, the markedly different heritabilities of signal rate for males reared collectively and individually (hs2=0.40 and 0.07, respectively) suggest that environmental variance may also be significant for this character.

Environmental variation is assumed to be largely attributable to differences in the extent of larval interactions in the collective and individual rearing conditions. However, the small difference in rearing temperature range between the two experiments may also have contributed to environmental variation.

The significant heritabilities found for peak amplitude and signal rate in A. grisella parallel findings for male signal characters in several other acoustic insects (e.g. Charalambous et al., 1994; Hedrick, 1994; but see Webb & Roff, 1992 for nonsignificant values). Because a significant additive genetic component of variance for the male trait(s) is prerequisite for indirect mechanisms of sexual selection, these findings suggest that either arbitrary or good-genes selection could be responsible for at least part of the observed patterns of female choice in these species. Operation of indirect sexual selection mechanisms also requires significant additive genetic variance for female preference and a genetic correlation between the male trait and female preference. In A. grisella, evidence for heritable variation in female preference has recently been found (Jang, 1997), but examination of associations between male signal characters and female preferences has thus far failed to reveal any positive evidence for genetic correlation (Jang, 1997). The absence of genetic correlation may represent a laboratory artifact, though, as the breeding protocol used in the experimental populations could have reduced linkage disequilibrium and therefore the likelihood of such correlation (see Bakker & Pomiankowski, 1995).

In light of the problems attending measurement of male signal × female preference correlation in the laboratory, can our determinations of the heritability of signal characters be considered as measures applicable to natural A. grisella populations? Traditionally, lower environmental variance is expected in laboratory than in natural populations; indeed, most laboratory studies attempt to minimize such variance. Because of greater phenotypic variability or genotype–environment interactions in the field, laboratory experiments may tend to overestimate heritability within natural populations (e.g. Aspi & Hoikkala, 1993). However, a recent survey (Weigensberg & Roff, 1996) revealed surprising concordance between laboratory and field-based determinations of heritability for various characters in a range of species. We have no reason to believe that A. grisella signal characters represents an exception.

Previous work on A. grisella showed that the male offspring of attractive sires developed faster, ceteris paribus, than those of nonattractive sires under favourable environmental conditions (Jia & Greenfield, 1997). Because rapid development may afford selective advantages if food is limiting or the season is short, the relationship between development time and attractiveness suggested that good-genes sexual selection may sometimes operate in A. grisella. Our current data, however, do not support this suggestion. Although development time was significantly heritable for males reared individually (Table 2) and showed significant phenotypic and environmental correlations with various signal characters for males reared collectively, there were no significant genetic correlations between development time and signal characters (Table 3). Thus, there was no indication that females who mate with males that produce attractive signals will necessarily receive genes for rapid offspring development. Clearly, more work on these relationships is needed under an expanded range of environmental conditions.

Directional selection exerted by female choice on male signal characters has traditionally been expected to restrict additive genetic variance for those characters. However, a recent review (Pomiankowski & Møller, 1995) indicates that additive genetic variance of sexually selected characters is relatively high in many species, and several theoretical treatments (Pomiankowski & Møller, 1995; Rowe & Houle, 1996) have been proffered to account for such values.

Our heritability estimates for several signal characters differed markedly between the two rearing conditions. Because the experiments used insects from the stock population at approximately the same time and included relatively large effective population sizes, it is unlikely that genetic differences accounted for differences in the significance of heritability values. If environmental sources of variance are sufficiently high, heritability may be negligible in some environments even though additive genetic variance is the same. Thus, additive genetic variance may be overwhelmed by environmental variance under some conditions, and selection by females for preferred male signal characters would be ineffective. Such changes in the strength and direction of selection imposed by females may explain, in part, the maintenance of genetic variance in male signals, and, ultimately, indirect female choice.

Although female choice in A. grisella may also result from direct mechanisms such as sensory bias, these mechanisms are not incompatible with indirect arbitrary or good-genes selection (Kirkpatrick, 1996). Moreover, they do not account for the maintenance of genetic variance for signal characters and the persistence of nonattractive males in populations. It is expected that investigations under an extensive range of environmental conditions may ultimately lead to a resolution of this dilemma.