Identification of a Novel Uromodulin-Like Gene Related to Predator-Induced Bulgy Morph in Anuran Tadpoles by Functional Microarray Analysis

Tadpoles of the anuran species Rana pirica can undergo predator-specific morphological responses. Exposure to a predation threat by larvae of the salamander Hynobius retardatus results in formation of a bulgy body (bulgy morph) with a higher tail. The tadpoles revert to a normal phenotype upon removal of the larval salamander threat. Although predator-induced phenotypic plasticity is of major interest to evolutionary ecologists, the molecular and physiological mechanisms that control this response have yet to be elucidated. In a previous study, we identified various genes that are expressed in the skin of the bulgy morph. However, it proved difficult to determine which of these were key genes in the control of gene expression associated with the bulgy phenotype. Here, we show that a novel gene plays an important role in the phenotypic plasticity producing the bulgy morph. A functional microarray analysis using facial tissue samples of control and bulgy morph tadpoles identified candidate functional genes for predator-specific morphological responses. A larger functional microarray was prepared than in the previous study and used to analyze mRNAs extracted from facial and brain tissues of tadpoles from induction-reversion experiments. We found that a novel uromodulin-like gene, which we name here pirica, was up-regulated and that keratin genes were down-regulated as the period of exposure to larval salamanders increased. Pirica consists of a 1296 bp open reading frame, which is putatively translated into a protein of 432 amino acids. The protein contains a zona pellucida domain similar to that of proteins that function to control water permeability. We found that the gene was expressed in the superficial epidermis of the tadpole skin.