A Bayesian Nonparametric Approach for Mapping Dynamic Quantitative Traits

In biology, many quantitative traits are dynamic in nature. They can often be described by some smooth functions or curves. A joint analysis of all the repeated measurements of the dynamic traits by functional quantitative trait loci (QTL) mapping methods has the benefits to (1) understand the genetic control of the whole dynamic process of the quantitative traits and (2) improve the statistical power to detect QTL. One crucial issue in functional QTL mapping is how to correctly describe the smoothness of trajectories of functional valued traits. We develop an efficient Bayesian nonparametric multiple-loci procedure for mapping dynamic traits. The method uses the Bayesian P-splines with (nonparametric) B-spline bases to specify the functional form of a QTL trajectory and a random walk prior to automatically determine its degree of smoothness. An efficient deterministic variational Bayes algorithm is used to implement both (1) the search of an optimal subset of QTL among large marker panels and (2) estimation of the genetic effects of the selected QTL changing over time. Our method can be fast even on some large-scale data sets. The advantages of our method are illustrated on both simulated and real data sets.     

A semiparametric approach for composite functional mapping of dynamic quantitative traits

Functional mapping has emerged as a powerful tool for mapping quantitative trait loci (QTL) that control developmental patterns of complex dynamic traits. Original functional mapping has been constructed within the context of simple interval mapping, without consideration of separate multiple linked QTL for a dynamic trait. In this article, we present a statistical framework for mapping QTL that affect dynamic traits by capitalizing on the strengths of functional mapping and composite interval mapping. Within this so-called composite functional-mapping framework, functional mapping models the time-dependent genetic effects of a QTL tested within a marker interval using a biologically meaningful parametric function, whereas composite interval mapping models the time-dependent genetic effects of the markers outside the test interval to control the genome background using a flexible nonparametric approach based on Legendre polynomials. Such a semiparametric framework was formulated by a maximum-likelihood model and implemented with the EM algorithm, allowing for the estimation and the test of the mathematical parameters that define the QTL effects and the regression coefficients of the Legendre polynomials that describe the marker effects. Simulation studies were performed to investigate the statistical behavior of composite functional mapping and compare its advantage in separating multiple linked QTL as compared to functional mapping. We used the new mapping approach to analyze a genetic mapping example in rice, leading to the identification of multiple QTL, some of which are linked on the same chromosome, that control the developmental trajectory of leaf age.     

Rank-based statistical methodologies for quantitative trait locus mapping

This article addresses the identification of genetic loci (QTL and elsewhere) that influence nonnormal quantitative traits with focus on experimental crosses. QTL mapping is typically based on the assumption that the traits follow normal distributions, which may not be true in practice. Model-free tests have been proposed. However, nonparametric estimation of genetic effects has not been studied. We propose an estimation procedure based on the linear rank test statistics. The properties of the new procedure are compared with those of traditional likelihood-based interval mapping and regression interval mapping via simulations and a real data example. The results indicate that the nonparametric method is a competitive alternative to the existing parametric methodologies.     

Bayesian functional mapping of dynamic quantitative traits

Without consideration of other linked QTLs responsible for dynamic trait, original functional mapping based on a single QTL model is not optimal for analyzing multiple dynamic trait loci. Despite that composite functional mapping incorporates the effects of genetic background outside the tested QTL in mapping model, the arbitrary choice of background markers also impact on the power of QTL detection. In this study, we proposed Bayesian functional mapping strategy that can simultaneously identify multiple QTL controlling developmental patterns of dynamic traits over the genome. Our proposed method fits the change of each QTL effect with the time by Legendre polynomial and takes the residual covariance structure into account using the first autoregressive equation. Also, Bayesian shrinkage estimation was employed to estimate the model parameters. Especially, we specify the gamma distribution as the prior for the first-order auto-regressive coefficient, which will guarantee the convergence of Bayesian sampling. Simulations showed that the proposed method could accurately estimate the QTL parameters and had a greater statistical power of QTL detection than the composite functional mapping. A real data analysis of leaf age growth in rice is used for the demonstration of our method. It shows that our Bayesian functional mapping can detect more QTLs as compared to composite functional mapping.     

