Phylogenetic Logistic Regression Research Articles

Citizen science data reveal altitudinal movement and seasonal ecosystem use by hummingbirds in the Andes Mountains

Ensuring connectivity is crucial to protect landscapes but it requires knowledge about how animals use ecosystems throughout the year. However, animal movements remain largely unknown in biodiversity hotspots, even for species that fulfill key ecological roles, as is the case of hummingbirds in the Andes. In the complex topography of mountain slopes, movement of these avian pollinators may occur either between habitat patches with asynchronous plant blooms or across ecosystems that are located within the same elevation bands or along altitudinal gradients. Here, we used two decades (2000–2020) of records from citizen science data and boosted regression trees to predict monthly distributions for 55 hummingbird species in the Andes. We identified shifts in altitudinal distribution between contiguous months and calculated changes in the proportion of predicted distributions occupied by ecosystem types. Our findings reveal substantial altitudinal movement and differences in the proportion of ecosystem types utilized throughout the year that had not been previously reported for several species. Yet the magnitude of altitudinal and ecosystem shifts varies between hummingbird clades, and in some cases changes in the proportion of ecosystem types within estimated distributions occurs with little variation in altitude. All ecosystems across the Andes show temporal changes in hummingbird occurrence, but these are higher in natural landscapes compared to croplands or urban areas. Finally, we used phylogenetic logistic regression to test whether altitudinal and ecosystem shifts affect population trends. We found that higher ecosystem seasonality is more strongly associated with decreasing populations in comparison to altitudinal shifts. Altogether, our study reveals complex patterns of movement in hummingbirds and highlights the importance of ecological connectivity across different ecosystem types. More generally, it demonstrates the opportunity of using citizen science data to increase understanding about species' seasonal occurrences, so that landscapes can be better managed to protect animal movement.Keywords: boosted regression trees, eBird, ecological connectivity, species distribution models

Evolutionary Patterns of Intersexual Power

Dominance and leverage are both possible causes of social inequality. If sexual dimorphism influences patterns of intersexual dominance, we predicted that highly dimorphic species are constrained to exhibit male-biased power (MP), but species with low sexual dimorphism are free to demonstrate a broader range of intersexual power relationships. If market effects influence intersexual leverage, we predicted that females have more power when group composition is more male-biased and estrus is asynchronous. We analyzed data on intersexual power, sexual dimorphism, expected estrous overlap, and sex ratio for 79 extant primate species using phylogenetic logistic regression and ancestral state reconstructions. Although MP is more common, every major primate clade includes non-MP species. MP was associated with greater body mass and canine length dimorphism and with female-biased sex-ratios. Low estrous overlap was associated with non-MP. Although MP was reconstructed as likely ancestral for anthropoids, the last common ancestor of this clade probably did not exhibit high sexual dimorphism. The last common ancestor of catarrhines was probably highly dimorphic, potentially constraining intersexual power relationships. Non-MP probably evolved multiple times in primates and may be less common because multiple traits are linked to MP while fewer traits are associated with female-biased power or equality.

Inferring the lifestyles of extinct Crocodyliformes using osteoderm ornamentation.

