Variabilidad genética del membrillo cimarrón (Malacomeles denticulata [Kunth] Jones) obtenida mediante marcadores Inter Secuencias Simples Repetidas o Intermicrosatélites (ISSR)
Genetic variability of Mexican serviceberry (Malacomeles denticulata [Kunth] Jones) obtained by inter simple sequence repeated or intermicrosatellites (ISSR) markers
RESUMEN
Malacomeles denticulata es un fruto nativo de México al que recientemente se le ha encontrado características funcionales para proponerlo como una opción frutal. Esta investigación tuvo como objetivo elucidar la variabilidad de doce poblaciones de M. denticulata mediante marcadores Inter Secuencias Simples Repetidas o Intermicrosatélites (ISSR). Todos los ISSR presentaron altos valores en el Contenido de Información Polimórfica (PIC, por sus siglas en inglés) y el índice de diferenciación poblacional de Nei (GST), así como altos porcentajes de polimorfismo. Se crearon tres grupos de variabilidad, donde siete poblaciones de Guanajuato y la población del Tepozán, Querétaro, conformaron un grupo; mientras que las poblaciones de Agua Zarca (Guanajuato), La Joya (Querétaro) y Santa Catarina del Monte (México) conformaron otro. Finalmente, la población de San Miguel Tlaixpan (México) quedó separada como un grupo atípico. De acuerdo con el Análisis Molecular de Varianza (Amova), la variabilidad intrapoblacional representa el 47% y la interpoblacional el 53% de la variabilidad total, lo que concuerda con la variabilidad en caracteres de semilla, pero difiere a lo reportado para marcadores SSR.
ABSTRACT
Malacomeles denticulata is a native fruit of Mexico that recently had been reported functional proprieties to be proposed as a fruit option. This research aims to elucidate the variability of twelve populations of M. denticulata by Inter Simple Sequence Repeated or Inter-microsatellites (ISSR) markers. All ISSR showed high values of Polymorphic Information Content (PIC) and Nei’s Index of Population Differentiation (GST) as well as high percentage of polymorphism. Three groups of variability were conformed, where seven populations of Guanajuato and the population of El Tepozán, Querétaro conformed the first group; while the populations of Agua Zarca (Guanajuato), La Joya (Querétaro), and Santa Catarina del Monte (México) conformed the other group. Finally, the population of San Miguel Tlaixpan (México) was separated as outlier group. According to Analysis of Molecular Variance (Amova), the within population variability was 47% and among population variability was 53% of total variability, that agree on the variability of seed traits but disagree on the variability of SSR markers.
INTRODUCCIÓN
El membrillo cimarrón (Malacomeles denticulata [Kunth] Jones) es una planta perteneciente a la familia de las Rosaceae subtribu Pyrinae (anteriormente la subfamilia Maloideae) (Campbell, Evans, Morgan, Dickinson & Arsenault, 2007; Turner, 2011). Malacomeles denticulata se encuentra distribuida desde el sur de Texas (Estados Unidos de América) hasta Guatemala y Honduras, teniendo en México su mayor distribución; por las características climáticas, las regiones de crecimiento óptimas están en los estados de Hidalgo, Estado de México, Querétaro, Guanajuato, Tlaxcala, Puebla, Nuevo León y Oaxaca (Núñez-Colín, 2010).
El género Malacomeles tiene una estrecha relación con los géneros Amelanchier y Peraphyllum. En este sentido, especies del género Amelanchier son explotadas de manera comercial en Estados Unidos de América y Canadá, tanto para arquitectura del paisaje como para producción de fruta (Stushnoff, 1991).
Para hacer un primer acercamiento a la variabilidad de M. denticulata, Hernández-Martínez et al. (2011) utilizaron características de semillas de membrillo cimarrón como patrón de variabilidad. En este estudio determinaron una alta variación tanto intra como interpoblacional, lo cual sugiere que esta especie puede ser de reproducción sexual.
Torres-Hernández et al. (2013) reportaron que a partir de 23 SSR de manzana (Malus domestica) aplicados a nueve poblaciones de membrillo cimarrón determinaron una alta heterocigocidad y falta de equilibrio Hardy-Weinberg en todas las poblaciones. Ello sugiere que esto podría estar relacionado con autoincompatibilidad gametofítica (AIG).
Debido a estos antecedentes, en este trabajo se pretendió determinar y validar la variabilidad de M. denticulata con marcadores moleculares Inter Secuencias Simples Repetidas o Intermicrosatélites (ISSR) que, a diferencia de los SSR, son marcadores semialeatorios (González & Aguirre, 2007) que brindan información complementaria para ampliar el conocimiento sobre la variabilidad genética del membrillo cimarrón.
MATERIALES Y MÉTODOS
Material vegetal
Durante el otoño de 2011 y primavera y otoño de 2012, una serie de viajes de colecta fueron hechos en las zonas con probable presencia de membrillo cimarrón basado en Núñez-Colín (2010). En el transcurso de esos viajes se obtuvieron once poblaciones silvestres de membrillo cimarrón en las comunidades de La Valenciana (dos poblaciones), El Derramadero, Victoria, Puerto de Ardillas, Agua Zarca y La Huerta en el estado de Guanajuato, así como en El Tepozán y La Joya en el estado colindante de Querétaro, y también en Santa Catarina del Monte y San Miguel Tlaixpan en el Estado de México. Adicionalmente, se consideró una población establecida en el Arboretum de Plantas Nativas de Guanajuato del Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP), la cual originalmente fue colectada de La Calera, municipio de San Miguel de Allende, Guanajuato (figura 1). La muestra utilizada de cada población fue de diez individuos colectados de manera aleatoria en cada caso. Las muestras colectadas se colocaron en sílica gel hasta su traslado al laboratorio, donde fueron almacenadas a –20 °C, para después liofilizarlas y molerlas en tubos tipo falcón; finalmente se almacenaron a –20 °C.
Figura 1. Ubicación de las poblaciones de Malacomeles denticulata estudiadas.
Fuente: Elaboración propia.
Figura 1. Ubicación de las poblaciones de Malacomeles denticulata estudiadas.
Fuente: Elaboración propia. Close
Aislamiento de ADN
El aislamiento de Ácido Desoxirribonucleico (ADN) se realizó con base en el protocolo Bromuro de hexadeciltrimetilamonio (CTAB) modificado con Cloruro de Sodio-Tris-EDTA (STE) (Doyle & Doyle, 1987); se utilizaron hojas jóvenes, sin manchas y sin evidencias de plagas. Su calidad y pureza fueron determinados utilizando el espectrofotómetro Thermo scientific Nanodrop 8000 version 2.1.0.
Amplificación por PCR de los ISSR
Se seleccionaron 17 iniciadores ISSR utilizados comúnmente en organismos vegetales (tabla 1). Para la amplificación por Reacción en cadena de la Polimerasa (PCR, por sus siglas en inglés) se utilizó la siguiente mezcla de reacción: 10 µL de buffer 1 x, 3 mM MgCl2, 0.2 mM dNTP, 1 pMol de primer y 1 U Taq Polimerasa recombinante (Invitrogen®), donde se usó un volumen de 20 µL con aproximadamente 20 ng de ADN. El perfil de temperatura utilizado fue el siguiente: 94 °C por 5 min, 35 x (94 °C por 60 s, Tm °C [tabla 1] por 60 s, 72 °C por 60 s), 72 °C por 10 min.
Tabla 1.
Lista de iniciadores ISSR utilizados en doce poblaciones de Malacomeles denticulata.
Nota: los iniciadores degenerados son mezclas de iniciadores PCR, donde los símbolos en la respectiva fórmula representan: Y = C o T.
Fuente: Elaboración propia.
Abrir
Nota: los iniciadores degenerados son mezclas de iniciadores PCR, donde los símbolos en la respectiva fórmula representan: Y = C o T.
Fuente: Elaboración propia. Close
Los productos de PCR fueron desnaturalizados con una solución de carga de formamida (98% formamida, 10 mM EDTA pH 8.0, ~40 mg azul de bromofenol) a 95 ºC por 5 min, después fueron puestos en hielo. El gel de poliacrilamida al 8.5% fue hecho con 40 mL de solución de poliacrilamida (19:1, Urea 8 M), 200 µL de persulfato de amonio 25% y 40 µL de TEMED. 12 µL de cada producto de PCR fue cargado en el gel. La electroforesis fue hecha usando el buffer TBE 1 x en una cámara de electroforesis vertical a 250 volts por 90 min. Finalmente, los geles de poliacrilamida fueron teñidos con bromuro de etidio para visualizar los fragmentos mediante luz UV.
Análisis estadísticos
El índice de información de Shannon (I), los estadísticos de diferenciación de subpoblaciones de Nei (HS, HT, GST), la prueba de flujo genético (NM), el porcentaje de polimorfismo (%P) y la distancia genética de Nei (Nei, 1972) fueron calculados con el programa POPGENE versión 1.32 (Yeh, Yang & Boyle, 1999). Además, a partir de los cálculos de POPGENE se dibujó un dendrograma con la distancia genética de Nei y el método de agrupación en pares por media aritmética no ponderada (UPGMA, por sus siglas en inglés) con TreeView versión 1.6.6
donde PIC es el contenido de información polimórfica del marcador $i$, $f_{i}$ la frecuencia de los fragmentos del marcador que estaban presentes y 1 – $f_{i}$ la frecuencia de los fragmentos del marcador que estaban ausentes. El PIC fue un promedio de los fragmentos para cada ISSR.
El PIC toma valores entre 0 y 0.5 para datos dominantes, como los ISSR, y se considera que valores menores a 0.15 es no informativo, valores entre 0.15 y 0.25 como informativo y mayor a 0.25 como altamente informativo (Laurentin & Karlovsky, 2007). También se calculó el análisis molecular de varianza (Amova) y un análisis de coordenadas principales (ACooP) utilizando el programa GenAlEx versión 6.5 (Peakall & Smouse, 2012).
RESULTADOS
Variabilidad intrapoblacional
Se encontró que las poblaciones presentaron diversos valores de HT, I y %P en las diferentes poblaciones (tabla 2).
El ISSR Pv02 fue monomórfico dentro de todas las poblaciones, Pv10 en las poblaciones de La Valenciana 2, La Joya y San Miguel Tlaixpan, el Pv09 en La Joya y Pv05 e ISSR19 en San Miguel Tlaixpan (tabla 2).
Se encontró que la población más homogénea fue la de San Miguel Tlaixpan, y la más heterogénea fue la de La Valenciana 1 (tabla 2). Las poblaciones variaron del 25.12% hasta el 55.02% de polimorfismo intrapoblacional.
Tabla 2.
Datos poblaciones del índice de heterogeneidad de Nei (HT), índice de Shannon (I) y porcentaje de polimorfismo (%P) en doce poblaciones de Malacomeles denticulata.
ISSR
M
El Derramadero
INIFAP
La Valenciana 2
Victoria
La Valenciana 1
Puerto de Ardillas
HT
I
% P
HT
I
% P
HT
I
% P
HT
I
% P
HT
I
% P
HT
I
% P
Pv01
10
0.1881
0.3019
0.7692
0.1385
0.2249
0.5769
0.0345
0.0493
0.0769
0.1646
0.2556
0.5769
0.1324
0.1963
0.3462
0.1556
0.2319
0.4615
Pv02
10
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
Pv03
10
0.2137
0.3218
0.6250
0.1921
0.2931
0.6250
0.1787
0.2704
0.5625
0.2689
0.3866
0.6250
0.1965
0.2949
0.6250
0.2057
0.3041
0.5625
Pv04
10
0.0530
0.0840
0.2308
0.1795
0.2611
0.4615
0.3211
0.4666
0.7692
0.1632
0.2517
0.5385
0.1481
0.2298
0.5385
0.1328
0.1916
0.3077
Pv05
10
0.2162
0.3250
0.6364
0.2162
0.3250
0.6364
0.2662
0.3945
0.7727
0.2894
0.4289
0.7727
0.3130
0.4561
0.8182
0.2544
0.3782
0.7273
Pv06
10
0.2113
0.3203
0.6500
0.2113
0.3203
0.6500
0.1669
0.2541
0.5000
0.1104
0.1626
0.3000
0.1918
0.2896
0.6000
0.1814
0.2616
0.4500
Pv07
10
0.0935
0.1413
0.2667
0.0935
0.1413
0.2667
0.1044
0.1541
0.2667
0.2043
0.3030
0.5333
0.2097
0.3388
0.8667
0.1331
0.1914
0.3333
Pv08
10
0.1651
0.2494
0.5000
0.1651
0.2494
0.5000
0.2050
0.3030
0.5500
0.2158
0.3161
0.6000
0.2731
0.4058
0.7500
0.2729
0.3946
0.7000
Pv09
10
0.2272
0.3486
0.8095
0.2272
0.3486
0.8095
0.1982
0.3075
0.6667
0.1710
0.2483
0.4286
0.1934
0.2788
0.4762
0.1051
0.1507
0.2381
Pv10
10
0.1987
0.2811
0.4286
0.1987
0.2811
0.4286
0.0000
0.0000
0.0000
0.0249
0.0381
0.0714
0.1677
0.2365
0.3571
0.0527
0.0788
0.1429
ISSR7
10
0.1887
0.2877
0.5882
0.1887
0.2877
0.5882
0.1164
0.1618
0.2353
0.1626
0.2304
0.3529
0.1788
0.2591
0.4706
0.0460
0.0681
0.1176
ISSR16
10
0.1445
0.2229
0.4615
0.3674
0.2229
0.4615
0.1318
0.1973
0.3846
0.0601
0.0944
0.2308
0.0715
0.1012
0.1538
0.0433
0.0662
0.1538
ISSR19
10
0.2314
0.3600
0.8000
0.1540
0.2320
0.4667
0.1936
0.2867
0.5333
0.2021
0.3048
0.6000
0.2480
0.3713
0.7333
0.1760
0.2584
0.4667
IS03
10
0.2762
0.4125
0.8000
0.2617
0.4058
0.9000
0.2820
0.4129
0.7500
0.2193
0.3486
0.8500
0.2422
0.3719
0.8000
0.2205
0.3239
0.5500
UBC830
10
0.3002
0.4333
0.7000
0.2343
0.3545
0.7000
0.2606
0.3740
0.6000
0.2939
0.4355
0.8000
0.2540
0.3610
0.6000
0.1760
0.2742
0.6000
UBC835
10
0.0801
0.1178
0.2174
0.1636
0.2473
0.4783
0.2203
0.3161
0.5217
0.1970
0.2935
0.5652
0.1488
0.2275
0.4783
0.0582
0.0830
0.1304
UBC841
10
0.2755
0.3993
0.6957
0.1352
0.2099
0.4783
0.1393
0.2115
0.4348
0.3673
0.5264
0.8696
0.2278
0.3504
0.7391
0.1624
0.2293
0.3478
Prom
0.1802
0.2710
0.5399
0.1839
0.2591
0.5310
0.1658
0.2447
0.4485
0.1832
0.2720
0.5126
0.1880
0.2805
0.5502
0.1398
0.2050
0.3700
ISSR
M
El Tepozán
La Huerta
La Joya
Santa Catarina del Monte
Agua Zarca
San Miguel Tlaixpan
HT
I
% P
HT
I
% P
HT
I
% P
HT
I
% P
HT
I
% P
HT
I
% P
Pv01
10
0.1471
0.2290
0.5385
0.1460
0.2227
0.4615
0.1224
0.1977
0.5000
0.1877
0.2753
0.5000
0.1857
0.2694
0.4615
0.1537
0.2316
0.4615
Pv02
10
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
Pv03
10
0.2179
0.3153
0.5625
0.1243
0.1882
0.3750
0.2729
0.3924
0.6250
0.2360
0.3387
0.5625
0.3509
0.5078
0.8750
0.2297
0.3258
0.5000
Pv04
10
0.1671
0.2395
0.3846
0.1191
0.1748
0.3077
0.0971
0.1400
0.2308
0.0150
0.0311
0.1538
0.0420
0.0685
0.1538
0.0394
0.0621
0.1538
Pv05
10
0.2110
0.3233
0.7273
0.3044
0.4550
0.8636
0.1389
0.2087
0.4091
0.1347
0.2090
0.5000
0.1304
0.1944
0.3636
0.0000
0.0000
0.0000
Pv06
10
0.2476
0.3644
0.6500
0.1695
0.2650
0.6000
0.1199
0.1732
0.3000
0.2141
0.3077
0.5000
0.1847
0.2660
0.4500
0.2880
0.4249
0.7500
Pv07
10
0.1963
0.2852
0.4667
0.1022
0.1570
0.3333
0.1834
0.2603
0.4000
0.0353
0.0551
0.1333
0.1581
0.2323
0.4000
0.1013
0.1495
0.2667
Pv08
10
0.0846
0.1226
0.2000
0.1417
0.2097
0.4000
0.1765
0.2595
0.4500
0.1095
0.1623
0.3000
0.1038
0.1506
0.2500
0.0665
0.0953
0.1500
Pv09
10
0.0988
0.1406
0.2381
0.0647
0.0922
0.1429
0.0000
0.0000
0.0000
0.1772
0.2571
0.4286
0.0647
0.0922
0.1429
0.0090
0.0161
0.0476
Pv10
10
0.1104
0.1640
0.2857
0.0716
0.1092
0.2143
0.0000
0.0000
0.0000
0.0877
0.1353
0.2857
0.0698
0.1116
0.2857
0.0000
0.0000
0.0000
ISSR7
10
0.0772
0.1110
0.1765
0.0273
0.0387
0.0588
0.0496
0.0718
0.1176
0.0827
0.1224
0.2353
0.1956
0.3035
0.6471
0.1835
0.2701
0.4706
ISSR16
10
0.1194
0.1768
0.3077
0.0934
0.1371
0.2308
0.0455
0.0685
0.1538
0.0934
0.1371
0.2308
0.0794
0.1199
0.2308
0.1261
0.1916
0.3846
ISSR19
10
0.1810
0.2775
0.6000
0.0619
0.0875
0.1333
0.0911
0.1366
0.2667
0.2843
0.4242
0.8000
0.2597
0.3787
0.6667
0.0000
0.0000
0.0000
IS03
10
0.3369
0.4946
0.9000
0.1472
0.2280
0.5000
0.1994
0.3056
0.7000
0.2501
0.3665
0.6500
0.1375
0.2016
0.3500
0.1286
0.1893
0.3500
UBC830
10
0.1772
0.2662
0.5000
0.2504
0.3752
0.7000
0.2239
0.3366
0.7000
0.3318
0.4808
0.8000
0.1979
0.2751
0.4000
0.1455
0.2034
0.3000
UBC835
10
0.0828
0.1250
0.2609
0.1144
0.1685
0.3043
0.1364
0.1991
0.3478
0.1284
0.1895
0.3478
0.1547
0.2212
0.3478
0.0460
0.0673
0.1304
UBC841
10
0.2288
0.3418
0.6522
0.2119
0.3045
0.5217
0.3134
0.4468
0.6957
0.2492
0.3601
0.6087
0.3717
0.5378
0.8696
0.1209
0.1780
0.3043
Prom
0.1579
0.2339
0.4383
0.1265
0.1890
0.3616
0.1277
0.1880
0.3469
0.1539
0.2266
0.4139
0.1580
0.2312
0.4056
0.0964
0.1415
0.2512
Prom = Promedio, M = Muestra.
Fuente: Elaboración propia. Abrir
Prom = Promedio, M = Muestra.
Fuente: Elaboración propia. Close
Variabilidad interpoblacional e índices de información de los marcadores ISSR
Todos los ISSR fueron polimórficos en las doce poblaciones, donde los porcentajes de polimorfismo fueron mayores del 85% y valores de PIC altamente informativos (> 0.25), con excepción del Pv02 que fue del 22.22 % y PIC de 0.0478 (tabla 3).
Los valores del índice de diferenciación entre poblaciones (GST) fueron altos (tabla 3) y destaca Pv02, ya que el polimorfismo que presenta en dos de las nueve bandas de este ISSR son específicos de poblaciones, lo cual permite identificarlos perfectamente; y el número de migrantes promedio dado por el índice de flujo genético (NM) varía entre 0 y 2.76 con un promedio de 1.0840, el cual es bajo, por tal motivo se traduce en que poca de la variación de las poblaciones puede explicarse por efectos de migración.
Tabla 3.
Datos de número de bandas, bandas polimórficas, porcentaje de polimorfismo (%P), el contenido de información polimórfica (PIC), índice de Shannon (I), índice de heterogeneidad de Nei (HT), índice de endogamia, índice de diferenciación entre poblaciones (GST) e índice de flujo genético (NM) de 17 ISSR aplicados a doce poblaciones de Malacomeles denticulata.
ISSR
M
Número de bandas
Bandas polimórficas
% P
PIC
I
HT
HS
GST
NM
Pv01
120
26
26
1.0000
0.2991
0.3905
0.2413
0.1464
0.3767
1.5540
Pv02
120
9
2
0.2222
0.0478
0.0819
0.0478
0.0000
1.0000
0.0000
Pv03
120
16
16
1.0000
0.3611
0.5483
0.3693
0.2239
0.4104
1.0785
Pv04
120
13
13
1.0000
0.3439
0.5177
0.3432
0.1231
0.6133
0.6134
Pv05
120
22
22
1.0000
0.4081
0.5704
0.3904
0.2038
0.4581
1.2893
Pv06
120
20
19
0.9500
0.3379
0.4893
0.3273
0.1864
0.4268
0.9126
Pv07
120
15
15
1.0000
0.2516
0.3890
0.2551
0.1412
0.3578
2.7605
Pv08
120
20
20
1.0000
0.3583
0.5261
0.3566
0.1535
0.5304
1.0521
Pv09
120
21
21
1.0000
0.3075
0.4720
0.3169
0.1248
0.5093
1.8498
Pv10
120
14
12
0.8571
0.2787
0.4544
0.3067
0.0679
0.7334
0.2636
ISSR7
120
17
16
0.9412
0.2807
0.4493
0.2939
0.1187
0.5487
0.7886
ISSR16
120
13
12
0.9231
0.3112
0.4799
0.3198
0.0926
0.6783
0.3162
ISSR19
120
15
15
1.0000
0.3503
0.5278
0.3566
0.1736
0.4371
1.2928
IS03
120
20
20
1.0000
0.4119
0.5772
0.3949
0.2251
0.4275
0.7877
UBC830
120
10
10
1.0000
0.3784
0.5350
0.3725
0.2371
0.3222
2.3320
UBC835
120
23
21
0.9130
0.3264
0.4601
0.3078
0.1276
0.5291
0.6971
UBC841
120
23
23
1.0000
0.3762
0.5814
0.3965
0.2336
0.4239
0.8397
Prom
17.4706
16.6471
0.9298
0.3194
0.4735
0.3175
0.1517
0.5167
1.0840
Prom = Promedio, M = Muestra.
Fuente: Elaboración propia.
Abrir
Prom = Promedio, M = Muestra.
Fuente: Elaboración propia. Close
Distancias genéticas y representación de agrupamiento de las poblaciones
La distancia genética de Nei (Nei, 1972) arrojó que las diferencias entre las distintas poblaciones de membrillo cimarrón variaron entre 0.0932 (9.32%) y 0.4194 (41.94%), lo que dio identidades genéticas entre el 0.6574 (65.74%) y 0.9110 (91.10%) (tabla 4).
La población de San Miguel Tlaixpan fue la que presentó la mayor distancia genética con todas las otras poblaciones, que entre ellas presentaron diferencias menores al 30% (tabla 4).
Tabla 4.
Distancia e identidad genética de Nei (1972) entre doce poblaciones de Malacomeles denticulata.
En el dendrograma se observa que la población de San Miguel Tlaixpan es diferente al resto, mientras que las poblaciones de Guanajuato, a excepción de Agua Zarca, están agrupadas en un mismo grupo, incluyendo a la población de El Tepozán (Querétaro), y quedaron en otro grupo las poblaciones de La Joya (Querétaro), Agua Zarca (Guanajuato) y Santa Catarina del Monte (Estado de México) (figura 2).
Figura 2. Dendrograma utilizando la distancia genética de Nei (1972) y método de aglomeración UPGMA de doce poblaciones de Malacomeles denticulata utilizando marcadores ISSR.
Fuente: Elaboración propia.
Figura 2. Dendrograma utilizando la distancia genética de Nei (1972) y método de aglomeración UPGMA de doce poblaciones de Malacomeles denticulata utilizando marcadores ISSR.
Fuente: Elaboración propia. Close
En el ACooP, el cual acumuló el 54.01% de la variabilidad total al tercer componente, se corroboran las agrupaciones visualizadas en el dendrograma (figura 3), donde se observa que la población de San Miguel Tlaixpan está separada del resto, principalmente por la Coordenada Principal 1 (CooP1), mientras que el grupo de La Joya, Agua Zarca y Santa Catarina del Monte está separado del grupo de poblaciones de Guanajuato y El Tepozán, principalmente por la Coordenada Principal 2 (CooP2) (figura 3).
En el grupo de las poblaciones de Guanajuato y en Tepozán se presentaron ligeras disimilitudes entre El Tepozán, Puerto de Ardillas y La Huerta con el resto del grupo en la Coordenada Principal 3 (CooP3). No obstante, por su cercanía en el CooP1 y CooP2 se corrobora que se trata de un mismo grupo.
Figura 3. Representación tridimensional de doce poblaciones de Malacomeles denticulata utilizando marcadores ISSR obtenido del análisis en coordenadas principales.
Fuente: Elaboración propia.
Figura 3. Representación tridimensional de doce poblaciones de Malacomeles denticulata utilizando marcadores ISSR obtenido del análisis en coordenadas principales.
Fuente: Elaboración propia. Close
Análisis Molecular de la Varianza (Amova)
En el Amova, la variabilidad intrapoblacional representa el 47% de la variación general, y el 53% corresponde a la variabilidad interpoblacional (tabla 5). Por lo anterior, aunque es ligeramente mayor la variabilidad interpoblacional, la variabilidad de esta especie se debe tanto a su variación intra como interpoblacional.
Tabla 5.
Análisis molecular de la varianza de doce poblaciones de Malacomeles denticulata usando 17 ISSR.
Los resultados obtenidos en este estudio presentaron disimilitudes con el de Torres-Hernández et al. (2013), debido principalmente a que las poblaciones de La Valenciana y el Tepozán fueron más similares en esta investigación a las poblaciones de Guanajuato, mientras que en el trabajo de Torres-Hernández et al. (2013) estaban como diferentes grupos de variación; caso contrario fue la población de La Joya, que en este caso queda separada del resto y que estaba agrupada con poblaciones de Guanajuato. Esto puede deberse a que los SSR, marcadores más específicos, definen las agrupaciones con diferencias genéticas puntuales (Nybom & Weising, 2010), mientras que los ISSR, marcadores que exploran todo el genoma, pueden ayudar a encontrar polimorfismos aleatorios en diferentes regiones del mismo para diferenciar las poblaciones de manera general o agrupar poblaciones parecidas genéticamente, pero con diferencias que sólo pueden ser elucidadas con marcadores específicos, como los SSR (Weising, Nybom, Wolff & Kahl, 2005).
Los polimorfismos en los ISSR tienden a ser altamente variables entre individuos debido a las altas tasas de mutación que experimentan, ya que cuando el ADN se replica durante la meiosis, la ADN polimerasa puede “tartamudear” hacia adelante o hacia atrás en las unidades repetidas, eliminando o agregando unidades a la cadena microsatétlite (González & Aguirre, 2007). Las cadenas resultantes pueden, entonces, presentar menos o más unidades de repetición (o pares de bases) que las cadenas parentales (Zietkiewicz, Rafalski & Labuda, 1994). Además de que su codificación y análisis es como marcadores dominantes, donde cada banda se considera un solo locus con dos alelos, por lo que no se sabe si la presencia de la banda se debe a un homocigoto dominante o un heterocigoto. SSR e ISSR, por tanto, son marcadores complementarios que permiten elucidar de una mejor manera la variabilidad real de una especie.
A diferencia de Torres-Hernández et al. (2013), con los ISSR la variabilidad interpoblacional es mayor y es casi igual la variabilidad interpoblacional a la intrapoblacional, lo que concuerda más con el estudio de Hernández-Martínez et al. (2011), con características morfológicas de semilla, quienes reportaron una alta variabilidad tanto intrapoblacional como interpoblacional, asunto comprobado en este trabajo.
No obstante, en este trabajo se consideraron poblaciones del Estado de México no consideradas por trabajos anteriores, así como una población atípica morfológicamente de Guanajuato, como la población de Agua Zarca, principalmente por tener las hojas más redondas que las otras poblaciones de Guanajuato, las cuales en este estudio se encuentran en grupos de variabilidad diferentes, en particular la población de San Miguel Tlaixpan, identificada como un grupo atípico (outlier). Por ello, estas poblaciones junto con la población de La Joya podrían ser una especie diferente a las de Guanajuato, como lo reportó Turner (2011), quien reclasificó a M. denticulata var. psilantha como M. psilantha y que la distribución reportada por este autor para dicha especie corresponde con el de las poblaciones del estado de Guanajuato, con excepción de Agua Zarca, y el Tepozán, Querétaro; mientras que el resto de las poblaciones podría tratarse de M. denticulata sensu stricto. Sin embargo, esto da cabida a que se haga un trabajo taxonómico sobre el género Malacomeles para asegurar si se trata de dos especies o solamente una, sobre todo cuando una de ellas era considerada una variedad botánica de la otra. Por lo tanto, este estudio podría servir de base como primer acercamiento para poder distinguirlas, por lo menos como grupos de variabilidad distintos. No obstante, un estudio como el realizado por Lo, Stefanovi$\grave{\mathrm{c}}$ & Dickinson (2007), quienes secuenciaron genes conservados nucleares y de cloroplasto en el género Crataegus, sería más apropiado para concluir si se trata de dos especies diferentes.
CONCLUSIÓN
Los marcadores ISSR aplicados a M. denticulata presentaron altos valores de PIC y de polimorfismo en general, motivo por el cual se consideran adecuados para estudios de variabilidad para esta especie.
Se encontró que la agrupación de las poblaciones obtenida con estos marcadores fue diferente a la reportada por Torres-Hernández et al. (2013), así como la distribución de la variabilidad inter e intrapoblacional; sin embargo, por la diferencia en el tipo de marcador se pueden considerar complementarios para conocer la diversidad real de la especie.
AGRADECIMIENTOS
Este trabajo fue financiado parcialmente por el Fondo de Ciencia Básica del Consejo Nacional de Ciencia y Tecnología (Conacyt) por el proyecto CB2009/134193, así como por los fondos fiscales del Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP) por el proyecto SIGI 1145232069.