Comprehensive phylogenetic analysis of all species of swordtails and platies (Pisces: Genus Xiphophorus) uncovers a hybrid origin of a swordtail fish, Xiphophorus monticolus, and demonstrates that the sexually selected sword originated in the ancestral lineage of the genus, but was lost again secondarily

Segregation Between an Ornamental and a Disease Driver Gene Provides Insights Into Pigment Cell Regulation

Genetic interactions are adaptive within a species. Hybridization can disrupt such species-specific genetic interactions and creates novel interactions that alter the hybrid progeny overall fitness. Hybrid incompatibility, which refers to degenerative genetic interactions that decrease the overall hybrid survival and sterility, is one of the results from combining two diverged genomes in hybrids. The discovery of spontaneous lethal tumorigenesis and underlying genetic interactions in select hybrids between diverged Xiphophorus species showed that lethal pathological process can result from degenerative genetic interactions. Such genetic interactions leading to lethal phenotype are thought to shield gene flow between diverged species. However, hybrids between certain Xiphophorus species do not develop such tumors. Here we report the identification of a locus residing in the genome of one Xiphophorus species that represses an oncogene from a different species. Our finding provides insights into normal and pathological pigment cell development, regulation and a molecular mechanism in hybrid incompatibility.

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

Segregation between an ornamental and a disease driver gene provides insights into pigment cell regulation

Genetic interactions are adaptive within a species. Hybridization can disrupt such species-specific genetic interactions and creates novel interactions that alter the hybrid progeny overall fitness. Hybrid incompatibility, which refers to degenerative genetic interactions that decrease the overall hybrid survival, is one of the results from combining two diverged genomes in hybrids. The discovery of spontaneous lethal tumorigenesis and underlying genetic interactions in select hybrids between diverged Xiphophorus species showed that lethal pathological process can result from degenerative genetic interactions. Such genetic interactions leading to lethal phenotype are thought to shield gene flow between diverged species. However, hybrids between certain Xiphophorus species do not develop such tumors. Here we report the identification of a locus residing in the genome of one Xiphophorus species that represses an oncogene from a different species. Our finding provides insights into normal and pathological pigment cell development, regulation and molecular mechanism in hybrid incompatibility.

Significance

The Dobzhansky-Muller model states epistatic interactions occurred between genes in diverged species underlies hybrid incompatibility. There are a few vertebrate interspecies hybrid cases that support the Dobzhansky-Muller model. This study reports a fish hybrid system where incompatible genetic interactions are involved in neuronal regulation of pigment cell biology, and also identified a novel point of regulation for pigment cells.

Disentangling historical relationships within Poeciliidae (Teleostei: Cyprinodontiformes) using ultraconserved elements

Poeciliids (Cyprinodontiformes: Poeciliidae), commonly known as livebearers, are popular fishes in the aquarium trade (e.g., guppies, mollies, swordtails) that are widely distributed in the Americas, with 274 valid species in 27 genera. This group has undergone various taxonomic changes recently, spurred by investigations using traditional genetic markers. Here we used over 1,000 ultraconserved loci to infer the relationships within Poeciliidae in the first attempt at understanding their diversification based on genome-scale data. We explore gene tree discordance and investigate potential incongruence between concatenation and coalescent inference methods. Our aim is to examine the influence of incomplete lineage sorting and reticulate evolution on the poeciliids' evolutionary history and how these factors contribute to the observed gene tree discordace. Our concatenated and coalescent phylogenomic inferences recovered four major clades within Poeciliidae. Most supra-generic level relationships we inferred were congruent with previous molecular studies, but we found some disagreements; the Middle American taxa Phallichthys and Poecilia (Mollienesia) were recovered as non-monophyletic, and unlike other recent molecular studies, we recovered Brachyrhaphis as monophyletic. Our study is the first to provide signatures of reticulate evolution in Poeciliidae at the family level; however, continued finer-scale investigations are needed to understand the complex evolutionary history of the family along with a much-needed taxonomic re-evaluation.

Dimerization of melanocortin 4 receptor controls puberty onset and body size polymorphism

Xiphophorus fish exhibit a clear phenotypic polymorphism in puberty onset and reproductive strategies of males. In X. nigrensis and X. multilineatus, puberty onset is genetically determined and linked to a melanocortin 4 receptor (Mc4r) polymorphism of wild-type and mutant alleles on the sex chromosomes. We hypothesized that Mc4r mutant alleles act on wild-type alleles by a dominant negative effect through receptor dimerization, leading to differential intracellular signaling and effector gene activation. Depending on signaling strength, the onset of puberty either occurs early or is delayed. Here, we show by Förster Resonance Energy Transfer (FRET) that wild-type Xiphophorus Mc4r monomers can form homodimers, but also heterodimers with mutant receptors resulting in compromised signaling which explains the reduced Mc4r signaling in large males. Thus, hetero- vs. homo- dimerization seems to be the key molecular mechanism for the polymorphism in puberty onset and body size in male fish.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Complete mitochondrial genome of a livebearing freshwater fish (Cyprinodontiformes: Poeciliidae): Poecilia parae

