Unisexual salamanders (genus Ambystoma) present a new reproductive mode for eukaryotes

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomics resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomics resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, antipredator strategies, and resilience and adaptive responses. They also serve as critical models for understanding widespread genomic characteristics, including evolutionary genome expansions and contractions given they have the largest range in genome sizes of any animal taxon and multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The advent of long-read sequencing technologies, along with computational techniques that enhance scaffolding capabilities and streamline computational workload is now enabling the ability to overcome some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC) in early 2023. This burgeoning community already has more than 282 members from 41 countries (6 in Africa, 131 in the Americas, 27 in Asia, 29 in Australasia, and 89 in Europe). The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and outline how the AGC can enable amphibian genomics research to "leap" to the next level.

A cyclical switch of gametogenic pathways in hybrids depends on the ploidy level

The cellular and molecular mechanisms governing sexual reproduction are conserved across eukaryotes. Nevertheless, hybridization can disrupt these mechanisms, leading to asexual reproduction, often accompanied by polyploidy. In this study, we investigate how ploidy level and ratio of parental genomes in hybrids affect their reproductive mode. We analyze the gametogenesis of sexual species and their diploid and triploid hybrids from the freshwater fish family Cobitidae, using newly developed cytogenetic markers. We find that diploid hybrid females possess oogonia and oocytes with original (diploid) and duplicated (tetraploid) ploidy. Diploid oocytes cannot progress beyond pachytene due to aberrant pairing. However, tetraploid oocytes, which emerge after premeiotic genome endoreplication, exhibit normal pairing and result in diploid gametes. Triploid hybrid females possess diploid, triploid, and haploid oogonia and oocytes. Triploid and haploid oocytes cannot progress beyond pachytene checkpoint due to aberrant chromosome pairing, while diploid oocytes have normal pairing in meiosis, resulting in haploid gametes. Diploid oocytes emerge after premeiotic elimination of a single-copied genome. Triploid hybrid males are sterile due to aberrant pairing and the failure of chromosomal segregation during meiotic divisions. Thus, changes in ploidy and genome dosage may lead to cyclical alteration of gametogenic pathways in hybrids.

Evolutionary mechanisms and practical significance of reproductive success and clonal diversity in unisexual vertebrate polyploids

Unisexual reproduction is generally relevant to polyploidy, and unisexual vertebrates are often considered an evolutionary "dead end" due to the accumulation of deleterious mutations and absence of genetic diversity. However, some unisexual polyploids have developed strategies to avoid genomic decay, and thus provide ideal models to unveil unexplored evolutionary mechanisms, from the reproductive success to clonal diversity creation. This article reviews the evolutionary mechanisms for overcoming meiotic barrier and generating genetic diversity in unisexual vertebrates, and summarizes recent research advancements in the polyploid Carassius complex. Gynogenetic gibel carp (Carassius gibelio) is a unique amphitriploid that has undergone a recurrent autotriploidy and has overcome the bottleneck of triploid sterility via gynogenesis. Recently, an efficient strategy in which ploidy changes, including from amphitriploid to amphitetraploid, then from amphitetraploid to novel amphitriploid, drive unisexual-sexual-unisexual reproduction transition and clonal diversity has been revealed. Based on this new discovery, multigenomic reconstruction biotechnology has been used to breed a novel strain with superior growth and stronger disease resistance. Moreover, a unique reproduction mode that combines both abilities of ameiotic oogenesis and sperm-egg fusion, termed as ameio-fusiongensis, has been discovered, and it provides an efficient approach to synthesize sterile allopolyploids. In order to avoid ecological risks upon escape and protect the sustainable property rights of the aquaculture seed industry, a controllable fertility biotechnology approach for precise breeding is being developed by integrating sterile allopolyploid synthesis and gene-editing techniques. This review provides novel insights into the origin and evolution of unisexual vertebrates and into the attempts being made to exploit new breeding biotechnologies in aquaculture.

Nuclear gene introgressions in hybrid populations of water frog Pelophylax esculentus complex: geographical analysis of the phenomenon and its interpretation

Reproduction of water frog hybrids Pelophylax esculentus (Pelophylax ridibundus x Pelophylax lessonae) is associated with hemiclonal reproduction and backcrossing. The hemiclonal mode of reproduction occurs within P. esculentus allodiploids. In this case, the unrecombined genome of one parental species is transmitted to the offspring after premeiotic elimination of the chromosome set of the second parental species. Usually, the chromosome set of P. lessonae is eliminated, and the altered genome of P. ridibundus is passed on to the progeny. The hemiclonal inheritance within diploid Pelophylax esculentus hybrids may be accompanied by certain aberrations of premeiotic elimination. As a result, the formation of P. ridibundus specimens with introgressions of the P. lessonae genetic material, or the formation of recombinant hybrids occurs, depending on which of the parental species backcrossing takes place. The aim of our study is to describe the aberration of premeiotic elimination within the water frog P. esculentus complex detected by the nuclear gene Ldh-B inheritance, with an attempt to find out the causes of this phenomenon. It has been established that aberrations of premeiotic elimination are widespread, but only within populations of water frog from the river system of Upper Dnieper within Ukraine. The highest level of introgression takes place in the water frog populations within Kiev metropolis under conditions of expressed anthropogenization, while the maximum frequency of recombinants was detected within populations from the basin of Desna River, that has preserved native ecosystems. It was demonstrated that the frequency of premeiotic aberrations does not correlate with the intensity of interspecific water frog hybridization. Populations with introgressions are more common than populations with recombinants, however, within the latter, the frequency of recombination events is higher. The primary factor of gametogenesis aberrations, most likely, is the genetic characteristics of the local populations of parental species, since unambiguous explanations of this phenomenon based on the action of environmental stress (pollution of water systems) are not obvious.

Cytogenetics of the Hybridogenetic Frog Pelophylax grafi and Its Parental Species Pelophylax perezi

