Correlated evolution between repertoire size and song plasticity predicts that sexual selection on song promotes open-ended learning

  1. Cristina M Robinson
  2. Kate T Snyder
  3. Nicole Creanza  Is a corresponding author
  1. Vanderbilt University, United States

Abstract

Some oscine songbird species modify their songs throughout their lives (‘adult song plasticity’ or ‘open-ended learning’), while others crystallize their songs around sexual maturity. It remains unknown whether the strength of sexual selection on song characteristics, such as repertoire size, affects adult song plasticity, or whether adult song plasticity affects song evolution. Here, we compiled data about song plasticity, song characteristics, and mating system and then examined evolutionary interactions between these traits. Across 67 species, we found that lineages with adult song plasticity show directional evolution toward increased syllable and song repertoires, while several other song characteristics evolved faster, but in a non-directional manner. Song plasticity appears to drive bi-directional transitions between monogamous and polygynous social mating systems. Notably, our analysis of correlated evolution suggests that extreme syllable and song repertoire sizes drive the evolution of adult song plasticity or stability, providing novel evidence that sexual selection may indirectly influence open- versus closed-ended learning.

https://doi.org/10.7554/eLife.44454.001

eLife digest

Every morning, people all over the world are greeted by the sound of songbirds singing to attract mates and defend their homes. Each type of songbird has its own unique song that it learns from other birds in the same species. Maintaining these signature songs is important, because songbirds that sing the wrong tune are unlikely to succeed in breeding. This type of vocal learning is rare in the animal kingdom and is similar to how humans learn to speak.

The time it takes for songbirds to learn their unique song varies between species: some species pick up their song within a few months after hatching, while others may continue learning for years or even the rest of their lives. It remains unclear, however, why certain species spend a longer time learning, and how changing this length of time affects the song they sing. One possibility is that differences in the amount of time spent learning evolved from female songbirds preferring mates with more impressive songs.

To address this possibility, Robinson et al. used previously published data to compare the song characteristics, learning periods, and mating strategies of 67 songbird species. This confirmed that longer learning windows allowed songbirds to develop a larger repertoire of songs with more unique sounds. Robinson et al. also found songbirds that learn elaborate songs over a short period of time quickly evolve either simpler songs or longer phases of learning. Meanwhile, songbirds that learn simple songs over a longer period tended to evolve more elaborate songs or shorten the amount of time they spent learning. Furthermore, species with longer learning were more likely to switch between having one mating partner and multiple partners over several generations.

These findings suggest that by choosing mates with elaborate songs and larger repertoires, female songbirds can end up favoring the evolution of longer learning. Studying the conditions that led to longer or shorter learning in songbirds could help scientists understand why some things, like human language, are easier to learn early in life.

https://doi.org/10.7554/eLife.44454.002

Introduction

Song is a learned behavior with a complex evolutionary history in the oscine songbirds. Birds’ songs have multiple functions, including species recognition, territory defense, and mate attraction (Catchpole and Slater, 2003). All species studied in this expansive clade have certain life stages during which they are more likely to learn and acquire songs, termed sensitive periods (Brenowitz and Beecher, 2005; Marler, 1990; Marler and Peters, 1987; Rauschecker and Marler, 1987). In some species, learning is restricted to a short sensitive period early in life, also called a ‘critical period’ (e.g. ~day 25–90 in zebra finches), after which no new song elements are acquired (Böhner, 1990; Immelman, 1969; Nottebohm, 1984). Other species appear to delay song crystallization until some time in adulthood (Dowsett-Lemaire, 1979; Kipper and Kiefer, 2010; Martens and Kessler, 2000); for example, chipping sparrows appear to have a second sensitive period immediately after their first migration, following which their song crystallizes (Liu and Kroodsma, 2006; Liu and Nottebohm, 2007). Still other species can continue to acquire new syllables or songs throughout their lives (Adret-Hausberger et al., 2010; Espmark and Lampe, 1993; Gil et al., 2001; Hausberger et al., 1991; Mountjoy and Lemon, 1995; Price and Yuan, 2011). Typically, this spectrum of variation in the timing of the sensitive period is simplified into a dichotomy of ‘open-ended learning’ and ‘closed-ended learning.’ While these temporal differences in the song-learning window have been studied for decades, it is unknown how they interact with the evolution of song itself. Previous hypotheses have suggested that seasonal factors, such as environmental variation and breeding season length, play a role in shaping adult song learning (Nottebohm et al., 1986; Smith et al., 1997; Tramontin et al., 2001). However, evidence from a small-scale comparative analysis suggests that a longer learning window in a species may be associated with larger average syllable repertoire sizes (Creanza et al., 2016).

Birdsong is composed of both culturally and genetically inherited features, any of which may be subject to evolutionary pressures. Two key modes of selection on song might act in conjunction: female choice can favor certain song characteristics, such as superior repertoire size, learning quality, or song performance (Beecher and Brenowitz, 2005; Gil and Gahr, 2002; Searcy and Marler, 1984; Searcy and Andersson, 1986), while the inherent metabolic cost of neuroplasticity should theoretically favor a shorter song-learning window and thus reduce the a for a bird to alter its song in adulthood (Garamszegi and Eens, 2004; Nottebohm et al., 1986; Tramontin and Brenowitz, 2000). Therefore, while learning in adulthood or elongated sensitive periods have not been shown to be directly under positive selection (sexual or natural) or to play an explicit role in female preferences, sexual selection acting on certain song features could indirectly favor longer or shorter song-learning windows. However, this theorized connection between sexual selection and adult learning hinges on establishing the evolutionary relationship between song as the target of sexual selection and the neurobehavioral phenotype of song learning, which has not yet been done.

Furthermore, sexual selection is hypothesized to be amplified in species with polygynous social mating systems or high rates of extra-pair paternity (EPP). A recent large-scale study found that polygyny drives faster, but non-directional, evolution of syllable repertoire size, and that syllable repertoire size is negatively correlated with the rate of EPP (Snyder and Creanza, 2019). Because of the evolutionary links between these mating strategies and song features, higher rates of EPP and polygyny could potentially have an effect on learning windows. We therefore investigate the evolutionary relationship between the critical learning period and mating strategies.

Here, we take a comparative, computational approach to the evolutionary history of open- and closed-ended song learning. We mined the literature for longitudinal field and laboratory observations to classify species as exhibiting ‘adult song stability’ or ‘adult song plasticity’. This classification is a quantifiable proxy for closed- and open-ended learning, as the true length of the song-learning window is difficult to assess directly in nature; ultimately, we obtained data for the classification of 67 species. For these species, we also compiled a database of seven species-level song characteristics that can represent either song complexity (syllable repertoire, syllables per song, and song repertoire) or song performance (song duration, inter-song interval, song rate, and song continuity). We then performed phylogenetically controlled analyses to evaluate the evolution of song and mating strategies alongside the relative plasticity of song over time. We find that adult song plasticity has evolved numerous times in bird species. Further, we find evidence of correlated evolution between adult song stability and plasticity and social mating system, with shifts in social mating system occurring more rapidly in lineages with adult song plasticity. In addition, we find a significant evolutionary pattern: species with plastic songs generally have larger repertoires than species with stable songs. Specifically, the evolution of larger syllable and song repertoires appears to drive an evolutionary transition toward open-ended learning.

Results

In this study, we analyzed the length of the song-learning window using three classification schemes. First, we classified 67 species as having either adult song plasticity (i.e. those that change their syllable repertoire after their first breeding season ends) or adult song stability. Second, when possible, we reclassified the species into three categories: 1) those that stop changing their song before their first breeding season, 2) those that modify their songs during their first breeding season but not after, and 3) those that modify their songs after their first breeding season has ended. This three-state classification was possible for 59 out of the 67 species. Third, for the same 59 species, we used reported estimates of the ages at which song stabilizes to create a continuous measure of the song-learning window. There were exceedingly few studies that examined song changes after the second breeding season, so our continuous metric ranges from 0 to 2 years. We then examined the evolutionary patterns of adult song stability and plasticity, as well as their interactions with species-level song characteristics and mating behaviors.

How has adult song plasticity evolved across clades?

We were first interested in examining the rate of evolution of adult song stability versus adult song plasticity, as well as when and where evolutionary transitions in these traits occurred on a phylogenetic tree (Jetz et al., 2012) using ancestral state reconstruction. As with any reconstruction of evolutionary history, these simulations cannot exactly predict the ancestral states but aim to approximate them. Furthermore, we note that only a subset of oscine families were represented in our analysis, and most of the early branching lineages that would be required to assess the ancestral state for all oscine species were not included in our dataset. Ultimately, we could not make a conclusion about whether the last common ancestor for the species included in this study had adult song plasticity, but our results hint that there might have been several early transitions in this trait, leading to clades that predominantly have adult song stability or plasticity, coupled with a number of more recent transitions (see pie charts in Figure 1A for the predicted likelihood of each state at each node). We found that a model allowing the transition rate from song stability to plasticity to be different from the transition rate from plasticity to stability (all-rates-different model [ARD]) did not fit the data significantly better than a simpler model allowing for only one rate of transition back and forth between song stability and plasticity (equal rates model [ER]) (LogLikelihoodER = −38.22, LogLikelihoodARD = −38.21, p=0.87). At least 14 transitions were required to explain the current binary song-stability states of our subset of bird species. Explaining the distribution of song plasticity in our subset of species most parsimoniously requires at least nine transitions to adult song plasticity if the last common ancestor was song-stable and seven transitions to song stability if the common ancestor was song-plastic (Figure 1—figure supplement 1).

Figure 1 with 8 supplements see all
Syllable repertoire size is larger in species with adult song plasticity even when controlling for phylogenetic relationships.

These phylogenies show the calculated evolution of natural-log transformed syllable repertoire size and either (A) stable and plastic song stability states or (B) early song-stable, delayed song-stable, and song-plastic states. Dots at the tips of branches represent the current song-stability state. Pie charts represent the likelihood that the common ancestor at that node was in each song-stability state. Dark purple colors represent small syllable repertoires while white represents large repertoires. For the sake of visualization, the color range was truncated based on the distribution of the data, such that the lowest value was the 25th percentile minus the range of the 25th to 50th percentile and the highest value was the 75th percentile plus the range of the 50th to 75th percentile. See Table 1 for PhylANOVA results.

https://doi.org/10.7554/eLife.44454.003
Figure 1—source data 1

Predicted values for internal nodes for each trait.

There are three plots with associated tables for each song characteristic: The first plot and table set gives the values when the binary categorization of song stability (song-stable and song-plastic) is used. The second plot and table set gives the values when the ternary categorization of early song-stable, delayed song-stable, and song-plastic is used. The final plot and table set gives values when the continuous categorization is used. In the table, ‘Node’ corresponds to the numbered boxes on the tree nodes. ‘State Likelihood’ gives the probability that the common ancestor at a node was in a given leaning state when discrete categorizations are used. S = song stable, p=song plastic, E = early song-stable, D = delayed song-stable. ‘Length of Plasticity’ predicts the value along a continuous spectrum. ‘Trait’ is the predicted value of the song trait being examined in the associated tree for each internal node.

https://doi.org/10.7554/eLife.44454.012

Which song traits differ between species with stable versus plastic songs?

We next tested whether song characteristics were affected by the length of the song-learning window on an evolutionary scale. Intuitively, it makes sense that a species that has a longer time-window to learn might be able to accumulate a larger repertoire. Indeed, this relationship is consistent with the pattern of song stability and repertoire size in several clades, such as the Phylloscopus species (Figure 2). However, many individual species do not follow this prediction: for example, Acrocephalus palustris appears to learn a large repertoire in a single year (Dowsett-Lemaire, 1979), and Philesturnus rufusater modifies its song for multiple years but maintains a small repertoire (Jenkins, 1978). Further, numerous species with adult song plasticity do not significantly increase their repertoire sizes over time (Eriksen et al., 2011; Galeotti et al., 2001; Garamszegi et al., 2005; Nicholson et al., 2007). Thus, an evolutionary link between adult song plasticity and larger repertoire sizes cannot be assumed. Using a phylogenetically controlled ANOVA (Garland et al., 1993; Revell, 2012), we found that species with adult song plasticity did possess significantly larger syllable repertoires than species with adult song stability (Figures 1A and 3A, Table 1). This concurs with a previous analysis using a smaller dataset (Creanza et al., 2016). Similarly, we found that song-plastic species had significantly larger song repertoires than song-stable species (Figure 3B, Figure 1—figure supplement 2, and Table 1).

Table 1
PhylANOVA results for all song traits when birds are divided into species with adult song stability or adult song plasticity.

Song traits are sorted from most to least significant. Song-stable and song-plastic columns show mean values of each log-transformed song trait. Corrected α indicates the threshold for significance with the Holm-Bonferroni correction.

https://doi.org/10.7554/eLife.44454.013
Song traitSong-stableSong-plasticF-ValueCorrected αp-value
Syllable repertoire1.88073.94641.50640.0071<0.001*
Song repertoire1.10553.868833.83340.0083<0.001*
Syllables/song1.25562.29629.26580.010.094
Duration0.77361.29272.07830.01250.42
Continuity-1.3453-1.02862.15370.01670.474
Interval1.60751.2181.38790.0250.567
Song rate1.89692.09710.60790.050.713
  1. *Denotes traits with significantly different groups.

Comparison of song stability state and song examples in Phylloscopus species.

P. collybita, a species with adult song stability, has a smaller syllable repertoire size than P. trochilus and P. fuscastus, two species with adult song plasticity. Colors of branches and nodes correspond with Figure 1. Sonograms were generated from recordings obtained from xeno-canto.org: XC340281 recorded by Tom Wulf (P. fuscatus), XC414221 recorded by Frank Lambert (P. collybita), and XC402265 recorded by Hans Matheve (P. trochilus). Sonograms are used only to demonstrate comparative repertoire size from one individual for each species and were stretched horizontally to fit the allotted space.

https://doi.org/10.7554/eLife.44454.014
Distribution of repertoire sizes in species with different song stability states.

(A) shows the distribution of syllable repertoires and (B) shows the distribution of song repertoires when species are broken into two groups based on song stability. (C) shows the distribution of syllable repertoires when species are broken into three groups based on song stability. Boxes indicate the 25th, 50th, and 75th percentile. The lower whisker is either the minimum value or the 25th percentile minus 1.5 times the interquartile range, whichever was larger. The upper whisker is either the maximum value or the 75th percentile plus 1.5 times the interquartile range, whichever was smaller. Dots are the raw values as a scatter plot. Solid dots are within the range of the box and whiskers, while open dots are outliers. (A-B) Species with adult song plasticity had significantly larger syllable and song repertoires. See Table 1 for full PhylANOVA results. (C) Species with adult song plasticity had significantly larger syllable repertoires than early song stable and delayed song stable species, but there was no significant difference between early and delayed song stable species (p=0.659). See Tables 2 and 3 for full PhylANOVA results. (D) shows the continuous relationship between syllable repertoire size and years spent learning, where song plasticity is truncated at 2 years due to lack of data from subsequent years.

https://doi.org/10.7554/eLife.44454.015

There were no significant differences between song-plastic and song-stable species for the other song characteristics that we tested: syllables per song, inter-song interval, song duration, song rate (calculated as 60/(interval + duration)), or song continuity (calculated as duration/(duration + interval)) (Table 1, Figure 1—figure supplements 37). When we used the classification scheme with three states, we could only test for differences in syllable repertoire, song repertoire, and syllables per song between groups, as there were very few early song-stable species for which we had data on the other song traits. We found no significant differences between early song-stable and delayed song-stable species for any tested traits, but both of these groups had significantly smaller syllable and song repertoires compared to song-plastic species (Figures 1B and 3C, Tables 2 and 3). When performing a phylogenetic generalized least squares (PGLS) analysis using continuous estimates of the duration of song plasticity, we found similar results; both syllable repertoire and song repertoire were correlated with duration of song plasticity, such that repertoire size increased with the song-plasticity duration (Figure 3D, Table 4).

Table 2
PhylANOVA results for all song traits when birds are divided into early song stability, delayed song stability, and song plasticity.

Song traits are sorted from most to least significant. Early, delayed, and plastic columns show mean values of each log-transformed song trait. Corrected α indicates the threshold for significance with the Holm-Bonferroni correction.

https://doi.org/10.7554/eLife.44454.016
Song traitEarlyDelayedPlasticF-ValueCorrected αp-value
Syllable repertoire1.64362.00623.94617.10990.00710.003*
Song repertoire0.67881.48193.868812.880.00830.011*
Syllables/song1.28521.24672.29623.68770.010.252
  1. *Denotes traits with significantly different groups.

Table 3
Post-hoc pairwise phylANOVA tests for significant song traits when birds are divided into early song stability, delayed song stability, and song plasticity.
https://doi.org/10.7554/eLife.44454.017
Song traitState 1State 2T-Valuep-value
Syllable repertoirePlasticDelayed4.89950.012*
Syllable repertoireEarlyPlastic4.60910.003*
Syllable repertoireEarlyDelayed0.68720.659
Song repertoirePlasticDelayed4.02680.044*
Song repertoireEarlyPlastic4.30740.015*
Song repertoireEarlyDelayed1.04440.55
  1. *Denotes traits with significantly different groups.

Table 4
Results of PGLS analysis between song characteristics and continuous song stability.

Test performed on the natural-log scaled values of song characteristics. λ is the value by which off-diagonal elements in the Brownian motion model are multiplied to make the correlation structure. Corrected α indicates the threshold for significance with the Holm-Bonferroni correction. Song traits are sorted from most to least significant.

https://doi.org/10.7554/eLife.44454.018
Song traitSlopeStd errorλT-ValueCorrected αp-value
Syllable repertoire0.90670.24490.89133.70210.0071<0.001*
Song repertoire1.10130.31230.83163.52630.0083<0.001*
Syllables/song0.37010.22240.46991.66420.010.1029
Interval0.42210.26460.88231.59530.01250.1215
Continuity-0.21350.14390.8832-1.48380.01670.1486
Duration0.37020.25691.01631.4410.0250.1578
Song rate-0.21130.250.7307-0.84530.050.4048
  1. *Denotes significant slopes.

Do song characteristics evolve at different rates in song-stable or song-plastic species?

Our result that species with adult song plasticity had significantly larger syllable and song repertoires raised the question of whether song stability versus plasticity also affected the rate of evolution for any of the song characteristics. To examine this possibility, we used the Brownie algorithm (O'Meara et al., 2006), which tests whether a model with two rates of evolution for each song characteristic—one rate for ancestral periods of song stability and another rate for song plasticity—fits the data significantly better than a model that allows for only a single rate of evolution of each song characteristic regardless of the ancestral states of song stability. Each calculation of the two-rate model is based on one stochastic projection of the ancestral traits across the phylogenetic tree, so we generated 1300 different stochastic simulation maps to use with Brownie. We plotted the distribution of potential rates (Figures 45) and compared the average log likelihood of the two-rate models to the log likelihood of the one-rate model (Table 5).

Table 5
Brownie results for song traits when birds are divided into species with adult song stability or adult song plasticity.

Rate columns show mean log likelihood. Song traits are sorted from most to least significant.

https://doi.org/10.7554/eLife.44454.019
Song traitOne rateTwo ratesp-value
Syllables/song-110.6482-100.7673<0.001*
Song rate-43.4397-38.49380.002*
Interval-45.2842-40.50040.002*
Duration-71.2042-66.31220.002*
Continuity-25.6471-24.72850.175
Syllable repertoire-120.2983-120.06950.499
Song repertoire-113.5829-113.37060.515
  1. *Denotes traits where the more complex model fit the data significantly better than the simpler model.

Distributions of rates for natural-log transformed song traits related to complexity.

The boxes at the top illustrate how we grouped the species for each model. Column 1 (A-C): Blue traces are song-stable, while red traces are song-plastic. Column 2 (D-F): Blue traces are early song-stable, purple traces are delayed song stable, and red traces are song-plastic. Column 3 (G-H): Blue traces are early song-stable, while red traces are delayed song-stable and song-plastic combined. The black line shows the rate value for the one-rate model in all columns. Asterisks indicate that the rate of evolution of that song characteristic significantly differed between groups. Lowercase zeta (ζ) the multi-rate model that best fit the data while using the fewest number of rates. In the case of syllable repertoire, the multi-rate models were not significantly better than the one-rate model. See Tables 59 for chi-square test results.

https://doi.org/10.7554/eLife.44454.020
Distributions of rates of evolution for natural-log transformed song traits related to performance.

Blue traces song-stable, while red traces are song-plastic. The black line shows the rate for the one-rate model. Asterisks indicate that the rate of evolution of that song characteristic significantly differed between song-stable and song-plastic lineages. See Table 5 for chi-square test results.

https://doi.org/10.7554/eLife.44454.021

We found that allowing for two different rates of song trait evolution depending on song stability or plasticity did not lead to a significantly better fit model than using one Brownian-motion rate for either syllable repertoire size or song repertoire size, even though syllable repertoires and song repertoires were both significantly larger in species with adult song plasticity (Figure 4A and B and Table 5). In contrast, the two-rate model led to a significantly better fit for syllables per song, song rate, inter-song interval, and song duration (Figures 4C and 5 and Table 5), indicating that evolution of these song characteristics was faster in song-plastic lineages (Figures 4C and 5, red traces).

We repeated this analysis with the three-state categorization for syllable repertoire, song repertoire, and syllables per song; for other song characteristics, we did not have enough species in the early song stability group. We found that the three-rate model was significantly better than the one-rate model for syllables per song and song repertoire, but not for syllable repertoire (Figure 4D–F and Table 6). However, the three-rate model was only significantly better than the two-rate model for song repertoire (Table 7). Thus, the two-rate model sufficiently approximated the evolution of syllables per song. We noticed that for both song repertoire and syllable repertoire, the rate of evolution in delayed song-stable lineages (purple traces in Figure 4D,E) was very similar to the rate in song-plastic lineages (corresponding red traces). We tested one more set of models where we combined delayed song-stable species with song-plastic species to create a ‘longer learning’ group, while early song-stable species were assigned to a ‘shorter learning’ group. For this comparison of shorter versus longer learning, the two-rate model was significantly better than the one-rate model for song repertoire and trending in that direction for syllable repertoire (Figure 4G,H and Table 8). The three-rate model was not significantly better than the longer/shorter-learning two-rate model for either syllable or song repertoire (Table 9). Taken together with our phylANOVA results, this pattern suggests that species with early song stability evolve their song repertoires and potentially their syllable repertoires at a slower rate than delayed song-stable and song-plastic species; however, only song-plastic species directionally evolve towards larger song and syllable repertoires.

Table 6
Brownie results for song traits when birds are divided into early song stability, delayed song stability, and song plasticity.

Rate columns show mean log likelihood. Song traits are sorted from most to least significant.

https://doi.org/10.7554/eLife.44454.022
Song traitOne rateThree ratesp-value
Syllables/song-97.8349-86.3206<0.001*
Song repertoire-100.812-97.76470.014*
Syllable repertoire-107.3206-105.58950.063
  1. *Denotes traits where the more complex model fit the data significantly better than the simpler model.

Table 7
Brownie results for song traits when birds are divided into either song stability (early plus delayed) and song plasticity (Two Rates) or early song stability, delayed song stability, and song plasticity (Three Rates).

Rate columns show mean log likelihood. Song traits are sorted from most to least significant.

https://doi.org/10.7554/eLife.44454.023
Song traitTwo ratesThree ratesp-value
Song repertoire-100.691-97.71480.015*
Syllable repertoire-107.1332-105.55320.075
Syllables/song-86.3125-86.34471
  1. *Denotes traits where the more complex model fit the data significantly better than the simpler model.

Table 8
Brownie results for song traits when birds are divided into shorter learning (early song stability) and longer learning (delayed song stability plus song plasticity).

Rate columns show mean log likelihood. Song traits are sorted from most to least significant.

https://doi.org/10.7554/eLife.44454.024
Song traitOne rateTwo ratesp-value
Song repertoire-100.812-97.99180.018*
Syllable repertoire-107.3206-105.84880.086
  1. *Denotes traits where the more complex model fit the data significantly better than the simpler model.

Table 9
Brownie results for song traits when birds are divided into either shorter learning (early song stability) and longer learning (delayed song stability plus song plasticity) (Two Rates) or early song stability, delayed song stability, and song plasticity (Three Rates).

Rate columns show mean log likelihood. Song traits are sorted from most to least significant.

https://doi.org/10.7554/eLife.44454.025
Song traitTwo ratesThree ratesp-value
Syllable repertoire-105.8156-105.55320.469
Song repertoire-97.9372-97.71480.505
  1. *Denotes traits where the more complex model fit the data significantly better than the simpler model.

Is song stability or plasticity influenced by the evolution of song traits and mating strategies, or vice versa?

While the Brownie algorithm tested whether adult song plasticity affected the rate of evolution for song characteristics, it did not address whether the opposite might be true. We used BayesTraits (Pagel, 1994; Pagel and Meade, 2006) to test whether the rate and order of evolutionary transitions in one trait is dependent on the state of another trait. Because the song features were continuous variables, we binarized them by setting a series of threshold values to delineate ‘low’ and ‘high’ categories, using each observed song feature value as a threshold in turn. We then tested whether there was correlated evolution between the binary classifications of adult song plasticity versus stability and each of the seven song characteristics.

In the lowest third of syllable repertoire thresholds, adult song plasticity with small syllable repertoires was an evolutionarily unstable state, with rapid transitions primarily toward a song-stable state and secondarily toward larger repertoires (82% of runs significant in this range, Figure 6). In the middle third of syllable repertoire thresholds, song stability with smaller syllable repertoires is an evolutionarily stable attractor state, with high rates of transition observed from large to small syllable repertoires in song-stable species and from plasticity to stability with a small syllable repertoire. These rate differences are highly significant (100% of runs significant in this range). In the highest third of syllable repertoire thresholds, adult song stability with a large syllable repertoire is an evolutionarily unstable state, transitioning primarily toward adult song plasticity (86% of runs significant in this range, Figure 6). We found similar trends when using two, four, and five bins for the song characteristic threshold values, with subtle differences. When using four or five bins, we still observe that song stability with larger syllable repertoires is an unstable combination. However, for the highest bin of threshold values, the transition rates are faster when changing to song plasticity, whereas for the second-highest bin, we observe faster transition rates toward repertoire size increases (Figure 6—source data 1). To rule out the possibility that syllable repertoire size evolution is faster in species with larger repertoire sizes regardless of learning program, we tested the rates of evolution of syllable repertoire size across monophyletic species pairs in our dataset. We found that lineages with larger syllable repertoire sizes do not systematically undergo faster or slower syllable repertoire size evolution (Figure 1—figure supplement 8).

Figure 6 with 6 supplements see all
Analysis of correlated evolution between adult song stability/plasticity and syllable repertoire size.

We repeated the BayesTraits analysis using each value of the continuous song trait as the threshold delineating the larger and smaller syllable repertoire groups. We performed a total of 100 runs per threshold. We pooled the results of all the runs into three groups based on whether the threshold was in the lowest, middle, or highest third of the unique trait values. Within these groups, we computed the mean percentage of runs that were significant at p<0.05 at each threshold. (A-C) Rate of transition plots when the lowest (red), middle (yellow), and highest (blue) thirds of the unique syllable repertoire values in the dataset were used as the threshold. Rates are the average across all runs when the threshold denoting small/large repertoire sizes was defined as each value within each segment. Arrows are labeled with the mean rate and the 95% confidence interval. Arrow weights are scaled to the mean rate values. (D) p-values from the 100 runs per threshold, plotted against threshold. Colored bars denote low, middle, and high threshold segments. Blue line denotes p=0.05.

https://doi.org/10.7554/eLife.44454.026
Figure 6—source data 1

Mean transition rates from BayesTraits analysis between song stability and song features averaged over 2, 4, and 5 bins.

On page 1, the top two panels depict the mean rates from all runs across the values of syllable repertoire in the lower and upper half of thresholds, respectively. The bottom panel shows the p-values from the 100 runs per threshold, plotted against threshold. Colored bars denote low and high threshold segments. The blue line denotes p=0.05. On page 2, the top four panels depict the mean rates from all runs across the values of syllable repertoire in the first (top left), second (top right), third (middle left), and fourth (middle right) quartiles of the thresholds. The bottom panel follows the pattern of the bottom panel from page 1. On page 3, the five panels (top left and right, middle left and right, and bottom left) depict the mean rates from all runs across the values of syllable repertoire in the five quintiles of the of the thresholds, respectively. The bottom panel follows the pattern of the bottom panels from pages 1 and 2. The three-page pattern repeats for each song feature.

https://doi.org/10.7554/eLife.44454.027
Figure 6—source data 2

BayesTraits analysis of syllable repertoire and learning window, jackknifed across families.

Labeling the same as in Figure 6 and Figure 6—figure supplement 1.

https://doi.org/10.7554/eLife.44454.034
Figure 6—source data 3

BayesTraits analysis of syllables per song and learning window, jackknifed across families.

Labeling the same as in Figure 6 and Figure 6—figure supplement 1.

https://doi.org/10.7554/eLife.44454.035

At low song repertoire threshold values, song plasticity with a small repertoire is an evolutionarily unstable state; there is rapid transition away from this combination, predominantly trending toward song stability but also transitioning secondarily to a larger song repertoire, with very high significance (100% of runs significant in this range). At moderate song repertoire thresholds, the highest rate is observed for species with small repertoires transitioning from song plasticity to song stability, also with very high significance (100% of runs significant in this range). At high song repertoire thresholds, the primary shift is from song stability to song plasticity in species with large song repertoires (89% of runs significant in this range) (Figure 7). When analyzing the results using five bins, transitions in the upper range of song repertoire values becomes more nuanced; in the highest bin, song stability with a larger song repertoire is very unstable, but is relatively stable in the second-highest bin. In the lower four bins, the dominant transition is from plastic to stable song with smaller song repertoires (Figure 6—source data 1). The results for syllables per song show some general trends that are complicated by the strong effect of the mimid species (see results of jackknife analysis below). We did not find evidence for correlated evolution between adult song stability or plasticity and any of the other song characteristics (Figure 6—figure supplements 15).

Analysis of correlated evolution between adult song stability/plasticity and song repertoire size.

Labeling is the same as in Figure 6.

https://doi.org/10.7554/eLife.44454.036
Figure 7—source data 1

BayesTraits analysis of song repertoire and learning window, jackknifed across families.

Labeling the same as in Figure 6 and Figure 6—figure supplement 2.

https://doi.org/10.7554/eLife.44454.037

In addition, it has been proposed that polygyny and extra-pair paternity (EPP) may increase sexual selection pressures on sexually selected traits, including song (Emlen and Oring, 1977; Payne, 1984), and increased selection pressure due to polygyny was theorized to accelerate the evolution of song learning in a mathematical model (Aoki, 1989). We tested for correlated evolution between adult song plasticity versus stability and both social mating system (polygyny vs. monogamy) and extra-pair paternity (low vs. high EPP), with the caveat that many species in our dataset lacked mating behavior classifications (57 species with social mating system data, 41 with EPP data). We did not find evidence for correlated evolution between song stability and EPP. There was, however, evidence for correlated evolution between polygynous/monogamous mating systems and song plasticity (100% of runs significant), with elevated rates of transition between polygyny and monogamy in the song-plastic state (Figure 6—figure supplement 6).

Do our results depend on the particular values or families included in analyses?

In many cases, there were multiple studies that gave different estimates for a given song trait in one species, so we used the median values across studies for our main analysis. To test whether our results depended on the particular values we used, we repeated the PhylANOVA and Brownie analysis using the either the maximum or minimum values reported in the literature. This did not significantly alter our PhylANOVA results for any song feature when species were split into song-stable and song-plastic groups (Supplementary file 1 Table S1), though when the early song-stable, delayed song-stable, and song-plastic dataset was used, using the maximum values for song repertoire led to non-significant results (Supplementary file 1 Table S2-3). When species were split into those with adult song plasticity and adult song stability, the Brownie analysis suggested that syllables per song was evolving significantly faster in song-plastic lineages when we used the median value, but using the minimum values led to non-significant results (Supplementary file 1 Table S4). When species were split into three states (early song stability, delayed song stability, and song plasticity), we found that the three-rate model was not better than the one-rate model when the minimum values for song repertoire were used, while the three-rate models for syllable repertoire became significant when either the maximum or the minimum values were used (Supplementary file 1 Table S5).

In our dataset, 24 songbird families were represented by 1 to 11 species each. This meant that families had unequal influence on the outcomes of the analyses. We performed a jackknife analysis to examine whether our results were affected by excluding individual families represented by four or more species in the full dataset. Exclusion of individual families did not significantly alter any of the phylANOVA results (Supplementary file 1 Tables S6-11). In general, exclusion of individual families did not affect the Brownie results (Supplementary file 1 Tables S12-17) or the BayesTraits results (Figure 6—source datas 23, Figure 7—source data 1). However, there were several notable exceptions, detailed below.

In several cases, removing a single family altered the significance of our findings. Removal of the Muscicapidae (three species) from the Brownie analysis of inter-song interval created a two-rate model that did not fit the data significantly better than the one-rate model (Supplementary file 1 Table S13). For song duration, removal of Fringillidae (three species) led to a two-rate model that was not a significantly better fit than the one-rate model (Supplementary file 1 Table S14). Removal of Mimidae (four species) from the Brownie analysis of syllables per song drastically changed the results, such that the two-rate model was no longer a significantly better fit than the one-rate model (Supplementary file 1 Table S15). All of the included Mimidae species are song-plastic, so these results suggest that mimids may be driving the the fit of the two-rate model for syllables per song. In addition, removal of Mimidae from the BayesTraits analysis of song plasticity and syllables per song altered the observed trends. At low threshold values of syllables per song, song plasticity with low syllables per song is an unstable state, with high rates of transition toward either higher syllables per song or toward song stability (77% of runs significant in this range); however, when mimid species are removed, gaining more syllables per song in the song-plastic state becomes far more likely. At moderate values of syllables per song, there are elevated rates of transition both toward more syllables per song and towards fewer syllables per song in song-plastic species (100% of runs significant in this range); there is a comparable trend when mimid species are removed. At high threshold values, more syllables per song with song stability is an unstable state (63% of runs significant in this range, Figure 6—figure supplement 2); however, with mimids removed, more syllables per song with song plasticity becomes the most unstable state (Figure 6—source data 3).

Due to the strong effect that the inclusion of Mimidae had on the Brownie analysis of syllables per song, we performed a second jackknife analysis, in which each of the four mimid species was removed in turn. Exclusion of Toxostoma rufum or Dumetella carolinensis had little effect on the results (Supplementary file 1 Table S18), leading to significant support for the two-rate model. However, exclusion of Mimus polyglottos (p=0.074) or Mimus gilvus (p=0.509) led to a two-rate model that did not fit the data significantly better than the one-rate model (Supplementary file 1 Table S18). Therefore, we concluded that these two Mimus species drove the estimated evolutionary rate of syllables per song in song-plastic species to be much greater than in song-stable species, and that faster evolution of syllables per song may not necessarily be a universal trend for song-plastic species.

Members of Mimidae are renowned for their vocal mimicry, frequently exhibiting improvisation and invention of syllables beyond simple imitation, and thus they lack the generally stereotyped song structure shown in other oscine families. Furthermore, mimids often have periods of continuous singing with minimal repetition of elements and irregular syllable spacing. Thus, quantification of song duration or number of syllables per song for mimids could be highly susceptible to listener perception (Wildenthal, 1965). Therefore, although we acknowledge that mimids are an important case study in extended learning durations, our results for the evolutionary rate of syllables per song might be more meaningful across all bird species when mimids are excluded, in which case we find the rate of evolution of syllables per song is be independent of adult song plasticity.

Discussion

Previously, it was unknown whether the song-learning window evolved in concert with song features associated with sexual selection, as predicted by a computational model of song learning (Creanza et al., 2016). This is a critical missing piece of the puzzle of song learning evolution, as previous evidence has suggested that sexual selection only acts upon the features of song and not the length of the song-learning window or maintenance of the song-learning pathways in the brain. Here, we performed phylogenetically controlled analyses to assess the interactions between the length of the song-learning window—using adult song stability versus plasticity as a proxy—and the evolution of song characteristics. Interestingly, we noted that several evolutionary events relatively early in passerine evolution accounted for much of the diversity in the song-plasticity period in our sample of species. We show that a bird’s ability to modify its song as an adult affects the characteristics of its species' song: adult song plasticity corresponds to larger syllable and song repertoires. Further, our results suggest that sexual selection for large repertoires could indirectly favor individuals with longer learning windows, driving the evolution of increased song plasticity.

We found two key trends in the trait correlation (phylANOVA) and evolutionary rate (Brownie) analyses. First, song plasticity affected the direction of evolution in traits that can be considered metrics of song complexity (syllable and song repertoire size, Figures 1 and 3, Table 1), leading to larger repertoires in species with adult song plasticity. Further, species with early song stability evolved their repertoires at a slower rate than species with longer learning (delayed song stability and song plasticity, Table 5), but song-plastic species did not evolve their repertoires at a faster rate than species with early or delayed song stability combined (Tables 69). Thus, while repertoires only evolve directionally in song-plastic species (Tables 23), our results suggest that extended learning through the first breeding season allows for faster, but nondirectional, evolution of repertoire size. A possible explanation is that delayed song learning allows individuals to modify their songs after migration and thus adapt their song to their new surroundings once they establish a territory, without necessarily corresponding with sexual selection for larger repertoires. This ability for an individual to adapt to a new local song might be beneficial, particularly when species have local dialect structure; however, this would not lead to directional evolution for any particular song feature, consistent with our findings. Second, song plasticity increased the rate of evolution primarily in traits that can be considered metrics of song performance (song duration, intersong interval, and song rate Figure 5, Table 2). While these performance-related song traits evolved faster in song-plastic lineages, this rapid evolution did not lead to significant differences in those song characteristics compared to song-stable lineages. A possible explanation may be that increases in repertoire size necessitate changes to song structure, but multiple structural aspects of song can be altered to accommodate these changes. Thus, there is no overall pattern of directional evolution in these other song characteristics. Alternatively, bird species may be required to adapt to increasing repertoire sizes while maintaining species-specific constraints imposed by innate aspects of song structure or female preferences for different performance characteristics. In the latter case, if information about innate characteristics and female preferences are known, it may be possible to predict how song traits will change in response to increasing repertoire sizes and greater adult song plasticity.

With our analyses of correlated evolution, we aimed to detect whether the state of the repertoire size or adult song stability versus plasticity consistently changes first in evolutionary history, facilitating a change in the other trait. Overall, our results suggest that there is not a consistent order of evolutionary transitions (Figures 67). For example, song stability with very large syllable or song repertoires, and song plasticity with very small syllable or song repertoires both formed evolutionarily unstable states, with high evolutionary rates of transition in both repertoire size and song-learning window. However, we do note that the fastest rates of transition in our analyses were those switching between song stability and song plasticity to leave those unstable states. This trend suggests that the magnitude of a species’ repertoire may be more likely to drive the evolution of learning window than vice versa (Figures 67). This is consistent with the idea that selection acting upon song features could indirectly place selective pressures on the learning window. We propose several hypotheses that could explain these evolutionary dynamics: 1) it may be disproportionately costly to maintain song plasticity when syllable or song repertoire sizes are very small, perhaps because the benefit of extra time to learn does not outweigh the metabolic cost of maintaining plasticity, or 2) species with small syllable or song repertoires may have highly stereotyped songs which are selected for based on accuracy of learning and consistency of song production, favoring males that only learn from their fathers or early-life neighbors. Alternatively, in species where females prefer larger repertoires, 3) it may simply require more time to learn a large song or syllable repertoire than is available with a short learning window, or 4) learning large repertoires may require too much energy devoted to song practice during the crucial period of development before and during fledging, favoring birds that can learn for longer periods. Further research into the physiological or reproductive costs of song plasticity is needed.

