Give ‘til it hurts: trade‐offs between immunity and male reproductive effort in the decorated cricket,Gryllodes sigillatus

Trade‐offs between life‐history variables can be manifested at either the phenotypic or genetic level, with vastly different evolutionary consequences. Here, we examined whether male decorated crickets (Gryllodes sigillatus) from eight inbred lines and the outbred founder population from which they were derived, trade‐off immune effort [lytic activity, phenoloxidase (PO) activity or encapsulation] to produce spermatophylaxes: costly nuptial food gifts essential for successful sperm transfer. Canonical correlation analysis of the outbred population revealed a trade‐off between spermatophylax mass and lytic activity. Analysis of our inbred lines, however, revealed that although PO activity, encapsulation, body mass, spermatophylax mass and ampulla (sperm capsule) mass were all highly heritable, lytic activity was not, and there was, therefore, no negative genetic correlation between lytic activity and spermatophylax mass. Thus, males showed a phenotypic but not a genetic trade‐off between spermatophylax mass and lytic activity, suggesting that this trade‐off is mediated largely by environmental factors.     

A trade-off between sexual signalling and immune function in a natural population of the drumming wolf spider Hygrolycosa rubrofasciata

The field of ecological immunology is ultimately seeking to address the question ‘Why is there variation in immune function?’ Here, we provide experimental evidence that costs of ubiquitous sexual signals are a significant source of variation in immune function. In the mating season, males of the wolf spider Hygrolycosa rubrofasciata drum against dry leaves while wandering around the habitat searching for receptive females. According to a previous study, the male metabolic rate during the drumming increases 22‐fold compared to the resting metabolic rate. In the present study, we examined whether investment in costly courtship drumming decreases male immune function in a wild population of H. rubrofasciata. We induced males to increase their drumming rate by introducing females in proximity. As estimates of male immune function, we used lytic activity and encapsulation rate. Lytic activity estimates the concentration of antimicrobial peptides in haemolymph, which have been shown to play an important role in defence against bacteria, viruses and fungi. Encapsulation is an important defence mechanism against nematodes and insect parasitoids, but it also plays a role in defence against viruses. Our results show that males with nonarbitrarily increased investment in drumming rate had considerably lower lytic activities than control males. Also, there was a tendency for males with nonarbitrarily increased investment in drumming rate to have lower encapsulation rates than control males. This study provides experimental evidence for the first time, to our knowledge, that there are direct immunological costs of sexual signalling in natural populations. Therefore, immunological costs of sexual signals may provide significant phenotypic variation to parasite‐mediated sexual selection.     

THE EVOLUTION OF THRESHOLD TRAITS: A QUANTITATIVE GENETIC ANALYSIS OF THE PHYSIOLOGICAL AND LIFE-HISTORY CORRELATES OF WING DIMORPHISM IN THE SAND CRICKET

Many traits are phenotypically discrete but polygenically determined. Such traits can be understood using the threshold model of quantitative genetics that posits a continuously distributed underlying trait, called the liability, and a threshold of response, individuals above the threshold displaying one morph and individuals below the threshold displaying the alternate morph. For many threshold traits the liability probably consists of a hormone or a suite of hormones. Previous experiments have implicated juvenile hormone esterase (JHE), a degratory enzyme of juvenile hormone, as a physiological determinant of wing dimorphism in the crickets Gryllus rubens and G. firmus. The present study uses a half-sib experiment to measure the heritability of JHE in the last nymphal stadium of G. firmus and its genetic correlation with fecundity, a trait that is itself genetically correlated with wing morph. The phenotypic and genetic parameters are consistent with the hypothesis that JHE is a significant component of the liability. Comparison of sire and dam estimates suggest that nonadditive effects may be important. Two models have been proposed to account for the fitness differences between morphs: the dichotomy model, which assumes that each morph can be characterized by a particular suite of traits, and the continuous model, which assumes that the associated fitness traits are correlated with the liability rather than the morphs themselves. The latter model predicts that the fitness differences will not be constant but change with the morph frequencies. Variation in fecundity and flight muscle histolysis are shown to be more consistent with the continuous model. Data from the present experiment on JHE are inconclusive, but results from a previous selection experiment also suggest that variation in JHE is consistent only with the continuous model.     

Good genes, genetic compatibility and the evolution of polyandry: use of the diallel cross to address competing hypotheses

Genetic benefits can enhance the fitness of polyandrous females through the high intrinsic genetic quality of females' mates or through the interaction between female and male genes. I used a full diallel cross, a quantitative genetics design that involves all possible crosses among a set of genetically homogeneous lines, to determine the mechanism through which polyandrous female decorated crickets (Gryllodes sigillatus) obtain genetic benefits. I measured several traits related to fitness and partitioned the phenotypic variance into components representing the contribution of additive genetic variance ('good genes'), nonadditive genetic variance (genetic compatibility), as well as maternal and paternal effects. The results reveal a significant variance attributable to both nonadditive and additive sources in the measured traits, and their influence depended on which trait was considered. The lack of congruence in sources of phenotypic variance among these fitness-related traits suggests that the evolution and maintenance of polyandry are unlikely to have resulted from one selective influence, but rather are the result of the collective effects of a number of factors.     

Quantitative genetic variation of metabolism in the nymphs of the sand cricket, Gryllus firmus, inferred from an analysis of inbred-lines

Compared with morphological and life history traits, quantitative genetic variation of metabolic and related traits in animals has been poorly studied. We used flow-through VCO(2) respirometry and simultaneous activity measurement on nymphs of the sand cricket (Gryllus firmus) from inbred lines to estimate broad-sense heritability of four metabolic variables. In addition, we measured a number of linear dimensions in the adults from the same inbred lines. There were significant multivariate effects of inbred lines for all traits and broad-sense heritability for physiological traits was 4.5%, 5.2%, 10.3% and 8.5% for average, resting, minimum and maximum CO(2) production in nymphs, respectively. Though the MANOVA indicated significant genetic variation among inbred lines in adult morphology, the broad-sense heritabilities were relatively low ranging from 0-18%. Our results indicate that the heritabilities of metabolic measures are large enough to potentially respond to selection.     

GENETIC CORRELATIONS AMONG TRAITS DETERMINING MIGRATORY TENDENCY IN THE SAND CRICKET,GRYLLUS FIRMUS

Migration by flight is an important component of the life cycles of most insects. The probability that a given insect will migrate by flight is influenced by many factors, most notably the presence or absence of fully-developed wings and functional flight musculature. Considerable variation has also been reported in the flight propensity of fully-winged individuals with functional flight musculature. We test the hypothesis that these components of migratory tendency are genetically correlated in a wing-dimorhic cricket, Gryllus firmus. Flight propensity and condition of the dorsal longitudinal flight muscles (DLM) are examined in fully-winged (LW) crickets from lines selected for increasing and for decreasing %LW, as well as from unselected control lines. Increased %LW is found to be associated with increased flight propensity among individuals with intact DLM, and with retention of functional DLM. The opposite is true for lines selected for decreased %LW. These results indicate both phenotypic and genetic correlations among behavioral, physiological, and morphological traits determining migratory tendency. We propose that these correlations may result from the multifunctional role of juvenile hormone, which has been reported to influence wing development, flight muscle development and degeneration, and flight propensity. Finally, we discuss the potential influence of genetic correlations for migratory traits on the evolution and maintenance of migratory polymorphisms in insects.     

Mapping determinants of variation in energy metabolism, respiration and flight in Drosophila

We employed quantitative trait locus (QTL) mapping to dissect the genetic architecture of a hierarchy of functionally related physiological traits, including metabolic enzyme activity, metabolite storage, metabolic rate, and free-flight performance in recombinant inbred lines of Drosophila melanogaster. We identified QTL underlying variation in glycogen synthase, hexokinase, phosphoglucomutase, and trehalase activity. In each case variation mapped away from the enzyme-encoding loci, indicating that trans-acting regions of the genome are important sources of variation within the metabolic network. Individual QTL associated with variation in metabolic rate and flight performance explained between 9 and 35% of the phenotypic variance. Bayesian QTL analysis identified epistatic effects underlying variation in flight velocity, metabolic rate, glycogen content, and several metabolic enzyme activities. A region on the third chromosome was associated with expression of the glucose-6-phosphate branchpoint enzymes and with metabolic rate and flight performance. These genomic regions are of special interest as they may coordinately regulate components of energy metabolism with effects on whole-organism physiological performance. The complex biochemical network is encoded by an equally complex network of interacting genetic elements with potentially pleiotropic effects. This has important consequences for the evolution of performance traits that depend upon these metabolic networks.     

Fluctuating asymmetry and immune function in a field cricket

Recently, fluctuating asymmetry (FA) of morphological traits has attracted great attention as a short-cut measure of individual quality. Whereas there is some evidence that FA of sexual ornaments is negatively associated with immune function, studies concerning FA and immune function in non-ornamental traits are absent. Here, we tested whether FA of three non-ornamental traits in hind limbs is related to male immune function in a population of the Mediterranean field cricket, Gryllus bimaculatus. As different measures of male immune function, we used encapsulation rate and lytic activity. We found that a composite measure of FA (cFA) was negatively related to encapsulation rate. However, lytic activity was not related to cFA, but there was a tendency that males with higher body mass had higher lytic activity than males with lower body mass. Our results suggest that FA in non-ornamental traits indicates male immunocompetence in G. bimaculatus .     

Genetic correlation between temperature-induced plasticity of life-history traits in a soil arthropod

Temperature is considered one of the most important mediators of phenotypic plasticity in ectotherms. However, the costs and benefits shaping the evolution of different thermal responses are poorly elucidated. One of the possible constraints to phenotypic plasticity is its intrinsic genetic cost, such as genetic linkage or pleiotropy. Genetic coupling of the thermal response curves for different life history traits may significantly affect the evolution of thermal sensitivity in thermally fluctuating environments. We used the collembolan Orchesella cincta to study if there is genetic variation in temperature-induced phenotypic plasticity in life history traits, and if the degree of temperature-induced plasticity is correlated across traits. Egg development rate, juvenile growth rate and egg size of 19 inbred isofemale lines were measured at two temperatures. Our results show that temperature was a highly significant factor for all three traits. Egg development rate and juvenile growth rate increased with increasing temperature, while egg size decreased. Line by temperature interaction was significant for all traits tested; indicating that genetic variation for temperature-induced plasticity existed. The degree of plasticity was significantly positively correlated between egg development rate and growth rate, but plasticity in egg size was not correlated to the other two plasticity traits. The findings suggest that the thermal plasticities of egg development rate and growth rate are partly under the control of the same genes or genetic regions. Hence, evolution of the thermal plasticity of traits cannot be understood in isolation of the response of other traits. If traits have similar and additive effects on fitness, genetic coupling between these traits may well facilitate the evolution of optimal phenotypes. However, for this we need to know the selective forces under field conditions.     

A test of the relationship between cuticular melanism and immune function in wild‐caught Mormon crickets

Cuticular melanism and innate immune parameters can share common physiological pathways in insects, and this functional connection may contribute to the maintenance of insect colour polymorphisms. However, evidence linking colouration and immune function has been equivocal, particularly when tested in wild populations. The present study investigates phenotypic links between colouration and immune function in migratory Mormon crickets (Anabrus simplex, Haldeman), in which juveniles occur in conspicuous colour variants but mature to become uniformly melanic adults. Wild‐caught insects are used to evaluate the relationship between juvenile colouration and three immune parameters: encapsulation ability, lysozyme‐like activity and phenoloxidase activity. As nymphs, brown crickets are better able to encapsulate an inert implant introduced into the haemocoel than green crickets, although the difference is slight and ceases after they all become darkly‐coloured adults. By contrast, adults that develop from brown nymphs have a higher basal phenoloxidase activity than those that develop from green nymphs, regardless of the fact that all adults are brown. Intrinsic factors other than colouration exert larger effects on immunity: males show stronger encapsulation responses but lower phenoloxidase activity than females, suggesting a sex‐specific trade‐off between these two immune parameters, and adults exhibit higher immune function than nymphs. In summary, modest support is found for a correlation between cuticular melanism and increased immune function in wild Mormon crickets. Additional intrinsic factors such as developmental stage and sex appear to interact with colouration and have a more substantial connection to immune function in the wild.     

PATERNAL EFFECTS ON THE EXPRESSION OF A MALE POLYPHENISM

