Given that selection removes genetic variance from evolving populations, thereby reducing exploration opportunities, it is important to find mechanisms that create genetic variation without the disruption of adapted genes and genomes caused by random mutation. Just such an alternative is offered by random epigenetic error, a developmental process that acts on materials and parts expressed by the genome. In this system of embodied computational evolution, simulated within a physics engine, epigenetic error was instantiated in an explicit genotype-to-phenotype map as transcription error at the initiation of gene expression. The hypothesis was that transcription error would create genetic variance by shielding genes from the direct impact of selection, creating, in the process, masquerading genomes. To test this hypothesis, populations of simulated embodied biorobots and their developmental systems were evolved under steady directional selection as equivalent rates of random mutation and random transcriptional error were covaried systematically in an 11 × 11 fully factorial experimental design. In each of the 121 different experimental conditions (unique combinations of mutation and transcription error), the same set of 10 randomly created replicate populations of 60 individuals were evolved. Selection for the improved locomotor behavior of individuals led to increased mean fitness of populations over 100 generations at nearly all levels and combinations of mutation and transcription error. When the effects of both types of error were partitioned statistically, increasing transcription error was shown to increase the final genetic variance of populations, incurring a fitness cost but acting on variance independently and differently from genetic mutation. Thus, random epigenetic errors in development feed back through selection of individuals with masquerading genomes to the population’s genetic variance over generational time. Random developmental processes offer an additional mechanism for exploration by increasing genetic variation in the face of steady, directional selection.

To fully understand the evolution of complex morphologies, analyses cannot stop at selection: It is essential to investigate the roles and interactions of multiple processes that drive evolutionary outcomes. The challenges of undertaking such analyses have affected both evolutionary biologists and evolutionary roboticists, with their common interests in complex morphologies. In this paper, we present analytical techniques from evolutionary biology, selection gradient analysis and morphospace walks, and we demonstrate their applicability to robot morphologies in analyses of three evolutionary mechanisms: randomness (genetic mutation), development (an explicitly implemented genotype-to-phenotype map), and selection. In particular, we applied these analytical techniques to evolved populations of simulated biorobots—embodied robots designed specifically as models of biological systems, for the testing of biological hypotheses—and we present a variety of results, including analyses that do all of the following: illuminate different evolutionary dynamics for different classes of morphological traits; illustrate how the traits targeted by selection can vary based on the likelihood of random genetic mutation; demonstrate that selection on two selected sets of morphological traits only partially explains the variance in fitness in our biorobots; and suggest that biases in developmental processes could partially explain evolutionary dynamics of morphology. When combined, the complementary analytical approaches discussed in this paper can enable insight into evolutionary processes beyond selection and thereby deepen our understanding of the evolution of robotic morphologies.