Members of the fish family Poeciliidae (livebearing 'tooth-carps') have historically been used as models in medical research, behavior ecology, and biological control. This group of primarily freshwater fishes is highly tolerant to environmental factors such as salinity and warm temperatures and includes some invasive species. Here, we present the mitochondrial genome of Poecilia parae. A representative of this species was obtained from Suriname. The complete mitochondrial genome was sequenced using Oxford Nanopore technology and is 16,559 bp long. The genome contains 13 protein-coding genes, two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and one control region (D-loop). Phylogenetic analysis yielded topologies similar to those previously published. The data generated here will be useful in future studies of comparative biology and those utilizing environmental DNA (eDNA).

Assessment of Various Standard Fish Diets on Growth and Fecundity of Platyfish (Xiphophorus maculatus) and Medaka (Oryzias latipes)

Several small freshwater fish species are utilized as models for human conditions and disease in biomedical research. Research animal diets are generally tailored to optimize growth, fecundity, and produce healthy research animals. However, a lack of reference diets presents a barrier in comparative studies between aquatic animal models and even among laboratories using the same species. Therefore, the objective of this study was to determine feeding regime and dietary effects on growth and fecundity in two commonly used freshwater fish, platyfish and medaka. From 1 through 6 months of age, platyfish and medaka were fed one of three feeding regime/diets: (1) our custom feeding regime consists of commercial flake food, beef liver paste, and live brine shrimp (CON); (2) a commercially available zebrafish diet, Gemma (GEM); and (3) a laboratory defined reference feeding regime (WAT). Weight, size, brood numbers, and survival rates for both species were measured monthly. Numbers of platyfish fry and hatch rate of medaka embryos were also determined. We observed that custom feeding regime (CON) fed platyfish and medaka grew larger, exhibited a higher survival rate, and had higher fecundity than WAT or GEM fed fish. These observations suggest that diets and regimes designed for zebrafish are not optimal to maintain platyfish or medaka. Thus, base diets, with clearly defined components and regimes, need to be developed with compositions that can be adjusted in a species-specific manner.

Secondary reversion to sexual monomorphism associated with tissue-specific loss of doublesex expression

Animal evolution is characterized by frequent turnover of sexually dimorphic traits-new sex-specific characters are gained, and some ancestral sex-specific characters are lost, in many lineages. In insects, sexual differentiation is predominantly cell autonomous and depends on the expression of the doublesex (dsx) transcription factor. In most cases, cells that transcribe dsx have the potential to undergo sex-specific differentiation, while those that lack dsx expression do not. Consistent with this mode of development, comparative research has shown that the origin of new sex-specific traits can be associated with the origin of new spatial domains of dsx expression. In this report, we examine the opposite situation-a secondary loss of the sex comb, a male-specific grasping structure that develops on the front legs of some drosophilid species. We show that while the origin of the sex comb is linked to an evolutionary gain of dsx expression in the leg, sex comb loss in a newly identified species of Lordiphosa (Drosophilidae) is associated with a secondary loss of dsx expression. We discuss how the developmental control of sexual dimorphism affects the mechanisms by which sex-specific traits can evolve.

Species delineation and systematics of a hemiclonal hybrid complex in Australian freshwaters (Gobiiformes: Gobioidei: Eleotridae: Hypseleotris)

The rivers of southeastern Australia host a species complex within the carp gudgeon genus Hypseleotris that includes parental species and hemiclonal hybrid lineages. These hemiclones can be difficult to distinguish from their parent taxa, making delineation of species unusually difficult. We approach this historical taxonomic problem by using single nucleotide polymorphism (SNP) genotyping to distinguish individuals of each species and hemiclones, enabling us to quantify the variation among evolutionary lineages and assign names to the species. Hypseleotris klunzingeri remains valid and does not have any hemiclones. We describe Hypseleotris bucephala and Hypseleotris gymnocephala from the Murray-Darling Basin and Hypseleotris acropinna from the Murray-Darling as well as eastern coastal streams north of the Mary River, part of the range attributed to H. galii. We further split H. galii to distinguish a species from the Mary River, Hypseleotris moolooboolaensis. We designate a neotype and redescribe H. galii due to uncertainty about the source and species identity of specimens used in the original description. We reconcile previous taxonomies, provide new common names for parental species, and advocate using the scientific names of both parents when referring to the hemiclone hybrids to avoid confusion with previous common names that did not distinguish parental taxa and hemiclones.

The complete mitochondrial genome of the Variable Platyfish Xiphophorus variatus

We present the complete mitochondrial genome sequence of the Variable Platyfish, Xiphophorus variatus (Meek 1904) (Cyprinodontiformes: Poeciliidae). The genome consists of 16,624 bp which encodes 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and 1 control region. Genome-wide nucleotide composition is 27.79% adenine, 31.11% cytosine, 15.63% guanine, and 25.48% thymine. The X. variatus mitochondrial genome shares similar GC content and identical gene order and gene strand location with other members of Poeciliidae. The sequence presented herein will be of utility for future phylogenetic and biomedical research and for designing primers for species detection from environmental DNA samples.

Rampant Genome-Wide Admixture across the Heliconius Radiation

