Evolution of phenotypic polymorphism in symbiont-pairing in plant-fungal symbiosis

This is Mathematica code to model the evolution of polyphenism in plants (i.e., the ability of a single genotype to randomly develop into individuals with different phenotypes). We are interested in the trait of roots allocated to two different fungal partners when the benefits of a particular ratio of partners vary with the environment. We use both an analytical model and a simulation to investigate the conditions for polyphenism to evolve and evolutionary endpoints -- the root allocation phenotypes that evolved to be expressed and the probability that each one is expressed. Along with the analytical model and simulation code, we also include saved simulation results as .csv files so that the results can be analyzed without needing to re-run the simulations. The saved simulation results include the evolved traits (root allocation and phenotype frequency) and the time it took for populations to evolve away from monophenism.
, , # Evolution of phenotypic polymorphism in symbiont-pairing in plant-fungal symbiosis

Dryad DOI: https://doi.org/10.5061/dryad.6m905qgcs

Authors: Alexandra L. Brown and Stephen Proulx
Contact: alexandra_brown (berkeley.edu)

This is Mathematica code that models the evolution of polyphenism in plants, in particular, polyphenism in the ratio of roots allocated to two different fungal partners in a variable environment. We use both an analytical model and a simulation to investigate the conditions for polyphenism to evolve and the endpoints of evolution. Along with the analytical model and simulation code, there are saved simulation results as .wdx files, which can be loaded from within `evolution_simulations.wls` and `evolution_simulations_alternate_pfix.wls` instead of running the simulations, if desired.

To produce the results in the manuscript, the code was run in Mathematica version 13.

## Analytical model

* `model_and_analysis.wls` contains the analytical model and its analysis....,