No abstract

Our knowledge of the conservation status of reptiles, the most diverse class of terrestrial vertebrates, has improved dramatically over the past decade, but still lags behind that of the other tetrapod groups. Here, we conduct the first comprehensive evaluation (~92% of the world's ~1714 described species) of the conservation 1 Joint senior authors. D.G. Chapple et al.

Aim The Lesser Sunda Islands are situated between the Sunda and Sahul Shelves, with a linear arrangement that has functioned as a two‐way filter for taxa dispersing between the Asian and Australo‐Papuan biogeographical realms. Distributional patterns of many terrestrial vertebrates suggest a stepping‐stone model of island colonization. Here we investigate the timing and sequence of island colonization in Asian‐origin fanged frogs from the volcanic Sunda Arc islands with the goal of testing the stepping‐stone model of island colonization. Location The Indonesian islands of Java, Lombok, Sumbawa, Flores and Lembata. Taxon Limnonectes dammermani and L. kadarsani (Family: Dicroglossidae) Methods Mitochondrial DNA was sequenced from 153 frogs to identify major lineages and to select samples for an exon‐capture experiment. We designed probes to capture sequence data from 974 exonic loci (1,235,981 bp) from 48 frogs including the outgroup species, L. microdiscus. The resulting data were analysed using phylogenetic, population genetic and biogeographical model testing methods. Results The mtDNA phylogeny finds L. kadarsani paraphyletic with respect to L. dammermani, with a pectinate topology consistent with the stepping‐stone model. Phylogenomic analyses of 974 exons recovered the two species as monophyletic sister taxa that diverged ~7.6 Ma with no detectable contemporary gene flow, suggesting introgression of the L. dammermani mitochondrion into L. kadarsani on Lombok resulting from an isolated ancient hybridization event ~4 Ma. Within L. kadarsani, the Lombok lineage diverged first while the Sumbawa and Lembata lineages are nested within a Flores assemblage composed of two parapatrically distributed lineages meeting in central Flores. Biogeographical model comparison found strict stepping‐stone dispersal to be less likely than models involving leap‐frog dispersal events. Main conclusions These results suggest that the currently accepted stepping‐stone model of island colonization might not best explain the current patterns of diversity in the archipelago. The high degree of genetic structure, large divergence times, and absent or low levels of migration between lineages suggests that L. kadarsani represents five distinct species.

Marker selection has emerged as an important component of phylogenomic study design due to rising concerns of the effects of gene tree estimation error, model misspecification, and data-type differences. Researchers must balance various trade-offs associated with locus length and evolutionary rate among other factors. The most commonly used reduced representation data sets for phylogenomics are ultraconserved elements (UCEs) and Anchored Hybrid Enrichment (AHE). Here, we introduce Rapidly Evolving Long Exon Capture (RELEC), a new set of loci that targets single exons that are both rapidly evolving (evolutionary rate faster than RAG1) and relatively long in length (>1,500 bp), while at the same time avoiding paralogy issues across amniotes. We compare the RELEC data set to UCEs and AHE in squamate reptiles by aligning and analyzing orthologous sequences from 17 squamate genomes, composed of 10 snakes and 7 lizards. The RELEC data set (179 loci) outperforms AHE and UCEs by maximizing per-locus genetic variation while maintaining presence and orthology across a range of evolutionary scales. RELEC markers show higher phylogenetic informativeness than UCE and AHE loci, and RELEC gene trees show greater similarity to the species tree than AHE or UCE gene trees. Furthermore, with fewer loci, RELEC remains computationally tractable for full Bayesian coalescent species tree analyses. We contrast RELEC to and discuss important aspects of comparable methods, and demonstrate how RELEC may be the most effective set of loci for resolving difficult nodes and rapid radiations. We provide several resources for capturing or extracting RELEC loci from other amniote groups.

The gekkonid genus Cyrtodactylus is a highly diverse group of lizards (280 + species), which covers an expansive geographic range. Although this genus has been the focus of many taxonomic and molecular systematic studies, species on the Southeast Asian island of Borneo have remained understudied, leading to an unclear evolutionary history with cascading effects on taxonomy and biogeographic inferences. We assembled the most comprehensive multilocus Bornean dataset (one mitochondrial and three nuclear loci) that included 129 novel sequences and representatives from each known Cyrtodactylus species on the island to validate taxonomic status, assess species diversity, and elucidate biogeographic patterns. Our results uncovered a high proportion of cryptic diversity and revealed numerous taxonomic complications, especially within the C. consobrinus, C. malayanus, and C. pubisulcus groups. Comparisons of pairwise genetic distances and a preliminary species delimitation analysis using the Automatic Barcode Gap Discovery (ABGD) method demonstrated that some wide-ranging species on Borneo likely comprise multiple distinct and deeply divergent lineages, each with more restricted distributional ranges. We also tested the prevailing biogeographic hypothesis of a single invasion from Borneo into the Philippines. Our analyses revealed that Philippine taxa were not monophyletic, but were likely derived from multiple separate invasions into the geopolitical areas comprising the Philippines. Although our investigation of Bornean Cyrtodactylus is the most comprehensive to-date, it highlights the need for expanded taxonomic sampling and suggests that our knowledge of the evolutionary history, systematics, and biogeography of Bornean Cyrtodactylus is far from complete.

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.

Episodic sea level changes that repeatedly exposed and inundated the Sunda Shelf characterize the Pleistocene. Available evidence points to a more xeric central Sunda Shelf during periods of low sea levels, and despite the broad land connections that persisted during this time, some organisms are assumed to have faced barriers to dispersal between land-masses on the Sunda Shelf. Eutropis rugifera is a secretive, forest adapted scincid lizard that ranges across the Sunda Shelf. In this study, we sequenced one mitochondrial (ND2) and four nuclear (BRCA1, BRCA2, RAG1, and MC1R) markers and generated a time-calibrated phylogeny in BEAST to test whether divergence times between Sundaic populations of E. rugifera occurred during Pleistocene sea-level changes, or if they predate the Pleistocene. We find that E. rugifera shows pre-Pleistocene divergences between populations on different Sundaic land-masses. The earliest divergence within E. rugifera separates the Philippine samples from the Sundaic samples approximately 16 Ma; the Philippine populations thus cannot be considered conspecific with Sundaic congeners. Sundaic populations diverged approximately 6 Ma, and populations within Borneo from Sabah and Sarawak separated approximately 4.5 Ma in the early Pliocene, followed by further cladogenesis in Sarawak through the Pleistocene. Divergence of peninsular Malaysian populations from the Mentawai Archipelago occurred approximately 5 Ma. Separation among island populations from the Mentawai Archipelago likely dates to the Pliocene/Pleistocene boundary approximately 3.5 Ma, and our samples from peninsular Malaysia appear to coalesce in the middle Pleistocene, about 1 Ma. Coupled with the monophyly of these populations, these divergence times suggest that despite consistent land-connections between these regions throughout the Pleistocene E. rugifera still faced barriers to dispersal, which may be a result of environmental shifts that accompanied the sea-level changes.

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 non-sister 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 over-water 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.