Introducción
⌅El caucho natural (Hevea brasiliensis [Willd. ex A. Juss.] Müll. Arg.), constituye una materia prima
estratégica a nivel global, esencial para la manufactura de más de
50,000 productos industriales (Priyadarshan, 2011Priyadarshan, P.M. 2011. Biology of Hevea Rubber. Wallingford, UK: CABI, 350 p.
).
La demanda actual exige un incremento sustancial en la productividad
por unidad de área; no obstante, el mejoramiento genético de esta
especie enfrenta desafíos críticos derivados de su biología perenne,
dado que los ciclos de selección convencionales pueden extenderse entre
25 y 30 años antes de la liberación comercial de un clon (Priyadarshan, 2017Priyadarshan, P.M. 2017. “Refinements to Hevea rubber breeding”. Tree Genetics & Genomes, 13(1): 1-17, DOI: http://doi.org/10.1007/s11295-017-1101-8.
).
En este contexto, la implementación de pruebas tempranas surge como una
herramienta indispensable para identificar genotipos sobresalientes en
etapas juveniles, permitiendo reducir significativamente los periodos de
evaluación (Kalpani et al., 2020Kalpani, K.L.M., Withange, S.P., Palihakkara, I.R. 2020. “Selection of superior genotypes at early stage of the rubber (Hevea brasiliensis) breeding cycle”. International Journal for Research in Applied Sciences and Biotechnology, 7(4): 146-154, DOI: http://doi.org/10.31033/ijrasb.7.4.17.
; Chandrasekhar et al., 2007Chandrasekhar,
T.R., Marattukalam, J.G., Mercykutty, V.C., Priyadarshan, P.M. 2007.
“Age of yield stabilization and its implications for optimising
selection and shortening breeding cycle in rubber (Hevea brasiliensis Muell. Arg.)”. Euphytica, 156(1-2): 67-75, DOI: http://doi.org/10.1007/s10681-006-9352-8.
).
Esta
estrategia se fundamenta en la existencia de correlaciones genéticas
positivas entre los caracteres evaluados en fases inmaduras y el
rendimiento en la etapa adulta, premisa ampliamente documentada en
diversos programas de mejoramiento (Gonçalves et al., 2005aGonçalves,
P. de S., Bortoletto, N., Cardinal, Á.B.B., Gouvêa, L.R.L., Costa,
R.B., Moraes, M.L.T. 2005. “Age-age correlation for early selection of
rubber tree genotypes in São Paulo State, Brazil”. Genetics and Molecular Biology, 28(4): 758-764, DOI: http://doi.org/10.1590/S1415-47572005000500018.
, 2005bGonçalves,
P. de S., Moraes, M.L.T., Bortoletto, N., Costa, R.B., Gonçalves,
E.C.P. 2005. “Genetic variation in growth traits and yield of rubber
trees (Hevea brasiliensis) growing in the Brazilian state of São Paulo”. Genetics and Molecular Biology, 28(4): 765-772, DOI: http://doi.org/10.1590/S1415-47572005000500019.
).
Tradicionalmente, la selección se ha basado en pruebas de rendimiento a
corto plazo; sin embargo, este enfoque debe integrarse con el análisis
dinámico del desarrollo ontogénico. Estudios contemporáneos confirman
que caracteres como el incremento del perímetro del tallo y el espesor
de corteza poseen una heredabilidad moderada a alta, asegurando que la
selección precoz genere un progreso genético predecible y eficiente (Peres Silva et al., 2014Peres
Silva, G.A., Gouvêa, L.R.L., Verardi, C.K., Bombonato de Oliveira,
A.L., Gonçalves, P.S. 2014. “Annual growth increment and stability of
rubber yield in the tapping phase in rubber tree clones: Implications
for early selection”. Industrial Crops and Products, 52: 801-808, DOI: http://doi.org/10.1016/j.indcrop.2013.12.010.
; Suwanchana et al., 2020Suwanchana,
S., Phumichai, C., Jenweerawat, S., Koonawootrittriron, S., Phumichai,
T., Teerawatnasuk, K. 2020. “Estimation of genetic parameters for rubber
yield and girth growth in a synthetic population”. Agriculture and Natural Resources, 54(6): 603-608, DOI: http://doi.org/10.34044/j.anres.2020.54.6.05.
; Gouvêa et al., 2013bGouvêa,
L.R.L., Silva, G.A.P., Verardi, C.K., Oliveira, A.L.B., Gonçalves, P.S.
2013b. “Simultaneous selection of rubber yield and girth growth in
young rubber trees”. Industrial Crops and Products, 50: 39-43, DOI: http://doi.org/10.1016/j.indcrop.2013.06.040.
).
Asimismo, la eficiencia productiva depende de la maduración temprana
del sistema laticífero, factor que puede monitorearse mediante
parámetros bioquímicos y fisiológicos desde la fase de vivero (Madigasekara et al., 2025Madigasekara,
M.M.S.K., Withanage, S.P., Liyanage, K.K., Dharmarathne, T.T.D.,
Meegahakumbura, M.K. 2025. “Precise selection of superior Hevea genotypes from 2014 hand-pollinated progeny”. Tropical Agricultural Research & Extension, 28(2), DOI: http://doi.org/10.4038/tare.v28i2.5779.
).
Complementariamente,
la eficacia de la selección precoz no depende exclusivamente del
rendimiento absoluto, sino de la homeostasis fisiológica frente a la
interacción genotipo-ambiente (G × A). La estabilidad de los caracteres
productivos es un componente crítico, tal como lo han señalado diversos
estudios sobre la estabilidad clonal de la producción de látex y la
susceptibilidad al secamiento (Omokhafe y Alika, 2003Omokhafe, K.O., Alika, J.E. 2003. “Clonal stability of latex yield in eleven clones of Hevea brasiliensis (Muell. Arg.)”. Genetics and Molecular Biology, 26(3): 313-317, DOI: http://doi.org/10.1590/S1415-47572003000300016.
; Omokhafe, 2004Omokhafe, K.O. 2004. “Clonal stability of tree dryness in Hevea brasiliensis Muell. Arg.”. Genetics and Molecular Biology, 27(2): 242-244, DOI: http://doi.org/10.1590/S1415-47572004000200019.
).
En regiones con marcada estacionalidad, los materiales genéticos pueden
exhibir una eficiencia sobresaliente bajo condiciones óptimas, pero
carecer de resiliencia ante eventos de estrés hídrico o térmico (Ajith et al., 2022Ajith, S.V., Nair, N.N., Sathik, M.B.M., Meenakumari, T. 2022. “Breeding for low-temperature stress tolerance in Hevea brasiliensis: screening of newly developed clones using latex biochemical parameters”. Journal of Rubber Research, 24(1): 63-76, ISSN: 2510-1498, DOI: http://doi.org/10.1007/s42464-021-00143-2.
; Jacob et al., 2021Jacob, J., Gitz, V., Gohet, E., et al., 2022. “Natural rubber contributions to adaptation to climate change”. In: XV World Forestry Congress, 2-6 de mayo, Seúl, Corea del Sur. (Nota: En congresos suele aceptarse et al. si el comité es muy extenso).
).
Por ello, el empleo de modelos lineales mixtos (REML/BLUP), es
fundamental para desglosar la varianza genética de la ambiental,
permitiendo una predicción más precisa de los valores genéticos (Gouvêa et al., 2013Gouvêa,
L.R.L., Silva, G.A.P., Verardi, C.K., Oliveira, A.L.B., Gonçalves, P.S.
2013b. “Simultaneous selection of rubber yield and girth growth in
young rubber trees”. Industrial Crops and Products, 50: 39-43, DOI: http://doi.org/10.1016/j.indcrop.2013.06.040.
).
La validación local de series clonales representa, por tanto, el
vínculo necesario entre la investigación básica y la recomendación
técnica fundamentada para la heveicultura moderna.
Bajo estas premisas, el presente estudio tuvo como objetivo evaluar la eficacia de las pruebas tempranas como criterio de selección mediante el seguimiento multianual de genotipos de las series IAN e IRCA. Se buscó determinar el vigor vegetativo y la estabilidad productiva de dichos clones en respuesta a la variabilidad climática estacional, a fin de identificar materiales con un balance óptimo entre potencial de rendimiento y resiliencia ambiental que sustente el uso de las pruebas tempranas como herramientas de predicción robusta en el trópico húmedo mexicano.
Materiales y Métodos
⌅Área de estudio y caracterización ambiental
⌅La investigación se realizó en el Campo Experimental El Palmar del INIFAP, ubicado en Tezonapa, Veracruz, México (18°32’ N, 96°46’ O; 180 msnm). La región presenta un clima cálido húmedo (Am), con una precipitación media anual de 2,800 mm y una temperatura media de 25 °C. Los suelos predominantes son Acrisoles (FAO/WRB), de naturaleza ácida y textura franco-arcillosa, característicos de las zonas productoras de caucho en el trópico húmedo mexicano.
Para asegurar la precisión del
análisis, se efectuó una caracterización ambiental detallada mediante el
software ClimateNA v7.60 (Wang et al., 2016Wang,
T., Hamann, A., Spittlehouse, D.L., Carroll, C. 2016. “Locally
downscaled and spatially customizable climate data for historical and
future periods for North America”. PLOS ONE, 11(6): e0156720, DOI: http://doi.org/10.1371/journal.pone.0156720.
).
Se extrajeron variables mensuales de temperatura (máxima, mínima y
media), precipitación (PPT), radiación (RAD) e índice térmico (DD5). Se
puso especial énfasis en el déficit hídrico mensual (DEFH), calculado
como balance hídrico acumulado, para cuantificar el impacto del estrés
hídrico estacional sobre la fisiología de los clones, factor ampliamente
documentado como limitante en la expresión del potencial genético (Priyadarshan y Clément-Demange, 2004Priyadarshan, P.M., Clément-Demange, A. 2004. “Breeding Hevea rubber: Formal and molecular genetics”. In: Hall, J. C., Dunlap, J. C., & Friedmann, T. (eds.), Advances in Genetics, vol. 52, Elsevier, pp. 51-115, DOI: http://doi.org/10.1016/S0065-2660(04)52003-8.
).
Material biológico y manejo experimental
⌅Se evaluaron cinco clones de Hevea brasiliensis seleccionados por su potencial genético: IAN 710 e IAN 873 (origen
brasileño), e IRCA 41, IRCA 19 e IRCA 111 (origen africano/francés,
introducidos mediante convenio INIFAP-CIRAD). Las plantas se produjeron
en bolsas de polietileno (20 x 40 cm) e injertadas sobre patrones de
seis meses de edad. El experimento se estableció bajo un diseño de
bloques completos al azar con cuatro repeticiones. Cada unidad
experimental consistió en 20 árboles, con un marco de plantación de 3 x 3
m (densidad de 1,111 árboles ha⁻¹). El manejo agronómico incluyó el
control sistemático de malezas, podas de formación y un programa de
fertilización semestral, ajustado a los requerimientos nutricionales
específicos para el desarrollo de la especie en la región (Ortiz, 2017Ortiz,
E.C. 2017. “Establecimiento y manejo de plantaciones de hule”. In:
Zamora-Martínez, M. C., Sánchez C., R., Zamora M., B. P., & Nieto de
Pascual P., Ma. C. del C. (comps.), Paquete Tecnológico para el Cultivo del Hule, México: Cenid-Comef, INIFAP, pp. 11-45.
).
Variables de respuesta y dinámica ontogénica
⌅El
seguimiento fenológico y productivo se realizó durante 48 meses,
segmentados en etapas anuales (P1 a P4), para facilitar el análisis de
estabilidad temporal, siguiendo los protocolos propuestos para la
selección temprana de genotipos (Gonçalves et al., 2005aGonçalves,
P. de S., Bortoletto, N., Cardinal, Á.B.B., Gouvêa, L.R.L., Costa,
R.B., Moraes, M.L.T. 2005. “Age-age correlation for early selection of
rubber tree genotypes in São Paulo State, Brazil”. Genetics and Molecular Biology, 28(4): 758-764, DOI: http://doi.org/10.1590/S1415-47572005000500018.
). Las variables evaluadas fueron:
- Vigor
vegetativo: Se midió mensualmente el diámetro a la altura del pecho
(DAP), a 1,30 m y la altura total (ALT). Con estos datos, se calcularon
las tasas de incremento anual (I-DAP e I-ALT), como indicadores de
precocidad genotípica, técnica validada en diversos estudios de
selección precoz (Suwanchana et al., 2020Suwanchana,
S., Phumichai, C., Jenweerawat, S., Koonawootrittriron, S., Phumichai,
T., Teerawatnasuk, K. 2020. “Estimation of genetic parameters for rubber
yield and girth growth in a synthetic population”. Agriculture and Natural Resources, 54(6): 603-608, DOI: http://doi.org/10.34044/j.anres.2020.54.6.05.
). - Rendimiento de caucho seco: Se determinó mediante la prueba Hamaker-Morris-Mann (H-M-M), modificada (Lopes et al., 1999Lopes, U.V., Marques, J.R.B., Monteiro, W.R. 1999. “Estimates of genetic parameters in rubber tree (Hevea spp). 1. General combining ability in juvenile seedlings”. Agrotrópica, 11(3): 163-168.
). Se empleó un sistema de pica en semi-espiral descendente (S/2), a una altura de 60 cm del suelo, con una frecuencia de día alterno (d/2). El látex recolectado cada cinco picas (serie), por árbol fue coagulado y deshidratado a 70 °C hasta peso constante para determinar el rendimiento por serie (RS, g) y, mediante el acumulado de las series mensuales, se obtuvo el rendimiento mensual (REN, g).
Análisis estadístico y modelos de estabilidad
⌅El análisis de los datos se realizó en SAS® v9.4 mediante un enfoque de dos fases que garantiza la precisión de los estimadores genéticos:
- Fase
1: Modelos Lineales Mixtos: Para evaluar el vigor y rendimiento, se
utilizó el procedimiento PROC MIXED mediante máxima verosimilitud
restringida (REML). Los clones y los periodos se definieron como efectos
fijos, mientras que los bloques e individuos dentro de bloques se
trataron como efectos aleatorios. La autocorrelación temporal de los
datos longitudinales se modeló con una estructura de simetría compuesta
(DDFM=BETWITHIN), enfoque ampliamente utilizado para maximizar la
eficiencia en la estimación de parámetros genéticos (Bombonato et al., 2015Bombonato,
T., Silva, P.H.S., Pisaroglo de Carvalho, M., Gouvêa, L.R.L., Verardi,
C.K., Gonçalves, P.S. 2015. “Rubber tree ortet-ramet genetic correlation
and early selection efficiency to reduce rubber tree breeding cycle”. Industrial Crops and Products, 77: 855-860, DOI: http://doi.org/10.1016/j.indcrop.2015.09.020.
; Gouvêa et al., 2013aGouvêa, L.R.L., Silva, G.A.P., Verardi, C.K., Oliveira, A.L.B., Gonçalves, P.S. 2013a. “Rubber tree early selection for yield stability in time and among locations”. Euphytica, 189(3): 343-350, DOI: http://doi.org/10.1007/s10681-013-0874-6.
). La comparación de medias se realizó mediante la prueba de Tukey-Kramer (α = 0,05). - Fase 2: Estabilidad y alometría: Se aplicó el modelo de Finlay y Wilkinson (1963)Finlay, K.W., Wilkinson, G.N. 1963. “The analysis of adaptation in a plant-breeding programme”. Australian Journal of Agricultural Research, 14(6): 742-754, DOI: http://doi.org/10.1071/AR9630742.
, para estimar la respuesta genotípica (β) y la predecibilidad (R²). Asimismo, se evaluó la eficiencia alométrica mediante regresión lineal (RS vs. DAP), para analizar la partición de carbono hacia la síntesis de látex. Finalmente, se calcularon coeficientes de correlación de Pearson (r) y Spearman (ρ), para validar la asociación entre las variables biológicas y los factores ambientales, complementando los criterios de selección para estabilidad productiva (Gonçalves et al., 2011Gonçalves, P. de S., Scaloppi Júnior, E.J., Martins, M.A., Moreno, R.M.B., Branco, R.B.F., Gonçalves, E.C.P. 2011. “Assessment of growth and yield performance of rubber tree clones of the IAC 500 series”. Pesquisa Agropecuária Brasileira, 46(12): 1643-1649, DOI: http://doi.org/10.1590/S0100-204X2011001200009.
).
Resultados
⌅Vigor estructural y precocidad vegetativa
⌅El análisis de varianza (Tabla 1), reveló diferencias altamente significativas (P < 0,001), para la tasa de incremento diamétrico (I-DAP) y altura (I-ALT), entre los genotipos evaluados. El clon IAN 873 presentó la mayor tasa de incremento diamétrico, con una media de 2,244 mm mes⁻¹, posicionándose significativamente por encima de los materiales de la serie IRCA. En términos de proyección temporal, los datos indican que el IAN 873 reduce la fase juvenil entre 10 y 12 meses en comparación con el clon IRCA 111. Como se observa en la figura 1, las trayectorias de crecimiento muestran una divergencia clara desde el inicio del ciclo, consolidando al IAN 873 como el genotipo con mayor dinamismo estructural.
| Clon | Vigor (Incrementos) | Desarrollo Final | Rendimiento Mensual | |||
|---|---|---|---|---|---|---|
| I-DAP (mm/mes) ¹ | I-ALT (m/mes) ² | DAP_F (cm)³ | ALT_F (m)⁴ | REN (g/mes) ⁵ | R_GAS (g/serie) ⁶ | |
| IRCA 41 | 1,982b | 0,237b | 10,32b | 11,54bc | 39,18a | 21,31a |
| IAN 873 | 2,244a | 0,252a | 11,94a | 12,41a | 34,55b | 18,66b |
| IAN 710 | 2,155a | 0,233b | 11,49a | 11,65b | 29,48c | 15,77c |
| IRCA 19 | 1,928bc | 0,222c | 10,18b | 10,84c | 26,23d | 13,84d |
| IRCA 111 | 1,833c | 0,224c | 9,61b | 11,07bc | 24,60d | 13,64d |
| n | 18,609 | 16,953 | 399 | 400 | 6,410 | 6,403 |
| Valor F | 27,60*** | 13,96*** | 24,85*** | 15,00*** | 61,95*** | 53,84*** |
| Varianza (σ2) | 3,81 | 0,035 | 301,16 | 1,95 | 750,00 | 258,41 |
Medias
con letras distintas en la misma columna indican diferencias
significativas (Tukey-Kramer, P≤ 0,05). Asteriscos (***) indican P <
0,001. Errores estándar (EE): ¹0,032, ²0,003, ³0,19, ⁴0,15, ⁵0,77,
⁶0,45.
Dinámica de rendimiento y eficiencia biofísica
⌅La evaluación productiva mostró una divergencia clara entre el vigor vegetativo y el rendimiento de caucho seco (Tabla 1). Aunque el IAN 873 lideró el vigor estructural, el clon IRCA 41 presentó el rendimiento más alto, con un promedio de 39,18 g mes⁻¹ (P < 0,0001). Esta diferencia evidencia una especialización funcional distinta: mientras los clones IAN privilegiaron la biomasa leñosa, los clones IRCA, particularmente el IRCA 41, maximizaron la producción de látex por unidad de tiempo.
Análisis de correlaciones
⌅La matriz de correlación (Tabla 2), indicó una asociación positiva significativa entre el DAP y el rendimiento acumulado (r = 0,657***). Por el contrario, se detectó una correlación negativa entre el incremento en altura (I-ALT) y el rendimiento (r = -0,105***). Respecto a las variables ambientales, el déficit hídrico (DEFH), presentó una correlación negativa y significativa con la productividad (r = -0,186***), estableciéndose como el modulador ambiental con mayor impacto negativo sobre la turgencia y el flujo de látex. Esta relación entre las variables climáticas y la eficiencia productiva se visualiza en la figura 2, donde el estrés hídrico limita la respuesta productiva de los genotipos más sensibles.
| Variable de Respuesta | DEFH (mm) | TMAX (°C) | PTME (mm) | DAP (cm) | ALT (m) | I-DAP |
|---|---|---|---|---|---|---|
| REN (g/mes) | -0,116*** | 0,132*** | 0,231*** | 0,584*** | 0,466*** | -0,165*** |
| R_GAS (g/serie) | -0,186*** | 0,113*** | 0,322*** | 0,657*** | 0,604*** | -0,105*** |
| I-DAP (mm/mes) | -0,194*** | 0,042*** | 0,151*** | -0,035*** | -0,070*** | - |
| I-ALT (m/mes) | -0,103*** | 0,116*** | 0,133* | -0,086*** | -0,043*** | 0,156*** |
DAP:
diámetro; ALT: altura; I: incrementos; REN: rendimiento mensual
promedio por mes; R_GAS: rendimiento mensual promedio por serie por
árbol. DEFH: déficit hídrico; TMAX: temperatura máxima; PTME:
precipitación. Significancia: *** P < 0,001; * P < 0,05; ns: no significativo.
Estabilidad y eficiencia biológica (β)
⌅El modelo de Finlay-Wilkinson clasificó a los genotipos según su respuesta ambiental (β). El clon IRCA 41 obtuvo un coeficiente β = 1,226, indicando una alta sensibilidad y maximización productiva en condiciones óptimas. El clon IAN 873 mostró un comportamiento más estable (β = 1,118, R² = 0,751), manteniendo un desempeño consistente ante la variabilidad del entorno. En contraste, el clon IRCA 111 presentó una respuesta negativa ante el estrés hídrico y térmico (β = -0,060***, β2 = -6,956***), denotando una alta vulnerabilidad ambiental (Tabla 3). La figura 3 permite observar la dispersión de los genotipos bajo el modelo de regresión conjunta, confirmando la posición superior del IRCA 41 frente a la estabilidad mostrada por el IAN 873.
| Clon | Estabilidad (F-W) | Sensibilidad ambiental (Betas) | ||||
|---|---|---|---|---|---|---|
| Coef. β | R2 | DEFH (β1) | TMAX (β2) | EVPR (β3) | Caracterización | |
| IRCA 41 | 1,226 | 0,599 | 0,000 (ref) | 0,000 (ref) | 0,000 (ref) | Especialista Élite |
| IAN 873 | 1,118 | 0,751 | -0,004ns | -3,665* | 0,329** | Homeostasis |
| IAN 710 | 0,889 | 0,691 | -0,014ns | -4,230** | 0,354** | Intermedio |
| IRCA 111 | 0,921 | 0,560 | -0,060*** | -6,956*** | 0,505*** | Vulnerable |
| IRCA 19 | 0,859 | 0,606 | -0,091*** | -3,589* | 0,279* | Sensible al estrés |
| Valor F (ANOVA) | 37,57* | - | 61,77* | 35,29* | 7,85* | Interacción GxA* |
β
> 1,0: Adaptabilidad específica a ambientes favorables; β≈1,0:
Adaptabilidad general; β < 1,0: Adaptabilidad específica a ambientes
desfavorables (bajo dinamismo). Valores negativos indican
susceptibilidad ante el déficit hídrico (DEFH) y temperatura máxima
(TMAX). EVPR: evaporación. Significancia: *** P < 0,001; ** P <
0,01; * P < 0,05; ns: no significativo.
Validación de la selección temprana
⌅El análisis de la interacción Genotipo x Periodo (F = 7,80***), mostró una variabilidad inicial alta en el rendimiento. No obstante, el coeficiente de correlación de Spearman (ρ), indicó una estabilización del ranking de genotipos a partir del tercer año (ρ = 0,679***; Tabla 4), confirmando que la jerarquía productiva de los clones se consolida tras los 36 meses de evaluación. Dicha estabilización y el comportamiento del ranking a través del tiempo se presentan en la figura 4, la cual corrobora que la selección temprana a los 36 meses es un predictor confiable del potencial productivo final.
| Clon | P1 (Año 1) | P2 (Año 2) | P3 (Año 3) | P4 (Año 4) | ρ (P vs P4) |
|---|---|---|---|---|---|
| IRCA 41 | 3,05b | 15,18a | 28,03a | 41,46a | 0,291*** (P1 vs P4) |
| IAN 873 | 3,31b | 13,15c | 23,46a | 37,83ab | |
| IRCA 111 | 1,48b | 8,99d | 16,16c | 31,70bc | 0,335*** (P2 vs P4) |
| IAN 710 | 3,39b | 12,87c | 18,88b | 30,34c | |
| IRCA 19 | 1,81b | 12,77c | 16,65c | 25,36d | 0,679*** (P3 vs P4) |
| EE (promedio) | 1,40 | 1,27 | 1,27 | 1,28 | |
| Valor F (ANOVA) | Clon: 13,46*** | Periodo: 677,72*** | Interacción: 7,80*** |
Medias
de mínimos cuadrados (LS-Means). Letras distintas en la misma columna
indican diferencias significativas (Tukey-Kramer, P ≤ 0,05). ρ de
Spearman: Coeficientes de correlación de rangos respecto al año final.
Efectos fijos del ANOVA (F): Clon (13,46***), Periodo (677,72***), Interacción (7,80***). Significancia: *** P < 0,001.
Discusión
⌅Plasticidad fenotípica y compensación metabólica
⌅Los
resultados confirman que el vigor vegetativo temprano es un predictor
sólido de la adaptabilidad al sitio. La superioridad del IAN 873 en
I-DAP corrobora la alta heredabilidad del crecimiento diamétrico
reportada por Gonçalves et al. (2005a)Gonçalves,
P. de S., Bortoletto, N., Cardinal, Á.B.B., Gouvêa, L.R.L., Costa,
R.B., Moraes, M.L.T. 2005. “Age-age correlation for early selection of
rubber tree genotypes in São Paulo State, Brazil”. Genetics and Molecular Biology, 28(4): 758-764, DOI: http://doi.org/10.1590/S1415-47572005000500018.
,
validando la viabilidad de la selección precoz mediante el seguimiento
de tasas de incremento. La divergencia observada entre el crecimiento
(IAN 873) y la productividad (IRCA 41), sustenta la hipótesis de un trade-off o compensación funcional: una canalización preferencial de carbono, ya
sea hacia biomasa leñosa o hacia la biosíntesis de polímeros (látex),
tal como lo ha señalado Priyadarshan (2011)Priyadarshan, P.M. 2011. Biology of Hevea Rubber. Wallingford, UK: CABI, 350 p.
.
Esta eficiencia biofísica del IRCA 41 lo posiciona como un material de
élite para esquemas de pica intensiva, reforzando la importancia de los
modelos de selección simultánea propuestos por Gouvêa et al. (2013b)Gouvêa,
L.R.L., Silva, G.A.P., Verardi, C.K., Oliveira, A.L.B., Gonçalves, P.S.
2013b. “Simultaneous selection of rubber yield and girth growth in
young rubber trees”. Industrial Crops and Products, 50: 39-43, DOI: http://doi.org/10.1016/j.indcrop.2013.06.040.
.
Respuesta al estrés hídrico y estabilidad clonal
⌅La correlación negativa observada entre el déficit hídrico (DEFH) y la producción corrobora los hallazgos de Ajith et al. (2022)Ajith, S.V., Nair, N.N., Sathik, M.B.M., Meenakumari, T. 2022. “Breeding for low-temperature stress tolerance in Hevea brasiliensis: screening of newly developed clones using latex biochemical parameters”. Journal of Rubber Research, 24(1): 63-76, ISSN: 2510-1498, DOI: http://doi.org/10.1007/s42464-021-00143-2.
, y Jacob et al. (2021)Jacob, J., Gitz, V., Gohet, E., et al., 2022. “Natural rubber contributions to adaptation to climate change”. In: XV World Forestry Congress, 2-6 de mayo, Seúl, Corea del Sur. (Nota: En congresos suele aceptarse et al. si el comité es muy extenso).
.
La sensibilidad de la turgencia laticífera al potencial hídrico es un
factor limitante crítico en el trópico húmedo; bajo estrés, la
viscosidad del látex aumenta, lo que incrementa la tasa de obturación
del panel de pica. La estabilidad clonal, por tanto, no es solo un valor
estadístico, sino una propiedad fisiológica intrínseca; autores como Omokhafe y Alika (2003)Omokhafe, K.O., Alika, J.E. 2003. “Clonal stability of latex yield in eleven clones of Hevea brasiliensis (Muell. Arg.)”. Genetics and Molecular Biology, 26(3): 313-317, DOI: http://doi.org/10.1590/S1415-47572003000300016.
, y Omokhafe (2004)Omokhafe, K.O. 2004. “Clonal stability of tree dryness in Hevea brasiliensis Muell. Arg.”. Genetics and Molecular Biology, 27(2): 242-244, DOI: http://doi.org/10.1590/S1415-47572004000200019.
,
han señalado que la susceptibilidad al secamiento es una variable
determinante para la estabilidad de la producción de látex. El
comportamiento "Especialista Élite" del IRCA 41 versus la homeostasis
ambiental del IAN 873 subraya la necesidad de diversificar las
recomendaciones técnicas: el IAN 873 reduce el riesgo en regiones con
marcada estacionalidad, mientras que el IRCA 41 es superior en ambientes
controlados o con menor estrés hídrico.
Implicaciones para el mejoramiento genético
⌅El
hallazgo central de este estudio es la validación del umbral de
selección a los 36 meses. El "ruido" fisiológico observado en la etapa
P1, derivado de la dependencia del vigor del sitio, disminuye
drásticamente a medida que el sistema laticífero madura. La
consolidación del ranking productivo (ρ = 0,679***), coincide con lo planteado por Chandrasekhar et al. (2007)Chandrasekhar,
T.R., Marattukalam, J.G., Mercykutty, V.C., Priyadarshan, P.M. 2007.
“Age of yield stabilization and its implications for optimising
selection and shortening breeding cycle in rubber (Hevea brasiliensis Muell. Arg.)”. Euphytica, 156(1-2): 67-75, DOI: http://doi.org/10.1007/s10681-006-9352-8.
,
respecto a la criticidad del tercer año de pica. Se concluye que la
integración de modelos lineales mixtos (REML/BLUP), en el monitoreo
multianual permite reducir sustancialmente los ciclos de mejoramiento,
proporcionando una herramienta robusta y estadísticamente confiable para
la selección de genotipos en condiciones agroecológicas del trópico
húmedo mexicano (Gouvêa et al., 2013aGouvêa,
L.R.L., Silva, G.A.P., Verardi, C.K., Oliveira, A.L.B., Gonçalves, P.S.
2013a. “Rubber tree early selection for yield stability in time and
among locations”. Euphytica, 189(3): 343-350, DOI: http://doi.org/10.1007/s10681-013-0874-6.
).
Conclusiones
⌅Esta investigación demuestra que el origen genético es el determinante principal de la precocidad estructural en Hevea brasiliensis; genotipos como IAN 873 e IAN 710 permiten reducir el periodo improductivo hasta en 12 meses, acelerando significativamente la viabilidad económica de la plantación.
Se identificó una divergencia metabólica fundamental: mientras los materiales de la serie brasileña priorizan el vigor vegetativo, el clon IRCA 41 canaliza los asimilados preferencialmente hacia la biosíntesis de látex, consolidándose como el ideotipo de alta eficiencia productiva bajo condiciones óptimas. Esta capacidad está supeditada a la resiliencia ambiental, donde el déficit hídrico actúa como el modulador más crítico; ante este factor, la homeostasis observada en el IAN 873 lo posiciona como un material robusto frente a la estacionalidad climática.
Finalmente, la validación estadística de la jerarquía productiva a partir del tercer año (36-48 meses), constituye un hito metodológico para la región. Este hallazgo permite optimizar los programas de mejoramiento al reducir drásticamente los ciclos de evaluación, pasando de protocolos tradicionales de 15 años a esquemas de selección precoz altamente confiables y escalables para el trópico húmedo.