Functional mapping of dynamic traits with robust t-distribution

Functional mapping has been a powerful tool in mapping quantitative trait loci (QTL) underlying dynamic traits of agricultural or biomedical interest. In functional mapping, multivariate normality is often assumed for the underlying data distribution, partially due to the ease of parameter estimation. The normality assumption however could be easily violated in real applications due to various reasons such as heavy tails or extreme observations. Departure from normality has negative effect on testing power and inference for QTL identification. In this work, we relax the normality assumption and propose a robust multivariate t-distribution mapping framework for QTL identification in functional mapping. Simulation studies show increased mapping power and precision with the t distribution than that of a normal distribution. The utility of the method is demonstrated through a real data analysis.     

家畜抗性等级性状的QTL定位方法

在广义线性模型的框架内模拟研究了家畜抗性等级性状的QTL定位方法,QTL参数的估计采用最大似然方法,比较了阈模型方法与一般线性方法的QTL定位效率,并对影响等级性状QTL定位效率的主要因素（QTL效应、性状的遗传力）进行了模拟研究,实验设计为多个家系的女儿设计,资源群体大小为500头。研究结果表明：在QTL位置参数估计及检验功效方面,阈模型方法具有一定的优势,对抗性等级性状QTL定位的功效也高于线性方法。另外,性状遗传力和QTL效应的大小对QTL定位的准确度也有直接的影响,随着性状遗传力QTL效应的     

High-density genotyping: an overkill for QTL mapping? Lessons learned from a case study in maize and simulations

High-density genotyping is extensively exploited in genome-wide association mapping studies and genomic selection in maize. By contrast, linkage mapping studies were until now mostly based on low-density genetic maps and theoretical results suggested this to be sufficient. This raises the question, if an increase in marker density would be an overkill for linkage mapping in biparental populations, or if important QTL mapping parameters would benefit from it. In this study, we addressed this question using experimental data and a simulation based on linkage maps with marker densities of 1, 2, and 5 cM. QTL mapping was performed for six diverse traits in a biparental population with 204 doubled haploid maize lines and in a simulation study with varying QTL effects and closely linked QTL for different population sizes. Our results showed that high-density maps neither improved the QTL detection power nor the predictive power for the proportion of explained genotypic variance. By contrast, the precision of QTL localization, the precision of effect estimates of detected QTL, especially for small and medium sized QTL, as well as the power to resolve closely linked QTL profited from an increase in marker density from 5 to 1 cM. In conclusion, the higher costs for high-density genotyping are compensated for by more precise estimates of parameters relevant for knowledge-based breeding, thus making an increase in marker density for linkage mapping attractive.     

Estimating the genetic architecture of quantitative traits

Understanding and estimating the structure and parameters associated with the genetic architecture of quantitative traits is a major research focus in quantitative genetics. With the availability of a well-saturated genetic map of molecular markers, it is possible to identify a major part of the structure of the genetic architecture of quantitative traits and to estimate the associated parameters. Multiple interval mapping, which was recently proposed for simultaneously mapping multiple quantitative trait loci (QTL), is well suited to the identification and estimation of the genetic architecture parameters, including the number, genomic positions, effects and interactions of significant QTL and their contribution to the genetic variance. With multiple traits and multiple environments involved in a QTL mapping experiment, pleiotropic effects and QTL by environment interactions can also be estimated. We review the method and discuss issues associated with multiple interval mapping, such as likelihood analysis, model selection, stopping rules and parameter estimation. The potential power and advantages of the method for mapping multiple QTL and estimating the genetic architecture are discussed. We also point out potential problems and difficulties in resolving the details of the genetic architecture as well as other areas that require further investigation. One application of the analysis is to improve genome-wide marker-assisted selection, particularly when the information about epistasis is used for selection with mating.     

远交群体动态性状基因定位的似然分析Ⅱ.模拟研究