How frequent is gene flow between species? The pattern of evolution is typically portrayed as a phylogenetic tree, yet gene flow between good species may be an important mechanism in diversification, spreading adaptive traits and leading to a complex pattern of phylogenetic incongruence. This process has thus far been studied mainly among a few closely related species, or in geographically restricted areas such as islands, but not on the scale of a continental radiation. Using a genomic representation of 40 out of 47 species in the genus, we demonstrate that admixture has played a role throughout the evolution of the charismatic Neotropical butterflies Heliconius. Modeling of phylogenetic networks based on the exome uncovers up to 13 instances of interspecific gene flow. Admixture is detected among the relatives of Heliconius erato, as well as between the ancient lineages leading to modern clades. Interspecific gene flow played a role throughout the evolution of the genus, although the process has been most frequent in the clade of Heliconius melpomene and relatives. We identify Heliconius hecalesia and relatives as putative hybrids, including new evidence for introgression at the loci controlling the mimetic wing patterns. Models accounting for interspecific gene flow yield a more complete picture of the radiation as a network, which will improve our ability to study trait evolution in a realistic comparative framework.

The Fish Family Poeciliidae as a Model to Study the Evolution and Diversification of Regenerative Capacity in Vertebrates

The capacity of regenerating a new structure after losing an old one is a major challenge in the animal kingdom. Fish have emerged as an interesting model to study regeneration due to their high and diverse regenerative capacity. To date, most efforts have focused on revealing the mechanisms underlying fin regeneration, but information on why and how this capacity evolves remains incomplete. Here, we propose the livebearing fish family Poeciliidae as a promising new model system to study the evolution of fin regeneration. First, we review the current state of knowledge on the evolution of regeneration in the animal kingdom, with a special emphasis on fish fins. Second, we summarize recent advances in our understanding of the mechanisms behind fin regeneration in fish. Third, we discuss potential evolutionary pressures that may modulate the regenerative capacity of fish fins and propose three new theories for how natural and sexual selection can lead to the evolution of fin regeneration: (1) signaling-driven fin regeneration, (2) predation-driven fin regeneration, and (3) matrotrophy-suppressed fin regeneration. Finally, we argue that fish from the family Poeciliidae are an excellent model system to test these theories, because they comprise of a large variety of species in a well-defined phylogenetic framework that inhabit very different environments and display remarkable variation in reproductive traits, allowing for comparative studies of fin regeneration among closely related species, among populations within species or among individuals within populations. This new model system has the potential to shed new light on the underlying genetic and molecular mechanisms driving the evolution and diversification of regeneration in vertebrates.

The Developmental and Genetic Architecture of the Sexually Selected Male Ornament of Swordtails

Sexual selection results in sex-specific characters like the conspicuously pigmented extension of the ventral tip of the caudal fin-the "sword"-in males of several species of Xiphophorus fishes. To uncover the genetic architecture underlying sword formation and to identify genes that are associated with its development, we characterized the sword transcriptional profile and combined it with genetic mapping approaches. Results showed that the male ornament of swordtails develops from a sexually non-dimorphic prepattern of transcription factors in the caudal fin. Among genes that constitute the exclusive sword transcriptome and are located in the genomic region associated with this trait we identify the potassium channel, Kcnh8, as a sword development gene. In addition to its neural function kcnh8 performs a known role in fin growth. These findings indicate that during evolution of swordtails a brain gene has been co-opted for an additional novel function in establishing a male ornament.

Melanocortin 4 receptor signaling and puberty onset regulation in Xiphophorus swordtails

Fish of the genus Xiphophorus provide a prominent example of genetic control of male body size and reproductive tactics. In X.nigrensis and X.multilineatus, puberty onset and body length are determined by melanocortin 4 receptor (Mc4r) allelic and copy number variations which were proposed to fine-tune the signaling output of the system. Accessory protein Mrap2 is required for growth across species by affecting Mc4r signaling. The molecular mechanism how Mc4r signaling controls puberty regulation in Xiphophorus and whether the interaction with Mrap2 is also involved was so far unclear. Hence, we examined Mc4r and Mrap2 in X.nigrensis and X.multilineatus, in comparison to a more distantly related species, X.hellerii. mc4r and mrap2 transcripts co-localized in the hypothalamus and preoptic regions in large males, small males and females of X.nigrensis, with similar signal strength for mrap2 but higher expression of mc4r in large males. This overexpression is constituted by wild-type and one subtype of mutant alleles. In vitro studies revealed that Mrap2 co-expressed with Mc4r increased cAMP production but did not change EC50. Cells co-expressing the wild-type and one mutant allele showed lower cAMP signaling than Mc4r wild-type cells. This indicates a role of Mc4r alleles, but not Mrap2, in puberty signaling. Different from X.nigrensis and X.multilineatus, X.hellerii has only wild-type alleles, but also shows a puberty onset and body length polymorphism, despite the absence of mutant alleles. Like in the two other species, mc4r and mrap2 transcripts colocalized and mc4r is expressed at substantially higher levels in large males. This demonstrates that puberty and growth regulation mechanism may not be identical even within same genus.

Is population subdivision different from speciation? From phylogeography to species delimitation

Species-level diversity and the underlying mechanisms that lead to the formation of new species, that is, speciation, have often been confounded with intraspecific diversity and population subdivision. The delineation between intraspecific and interspecific divergence processes has received much less attention than species delimitation. The ramifications of confounding speciation and population subdivision are that the term speciation has been used to describe many different biological divergence processes, rendering the results, or inferences, between studies incomparable. Phylogeographic studies have advanced our understanding of how spatial variation in the pattern of biodiversity can begin, become structured, and persist through time. Studies of species delimitation have further provided statistical and model-based approaches to determine the phylogeographic entities that merit species status. However, without a proper understanding and delineation between the processes that generate and maintain intraspecific and interspecific diversity in a study system, the delimitation of species may still not be biologically and evolutionarily relevant. I argue that variation in the continuity of the divergence process among biological systems could be a key factor leading to the enduring contention in delineating divergence patterns, or species delimitation, meriting future comparative studies to help us better understand the nature of biological species.

