The host phylogeny determines viral infectivity and replication across Staphylococcus host species.

Virus host shifts, where a virus transmits to and infects a novel host species, are a major source of emerging infectious disease. Genetic similarity between eukaryotic host species has been shown to be an important determinant of the outcome of virus host shifts, but it is unclear if this is the case for prokaryotes where anti-virus defences can be transmitted by horizontal gene transfer and evolve rapidly. Here, we measure the susceptibility of 64 strains of Staphylococcaceae bacteria (48 strains of Staphylococcus aureus and 16 non-S. aureus species spanning 2 genera) to the bacteriophage ISP, which is currently under investigation for use in phage therapy. Using three methods-plaque assays, optical density (OD) assays, and quantitative (q)PCR-we find that the host phylogeny explains a large proportion of the variation in susceptibility to ISP across the host panel. These patterns were consistent in models of only S. aureus strains and models with a single representative from each Staphylococcaceae species, suggesting that these phylogenetic effects are conserved both within and among host species. We find positive correlations between susceptibility assessed using OD and qPCR and variable correlations between plaque assays and either OD or qPCR, suggesting that plaque assays alone may be inadequate to assess host range. Furthermore, we demonstrate that the phylogenetic relationships between bacterial hosts can generally be used to predict the susceptibility of bacterial strains to phage infection when the susceptibility of closely related hosts is known, although this approach produced large prediction errors in multiple strains where phylogeny was uninformative. Together, our results demonstrate the ability of bacterial host evolutionary relatedness to explain differences in susceptibility to phage infection, with implications for the development of ISP both as a phage therapy treatment and as an experimental system for the study of virus host shifts.

Directional selection and the evolution of breeding date in birds, revisited: Hard selection and the evolution of plasticity.

The mismatch between when individuals breed and when we think they should breed has been a long-standing problem in evolutionary ecology. Price et al. is a classic theory paper in this field and is mainly cited for its most obvious result: if individuals with high nutritional condition breed early, then the advantage of breeding early may be overestimated when information on nutritional condition is absent. Price at al.'s less obvious result is that individuals, on average, are expected to breed later than the optimum. Here, we provide an explanation of their non-intuitive result in terms of hard selection, and go on to show that neither of their results are expected to hold if the relationship between breeding date and nutrition is allowed to evolve. By introducing the assumption that the advantage of breeding early is greater for individuals in high nutritional condition, we show that their most cited result can be salvaged. However, individuals, on average, are expected to breed earlier than the optimum, not later. More generally, we also show that the hard selection mechanisms that underpin these results have major implications for the evolution of plasticity: when environmental heterogeneity becomes too great, plasticity is selected against, prohibiting the evolution of generalists.

Decomposing phenotypic skew and its effects on the predicted response to strong selection.

The major frameworks for predicting evolutionary change assume that a phenotype's underlying genetic and environmental components are normally distributed. However, the predictions of these frameworks may no longer hold if distributions are skewed. Despite this, phenotypic skew has never been decomposed, meaning the fundamental assumptions of quantitative genetics remain untested. Here we demonstrate that the substantial phenotypic skew in the body size of juvenile blue tits (Cyanistes caeruleus) is driven by environmental factors. Although skew had little impact on our predictions of selection response in this case, our results highlight the impact of skew on the estimation of inheritance and selection. Specifically, the nonlinear parent-offspring regressions induced by skew, alongside selective disappearance, can strongly bias estimates of heritability. The ubiquity of skew and strong directional selection on juvenile body size imply that heritability is commonly overestimated, which may in part explain the discrepancy between predicted and observed trait evolution.

Contemporary selection on MHC genes in a free-living ruminant population.

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate that contemporary natural selection is acting on MHC class II genes in Soay sheep and that the mode of selection on specific fitness components can be different mode from selection on total fitness.

The relative importance of plasticity versus genetic differentiation in explaining between population differences; a meta-analysis.

Both plasticity and genetic differentiation can contribute to phenotypic differences between populations. Using data on non-fitness traits from reciprocal transplant studies, we show that approximately 60% of traits exhibit co-gradient variation whereby genetic differences and plasticity-induced differences between populations are the same sign. In these cases, plasticity is about twice as important as genetic differentiation in explaining phenotypic divergence. In contrast to fitness traits, the amount of genotype by environment interaction is small. Of the 40% of traits that exhibit counter-gradient variation the majority seem to be hyperplastic whereby non-native individuals express phenotypes that exceed those of native individuals. In about 20% of cases plasticity causes non-native phenotypes to diverge from the native phenotype to a greater extent than if plasticity was absent, consistent with maladaptive plasticity. The degree to which genetic differentiation versus plasticity can explain phenotypic divergence varies a lot between species, but our proxies for motility and migration explain little of this variation.