Beginning with Darwin (1872), numerous researchers have proposed that polygynous mating systems could lead to amplified sexual selection (Emlen and Oring, 1977; Payne, 1984). The elevated rates of transition that we observed between social mating systems in song-plastic lineages suggest an interesting hypothesis for further investigation: perhaps having a plastic song-learning program facilitates evolutionary transitions in mating systems.

Our results provide key evidence that sexual selection upon song characteristics might indirectly act upon the song-learning window. We do not fully understand the mechanisms underlying the maintenance or reopening of the song learning window in adulthood, but genetic, environmental, hormonal, and social factors are likely contributors (Eales, 1987; Eales, 1985; Kroodsma and Pickert, 1980; Morrison and Nottebohm, 1993; Nottebohm, 1969). For example, when zebra finches were reared in isolation, their sensitive periods were lengthened. These isolated birds maintained both gene expression profiles associated with song learning in the song system and high levels of neuronal addition to the HVC (a key region in the song system of the songbird brain) for longer than birds reared with an adult male tutor (Kelly et al., 2018; Wilbrecht et al., 2006), linking the neural underpinnings of song learning to the length of the song-learning window (Brenowitz and Beecher, 2005; Nordeen and Nordeen, 1990; Nordeen, 1997; Nottebohm, 1992). Furthermore, a positive association between HVC volume and song repertoire size has been demonstrated both intraspecifically (in song sparrows; Pfaff et al., 2007) and interspecifically (Devoogd et al., 1993). In light of our findings that adult song plasticity correlates with an increase in song repertoire size, there is a logical prediction that extended song learning may be associated with increased HVC volume across species. This is an important avenue for future research.

Although our dataset includes species from 24 different songbird families, many families are not represented due to a lack of data about song stability. It will be important to expand this dataset in future studies. It would also be interesting to explore the evolutionary interactions between adult song plasticity and mimicry of heterospecific sounds, which has been observed in Mimidae and numerous other clades across the songbird lineage (Goller and Shizuka, 2018). With our current dataset, we could not adequately explore the effects of mimicry on the evolution of song learning outside of the mimids, but the repeated evolution of mimicry makes it a particularly interesting topic for follow-up studies on the length of the song-learning window. In addition, different song metrics that are tailored to mimicry would be important in studying the evolution of vocal mimics and the dynamics of their unique vocal patterns. Furthermore, there is increasing interest in the importance of female song in species, which is more common than previously thought (Langmore, 1998; Odom et al., 2014). Our dataset includes numerous species wherein females are known to sing at least occasionally (Langmore, 1998), but the length of the song-learning window in females has not been assessed in any of these species. There is, however, some evidence that female birds can modify their song preferences in adulthood (Nagle and Kreutzer, 1997). Thus, it remains an open question whether song plasticity affects the evolution of female song in the same way it affects male song, and whether species with adult song plasticity in males also have adult song plasticity in females.

Our findings shed new light on the broader subject of song evolution, specifically the evolution of the process of song learning. We hypothesized that sexual selection on certain aspects of song could in turn alter the selection pressures on the length of the song-learning window. Here, we performed phylogenetically controlled analyses across 67 songbird species to assess the evolutionary interactions between song characteristics and song plasticity in adulthood. With these analyses, we show the first evidence for this evolutionary relationship. Adult song stability versus plasticity may be evolutionarily dependent upon the properties of the song itself: large syllable and song repertoires appear to drive the evolution of adult song plasticity and thus open-ended song learning. This provides context for the remarkable interspecific variation in song-learning windows across the songbird lineage and suggests an evolutionary mechanism by which sexual selection might have influenced the evolution of songbird brains.

Materials and methods

Data collection

Request a detailed protocol

It is difficult to precisely measure the length of the song-learning window in both field and lab studies. In field studies, if a bird is recorded singing a new song element in its second year, researchers often cannot rule out the possibility that it learned that element during its first year but did not incorporate it into its song until later (Marler and Peters, 1981; Vargas-Castro et al., 2015). Likewise, if a bird learns a new element in its second year but elects not to produce it, a recordist would be unlikely to capture it. In contrast, raising birds in the laboratory with known song exposure enables researchers to assess when a bird learned a particular song element (Chaiken et al., 1994; Marler and Peters, 1988; Nelson (1998); Nottebohm, 1969) but also raises questions about whether the lab-reared birds are exhibiting their normal behaviors (Baptista and Petrinovich, 1984; Kroodsma and Pickert, 1984). In addition, these lab-rearing procedures have only been performed in a handful of species. In this paper, we examined studies that include longitudinal measures of adult song stability versus plasticity (e.g. Nordby et al., 2002) to determine whether song is modified across time as a proxy for the length of the song learning window.

We performed a literature search to gather information on the stability or plasticity of male song in oscine species. Studies with information about learning style were found via Google Scholar using the following search terms: [species name] or [common name] in combination with ‘open-ended,’ ‘close-ended,’ ‘closed-ended,’ ‘age-limited,’ ‘crystal*,” ‘adult learning,’ and ‘song changes.’ We used three strategies to assign the song stability classification for a species. We first defined a species as having adult song stability (‘song-stable species’) if the literature indicated that males did not modify their songs after the first breeding season. Species in which males modified their syllable repertoires after their first breeding season were classified as having adult song plasticity (‘song-plastic species’). This strategy was meant to approximate the dichotomy of open- and closed-ended learners often used in the birdsong field. We made two exceptions to this for the species Cacicus cela and Phoenicurus ochruros, in which males do not gain their mature plumage until their second breeding season and may therefore be delayed in reaching sexual maturity relative to other bird species (Draganoiu et al., 2014; Trainer and Parsons, 2002). Because these birds cease modifying their repertoires before reaching their second breeding season with mature plumage, they were considered song-stable (Source data 1). Additionally, past research defined Melospiza lincolnii as an ‘open-ended improviser.’ It is unclear whether improvisation throughout the lifespan is equivalent to learning throughout the lifespan, however it does fit our definition for adult song plasticity. For the main analysis, we considered this species to have adult song plasticity, but we repeated the main analysis with this species reclassified as having adult song stability. This reclassification produced a negligible effect on our results (Supplementary file 1 Tables S19-20). For our second strategy, we separated song-stable species into two sub-groups: early song-stable (species that cease modifying their songs before the first breeding season) and delayed song-stable (species that modify their songs during their first breeding season but not after). There was not enough information to make this determination for some species, so our dataset was reduced to 59 species. We used these same 59 species to generate a continuous measure of song stability for our third strategy. Because no information was available about the prevalence of song changes beyond the second breeding season for most of the species in our dataset, this measure only ranges from 0 to 2 years, and all song-plastic species were assigned a value of 2. Furthermore, it was not clear exactly when most of the delayed song-stable species stopped learning, so they were given a value of 1.33, at which point the first breeding season should have ended. The two species mentioned above that display delayed song and plumage maturation were assigned a value of 2.

To gather data on the song characteristics, we performed a literature search via Google Scholar and Web of Science using the search terms ‘Passeriformes’ or [species name] in combination with ‘syllable repertoire’ or ‘song repertoire.’ This yielded a mix of primary sources and studies that had previously aggregated repertoire size data. We also gathered data from the curated field guides Birds of North America (Rodewald, 2015) and Handbook of Birds of the World (del Hoyo et al., 1992). We did not perform explicit searches for any of the other included song characteristics, but we collected this data whenever we encountered it. Song characteristic nomenclature is variable across studies, so, when possible, we read the methods of the primary sources to ensure that the data collected were consistent with our definitions for song characteristics.

We utilized the following definitions for song characteristics:

  1. Syllable repertoire - the mean number of unique syllables sung per individual.

  2. Syllables per song - the mean number of unique syllables per song.

  3. Song repertoire - the mean number of unique songs produced per individual.

  4. Song duration - the mean length of a song, delineated by pauses or beginnings of a repeated motif (seconds).

  5. Intersong interval - the mean pause length between songs (seconds).

  6. Song rate - the mean number of songs produced per minute (calculated value; Rate = 60/(Duration + Interval)).

  7. Continuity - the fraction of singing time spent actively producing song (calculated value; Continuity = Duration/(Duration + Interval)).

For three species with song-learning window data, syllable repertoire size estimates were not available in the literature, so we estimated these repertoire sizes from published sources or song recordings. For Philesturnus rufusater (Supplementary file 1 Table S21), (Jenkins, 1978) coded songs into different types and gave the song repertoire types of each male studied in Table 2 of that paper. Information from one male was missing from this table. We deduced the repertoire of the missing male by first looking at Table 1 from that paper (Jenkins’s Table 1), which showed that 16 males had a repertoire size of one song. Only 15 of the males in Jenkins’s Table 2 had a repertoire size of one song. Thus, the missing male had a repertoire size of one song. We then compared the bands of males present in Jenkins’s Table 2 to the territory map in Figure 7 of that paper, and A_RW was the only male missing from Jenkins’s Table 2. A_RW was located in the DC region of Jenkins’s Figure 7, so we assigned that as his repertoire. Jenkins notes that neighboring males share song types, so the only other song A_RW could have known instead of DC was ZZ, which has the same number of unique syllables as DC. For Geospiza fortis (Supplementary file 1 Table S22) and G. scandens (Supplementary file 1 Table S23), we used recordings from the Macaulay Library at the Cornell Lab of Ornithology (Cornell Lab of Ornithology, 2019) to estimate the syllable repertoire size and syllables per song. For G. fortis, we also used the sonogram examples present in Grant and Grant (1996).

