Controlling for non-independence in comparative analysis of patterns across populations within species

New horizons for comparative studies and meta-analyses

Comparative analyses and meta-analyses are key tools to elucidate broad biological principles, yet the two approaches often appear different in purpose. We propose an integrated approach that can generate deeper insights into ecoevolutionary processes. Marrying comparative and meta-analytic approaches will allow for (i) a more accurate investigation of drivers of biological variation, (ii) a greater ability to account for sources of non-independence in experimental data, (iii) more effective control of publication bias, and (iv) improved transparency and reproducibility. Stronger integration of meta-analytic and comparative studies can also broaden the scope from species-centric investigations to community-level responses and function-valued traits (e.g., reaction norms). We illuminate commonalities, differences, and the transformative potential of combining these methodologies for advancing ecology and evolutionary biology.

Parallel shifts in trout feeding morphology suggest rapid adaptation to alpine lake environments

Eco-evolutionary interactions following ecosystem change provide critical insight into the ability of organisms to adapt to shifting resource landscapes. Here we explore evidence for the rapid parallel evolution of trout feeding morphology following eco-evolutionary interactions with zooplankton in alpine lakes stocked at different points in time in the Wind River Range (Wyoming, USA). In this system, trout predation has altered the zooplankton species community and driven a decrease in average zooplankton size. In some lakes that were stocked decades ago, we find shifts in gill raker traits consistent with the hypothesis that trout have rapidly adapted to exploit available smaller-bodied zooplankton more effectively. We explore this morphological response in multiple lake populations across two species of trout (cutthroat trout, Oncorhynchus clarkii, and golden trout Oncorhynchus aguabonita) and examine the impact of resource availability on morphological variation in gill raker number among lakes. Furthermore, we present genetic data to provide evidence that historically stocked cutthroat trout populations likely derive from multiple population sources, and incorporate variation from genomic relatedness in our exploration of environmental predictors of feeding morphology. These findings describe rapid adaptation and eco-evolutionary interactions in trout and document an evolutionary response to novel, contemporary ecosystem change.

Variation in ontogenetic trajectories of limb dimensions in humans is attributable to both climatic effects and neutral evolution

Previous studies showed that there is variation in ontogenetic trajectories of human limb dimensions and proportions. However, little is known about the evolutionary significance of this variation. This study used a global sample of modern human immature long bone measurements and a multivariate linear mixed-effects model to study 1) whether the variation in ontogenetic trajectories of limb dimensions is consistent with ecogeographic predictions and 2) the effects of different evolutionary forces on the variation in ontogenetic trajectories. We found that genetic relatedness arising from neutral (nonselective) evolution, allometric variation associated with the change in size, and directional effects from climate all contributed to the variation in ontogenetic trajectories of all major long bone dimensions in modern humans. After accounting for the effects of neutral evolution and holding other effects considered in the current study constant, extreme temperatures have weak, positive associations with diaphyseal length and breadth measurements, while mean temperature shows negative associations with diaphyseal dimensions. The association with extreme temperatures fits the expectations of ecogeographic rules, while the association with mean temperature may explain the observed among-group variation in intralimb indices. The association with climate is present throughout ontogeny, suggesting an explanation of adaptation by natural selection as the most likely cause. On the other hand, genetic relatedness among groups, as structured by neutral evolutionary factors, is an important consideration when interpreting skeletal morphology, even for nonadult individuals.

Does genome size increase with water depth in marine fishes?

A growing body of research suggests that genome size in animals can be affected by ecological factors. Half a century ago, Ebeling et al. proposed that genome size increases with depth in some teleost fish groups and discussed a number of biological mechanisms that may explain this pattern (e.g., passive accumulation, adaptive acclimation). Using phylogenetic comparative approaches, we revisit this hypothesis based on genome size and ecological data from up to 708 marine fish species in combination with a set of large-scale phylogenies, including a newly inferred tree. We also conduct modeling approaches of trait evolution and implement a variety of regression analyses to assess the relationship between genome size and depth. Our reanalysis of Ebeling et al.'s dataset shows a weak association between these variables, but the overall pattern in their data is driven by a single clade. Although new analyses based on our "all-species" dataset resulted in positive correlations, providing some evidence that genome size evolves as a function of depth, only one subclade consistently yielded statistically significant correlations. By contrast, negative correlations are rare and nonsignificant. All in all, we find modest evidence for an increase in genome size along the depth axis in marine fishes. We discuss some mechanistic explanations for the observed trends.

Plasmids do not consistently stabilize cooperation across bacteria but may promote broad pathogen host-range

Horizontal gene transfer via plasmids could favour cooperation in bacteria, because transfer of a cooperative gene turns non-cooperative cheats into cooperators. This hypothesis has received support from theoretical, genomic and experimental analyses. In contrast, we show here, with a comparative analysis across 51 diverse species, that genes for extracellular proteins, which are likely to act as cooperative ‘public goods’, were not more likely to be carried on either: (i) plasmids compared to chromosomes; or (ii) plasmids that transfer at higher rates. Our results were supported by theoretical modelling which showed that while horizontal gene transfer can help cooperative genes initially invade a population, it has less influence on the longer-term maintenance of cooperation. Instead, we found that genes for extracellular proteins were more likely to be on plasmids when they coded for pathogenic virulence traits, in pathogenic bacteria with a broad host-range.

Ecological Interactions, Environmental Gradients, and Gene Flow in Local Adaptation

Despite long-standing interest in local adaptation of plants to their biotic and abiotic environment, existing theory, and many case studies, little work to date has addressed within-species evolution of concerted strategies and how these might contrast with patterns across species. Here we consider the interactions between pollinators, herbivores, and resource availability in shaping plant local adaptation, how these interactions impact plant phenotypes and gene flow, and the conditions where multiple traits align along major environmental gradients such as latitude and elevation. Continued work in emerging model systems will benefit from the melding of classic experimental approaches with novel population genetic analyses to reveal patterns and processes in plant local adaptation.

Temperature heterogeneity correlates with intraspecific variation in physiological flexibility in a small endotherm

Phenotypic flexibility allows individuals to reversibly modify trait values and theory predicts an individual's relative degree of flexibility positively correlates with the environmental heterogeneity it experiences. We test this prediction by integrating surveys of population genetic and physiological variation with thermal acclimation experiments and indices of environmental heterogeneity in the Dark-eyed Junco (Junco hyemalis) and its congeners. We combine field measures of thermogenic capacity for 335 individuals, 22,006 single nucleotide polymorphisms genotyped in 181 individuals, and laboratory acclimations replicated on five populations. We show that Junco populations: (1) differ in their thermogenic responses to temperature variation in the field; (2) harbor allelic variation that also correlates with temperature heterogeneity; and (3) exhibit intra-specific variation in thermogenic flexibility in the laboratory that correlates with the heterogeneity of their native thermal environment. These results provide comprehensive support that phenotypic flexibility corresponds with environmental heterogeneity and highlight its importance for coping with environmental change.

Replications, Comparisons, Sampling and the Problem of Representativeness in Animal Cognition Research

Animal cognition research often involves small and idiosyncratic samples. This can constrain the generalizability and replicability of a study's results and prevent meaningful comparisons between samples. However, there is little consensus about what makes a strong replication or comparison in animal research. We apply a resampling definition of replication to answer these questions in Part 1 of this article, and, in Part 2, we focus on the problem of representativeness in animal research. Through a case study and a simulation study, we highlight how and when representativeness may be an issue in animal behavior and cognition research and show how the representativeness problems can be viewed through the lenses of, i) replicability, ii) generalizability and external validity, iii) pseudoreplication and, iv) theory testing. Next, we discuss when and how researchers can improve their ability to learn from small sample research through, i) increasing heterogeneity in experimental design, ii) increasing homogeneity in experimental design, and, iii) statistically modeling variation. Finally, we describe how the strongest solutions will vary depending on the goals and resources of individual research programs and discuss some barriers towards implementing them.

Free-Living, Psychrotrophic Bacteria of the Genus Psychrobacter Are Descendants of Pathobionts

Host-adapted microorganisms are generally assumed to have evolved from free-living, environmental microorganisms, as examples of the reverse process are rare. In the phylum Gammaproteobacteria, family Moraxellaceae, the genus Psychrobacter includes strains from a broad ecological distribution including animal bodies as well as sea ice and other nonhost environments. To elucidate the relationship between these ecological niches and Psychrobacter's evolutionary history, we performed tandem genomic analyses with phenotyping of 85 Psychrobacter accessions. Phylogenomic analysis of the family Moraxellaceae reveals that basal members of the Psychrobacter clade are Moraxella spp., a group of often-pathogenic organisms. Psychrobacter exhibited two broad growth patterns in our phenotypic screen: one group that we called the "flexible ecotype" (FE) had the ability to grow between 4 and 37°C, and the other, which we called the "restricted ecotype" (RE), could grow between 4 and 25°C. The FE group includes phylogenetically basal strains, and FE strains exhibit increased transposon copy numbers, smaller genomes, and a higher likelihood to be bile salt resistant. The RE group contains only phylogenetically derived strains and has increased proportions of lipid metabolism and biofilm formation genes, functions that are adaptive to cold stress. In a 16S rRNA gene survey of polar bear fecal samples, we detect both FE and RE strains, but in in vivo colonizations of gnotobiotic mice, only FE strains persist. Our results indicate the ability to grow at 37°C, seemingly necessary for mammalian gut colonization, is an ancestral trait for Psychrobacter, which likely evolved from a pathobiont.IMPORTANCE Host-associated microbes are generally assumed to have evolved from free-living ones. The evolutionary transition of microbes in the opposite direction, from host associated toward free living, has been predicted based on phylogenetic data but not studied in depth. Here, we provide evidence that the genus Psychrobacter, particularly well known for inhabiting low-temperature, high-salt environments such as sea ice, permafrost soils, and frozen foodstuffs, has evolved from a mammalian-associated ancestor. We show that some Psychrobacter strains retain seemingly ancestral genomic and phenotypic traits that correspond with host association while others have diverged to psychrotrophic or psychrophilic lifestyles.

Phylogenetic Permulations: A Statistically Rigorous Approach to Measure Confidence in Associations in a Phylogenetic Context

Many evolutionary comparative methods seek to identify associations between phenotypic traits or between traits and genotypes, often with the goal of inferring potential functional relationships between them. Comparative genomics methods aimed at this goal measure the association between evolutionary changes at the genetic level with traits evolving convergently across phylogenetic lineages. However, these methods have complex statistical behaviors that are influenced by nontrivial and oftentimes unknown confounding factors. Consequently, using standard statistical analyses in interpreting the outputs of these methods leads to potentially inaccurate conclusions. Here, we introduce phylogenetic permulations, a novel statistical strategy that combines phylogenetic simulations and permutations to calculate accurate, unbiased P values from phylogenetic methods. Permulations construct the null expectation for P values from a given phylogenetic method by empirically generating null phenotypes. Subsequently, empirical P values that capture the true statistical confidence given the correlation structure in the data are directly calculated based on the empirical null expectation. We examine the performance of permulation methods by analyzing both binary and continuous phenotypes, including marine, subterranean, and long-lived large-bodied mammal phenotypes. Our results reveal that permulations improve the statistical power of phylogenetic analyses and correctly calibrate statements of confidence in rejecting complex null distributions while maintaining or improving the enrichment of known functions related to the phenotype. We also find that permulations refine pathway enrichment analyses by correcting for nonindependence in gene ranks. Our results demonstrate that permulations are a powerful tool for improving statistical confidence in the conclusions of phylogenetic analysis when the parametric null is unknown.

No impact of a short-term climatic "El Niño" fluctuation on gut microbial diversity in populations of the Galápagos marine iguana (Amblyrhynchus cristatus)

Gut microorganisms are crucial for many biological functions playing a pivotal role in the host's well-being. We studied gut bacterial community structure of marine iguana populations across the Galápagos archipelago. Marine iguanas depend heavily on their specialized gut microbiome for the digestion of dietary algae, a resource whose growth was strongly reduced by severe "El Niño"-related climatic fluctuations in 2015/2016. As a consequence, marine iguana populations showed signs of starvation as expressed by a poor body condition. Body condition indices (BCI) varied between island populations indicating that food resources (i.e., algae) are affected differently across the archipelago during 'El Niño' events. Though this event impacted food availability for marine iguanas, we found that reductions in body condition due to "El Niño"-related starvation did not result in differences in bacterial gut community structure. Species richness of gut microorganisms was instead correlated with levels of neutral genetic diversity in the distinct host populations. Our data suggest that marine iguana populations with a higher level of gene diversity and allelic richness may harbor a more diverse gut microbiome than those populations with lower genetic diversity. Since low values of these diversity parameters usually correlate with small census and effective population sizes, we use our results to propose a novel hypothesis according to which small and genetically less diverse host populations might be characterized by less diverse microbiomes. Whether such genetically depauperate populations may experience additional threats from reduced dietary flexibility due to a limited intestinal microbiome is currently unclear and calls for further investigation.

The evolution of sex chromosome dosage compensation in animals

The evolution of heteromorphic sex chromosomes shall lead to gene expression dosage problems, as in at least one of the sexes, the sex-linked gene dose has been reduced by half. It has been proposed that the transcriptional output of the whole X or Z chromosome should be doubled for complete dosage compensation in heterogametic sex. However, owing to the variability of the existing methods to determine the transcriptional differences between sex chromosomes and autosomes (S:A ratios) in different studies, we collected more than 500 public RNA-Seq data set from multiple tissues and species in major clades and proposed a unified computational framework for unbiased and comparable measurement of the S:A ratios of multiple species. We also tested the evolution of dosage compensation more directly by assessing changes in the expression levels of the current sex-linked genes relative to those of the ancestral sex-linked genes. We found that in mammals and birds, the S:A ratio is approximately 0.5, whereas in insects, fishes, and flatworms, the S:A ratio is approximately 1.0. Further analysis showed that the fraction of dosage-sensitive housekeeping genes on the X/Z chromosome is significantly correlated with the S:A ratio. In addition, the degree of degeneration of the Y chromosome may be responsible for the change in the S:A ratio in mammals without a dosage compensation mechanism. Our observations offer unequivocal support for the sex chromosome insensitivity hypothesis in animals and suggest that dosage sensitivity states of sex chromosomes are a major factor underlying different evolutionary strategies of dosage compensation.

The link between selection for function and human-directed play behaviour in dogs