Polyphenic traits are widespread, but compared to other traits, relatively few studies have explored the mechanisms that influence their inheritance. Here we investigated the relative importance of additive, nonadditive genetic, and parental sources of variation in the expression of polyphenic male dimorphism in the mite Rhizoglyphus echinopus, a species in which males are either fighters or scramblers. We established eight inbred lines through eight generations of full‐sibling matings, and then crossed the inbred lines in a partial diallel design. Nymphs were isolated and raised to adulthood with ad libitum food. At adulthood, male morph was recorded for all male offspring. Using a Cockerham–Weir model, we found strong paternal effects for this polyphenic trait that could be either linked to the Y chromosome of males or an indirect genetic effect that is environmentally transmitted. In additional analyses, we were able to corroborate the paternal effects but also detected significant additive effects questioning the Cockerham–Weir analysis. This study reveals the potential importance of paternal effects on the expression of polyphenic traits and sheds light on the complex genetic architecture of these traits.     

Variation in sleep and metabolic function is associated with latitude and average temperature in Drosophila melanogaster

Regulation of sleep and metabolic homeostasis is critical to an animal's survival and under stringent evolutionary pressure. Animals display remarkable diversity in sleep and metabolic phenotypes; however, an understanding of the ecological forces that select for, and maintain, these phenotypic differences remains poorly understood. The fruit fly, Drosophila melanogaster, is a powerful model for investigating the genetic regulation of sleep and metabolic function, and screening in inbred fly lines has led to the identification of novel genetic regulators of sleep. Nevertheless, little is known about the contributions of naturally occurring genetic differences to sleep, metabolic phenotypes, and their relationship with geographic or environmental gradients. Here, we quantified sleep and metabolic phenotypes in 24 D. melanogaster populations collected from diverse geographic localities. These studies reveal remarkable variation in sleep, starvation resistance, and energy stores. We found that increased sleep duration is associated with proximity to the equator and elevated average annual temperature, suggesting that environmental gradients strongly influence natural variation in sleep. Further, we found variation in metabolic regulation of sleep to be associated with free glucose levels, while starvation resistance associates with glycogen and triglyceride stores. Taken together, these findings reveal robust naturally occurring variation in sleep and metabolic traits in D. melanogaster, providing a model to investigate how evolutionary and ecological history modulate these complex traits.     

Genetic analyses of correlated solids, flavor, and health-enhancing traits in onion (Allium cepa L.)

Onion possesses organosulfur compounds and carbohydrates that provide unique flavor and health-enhancing characteristics. Significant phenotypic correlations have been reported among soluble solids content (SSC), total dry matter, pungency, and onion-induced in vitro antiplatelet activity. A genetic map and segregating F3M families derived from a cross between two inbred populations were used to identify and estimate the effects of quantitative trait loci (QTLs) controlling these traits at 30 and 90 days postharvest. In vitro antiplatelet activities among different onion populations were consistent across six human blood donors. Most of the populations showed in vitro antiplatelet activities; however, for some donors, one of the parental lines and two F3M families had pro-aggregatory effects under our experimental conditions. SSC, dry matter, pungency, and in vitro antiplatelet activity showed significant positive phenotypic and genetic correlations. A chromosome region on linkage group E accounted for a significant amount of the phenotypic variation for all of these traits. The correlations among these traits may be due to linkage or pleiotropy of genes controlling solids content. Our results indicate that it will be difficult to develop onion populations with lower pungency and high in vitro antiplatelet activity; however, the strong genetic and phenotypic correlations between high in vitro antiplatelet activity and high SSC are beneficial for the health functionality of onion.     

Efficiency of selection, as measured by single nucleotide polymorphism variation, is dependent on inbreeding rate in Drosophila melanogaster

It is often hypothesized that slow inbreeding causes less inbreeding depression than fast inbreeding at the same absolute level of inbreeding. Possible explanations for this phenomenon include the more efficient purging of deleterious alleles and more efficient selection for heterozygote individuals during slow, when compared with fast, inbreeding. We studied the impact of inbreeding rate on the loss of heterozygosity and on morphological traits in Drosophila melanogaster. We analysed five noninbred control lines, 10 fast inbred lines and 10 slow inbred lines; the inbred lines all had an expected inbreeding coefficient of approximately 0.25. Forty single nucleotide polymorphisms in DNA coding regions were genotyped, and we measured the size and shape of wings and counted the number of sternopleural bristles on the genotyped individuals. We found a significantly higher level of genetic variation in the slow inbred lines than in the fast inbred lines. This higher genetic variation was resulting from a large contribution from a few loci and a smaller effect from several loci. We attributed the increased heterozygosity in the slow inbred lines to the favouring of heterozygous individuals over homozygous individuals by natural selection, either by associative over‐dominance or balancing selection, or a combination of both. Furthermore, we found a significant polynomial correlation between genetic variance and wing size and shape in the fast inbred lines. This was caused by a greater number of homozygous individuals among the fast inbred lines with small, narrow wings, which indicated inbreeding depression. Our results demonstrated that the same amount of inbreeding can have different effects on genetic variance depending on the inbreeding rate, with slow inbreeding leading to higher genetic variance than fast inbreeding. These results increase our understanding of the genetic basis of the common observation that slow inbred lines express less inbreeding depression than fast inbred lines. In addition, this has more general implications for the importance of selection in maintaining genetic variation.     

