SciELO - Scientific Electronic Library Online

 
vol.23 issue3Presence of Non-African Haplotypes Increase Genetic Diversity in Sickle Cell Anemia Patients in ColombiaNew Records and Range Extensions for Birds from the Western Slope of the Eastern Andes, Santander, Colombia author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Acta Biológica Colombiana

Print version ISSN 0120-548X

Acta biol.Colomb. vol.23 no.3 Bogotá Sep./Dec. 2018

https://doi.org/10.15446/abc.v23n3.65407 

Articles

ANÁLISIS DE ADN MITOCONDRIAL EN UNA MUESTRA DE RESTOS ÓSEOS PRECOLOMBINOS DE NORTE DE SANTANDER, COLOMBIA (ÁREA CULTURAL CHITARERA)

Mitochondrial DNA Analysis of a Sample of Pre-Columbian Acient Remains from Norte de Santander, Colombia (Chitarera Cultural Area)

Andrea CASAS-VARGAS1  * 

Liza Maria ROMERO1 

José Vicente RODRÍGUEZ2 

William USAQUÉN1 

1 Grupo de Genética de Poblaciones e Identificación, Instituto de Genética, Universidad Nacional de Colombia, Sede Bogotá. Bogotá, Colombia.

2 Laboratorio de Antropología Física, Universidad Nacional de Colombia, Sede Bogotá. Bogotá, Colombia.


RESUMEN

Los análisis de ADN antiguo (ADNa) han incrementado en los últimos años permitiendo conocer la diversidad genética de las poblaciones precolombinas. En Colombia, existen pocos registros arqueológicos de la población prehispánica del Norte de Santander habitada en el siglo XVI por el grupo Chitarero. Por este motivo, nos propusimos analizar la diversidad genética a partir de secuencias de la región HVRI del ADNmt y determinar sus posibles relaciones con otras comunidades tanto antiguas como contemporáneas. Se analizaron siete individuos precolombinos asociados a este grupo prehispánico, recuperados en los municipios de Cácota y Silos en el departamento de Norte de Santander de los Andes Orientales colombianos, siguiendo criterios estrictos de autenticidad para el ADNa. En todos los individuos se logró identificar el haplogrupo B caracterizado por el polimorfismo en la posición 16217C, siendo éste uno de los más frecuentes en comunidades precolombinas y contemporáneas de los Andes Suramericanos. Este hallazgo indica que este grupo poblacional se encuentra estrechamente emparentado por línea materna, con posibles índices de endogamia, con una probable densidad demográfica baja y una baja diversidad genética, similares a lo observado en comunidades pertenecientes a periodos anteriores como el Formativo. Este grupo precolombino exhibe una de las diversidades genéticas más bajas reportadas en las poblaciones pertenecientes a la familia lingüistica Chibcha. Estos resultados genéticos coinciden con los planteamientos sobre el grupo Chitarero de pertenecer a comunidades pequeñas independientes, con asentamientos dispersos, apartados unos de otros.

Palabras clave: ADN antiguo; ADNmt; Andes; Chitareros; Colombia; HVRI

ABSTRACT

In the last few years there has been an increase in ancient DNA (aDNA) analyses that has allowed shedding light on the diversity of pre-Columbian populations. In Colombia, there are few archaeological records belonging to the prehispanic population from Norte de Santander inhabited in the XVI century by the Chitarero. For this reason, we performed a genetic diversity analysis of the HVRI region of mtDNA in order to determine their possible relationships with other communities both ancient and contemporary. We analyzed seven pre-Columbian individuals belonging to this pre-Hispanic group, recovered from the municipalities of Cácota and Silos in the department of Norte de Santander located at the Colombian Andes, following strict authenticity criteria for aDNA. All individuals were identified as belonging to haplogroup B, characterized by the polymorphism found at position 16217C which is one of the most frequent haplogroups in pre-Columbian and contemporary communities of the South American Andes. This finding suggests that this population group was closely related through its maternal lineage, with possible inbreeding indexes, low population density and therefore low genetic diversity, similar to what is observed in communities belonging to previous periods such as the Formative period. This pre-Columbian group exhibits one of the lowest genetic diversities reported in populations belonging to the Chibcha linguistic family. These genetic results coincide with the views on the Chitarero group as belonging to small independent communities, with dispersed settlements separated among them.

Keywords: ancient DNA; Andes; Chitareros; Colombia; HVRI; mtDNA

INTRODUCCIÓN

Los estudios de la variación del ADN mitocondrial (ADNmt) han mostrado que los nativos americanos pertenecen a uno de los cinco haplogrupos maternos A, B, C, D y X (Schurr et al., 1990; Torroni et al., 1992; Brown et al., 1998). Análisis recientes de mitogenomas han permitido establecer la presencia de linajes más específicos, algunos de distribución Pan-Americana (A2, B2, C1b, C1c, C1d, D1 y D4h3), descritos como linajes fundadores, y otros presentes sólo en algunas zonas del Norte y Sur de América (A2a, A2b, C4c, D2A, D3, D4E1, X2a, X2g, A2ad, A2af1 B2i2, B2J, B2K, B2L, C1b13, D1G, D1J) (Tamm et al., 2007; Perego et al., 2009; De Saint Pierre et al., 2012a; De Saint Pierre et al., 2012b; Gómez-Carballa et al., 2012; Perego et al., 2012; Gómez-Carballa et al., 2015). Así mismo, el análisis de la variación genética de la región control del ADNmt ha permitido identificar sitios polimórficos que comparten los individuos de un mismo haplogrupo fundador, los cuales definen haplotipos (sublinajes) específicos en cada linaje. Los sublinajes pueden ser específicos de grupos poblacionales (Torroni etal., 1993; Malhi etal., 2002) o de familias lingüísticas (Derbeneva et al., 2002), por lo que se utilizan para estudiar con mayor precisión las relaciones genéticas entre las poblaciones.

