Dynamic formation of asexual diploid and polyploid lineages: multilocus analysis of Cobitis reveals the mechanisms maintaining the diversity of clones

Genetic and karyotype divergence between parents affect clonality and sterility in hybrids

Asexual reproduction can be triggered by interspecific hybridization, but its emergence is supposedly rare, relying on exceptional combinations of suitable genomes. To examine how genomic and karyotype divergence between parental lineages affect the incidence of asexual gametogenesis, we experimentally hybridized fishes (Cobitidae) across a broad phylogenetic spectrum, assessed by whole exome data. Gametogenic pathways generally followed a continuum from sexual reproduction in hybrids between closely related evolutionary lineages to sterile or inviable crosses between distant lineages. However, most crosses resulted in a combination of sterile males and asexually reproducing females. Their gametes usually experienced problems in chromosome pairing, but females also produced a certain proportion of oocytes with premeiotically duplicated genomes, enabling their development into clonal eggs. Interspecific hybridization may thus commonly affect cell cycles in a specific way, allowing the formation of unreduced oocytes. The emergence of asexual gametogenesis appears tightly linked to hybrid sterility and constitutes an inherent part of the extended speciation continuum.

Why do hybrids turn down sex?

Asexual reproduction is ancestral in prokaryotes; the switch to sexuality in eukaryotes is one of the major transitions in the history of life. The study of the maintenance of sex in eukaryotes has raised considerable interest for decades and is still one of evolutionary biology's most prominent question. The observation that many asexual species are of hybrid origin has led some to propose that asexuality in hybrids results from sexual processes being disturbed because of incompatibilities between the two parental species' genomes. However, in some cases, failure to produce asexual F1s in the lab may indicate that this mechanism is not the only road to asexuality in hybrid species. Here, we present a mathematical model and propose an alternative, adaptive route for the evolution of asexuality from previously sexual hybrids. Under some reproductive alterations, we show that asexuality can evolve to rescue hybrids' reproduction. Importantly, we highlight that when incompatibilities only affect the fusion of sperm and egg's genomes, the two traits that characterize asexuality, namely unreduced meiosis and the initiation of embryogenesis without the incorporation of the sperm's pronucleus, can evolve separately, greatly facilitating the overall evolutionary route. Taken together, our results provide an alternative, potentially complementary explanation for the link between asexuality and hybridization.

Sperm-dependent asexual species and their role in ecology and evolution

Sexual reproduction is the primary mode of reproduction in eukaryotes, but some organisms have evolved deviations from classical sex and switched to asexuality. These asexual lineages have sometimes been viewed as evolutionary dead ends, but recent research has revealed their importance in many areas of general biology. Our review explores the understudied, yet important mechanisms by which sperm-dependent asexuals that produce non-recombined gametes but rely on their fertilization, can have a significant impact on the evolution of coexisting sexual species and ecosystems. These impacts are concentrated around three major fields. Firstly, sperm-dependent asexuals can potentially impact the gene pool of coexisting sexual species by either restricting their population sizes or by providing bridges for interspecific gene flow whose type and consequences substantially differ from gene flow mechanisms expected under sexual reproduction. Secondly, they may impact on sexuals' diversification rates either directly, by serving as stepping-stones in speciation, or indirectly, by promoting the formation of pre- and postzygotic reproduction barriers among nascent species. Thirdly, they can potentially impact on spatial distribution of species, via direct or indirect (apparent) types of competition and Allee effects. For each such mechanism, we provide empirical examples of how natural sperm-dependent asexuals impact the evolution of their sexual counterparts. In particular, we highlight that these broad effects may last beyond the tenure of the individual asexual lineages causing them, which challenges the traditional perception that asexual lineages are short-lived evolutionary dead ends and minor sideshows. Our review also proposes new research directions to incorporate the aforementioned impacts of sperm-dependent asexuals. These research directions will ultimately enhance our understanding of the evolution of genomes and biological interactions in general.

Evidence for cryptic sex in parthenogenetic stick insects of the genus Timema

Obligately parthenogenetic species are expected to be short lived since the lack of sex and recombination should translate into a slower adaptation rate and increased accumulation of deleterious alleles. Some, however, are thought to have been reproducing without males for millions of years. It is not clear how these old parthenogens can escape the predicted long-term costs of parthenogenesis, but an obvious explanation is cryptic sex. In this study, we screen for signatures of cryptic sex in eight populations of four parthenogenetic species of Timema stick insects, some estimated to be older than 1 Myr. Low genotype diversity, homozygosity of individuals and high linkage disequilibrium (LD) unaffected by marker distances support exclusively parthenogenetic reproduction in six populations. However, in two populations (namely, of the species Timema douglasi and T. monikensis) we find strong evidence for cryptic sex, most likely mediated by rare males. These populations had comparatively high genotype diversities, lower LD, and a clear LD decay with genetic distance. Rare sex in species that are otherwise largely parthenogenetic could help explain the unusual success of parthenogenesis in the Timema genus and raises the question whether episodes of rare sex are in fact the simplest explanation for the persistence of many old parthenogens in nature.

Sequence capture: Obsolete or irreplaceable? A thorough validation across phylogenetic distances and its applicability to hybrids and allopolyploids

