Transcriptome analysis of predator- and prey-induced phenotypic plasticity in the Hokkaido salamander (Hynobius retardatus)

Predator- and prey-induced phenotypic plasticity is widely observed among amphibian species. Although ecological factors inducing diverse phenotypic responses have been extensively characterized, we know little about the molecular bases of variation in phenotypic plasticity. Larvae of the Hokkaido salamander, Hynobius retardatus, exhibit two distinct morphs: the presence of their prey, Rana pirica tadpoles, induces a broad-headed attack morph, and the presence of predatory dragonfly nymphs (Aeshna nigroflava) induces a defence morph with enlarged external gills and a high tail. To compare the genes involved in predator- and prey-induced phenotypic plasticity, we carried out a de novo transcriptome analysis of Hokkaido salamander larvae exposed to either prey or predator individuals. First, we found that the number of genes involved in the expression of the defence morph was approximately five times the number involved in the expression of the attack morph. This result is consistent with...