Past and Recent Effects of Livestock Activity on the Genetic Diversity and Population Structure of Native Guanaco Populations of Arid Patagonia

History of Diversification and Adaptation from North to South Revealed by Genomic Data: Guanacos from the Desert to Sub-Antarctica

The increased availability of quality genomic data has greatly improved the scope and resolution of our understanding of the recent evolutionary history of wild species adapted to extreme environments and their susceptibility to anthropogenic impacts. The guanaco (Lama guanicoe), the largest wild ungulate in South America, is a good example. The guanaco is well adapted to a wide range of habitats, including the Sechura Desert, the high Andes Mountains to the north, and the extreme temperatures and conditions of Navarino Island to the south. Guanacos also have a long history of overexploitation by humans. To assess the evolutionary impact of these challenging habitats on the genomic diversity, we analyzed 38 genomes (∼10 to 16×) throughout their extensive latitudinal distribution from the Sechura and Atacama Desert to southward into Tierra del Fuego Island. These included analyses of patterns of unique differentiation in the north and geographic region further south with admixture among L. g. cacsilensis and L. g. guanicoe. Our findings provide new insights on the divergence of the subspecies ∼800,000 yr BP and document two divergent demographic trajectories and to the initial expansion of guanaco into the more southern portions of the Atacama Desert. Patagonian guanacos have experienced contemporary reductions in effective population sizes, likely the consequence of anthropogenic impacts. The lowest levels of genetic diversity corresponded to their northern and western limits of distribution and some varying degrees of genetic differentiation. Adaptive genomic diversity was strongly linked with environmental variables and was linked with colonization toward the south followed by adaptation.

Microsatellite-Based Genetic Diversity and Population Structure of Huacaya Alpacas (Vicugna pacos) in Southern Peru

The alpaca population mostly consists of the Huacaya phenotype and is widely distributed in Southern Peru. This study aimed to estimate the genetic diversity and population structure of two Huacaya alpaca populations (Ajoyani and Quimsachata) using fourteen and twelve microsatellite markers for each population, respectively. A total of 168 alpaca biological samples were outsourced to Peruvian laboratories for DNA extraction and genotyping. For genetic diversity, observed heterozygosity (H_o), expected heterozygosity (H_e), polymorphism information content (PIC), and fixation indices values were estimated. An admixture analysis was performed for the population structure analysis. Different programs were used for these estimations. In total, 133 (Ajoyani) and 129 (Quimsachata) alleles were found, with a range of 4 to 17 by locus. The mean H_O, H_E, and PIC per marker for Ajoyani were 0.764 ± 0.112, 0.771 ± 0.1, and 0.736; for Quimsachata, they were 0.783 ± 0.087, 0.773 ± 0.095, and 0.738, respectively. The population structure showed no structure with K = 2. This study provides useful indicators for the creation of appropriate alpaca conservation programs.