Por otra parte, el desarrollo de técnicas moleculares en la recuperación de ADN en muestras antiguas como restos óseos, tejidos momificados, dientes y coprolitos humanos, han permitido estudiar de manera directa la genética de las poblaciones antiguas (Páábo, 1985; De La Cruz et al., 2008; Kemp et al., 2010). En Colombia este tipo de análisis en poblaciones precolombinas han ido en aumento, especialmente en los Andes Orientales. Se ha logrando amplificar el ADNmt en diferentes poblaciones, principalmente en comunidades chibchas tales como Muiscas, Guanes y Laches (Uricoechea, 2010; Casas-Vargas et al., 2011; Casas-Vargas et al., 2017); a su vez en grupos precerámicos como Checua (Cundinamarca) (Díaz-Matallana et al., 2016) y del período Herrera o Formativo en Madrid (Cundinamarca) (Silva et al., 2008).

A la llegada de los españoles otra comunidad chibcha importante fueron los Chitareros, por referencias históricas se sabe que se ubicaron hacia el nororiente del país en los departamentos de Norte de Santander, extendiéndose hasta San Cristobal y el Estado de Táchira en Venezuela, hacia el oriente se asentaron en los valles del Chitagá, Silos, Labateca y Toledo. La denominación Chitarero y de otros grupos étnicos, resulta problemática puesto que los españoles generalizaron algunas denominaciones a vastas regiones, cuando en realidad los indígenas se circunscribían a pequeños valles y altiplanos; para nuestro caso los indígenas de la provincia de Pamplona fueron denominados "chitareros" (Aguado [1581]1956; Langebaek, 1992; Rodríguez, 2011) (Fig. 1), porque usaban calabazos atados a la cintura en los cuales almacenaban chicha de maíz y yuca. Se ha descrito que desde 1530 los Chitareros fueron sometidos a un proceso acelerado de conquista y mestizaje por parte de los españoles, cuya consecuencia inmediata fue la rápida extinción de la etnia (Ruge, 2014). Desafortunadamente, no existe un contexto histórico y arqueológico completo que permita inferir acerca del origen y dinámicas poblacionales de este grupo.

Figura 1 Caracterización geográfica del Área Chitarera. Contorno gris: Territorio ocupado por el grupo étnico Chitarero en Colombia y Venezuela en el tiempo de la conquista española. Tomado de (Ruge, 2014). 

El objetivo de esta investigación es caracterizar la composición genética de la comunidad precolombina Chitarera de Norte de Santander, a partir del análisis del ADN mitocondrial en la región hipervariable I (HVRI) y establecer sus posibles relaciones biológicas con otras poblaciones para así contribuir al conocimiento de sus posibles orígenes.

MATERIALES Y MÉTODOS

Población de estudio

Se analizaron siete individuos que hacen parte de la colección arqueológica del Laboratorio de Antropología Física de la Universidad Nacional de Colombia. Todos los individuos fueron encontrados en abrigos rocosos de los municipios de Silos y Cácota (Norte de Santander, Colombia) en los años 2002 y 2006 en marco del Proyecto Inventario Arqueológico, Paleontológico y de Arte Rupestre de Norte de Santander (Ruge, 2014). Se obtuvieron muestras de dientes y huesos largos para el análisis molecular.

Desde el punto de vista arqueológico los lugares donde fueron encontrados los restos óseos corresponden en su mayoría a abrigos rocosos que fueron utilizados como recintos funerarios, similar a lo encontrado en la Mesa de Los Santos (Santander) y en Silos (Norte de Santander) (Silva, 1947). Las estimaciones de sexo y edad se obtuvieron de los análisis bioarqueológicos previos. No se cuenta con una datación de radiocarbono de las muestras que permita situar cronológicamente a estos individuos, pero a juzgar por el ajuar corresponden al período Reciente (siglos X-XVI d. C.) (Tabla 1).

Tabla 1 Muestras de restos óseos analizadas. Codificación de las muestras según ficha, localización del yacimiento y estimación de sexo y edad a partir de datos antropológicos. 

Precauciones de contaminación para el análisis de ADN antiguo

Para garantizar los resultados y evitar cualquier contaminación exógena, se siguieron algunos de los protocolos propuestos por Cooper y Poinar (2000). Las muestras fueron analizadas en un laboratorio ubicado en el Instituto de Genética de la Universidad Nacional de Colombia que es exclusivo para el análisis de los restos óseos antiguos. Se contó con áreas de trabajo separadas, que garantizaron todas las condiciones de asepsia para el análisis de cada una de las muestras, material desechable (batas, guantes, tapabocas y gorro) y limpieza e irradiación con luz UV de los materiales antes de su uso.

Para detectar la posible contaminación por ADN exógeno moderno, blancos de extracción y amplificación se utilizaron como controles negativos.

Para cada una de las muestras procesadas, se realizaron réplicas de los análisis (de dos a tres aislamientos y amplificaciones por cada muestra). Adicionalmente, se realizó reproducibilidad de los resultados de una muestra escogida al azar en otro laboratorio de ADNa ubicado en la Universidad de la Sabana (Bogotá, Colombia). Por último, todas las personas involucradas en el proceso de las muestras, se tipificaron genéticamente (secuenciación HVR I) y fueron comparados con los resultados obtenidos.

Extracción de ADN

Cada muestra fue sometida a una limpieza antes de su procesamiento como se ha descrito previamente (Casas-Vargas et al., 2011; Casas-Vargas et al., 2017). Brevemente, la descontaminación se realizó removiendo completamente la capa externa del hueso empleando un mototool (para las piezas dentales no fue necesario realizar este primer proceso). Posteriormente las muestras fueron lavadas con hipoclorito al 5 % por diez minutos, seguido de etanol al 70 % por cinco minutos y finalmente exposición a luz ultravioleta por cada uno de los lados durante 15 minutos. Finalmente, se procedió a obtener aproximadamente dos gramos de pulverizado de hueso o diente usando el molino MM400 de Retsch.

Aproximadamente un gramo de pulverizado fue incubado a 56 °C durante 18 horas en una solución de 15 mL que contenía EDTA 0.5M, N-Lauryl-Sarcosyl al 1 %, y proteinasa-K (20 mg/mL). Luego se realizó la extracción del ADN empleando el kit Qiaquick PCR Purification (Qiagen) siguiendo las recomendaciones del fabricante y las modificaciones realizadas por la Comisión Internacional de Personas Desaparecidas (International Commission on Missing Persons [ICMP] (Huffine et al., 2001; Parsons et al., 2007; Huel et al., 2012). Brevemente, el lisado obtenido se centrifuga a 4000 rpm durante cinco min, el sobrenadante se transfiere a filtros Amicon y se centrifuga a 4000 rpm para concentrar la muestra hasta que quede un volumen final de 250 uL. El volumen recuperado se transfiere a tubos de 1,5 mL que contienen 1,250 uL de buffer PB y alícuotas de 720 uL se centrifugan en filtros QIAquick a 14000 rpm por un min. A continuación, se añaden 720uL de buffer PE al filtro y se centrifugan a 14000 rpm por un min. El ADN se eluye mediante la adición de buffer EB, se incuba durante dos min, y se centrifuga a 14000 rpm durante un minuto.

Este proceso se repitió durante dos ocasiones para algunas de las muestras escogidas al azar en momentos separados para confirmar los resultados obtenidos. Los aislados de ADN fueron cuantificados empleando el equipo Qubit® 2.0 Fluorometer de Invitrogen, siguiendo las recomendaciones del fabricante. Para la reacción de PCR fueron empleados de 3 a 7 ng de ADN.

Como control de contaminación, se realizó la extracción de ADN de los investigadores que intervinieron en el proceso, a partir de sangre empleando tarjetas FTA® (Whatman Bioscience, Cambridge, United Kingdom), las cuales están diseñadas para la recolección, almacenamiento y posterior purificación del ADN, mediante el uso de un buffer FTA ®.

Análisis del ADN mitocondrial

Los aislados de ADN fueron amplificados para cuatro fragmentos de PCR de tamaños entre los 160 y 180 pb que amplifican la región hipervariable I del ADN mitocondrial, (Tabla S1). El ADN fue amplificado mediante PCR utilizando master mix HotStarTaq Master Mix (Qiagen®) 2X, 1X de enhancer (Epicentre®),100 ug de BSA (Roche®) y 0,8uM de cada cebador o primer. El programa de PCR utilizado considera: denaturación inicial, 95 °C por 15 min, 50 ciclos de: denaturación, 95 °C por 30 seg; anillamiento, 30 seg a una temperatura de 57- 58 °C dependiendo de los primers empleados; elongación, 72 °C por 30 seg; y elongación final a 72 °C por diez minutos.

La presencia e integridad de los amplificados obtenidos en las PCRs fueron verificadas por electroforesis en gel de agarosa al 1,5 % teñido con Syber® Safe. Se confirmó la concentración y altura de la banda esperada por comparación con el patrón de peso molecular HyperLadderTM II (Bioline). La purificación de cada fragmento amplificado se realizó con el Kit WizardTM PCR Prep System (Promega), siguiendo las recomendaciones del fabricante. Los productos de purificados de PCR fueron enviados a MACROGEN (Seoul, South Korea). Se realizaron secuencias bidireccionales empleando los mismos primers de la amplificación de por lo menos dos-productos de PCR por cada muestra para cada uno de los fragmentos (con el fin de descartar errores). Para las muestras control de los investigadores se emplearon los primers L15,997 y H16,391 que amplifican un fragmento de 433 pb de la HVR I del ADN mitocondrial (Sánchez, 2007).

Análisis de la deleción de 9pb

La confirmación de los resultados de la secuenciación se hizo mediante el análisis de la deleción de 9 pb en COII/ tARNlys, característica del haplogrupo B. El ADN fue amplificado mediante PCR utilizando master mix GoTaq Hot Start® ADN polimerasa (Qiagen®) 2X, 1X de enhancer (Epicentre®), y 100 ug de BSA (Roche®) y 0,5uM de cada cebador o primer (Tabla S1). El programa de PCR utilizado considera: denaturación inicial, 95 °C por 15 min, 50 ciclos de: denaturación, 95 °C por 45 seg; anillamiento, 55 °C por 45 seg; elongación, 72 °C por 45 seg.; y elongación final a 72 °C por tres minutos. La deleción fue evidenciada en gel de poliacrilamida:bisacrilamida en concentración 40:1 teñido con Syber® Safe.

Amplificación del ADN para marcadores autosómicos y amelogenina

La amplificación se realizó empleando el kit comercial multiplex AmpFiSTR® SEfilerTM (Applied Biosystems) el cual cuenta con 11 marcadores STRs autosómicos (D2S1338, D3S1358, D8S1179, D16S539, D18S51, D19S433, D21S11, FGA, SE33, TH01 y vWA) y el marcador sexual Amelogenina. En un volumen final de 25 uL se mezclaron los siguientes componentes: mezcla de reacción para PCR (AmpFiSTR® SEfilerTM MIX: MgCl2, dNTPs, albúmina de suero bovino en buffer con azida de sodio 0.05 %) 10 uL; set de primers (forward y reverse) 5 uL; enzima ADN polimerasa (AmplitaqGoldTM DNA Polimerase 5 U/uL) 0,5 uL; ADN (1,0-2,5 ng/uL) 5 uL y 4,5 uL de agua destilada desionizada. La PCR se llevó a cabo en un termociclador Thermal cycler ®PCRSystem 2720 de Applied-Biosystems®, siguiendo las recomendaciones del fabricante. Se utilizó un Analizador Genético ABI PRISM® 310 (Applied Biosystems) para la tipificación molecular siguiendo las recomendaciones del fabricante. Se empleó el software Genemapper ID v3.2, el cual analiza el estándar de tamaño de control interno GeneScan-500 LIZ® en cada muestra y la escalera alélica para asegurar la precisión en cada inyección del capilar. La asignación alélica de las muestras analizadas se realizó por comparación de los fragmentos obtenidos en pares de bases para cada pico del alelo con los picos de la escalera alélica provista por la casa comercial para el kit SEfi ler ® (Applied Biosystems®, Foster City, CA).

Análisis de datos

Para la asignación de los haplotipos, las secuencias fueron comparadas con la Secuencia de Referencia Revisada de Cambridge (rCRS) (Andrews et al., 1999), en el software Sequencher 5.2.4 (Gene Codes Corporation, Ann Arbor, MI, USA). La asignación del haplogrupo fue acorde a Phylotree build 17 (Van Oven y Kayser, 2009) (Fig. S1).

La frecuencia de los haplogrupos se calculó dividiendo el número de individuos de cada haplogrupo entre el número total de individuos analizados.

Para establecer la posible relación biológica entre los individuos precolombinos encontrados en el área Chitarera con otras comunidades nativas antiguas y contemporáneas se compararon los datos obtenidos con 31 comunidades nativas americanas contemporáneas, y con 34 poblaciones antiguas precolombinas suramericanas (Tabla S2). Se realizó análisis de UPGMA basado en las frecuencias de haplogrupo de las diferentes poblaciones amerindias colombianas empleando el software MVSP® v3.1 (Kovach, 2007).

Las distancias genéticas se estimaron con base en las frecuencias de cada haplogrupo del ADNmt, calculando los datos del índice de fijación (F) con el programa Arlequín® (Excoffier et al. 2005). Se compararon las poblaciones por medio del análisis de escalamiento multidimensional basado en los datos del FST, con el programa MVSP® , versión 3.1 (Kovach, 2007).

RESULTADOS

Marcadores autosómicos y amelogenina

Para las siete muestras de restos óseos analizadas, se realizó la tipificación de 12 marcadores autosómicos con el kit comercial multiplex AmpFiSTR® SefilerTM (Applied Biosystems®). Este análisis se realizó con el fin de garantizar que se trataran de individuos únicos y no se tomaran dos piezas óseas de un mismo individuo. La determinación del sexo con el marcador de la amelogenina fue utilizada para comparar la evidencia genética con los datos arqueológicos según sus características morfológicas, encontrando concordancia en todas ellas. En tres muestras se logró obtener perfiles en 11 STRs (individuos B15, B17 y B49); en la muestra B47 no se obtuvo amplificación debido al estado de degradación del ADN, lo cual es característico de las muestras antiguas (Tabla S3), para el individuo B48 fue posible determinar el sexo femenino (XX), aportando datos genéticos a los análisis antropológicos previos. Finalmente, los perfiles de ADN encontrados en los restos óseos no coincidieron con los de los investigadores que manipularon el material durante la experimentación (Tabla S3).

Análisis de haplogrupos mitocondriales

La metodología de extracción de ADNa implementada en este trabajo presentó una reproducibilidad del 100 % la cual es bastante alta dada la complejidad del procedimiento realizado, es decir, que se logró extraer y cuantificar ADN de las siete muestras analizadas (Tabla S3). Estos resultados mostraron cantidades suficientes de ADN para su amplificación y asegurar la repetibilidad de los mismos.

La asignación del haplogrupo mitocondrial se realizó siguiendo los lineamientos de Phylotree build 17 (Van Oven y Kayser, 2009). Los siete individuos presentaron el mismo haplotipo con la presencia de un solo polimorfismo en la posición 16217C característicos del haplogrupo B2 (Tabla 2).

Tabla 2 Haplotipos mitocondriales región HVRI de las muestras analizadas. investigadores involucrados en la manipulación del material óseo. 

La amplificaron y secuenciación de HVRI de los tres fragmentos correspondientes a las regiones 1 (16,011 - 16,131), 3 (16,210-16,330) y 4 (16,250-16,382), fueron exitosos. Sin embargo, la secuenciación del fragmento 2 (16,127-16,229) provoca un deslizamiento de la polimerasa, lo cual, no hace fácil la interpretación de los polimorfismos en esta región. Por este motivo el polimorfismo en la posición T16,189C, característico del haplogrupo B4 no fue detectado (Tabla 2).

Con el fin de confirmar este haplogrupo se realizó el análisis de la delecion de 9pb, en COII/tARNlys característica del haplogrupo B; todas las muestras presentaron esta deleción.

Teniendo en cuenta que en poblaciones Nativo Americanas el haplogrupo B predominante es el B2, caracterizado por la deleción de nueve pares de bases en la región COII/ tARNlys y presencia del polimorfismo en la posición 16217 (De Saint Pierre et al., 2012a; De Saint Pierre et al., 2012b), designaremos las muestras como haplogrupo B2.

Relación genética de los Chitareros con otros grupos étnicos colombianos: un análisis regional

En Colombia, los estudios de ADNa realizados con restos humanos han involucrado principalmente poblaciones andinas del Altiplano Cundiboyacense (Silva et al., 2008; Casas-Vargas et al., 2017), Santander (Casas-Vargas et al., 2011) y el Valle del Cauca (Barreto et al., 2010). Para Norte de Santander no se tienen registros, siendo una zona culturalmente importante para el país debido a que fue un área de asentamiento de poblaciones prehispánicas de diversa filiación étnico-cultural y lingüística (Rodríguez, 2007).

Al construir un dendograma UPGMA basado en las frecuencias de los haplogrupos de las diferentes poblaciones amerindias antiguas y actuales colombianas, se observó una estrecha relación genética de los individuos del estudio (Chit-COL-ANT), con tres poblaciones antiguas pertenecientes a los tres periodos de ocupación precolombina: Precerámico (Precer- COL-ANT), Formativo (Herr-COL-ANT) y Reciente (Quimb- COL-ANT) (Fig. 2).

Figura 2 Dendograma NJ construido a partir de las frecuencias de haplogrupos mitocondriales de las muestras del presente proyecto comparado con muestras de otras poblaciones amerindias antiguas y actuales de Colombia. Los nombres claves de las abreviaturas están indicados en paréntesis: ChibCOL (ChibchaCOL), CubCOL (CubeoCOL), Paez COL (PaezCOL), DesCOL (DesanoCOL), CurrCOL (CurripacoCOL), WayCOL (WayuuCOL), WauCOL (WaunanaCOL), IngCOL (InganoCOL), Tu-CuCOL (Tule-CunaCOL), TucCOL (TucanoCOL), EmbCOL (EmberaCOL), AntCOL (AntioquiaCOL), TatCOL (TatuyoCOL), ZenuCOL (ZenuCOL), PiuCOL (PuinaveCOL), CauCOL (CaucaCOL), Guan-ButarCOL (Guane-ButareguaCOL), GuamCOL (GuambianoCOL), ArsCOL (ArsarioCOL), KunCOL (KunaCOL), KogCOL (KoguiCOL), ArhCOL (ArhuacoCOL), IjkaCOL (IjkaCOL), ChimCOL (ChimilaCOL), PiaCOL (PiapocoCOL), PirCOL (PiratapuyoCOL), YerCOL (YeralCOL), TuyCOL (TuyucaCOL), BarCOL (BarasanoCOL). 

Chitareros y la dinámica poblacional precolombina en Suramérica

Se construyó una matriz de distancias a partir de datos del FST y se analizaron los componentes principales (PCA) (Fig. 3). Se observa la cercanía de los chitareros (No. 28) con la población antigua de Madrid, Cundinamarca (No. 30) clasificados en el período Formativo o Herrera (Fig. 3).

Figura 3 Análisis de componentes principales (PCA) basado en datos del FST de poblaciones de Precolombinas de Colombia y Sur América para haplogrupos de ADNmt. Las poblaciones enumeradas se encuantran descritas en la Tabla S2. 

DISCUSIÓN

Este estudio constituye el primer esfuerzo en la realización de estudios de ADN antiguo en el Nororiente colombiano, suministrando información complementaria a las reseñas históricas y los datos antropológicos recolectados previamente para este grupo étnico. El análisis molecular permitió clasificar a todos los individuos chitareros estudiados como pertenecientes al haplogrupo B2 por la presencia del polimorfismo en la posición T16, 217C respecto a la secuencia de referencia rCRS. Todos presentaron el mismo haplotipo indicando posiblemente que se tratara de un grupo humano estrechamente emparentado, con un alto nivel de endogamia, teniendo en cuenta una probable densidad demográfica baja y por lo tanto una baja diversidad genética, similar a lo encontrado en otras poblaciones precolombinas como en Madrid (Cundinamarca-Colombia) (Periodo Formativo o Herrera) (Silva et al., 2008), lo cual sería un dato valioso y aportaría a la información etno-histórica disponible hasta la fecha. Un factor importante que debemos tener en cuenta es que la muestra no sea realmente representativa de la población conocida como Chitareros asentada en los municipios de Santo Domingo de Cácota y Silos en Norte de Santander durante el periodo Reciente o chibcha; a pesar de ello, estos datos representan los primeros hallazgos moleculares que complementarían los análisis craneométricos y dentales que indican una gran afinidad con el grupo Muisca del Norte (Rodríguez, 2011).

Según la información de los cronistas sobre los grupos chitareros del siglo XVI, ubicados entre los ríos Guaca, Servitá y Suratá (Norte de Santander), los describen como comunidades independientes, con asentamientos dispersos, pequeños, apartados unos de otros, en distintos pisos térmicos, que no superaban la cifra de 30.000 habitantes (Aguado [1581] 1956; Langebaek, 1992; Rodríguez, 2001). Por el contrario, los tamaños poblaciones de las comunidades precolombinas aledañas son superiores a estos valores como es el caso de los muiscas en donde se ha calculado una cifra de 600.000-900.000 habitantes para el siglo XVI, para los guanes 100.000-150.000 habitantes (Rodríguez, 2011), esto afectaría la diversidad genética entre cada uno de los grupos poblacionales.

El patrón de distribución de frecuencias de la población precolombina Chitarera es muy similar al encontrado en otro estudio realizado en restos óseos antiguos del periodo Herrera ó Formativo de la sabana de Bogotá en donde se analizaron seis individuos y todos ellos presentaron el mismo haplogrupo mitocondrial (B) (Silva et al., 2008). Caso contrario sucede cuando se compara con otros grupos cercanos geográfica y culturalmente como los Muiscas y los Guanes (Casas-Vargas et al., 2011; Casas-Vargas et al., 2017), donde se ha reportado una mayor variabilidad genética con altas frecuencias de haplogrupo A, frecuencias intermedias de haplogrupo B, y en menores frecuencias los linajes Cy. Los hallazgos de esta distribución de frecuencias en poblaciones antiguas colombianas ha sido atribuida a un arribo inicial de individuos que traían consigo los haplogrupos B y D, y posteriormente un flujo génico más reciente de pobladores que presentaban una alta frecuencia del haplogrupo A (Jara, 2007).

Moraga et al. (2005) han descrito una frecuencia del haplogrupo B en un 42 % en población antigua (1800-1300 AP) del norte de Chile. En muestras precolombinas de la región de la Pampa Grande (Salta, Argentina), se reportó este haplogrupo en un 42 % de las muestras (Carnese et al., 2010), reflejando un flujo genético entre las tierras altas a los valles de la región del sur del complejo Andino. Altas frecuencias del haplogrupo B han sido detectadas en poblaciones modernas de los Andes Centrales, como los Aymara, Quechua, Tupe, Yungay, Arequipa (Fuselli et al., 2003; Lewis et al., 2007). Este haplogrupo ha sido encontrado también en todos los periodos arqueológicos colombianos (Precerámico, Herrera, Reciente), sugiriendo una continuidad poblacional a lo largo del tiempo (Fernández, 1999; Silva et al., 2008; Casas-Vargas et al., 2011; Casas-Vargas et al., 2017).

El estudio de la variación genética en el ADNmt en poblaciones de nativos americanos, ha dado lugar a una extensa discusión acerca del poblamiento de América (Powell y Neves 1999; Perez et al., 2009; Chatters et al., 2014). Cuando se compara la diversidad y la estructura genética de poblaciones andinas y amazónicas se observa que las primeras tienen una mayor diversidad genética dentro de las poblaciones, pero pequeñas distancias genéticas, en contraste con las amazónicas (Luiselli et al., 2000; Tarazona-Santos et al., 2001; Fuselli et al., 2003; Lewis et al., 2005). Adicionalmente, la distribución de las frecuencias de los haplogrupos del ADNmt (A-D) en poblaciones actuales, ha mostrado que las Andinas tienen un alto porcentaje del haplogrupo B (63 %), mientras que los otros haplogrupos se encuentran en menor proporción, en comparación con las amazónicas, en las cuales el haplogrupo C, tiene la mayor frecuencia (Rodriguez-Delfin et al., 2001).

Los datos de este estudio muestran la importancia del ADNa en los estudios arqueológicos, y constituyen una primera aproximación genética en esta comunidad precolombina del Norte de Santander. Es necesario ampliar el rango de lectura en el ADN mitocondrial, con el fin de dar mayor poder de resolución y lograr retratar sublinajes mas informativos, profundizando de esta manera en el conocimiento de los origenes y dinámicas poblacionales de estas comunidades.

Finalmente, uno de los principales problemas en los análisis de ADNa humano es la abundancia del ADN humano moderno, debido a esto es necesario adoptar ciertos criterios de autenticidad. Sin embargo, muchos de sus criterios no están al alcance de los investigadores. Gilbert et al. (2005) argumentan que los investigadores en este campo deben tomar una aproximación más cognoscitiva y autocrítica, especialmente, en el momento de comprobar en detalle todos los criterios de la lista, los investigadores pueden explicar cómo los datos fueron obtenidos y por qué ellos creen que son auténticos (Gilbert et al., 2005). Por lo anterior, los resultados que se obtuvieron fueron realizados teniendo en cuenta por lo menos siete de los criterios sugeridos por Cooper y Poinar (2000).

Esto conllevó a ser más rigurosos en los análisis con el fin de descartar cualquier tipo de contaminación cruzada entre las muestras o con otro tipo de ADN exógeno. Dado que en muchas ocasiones los haplogrupos de ADN mitocondrial no permiten individualizar, fue necesario emplear marcadores autosómicos de tipo STR's para garantizar que se trataran de individuos únicos y que los perfiles genéticos no correspondieran a algunos de los investigadores que intervinieron en el proceso. Así mismo los resultados fueron evaluados bajo criterios de repetibilidad (cambio de analista) y reproducibilidad (laboratorio alterno de ADNa). En ambas situaciones se presentaron los mismos resultados brindándole un mayor soporte a nuestros resultados.

CONCLUSIONES

Nuestros datos resaltan la importancia de la contribución del ADN antiguo en los estudios arqueológicos. Durante este estudio se implementó una metodología con alta reproducibilidad para la extracción de ADNa, que permitió un análisis del ADN mitocondrial y la confirmación molecular del sexo antropológico. De los cuatro haplogrupos mitocondriales reportados en la literatura se determinó la presencia del linaje B en los siete individuos precolombinos con fi liación cultural Chitarera. Adicionalmente, se encontró una fuerte asociación de este grupo con otras poblaciones antiguas andinas de diferentes periodos antropológicos sugiriendo una afinidad poblacional a lo largo del tiempo.

Para comprender mejor la dinámica poblacional de esta comunidad, se hace necesario ampliar el tamaño de la muestra, explorar un mayor territorio y analizar marcadores genéticos adicionales. Así mismo en este estudio se resalta la importancia de que exista una contextualización arqueológica previa al análisis molecular, ya que sin ella las inferencias genéticas acerca de sus relaciones biológicas y comportamientos culturales (endogamia o rituales funerarios) pueden ser imprecisas e ineficientes.

La información genética obtenida en el presente estudio puede contribuir a los arqueólogos y antropólogos interesados en estudiar esta zona del país y que deseen contextualizar mejor a este grupo étnico, en conjunto con los datos obtenidos se aumentará la consistencia de los resultados y aportarán mucha más información sobre la historia poblacional de este grupo étnico.

AGRADECIMIENTOS

Los autores expresan un agradecimiento especial al Laboratorio de Antropología Física de la Universidad Nacional de Colombia por facilitar el acceso a las muestras que hicieron posible este trabajo. También al ICANH, por autorizar el acceso a los restos óseos (Licencia No 4214 de 2014) en el marco del proyecto financiado por Colciencias (ref. 110156934771). Finalmente, a la Universidad de la Sabana, especialmente al Dr. Ignacio Briceño y María Claudia Noguera por realizar la reproducibilidad del estudio.

REFERENCIAS

Aguado P. Recopilación historial (1581) (Introducción, notas comentarios de Juan Friede). 4 vols. Bogotá D.C.: Biblioteca de la Presidencia de la Repúplica de Colombia. 1956. [ Links ]

Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999;23(2):147. Doi:10.1038/13779. [ Links ]

Barreto G, Rodríguez C, Braga Y, Quintero M, Jurado L, Arias L, et al. Extracción y amplificación de ADN arqueológico de dos poblaciones prehispánicas del suroccidente colombiano. Int J South Am Archaeol. 2010;6:46-57. [ Links ]

Brown MD, Hosseini SH, Torroni A, Bandelt HJ, Allen JC, Schurr TG, et al. mtDNA haplogroup X: An ancient link between Europe/Western Asia and North America? Am J Hum Genet. 1998;63(6):1852-1861. Doi:10.1086/302155. [ Links ]

Carnese FR, Mendisco F, Keyser C, Dejean CB, DugoujonJ-M, Bravi CM, et al. Paleogenetical study of pre-Columbian samples from Pampa Grande (Salta, Argentina). Am J Phys Anthropol. 2010;141(3):452-462. Doi:10.1002/ajpa.21165. [ Links ]

Casas-Vargas A, Gómez A, Briceño I, Díaz-Matallana M, Bernal JE, Rodríguez JV. High genetic diversity on a sample of pre-Columbian bone remains from Guane territories in northwestern Colombia. Am J Phys Anthropol . 2011;649:637-649. Doi:10.1002/ajpa.21626. [ Links ]

Casas-Vargas A, Romero L, Usaquen W, Zea S, Silva M, Briceño I, et al. Diversidad del ADN mitocondrial en restos óseos prehispánicos asociados al templo del sol en los andes orientales colombianos. Biomédica. 2017;37(4):1-9. [ Links ]

Chatters JC, Kennett DJ, Asmerom Y, Kemp BM, Polyak V, Blank AN, et al. Late Pleistocene human skeleton and mtDNA link Paleoamericans and modern Native Americans. Science. 2014;344(6185):750-754. Doi:10.1126/science.1252619. [ Links ]

Cooper A, Poinar HN. Ancient DNA: do it right or not at all. Science. 2000 p. 1139. Doi:10.1126/science.289.5482.1139b. [ Links ]

De La Cruz I, Gonzalez-Oliver A, Kemp BM, Roman JA, Smith DG, Torre-Blanco A. Sex identification of children sacrificed to the ancient Aztec rain gods in Tlatelolco. Curr Anthropol. 2008;49(3):519-526. Doi:10.1086/587642. [ Links ]

De Saint Pierre M, Bravi CM, Motti JMB, Fuku N, Tanaka M, Llop E, et al. An alternative model for the early peopling of southern South america revealed by analyses of three mitochondrial DNA haplogroups. PLoS One. 2012a;7(9):e43486. Doi:10.1371/journal.pone.0043486. [ Links ]

De Saint Pierre M, Gandini F, Perego UA, Bodner M, Gómez-Carballa A, Corach D, et al. Arrival of Paleo-Indians to the southern cone of South America: new clues from mitogenomes. PLoS One . 2012b;7(12):e51311. Doi:10.1371/journal.pone.0051311. [ Links ]

Derbeneva OA, Sukernik RI, Volodko N V., Hosseini SH, Lott MT, Wallace DC. Analysis of mitochondrial DNA diversity in the Aleuts of the Commander Islands and its implications for the genetic history of Beringia. Am J Hum Genet . 2002;71(2):415-421. Doi:10.1086/341720. [ Links ]

Díaz-Matallana M, Gómez A, Briceño I, Rodríguez JV. Genetic analysis of Paleo-Colombians from Nemocón, Cundinamarca provides insights on the early peopling of Northwestern South America. Rev Acad Colomb Cienc. 2016;40(156):461-483. [ Links ]

Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1:47-50. [ Links ]

Fernández C. La arqueología molecular aplicada a la solución de problemas prehistoricos: análisis de ADN mitocondrial en momias y restos óseos prehispánicos (Tesis de pregrado). Bogotá, Departemento de Antropología, Facultad de Ciencias Humanas, Universidad Nacional de Colombia; 1999. 196 p. [ Links ]

Fuselli S, Tarazona-Santos E, Dupanloup I, Soto A, Luiselli D, Pettener D. Mitochondrial DNA diversity in South America and the genetic history of andean highlanders. Mol Biol Evol. 2003;20(10):1682-1691. Doi:10.1093/molbev/msg188. [ Links ]

Gilbert MTP, Hofreiter M, Barnes I. Assessing ancient DNA studies. 2005;20(10):541-544. Doi:10.1016/j.tree.2005.07.005. [ Links ]

Gómez-Carballa A, Ignacio-Veiga A, Alvarez-Iglesias V, Pastoriza-Mourelle A, Ruíz Y, Pineda L, et al. A melting pot of multicontinental mtDNA lineages in admixed Venezuelans. Am J Phys Anthropol . 2012;147:78-87. Doi:10.1002/ajpa.21629. [ Links ]

Gómez-Carballa A, Catelli L, Pardo-Seco J, Martinón-Torres F, Roewer L, Vullo C, et al. The complete mitogenome of a 500-year-old Inca child mummy. Sci Rep. Nature Publishing Group; 2015;5:16462. Doi:10.1038/srep16462. [ Links ]

Huel R, Amory S, Bilic A, Vidovic S, Jasaragic E, Parsons T. DNA extraction from aged skeletal samples for STR typing by capillary electrophoresis. Methods Mol Biol. 2012;830:185-198. [ Links ]

Huffine E, Crews J, Kennedy B, Bomberger K, Zinbo A. Mass identification of persons missing from the break-up of the former Yugoslavia: structure, function, and role of the International Commission on Missing Persons. Croat Med J. 2001;42(3):271-275. [ Links ]

Jara N. Determinación de la estructura genética en un grupo poblacional Muisca mediante el análisis de polimorfismos en el ADNmt (tesis de maestría). Bogotá: Instituto de Genética, Facultad de Medicina, Pontificia Universidad Javeriana; 2007. 152 p. [ Links ]

Kemp BM, González-Ol iver A, Malhi RS, Monroe C, Schroeder KB, McDonough J, et al. Evaluating the Farming/Language Dispersal Hypothesis with genetic variation exhibited by populations in the Southwest and Mesoamerica. Proc Natl Acad Sci. 2010;107(15):6759-6764. Doi:10.1073/pnas.0905753107. [ Links ]

Kovach W. MVSP-A Multivariate Statistical Package for Windows, ver. 3.1. Pentraeth, Wales, U.K.: Kovack Computing Services; 2007. [ Links ]

Langebaek C. Noticias de caciques muy mayores: origen y desarrollo de sociedades complejas en el nororiente de Colombia y norte de Venezuela. 1a ed. Bogotá: Ediciones Uniandes-Editorial Universidad de Antioquia; 1992. 256 p. [ Links ]

Lewis CM, Tito RY, Lizárraga B, Stone AC. Land, language, and loci: mtDNA in Native Americans and the genetic history of Peru. Am J Phys Anthropol. Wiley Subscription Services, Inc., A Wiley Company; 2005;127(3):351-360. Doi:10.1002/ajpa.20102. [ Links ]

Lewis CM, Lizárraga B, Tito RY, López PW, Iannacone GC, Medina A, et al. Mitochondrial DNA and the peopling of South America. Hum Biol. 2007;79(2):159-178. [ Links ]

Luiselli D, Simoni L, Tarazona-Santos E, Pastor S, Pettener D. Genetic structure of Quechua-speakers of the Central Andes and geographic patterns of gene frequencies in South Amerindian populations. Am J Phys Anthropol . 2000;113:5-17. Doi:10.1002/1096-8644(200009)113:1<5::AID-AJPA2>3.0.CO;2-U. [ Links ]

Malhi RS, Eshleman JA, Greenberg JA, Weiss DA, Schultz Shook BA, Kaestle FA, et al. The structure of diversity within New World mitochondrial DNA haplogroups: implications for the prehistory of North America. Am J Hum Genet . 2002;70(4):905-919. [ Links ]

Moraga M, Santoro CM, Standen VG, Carvallo P, Rothhammer F. Microevolution in prehistoric Andean populations: chronologic mtDNA variation in the desert valleys of northern Chile. Am J Phys Anthropol . 2005;127(2):170-181. Doi:10.1002/ajpa.10438. [ Links ]

Pääbo S. Molecular cloning of Ancient Egyptian mummy DNA. Nature. 1985;314(6012):644-545. Doi:10.1038/314644a0. [ Links ]

Parsons TJ, Huel R, Davoren J, Katzmarzyk C, Milos A, Selmanovic A, et al. Application of novel "mini-amplicon" STR multiplexes to high volume casework on degraded skeletal remains. Forensic Sci Int Genet. 2007;1:175-179. Doi:10.1016/j.fsigen.2007.02.003. [ Links ]

Perego UA, Achilli A, Angerhofer N, Accetturo M, Pala M, Olivieri A, et al. Distinctive Paleo-Indian migration routes from beringia marked by two rare mtDNA haplogroups. Curr Biol. 2009:1-8. Doi:10.1016/j.cub.2008.11.058. [ Links ]

Perego U, Lancioni H, Tribaldos M, Angerhofer N, Ekins JE, Olivieri A, et al. Decrypting the mitochondrial gene pool of modern Panamanians. PLoS One . 2012;7(6):e38337. Doi:10.1371/journal.pone.0038337. [ Links ]

Perez SI, Bernal V, Gonzalez PN, Sardi M, Politis GG. Discrepancy between cranial and DNA data of early Americans: implications for American peopling. PLoS One . 2009;4(5):e5746. Doi:10.1371/journal.pone.0005746. [ Links ]

Powell JF, Neves W. Craniofacial morphology of the first Americans: Pattern and process in the peopling of the New World. Am J Phys Anthropol . 1999;Suppl 29:153-188. [ Links ]

Rodriguez-Delfin LA, Rubin-de-Celis VE, Zago MA. Genetic diversity in an Andean population from Peru and regional migration patterns of Amerindians in South America: Data from Y chromosome and mitochondrial DNA. Hum Hered. 2001;51(1-2):97-106. Doi:10.1159/000022964. [ Links ]

Rodríguez JV. Los chibchas. 1st ed. Rodríguez C JV, editor. Bogotá D.C.: Editora Guadalupe Ltda; 2001. 311 p. [ Links ]

Rodríguez JV. La diversidad poblacional de Colombia en el tiempo y el espacio: Estudio craneométrico. Rev Acad Colomb Cienc . 2007;31(120):321-346. [ Links ]

Rodríguez JV. Los chibchas: hijos del sol, la luna y los Andes. Orígenes de su diversidad. Bogotá D.C.: Instituto de Desarrollo Urbano: Universidad Nacional de Colombia; 2011. 282 p. [ Links ]

Ruge L. Condiciones de vida del grupo prehispánico Chitarero. Evidencias bioarqueológicas de violencia interpersonal (Tesis de pregrado). Bogotá: Antropología, Faculdad de Ciencias Humanas, Universidad Nacional de Colombia; 2014. 95 p. [ Links ]

Sánchez C. Secuenciación de ADN mitocondrial a partir de fragmentos óseos prehispánicos hallados en el sector de Candelaria La Nueva en Bogotá (Tesis de maestria). Bogotá: Instituto de Genética, Facultad de Medicina, Pontificia Universidad Javeriana; 2007. 174 p. [ Links ]

Schurr TG, Ballinger SW, Gan YY, Hodge J, Merriwether D, Lawrence DN, et al. Amerindian mitochondrial DNAs have rare Asian mutations at high frequencies, suggesting they derived from four primary maternal lineages. Am J Hum Genet . 1990;46(3):613-623. [ Links ]

Silva A, Briceño I, Burgos J, Torres D, Villegas V, Gómez A, et al. Análisis de ADN mitocondrial en una muestra de restos óseos arcaicos del periodo Herrera en la sabana de Bogotá. Biomédica. 2008:569-577. [ Links ]

Silva E. Sobre arqueología y antropología chibcha. Rev Univ Nac. 1947;8:233-253. [ Links ]

Tamm E, Kivisild T, Reidla M, Metspalu M, Smith DG, Mulligan CJ, et al. Beringian standstill and spread of native American founders. PLoS One . 2007;2(9):1-6. Doi:10.1371/journal.pone.0000829. [ Links ]

Tarazona-Santos E, Carvalho-Silva DR, Pettener D, Luiselli D, De Stefano GF, Labarga CM, et al. Genetic differentiation in South Amerindians is related to environmental and cultural diversity: evidence from the Y chromosome. Am J Hum Genet . 2001;68(6):1485-1496. Doi:10.1086/320601. [ Links ]

Torroni A, Schurr TG, Yang CC, Szathmary EJ, Williams RC, Schanfield MS, et al. Native American mitochondrial DNA analysis indicates that the Amerind and the Nadene populations were founded by two independent migrations. Genetics. 1992;130:153-162. [ Links ]

Torroni A, Schurr TG, Cabell MF, Brown MD, Neel J V, Larsen M, et al. Asian affinities and continental radiation of the four founding Native American mtDNAs. Am J Hum Genet . 1993;53(3):563-590. [ Links ]

Uricoechea D. Análisis de ADN mitocondrial en una muestra de restos humanos prehispánicos encontrados en Jericó, Boyacá, Colombia (Tesis de pregrado). Bogotá: Departamento de Biología, Facultad de Ciencias, Universidad El Bosque; 2010. 87 p. [ Links ]

Van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009;30(2):386-394. Doi:10.1002/humu.20921. [ Links ]

Associate Editor: Graciela Bailliet.

Citation/Citar este artículo como: Casas-Vargas A, Romero LM, Rodríguez JV, Usaquén W. Análisis de ADN mitocondrial en una muestra de restos óseos precolombinos de Norte de Santander, Colombia (área cultural Chitarera). Acta biol. Colomb. 2018;23(3):263-273. DOI: http://dx.doi.org/10.15446/abc.v23n3.65407

CONFLICTOS DE INTERÉS Los autores declaramos que no tenemos conflictos de interés.

Recibido: 01 de Junio de 2017; Revisado: 08 de Agosto de 2018; Aprobado: 14 de Junio de 2018

* For correspondence. lacasasv@unal.edu.co

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons