Y-chromosome and mtDNA variation confirms independent domestications and directional hybridization in South American camelids

Whole genome sequencing analysis of alpaca suggests TRPV3 as a candidate gene for the suri phenotype

BACKGROUND

Alpaca is a domestic South American camelid probably arising from the domestication of two wild camelids, the vicugna and the guanaco. Two phenotypes are described for alpaca, known as huacaya and suri. Huacaya fleece is characterized by compact, soft, and highly crimped fibers, while suri fleece is longer, straight, less crimped, and lustrous. The gene variants determining these phenotypes are still unknown, although previous studies suggested a dominant inheritance of the suri. Based on that, the aim of this study was the identification of the gene variants determining alpaca coat phenotypes through whole genome sequencing (WGS) analysis.

RESULTS

The sample used includes two test-cross alpaca families, suri × huacaya, which produced two offspring, one with the suri phenotype and one with the huacaya phenotype. The analyzed sample was expanded through the addition of WGS data from six vicugnas and six guanacos; this because we assumed the absence of the gene variants linked to the suri phenotype in these wild species. The analysis of gene variant segregation with the suri phenotype, coupled with the filtering of gene variants present in the wild species, disclosed the presence in all the suri samples of a premature termination codon (PTC) in TRPV3 (transient receptor potential cation channel subfamily V member 3), a gene known to be involved in hair growth and cycling, thermal sensation, cold tolerance and adaptation in several species. Mutations in TRPV3 were previously associated with the alteration of hair structure leading to an impaired formation of the hair canal and the hair shaft in mouse. This PTC in TRPV3, due to a G > T substitution (p.Glu475*), results in a loss of 290 amino acids from the canonical translated protein, plausibly leading to a physiological dysfunction.

CONCLUSION

The present results suggest that the suri phenotype may arise from a TRPV3 gene variant which may explain some of the suri features such as its longer hair fibre with lower number of cuticular scales compared to huacaya.

Genome-wide scan for runs of homozygosity in South American Camelids

BACKGROUND

Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection.

RESULTS

A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively.

CONCLUSIONS

The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

Genotyping-by-sequencing (GBS) as a tool for interspecies hybrid detection

Genotyping-by-sequencing (GBS) is an extremely useful, modern and relatively inexpensive approach to discovering high-quality single-nucleotide polymorphisms (SNPs), which seem to be the most promising markers for identifying hybrid individuals between different species, especially those that can create backcrosses. In addition, GBS could become an invaluable tool in finding backcrosses, even several generations back. Its potential for the use of restriction enzymes and species is almost unlimited. It can also be successfully applied to species for which a reference genome is not established. In this paper, we describe the GBS technique, its main advantages and disadvantages, and the research carried out using this method concerning interspecies hybridisation and the identification of fertile hybrids. We also present future approaches that could be of interest in the context of the GBS method.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

Microsatellite-Based Genetic Structure and Hybrid Detection in Alpacas Bred in Poland

This study aimed to characterize the population structure and genetic diversity of alpacas maintained in Poland using 17 microsatellite markers recommended by the International Society for Animal Genetics. The classification of llamas, alpacas, and hybrids of both based on phenotype is often difficult due to long-term admixture. Our results showed that microsatellite markers can distinguish alpacas from llamas and provide information about the level of admixture of one species in another. Alpacas admixed with llamas constituted 8.8% of the tested individuals, with the first-generation hybrid displaying only 7.4% of llama admixture. The results showed that Poland hosts a high alpaca genetic diversity as a consequence of their mixed origin. More than 200 different alleles were identified and the average observed heterozygosity and expected heterozygosity values were 0.745 and 0.768, respectively, the average coefficient of inbreeding was 0.034, and the average polymorphism information content value was 0.741. The probability of exclusion for one parent was estimated at 0.99995 and for two parents at 0.99999.

Ancient DNA reveals the lost domestication history of South American camelids in Northern Chile and across the Andes

