Menéndez, López Sosa, and Solari: Variación morfológica y afinidades biológicas de esqueletos del Holoceno temprano/medio de Norteamérica y Sudamérica:

Introducción

Los estudios sobre la diversidad biológica de los individuos más tempranos de las Américas se remontan a casi 200 años atrás, cuando se encontraron y describieron los primeros esqueletos humanos de Sudamérica considerados antiguos (Lund, 1840, 1844; Ameghino, 1880; Hrdlička, Holmes, Willis, Wright y Fenner, 1912). Desde entonces, se han planteado varios interrogantes, tales como: dónde se originaron los ancestros de los individuos que llegaron a las Américas, cuándo llegaron, y cuáles fueron las dinámicas de expansión poblacional, por mencionar sólo algunos (Neves y Pucciarelli, 1989; Neves y Hubbe, 2005; Fagundes et al., 2008; González-José, Bortolini, Santos y Bonatto, 2008; Pucciarelli, Perez y Politis, 2010; de Azevedo et al., 2011; de Saint Pierre et al., 2012; Hubbe, Strauss, Hubbe y Neves, 2015; Raghavan et al., 2015; Skoglund et al., 2015; Posth et al., 2018; entre otros). Actualmente, gracias al aporte de distintas disciplinas, al creciente número de muestras arqueológicas disponibles, y a los avances tecnológicos de las últimas décadas, muchas de las preguntas pueden ser abordadas con mayor profundidad. Esto ha permitido avanzar en aspectos tales como el origen de la(s) población(es) ancestral(es) de las Américas, para el cual se considera que la mayor parte proviene del este de Asia (Moreno-Mayar et al., 2018a; Posth et al., 2018); y la antigüedad de los humanos en las Américas, que gracias al avance de estudios recientes se sostiene la presencia humana en el centro sur de Norteamérica antes del último máximo glacial (Ardelean et al., 2020; Boëda et al., 2021; Bennett et al., 2021; Madsen et al., 2022). Sin embargo, más allá de la gran atención que esta región ha recibido a lo largo de las últimas décadas, aún quedan muchos aspectos por ser estudiados e interpretados en mayor profundidad.

Uno de los aspectos que aún queda por dilucidar es aquel relacionado con las dinámicas de expansión de las poblaciones humanas a lo largo y a lo ancho de ambos continentes. Una vez llegados a Beringia, múltiples caminos podrían haber sido adoptados en la dispersión hacia el sur (Miotti, 2006; Achilli et al., 2008; Goebel, Waters y O’Rourke, 2008; Perego et al., 2010; Bodner et al., 2012; Reich et al., 2012; Achilli et al., 2013; Battaglia et al., 2013; Llamas et al., 2016). En Norteamérica, gracias al avance de estudios paleoambientales, arqueológicos y genéticos, existen dos hipótesis principales que explican la dispersión humana de esta amplia región: la ruta costera, a lo largo de la cual los primeros humanos se desplazaron a lo largo de la costa Pacífica y desde allí se expandieron tanto hacia al interior de Norteamérica como hacia Centroamérica y Sudamérica (Erlandson et al., 2007; Tamm et al., 2007; Fagundes et al., 2008; Bodner et al., 2012); y la ruta interior, según la cual se plantea que los desplazamientos se realizaron mediante el corredor libre de hielo entre las masas de hielo Laurentino y Cordillera (que cubrían lo que son actualmente las provincias de Alberta y el este de la Columbia Británica), cuya disponibilidad como pasaje se estima habría empezado hace 13.000 años AP (Pedersen et al., 2016; Potter et al., 2017). No se descarta la utilización de ambas rutas, lo cual inevitablemente queda supeditado a las condiciones climáticas que prevalecieron en diferentes momentos.

En base a expectativas arqueológicas y antropológicas previas, estudios genéticos recientes plantean que las expansiones humanas iniciales a Sudamérica habrían comenzado hace aproximadamente 16.000 años AP, luego de la divergencia entre los linajes que fueron descriptos como NNA/SNA (North Native Americans, South Native Americans), seguido de una expansión desde Mesoamérica que habría comenzado hace 8.000 años AP y que alcanzó Patagonia hace aproximadamente 5.000 años AP (Moreno-Mayar et al., 2018a; Posth et al., 2018). De acuerdo con interpretaciones de la evidencia arqueológica, lingüística, genética y osteológica, existen al menos dos modelos que explican la dispersión humana temprana en Sudamérica: a) una divergencia temprana en el noroeste de Sudamérica, a partir de la cual las poblaciones humanas se habrían bifurcado siguiendo, por un lado, una ruta andina a lo largo del océano Pacífico; y por otro lado, poblaciones más pequeñas habrían seguido un trayecto a lo largo de llanuras aluviales y/o a lo largo del océano Atlántico (Gruhn, 1994; Tarazona-Santos et al., 2001; Miotti, 2006; Pucciarelli et al., 2006; Wang et al., 2007; Rothhammer y Dillehay, 2009; Yang et al., 2010; Harris et al., 2018). La alternativa más fuerte a dicho modelo es b) una expansión inicial mediante una ruta costera a lo largo del Pacífico, a partir de la cual se habrían producido migraciones trasandinas hacia el este del continente (Wang et al., 2007; Bodner et al., 2012; Rothhammer, Fehren-Schmitz, Puddu y Capriles, 2017; Gómez-Carballa et al., 2018); aunque no se debería descartar c) una combinación de ambos modelos (Miotti, 2006; Menéndez, Rademaker y Harvati, 2019) (Figura 1). Luego, durante fines del Holoceno medio y comienzos del Holoceno tardío se estableció otra etapa de dinámicas expansivas que no serán abordadas en este trabajo, tales como el arribo de individuos de Mesoamérica desde hace aproximadamente 5.000 años (Moreno-Mayar et al., 2018a; Posth et al., 2018) y los movimientos de escala continental de grupos Tupí desde el Holoceno medio al tardío (Ramallo et al., 2013; Iriarte et al., 2017).

Figura 1

Modelos que explican la dispersión humana temprana en Sudamérica: (a) divergencia temprana en el noroeste de Sudamérica; (b) ruta costera a lo largo del Pacífico y migraciones trasandinas; (c) una combinación de los anteriores.

El análisis de la variación morfológica de los esqueletos humanos del Holoceno temprano y medio de las Américas puede contribuir al debate sobre las más antiguas dinámicas de expansión en el continente (Neves y Pucciarelli, 1989, 1991; Neves y Hubbe, 2005; Sardi, Rozzi, González-José y Pucciarelli, 2005; González-José et al., 2008; Perez, Bernal, Gonzalez, Sardi y Politis, 2009; Pucciarelli et al., 2010; de Azevedo et al., 2011; Manríquez et al., 2011; Hubbe et al., 2015; Menéndez et al., 2015; Delgado Burbano, 2017; von Cramon Taubadel, Strauss y Hubbe, 2017; Kuzminsky et al., 2018, entre otros). A partir del estudio de la variación craneométrica de los esqueletos humanos más tempranos del continente hallados a la fecha pueden discutirse posibles vías de dispersión con base en las similitudes y diferencias morfométricas entre individuos de distintas áreas (Pucciarelli et al., 2006, 2008). Se ha planteado previamente la existencia de una estructura geográfica este-oeste principalmente para poblaciones humanas recientes, donde las poblaciones del este fueron descritas presentando los niveles más bajos de variación intrapoblacional, neurocráneo de tamaño pequeño y grandes similitudes con los individuos más tempranos de Sudamérica (Tarazona-Santos et al., 2001; Pucciarelli et al., 2006). Asimismo, trabajos previos han demostrado que muchos de los individuos tempranos procedentes de distintas regiones presentan grandes similitudes entre sí, e.g., Sabana de Bogotá y Lagoa Santa (Neves et al., 2013; Hubbe et al., 2015), mientras que otros investigadores destacan la gran variación intrapoblacional e interpoblacional de las primeras poblaciones humanas que se dispersaron en las Américas (Gonzalez-José et al., 2008; Hubbe et al., 2020, entre otros). Sin embargo, aún quedan por analizar las causas de dicha variación morfológica, así como las similitudes y diferencias de individuos del Holoceno temprano/medio con el fin de interpretar dinámicas poblacionales de expansión.

El objetivo de este trabajo es analizar las afinidades morfológicas de algunos de los esqueletos más antiguos de Sudamérica con el fin de evaluar los modelos que describen patrones de expansión en el continente durante el Holoceno temprano y medio, es decir, momentos inmediatamente posteriores a las primeras expansiones humanas. Mediante el análisis de la variación morfológica los esqueletos humanos del Holoceno temprano y Holoceno medio inicial, discutiremos las dinámicas de expansión ocurridas previamente al cambio poblacional que tuvo lugar luego del Holoceno medio (“mid-Holocene population turnover”, sensuPosth et al., 2018). Pondremos a prueba el modelo según el cual se habría producido una divergencia temprana en el noroeste de Sudamérica, a partir de la cual las poblaciones humanas se habrían escindido siguiendo, ya sea una ruta andina o costera a lo largo del océano Pacífico y, por otro lado, un trayecto hacia el interior del continente y/o a lo largo del océano Atlántico. Si este modelo permite explicar la dispersión de las poblaciones humanas en Sudamérica, entonces esperamos encontrar una estructuración poblacional con muestras del oeste y este del continente que presenten las mayores diferencias entre sí, mientras que las muestras al interior de cada una de estas regiones deberían presentar mayores similitudes entre sí. Si, por el contrario, los resultados muestran similitudes entre individuos del oeste y este de Sudamérica, entonces dicho modelo no es suficiente para explicar las dispersiones humanas tempranas. Si este fuese el caso, el modelo que plantea una expansión inicial mediante una ruta costera a lo largo del Pacífico, a partir de la cual se habrían producido migraciones trasandinas hacia el este del continente, explicaría las primeras expansiones humanas en Sudamérica. Cabe destacar que aquí estamos evaluando el rango de ocupaciones inmediatamente posteriores a las expansiones más tempranas del continente, así como también los posibles movimientos ocurridos entre el Holoceno temprano y medio. Por lo tanto, no estamos considerando la existencia de migraciones más recientes ni de flujo génico entre regiones durante el Holoceno tardío, lo cual ha sido evaluado en detalle en otros trabajos, y dada su complejidad requiere un análisis pormenorizado entre las distintas regiones (García, Pauro, Nores, Bravi y Demarchi, 2012; Ramallo et al., 2013; Posth et al., 2018; Barberena et al., 2020; Postillone, Martínez, Flensborg y Dejean, 2020; entre otros).

