Introducción
Los granos de polen son estructuras microscópicas que se encuentran en las anteras de los estambres de las plantas angiospermas. Esos granos se adhieren a los cuerpos de las abejas obreras que visitan las flores por medio de fuerzas electrostáticas débiles (Clarke et al. 2017), varios granos de polen son aglutinados con las patas traseras y son humedecidos con secreciones salivales y néctar para formar gránulos de 1 a 4 mm, los cuales son transportados en las curbículas de la abeja y almacenados en alveolos del panal para su posterior consumo (Thakur y Nanda 2020).
El polen apícola se puede recolectar colocando una trampa en la entrada de las colmenas, la cual permite desprender el polen de las patas traseras de las abejas que regresan a la colmena. El polen apícola debe consumirse fresco o secarse cuidadosamente para que conserve sus propiedades nutricionales (Kieliszek et al. 2018). La composición del polen apícola varía dependiendo de factores como la estacionalidad, la ubicación geográfica, la humedad, entre otros. Estas variaciones le confieren diferentes propiedades nutricionales, fisicoquímicas y funcionales, no solo al polen sino también a los otros productos apícolas tales como la jalea real, la miel y el propóleo (Galimberti et al. 2014). El polen de abeja es rico en carbohidratos, proteínas y lípidos, contiene todos los aminoácidos y minerales esenciales, ácidos grasos y vitaminas principalmente del complejo B, carotenoides y flavonoides (Thakur y Nanda 2018, 2020). Actualmente se consume principalmente como suplemento dietético, por lo cual ha suscitado gran interés la investigación sobre este alimento, desde poder caracterizarlo fisicoquímicamente hasta la determinación de su origen botánico (Thakur y Nanda 2018, Végh et al. 2021).
El origen botánico del polen afecta el color de los gránulos de polen el cual puede variar desde blanco hasta café oscuro (Thakur y Nanda 2020). El polen se considera monofloral si presenta un único origen botánico en una frecuencia mayor al 90%, o bien no contar con polen accesorio en una frecuencia mayor al 60%, cuando no hay un polen predominante se considera polen multifloral (Thakur y Nanda 2020). Conocer cuáles son las plantas que son visitadas por las abejas es de gran importancia para los apicultores, para asegurar la calidad y seguridad de sus productos. Esta información proporciona datos valiosos sobre cuáles son las fuentes relevantes de néctar y la preferencia alimenticia de las abejas en diferentes épocas del año (Tanaka et al. 2020). Lo anterior permite elaborar estrategias para preservar la salud y productividad de las colmenas, determinar las características florales de una región y sus cambios fenológicos, detectar la presencia de especies invasivas, derivado de lo anterior, funciona como herramienta de monitoreo de la vegetación (Milla et al. 2022). Tradicionalmente el origen botánico del polen ha sido identificado por medio de análisis microscópico del tamaño, forma y superficie de los granos de polen, para asociarlo con una especie de plantas. Sin embargo, este método demanda un amplio conocimiento de las características distintivas del polen de cada especie, y frecuentemente este método resulta en una identificación imprecisa (Galimberti et al. 2014, Végh et al. 2021).
Las técnicas basadas en biología molecular muestran un gran potencial para superar esas limitaciones, ya que el análisis de secuencias estandarizadas como los códigos de barras de ADN permiten identificar y discriminar exitosamente plantas, incluso de taxones cercanos. En plantas, los marcadores propuestos por el Consortium for the Barcode of Life (CBOL) son las regiones del gen de la Maturasa K (matK), y de la ribulosa bifosfato carboxilasa (rbcL), adicionalmente se agregó el espaciador interno transcrito (ITS2) debido a su alto poder de discriminación de especies (Chen et al. 2010). Aunque estos fragmentos estandarizados permiten la discriminación entre especies, en muchos casos prácticos, no son capaces de identificar a todas las especies de plantas o incluso diferenciar entre plantas de un mismo grupo taxonómico (Shneyer y Rodionov, 2019). Un punto clave para obtener exitosamente códigos de barras de ADN, es la extracción del material genético, donde se requiere un método rápido, fácil y económico, del cual resulte un alto rendimiento y calidad de ADN. Para garantizar la reproducibilidad de los resultados es necesario descartar la presencia de inhibidores de la PCR (Safeena et al. 2021). Un ADN libre de inhibidores permite realizar no solo la amplificación de los códigos de barras, si no también hace posible la detección de polen de plantas transgénicas, alergénicas o bien de plantas invasoras (Guertler et al. 2014). En el presente estudio, se evaluó el efecto del tiempo de incubación, numero de lavados y el solvente utilizado sobre la cantidad y calidad del ADN genómico extraído de polen extraído por medio de la técnica con bromuro de hexadecil trimetil amonio (CTAB). Además, se determinó la utilidad del ADN extraído para identificar molecularmente algunas de las plantas que son visitadas por las abejas para colectar polen.
Materiales y métodos
Purificación del ADN genómico de polen apícola
Se utilizó polen apícola multifloral colectado en la comunidad de Las Rosas, Chiapas (16°22’ 00” LN 92° 21’ 00” LO). Para la extracción del ADN, se llevaron a cabo 12 variantes del protocolo previamente documentado para extraer ADN de plantas (Sharma et al. 2003). Las variaciones fueron en el tiempo de incubación inicial (tres variaciones), lavados (dos variaciones) y el solvente utilizado (dos variantes) (Tabla 1). El protocolo consistió en que 1 g de polen fue sometido a diferentes tiempos de incubación a 50 °C (0, 10 y 30 minutos) en 12 mL de tampón de fosfatos salino (PBS). Posterior a la incubación se centrifugó a 4 000 rpm por 5 min, el sobrenadante se decantó lentamente cuidando de no desprender la pastilla. Se realizaron lavados agregando 8 mL de agua ultrapura a la pastilla, se agitó durante 30 s con un agitador vortex y se centrifugó a 4 000 rpm por 20 min y se descartó el sobrenadante, este proceso se realizó de una o tres veces según el tratamiento. Posterior al lavado, se agregaron 500 µL de amortiguador de extracción CTAB comercial (Promega Corporation, Madison, WI, USA) y se mezcló durante 40 s con un agitador vortex. La solución se transfirió a un microtubo de 1.5 mL estéril y los tubos fueron incubados a 65°C a 200 rpm, durante 30 min. Luego se agregaron 0.6 volúmenes del solvente cloroformo-alcohol isoamílico (C-I) en proporción 24:1 o fenolcloroformo-alcohol isoamílico (F-C-I) en proporción 25:24:1 según el tratamiento. Se agitó durante 10 s en un agitador vortex y se centrifugó cada tubo a 15 000 rpm durante 10 min. Con la ayuda de una micropipeta se tomó la fase acuosa y se depositó en un microtubo de 1.5 mL estéril y se agregaron dos volúmenes de etanol absoluto para posteriormente centrifugar a 15 000 rpm por 5 min. Se decantó el sobrenadante y para lavar la pastilla formada se agregaron 400 µL de etanol al 70%, se volvió a centrifugar a 15 000 rpm por 5 min y se descartó el exceso de etanol y dejó secar a temperatura ambiente. La pastilla fue rehidratada con 100 µL de agua ultrapura libre de nucleasas y se incubó a 50 °C a 200 rpm por 15 min. La integridad del ADN genómico extraído fue verificada por medio de una electroforesis en gel de agarosa (Tris-acetato-EDTA 1x, 0.8% agarosa, 90V).
Tabla 1. Variantes de los tratamientos analizados.
Tiempo de Incubación (min) | Numero de Lavados | Solvente | Tratamiento |
---|---|---|---|
0 | 1 | C-I | T1 |
F-C-I | T2 | ||
2 | C-I | T3 | |
F-C-I | T4 | ||
10 | 1 | C-I | T5 |
F-C-I | T6 | ||
2 | C-I | T7 | |
F-C-I | T8 | ||
30 | 1 | C-I | T9 |
F-C-I | T10 | ||
2 | C-I | T11 | |
F-C-I | T12 |
C-I: cloroformo-alcohol isoamílico en proporción 24:1, F-C-I: fenol-cloroformo-alcohol isoamílico en proporción 25:24:1
Cuantificación del ADN genómico
Con el equipo Nanodrop (NANODROP 2000/2000c Spectrophotometer, Thermo Scientific) se realizó por triplicado la cuantificación del rendimiento de ADN genómico y las proteínas obtenidas, utilizando la aplicación para ácidos nucleicos NANODROP 2000/2000c v.1.0 la cual arroja como resultados la cantidad de ADN en ng y las absorbancias a 280, 260 y 230 nm. Para esto se utilizó 1 µL agua ultrapura como solución de referencia (blanco) y de las muestras de ADN, las cuales fueron depositadas individualmente en la superficie de medición. La intensidad de luz de la muestra y del blanco son requeridas para calcular la absorbancia de una longitud de onda determinada (Thermo Scientific; v1.0 User Manual, NanoDrop 2000/2000c Spectrophotometer, 2009).
Índices A260/A280 y A260/A230
Los índices A260/A280 y A260/A230 de cada muestra fueron determinados automáticamente por la aplicación NANODROP 2000/2000c v.1.0. Se calculan dividiendo el valor de absorbancia a 260 nm entre el valor a 280 nm y, dividiendo el valor de absorbancia a 260 nm entre el valor a 230 nm, respectivamente.
Amplificación del marcador ITS2
El ADN obtenido en cada tratamiento se utilizó para la amplificación por PCR del fragmento ITS2 usando cebadores específicos (ITS-2F:5’-ATGCGATACTTGGTGTGAAT-3’, ITS-2R:5’-GACGCTTCTCCAGACTACAAT-3’) previamente documentados (Chen et al. 2010). Para la PCR se utilizó el kit comercial GoTaq Green Master Mix, 2XTM (Promega Corporation, Madison, WI, USA), con reacciones de 25 µL de mezcla según las especificaciones del fabricante, las cuales se colocaron bajo un programa de amplificación (95°C, 4min; 35x (94°C, 45 s; 56°C, 1 min; 72°C, 1 min); 72°C, ∞ min) en un termociclador (Axygen® MaxyGene™ Thermal Cycler II), de cada tratamiento se establecieron reacciones por triplicado. Los productos de PCR fueron visualizados por medio de electroforesis en gel de agarosa (Tris-acetato-EDTA 1x, 0.8% agarosa, 90V), la banda putativa correspondiente al fragmento ITS2, amplificada utilizando el DNA genómico del tratamiento T1 (Tabla 1), fue purificada utilizando el kit Wizard® SV Gel y PCR Clean-Up System (Promega Corporation, Madison, WI, USA), de acuerdo con las instrucciones del fabricante.
Clonación y secuenciación de los fragmentos ITS2
Los amplicones purificados fueron clonados usando el kit comercial pGEM-T Easy Vector System (Promega Corporation, Madison, WI, USA), siguiendo las especificaciones del fabricante. El resultado de la ligación se utilizó para transformar cepas de E. coli DH5α calcio competentes. Las cepas transformantes resultantes fueron cultivadas en placas que contenían medio agar LB con 50 mg mL-1 de ampicilina. Las colonias que resultaron resistentes a la ampicilina fueron cultivadas en medio LB liquido con 50 mg mL-1 de ampicilina con la finalidad de realizar la extracción de los vectores recombinantes utilizando el kit comercial Wizard Plus SV Minipreps DNA Purification Systems (Promega Corporation, Madison, WI, USA), según las especificaciones del fabricante. El ADN plásmidico fue verificado por medio de electroforesis en gel de agarosa y utilizado como sustrato para realizar la PCR utilizando cebadores universales M13. Un total de 15 vectores presentaron un amplicon del tamaño esperado los cuales fueron secuenciados de forma unidireccional con el cebador M13F en Langebio-CINVESTAV (Irapuato, México) según las especificaciones del servicio.
Análisis estadístico
Se realizó un análisis de varianza de tres vías (factores: tiempo de incubación * número de lavados * tipo de solvente, p ≤ 0.05). Para determinar si hubo diferencias entre los rendimientos de ADN, proteínas, índices A260/A280 y A260/A230 y, además determinar si hay interacción entre los factores. Se realizaron graficas comparativas de las medias (+/- error estándar) de cada variable evaluada (rendimiento de ADN, rendimiento de proteínas, índice 260/280 e índice 260/230) con respecto a cada factor (tiempo de incubación * número de lavados * tipo de solvente) y sus variantes. Utilizando el software R-studio 4.1.2. En las secuencias resultantes se delimitó el inserto por medio de la comparación de la secuencia obtenida con la secuencia del pGEM-T Easy Vector, las secuencias obtenidas fueron analizadas con la herramienta BLASTn (NCBI 2022) considerando un valor de E <1 × 10-5 y aciertos máximos de 98 a 100% para considerarse perteneciente a una especie en particular (Ross et al. 2008).
Resultados
Rendimiento de ADN
La electroforesis en gel de agarosa del ADN genómico extraído permitió verificar la integridad del material genético en todos los tratamientos realizados (Figura 1). El análisis de varianza (p ≤ 0.05) indica que el tiempo de incubación influye de forma significativa sobre el rendimiento del ADN obtenido, con diferencias altamente significativas derivadas del número de lavados realizados y del tipo de solvente utilizado (Tabla 2). Por otra parte, los tratamientos se tradujeron en diferencias en rendimiento de proteínas (p ≤ 0.05), jugando un papel relevante el tiempo de incubación, al igual que el número de lavados y el solvente en la proteína obtenida. El rendimiento del ADN y de proteínas, no tuvieron evidencias de interacción entre los tratamientos, por lo que se puede inferir que el tiempo de incubación, lavados y tipo de solvente no tienen influencia entre sí (Tabla 2).
Tabla 2 Resultados del ANOVA de tres vías de cada una de las variables analizadas.
Factor | ADN p-value | Proteínas p-value | Índice 260/280 p-value | Índice 260/230 p-value |
---|---|---|---|---|
Incubación | 0.07165* | 0.046* | 0.078 | 0.007** |
Lavados | 2.19x10-04** | 1.69x10-4** | 0.066 | 0.006** |
Solvente | 0.0022** | 7.80x10-4** | 0.020* | 0.017* |
Incubación*lavados | 0.0521 | 0.093 | 0.0018** | 0.002** |
Incubación*solvente | 0.22943 | 0.169 | 0.524 | 0.588 |
Lavados*solvente | 0.29819 | 0.394 | 0.49 | 0.573 |
Incubación*lavados* Solvente | 0.7463 | 0.770 | 0.129 | 0.661 |
*,** significativo y altamente significativo, respectivamente.
Los resultados de la cuantificación, señalan que el número de lavados reduce el rendimiento de ADN extraído y, a su vez el uso de C:I como solvente resulta en mayores rendimientos del ADN, aunque no hay diferencias estadísticas. Los rendimientos obtenidos en los diferentes tiempos de incubación fueron mayores en ausencia de incubación y más bajos con 10 y 30 min de incubación. Se evidencia que el mejor protocolo para obtener mayor rendimiento de ADN, no incluye tiempo de incubación a 50°C, un solo lavado del polen y el uso de C:I como solvente (Figura 2A). El uso de C:I tuvo mayor rendimiento de las proteínas comparado con el rendimiento obtenido usando F:C:I (Figura 2B).
![](/img/revistas/era/v10nspe3//2007-901X-era-10-iii-e3656-gf2.jpg)
Figura 2 A) Rendimiento de ADN y B) Rendimiento de proteínas promedios obtenidos (+/- error estándar).
Los tratamientos con mayores rendimientos del ADN permitieron obtener valores de hasta 180 ng µL-1, y los tratamientos con menores rendimientos registraron rendimientos del ADN de hasta 120 ng µL-1, es decir, hasta un 33% menos. El rendimiento de proteínas exhibió un comportamiento similar al del rendimiento del ADN; es decir, los valores más altos del rendimiento de proteínas se presentaron en aquellos tratamientos que tuvieron mayores rendimientos del ADN (Figura 2 A y B), siendo ambos afectados por los mismos factores, indicando una fuerte asociación entre ambos tipos de moléculas. Los mayores valores de rendimiento de proteínas alcanzar hasta 2 mg mL-1 y los menores que se registraron fueron de hasta 1.32 mg mL-1, es decir, hasta un 35% menos.
Índices 260/280 y 260/230
Se encontraron diferencias significativas en el índice 260/280 relacionadas a los diferentes solventes utilizados. También se detectó interacción entre los tiempos de incubación y el número de lavados (Tabla 2). En cuanto al índice 260/230, se encontró diferencia significativa relacionada con el tiempo de incubación, número de lavados del polen, y solvente utilizado; además de interacción entre los tiempos de incubación y número de lavados (Tabla 2). Los índices 260/280 más bajos se obtuvieron en aquellos tratamientos que no contaron con tiempo de incubación, alcanzando un valor medio de 1.8 cuando se incluyó un tiempo de incubación de 30 min, el cual es considerado óptimo (Figura 3A). Los índices 260/230 obtenidos fueron bajos, oscilando entre 1 y 1.3, aunque ninguna muestra obtuvo un índice 260/230 mayor a 1.8 para ser considerado aceptable (Figura 3B).
Amplificación del fragmento ITS2 y secuenciación
La amplificación de los fragmentos ITS2 por PCR ocurrió en todas las muestras de ADN obtenidas. De las 15 muestras enviadas a secuenciación, se obtuvieron secuencias cuyo tamaño estuvo entre 365 y 482 pares de base, de las cuales fue posible determinar la especie en 11 secuencias. En el análisis por medio de Blastn indicó que cuatro muestras presentaron un porcentaje de identificación menor al 98%, lo cual permite determinar únicamente el género, en la Tabla 3 se indica la especie más probable. Tres muestras secuenciadas presentaron dos especies probables y, en cuatro muestras se identificó una misma especie (Melampodium divericatum).
Tabla 3 Especies identificadas por medio de Blastn.
Género | Especie | Nucleótidos idénticos | Porcentaje de identidad |
---|---|---|---|
Helianthus | annus/maximilani | 479/495 | 96.77 |
Bidens | alba/pilosa | 487/492 | 98.98 |
Juglans | Regia | 437/453 | 96.47 |
Brassica | Juncea | 441/445 | 99.1 |
Alnus | Cordata | 482/498 | 96.79 |
Melampodium | Divaricatum | 365/365 | 100 |
Sigesbeckia | Orientalis | 423/468 | 90.4 |
Melampodium | Divaricatum | 368/368 | 100 |
Melampodium | Divaricatum | 383/383 | 100 |
Melampodium | Divaricatum | 368/369 | 99.7 |
Ageratum | conyzoides/houstonianum | 378/378 | 100 |
Discusión
Rendimiento de ADN
La identificación molecular del origen botánico del polen apícola, requiere un método de extracción de ADN que permita obtener buenos rendimientos y la pureza necesaria para que no interfiera con los análisis posteriores. Uno de los principales retos a vencer en la extracción del ADN de polen apícola es lograr el rompimiento de la exina, la cual es una estructura dura que cubre los granos de polen (Thakur y Nanda 2020, Swenson y Gemeinholzer 2021).
Los protocolos que se han desarrollado, en su mayoría, incluyen mecanismos físicos para la ruptura de la exina, incluyendo el uso de proteinasas o de nitrógeno líquido junto al uso de kits comerciales (Galimberti et al. 2014, Tanaka et al. 2020). Además de los que incluyen el uso de DNAzol (Thermo Fischer Scientific) y amortiguador de extracción CTAB (125 mM Tris-HCL, 50mM EDTA, 50 mM NaCL) (Torricelli et al. 2016, Manivanan et al. 2018). La mayoría de los reportes previos enfocados al estudio del polen apícola, analizan el polen que está contenido en la miel, lo cual hace necesario adicionar una metodología para extraer los granos de polen de la miel. Los protocolos que se han documentado para este objetivo incluyen ciclos de disolución de la miel en agua y una posterior centrifugación para sedimentar los granos de polen (Guertler et al. 2014, Bell et al. 2016, Manivanan et al. 2018). En este estudio, la extracción se realizó directamente con polen apícola, utilizando un protocolo basado en el uso del amortiguador de extracción CTAB.
El rendimiento del ADN obtenido en este trabajo oscila entre 120-180 ng µL-1, a la fecha existen pocos estudios que permitan establecer una comparativa real. La mayoría los estudios palinológicos con técnicas moleculares, parten de una muestra de miel a la cual se le extrae el polen y posteriormente realizan la extracción de ADN genómico, obteniendo rendimientos inferiores a 50 ng µL-1 (Galimberti et al. 2014, Bruni et al. 2015, Torricelli et al. 2016). En este estudio se utilizó polen apícola, la cual es una mezcla compleja que contiene moléculas de distinta naturaleza. Un trabajo similar fue realizado por Galimberti et al. (2014) quienes extrajeron ADN genómico de polen apícola, utilizando un kit comercial (Plant DNeasy Isolationy Purification kit, Qiagen, Italy) y nitrógeno líquido permitió obtener rendimientos de hasta 20 ng µL-1.
Índices 260/280 y 260/230
Los índices 260/280 y 260/230 son útiles en la detección de como proteínas, sales y polisacáridos en muestras de ADN; los cuales pueden inhibir e interferir con los análisis posteriores como la amplificación o secuenciación. La pureza del ADN puede ser afectada por varios componentes presentes en el polen apícola y las plantas en general, tales como polisacáridos, lípidos, polifenoles o bien compuestos químicos utilizados en la extracción. El índice 260/280 de 1.8 es indicativo que una muestra de ADN tiene alta pureza y es libre de proteínas y fenoles (Lucena-Aguilar et al. 2016, Abdel-Latif y Osman 2017). Durante la extracción del ADN, se liberan diversas sustancias endógenas contaminantes, sin embargo, el problema más común es la degradación del ADN genómico por la presencia de formas oxidadas de polifenoles, que se unen covalentemente al ADN dañándolo. El β-mercaptoetanol, que forma parte del amortiguador de extracción CTAB, es un fuerte agente reductor que previene la oxidación de los polifenoles en las muestras vegetales (Youssef et al. 2015).
En este estudio, el índice 260/280 se ubicó en el rango de 1.76 a 1.82, pudiendo resaltar la pureza del ADN extraído sin el uso de nitrógeno líquido. Se considera que un índice 260/280 de 1.93 a 2.27 indica niveles insignificantes de contaminación por proteínas, polisacáridos y ARN, y un nivel menor de 1.7 indica la presencia de fenol y otros contaminantes (Pervaiz et al. 2011, Abdel-Latif y Osman 2017, Safeena et al. 2021). Las variaciones en el índice 260/280 son atribuibles al método de extracción. Este índice resultó más cercano a 1.8 cuando se realiza una incubación de 30 min a 50°C, lo cual disminuye los rendimientos de proteínas, mismas que aumentaron en caso de no realizar la incubación.
En cuanto al índice 260/230, es considerado por diversos autores como una medida secundaria de la pureza del ADN (Lucena-Aguilar et al. 2016), las muestras aquí analizadas presentaron valores entre 1.0 y 1.3. Estudios previos consideran que los valores aceptables del índice 260/230 oscilan entre 2.0 y 2.2, mientras que valores menores a 1.8 indican la contaminación con fenoles, compuestos orgánicos, EDTA, carbohidratos y polisacáridos (Lucena-Aguilar et al. 2016, Safeena et al. 2021). Existen cuestionamientos sobre la confiabilidad de este índice, debido a que el amortiguador de elución incrementa la concentración de sales en la muestra, reduciendo el índice 260/230 (Lucena-Aguilar et al. 2016). En este trabajo, la rehidratación del ADN se realizó con agua ultrapura, por lo que los bajos índices 260/230 son atribuibles a otros pasos de la extracción. A pesar del índice 260/230 bajo, la calidad del ADN genómico extraído fue suficiente para amplificar el fragmento ITS2. Estos resultados apoyan los cuestionamientos sobre la confiablidad de este índice, a la fecha existen diferentes reportes que indican que un índice 260/230 bajo no interfiere en análisis moleculares basados en amplificación (Shepherd et al. 2015).
Amplificación del fragmento ITS2 y secuenciación
Los códigos de barras del ADN de origen plastídico (como rbcL y matK), presentan cierta problemática, debido al número variable de plástidos en diferentes especies y a los casos de herencia biparental, lo cual interfiere en la identificación (Bell et al. 2016). Sin embargo, han sido previamente utilizados de modo generalizado en estudios de identificación de polen apícola basado en métodos moleculares (Galimberti et al. 2014, Guertler et al. 2014, Manivanan et al. 2018, Tanaka et al. 2020, Milla et al. 2022).
A la fecha se reportan cinco marcadores como códigos de barras de ADN para la identificación de plantas, tres de ellos (rbcL, matK e ITS2) tienen amplia aceptación en el Barcode of Life Database (BOLDSYSTEMS 2022). No obstante, el marcador y cebadores seleccionados pueden representar fuentes de sesgo durante la amplificación por PCR. Todos los marcadores presentan diferente tasa de amplificación en múltiples especies, frecuentemente es necesario rediseñar cebadores para un grupo taxonómico especifico, lo cual conlleva a un sesgo en los amplicones producidos, donde el ADN blanco tiene fuentes muy diversas (Kraaijeveld et al. 2015, Bell et al. 2016). En estudios que involucran una especie o un número limitado de individuos, la combinación de múltiples locus puede resultar en mayor poder de discriminación. Sin embargo, en el caso del polen un análisis multilocus debe ir acompañado de un análisis bioinformático para determinar con mayor exactitud las especies presentes en la muestra.
En este trabajo se utilizó el marcador molecular ITS2 para la identificación de las especies de plantas. Este marcador no plastídico, puede presentar problemas debido a la ploidía de la especie en cuestión, pero al mismo tiempo ha demostrado tener mayor resolución taxonómica (Chen et al. 2010) y ser un discriminador eficiente entre especies en innumerables familias de plantas (Timpano et al. 2020). A pesar de su alta resolución especifica en plantas, este marcador tiene la desventaja de potencialmente poder amplificar ADN fúngico (Bell et al. 2016, Timpano et al. 2020).
Los cebadores utilizados permitieron la amplificación del fragmento ITS2 en las muestras de ADN de cada uno de los tratamientos, lo cual resalta la buena calidad del ADN obtenido. Con relación a la universalidad de los cebadores ITS2, debido a que el polen apícola es una mezcla con diversos taxones representados, es de esperarse que algunos de ellos permitan la amplificación, no obstante, esto no implica mayor abundancia de una especie en particular, aunque si resalta que cuatro muestras resultaran en una misma especie (M. divericatum). Entre las especies identificadas destaca Alnus cordata, considerada endémica de Italia (Smeriglio et al. 2019, Villani et al. 2021), junto con Juglans regia fueron las únicas especies arbóreas identificadas. Mientras que el porcentaje de identidad de A. cordata es marcadamente bajo, por lo cual hay duda sobre su identidad a nivel especifico. La mayoría de las especies identificadas son consideradas herbáceas no perennes, entre las que se identificó a Ageratum conyzoides especie considerada tóxica (Kotta et al. 2020).
La secuenciación de las muestras se realizó con la técnica de Sanger, la cual es útil para hacer inferencias en una escala detallada, pero es impráctica en la aplicación de estudios metagenómicos. La tecnología de secuenciación de próxima generación (NGS) permite la secuenciación de alto rendimiento para realizar análisis metagenómicos, permitiendo el estudio a nivel bioma, incluyendo a los organismos de un hábitat especifico. Otros autores han reportado estudios metagenómicos donde la identificación a nivel especie es limitada o nula (Richardson et al. 2015, Tanaka et al. 2020). A la fecha diversos estudios con ADN genómico de polen demuestran que una herramienta eficiente para el monitoreo ambiental, ya que permite detectar no solo de la flora relevante en la alimentación de las abejas, sino también la identificación de plantas tóxicas (Bruni et al. 2015), alergénicas (Kraaijeveld et al. 2015), invasoras (Noël et al. 2023) y transgénicas (Wu et al. 2017, Bayram 2021).
Conclusiones
Las modificaciones realizadas al protocolo de extracción del ADN basado en CTAB, tuvieron efectos en el rendimiento de ADN genómico y la contaminación por proteínas. En términos generales, realizar un lavado y el uso de cloroformo-alcohol isoamílico resultaron en mayor rendimiento del ADN genómico con calidad suficiente para realizar la amplificación del marcador ITS2 y la identificación de especies de plantas. La información es importante para hacer eficiente el uso de recursos y tiempo para extraer ADN genómico de polen apícola. Este método permite la identificación de especies de relevancia ecológica para las abejas y demás polinizadores. Además, su utilidad puede extrapolarse al monitoreo de hábitats para la detección de especies de plantas invasoras, transgénicas, alergénicas, entre otras.