INTRODUCCIÓN

La vinculación entre poblaciones de regiones distantes constituyó parte integral de la forma de vida prehispánica en el Área Centrosur Andina. Las relaciones incluían las esferas económicas (intercambio de bienes), políticas (alianzas, conquistas) y socioculturales (prácticas culturales compartidas o resignificadas a partir de las de otros pueblos).

Dentro de este universo de posibles vínculos, el tráfico de bienes ha sido ampliamente estudiado desde la arqueología. Las evidencias de intercambio se han hallado en sitios correspondientes a todos los períodos de ocupación, incluso de momentos tempranos como el Arcaico (11.000-2500 años AP) (Aschero 1984; Fernández Distel 1986) y el Formativo (2500-1050 años AP) (Yacobaccio 2012; Martel 2014), y corresponden a materiales procedentes de otras regiones. El tráfico habría implicado tanto el intercambio de objetos de prestigio, como de bienes para consumo cotidiano (Yacobaccio 2012).

El intercambio de bienes de consumo se encuentra ligado a la complementariedad de recursos entre distintas regiones puesto que, en el área andina, las diferentes zonas productivas se distribuyen a lo largo de grandes distancias (Nielsen 2000). Estas zonas, a su vez, se despliegan a lo largo de diferentes fajas altitudinales correspondientes a pisos ecológicos con recursos particulares. Las sociedades prehispánicas habrían desarrollado una diversidad de prácticas que les permitieron el acceso y la explotación de esas diferentes ecorregiones. La explotación podía realizarse en forma directa, con individuos de un grupo social o étnico involucrados directamente en la producción o extracción de recursos, o en forma indirecta, para obtener los recursos por comercio o intercambio. Este comercio también podía ser directo o estar mediado por grupos especializados en el tráfico, como los caravaneros del altiplano, que desarrollaron esta práctica para el transporte de bienes utilizando tropillas de animales (Nielsen 2000; Llagostera 2010).

Las diferentes estrategias pueden pensarse como parte de un sistema político sostenido a lo largo del tiempo. En este sentido, habría sido importante la definición de los recursos que se habrían de comerciar de acuerdo con las necesidades, la frecuencia de los viajes, la forma de organizar el tráfico (directo, indirecto mediante caravanas, etc.) y las rutas establecidas para la circulación. Es decir que el intercambio fue consecuencia de una estructura social que garantizaba su efectividad y continuidad (Nielsen 2000; Yacobaccio 2012).

Las evidencias de un sistema organizado de intercambio de bienes permiten adjudicarles a las sociedades andinas una gran movilidad, incluso luego de haberse producido el cambio hacia una forma de subsistencia agroalfarera (Olivera 2001). No obstante, es necesario diferenciar entre la movilidad asociada al tráfico de recursos y aquella relacionada con eventos migratorios. En este sentido, es importante tener en cuenta que el estudio de la movilidad puede ser abordado desde diversas perspectivas. Por ejemplo, la presencia de individuos migrantes en un sitio arqueológico puede evaluarse mediante la identificación de prácticas culturales particulares de regiones específicas, como las culinarias o los rituales funerarios (e.g., Goldstein y Owen 2001), etc. Sin embargo, estos indicadores constituyen evidencias indirectas de migración en una comunidad, puesto que las prácticas culturales pueden tener diferentes interpretaciones, que incluyen otros procesos además de los migratorios.

Algunas evidencias más directas pueden hallarse a partir de análisis bioarqueológicos. Por ejemplo, mediante el estudio de prácticas corporales ha podido relacionarse la deformación craneana con la procedencia étnica de individuos en zonas donde habría ocurrido la articulación entre diferentes grupos (Blom 2005). Otra forma de evaluar con mayor precisión la procedencia de un individuo es analizar isótopos de estroncio, puesto que la señal isotópica de una región geográfica determinada se incorpora a los tejidos óseos y dentales de los individuos mediante su dieta y, por lo tanto, brinda información sobre el lugar de origen o sobre los lugares de residencia en los años previos a la muerte, según el elemento analizado (Manzanilla 2007; Tung y Knudson 2008). Además, algunos estudios han intentado establecer si las relaciones de intercambio podrían haber favorecido el flujo génico entre las poblaciones implicadas a partir de análisis craneométricos (e.g., Blom et al. 2004; Varela et al. 2014).

Una de las formas más directas de evaluar el aporte de los migrantes al acervo génico de una población, es decir, el flujo génico, es por medio de estudios de ADN antiguo (ADNa). Algunos procesos socioculturales de las poblaciones humanas pueden tener influencia sobre sus características a nivel genético (Jobling 2012). En particular, los patrones de movilidad pueden relacionarse con el flujo génico y la diferenciación entre distintos grupos. La posibilidad de que el sistema de intercambios también haya implicado el movimiento de personas de una población a otra a partir de nexos matrimoniales ha sido planteada para los Andes circumpuneños (Llagostera 2010: 292). Sin embargo, es importante resaltar que el flujo génico no necesariamente debió ocurrir mediante los mismos mecanismos involucrados en otros tipos de movilidad, como el intercambio de bienes. Por ejemplo, en términos hipotéticos, una interacción como el caravaneo, efectuado mayormente por individuos masculinos (Nielsen 2000), pudo haber conectado distintas poblaciones que, al mismo tiempo, mantenían flujo génico mediante la migración de mujeres.

Los avances realizados en los estudios de ADNa han permitido obtener datos de diversos sitios arqueológicos del Área Centrosur Andina, principalmente de ADN mitocondrial (ADNmt). A pesar de que este marcador permite evaluar únicamente la línea materna, debido a la gran cantidad de secuencias de ADNmt disponibles, los Andes Centrosur constituyen un área apropiada para intentar estudiar la problemática de la movilidad y el flujo génico desde una aproximación genética. Hasta el momento se cuenta con datos para sitios arqueológicos de los actuales territorios de Perú, altiplano de Bolivia, norte de Chile y noroeste argentino (e.g., Luciani et al. 2006; Kemp et al. 2009; Carnese et al. 2010; Mendisco et al. 2011, 2014; Baca et al. 2014; Fehren-Schmitz et al. 2014; Llamas et al. 2016; Russo et al. 2017). En algunas de estas poblaciones también se han analizado marcadores del cromosoma Y (línea paterna) y algunos autosómicos (e.g., Carnese et al. 2010; Baca et al. 2012; Mendisco et al. 2014, 2018; Fehren-Schmitz et al. 2015; Russo et al. 2016, 2018), pero esta información aún resulta escasa comparada con la disponible para ADNmt.

