Genetic network analysis uncovers spatial variation in diversity and connectivity of a species presenting a continuous distribution

The conservation of genetic diversity and connectivity is essential for the long-term persistence and adaptive ability of a species. Recent calls have been made for the inclusion of genetic diversity and differentiation measures in the assessment, management, and conservation of species. However, the literature often lacks direction on how to do so for species with continuous distributions or no distinct breaks in genetic connectivity. There are many considerations to overcome when investigating genetic diversity and connectivity of such species. We combine multiple genetic network methodologies with more traditional population genetic analyses within a single framework to address the challenges of investigating population structure and quantifying variation in genetic diversity and connectivity of wide-ranging species with continuous distributions. We demonstrate the efficacy and applicability of our framework through a study on woodland caribou (Rangifer tarandus) occupying the boreal forest of Canada; a species of significant conservation concern. The dataset consisted of 4911 unique individuals genotyped at 9 microsatellite loci, which were subsequently partitioned into 103 spatial nodes to create a population-based genetic network. The Walktrap community detection algorithm was used to detect hierarchical population genetic structure within the study area and node-based network metrics such as mean inverse edge weight and clustering coefficient were used to quantify the variation in genetic connectivity across the range. Lastly, genetic diversity was assessed by calculating allelic richness and heterozygosity of the nodes making up the network. The community detection analysis identified two communities at the coarsest scale to nine communities at the optimal partition. A strong pattern of Isolation by Distance (IBD) was found across the range at multiple scales. Furthermore, signs of genetic erosion along the study areas southern boundaries were depicted by nodes presenting low genetic diversity and low centrality values. These results are important to the species status assessments in providing previously unavailable information on connectivity and diversity within and beyond the current local population units used in management. Our approach to quantify the patterns and extent of connectivity across the boreal range is comprehensive and could easily be adapted to other species. The results are robust and provide a solid foundation for the continued monitoring and recovery of the species.