Gradients in richness and turnover of a forest passerine's diet prior to breeding: A mixed model approach applied to faecal metabarcoding data.

Rather little is known about the dietary richness and variation of generalist insectivorous species, including birds, due primarily to difficulties in prey identification. Using faecal metabarcoding, we provide the most comprehensive analysis of a passerine's diet to date, identifying the relative magnitudes of biogeographic, habitat and temporal trends in the richness and turnover in diet of Cyanistes caeruleus (blue tit) along a 39 site and 2° latitudinal transect in Scotland. Faecal samples were collected in 2014-2015 from adult birds roosting in nestboxes prior to nest building. DNA was extracted from 793 samples and we amplified COI and 16S minibarcodes. We identified 432 molecular operational taxonomic units that correspond to putative dietary items. Most dietary items were rare, with Lepidoptera being the most abundant and taxon-rich prey order. Here, we present a statistical approach for estimation of gradients and intersample variation in taxonomic richness and turnover using a generalised linear mixed model. We discuss the merits of this approach over existing tools and present methods for model-based estimation of repeatability, taxon richness and Jaccard indices. We found that dietary richness increases significantly as spring advances, but changes little with elevation, latitude or local tree composition. In comparison, dietary composition exhibits significant turnover along temporal and spatial gradients and among sites. Our study shows the promise of faecal metabarcoding for inferring the macroecology of food webs, but we also highlight the challenge posed by contamination and make recommendations of laboratory and statistical practices to minimise its impact on inference.

The evolution of phenotypic plasticity when environments fluctuate in time and space.

Most theoretical studies have explored the evolution of plasticity when the environment, and therefore the optimal trait value, varies in time or space. When the environment varies in time and space, we show that genetic adaptation to Markovian temporal fluctuations depends on the between-generation autocorrelation in the environment in exactly the same way that genetic adaptation to spatial fluctuations depends on the probability of philopatry. This is because both measure the correlation in parent-offspring environments and therefore the effectiveness of a genetic response to selection. If the capacity to genetically respond to selection is stronger in one dimension (e.g., space), then plasticity mainly evolves in response to fluctuations in the other dimension (e.g., time). If the relationships between the environments of development and selection are the same in time and space, the evolved plastic response to temporal fluctuations is useful in a spatial context and genetic differentiation in space is reduced. However, if the relationships between the environments of development and selection are different, the optimal level of plasticity is different in the two dimensions. In this case, the plastic response that evolves to cope with temporal fluctuations may actually be maladaptive in space, resulting in the evolution of hyperplasticity or negative plasticity. These effects can be mitigated by spatial genetic differentiation that acts in opposition to plasticity resulting in counter-gradient variation. These results highlight the difficulty of making space-for-time substitutions in empirical work but identify the key parameters that need to be measured in order to test whether space-for-time substitutions are likely to be valid.

No evidence for sibling or parent-offspring coadaptation in a wild population of blue tits, despite high power.

Parent and offspring behaviors are expected to act as both the agents and targets of selection. This may generate parent-offspring coadaptation in which parent and offspring behaviors become genetically correlated in a way that increases inclusive fitness. Cross-fostering has been used to study parent-offspring coadaptation, with the prediction that offspring raised by non-relatives, or parents raising non-relatives, should suffer fitness costs. Using long-term data from more than 400 partially crossed broods of blue tits (Cyanistes caeruleus), we show that there is no difference in mass or survival between crossed and non-crossed chicks. However, previous studies for which the evidence for parent-offspring coadaptation is strongest compare chicks from fully crossed broods with those from non-crossed broods. When parent-offspring coadaptation acts at the level of the brood then partial cross-fostering experiments are not expected to show evidence of coadaptation. To test this, we performed an additional experiment (163 broods) in which clutches were either fully crossed, non-crossed, or partially crossed. In agreement with the long-term data, there was no evidence for parent-offspring coadaptation on offspring fitness despite high power. In addition there was no evidence of effects on parental fitness, nor evidence of sibling coadaptation, although the power of these tests was more modest.

Changes in temperature alter the potential outcomes of virus host shifts.

