Next Article in Journal
Complex Analysis of Retroposed Genes’ Contribution to Human Genome, Proteome and Transcriptome
Next Article in Special Issue
Genome-Wide Diversity, Population Structure and Demographic History of Dromedaries in the Central Desert of Iran
Previous Article in Journal
Genome-Wide Transcriptome and Metabolome Analyses Provide Novel Insights and Suggest a Sex-Specific Response to Heat Stress in Pigs
Previous Article in Special Issue
Cytogenetic Mapping of 35 New Markers in the Alpaca (Vicugna pacos)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity and Population Structure of Llamas (Lama glama) from the Camelid Germplasm Bank—Quimsachata

1
Laboratorio de Biología Molecular y Genómica, Dirección de Recursos Genéticos y Biotecnología, Instituto Nacional de Innovación Agraria (INIA). Av. La Molina 1981, 15024 Lima, Peru
2
Facultad de Zootecnia, Universidad Nacional Agraria La Molina (UNALM), Av. La Molina s/n, 15024 Lima, Peru
*
Authors to whom correspondence should be addressed.
Current address: Environmental Cell Biology Group, Department of Functional and Evolutionary Ecology, University of Vienna, 1010 Vienna, Austria.
Genes 2020, 11(5), 541; https://doi.org/10.3390/genes11050541
Submission received: 6 April 2020 / Revised: 30 April 2020 / Accepted: 4 May 2020 / Published: 12 May 2020
(This article belongs to the Special Issue Camelid Genomics)

Abstract

:
Llamas (Lama glama) are invaluable resources of Peru. Despite their importance, their population is decreasing. The Camelid Germplasm Bank—Quimsachata was created as a guardian of this South American camelid (SAC) species and established a bank of llamas from their two types, Ch’aku and Q’ara. However, these populations need to present high genetic diversity to be considered suitable conservation stocks. Thus, in the present study, 13 microsatellites specific for the SAC were used to assess the current genetic variability and differentiation of the llama population from the Bank. The global population showed high genetic diversity with a total of 157 different alleles, with an average of 12.08 alleles per microsatellite, an expected and observed heterozygosity of 0.758 and 0.707, respectively, and an average polymorphic information content (PIC) of 0.723. Although considered as two different breeds and managed separately, the genetic differentiation between Ch’aku and Q’ara was low (FST = 0.01). Accordingly, the gene flow value was high (Nm = 30.5). Overall, our results indicate the existence of high genetic variation among individuals, and thus, this llama population could be considered a suitable genetic stock for their conservation and for sustainability programs. Additionally, the 13 microsatellites can be used to study other Peruvian llama populations and monitor the genetic variability of llamas from the Camelid Germplasm Bank—Quimsachata.

1. Introduction

Llamas (Lama glama) are the largest South American camelids (SACs) and the best adapted to a wide range of environmental conditions [1]. They can be mainly found in the Andean region between 2300 and 4000 m above sea level in Peru, Bolivia, Ecuador, northwest Argentina, and central Chile [2]. After Bolivia, Peru is the country with the second-largest population of llamas in the world [3], being primarily found in the Department of Puno, with almost 35% of the total Peruvian llama population [4]. For the local economy, llamas are seen as multipurpose animals; given their low-fat and low-cholesterol but high-protein content, their meat is consumed [2], and occasionally their intestines are used to make string and drums, while their excrement is used as fuel. Moreover, they are pack animals and their fiber is often used for clothing [5]. Two main types of llamas are recognized: Q’ara and Ch’aku (or T’hampulli). The Q’ara is a light-wool type with a long and slim body and has lower quality fiber, but it possesses a greater aptitude for meat production [6,7,8]; the heavy-wool type, Ch’aku, has a shorter body but has a higher potential for fiber production [9].
Despite being used for different utilities, the population of llamas is decreasing, with around 300,000 fewer llamas reported in the last national census [10] compared to the one in 1994 [11,12]. The diminishing llama population in Peru is due to various factors, such as the intense selection of the white fiber of alpacas, which leads to the decrease of colored (non-white) alpacas and llamas in order to increase the rearing of more white alpacas [13]. Furthermore, llama keepers tend to be neglected, as project funding greatly benefits the production system of other SAC species, mainly alpacas [12,14].
Therefore, before the imminent loss of their phenotypic and genetic diversity, in 1987 the National Institute of Agricultural Innovation (INIA) established the “Camelid Germplasm Bank—Quimsachata” in Puno, Peru, oriented to the recovery and conservation of colored alpacas and the two llama types, Ch’aku and Q’ara [13]. However, to be considered adequate genetic stocks, these populations must present high genetic diversity to guarantee the conservation and augment the variability and productivity of these animals [15]. However, to date, there are no genetic characterization studies of the llama population from the Bank that could help to develop appropriate management strategies, detect rare alleles that indicate the presence of unique genetic variants [16], and avoid inbreeding depression [17].
Microsatellite markers (hereafter referred to as microsatellites), due to their high polymorphism and abundant distribution throughout the genome [18], are important tools for the assessment of genetic diversity and the structure of livestock populations [19,20,21]. There are over 70 nuclear microsatellites reported for the SAC [22,23,24,25,26,27,28,29] that can be potentially used for the genetic characterization of SAC populations—a first step in making conservation strategies. In Peru, most studies (some involving microsatellites) have been carried out in alpaca and vicuña populations [21,30,31,32], whereas in llama populations these were either focused on studying their origin [29,33,34] or on analyzing quantitative parameters related to the physiology, production, growth, or fiber of the animal [35,36,37,38].
Therefore, this is the first population genetics study employing microsatellites carried out on the llama population of the Camelid Germplasm Bank—Quimsachata, whose biodiversity represents the natural, economic, cultural, and historical patrimony of Peru, and it aims to (i) assess its genetic diversity and structure and (ii) estimate the genetic differentiation between the two llama subpopulations, Ch’aku and Q’ara.