As whole-genome sequencing has become pervasive, some have suggested that reduced genomic representation approaches, for example, sequence capture, are becoming obsolete. In the present study, we argue that these techniques still provide excellent tools in terms of price and quality of data as well as in their ability to provide markers with specific features, as required, for example, in phylogenomics. A potential drawback of the wide-scale application of reduced representation approaches could be their drop in efficiency with increasing phylogenetic distance from the reference species. While some studies have focused on the degree and performance of reduced representation techniques in such situations, to our knowledge, none of them evaluated their applicability to inter-specific hybrids and polyploids. This highlights a significant gap in current knowledge since there is increasing evidence for the frequent occurrence of natural hybrids and polyploids, as well as for the major importance of both phenomena in evolution. The main aim of the present study was to carry out a thorough validation of SEQcap applicability to (1) a set of non-model taxa with a wide range of phylogenetic relatedness and (2) inter-specific hybrids of various ploidies and genomic compositions. Considering the latter point, we especially focused on mechanisms causing allelic bias and consequent allelic dropout, as these could have confounding effects with respect to the evolutionary genomic dynamics of hybrids, especially in asexuals, which virtually reproduce as a frozen F1 generation.

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.

Clonal gametogenesis is triggered by intrinsic stimuli in the hybrid's germ cells but is dependent on sex differentiation

Interspecific hybridization may trigger the transition from sexual reproduction to asexuality, but mechanistic reasons for such a change in a hybrid's reproduction are poorly understood. Gametogenesis of many asexual hybrids involves a stage of premeiotic endoreduplication (PMER), when gonial cells duplicate chromosomes and subsequent meiotic divisions involve bivalents between identical copies, leading to production of clonal gametes. Here, we investigated the triggers of PMER and whether its induction is linked to intrinsic stimuli within a hybrid's gonial cells or whether it is regulated by the surrounding gonadal tissue. We investigated gametogenesis in the Cobitis taenia hybrid complex, which involves sexually reproducing species (Cobitis elongatoides and C. taenia) as well as their hybrids, where females reproduce clonally via PMER while males are sterile. We transplanted spermatogonial stem cells (SSCs) from C. elongatoides and triploid hybrid males into embryos of sexual species and of asexual hybrid females, respectively, and observed their development in an allospecific gonadal environment. Sexual SSCs underwent regular meiosis and produced normally reduced gametes when transplanted into clonal females. On the other hand, the hybrid's SSCs lead to sterility when transplanted into sexual males, but maintained their ability to undergo asexual development (PMER) and production of clonal eggs, when transplanted into sexual females. This suggests that asexual gametogenesis is under complex control when somatic gonadal tissue indirectly affects the execution of asexual development by determining the sexual differentiation of stem cells and once such cells develop to female phenotypes, hybrid germ cells trigger the PMER from their intrinsic signals.

The Role of Interspecific Hybridisation in Adaptation and Speciation: Insights From Studies in Senecio

Hybridisation is well documented in many species, especially plants. Although hybrid populations might be short-lived and do not evolve into new lineages, hybridisaiton could lead to evolutionary novelty, promoting adaptation and speciation. The genus Senecio (Asteraceae) has been actively used to unravel the role of hybridisation in adaptation and speciation. In this article, we first briefly describe the process of hybridisation and the state of hybridisation research over the years. We then discuss various roles of hybridisation in plant adaptation and speciation illustrated with examples from different Senecio species, but also mention other groups of organisms whenever necessary. In particular, we focus on the genomic and transcriptomic consequences of hybridisation, as well as the ecological and physiological aspects from the hybrids' point of view. Overall, this article aims to showcase the roles of hybridisation in speciation and adaptation, and the research potential of Senecio, which is part of the ecologically and economically important family, Asteraceae.

Genome Fractionation and Loss of Heterozygosity in Hybrids and Polyploids: Mechanisms, Consequences for Selection, and Link to Gene Function

Hybridization and genome duplication have played crucial roles in the evolution of many animal and plant taxa. The subgenomes of parental species undergo considerable changes in hybrids and polyploids, which often selectively eliminate segments of one subgenome. However, the mechanisms underlying these changes are not well understood, particularly when the hybridization is linked with asexual reproduction that opens up unexpected evolutionary pathways. To elucidate this problem, we compared published cytogenetic and RNAseq data with exome sequences of asexual diploid and polyploid hybrids between three fish species; Cobitis elongatoides, C. taenia, and C. tanaitica. Clonal genomes remained generally static at chromosome-scale levels but their heterozygosity gradually deteriorated at the level of individual genes owing to allelic deletions and conversions. Interestingly, the impact of both processes varies among animals and genomic regions depending on ploidy level and the properties of affected genes. Namely, polyploids were more tolerant to deletions than diploid asexuals where conversions prevailed, and genomic restructuring events accumulated preferentially in genes characterized by high transcription levels and GC-content, strong purifying selection and specific functions like interacting with intracellular membranes. Although hybrids were phenotypically more similar to C. taenia, we found that they preferentially retained C. elongatoides alleles. This demonstrates that favored subgenome is not necessarily the transcriptionally dominant one. This study demonstrated that subgenomes in asexual hybrids and polyploids evolve under a complex interplay of selection and several molecular mechanisms whose efficiency depends on the organism's ploidy level, as well as functional properties and parental ancestry of the genomic region.

Challenges and Costs of Asexuality: Variation in Premeiotic Genome Duplication in Gynogenetic Hybrids from Cobitis taenia Complex

