In this article we are presenting method of construction of three series of optimal partial triallel cross designs through complete diallel cross designs which have been obtained by nested balanced incomplete block designs.
Mating design play very important role to study the genetic properties of a set of inbred lines in plant breeding experiments. Most commonly used mating designs are diallel or a two-way cross . Suppose there are inbred lines and it is desired to perform a diallel cross experiment involving crosses of the type for , this type of mating design is called complete diallel cross (CDC) method (4) of Griffing (1956).
Triallel crosses form an important class of mating designs, which are used for studying the genetic properties of a set of inbred lines in plant breeding experiments. For inbred lines, the number of different crosses for a complete triallel experiment is of the type . Rawlings and Cockerham (1962) were the first to introduce mating designs for triallel crosses. Let denote the total number of crosses (experimental units) involved in a triallel experiment. It is desired to compare the lines with respect to their general combining abilities, the specific combining abilities being not included in the model.
Triallel cross (TC) experiments are generally conducted using a completely randomized design (CRD) or a randomized complete block (RCB) design as environmental design involving crosses. Even with a moderate number of parents, say , in a TC experiment; the number of crosses becomes unmanageable to be accommodated in homogeneous blocks. For such situations, Hinkelmann (1965) developed partial triallel crosses (PTC) involving only a sample of all possible crosses by establishing a correspondence between PTC and generalized partially balanced incomplete block designs (GPBIBD). Ponnuswamy and Srinivasan (1991) and Subbarayan (1992) obtained PTC using a class of balanced incomplete block (BIB) designs.
Other research workers who contributed in this area are Arora and Aggarwal (1984, p. 89), Ceranka et al. (1990). Recently Harun et al. (2016a, 2016b) constructed efficient PTC designs using mutually orthogonal latin squares and two associate partially balanced incomplete block designs and also gave a simple method of construction of variance – balanced PTC designs for comparing set of test treatments with a control line. More details on TC experiments can be found in Hinkelmann (1975) and Narain (1990).
Following Gupta and Kageyama (1994) and Dey and Midha (1996), Das and Gupta (1997) constructed block designs for TC by using the nested balanced block design (NBIBD) with parameters . Their method yields designs which are universally optimal in , the class of connected block designs for TC in lines with blocks each of size such that the total number of experimental units are . Sharma et al. (2011) and Sharma and Fanta (2012) obtained optimal designs for PTC by using mutually orthogonal Latin squares and gave their analysis along with first line effects and second line effects.
In this paper we are presenting methods of construction of block designs for TC experiments through block design of CDC experiment method (4). These designs are found to be optimal in the sense of Das and Gupta (1997). The paper is structured as: in Section 2 we gave some definitions and in Section 3 we discussed optimality tools. In Section 4 we discussed the construction and optimality of these designs. In Section 5 we have discussed the utility of these designs.
Some definitions
Definition: The triallel cross (T.C.) has been defined by Rawlings and Cockerham (1962) as a set of all possible three-way hybrids among a group of (inbred) lines. Given three lines and , there are distinct triallel crosses, namely and involving these three lines.
Thus given a set of lines, the TC will consist of a set of three way crosses.
Hinkelmann (1965) proposed the definition of PTC as given below:
Suppose we have lines which are denoted by . A three way cross is then represented by a triplet , where stands for an offspring of the single cross . We shall call and half-parents and full-parent. The crosses , and are considered to be identical in three way crosses. Then PTC can be defined as follows:
A set of matings is said to be a PTC if it satisfies the following conditions:
Each line occurs exactly times as half-parent and times as full parent.
Each cross occurs either once or not at all.
The total number of crosses is and . Let , whence .
Optimality tool
Following Parsad et al. (2005), let be a block design for an -allel cross experiment involving inbred lines, blocks each of size , where or 3 in our context. This means that there are crosses in each of the blocks of . Further, let and denote the number of replication of the cross and the number of replications of the line in different crosses, respectively, in , where means number of crosses. Evidently, and , the total number of observations. In a TC experiment, the genotypic effect of the hybrid consists of single line effects, two line effects and three line specific effects. However, if we assume that for a PTC experiment (in which every line appears as half parent an equal number of times, say , and every line appears as full parent an equal number of times, say , and each of the crosses appears at most once) the two line specific effects and three line specific effects are not of importance, still the line effects are of two types viz effects as half parent and effect as full parent i.e. the ordering of lines in a triallel cross is important. Some plant breeders argue that these ordering effects can also be averaged over line effects. Das and Gupta (1997) considered the situations where ordering of lines in a triallel cross is not of importance. We also considered this and we took the following additive model for the observations obtained from design .
where is vector of observations, is the vector of ones, is the design matrix for lines and is an design matrix for blocks, that is, the element of (respectively, of ) is 1 if the observation pertains to the line (respectively, of block) and is zero otherwise. is a general mean, is a vector of line parameters, is a vector of block parameters and is an vector of residuals. It is assumed that vector is fixed and is normally distributed with and , where is the identity matrix of conformable order.
The method of least squares for the analysis of proposed design leads to the following reduced normal equations for estimating the linear functions of the general combining effects of lines under model Eq. (1).
where , and for , is the number of crosses in in which the lines and appear together. is the number of times the line occurs in block of and is the diagonal matrix of block sizes.
A design will be called connected if and only if rank, or equivalently, if and only if all elementary comparisons among general combining ability (gca) effects are estimable using . We denote by , the class of all such connected block design with lines, blocks each of size . To prove optimality of design , we need the following well known lemma [(see, e.g., Cheng, 1978, p1246).]
Lemma 1. For given positive integers and , the minimum of subject to , where ’s are non-negative integers, is obtained when of the ’s are equal to and , where denotes the largest integer not exceeding . The corresponding minimum of is .
We then have:
Theorem 1. For any design
where and for a square matrix , stands for the sum of the diagonal elements of .
Proof. Follows on the lines similar to that of Das and Gupta (1997), Das, Dey and Dean (1998) and Parsad, Gupta and Srivastava (2005).
For any design , we have
Now, since , using Lemma 1
Hence
By Lemma 1, the above equality is attained if and only if or , for ; . Further, using proposition 1 of Kiefer (1975), we have following theorem:
Theorem 2. Let be a block design for -allel crosses, satisfying
and
is completely symmetric.
Then is universally optimal in
where and for a square matrix , stands for the sum of the diagonal elements of .
Kiefer (1975) showed that a design is universally optimal in a relevant class of competing design if (i) the information matrix (the matrix) of the design is completely symmetric in the sense that has all its diagonal elements equal and all of its off-diagonal elements equal, and (ii) the matrix has maximum trace in the class of competing designs, that is, such a design minimize the average variance of the best linear unbiased estimators of all elementary contrasts among the parameters of interest i.e. the general combining ability in our context. Such a design has
where is an identity matrix of order and is a matrix of all ones. Furthermore, using all elementary contrasts among gca effects are estimated with variance
Now we will show a connection between CDC designs method (4) and PTC designs through NBIBDs of Preece (1967).
Definition 1. A NBIBD with parameters is a design for treatments, each replicated times with two systems of blocks such that:
The second system is nested within the first, with each block from the first system, called henceforth as ‘block’ containing exactly blocks from the second system, called hereafter as ‘sub-blocks’;
Ignoring the first system leaves a balanced incomplete block design with parameters .
The following parametric relations hold for a nested balanced incomplete block design:
Since we are using diallel cross designs for PTC experiments using NBIBD, so first we will establish relation of NBIBD with diallel cross designs. Consider now a NBIBD with parameters . If we identify the treatments of as lines of a CDC method (4) experiment and perform crosses among the lines appearing in the same sub block of , we get a block design for a CDC method (4) experiment involving lines with crosses, each replicated times, and blocks, each of size . Such a design ; also, for such a design 0 or 1 for . and
where is an identity matrix of order and is a unit column vector of order . Clearly given by Eq. (3) is completely symmetric and which equals the equality given in Theorem 2. Thus the design is universally optimal in and using each elementary contrast among general combining ability effects is estimated with a variance
Since 0 or 1 for , a PTC triallel cross designs can be derived from by attaching line with each cross in block provided line does not appear in block, where . Hence we get a block design for a PTC experiment involving lines with triallel crosses, each replicated times, and blocks, each of size . Such a design ; and, for a design for . and
The given by Eq. (5) is completely symmetric and which equals the equality given in Theorem 2. Thus the design is universally optimal in and using each elementary contrast among gca effects is estimated with a variance
Hence we can state the following theorem:
Theorem 3. The existence of a NBIBD with parameters implies the existence of two types of designs
A universally optimal incomplete block design for CDC experiment method (4);
A universally optimal PTC design .
Method of construction of PTC through diallel cross designs
Series 1: Let be a prime or a prime power and be a primitive element of the GF (). Consider the following initial blocks.
As shown by Dey et al. (1986), these initial blocks, when developed in the sense of Bose (1939), give rise to a nested balanced incomplete block design with parameters .
Following Sharma and Tadesse (2017) we can put the above blocks into single block and developing the single block over mod (), we obtained an optimal block design for diallel crosses with minimal number of experimental units with parameters , and 1. This diallel cross design can be converted into optimal partial triallel cross design with parameters , by attaching line with the crosses in block in which th line does not appear at all, where and considering rows as blocks,
Example 1. Let 2. We get the following two columns.
where is a primitive element of GF and the elements of GF are .
We can convert both columns in single column and adding successively the non-zero elements of GF to the contents of the single column, the CDC design method (4) is obtained as given below, where the lines have been relabeled 1–9, using the correspondence :
CDC design method (4)
The above design is an optimal diallel cross design for method (4) with parameters and 1.
Example 2. Now attaching the 1, 2,…, 9 elements, respectively, with the crosses of the 9 blocks because these elements do not appear in the respective blocks of the above design. Considering rows as blocks, we obtain optimal block design for triallel cross with parameters and 1, which fulfill the all conditions for PTC design. The design is given below.
Partial Triallel Cross Design Series 1
Series 2: Let be a prime or a prime power and be a primitive element of the Galois field of order , GF (p). Consider the initial blocks
Dey et al. (1986) showed that these initial blocks, when developed give a solution of a nested incomplete block design with parameters .
Similarly arranging the above initial blocks into single block as in series 1 and developing , will yield an optimal CDC design method (4) with parameters , and 1.
This design can be converted into partial triallel cross design with parameters , , and 1 by the procedure described above in Series 1.
Example 3. Let 2. Then we get the following two initial blocks.
Now we arrange these two blocks in a single block as given below.
Now developing the above block mod (13), we obtain optimal CDC design with parameters 13, 6, 13 and 1. Following the procedure of example 2, we can obtain block design for PTC with parameters 13, 13, 7 and 1.
Series 3: Let be a prime or a prime power, then cyclically developing the following columns
yields an optimal CDC design method (4) with parameters . A PTC design with parameters , and can be obtained by the procedure described in Series1.
Example 4. Let 3. Then 7 and developing the following columns mod (7)
yields optimal CDC design with parameters 7, 3, and 7 and 1. A PTC design with parameters 7, 3, 7 and 1 can be obtained by the procedure given in Example 2.
Partial Triallel Cross Design Series 2
Note: The columns form a nested balanced incomplete block design with parameters
Universally optimal block designs, with 13, obtained by the above method from NBIB designs of Morgan et al. (2001), are listed below in Table. Using these designs optimal partial triallel cross designs can be constructed.
Universally optimal block design for complete diallel crosses with 13 generated by using NBIB designs of Morgan, Preece and Rees (2001). The following PTC designs can be obtained using the procedures given in Section 4
S. No.
Source
1.
7
7
6
MPR2
2.
9
18
4
MPR5w
3.
9
9
8
MPR8
4.
11
11
5
MPR 14
5.
13
39
4
MPR 20w
6.
13
26
6
MPR21
7.
13
13
12
MP23
Remark: Sources can be found in Morgan, Preece and Rees (2001).
Conclusion
We have obtained three series of optimal PTC design using optimal CDC designs method (4). These designs require only experimental units i.e. th crosses in comparison to complete TC designs. According to Shunmugathai and Srinivasan (2012), these designs can be used to improve the quantitative traits of economic and nutritional importance in crops and animals. It has been established that the three way hybrids are more stable than pure lines and single cross hybrids and exhibits individual as well as population buffering mechanisms because of the genetic base are replica of lines are available.
References
1.
AroraB. S., & AggarwalK. R. (1984). Confounded trialllel experiments and their applications. Sankhya, B46, 54-63.
2.
AroraB. S., & AggarwalK. R. (1989). Triallel experiments with reciprocal effects. Jour. Ind. Soc. Agril. Stat., 41, 91-103.
3.
BoseR. C. (1939). On the construction of balanced incomplete block designs. Ann. Eugen., 9, 353-399.
4.
CerankaB.ChudzikH.DobekA., & KielczewskaH. (1990). Estimation of parameters for triallel crosses compared in block designs. Statist. Appl., 2, 27-35.
5.
ChengC. S. (1978). Optimality of certain asymmetrical experimental designs. Ann. Statist., 1239-1261.
6.
DasA., & GuptaS. (1997). Optimal block designs for triallel cross experiments. Commu. Stat.-theory and Methods, 26, 1767-1777.
7.
DasA.DeyA., and DeanA. M. (1998). Optimal designs for diallel cross experiments. Statist. Prob. Letters., 36, 427-436.
8.
DeyA.DasU. S., & BanerjeeA. K. (1986). On nested balanced incomplete block designs. Calcutta Statist. Assoc. Bull., 35, 161-167.
9.
DeyA., & MidhaC. K. (1996). Optimal designs for diallel crosses. Biometrika, 83(2), 484-489.
10.
DasA.DeyA., & DeanA. M. (1998). Optimal designs for diallel cross experiments. Statistics and Probabilty Letters, 36, 427-436.
11.
GriffingB. (1956). Concepts of general and specific combining ability in relation to diallel crossing system. Aus. J. Biol. Sci., 9, 463-493.
HarunM.VargheseC.VargheseE., & JaggiS. (2016a). Three-way cross designs for animal breeding experiments. The Indian Journal of Animal Sciences, 86(6), 99-00.
14.
HarunM.VargheseC.VargheseE., & JaggiS. (2016b). Triallel cross designs for comparing a set of test lines with a control line. The Indian Journal of Genetics and plant Breeding, 76(1), 84-87.
HinkelmannK. (1975). A Survey of Statistical Design and Linear Models, 243-269, Amsterdam; North Hollond.
17.
KieferJ. (1975). Construction and optimality of generalized Youden designs. In a Survey of Statistical Design and Linear Models, ED. J. N. Srivastava, (1975), pp. 333-353. Amsterdam: North Hollond.
18.
NarainP. (1990). Statistical Genetics. Wiley Eastern Limited, New Delhi.
19.
ParsadR.GuptaV. K., & GuptaS. (2005). Optimal designs for diallel and double cross experiments. Utilitas Mathematica, 68, 11-32.
20.
PonnuswamyK. N., & SrinivasanM. R. (1991). Construction of partial triallel crosses (PTC) using a class of balanced incomplete block designs (BIBD). Commu. Stat-Theory and Methods, A20, 3315-3.
SharmaM. K.AsedewoinA., & Misganaw (2011). Optimal partial triallel cross designs. Model Assisted Statistics and Applications, 6, 1-8.
24.
SharmaM. K., & FantaS. (2012). Optimal partial triallel cross. Jour. of Indian Society of Ag. Stat., 66(2), 289-294.
25.
SharmaM. K., & TadesseM. (2017). Optimal block and row-column designs for CDC methods. Sankhya Ser. B, 1-14. DOI 10.1007/s13571-016-0126-6.
26.
ShunmugathaiR., & SrinivasanM. R. (2012). Robustness of triallel cross experiments using NBIB mating designs against interchange of a cross. Internaional Journal of Mathematics and Scientific Computing, 2(2), 86-92.
27.
SubbarayanA. (1992). On the applications of pure cyclic triple system for plant breeding experiments. Journal of Applied Statistics, 19(4), 489-500.