Osteoderms are bony plates formed within the dermis of diverse vertebrate groups. They are present in all crocodylomorphs but Metriorhynchidae. Most of them show typical bone ornamentation consisting of pits and ridges on their outer surface. The most widely discussed functional hypothesis suggests that the ornamentation of osteoderms influences heat exchange with the environment through the adjacent vascular network, facilitating the absorption of solar radiation. This process allows semiaquatic crocodiles to compensate for heat loss resulting from the high thermal conductivity of surrounding water. In order to test this assertion, we conducted a phylogenetic logistic regression analysis to evaluate the relationship between osteoderm relative area of pits (RAP) and lifestyle (terrestrial versus aquatic) in a sample of crocodyliforms. Our results revealed that lifestyle is significantly explained by RAP: the lower the degree of ornamentation (RAP), the higher the probability of a terrestrial lifestyle. We used this model to infer the lifestyle of two extinct taxa, Peirosaurus torminni and Microsuchus schilleri. We concluded that terrestrial notosuchians may have lost osteoderm ornamentation due to the lower thermal conductivity of air and reduced heat loss in a terrestrial environment compared to what happens in water. Among these notosuchians, we hypothesize that large terrestrial baurusuchids maintained a stable body temperature due to thermal inertia, whereas small notosuchians took advantage of the early morning sun exposure to warm up and stayed in terrestrial burrows during periods of intense solar radiation. Finally, unlike the almost motionless behavior of freshwater crocodiles, fully marine Metriorhynchidae probably lost osteoderms because they constantly swim, generating heat by muscular contraction, so osteoderms with a thermoregulatory function for heat absorption were no longer positively selected.

Apomixis occurs frequently along the entire American Cordillera

Abstract Apomictic plants mostly occur at higher latitudes and higher elevations, where they tend to occupy deglaciated areas. By analysing samples of the alpine floras of different latitudinal zones of the Americas (Argentina, Bolivia, Ecuador, USA), we investigated whether the proportion of gametophytic apomixis is comparable between the temperate alpine zones of the Northern Hemisphere and the Southern Hemisphere and whether the proportion of apomictic species increases from the alpine tropics towards the temperate alpine zone. The reproductive systems of 196 species and 275 individuals were determined by a flow-cytometric seed screening of 999 seeds. Phylogenetic logistic regression was used to evaluate differences in the proportion of apomictic species in samples of regional alpine floras. The apomictic embryo:endosperm ratio was determined for 33 species (two pseudogamous species and 31 species with autonomous apomixis) belonging to 28 genera of 13 families. The probability of a species being apomictic did not significantly differ between the regions. The highest probability of a species being apomictic was estimated for the northern temperate zone (29.6%), followed by the tropical zone (12.5%), and the southern temperate zone (11.9%). Our results show that asexual plant reproduction by apomixis is not restricted to the alpine zone of the Northern Hemisphere and its representation in alpine floras of the Southern Hemisphere, including the alpine tropics, has been substantially underestimated.

Why are animals conspicuously colored? Evolution of sexual versus warning signals in land vertebrates.

Conspicuous colors (e.g., red, yellow, blue) have evolved numerous times across animals. But the function of this coloration can differ radically among species. Many species use this coloration as a sexual signal to conspecifics, whereas others use it as a warning signal to predators. Why do different species evolve conspicuous coloration in association with one function as opposed to the other? We address this question in terrestrial vertebrates (tetrapods) using phylogenetic approaches, and test whether day-night activities of species help determine these patterns. Using phylogenetic logistic regression, we found that conspicuous, sexually dimorphic coloration is significantly associated with diurnal lineages (e.g., many birds and lizards). By contrast, the evolution of warning signals was significantly associated with large-scale clades that were ancestrally nocturnal (e.g., snakes, amphibians), regardless of the current diel activity of species. Overall, we show that the evolution of conspicuous coloration as warning signals or sexual signals is influenced by the ecology of species, both recently and in the ancient past.

Neotropical Birds Respond Innately to Unfamiliar Acoustic Signals.

