The genetic signature left by the range expansion of red foxes in Australia is detectable after more than 80 years of population stability.

Reconstructing biological invasions from historical sources can provide insights into how they occur but are difficult to do when invasions are poorly documented. Genetic signatures left by invaders can also offer insights into invasion routes, points of origin and general biology but often present conclusions that are contradictory to expectations. Here, we test the ability of continental-wide microsatellite genotype data from 29 loci and 3122 samples to reconstruct the well-documented invasion of red foxes Vulpes vulpes from the United Kingdom into Australia over 150 years ago, an invasion that has led to the extinction of many native species. Our analysis reveals several key signals of invasion evident in Australian foxes. They display lower levels of diversity than foxes sampled from the UK, exhibit clines in diversity from the point of introduction (south-east Australia) to the edge of their range, and show strong evidence of allele surfing in westerly and north-easterly directions. These characteristics are consistent with a single point of origin followed by rapid expansion in westerly and north-easterly directions as suggested by historical records. We also find little genetic structure in foxes across Australia with only the vast Nullarbor Plains and Great Victoria Desert region presenting a detectable barrier to their dispersal. As such, no mainland region within the current range of foxes can be considered genetically isolated and therefore appropriate for localised eradication efforts. Overall, our analyses demonstrate the ability of comprehensive population genetic studies to reconstruct invasion histories even after more than 80 years since colonisation was stabilised.

dartr: An r package to facilitate analysis of SNP data generated from reduced representation genome sequencing.

Although vast technological advances have been made and genetic software packages are growing in number, it is not a trivial task to analyse SNP data. We announce a new r package, dartr, enabling the analysis of single nucleotide polymorphism data for population genomic and phylogenomic applications. dartr provides user-friendly functions for data quality control and marker selection, and permits rigorous evaluations of conformation to Hardy-Weinberg equilibrium, gametic-phase disequilibrium and neutrality. The package reports standard descriptive statistics, permits exploration of patterns in the data through principal components analysis and conducts standard F-statistics, as well as basic phylogenetic analyses, population assignment, isolation by distance and exports data to a variety of commonly used downstream applications (e.g., newhybrids, faststructure and phylogeny applications) outside of the r environment. The package serves two main purposes: first, a user-friendly approach to lower the hurdle to analyse such data-therefore, the package comes with a detailed tutorial targeted to the r beginner to allow data analysis without requiring deep knowledge of r. Second, we use a single, well-established format-genlight from the adegenet package-as input for all our functions to avoid data reformatting. By strictly using the genlight format, we hope to facilitate this format as the de facto standard of future software developments and hence reduce the format jungle of genetic data sets. The dartr package is available via the r CRAN network and GitHub.