INTRODUCCIÓN

El frijol ayocote (Phaseolus coccineus L.) es una leguminosa que tiene su centro de origen y domesticación en México (^{Llaca et al., 1994}; ^{Sicard et al., 2005}) el cual, de acuerdo con ^{Salinas (1988)} se ubica en las zonas altas de los estados de Puebla, Oaxaca, y Chiapas. Sin embargo, el frijol ayocote en México es poco cultivado y consumido no así en regiones de Europa principalmente España, Holanda y Reino Unido, donde su consumo ha sustituido al frijol común (P. vulgaris L.) y se consume principalmente en su presentación de ejote y de grano seco en el sur de Italia (^{Santalla et al., 2004}; ^{Rodiño et al., 2007}). Por otra parte, se ha reportado su consumo en Japón donde se cultiva principalmente en zonas montañosas y su semilla es muy apreciada por su excelente calidad culinaria (^{Watanabe & Tojo, 2006}). De las especies domesticadas del género Phaseolus spp., el frijol ayocote es la segunda especie en importancia económica (^{Saburido-Álvarez & Herrera-Estrella, 2015}). El cultivo tradicional ocurre en pequeña escala (^{Schwember et al., 2017}). En Centroamérica, el frijol ayocote se cultiva casi exclusivamente en asociación con maíz (Zea mays L.) (^{Vargas-Vázquez et al., 2007}). Dado que los frijoles son un elemento importante en el autoabasto de poblaciones rurales, además de ser la leguminosa de grano más importante constituye la fuente principal de proteína de calidad a bajo costo, con diversas variaciones en cuanto a formas de cultivo, uso y maneras de preparación para consumo humano (^{Svetleva et al., 2003}).

El frijol ayocote en México se siembra con variedades nativas o criollas y su cultivo se ha reducido de manera alarmante en los últimos 50 años (^{Vargas-Vázquez et al., 2007}). Además, el manejo agronómico de esta especie es diferente al del frijol común porque presenta aproximadamente 14.7% de cruzamiento natural (^{Muruaga et al., 1992}), lo que dificulta conservar la pureza genética de las variedades criollas. Lo anterior, aunado al manejo en la selección e intercambio de semilla por parte de los productores, impacta directamente en la estructura poblacional de esta especie (^{Mercati et al., 2015}). Sin embargo, se han conducido pocos trabajos que comparen el germoplasma de P. coccineus y de P. vulgaris con respecto a los rasgos morfoagronómicos y genéticos útiles para el mejoramiento (^{Acampora et al., 2007}).

La diversidad genética del frijol ayocote se ha estudiado previamente aplicando diversas estrategias (^{Nowosielski et al., 2002}; ^{Sicard et al., 2005}; ^{Boczkowska et al., 2012}), principalmente marcadores moleculares de ADN y marcadores morfoagronómicos. Por su parte, ^{Spataro et al. (2011)} desarrollaron una colección núcleo de frijol ayocote con 228 accesiones (52 de Mesoamérica y 148 de Europa) y se ha confirmado la diferenciación genética entre ambos acervos genéticos (^{Rodriguez et al., 2013}; ^{Mercati et al., 2015}) y diferenciación entre poblaciones debido a su estructuración genética (^{Catarcione et al., 2013}; ^{Mercati et al., 2015}). En Mexico también se estructuró una colección nacional de P. coccineus que incluye 798 accesiones, y donde alrededor del 80% proviene de las regiones del eje Neovolcánico y la Sierra Madre Oriental (^{Vargas-Vázquez et al., 2012}). El frijol ayocote del Carso Huasteco se subdivide en dos grupos: uno que incluye germoplasma con ciclo de vida tardío y regiones con temperaturas mínimas de 2-5°C con semillas y vainas grandes y otro grupo con genotipos precoces que se colectaron en sitios con temperaturas mínimas de 0 a -2°C, así como semillas y vainas pequeñas (^{Ruíz-Salazar et al., 2010}; ^{Vargas-Vázquez et al., 2011}; ²⁰¹³; ²⁰¹⁴).

Por lo anterior, el objetivo de la presente investigación fue determinar, por primera vez, la variabilidad y la estructura genética de germoplasma de frijol ayocote de la Subprovincia Carso Huasteco de Puebla, México, mediante marcadores moleculares AFLP.

MATERIALES Y MÉTODOS

RESULTADOS

