Unidirectional response to bidirectional selection on body size. I. Phenotypic, life-history, and endocrine responses

Evolutionary implications of size-selective mortality on the ontogenetic development of shoal cohesion: a neurochemical approach using a zebrafish, Danio rerio, harvest selection experiment

Abstract

Size-selective mortality may evolutionarily alter life-history as well as individual behavioral and physiological traits. Moreover, size-selective mortality can affect group behavioral traits, such as shoaling and collective properties (e.g., shoal cohesion), which are relevant for finding food and reducing risk of predation. Here, we present experimental evidence using selection lines of zebrafish (Danio rerio) that were exposed to positive (large-harvested), negative (small-harvested), and random (control) size-selective mortality for five generations, followed by eight generations during which harvesting was halted to remove maternal effects and to study evolutionarily fixed outcomes. We investigated changes in shoal cohesion and turnover in monoamines in zebrafish through ontogeny. To that end, we repeatedly measured inter-individual distance in groups of eight fish and the turnovers of dopamine and serotonin in brains of fish from juvenile to the adult stage at 40-day intervals. We, firstly, found that shoal cohesion was overall consistent through ontogeny at group levels suggesting the presence of collective personality. Secondly, we found a decrease in shoal cohesion through ontogeny in the small-harvested and control lines, while the large-harvested line did not show any ontogenetic change. Thirdly, the selection lines did not differ among each other in shoal cohesion at any ontogenetic stage. Fourthly, dopamine turnover increased through ontogeny in a similar way for all lines while the serotonin turnover decreased in the large-harvested and control lines, but not in the small-harvested line. The large-harvested line also had higher serotonin turnover than controls at specific time periods. In conclusion, intensive size-selective mortality left an evolutionary legacy of asymmetric selection responses in the ontogeny of shoal cohesion and the underlying physiological mechanisms in experimentally harvested zebrafish in the laboratory.

Significant statement

The evolution of animal behavior can be affected by human activities both at behavioral and physiological levels, but causal evidence is scarce and mostly focusing on single life-stages. We studied whether and to what extent size-selective harvesting, a common selection pattern in fisheries, can be an evolutionary driver of the development of shoal cohesion during ontogeny. We used a multi-generation experiment with zebrafish to study cause-and-effects of opposing size-selection patterns. We quantified shoal cohesion, and serotonin and dopamine turnover in the brain. We found that shoal cohesion emerged as a collective personality trait and that behavioral and physiological responses were asymmetrical with respect to the opposing selection patterns.

Ecological ramifications of adaptation to size-selective mortality

Size-selective mortality due to harvesting is a threat to numerous exploited species, but how it affects the ecosystem remains largely unexplored. Here, we used a pond mesocosm experiment to assess how evolutionary responses to opposite size-selective mortality interacted with the environment (fish density and light intensity used as a proxy of resource availability) to modulate fish populations, prey community composition and ecosystem functions. We used medaka (Oryzias latipes) previously selected over 10 generations for small size (harvest-like selection; small-breeder line) or large size (large-breeder line), which displayed slow somatic growth and early maturity or fast somatic growth and late maturity, respectively. Large-breeder medaka produced more juveniles, which seemed to grow faster than small-breeder ones but only under high fish density. Additionally, large-breeder medaka had an increased impact on some benthic prey, suggesting expanded diet breadth and/or enhanced foraging abilities. As a consequence, increased light stimulated benthic algae biomass only in presence of large-breeder medaka, which were presumably better at controlling benthic grazers. Aggregated effect sizes at the community and ecosystem levels revealed that the ecological effects of medaka evolution were of similar magnitude to those induced by the environment and fish introduction. These findings indicate the important environmental dependency of evolutionary response to opposite size-selective mortality on higher levels of biological organizations.

Density-dependent natural selection mediates harvest-induced trait changes

