Transcriptome variation in banded newt (Ommatotriton vittatus) during its life cycle and habitat transition

Israel represents the southern limit of the distribution of the banded newt (Ommatotriton vittatus). The life cycle of O. vittatus includes several distinct phases: eggs, aquatic larvae, a terrestrial phase and an aquatic reproductive phase. We investigated differences in gene expression during the life cycle and transition of banded newts between terrestrial and aquatic habitats using mRNA-seq. We identified ∼10 k genes that were differentially expressed (DE) in one of the pairwise comparisons between 3 groups: 1 - terrestrial newts (males and females), 2 - aquatic newts (males and females), 3 - aquatic larvae before metamorphosis. The groups were clearly defined by Principal Components Analysis (PCA). The greatest difference was between aquatic newts (males and females) and aquatic larvae: ∼7.4 k DE genes. Of special interest were the ∼2.4 k genes DE between the aquatic and terrestrial phenotypes. These included prominent candidates with known roles in kidney function (uromodulin homologs were strongly associated with aquatic lifestyle), tissue structure (keratins), and the thyroid hormone signaling modulator DUOXA1. Additional developmental and metabolic pathways overrepresented among the identified DE genes included "epidermis development", "nervous system development", "nucleotide-sugar biosynthesis". Overall, both metamorphosis and environmental transition of banded newts involve extensive transcriptomic remodeling involving developmental, metabolic, and cellular pathways. Understanding the roles of these pathways and individual genes is instrumental for studies of transition between habitats, especially those affected by climate change. Furthermore, the phenotypic flexibility of the newt and the underlying regulation of gene expression can shed light on the evolution of terrestrial vertebrates.