Introducción
En la actualidad, la agricultura enfrenta las modificaciones al ambiente por efecto del cambio climático, que, de manera particular, modifica la temperatura, la radiación y la precipitación pluvial, los cuales afectan negativamente la producción (Arce-Romero et al., 2018). Se considera que una alternativa para mitigar los efectos del cambio climático es el uso de la diversidad genética de las variedades nativas, debido a que están adaptadas a condiciones adversas, y comúnmente enfrentan diversos tipos de estrés, ocasionados por factores bióticos y abióticos (Rincón-Sánchez y Ruiz-Torres, 2017).
En el desarrollo de variedades mejoradas para incrementar el potencial de rendimiento, productividad y adaptación a las condiciones ambientales se requiere identificar fuentes de variación genética para ser incorporada en los procesos de selección. Los maíces nativos cuentan con información genética de importancia agronómica, como la capacidad de adaptación a condiciones adversas, precocidad, potencial de rendimiento y estabilidad en diferentes ambientes, los cuales pueden ser utilizados en programas de mejoramiento genético (González-Martínez et al., 2020). Además de la evaluación agronómica, se requiere una alternativa para identificar germoplasma con potencial bajo condiciones de estrés, Farooq et al. (2015) proponen realizar ensayos donde los genotipos se sometan a condiciones de salinidad.
En un estudio de la respuesta del sistema radical en maíz a condiciones de estrés por sequía Li et al. (2015) reportan variación genética de la raíz, la cual puede ser utilizada como criterio de selección en estado de plántula. En maíz, se pueden realizar pruebas masivas para identificar genotipos tolerantes a salinidad con base en el peso de las plántulas y tasa de crecimiento bajo condiciones de estrés salino (Farooq et al., 2015). En ensayos realizados durante las primeras etapas de crecimiento en maíz, Khan et al. (2003) encontraron que el crecimiento de las raíces en función de las concentraciones de sal es útil para seleccionar de forma preliminar genotipos tolerantes a salinidad.
Asimismo, un buen desarrollo del sistema radical en la emergencia y etapa temprana de desarrollo es de gran importancia para que el cultivo de maíz pueda establecerse en condiciones de estrés hídrico (Leach et al., 2011). Los objetivos de la presente investigación fueron analizar la producción de materia seca en etapas tempranas de desarrollo y selección de poblaciones de maíz bajo condiciones de estrés por salinidad. El trabajo se realizó bajo los supuestos que en las poblaciones de maíz existe variación genética y respuesta diferencial en la producción de materia seca que puede ser utilizada como criterio de selección temprana.
Materiales y métodos
El material genético consistió en 118 poblaciones de maíz clasificadas como raza Ratón procedentes de los estados de Coahuila (83), Nuevo León (20) y Tamaulipas (15), resguardados en el Centro de Conservación de Semillas Ortodoxas Región Norte (CCSONorte), ubicado en la Universidad Autónoma Agraria Antonio Narro (UAAAN). En el experimento se incluyeron dos híbridos comerciales como testigo: Rx715 de Asgrow y 6818 de Syngenta. El ensayo se estableció en el invernadero número dos de la UAAAN, Saltillo, Coahuila, México. El invernadero cuenta con control pasivo de la temperatura, la cual fue registrada cada cuatro horas utilizando un Hobo® data logger. La temperatura promedio registrada durante el desarrollo del experimento fue de 26.5 y 18.7 °C y humedad relativa promedio de 54.1 y 70.4%, durante el día y noche, respectivamente.
El estudio se realizó en dos ambientes contrastantes con y sin estrés por salinidad. En el ensayo se establecieron 120 genotipos (118 poblaciones y dos testigos) en los meses de junio, julio y agosto de 2021 bajo el diseño experimental de parcelas divididas en bloques completos al azar con cuatro repeticiones, en donde la parcela grande A fueron los ambientes (con/sin sal) y la parcela chica B, los genotipos. Por cuestiones de manejo del material experimental y toma de datos, cada bloque fue establecido en diferentes fechas de siembra, con separación de una semana y una duración del experimento de cinco semanas a partir de la fecha de siembra.
La unidad experimental (UE) consistió en macetas d de polietileno negro de 6 x 20 cm, con un sustrato compuesto de peat-moss, Perlita y Vermiculita en relación 2:1:1, respectivamente. Se sembraron tres semillas por maceta a tres centímetros de profundidad y posterior a la emergencia se eliminó una planta, dejando solo dos plantas por UE, sin fertilización adicional. Los riegos para el ambiente de estrés consistieron en una solución salina a base de NaCl (18.57 mS/cm) y aplicación de agua destilada para el ambiente sin estrés. En cada caso se aplicaron 100 ml por maceta con un intervalo de riego de siete días. Para apoyar a las plantas en la etapa de germinación y emergencia, se aplicó un riego de auxilio al tercer día después de la siembra.
Después de cada riego se realizaron lecturas de conductividad eléctrica de forma directa en el sustrato, con el medidor portátil de conductividad para suelo Hanna® HI993310. Al finalizar el ensayo, se realizó la estimación del contenido de clorofila (CLO) en unidades SPAD (Índice de contenido relativo de clorofila) con ayuda de un medidor portátil SPAD 501 Konica Minolta (Konica Minolta Inc, 2009). La muestra se realizó en la última hoja ligulada de una planta de la unidad experimental y se realizaron tres lecturas, registrando el valor promedio. Una vez transcurridas las cinco semanas a partir de la siembra, las raíces fueron extraídas del sustrato con cuidado.
Cada planta fue seccionada en dos partes (raíz y vástago) y luego se depositaron en bolsas de papel de estraza, para someterlas a secado en una estufa a 70 °C por 48 h. Después, se pesaron las muestras en una balanza analítica marca AMD HR-200 y el resultado fue expresado como peso seco de raíz (PSR) y peso seco de vástago (PSV) en mg planta-1.
Se realizó análisis de varianza de acuerdo con el diseño experimental utilizando y una comparación múltiple de medias de Tukey (p≤ 0.05) con el procedimiento PROC GLM de SAS (SAS Institute Inc, 2018). Se computó el valor medio del peso seco total (PST) con las fracciones del peso seco de raíz (PSR) y peso seco de vástago (PSV). Con los valores medios se realizó un análisis de correlación simple de Pearson con el procedimiento PROC CORR de SAS (SAS Institute Inc, 2018). Los ambientes, genotipos y la interacción genotipo x ambiente fueron considerados como efectos fijos, el resto de los efectos en el modelo como aleatorios. Se calculó un valor de decisión para identificar genotipos superiores, definido por la media más uno y dos veces el error estándar de la media (𝜇 + ɛɛ; 𝜇 + 2ɛɛ).
Finalmente, se utilizó un modelo aleatorio para estimar los componentes de varianza con el procedimiento PROC MIXED de SAS (SAS Institute Inc, 2018). Con los componentes de varianza se estimó el coeficiente de determinación genética o heredabilidad en sentido amplio para los caracteres en estudio (Hallauer et al., 2010). Las medias de genotipos x caracteres en ambientes fueron utilizadas para realizar el análisis de dispersión gráfica de estos, basado en el análisis de componentes principales (ACP), el cual fue obtenido con el modelo GGEbiplot (Yan, 2014).
Resultados y discusión
Los cuadrados medios del análisis de varianza y la estimación de los componentes de varianza para los caracteres peso seco de raíz (PSR), peso seco de vástago (PSV) y contenido de clorofila (CLO), obtenidas en los ensayos de invernadero en 2021 se presentan en el (Cuadro 1).
FV | GL | Peso seco de raíz (mg planta-1) | Peso seco de vástago (mg planta-1) | >Contenido de clorofila (SPAD) | |||||
CM | σ 2 | CM | σ 2 | CM | σ 2 | ||||
Bloques (BLO) | 3 | 55 258.4 | 0 | 424 974.9 | 0 | 2 075.8 | 6.3 | ||
Ambientes (AMB) | 1 | 447 912.7 | 732.1 | 17 893 944.1** | 36 409 | 10 881* | 21.5 | ||
BLO × AMB | 3 | 135 315.3** | 775.8 | 475 056.7** | 3 652.3 | 589.5** | 4.8 | ||
Genotipos (GEN) | 119 | 5 293.1** | 397.7 | 40 057.8** | 3 150.6 | 11 | 0 | ||
AMB × GEN | 119 | 2 143.7 | 78.8 | 15 098.2 | 295 | 12.4 | 0 | ||
Error | 712 | 1 826.6 | 1 827.9 | 13 994.1 | 13 989 | 11.8 | 11.8 | ||
CV (%) | 29.6 | 27 | 11.7 | ||||||
H2 | 0.6 | 0.62 | 0 |
*, **= significativo al 0.05 y 0.01 niveles de probabilidad, respectivamente; FV= fuentes de variación; GL= grados de libertad; CM= cuadrados medios; CV= coeficiente de variación; H2= coeficiente de determinación genética o heredabilidad en sentido amplio.
En los ambientes de evaluación (con y sin aplicación salina) se encontraron diferencias estadísticas en PSV y CLO, en tanto que en PSR no hubo evidencia diferencias estadísticas. En los genotipos (GEN) (poblaciones y testigos) se encontraron diferencias estadísticas (p≤ 0.01) para las variables de peso seco, excepto en CLO, lo cual se debe a la variación existente en el material experimental en estudio. No se encontraron diferencias estadísticas en la interacción de genotipos × ambientes en ninguna de las variables en estudio, lo cual indica que los genotipos mostraron un comportamiento similar en los dos ambientes de evaluación.
Considerando que en el ensayo no se aplicó fertilización, Bänziger et al. (2012) comparten un punto de vista que permite analizar el fenómeno de la relación raíz-vástago y mencionan que el nitrógeno es el nutriente que la planta requiere en más grandes cantidades y éste está estrechamente relacionado con la producción de biomasa vegetativa y el contenido de clorofila; además, mencionan que en condiciones de estrés por deficiencia de nitrógeno, las plantas favorecen el crecimiento de raíces más que el crecimiento de vástago y, por consecuencia, la proporción raíz: vástago aumenta.
En el Cuadro 2 se presentan los valores medios de los ambientes de evaluación (con y sin sal), identificados como las condiciones de estrés y testigo, respectivamente, que representan los ensayos realizados en invernadero para las variables en estudio. De acuerdo con este, las condiciones de estrés salino (estrés) tuvieron un efecto en PSR y PSV con respecto a las condiciones sin estrés (testigo), lo cual se manifestó en una reducción de 25.9 y 47.5%, respectivamente. En el caso del contenido de clorofila (CLO), la respuesta fue opuesta, pues se encontró un incremento promedio de 26% en las condiciones de estrés con respecto al testigoBALANZA. Lo anterior debido a que los materiales que no estuvieron expuestos a sal, mostraron una mayor producción de biomasa, lo cual provocó un agotamiento de nutrientes, debido a que ellos nunca recibieron ningún tipo de fertilización.
Ambientes | Peso seco de raíz (mg planta-1) | Peso seco de vástago (mg planta-1) | Contenido de clorofila (SPAD) |
Testigo (sin sal) | 165.6 a† | 574.66 a | 25.99 b |
Estrés (con sal) | 122.7 a | 301.75 b | 32.74 a |
Media | 144.15 | 438.21 | 29.36 |
Tukey (α= 0.05) | 75.65 | 141.74 | 5 |
†= valores con la misma letra en la columna no son estadísticamente diferentes.
A este respecto, Correndo y García (2012) mencionan que existe un fenómeno en donde el aumento de producción de materia seca provoca una dilución de nutrientes en los tejidos cuando no existe aportación de nutrientes. Durante el desarrollo del experimento, en las condiciones de estrés salino, las plantas de maíz fueron expuestas a una conductividad eléctrica promedio de 3.45 mS cm-1, con una amplitud entre bloques (fechas de siembra) de 3.13 a 3.88 mS cm-1. Leland y Mass (1999), refieren que el cultivo de maíz comienza a ser afectado cuando las conductividades eléctricas son superiores a 1.7 mS cm-1. La raíz es el primer órgano expuesto al estrés por sal, los brotes subsecuentes son más sensibles en la respuesta y expresión en relación con las raíces en el cultivo del maíz (Iqbal et al., 2020).
El coeficiente de determinación genética o heredabilidad en sentido amplio (fracción de la varianza fenotípica debida a los efectos genéticos aditivos y dominancia) (Falconer y Mackay, 1996) resalta la importancia de los genotipos con relación a las diferencias fenotípicas para un carácter de interés en el proceso de selección. En este estudio, se estimó el coeficiente de heredabilidad en sentido amplio obteniendo valores de 0.6 y 0.62 para los caracteres peso seco de raíz (PSR) y peso seco de vástago (PSV), respectivamente.
En el caso del contenido de clorofila (CLO), la mayor parte de la variación se debió a efectos ambientales en relación con los efectos genéticos, por lo que este carácter, bajo las condiciones del presente estudio, no es de utilidad para efectos de selección en la etapa temprana de desarrollo. En un estudio de vigor híbrido para tolerancia a sequía, Ali et al. (2013) reportan valores de heredabilidad de 0.9 para el peso seco de raíz y 0.76 peso seco de vástago. También, en la evaluación de cuatro genotipos de maíz para tolerancia a sequía, Haseeb et al. (2020) reportan heredabilidades de 0.89 y 0.9 para el peso seco de raíz y peso seso de vástago, respectivamente. Sin embargo, en la evaluación de líneas endogámicas bajo condiciones de estrés hídrico, Masood et al. (2020) reportan valores de heredabilidad para PSR y PSV de 0.49 y 0.47, y de 0.34 y 0.46 para las condiciones de ambientes normal y de estrés, respectivamente.
En la evaluación de líneas endogámicas de maíz bajo diferentes niveles de nitrógeno, Abdel-Ghani et al. (2013), reportan valores estimados de heredabilidad de 89.6 y 0.83 para el peso seso de raíz y peso seco de vástago, respectivamente. El análisis de los resultados en el presente estudio indica la importancia de los caracteres de PSR y PSV en la selección de genotipos de maíz en etapas tempranas de desarrollo. Debido a la ausencia de interacción genotipo × ambiente, se computó el valor medio del peso seco total (PST) con las fracciones del peso seco de raíz (PSR) y peso seco de vástago (PSV). Con los valores medios del PSR, PSV y contenido de clorofila (CLO) en los ambientes de evaluación, se realizó un análisis de correlación simple de Pearson. El PST correlacionó positiva y significativamente (p≤ 0.01) con las fracciones PSR y PSV en las condiciones ambientales contrastantes (estrés con salinidad y testigo).
Sin embargo, no se encontró asociación entre los caracteres de peso seco con CLO en los ambientes de evaluación. Las asociaciones entre los caracteres en estudio se presentan de manera gráfica más adelante. En el Cuadro 3 se presentan las medias de peso seco de los 25 genotipos superiores identificados con base en el peso seco de la planta en los ensayos en invernadero en 2021.
Genotipos | Grupos‡ | PST | PSR-E | PSV-E | PSR-T | PSV-T | |||||
Poblaciones† | |||||||||||
C-25 | C | 762.06 | ** | 147.13 | * | 440.13 | ** | 190.63 | 746.25 | ** | |
C-14 | A | 739.94 | ** | 113.75 | 297.88 | 265.63 | ** | 802.63 | ** | ||
C-77 | A | 726.44 | ** | 105.88 | 374 | * | 220 | * | 753 | ** | |
C-79 | A | 726.06 | ** | 152 | * | 324.5 | 203.38 | * | 772.25 | ** | |
C-55 | A | 724.5 | ** | 117.75 | 308.75 | 202.88 | * | 819.63 | ** | ||
C-8 | B | 704.94 | ** | 180 | ** | 417.75 | ** | 182.13 | 630 | ||
C-23 | C | 703.81 | ** | 96.5 | 307 | 239.38 | ** | 764.75 | ** | ||
T-106 | C | 698.19 | ** | 153.5 | * | 401.13 | * | 170 | 671.75 | * | |
C-17 | A | 696.75 | ** | 143.25 | 341.13 | 166 | 743.13 | ** | |||
C-19 | B | 696.25 | ** | 165.88 | ** | 386.5 | * | 199.5 | * | 640.63 | |
N-36 | A | 695.94 | ** | 111.38 | 346 | 207.13 | * | 727.38 | * | ||
C-24 | B | 692.31 | * | 167.5 | ** | 477.5 | ** | 172.75 | 566.88 | ||
N-50 | A | 692.19 | * | 156 | * | 302.38 | 183.75 | 742.25 | ** | ||
T-103 | A | 691 | * | 162.63 | * | 403.38 | * | 234.5 | ** | 581.5 | |
T-104 | B | 686.13 | * | 194 | ** | 363.13 | * | 182.88 | 632.25 | ||
C-4 | C | 679.56 | * | 152.75 | * | 376.75 | * | 185.38 | 644.25 | ||
C-11 | A | 678.94 | * | 123.63 | 300.88 | 178.13 | 755.25 | ** | |||
T-99 | C | 677.19 | * | 170.38 | ** | 360.13 | * | 228.88 | ** | 595 | |
N-42 | C | 676.88 | * | 114.13 | 364.13 | * | 185.25 | 690.25 | * | ||
C-16 | C | 676.75 | * | 126.38 | 405.75 | ** | 180 | 641.38 | |||
C-30 | C | 676.5 | * | 123.25 | 360.38 | * | 196.13 | * | 673.25 | * | |
C-21 | B | 674.19 | * | 145.63 | * | 368.13 | * | 187.38 | 647.25 | ||
C-31 | C | 671.06 | * | 144.63 | * | 381.5 | * | 190.13 | 625.88 | ||
T-109 | B | 669.75 | * | 136.5 | 387.5 | * | 180.25 | 635.25 | |||
C-29 | A | 669.38 | * | 119.25 | 292 | 241.13 | ** | 686.38 | * | ||
Testigos | |||||||||||
6818 | 549.31 | 131.13 | 288.25 | 157.5 | 521.75 | ||||||
Rx715 | 296.25 | 54.25 | 112.63 | 95.88 | 329.75 | ||||||
Media | 581.92 | 122.49 | 301.29 | 165.6 | 574.66 | ||||||
EE | 55.9 | 21.64 | 52.12 | 28.92 | 78.19 |
*, **= selección con base μ+EE y μ+2EE, respectivamente; PST, PSR y PSV = peso seco total, peso seco de raíz y peso seco de vástago en mg planta-1; CLO= contenido de clorofila; E y T= estrés por salinidad y testigo (sin sal); †= primer dígito indica la procedencia; C= Coahuila, N= Nuevo León, T= Tamaulipas, seguido del número de entrada; ‡grupos: A= selección en ambiente favorable; B= selección en condiciones de estrés; C= selección promedio; EE= error estándar de la media.
En la producción promedio total de materia seca (PST) se encontró variación entre los genotipos evaluados, con una amplitud de 296.25 a 762.06, con una media de 581.92 mg planta-1, por lo que es factible de realizar selección (Cuadro 3). En el Cuadro 3 sólo se presentan los valores medios de los 25 genotipos superiores, donde, de manera particular, se puede dilucidar la variación de estos, con base en las fracciones de peso seco de raíz (PSR) y vástago (PSV) en los ambientes de la evaluación (condiciones de estrés salino y testigo).
Es importante resaltar que, de los 120 genotipos utilizados sólo sobresalen por su comportamiento, en PSR y PSV las poblaciones de la raza Ratón; es decir, ninguno de los testigos utilizados (híbridos), mostraron producción de materia seca significativa en PST o la expresión en las condiciones ambientales contrastantes. El análisis del comportamiento de los genotipos en los ambientes de evaluación y la asociación de los caracteres se presenta de manera gráfica (Figura 1), para explorar la interacción de los genotipos × caracteres en ambientes, obtenida a partir del análisis de componentes principales utilizando el modelo GGEbiplot (Yan, 2014).
En la Figura 1 se consideran los dos primeros componentes principales (PC1 y PC2), los cuales en conjunto explican 66.6% de la variación total acumulada en los valores medios de los genotipos y los caracteres × ambientes de evaluación. En el modelo de dispersión gráfica GGEbiplot se analiza la interacción de factores, en este caso, los caracteres en los ambientes de evaluación se usan como probadores de los genotipos en estudio, por lo tanto, se analiza la interacción genotipos × caracteres en ambientes.
En la Figura 1, el coseno del ángulo que forman los vectores de la combinación de caracteres en ambientes que se producen a partir del origen se aproxima al coeficiente de correlación entre ellos, un ángulo menor a 30º indica una relación muy estrecha, en tanto que un ángulo aproximadamente de 90º, una relación independiente, un ángulo de 180º estará representando una asociación negativa (Yan, 2014). Similarmente, cada genotipo proyecta un vector a partir del origen y, por lo tanto, la relación con los caracteres dependerá de la ubicación en el cuadrante.
De acuerdo con la proyección y ángulos de los vectores, los caracteres del peso seco están positivamente correlacionados (PSR y PSV), con una asociación muy estrecha dentro de cada uno de los ambientes. En este sentido, la relación del contenido de clorofila (CLO) con PSR y PSV en ambientes es relativamente baja, y con el PST, una correlación independiente. En el eje de la abscisa (PC1) de la (Figura 1) se proyecta gráficamente las medias del peso seco de los genotipos con la amplitud señalada en el (Cuadro 3), donde Rx715 obtuvo el valor mínimo de 296.25 mg planta-1 y el genotipo C-25 el valor máximo de 762.06 mg planta-1.
En el caso del eje de la ordenada (PC2) representa la interacción de los genotipos con los caracteres en ambientes, que de acuerdo con el modelo de interacción GGEbiplot, sufre una inclinación para efectos de la interpretación. Se puede dilucidar en la Figura 1 la dispersión amplia de los genotipos, como interpretación de la variación genética de estos con la producción de materia seca representada por los caracteres de este estudio. Asimismo, la Figura 1 es de utilidad ya que es posible identificar genotipos con características particulares de acuerdo con los patrones de variación en el estudio.
De esta manera, se puede identificar a los genotipos sobresalientes en el ambiente favorable (testigo) identificado como grupo A, aquellos sobresalientes en el ambiente de estrés por salinidad (grupo B), y los genotipos con un comportamiento promedio a través de ambientes, identificados con el grupo C. Por lo tanto, de los 25 genotipos superiores (Cuadro 3), se identificaron 10 en el ambiente favorable (testigo), seis en el ambiente de estrés (tolerantes a salinidad) y nueve con un comportamiento promedio a través de los ambientes de evaluación.
El peso seco de raíz (PSR) y el peso seco de vástago (PSV) con y sin estrés por salinidad fueron considerados en la evaluación temprana para determinar diferencias en el potencial fisiológico y como alternativa como criterio de selección en el cultivo de maíz con tolerancia a la salinidad como ha sido señalado por diversos autores (Giaveno et al., 2007; Marcos-Filho, 2015). Estos caracteres han sido estudiados para identificar genotipos bajo condiciones de estrés por sequía (Li et al., 2015) y con tolerancia a salinidad (Khan et al., 2003; Farooq et al., 2015), los cuales, se ha sugerido, pueden ser utilizados como criterios de selección en etapas tempranas de desarrollo en maíz.
Conclusiones
Se determina que existe variación genética en las poblaciones estudiadas de la raza Ratón para el peso seco de raíz y de vástago, con una estimación de heredabilidad en sentido amplio de 0.6 y 0.62, respectivamente, lo que permite realizar selección de genotipos en la etapa temprana de desarrollo. De los 25 genotipos superiores se identificaron 10 en el ambiente favorable (testigo), seis en el ambiente de estrés (tolerantes a salinidad) y nueve con un comportamiento promedio a través de ambientes.