The transition from sexual reproduction to asexuality is often triggered by hybridization. The gametogenesis of many hybrid asexuals involves premeiotic genome endoreplication leading to bypass hybrid sterility and forming clonal gametes. However, it is still not clear when endoreplication occurs, how many gonial cells it affects and whether its rate differs among clonal lineages. Here, we investigated meiotic and premeiotic cells of diploid and triploid hybrids of spined loaches (Cypriniformes: Cobitis) that reproduce by gynogenesis. We found that in naturally and experimentally produced F1 hybrids asexuality is achieved by genome endoreplication, which occurs in gonocytes just before entering meiosis or, rarely, one or a few divisions before meiosis. However, genome endoreplication was observed only in a minor fraction of the hybrid's gonocytes, while the vast majority of gonocytes were unable to duplicate their genomes and consequently could not proceed beyond pachytene due to defects in bivalent formation. We also noted that the rate of endoreplication was significantly higher among gonocytes of hybrids from natural clones than of experimentally produced F1 hybrids. Thus, asexuality and hybrid sterility are intimately related phenomena and the transition from sexual reproduction to asexuality must overcome significant problems with genome incompatibilities with a possible impact on reproductive potential.

Parthenogenesis is self-destructive for scaled reptiles

Parthenogenesis is rare in nature. With 39 described true parthenogens, scaled reptiles (Squamata) are the only vertebrates that evolved this reproductive strategy. Parthenogenesis is ecologically advantageous in the short term, but the young age and rarity of parthenogenetic species indicate it is less advantageous in the long term. This suggests parthenogenesis is self-destructive: it arises often but is lost due to increased extinction rates, high rates of reversal or both. However, this role of parthenogenesis as a self-destructive trait remains unknown. We used a phylogeny of Squamata (5388 species), tree metrics, null simulations and macroevolutionary scenarios of trait diversification to address the factors that best explain the rarity of parthenogenetic species. We show that parthenogenesis can be considered as self-destructive, with high extinction rates mainly responsible for its rarity in nature. Since these parthenogenetic species occur, this trait should be ecologically relevant in the short term.

Sex and Asex: A Clonal Lexicon

Organisms across the tree of life have complex life cycles that include both sexual and asexual reproduction or that are obligately asexual. These organisms include ecologically dominant species that structure many terrestrial and marine ecosystems, as well as many pathogens, pests, and invasive species. We must consider both the evolution and maintenance of these various reproductive modes and how these modes shape the genetic diversity, adaptive evolution, and ability to persist in the species that exhibit them. Thus, having a common framework is a key aspect of understanding the biodiversity that shapes our planet. In the 2019 AGA President's Symposium, Sex and Asex: The genetics of complex life cycles, researchers investigating a wide range of taxonomic models and using a variety of modes of investigation coalesced around a common theme-understanding not only how such complex life cycles may evolve, but how they are shaped by the evolutionary and ecological forces around them. In this introduction to the Special Issue from the symposium, we give an overview of some of the key ideas and areas of investigation (a common clonal lexicon, we might say) and introduce the breadth of work submitted by symposium participants.

No evidence for accumulation of deleterious mutations and fitness degradation in clonal fish hybrids: Abandoning sex without regrets

Despite its inherent costs, sexual reproduction is ubiquitous in nature, and the mechanisms to protect it from a competitive displacement by asexuality remain unclear. Popular mutation-based explanations, like the Muller's ratchet and the Kondrashov's hatchet, assume that purifying selection may not halt the accumulation of deleterious mutations in the nonrecombining genomes, ultimately leading to their degeneration. However, empirical evidence is scarce and it remains particularly unclear whether mutational degradation proceeds fast enough to ensure the decay of clonal organisms and to prevent them from outcompeting their sexual counterparts. To test this hypothesis, we jointly analysed the exome sequences and the fitness-related phenotypic traits of the sexually reproducing fish species and their clonal hybrids, whose evolutionary ages ranged from F1 generations to 300 ky. As expected, mutations tended to accumulate in the clonal genomes in a time-dependent manner. However, contrary to the predictions, we found no trend towards increased nonsynonymity of mutations acquired by clones, nor higher radicality of their amino acid substitutions. Moreover, there was no evidence for fitness degeneration in the old clones compared with that in the younger ones. In summary, although an efficacy of purifying selection may still be reduced in the asexual genomes, our data indicate that its efficiency is not drastically decreased. Even the oldest investigated clone was found to be too young to suffer fitness consequences from a mutation accumulation. This suggests that mechanisms other than mutation accumulation may be needed to explain the competitive advantage of sex in the short term.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Cytogenetic Characterization of Seven Novel satDNA Markers in Two Species of Spined Loaches (Cobitis) and Their Clonal Hybrids

Interspecific hybridization is a powerful evolutionary force. However, the investigation of hybrids requires the application of methodologies that provide efficient and indubitable identification of both parental subgenomes in hybrid individuals. Repetitive DNA, and especially the satellite DNA sequences (satDNA), can rapidly diverge even between closely related species, hence providing a useful tool for cytogenetic investigations of hybrids. Recent progress in whole-genome sequencing (WGS) offers unprecedented possibilities for the development of new tools for species determination, including identification of species-specific satDNA markers. In this study, we focused on spined loaches (Cobitis, Teleostei), a group of fishes with frequent interspecific hybridization. Using the WGS of one species, C. elongatoides, we identified seven satDNA markers, which were mapped by fluorescence in situ hybridization on mitotic and lampbrush chromosomes of C. elongatoides, C. taenia and their triploid hybrids (C. elongatoides × 2C. taenia). Two of these markers were chromosome-specific in both species, one had centromeric localization in multiple chromosomes and four had variable patterns between tested species. Our study provided a novel set of cytogenetic markers for Cobitis species and demonstrated that NGS-based development of satDNA cytogenetic markers may provide a very efficient and easy tool for the investigation of hybrid genomes, cell ploidy, and karyotype evolution.

