Tesis de Soya
Tesis de Soya
Tesis de Soya
Directores:
GUSTAVO LIGARRETO M., Ph.D.
Universidad Nacional Sede Bogot
JOS DLMER MORENO M., Ph.D
Corpoica, C.I. Tibaitat.
AGRADECIMIENTOS
A Corpoica, con la direccin ejecutiva del Dr. Arturo Vega por haberme dado la oportunidad de
realizar los estudios doctorales y financiar los mismos.
Al Dr. Jaime Triana R., director del C.I La Libertad de Corpoica por su confianza y apoyo.
A los Drs. Gustavo Ligarreto y Jos Dlmer Moreno por su direccin y asesora en mis estudios
doctorales y trabajo de tesis.
A los Drs. Jorge Franco y Alia Rodrguez por ser parte del comit asesor de mis estudios
doctorales.
A mis compaeros M.Sc. Carmen Rosa Salamanca S. y Samuel Caicedo G. por su
colaboracin y ayuda en el desarrollo de las actividades de campo.
Al M.Sc. Victor Nnez, a las I.A Linda Gmez y Erika Snchez B., por su colaboracin y apoyo
en el desarrollo de est investigacin en el laboratorio de Genetica Molecular Vegetal de
Corpoica, C.I. Tibaitat.
A las M.Sc. Margarita Ramrez y Ruth Bonilla y a la biloga Lina Moreno del laboratorio de
microorganismos de Corpoica por su apoyo en la consecucin y aislamiento de las cepas de B.
japonicum.
Al personal tcnico y de apoyo de Corpoica del C.I. La Libertad, en particular a Luis Vega,
Graciela Casiano, Ivn Chavarro, Mesas Moreno, Maritza Prado, a los pasantes Andrs
Snchez y Julio Quiroga.
Al M.Sc. Oscar Pardo y sus colaboradores del laboratorio de nutricin animal del C.I. La
Libertad, por su apoyo en la determinacin de protena de grano.
A M.Sc. Mauricio lvarez e
I.A
georeferenciada.
Al Dr. Pengyin Chen por haberme dado la oportunidad de realizar la intership en el
departamento de mejoramiento gentico de soya de la Universidad de Arkansas y por compartir
conmigo sus conocimientos.
A la Universidad Nacional de Colombia, por admitirme en sus claustros de estudio, orientarme y
fortalecer mi formacin acadmica.
A todos aquellos que de una u otra manera contribuyeron con la realizacin de mis estudios y
trabajo de tesis doctoral.
DEDICATORIA
A mis
padres Salomn Valencia y Marina Ramrez, por su amor, sus bendiciones y por
haberme hecho un hombre de bien. Gratitud por siempre para mis amados padres.
A mi esposa Elizabeth y a mis hijas Lizth Alexandra, Diana Vanessa y Jessica a quienes amo,
por su comprensin, paciencia, tolerancia y apoyo incondicional.
A mis hermanos que tanto quiero, por hacer de los Valencia una familia de hombres y mujeres
de grandes principios y valores, familia pujante y emprendedora, guiada por unos maravillosos
padres.
CONTENIDO
Pg
RESUMEN
ABSTRACT
1. INTRODUCCIN
10
12
14
especificidad
1.4.3. Diversidad de rizobios a nivel molecular
16
17
19
21
21
24
26
ambiente
29
analysis).
1.6 Propsitos de la investigacin en soya.
2. OBJETIVOS
33
35
35
35
3. MATERIALES Y MTODOS
36
3.1 Localizacin.
36
39
40
41
42
44
45
46
48
51
51
microsatlites.
3.7.2. Caracterizacin molecular de cepas de B. japonicum mediante
53
rep-PCR.
4. RESULTADOS Y DISCUSIN
4.1 Efecto de la asociacin simbitica de cepas de B. japonicum sobre el
55
55
55
58
60
63
vainas.
4.1.5. Eficiencia de la fijacin biolgica del nitrgeno.
69
74
76
87
87
90
95
de grano.
4.3.1. Respuesta de variedades de soya por cepa de B. japoniucum, a
95
99
103
104
109
111
113
114
116
118
119
agronmicas.
4.3.4.2.
Asociacin
de
PC1
ambientales
las
covariables
121
125
125
edafoclimticas.
mediante microsatlites.
4.4.2 Caracterizacin molecular de las cepas de B. japonicum J-01, J-
128
132
6 RECOMENDACIONES
134
7 BIBLOGRAFA
135
8 ANEXOS
154
NDICE DE TABLAS
Tabla 1
Tabla 2
Tabla 3
37
Tabla 4
38
Tabla 5
39
Tabla 6
47
Tabla 7
50
Tabla 8
61
Tabla 9
64
Tabla 10
66
Tabla 11
67
Tabla 12
68
Tabla 13
69
Tabla 14
76
Tabla 15
78
Tabla 16
79
Tabla 17
82
Tabla 18
84
Tabla 19
85
Tabla 20
86
Tabla 21
88
Tabla 22
89
Tabla 23
91
Tabla 24
94
Tabla 25
96
Tabla 26
99
Tabla 27
104
Tabla 28
109
Tabla 29
113
Tabla 30
114
Tabla 31
118
Tabla 32
120
Tabla 33
125
NDICE DE FIGURAS
Figura 1
Figura 2
10
Figura 3
33
Figura 4
43
Figura 5
46
Figura 6
57
Figura 7
59
Figura 8
60
Figura 9
62
Figura 10
70
Figura 11
71
Figura 12
72
Figura 13
81
Asociacin del nmero de vainas con 3 granos (V3), con vainas totales
(VT) y rendimiento de grano, a travs de ambientes y cepas de B.
japonicum -N.
92
Figura 15
97
Figura 16
Anlisis biplot GGA del rendimiento de grano de una variedad de soya 101
inoculada con diferentes cepas de B. japonicum y evaluadas en 8
ambientes de la Orinoquia colombiana. A. Soyica P-34, B. Orinoquia 3,
C. Corpoica Libertad 4, D. Corpoica Taluma 5, E. Corpoica superior 6 y
F. Corpoica Sabana 7.
Figura 17
Figura 18
Figura 19
Figura 20
Figura 21
Biplot GGA con vector virtual CP1v y CP2v promedio. Los Vc en el 115
crculo representan los ms estables y de mayor rendimiento. Las lneas
punteadas paralelas al eje virtual Cp2v, representan la proyeccin en
estabilidad.
Figura 22
Biplot GGA con vector virtual CP1v y CP2v sobre el ambiente Lib08a 117
ms discriminante y representativo. Los Vc ordenados de derecha a
izquierda sobre el eje x representan el ranqueo de mayor a menor por
rendimiento de grano sobre este ambiente.
Figura 23
con el 106
diferencias genticas que pueden ser de utilidad para futuros programas de identificacin de
individuos, estudios poblacionales y seleccin asistida. Las cepas de B. japonicum J-96 y J-98
se constituyen en nuevas opciones de fijacin biolgica de nitrgeno altamente efectivas para
las variedades de soja, en la Orinoquia colombiana.
Palabras claves: Oxisoles, soya, FBN, mejoramiento gentico.
, in eight environments. Strain J-96 presented greater number and total dry weight of nodules
followed by the J-98 strain, both strains contributing in 93 and 92% of N required, respectively.
The J-01 strain had the lowest fixation of N, with 84% of N requirement. Grain yield varied
across environments and was influenced by the strain and soybean variety. The main effects of
variety x strain (CP1) of SRGE model were high and positively correlated with grain yield (r =
0.97, P <0.01). The varieties more stable and productive were Corpoica Superior 6 and
Orinoquia 3 with J-96 strain. A consistent association was observed between CP1 environments
with water balance (r=0.92, P<0.01) in the period R3-R6. The environmental most discriminating
and representative for selection process was C.I. La Libertad in the Piedemonte. Varietal
molecular patterns defined from microsatellite (SSR) and strains patterns derived through BoxPCR and Eric-PCR primers, allowed the establishment of genetic differences that may be useful
for future individual identification program, genetic diversity studies and assisted selection. The
J-96 and J-98 strains are new options for its high infective capacity and its high biological
nitrogen fixation affectivity in soybean varieties in the colombian Orinoquia natural region.
Key words: Oxisoils, soybean, NBF, plant breeding.
1. INTRODUCCIN
1.1 Generalidades y potencialidades de la Orinoquia colombiana
El incremento creciente en la demanda de alimentos alrededor del mundo es un gran reto para
los fitomejoradores en busca de plantas de alto potencial gentico que mejoren la productividad
de los suelos con limitaciones para la produccin y que puedan ser vinculadas en sistemas
agrcolas sostenibles.
En Suramrica, cerca de 250 millones de hectreas en sabanas tropicales presentan problemas
de toxicidad de Al, incluidas las sabanas de la Orinoquia colombiana (Vera, 2000). La regin
natural de la Orinoquia colombiana representa el 23% del rea total del territorio nacional
equivalente a 26 millones de hectreas (Corpoica, 2009). En esta vasta regin, se han definido
dos subregiones: el piedemonte llanero y las sabanas no inundables, de alto inters para el
desarrollo agropecuario sostenible. Las sabanas no inundables, se encuentran ubicadas al sur
oriente de los ros Meta y Metica hasta hacer contacto con el Andn del Orinoco, la cual a su
vez, se subdivide en dos unidades: la Altillanura plana y la Altillanura disectada. La Altillanura
plana, como paisaje particular, se sita al sur del ro Meta y comprende buena parte del
departamento del Vichada (Molano, 1998).
La Altillanura plana y las terrazas del piedemonte son de gran inters para el sector
agropecuario por presentar ventajas comparativas como topografa de fcil mecanizacin,
excelente rgimen de lluvias y cercana a los principales centros de consumo del pas
(Corpoica, 2009). Los suelos oxisoles (Tropeptic Haplustox isohipertermic) de estos
agroecosistemas presentan un pH promedio de 4.5, baja disponibilidad
de P, Ca, Mg, K
(Amzquita et al., 2005; Rivas et al., 2004; Molina et al., 2003), tienen baja capacidad de
intercambio catinico, la materia orgnica tiende a ser baja y alta fijacin de los fosfatos
(Borkert y Sfredo, 1995) y poseen serios problemas de fitotoxicidad de aluminio (Bianchi et al.,
2000), que limita el crecimiento de las plantas y la productividad agrcola (Von y Mutert, 1995).
Sin embargo, estudios realizados por Rivas et al. (2004), resaltan los beneficios econmicos y
ambientales de la produccin en la Altillanura colombiana.
desde 1954 lideran la produccin mundial (Ridner, 2006). En Brasil fue introducida en 1882,
pero su difusin se produjo a principios del siglo XX (Purseglove, 1984). En 1929 la soya fue
introducida a Colombia a nivel experimental y su produccin comercial se inici en el siglo XX,
dcada de los 50's, con variedades introducidas de Estados Unidos (Bastidas y Agudelo, 1994).
En Colombia, el grano de soya (Glycine max (L) Merril) es fuente fundamental de protena como
materia prima para la elaboracin de alimentos balanceados, particularmente para la cadena
avcola-porccola, con usos alternativos como ensilaje en la cadena crnica-lctea y, alimento
rico en aminocidos esenciales, con grandes beneficios para la salud y nutricin humana. La
soya es el cultivo por excelencia de los sistemas de rotacin con gramneas como el arroz, el
maz y el sorgo, base de la seguridad alimentara del pas (Valencia et al., 2006).
1.2.1
(2009), fue de 90,2 millones de hectreas, para una produccin global de 220,5 millones de
toneladas y un rendimiento promedio de grano de 2,4 t.ha-1. Estados Unidos con el 33 % de la
produccin total ocupa el primer lugar, seguido por Brasil con el 26,2% y luego Argentina con
21,5%. Colombia solo alcanza a producir cerca del 0,023% del total mundial (Tabla 1). En el
mercado internacional, el mayor importador es China con 34 millones de toneladas seguido por
la Unin Europea con 15,45 millones, Japn 4,15 millones y Mxico con 3,95 millones de
toneladas (USDA 2007)
Tabla 1. rea, produccin y rendimiento de los principales productores de soya, en el mundo.
2007.
Pas
Estados Unidos
Brasil
Argentina
China
India
Paraguay
Canad
Colombia
Total
Produccin (t.)
72.860.400
57.857.200
47.482.784
13.800.147
10.968.000
5.856.000
2.695.700
51.471
220.532.612
Rendimiento t.ha-1
2,80
2,80
2,97
1,55
1,24
2,41
2,30
2,20
2,44
La soya es el sexto cultivo con mayor produccin en el mundo superado por caa de azcar,
maz, trigo, arroz y papa. Es reconocida como un cultivo sensible al fotoperodo y que
originalmente est adaptado a latitudes templadas con ms horas de luz (14-16 h); en este
sentido el mejoramiento gentico ha tenido una importante participacin, ya que ha permitido
aumentar el rea de siembra en los trpicos con la identificacin e incorporacin del carcter
juvenil, que confiere insensibilidad al fotoperodo (Alezones y Zocco, 2007).
Colombia a pesar de tener condiciones agroecolgicas ideales para la produccin de soya,
importa cerca del 93% de sus necesidades, cercanas al milln de toneladas anuales. Los bajos
rendimientos de grano y los altos costos de produccin son los dos elementos que coexisten
para que el panorama de la produccin de soya en el pas sea desalentador y la necesidad de
importar el grano y torta (MADR, 2008), para suplir parcialmente las necesidades internas. Las
importaciones de soya en el 2007 fueron de 332.064 toneladas de grano y 810.605 toneladas
de torta (FAO, 2009).
La reduccin del rea de siembra ha variado en los ltimos 20 aos en Colombia, pasando de
116.150 hectreas en 1990, a 26.344 hectreas en 2008, con una produccin de 55.950
toneladas (Agronet, 2009). En la Figura 1, se puede observar el comportamiento promedio del
rea y la produccin de soya en Colombia desde 1987, donde se visualiza una drstica cada
del rea de siembra y produccin de soya desde 1990. Los altos costos de produccin e el
Valle del Cauca convirtieron al departamento del Meta en el mayor productor del grano en
Colombia al aportar cerca del 76.1% de la produccin nacional (Agronet, 2009), por su ventajas
comparativas. El rea de siembra de soya en la Orinoquia colombiana en el 2009 fue de 23.480
hectreas, con una participacin de la Altillanura del 53 %
publicada Corpoica y Cadenas productivas del Meta).
Rendimiento (Kg.ha)
2600
2400
2200
2000
1800
1600
1400
1200
1000
rea (ha)
Produccin (t)
Area (has)
120000
250000
200000
100000
80000
150000
60000
100000
40000
50000
20000
19
8
19 7
88
19
8
19 9
90
19
9
19 1
92
19
9
19 3
94
19
9
19 5
96
19
9
19 7
98
19
9
20 9
00
20
0
20 1
02
20
0
20 3
0
20 4
0
20 5
06
20
0
20 7
08
Produccin (t)
140000
Ao
Semestre A
Semestre B
Total
2000
57
57
2001
50
630
680
2002
175
1.625
1.800
2003
1.270
2.822
4.092
2004
3.495
6.300
9.795
2005
1.471
3.100
4.571
2006
2.100
3.800
5.900
2007
1.050
3.150
4.200
2008
3.250
6.415
9.665
2009
5.953
6.481
12.434
Mercado y precios de la soya. Los precios de esta materia prima han sido muy
tonelada, en promedio, durante el ltimo trimestre de 2005, para igual perodo en el 2007 fue de
427. Por su parte, la torta de soya pas, en los mismos perodos de la referencia, de 254 a 333
dlares en el 2007 (Portafolio, 2008). El precio ms alto registrado para la soya se dio en el
2008 de hasta 565 dlares (Figura 2), desbordando su fluctuacin histrica que bordeaba los
250 dlares. Aunque esta fue una buena seal de precios para los productores de materias
primas, los avicultores vieron amenazados sus negocios. En este escenario, las condiciones
reales de mercado perdieron sustento en trminos de oferta y demanda, donde primaron los
resultados de corto plazo. En el primer semestre de 2009, los precios de soya se redujeron en
un 23,6% en la bolsa de Chicago y los fletes en un 50%, aunque con una revaluacin fluctuante
del peso frente al dlar. En promedio, en el primer semestre de 2009, la la torta de soya registr
una variacin de 3% y un incremento de 16% (Fenavi, 2009).
600
500
400
300
200
100
0
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
Ao
Figura 2. Precio promedio (US$ -CIF) por tonelada de soya primer semestre (2000-2009).
Fuente: Fenavi, 2009.
promisorias se realizaron con utilizacin de la cepa ICA J-01 desde su liberacin en 1989,
desconciendo la importancia de la interaccin especfica rizobio-variedad-ambiente. Salamanca
et al. (2006), referencia de manera prelimiar el comportamiento de algunas variedades de soya
a diferentes cepas de B. japonuicum.
crecimiento, productoras de lcali. Las colonias son circulares, rara vez translcidas, blancas y
convexas con un dimetro menor a 1 mm despus de 5-7 das de incubacin a una temperatura
ptima de 25 a 30C (Jordan, 1982). El gnero Bradyrhizobium tiene un amplio rango de
plantas huspedes, incluyendo muchas leguminosas del trpico y de la zona templada (Ba et
al., 2004). En el mundo, han sido reportadas las bacterias B. japonicum y B. elkanii, fijadoras de
nitrgeno en soya (Glycine max (L.) Merril), como las ms predominantes (Chao et al., 2008).
La soya tiene la habilidad de asociarse en forma simbitica con estas bacterias para obtener a
travs de la FBN gran parte del N que requiere para su crecimiento (Hungria et al., 1998;
Perticari et al., 2003), con aportes hasta del 94% del N, en los cultivares ms productivos
(Hungria et al., 2006).
La soya est genticamente diseada para acumular protenas en sus granos, razn por la cual
se lo considera como cultivo altamente demandante de nitrgeno. Las protenas de soya tienen
un 5,75 % de nitrgeno en su composicin, por lo cual el cultivo debe acumular alrededor de 80
Kg. de este nutriente en su biomasa vegetal para producir cada tonelada de grano (Gonzlez,
2002).
Los cultivos de soya con un rendimiento de 2.500 Kg.ha-1 acumulan cerca de 200 Kg. de N y del
67 al 75% de este N se encuentra almacenado en el grano (Boddey y Hungra, 1997). Las
fuentes de N requerido para el crecimiento de la planta proviene de la fijacin simbitica del N
atmosfrico (N2), de la asimilacin directa de N inorgnico del suelo (Cattelan y Hungria, 1995),
y aqul movilizado desde rganos de acumulacin temporaria en la propia planta (Luyindula &
Weaver, 1989). La relacin entre el nitrgeno absorbido desde el suelo y el proveniente de la
simbiosis vara con la especie de leguminosa, las condiciones ambientales y las caractersticas
genticas del microsimbionte (Shabaev et al., 1996). En primera instancia, mientras se genera
el sistema nodular, la planta utiliza nitrgeno del suelo. Si hay mucho nitrgeno en el ambiente
de la raz, el nmero de ndulos formados sera menor (Gonzlez, 2002). Las plantas absorben
el N del suelo en forma de nitrato o amonio, donde el exceso de nitrgeno en el suelo en forma
de nitratos, tiene un efecto inhibitorio sobre la simbiosis en todos los pasos, desde la infeccin,
formacin del ndulo y la fijacin de N2 Fernndez (2003). As, el aporte global por FBN para la
soya es menor en suelos bien provistos de N que en aqullos en que el nutriente es deficitario
(Gonzlez, 2002; Pietrarelli et al., 2008).
Por lo anterior, la FBN, tiene grandes repercusiones sobre la produccin y calidad del grano,
reduccin de costos y sobre la sostenibilidad del sistema de produccin. La herencia del
proceso de nodulacin puede ser relativamente simple. Segn Tsyganov et al. (1998), los
genes involucrados en el control del proceso de nodulacin son recesivos y pueden no ser
numerosos.
1.4.1 Proceso de infeccin y nodulacin. En la simbiosis, la planta husped tiene la
informacin gentica para la infeccin simbitica y para la nodulacin. El papel de la bacteria es
slo el de disparar el proceso. Este proceso de asociacin simbitica centra su importancia en
dos aspectos fundamentales de la relacin
Nod estimulan la divisin de las clulas vegetales, produciendo finalmente el ndulo. 6. Ingreso
de las bacterias a las clulas de la raz y diferenciacin de los bacteroides. 7. Establecimiento
del ndulo funcional maduro: El ndulo comienza a fijar nitrgeno slo cuando termina el
desarrollo de los bacteroides. El sistema vascular de la planta se extiende dentro del ndulo y
transporta nutrientes hacia y desde el ndulo. Se sintetizan enzimas especficas que participan
en este proceso simbitico, entre ellas:
-Nitrogenasa: enzima es responsable de convertir o reducir el nitrgeno atmosfrico (N2) en
NH4+, es sintetizada en el citosol por los bacteroides. Las leguminosas utilizan el NH4+, para
convertir precursores metablicos en aminocidos que son convertidos en protenas. Esta
protena de gran tamao que consiste de dos componentes, la protena homodimrica que
contiene Fe y es codificada por nifH, y la protena tetramrica que contiene Fe y molibdeno
(Mo), codificada por los genes nifD y nifK (Wang et al., 2001).
-Leghemoglobina: es una protena donde la subunidad hemina es formada por los
bacteroides, mientras que la globina es formada por la planta. La leghemoglobina es
responsable del color rojo de los ndulos activos. Sirve para limitar la concentracin de oxgeno
dentro del ndulo. De esta forma los bacteroides reciben suficiente oxgeno para sobrevivir y se
evita que la nitrogenasa pueda ser inactivada por el oxgeno (Olivares y Barea, 1995).
-Hidrogenasa: La nitrogenasa cataliza la reduccin de nitrgeno atmosfrico a amonio y este
proceso conlleva la produccin de hidrgeno. Este hidrgeno se pierde como gas por lo que
constituye una prdida de energa en el proceso de fijacin de nitrgeno (Ruz et al., 2000). El
reciclado de hidrgeno se ha asociado a un aumento de la productividad de las plantas debido
al incremento de la eficiencia energtica en la simbiosis rizobios-leguminosas. La hdrogenasa
es una enzima heterodimrica que esta unida a la membrana y contiene un centro activo
bimetlico NiFe (Brito et al., 1994).
Aunque el proceso de nodulacin es variable entre genotipos, ambientes y cepas, segn
Imsande (1998), la nodulacin puede ser observada de 3 a 5 das de la emergencia de la
planta cuando las condiciones son ideales para la expresin de la simbiosis y la actividad de
fijacin de 7 a 10 das despus de la visualizacin de los ndulos. Los valores de N fijado son
bajos en los estados vegetativos hasta comienzos de floracin, con valores entre 0,5 y 1 Kg.
de N fijado.ha-1.da-1. De ah en adelante, se registra la mayor actividad. Las tasas mximas
de fijacin se sitan entre los estados reproductivos R5-R6 con valores promedio de 3 y
mximos de 5 Kg. de N fijado.ha-1.da-1. Luego de esta etapa el proceso cae en forma abrupta.
1.4.2 Seales moleculares de la fijacin biolgica del N y su especificidad. Las seales
que intercambian los rizobio con la planta y su especificidad, son ampliamente descritos por
Schultze y Kondorosi (1998), Spaink (2000) y
Estas
asociaciones
nodABC (genes nod comunes) y de otros genes nod esenciales. Para que se produzca la
infeccin, la protena NodD tiene que ser activada y para esto tiene que interaccionar con el
flavonoide especfico. Por esto, los factores NodD son determinantes de la especificidad de
husped. Se han podido establecer que los genes relacionados con la fijacin del nitrgeno a
nivel de la bacteria, presentan cambios mutacionales o transferencia gentica horizontal dando
origen a una gran diversidad intraespecfica. Son varios los genes nod que participan en la
simbiosis y su funcin es descrita por Hartwig et al. (1990), Pandya et al. (1998), Roche et al.
(1996) y Tsvetkova et al. (2006):
Genes nod comunes: (nodABC), son genes conservados en todos los rizobios y pueden
intercambiarse entre especies y gneros. Su ausencia impide el proceso de infeccin. Aunque
en estudios recientes han encontrado polimorfismos en nodC con series allicas que generan
variabilidad intraespecfica.
Genes nod especficos: (nodFE, nodH, nodM, nodPQ), no necesariamente presentes en todos
los rizobios. Son los responsables de la especificidad del husped. Los genes nodEF, nodM y
nodPQ modifican el factor Nod hacindolo especfico. Las mutaciones en ellos alteran o
amplan el rango de husped puesto que hay genes nod especficos que amplan el rango de
especificidad y otros que lo reducen impidiendo que un rizobio determinado infecte una planta
dada.
Genes responsables de la sntesis de polisacridos: Los productos de estos genes son
importantes para formar los canales de infeccin. Estos constituyentes de la superficie de los
rizobios como los polisacridos, son necesarios para la formacin de ndulos simbiticamente
efectivos. Los polisacridos que han sido estudiados en relacin a su funcin en simbiosis son:
exopolisacridos (EPSs), lipopolisacridos (LPSs), antgenos de tipo K (KPSs) y glucanos
cclicos (Spaink, 2000), que le dan especificidad a una cepa determinada. Los genes de
lipopolisacridos (LPS), son hipervariables.
Un aspecto importante de considerar sobre los factores Nod, fue reportado por Giraud et al.
(2007), quienes
descubrieron
Los costos econmicos y ambientales del uso intensivo de qumicos fertilizantes nitrogenados
en la agricultura son una preocupacin mundial. Una de las mayores amenazas para la calidad
del agua es la escorrenta de N en los sistemas agrcolas donde se utilizan fertilizantes ricos en
N mineral. Por lo anterior, es necesario promover el uso de mejores cultivares y de prcticas
de manejo conservacionistas, donde se haga uso eficiente de las fuentes biolgicas de fijacin
de nitrgeno para la produccin de los cultivos (Bohlool et al., 1992)
La fijacin biolgica del nitrgeno (FBN), es muy importante tanto desde una perspectiva
econmica como ambiental. La energa total utilizada para el cultivo de la soya cuando el
inculo se aplica como oferta de nitrgeno necesario para la cosecha es menor que la
necesaria para mantener el cultivo cuando los fertilizantes qumicos son la fuente de nitrgeno.
Adems, la contribucin del inculo a la energa total utilizada en el sistema (alrededor de
0,23%) es muy baja, en comparacin con la contribucin dada por el nitrgeno de los
fertilizantes empleados (12,9%). La fraccin del total de la energa no renovable necesaria para
mantener el cultivo es mayor con el uso de fertilizantes que utilizando el inculo (Panzieri et al.,
2000).
La inoculacin de soya genera un gran impacto sobre el rendimiento de grano, en suelos
donde se cultiva por primera vez. El rendimiento se puede optimizar mediante el empleo de
tecnologa ajustada de fecha se siembra, espaciamiento, control de malezas y plagas y
fertilizacin fosforada, con incrementos experimentales de 700 a 1.000 Kg. ha-1 sobre el testigo
no inoculado, cuando se utilizan inoculantes con base en cepas eficientes de B. japonicum.
Este valor disminuir si el cultivo est sometido a estrs hdrico, durante el perodo de
formacin de ndulos y/o durante el llenado de granos (Gonzlez, 2002; Gonzlez et al., 1998).
En trabajos realizados por Pietrarelli et al. (2008), se encontr que una condicin de mayor
disponibilidad hdrica determin un promedio de peso seco de ndulos totales por planta 2,2
veces superior a los registrados en la campaa con estrs hdrico.
Algunos estudios indican que en total, la FBN puede suplir hasta 300 Kg. ha-1 de N en todo el
ciclo del cultivo (Keyser y Li, 1992; Hungria y Campo, 2005). La evolucin de la fijacin
biolgica del nitrgeno (FBN), est relacionada con la tasa de acumulacin de carbono (C)
producto de la fotosntesis, por lo tanto, las limitaciones nutricionales que afecten el
crecimiento del cultivo afectarn la tasa de acumulacin de N. Por otro lado, numerosos
nutrientes intervienen directamente en el proceso de fijacin, como el magnesio Mg, Mo, Fe y
Co (Purcell, 1999).
1.4.5. Factores que afectan la asociacin simbitica. Las altas temperaturas, la sequa y la
acidez del suelo limitan la formacin y la funcin de ndulos en las races de leguminosas. En la
mayora de los casos, el microsimbionte es el ms afectado. Las plantas que crecen sobre N
mineral por lo general son menos sensibles a estos estreses. Ha sido posible incrementar los
rendimientos de grano en soya con la seleccin y utilizacin de cepas tolerantes al estrs.
Adicionalmente, las prcticas adecuadas de manejo del suelo como la no labranza, son otros
enfoques que pueden reducir la temperatura del suelo y conservar la humedad, lo que
aumentar la fijacin de N (Hungria y Vargas, 2000).
Las seales moleculares entre planta-rizobio se pueden ver afectadas a 39oC donde la
liberacin de los inductores de los genes nod de la soya, decrece (Hungria, 1995). La alta
temperatura inhibe la formacin de pelos de la raz, reduciendo el nmero de sitios para la
nodulacin y acelera su senescencia (Hungria y Franco, 1993). Las temperaturas bajas
retardan el desarrollo de la planta, la formacin de ndulos y consecuentemente, disminuyen
las tasas de fijacin de N. En soya, la nodulacin es inhibida por temperaturas superiores a 42 a
45C, aunque la fijacin de N en soya muestra ser ms tolerante al calor que en otras
leguminosas (Fernndez, 2003)
En el trpico, la disponibilidad de agua est determinada principalmente por la capacidad de
retencin de humedad del sustrato. Por lo general, a potenciales de agua -0,5 a -1,5 MPa, el
crecimiento y la fijacin de nitrgeno se ven afectadas en las zonas tropicales (Hungria y
Vargas, 2000). El proceso de FBN es altamente sensible al estrs hdrico (Sinclair et al., 2000),
comprometindose cuando la humedad del suelo es menor al 50% del agua til (Serraj et al.,
1999). El estrs hdrico afecta la supervivencia de rizobios, el crecimiento y estructura de la
poblacin en el suelo, la formacin y la longevidad de los ndulos, la sntesis de
leghemoglobina y la funcin de los ndulos. La sequa implica menos agua disponible para el
transporte de compuestos nitrogenados procedentes de ndulos. La fijacin de N con estrs
hdrico es limitada por la disponibilidad de fotosintatos (Hungria y Vargas, 2000). Frente a dficit
hdricos crecientes, la primera respuesta de la planta a nivel del sistema nodular, es inactivar la
nitrogenasa. Si la humedad edfica se restablece, tambin lo har la FBN, pero si esta se
profundiza, se producir la inactivacin definitiva de la enzima, seguida de muerte y
desprendimiento nodular. La planta crecer, a partir de entonces, a expensas del N edfico, sin
aporte adicional de N atmosfrico (Racca, R. y Pollino, D. 2005).
El exceso de agua o anegamiento, tambin afecta tanto al rizobio como a la planta, al generar
falta de O2 en el suelo. El encharcamiento llevar a la muerte de las bacterias del suelo, y a un
detrimento de las condiciones de crecimiento para la planta, que impedir la nodulacin y la
fijacin biolgica, hasta la eliminacin de los ndulos ya establecidos. El efecto de la inundacin
es ms severo en las plantas que crecen a expensas del N proveniente de la fijacin biolgica
que en las que usan el N edfico, indicando que la fijacin de N es ms sensible que la planta
misma (Fernndez, 2003).
La acidez del suelo afecta todos los aspectos de la simbiosis, desde la supervivencia y
multiplicacin de los rizobios en el suelo, la infeccin y nodulacin hasta la fijacin de N
(Fernndez, 2003). El pH bajo del suelo es a menudo asociado con toxicidad de aluminio (Al) y
de manganeso (Mn), y reducido suministro de calcio (Ca) y fsforo (P). Estas condiciones
afectan el crecimiento de rizobios. La liberacin de inductores de los genes nod de las races de
la soya fue menor a pH 4,5 que a pH 5,8 (Hungria y Stacey, 1997). Algunas especies de
rizobios pueden tolerar mejor la acidez que otros, y la tolerancia puede variar entre las cepas
dentro de una especie (Hungria et al., 1997b).
Olivera et al. (2004), sostienen que las leguminosas fijadoras de nitrgeno tienen mayores
requerimientos de fsforo que las no fijadoras, de hecho se atribuye a la aplicacin de fsforo el
incremento de nitrgeno total condicionado a la alta actividad de fijacin simbitica de
nitrgeno. Por lo que una deficiencia de P afectar negativamente el proceso simbitico.
Otro factor que afecta la nodulacin es la presencia de cepas naturalizadas que son por lo
general, ms competitivas y menos efectivas. Por ello, la inoculacin en cada ciclo de cultivo
favorece el desarrollo de ndulos implantados en la raz primaria cercanos al cuello de la planta,
los cuales tienen una actividad nitrogenasa especfica varias veces superior a los de las races
Las cepas
P=G+A+GA, modelo que puede ser descrito estadsticamente, as (Yan y Kang, 2002):
Pij = + Gi + Aj + (GA)ij
Donde,
Pij = Cada uno de las observaciones de la variable respuesta o fenotipo
comportamiento de los cultivos (Butron et al., 2004; Crossa et al., 2002; Dehghani et al., 2006;
Yan y Hunt, 2002; Kaya et al., 2006; Samonte et al., 2005; Yan y Tinker, 2005; Yan et al.,
2000; Yan, 2002; Yan y Kang, 2003; Yan y Ma, 2006).
Modelo AMMI (Additive main effects and multiplicative interaction model). Es el modelo
estadstico en el que los efectos principales de genotipos y de ambientes, considerados
trminos lineales, se explican mediante un anlisis de varianza convencional; el componente no
aditivo se atribuye a la interaccin genotipo-ambiente (GA) y se analiza mediante la tcnica de
componentes principales (CP). Se puede emplear el programa SAS para obtener el modelo
AMMI segn al algoritmo computacional creado por Vargas y Crossa (2000), donde los dos
primeros componentes principales explican en el biplot gran parte de la variabilidad generada
por la interaccin genotipo-ambiente.
Modelo SREG (Sites Regression Model). Este modelo permite obtener un biplot GGA, donde,
la varianza total es explicada ampliamente por el efecto ambiental (A) y usualmente una menor
participacin del genotipo (G) y de la interaccin (Yan, 2002; Ibaez et al., 2006). Este modelo
representa simultneamente la variabilidad de genotipos y ambientes, basada en el anlisis de
componentes principales (Yan et al., 2000), donde incluye G+GA en el termino bilineal para
producir el biplot (Crossa y Cornelius, 1997; Burgueo et al., 2002, Crossa, et al., 2002).
El biplot GGA se refiere al afecto principal genotpico (G) de la interaccin genotipo-ambiente
(GA), las cuales son las dos fuentes de variacin del modelo SREG. Aunque la medida del
rendimiento es el resultado combinado de los efectos del genotipo (G), ambiente (A) e
interaccin genotipo x ambiente (GA), solo G y GA son relevantes para la evaluacin de
cultivares. Con esta premisa Yan et al. (2000) combin los efectos de G y GA para crear la
conocida GGA, siguiendo la propuesta de Gabriel (1971) sobre biplot, como una herramienta
grfica para presentar resultados del anlisis de componentes principales. Es una tcnica de
anlisis multivariado que despliega grficamente los datos en doble va que permite la
visualizacin de la interrelacin entre ambientes (A), genotipos (G) e interrelacin entre
genotipos y ambientes (GA).
en la comprensin (Yan et al., 2000). Por lo general, el biplot GGA del modelo SREG se
construye a partir de los componentes principales CP1 y CP2 derivados de la descomposicin
del valor singular (SVD).
correlacionado con el efecto principal de genotipo, este representa la proporcin del rendimiento
que se debe slo al genotipo (G), mientras que, el CP2 representa la parte del rendimiento
debida a la interaccin genotipo-ambiente (GA). La asociacin entre ambientes est dada por el
ngulo que se forma entre sus vectores en el biplot GGA, de tal manera, que en los ambientes
cercanos con asociacin positiva se conformar un ngulo agudo entre sus vectores, mientras
que, un ngulo recto entre vectores determinar ausencia de asociacin entre ambientes y el
ngulo obtuso, una asociacin negativa u opuesta (Yan et al., 2000). Cuando los CP1 tienen
una correlacin casi perfecta con el efecto principal genotpico, los cultivares ideales podran
tener un alto CP1 o alta habilidad productiva y un bajo coeficiente CP2 o alta estabilidad. As, el
ambiente ideal de evaluacin podra tener un alto CP1 y ser ms discriminante de los genotipos
en trminos
que la correlacin de un carcter particular con otros caracteres es de gran importancia para la
seleccin indirecta de genotipos deseables.
Con frecuencia la informacin es analizada estableciendo mltiples relaciones entre variables a
travs de matrices de correlacin. En estas matrices se incluyen todas las relaciones posibles
sin establecer el sentido de la asociacin. As, los coeficientes de correlacin a pesar de tener
utilidad en la cuantificacin de la magnitud y direccin de los efectos de factores en la
determinacin de caracteres complejos, estos coeficientes no indican la exacta importancia que
tienen los efectos directos e indirectos de esos caracteres sobre el rendimiento (Abbott et al.,
2007). Un coeficiente de correlacin alto o bajo entre dos variables se puede deber al efecto de
una tercera variable o grupo de variables. Por ello, segn Seker y Serin (2004), la correlacin
entre dos caractersticas no predice el xito de la seleccin.
Como alternativa para superar este inconveniente, Wright (1921) desarroll un mtodo analtico
denominado anlisis de coeficientes de sendero (path coefficient analysis) o coeficientes de
trayectoria, creado como alternativa estadstica para investigar la ramificacin de varios
modelos causales en poblaciones genticas y hoy da, utilizado en diferentes reas del
conocimiento. El mtodo asume que las variables se asocian entre si mediante relaciones
lineales (Roehrig, 1996) y que sus coeficientes son estandarizados por cuanto son estimados a
partir de los coeficientes de correlacin. El anlisis permite descomponer las correlaciones
entre dos variables en efectos directos e indirectos (Mitchell, 1992) y facilita identificar las
posibles explicaciones causales de las correlaciones observadas entre una variable respuesta
(endgena o dependiente) y una serie de variables predictoras (exgenas o independientes)
(http: luna.cas.usf.edu/// regression/Pathan.html). Este mtodo ha sido de gran utilidad para el
mejoramiento de plantas donde se hace descomposicin de la correlacin en efectos directos e
indirectos identificando cadenas de causalidad. Diversas investigaciones (Whittaker et al., 2009;
Oz et al., 2009; Gaikwad et al., 2007), han demostrado la importancia de este tipo de anlisis
para la formulacin de estrategias apropiadas de seleccin. Este mtodo difiere de los modelos
de regresin mltiple en que permite establecer las relaciones causales entre las variables
(Carrascal, 2004; Mitchell, 1992).
En el diagrama simple del anlisis de sendero (Figura 3), se observan los coeficientes de
Wright o coeficientes de causalidad, que corresponden a la relacin entre dos variables. Se
denotan como pij, el subndice i representa la variable exgena y el subndice j representa a la
variable endgena. Para las variables exgenas el rxy representa el coeficiente de correlacin
simple entre ellas.
Figura 3. Diagrama simple del anlisis de
X1
p13
Z3
rxy
Y2
p23
Residual
Los coeficientes de Wright o coeficientes de sendero, han sido utilizado en soya para diversos
tipos de estudios. Para identificar criterios de seleccin indirecta en siembras tardas (Board et
al., 1997), efectos del ambiente sobre el rendimiento (Board et al., 1999), componentes de
rendimiento en diferentes tipos de siembra en soya (Pandey y Torrie, 1973), efecto de las
densidades de poblacin sobre el rendimiento (Ball et al., 2001), componentes de rendimiento
en hbridos de soya (Taware et al., 1997), interferencia de malezas sobre el crecimiento de la
planta (Jordan, 1992),
y de
cultivo, dominantes que difieren de un regin a otra y que pueden ser ms marcadas entre
condiciones tropicales y templadas. Por lo anterior, cualquier programa de mejoramiento que
pretenda desarrollar con xito nuevas alternativas varietales para un ecosistema especfico,
deber tener en consideracin el conocimiento previo de la interrelacin de caracteres
genotpicos en su zona de influencia.
1.6 Propsitos de la investigacin en soya.
Para abastecer parte de la demanda interna de grano de soya y competir con los mercados
internacionales, es necesario que Colombia afiance los procesos de investigacin y desarrollo
tecnolgico en busca de una produccin competitiva. Una alternativa viable, est en la
generacin de nuevas variedades con adaptacin especfica, de alto potencial de rendimiento,
de alta FBN y eficientes en la utilizacin de los recursos de produccin, enmarcadas en el
contexto conservacionista de los agroecosistemas.
Aunque los avances han sido significativos en la generacin metodos de seleccin rpida de
genotipos
tolerantes
aluminio,
generacin
de
nuevas
variedades
estables
Orinoquia colombiana.
El propsito de la presente investigacin fue determinar e interpretar el fenmeno de la
interaccin entre variedades, cepas de B. japonicum y ambientes, por medio del anlisis de los
efectos principales adictivos y multiplicativos de los modelos AMMI y SREG, para definir los
ambientes de mayor discriminacin y la mejor combinacin de variedades y cepas para
maximizar el rendimiento de grano, estimar las asociaciones directas e indirectas entre
caractersticas agronmicas, de fijacin biolgica del N y los componentes de rendimiento sobre
el rendimiento de grano, en oxisoles de la Orinoquia colombiana, como estrategia para futuros
programas de mejoramiento.
2. OBJETIVOS
2.1 General:
Determinar la interaccin existente entre variedades de soya y cepas de bacterias simbiticas
de Bradyrhizobium japonicum, en oxisoles de la Orinoquia colombiana, para potenciar el efecto
FBN en la productividad del cultivo.
2.2
Especficos
3.
MATERIALES Y MTODOS
3.1 Localizacin.
Los ensayos experimentales multiambientes (EMA), de variedades de soya y cepas de B.
japonicum, se realizaron durante el 2007 y 2008, en oxisoles de la Orinoquia colombiana, en
latitudes entre 0403' norte y 7329' oeste, con altitudes entre 150 y 336 msnm. Las Localidades
o ambientes fueron identificados como: Lib: Centro de Investigaciones La Libertad, ubicado en
el municipio de Villavicencio, Tal: Estacin experimental Taluma, localizada en Pueto Lpez
en cercana al municipio de Puerto Gaitn y los ambientes de Scruz y Hor corresponden a
las fincas Santa Cruz y Hortensias del municipio de Puerto Lpez.
fue
estimada
mediante
el
software
Daily
versin
(Granfield
University,
Municipio-Dpto.
Localidades experimentales
Altitud msnm
La Libertad
Villavicencio-Meta
336
Puerto Lpez-Meta
Tal08A y Tal08B
150
Puerto Lpez-Meta
Scruz08A y Hor08B
195
0403' N 7329' W
Las Margaritas
0421' N 7210' W
La Palomera
4 15' N 72 33' W
Tabla 4. Datos climatolgicos acumulados por ciclo del cultivo y por semestre. Piedemonte
llanero y Altillanura colombiana (2007-2008).
Subregin
Piedemonte 2007
Piedemonte 2008
Altillanura 2008
Sem.
Precipitacin
ETP
Temperatura
Humedad
Brillo
acumulada
acumulada
media
relativa
Solar
(mm)
(mm)
( C)
(%)
(h.da-1)
1370
325
25,3
87
4,5
1122
355
25,9
86
5,5
1566
293
25,0
88
4,3
976
358
25,7
84
6,2
1238
332
26,5
80
4,4
634
395
26,0
85
6,5
para Scruz08A; Maz, para Hor08B y sabana nativa, para Tal08A y Tal08B. En ninguno de los
lotes hubo siembras previas de soya en su historial.
Tabla 5. Caractersticas fsico-qumicas de los suelos de las localidades en estudio.
Variable/ loc
Lib07a
Na (cmol.Kg-1)
CICE (cmol.Kg-1)
Sat Al (%)
Sat bases (%)
S (%)
Fe (ppm)
B (ppm)
Cu (ppm)
Mn (ppm)
Zn (ppm)
4,6
2,4
8,5
1,6
1,3
0,8
0,2
0,2
0,1
2,9
44,9
44,3
10,0
23,5
0,1
0,3
5,0
0,2
Arena
Limo
Arcilla
Textura
Densidad Ap. 10 cm
Densidad Ap. 20 cm
58
22
21
FArA
1,50
1,61
PH
M.O (%)
P(ppm)
Acidez interc. (cmol.Kg-1)
-1
Al intercamb. (cmol.Kg )
-1
Ca (cmol.Kg )
Mg (cmol.Kg-1)
K (cmol.Kg-1)
Tal08b
Scruz08a
Hor08b
5,0
1,5
1,8
0,7
0,5
0,4
0,2
0,0
0,0
1,3
41,6
44,5
2,0
33,3
0,2
0,2
1,9
0,3
4,5
3,1
4,0
1,9
1,4
0,7
0,3
0,1
0,1
3,2
44,9
37,1
1,0
46,5
0,2
0,8
2,3
0,8
4,7
2,2
2,0
2,0
1,4
0,3
0,1
0,1
0,1
2,5
57,4
18,0
1,0
10,0
0,1
0,2
0,5
0,3
59
8
33
FA
1,62
1,77
53
30
17
FA
1,43
1,40
34
32
35
Far
1,40
1,50
La siembra de los experimentos se realiz bajo el sistema de labranza reducida, con un pase de
rastra e incorporacin de la cal dolomita un mes antes de la siembra, en dosis basada en la
saturacin de bases segn la frmula de Raij et al. (1985) y ajustada al 50%, un pase de cincel
rgido y un pase pulidor.
Cal dolomita t.ha-1 = [(SB2 SB1) x CICE] / PRNT
-1
urea. Este nivel de N fue establecido segn estudios de Snchez et al. (1990), realizados en
oxisoles de la Orinoquia colombiana
La inoculacin de la soya se realiz aplicando 10 g. Kg-1 de semilla, de una suspensin de
bacterias >108 ucf/ml en mezcla con agua azucarada al 10%. Para el tratamiento de cepas (J01+J-96), se utiliz 5 g. Kg-1 de semilla, de cada una. Las variedades de soya fueron inoculadas
antes de la siembra y en condiciones de sombra, con las cepas de
referencia en forma
individual. Se dejaron testigos de cada variedad con la aplicacin de nitrgeno mineral (150
Kg.ha-1). La dosis de N (urea), fue fraccionada en dos tiempos (10 y 25 das despus de
emergencia) y aplicada en franja paralelas al surco de plantas. La dosis de potasio fue
igualmente fraccionada, con una primera aplicacin al fondo del surco al momento de la
siembra y una segunda aplicacin a los 10 das con la primera dosis de N. La fertilizacin con
elementos mayores y menores se realiz al momento de la siembra, dosificada de acuerdo con
el anlisis de suelo. Las fuentes de fertilizantes fueron: superfosfato triple, cloruro de potasio,
borocinco y sulcamag.
Los efectos y la especificidad de la accin simple o mezcla de cepas de B. japonicum sobre el
comportamiento agronmico y contenido proteico de la soya, se realiz a travs de la medicin
de variables agronmicas y de fijacin biolgica del N.
Se tomaron
Estado de evaluacin
Perodo comprendido entre la emergencia y la aparicin de una
flor en el 50% de las plantas de la parcela.
Perodo comprendido entre la emergencia y cuando el 95% de
las vainas estn secas.
Altura en cm. planta -1 tomada desde la superficie del suelo hasta
el ltimo nudo.
Nmero de nudos desde la hoja simple al nudo superior donde la
hoja trifoliada este completamente abierta.
Evaluacin a los 65 das despus de la siembra: Bacteriosis,
Pstula bacterial, Cercospora, Peronospora, roya y virosis. La
evaluacin de severidad se midi en porcetaje del rea afectada.
Nmero de vainas totales por planta, como promedio de 5
plantas a cosecha.
Nmero de vainas sin granos, promedio de 5 plantas, a cosecha.
Nmero de compartimentos de la vaina sin grano, promedio de 5
plantas, a cosecha.
Nmero de vainas con 1, 2, 3 y 4 granos, promedio de 5
plantas, a cosecha.
Peso (g) de 100 granos por parcela, a cosecha.
Anlisis de los contenidos de protena (%) del grano, a cosecha.
Peso de grano por parcela y transformado a t.ha-1 con 13% de
humedad.
separadas en parte rea y races, y las races fueron lavadas con agua de grifo sobre zarandas
para extraer los ndulos. Se cuantific el nmero de ndulos de las races (NN) por planta y
luego, secados en horno a 60C durante 24 h, para registrar su peso seco en gramos (PSN). En
cada cosecha, la parte area y races se secaron en estufa a 70 C durante 48 h para obtener
los correspondientes pesos secos areos (PST) y de raz (PSR).
VARIEDADES
V1: Soyica P-34
CEPAS
C1: J-01
V2:
V3:
V4:
V5:
V6:
C2:
C3:
C4:
N :
Orinoquia 3
C. La Libertad 4
C. Taluma 5
C. Superior 6
C. Sabana 7
J-96
J-98
J-96 + J-01
Test. 150 Kg.ha-1
Las variables
NN,
componentes de rendimiento (V1, V2, V3, V4, VV, AE, VT), a raz cuadrada de (x+1). Para el
anlisis del rendimiento de grano, se consider como covariable el nmero de plantas por
parcela. Las medias de estas variables fueron ajustadas utilizando el procedimiento LSMEANS.
Los coerfcientes de correlaciones de Pearson fueron estimados con el procedimiento PROC
CORR del paquete estadstico SAS.
Los anlisis de interaccin variedad-cepa-ambiente para las variables agronmicas con
significancia estadstica (P<0.05), se realizaron a travs del modelo SREG (Burgueo et al.
2002), desagregando los componentes principales de mayor significancia y graficando sus
tendencias a travs del biplot de interaccin GGA.
Los componentes de varianza, las correlaciones genticas y fenotpicas de variedades de soya
inoculadas con diferentes cepas de rizobio en suelos de la altillanura colombiana, fueron
obtenidos a travs de los anlisis simples y combinados de la informacin registrada en las
pruebas de campo por localidad y las pruebas de laboratorio. Se aplic el anlisis de sendero
(Wright, 1921), para particionar los coeficientes de correlacin en efectos directos e indirectos
de variables agronmicas, de fijacin biolgica, rendimiento de grano y sus compenentes.
3.6.1. Anlisis de correlacin y de sendero (Path analysis). Las correlaciones simples (r)
fueron obtenidas para todas las posibles combinaciones de las caractersticas relacionadas con
el rendimiento de grano y su significancia fue establecida a P<0.05 y P<0.01, como lo describen
Steel y Torrie (1989). La correlaciones genticas fueron estimadas con el programa Genes
(Cruz, 2007). Este programa, aplica las formulas de correlacin descritas por Falconer (1986).
- Correlacin fenotpica
- Correlacin gentica
- Correlacin ambiental
Donde,
(x)
(y)
son la desviaciones
p15
p25
X5
p35
p45
X1
Los
r12
X2
r14
r24
r34
Residual
Figura 5.
de
sendero
fueron
p 45,
r13
r23
X3
coeficientes
X4
VARIEDAD
CEPA
Vc
Soyica P-34
J-01
C. Taluma-5
J-98
16
Orinoquia 3
J-01
C. Superior-6
J-98
17
C. Libertad-4
J-01
C. Sabana-7
J-98
18
C. Taluma-5
J-01
C. Sabana-7
J-01 + J-96
19
C. Superior-6
J-01
C. Superior-6
J-01 + J-96
20
C. Sabana-7
J-01
C. Taluma-5
J-01 + J-96
21
C. Sabana-7
J-96
C. Libertad-4
J-01 + J-96
22
C. Superior-6
J-96
Orinoquia 3
J-01 + J-96
23
C. Taluma-5
J-96
Soyica P-34
J-01 + J-96
24
Soyica P-34
-1
25
-1
26
-1
27
-1
28
-1
29
-1
30
C. Libertad-4
Orinoquia 3
Soyica P-34
Soyica P-34
Orinoquia 3
C. Libertad-4
J-96
J-96
J-96
J-98
J-98
J-98
10
11
12
13
14
15
VARIEDAD
Orinoquia 3
C. Libertad-4
C. Taluma-5
C. Superior-6
C. Sabana-7
CEPA
150 Kg.ha N
150 Kg.ha N
150 Kg.ha N
150 Kg.ha N
150 Kg.ha N
150 Kg.ha N
Vc
Los estudios de interaccin genotipo-ambiente fueron realizados mediante los modelos AMMI
(Additive Main effects and Multiplicative Interaction model) y SREG (Sites Regression model).
Los anlisis de varianza de estos modelos permitieron determinar la contribucin relativa del G,
A e interaccin GA de la varianza total.
Se emple el programa SAS para obtener el modelo AMMI segn al algoritmo computacional
creado por Vargas y Crossa (2000) y del modelo SREG (Burgueo et al. 2002), los cuales
permitieron generar cada biplot donde GA hace referencia al grfico de la interaccin
genotipo-ambiente obtenido del modelo AMMI y el biplot GGA del modelo SREG. Con el
modelo SREG se realizaron los anlisis de estabilidad fenotpica e identificacin de ambientes
ms discriminatorios y representativos.
3.6.3. Modelos estadsticos. Los anlisis estadsticos de las variables agronmicas de las
variedades de soya, como las variables de rizobio, fueron analizadas con el siguiente modelo
estadstico para franjas divididas (Steel y Torrie, 1989).
variedad
Yge = + e + N n gn en + ge
Yge = Rendimiento promedio de un genotipo g en el ambiente e.
= Media general.
e = Desviaciones de las medias de ambientes.
N = Nmero de CP retenidos en el modelo.
n = Valor singular para el CP.
Yge = + g + e + N n gn en + ge
Yge = Rendimiento promedio de un genotipo g en el ambiente e.
= Media general.
g = Desviaciones de las medias de genotipos.
e = Desviaciones de las medias de ambientes.
N = Nmero de CP retenidos en el modelo.
n = Valor singular para el CP.
gn = Valores de los vectores de los genotipos para cada CP.
en = Valores de los vectores de los ambientes para cada CP.
ge = Residual.
El paquete estadstico y grfico empleado para procesar ambos anlisis multiplicativos fue el
propuesto por Burgueo et al., (2002), realizado en SAS y denominado Graphing GE and GGE
Biplot.
El nmero de trminos AMMI y SREG necesarios para una descripcin adecuada de la
interaccin genotipo por ambiente, se evalu con el mtodo de Gollob (1968). En la Tabla 7 se
observa la particin de la varianza de la interaccin en sus componentes. El cuadrado medio
de cada trmino se obtuvo a partir de los eigenvalues o valores propios, la suma de cuadrados
de estos valores caractersticos se dividi por los grados de libertad, calculados segn el
procedimiento propuesto por Gollob (1968): gl= g + a - 1- 2n (g, nmero de genotipos; a:
G.L Gollob
Suma de cuadrados
Genotipos (G)
G -1
SC (G)
Ambientes (A)
A1
SC (A)
(G-1)(A-1)
SC (GxA)
CP1
G+A-1-(2x1)
21
CP2
G+A-1-(2x2)
22
CPn
G+A-1-(2xn)
2n
GA(n-1)
SC(error medio)
GAr-1
SC (total)
Interaccin (GxA)
Error medio
Total
coeficientes del CP2 (eje y). Los ambientes fueron presentados como vectores que se originan
en la coordenada (0,0) del biplot y se extienden hasta el marcador correspondiente.
Con el anlisis de la informacin, se identificaron los mejores tratamientos en trminos de
produccin de grano y de estabilidad fenotpica, para generar la respectiva recomendacin
tcnica. Los modelos AMMI y SREG no permiten establecer las diferentas estadsticas entre
genotipos o tratamientos en los biplot, razn por la cual, mediante anlisis de contrastes (Steel
y Torrie, 1989), utilizando el comando CONTRAST del SAS, se seleccionaron las mejores
combinaciones de la interaccin entre variedades y cepas de B. japonicum.
Extraccin de ADN: El ADN total se aisl a partir de hojas jvenes de las variedades de soya
siguiendo el protocolo del Laboratorio de Gentica Molecular Vegetal de Corpoica. En ste
procedimiento se utiliz 500 mg de hojas maceradas con nitrgeno lquido. Se adicion buffer
de extraccin (0,15 M de Tris-HCl pH 8,0; 0,015 M de EDTA pH 8,0; 1 M de NaCl; 5% de mercaptoetanol; 1% de CTAB y 1% de PVP 360000) y se incub a 65C por 30 minutos. Luego,
se utilizaron 400 L de fenol equilibrado y 400 L de cloroformo: isoamlico en un mismo paso,
para eliminar contaminantes como protenas y fenoles. Se realizaron dos precipitaciones, la
primera con 1 mL de 2-propanol seguida por una limpieza con 400 L de cloroformo: isoamlico
del pellet diluido en agua destilada; la segunda precipitacin se hizo con 1/10 de volumen de
acetato de sodio 3 M, pH 5.2 fro (-20 C) y 6/10 de volumen de isopropanol fro. Por ltimo, se
realizaron dos lavados, cada uno con 1 mL de etanol al 70%. El ADN obtenido se diluy en 100
L de solucin TE/RNAsa (10 mM de Tris HCl, pH 8.0; 1 mM EDTA, pH 8.0 y 25 mg/mL de
RNAsa).
La calidad del ADN se verific en geles de agarosa al 0,8% en cmara de electroforesis
horizontal Horizon 58 (Life Technologies), y teidos con bromuro de etidio a 0,1 g/ml, se corri
a 100V por 1 hora usando la fuente de poder Bio-Rad 1000/500. La concentracin del ADN se
determin por cuantificacin a partir de un marcador de masa (High DNA Mass Ladder,
Invitrogen), usando el software Gene-Tool. Este programa compara el patrn de bandas
generado y la banda obtenida por cada muestra y da un estimativo de la concentracin obtenida
Sentido
Satt 131
Forward
Secuencia 5 - 3
GGC TTG TGG TAA GGG GAA GAC TTC
Forward
GCG AAT ACA TAA AAC TCA AAT TCA AAT CAT A
Reverse GCG TTC TAT AAA TTT CAT TCA TAG TTT CAA T
Sat_039
Forward
Forward
Forward
Forward
se coloc el vidrio con el gel polimerizado en la solucin fijadora (1% de cido actico glacial v/v
+ 10% etanol 95%) durante 10 minutos; segundo, se hizo un lavado con agua desionizada
durante 1 min; tercero, se introdujo en solucin oxidante (1,5% de cido ntrico 65% v/v) durante
3 min; cuarto, lavado con agua desionizada durante 1 min; quinto, se sumergi el vidrio en
solucin de tincin (1% de nitrato de plata p/v + 0,12% de formaldehido 37% v/v) por 20 min;
sexto, se realiz un lavado con agua desionizada durante 10 segundos; sptimo, se pas a la
solucin de revelado (3% carbonato de sodio p/v + 0,054% de formaldehido 37% v/v) con
agitacin constante hasta que las bandas se hicieron visibles; octavo, una vez reveladas las
bandas de ADN, se llev el gel a una solucin de parada (1% de cido actico glacial v/v)
durante 3 min; noveno, se realiz un lavado con agua desionizada por 1 min. El gel se dej
secar durante la noche. Las imgenes se registraron con un escner HP ScanJet XPA en
formato JPEG. El peso de las bandas se determin utilizando el programa Gene Tools y como
patrn el marcador 25pb (Invitrogen).
3.7.2
rep-PCR. Las reacciones PCR con secuencias repetidas altamente conservadas se realizaron
con los primers o cebadores ERIC1R (3'- CACTTAGGGGTCCTCGAATGTA-5'), ERIC2 (5'AAGTAAGTGACTGGGGTGAGCG-3'),
BOXA1
(CTACGGCAAGGCG
ACGCTGACG),
descritos por Versalovic, et al. (2004). Las reacciones fueron realizadas en un volumen final de
25 l, conteniendo 50 ng de ADN, 1U de Taq ADN polimerasa,
dNTPs y 0,1 M de cada cebador. El proceso de amplificacin con los cebadors BOX y ERIC
consisti de 1 ciclo inicial de denaturacin a 95oC por 6 min, 35 ciclos de denaturacin con la
siguiente secuencia: 94C for 1 min, 50C for 1 min y 65C for 8 min., y un 1 ciclo final de
extensin de 65 oC por 16 min. Esta PCR fue corrida en un termociclador PT100 MJ-Research y
los fragmentos fueron visualizados por electroforesis sobre un gel de agarosa al 1,5% teido
con bromuro de etidio al 0,1 g/ml y corrido a 100 V por 3 h. Posteriormente, el tamao de las
bandas que gener cada cebador fue establecido con base en el marcador de peso molecular 1
kb Plus DNA ladder (Gibco BRL). Para cada cebador se hicieron ms de tres amplificaciones y
corrido para verificar su reproducibilidad y la consistencia del patrn de bandas.
Anlisis de la informacin: Las bandas obtenidas por cada cebador fueron transformadas en
una matriz binaria basada en la presencia (1) o ausencia (0) de cada fragmento amplificado.
Con el paquete estadstico NTSY-pc (Rohlf, 1997), se determin el ndice de similitud de Dice
(1945) y el dendograma final de variedades de soya se construy utilizando el mtodo de
agrupamiento UPGMA (Unweighted Pair Group Method with Arithmetic Mean). El contenido de
informacin polimrfica para SSR fue determinada por la frmula CIP=1-p2, donde p es la
frecuencia de los alelos en estado homocigoto.
4. RESULTADOS Y DISCUSIN
4.1
distintas variedades de soya. Investigaciones realizadas por Ellis et al. (1992), demostraron
que la sensibilidad al fotoperodo en soya no es uniforme durante el perodo comprendido entre
la siembra y la floracin.
Las variedades de soya difieren en los tiempos de cada fase de induccin de la floracin y su
duracin puede estar afectada por la temperatura. La interaccin compleja del fotoperodo con
la temperatura controlan los procesos de floracin, razn por la cual, es difcil mediante
modelos predecir la floracin y la madurez de este cultivo (Acock y Acock, 1995; Zhang et al.,
2001). Sin embargo, las variedades colombianas en estudio, presentan alta adaptacin al
trpico con la presencia de un perodo juvenil prolongado que evita que la soya se florezca
tempranamente y comprometa el rendimiento de grano, por lo que, la temperatura se constituye
en el factor de mayor inters. La mayor temperatura mxima, promedia y mnima se observaron
en el segundo semestre y la localidad con mayores temperaturas se dieron en la localidad
Tal08B (Tem mxima 31.7oC, media 26oC y mnima 23oC). Sin embargo, el efecto de la
temperatura y el dficit hdrico se hizo ms marcada en la madurez que en la floracin. Las
localidades con mayores temperaturas y menor disponibilidad de agua en el segundo semestre
se asociaron con mayor precocidad, hasta con 18 das de diferencia con repecto al primer
semestre, aunque la interaccin GA fue significativa.
En contraste con la floracin, la variedad ms precoz a madurez fisiolgica fue C. superior 6
con un promedio de 86,7 das, seguida por C. Libertad 4 con 89,1 das. Las variedades ms
tardas fueron Soyica P-34 y C.Taluma 5 con 93,1 y 94,7 das, respectivamente (Figura 6). Esto
significa que la variedad C. Superior 6 aunque tarda ms en florecer que Soyica P-34, se
madura precoz, al reducir su periodo de llenado de grano y por ende, produciendo granos de
menor peso.
D a s a f lo ra c i n y m a d u re z f is io l g ic a
100,0
ab
cd
90,0
bc
DF
DM
80,0
70,0
60,0
50,0
40,0
ab
ab
C. Sabana 7
C. Taluma 5
30,0
20,0
10,0
0,0
Soyica P-34
Orinoquia 3
C. Superior 6
C. Libertad 4
Variedades
Figura 6. Das a floracin (DF) y a madurez fisiolgica (DM) de seis variedades de soya, en
oxisoles de la Orinoquia colombiana. Barras con la misma letra no difieren estadsticamente
Los efectos ambientales sobre DM en las variedades, fueron marcados para las variedades
Soyica P-34 y C. Talma 5. La primera variedad present alta frecuencia de tallos verdes y
retencin foliar cuando su perodo de madurez coincidi con altos niveles de precipitacin,
situacin muy frecuente en el primer semestre del ao en el piedemonte llanero y en algunas
ocasiones, en la altillanura. Este fenmeno de retencin foliar pudo estar asociado con el
denominado sndrome de tallo verde (STV), que afecta negativa y en forma severa la cantidad
y calidad de los granos. Segn Hill et al. (2006), el STV es un trastorno de la soya que hace
que los tallos permanecen verdes, no senescentes y hmedos, aunque las vainas y semillas
pueden estar maduras o secas. Estos autores mencionan que existe variabilidad gentica entre
cultivares para este carcter, que puede ser de utilidad para desarrollar variedades con baja
sensibilidad a este sndrome.
Las variedades colombianas en estudio difieren en este carcter STV, donde Soyica P-34 es la
ms sensible, razn por la cual no se recomienda para siembras de primer semestre en la
Orinoquia colombiana (Valencia et al., 2006). Varios trabajos de investigacin en soya
(Schwenk y Nickell (1980), Hobbs et al., 2006; Villar et al., 2006; Hill et al., 2006; Basanta et al.,
2007), se han realizado con el propsito de determinar las causas del STV y aunque algunos
han establecido su origen en factores biticos como la presencia de plagas como el chinche,
algunos patgenos y factores abiticos de estrs, en general, han concluido que el STV es un
fenmeno muy complejo.
En el caso de C. Taluma 5, su mayor limitante en DM est en el segundo semestre, que ante el
dficit hdrico frecuente en la altillanura, se madura precoz. Esta variedad es la ms tarda y por
lo tanto, necesita un perodo ms prolongado de lluvias en el segundo semestre en la altillanura
con respecto a las otras variedades, para alcanzar el llenado completo del grano. En el primer
semestre en el piedemonte C. Taluma 5 presenta un ciclo ms tardo con madurez
desuniforme, por ello, solo se recomienda como forraje para corte entre 65 y 70 das despus
de emergencia (Valencia et al., 2006b).
4.1.2. Altura de planta y nmero de nudos a madurez fisiolgica. La altura de planta (ALM)
y el nmero de nudos a madurez (NM), presentaron diferencias significativas (P<0.01), entre
ambientes (A), variedades (G), cepas-N(C) e interacciones AG, AC y GC. La interaccin triple
AGC no fue estadsticamente significativa. Esto significa que, la ALM y NM es afectada por el
ambiente y que su comportamiento es variable dependiendo de la variedad y de la cepa-N
utilizada.
Las variedades que en promedio presentaron la mayor altura de planta a madurez fueron, en su
orden: C. Taluma 5, Orinoquia 3 y C. Sabana 7, con valores de 80, 74,3 y 73,3 cm., a diferencia
de Soyica P-34 y C. Superior 6 que obtuvieron un altura promedia de 59 y 56,7 cm.,
respectivamente. La menor altura fue para la variedad Corpoica Libertad 4, la cual difiri
estadsticamente de las dems con todas las cepas. Es de anotar que esta variedad es de
hbito de crecimiento determinado, que por lo general, son de menor porte que los
indeterminados. Con excepcin de C. Libertad 4 y Soyica P-34, la altura de planta fue
significativamente ms baja con la cepa J-01 (Figura 7) y la altura fue superior en todas las
variedades al aplicar 150 Kg.ha-1 de N.
100
Lib4
90
Ori3
P34
Sab7
Sup6
Tal5
80
70
60
50
40
30
20
10
0
J01
J96
J98
J96+J01
Fuente de nitrgeno
20
Lib4
Ori3
P34
Sab7
Sup6
Tal5
18
16
14
12
10
8
6
4
2
0
J01
J96
J98
J96+J01
Fuente de nitrgeno
Figura 8. Nmero de nudos por planta (NM) promedio y error estndar de 6 variedades de
soya inoculadas con diferentes fuentes de nitrgeno (cepas de B. japonicum 150 Kg.ha-1 de
N).
La variedad C. Taluma 5 obtuvo el mayor nmero de nudos asociado a su periodo vegetativo
ms prolongado. El nmero de nudos por planta fue estadsticamente superior en la localidad
Lib07A con un promedio de 16,4 pero con entrenudos cortos, en contraste con la localidad
Lib08A donde las plantas en promedio alcanzaron el menor nmero de nudos y entrenudos ms
largos. Segn FAO 1995, el nmero de nudos en el tallo principal depende de la reaccin de los
individuos al fotoperodo durante el cual crecen y del tipo de hbito de crecimiento. En general,
la altura de planta depende del nmero de nudos formados y el nmero de nudos depende de
la duracin del da o fotoperodo (Bastidas y Agudelo, 1994).
4.1.3 Nmero y peso seco de ndulos por fase fenolgica. El nmero de ndulos (NN) y su
peso seco (PSN) por planta, fueron estadsticamente diferentes (P<0.01) entre la fases
fenolgicas de floracin, formacin y llenado de vainas, y entre cepas-N, sin interaccin
significativa de fase fenolgica x cepa. Con los resultados descritos en la Tabla 8, se puede
evidenciar el crecimiento progresivo de los rizobios al incrementar significativamente el NN y
PSN de floracin a llenado de grano. Los mayores valores para NN y PSN se alcanzaron en la
fase de llenado de grano (R5-R6), en un relacin 3:1 con respecto a la fase de floracin (R1R2), mientras que con la fase de formacin de vainas (R3-R4), esta relacin fue de 1,2:1 para
NN y 1,4:1 para PSN.
Tabla 8. Nmero de ndulos (NN) y peso seco de ndulos (PSN), en tres fases fenolgicas y su
correlacin con rendimiento de grano, en oxisoles de la Orinoquia colombiana.
Nmero de
Peso seco de
r
ndulos
ndulos/planta (g)
0,50
Floracin (R1/R2)
11,21 c
0,10 c
0,54
Formacin de vainas (R3/R4)
25,70 b
0,21 b
0,44
Llenado de grano (R5/R6)
29,80 a
0,30 a
*Letras diferentes representan significancia estadstica (Tukey 5%).
r= correlacin de la variable con rendimiento de grano.
Fase fenolgica
Las
cepas
J-96
mezcla
(J-01+J-96),
presentaron
un
r
0,52
0,59
0,39
comportamiento
similar
estadsticamente diferentes (P<0.01) de J-98 y J-01, como se puede observar en las Figuras 9A
y 9B. Aunque la competitividad en estudios de B. japonicum tradicionalmente hace referencia a
la poblacin de rizobio formados de un determinado inoculante comparado contra los ndulos
formados por cepas naturalizadas (Simanungkalit et al., 1996), en el presente experimento se
pude considerar la competitividad como la diferencia estadstica entre las cepas evaluadas, por
cuanto las siembras se realizaron en lotes sin siembras previas de soya.
Las cepas de B. japonicum son competitivas por tener una mejor adaptacin a las condiciones
adversas, pero no necesariamente son ms efectivas (Fernndez, 2003). Con este criterio y
como se pude observar en la Figura 9, la cepa J-96 y la mezcla presentaron una tasa de
crecimiento mayor (1,11 a 1,22 ndulos por da, respectivamente), mientras que en J-01 y J-98
esta estuvo entre 0,67 y 0,87 ndulos por da, hasta la fase de formacin de vainas. Entre
formacin de vainas y llenado, la tasa de crecimiento de ndulos se redujo en todos los casos,
aunque de manera ms marcada con J-96 y mezcla. Investigaciones realizadas por el
laboratorio de germoplasma de microorganismos de Corpoica, corroboran esta observacin, al
establecer en medio lquido que la cepa J-96 crece a una tasa superior (aprox. dos veces), que
J-01 (informacin aun no publicada).
J01
J98
35,0
J96
Mez
30,0
25,0
20,0
15,0
10,0
5,0
0,0
Floracin
Formacn vainas
Etapa fenlogica
Llenado vainas
P e s o s e c o d e n o d u l o s / p l a n ta (g )
N m e ro d e n o d u l o s / p l a n ta
40,0
0,40
0,35
J01
J98
0,30
J96
Mez
0,25
0,20
0,15
0,10
0,05
0,00
Floracin
Formacn vainas
Llenado vainas
Etapa fenlogica
Figura 9. A. Nmero (NN) y B. peso secos de ndulos (PSN) de cepas de B. japonicum en tres
fases fenolgicas de la soya. Las letras a y b representan diferencias estadsticas promedias
entre cepas.
Al asociar la nodulacin con la FBN, estos resultados sugieren que la fijacin de N puede ser
variable a lo largo del ciclo. Segn Fernndez (2003), esta es muy baja en los estadios
vegetativos, donde son ms importantes los aportes de nitrgeno desde el suelo. No obstante,
es la etapa crtica para la formacin de los ndulos con cepas de calidad seleccionada provistas
por los inoculantes. La mayor tasa de fijacin biolgica en soya se produce a partir del
comienzo de las fases reproductivas de fin de floracin y llenado de grano, que coincide con la
mayor demanda de nitrgeno. Luego de R5, la tasa de fijacin disminuye hasta hacerse nula
antes de la madurez fisiolgica. Este fenmeno ocurre porque la planta destina gran cantidad
de hidratos de carbono (fuentes de energa para el rizobio) hacia los granos, desactivando as
la fijacin biolgica de N2. As, los ndulos mueren y se desprenden.
Con los resultados previos, y aunque los coeficientes de correlacin de NN y PSN en las tres
fases fenolgicas con rendimiento de grano fueron positivos y significativos, los anlisis de las
variables respuesta del comportamiento agronmico fueron interpretados con base en los
promedios obtenidos en la fase fenolgica de formacin de vainas por tener la mayor
correlacin con rendimiento (r=0,54 y r=0,59, con NN y PSN, respectivamente).
Estudios
realizados por Greder et al. (1987), tambin encontraron relacin positiva entre PSN y
rendimiento de grano y sugirieron el PSN como criterio de seleccin, mientras que Pietrarelli et
al. (2008) y Lzaro et al. (2004), no hallaron asociacin entre rendimientos en grano y el patrn
de nodulacin.
4.1.4 Nodulacin y produccin de biomasa, en estado de formacin de vainas. Al evaluar
el comportamiento de las cepas de B. japonicum en las diferentes variedades de soya, se
obtuvieron resultados que evidencian la especificidad cepa-variedad y la habilidad de los
ndulos para fijar nitrgeno.
El anlisis combinado de varianza de ambientes, variedades y cepas-N, permiti establecer
diferencias estadsticas a nivel de P<0.01, entre ambientes (A), variedades (G), cepas-N (C) e
interacciones AG, AC, GC y AGC para nmero ndulos por planta (NN), peso seco de ndulos
por planta (PSN), peso seco de raz (PSR) y peso seco de parte area (PST), en fase de
formacin de vainas, con excepcin de PSR que no present diferencias estadsticas en las
interacciones AC y AGC.
En la Tabla 9, se puede observar la comparacin de medias de NN y PSN,
a nivel de
ambientes y cepas. En general, la cepa J-01 obtuvo los valores ms bajos en NN y PSN,
mientras que las otras cepas y la mezcla mantuvieron valores similares. Estos resultados hacen
notorias las ventajas comparativas promedias de J-96 y J-98 sobre la cepa comercial J-01. Es
razonable pensar que la cepa J-96 realiza la mayor funcin simbitica cuando sta se mezcla
con J-01, por cuanto en todos los casos, el comportamiento de la mezcla fue similar a J-96.
Estos resultados coinciden con lo expuesto por Kober et al. (2004), donde algunas mezclas de
cepas pueden resultar menos infectivas y efectivas en la fijacin del N, o en ocasiones una de
ellas pude predominar.
Segn Neyra (1995), el peso de los ndulos de una planta est relacionada en forma directa
con la actividad de la fijacin del nitrgeno y el nmero de ndulos es un indicativo de la
condicin de la raz para la infeccin de los rizobios. Sin embargo, estos dos factores no son
completamente independientes. Cuando un planta tiene pocos ndulos, en general estos son
ms grandes que cuando los ndulos son abundantes. El grado de asociacin entre NN y PSN
fue positiva y significativa (r=0,67). En contraste, Iribarne et al. (1998), reportaron una baja
correlacin de estas variables, al encontrar que el mayor peso seco de ndulos no fue el
resultado del desarrollo de un nmero significativamente mayor de ndulos.
Tabla 9. Nmero y peso de ndulos en fase de formacin de vainas, en 8 ambientes de la
Orinoquia colombiana.
Loc/cepa
HOR08B
LIB07A
LIB07B
LIB08A
LIB08B
SCRUZ08A
TAL08A
TAL08B
Promedio
CV (%)
J-01
8,1
3,0
26,6
5,1
0,8
15,5
17,2
45,4
15,2 b
26,19
HOR08B
LIB07A
LIB07B
LIB08A
LIB8B
SCRUZ08A
TAL08A
TAL08B
Promedio
CV (%)
J-01
0,06
0,04
0,22
0,04
0,03
0,39
0,26
0,16
0,15 b
25,19
Nmero de ndulos/planta
J-96
J-96+J-01
48,7
46,4
14,2
19,5
45,5
41,4
31,7
25,7
25,0
22,2
15,2
20,3
52,7
50,2
70,3
76,4
37,9 a
37,8 a
R2
0,88
Peso seco de ndulos (g/planta)
J-96
J-96+J-01
0,29
0,26
0,16
0,19
0,33
0,32
0,18
0,17
0,34
0,24
0,29
0,45
0,22
0,23
0,21
0,21
0,25 a
0,26 a
R2
0,85
J-98
52,3
8,1
19,4
25,1
31,2
23,9
62,0
77,5
37,4 a
Promedio
38,9 b
11,2 e
33,2 bc
21,9 cd
19,8 ed
18,7 de
45,5 b
67,4 a
J-98
0,39
0,10
0,17
0,17
0,22
0,30
0,30
0,22
0,23 a
Promedio
0,25 ac
0,12 d
0,26 ab
0,14 cd
0,21 bd
0,36 a
0,25 ab
0,20 bd
*Letras diferentes entre promedios representan significancia estadstica (segn pueba de Tukey P<0,05).
En la Tabla 10, se presentan los promedios por localidad y cepa para PSR, PST y Rend. El
peso seco de raz (PSR), en promedio por localidades fue similar estadsticamente entre cepas
y con diferencias sobre la aplicacin de N, con el cual se obtuvo mayor peso de races. El PSR
correlacion positiva y significativamente con NN con un valor muy bajo (r=0,14; P<0.01), con
PST y Rend, esta correlacin fue de 0,47 y 0,11, respectivamente, razn por la cual, PSR se
considera una variable de bajo inters para los propsitos de evaluacin y seleccin de cepas.
La asociacin del PST con NN y PSN fue de r=0,20 (P<0.01) y con rendimiento de grano
(Rend), de r=0,33 (P<0.01). Estos resultados sugieren que, la variable PST aunque de
correlacion positiva con Rend,
evaluacin y seleccin por cuanto los ndulos solo explican el 20% del PST y esta ltima solo
explica el 33% del rendimiento de grano. En contraste, Iribarne et al. (1998) encontraron
correlacin positiva entre el peso seco de ndulos y el peso seco de la parte area, con valores
entre 0,65 y 0,91 (P<0.05). Kumudini et al. (2001), reportaron en soya incrementos en
rendimiento de grano asociados con incrementos en materia seca. De igual manera,
Koutroubas et al. (1998), demostaron que las plantas de soya inoculadas acumulan ms
materia seca y N en los tejidos vegetativos en la R2 y R5.
rendimiento de grano correlacion positivamente con la materia seca en R2, la materia vegetal
seca y el contenido de N, en R5. Estos autores sugieren que la acumulacin de materia seca y
de N en etapas tempranas de crecimiento son factores importantes en la consecucin de altos
rendimientos. Similares resultados fueron reportados por Iribarne et al. (1998), donde la
inoculacin con diversas cepas result en incrementos variables de la produccin de materia
seca que sugiere una capacidad diferencial de los rizobios para fijar N. En general, se puede
establecer que la produccin de una mayor cantidad de biomasa es una caracterstica deseable
que contribuye con una mayor disponibilidad de N. En la Tabla 10, se puede observar que en
promedio de localidades, el PST obtenido con la cepa J-01 fue estadsticamente inferior a las
otras cepas y a la mezcla, mientras que con la aplicacin de N se logr la mayor acumulacin
de biomasa seca. En general, la localidad Lib07A tuvo los valores ms bajos en PSN y PST,
con bajo PSR, el cual a su vez, coincidi con el rendimientos de grano ms bajo.
Para determinar el efecto de las siembras de primer y segundo semestre sobre el PSN, se
realiz el anlisis combinado por semestre. No se encontraron diferencias estadsticas entre
semestres ni en la interaccin semestre por cepa, mientras que las diferencias entre cepas
fueron estadsticamente diferentes (P<0.01). Sin embargo, en el anlisis combinado de varianza
de PSN excluyendo el factor semestre y considerando todas las localidades, se observaron
diferencias significativas (P<0.01), entre ambientes, cepas y variedades y sus respectivas
interacciones. Esto sugiere una adaptacin varietal diferencial por semestre que enmascara el
efecto de semestres sobre las cepas. Para determinar las diferencias varietales con cada cepa,
Localidad
HOR08B
LIB07A
LIB07B
LIB08A
LIB8B
SCRUZ08A
TAL08A
TAL08B
Promedio
J-01
1,13
0,67
1,57
0,56
0,91
0,60
0,65
0,84
0,87 b
CV (%)
13,34
Localidad
HOR08B
LIB07A
LIB07B
LIB08A
LIB08B
SCRUZ08A
TAL08A
TAL08B
Promedio
J-01
5,47
3,06
6,74
4,00
4,17
6,98
4,30
5,51
5,03 c
CV (%)
11,49
Localidad
HOR08B
LIB07A
LIB07B
LIB08A
LIB08B
SCRUZ08A
TAL08A
TAL08B
Promedio
J-01
1,80
0,91
2,20
1,90
0,84
2,17
2,14
1,78
1,72 b
N
1,35
0,65
1,72
0,62
1,02
1,04
0,84
1,12
1,05 a
Promedio
1,10 b
0,63 d
1,52 a
0,58 d
0,96 bc
0,72 cd
0,74 cd
0,96 bc
R2
0,84
Peso seco de parte area (g/planta)
J-96
J-96+J-01
J-98
9,18
6,87
6,85
2,96
3,02
2,79
7,07
6,86
5,87
4,92
5,08
5,57
4,60
4,66
4,38
9,01
8,09
8,38
5,54
6,07
6,20
6,53
6,07
6,41
6,23 b
5,84 b
5,81 b
N
10,74
3,50
10,59
7,07
6,66
9,27
6,02
7,68
7,69 a
Promedio
7,82 a
3,07 d
7,43 ab
5,33 c
4,89 c
8,35 a
5,63 bc
6,44 ac
R2
0,87
Rendimiento de grano (t/ha)
J-96
J-96+J-01
J-98
2,20
1,98
2,07
1,50
1,83
1,72
2,55
2,25
2,56
2,58
2,70
2,52
1,46
1,43
1,47
2,31
2,19
2,01
2,27
2,24
2,32
1,76
1,78
1,76
2,08 a
2,05 a
2,05 a
N
2,32
1,02
2,70
1,92
1,32
2,09
2,11
1,82
1,91 a
Promedio
2,07 ab
1,39 c
2,45 a
2,32 ab
1,31 c
2,15 ab
2,22 ab
1,78 bc
CV (%)
17,90
R2
0,89
* Letras diferentes entre promedios representan significancia estadstica (segn pueba de Tukey P<0,05).
Tabla 11. Comparacin de medias para PSN de variedades y cepas, combinado de ocho
localidades, en oxisoles de la Orinoquia colombiana.
Var/cepa
J-01
J-98
J-96
J-01+J-96
C. Libertad 4
0,11
0,22
0,07
0,12
0,13 c
Orinoquia 3
0,10
0,21
0,26
0,26
0,21 b
Soyca P-34
0,16
0,22
0,28
0,30
0,24 ab
C. Sabana 7
0,14
0,21
0,29
0,29
0,23 ab
C. Superior 6
0,14
0,30
0,34
0,32
0,27 a
C. Taluma 5
0,12
0,22
0,30
0,25
0,22 ab
0,13 b
0,23 a
0,26 a
0,26 a
Promedio
Promedio
* Letras diferentes entre promedios representan significancia estadstica (segn pueba de Tukey P<0,05).
En promedio, la cepa J-01 alcanz el menor PSN difiriendo estadsticamente de las otras
cepas. Las cepas J-98, J-96 y mezcla, superaron a J-01 en PSN entre 76,9 y 100%. En
general, las variedades respondieron de forma similar en PSN con la cepa J-98, mientras que
con J-96 y Mezcla, se hicieron evidentes las diferencias entre variedades particularmente por el
menor PSN de la variedad C. Libertad 4. Al comparar el peso seco promedio por ndulo
(PN=NN/PSN), entre variedades y cepas, se observaron diferencias estadsticas entre
variedades, cepas y localidades, con interaccin significativa solamente entre variedades por
localidad.
En promedio, la cepa J-01 tuvo el mayor PN con 13,2 miligramos por ndulo,
mientras que la variedad con mayor PN fue C. Libertad 4 (Tabla 12), aunque esta no difiri de
Orinoquia 3. King y Purcell (2001), encontraron que en general, los ndulos ms grandes tienen
un mayor suministro de floema y por lo tanto, tienen proporcionalmente mayor capacidad de
fijar N en sus tejidos que los ndulos pequeos. Ellos sostienen que entre ms grandes y
numerosos sean los ndulos mayor tolerancia tendr la planta al estrs hdrico. Champion et al.
(1992), aunque no encontraron diferencias en nmero de nudos entre cepas, obtuvieron un alta
relacin positiva entre el ndulos grandes y eficiencia simbitica, atribuido a una efectiva
interaccin cepa-cultivar. Sin embargo, en nuestro caso a pesar de que la cepa J-01 logr
mayor PN, su reducido NN promedio le confiri desventajas notorias sobre las dems cepas.
Tabla 12. Peso seco relativo de ndulos en mg (PN) por variedad y cepa, combinado de ocho
ambientes.
cepa/var
J-01
J-96
J-98
J-96+J-01
Libert 4
14,5
14,6
7,6
12,5
Orinoq. 3
15,2
8,5
10,9
6,9
Soy P-34
12,6
8,1
8,5
6,7
Sabana 7
11,0
7,4
8,0
6,6
Super 6
11,6
6,8
7,2
5,4
Taluma 5
14,7
5,0
6,4
4,5
Promedio
12,3 a
10,4 ab
9,0 b
8,2 b
7,8 b
7,7 b
Promedio
13,2 a
8,4 b
8,1 b
7,1 b
* Letras diferentes entre promedios representan significancia estadstica (segn pueba de Tukey P<0,05).
Las diferencias estadsticas para PSN fueron tambin evidentes en la interaccin VcA (var-cepa
x ambiente). Para establecer esta interaccin, se emple el modelo SREG, el cual permiti
integrar los factores aditivos de ambientes (A), tratamientos (Vc) (combinacin var-cepa) y el
factor multiplicativo de interaccin VcA. Los dos primeros componentes principales para PSN
fueron estadsticamente significativos. Los componentes CP1 y CP2 explicaron el 78% de la
variacin total (CP1=61,9 y CP2=16,1%, respectivamente). Los valores de los dos primeros
componentes principales del anlisis SREG por variedad y cepa, se presentan en la Tabla 13.
La alta correlacin r=0,97 entre PSN promedio y su CP1, permitio confirmar con la menor
respuesta general de las variedades a la inoculacin con la cepa J-01, en los diferentes
ambientes. Estos componentes que definen la interaccin var-cepa por ambiente (VcA),
permiti deducir que, las variedades presentaron mejor asociacin con las cepas J-96, mezcla y
J-98 para PSN por tener valores CP1 positivos de magnitud variable, con excepcin de la
variedad C. Libertad 4 que present valores CP1 negativos, y solo respondi favorablemente a
la cepa J-98. Sin embargo, con la cepa J-96 se present la mayor proyeccin de los vectores
CP1, es decir mayor PSN y con un menor efecto CP2, es decir una mayor estabilidad, mientras
que con J-98 la respuesta fue ms inestable. El mejor comportamiento para la cepa J-01 la
obtuvo la variedad Soyica P-34 y la mejor respuesta general en PSN, se observ en la variedad
C. Superior 6 con las cepas J-96, mezcla y J-98.
Tabla 13. Componentes principales del anlisis SREG para peso seco de ndulos (PSN) por
variedad y cepa.
Variedad
Soyica P-34
Orinoquia 3
C. Libertad-4
C. Taluma-5
C. Superior-6
C. Sabana-7
Cepa
J-01
J-01
J-01
J-01
J-01
J-01
PSN
promedio
0,16
0,10
0,11
0,12
0,14
0,14
CP1
-0,179
-0,362
-0,332
-0,361
-0,348
-0,355
CP2
-0,271
0,013
0,073
0,173
0,433
0,018
C. Sabana-7
C. Superior-6
C. Taluma-5
C. Libertad-4
Orinoquia 3
Soyica P-34
J-96
J-96
J-96
J-96
J-96
J-96
0,29
0,34
0,30
0,07
0,26
0,28
0,231
0,350
0,331
-0,363
0,152
0,220
0,060
0,163
0,056
-0,253
0,075
0,003
Soyica P-34
Orinoquia 3
C. Libertad-4
C. Taluma-5
C. Superior-6
C. Sabana-7
J-98
J-98
J-98
J-98
J-98
J-98
0,22
0,21
0,22
0,22
0,29
0,21
0,064
0,049
0,031
-0,009
0,225
0,091
-0,249
-0,275
-0,298
-0,042
-0,012
-0,252
C. Sabana-7
C. Superior-6
C. Taluma-5
C. Libertad-4
Orinoquia 3
Soy P-34
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
0,29
0,32
0,25
0,12
0,26
0,30
0,214
0,268
0,146
-0,347
0,109
0,175
0,110
0,120
0,188
0,068
0,099
0,000
variedades y cepas puede ser establecida con la comparacin de medias en cada variable
independiente. Las cepas en todos los casos fueron inferiores al control (150 Kg. de N.ha-1),
para la fase fenolgica de formacin de vainas (Figura 10). Las cepas J-01, J-98, J-96 y Mezcla,
fueron superadas significativamente por el control con N en 37, 24, 22 y 20 %, respectivamente,
en el contenido de nitrgeno de la parte area de la planta, lo que significa que la cepa J-01
aport la menor cantidad de nitrgeno a la planta, mientras que las otras cepas y la mezcla
hicieron aportes similares. En promedio, las cepas evaluadas no fueron capaces de suministrar
tanto N como lo hizo el testigo control (150 Kg.ha-1 de N). Segn Racca y Collino (2005), el
suministro de N es factor clave en el rendimiento de grano por cuanto cada Kg. de N acumulado
en biomasa produce un incremento de 8-10 kg de granos.ha-1. La asociacin entre rendimiento
de grano en soya y acumulacin de N en la planta, sugiere que el rendimiento podra ser
incrementado por nodulacin de cepas superiores de B. japonicum para fijar N (Bruce y
Fuhrmann, 1993), como sucede con las cepas J-96 y J-98, al compararlas con la ICA J-01
(cepa comercial de alta adopcin en la Orinoquia) .
0,40
a
Contenido de nitrgeno
0,35
0,30
0,25
0,20
0,15
0,10
0,05
0,00
J01
J96
J98
J96+J01
Fuente de nitrgeno
150 Kg.ha-1
Figura 10. Contenido de nitrgeno promedio (g.planta-1), en soya, aportado por las cepas de B.
japonicum vs control con 150 Kg de N.ha-1, en oxisoles de la Orinoquia colombiana.
El contenido de N entre variedades fue estadsticamente diferente (P<0.01). La variedad C. La
Libertad 4 obtuvo el menor contenido de N, aunque no difiri de C. Taluma 5.
variedades fueron estadisticamnte superiores a las dos anteriores (Figura 11).
las otras
0,35
Contenido de nitrgeno
0,30
ab
0,25
0,20
0,15
0,10
0,05
0,00
C. Libertad 4
Orinoquia 3
Soyica P-34
C. Sabana 7
C. Superior 6
C. Taluma 5
Variedades
Sin
embargo, aunque las plantas de soya pueden adquirir una parte sustancial de su N a partir de la
asociacin simbitica con B. japonicum, el manejo y las variables ambientales afectan esta
simbiosis. Cualquier factor que afecte la disponibilidad y absorcin del N del suelo, puede influir
de manera directa en la fijacin de N y la produccin de soya (George et al., 1988).
Los resultados demostraron que las cepas en evaluacin no suplieron los requerimientos totales
de N de las variedades, en las condiciones ambintales bajo estudio. S asumimos que los
niveles de N del suelo fueron muy bajos, las cepas aportaron entre 126 y 139,5 Kg. de N.ha-1.
Por lo anterior, se justifica la aplicacin de N mineral en niveles que no inhiban el desempeo
de los rizobios, ms aun en los oxisoles de la altillanura donde la disponibilidad de N del suelo
es muy baja. Segn Silvester (1983), la aplicacin de 30 Kg. de N.ha-1 puede favorecer el
crecimiento inicial de las plantas sin afectar la nodulacin.
0,94
J96
J98
0,92
0,90
0,88
0,86
0,84
0,82
0,80
0,78
J01
J96+J01
Cepas de B. japonicum
Figura 12. ndice relativo de fijacin biolgica de nitrgeno (IRFBN) de cepas de B. japonicum
sobre soya, en oxisoles de la Orinoquia colombiana. Fase fenolgica de formacin de vainas
(C.V (%)=15.1 y R2 = 0.81).
Estos aportes de N tambin estn en el rango reportado por Agha et al. (2004), quienes
sostienen que la soya puede fijar entre 14 y 300 Kg. de N.ha-1, dependiendo del potencial de
produccin, la disponibilidad de N del suelo y la interaccin gentica entre el genotipo husped
y las cepas de B. japonicum. Adems, segn Racca y Collino (2005), la fijacin biolgica del
nitrgeno (FBN) en soya, no es constante a lo largo de la ontogenia de la planta: es escasa y
lenta del comienzo del cultivo hasta R3, crece luego de manera exponencial hasta R5-R6, para
luego disminuir y anularse como consecuencia de la competencia que ejercen los granos por
los fotoasimilados.
Es de resaltar que el comportamiento de la cepa J-96 y la mezcla (J-96+J-01) fueron similares y
la efectividad de la J-96 no fue deprimida por la mezcla con un cepa menos efectiva como la J01, mientras que la mezcla parece haber inhibido a la J-01. Segn Daramola et al. (1994), en
una mezcla de cepas, una puede deprimir significativamnte a otra. En nuestro caso, la J-01 con
crecimiento ms lento que J-96 pudo haber favorecido la infectividad y efectividad de esta
ltima.
En general, los factores que controlan la cantidad de N fijado incluyen: la disponibilidad de N en
el suelo, los determinantes genticos de la compatibilidad entre planta y cepa, el potencial de
rendimiento de la planta y la ausencia de otros factores que limitan el rendimiento (Keyser y Li,
1992). En nuestro caso, un factor ambiental determinante en el IRFBN fue el estres hdrico. El
menor valor de IRFBN se obtuvo en la localidad LIB08B donde el dficit hdrico fue marcado en
llenado de grano y donde la cepa J-01 alcanz solo un IRFBN de 0.62. Estos resultados
ratifican las investigaciones realizadas por King y Purcell (2001), quienes encontraron que el
dficit hdrico reduce la fijacin simbitica y por ende, limita el rendimiento de grano y de
protena. Este dficit hdrico produce su mayor efecto en floracin y llenado de grano, por
cuanto, la mayor fijacin de N se produce en este perodo (Fabre y Planchon, 2000). Adems,
el establecimiento efectivo de la fijacin simbitica del nitrgeno en las zonas donde la soya no
ha sido cultivada, puede ser difcil, particularmente si la soya se planta en suelos tropicales
(Borkert y Sfredo 1995), donde, es necesario complementar la biofertilizacin con la fertilizacin
nitrogenada, como lo sugiere Ruz et al. (2008), para suplir los requerimientos de la planta.
mayor
Tabla 14. Contenido de protena (%) y protena total (t.ha-1) de variedades de soya inoculadas
con diferentes cepas-N, en oxisoles de la Orinoquia colombiana.
Cepa-N
Libertad 4
Orinoq. 3
Soy. P-34
Sabana 7
Superior 6
Taluma 5
Promedio
J-01
36,3 (0,59)
36,2 (0,48)
34,6 (0,48)
37,5 (0,57)
35,5 (0,55)
33,0 (0,53)
35,5 (0,53)
J-96
34,4 (0,45)
40,2 (0,76)
37,9 (0,66)
36,8 (0,72)
39,0 (0,72)
38,3 (0,81)
37,8 (0,69)
J-98
38,2 (0,68)
37,6 (0,64)
37,1 (0,60)
35,9 (0,59)
38,7 (0,64)
36,2 (0,73)
37,3 (0,65)
J-96+J-01
33,9 (0,41)
37,1 (0,70)
37,2 (0,64)
38,7 (0,70)
39,1 (0,69)
35,9 (0,79)
37,0 (0,65)
36,0 (0,75)
37,2 (0,68)
34,7 (0,60)
37,8 (0,64)
37,5 (0,69)
38,4 (0,71)
36,9 (0,68)
Promedio
35,8 (0,58)
37,7 (0,65)
36,3 (0,60)
37,3 (0,64)
38,0 (0,66)
36,4(0,71)
genotipos por alto contenido de protena resulta en decrementos del rendimiento de grano. Sin
emabargo, Yin y Vyn (2005), sugieren que la alta produccin de grano de soya puede estar
acompaada por altas concentraciones de isoflavonoides sin una sustancial reduccin de la
concentracin de la protena.
Egli y Bruening (2007), sugieren que la relacin negativa entre la concentracin de protena en
las semillas y el rendimiento no puede ser explicado solo por la semilla, al parecer porque se
trata de un fenmeno de toda la planta, que posiblemente impliquen reducciones en el
suministro de asimilados a las semillas cuando las plantas acumulan ms N. La concentracin
de N en granos maduros de soya es en parte determinado por las condiciones ambientales. En
promedio, los contenidos de protenas fueron bajos debido a la mala distribucin de las lluvias
durante el desarrollo de los granos, en los ambientes de estudio.
4.1.7. Componentes de rendimiento de grano.
El anlisis estadstico combinado de varianza para nmero de vainas por planta (VT), nmero
de granos por vaina (GV) y peso de grano (PG), exhibieron diferencias (P< 0.01), entre
ambientes(A), variedades (G) y en la interaccin GA. Entre cepas-N (C) e interaccin CA, las
diferencias fueron significativas para VT y GV, y no significativas para PG. Para la interaccin
GC estas diferencias fueron estadsticamente significativas para las tres variables. Los
componentes de rendimiento V3, VV y AE, difirieron estadsticamente (P< 0.01), entre
ambientes (A) y entre variedades (G). En la interaccin CG y CA, se presentaron diferencias
para V3 y AE. No se observaron diferencias en la interaccin GCA.
En VT, la variedad C. Libertad 4 present el menor nmero de vainas por planta promedio
(VT=27,7), diferente estadsticamente del resto de variedades. Aunque C. Sabana 7 super a
Orinoquia 3, estas fueron similares a las otras variedades. Entre cepas, las diferencias fueron
evidentes entre J-98 y mezcla con respecto a J-01, la cual fue superada en 14 y 12%,
respectivamente. Los promedios de VT con la cepa J-96 y N, fueron similares a las dems
(Tabla 15). C. Libertad 4, present la mayor proporcin de vainas de dos granos (44,1%),
mientras que en el resto de variedades, la mayor proporcin fue de tres granos por vaina (entre
51,4 y 57,7%). Con respecto al nmero de vainas de 4 granos, fue notoria la superioridad de
Soyica P-34 con un 13,6 % del total de vainas (Tabla 16). Estos resultados coinciden con lo
expuesto por la FAO (1995), quienes exponen que la forma de los foliolos y el nmero de
granos por vaina no son heredados independientemente. Las plantas con fololos normales y
ovales producen dos o tres semillas por vaina, las plantas con fololos estrechos producen tres
a cuatro granos por vaina como en la variedad Soyica P-34, y plantas con foliolos ovales con
uno o dos granos por vaina como en C. Libertad 4. La presencia de vainas de 4 granos de las
variedades con fololos estrechos o lanceolados, no ha demostrado tener relacin con mayor
productividad. Ms aun, este tipo de genotipos se los ha asociado con permitir la entrada de luz
hasta el suelo y favorecer la germinacin y desarrollo de malezas que competirn con la soya
por agua, radiacin solar y nutrientes.
El vaneamiento de vainas fue ms alta en las variedades C. Libertad 4 y Soyica P-34, atribuido
en gran medida a la presencia de Cercopora sojina en la primera variedad al reducir el rea
fotosinttica y por la presencia del sndrome de tallo verde en la segunda variedad (STV). En
las otras variedades, la proporcin de VV estuvo entre 4.7 y 4.9, lo que significa que factores
biticos y/o abiticos afectan el llenado de las vainas. Las variedades C. Superior 6 y Orinoquia
3 presentaron el menor vaneamiento, significativamente inferior a C. Libertad 4 y Soyica P-34.
En AE,
J-01
J-96
J-98
J-96+J-01
N
Promedio
Lib4
24,0
23,7
32,8
25,7
32,4
27,7 c
J-01
J-96
J-98
J-96+J-01
N
Promedio
2,43
2,44
2,50
2,42
2,41
2,44 d
J-01
J-96
J-98
J-96+J-01
N
Promedio
14,79
14,11
14,35
14,75
14,11
14,42 a
Tal5
31,5
34,2
36,4
37,2
31,3
34,1 ab
Promedio
30,6 b
33,4 ab
34,9 a
34,2 a
32,1 ab
2,57
2,62
2,66
2,63
2,64
2,63 b
2,57 b
2,60 b
2,61 a
2,61 a
2,60 ab
13,00
13,63
13,50
13,45
12,89
13,29 bc
13,55 ab
13,71 a
13,88 a
13,91 a
13,14 b
El nmero de vainas totales por planta (VT), present correlacin positiva con el rendimiento de
grano (r=0,81). Resultados similares fueron encontrados por Ponce et al. (2002) y Hernndez
et al. (2004). Los altos rendimientos de grano en soya sern posibles cuando se optimicen las
condiciones para el desarrollo de un mayor nmero de nudos sobre los cuales se desarrollaran
los frutos, lo que representara un aumento en el nmero de destinos o vainas totales, como lo
menciona Salvagiotti (2009).
Tabla 16. Proporcin de vainas con un grano (V1), dos (V2), tres (V3) y cuatro granos (V4) por
vaina del total de vainas por planta (VT).
Var/vainas
C. Libertad 4
Orinoquia 3
Soyica P-34
C. Sabana 7
C. Superiro 6
C. Taluma 5
V1
3,4
2,5
3,0
3,3
4,5
2,1
V2
44,1
34,0
25,6
31,0
34,5
31,3
V3
43,9
56,2
51,4
59,4
55,6
57,7
V4
1,3
2,5
13,6
1,4
0,8
3,9
En granos por vaina (GV), la variedad Soyica P-34 con un promedio de 2,75 super al resto de
variedades por su alta proporcin de vainas con 4 granos, mientras que el menor valor lo
obtuvo C. Libertad 4 con 2,44 granos por vaina por su alta proporcin de vainas de 2 granos. En
promedio, la cepa J-01 y J-96, obtuvieron el menor nmero de granos por vaina y el mayor fue
para J-98, aunque este ltimo no difiri de la mezcla. Ninguna cepa difiri en GV del testigo con
aplicacin de N. Estos resultados sugieren que los GV puede ser variable no solo
entre
variedades sino tambin por la cepa utilizada. En promedio, la mejor respuesta se logr con la
J-98 y la menor con J-01. Segn Bastidas y Agudelo (1994), el dficit hdrico reduce el nmero
de granos por vaina y el peso de grano.
El peso de grano (PG), fue superior en Soyica P-34 y C. Libertad 4, con respecto a las dems
variedades. Sin embargo, a pesar de que estas dos variedades se han caracterizado por
producir granos de mayor tamao que las otras variedades (Valencia et al., 2006), una
reduccin de PG fue observada
similar a J-01.
Es de resaltar que la comparacin de medias de las caractersticas descritas, son informativas
de las tendencias, por cuanto se presentaron diferentes niveles de interaccin entre variedades,
cepas y ambientes.
Por lo anterior,
variedades por cepas), para vincular al anlisis de manera conjunta los efectos de variedadcepa (Vc) y ambiente (A), para realizar el anlisis de interaccin (VcA)
La interaccin VcA fue establecida con la aplicacin del modelo SREG. Para la variable VT, los
dos primeros componentes principales fueron estadsticamente significativos. Los componentes
CP1 y CP2 explicaron el 63.1% de la variacin total. La mayor interaccin se present en los
ambientes Lib7A y Lib08A, donde se presentaron las mayores diferencias varietales. Los Vc 7,
18, 24, 11, 14, 23, 10,3, 28 y 1, que corresponden a los vrtices del polgono externo, son los
de mayor interaccin. Con valores positivos de interaccin se presenta C. Sabana7 con J-96 y
J-98, C. Superior 6 con mezcla, Soyica P-34 con la mezcla, Orinoquia 3 con J-96, J-98 y
mezcla. Mientras qued CP1 negativos fueron para C. Libertad 4 con J-96 y Soyica P-34 con J01 (Figura 13). Los tratamientos con la cepa J-01 se agruparon en el margen izquierdo de la
Figura y dan una connotacin de baja respuesta de las variedades a esta cepa en lo referente a
la produccin de vainas. En general, una situacin similar se present con la aplicacin de N,
con excepcin de C. Superior 6.
La correlacin del componente CP1 del nmero de vaina totales (VT) con peso seco de ndulos
(PSN), fue positiva y estadsticamente significativa con r=0,77 y con el rendimiento de grano fue
de r=0,89. As mismo, el CP1 de PSN correlacion con rendimiento de grano en r=0,79, lo que
significa que los factores que intervienen en interaccin de PSN y VT con el ambiente son
similares a los que afectan el rendimiento de grano.
Figura 13. Biplot GGA para vainas totales (VT) de variedades de soya inoculadas con
diferentes cepas de B, japonicum-N, evaluadas en diferentes ambientes de la Orinoquia. Los
valores en color rojo, corresponden a los tratamientos con la cepa J-01.
Para granos por vaina (GV), los primeros dos componentes fueron significativos y explicaron el
82,3% de la variacin total (65,1% del CP1 y 17,2 del CP2). Los valores obtenidos de cada
componente se presentan en la Tabla 17.
present alta interaccin en la localidad Lib07A, seguido por Tal08A. Los tratamientos Vc 7, 24
y 25 correspondientes, tuvieron valores de CP1 positivos, mientras que para los Vc 3, 22, 20 y
23, dichos valores fueron negativos. Los Vc positivos estuvieron asociados con: la variedad C.
Sabana 7 con la cepa J-96, Soyica P-34 con la J-96, mezcla y con N. Los Vc negativos
correspondieron a C. Libertad 4 con J-01 y Mezcla.
Tabla 17. Componentes principales del anlisis SREG para granos por vaina (GV) por variedad
y cepa.
Variedad
Soyica P-34
Orinoquia 3
Corpoica Libertad-4
Corpoica Taluma-5
Corpoica Superior-6
Corpoica Sabana-7
Cepa
J-01
J-01
J-01
J-01
J-01
J-01
GV
2,71
2,60
2,43
2,58
2,56
2,59
CP1
0,141
-0,007
-0,206
-0,031
-0,074
0,002
CP2
-0,141
-0,010
-0,120
-0,070
0,075
0,037
Corpoica Sabana-7
Corpoica Superior-6
Corpoica Taluma-5
Corpoica Libertad-4
Orinoquia 3
Soyica P-34
J-96
J-96
J-96
J-96
J-96
J-96
2,64
2,56
2,64
2,44
2,62
2,77
0,034
-0,057
0,029
-0,175
0,000
0,186
0,137
0,024
0,030
-0,123
0,123
-0,035
Soyica P-34
Orinoquia 3
Corpoica Libertad-4
Corpoica Taluma-5
Corpoica Superior-6
Corpoica Sabana-7
J-98
J-98
J-98
J-98
J-98
J-98
2,77
2,64
2,51
2,68
2,54
2,58
0,189
0,013
-0,117
0,084
-0,086
-0,034
-0,051
0,104
-0,005
0,044
0,058
0,080
Corpoica Sabana-7
Corpoica Superior-6
Corpoica Taluma-5
Corpoica Libertad-4
Orinoquia 3
Soyica P-34
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
2,65
2,55
2,65
2,42
2,63
2,80
0,055
-0,071
0,042
-0,208
0,003
0,226
0,054
0,073
0,043
-0,078
0,129
0,000
Soyica P-34
2,82
0,240
-0,187
Orinoquia 3
-1
2,61
0,008
-0,004
-1
2,42
-0,195
-0,118
-1
2,66
0,071
-0,035
-1
2,56
-0,052
0,034
-1
2,59
-0,007
-0,070
Corpoica Libertad-4
Corpoica Taluma-5
Corpoica Superior-6
Corpoica Sabana-7
150 Kg. N. ha
150 Kg. N. ha
150 Kg. N. ha
150 Kg. N. ha
150 Kg. N. ha
En sntesis, el nmero de granos por vaina promedio es afectado por la interaccin VcA, donde
los mayores efectos se presentan para las variedades C. Libertad 4 y Soyica P-34, aunque con
efectos contrarios. En la primera variedad todos los valores de CP1 fueron negativos para
todas las cepas e inferiores a las dems variedades, mientras que para Soyica P-34, todos los
coeficientes CP1 fueron positivos y superiores al resto de variedades. Esta variable GV aunque
con positiva correlacin de su CP1 con rendimiento de grano, VT y PSN, los valores fueron
inferiores a 0,41.
Para peso de granos (PG), los dos primeros componentes principales fueron significativos y
explicaron el 79,2 % de la variacin total (60 y 19,2 % entre CP1 y CP2, respectivamente). Los
vrtices de interaccin positiva fueron los Vc 3, 13, 12, 24, 1 y 23 asociados con las variedades
C. Libertad 4, Soyica P-34 y Orinoquia 3. Los CP1 negativos fueron para los VC 30, 8 y 17,
asociados con las variedades C. superior 6 y C. sabana 7 (Tabla 18). Estos resultados sugieren
que el peso de grano es afectado por las condiciones ambientales, la fuente de N aplicada y
est altamente asociado con la variedad de soya.
Anlisis de correlacin de los componentes asociados con rendimiento de grano. Los
anlisis de varianza mostraron diferencias (P<0.01) entre los tratamientos (Vc), ambientes (A) e
interaccin (VcA). Los resultados obtenidos de la evaluacin de 30 tratamientos (Vc) como
promedios, combinado de ambientes y por semestre, la desviacin estndar, los valores
mnimos y mximos, se presentan en la Tabla 19. El rendimiento promedio (Rend) del anlisis
combinado y por semestre estuvo en un rango de 0,81 a 2,67 t.ha-1. Los promedios ms altos
estuvieron asociados a las combinaciones de las variedades Orinoquia 3 y C. superior 6 con
las cepas J-96 y J-98, y los ms bajos, estuvieron asociados con la cepa J-01 y la variedad C.
Libertad 4. Los componentes de rendimiento estuvieron en un rango de: nmero de nudos
(NM) , entre 10 y 19; granos por vaina (GV), entre 2,24 y 2,97; vainas por planta (VT) entre 17 y
39, y peso de grano (PG), entre 11,63 y 20 g.
Los coeficientes de correlacin entre los caracteres evaluados se presentan en la Tabla 20. En
el anlisis combinado, una positiva y significativa relacin entre el rendimiento de grano y sus
componentes se observ con el nmero de vainas con 3 granos (V3), nmero de vainas por
planta o vainas totales (VT) y nmero de granos por vaina (GV), y fue negativa y significativa
estadsticamente con nmero de vainas vanas (VV) y peso de grano (PG). Los componentes,
nmero de vainas con 1 grano (V1), con 2 granos (V2), con 4 granos (V4) y abortos
embrionarios (AE), no correlacionaron estadsticamente con rendimiento de grano.
Tabla 18. Componentes principales del anlisis SREG para peso de grano (PG) por variedad y
cepa.
Variedad
Soyica P-34
Orinoquia 3
Corpoica Libertad-4
Corpoica Taluma-5
Corpoica Superior-6
Corpoica Sabana-7
Cepa
J-01
J-01
J-01
J-01
J-01
J-01
GV
13,3
12,9
13,9
11,9
11,5
11,9
CP1
0,293
0,081
0,135
-0,127
-0,350
-0,145
CP2
-0,131
-0,081
0,506
-0,103
-0,008
-0,137
Corpoica Sabana-7
Corpoica Superior-6
Corpoica Taluma-5
Corpoica Libertad-4
Orinoquia 3
Soyica P-34
J-96
J-96
J-96
J-96
J-96
J-96
12,6
11,5
12,7
13,4
12,6
13,9
-0,006
-0,415
0,075
0,095
-0,070
0,370
-0,001
0,033
-0,126
0,369
-0,041
-0,073
Soyica P-34
Orinoquia 3
Corpoica Libertad-4
Corpoica Taluma-5
Corpoica Superior-6
Corpoica Sabana-7
J-98
J-98
J-98
J-98
J-98
J-98
14,1
13,0
13,3
12,4
11,8
12,5
0,344
0,003
0,211
0,003
-0,344
-0,029
0,101
0,088
-0,042
-0,080
0,137
-0,026
Corpoica Sabana-7
Corpoica Superior-6
Corpoica Taluma-5
Corpoica Libertad-4
Orinoquia 3
Soyica P-34
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
J-01 + J-96
12,8
11,9
12,4
13,7
12,4
14,0
0,038
-0,296
-0,063
0,264
0,036
0,368
-0,026
0,075
0,021
0,167
-0,330
-0,089
Soyica P-34
Orinoquia 3
Corpoica Libertad-4
Corpoica Taluma-5
Corpoica Superior-6
Corpoica Sabana-7
150 Kg N.ha-1
150 Kg N.ha-1
150 Kg N.ha-1
150 Kg N.ha-1
150 Kg N.ha-1
150 Kg N.ha-1
13,0
12,1
13,2
12,2
11,7
11,6
0,131
-0,152
0,211
-0,163
-0,341
-0,157
-0,118
-0,038
0,058
0,111
0,104
-0,321
Tabla 19. Media, desviacin eatndar, valores mnimos y mximos para diferentes caracteres
de soya del anlisis combinado y por semestres, en la Orinoquia colombiana.
Combinado semestre A y B
Variable
Nmero de nudos (NM)
Nmero de granos (GV)
Vainas totales(VT)
Peso granos (PG)
Rendimiento (Rend)
Media
15,30
2,61
33,09
13,74
1,96
14,88
2,51
33,21
16,09
1,97
15,62
2,68
33,00
13,15
1,90
D.estndar
2,17
0,10
3,97
0,78
0,29
Semestre A
2,44
0,13
5,97
1,16
0,51
Semestre B
2,03
0,10
3,45
0,93
0,23
Mnimo
11,11
2,42
23,02
12,60
1,24
Mximo
18,25
2,83
39,93
15,08
2,38
10,17
2,24
17,14
13,67
0,81
19,08
2,78
41,42
20,00
2,67
11,50
2,51
25,38
11,63
1,42
18,69
2,97
38,81
14,81
2,24
Las caractersticas agronmicas, nmero de nudos a madurez (NM), nmero de ndulos (NN),
peso seco de ndulos (PSN), peso seco de raz (PSR), peso seco de parte rea (PST) y
contenido de protena, correlacionaron positiva y significativamente con rendimiento de grano y
negativamente con la presencia de Cercopora sojina (CS), a diferencia de das a floracin
completa (DF), das a madurez (DM) y la reaccin a las enfermedades como pstula bacteriana
(PB) y roya (RY), que no presentaron asociacin estadsticamente significativa con rendimiento
de grano.
La positiva asociacin entre rendimiento de grano y altura de planta madurez (ALM), nmero de
nudos por planta (NM), nmero de vainas con 3 granos (V3), nmero de vainas por planta (VT),
peso seco de raz (PSR) y peso seco de parte area (PST), fue consistente entre semestres y
el anlisis combinado. Similares resultados fueron reportados por Oz et al. (2009).
r
combinado
-0,07
0,01
0,63
0,67
0,11
0,18
0,19
0,84
0,02
0,44
-0,33
0,28
0,81
-0,44
-0,69
-0,13
-0,29
0,70
0,63
0,38
0,65
Prob
Prob
Sem A
0,730
0,980
<,0001
<,0001
0,550
0,350
0,310
<,0001
0,910
0,020
0,080
0,130
<,0001
0,010
<,0001
0,490
0,120
<,0001
<,0001
0,040
<,0001
-0,01
0,01
0,61
0,64
0,26
0,18
0,46
0,84
0,00
0,53
-0,44
0,28
0,79
-0,21
-0,65
*
*
0,76
0,70
0,30
0,68
Prob
Sem B
0,965
0,950
<,0001
<,0001
0,173
0,345
0,010
<,0001
0,999
0,003
0,016
0,141
<,0001
0,276
<,0001
*
*
<,0001
<,0001
0,101
<,0001
-0,01
0,02
0,56
0,52
0,08
-0,08
0,18
0,53
0,07
0,18
0,23
0,46
0,63
-0,18
-0,35
-0,07
-0,26
0,27
0,13
0,51
0,733
0,946
0,921
0,001
0,004
0,688
0,660
0,332
0,002
0,702
0,337
0,224
0,011
<,0001
0,334
0,062
0,705
0,172
0,143
0,491
0,004
<,0001
Con relacin a roya (Phakopsora sp.) y pstula bacteriana (Xanthomanas axonopodis var.
Campestres), estas estuvieron ausentes en el primer semestre y en el segundo semestre, roya
solo se present en la localidad Lib08 y en severidades de hasta 25% del rea foliar afectada,
con sintomatologa similar en todas las variedades evaluadas. Es de resaltar, que por
evaluaciones previas del programa de mejoramiento de soya de Corpoica,
en el C.I La
Libertad, se ha establecido que las variedades desarrolladas para la Orinoquia colombiana son
susceptibles al patgeno. En el mundo, no existen variedades comerciales resistentes a roya
(McGinnis y Suszkiw, 2007). La presencia de pstula bacteriana fue tangencial y afect en
particular a la variedad Corpoica Sabana 7, sin comprometer el rendimiento de grano.
Las variables NN, PSN y PST, tuvieron efectos directos positivos e indirectos
notorios sobre el rendimiento va nmero de nudos por planta (NM). En general, el mayor efecto
indirecto de las variables evaluadas se present con nmero de nudos por planta (NM). Esta
caracterstica por ser altamente consistente con un aporte directo superior al 50% de la
correlacin total, fue includa en el anlisis de sendero como variable exgena para valorar el
aporte de los componentes de rendimiento sobre la produccin de grano.
Tabla 21. Coeficientes de sendero y participacin (%) con efectos directos e indirectos de los caracteres agronmicos sobre el
rendimiento de grano de Vc (variedades de soya combinadas con cepas de B. japonicum), en oxisoles de la Orinoquia
colombiana.
Carcter
DF
DM
ALM
NM
CS
NN
PSN
PSR
PST
-0,131
28,1
0,027
2,6
0,008
0,50
0,020
1,3
-0,039
2,8
-0,026
2,9
-0,010
1,3
-0,020
1,8
0,005
0,5
0,079
16,9
-0,378
36,8
-0,216
14,80
-0,208
13,7
0,121
8,8
0,023
2,5
0,026
3,3
-0,182
16,7
-0,057
5,6
0,004
0,9
-0,037
3,6
-0,066
4,50
-0,063
4,2
0,047
3,4
-0,014
1,5
-0,009
1,1
-0,045
4,1
-0,036
3,5
-0,121
25,9
0,443
43,1
0,774
53,00
0,806
53,1
-0,677
49,2
0,242
26,7
0,161
20,5
0,500
45,8
0,371
36,5
Cercospora (CS)
0,049
10,5
-0,052
5,1
-0,116
7,90
-0,137
9,0
0,163
11,8
-0,060
6,6
-0,059
7,5
-0,078
7,1
-0,073
7,2
0,049
10,5
-0,015
1,5
0,054
3,70
0,074
4,9
-0,091
6,6
0,246
27,1
0,222
28,2
0,025
2,3
0,074
7,3
0,017
3,6
-0,015
1,5
0,028
1,90
0,043
2,8
-0,078
5,7
0,194
21,4
0,216
27,5
-0,004
0,4
0,054
5,3
-0,004
0,9
-0,012
1,2
-0,018
1,20
-0,016
1,1
0,012
0,9
-0,003
0,3
0,001
0,1
-0,026
2,4
-0,017
1,7
-0,013
2,8
0,049
4,8
0,181
12,40
0,151
9,9
-0,148
10,8
0,099
10,9
0,082
10,4
0,211
19,3
0,329
32,4
Correlacin total
-0,070
0,010
0,630
0,670
-0,690
0,700
Nmeros con negrilla corresponden a los efectos directos y el resto, son efectos indirectos.
0,630
0,380
0,650
Combinado
Semestre A
Semestre B
0,60**
0,55**
-0,48*
0,67**
0,65**
0,46*
-0,35 ns
0,64**
0,47*
0,55**
-0,40*
0,52**
0,48*
0,14 ns
0,44*
0,71**
-0,39*
0,53**
0,06 ns
0,30 ns
0,18 ns
-0,29 ns
0,81**
-0,28 ns
0,79**
-0,41*
0,63**
-0,44*
-0,21 ns
-0,18 ns
es el
resultado de la adicin del efecto directo de cada variable con sus efectos indirectos. En el
anlisis de sendero conjunto de los efectos combinados y por semestre, el mayor aporte
en efectos directos fue para nmero de vainas por planta (VT) con coeficientes que
estuvieron entre 0,628 y 0,836, con una participacin de la correlacin total entre 61,3 a
77,5 %. Estos resultados coinciden con Milik et al. (2007), Malik et al. (2006), Kantolic y
Slafer (2005), Rajanna et al., (2000), Chand (1999) y Jadhav et al. (1995) y Mishra et al.
(1994), quienes concluyeron que el nmero de vainas es un importante componente del
rendimiento del cultivo de la soya (Glycine max (L.) Merrill). Con relacin al nmero de
vainas por planta (VT), las investigaciones de Oz et al. (2009), concluyeron que la
asociacin positiva (r=0,458; P<0.01) con rendimiento de grano fue estable a travs de
aos y localidades, y que el nmero de vainas por planta ocurri desde el efecto directo
en un 50,4%. Para Board et al. (1997), el nmero de vainas por nudo reproductivo fue el
mejor criterio de seleccin indirecto para rendimiento de grano.
En segundo lugar, se ubic el nmero de nudos a madurez (NM), con efectos directos de
0,26 a 0,524 y aportes entre 38,8 y 46%. En esta sentido, Bizeti et al. (2004), encontraron
que el nmero de nudos correlacion significativamente con rendimiento de grano y fue el
de mayor efecto directo sobre rendimiento, por lo que, para estos investigadores NM
podra ser gran ayuda para la seleccin indirecta de genotipos con alto rendimiento.
Las variables VT y NM no solo presentaron los mayores efectos directos sino tambin, los
mayores aportes en efectos indirectos, por lo que, se constituyen en dos variables
altamente tiles como predictoras confiables del rendimiento de grano de las variedades
de soya, en oxisoles de la Orinoquia.
El nmero de granos por vaina (GV), no fue relevante para la produccin de grano por
tener efectos directos bajos y en general, sus aportes en efectos indirectos fueron tambin
bajos, en el anlisis combinado y semestre B, aunque para el semestre A, el efecto directo
sobre el rendimiento fue negativo con un aporte del 28,5%. Estos resultados contrastan
con los obtenidos por Ball et al. (2001), Iqbal et al. (2003) y Oz et al. (2009) quines a
travs del anlisis de sendero encontraron que, el nmero de granos por vaina ejerci un
efecto directo positivo sobre el rendimiento de grano. Estas discrepancias pueden ser
atribuidas al efecto ambiental y a que tradicionalmente la variable es medida como el
promedio de nmero de granos por vaina. En nuestro estudio, no se encontr correlacin
entre el rendimiento de grano con vainas de 1 grano (V1), de 2 granos (V2) ni con vainas
de 4 granos (V4), pero fue estadsticamente significativa con vainas de 3 granos. Esta
variable present consistencia a travs de localidades y combinacin de variedades con
cepas, con una alta asociacin con el nmero de vainas por planta (VT) y rendimiento de
grano (Figura 14), lo que permite deducir que V3 y VT son muy importantes, muy
predictivas del rendimiento de grano y de gran utilidad como criterios de seleccin en
programas de mejoramiento de la especie, en oxisoles la Orinoquia.
EFECTOS
Efectos directos
Efectos indirectos
-por nmero de granos (GV)
0,001
0,1
-0,248
-por vainas totales(VT)
0,345
51,5
0,385
-por peso granos (PG)
0,064
9,6
-0,021
Correlacin total
0,67
0,64
Nmero de vainas totales (VT)
Efectos directos
0,628
77,5
0,836
Efectos indirectos
-por nmero de nudos (NM)
0,143
17,6
0,241
-por nmero de granos (GV)
0,001
0,1
-0,271
-por peso granos (PG)
0,039
4,8
-0,016
Correlacin total
0,81
0,79
Nmero de granos por vaina (GV)
Efectos directos
0,002
0,3
-0,381
Efectos indirectos
-por nmero de nudos (NM)
0,156
32,7
0,340
-por vainas totales(VT)
0,301
63,1
0,594
-por peso granos (PG)
-0,019
3,9
-0,023
Correlacin total
0,44
0,53
Peso de 100 granos (PG)
Efectos directos
-0,133 30,3
0,059
Efectos indirectos
-por nmero de nudos (NM)
-0,125 28,3 -0,183
-por nmero de granos (GV)
0,000
0,0
0,149
-por vainas totales(VT)
-0,182 41,3 -0,234
Correlacin total
-0,44
-0,21
R2
Sem B
Coefic.
%
0,340
46,0
21,1
32,7
1,7
-0,033
0,289
-0,077
0,52
4,4
39,1
10,4
61,3
0,526
66,1
17,7
19,8
1,2
0,187
-0,004
-0,079
0,63
23,5
0,5
9,9
28,5
-0,069
21,7
25,4
44,4
1,7
0,160
0,032
0,058
0,18
50,2
9,9
18,1
9,4
0,193
34,1
29,3
23,8
37,5
-0,136
-0,021
-0,216
-0,18
24,1
3,7
38,2
50,00
3,00
Vainas 3 granos
Vainas Totales
45,00
Rendimiento
2,50
40,00
35,00
2,00
30,00
1,50
25,00
20,00
1,00
15,00
10,00
0,50
5,00
N
J0
1
J9
6
J9
8
M
ez
N
J0
1
J9
6
J9
8
M
ez
N
J0
1
J9
6
J9
8
M
ez
N
J0
1
J9
6
J9
8
M
ez
N
J0
1
J9
6
J9
8
M
ez
0,00
N
J0
1
J9
6
J9
8
M
ez
J0
1
J9
6
J9
8
M
ez
0,00
Figura 14. Asociacin del nmero de vainas con 3 granos (V3), con vainas totales (VT) y
rendimiento de grano, a travs de ambientes y cepas de B. japonicum -N.
Con respecto a peso de grano (PG), los efectos directos aunque algunos negativos,
fueron poco relevantes para afectar la variable endgena. Los efectos directos estuvieron
entre -0,133 a 0,193, considerados como valores muy bajos de interrelacin. Taware et al.
(1997), Shukla et al. (1998) y Board et al. (1999), no encontraron relacin significativa
entre peso de grano y rendimiento. En contraste, Malik et al. (2007), en sus
investigaciones con soya, revelaron con el anlisis de sendero que el peso de semillas
tuvo un efecto directo negativo sobre el rendimiento. Estas diferencias en los resultados
puede ser atribuidos a la influencia de factores ambientales e interacciones en esta
variable tan compleja, y su expresin depende de factores como: la variedad, las
condiciones adafoclimticas predominantes (agua, radiacin, nutricin, etc) y presencia de
factores biticos como patgenos y dems plagas.
Las correlaciones genotpicas y fenotpicas de los componentes de rendimiento se
presentan en la Tabla 24. En general, las correlaciones genotpicas fueron superiores a
sus respectivas fenotpicas, por lo que segn Iqbal et al. (2003), es un indicativo de una
mayor contribucin de factores genticos sobre la asociacin; sin embargo, varios
coeficientes fenotpicos con menor magnitud resultaron estadsticamente significativos a
diferencia de los coeficientes genotpicos, sugiriendo un alto efecto ambiental sobre los
Los
NM
1
0,847*
0,602
0,823*
-0,240
0,853*
V3
0,893
1
0,599
0,980**
-0,440
0,868*
NM
1,028**
V3
GV
0,906*
0,504
1
0,258
VT
PG
Rend
0,567
-0,369
0,802
0,835*
-0,747
0,743
NM
V3
GV
VT
PG
Rend
1
0,782
0,628
0,711
-0,643
0,610
0,832
1
0,582
0,938**
-0,938
0,870*
NM
V3
GV
VT
PG
Rend
Cepa J-96
GV
0,626
0,631
1
0,728
0,478
0,583
Cepa J-98
0,546
0,241
1
0,116
0,471
0,017
Cepa J-01
0,694
0,589
1
0,779
-0,113
0,323
VT
0,918
1,005**
0,848
1
-0,311
0,870*
PG
-0,255
-0,606
0,619
-0,447
1
-0,395
Rend
0,909
0,923
0,634
0,954
-0,616
1
1,000**
-0,432
1,146**
1,000**
1,000**
-1,129
0,611
1,092**
0,114
1
-0,658
0,314
1,000**
1
-0,691
1,000**
1,000**
1
0,777
0,962
0,870
1
-0,868
0,797
-0,730
-1,075
-0,120
-1,114
1
-0,522
1,338*
1,851*
0,627
1,777*
1,000**
1
e interaccin GA y su participacin de la
varianza total. La localidad o ambiente fue la fuente de mayor aporte en la variacin del
rendimiento de grano con una participacin del 35,3 al 73,6% de la variacin total. Auque
el ambiente present el mayor aporte a la varianza total, la magnitud de la interaccin
relativa del 21,5 al 34,5%, hace irrelevante el efecto ambiental como factor individual, y
justifica la consideracin de la interaccin GA como elemento clave en la evaluacin y
seleccin de las variedades de mejor comportamiento. Segn Yan et al. (2000) y
Dehghani et al. (2006), la alta magnitud de la interaccin con respecto a G, sugiere la
posible existencia de diferentes mega-ambientes, es decir, grupos de ambientes con
respuesta similar pero que difieren en su comportamiento con respecto a otros. Los
efectos de interaccin (GA) fueron superiores al efecto principal genotpico (G) en las
cepas-N, con excepcin de la J-96 y la mezcla, donde este efecto fue inverso. Con estos
resultados, qued demostrado que independiente de la cepa-N utilizada, el rendimiento de
grano de las variedades es afectada por las condiciones ambientales y que stas no
responden de manera similar a los cambios ambientales por efecto de la interaccin GA.
La mayor participacin del componente gentico sobre la varianza total se obtuvo en la
cepa J-96 y en mezcla de sta con J-01 (36,4 y 34,9%, respectivamente), aunque menor
interaccin GA se observ con la cepa J-96. Al ser menos afectada por el ambiente y por
la interaccin, con la J-96 las variedades manifestaron mayor estabilidad promedia,
mientras que con la cepa J-98 se obtuvo la mayor participacin de la interaccin (34,5%).
Con la aplicacin de N se present una menor interaccin GA, pero con el mayor aporte
de la varianza por el factor ambiental. Este ltimo al igual que J-01, tuvieron los menores
aportes en la suma de cuadrados entre variedades.
La contribucin relativa del efecto principal G permite una inferencia en sentido amplio o
comportamiento general de variedade a travs de ambientes, mientras que el efecto de la
interaccin GA est relacionado con una inferencia de tipo estrecha sugiriendo una
adaptacin especfica. Aunque para todos los casos, la interaccin fue estadsticamente
significativa, la cepa J-96 y mezcla, tuvieron una interaccin GA inferior al efecto G y por
lo tanto, una adaptacin amplia, mientras que con las dems cepas y la aplicacin de N,
los mayores efectos fueron ambintales. Desde el punto de vista fenotpico, este
fenmeno se evidenci con la presencia de sntomas de deficiencia de N marcada en la
variedad C. Libertad 4 con la cepa J-96 y mezcla, en todas las localidades. En general,
las cepas como elemento del ambiente generaron una respuesta diferencial entre
variedades.
Tabla 25. Suma de cuadrados por componente y proporcin de la varianza total de los
efectos principales de variedades (G), ambientes (A), e interaccin (GA), para rendimiento
de grano de variedades de soya por cepa de B. japonicum, en diferentes ambientes.
Componentes de varianza
SC IV
Cepa - N
% de la varianza total
GA
GA
J-01
43,76
5,25
17,78 i
65,5
7,9
26,6
J-98
23,27
5,65
15,25
52,7
12,8
34,5
J-96
29,65
26,47
16,60
40,8
36,4
22,8
20,62
20,40
17,42
35,3
34,9
29,8
45,30
3,00
13,21
73,6
4,9
21,5
J-01+J-96
150 kg de N. ha
-1
CP1: 54.6
CP2: 25.3
Acum : 79.9%
J-01
CP1: 68.6
CP2: 20.1
Acum : 88.7%
J-98
CP1: 85.4
CP2: 8.8
Acum : 94.2%
CP1: 63.0
CP2: 19.8
Acum : 82.8%
J-96
150 Kg N ha -1
Figura 15. Anlisis biplot GGA del rendimiento de grano de 6 variedades de soya por
cepa-N, evaluadas en ocho ambientes de la Orinoquia colombiana. A. cepa J-01, B. Cepa
J-96, C. Cepa J-98 y D. 150 Kg.de N. ha-1
En los biplot GGA, la cercana de los vectores ambientales es sinnimo de alta correlacin
entre ellos (Yan y Rajcan, 2002b). Un ngulo inferior de 90 o superior a 270 entre un
vector genotpico y un vector ambiental indica que el genotipo tiene una respuesta positiva
en ese ambiente (Yang et al., 2009). En la cepa J-96 estos ejes se encuentran ms
relacionados, conformando solo tres sectores de ambientes que comparten condiciones
similares. Esto ratifica el comportamiento ms estable de las variedades con la cepa J-96.
Los resultados del biplot coinciden con el anlisis de los componentes de varianza, donde
la varianza A fue inferior a G, comparativamente con el resto de cepas.
En los biplot, es posible identificar las variedades de menor interaccin GA con cada cepa,
por su menor CP1 y CP2. La variedad C. Sabana 7 con las cepas J-01 y J-98, Soyica
P-34 con J-96 y C. Taluma 5 con la aplicacin de N.
ambientes (A) e interaccin CA, para rendimiento de grano. Estos resultados sugieren una
respuesta variable de las cepas ante los cambios de las condiciones ambientales y ms
aun, un comportamiento diferencial atribuido a la interaccin Cepa-N x ambiente (CA), lo
que significa, que el nivel promedio de cada cepa-N no puede ser utilizado como criterio
de seleccin de la mejor cepa, puesto que est implcito que cada cepa-N responde
diferencialmente al ambiente y el orden de respuesta de las cepas-N varia de un ambiente
a otro. Al separar los efectos principales cepa-N (C), ambientes (A) e interaccin CA, el
mayor aporte a la varianza total lo hizo el ambiente en un rango de 56,1 a 75,6%, seguido
por el aporte de la interaccin CA con 15,6 a 26,4 %. Los menores aportes a la suma de
cuadrados se observ entre cepas-N (Tabla 26). Estos aportes significativos de la
interaccin CA, demuestran su importancia en el proceso de seleccin de cepas.
Independientemente de las variedades, las cepas-N son afectadas por las condiciones
ambientales variables y con respuesta diferencial a estos cambios.
SC IV
% de la varianza total SC
Soyica P-34
A
25,50
C
4,67
CA
9,01 i
A
65,1
C
11,9
CA
23,0
Orinoquia 3
31,30
10,93
13,61
56,1
19,6
24,4
C. Libertad 4
40,28
7,20
8,78
71,6
12,8
15,6
C. Taluma 5
17,27
3,91
7,61
60,0
13,6
26,4
C. superior 6
42,77
3,47
10,32
75,6
6,1
18,2
68,9
13,7
17,4
C. Sabana 7
25,72
5,12
6,52
Componente de varianza para ambientes.
Componente de varianza para cepas por variedad.
i Componente de varianza de la interaccin ambiente x cepa
Soyica P-34
Orinoquia 3
CP1: 70.4
CP2: 19.0
Acum : 89.4%
CP1: 70.0
CP2: 19.8 ns
Acum : 89.8%
Cepa
Vc
J-01
1
J-96
12
J-98
13
Mez
24
N
25
J-01
2
J-96
11
J-98
14
C. Libertad 4
C. Taluma 5
CP1: 67.1
CP2: 24.0
Acum : 91.1%
CP1: 67.0
CP2: 19.2
Acum : 86.2%
Cepa
Vc
J-01
3
J-96
10
J-98
15
Mez
22
N
27
J-01
4
J-96
9
J-98
16
Mez
23
Mez
21
N
26
N
28
C. Superior 6
C. Sabana 7
CP1: 72.1
CP2: 17.5 ns
Acum : 89.6%
CP1: 84.5
CP2: 8.3 ns
Acum : 92.8%
Cepa
Vc
J-01
5
J-96
8
J-98
17
Mez
20
N
29
J-01
6
J-96
7
J-98
18
Mez
19
N
30
Figura 16. Anlisis biplot GGA del rendimiento de grano de una variedad de soya
inoculada con diferentes cepas de B. japonicum y evaluadas en ocho ambientes de la
Orinoquia colombiana. A. Soyica P-34, B. Orinoquia 3, C. Corpoica Libertad 4, D.
Corpoica Taluma 5, E. Corpoica superior 6 y F. Corpoica Sabana 7.
En general, las variedades presentaron una respuesta positiva a las cepas J-96, mezcla y
J-98 y una respuesta negativa con J-01 y aplicacin de N, con excepcin de C. Libertad 4
que present efectos solo efectos positivos con J-98 y N. Los tratamientos con N mineral
y la utilizacin de J-01, en todas las variedades se ubicaron en el margen izquierdo de la
abscisa, lo que les da una connotacin de respuesta desfavorable ante las variedades. Si
comparamos los rendimientos de grano promedio de las variedades con cada cepa, los
resultados demuestran que los rendimientos de grano fueron superiores con la utilizacin
de cepas efectivas como J-96 y J-98 comparativamente con la cepa J-01 o aplicacin de
150. Kg.ha-1 de N. Estos resultados sugieren que la concentracin de N en el follaje
aunque es un elemento importante de la medicin de FBN, este debe estar acompaado
de otros elementos gentico-ambientales para mximizar el rendimiento de grano de una
variedad.
En la variedad C. Libertad 4, los ambientes se agruparon en dos sectores, separados en
forma clara por semestres A y B. Esto significa que las condiciones ambintales variables
entre
semestres
generaron
sobre
sta
variedad
marcadas
diferencias
en
su
Kaya et al.
3,00
Corpoica
Taluma 5
Rendimiento (t ha-1)
2,50
Corpoica
Superior 6
Corpoica
Sabana 7
Soyica
P-34
Corpoica
Libertad 4
Orinoquia 3
2,00
1,50
1,00
0,50
J96
J01
J98
J01
J98
J96
J01
J98
J96
J01
J98
J96
J01
J98
J96
J01
J98
J96
0,00
Figura 17. Promedios y error estndar del rendimiento de grano (t.ha-1) de 6 variedades
de soya, inoculadas con diferentes cepas de B. japonicum o aplicacin de N mineral.
4.3.3.1 Anlisis GA del modelo AMMI para el comportamiento combinado de
variedades-cepas-N (Vc) y ambientes (A). En el anlisis AMMI, la varianza combinada a
travs de los ambientes mostr diferencias estadsticamente significativas (P<0.01) entre
var-cepa (Vc), entre ambientes (A) y en la interaccin var-cepa x ambiente (VcA). En la
Figura 18 se presenta el biplot en el que se consideran los efectos de los primeros dos
componentes principales (CP), para el anlisis e interpretacin de los patrones de
respuesta de A, Vc e interaccin VcA y el poder discriminativo de los ambientes.
El efecto principal del ambiente (A), explic el 41,95% de la suma de cuadrados de los
tratamientos y el efecto principal de var-cepa (Vc) el 25,03% de la variacin y el 33,06%
restante, fue debido a la interaccin VcA. Estos resultados indican un alto efecto
ambiental sobre el comportamiento de las Vc y que en gran medida, las diferencias en
rendimiento de grano son explicadas por la diversidad de ambientes y su interaccin VcA.
La interaccin VcA fue superior al efecto principal de Vc, lo que demuestra que los Vc
(var-cepa), no responden de manera similar a los cambios ambientales. Los dos primeros
CP, empleados en la construccin del biplot VcA fueron estadsticamente significativos
(P<0.01) y explicaron el 63,64% de la variacin de la interaccin VcA con el 33,5% de los
grados de libertad de la interaccin. El CP1 explic un 47,23% y el CP2 el 16,41%, como
se observa en la Tabla 27.
Tabla 27. Efectos principales aditivos y multiplicativos del anlisis de varianza del
rendimiento de grano (t ha-1) de var-cepas (Vc) a travs de ambientes(A), en oxisoles de
la Orinoquia clombiana.
Fuente Variacin
Var-cepa (VC)
Ambientes (A)
Var-cepa x Amb (VCxA)
CP1
CP2
CP3
CP4
CP5
CP6
CP7
Error medio
C.v (%)
G.L
29
7
203
35
33
31
29
27
25
23
651
20,7
SC
79,193
132,637
104,603
49,4030
17,1704
13,0521
12,1607
7,2900
3,5239
2,0026
107,9385
R2 =0,78
Participacin (%)
25,03
41,92
33,06
47,23
16,41
12,48
11,63
6,97
3,37
1,91
CM
2,7308**
18,9481**
0.5153**
1,41151**
0.52031**
0.42103**
0.41933*
0.27000
0.14096
0.08707
Los coeficientes del CP1 y los promedios de rendimiento de grano de los Vc tuvieron una
alta correlacin (r=0,75; P<0.01), por lo tanto, el CP1 represent significativamente el
rendimiento de los Vc y as, el biplot pudo ser utilizado directamente para identificar los Vc
superiores, como lo sugieren Crossa et al. (2002). Segn Yan et al. (2001), una baja
correlacin entre los coeficientes de CP1 y los promedios de los genotipos puede ser
atribuida a una fuerte y compleja interaccin GA. El rendimiento promedio mostr un
amplio rango entre los Vc, de 1,24 a 2,38 t ha-1. Este rango de promedios de rendimiento
Una localidad discriminante segn Blanche y Myers (2006), es aquella que maximiza la
variacin fenotpica observada entre genotipos para determinada caracterstica. En estas
localidades, es alta la eficiencia y xito de la seleccin de cultivares en programas de
mejoramiento gentico. La habilidad discriminante de una localidad esta comprendida por
una diversidad de factores como tipo de suelo, presencia de patgenos, drenaje de
campo, temperatura, precipitacin, fertilidad del suelo y prcticas de manejo. Algunos de
estos factores son estticos y propios de cada localidad como la acidez de los suelos,
mientras que otros son dinmicos como las enfermedades. En la localidad Lib08a ms
discriminante, se obtuvo el segundo promedio mas alto de rendimiento de grano entre
localidades (2,29 t.ha-1) y el mejor promedio por variedad cepa para C. Superior 6 con la
mezcla de J-01+J-96 (3,64 t.ha-1), aunque sin diferencias con la J-96 y J-98 (3,38 y 3,32
t.ha-1, respectivamente), con la misma variedad. Por lo anterior, el Centro de
Investigaciones La Libertad en evaluaciones de primer semestre sera la condicin ideal
para procesos de seleccin de genotipos superiores.
En la Figura 18, segn Yan et al. (2000) y Yan (2002), la relacin de ambientes est dada
por el ngulo que forman sus vectores con respecto al origen (0,0). Un ngulo menor de
90 entre ambientes permite identificar los sectores donde los genotipos presentan un
comportamiento similar. Los ambientes que forman un ngulo cercano a 90 no guardan
relacin, mientras los que tienen un ngulo cercano a los 180 tienden a ordenar los
genotipos de manera inversa. En el biplot GA se puede apreciar que los ambientes Lib7a
y Lib08a conforman un primer sector, Scruz08a y Tal08a, el segundo sector, Lib7b, Lib8b
y Hor08b, un tercer sector y el ambiente Tal08b puede estar conformando un cuarto
sector. Este agrupamiento en sectores permitira prescindir de uno o ms ambientes que
tienden a jerarquizar de manera similar a las Vc, sin alterar significativamente los
resultados. Por otra parte, los vectores opuestos como los ambientes Lib7a y Lib08a
contra el ambiente Tal08b, o en el caso de Lib7b, Lib8b y Hor08b contra Scruz08a,
jerarquizan los genotipos de manera inversa.
Este anlisis GA permiti separar los semestres A y B del piedemonte llanero (localidad
del C.I La Libertad), donde los semestres discriminan y ordenan de manera diferente los
Vc en las evaluaciones de 2007 y 2008. Esto sugiere un alto efecto de las condiciones
climticas sobre el comportamiento diferencial de las Vc. Por ello, es importante realizar
evaluaciones espacio-temporales antes de generar cualquier recomendacin sobre
variedades.
En general, los tratamientos ms estables o aquellos con CP1 y CP2 cercanos a cero
fueron: Vc 29 y 5 (C. Superior 6 con N y J-01), Vc 16 y 9 (C.Taluma 5 con J-98 y J-96), Vc
1 y 24 (Soyica P-34 con J-01 y mezcla), Vc 7 (C. Sabana7 con J-96). Sin embargo, no
fueron los de mejor rendimiento de grano. Para propsitos de seleccin, se identificaron
los Vc con mayor CP1 (atributo de potencial de rendimiento) y bajo CP2 (atributo de
estabilidad). As, los mejores Vc fueron: 8, 11, 14, 17 y 20, en todos los casos con
rendimientos promedios superiores al promedio general. Estos tratamientos definen el
comportamiento de las variedades C. Superior 6 y Orinoquia 3 con las cepas J-96, J-98 y
mezcla, lo que se puede visualizar con mayor claridad al graficar el rendimiento en funcin
del CP1 (Figura 19).
G.L
SC
29
7
203
35
33
31
29
27
25
23
551
20,7
79,19
132,63
104,60
114,07
19,89
13,83
12,32
12,08
7,25
2,46
0,1658
R2 =0,78
Participacin
(%)
25,0
41,9
33,1
62,1
10,8
7,5
6,7
6,5
3,9
1,3
CM
2,73**
18,94**
0,51**
3,25**
0,60**
0,44**
0,42**
0,44**
0,29*
0,10
Esta interaccin de alta magnitud demuestra que las Vc no responden de manera similar
a los cambios ambientales. Para identificar los Vc superiores y su estabilidad, se elabor
el biplot GGA a partir de los dos primeros componentes principales, considerando G como
el comportamiento combinado de var-cepa (Vc).
Con este modelo, se gener un polgono que uni los puntos de los vectores de los Vc
(var-cepa) ms alejados del origen y se traz una lnea perpendicular por cada lado del
polgono con relacin al origen, ocasionando que los ambientes y Vc fueran separados
en subgrupos, donde los Vc ubicados en el vrtice de cada sector fueron los de mejor
desempeo en los ambientes incluidos, como lo describen Yan et al. (2000) y Yan e Yang
(2002). As, en el biplot se conformaron siete sectores, en cuatro de los cuales quedaron
incluidos todos los ambientes evaluados (Figura 20). Este polgono es la mejor manera de
visualizar el patrn de interaccin entre genotipos y ambientes y para la interpretacin
efectiva del biplot (Yan e Kang, 2003).
El biplot explic el 72,9% del GGA con los CP1 y CP2, donde el CP1 aport el 62,1%.
Los otros componentes presentaron una participacin inferior al 7,6%. Por lo anterior, los
dos primeros dos componentes fueron considerados para la generacin del biplot GGA,
como lo siguieren Dardanelli et al. (2006). Los efectos
principales de Vc (CP1),
presentaron correlacin alta y positiva con rendimiento de grano (r=0,97; P<0.01), muy
superior a la obtenida con el modelo AMMI (r=0,75; P<0.01). Yan et al. (2001), mencionan
que este valor de asociacin debe ser r>0.95 para identificar con certeza los genotipos
superiores. Este tipo de asociacin observado tambin por Yan y Hunt (2001a), Dardanelli
et al. (2006), Crossa et al. (2002) y Yan (2002), permite el uso del biplot para identificar
los Vc superiores y su estabilidad, as como, los ambientes ms discriminantes y de
respuesta similar.
El coeficiente CP2 de Vc y ambientes tomaron valores positivos y negativos; en
consecuencia, segn Yan y Hunt (2001), un genotipo con un alto y positivo CP2 en una
localidad, puede tener un valor grande y negativo en otra localidad, por ello, el CP2
presenta una respuesta desproporcionada para genotipos, el cual puede ser la mayor
fuente de algunas interacciones cruzadas.
- Mejor genotipo por ambiente: Los genotipos que caen sobre el vrtice del polgono en
un sector determinado, presentan los mejores rendimientos en los ambientes de dicho
sector (Yan, 2002; Yan y Hunt, 2002). En la Figura 20, los marcadores de los Vc ms
alejados en el espacio bidimensional (4, 2, 27, 23, 11, 20, 8 y 17), conectados mediante
lneas rectas, configuraron los vrtices del polgono. Estos Vc vrtices son los que ms
aportan a la interaccin, es decir, los de ms alto y ms bajo rendimiento. Los Vc
ubicados en el margen izquierdo del biplot sobre el eje de la abscisa, presentaron baja
adaptacin a los ambientes evaluados, tal como lo menciona Medina et al. (2001). Los del
lado derecho, presentan una relacin positiva con los ambientes considerados, es decir,
J01
J96
3
J98
1
5
66
2
N
Figura 20. Biplot GGA para rendimiento de grano (t.ha-1) de variedades de soya
inoculadas con diferentes cepas de B, japonicum-N, evaluadas en diferentes ambientes
de la Orinoquia.
que tienden a presentar una adaptacin amplia. Estos ltimos Vc se encuentran ubicados
en el ngulo formado por Scruz08a y Hor08b. Los Vc 23, 11, 20, 8 y 17 son los de mejor
comportamiento general en los ambientes evaluados. Los mejores promedios los obtuvo
el Vc 8 y 20 en el ambiente Lib08a, con rendimientos promedios de grano de 3,38 y 3,64
t.ha-1, respectivamente, y que corresponden a la variedad C. Superior 6 con la cepa J-96
y la mezcla. De forma similar, la variedad Orinoquia 3 con la cepa J-96 y mezcla (Vc 11 y
23), se destac por su alto rendimiento promedio general de 2,38 y 2,34 t.ha-1,
respectivamente, con su mejor respuesta en la localidad de Lib07b con 3,33 t.ha-1. Por
otra parte, los Vc 3, 10, 22 y 27, presentan en general bajos rendimientos promedios
entre 1,43 y 1,73 t.ha-1, en los ambientes evaluados, y que representan al comportamiento
de la variedad C. Libertad 4 con la cepa J-01, J-96, mezcla y N. Esta variedad tuvo su
mejor comportamiento en la localidad Hor08B con la cepa J-98 y la aplicacin de N (2,61 y
2,89 t.ha-1, respectivamente). En todos los casos, C. Libertad 4 se ubic en el margen
izquierdo de la abscisa, con una respuesta promedio inferior al promedio.
- Relacin entre ambientes: Los CP1 para ambientes fueron positivos y cercanos a cero.
Tal08b fue el nico ambiente con CP1 negativo pero muy cercano a cero (Figura 20). La
relacin entre ambientes se determin con la metodologa propuesta por Yan (2002),
para un biplot GGA, donde el ngulo formado entre los vectores de dos ambientes esta
relacionado con el coeficiente de correlacin entre ellos. El coseno de este ngulo
corresponde de manera cercana al coeficiente correlacin de los dos ambientes. Los
ngulos inferiores a 90o representan una correlacin positiva entre ambientes, los ngulos
cercanos a 90o pueden tener muy baja o nula relacin con coeficientes de correlacin
cercana a cero y los superiores a 90 presentan correlacin negativa. Esta relacin se
fundamenta en que, las distancias entre ambientes esta dada en grados, donde el coseno
de 90o= 0 (ambientes independientes), coseno de 180o = -1 (ambientes con relacin
opuesta o negativa) y coseno 0 = 1 (ambientes con relacionados positiva)
(www.ggebiplot).
Con estas consideraciones, los ambientes fueron separados de manera notable por
semestre; las localidades del primer semestre Lib07a, Lib08a, Tal08a y Scaruz08a
guardan una estrecha relacin entre si, de la misma manera los ambientes Lib07b,
Lib08b, Tal08b y Hor08b. La mayor correlacin se obtuvo entre las localidades Lib07a y
Lib08a, o Lib07b y Lib08b, mientras que los ms contrastantes fueron Tal08b y Scruz 08a,
como se puede observar en el biplot GGA. Esta observacin coincide con los coeficientes
de correlacin de la Tabla 29, donde las localidades de Lib07a y Lib08a
del primer
semestre tuvieron un r=0,75 (P<0.01), Lib08a con Tal08a r=0,67 (P<0.01), Lib07b con
Hor08b r= 0,45 (P<0.05). Scruza08a present alta correlacin con las dems localidades
del primer semestre, pero no correlacion con ninguna de las localidades del segundo
semestre. No se observaron correlaciones negativas estadsticamente significativas. Con
Tal08b no es posible mantener la relacin entre el biplot y la Tabla 29. Segn Yan, 2002,
se pueden esperar algunas inconsistencias porque el biplot no explica el 100% de la
variacin GGA. Por lo anterior, no sobra utilizar tambin el modelo AMMI. Los resultados
obtenidos con ambos modelos sugieren una diferencia marcada entre los ambientes,
atribuida en gran medida por la poca de siembra, ms relacionada con la disponibilidad
de agua lluvia o balance hdrico entre semestre.
Hor08b
0.074 ns
Lib07a
0.59 **
0.22 ns
Lib07b
0.25 ns
0.45 *
0.33 ns
Lib08a
0.63 **
0.21 ns
0.75 **
0.40 *
Lib08b
0.34 ns
0.38 *
0.64 **
0.45 *
0.56 **
Tal08a
0.60 **
-0.03 ns
0.50 **
0.35 ns
0.67 **
0.59 **
Tal08b
-0.08 ns
0.007 ns
-0.10 ns
-0.12 ns
-0.05 ns
0.38 *
0.34 ns
- Estabilidad fenotipica del rendimiento grano: En le biplot GGA, CP1 representa los
efectos principales del genotipo o comportamiento medio del genotipo (Vc) y el CP2 est
asociado con los efectos de interaccin VcA con cada genotipo (Vc), el cual es
considerado como medida de estabilidad (Yan et al., 2000), Es decir, la estabilidad en el
biplot es interpretada solo por la proyeccin del vector VcA o del vector virtual CP2v. En
el biplot, el efecto principal de genotipos es graficado en las abscisas, para la
interpretacin en trminos de comportamiento medio y estabilidad.
Para crear una coordenada promedia de ambientes y establecer el comportamiento
promedio y estabilidad de los Vc, los coeficientes CP1 y CP2 de los ambientes fueron
promediados para crear el punto medio que al pasar por el origen (0,0) conform el eje X
o CP1 virtual del nuevo biplot (CP1v) y su perpendicular que pasa por el origen (0,0)
conform el nuevo CP2 virtual (CP2v).
existe, estos valores sirven de referencia para la seleccin de cultivares. Los mejores
genotipos sern aquellos que se acerquen al cultivar ideal.
La clasificacin de Vc (var-cepa) extremos por rendimiento y estabilidad se describen en
la Tabla 30. Los Vc con CP1v negativos fueron 1, 5, 6, 15, 25, 28 y 30, asociados con
J-01 y aplicacin de N, y la variedad C. Libertad 4, que por sus bajos rendimientos deben
ser descartados en el proceso de seleccin. En el biplot GGA de rendimiento y estabilidad
(Figura 21), se puede observar con claridad la distribucin de los Vc sobre la abscisa y
ordenada, donde el mejorador puede realizar su proceso de seleccin acorde con sus
intereses. Sobre el margen derecho de la ordenada se ubican las mejores opciones para
la seleccin destacndose por su adaptacin amplia y alto rendimiento los Vc 8, 11 y 20,
y con adaptacin especfica y alto rendimiento los Vc 17 y 23. En contraste, los Vc 3, 10 y
22 obtuvieron los rendimientos promedios ms bajos a pesar de su alta estabilidad, lo
que significa que son malos en todos los ambientes. El comportamiento tpico de un Vc
con baja estabilidad y bajo rendimiento es el 27, por lo que se deduce que es un Vc malo
en algunos sitios.
Tabla 30. Clasificacin de los genotipos (Vc) en funcin del rendimiento de grano
promedio y estabilidad fenotpica.
Estabilidad/Rendimiento
Alta estabilidad
Baja estabilidad
Alto rendimiento
Bajo rendimiento
Adaptacin amplia
Adaptacin especfica
Los modelos AMMI y SREG no permiten establecer las diferencias estadsticas entre
genotipos o tratamientos en los biplot, razn por la cual, mediante anlisis de contrastes,
se seleccionaron las mejores combinaciones de la interaccin entre variedades y cepas
de B. japonicum, para generar la respectiva recomendacin tcnica. Las variedades
Orinoquia 3 y C. Superior 6 con la cepa J-96 fueron las ms estables y rendidoras,
Figura 21. Biplot GGA con vector virtual CP1v y CP2v promedio. Los Vc en el crculo
representan los ms estables y de mayor rendimiento. Las lneas punteadas paralelas al
eje virtual CP2v, representan la proyeccin en estabilidad (ms cercano al eje X, ser ms
estable).
- Ambiente ms discriminante y representativo: El CP1 del biplot GGA estuvo
correlacionado con el potencial de rendimiento y por su cercana al origen del segundo
eje, se le considera el ms estable, como lo sugiere Yan et al. (2000). Por lo tanto, el
sector de mayor inters por discriminar de forma efectiva los Vc a travs de ambientes,
seran aquellos con valores altos en CP1 y ms bajo en CP2. Este principio se cumple
con las localidades Lib08a y Lib07a con un CP1 de 1,16 y 1,13 y CP2 de -0,16 y 0,09,
respectivamente. Sin embargo, segn Yang et al. (2009), los vectores ambintales ms
alejados del origen discriminan mejor los genotipos que aquellos que estn cerca del
origen. En este sentido, el vector con mayor habilidad discriminante de los Vc corresponde
a la localidad de Lib08a, el cual al promediar los ambientes y obtener los ejes virtuales
CP1v y CP2v, present la mayor proyeccin y cercana al eje x. Situacin similar ocurri
con el modelo AMMI. Las localidades con valores PC1 cerca de cero tienen poca
interaccin y baja discriminacin a travs de genotipos (Crossa et al., 1991), como ocurre
con las localidades de Tal08a y Tal08b.
El menor ngulo de los vectores ambientales formado contra el eje virtual CP1v se
present con las localidades Lib07a y Lib08a, que constituyen los ambientes ms
representativos. En este proceso de evaluacin de ambientes, el ambiente ideal ms
discriminante y representativo en promedio fue Lib08a, donde se logr establecer con
mayor claridad el ranqueo de genotipos o Vc. Estos resultados sugieren que en el C.I La
Libertad en siembras de primer semestre es posible obtener los mayores beneficios en un
programa de mejoramiento gentico.
En el biplot GGE (Figura 22), se presenta el comportamiento de diferentes Vc sobre el
ambiente ideal (Lib08a). El eje X para este ambiente es generado por una lnea que
pasa por el origen (0,0) y llega hasta vrtice del ambiente, en este caso a Lib08a y que
conforma el CP1v. Los Vc son ranqueados de acuerdo con la proyeccin sobre este eje
X del ambiente ideal. La lnea perpendicular a este eje x o CP2v, permiti identificar los Vc
con rendimientos inferiores o superiores a la media sobre este ambiente. As, los Vc 11,
20, 8 y 17 se destacaron por su mejor comportamiento en el ambiente Lib08a. Aunque
estos mismos Vc sobresaliron con el modelo AMMI, este modelo no los ranque en forma
clara por su potencial de rendimiento como el SREG.
VC ideales
Vc
ideales
Mejor
ambiente
Figura 22. Biplot GGA con vector virtual CP1v y CP2v sobre el ambiente Lib08a ms
discriminante y representativo. Los Vc ordenados de derecha a izquierda sobre el eje x
representan el ranqueo de mayor a menor por rendimiento de grano sobre este ambiente.
4.3.4. Interpretacin de la interaccin genotipo-ambiente (VcA). El trmino de
interaccin de GA hace referencia en este caso, a la variacin en rendimiento que no
puede ser explicada por el genotipo efecto principal genotpico (G) y el efecto principal
ambiental (A) (Yan y Hunt, 2001b). El conocimiento de magnitud de la interaccin GA
carece de importancia si no se conocen sus causas. Segn Yan y Hunt (2001b), este
conocimiento puede ser utilizado para establecer los objetivos del mejoramiento,
identificar condiciones ideales de evaluacin y formular recomendaciones para reas de
adaptacin ptima de cultivares; adems, puede permitir la manipulacin de algunos
factores ambintales para facilitar la seleccin de los cultivares superiores. Las causas de
la interaccin de variedades de soya y cepas de B. japonicum (Vc) con el ambiente (A),
fueron establecidas mediante la asociacin entre variables agronmicos y los coeficientes
CP1 del efecto principal genotpico (Vc) del anlisis del modelo SREG para rendimiento
de grano, considerando la alta asociacin CP1 y rendimiento de grano de los Vc.
4.3.4.1 Asociacin de PC1 de var-cepa (Vc) y las covariables agronmicas. Una
asociacin positiva y significativa casi perfecta se present entre los coeficientes CP1
(primer componente principal ) y el efecto principal de Vc (rendimiento de grano de los Vc
a travs de ambientes) con r= 0,97 y significancia P<0.01, por ello, estos CP1 pueden ser
interpretados por su efecto principal de Vc, lo que permite deducir que las variables
agronmicas que correlacionaron con estos CP1, tambin estn asociadas e interactan
con la produccin de grano de las variedades a travs de los ambientes en estudio. Las
variables asociadas con valores positivos y significativos estadsticamente fueron: vainas
totales (VT) y dentro de estas, las vainas con 3 granos (V3), la altura de planta (ALM) y el
nmero de nudos (NM), el nmero de ndulos (NN), el peso seco de ndulos (PSN), el
contenido total de nitrgeno (CN) en el tejido areo y el ndice de fijacin biolgica de
nitrgeno (IRFBN) (Tabla 31). Estos resultados demuestran la importancia de la
asociacin simbitica en la produccin de grano de soya y donde, cualquier factor externo
que afecte su normal desempeo tambin afectar el rendimiento de grano y sus
componentes.
Tabla 31. Correlacin del primer componente (CP1) del modelo SREG con diversas
variables agronmicas y de asociacin simbitica.
Caracterstica
Rendimiento de grano
Vainas totales
Vainas de con 1 grano
Vainas con 2 granos
Vainas con 3 granos
Vainas con 4 granos
Vainas vanas
Abortos embrionarios
Peso de grano
Das a floracin
Das a madurez
r
0,97
0,76
0,29
0,23
0,83
0,01
-0,49
0,13
-0,44
-0,16
-0,13
prob
<0,0001
<0,0001
0,1172
0,2275
<0,0001
0,9733
0,0055
0,4785
0,016
0,397
0,4919
Caracterstica
Altura de planta
Nmero de nudos
Nmero de ndulos
Peso seco de ndulos
Peso seco de raz
Peso seco areo
Contenido de N total
ndice relativo IRFBN
Cecospora
Pstula
Roya
r
0,55
0,60
0,74
0,69
0,26
0,35
0,40
0,52
-0,69
-0,10
-0,18
prob
0,0018
0,0005
<0,0001
<0,0001
0,1591
0,0597
0,0278
0,0088
<0,0001
0,5819
0,3502
Los
sentido inverso, es decir, aumento del efecto para unos tratamientos y disminucin para
otros. Con relacin a la presencia de patgenos, la presencia de C. sojina estuvo siempre
asociada con reduccin en rendimientos de grano en los ambientes evaluados y los
factores ambientales favorables, siempre beneficiaron la expresin del rendimiento en
todos los tratamientos en diferente magnitud, lo que corrobora la presencia de la
interaccin cuantitativa.
4.3.4.2. Asociacin de PC1 ambientales y las covariables edafoclimticas. Para
establecer la interaccin GA, es importante descifrar los elementos que ocasionan la
diferenciacin ambiental (Kang, 2002), las cuales pueden ser de tipo bitico o abitico de
carcter esttico (ej. latitud, altitud, textura del suelo, acidez) o dinmico (ej. precipitacin,
fertilidad de suelos, manejo agronmico) y que pueden ser vinculadas al anlisis como
covariables para explicar la interaccin, permitiendo una mayor confiabilidad en la
prediccin del comportamiento de las variedades y cepas. Para establecer la relacin de
las variables edafoclimticas con la interaccin, se asumi el criterio expuesto por Yan
(2002), donde el ngulo formado entre los vectores de dos ambientes en un biplot GGA
esta relacionado con el coeficiente de correlacin entre ellos. De tal manera que dos
vectores ambientales cercanos comparten condiciones ambintales similares. Los
coeficientes de correlacin entre los CP1 de ambientes y las covariables edafoclimticas
se presentan en la Tabla 32.
Tabla 32. Correlacin de CP1 de ambientes para la variable rendimiento de grano con
diferentes covariables edafoclimticas, en oxisoles de la Orinoquia colombiana.
r
prob
Arena
0,43
0,292
Limo
0,21
Arcilla
pH
variable
prob.
BSLl
-0,66
0,151
0,619
RSF
-0,71
0,116
-0,72
0,043
RSLl
-0,56
0,250
-0,38
0,351
RP10
-0,55
0,154
variable
MO
0,40
0,329
RP20
-0,44
0,281
BHF
0,68
0,066
D.a 10
- 0,38
0,357
BHLl
0,92
0,001
D.a 20
- 0,19
0,649
BSF
-0,79
0,061
MO= materia orgnica, BHF=balance hdrico desde siembra hasta floracin, BHLl=balance hdrico
de formacin de vainas a llenado de granos, BS= brillo solar a floracin o a llenado de granos.
RS=radiacin solar a floracin o llenado de grano, RP= resistencia a la penetracin a los 10 y 20
cm de profundidad, D.a.= densidad aparente a 10 y 20 cm.
La precipitacin durante el ciclo del cultivo fue superior en el primer semestre del ao en
todas las localidades (1.238 a 1.566 mm), mientras que en segundo a pesar de
presentarse un buen volumen de agua (634 a 1.122 mm), se presentaron limitantes en la
frecuencia con que ocurri la precipitacin, generndose un dficit hdrico en las etapas
crticas del cultivo como floracin y llenado de grano.
localidad Lib07a, con bajo rendimiento de grano promedio (1,44 t.ha-1), al caer en forma
notoria la precipitacin del 1 al 14 de julio de 2007, donde el balance hdrico fue negativo
en 9 de los 14 das (anexo 1), para un balance negativo total de -7,8 mm., en este perodo
que coincidi con la fase fenolgica de llenado de grano. Este dficit fue acentuado por
una alta resistencia a la penetracin de la races (1,36 Mpa a 20 cm de profundidad).
En la medicin de la resistencia a la penetracin promedia a 10 y 20 cm de profundidad
(RP10 y RP20), aunque no se observaron diferencias estadsticas para los coeficientes de
correlacin, fue notoria la menor resistencia a la penetracin en las localidades Lib07b y
Lib08a con valores de 0,73 y 1,01 Mpa, respectivamente, a 20 cm de profundidad del
suelo, mientras que en las otras localidades, estos valores estuvieron entre 1,36 y 2,31
Mpa. Las dos localidades con la menor RP, presentaron los mayores rendimientos de
grano promedios (2,46 y 2,29 t.ha-1, respectivamente).
Bonel y Montico (2007), encontraron una reduccin en rendimiento asiciada con el
aumento de la resistencia a la penetracin. Ellos encontraron que incrementos de
resistencia de 0,6 Mpa a cerca de 1,2 Mpa, la reduccin en rendimiento podra llegar
hasta el 45% (R2=0,93). Aunque es muy discutido los valores crticos para soya, Beutler et
al. (2005), sugieren que con valores de 0,85 Mpa y un densidad aparente de 1,48 g.cm3
se inicia la reduccin del rendimiento de la soya. La densidad apartente aunque no
present asociacin con rendimiento de grano, sus valores promedios fueron iguales o
superiores a 1,4 g.cm3 en los ambientes de la Oronoquia. Asociado con la compactacin,
est la menor tasa de absorcin de N, segn Faltico et al. (1999). Una situacin similar
se observ en Lib08b donde la resistencia a la penetracin fue de 1,42 Mp a 20 cm de
profundidad y donde dficit hdrico fue marcado del 4 al 12 de noviembre y del 28 de
noviembre al 6 diciembre. Estos factores de dficit hdrico y compactacin de suelos,
fueron los factores que ms afectaron el rendimiento de grano. Segn Rangeon et al.
(2008), la resistencia a la penetracin y la humedad del suelo estn muy relacionadas, a
fertilidad diferencial del suelo y/o las prcticas culturales. Dehghani et al. (2006),
encontraron efectos significativos de la precipitacin, temperatura, estacin de siembra y
humedad relativa sobre la interaccin GA.
Para reducir el impacto de los factores de interaccin expuestos y potenciar el rendimiento
de grano en los oxisoles de la Orinoquia colombiana, se sugiere:
-
beneficios
para
Colombia,
pas
altamente
importador
de
fertilizantes
nmero variable de bandas sobre las variedades de soya, con una amplificacin de
patrones reproducibles sin presencia de bandas monomrficas. El nmero de alelos
detectados por locus estuvo entre 4 y 7, para un total de 32 alelos. Los oligonucletidos
fueron altamente informativos con un contenido polimrfico (CIP) entre 0,69 y 0.85, lo que
gener un promedio de 5,3 bandas por cada cebador. El oligonucletido con mayor CIP
fue Satt 165 de secuencias repetidas tipo (ATT)18, con la generacin de 7 alelos y el
menor para Satt 472, con 4 alelos de la secuencia (ATT)37 (Tabla 33). Grant y Shoemaker
(2009), asociaron estos cebadores con importantes QTLs, as: Satt131 con el QTL NitR7
1-7, responsable de la acumulacin de nitrgeno en estado R5; Satt 472 con los QTLs
CEW 9-2 y Sclero 5-9, para resistencia a la plaga Helicoverpa zea y al patgeno
Sclerotinia sp., respectivamente; Sat_039 asociado con el QTL Sd set 1-2,
para
Orinoquia 3
GMABAB
satt131
Satt472
C. Sabana 7
C. Superior 6
C. Libertad 4
Soyica P-34
C. Taluma 5
0.15
0.22
0.28
0.35
0.42
Figura 23. Grado de similitud entre seis variedades de soya. En el eje X aparece el
coeficiente de similitud de Dice. Se observan las bandas amplificadas con los cebadores
GMABAB, Satt 131 y Satt 472, con carriles ordenados asi: Orinoquia 3, C. Superior, C.
Sabana 7, C. Libertad, Soyica P-34 y C. Taluma 5, respectivamente. En el extremo
izquierdo aparece el marcador de peso 25pb
Los bajos valores en los ndices de similitud demuestran las diferencias marcadas entre
variedades, lo que ratifica la utilizacin de germoplasma divergente en la construccin de
las nuevas variedades de soya para oxisoles de la Orinoquia colombiana. En general, la
especie cultivada Glycine max (L) Merril, posee bajos niveles de diversidad gentica en
comparacin con otras especies autgamas (Pecina, et al., 2005). En este sentido, los
estudios con marcadores moleculares realizados por Brown-Guedira et al. (2000),
diferentes pruebas. Las diferencias moleculares entre cepas se hicieron evidentes con los
dos protocolos de rep-PCR.
El sistema ms complejo se obtuvo con ERIC al generar 16 fragmentos totales en las
cepas de referencia. Resultados similares fueron obtenidos por Ferreira y Hungria (2002),
quienes encontraron un promedio de 15 bandas en B. japonicum. Las bandas tuvieron
tamaos entre 420 y 3500 pb, con la presencia de una banda comn de 1200 bp, entre
cepas. Un mayor bandeo en las cepas de referencia se present en la J-96 y el menor, en
J-01 (Figura 24A). El mayor polimorfismo se observ entre 3000 y 3500 pb y entre 420 y
900 pb.
Con BOX-PCR las bandas presentaron un tamao entre 1900 pb y 3800 pb, en un patrn
molecular de 6 fragmentos. Las cepas J-01 y J-96 compartieron 4 de las 6 bandas,
mientras que J-98 solo comparti una banda de 2500 pb con las otras dos cepas (Figura
24B).
Aunque ERIC present mayor fragmentacin y polimorfismo que BOX, los dos sistemas
permitieron establecer diferencias genticas intraespecficas entre cepas de B. japonicum.
Sin embargo, los patrones moleculares ERIC-PCR obtenidos a partir de ndulos con
respecto a las cepas de referencia del banco de germoplasma fueron diferentes. Las
mayores discrepancias se observaron con la cepa J-96, seguida por la J-01, mientras que
la J-98 fue muy estable. Estos resultados concuerdan con los resultados obtenidos por
Batista et al. (2007), quien encontr que la variabilidad puede ser tal que incluso dentro de
un grupo genticamente ms estable (como la CPAC 15), slo el 6,4% de los aislamientos
mostraron una alta similitud a la cepa inoculante, mientras que ninguno fue similar a la
CPAC 7. La variabilidad gentica parece ser el resultado de una variedad y combinacin
de eventos, incluyendo la dispersin de la cepa, la recombinacin genmica y la
transferencia horizontal de genes. Adems, la variabilidad gentica de las cepas parece
que se asocia principalmente con la adaptacin, capacidad saprofita y competitividad, y no
con la eficacia simbitica, ya que la similitud de los genes simbiticos fue superior a la de
regiones conservadas del ADN.
Estos
protocolos han sido aplicados con xito para caracterizar en forma intraespecfica a B.
japonicum (Judd et al. ,1993; Xue et al., 1999; Ferreira y Hungria, 2002 y Giongo, et al.,
2008).
Estos resultados podran ser tiles para la identificacin de nuevas cepas y formulacin de
nuevos inoculantes con mayor capacidad para fijar N, lo que permitir incrementar la
productividad de la soja en oxisoles de la Orinoquia colombiana. Varios trabajos de
investigacin (Hungria et al., 1996; Santos et al., 1999), reportan la existencia de
diferencias entre cepas parentales y sus variantes, en: produccin de moco, sntesis de
IAA, perfil del factor Nod y el patrn molecular de ADN obtenido con cebadores arbitrarios
o especficos. Godoy et al. (2008), encontraron en estudios genmicos diferencias
genticas notorias entre las cepas CPAC 15 y USDA 110, de las cuales se derivaron las
cepas J-98 e ICA J-01, respectivamente, en aspectos de nodulacin, relacin planta
husped y adaptacin, y genes relacionados con el sistema de secrecin y transporte.
Es importante considerar que, no existen cepas indgenas o nativas en Brasil (Ferreira, y
Hungria, 2002), ni en Colombia, por no ser centros de origen. Sin embargo, la siembras
sucesivas de soyas inoculantes con rizobios introducidos, han generado en estos pases
cepas naturalizadas que han logrado adaptarse a las nuevas condiciones edafoclimticas
y que se constituyen en fuente de diversidad biolgica para procesos de seleccin de
cepas ms efectivas para fijar N.
5. CONCLUSIONES
La asociacin simbitica para fijar N atmosfrico fue diferencial entre cepas de B.
japonicum, variacin atribuida a la especificidad de las cepas por un determinado genotipo
varietal
las
condiciones
ambintales
predominantes.
La
cepa
J-01
fue
6. RECOMENDACIONES
El mejoramiento en la fijacin biolgica de nitrgeno puede contribuir con el incremento
del rendimiento de grano. Por ello, los esfuerzos de investigacin para mejorar la FBN se
deben orientar tanto en la generacin de variedades de alta adaptacin y potencial
gentico como en la seleccin de cepas altamente infectivas y efectivas para la FBN, que
interacten favorablemente para superar la brecha del rendimiento de grano en soya.
Nuevas alternativas genticas estn disponibles en el mercado de variedades de soya
para la Orinoquia colombiana. Estas opciones hacen necesaria la valoracin de la
interaccin especfica con las cepas de B. japonicum existentes para reducir los riesgos
de una baja asociacin simbitica que atente contra la rentabilidad del cultivo.
Aunque las cepas naturalizas en general, son ms competitivas y menos efectivas en el
proceso de la FBN que las cepas introducidas, la bsqueda de cepas naturalizas como
variantes recombinantes, infectivas y efectivas adaptadas a los suelos de sabana debera
ser de alta prioridad para contribuir al desarrollo competitivo y sostenible de la soya en
estos agroecosistemas de la Orinoquia colombiana, hacia donde se estn dirigiendo los
recursos financieros pblicos y privados para la produccin de materias primas de la
cadena avcola-porccola.
7. BILBIOGRAFIA
Abbott L., Pistorale, S. y Filippini, O. 2007 .Anlisis de coeficientes de sendero para el
rendimiento de semillas en Bromus catharticus. Cien. Inv. Agr. 34(2): 141-149.
Acock, M.C. and Acock, B. 1995. Photoperiod sensitivity during soybean flower development.
Biotronics, 24: 25-34.
Agronet. 2009. http:www.agronet.gov.co. Consultado 20 de septiembre de 2009
Akem, C. and Dashiell, K.1994. Effect of planting date on severity of frogeye leaf spot and
grain yield of soybeans. Crop Protec. 13: 607-610.
Alezones, J. y Zocco, J. 2007. Mejoramiento gentico y produccin de semilla de soya.
Fundacin Danac. Gerencia de Investigacin-Gestin de innovacin tecnolgica e
informacin. Yaracuy, Venezuela.
Almansa, E. 2006a. Manejo del recuro hdrico para el cultivo de la soya en la Orinoquia
colombiana. En: Soya (Glycine max (L.) Merril, alternativa para los sistemas de produccin de
la Orinoquia colombiana. Corpoica. Aditora Guadalupe, Ltda. Villavicencio, pp. 135-142.
Almansa, E. 2006b. Requerimientos hdricos del cultivo de la soya en la Altillanura. Plegable
divulgativo No. 50. Corpoica, C.I. La Libertad. Tipografa y Litografa Moderna. Villavicencio.
Amzquita, E., Hoyos, Ph. y Molina, D. 2005. Mejoramiento de la calidad del suelo en
rotacin maz-soya en la altillanura colombiana. Resmenes III Seminario Regional
Agrociencia y Tecnologa Siglo XXI Orinoquia Colombiana. Corpoica, Villavicencio, 45 p.
Ariyo, O. 1995. Correlations and path-coefficients analysis of components of seed yield in
soybeans. African crop science journal, 3(1): 29-33.
Arshad, M., Ali, N., and Ghafoor, A. 2006. Character correlation and path coefficient in
soybean [Glycine max (L.) Merrill]. Pak. J. Botany, 38(1): 121-130.
vila, M., Bracccini, A., Scapim, C., Mandarino, J., Albrechet, L. y Soares, P. 2007.
Componentes do rendimento, teores de isoflavonas, protenas, leo e qualidade de sementes
de soja. Revista brasileira de sementes, 29(3): 1-17.
Ba, A., Samba, R., Sylla, S., Roux, C., Neyra, M., Rousteau, A., Imbert, D., and Toribio, A.
2004. Characterization of the diversity of symbiotic microorganisms in Pterocarpus officinalis
in swamp forests of Guadeloupe and Martinique. Revue de Ecologie-la terre et la vie, 59: 163170.
Baca, B., Soto, L. y Pardo, M. 2000. Fijacin biolgica de nitrgeno. Elementos, 7(38): 43-49.
Ball, R., McNew, R., Vories, E., Keisling T., and Purcell, L. 2001. Path analyses of population
density effects on short-season soybean yield. Agron. J., 93: 187195.
Barcellos, F., Menna, P., da Silva Batista, J., and Hungria, M. 2007. Evidence of horizontal
transfer of symbiotic genes from a Bradyrhizobium japonicum inoculant strain to indigenous
diazotrophs Sinorhizobium (Ensifer) fredii and Bradyrhizobium elkanii in a Brazilian Savannah
soil. Appl. Environ. Microbiol. 73:26352643.
Barrientos, L., Higuera, M., Acua, H., Guerrero, J., Ortega, F. y Seguel I. 2002. Efectividad
simbitica de cepas naturalizadas de Mesorhizobium loti y Bradyrhizobium sp. (Lotus) en
plantas de tres especies del genero Lotus. Agric. Tc. 62(2): 226-236.
Basanta, M., Vega, C., Peirone, L., Alvarez, C. Lovera, E. 2007. Incidencia del sndrome de
tallo verde en soja (Glycine Max L. Merril) bajo siembra directa y labranza convencional en el
centro de la provincia de Crdoba. rea de Agronoma de INTA Manfredi, pp. 62-69.
Bastidas G. y Agudelo O. 1994. El Cultivo de la soya. Manual de asistencia tcnica No.60.
ICA-Corpoica. Palmira Valle del Cauca, Colombia. 460 p.
Batista, J., Hungria, M., Barcellos, F., Ferreira, M., and Mendes, I. 2007. Variability in
Bradyrhizobium japonicum and B. elkanii seven years after introduction of both the exotic
microsymbiont and the soybean host in a Cerrados soil. Microbial Ecol., 53: 270284.
Beutler, N., Centurio, J., and Silva, A. 2005. Soil resistance to penetration and least limiting
water range for soybean yield in a haplustox from Brazil. Braz. arch. biol. technol., 48(6):863871.
Bhardwaj, H.,. Reese, P. Jr., Kang, M., and Zhang M., 2005. Soybean cultivar evaluation via
GGE biplot analysis, J. New Seeds, 7: 37-50.
Bianchi, C., Carter, T., Bailey, Jr, M., Mian, M., Rufty, T., Ashley, D., Boerma, H., Arellano, C.,
Hussey, R., and Parrott W. 2000. Aluminum Tolerance Associated with Quantitative Trait
Loci Derived from Soybean PI 416937 in Hydroponics. Crop Sci., 40:538-545.
Bizeti, H., Carvalho, C., Souza, J., and Destro, D. 2004. Path analysis under multicollinearity
in soybean. Brazilian Archives of Biology and Technology, 47(5): 669-676.
Blanche, S. and Myers G. 2006. Identifying discriminating locations for cultivar selection in
Louisiana. Crop Sci., 46: 946-949.
Board, J., Kang, M., and Harville, B.1999. Path Analyses of the Yield Formation Process for
Late-Planted Soybean. Agron J., 91:128-135.
Board, J., Kang, M., and Harville, B. 1997. Path analyses identify indirect selection criteria for
yield of late-planted soybean. Crop Sci., 37: 879884.
Boddey, L. and Hungria, M. 1997. Phenotypic grouping of Brazilian Bradyrhizobium strains
which nodulate soybean. Biol Fertil Soils, 25: 407415.
Bohlool, B., Ladha, J., Garrity, D., and George, T. 1992. Biological nitrogen fixation for
sustainable agriculture: A perspective Plant and Soil, 141(1-2): 1-11.
Bonel, B. y Montico, S. 2007. Las posibilidades del escarificado en suelos del sur de Santa
66(11): 5099-5103.
CIAT. 1993. Catalogo de Cepas de Rizobio para Leguminosas Forrajeras Tropicales.
Documento de trabajo N 14. 5 Ed. Compilado por: Manuel Franco, Gloria Amparo Ocampo,
Edilfonso Melo y Richard Thomas .Cali, Colombia. 123 p.
Corpoica, 2009. Rutas tecnolgicas y productivas. Llanos Orientales de Colombia. Boletn
informativo.
Crossa J., Fox, P. Pfeiffer, W., Rajaram, S., and Gauch, H. Jr. 1991. AMMI adjustment for
statistical analysis of an international wheat yield trial.Theor. Appl Genet, 81: 27-37.
Crossa, J. and Cornelius, P. 1997, Sites regression and shifted multiplicative model clustering
of cultivar trial sites under heterogeneity of error variances, Crop Sci., 37: 406-415.
Crossa, J., Cornelius, P., and Yan, W. 2002. Biplots of linear-bilinear models for studying
crossover genotype x environment interaction. Crop Sci., 42: 619-633.
Cruz, C. D. 2004. Programa GENES verso windows - Anlises estatstica, biomtrica,
multivariada e diversidade gentica de dados quantitativos e moleculares. Editora UFV,
Universidade Federal de Viosa (Brasil). En: www.ufv.br/dbg/genes/ genes.htm; consulta:
febrero de 2009. verso atual 2007.
Daramola, D., Danso, S., and Hardarson, G. 1994. Nodulation, N2 fixation and dry matter
yield of soybean [Glycine max (L.) merrill] inoculated with effective and ineffective
Bradyrhizobium japonicum strains. Soil Biology and Biochemistry, 26(7):883-889.
Dardanelli, J., Balzarini, M., Martnez, M., Cuniberti, M., Resnik, S., Ramunda, S., Herrero, R.,
and Baigorri, H. 2006. Soybean maturity groups, environments, and their interaction define
mega-environments for seed composition in Argentina. Crop Sci., 46:1939-1947.
Date, R. 2000. Inoculated legumes in cropping systems of the tropics. Field Crops Research,
65(2-3): 123-136.
De la Vega, A. 2000. Interaccin Genotipo x Ambiente y mejoramiento de Girasol para las
regiones Central y Norte de Argentina. Tesis doctoral Facultad de Agronoma. Universidad de
Buenos Aires.
De Sousa Vieira, O. y Milligan, Scott. 2005. Interreltionships of cane yield components and
their utility in sugarcane family selection: path coefficient analysis. Interciencia, 30(2): 93-96.
Dehghani, H., Ebadi, A., and Yousefi, A. 2006. Biplot analysis of genotype by environment
interaction for Barley yield in Iran. Agron. J. 98: 388-393.
Delgado, M., Tresierra, A., Talbi, Ch., Bedmar, E. 2004. Los genes mod de Bradyrhizobium
japonicum codidfican un transportador de molideno de tipo ABC. Resumen VII Reunin
nacional del metabolismo del nitrgeno. Almonte, Huelva.
Ferreira, M. and Hungria, M. 2002. Recovery of soybean inoculant strains from uncropped
soils in Brazil. Field Crops Res. 79: 139152.
Finlay, K. and Wilkinson, G. 1963. The analysis of adaptation in a plant breeding programme.
Aust. J. Agri. Res., 14:742-754.
Fresoli, D. Beret, P. y Rojas, P. 2004. Desarrollo de cultivares de soja mejorados con alto
potencial de rendimiento para la regin litoral. Ciencias Agrcolas y Ciencias Mdicas.
Ciencia, Docencia y Tecnologa, 28: 247-252.
Gabriel, K. 1971. The biplot graphic of matrices with application to principal component
analysis. Biometrics, 58: 453-467.
Gaikwad, S., Bangar N., and Chavan, B. 2007. Correlation and path coefficient analysis in
soybean. J. Maharashtra Agric. Univ., 32: 276277.
Gauch, H. Jr. 2006. Statistical Analysis of Yield Trials by AMMI and GGE. Crop Sci., 46:
14881500.
George, T., Bohlool, B., and Singleton, P. 1987. Bradyrhizobium japonicum-Environment
Interactions: Nodulation and Interstrain Competition in Soils along an Elevational Transect.
American Society for Microbiology, 53(5): 1113-1117.
George,T., Singleton, P., and Ben, B. 1988. Yield, Soil Nitrogen Uptake, and Nitrogen Fixation
by Soybean from Four Maturity Groups Grown at Three Elevations Agron. J., 80: 563-567.
GGEbiplot Pattern Explorer. Complete solution to biplot analysis. Biplot analysis of genotype
by environment data. www.ggebiplot.com.
Giongo, A., Ambrosini, A., Vargas, L., Freire, J., Bodanese, M., Passaglia, L. 2008. Evaluation
of genetic diversity of bradyrhizobia strains nodulating soybean [Glycine max (L.) Merrill]
isolated from South Brazilian fields. Applied soil ecology, 38: 261-269.
Giraud, E., Moulin, L., Vallenet, D., Barbe, V., Cytryn, E., Avarre, J., Jaubert, M., Simon, D.,
Cartieaux, F., Prin, Y., Bena, G., Hannibal, L., Fardoux, J., Kojadinovic, M., Vuillet, L., Lajus,
L., Cruveiller, S., Rouy, Z., Mangenot, S., Segurens, B., Dossat, C., Franck, W., Chang, W.,
Saunders, E., Bruce, D., Richardson, P., Normand, P., Dreyfus, P., Pignol, D., Stacey, G.,
Emerich, D., Vermglio, A., Mdigue, C., and Sadowsky, M. 2007. Legumes Symbioses:
Absence of Nod Genes in Photosynthetic Bradyrhizobia. Science, 316(5829): 1307 1312.
Gizlice, Z., Carter, T. Jr., and Burton, J. 1994. Genetic base for North American public
soybean cultivars released between 1947 and 1988. Crop Sci. 34: 11431151.
Godoy, L, Vasconcelos, A., Chueire, L., Souza, R., Nicolas, M., Barcellos, F., and Hungria, M.
2008. Genomic panorama of Bradyrhizobium japonicum CPAC 15, a commercial inoculant
strain largely established in Brazilian soils and belonging to the same serogroup as USDA
123. Soil Biology & Biochemistry, 40: 2743-2753.
Gollob, H. 1968. A statistical model which combines features of factor analytic and analysis of
variances techniques. Psychometrika, 33: 73-115.
Kaya, Y., Akcura, M., and Taner, S. 2006. GGE-biplot analysis of multi-environment yield
trials in bread wheat. Turk J. Agric. For., 30: 325-337.
Kaya, Y., Palta, ., and Taner, S. 2002. Additive main effects and multiplicative interactions
analysis of yield performances in bread wheat genotypes across environments. Turk J Agric.,
26:275-279.
Keyser, H. and Li, F. 1992. Potential for increasing biological nitrogen fixation in soybean.
Plan and soil, 141(1-2):119-135.
King, C. and Purcell, L. 2001. Soybean Nodule Size and Relationship to Nitrogen Fixation
Response to Water Deficit. Crop Sci., 41:1099-1107.
Kober, M., Saccol de S, E., Freire, J., and Giongo, A. 2004. Characterization of variants of
Bradyrhizobium elkanii and B. japonicum and symbiotic behaviour in soybeans. Cincia Rural,
34(5):1459-1464.
Koutroubas, D., Papakosta, D., and Gagianas, A. 1998. The importance of early dry matter
and nitrogen accumulation in soybean yieldEuropean Journal of Agronomy, 9(1):1-10.
Krishnan, H., Lorio, J., Kim, W., Jiang, G., Kim, K., DeBoer, M., and Pueppke, S. 2003.
Extracellular Proteins Involved in Soybean Cultivar-Specific Nodulation Are Associated with
Pilus-Like Surface Appendages and Exported by a Type III Protein Secretion System in
Sinorhizobium fredii USDA257. MPMI, 16(7):617625.
Kumudini, S., Hume, D., and Chu, D. 2001. Genetic Improvement in Short Season Soybeans:
I. Dry Matter Accumulation, Partitioning, and Leaf Area Duration. Crop Sci., 41:391-398.
Lzaro, L.; Ressia, J., Lett, L., Mendivil, G., Agostini, M. De Pablo, C. y Balbuena, R. 2004.
Efectos de los sistemas de labranza y la inoculacin en soja en siembra tarda. Agriscientia,
21(2):59-66.
Lpez, S, Perticari, A., Piccinetti, C., Ventimiglia, L., Arias, A., De Battista, J., Althabegoiti,
M., Mongiardini, E., Prez, J., Quelas, J., and Lodeiro, A. 2009. In-Furrow Inoculation and
Selection for Higher Motility Enhances the Efficacy of Bradyrhizobium japonicum Nodulation
Agron. J., 101:357-363.
Louws, F., Rademaker, J., and Bruijn, J. 1999. The three Ds of PCR- based genomic analysis
of phytobacteria: diversity, detection, and disease diagnosis. Annu. Rev. Phytopathol, 37:81
125.
Luyindula, N. and Weaver, R. 1989. Nitrogen partitioning in cowpea as influenced by rhizobial
strain and mineral nitrogen. Agronomy Journal, 81:593-596.
Maccheroni, W., Jordao, H., De Gaspari, R., De Moura, G.L., and Matsuoka, S. 2009.
Development of a dependable microsatellite-based fingerprinting system for sugarcane.
Sugar Cane Int. 27(2): 4752.
Malik, M., Qureshi, A., and Ghafoor, A. 2006b. Utilization of Diverse Germplasm for Soybean
Yield Improvement. Asian Journal of Plant Sciences 5(4):663-667.
Malik, M., Ashraf, M., Qureshi, A., and Ghafoor, A. 2007. Assesment of genetic variability,
correlation and path analysis for yield and its components in soybean. Pakistan J. Bot., 39:
405413.
Malik, M., Qureshi, A., Ashraf, M., and Ghafoor, A. 2006a. Genetic variability of the main yield
related characters in soybean. Int. J. Agric. Biol., 8(6):815819.
McGinnis, L. y Suszkiw, J. 2007. En busca de una soya resistente a las enfermedades.
Oleaginosas, Comit Nacional Sistema- Producto.
http://www.oleaginosas.org/impr_144.shtml.
Medina, M, Abadie T., Vilar, D. y Ceretta, S. 2001. Estudio metodolgico de adaptacin de
cultivares de maz para silo a las condiciones de Uruguay. Agrociancia 5(1):23-31.
Minamisawa, K., Seki, T., Onodera, S., Kubota, M., and Asami, T. 1992. Genetic relatedness
of Bradyrhizobium japonicum field isolates as revealed by repeated sequences and various
other characteristics. Appl Environ Microbiol, 58:28322839.
MADR-Ministerio De Agricultura Y Desarrollo Rural. 2008. Estadsticas del sector
agropecuario. Principales indicadores agropecuarios.
Mishra, R., Baenziger, P., Russell, W.Graybosch, R., Baltensperger, D., and Eskridge, K.
2006. Crossover Interactions for Grain Yield in Multienvironmental Trials of Winter Wheat ,
Crop Sci., 46:1291-1298
Mishra, A., Ali, S., Tiwari R., and Raghuwanshi, R. 1994. Correlation and path analysis in
segregating populations of soybean. International journal of tropical agriculture, 12(3-4): 278281.
Mitchell, R. 1992. Testing evolutionary and ecological hypotheses using path analysis and
structural equation modeling. Funct. Ecol., 6:123-129.
Molano, J. 1998. Biogeografa de la Orinoquia colombiana. En: Colombia Orinoco. Fondo
para la Proteccin del Medio Ambiente Jos Celestino Mutis FEN-Colombia. Biblioteca
virtual Luis ngel Arango.
Molina, D., Amzquita, E. y Hoyos, Ph. 2003. Construccin de capas arables en suelos
oxisoles de la Altillanura colombiana. Evaluacin de parmetros y procesos hidrolgicos en el
suelo. VII Escuela Latinoamericana de Fsica de Suelos. La Serena, Chile. Pp. 113-117.
Mumns, D. 1987. Nitrogen fixation potential of bean (P. vulgaris) compared with other grains
legumes under controlled conditions. Plant and Soil, 98(2):169-182.
Njera, G. y Moreno, L. 2004. Comportamiento de familias S1 de maz en distintos pH del
Perticari, A., Arias, N., Baigorri, H., De Battista, J., Lett, L., Montecchia, M., Pacheco, J.,
Simonella, A., Toresani, S., Ventimiglia, L. y Vicentini, R. 2003. Inoculacin y fijacin biolgica
de nitrgeno en el cultivo de soja. En: El libro de la soja. Ed. SEMA. Captulo 7, pp. 69-76.
Pietrarelli, L., Zamar, J., Legua, H., Alessandria, E., Snchez, J., Arborno, M., and Luque, S.
2008. Effects of different management practices on nodulation and on the yield of soybean
crops. Agriscientia, 25(2):81-87.
Ppolo, A.E. 2002. Influncia da temperatura sobre as concentraes de protena e leo em
sementes de soja (Glycine max (L.) Merrill).. Tesis de Doctorado. Universidade de So Paulo,
Escola Superior de Agricultura Luiz de Queiroz. Piracicaba, So Paulo, Brasil. 128 p.
Ponce, M., Ortiz, R., De la F, C y Moya, C. 2002. Estudio comparativo de nuevas variedades
de soya (Glycine max L. Merrill) para las condiciones de primavera en Cuba. Cultivos
tropicales, 23 (2):55-58.
Portafolio. 2008. Negocios. Altos precios internacionales de insumos comenzaron a sentirse
en Colombia. Portafolio.com.co.
Prithiviraj, B., Souleimanov, A., Zhou, X., and Smith, D. 2000. Differential response of
soybean (Glycine max (L.) Merr.) genotypes to lipo-chito-oligosaccharide Nod Bj V (C(18:1)
MeFuc). J. Exp. Bot., 51:2045-2051.
Pueppke, S., Bolaos, M., Werner, D., Bec., M., Prom, J., and Krishnan, H. 1998. Release
of Flavonoids by the Soybean Cultivars McCall and Peking and Their Perception as Signals
by the Nitrogen-Fixing Symbiont Sinorhizobium fredii. Plant Physiol., 117: 599-606.
Purcell, L. 1999. Physiology, nutrition and fertilization of soybeans in the United States.
Anales Simposio "Monitoramiento nutricional para la recomendacin de fertilizacin de
cultivos". POTAFOS. Piracicaba, SP, Brasil.
Purseglove, J. 1984. Tropical Crops. Dicotyledons. Ed. Longman Group Ltd. pp. 265 271.
R program ver. 2.7.0. 2008. The R Foundation for Statistical Computing.
Racca, R. y Pollino, D. 2005. Bases fisiolgicas para el manejo de la fijacin biolgica del
nitrgeno en soja. Actas del Congreso Mundo Soja, Buenos Aires (Argentina), pp. 111-120.
Raij, B., da Silva, N. Bataglia, O., Quaggio, J., Hiroce, R., Cantarella, H., Belinazzi, R. Jr.,
Dechen, A. y Trani, P. 1985. Recomendaes de adubao e calagem para o Estado de So
Paulo. Bol. Tec. 100. Campinas, SP, Brasil.
Rakhi Agrawal, Shweta Bajoria, S., and Pareek, P. 2006. DNA Isolation from Rhizobium by
Phenol Chloroform Method. Birla Institute of Scientific Research, Jaipur (IND). Protocol
Online. http://www.protocol-online.org
Rajanna, M., Viswanatha, S., Kulkarni R., and Ramesh, S. 2000. Correlation and path
analysis in soybean (Glycine max (L.) Merrill). Crop Research Hisar, 20(2):244-247.
Rangeon, N., Aciar, L., Osinaga, R., Arzeno, J. y Snchez, C. 2008. Anlisis de la resistencia
a la penetracin y humedad de suelo como indicadores de calidad en distintos sistemas de
labranza. XXI Congreso argentino de ciencia del suelo. Potrero de los Funes (SL), Argentina.
Rengel, Z. 2002. Breeding for better simbiosis. Plant and Soil, 245:147162.
Ridner, E. 2006. Soya, propiedades nutricionales y su impacto en la salud. 1ra. ed. - Buenos
Aires. Grupo Q S.A.: Sociedad Argentina de Nutricin. 99 p.
Rivas, L., Hoyos Ph., Amezquita, E. y Molina, D. 2004. Manejo y usos de los suelos de la
altillanura colombiana. Proyecto de evaluacin de impacto, convenio MADR-Ciat.
Roche, P., Maillet, F., Plazanet, C., Debelle, F., Ferro, M., Truchet, G., Prome, J. C., and De
Narie, J. 1996. The common nodABC genes of Rhizobium meliloti are host-range
determinants. Proc Natl Acad Sci USA, 93:1530515310.
Roehrig, S. 1996.
16:55-66.
Rohlf, F. 1997. Software NTSYS-pc. Numerical Taxonomy and Multivariate Analysis System
versin 2.02.
Ruiz, T., Imperial, J., and Palacios, J. 2000. In Prokaryotic Nitrogen Fixation: A Model System
for Analysis of a Biological Process (Triplett, E.W., ed.), Horizon Scientific Press,
Wymondham, U.K. pp. 489507.
Ruz, R., Santos, E. y Torres, J. 2008. Influencia del bradyrhizobium japonicum y la
fertilizacion nitrogenada en los rendimientos de la soya (Glycine max (L.) Merril) desarrollada
en un suelo fersialitico pardo rojizo tpico. p. 7 Biblioteca virtual de ciencias en Cuba.
Sadowsky M., Raymond, T., Perry, C., and Keyser H. 1987. Genetic Diversity in
Bradyrhizobium japonicum Serogroup 123 and Its Relation to Genotype-Specific Nodulation
of Soybean. Applied And Environmental Microbiology, 53(11):2624-2630.
Sadowsky, M., Cregant, P., Gottfert, M., Sharma,A., Gerhold,D., Rodriguez, F., Keyserii,
H.,Hennecke, H., and Stacey, G. 1991. The Bradyrhizobium japonicum nolA gene and its
involvement in the genotype-specific nodulation of soybeans. Proc. Nati. Acad. Sci. 88: 637641.
Saghai, M., Soliman, K., Jorgensen, R., and Allard, R. 1984 Ribosomal DNA spacer-length
polymorphisms in barley: Mendelian inheritance, chromosomal location, and population
dynamics. Proc. Natl. Acad. Sci., 81:8014-8019.
Salamanca, C. y Ramrez, M. 2000. ICA J 01 Biofertilizante para soya. Inoculante de rizobio
para los Llanos Orientales. Plegable divulgativo N 16 CORPOICA-Ministerio de agricultura.
FOCC. Villavicencio, Meta, Colombia.
Shukla, S., Singh, K., and Pushpendra. 1998. Correlation and path coefficient analysis of yield
and its components in soybean (Glycine max L. Merrill). Soybean genetics newsletter (USA),
25:67-70.
Silvester, R. 1983. Fijacin biolgica de nitrgeno por leguminosas: aspectos agronmicos
relacionados con su nodulacin con Rhizobium. Suelos Ecuatoriales, 13(2):28- 35.
Silvester, R., Kipe, J., Harris, D. 1987. Simbiosis leguminosa-rizobio: Evaluacin, seleccin y
manejo. Centro Internacional de Agricultura Tropical. Cali, Colombia. 72 p.
Simanungkalit, R., Roughley, R., Hastuti, R., Indrasumunar, A., and Pratiwi, E. 1996.
Inoculation of soybean with selected strains of Bradyrhizobium japonicum can increase yield
on acid soils in Indonesia. Soil Biology and Biochemistry, 28(2):257-259.
Sinclair, T., Purcell, L., Vadez, V., Serraj, R., King, C., and Nelson, R. 2000. Identification of
Soybean Genotypes with N2 Fixation Tolerance to Water Deficits Crop Sci., 40(6):1803
1809.
Spaink, H. 2000. Root nodulation and infection factors produced by rhizobial bacteria. Ann.
Rev. Microbiol., 54:257-288.
Sprent, J. 1980. Root nodule anatomy, type of export product and evolutionary origen of so
me leguminosae. Plant Cell Environ., 3:35-43.
Stacey, G. 1995. Bradyrhizobium japonicum nodulation genetics. FEMS Microbiology Letters,
127(1-2):1-9.
Steel, R. and Torrie, J. 1989. Bioestadstica, principios y procedimientos. Segunda edicin
(primera en espaol). McGraw-Hill. Mxico. 622 p.
Taware, S., Halvankar, G., Raut, V., and Patil, V. 1997. Variability, correlation and path
analysis in soybean hybrids. Soybean genetics newsletter (USA), 24:96-98.
Thies, J. 1990. Modeling ecological determinants of the symbiotic performance of introduced
rhizobia in tropical soils. Tesis Ph.D. University Of Hawaii. 166 p.
Tiwari, S. and Bhatnagar, P. 1991. Pod shattering as related to other agronomic attributes in
soybean. Tropical agriculture, 68(1):102-103.
Tsvetkova, G., Teofilova, T., and Georgiev, G. 2006. Effect of naringein and quercetin on
Activity of nodabc genes of strain d293 and Following nodulation and nitrogen fixation
Response of inoculated pea plants (pisum sativum L.). Gen. Appl. Plant Physiology, Special
Issue, pp. 67-71.
Tsyganov, V., Morzhina, E., Stefanov, S., Borisov, A., Lebsky, V., and Tikhonovich, I. 1998.
The pea (Pisum sativum L.) genes sym33 and sym40 control infection thread formation and
root nodule function. Mol. Gen. Genet., 295:491503.
UPOV. 2007. Proyecto directrices para los perfiles de ADN: Seleccin de marcadores
moleculares
y
Creacin
de
una
base
de
datos
(directrices
BMY).
http://www.upov.int/es/documents/c/41/bmt_guidelines_proj_9.pdf
USDA. 2007. Soybeans and Oil Crops. Soybean Projections, 2007-16.
http://www.ers.usda.gov/Briefing/SoybeansOilCrops/2007baseline.htm.
Vincent, J. 1970. A Manual for the Practical Study of Root Nodule Bacteria. Blackwell
Scientific, Public. Oxford.UK
Valencia, R. and Leal, D. 2004. Genetic Alternatives form Production System in the Acid-Soil
Savannas of the Colombian Orinoquia. Agropastoral System for the Tropical Savannas of
Latin America. pp. 127-140.
Valencia, R., Carmen, H., Vargas, H. y Arrieta, G. 2006. Variedades mejoradas de soya para
zonas productoras actuales y potenciales de Colombia. Innovacin y cambio tecnolgicoCORPOICA. 4(2-3):7-15.
Valencia, R., Pardo, O., Rincn, A. y Bueno, G. 2006b. CORPOICA TALUAM-5. Variedad
de soya de uso forrajero para sistemas de produccin de la Orinoquia colombiana. Boletn
tcnico No. 48. NMC Publicidad e impresos. Villavicencio-Meta. 32 p.
Valencia, R. 2009. Semilla Corpoica Sueprior 6, variedad de soya. Oferta gentica para la
orinoquia colombiana. Plegable divulgativo No. 45. Publicacin Corpoica, C.I. La Libertad,
Villavicencio, Meta.
Vallejo, F. y Estrada, E. 2002. Mejoramiento gentico de plantas. Universidad Nacional de
Colombia. Impresin Feriva. S.A. pp. 183-201.
Vargas, M. and Crossa, J. 2000. The AMMI analysis and the graph of the Biplot in SAS.
CIMMYT, Int. Mxico. 42 pp. (www.cimmyt.org/biometrics).
Vencovsky, R. y Barriga, P. 1992. Gentica Biomtrica no Fitomelhorarnento de Plantas.
Sociedade Brasileira de Gentica. Riberao Preto, Sao Paulo, Brasil. 486p.
Vera, R. 2000. Investigacin en sistemas agropastoriles: antecedentes y estrategias. En:
Sistemas agropastoriles en sabanas tropicales de america latina. Ciat-Embrapa. pp. 1-6.
Versalovic, J., Schneider, M., de Bruijn, F., Lupski, J. 1994. Genomic fingerprinting of bacteria
using repetitive sequence-based polymerase chain reaction. Methods Mol. Cell Biol. 5: 2540.
Villalobos E., Barrantes, R. y Echandi, R. 1985. Asociacin de algunas caractersticas
cuantitativas con el rendimiento de la soya (Glycine max) en el trpico mediante dos tcnicas
de regresin. Agron., Castano, 9 (1):65-70.
Villar, J., Astegiano, E. y Cencig, G. 2006. Sndrome de tallo verde (stv) en el centro de la
provincia de santa fe: su asociacin con factores ambientales. Campaa 2004/05. In
Resmenes ACSOJA, ed. Tercer Congreso de Soja del Mercosur, Rosario, Argentina. pp. 1720.
Von, H. and Mutert, E. 1995. Global extent, development and economicimpact of acid soils.
Plant Soil, 171:1-15.
Wang, H., Man, Ch., Wang, E., and Chen, W. 2009. Diversity of rhizobia and interactions
among the host legumes and rhizobial genotypes in an agricultural-forestry ecosystem. Plant
Soil, 314:169182.
Wang, E., Martnez, J. y Lpez, I. 2001. Rhizobium y su destacada simbiosis con plantas.
En: Microbios. Capitulo 8. www.microbiologia.org.mx/microbiosenlinea.
Wang, H., Man Chao Xin, Wang En Tao, and Chen Wen Xin. 2009. Diversity of rhizobia and
interactions among the host legumes and rhizobial genotypes in an agricultural-forestry
Ecosystem. Plant Soil, 314:169182.
Whittaker, A., Vazzana, C., Vecchio, V., and Benedettelli, S. 2009. Evaluation of direct and
indirect effects of flavonoids, mineral elements and dry weight on antiradical scavenging
activity in leaf material of field-grown Trifolium pratense cultivars using Path Analysis. Field
Crops Research, 113(1):1-11.
Wright, S. 1921. Correlation and causation. Journal of Agriculture Research, 20:557-585.
Xue, Z., Muilenburg, M., Lie, T., Jun, Z., and Fu, L. 1999. Characterization of soybean
rhizobia at different levels using PCR based techniques. Biodiversity Science, 7(4): 277-284.
Yan, W. and Hunt, L. 2002, Biplot analysis of diallel data. Crop Sci., 42: 21-30.
Yan, W., Hunt, L., Sheng, Q., and Szlavnics, Z. 2000. Cultivar evaluation and megaenvironment investigation based on the GGE biplot. Crop Sci., 40:597-605.
Yan, W. 2002. Singular-value partitioning in biplot analysis of multienvironment trial data.
Agron. J., 94:990-996.
Yan, W. and Kang, M. 2003. GGE biplot analysis: A graphical tool for breeders, geneticists,
and agronomists. CRC Press, Boca Raton, FL.
Yan, W. and Ma, B. 2006. Model diagnosis and GGE biplot analysis, Eastern Cereal and
Oilseed Research Centre, Agriculture and Agri-food Canada. 39 p.
Yan, W. and Rajcan, I. 2002. Biplot evaluation of test sites and trait relations of soybean in
Ontario, Crop Sci., 42:11-20.
Yan, W. and Tinker, N. 2005. A biplot approach to the investigation of QTL-by environment
patterns, Molecular Breeding, 15:31-43.
Yan, W. and Tinker, N. 2005. An integrated biplot analysis system for displaying, interpreting,
and exploring genotype X environment interaction. Crop Sci., 45:1004-1016.
Yan, W. and. Kang. M. 2002. GGE biplot analysis: a graphical tool for breeders, Geneticists,
and Agronomists. CRC Press. 271 p.
Yan, W., Cornelius, P., Crossa, J., and Hunt, L. 2001. Two types of GGE biplots for analyzing
multi-environment trial data. Crop Sci., 41:656663.
Yan, W. and Hunt L. 2001. Interpretation of Genotype x Environment Interaction for Winter
Wheat Yield in Ontario. Crop Sci., 41:19-25.
Yang, R., Crossa, J., Cornelius, P., and Burgueo, J. 2009. Biplot Analysis of Genotype
Environment Interaction: Proceed with Caution. Crop Sci., 49:15641576.
Yang, R. 2007. Mixed-Model Analysis of Crossover GenotypeEnvironment Interactions, Crop
Sci., 47:10511062.
Yin, X. and Vyn, T. 2005. Relationships of Isoflavone, Oil, and Protein in Seed with Yield of
Soybean. Agron. J., 97(5):1314 1321.
Zdor, R. and Pueppke, S. 1990. Competition for nodulation of soybean by Bradyrhizobium
japonicum 123 & 138 in soil containing indigenous rhizobia. Soil Biol. Biochem, 22:607-613
Zhang, H., Charles, T., Drisco, B., Prithiviraj, B., and Smith, D. 2002. Low Temperature
Tolerant Bradyrhizobium japonicum Strains Allowing Improved Soybean Yield in ShortSeason Areas. Agronomy Journal 94:870-875.
Zhang, L., Wang, R., and Hesketh, J. 2001. Effects of Photoperiod on Growth and
Development of Soybean Floral Bud in Different Maturity. Agronomy Journal, 93:944-948.
Zobel, R., Wright, M., and Gauch, H. 1988. Statistical anlisis of a yield trial. Agron J., 80:388393.
8. ANEXOS
40
Floracin
Siembra
Precipitacin (mm)
ETP (mm)
Cosecha
60
Cosecha
80
100
Floracin
120
Siembra
140
20
M M M M M M M M J J J J J J J Jl Jl Jl Jl Jl Jl Jl Jl A
10
14
18
22
26
30
4
8
12
16
20
24
28
1
5
9
13
17
21
25
29
3
7
11
15
3
7
11
15
19
23
27
31
4
8
12
16
20
24
28
2
6
10
14
18
22
26
30
3
7
0
S S S S S S O O O O O O O N N N N N N N N D D D D
Anexo 1. Precipitacin y ETP acumulada (mm), estacin La Libertad 2007, localidades Lib07A y Lib07B.
60
Floracin
Floracin
Siembra
Cosecha
80
100
Precipitacin (mm)
ETP (mm)
Cosecha
120
140
Siembra
160
40
20
M M M M M M J J J J J J J Jl Jl Jl Jl Jl Jl Jl Jl A A A
21
25
29
3
7
11
15
19
23
27
31
4
8
12
16
20
24
28
2
6
10
14
18
22
0
11
15
19
23
27
31
4
8
12
16
20
24
28
2
6
10
14
18
22
26
30
3
7
11
15
180
S S SOOOOOOOO N N N N NN N D D D D D D
Anexo 2. Precipitacin y ETP acumulada (mm), estacin La Libertad 2008, localidades Lib08A y Lib08B.
Cosecha
Floracin
Siembra
60
80
ETP (mm)
Cosecha
100
120
Floracin
140
Siembra
Precipitacin (mm)
160
40
20
17
21
25
29
2
6
10
14
18
22
26
30
4
8
12
16
20
24
28
1
5
9
13
17
21
13
17
21
25
29
3
7
11
15
19
23
27
31
4
8
12
16
20
24
28
2
6
10
14
18
0
M M M M J J J J J J J J Jl Jl Jl Jl Jl Jl Jl A A A A A
S S S S S O O O O O O O O N N N N N N N D D D D D
Anexo 3. Precipitacin y ETP acumulada (mm), estacin Las Margaritas 2008, localidades Tal08A y Tal08B.
Cosecha
60
Floracin
Siembra
80
Cosecha
100
120
Siembra
140
Floracin
Precipitacin (mm)
ETP (mm) Taluma
40
20
M M M M J J J J J J J J J J J J J J J J A A A A A
6
10
14
18
22
26
30
4
8
12
16
20
24
28
1
5
9
13
17
21
25
29
3
7
11
16
20
24
28
1
5
9
13
17
21
25
29
3
7
11
15
19
23
27
31
4
8
12
16
20
24
0
S S S S S S S O O O O O O O N N N N N N N N D D D
Anexo 4. Precipitacin y ETP acumulada (mm), estacin La Palomera 2008, localidades Scruz08A y Hor08B.