When the song repertoire for a species equalled one, we assumed that its syllable repertoire was equal to its number of syllables per song. In many cases, there were multiple studies that gave different estimates for a given song variable in one species. To handle these discrepancies, we created three datasets (Source data 2). For one, the main dataset, we used the median value across studies. For the second dataset, we used the minimum value reported in the literature for each species, and for the third, we used the maximum value reported in the literature. We log normalized all song trait data. The three datasets revealed similar results; analyses using the maximum and minimum literature values are presented in Supplementary file 1 Tables S1-S5.

We also cataloged data on mating behavior at the species level. In particular, we assembled binary classifications of social mating system (monogamy vs polygyny) and extra-pair paternity (low EPP vs high EPP). We considered a species to be monogamous or polygynous when a source unambiguously categorized that species’ social mating system; that is, we did not assign a social mating system to species labeled ‘probably,’ ‘usually,’ ‘mostly,’ ‘normally,’ ‘typically,’ and ‘generally monogamous/polygynous,’ etc. unless quantitative measurements were also provided. When quantitative data were available, species were defined as polygynous when at least 5% of males had more than one social mate, as in Snyder and Creanza (2019). A review of extra-pair paternity studies estimated an average of ~11% of offspring per nest were attributable to extra-pair mates across species (Griffith et al., 2002). In accordance with this estimate and with previous studies that used a binary classification of EPP (Snyder and Creanza, 2019; Soma and Garamszegi, 2011), we used a 10% threshold for either extra-pair young or nests containing at least one extra-pair chick to estimate the frequency of extra-pair paternity in that species (<10% = low EPP; ≥10% = high EPP, Source data 2).

Assessing the evolutionary history of adult song plasticity

Request a detailed protocol

To predict the rate of transition between adult song stability and adult song plasticity, we used the ace function from the R package Phytools and a phylogeny from Jetz et al. (2012). We note that this phylogeny has broad coverage of oscine songbirds but is based on limited genetic data (often mitochondrial), so the topology could change as more avian genomes are sequenced (Jarvis et al., 2014; Lamichhaney et al., 2015). With this tree, we tested whether an all rates different (ARD) model fit the data significantly better than the equal rates (ER) model using an ANOVA. We then used the better-fit equal rates model to generate 10,000 trees with make.simmap (Phytools). This function uses the rate from ace and a phylogenetic tree with annotated tips to create stochastic simulation maps for the potential evolutionary transitions between the song-stable and song-plastic states. We found the predicted ancestral state for each of these 10,000 simulations and used countSimmap (Phytools) to count the total number of transitions that occurred in each map. The minimum number of predicted evolutionary transitions across these 10,000 simulations was considered to be the most parsimonious; we also compared this to a manual count of evolutionary transitions starting from either ancestral state.

Detecting differences in song characteristic evolution in song-stable vs. song-plastic species

Request a detailed protocol

To test whether there were significant differences between song-stable species and song-plastic species for the song traits, we performed a phylogenetically controlled ANOVA (phylANOVA, Phytools) for each song characteristic. We repeated this analysis with the subset of species we classified into early song-stable, delayed song-stable, and song-plastic. Because there were relatively few early song-stable species in this dataset, we only performed this re-analysis for song traits that had data for at least nine early song-stable species (syllable repertoire size: nine species with early song stability, song repertoire size: nine species, and syllables per song: ten species). In this paper, we visualize the predicted ancestral traits on the phylogenetic tree with color and pie graphs, however, the raw values are available in Figure 1—source data 1. To test for correlations between song characteristics and the continuous values for the duration of song plasticity, we performed a phylogenetic generalized least squares (PGLS) analysis. We used the function gls (R package: nlme), with the ‘correlation’ parameter lambda computed using the function corPagel (R package: ape).

To test whether adult song stability or plasticity affected the rate of evolution for the song characteristics, we used the function brownie.lite (R package: phytools). This function first calculates a one-rate model of evolution for a song characteristic using a phylogenetic tree and the current states of the tips for that song trait. This one-rate model assumes that change in the value of the song characteristic is random across evolutionary time and can be approximated by Brownian motion. Next, a model is generated wherein two different rates are calculated; this two-rate model assumes that the evolution of the song characteristic has one rate in the song-stable state and a different rate in the song-plastic state. This model requires estimations for the ancestral states of song stability for each branch of the phylogeny. To create these estimates, we used the function ace to calculate the rate of transition between the song-stable and song-plastic states for the full dataset. We then used these transition rates to generate 1300 different stochastic simulation maps (make.simmap) for the subset of species that had data for each song trait. For the Brownie analysis, we tested whether the two-rate model fit the data significantly better than the one-rate model by performing a chi-square test on the mean log likelihoods of the two models. We repeated this analysis for the set of species we classified as early song-stable, delayed song-stable, and song-plastic for traits for which we had data on at least nine early song-stable species. We compared the three-rate model to the one-rate model. We also reran the two-rate model in this reduced dataset by combining the early and delayed song-stable groups and testing whether the three-rate model was better than the two-rate model. Because the delayed song-stable trace peaked at a similar position to the song-plastic trace for syllable and song repertoire size (Figure 4D–E), we also compared the three-rate model to another version of the two-rate model, in which one group was early song-stable (shorter learning), and the other was delayed song-stable and song-plastic combined (longer learning).

We used BayesTraits to test for correlated evolution between song stability and song characteristics, or, in other words, whether the rate and direction of evolutionary transitions of one trait are dependent on the state of another trait, and whether an order of transition events can be inferred. Specifically, we tested the hypothesis that an evolutionary change in song stability increases the likelihood of an evolutionary change in certain song variables or mating behaviors, or vice versa. BayesTraits compares two models of discrete trait evolution for a pair of binary traits and a given phylogenetic tree: 1) an independent model (i.e. the evolution of one trait does not depend on the other trait) and 2) a dependent model (i.e. the evolutionary transitions of each trait depend on the state of the other trait, suggesting correlated evolution). Using the maximum likelihood method, BayesTraits reports marginal likelihoods for the computed dependent and independent models (function Discrete in package btw [BayesTraitsWrapper]), allowing us to determine whether the dependent model describes the data significantly better than the independent model. We used function LRtest (package: lmtest) to perform the likelihood ratio test. Since this model required both input traits to be binary, we classified the continuous song characteristics as binary groups (‘low’ or ‘high’) based on a delineating threshold. Instead of choosing the threshold arbitrarily, we used each unique value of the song characteristic data as the threshold and repeated the analysis 100 times at each threshold. This method of using a spectrum of thresholds to delineate the ‘low’ and ‘high’ value categories resulted in transition rates that varied dramatically depending on where the threshold was placed. In essence, when the threshold was set as a value in the bottom third of the unique trait values present in the data, the analysis evaluated the rate of transition from low-to-moderate and larger values for a song trait and vice versa. When the threshold was set as a value in the upper range of the unique trait values present in the data, the analysis calculated the rates of transition from higher song trait values to medium-to-low values. To account for this nuance, we binned the threshold data into two to five bins, with three bins as the default: low (bottom third of unique trait values), medium (middle third) and high (top third). We then calculated the mean of each state transition rate in each bin. In addition to the song characteristics, we also analyzed song stability versus social mating system (i.e. social monogamy or polygyny) and rate of EPP. These analyses were performed for 1000 runs each.

Jackknife analysis

Request a detailed protocol

Some families of birds were well represented in our sample, while others were only represented by one or two species. To test whether any well-represented family significantly skewed our results, we removed each family that was represented by four or more species in the full dataset in turn, and repeated the phylANOVA, brownie.lite, and BayesTraits analyses. Jackknife analyses were only performed when significant results were obtained in the main analysis. Thus, all song traits except continuity were tested in the phylANOVA and brownie.lite jackknife analysis, while only syllable repertoire, song repertoire, and syllables per song were tested in the BayesTraits jackknife analysis. Each Brownie analysis was run on 1300 unique stochastic character maps, and each BayesTraits analysis was repeated 20 times. We determined the family of each species based on its classification in the 2017 version of the eBird Clements Integrated Checklist (Clements, 2007). The family Locustella was combined with Acrocephalidae, as Acrocephalidae was paraphyletic when Locustella was considered to be a separate family. The Mimidae family alone had a large effect on the syllables per song metric, so we performed another jackknife analysis with phylANOVA and brownie.lite by removing each mimid species in turn.

Correction for multiple testing

Request a detailed protocol

We used a Holm-Bonferroni correction to control for testing multiple hypotheses with the same data (Holm, 1979). This correction is appropriate for data wherein the outcome of one test is likely to be related to the outcome of another test, as would be the case for song characteristics.

Data availability

All data are made available as supplementary information provided with this manuscript, and are also provided at https://github.com/CreanzaLab/SongLearningEvolution (copy archived at https://github.com/elifesciences-publications/SongLearningEvolution).

The following data sets were generated

References

  1. Book
    1. Catchpole CK
    2. Slater PJB
    (2003)
    Bird Song: Biological Themes and Variations
    Leiden: Cambridge University Press.
  2. Book
    1. Clements JF
    (2007)
    The Clements Checklist of Birds of the World
    Ithaca: Comstock Publishing Associates.
  3. Book
    1. Darwin C
    (1872)
    The Descent of Man, and Selection in Relation to Sex
    London: J. Murray.
  4. Book
    1. del Hoyo J
    2. Elliott A
    3. Sargatal J
    4. Christie DA
    5. de Juana E
    (1992)
    Handbook of the Birds of the World Alive
    Barcelona: Lynx Edicions.
    1. Hausberger M
    2. Jenkins PF
    3. Keene J
    (1991)
    Species-Specificity and mimicry in bird song: are they paradoxes? A reevaluation of song mimicry in the european starling
    Behaviour 117:53–81.
    1. Holm S
    (1979)
    A simple sequentially rejective multiple test procedure
    Scandinavian Journal of Statistics, Theory and Applications 6:65–70.
    1. Immelman K
    (1969)
    Song development in the zebra finch and other estrildid finches
    Bird Vocalizations pp. 61–74.
    1. Marler P
    (1990) Song learning: the interface between behaviour and neuroethology
    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 329:109–114.
    https://doi.org/10.1098/rstb.1990.0155
    1. O'Meara BC
    2. Ané C
    3. Sanderson MJ
    4. Wainwright PC
    (2006)
    Testing for different rates of continuous trait evolution using likelihood
    Evolution 60:922–933.
  5. Book
    1. Rauschecker JP
    2. Marler P
    (1987)
    Imprinting and Cortical Plasticity: Comparative Aspects of Sensitive Periods
    New York: Wiley-Interscience.
  6. Book
    1. Rodewald P
    (2015)
    The Birds of North America
    Ithaca: Cornell Laboratory of Ornithology.

Decision letter

  1. Patricia J Wittkopp
    Senior and Reviewing Editor; University of Michigan, United States
  2. Ofer Tchernichovski
    Reviewer; Hunter College, United States

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

Thank you for submitting your article "Correlated evolution of repertoire size and song plasticity predicts that sexual selection promotes open-ended learning" for consideration by eLife. Your article has been reviewed by three peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Patricia Wittkopp as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Ofer Tchernichovski (Reviewer #3).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

The authors provide some of the first evidence that an important and variable (but previously poorly understood) aspect of bird song learning – the song learning window – may be shaped by selection for large syllable and song repertoires. Thus, selection pressures on the song structure itself may be selecting for this variation, whereas, in the past researchers have often thought that environmental variation and breeding season length were leading factors shaping this trait. This is an exceptionally novel study and finding that is likely to spur much future research in this area and change how we think about how this aspect of song learning behavior evolves. Two of the reviewers had more minor comments, but a third had more serious concerns about presentation and data availability. I have summarized their comments (often sharing them verbatim) below.

Essential major revisions:

Please pay special attention to the comments below on:

1) The appropriateness of ancestral state reconstruction.

2) The methods on how birds that learn only into the beginning of their second breeding season were handled.

3) The fact that 10% EPP seems low for considering those species as having EPP.

4) Including example sonograms and links to audios.

5) Being more explicit that the behaviors are a long a continuum.

In addition, any additional evidence that can be added on the attractor models would also be appreciated. We recognized this might be challenging, however.

Additional specific comments: (from individual reviewers, lightly edited, in a combined list below)

