Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs.

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

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.

A pheromone by any other gene would smell as sweet.

Reproductive isolation is the result of either the inability to produce viable and fertile offspring or the avoidance of mating altogether. While these mechanisms can evolve either over time via genetic drift or natural selection, the genetic result is usually a complex set of traits that are often linked. Explaining how reproductive isolation proceeds from the initiation of divergence to the complete prevention of mating is often a difficult task, as the underlying genes for traits associated with reproductive isolation can change via molecular evolution and subsequent protein coding alterations or through alterations of gene expression regulation. In this issue of Molecular Ecology, Treer, Maex, VanBocxlaer, Proost, and Bossuyt () use transcriptomic, proteomic and phylogenetic analyses to show that species-specific sex pheromones are the result of gradual sequence divergence on the same set of proteins in two closely related newt species (Ichthyosaura alpestris and Lissotriton helveticus). This study shows that salamander pheromone systems provide an enticing opportunity to connect the evolution of reproductive isolation to the changes in genes that underlie a key phenotype.

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.

No safety in the trees: Local and species-level adaptation of an arboreal squirrel to the venom of sympatric rattlesnakes.

Within some species, squirrels respond to variable selection from venomous snake predators by showing population-level variation in resistance, while between species, some rattlesnakes possess venom that is more effective at overcoming venom resistance in different species of squirrels. A functional evaluation of resistance variation to venom within and between species of squirrels and snakes can link resistance variation to its evolutionary causes across these different evolutionary scales. To do this, we compared the effectiveness of squirrel sera in inhibiting rattlesnake (Crotalus spp.) venom metalloproteinase activity between populations and between species to test for a response to local variation in selection from a single rattlesnake predator and for specialization of two resistant squirrel species to each of their distinct sympatric snake predators. We found that Timber Rattlesnake (Crotalus horridus) venom inhibition by Eastern gray squirrels (Sciurus carolinensis) is higher at a site where the rattlesnakes are present, which suggests selection may maintain venom resistance in populations separated by short distances. Next, we performed a reciprocal cross of venoms and sera from two rattlesnake and two squirrel species. This showed that squirrel resistance is lower when tested against venom from allopatric compared to sympatric rattlesnake species, demonstrating that squirrel inhibitors are specialized to sympatric venom and suggesting a tradeoff in terms of specialization to the venom of a specific species of rattlesnake predator. This pattern can be explained if inhibitors must recognize venom proteins and resistance evolution tracks venom evolution.

Cryptic sex? Estimates of genome exchange in unisexual mole salamanders (Ambystoma sp.).

Cryptic sex has been argued to explain the exceptional longevity of certain parthenogenetic vertebrate lineages, yet direct measurements of genetic exchange between sexual and apparently parthenogenetic forms are rare. Female unisexual mole salamanders (Ambystoma sp.) are the oldest known unisexual vertebrate lineage (~5 million years), and one hypothesis for their persistence is that allopolyploid female unisexuals periodically exchange haploid genomes 'genome exchange' during gynogenetic reproduction with males from sympatric sexual species. We test this hypothesis by using genome-specific microsatellite DNA markers to estimate the rates of genome exchange between sexual males and unisexual females in two ponds in NE Ohio. We also test the prediction that levels of gene flow should be higher for 'sympatric' (sexual males present) genomes in unisexuals compared to 'allopatric' (sexual males absent) unisexual genomes. We used a model testing framework in the coalescent-based program MIGRATE-N to compare models where unidirectional gene flow is present and absent between sexual species and unisexuals. As predicted, our results show higher levels of gene flow between sexuals and sympatric unisexual genomes compared to lower (likely artefactual) levels of gene flow between sexuals and allopatric unisexual genomes. Our results provide direct evidence that genome exchange between sexual and unisexual Ambystoma occurs and demonstrate that the magnitude depends on which sexual species are present. The relatively high levels of gene flow suggest that unisexuals must be at a selective advantage over sexual forms so as to avoid extinction due to genetic swamping through genome exchange.

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.).

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species-specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri-specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure), bidirectional estimates of gene flow (migrate-n and IMa2) and phylogeny-based methods (*beast, bucky) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species-specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri-like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.