Host shifts-where a pathogen jumps between different host species-are an important source of emerging infectious disease. With on-going climate change there is an increasing need to understand the effect changes in temperature may have on emerging infectious disease. We investigated whether species' susceptibilities change with temperature and ask if susceptibility is greatest at different temperatures in different species. We infected 45 species of Drosophilidae with an RNA virus and measured how viral load changes with temperature. We found the host phylogeny explained a large proportion of the variation in viral load at each temperature, with strong phylogenetic correlations between viral loads across temperature. The variance in viral load increased with temperature, while the mean viral load did not. This suggests that as temperature increases the most susceptible species become more susceptible, and the least susceptible less so. We found no significant relationship between a species' susceptibility across temperatures, and proxies for thermal optima (critical thermal maximum and minimum or basal metabolic rate). These results suggest that whilst the rank order of species susceptibilities may remain the same with changes in temperature, some species may become more susceptible to a novel pathogen, and others less so.

RNA-Interference Pathways Display High Rates of Adaptive Protein Evolution in Multiple Invertebrates.

Conflict between organisms can lead to a reciprocal adaptation that manifests as an increased evolutionary rate in genes mediating the conflict. This adaptive signature has been observed in RNA-interference (RNAi) pathway genes involved in the suppression of viruses and transposable elements in Drosophila melanogaster, suggesting that a subset of Drosophila RNAi genes may be locked in an arms race with these parasites. However, it is not known whether rapid evolution of RNAi genes is a general phenomenon across invertebrates, or which RNAi genes generally evolve adaptively. Here we use population genomic data from eight invertebrate species to infer rates of adaptive sequence evolution, and to test for past and ongoing selective sweeps in RNAi genes. We assess rates of adaptive protein evolution across species using a formal meta-analytic framework to combine data across species and by implementing a multispecies generalized linear mixed model of mutation counts. Across species, we find that RNAi genes display a greater rate of adaptive protein substitution than other genes, and that this is primarily mediated by positive selection acting on the genes most likely to defend against viruses and transposable elements. In contrast, evidence for recent selective sweeps is broadly spread across functional classes of RNAi genes and differs substantially among species. Finally, we identify genes that exhibit elevated adaptive evolution across the analyzed insect species, perhaps due to concurrent parasite-mediated arms races.

Estimating the ability of plants to plastically track temperature-mediated shifts in the spring phenological optimum.

One consequence of rising spring temperatures is that the optimum timing of key life-history events may advance. Where this is the case, a population's fate may depend on the degree to which it is able to track a change in the optimum timing either via plasticity or via adaptation. Estimating the effect that temperature change will have on optimum timing using standard approaches is logistically challenging, with the result that very few estimates of this important parameter exist. Here we adopt an alternative statistical method that substitutes space for time to estimate the temperature sensitivity of the optimum timing of 22 plant species based on >200 000 spatiotemporal phenological observations from across the United Kingdom. We find that first leafing and flowering dates are sensitive to forcing (spring) temperatures, with optimum timing advancing by an average of 3 days °C-1 and plastic responses to forcing between -3 and -8 days °C-1 . Chilling (autumn/winter) temperatures and photoperiod tend to be important cues for species with early and late phenology, respectively. For most species, we find that plasticity is adaptive, and for seven species, plasticity is sufficient to track geographic variation in the optimum phenology. For four species, we find that plasticity is significantly steeper than the optimum slope that we estimate between forcing temperature and phenology, and we examine possible explanations for this countergradient pattern, including local adaptation.

Selection on parental performance opposes selection for larger body mass in a wild population of blue tits.

There is abundant evidence in many taxa for positive directional selection on body size, and yet little evidence for microevolutionary change. In many species, variation in body size is partly determined by the actions of parents, so a proposed explanation for stasis is the presence of a negative genetic correlation between direct and parental effects. Consequently, selecting genes for increased body size would result in a correlated decline in parental effects, reducing body size in the following generation. We show that these arguments implicitly assume that parental care is cost free, and that including a cost alters the predicted genetic architectures needed to explain stasis. Using a large cross-fostered population of blue tits, we estimate direct selection on parental effects for body mass, and show it is negative. Negative selection is consistent with a cost to parental care, mainly acting through a reduction in current fecundity rather than survival. Under these conditions, evolutionary stasis is possible for moderately negative genetic correlations between direct and parental effects. This is in contrast to the implausibly extreme correlations needed when care is assumed to be cost-free. Thus, we highlight the importance of accounting correctly for complete selection acting on traits across generations.

Transparency in Ecology and Evolution: Real Problems, Real Solutions.