1) The most critical reviewer writes: It is obvious that song learning takes time: the more syllable types and transitions a bird of a certain species 'happens' to learn, the longer it should take. The authors emphasize this as their principal result, but there is much more in this work than this generic conclusion, which by itself, is not very surprising. It is quite obvious to me that repertoire size should be correlated with open ended learning (adult song plasticity). On the other hand, the results about song evolution rate across categories, and about attractor and repeller states in evolutionary dynamics (Figures 4-5) are much more interesting. The problem is that, without some additional analysis, the case is not strong enough and not even clear enough. Throughout the paper no specific example are presented in a manner that could convince the reader that the evolutionary trajectories are interesting and robust – I think they are, but authors chose not to present any examples that could help me judge.

2) It is surprising that a study about the evolution of song complexity and performance does not show a single sonogram, and does not present even a single example of an evolutionary trajectory. Abstract statistical models with measures that are combined across many species can be used for presenting scientific results, but only after specific 'typical' examples are presented. Specifically, the phylogenetic tree presented in Figure 2 contains several trajectories that look interesting: going either from stability to plasticity, or the other way around over several transitions. Those cases should be analyzed separately to examine if adult plasticity evolved before, during or after changes in features such as syllables per song, etc. Those examples should be presented as sonograms (with attached media) so that readers can assess those changes. One way to do this would be to present at least two examples of each type in the main paper, and then the supplement can show several more. Those specific examples should lead the overall abstract analysis.

3) Evolution of adult song plasticity: I am a bit dissatisfied with the dichotomy between open ended and close ended learners in the input set, which could have affected the number of transitions required by model. Say that we had a continuous variable, such as the lifetime duration of song learning. It seems likely that an ancestor with a narrow song learning window (such as in a zebra finch) could gradually evolve an extended period, which only in specific, extreme cases, would show a phenotype of open-ended learning and adult plasticity. So, I would like the authors to look more carefully at the evolutionary trajectories, and within a small number of species assess the possibility of continuity. I am not sure if sufficient data to make such an assessment are available, but I would assume that song development should correlate with sexual maturity and growth rate, which may be easier to compare across species.

4) I would have thought that the rate of song feature changes should be slower when repertoire size is large (as it takes longer to turn a cargo ship than to turn a canoe). So maybe the rate of change analysis should take this issue into account?

5) The analysis of attractor and repeller states is potentially very interesting. But many arbitrary decisions about thresholds and input measures make me worried about this analysis. First of all, there is a need to look at all those results continuously – if there is an attractor state in low threshold and a repeller state in high threshed, it could mean that there is an intermediate state of equilibrium. In other words, the evolutionary trajectories might be 'hovering' around some stable 'prior'. More analysis is needed here across graphs. Here too, I would start by looking at specific trajectories (and please present sonograms).

6) Throughout the paper I kept wondering if there is a simple way of telling which came first: song learning duration or syllable per song? Can you look at specific trajectories to get some clues, even if not statistical?

7) Also, since mimids drive much of the correlations maybe it is a good idea to look at them more carefully as a post hoc analysis. I understand that authors want to rid of unrepresentative example, but this does not mean that the most convincing cases should be ignored…

8) Overall, it would be good to show, at least in the supplement, some examples of evolutionary of complexity and performance side by side with evolution of open-ended learning. I feel that without seeing what's going on in specific cases it is difficult to judge the quality of the data.

9) Regarding data quality of adult song plasticity: authors used a categorical measure (taken from multiple published papers) and there is no mention in the Materials and methods about assessment and quality control of integrating such a categorical assessment across studies. Song traits, on the other hand, are continuous variables. I think that this categorical, binary assessment of song plasticity has limitation that should be discussed.

10) Line numbers would make reviewing easier.

11) It is not clear whether the data included were all from males in each species, or whether some females were included. If a mixture of sexes, I think this factor should controlled for in one or more analyses. The authors mentioned females in the discussion for future studies, but never mentioned (maybe assumed) that they only analyzed studies with males in this work.

12) The authors should make a cautionary statement that simulations may not match the real biology, but approximate it.

13) The Introduction seems excessively long on background material and different hypotheses. I suggest dealing with the various hypothesis in the Results, and only mention two in the Introduction, so that the paper is less repetitive.

14) Also in the Introduction and other parts of the paper, stability and plasticity are treated as different variables (whether in written text or some of the figures). But they are along a continuum, that the authors results support. I suggest making this clearer.

15) On its first mention, the authors need to cite what phylogenetic tree of oscine families that they used. The accuracy of the tree could affect their results. I couldn't find any citation in the main text for the tree used.

16) About the tree results, it looks pretty clear from all the phylogenetic analyses that there was an early split or two splits in the songbird lineage between two predominantly plastic song clades and a predominantly stable song clade. And then some independent deviations from those early splits moving into the other clade's phenotype. Maybe the common songbird ancestor was in the middle of the continuum. This exciting finding of the major clades with plastic vs. stable song should be highlighted in the Abstract and Results.

17) Subsection “Do our results depend on the particular values or families included in analyses?”, first paragraph. I thought other song production factors, like song intervals, were not affected in the above analyses? But they are here, when using maximum or minimum values?

18) The Discussion is on the long side. It could be shortened without losing content, by not repeating the results, but simply moving into a discussion of them.

19) The following two sentences appear to contradict each other: "We did not find evidence for correlated evolution between song stability and EPP. Most of our simulations testing for correlated evolution between mating system and learning window were significant;"

20) Discussion, seventh paragraph. Explain what is HVC, to scientists unfamiliar with the avian brain.

21) Figure 1: Label y axes (log scale). Mention stat results in the legend

22) Figure 2. Legend text is hard to follow, in terms of grammar.

23) Please consider archiving the data in a proper repository so that they can be more readily built upon.

24) Abstract: Song-plastic species. This is not a term I am used to. Do you mean open ended learners? If so, can you use that term? Or at least consider "species with plastic song" – and perhaps briefly define what you mean.

25) Abstract: "Notably…": What is the evidence for this? Is this derived from the above result, or is it from other results. Please preface this with what kind of analysis demonstrated this result.

26) Introduction – "small-scale analysis of comparative evolution": more appropriate = "small-scale comparative analysis".

27) Introduction: "female choice can favor certain characteristics" – can you give a bit more description of process, or what is selected for (similar to the description for neuroplasticity). This is crucial to your findings.

28) Introduction, third paragraph – I would expect this before talking about song-learning windows – as it is the more standard dogma under which bird song is studied.

29) Results – first paragraph: this is all good info, but does not present any results. It seems better suited for the Materials and methods (and, in fact, this useful explanation of why you went with the novel, broad "song-plastic/stable" terminology is left out of the Materials and methods. I would recommend moving more of this text there and skip straight to your classification and Results (with a brief definition here)).

30) Results – second paragraph: the idea of examining the ancestral state of all songbirds for song plasticity is interesting, however, you don't really have an appropriate dataset for this. Because you have so many missing taxa (only 64 species out of ~5000 categorized) and they fall throughout the songbirds (you could predict the ancestral state for a particular family if you had sampled many representatives of that family, or all songbirds if you densely sampled representatives of lineages branching more directly from the common ancestor of all songbirds). However, your dataset is a good one for examining correlated evolution, as you have done, and you can also get an idea of transition rates and when/ how often transitions between song stability and plasticity occur (although, it is good to be aware that this may be a fairly rough estimate just because there are probably many species that are missing among the species sampled). I would consider removing/avoiding mention of ancestral states, and approach this in light of major transitions among these traits. You could mention that the ancestral state for these lineages could not be resolved.

31) Results – paragraph beginning "Unexpectedly" – I'm not sure this is that unexpected (I'm not certain what I would predict here). Perhaps skip this word.

32) Results – subheading "Does song stability or plasticity affect the direction of evolution for song traits and mating system?": seems to me the predicted direction of causality is backwards here. I would predict that song learning strategy is dictated by ecological traits or the characteristics of the songs that are being selected (as you find). Therefore, I suggest alternate subheading and approaching this discussion from the perspective of "Is song stability or plasticity influenced by the evolution of song traits or mating system?"

I think this is also approached in the same way in the Abstract/Introduction. There I also suggest pitching this the opposite way (that's how I automatically think of it).

33) Discussion – first paragraph: "in concert with sexually selected song features" – we don't know that these features are strictly sexually selected. Consider "song features often thought of as sexually selected" or "song features associated with sexual selection".

34) Discussion – third paragraph: the two sentences after "We observed that with song stability….". Can this be explained more simply? A lot of jargon here (repeller/attractor states). Can you explain these trends in terms of what is evolutionarily stable/unstable – e.g. song stability coupled with small syllables and song repertoires is an evolutionarily unstable state, so transitions away from it are high.

35) Discussion – fourth paragraph: I think you undersell the importance of mimics (not strictly mimids) in songbird song evolution here, and the implications of your result including the mimids. Mimicry is common throughout songbirds – including in the earliest-branching songbird lineage (the lyrebirds). Really the sample size of the study is small compared to all songbirds (understandably limited by our knowledge of the song learning window for most species). However, I think it would be nice to acknowledge here or elsewhere that this is a relatively small sample, but that mimicry is common in songbirds, so the pattern seen in the mimids may be meaningful and it will be interesting to see how it holds up as we learn more.

36) Discussion – overall comment – there are a lot of thoughtful points here, but the Discussion is rather long. I wonder if these can be delivered more briefly. Perhaps some points can be combined, or overall length cut?

37) Materials and methods – first paragraph: "adult song plasticity ('song-stable species')" – I think there is a typo here – doesn't this statement refer to species without song plasticity?

38) Materials and methods – definitions for 'Song-stable/plastic': there seems to be a gap in your definitions – which leaves birds that learn prior to their second breeding season but do not subsequently modify their song following the second breeding season. I am under the impression that this is not an uncommon strategy. How did you handle these species?

39) Materials and methods: defined high EPP when at least 10% of offspring had a genetic male parent other than social male. This seems low to me. What was the justification for this cutoff?

https://doi.org/10.7554/eLife.44454.050

Author response

Essential major revisions:

Please pay special attention to the comments below on:

1) The appropriateness of ancestral state reconstruction.

2) The methods on how birds that learn only into the beginning of their second breeding season were handled.

3) The fact that 10% EPP seems low for considering those species as having EPP.

4) Including example sonograms and links to audios.

5) Being more explicit that the behaviors are a long a continuum.

In addition, any additional evidence that can be added on the attractor models would also be appreciated. We recognized this might be challenging, however.

Additional specific comments: (from individual reviewers, lightly edited, in a combined list below)

1) The most critical reviewer writes: It is obvious that song learning takes time: the more syllable types and transitions a bird of a certain species 'happens' to learn, the longer it should take. The authors emphasize this as their principal result, but there is much more in this work than this generic conclusion, which by itself, is not very surprising. It is quite obvious to me that repertoire size should be correlated with open ended learning (adult song plasticity). On the other hand, the results about song evolution rate across categories, and about attractor and repeller states in evolutionary dynamics (Figures 4-5) are much more interesting. The problem is that, without some additional analysis, the case is not strong enough and not even clear enough. Throughout the paper no specific example are presented in a manner that could convince the reader that the evolutionary trajectories are interesting and robust – I think they are, but authors chose not to present any examples that could help me judge.

Thank you for this feedback! We agree that it is appealingly intuitive that larger repertoires would be correlated with adult song plasticity, but there is of course nuance in this relationship:

We note that (1) the distribution of repertoire sizes shows substantial overlap between species with adult song plasticity and adult song stability (Figure 1), (2) that there are species that do not conform to this trend (Acrocephalus palustris, Philesturnus rufusater, etc., Figure 2), and (3) that several species with adult song plasticity have been shown to modify which syllables they perform as they age rather than significantly increasing their repertoire sizes over time. Thus, it was important to establish the evolutionary correlation between larger repertoires and longer periods of song plasticity as the baseline for the rest of our study. Your point is well taken that, in an effort to conduct a broad evolutionary study, we did not highlight any specific examples or show individual songs. We now highlight in Figure 3 a genus that appears to show a relatively recent evolutionary trajectory that includes both a transition to closedended learning and a corresponding decrease in repertoire size. Since our repertoire estimates were primarily pulled from the literature (except in a few cases where we note that we estimated repertoire size ourselves from song data), we do not want to mistakenly give the impression that we routinely analyzed individual song recordings for our analysis, and we hope we balanced this appropriately.

