Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus)

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Complex genetic patterns and distribution limits mediated by native congeners of the worldwide invasive red-eared slider turtle.

Non-native (invasive) species offer a unique opportunity to study the geographical distribution and range limits of species, wherein the evolutionary change driven by interspecific interactions between native and non-native closely related species is a key component. The red-eared slider turtle, Trachemys scripta elegans (TSE), has been introduced and successfully established worldwide. It can coexist with its native congeners T. cataspila, T. venusta and T. taylori in Mexico. We performed comprehensive fieldwork, executed a battery of genetic analyses and applied a novel species distribution modelling approach to evaluate their historical lineage relationships and contemporary population genetic patterns. Our findings support the historical common ancestry between native TSE and non-native (TSEalien ), while also highlighting the genetic differentiation of the exotic lineage. Genetic patterns are associated with their range size/endemism gradient; the microendemic T. taylori showed significant reduced genetic diversity and high differentiation, whereas TSEalien showed the highest diversity and signals of population size expansion. Counter to our expectations, lower naturally occurring distribution overlap and little admixture patterns were found between TSE and its congeners, exhibiting reduced gene flow and clear genetic separation across neighbouring species despite having zones of contact. We demonstrate that these native Trachemys species have distinct climatic niche suitability, probably preventing establishment of and displacement by the TSEalien . Moreover, we found major niche overlap between TSEalien and native species worldwide, supporting our prediction that sites with closer ecological optima to the invasive species have higher establishment risk than those that are closer to the niche-centre of the native species.

Interspecific Gene Flow and Mitochondrial Genome Capture during the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae).

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, Anolis opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting (ILS). Here, we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene tree, species tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene tree recovers the same relationships previously inferred for this group, which is strikingly different from the species tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model (MSCM), so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group; however, coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent, adaptive radiations. [Adaptive radiation; hybridization; introgression; multispecies network coalescent; posterior predictive simulation.].

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome.

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.