Selection for small body size favours contrasting sex-specific life histories, boldness and feeding in medaka, Oryzias latipes

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.

Size-selective mortality fosters ontogenetic changes in collective risk-taking behaviour in zebrafish, Danio rerio

Size-selective mortality is common in fish populations and can operate either in a positive size-selective fashion by harvesting larger-than-average fish or be negatively size-selective by harvesting smaller-than-average fish. Through various mechanisms (like genetic correlations among behaviour and life-history traits or direct selection on behaviour co-varying with growth rate or size-at-maturation), size-selection can result in evolutionary changes in behavioural traits. Theory suggests that both positive and negative size-selection without additional selection on behaviour favours boldness, while evolution of shyness is possible if the largest fish are harvested. Here we examined the impact of size-selective mortality on collective boldness across ontogeny using three experimental lines of zebrafish (Daniorerio) generated through positive (large-harvested), negative (small-harvested) and random (control line) size-selective mortality for five generations and then relaxed selection for 10 generations to examine evolutionarily fixed outcomes. We measured collective risk-taking during feeding (boldness) under simulated aerial predation threat, and across four contexts in presence/absence of a cichlid. Boldness decreased across ontogeny under aerial predation threat, and the small-harvested line was consistently bolder than controls. The large and small-harvested lines showed higher behavioural plasticity as larvae and developed personality earlier compared to the controls. The large-harvested line showed increased variability and plasticity in boldness throughout ontogeny. In the presence of a live predator, fish did not differ in boldness in three contexts compared to the controls, but the large-harvested line showed reduced behavioural plasticity across contexts than controls. Our results confirmed theory by demonstrating that size-selective harvesting evolutionarily alters collective boldness and its variability and plasticity.

Environmental temperature during early life affects the personality of mosquitofish in adulthood

Personality has been observed in a variety of animal taxa with important implications in ecology and evolution. Exploring the influence of environmental temperature during early life on personality could help to understand the ontogeny of this phenotypic trait in animals. In this study, we reared newborn mosquitofish Gambusia affinis at high (30°C) and low (25°C) water temperatures and measured their shyness and exploration upon sexual maturity. We tested the repeatability of each behavioral trait; the correlation between them; and the effects of rearing temperature, sex, and body length on the behaviors. When growing up at low temperatures, female fish exhibited repeatability in shyness and exploration, and males exhibited marginal repeatability in shyness. However, neither of the 2 behaviors were repeatable when the fish were reared at high temperatures. There was a negative correlation between shyness and exploration, indicating that the 2 behaviors comprise a behavioral syndrome in this species. Mosquitofish reared at high temperatures were more explorative than those reared at low temperatures, while there was no difference in shyness between the 2 treatments. Body length and sex had no significant effects on the average values of the 2 behaviors. The results indicate that environmental temperature during early life could shape the personality of mosquitofish and modify the average of the behavioral traits. These findings might provide insights to understand the ontogeny of animal personality and how changes in environmental temperature influence animal dispersal by shaping their personality.

Density influences the heritability and genetic correlations of fish behaviour under trawling-associated selection

Fishing-associated selection is one of the most important human-induced evolutionary pressures for natural populations. However, it is unclear whether fishing leads to heritable phenotypic changes in the targeted populations, as the heritability and genetic correlations of traits potentially under selection have received little attention. In addition, phenotypic changes could arise from fishing-associated environmental effects, such as reductions in population density. Using fish reared at baseline and reduced group density and repeatedly harvested by simulated trawling, we show that trawling can induce direct selection on fish social behaviour. As sociability has significant heritability and is also genetically correlated with activity and exploration, trawling has the potential to induce both direct selection and indirect selection on a variety of fish behaviours, potentially leading to evolution over time. However, while trawling selection was consistent between density conditions, the heritability and genetic correlations of behaviours changed according to the population density. Fishing-associated environmental effects can thus modify the evolutionary potential of fish behaviour, revealing the need to use a more integrative approach to address the evolutionary consequences of fishing.

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.

Size-selective mortality induces evolutionary changes in group risk-taking behaviour and the circadian system in a fish

Intensive and trait-selective mortality of fish and wildlife can cause evolutionary changes in a range of life-history and behavioural traits. These changes might in turn alter the circadian system due to co-evolutionary mechanisms or correlated selection responses both at behavioural and molecular levels, with knock-on effects on daily physiological processes and behavioural outputs. We examined the evolutionary impact of size-selective mortality on group risk-taking behaviour and the circadian system in a model fish species. We exposed zebrafish Danio rerio to either large or small size-selective harvesting relative to a control over five generations, followed by eight generations during which harvesting was halted to remove maternal effects. Size-selective mortality affected fine-scale timing of behaviours. In particular, small size-selective mortality, typical of specialized fisheries and gape-limited predators targeting smaller size classes, increased group risk-taking behaviuor during feeding and after simulated predator attacks. Moreover, small size-selective mortality increased early peaks of daily activity as well as extended self-feeding daily activity to the photophase compared to controls. By contrast large size-selective mortality, typical of most wild capture fisheries, only showed an almost significant effect of decreasing group risk-taking behaviour during the habituation phase and no clear changes in fine-scale timing of daily behavioural rhythms compared to controls. We also found changes in the molecular circadian core clockwork in response to both size-selective mortality treatments. These changes disappeared in the clock output pathway because both size-selected lines showed similar transcription profiles. This switch downstream to the molecular circadian core clockwork also resulted in similar overall behavioural rhythms (diurnal swimming and self-feeding in the last hours of darkness) independent of the underlying molecular clock. To conclude, our experimental harvest left an asymmetrical evolutionary legacy in group risk-taking behaviour and in fine-scale daily behavioural rhythms. Yet, the overall timing of activity showed evolutionary resistance probably maintained by a molecular switch. Our experimental findings suggest that size-selective mortality can have consequences for behaviour and physiological processes.

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).

Stoichiometric consequences of size-selective mortality: An experimental test using the Japanese medaka (Oryzias latipes)

The determinants of intraspecific stoichiometric variation remain difficult to elucidate due to their multiple origins (e.g. genetic vs. environmental) and potential interactive effects. We evaluated whether two size-selected lines of medaka (Oryzias latipes) with contrasted life-history strategies (small- and large-breeder lines with slow growth and early maturity vs. fast growth and late maturity) differed in their organismal stoichiometry (percentage and ratios of carbon [C], nitrogen [N] and phosphorus [P]) in a mesocosm experiment. We also tested how size-selection interacted with environmental conditions (i.e. two levels of fish density and light intensity), body condition and sex. Results showed that large-breeder fish were significantly N-enriched compared to small-breeders, while the two size-selected lines did not differ in body P composition. Size-selection interacted with density - high density only affected small-breeders leading to decreasing %C and C: N - and with sex - large-breeder females had higher %C and C:N values than large-breeder males. Finally, C:P and N:P ratios increased with body condition due to decreasing %P. Overall, our results show that the ecological consequences of size-selective mortality extend to organismal stoichiometry and may, from there, change nutrient cycling and ecosystem functioning.

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.