Courtship song and immune function in the field cricket Gryllus bimaculatus

It has been assumed that sexual ornaments have evolved to reveal males' health and vigour for females. Choosy females may indirectly use ornaments as an indicator of the presence and effectiveness of genes for resistance against parasites. In this study we tested whether females of the Mediterranean field cricket, Gryllus bimaculatus, can use courtship song as a cue for choosing males with high immunocompetence, measured as encapsulation rate of nylon implants and lytic activity of haemolymph. We found that female crickets preferred courtship songs from males with a high encapsulation rate. Female crickets also had a tendency to prefer courtship songs with high tick rate and long high-frequency tick duration. These preferred song components were positively correlated with encapsulation rate, but negatively correlated with lytic activity of the male. In contrast to previous studies of crickets, there was no correlation between male weight and encapsulation rate or lytic activity. There is some evidence in another cricket species that the ability to encapsulate pathogens is heritable. Thus, in light of this study it seems possible that by preferring males according to their courtship song, females might benefit by increasing the parasite resistance of their offspring.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 79 , 503–510.     

INBREEDING: ITS EFFECT ON RESPONSE TO SELECTION FOR PUPAL WEIGHT AND THE HERITABLE VARIANCE IN FITNESS IN THE FLOUR BEETLE,TRIBOLIUM CASTANEUM

We report our studies of the effect of inbreeding on the response to selection for increased pupal weight in the flour beetle, Tribolium castaneum. We also report the effects of inbreeding and selection for pupal weight on the heritable variation in fitness and fitness components. We created replicate and independent inbred lines with F-values of 0.00, 0.375, and 0.672, by 0, 2, and 5 generations, respectively, of brother-sister mating of adult beetles from an outbred stock population. Subsequently, we imposed artificial within-family selection for increased pupal weight in each of 15 inbred lines for eight generations; each line had its own paired, unselected control. We compared the response to selection across the three levels of inbreeding with theoretical expectation, and investigated the effects of inbreeding and selection on fitness variation among families within all 30 selected and control lines. Among-line variation in pupal weight increased with increased inbreeding prior to selection but diminished with directional selection. Inbreeding reduced the realized heritability of pupal weight concordant with quantitative predictions of additive theory. Mean fitness, measured in several ways, declined with inbreeding and declined further with selection. In contrast, the genetic variation for fitness in the inbred and selected lines lines equalled or exceeded that of the outbred controls. Our results suggest that inbreeding and selection may affect traits in different ways depending on the relative amounts of additive and nonadditive genetic variation.     

Quantitative genetics of oviposition behaviour and interactions among oviposition traits in the sand cricket

We performed a quantitative genetic study of oviposition behaviours and oviposition traits in the sand cricket Gryllus firmus. Egg survival in crickets depends on the depth at which they are inserted into the soil with the ovipositor. We examined whether egg depth depends on ovipositor length alone, or on both morphological and behavioural traits associated with oviposition. Heritability estimates were high (h² >0.5) for ovipositor length and small (h²=0.2) for oviposition behaviours. Negative genetic correlations between ovipositor length and some behavioural traits (digging depth and the behavioural component of egg depth) indicated compensation between oviposition traits on egg depth. Because of behavioural compensation, females with different ovipositor lengths subject to stabilizing selection on egg depth could have equal fitnesses. Females laid their eggs deeper, and their eggs were marginally more evenly distributed in dry than in wet sand. This suggests adaptive phenotypic plasticity in laying behaviour, but may also result from physical constraints of the substrate on the insertion of the ovipositor. The absence of significant between-family variation in oviposition traits in response to sand moisture indicates low evolutionary potential for phenotypic plasticity in oviposition traits according to soil moisture. In highly unpredictable environments, females could spread the risk of desiccation by laying eggs at different depths independently of environmental conditions (bet hedging). Our results show significant additive genetic influences on the ability of a female to spread risks as measured by genetic variation in egg distribution, suggesting that a bet-hedging strategy of egg laying has the potential to evolve in this population. Copyright 2002 The Association for the Study of Animal Behaviour. Published by Elsevier Science Ltd. All rights reserved.     