Con el objetivo de aportar al estudio de las interacciones entre poblaciones prehispánicas del Área Centrosur Andina, en este trabajo se analizaron la variabilidad y diferenciación genética mitocondrial, así como su asociación con la distancia geográfica entre grupos. La principal hipótesis subyacente a este objetivo fue que las poblaciones de diversas regiones conectadas por relaciones de intercambio han mantenido también flujo génico. Esto implica que aquellos grupos más conectados con otros tendrían mayor variabilidad y menor diferenciación genética. La hipótesis general puede dividirse en dos aspectos particulares: temporal y geográfico. En el primer caso, debido a que las dinámicas de intercambio se incrementaron hacia el período Tardío o de Desarrollos Regionales (1050 AP aproximadamente), es esperable que las poblaciones más tempranas hayan tenido un mayor grado de diferenciación y una menor variabilidad debido a la menor influencia de migrantes. En el segundo caso, se espera que las poblaciones más cercanas geográficamente hayan tenido un mayor grado de contacto y, por lo tanto, una menor diferenciación genética, siempre que ningún otro factor haya influenciado el flujo génico.

Materiales y Métodos

Grupos poblacionales estudiados

Se recopilaron secuencias de la Región Hipervariable I (RHVI) del ADNmt de individuos procedentes de 50 sitios arqueológicos del Área Centrosur Andina (Tabla 1). Además, se sumó una secuencia no publicada previamente (pero ver Russo 2016) de un individuo proveniente del sitio El Tero (Valle Calchaquí, Noroeste Argentino-NOA) del período de Desarrollos Regionales (PDR) (1050-520 años AP). El ADN de dicho individuo se extrajo a partir de una pieza dental y se amplificó por PCR la RHVI en tres fragmentos solapados que posteriormente fueron secuenciados. Las metodologías de extracción, amplificación y secuenciación, así como las medidas de prevención de la contaminación, se emplearon siguiendo un protocolo ya descripto (Russo et al. 2017).

Tabla 1

Grupos poblacionales prehispánicos del Área Centrosur Andina analizados en el presente trabajo.

La base de datos final comprendió un total de 415 individuos (mayormente adultos, de ambos sexos) provenientes de poblaciones que habitaron los actuales territorios del NOA, norte de Chile, Bolivia y Perú (Figura 1). Los individuos fueron agrupados en 18 poblaciones teniendo en cuenta la cercanía geográfica y temporal entre los sitios arqueológicos incluidos, e intentando que el tamaño muestral fuera lo más homogéneo posible (Tabla 1 y Figura 1).

Figura 1
Ubicación de los grupos poblacionales analizados en el presente trabajo. A modo de ejemplo, se nombran algunos de los sitios arqueológicos incluidos en cada grupo.

Si bien se procuró que los grupos poblacionales tuvieran una consistencia geográfica y temporal, en algunos casos quedaron conformados por individuos de sitios con temporalidades variadas. El principal motivo fue la poca disponibilidad de datos genéticos, sumado a la falta de información en bibliografía para asignar una temporalidad a cada individuo. Cabe mencionar que el grupo con el rango temporal más amplio (norte de Chile) posee solo un individuo de un período muy temprano, que fue incluido a fin de priorizar un tamaño muestral mayor (Tabla 1) y habiendo comprobado previamente que su exclusión no modificaba de forma sustancial los resultados obtenidos (Tablas S1 y S2).

Análisis estadísticos

Para cada grupo (Tabla 1 y Figura 1) se calcularon medidas de variabilidad ampliamente utilizadas, tales como el número de haplotipos (h), la diversidad haplotípica (Hd) (Nei 1987) y la diversidad nucleotídica (π) (Nei y Li 1979), que son proporcionales a la diversidad genética de la población analizada. Estos estadísticos fueron calculados con el programa DnaSP v. 5.10.01 (Librado y Rozas 2009).

La diferenciación genética entre los grupos estudiados se evaluó mediante un análisis molecular de la varianza (AMOVA) con el programa Arlequin v. 3.5.2.2 (Excoffier y Lischer 2010). Se calcularon los estadísticos FST entre pares de poblaciones seleccionando 1000 permutaciones para obtener la significación (α = 0,05). Se utilizó el modelo de evolución molecular de Tamura y Nei (1993) con un valor de 0,311 para el parámetro de la distribución Gamma, elegido con el programa MEGA6 (Tamura et al. 2013) acorde con el criterio BIC (Schwarz 1978).

Además, se analizó la correlación entre el número de poblaciones con las que se diferenció un determinado grupo y su diversidad nucleotídica utilizando el programa Infostat v. 2015 (Di Rienzo et al. 2015).

