Epistasis interaction of QTL effects as a genetic parameter influencing estimation of the genetic additive effect

RESEARCH ARTICLE

Epistasis interaction of QTL effects as a genetic parameter influencing estimation of the genetic additive effect

Jan Bocianowski

Department of Mathematical and Statistical Methods, Poznan University of Life Sciences, Poznan, Poland.

^{Send correspondence to}

ABSTRACT

Epistasis, an additive-by-additive interaction between quantitative trait loci, has been defined as a deviation from the sum of independent effects of individual genes. Epistasis between QTLs assayed in populations segregating for an entire genome has been found at a frequency close to that expected by chance alone. Recently, epistatic effects have been considered by many researchers as important for complex traits. In order to understand the genetic control of complex traits, it is necessary to clarify additive-by-additive interactions among genes. Herein we compare estimates of a parameter connected with the additive gene action calculated on the basis of two models: a model excluding epistasis and a model with additive-by-additive interaction effects. In this paper two data sets were analysed: 1) 150 barley doubled haploid lines derived from the Steptoe x Morex cross, and 2) 145 DH lines of barley obtained from the Harrington x TR306 cross. The results showed that in cases when the effect of epistasis was different from zero, the coefficient of determination was larger for the model with epistasis than for the one excluding epistasis. These results indicate that epistatic interaction plays an important role in controlling the expression of complex traits.

Keywords: additive gene action effect, barley, doubled haploid lines, epistasis, QTL.

Introduction

Epistasis (an additive-by-additive interaction between quantitative trait loci or a nonallelic interaction of homozygous loci) was recognized as early as 100 years ago by Bateson (1909) to describe a situation where the action of one locus masks the allelic effects at another locus. However, in classical statistical genetics, epistasis has been used as a statistical abstraction, so that less attention has been paid to the molecular and physiological nature of the gene interaction involved (Tachida and Cockerham, 1989).

Epistatic interactions of major genes, denoted as aa (Kearsey and Pooni, 1998), have often been inferred from modified segregation ratios (Eshed and Zamir, 1996). For genes affecting quantitative traits epistasis has been defined as a deviation from the sum of independent effects of individual genes (Falconer, 1989). Epistasis between quantitative trait loci (QTLs) assayed in populations segregating for an entire genome has been found at a frequency close to that expected by chance alone (Edwards et al., 1987; Doebley and Stec, 1991; Paterson et al., 1991; Stuber et al., 1992; De-Vicente and Tanksley, 1993; Lin et al., 1995; Xiao et al., 1995). Recently epistatic effects have been considered by many researchers as important for complex traits (Lark et al., 1995; Eshed and Zamir, 1996; Cockerham and Zeng, 1996; Yu et al., 1997; Conti et al., 2011; Gowda et al., 2011; Jiang et al., 2011; Li et al., 2011; Mao et al., 2011; Upadhyaya et al., 2011; Bocianowski, 2012a,b,c). Hence, genetic models for QTL mapping assuming no epistasis could lead to a biased estimation of QTL parameters.

A common problem reported so far as associated with the analyses of data is that analyses of single-locus QTLs and epistatic interactions were conducted separately using different analytical tools (Xing et al., 2002; Bocianowski, 2008, 2012a,b,c; Krajewski et al., 2012). Although both of the analytical tools can provide statistical estimates for the amount of the effects and the proportions of variance explained, it is necessary for a joint estimation to evaluate the relative importance of individual QTLs and epistatic interactions in determining the performance of these traits.

The aim of the current study was to compare estimates of a parameter connected with the additive gene action calculated on the basis of two models: a model excluding epistasis and a model with additive-by-additive interaction effects. To this end, two data sets were analysed: 1) 150 barley doubled haploid lines derived from the Steptoe x Morex cross, and 2) 145 doubled haploid lines of barley obtained from the HarringtonTR x 306 cross.

Material and Methods

Estimation method

If in an experiment n biparental homozygous (doubled haploid, DH) plant lines were observed, an n-vector of phenotypic mean observations y =[y₁ y₂... y_n]' and qn-vectors of marker genotype observations m_l, l = 1, 2, ..., q were obtained. The i-th element (i = 1, 2, ..., n) of vector m_lis equal to -1 or 1, depending on the parent's genotype exhibited by the i-th line.

The estimation of an additive effect of a QTL (denoted by a) is based on the assumption that genes responsible for the trait are closely linked to the observed molecular marker. By choosing p from all the observed markers we can explain the variability of the trait, and model observations for the lines as:

or

where 1 denotes the n-dimensional vector of ones, µ the general mean, X a (n × p)-dimensional matrix of the form X = [

... ], l₁, l₂, ..., l_pE {1, 2, ..., q} and β the p-dimensional vector of unknown parameters of the form β^' = [

... ]. Finally, Z denotes a matrix whose columns are products of certain columns of matrix X, γ is the vector of unknown parameters of the form γ^' = [

... ], and e the n-dimensional vector of random variables, such that E(e_i)=0, Var(e_i) = σ², and Cov(e_i, e_j)=0for i ≠ j, i, j = 1, 2, ..., n. The parameters, , ..., are the additive effects of the genes controlling the trait and the parameters , , ..., are the additive × additive interaction effects. We assume that the epistatic interaction effects reveal only loci with significant additive gene action effects. This assumption significantly reduces the number of potential significant effects and improves the usefulness of the regression model.

The total additive effect of genes influencing the traits, defined as the sum of absolute values of individual additive effects, can be found in Bocianowski and Krajewski (2009) as:

The total epistasis effect of gene pairs influencing the trait, defined as the sum of values of individual pairs effects, is given in Bocianowski (2012b) as:

Markers chosen for models (1) and (2) may, for example, be selected, by a stepwise regression procedure (Charcosset et al., 2001). Here we used a three-stage algorithm, in which selection was first made by a backward stepwise search conducted independently within all linkage groups and then markers selected in this way were placed in one group and subjected to the second backward selection (see Jansen and Stam, 1994). Finally, at the third stage, we considered situations in which selected markers were located on the chromosome very close to each other (closer than 5 cM). Because these markers are probably linked to one QTL, only the marker with the largest value of the statistic test was retained in the set. At the first and second stages the Bonferroni correction was applied to control type I error for multiple tests (Province, 2001). For epistasis only markers showing significant association with traits were tested.

The coefficient of determination (R²) was used to measure how the model fits the data and, in this study, how the amount of the phenotypic variance is explained by the marker effects and marker interaction effects.

Examples

To compare the estimates of a obtained by the two methods (excluding and including epistasis), the following data sets were used.

Example 1

The data concern 150 doubled haploid (DH) lines of barley obtained from the Steptoe × Morex cross, used in the North American Barley Genome Mapping (NABGM) project and tested in sixteen environments (Kleinhofs et al., 1993; Romagosa et al., 1996; GrainGenes database, Steptoe x Morex cross). The linkage map used consisted of 223 molecular markers, mostly RFLP, with a mean distance between markers of 5.66 cM. The lines were analysed for eight phenotypic traits (alpha amylase, AA; diastatic power, DP; grain protein, GP; grain yield, GY; height, H; heading date, HD; lodging, L; malt extract, ME; Hayes et al., 1993). Grain protein, lodging and malt extract were transformed by arcsin . Missing marker data were estimated by the method of Martinez and Curnow (1994), using non-missing data of flanking markers.

Example 2

The data also come from the NABGM project (Tinker et al., 1996; GrainGenes database, Hordeum) and concern 145 doubled haploid (DH) lines of barley obtained from the Harrington × TR306 cross. The lines were analysed for seven phenotypic traits (weight of grain harvested per unit area, GY; number of days from planting until emergence of 50% of heads on main tillers, HD; number of days from planting until physiological maturity, NM; plant height, H; lodging transformed by arcsin , L; 1000 kernel weight, KW; and test weight, TW). We used a map composed of 127 molecular markers (mostly RFLP) with the mean distance between markers of 10.62 cM. Results shown below concern observations from five environments (in four environments observations were made over two years).

The total additive effects of QTLs were estimated for model (1) -without epistasis and for model (2) -with epistasis, for each environment independently for both data sets.

Results

Table 1 presents estimates of the total additive effects for the 150 doubled haploid lines of barley obtained from the Steptoe × Morex cross calculated by an assumption of a lack of epistasis effect (model 1), as well as by an assumption of the effect of epistasis interactions of genes (model 2). The obtained results showed that in 27 cases (30%) no statistically significant epistasis interaction effects were found. In 24 cases the values of additive effects were lower after epistasis had been incorporated in the model, compared to the case when this effect was excluded. The largest decrease of values of additive effects was observed for AA in MN92 -by 49.97% (Table 1). In 39 cases, allowance for the epistasis effect caused an increase in the value of additive effects (Table 1). The largest increase of a amounted to 142.86% (for GP in MTi92). The percentage phenotypic variance explained by QTL effects and their epistasis effects was larger than R²100 for model (1), except for one case, ME in MN92, where the decrease amounted to 0.3% (Table 2). The maximal increase of the observed phenotypic variation was 16.6% (from 25.2% to 41.8% for GY in WA92). In ten cases, an increase in R² was observed, amounting to at least 10% (Table 2).

For the second dataset (146 doubled haploid lines of barley obtained from the HarringtonTR x 306 cross), no epistasis effects were found in 38 cases (Table 3). With the incorporation of the effect of epistasis interaction in 14 cases, an increase of additive effect value and a reduction in the value of a was denoted in 11 cases (Table 3). The largest reduction in the value of parameter a was 46.10% (for H in MB93), whereas the largest increase of an additive effect was 28.88% (for WG in QC93). In all the cases, when the effect of epistasis was different from zero, the coefficient of determination was larger for model (2) than for model (1).

Discussion

The identification of QTLs and the elucidation of their genetic control (main effects and their epistatic effects) are essential for the development of efficient marker-assisted selection (MAS), aimed at improving breeding efficiency (Govindaraj et al., 2009). A direct implication of epistasis, especially the involvement of QTLs in epistatic interactions, is that the effects of single-locus QTLs are mostly dependent on the genotypes of other loci, and the effect of a QTL can sometimes be negated by the genotypes of a second locus. Thus any attempt to utilize QTLs in breeding programs has to take into account such epistatic effects. It is worthy of note that, although this study revealed a large number of epistases events through statistical genetic analysis, many further studies are needed before we can fully understand the biological meaning of these phenomena.

The most-important results of this study relate to the statistical characterization of the genetic components that control the expression of the traits, including additive-byadditive epistatic interactions. Ma et al. (2007) observed that 37% of the main-effect QTLs were involved in the epistatic interactions in maize grain yield and its components. This indicated that many loci in epistatic effects might not have significant effects for studied trait alone but might affect its expression by epistatic effects with the other loci. The results obtained herein also suggest that some of the additive QTLs might be detected with effects confounded by epistatic effects, if the epistatic effects were ignored in QTL mapping.

Incorporation of epistasis interaction of QTLs provided a more comprehensive characterization of the analyzed DH lines. This is evidenced by the higher R² values for model (2) than for model (1), i.e. the one excluding epistasis (Tables 2 and ⁴). Thus it may be concluded that QTL epistasis is a significant component for the understanding of the genetic control of determined phenotypic values, while failure to include epistasis may result in an incomplete or even erroneous characterization of the analyzed. In the presence of epistasis, however, the control of only main-effect markers is insufficient, because the epistatic effects of QTLs will also show influences, particularly in the case of complex phenotypes (Li et al., 1997a,b). Thus, inclusion of interaction markers closely linked to epistatic QTLs in the statistical models is expected to improve the power and accuracy of QTL mapping. The phenomenon of a biased estimation of additive effects in the absence of important interaction effects has already been addressed by Zeng et al. (2005) in the analysis of simulation data. A significant proportion of the identified additive effect QTLs were involved in digenic interactions with background loci. Thus, the usual estimates of additive effects of a QTL can be confounded by interactions, which may change according to genetic backgrounds, environments, and other factors. This means that QTLs and the epistatic loci are interchangeable, depending on the genetic backgrounds and probably environments where they are identified. This study showed that, besides the main (additive) effect QTLs, epistatic QTLs also play a crucial role in determining phenotypic values. Even if the epistatic interactions of main effect QTLs limit their usefulness in MAS programmes (Tan et al., 2001), the pronounced individual additive effects of these QTLs are sufficient enough to recruit them for MAS (Govindaraj et al., 2009). Because of the interaction between different loci, the offspring phenotype will be largely influenced by the genetic background of the receptor line when marker-directed selection is carried out (Tan et al., 2001).

The results obtained herein reinforce the importance of epistasis investigations in marker trait association studies, as the individual effect of a marker as locus depends on the marker genotype at other interacting loci. In fact, a favorable allele at one locus may be an unfavorable one in a different genetic background, and vice versa (Holland, 2001). Thus, this has to be taken in consideration, especially for sugarcane, due to the several possible interactions between the multiple alleles from different loci.

Furthermore, the results indicate that epistatic interaction plays an important role in controlling the expression of complex traits. Yu et al. (1997) and Rahman et al. (2007) also identified a number of epistatic QTLs influencing yield and yield components. Thus, the utilization of marker assisted selection in different plant breeding programs has to take epistatic effects into consideration.

Internet Resources

Received: July 16, 2012

Accepted: December 12, 2012.

Associate Editor: Dario Grattapaglia

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

• Bateson W (1909) Mendel's Principles of Heredity. Cambridge University Press, Cambridge, 396 pp.
• Bocianowski J (2008) Comparison of two methods of estimation of nonallelic interaction of QTL effects on the basis of doubled haploid lines in barley. Agric Conspec Sci 73:183-187.
• Bocianowski J (2012a) A comparison of two methods to estimate additive-by-additive interaction of QTL effects by a simulation study. J Theor Biol 308:20-24.
• Bocianowski J (2012b) Analytical and numerical comparisons of two methods of estimation of additive x additive interaction of QTL effects. Sci Agric 69:240-246.
• Bocianowski J (2012c) Using weighted multiple linear regression in estimation of QTL-by-QTL epistasis effects. Genet Mol Biol 35:802-809.
• Bocianowski J and Krajewski P (2009) Comparison of the genetic additive effect estimators based on phenotypic observations and on molecular marker data. Euphytica 165:113-122.
• Charcosset A, Mangin B, Moreau L, Combes L, Jourjon M-F and Gallais A (2001) Heterosis in maize investigated using connected RIL populations. In: Gallais A, Dillmann C and Goldringer I (eds) Quantitative Genetics and Breeding Methods: The Way Ahead. Proceedings of the Eleventh Meeting of the EUCARPIA, Section Biometrics in Plant Breeding, Institut National de la Recherche Agronomique, Paris, France, Les Colloques 96, pp 89-98.
• Cockerham CC and Zeng ZB (1996) Design with marker loci. Genetics 143:1437-1456.
• Conti V, Roncallo PF, Beaufort V, Cervigni GL, Miranda R, Jensen CA and Echenique VC (2011) Mapping of main and epistatic effect QTLs associated to grain protein and gluten strength using a RIL population of durum wheat. J Appl Genet 52:287-298.
• De-Vicente MC and Tanksley SD (1993) QTL analysis of transgressive segregation in an interspecific tomato cross. Genetics 134:585-596.
• Doebley J and Stec A (1991) Genetic analysis of the morphological differences between maize and teosinte. Genetics 129:285-295.
• Edwards MD, Stuber CW and Wendel JF (1987) Molecularmarker-facilitated investigations of quantitative-trait loci in maize. I. Numbers, genomic distribution and types of gene action. Genetics 116:113-125.
• Eshed Y and Zamir D (1996) Less-than-additive epistatic interactions of quantitative trait loci in tomato. Genetics 143:1807- 1817.
• Falconer DS (1989) Introduction to Quantitative Genetics. John Wiley and Sons, New York, 340 pp.
• Govindaraj P, Vinod KK, Arumugachamy S and Maheswaran M (2009) Analysing genetic control of cooked grain traits and gelatinization temperature in a double haploid population of rice by quantitative trait loci mapping. Euphytica 166:165-176.
• Gowda SJM, Radzika P, Mhase LB, Jamadagni BM, Gupta VS and Kadro NY (2011) Mapping of QTLs governing agronomic and field traits in chickpea. J Appl Genet 52:9-21.
• Hayes PM, Liu BH, Knapp SJ, Chen F, Jones B, Blake T, Franckowiak J, Rasmusson D, Sorrells M, Ullrich SE, et al. (1993) Quantitative trait locus effects and environmental interaction in a sample of North American barley germ plasm. Theor Appl Genet 87:392-401.
• Holland J (2001) Epistasis and plant breeding. In: Janick J (ed) Plant Breeding Reviews. Vol. 21. John Wiley & Sons, Inc., New York, pp 27-112.
• Jansen RC and Stam P (1994) High resolution of quantitative traits into multiple loci via interval mapping. Genetics 136:1447-1455.
• Jiang Z, Zhang B, Teng W, Han Y, Zhao X, Sun D, Zhang Z and Li W (2011). Impact of epistasis and QTL x environmental interaction on the oil filling rate of soybean seed at different developmental stages. Euphytica 177:431-442
• Kearsey MJ and Pooni HS (1998) The Genetical Analysis of Quantitative Traits. Chapman and Hall, London, 381 pp.
• Kleinhofs A, Kilian A, Saghai Maroof MA, Biyashev RM, Hayes P, Chen FQ, Lapitan N, Fenwick A, Blake TK, Kanazin V, et al. (1993) A molecular, isozyme and morphological map of the barley (Hordeum vulgare) genome. Theor Appl Genet 86:705-712.
• Krajewski P, Bocianowski J, Gawlowska M, Kaczmarek Z, Pniewski T, Swiecicki W and Wolko B (2012) QTL for field components and protein content: A multienvironmental study of two pea (Pisum sativum L.) populations. Euphytica 183:323-336.
• Lark KG, Chase K, Adler FR, Mansur LM and Orf JJ (1995) Interactions between quantitative trait loci in soybean in which trait variation at one locus is conditional upon a specific allele at another. Proc Natl Acad Sci USA 92:4656-4660.
• Li ZK, Pinson SRM, Park WD, Paterson AH and Stansel JW (1997a) Epistasis for three grain yield components in rice (Oryza sativa L.). Genetics 145:453-465.
• Li ZK, Pinson SRM, Paterson AH, Park WD and Stansel JW (1997b) Genetics of hybrid sterility and hybrid breakdown in an inter-subspecific rice (Oryza sativa L.) population. Genetics 145:1139-1148.
• Li Z-M, Ding J-Q, Wang R-X, Chen J-F, Sun X-D, Chen W, Song W-B, Dong H-F, Dai X-D, Xia Z-L, et al. (2011) A new QTL for resistance to Fusarium ear rot in maize. J Appl Genetics 52:403-406.
• Lin Y, Schertz KF and Paterson AH (1995) Comparative analysis of QTLs affecting plant height and maturity across Poacceae, in reference to an interspecific sorghum population. Genetics 141:391-411.
• Ma XQ, Tang JH, Teng WT, Yan JB, Meng YJ and Li JS (2007) Epistatic interaction is an important genetic basis of grain yield and its components in maize. Mol Breed 20:41-51.
• Mao D, Liu T, Xu C, Li X and Xing Y (2011). Epistasis and complementary gene action adequately account for the genetic bases of transgressive segregation of kilo-grain weight in rice. Euphytica 180:261-271.
• Martinez O and Curnow RN (1994) Missing markers when estimating quantitative trait loci using regression mapping. Heredity 73:198-206.
• Paterson AH, Damon S, Hewitt JD, Zamir D, Rabinowitch HD, Lincoln SE, Lander ES and Tanksley SD (1991) Mendelian factors underlying quantitative traits in tomato: Comparison across species, generations and environments. Genetics 127:181-197.
• Province MA (2001) 30 Sequential methods of analysis for genome scan. Advance in Genetics 42:499-514.
• Rahman ML, Chu SH, Choi MS, Qiao YL, Jiang W, Piao R, Khanam S, Cho YI, Jeung J, Jena KK, et al. (2007) Identification of QTLs for some agronomic traits in rice using an introgression line from Oryza minuta Mol Cells 24:16-26.
• Romagosa I, Ullrich SE, Han F and Hayes PM (1996) Use of the additive main effects and multiplicative interaction model in QTL mapping for adaptation in barley. Theor Appl Genet 93:30-37.
• Stuber CW, Lincoln SE, Wolff DW, Helentjaris T and Lander ES (1992) Identification of genetic factors contributing to heterosis in a hybrid from two elite maize inbred lines using molecular markers. Genetics 132:823-839.
• Tachida H and Cockerham CC (1989) A building block model for quantitative genetics. Genetics 121:839-844.
• Tan YF, Sun M, Xing YZ, Hua JP, Sun XL, Zhang QF and Corke H (2001) Mapping quantitative trait loci for milling quality, protein content and color characteristics of rice using a recombinant inbred line population derived from an elite rice hybrid. Theor Appl Genet 103:1037-1045.
• Tinker NA, Mather DE, Rossnagel BG, Kasha KJ, Kleinhofs A, Hayes PM, Falk DE, Ferguson T, Shugar LP, Legge WG, et al. (1996) Regions of the genome that affect agronomic performance in two-row barley. Crop Sci 36:1053-1062.
• Upadhyaya HD, Sharma S, Singh S and Singh M (2011). Inheritance of drought resistance related traits in two crosses of groundnut (Arachis hypogaea L.). Euphytica 177:55-66.
• Xiao J, Li J, Yuan L and Tanksley SD (1995) Dominance is the major genetic basis of heterosis in rice as revealed by QTL analysis using molecular markers. Genetics 140:745-754.
• Xing YZ, Tan YF, Hua JP, Sun XL and Zhang Q (2002) Characterization of the main effects, epistatic effects and their environmental interactions of QTLs on the genetic basis of yield traits in rice. Theor Appl Genet 105:248-257.
• Yu SB, Li JX, Xu CG, Tan YF, Gao YJ, Li XH, Zhang Q and Maroof MAS (1997) Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci USA 94:9226-9231.
• Zeng Z-B, Wang T and Zou W (2005) Modeling quantitative trait loci and interpretation of models. Genetics 169:1711-1725.
• GrainGenes, a database for Triticeae and Avena: The Steptoe x Morex Barley Mapping Population, http://wheat.pw.usda.gov/ggpages/SxM (July 16, 2012).
• GrainGenes, a database for Triticeae and Avena: Hordeum, http://wheat.pw.usda.gov/ggpages/maps/ Hordeum (July 16, 2012).

Publication Dates

• Publication in this collection
  28 Feb 2013
• Date of issue
  2013

History

• Received
  16 July 2012
• Accepted
  12 Dec 2012