The study of South American camelids and their domestication is a highly debated topic in zooarchaeology. Identifying the domestic species (alpaca and llama) in archaeological sites based solely on morphological data is challenging due to their similarity with respect to their wild ancestors. Using genetic methods also presents challenges due to the hybridization history of the domestic species, which are thought to have extensively hybridized following the Spanish conquest of South America that resulted in camelids slaughtered en masse. In this study, we generated mitochondrial genomes for 61 ancient South American camelids dated between 3,500 and 2,400 years before the present (Early Formative period) from two archaeological sites in Northern Chile (Tulán-54 and Tulán-85), as well as 66 modern camelid mitogenomes and 815 modern mitochondrial control region sequences from across South America. In addition, we performed osteometric analyses to differentiate big and small body size camelids. A comparative analysis of these data suggests that a substantial proportion of the ancient vicuña genetic variation has been lost since the Early Formative period, as it is not present in modern specimens. Moreover, we propose a domestication hypothesis that includes an ancient guanaco population that no longer exists. Finally, we find evidence that interbreeding practices were widespread during the domestication process by the early camelid herders in the Atacama during the Early Formative period and predating the Spanish conquest.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Were domestic camelids present on the prehispanic South American agricultural frontier? An ancient DNA study

The southern boundary of prehispanic farming in South America occurs in central Mendoza Province, Argentina at approximately 34 degrees south latitude. Archaeological evidence of farming includes the recovery of macrobotanical remains of cultigens and isotopic chemistry of human bone. Since the 1990s, archaeologists have also hypothesized that the llama (Lama glama), a domesticated South American camelid, was also herded near the southern boundary of prehispanic farming. The remains of a wild congeneric camelid, the guanaco (Lama guanicoe), however, are common in archaeological sites throughout Mendoza Province. It is difficult to distinguish bones of the domestic llama from wild guanaco in terms of osteological morphology, and therefore, claims that llama were in geographic areas where guanaco were also present based on osteometric analysis alone remain equivocal. A recent study, for example, claimed that twenty-five percent of the camelid remains from the high elevation Andes site of Laguna del Diamante S4 were identified based on osteometric evidence as domestic llama, but guanaco are also a likely candidate since the two species overlap in size. We test the hypothesis that domesticated camelids occurred in prehispanic, southern Mendoza through analysis of ancient DNA. We generated whole mitochondrial genome datasets from 41 samples from southern Mendoza late Holocene archaeological sites, located between 450 and 3400 meters above sea level (masl). All camelid samples from those sites were identified as guanaco; thus, we have no evidence to support the hypothesis that the domestic llama occurred in prehispanic southern Mendoza.

Untangling the fibre ball: Proteomic characterization of South American camelid hair fibres by untargeted multivariate analysis and molecular networking

The proteomic analysis of hairs, yarns or textiles has emerged as a powerful method to determine species of origin, mainly used in archaeozoological research and fraud control. Differentiation between the South American camelid (SAC) species (the wild guanaco and vicuña and their respective domesticates the llama and alpaca) is particularly challenging due to poor database information and significant hybridization between species. In this study, we analysed 41 modern and 4 archaeological samples from the four SACs species. Despite strong similarities with Old World Camelidae, we identified 7 peptides specific to SACs assigned to keratin K86 and the keratin-associated proteins KAP13-1 and KAP11-1. Untargeted multivariate analysis of the LC-MS data permitted to distinguish SAC species and propose discriminant features. MS/MS-based molecular networking combined with database-assisted de novo sequencing permitted to identify 5 new taxonomic peptides assigned to K33a, K81 and/or K83 keratins and KAP19-1. These peptides differentiate the two wild species, guanaco and vicuña. These results show the value of combining database search and untargeted metabolomic approaches for paleoproteomics, and reveal for the first time the potential of molecular networks to highlight deamidation related to diagenesis and cluster highly similar peptides related to interchain homologies or intra- or inter-specific polymorphism. SIGNIFICANCE: This study used an innovative approach combining multivariate analysis of LC-MS data together with molecular networking and database-assisted de novo sequencing to identify taxonomic peptides in palaeoproteomics. It constitutes the first attempt to differentiate between hair fibres from the four South American camelids (SACs) based on proteomic analysis of modern and archaeological samples. It provides different proteomic signatures for each of the four SAC species and proposes new SAC taxonomic peptides of interest in archaeozoology and fraud control. SACs have been extensively exploited since human colonization of South America but have not been studied to the extent of their economic, cultural and heritage importance. Applied to the analysis of ancient Andean textiles, our results should permit a better understanding of cultural and pastoral practices in South America. The wild SACs are endangered by poaching and black-market sale of their fibre. For the first time, our results provide discriminant features for the determination of species of origin of contraband fibre.

