Advances in Anthropology
Vol.06 No.01(2016), Article ID:63319,10 pages
10.4236/aa.2016.61001

Kinship and Phenotypic Divergence in the Ancient Population of the Puna Plateau of Northwestern Argentina

María L. Fuchs*, Héctor H. Varela, José A. Cocilovo

Departamento de Ciencias Naturales, Facultad de Ciencias Exactas, Físico-Química y Naturales, Universidad Nacional de Río Cuarto, Río Cuarto, Argentina

Copyright © 2016 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 1 December 2015; accepted 1 February 2016; published 4 February 2016

ABSTRACT

The Puna de Jujuy is a characteristic region of Northwest Argentina, and archaeological evidence shows remarkable development in this area during the period of Regional Development (PRD) or Late period. The aim of this paper is to study the genetic structure of the late population of the Puna de Jujuy (Argentina) using cranial metric traits. The sample comprised 302 adult individuals belonging to six different localities: Agua Caliente, Casabindo, Doncellas, Queta, Río Negro and Sorcuyo. The research included 27 cranial metric variables, a discriminant analysis and the calculation of Mahalanobis’ D2 distance and the fixation or genetic divergence index. The results showed little phenotypic differentiation between six locations; 1% of the genetic variation explaind the differences between groups, and the rest represented variations within the group. The sites of Queta and Sorcuyo had lower observed variance values than expected in the model of balance between genetic drift and gene flow, and Agua Caliente showed higher observed variance values than expected. The remaining sites showed values close to the ones expected in the model. The calculated distances suggested that the differentiation between localities was moderate. These results were compared to FST calculated for other local populations. A differential distribution of morphological variability between sites was observed. Agua Caliente, Casabindo, Doncellas and Río Negro showed higher morphological variability due to higher gene flow values, whereas Queta and Sorcuyo showed lower morphological variability.

Keywords:

Biological Affinity, Late period, Genetic Drift, Gene Flow

1. Introduction

The Puna de Jujuy is a characteristic region of Northwest Argentina; it is a high plateau with altitudes varying between 3500 and 4500 masl, cold and dry climate and wide temperature variations during day and night (Paoli et al., 2002) . The archaeological information available shows the presence of hunter-gatherers since the Early and Middle Holocene, but the evidence is most apparent in the Late Holocene (Aschero, 2011; Yacobaccio, 2011) . The emergence of farmer-pastoral societies constitutes the basis for the development of farmer-potter societies, which reach their peak in the period of Regional Development or Late period (Albeck & Ruiz, 2003; Albeck, 2007) .

The structure of a population may be determined by historical, political, social, economic, biological and geographical factors, all of which interact in complex ways with the processes of genetic drift, natural selection and gene flow, determining the temporal and spatial differentiation of subpopulations. We can approach the analysis of this structure by applying direct and indirect models, also known as bound and free models (Howells, 1973; Relethford, 1980; Relethford & Lees, 1982; Relethford, 2007) . These types of studies have been performed using gene frequencies, surnames, pedigrees and quantitative traits (Morton et al., 1968; Relethford et al., 1981; Varela, 1997) .

Several authors have made contributions to the study and estimation of different parameters in order to understand the structure of the population from quantitative traits; among these authors: Morton et al., 1968, Morton et al., 1971, Morton, 1973, 1975; Relethford, 1980, 1988, 1991a, 1991b, 1994; Relethford et al., 1981, Relethford & Lees, 1982; Relethford & Crawford, 1995 .

Following the explanation above, the aim of this paper is to study the genetic structure of the late population of the Puna de Jujuy using cranial metric traits. This involves knowing the existing biological variation within and among different subpopulations and the evolutionary factors that explain this variation.

One of the first contributions to establish the microevolutionary processes in the South Central Andes region consists of studying the temporal variation of five prehistoric groups from Northern coastal of Chile through cranial metric traits and the calculation of Mahalanobis’ D2 distance (Rothhammer et al., 1982) . The authors use a model that predicts an exponential increase of kinship over time. The model is applied to both metric and non-metric traits in populations from the Valley of Azapa and Pisagua, Chile (Cocilovo, 1995; Cocilovo & Rothhammer, 1996a, 1996b, 1999; Cocilovo & Varela, 1998) . For the prehispanic population of San Pedro de Atacama, the authors study the population structure using free and bound models in order to estimate kinship based phenotypic and chronological distance (Varela, 1997; Varela & Cocilovo, 2000) .

The method suggested by Relethford & Blangero (1990) for quantitative traits proposes a balance between gene flow and genetic drift, and allows us to estimate Wright’s FST as an indicator of genetic divergence between groups. The key assumptions of this model are to consider characters to be selectively neutral, that there is a similar mutation rate for all populations, that migration is constant and that the number of migrants is similar between any pair of populations. The deviations in relation to the model may be accounted for by variations in the effective size or by differences in the gene flow rate (Relethford & Blangero, 1990; Relethford, 1994, 1996; Relethford & Harpending, 1994; Relethford et al., 1997) . The model is applied to prehispanic populations from San Pedro de Atacama (Varela & Cocilovo, 2009, 2011) and Arica (Cocilovo et al., 2001; Varela & Cocilovo, 2002; Varela et al., 2004; Varela et al., 2006) in Northern Chile, with the aim of studying the temporal and spatial genetic divergence between the different periods of cultural development.

For Northwest Argentina, a research was carried out in the Calchaquí Valley which proposes the existence of a biological unit in the Late period, as a result of a common biological history and a shared sociopolitical and cultural organization (Baffi & Cocilovo, 1989-1990) . In the Quebrada de Humahuaca, it was possible to prove the existence of a biological heterogeneity higher than expected, which reflects the intervention of different entities in the configuration of groups that inhabited the area (Bordach & Cocilovo, 1991; Cocilovo et al., 1999; Varela et al., 1999; Cocilovo et al., 2001) . Finally, in the eastern sector of the Puna de Jujuy, what is proposed is an emerging structuring of the late population and relations with the Quebrada de Humahuaca and San Pedro de Atacama (Mendonça et al., 1990-1991) .

2. Materials and Methods

The skulls sample of the Puna de Jujuy is stored at the Museum of La Plata (University of La Plata) and in the Ethnographic Museum (University of Buenos Aires). The material belongs to different archaeological sites: Doncellas, Agua Caliente, Casabindo, Río Negro, Sorcuyo and Queta (see Figure 1) and their radiocarbon dating places the population of the Puna within the Late period, 1000 - 1450 A.D (Fuchs & Varela, 2013; Fuchs, 2014) .

The sex of the individuals was determined by observing the morphological traits of the skull (Acsádi & Nemeskéri, 1970; Bass, 1981; Buisktra & Ubelaker, 1994; White & Folkens, 2005) . The contribution of ancient DNA allowed us to certify that the sex of the individuals was the same as that determined by morphological traits (Postillone et al., 2013) . The age determination was performed according to Milner et al. (2000) , Meindl & Lovejoy (1985) , Molnar (1971) , Ubelaker (1984) and White & Folkens (2005) . The artificial cranial deformation was determined according to Dembo & Imbelloni (1938) , and the following types were considered: erect tabular, oblique tabular, erect circular, oblique circular and non-deformed. The observation of the metric traits was performed according to Bass (1981) , Buikstra & Ubelaker (1994) and Comas (1966) . In this study, 27 cranial metric variables were used (Table 1).

The sample included 302 adults of both sexes, with different types of artificial deformation, and from six sites (Agua Caliente, Casabindo, Doncellas, Queta, Río Negro and Sorcuyo) (see Table 2). In order to study the structure of the local population, the data were corrected for sex and deformation, using the general linear model (Seber, 1984) , i.e., we use the residuals of the regression between the dependent variables (cranial metric characters) and independent categorical variables (sex and artificial deformation). Later, the discriminant analysis was performed and Mahalanobis’ D2 distances between groups were calculated. Finally, a dendrogram was made with these distances using the Ward’s method for the six localities.

Furthermore, the genetic divergence index was calculated (FST) applying the model proposed for quantitative traits by Relethford & Blangero (1990) . The model can be expressed as: E (VPi) = VPw (1 − rii)/(1 − FST), where E (VPi) is the expected average phenotypic variance of the population i, VPi is the pooled average within-group phenotypic variance, rii is the genetic distance of population i to the centroid and FST is the average genetic distance to the centroid. For this study, we used the Rmet 5.0 software to calculate the value of FST, assuming equal relative population size and heritabilities of 0.55 and 0.78. Heritability (h2) of 0.55 was used since it is the value most used in the specific literature and an h2 of 0.78 which represents the estimated value for the population of the Puna de Jujuy based on the calculation of repeatability (Fuchs et al., 2014) .

Figure 1. Sites in the Puna de Jujuy, Argentine.

Table 1. Craniometrical variables.

Table 2. Sample used to analyze the population structure.

References. TE: erect tabular, TO: oblique tabular, CO: oblique circular and ND: non-deformed.

3. Results

The upper hemimatrix of Table 3 displays Mahalanobis distances and the lower hemimatrix shows the probability values for the comparison between pairs of localities. Figure 2 shows the distribution of the groups based on the first two canonical variables which account for 67% of the variability between groups. The first canonical variable separates mainly Agua Caliente and Casabindo, while the second canonical variable separates Queta and Sorcuyo.

Figure 3 shows the results of cluster analysis where Agua Caliente, Doncellas and Río Negro make up a group, Casabindo is associated later, then Sorcuyo, and finally Queta.

The results of the fixation index for h2 = 0.55 shows FST biased of 0.047, FST unbiased of 0.025 and standard error of 0.005, whereas for h2 = 0.78, values were 0.034 (FST biased), 0.013 (FST unbiased) and standard error of 0.005.

Table 4 shows the results of applying the model of balance between genetic drift and migration proposed by Relethford & Blangero (1990) : unbiased genetic distance of population i to the centroid (r(ii)), standard error (se), observed variance (VObs) and expected variance (VEsp) and the difference between both variances for heritability of 0.55 and 0.78. The sites of Queta and Sorcuyo showed less observed variance than expected, Agua Caliente

Figure 2. Locality distribution according to the centroids of the first two can- onical variables.

Figure 3. Cluster diagram among sites of the Puna de Jujuy.

Table 3. Mahalanobis’ D2 distance and values of probability for the Puna de Jujuy.

Wilks’ Lambda: 0.49313 approx. F (135, 1336) = 1.5261, p < 0.0002

showed greater observed variance than expected, and the rest of the sites presented similar values between observed and expected. Figure 4 shows the connection between the phenotypic variance and the distance of each group to the population centroid.

4. Discussion and Conclusion

The calculated distances suggest that the differentiation between localities is moderate, 60% of them are not statistically significant (Table 3). FST calculated with values of h2 = 0.78 indicates that 1% of the total genetic

Figure 4. Graphic representation of the drift/migration model in the Puna de Jujuy with h2 = 0.55 (a = left) and h2 = 0.78 (b = right). References: r(ii): genetic distance of population i to the centroid. Black points: observed values; gray points: expected values. 1: Agua Caliente; 2: Casabindo; 3: Doncellas; 4: Queta; 5: Río Negro; 6: Sorcuyo.

Table 4. Results for the model of balance drift/migration.

References: unbiased genetic distance of population i to the centroid (r(ii)), standard error (se), observed variance (VObs) and expected variance (VEsp) and the difference between both variances (Resi) for heritability (h2) of 0.55 and 0.78.

variation explains the difference between localities, and the remaining 99% constitutes the existing variability within localities. These results are compared to FST calculated for other local populations, for example, the Quitor Cemetery in San Pedro de Atacama, where there is a study on the temporal genetic divergence between the Formative, Tiwanaku and Regional Development periods (Varela & Cocilovo, 2009) . For this population, results show FST values of 0.027 (with h2 = 0.8), 0.046 (with h2 = 0.55) and 0.018 (with h2 = 1). Estimated values for h2 = 0.8 and 1 in San Pedro de Atacama are very similar to the values of the Puna de Jujuy, where FST values are 0.025 (with h2 = 0.55) and 0.013 (with h2 = 0.78). Later, the phenotypic differentiation between the Early, Middle and Late periods of the ayllus of San Pedro de Atacama is analyzed. Results show a moderate process of genetic divergence between periods, with an average FST value of 0.0481. The FST for the Early, Middle and Late periods is 0.114, 0.005 and 0.026 respectively (Varela & Cocilovo, 2011) . This last FST value for the Late period is very similar to the value obtained in the population of the Puna de Jujuy.