Hybrid taxa from the genus Pelophylax can propagate themselves in a modified way of sexual reproduction called hybridogenesis ensuring the formation of clonal gametes containing the genome of only one parental (host) species. Pelophylax grafi from South-Western Europe is a hybrid composed of P. ridibundus and P. perezi genomes and it lives with a host species P. perezi (P-G system). Yet it is unknown, whether non-Mendelian inheritance is fully maintained in such populations. In this study, we characterize P. perezi and P. grafi somatic karyotypes by using comparative genomic hybridization, genomic in situ hybridization, fluorescent in situ hybridization, and actinomycin D-DAPI. Here, we show the homeology of P. perezi and P. grafi somatic karyotypes to other Pelophylax taxa with 2n = 26 and equal contribution of ridibundus and perezi chromosomes in P. grafi which supports F1 hybrid genome constitution as well as a hemiclonal genome inheritance. We show that ridibundus chromosomes have larger regions of interstitial (TTAGGG)n repeats flanking the nucleolus organizing region on chromosome no. 10 and a high quantity of AT pairs in the centromeric regions. In P. perezi, we found species-specific sequences in metaphase chromosomes and marker structures in lampbrush chromosomes. Pericentromeric RrS1 repeat sequence was present in perezi and ridibundus chromosomes, but the blocks were stronger in ridibundus. Various cytogenetic techniques applied to the P-G system provide genome discrimination between ridibundus and perezi chromosomal sets. They could be used in studies of germ-line cells to explain patterns of clonal gametogenesis in P. grafi and broaden the knowledge about reproductive strategies in hybrid animals.

Diversity and composition of mixed-ploidy unisexual salamander assemblages reflect the key influence of host species

Understanding processes that govern and sustain biological diversity is a central goal of community ecology. Unisexual complexes, where reproduction depends on sperm from males of one or more bisexual host species, are rare and the processes driving their diversity and structure remain poorly understood. Unisexual Ambystoma salamanders produce distinct biotypes ('genomotypes') depending on which bisexual species they 'steal' sperm from. This reproductive mode should generate distinct assemblages depending on the locally available bisexual host species. Yet, how availability and relative abundance of multiple bisexual hosts influences composition and diversity of natural unisexual assemblages at local or regional scales remains unknown. We hypothesize that host identity most directly drives local assemblage composition, with host variation associated with increased beta and gamma diversity within unisexuals. We collected genetic samples from Ambystoma salamanders across Pelee Island, Ontario, Canada (2015-2022). Two host species were identified (A. texanum and A. laterale) with nine sites having a single host and one site having both. Unisexual assemblages were grouped into four clusters by similarity, with host identity being a key determinant. Gamma diversity increased as a result of distinct host-specific assemblages forming at different sites on the island (i.e., high beta diversity). Assemblage composition, but not diversity, was correlated with relative host abundance, which may reflect matching niche requirements between host and unisexual forms they produce. Our results demonstrate that diversity and structure of unisexual assemblages are clearly shaped by their host(s) and such systems may serve as models for studying how biotic interactions shape ecological communities.

An efficient approach to synthesize sterile allopolyploids through the combined reproduction mode of ameiotic oogenesis and sperm-egg fusion in the polyploid Carassius complex

The association between polyploidy and reproduction transition, which is an intriguing issue in evolutionary genetics, can also be exploited as an approach for genetic improvement in agriculture. Recently, we generated novel amphitriploids (NA3n) by integrating the genomes of the gynogenetic Carassius gibelio and sexual C. auratus, and found gynogenesis was recovered in most NA3n females (NA3n♀I). Here, we discovered a unique reproduction mode, termed ameio-fusiongenesis, which combines the abilities of both ameiotic oogenesis and sperm-egg fusion, in a few NA3n females (NA3n♀II). These females inherited ameiotic oogenesis to produce unreduced eggs from gynogenetic C. gibelio and sperm-egg fusion from sexual C. auratus. Subsequently, we utilized this unique reproduction mode to generate a group of synthetic alloheptaploids by crossing NA3n♀II with Megalobrama amblycephala. They contained all chromosomes of maternal NA3n♀II and a chromosomal set of paternal M. amblycephala. Intergenomic chromosome translocations between NA3n♀II and M. amblycephala were also observed in a few somatic cells. Primary oocytes of the alloheptaploid underwent severe apoptosis owing to incomplete double-strand break repair at prophase I. Although spermatocytes displayed similar chromosome behavior at prophase I, they underwent apoptosis due to chromosome separation failure at metaphase I. Therefore, the alloheptaploid females and males were all sterile. Finally, we established a sustainable clone for the large-scale production of NA3n♀II and developed an efficient approach to synthesize diverse allopolyploids containing genomes of different cyprinid species. These findings not only broaden our understanding of reproduction transition but also offer a practical strategy for polyploidy breeding and heterosis fixing.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

LAY SUMMARY

- What factors are common among bird species with cryptic diversity in Panama? What role do geography, ecology, phylogeographic history, and other factors play in generating bird diversity?- 19% of widely-sampled bird species form two or more distinct DNA barcode clades, suggesting widespread unrecognized diversity.- Traits associated with reduced dispersal ability, such as use of forest understory, high territoriality, low hand-wing index, and insectivory, were more common in taxa with cryptic diversity. Filogeografía comparada revela amplia diversidad críptica causada por la ecología en las aves de Panamá.

RESUMEN

Especies extendidas frecuentemente tiene diversidad genética no reconocida, y investigando los factores asociados con esta variación críptica puede ayudarnos a entender las fuerzas que impulsan la diversificación. Aquí, identificamos especies crípticas potenciales basadas en un conjunto de datos de códigos de barras de ADN mitocondrial de 2,333 individuos de aves de Panama en 429 especies, representando 391 (59%) de las 659 especies de aves terrestres residentes del país, además de algunas aves acuáticas muestreada de manera oportunista. Adicionalmente, complementamos estos datos con secuencias mitocondriales disponibles públicamente de otros loci, tal como ND2 o citocroma b, obtenidos de los genomas mitocondriales completos de 20 taxones. Utilizando los números de identificación de código de barras (en ingles: BINs), un sistema taxonómico numérico que proporcina una estimación imparcial de la diversidad potencial a nivel de especie, encontramos especies crípticas putativas en 19% de las especies de aves terrestres, lo que destaca la diversidad oculta en la avifauna bien descrita de Panamá. Aunque algunos de estos eventos de divergencia conciden con características geográficas que probablemente aislaron las poblaciones, la mayoría (74%) de la divergencia en las tierras bajas se encuentra entre las poblaciones orientales y occidentales. El tiempo de esta divergencia no coincidió entre los taxones, sugiriendo que eventos históricos tales como la formación del Istmo de Panamá y los ciclos climáticos del pleistoceno, no fueron los principales impulsores de la especiación. En cambio, observamos asociaciones fuertes entre las características ecológicas y la divergencia mitocondriale: las especies del bosque, sotobosque, con una dieta insectívora, y con territorialidad fuerte mostraton múltiple BINs probables. Adicionalmente, el índice mano-ala, que está asociado a la capacidad de dispersión, fue significativamente menor en las especies con BINs multiples, sugiriendo que la capacidad de dispersión tiene un rol importamente en la generación de la diversidad de las aves neotropicales. Estos resultos demonstran la necesidad de que estudios evolutivos de las comunidades de aves tropicales consideren los factores ecológicos en conjunto con las explicaciones geográficos. Palabras clave: biodiversidad tropical, biogeografía, códigos de barras, dispersión, especies crípticas.

Environmental Drivers of Amphibian Breeding Phenology across Multiple Sites

A mechanistic understanding of phenology, the seasonal timing of life history events, is important for understanding species’ interactions and the potential responses of ecological communities to a rapidly changing climate. We present analysis of a seven-year dataset on the breeding phenology of wood frogs (Rana sylvatica), tiger salamanders (Ambystoma tigrinum), blue-spotted salamanders (Ambystoma laterale), and associated unisexual Ambystoma salamanders from six wetlands in Southeast Michigan, USA. We assess whether the ordinal date of breeding migrations varies among species, sexes, and individual wetlands, and we describe the specific environmental conditions associated with breeding migrations for each species/sex. Breeding date was significantly affected by species/sex identity, year, wetland, and the interactions between species/sex and year as well as wetland and year. There was a great deal of variation among years, with breeding occurring nearly synchronously among groups in some years but widely spaced between groups in other years. Specific environmental triggers for movement varied for each species and sex and changed as the breeding season progressed. In general, salamanders responded to longer temperature lags (more warmer days in a row) than wood frogs, whereas wood frogs required longer precipitation lags (more rainy days in a row) than salamanders. Wood frogs were more likely to migrate around the time of a new moon, whereas in contrast, Ambystoma salamander migration was not associated with a moon phase. Ordinal day was an important factor in all models, suggesting that these amphibians require a latency period or similar mechanism to avoid breeding too early in the year, even when weather conditions appear favorable. Male wood frogs migrated earlier than female wood frogs, and male blue-spotted salamanders migrated earlier than female A. laterale and associated unisexual females. Larger unisexual salamanders migrated earlier than smaller individuals. Differences in species’ responses to environmental cues led to wood frogs and A. laterale breeding later than tiger salamanders in colder years but not in warmer years. This suggests that, as the climate warms, wood frog and A. laterale larvae may experience less predation from tiger salamander larvae due to reduced size differences when they breed simultaneously. Our study is one of few to describe the proximate drivers of amphibian breeding migrations across multiple species, wetlands, and years, and it can inform models predicting how climate change may shift ecological interactions among pond-breeding amphibian species.

Origin of scarlet gynogenetic triploid Carassius fish: Implications for conservation of the sexual-gynogenetic complex

Conservation of sperm-dependent asexual (gynogenetic) species is challenging due to their complicated ecological dynamics, which requires the stable coexistence with their sperm-providing sexual relatives, who often share similar niches. A symbolic but vulnerable gynogenetic animal is the scarlet Carassius fish, or Hibuna, which is mainly found in Lake Harutori on Hokkaido, Japan. Although Hibuna in Lake Harutori has been protected as a symbol of the Natural Monument of Japan, it has recently suffered population decline. To establish an effective conservation strategy for Hibuna, we investigated its origin, reproductive mode, and genetic diversity, with reference to the surrounding wild populations, using nuclear microsatellites and mitochondrial gene sequences. Our genetic analyses revealed that the main ploidy of Hibuna was triploid or tetraploid, and it reproduces gynogenetically. However, no co-existing sexual diploid Carassius was detected among our samples, suggesting that the sexual diploids and the gynogenetic population including Hibuna would be at risk of co-extirpation. In addition, Hibuna showed high genetic/clonal diversity and most Hibuna had nonindigenous mitochondrial haplotypes that are mostly identical to those reported from goldfish. These results indicate that Hibuna most probably originated from hybridization between indigenous gynogenetic triploids and goldfish introduced about 100 years ago, involving rare sexual reproduction. This spontaneous long-term field experiment exemplifies the recently documented diversification process of gynogenetic Carassius via complex interploidy gene flow. Although the priority to be placed on the conservation of Hibuna is controversial, the maintenance of gynogenetic Carassius, including Hibuna, requires strategic conservation of sexual populations.

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.

Conservation genomics of urban populations of Streamside Salamander (Ambystoma barbouri)

In Tennessee, populations of the state endangered Streamside Salamander (Ambystoma barbouri) are in decline as their distribution lies mostly within rapidly developing areas in the Nashville Basin. Information regarding the partitioning of genetic variation among populations of A. barbouri and the taxonomic status of these populations relative to northern populations and their congener, the Small-mouthed Salamander (A. texanum), have important implications for management and conservation of this species. Here we combined mitochondrial sequencing and genome-wide single nucleotide polymorphism (SNP) data generated using Genotyping-by-Sequencing (GBS) to investigate patterns of genetic variation within Tennessee populations of A. barbouri, to assess their relationship to populations in Kentucky, and to examine their phylogenetic relationship to the closely related A. texanum. Results from phylogenetic reconstructions reveal a complex history of Tennessee A. barbouri populations with regards to northern populations, unisexual A. barbouri, and A. texanum. Patterns of mitochondrial sequence variation suggest that A. barbouri may have originated within Tennessee and expanded north multiple times into Kentucky, Ohio, Indiana, and West Virginia. Phylogenetic reconstructions based on genome-wide SNP data contradict results based on mitochondrial DNA and correspond to geographic and taxonomic boundaries. Variation in allele frequencies at SNP genotypes, as identified by multivariate analyses and Bayesian assignment tests, identified three evolutionary significant units (ESUs) for A. barbouri within Tennessee. Collectively, these results emphasize the need for prioritizing conservation needs for Tennessee populations of A. barbouri to ensure the long-term persistence of this species.

Gigantic genomes of salamanders indicate that body temperature, not genome size, is the driver of global methylation and 5-methylcytosine deamination in vertebrates

Transposable elements (TEs) are sequences that replicate and move throughout genomes, and they can be silenced through methylation of cytosines at CpG dinucleotides. TE abundance contributes to genome size, but TE silencing variation across genomes of different sizes remains underexplored. Salamanders include most of the largest C-values - 9 to 120 Gb. We measured CpG methylation levels in salamanders with genomes ranging from 2N = ∼58 Gb to 4N = ∼116 Gb. We compared these levels to results from endo- and ectothermic vertebrates with more typical genomes. Salamander methylation levels are approximately 90%, higher than all endotherms. However, salamander methylation does not differ from other ectotherms, despite an approximately 100-fold difference in nuclear DNA content. Because methylation affects the nucleotide compositional landscape through 5-methylcytosine deamination to thymine, we quantified salamander CpG dinucleotide levels and compared them to other vertebrates. Salamanders and other ectotherms have comparable CpG levels, and ectotherm levels are higher than endotherms. These data show no shift in global methylation at the base of salamanders, despite a dramatic increase in TE load and genome size. This result is reconcilable with previous studies that considered endothermy and ectothermy, which may be more important drivers of methylation in vertebrates than genome size.

Premeiotic endoreplication is essential for obligate parthenogenesis in geckos

Obligate parthenogenesis evolved in reptiles convergently several times, mainly through interspecific hybridization. The obligate parthenogenetic complexes typically include both diploid and triploid lineages. Offspring of parthenogenetic hybrids are genetic copies of their mother; however, the cellular mechanism enabling the production of unreduced cells is largely unknown. Here, we show that oocytes go through meiosis in three widespread, or even strongly invasive, obligate parthenogenetic complexes of geckos, namely in diploid and triploid Lepidodactylus lugubris, and triploid Hemiphyllodactylus typus and Heteronotia binoei. In all four lineages, the majority of oocytes enter the pachytene at the original ploidy level, but their chromosomes cannot pair properly and instead form univalents, bivalents and multivalents. Unreduced eggs with clonally inherited genomes are formed from germ cells that had undergone premeiotic endoreplication, in which appropriate segregation is ensured by the formation of bivalents made from copies of identical chromosomes. We conclude that the induction of premeiotic endoreplication in reptiles was independently co-opted at least four times as an essential component of parthenogenetic reproduction and that this mechanism enables the emergence of fertile polyploid lineages within parthenogenetic complexes.

Delete and survive: strategies of programmed genetic material elimination in eukaryotes

Genome stability is a crucial feature of eukaryotic organisms because its alteration drastically affects the normal development and survival of cells and the organism as a whole. Nevertheless, some organisms can selectively eliminate part of their genomes from certain cell types during specific stages of ontogenesis. This review aims to describe the phenomenon of programmed DNA elimination, which includes chromatin diminution (together with programmed genome rearrangement or DNA rearrangements), B and sex chromosome elimination, paternal genome elimination, parasitically induced genome elimination, and genome elimination in animal and plant hybrids. During programmed DNA elimination, individual chromosomal fragments, whole chromosomes, and even entire parental genomes can be selectively removed. Programmed DNA elimination occurs independently in different organisms, ranging from ciliate protozoa to mammals. Depending on the sequences destined for exclusion, programmed DNA elimination may serve as a radical mechanism of dosage compensation and inactivation of unnecessary or dangerous genetic entities. In hybrids, genome elimination results from competition between parental genomes. Despite the different consequences of DNA elimination, all genetic material destined for elimination must be first recognised, epigenetically marked, separated, and then removed and degraded.

Interploidy gene flow involving the sexual-asexual cycle facilitates the diversification of gynogenetic triploid Carassius fish

Asexual vertebrates are rare and at risk of extinction due to their restricted adaptability through the loss of genetic recombination. We explore the mechanisms behind the generation and maintenance of genetic diversity in triploid asexual (gynogenetic) Carassius auratus fish, which is widespread in East Asian fresh waters and exhibits one of the most extensive distribution among asexual vertebrates despite its dependence on host sperm. Our analyses of genetic composition using dozens of genetic markers and genome-wide transcriptome sequencing uncover admixed genetic composition of Japanese asexual triploid Carassius consisting of both the diverged Japanese and Eurasian alleles, suggesting the involvement of Eurasian lineages in its origin. However, coexisting sexual diploid relatives and asexual triploids in Japan show regional genetic similarity in both mitochondrial and nuclear markers. These results are attributed to a unique unidirectional gene flow from diploids to sympatric triploids, with the involvement of occasional sexual reproduction. Additionally, the asexual triploid shows a weaker population structure than the sexual diploid, and multiple triploid lineages coexist in most Japanese rivers. The generated diversity via repeated interploidy gene flow as well as an increased establishment of immigrants is assumed to offset the cost of asexual reproduction and might contribute to the successful broad distribution of this asexual vertebrate.

Sex chromosomes in meiotic, hemiclonal, clonal and polyploid hybrid vertebrates: along the 'extended speciation continuum'

We review knowledge about the roles of sex chromosomes in vertebrate hybridization and speciation, exploring a gradient of divergences with increasing reproductive isolation (speciation continuum). Under early divergence, well-differentiated sex chromosomes in meiotic hybrids may cause Haldane-effects and introgress less easily than autosomes. Undifferentiated sex chromosomes are more susceptible to introgression and form multiple (or new) sex chromosome systems with hardly predictable dominance hierarchies. Under increased divergence, most vertebrates reach complete intrinsic reproductive isolation. Slightly earlier, some hybrids (linked in 'the extended speciation continuum') exhibit aberrant gametogenesis, leading towards female clonality. This facilitates the evolution of various allodiploid and allopolyploid clonal ('asexual') hybrid vertebrates, where 'asexuality' might be a form of intrinsic reproductive isolation. A comprehensive list of 'asexual' hybrid vertebrates shows that they all evolved from parents with divergences that were greater than at the intraspecific level (K2P-distances of greater than 5-22% based on mtDNA). These 'asexual' taxa inherited genetic sex determination by mostly undifferentiated sex chromosomes. Among the few known sex-determining systems in hybrid 'asexuals', female heterogamety (ZW) occurred about twice as often as male heterogamety (XY). We hypothesize that pre-/meiotic aberrations in all-female ZW-hybrids present Haldane-effects promoting their evolution. Understanding the preconditions to produce various clonal or meiotic allopolyploids appears crucial for insights into the evolution of sex, 'asexuality' and polyploidy. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Testis Development and Differentiation in Amphibians

