nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
El maíz es el cereal más producido en todo el mundo según los últimos datos de la Organización para la Agricultura y la Alimentación, con una superficie cosechada de 205 ×106 ha y una producción de 1210 × 106 Mg en el año 2021, superando al arroz y al trigo en cuanto a productividad (FAOSTAT, 2021), debido a que entre los cereales posee el más alto valor nutricional (Singh, Singh y Singh, 2020), es considerado un componente importante para la alimentación humana y animal. En México el cultivo del maíz es el más representativo por su importancia cultural, económica y social (SAGARPA, 2017). Dentro de la clasificación del maíz se encuentra el criollo que durante años ha sido mejorado por los campesinos, al realizar un intercambio de semillas mediante la selección empírica y con ello, intercambio de genes (Sánchez-Hernández, De la Cruz y Sánchez, 2015). Como resultado de esta selección existen diferentes razas de maíz, clasificadas según sus características isoenzimáticas y morfológicas (Sánchez, Goodman y Stuber, 2000). Navarro-Garza, Hernández, Castillo y Pérez (2012) destacan que el término raza, en las plantas, se refiere al conjunto de características o fenotipo y que se le asigna un nombre el cual evidencia alguna asociación en especial. En México existe un reporte de 64 razas de maíz, de las cuales 59 se consideran nativas y cinco que fueron descritas inicialmente en otras regiones (CONABIO, 2020).
En cuanto a rendimiento, el maíz criollo se ve superado considerablemente por las variedades mejoradas, no obstante, los pequeños agricultores se inclinan por las variedades nativas dado que poseen algunas ventajas, entre las cuales se encuentran: a) la capacidad para sembrarse en terrenos edáfica y climáticamente más limitados (Turrent, Wise y Garvey, 2012), reportándose razas que sobreviven donde los maíces mejorados no lo hacen (Vázquez-Carrillo, Pérez, Hernández, Marrufo y Martínez, 2010); b) aptitud para elaborar platillos tradicionales; c) estabilidad y adaptabilidad frente a climas locales y fluctuantes; d) reducción de costos en insumos para su producción; y e) mejor manejo del riesgo agrícola (Guillén-Pérez, Sánchez, Mercado y Navarro, 2002; Turiján-Altamirano, Damián, Ramírez, Juárez y Estrella, 2012; Turrent et al., 2012). La Secretaría de Agricultura y Desarrollo Rural (SADER), informó que en 2020 la producción de maíz grano fue de 27 707 775 toneladas, mientras que en 2021 se reportó un crecimiento de 2.6% en la producción con una estimación de 28 427 436 toneladas (SADER, 2021).
El factor edáfico comprende las propiedades físicas, químicas y biológicas, las cuales influyen en el desarrollo de las plantas (Rajakaruna y Boyd, 2008). Dentro de las características biológicas del suelo, las comunidades microbianas son elementales en el desempeño de las plantas, al controlar su fisiología y crecimiento, siendo muy grande la diversidad microbiana, asociada con las raíces de las plantas (Karnwal, 2022). La rizósfera se define como el suelo más próximo a las raíces de una planta, y el término “rizobacteria” engloba aquellas bacterias que están asociadas a dichas raíces en la rizósfera, las cuales pueden tener un efecto benéfico, neutral o dañino para la misma planta (Odelade y Babalola, 2019). Entre las rizobacterias se encuentran las bacterias promotoras del crecimiento vegetal (BPCV), grupos de diversas especies que mejoran el crecimiento de las plantas (González y Fuentes, 2017). Estas bacterias simbióticas que habitan la rizósfera de ciertas especies del maíz, además de favorecer su productividad, también participan en la defensa de la salud del cultivo sin dañar al ambiente (Akhtar, Qureshi, Iqbal, Ahmad y Khan, 2012; Raza, Ling, Yang, Huang y Shen, 2016).
Cada vez hay más estudios que demuestran que el genoma de la planta tiene un impacto dinámico y directo en la estructura de su microbioma (Aira, Gómez, Lazcano, Bååth y Domínguez, 2010; Cregger et al., 2018), este último definido como “…todo el hábitat, incluidos los microorganismos (bacterias, arqueas, eucariotas inferiores, superiores y virus), sus genomas (es decir, genes) y las condiciones ambientales circundantes…” (Marchesi y Ravel, 2015). Considerando la asociación existente entre el genoma y el microbioma de la planta, es indispensable sopesar de manera conjunta al hospedero y el microbioma para fortalecer el conocimiento sobre el funcionamiento de la planta como un holobionte (Ravanbakhsh, Kowalchuk y Jousset, 2020). Por lo tanto, el estudio del microbioma bacteriano rizosférico es fundamental para comprender el rendimiento de la planta. En el caso del maíz, el rendimiento tiene una relevancia considerable debido a que, a nivel regional y nacional, se espera un aumento continuo de la población, de ahí que la demanda de este cereal se incremente para ser destinado al consumo humano y animal. El presente estudio fue realizado sobre una transición agroecológica, cuestión imprescindible en la actualidad, ya que las técnicas agroecológicas aseguran un uso sostenible de los recursos naturales, el cuidado del ambiente y de las personas, permitiendo que se pueda garantizar la alimentación de cada familia a largo plazo. En atención a lo descrito anteriormente, el objetivo del presente estudio fue, describir el microbioma rizosférico de bacterias y su impacto sobre el rendimiento en tres razas de maíz criollo de grano, utilizadas en la producción de forraje, grano y doble propósito.
Materiales y Métodos
Ubicación y condiciones ambientales del área de estudio
La Comarca Lagunera presenta una precipitación pluvial promedio anual de 258 mm, con temperatura media anual de 21 °C, con máxima y mínima de 33.7 y 7.5 ºC, respectivamente (Montemayor-Trejo et al., 2012). Se recolectaron muestras del microbioma rizosférico de suelo (franco-arenoso, pH 8.1, C.E. 1.8 mS cm-1, materia orgánica de 1.1%) del Campo Experimental El Retiro, de la Universidad Autónoma Agraria Antonio Narro, ubicado en el municipio de Francisco I. Madero, Coahuila de Zaragoza, México (25º 49’ 53.45” N y 103º 7’ 4.03” O), durante el desarrollo de tres razas de maíz criollo, para las cuales se ocuparon semillas del banco de germoplasma del Departamento de Agroecología de la propia institución: C1, raza forrajera; C2, raza doble propósito (grano y forraje) y C3, raza para grano, cada una con tres repeticiones etiquetadas como C1-1, C1-2, C1-3, C2-1, C2-2, C2-3, C3-1, C3-2 y C3-3.
Establecimiento del cultivo y manejo de las muestras
La siembra del cultivo, de las tres razas, de maíz criollo se efectuó el 13 de junio del 2021 en melgas de 120 m de largo y 10 m de ancho. Los surcos presentaron una separación de 0.75 m y una distancia entre plantas de 0.20 m. Las razas se sembraron en dos surcos de 60 m con cuatro repeticiones, y dos surcos de bordo entre éstas. La preparación del suelo consistió en doble rastreo, bordeo, con la aplicación de un riego de pre siembra de 0.20 m y se aplicaron cuatro riegos de auxilio de 0.15 m de lámina cada uno, utilizando un sistema de multi-compuertas. Se aplicó biol como repelente y fertilizante foliar cada ocho días, a los 40 días después de sembrado. El manejo de maleza se realizó de manera manual en una ocasión con dos jornales. A los 51 días se realizó el muestreo de la rizósfera. Se extrajo la planta de maíz y se obtuvieron 100 mg de suelo asociado a las raíces. Cada muestra fue colocada en un tubo bashing bead Zymo ResearchTM con 750 µL de buffer lisante/estabilizador XpeditionTM; estos tubos se procesaron en un disruptor celular TerraLyzerTM, durante 45 segundos, para la lisis celular y la conservación del DNA. Posteriormente, las muestras se transportaron al laboratorio en hielera portátil.
Extracción de DNA del microbioma del suelo
El DNA de las muestras se extrajo utilizando el kit DNA MiniPrep de Zymo ResearchTM y se cuantificó con un nanodrop ThermoScientific. Para efectuar la amplificación de la región V3-V4 del gen 16S rRNA se utilizaron los primeros propuestos por Klindworth et al. (2012): S-D-Bact-0341-b785-a-A-21, 5´CCTACGGGNGGCWGCAG-3´ y S-D-Bact-0785-a-A-21, 5´-GACTACHVGGGTATCTAATCC-3´, con los cuales se generó un amplicón de ~460 pb. Se empleó el protocolo de PCR Illumina (2021a), suministrando 12.5 µL de MyTaqTM Ready Mix 1X (Bioline®), 1 µL de cada primer (10 uM), 5 µL de DNA (50 ng totales) y 5.5 µL de H2O grado molecular. Se utilizó el siguiente ciclo: 95 °C, por 3 minutos; 25 ciclos de 95 °C, durante 30 segundos, 55 °C por 30 segundos, 72 °C durante 30 segundos; 72 °C por 5 minutos, en un termociclador Labnet MultigeneTM Gradient PCR. Se efectuó la purificación de los amplicones con perlas Agentcourt® AMpure® XP al 0.8%. Posteriormente, los amplicones se etiquetaron utilizando el Nextera XT Index KitTM con ayuda del protocolo de Illumina (2021b), aplicando 25 µL de MyTaqTM Ready Mix 1X (Bioline®) 5 µL de cada primer (N7xx y S5xx), 5 µL de DNA y 10 µL de H2O, grado molecular; se empleó el siguiente ciclo: 95 °C por 3 minutos; 10 ciclos de 95 °C durante 30 segundos; 55 °C por 30 segundos, 72 °C durante 30 segundos; 72 °C por 5 minutos. Se efectuó la purificación de las bibliotecas con perlas Agentcourt® AMpure® XP al 1.2%. Para concluir, se llevó a cabo la cuantificación, normalización (equimolaridad) y secuenciación masiva de siguiente generación (MiSeq Illumina® de 2 × 250 lecturas de final pareado) empleando el protocolo para metagenómica 16S (Illumina, 2021a).
Análisis bioinformático
El análisis de las secuencias se efectuó en máquina virtual Oracle® VM VirtualBox® 6.1.34 en plataforma Linux, por medio del software bioinformático Quantitative Insights Into Microbial Ecology (QIIME) v.1.46 (Caporaso et al., 2010). El proceso comenzó con el ensamblaje de las secuencias, forward y reverse, de las muestras, utilizando el programa PEAR (Zhang, Kobert, Flouri y Stamatakis, 2014) con un traslape de 50 pb, una mínima longitud por lectura de 430 pb y máxima de 470 pb, con un criterio de calidad Q30 (una base errónea por cada 1000 bases) de acuerdo con lo recomendado por García-de la Peña et al. (2019). Posteriormente, se transformaron los archivos a formato FASTA y se realizó la eliminación de secuencias quiméricas de las muestras con USEARCH (Edgar, 2010).
Para la selección de las Unidades taxonómicas operacionales (OTUs) se aplicó el método UCLUST (Edgar, 2010) al 97% de similitud; se seleccionó al azar una secuencia representativa para cada OTU y se realizó la asignación taxonómica, utilizando como referencia la base de datos EzBioCloud (Yoon et al., 2017). Se elaboró el cuadro de OTUs en formato Biological observation matrix (Biom) (McDonald et al., 2012), se realizó la separación de los dominios y se filtraron los singletons (OTUs que solo tuvieron una observación) de acuerdo con Navas-Molina et al. (2013). Con los datos registrados se obtuvo la abundancia absoluta de OTUs y con esta información se graficó el número de secuencias por número de taxa para visualizar la profundidad de cobertura (curvas de tendencia asíntota); aplicando el programa PAST ver. 4.10 (Hammer, Harper y Ryan, 2001).
Para los análisis ecológicos y estadísticos se realizó un proceso de rarefacción aleatoria simple (Weiss et al., 2017) empleando como referencia el menor valor de las secuencias bacterianas, después de la eliminación de los singletons. Con el archivo estandarizado generado se calculó la diversidad alfa mediante los índices de Shannon y Simpson; para observar diferencia significativa en estos índices, entre las tres razas de maíz se emplearon pruebas no paramétricas de Kruskal-Wallis (P < 0.05) y pruebas de comparación de medias de índices para cada raza. Para la diversidad beta se utilizó el índice de Bray-Curtis (Beals, 1984), cuya matriz se empleó para efectuar una prueba PERMANOVA (P < 0.05) buscando diferencia significativa entre las microbiotas de las tres razas de maíz criollo. La visualización de la diversidad beta se realizó mediante una gráfica de análisis de coordenadas principales (PCoA) utilizando el programa Emperor (Vázquez-Baeza, Pirrung, Gonzalez y Knight, 2013). Las abundancias relativas para phylum, clase, orden, familia y género se analizaron en Microsoft® Excel® 2019. Los géneros con una abundancia ≥ 1% fueron representados en un mapa de calor; para la elaboración de dendrograma, de las muestras, se utilizó el método de conglomerados jerárquicos con medida euclidiana en el programa de libre acceso en línea Morpheus, Broad Institute. Finalmente, se llevó a cabo un análisis SIMPER y pruebas no paramétricas de Kruskal-Wallis a las rizobacterias encontradas entre las tres razas de maíz criollo, para cuantificar la contribución significativa de los taxa en las diferencias observadas con el índice de Bray-Curtis.
Se realizó la selección de las 100 mejores mazorcas, de las cuales se eligieron 10, basadas en el tamaño y peso, consideradas como las mejores en cada raza; para cada una se determinó su peso, utilizando una báscula digital marca Santul®, además, se registraron datos de número de hileras, número de granos por hilera, peso del grano y olote. Con estos datos se determinó el rendimiento, contemplando la densidad de plantas (67 000 plantas ha-1).
Resultados y Discusión
La media del total de secuencias obtenidas antes del ensamblaje fue de 35 607; las secuencias ensambladas fueron 33 142 y se descartaron 2465 secuencias. Se eliminaron en promedio 3608 secuencias quiméricas y la media de secuencias de calidad fue 29 523. Después de la asignación taxonómica y de la eliminación de los singletons se obtuvo un promedio de 21 134 secuencias bacterianas y finalmente se registraron 6853 OTUs bacterianos sin singletons (Cuadro 1). Las OTU en la raza C1 fueron 31 phyla, 88 clases, 171 órdenes, 380 familias, 977 géneros y 1059 especies; en la raza C2, 27 phyla, 86 clases, 164 órdenes, 364 familias, 973 géneros y 1040 especies; y la raza C3, 28 phyla, 84 clases, 163 órdenes, 351 familias, 937 géneros y 1012 especies. Se presentó una apropiada profundidad de cobertura referente al número de OTUs bacterianos de las muestras, dado que todas las curvas lograron la asíntota alrededor de las 13 000 secuencias (Figura 1).
Table 1: Information of the bacterial sequences analyzed in this study.
| Muestra | ST | SE | SD | QE | SC | SBSS | OTUs |
| C1-1 | 42 542 | 39 505 | 3 037 | 3162 | 36 336 | 26 136 | 7711 |
| C1-2 | 36 579 | 34 209 | 2 370 | 4050 | 30 148 | 21 937 | 7265 |
| C1-3 | 41 190 | 38 072 | 3 118 | 4048 | 34 007 | 24 301 | 7475 |
| C2-1 | 27 361 | 25 342 | 2 019 | 2739 | 22 599 | 14 529 | 5775 |
| C2-2 | 37 218 | 34 527 | 2 691 | 3986 | 30 532 | 22 146 | 7155 |
| C2-3 | 37 860 | 35 374 | 2 486 | 3721 | 31 634 | 23 475 | 7198 |
| C3-1 | 32 456 | 30 373 | 2 083 | 3609 | 26 753 | 18 891 | 6010 |
| C3-2 | 32 765 | 30 564 | 2 201 | 3827 | 26 730 | 19 150 | 6547 |
| C3-3 | 32 489 | 30 310 | 2 179 | 3334 | 26 964 | 19 639 | 6545 |
| Media | 35 607 | 33 142 | 2 465 | 3608 | 29 523 | 21 134 | 6853 |
C1 = criollo 1; C2 = criollo 2; C3 = criollo 3; 1 a 3 = repeticiones; ST = secuencias totales; SE = secuencias ensambladas; SD = secuencias descartadas; QE = quimeras eliminadas; SC = secuencias de calidad después de la eliminación de quimeras; SBSS = Secuencias bacterianas después de la eliminación de los singletons; OTUs = unidades taxonómicas operacionales.
C1 = creole 1; C2 = creole 2; C3 = creole 3; 1 to 3 = replicates; ST = total sequences; SE = assembled sequences; SD = discarded sequences; QE = chimeras removed; SC = quality sequences after chimeras removal; SBSS = bacterial sequences after singletons removal; OTUs = operational taxonomic Units.
En un estudio realizado por Peiffer et al. (2013a) demostraron que el genotipo Z. mays, spp., mays influye significativamente en la diversidad bacteriana alfa y beta. En cuanto a la diversidad alfa se observó que no existió diferencia significativa entre las tres razas de maíz criollo, con respecto a los índices de Shannon y Simpson (prueba de Kruskal-Wallis y medias de índices para cada raza de maíz). Es probable que la significancia en el índice de Simpson para la raza C3 se deba a los exudados de las raíces y al genotipo (Peiffer y Ley, 2013b), lo cual coincide con lo que afirma Nazir et al. (2016) sobre los exudados de las raíces, que éstas están estrechamente relacionadas al reclutamiento de los microorganismos como son las BPCV. Se observó diferencia significativa en la diversidad beta entre las tres razas de maíz criollo (PERMANOVA: pseudo-F = 1.178, p = 0.007). En la gráfica de coordenadas principales (Figura 2) se aprecia la agrupación compacta del grupo C2-1, C2-2 y C2-3; el grupo C1-1, C1-2 y C1-3 aparecen como un grupo más heterogéneo que el anterior, pero aun formando un grupo distinguible; por otra parte, las muestras C3-1, C3-2 y C3-3 registraron la mayor heterogeneidad individual y no formaron un grupo definido. En general, se aprecia que las tres razas de maíz criollo fueron distintas en cuanto a su diversidad bacteriana (Figura 2).
Investigaciones previas que han utilizado herramientas moleculares para estudiar la composición de la comunidad bacteriana del suelo in situ demostraron que, aunque las bacterias se han subdividido en más de 100 phyla, menos de 10 son abundantes en el suelo (Janssen, 2006). En el presente estudio se reportaron 31 phyla: los tres más abundantes para C1, C2 y C3, respectivamente, fueron Proteobacteria ( = 30, 33 y 35%), seguido de Actinobacteria ( = 19, 18 y 17%) y Acidobacteria ( = 13, 13 y 12%); mientras tanto, Gemmatimonadetes ( = 10, 11 y 10%), Chloroflexi ( = 9, 7 y 8%), Bacteroidetes ( = 5, 6 y 6%) y Planctomycetes ( = 5%) resultaron ser los menos frecuentes (Figuras 3a, 3b). Una tendencia similar fue reportada por Qaisrani et al. (2019), quienes concluyeron, a través de diversas investigaciones, que es claro que estos phyla suelen ser las bacterias más dominantes en la rizósfera del maíz. De la misma manera, Nemergut et al. (2011) destacan, además, que la abundancia relativa estimada de los principales phyla varía entre diferentes suelos o muestras.
Para las tres razas de maíz C1, C2, C3, se registraron 96 clases taxonómicas de las cuales Alphaproteobacteria ( = 18, 19 y 20%), Actinobacteria_c ( = 11%), Vicinamibacter_c ( = 8.8 y 7%), Gammaproteobacteria ( = 5.6 y 7%), Gemmatimonadetes_c ( = 6.6 y 5%), Longimicrobia ( = 5%) y Deltaproteobacteria ( = 4%) fueron las más abundantes. A nivel de orden se obtuvieron 185 taxa de los cuales Sphingomonadales ( = 7.8 y 9%), Vicinamibacter_o ( = 8.8 y 7%), Gemmatimonadales ( = 6.6 y 5%), Rhizobiales ( = 5%), Rhodospirillales ( = 4.5 y 4%), Myxococcales ( = 3%) y Gaiellales ( = 3%) mostraron los porcentajes más altos. En cuanto a familias se reportaron 414 taxa donde predominaron Shingomonadaceae (= 6.7 y 8%), Vicinamibacter_f ( = 7.7 y 6%), Gemmatimonadaceae ( = 4.3 y 3%), Gaiellaceae ( = 3%), GQ263235_f ( = 3%), Rhodospirillaceae ( = 3%) y Tepidisphaeraceae ( = 3.2 y 3%). Asimismo, se obtuvieron 1164 géneros de los cuales los más abundantes fueron: Sphingomonas ( = 6.6 y 7 %), Gaiella ( = 3%), PAC001846_g ( = 3.3 y 2%), PAC001874_g ( = 3.2 y 2%), Tepidisphaera ( = 3.2 y 3%), PAC000624_g ( = 2%) y Steroidobacter ( = 2%).
De las 1283 especies registradas solo 14 de ellas poseen nombre taxonómico binomial: Allorhizobium borbori, Amycolatopsis keratiniphila, Bacillus subtilis, Enterobacter cloacae, Goodfellowiella coeruleoviolacea, Kallotenue papyrolyticum, Kibdelosporangium aridum, Klebsiella pneumonia, Neochlamydia hartmannellae, Nonomuraea roseoviolacea, Parachlamydia acanthamoebae , Photorhabdus luminescens, Roseomonas gilardii y Streptomyces thermoviolaceus; sin embargo, ninguna de ellas presentó una abundancia relativa mayor al 1%. En la Figura 4 se presenta un mapa de calor de las rizobacterias a nivel género con una abundancia relativa mayor al 1% de cada muestra edáfica de las tres razas de maíz criollo.
Sphingomonas (Proteobacteria) fue el género más abundante en los tres maíces criollos estudiados (Figura 4); este género posee diversas cualidades entre las cuales se encuentran: a) ser en un buen candidato para utilizarse como fertilizante biológico dado sus aptitudes para fijar N2 (Sun, Zhang, Xu y Hu, 2010); b) ser una BPCV que solubiliza potasio (Etesami, Emami y Alikhani, 2017), y por ser una alfa-proteobacteria posee una sorprendente habilidad para degradar contaminantes ambientales recalcitrantes (He et al., 2017) como pesticidas (Manickam, Reddy, Saini y Shanker, 2008) o herbicidas (Keum, Lee y Kim, 2008), considerándose como una excelente opción para la biorremediación (Peng et al., 2008); y por último, c) algunas de sus cepas producen ciertos compuestos antibióticos que pueden ser responsables, de forma parcial, del control de algunas infecciones fúngicas que suelen infectar los maíces domesticados (Chagas, Pessotti, Caraballo y Pupo, 2018).
Gaiella fue el segundo género más abundante en las tres razas de maíz. A pesar de que sus funciones son desconocidas, al pertenecer al phylum Actinobacteria posee diversas aptitudes que corresponden a una BPCV (Yadav et al., 2018), favoreciendo la concentración de N, la materia orgánica y los micro y macroelementos esenciales en el suelo, así como el metabolismo del carbono y el desarrollo de las plantas (AbdElgawad et al., 2020); adicionalmente, es un phylum que suele ser persistente en suelos áridos (Santos-Medellin, Edwards, Liechty, Nguyen y Sundaresan, 2017) ya que puede sobrevivir en temperaturas mesófilas y algunos de sus miembros en temperaturas termófilas que alcanzan los 60 °C (Edwards, 1993), protegiendo a los cultivos de los efectos negativos de la sequía (Chukwuneme, Babalola, Kutu y Ojuederie, 2020).
Le prosiguen en abundancia los géneros no cultivables PAC001846_g y PAC001874_g, miembros de Acidobacteria, phylum que colabora con los miembros de Proteobacteria para degradar los polisacáridos fúngicos y vegetales (Lladó, Žifčáková, Větrovský, Eichlerová y Baldrian, 2016; Zimmermann, Gonzalez y Saiz, 2006), describiéndose como el que tiende a incrementar su abundancia cuando el pH del suelo decrece (Männistö, Tiirola y Häggblom, 2007; Jones et al., 2009; Dimitriu y Grayston, 2010).
Sphingomonas, Steroidobacter, Microvirga y EF125410_g pertenecen a Proteobacteria, el phylum más abundante en el estudio, relevante para el ciclo global del nitrógeno, del azufre y del carbono con el fin de asegurar los procesos de los ciclos biogeoquímicos sostenibles (Itävaara, Salavirta, Marjamaa y Ruskeeniemi, 2016); de la misma manera, tienen un papel fundamental en la degradación de la materia orgánica (Coates, Cole, Chakraborty, O’Connor y Achenbach, 2002; Chaudhry, Rehman, Mishra, Chauhan y Nautiyal, 2012; Lladó et al., 2016) para producir una gran cantidad de alcoholes aromáticos y oligosacáridos que suelen ser utilizados como fuente de carbono por las acidobacterias (Wei et al., 2017), aumentando la abundancia de estas últimas (PAC001846_g, PAC001874_g y PAC002019_g).
Tepidisphaera (Planctomycetes) es un phylum que se encuentra, comúnmente en las aguas termales terrestres y en diversos suelos de forma natural, inclusive encontrándose en diferentes sistemas productivos agrícolas e industriales (Kovaleva et al., 2019); en el presente estudio respecto a las tres razas de maíz, obtuvo el quinto lugar en abundancia relativa. Los géneros PAC000624_g, EF540396_g y “other” pertenecen a Gemmatimonadetes, phylum que suele representar aproximadamente del 1-2% de las bacterias del suelo a nivel mundial (Delgado-Baquerizo et al., 2018), siendo uno de los grupos menos estudiados; asimismo, al ser cosmopolitas en suelos, son capaz de acoplarse a un amplio abanico de nutrientes al tener un metabolismo versátil (Mujakic, Piwosz y Koblizek, 2022), con la propiedad de adaptarse a ambientes secos, incrementando su abundancia en suelos áridos y semiáridos (Acosta-Martínez, Dowd, Sun y Allen, 2008; Kim, Dungan y Crowley, 2008; Cary, McDonald, Barrett y Cowan, 2010; Neilson et al., 2017; Ahmed, Verma, Gupta, Mandhan y Chauhan, 2018; Zhao et al., 2020); esto tiene gran trascendencia debido a que la Comarca Lagunera se ubica dentro de la zona árida de México.
Por su parte, Flavisolibacter pertenece al phylum Bacteroidetes, este último caracterizado por estar involucrado en la degradación de moléculas orgánicas complejas de forma aeróbica, siendo copiótrofo, ya que tiende a incrementar su abundancia en el suelo después de que se incorpora carbono (Fierer, Bradford y Jackson, 2007; Eilers, Lauber, Knight y Fierer, 2010) y después de incrementar el pH del suelo (Lauber, Hamady, Knight y Fierer, 2009) obteniendo el décimo lugar en el presente estudio. Aunque Bacillus perteneciente a Firmicutes, es una rizobacteria solubilizadora de potasio o BPCV (Etesami et al., 2017) en el presente estudio se registró una abundancia baja.
El análisis de similitud porcentual (SIMPER) realizado para cuantificar la contribución de los taxa en términos de abundancia relativa y las diferencias observadas con el índice de Bray-Curtis se presentan en el Cuadro 2. Solo se incluyen las rizobacterias que obtuvieron un promedio de disimilitud > 0.10 y que obtuvieron un valor P < 0.05, explicando la mayor parte de las diferencias.
Table 2: Percentage similarity analysis (SIMPER) and Kruskal-Wallis test.
| Taxa | Disimilitud Promedio | Contrib. | Acumulado | Media C1 | Media C2 | Media C3 | H | g.l | p | |
| % | % | |||||||||
| Phylum | Proteobacteria | 1.668 | 26.110 | 26.110 | 0.303 | 0.329 | 0.351 | 5.956 | 2 | 0.050 |
| Clase | Phycisphaerae | 0.132 | 1.485 | 80.340 | 0.031 | 0.027 | 0.029 | 6.489 | 2 | 0.038 |
| Orden | Sphingomonadales | 0.696 | 5.676 | 5.676 | 0.068 | 0.078 | 0.088 | 5.956 | 2 | 0.050 |
| Micrococcales | 0.362 | 2.955 | 21.820 | 0.012 | 0.017 | 0.023 | 7.200 | 2 | 0.026† | |
| Phycisphaerales | 0.129 | 1.055 | 69.430 | 0.031 | 0.028 | 0.029 | 6.489 | 2 | 0.038 | |
| Familia | Moraxellaceae | 0.384 | 2.518 | 11.390 | 0.000 | 0.003 | 0.012 | 7.200 | 2 | 0.027 |
| Micrococcaceae | 0.278 | 1.823 | 21.770 | 0.007 | 0.011 | 0.015 | 7.200 | 2 | 0.027† | |
| Chitinophagaceae | 0.265 | 1.737 | 23.500 | 0.018 | 0.022 | 0.021 | 10.110 | 2 | 0.001 | |
| Sphingobacteriaceae | 0.194 | 1.275 | 35.650 | 0.002 | 0.004 | 0.007 | 6.489 | 2 | 0.038 | |
| Enterobacteriaceae | 0.119 | 0.783 | 58.210 | 0.003 | 0.002 | 0.004 | 6.489 | 2 | 0.037 | |
| Género | Acinetobacter | 0.385 | 1.899 | 4.761 | 0.000 | 0.002 | 0.012 | 7.200 | 2 | 0.027† |
| Chryseobacterium | 0.194 | 0.957 | 21.07 | 0.001 | 0.002 | 0.006 | 5.956 | 2 | 0.050 | |
| Sphingobacterium | 0.135 | 0.664 | 28.68 | 0.000 | 0.001 | 0.004 | 6.489 | 2 | 0.037 |
Contrib. % = porcentaje de contribución; H = valor Kruskal-Wallis; g.l. = grados de libertad; P = valor de significancia; †bacteria promotora del crecimiento vegetal.
Contrib. % = contribution percentage; H = Kruskal-Wallis value; g.l. = degrees of freedom; p = significance value; †plant growth promoting bacteria.
Entre las abundancias relativas de las razas estudiadas, Sphingomonadales miembro de Proteobacteria, y Chryseobacterium perteneciente al phylum Flavobacteria no registraron diferencia significativa (P < 0.05); Cuadro 2. Por el contrario, los taxa que presentaron diferencia significativa se describen a continuación: Enterobacteriaceae miembro de Proteobacteria; los miembros de Enterobacteriaceae se pueden encontrar en suelo, granos y plantas, muchas de sus especies causan tizón, marchitez y podredumbre blanda en el maíz, a menudo destruyendo cantidades sustanciales de la pared celular vegetal de las plantas (Brenner et al., 2005), obteniendo en abundancia de 0.2% entre las tres razas; Acinetobacter es miembro de la familia Moraxellaceae que suele encontrarse de forma natural en los suelos (Brenner et al., 2005); es un género BPCV que promueve el crecimiento vegetal produciendo fitohormonas, sideróforos y solubilizando el fosfato tricálcico insoluble (Sachdev, Nema, Dhakephalkar, Zinjarde y Chopade, 2010).
El orden Phycisphaerales es miembro de la clase Phycisphaera, phylum Planctomycetes; este último suele ser afectado por la distribución de oxígeno en el suelo ya que su abundancia es menor en el suelo cercano a la rizósfera óxica que en el suelo anóxico (Derakshani, Lukow y Liesack, 2001) registrando un 5% de abundancia en el presente estudio. La familia Micrococcaceae pertenece al orden Micrococcales y al phylum Actinobacteria, el cual se caracteriza por ser BPCV, capaz de inhibir patógenos, sintetizar sustancias reguladoras del desarrollo vegetal y la disponibilidad de minerales y nutrientes (Franco-Correa y Cahavarro, 2016).
El género Sphingobacterium pertenece a la familia Sphingobacteriaceae; este género, en determinados cultivos de gramíneas, atenúa la respuesta al estrés por sequía (Khan, Glielich, Phan, Redman y Doty, 2012; Comas, Becker, Cruz, Byrne y Dierig, 2013); algunas de sus especies tienen aptitudes para incrementar el crecimiento, además de que se han aislado de las raíces del maíz (Mehnaz, Weselowski y Lazarovits, 2007). Proteobacteria es uno de los phylum que suele ser predominante en la rizósfera de maíz; es de crecimiento rápido, ya que participa en la liberación del carbono lábil; no obstante, en los suelos no rizosféricos suele enriquecerse con los microbios oligotróficos de crecimiento lento, tales como Planctomycetes, Actidobacteria y Chlorofexi. En estudios previos se ha demostrado que las diversidades microbianas contribuyen a la funcionalidad del maíz (Mehta et al., 2021). La familia Chitinophagaceae es miembro de Bacteroidetes, phylum involucrado en la degradación de moléculas orgánicas; en el suelo degradan material vegetal y su abundancia relativa aumenta con la adición de carbono en la rizósfera (Fierer et al., 2007; Eilers et al., 2010).
En el Cuadro 3 se muestran los rendimientos en peso de grano y peso seco de cada raza de maíz criollo junto con una estimación aproximada en mega gramos por hectárea. Asimismo, se presentan los resultados del peso de grano evaluado con respecto a la selección de las 10 mejores mazorcas de cada tratamiento, el peso seco calculado (considerando tres plantas por tratamiento) y la estimación muestran también.
Cuadro 3: Rendimiento de peso en grano, peso seco y estimación aproximada de las toneladas de peso seco por hectárea.
Table 3: Grain weight yield, dry weight and approximate estimate of dry weight tons per hectare.
| Criollo | PG† | PS | Rendimiento PS | Rendimiento grano |
| g mazorca-1 | g planta-1 | - - - - - - - - - - Mg ha-1 - - - - - - - - - | ||
| C1 | 116.2±24.02 b | 84.23±12.38 b | 7.41±12.38 b | 7.78±24.02 b |
| C2 | 232.3±39.60 a | 120.6±10.44 a | 10.61±10.44 a | 15.56±39.6 a |
| C3 | 126.8±31.73 b | 83.57±4.65 b | 7.35±4.56 b | 8.49±31.73 b |
PG†: peso promedio de 10 mazorcas; PS = peso seco de plantas.
PG†: mean weight of 10 corncobs; PS = dry weight of plants.
La comparación múltiple de medias respecto a los datos obtenidos del peso de grano dio como resultado una diferencia notable entre la media del C2, destinado para grano, contra el resto de las razas. De acuerdo con los resultados de la estimación de producción en peso seco por hectárea, el C2 presentó la mayor cantidad en toneladas con respecto a las otras dos razas de maíz, mientras que los criollos 1 y 3 no presentaron una diferencia significativa entre ambas. En cuanto al peso seco estimado por hectárea, el C2 obtuvo la mayor producción (10.61 Mg ha-1), de la misma forma para el rendimiento en grano (15.56 Mg ha-1) con respecto a las otras dos razas (Cuadro 3).
Conclusiones
Entre las tres razas de maíz criollo, el C3, destinado para grano, registro la mayor abundancia relativa entre las bacterias con capacidad BPCV (Orden Micrococcales, Familia Micrococcales, género Acinetobacter y Sphingobacterium) que demostraron tener un promedio de disimilitud > 0.10 y un valor P < 0.05 en la prueba de Kruskal-Wallis, con una correlación positiva con el rendimiento de los cultivos. Por tanto, el genotipo de la raza C2 cumple funciones de doble propósito (grano y forraje) y puede ser considerado como candidato para preservar y evaluar, bajo condiciones diferentes, su rendimiento en grano; de igual manera C2 fue la raza criolla con mayor rendimiento estimado en el presente estudio con 15.57 toneladas por hectárea.
Declaración de Ética
No aplicable. En virtud de que no se utilizó ningún agente químico que pudiera dañar el ambiente del suelo, ni organismos transgénicos. Se cumple con las normas establecidas por Terra Latinoamericana.
Consentimiento para Publicación
Todos los autores han leído y aprobado el manuscrito final y han dado su consentimiento para su publicación.
Contribución de los Autores
Conceptualización de la investigación, diseño metodológico de la misma e interpretación de datos, redacción y revisión de manuscrito: J.V.A. Realización de trabajo experimental, análisis metodológico metagenómico, revisión y redacción de artículo: M.L.A. Trabajo de campo y rendimientos en grano de las razas de maíz: Y.D.C. Análisis metodológico metagenómico, revisión y redacción de artículo: E.M.M.M. Diseño, análisis de los resultados y la revisión del manuscrito: E.B.C., G.Z.C., y A.C.M. Revisión de datos y manuscrito: R.R.G.E. Revisión, supervisión de avances del proyecto, redacción y responsable de la publicación del manuscrito: A.M.R. Conceptualización de la investigación, diseño metodológico, análisis metagenómico, revisión de manuscrito: C.G.P.
