Mapping quantitative trait loci underlying triploid endosperm traits

Mixed linear model approach for mapping quantitative trait loci underlying crop seed traits

The crop seed is a complex organ that may be composed of the diploid embryo, the triploid endosperm and the diploid maternal tissues. According to the genetic features of seed characters, two genetic models for mapping quantitative trait loci (QTLs) of crop seed traits are proposed, with inclusion of maternal effects, embryo or endosperm effects of QTL, environmental effects and QTL-by-environment (QE) interactions. The mapping population can be generated either from double back-cross of immortalized F2 (IF2) to the two parents, from random-cross of IF2 or from selfing of IF2 population. Candidate marker intervals potentially harboring QTLs are first selected through one-dimensional scanning across the whole genome. The selected candidate marker intervals are then included in the model as cofactors to control background genetic effects on the putative QTL(s). Finally, a QTL full model is constructed and model selection is conducted to eliminate false positive QTLs. The genetic main effects of QTLs, QE interaction effects and the corresponding P-values are computed by Markov chain Monte Carlo algorithm for Gaussian mixed linear model via Gibbs sampling. Monte Carlo simulations were performed to investigate the reliability and efficiency of the proposed method. The simulation results showed that the proposed method had higher power to accurately detect simulated QTLs and properly estimated effect of these QTLs. To demonstrate the usefulness, the proposed method was used to identify the QTLs underlying fiber percentage in an upland cotton IF2 population. A computer software, QTLNetwork-Seed, was developed for QTL analysis of seed traits.

Genetic dissection of complex endosperm traits

The endosperm of plants is a major source of food, feed and industrial raw materials. The genetic analysis of endosperm traits poses numerous challenges due to the endosperm's complex genetic composition and unique physical and developmental properties. Modern molecular techniques and statistical methods have greatly improved the mapping of quantitative trait loci underlying endosperm traits and have led to revolutionary insights regarding epistatic and epigenetic effects. This article describes the current state of the methodologies used in the genetic dissection of endosperm traits and highlights practical issues and statistical concepts and procedures.

Improved genetic mapping of endosperm traits using NCIII and TTC designs

The endosperm in plants is a major source of human nutrition and industrial raw material. The genetic study of endosperm poses a great challenge due to its complex genetic composition and unique physical and developmental properties. In this note, we shall revisit 2 classic mating designs-North Carolina Design III (NCIII) and triple test cross (TTC)-and demonstrate their efficiency in detecting quantitative trait loci underlying endosperm traits.

Mapping quantitative trait loci for binary trait in the F2:3 design

In the analysis of inheritance of quantitative traits with low heritability, an F(2:3) design that genotypes plants in F(2) and phenotypes plants in F(2:3) progeny is often used in plant genetics. Although statistical approaches for mapping quantitative trait loci (QTL) in the F(2:3) design have been well developed, those for binary traits of biological interest and economic importance are seldom addressed. In this study, an attempt was made to map binary trait loci (BTL) in the F(2:3) design. The fundamental idea was: the F(2) plants were genotyped, all phenotypic values of each F(2:3) progeny were measured for binary trait, and these binary trait values and the marker genotype informations were used to detect BTL under the penetrance and liability models. The proposed method was verified by a series of Monte-Carlo simulation experiments. These results showed that maximum likelihood approaches under the penetrance and liability models provide accurate estimates for the effects and the locations of BTL with high statistical power, even under of low heritability. Moreover, the penetrance model is as efficient as the liability model, and the F(2:3) design is more efficient than classical F(2) design, even though only a single progeny is collected from each F(2:3) family. With the maximum likelihood approaches under the penetrance and the liability models developed in this study, we can map binary traits as we can do for quantitative trait in the F(2:3) design.

Mapping epistatic quantitative trait loci underlying endosperm traits using all markers on the entire genome in a random hybridization design

Triploid endosperm is of great economic importance owing to its nutritious quality. Mapping endosperm trait loci (ETL) can provide an efficient way to genetically improve grain quality. However, most triploid ETL mapping methods do not produce unbiased estimates of the two dominant effects of ETL. A random hybridization design is an alternative method that may be used to overcome this problem. However, epistasis has an important role in the dissection of genetic architecture for complex traits. In this study, therefore, an attempt was made to map epistatic ETL (eETL) under a triploid genetic model of endosperm traits in a random hybridization design. The endosperm trait means of random hybrid lines, together with known marker genotype information from their corresponding parental F(2) plants, were used to estimate, efficiently and without bias, the positions and all of the effects of eETL using a penalized maximum likelihood method. The method proposed in this article was verified by a series of Monte Carlo simulation experiments. Results from the simulated studies show that the proposed method provides accurate estimates of eETL parameters with a low false-positive rate and a relatively short running time. This new method enables us to map triploid eETL in the same way as diploid quantitative traits.

The QTL analysis on maternal and endosperm genome and their environmental interactions for characters of cooking quality in rice (Oryza sativa L.)

Investigations to identify quantitative trait loci (QTLs) governing cooking quality traits including amylose content, gel consistency and gelatinization temperature (expressed by the alkali spread value) were conducted using a set of 241 RIL populations derived from an elite hybrid cross of "Zhenshan 97"x"Minghui 63" and their reciprocal backcrosses BC1F1 and BC2F1 populations in two environments. QTLs and QTLxenvironment interactions were analyzed by using the genetic model with endosperm and maternal effects and environmental interaction effects on quantitative traits of seed in cereal crops. The results suggested that a total of seven QTLs were associated with cooking quality of rice, which were subsequently mapped to chromosomes 1, 4 and 6. Six of these QTLs were also found to have environmental interaction effects.

Interval Mapping of Quantitative Trait Loci Underlying Triploid Endosperm Traits Using F3 Seeds

A statistical method for mapping quantitative trait loci (QTLs) underlying endosperm traits is proposed. The method is based on a genetic model containing both the direct effects and maternal effects of an endosperm QTL and on an experimental design termed two-stage hierarchical design, in which the trait information is obtained from F(3) endosperms and molecular marker information is obtained from F(2) plants and F(3) embryos (plants). Results of computer simulations indicate that the method can efficiently map endosperm QTLs and precisely estimate both the direct and maternal effects of endosperm QTLs when the sample size is sufficiently large.

A mixture model approach to the mapping of QTL controlling endosperm traits with bulked samples

Endosperm traits are of triploid inheritance and have become a focus of breeding effort for their close relations with the grain quality. Current methods for mapping quantitative trait loci (QTL) underlying endosperm traits are restricted to the use of the phenotypes of single grain samples as input data set, which are often not available in practice due to the small size of the cereal seeds. This paper proposed a statistical model for one specially tailored mapping strategy, where the marker genotypes are obtained from the maternal plants in the segregation population and the phenotypic responses are replaced by the trait means of composite endosperm samples pooled from each plant. It should therefore be more practical and have wide applicability in mapping endosperm traits. The method was implemented by fitting the phenotypic means of endosperms into a Gaussian mixture model. Both the exact and approximate Expectation-Maximization algorithms were proposed to estimate the model parameters. The presence of the QTL was determined by likelihood ratio test statistics. Statistical power and other properties of the new method were investigated and compared to the current single-seed method under a variety of scenarios through simulation studies. The simulations suggest a reasonable sample size should be used to ensure reliable results. The proposed method was also applied to a simulated genome data for further evaluation. As an illustration, a real data of maize was analyzed to find the loci responsible for the popping expansion volume.

Mapping quantitative trait loci using the experimental designs of recombinant inbred populations

In the data collection of the QTL experiments using recombinant inbred (RI) populations, when individuals are genotyped for markers in a population, the trait values (phenotypes) can be obtained from the genotyped individuals (from the same population) or from some progeny of the genotyped individuals (from the different populations). Let Fu be the genotyped population and Fv (v>or=u) be the phenotyped population. The experimental designs that both marker genotypes and phenotypes are recorded on the same populations can be denoted as (Fu/Fv, u=v) designs and that genotypes and phenotypes are obtained from the different populations can be denoted as (Fu/Fv, v>u) designs. Although most of the QTL mapping experiments have been conducted on the backcross and F2(F2/F2) designs, the other (Fu/Fv, v>or=u) designs are also very popular. The great benefits of using the other (Fu/Fv, v>or=u) designs in QTL mapping include reducing cost and environmental variance by phenotyping several progeny for the genotyped individuals and taking advantages of the changes in population structures of other RI populations. Current QTL mapping methods including those for the (Fu/Fv, u=v) designs, mostly for the backcross or F2/F2 design, and for the F2/F3 design based on a one-QTL model are inadequate for the investigation of the mapping properties in the (Fu/Fv, uor=u) designs. In addition, the QTL mapping properties of the proposed and approximate methods in different designs are discussed. Simulations were performed to evaluate the performance of the proposed and approximate methods. The proposed method is proven to be able to correct the problems of the approximate and current methods for improving the resolution of genetic architecture of quantitative traits and can serve as an effective tool to explore the QTL mapping study in the system of RI populations.

A Method for Identification of the Expression Mode and Mapping of QTL Underlying Embryo-Specific Characters

Embryos of crop seeds are one of the major sources of the plant protein and lipid for human nutrition. The genetic expression for embryo-specific characters in crop seeds can be controlled exclusively by the embryo or the maternal genotypes and sometimes by both simultaneously. However, current methods for mapping quantitative trait loci (QTLs) underlying characters of maternal plants have not been effective in dealing with the QTL analysis of embryo characters. On the basis of the expression feature of embryo, a statistical method was proposed for the identification of expression mode and mapping of QTL controlling embryo traits. The maximum likelihood method implemented via the expectation maximization algorithm was used to estimate parameters of a putative embryo-specific QTL. The QTL expression mode was identified by the likelihood ratio test statistic. Statistical power and other properties of the proposed method were investigated under a variety of scenarios through simulation studies. The results showed that the mapping method neglecting the effects of embryo genotype or maternal effects could neither identify the expression mode of QTL nor estimate its genetic effects accurately, whereas the proposed method could effectively map the embryo-specific QTL of various expression modes.

Multiple-interval mapping for quantitative trait loci controlling endosperm traits

Endosperm traits are trisomic inheritant and are of great economic importance because they are usually directly related to grain quality. Mapping for quantitative trait loci (QTL) underlying endosperm traits can provide an efficient way to genetically improve grain quality. As the traditional QTL mapping methods (diploid methods) are usually designed for traits under diploid control, they are not the ideal approaches to map endosperm traits because they ignore the triploid nature of endosperm. In this article, a statistical method considering the triploid nature of endosperm (triploid method) is developed on the basis of multiple-interval mapping (MIM) to map for the underlying QTL. The proposed triploid MIM method is derived to broadly use the marker information either from only the maternal plants or from both the maternal plants and their embryos in the backcross and F2 populations for mapping endosperm traits. Due to the use of multiple intervals simultaneously to take multiple QTL into account, the triploid MIM method can provide better detection power and estimation precision, and as shown in this article it is capable of analyzing and searching for epistatic QTL directly as compared to the traditional diploid methods and current triploid methods using only one (or two) interval(s). Several important issues in endosperm trait mapping, such as the relation and differences between the diploid and triploid methods, variance components of genetic variation, and the problems if effects are present and ignored, are also addressed. Simulations are performed to further explore these issues, to investigate the relative efficiency of different experimental designs, and to evaluate the performance of the proposed and current methods in mapping endosperm traits. The MIM-based triploid method can provide a powerful tool to estimate the genetic architecture of endosperm traits and to assist the marker-assisted selection for the improvement of grain quality in crop science. The triploid MIM FORTRAN program for mapping endosperm traits is available on the worldwide web (http://www.stat.sinica.edu.tw/chkao/).

Mapping genome-genome epistasis: a high-dimensional model

MOTIVATION

The proper development of any organ or tissue requires the coordinated expression of its underlying genes that can be located on different genomes present in an organism. For instance, each step in the development of seed for a higher plant is the consequence of gene interactions from the maternal, embryo and endosperm genomes.

RESULTS

We present a multivariate statistical model for mapping quantitative trait loci (QTL) by incorporating two important aspects of seed development in plants-QTL interactions derived from different genomes, the maternal, embryo and endosperm, and genetic correlations among phenotypic traits expressed in different genome-specific tissues. This model, which has a high dimensionality, is constructed within the maximum-likelihood context based on a finite mixture model. The implementation of the expectation-maximization algorithm allows for the efficient estimation of QTL positions, their action and interaction effects and pleiotropic effects. The application of this high-dimensional model to a real rice dataset has validated its usefulness.

CONCLUSIONS

Our model was derived for self-pollinated plants, but it can be extended to cross-pollinated plants and to animals. With the burgeoning of genetic and genomic data, this high-dimensional model will have many implications for agricultural and evolutionary genetic research.

AVAILABILITY

A package of software will be provided from the corresponding author upon request.

Mapping quantitative trait loci in F2 incorporating phenotypes of F3 progeny

In plants and laboratory animals, QTL mapping is commonly performed using F(2) or BC individuals derived from the cross of two inbred lines. Typical QTL mapping statistics assume that each F(2) individual is genotyped for the markers and phenotyped for the trait. For plant traits with low heritability, it has been suggested to use the average phenotypic values of F(3) progeny derived from selfing F(2) plants in place of the F(2) phenotype itself. All F(3) progeny derived from the same F(2) plant belong to the same F(2:3) family, denoted by F(2:3). If the size of each F(2:3) family (the number of F(3) progeny) is sufficiently large, the average value of the family will represent the genotypic value of the F(2) plant, and thus the power of QTL mapping may be significantly increased. The strategy of using F(2) marker genotypes and F(3) average phenotypes for QTL mapping in plants is quite similar to the daughter design of QTL mapping in dairy cattle. We study the fundamental principle of the plant version of the daughter design and develop a new statistical method to map QTL under this F(2:3) strategy. We also propose to combine both the F(2) phenotypes and the F(2:3) average phenotypes to further increase the power of QTL mapping. The statistical method developed in this study differs from published ones in that the new method fully takes advantage of the mixture distribution for F(2:3) families of heterozygous F(2) plants. Incorporation of this new information has significantly increased the statistical power of QTL detection relative to the classical F(2) design, even if only a single F(3) progeny is collected from each F(2:3) family. The mixture model is developed on the basis of a single-QTL model and implemented via the EM algorithm. Substantial computer simulation was conducted to demonstrate the improved efficiency of the mixture model. Extension of the mixture model to multiple QTL analysis is developed using a Bayesian approach. The computer program performing the Bayesian analysis of the simulated data is available to users for real data analysis.

Mapping Quantitative Trait Loci in F2 Incorporating Phenotypes of F3 Progeny

In plants and laboratory animals, QTL mapping is commonly performed using F2 or BC individuals derived from the cross of two inbred lines. Typical QTL mapping statistics assume that each F2 individual is genotyped for the markers and phenotyped for the trait. For plant traits with low heritability, it has been suggested to use the average phenotypic values of F3 progeny derived from selfing F2 plants in place of the F2 phenotype itself. All F3 progeny derived from the same F2 plant belong to the same F2:3 family, denoted by F2:3. If the size of each F2:3 family (the number of F3 progeny) is sufficiently large, the average value of the family will represent the genotypic value of the F2 plant, and thus the power of QTL mapping may be significantly increased. The strategy of using F2 marker genotypes and F3 average phenotypes for QTL mapping in plants is quite similar to the daughter design of QTL mapping in dairy cattle. We study the fundamental principle of the plant version of the daughter design and develop a new statistical method to map QTL under this F2:3 strategy. We also propose to combine both the F2 phenotypes and the F2:3 average phenotypes to further increase the power of QTL mapping. The statistical method developed in this study differs from published ones in that the new method fully takes advantage of the mixture distribution for F2:3 families of heterozygous F2 plants. Incorporation of this new information has significantly increased the statistical power of QTL detection relative to the classical F2 design, even if only a single F3 progeny is collected from each F2:3 family. The mixture model is developed on the basis of a single-QTL model and implemented via the EM algorithm. Substantial computer simulation was conducted to demonstrate the improved efficiency of the mixture model. Extension of the mixture model to multiple QTL analysis is developed using a Bayesian approach. The computer program performing the Bayesian analysis of the simulated data is available to users for real data analysis.