Rapid life-history changes caused by size-selective harvesting are often interpreted as a response to direct harvest selection against a large body size. However, similar trait changes may result from a harvest-induced relaxation of natural selection for a large body size via density-dependent selection. Here, we show evidence of such density-dependent selection favouring large-bodied individuals at high population densities, in replicated pond populations of medaka fish. Harvesting, in contrast, selected medaka directly against a large body size and, in parallel, decreased medaka population densities. Five years of harvesting were enough for harvested and unharvested medaka populations to inherit the classically predicted trait differences, whereby harvested medaka grew slower and matured earlier than unharvested medaka. We show that this life-history divergence was not driven by direct harvest selection for a smaller body size in harvested populations, but by density-dependent natural selection for a larger body size in unharvested populations.

Unidirectional response to bidirectional selection on body size II. Quantitative genetics

Anticipating the genetic and phenotypic changes induced by natural or artificial selection requires reliable estimates of trait evolvabilities (genetic variances and covariances). However, whether or not multivariate quantitative genetics models are able to predict precisely the evolution of traits of interest, especially fitness-related, life history traits, remains an open empirical question. Here, we assessed to what extent the response to bivariate artificial selection on both body size and maturity in the medaka Oryzias latipes, a model fish species, fits the theoretical predictions. Three lines (Large, Small, and Control lines) were differentially selected for body length at 75 days of age, conditional on maturity. As maturity and body size were phenotypically correlated, this selection procedure generated a bi-dimensional selection pattern on two life history traits. After removal of nonheritable trends and noise with a random effect ("animal") model, the observed selection response did not match the expected bidirectional response. For body size, Large and Control lines responded along selection gradients (larger body size and stasis, respectively), but, surprisingly, the Small did not evolve a smaller body length and remained identical to the Control line throughout the experiment. The magnitude of the empirical response was smaller than the theoretical prediction in both selected directions. For maturity, the response was opposite to the expectation (the Large line evolved late maturity compared to the Control line, while the Small line evolved early maturity, while the opposite pattern was predicted due to the strong positive genetic correlation between both traits). The mismatch between predicted and observed response was substantial and could not be explained by usual sources of uncertainties (including sampling effects, genetic drift, and error in G matrix estimates).

Ecosystem consequences of multi-trait response to environmental changes in Japanese medaka, Oryzias latipes

Intraspecific trait variation has large effects on the ecosystem and is greatly affected by human activities. To date, most studies focused on single-trait analyses, while considering multiple traits is expected to better predict how an individual interacts with its environment. Here, we used a mesocosm experiment with fish Oryzias latipes to test whether individual growth, boldness and functional traits (feeding rate and stoichiometric traits) formed one functional pace-of-life syndrome (POLS). We then tested the effects of among-individual mean and variance of fish functional POLSs within mesocosms on invertebrate community (e.g. zoobenthos and zooplankton abundances) and ecosystem processes (e.g. ecosystem metabolism, algae stock, nutrient concentrations). Stoichiometric traits correlated with somatic growth and behaviours, forming two independent functional POLS (i.e. two major covariance axes). Mean values of the first syndrome were sex- and environment-dependent and were associated with (i) long-term (10 generations; 4 years) selection for small or large body size resulting in contrasting life histories and (ii) short-term (6 weeks) effects of experimental treatments on resource availability (through manipulation of light intensity and interspecific competition). Specifically, females and individuals from populations selected for a small body size presented fast functional POLS with faster growth rate, higher carbon body content and lower boldness. Individuals exposed to low resources (low light and high competition) displayed a slow functional POLS. Higher mesocosm mean and variance values in the second functional POLS (i.e. high feeding rate, high carbon:nitrogen body ratio, low ammonium excretion rate) were associated to decreased prey abundances, but did not affect any of the ecosystem processes. We highlighted the presence of functional multi-trait covariation in medaka, which were affected by sex, long-term selection history and short-term environmental conditions, that ultimately had cascading ecological consequences. We stressed the need for applying this approach to better predict ecosystem response to anthropogenic global changes.