In the prehistoric groups of the coast and the Azapa Valley (Northern Chile), evidence shows that gene flow is lower on the coast than in the valley, minimal estimated FST for the total population (coast and valley) being 0.02, 0.01 for the coastal populations, and 0.006 for valley populations (Varela & Cocilovo, 2002) . FST values for the total population and for the population of the coast of Azapa are similar to the values obtained for the Puna de Jujuy.

The results of applying the drift/migration model (Relethford & Blangero, 1990) show that the sites of Queta and Sorcuyo present observed variance values lower than expected for those localities; this can be due to a gene flow that is lower than the average. In Agua Caliente, a higher variance than expected is observed, which may reflect a higher migratory effect than the expected average. The rest of the sites show values close to the expected by the model of balance between gene flow and genetic drift.

The site of Agua Caliente (Alfaro, 1983) , because of its greater gene flow, will show less kinship within-group and minor isolation with other sites of sub-region. This site is located on the same line than Doncellas but on the opposite edge of the river, it is a conglomerate or semi-conglomerate undefended, with rectangular and quadrangular enclosures forming simple and compound units.

By contrast, the lower gene flow of Queta and Sorcuyo is related to increase within-group kinship and isolation with other groups in the sub-region. Queta has been cited by archaeologists who carry out research or visit the site, but the only available information is Boman (1908) and Madrazo & Ottonello (1966) , who define it as a conglomerate village without defense. In a later contribution, Alfaro (1988) walks around the site and can only find a partially-opened cist and three utilitarian urns with the remains of infants.

The archaeological site of Sorcuyo (Casanova, 1938) or Pueblo Viejo de Tucute (Albeck, 1999; Albeck & Ruiz, 2003) is a large semiconglomerate village located in two hills on both sides of the Tucute stream, and appears to be inhabited as from late 900 A.D. until early 1500 A.D. (Albeck & Zaburlín, 2008) . Some years ago, Albeck (2007) proposed the entry of Highland tradition groups, which settled in the Casabindo area. The estimated date of this event is previous to the year 1300 A.D., when the village of Pueblo Viejo appears to be settled permanently. These two groups may have caused greater isolation as a consequence of the decrease in genetic exchange when compared to the rest of the population, or of practicing a reproductive system with more frequent mating within-group than between groups.

The distances calculated between the six locations show a moderate phenotypic differentiation between them. The results of the model application drift/migration (Relethford & Blangero, 1990) , reveals a differential distribution of morphological variability between sites. Agua Caliente, Casabindo, Doncellas and Río Negro show higher morphological variability, product of a higher gene flow, whereas Queta and Sorcuyo show lower variability. This differential distribution of the morphological variation could be explained by the entry of groups possibly of altiplano origin, or other regions in the Late Period. The greater isolation of these two groups may be the result of a decrease in the genetic exchange regarding the rest of the population, or by practicing a reproductive system that promotes mating within the group.

Acknowledgements

We thank the institutions that financed this research: Consejo Nacional de Investigaciones Científicas, Secretaría de Ciencia y Técnica-Universidad Nacional de Río Cuarto and Agencia Nacional de Promoción Científica y Técnica. We also want to thank the following people: Myriam Tarragó, Claudia Aranda, Karina Zucala, Héctor Pucciarelli, Mariano Del Papa, Julián Valetti and Maia Upstein.

Cite this paper

María L.Fuchs,Héctor H.Varela,José A.Cocilovo, (2016) Kinship and Phenotypic Divergence in the Ancient Population of the Puna Plateau of Northwestern Argentina. Advances in Anthropology,06,1-10. doi: 10.4236/aa.2016.61001