Materiales y métodos

Muestra

Compilamos una muestra de 92 individuos adultos del Holoceno temprano (e.g., 11.700-8.200 años AP; Walker et al., 2012) y Holoceno medio (e.g., 8.200-4.200 años AP; Walker et al., 2012) procedentes del norte de Chile, Perú, pampas argentinas, Colombia, este de Brasil, México y Estados Unidos (Figura 2; Tabla 1). Los datos fueron obtenidos principalmente de la base de datos abierta recolectada por el Dr. Héctor Pucciarelli (Pucciarelli et al., 2006, 2008; Paschetta, González-José y Lanata, 2017), y la base de datos publicada por el Dr. Richard Jantz y el Dr. Douglas Owsley (Owsley y Jantz, 1999). Los nombres de las medidas (e.g., distancias lineales) registradas en ambos casos difieren dado que dichos investigadores tenían distintas preguntas de trabajo. Mientras que Jantz y Owsley estaban interesados en inferir la historia poblacional a partir de las muestras que estudiaron -para lo cual emplearon las variables definidas por W. W. Howells (Howells, 1995; Jantz y Owsley, 2001)-, Pucciarelli definió variables craneofuncionales, ya que su interés era evaluar la influencia de factores ambientales en la variación craneofacial (Pucciarelli, 2008; Sardi, 2017). Sin embargo, visto que, en algunos casos, las medidas morfométricas fueron definidas según los mismos puntos anatómicos, estas pueden considerarse homólogas a pesar de que tienen distintos nombres.

Tabla 1

Composición de las muestras utilizadas para el presente estudio e información cronológica asociada. La cronología se presenta de dos maneras: rangos cronológicos relativos basados en la estratigrafía indicados al comienzo con “~”, y fechados radiocarbónicos directos calibrados con un rango de dos sigmas indicados con “*”. En este último caso, se presenta el fechado entre paréntesis.

País/Región	Sitio arqueológico	Número de individuos	Sexo (M/F/I)	Cronología en años calibrados con dos sigma (cronología en años ¹⁴C AP)	Referencia bibliográfica
Estados Unidos	Plattsmouth Oss, Nebraska	2	1/1/0	~10,000-9,000	Owsley y Jantz (1999); Jantz y Owsley (2001)
	Wet Gravel Pit, Nebraska	2	1/1/0	~10,000-9,000	Key (1982); Owsley y Jantz (1999); Jantz y Owsley (2001)
	Gilder Mound, Nebraska	1	0/1/0	~10,000-9,000	Key (1982); Owsley y Jantz (1999); Jantz y Owsley (2001)
	Lime Creek, Nebraska	1	1/0/0	~10,000-9,000	Key (1982); Owsley y Jantz (1999); Jantz y Owsley (2001)
	Buhl, Idaho	1	0/1/0	12,772-12,474 (10,675 ± 95)*	Green et al. (1998); Lepper (2014)
	Wizards Beach, Nevada	1	1/0/0	10,516-10,245 (9,225 ± 60)*	Tuohy y Dansie (1997); Owsley y Jantz (1999); Lepper (2014)
	Spirit Cave 2, Nevada	1	1/0/0	10,720-10,571 (9,415 ± 25)*	Tuohy y Dansie (1997); Jantz y Owsley (2001); Lepper (2014)
México	Peñón de los Baños III	1	0/1/0	10,720-10,571 (10,755 ± 75)*	Gonzalez et al. (2003, 2015); Pucciarelli, González-José, Neves, Sardi y Rozzi (2008)
	Tlapacoya 1	1	1/0/0	12,104-11,606 (10,200 ± 65)*	Gonzalez et al. (2003, 2015); Pucciarelli, González-José, Neves, Sardi y Rozzi (2008)
	Metro Balderas	1	1/0/0	~10,500	Pucciarelli González-José, Neves, Sardi y Rozzi (2008); Gonzalez et al. (2015)
	Paleoméxico Cuenca 10-613652	1	1/0/0	~10,000-9,000	Pucciarelli González-José, Neves, Sardi y Rozzi (2008)
	Chimalhuacán	1	1/0/0	~10,500	Gonzalez et al. (2015); Pucciarelli González-José, Neves, Sardi y Rozzi (2008)
	Cueva del Tecolote	1	1/0/0	~9,000-7,000	Gonzalez et al. (2015); Pucciarelli González-José, Neves, Sardi y Rozzi (2008)
Colombia	Checua	9	7/2/0	9443-8307 (8,200 + 110 -6800 +)*	Groot (1992); Correal Urrego y Ardila (1979); Pucciarelli et al. (2006); Delgado Burbano (2017)
Colombia	Tequendama	10	6/4/0	8182-6480 (7269 + 65 - 5848 + 56)*	Correal Urrego y van der Hammen (1977); Pucciarelli et al. (2006); Delgado Burbano (2015, 2017)
Chile (Norte)	Camarones 14	5	3/2/0	8724-6621(7420 + 225 - 6615 + 390)*	Schiappacasse y Niemeyer (1984); Arriaza (1995); Pucciarelli et al. (2006)
Perú	Lauricocha	4	3/1/0	8530-5748 (7871 + 30 -5158 + 27)*	Cardich (1964); Pucciarelli et al. (2006); Fehren-Schmitz et al. (2015)
Perú	Paiján	1	1/0/0	11,945-8,452 (9,020 + 650)*	Chauchat (1975); Chauchat y Lacombe (1984); Chauchat, Ivez, Briceno y Uceda (1998); Pucciarelli et al. (2006)
Brasil (Este/Noreste)	Serra da Capivara	1	0/0/1	9526-9681 (8640 + 30)*	Menéndez et al. (2022)
Brasil (Este/Noreste)	Lagoa Santa	39	27/12/0	~10000-7000	Pucciarelli et al. (2006); Da-Gloria, Neves y Hubbe (2017)
Argentina (Pampas)	La Tigra	1	1/0/0	8177-7941 (7270 + 60)*	Politis, Barrientos y Stafford (2011); Menéndez et al. (2015)
	Laguna de los Pampas	1	1/0/0	10,238-9,762 (8971 + 77)*	Politis, Messineo, González, Álvarez y Favier Dubois (2012); Menéndez et al. (2015)
	Necochea	1	1/0/0	7938-7676 (7013 + 67)*	Politis y Bonomo (2011); Menéndez et al. (2015)
	Arroyo del Moro	2	1/1/0	7845-7571 (6885 + 73)*	Politis y Bonomo (2011); Menéndez et al. (2015)
	Arroyo Seco 2	2	1/1/0	8554-8370 (7700 + 60)*	Politis, Gutiérrez y Scabuzzo (2014); Menéndez et al. (2015)
	Arroyo Chocorí	1	1/0/0	7882-7681 (7010 + 60)*	Politis, Barrientos y Stafford (2011); Menéndez et al. (2015)
TOTAL		92	63/28/1

Figura 2

Mapa con la localización de las muestras analizadas.

En la Tabla 1 se detalla el número de individuos, información cronológica asociada, y la estimación del sexo para cada una de las muestras. Mientras que algunas muestras presentan fechados radiocarbónicos directos, otras presentan información del contexto arqueológico que permite asociarlas a un rango cronológico. En el presente trabajo, las muestras fueron agrupadas por sitio arqueológico y también por región geográfica (correspondiente con el territorio político del país donde hoy se encuentran), siguiendo el agrupamiento utilizado previamente por Menéndez et al. (2015) para esta misma base de datos. Dado el carácter fragmentario de las muestras del Holoceno temprano, hemos decidido analizar todo el conjunto disponible, a pesar de que los distintos tamaños muestrales puedan afectar la interpretación de los resultados.

Variables morfométricas

Se seleccionó un total de diez medidas lineales que describen tanto rasgos morfométricos de la bóveda craneana como del esqueleto facial (Tabla 2). Estas diez medidas se encontraban presentes en ambas bases de datos morfométricas, y fueron registradas siguiendo a Pucciarelli (2008) y, en el caso de Jantz y Owsley (2001), siguiendo a Howells (1973). Las medidas fueron seleccionadas buscando encontrar un equilibro para maximizar tanto la inclusión de la mayor cantidad de individuos del Holoceno temprano en la muestra, como el número de variables para describir de manera más completa la variación morfológica de esta. Sin embargo, en muchos casos nos enfrentamos con dos dificultades: 1) los cráneos estudiados no se encuentran completos y, por ende, no permiten registrar un mayor número de medidas; 2) la cantidad de variables equivalentes entre ambas bases de datos es baja. De todas maneras, consideramos que las diez medidas lineales seleccionadas representan una solución intermedia cuando se compara con otros trabajos que utilizaron una muestra similar, algunos empleando solo seis variables (Steele y Powell, 1992), y otros, 22 variables (Jantz y Owsley, 2001).

Tabla 2

Lista de las medidas lineales craneofaciales usadas para este estudio.

Abreviatura	Nombre de la medida (puntos craneométricos que la delimitan)
LN	Longitud neurocraneana (nasion-opistocraneo)
AN	Ancho neurocraneano (eurion-eurion)
AF	Ancho facial (cigion-cigion)
HNM	Altura neurocraneana media (basion-bregma)
HF	Altura facial (nasion-prostion)
HR	Altura respiratoria (nasion-nasoespinal)
AR	Ancho respiratorio (alar-alar)
HO	Altura óptica (punto supraorbitario medio-punto infraorbitario medio)
AO	Ancho óptico (dacrion-ectoconquio)
AA	Ancho alveolar (ectomolar-ectomolar)

No hemos realizado análisis del error interobservador dado que no existen individuos en común en ambas bases de datos para los cuales se pueda calcular si existen diferencias en el registro de datos. Sin embargo, a pesar de que la base de datos de los individuos de Estados Unidos fue registrada por otros observadores distintos al resto de la muestra, consideramos relevante incluir tales muestras como comparativas en nuestro trabajo. Por esto es por lo que seguimos el criterio previamente aplicado por Delgado Burbano (2015), quien también utilizó dicha base de datos comparativa después de haber constatado que se siguieron los mismos procedimientos de registro de las variables aquí analizadas. Asimismo, a pesar de que las muestras no se encuentran balanceadas según el sexo de los individuos (63 individuos masculinos, 28 femeninos, 1 indeterminado), seguimos trabajos previos que han sido analizado individuos de esta base de datos sin ningún tipo de corrección preliminar por sexo (Pucciarelli et al., 2006, 2010; Menéndez et al., 2015, 2019). Esta base de datos contiene la mayoría de los individuos más antiguos recuperados en las Américas, por lo que dicho sesgo se encuentra asociado con la visibilidad arqueológica en relación con los entierros del Pleistoceno final y del Holoceno temprano (Dillehay, 1997).