Variable occurrence of apoptosis in the testes of diploid and sterile allotetraploid Cobitis (Teleostei, Cobitidae) males during the reproductive cycle

Cobitis species exist in both diploid and diploid-polyploid (d-p) populations, but mostly occur in the latter. They are considered an important model organism to study the biology and physiology of natural hybrid and polyploid vertebrates. Indeed, polyploidization causes a huge stress for in terms of cell physiology and alters spermatogenesis in polyploid fish. The most extensively studied mode of germ cell death during spermatogenesis in vertebrates is apoptosis. The aim of this study was to examine caspase-3 immunoexpression in the testes of Cobitis taenia from a diploid population as well as C. taenia and sterile tetraploid Cobitis from d-p populations before, during and after spawning. The obtained results suggest a different performance of apoptosis in the testes of C. taenia from the two studied populations and seems to be conditioned by their role as the only sperm donors in d-p populations. Moreover, apoptosis was an active cell death process in the testes of tetraploid Cobitis.

The Legacy of Sexual Ancestors in Phenotypic Variability, Gene Expression, and Homoeolog Regulation of Asexual Hybrids and Polyploids

Hybridization and polyploidization are important evolutionary processes whose impacts range from the alteration of gene expression and phenotypic variation to the triggering of asexual reproduction. We investigated fishes of the Cobitis taenia-elongatoides hybrid complex, which allowed us to disentangle the direct effects of both processes, due to the co-occurrence of parental species with their diploid and triploid hybrids. Employing morphological, ecological, and RNAseq approaches, we investigated the molecular determinants of hybrid and polyploid forms. In contrast with other studies, hybridization and polyploidy induced relatively very little transgressivity. Instead, Cobitis hybrids appeared intermediate with a clear effect of genomic dosing when triploids expressed higher similarity to the parent contributing two genome sets. This dosage effect was symmetric in the germline (oocyte gene expression), interestingly though, we observed an overall bias toward C. taenia in somatic tissues and traits. At the level of individual genes, expression-level dominance vastly prevailed over additivity or transgressivity. Also, trans-regulation of gene expression was less efficient in diploid hybrids than in triploids, where the expression modulation of homoeologs derived from the "haploid" parent was stronger than those derived from the "diploid" parent. Our findings suggest that the apparent intermediacy of hybrid phenotypes results from the combination of individual genes with dominant expression rather than from simple additivity. The efficiency of cross-talk between trans-regulatory elements further appears dosage dependent. Important effects of polyploidization may thus stem from changes in relative concentrations of trans-regulatory elements and their binding sites between hybridizing genomes. Links between gene regulation and asexuality are discussed.

Chromosomal inheritance of parental rDNAs distribution pattern detected by FISH in diploid F1 hybrid progeny of Cobitis (Teleostei, Cobitidae) species has non-Mendelian character

This study was conducted to describe the major and the minor rDNA chromosome distribution in the spined loach Cobitis taenia (2n = 48) and the Danubian loach Cobitis elongatoides (2n = 50), and their laboratory-produced diploid reciprocal F₁ hybrid progeny. It was tested by fluorescence in situ hybridisation (FISH) whether the number of 28s and 5s rDNA sites in the karyotypes of diploid hybrids corresponds to the expectations resulting from Mendelian ratio and if nucleolar organiser regions (NOR)were inherited from both parents or nucleolar dominance can be observed in the induced F₁ hybrid progeny. Ten (females) or twelve (males) 28s rDNA loci were located in nine uniarm chromosomes of C. taenia. Two of such loci terminally bounded on one acrocentric chromosome were unique and indicated as specific for this species. Large 5s rDNA clusters were located on two acrocentric chromosomes. In C. elongatoides of both sexes, six NOR sites in terminal regions on six meta-submetacentric chromosomes and two 5s rDNA sites on large submetacentrics were detected. The F₁ hybrid progeny (2n = 49) was characterised by the intermediate karyotype with the sites of ribosome synthesis on chromosomes inherited from both parents without showing nucleolar dominance. 5s rDNA sites were detected on large submetacentric and two acrocentric chromosomes. The observed number of both 28s and 5s rDNAs signals in F₁ diploid Cobitis hybrids was disproportionally inherited from the two parental species, showing inconsistency with the Mendelian ratios. The presented rDNA patterns indicate some marker chromosomes that allow the species of the parental male and female to be recognised in hybrid progeny. The 5s rDNA was found to be a particularly effective diagnostic marker of C. elongatoides to partially discern genomic composition of diploid Cobitis hybrids and presumably allopolyploids resulting from their backcrossing with one of the parental species. Thus, the current study provides insight into the extent of rDNA heredity in Cobitis chromosomes and their cytotaxonomic character.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Diversity of spined loaches from Asia Minor in a phylogenetic context (Teleostei: Cobitidae)

