DATOS GENÉTICOS E ISOTÓPICOS DE INDIVIDUOS VINCULADOS A LA REDUCCIÓN FRANCISCANA DE BARADERO (SIGLO XVII)

GENETIC AND ISOTOPIC PROFILE OF INDIVIDUALS LINKED TO BARADERO`S FRANCISCAN REDUCTION (XVII CENTURY)

DADOS GENÉTICOS E ISOTÓPICOS DE INDIVÍDUOS LIGADOS À REDUÇÃO FRANCISCA DE BARADERO (SÉC. XVII)

INTRODUCCIÓN

Durante la conquista de América los europeos fueron ocupando territorios mediante la instalación de diversos asentamientos con poblaciones estables, ya sea con fines militares, urbanísticos, religiosos o mercantiles. El propósito consistía en colonizar espacios donde se esperaba hallar nuevas riquezas y, por lo tanto, evitar que fueran susceptibles de ser invadidos por otros países conquistadores. En las colonias establecidas por la Corona Española, los asentamientos que tuvieron un marcado impacto sobre la población indígena fueron las encomiendas y las reducciones o pueblos de indios. Ambos tipos de instalaciones fueron parte de misma empresa colonial, que buscaba controlar la fuerza de trabajo indígena y transformar sus prácticas sociales, económicas, políticas y religiosas (Assadourian 1983; Stern 1986; Bethell 1990). A partir de estas estrategias de reducción de la población indígena en espacios controlados se intentó crear sirvientes confiables, que pudiesen imitar en alguna medida las prácticas culturales (lenguaje, vestimenta, dieta y ética laboral) de la clase dominante. El sistema reduccional facilitó el adoctrinamiento y la defensa contra los indígenas hostiles o de otros agentes antagonistas de las colonias; también permitió a los misioneros tener bajo control a los nuevos indígenas que nacían en los pueblos de indios y a aquellos otros que iban ingresando en diferentes momentos y circunstancias. Las prácticas específicas de cada misión variaron en el tiempo y el espacio, pero básicamente implicaban largas jornadas de prácticas agrícolas y otras actividades destinadas a sostener las misiones y el trabajo en las encomiendas, bajo un régimen de trabajo comunal forzado, y actividades religiosas obligatorias (Lozano Pedro 1873; Millé 1961; Palomeque 2000; Lightfoot 2005; Stein 2005).

Diversas fuentes documentales de fines del siglo XVI y del siglo XVII proporcionan información sobre las estrategias de dominación implementadas durante la colonización hispánica en el Río de La Plata (Archivo General de Indias- Archivo de Charcas [en adelante AGI-ACH] 1609, 1610, 1613, 1616 a 1622, 1626, 1628, 1663, 1664, 1673, 1678, 1690, 1714, 1718, 1748, 1749; AB 1679, 1691, 1694; Schmidel 1903; Garay [1582] 1969; Díaz de Guzmán [1612] 2012). En ellas se describen aspectos vinculados a la relocalización de los pueblos indígenas que estaban asentados en las costas, los riachos del Paraná Inferior y el noreste de la llanura bonaerense. Provenientes de esos espacios, algunos representantes de los grupos querandí, guaraní y chaná fueron trasladados hasta las primeras encomiendas y pueblos de indios, que se establecieron en el área, con el objetivo de controlar los conflictos interétnicos, aprovechar su fuerza de trabajo y evangelizarlos.

Los acontecimientos ocurridos en esas primeras instalaciones hispánicas fueron narrados desde la perspectiva hegemónica de los funcionarios coloniales, especialmente para informar a las autoridades españolas sobre la situación que atravesaban los gobiernos que se sucedieron en el Río de la Plata durante el siglo XVII; a veces sobredimensionando los éxitos y disminuyendo el número de fracasos o a la inversa, según el arbitrio de cada funcionario o superioridad eclesiástica y según las circunstancias sociohistóricas particulares. No obstante, ya sea de manera implícita o explícita, en las fuentes documentales también se da cuenta de las consecuencias que tuvieron las políticas implementadas por la Corona española para asegurar el funcionamiento de las encomiendas y las reducciones indígenas. En especial, a pesar del sesgo ideológico del relato histórico, se evidencian sucesos que condujeron al desarrollo de procesos de desestructuración étnica, de etnogénesis y de resistencia cultural en las primeras reducciones franciscanas organizadas en el Río de la Plata a comienzos del siglo XVII (AGI-ACH 1620, 1626, 1628, 1663, 1678, 1679; Buscaglia 2011; Oliva y Lisboa 2009; Tapia et al. 2015).

Además de la documentación escrita, la complejidad de dichos procesos puede abordarse a través de los conocimientos que proporciona el registro arqueológico. Al respecto, para evaluar el impacto de la conquista en el área se han venido desarrollando varios proyectos de investigación planteados desde una perspectiva interdisciplinaria, que vincula los campos de la Arqueología Histórica, la Antropología Biológica y la Antropología Genética. Desde el punto de vista teórico-metodológico las investigaciones incluyen la interrelación entre los datos generados desde estos tres enfoques disciplinarios y los que proporcionan las fuentes documentales. Como caso paradigmático de estudio, las investigaciones se basan en la información histórica y el registro arqueológico vinculado con la reducción franciscana Santiago del Baradero, que Hernandarias fundó en 1615, entre otras reducciones instaladas contemporáneamente en el Río de la Plata (Tubichaminí, San José o de El Bagual) y otros sitios de la costa del Paraná Inferior y Medio (San Bartolomé de los Chaná, San Miguel de los Calchines Boca del Monje y San Lorenzo de los Mocoretá). La particularidad de esa reducción reside en que fue la única que prosperó desde su fundación hasta mediados del siglo XVIII (AGI-ACH 1616 a 1622, 1626, 1628, 1714, 1690, 1691, 1694, 1708, 1714, 1749; Tapia 2002, 2014).

En Santiago del Baradero fueron reducidos 250 indígenas pertenecientes a los grupos chana y guaraní, que si bien habitaban en espacios similares de costas fluviales e islas del Paraná Inferior, habían practicado formas de subsistencia tradicional diferenciadas: cazarecolección-pesca para el primer grupo y caza-recolección-horticultura para el último (Schmidel 1903; Lothrop 1932; Ramírez [1528] 1939; Lopes de Souza [1531] 1969; Díaz de Guzmán [1612] 2012). Con el transcurso del tiempo a la población inicial se fueron sumando individuos provenientes de otros grupos indígenas y de otras regiones del país. Estos nuevos aportes poblacionales fueron registrados en un padrón de indios de Santiago del Baradero, realizado entre 1684 y 1686 por orden del gobernador Herrera y Sotomayor: “mando se haga padrón de todas las dichas familias de indios con distinción de sus naturalezas, sexos y edades poniendo con separación los que han sido descendientes de la misma reducción y los que son agregados expresando las familias de los unos y los otros” (AGI-ACH, 1690, Leg. 282:4). En el padrón se registraron 99 indígenas de los cuales 9 varones adultos eran foráneos que procedían de otras regiones del nordeste y del noroeste, y residían allí porque estaban casados con mujeres nacidas en la reducción.

Si se consideran los 250 indígenas chaná y guaraníes iniciales -con los cuales Hernandarias organizó la reducción en 1615- y se compara con los 99 registrados en el padrón, se observa que la población de Santiago del Baradero después de aproximadamente 70 años tenía 151 individuos menos. Por otra parte, los indígenas que se empadronaron hacia fines del siglo XVII, ya no pueden asignarse a etnías específicas dado que -como también refieren los documentos de años anteriores-, la dispersión y relocalización de indios produjo procesos de mestización muy marcados. Estas estrategias de dominación colonial también generaron reagrupaciones sociales multiétnicas o re-etnificaciones, en términos de Boccara (2003), que respondieron a objetivos de supervivencia, lucha y resistencia para hacer frente al colonialismo hispánico (Tapia 2014).

El registro arqueológico utilizado como evidencia del funcionamiento de la reducción Santiago del Baradero procede del sitio Cementerio Indígena, Baradero, Provincia de Buenos Aires (Figura 1). Este sitio fue excavado por Salvador Debenedetti en 1910, oportunidad en la que realizó el hallazgo de catorce enterratorios asociados con materiales tanto de procedencia europea como indígena (Debenedetti 1911). Teniendo en cuenta los materiales de metal y las cuentas de vidrio consideró que la formación del sitio debía ubicarse en tiempos poshispánicos y que, según los documentos escritos, podría vincularse con dicha reducción.

Las investigaciones arqueológicas actuales han permitido corroborar la atribución cronológica de Debenedetti a partir de cuatro fechados radiocarbónicos: dos fechados de C14 convencional efectuados en el LATYR (cal 1578 - 1621 DC 1σ y cal 1646 - 1697 DC 1σ) y dos fechados AMS obtenidos en el Laboratorio de la Universidad de Georgia (cal 1598 - 1609 DC 1σ y cal 1594 – 1613 DC 2σ), ya presentados en publicaciones anteriores (Kozameh et al. 2015; Rizzo et al. 2015).

Asimismo, a partir de los materiales arqueológicos recuperados en ese sitio -que actualmente integran una colección depositada en el Museo Histórico A. Barbich de Baradero, Pcia de Buenos Aires-, se han efectuado estudios morfológico-funcionales y arqueométricos de los materiales de metal y de las cuentas de vidrio (Tapia et al. 2009; Tapia y Pineau 2011, 2013), así como también de los restos óseos humanos desde un abordaje bioantropológico (Kozameh et al. 2015). Los resultados obtenidos han permitido realizar inferencias en torno de los procesos de cambio cultural y biológico que habrían afectado a los individuos reducidos (Tapia 2014).

En esta presentación, se propone contribuir con nuevos datos bioantropológicos al conocimiento del impacto que produjo la colonización en las reducciones franciscanas en el área de estudio. A partir del análisis de los restos óseos humanos recuperados en el sitio Cementerio Indígena, se han abierto dos nuevas líneas de indagación (Cardozo et al. 2013, 2014, Rizzo et al. 2015): 1- la identificación de posibles variaciones genéticas ocasionadas por la relocalización y coexistencia de diferentes grupos indígenas en Santiago del Baradero, mediante el análisis de ADN mitocondrial (en adelante ADNmt) antiguo; y 2- la caracterización del tipo de dieta y sus cambios respecto de las poblaciones prehispánicas, a partir del análisis de isótopos estables δ13C, δ15N.

En primer lugar, bajo el encuadre teórico de la genética de poblaciones se efectúan los análisis del ADNmt el cual está relacionado con los linajes maternos. Este tipo de investigaciones realizadas tanto en restos arqueológicos como en poblaciones autóctonas contemporáneas, ha permitido definir cinco haplogrupos (en adelante Hg) amerindios fundadores de la población del continente americano, llamados A2, B2, C1 y D1 correspondientes al 95% de los Hg amerindios y X2a solo descripto para Norteamérica (Tamm et al. 2007; Achilli et al. 2008; Kemp y Schurr 2010; Perego et al. 2010; Crespo 2011). Estos Hg a su vez se encuentran subdivididos en haplotipos según mutaciones puntuales adicionales a lo largo de todo el genoma mitocondrial (Tamm et al. 2007; Achilli et al. 2008; Perego et al. 2010). De esta manera en la genética de poblaciones se puede incorporar el componente temporal al estudio de la variabilidad genética y, en consecuencia, realizar el análisis filogenético de las poblaciones pasadas y las especies extintas. A su vez, la información obtenida permite responder a interrogantes sobre los procesos migratorios, patrones de residencia y relaciones genéticas entre poblaciones humanas antiguas, reconstruyendo así la historia poblacional en tiempo y espacio (Stone y Stoneking 1993; Hermann y Hummel 1994).

En segundo lugar, dado que las prácticas dietarias constituyen una de las actividades culturales fundamentales en la vida de los seres humanos, desde el abastecimiento de las necesidades calóricas básicas hasta las variables socioeconómicas que afectan la producción de los alimentos, adquisición, preparación y consumo; es posible considerar que el proceso colonizador haya impactado en gran medida sobre la forma de alimentarse de los grupos aborígenes. En los establecimientos coloniales entre algunos de los factores que facilitaron la reducción y adoctrinamiento de los indígenas se encontraba la oferta de alimentos, en especial los hidratos de carbono basados en harina de maíz y trigo, vegetales cultivados de manera extensiva y procesados en las mismas reducciones (Larsen 1987, 2002; Tapia 2002).

Dada la importancia que tienen las prácticas alimenticias para el conocimiento de las formas de vida de las poblaciones del pasado, se han utilizado diversas fuentes de información para su estudio. Sin embargo, en muchos casos los resultados obtenidos consisten simplemente en determinar el menú principal, y solo alcanzan para proponer estimaciones semi cuantitativas de la porción animal que se incluye en la dieta, ya que el componente vegetal resulta más difícil de determinar debido a que es más susceptible a la degradación (Tykot et al. 2014). A principios de los años ´70 se intentó adquirir un conocimiento más detallado sobre el tipo de dieta que habrían tenido los individuos en el pasado. Esto se logró mediante la implementación de una técnica que permite analizar, en tejidos biológicos, la relación existente entre los átomos de Carbono (12C) y Nitrógeno (14N) con sus isótopos estables (13C y 15N respectivamente). La misma se sustenta en el hecho de que estos elementos, adquiridos a través de la alimentación, constituyen dos de los átomos principales que componen a los seres vivos (Schoeninger 1995; Tykot 2006). Es importante destacar que en la naturaleza estas variantes se encuentran presentes en una proporción determinada con respecto al átomo más abundante, por ejemplo 12C se encuentra presente en un 98,93% mientras que su isótopo 13C en un 1,07%.

Los isótopos estables son variantes atómicas de distinto peso molecular que no decaen radioactivamente. Se caracterizan por poseer un núcleo con distinto número de neutrones pero igual número de protones y electrones, a modo de ejemplo 12C posee 6 protones y 6 neutrones, mientras que 13C posee 6 protones y 7 neutrones. Esta variación en el número de neutrones no modifica sus propiedades químicas, determinadas por los electrones, pero sí su velocidad de reacción debido a los distintos pesos moleculares de los isótopos (Schoeninger 1995). Las diferencias ocasionadas en las velocidades de reacción se producen en los distintos procesos bioquímicos como la fotosíntesis y vías metabólicas, ocasionando un fraccionamiento. Este se define como la diferencia en las relaciones isotópicas entre el sustrato (la comida ingerida) y el producto (colágeno, músculo, grasa, etc.). Estas diferencias son transmitidas a los largo de las redes tróficas provocando fraccionamientos particulares que pueden ser detectados en un espectrómetro de masa. Estos procesos presentan variabilidad en función del isótopo estable que se considere, del tipo de especie elegida y del material que se utilice en un mismo individuo (Panarello et al. 2009). A modo de ejemplo, según los valores de los isótopos estables de carbono se puede distinguir si la proteína de la dieta proviene de animales marinos o terrestres o bien de algún tipo de planta, dependiendo de la vía fotosintética que estas utilicen (C3, C4 o CAM). Por su parte, en los valores de los isótopos estables de nitrógeno se puede distinguir si la proteína proviene de plantas o animales y, de esta manera, definir el nivel trófico que un organismo ocupa en la cadena alimenticia (Wilson et al. 2007).

Otra característica importante de las investigaciones realizadas en isótopos estables corresponde al recambio atómico que ocurre en los tejidos de los individuos, por ejemplo el colágeno de los huesos posee un recambio atómico lento por lo que los resultados obtenidos al respecto constituirán el promedio de la dieta de ese individuo en los últimos 7 a 10 años de su vida; por su parte las piezas dentales una vez formadas no poseen un recambio atómico (Tykot 2006). Estos datos son relevantes y deben tenerse presentes en las interpretaciones, ya que los huesos nos darán el promedio de la dieta de los últimos años de la vida de un individuo mientras que los datos obtenidos a partir de piezas dentales aportaran información de los primeros años en los cuales se formaron las piezas analizadas.

Se considera que la información recabada a partir de los análisis bioarqueológicos y arqueogenéticos pueden aportar información crítica para entender la historia adaptativa de nuestra especie, el modo de vida, la variabilidad genética de los individuos y su relación con otras poblaciones, además de permitir interpretar los acontecimientos ocurridos durante la conquista hispánica a partir de datos que provienen directamente de los restos de los individuos involucrados. Si bien este tipo de información también se encuentra expuesta a problemas de interpretación, estos son diferentes de aquellos que surgen al analizar los documentos escritos y la cultura material, debido a los sesgos particulares, característicos de ambos tipos de registro (Larsen y Walker 2005).

MATERIALES Y MÉTODOS

Tipo de muestras y tratamiento utilizado en los análisis genéticos

El tratamiento de las muestras y la extracción de ADNmt antiguo se efectuaron en el Laboratorio del Centro de Estudios Biomédicos, Biotecnológicos, Ambientales y Diagnóstico (CEBBAD) de la Universidad Maimónides y la secuenciación de las muestras con resultado positivo se realizó en la Unidad de Genómica, Nodo Plataforma Nacional de Genómica CATG del INTA sede Hurlingham, Buenos Aires. Para efectuar los análisis se seleccionaron 20 muestras que proceden del sitio poshispánico Cementerio Indígena (33º 48` 0,4`` Lat. Sur, 59º 30` 6,18`` Long. Oeste, Figura 1); 13 de ellas corresponden a restos óseos y 7 a piezas dentales (Tabla 1). Teniendo en cuenta que los restos forman parte de una colección con elementos anatómicos donde las mandíbulas y maxilares no se encuentran integradas con los restos postcraneales, no se consideró a estos últimos para el cálculo del total de individuos muestreados. Por ello, solo resulta confiable establecer un número mínimo de 7 individuos a partir de las piezas dentarias (CIB15 a CIB21), que fueron extraídas de mandíbulas o de maxilares de diferentes individuos.

Para el tratamiento de las muestras se aplicaron los procedimientos habituales según se trate de piezas dentarias o de hueso. En el caso de las muestras de molares, en primer lugar se confeccionó un molde de la pieza, dado que por el estado de conservación esta podía deteriorarse durante su tratamiento, siendo luego imposible remontarla. En el caso de los restos óseos se procedió a la limpieza de la superficie con una gasa humedecida con hipoclorito de sodio (lavandina) al 10%, a fin de eliminar la carga de ADN exógeno. Posteriormente con un taladro Dremel® provisto de un disco de corte fino se obtuvo un pequeño fragmento del hueso, preferentemente de la región más compacta y gruesa del mismo. Al fragmento óseo obtenido se le eliminaron unos milímetros de la capa superficial utilizando el taladro Dremel® provisto de una fresa de pulido de alta velocidad.

Tanto el resto óseo limpio como la pieza dental se sumergieron en lavandina al 10% y se los enjuagó dos veces con agua para biología molecular, bidestilada y esteril. Finalmente se irradió la muestra con luz ultravioleta (UV) durante 15 minutos de cada lado. Una vez eliminada la contaminación superficial se procedió a la obtención de polvo de hueso o diente para permitir una mayor superficie de contacto entre los reactivos y la muestra. En el caso de los fragmentos óseos y de las piezas dentales pequeñas o en mal estado de conservación se procedió a la molienda dentro de un molino Ika yellow line A10. Por su parte, para el resto de las piezas dentales se seccionaron transversalmente desbastando su interior con un minitorno Dremel®, con el objetivo de obtener dentina a partir de la pulpa; finalmente la pieza fue reconstituida.

Una alícuota del polvo obtenido se trató con distintos reactivos, dependiendo si se trataba de extraer ADN mediante equipos comerciales o de protocolos diseñados especialmente para muestras antiguas, ya que estos últimos permiten trabajar con mayor cantidad de muestra y posibilitan la recuperación de material genético a partir de muestras difíciles. En el primer caso, utilizamos distintos equipos como Gene Clean for Ancient DNA (MPBiomedicals) e Investigator (Qiagen), en todos los casos siguiendo las especificaciones del fabricante. En el caso de protocolos diseñados para ADN antiguo, se utilizó la extracción orgánica por fenol-cloroformo (Carnesse et al., 2010) seguido de una purificación utilizando el kit comercial Bioneer.

Reacción en Cadena de la Polimerasa (PCR): Haplogrupos Mitocondriales

La reacción que amplifica el ADN, denominada PCR, se llevó a cabo utilizando los cebadores especificados en la Tabla 2 y los siguientes reactivos: MgCl2, Oligonucleótidos trifosfatos, GoTaq PCR Buffer, Agua calidad biología molecular, muestra de ADN antiguo y GoTaq ADN polimerasa, en volúmenes dependientes de las especificaciones del fabricante. Las condiciones del termociclador para realizar la PCR comprende un primer paso de 5 minutos donde se produce la apertura de las cadenas de ADN a 94°C, posteriormente 45 ciclos que incluyen tres pasos: apertura de la cadena 45 segundos a 94°, emparejamiento de los cebadores a las cadenas 45 segundos a 53°C y elongación de la cadena 45 segundos a 73°C; finalmente una extensión durante 5 minutos a 73°C.

Para la determinación de los Hg se utilizó la técnica de Polimorfismo de Longitud de Fragmento de Restricción (en adelante RFLP por sus siglas en ingles) en el cual los productos de amplificación de los Hg A2, C1 y D1 son tratados con la enzima de restricción correspondiente (Tabla 2) según indicaciones del fabricante. Los fragmentos se analizaron en un gel de agarosa o poliacrilamida para verificar la existencia o no del corte enzimático, o en el caso del análisis del Hg B2 la deleción de nueve pares de bases características (Tabla 2). En todos los casos los datos obtenidos corresponden a resultados concordantes de al menos dos amplificaciones provenientes de extracciones diferentes Tipo de muestras y tratamiento utilizado en el análisis de isótopos estables De cada una de las muestras CIB16 y CIB20 se seleccionó un fragmento óseo de la mandíbula y una pieza dental (segundo molar y tercer molar inferior derecho respectivamente) con el fin de analizar el cambio isotópico ocurrido en los primeros y últimos años de vida. En estos fragmentos se analizó la proporción δ13C y δ15N en la fracción Hg orgánica (colágeno en hueso y dentina en el molar). Los análisis se llevaron a cabo siguiendo el protocolo propuesto por Cherkinsky (2009) en el Center for Applied Isotope Studies de la Universidad de Georgia, EE.UU. La proporción 13C/12C del carbono se midió en relación a PDB y la 15N/14N con respecto al aire.

Considerando que este trabajo corresponde a una primera aproximación exploratoria de la dieta de los individuos que se encontraban instalados en la misión Santiago del Baradero, mayormente aborígenes, se recurrió a la búsqueda bibliográfica de datos que pudieran utilizarse para realizar la reconstrucción de la dieta isotópica. Los valores de peces (Loponte y Acosta 2007) se corresponden a muestras arqueológicas de la región de estudio, mientras que los provenientes de novillos (Balasse et al. 2001), maíz (Rossi et al. 2007), porotos y habas (Fraser et al. 2013) pertenecen a muestras de otros países, las cuales fueron elegidas debido a que corresponden a experimentos controlados donde las condiciones de cultivo o pastoreo se asemejan a las que se habrían adoptado durante el siglo XVI y XVII en el área de estudio. Estos alimentos fueron específicamente elegidos ya que se los menciona en las fuentes escritas: “Entro en otra reducción nombrada Santiago del Varadero…tiene los indios que estaban en las islas que serán hasta ciento cuarenta y su cacique general se llama Don Bartholome, y que siembran maíz, frijoles y habillas y otras semillas con lo cual y carne de yeguas cimarronas y pescado que pescan en el rio se sustentan” (AGI –ACH legajo 27, 1619:5).

La interpretación de los datos isotópicos obtenidos se efectuó utilizando el paquete Stable Isotope Analysis in R (SIAR V4, 2011), que se encuentra diseñado para resolver modelos de mezcla para la información isotópica dentro de un marco bayesiano (Inger et al. 2010). Los valores de δ13C de los individuos de Cementerio Indígena y aquellos valores de fauna arqueológica fueron corregidos a valores de dieta según las ecuaciones propuestas por Newsome et al. (2004); además se corrigieron los valores de las muestras modernas a los cuales se les sumo 1,5 debido al efecto Suess (Panarello et al. 2009).

RESULTADOS Y DISCUSION

Datos genéticos

De las 20 muestras analizadas solo obtuvimos resultados de RFLP positivo en 10 de ellas (lo cual representa un porcentaje de recuperación del 50%). De esas 10 muestras positivas, 4 pertenecen a restos óseos mientras que las 6 restantes son piezas dentales. Por lo tanto, al comparar el resultado de los diferentes tipos de muestras observamos que el porcentaje de recuperación de ADNmt antiguo en hueso fue solo del 30,76% mientras que el de las piezas dentales fue del 85,7%. Estos materiales poseen distintas características estructurales: el hueso es más poroso que los dientes y estos últimos además poseen una capa externa de esmalte la cual, sino ha sido dañada por acción mecánica o bacteriana, mantiene a la dentina aislada del medio exterior. Esta característica de las piezas dentales permite una mayor recuperación de material genético.

Nueve de los resultados obtenidos por RFLP corresponden a haplogrupos amerindios (Tabla 3): 2 muestras corresponden a Hg A (CIB8 y CIB15), 5 a Hg B (CIB 11, CIB16, CIB17, CIB18 y CIB 21), 1 a Hg C (CIB10) y 1 a Hg D (CIB20), en tanto que la muestra CIB3 hasta elmomento no ha sido identificada como amerindia (aún resta determinar su tipo de Hg). El número mínimo de individuos asciende a 7, calculando las 6 piezas dentarias que dieron positivo y sumando el Hg C solamente identificado en una piezas ósea. Hasta el momento la información genética obtenida por RFLP ha permitido determinar que la mayoría de los haplogrupos amerindios encontrados pertenecen a B2.

Para la discusión de los resultados se tienen en cuenta los datos genéticos de las poblaciones modernas del litoral mesopotámico, que fueron dados a conocer por varios investigadores. Los análisis genéticos de individuos Mbyá-Guaraní actuales del sur de Brasil indican que los mismos pertenecen a los Hg A y D (Marrero et al. 2007). Sin embargo, en un estudio llevado a cabo por Sala et al. (2010) sobre poblaciones de esa misma etnia pero en la provincia de Misiones, se identificaron los cuatro haplogrupos mitocondriales amerindios, siendo en este caso el más abundante el Hg A (40%) seguido por los Hg D (36%), B (17%) y por último el C. Por otra parte, de acuerdo con los fechados obtenidos en investigaciones arqueológicas de grupos prehispánicos, se estima que los individuos guaraní arribaron a Misiones hace aproximadamente 1000 años siguiendo el curso de los ríos Paraná y Paraguay (Loponte 2008; Loponte y Acosta 2008; Loponte et al. 2011), mientras que a la región del Delta del Río Paraná habrían arribado unos 300 a 200 años antes del arribo de los conquistadores hispánicos (Bonomo et al. 2015).

Los datos genéticos obtenidos hasta el momento de los restos humanos del sitio Cementerio Indígena, indican que la variabilidad de Hg –con un mayor porcentaje de Hg B2- podría ser el resultado de dos situaciones, según los individuos muestrados hayan sido de la etnía guaraní o de los chaná. En primer lugar, considerando que los grupos Mbyá-guaraní del sur brasilero solo presentan dos de los cuatro haplogrupos amerindios mientras que en Misiones se han identificado a todos ellos, se podría inferir que la alta movilidad de esa etnia desde el norte del litoral hasta el río de la Plata podría haber generado una intensificación de la interacción reproductiva con otros grupos étnicos, derivando en una mayor variabilidad genética al instalarse en el Río de la Plata.

En segundo lugar, la identificación preponderante del Hg B es relevante dado que en un estudio realizado por Alechine et al. (2009) sobre poblaciones modernas de Chaco y Formosa, se observó que el Hg B2 era mayoritario entre los grupos Pilagá y Toba del Chaco mientras que D1 era mayoritario en Pilagá y Toba de Formosa. Por ello, es plausible considerar que otras poblaciones del pasado que habitaron en el área de Paraná Inferior como la etnia chaná, también haya presentado mayoritariamente el Hg B. La presencia de los Hg A, C y D en parte pueden dar cuenta de la variabilidad genética de los individuos encontrados en el sitio Cementerio Indígena, hecho esperable en los espacios reduccionales. En efecto, para mantener un número aceptable de nativos que pudiesen sostener el funcionamiento de los pueblos de indios, tanto laboral como económicamente, también se captaba individuos provenientes de grupos indígenas foráneos, procedentes de lugares lejanos como el norte o centro oeste del país.

Otro de los análisis que se llevaron a cabo corresponde a la secuenciación de la Región Hipervariable 1 del ADN mitocondrial para poder determinar el haplotipo específico de estas muestras, sin embargo aún no se han conseguido datos concluyentes. La razón por la cual existen estas complicaciones a la hora de secuenciar puede deberse a las condiciones de preservación. Aquellos restos que se encuentran depositados en las cercanías al río (como es el caso de Cementerio Indígena) están expuestos a situaciones ambientales adversas, como la acidez del suelo y la elevada humedad, que favorecen los procesos diagenéticos y la degradación de las moléculas orgánicas. El ADN es particularmente sensible a estas condiciones de temperatura y humedad, siendo más favorecida su preservación en climas fríos y secos. Considerando estas características es de esperar que existan pocas moléculas de ADN de mayor longitud (como las requeridas para la secuenciación) en aquellos restos obtenidos en zonas de humedales (Hofreiter et al. 2015).

Datos Isotópicos

Los resultados isotópicos de las muestras CIB16 y CIB20 utilizadas para el análisis -dos mandíbulas, de las cuales se extrajeron tanto un fragmento óseo como un molar- indican, en ambos casos, un leve incremento en los valores de δ13C cuando se comparan los primeros años de vida de los individuos con los últimos antes de su muerte y un descenso de los valores de δ15N. Este último valor se encuentra vinculado al consumo de proteínas animales, lo que estaría indicando una disminución de este tipo de alimentos en la dieta de los individuos en los últimos años de vida en la reducción Santiago del Baradero, mientras que el aumento en los valores de δ13C está ligado al consumo de plantas de patrón fotosintético C4 como es el maíz (Tabla 4, Figura 2).

Para establecer la proporción de los distintos alimentos que fueron consumidos por los individuos muestreados, se tomaron los valores isotópicos determinados por otros investigadores. De esta manera se registraron datos para el maíz como representante de las plantas C4 (Rossi et al. 2007) y porotos y habas como representantes de plantas C3 (Fraser et al. 2013), el surubí como representante de los peces (Loponte y Acosta 2007) y de los novillos (Balasse et al. 2001). Utilizando el programa SIAR, que se basa en una metodología bayesiana, se obtuvieron las proporciones de los distintos alimentos seleccionados para el muestreo que habrían contribuido a la dieta de estos individuos.

Como se observa en las Figuras 3 y 4 tanto en los primeros como en los últimos años de vida de estos individuos la dieta fue variable y no se observa la predominancia marcada de una de las posibles fuentes, aunque sí en ambos casos se distingue una leve preferencia hacia el consumo de peces (principalmente) y del maíz, particularmente en las muestras de molares. Estos resultados son contrastantes con la información que proporcionan tres cartas escritas por el Gobernador Góngora al rey de España en las que detalla la situación de la reducción de Santiago del Baradero, registrada durante su visita en 1619 y 1620. En esos documentos se indica que los indígenas consumían vegetales como el maíz, frijoles y habillas y en menor medida el pescado: “por ser malo el sitio, pantanoso y anegadizo sin leña mala agua y tener desviado el pescadero mas de 3 leguas” (AGI-ACH legajo 162:14). Si bien la descripción de Góngora coincide con los valores isotópicos que indican consumo de plantas C4 como el maíz (Figura 5), es contraria respecto del consumo de peces: los valores isotópicos señalan que la diferencia encontrada entre los primeros años de vida de los individuos (antes de entrar en la reducción) y los últimos, no es tan significativa. En consecuencia, aunque los datos de las piezas óseas muestran una leve disminución en el consumo de peces hacía los últimos años de vida, estos recursos habrían continuado siendo un aporte importante en la dieta a lo largo del tiempo (Figura 4 y 5). Por su parte en las muestras de restos óseos puede observarse un leve incremento en el consumo de novillo con respecto a los molares. En ambos casos es destacable la escasa proporción de alimentos con patrón fotosintético C3 caracterizados en este caso por porotos y habas (PlntC3).

Con respecto al análisis de la dieta, se debe tener en cuenta que los individuos chaná y guaraníes que ingresaron a la reducción, tenían una forma de subsistencia tradicional parcialmente diferenciable. Considerando la variedad de recursos disponibles en la zona, las fuentes históricas y los trabajos arqueológicos previos realizados por otros investigadores sobre grupos prehispánicos, se infiere que un componente importante de la dieta para ambas etnias fue aportado por peces y mamíferos terrestres (Loponte y Acosta 2007; Scabuzzo y Gonzalez 2007; Politis et al. 2009; Bo y Quintana 2013). Respecto de los vegetales, para los chaná se ha mencionado la recolección estacional de algunos frutos como los de la palmera pindo (Syagrus ramanzoffiana) y, dada las prácticas hortícolas los guaraníes habrían incluido el maíz y otros recursos como la calabaza. Una vez reducidos en Santiago del Baradero, los individuos de ambas etnías habrían consumido los mismos recursos homogenizando sus dietas; por ello, para los últimos años de vida cabría esperar valores isotópicos similares.

A MODO DE CONCLUSIÓN

Hasta el momento no existían datos genéticos de poblaciones nativas post hispánicas del Paraná Inferior, por lo tanto los datos de ADNmt antiguo obtenidos son los primeros en darse a conocer para el área. Si bien se debe profundizar el análisis con un mayor número de muestras (óseas o dentarias), es posible realizar algunas inferencias iniciales a partir de la identificación de la preponderancia del Hg B2. Este haplogrupo se ha registrado entre algunos grupos amerindios actuales del litoral mesopotámico y quizá también fue preponderante en otros grupos más tempranos como los chaná o bien los descendientes de los Mbayá guaraní.

No obstante, recién se ha empezado a visualizar “la punta del iceberg” de una compleja problemática, ya que la determinación de otros haplogrupos amerindios también indica la existencia de variabilidad genética, un resultado esperado para el caso de todos los “pueblos de indios” donde fue frecuente la relocalización de individuos de diferentes grupos culturales como parte de las estrategias de dominación colonial. De hecho, algunas fuentes mencionan la presencia de indígenas “serranos” y/o “pueblos de la cordillera” (mapuches) de lengua mapudungun, que a mitad del siglo XVII se acercaban a las pampas a intercambiar objetos por ganado: “este gentío de indios (…) con la libertad que goza se retira por tiempos a las cordilleras de Chile, y se comunica con los indios Serranos, y enemigos de aquel Reino, y los proveen de caballos y algunas armas de alfanjes, y espadas anchas y otros géneros” (AGI-ACH 1673:5).

Precisamente, esta referencia documental se relaciona con datos del registro arqueológico del sitio Cementerio Indígena, donde se han recuperado y analizado discos metálicos con punzonado en torno al borde y orificio central (Tapia et al. 2009) los cuales presentan similitudes con otros objetos post hispánicos similares, encontrados en sitios lejanos a Baradero como Mendoza, Neuquén y Río Negro (Tucker et al. 2012; De Rosa y Tapia 2014). El hallazgo de estos discos, que formaron parte de la vestimenta y los adornos, abren interrogantes acerca de si los mismos fueron adquiridos por intercambio con otros grupos aborígenes o si en la Misión también se habrían reducido individuos pertenecientes a grupos étnicos de aquellas zonas (Sanchez Labrador 1936; Schindler 1978). Al respecto, recientemente se iniciaron análisis de ADNmt de muestras de restos humanos encontrados en el sitio Cueva Epullán Grande, ubicado en la provincia de Neuquén (ocupada desde los 2200 AP al periodo poshispánico).

Las dos muestras analizadas para obtener valores isotópicos, también proporcionaron datos genéticos que se corresponden con los haplogrupos amerindios B y D (Tabla 3) y datos cronológicos que los ubican a ambos en los primeros años de vida de la reducción Santiago del Baradero. En estos dos individuos amerindios, se observa que la dieta no habría dependido exclusivamente de un tipo de alimento particular. Si bien en ambas muestras analizadas se observa una disminución en el consumo de proteínas animales y un leve incremento en el consumo de maíz, el cambio de dieta no es el que cabría inferir según las fuentes documentales. Según éstas, se podría esperar una mayor proporción de plantas C3 y por lo tanto un mayor consumo de mamíferos herbívoros terrestres domesticados (novillos, caballos, ovejas). Sin embargo, los datos isotópicos señalan semejanza con la dieta de individuos guaraníes prehispánicos, que incluían el consumo del maíz y la pesca en el modo de subsistencia (Loponte y Acosta 2007). Aunque durante los últimos años de vida en la misión, algunos individuos de tradición cazadora recolectora como los chaná, podrían haber alcanzado una dieta similar.

Para poder interpretar los valores obtenidos se consideran dos posibles explicaciones alternativas. Una de ellas se fundamenta en el hecho que el tiempo de recambio atómico en hueso se estima en aproximadamente 7-10 años, por lo que los resultados isotópicos obtenidos representan un promedio de la dieta a lo largo de este tiempo (Tykot 2014). Por este motivo, los individuos analizados no habrían vivido en la Misión el tiempo suficiente para que se produjera el cambio de registro isotópico de acuerdo al tipo de dieta que menciona el gobernador Góngora, este hecho se corrobora con los dos fechados AMS del Laboratorio de la Universidad de Georgia que ubican los restos para aproximadamente la primera década de vida en Santiago del Baradero. Otra explicación podría apoyarse en el sesgo subjetivo de la observación de Góngora acerca de lo que realmente ocurría en la reducción en esos momentos; después de todo él solo estuvo allí no más de dos días. Por lo tanto, la pesca pudo haber tenido más importancia de lo que se infiere en su relato al rey de España, particularmente si consideramos las frecuentes hambrunas que pasaron los habitantes de esa misión franciscana (AGIACH 1620, 1621-1622).

Estos primeros aportes genéticos e isotópicos abren nuevos interrogantes a responder mediante futuras indagaciones. En especial, se requiere profundizar la reconstrucción de la ecología isotópica de los tiempos de la colonia en el Río de La Plata y ampliar los análisis de ADNmt antiguo, incluyendo otras poblaciones prehispánicas y post hispánicas del área del Paraná Inferior y del centro del país como referentes comparativos.

AGRADECIMIENTOS

Las investigaciones se realizaron mediante los subsidios otorgados por la AGENCIA-FONCYT al proyecto PICT 2008-1981 y por la SECYT-UNLU a los proyectos de investigación arqueológica de las Programaciones científicas 2010-2012, 20132015 y 2016-2017, dirigidos por la Dra. Alicia H. Tapia. Asimismo se ha recibido el apoyo financiero de CONICET mediante el otorgamiento de una beca doctoral al Lic. Darío Cardozo y de la Fundación Científica Felipe Fiorellino para la compra de insumos en el Laboratorio del CEBBAD de la Universidad Maimónides, donde se realizaron los análisis genéticos. Agradecemos a Raquel Negrão por su ayuda con la traducción al portugués.