以远交设计群体为例,在推导出动态性状基因定位的似然法分析过程的基础上,选择3阶Legendre多项式为子模型,采用Monte-Carlo方法模拟研究了不同个体数、测定日频数、标记密度和QTL遗传贡献率对两种分析方法检测QTL效率的影响。每个因素都取高、中和低3个水平,利用正交设计安排模拟因素试验组合。模拟试验结果表明：高QTL遗传贡献率要比低QTL遗传贡献率的QTL在检测时需要较少个体数和测定日抽样;但不论QTL遗传贡献率多大,300以上的群体大小和5％以上的测定日频数都可以保证足够高的检测效率。个体数和测定日频数对动态性状QTL的分析和检测具有几乎相同的作用,而且相同样本含量条件下两者呈现互补的关系。就某个动态点的QTL检测而言,模拟试验也同时证明：提出的这种以整个动态过程为定位目标的动态性状基因定位方法明显优于传统的逐个动态点的定位分析方法。     

Functional mapping of quantitative trait loci underlying the character process: a theoretical framework

Unlike a character measured at a finite set of landmark points, function-valued traits are those that change as a function of some independent and continuous variable. These traits, also called infinite-dimensional characters, can be described as the character process and include a number of biologically, economically, or biomedically important features, such as growth trajectories, allometric scalings, and norms of reaction. Here we present a new statistical infrastructure for mapping quantitative trait loci (QTL) underlying the character process. This strategy, termed functional mapping, integrates mathematical relationships of different traits or variables within the genetic mapping framework. Logistic mapping proposed in this article can be viewed as an example of functional mapping. Logistic mapping is based on a universal biological law that for each and every living organism growth over time follows an exponential growth curve (e.g., logistic or S-shaped). A maximum-likelihood approach based on a logistic-mixture model, implemented with the EM algorithm, is developed to provide the estimates of QTL positions, QTL effects, and other model parameters responsible for growth trajectories. Logistic mapping displays a tremendous potential to increase the power of QTL detection, the precision of parameter estimation, and the resolution of QTL localization due to the small number of parameters to be estimated, the pleiotropic effect of a QTL on growth, and/or residual correlations of growth at different ages. More importantly, logistic mapping allows for testing numerous biologically important hypotheses concerning the genetic basis of quantitative variation, thus gaining an insight into the critical role of development in shaping plant and animal evolution and domestication. The power of logistic mapping is demonstrated by an example of a forest tree, in which one QTL affecting stem growth processes is detected on a linkage group using our method, whereas it cannot be detected using current methods. The advantages of functional mapping are also discussed.     

The effect of an imprecise map on interval mapping QTLs

The statistical analysis of quantitative trait locus (QTL) experiments relies on the use of a linkage map of the markers genotyped. Such a map is, at best, a good estimate of the true map. Resources might be diverted into developing better marker maps or improved maps become available after the analysis, raising concerns over the original analysis. It is therefore important to understand the sensitivity of QTL analysis to map inaccuracy. We have used simulation methods to investigate the consequences of an incorrect map on the results of a QTL analysis using interval mapping. Backcross data sets were generated with a particular map and then analysed with both the correct map and incorrect maps. If the incorrect maps maintained the true linkage groups (i.e. no markers were incorrectly assigned to another linkage group), the accuracy of the map had little or no impact on the ability to detect QTLs, the true significance levels of the tests or the relative placement of QTLs. When a marker was incorrectly placed on another linkage group, there was a small increase in the level of the test. After adjusting for this increase, there was a decrease in power to detect a QTL near the misplaced marker. This decrease was of a similar magnitude to that found when using a single-marker analysis compared with interval mapping. These results mean that QTL analyses can proceed without the need for very accurate marker maps, and that estimated QTL positions can be translated onto updated maps without the need for reanalysis.     

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 Trait Analysis of Genetic Mapping for Quantitative Trait Loci

We present in this paper models and statistical methods for performing multiple trait analysis on mapping quantitative trait loci (QTL) based on the composite interval mapping method. By taking into account the correlated structure of multiple traits, this joint analysis has several advantages, compared with separate analyses, for mapping QTL, including the expected improvement on the statistical power of the test for QTL and on the precision of parameter estimation. Also this joint analysis provides formal procedures to test a number of biologically interesting hypotheses concerning the nature of genetic correlations between different traits. Among the testing procedures considered are those for joint mapping, pleiotropy, QTL by environment interaction, and pleiotropy vs. close linkage. The test of pleiotropy (one pleiotropic QTL at a genome position) vs. close linkage (multiple nearby nonpleiotropic QTL) can have important implications for our understanding of the nature of genetic correlations between different traits in certain regions of a genome and also for practical applications in animal and plant breeding because one of the major goals in breeding is to break unfavorable linkage. Results of extensive simulation studies are presented to illustrate various properties of the analyses.     

