Genomic signatures associated with transitions to viviparity in Cyprinodontiformes

The transition from oviparity to viviparity has occurred independently over a hundred times across vertebrates, presenting a compelling case of phenotypic convergence. However, whether repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common genetic signature in genomic divergence remains unknown. Whilst investigations into the evolution of viviparity have demonstrated striking similarity among the genes and pathways involved across vertebrate groups, quantitative tests for genome-wide convergence provide ambivalent answers. Here, we investigate molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo and publicly-available genomes of viviparous and oviparous species to quantify molecular convergence across coding and non-coding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, statistical power and biological confounds (hemiplasy and introgression) could constrain our ability to detect correlated evolution. We therefore also identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. While we detected no evidence of positive or relaxed selection for these genes in branches associated with the evolution of viviparity in Cyprinodontiformes, motif-enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development and immunity occurred alongside the evolution of viviparity. Overall, our findings indicate that an excess of molecular convergence, at any level, is not strongly associated with independent transitions to viviparity in these fish.