Rapid Radiation and Rampant Reticulation: Phylogenomics of South American Liolaemus Lizards.

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction site-associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.[Andes; ddRADSeq; introgression; lizards; mtDNA; reptiles; SNPs.].

Polyploidy breaks speciation barriers in Australian burrowing frogs Neobatrachus.

Polyploidy has played an important role in evolution across the tree of life but it is still unclear how polyploid lineages may persist after their initial formation. While both common and well-studied in plants, polyploidy is rare in animals and generally less understood. The Australian burrowing frog genus Neobatrachus is comprised of six diploid and three polyploid species and offers a powerful animal polyploid model system. We generated exome-capture sequence data from 87 individuals representing all nine species of Neobatrachus to investigate species-level relationships, the origin and inheritance mode of polyploid species, and the population genomic effects of polyploidy on genus-wide demography. We describe rapid speciation of diploid Neobatrachus species and show that the three independently originated polyploid species have tetrasomic or mixed inheritance. We document higher genetic diversity in tetraploids, resulting from widespread gene flow between the tetraploids, asymmetric inter-ploidy gene flow directed from sympatric diploids to tetraploids, and isolation of diploid species from each other. We also constructed models of ecologically suitable areas for each species to investigate the impact of climate on differing ploidy levels. These models suggest substantial change in suitable areas compared to past climate, which correspond to population genomic estimates of demographic histories. We propose that Neobatrachus diploids may be suffering the early genomic impacts of climate-induced habitat loss, while tetraploids appear to be avoiding this fate, possibly due to widespread gene flow. Finally, we demonstrate that Neobatrachus is an attractive model to study the effects of ploidy on the evolution of adaptation in animals.

How mountains shape biodiversity: The role of the Andes in biogeography, diversification, and reproductive biology in South America's most species-rich lizard radiation (Squamata: Liolaemidae).

Testing hypotheses on drivers of clade evolution and trait diversification provides insight into many aspects of evolutionary biology. Often, studies investigate only intrinsic biological properties of organisms as the causes of diversity, however, extrinsic properties of a clade's environment, particularly geological history, may also offer compelling explanations. The Andes are a young mountain chain known to have shaped many aspects of climate and diversity of South America. The Liolaemidae are a radiation of South American reptiles with over 300 species found across most biomes and with similar numbers of egg-laying and live-bearing species. Using the most complete dated phylogeny of the family, we tested the role of Andean uplift in biogeography, diversification patterns, and parity mode of the Liolaemidae. We find that the Andes promoted lineage diversification and acted as a species pump into surrounding biomes. We also find strong support for the role of Andean uplift in boosting the species diversity of these lizards via allopatric fragmentation. Finally, we find repeated shifts in parity mode associated with changing thermal niches, with live-bearing favored in cold climates and egg-laying favored in warm climates. Importantly, we find evidence for possible reversals to oviparity, an evolutionary transition believed to be extremely rare.

Positive Darwinian selection results in resistance to cardioactive toxins in true toads (Anura: Bufonidae).

Members of the Family Bufonidae, true toads, are famous for their endogenously synthesized cardioactive steroids that serve as defensive toxins. Evolution of resistance to these toxins is not understood. We sequenced a key region of the toxin's binding site in the Na(+)/K(+) ATPase for relevant taxa representing Hyloidea (including bufonids), Ranoidea and Archaeobatrachia and tested for positive selection in a phylogenetic context. Bufonidae were distinct from other Hyloidea at 4-6 of 12 sites and, with one exception, had a homologous amino acid sequence. Melanophryniscus stelzneri had a distinct sequence, consistent with other independent evidence for a differentiated toxin. Tests within Bufonidae detected positive selection within the binding region, providing, to our knowledge, the first evidence of this type for positive selection within Amphibia. There was no evidence for positive selection on Bufonidae or M. stelzneri lineages. Sequence change in Leptodactylus ocellatus, a leptodactylid predator of Bufonidae, provides a molecular basis for predator resistance possibly associated with gene duplication.

Molecular phylogeny and divergence dates for Australasian elapids and sea snakes (hydrophiinae): evidence from seven genes for rapid evolutionary radiations.

One of the most prolific radiations of venomous snakes, the Australo-Melanesian Hydrophiinae includes approximately 100 species of Australasian terrestrial elapids plus all approximately 60 species of viviparous sea snakes. Here, we estimate hydrophiine relationships based on a large data set comprising 5800 bp drawn from seven genes (mitochondrial: ND4, cytb, 12S, 16S; nuclear: rag1, cmos, myh). These data were analysed using parsimony, likelihood and Bayesian methods to better resolve hydrophiine phylogeny and provide a timescale for the terrestrial and marine radiations. Among oviparous forms, Cacophis, Furina and Demansia are basal to other Australian elapids (core oxyuranines). The Melanesian Toxicocalamus and Aspidomorphus group with Demansia, indicating multiple dispersal events between New Guinea and Australia. Oxyuranus and Pseudonaja form a robust clade. The small burrowing taxa form two separate clades, one consisting of Vermicella and Neelaps calanotus, and the other including Simoselaps, Brachyurophis and Neelaps bimaculatus. The viviparous terrestrial elapids form three separate groups: Acanthophis, the Rhinoplocephalus group and the Notechis-Hemiaspis group. True sea snakes (Hydrophiini) are robustly united with the Notechis-Hemiaspis group. Many of the retrieved groupings are consistent with previous molecular and morphological analyses, but the polyphyly of the viviparous and burrowing groups, and of Neelaps, are novel results. Bayesian relaxed clock analyses indicate very recent divergences: the approximately 160 species of the core Australian radiation (including sea snakes) arose within the last 10 Myr, with most inter-generic splits dating to between 10 and 6 Ma. The Hydrophis sea snake lineage is an exceptionally rapid radiation, with > 40 species evolving within the last 5 Myr.