植物QTL定位方法的研究进展

本文系统地介绍了QTL定位的单一标记分析法、区间作图法以及复合区间作图法、混合显性模型的分析方法,概述了一些主要定位方法的分析原理、存在的主要优缺点。单一标记分析法可以采用方差分析、回归分析或似然比检验的方法分析。区间作图法和复合区间作图法是基于两个相邻标记的QTL定位方法,可采用回归分析或最大似然法分析。复合区间作图法在模型中包括了与其他QTL连锁的标记,可以提高作图的精度和效率。混合线性模型的QTL定位方法可以包括复杂的遗传效应及QTL与环境的互作效应,具有更广阔的应用前景。 Abstract:QTL mapping methods are reviewed for single-marker mapping,interval mapping,composite interval mapping,and mixed-model based method.Statistical approaches along with their properties are discussed for the mapping methods.ANOVA,regression method and likelihood ratio test can be applied in single-marker mapping.Interval mapping and composite interval mapping can be conducted,based on two interval markers,by regression method and maximum likelihood method.Since markers linked with other QTLs are include in the model,composite interval mapping is more precision and powerful.Mapping QTL by mixed-model approaches is more applicable when complicated QTL effects as well as QTL by environment interaction are analyzed.     

Semiparametric functional mapping of quantitative trait loci governing long-term HIV dynamics

MOTIVATION: Functional mapping has proven to be powerful for characterizing quantitative trait loci (QTL) that control complex dynamic traits. More recently, functional mapping has been extended to identify the host QTL responsible for HIV dynamics by incorporating a parametric bi-exponential function for earlier stages of viral load trajectories. However, existing functional mapping cannot be used to map long-term HIV dynamics because no mathematical functions are available for later stages of HIV dynamic changes. RESULTS: We derived a statistical model for functional mapping of dynamic QTL through characterizing HIV load trajectories during a long-term period semiparametrically. The new model was constructed within the maximum likelihood framework and implemented with the EM-simplex algorithm. It allows for the test of differences in the genetic control of short- and long-term HIV dynamics and the characterization of the effects of viral-host genome interaction. Extensive simulation studies have been performed to test the statistical behavior of this model. The new model will provide an important tool for genetic and genomic studies of human complex diseases like HIV/AIDS and their pathological progression. AVAILABILITY: Available on request from the corresponding author.     

QTL mapping of grain length in rice (Oryza sativa L.) using chromosome segment substitution lines

Chromosome segment substitution (CSS) lines have the potential for use in QTL fine mapping and map-based cloning. The standard t-test used in the idealized case that each CSS line has a single segment from the donor parent is not suitable for non-idealized CSS lines carrying several substituted segments from the donor parent. In this study, we present a likelihood ratio test based on stepwise regression (RSTEP-LRT) that can be used for QTL mapping in a population consisting of non-idealized CSS lines. Stepwise regression is used to select the most important segments for the trait of interest, and the likelihood ratio test is used to calculate the LOD score of each chromosome segment. This method is statistically equivalent to the standard t-test with idealized CSS lines. To further improve the power of QTL mapping, a method is proposed to decrease multicollinearity among markers (or chromosome segments). QTL mapping with an example CSS population in rice consisting of 65 non-idealized CSS lines and 82 chromosome segments indicated that a total of 18 segments on eight of the 12 rice chromosomes harboured QTLs affecting grain length under the LOD threshold of 2.5. Three major stable QTLs were detected in all eight environments. Some minor QTLs were not detected in all environments, but they could increase or decrease the grain length constantly. These minor genes are also useful in marker-assisted gene pyramiding.     

Multiple interval mapping for quantitative trait loci.

