|
|
||||||||
a Dep. of Agron., Iowa State Univ., Ames, IA 50011-1010
b Dep. of Agron. and Environ. Manage., Louisiana State Univ. Agric. Cent., Baton Rouge, LA 70803-2110
* Corresponding author (mkang{at}agcenter.lsu.edu)
Received for publication October 8, 2004.
| ABSTRACT |
|---|
|
|
|---|
2g) and SCA (
2s) variances, which can be used to estimate additive (
2A) and dominance (
2D) variances, and ultimately narrow-sense heritability (h2). DIALLEL-SAS05 also provides information on GCA x ENV, SCA x ENV, reciprocal x ENV, maternal x ENV, and nonmaternal x ENV interactions, when applicable. DIALLEL-SAS05 should greatly improve researchers' efficiency in analyzing and interpreting diallel-cross data. The program code is available on a CD from the corresponding author.
Abbreviations: GCA, general combining ability GEAN, GardnerEberhart analysis SCA, specific combining ability
| INTRODUCTION |
|---|
|
|
|---|
Gardner and Eberhart (1966) proposed alternative analyses of data from diallel crosses produced from heterogeneous parents/populations ("varieties"). GardnerEberhart Analysis (GEAN) I is resource-intensive, requiring the evaluation of n parents, n(n 1)/2 F1 crosses, and inbred progeny of parents and crosses, but provides information on additive and dominance gene action, heterosis, and inbreeding depression (Gardner and Eberhart, 1966; Murray et al., 2003). The GEAN II is useful in evaluating n populations (varieties) and their n(n 1)/2 F1 crosses; variation among populations is partitioned into populations and midparent heterosis (Gardner and Eberhart, 1966; Hallauer and Miranda, 1988; Murray et al., 2003). Heterosis is further partitioned into average, variety, and specific heterosis, but additive and dominance parameters cannot be determined in Analysis II because they are confounded with the source of variation labeled "variety" (Murray et al., 2003). The GEAN III provides estimates of variety and GCA effects from an analysis that contains the following sources of variation: parents, parents vs. F1 crosses, and F1 crosses; GCA effects are estimated in a manner similar to Griffing's Method 4, Model 1 (Murray et al., 2003). Both GEAN II and GEAN III provide estimates of average heterosis and SCA (Gardner and Eberhart, 1966; Murray et al., 2003). Many researchers have fruitfully used GEAN for studying heterosis and estimating GCA and SCA in variety diallels (Crossa et al., 1987; Ali et al., 2001; Lee et al., 2003).
Murray et al. (2003) pointed out that formulas for the effects for GEAN II and GEAN III were nonintuitive both biologically and genetically and incorporated the number of parents as multipliers. They further opined that (i) the variety effects obtained in GEAN III were "unconstrained" estimates, whereas those from GEAN II were constrained estimates because of the assumption of no heterosis; and (ii) the results had implications for the use and interpretation of such effects. For further specific details about GEAN I, II, and III, the reader is referred to Murray et al. (2003).
Naturally because of cumbersome calculations, researchers have felt the need to develop computer programs for analyzing diallel data. Program codes have been developed for use in crop breeding (Burrow and Coors, 1994; Kang, 1994; Magari and Kang, 1994; Zhang and Kang, 1997) and tree breeding (Wu and Matheson, 2000, 2001; Xiang and Li, 2001). Murray et al. (2003) reported a SAS computer program for GEAN II and GEAN III; the program used the general linear model with cell means to estimate various genetic effects.
The DIALLEL-SAS program of Zhang and Kang (1997) has been extensively usedfor example, to analyze diallel-cross data in maize (Zea mays L.) (Kang et al., 1999, 2005; Naidoo et al., 2002; Lee et al., 2003), rapeseed (Brassica napus L.) (Goffman and Becker, 2001), wheat (Triticum aestivum L.) (Le Gouis et al., 2002; Hakizimana et al., 2004), durum wheat (Triticum turgidum L. var. durum) (Solomon and Labuschagne, 2004), alfalfa (Medicago sativa L.) (Guines et al., 2002; Riday and Brummer, 2002a, 2002b), and white lupin (Lupinus albus L.) (Lagunes-Espinoza et al., 2003). The usefulness of DIALLEL-SAS can be enhanced, however. The limitations of DIALLEL-SAS (Zhang and Kang, 1997) are (i) while it can be easily adapted for a particular diallel set (Griffing's Methods 1 to 4) for two environments, the user must modify the code to analyze data from more than two environments, which can be cumbersome; (ii) it provided results only for a fixed-effects model; (iii) it required extensive and laborious program modifications for parent number other than n = 5 used in the example programs in Zhang and Kang (1997); and (iv) it did not include Gardner and Eberhart (1966) Analyses II and III.
Currently, there is no known computer program that handles analyses of both Griffing's four diallel methods and GardnerEberhart Analyses II and III. Thus, a need exists for a comprehensive program that computes not only Griffing's diallel analysesboth fixed- and random-effects modelsbut also GardnerEberhart Analyses II and III. We report here the development of such a program called DIALLEL-SAS05.
| PROGRAM DESCRIPTION |
|---|
|
|
|---|
In the current program, DIALLEL-SAS05, statistical models used for GEAN2 and GEAN3 are as follows (Gardner and Eberhart, 1966; Murray et al., 2003):
![]() |
![]() |
j'; µv = mean of all parents; vj and vj' are effects for variety j and j', respectively, when parents are included in the analysis (restriction:
vj = 0);
= average heterosis;
= 0 when j = j' and
= 1 when j
j'; hj and hj' = variety heterosis for variety j and j', respectively; sjj' = specific heterosis in GEAN2 and SCA in GEAN3 (restriction
sjj' =
sjj = 0); xjj' = the cross effect for the mating between varieties j and j'; and gj and gj' are GCA effects for variety j and j', respectively.
Program Components
The DIALLEL-SAS05 program consists of two major parts. The first part is a user-interface program called "Diallel.sas." This is where the user inputs data and invokes necessary macro functions. The second part is called "DialAll.sas," which includes seven major macro functions and other utility macros. A flow chart for functional relationship among different macros is given in Fig. 1
.
|
| AN EXAMPLE |
|---|
|
|
|---|
OPTIONS PS=56 LS=100 mLOGIC PAGENO=1;
%INCLUDE 'C:\DIALLEL\DIALALL.SAS';
DATA DT1;
INPUT I J REP ENTRY YIELD ENV;
IF I=J THEN AVGH=0;ELSE AVGH=1;
CARDS;
1 1 1 1 10.5 1
1 1 2 1 10.7 1
1 2 1 2 11.9 1
1 2 2 2 12.0 1
1 3 1 3 14.5 1
1 3 2 3 14.2 1
1 4 1 4 9.0 1
1 4 2 4 8.5 1
1 5 1 5 13.5 1
1 5 2 5 14.2 1
.
more data records
.
5 1 1 21 22.1 2
5 1 2 21 21.0 2
5 2 1 22 21.9 2
5 2 2 22 20.2 2
5 3 1 23 20.2 2
5 3 2 23 20.8 2
5 4 1 24 21.2 2
5 4 2 24 20.6 2
5 5 1 25 20.2 2
5 5 2 25 21.0 2
;
RUN;
/* FOR GRIFFING'S METHOD 2, DO THE FOLLOWING TO DELETE RECIPROCAL CROSSES */
DATA DT2;SET DT1;
IF I>J THEN DELETE;
RUN;
/* FOR GRIFFING'S METHOD 3, DELETE INBRED LINES AS FOLLOWS */
DATA DT3;SET DT1;
IF I=J THEN DELETE;
RUN;
/* FOR GRIFFING'S METHOD 4, DELETE INBRED LINES AND RECIPROCAL CROSSES AS FOLLOWS */
DATA DT4;SET DT1;
IF I>=J THEN DELETE;
RUN;
/* TO MAKE DATASET FOR GARDNER-EBERHART ANALYSIS II AND III WHEN PLOT MEAN DATA ARE USED, DO THE FOLLOWING (NOTE: YOU WILL NEED TO REPLACE YIELD IN AVG(YIELD) WITH APPROPRIATE VARIABLE BEING ANALYZED): */
PROC SQL;
CREATE TABLE DT5 AS
SELECT I,J,ENTRY,AVG(YIELD) AS AVG_Y,AVGH
FROM DT1
GROUP BY I,J,ENTRY;
/* GENERATING DESIGN MATRICES ONLY */
%MGENERATOR(NUM_P=5,DSN=DT1, METHOD=1);
%GARDNERANAL2(NUM_P=5, DSN=DT2);
%GARDNERANAL3(NUM_P=5, DSN=DT2);
/* FOR GRIFFINGS FIXED-EFFECTS MODEL */
%DIALANALFIXMODEL(NUM_P=5,METHOD=
1,YVAR=YIELD,ENV=2,REP=2,DSN=DT1);
%DIALANALFIXMODEL(NUM_P=5,METHOD=
2,YVAR=YIELD,ENV=2,REP=2,DSN=DT2);
%DIALANALFIXMODEL(NUM_P= 5,METHOD=
3,YVAR=YIELD,ENV=2,REP=2,DSN=DT3);
%DIALANALFIXMODEL(NUM_P= 5,METHOD=
4,YVAR=YIELD,ENV=2, REP=2, DSN=DT4);
/* FOR GRIFFINGS RANDOM-EFFECTS MODEL */
%DIALANALRANMODEL(NUM_P=5, METHOD=
1,YVAR=YIELD, ENV=2, REP=2, DSN=DT1);
%DIALANALRANMODEL(NUM_P=5, METHOD=
2,YVAR=YIELD, ENV=2, REP=2, DSN=DT2);
%DIALANALRANMODEL(NUM_P= 5, METHOD=
3,YVAR=YIELD, ENV=2, REP=2, DSN=DT3);
%DIALANALRANMODEL(NUM_P=5, METHOD=
4,YVAR=YIELD, ENV=2, REP=2, DSN=DT4);
/* GARDNER AND EBERHART ANALYSES */
/* FOR GARDNER-EBERHART ANALYSIS II AND ANALYSIS III, RECIPROCAL CROSSES ARE DELETED(SEE DATA DT2 ABOVE) */
TITLE "GARDNER AND EBERHART ANALYSIS II";
%DIALANALYSISII(NUM_P=5, YVAR=YIELD,
ENV=2,REP=2, DSN=DT2, MSE=1,
DFE=0);
TITLE "ANALYSIS II, MSE=1.2";
%DIALANALYSISII(NUM_P=5, YVAR=AVG_Y,
ENV=1,REP=1, DSN=DT5, MSE=1.2,
DFE=28);
TITLE "GARDNER AND EBERHART ANALYSIS III";
%DIALANALYSISIII(NUM_P=5, YVAR=YIELD,
ENV=2,REP=2,DSN=DT2, MSE=1,
DFE=0);
TITLE "ANALYSIS III, MSE=1.2";
%DIALANALYSISIII(NUM_P=5, YVAR=AVG_Y,
ENV=1,REP=1,DSN=DT5, MSE=1.2,
DFE=28);
RUN;
QUIT;
In the above SAS statements relative to the GardnerEberhart analyses, if MSE
0, the program recognizes that the data have multiple replications. If MSE is >0 (actual value to be provided by user), the program recognizes that data represent means across replications.
For brevity, only sample outputs for Griffing's Method 1 for a fixed-effects model or Model 1 (Fig. 2) and a random-effects model or Model 2 (Fig. 3) are shown. The results for GEAN II and III are given in Fig. 4 and 5 , respectively. That part of the output that is identical for GEAN II and GEAN III has been omitted from Fig. 5. Only essential information of various outputs has been shown in Fig. 2 through 5.
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
2g,
2s, and reciprocals (
2r). Kuehl et al. (1968) used a diallel design to estimate
2A,
2D, and
2AA. Kang (1994) provided computations for estimating
2A,
2D, and narrow-sense heritability (h2) for all four Griffing's methods. DIALLEL-SAS05 is more user-friendly, comprehensive, and less cumbersome to modify than the DIALLEL-SAS program of Zhang and Kang (1997). The major features of the new program are
2g,
2s,
2r, and their variances are computed for the random-effects model to estimate
2A,
2D, and h2. To compute
2A and
2D (e.g., from the information given in Fig. 3), one would do the following:
When inbreeding coefficient (F) of parents = 0 (no inbreeding)
![]() |
When inbreeding coefficient (F) of parents = 1 (100% inbreeding)
![]() |
Software Installation and Support
The DIALLEL-SAS05 program was developed and tested on SAS 9.0 on workstation with Windows XP (SAS Inst., 2002). It can run on any SAS version 8.0 or above, however. The program is expected to run on other SAS versions without any significant problems/modifications, but we did not test it on older SAS versions.
The software is contained in two files: Diallel.sas and DialAll.sas. The DialAll.sas should be saved in a directory in the user's local PC or a server that the user can access. For example, the user may save it in C:\Documents\DialAll.sas. The other file, Diallel.sas, may be saved in the same folder or in a different folder in user's personal computer. Once the two files have been saved, the user would open Diallel.sas with the SAS software. Upon opening Diallel.sas, the user will need to change the %include file path to the user's own specific path (%include C:\Documents\DialAll.sas;). Now, the user is ready to input data and run any of the seven macro functions. To analyze data for more than one trait, the user would need to have as many statements as the number of traits. For example, to analyze five-parent diallel data on yield and testwt for Griffing's Method 1, Model 1, user must include the following statements:
%DIALANALFIXMODEL(NUM_P=5,METHOD=
1,YVAR=YIELD,ENV=2,REP=2,DSN=DT1);
%DIALANALFIXMODEL(NUM_P=5,METHOD=
1,YVAR=TESTWT,ENV=2,REP=2,DSN=DT1);
In summary, DIALLEL-SAS05 has the following major functions: (i) it deals with both fixed- and random-effects models relative to all four Griffing's diallel methods, (ii) it computes GardnerEberhart's Analyses II and III, (iii) it is effective for analyzing diallel experiments with number of parents ranging from 4 to 12, (iv) it can analyze diallel data for any number of environments, and (v) it provides information on GCA x ENV, SCA x ENV, reciprocal x ENV, maternal x ENV, and nonmaternal x ENV interactions, when applicable. The user can generate a linear matrix for any number of parents in a diallel cross by easily extending the program code to diallel designs with more than 12 parents. DIALLEL-SAS05 is expected to greatly enhance researchers' efficiency in analyzing and interpreting diallel-cross data.
The DIALLEL-SAS05 software will be provided free of charge to those affiliated with academic institutes. The user may request a copy of the software on a CD, along with documentation/instructions regarding DIALLEL-SAS05, by sending US$10 to the corresponding author to defray the cost of a CD and shipping and handling. Outputs for the five-parent example for all four methods of Griffing (both fixed and random effects) and for GEAN II and GEAN III will be provided along with other documentation for the software.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S.M. Smith, J.W. Scott, J.A. Bartz, and S.A. Sargent Diallel Analysis of Fruit Water Absorption in Tomato, a Contributing Factor in Postharvest Decays J. Amer. Soc. Hort. Sci., January 1, 2008; 133(1): 55 - 60. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. I. Ragsdale and C. W. Smith Germplasm Potential for Trait Improvement in Upland Cotton: Diallel Analysis of Within-Boll Seed Yield Components Crop Sci., May 31, 2007; 47(3): 1013 - 1017. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Sambandan, A. Yamamoto, J.-J. Fanara, T. F. C. Mackay, and R. R. H. Anholt Dynamic Genetic Interactions Determine Odor-Guided Behavior in Drosophila melanogaster Genetics, November 1, 2006; 174(3): 1349 - 1363. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. de la Vega and S. C. Chapman Multivariate Analyses to Display Interactions between Environment and General or Specific Combining Ability in Hybrid Crops Crop Sci., February 24, 2006; 46(2): 957 - 967. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Crop Science | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||