Journal of Heredity Advance Access originally published online on July 13, 2005
Journal of Heredity 2005 96(5):586-592; doi:10.1093/jhered/esi084
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Brief Communication |
Epistasis in the Expression of Relevant Traits in Cassava (Manihot esculenta Crantz) for Subhumid Conditions
From the Hue University of Agriculture and Forestry, 24 Phuong Hung Street, Hue City, Vietnam (Cach); International Center for Tropical Agriculture (CIAT), Apartado Aéreo 6713, Cali, Colombia (Lenis, Perez, Calle, Morante, and Ceballos); and Universidad Nacional de Colombia, Carrera 32, Chapinero vía Candelaria, Palmira, Colombia (Ceballos)
Address correspondence to H. Ceballos at CIAT at the address above, or e-mail: h.ceballos{at}cgiar.org.
| Abstract |
|---|
|
|
|---|
There is limited knowledge on the inheritance of agronomic traits in cassava and the importance of epistasis for most crops. A nine-parent diallel study was conducted in subhumid environments. Thirty clones were obtained from each F1 cross. Each clone was represented by six plants, which were distributed in three replications at two locations. Therefore the same 30 genotypes of each F1 cross were planted in the three replications at the two locations. Analysis of variance suggested significant genetic effects for all variables analyzed (reaction to thrips, fresh root and foliage yields, harvest index, dry matter content, and root dry matter yield). Significant epistatic effects were observed for all variables, except harvest index. Dominance variance was always significant, except for dry matter content and dry matter yield. Additive variance was significant only for reaction to thrips. Results suggested that dominance plays an important role in complex traits such as root yield. The significance of epistasis can help us understand the difficulties of quantitative genetics models and QTLs in satisfactorily explaining phenotypic variation in traits with complex inheritance. Significant epistasis would justify the production of inbred parental lines to fix favorable allele combinations in the production of hybrid cassava cultivars.
| Introduction |
|---|
|
|
|---|
Cassava (Manihot esculenta Crantz), along with maize, sugarcane, and rice constitute the most important sources of energy in the diet of most tropical countries of the world. Cassava is the fourth most important basic food after rice, wheat, and maize and is a fundamental component in the diet of millions of people (FAO/FIDA 2000). Scott et al. (2000) estimated that for the 199597 period, annual production of cassava was about 165.3 million tons, with a value of
$8.8 billion (US$). Little progress in understanding the inheritance of agronomic traits in cassava has been achieved. Few articles regarding the inheritance of quantitative traits have been published (Easwari et al. 1995; Easwari and Sheela 1998; Losada 1990). Cassava is perhaps unique in that a molecular map has been already developed (Cortes et al. 2002; Fregene et al. 1997; Jorge et al. 2000, 2001; Mba et al. 2001; Okogbenin and Fregene 2003), but it is complemented with limited traditional genetics knowledge. Cassava is also an interesting crop because its vegetative propagation allows the estimation of within-family genetic variation and, indirectly, the relative importance of epistatic effects. Genetic studies analyzing the importance of epistatic effects are not very common, particularly in annual crops.
Accurate measurement of epistatic effects for complex traits, such as yield, is difficult and expensive. Reports in the literature on the relevance of epistasis are not as frequent as those estimating additive and dominance variances or effects and generally take advantage of the vegetative multiplication that some species offer (Comstock et al. 1958; Foster and Shaw 1988; Isik et al. 2003; Rönnberg-Wästljung and Gullberg 1999; Rönnberg-Wästljung et al. 1994; Stonecypher and McCullough 1986). In many cases, these reports are on forest trees. Because of the complexities of these analyses and the costs involved, reports in the literature related to epistatic effects are frequently based on a limited number of genotypes.
Holland (2001) published a comprehensive review on epistasis and plant breeding. Several cases of significant epistasis have been reported in self-pollinated (Brim and Cockerham 1961; Busch et al. 1974; Gravois 1994; Hanson and Weber 1961; Pixley and Frey 1991; Orf et al. 1999) and cross-pollinated (Ceballos et al. 1998; Eta-Ndu and Openshaw 1999; Lamkey et al. 1995; Melchinger et al. 1986; Wolf and Hallauer 1997) crops. According to Holland (2001) finding significant epistasis seems to be easier in self- than in cross-pollinated species and in designs based in the contrasts of means rather than the analysis of variances.
The objective of this study was to analyze the within-family variation in a diallel study conducted in two subhumid environments and assess the relative importance of additive, dominance, and epistatic genetic effects on the expression of several relevant traits of cassava.
| Materials and Methods |
|---|
|
|
|---|
A diallel mating design was used to generate F1 crosses among nine parents. Inbreeding level of parental lines was considered zero because no self-pollination has been involved in cassava breeding, and crosses among related clones are generally avoided. Kawano et al. (1978) provided evidence that cassava is a highly heterozygous species. Controlled pollinations were performed following the standard procedures described by Kawano (1980). Many parental clones were initially involved, but the parents ultimately used (as well as the number of parents involved) were those that allowed for as balanced a set of crosses as possible. Botanical seed were germinated and grown in a screen house until the seedlings were 2 months old, when they were transplanted to the field at CIAT experimental station in Palmira, Valle del Cauca, Colombia. F1 plants were grown in the field for 10 months. Among the many genotypes (>30) from a given F1 cross, 30 were randomly chosen for this study based solely on their capacity to produce at least six vegetative cuttings. Each of these stakes was planted in one of three replications at one of two locations.
Trials were planted during July 2001 in two subhumid locations in Colombia (Cach et al. forthcoming). A randomized complete block design was used. The evaluation was similar to a split-plot design. Each replication contained 36 main plots, one for each of the 36 F1 crosses of the diallel. Each F1 cross was therefore randomly allocated within each replication. Main plots contained eight rows with seven plants per row. The first and last rows and the first and last plant within each row were filled with border plants. The rest of the plot (6 x 5 = 30 subplots) was used to plant the experimental material. The 30 clones constituting each F1 cross were planted together in the respective main plots of each replication. The experimental design, therefore, offered two types of error: (1) associated with the main plots or F1 averages, and (2) the error associated with the subplots or within-F1 variation. Row-to-row distances and separation of plants within row were 1 m for a final plant density of 10,000 plants ha1.
The six vegetative cuttings obtained from each plant in the nursery at Palmira were distributed in the three replications at the two locations for the evaluation trials. Therefore for each F1 cross, the same group of 30 genotypes was used in each experimental plot. Trials were harvested in May 2002, 10 months after planting (the usual age for harvesting cassava in this environment). One month after planting, 330 kg ha1 of a 15-15-15 NPK fertilizer was applied to the soil, following the standard recommendations for cassava grown in this region of Colombia.
Plants were hand harvested individually. The roots produced by each plant were weighed as well as the above-ground biomass (stem and foliage). Harvest index was measured as the ratio between root weight and total biomass. Root dry matter content was estimated using the specific gravity methodology (Kawano et al. 1987). Approximately 35 kg of roots were weighed in a hanging scale (WA) and then, the same sample was weighed with the roots submerged in water (WW). Dry matter content of the roots produced from each plant was estimated individually utilizing the following formula:
![]() |
Reaction to thrips (Frankliniella williamsi), plant type architecture, and general root appearance were scored using a 15 scale where 1 = resistant or excellent and 5 = susceptible or very poor (CIAT 2002). Plant type score took into consideration several important characteristics, such as plant vigor, erect architecture with few branches and reduced branching angle, adequate capacity to produce vegetative cuttings, amount of foliage present at harvest time, and absence of foliar diseases (which in this particular environment are not frequent).
Statistical Model
The analysis of variance was conducted following the expectations for each mean square described in Table 1. The analysis takes advantage of the full- (FS) and half-sib (HS) families that the diallel mating design creates. As is commonly the case, a few plants died or failed to develop normally to be harvested. Therefore in a few F1 crosses, fewer that 30 clones were actually evaluated in the field in each of the three replications at the two locations. To take into consideration this lack of uniformity, the harmonic (not the arithmetic) mean was used as k in the expected mean squares formulas (Vencovsky and Barriga 1992; see bottom of Table 2). The total genetic variance was partitioned into between-family variation (
) and the within-family variation (
). The between-family variation, in turn, was partitioned into the well-known variances related to general (
) and specific (
) combining ability, which in turn allow the estimation of
and
(Griffing 1956; Hallauer and Miranda 1988):
![]() | (1a) |
![]() | (1b) |
![]() | (2a) |
![]() | (2b) |
![]() | (3a) |
![]() | (3b) |
), genotype by environment (
) and environmental (
) components, as illustrated in Table 1.
|
|
The within-family analysis allows estimation of the relative importance of epistatic effects. In the absence of epistasis the following equation holds true (Hallauer and Miranda 1988):
![]() | (4) |
![]() | (5) |
) = 0 and 4 Cov(
) = 0, the formula can be simplified:
![]() | (6) |
![]() |
![]() |
![]() |
![]() |
![]() |
Equation 6 can now be written as follows:
![]() |
The analysis of between-family variation was published elsewhere (Cach et al. forthcoming). In that article, genetic effects, rather than genetic variances, were of interest, and they were considered fixed effects. In the present study, however, the analysis of within-family variance and the relative importance of epistatic effects are of prime interest. All effects, therefore, were considered random and normally distributed. The 30 genotypes representing each F1 cross are clearly a random sample of all possible genotypes that could possibly be derived from the respective parents. The only criterion defining which genotype would be used was the capacity to produce six stakes in an environment different from the target environment where the evaluation was conducted. The parents involved in this study were among a group of 2530 clones characterized by their adaptation to subhumid conditions: long periods without rain, tolerance or resistance to insect and arthropod pests typical for these environments (particularly thrips and different species of mites), and a frequent susceptibility to foliar diseases (because they are not common in this kind of environment). Eight of the parents evaluated come from CIAT's cassava-breeding project in Colombia, and the remaining clone was a cultivar released many years ago in Thailand. These parents are considered to be part of a reference population of clones adapted to the subhumid, lowland, tropical environment.
The actual nine parents eventually included were those that allowed for a balanced set of progenies for the study. Therefore, the main criterion for the selection of the parental lines was their capacity to flower and produce adequate samples of botanical seed from many different crosses. It is difficult to assess the impact (if any) of this selection because crossings are made in the midaltitude valleys environment where CIAT headquarters are located, but the evaluation was conducted in a completely different environment. This is important because the flowering habit, which profoundly affects plant architecture, varies drastically from one environment to the other. A nonbranching, erect type in the subhumid environment may be bushy and flower profusely at Palmira. Because of this situation, it can be assumed that the effect of selection of parents at Palmira had a neutral impact on the general performance of the progenies selected and evaluated for this study.
The analysis of variance for the between-family variation follows Method 4 proposed by Griffing (1956). The usual assumptions for Method 4 analysis are: regular diploid behavior during meiosis, absence of cytoplasmic effects, linkage equilibrium, relatives are random members of a specified population, and, because of the vegetative propagation of cassava, negligible C-effects (Libby and Jund 1962). In the case of cassava, C-effects would result from differences in the physiological/sanitary status between F1 mother plants and/or among the six stakes used to clone each genotype and these differences would be confounded with the environmental and/or genotype x environment interactions components of variation. Because the F1 plants from which the six stakes were taken had been grown in Palmira under excellent management practices, differences (if any) in the physiological/sanitary status of these vegetative cuttings are reasonably expected to be small and negligible.
| Results |
|---|
|
|
|---|
The coefficients of variability indicated that the experimental error involved in this study was relatively low. Results, therefore are reliable and the precision of the analysis, acceptable (Cach et al. forthcoming). The two locations used in the evaluation showed statistical differences only for foliage yield and harvest index (Table 2). On the other hand statistical differences among crosses were found for all the variables analyzed. With the exception of the thrips score, the crosses by environment interactions were also significant. GCA mean squares were significant for thrips score and harvest index (Table 2). SCA mean squares were also significant for all variables except dry matter content.
Because individual clone data have been included, the degrees of freedom involved are considerably larger (Table 2) than those reported in the between family analysis (Cach et al. forthcoming). In every case, within-family genetic variation (
) was statistically significant. The interaction between environment and the within-family genetic variation also proved to be statistically significant. From the mean squares presented in Table 2 the estimates for
and the test for epistasis were obtained as already described.
Variance components were considered important if the standard errors were less than half of the component estimates (Isik et al. 2003). The estimate for
was larger than that for
for fresh root and foliage yields, harvest index, and dry matter yield and smaller for reaction to thrips and dry matter content (Table 3). Epistasis was highly significant for all variables (test values > two times the value of their respective standard errors) except harvest index (Table 3).
|
| Discussion |
|---|
|
|
|---|
Based on the magnitude of the estimates for between- and within-family genetic variances, a large proportion of the genetic variability (7993%) remained as within family variation (Table 3). These results agree with observations during the selection in evaluation trials where large numbers of crosses among elite parental lines are represented by several clones. As expected, the lowest within-family variation (79% of total genetic variance) was measured for a relatively simply inherited trait, such as the reaction to thrips (Bellotti 2002), which showed the only statistically significant additive variance. The tolerance/resistance in outstanding parents transmitted to the progeny tended to accentuate differences among families and reduce the variability among sister clones. A similar situation was observed in a similar study for the midaltitude valleys environment (Pérez et al. forthcoming). However, it is clear that a considerable within-family variation still remained even for the reaction to thrips. On the other hand, complex traits such as root and foliage yields showed a larger partitioning of the total genetic variance (>90%) into the within-family variation, suggesting that there were, comparatively, smaller differences in the breeding values of the progenitors.
The within-family variation suggested not only important genetic effects but also significant genotype-by-environment variation for all variables analyzed. This interaction implies that reliable selection can only be made when enough planting material for replicated trials at more than one location has been produced. In practice, this means the third or fourth stage in the selection process (Ceballos et al. 2004). One alternative for overcoming this problem would be to modify the clonal evaluation trials (first stage in the selection process), which currently is conducted as an unreplicated trial at a single location, with seven plants per genotype (Ceballos et al. 2004). The total number of plants per genotype can be raised to eight, so that two trials at two different locations and with four plants per genotype at each location can be planted. Although the costs related to this change are large, and the logistic complications considerable, the data provided by this experiment (and other similar studies) suggest that they may be justifiable.
Dominance effects were very important for thrips, harvest index, and root and foliage yields, with variance estimates significantly different from zero (estimates two times or more the size of the respective standard error). Only the score for thrips and dry matter content showed larger estimates for the additive compared with the dominance variance (Table 3). This highlights the importance of heterosis in cassava breeding for many relevant traits, which in turn justifies the implementation of a reciprocal recurrent selection scheme for cassava genetic improvement.
Epistatic effects were significant for all variables, except harvest index, based on the test for epistasis (Table 3). It was surprising to see the size and generalized significance of epistatic effects. In many cases reported in the literature, epistatic effects may have been large but failed to reach statistical significance, in part because of the size of the standard errors typical for complex linear functions (Hallauer and Miranda 1988; Hinze and Lamkey 2003; Holland 2001). In this study, however, this was not the case. To a large extent this may be the result of the large size of this experiment, which resulted in large degrees of freedom for the overall analysis, including the number of clones within family and the number of replications and environments employed. However, the large and frequent epistasis found in this study may also be the result of the evolutionary history of this species that can multiply both sexually or clonally. It is feasible that cassava has evolved to take advantage of favorable gene combinations resulting from dominance and epistatic relationships by fixing them through the vegetative mode of reproduction. The results of this study reveal the limitation of most quantitative genetic studies based on the assumption of negligible epistasis. These results would also help explain the difficulties in finding QTLs that satisfactorily explain the phenotypic variation observed in complex traits, such as yield (Kao and Zeng 2002).
The phenotypic clonal selection used for cassava breeding takes advantage of the vegetative reproduction of the crop. In selecting outstanding clones all genetic effects (additive, dominance, and epistatic) are exploited (Ceballos et al. 2004; Mullin and Park 1992). However, the current recurrent selection system lacks the capacity to direct genetic improvement in such a way that the frequency of favorable (within or between loci) genetic combinations is maximized. To achieve this, special efforts to design parental clones that produce better crosses are required.
CIAT has recently introduced modifications that allow for the estimation of GCA effects in early stages of the selection process (Ceballos et al. 2004). This, in turn, allows the implementation of the Backward GCA Selection described by Mullin and Park (1992). Results from this study suggest that this approach would be ideal for traits such as the reaction to thrips, given the importance of GCA effects and the comparatively low relevance of dominance and epistatic effects. For complex traits such as fresh-root yield, however, the prevalence of nonadditive effects suggested by this study would require a different approach. The development of clones specifically designed for their use as parents in breeding nurseries would be one alternative that offers interesting advantages. Introduction of inbreeding in these parental clones would facilitate the gradual and consistent assembly of favorable gene combinations, which in the current system occur just by chance. Inbreeding would also facilitate the reduction of the genetic load of this crop, which is expected to be relatively large at this point in time.
One major constraint for the introduction of inbreeding in cassava is the time required for it. The production of doubled haploids through anther or microspore culture is an interesting approach that would reduce the time required to obtain homozygous genotypes. This, in turn, will maximize the exploitation of dominance and epistatic genetic variation (in crosses), which have been found to be significant in this study. CIAT is currently executing a project financed by the Rockefeller Foundation to develop the protocol for the production of doubled haploids in cassava.
| Acknowledgments |
|---|
The technical assistance of Drs. B. Li and R. Macchiavelli for the statistical analysis of the within family variation is greatly appreciated. Any mistake in the analysis, however, is solely responsibility of the authors. This research was conducted with resources provided by the Ministry of Agriculture and Rural Development of Colombia.
| Footnotes |
|---|
Corresponding Editor: William Tracy
Received November 10, 2004
Accepted April 11, 2005
| References |
|---|
|
|
|---|
-
Becker WA, 1985. Manual of quantitative genetics, 4th ed Pullman, WA: Academic Enterprises; 9596.
Bellotti AC, 2002. Arthropod pests In: Cassava: biology, production and utilization (Hillocks RJ, Thresh JM, and Bellotti AC, eds). Wallingford, NY: CABI Publishing; 209235.
Brim CA and Cockerham CC, 1961. Inheritance of quantitative characters in soybeans. Crop Sci 1:187190.
Busch RH, Janke JC, and Frohberg RC, 1974. Evaluation of crosses among high and low yielding parents of spring wheat (Triticum aestivum L.) and bulk prediction of line performance. Crop Sci 14:4750.
Cach NT, Lenis JI, Perez JC, Morante N, Calle F, and Ceballos H, forthcoming. Inheritance of agronomically relevant traits in cassava (Manihot esculenta Crantz) for sub-humid conditions Plant Breedi (in press).
Ceballos H, Iglesias CA, Pérez JC, and Dixon AGO, 2004. Cassava breeding: opportunities and challenges. Plant Mol Biol 56:503515.[CrossRef][ISI][Medline]
Ceballos H, Pandey S, Narro L, and Perez JC, 1998. Additive, dominance, and epistatic effects for maize grain yield in acid and non-acid soils. Theor Appl Genet 96:662668.[CrossRef]
CIAT, 2002. Annual report from IP3 project: improved cassava for the developing world Cali, Colombia: CIAT, Apdo Aéreo 6713.
Comstock RE, Kelleher T, and Morrow EB, 1958. Genetic variation in an asexual species, the garden strawberry. Genetics 43:634646.
Cortes DF, Reily K, Okogbenin E, Beeching JR, Iglesias C, and Tohme J, 2002. Mapping wound-response genes involved in post-harvest physiological deterioration (PPD) of cassava (Manihot esculenta Crantz). Euphytica 128:4753.
Easwari Amma CS and Sheela MN, 1998. Genetic analysis in a diallel cross of inbred lines of cassava. Madras Agric J 85(5,6):264268.
Easwari Amma CS, Sheela MN, and Thankamma Pillai PK, 1995. Combining ability analysis in cassava. J Root Crops 21(2):6571.
Eta-Ndu JT and Openshaw SJ, 1999. Epistasis for grain yield in two F2 populations of maize. Crop Sci 39:346352.
FAO/FIDA, 2000. La economía mundial de la yuca. Hechos, tendencias y perspectivas. Fondo Internacional de Desarrollo Agrícola Organización de las Naciones Unidas para la Agricultura y la Alimentación. Roma, Italy.
Foster GS and Shaw DV, 1988. Using clonal replicates to explore genetic variation in a perennial plant species. Theor Appl Genet 76:788794.[CrossRef]
Fregene M, Angel F, Gomez G, Rodriguez F, Chavarriaga P, Roca W, Tohme J, and Bonierbale M, 1997. A molecular genetic map of cassava. Theor Appl Genet 95:431441.[CrossRef][ISI]
Griffing B, 1956. Concept of general and specific combining ability in relation to diallel crossing systems. Aust J Biol Sci 9:46393.
Gravois KA, 1994. Diallel analysis of head rice percentage, total milled rice percentage, and rough rice yield. Crop Sci 34:4245.
Hallauer AR and Miranda Fo JB, 1988. Quantitative genetics in maize breeding, 2nd ed Ames: Iowa State University Press; 45114.
Hanson WD and Weber CR, 1961. Resolution of genetic variability in self-pollinated species with an application to the soybean. Genetics 46:14251434.
Hinze LL and Lamkey KR, 2003. Absence of epistasis for grain yield in elite maize hybrids. Crop Sci 43:4656.
Holland JB, 2001. Epistasis and plant breeding. Plant Breed Rev 21:2792.
Isik F, Li B, and Frampton J, 2003. Estimates of additive, dominance and epistatic genetic variances from a clonally replicated test of loblolly pine. Forest Sci 49(1):7788.
Jorge V, Fregene M, Duque MC, Bonierbale MW, Tohme J, and Verdier V, 2000. Genetic mapping of resistance to bacterial blight disease in cassava (Manihot esculenta Crantz). Theor Appl Genet 101:865872.[CrossRef]
Jorge V, Fregene M, Velez CM, Duque MC, Tohme J, and Verdier V, 2001. QTL analysis of field resistance to Xanthomonas axonopodis pv manihotis in cassava. Theor Appl Genet 102:564571.[CrossRef]
Kao CH and Zeng ZB, 2002. Modeling epistasis of quantitative trait loci using Cockerham's model. Genetics 160(3):12431261.
Kawano K, 1980. Cassava In: Hybridization of crop plants (Fehr WR and Hadley HH, eds). Madison, WI: ASA, CSSA; 225233.
Kawano K, Gonçalvez Fukuda WM, and Cenpukdee U, 1987. Genetic and environmental effects on dry matter content of cassava root. Crop Sci 27:6974.
Lamkey KR, Schnicker BJ, and Melchinger AE, 1995. Epistasis in an elite maize hybrid and choice of generation for inbred line development. Crop Sci 35:12721281.
Libby WJ and Jund E, 1962. Variance associated with cloning. Heredity 17:533540.
Losada VT, 1990. Cruzamentos dialélicos em mandioca (Manihot esculenta Crantz) Ph.D. dissertation. Piracicaba, SP, Brazil. Escola Superior de Agricultura Luiz de Queiroz. Universidade de São Paulo.
Lynch M and Walsh B, 1998. Genetics and analysis of quantitative traits Sinauer Associates, esp. pp 558563 and 813816.
Mba REC, Stephenson P, Edwards K, Melzer S, Mkumbira J, Gullberg U, Apel K, Gale M, Tohme J, and Fregene M, 2001. Simple sequence repeat (SSR) markers survey of the cassava (Manihot esculenta Crantz) genome: towards an SSR-based molecular genetic map of cassava. Theor Appl Genet 102:2131.[CrossRef]
Melchinger AE, Geiger HH, and Schnell FW, 1986. Epistasis in maize (Zea mays L.). 2. Genetic effects in crosses among early flint and dent inbred lines determined by three methods. Theor Appl Genet 72:231239.[CrossRef][ISI]
Mullin TJ and Park YS, 1992. Estimating genetic gains from alternative breeding stratregies for clonal forestry. Can J For Res 22:1423.
Okogbenin E and Fregene M, 2003. Genetic mapping of QTLs affecting productivity and plant architecture in a full-sib cross from non-imbred parents in cassava (Manihot esculenta Crantz). Theor Appl Genet 107:14521462.[Medline]
Orf JH, Chase K, Adler FR, Mansur LM, and Lark KG, 1999. Genetics of soybean agronomic traits: II. Interactions between yield quantitative trait loci in soybean. Crop Sci 39:16521657.
Pérez JC, Ceballos H, Jaramillo G, Morante N, Calle F, Arias B, and Bellotti AC, forthcoming. Epistasis in cassava (Manihot esculenta Crantz) adapted to mid-altitude valley environment Crop Sci (in press).
Pixley KV and Frey KJ, 1991. Combining ability for test weight and agronomic traits of oat. Crop Sci 31:14481451.
Rönnberg-Wästljung AC and Gullberg U, 1999. Genetics of breeding characters with possible effects on biomass production in Salix viminalis (L.). Theor Appl Genet 98:531540.[CrossRef]
Rönnberg-Wästljung AC, Gullberg U, and Nilsson C, 1994. Genetic parameters of growth characters in Salix viminalis grown in Sweden. Can J For Res 24:19601969.
Scott GJ, Rosegrant MK, and Ringler C, 2000. Global projections for root and tuber crops to the year 2020. Food Policy 25:561597.[CrossRef]
Stonecypher RW and McCullough RB, 1986. Estimates of additive and non-additive genetic variance from a cloneal diallel of Douglas-fir Pseudotsuga menziesii (Mirb.) Franco Proc Int Union For Res Org, Join Mtg Working Parties Breed Theor, Prog Test, Seed Orch. Williamsburg, VA: NCSU-Industry Coop Tree Imp Program; 211227.
Vega OPC, 1987. Introducción a la teoría de genética cuantitativa Caracas: Universidad Central de Venezuela Press.
Vencovsky R and Barriga P, 1992. Genética Biométrica no Fitomelhoramento Ribeirão Preto, Brasil: Sociedade Brasileira de Genética.
Wolf DP and Hallauer AR, 1997. Triple test cross analysis to detect epistasis in maize. Crop Sci 37:763770.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||