La amplificación de cuatro combinaciones de oligonucleótidos AFLP arrojó 256 bandas, de las cuales 224 fueron polimórficas (87.4%) y 32 monomórficas (12.6%) con estos datos se construyó la matriz de presencia/ausencia de bandas (cuadro 2). El AMOVA presentó un valor F_ST = 0.067; además, se observó que la mayor proporción de la varianza molecular explicada se encuentra dentro de las accesiones (93.3%) (cuadro 3). El análisis de conglomerados con datos AFLP formó dos grupos principales (A y B) sin embargo no se observaron agrupamientos definidos con base en el origen del germoplasma pues los testigos Pinto Villa, Pinto Zapata (P. vulgaris) L. y P. glabellus L. (grupo B) fueron distintos a P. coccineus (grupo A), donde se observó una mezcla de poblaciones (fig. 1).

Fig. 1 Dendrograma de 122 accesiones de frijol ayocote mediante el metodo UPGMA. 

El análisis de la estructura poblacional de accesiones de frijol ayocote determinó un valor de delta K = 3. El primer grupo (K1) incluye accesiones de once localidades: Tlatlauquitepec, M. Serdán, Zacapoaxtla, Nauzontla, M. Cuetzalan, M. Zacapoaxtla, Ahuacatlán, Xochiapulco, Atempan y Zoquiapan y los testigos; el segundo conglomerado (K2) incluye germoplasma de diez localidades (Nauzontla, Zacapoaxtla, Atempan, M. Zacapoaxtla, M. Tlatlauquitepec, Zoquiapan, Chignahuapan, Tlatlauquitepec, Zacatlán y Tételes de Ávila Castillo). Finalmente, el tercer grupo (K3) muestra genotipos de 12 localidades: Tlatlauquitepec, M. Serdán, Zacapoaxtla, Zacatlán, Atempan, Zoquiapan, Chignahuapan, Huauchinango, M. Cuetzalan, M. Tlatlauquitepec, Nauzontla y Ahuacatlán. Lo anterior confirmó los resultados del análisis de conglomerados dado que no se observaron agrupamientos definidos con base en el origen de colecta de cada accesión. Se realizó un segundo análisis donde se consideró un valor de delta K = 4 y se observó que el grupo testigo que incluyó dos accesiones de de P. vulgaris (Pinto Villa y P. Zapata) y una de P. glabellus se agruparon de manera conjunta en el clóster cuatro (K4) (fig. 2).

Fig. 2 Determinación de la estructura poblacional mediante AFLP de 122 accesiones de Phaseolus coccineus de la Subprovincia Carso Huasteco. 

DISCUSIÓN

En este trabajo se llevó a cabo el análisis de la variabilidad y la estructura genética de accesiones de P. coccineus provenientes del Carso Huasteco de México, una de las regiones más importantes de origen y diversidad de frijol ayocote en Mesoamérica, aplicando la estrategia de marcadores moleculares tipo AFLP. La estrategia molecular fue eficiente para el análisis genético en virtud de que detectó más del 87% de productos polimórficos en el germoplasma analizado y a su condición de estrategia de 'alta eficiencia multiplex al generar alta cantidad de productos amplificados por ensayo o reacción de amplificación (^{Bitocchi et al., 2015}).

El análisis molecular separó claramente al germoplasma de P. coccineus de los testigos pertenecientes a las especies P. vulgaris y P. glabellus (^{Nowosielski et al., 2002}); aunque la separación entre el germoplasma del Carso Huasteco y P. coccineus tipo y la variedad de frijol ayocote Blanco Tlaxcala ya no fue tan clara. Sin embargo, la formación de al menos tres sub-grupos de genotipos de frijol ayocote con diferentes orígenes demuestra la existencia de poblaciones genéticamente similares y adaptadas a rangos climáticos específicos y en constante intercambio debido a la 'captura' de nuevos alelos dada la posible polinización cruzada que ocurre naturalmente en sus zonas de distribución (^{Spataro et al., 2011}; ^{Bitocchi et al., 2015}).

Éste es el primer trabajo en México que analiza genéticamente P. coccineus con base en una estrategia molecular altamente confiable e informativa y donde, además, se incluyen poblaciones representativas de la especie. Los estudios anteriores han enfatizado el estudio de P. coccineus con base en rasgos fisiológicos, fenológicos, morfológicos o productivos (^{Vargas-Vázquez et al., 2011}; ²⁰¹³; ²⁰¹⁴). Destaca el hecho de encontrar la mayor proporción de la varianza genética molecular dentro de las accesiones más que entre accesiones o grupos de accesiones con base en el origen (^{Spataro et al., 2011}) mientras que en frijol común se reportan proporciones diferentes (^{Gill-Langarica et al., 2011}). Lo anterior sugiere genotipos de frijol ayocote diversos. Además, esto indica que, aunque el análisis genético de plantas basado en la formación de 'bulks' de DNA reduce el trabajo y costos, puede subestimar la real naturaleza genética de cada accesión o población de genotipos de frijol. Las futuras estrategias de conservación del germoplasma de frijol ayocote deben considerar mecanismos de multi plicación que garanticen preservar la variabilidad genética de cada genotipo.

