Across the tree of life, species have repeatedly evolved similar phenotypes. While well-studied for ecological traits, there is also evidence for recurrent evolution of sexually selected traits. Swordtail fish (Xiphophorus) are a classic model system for studying sexual selection, and female Xiphophorus exhibit strong mate preferences for large male body size and a range of sexually dimorphic ornaments. Interestingly, sexually selected traits have also been lost multiple times in the genus. However, there has been uncertainty over the number of losses of ornamentation and large body size because phylogenetic relationships between species in this group have historically been controversial, partially due to prevalent gene flow. Here, we use whole-genome sequencing approaches to re-examine phylogenetic relationships within a Xiphophorus clade that varies in the presence and absence of sexually selected traits. Using wild-caught individuals, we determine the phylogenetic placement of a small, unornamented species, X. continens, confirming an additional loss of ornamentation and large body size in the clade. With these revised phylogenetic relationships, we analyze evidence for coevolution between body size and other sexually selected traits using phylogenetic comparative methods. These results provide insights into the evolutionary pressures driving the recurrent loss of suites of sexually selected traits.

The remarkable diversity of sex determination mechanisms known in fish may be fuelled by exceptionally high rates of sex chromosome turnovers or transitions. However, the evolutionary causes and genomic mechanisms underlying this variation and instability are yet to be understood. Here we report on an over 30-year evolutionary experiment in which we tested the genomic consequences of hybridisation and selection between two Xiphophorus fish species with different sex chromosome systems. We find that introgression and imposing selection for pigmentation phenotypes results in the retention of an unexpectedly large maternally derived genomic region. During the hybridisation process, the sex-determining region of the X chromosome from one parental species was translocated to an autosome in the hybrids leading to the evolution of a new sex chromosome. Our results highlight the complexity of factors contributing to patterns observed in hybrid genomes, and we experimentally demonstrate that hybridisation can catalyze rapid evolution of a new sex chromosome.

Sex determination systems are diverse, especially among fish, and include genetic and/or environmental components. Unexpectedly for such a basic aspect of development, sex determination systems change rapidly during evolution and gonadal fate is not ultimate, being actively maintained lifelong. Here, sequences of expressed genes involved in maintenance of gonad identity and reproduction processes were obtained through transcriptome assembly of the brain-gonadal axis tissues of a freshwater fish inhabiting highly variable environments, the gonochoristic Iberian fish Squalius pyrenaicus. Through Illumina total RNA-sequencing, male and female transcriptomes of brain and gonad tissues were assembled with Trans-ABySS software and merged to produce a more comprehensive S. pyrenaicus transcriptome. Coding sequences (CDS) predicted by TransDecoder were annotated using blastx. By means of read mapping against the reference transcriptome and CDS datasets, using Bowtie2, the accuracy of read mapping was assessed. This first endemic Iberian cyprinid transcriptome of organs involved in reproduction processes may serve as a valuable genomic resource for studying sexual mechanisms and other aspects of evolution, such as speciation and responses to environmental changes, and may be a useful tool for conservation studies since S. pyrenaicus is an endangered species.

Next-generation sequencing (NGS) techniques are now key tools in the detection of population genomic and gene expression differences in a large array of organisms. However, so far few studies have utilized such data for phylogenetic estimations. Here, we use NGS data obtained from genome-wide restriction site-associated DNA (RAD) (∼66000 SNPs) to estimate the phylogenetic relationships among all 26 species of swordtail and platyfish (genus Xiphophorus) from Central America. Past studies, both sequence and morphology-based, have differed in their inferences of the evolutionary relationships within this genus, particularly at the species-level and among monophyletic groupings. We show that using a large number of markers throughout the genome, we are able to infer the phylogenetic relationships with unparalleled resolution for this genus. The relationships among all three major clades and species within each of them are highly resolved and consistent under maximum likelihood, Bayesian inference and maximum parsimony. However, we also highlight the current cautions with this data type and analyses. This genus exhibits a particularly interesting evolutionary history where at least two species may have arisen through hybridization events. Here, we are able to infer the paternal lineages of these putative hybrid species. Using the RAD-marker-based tree we reconstruct the evolutionary history of the sexually selected sword trait and show that it may have been present in the common ancestor of the genus. Together our results highlight the outstanding capacity that RAD sequencing data has for resolving previously problematic phylogenetic relationships, particularly among relatively closely related species.