Finalmente, se estudió la relación entre la diferenciación genética y la distancia geográfica. Para ello, se obtuvo una matriz de distancias genéticas entre grupos (Tabla S3), medida como el número promedio de sustituciones por sitio entre pares de secuencias. Estas distancias se obtuvieron con el programa MEGA6 utilizando el mencionado modelo de evolución molecular. Por otro lado, se calcularon las distancias en kilómetros a partir de las coordenadas geográficas. A los grupos poblacionales que incluyeron distintos sitios se les asignó una coordenada geográfica correspondiente a un punto intermedio entre ellos (Tabla S4). Aquellos grupos que poseían iguales coordenadas geográficas pero diferente temporalidad (por ejemplo, los conjuntos Valles Altos y Costa del Sur de Perú; Tabla 1) fueron agrupados o analizados por separado dependiendo de si se encontraron diferencias genéticas entre ellos a partir del AMOVA. En efecto, debido a que los grupos del Horizonte Medio (HM) y del período Intermedio Tardío (PITar) de la costa sur de Perú se diferenciaron, no pudieron ser agrupados en un solo grupo llamado “Costa Sur de Perú”. Por lo tanto, se contó con dos matrices de distancias genéticas: en una de ellas se incluyó dentro del grupo “Costa Sur de Perú” a los individuos de todas las temporalidades menos a los del PITar, mientras que en la otra, ese grupo incluyó individuos de todos los períodos exceptuando el HM (Tabla S3). Utilizando la latitud y la longitud, la matriz de distancias geográficas entre pares de grupos (Tabla S5) se calculó con las funciones distm y distVincentyEllipsoid (Vincenty 1975) disponibles en el paquete geosphere (Hijmans 2015) de R v. 3.0.3 (R Core Team 2014). La correlación entre las matrices de distancias genéticas y geográficas se evaluó mediante la prueba de Mantel (Mantel 1967) utilizando Infostat.

RESULTADOS

Los estadísticos de variabilidad genética indicaron que, en líneas generales, los grupos con mayor tamaño muestral tuvieron más haplotipos (h) y aquellos con alta diversidad haplotípica también presentaron elevada diversidad nucleotídica (Figuras 2 y S1; Tabla S6). Es importante destacar que la correlación hallada entre el tamaño muestral y el número de haplotipos puede explicarse por razones probabilísticas y es esperable, al igual que la no correlación entre los estadísticos de variabilidad (Hd y π) y el tamaño muestral (Tabla S7). Este hecho permite analizar los resultados de diversidad genética obtenidos para cada grupo poblacional independientemente de su tamaño muestral.

Entre los grupos poblacionales con mayor variabilidad genética se encontraron las Tierras Altas del sur de Perú y Tiwanaku (Figura 2). Estos sitios están asociados a grandes centros ceremoniales o capitales estatales de los Imperios inca (520–414 años AP) y Tiwanaku (1450-800 años AP), respectivamente. Asimismo, los sitios incluidos en el norte de Chile y el sector sur de la Quebrada de Humahuaca (QH) presentaron valores relativamente altos de variabilidad genética (Figura 2).

Figura 2
Diversidad nucleotídica (π) en función de la diversidad haplotípica (Hd) para cada grupo poblacional estudiado en el presente trabajo. El tamaño de los puntos es proporcional al número de haplotipos (h). Los grupos Quebrada de Humahuaca Norte y Valle Calchaquí no se muestran por presentar valores muy extremos y para facilitar la visualización (Figura S1). HM: Horizonte medio. PITar: Período Intermedio tardío. HT: Horizonte temprano. PITem: Período Intermedio temprano.

En el extremo opuesto, dentro de los grupos con menor variabilidad, se encontraron el sitio Tompullo 2 y el sector norte de la QH, así como los sitios más tempranos de la costa sur de Perú.

Con respecto a la diferenciación genética, se observó un alto grado de estructuración considerando toda el área andina analizada (p < 0,00001; Tabla S8) e incluso se hallaron diferencias significativas entre pares de poblaciones dentro de una misma región (Figura 3). Por ejemplo, dentro del NOA, el sector norte de la QH se diferenció tanto del grupo sur como de Las Pirguas, si bien no se hallaron evidencias de diferenciación entre los grupos QH Sur, Las Pirguas y Valle Calchaquí (Figura 3).

Figura 3
Comparaciones de a pares (FST) entre los grupos poblacionales prehispánicos estudiados en el presente trabajo. *p<0,05. HM: Horizonte medio. PITar: Período Intermedio tardío. HT: Horizonte temprano. PITem: Período Intermedio temprano.

Cabe destacar que tanto el norte de Chile como Tiwanaku no se diferenciaron entre sí, ni de ninguno de los grupos del NOA (a excepción de QH Norte, que difirió de Tiwanaku). Sin embargo, ambos presentaron diferencias con respecto a muchos grupos del actual territorio de Perú (Figura 3), aunque Tiwanaku se diferenció de un menor número de poblaciones (Figura 3).

En líneas generales, la mayoría de los grupos localizados en el actual Perú se diferenciaron entre sí (Figura 3). Cabe destacar que en la costa sur del actual Perú se hallaron diferencias entre grupos muy cercanos pertenecientes a distintas temporalidades (Fehren-Schmitz et al. 2014) (Figura 3).

Las poblaciones más diferenciadas del resto fueron Tompullo 2, costa sur de Perú del Horizonte temprano (HT) y costa sur de Perú del período Intermedio temprano (PITem), seguidos de QH Norte y Las Pirguas (Figura 3). En el extremo opuesto se encontraron los grupos QH Sur, Valle Calchaquí y Tiwanaku (Figura 3). La cantidad de poblaciones de las cuales se diferenció un determinado grupo se correlacionó negativamente con su diversidad nucleotídica (r = -0,714; p = 0,00087). Es decir que, en líneas generales, los grupos que más se diferenciaron del resto presentaron valores más bajos de diversidad genética, y viceversa (Figura 4).

Figura 4
Cantidad de poblaciones de las cuales se diferencia (eje de la izquierda) y diversidad nucleotídica (eje de la derecha) de cada grupo poblacional analizado en el presente trabajo. HM: Horizonte medio. PITar: Período Intermedio tardío. HT: Horizonte temprano. PITem: Período Intermedio temprano.

Finalmente, no se hallaron evidencias de una relación entre las distancias genéticas y geográficas de los grupos poblacionales estudiados. Ninguna de las dos matrices de distancias genéticas (generadas para analizar por separado a las poblaciones de la costa sur de Perú HM y PITar) se correlacionó significativamente con la matriz de distancias geográficas (p = 0,08 y p = 0,07). Además, presentaron coeficientes de correlación lineal muy bajos (0,187 y 0,215).

DISCUSIÓN