Accurate determination of species diversity in areas of high endemicity, particularly those lacking comprehensive systematic knowledge, represents a challenge for both taxonomists and conservationists. This need is particularly evident in areas greatly affected by anthropogenic disturbances such as the Eastern Mediterranean and its freshwater environments. To improve our knowledge of Eastern Mediterranean freshwater fishes, we phylogenetically studied Western Palearctic Cobitis species, focusing on those found in Turkey. Overall, our results provide a robust framework to assess the number of species of Cobitis. Phylogenetic reconstructions based on mitochondrial (cyt b) and nuclear (RAG1) sequences show seven major clades (Clades 1-7) grouping all Western Palearctic Cobitis species, except C. melanoleuca. In general, each major clade comprises Cobitis species that inhabit geographically close areas and have similar secondary sexual characters. Multiple divergent lineages were identified in our analyses, some of which were highly divergent such as the ones inhabiting Turkish freshwaters. Moreover, in some analyses, several of the identified lineages were incongruent with a priori defined species. Furthermore, our analyses identified eight potentially new candidate species, six that had been suggested in previous studies and two that are reported here for the first time. Our results reveal Turkey as the area with the greatest diversity of spined loaches in the Mediterranean.

Fine-scale genetic structure of the European bitterling at the intersection of three major European watersheds

BACKGROUND

Anthropogenic factors can have a major impact on the contemporary distribution of intraspecific genetic diversity. Many freshwater fishes have finely structured and locally adapted populations, but their natural genetic structure can be affected by river engineering schemes across river basins, fish transfers in aquaculture industry and conservation management. The European bitterling (Rhodeus amarus) is a small fish that is a brood parasite of freshwater mussels and is widespread across continental Europe. Its range recently expanded, following sharp declines in the 1970s and 1980s. We investigated its genetic variability and spatial structure at the centre of its distribution at the boundary of three watersheds, testing the role of natural and anthropogenic factors in its genetic structure.

RESULTS

Sequences of mitochondrial cytochrome B (CYTB) revealed that bitterling colonised central Europe from two Ponto-Caspian refugia, which partly defines its contemporary genetic structure. Twelve polymorphic microsatellite loci revealed pronounced interpopulation differentiation, with significant small-scale differentiation within the same river basins. At a large scale, populations from the Baltic Sea watershed (middle Oder and Vistula basins) were distinct from those from the Black Sea watershed (Danube basin), while populations from rivers of the North Sea watershed (Rhine, Elbe) originated from the admixture of both original sources. Notable exceptions demonstrated the potential role of human translocations across watersheds, with the upper River Oder (Baltic watershed) inhabited by fish from the Danube basin (Black Sea watershed) and a population in the southern part of the River Elbe (North Sea watershed) basin possessing a signal of admixture from the Danube basin.

CONCLUSIONS

Hydrography and physical barriers to dispersal are only partly reflected in the genetic structure of the European bitterling at the intersection of three major watersheds in central Europe. Drainage boundaries have been obscured by human-mediated translocations, likely related to common carp, Cyprinus carpio, cultivation and game-fish management. Despite these translocations, populations of bitterling are significantly structured by genetic drift, possibly reinforced by its low dispersal ability. Overall, the impact of anthropogenic factors on the genetic structure of the bitterling populations in central Europe is limited.

All-male hybrids of a tetrapod Pelophylax esculentus share its origin and genetics of maintenance

Background

Sexual parasites offer unique insights into the reproduction of unisexual and sexual populations. Because unisexuality is almost exclusively linked to the female sex, most studies addressed host-parasite dynamics in populations where sperm-dependent females dominate. Pelophylax water frogs from Central Europe include hybrids of both sexes, collectively named P. esculentus. They live syntopically with their parental species P. lessonae and/or P. ridibundus. Some hybrid lineages consist of all males providing a chance to understand the origin and perpetuation of a host-parasite (egg-dependent) system compared to sperm-dependent parthenogenesis.

Methods

We focused on P. ridibundus-P. esculentus populations where P. ridibundus of both sexes lives together with only diploid P. esculentus males. Based on 17 microsatellite markers and six allozyme loci, we analyzed (i) the variability of individual genomes, (ii) the reproductive mode(s) of all-male hybrids, and (iii) the genealogical relationships between the hybrid and parental genomes.

Results

Our microsatellite data revealed that P. esculentus males bear Mendelian-inherited ridibundus genomes while the lessonae genome represents a single clone. Our data indicate that this clone did not recently originate from adjacent P. lessonae populations, suggesting an older in situ or ex situ origin.

Conclusions

Our results confirm that also males can perpetuate over many generations as the unisexual lineage and successfully compete with P. ridibundus males for eggs provided by P. ridibundus females. Natural persistence of such sex-specific hybrid populations allows to studying the similarities and differences between male and female reproductive parasitism in many biological settings.

Electronic supplementary material

The online version of this article (10.1186/s13293-018-0172-z) contains supplementary material, which is available to authorized users.

Profound genetic divergence and asymmetric parental genome contributions as hallmarks of hybrid speciation in polyploid toads