South American Camelids: their values and contributions to people

South American Camelids (SACs) make several material and non-material contributions to people and are a key component of the Andean biocultural heritage. From the perspective of the IPBES' Conceptual Framework, SACs constitute the "nature" component in the complex system of interactions between human societies and the Andean mountain environment. There are four SAC living species today, two of which are wild, or Salqa, in the indigenous cosmovision: guanaco (Lama guanicoe) and vicuña (Vicugna vicugna). Llama (Lama glama) and alpaca (Vicugna pacos) were domesticated 5000 years ago, and are therefore Uywa, in the indigenous cosmovision. Both wild and domestic camelids were, and in several cases still are, the most highly appreciated resource for Andean livelihoods. Historically, camelids and their contributions have been used by Andean people since the peopling of the Americas over 11,000 years ago. In this paper, we present three case studies (chakus for vicuña management, llama caravans, and llama nanobodies) to bring attention to the essential role of vicuñas and llamas for Andean communities today, their intercultural linkages with the Western world, and telecoupling interactions.

Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca

BACKGROUND

Despite their regional economic importance and being increasingly reared globally, the origins and evolution of the llama and alpaca remain poorly understood. Here we report reference genomes for the llama, and for the guanaco and vicuña (their putative wild progenitors), compare these with the published alpaca genome, and resequence seven individuals of all four species to better understand domestication and introgression between the llama and alpaca.

RESULTS

Phylogenomic analysis confirms that the llama was domesticated from the guanaco and the alpaca from the vicuña. Introgression was much higher in the alpaca genome (36%) than the llama (5%) and could be dated close to the time of the Spanish conquest, approximately 500 years ago. Introgression patterns are at their most variable on the X-chromosome of the alpaca, featuring 53 genes known to have deleterious X-linked phenotypes in humans. Strong genome-wide introgression signatures include olfactory receptor complexes into both species, hypertension resistance into alpaca, and fleece/fiber traits into llama. Genomic signatures of domestication in the llama include male reproductive traits, while in alpaca feature fleece characteristics, olfaction-related and hypoxia adaptation traits. Expression analysis of the introgressed region that is syntenic to human HSA4q21, a gene cluster previously associated with hypertension in humans under hypoxic conditions, shows a previously undocumented role for PRDM8 downregulation as a potential transcriptional regulation mechanism, analogous to that previously reported at high altitude for hypoxia-inducible factor 1α.

CONCLUSIONS

The unprecedented introgression signatures within both domestic camelid genomes may reflect post-conquest changes in agriculture and the breakdown of traditional management practices.

Evaluation of the short- and long-term complications and outcomes of phacoemulsification surgery in alpacas

OBJECTIVE

To describe the complications and determine the visual outcomes associated with phacoemulsification surgery in alpacas.

ANIMALS

Five alpacas (nine operated eyes).

PROCEDURES

Medical records of alpacas with a clinical diagnosis of cataracts that received elective phacoemulsification surgery at Cornell University over a 15 year period were reviewed.

RESULTS

The mean (± standard deviation) age of the five operated alpacas was 104 days (±64 days). The cataracts were suspected to be congenital in all alpacas. Concurrent congenital ocular abnormalities were common in the operated eyes and included microphakia, persistent patent hyaloid arteries, and persistent patent tunica vasculosa lentis. No significant intraoperative complications were identified in any cases. In the immediate postoperative period, uveitis was mild-to-moderate in most operated eyes, ocular hypertension was not detected, and clinically significant corneal edema did not occur in any eyes. The mean (± standard deviation) duration of clinical follow-up after surgery was 32 months (±25 months). A peri-incisional corneal stromal ulcer associated with fibrinous uveitis was identified in a single eye during recheck examination 4 weeks postoperative that resolved with medical therapy. Long-term complications included unilateral glaucoma in two alpacas that resulted in vision loss. On the final recheck ophthalmic examination, 78% (7 of 9) of operated eyes were visual and all operated alpacas had vision in at least one operated eye.