Ichthyological Differentiation and Homogenization in the Pánuco Basin, Mexico

Species introductions and extirpations are key aspects of aquatic ecosystem change that need to be examined at large geographic and temporal scales. The Pánuco Basin (Eastern Mexico) has high ichthyological diversity and ecological heterogeneity. However, freshwater fish (FWF) introductions and extirpations since the mid-1900s have modified species range and distribution. We examine changes in FWF species composition in and among four sub-basins of the Pánuco by comparing fish collection records pre-1980 to 2018. Currently, the FWF of the Pánuco includes 95 species. Fishes in the Poeciliidae, Cyprinidae, and Cichlidae, respectively, comprised most records over time. Significant differences in species composition were found between the first (pre-1980) and last (2011–2018) study periods, but not for periods in-between. Eight independent species groups were key for explaining changes in Pánuco river ichthyofauna; one group was dominated by invasive species, and saw increases in the number of records across study periods (faunal homogenization). Another group was formed by species with conservation concern with a declining number of records over time. Thirteen (2 native and 11 non-native) species were responsible for temporal turnover. These results strongly suggest high rates of differentiation over time (via native species loss) following widespread non-native species introductions.

Independent Origin of XY and ZW Sex Determination Mechanisms in Mosquitofish Sister Species

Fish are known for the outstanding variety of their sex determination mechanisms and sex chromosome systems. The western (Gambusia affinis) and eastern mosquitofish (G. holbrooki) are sister species for which different sex determination mechanisms have been described: ZZ/ZW for G. affinis and XX/XY for G. holbrooki Here, we carried out restriction-site associated DNA (RAD-) and pool sequencing (Pool-seq) to characterize the sex chromosomes of both species. We found that the ZW chromosomes of G. affinis females and the XY chromosomes of G. holbrooki males correspond to different linkage groups, and thus evolved independently from separate autosomes. In interspecific hybrids, the Y chromosome is dominant over the W chromosome, and X is dominant over Z. In G. holbrooki, we identified a candidate region for the Y-linked melanic pigmentation locus, a rare male phenotype that constitutes a potentially sexually antagonistic trait and is associated with other such characteristics, e.g., large body size and aggressive behavior. We developed a SNP-based marker in the Y-linked allele of GIPC PDZ domain containing family member 1 (gipc1), which was linked to melanism in all tested G. holbrooki populations. This locus represents an example for a color locus that is located in close proximity to a putative sex determiner, and most likely substantially contributed to the evolution of the Y.

Intergeneric hybrids inform reproductive isolating barriers in the Antarctic icefish radiation

Interspecific hybridization or barriers to hybridization may have contributed to the diversification of Antarctic icefishes (Channichthyidae), but data supporting these hypotheses is scarce. To understand the potential for hybridization and to investigate reproductive isolating mechanisms among icefish species, we performed in vitro fertilization experiments using eggs from a female blackfin icefish Chaenocephalus aceratus and sperm from a male of another genera, the ocellated icefish Chionodraco rastrospinosus. Sequencing of genomic and mitochondrial DNA confirmed the intergeneric hybrid nature of resulting embryos which successfully developed and hatched as active larvae at about four and a half months during the Antarctic winter. This result demonstrates the compatibility of gametes of these two species and the viability of resulting zygotes and larvae. Due to logistic constraints and the slow developmental rate of icefishes, we could not test for long-term hybrid viability, fertility, fitness, or hybrid breakdown. Analysis of our fishing records and available literature, however, suggests that the strongest barriers to hybridization among parapatric icefish species are likely to be behavioral and characterized by assortative mating and species-specific courtship and nesting behaviors. This conclusion suggests that, in long-lived fish species with late sexual maturity and high energetic investment in reproduction like icefishes, pre-mating barriers are energetically more efficient than post-mating barriers to prevent hybridization.

Standing geographic variation in eclosion time and the genomics of host race formation in Rhagoletis pomonella fruit flies

Taxa harboring high levels of standing variation may be more likely to adapt to rapid environmental shifts and experience ecological speciation. Here, we characterize geographic and host-related differentiation for 10,241 single nucleotide polymorphisms in Rhagoletis pomonella fruit flies to infer whether standing genetic variation in adult eclosion time in the ancestral hawthorn (Crataegus spp.)-infesting host race, as opposed to new mutations, contributed substantially to its recent shift to earlier fruiting apple (Malus domestica). Allele frequency differences associated with early vs. late eclosion time within each host race were significantly related to geographic genetic variation and host race differentiation across four sites, arrayed from north to south along a 430-km transect, where the host races co-occur in sympatry in the Midwest United States. Host fruiting phenology is clinal, with both apple and hawthorn trees fruiting earlier in the North and later in the South. Thus, we expected alleles associated with earlier eclosion to be at higher frequencies in northern populations. This pattern was observed in the hawthorn race across all four populations; however, allele frequency patterns in the apple race were more complex. Despite the generally earlier eclosion timing of apple flies and corresponding apple fruiting phenology, alleles on chromosomes 2 and 3 associated with earlier emergence were paradoxically at lower frequency in the apple than hawthorn host race across all four sympatric sites. However, loci on chromosome 1 did show higher frequencies of early eclosion-associated alleles in the apple than hawthorn host race at the two southern sites, potentially accounting for their earlier eclosion phenotype. Thus, although extensive clinal genetic variation in the ancestral hawthorn race exists and contributed to the host shift to apple, further study is needed to resolve details of how this standing variation was selected to generate earlier eclosing apple fly populations in the North.