2. Materials and Methods

2.1. Sample Collection

Blood samples from 251 adult llamas of types Ch’aku (n = 92, 67 females and 25 males) and Q’ara (n = 159, 109 females and 50 males) were collected at the Camelid Germplasm Bank—Quimsachata, located in the district of Santa Lucia in the Department of Puno, Peru (Figure 1). The protocol used for the blood sample collection agreed with the requirements of the National Law No. 30407 “Ley de Protección y Bienestar Animal (Animal Protection and Welfare Law)”. To measure the genetic differentiation between both phenotypes, Ch’aku and Q’ara llamas were considered as two different subpopulations. Llamas were selected after analyzing their pedigree records, and only unrelated animals were sampled. The difference in the number of animals between the two populations and in the number of males and females that were enrolled in the study were due to the fact that the Camelid Germplasm Bank—Quimsachata has more Q’ara than Ch’aku phenotypes and a male to female ratio of approximately 30 to 70. Thus, the sampling was designed to keep this ratio and to have a representative sample of the population. Total genomic DNA was extracted using a standard phenol–chloroform and ethanol precipitation protocol [39]. The DNA pellet was resuspended in TE buffer and stored at −20 °C.

2.2. DNA Amplification and Microsatellite Genotyping

Thirteen microsatellites specific to llamas and alpacas, namely LCA82, LCA54, LCA65, LCA83, LCA77, LCA85, YWLL08, YWLL44, YWLL59, LAB1, GLM4, Lgu76, and VOLP03 (Table S1), were chosen due to their high polymorphic information content (PIC), elevated heterozygosity, and genetic diversity [22,24,25,26,28,40] and were genotyped across all samples. Genotyping was performed following the procedure described by De Arruda et al. (2010) [41]. The forward primer carried an extension sequence of 19 bp attached to its 5′ end (M13 sequence), which allowed labeling with three different fluorescent dyes (6-FAM, NED, and HEX; Table S1). The reverse primer remained unaltered.
Genomic DNA was amplified by polymerase chain reaction (PCR) using a Mastercycler Pro S (Eppendorf, Hauppauge, NY, USA). Microsatellites with similar PCR conditions were co-amplified using a multiplex PCR, and the fluorescent labeling of the forward primers allowed for the design of multiloading panels. Each PCR run was performed in a total volume of 10 μL, containing 50 ng of DNA. Reaction mixtures contained 1x PCR buffer, 2 mM of MgCl2, 0.2 μM of each primer, 0.2 mM of dNTPs, and 0.5U of Taq polymerase. The amplification conditions included an initial denaturation step of 95 °C for 5 min, followed by 35 cycles of 95 °C for 30 s, 90 s at 56 °C or 58 °C, 72 °C for 1 min, and a final extension at 72 °C for 30 min. The PCR products were separated by capillary electrophoresis in an automatic ABI Prism 3130XL Genetic Analyzer® (Applied Biosystems, Foster City, CA, USA). Genotyping was performed using GeneMapper v.4.0 software (Applied Biosystems).

2.3. Statistical Analysis

The most used genetic diversity parameters, such as the number of alleles per microsatellite, allelic frequencies, mean number of alleles per microsatellite, private alleles, expected and observed heterozygosity (He and Ho, respectively), and polymorphic information content (PIC), were calculated for each of the 13 microsatellite markers, using the CERVUS 3.0.3 [42] and GENETIX 4.0.5 [43] software. Possible deviations from the Hardy–Weinberg equilibrium (HWE), either due to an excess or to a deficit of heterozygous in the total population and within subpopulations, were estimated using Fisher’s exact test implemented in the GENEPOP 4.0.11 software [44]. Furthermore, to guarantee the quality of the results, the null allele frequencies were also calculated using Micro-Checker 2.2.3 [45]. The level of genetic differentiation among individuals, within and between subpopulations, was calculated by the analysis of molecular variance (AMOVA) test using the ARLEQUIN 3.1 software [46]. The extent of genetic differentiation between the two llama subpopulations, Ch’aku and Q’ara, was quantified using the F-statistics (FIS, FIT, and FST; [47]) using GENEPOP 4.0.11 and corroborated with FSTAT 2.9.3.2 [48]. The effective number of migrants per generation (Nm, [19]) was estimated by GENEPOP 4.0.11.
The genetic structure was determined through a grouping analysis of the individuals in a different number of inferred clusters (K) using an analysis based on the ‘admixture’ ancestry model implemented in the STRUCTURE 2.3 software [49]. The burn-in period was set to 50,000 followed by 500,000 Markov chain Monte Carlo (MCMC) iterations. Independent runs of K were performed from 1 to 7 clusters and were repeated 4 times to check the consistency of the results. Finally, a factorial correspondence analysis was performed with GENETIX 3.0.3 to further assess the genetic relationships between the llama types, describing the association of qualitative variables in which each individual is represented just once for the value of each modality (microsatellites) and variable (alleles per microsatellite).