2) It is surprising that a study about the evolution of song complexity and performance does not show a single sonogram, and does not present even a single example of an evolutionary trajectory. Abstract statistical models with measures that are combined across many species can be used for presenting scientific results, but only after specific 'typical' examples are presented. Specifically, the phylogenetic tree presented in Figure 2 contains several trajectories that look interesting: going either from stability to plasticity, or the other way around over several transitions. Those cases should be analyzed separately to examine if adult plasticity evolved before, during or after changes in features such as syllables per song, etc. Those examples should be presented as sonograms (with attached media) so that readers can assess those changes. One way to do this would be to present at least two examples of each type in the main paper, and then the supplement can show several more. Those specific examples should lead the overall abstract analysis.

We originally did not include any sonograms because we gathered our values regarding song complexity and performance from specieslevel estimates from previously published works. As mentioned above, we have now highlighted one example of song stability changing within a genus and included sonograms from xenocanto for visual inspection.

As far as analyzing some of these interesting cases separately to test for order of evolution, we now clarify in the text that this was the hypothesis we were testing with the Bayes Traits analysis. This test for correlated evolution can sometimes reveal a predominant order of transitional events, though it did not conclusively do so in our study. We did find that a transition in learning window was often more likely than a transition in repertoire size when comparing the predicted rate of evolutionary transitions out of unstable states, and we posit in the text that this change in learning window might often come first.

Since you specifically mentioned syllables per song, we also want to clarify those results: since all of the studied mimids were song plastic, and the change in the rate of evolution depended almost entirely on one species (Mimus gilvus), we would hypothesize that this change in syllables per song occurred after the transition to songplasticity, though more data on mimids would be needed to be sure.

3) Evolution of adult song plasticity: I am a bit dissatisfied with the dichotomy between open ended and close ended learners in the input set, which could have affected the number of transitions required by model. Say that we had a continuous variable, such as the lifetime duration of song learning. It seems likely that an ancestor with a narrow song learning window (such as in a zebra finch) could gradually evolve an extended period, which only in specific, extreme cases, would show a phenotype of open-ended learning and adult plasticity. So, I would like the authors to look more carefully at the evolutionary trajectories, and within a small number of species assess the possibility of continuity. I am not sure if sufficient data to make such an assessment are available, but I would assume that song development should correlate with sexual maturity and growth rate, which may be easier to compare across species.

Thank you for this helpful feedback. In response, we have implemented two new sets of analyses that we think have strengthened the manuscript.

(As a side note: for this revision, we were able to add a few more species to our analysis. We had the great honor of hosting Rosemary and Peter Grant for our studentinvited seminar series, and cofirstauthors Cristina and Kate described this project to them. They were excited to see our results, and they also helped us add two of Darwin’s finches to our dataset. We were also able to determine repertoire size data for another species, bringing our total to 67 species.)

To address the gap implied by our original framing, we modified our definition of songplastic to be birds that modify their repertoires after their first breeding season. We then revisited the literature on each of our species to gather more detailed information about the length of the songlearning window if available.

First, when possible, we reclassified the species into three categories: 1) those that stop changing their song before their first breeding season (early songstable), 2) those that modify their songs during their first breeding season but not after (delayed songstable), and 3) those that modify their songs after their first breeding season is over (songplastic). Species that showed delayed learning until their second year and also delayed plumage maturation were added to the delayed songstable group. This threestate classification was possible for 59 out of the 67 species.

Second, for the same 59 species, we used reported estimates of the ages at which song stabilizes in a species to create a continuous measure of the songlearning window. There were exceedingly few studies that examined song changes after the second breeding season, so our continuous metric ranges from 0 to 2 years.

Thus, we can now analyze the phenomenon of song stability/plasticity in three ways; in addition to the previous binary categorization (adult song stability vs. adult song plasticity) we have classified the data as a threestate discrete trait and as a continuous trait.

We present results using these two new classification schemes in the main text. As described in our response to the previous comment, we also highlight a specific example with spectrograms to illustrate interesting evolutionary transitions.

4) I would have thought that the rate of song feature changes should be slower when repertoire size is large (as it takes longer to turn a cargo ship than to turn a canoe). So maybe the rate of change analysis should take this issue into account?

We gave a lot of thought to this question. (A limitation, of course, is that archived recordings do not go back far enough to truly know the rate of song change on an evolutionary scale.) One implication of the canoe/cargo ship metaphor is that species with large and small repertoires should have a similar scale of variation in song traits: for example, if, on average, syllables repertoires vary by a few syllable per generation, then a small repertoire would indeed change faster than a large repertoire. We analyzed data across species to see whether this metaphor holds.

We plotted mean repertoire size versus standard deviation in repertoire size (data from Garamszegi et al., 2005).

From this plot, we notice that the standard deviation of syllable repertoire size increases with the mean syllable repertoire size. Thus, it appears that a cargo ship might be able to change at the same rate as a canoe, so to speak. We might hypothesize that other song features would show a similar pattern: the standard deviation of song duration, rate, etc. might also increase as syllable repertoire size increases.

5) The analysis of attractor and repeller states is potentially very interesting. But many arbitrary decisions about thresholds and input measures make me worried about this analysis. First of all, there is a need to look at all those results continuously – if there is an attractor state in low threshold and a repeller state in high threshed, it could mean that there is an intermediate state of equilibrium. In other words, the evolutionary trajectories might be 'hovering' around some stable 'prior'. More analysis is needed here across graphs. Here too, I would start by looking at specific trajectories (and please present sonograms).

Thank you we have addressed this comment in several ways. First, we have recoded the songwindow data as a continuous metric (amount of time that the song is plastic, described above), and performed a PGLS (Phylogenetic Generalized Least Squares regression) test to detect evolutionary correlations in continuous characters. For syllable repertoire and song repertoire, we found a significant correlation between the length of the learning window and repertoire size when controlling for phylogenetic relatedness. We also plotted the continuous data in Figure 1D, so that the reader can evaluate it. This figure illustrates that learning windows of less than one year are primarily associated with repertoires of 10 or fewer syllables; repertoire sizes around one year are compatible with a range of repertoire sizes below ~50 syllables, and learning during and after the second breeding season corresponds to a wide range of possible repertoire sizes, mostly between ~20 and ~500 syllables.

We also examined whether the decision to break species into 3 bins was influencing our results by also dividing the data into 2, 4, and 5 bins. We saw similar results regardless of the number of bins used, where the smallest repertoire sizes with song plasticity were evolutionarily unstable, and the largest repertoires with song stability were evolutionarily unstable. For bins in which the threshold was set in the intermediate range, we also saw similar patterns for syllable repertoire, though it does appear that in the range of ~1460 syllables, there is little evolution between song stability states. Perhaps this range represents an equilibrium point where birds can evolve their repertoire sizes independently of their song stability. This concurs with the results of the continuous analysis noted above: it would make sense that we do not see a high rate of evolution of the length of the learning window in this threshold range, since repertoires of ~1460 syllables are common in both species with learning windows lasting around one year (delayed stability) and those with learning windows lasting two or more years (song plasticity) These results suggest that the evolutionary transitions in repertoire size within this range does not necessarily require alteration of the songlearning window length. For song repertoire, the intermediate bins show that species with smaller song repertoires tend to transition from song plasticity to song stability but there is little change in song repertoire sizes.

We also conducted another test (Author response image 2) and can include in the main text if desired.

Regarding the possibility of an equilibrium state, we made note of an interesting change in repertoire size stability when we repeated the BayesTraits analyses with a different grouping, separating species with early song stability (shorter learning windows) from those with either delayed song stability or adult song plasticity (longer learning windows). In these analyses, we observed that the rate of transition between shorter and longer learning windows was reduced

compared to when species were grouped as songstable and songplastic, hinting that the predominant evolutionary transitions might be between species with delayed song stability and adult song plasticity. The “longer learning” grouping helped stabilize species with the smallest repertoires and longer leaning by reducing the transition to early song stable. Thus, we hypothesize that delayed song stability can allow species with small repertoires to maintain some song flexibility, but learning extending across multiple years is unstable.

6) Throughout the paper I kept wondering if there is a simple way of telling which came first: song learning duration or syllable per song? Can you look at specific trajectories to get some clues, even if not statistical?

If there were a consistent pattern of which came first, song learning duration or song complexity, we should see a clear trend with the BayesTraits analysis, which would appear as a single directional arrow away from an unstable combination of traits. Our results were more nuanced: evolutionarily unstable combinations, such as very small repertoires with adult song plasticity or very large repertoires with adult song stability, were able to transition in either trait; one transition did not always precede the other. However, we do report a trend in these rates: the transitions between adult song stability and plasticity appear to occur at a higher rate than transitions in repertoire size, perhaps indicating that the song learning window is more likely to change first, and song complexity changes in response. We now note this trend in the text.

7) Also, since mimids drive much of the correlations maybe it is a good idea to look at them more carefully as a post hoc analysis. I understand that authors want to rid of unrepresentative example, but this does not mean that the most convincing cases should be ignored.

We would like to emphasize that mimids only affected the patterns in syllables per song, not in syllable repertoire size, song repertoire size, or any song performance trait. We note in the text that syllables per song might be a metric that is especially susceptible to differences in listener perception: where does one mockingbird song end and another begin?

We also note that with mimids included, we found an increased rate of evolution of syllables per song, but there was no change in the direction of evolution, and the results regarding the rate of evolution relied on the inclusion of two species in particular, Mimus gilvus and M. polyglottis. Removing M. polyglottis made the rate of evolution nonsignificantly different between songstable and songplastic lineages (p = 0.074), while removal of M. gilvus made the difference very nonsignificant (p = 0.504). Thus, more research is required before knowing which aspects of this result reflect real biological differences and which aspects hinge on the coding of syllables per song for particular mimids. Furthermore, in our new continuous analysis (PGLS), syllables per song was not significantly correlated with the length of the songlearning window.

In total, this suggests that while mimids warrant further investigation, they are acting differently than other openended learners, and it may be that M. gilvus specifically is acting differently from other bird species or that its song is particularly difficult to parse accurately into syllables per song.

8) Overall, it would be good to show, at least in the supplement, some examples of evolutionary of complexity and performance side by side with evolution of open-ended learning. I feel that without seeing what's going on in specific cases it is difficult to judge the quality of the data.

Thank you for this suggestion. First, we apologize if our datasets with all of the specieslevel data were not transmitted properly we fully intended for the reviewers to be able to judge the data on song characteristics and learning window and the sources where we obtained those data. As a precaution, we have uploaded all data and code to github (github.com/CreanzaLab/SongLearningEvolution) in addition to the journal website.

In addition, so that readers can evaluate any specific evolutionary trajectory, we now have a supplemental section in which we provide the estimated ancestral values of both song characteristics and learning window. To do this, we numbered every node in the phylogenies and provide a table with song trait and plasticity/stability estimates for each node (Figure 2—source data 1). Thus, the evolution of complexity and performance is shown side by side with the evolution of openended learning for each of the song characteristics. We think that this complements the colorcoding of the tree branches in the main text, providing a quantitative estimate of each evolutionary trajectory.

9) Regarding data quality of adult song plasticity: authors used a categorical measure (taken from multiple published papers) and there is no mention in the Materials and methods about assessment and quality control of integrating such a categorical assessment across studies. Song traits, on the other hand, are continuous variables. I think that this categorical, binary assessment of song plasticity has limitation that should be discussed.

We have addressed this comment both by explicitly describing our categorization methods in the text and by reanalyzing our data in two new ways, multistate discrete and continuous, as described above.

10) Line numbers would make reviewing easier.

We apologize for this oversight! The current version has line numbers.

11) It is not clear whether the data included were all from males in each species, or whether some females were included. If a mixture of sexes, I think this factor should controlled for in one or more analyses. The authors mentioned females in the discussion for future studies, but never mentioned (maybe assumed) that they only analyzed studies with males in this work.

Thank you – we analyzed male song and song learning without explicitly stating that we did so. We have added text to the Materials and methods to clarify this point.

12) The authors should make a cautionary statement that simulations may not match the real biology, but approximate it.

We now include a sentence to this effect in the Results section.

13) The Introduction seems excessively long on background material and different hypotheses. I suggest dealing with the various hypothesis in the Results, and only mention two in the Introduction, so that the paper is less repetitive.

We have streamlined the Introduction and integrated some of the text into the Results section.

14) Also in the Introduction and other parts of the paper, stability and plasticity are treated as different variables (whether in written text or some of the figures). But they are along a continuum, that the authors results support. I suggest making this clearer.

Please see our third response above.

15) On its first mention, the authors need to cite what phylogenetic tree of oscine families that they used. The accuracy of the tree could affect their results. I couldn't find any citation in the main text for the tree used.

Apologies for this missing citation the tree citation (Jetz et al., 2012) is now mentioned in the Results and Materials and methods sections.