EFFECT OF AN EXPERIMENTAL BOTTLENECK ON MORPHOLOGICAL INTEGRATION IN THE HOUSEFLY

Three measures of multivariate integration were derived from both additive genetic covariance and correlation matrices estimated from parent-offspring covariances to investigate the effect of bottlenecks of different sizes on genetic integration of morphological traits in the housefly, Musca domestica L. Bottleneck lines were initiated with one, four, or 16 pairs of flies sampled from a natural outbred (control) population. Bottlenecks of intermediate size significantly increased the average genetic correlation among traits, resulting in nearly isomorphic variation among all traits in these lines. Single-pair bottlenecks significantly disrupted the trait interrelationships, and the suites of traits identified by principal components of the additive genetic correlation and covariance matrices for the control population were no longer evident in these bottleneck lines. The alteration of the genetic relationships among traits as a result of a bottleneck suggests that nonadditive components of genetic variation affecting these traits were present in the control line. We discuss the implications of nonadditive gene action, particularly epistasis, for speciation via bottlenecks.     

Naturally occurring genetic variation in the age-specific immune response of Drosophila melanogaster

Immunosenescence, the age‐related decline in immune response, is a well‐known consequence of aging. To date, most studies of age‐related changes in immune response focused on the cellular and physiological bases of this decline; we have virtually no understanding of the genetic basis of age‐related changes in the immune system or if indeed such control exists. We used 25 chromosome substitution lines of Drosophila melanogaster derived from a natural population to address three questions: (i) How is the function of the innate immune system influenced by age? (ii) Is there a genetic basis for phenotypic variation in immune response at different ages? (iii) Is there a genetic basis for differences in the way that age influences the immune function? Virgin females from each line were assayed for immune response using clearance of infection with Escherichia coli at 1 and 4 weeks of age. We found significant genetic variation among lines in immune response at each age. Unexpectedly, when averaged across all lines, the immune response actually improved with age. However, there was significant variation in the effect of age on immune response with 11 lines showing improvement, nine lines showing no change and five exhibiting a decline with age. There was no genetic correlation of immune response across ages suggesting that different loci contribute to variation in immune response at each age. The genetic component of the variation in immune response increased with age, a pattern predicted by the mutation accumulation model of senescence. However, this increase in variation resulted in part from the improvement of the immune response in some lines with age. Thus the observed changes in genetic variation in immune function with age are not entirely explained by the mutation accumulation model.     

Phenotypic variation and covariation indicate high evolvability of acoustic communication in crickets

Studying the genetic architecture of sexual traits provides insight into the rate and direction at which traits can respond to selection. Traits associated with few loci and limited genetic and phenotypic constraints tend to evolve at high rates typically observed for secondary sexual characters. Here, we examined the genetic architecture of song traits and female song preferences in the field crickets Gryllus rubens and Gryllus texensis. Song and preference data were collected from both species and interspecific F1 and F2 hybrids. We first analysed phenotypic variation to examine interspecific differentiation and trait distributions in parental and hybrid generations. Then, the relative contribution of additive and additive‐dominance variation was estimated. Finally, phenotypic variance–covariance ( P ) matrices were estimated to evaluate the multivariate phenotype available for selection. Song traits and preferences had unimodal trait distributions, and hybrid offspring were intermediate with respect to the parents. We uncovered additive and dominance variation in song traits and preferences. For two song traits, we found evidence for X‐linked inheritance. On the one hand, the observed genetic architecture does not suggest rapid divergence, although sex linkage may have allowed for somewhat higher evolutionary rates. On the other hand, P matrices revealed that multivariate variation in song traits aligned with major dimensions in song preferences, suggesting a strong selection response. We also found strong covariance between the main traits that are sexually selected and traits that are not directly selected by females, providing an explanation for the striking multivariate divergence in male calling songs despite limited divergence in female preferences.