text new page (beta)
Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Similars in
SciELO 
Permalink
Introducción
Los cereales, como el maíz (Zea mays L.), son el principal alimento y fuente de energía para los habitantes de la mayoría de los países en desarrollo. Cerca del 60 % de las calorías se derivan directamente de ellos y en los países más pobres se llega hasta el 80 % (Awika, 2011). Entre los tres principales cereales (maíz, trigo y arroz), el maíz contribuye con más del 40 % de la producción mundial con 1,038 millones de toneladas (FAOSTAT, 2015); además, por su alto potencial de rendimiento, es reconocido como un cultivo que puede contribuir en alto grado a la seguridad alimentaria (Almeida et al., 2014).
Se ha comprobado que el cambio climático se está presentando con un aumento en los fenómenos meteorológicos extremos como la temperatura y la variabilidad de las lluvias (Olivera, 2013). En el informe del Panel Intergubernamental del Cambio Climático se menciona que las sequías se intensificarán durante este siglo en algunas zonas de Europa Central, América Central, México, Noreste de Brasil y Sur de África debido a la reducción de las precipitaciones y al aumento de la evapotranspiración (IPCC, 2012). La escasez de agua es el factor más limitante del medio ambiente para la productividad del maíz (Inzunza-Ibarra et al., 2018) en las regiones tropicales y subtropicales (Messmer et al., 2011). Se estima que el rendimiento se reduce hasta en 39 % cuando hay una escasez de agua de 40 % aproximadamente (Daryanto et al., 2016).
En la historia del mejoramiento por hibridación se ha demostrado que la tolerancia al estrés está más asociada a la ganancia genética que a la heterosis en los genotipos más recientes (Duvick, 2005; Reyes et al., 2015) debido a que se ha mejorado la eficiencia del uso del agua y el índice de cosecha. El mejoramiento genético por selección para la tolerancia a la sequía mediante rendimiento de grano es ineficiente debido a la baja heredabilidad de ambos caracteres, por lo que se sugiere el uso de variables secundarias, que proporcionan información de la forma en que el rendimiento cambia en condiciones de sequía. Idealmente, un carácter secundario debe de estar correlacionado con el rendimiento de grano en condiciones de sequía, con variación genética y altamente heredable; también debe ser fácil, barato y rápido de observar o medir de forma no destructiva, estable durante el periodo de medición y no asociado con la pérdida de rendimiento en condiciones sin estrés (Lafitte et al., 2003).
Debido a que la tolerancia a la sequía es de herencia cuantitativa, controlada por muchos genes, resulta apropiado apoyarse en técnicas moleculares para identificar los principales procesos fisiológicos asociados con la mejora del rendimiento y la determinación de las asociaciones gen-fenotipo que permitan aumentar la eficiencia de la selección. El mejoramiento asistido por marcadores moleculares es una herramienta que agiliza el avance genético (Ribaut y Ragot, 2007), apoyado por loci de carácter cuantitativo. Los QTLs permite conocer la importancia de una región del genoma asociada a un carácter. En maíz se han descrito varios QTLs asociados a condiciones de estrés hídrico, como intervalo de floración, enrollamiento de la hoja, senescencia, vigor, entre otros (Almeida et al., 2014; Lu et al., 2012; Messmer et al., 2011; Trachsel et al., 2016).
La detección de QTLs en grupos de líneas puede ser no del todo buena debido a la carga genética para el caso del maíz, por lo que una solución es la evaluación de éstas en combinación con un probador (Frascaroli et al., 2009; Ramírez-Díaz et al., 2019) y así detectar QTLs consistentes. El objetivo de este trabajo fue identificar QTLs para la tolerancia a sequía en una población de maíz de trópico. La consecución de dicho objetivo incluyó i) estimación de la heredabilidad y la relación de los rasgos secundarios con el rendimiento de grano bajo ambientes de sequía y riego, ii) identificación de las regiones genómicas que influyen en los caracteres secundarios de los regímenes hídricos, y iii) detección de regiones del genoma con QTLs clave para posteriores trabajos de introgresión asistida por marcadores de tolerancia a la sequía en maíz tropical.
Materiales y métodos
Material vegetal
El material utilizado fue un grupo de 166 líneas de maíz biparentales obtenidas por la técnica de dobles haploides a partir de las líneas LPSC7F64 (La Posta Sequía C7 F64-2-6-2-2-B-B-B) y CML495, además de 196 cruzas de prueba obtenidas utilizando la línea CML494. La LPSC7F64 es una línea que fue derivada de la población La Posta Sequía, de grano blanco dentado del patrón heterótico B; es un sintético de la raza Tuxpeño, adaptado a la condición del trópico bajo, desarrollado mediante siete ciclos de selección recurrente bajo condiciones de sequía. La línea CML495 es de grano blanco cristalino desarrollada para condiciones de trópicos bajos del patrón heterótico A. El probador CML494 es una línea de grano blanco dentado obtenido de la población La Posta Sequía Ciclo 4 Familia 7, de grupo heterótico AB (http://apps.cimmyt.org).
Sitios experimentales
El experimento se estableció en la estación experimental del Centro International de Mejoramiento de Maíz y Trigo (CIMMYT) en Tlaltizapán, Morelos (18º 41’ N, 99º 08’ O, 940 msnm) e Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP) de Iguala, Guerrero (18º 20’ N, 99º 30’ O, 732 msnm). La siembra se realizó durante la temporada de invierno seco 2013 y 2014 cuando normalmente no hay lluvias y se depende de la irrigación.
Diseño y unidad experimental
Se utilizó un diseño experimental α-látice (0,1) con bloques incompletos y dos repeticiones. La unidad experimental fue de un surco de 4.5 m y 0.75 m entre surcos, un espaciamiento entre plantas de 0.2 m, para una densidad de 66.6 plantas m-2.
Manejo experimental
Se realizaron prácticas culturales convencionales en las parcelas de ensayo y los riegos fueron mediante goteo. El estrés por sequía se indujo de acuerdo con el protocolo de CIMMYT (Bänziger et al., 2012) a los 12-15 días antes de la antesis (~190 UC). Para el caso de un estrés hídrico severo indicado por un intervalo de floración (IF) mayor de 5 d en promedio del ensayo, el riego se aplicó 7 d después de finalizada la antesis, mientras que un segundo riego se aplicó tres semanas después de la finalización de la antesis. En el caso de estrés hídrico moderado, se aplicó un solo riego durante el periodo de llenado del grano dos semanas después de la finalización de la antesis. El monitoreo de la humedad se realizó tres veces por semana con una sonda de humedad del suelo PR2 /6 (Delta-T Devices, Cambridge, Reino Unido) a profundidades de 10, 20, 30, 40, 60 y 100 cm para la programación de riegos en el tratamiento de estrés por sequía.
Datos fenotípicos
Las variables fenotípicas se obtuvieron a lo largo del ciclo del cultivo. El índice de vegetación de diferencia normalizada (NDVI) se capturó a las dos, cuatro, cinco y seis semanas después de la siembra con un RT-Greenseeker (Trimble, Ukia, CA, EUA). El NDVI fue medido en cada unidad experimental a una altura de 0.8 m caminando con el sensor en medio de los surcos. El NDVI fue calculado de la siguiente manera:
donde: IRC es la reflectividad en el infrarrojo cercano y R es la reflectividad en el rojo.
Los días a antesis (DA) se obtuvieron con el número de días desde la siembra hasta el 50 % de las plantas con dehiscencia de polen. El intervalo de floración fue medido por el intervalo de tiempo en días entre la floración masculina y la femenina. La altura de planta (AP) se midió en cm al momento de la floración desde el suelo hasta la última hoja. La senescencia se registró a las dos, cuatro y seis semanas después de la floración mediante una calificación visual con una escala de 1 (sin senescencia) a 9 (senescencia completa). Para NDVI y senescencia, el área bajo la curva (AUC) se calculó mediante la integración de una función polinomial de segundo grado ajustada a mediciones individuales tomadas antes (por NDVI) o después de la floración (por senescencia). Se adicionó un índice que integra el vigor temprano (NDVI), vigor final (altura de planta) y senescencia, denominado COMBO, que se calculó de la siguiente manera:
donde: AP es la altura de planta a floración, AUCNDVI es el área bajo la curva de NDVI y AUCSEN es el área bajo la curva de la senescencia.
El rendimiento de grano se calculó con el peso de grano ajustado a una humedad de 12.5 %.
Análisis estadístico
Los componentes de varianza fueron estimados con máxima verosimilitud lineal restringida mediante META-R 2.5 (Alvarado et al., 2015). El modelo estándar lineal mixto para la variable respuesta estuvo representado por:
En el que µ es la media general, G i es el efecto del i-ésimo genotipo (i = 166 y 196 líneas y cruzas de prueba), E j es el efecto del j-ésimo ambiente (CP_R: j = 2, CP_S: j = 4, Líneas DH_R j = 2 y Líneas DH_S j = 2), GE ij es el efecto de interacción del i-ésimo genotipo con el j-ésimo ambiente, R(E) kj es el efecto de la k-ésima repetición dentro del j-ésimo ambiente, BI(RE) bkj es el efecto de los bloques incompletos dentro de la k-ésima repetición en el j-ésimo ambiente. El efecto ambiente no fue considerado aleatorio. Los mejores predictores lineales insesgados (BLUP) se calcularon para las líneas y cruzas de prueba en las dos condiciones de riego por separado.
La heredabilidad en sentido amplio fue estimada como:
donde: σ 2 g es la varianza genotípica, σ 2 g x e es la varianza de la interacción genotipo × ambiente, σ 2 e es la varianza residual, e es el número de ambientes y r es el número de repeticiones por ambiente.
La correlación genética (r g ) entre las características se calculó con la siguiente ecuación:
donde: σ g(ii') es la covarianza genotípica entre cada par de caracteres y σ g(i) σ g(i') es el producto de las desviaciones estándar genotípicas de cada carácter.
Las correlaciones fenotípicas entre las características se calcularon con base en los BLUP como los coeficientes de correlación de Pearson.
Genotipificación y construcción de mapa de ligamiento
Se tomaron 12 muestras de hoja por cada una de las líneas usadas en el estudio. Para el aislamiento del ADN se usó el procedimiento de CTAB (CIMMYT, 2005). El ADN de todas las muestras se envió a Cornell University Biotechnology Resource Center (Ithaca, NY, EUA) para genotipeo por secuenciación (GBS), donde se aplicó el protocolo de Elshire et al. (2011). El ADN genómico se digirió con la enzima de restricción ApaK1. Las bibliotecas de GBS fueron construidas en 96-plex y secuenciadas en un equipo Illumina HiSeq 2000 (Illumina, Inc., San Diego, California, EUA) (Glaubitz et al., 2014; Zhang et al., 2015). Los marcadores de polimorfismo de nucleótido individual (SNP) se filtraron e imputaron en el programa TASSEL GBS y se usó como referencia el genoma de la línea B73 (Glaubitz et al., 2014). El mapa de ligamiento se construyó en ICIMv4 con 476 SNP; para estimar la distancia del mapa se usó la función de Kosambi (Kosambi, 1943).
Análisis de QTLs
El mapeo de intervalo compuesto fue realizado con ICIMv4 (Li et al., 2007). Se incluyeron 174 genotipos en el mapa con una convergencia de 1 cM. Un QTL con valor LOD mayor o igual que 2.5 fue considerado como significativo (P ≤ 0.05). Un QTL para diferentes caracteres fue declarado común cuando sus posiciones con altas puntuaciones LOD (pico) se ubicaban en la misma posición. El efecto aditivo se obtuvo con ICIMv4. El signo de efectos aditivos fue usado para identificar el origen del alelo favorable. Aditividad positiva indica un efecto de sustitución proveniente del progenitor recurrente (CML495) y aditividad negativa indica el efecto de sustitución proviene del donador (LPSC7F64). R2 representa el porcentaje de la varianza fenotípica explicada por todos los QTLs detectados para cada carácter (Zeng et al., 1999).
Resultados y discusión
Las estimaciones de medias, componentes de varianza y heredabilidad de las líneas DH se muestran en el Cuadro 1 y de las cruzas de prueba en el Cuadro 2. El rendimiento promedio de la población de líneas DH bajo riego y sequia fue de 1.9 y 0.7 t ha-1, mientras que en las cruzas de prueba fue de 4.11 y 3.11 t ha-1 para riego y sequía, respectivamente. La heredabilidad fue de 0.40 y 0.56 para líneas bajo riego y sequía, respectivamente (Figura 1), mientras que en cruzas de prueba (CP) ésta fue de 0.61 y 0.28 para riego y sequía, respectivamente. La reducción del rendimiento por efecto de la sequía fue de 63.16 % para las líneas DH y de 24.33 % para las cruzas de prueba.
Cuadro 1 Estimación de medias, componentes de la varianza y heredabilidad del rendimiento y caracteres secundarios de líneas DH bajo riego (R) y sequía (S).
| Característica | Condición | Media | σ2 g | σ2 g x e | σ2 e | H 2 | CV (%) |
|---|---|---|---|---|---|---|---|
| RG (t ha-1) | R | 1.9 | 0.15 | 0.29** | 0.29 | 0.40 | 28.1 |
| S | 0.7 | 0.04 | 0.02* | 0.09 | 0.56 | 41.8 | |
| DA (d) | R | - | - | - | - | - | - |
| S | 89.1 | 8.96 | 1.57** | 4.76 | 0.82 | 2.45 | |
| IF (d) | R | 0.28 | 3.6 × 10-13 | 3.85** | 13.05 | 7.0 × 10-14 | 1292 |
| S | -1.1 | 2.91 | 0.58* | 4.37 | 0.68 | 199.9 | |
| AP (cm) | R | 111.3 | 108.63 | 1.90 | 146.35 | 0.74 | 10.9 |
| S | 100.5 | 69.36 | 0.00 | 124.48 | 0.69 | 11.1 | |
| SEN1 (1-10) | R | 1.4 | 0.02 | 0.03** | 0.21 | 0.24 | 31.9 |
| S | 2.4 | 0.03 | 0.06* | 0.33 | 0.23 | 23.7 | |
| SEN2 (1-10) | R | 2.8 | 0.04 | 1.6 × 10-3 | 0.31 | 0.34 | 19.9 |
| S | 4.8 | 0.06 | 0.05** | 0.25 | 0.42 | 10.4 | |
| SEN3 (1-10) | R | 5.5 | 0.31 | 0.26** | 0.66 | 0.51 | 15 |
| S | 7.8 | 0.35 | 0.20* | 1.27 | 0.46 | 14.5 | |
| AUCSEN | R | 91.5 | 58.18 | 37.36** | 173.35 | 0.48 | 14.4 |
| S | 142.1 | 85.82 | 43.01* | 239.70 | 0.51 | 10.9 | |
| NDVI1 | R | 0.2 | 2.1 × 10-4 | 2.0 × 10-4** | 5.1 × 10-4 | 0.49 | 9.64 |
| S | - | - | - | - | - | ||
| NDVI2 | R | 0.4 | 8.7 × 10-4 | 4.0 × 10-4* | 1.9 × 10-3 | 0.56 | 10.72 |
| S | - | - | - | - | - | ||
| NDVI3 | R | 0.4 | 1.4 × 10-3 | 3.7 × 10-4* | 3.2 × 10-3 | 0.60 | 13.69 |
| S | - | - | - | - | - | ||
| NDVI4 | R | 0.6 | 1.4 × 10-3 | 2.6 × 10-4* | 2.7 × 10-3 | 0.63 | 9.09 |
| S | - | - | - | - | - | ||
| AUCNDVI | R | 11.3 | 0.72 | 0.21* | 1.25 | 0.63 | 9.85 |
| S | - | - | - | - | - | ||
| COMBO | R | 15.2 | 7.24 | 4.91** | 9.68 | 0.60 | 20.43 |
| S | - | - | - | - | - |
RG: rendimiento de grano, DA: días a antesis, IF: intervalo de floración, AP: altura de planta, SEN: senescencia, NDVI: índice de vegetación de diferencia normalizada, AUC (AUCSEN y AUCNDVI): área bajo la curva, σ2 g: varianza genética, σ2 g x e: varianza genotipo × ambiente, σ2 e: varianza residual, H2: heredabilidad en sentido amplio, CV: coeficiente de variación. -; información no disponible, líneas que no prosperaron en una localidad se eliminaron del análisis combinado. *: significancia P ≤ 0.05, **: significancia P ≤ 0.01.
Cuadro 2 Estimación de medias, componentes de la varianza y heredabilidad del rendimiento y caracteres secundarios de cruzas de prueba bajo riego (R) y sequía (S).
| Característica | Condición | Media | σ2 g | σ2 g x e | σ2 e | H 2 | CV (%) |
|---|---|---|---|---|---|---|---|
| RG (t ha-1) | R | 4.11 | 0.32 | 0 | 0.82 | 0.61 | 22.08 |
| S | 3.11 | 0.05 | 0.15** | 0.68 | 0.28 | 26.5 | |
| DA (d) | R | 84.8 | 1.45 | 0.37 | 3.82 | 0.56 | 2.3 |
| S | 81.0 | 2.08 | 0.29 | 2.93 | 0.83 | 2.1 | |
| IF (d) | R | 1.1 | 0.30 | 0.09 | 3.15 | 0.26 | 158.5 |
| S | 1.8 | 0.13 | 0.48** | 2.79 | 0.21 | 92.8 | |
| AP (cm) | R | 181.6 | 16.33 | 6.54 | 99.64 | 0.37 | 5.5 |
| S | 211.4 | 101.49 | 15.47* | 84 | 0.84 | 4.3 | |
| SEN1 (1-10) | R | 1.7 | 0.01 | 0.02 | 0.16 | 0.1 | 23.6 |
| S | 2.7 | 0.01 | 0.01 | 0.23 | 0.19 | 17.9 | |
| SEN2 (1-10) | R | - | - | - | - | - | - |
| S | 3.9 | 0.00 | 0.11** | 0.34 | 0.00 | 14.7 | |
| SEN3 (1-10) | R | 3.2 | 0.05 | 0.001 | 0.37 | 0.37 | 19.2 |
| S | 7.3 | 0.07 | 0.14** | 0.50 | 0.40 | 9.5 | |
| AUCSEN | R | - | - | - | - | - | - |
| S | 124.2 | 10.82 | 31.7** | 155.99 | 0.28 | 10.1 | |
| NDVI1 | R | 0.31 | 4.3 × 10-4 | 1.1 × 10-3** | 2.0 × 10-3 | 0.29 | 14.8 |
| S | 0.24 | 1.8 × 10-4 | 8.0 × 10-4** | 7.2 × 10-4 | 0.39 | 11.6 | |
| NDVI2 | R | 0.46 | 2.0 × 10-3 | 9.0 × 10-6 | 3.6 × 10-3 | 0.69 | 12.9 |
| S | 0.41 | 1.8 × 10-3 | 7.7 × 10-4** | 2.7 × 10-3 | 0.71 | 13.1 | |
| NDVI3 | R | 0.53 | 2.8 × 10-3 | 1.4 × 10-18 | 0.01 | 0.69 | 13.6 |
| S | 0.49 | 2.7 × 10-3 | 2.1 × 10-3** | 4.1 × 10-3 | 0.72 | 13.4 | |
| NDVI4 | R | - | - | - | - | - | - |
| S | 0.59 | 2.8 × 10-3 | 2.0 × 10-3** | 3.9 × 10-3 | 0.74 | 10.8 | |
| AUCNDVI | R | 12.15 | 1.47 | 0.00 | 1.95 | 0.75 | 11.5 |
| S | 12.00 | 1.22 | 1.03** | 1.72 | 0.72 | 11.2 | |
| COMBO | R | - | - | - | - | - | |
| S | 19.71 | 2.90 | 1.77** | 8.96 | 0.65 | 15.4 |
RG: rendimiento de grano, DA: días a antesis, IF: intervalo de floración, AP: altura de planta, SEN: senescencia, NDVI: índice de vegetación de diferencia normalizada, AUC (AUCSEN y AUCNDVI): área bajo la curva, σ2 g: varianza genética, σ2 g x e: varianza genotipo × ambiente, σ2 e: varianza residual, H 2 : heredabilidad en sentido amplio, CV: coeficiente de variación, *: significancia P ≤ 0.05, **: significancia P ≤ 0.01.
En las CP la antesis fue más temprana que en las líneas DH. Las heredabilidades fueron moderadas y altas (H 2 > 0.56) y los coeficientes de variación bajos (< 2.5 %). La antesis se presentó antes en las cruzas de prueba bajo sequía que en las de riego (Figura 1), con lo que se presenta un escape a la sequía de aproximadamente 4 d.
Figura 1 Diagramas de caja donde se muestran el primer cuartil, mediana, tercer cuartil y valores extremos de las mediciones de líneas DH y cruzas de prueba bajo sequía y riego. LS: líneas DH bajo sequía, LR: líneas DH bajo riego, CPS: cruzas de prueba bajo sequía, CPR: cruzas de prueba bajo riego. A) rendimiento de grano, B) días a antesis, C) intervalo de floración, D) altura de planta, E) senescencia SEN1, F) SEN2, G) SEN3, H) área bajo la curva AUCSEN, I) NDVI1, J) NDVI2, K) NDVI3, L) NDVI4, M) AUCNDVI, N) COMBO.
El intervalo de floración presentó diferencias entre genotipos únicamente en líneas DH bajo sequía, con un intervalo de -8 a 2.3 d mientras que en las CP en sequía el intervalo fue reducido, de 1.3 a 2.3 d. Las líneas bajo sequía presentaron una heredabilidad intermedia (H 2 = 0.68), mientras que las cruzas de prueba presentaron baja heredabilidad (H 2 = 0.21).
La altura de planta presentó diferencias significativas entre las poblaciones en las dos condiciones, donde las poblaciones bajo sequía presentaron valores mayores que las de riego para las CP. La diferencia en altura promedio de la población de líneas DH y CP bajo riego (H 2 = 0.69 y 0.37) y sequía (H 2 = 0.74 y 0.84) fue de 10.8 y 29.8 cm, respectivamente.
La calificación de SEN1 no presentó diferencias entre genotipos para la población CP, y las heredabilidades fueron bajas (H 2 < 0.24). Mientras que en las calificaciones de SEN2 sólo en líneas DH bajo riego y sequía se presentaron diferencias con medias de calificación de 2.8 para riego (H 2 = 0.34) y sequía de 4.8 (H 2 = 0.42). En SEN3 se presentaron diferencias entre genotipos en las diferentes poblaciones de riego y sequía con heredabilidades intermedias, en líneas DH la media bajo riego fue de 5.5 (H 2 = 0.51) y para sequia de 7.8 (H 2 = 0.46), en tanto que en cruzas de prueba en riego fue de 3.2 (H 2 = 0.37) y en sequía 7.3 (H 2 = 0.40). Los resultados en AUCSEN revelaron diferencias entre genotipos en las líneas DH bajo riego (H 2 = 0.48) y sequía (H 2 = 0.51), mientras que las CP bajo sequía no presentaron diferencias entre genotipos (H 2 = 0.28).
Los valores de NDVI en las CP bajo riego fueron mayores, debido quizá a una mejor fertilidad del suelo. Las heredabilidades van de forma creciente de NDVI1 a NDVI4 (H 2 = 0.29 a 0.74). AUCNDVI presenta heredabilidad alta (H 2 = 0.63 a 0.75).
Los resultados de COMBO sólo se presentan para líneas DH bajo riego (H 2 = 0.60) y cruzas de prueba en sequía (H 2 = 0.65), donde se obtuvieron diferencias entre genotipos.
y condiciones de riego
En el Cuadro 3 se presentan las correlaciones de las características secundarias con el rendimiento. En las líneas DH bajo riego la mayoría de las características tuvieron asociaciones altas con el rendimiento, con excepción del intervalo de floración y senescencia. Para las CP y DH en riego todas las calificaciones de senescencia, AUCSEN, NDVI (todas las capturas), AUCNDVI y COMBO presentaron correlaciones altas y positivas; IF presentó correlación intermedia, mientras que el resto no presentó correlación. En las Líneas DH en sequía la altura de planta presentó la mayor asociación con el rendimiento con una correlación genotípica de 0.54. El resto de las características tuvieron correlaciones bajas y negativas. Para las cruzas de prueba en riego, los caracteres IF, SEN1, NDVI (todas las tomas), AUCNDVI y COMBO presentaron correlaciones genéticas altas y positivas.
Cuadro 3 Correlación genética (r g ) y fenotípica (r f ) del rendimiento de grano con caracteres secundarios de líneas DH y cruzas de prueba (CP) bajo riego y sequía.
| Carácter | Líneas DH Riego | CP Riego | Líneas DH Sequía | CP Sequía | |||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| rg | r f | r g | r f | r g | r f | r g | r f | ||||
| DA | 0.78** | -0.01 | -0.01 | -0.07 | -0.19* | -0.17* | -0.17* | -0.09 | |||
| IF | 0 | -0.24** | 0.55** | 0.15 | -0.26** | -0.17* | 0.73** | 0.07 | |||
| AP | 0.30** | 0.26** | -0.09 | 0.04 | 0.54** | 0.41** | 0.42** | 0.34** | |||
| SEN1 | -0.45** | -0.26** | 0.99** | 0.42** | -0.10 | -0.20* | 0.99** | 0.19** | |||
| SEN2 | -0.12 | -0.22** | 0.99** | 0.42** | -0.11 | -0.27** | 0.00 | -0.06 | |||
| SEN3 | -0.33** | -0.33** | 0.59** | 0.27** | -0.12 | -0.23** | -0.26** | -0.12 | |||
| AUCSEN | -0.29** | -0.33** | 0.99** | 0.46** | -0.14 | -0.29** | -0.58** | -0.06 | |||
| NDVI1 | 0.77** | 0.59** | 0.99** | 0.67** | - | - | 0.99** | 0.30** | |||
| NDVI2 | 0.85** | 0.61** | 0.85** | 0.70** | - | - | 0.99** | 0.44** | |||
| NDVI3 | 0.89** | 0.63** | 0.84** | 0.72** | - | - | 0.99** | 0.43** | |||
| NDVI4 | 0.99** | 0.66** | - | - | - | - | 0.90** | 0.42** | |||
| AUCNDVI | 0.85** | 0.65** | 0.86** | 0.75** | - | - | 0.94** | 0.44** | |||
| COMBO | 0.68** | 0.60** | - | - | - | - | 0.96** | 0.55** | |||
DA: días a antesis, IF: intervalo de floración, AP: altura de planta, SEN: senescencia, NDVI: índice de vegetación de diferencia normalizada, AUC (AUCSEN y AUCNDVI): área bajo la curva, *: significancia P ≤ 0.05, **: significancia P ≤ 0.01.
con tolerancia a sequía
QTLs para rendimiento
En la población de líneas DH con riego se ubicó un QTL para rendimiento en el Bin 4.1 que explica el 17.8 % de la variación fenotípica, el efecto de sustitución de esta región contribuye con un incremento de 0.12 t ha-1 aportado por LPS-C7-F64 (Cuadro 4). En la población de líneas DH bajo sequía se localizaron dos QTLs en los Bin 1.07 y 4.09 que participan con el 7.9 y 9.8 % de la variación fenotípica con 0.2 y 0.6 t ha-1, el primero proveniente de la línea CML495 y el segundo de LPS-C7-F64.
Cuadro 4 QTLs obtenidos a través mapeo del intervalo compuesto (ICIM) de los rasgos de las líneas DH y cruzas de prueba bajo riego y sequía.
| Carácter | Cromosoma | Bin† | Posición (cM)†† | Intervalo del marcador (cM) | Intervalo físico (Mb)¶ | LOD¶¶ | PVFx | Addxx |
|---|---|---|---|---|---|---|---|---|
| SEN3_LR | 1 | 1.01/1.02 | 131 | 128.25 - 131.25 | 14.26 - 12.24 | 3.11 | 8.66 | -0.11 |
| Combo_CPS | 1 | 1.02/1.03 | 118 | 116.51 - 118.28 | 32.58 - 28.66 | 2.57 | 7.49 | 0.28 |
| DA_LS | 1 | 1.03 | 113 | 112.42 - 113.02 | 38.49 - 36.77 | 3.94 | 12.91 | 0.88 |
| NDVI2_CPS | 1 | 1.03 | 107 | 104.12 - 107.04 | 46.37 - 42.29 | 3.48 | 9.99 | 0.01 |
| AP_LR | 1 | 1.03 | 104 | 103.53 - 104.12 | 49.77 - 46.37 | 2.83 | 6.28 | 2.28 |
| NDVI1_CPS | 1 | 1.03 | 104 | 103.53 - 104.12 | 49.77 - 46.37 | 2.72 | 7.2 | 0.002 |
| NDVI3_CPS | 1 | 1.03 | 104 | 103.53 - 104.12 | 49.77 - 46.37 | 3.38 | 9.76 | 0.01 |
| NDVI4_CPS | 1 | 1.03 | 104 | 103.53 - 104.12 | 49.77 - 46.37 | 3.37 | 8.89 | 0.01 |
| AUCNDVI_CPS | 1 | 1.03 | 104 | 103.53 - 104.12 | 49.77 - 46.37 | 3.82 | 10.94 | 0.26 |
| SEN1_LR | 1 | 1.06 | 75 | 73.88 - 75.65 | 189.71 - 182.01 | 2.83 | 8.78 | 0.02 |
| RG_CPS | 1 | 1.07 | 55 | 53.86 - 55.06 | 227.17 - 221.20 | 2.72 | 7.92 | 0.02 |
| AP_LR | 2 | 2.07 | 36 | 35.35 - 36.52 | 204.30 - 201.75 | 4.06 | 9.2 | -2.76 |
| SEN3_CPS | 3 | 3.04 | 88 | 86.9 - 88.06 | 95.83 - 49.74 | 4.94 | 13.26 | 0.05 |
| NDVI4_CPS | 3 | 3.04 | 87 | 86.9 - 88.06 | 95.83 - 49.74 | 2.53 | 6.62 | 0.01 |
| SEN3_CPS | 3 | 3.07 | 52 | 51.95 - 53.11 | 196.10 - 191.58 | 2.96 | 7.71 | -0.04 |
| AP_CPS | 3 | 3.09 | 21 | 19.93 - 22.3 | 223.35 - 221.21 | 2.83 | 7.63 | 1.56 |
| RG_LR | 4 | 4.10 | 125 | 123.32 - 126.89 | 239.21 - 239.97 | 4.96 | 17.84 | -0.12 |
| AUCSEN_LR | 4 | 4.05 | 43 | 42.97 - 43.54 | 61.37 - 134.00 | 3.2 | 9.86 | 1.62 |
| SEN1_LS | 4 | 4.05 | 43 | 42.97 - 43.54 | 61.37 - 134.00 | 3.4 | 10.01 | 0.03 |
| SEN3_LS | 4 | 4.05 | 43 | 42.97 - 43.54 | 61.37 - 134.00 | 3.87 | 12 | 0.13 |
| AUCSEN_LS | 4 | 4.05 | 43 | 42.97 - 43.54 | 61.37 - 134.00 | 4.47 | 13.81 | 2.27 |
| Combo_LR | 4 | 4.05 | 45 | 43.54 - 46.44 | 110.49 - 148.67 | 5 | 12.99 | -0.72 |
| AP_LR | 4 | 4.05/4.06 | 47 | 46.44 - 47.02 | 148.67 - 155.30 | 5.13 | 11.64 | -3.1 |
| SEN1_LR | 4 | 4.06 | 49 | 47.02 - 49.38 | 155.30 - 156.72 | 3.23 | 9.98 | 0.02 |
| AP_LS | 4 | 4.06 | 74 | 66.57 - 4.25 | 170.14 - 167.07 | 4.19 | 12.57 | -2.52 |
| IF_CPS | 4 | 4.08 | 95 | 92.88 - 95.39 | 180.43 - 181.36 | 2.53 | 6.78 | 0.04 |
| RG_LS | 4 | 4.09 | 108 | 107.29 - 108.46 | 209.56 - 222.18 | 2.95 | 9.79 | -0.06 |
| RG_LS | 4 | 4.09 | 115 | 112.11 - 115.04 | 233.93 - 235.20 | 2.62 | 8.75 | -0.04 |
| SEN1_CPS | 5 | 5.01 | 84 | 84 - 87.53 | 5.29 - 4.54 | 2.51 | 7.78 | -0.01 |
| AUCSEN_LR | 5 | 5.03 | 62 | 61.63 - 62.88 | 69.09 - 62.76 | 2.68 | 8.07 | 1.47 |
| NDVI1_CPS | 5 | 5.03 | 60 | 59.89 - 61.63 | 77.67 - 69.09 | 3.43 | 9.21 | 0.002 |
| IF_LS | 5 | 5.03/5.04 | 59 | 58.14 - 59.89 | 89.17 - 77.67 | 3.32 | 10.91 | 0.48 |
| SEN3_LR | 5 | 5.04 | 50 | 49.92 - 51.69 | 168.23 - 165.61 | 2.98 | 8.22 | 0.11 |
| Combo_LR | 5 | 5.05 | 42 | 41.57 - 43.98 | 174.86 - 172.70 | 4.56 | 11.61 | -0.68 |
| AP_LS | 6 | 6.01 | 61 | 60.72 - 63.61 | 60.18 - 30.46 | 3.76 | 11.12 | 2.35 |
| SEN3_LR | 6 | 6.04 | 44 | 43.44 - 44.03 | 110.16 - 106.67 | 2.5 | 6.85 | 0.09 |
| IF_LS | 6 | 6.07/6.08 | 1 | 0 - 1.18 | 167.57 - 166.33 | 3.24 | 9.77 | 0.44 |
| AP_LR | 7 | 7.03 | 49 | 48.07 - 49.26 | 136.38 - 137.35 | 5.87 | 13.59 | 3.34 |
| Combo_LR | 7 | 7.03 | 49 | 48.07 - 49.26 | 136.38 - 137.35 | 3.94 | 9.89 | 0.62 |
| AP_CPS | 8 | 8.04/8.05 | 73 | 72.61 - 74.46 | 116.33 - 123.33 | 2.73 | 7.35 | 1.53 |
| SEN1_LS | 9 | 9.02/9.03 | 30 | 29 - 30.75 | 23.11 - 32.95 | 4.33 | 13.16 | -0.03 |
| SEN3_LS | 9 | 9.02/9.03 | 29 | 29 - 30.75 | 23.11 - 32.95 | 2.60 | 7.85 | -0.10 |
| AUCSEN_CPS | 9 | 9.02/9.03 | 29 | 29 - 30.75 | 23.11 - 32.95 | 2.70 | 7.83 | -0.35 |
| IF_CPS | 9 | 9.06 | 62 | 61.79 - 62.98 | 146.48 - 147.24 | 2.95 | 8.04 | -0.04 |
LR: líneas DH bajo riego, LS: líneas bajo sequía, CPR: cruzas de prueba bajo riego; CPS: cruzas de prueba bajo sequía, RG: rendimiento de grano, DA: días a antesis, IF: intervalo de floración, AP: altura de planta, SEN: senescencia, AUC (AUCSEN y AUCNDVI) área bajo la curva, NDVI: índice de vegetación de diferencia normalizada. †Bin: región referida en el cromosoma (https://www.maizegdb.org/bin_viewer), ††centiMorgans; ¶Millones de pares de bases, ¶¶valor LOD, xPorción de la varianza fenotípica del QTL, xxEfecto aditivo de cada QTL.
QTL para días a antesis
Se ubicó un QTL para días a antesis en las líneas DH bajo sequía, localizado en el Bin 1.03 (LOD 3.94) que explica el 12.9 % de la variación fenotípica y que puede modificar los días a antesis en hasta un día (aditividad +0.88).
QTLs para intervalo de floración
En el intervalo de floración la línea CML495 contribuyó con dos QTLs favorables ubicados en los Bin 5.02/5.03 y 6.07/6.08 que explicaron 10.9 y 9.8 % de la variación fenotípica, con efectos de asincronía de 0.48 y 0.44 d bajo condiciones de sequía.
QTLs para senescencia
En el Bin 4.5 se ubicó un QTL constitutivo de senescencia que explica el 13.8 % de la variación fenotípica en AUCSEN_LS y 10 % en AUCSEN_LR, donde el alelo favorable proviene de la línea LPSC7F64. Otros QTLs importantes son los localizados en el Bin 9.02/9.03 provenientes de las líneas DH y cruzas de prueba (SEN1_LS, SEN3_LS y AUCSEN_CPS) con una varianza fenotípica entre 7.8 y 13.2 %, donde los alelos favorables provienen de la línea CML495.
QTLs para vigor
Se localizaron siete QTLs para NDVI y ocho para altura de planta. Cinco QTLs de altura de planta y todos los de NDVI provienen de la línea CML495, lo que confiere un desarrollo vigoroso en la etapa vegetativa, que está altamente correlacionado con un incremento del rendimiento. En el Bin 1.03 se ubicaron los QTL de NDVI y altura de planta que explican hasta el 10.9 % (AUCNDVI_CPS) y 6.3 % (AP_LR) de la variación fenotípica y confiere un incremento en altura de planta de 2.28 cm. Tal alelo podría implicar un QTL constitutivo proveniente de la línea CML495. La región 5.03 de NDVI1_CPS explica el 9.2 % de la variación. Estas dos importantes regiones también fueron localizadas en el trabajo de Trachsel et al. (2016).
QTLs para altura de planta
Para condiciones de sequía se ubicaron en los Bin 3.09 4.06, 6.01 y 8.04/8.05, donde la variación fenotípica explicada fue de 7.4 y 12.6 % y valores aditivos de 1.5 y 2.5 cm. Otro QTL importante aportado por la línea CML495 es de AP_LR en el Bin 7.03 que explica el 13.6 % de la variación fenotípica y un incremento de 3.34 cm. Esta región también fue localizada por Almeida et al. (2013) para altura de planta.
QTLs para COMBO
Los QTLs que integran la relación vigor/senescencia para riego se encontraron junto con las de altura de planta en los Bin 4.05/4.06 y 7.03 y de manera individual en el 5.05, que juntos explicaron el 34.5 % de la variación fenotípica. Para sequía se ubicó un QTL en el Bin 1.02/1.03, que proporciona el 7.5 % de la variación fenotípica.
La heredabilidad del rendimiento se considera moderada para las líneas DH bajo condiciones de riego y sequía, en tanto que en las CP bajo riego la heredabilidad resultante es considerada como alta, y para sequía se comportó como baja, lo que pudo deberse a las altas varianzas residual y de interacción G × A, como lo reportan Lafitte et al. (2003). Las heredabilidades de los caracteres secundarios en las cruzas de prueba en su mayoría fueron altas bajo sequía y altamente correlacionadas con el rendimiento. En esta investigación se pusieron de manifiesto correlaciones fenotípicas y genotípicas importantes entre el rendimiento, el intervalo de floración (líneas DH), altura de planta, SEN3, AUCSEN (CP), NDVI2, NDVI3, NDVI4, AUCNDVI y COMBO, los que también se muestran en otros trabajos (Almeida et al., 2013; Lu et al., 2011; Messmer et al., 2011; Trachsel et al., 2016).
En el presente estudio los caracteres secundarios NDVI y altura de planta fueron los que mejor correlación tuvieron con el rendimiento y con alta heredabilidad, tanto en estrés hídrico como en condiciones óptimas de humedad; la misma tendencia se presentó en la relación vigor y senescencia, referida como COMBO. El contenido de clorofila y la senescencia retrasada están correlacionados con NDVI y altura de planta, y éstos a su vez con el rendimiento, como lo describe Lu et al. (2011), por lo que, utilizar a COMBO como carácter secundario ofrece una mejor identificación de genotipos tolerantes a sequía. De los QTLs de rendimiento y caracteres secundarios de rendimiento con alto valor LOD y porcentaje de variación fenotípica se ubicaron varios segmentos contenidos en los Bin 1.02/1.03, 1.07, 4.05, 4.10 y 7.03.
QTL rendimiento bajo sequía Bin 1.07
Este QTL para rendimiento se ha identificado en otros estudios bajo condiciones normales (Frascaroli et al., 2007) y en condiciones de estrés (Zaidi et al., 2015). En la región del cromosoma 1 se encuentra el gen agpsl1, que está relacionado con el incremento de almidón en el grano (Huang et al., 2014; Liu et al., 2016).
QTL rendimiento bajo riego Bin 4.10
Este QTL es un candidato para ser introducido en líneas elite de maíz. Almeida et al. (2013) lo ubicaron en la línea CML504 en la posición 239.24 - 242.02 Mb. En esta región se ubica el gen NRT2.1. El estudio de Sorgonà et al. (2011) sugiere que NRT2.1 es un componente importante del sistema transportador de alta afinidad (iHATS) para la absorción de NO3 - a lo largo del eje de la raíz del maíz.
QTL Bin 1.03 altura de planta, NDVI (vigor) y COMBO_CPS
Esta región está relacionada con el rendimiento en condiciones de sequía (Millet et al., 2016; Zhu et al., 2011) y en condiciones normales (Semagn et al., 2013) y también se han ubicado QTLs para altura de planta (Frascaroli et al., 2007), stay Green (Almeida et al., 2014) y senescencia (Trachsel et al., 2016). Los efectos de esta región indican una favorable capacidad de crecimiento y producción bajo condiciones de sequía.
QTL 4.05 COMBO y senescencia
En el Bin 4.05 se han ubicado varios QTLs asociados a rendimiento bajo sequia (Zhu et al., 2011), rendimiento en condiciones normales (Chen et al., 2016), altura de planta (Lu et al., 2012) y senescencia (Trachsel et al., 2016). En esta región se ubica el gen nnr1 que está relacionado con la absorción de nitrógeno en maíz (Liseron-Monfils et al., 2013).
QTL 7.03 COMBO
Se ha localizado un QTL para altura de planta tanto en riego como en sequía (Almeida et al., 2014; Peiffer et al., 2014), rendimiento en sequía (Almeida et al., 2013) en el que se encuentra el gen glutation transferasa 23 (gst23) (Chen et al., 2012; Varga et al., 2012) y fitoeno sintasa 3 (psy3), donde se demostró que la acumulación de ABA en raíz está limitada por la expresión de PSY3 en condiciones de estrés por sequía (Li et al., 2009).
Conclusiones
Los caracteres secundarios NDVI, altura de planta, senescencia (Sen3 y AUCSEN) y COMBO proporcionaron mejor información y permitieron identificar genotipos con buen rendimiento de líneas DH y cruzas de prueba en etapas tempranas bajo condiciones de sequía y régimen de humedad adecuado. Los QTLs ubicados en los cromosomas 1, 4, y 7 son los de mayor relevancia y están relacionados con vigor y rendimiento; éstos pueden ser útiles en la introgresión asistida por marcadores en líneas e híbridos de trópico para aumentar el rendimiento bajo condiciones de sequía y de buena humedad. Los resultados de rendimiento y los caracteres secundarios son menos sesgados al utilizar las mediciones de líneas y cruzas de prueba, lo que reduce la interacción con alelos deletéreos.