The evolutionary causes and consequences of allopolyploidization, an exceptional pathway to instant hybrid speciation, are poorly investigated in animals. In particular, when and why hybrid polyploids versus diploids are produced, and constraints on sources of paternal and maternal ancestors, remain underexplored. Using the Palearctic green toad radiation (including bisexually reproducing species of three ploidy levels) as model, we generate a range-wide multi-locus phylogeny of 15 taxa and present four new insights: (i) at least five (up to seven) distinct allotriploid and allotetraploid taxa have evolved in the Pleistocene; (ii) all maternal and paternal ancestors of hybrid polyploids stem from two deeply diverged nuclear clades (6 Mya, 3.1-9.6 Mya), with distinctly greater divergence than the parental species of diploid hybrids found at secondary contact zones; (iii) allotriploid taxa possess two conspecific genomes and a deeply diverged allospecific one, suggesting that genomic imbalance and divergence are causal for their partly clonal reproductive mode; (iv) maternal versus paternal genome contributions exhibit asymmetry, with the maternal nuclear (and mitochondrial) genome of polyploids always coming from the same clade, and the paternal genome from the other. We compare our findings with similar patterns in diploid/polyploid vertebrates, and suggest deep ancestral divergence as a precondition for successful allopolyploidization.

Hybrid asexuality as a primary postzygotic barrier between nascent species: On the interconnection between asexuality, hybridization and speciation

Although sexual reproduction is ubiquitous throughout nature, the molecular machinery behind it has been repeatedly disrupted during evolution, leading to the emergence of asexual lineages in all eukaryotic phyla. Despite intensive research, little is known about what causes the switch from sexual reproduction to asexuality. Interspecific hybridization is one of the candidate explanations, but the reasons for the apparent association between hybridization and asexuality remain unclear. In this study, we combined cross-breeding experiments with population genetic and phylogenomic approaches to reveal the history of speciation and asexuality evolution in European spined loaches (Cobitis). Contemporary species readily hybridize in hybrid zones, but produce infertile males and fertile but clonally reproducing females that cannot mediate introgressions. However, our analysis of exome data indicates that intensive gene flow between species has occurred in the past. Crossings among species with various genetic distances showed that, while distantly related species produced asexual females and sterile males, closely related species produce sexually reproducing hybrids of both sexes. Our results suggest that hybridization leads to sexual hybrids at the initial stages of speciation, but as the species diverge further, the gradual accumulation of reproductive incompatibilities between species could distort their gametogenesis towards asexuality. Interestingly, comparative analysis of published data revealed that hybrid asexuality generally evolves at lower genetic divergences than hybrid sterility or inviability. Given that hybrid asexuality effectively restricts gene flow, it may establish a primary reproductive barrier earlier during diversification than other "classical" forms of postzygotic incompatibilities. Hybrid asexuality may thus indirectly contribute to the speciation process.

Reproductive dynamics shapes genomotype composition in an allopolyploid complex

Hybrid complexes are composed of organisms with multiple combinations of parental genomes (genomotypes) that interconnect through nets of crosses. Although several such complexes are well established without speciation or extinction, mechanisms shaping their dynamics remain poorly understood. In this study, we quantified the reproductive success of the allopolyploid Iberian fish Squalius alburnoides in experimental free-access and directional crosses involving the most common genomotypes. Specifically, we analysed the paternity of the offspring produced when females had free access to male genomotypes and quantified variations in egg allocation, fertilization rate, and offspring survival among crosses involving each male genomotype. The composition of the offspring produced from free-access crosses varied significantly from that expected from random mating, suggesting that offspring production and viability are not independent of parental male genomotype. Moreover, directional crosses producing the genomotype most commonly found in wild populations appeared to be the most successful, with females laying more eggs, and fertilization rate and offspring survival being the highest. These results suggest that reproductive dynamics plays a relevant role in structuring the genomotype composition of populations and opens a path to future research on the ecology and evolutionary biology of allopolyploids and their multiplicity of possible evolutionary pathways.

Asexual Reproduction Does Not Apparently Increase the Rate of Chromosomal Evolution: Karyotype Stability in Diploid and Triploid Clonal Hybrid Fish (Cobitis, Cypriniformes, Teleostei)

Interspecific hybridization, polyploidization and transitions from sexuality to asexuality considerably affect organismal genomes. Especially the last mentioned process has been assumed to play a significant role in the initiation of chromosomal rearrangements, causing increased rates of karyotype evolution. We used cytogenetic analysis and molecular dating of cladogenetic events to compare the rate of changes of chromosome morphology and karyotype in asexually and sexually reproducing counterparts in European spined loach fish (Cobitis). We studied metaphases of three sexually reproducing species and their diploid and polyploid hybrid clones of different age of origin. The material includes artificial F1 hybrid strains, representatives of lineage originated in Holocene epoch, and also individuals of an oldest known age to date (roughly 0.37 MYA). Thereafter we applied GISH technique as a marker to differentiate parental chromosomal sets in hybrids. Although the sexual species accumulated remarkable chromosomal rearrangements after their speciation, we observed no differences in chromosome numbers and/or morphology among karyotypes of asexual hybrids. These hybrids possess chromosome sets originating from respective parental species with no cytogenetically detectable recombinations, suggesting their integrity even in a long term. The switch to asexual reproduction thus did not provoke any significant acceleration of the rate of chromosomal evolution in Cobitis. Asexual animals described in other case studies reproduce ameiotically, while Cobitis hybrids described here produce eggs likely through modified meiosis. Therefore, our findings indicate that the effect of asexuality on the rate of chromosomal change may be context-dependent rather than universal and related to particular type of asexual reproduction.

Genetic diversity and distribution patterns of diploid and polyploid hybrid water frog populations (Pelophylax esculentus complex) across Europe

Polyploidization is a rare yet sometimes successful way for animals to rapidly create geno- and phenotypes that may colonize new habitats and quickly adapt to environmental changes. In this study, we use water frogs of the Pelophylax esculentus complex, comprising two species (Pelophylax lessonae, genotype LL; Pelophylax ridibundus, RR) and various diploid (LR) and triploid (LLR, LRR) hybrid forms, summarized as P. esculentus, as a model for studying recent hybridization and polyploidization in the context of speciation. Specifically, we compared the geographic distribution and genetic diversity of diploid and triploid hybrids across Europe to understand their origin, maintenance and potential role in hybrid speciation. We found that different hybrid and parental genotypes are not evenly distributed across Europe. Rather, their genetic diversity is structured by latitude and longitude and the presence/absence of parental species but not of triploids. Highest genetic diversity was observed in central and eastern Europe, the lowest in the northwestern parts of Europe. This gradient can be explained by the decrease in genetic diversity during postglacial expansion from southeastern glacial refuge areas. Genealogical relationships calculated on the basis of microsatellite data clearly indicate that hybrids are of multiple origin and include a huge variety of parental genomes. Water frogs in mixed-ploidy populations without any parental species (i.e. all-hybrid populations) can be viewed as evolutionary units that may be on their way towards hybrid speciation. Maintenance of such all-hybrid populations requires a continuous exchange of genomes between diploids and triploids, but scenarios for alternative evolutionary trajectories are discussed.

Mate Choice Drives Evolutionary Stability in a Hybrid Complex

Previous studies have shown that assortative mating acts as a driver of speciation by countering hybridization between two populations of the same species (pre-zygotic isolation) or through mate choice among the hybrids (hybrid speciation). In both speciation types, assortative mating promotes speciation over a transient hybridization stage. We studied mate choice in a hybrid vertebrate complex, the allopolyploid fish Squalius alburnoides. This complex is composed by several genomotypes connected by an intricate reproductive dynamics. We developed a model that predicts the hybrid complex can persist when females exhibit particular mate choice patterns. Our model is able to reproduce the diversity of population dynamic outcomes found in nature, namely the dominance of the triploids and the dominance of the tetraploids, depending on female mate choice patterns and frequency of the parental species. Experimental mate choice trials showed that females exhibit the preferences predicted by the model. Thus, despite the known role of assortative mating in driving speciation, our findings suggest that certain mate choice patterns can instead hinder speciation and support the persistence of hybrids over time without speciation or extinction.

Contrasting reproductive strategies of triploid hybrid males in vertebrate mating systems

The scarcity of parthenogenetic vertebrates is often attributed to their 'inferior' mode of clonal reproduction, which restricts them to self-reproduce their own genotype lineage and leaves little evolutionary potential with regard to speciation and evolution of sexual reproduction. Here, we show that for some taxa, such uniformity does not hold. Using hybridogenetic water frogs (Pelophylax esculentus) as a model system, we demonstrate that triploid hybrid males from two geographic regions exhibit very different reproductive modes. With an integrative data set combining field studies, crossing experiments, flow cytometry and microsatellite analyses, we found that triploid hybrids from Central Europe are rare, occur in male sex only and form diploid gametes of a single clonal lineage. In contrast, triploid hybrids from north-western Europe are widespread, occur in both sexes and produce recombined haploid gametes. These differences translate into contrasting reproductive roles between regions. In Central Europe, triploid hybrid males sexually parasitize diploid hybrids and just perpetuate their own genotype--which is the usual pattern in parthenogens. In north-western Europe, on the other hand, the triploid males are gamete donors for diploid hybrids, thereby stabilizing the mixed 2n-3n hybrid populations. By demonstrating these contrasting roles in male reproduction, we draw attention to a new significant evolutionary potential for animals with nonsexual reproduction, namely reproductive plasticity.

Ploidy-dependent survival of progeny arising from crosses between natural allotriploid Cobitis females and diploid C. taenia males (Pisces, Cobitidae)

Crosses between 21 triploid hybrid Cobitis females and 19 C. taenia (2n = 48) males led to viable progeny; whereas no embryonic development was observed in crosses with tetraploid males (4n = 98). The ploidy status of 491 progenies randomly selected with flow cytometry (316) or chromosome analysis (175) revealed an average of 55.2 % triploids and 44.8 % tetraploids, but the ratio of 3n versus 4n fish did change during development. In the first 2 days after hatching, approximately 65.1 % of tetraploid larvae were observed. Their number decreased significantly to 30.8 and 6.2 % on average during 2-5 and 10-15 months of life, respectively. The karyotype of tetraploid progeny (4n = 98) included 3n = 74 chromosomes of the parental female and n = 24 of C. taenia male. The number of tetraploid progeny indicated indirectly that about 66 % of eggs from 3n females were fertilized with C. taenia. The rest of the eggs developed clonally via gynogenesis or hemiclonally via hybridogenesis into triploids of the same karyotype structure as parental females. We have documented for the first time that (at least under experimental conditions) tetraploids are commonly formed, but are less viable than triploids, and a ratio similar to what is found under natural conditions is finally attained. The current explanation concerning the ploidy and karyotype structure of the progeny confirms that the eggs of 3n Cobitis females are not only capable of maintaining all chromosomes but are also capable of incorporating the sperm genome, thus creating the potential to produce tetraploids.

Distinguishing between incomplete lineage sorting and genomic introgressions: complete fixation of allospecific mitochondrial DNA in a sexually reproducing fish (Cobitis; Teleostei), despite clonal reproduction of hybrids

Distinguishing between hybrid introgression and incomplete lineage sorting causing incongruence among gene trees in that they exhibit topological differences requires application of statistical approaches that are based on biologically relevant models. Such study is especially challenging in hybrid systems, where usual vectors mediating interspecific gene transfers--hybrids with Mendelian heredity--are absent or unknown. Here we study a complex of hybridizing species, which are known to produce clonal hybrids, to discover how one of the species, Cobitis tanaitica, has achieved a pattern of mito-nuclear mosaic genome over the whole geographic range. We appplied three distinct methods, including the method using solely the information on gene tree topologies, and found that the contrasting mito-nuclear signal might not have resulted from the retention of ancestral polymorphism. Instead, we found two signs of hybridization events related to C. tanaitica; one concerning nuclear gene flow and the other suggested mitochondrial capture. Interestingly, clonal inheritance (gynogenesis) of contemporary hybrids prevents genomic introgressions and non-clonal hybrids are either absent or too rare to be detected among European Cobitis. Our analyses therefore suggest that introgressive hybridizations are rather old episodes, mediated by previously existing hybrids whose inheritance was not entirely clonal. Cobitis complex thus supports the view that the type of resulting hybrids depends on a level of genomic divergence between sexual species.

Comparative analysis of mitochondrial genomes in distinct nuclear ploidy loach Misgurnus anguillicaudatus and its implications for polyploidy evolution

Misgurnus anguillicaudatus has several natural ploidy types. To investigate whether nuclear polyploidy have an impact on mitochondrial DNA (mtDNA), the complete mitochondrial genomes (mitogenomes) of five distinct ploidy M. anguillicaudatus (natural diploid, triploid, tetraploid, pentaploid and hexaploid), which were collected in central China, were sequenced and analyzed. The five mitogenomes share the same gene arrangement and have similar gene size, base composition and codon usage pattern. The most variable regions of the mitogenome were the protein-coding genes, especially the ND4L (5.39% mutation rate). Most variations occurred in tetraploids. The phylogenetic tree showed that the tetraploid M. anguillicaudatus separated early from other ploidy loaches. Meanwhile, the mitogenomes from pentaploids, and hexaploids have the closest phylogenetic relations, but far from that of tetraploids, implying that pentaploids and hexaploids could not be formed from tetraploids, possibly from the diploids and triploids. The genus Misgurnus species were divided into two divergent inter-genus clades, and the five ploidy M. anguillicaudatus were monophyletic, which support the hypotheses about the mitochondrial introgression in loach species.

When a clonal genome finds its way back to a sexual species: evidence from ongoing but rare introgression in the hybridogenetic water frog complex

Besides several exceptions, asexual metazoans are usually viewed as ephemeral sinks for genomes, which become 'frozen' in clonal lineages after their emergence from ancestral sexual species. Here, we investigated whether and at what rate the asexuals are able to introgress their genomes back into the parental sexual population, thus more or less importantly affecting the gene pools of sexual species. We focused on hybridogenetic hybrids of western Palaearctic water frogs (Pelophylax esculentus), which originate through hybridization between P. ridibundus and P. lessonae, but transmit only clonal ridibundus genome into their gametes. Although usually mating with P. lessonae, P. esculentus may upon mating with P. ridibundus or another hybrid produce sexually reproducing P. ridibundus offspring with the introgressed ex-clonal genome. We compared the rate of nuclear amplified fragment length polymorphism (AFLP) and mitochondrial introgression in two types of populations, that is, those where P. ridibundus occurs in isolation and those where it lives with the hybridogens. Although significant differentiation (Φpt) between sexual and clonal ridibundus genomes suggested limited gene flow between sexuals and hybridogens, a non-negligible (~5%) proportion of P. ridibundus bore introgressed mtDNA and AFLP markers. Whereas transfer of mtDNA was exclusively unidirectional, introgression of nuclear markers was bidirectional. The proportion of introgressed P. ridibundus was highest in syntopic populations with P. esculentus, proving an ongoing and site-specific interspecific genetic transfer mediated by hybridogenetic hybrids. It turns out that asexual hybrids are not just a sink for genes of sexual species, but may significantly influence the genetic architecture of their sexual counterparts.

Let us not be unfair to asexuals: their ephemerality may be explained by neutral models without invoking any evolutionary constraints of asexuality

Phylogenetic studies typically demonstrate lower evolutionary ages of clones, relative to their sexual ancestors. This has often been attributed to heightened extinction risk of asexual organisms. We previously criticized such interpretations and demonstrated that the life span of clones is ultimately limited by neutral drift depending on the rate at which new clones are spawned into an asexual community of a finite size. Therefore, it is important to investigate whether the natural rates of such influxes are sufficiently high to account for the relative ephemerality of clones without assuming their increased extinction rate. I applied the neutral clonal turnover model to phylogenies of polyploid asexual ferns and simulated the coalescent trees over a wide range of demographic structures and sampling schemes. On parameterizing the model with biologically relevant estimates of population sizes and plant polyploidization rates, simulated clonal assemblages appeared younger than their sexual counterparts even in the absence of selection against clones. Therefore, differences observed between the ages of sexual and clonal lineages may be explained by the neutral clonal turnover. Researchers should consider the possibility that natural clones may get lost by neutral drift before their fate could eventually be affected by any long-term constraints of asexuality.