Las distancias lineales que fueron registradas por los observadores de ambas bases de datos constituyen variables de forma (form). La forma (form en inglés) está compuesta por forma (shape en inglés) y tamaño (centroid size) (Zelditch, Swiderski, Sheets y Fink, 2004). Dado que el objetivo de este trabajo es analizar la variación en forma (shape), calculamos las razones de Mosimann para poder obtener variables que excluyan la información de tamaño geométrico (i.e., media geométrica). Estas razones fueron calculadas como el cociente entre cada una de las distancias y la media geométrica de la totalidad de las distancias por cada individuo (Darroch y Mosimann, 1985; Jungers, Falsetti y Wall, 1995). Todos los siguientes análisis fueron realizados empleando las variables de forma (shape).

Análisis estadísticos

La variación morfológica de los individuos de las Américas fue explorada por medio de un análisis de componentes principales (ACP). Este método es ampliamente utilizado en estudios morfométricos, ya que permite reducir la dimensionalidad de los datos estudiados, al mismo tiempo que resume y ordena el porcentaje de la variación descrita por las diferentes variables (Mitteroecker y Gunz, 2009). Mediante el ACP se transforma el conjunto original de variables (aquí, las diez medidas craneométricas) en un nuevo conjunto de variables no correlacionadas entre sí que representan la mayor variación del conjunto de datos (Zar, 2010). Se realizó un ACP utilizando la totalidad de los individuos y agrupando las muestras por región geográfica, lo cual se corresponde con el actual territorio político de cada país (ver Tabla 1). Complementariamente, utilizando los scores de los diez componentes principales (CP) que resumen la variación total en la muestra se realizó un análisis de funciones discriminantes (AFD). El AFD es un análisis estadístico multivariado que se diferencia del ACP porque su propósito es encontrar un eje de combinación de las variables que permita maximizar la diferencia existente entre grupos definidos a priori (Zar, 2010). Asimismo, el AFD estima la distancia que hay entre cada individuo y el centroide de cada grupo usando distancias Mahalanobis, que permiten clasificar las muestras según las variables dentro de los grupos previamente establecidos. Ya que las muestras son asignadas al grupo más cercano, dicha clasificación nos permite evaluar afinidades entre individuos. Tomando en cuenta la heterogeneidad que existe en el tamaño muestral, llevamos a cabo el AFD especificando las probabilidades a priori para cada grupo proporcional a la frecuencia de individuos por grupo. De esta manera, presentamos los gráficos que muestran la dispersión de las muestras a lo largo de los dos primeros componentes principales, y de los dos primeros ejes discriminantes, así como también tablas que muestran la contribución de las variables utilizadas en los dos primeros componentes y ejes principales. Complementariamente, reportamos la matriz de la tasa de clasificación predicha para cada grupo, la cual resume las probabilidades posteriores esperadas para cada muestra.

Las afinidades biológicas entre los individuos tempranos de las Américas fueron evaluadas mediante distancias biológicas. Para esto calculamos distancias Mahalanobis, mismas que son ampliamente utilizadas para evaluar similitudes entre poblaciones, ya que consideran la posible correlación entre las variables aleatorias, además de no verse afectadas por la escala y la unidad de las variables (Mahalanobis, 1930). Las distancias Mahalanobis fueron representadas mediante un análisis de coordenadas principales (ACoP), cuyos resultados se presentan mediante un gráfico de dispersión de los dos primeros ejes (Davis, 1986). Para el cálculo de distancias y su representación visual, se agruparon las muestras por sitio arqueológico usando la media de las medidas morfométricas de los individuos agrupados según el sitio arqueológico, aunque para México y Estados Unidos utilizamos las medias del conjunto. Dado que el presente estudio apunta a discutir dinámicas de expansión en Sudamérica, las muestras de México y Estados Unidos son utilizadas como comparativas, y por eso decidimos mantener a los individuos agrupados en estas dos regiones.

El análisis de componentes principales, el análisis de coordenadas principales y el análisis discriminante fueron realizados con los programas PAST (v. 2.17, Hammer, Harper y Ryan, 2001) y JMP (v. 16, by the SAS Institute Inc., 2021), mientras que las distancias Mahalanobis fueron calculadas en R (R Core Team, 2021).

Resultados

Variación morfológica de los individuos tempranos de las Américas

El ACP realizado con la totalidad de los individuos agrupados por región geográfica muestra que los dos primeros componentes principales (CP) explican el 50% de la variación total de la muestra (Figura 3A). A pesar de que existe una gran superposición en el rango de variación, al observar la distribución de las muestras a lo largo del CP2 se puede ver que los individuos del este/noreste de Sudamérica (Brasil) se encuentran hacia el extremo positivo, mientras que los individuos provenientes de la región andina o el oeste del continente (Colombia, Perú, Chile) se encuentran hacia el extremo negativo de dicho rango (Figura 3A). Adicionalmente, las muestras de Norteamérica (México, Estados Unidos) se ubican hacia el centro de la distribución, aunque más cercanas a las del oeste de Sudamérica, mientras que las de Argentina se distribuyen a lo largo de toda la variación del CP2, y ocupan un lugar central en la distribución. El CP1 no muestra diferencias que puedan explicarse según los criterios aquí seguidos; otros factores que exceden a los objetivos de este trabajo contribuirían a explicar la variación a lo largo de este componente. Sin embargo, se destaca que algunas muestras presentan una amplia variación a lo largo del CP2 (Brasil, Colombia, Perú), mientras que otras presentan una variación más acotada (Argentina, Chile). Las muestras de Norteamérica (México, Estados Unidos) se encuentran hacia el extremo negativo del CP1, presentan un rango acotado de variación y una fuerte superposición. La Tabla 3 muestra la contribución que tiene cada una de las variables morfométricas en los componentes principales. Para el CP1, puede verse que las variables que mayormente contribuyen son el ancho facial (AF) y la altura neurocraneana media (HNM). Los individuos del este de Sudamérica (Brasil, y parte de los individuos de las pampas argentinas) presentan esqueletos faciales más angostos que los individuos del noroeste del continente (Colombia, Perú). Complementariamente, aunque en menor medida, los individuos del noroeste del continente presentan bóvedas craneanas más bajas en sentido superior-inferior, mientras que, en comparación, los individuos de Brasil, y parte de los de las pampas presentan bóvedas más altas.

Figura 3

Análisis de la variación morfológica de los individuos tempranos de las Américas, agrupados por región geográfica: A) Dispersión de las muestras a lo largo de los dos primeros componentes principales; B) Dispersión de las muestras a lo largo de los dos primeros ejes discriminantes.

Tabla 3

Contribución de las variables utilizadas para este estudio en cada componente principal (CP) o eje discriminante (ED); las dos primeras columnas refieren al ACP de todos los individuos agrupados por región; las columnas tercera y cuarta, al AD de todos los individuos agrupados por región. En negrita se indican las variables que mayormente contribuyen a maximizar la variación de la muestra (ACP) o de los grupos (AD).

	PCA todos los individuos por región		Análisis discriminante todos los individuos por región
Componente o eje Principal	1	2	1	2
LN	0.256	0.047	4.020	-0.561
AN	-0.111	0.234	-3.804	1.523
AF	-0.847	0.082	-0.914	-3.147
HNM	0.400	0.515	1.585	-0.924
HF	0.091	-0.757	-21.675	-9.869
HR	0.175	-0.293	8.071	22.379
AR	0.032	0.012	2.638	9.856
HO	0.047	0.002	-3568	7.995
AO	-0.003	0.091	0.666	-17.330
AA	-0.041	0.066	12.980	-9.923

La variación a lo largo de los dos primeros ejes discriminantes (ED) resume un 75% de la variación de la muestra, y presenta una moderada estructuración geográfica (Figura 3B). El ED1 diferencia las muestras de Sudamérica en sentido oeste-este. Las muestras de Colombia, Chile, y Perú se encuentran en el extremo negativo del ED2, mientras que, en el extremo positivo, se ubican las de Brasil (Figura 3B). Las variables que mayormente contribuyen a las diferencias a lo largo del ED1 son la altura facial, altura respiratoria, y ancho alveolar (Tabla 3). Las muestras de Sudamérica, a diferencia de las de Norteamérica, presentan esqueletos faciales más altos en sentido supero-inferior, maxilares más amplios y narices más altas. El ED2 separa a las muestras en sentido norte-sur; las muestras de Norteamérica (Estados Unidos, México) se ubican hacia el extremo negativo, mientras las muestras de Sudamérica se encuentran hacia el extremo positivo (Colombia, Chile, Perú). Algunos individuos de Perú, Argentina, y Brasil presentan afinidad con las muestras de Norteamérica. Las muestras de Argentina se encuentran en una posición central considerando ambos EDs; algunos individuos son similares a aquellos de Brasil y otros a las muestras del oeste del continente (Figura 3B). Algunos individuos de Colombia y México se encuentran más cercanos a la distribución de Brasil. Las variables que mayormente contribuyen a dicha diferenciación son la altura respiratoria y la altura óptica, de manera que las muestras del oeste del continente presentan narices y órbitas más bajas que las de los individuos de Brasil (Tabla 3). La matriz de clasificación muestra que los grupos de Brasil, Colombia, Chile y Estados Unidos son los más homogéneos; los dos últimos son los que muestran tasas de clasificación correcta más altas (1.000 y 0.889, respectivamente). En contraste, las muestras de México, Perú, y Argentina presentan una mayor variación intramuestral, con valores entre 0.667 y 0.375 (Tabla 4). Esto confirma los resultados previos, de acuerdo con los cuales algunos individuos de Argentina presentan grandes afinidades con individuos de Colombia y Brasil, mientras que las muestras de México presentan afinidades morfológicas con las de Brasil.

Tabla 4

Matriz clasificatoria del análisis discriminante. Resumen de la tasa de clasificación predicha para cada grupo. Se destacan en negrita los valores mencionados en el texto.

Real	Tasa predicha
	Argentina	Brasil	Chile	Colombia	EUA	México	Perú
