Introducción
Las especies Theobroma cacao L. y T. bicolor Humb. & Bonpl. pertenecen a la familia Sterculiaceae, alternativamente Malvaceae (Southgate, 2000; Zequeira y Beauregard, 2024). Ambas especies son árboles con hojas siempre verdes que alcanzan altura de 3 - 12 m, y se desarrollan en la franja geográfica entre los 18( LN y 15( LS, con temperaturas de 20º a 30º C y más de 2,000 mm de lluvia distribuida en todo el año (Cuatrecasas, 1964). Por otro lado, Southgate (2000) menciona que T. cacao se cultivó por primera vez durante los años 250 - 900 de la era actual, aunque Cornejo et al., (2018) señalan que el cacao fue domesticado hace 3,600 años de donde surgió la variedad ‘Criollo’ y posteriormente fue introducido a Europa a mediados de 1550 (Dillinger et al., 2000; Coe & Coe, 2018). El cacao es una planta cuyos frutos tienen gran capacidad antioxidante debido a los polifenoles que contiene (Kim et al., 2011; Córdova-Avalos et al., 2023; Gómez-Rivera et al., 2024). A nivel mundial se reconocen los tipos de cacao Criollo, Forastero y Trinitario, este último es una variedad híbrida resultado del cruzamiento de las dos primeras. En México el cacao criollo es apreciado por su perfil sensorial y también es un cultivo importante para la economía familiar del sureste de la república Mexicana (Córdova-Lázaro et al., 2018; Pérez-Obrador et al., 2025). Con respecto al centro de origen del cacao, Zarrillo et al., (2018) mencionan que el cacao es originario de la zona amazónica de Santa Ana-La Florida al sureste de Ecuador y refuerzan su hipótesis con el descubrimiento de semillas de 5,300 años de antigüedad y el correspondiente análisis molecular de DNA (Pérez-Obrador et al., 2025).
Respecto al pataxte (T. bicolor), éste es una especie con alto potencial productivo que se encuentra en proceso de domesticación y su distribución se restringe al Continente Americano (Tinajero-Carrizales et al., 2021). En México, al pataxte se le encuentra en los estados de Chiapas, Oaxaca, Tabasco y con pequeñas poblaciones en el estado de Guerrero. Sin embargo, también crece y se desarrolla en Colombia, Venezuela y Brasil, de manera general al pataxte se le puede encontrar asociado a cacao o como cultivo de traspatio (Rendón-Aguilar et al., 1998; Gálvez-Marroquín et al., 2016).
La diversidad genética y estructura poblacional del cacao ha sido estudiada con diferentes tipos de marcadores moleculares de ADN como los polimorfimos de secuencia simple (SNP) (Martínez et al., 2017; Osorio-Guarín et al., 2017; Cornejo et al., 2018; Lindo et al., 2018; De Wever et al., 2019; Bustamante et al., 2022), polimorfismo en longitud de los fragmentos de restricción (RFLP) (Motamayor & Lanaud, 2002; Motamayor et al., 2003), secuencias simples repetidas (SSRs) (Zhang et al., 2006; Ruíz-Erazo et al., 2015) y con marcadores de cloroplasto (Nieves-Orduña et al., 2021). Sin embargo, hay pocos estudios en los que se haya realizado análisis de la diversidad genética del pataxte con este tipo de marcadores, dentro de ellos se puede mencionar el uso de los RAPDs o amplificación aleatoria de ADN polimórfico (Moreno et al., 2004). Arteaga-Voigt et al. (2016) estudiaron con nueve marcadores SSRs la estructura genética y relación filogenética de poblaciones de Cacao Nacional Silvestre y Cultivado procedente de 19 zonas productoras del Norte de La Paz, Bolivia. De sus resultados los autores citados reportan que el número de polimorfismo por locus fue de 4 a 13 alelos, el contenido de información polimórfica (PIC), presentó patrones de distribución semejantes entre sí. En todas las poblaciones estudiadas la heterocigosidad esperada (He) y heterocigosidad observada (Ho) presentó un 99.98%. Por otro lado, Ricaño-Rodríguez et al. (2019) genotipificaron muestras de cacao de Tabasco y Chiapas, y encontraron que T. cacao posee porcentajes entre 2-11% de SNPs. Los autores referidos observaron secuencias conservadas en las regiones codificantes de los genomas de las variedades de cacao estudiadas. Por su parte, Pérez-Pérez et al. (2025) realizaron un estudio morfológico y molecular en pataxte de los cuatro estados en los que se le encuentra a esta especie en México, de sus principales resultados reportan no haber encontrado variación genética en las regiones de ADN analizadas, pero sí diferencias morfológicas, las que deben ser modificadas por los diferentes ambientes en dónde se tienen establecidas las plantaciones de pataxte referidas. Por la relevancia e importancia que tiene el cultivo de cacao y pataxte en el sureste de la república Mexicana, el objetivo del presente estudio fue determinar y caracterizar la estructura poblacional de 56 accesiones de Theobroma spp. mediante marcadores microsatélites.
Materiales y métodos
Material vegetal
La colecta del material experimental de esta investigación se realizó durante septiembre - diciembre 2018 y enero - agosto 2019 en los estados de Tabasco y Chiapas, México, se obtuvieron 56 accesiones, de las cuales, 43 pertenecen a cacao (T. cacao L.) representadas por cinco grupos genéticos: Guayaquil, Criollo, Calabacillo, Patastillo y Guayaquil/rojo y uno de pataxte (T. bicolor) representado por 13 colectas (Cuadro 1). Se tomaron de 2 a 3 mazorcas y de una a cinco plantas por variedad en manera aleatoria. La edad de los árboles muestreados oscila entre 15 y 25 años. Los suelos son profundos, temperatura apta para el desarrollo de los árboles y precipitación mayor a 2,000 mm durante el año.
Cuadro 1 Datos de pasaporte de 56 colectas que incluyen 43 accesiones de cacao y 13 de pataxte procedentes de Tabasco y Chiapas, MéxicoTable 1. Passport data for 56 collections, including 43 cacao accessions and 13 pataxte accessions from Tabasco and Chiapas, Mexico
| No. | Variedades | Estado | Municipio | Localidad | Tipo de tejido | LN | LO |
|---|---|---|---|---|---|---|---|
| 1 | Guayaquil | Tabasco | Centro | González 1ͣ | Fresco | 17˚10’05 | 93˚04’58” |
| 2 | Guayaquil | Chiapas | Acapetahua | Ejido Arenal | Fresco | 17˚58’38.21 | 92˚18’18.0” |
| 3 | Criollo | Chiapas | Acapetahua | Ejido Arenal | Fresco | 17˚58’38.21 | 92˚18’18.0” |
| 4 | Guayaquil | Chiapas | Chicomucelo | Zacoalpa | Fresco | 15˚10’23.27” | 92˚41’59.93” |
| 5 | Guayaquil | Tabasco | Paraíso | Oriente 1ͣ | Fresco | 18˚21’12.2” | 93˚11’54.6” |
| 6 | Calabacillo | Chiapas | Acapetahua | Ejido Arenal | Fresco | 17˚58’38.21 | 92˚18’18.0” |
| 7 | Criollo | Tabasco | Paraíso | Oriente 1ͣ | Fresco | 18˚21’12.2” | 93˚11’54.6” |
| 8 | Guayaquil | Chiapas | Chicomucelo | Piedra Labrada | Fresco | 15˚10’23.27” | 92˚41’59.93” |
| 9 | Calabacillo | Tabasco | Paraíso | Oriente 1ͣ | Fresco | 18˚21’12.2” | 93˚11’54.6” |
| 10 | Guayaquil | Tabasco | Teapa | Cerro de Puyacatengo | Fresco | 18˚07’27.7” | 92˚55’86” |
| 11 | Guayaquil | Tabasco | Cunduacán | Huapacal | Fresco | 18˚3’2.52” | 93˚9’21.24” |
| 12 | Calabacillo | Tabasco | Teapa | Cerro de Puyacatengo | Fresco | 18˚07’27.7” | 92˚55’86” |
| 13 | Patastillo | Tabasco | Cunduacán | Las Cruces | Fresco | 18˚07’36.5” | 93˚11’29” |
| 14 | Calabacillo | Tabasco | Cunduacán | Rio Seco | Fresco | 18˚07’39.1” | 93˚18’02.9” |
| 15 | Patastillo | Tabasco | Cunduacán | Rio Seco | Fresco | 18˚07’39.1” | 93˚18’02.9” |
| 16 | Calabacillo | Chiapas | Reforma | Macayo | Fresco | 17˚8’20.6” | 93˚16’27.2” |
| 17 | Calabacillo | Tabasco | Cunduacán | Los Cerros | Seco | 18˚02’14.7” | 93˚06’04.9” |
| 18 | Pataxte | Tabasco | Huimanguillo | Los Naranjos | Seco | 17˚54’05” | 93˚35’11” |
| 19 | Guayaquil | Tabasco | Cunduacán | Rio Seco | Fresco | 18˚07’36.5” | 93˚11’29” |
| 20 | Patastillo | Tabasco | Cunduacán | Pechucalco | Fresco | 18˚06’42.2” | 93˚11˚13.6” |
| 21 | Patastillo | Tabasco | Tacotalpa | Zunup y Patastal | Fresco | 17˚28’22.1” | 92˚48’48.1” |
| 22 | Calabacillo | Tabasco | Cunduacán | La Trinidad | Fresco | 18˚09’12” | 93˚16’11.7” |
| 23 | Pataxte | Tabasco | Pichucalco | Los Cascabeles | Seco | 17˚50’16” | 93˚11˴33” |
| 24 | Patastillo | Tabasco | Cunduacán | Los Cerros | Fresco | 18˚02’14.7” | 93˚06’04.9” |
| 25 | Pataxte | Tabasco | Cunduacán | Las Cruces | Seco | 18˚07’36.5” | 93˚11’29” |
| 26 | Pataxte | Tabasco | Teapa | Cerro de Puyacatengo | Fresco | 15˚41’12.26” | 92˚18’40.08” |
| 27 | Patastillo | Chiapas | Macayo | Reforma | Fresco | 17˚58’20.6” | 93˚16’27.2” |
| 28 | Calabacillo | Tabasco | Cunduacán | Las Cruces | Fresco | 18˚07’36.5” | 93˚11’29” |
| 29 | Criollo | Tabasco | Teapa | Cerro de Puyacatengo | Fresco | 17˚32’45.1” | 92˚55’04.6” |
| 30 | Patastillo | Tabasco | Tacotalpa | Agua Escondida | Fresco | 17˚29’14.3” | 92˚46’37.5” |
| 31 | Calabacillo | Tabasco | Tacotalpa | Zunup y Patastal | Fresco | 17˚28’22.1” | 92˚46’37.5” |
| 32 | Calabacillo | Tabasco | Tacotalpa | Agua Escondida | Fresco | 17˚29’14.3” | 92˚46’37.5” |
| 33 | Patastillo | Tabasco | Paraíso | Oriente 1ͣ | Fresco | 18˚21’12.2” | 93˚11’54.6” |
| 34 | Guayaquil | Tabasco | Comacalco | Ra. Reyes Hernández | Seco | 18˚13’52.0” | 93˚15’09.5” |
| 35 | Trinitario | Tabasco | Comacalco | Ra. Reyes Hernández | Seco | 18˚13’52.0” | 93˚15’09.5” |
| 36 | Pataxte | Tabasco | Cárdenas | Nueva Zelandia | Seco | 17˚54’05” | 93˚35’11” |
| 37 | Criollo Blanco | Tabasco | Paraíso | Oriente 1ͣ | Seco | 18˚21’12.2” | 93˚11’54.6” |
| 38 | Criollo Rojo | Tabasco | Paraíso | Oriente 1ͣ | Seco | 18˚21’12.2” | 93˚11’54.6” |
| 39 | Pataxte | Tabasco | Pichucalco | Los Cascabeles | Seco | 17˚50’16” | 93˚11’33” |
| 40 | Pataxte | Tabasco | Paraíso | Oriente 1ͣ | Seco | 18˚21’12.2” | 93˚11’54.6” |
| 41 | Patastillo | Tabasco | Tenosique | El Palmar | Seco | 17˚28’00” | 91˚25’00” |
| 42 | Calabacillo | Tabasco | Cunduacán | Pechucalco | Seco | 18˚06’42.2” | 93˚11’13.6” |
| 43 | Guayaquil Rojo | Tabasco | Cunduacán | Pechucalco | Seco | 18˚06’42.2” | 93˚11’13.6” |
| 44 | Pataxte | Tabasco | Comalcalco | Ra. Reyes Hernández | Seco | 18˚13’52.0” | 93˚15’09.5” |
| 45 | Criollo/ Calabacillo | Tabasco | Comalcalco | Ra. Reyes Hernández | Seco | 18˚13’52.0” | 93˚15’09.5” |
| 46 | Guayaquil | Tabasco | Centro | Plátano y Cacao 1ͣ Sección | Seco | 17˚57’47.5” | 93˚10’23.0” |
| 47 | Guayaquil | Chiapas | Reforma | Macayo | Seco | 17˚58’20.6” | 93˚16’27.2” |
| 48 | Pataxte | Chiapas | Reforma | Macayo | Seco | 17˚58’20.6” | 93˚16’27.2” |
| 49 | Patastillo | Tabasco | Cunduacán | La Trinidad | Seco | 18˚09’12.0” | 93˚16’11.7” |
| 50 | Pataxte | Tabasco | Nacajuca | Bosques de Saloya | Seco | 18˚05’49” | 92˚56’41” |
| 51 | Pataxte | Tabasco | Comalcalco | R/A Norte 1ͣ Sección | Seco | 18˚21’12.2” | 93˚15’09.5” |
| 52 | Guayaquil | Tabasco | Paraiso | Oriente 1ͣ | Fresco | 18˚21’12.2” | 93˚11’54.4” |
| 53 | Criollo | Tabasco | Paraiso | Oriente 1ͣ | Fresco | 18˚21’12.2” | 93˚11’54.4” |
| 54 | Calabacillo | Tabasco | Cunduacán | Huapacal | Fresco | 18˚3’2.52” | 93˚9’21.25” |
| 55 | Pataxte | Tabasco | Teapa | Cerro de Puyacatengo | Seco | 18˚07’27.7” | 92˚55’86” |
| 56 | Guayaquil Rojo | Tabasco | Teapa | Cerro de Puyacatengo | Fresco | 18˚07’27.7” | 92˚55’86” |
Extracción de ADN
La caracterización molecular de las accesiones de cacao y pataxte se realizó en el Laboratorio de Biotecnología de Alimentos de la Unidad Académica Multidisciplinaria Reynosa-Aztlán de la Universidad Autónoma de Tamaulipas (UAMRA - UAT). La semilla de cada una de las 56 colectas se hizo germinar en charolas de 21 cavidades con una mezcla de “Perlita y Peatmoss” en proporción 40/60. Después de la germinación se colectaron las primeras dos hojas de cada plántula, de las cuales se extrajo el DNA con el kit comercial ‘DNeasy Plant Mini Kit’ de Quiagen( sin modificaciones del instructivo. Después de la extracción, el ADN genómico obtenido de cada una de las accesiones se cuantificó mediante un equipo NanoDrop( de Thermo Fisher, posterior a ello, se realizó una mezcla de colorante para tinción de ácidos nucléicos SYBR Gold en dilución 1:100,000 con 2 (L de DNA genómico, la mezcla se colocó en gel de agarosa al 1% en una cámara de electroforesis cuyas condiciones de electroforesis fueron 95 Volts y 50 mA por una hora, posterior a esto, las muestras se visualizaron mediante un fotodocumentador de luz UV ‘Gel DocTM EZ Imager’ de Bio-Rad( .
Condiciones de amplificación y visualización de los microsatélites
Los cuatro marcadores microsatélites empleados en el estudio se presentan en el Cuadro 2 y fueron seleccionados por ser informativos de acuerdo con el trabajo realizado por Saunders et al. (2004). La amplificación del ADN de las muestras colectadas se realizó mediante la técnica de Reacción en Cadena de la Polimerasa (PCR) en un termociclador ‘T100’ de Bio-Rad(, con temperaturas acordes a las condiciones de amplificación indicadas por Saunders et al. (2004). La PCR se realizó en un volumen final de 25 µL compuesto por: 70 ng de ADNg (molde), 50 pmol de cada iniciador, 200 µM de dNTP´s, búfer PCR 5X (Tris-HCl a 10 mM, pH 9; KCl a 50 mM; Triton X-100 al 0.1%), MgCl2 a 4 mM y 1 U de Taq DNA polimerasa de Promega (Madison, WI, EE.UU.). Los productos amplificados se visualizaron y cuantificaron en geles de agarosa al 1.2% ‘NuSieve’ de Lonza( con búfer SB (Borato de sodio decahidratado 5 mM) para alta resolución de los alelos producto de la amplificación y se analizaron mediante el software GelAnalizer 19.1 con marcadores de peso estándar para determinar los tamaños de los fragmentos amplificados por SSR.
Cuadro 2 Descripción de los marcadores microsatélites empleados en la caracterización de accesiones de cacao y pataxteTable 2. Description of the microsatellite markers used in the characterization of cacao and pataxte accessions
| SSR | Secuencias | Tm (C | Pb | Motivo |
|---|---|---|---|---|
| mTcCIR11 | F-5’TTTGGTGATTATTAGCAG3’ R-3’GATTCGATTTGATGTGAG5’ | 46 | 286 - 321 | (TC)13 |
| mTcCIR12 | F-5’TCTGACCCCAAACCTGTA3’ R-3’ATTCCAGTTAAAGCACAT5’ | 46 | 165 - 256 | (CATA)4N18(TG)6 |
| mTcCIR18 | F-5’GATAGCTAAGGGGATTGAGGA3’ R-3’GGTAATTCAATCATTTGAGGATA5’ | 51 | 331 - 354 | (GA)12 |
| mTcCIR58 | F-5’TTTTTGGTGATGGAACTAT3’ R-3’TGGTTAAGCAACACTAAACT5’ | 52 | 254 - 274 | (GT)40 |
Tm= temperatura de alineamiento; Pb= pares de bases; Fuente: Saunders et al. (2004).
Análisis estadístico
Durante la visualización de los geles, se anotó el peso molecular del producto amplificado para cada alelo de cada colecta evaluada donde se tomó como referencia el un marcador de peso molecular de 50 a 1,000 pares de bases Promega( G3161 y con ello, se estimó la frecuencia alélica para cada microsatélite probado en el presente estudio. La estimación de la diferenciación genética entre colectas, se realizó mediante un Análisis molecular de varianza (AMOVA), con el software GenAlEx 6.5.1 (Peakall & Smouse, 2012), y de acuerdo con De Vicente et al. (2004), los parámetros de diversidad estimados en esta investigación son:
HT = heterocigosidad esperada en la población total, estimada a partir de las frecuencias alélicas combinadas
Ho = heterocigosidad promedio observada en un grupo de poblaciones
He = heterocigosidad esperada, estimada a partir de cada subpoblación
FIS = 1 -
FST = 1 -
FIT =1 -
P=
Na= número de alelos por locus
Ne=número promedio de alelos por locus
Los datos de los pesos moleculares fueron transformados a una matriz de disimilitud alélica con la fórmula siguiente:
Donde:
d ij = disimilitud entre las unidades de i y j
L = número de loci
π = ploidía
m i = número de alelos pares para el locus l
Con la matriz de disimilitud generada se construyó un dendrograma mediante el método ‘UnWeighted Neighbor Joining’ “Vecinos no ponderados con el software DARWin 6.0.014 (Perrier et al., 2003). La estructura de las poblaciones analizadas se determinó con el software STRUCTURE V 2.3.4 (Chumacero et al., 2013) mediante ajuste del valor de ΔK que fue de uno a ocho con 200,000 ‘burn-ins periods’ “Periodos de alineamiento” y 800,000 repeticiones a través de veinte simulaciones (iteraciones) independientes para cada valor de ΔK, se empleó un modelo mezclado y correlacionado en cada caso. El número óptimo de poblaciones se estimó con la prueba de ΔK propuesta por Evanno et al. (2005) con el software disponible en línea STRUCTURE HARVESTER V 0.6.94 (Earl & vonHoldt, 2012).
Resultados
El rango alélico observado para cada marcador fue el siguiente: mTcCIR11 (121-313 pb), mTcCIR12 (190-232 pb), mTcCIR18 (333-357 pb) y mTcCIR58 (254-274 pb). Se observó exceso de homocigotos dado el valor de FIS=0.861 lo cual es derivado de la presencia de alelos nulos en las colectas evaluadas. Se encontraron 10 alelos para la variedad de cacao ‘Guayaquil’ y 10 para pataxte, el valor más alto para el número de alelos promedio (Ne=4), lo mostró la variedad ‘Guayaquil’, a diferencia del pataxte, que presentó tres alelos, pero superó al Ne=2 observado en la variedad Criollo. El valor más alto de heterocigosidad (He) lo presentó la variedad Guayaquil (0.595), pero fue inferior al 0.602 observado en Pataxte (Cuadro 3).
Cuadro 3 Parámetros obtenidos en la evaluación de la diversidad genética de poblaciones de cacao y pataxte colectadas en Tabasco y Chiapas, México.Table 3. Parameters obtained in the evaluation of the genetic diversity of cacao and pataxte populations collected in Tabasco and Chiapas, Mexico.
| Variedades | Na | Ne | He | Ho | Nivel de significancia |
|---|---|---|---|---|---|
| Variedades | Na | Ne | He | Ho | Nivel de significancia |
| Guayaquil | 10 | 4.00 | 0.595 | 0.100 | 0.000 |
| Criollo | 3 | 2.00 | 0.278 | 0.167 | 0.039 |
| Calabacillo | 6 | 2.50 | 0.458 | 0.083 | 0.000 |
| Patastillo | 8 | 3.00 | 0.543 | 0.188 | 0.000 |
| Guayaquil/Rojo | 2 | 1.75 | 0.281 | 0.125 | 0.111 |
| Pataxte | 10 | 3.00 | 0.602 | 0.165 | 0.000 |
| Promedio | 6.5 | 2.5 | 0.460 | 0.138 |
Na= Número de alelos por locus; Ne= Número promedio de alelos por locus; He= Heterocigosidad esperada; Ho= Heterocigosidad observada.
De los resultados arrojados por el Análisis molecular de varianza (AMOVA), se puede observar (Cuadro 4) que el 1.0% de la variación fue explicada por la moderada diferenciación genética entre las poblaciones de cacao y pataxte estudiadas. Sin embargo, se observó que el 85.0% de la varianza estimada se atribuye a diferencias entre los individuos, y 14.0% a la variación entre los individuos dentro de las poblaciones.
Con respecto a los índices de fijación o estadísticos de Wright, se observa que el FST presentó un valor de 0.011, lo que se puede interpretar como poca diferenciación genética debido a las frecuencias alélicas entre las accesiones evaluadas. El FIT de 0.861 y FIS también con un valor de 0.861, señalan exceso promedio de heterocigotos en promedio en un grupo de poblaciones y en cada una de las poblaciones evaluadas.
Cuadro 4 Análisis molecular de varianza de 56 accesiones de Theobroma spp. colectados en los estados de Tabasco y Chiapas, México.Table 4. Molecular analysis of variance of 56 Theobroma spp. accessions collected in the states of Tabasco and Chiapas, Mexico.
| Fuente de variación | Grados de libertad |
Suma de cuadrados |
Cuadrados Medios |
% de varianza explicada |
Índices de fijación |
|---|---|---|---|---|---|
| Entre poblaciones | 4 | 11.211 | 2.803 | 1.0 | FST= 0.011 |
| Entre individuos | 51 | 126.039 | 2.471 | 85.0 | FIT= 0.861 |
| Entre individuos dentro de poblaciones | 56 | 10.500 | 0.188 | 14.0 | FIS= 0.861 |
| Total | 111 | 147.750 | 100.0 |
El dendrograma (Figura 1), muestra las relaciones genéticas entre las accesiones evaluadas en el presente estudio, obsérvese que se formaron cuatro grupos o clústers donde se consideró una distancia de similitud igual 0.2. En el primer grupo (clúster 1) se agruparon 18 accesiones, entre ellas, se encuentran dos de Chiapas (Chis_Pataxte-1 y Chis_Patastillo-55), y el resto es germoplasma colectado en el estado de Tabasco. Este agrupamiento es posible que se haya debido a que presentaron similitud en el número de bandas identificadas por los marcadores en cada accesión, además de su origen geográfico. Dentro de este primer grupo (clúster 1) predominaron Pataxte (5) y Patastillo (4). El segundo grupo (clúster 2) se conformó por 13 accesiones y al igual que en el primer clúster, en éste, también predominaron los tipos Pataxte (4) y Patastillo (3). El grupo 3 (clúster 3) agrupó a 18 variedades de cacao y fue la variedad Guayaquil la que sobresalió en mayor número, cuatro de Chiapas, y tres de Tabasco, luego, el Pataxte y Patastillo con cuatro colectas cada uno, y por último el Calabacillo con tres colectas. El último grupo (clúster 4) fue en el que se agrupó el menor número de colectas (7) donde predominaron en número, los tipos Calabacillo con tres colectas y Criollo con dos.
El análisis de estructura poblacional (Figura 2) clasificó a las variedades de cacao y pataxte evaluadas en tres grupos (ΔK=3). El grupo 2 lo conformaron 23 accesiones, el grupo 3 se integró por 19 accesiones y el grupo 1 fue el de menor número de accesiones ya que solamente incluyó 14 accesiones. Obsérvese (Figura 2) que los pataxtes se encuentran distribuidos en cada uno de los tres grupos en cantidad más o menos proporcional por grupo.
Discusión
Los marcadores moleculares tipo microsatélite (SSR) en el presente estudio resultaron ser una herramienta útil y eficiente para determinar la diversidad de las colectas de cacao y pataxte evaluadas, lo cual quedó demostrado debido al número de alelos por locus amplificado. Cabe mencionar el trabajo de Chumacero et al., (2013) quienes emplearon los microsatélites mTcCir11, mTcCIR12 y mTcCIR18 entre otros, en reacciones de PCR multiplex, con lo cual pudieron detectar el flujo génico y la diversidad genética entre variedades de cacao silvestres y cultivadas. El microsatélite, que por los resultados arrojados, mostró más eficiencia en la presente investigación, además de ser significativo es el mTcCIR12, sin embargo, Smulders et al., (2012) en su investigación realizada en cacao en Trinidad y Tobago, reportaron que el microsatélite mTcCIR18 resultó ser el más eficiente en cuanto al número de alelos amplificados, a diferencia del presente trabajo, pero que coincide con los resultados del marcador mTcCIR12 que reportó en cacao Londoño-Murillo et al. (2011), y menciona que este marcador es uno de los más informativos y polimórficos, esta discrepancia puede deberse a que los autores referidos emplearon plata para teñir sus geles, y en la presente investigación se usó SYBR Gold lo que permitió una tinción de banda muy eficiente en todos los microsatélites evaluados.
Al comparar los resultados obtenidos en la presente investigación con los encontrados por (Aikpokpodion et al., 2010; Londoño-Murillo et al., 2011; Loor Solorzano et al., 2012; Vázquez-Ovando et al., 2014) quienes en sus trabajos emplearon los mismos microsatélites, se puede establecer que el número de alelos y rangos fueron consistentes entre sí, lo que se corrobora también por lo reportado por Saunders et al. (2004). En el presente estudio, se identificaron pequeñas variaciones en la estimación de los pesos moleculares de los alelos, pero esto puede ser un efecto de la sensibilidad que se presenta al realizar la lectura del programa GelAnalyzer 19.1 (Lazar et al., 2010). O también del ambiente en que se desarrollan las diferentes colecciones de cacao o pataxte de cada sitio experimental.
En relación al número de alelos por locus reportados en algunos trabajos con cacao se tienen los de Aikpokpodion et al. (2010) con 3.25, Aikpokpodion et al. (2009) con 4.39 alelos, Sereno et al. (2006) con 4.45 alelos, Boza et al. (2014) que reportó 3.14 alelos, y Velayutham et al., (2013) encontraron 3 alelos, que son valores menores al 6.50 de alelos promedio por locus detectados en este trabajo para las variedades cultivadas de cacao y pataxte, y similar al promedio de 6 alelos por locus reportado por Vázquez-Ovando et al. (2014).
Es posible que el número variable de alelos por locus reportados en los estudios previamente mencionados se deba a que se incluyeron materiales cultivados procedentes de distintas poblaciones, y por ello, cuando se comparan los resultados de este trabajo con estudios realizados en una única población se obtienen polimorfismos similares en cuanto al número de alelos amplificados.
Cabe mencionar que los resultados obtenidos se basan en la comparativa de las diferencias en cuanto a los tamaños y número de alelos que se pueden obtener, los cuales, en gran medida se deben a la técnica de detección empleada (Smulders et al., 2010). Mientras que la detección por radioactividad (NUCTECH) y por electroforesis capilar (poliacrilamida) permite la identificación de un mayor número de alelos, el uso de nitrato de plata es efectivo en menor proporción para la detección de alelos, pero con la ventaja de ser de bajo costo, mientras que en el presente trabajo, se optó por emplear la tinción de SYBR Gold para visualizar ácidos nucleicos ello, por su alta similitud con la tinción de bromuro de etidio (EtBr) una vez que se expone a luz ultravioleta.
En relación a la varianza encontrada en el AMOVA, la variación de 1% de varianza explicada entre poblaciones puede deberse a que las accesiones evaluadas presentan poca diferenciación, la cual puede deberse a la similitud génica o alélica que existe ente ellas. Lo anterior se refuerza por lo encontrado por Aikpokpodion et al. (2010), quienes al evaluar 11 poblaciones (híbridos, poblaciones, y colecciones) del banco de germoplasma de Nigeria, encontraron 37% de varianza entre poblaciones y 67% dentro de poblaciones.
Con respecto a los estimadores de la estructura poblacional, Opoku et al. (2007) obtuvieron un valor de FST = 0.076, no obstante que este índice de fijación indica moderada diferenciación genética entre las subpoblaciones evaluadas (no existe divergencia genética), los autores mencionan haber observado grandes diferencias genéticas entre las accesiones evaluadas, con lo cual separaron el germoplasma mejorado del criollo de los agricultores. En los resultados de sus trabajos de investigación Aikpokpodion et al. (2009); Loor Solorzano et al. (2012); Ruiz et al. (2011); y Vázquez-Ovando et al. (2014), encontraron valores de FST que fluctuaron de 0.123, 0.209, 0.4979 y 0.230, los que son superiores al FST= 0.011 estimado en la presente investigación, es decir, las poblaciones evaluadas de los investigadores citados mostraron diferenciación genética de moderada a muy grande, mientras que las accesiones evaluadas en el presente estudio, presentan poca diferenciación o pequeña, misma que es posible que sea producto del intercambio de germoplasma entre agricultores de cacao de las localidades estudiadas en cada estado. Con respecto al valor FIS=0.861 se encontró más exceso de homocigotos que el reportado (FIS=0.2402) por Loor Solorzano et al. (2012) oeste índice de fijación, lo cual pone de manifiesto la poca diferenciación genética que existe entre las accesiones evaluadas.
En cuanto a la heterocigosidad esperada (He=0.460), y la observada (Ho=0.138) el valor promedio encontrado para He, superó al reportado por Sereno et al. (2006); Chumacero et al. (2013); Vázquez-Ovando et al. (2014) y Gopaulchan et al. (2020), pero inferior al estimado por Aikpokpodion et al. (2010); y menor para Ho del trabajo de Ruiz et al. (2011) y Batista et al. (2016) . Resultados que pueden deberse, en parte a la constitución genética del germoplasma evaluado en cada trabajo.
Con relación al agrupamiento encontrado en las accesiones evaluadas, obsérvese (Figura 1 y 2), que el pataxte se mezcló con los demás tipos o accesiones de cacao, esto se contrapone con lo asegurado por Azurdia, (2014) y Avalos et al., (2012) de que no existe relación filogenética entre el cacao criollo y el pataxte, y que por ello la hibridación natural entre ambas especies no es posible que se lleve a cabo.
Conclusiones
La investigación realizada en la que se emplearon cuatro microsatélites, generó información con la que se puede establecer la hipótesis de que es posible que haya alguna relación filogenética entre el cacao y pataxte que se cultivan en los estados de Tabasco y Chiapas. Lo anterior se sustenta en el hecho de que es común que los agricultores de los estados referidos cultivan juntos en sus plantaciones al cacao y pataxte, y seleccionen semilla de las mejores mazorcas de cacao para tener nuevas plantas en sus predios, las que es posible que pudieran ser producto de posible hibridación, ya sea entre variedades de cacao o entre cacao y pataxte.
El mayor número efectivo de alelos (Ne) se encontró en las accesiones de cacao Guayaquil, que puede ser resultado del posible entrecruzamiento entre las variedades mencionadas.










nueva página del texto (beta)