Sex is determined genetically in amphibians; however, little is known about the sex chromosomes, testis-determining genes, and the genes involved in testis differentiation in this class. Certain inherent characteristics of the species of this group, like the homomorphic sex chromosomes, the high diversity of the sex-determining mechanisms, or the existence of polyploids, may hinder the design of experiments when studying how the gonads can differentiate. Even so, other features, like their external development or the possibility of inducing sex reversal by external treatments, can be helpful. This review summarizes the current knowledge on amphibian sex determination, gonadal development, and testis differentiation. The analysis of this information, compared with the information available for other vertebrate groups, allows us to identify the evolutionarily conserved and divergent pathways involved in testis differentiation. Overall, the data confirm the previous observations in other vertebrates-the morphology of the adult testis is similar across different groups; however, the male-determining signal and the genetic networks involved in testis differentiation are not evolutionarily conserved.

Skipping sex: A nonrecombinant genomic assemblage of complementary reproductive modules

The unusual occurrence and developmental diversity of asexual eukaryotes remain a puzzle. De novo formation of a functioning asexual genome requires a unique assembly of sets of genes or gene states to disrupt cellular mechanisms of meiosis and gametogenesis, and to affect discrete components of sexuality and produce clonal or hemiclonal offspring. We highlight two usually overlooked but essential conditions to understand the molecular nature of clonal organisms, that is, a nonrecombinant genomic assemblage retaining modifiers of the sexual program, and a complementation between altered reproductive components. These subtle conditions are the basis for physiologically viable and genetically balanced transitions between generations. Genomic and developmental evidence from asexual animals and plants indicates the lack of complementation of molecular changes in the sexual reproductive program is likely the main cause of asexuals' rarity, and can provide an explanatory frame for the developmental diversity and lability of developmental patterns in some asexuals as well as for the discordant time to extinction estimations.

Genotypic similarities among the parthenogenetic Darevskia rock lizards with different hybrid origins

BACKGROUND

The majority of parthenogenetic vertebrates derive from hybridization between sexually reproducing species, but the exact number of hybridization events ancestral to currently extant clonal lineages is difficult to determine. Usually, we do not know whether the parental species are able to contribute their genes to the parthenogenetic vertebrate lineages after the initial hybridization. In this paper, we address the hypothesis, whether some genotypes of seven phenotypically distinct parthenogenetic rock lizards (genus Darevskia) could have resulted from back-crosses of parthenogens with their presumed parental species. We also tried to identify, as precise as possible, the ancestral populations of all seven parthenogens.

RESULTS

We analysed partial mtDNA sequences and microsatellite genotypes of all seven parthenogens and their presumed ansectral species, sampled across the entire geographic range of parthenogenesis in this group. Our results confirm the previous designation of the parental species, but further specify the maternal populations that are likely ancestral to different parthenogenetic lineages. Contrary to the expectation of independent hybrid origins of the unisexual taxa, we found that genotypes at multiple loci were shared frequently between different parthenogenetic species. The highest proportions of shared genotypes were detected between (i) D. sapphirina and D. bendimahiensis and (ii) D. dahli and D. armeniaca, and less often between other parthenogens. In case (ii), genotypes at the remaining loci were notably distinct.

CONCLUSIONS

We suggest that both observations (i-ii) can be explained by two parthenogenetic forms tracing their origin to a single initial hybridization event. In case (ii), however, occasional gene exchange between the unisexual and the parental bisexual species could have taken place after the onset of parthenogenetic reproduction. Indeed, backcrossed polyploid hybrids are relatively frequent in Darevskia, although no direct evidence of recent gene flow has been previously documented. Our results further suggest that parthenogens are losing heterozygosity as a result of allelic conversion, hence their fitness is expected to decline over time as genetic diversity declines. Backcrosses with the parental species could be a rescue mechanism which might prevent this decline, and therefore increase the persistance of unisexual forms.

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.

Unexpected Oogenic Pathways for the Triploid Fish Chrosomus eos-neogaeus

Triploid vertebrates from unisexual complexes often perpetuate themselves asexually. In the fish Chrosomus eos×eos-neogaeus, triploids are continuously produced by diploid hybrids. However, they are not expected to perpetuate since C. eos are their only known progeny. This study aims to investigate the oogenesis of these triploid hybrids through experimental crosses. A total of 337 larvae from 12 female triploids and 3 2n/3n mosaics fertilized with C. eos sperm were genetically characterized. The detection of C. eos as progeny of triploid hybrids confirmed the occurrence of a pathway similar to meiotic hybridogenesis but only for half of the tripoids. The presence of tetraploid offspring for all these females revealed the formation of unreduced triploid eggs as a probable failure of meiotic hybridogenesis. The remaining female triploids and all mosaics produced diploid and triploid hybrids. Triploids excluded the haplome from paternal leakage and produced eggs with the diploid hybrid genome through an ameiotic hybridogenesis. Both types of hybridogenesis occurred in a mutually exclusive manner. This leads us to consider 2 hypothetical scenarios: First, any female triploids can perform either type of hybridogenesis, allowing the long-term persistence of triploid hybrids by a fraction of the population. Alternatively, ameiotic hybridogenesis occurs in triploids of the first generation (from diploid mothers), while meiotic hybridogenesis occurs in triploids of the second generation (from triploid mothers); triploid hybrids then are not perpetuating lineages. The population dynamics of the C. eos-neogaeus complex appears a step more complicated than previously expected.

Micronuclei in germ cells of hybrid frogs from Pelophylax esculentus complex contain gradually eliminated chromosomes

In most organisms, cells typically maintain genome integrity, as radical genome reorganization leads to dramatic consequences. However, certain organisms, ranging from unicellular ciliates to vertebrates, are able to selectively eliminate specific parts of their genome during certain stages of development. Moreover, partial or complete elimination of one of the parental genomes occurs in interspecies hybrids reproducing asexually. Although several examples of this phenomenon are known, the molecular and cellular processes involved in selective elimination of genetic material remain largely undescribed for the majority of such organisms. Here, we elucidate the process of selective genome elimination in water frog hybrids from the Pelophylax esculentus complex reproducing through hybridogenesis. Specifically, in the gonads of diploid and triploid hybrids, but not those of the parental species, we revealed micronuclei in the cytoplasm of germ cells. In each micronucleus, only one centromere was detected with antibodies against kinetochore proteins, suggesting that each micronucleus comprises a single chromosome. Using 3D-FISH with species-specific centromeric probe, we determined the role of micronuclei in selective genome elimination. We found that in triploid LLR hybrids, micronuclei preferentially contain P. ridibundus chromosomes, while in diploid hybrids, micronuclei preferentially contain P. lessonae chromosomes. The number of centromere signals in the nuclei suggested that germ cells were aneuploid until they eliminate the whole chromosomal set of one of the parental species. Furthermore, in diploid hybrids, misaligned P. lessonae chromosomes were observed during the metaphase stage of germ cells division, suggesting their possible elimination due to the inability to attach to the spindle and segregate properly. Additionally, we described gonocytes with an increased number of P. ridibundus centromeres, indicating duplication of the genetic material. We conclude that selective genome elimination from germ cells of diploid and triploid hybrids occurs via the gradual elimination of individual chromosomes of one of the parental genomes, which are enclosed within micronuclei.

Triad hybridization via a conduit species

Introgressive hybridization can affect the evolution of populations in several important ways. It may retard or reverse divergence of species, enable the development of novel traits, enhance the potential for future evolution by elevating levels of standing variation, create new species, and alleviate inbreeding depression in small populations. Most of what is known of contemporary hybridization in nature comes from the study of pairs of species, either coexisting in the same habitat or distributed parapatrically and separated by a hybrid zone. More rarely, three species form an interbreeding complex (triad), reported in vertebrates, insects, and plants. Often, one species acts as a genetic link or conduit for the passage of genes (alleles) between two others that rarely, if ever, hybridize. Demographic and genetic consequences are unknown. Here we report results of a long-term study of interbreeding Darwin's finches on Daphne Major island, Galápagos. Geospiza fortis acted as a conduit for the passage of genes between two others that have never been observed to interbreed on Daphne: Geospiza fuliginosa, a rare immigrant, and Geospiza scandens, a resident. Microsatellite gene flow from G. fortis into G. scandens increased in frequency during 30 y of favorable ecological conditions, resulting in genetic and morphological convergence. G. fortis, G. scandens, and the derived dihybrids and trihybrids experienced approximately equal fitness. Especially relevant to young adaptive radiations, where species differ principally in ecology and behavior, these findings illustrate how new combinations of genes created by hybridization among three species can enhance the potential for evolutionary change.

Stable Genome Incorporation of Sperm-derived DNA Fragments in Gynogenetic Clone of Gibel Carp

How unisexual animals eliminate deleterious mutations to avoid dead ends is one of the most interesting puzzles in evolutionary genetics. Incorporation of microchromosomes derived from exogenous sperm had been observed in gynogenetic animals, but little is known about their detailed process and hereditary fate. Here, we show a stable genome incorporation case in an artificial clone F of gynogenetic gibel carp (Carassius gibelio). A total of 12 exogenous DNA fragments were screened through a read depth-dependent comparison strategy and confirmed to be specific to the clone F and the paternal blunt snout bream (Megalobrama amblycephala Yin) by SCAR (sequence characterized amplified regions) marker detection. Moreover, these sperm-derived DNA fragments were not detected in some samples in early gynogenetic generations, but they were found to exist in all examined individuals through artificial gynogenetic selections of 13 generations, implying that they might have stably incorporated into the genome of clone F. Furthermore, chromosome localization and sequence characterization indicate that the largest fragment CgA22_34 is derived from blunt snout bream non-LTR retrotransposon and durably incorporated into only one of three homologous chromosomes of gibel carp clone F. Our results suggest that the incorporated sperm-derived DNA fragments by allogynogenesis should increase genetic diversity and introduce new traits into unisexual animals which will benefit genetic breeding of gibel carp. During the process, transposable elements (TEs) may play significant roles in shaping the genome structures. Simultaneously, the incorporated DNA fragments are able to be used as genetic markers to perform selective breeding programs in aquaculture practices of gibel carp.

Estimating critical habitat based on year-round movements of the endangered Jefferson Salamander (Ambystoma jeffersonianum) and their unisexual dependents

Habitat protection is a key component of endangered species conservation, but critical habitat designations are often based on limited data or habitat use during only a portion of a species’ life cycle. Protected habitat around breeding pools for the endangered Jefferson Salamander (Ambystoma jeffersonianum (Green, 1827)) and their unisexual dependents (Ambystoma laterale – (2) jeffersonianum) is based upon limited movement data from primarily spring and summer seasons. Furthermore, despite their treatment as distinct species under Canada’s Species at Risk Act, differences in habitat use have not been directly compared in areas where they co-occur. We used radiotelemetry to track A. jeffersonianum (JJ) and A. laterale – jeffersonianum (LJJ and LJJJ) during fall migrations to overwintering sites. We used these data and summarized available movement data from past studies that tracked movements in other periods of the annual cycle to estimate year-round critical habitat. Ambystoma jeffersonianum travelled significantly longer distances to overwintering locations than unisexuals. Individuals were more likely to overwinter next to a similar genomotype individual than not. Critical habitat encompassing the entire annual cycle of A. jeffersonianum extends up to 400–450 m from breeding ponds indicating existing regulatory habitat protections in Canada do not currently protect sufficient habitat.

Survival of Polyploid hybrid salamander embryos

Background

Animals with polyploid, hybrid nuclei offer a challenge for models of gene expression and regulation during embryogenesis. To understand how such organisms proceed through development, we examined the timing and prevalence of mortality among embryos of unisexual salamanders in the genus Ambystoma.

Results

Our regional field surveys suggested that heightened rates of embryo mortality among unisexual salamanders begin in the earliest stages of embryogenesis. Although we expected elevated mortality after zygotic genome activation in the blastula stage, this is not what we found among embryos which we reared in the laboratory. Once embryos entered the first cleavage stage, we found no difference in mortality rates between unisexual salamanders and their bisexual hosts. Our results are consistent with previous studies showing high rates of unisexual mortality, but counter to reports that heightened embryo mortality continues throughout embryo development.

Conclusions

Possible causes of embryonic mortality in early embryogenesis suggested by our results include abnormal maternal loading of RNA during meiosis and barriers to insemination. The surprising survival rates of embryos post-cleavage invites further study of how genes are regulated during development in such polyploid hybrid organisms.

A family study to examine clonal diversity in unisexual salamanders (genus Ambystoma)

Unisexual Ambystoma are the oldest known unisexual vertebrates and comprise a lineage of eastern North American all female salamanders that reproduce by stealing sperm from as many as five normally bisexual congeneric species. The sperm may be used to only stimulate egg development by gynogenesis but can be incorporated in the zygote to elevate the ploidy level or to replace one of the female’s haploid genomes. This flexible and unique reproductive system, termed kleptogenesis, is investigated using a microsatellite examination of 988 offspring from 14 unisexual mothers. All mothers produced clonal and ploidy-elevated offspring. Genome replacement and multiple paternity are confirmed for the first time in unisexual Ambystoma. Microsatellite mutations were found in all five microsatellite loci and the estimated microsatellite mutation rate varied by locus and by genome. Clonal variation is attributed to the inclusion of sperm donors’ haploid genomes for ploidy elevation, genome replacement, mutations, and natural selection.

Naturally clonal vertebrates are an untapped resource in ecology and evolution research

Science requires replication. The development of many cloned or isogenic model organisms is a testament to this. But researchers are reluctant to use these traditional animal model systems for certain questions in evolution or ecology research, because of concerns over relevance or inbreeding. It has largely been overlooked that there are a substantial number of vertebrate species that reproduce clonally in nature. Here we highlight how use of these naturally evolved, phenotypically complex animals can push the boundaries of traditional experimental design and contribute to answering fundamental questions in the fields of ecology and evolution.

Sperm-dependent asexual hybrids determine competition among sexual species

Interspecific competition is a fundamental process affecting community structure and evolution of interacting species. Besides direct competition, this process is also mediated by shared enemies, which can change the outcome of competition dramatically. However, previous studies investigating interactions between competing species and their parasites (parasite-mediated competition) completely overlooked the effect of ‘sperm’ parasites (i.e. sperm-dependent parthenogens or pseudogams) on competition. These organisms originate by interspecific hybridization, produce clonal gametes, but exploit parental species for their own reproduction, being therefore analogous to classical parasites. Here we use the reaction-diffusion model and show that pseudogams alter the outcome of interspecific competition significantly. They may either slow down competitive exclusion of the inferior competitor or even turn the outcome of competition between the species. Asexual organisms may thus have unexpectedly strong impact on community structure, and have more significant evolutionary potential than was previously thought.

The Biogeography of Deep Time Phylogenetic Reticulation

Most phylogenies are typically represented as purely bifurcating. However, as genomic data have become more common in phylogenetic studies, it is not unusual to find reticulation among terminal lineages or among internal nodes (deep time reticulation; DTR). In these situations, gene flow must have happened in the same or adjacent geographic areas for these DTRs to have occurred and therefore biogeographic reconstruction should provide similar area estimates for parental nodes, provided extinction or dispersal has not eroded these patterns. We examine the phylogeny of the widely distributed New World kingsnakes (Lampropeltis), determine if DTR is present in this group, and estimate the ancestral area for reticulation. Importantly, we develop a new method that uses coalescent simulations in a machine learning framework to show conclusively that this phylogeny is best represented as reticulating at deeper time. Using joint probabilities of ancestral area reconstructions on the bifurcating parental lineages from the reticulating node, we show that this reticulation likely occurred in northwestern Mexico/southwestern US, and subsequently, led to the diversification of the Mexican kingsnakes. This region has been previously identified as an area important for understanding speciation and secondary contact with gene flow in snakes and other squamates. This research shows that phylogenetic reticulation is common, even in well-studied groups, and that the geographic scope of ancient hybridization is recoverable.

Effect of a locally adapted genome on environmentally induced epigenetic variation

Both genetic variation and environmentally induced epigenetic changes allow organisms to persist through the heterogeneity of their habitats. Selection on genetic variation can promote local adaptation of populations. However, in absence of genetic variation, clonal organisms mostly rely on epigenetics to respond to environmental heterogeneity. We used the potential of unisexual organisms in incorporating their host genome, to empirically assess whether the presence of a locally adapted genome affects environmentally induced epigenetic changes in clonal organisms. We addressed this problematic by using unisexual lineages of the kleptogen vertebrate Ambystoma laterale-jeffersonianum complex that can optionally incorporate genetic material from locally adapted sexual hosts through genomic exchanges. More specifically, we compared environmentally induced epigenetic changes between lineages strictly reproducing clonally vs. those incorporating a locally adapted genome. The results revealed that both lineage and sample site components, as well as their interaction, affected epigenetic variation. When lineages were analysed separately, differences among sample sites were only detected in lineages impervious to genomic exchanges. Sample sites had no significant effect on the epigenetic variation of lineages that performed genomic exchanges. These results suggest that environmentally induced epigenetic variation among sites depends more on the lack of locally adapted alleles than on the level of genetic variation.

Genome-specific histories of divergence and introgression between an allopolyploid unisexual salamander lineage and two ancestral sexual species

Quantifying introgression between sexual species and polyploid lineages traditionally thought to be asexual is an important step in understanding what drives the longevity of putatively asexual groups. Here, we capitalize on three recent innovations-ultraconserved element (UCE) sequencing, bioinformatic techniques for identifying genome-specific variation in polyploids, and model-based methods for evaluating historical gene flow-to measure the extent and tempo of introgression over the evolutionary history of an allopolyploid lineage of all-female salamanders and two ancestral sexual species. Our analyses support a scenario in which the genomes sampled in unisexual salamanders last shared a common ancestor with genomes in their parental species ∼3.4 million years ago, followed by a period of divergence between homologous genomes. Recently, secondary introgression has occurred at different times with each sexual species during the last 500,000 years. Sustained introgression of sexual genomes into the unisexual lineage is the defining characteristic of their reproductive mode, but this study provides the first evidence that unisexual genomes have undergone long periods of divergence without introgression. Unlike other sperm-dependent taxa in which introgression is rare, the alternating periods of divergence and introgression between unisexual salamanders and their sexual relatives could explain why these salamanders are among the oldest described unisexual animals.

Influence of genome and bio-ecology on the prevalence of genome exchange in unisexuals of the Ambystoma complex

Background

Unisexuals of the blue-spotted salamander complex are thought to reproduce by kleptogenesis. Genome exchanges associated with this sperm-dependent mode of reproduction are expected to result in a higher genetic variation and multiple ploidy levels compared to clonality. However, the existence of some populations exclusively formed of genetically identical individuals suggests that factors could prevent genome exchanges. This study aimed at assessing the prevalence of genome exchange among unisexuals of the Ambystoma laterale-jeffersonianum complex from 10 sites in the northern part of their distribution.

Results

A total of 235 individuals, including 207 unisexuals, were genotyped using microsatellite loci and AFLP. Unisexual individuals could be sorted in five genetically distinct groups, likely derived from the same paternal A. jeffersonianum haplome. One of these groups exclusively reproduced clonally, even when found in sympatry with lineages presenting signature of genome exchange. Genome exchange was site-dependent for another group. Genome exchange was detected at all sites for the three remaining groups.

Conclusion

Prevalence of genome exchange appears to be associated with ecological conditions such as availability of effective sperm donors. Intrinsic genomic factors may also affect this process, since different lineages in sympatry present highly variable rate of genome exchange. The coexistence of clonal and genetically diversified lineages opens the door to further research on alternatives to genetic variation.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1200-7) contains supplementary material, which is available to authorized users.

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.

General environmental heterogeneity as the explanation of sexuality? Comparative study shows that ancient asexual taxa are associated with both biotically and abiotically homogeneous environments

Ecological theories of sexual reproduction assume that sexuality is advantageous in certain conditions, for example, in biotically or abiotically more heterogeneous environments. Such theories thus could be tested by comparative studies. However, the published results of these studies are rather unconvincing. Here, we present the results of a new comparative study based exclusively on the ancient asexual clades. The association with biotically or abiotically homogeneous environments in these asexual clades was compared with the same association in their sister, or closely related, sexual clades. Using the conservative definition of ancient asexuals (i.e., age >1 million years), we found eight pairs of taxa of sexual and asexual species, six differing in the heterogeneity of their inhabited environment on the basis of available data. The difference between the environmental type associated with the sexual and asexual species was then compared in an exact binomial test. The results showed that the majority of ancient asexual clades tend to be associated with biotically, abiotically, or both biotically and abiotically more homogeneous environments than their sexual controls. In the exploratory part of the study, we found that the ancient asexuals often have durable resting stages, enabling life in subjectively homogeneous environments, live in the absence of intense biotic interactions, and are very often sedentary, inhabiting benthos, and soil. The consequences of these findings for the ecological theories of sexual reproduction are discussed.

Mutual maintenance of di- and triploid Pelophylax esculentus hybrids in R-E systems: results from artificial crossings experiments

Background

Interspecies animal hybrids can employ clonal or hemiclonal reproduction modes where one or all parental genomes are transmitted to the progeny without recombination. Nevertheless, some interspecies hybrids retain strong connection with the parental species needed for successful reproduction. Appearance of polyploid hybrid animals may play an important role in the substitution of parental species and in the speciation process.

Results

To establish the mechanisms that enable parental species, diploid and polyploid hybrids coexist we have performed artificial crossing experiments of water frogs of Pelophylax esculentus complex. We identified tadpole karyotypes and oocyte genome composition in all females involved in the crossings. The majority of diploid and triploid hybrid frogs produced oocytes with 13 bivalents leading to haploid gametes with the same genome as parental species hybrids usually coexist with. After fertilization of such gametes only diploid animals appeared. Oocytes with 26 bivalents produced by some diploid hybrid frogs lead to diploid gametes, which give rise to triploid hybrids after fertilization. In gonads of all diploid and triploid hybrid tadpoles we found DAPI-positive micronuclei (nucleus-like bodies) involved in selective genome elimination. Hybrid male and female individuals produced tadpoles with variable karyotype and ploidy even in one crossing owing to gametes with various genome composition.

Conclusions

We propose a model of diploid and triploid hybrid frog reproduction in R-E population systems. Triploid Pelophylax esculentus hybrids can transmit genome of parental species they coexist with by producing haploid gametes with the same genome composition. Triploid hybrids cannot produce triploid individuals after crossings with each other and depend on diploid hybrid females producing diploid eggs. In contrast to other population systems, the majority of diploid and triploid hybrid females unexpectedly produced gametes with the same genome as parental species hybrids coexist with.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1063-3) contains supplementary material, which is available to authorized users.

Genome Expression Balance in a Triploid Trihybrid Vertebrate

Polyploidy is increasingly recognized as a driver of biological diversity. How and why polyploidization affects gene expression is critical to understanding the link between ploidy elevation and diversification. In polyploid plants, multiple studies have demonstrated that ploidy elevation can confer major but variable consequences for gene expression, ranging from gene-by-gene alterations to entirely silenced genomes. By contrast, animal polyploids remain largely uncharacterized. Accordingly, how animals respond to and manage polyploidy events is not understood. Here, we address this important knowledge gap by analyzing transcriptomes from a triploid hybrid animal, a unisexual Ambystoma salamander, and three sexual Ambystoma species that represent all three parental genomes in the unisexual. We used a novel bioinformatics pipeline that includes competitively mapping triploid sequences to a reference set of orthologous genes in the sexual species to evaluate subgenome expression. Our comparisons of gene expression levels across the three parental genomes revealed that the unisexual triploid displays a pattern of genome balance, where 72% of the genes analyzed were expressed equally among the subgenomes. This result is strikingly different from the genome imbalance typically observed in hybrid polyploid plants. Our analyses represent the first to address gene expression in a triploid hybrid animal and introduce a novel bioinformatic framework for analyzing transcriptomic data.