Long-term experimental hybridisation results in the evolution of a new sex chromosome in swordtail fish

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.

Fish have a high diversity of sex-determining systems, but the mechanisms responsible for this are not well understood. Here, Franchini et al. show how hybridization and backcrossing have led to the evolution of a new sex chromosome in swordtail fish during 30 years of experimental evolution.

In search of the Goldilocks zone for hybrid speciation

Hybridization has recently gained considerable interest both as a unique opportunity for observing speciation mechanisms and as a potential engine for speciation. The latter remains a controversial topic. It was recently hypothesized that the reciprocal sorting of genetic incompatibilities from parental species could result in hybrid speciation, when the hybrid population maintains a mixed combination of the parental incompatibilities that prevents further gene exchange with both parental populations. However, the specifics of the purging/sorting process of multiple incompatibilities have not been examined theoretically. We here investigate the allele-frequency dynamics of an isolated hybrid population that results from a single hybridization event. Using models of two or four loci, we investigate the fate of one or two genetic incompatibilities of the Dobzhansky-Muller type (DMIs). We study how various parameters affect both the sorting/purging of the DMIs and the probability of observing hybrid speciation by reciprocal sorting. We find that the probability of hybrid speciation is strongly dependent on the linkage architecture (i.e. the order and recombination rate between loci along chromosomes), the population size of the hybrid population, and the initial relative contributions of the parental populations to the hybrid population. We identify a Goldilocks zone for specific linkage architectures and intermediate recombination rates, in which hybrid speciation becomes highly probable. Whereas an equal contribution of parental populations to the hybrid population maximizes the hybrid speciation probability in the Goldilocks zone, other linkage architectures yield unintuitive asymmetric maxima. We provide an explanation for this pattern, and discuss our results both with respect to the best conditions for observing hybrid speciation in nature and their implications regarding patterns of introgression in hybrid zones.

Hybridization is observed ubiquitously in nature. Its outcome can range from extinction to the creation of new species. With respect to the latter, the probability of homoploid hybrid speciation, i.e. the formation of a new species as a result of a hybridization event without changes in the ploidy of the organism, is a hotly debated topic. Here, we analyze a minimal model for homoploid hybrid speciation, in which reproductive isolation is achieved by means of (postzygotic) Dobzhansky-Muller incompatibilities. When these postzygotic genetic incompatibilities are resolved in the hybrid population, their reciprocal sorting can result in reproductive isolation from both parental populations, thus creating a hybrid species. We show that, in accordance with the current literature, hybrid speciation tends to be rare. However, specific arrangements of the genes responsible for reproductive isolation can make reciprocal sorting almost unavoidable and thus create barriers to the parental population in an almost deterministic matter. We discuss the implications of these results for hybrid speciation and patterns of introgression in nature.

Female Choice Undermines the Emergence of Strong Sexual Isolation between Locally Adapted Populations of Atlantic Mollies (Poecilia mexicana)

Divergent selection between ecologically dissimilar habitats promotes local adaptation, which can lead to reproductive isolation (RI). Populations in the Poecilia mexicana species complex have independently adapted to toxic hydrogen sulfide and show varying degrees of RI. Here, we examined the variation in the mate choice component of prezygotic RI. Mate choice tests across drainages (with stimulus males from another drainage) suggest that specific features of the males coupled with a general female preference for yellow color patterns explain the observed variation. Analyses of male body coloration identified the intensity of yellow fin coloration as a strong candidate to explain this pattern, and common-garden rearing suggested heritable population differences. Male sexual ornamentation apparently evolved differently across sulfide-adapted populations, for example because of differences in natural counterselection via predation. The ubiquitous preference for yellow color ornaments in poeciliid females likely undermines the emergence of strong RI, as female discrimination in favor of own males becomes weaker when yellow fin coloration in the respective sulfide ecotype increases. Our study illustrates the complexity of the (partly non-parallel) pathways to divergence among replicated ecological gradients. We suggest that future work should identify the genomic loci involved in the pattern reported here, making use of the increasing genomic and transcriptomic datasets available for our study system.

Rapid hybrid speciation in Darwin’s finches

Galapagos finches have driven hypotheses of how speciation occurs. Most commonly, it is assumed that natural selection separates species originating from a single population on the basis of variation in traits that confer advantages for survival and reproduction. Lamichhaney et al. document a case where cross-species hybridization established a reproductively isolated lineage, which demonstrates a process known as homoploid hybrid speciation in action (see the Perspective by Wagner). The authors used genetic markers and phenotypic analyses to create a pedigree that revealed how a cross-island migrant bred with a native species to form a self-perpetuating hybrid population that was reproductively isolated from both parental species.

Science , this issue p. 224 ; see also p. 157

Characterization and phylogenetic analysis of complete mitochondrial genomes for two desert cyprinodontoid fishes, Empetrichthys latos and Crenichthys baileyi

The Pahrump poolfish (Empetrichthys latos) and White River springfish (Crenichthys baileyi) are small-bodied teleost fishes (order Cyprinodontiformes) endemic to the arid Great Basin and Mojave Desert regions of western North America. These taxa survive as small, isolated populations in remote streams and springs and evolved to tolerate extreme conditions of high temperature and low dissolved oxygen. Both species have experienced severe population declines over the last 50-60years that led to some subspecies being categorized with protected status under the U.S. Endangered Species Act. Here we report the first sequencing of the complete mitochondrial DNA genomes for both E. l. latos and the moapae subspecies of C. baileyi. Complete mitogenomes of 16,546bp nucleotides were obtained from two E. l. latos individuals collected from introduced populations at Spring Mountain Ranch State Park and Shoshone Ponds Natural Area, Nevada, USA, while a single mitogenome of 16,537bp was sequenced for C. b. moapae. The mitogenomes of both species contain 13 protein-encoding genes, twenty-two tRNAs, and two rRNAs (12S and 18S) following the syntenic arrangement typical of Actinopterygiian fish mitogenomes, as well as D-loop control regions of 858bp for E. latos and 842bp for C. baileyi moapae. The two E. latos individuals exhibited only 0.0181% nucleotide sequence divergence across the entire mitogenome, implying little intraspecific mtDNA genetic variation. Comparative phylogenetic analysis of the poolfish and springfish mitochondrial genomes to available mitogenomes of other Cyprinodontoid fishes confirmed the close relationship of these oviparous Empetrichthys and Crenichthys genera to the viviparous goodeid fishes of central Mexico, and showed the combined clade of these fishes to be a sister group to the Profundulidae killifishes. Despite several significant life history and morphological differences between the Empetrichthyinae and Goodienae, estimates of evolutionary genetic distances using two partial regions of mtDNA point to inclusion of the Empetrichthys and Crenichthys genera within the family Goodeidae along with the goodeid fishes of central Mexico.

Characterization of appetite-regulating factors in platyfish, Xiphophorus maculatus (Cyprinodontiformes Poeciliidae)

The regulation of energy in fish, like most vertebrates, is a complex process that involves a number of brain and peripheral hormones. These signals include anorexigenic (e.g. cholecystokinin (CCK) and cocaine- and amphetamine-regulated transcript (CART)) as well as orexigenic (e.g. orexin and neuropeptide Y (NPY)) peptides. Platyfish, Xiphophorus maculatus, are freshwater viviparous fish for which little is known about the endocrine mechanisms regulating feeding. In order to elucidate the role of these peptides in the regulation of feeding of platyfish, we examined the effects of peripheral injections of CCK and orexin on feeding behavior and food intake. Injections of CCK decreased both food intake and searching behavior, while injections of orexin increased searching behavior but did not affect food consumption. In order to better characterize these peptides, we examined their mRNA tissue distribution and assessed the effects of a 10-day fast on their brain and intestine expressions in both males and females. CCK, CART, NPY and orexin all show widespread distributions in brain and several peripheral tissues, including intestine and gonads. Fasting induced decreases in both CCK and CART and an increase in orexin mRNA expressions in the brain and a decrease in CCK expression in the intestine, but did not affect either expressions of NPY. There were no significant sex-specific differences in either the behavioral responses to injections or the expression responses to fasting. The widespread distribution and the fasting-induced changes in expression of these peptides suggest that they might have several physiological roles in platyfish, including the regulation of feeding.

Do relaxed selection and habitat temperature facilitate biased mitogenomic introgression in a narrowly endemic fish?

Introgression might be exceptionally common during the evolution of narrowly endemic species. For instance, in the springs of the small and isolated Cuatro Ciénegas Valley, the mitogenome of the cichlid fish Herichthys cyanoguttatus could be rapidly introgressing into populations of the trophically polymorphic H. minckleyi. We used a combination of genetic and environmental data to examine the factors associated with this mitochondrial introgression. A reduced representation library of over 6220 single nucleotide polymorphisms (SNPs) from the nuclear genome showed that mitochondrial introgression into H. minckleyi is biased relative to the amount of nuclear introgression. SNP assignment probabilities also indicated that cichlids with more hybrid ancestry are not more commonly female providing no support for asymmetric backcrossing or hybrid‐induced sex‐ratio distortion in generating the bias in mitochondrial introgression. Smaller effective population size in H. minckleyi inferred from the SNPs coupled with sequences of all 13 mitochondrial proteins suggests that relaxed selection on the mitogenome could be facilitating the introgression of “H. cyanoguttatus” haplotypes. Additionally, we showed that springs with colder temperatures had greater amounts of mitochondrial introgression from H. cyanoguttatus. Relaxed selection in H. minckleyi coupled with temperature‐related molecular adaptation could be facilitating mitogenomic introgression into H. minckleyi.

Ecological divergence and conservatism: spatiotemporal patterns of niche evolution in a genus of livebearing fishes (Poeciliidae: Xiphophorus)

Background

Ecological factors often have a strong impact on spatiotemporal patterns of biodiversity. The integration of spatial ecology and phylogenetics allows for rigorous tests of whether speciation is associated with niche conservatism (constraints on ecological divergence) or niche divergence. We address this question in a genus of livebearing fishes for which the role of sexual selection in speciation has long been studied, but in which the potential role of ecological divergence during speciation has not been tested.

Results