Argentina	0.375	0.250	0.125	0.250	0.000	0.000	0.000
Brasil	0.025	0.875	0.000	0.000	0.025	0.050	0.025
Chile	0.000	0.000	1.000	0.000	0.000	0.000	0.000
Colombia	0.053	0.105	0.000	0.842	0.000	0.000	0.000
EEUU	0.000	0.000	0.000	0.000	0.889	0.111	0.000
México	0.000	0.167	0.000	0.167	0.000	0.667	0.000
Perú	0.000	0.200	0.000	0.200	0.200	0.000	0.400

Cabe destacar que las interpretaciones de los resultados del análisis de componentes principales y análisis discriminante pueden verse afectadas por el tamaño muestral diferencial.

Afinidades biológicas entre los individuos tempranos de las Américas

Los resultados de la evaluación de las afinidades biológicas entre individuos procedentes de sitios arqueológicos del Holoceno temprano y medio muestran que existe cierta asociación entre algunas muestras en función de cercanía geográfica, cronología, o las rutas de dispersión utilizadas (Tabla 5). Las distancias Mahalanobis permiten vislumbrar cierto patrón geográfico y cronológico. Los individuos más similares entre sí son aquellos de Arroyo Seco 2/Tequendama, Camarones/Arroyo del Moro, Lagoa Santa/Arroyo Seco 2, Tequendama/Checua; Tequendama/Camarones; Estados Unidos/México; Lagoa Santa/Tequendama; Tequendama/Arroyo del Moro, y Lauricocha/Estados Unidos (Tabla 5). Mientras que en la mayoría de los casos estas muestras presentan distancias geográficas medias (salvo Checua/Tequendama, que son muy próximas entre sí), las muestras de Lagoa Santa/Arroyo Seco 2 se encuentran en el este de Sudamérica; las de Checua/Tequendama, Camarones/Tequendama en el oeste del continente; y las muestras de Camarones/Arroyo del Moro se encuentran ambas en el Cono Sur (Tabla 5). No existen similitudes entre muestras muy distantes geográficamente, como podría ser Estados Unidos/Arroyo Seco 2. Asimismo, la muestra de Camarones/Arroyo del Moro y Tequendama/Arroyo del Moro presenta una superposición cronológica ~7800-7500 años AP, mientras Camarones/Tequendama presentan una superposición de ~8100-6600 años AP (Tabla 1); sin embargo, es interesante observar que, a pesar de que la muestra de Tequendama presenta similitudes con la de Arroyo Seco 2, estas no presentan un rango cronológico similar (Tequendama es más reciente). Asimismo, Lagoa Santa/Arroyo Seco 2 presentan una superposición cronológica alrededor de ~8500-8300 años AP (Tabla 1) y ambas se encuentran ubicadas en el este del continente, pero los individuos no parecen presentar similitudes morfológicas. Una gran cantidad de pares de muestras presentan diferencias significativas entre sí. Entre las muestras que más difieren entre sí se pueden destacar las de Laguna de los Pampas y Paiján, que presentan grandes diferencias con la mayoría de las otras. La muestra de Laguna de los Pampas presenta grandes diferencias con otras de la misma región: Necochea, Chocorí, Arroyo del Moro, Arroyo Seco 2, pero también con Camarones y Paiján. La muestra de Paiján presenta grandes diferencias con las muestras de Toca dos Coqueiros y Estados Unidos (Tabla 5). Dichas muestras no presentan superposiciones en los rangos cronológicos; algunas son del Holoceno medio, y otras, del temprano (Tabla 1). A lo largo de la primera coordenada del PCoA se puede ver cierta estructuración geográfica de acuerdo con la cual las muestras del oeste del continente (Paiján, Lauricocha, Tequendama, Checua) y las de Norteamérica (Estados Unidos, México) se encuentran hacia el extremo positivo, y las del este, hacia el extremo negativo (Laguna de los Pampas, Arroyo Seco 2, Lagoa Santa, Arroyo del Moro, Necochea, Toca dos Coqueiros), aunque existen algunas excepciones a este patrón (Camarones, La Tigra, Chocorí) (Figura 4). Particularmente, Paiján, Toca dos Coqueiros y Necochea son las más disímiles de toda la muestra, y se ubican en los extremos de la distribución (Figura 4).

Tabla 5

Distancias Mahalanobis (D2) de los individuos según el sitio arqueológico o región. Las D2 se indican por debajo de la diagonal, mientras que los valores de p asociados se encuentran por encima de ella. EEUU=Estados Unidos; MX=México; CAM=Camarones 14; CHE=Checua; LAU=Lauricocha; COQ=Toca dos Coqueiros; LS=Lagoa Santa; LT=La Tigra; LLP=Laguna de los Pampas; NE=Necochea; AM=Arroyo del Moro; AS=Arroyo Seco 2; CHO=Chocorí; PJ=Paiján; TQ=Tequendama. Se destacan en negrita los valores de distancias más altos y más bajos que se describen en el texto.

	EEUU	MX	CAM	CHE	LAU	COQ	LS	LT	LLP	NE	AM	AS	CHO	PJ	TQ
EEUU	0	0.862	0.028	0.143	0.753	0.572	0.196	0.002	0.000	0.000	0.004	0.064	0.000	0.000	0.222
MX	05.40	0	0.080	0.165	0.554	0.118	0.709	0.015	0.000	0.000	0.050	0.364	0.009	0.000	0.667
CAM	20.03	16.71	0	0.295	0.039	0.000	0.100	0.000	0.000	0.000	0.945	0.456	0.110	0.000	0.825
CHE	14.70	14.16	11.84	0	0.578	0.032	0.264	0.000	0.000	0.000	0.319	0.556	0.000	0.000	0.861
LAU	06.70	08.76	19.07	8.52	0	0.113	0.283	0.000	0.000	0.002	0.083	0.222	0.000	0.000	0.553
COQ	08.57	15.39	36.59	19.71	15.52	0	0.070	0.000	0.002	0.000	0.000	0.022	0.000	0.000	0.006
LS	13.51	7.16	15.98	12.31	12.02	17.17	0	0.001	0.000	0.000	0.143	0.977	0.000	0.000	0.772
LT	26.79	21.87	37.04	44.35	32.28	45.45	29.29	0	0.000	0.000	0.000	0.000	0.000	0.000	0.000
LLP	32.76	40.20	72.32	50.15	53.03	27.32	41.40	61.10	0	0.000	0.000	0.000	0.000	0.000	0.000
NE	41.91	30.03	30.92	32.51	26.78	74.78	40.57	45.71	105.93	0	0.046	0.000	0.000	0.000	0.005
AM	25.39	18.27	4.04	11.49	16.59	43.48	14.68	37.81	78.39	18.53	0	0.384	0.026	0.001	0.773
AS	17.46	10.91	9.81	8.74	13.02	20.75	3.16	35.53	56.51	40.83	10.65	0	0.011	0.000	0.973
CHO	37.43	23.46	15.63	35.82	38.31	56.52	30.86	30.58	96.79	38.90	20.26	22.76	0	0.000	0.035
PJ	62.67	47.98	36.35	34.83	46.18	65.01	43.60	79.95	106.96	51.49	28.12	32.37	34.83	0	0.000
TQ	13.01	7.59	5.87	5.41	8.77	24.49	6.48	31.89	56.62	24.95	6.47	3.30	19.39	33.48	0

Figura 4

Análisis de coordenadas principales que muestra las afinidades biológicas entre los individuos de las Américas en base a distancias Mahalanobis entre las muestras a lo largo de las dos primeras dimensiones.

Discusión

En este artículo se analizó la variación morfológica de un conjunto de individuos del Holoceno temprano y medio de las Américas, con el fin de evaluar las afinidades biológicas entre ellos para discutir posibles modos de expansión de los descendientes de los primeros humanos en Sudamérica. El linaje de grupos Nativos Sudamericanos (SNA: South Native Americans) divergió hace ~16,000 años AP, lo cual coincide con estimaciones de la llegada de los primeros humanos a Sudamérica basadas en la distribución de fechados directos calibrados (Prates, Politis y Perez, 2020; Moreno-Mayar et al., 2018b). Sin embargo, cabe destacar que, actualmente, el esqueleto más antiguo de Sudamérica reportado con fechado radiocarbónico directo proviene de la Cueva de Candonga, en Córdoba, Argentina (media 2 sigma: 12.272 años calendáricos AP; Cornero, Neves y Rivero, 2014), lo que significa que los esqueletos del Holoceno temprano que han sido hasta ahora recuperados -muchos de los cuales forman parte de las bases de datos que fueron aquí analizadas- no representan los descendientes directos de los primeros migrantes, ya que son ~200 generaciones más recientes. Es decir, los esqueletos más antiguos disponibles para ser estudiados pertenecen a individuos que ya estaban asentados en Sudamérica, por lo que su variación morfológica brinda información acerca de posibles dinámicas de expansión en el continente durante los primeros milenios. En este trabajo en particular, nos focalizamos en la variación morfológica de los individuos recuperados que vivieron durante el Holoceno temprano y medio.

Nuestros resultados muestran que existe una relativa estructuración geográfica este-oeste de las muestras estudiadas, lo cual en parte sustenta el modelo de divergencia poblacional temprana en el noroeste de Sudamérica como resultado de una separación en al menos dos grupos que siguieron distintas direcciones (Gruhn, 1994; Tarazona-Santos et al., 2001; Miotti, 2006; Pucciarelli et al., 2006; Wang et al., 2007; Rothhammer y Dillehay, 2009; Yang et al., 2010; Harris et al., 2018). Las similitudes morfológicas entre individuos del este de Sudamérica (Brasil, Argentina) y sus diferencias con muestras del oeste del continente (Perú, Colombia) -lo cual ha sido ilustrado con los resultados del ADF, ACP, y PCoA- provee sustento al modelo según el cual los primeros grupos humanos que llegaron al noroeste de Sudamérica se separaron en dos grupos que migraron de manera independiente, uno siguiendo la costa Pacífica y el otro siguiendo la costa Atlántica o el interior del continente. Es decir, la diversificación morfológica basal entre poblaciones del este y oeste de Sudamérica se remontaría a ~16,000 años AP (Moreno-Mayar et al., 2018a; Prates et al., 2020). En este trabajo hemos asumido la existencia de un pulso de expansión antes del Holoceno medio (Rothhammer y Dillehay, 2009; Moreno-Mayar et al., 2018a; Posth et al., 2018), pero no debe descartarse que la descripta diferenciación entre muestras del este y oeste del continente puede reflejar una contribución diferencial de varios linajes ancestrales (Pucciarelli et al., 2006). Asimismo, la gran variación morfológica de la muestra -fundamentalmente evidenciada a través de los individuos de Argentina- apoya trabajos previos que plantean que las poblaciones del este de Sudamérica presentan mayor variación intragrupal que las poblaciones del oeste del continente (Tarazona-Santos et al., 2001; Pucciarelli et al., 2006; Menéndez et al., 2015), y trabajos que proponen una diversidad morfológica alta para los individuos del Holoceno temprano de las Américas (Sardi et al., 2005; González-José et al., 2008; Hubbe et al., 2020) y particularmente para Brasil (Menéndez et al., 2022). Sin embargo, esta afirmación no se mantiene en todos los casos ya que las muestras de Brasil se caracterizan por una gran homogeneidad intrapoblacional, lo cual concuerda con lo descripto anteriormente (Hubbe et al., 2015).

La amplia variación morfológica que muestran los individuos provenientes de las pampas argentinas, las similitudes que presentan con individuos de áreas geográficas distantes e incluso del oeste del continente (Camarones, Tequendama), así como las grandes diferencias entre individuos de esta región (Laguna de los Pampas con Necochea, Arroyo del Moro, Chocorí, Arroyo Seco 2) complejiza la interpretación de dicho modelo. Es decir, si bien el modelo de divergencia temprana en el noroeste de Sudamérica explica la variación de la mayoría de los individuos analizados, no explica completamente la variación morfológica de los individuos del Holoceno temprano y medio aquí analizados. La gran variación que presentan las muestras de las pampas argentinas, ilustradas por los resultados del ACP y el PCoA donde algunos individuos muestran mayores afinidades entre sí (Arroyo del Moro/Arroyo Seco 2), otros con individuos de Brasil (Arroyo Seco 2/Lagoa Santa), mientras que otros con muestras de Chile (Arroyo del Moro/Camarones), apoya la hipótesis de que algunas poblaciones del Cono Sur derivan de dos linajes que se habrían expandido independientemente, uno Pacífico/Andino y otro Amazónico/Atlántico (sensuGómez-Carballa et al., 2018). Estos resultados coinciden con resultados craneométricos previos que muestran una gran variación morfológica en individuos del Holoceno temprano de las pampas argentinas, y su asociación con individuos tempranos de otras regiones distantes tales como Lagoa Santa (Pucciarelli et al., 2010; Menéndez et al., 2015, 2017), y estudios genéticos que muestran flujo génico entre individuos de las pampas y de los Andes Centrales (Muzzio et al., 2018). La variación morfológica que ilustran los resultados del ACP y ADF, donde las muestras de Pampa se encuentran en el centro de la distribución, y las muestras de Brasil y Colombia en los extremos, apoyan esta interpretación. Sin embargo, debe tenerse en cuenta que esta puede estar influenciada por el bajo tamaño muestral de la muestra aquí analizada de las pampas argentinas, que no permite poner a prueba estadísticamente dicha hipótesis ni llegar a interpretaciones conclusivas. Asimismo, otro factor que también debe tenerse en cuenta cuando se analizan los cambios en la morfología craneofacial es el impacto de factores no-aleatorios, los cuales aquí no fueron evaluados, pero ya han sido planteados previamente en muestras del Cono Sur de Sudamérica (Perez y Monteiro, 2009; Fabra y Demarchi, 2011; Perez et al., 2011; Menéndez, 2018).

Asimismo, las grandes diferencias morfológicas que encontramos entre el individuo de Laguna de los Pampas respecto de otros de la misma región pueden explicarse como resultado de que aquel representa un linaje basal aún más antiguo, mientras los otros individuos representan linajes más recientes o mestizajes entre unos y otros (Roca-Rada et al., 2021). El fechado más antiguo que presenta este sitio (es el cráneo más antiguo con fechado directo en Sudamérica), en comparación con los otros, provee sustento a dicha interpretación (media 2 sigma: 9997 años calendáricos AP; Politis et al., 2012). Complementariamente, nuestros resultados muestran que el individuo de Arroyo Seco 2 se encuentra en el centro de la distribución en cuanto a la variación morfológica de las muestras de las pampas argentinas, lo cual apoya estudios previos que describen a las muestras de Arroyo Seco 2 como mixtas con un componente de Lagoa Santa/Anzick, pero también con un componente común a las otras poblaciones sudamericanas del Holoceno medio (Posth et al., 2018; Roca-Rada et al., 2021). Esto sugiere que durante el Pleistoceno final al Holoceno medio habrían confluido en las pampas argentinas poblaciones provenientes de linajes independientes que habrían llegado a través del Atlántico, pero también poblaciones que habrían migrado como producto de expansiones trasandinas. Es decir, consideramos que en las pampas argentinas habría existido continuidad poblacional durante el Holoceno temprano y medio con un moderado influjo de individuos asociados a linajes diferentes (Politis, 2008; Messineo et al., 2019; Postillone et al., 2020). De esta manera, los resultados de nuestro trabajo podrían interpretarse según un modelo de expansiones de linajes ancestrales de distintas procedencias (Figura 1c) para el caso de las pampas argentinas.

Con relación a la variación de las muestras del oeste de Sudamérica, las grandes similitudes entre algunas muestras (Checua/Paiján, Tequendama/Lauricocha), y también con las muestras de Norteamérica (Estados Unidos/Lauricocha, Estados Unidos/Tequendama), puede deberse a que estas derivarían de un mismo linaje ancestral. Esto concuerda con trabajos previos, en los cuales, a partir de análisis integrales de datos genéticos y morfológicos, se planteó que la falta de diferenciación genética y la moderada diversificación morfológica de muestras de la sabana de Bogotá puede deberse a un proceso de asimilación entre poblaciones locales y poblaciones de regiones adyacentes (Delgado Burbano, Rodríguez, Kassadjikova y Fehren-Schmitz, 2021). Además, estudios genéticos muestran una asociación estrecha entre muestras de Sudamérica y el individuo de Anzick-1 del noroeste de Estados Unidos (Posth et al., 2018), lo cual indica que el individuo de Anzick-1 formaba parte del mismo linaje ancestral del cual divergieron las poblaciones de Sudamérica. Respecto del centro y sur de los Andes, en Perú, las moderadas diferencias que encontramos entre los individuos de Lauricocha y el individuo de Paiján pueden explicarse por el relativo aislamiento geográfico de los grupos humanos que vivían en la sierra y en la costa. Según estudios genéticos recientes, dicha diversificación entre poblaciones de la sierra y de la costa tiene una profundidad temporal que se remonta a, al menos, ocho milenios (Lindo et al., 2018; Harris et al., 2018; Nakatsuka et al., 2020). Las similitudes entre las muestras del oeste de Sudamérica pueden interpretarse como resultado de la expansión a lo largo de la costa del Pacífico (Kuzminsky, 2013; Kuzminsky et al., 2018), mientras que las similitudes de la muestra de Chile con individuos de Pampa (Arroyo del Moro/Camarones) apoyan la existencia de migraciones trasandinas que se describieron más arriba, desde la región del Pacífico hacia el interior del continente (Harris et al., 2018; Gómez-Carballa et al., 2018).

Finalmente, las muestras de Norteamérica presentan una variación morfológica moderada, con fuertes afinidades y una gran superposición entre sí. Además, en general, se encuentran mayormente asociadas a las muestras del oeste de Sudamérica. Creemos que este patrón de afinidades podría explicarse como resultado de la expansión ocurrida de Mesoamérica hacia Sudamérica durante el Holoceno medio (Moreno-Mayar et al., 2018a). Asimismo, la ausencia de la señal Austromelanesia entre los individuos del oeste de Sudamérica (de acuerdo con Moreno-Mayar et al., 2018a, por el momento únicamente presente en individuos de Lagoa Santa; aunque según Posth et al., 2018, probablemente atribuible a un sesgo metodológico en los análisis realizados por Moreno-Mayar et al., 2018a), y la divergencia temprana de los individuos del oeste de Sudamérica (Harris et al., 2018) podrían explicar la menor variación de los individuos de Norteamérica y del oeste de Sudamérica. Sin embargo, estas muestras son de tamaño pequeño y han sido agrupadas para fines comparativos, por lo que se espera que el análisis de la variación de los individuos de Norteamérica pueda explorarse en más detalle en futuros trabajos.

Para concluir, nuestros resultados apoyan que la morfología craneofacial es una aproximación útil para discutir dinámicas de dispersión humana en Sudamérica. Asimismo, muestran que la dispersión humana en Sudamérica es más compleja que lo que ha sido planteado en los recientes modelos. Si bien el modelo de divergencia temprana en el noroeste de Sudamérica explica gran parte de la estructuración geográfica entre muestras del este y oeste del continente, hay ciertas regiones como las pampas argentinas para las cuales dicho modelo no resulta válido. Cabe destacar que para esta región existe una gran cantidad de sitios arqueológicos con pocos individuos por sitio (a excepción de Arroyo Seco 2, para el cual, de todas maneras, analizamos un individuo solamente). Para las otras regiones (e.g., Lagoa Santa, Brasil; Tequendama, Colombia), existe una gran cantidad de individuos provenientes de un mismo sitio arqueológico. Es decir, la gran variación que encontramos en la muestra de las pampas argentinas puede deberse a que no contamos con una representación de una gran cantidad de individuos de un solo sitio arqueológico (en términos de conservación de cráneos completos), o a que no haya una gran cantidad de sitios representados en las muestras de las otras regiones (a excepción de Colombia, representada por dos sitios). Estudios futuros deberán mejorar el diseño muestral e incluir, por ejemplo, para Brasil, individuos de varios sitios arqueológicos de distintas regiones para poner a prueba la hipótesis de un modelo mixto (e.g., divergencia temprana y migraciones trasandinas) en pos de explicar la diversificación de los humanos en Sudamérica durante el Holoceno temprano y medio.

Agradecimientos:

Queremos agradecer a Richard Jantz por amablemente facilitarnos información asociada a la base de datos de Norteamérica que publicó previamente, a Mark Hubbe por compartir con nosotros un script de R que programó para el cálculo de las distancias Mahalanobis, y a dos revisores anónimos por sus observaciones que contribuyeron a mejorar una versión previa de este trabajo.

Referencias bibliográficas

Achilli, A., Perego, U. A., Bravi, C. M., Coble, M. D., Kong, Q. P., Woodward, S. R., Salas, A., Torroni, A., y Bandelt, H. J. (2008). The phylogeny of the four pan-American MtDNA haplogroups: implications for evolutionary and disease studies. PLoS One, 3(3), e1764.

Achilli, A., Perego, U. A., Lancioni, H., Olivieri, A., Gandini, F., Kashani, B. H., Battaglia, V., Grugni, V., Angerhofer, N., Rogers, M. P. y Herrera, R. J. (2013). Reconciling migration models to the Americas with the variation of North American native mitogenomes. Proceedings of the National Academy of Sciences, 110(35), 14308-14313.

Ameghino, F. (1880). La Antigüedad del Hombre en el Plata, II. París-Buenos Aires: Masson-Igon Hnos.

Ardelean, C. F., Becerra-Valdivia, L., Pedersen, M. W., Schwenninger, J. L., Oviatt, C. G., Macías-Quintero, J. I., Arroyo-Cabrales, J., Sikora, M., Ocampo-Díaz, Y. Z. E., Rubio-Cisneros, I. I., y Watling, J. G. (2020). Evidence of human occupation in Mexico around the Last Glacial Maximum. Nature, 584(7819), 87-92.

Arriaza, B. T. (1995). Beyond death: the Chinchorro mummies of ancient Chile Washington, D. C.: Smithsonian Institution Press.

Barberena, R., Menéndez, L., Le Roux, P. J., Marsh, E. J., Tessone, A., Novellino, P., Lucero, G., Luyt, J., Sealy, J., Cardillo, M. y Gasco, A. (2020). Multi-isotopic and morphometric evidence for the migration of farmers leading up to the Inka conquest of the southern Andes. Scientific Reports, 10(1), 1-12.

Battaglia, V., Grugni, V., Perego, U. A., Angerhofer, N., Gomez-Palmieri, J. E., Woodward, S. R., Achilli, A., Myres, N., Torroni, A. y Semino, O. (2013). The first peopling of South America: new evidence from Y-chromosome haplogroup Q. PloS One, 8(8), e71390.

Bennett, M. R., Bustos, D., Pigati, J. S., Springer, K. B., Urban, T. M., Holliday, V. T., Reynolds, S. C., Budka, M., Honke, J. S., Hudson, A. M. y Fenerty, B. (2021). Evidence of humans in North America during the Last Glacial Maximum. Science, 373(6562), 1528-1531.

Bodner, M., Perego, U. A., Huber, G., Fendt, L., Röck, A. W., Zimmermann, B., Olivieri, A., Gómez-Carballa, A., Lancioni, H., Angerhofer, N. y Bobillo, M. C. (2012). Rapid coastal spread of First Americans: Novel insights from South America’s Southern Cone mitochondrial genomes. Genome Research, 22(5), 811-820.

Boëda, E., Gruhn, R., Vialou, A. V., Aschero, C., Vialou, D., Pino, M., Gluchy, M., Pérez, A. y Ramos, M. P. (2021). The Chiquihuite Cave, a Real Novelty? Observations about the Still-ignored South American Prehistory. PaleoAmerica, 7(1), 1-7.

Cardich, A. (1964). Lauricocha. Fundamentos para una Prehistoria de los Andes Centrales. Studia Praæhistorica. Revista Andina 1: 151-173. Centro Argentino de Estudios Prehistóricos, Buenos Aires.

Chauchat, C. (1975). The paijan complex, pampa de Cupisnique, Peru. Nawpa Pacha, 13(1), 85-96.

Chauchat, C. y J. -P. Lacombe. (1984). El Hombre de Paiján: el más antiguo peruano? Gaceta Arqueológica Andina 11, 4-6, 12.

Chauchat, C., lvez, C. G., Briceno, J. y Uceda, S. (1998). Sitios arqueológicos de la zona de Cupisnique y margen derecha del valle de Chicama, Patrimonio Arqueológico Zona Norte/4 Travaux de L‘lnstitut Francais d‘Etudes Andines, 113. Trujillo y Lima: Instituto Nacional de Cultura La Libertad, e Instituto Francés de Estudios Andinos.

Cornero, S., Neves, W. y Rivero, D. (2014). Nuevos aportes a la cronología de las ocupaciones tempranas en las sierras de Córdoba. Relaciones de la Sociedad Argentina de Antropología, 39, 285-292.

Correal Urrego, G. y Ardila, G. (1979). Investigaciones arqueológicas en abrigos rocosos de Nemocón y Sueva (Vol. 3). Bogotá: Fundación de Investigaciones Arqueológicas Nacionales; Banco de la República.

Correal Urrego, G. y Van der Hammen, T. (1977). Investigaciones arqueológicas en los abrigos de Tequendama Bogotá: Banco Popular.

Da-Gloria, P., Neves, W. A. y Hubbe, M. (Eds.). (2017). Archaeological and paleontological research in Lagoa Santa: the quest for the first Americans Chicago: Springer.

Darroch, J. N. y Mosimann, J. E. (1985). Canonical and principal components of shape. Biometrika, 72(2), 241-252.

Davis, J.C. (1986). Statistics and Data Analysis in Geology Nueva York: John Wiley & Sons,

De Azevedo, S., Nocera, A., Paschetta, C., Castillo, L., González, M. y González‐José, R. (2011). Evaluating microevolutionary models for the early settlement of the New World: the importance of recurrent gene flow with Asia. American Journal of Physical Anthropology, 146(4), 539-552.

de Saint Pierre, M., Bravi, C. M., Motti, J. M., Fuku, N., Tanaka, M., Llop, E., Bonatto, S. L. y Moraga, M. (2012). An alternative model for the early peopling of southern South America revealed by analyses of three mitochondrial DNA haplogroups. PLoS One , 7(9): e43486. doi: 10.1371/journal.pone.0043486

Delgado Burbano, M.E. (2015). Variación dental y craneofacial en el Norte de los Andes durante el Pleistoceno y el Holoceno (tesis de doctorado). Universidad Nacional de La Plata, La Plata, Argentina.

Delgado Burbano, M. (2017). Holocene population history of the Sabana de Bogotá region, Northern South America: An assessment of the craniofacial shape variation. American Journal of Physical Anthropology , 162(2), 350-369.

Delgado Burbano, M., Rodríguez, F., Kassadjikova, K. y Fehren-Schmitz, L. (2021). A paleogenetic perspective of the Sabana de Bogotá (Northern South America) population history over the Holocene (9000-550 cal BP). Quaternary International, 578, 73-86.

Dillehay, T.D., 1997. ¿Dónde están los restos óseos humanos del período pleistoceno tardío? Problemas y perspectivas en la búsqueda de los primeros americanos. Boletín de Arqueología PUCP, (1), 55-63.

Erlandson, J. M., Graham, M. H., Bourque, B. J., Corbett, D., Estes, J. A. y Steneck, R. S. (2007). The kelp highway hypothesis: marine ecology, the coastal migration theory, and the peopling of the Americas. The Journal of Island and Coastal Archaeology, 2(2), 161-174.

Fabra, M. y Demarchi, D. A. (2011). Geographic patterns of craniofacial variation in pre-hispanic populations from the southern cone of South America. Human Biology, 83(4), 491-507.

Fagundes, N. J., Kanitz, R., Eckert, R., Valls, A. C., Bogo, M. R., Salzano, F. M., Smith, D. G., Silva Jr, W. A., Zago, M. A., Ribeiro-dos-Santos, A. K. y Santos, S. E. (2008). Mitochondrial population genomics supports a single pre-Clovis origin with a coastal route for the peopling of the Americas. The American Journal of Human Genetics, 82(3), 583-592.

Fehren-Schmitz, L., Llamas, B., Lindauer, S., Tomasto-Cagigao, E., Kuzminsky, S., Rohland, N., Santos, F. R., Kaulicke, P., Valverde, G., Richards, S. M. y Nordenfelt, S. (2015). A re-appraisal of the early Andean human remains from Lauricocha in Peru. PloS One, 10(6), e0127141.

García, A., Pauro, M., Nores, R., Bravi, C. M. y Demarchi, D. A. (2012). Phylogeography of mitochondrial haplogroup D1: An early spread of subhaplogroup D1j from Central Argentina. American Journal of Physical Anthropology , 149(4), 583-590.

Goebel, T., Waters, M.R. y O’Rourke, D.H. (2008). The late Pleistocene dispersal of modern humans in the Americas. science, 319(5869), 1497-1502.

Gómez-Carballa, A., Pardo-Seco, J., Brandini, S., Achilli, A., Perego, U. A., Coble, M. D., Diegoli, T. M., Álvarez-Iglesias, V., Martinón-Torres, F., Olivieri, A. y Torroni, A. (2018). The peopling of South America and the trans-Andean gene flow of the first settlers. Genome Research , 28(6), 767-779.

Gonzalez, S., Jiménez-López, J. C., Hedges, R., Huddart, D., Ohman, J. C., y Turner, A. (2003). Earliest humans in the Americas: new evidence from Mexico. Journal of Human Evolution, 44, 379-387.

Gonzalez, S., Huddart, D., Israde-Alcántara, I., Domínguez-Vázquez, G., Bischoff, J. y Felstead, N. (2015). Paleoindian sites from the Basin of Mexico: evidence from stratigraphy, tephrochronology and dating. Quaternary International, 363(30), 4-19.

González‐José, R., Bortolini, M. C., Santos, F. R. y Bonatto, S. L. (2008). The peopling of America: craniofacial shape variation on a continental scale and its interpretation from an interdisciplinary view. American Journal of Physical Anthropology , 137(2), 175-187.

Green, T. J., Cochran, B., Fenton, T. W., Woods, J. C., Titmus, G. L., Tieszen, L., Davis, M. A. y Miller, S. J. (1998). The Buhl burial: a Paleoindian woman from southern Idaho. American Antiquity, 63(3), 437-456.

Groot de Mahecha, A. M. (1992). Checua: una secuencia cultural entre 8500 y 3000 años antes del presente Bogotá: Fundación de Investigaciones Arqueológicas Nacionales .

Gruhn, R. (1994). The Pacific Coast route of initial entry: an overview. En R. Bonnichsen y D. G. Steele (Eds.). Methods and Theory for Investigating the Peopling of the Americas (pp. 249-256). Corvallis: Oregon State University.

Hammer, Ø., Harper, D. A. y Ryan, P. D. (2001). PAST: Paleontological statistics software package for education and data analysis. Palaeontologia electronica, 4(1), 1-9.

Harris, D. N., Song, W., Shetty, A. C., Levano, K. S., Cáceres, O., Padilla, C., Borda, V., Tarazona, D., Trujillo, O., Sanchez, C. y Kessler, M.D. (2018). Evolutionary genomic dynamics of Peruvians before, during, and after the Inca Empire. Proceedings of the National Academy of Sciences , 115(28), E6526-E6535.

Howells, W. W. (1973). Cranial variation in man Volume 67. Cambridge, MA: Papers of the Peabody Museum of Archaeology and Ethnology, Harvard University.

Howells, W. W. (1995). Who’s who in skulls Volume 82. Cambridge, MA: Peabody Museum of Archaeology and Ethnology, Harvard University.

Hrdlička, A., Holmes, W. H., Willis, B., Wright, F. y Fenner, C. (1912). Early Man in South America Washington, DC: Bureau of American Ethnology, Bulletin no 52.

Hubbe, M., Strauss, A., Hubbe, A. y Neves, W.A. (2015). Early South Americans cranial morphological variation and the origin of American biological diversity. PLoS One , 10(10), e0138090.

Hubbe, M., Terrazas Mata, A., Herrera, B., Benavente Sanvicente, M. E., González González, A., Rojas Sandoval, C., Avilés Olguín, J., Acevez Núñez, E. y von Cramon-Taubadel, N., (2020). Morphological variation of the early human remains from Quintana Roo, Yucatán Peninsula, Mexico: Contributions to the discussions about the settlement of the Americas. PloS One, 15(1), e0227444.

Iriarte, J., Smith, R.J., Gregorio de Souza, J., Mayle, F. E., Whitney, B. S., Cárdenas, M. L., Singarayer, J., Carson, J. F., Roy, S. y Valdes, P. (2017). Out of Amazonia: Late-Holocene climate change and the Tupi-Guarani trans-continental expansion. The Holocene, 27(7), 967-975.

Jantz, R. L. y Owsley, D. W. (2001). Variation among early North American crania. American Journal of Physical Anthropology , 114(2), 146-155.

Jungers, W. L., Falsetti, A. B. y Wall, C. E. (1995). Shape, relative size, and size‐adjustments in morphometrics. American Journal of Physical Anthropology , 38(S21), 137-161.

Key, P.J. (1982). Craniometric relationships among Plains Indians: Culture-historical and evolutionary implications (tesis doctoral). The University of Tennessee, EEUU.

Kuzminsky, S. C. (2013). Human craniofacial variation among Pacific Rim Populations through the Holocene: A test of the coastal migration hypothesis using 3D morphometric methods (tesis de doctorado). University of California, Santa Barbara, EE.UU.

Kuzminsky, S. C., Reyes Báez, O., Arriaza, B., Méndez, C., Standen, V. G., San Román, M., Muñoz, I., Durán Herrera, A. y Hubbe, M. (2018). Investigating cranial morphological variation of early human skeletal remains from Chile: A 3D geometric morphometric approach. American Journal of Physical Anthropology , 165(2), 223-237.

Lepper, B. T. (2014). The people who peopled America. En D. W. Owsley y R. L. Jantz (Eds.). Kennewick Man: The Scientific Investigation of an Ancient American Skeleton (pp. 7-29). College Station: Texas A&M University Press.

Lindo, J., Haas, R., Hofman, C., Apata, M., Moraga, M., Verdugo, R. A., Watson, J. T., Llave, C. V., Witonsky, D., Beall, C. y Warinner, C. (2018). The genetic prehistory of the Andean highlands 7000 years BP though European contact. Science Advances, 4(11), eaau4921.

Llamas, B., Fehren-Schmitz, L., Valverde, G., Soubrier, J., Mallick, S., Rohland, N., Nordenfelt, S., Valdiosera, C., Richards, S. M., Rohrlach, A. y Romero, M. I. B. (2016). Ancient mitochondrial DNA provides high-resolution time scale of the peopling of the Americas. Science Advances, 2(4), e1501385.

Lund, P. W. (1840). View of the fauna of Brazil previous to the last geological revolution. Magazine of Natural History, 4(1-8), 1-389.

Lund, P. W. (1844). Notice sur des ossements humains fossiles, trouvés dans une caverne du Brésil. Mémoires de la Société Royale des Antiquaires du Nord, 49-77.

Madsen, D. B., Davis, L. G., Rhode, D. y Oviatt, C. G. (2022). Comment on “Evidence of humans in North America during the Last Glacial Maximum”. Science, 375(6577), p.eabm4678.

Mahalanobis, P.C. (1930). On tests and measures of group divergence. Journal and Proceedings of Asiatic Society of Bengal (New series) 26, 541-588.

Manríquez, G., Moraga, M., Santoro, C., Aspillaga, E., Arriaza, B. T. y Rothhammer, F. (2011). Morphometric and mtDNA analyses of archaic skeletal remains from southwestern South America. Chungara, Revista de Antropología Chilena, 43(2), 283-292.

Menéndez, L. P. (2018). Moderate climate signature in cranial anatomy of late Holocene human populations from Southern South America. American Journal of Physical Anthropology , 165(2), 309-326.

Menéndez, L. P., Perez, S. I., Pucciarelli, H. M., Bonomo, M., Messineo, P. G., Gonzalez, M. E. y Politis, G. G. (2015). Early Holocene human remains from the Argentinean Pampas: cranial variation in South America and the American peopling. PaleoAmerica, 1(3), 251-265.

Menéndez, L. P., Bonomo, M., Messineo, P. G., González, M. E., Politis, G. G. y Pérez, S. I. (2017). Craniofacial Variation of Early Americans from Argentinean Pampas: Implications for the Settlement of the New World. En C. Paschetta, R. González-José, y J. L. Lanata (Eds.). De cómo cruzar fronteras en la ciencia. Homenaje a Héctor Pucciarelli (pp. 87-109). Puerto Madryn y Bariloche: CONICET.

Menéndez, L. P., Rademaker, K. y Harvati, K. (2019). Revisiting east-west skull patterns and the role of random factors in South America: Cranial reconstruction and morphometric analysis of the facial skeleton from Cuncaicha Rockshelter (southern Peru). PaleoAmerica, 5(4), 315-334.

Menéndez, L. P., López-Sosa, M. C., da Silva, S. F. S. M., Martin, G., Pessis, A. M., Guidon, N. G. y Solari, A. (2022). Morphometric affinities and first radiocarbon dating of the Toca dos Coqueiros’ skull (Serra da Capivara, Brazil). Scientific Reports, 12 , 7807.

Messineo, P. G., Scheifler, N. A., Álvarez, M. C., González, M. E., Pal, N., Barros, M. P. y Politis, G. G. (2019). Was the Central Pampean Dunefields of Argentina Occupied during the Late Pleistocene? A Reappraisal of the Evidence. PaleoAmerica, 5(4), 378-391.

Miotti, L. (2006). La fachada atlántica, como puerta de ingreso alternativa de la colonización humana de América del Sur durante la transición Pleistoceno/Holoceno. En II Simposio Internacional El Hombre Temprano en América (pp. 155-188). México: Instituto Nacional de Antropología e Historia.

Mitteroecker, P. y Gunz, P. (2009). Advances in geometric morphometrics. Evolutionary Biology, 36(2), 235-247.

Moreno-Mayar, J.V., Vinner, L., de Barros Damgaard, P., De La Fuente, C., Chan, J., Spence, J. P., Allentoft, M. E., Vimala, T., Racimo, F., Pinotti, T. y Rasmussen, S. (2018a). Early human dispersals within the Americas. Science, 362(6419), eaav2621.

Moreno-Mayar, J.V., Potter, B. A., Vinner, L., Steinrücken, M., Rasmussen, S., Terhorst, J., Kamm, J. A., Albrechtsen, A., Malaspinas, A. S., Sikora, M. y Reuther, J. D. (2018b). Terminal Pleistocene Alaskan genome reveals first founding population of Native Americans. Nature, 553(7687), 203-207.

Muzzio, M., Motti, J. M. B., Paz Sepulveda, P. B., Yee, M. C., Cooke, T., Santos, M. R., Ramallo, V., Alfaro, E. L., Dipierri, J. E., Baillet, G., Bravi, C. M., Bustamante, C. D. y Kenny, E. E. (2018). Population structure in Argentina. PLoS ONE, 13, e0196325.

Nakatsuka, N., Lazaridis, I., Barbieri, C., Skoglund, P., Rohland, N., Mallick, S., Posth, C., Harkins-Kinkaid, K., Ferry, M., Harney, É. y Michel, M. (2020). A paleogenomic reconstruction of the deep population history of the Andes. Cell, 181(5), 1131-1145.

Neves, W. A. y Hubbe, M., 2005. Cranial morphology of early Americans from Lagoa Santa, Brazil: Implications for the settlement of the New World. Proceedings of the National Academy of Sciences , 102(51), 18.309-18.314.

Neves, W. A. y Pucciarelli, H. M. (1989). Extra-continental biological relationships of early South American human remains: A multivariate analysis. Ciência e Cultura, 41, 566-575.

Neves, W. A. y Pucciarelli, H. M. (1991). Morphological affinities of the first Americans: an exploratory analysis based on early South American human remains. Journal of Human Evolution , 21(4), 261-273. https://doi.org/10.1016/0047-2484(91)90107-7

Neves, W. A., Hubbe, M., Bernardo, D., Strauss, A., Araujo, A. y Kipnis, R. (2013). Early human occupation of Lagoa Santa, Eastern Central Brazil: craniometric variation of the initial settlers of South America. En K. E. Graf, C. V. Ketron y M. R. Waters (Eds.), Paleoamerican Odyssey (pp. 397-412). Texas: Texas A&M University Press.

Owsley, D. W. y Jantz, R. L. (1999). Databases for Paleo-American skeletal biology research. En R. Gruhn y R. Bonnichsen (Eds.). Who were the first Americans: Proceedings of the 58th Annual Biology Colloquium (pp. 79-96). Texas: Texas A & M University Press.

Paschetta, C., González-José, R. y Lanata, J. L. (Comps.) (2017). De cómo cruzar fronteras en la ciencia. Homenaje a Héctor M. Pucciarelli Puerto Madryn y San Carlos de Bariloche: Instituto Patagónico de Ciencias Sociales y Humanas e Instituto de Investigaciones en Diversidad Cultural y Procesos de Cambio.

Pedersen, M. W., Ruter, A., Schweger, C., Friebe, H., Staff, R. A., Kjeldsen, K. K., Mendoza, M. L., Beaudoin, A. B., Zutter, C., Larsen, N. K. y Potter, B. A. (2016). Postglacial viability and colonization in North America’s ice-free corridor. Nature, 537(7618), 45-49.

