Next-generation sequencing (NGS) approaches are increasingly being used to generate multi-locus data for phylogeographic and evolutionary genetics research. We detail the applicability of a restriction enzyme-mediated genome complexity reduction approach with subsequent NGS (DArTseq) in vertebrate study systems at different evolutionary and geographical scales. We present two case studies using SNP data from the DArTseq molecular marker platform. First, we used DArTseq in a large phylogeographic study of the agamid lizard Ctenophorus caudicinctus, including 91 individuals and spanning the geographical range of this species across arid Australia. A low-density DArTseq assay resulted in 28 960 SNPs, with low density referring to a comparably reduced set of identified and sequenced markers as a cost-effective approach. Second, we applied this approach to an evolutionary genetics study of a classic frog hybrid zone (Litoria ewingii–Litoria paraewingi) across 93 individuals, which resulted in 48 117 and 67 060 SNPs for a low- and high-density assay, respectively. We provide a docker-based workflow to facilitate data preparation and analysis, then analyse SNP data using multiple methods including Bayesian model-based clustering and conditional likelihood approaches. Based on comparison of results from the DArTseq platform and traditional molecular approaches, we conclude that DArTseq can be used successfully in vertebrates and will be of particular interest to researchers working at the interface between population genetics and phylogenetics, exploring species boundaries, gene exchange and hybridization.

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 Using the rock‐specialist agamid Ctenophorus caudicinctus as a model, we test hypothesized biogeographical dispersal corridors for lizards in the Australian arid zone (across the western sand deserts), and assess how these dispersal routes have shaped phylogeographical structuring. Location Arid and semi‐arid Australia. Methods We sequenced a c. 1400 bp fragment of mtDNA (ND2) for 134 individuals of C. caudicinctus as well as a subset of each of the mtDNA clades for five nuclear loci (BDNF, BACH1, GAPD, NTF3, and PRLR). We used phylogenetic methods to assess biogeographical patterns within C. caudicinctus, including relaxed molecular clock analyses to estimate divergence times. Ecological niche modelling (Maxent) was employed to estimate the current distribution of suitable climatic envelopes for each lineage. Results Phylogenetic analyses identified two deeply divergent mtDNA clades within C. caudicinctus – an eastern and western clade – separated by the Western Australian sand deserts. However, divergences pre‐date the Pleistocene sand deserts. Phylogenetic analyses of the nuclear DNA data sets generally support major mtDNA clades, suggesting past connections between the western C. c. caudicinctus populations in far eastern Pilbara (EP) and the lineages to the east of the sand deserts. Ecological niche modelling supports the continued suitability of climatic conditions between the Central Ranges and the far EP for C. c. graafi. Main conclusions Estimates of lineage ages provide evidence of divergence between eastern and western clades during the Miocene with subsequent secondary contact during the Pliocene. Our results suggest that this secondary contact occurred via dispersal between the Central Ranges and the far EP, rather than the more southerly Giles Corridor. These events precede the origins of the western sand deserts and divergence patterns instead appear associated with Miocene and Pliocene climate change.

No abstract