Introducción

La genética de poblaciones se mantiene más vigente que nunca y con mayor fuerza debido a los avances que ha tenido la biología molecular en la secuenciación e identificación de genes y regiones no codificantes en los seres humanos. A través de los marcadores encontrados en el ADN y su variación, podemos dilucidar los orígenes del hombre y la dinámica que han presentado nuestras poblaciones. En el caso de Venezuela, su población se caracteriza por presentar una gran variabilidad genética y cultural como resultado de las numerosas migraciones que ha recibido a lo largo de la historia y, por lo tanto, su composición étnica constituye una mezcla entre indígenas, africanos y en mayor medida europeos. A su vez, cada región del país posee su propia estructura genética, que refleja su historia particular. En la actualidad, según los resultados del censo poblacional realizado por el Instituto Nacional de Estadísticas (INE) del año 2011, 725 141 individuos se declararon indígenas y se encuentran distribuidos a lo largo de Venezuela. En la población total indígena registrada, 50.6% son hombres (365 920) y 49.9% son mujeres (359 208); además, se reportó un aumento significativo de población indígena respecto al total de 41.8% entre 2001 y 2011 (^{Instituto Nacional de Estadística 2013}). En el estado Zulia se encuentra la mayor concentración de indígenas con 443 673 lo que representa 61.18% de la población aborigen de Venezuela.

La ^{Ley de Demarcación y Garantía del Hábitat y Tierras de los Pueblos Indígenas (2001)} indica la existencia de 35 grupos étnicos en el país y entre ellos reconoce cinco en el estado Zulia: wayuu, bari, yukpa, japreria y los añu o paraujanos quienes mantienen un arraigo sociocultural e histórico importante. Los añu o paraujanos son el segundo grupo de indígenas desde el punto de vista numérico en el estado Zulia, después de los wayuu que ocupan el primer lugar y los quintos a nivel nacional, representando 2.9% de toda la población indígena; siguiendo a los wayuu se encuentran los waraos, kariña y pemones 2011 (^{Instituto Nacional de Estadística 2013}).

Los indígenas de la etnia añu son conocidos también como la “población de agua”; habitan el noroeste del estado Zulia, tanto en palafitos como en tierra firme, desde la laguna de Sinamaica, ciénagas vecinas y el río Limón, hasta Carrasquero, Campo Mara, San Rafael del Moján e islas de la Bahía de Urubá; el barrio Santa Rosa de Agua y barrios vecinos en Maracaibo; y en la costa noroeste del Lago de Maracaibo, desde el Curarire hasta la desembocadura del Río Palmar; así mismo, se conoce que un grupo de ellos habita al sur del lago específicamente en el Municipio Colón (^{Bermúdez et al. 2009}; ^{Novelli 2012}). Los añú se autodenominan “laguneros” para reconocerse como habitantes de la laguna de Sinamaica y diferenciarse de otras personas que allí viven y que no pertenecen a esta etnia indígena. Asimismo, han presentado una alta tasa de endogamia por ser una sociedad culturalmente muy cerrada, similar a los grupos étnicos barí y yukpa. Para algunos investigadores la etnia zuliana se encuentra en peligro de extinción debido a la occidentalización, de ahí la importancia del estudio de su estructura genética y establecer cuánto pueden semejarse o diferenciarse de otras etnias de la región y obtener información acerca de su mestizaje (^{Amodio et al. 2005}). En este sentido se puede resaltar desde el punto de vista genético la importancia de conocer la estructura de estos pueblos y compararlas con las ya estudiadas en el continente americano, para conocer sus orígenes y, de igual manera, identificar el grado de mezcla que pudo haberse producido tras su incorporación cada vez mayor a las ciudades. Lo dicho anteriormente es de gran relevancia en el estado Zulia donde se encuentra concentrada la mayor parte de la población indígena del país; y más aún la comunidad añu o paraujana elegida para este estudio, que representa el quinto grupo en frecuencia después de los wayuu, los waraos, los kariñas y los pemones, siendo además la segunda población indígena con mayor frecuencia en el estado.

Materiales y métodos

Resultados y discusión

La población añú se encuentra en equilibrio de Hardy-Weinberg (HW) para 14 de los marcadores autosómicos estudiados, lo que permite inferir que es una población en donde los individuos tienen un apareamiento aleatorio y que los eventos migratorios que han podido haber marcado una dinámica poblacional, han tenido poco efecto en la variación de las proporciones genotípicas de la misma. Sin embargo, se observa que el marcador D8S1179 se encuentra en desequilibrio de Hardy-Weinberg (P < 0,003) puesto que la Ho es significativamente inferior a la He. El análisis de Hardy-Weinberg busca establecer si en una población las frecuencias génicas y genotípicas permanecen constantes de generación en generación, cuando la población cumple con apareamiento aleatorio, herencia mendeliana y ausencia de selección, mutación o migración (^{Falconer y Mackay 1996}). En un porcentaje elevado de estos loci, el desvió podría explicarse a causa de un déficit de heterocigotos, presumiblemente dado por un mayor nivel de endogamia, ya que los loci con mayor frecuencia se desvían del equilibrio HW.

De los 15 marcadores autosómicos estudiados, catorce presentaron heterocigosis mayores a 58% con excepción de D7S820 con 48%, siendo los loci PENTA E y D13S317 los que presentaron la mayor heterocigosis con 89.6%. En relación con la variabilidad de los marcadores en la población añú, de los 15 marcadores STR’S analizados se observó que los marcadores PENTA E, D18S51, FGA y PENTA D fueron los que presentaban una mayor variabilidad en cuanto al número de alelos. PENTA E resultó ser el más polimórfico en esta investigación, al tener presente 16 alelos que van desde el 7 hasta el 21, el de mayor frecuencia fue el alelo 14 con 24.1%, seguido del marcador FGA, el cual presentó diez alelos que van desde el 19 al 28 siendo el de mayor frecuencia el alelo 24 con 39.6%; D18S51 presentó 9 alelos del 12 al 23, siendo el de mayor frecuencia el 12 con 25.8%, y penta D mostró 8 alelos distintos desde el 8 al 15 siendo el de mayor frecuencia el alelo 9 con 32.8%.

En cambio, los marcadores TH01 y D7S820 fueron los que presentaron menos variabilidad con 5 y 4 alelos respectivamente. El marcador con mayor frecuencia alélica independientemente del número de alelos presentado en la población añú fue en primer lugar el alelo 16 de VWA con 62%, en segundo lugar el alelo 11 de D5S818 con 60% y, en tercer lugar los alelos 6 de TH01 y 11 de D7S820, ambos con una frecuencia de 58%. No se obtuvo presencia de alelos raros (frecuencia alélica < 0,005%) para la población indígena añu (cuadro 1).

Pineda y colaboradores estudiaron 15 marcadores STR’S autosómicos en la población de Maracaibo y encontraron que el marcador con mayor variabilidad fue FGA con 14 alelos siendo los más frecuente los alelos 21, 23 y 24 todos con 16%, y los de menor variabilidad D3S1358 y TH01 ambos con 7 alelos, mostrando el alelo 15 con la mayor frecuencia; para D3S1358 con 30% y el alelo 6 del marcador TH01 también con 30% (Pineda et al. 2006).

Al comparar la variabilidad de los alelos encontrados en nuestro estudio con el llevado a cabo por Filho Rigotti en la población brasileña, observamos que el marcador con mayor variabilidad en cuanto al número de alelos fue D21S11, el cual presentó 21 alelos distintos, siendo el de mayor frecuencia el alelo 30, el FGA presentaba 20 alelos siendo el de mayor frecuencia el alelo 21 y el D18S51 presentaba 16 alelos del 9 al 23 siendo el más frecuente el alelo 11, en cambio, los marcadores D7S820, CSF1PO, D16S539, TPOX y el D5S818 presentaban menor variabilidad al mostrar en todos los casos observados 9 alelos (^{Filho 2012}). Por su parte, en un estudio realizado por Ossa y colaboradores en la población de Antioquía, Colombia, de 14 marcadores STR’S autosómicos, se encontró que el sistema PENTA E fue el marcador que presentó un mayor número de alelos entre los individuos con 16, siendo el de mayor frecuencia el alelo 11 con 19% y D3S1358 el marcador que presentó un menor número de alelos con 5, entre los cuales, el más frecuente fue el alelo 14 (^{Ossa et al. 2010}).

En otro estudio llevado a cabo por Ramírez en la Ciudad de México, mediante el uso de 22 loci STR’S autosómicos, se determinó que el loci con mayor variabilidad fue el PENTA E con 23 alelos siendo el más frecuente el 12 con 16%, D18S51 con 15 y el de mayor frecuencia el alelo 17 con 18% y FGA con 14, lo cual mostró que el de mayor frecuencia fue el alelo 24 con 16%; por el contrario, el de menor variabilidad fue TPOX con 5 alelos y el más frecuente fue el alelo 8 (^{Ramírez 2014}).

Del mismo modo, en un estudio llevado a cabo por Flores y colaboradores en la población de Belice, se determinaron las frecuencias alélicas de 15 marcadores STR’S autosómicos y encontraron la mayor variabilidad para el locus D21S11 con 17 alelos siendo el más frecuente el 30 con 22% y el locus D18S51 también con 17, donde el alelo 17 fue el más frecuente para este marcador con 21%; los que presentaron menos variabilidad en esta población fueron CSF1PO, D3S1358, TH01 y D16S539 todos con 8 alelos (^{Flores et al. 2015}).

En el caso de una investigación realizada por Hernández y Trejo en la ciudad de Zacatecas, México, mediante el empleo de 15 marcadores STR’S autosómicos se encontró que los marcadores con mayor variabilidad fueron D18S51 y D21S11 ambos con 15 alelos, siendo más frecuente el alelo 14 con 18% para D18S51 y el alelo 30 con 27% para D21S11 y el de menor variabilidad fue TH01 con 6 alelos, siendo el alelo 7 el más frecuente con 33% (^{Hernández-Rodríguez y Trejo-Medinilla 2014}).

Con esto podemos observar que los marcadores PENTA E, D18S51 y FGA, son los que presentan mayor variabilidad en cuanto al número de alelos, tanto en la población añú como en otras. Por su parte, la elaboración de la matriz de distancias genéticas entre poblaciones se estimó con base en ^{Nei (1974)} tomando en cuenta las frecuencias alélicas de 7 marcadores de la población añú: CSF1PO, TPOX, TH01, vWA, D7S820, D13S317 y D5S818 y las reportadas para diversas poblaciones del estado Zulia tales como las poblaciones indígenas barí y yukpa, la población mestiza de Maracaibo, la población afro-descendiente de San José de Heras e Isla de Toas, conjuntamente con otras poblaciones de Latinoamérica, como los indígenas waorani de Ecuador, mestizos de Zacatecas, México, mestizos de Antioquía, Colombia, mestizos del sur de Brasil y de Buenos Aires, Argentina. Los resultados de la misma pueden observarse en el cuadro 2.

Las variaciones biológicas entre poblaciones humanas pueden ser investigadas a través de distintas herramientas metodológicas, sin embargo la bibliografía internacional recomienda el uso de las distancias genéticas y la representación gráfica de las matrices a través de dendrogramas o mapas genéticos que son usados frecuentemente para hacer inferencias con respecto a patrones geográficos, lingüísticos o etnohistóricos (^{Demarchi 2009}). Existen diferentes tipos de distancias genéticas que se pueden aplicar para estudios de diferenciación génica de las cuales la más utilizada para polimorfismos del tipo microsatélite es la distancia genética del tipo DA de ^{Saitou y Nei (1987)}. Ésta ha demostrado tener ciertas ventajas sobre las otras, ya que es más discriminativa entre poblaciones cercanamente relacionadas como las humanas, y aunque no es lineal con el tiempo evolucionario, se observa que es más eficiente para obtener correctas tipologías a nivel de los árboles, por lo cual se seleccionó para esta investigación (^{Zabala et al. 2005})

Según Cavalli-Sforza, cuando las distancias genéticas son estudiadas bajo la construcción de un árbol filogenético, las poblaciones tenderán a unirse o agruparse con aquellas a las que tengan mayor afinidad. En el caso de poblaciones mezcladas, formarán grupos con las poblaciones que hayan aportado el mayor componente a su acervo genético (^{Cavalli-Sforza 2007}).

De acuerdo con hipótesis teóricas y algunos resultados genéticos, los pueblos indígenas de América provienen de Asia, entraron a América aproximadamente hace 18 000 años al atravesar el estrecho de Bering; a medida que descendían por el continente se iban asentando en diferentes grupos. El territorio venezolano fue alcanzado por estos pobladores aproximadamente hace 15 000 años en la época Paleoindia; se asentaron en la zona noroccidental del país y se caracterizaban por ser recolectores, cazadores y pescadores (^{Sans 2000}; ^{Leal-Jiménez 2008}). Al considerar lo anteriormente expuesto, sería lógico plantear que los pueblos indígenas americanos comparten un mapa genético muy semejante, sobre todo aquellos que se encuentran geográficamente cercanos, como los indígenas añú, barí y yukpa que están ubicados en el estado Zulia, sin embargo, en este estudio queda demostrado que no es así.

Cuando analizamos el dendrograma construido al medirse por frecuencias alélicas de los microsatelites estudiados: CSF1P0, TPOX, TH01, vWA, D7S820, D13S317 y D5S818, encontramos que la población indígena añú se ubicó en una rama totalmente aislada, lo cual nos indica que no comparten su acervo genético con otros grupos indígenas zulianos, ni con estas poblaciones, por lo que se podría suponer que este grupo étnico tiene su propio origen ancestral, además de poco o ningún mestizaje (figura 1).

Figura 1 Dendrograma construido con el método de Distancia Da Nei con NJ a partir de los microsatelites estudiados: CSF1P0, TPOX, TH01, vWA, D7S820, D13S317 y D5S818 en la población añú y diversas poblaciones americanas. 

Los grupos indígenas barí, yukpa y añú pertenecen a familias etnolingüísticas diferentes. En el caso de los barí provienen del tronco de los chibchas, los yukpas provienen de la rama norte de la familia lingüística caribe y los añú de la familia lingüística arawak, por lo que no solo existe un distanciamiento geográfico sino cultural que contribuye a que sean poblaciones cerradas y genéticamente diferentes. A continuación se encuentran ubicadas las poblaciones de América Central y del Sur, las cuales, a pesar de tener historias y culturas diferentes, tienen características comunes. En primer lugar, la población de Zacatecas (México), podría explicar su mayor cercanía con las poblaciones indígenas debido a que se trata de un país que tenía una alta densidad de indios americanos en la época precolombina y que, a pesar de haber pasado también por procesos colonizadores y muchas migraciones, aún mantienen una fuerte contribución amerindia, mostrando un grado más alto de linaje americano nativo. A continuación se encuentra la población del sur (Brasil) y del departamento de Antioquía, Colombia, las cuales a su vez se agruparon por su cercanía genética. Luego encontramos la población de Buenos Aires (Argentina). Por último y más alejadas genéticamente de la población indígena añú, a pesar de su cercanía geográfica, se encuentran las poblaciones del estado Zulia como Maracaibo y San José de Heras e Isla de Toas, quienes poseen un alto componente afrodescendiente.

Algunos investigadores encuentran correlación entre genes y lenguaje; contrario a los resultados que se encontraron en este trabajo, en el dendrograma NJ se muestra que la etnia yukpa y barí del estado Zulia y waoranis de Ecuador están próximos genéticamente a pesar de que sus lenguas son radicalmente diferentes. Los yukpa son de la filiación lingüística caribe, los barí de la familia lingüística chibcha y waoranis de Ecuador cuyo idioma es el huao terero, no derivado de ninguna de estas ramas lingüísticas. Por su parte, los indígenas añú se encuentran en una rama separada, lo que demuestra que están distantes de estos grupos desde el punto de vista genético y del lingüístico ya que originalmente éstos hablaban un idioma propio, “el añú”, de la filiación lingüística arawaca, que es la más extendida e importante familia lingüística de América del sur. Actualmente, es bien sabido que ha cesado la transmisión intergeneracional del idioma, pues han adoptado en su mayoría el castellano para comunicarse.

Otro factor que podemos considerar al observar la diferenciación de la población añú con respecto a otras estudiadas, es su ubicación geográfica. La laguna de Sinamaica está ubicada al noroeste del estado Zulia, una zona de difícil acceso incluso en la actualidad, por lo que se ha dificultado el flujo genético entre poblaciones. La cultura es un elemento de gran influencia pues la población añú es una comunidad cerrada cuya estructura social es muy importante y siempre buscan mantener la identidad de su pueblo, por lo tanto, prevalece la unión consanguínea entre primos.

Por todo lo expuesto anteriormente, podemos suponer que condiciones como las barreras lingüísticas, la distancia geográfica y los patrones culturales han traído como consecuencia el aislamiento y cierre de la población añú, evitando así el intercambio genético con otros grupos étnicos.

Estos resultados, muestran la gran utilidad que tienen los marcadores moleculares para establecer las diferencias genéticas entre grupos que cohabitan dentro de una misma área geográfica a la par que proveen información que apoya la evolución demográfica y destaca la importancia de disponer de bases de datos locales para utilidad de la genética forense y la antropología molecular.

Conclusiones

Éste es el primer reporte acerca de la estructura genética de la población indígena añú de la laguna de Sinamaica; aporta datos importantes que contribuyen al conocimiento del Acervo Genético de la región zuliana. Las frecuencias alélicas obtenidas son importantes por el establecimiento de bases de datos locales para su aplicación en estudios de predisposición a diversas enfermedades, investigaciones en el área de la Farmacogenética, entre otros aplicados a la población indígena añú. A pesar de que los indígenas añú constituyen un grupo poblacional cerrado, el alto nivel de heterocigosis demostrado en algunos marcadores, favorece los estudios de genética forense. Con base en los resultados de las distancias genéticas obtenidas se pudo demostrar que la etnia indígena añú posee identidad genética propia que permite diferenciarla del resto de las poblaciones estudiadas.