Perego, U. A., Angerhofer, N., Pala, M., Olivieri, A., Lancioni, H., Kashani, B. H., Carossa, V., Ekins, J. E., Gómez-Carballa, A., Huber, G. y Zimmermann, B. (2010). The initial peopling of the Americas: a growing number of founding mitochondrial genomes from Beringia. Genome Research , 20(9), 1174-1179.

Perez, S. I. y Monteiro, L. M. (2009) Nonrandom factors in modern human morphological diversification: a study of craniofacial variation in southern South American populations. Evolution, 63, 978-993.

Perez, S. I., Bernal, V., Gonzalez, P.N., Sardi, M. y Politis, G.G. (2009). Discrepancy between cranial and DNA data of early Americans: implications for American peopling. PLoS One , 4(5), e5746.

Perez, S. I., Lema, V., Diniz‐Filho, J. A. F., Bernal, V., Gonzalez, P.N., Gobbo, D. y Pucciarelli, H.M. (2011). The role of diet and temperature in shaping cranial diversification of South American human populations: an approach based on spatial regression and divergence rate tests. Journal of Biogeography, 38(1), 148-163.

Politis, G. (2008). The pampas and campos of South America. En H. Silverman y W. Isbell (Eds.), The Handbook of South American Archaeology (pp. 235-260). Nueva York: Springer.

Politis, G. y Bonomo, M. (2011). Nuevos datos sobre el “Hombre Fósil de Ameghino”. Publicación Especial APA, 12, 101-135.

Politis, G., Barrientos, G. y Stafford, T. (2011). Revisiting Ameghino: new 14c dates from ancient human skeletons from the Argentine Pampas. En D. Vilaou (Ed.), Peuplements et préhistorie en Amériques (pp. 43-54). París: Editorial du CTHS.

Politis, G. G., Messineo, P. G., González, M. E., Álvarez, M. C. y Favier Dubois, C. (2012). Primeros resultados de las investigaciones en el sitio Laguna de los Pampas (partido de Lincoln, provincia de Buenos Aires). Relaciones de la Sociedad Argentina de Antropología , 37(2), 463-472.

Politis, G. G., Gutiérrez, M. y Scabuzzo, C. (2014). Estado actual de las Investigaciones en el sitio Arroyo Seco 2 (Partido de Tres Arroyos, Provincia de Buenos Aires, Argentina) Olavarría: Serie Monográfica N° 5. INCUAPA-CONICET, UNICEN.

Posth, C., Nakatsuka, N., Lazaridis, I., Skoglund, P., Mallick, S., Lamnidis, T. C., Rohland, N., Nägele, K., Adamski, N., Bertolini, E. y Broomandkhoshbacht, N. (2018). Reconstructing the deep population history of Central and South America. Cell, 175(5), 1185-1197.

Postillone, M. B., Martínez, G., Flensborg, G. y Dejean, C. B. (2020). First analysis of mitochondrial lineages from the eastern Pampa-Patagonia transition during the final late Holocene. American Journal of Physical Anthropology , 171(4), 659-670.

Potter, B. A., Reuther, J. D., Holliday, V. T., Holmes, C. E., Miller, D. S. y Schmuck, N. (2017). Early colonization of Beringia and Northern North America: Chronology, routes, and adaptive strategies. Quaternary International, 444, 36-55.

Prates, L., Politis, G. G. y Perez, S. I. (2020). Rapid radiation of humans in South America after the last glacial maximum: A radiocarbon-based study. PLoS One , 15(7), e0236023.

Pucciarelli, H. M. (2008). Evolución y diversificación biológica humana desde la perspectiva craneofuncional México: Universidad Nacional Autónoma de México (UNAM), Instituto de Investigaciones Antropológicas-Escuela Nacional de Antropología e Historia (ENAH).

Pucciarelli, H. M., Neves, W. A., González-José, R., Sardi, M. L., Rozzi, F. R., Struck, A. y Bonilla, M. Y. (2006). East-west cranial differentiation in pre-Columbian human populations of South America. Homo, 57(2), 133-150.

Pucciarelli, H. M., González-José, R., Neves, W. A., Sardi, M. L. y Rozzi, F. R. (2008). East-West cranial differentiation in pre-Columbian populations from Central and North America. Journal of Human Evolution , 54(3), 296-308.

Pucciarelli, H. M., Perez, S. I. y Politis, G. G. (2010). Early Holocene human remains from the Argentinean Pampas: additional evidence for distinctive cranial morphology of early South Americans. American Journal of Physical Anthropology , 143(2), 298-305.

R Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Recuperado de https://www.R-project.org/

Raghavan, M., Steinrücken, M., Harris, K., Schiffels, S., Rasmussen, S., DeGiorgio, M., Albrechtsen, A., Valdiosera, C., Ávila-Arcos, M.C., Malaspinas, A. S. y Eriksson, A. (2015). Genomic evidence for the Pleistocene and recent population history of Native Americans. Science, 349(6250), aab3884.

Ramallo, V., Bisso‐Machado, R., Bravi, C., Coble, M. D., Salzano, F. M., Hünemeier, T. y Bortolini, M. C. (2013). Demographic expansions in South America: enlightening a complex scenario with genetic and linguistic data. American Journal of Physical Anthropology , 150(3), 453-463.

Reich, D., Patterson, N., Campbell, D., Tandon, A., Mazieres, S., Ray, N., Parra, M. V., Rojas, W., Duque, C., Mesa, N. y García, L.F. (2012). Reconstructing native American population history. Nature, 488(7411), 370-374.

Roca-Rada, X., Politis, G. G., Messineo, P. G., Scheifler, N., Scabuzzo, C., González, M., Harkins, K. M., Reich, D., Souilmi, Y., Teixeira, J. C. y Llamas, B. (2021). Ancient mitochondrial genomes from the Argentinian Pampas inform the early peopling of the Southern Cone of South America. Iscience, 24(6), 102553.

Rothhammer, F. y Dillehay, T. D. (2009). The late Pleistocene colonization of South America: an interdisciplinary perspective. Annals of Human Genetics, 73(5), 540-549.

Rothhammer, F., Fehren‐Schmitz, L., Puddu, G. y Capriles, J. (2017). Mitochondrial DNA haplogroup variation of contemporary mixed South Americans reveals prehistoric displacements linked to archaeologically‐derived culture history. American Journal of Human Biology, 29(6), e23029.

Sardi, M.L. (2017). El método craneofuncional. En C. Paschetta, R. González-José, y J. L. Lanata (Eds.). De cómo cruzar fronteras en la ciencia. Homenaje a Héctor Pucciarelli (pp. 120-135). Puerto Madryn: CONICET.

Sardi, M. L., Rozzi, F. R., González‐José, R. y Pucciarelli, H. M. (2005). South Amerindian craniofacial morphology: diversity and implications for Amerindian evolution. American Journal of Physical Anthropology , 128(4), 747-756.

Schiappacasse, V. y H. Niemeyer. (1984). Descripción y Análisis Interpretativo de un Sitio Arcaico Temprano en la Quebrada de Camarones Santiago: Museo Nacional de Historia Natural. Publicación Ocasional 14.

Skoglund, P., Mallick, S., Bortolini, M. C., Chennagiri, N., Hünemeier, T., Petzl-Erler, M. L., Salzano, F. M., Patterson, N. y Reich, D. (2015). Genetic evidence for two founding populations of the Americas. Nature, 525(7567), 104-108.

Steele, D. G. y Powell, J. F. (1992). Peopling of the Americas: paleobiological evidence. Human biology, 64(3), 303-336.

Tamm, E., Kivisild, T., Reidla, M., Metspalu, M., Smith, D. G., Mulligan, C. J., Bravi, C. M., Rickards, O., Martinez-Labarga, C., Khusnutdinova, E. K. y Fedorova, S. A. (2007). Beringian standstill and spread of Native American founders. PloS One, 2(9), e829.

Tarazona-Santos, E., Carvalho-Silva, D. R., Pettener, D., Luiselli, D., De Stefano, G. F., Labarga, C. M., Rickards, O., Tyler-Smith, C., Pena, S. D. y Santos, F. R. (2001). Genetic differentiation in South Amerindians is related to environmental and cultural diversity: evidence from the Y chromosome. The American Journal of Human Genetics , 68(6), 1485-1496.

Tuohy, D. R. y Dansie, A. J. (1997). New information regarding early Holocene manifestations in the western Great Basin. Nevada Historical Society Quarterly, 40(1), 24-53.

von Cramon-Taubadel, N., Strauss, A. y Hubbe, M. (2017). Evolutionary population history of early Paleoamerican cranial morphology. Science Advances, 3(2), e1602289.

Walker, M. J., Berkelhammer, M., Björck, S., Cwynar, L. C., Fisher, D. A., Long, A. J., Lowe, J. J., Newnham, R. M., Rasmussen, S. O. y Weiss, H. (2012). Formal subdivision of the Holocene Series/Epoch: a Discussion Paper by a Working Group of INTIMATE (Integration of ice‐core, marine and terrestrial records) and the Subcommission on Quaternary Stratigraphy (International Commission on Stratigraphy). Journal of Quaternary Science, 27(7), 649-659.

Wang, S., Lewis Jr, C. M., Jakobsson, M., Ramachandran, S., Ray, N., Bedoya, G., Rojas, W., Parra, M. V., Molina, J. A., Gallo, C. y Mazzotti, G. (2007). Genetic variation and population structure in Native Americans. PLoS Genetics, 3(11), e185.

Yang, N. N., Mazieres, S., Bravi, C., Ray, N., Wang, S., Burley, M. W., Bedoya, G., Rojas, W., Parra, M. V., Molina, J. A., et al (2010). Contrasting patterns of nuclear and mtDNA diversity in Native American populations. Annals of Human Genetics, 74, 525-538.

Zar, J.H. (2010). Biostatistical analysis Upper Saddle River: Pearson Prentice-Hall.

Zelditch, M. L., Swiderski, D. L., Sheets, H. D. y Fink, W. L. (2004). Geometric morphometrics for biologists Londres: Elsevier Academic Press.

Notas

[1] Financiamiento: Este trabajo se encuentra financiado por la Fundación Alemana para la Investigación Científica (DFG: Deutsche Forschungsgemeinschaft) y fue realizado en el marco del proyecto “Human morphological diversification in the Argentinean Pampas: Implications for the Peopling of South America” (N° de Proyecto: 415489479).