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

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

The Lesser Sunda Archipelago offers exceptional potential as a model system for studying the dynamics of dispersal-driven diversification. The geographic proximity of the islands suggests the possibility for successful dispersal, but this is countered by the permanence of the marine barriers and extreme intervening currents that are expected to hinder gene flow. Phylogenetic and species delimitation analyses of flying lizards (genus Draco) using single mitochondrial genes, complete mitochondrial genomes, and exome-capture data sets identified 9-11 deeply divergent lineages including single-island endemics, lineages that span multiple islands, and parapatrically distributed nonsister lineages on the larger islands. Population clustering and PCA confirmed these genetic boundaries with isolation-by-distance playing a role in some islands or island sets. While gdi estimates place most candidate species comparisons in the ambiguous zone, migration estimates suggest 9 or 10 species exist with nuclear introgression detected across some intra-island contact zones. Initial entry of Draco into the archipelago occurred at 5.5-7.5 Ma, with most inter-island colonization events having occurred between 1-3 Ma. Biogeographical model testing favors scenarios integrating geographic distance and historical island connectivity, including an initial stepping-stone dispersal process from the Greater Sunda Shelf through the Sunda Arc as far eastward as Lembata Island. However, rather than reaching the adjacent island of Pantar by dispersing over the 15-km wide Alor Strait, Draco ultimately reached Pantar (and much of the rest of the archipelago) by way of a circuitous route involving at least five overwater dispersal events. These findings suggest that historical geological and oceanographic conditions heavily influenced dispersal pathways and gene flow, which in turn drove species formation and shaped species boundaries. [Biogeography; genomics, Indonesia; lizards; phylogeography; reptiles].

Venom Composition in a Phenotypically Variable Pit Viper ( Trimeresurus insularis) across the Lesser Sunda Archipelago.

The genus Trimeresurus comprises a group of venomous pitvipers endemic to Southeast Asia and the Pacific Islands. Of these, Trimeresurus insularis, the White-lipped Island Pitviper, is a nocturnal, arboreal species that occurs on nearly every major island of the Lesser Sunda archipelago. In the current study, venom phenotypic characteristics of T. insularis sampled from eight Lesser Sunda Islands (Flores, Lembata, Lombok, Pantar, Sumba, Sumbawa, Timor, and Wetar) were evaluated via SDS-PAGE, enzymatic activity assays, fibrinogenolytic assays, gelatin zymography, and RP-HPLC, and the Sumbawa sample was characterized by venomic analysis. For additional comparative analyses, venoms were also examined from several species in the Trimeresurus complex, including T. borneensis, T. gramineus, T. puniceus, T. purpureomaculatus, T. stejnegeri, and Protobothrops flavoviridis. Despite the geographical isolation, T. insularis venoms from all eight islands demonstrated remarkable similarities in gel electrophoretic profiles and RP-HPLC patterns, and all populations had protein bands in the mass ranges of phosphodiesterases (PDE), l-amino acid oxidases (LAAO), P-III snake venom metalloproteinases (SVMP), serine proteases, cysteine-rich secretory proteins (CRISP), phospholipases A2 (PLA2), and C-type lectins. An exception was observed in the Lombok sample, which lacked protein bands in the mass range of serine protease and CRISP. Venomic analysis of the Sumbawa venom also identified these protein families, in addition to several proteins of lesser abundance (<1%), including glutaminyl cyclase, aminopeptidase, PLA2 inhibitor, phospholipase B, cobra venom factor, 5'-nucleotidase, vascular endothelial growth factor, and hyaluronidase. All T. insularis venoms exhibited similarities in thrombin-like and PDE activities, while significant differences were observed for LAAO, SVMP, and kallikrein-like activities, though these differences were only observed for a few islands. Slight but noticeable differences were also observed with fibrinogen and gelatin digestion activities. Trimeresurus insularis venoms exhibited overall similarity to the other Trimeresurus complex species examined, with the exception of P. flavoviridis venom, which showed the greatest overall differentiation. Western blot analysis revealed that all major T. insularis venom proteins were recognized by Green Pitviper ( T. albolabris) antivenom, and reactivity was also seen with most venom proteins of the other Trimeresurus species, but incomplete antivenom-venom recognition was observed against P. flavoviridis venom proteins. These results demonstrate significant conservation in the venom composition of T. insularis across the Lesser Sunda archipelago relative to the other Trimeresurus species examined.

Habitat preference modulates trans-oceanic dispersal in a terrestrial vertebrate.

The importance of long-distance dispersal (LDD) in shaping geographical distributions has been debated since the nineteenth century. In terrestrial vertebrates, LDD events across large water bodies are considered highly improbable, but organismal traits affecting dispersal capacity are generally not taken into account. Here, we focus on a recent lizard radiation and combine a summary-coalescent species tree based on 1225 exons with a probabilistic model that links dispersal capacity to an evolving trait, to investigate whether ecological specialization has influenced the probability of trans-oceanic dispersal. Cryptoblepharus species that occur in coastal habitats have on average dispersed 13 to 14 times more frequently than non-coastal species and coastal specialization has, therefore, led to an extraordinarily widespread distribution that includes multiple continents and distant island archipelagoes. Furthermore, their presence across the Pacific substantially predates the age of human colonization and we can explicitly reject the possibility that these patterns are solely shaped by human-mediated dispersal. Overall, by combining new analytical methods with a comprehensive phylogenomic dataset, we use a quantitative framework to show how coastal specialization can influence dispersal capacity and eventually shape geographical distributions at a macroevolutionary scale.

Toxic toad invasion of Wallacea: A biodiversity hotspot characterized by extraordinary endemism.

(a) A map of Wallacea showing islands invaded by Duttaphrynus melanostictus in red, islands inhabited by Varanus komodoensis in blue, and localities of genetic samples in yellow points. (b) A D. melanostictus from Lombok Island. (c) Environmental niche model for the Sunda Islands clade of D. melanostictus projected into Wallacea. Green color indicates very high suitability, yellow color indicates high suitability, and orange color indicates moderate suitability.

An integrative approach to phylogeography: investigating the effects of ancient seaways, climate, and historical geology on multi-locus phylogeographic boundaries of the Arboreal Salamander (Aneides lugubris).

BACKGROUND: Phylogeography is an important tool that can be used to reveal cryptic biodiversity and to better understand the processes that promote lineage diversification. We studied the phylogeographic history of the Arboreal Salamander (Aneides lugubris), a wide-ranging species endemic to the California floristic province. We used multi-locus data to reconstruct the evolutionary history of A. lugubris and to discover the geographic location of major genetic breaks within the species. We also used species distribution modeling and comparative phylogeography to better understand the environmental factors that have shaped the genetic history of A. lugubris. RESULTS: We found six major mitochondrial clades in A. lugubris. Nuclear loci supported the existence of at least three genetically distinct groups, corresponding to populations north of the San Francisco Bay and in the Sierra Nevada, in the Santa Cruz Mountains, and in the central coast and southern California. All of the genetic breaks in mitochondrial and nuclear loci corresponded to regions where historical barriers to dispersal have been observed in other species. Geologic or water barriers likely were the most important factors restricting gene flow among clades. Climatic unsuitability during glacial maximum may have contributed to the isolation of the mitochondrial clades in the central coast and southern California. A projection of our species distribution model to a future scenario with a moderate amount of climate change suggests that most of the range of A. lugubris will remain climatically suitable, but climatic conditions in the Sierra Nevada and low elevation areas in Southern California are likely to deteriorate. CONCLUSIONS: Aneides lugubris contains substantial cryptic genetic diversity as a result of historical isolation of populations. At least two (and perhaps three) evolutionarily significant units in A. lugubris merit protection; all six mitochondrial clades should be considered as management units within the species.

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.

A new species of terrestrially-nesting fanged frog (Anura: Dicroglossidae) from Sulawesi Island, Indonesia.

Herein, we describe a new species of terrestrially-nesting fanged frog from Sulawesi Island, Indonesia. Though male nest attendance and terrestrial egg deposition is known in one other Sulawesi fanged frog (Limnonectes arathooni), the new species exhibits a derived reproductive mode unique to the Sulawesi assemblage; male frogs guard one or more clutches of eggs festooned to leaves or mossy boulders one to two meters above small slow-moving streams, trickles, or seeps. This island endemic has thus far been collected at three sites on Sulawesi: one in the Central Core of the island, and two on the Southwest Peninsula-south of the Tempe Depression (a major biogeographical boundary). The new Limnonectes has the smallest adult body size among its Sulawesi congeners-with a maximum snout-vent length of about 30 millimeters. Beyond its unique reproductive behavior and body size, the species is further diagnosed on the basis of advertisement call and genetic distance from sympatric fanged frogs. The discovery and description of the new species highlights the remarkable reproductive trait diversity that characterizes the Sulawesi fanged frog assemblage despite that most species in this radiation have yet to be formally described.

Defining evolutionary boundaries across parapatric ecomorphs of Black Salamanders (Aneides flavipunctatus) with conservation implications.

The accurate delimitation of evolutionary population units represents an important component in phylogeographic and conservation genetic studies. Here, we used a combined population assignment and historical demographic approach to study a complex of ecomorphologically distinctive populations of Black Salamanders (Aneides flavipunctatus) that are parapatrically distributed and meet at a three-way contact zone in north-western California. We used mitochondrial tree-based and multilocus clustering methods to evaluate a priori two- (Northern and Southern) and three (Northern, Coast and Inland) population hypotheses derived from previous studies. Mitochondrial results were consistent with the two- and three-population hypotheses, while the nDNA clustering results supported only the two-population hypothesis. Historical demographic analyses and mtDNA gene divergence estimates revealed that the Northern and Southern populations split during the Pliocene (2-5 Ma). Subdivision of the Southern population into Coast and Inland populations was estimated to be late Pleistocene (0.24 Ma), although our mtDNA results suggested a Pliocene divergence. Effective gene flow estimates (2N(e)m) suggest that either the two- or three-population hypotheses remain valid. However, our results unexpectedly revealed that the Northern population might instead represent two parapatric populations that separated nearly 4 Ma. These results are surprising because the Pliocene divergence between these ecomorphologically conservative forms is similar or older than for the ecomorphologically divergent Coast and Inland sister populations. We conclude that Black Salamanders in north-western California belong to at least three or four populations or species, and these all meet criteria for being Evolutionary Significant Units or 'ESUs' and therefore warrant conservation consideration.

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.

Phylogenetic relationships of southern Wallacean ranid frogs (Anura: Ranidae: Hylarana).

Frogs in the family Ranidae are diverse in Asia and are thought to have dispersed to the Sahul Shelf approximately 10 million years ago, where they radiated into more than a dozen species. Ranid species in the intervening oceanic islands of Wallacea, such as Hylarana florensis and H. elberti from the Lesser Sundas and H. moluccana from eastern Wallacea, are assumed to belong to the subgenus Papurana, yet this has not been confirmed with molecular data. We analyzed mitochondrial DNA of Hylarana species from five islands spanning the reported ranges of H. florensis and H. elberti and compared them to confirmed Papurana species and closely related subgenera within Hylarana. We find that the Lesser Sunda H. florensis and H. elberti form a clade that is sister to the rest of the Australo-Papuan Papurana assemblage. Species delimitation analyses and divergence time estimates suggest that populations of H. florensis on Lombok may be distinct from those on Flores at the species level. Likewise, populations of H. elberti on Sumba and Timor may be distinct from each other and from those on Wetar, tshe type locality of H. elberti. Samples from Babar Island thought to be members of H. elberti in fact belong to the wide-ranging H. daemeli, which occurs in northern Australia, across New Guinea, and on the neighboring island of Tanimbar. These results suggest that the Lesser Sundas may have served as a stepping-stone for colonization of the Sahul Shelf and that species diversity of Papurana frogs is underestimated in the Lesser Sundas.

Taxonomic revision of black salamanders of the Aneides flavipunctatus complex (Caudata: Plethodontidae).

We present a taxonomic revision of the black salamander (Aneides flavipunctatus) complex of northwestern California and extreme southeastern Oregon. The revision is based on a number of published works as well as new molecular and morphological data presented herein. The subspecies Aneides flavipunctatus niger Myers & Maslin 1948 is raised in rank to a full species. It is isolated far to the south of the main range on the San Francisco Peninsula, south and west of San Francisco Bay. Another geographically isolated set of populations occurs well inland in Shasta County, northern CA, mainly in the vicinity of Shasta Lake. It is raised from synonymy and recognized as Aneides iecanus (Cope 1883). The remaining taxa occur mainly along and inland from the coast from the vicinity of the Russian River and Lake Berryessa/Putah Creek, north to the vicinity of the Smith River near the Oregon border and more inland along the Klamath and Trinity Rivers and tributaries into Oregon. The northern segment of this nearly continuous range is named Aneides klamathensis Reilly and Wake 2019. We use molecular data to provide a detailed examination of a narrow contact zone between the northern A. klamathensis and the more southern A. flavipunctatus in southern Humboldt County in the vicinity of the Van Duzen and main fork of the Eel rivers. To the south is the remnant of the former species and it takes the name Aneides flavipunctatus (Strauch 1870). It is highly diversified morphologically and genetically and requires additional study.

Squeezing water from a stone: high-throughput sequencing from a 145-year old holotype resolves (barely) a cryptic species problem in flying lizards.

We used Massively Parallel High-Throughput Sequencing to obtain genetic data from a 145-year old holotype specimen of the flying lizard, Draco cristatellus. Obtaining genetic data from this holotype was necessary to resolve an otherwise intractable taxonomic problem involving the status of this species relative to closely related sympatric Draco species that cannot otherwise be distinguished from one another on the basis of museum specimens. Initial analyses suggested that the DNA present in the holotype sample was so degraded as to be unusable for sequencing. However, we used a specialized extraction procedure developed for highly degraded ancient DNA samples and MiSeq shotgun sequencing to obtain just enough low-coverage mitochondrial DNA (721 base pairs) to conclusively resolve the species status of the holotype as well as a second known specimen of this species. The holotype was prepared before the advent of formalin-fixation and therefore was most likely originally fixed with ethanol and never exposed to formalin. Whereas conventional wisdom suggests that formalin-fixed samples should be the most challenging for DNA sequencing, we propose that evaporation during long-term alcohol storage and consequent water-exposure may subject older ethanol-fixed museum specimens to hydrolytic damage. If so, this may pose an even greater challenge for sequencing efforts involving historical samples.

A new species of Lepidodactylus (Reptilia: Squamata: Gekkonidae) from the Kei Islands, Maluku, Indonesia.

Lepidodactylus pantai is a new species of gecko from the Kei Islands, Maluku, Indonesia that is closely associated with intertidal habitats. This species does not fit cleanly into any of the three species groups described for the genus because it possesses the unique combination of both divided terminal scansors on all toes and a nearly completely cylindrical tail without fringes or evidence of dorsoventral compression. A phylogenetic analysis including this species demonstrates that it is the sister taxon of a population from Palau, and that this clade is sister to the clade containing Group III species for which we have molecular data.