3. Results

3.1. Genetic Diversity Assessment

The 13 microsatellites used in this study were polymorphic and revealed a high level of genetic diversity in the population of llamas from the Camelid Germplasm Bank—Quimsachata. The number of alleles, He, Ho, and PIC values for each microsatellite for Ch’aku and Q’ara as separated populations and as a global population are shown in Table 1.
In total, 157 alleles were found across the 251 individuals, with an average of 12.08 alleles per microsatellite (ranging from 8 for LCA54, LCA83, and LCA85 up to 19 for YWLL08 and YWLL59). The He value was higher than the Ho value, with an average of 0.707 in the global population. Additionally, the He and Ho values were always higher in the Ch’aku subpopulation compared to the Q’ara. The PIC value ranged from 0.517 (LCA54) to 0.883 (YWLL08).
For the global population, 7 out of the 13 microsatellites departed from the HWE (p < 0.05). The inbreeding coefficient (FIS) was estimated for each microsatellite for the global population and for the Ch’aku and Q’ara subpopulations, and they were positive and significantly different from zero (p > 0.05) (Table 2). Moreover, we observed a deficit of heterozygotes and 7.1% more homozygotes than would be expected under the HWE (Table 2). The existence of null alleles was assessed with Micro-Checker, and only the microsatellites LCA82A, YWLL59A, LCA85A, and GLM4 showed signs of having them (Table S2). On the other hand, a total of 35 private alleles were observed (15 for Ch’aku and 20 for Q’ara subpopulations), but their frequencies were very low (<0.05) (Table 3).

3.2. Genetic Differentiation between Ch’aku and Q’ara

The genetic differentiation coefficient value (FST) between the Ch’aku and Q’ara types showed very low genetic differentiation (FST = 0.01) (Table 4), whereas a high gene flow value was observed (Nm = 30.9 migrants per generation). Additionally, the analysis of molecular variance indicated that the highest variance was due to variation within the populations (93.95%), while only 1.02% variance was observed between the two subpopulations (Table S3). The FST value calculated by Arlequin 3.1 was 0.0102 and coincided with the value given by FSTAT 2.9.3.2 [48].
Bayesian analysis carried out by the STRUCTURE software [49], using seven independent runs (K = 1, 2, 3, 4, 5, 6, 7) and each one repeated four times, did not show any evidence of genetic differentiation or population subdivisions. The highest likelihood was obtained when K = 2 and the individuals of both phenotypes were assigned to two clusters, but there was not a clear separation between the two llama subpopulations (Figure 2). Thus, there was no evidence of population structure in the Camelid Germplasm Bank—Quimsachata. Additionally, the correspondence factorial analysis showed a graphic representation of the genetic relationship between the Ch’aku and Q’ara subpopulations (Figure 3). This analysis indicated once again the low genetic differentiation between the subpopulations, given that both overlap with each other and do not form clear independent groups.

4. Discussion