Human-directed play behaviour is a distinct behavioural feature of domestic dogs. But the role that artificial selection for contemporary dog breeds has played for human-directed play behaviour remains elusive. Here, we investigate how human-directed play behaviour has evolved in relation to the selection for different functions, considering processes of shared ancestry and gene flow among the different breeds. We use the American Kennel Club (AKC) breed group categorization to reflect the major functional differences and combine this with observational data on human-directed play behaviour for over 132 breeds across 89 352 individuals from the Swedish Dog Mentality Assessment project. Our analyses demonstrate that ancestor dogs already showed intermediate levels of human-directed play behaviour, levels that are shared with several modern breed types. Herding and Sporting breeds display higher levels of human-directed play behaviour, statistically distinguishable from Non-sporting and Toy breeds. Our results suggest that human-directed play behaviour played a role in the early domestication of dogs and that subsequent artificial selection for function has been important for contemporary variation in a behavioural phenotype mediating the social bond with humans.

Monotocy and the evolution of plural breeding in mammals

In many mammals, breeding females are intolerant of each other and seldom associate closely but, in some, they aggregate in groups that vary in size, stability, and kinship structure. Aggregation frequently increases competition for food, and interspecific differences in female sociality among mammals are commonly attributed to contrasts in ecological parameters, including variation in activity timing, the distribution of resources, as well as the risk of predation. However, there is increasing indication that differences in female sociality are also associated with phylogenetic relationships and with contrasts in life-history parameters. We show here that evolutionary transitions from systems where breeding females usually occupy separate ranges ("singular breeding") to systems where breeding females usually aggregate ("plural breeding") have occurred more frequently in monotocous lineages where females produce single young than in polytocous ones where they produce litters. A likely explanation of this association is that competition between breeding females for resources is reduced where they produce single young and is more intense where they produce litters. Our findings reinforce evidence that variation in life-history parameters plays an important role in shaping the evolution of social behavior.

Modern morphometrics and the study of population differences: Good data behind clever analyses and cool pictures?

The study of phenotypic variation in time and space is central to evolutionary biology. Modern geometric morphometrics is the leading family of methods for the quantitative analysis of biological forms. This set of techniques relies heavily on technological innovation for data acquisition, often in the form of 2D or 3D digital images, and on powerful multivariate statistical tools for their analysis. However, neither the most sophisticated device for computerized imaging nor the best statistical test can produce accurate, robust and reproducible results, if it is not based on really good samples and an appropriate use of the 'measurements' extracted from the data. Using examples mostly from my own work on mammal craniofacial variation and museum specimens, I will show how easy it is to forget these most basic assumptions, while focusing heavily on analytical and visualization methods, and much less on the data that generate potentially powerful analyses and visually appealing diagrams.

Studying convergent evolution to relate genotype to behavioral phenotype

Neuroscience has a long, rich history in embracing unusual animals for research. Over the past several decades, there has been a technology-driven bottleneck in the species used for neuroscience research. However, an oncoming wave of technologies applicable to many animals hold promise for enabling researchers to address challenging scientific questions that cannot be solved using traditional laboratory animals. Here, we discuss how leveraging the convergent evolution of physiological or behavioral phenotypes can empower research mapping genotype to phenotype interactions. We present two case studies using electric fish and poison frogs and discuss how comparative work can teach us about evolutionary constraint and flexibility at various levels of biological organization. We also offer advice on the potential and pitfalls of establishing novel model systems in neuroscience research. Finally, we end with a discussion on the use of charismatic animals in neuroscience research and their utility in public outreach. Overall, we argue that convergent evolution frameworks can help identify generalizable principles of neuroscience.

The role of common ancestry and gene flow in the evolution of human-directed play behaviour in dogs

Among-population variance of phenotypic traits is of high relevance for understanding evolutionary mechanisms that operate in relatively short timescales, but various sources of nonindependence, such as common ancestry and gene flow, can hamper the interpretations. In this comparative analysis of 138 dog breeds, we demonstrate how such confounders can independently shape the evolution of a behavioural trait (human-directed play behaviour from the Dog Mentality Assessment project). We combined information on genetic relatedness and haplotype sharing to reflect common ancestry and gene flow, respectively, and entered these into a phylogenetic mixed model to partition the among-breed variance of human-directed play behaviour while also accounting for within-breed variance. We found that 75% of the among-breed variance was explained by overall genetic relatedness among breeds, whereas 15% could be attributed to haplotype sharing that arises from gene flow. Therefore, most of the differences in human-directed play behaviour among breeds have likely been caused by constraints of common ancestry as a likely consequence of past selection regimes. On the other hand, gene flow caused by crosses among breeds has played a minor, but not negligible role. Our study serves as an example of an analytical approach that can be applied to comparative situations where the effects of shared origin and gene flow require quantification and appropriate statistical control in a within-species/among-population framework. Altogether, our results suggest that the evolutionary history of dog breeds has left remarkable signatures on the among-breed variation of a behavioural phenotype.

Population genomics perspectives on convergent adaptation

Convergent adaptation is the independent evolution of similar traits conferring a fitness advantage in two or more lineages. Cases of convergent adaptation inform our ideas about the ecological and molecular basis of adaptation. In judging the degree to which putative cases of convergent adaptation provide an independent replication of the process of adaptation, it is necessary to establish the degree to which the evolutionary change is unexpected under null models and to show that selection has repeatedly, independently driven these changes. Here, we discuss the issues that arise from these questions particularly for closely related populations, where gene flow and standing variation add additional layers of complexity. We outline a conceptual framework to guide intuition as to the extent to which evolutionary change represents the independent gain of information owing to selection and show that this is a measure of how surprised we should be by convergence. Additionally, we summarize the ways population and quantitative genetics and genomics may help us address questions related to convergent adaptation, as well as open new questions and avenues of research. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Evolutionary potential varies across populations and traits in the neotropical oak Quercus oleoides

Heritable variation in polygenic (quantitative) traits is critical for adaptive evolution and is especially important in this era of rapid climate change. In this study, we examined the levels of quantitative genetic variation of populations of the tropical tree Quercus oleoides Cham. and Schlect. for a suite of traits related to resource use and drought resistance. We tested whether quantitative genetic variation differed across traits, populations and watering treatments. We also tested potential evolutionary factors that might have shaped such a pattern: selection by climate and genetic drift. We measured 15 functional traits on 1322 1-year-old seedlings of 84 maternal half-sib families originating from five populations growing under two watering treatments in a greenhouse. We estimated the additive genetic variance, coefficient of additive genetic variation and narrow-sense heritability for each combination of traits, populations and treatments. In addition, we genotyped a total of 119 individuals (with at least 20 individuals per population) using nuclear microsatellites to estimate genetic diversity and population genetic structure. Our results showed that gas exchange traits and growth exhibited strikingly high quantitative genetic variation compared with traits related to leaf morphology, anatomy and photochemistry. Quantitative genetic variation differed between populations even at geographical scales as small as a few kilometers. Climate was associated with quantitative genetic variation, but only weakly. Genetic structure and diversity in neutral markers did not relate to coefficient of additive genetic variation. Our study demonstrates that quantitative genetic variation is not homogeneous across traits and populations of Q. oleoides. More importantly, our findings suggest that predictions about potential responses of species to climate change need to consider population-specific evolutionary characteristics.

Rapid phenotypic evolution following shifts in life cycle complexity

Life cycle strategies have evolved extensively throughout the history of metazoans. The expression of disparate life stages within a single ontogeny can present conflicts to trait evolution, and therefore may have played a major role in shaping metazoan forms. However, few studies have examined the consequences of adding or subtracting life stages on patterns of trait evolution. By analysing trait evolution in a clade of closely related salamander lineages we show that shifts in the number of life cycle stages are associated with rapid phenotypic evolution. Specifically, salamanders with an aquatic-only (paedomorphic) life cycle have frequently added vertebrae to their trunk skeleton compared with closely related lineages with a complex aquatic-to-terrestrial (biphasic) life cycle. The rate of vertebral column evolution is also substantially lower in biphasic lineages, which may reflect the functional compromise of a complex cycle. This study demonstrates that the consequences of life cycle evolution can be detected at very fine scales of divergence. Rapid evolutionary responses can result from shifts in selective regimes following changes in life cycle complexity.

Intraspecific variation in lizard heat tolerance alters estimates of climate impact

Research addressing the effects of global warming on the distribution and persistence of species generally assumes that population variation in thermal tolerance is spatially constant or overridden by interspecific variation. Typically, this rationale is implicit in sourcing one critical thermal maximum (CT_max ) population estimate per species to model spatiotemporal cross-taxa variation in heat tolerance. Theory suggests that such an approach could result in biased or imprecise estimates and forecasts of impact from climate warming, but limited empirical evidence in support of those expectations exists. We experimentally quantify the magnitude of intraspecific variation in CT_max among lizard populations, and the extent to which incorporating such variability can alter estimates of climate impact through a biophysical model. To do so, we measured CT_max from 59 populations of 15 Iberian lizard species (304 individuals). The overall median CT_max across all individuals from all species was 42.8°C and ranged from 40.5 to 48.3°C, with species medians decreasing through xeric, climate-generalist and mesic taxa. We found strong statistical support for intraspecific differentiation in CT_max by up to a median of 3°C among populations. We show that annual restricted activity (operative temperature > CT_max ) over the Iberian distribution of our study species differs by a median of >80 hr per 25-km² grid cell based on different population-level CT_max estimates. This discrepancy leads to predictions of spatial variation in annual restricted activity to change by more than 20 days for six of the study species. Considering that during restriction periods, reptiles should be unable to feed and reproduce, current projections of climate-change impacts on the fitness of ectotherm fauna could be under- or over-estimated depending on which population is chosen to represent the physiological spectra of the species in question. Mapping heat tolerance over the full geographical ranges of single species is thus critical to address cross-taxa patterns and drivers of heat tolerance in a biologically comprehensive way.

Feeding ecology and sexual dimorphism in a speciose flower beetle clade (Hopliini: Scarabaeidae)

The relationship between feeding ecology and sexual dimorphism is examined in a speciose South African monkey beetle clade. We test whether feeding and mating at a fixed site (embedding guild) is associated with greater levels of sexual dimorphism and possibly sexual selection than species using unpredictable feeding resources (non-embedding guild). Sexual dimorphism was measured using a point scoring system for hind leg and colour across the two feeding guilds for >50% of the regional fauna. Quantification of hind leg dimorphism using a scoring system and allometric scaling were used to identify traits subject to sexual selection. Feeding guild had a significant effect on hind leg dimorphism, with embedders having high and non-embedders low scores. The sessile and defendable distribution of females on stable platform flowers may favour contests and associated hind leg weaponry. In contrast, degree of colour dimorphism between the sexes was not associated with any particular feeding guild, and may serve to reduce male conflict and combat. Embedder males had high proportions (∼76%) of species with positive allometric slopes for almost all hind leg traits. For male non-embedders, only ∼37% of species showed positive scaling relationships. Phylogenetic data, in conjunction with behavioural data on the function of leg weaponry and visual signalling among males is needed to better understand the link between sexual dimorphism and sexual selection in the radiation of the monkey beetles.

Resilience to fire of phylogenetic diversity across biological domains

Fire alters the structure and composition of above- and belowground communities with concurrent shifts in phylogenetic diversity. The inspection of postfire trends in the diversity of ecological communities incorporating phylogenetic information allows to better understand the mechanisms driving fire resilience. While fire reduces plant phylogenetic diversity based on the recruitment of evolutionarily related species with postfire seed persistence, it increases that of soil microbes by limiting soil resources and changing the dominance of competing microbes. Thus, during postfire community reassembly, plant and soil microbes might experience opposing temporal trends in their phylogenetic diversity that are linked through changes in the soil conditions. We tested this hypothesis by investigating the postfire evolution of plant and soil microbial (fungi, bacteria and archaea) communities across three 20-year chronosequences. Plant phylogenetic diversity increased with time since fire as pioneer seeders facilitate the establishment of distantly related late-successional shrubs. The postfire increase in plant phylogenetic diversity fostered plant productivity, eventually recovering soil organic matter. These shifts over time in the soil conditions explained the postfire restoration of fungal and bacterial phylogenetic diversity, which decreased to prefire levels, suggesting that evolutionarily related taxa with high relative fitness recover their competitive superiority during community reassembly. The resilience to fire of phylogenetic diversity across biological domains helps preserve the evolutionary history stored in our ecosystems.

Improving generation length estimates for the IUCN Red List

The International Union for the Conservation of Nature (IUCN) Red List classifies species according to their risk of extinction, informing local to global conservation decisions. Here we look to advance the estimation of generation length, which is used as a time-scalar in the Red List as a way of accounting for differences in species’ life-histories. We calculated or predicted generation length for 86 species of antelope following the Rspan approach. We also tested the importance of both allometry (body-mass) and phylogeny (phylogenetic eigenvectors) as predictors of generation length within a Phylogenetic Eigenvector Map (PEM) framework. We then evaluated the predictive power of this PEM and two binning approaches, following a leave-one-out cross-validation routine. We showed that captive and wild longevity data are nonequivalent and that both body-mass and phylogeny are important predictors for generation length (body-mass explained 64% and phylogeny 36% of the partitioned explained variance). Plus, both the PEM, and the binning approach that included both taxonomic rank and body-mass, had good predictive power and therefore are suitable for extrapolating generation length to missing-data species. Therefore, based on our findings, we advise separating captive and wild data when estimating generation length, and considering the implications of wild and captive data more widely in life-history analyses. We also recommend that body-mass and phylogeny should be used in combination, preferably under a PEM framework (as it was less reliant on available reference species and more explicitly accounts for phylogenetic relatedness) or a binning approach if a PEM is not feasible, to extrapolate generation length to missing-data species. Overall, we provide a transparent, consistent and transferable workflow for improving the use of the Rspan method to calculate generation length for the IUCN Red List.

Using the Ornstein-Uhlenbeck process to model the evolution of interacting populations

The Ornstein-Uhlenbeck (OU) process plays a major role in the analysis of the evolution of phenotypic traits along phylogenies. The standard OU process includes random perturbations and stabilizing selection and assumes that species evolve independently. However, evolving species may interact through various ecological process and also exchange genes especially in plants. This is particularly true if we want to study phenotypic evolution among diverging populations within species. In this work we present a straightforward statistical approach with analytical solutions that allows for the inclusion of adaptation and migration in a common phylogenetic framework, which can also be useful for studying local adaptation among populations within the same species. We furthermore present a detailed simulation study that clearly indicates the adverse effects of ignoring migration. Similarity between species due to migration could be misinterpreted as very strong convergent evolution without proper correction for these additional dependencies. Finally, we show that our model can be interpreted in terms of ecological interactions between species, providing a general framework for the evolution of traits between "interacting" species or populations.

Factors contributing to the accumulation of reproductive isolation: A mixed model approach

The analysis of large datasets describing reproductive isolation between species has been extremely influential in the study of speciation. However, the statistical methods currently used for these data limit the ability to make direct inferences about the factors predicting the evolution of reproductive isolation. As a result, our understanding of iconic patterns and rules of speciation rely on indirect analyses that have clear statistical limitations. Phylogenetic mixed models are commonly used in ecology and evolution, but have not been applied to studies of reproductive isolation. Here I describe a flexible framework using phylogenetic mixed models to analyze data collected at different evolutionary scales, to test both categorical and continuous predictor variables, and to test the effect of multiple predictors on rates and patterns of reproductive isolation simultaneously. I demonstrate the utility of this framework by re-analyzing four classic datasets, from both animals and plants, and evaluating several hypotheses that could not be tested in the original studies: In the Drosophila and Bufonidae datasets, I found support for more rapid accumulation of reproductive isolation in sympatric species pairs compared to allopatric species pairs. Using Silene and Nolana, I found no evidence supporting the hypothesis that floral differentiation elevates postzygotic reproductive isolation. The faster accumulation of postzygotic isolation in sympatry is likely the result of species coexistence determined by the level of postzygotic isolation between species. In addition, floral trait divergence does not appear to translate into pleiotropic effects on postzygotic reproductive isolation. Overall, these methods can allow researchers to test new hypotheses using a single statistical method, while remedying the statistical limitations of several previous methods.

Temporal shifts and temperature sensitivity of avian spring migratory phenology: a phylogenetic meta-analysis

There are wide reports of advances in the timing of spring migration of birds over time and in relation to rising temperatures, though phenological responses vary substantially within and among species. An understanding of the ecological, life-history and geographic variables that predict this intra- and interspecific variation can guide our projections of how populations and species are likely to respond to future climate change. Here, we conduct phylogenetic meta-analyses addressing slope estimates of the timing of avian spring migration regressed on (i) year and (ii) temperature, representing a total of 413 species across five continents. We take into account slope estimation error and examine phylogenetic, ecological and geographic predictors of intra- and interspecific variation. We confirm earlier findings that on average birds have significantly advanced their spring migration time by 2·1 days per decade and 1·2 days °C-1 . We find that over time and in response to warmer spring conditions, short-distance migrants have advanced spring migratory phenology by more than long-distance migrants. We also find that larger bodied species show greater advance over time compared to smaller bodied species. Our results did not reveal any evidence that interspecific variation in migration response is predictable on the basis of species' habitat or diet. We detected a substantial phylogenetic signal in migration time in response to both year and temperature, suggesting that some of the shifts in migratory phenological response to climate are predictable on the basis of phylogeny. However, we estimate high levels of species and spatial variance relative to phylogenetic variance, which is consistent with plasticity in response to climate evolving fairly rapidly and being more influenced by adaptation to current local climate than by common descent. On average, avian spring migration times have advanced over time and as spring has become warmer. While we are able to identify predictors that explain some of the true among-species variation in response, substantial intra- and interspecific variation in migratory response remains to be explained.

The ultimate and proximate mechanisms driving the evolution of long tails in forest deer mice

Understanding both the role of selection in driving phenotypic change and its underlying genetic basis remain major challenges in evolutionary biology. Here, we use modern tools to revisit a classic system of local adaptation in the North American deer mouse, Peromyscus maniculatus, which occupies two main habitat types: prairie and forest. Using historical collections, we find that forest‐dwelling mice have longer tails than those from nonforested habitat, even when we account for individual and population relatedness. Using genome‐wide SNP data, we show that mice from forested habitats in the eastern and western parts of their range form separate clades, suggesting that increased tail length evolved independently. We find that forest mice in the east and west have both more and longer caudal vertebrae, but not trunk vertebrae, than nearby prairie forms. By intercrossing prairie and forest mice, we show that the number and length of caudal vertebrae are not correlated in this recombinant population, indicating that variation in these traits is controlled by separate genetic loci. Together, these results demonstrate convergent evolution of the long‐tailed forest phenotype through two distinct genetic mechanisms, affecting number and length of vertebrae, and suggest that these morphological changes—either independently or together—are adaptive.

Navigating the currents of seascape genomics: how spatial analyses can augment population genomic studies

Population genomic approaches are making rapid inroads in the study of non-model organisms, including marine taxa. To date, these marine studies have predominantly focused on rudimentary metrics describing the spatial and environmental context of their study region (e.g., geographical distance, average sea surface temperature, average salinity). We contend that a more nuanced and considered approach to quantifying seascape dynamics and patterns can strengthen population genomic investigations and help identify spatial, temporal, and environmental factors associated with differing selective regimes or demographic histories. Nevertheless, approaches for quantifying marine landscapes are complicated. Characteristic features of the marine environment, including pelagic living in flowing water (experienced by most marine taxa at some point in their life cycle), require a well-designed spatial-temporal sampling strategy and analysis. Many genetic summary statistics used to describe populations may be inappropriate for marine species with large population sizes, large species ranges, stochastic recruitment, and asymmetrical gene flow. Finally, statistical approaches for testing associations between seascapes and population genomic patterns are still maturing with no single approach able to capture all relevant considerations. None of these issues are completely unique to marine systems and therefore similar issues and solutions will be shared for many organisms regardless of habitat. Here, we outline goals and spatial approaches for landscape genomics with an emphasis on marine systems and review the growing empirical literature on seascape genomics. We review established tools and approaches and highlight promising new strategies to overcome select issues including a strategy to spatially optimize sampling. Despite the many challenges, we argue that marine systems may be especially well suited for identifying candidate genomic regions under environmentally mediated selection and that seascape genomic approaches are especially useful for identifying robust locus-by-environment associations.

Animal choruses emerge from receiver psychology

Synchrony and alternation in large animal choruses are often viewed as adaptations by which cooperating males increase their attractiveness to females or evade predators. Alternatively, these seemingly composed productions may simply emerge by default from the receiver psychology of mate choice. This second, emergent property hypothesis has been inferred from findings that females in various acoustic species ignore male calls that follow a neighbor's by a brief interval, that males often adjust the timing of their call rhythm and reduce the incidence of ineffective, following calls, and from simulations modeling the collective outcome of male adjustments. However, the purported connection between male song timing and female preference has never been tested experimentally, and the emergent property hypothesis has remained speculative. Studying a distinctive katydid species genetically structured as isolated populations, we conducted a comparative phylogenetic analysis of the correlation between male call timing and female preference. We report that across 17 sampled populations male adjustments match the interval over which females prefer leading calls; moreover, this correlation holds after correction for phylogenetic signal. Our study is the first demonstration that male adjustments coevolved with female preferences and thereby confirms the critical link in the emergent property model of chorus evolution.

Macroscale intraspecific variation and environmental heterogeneity: analysis of cold and warm zone abundance, mortality, and regeneration distributions of four eastern US tree species

I test for macroscale intraspecific variation of abundance, mortality, and regeneration of four eastern US tree species (Tsuga canadensis,Betula lenta,Liriodendron tulipifera, and Quercus prinus) by splitting them into three climatic zones based on plant hardiness zones (PHZs). The primary goals of the analysis are to assess the differences in environmental heterogeneity and demographic responses among climatic zones, map regional species groups based on decision tree rules, and evaluate univariate and multivariate patterns of species demography with respect to environmental variables. I use the Forest Inventory Analysis (FIA) data to derive abundance, mortality, and regeneration indices and split the range into three climatic zones based on USDA PHZs: (1) cold adapted, leading region; (2) middle, well‐adapted region; and (3) warm adapted, trailing region. I employ decision tree ensemble methods to assess the importance of environmental predictors on the abundance of the species between the cold and warm zones and map zonal variations in species groups. Multivariate regression trees are used to simultaneously explore abundance, mortality, and regeneration in tandem to assess species vulnerability. Analyses point to the relative importance of climate in the warm adapted, trailing zone (especially moisture) compared to the cold adapted, leading zone. Higher mortality and lower regeneration patterns in the warm trailing zone point to its vulnerability to growing season temperature and precipitation changes that could figure more prominently in the future. This study highlights the need to account for intraspecific variation of demography in order to understand environmental heterogeneity and differential adaptation. It provides a methodology for assessing the vulnerability of tree species by delineating climatic zones based on easily available PHZ data, and FIA derived abundance, mortality, and regeneration indices as a proxy for overall growth and fitness. Based on decision tree rules, ecologically meaningful variations in species abundance among the climatic zones can be related to environmental variability and mapped.

Impacts of local adaptation of forest trees on associations with herbivorous insects: implications for adaptive forest management

Disruption of species interactions is a key issue in climate change biology. Interactions involving forest trees may be particularly vulnerable due to evolutionary rate limitations imposed by long generation times. One mitigation strategy for such impacts is Climate matching - the augmentation of local native tree populations by input from nonlocal populations currently experiencing predicted future climates. This strategy is controversial because of potential cascading impacts on locally adapted animal communities. We explored these impacts using abundance data for local native gallwasp herbivores sampled from 20 provenances of sessile oak (Quercus petraea) planted in a common garden trial. We hypothesized that non-native provenances would show (i) declining growth performance with increasing distance between provenance origin and trial site, and (ii) phenological differences to local oaks that increased with latitudinal differences between origin and trial site. Under a local adaptation hypothesis, we predicted declining gallwasp abundance with increasing phenological mismatch between native and climate-matched trees. Both hypotheses for oaks were supported. Provenance explained significant variation in gallwasp abundance, but no gall type showed the relationship between abundance and phenological mismatch predicted by a local adaptation hypothesis. Our results show that climate matching would have complex and variable impacts on oak gall communities.