Phylogenomic Analysis Reveals Dispersal-Driven Speciation and Divergence with Gene Flow in Lesser Sunda Flying Lizards (Genus Draco)

To and fro in the archipelago: Repeated inter-island dispersal and New Guinea's orogeny affect diversification of Delias, the world's largest butterfly genus

The world's largest butterfly genus Delias, commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15-16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene-likely facilitated by frequent sea level change-sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification.

Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group)

Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.

Integrative methods reveal multiple drivers of diversification in rice paddy snakes

Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae: Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis that Hypsiscopus diversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggest Hypsiscopus diversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species of Hypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets.

The natural and human-mediated expansion of a human-commensal lizard into the fringes of Southeast Asia

BACKGROUND

Human-commensal species often display deep ancestral genetic structure within their native range and founder-effects and/or evidence of multiple introductions and admixture in newly established areas. We investigated the phylogeography of Eutropis multifasciata, an abundant human-commensal scincid lizard that occurs across Southeast Asia, to determine the extent of its native range and to assess the sources and signatures of human introduction outside of the native range. We sequenced over 350 samples of E. multifasciata for the mitochondrial ND2 gene and reanalyzed a previous RADseq population genetic dataset in a phylogenetic framework.

RESULTS

Nuclear and mitochondrial trees are concordant and show that E. multifasciata has retained high levels of genetic structure across Southeast Asia despite being frequently moved by humans. Lineage boundaries in the native range roughly correspond to several major biogeographic barriers, including Wallace's Line and the Isthmus of Kra. Islands at the outer fringe of the range show evidence of founder-effects and multiple introductions.

CONCLUSIONS

Most of enormous range of E. multifasciata across Southeast Asia is native and it only displays signs of human-introduction or recent expansion along the eastern and northern fringe of its range. There were at least three events of human-introductions to Taiwan and offshore islands, and several oceanic islands in eastern Indonesia show a similar pattern. In Myanmar and Hainan, there is a founder-effect consistent with post-warming expansion after the last glacial maxima or human introduction.

Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

Bewildering biogeography: waves of dispersal and diversification across southern Wallacea by bent-toed geckos (genus: Cyrtodactylus)

Bent-toed Geckos, genus Cyrtodactylus, are one of the most diverse terrestrial vertebrate groups, and their range extends from South Asia into Australo-Papua and adjacent Pacific islands. Given the generally high faunal endemism on Wallacean islands, it is rather paradoxical that the diversity in these geckos appears to be so low (21 species in Wallacea, 15 in the Philippines) compared with continental shelf assemblages (>300 species on Sunda + Sahul Shelves + adjacent islands). To determine whether this shortfall was real or an artifact of historical undersampling, we analyzed mitochondrial DNA sequences of hundreds of southern Wallacean samples (Lesser Sundas + southern Maluku). After screening to guide sample selection for target capture data collection, we obtained a 1150-locus genomic dataset (1,476,505 bp) for 119 samples of southern Wallacean and closely related lineages. The results suggest that species diversity of Cyrtodactylus in southern Wallacea is vastly underestimated, with phylogenomic and clustering analyses suggesting as many as 25 candidate species, in contrast to the 8 currently described. Gene exchange between adjacent candidate species is absent or minimal across the archipelago with only one case of >0.5 migrants per generation. Biogeographical analysis suggests that the hitherto unrecognized diversity is the result of at least three independent dispersals from Sulawesi or its offshore islands into southern Wallacea between 6-14 Ma, with one invasion producing small-bodied geckos and the other two or three producing larger-bodied geckos. The smaller-bodied laevigatus group appears to be able to coexist with members of either larger-bodied clade, but we have yet to find members of the two larger-bodied clades occurring in sympatry, suggesting that ecological partitioning or competitive exclusion may be shaping individual island assemblages.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Diverge and Conquer: Phylogenomics of southern Wallacean forest skinks (Genus: Sphenomorphus) and their colonization of the Lesser Sunda Archipelago

The archipelagos of Wallacea extend between the Sunda and Sahul Shelves, serving as a semipermeable two-way filter influencing faunal exchange between Asia and Australo-Papua. Forest skinks (Genus Sphenomorphus) are widespread throughout southern Wallacea and exhibit complex clinal, ontogenetic, sexual, and seasonal morphological variation, rendering species delimitation difficult. We screened a mitochondrial marker for 245 Sphenomorphus specimens from this area to inform the selection of 104 samples from which we used targeted sequence capture to generate a dataset of 1154 nuclear genes (∼1.8 Mb) plus complete mitochondrial genomes. Phylogenomic analyses recovered many deeply divergent lineages, three pairs of which are now sympatric, that began to diversify in the late Miocene shortly after the oldest islands are thought to have become emergent. We infer a complex and nonstepping-stone pattern of island colonization, with the group having originated in the Sunda Arc islands before using Sumba as a springboard for colonization of the Banda Arcs. Estimates of population structure and gene flow across the region suggest total isolation except between two Pleistocene Aggregate Island Complexes that become episodically land-bridged during glacial maxima. These historical processes have resulted in at least 11 Sphenomorphus species in the region, nine of which require formal description. This fine-scale geographic partitioning of undescribed species highlights the importance of utilizing comprehensive genomic studies for defining biodiversity hotspots to be considered for conservation protection.