By combining reconstruction of ancestral climate tolerances and disparity indices, we show that the earliest evolutionary split in Xiphophorus was associated with significant divergence for temperature variables. Niche evolution and present day niches were most closely associated with each species’ geographic distribution relative to a biogeographic barrier, the Trans-Mexican Volcanic Belt. Tests for similarity of the environmental backgrounds of closely related species suggested that the relative importance of niche conservatism and divergence during speciation varied among the primary clades of Xiphophorus. Closely related species in the two swordtail clades exhibited higher levels of niche overlap than expected given environmental background similarity indicative of niche conservatism. In contrast, almost all species of platyfish had significantly divergent niches compared to environmental backgrounds, which is indicative of niche divergence.

Conclusion

The results suggest that the relative importance of niche conservatism and divergence differed among the clades of Xiphophorus and that traits associated with niche evolution may be more evolutionarily labile in the platyfishes. Our results ultimately suggest that the taxonomic scale of tests for conservatism and divergence could greatly influence inferences of their relative importance in the speciation process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0593-4) contains supplementary material, which is available to authorized users.

X. couchianus and X. hellerii genome models provide genomic variation insight among Xiphophorus species

BACKGROUND

Xiphophorus fishes are represented by 26 live-bearing species of tropical fish that express many attributes (e.g., viviparity, genetic and phenotypic variation, ecological adaptation, varied sexual developmental mechanisms, ability to produce fertile interspecies hybrids) that have made attractive research models for over 85 years. Use of various interspecies hybrids to investigate the genetics underlying spontaneous and induced tumorigenesis has resulted in the development and maintenance of pedigreed Xiphophorus lines specifically bred for research. The recent availability of the X. maculatus reference genome assembly now provides unprecedented opportunities for novel and exciting comparative research studies among Xiphophorus species.

RESULTS

We present sequencing, assembly and annotation of two new genomes representing Xiphophorus couchianus and Xiphophorus hellerii. The final X. couchianus and X. hellerii assemblies have total sizes of 708 Mb and 734 Mb and correspond to 98 % and 102 % of the X. maculatus Jp 163 A genome size, respectively. The rates of single nucleotide change range from 1 per 52 bp to 1 per 69 bp among the three genomes and the impact of putatively damaging variants are presented. In addition, a survey of transposable elements allowed us to deduce an ancestral TE landscape, uncovered potential active TEs and document a recent burst of TEs during evolution of this genus.

CONCLUSIONS

Two new Xiphophorus genomes and their corresponding transcriptomes were efficiently assembled, the former using a novel guided assembly approach. Three assembled genome sequences within this single vertebrate order of new world live-bearing fishes will accelerate our understanding of relationship between environmental adaptation and genome evolution. In addition, these genome resources provide capability to determine allele specific gene regulation among interspecies hybrids produced by crossing any of the three species that are known to produce progeny predisposed to tumor development.

Genes as cues: phenotypic integration of genetic and epigenetic information from a Darwinian perspective

The development of multicellular organisms involves a delicate interplay between genetic and environmental influences. It is often useful to think of developmental systems as integrating available sources of information about current conditions to produce organisms. Genes and inherited physiology provide cues, as does the state of the environment during development. The integration systems themselves are under genetic control and subject to Darwinian selection, so we expect them to evolve to produce organisms that fit well with current ecological (including social) conditions. We argue for the scientific value of this explicitly informational perspective by providing detailed examples of how it can elucidate taxonomically diverse phenomena. We also present a general framework for linking genetic and phenotypic variation from an informational perspective. This application of Darwinian logic at the organismal level can elucidate genetic influences on phenotypic variation in novel and counterintuitive ways.

Automatic selection of partitioning schemes for phylogenetic analyses using iterative k-means clustering of site rates

BACKGROUND

Model selection is a vital part of most phylogenetic analyses, and accounting for the heterogeneity in evolutionary patterns across sites is particularly important. Mixture models and partitioning are commonly used to account for this variation, and partitioning is the most popular approach. Most current partitioning methods require some a priori partitioning scheme to be defined, typically guided by known structural features of the sequences, such as gene boundaries or codon positions. Recent evidence suggests that these a priori boundaries often fail to adequately account for variation in rates and patterns of evolution among sites. Furthermore, new phylogenomic datasets such as those assembled from ultra-conserved elements lack obvious structural features on which to define a priori partitioning schemes. The upshot is that, for many phylogenetic datasets, partitioned models of molecular evolution may be inadequate, thus limiting the accuracy of downstream phylogenetic analyses.

RESULTS

We present a new algorithm that automatically selects a partitioning scheme via the iterative division of the alignment into subsets of similar sites based on their rates of evolution. We compare this method to existing approaches using a wide range of empirical datasets, and show that it consistently leads to large increases in the fit of partitioned models of molecular evolution when measured using AICc and BIC scores. In doing so, we demonstrate that some related approaches to solving this problem may have been associated with a small but important bias.

CONCLUSIONS

Our method provides an alternative to traditional approaches to partitioning, such as dividing alignments by gene and codon position. Because our method is data-driven, it can be used to estimate partitioned models for all types of alignments, including those that are not amenable to traditional approaches to partitioning.

The effects of partitioning on phylogenetic inference

Partitioning is a commonly used method in phylogenetics that aims to accommodate variation in substitution patterns among sites. Despite its popularity, there have been few systematic studies of its effects on phylogenetic inference, and there have been no studies that compare the effects of different approaches to partitioning across many empirical data sets. In this study, we applied four commonly used approaches to partitioning to each of 34 empirical data sets, and then compared the resulting tree topologies, branch-lengths, and bootstrap support estimated using each approach. We find that the choice of partitioning scheme often affects tree topology, particularly when partitioning is omitted. Most notably, we find occasional instances where the use of a suboptimal partitioning scheme produces highly supported but incorrect nodes in the tree. Branch-lengths and bootstrap support are also affected by the choice of partitioning scheme, sometimes dramatically so. We discuss the reasons for these effects and make some suggestions for best practice.

Transcriptomics of two evolutionary novelties: how to make a sperm-transfer organ out of an anal fin and a sexually selected "sword" out of a caudal fin

Swords are exaggerated male ornaments of swordtail fishes that have been of great interest to evolutionary biologists ever since Darwin described them in the Descent of Man (1871). They are a novel sexually selected trait derived from modified ventral caudal fin rays and are only found in the genus Xiphophorus. Another phylogenetically more widespread and older male trait is the gonopodium, an intromittent organ found in all poeciliid fishes, that is derived from a modified anal fin. Despite many evolutionary and behavioral studies on both traits, little is known so far about the molecular mechanisms underlying their development. By investigating transcriptomic changes (utilizing a RNA-Seq approach) in response to testosterone treatment in the swordtail fish, Xiphophorus hellerii, we aimed to better understand the architecture of the gene regulatory networks underpinning the development of these two evolutionary novelties. Large numbers of genes with tissue-specific expression patterns were identified. Among the "sword genes" those involved in embryonic organ development, sexual character development and coloration were highly expressed, while in the gonopodium rather more morphogenesis-related genes were found. Interestingly, many genes and genetic pathways are shared between both developing novel traits derived from median fins: the sword and the gonopodium. Our analyses show that a larger set of gene networks was co-opted during the development and evolution of the "older" gonopodium than in the "younger," and morphologically less complex trait, the sword. We provide a catalog of candidate genes for future efforts to dissect the development of those sexually selected exaggerated male traits in swordtails.

Interspecific variation in the repair of UV damaged DNA in the genus Xiphophorus as a factor in the decline of the Rio Grande Platyfish

The fish genus Xiphophorus consists of 26 species distributed along the eastern slopes of mountain ranges extending from northern Mexico to Belize and Nicaragua. We analyzed light-dependent repair of UV-induced DNA damage in at least two species from each of the four monophyletic Xiphophorus groups. We found that the northern platyfish had significantly reduced photoenzymatic repair compared to the other three groups, including the northern swordtails, southern platyfish and southern swordtails. All of the species of the northern platyfish, including the Marbled (meyeri), Northern (gordoni) and Monterrey Platyfish (couchianus) are the northernmost species in the genus and are the only three species in the genus that are currently found on the IUCN Red List of Threatened Species. Satellite data from the past 30 years (1979-2008) correlate greater increases in shorter wavelength UVB with higher latitudes within the Xiphophorus range. We suggest that, combined with other consequences of human population growth, anthropogenic deozonation resulting in a disproportionate increase in UVB in temperate latitudes may be a contributing factor in the decline and extirpation of the northern platyfish.

Characterization and differential expression of CPD and 6-4 DNA photolyases in Xiphophorus species and interspecies hybrids

Among the many Xiphophorus interspecies hybrid tumor models are those that exhibit ultraviolet light (UVB) induced melanoma. In previous studies, assessment of UVB induced DNA damage and nucleotide excision DNA repair has been performed in parental lines and interspecies hybrids. Species and hybrid specific differences in the levels of DNA damage induced and the dark repair rates for cyclobutane pyrimidine dimers (CPDs) and 6-4 pyrimidine pyrimidine photoproducts (6-4PPs) have been reported. However, UVB induced DNA lesions in Xiphophorus fishes are thought to primarily be repaired via light dependent CPD and 6-4PP specific photolyases. Photolyases are of evolutionary interest since they are ancient and presumably function solely to ameliorate the deleterious effects of UVB exposure. Herein, we report results from detailed studies of CPD and 6-4PP photolyase gene expression within several Xiphophorus tissues. We determined photolyase gene expression patterns before and after exposure to fluorescent light in X. maculatus, X. couchianus, and for F1 interspecies hybrids produced from crossing these two parental lines (X. maculatus Jp 163 B×X. couchianus). We present novel results showing these two photolyase genes exhibit species, tissue, and hybrid-specific differences in basal and light induced gene expression.

The evolution of the sexually selected sword in Xiphophorus does not compromise aerobic locomotor performance

Sexual selection can increase morphological diversity within and among species. Little is known regarding how interspecific variation produced through sexual selection affects other functional systems. Here, we examine how morphological diversity resulting from sexual selection impacts aerobic locomotor performance. Using Xiphophorus (swordtail fish) and their close relatives (N = 19 species), we examined whether the evolution of a longer sexually selected sword affects critical swimming speed. We also examined the effect of other suborganismal, physiological, and morphological traits on critical swimming speed, as well as their relationship with sword length. In correlation analyses, we found no significant relationship between sword length and critical swimming speed. Unexpectedly, we found that critical swimming speed was higher in species with longer swords, after controlling for body size in multiple regression analyses. We also found several suborganismal and morphological predictors of critical swimming speed, as well as a significant negative relationship between sword length and heart and gill mass. Our results suggest that interspecific variation in sword length is not costly for this aspect of swimming performance, but further studies should examine potential costs for other types of locomotion and other components of Darwinian fitness (e.g., survivorship, life span).