The causes of craniofacial variation among human populations have been the subject of controversy. In this work, we studied aboriginal populations from southern South America, the last continental region peopled by humans and with a wide range of ecological conditions. Because of these characteristics, southern South America provides a unique opportunity to study the relative importance of random and nonrandom factors in human diversification. Previous craniometric studies recognized remarkable differences among populations from this region, usually resorting to random factors as the main explanation. In contrast, here we suggest, using tests based on quantitative genetic models, that: (1) the rate of craniofacial divergence among these populations is too high and (2) the patterns of variation within and between populations are too different to be explained by genetic drift alone. In addition, the among-sample craniofacial variation is correlated with climate and diet but not with mtDNA variation. We suggest that the influence of nonrandom factors (e.g., plasticity, selection) on human craniofacial diversification in regions with large ecological variation is more important than generally acknowledged and capable to generate large craniofacial divergence in a short period of time. These results bring nonrandom factors into focus for the interpretation of human craniofacial variation. © 2009 The Author(s).634978993Ackermann, R.R., Cheverud, J.M., Discerning evolutionary processes in patterns of tamarin (genus Saguinus) craniofacial variation (2002) Am. J. Phys. Anthropol., 117, pp. 260-271Ackermann, R.R., Cheverud, J.M., Detecting genetic drift versus selection in human evolution (2004) Proc. Natl. Acad. Sci. USA, 101, pp. 17946-17951Adams, D.C., Rohlf, F.J., Slice, D.E., Geometric morphometrics: Ten years of progress following the 'revolution.' (2004) Ital. J. Zool., 71, pp. 5-16Anderson, D.G., Gillam, J.C., Paleoindian interaction and mating networks: Reply to Moore and Moseley (2001) Am. Antiq., 66, pp. 530-535Barton, N.H., Briggs, D.E.G., Eisen, J.A., Goldstein, D.B., Patel, N.H., (2007) Evolution., , and. Cold Spring Harbor Laboratory Press, New YorkBeals, K.L., Smith, C.L., Dodd, S.M., Climate and the evolution of brachycephalization (1983) Am. J. Phys. Anthropol., 62, pp. 425-437Berberián, E.E., Nielsen, A.E., (2001) Historia Argentina Prehispánica., , and. Brujas, Córdoba, ArgentinaBernal, V., Procesos de diferenciación biológica entre poblaciones humanas del Holoceno tardío de Patagonia. Una aproximación desde la variación métrica dental (2008) Ph.D. Dissertation., , Universidad Nacional de La Plata, La Plata, ArgentinaBernal, V., Perez, S.I., Gonzalez, P.N., Variation and causal factors of craniofacial robusticity in Patagonian hunter-gatherers from the late Holocene (2006) Am. J. Hum. Biol., 18, pp. 748-765Binford, L.R., Constructing frames of reference (2000) An Analytical Method for Archaeological Theory Building Using Hunter-gatherer and Environmental Data Sets., , Univ. of California Press, Berkeley, CABookstein, F.L., Morphometric tools for landmark data (1991) Geometry and Biology., , Cambridge Univ. Press, CambridgeBookstein, F.L., Landmark methods for forms without landmarks: Localizing group differences in outline shape (1997) Med. Image Anal., 1, pp. 225-243Bookstein, F.L., Streissguth, A.P., Sampson, P.D., Connor, P.D., Barr, H.M., Corpus callosum shape and neuropsychological deficits in adult males with heavy fetal alcohol exposure (2002) Neuroimage, 15, pp. 233-251Borrero, L.A., The prehistoric exploration and colonization of Fuego-Patagonia (1999) J. World Prehistory, 13, pp. 321-355Brandt, M., Siegel, M.I., The effects of stress on cortical bone thickness in rodents (1978) Am. J. Phys. Anthropol., 49, pp. 31-34Buikstra, J., Ubelaker, D., (1994) Standards for Data Collection from Human Skeletal Remains, , and. Arkansas Archaeological Survey, Research Series 44, FayettevilleCadien, J., Harris, E., Jones, W., Mandarino, L., Biological lineages, skeletal populations and microevolution (1976) Yearb. Phys. Anthropol., 18, pp. 194-201Carroll, S.P., Hendry, A.P., Reznick, D.N., Fox, C.W., Evolution on ecological time-scales (2007) Funct. Ecol., 21, pp. 387-393Carson, E.A., Maximum likelihood estimation of human craniometric heritabilities (2006) Am. J. Phys. Anthropol., 131, pp. 169-180Cavalli-Sforza, L.L., Bodmer, W.F., (1971) The Genetics of Human Populations., , and. General Publishing Company, TorontoCavalli-Sforza, L.L., Menozzi, P., Piazza, A., (1994) The History and Geography of Human Genes., , and. Princeton Univ. Press, Princeton, NJCheverud, J.M., Phenotipic, genetic and environmental morphological integration in the cranium (1982) Evolution, 36, pp. 499-516Cheverud, J.M., A comparison of genetic and phenotypic correlations (1988) Evolution, 42, pp. 958-968Cheverud, J.M., Morphological integration in the saddle-back tamarin (Saguinus fuscicollis) cranium (1995) Am. Nat., 145, pp. 63-89Cheverud, J.M., The relationship between development and evolution through heritable variation (2007) Tinkering: The Microevolution of Development., pp. 55-65. , Pp. in. G. Bock. and. J. Goode, eds. John Wiley & Sons, Ltd., West Sussex, UKCurtis, N., Kupczik, K., O'Higgins, P., Moazen, M., Fagan, M., Predicting skull loading: Applying multibody dynamics analysis to a macaque skull (2008) Anat. Rec., 291, pp. 491-501Dejean, C.B., Lanata, J.L., Martino, L., Carnese, F.R., Osella, A., Demografía y distribución de haplogrupos mitocondriales durante la dispersión inicial en las Américas (2007) Revista Argentina de Antropología Biológica., 9, p. 51Fenner, J., Cross-cultural estimation of the human generation interval for use in genetics-based population divergence studies (2005) Am. J. Phys. Anthropol., 128, pp. 415-423Frisancho, A.R., (1996) Human Adaptation and Accommodation., , The Univ. of Michigan Press, Ann. Arbor, MIGarcía-Bour, J., Pérez-Pérez, A., Álvarez, S., Fernández, E., López-Parra, A.M., Arroyo-Pardo, E., Turbón, D., Early population differentiation in extinct aborigines from Tierra del Fuego-Patagonia: Ancient mtDNA sequences and Y-chromosome STR characterization (2004) Am. J. Phys. Anthropol., 123, pp. 361-370Gil, A., Shelnut, N., Neme, G., Tykot, R., Michieli, C., Isótopos estables y dieta humana en el centro oeste: Datos muestras de San Juan (2006) Revista Cazadores Y Recolectores del Cono Sur., 1, pp. 149-161Gower, J.C., Statistical methods of comparing different multivariate analyses of the same data (1971) Mathematics in the Archaeological and Historical Sciences., pp. 138-149. , Pp. in. F. R. Hodson, D. G. Kendall. and. P. Tautu, eds. Edinburgh Univ. Press, EdinburghHair, J.F., Tatham, R.L., Anderson, R.E., Black, W., (1998) Multivariate Data Analysis., , and. Prentice Hall, New YorkHallgrímsson, B.K., Willmore, C.D., Cooper, D.M.L., Craniofacial variability and modularity in macaques and mice (2004) J. Exp. Zool. (Mol. Dev. Evol.), 302, pp. 207-225Hallgrímsson, B., Lieberman, D.E., Liu, W., Hutchinson, A.F., Jirik, F.R., Epigenetic interactions and the structure of phenotypic variation in the cranium (2007) Evol. Dev., 9, pp. 76-91Hanken, J., Hall, B., (1993) The Skull., , and. Univ. of Chicago Press, ChicagoHarvati, K., Weaver, T.D., Human cranial anatomy and the differential preservation of population history and climate signatures (2006) Anat. Rec. Part a, 288, pp. 1225-1233Hernández, M., Lalueza Fox, C., García-Moro, C., Fueguian cranial morphology: The adaptation to a cold, harsh environment (1997) Am. J. Phys. Anthropol., 103, pp. 103-117Hendry, A.P., Kinnison, M.T., The pace of modern life: Measuring rates of contemporary microevolution (1999) Evolution, 53, pp. 1637-1653Holliday, T.W., Body proportions in Late Pleistocene Europe and modern human origins (1997) J. Hum. Evol., 32, pp. 423-447Hunt, G., Phenotypic variance inflation in fossil samples: An empirical assessment. (2004) Paleobiology, 30, pp. 487-506Hunt, G., Bell, M.A., Travis, M.P., Evolution toward a new adaptive optimum: Phenotypic evolution in a fossil stickleback lineage (2007) Evolution, 623, pp. 700-710Jackson, D.A., PROTEST: A Procrustean randomization test of community environment concordance (1995) Écoscience, 2, pp. 297-303Katzmarzyk, P.T., Leonard, W.R., Climatic influences on human body size and proportions: Ecological adaptations and secular trends (1998) Am. J. Phys. Anthropol., 106, pp. 483-503Kelly, J., Grammer, S., Tykot, R.H., Belardi, J., Borrero, L., Gomez-Otero, J., Guichón, R., Integrating archaeological, ethnographic and analytical subsistence data: A case study from Patagonia, South America. (2001) 66th Annual Meeting of the Society for American Archaeology., , New Orleans, LouisianaKennett, D.J., Winterhalder, B., (2006) Behavioral Ecology and the Transition to Agriculture., , and. Univ. of California Press, Berkeley, CAKlingenberg, C.P., Monteiro, L.R., Distances and directions in multidimensional shape spaces: Implications for morphometric applications (2005) Syst. Biol., 54, pp. 678-688Lahr, M.M., Patterns of modern human diversification: Implications for Amerindian origins (1995) Yearb. Phys. Anthropol., 38, pp. 163-198Lalueza, C., Pérez-Pérez, A., Prats, E., Cornudella, L., Turbón, D., Lack of founding Amerindian mitochondrial DNA lineages in extinct aborigines from Tierra del Fuego-Patagonia (1997) Hum. Mol. Genet., 6, pp. 41-46Lanata, J.L., Martino, L., Osella, A., Garcia-Herbst, A., Ambiente y demografía durante la dispersión humana inicial en Sudamérica. (2006) Ecología Histórica: Interacciones Sociedad-ambiente en Distintas Escalas Espacio Temporales., , and. In. C. López. and. G. Ospina. eds. Universidad Tecnológica de Pereira, Pereira, Colombia. In pressLande, R., Statistical tests for natural selection on quantitative characters (1977) Evolution, 31, pp. 442-444Lande, R., Quantitative genetic analysis of multivariate evolution, applied to brain: Body size allometry (1979) Evolution, 33, pp. 402-416Legendre, P.F.J., Lapointe, Casgrain, P., Modeling brain evolution from behavior: A permutational regressio approach (1994) Evolution, 48, pp. 1487-1499Leonard, W.R., Snodgrass, J.J., Sorensen, M.V., Metabolic adaptation in indigenous Siberian populations (2005) Annu. Rev. Anthropol., 34, pp. 451-471Lemos, B., Marroig, G., Cerqueira, R., Evolutionary rates and stabilizing selection in large-bodied opossum skulls (Didelphimorphia: Didelphidae) (2001) J. Zool. (Lond.), 255, pp. 181-189Levin, S.A., The problem of pattern and scale in ecology: The Robert H. MacArthur award lecture (1992) Ecology, 73, pp. 1943-1967Lieberman, D.E., Pearson, O.M., Mowbray, K.M., Basicranial influence on overall cranial shape (2000) J. Hum. Evol., 38, pp. 291-315Lockwood, C.A., Kimbel, W.H., Lynch, J.M., Morphometrics and hominoid phylogeny: Support for a chimpanzee-human clade and differentiation among great ape subspecies (2004) Proc. Natl. Acad. Sci. USA, 101, pp. 4356-4360Losos, J.B., Uncertainty in the reconstruction of ancestral character states and limitations on the use of phylogenetic comparative methods (1999) Anim. Behav., 58, pp. 1319-1324Losos, J.B., U'Arheit, K.I., Schoener, T.M., Adaptive differentiation following experimental island colonization in Anolis lizards (1997) Nature, 387, pp. 70-73Lynch, M., The rate of morphological evolution in mammals from the standpoint of the neutral expectation (1990) Am. Nat., 136, pp. 727-741Marriog, G., Cheverud, J., Cranial evolution in Sakis (Pithecia, Platyrrhini) I: Interspecific differentiation and allometric patterns (2004) Am. J. Phys. Anthropol., 125, pp. 266-278Mitteroecker, P., Bookstein, F., The evolutionary role of modularity and integration in the hominoid cranium (2008) Evolution, 62, pp. 943-958Monteiro, L.R., Gomes Jr., J.L., Morphological divergence rate tests for natural selection: Uncertainty of parameter estimation and robustness of results (2005) Genet. Mol. Biol., 28, pp. 345-355Moraga, M., Rocco, P., Miquel, J.F., Nervi, F., Llop, E., Chakraborty, R., Rothhammer, F., Carvallo, P., Mitochondrial DNA polymorphisms in Chilean aboriginal populations: Implications for the peopling of the southern cone of the continent (2000) Am. J. Phys. Anthropol., 113, pp. 19-29Morriss-Kay, G.M., Derivation of the mammalian skull vault (2001) J. Anat., 199, pp. 143-151Moss, M.L., Young, R.W., A functional approach to craniology (1960) Am. J. Phys. Anthropol., 18, pp. 281-292Neves, W.A., Powell, J.F., Ozolins, E.G., Extra-continental morphological affinities of Palli Aike, Southern Chile (1999) Interciencia, 24, pp. 258-263Oksanen, J., Kindt, R., Legendre, P., O'Hara, B., Simpson, G.L., Stevens, M.H.H., (2008) Vegan: Community Ecology Package, , http://cran.r-project.org, and. R package version 1.11-4Opperman, L.A., Cranial sutures as intramembranous bone growth sites (2000) Dev. Dyn., 219, pp. 472-485Opperman, L.A., Gakunga, P.T., Carlson, D.S., Genetic factors influencing morphogenesis and growth of sutures and synchondroses in the craniofacial complex (2005) Semin. Orthod., 11, pp. 199-208Pearson, O.M., Activity, climate, and postcranial robusticity implications for modern human origins and scenarios of adaptive change (2000) Curr. Anthropol., 41, pp. 569-207Pearson, O.M., Millones, M., Rasgos esqueletales de adaptación al clima y a la actividad entre los habitantes aborígenes de Tierra del Fuego (2005) Magallania, 33, pp. 37-50Peres-Neto, P.R., Jackson, D.A., How well do multivariate data sets match? the advantages of a Procrustean superimposition approach over the Mantel test (2001) Oecologia (Berl.)., 129, pp. 169-178Perez, S.I., Bernal, V., Gonzalez, P., Differences between sliding semilandmarks methods: Implications for shape analyses of human populations (2006) J. Anat., 208, pp. 769-784Perez, S.I., Bernal, V., Gonzalez, P., Morphological differentiation of aboriginal human populations from Tierra del Fuego (Patagonia): Implications for South American peopling (2007) Am. J. Phys. Anthropol., 133, pp. 1067-1079Perez, S.I., Bernal, V., Gonzalez, P., Evolutionary relationships among prehistoric human populations: An evaluation of facial morphometric data employing molecular based genealogies (2007) Hum. Biol., 79, pp. 25-50Pucciarelli, H.M., The effects of race, sex, and nutrition on craniofacial differentiation in rats. a multivariate analysis (1980) Am. J. Phys. Anthropol., 53, pp. 359-368Pucciarelli, H.M., Growth of the functional components of the rat skull and its alteration by nutritional effects. a multivariate analysis (1981) Am. J. Phys. Anthropol., 56, pp. 33-41Pucciarelli, H.M., Oyhenart, E.E., Effects of maternal food restriction during lactation on craniofacial growth in weanling rats (1987) Am. J. Phys. Anthropol., 72, pp. 67-75Rae, T.C., Strand Viearsdóttir, U., Jeffery, N., Steegmann Jr., A.T., Developmental response to cold stress in cranial morphology of Rattus: Implications for the interpretation of climatic adaptation in fossil hominins (2006) Proc. R. Soc. Lond. B, 273, pp. 2605-2610(2008) R: A Language and Environment for Statistical Computing., , http://www.R-project.org, R Development Core Team. R Foundation for Statistical Computing, Vienna, AustriaRelethford, J.H., Craniometric variation among modern human populations (1994) Am. J. Phys. Anthropol., 95, pp. 53-62Relethford, J.H., Apportionment of global human genetic diversity based on craniometrics and skin color (2002) Am. J. Phys. Anthropol., 118, pp. 393-398Relethford, J.H., Global patterns of isolation by distance based on genetic and morphological data (2004) Hum. Biol., 76, pp. 499-513Reznick, D.N., Shaw, F.H., Rodd, F.H., Shaw, R.G., Evaluation of the rate of evolution in natural populations of guppies (Poecilia reticlrlatri) (1997) Science, 275, pp. 1934-1937Riesenfeld, A., The effect of extreme temperatures and starvation on the body proportions of the rat (1973) Am. J. Phys. Anthropol., 39, pp. 427-460Riesenfeld, A., The role of body mass in thermoregulation (1981) Am. J. Phys. Anthropol., 55, pp. 95-99Roberts, D.F., Body weight, race and climate (1953) Am. J. Phys. Anthropol., 11, pp. 533-558Roff, D.A., The evolution of genetic correlations: An analysis of patterns (1996) Evolution, 50, pp. 1392-1403Rohlf, F.J., Relative warps analysis and an example of its application to Mosquito wings (1993) Contributions to Morphometrics., pp. 132-159. , Pp. in. L. F. Marcus, E. Bello. and. A. García-Valdecasas, eds. Monografías del Museo Nacional de Ciencias Naturales, MadridRohlf, F.J., (2007) Tps Serie Softwares, , http://life.bio.sunysb.edu/morph/, Available atRohlf, F.J., Slice, D.E., Extensions of the Procrustes Method for the optimal superimposition of landmarks (1990) Syst. Zool., 39, pp. 40-59Roseman, C.C., Detecting interregionally diversifying natural selection on modern human cranial form by using matched molecular and morphometric data (2004) Proc. Natl. Acad. Sci. USA, 101, pp. 12824-12829Rothhammer, F., Silva, C., Craniometrical variation among South American prehistoric populations: Climatic, altitudinal, chronological and geographic contributions (1990) Am. J. Phys. Anthropol., 82, pp. 9-17Ruff, C.B., Morphological adaptation to climate in modern and fossil hominids (1994) Yearb. Phys. Anthropol., 37, pp. 65-107Sardi, M.L., Ramirez-Rozzi, F.V., A cross-sectional study of human craniofacial growth (2005) Ann. Hum. Biol., 32, pp. 390-396Sardi, M.L., Ramírez-Rozzi, F., González-José, R., Pucciarelli, H.M., South Amerindian craniofacial morphology: Diversity and implications for Amerindian evolution (2005) Am. J. Phys. Anthropol., 128, pp. 747-756Scheuer, L., Black, S., (2000) Developmental Juvenile Osteology., , and. Academic Press, LondonSchluter, D., (2000) The Ecology of Adaptative Radiation., , Oxford Univ. Press, New YorkSheets, H.D., Mitchell, C.E., Why the null matters: Statistical tests, random walks and evolution (2001) Genetica, 112-113, pp. 105-125Slatkin, M., Gene flow and the geographic structure of natural populations (1987) Science, 236, pp. 787-792(1961) Climatic Database of Servicio Meteorológico Nacional, , http://www.smn.gov.ar/, SMN. Available atSober, E., Reconstructing the past (1988) Parsimony, Evolution and Inference., , The MIT Press, Cambridge, MASperber, G.H., (2001) Craniofacial Development., , BC Decker Inc., HamiltonSpicer, G.S., Morphological evolution of the Drosophila virilis species group as assessed by rate tests for natural selection on quantitative characters (1993) Evolution, 47, pp. 1240-1254Steegmann, A.T., Platner, W.S., Experimental cold modification of cranio-facial morphology (1968) Am. J. Phys. Anthropol., 28, pp. 17-30Steward, J.H., (1950) Handbook of South American Indians., , Pub. Official, WashingtonStynder, D.D., Ackermann, R.R., Sealy, J.C., Craniofacial variation and population continuity during the South African Holocene (2007) Am. J. Phys. Anthropol., 134, pp. 489-500Turelli, M., Gillespie, J.H., Lande, R., Rate tests for selection on quantitative characters during macroevolution and microevolution (1988) Evolution, 42, pp. 1085-1089Vitzthum, V.J., A number no greater than the sum of its parts: The use and abuse of heritability (2003) Hum. Biol., 75, pp. 539-558Weaver, T.D., Roseman, C.C., Stringer, C.B., Were neandertal and modern human cranial differences produced by natural selection or genetic drift? (2007) J. Hum. Evol., 53, pp. 135-145Wobst, M., Boundary conditions for paleolithic social systems: A simulation approach (1974) Am. Antiq., 39, pp. 147-178