To make progress scientists need to know what other researchers have found and how they found it. However, transparency is often insufficient across much of ecology and evolution. Researchers often fail to report results and methods in detail sufficient to permit interpretation and meta-analysis, and many results go entirely unreported. Further, these unreported results are often a biased subset. Thus the conclusions we can draw from the published literature are themselves often biased and sometimes might be entirely incorrect. Fortunately there is a movement across empirical disciplines, and now within ecology and evolution, to shape editorial policies to better promote transparency. This can be done by either requiring more disclosure by scientists or by developing incentives to encourage disclosure.

Passerines may be sufficiently plastic to track temperature-mediated shifts in optimum lay date.

Projecting the fates of populations under climate change is one of global change biology's foremost challenges. Here, we seek to identify the contributions that temperature-mediated local adaptation and plasticity make to spatial variation in nesting phenology, a phenotypic trait showing strong responses to warming. We apply a mixed modeling framework to a Britain-wide spatiotemporal dataset comprising >100 000 records of first egg dates from four single-brooded passerine bird species. The average temperature during a specific time period (sliding window) strongly predicts spatiotemporal variation in lay date. All four species exhibit phenological plasticity, advancing lay date by 2-5 days °C(-1) . The initiation of this sliding window is delayed further north, which may be a response to a photoperiod threshold. Using clinal trends in phenology and temperature, we are able to estimate the temperature sensitivity of selection on lay date (B), but our estimates are highly sensitive to the temporal position of the sliding window. If the sliding window is of fixed duration with a start date determined by photoperiod, we find B is tracked by phenotypic plasticity. If, instead, we allow the start and duration of the sliding window to change with latitude, we find plasticity does not track B, although in this case, at odds with theoretical expectations, our estimates of B differ across latitude vs. longitude. We argue that a model combining photoperiod and mean temperature is most consistent with current understanding of phenological cues in passerines, the results from which suggest that each species could respond to projected increases in spring temperatures through plasticity alone. However, our estimates of B require further validation.

The spatial scale of local adaptation in a stochastic environment.

The distribution of phenotypes in space will be a compromise between adaptive plasticity and local adaptation increasing the fit of phenotypes to local conditions and gene flow reducing that fit. Theoretical models on the evolution of quantitative characters on spatially explicit landscapes have only considered scenarios where optimum trait values change as deterministic functions of space. Here, these models are extended to include stochastic spatially autocorrelated aspects to the environment, and consequently the optimal phenotype. Under these conditions, the regression of phenotype on the environmental variable becomes steeper as the spatial scale on which populations are sampled becomes larger. Under certain deterministic models - such as linear clines - the regression is constant. The way in which the regression changes with spatial scale is informative about the degree of phenotypic plasticity, the relative scale of effective gene flow and the environmental dependency of selection. Connections to temporal models are discussed.

Heritability of heterozygosity offers a new way of understanding why dominant gene action contributes to additive genetic variance.

Whenever allele frequencies are unequal, nonadditive gene action contributes to additive genetic variance and therefore the resemblance between parents and offspring. The reason for this has not been easy to understand. Here, we present a new single-locus decomposition of additive genetic variance that may give greater intuition about this important result. We show that the contribution of dominant gene action to parent-offspring resemblance only depends on the degree to which the heterozygosity of parents and offspring covary. Thus, dominant gene action only contributes to additive genetic variance when heterozygosity is heritable. Under most circumstances this is the case because individuals with rare alleles are more likely to be heterozygous, and because they pass rare alleles to their offspring they also tend to have heterozygous offspring. When segregating alleles are at equal frequency there are no rare alleles, the heterozygosities of parents and offspring are uncorrelated and dominant gene action does not contribute to additive genetic variance.

Are molecular markers useful predictors of adaptive potential?

Estimates of molecular genetic variation are often used as a cheap and simple surrogate for a population's adaptive potential, yet empirical evidence suggests they are unlikely to be a valid proxy. However, this evidence is based on molecular genetic variation poorly predicting estimates of adaptive potential rather than how well it predicts true values. As a consequence, the relationship has been systematically underestimated and the precision with which it could be measured severely overstated. By collating a large database, and using suitable statistical methods, we obtain a 95% upper bound of 0.26 for the proportion of variance in quantitative genetic variation explained by molecular diversity. The relationship is probably too weak to be useful, but this conclusion must be taken as provisional: less noisy estimates of quantitative genetic variation are required. In contrast, and perhaps surprisingly, current sampling strategies appear sufficient for characterising a population's molecular genetic variation at comparable markers.