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

Sexual Dichromatism Drives Diversification within a Major Radiation of African Amphibians.

Theory predicts that sexually dimorphic traits under strong sexual selection, particularly those involved with intersexual signaling, can accelerate speciation and produce bursts of diversification. Sexual dichromatism (sexual dimorphism in color) is widely used as a proxy for sexual selection and is associated with rapid diversification in several animal groups, yet studies using phylogenetic comparative methods to explicitly test for an association between sexual dichromatism and diversification have produced conflicting results. Sexual dichromatism is rare in frogs, but it is both striking and prevalent in African reed frogs, a major component of the diverse frog radiation termed Afrobatrachia. In contrast to most other vertebrates, reed frogs display female-biased dichromatism in which females undergo color transformation, often resulting in more ornate coloration in females than in males. We produce a robust phylogeny of Afrobatrachia to investigate the evolutionary origins of sexual dichromatism in this radiation and examine whether the presence of dichromatism is associated with increased rates of net diversification. We find that sexual dichromatism evolved once within hyperoliids and was followed by numerous independent reversals to monochromatism. We detect significant diversification rate heterogeneity in Afrobatrachia and find that sexually dichromatic lineages have double the average net diversification rate of monochromatic lineages. By conducting trait simulations on our empirical phylogeny, we demonstrate that our inference of trait-dependent diversification is robust. Although sexual dichromatism in hyperoliid frogs is linked to their rapid diversification and supports macroevolutionary predictions of speciation by sexual selection, the function of dichromatism in reed frogs remains unclear. We propose that reed frogs are a compelling system for studying the roles of natural and sexual selection on the evolution of sexual dichromatism across micro- and macroevolutionary timescales.

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.

Discordance between genomic divergence and phenotypic variation in a rapidly evolving avian genus (Motacilla).

Generally, genotypes and phenotypes are expected to be spatially congruent; however, in widespread species complexes with few barriers to dispersal, multiple contact zones, and limited reproductive isolation, discordance between phenotypes and phylogeographic groups is more probable. Wagtails (Motacilla) are a genus of birds with striking plumage pattern variation across the Old World. Up to 13 subspecies are recognized within a single species, yet previous studies using mitochondrial DNA have supported polyphyletic phylogeographic groups that are inconsistent with subspecies plumage characteristics. In this study, we investigate the link between phenotypes and genotype by taking a phylogenetic approach. We use genome-wide SNPs, nuclear introns, and mitochondrial DNA to estimate population structure, isolation by distance, and species relationships. Together, our genetic sampling includes complete species-level sampling and comprehensive coverage of the three most phenotypically diverse Palearctic species. Our study provides strong evidence for species-level patterns of differentiation, however population-level differentiation is less pronounced. SNPs provide a robust estimate of species-level relationships, which are mostly corroborated by a combined analysis of mtDNA and nuclear introns (the first time-calibrated species tree for the genus). However, the mtDNA tree is strongly incongruent and is considered to misrepresent the species phylogeny. The extant wagtail lineages originated during the Pliocene and the Eurasian lineage underwent rapid diversification during the Pleistocene. Three of four widespread Eurasian species exhibit an east-west divide that contradicts both subspecies taxonomy and phenotypic variation. Indeed, SNPs fail to distinguish between phenotypically distinct subspecies within the M. alba and M. flava complexes, and instead support geographical regions, each of which is home to two or more different looking subspecies. This is a major step towards our understanding of wagtail phylogeny compared to previous analyses of fewer species and considerably less sequence data.

Phylogenomic analysis of the Chilean clade of Liolaemus lizards (Squamata: Liolaemidae) based on sequence capture data.

The genus Liolaemus is one of the most ecologically diverse and species-rich genera of lizards worldwide. It currently includes more than 250 recognized species, which have been subject to many ecological and evolutionary studies. Nevertheless, Liolaemus lizards have a complex taxonomic history, mainly due to the incongruence between morphological and genetic data, incomplete taxon sampling, incomplete lineage sorting and hybridization. In addition, as many species have restricted and remote distributions, this has hampered their examination and inclusion in molecular systematic studies. The aims of this study are to infer a robust phylogeny for a subsample of lizards representing the Chilean clade (subgenus Liolaemus sensu stricto), and to test the monophyly of several of the major species groups. We use a phylogenomic approach, targeting 541 ultra-conserved elements (UCEs) and 44 protein-coding genes for 16 taxa. We conduct a comparison of phylogenetic analyses using maximum-likelihood and several species tree inference methods. The UCEs provide stronger support for phylogenetic relationships compared to the protein-coding genes; however, the UCEs outnumber the protein-coding genes by 10-fold. On average, the protein-coding genes contain over twice the number of informative sites. Based on our phylogenomic analyses, all the groups sampled are polyphyletic. Liolaemus tenuis tenuis is difficult to place in the phylogeny, because only a few loci (nine) were recovered for this species. Topologies or support values did not change dramatically upon exclusion of L. t. tenuis from analyses, suggesting that missing data did not had a significant impact on phylogenetic inference in this data set. The phylogenomic analyses provide strong support for sister group relationships between L. fuscus, L. monticola, L. nigroviridis and L. nitidus, and L. platei and L. velosoi. Despite our limited taxon sampling, we have provided a reliable starting hypothesis for the relationships among many major groups of the Chilean clade of Liolaemus that will help future work aimed at resolving the Liolaemus phylogeny.

Phylogenomics of a rapid radiation: is chromosomal evolution linked to increased diversification in north american spiny lizards (Genus Sceloporus)?

BACKGROUND: Resolving the short phylogenetic branches that result from rapid evolutionary diversification often requires large numbers of loci. We collected targeted sequence capture data from 585 nuclear loci (541 ultraconserved elements and 44 protein-coding genes) to estimate the phylogenetic relationships among iguanian lizards in the North American genus Sceloporus. We tested for diversification rate shifts to determine if rapid radiation in the genus is correlated with chromosomal evolution. RESULTS: The phylogenomic trees that we obtained for Sceloporus using concatenation and coalescent-based species tree inference provide strong support for the monophyly and interrelationships among nearly all major groups. The diversification analysis supported one rate shift on the Sceloporus phylogeny approximately 20-25 million years ago that is associated with the doubling of the speciation rate from 0.06 species/million years (Ma) to 0.15 species/Ma. The posterior probability for this rate shift occurring on the branch leading to the Sceloporus species groups exhibiting increased chromosomal diversity is high (posterior probability = 0.997). CONCLUSIONS: Despite high levels of gene tree discordance, we were able to estimate a phylogenomic tree for Sceloporus that solves some of the taxonomic problems caused by previous analyses of fewer loci. The taxonomic changes that we propose using this new phylogenomic tree help clarify the number and composition of the major species groups in the genus. Our study provides new evidence for a putative link between chromosomal evolution and the rapid divergence and radiation of Sceloporus across North America.