References

  1. 1. Acsádi, G., & Nemeskéri, J. (1970). History of Human Life Span and Mortality. Budapest: Akadémi Kiadó.

  2. 2. Albeck, M. E. (1999). El recinto R1 de Pueblo Viejo de Tucute: Sorcuyo revisitado. Chungara Revista de Antropología Chilena, 130, 143-159.

  3. 3. Albeck, M. E. (2007). El Intermedio Tardío: Interacciones económicas y políticas en la Puna de Jujuy. In V. Williams, B. Ventura, A. Callegari, & H. Yacobaccio (Eds.), Sociedades Precolombinas Surandinas. Temporalidad, Interacción y Dinámica cultural del NOA en el ámbito de los Andes Centro-Sur (pp. 125-146). Buenos Aires: Taller Internacional de Arqueología del NOA y Andes Centro Sur.

  4. 4. Albeck, M. E., & Ruíz, M. (2003). El tardío en la Puna de Jujuy: Poblados, etnias y territorios. Cuadernos FHyCS-UNJu, 20, 199-221.

  5. 5. Albeck, M. E., & Zaburlín, M. A. (2008). Aportes a la cronología de los asentamientos agropastoriles de la Puna de Jujuy. Relaciones de la Sociedad Argentina de Antropología, 33, 155-180.

  6. 6. Alfaro, L. C. (1983). Investigación arqueológica en la cuenca del Río Doncellas (Provincia de Jujuy). Integración de la Puna Jujena a los centros cúlticos andinos. Relaciones de la Sociedad Argentina de Antropología, 15, 25-47.

  7. 7. Alfaro, L. C. (1988). Investigación Arqueológica en la Cuenca del río Doncellas. Dpto. Cochinoca, Jujuy. Reconstrucción de una cultura olvidada en la Puna Jujena. San Salvador de Jujuy: Talleres Gráficos del Boletín Oficial e Imprenta del Estado de la Provincia de Jujuy.

  8. 8. Aschero, C. A. (2011). Holoceno Medio en la Puna Argentina: Dos puntos de observación para la cultura material y una perspectiva regional. In M. Mondini, J. G. Martínez, H. Muscio, & M. B. Marconetto (Eds.), Poblaciones Humanas y Ambientes en el Noroeste Argentino Durante el Holoceno Medio (pp. 33-43). Córdoba: Talleres Gráficos de Corintios 13.

  9. 9. Baffi, E. I., & Cocilovo, J. A. (1989-1990). La población prehispánica tardía del sector septentrional del valle Calchaquí. Runa, 19, 11-26.

  10. 10. Bass, W. M. (1981). Human Osteology: A Laboratory and Field Manual of the Human Skeleton. Springfield: University of Missouri Press.

  11. 11. Boman, E. (1908). Antiquitiés de la région andine de la République Argentine et du désert d’Atacama (Vol. 1). Paris: Imprimerie nationale, H. Le Soudier.

  12. 12. Bordach, M. A., & Cocilovo, J. A. (1991). Composición y estructura de la población prehispánica de la quebrada de Humahuaca, primera aproximación. Antropología Biológica, 1, 15-31.

  13. 13. Buikstra, J. E., & Ubelaker, D. H. (1994). Standards for Data Collection from Human Skeletal Remains. Research Series No. 44, Fayetteville: Arkansas Archeological Survey.

  14. 14. Casanova, E. (1938). Investigaciones arqueológicas en Sorcuyo, Puna de Jujuy. Anales del Museo Argentino de Ciencias Naturales, 80, 423-462.

  15. 15. Cocilovo, J. A. (1995). La Población Prehistórica de Pisagua. Continuidad y Cambio Biocultural en el Norte de Chile. Tesis doctoral inédita, Córdoba: Facultad de Ciencias Exactas, Físicas y Naturales, Universidad Nacional de Córdoba.

  16. 16. Cocilovo, J. A., & Rothhammer, F. (1996a). El valle de Azapa en el norte de Chile. Análisis de un modelo biocronológico para establecer la estructura de la población. Revista Argentina de Antropología Biológica, 1, 232-245.

  17. 17. Cocilovo, J. A., & Rothhammer, F. (1996b). Methodological Approaches for the Solution of Ethnohistorical Problems: Bioassay of Kinship in Prehistoric Populations of Arica, Chile. Homo, 47, 117-190.

  18. 18. Cocilovo, J. A., & Rothhammer, F. (1999). Microevolución morfológica y extinción del parentesco en asentamientos humanos prehistóricos del Valle de Azapa, Chile. Revista Chilena de Historia Natural, 72, 213-218.

  19. 19. Cocilovo, J. A., & Varela, H. H. (1998). Variación morfológica, estructura canónica y bioensayo de parentesco en poblaciones prehistóricas del norte de Chile. Chungara Revista de Antropología Chilena, 30, 75-85.
    http://dx.doi.org/10.4067/s0717-73561998000100006

  20. 20. Cocilovo, J. A., Varela, H. H., & Valdano, S. G. (2001). Estructura de la población antigua de la Quebrada de Humahuaca. In E. E. Berberián, & A. E. Nielsen (Eds.), Historia Argentina Prehispánica (pp. 265-287). Córdoba: Editorial Brujas.

  21. 21. Cocilovo, J. A., Varela, H. H., Baffi, E. I., Valdano, S. G., & Torres, M. F. (1999). Estructura y composición de la población antigua de la Quebrada de Humahuaca. Análisis exploratorio. Relaciones de la Sociedad Argentina de Antropología, 24, 59-80.

  22. 22. Cocilovo, J. A., Varela, H. H., Espoueys, O., & Standen, V. G. (2001). El proceso microevoulutivo de la población nativa antigua de Arica. Chungara Revista de Antropología Chilena, 33, 13-20.

  23. 23. Comas, J. (1966). Manual de Antropología Física. México: Universidad Nacional Autónoma de México, Instituto de Investigaciones Históricas, Sección de Antropología.

  24. 24. Dembo, A., & Imbelloni, J. (1938). Deformaciones intencionales del cuerpo humano de carácter étnico. Buenos Aires: Humanior.

  25. 25. Fuchs, M. L. (2014). Caracterización bioestructural de la población antigua de la Puna de Jujuy (Argentina). Tesis doctoral inédita, Río Cuarto: Facultad de Ciencias Exactas, Físico-Químicas y Naturales, Universidad Nacional de Río Cuarto.

  26. 26. Fuchs, M. L., & Varela, H. H. (2013). Fechados radiocarbónicos de colecciones osteológicas de la Puna de Jujuy, Argentina. Relaciones de la Sociedad Argentina de Antropología, 38, 553-558.

  27. 27. Fuchs, M. L., Valetti, J. A., Varela, H. H., & Cocilovo, J. A. (2014). Varianza ambiental y estabilidad del desarrollo en la población prehispánica de la Puna de Jujuy. In V. Seldes, & M. S. Gheggi (Eds.), Antropología Biológica y estudios del comportamiento mortuorio de los pueblos prehispánicos del Noroeste Argentino (pp. 23-35). Buenos Aires: Fundación Azara.

  28. 28. Howells, W. W. (1973). Cranial Variation in Man: A Study by Multivariate Analysis of Patterns of Difference among Recent Human Population. Papers of the Peabody Museum of Archaeology and Ethnology, 67, 1-259.

  29. 29. Madrazo, G. B., & Ottonello de García Reynoso, M. M. (1966). Tipos de instalación prehispánica en la región de la Puna y su borde (No. 1). Olavarría: Museo Etnográfico Municipal “Dámaso Arce”.

  30. 30. Meindl, R. S., & Lovejoy, C. O. (1985). Ectocranial Suture Closure: A Revised Method for the Determination of Skeletal Age at Death Based on the Lateral-Anterior Sutures. American Journal of Physical Anthropology, 68, 57-66.
    http://dx.doi.org/10.1002/ajpa.1330680106

  31. 31. Mendonca, O. J., Cocilovo, J. A., & Valdano, S. (1990-1991). La población prehistórica de la cuenca de Miraflores-Guayatayoc-Salinas Grandes en el sector Oriental de la Puna Jujena. Avances en Arqueología, 1, 59-80.

  32. 32. Milner, G. R., Wood, J. W., & Boldsen, J. L. (2000). Paleodemography. In M. A. Katzenberg, & S. R. Saunders (Eds.), Biological Anthropology of Human Skeleton (pp. 467-497). New York: Willey-Liss.

  33. 33. Molnar, S. (1971). Human Tooth Wear. Tooth Function and Cultural Variability. American Journal of Physical Anthropology, 34, 175-189.
    http://dx.doi.org/10.1002/ajpa.1330340204

  34. 34. Morton, N. E. (1973). Kinship and Population Structure. In N. E. Morton (Ed.), Genetic Structure of Population (pp. 66-71). Honolulu: University of Hawaii Press.

  35. 35. Morton, N. E. (1975). Kinship, Information and Biological Distance. Theoretical Population Biology, 7, 246-255.
    http://dx.doi.org/10.1016/0040-5809(75)90019-2

  36. 36. Morton, N. E., Miki, C., & Yee, S. (1968). Bioassay of Population Structure under Isolation by Distance. American Journal of Human Genetics, 20, 411-419.

  37. 37. Morton, N. E., Yee, S., Harris, D. E., & Lew, R. (1971). Bioassay of Kinship. Theoretical Population Biology, 2, 507-524.
    http://dx.doi.org/10.1016/0040-5809(71)90038-4

  38. 38. Paoli, H., Bianchi, A. R., Yanez, C. E., Volante, J. N., Fernández, D. R., Mattalía, M. C., & Noé, Y. E. (2002). Recursos hídricos de la Puna, valles y bolsones áridos del Noroeste Argentino. Salta: Convenio INTA EEA Salta-CIED.

  39. 39. Postillone, M. B., Avena, S. A., Carnese, F. R., & Dejean, C. B. (2013). Análisis genético de una muestra precolombina del sitio Agua Caliente (Jujuy, Argentina). In J. R. Bárcena, & S. E. Martin (Eds.), Actas XVIII Congreso Nacional de Arqueología Argentina (p. 47). La Rioja: XVIII Congreso Nacional de Arqueología Argentina.

  40. 40. Relethford, J. H. (1980). Bioassay of Kinship from Continuous Traits. Human Biology, 52, 689-700.

  41. 41. Relethford, J. H. (1988). Estimation of Kinship and Genetic Distance from Surnames. Human Biology, 60, 475-492.

  42. 42. Relethford, J. H. (1991a). Effect of Population Size on Marital Migration Distance. Human Biology, 63, 95-98.

  43. 43. Relethford, J. H. (1991b). Genetic Drift and Anthropometric Variation Ireland. Human Biology, 63, 155-165.

  44. 44. Relethford, J. H. (1994). Craniometric Variation among Modern Human Populations. American Journal of Physical Anthropology, 95, 53-62.
    http://dx.doi.org/10.1002/ajpa.1330950105

  45. 45. Relethford, J. H. (1996). Genetic Drift Can Obscure Population History: Problem and Solution. Human Biology, 68, 29-44.

  46. 46. Relethford, J. H. (2007). The Use of Quantitative Traits in Anthropological Genetic Studies of Population Structure and History. In M. H. Crawford (Ed), Anthropological Genetics: Theory, Methods and Applications (pp. 187-209). Cambridge: Cambridge University Press.

  47. 47. Relethford, J. H., & Lees, F. C. (1982). The Use of Quantitative Traits in the Study of Human Population Structure. American Journal of Physical Anthropology, 25, 113-132.
    http://dx.doi.org/10.1002/ajpa.1330250507

  48. 48. Relethford, J. H., & Blangero, J. (1990). Detection of Differential Gene Flow from Patterns of Quantitative Variation. Human Biology, 62, 5-25.

  49. 49. Relethford, J. H., & Crawford, M. H. (1995). Anthropometric Variation and the Population History of Ireland. American Journal of Physical Anthropology, 96, 25-38.
    http://dx.doi.org/10.1002/ajpa.1330960104

  50. 50. Relethford, J. H., & Harpending, H. C. (1994). Craniometric Variation, Genetic Theory, and Modern Human Origins. American Journal of Physical Anthropology, 95, 249-270.
    http://dx.doi.org/10.1002/ajpa.1330950302

  51. 51. Relethford, J. H., Lees, F. C., & Crawford, M. H. (1981). Population Structure and Anthropometric Variation in Rural Western Ireland: Isolation by Distance and Analysis of the Residuals. American Journal of Physical Anthropology, 55, 233-245.
    http://dx.doi.org/10.1002/ajpa.1330550212

  52. 52. Relethford, J. H., Crawford, M. H., & Blangero, J. (1997). Genetic Drift and Gene Flow in Post-Famine Ireland. Human Biology, 69, 443-465.

  53. 53. Rothhammer, F., Cocilovo, J. A., Quevedo, S., & Llop, E. (1982). Microevolution in Prehistoric Andean Populations: I. Chronology Craniometric Variation. American Journal of Physical Anthropology, 58, 391-396.
    http://dx.doi.org/10.1002/ajpa.1330580406

  54. 54. Seber, G. A. F. (1984). Multivariate Observations. New York: John Wiley & Sons.
    http://dx.doi.org/10.1002/9780470316641

  55. 55. Ubelaker, D. (1984). Human Skeletal Remains. Washington DC: Taraxacum.

  56. 56. Varela, H. H. (1997). La población prehistórica de San Pedro de Atacama. Composición. Estructura y Relaciones Biológicas. Tesis doctoral inédita, Río Cuarto: Facultad de Ciencias Exactas, Físico-Químicas y Naturales, Universidad Nacional de Río Cuarto.

  57. 57. Varela, H. H., & Cocilovo, J. A. (2000). Structure of the Prehistoric Population of San Pedro de Atacama. Current Anthropology, 41, 125-132.
    http://dx.doi.org/10.1086/300112

  58. 58. Varela, H. H., & Cocilovo, J. A. (2002). Genetic Drift and Gene Flow in to Prehistoric Population of the Azapa Valley and Coast, Chile. American Journal of Physical Anthropology, 118, 259-267.
    http://dx.doi.org/10.1002/ajpa.10075

  59. 59. Varela, H. H., & Cocilovo, J. A. (2009). Microevolución en San Pedro de Atacama (Norte de Chile): El cementerio de Quitor. Latin American Antiquity, 20, 333-342.

  60. 60. Varela, H. H., & Cocilovo, J. A. (2011). Divergencia fenotípica en los oasis de san Pedro de Atacama, norte de Chile. Estudios Atacamenos, 42, 101-112.
    http://dx.doi.org/10.4067/S0718-10432011000200006

  61. 61. Varela, H. H., Cocilovo, J. A., Baffi, E. I., & Valdano, S. G. (1999). La población antigua de la Quebrada de Humahuaca y sus relaciones biológicas con áreas aledanas. Revista Espanola de Antropología Biológica, 20, 7-24.

  62. 62. Varela, H. H., Cocilovo, J. A., Quevedo, S., & Costa, M. A. (2004). La estructura de la población de pescadores del Periodo Arcaico Tardío y Formativo de Arica. Boletín del Museo Nacional de Historia Natural, 53, 149-160.

  63. 63. Varela, H. H., Cocilovo, J. A., Santoro, C. M., & Rothhammer, F. (2006). Microevolution of Human Archaic Groups of Arica, Northern Chile, and Its Genetic Contribution to Populations from the Formative Period. Revista Chilena de Historia Natural, 79, 185-193.
    http://dx.doi.org/10.4067/S0716-078X2006000200005

  64. 64. White, T. D., & Folkens, P. A. (2005). The Human Bone Manual. London: Elsevier.

  65. 65. Yacobaccio, H. D. (2011). Hacia una ecología humana para el Holoceno Medio. In M. Mondini, J. G. Martínez, H. Muscio, & M. B. Marconetto (Eds.), Poblaciones humanas y ambientes en el Noroeste Argentino durante el Holoceno Medio (pp. 23-31). Córdoba: Talleres Gráficos de Corintios 13.

NOTES

*Corresponding author.