Comparative Transcriptomic Analyses Identify Candidate Genes for Convergent Reproductive Shifts in a Bimodal...

Shifts in reproductive mode represent key evolutionary innovations that shape species' life histories and evolutionary trajectories. Species showing bimodal reproductive strategies with multiple independent origins offer a rare opportunity to gain insights into the adaptive processes and mechanisms underlying convergent traits. The fire salamander, Salamandra salamandra, is the only amphibian exhibiting intraspecific variation in reproductive mode across multiple independent reproductive shifts, enabling investigation of the transition between larviparity (females give birth to aquatic larvae) and pueriparity (females give birth to fully developed terrestrial juveniles) within a single species and across different timescales. Pueriparity is an adaptive innovation that skips the aquatic larval stage, allowing individuals to exploit habitats with no available water bodies. The fire salamander is larviparous across most of its range, but pueriparity has evolved independently at least three times: once in the early Pleistocene within S. s. bernardezi in the mountains of northern Spain, and more recently on two land-bridge islands (NW Spain) inhabited by S. s. gallaica. To identify candidate genes associated with these distinct reproductive modes, we compared gene expression profiles of the uterus and oviduct of pregnant females across two independent evolutionary transitions using RNA-sequencing. We detected shared changes in maternal gene expression among pueriparous S. s. bernardezi and S. s. gallaica relative to their larviparous counterparts, in addition to differences unique to each independent evolutionary transition. Functional enrichment analyses indicated that differentially expressed genes were associated with reproductive timing, angiogenesis, and maternal signalling, consistent with the phenotypic differences observed in the uterine environment and embryonic development between the two reproductive modes. This study represents an important first step towards understanding the genomic basis of the evolution of pueriparity in a remarkable bimodal reproductive system, and provides transcriptomic resources and candidate genes for future research into the genomic architecture underlying this poorly understood adaptive trait.