Population Genomics of the Foothill Yellow-Legged Frog (Rana boylii) and RADseq Parameter Choice for Large-Genome Organisms

Genomic data are useful for attaining high resolution in population genetic studies and have become increasingly available for answering questions in biological conservation. We analyzed RADseq data for the protected foothill yellow-legged frog (Rana boylii) throughout its native range in California and Oregon, including many of the same localities included in an earlier study based on mitochondrial DNA. We recovered five primary clades that correspond to geographic regions within California and Oregon, with better resolution and more spatially consistent patterns than the previous study, confirming the increased resolving power of genomic approaches compared to single-locus analyses. Bayesian clustering, PCA and population differentiation with admixture analyses all indicated that approximately half the range of R. boylii consists of a single, relatively uniform population, while regions in the Sierra Nevada and Central Coast Range of California are deeply differentiated genetically. Additionally, a major methodological challenge for large genome organisms, including many amphibians, is deciding on sequence similarity clustering thresholds for population genetic analyses using RADseq data, and we develop a novel set of metrics that allow researchers to set a sequence similarity threshold that maximizes the separation of paralogous regions while minimizing the oversplitting of naturally occurring allelic variation within loci.