Figure Legends
Fig. 1. Photograph of adult Desmognathus aeneus from northern Georgia (courtesy Todd W. Pierson, KSU) and map of 90 sample localities with SDM results of low, medium, and high habitat suitability based on 10 BIOCLIM, ENVIREM, and ecoregion variables (see SI text).
Fig. 2. Linear morphometric analysis of size-corrected phenotypic traits using LDA based on 89 individuals from 31 populations, analysis of which do not generally separate the three geographic genetic lineages. Axilla-Groin Length (AG) was the most substantial variable distinguishing lineages but was not significantly different between them.
Fig. 3. The best-fit three-population demographic model from GADMA estimated using the ‘moments’ engine, showing topology, divergence times, population sizes, and migration rates (a), and estimates of individual ancestry coefficients from sNMF (b, d) along with a PCA of the SNP matrix and the relationship between geographic and genetic distance from the GDM (c).
Fig, 4. Analysis of IBD based on all individuals, between lineages, and within lineages, showing linear (all) and LOESS (a) relationship between geographic distance (km) and linearizedFst .
Figure 5. Results from RDA showing genotype-environment and genotype-phenotype association by sites (a, b) and SNPs (c, d). The 18 SNPs shared in both are outlined in red.