En este trabajo se estudió la variabilidad genética de poblaciones humanas que habitaron el Área Centrosur Andina en tiempos prehispánicos a partir del análisis de secuencias de la RHVI del ADNmt, teniendo como hipótesis que la conexión entre diversas regiones por medio de relaciones de intercambio habría favorecido el flujo génico entre ellas. Bajo esta hipótesis general, las diferencias en los niveles hallados de diversidad y diferenciación genética entre las poblaciones estudiadas pueden explicarse por distintas causas.

En primer lugar, con respecto al factor temporal, se observó una tendencia a que las poblaciones más antiguas presentaran menos variabilidad que las tardías y se encontraran más diferenciadas del resto, como era esperable. Por ejemplo, la costa sur de Perú (períodos Intermedio temprano y Horizonte temprano, entre 2790 y 1310 años AP) podría incluir sitios de pequeña envergadura y que posiblemente habrían estado más aislados. Es decir que, en los momentos más tempranos, las poblaciones habrían tenido menor movilidad entre diferentes regiones y, por consiguiente, escaso flujo génico, al menos en comparación con períodos posteriores. Es importante destacar que la costa sur de Perú representa una de las muestras más numerosas en su conjunto, por lo que esta interpretación no estaría afectada por el tamaño muestral. En el NOA, si bien las relaciones interregionales se evidencian desde tiempos tempranos (e.g., Aschero 1984; Fernández Distel 1986; Yacobaccio 2012), es en el PDR cuando se despliega un importante tráfico de bienes, mediante, por ejemplo, las caravanas de llamas (Nielsen 2001). Este incremento en el contacto entre poblaciones distantes podría haber aumentado el flujo génico. Puede pensarse entonces que las sociedades agroalfareras tempranas poseían una baja densidad demográfica, una mayor probabilidad de estar bajo los efectos de la deriva génica y menores niveles de flujo génico entre regiones distantes, lo cual se vería reflejado en los bajos niveles de diversidad genética y altos de diferenciación hallados en relación con los grupos más tardíos (Figura 4).

Otros estudios han asociado el factor temporal a la composición genética de los grupos prehispánicos. Por ejemplo, Fehren-Schmitz et al. (2014) registraron cambios en el flujo génico mantenido a lo largo del tiempo entre regiones cercanas debidos a fluctuaciones climáticas que afectarían la disponibilidad de recursos. Además, la estructura genética de algunos grupos andinos se habría modificado luego de la conformación del Imperio inca (Russo et al. 2018).

Sin embargo, existen excepciones a este patrón general. Por ejemplo, es necesaria otra explicación para la menor variabilidad hallada en el sector norte de la QH, con ocupaciones durante el PDR. En este caso, podría pensarse que la ubicación del asentamiento Los Amarillos (de donde proviene la mayoría de los individuos de este grupo), a 6 km aproximadamente de la quebrada troncal (Figura S2), pudo haber provocado una disminución en su capacidad de mantener relaciones interpoblacionales, principalmente en comparación con los grupos localizados a lo largo de la principal ruta de acceso y comercio de la región. En este sentido, sería importante realizar estudios moleculares en otros sitios ubicados en las quebradas laterales (Figura S2).

Por otro lado, respecto del factor geográfico, es importante destacar el rol de las interacciones interregionales en relación con los procesos migratorios. En líneas generales, los resultados hallados en este trabajo son consistentes con la existencia de flujo génico entre las poblaciones del NOA analizadas, puesto que no se encontró diferenciación entre la mayoría de ellas (Figura 3). Asimismo, no se hallaron evidencias de diferenciación entre el NOA y los actuales territorios de Bolivia y norte de Chile (Figura 3). Estos resultados indican una posible vinculación genética entre estas poblaciones que puede deberse a diversos motivos no necesariamente excluyentes. Es posible que los grupos prehispánicos de los actuales territorios del NOA, norte de Chile y Bolivia (zona andina) hayan tenido un origen poblacional común, lo cual se vería reflejado en la existencia de un acervo génico compartido y, por lo tanto, en la falta de diferenciación genética entre ellos. Sin embargo, el tiempo transcurrido desde el poblamiento inicial de estas regiones a partir de un sustrato común pudo haber generado diferencias. En efecto, en este trabajo se observó, en general, una diferenciación entre las poblaciones andinas tempranas y tardías, que coincide con otras evidencias de cambios temporales mencionadas con anterioridad. Teniendo en cuenta esto, es probable que las relaciones que mantuvieron las poblaciones prehispánicas del actual NOA, norte de Chile y Bolivia, conocidas desde la arqueología gracias a los mecanismos de intercambio de bienes (Nielsen 2003), hayan involucrado también flujo génico.

A pesar de que para estas poblaciones de regiones cercanas, como los territorios del NOA, norte de Chile y Bolivia, no se encontró diferenciación genética, no fue posible definir un patrón claro de distancia genética que se asociara con la distribución geográfica de los grupos a lo largo del Área Centrosur Andina. Es decir que no se hallaron evidencias de correlación entre la diferenciación genética y la distancia geográfica entre los grupos estudiados. Esto puede deberse a que algunas poblaciones muy distanciadas no presentaron diferencias (por ejemplo, QH Sur y varios grupos de Perú), mientras que otras muy próximas, sí (como Puca y Tompullo 2, los sitios más cercanos del conjunto analizado, que se encuentran a escasos 19 km) (Figuras 1 y 3). Idealmente, se espera que, a medida que aumenta la distancia geográfica, se incremente también la diferenciación genética entre las poblaciones. En el marco de la genética de poblaciones, este efecto es conocido como aislamiento por distancia (Wright 1943). El hecho de que los grupos humanos analizados en este trabajo se alejen de este modelo podría explicarse por diversos factores.

Por ejemplo, podría pensarse que algunos asentamientos prehispánicos diferían en cuanto a su variabilidad genética debido a niveles de endogamia variables, relacionados con pautas socioculturales, políticas y/o económicas que habrían determinado su grado de vinculación con otros. Por ejemplo, de acuerdo con Baca et al. (2014), el sitio Tompullo 2 (siglos XV y XVI) habría sido habitado por un pequeño grupo de pastores posiblemente trasladados a esa región como mitmaqkunas o mitimaes, práctica habitual durante el Imperio incaico (e.g., Williams 2000). En este caso, la escasa variabilidad observada se explicaría como una consecuencia de que dicho grupo permaneció relativamente aislado y durante un breve período en ese lugar, lo que habría dificultado el contacto con poblaciones vecinas.

En el extremo opuesto, la mayor variabilidad y el menor grado de diferenciación con otras poblaciones pueden explicarse por procesos vinculados también a los patrones demográficos y la organización social de cada asentamiento. Por un lado, los grupos Valle Calchaquí y QH Sur incluyeron sitios del PDR, como algunos pucarás y El Tero (Tabla 1), donde se habría concentrado gran cantidad de población. En efecto, los lugares con alta densidad demográfica habrían aumentado durante el PDR como consecuencia de un abandono de pequeñas comunidades y aglutinamiento alrededor de grandes asentamientos (Nielsen 2001). Este proceso de concentración habría implicado un aumento del tamaño poblacional en cada sitio, lo que habría incrementado, por consiguiente, la diversidad y variabilidad genética, sobre todo suponiendo que las poblaciones del período anterior eran más pequeñas, más homogéneas internamente y poseían menor flujo génico entre grupos, incluso dentro de una misma región.

Por otro lado, mayores niveles de variabilidad y menor diferenciación también se hallaron en sitios correspondientes a las capitales de dos Estados consolidados como fueron Tiwanaku (1450-800 años AP) y el Imperio incaico (520-414 años AP). Estos asentamientos, además de poseer gran densidad poblacional, constituían grandes centros receptores de migrantes que habrían aportado diversidad al acervo génico.

En el caso del Estado Tiwanaku (1450-800 años AP), se establecieron colonias compuestas por individuos de diversas regiones, con lo que se conformaron asentamientos con diferentes grupos sociales y étnicos (Janusek 2004). A su vez, la capital imperial constituía un centro ceremonial y de peregrinación desde diferentes regiones controladas por Tiwanaku (Smith y Schreiber 2006). La alta variabilidad genética y la baja diferenciación registradas para este sitio (Figuras 2, 3 y 4) se explicarían entonces como consecuencia de la existencia de espacios donde convivían distintas comunidades.

Por su parte, las prácticas imperiales incaicas también implicaron un aumento en la movilidad de las poblaciones y, con esto, una reconfiguración de las relaciones sociales. Dentro del Tawantinsuyu, el traslado de poblaciones formó parte del proceso de conquista e implicó la relocalización de grupos en otras regiones para cumplir con diversas tareas bajo los sistemas de mitmaqkunas, yanakunas (Williams 2000; Moseley 2008) y mitas (Acuto 2012). Además, es importante considerar que, bajo este Imperio, algunos asentamientos constituyeron centros ceremoniales. Por ejemplo, Acchaymarca (ocupado entre los siglos XI y XVI) (Baca et al. 2014) habría sido un lugar de peregrinaje desde diferentes regiones, y habría representado un espacio multiétnico. Consistente con esto, este sitio fue uno de los menos diferenciados del resto (Figuras 3 y 4), y su diversidad haplotípica (o de linajes) se encontró entre los valores más altos (Figura 2 y Tabla S6), lo cual confirma hallazgos anteriores (Baca et al. 2014).

En síntesis, el auge y la caída de Estados e Imperios ocurridos durante el último milenio en la región surandina dio lugar a intensos movimientos poblacionales que habrían afectado la composición genética de diversos grupos. Asimismo, el hecho de que los grupos humanos analizados en este trabajo se alejen del modelo esperado de aislamiento por distancia refuerza la idea de que los patrones de flujo génico y diferenciación habrían estado también bajo la influencia de otro tipo de factores. En particular, sumado a la temporalidad, el grado de aislamiento o “apertura” (recepción de migrantes) de un determinado grupo o comunidad habría estado altamente influenciado por componentes socioculturales, políticos y económicos.

CONCLUSIONES

En las poblaciones prehispánicas andinas se han registrado interacciones y procesos migratorios entre diferentes comunidades a partir de distintas líneas de evidencias, como arqueológicas (e.g., Nielsen 2001, 2003; Janusek 2004; Yacobaccio 2012), bioarqueológicas (e.g., Blom et al. 2004) y métodos químicos (e.g., Tung y Knudson 2011). En el presente trabajo se evaluó la evidencia genética como otra fuente valiosa de información para aportar a la discusión.

Al analizar secuencias de ADNmt en poblaciones del Área Centrosur Andina correspondientes a distintas temporalidades, se halló que diversos aspectos socioculturales habrían influenciado la estructuración genética, la distribución de la variabilidad y las relaciones que se establecieron entre grupos.

En primer lugar, la temporalidad sería un factor importante a tener en cuenta al analizar la diferenciación genética entre poblaciones prehispánicas. En este trabajo se registraron diferencias entre los grupos más antiguos, que además poseían menos variabilidad, y los de periodos posteriores.

Por otro lado, al estudiar grupos humanos, la distancia geográfica no es el único factor que puede determinar las interacciones y el flujo génico. En efecto, se observó que los grupos más alejados geográficamente no fueron necesariamente los más diferenciados, y que algunas comunidades muy cercanas no habrían mantenido flujo génico. Esto concuerda con una prevalencia de factores socioculturales, políticos y económicos en el establecimiento de vínculos entre las poblaciones humanas prehispánicas. Tanto el flujo génico como otros tipos de relaciones entre grupos andinos de diferentes regiones pudieron estar determinados por la disponibilidad de recursos, por ejemplo, y ser relativamente independientes de la distancia geográfica.

Además, la diversidad genética también se halló relacionada con características sociopolíticas. Por ejemplo, aquellos sitios que fueron centros ceremoniales o capitales de Estados, con alta influencia de migrantes de distintas regiones, presentaron mayores niveles de variabilidad. En contraposición, una menor diversidad genética se encontró en asentamientos de menor envergadura y relativamente más aislados o con menor flujo génico con otras poblaciones.

Finalmente, se destaca la importancia de incorporar a estos estudios otro tipo de marcadores moleculares, como los autosómicos. Su análisis permitiría, entre otras cosas, obtener información acerca de las relaciones de parentesco genético entre individuos (e.g., Russo et al. 2016; Mendisco et al. 2018). Evaluar la importancia del parentesco en la organización social andina resulta relevante porque, por ejemplo, las pautas que definen las uniones entre individuos pueden influenciar la diversidad genética poblacional (Guillot et al. 2015). Además, los marcadores autosómicos permiten analizar los patrones de migración diferencial entre hombres y mujeres (e.g., Baca et al. 2014), relacionados con los sistemas de residencia y la movilidad. En tal sentido, los datos genéticos podrían profundizar los conocimientos sobre fenómenos como el intercambio y las situaciones de conflicto, y su relación con los mecanismos de flujo génico.

Agradecimientos

Los autores deseamos agradecer a los revisores anónimos, cuyas sugerencias ayudaron a mejorar sustancialmente el manuscrito. El presente trabajo se llevó a cabo gracias a una beca doctoral cofinanciada por CONICET y Universidad Maimónides. Agradecemos además el apoyo financiero de ANPCyT (PICT 2014-3012), UBACYT (20020150200233BA-2016), Fundación Científica Felipe Fiorellino y Fundación de Historia Natural Félix de Azara.

REFERENCIAS

Acuto, F. (2012). Landscapes of Inequality, Spectacle and Control: Inka Social Order in Provincial Contexts. Revista Chilena de Antropología 25: 9-64. doi: 10.5354/0719-1472.2012.20256

Aschero, C. (1984). El sitio ICc-4: Un asentamiento precerámico en la Quebrada de Inca Cueva (Jujuy. Argentina). Estudios Atacameños 7: 62-72. doi: 10.22199/S07181043.1984.0007.00006

Baca, M., K. Doan, M. Sobczyk, A. Stankovic y P. Węgleński. (2012). Ancient DNA reveals kinship burial patterns of a pre-Columbian Andean community. BMC Genetics 13 (1): 30. doi: 10.1186/1471-2156-13-30

Baca, M., M. Molak, M. Sobczyk, P. Węgleński y A. Stankovic. (2014). Locals, resettlers, and pilgrims: a genetic portrait of three pre-Columbian Andean populations. American Journal of Physical Anthropology 154 (3): 402-412. doi: 10.1002/ajpa.22524

Blom, D. E. (2005). Embodying borders: human body modification and diversity in Tiwanaku society. Journal of Anthropological Archaeology 24 (1): 1-24. doi: 10.1016/j.jaa.2004.10.001

Blom, D. E., M. C. Lozada Cerna y B. Hallgrimson. (2004). Diversidad interna en Tiwanaku y Chiribaya: estructuras sociopolíticas serranas y costeras. En Tiwanaku: Aproximaciones a sus contextos históricos y sociales, compilado por M. A. Rivera y A. L. Kolata, pp. 165-177. Editorial Universidad Bolivariana, Santiago de Chile.

Carnese, F., F. Mendisco, C. Keyser, C. B. Dejean, J. M. Dugoujon, C. M. Bravi, B. Ludes y E. Crubézy. (2010). Paleogenetical study of pre-Columbian samples from Pampa Grande (Salta, Argentina). American Journal of Physical Anthropology 141 (3): 452-462. doi: 10.1002/ajpa.21165

Di Rienzo, J.A., F. Casanoves, M. G. Balzarini, L. González, M. Tablada y C. W. Robledo. (2015). InfoStat versión 2015. Grupo InfoStat, FCA, Universidad Nacional de Córdoba, Argentina. http://www.infostat.com.ar

Excoffier, L. y H. E. L. Lischer. (2010). Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources 10 (3): 564-567. doi: 10.1111/j.1755-0998.2010.02847.x

Fehren-Schmitz, L., M. Reindel, E. Tomasto Cagigao, S. Hummel y B. Herrmann. (2010). Pre-Columbian population dynamics in coastal southern Peru: A diachronic investigation of mtDNA patterns in the Palpa region by ancient DNA analysis. American Journal of Physical Anthropology 141 (2): 208-221. doi: 10.1002/ajpa.21135

Fehren-Schmitz, L., O. Warnberg, M. Reindel, V. Seidenberg, E. Tomasto-Cagigao, J. Isla-Cuadrado, S. Hummel y B. Herrmann. (2011). Diachronic investigations of mitochondrial and Y-chromosomal genetic markers in pre-Columbian Andean highlanders from south Peru. Annals of Human Genetics 75 (2): 266-283. doi: 10.1111/j.1469-1809.2010.00620.x

Fehren-Schmitz, L., W. Haak, B. Mächtle, F. Masch, B. Llamas, E. Tomasto Cagigao, V. Sossna, K. Schittek, J. Isla Cuadrado, B. Eitel y M. Reindel. (2014). Climate change underlies global demographic, genetic, and cultural transitions in pre-Columbian southern Peru. Proceedings of the National Academy of Sciences USA 201403466. doi: 10.1073/pnas.1403466111

Fehren-Schmitz, L., B. Llamas, S. Lindauer, E. Tomasto-Cagigao, S. Kuzminsky, N. Rohland, F. R. Santos, P. Kaulicke, G. Valverde, S. M. Richards, S. Nordenfelt, V. Seidenberg, S. Mallick, A. Cooper, D. Reich y W. Haak. (2015). A Re-Appraisal of the Early Andean Human Remains from Lauricocha in Peru. PLoS One 10 (6): e0127141. doi: 10.1371/journal.pone.0127141

Fernández Distel, A. (1986). Las cuevas de Huachichocana, su posición dentro del precerámico con agricultura incipiente del Noroeste Argentino. Beiträge zur Allgemeinen und Vergleichenden Archäologie 8: 353-430.

Goldstein, P. S. y B. Owen. (2001). Tiwanaku en Moquegua: las colonias altiplánicas. Boletín de Arqueología, Pontificia Universidad Católica del Perú 5: 139-168.

Guillot, E. G., M. L. Hazelton, T. M. Karafet, J. S. Lansing, H. Sudoyo y M. P. Cox. (2015). Relaxed observance of traditional marriage rules allows social connectivity without loss of genetic diversity. Molecular Biology and Evolution 32 (9): 2254-2262. doi: 10.1093/molbev/msv102

Hijmans, R. J. (2015). geosphere: Spherical Trigonometry. R package version 1.3-13. http://CRAN.R-project.org/package=geosphere

Janusek, J. W. (2004). Tiwanaku and Its Precursors: Recent Research and Emerging Perspectives. Journal of Archaeological Research 12 (2): 121-183. doi: 10.1023/B:JARE.0000023711.96664.1b

Jobling, M. A. (2012). The impact of recent events on human genetic diversity. Philosophical Transactions of the Royal Society of London B Biological Sciences 367 (1590): 793-799. doi: 10.1098/rstb.2011.0297

Kemp, B. M., T. A. Tung y M. L. Summar. (2009). Genetic continuity after the collapse of the Wari Empire: mitochondrial DNA profiles from Wari and post-Wari populations in the ancient Andes. American Journal Physical Anthropology 140 (1): 80-91. doi: 10.1002/ajpa.21037

Llagostera, A. (2010). Retomando los límites y las limitaciones del “archipiélago vertical”. Chungara (Arica) 42 (1): 283-295. doi: 10.4067/S0717-73562010000100036

Librado, P. y J. Rozas. (2009). DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451-1452. doi: 10.1093/bioinformatics/btp187

Llamas, B., L. Fehren-Schmitz, G. Valverde, J. Soubrier, S. Mallick, N. Rohland, S. Nordenfelt, C. Valdiosera, S. M. Richards, A. Rohrlach, M. I. Romero, I. F. Espinoza, E. T. Cagigao, L. W. Jiménez, K. Makowski, I. S. Reyna, J. M. Lory, J. A. Torrez, M. A. Rivera, R. L. Burger, M. C. Ceruti, J. Reinhard, R. S. Wells, G. Politis, C. M. Santoro, V. G. Standen, C. Smith, D. Reich, S. Y. Ho, A. Cooper y W. Haak. (2016). Ancient mitochondrial DNA provides high-resolution time scale of the peopling of the Americas. Science Advances 2 (4): e1501385. doi: 10.1126/sciadv.1501385

Luciani, S., G. Fornaciari, O. Rickards, C. Martínez Labraga y F. Rollo. (2006). Molecular characterization of a pre-columbian mummy and in situ coprolite. American Journal of Physical Anthropology 129 (4): 620-629. doi: 10.1002/ajpa.20314

Mantel, N. A. (1967). The detection of disease clustering and a generalized regression approach. Cancer Research 27 (2 part 1): 209-220.

Manzanilla, L. (2007). Comentarios finales: el estudio de la complejidad, las identidades y las interacciones. Una propuesta metodológica. En Sociedades precolombinas surandinas: temporalidad, interacción y dinámica cultural del NOA en el ámbito de los Andes Centro-Sur, editado por V. I. Williams, B. N. Ventura, A. B. M. Callegari y H. D. Yacobaccio, pp. 393-418. Taller Internacional de Arqueología del NOA y Andes Centro Sur, Buenos Aires.

Martel, A. (2014). Aguas Calientes. Evidencias directas de tráfico caravanero entre la Puna meridional y el valle Calchaquí. Estudios Sociales del NOA: nueva serie 13: 103-124.

Mendisco, F., C. Keyser, C. Hollard, V. Seldes, A. Nielsen, E. Crubézy y B. Ludes. (2011). Application of the iPLEXTM Gold SNP genotyping method to the analysis of Amerindian ancient DNA samples: benefits for ancient population-based studies. Electrophoresis 32 (3-4): 386-393. doi: 10.1002/elps.201000483

Mendisco, F., C. Keyser, V. Seldes, C. Rivolta, P. Mercolli, P. Cruz, A. E. Nielsen, E. Crubezy y B. Ludes. (2014). Genetic diversity of a late prehispanic group of the Quebrada de Humahuaca, Northwestern Argentina. Annals of Human Genetics 78 (5): 367-380. doi: 10.1111/ahg.12075

Mendisco, F., C. Keyser, V. Seldes, A. E. Nielsen, M. G. Russo, E. Crubézy y B. Ludes. (2018). An insight into the burial practices of the late pre-Hispanic Los Amarillos community (northwestern Argentina) through the study of ancient DNA. Journal of Archaeological Science 91: 12-19. doi: 10.1016/j.jas.2018.01.005

Moseley, M. (2008). The Incas and their Ancestors: The Archaeology of Peru. Revised Edition. Thames & Hudson, Londres.

Nei, M. (1987). Molecular evolutionary genetics. Columbia University Press, Nueva York.

Nei, M. y W. H. Li. (1979). Mathematical model for studying genetic variation in terms of restriction endonucleases. Proceedings of the National Academy of Sciences USA 76 (10): 5269-5273. doi: 10.1073/pnas.76.10.5269

Nielsen, A. E. (2000). Andean Caravans: an ethnoarchaeology. Tesis Doctoral inédita. Department of Anthropology, The University of Arizona, Tucson.

Nielsen, A. E. (2001). Evolución social en Quebrada de Humahuaca (AD 700-1536). En Historia argentina prehispánica, t. 1, editado por E. Berberián y A. E. Nielsen, pp. 171-264. Brujas, Córdoba.

Nielsen, A. E. (2003). Por las rutas del Zenta: evidencias directas del tráfico prehispánico entre Humahuaca y las Yungas. En La mitad verde del mundo andino. Investigaciones arqueológicas en la vertiente oriental de los Andes y las tierras bajas de Bolivia y Argentina, editado por B. Ventura y G. Ortiz, pp. 261-284. EdiUnju, Jujuy.

Olivera, D. (2001). Sociedades agropastoriles tempranas: el Formativo Inferior del Noroeste Argentino. En Historia Argentina Prehispánica, t. 1, editado por E. Berberian y A. Nielsen, pp. 83-127. Brujas, Córdoba.

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

Raghavan, M., M. Steinrücken, K. Harris, S. Schiffels, S. Rasmussen, M. DeGiorgio, A. Albrechtsen, C. Valdiosera, M. C. Ávila-Arcos, A. S. Malaspinas, A. Eriksson, I. Moltke, M. Metspalu, J. R. Homburger, J. Wall, O. E. Cornejo, J. V. Moreno-Mayar, T. S. Korneliussen, T. Pierre, M. Rasmussen, P. F. Campos, P. de Barros Damgaard, M. E. Allentoft, J. Lindo, E. Metspalu, R. Rodríguez-Varela, J. Mansilla, C. Henrickson, A. Seguin-Orlando, H. Malmström, T. Stafford Jr., S. S. Shringarpure, A. Moreno-Estrada, M. Karmin, K. Tambets, A. Bergström, Y. Xue, V. Warmuth, A. D. Friend, J. Singarayer, P. Valdes, F. Balloux, I. Leboreiro, J. L. Vera, H. Rangel-Villalobos, D. Pettener, D. Luiselli, L. G. Davis, E. Heyer, C. P. E. Zollikofer, M. S. Ponce de León, C. I. Smith, V. Grimes, K. A. Pike, M. Deal, B. T. Fuller, B. Arriaza, V. Standen, M. F. Luz, F. Ricaut, N. Guidon, L. Osipova, M. I. Voevoda, O. L. Posukh, O. Balanovsky, M. Lavryashina, Y. Bogunov, E. Khusnutdinova, M. Gubina, E. Balanovska, S. Fedorova, S. Litvinov, B. Malyarchuk, M. Derenko, M. J. Mosher, D. Archer, J. Cybulski, B. Petzelt, J. Mitchell, R. Worl, P. J. Norman, P. Parham, B. M. Kemp, T. Kivisild, C. Tyler-Smith, M. S. Sandhu, M. Crawford, R. Villems, D. Glenn Smith, M. R. Waters, T. Goebel, J. R. Johnson, R. S. Malhi, M. Jakobsson, D. J. Meltzer, A. Manica, R. Durbin, C. D. Bustamante, Y. S. Song, R. Nielsen y E. Willerslev. (2015). Genomic evidence for the Pleistocene and recent population history of Native Americans. Science 349 (6250): aab3884. doi: 10.1126/science.aab3884

Russo, M. G. (2016). Análisis genéticos aplicados al estudio de la dinámica de poblaciones humanas prehispánicas del Noroeste Argentino. Tesis Doctoral inédita. Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires.

Russo, M. G., F. Mendisco, S. Avena, C. B. Dejean y V. Seldes. (2016). Pre-Hispanic Mortuary Practices in Quebrada de Humahuaca (North-western Argentina): Genetic Relatedness among Individuals Buried in the Same Grave. Annals of Human Genetics 80 (4): 210-220. doi: 10.1111/ahg.12159

Russo, M. G., M. S. Gheggi, S. A. Avena, C. B. Dejean y M. B. Cremonte. (2017). Linajes mitocondriales en muestras de Esquina de Huajra (Jujuy, Argentina). Aportes al estudio de la ocupación Incaica en la región y la procedencia de sus habitantes. Revista Argentina de Antropología Biológica 19 (1): 15. doi: 10.17139/raab.19.1.15

Russo, M. G., F. Mendisco, S. A. Avena, C. M. Crespo, V. Arencibia, C. B. Dejean y V. Seldes. (2018). Ancient DNA reveals temporal population structure within the Southcentral Andes Area. American Journal of Physical Anthropology 166 (4): 851-860. doi: 10.1002/ajpa.23475

Schwarz, G. (1978). Estimating the dimension of a model. The Annals of Statistics 6 (2): 461-464. doi: 10.1214/aos/1176344136

Smith, M. E. y K. J. Schreiber. (2006). NewWorld States and Empires: Politics, Religion, and Urbanism. Journal of Archaeological Research 14 (1): 1-52. doi: 10.1007/s10814-005-9000-1

Tamura, K. y M. Nei. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular Biology and Evolution 10 (3): 512-526. doi: 10.1093/oxfordjournals.molbev.a040023

Tamura, K., G. Stecher, D. Peterson, A. Filipski y S. Kumar. (2013). MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Molecular Biology and Evolution 30 (12): 2725-2729. doi: 10.1093/molbev/mst197

Tung, T. A. y K. J. Knudson. (2008). Social identities and geographical origins of Wari trophy heads from Conchopata, Peru. Current Anthropology 49 (5): 915-925. doi: 10.1086/591318

Tung, T. A. y K. J. Knudson. (2011). Identifying locals, migrants, and captives in the Wari Heartland: A bioarchaeological and biogeochemical study of human remains from Conchopata, Peru. Journal of Anthropological Archaeology 30 (3): 247-261. doi: 10.1016/j.jaa.2011.06.005

Varela, H. H., J. A. Cocilovo, M. L. Fuchs y T. G. O’Brien. (2014). Coastal, valley, and oasis interaction: impact on the evolution of ancient populations in the South Central Andes. American Journal of Physical Anthropology 155 (4): 591-599. doi: 10.1002/ajpa.22614

Vincenty, T. (1975). Direct and inverse solutions of geodesics on the ellipsoid with application of nested equations. Survey Review 23 (176): 88-93. doi: 10.1179/sre.1975.23.176.88

Williams, V. (2000). El imperio Inka en la provincia de Catamarca. Intersecciones en Antropología 1 (1): 55-78.

Wright, S. (1943). Isolation by distance. Genetics 28 (2): 114-138.

Yacobaccio, H. D. (2012). Intercambio y caravanas de llamas en el sur andino (3000-1000 AP). Comechingonia 16 (1): 31-51.

Notas