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.].

Phylogenomic evidence for a recent and rapid radiation of lizards in the Patagonian Liolaemus fitzingerii species group.

Rapid evolutionary radiations are difficult to resolve because divergence events are nearly synchronous and gene flow among nascent species can be high, resulting in a phylogenetic "bush". Large datasets composed of sequence loci from across the genome can potentially help resolve some of these difficult phylogenetic problems. A suitable test case is the Liolaemus fitzingerii species group of lizards, which includes twelve species that are broadly distributed in Argentinean Patagonia. The species in the group have had a complex evolutionary history that has led to high morphological variation and unstable taxonomy. We generated a sequence capture dataset for 28 ingroup individuals of 580 nuclear loci, alongside a mitogenomic dataset, to infer phylogenetic relationships among species in this group. Relationships among species were generally weakly supported with the nuclear data, and along with an inferred age of ∼2.6 million years old, indicate either rapid evolution, hybridization, incomplete lineage sorting, non-informative data, or a combination thereof. We inferred a signal of mito-nuclear discordance, indicating potential hybridization between L. melanops and L. martorii, and phylogenetic network analyses provided support for 5 reticulation events among species. Phasing the nuclear loci did not provide additional insight into relationships or suspected patterns of hybridization. Only one clade, composed of L. camarones, L. fitzingerii, and L. xanthoviridis was recovered across all analyses. Genomic datasets provide molecular systematists with new opportunities to resolve difficult phylogenetic problems, yet the lack of phylogenetic resolution in Patagonian Liolaemus is biologically meaningful and indicative of a recent and rapid evolutionary radiation. The phylogenetic relationships of the Liolaemus fitzingerii group may be best modeled as a reticulated network instead of a bifurcating phylogeny.

Cryptic diversity and morphological convergence in threatened species of fossorial skinks in the genus Scelotes (Squamata: Scincidae) from the Western Cape Coast of South Africa: implications for species boundaries, digit reduction and conservation.

We investigated the evolutionary relationships among populations of two threatened Red Data Book fossorial skinks, Scelotes gronovii and Scelotes kasneri, along the Western Cape Coast of South Africa. The genus Scelotes shows considerable variation in limb and digit reduction. We sampled four localities purported to contain S. gronovii and seven of S. kasneri, encompassing all of each species' limited distribution. Each of these species lack forelimbs, and differ by the number of digits on the hind limbs, among other morphological characters; S. gronovii bears a single digit and S. kasneri bears two digits on the hind limbs. Sequence data obtained from three mtDNA (16S ribosomal RNA, cytochrome b, and nicotinamide adenine dinucleotide dehydrogenase 1 unit; 2035 bp ttl.) and two nuclear (dynein axonemal heavy chain 3 and the natural killer tumor recognition; 1848 bp ttl.) gene regions were used to reconstruct the evolutionary relationships among the two focal species and several other co-distributed species (Scelotes bipes, Scelotes montispectus, and Scelotes sexlineatus). Phylogenetic results (Bayesian and parsimony) revealed that several populations previously considered S. kasneri actually belong to other species, and others are paraphyletic with respect to one another. Additionally, populations of S. gronovii were also found to be paraphyletic, with populations south of the Berg River supported as sister to S. bipes, and populations north of the Berg River sister the remaining sampled species. Our results require a redefinition of S. sexlineatus to encompass populations morphologically convergent with S. kasneri and restrict the ranges of the already threatened S. kasneri and S. gronovii even further. The paraphyly of S. gronovii and the placement of each clade as sister to clades of species bearing two digits on the hind limbs suggests that digit loss has occurred at least twice in this group.

Phylogeny of iguanian lizards inferred from 29 nuclear loci, and a comparison of concatenated and species-tree approaches for an ancient, rapid radiation.

Iguanian lizards form a diverse clade whose members have been the focus of many comparative studies of ecology, behavior, and evolution. Despite the importance of phylogeny to such studies, interrelationships among many iguanian clades remain uncertain. Within the Old World clade Acrodonta, Agamidae is sometimes found to be paraphyletic with respect to Chamaeleonidae, and recent molecular studies have produced conflicting results for many major clades. Within the largely New World clade Pleurodonta, relationships among the 12 currently recognized major subclades (mostly ranked as families) have been largely unresolved or poorly supported in previous studies. To clarify iguanian evolutionary history, we first infer phylogenies using concatenated maximum-likelihood (ML) and Bayesian analyses of DNA sequence data from 29 nuclear protein-coding genes for 47 iguanian and 29 outgroup taxa. We then estimate a relaxed-clock Bayesian chronogram for iguanians using BEAST. All three methods produce identical topologies. Within Acrodonta, we find strong support for monophyly of Agamidae with respect to Chamaeleonidae, and for almost all relationships within agamids. Within Pleurodonta, we find strong Bayesian support for almost all relationships, and strong ML support for some interfamilial relationships and for monophyly of almost all families (excepting Polychrotidae). Our phylogenetic results suggest a non-traditional biogeographic scenario in which pleurodonts originated in the Northern Hemisphere and subsequently spread southward into South America. The pleurodont portion of the tree is characterized by several very short, deep branches, raising the possibility of deep coalescences that may confound concatenated analyses. We therefore also use 27 of these genes to implement a coalescent-based species-tree approach for pleurodonts. Although this analysis strongly supports monophyly of the pleurodont families, interfamilial relationships are generally different from those in the concatenated tree, and support is uniformly poor. However, a species-tree analysis using only the seven most variable loci yields higher support and more congruence with the concatenated tree. This suggests that low support in the 27-gene species-tree analysis may be an artifact of the many loci that are uninformative for very short branches. This may be a general problem for the application of species-tree methods to rapid radiations, even with phylogenomic data sets. Finally, we correct the non-monophyly of Polychrotidae by recognizing the pleurodont genus Anolis (sensu lato) as a separate family (Dactyloidae), and we correct the non-monophyly of the agamid genus Physignathus by resurrection of the genus Istiurus for the former Physignathus lesueurii.

Phylogenetic relationships and limb loss in sub-Saharan African scincine lizards (Squamata: Scincidae).

Skinks are the largest family of lizards and are found worldwide in a diversity of habitats. One of the larger and more poorly studied groups of skinks includes members of the subfamily Scincinae distributed in sub-Saharan Africa. Sub-Saharan African scincines are one of the many groups of lizards that show limb reduction and loss, and the genus Scelotes offers an excellent opportunity to look at limb loss in a phylogenetic context. Phylogenetic relationships were reconstructed for a total of 52 taxa representing all subfamilies of skinks as well as other Autarchoglossan families using sequence from six gene regions including; 12S, 16S, and cytochrome b (mitochondrial), as well as alpha-Enolase, 18S, and C-mos (nuclear). The family Scincidae is recovered as monophyletic and is the sister taxon to a (Cordylidae+Xantusiidae) clade. Within skinks the subfamily Acontinae is monophyletic and sister group to all remaining skinks. There is no support for the monophyly of the subfamilies Lygosominae and Scincinae, but sub-Saharan African scincines+Feylinia form a well supported monophyletic group. The monophyly of Scelotes is confirmed, and support is found for two geographic groups within the genus. Reconstructions of ancestral states for limb and digital characters show limited support for the reversal or gain of both digits and limbs, but conservative interpretation of the results suggest that limb loss is common, occurring multiple times throughout evolutionary history, and is most likely not reversible.