AbstractMany animals respond to heterospecific signals that indicate the presence of food or predators. Although the benefits of responding are clear, the behavioral and cognitive mechanisms underlying responses are not. Whether responses are learned, innate, or an epiphenomenon created by following other species as they respond to signals remains unknown because most studies have involved respondents that are sympatric with their heterospecific signalers and that have therefore had opportunities to learn their signals. In this study, we tested the mechanisms underlying avian responses to heterospecific chick-a-dee calls. All North American parids produce chick-a-dee calls in response to arousing stimuli, such as food and predators, and diverse species respond by approaching the caller and consuming the food or mobbing the predator. We broadcast chick-a-dee calls plus two control stimuli in Costa Rica, Colombia, and Brazil, where no parids ever occur. We conducted our trials in the winter, when Neotropical migrants that might be familiar with chick-a-dee calls were present, and in the temperate breeding season, when migrants were absent. Across 138 trials, 38 resident species from 14 families and four orders responded to chick-a-dee calls by approaching to within 5 m of the playback speaker. A phylogenetic logistic regression showed that whether a species responded was not significantly associated with the species' mean body mass or the structural similarity between its calls and chick-a-dee calls. Residents were significantly more likely to approach chick-a-dee calls than either control stimulus. This pattern was unaffected by the presence of migrants, thus demonstrating that the observed responses are innate. Our study shows that learning cannot fully explain responses to heterospecific chick-a-dee calls and that structural features distinguishing these calls from other vocalizations are important.

Sociality and migration predict hybridization across birds.

Hybridization represents a natural experiment that can provide insight into processes of speciation and diversification. Though considerable research has focused on hybrid zone dynamics, macroevolutionary investigations of the factors that influence hybridization are few. Here, we compile a database of avian hybrids and perform comparative analyses to determine whether several social and life-history variables influence broad patterns of hybridization. We perform three main analyses: phylogenetic logistic regression to examine variables that are associated with the presence of hybridization, phylogenetic Poisson regression of only those species exhibiting hybridization to determine which variables are associated with the extent of hybridization, and a phylogenetic logistic regression on a subset of data to assess potential pseudoreplication. After testing several social and life-history variables, we found that social bond duration is associated with the presence and extent of hybridization. Specifically, lengthy social bonds are negatively associated with the presence and extent of hybridization. In addition to social bond length, migration is positively linked with a greater likelihood of hybridization. The broad-scale correlations between species-specific traits and hybridization across diverse avian lineages suggest commonalities in the fine-scale processes involved in mating with heterospecifics, which in turn has implications for how we think about, study and understand hybridization processes and their influence on evolutionary trajectories.

Invertebrate‐mediated dispersal plays an important role in shaping the current distribution of a herbaceous monocot

AbstractAimInvertebrate‐mediated dispersal has previously been proposed to promote angiosperm diversification and distribution. However, little is known about the specific impact of invertebrate‐mediated dispersal on the biogeography and current distribution of plants. We aim to infer the influence of vespicochorous (hornet) and myrmecochorous (ant) dispersal on the historical biogeography of herbaceous monocot species.LocationSoutheast Asia, East Asia, Australia, North America.TaxonFamily Stemonaceae.MethodWe sampled ca. 75% of the species diversity in Stemonaceae (28 out of 37 species), covering the entire distribution range of the family, to reconstruct the biogeographic history of this family. Using phylogenetic logistic regression analyses, we then tested the relationship between dispersal modes and geographic distributions.ResultsStemonaceae originated on the Asian mainland during the late Cretaceous and then dispersed to North America, the western Malay Archipelago and eastern Malay Archipelago and Australia between the late Cretaceous and Pliocene. Geographical ranges of ant‐ versus hornet‐dispersed Stemonaceae species are significantly different, with vespicochorous species having broader distribution ranges than myrmecochorous species.Main conclusionsInvertebrate‐mediated dispersal in Stemonaceae may promote narrow endemism and play an important role in shaping the current distribution of species. Most lineages dispersed by ants failed to cross biogeographic barriers and exhibit limited range expansion overland. Vespicochorous lineages were able to cross oceanic barriers and occupy larger continental areas and/or occur on oceanic islands.

Comparative genomic analysis of Vibrios yields insights into genes associated with virulence towards C. gigas larvae

BackgroundVibriosis has been implicated in major losses of larvae at shellfish hatcheries. However, the species of Vibrio responsible for disease in aquaculture settings and their associated virulence genes are often variable or undefined. Knowledge of the specific nature of these factors is essential to developing a better understanding of the environmental and biological conditions that lead to larvae mortality events in hatcheries. We tested the virulence of 51 Vibrio strains towards Pacific Oyster (Crassostreae gigas) larvae and sequenced draft genomes of 42 hatchery-associated vibrios to determine groups of orthologous genes associated with virulence and to determine the phylogenetic relationships among pathogens and non-pathogens of C. gigas larvae.ResultsV. coralliilyticus strains were the most prevalent pathogenic isolates. A phylogenetic logistic regression model identified over 500 protein-coding genes correlated with pathogenicity. Many of these genes had straightforward links to disease mechanisms, including predicted hemolysins, proteases, and multiple Type 3 Secretion System genes, while others appear to have possible indirect roles in pathogenesis and may be more important for general survival in the host environment. Multiple metabolism and nutrient acquisition genes were also identified to correlate with pathogenicity, highlighting specific features that may enable pathogen survival within C. gigas larvae.ConclusionsThese findings have important implications on the range of pathogenic Vibrio spp. found in oyster-rearing environments and the genetic determinants of virulence in these populations.

The relative importance of phylogeny and habitat in determining the presence and prominence of a granula iridica in hooved mammals.

To investigate the relationship between phylogeny and amount of shade in a species' habitat regarding the presence or absence of an iridal granula iridica (GI) in a large sample of Artiodactyl and Perissodactyl clades and using online resources. The Comparative Ocular Pathology Laboratory of Wisconsin (COPLOW) archives were searched for glass slide material from Artiodactyl (even-toed) and Perissodactyl (odd-toed) ungulates. The slides were examined, and the presence or absence of the GI was noted. The phylogenetic tree of the ungulate species was inferred using TimeTree (http://www.timetree.org), and the habitat data are derived from Animal Diversity Web (https://animaldiversity.org/). We assessed the probability of the presence of GI occurring given the amount of shade in a species' environment using phylogenetic logistic regression. Forty-eight artiodactyl species were able to be evaluated and tabulated. Nine perissodactyl species were able to be evaluated. The phylogenetic logistic regression showed that the probability of GI presence was lower in artiodactyl species that inhabited shaded environments (βshaded =-1.774). Arctiodacyl species inhabiting a nonshaded environment were slightly more probable to have the GI present (βnonshaded =0.023), with species inhabitating ambiguously shaded environments having a high probability of GI presence (βambiguous =2.214). Our results suggest that the GI may be a common morphological feature to shade the pupil in nonshaded environments, and, in its absence, increase the amount of light reaching the retina to improve vision in shaded environments for hooved mammals. Further research on the functional optics of the GI and studies that include additional ungulate species would further elucidate phylogenetic and ecological factors influencing the occurrence of GI in hooved mammals.

Expansion of vomeronasal receptor genes (OlfC) in the evolution of fright reaction in Ostariophysan fishes

Ostariophysans are the most diverse group of freshwater fishes and feature a pheromone-elicited fright reaction. However, the genetic basis of fright reaction is unclear. Here, we compared vomeronasal type 2 receptor-like (OlfC) genes from fishes having and lacking fright reaction, to provide insight into evolution of pheromonal olfaction in fishes. We found OlfC genes expanded remarkably in ostariophysans having fright reaction compared with fishes lacking fright reaction. Phylogenetic analysis indicates OlfC subfamily 9 expanded specifically in ostariophysans having fright reaction. Principle component and phylogenetic logistic regression analysis partitioned fishes by ecotype (having or lacking fright reaction) and identified OlfC subfamily 9 as being an important factor for fright reaction. Expression levels of expanded OlfC subfamily genes after fright reaction in zebrafish changed more than did genes that had not expanded. Furthermore, evidence of positive selection was found in the expanded OlfC proteins in ostariophysan fishes having fright reaction. These results provide new insight into the genetic basis of fright reaction in ostariophysan fish and will enable future research into the mechanism of action of OlfC proteins.

Sodium-powered stators of the bacterial flagellar motor can generate torque in the presence of phenamil with mutations near the peptidoglycan-binding region.

The bacterial flagellar motor powers the rotation that propels the swimming bacteria. Rotational torque is generated by harnessing the flow of ions through ion channels known as stators which couple the energy from the ion gradient across the inner membrane to rotation of the rotor. Here, we used error-prone PCR to introduce single point mutations into the sodium-powered Vibrio alginolyticus/Escherichia coli chimeric stator PotB and selected for motors that exhibited motility in the presence of the sodium-channel inhibitor phenamil. We found single mutations that enable motility under phenamil occurred at two sites: (i) the transmembrane domain of PotB, corresponding to the TM region of the PomB stator from V. alginolyticus and (ii) near the peptidoglycan binding region that corresponds to the C-terminal region of the MotB stator from E. coli. Single cell rotation assays confirmed that individual flagellar motors could rotate in up to 100µM phenamil. Using phylogenetic logistic regression, we found correlation between natural residue variation and ion source at positions corresponding to PotB F22Y, but not at other sites. Our results demonstrate that it is not only the pore region of the stator that moderates motility in the presence of ion-channel blockers.

Correlated Evolution of White Spots on Ears and Closed Habitat Preferences in Felids

The pigmentation patterns of many carnivorous mammals comprise contrasting white patches of hair in different parts of the body whose evolution remains largely misunderstood. Some felids (Felidae) exhibit conspicuous white spots on the posterior part of the ears, while the ear color of others is uniform. On the basis that ear movement in felids has a role in intraspecific communication and that color contrast enhances detection, here it is hypothesized that white spots on ears may be particularly adaptive under conditions of poor visibility and thus be associated with the occupancy of closed habitats. This prediction was tested using phylogenetic logistic regression models with all species of extant felids. Results show a clear association between the occurrence of white spots on ears and preference for closed habitats, and this is independent of body size and whether species that occupy both closed and open habitats are considered as closed- or open-habitat specialists. Phylogenetic signal analyses indicate that the occurrence of white spots on ears is a highly conserved trait while habitat preferences are evolutionarily labile, suggesting that the presence of white spots may have partly contributed to the adaptation of felids to closed habitats. These findings indicate that some subtle pigment traits have fulfilled a significant role in determining the success of habitat occupancy by felids and possibly other mammals, which in turn has driven the evolutionary maintenance of such traits.

R$^{2}$s for Correlated Data: Phylogenetic Models, LMMs, and GLMMs

Many researchers want to report an $R^{2}$ to measure the variance explained by a model. When the model includes correlation among data, such as phylogenetic models and mixed models, defining an $R^{2}$ faces two conceptual problems. (i) It is unclear how to measure the variance explained by predictor (independent) variables when the model contains covariances. (ii) Researchers may want the $R^{2}$ to include the variance explained by the covariances by asking questions such as "How much of the data is explained by phylogeny?" Here, I investigated three $R^{2}$s for phylogenetic and mixed models. $R^{2}_{resid}$ is an extension of the ordinary least-squares $R^{2}$ that weights residuals by variances and covariances estimated by the model; it is closely related to $R^{2}_{glmm}$ presented by Nakagawa and Schielzeth (2013. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 4:133-142). $R^{2}_{pred}$ is based on predicting each residual from the fitted model and computing the variance between observed and predicted values. $R^{2}_{lik}$ is based on the likelihood of fitted models, and therefore, reflects the amount of information that the models contain. These three $R^{2}$s are formulated as partial $R^{2}$s, making it possible to compare the contributions of predictor variables and variance components (phylogenetic signal and random effects) to the fit of models. Because partial $R^{2}$s compare a full model with a reduced model without components of the full model, they are distinct from marginal $R^{2}$s that partition additive components of the variance. I assessed the properties of the $R^{2}$s for phylogenetic models using simulations for continuous and binary response data (phylogenetic generalized least squares and phylogenetic logistic regression). Because the $R^{2}$s are designed broadly for any model for correlated data, I also compared $R^{2}$s for linear mixed models and generalized linear mixed models. $R^{2}_{resid}$, $R^{2}_{pred}$, and $R^{2}_{lik}$ all have similar performance in describing the variance explained by different components of models. However, $R^{2}_{pred}$ gives the most direct answer to the question of how much variance in the data is explained by a model. $R^{2}_{resid}$ is most appropriate for comparing models fit to different data sets, because it does not depend on sample sizes. And $R^{2}_{lik}$ is most appropriate to assess the importance of different components within the same model applied to the same data, because it is most closely associated with statistical significance tests.

