Experimental and theoretical support for costs of plasticity and phenotype in a nematode cannibalistic trait

Developmental plasticity is the ability of a genotype to express multiple phenotypes under different environmental conditions and has been shown to facilitate the evolution of novel traits. However, while the associated cost of plasticity, i.e., the loss in fitness due to the ability to express plasticity in response to environmental change, and the cost of phenotype, i.e., the loss of fitness due to expressing a fixed phenotype across environments, have been theoretically predicted, empirically such costs remain poorly documented and little understood. Here, we use a plasticity model system, hermaphroditic nematode Pristionchus pacificus, to experimentally measure these costs in wild isolates under controlled laboratory conditions. P. pacificus can develop either a bacterial feeding or predatory mouth morph in response to different external stimuli, with natural variation of mouth-morph ratios between strains. We first demonstrated the cost of phenotype by analyzing fecundity and developmental speed in relation to mouth morphs across the P. pacificus phylogenetic tree. Then, we exposed P. pacificus strains to two distinct microbial diets that induce strain-specific mouth-form ratios. Our results indicate that the plastic strain does shoulder a cost of plasticity, i.e., the diet-induced predatory mouth morph is associated with reduced fecundity and slower developmental speed. In contrast, the non-plastic strain suffers from the cost of phenotype since its phenotype does not change to match the unfavorable bacterial diet but shows increased fitness and higher developmental speed on the favorable diet. Furthermore, using a stage-structured population model based on empirically derived life history parameters, we show how population structure can alleviate the cost of plasticity in P. pacificus. The results of the model illustrate the extent to which the costs associated with plasticity and its effect on competition depend on ecological factors. This study provides support for costs of plasticity and phenotype based on empirical and modeling approaches.

Methods
This dataset contains the following scripts:

digital_worm.py: The model written in Python 3.10
main_figs_model.ipynb: The code with instructions to generate the main modeling figures.
main_figs_exp.ipynb: The code with instructions to generate the main experimental figures, including the Bayesian statistical analyses associated with these figures.
raw_data.zip: The raw data required to generate main and supp figures. 
supp_figs.ipynb: The code with instructions to generat the supp figures
digital_worm_n_dev_dif.py: A verison of the model written in Python 3.10 in which the developmental speeds remain unchanged upon switching the diet. Used in supp_figs.ipynb. 

The dataset can be found in raw_data.zip