16) About the tree results, it looks pretty clear from all the phylogenetic analyses that there was an early split or two splits in the songbird lineage between two predominantly plastic song clades and a predominantly stable song clade. And then some independent deviations from those early splits moving into the other clade's phenotype. Maybe the common songbird ancestor was in the middle of the continuum. This exciting finding of the major clades with plastic vs. stable song should be highlighted in the Abstract and Results.

Thank you for this suggestion! We now highlight this point in Results and Discussion.

17) Subsection “Do our results depend on the particular values or families included in analyses?”, first paragraph. I thought other song production factors, like song intervals, were not affected in the above analyses? But they are here, when using maximum or minimum values?

Song production/performance characteristics were not significantly different when we compared species with adult song plasticity versus species with adult song stability, but some of these characteristics did appear to evolve significantly faster in songplastic lineages. To be stringent, we repeated all analyses with the minimum and maximum values that we found in the literature to ensure that our results did not hinge on the exact value we used (as in Moore et al., 2011).

18) The Discussion is on the long side. It could be shortened without losing content, by not repeating the results, but simply moving into a discussion of them.

We have revised the Discussion to shorten it, by eliminating repeated results, combining semirepetitive and by streamlining the rest of the text.

19) The following two sentences appear to contradict each other: "We did not findevidence for correlated evolution between song stability and EPP. Most of our simulations testing for correlated evolution between mating system and learning window were significant;"

We apologize for the confusing wording here. We were using “mating system” as a general term for monogamy and polygyny, not including EPP. We have rephrased this section and now refer to mating system as “social mating system” to clarify.

20) Discussion, seventh paragraph. Explain what is HVC, to scientists unfamiliar with the avian brain.

We have clarified the wording of this sentence to define HVC.

21) Figure 1: Label y axes (log scale). Mention stat results in the legend.

Thank you, we have made the suggested changes to the figure and its caption.

22) Figure 2. Legend text is hard to follow, in terms of grammar.

We have revised the caption for this figure.

23) Please consider archiving the data in a proper repository so that they can be more readily built upon.

We apologize that the raw data was not available to the reviewers. For our initial submission of the manuscript, we uploaded the database of information that we collected and analyzed as supplementary tables (Supplementary Datasets S1 and S2), and we uploaded the raw data and code as a zipped archive. We intend for both of these (Supplemental Datasets and Supplemental Code) to be supplemental information that would be available alongside our paper on the journal’s website at publication. In case this supplemental data and code are not provided to reviewers, we have also uploaded them to github.com/CreanzaLab/SongLearningEvolution. In response to this comment, we have also added several files to our supplementary information (e.g. Figure 2—source data 1) so that the raw reconstructed evolutionary states at each node are available in conjunction with the phylogenies.

24) Abstract: Song-plastic species. This is not a term I am used to. Do you mean open ended learners? If so, can you use that term? Or at least consider "species with plastic song" – and perhaps briefly define what you mean.

25) Abstract: "Notably…": What is the evidence for this? Is this derived from the above result, or is it from other results. Please preface this with what kind of analysis demonstrated this result.

Thank you for these suggestions; we have rephrased the Abstract accordingly.

26) Introduction – "small-scale analysis of comparative evolution": more appropriate = "small-scale comparative analysis".

We have made this change.

27) Introduction: "female choice can favor certain characteristics" – can you give a bit more description of process, or what is selected for (similar to the description for neuroplasticity). This is crucial to your findings.

We have reworded this section with more details to better parallel the description of neuroplasticity.

28) Introduction, third paragraph – I would expect this before talking about song-learning windows – as it is the more standard dogma under which bird song is studied.

We have edited the Introduction and we hope that it flows better.

29) Results – first paragraph: this is all good information, but does not present any results. It seems better suited for the Materials and methods (and, in fact, this useful explanation of why you went with the novel, broad "song-plastic/stable" terminology is left out of the Materials and methods. I would recommend moving more of this text there and skip straight to your classification and Results (with a brief definition here)).

Thank you! We moved this section to the Materials and methods and instead start the Results with a short definition of our classification systems.

30) Results – second paragraph: the idea of examining the ancestral state of all songbirds for song plasticity is interesting, however, you don't really have an appropriate dataset for this. Because you have so many missing taxa (only 64 species out of ~5000 categorized) and they fall throughout the songbirds (you could predict the ancestral state for a particular family if you had sampled many representatives of that family, or all songbirds if you densely sampled representatives of lineages branching more directly from the common ancestor of all songbirds). However, your dataset is a good one for examining correlated evolution, as you have done, and you can also get an idea of transition rates and when/ how often transitions between song stability and plasticity occur (although, it is good to be aware that this may be a fairly rough estimate just because there are probably many species that are missing among the species sampled). I would consider removing/avoiding mention of ancestral states, and approach this in light of major transitions among these traits. You could mention that the ancestral state for these lineages could not be resolved.

We no longer focus on the ancestral state of all birds other than to mention that we could not resolve it with our current dataset.

31) Results – paragraph beginning "Unexpectedly" – I'm not sure this is that unexpected (I'm not certain what I would predict here). Perhaps skip this word.

We have altered this phrasing.

32) Results – subheading "Does song stability or plasticity affect the direction of evolution for song traits and mating system?": seems to me the predicted direction of causality is backwards here. I would predict that song-learning strategy is dictated by ecological traits or the characteristics of the songs that are being selected (as you find). Therefore, I suggest alternate subheading and approaching this discussion from the perspective of "Is song stability or plasticity influenced by the evolution of song traits or mating system?"

I think this is also approached in the same way in the Abstract/Introduction. There I also suggest pitching this the opposite way (that's how I automatically think of it).

We agree and have rephrased the text in the Abstract and Results to reflect this idea.

33) Discussion – first paragraph: "in concert with sexually selected song features" – we don't know that these features are strictly sexually selected. Consider "song features often thought of as sexually selected" or "song features associated with sexual selection".

We have changed this phrasing as suggested; thank you.

34) Discussion – third paragraph: the two sentences after "We observed that with song stability….". Can this be explained more simply? A lot of jargon here (repeller/attractor states). Can you explain these trends in terms of what is evolutionarily stable/unstable – e.g. song stability coupled with small syllables and song repertoires is an evolutionarily unstable state, so transitions away from it are high.

We have rephrased the relevant sections using the suggested terms. When drafting the manuscript, we originally used “evolutionarily stable/unstable,” but then we grew concerned that talking about the evolutionary stability of song stability would become confusing.

35) Discussion – fourth paragraph: I think you undersell the importance of mimics (not strictly mimids) in songbird song evolution here, and the implications of your result including the mimids. Mimicry is common throughout songbirds – including in the earliest-branching songbird lineage (the lyrebirds). Really the sample size of the study is small compared to all songbirds (understandably limited by our knowledge of the song learning window for most species). However, I think it would be nice to acknowledge here or elsewhere that this is a relatively small sample, but that mimicry is common in songbirds, so the pattern seen in the mimids may be meaningful and it will be interesting to see how it holds up as we learn more.

We now address explicitly in the Discussion the interesting and important questions surrounding the evolution of mimicry and its interaction with the songlearning window. To shorten the Discussion, though, we moved some of the previous discussion of our Mimidaespecific findings back to the relevant section in the Results.

36) Discussion – overall comment – there are a lot of thoughtful points here, but the Discussion is rather long. I wonder if these can be delivered more briefly. Perhaps some points can be combined, or overall length cut?

Please see our response to comment #18.

37) Materials and methods – first paragraph: "adult song plasticity ('song-stable species')" – I think there is a typo here – doesn't this statement refer to species without song plasticity?

Thank you for catching this confusing typo! We have corrected this in the text.

38) Materials and methods – definitions for 'Song-stable/plastic': there seems to be a gap in your definitions – which leaves birds that learn prior to their second breeding season but do not subsequently modify their song following the second breeding season. I am under the impression that this is not an uncommon strategy. How did you handle these species?

Please see our response to comment #3.

39) Materials and methods: defined high EPP when at least 10% of offspring had a genetic male parent other than social male. This seems low to me. What was the justification for this cutoff?

We now provide the justification for this cut-off in the text: first, a review of extrapair paternity studies (Griffith et al., 2002) estimated the crossspecies average to be ~11% of offspring per nest to be attributable to extrapair mates, so we wanted ‘low’ EPP to be below average. In addition, previous analyses used a 10% threshold (e.g. Soma and Garamszegi, 2011), and we wanted our study to be in line with previous efforts to delineate low/high EPP.

https://doi.org/10.7554/eLife.44454.051

Article and author information

Author details

  1. Cristina M Robinson

    Department of Biological Sciences, Vanderbilt University, Nashville, United States
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing
    Contributed equally with
    Kate T Snyder
    Competing interests
    No competing interests declared
  2. Kate T Snyder

    Department of Biological Sciences, Vanderbilt University, Nashville, United States
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing
    Contributed equally with
    Cristina M Robinson
    Competing interests
    No competing interests declared
  3. Nicole Creanza

    Department of Biological Sciences, Vanderbilt University, Nashville, United States
    Contribution
    Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Validation, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing
    For correspondence
    nicole.creanza@vanderbilt.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8821-7383

Funding

No external funding was received for this work

Senior and Reviewing Editor

  1. Patricia J Wittkopp, University of Michigan, United States

Reviewer

  1. Ofer Tchernichovski, Hunter College, United States

Version history

  1. Received: December 17, 2018
  2. Accepted: July 16, 2019
  3. Version of Record published: September 3, 2019 (version 1)
  4. Version of Record updated: September 13, 2019 (version 2)

Copyright

© 2019, Robinson et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 2,462
    Page views
  • 273
    Downloads
  • 24
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Cristina M Robinson
  2. Kate T Snyder
  3. Nicole Creanza
(2019)
Correlated evolution between repertoire size and song plasticity predicts that sexual selection on song promotes open-ended learning
eLife 8:e44454.
https://doi.org/10.7554/eLife.44454

Further reading

    1. Computational and Systems Biology
    Zeyneb Kurt, Jenny Cheng ... Xia Yang
    Research Article

    Mouse models have been used extensively to study human coronary artery disease (CAD) or atherosclerosis and to test therapeutic targets. However, whether mouse and human share similar genetic factors and pathogenic mechanisms of atherosclerosis has not been thoroughly investigated in a data-driven manner. We conducted a cross-species comparison study to better understand atherosclerosis pathogenesis between species by leveraging multiomics data. Specifically, we compared genetically driven and thus CAD-causal gene networks and pathways, by using human GWAS of CAD from the CARDIoGRAMplusC4D consortium and mouse GWAS of atherosclerosis from the Hybrid Mouse Diversity Panel (HMDP) followed by integration with functional multiomics human (STARNET and GTEx) and mouse (HMDP) databases. We found that mouse and human shared >75% of CAD causal pathways. Based on network topology, we then predicted key regulatory genes for both the shared pathways and species-specific pathways, which were further validated through the use of single cell data and the latest CAD GWAS. In sum, our results should serve as a much-needed guidance for which human CAD-causal pathways can or cannot be further evaluated for novel CAD therapies using mouse models.

    1. Computational and Systems Biology
    2. Microbiology and Infectious Disease
    Vanessa Dumeaux, Samira Massahi ... Michael T Hallett
    Research Article Updated

    Candida albicans, an opportunistic human pathogen, poses a significant threat to human health and is associated with significant socio-economic burden. Current antifungal treatments fail, at least in part, because C. albicans can initiate a strong drug tolerance response that allows some cells to grow at drug concentrations above their minimal inhibitory concentration. To better characterize this cytoprotective tolerance program at the molecular single-cell level, we used a nanoliter droplet-based transcriptomics platform to profile thousands of individual fungal cells and establish their subpopulation characteristics in the absence and presence of antifungal drugs. Profiles of untreated cells exhibit heterogeneous expression that correlates with cell cycle stage with distinct metabolic and stress responses. At 2 days post-fluconazole exposure (a time when tolerance is measurable), surviving cells bifurcate into two major subpopulations: one characterized by the upregulation of genes encoding ribosomal proteins, rRNA processing machinery, and mitochondrial cellular respiration capacity, termed the Ribo-dominant (Rd) state; and the other enriched for genes encoding stress responses and related processes, termed the Stress-dominant (Sd) state. This bifurcation persists at 3 and 6 days post-treatment. We provide evidence that the ribosome assembly stress response (RASTR) is activated in these subpopulations and may facilitate cell survival.