A new statistical method for mapping quantitative trait loci (QTL), called multiple interval mapping (MIM), is presented. It uses multiple marker intervals simultaneously to fit multiple putative QTL directly in the model for mapping QTL. The MIM model is based on Cockerham's model for interpreting genetic parameters and the method of maximum likelihood for estimating genetic parameters. With the MIM approach, the precision and power of QTL mapping could be improved. Also, epistasis between QTL, genotypic values of individuals, and heritabilities of quantitative traits can be readily estimated and analyzed. Using the MIM model, a stepwise selection procedure with likelihood ratio test statistic as a criterion is proposed to identify QTL. This MIM method was applied to a mapping data set of radiata pine on three traits: brown cone number, tree diameter, and branch quality scores. Based on the MIM result, seven, six, and five QTL were detected for the three traits, respectively. The detected QTL individually contributed from approximately 1 to 27% of the total genetic variation. Significant epistasis between four pairs of QTL in two traits was detected, and the four pairs of QTL contributed approximately 10.38 and 14.14% of the total genetic variation. The asymptotic variances of QTL positions and effects were also provided to construct the confidence intervals. The estimated heritabilities were 0.5606, 0.5226, and 0. 3630 for the three traits, respectively. With the estimated QTL effects and positions, the best strategy of marker-assisted selection for trait improvement for a specific purpose and requirement can be explored. The MIM FORTRAN program is available on the worldwide web (http://www.stat.sinica.edu.tw/chkao/).     

Mapping and analysis of quantitative trait loci in experimental populations

Simple statistical methods for the study of quantitative trait loci (QTL), such as analysis of variance, have given way to methods that involve several markers and high-resolution genetic maps. As a result, the mapping community has been provided with statistical and computational tools that have much greater power than ever before for studying and locating multiple and interacting QTL. Apart from their immediate practical applications, the lessons learnt from this evolution of QTL methodology might also be generally relevant to other types of functional genomics approach that are aimed at the dissection of complex phenotypes, such as microarray assessment of gene expression.     

Semiparametric variance components models for genetic studies with longitudinal phenotypes

In a family-based genetic study such as the Framingham Heart Study (FHS), longitudinal trait measurements are recorded on subjects collected from families. Observations on subjects from the same family are correlated due to shared genetic composition or environmental factors such as diet. The data have a 3-level structure with measurements nested in subjects and subjects nested in families. We propose a semiparametric variance components model to describe phenotype observed at a time point as the sum of a nonparametric population mean function, a nonparametric random quantitative trait locus (QTL) effect, a shared environmental effect, a residual random polygenic effect and measurement error. One feature of the model is that we do not assume a parametric functional form of the age-dependent QTL effect, and we use penalized spline-based method to fit the model. We obtain nonparametric estimation of the QTL heritability defined as the ratio of the QTL variance to the total phenotypic variance. We use simulation studies to investigate performance of the proposed methods and apply these methods to the FHS systolic blood pressure data to estimate age-specific QTL effect at 62cM on chromosome 17.     

The Statistical Power of Inclusive Composite Interval Mapping in Detecting Digenic Epistasis Showing Common F2 Segregation Ratios

Epistasis is a commonly observed genetic phenomenon and an important source of variation of complex traits,which could maintain additive variance and therefore assure the long-term genetic gain in breeding.Inclusive composite interval mapping（ICIM） is able to identify epistatic quantitative trait loci（QTLs） no matter whether the two interacting QTLs have any additive effects.In this article,we conducted a simulation study to evaluate detection power and false discovery rate（FDR） of ICIM epistatic mapping,by considering F₂ and doubled haploid（DH） populations,different F₂ segregation ratios and population sizes.Results indicated that estimations of QTL locations and effects were unbiased,and the detection power of epistatic mapping was largely affected by population size,heritability of epistasis,and the amount and distribution of genetic effects.When the same likelihood of odd（LOD） threshold was used,detection power of QTL was higher in F₂ population than power in DH population;meanwhile FDR in F₂ was also higher than that in DH.The increase of marker density from 10 cM to 5 cM led to similar detection power but higher FDR.In simulated populations,ICIM achieved better mapping results than multiple interval mapping（MIM） in estimation of QTL positions and effect.At the end,we gave epistatic mapping results of ICIM in one actual population in rice（Oryza sativa L.）.