Breeding biology and the evolution of dynamic sexual dichromatism in frogs

Dynamic sexual dichromatism is a temporary colour change between the sexes and has evolved independently in a wide range of anurans, many of which are explosive breeders wherein males physically compete for access to females. Behavioural studies in a few species indicate that dynamic dichromatism functions as a visual signal in large breeding aggregations; however, the prevalence of this trait and the social and environmental factors underlying its expression are poorly understood. We compiled a database of 178 anurans with dynamic dichromatism that include representatives from 15 families and subfamilies. Dynamic dichromatism is common in two of the three subfamilies of hylid treefrogs. Phylogenetic comparative analyses of 355 hylid species (of which 95 display dynamic dichromatism) reveal high transition rates between dynamic dichromatism, ontogenetic (permanent) dichromatism and monochromatism reflecting the high evolutionary lability of this trait. Correlated evolution in hylids between dynamic dichromatism and forming large breeding aggregations indicates that the evolution of large breeding aggregations precedes the evolution of dynamic dichromatism. Multivariate phylogenetic logistic regression recovers the interaction between biogeographic distribution and forming breeding aggregations as a significant predictor of dynamic dichromatism in hylids. Accounting for macroecological differences between temperate and tropical regions, such as seasonality and the availability of breeding sites, may improve our understanding of ecological contexts in which dynamic dichromatism is likely to arise in tropical lineages and why it is retained in some temperate species and lost in others.

Biogeography of boreal passerine range dynamics in western North America: past, present, and future

Many of the Neotropical migrant bird species that breed throughout the Canadian boreal region are not found in the Alaskan boreal region, separated by the northwestern cordilleran mountains, despite the presence of climatically suitable habitat. We asked whether biological or climatic factors constrain certain species from crossing this geographic barrier. Analyzing a comprehensive dataset for 80 boreal passerine species, we used phylogenetic logistic regression to evaluate the relative importance of physical, migratory and competition metrics versus current and paleoclimatic suitability factors. Controlling for current climatic suitability within boreal Alaska, we found that species with the greatest climatic suitability across the northwestern cordillera, presently and also during the mid‐Holocene period, were most likely to be regular breeders in the Alaskan boreal region. Migratory strategy also played a role, but could not be disentangled from its strong phylogenetic basis. Our analysis suggests that the perceived barrier of the northwestern cordillera may be easily weakened as climate change improves conditions there for many forest species. The weakening of this barrier may lead to relatively rapid range expansions and the reshuffling of species communities. Species’ realized distributional shifts will be a function of the interplay between a changing climate and static topographic features.

Evolutionary correlation of heterodichogamy and floral longevity

Heterodichogamy, the coexistence of morphs differing in the time when the shift between male and female phases occurs, is a rare sexual dimorphism derived from synchronous dichogamy. In dichogamous populations, plants need to adjust the temporal sex role to deal with the negative frequency dependent selection imposed by population sex ratios. However, such temporal adjustment may be constrained in dichogamous species with bisexual flowers and short floral longevities (<2 days) and therefore the reciprocal floral strategy can easily invade the population, ultimately leading to heterodichogamy. We predict that heterodichogamy has evolved in species with short floral longevities. To test this prediction, we compiled data on floral longevity for 377 dichogamous and heterodichogamous species with bisexual flowers and ran a phylogenetic logistic regression between floral strategy (dichogamy/heterodichogamy) and floral longevity. The results showed that heterodichogamous species had significantly shorter floral longevities than dichogamous species. This result together with the fact that evolutionary transitions always occurred from dichogamy to heterodichogamy supports the conclusion of the existence of a significant and negative evolutionary correlation between the floral longevity and the probability to evolve heterodichogamy from dichogamy in bisexual species. We discuss that the rarity of heterodichogamy in nature may be the result of other selective pressures towards long floral longevity in dichogamous species.

Resource predictability and specialization in avian malaria parasites.

We tested the hypothesis that avian haemosporidian (malaria) parasites specialize on hosts that can be characterized as predictable resources at a site in Amazonian Ecuador. We incorporated host phylogenetic relationship and relative abundance in assessing parasite specialization, and we examined associations between parasite specialization and three host characteristics - abundance, mass and longevity - using quantile regression, phylogenetic logistic regression and t-tests. Hosts of specialist malaria parasite lineages were on average more abundant than hosts of generalist parasite lineages, but the relationship between host abundance and parasite specialization was not consistent across analyses. We also found support for a positive association between parasite specialization and host longevity, but this also was not consistent across analyses. Nonetheless, our findings suggest that the predictability of a host resource may play a role in the evolution of specialization. However, we also discuss two alternative explanations to the resource predictability hypothesis for specialization: (i) that interspecific interactions among the parasites themselves might constrain some parasites to a specialist strategy, and (ii) that frequent encounters with multiple host species, mediated by blood-sucking insects, might promote generalization within this system.

The ecological drivers of nuptial color evolution in darters (Percidae: Etheostomatinae).

Closely related animal lineages often vary in male coloration, and ecological selection is hypothesized to shape this variation. The role of ecological selection in inhibiting male color has been documented extensively at the population level, but relatively few studies have investigated the evolution of male coloration across a clade of closely related species. Darters are a diverse group of fishes that vary in the presence of elaborate male nuptial coloration, with some species exhibiting vivid color patterns and others mostly or entirely achromatic. We used phylogenetic logistic regression to test for correlations between the presence/absence of color traits across darter species and the ecological conditions in which these species occur. Environmental variables were correlated with the presence of nuptial color in darters with colorful species tending to inhabit environments that would support fewer predators and potentially transmit a broader spectrum of natural light compared to species lacking male coloration. We also tested the color preferences of a common darter predator, largemouth bass, and found that it exhibits a strong preference for red, providing further evidence of predation as a source of selection on color evolution in darters. Ecological selection therefore appears to be an important factor in dictating the presence or absence of male coloration in this group of fishes.

Evolutionary models for the retention of adult–adult social play in primates: The roles of diet and other factors associated with resource acquisition

What factors in animal life history facilitate or reduce the probability that a species will perform play behavior? While some relationships are known within species and across individuals, it is not obvious that such relationships can be used to explain differences and similarities in amount and type of play across large taxonomic groupings of animals, let alone transitions among them. Primates encompass a relatively large assemblage of species that differ in numerous dietary, habitat, reproductive, and physiological processes. While all juvenile primates engage in social play, far fewer primate species engage in social play as adults. Here, derived from theory and more small-scale comparisons, we explore several biological and behavioral phenomena that differ among nonhuman primates and which may explain differences in the occurrence of adult–adult social play. We used phylogenetic logistic regression to assess the correlation of adult play with various life-history, metabolic and socioecological variables. Although the main contribution of our paper is demonstrating use of phylogenetic methods in the context of play evolution, suggestive but not significant evidence was found that adult sexual and nonsexual social play is influenced by diet, habitat, reproductive, and metabolic factors, sometimes in opposite directions, that we discuss along with needed future analyses involving improved data and interactions among traits.