El análisis genético de P. coccineus ha puesto énfasis en germoplasma de Europa y los resultados coinciden en reportar diversidad genética reducida (^{Sicard et al., 2005}; ^{Boczkowska et al., 2012}). Incluso, se han podido diferenciar genotipos de frijol ayocote con base en acervos genéticos (^{Sicard et al., 2005}; ^{Acampora et al., 2007}). El análisis de una colección 'mundial' de P. coccineus (228 genotipos de Mesoamérica -52-, Europa -148- y silvestres -28-) volvió a recalcar la diversidad genética reducida del frijol ayocote (^{Spataro et al., 2011}). En este sentido se enfatiza que el germoplasma analizado fue introducido a Europa proveniente de América (^{Rodriguez et al., 2013}; ^{Mercati et al., 2015}) sufrió el efecto del 'cuello de botella' que precedió a la expansión de la especie en el continente Europeo (^{Rodríguez et al., 2013}).

Sin embargo, ^{Rodríguez et al. (2013)} sugieren la diferenciación genética adaptativa que luego condujo a la diversificación de P. coccineus en Europa, con la consiguiente estructuración genética de poblaciones (^{Catarcione et al., 2013}; (^{Mercati et al., 2015})en virtud a factores tales como la polinización abierta (≈14% según ^{Muruaga et al., 1992}) y el intercambio de germoplasma (semilla) y la selección empírica de genotipos por los agricultores locales (^{Castillo-Mendoza et al., 2006}; ^{Vargas-Vázquez et al., 2011}; ²⁰¹²; ²⁰¹³; ²⁰¹⁴); Así, los estudios previos concluyen que Europa debe considerarse un centro de diversificación secundario de P. coccineus (^{Mercati et al., 2015}).

La variabilidad genética de P. coccineus podrá ser aprovechada para el mejoramiento genético de la propia especie o de otras cercanas como P. vulgaris, en la medida que se intensifiquen los esfuerzos en caracterizar el germoplasma disponible de la especie y se identifiquen así genotipos que puedan servir como fuentes de resistencia a factores como las enfermedades y el frío. También, se propone apoyarse en tecnologías de secuenciación de nueva generación para el análisis de la diversidad genética, el mapeo genético y estudios de sintenia para identificar SNPs o genes funcionales e 'indels' asociados con características agronómicas importantes (^{Schwember et al., 2017}).

El análisis con 42 548 SNPs generados por secuenciación de 242 individuos de P. coccineus y las formas domesticadas de especies cercanas a P. vulgaris (20) y P. dumosus (35) indicó la formación de ocho conglomerados o grupos de genotipos, donde la mitad correspondió a poblaciones silvestres y el resto, a domesticadas. Dado que las poblaciones domesticadas formaron un solo clado monofilético se asume que P. coccineus de Méxio sólo sufrió un evento de domesticación mismo que posiblemente ocurrió en germoplasma del Cinturón Volcánico Trans-Mexicano. Contrario a lo observado previamente, no se observaron evidencias de flujo de genes recientes entre poblaciones silvestres y cultivadas. Finalmente, se identificaron 24 SNPs relacionados con domesticación, 12 con diversificación genética y ocho con selección natural. De los pocos SNPs identificados y que están anotados, se expresan en flores y vainas (^{Guerra-García et al., 2017}).

El germoplasma silvestre de P. coccineus se dividió en cuatro grupos genéticos que exhiben diferenciación genética significativa (^{Guerra-García et al., 2017}). Así como la diversidad contenida en los parientes silvestres es un determinante importante en la evolución de cultivos, el de las poblaciónes domesticadas también lo es. La domesticación reduce la diversidad genética (^{Schmutz et al., 2014}) debido a la deriva genética resultante de los 'cuellos de botella' poblacionales y la selección artificial (^{Schmutz et al., 2014}); además de los tamaños efectivos de las poblaciones o la alta tasa de polinización cruzada (^{Muruaga et al., 1992}).

CONCLUSIONES

El análisis genético del germoplasma de P. coccineus del Carso Huasteco de México indicó la formación de tres grupos de accesiones, aunque no hubo relación con el origen del germoplasma; así como mayor varianza genética molecular dentro de accesiones que entre ellas.