There has been an upsurge in attempts to conserve natural resources, since it is well known that the loss of genetic variability diminishes the ability to recover endangered species [50] and decreases the chance to improve the performance of animals involved in production systems [51]. In this scenario, the Camelid Germplasm Bank—Quimsachata plays a crucial role in conserving a natural, economic, cultural, and historical resource of the SAC populations in Peru. The population of llamas in this Bank must have high genetic diversity to be considered an appropriate genetic stock that could contribute to ensuring its conservation, aid in implementing future strategies to face the loss of diversity, and increase the viability and productivity of llama populations from other regions of Peru. This study is the first to report on the genetic diversity and population structure of the llama population in the Camelid Germplasm Bank—Quimsachata using microsatellite markers.
A representative sample of llama individuals (n = 251) was analyzed by using 13 microsatellites specific to the SAC. The results showed a high level of genetic diversity in the population of llamas, with an average of 12.08 alleles per microsatellite and an expected heterozygosity of 0.758. The Camelid Germplasm Bank—Quimsachata keeps the birth and mating records of each individual and that information is used to avoid mating between close relatives [13]. Thus, our results might indicate that this management strategy has contributed to maintaining high genetic variability, and it will be complemented with the routine assessment of the genetic variability by means of molecular markers such as microsatellites or single-nucleotide polymorphisms (SNPs). Furthermore, our results are similar to those from Argentinean (7.33–8.33 mean alleles per microsatellite and He of 0.47–0.9; [1,52]) and Bolivian (12.04 alleles per microsatellite and He = 0.68) llama populations [2]. Likewise, in those studies, the He was always higher than the Ho, and a higher value of He was reported in the Quimsachata llama population compared to the Bolivian llama population. Nonetheless, it is important to mention that they used different sets of microsatellite markers (many of which were individually used in this study), and in the case of Barreta et al. (2013), they were not necessarily specific for SACs. Thus, to make a more adequate comparison, we would ultimately have to use the same set of microsatellites used in those studies.
All the microsatellites analyzed were highly polymorphic (PIC > 0.5), especially YWLL08 (0.883), YWLL44 (0.863), YWLL59 (0.819), and LCA85 (0.8). Interestingly, we found alleles previously unrecorded in their original studies; for instance, YWLL08, YWLL44, YWLL59, LCA82, LCA83, GLM4, and LGU76 reported 13, 11, 10, 5, 7, 9, and 8 alleles in the first studies, respectively [22,26,27,28,40], whereas we recorded 19, 17, 19, 9, 8, 10, and 15, respectively (Table 1). However, this difference could be due to the distinctive population sample and size used by the initial investigators.
The FIS media value for the total population (0.063) and the microsatellites showing deviations from the HWE were explained by heterozygote deficiency. The excess of homozygotes in a domestic population indicates loss of genetic variability and could be explained by a lack of random mating, which occurs during the artificial selection of herds [1], population subdivisions (Wahlund effect), gene flow, or the existence of null alleles [53]. Indeed, we observed the presence of null alleles in four microsatellite markers (LCA82A, YWLL59A, LCA85A, and GLM4; Table S2), which could have potentially increased the observed homozygosity value. Additionally, other possible explanations for homozygote excess could be due to the evolutionary history of llamas, consisting of polygynous behavior in which herds contain an α male who controls the access of other males to its territory and expels young males while retaining the female individuals. This behavior of a young male’s exclusion is still present in managed populations [1,54]. A departure from the HWE could be due not only to a genotyping error (null alleles) in some of the microsatellites but also to the repetitive mating of individuals within the same herd and to the low quantity of breeding males.
Regarding the llama population structure from the Camelid Germplasm Bank—Quimsachata, although we found a large number of private alleles in the Ch’aku and Q’ara subpopulations (15 and 20, respectively; Table 3), the AMOVA analysis showed that most of the variation came from the variability within the Ch’aku and Q’ara subpopulations (93.95%), whereas only 1.02% was due to the variability between the subpopulations (Table S3). Furthermore, their genetic differentiation index (FST) was 0.0102 along with a high gene flow value (Nm = 30.5). These values, according to [47], are indicators of low genetic differentiation between the two subpopulations and a weak genetic structure. These results are confirmed by the cluster analysis and the correspondence factorial analysis (Figure 3) and are similar to the values estimated by [2], among and within regional groups of Bolivian llamas.
Importantly, the breeding management of the Bank has been highly controlled since its beginnings, that is, all Ch’aku and Q’ara males are reared together but kept separately from Ch’aku and Q’ara females (who are also reared together) and the males and females of the same phenotype are only brought together for breeding. Therefore, the low genetic structure is most likely a result of their own ancestral domestication process, which involved frequent exchanges of reproductive males and, hence, crossbreeding between both types. For future studies to elucidate a potential genetic differentiation between these phenotypes, we propose the use of genes related to observed phenotypic traits (e.g., the diameter of fiber), which could help to identify allelic variants related to the studied phenotypes [2], or the use of SNPs, which cover the whole llama genome [12].

5. Conclusions

This study is the first assessment of the genetic diversity and structure of the Peruvian llama (L. glama) population from the Camelid Germplasm Bank—Quimsachata utilizing microsatellites. The set of microsatellites used was highly polymorphic and, hence, can be utilized to track the genetic variability of these animals to avoid a reduction in the effective population size. Overall, our results show that the llama population of the Camelid Germplasm Bank—Quimsachata presents a high level of genetic diversity; thus, it can be considered as an adequate stock for the conservation of this natural resource. However, it would be important to identify other Peruvian llama populations that could help to increase the diversity of the Bank and enable a higher representation of Peruvian llamas.
On the other hand, although bred separately and managed as two distinct subpopulations, we observed low genetic differentiation between the Ch’aku and Q’ara phenotypes. Therefore, additional analysis with genes related to phenotypic observed differences should be carried out to assess potential genetic differentiation between these two phenotypes.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/5/541/s1, Table S1: Primer sets used for genotyping the microsatellites, Table S2: Existence and Frequencies of null alleles found in each microsatellite marker analyzed with the software Micro Checker 2.2.3 and CERVUS 3.0.3, respectively Table S3: Analysis of molecular variance (AMOVA) between the subpopulations of llamas Ch’aku and Q’ara.

Author Contributions

Conceptualization, E.A.V.-R.; investigation, G.F.P. and C.E.Y.-M.; methodology, E.A.V.-R.; visualization, G.F.P., C.E.Y.-M., G.A.G., and E.A.V.-R.; writing—original draft, G.F.P., C.E.Y.-M., G.A.G., and E.A.V.-R. All authors have read and agreed to the published version of the manuscript.

Funding

The Division of Genetic Resources and Biotechnology for its project on the use of molecular tools to characterize zoo genetic resources received the funding from the National Institute of Agricultural Innovation (INIA). The National Program of Agricultural Innovation (PNIA) funded this publication.

Acknowledgments

We would like to express our sincere gratitude to the Camelid Germplasm Bank—Quimsachata, especially to Teodosio Huanca and his team for helping to collect the samples; Gerardo Mamani and Evelyn Farfán for their insightful comments on the manuscript; and Dina L. Gutierrez and Pamela Burger for reading and commenting on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bustamante, A.V.; Maté, M.L.; Zambelli, A.; Vidal-Rioja, L.; Lamas, H.E.; Giovambattista, G. Análisis de diversidad genética en tres poblaciones de llamas (Lama glama) del noroeste argentino. Rev. Chil. Hist. Nat. 2006, 79, 175–184. [Google Scholar] [CrossRef]
  2. Barreta, J.; Gutiérrez-Gil, B.; Iñiguez, V.; Romero, F.; Saavedra, V.; Chiri, R.; Rodríguez, T.; Arranz, J.J. Analysis of genetic diversity in Bolivian llama populations using microsatellites. J. Anim. Breed. Genet. 2013, 130, 321–331. [Google Scholar] [CrossRef] [PubMed]
  3. Fernández, B.S. Situación Actual de los Camélidos Sudamericanos en Perú. Proyecto de Cooperación Técnica en Apoyo de la Crianza y Aprovechamiento de los Camélidos Sudamericanos en la Región Andina Tcp/Rla/2914; Organización de las Naciones Unidas para la Agricultura y la Alimentación (FAO): Roma, Italy, 2005. [Google Scholar]
  4. Llacsa, J. Determinación de Indicadores Biométricos para la Selección de Llamas (Lama glama) Productoras de Carne. Ph.D. Thesis, Tesis de Médico Veterinario y Zootecnista, Universidad Nacional del Altiplano, Puno, Peru, 2006. [Google Scholar]
  5. Wilson, D.E.; Reeder, D.M. Mammal Species of the World: A Taxonomic and Geographic Reference, 3rd ed.; Johns Hopkins University Press: Baltimore, MD, USA, 2005; Volume 1, p. 645, ISBN 0801882214, 9780801882210. [Google Scholar]
  6. Wheeler, J.C. Evolution and present situation of the South American Camelidae. Biol. J. Linn. Soc. 1995, 54, 271–295. [Google Scholar] [CrossRef]
  7. Iñiguez, L.C.; Alem, R.; Wauer, A.; Mueller, J. Fleece types, fiber characteristics and production system of an outstanding llama population from Southern Bolivia. Small Rumin. Res. 1998, 30, 57–65. [Google Scholar] [CrossRef]
  8. Stemmer, A.; Wurzinger, M.; Soelkner, J.; Valle Zárate, A.; Nuernberg, M.; Delgado, J. La llama de Ayopaya: Descripción de un recurso genético autóctono. Arch. Zootec. 2005, 54, 253–259. [Google Scholar]
  9. Maquera, F.E. Caracterización y Persistencia Fenotípica en Llamas Karas y Lanudas del Centro Experimental La Raya-Puno. Master’s Thesis, Universidad Nacional Agraria La Molina, Lima, Peru, 1991. [Google Scholar]
  10. INEI (Instituto Nacional de Estadística e Informática). IV Censo Nacional Agropecuario 2012, Lima, Peru. Available online: http://proyectos.inei.gob.pe/web/DocumentosPublicos/ResultadosFinalesIVCENAGRO.pdf (accessed on 25 July 2019).
  11. INEI (Instituto Nacional de Estadística e Informática). III Censo Nacional Agropecuario 1994, Lima, Peru. Available online: http://censos.inei.gob.pe/bcoCuadros/IIIcenagro.htm (accessed on 10 July 2019).
  12. Gutierrez, G.A.; Gutierrez, J.P.; Huanca, T.; Wurzinger, M. Challenges and opportunities of genetic improvement in alpacas and llamas in Peru. In Proceedings of the World Congress on Genetics Applied to Livestock Production, Auckland, New Zealand, 12–16 February 2018; p. 762. [Google Scholar]
  13. Huanca, T.; Apaza, N.; Gonzales, M. Experiencia del INIA en el fortalecimiento del Banco de germoplasma de camélidos domésticos. Arch. Latinoam. Prod. Anim. 2007, 15, 186–194. [Google Scholar]
  14. Wurzinger, M.; Gutiérrez, G. Analysis of a multi-stakeholder process during the start-up phase of two community-based llama breeding programs in Peru. Livest. Res. Rural. Dev. 2017, 29, 10. [Google Scholar]
  15. Hammond, K. Conservation of domestic animal diversity: Global overview. In Proceedings of the 5th World Congress on Genetics Applied to Livestock Production, Guelph, Canada, 7–12 August 1994; Volume 21, pp. 423–439. [Google Scholar]
  16. Luikart, G.; Gielly, L.; Excoffier, L.; Vigne, J.D.; Bouvet, J.; Taberlet, P. Multiple maternal origins and weak phylogeographic structure in domestic goats. Proc. Natl. Acad. Sci. USA 2001, 98, 5927–5932. [Google Scholar] [CrossRef] [Green Version]
  17. Falconer, D.S.; Mackay, T.F.C. Introduction to Quantitative Genetics, 4th ed.; Logmans Green: Harlow, UK, 1996; ISBN-13 978-0582243026. [Google Scholar]
  18. Ihara, N.; Takasuga, A.; Mizoshita, K.; Takeda, H.; Sugimoto, M.; Mizoguchi, Y.; Hirano, T.; Itoh, T.; Watanabe, T.; Reed, K.M.; et al. A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Res. 2004, 14, 1987–1998. [Google Scholar] [CrossRef] [Green Version]
  19. Slatkin, M. A measure of population subdivision based on microsatellite allele frequencies. Genetics 1995, 139, 457–462. [Google Scholar]
  20. Maté, M.L.; Bustamante, A.; Giovambattista, G.; De Lamo, D.; Von Thüngen, J.; Zambelli, A.; Vidal-Rioja, L. Genetic diversity and differentiation of guanaco populations from Argentina inferred from microsatellite data. Anim. Genet. 2005, 36, 316–321. [Google Scholar] [CrossRef] [PubMed]
  21. Paredes, M.; Machaca, J.; Azor, P.J.; Alonso-Moraga, A.; Membrillo, A.; Muñoz-Serrano, A. Genetic Differentiation of Six Peruvian Alpaca Populations. In Fibre Production in South American Camelids and Other Fibre Animals; Wageningen Academic Publishers: Wageningen, The Netherlands, 2011; pp. 161–166. [Google Scholar]
  22. Lang, K.D.M.; Wang, Y.; Plante, Y. Fifteen polymorphic dinucleotide microsatellites in llamas and alpacas. Anim. Genet. 1996, 27, 293. [Google Scholar] [CrossRef]
  23. Penedo, M.C.T.; Caetano, A.R.; Cordova, K.I. Microsatellite markers for South American camelids. Anim. Genet. 1998, 29, 411–412. [Google Scholar] [PubMed]
  24. Obreque, V.; Coogle, L.; Henney, P.J.; Bailey, E.; Mancilla, R.; García-Huidobro, J.; Hinrichsen, P.; Cothran, E.G. Characterization of 10 polymorphic alpaca dinucleotide microsatellites. Anim. Genet. 1998, 29, 461–462. [Google Scholar] [PubMed]
  25. Obreque, V.; Mancilla, R.; García-Huidobro, J.; Cothran, E.G.; Hinrichsen, P. Thirteen new dinucleotide microsatellites in Alpaca. Anim. Genet. 1999, 30, 397–398. [Google Scholar] [CrossRef] [PubMed]
  26. Penedo, M.C.T.; Caetano, A.R.; Cordova, K. Eight microsatellite markers for South American camelids. Anim. Genet. 1999, 30, 166–167. [Google Scholar] [CrossRef]
  27. Penedo, M.C.T.; Caetano, A.R.; Cordova, K.I. Six microsatellite markers for South American camelids. Anim. Genet. 1999, 30, 399. [Google Scholar] [CrossRef]
  28. Sarno, R.J.; David, V.A.; Franklin, W.L.; O’Brien, S.J.; Johnson, W.E. Development of microsatellite markers in the guanaco, Lama guanicoe: Utility for South American camelids. Mol. Ecol. 2000, 9, 1922–1924. [Google Scholar] [CrossRef]
  29. Wheeler, J.C.; Chikhi, L.; Bruford, M.W. Genetic analysis of the origins of domestic South American camelids. In Archaeology and Animal Domestication: New Genetic and Archaeological Paradigms: New Genetic and Archaeological Paradigms; Zeder, M., Ed.; University of California Press: Berkeley, CA, USA, 2006; pp. 329–341. [Google Scholar]
  30. Rodríguez, J.B.; Wheeler, J.C.; Dodd, C.S.; Bruford, M.W.; Rosadio, R.A. Determinación de parentesco en alpacas (Vicugna pacos) por medio del análisis de ADN microsatélite. Rev. Investig. Vet. del Peru 2004, 15, 113–119. [Google Scholar] [CrossRef] [Green Version]
  31. Marín, J.C.; Casey, C.S.; Kadwell, M.; Yaya, K.; Hoces, D.; Olazabal, J.; Rosadio, R.; Rodriguez, J.; Spotorno, A.; Bruford, M.W.; et al. Mitochondrial phylogeography and demographic history of the Vicuña: Implications for conservation. Heredity 2007, 99, 70–80. [Google Scholar] [CrossRef] [Green Version]
  32. Melo, C.; Manunza, A.; Melo, M.; Olivera, L.; Amills, M. Analysis of the mitochondrial diversity of alpacas in eight farming areas of the south of Peru. In Fibre Production in South American Camelids and Other Fibre Animals; Wageningen Academic Publishers: Wageningen, The Netherlands, 2011; pp. 87–91. [Google Scholar] [CrossRef]
  33. Kadwell, M.; Fernandez, M.; Stanley, H.F.; Baldi, R.; Wheeler, J.C.; Rosadio, R.; Bruford, M.W. Genetic analysis reveals the wild ancestors of the llama and the alpaca. Proc. R. Soc. B Biol. Sci. 2001, 268, 2575–2584. [Google Scholar] [CrossRef] [Green Version]
  34. Cano, Z.L.; Rosadio, A.R.; Maturrano, H.L.; Dávalos, R.R.; Wheeler, J.C. Caracterización fenotípica y análisis de ADN mitocondrial de llamas de marcapomacocha, Peru. Rev. Investig. Vet. Peru. 2012, 23, 388–398. [Google Scholar] [CrossRef] [Green Version]
  35. Heath, D.; Smith, P.; Williams, D.; Harris, P.; Arias-Stella, J.; Krüger, H. The heart and pulmonary vasculature of the llama (Lama glama). Thorax 1974, 29, 463–471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Cristofanelli, S.; Antonini, M.; Torres, D.; Polidori, P.; Renieri, C. Meat and carcass quality from Peruvian llama (Lama glama) and alpaca (Lama pacos). Meat Sci. 2004, 66, 589–593. [Google Scholar] [CrossRef]
  37. Cristofanelli, S.; Antonini, M.; Torres, D.; Polidori, P.; Renieri, C. Carcass characteristics of peruvian llama (Lama glama) and alpaca (Lama pacos) reared in the Andean highlands. Small Rumin. Res. 2005, 58, 219–222. [Google Scholar] [CrossRef]
  38. Quispe, E.C.; Rodríguez, T.C.; Iñiguez, L.R.; Mueller, J.P. Producción de fibra de alpaca, llama, vicuña y guanaco en Sudamérica. Anim. Genet. Resour. 2009, 45, 1–14. [Google Scholar] [CrossRef]
  39. Sambrook, J.; Russell, W. Molecular Cloning: A Laboratory Manual, 3rd ed.; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 2001. [Google Scholar]
  40. Bustamante, A.V.; Maté, M.L.; Zambelli, A.; Vidal-Rioja, L. Isolation and characterization of 10 polymorphic dinucleotide microsatellite markers for llama and guanaco. Mol. Ecol. Notes 2003, 3, 68–69. [Google Scholar] [CrossRef]
  41. De Arruda, M.P.; Gonçalves, E.C.; Schneider, M.P.C.; Da Costa Da Silva, A.L.; Morielle-Versute, E. An alternative genotyping method using dye-labeled universal primer to reduce unspecific amplifications. Mol. Biol. Rep. 2010, 37, 2031–2036. [Google Scholar] [CrossRef]
  42. Marshall, T.C.; Slate, J.; Kruuk, L.E.B.; Pemberton, J.M. Statistical confidence for likelihood-based paternity inference in natural populations. Mol. Ecol. 1998, 7, 639–655. [Google Scholar] [CrossRef] [Green Version]
  43. Belkhir, K.; Borsa, P.; Chikhi, L.; Raufaste, N.; Bonhomme, F. GENETIX 4.05, Logiciel sous WindowsTM pour la génétique des populations. In Laboratoire Génome, Populations, Interactions, CNRS UMR 5171; Université de Montpellier II: Montpellier, France, 2004. [Google Scholar]
  44. Raymond, M.; Rousset, F. GENEPOP (Version 1.2): Population Genetics Software for Exact Tests and Ecumenicism. J. Hered 1995, 86, 248–249. [Google Scholar] [CrossRef]
  45. Van Oosterhout, C.; Hutchinson, W.F.; Wills, D.P.; Shipley, P. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 2004, 4, 535–538. [Google Scholar] [CrossRef]
  46. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinforma 2005, 1. [Google Scholar] [CrossRef] [Green Version]
  47. Wright, S. Evolution and the Genetics of Populations, Volume 4, Variability within and among Natural Populations; University of Chicago Press: Chicago, IL, USA, 1984; p. 590, ISBN-13: 978-0226910529. [Google Scholar]
  48. FSTAT A Program to Estimate and Test Gene Diversities and Fixation Indices (Version 2.9.3). 2001. 149, 507–526. Available online: https://www2.unil.ch/popgen/softwares/fstat.htm (accessed on 20 February 2019).
  49. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar]
  50. Moritz, C. Conservation units and translocations: Strategies for conserving evolutionary processes. Hereditas 1999, 130, 217–228. [Google Scholar] [CrossRef]
  51. Aranguren-Méndez, J.; Jordana, J.; Gomez, M. Genetic diversity in Spanish donkey breeds using microsatellite DNA markers. Genet. Sel. Evol. 2001, 33, 433–442. [Google Scholar] [CrossRef] [PubMed]
  52. Bustamante, A.V.; Zambelli, A.; De Lamo, D.A.; Von Thungen, J.; Vidal-Rioja, L. Genetic variability of guanaco and llama populations in Argentina. Small Rumin. Res. 2002, 44, 97–101. [Google Scholar] [CrossRef]
  53. Callen, D.F.; Thompson, A.D.; Shen, Y.; Phillips, H.A.; Richards, R.I.; Mulley, J.C.; Sutherland, G.R. Incidence and origin of “null” alleles in the (AC)n microsatellite markers. Am. J. Hum. Genet. 1993, 52, 922. [Google Scholar] [CrossRef]
  54. Franklin, W.L. Biology, ecology, and relationship to man of the South American camelids. Mamm. Biol. S. Am. 1982, 6, 457–489. [Google Scholar]
Figure 1. Geographical localization of the Camelid Germplasm Bank—Quimsachata in the Department of Puno, Peru at around 4200 m above sea level (map created with DIVA-GIS software).
Figure 1. Geographical localization of the Camelid Germplasm Bank—Quimsachata in the Department of Puno, Peru at around 4200 m above sea level (map created with DIVA-GIS software).
Genes 11 00541 g001
Figure 2. Structure analysis of the global llama population by cluster assignment using the Bayesian analysis (K = 2). The colored bars indicate the probability of assignment to either Cluster 1 (red) or 2 (green). The black line divides the Ch’aku and Q’ara subpopulations.
Figure 2. Structure analysis of the global llama population by cluster assignment using the Bayesian analysis (K = 2). The colored bars indicate the probability of assignment to either Cluster 1 (red) or 2 (green). The black line divides the Ch’aku and Q’ara subpopulations.
Genes 11 00541 g002
Figure 3. Factorial correspondence analysis between the Ch’aku and Q’ara subpopulations. Yellow, Ch’aku; Blue, Q’ara.
Figure 3. Factorial correspondence analysis between the Ch’aku and Q’ara subpopulations. Yellow, Ch’aku; Blue, Q’ara.
Genes 11 00541 g003
Table 1. Genetic diversity parameters of each llama subpopulation and the global population.
Table 1. Genetic diversity parameters of each llama subpopulation and the global population.
Micro SatelliteCh’akuQ’araGlobal Population
AHoHePICAHoHePICAHoHePIC
GLM490.660.780.75100.760.780.66100.720.780.71
LAB1130.720.790.53130.700.730.5140.710.760.52
LCA5440.60.610.7580.530.560.7580.560.590.75
LCA6580.610.650.8280.610.640.7790.610.640.79
LCA7780.660.640.86100.670.680.86110.660.670.86
LCA8290.680.790.8270.520.720.8190.570.750.82
LCA8380.790.780.7770.720.740.7180.750.760.74
LCA8580.760.850.9070.740.810.8680.750.820.88
LGU76140.830.840.83140.750.790.78150.780.810.8
YWLL08170.900.910.59170.830.870.63190.860.890.62
YWLL44160.840.890.58170.850.870.59170.850.8750.59
YWLL59150.730.840.74170.750.840.70190.740.840.73
VOLP0390.740.680.6280.590.640.57100.640.660.59
Average10.50.730.770.7410.90.690.740.7112.080.710.760.72
A, allele number; Ho, observed heterozygosity; He, expected heterozygosity; PIC, polymorphic information content.
Table 2. Analysis of the Hardy–Weinberg equilibrium for each microsatellite of the global llama population of the Camelid Germplasm Bank—Quimsachata.
Table 2. Analysis of the Hardy–Weinberg equilibrium for each microsatellite of the global llama population of the Camelid Germplasm Bank—Quimsachata.
Global PopulationCh’aku PopulationQ’ara Population
MicrosatelliteDeficit of Heterozygotes (p-Value)Standard DeviationFISMicrosatelliteDeficit of Heterozygotes (p-Value)Standard DeviationFISMicrosatelliteDeficit of Heterozygotes (p-Value)Standard DeviationFIS
GLM4 0.00620.00140.073GLM4 0.00160.00040.145GLM4 0.28510.0120.0335
LAB1 0.00020.00020.056LAB1 000.087LAB1 0.1760.01280.0383
LCA54 0.02170.00510.04LCA54 0.083100.0211LCA54 0.04750.00540.0512
LCA65 0.0540.01340.054LCA65 0.24890.01670.0676LCA65 0.18580.01270.0467
LCA77 0.2590.0150.004LCA77 0.63920.0144−0.0216LCA77 0.21130.02060.0171
LCA82 0.00010.00020.223LCA82 0.01420.00180.1327LCA82 0.00240.00060.2788
LCA83 0.04630.00660.01LCA83 0.03020.0025−0.0152LCA83 0.29070.00890.0251
LCA85 0.00020.00010.089LCA85 0.00380.00060.1191LCA85 0.01550.00130.0788
LGU76 0.1390.02270.039LGU76 0.19390.01280.0149LGU76 0.26820.01550.0537
YWLL08 0.1820.04580.034YWLL08 0.47790.02050.0127YWLL08 0.11050.01130.0473
YWLL44 0.1180.01330.033YWLL44 0.17220.01330.0372YWLL44 0.22470.01530.0335
YWLL59 0.00010.00010.118YWLL59 000.1255YWLL59 0.03730.00680.109
VOLP03 0.4340.05110.027VOLP03 0.94120.0091−0.0785VOLP03 0.15080.01130.0905
FIS, coefficient of inbreeding.
Table 3. Private allele size of the 13 analyzed microsatellite markers of the Ch’aku and Q’ara llama subpopulations.
Table 3. Private allele size of the 13 analyzed microsatellite markers of the Ch’aku and Q’ara llama subpopulations.
MicrosatelliteCh’akuQ’ara
GLM4*214
LAB1177203
LCA54*159, 171, 173, 179
LCA65202200
LCA77264256, 262, 278
LCA82126, 132*
LCA83216*
LCA85214*
LGU76285279
YWLL08181, 201157, 189
YWLL44127133, 143
YWLL59102, 114138, 142, 148, 150
VOLP03153, 171191
Total1520
* No private alleles found.
Table 4. FST values of each microsatellite marker in the global population of llamas from the Camelid Germplasm Bank—Quimsachata.
Table 4. FST values of each microsatellite marker in the global population of llamas from the Camelid Germplasm Bank—Quimsachata.
MicrosatellitesFST
GLM40.005
LAB10.009
LCA540.033
LCA650.001
LCA770.005
LCA820.017
LCA830.022
LCA850.009
LGU760.009
YWLL080.016
YWLL440.005
YWLL590.001
VOLP030.002
Average0.01

Share and Cite

MDPI and ACS Style

Paredes, G.F.; Yalta-Macedo, C.E.; Gutierrez, G.A.; Veli-Rivera, E.A. Genetic Diversity and Population Structure of Llamas (Lama glama) from the Camelid Germplasm Bank—Quimsachata. Genes 2020, 11, 541. https://doi.org/10.3390/genes11050541

AMA Style

Paredes GF, Yalta-Macedo CE, Gutierrez GA, Veli-Rivera EA. Genetic Diversity and Population Structure of Llamas (Lama glama) from the Camelid Germplasm Bank—Quimsachata. Genes. 2020; 11(5):541. https://doi.org/10.3390/genes11050541

Chicago/Turabian Style

Paredes, Gabriela F., Claudia E. Yalta-Macedo, Gustavo A. Gutierrez, and Eudosio A. Veli-Rivera. 2020. "Genetic Diversity and Population Structure of Llamas (Lama glama) from the Camelid Germplasm Bank—Quimsachata" Genes 11, no. 5: 541. https://doi.org/10.3390/genes11050541

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop