A blocking Gibbs sampling method to detect major genes with phenotypic data from a diallel mating

Diallel mating is a frequently used design for estimating the additive and dominance genetic (polygenic) effects involved in quantitative traits observed in the half- and full-sib progenies generated in plant breeding programmes. Gibbs sampling has been used for making statistical inferences for a mixed-inheritance model (MIM) that includes both major genes and polygenes. However, using this approach it has not been possible to incorporate the genetic properties of major genes with the additive and dominance polygenic effects in a diallel mating population. A parent block Gibbs sampling method was developed in this study to make statistical inferences about the major gene and polygenic effects on quantitative traits for progenies derived from a half-diallel mating design. Using simulated data sets with different major and polygenic effects, the proposed method accurately estimated the major and polygenic effects of quantitative traits, and possible genotypes of parents and progenies. The impact of specifying different prior distributions was examined and was found to have little effect on inference on the posterior distribution. This approach was applied to an experimental data set of Loblolly pine (Pinus taeda L.) derived from a 6-parent half-diallel mating. The result indicated that there might be a recessive major gene affecting height growth in this diallel population.     

Multivariate segregation analysis for quantitative traits in line crosses

Segregation analysis is a method of detecting major genes for quantitative traits without using marker information. It serves as an important tool in helping investigators to plan further studies such as quantitative trait loci mapping or more sophisticated genomic analyses. However, current methods of segregation analysis for a single trait typically have low statistical power. We propose a multivariate segregation analysis (MSA) that takes advantage of the correlation structure of multiple quantitative traits to detect major genes. This method not only increases the statistical power, but allows dissection of the genetic architecture underlying the trait complex. In MSA the observed phenotypes of multiple correlated traits are fitted to a multivariate Gaussian mixture model. Model parameters are estimated under the maximum likelihood framework via the expectation-maximization algorithm. The presence of major genes is tested using likelihood ratio test statistics. Pleiotropy is distinguished from close linkage by comparing three possible models using the Bayesian information criterion. Two simulation experiments were performed based on the F(2) mating design. In the first, the statistical properties of MSA under varying heritabilities and sample sizes were investigated and the results compared with those obtained from single-trait analysis. In the second simulation the efficacy of MSA in separating pleiotropy from close linkage was demonstrated. Finally, the new method was applied to real data and detected a major gene responsible for both plant height and tiller number in rice.     

A quantitative genetic model for analyzing species differences in outcrossing species

A genetic model based on a two-level intra- and interspecific mating design is proposed to estimate the genetic architecture of species differences and heterosis for outcrossing species. The underlying genetic analyses make use of classical quantitative genetic theories and recent results from molecular genetic studies. Gene effects across different quantitative trait loci (QTL) can be approximated by a geometric series. Under natural selection, gene effects are often associated with allele frequencies in a particular way, which can be approximated by the gamma distribution. By incorporating these approximations into family structural analyses in the mating design, we are able to estimate a number of genetic parameters that contribute to quantitative genetic variation based on a nonlinear optimization approach. These parameters include the number of QTL, their gene effects, and their allele frequencies in the parental populations. We perform simulation studies and illustrate an example to demonstrate the statistical property and procedure of the method.     

Multifactorial genetic models for quantitative traits in humans.

Quantitative traits measured in human families can be analyzed to partition the total population variance into genetic and environmental components, or to elucidate the genetic mechanism involved. We review the estimation of variance components directly from human pedigree data, or in the form of path coefficients from correlations between pairs of relatives. To elucidate genetic mechanisms, a mixed model that allows for segregation at a major locus, a polygenic effect and a sibling environmental correlation is described for nuclear families. In each case appropriate likelihoods are derived as a basis, using numerical maximum likelihood methods, for parameter estimation and hypothesis testing. A general model is then described that allows for several familial sources of environmental variation, assortative mating, and both major gene and polygenic effects; and an algorithm for calculating the likelihood of a pedigree under this model is indicated. Finally, some of the remaining problems in this area of biometric analysis are pointed out.     

QTL detection power of multi-parental RIL populations in Arabidopsis thaliana

A major goal of today's biology is to understand the genetic basis of quantitative traits. This can be achieved by statistical methods that evaluate the association between molecular marker variation and phenotypic variation in different types of mapping populations. The objective of this work was to evaluate the statistical power of quantitative trait loci (QTL) detection of various multi-parental mating designs, as well as to assess the reasons for the observed differences. Our study was based on an empirical data of 20 Arabidopsis thaliana accessions, which have been selected to capture the maximum genetic diversity. The examined mating designs differed strongly with respect to the statistical power to detect QTL. We observed the highest power to detect QTL for the diallel cross with random mating design. The results of our study suggested that performing sibling mating within subpopulations of joint-linkage mapping populations has the potential to considerably increase the power for QTL detection. Our results, however, revealed that using designs in which more than two parental alleles segregate in each subpopulation increases the power even more.     

Estimating effects of a single gene and polygenes on quantitative traits from a diallel design

A genetic model is developed with additive and dominance effects of a single gene and polygenes as well as general and specific reciprocal effects for the progeny from a diallel mating design. The methods of ANOVA, minimum norm quadratic unbiased estimation (MINQUE), restricted maximum likelihood estimation (REML), and maximum likelihood estimation (ML) are suggested for estimating variance components, and the methods of generalized least squares (GLS) and ordinary least squares (OLS) for fixed effects, while best linear unbiased prediction, linear unbiased prediction (LUP), and adjusted unbiased prediction are suggested for analyzing random effects. Monte Carlo simulations were conducted to evaluate the unbiasedness and efficiency of statistical methods involving two diallel designs with commonly used sample sizes, 6 and 8 parents, with no and missing crosses, respectively. Simulation results show that GLS and OLS are almost equally efficient for estimation of fixed effects, while MINQUE (1) and REML are better estimators of the variance components and LUP is most practical method for prediction of random effects. Data from a Drosophila melanogaster experiment (Gilbert 1985a, Theor appl Genet 69:625–629) were used as a working example to demonstrate the statistical analysis. The new methodology is also applicable to screening candidate gene(s) and to other mating designs with multiple parents, such as nested (NC Design I) and factorial (NC Design II) designs. Moreover, this methodology can serve as a guide to develop new methods for detecting indiscernible major genes and mapping quantitative trait loci based on mixture distribution theory. The computer program for the methods suggested in this article is freely available from the authors.     

Multilocus genomics of outcrossing plant populations

The structure and organization of natural plant populations can be understood by estimating the genetic parameters related to mating behavior, recombination frequency, and gene associations with DNA-based markers typed throughout the genome. We developed a statistical and computational model for estimating and testing these parameters from multilocus data collected in a natural population. This model, constructed by a maximum likelihood approach and implemented within the EM algorithm, is shown to be robust for simultaneously estimating the outcrossing rate, recombination frequencies and linkage disequilibria. The algorithm built with three or more markers allows the characterization of crossover interference in meiosis and high-order disequilibria among different genes, thus providing a powerful tool for illustrating a detailed picture of genetic diversity and organization in natural populations. Computer simulations demonstrate the statistical properties of the proposed model. This multilocus model will be useful for studying the pattern and amount of genetic variation within and among populations to further infer the evolutionary history of a plant species.     

Selecting a minimal number of relevant genes from microarray data to design accurate tissue classifiers

It is essential to select a minimal number of relevant genes from microarray data while maximizing classification accuracy for the development of inexpensive diagnostic tests. However, it is intractable to simultaneously optimize gene selection and classification accuracy that is a large parameter optimization problem. We propose an efficient evolutionary approach to gene selection from microarray data which can be combined with the optimal design of various multiclass classifiers. The proposed method (named GeneSelect) consists of three parts which are fully cooperated: an efficient encoding scheme of candidate solutions, a generalized fitness function, and an intelligent genetic algorithm (IGA). An existing hybrid approach based on genetic algorithm and maximum likelihood classification (GA/MLHD) is proposed to select a small number of relevant genes for accurate classification of samples. To evaluate the performance of GeneSelect, the gene selection is combined with the same maximum likelihood classification (named IGA/MLHD) for convenient comparisons. The performance of IGA/MLHD is applied to 11 cancer-related human gene expression datasets. The simulation results show that IGA/MLHD is superior to GA/MLHD in terms of the number of selected genes, classification accuracy, and robustness of selected genes and accuracy.     

A score test for the linkage analysis of qualitative and quantitative traits based on identity by descent data from sib-pairs

We propose a general likelihood-based approach to the linkage analysis of qualitative and quantitative traits using identity by descent (IBD) data from sib-pairs. We consider the likelihood of IBD data conditional on phenotypes and test the null hypothesis of no linkage between a marker locus and a gene influencing the trait using a score test in the recombination fraction theta between the two loci. This method unifies the linkage analysis of qualitative and quantitative traits into a single inferential framework, yielding a simple and intuitive test statistic. Conditioning on phenotypes avoids unrealistic random sampling assumptions and allows sib-pairs from differing ascertainment mechanisms to be incorporated into a single likelihood analysis. In particular, it allows the selection of sib-pairs based on their trait values and the analysis of only those pairs having the most informative phenotypes. The score test is based on the full likelihood, i.e. the likelihood based on all phenotype data rather than just differences of sib-pair phenotypes. Considering only phenotype differences, as in Haseman and Elston (1972) and Kruglyak and Lander (1995), may result in important losses in power. The linkage score test is derived under general genetic models for the trait, which may include multiple unlinked genes. Population genetic assumptions, such as random mating or linkage equilibrium at the trait loci, are not required. This score test is thus particularly promising for the analysis of complex human traits. The score statistic readily extends to accommodate incomplete IBD data at the test locus, by using the hidden Markov model implemented in the programs MAPMAKER/SIBS and GENEHUNTER (Kruglyak and Lander, 1995; Kruglyak et al., 1996). Preliminary simulation studies indicate that the linkage score test generally matches or outperforms the Haseman-Elston test, the largest gains in power being for selected samples of sib-pairs with extreme phenotypes.     

两个位点主基因控制的质量—数量性状的遗传分析

应用极大似然法和EM算法提出了关于两个位点主基因控制的质量．数量性状的遗传分析方法,参照质量性状两位点互作在F_2代的分离比率建立了7种遗传假设及其似然比测验的程序,讨论了应用这一方法时应注意的几个问题．     

Spatial auto-correlation and auto-regressive models estimation from sample survey data

Maximum likelihood estimation of the model parameters for a spatial population based on data collected from a survey sample is usually straightforward when sampling and non-response are both non-informative, since the model can then usually be fitted using the available sample data, and no allowance is necessary for the fact that only a part of the population has been observed. Although for many regression models this naive strategy yields consistent estimates, this is not the case for some models, such as spatial auto-regressive models. In this paper, we show that for a broad class of such models, a maximum marginal likelihood approach that uses both sample and population data leads to more efficient estimates since it uses spatial information from sampled as well as non-sampled units. Extensive simulation experiments based on two well-known data sets are used to assess the impact of the spatial sampling design, the auto-correlation parameter and the sample size on the performance of this approach. When compared to some widely used methods that use only sample data, the results from these experiments show that the maximum marginal likelihood approach is much more precise.     

A quantitative genetic model for mixed diploid and triploid hybrid progenies in tree breeding and evolution

Interspecific hybridization has played a critical role in tree evolution and breeding. The findings of triploidy in forest trees stimulate the development of a quantitative genetic model to estimate the nature of gene action. The model is based on clonally replicated triploid progenies derived from a two-level population and individual-within-population mating design in which offspring have a double dose of alleles from the parent and a single dose of alleles from the other parent. With the same genetic assumptions of a diploid model, except non-Mendelian behavior at meiosis, and the experimental variances estimated from a linear statistical model, total genetic variances in the triploid progenies are separated into additive, dominance, and epistatic components. In addition, by combining the new model with the already existing model based on disomic expression, the partitioning of additive, dominant, and epistatic variances can be obtained for a mixed diploid/triploid F₁ progeny population. This paper provides an alternative technique to study the modes of quantitative inheritance for outcrossing, long-lived forest trees in which inbred lines cannot be easily generated. The accuracy for estimating gene action using this technique is discussed.     

Inheritance of zingiberene in Lycopersicon

Summary The simple mating designs provide unbiased estimates for genetic components of variance (additive genetic variance and dominance variance) under the assumption of no epistatic effect. There is empirical evidence, however, that suggests the existence of epistatic gene effects. The triallel and double cross mating designs permit the estimation of epistatic gene effects. A systematic and mathematical approach is suggested for the estimation of variance components based on the alternate model for triallel mating design.     

Directional genetic differentiation and relative migration

Understanding the population structure and patterns of gene flow within species is of fundamental importance to the study of evolution. In the fields of population and evolutionary genetics, measures of genetic differentiation are commonly used to gather this information. One potential caveat is that these measures assume gene flow to be symmetric. However, asymmetric gene flow is common in nature, especially in systems driven by physical processes such as wind or water currents. As information about levels of asymmetric gene flow among populations is essential for the correct interpretation of the distribution of contemporary genetic diversity within species, this should not be overlooked. To obtain information on asymmetric migration patterns from genetic data, complex models based on maximum‐likelihood or Bayesian approaches generally need to be employed, often at great computational cost. Here, a new simpler and more efficient approach for understanding gene flow patterns is presented. This approach allows the estimation of directional components of genetic divergence between pairs of populations at low computational effort, using any of the classical or modern measures of genetic differentiation. These directional measures of genetic differentiation can further be used to calculate directional relative migration and to detect asymmetries in gene flow patterns. This can be done in a user‐friendly web application called divMigrate‐online introduced in this study. Using simulated data sets with known gene flow regimes, we demonstrate that the method is capable of resolving complex migration patterns under a range of study designs.     

An analysis of genetic differentiation among assortatively mating Drosophila melanogaster in Zimbabwe

African Drosophila melanogaster populations, and those from Zimbabwe in particular, have attracted much interest recently. African flies differ genetically from 'cosmopolitan' populations and were found to exhibit discriminative mating behaviour against individuals from 'cosmopolitan' populations. It has therefore been proposed that Zimbabwean and some other African populations are in an 'incipient stage of speciation'. However, whether the mating behaviour is an effective barrier against gene flow from other populations, and whether intra-population genetic differentiation has already evolved in sympatry is not known. Here, we took a population-based approach to test whether the well-characterized mating type differences have resulted in a genome-wide differentiation at the population level. Using 122 polymorphic microsatellite loci mapping to the third chromosome, we demonstrate a significant genetic differentiation between Zimbabwean flies differing in their mating behaviour. We also provide evidence to suggest that this difference is unlikely to be attributable to population structure within Zimbabwe. However, the analysis of individual microsatellite loci did not indicate more loci differentiating these two groups than expected by chance. Our data suggest that the 'Z'-'M' mating behaviour is strong enough to result in a small but significant genetic differentiation. Thus, future studies based on a larger population sample of flies characterized for their mating behaviour and using more markers are expected to provide more information on the genetic basis of the mating traits in the Zimbabwe flies.     

Mapping genome-wide QTL of ratio traits with Bayesian shrinkage analysis for its component traits

The ratio trait is defined as a ratio of two regular quantitative traits with normal distribution, which is distinguished from regular quantitative traits in the genetic analysis because it does not follow the normal distribution. On the basis of maximum likelihood method that uses a special linear combination of the two component traits, we develop a Bayesian mapping strategy for ratio traits, which firstly analyzes the two component traits by Bayesian shrinkage method, and then generates a new posterior sample of genetic effects for a ratio trait from ones of population means and genetic effects for the two component traits, finally, infers QTL for the ratio trait via post MCMC analysis for the new posterior sample. A simulation study demonstrates that the new method has higher detecting power of the QTL than maximum likelihood method. An application is illustrated to map genome-wide QTL for relative growth rate of height on soybean.     

EVOLUTION AND EXTINCTION IN A CHANGING ENVIRONMENT: A QUANTITATIVE-GENETIC ANALYSIS

Because of the ubiquity of genetic variation for quantitative traits, virtually all populations have some capacity to respond evolutionarily to selective challenges. However, natural selection imposes demographic costs on a population, and if these costs are sufficiently large, the likelihood of extinction will be high. We consider how the mean time to extinction depends on selective pressures (rate and stochasticity of environmental change, and strength of selection), population parameters (carrying capacity, and reproductive capacity), and genetics (rate of polygenic mutation). We assume that in a randomly mating, finite population subject to density-dependent population growth, individual fitness is determined by a single quantitative-genetic character under Gaussian stabilizing selection with the optimum phenotype exhibiting directional change, or random fluctuations, or both. The quantitative trait is determined by a finite number of freely recombining, mutationally equivalent, additive loci. The dynamics of evolution and extinction are investigated, assuming that the population is initially under mutation-selection-drift balance. Under this model, in a directionally changing environment, the mean phenotype lags behind the optimum, but on the average evolves parallel to it. The magnitude of the lag determines the vulnerability to extinction. In finite populations, stochastic variation in the genetic variance can be quite pronounced, and bottlenecks in the genetic variance temporarily can impair the population's adaptive capacity enough to cause extinction when it would otherwise be unlikely in an effectively infinite population. We find that maximum sustainable rates of evolution or, equivalently, critical rates of environmental change, may be considerably less than 10% of a phenotypic standard deviation per generation.     

Partitioning of population genetic variance under multiplicative-epistatic gene action

Determining the way in which different QTLs interact (epistasis) in their effects on the phenotype is crucial to many areas in population genetics and evolutionary biology. For example, in the founder event, a separated population readapts to a new environment through the release of cryptic gene-gene interactions. In hybrid zones, hybrid speciation must be subjected to natural selection for epistasis resulting from genomic recombinations between different species. However, there is a severe shortage of relevant methodologies to estimate epistatic genetic effects and variances. A statistical model has recently been proposed to estimate the number of QTLs, their genetic effects and allelic frequencies in segregating populations. This model is based on multiplicative gene action and derived from a two-level intra- and interspecific mating design. In this paper, we formulate a statistical procedure for partitioning the genetic variance into additive, dominant and various kinds of epistatic components in an intra- or mixed intra- and interspecific hybrid population. The procedure can be used to study the genetic architecture of fragmented populations and hybrid zones, thus allowing for a better recognition of the role of epistasis in evolution and hybrid speciation. A real example for two Populus species, P. tremuloides and P. tremula, is provided to illustrate the procedure. In this example, we found that considerable new genetic variation is formed through genomic recombination between two aspen species. Received: 1 May 1999 / Accepted: 27 July 1999     

A non-stationary model for functional mapping of complex traits

SUMMARY: Understanding the genetic control of growth is fundamental to agricultural, evolutionary and biomedical genetic research. In this article, we present a statistical model for mapping quantitative trait loci (QTL) that are responsible for genetic differences in growth trajectories during ontogenetic development. This model is derived within the maximum likelihood context, implemented with the expectation-maximization algorithm. We incorporate mathematical aspects of growth processes to model the mean vector and structured antedependence models to approximate time-dependent covariance matrices for longitudinal traits. Our model has been employed to map QTL that affect body mass growth trajectories in both male and female mice of an F2 population derived from the Large and Small mouse strains. The results from this model are compared with those from the autoregressive-based functional mapping approach. Based on results from computer simulation studies, we suggest that these two models are alternative to one another and should be used simultaneously for the same dataset.