CONCLUSIONS AND CLINICAL RELEVANCE

In this series, the visual outcomes associated with elective phacoemulsification surgery in alpacas were good. Glaucoma was the most common vision-threatening long-term postoperative complication.

Utility of genetic variation in coat color genes to distinguish wild, domestic and hybrid South American camelids for forensic and judicial applications

A molecular genetic protocol for distinguishing pure and hybrid South American camelids was developed to provide strong, quantifiable, and unbiased species identification. We detail the application of the approach in the context of a criminal case in the Andes Mountains of central Chile where the defendants were alleged to have illegally hunted three wild guanacos (Lama guanicoe), as opposed to hybrid domestic llama (Lama glama)/wild guanaco crosses, which are unregulated. We describe a workflow that differentiates among wild, domestic and hybrid South American camelids (Lama versus Vicugna) based on mitochondrial cytochrome b genetic variation (to distinguish between Lama and Vicugna), and MC1R and exon 4 variation of the ASIP gene (to differentiate wild from domestic species). Additionally, we infer the population origin and sex of each of the three individuals from a panel of 15 autosomal microsatellite loci and the presence or absence of the SRY gene. Our analyses strongly supported the inference that the confiscated carcasses corresponded with 2 male and 1 female guanacos that were hunted illegally. Statistical power analyses suggested that there was an extremely low probability of misidentifying domestic camelids as wild camelids (an estimated 0 % Type I error rate), or using more conservative approached a 1.17 % chance of misidentification of wild species as domestic camelids (Type II error). Our case report and methodological and analytical protocols demonstrate the power of genetic variation in coat color genes to identify hybrids between wild and domestic camelid species and highlight the utility of the approach to help combat illegal wildlife hunting and trafficking.

Phylogeography and Population Genetics of Vicugna vicugna: Evolution in the Arid Andean High Plateau

The vicuña (Vicugna vicugna) is the most representative wild ungulate of the high Andes of South America with two recognized morphological subspecies, V. v. mensalis in the north and V. v. vicugna in the south of its distribution. Current vicuña population size (460,000-520,000 animals) is the result of population recovery programs established in response to 500 years of overexploitation. Despite the vicuña's ecosystemic, economic and social importance, studies about their genetic variation and history are limited and geographically restricted. Here, we present a comprehensive assessment of the genetic diversity of vicuña based on samples collected throughout its distribution range corresponding to eleven localities in Peru and five in Chile representing V. v. mensalis, plus four localities each in Argentina and Chile representing V. v. vicugna. Analysis of mitochondrial DNA and microsatellite markers show contrasting results regarding differentiation between the two vicuña types with mitochondrial haplotypes supporting subspecies differentiation, albeit with only a few mutational steps separating the two subspecies. In contrast, microsatellite markers show that vicuña genetic variation is best explained as an isolation by distance pattern where populations on opposite ends of the distribution present different allelic compositions, but the intermediate populations present a variety of alleles shared by both extreme forms. Demographic characterization of the species evidenced a simultaneous and strong reduction in the effective population size in all localities supporting the existence of a unique, large ancestral population (effective size ∼50,000 individuals) as recently as the mid-Holocene. Furthermore, the genetic variation observed across all localities is better explained by a model of gene flow interconnecting them rather than only by genetic drift. Consequently, we propose space "continuous" Management Units for vicuña as populations exhibit differentiation by distance and spatial autocorrelation linked to sex biased dispersal instead of population fragmentation or geographical barriers across the distribution.

Energy expenditure and body temperature variations in llamas living in the High Andes of Peru

Some large herbivores exhibit seasonal adjustments in their energy metabolism. Therefore, our aim was to determine if the llama (one of the most extensively kept livestock breeds) exhibits seasonal adjustment of their energy expenditure, body temperature and locomotion, under its natural high altitude Andean habitat. For this purpose, energy expenditure, body temperature and locomotion were measured in seven non-pregnant llama dams for ten months on the Andean High Plateau (4400 m above sea level). Daily energy expenditure was measured as field metabolic rate using the doubly labelled water method at four different measurement times. Additionally, a telemetry system was used to continuously record activity, body temperature (3 min intervals) as well as the position (hourly) of each animal. The results show that llamas adjusted their body temperature and daily energy expenditure according to environmental conditions. Furthermore, llamas under high altitude Andean climatic conditions exhibited a pronounced daily rhythm in body temperature and activity, with low values at sunrise and increasing values towards sunset. Llamas also had remarkably low energy expenditure compared to other herbivores. Thus, despite the domestication process, llamas have not lost the ability to adjust their body temperature and daily energy expenditure under adverse environmental conditions, similar to some wild herbivores.

Genetic Variation in Coat Colour Genes MC1R and ASIP Provides Insights Into Domestication and Management of South American Camelids

The domestication of wild vicuña and guanaco by early pre-Inca cultures is an iconic example of wildlife management and domestication in the Americas. Although domestic llamas and alpacas were clearly selected for key, yet distinct, phenotypic traits, the relative patterns and direction of selection and domestication have not been confirmed using genetic approaches. However, the detailed archaeological records from the region suggest that domestication was a process carried out under significant control and planning, which would have facilitated coordinated and thus extremely effective selective pressure to achieve and maintain desired phenotypic traits. Here we link patterns of sequence variation in two well-characterised genes coding for colour variation in vertebrates and interpret the results in the context of domestication in guanacos and vicuñas. We hypothesise that colour variation in wild populations of guanacos and vicunas were strongly selected against. In contrast, variation in coat colour variation in alpaca was strongly selected for and became rapidly fixed in alpacas. In contrast, coat colour variants in llamas were of less economic value, and thus were under less selective pressure. We report for the first time the full sequence of MC1R and 3 exons of ASIP in 171 wild specimens from throughout their distribution and which represented a range of commonly observed colour patterns. We found a significant difference in the number of non-synonymous substitutions, but not synonymous substitutions among wild and domestics species. The genetic variation in MC1R and ASIP did not differentiate alpaca from llama due to the high degree of reciprocal introgression, but the combination of 11 substitutions are sufficient to distinguish domestic from wild animals. Although there is gene flow among domestic and wild species, most of the non-synonymous variation in MC1R and ASIP was not observed in wild species, presumably because these substitutions and the associated colour phenotypes are not effectively transmitted back into wild populations. Therefore, this set of substitutions unequivocally differentiates wild from domestic animals, which will have important practical application in forensic cases involving the poaching of wild vicuñas and guanacos. These markers will also assist in identifying and studying archaeological remains pre- and post-domestication.

Change of niche in guanaco (Lama guanicoe): the effects of climate change on habitat suitability and lineage conservatism in Chile

Background

The main goal of this contribution was to define the ecological niche of the guanaco (Lama guanicoe), to describe potential distributional changes, and to assess the relative importance of niche conservatism and divergence processes between the two lineages described for the species (L.g. cacsilensis and L.g. guanicoe).

Methods

We used maximum entropy to model lineage’s climate niche from 3,321 locations throughout continental Chile, and developed future niche models under climate change for two extreme greenhouse gas emission scenarios (RCP2.6 and RCP8.5). We evaluated changes of the environmental niche and future distribution of the largest mammal in the Southern Cone of South America. Evaluation of niche conservatism and divergence were based on identity and background similarity tests.

Results

We show that: (a) the current geographic distribution of lineages is associated with different climatic requirements that are related to the geographic areas where these lineages are located; (b) future distribution models predict a decrease in the distribution surface under both scenarios; (c) a 3% decrease of areal protection is expected if the current distribution of protected areas is maintained, and this is expected to occur at the expense of a large reduction of high quality habitats under the best scenario; (d) current and future distribution ranges of guanaco mostly adhere to phylogenetic niche divergence hypotheses between lineages.

Discussion

Associating environmental variables with species ecological niche seems to be an important aspect of unveiling the particularities of, both evolutionary patterns and ecological features that species face in a changing environment. We report specific descriptions of how these patterns may play out under the most extreme climate change predictions and provide a grim outlook of the future potential distribution of guanaco in Chile. From an ecological perspective, while a slightly smaller distribution area is expected, this may come with an important reduction of available quality habitats. From the evolutionary perspective, we describe the limitations of this taxon as it experiences forces imposed by climate change dynamics.