Functional mapping of quantitative trait loci underlying growth trajectories using a transform-both-sides logistic model

The incorporation of developmental control mechanisms of growth has proven to be a powerful tool in mapping quantitative trait loci (QTL) underlying growth trajectories. A theoretical framework for implementing a QTL mapping strategy with growth laws has been established. This framework can be generalized to an arbitrary number of time points, where growth is measured, and becomes computationally more tractable, when the assumption of variance stationarity is made. In practice, however, this assumption is likely to be violated for age-specific growth traits due to a scale effect. In this article, we present a new statistical model for mapping growth QTL, which also addresses the problem of variance stationarity, by using a transform-both-sides (TBS) model advocated by Carroll and Ruppert (1984, Journal of the American Statistical Association 79, 321-328). The TBS-based model for mapping growth QTL cannot only maintain the original biological properties of a growth model, but also can increase the accuracy and precision of parameter estimation and the power to detect a QTL responsible for growth differentiation. Using the TBS-based model, we successfully map a QTL governing growth trajectories to a linkage group in an example of forest trees. The statistical and biological properties of the estimates of this growth QTL position and effect are investigated using Monte Carlo simulation studies. The implications of our model for understanding the genetic architecture of growth are discussed.     

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.     

Functional mapping of growth and development

Understanding how an organism develops into a fully functioning adult from a mass of undifferentiated cells may reveal different strategies that allow the organism to survive under limiting conditions. Here, we review an analytical model for characterizing quantitative trait loci (QTLs) that underlie variation in growth trajectories and developmental timing. This model, called functional mapping, incorporates fundamental principles behind biological processes or networks that are bridged with mathematical functions into a statistical mapping framework. Functional mapping estimates parameters that determine the shape and function of a particular biological process, thus providing a flexible platform to test biologically meaningful hypotheses regarding the complex relationships between gene action and development.     

A general framework for analyzing the genetic architecture of developmental characteristics

The genetic architecture of growth traits plays a central role in shaping the growth, development, and evolution of organisms. While a limited number of models have been devised to estimate genetic effects on complex phenotypes, no model has been available to examine how gene actions and interactions alter the ontogenetic development of an organism and transform the altered ontogeny into descendants. In this article, we present a novel statistical model for mapping quantitative trait loci (QTL) determining the developmental process of complex traits. Our model is constructed within the traditional maximum-likelihood framework implemented with the EM algorithm. We employ biologically meaningful growth curve equations to model time-specific expected genetic values and the AR(1) model to structure the residual variance-covariance matrix among different time points. Because of a reduced number of parameters being estimated and the incorporation of biological principles, the new model displays increased statistical power to detect QTL exerting an effect on the shape of ontogenetic growth and development. The model allows for the tests of a number of biological hypotheses regarding the role of epistasis in determining biological growth, form, and shape and for the resolution of developmental problems at the interface with evolution. Using our newly developed model, we have successfully detected significant additive x additive epistatic effects on stem height growth trajectories in a forest tree.     

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.     

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 - how to map and study the genetic architecture of dynamic complex traits

The development of any organism is a complex dynamic process that is controlled by a network of genes as well as by environmental factors. Traditional mapping approaches for analysing phenotypic data measured at a single time point are too simple to reveal the genetic control of developmental processes. A general statistical mapping framework, called functional mapping, has been proposed to characterize, in a single step, the quantitative trait loci (QTLs) or nucleotides (QTNs) that underlie a complex dynamic trait. Functional mapping estimates mathematical parameters that describe the developmental mechanisms of trait formation and expression for each QTL or QTN. The approach provides a useful quantitative and testable framework for assessing the interplay between gene actions or interactions and developmental changes.     

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.     

Functional mapping imprinted quantitative trait loci underlying developmental characteristics

Background

Genomic imprinting, a phenomenon referring to nonequivalent expression of alleles depending on their parental origins, has been widely observed in nature. It has been shown recently that the epigenetic modification of an imprinted gene can be detected through a genetic mapping approach. Such an approach is developed based on traditional quantitative trait loci (QTL) mapping focusing on single trait analysis. Recent studies have shown that most imprinted genes in mammals play an important role in controlling embryonic growth and post-natal development. For a developmental character such as growth, current approach is less efficient in dissecting the dynamic genetic effect of imprinted genes during individual ontology.

Results

Functional mapping has been emerging as a powerful framework for mapping quantitative trait loci underlying complex traits showing developmental characteristics. To understand the genetic architecture of dynamic imprinted traits, we propose a mapping strategy by integrating the functional mapping approach with genomic imprinting. We demonstrate the approach through mapping imprinted QTL controlling growth trajectories in an inbred F₂ population. The statistical behavior of the approach is shown through simulation studies, in which the parameters can be estimated with reasonable precision under different simulation scenarios. The utility of the approach is illustrated through real data analysis in an F₂ family derived from LG/J and SM/J mouse stains. Three maternally imprinted QTLs are identified as regulating the growth trajectory of mouse body weight.

Conclusion

The functional iQTL mapping approach developed here provides a quantitative and testable framework for assessing the interplay between imprinted genes and a developmental process, and will have important implications for elucidating the genetic architecture of imprinted traits.     

A dynamic model for functional mapping of biological rhythms

Functional mapping is a statistical method for mapping quantitative trait loci (QTLs) that regulate the dynamic pattern of a biological trait. This method integrates mathematical aspects of biological complexity into a mixture model for genetic mapping and tests the genetic effects of QTLs by comparing genotype-specific curve parameters. As a way of quantitatively specifying the dynamic behavior of a system, differential equations have proven to be powerful for modeling and unraveling the biochemical, molecular, and cellular mechanisms of a biological process, such as biological rhythms. The equipment of functional mapping with biologically meaningful differential equations provides new insights into the genetic control of any dynamic processes. We formulate a new functional mapping framework for a dynamic biological rhythm by incorporating a group of ordinary differential equations (ODE). The Runge-Kutta fourth order algorithm was implemented to estimate the parameters that define the system of ODE. The new model will find its implications for understanding the interplay between gene interactions and developmental pathways in complex biological rhythms.     

Bayesian inference for genomic imprinting underlying developmental characteristics

The identification of imprinted genes is becoming a standard procedure in searching for quantitative trait loci (QTL) underlying complex traits. When a developmental characteristic such as growth or drug response is observed at multiple time points, understanding the dynamics of gene function governing the underlying feature should provide more biological information regarding the genetic control of an organism. Recognizing that differential imprinting can be development-specific, mapping imprinted genes considering the dynamic imprinting effect can provide additional biological insights into the epigenetic control of a complex trait. In this study, we proposed a Bayesian imprinted QTL (iQTL) mapping framework considering the dynamics of imprinting effects and model multiple iQTLs with an efficient Bayesian model selection procedure. The method overcomes the limitation of likelihood-based mapping procedure, and can simultaneously identify multiple iQTLs with different gene action modes across the whole genome with high computational efficiency. An inference procedure using Bayes factors to distinguish different imprinting patterns of iQTL was proposed. Monte Carlo simulations were conducted to evaluate the performance of the method. The utility of the approach was illustrated through an analysis of a body weight growth data set in an F(2) family derived from LG/J and SM/J mouse stains. The proposed Bayesian mapping method provides an efficient and computationally feasible framework for genome-wide multiple iQTL inference with complex developmental traits.     

A mechanistic model for genetic machinery of ontogenetic growth

Two different genetic mechanisms can be proposed to explain variation in growth trajectories. The allelic sensitivity hypothesis states that growth trajectory is controlled by the time-dependent expression of alleles at the deterministic quantitative trait loci (dQTL) formed during embryogenesis. The gene regulation hypothesis states that the differentiation in growth process is due to the opportunistic quantitative trait loci (oQTL) through their mediation with new developmental signals. These two hypotheses of genetic control have been elucidated in the literature. Here, we propose a new statistical model for discerning these two mechanisms in the context of growth trajectories by integrating growth laws within a QTL-mapping framework. This model is developed within the maximum-likelihood context, implemented with a grid approach for estimating the genomic positions of the deterministic and opportunistic QTL and the simplex algorithm for estimating the growth curve parameters of the genotypes at these QTL and the parameters modeling the residual (co)variance matrix. Our model allows for extensive hypothesis tests for the genetic control of growth processes and developmental events by these two types of QTL. The application of this new model to an F(2) progeny in mice leads to the detection of deterministic and opportunistic QTL on chromosome 1 for mouse body mass growth. The estimates of QTL positions and effects from our model are broadly in agreement with those by traditional interval-mapping approaches. The implications of this model for biological and biomedical research are discussed.     

Nonparametric functional mapping of quantitative trait loci

Summary .  Functional mapping is a useful tool for mapping quantitative trait loci (QTL) that control dynamic traits. It incorporates mathematical aspects of biological processes into the mixture model-based likelihood setting for QTL mapping, thus increasing the power of QTL detection and the precision of parameter estimation. However, in many situations there is no obvious functional form and, in such cases, this strategy will not be optimal. Here we propose to use nonparametric function estimation, typically implemented with B-splines, to estimate the underlying functional form of phenotypic trajectories, and then construct a nonparametric test to find evidence of existing QTL. Using the representation of a nonparametric regression as a mixed model, the final test statistic is a likelihood ratio test. We consider two types of genetic maps: dense maps and general maps, and the power of nonparametric functional mapping is investigated through simulation studies and demonstrated by examples.     

A nonlinear mixed-effect mixture model for functional mapping of dynamic traits

Functional mapping has emerged as a next-generation statistical tool for mapping quantitative trait loci (QTL) that affect complex dynamic traits. In this article, we incorporated the idea of nonlinear mixed-effect (NLME) models into the mixture-based framework of functional mapping, aimed to generalize the spectrum of applications for functional mapping. NLME-based functional mapping, implemented with the linearization algorithm based on the first-order Taylor expansion, can provide reasonable estimates of QTL genotypic-specific curve parameters (fixed effect) and the between-individual variation of these parameters (random effect). Results from simulation studies suggest that the NLME-based model is more general than traditional functional mapping. The new model can be useful for the identification of the ontogenetic patterns of QTL genetic effects during time course.     

An algorithm for molecular dissection of tumor progression

The volumetric growth of tumor cells as a function of time is most often likely to be a complex trait, controlled by the combined influences of multiple genes and environmental influences. Genetic mapping has proven to be a powerful tool for detecting and identifying specific genes affecting complex traits, i.e., quantitative trait loci (QTL), based on polymorphic markers. In this article, we present a novel statistical model for genetic mapping of QTL governing tumor growth trajectories in humans. In principle, this model is a combination of functional mapping proposed to map function-valued traits and linkage disequilibrium mapping designed to provide high resolution mapping of QTL by making use of recombination events created at a historic time. We implement an EM-simplex hybrid algorithm for parameter estimation, in which a closed-form solution for the EM algorithm is derived to estimate the population genetic parameters of QTL including the allele frequencies and the coefficient of linkage disequilibrium, and the simplex algorithm incorporated to estimate the curve parameters describing the dynamic changes of cancer cells for different QTL genotypes. Extensive simulations are performed to investigate the statistical properties of our model. Through a number of hypothesis tests, our model allows for cutting-edge studies aimed to decipher the genetic mechanisms underlying cancer growth, development and differentiation. The implications of our model in gene therapy for cancer research are discussed.     

Regression-based multi-trait QTL mapping using a structural equation model

Quantitative trait loci (QTL) mapping often results in data on a number of traits that have well-established causal relationships. Many multi-trait QTL mapping methods that account for the correlation among multiple traits have been developed to improve the statistical power and the precision of QTL parameter estimation. However, none of these methods are capable of incorporating the causal structure among the traits. Consequently, genetic functions of the QTL may not be fully understood. Structural equation modeling (SEM) allows researchers to explicitly characterize the causal structure among the variables and to decompose effects into direct, indirect, and total effects. In this paper, we developed a multi-trait SEM method of QTL mapping that takes into account the causal relationships among traits related to grain yield. Performance of the proposed method is evaluated by simulation study and applied to data from a wheat experiment. Compared with single trait analysis and the multi-trait least-squares analysis, our multi-trait SEM improves statistical power of QTL detection and provides important insight into how QTLs regulate traits by investigating the direct, indirect, and total QTL effects. The approach also helps build biological models that more realistically reflect the complex relationships among QTL and traits and is more precise and efficient in QTL mapping than single trait analysis.     

A dynamic model for functional mapping of biological rhythms

Functional mapping is a statistical method for mapping quantitative trait loci (QTLs) that regulate the dynamic pattern of a biological trait. This method integrates mathematical aspects of biological complexity into a mixture model for genetic mapping and tests the genetic effects of QTLs by comparing genotype-specific curve parameters. As a way of quantitatively specifying the dynamic behaviour of a system, differential equations have proved to be powerful for modelling and unravelling the biochemical, molecular, and cellular mechanisms of a biological process, such as biological rhythms. The equipment of functional mapping with biologically meaningful differential equations provides new insights into the genetic control of any dynamic processes. We formulate a new functional mapping framework for a dynamic biological rhythm by incorporating a group of ordinary differential equations (ODE). The Runge–Kutta fourth-order algorithm was implemented to estimate the parameters that define the system of ODE. The new model will find its implications for understanding the interplay between gene interactions and developmental pathways in complex biological rhythms.     

A hyperspace model to decipher the genetic architecture of developmental processes: allometry meets ontogeny

To better utilize limited resources for their survival and reproduction, all organisms undergo developmental changes in both body size and shape during ontogeny. The genetic analysis of size change with increasing age, i.e., growth, has received considerable attention in quantitative developmental genetic studies, but the genetic architecture of ontogenetic changes in body shape and its associated allometry have been poorly understood partly due to the lack of analytical tools. In this article, we attempt to construct a multivariate statistical framework for studying the genetic regulation of ontogenetic growth and shape. We have integrated biologically meaningful mathematical functions of growth curves and developmental allometry into the estimation process of genetic mapping aimed at identifying individual quantitative trait loci (QTL) for phenotypic variation. This model defined with high dimensions can characterize the ontogenetic patterns of genetic effects of QTL over the lifetime of an organism and assess the interplay between genetic actions/interactions and phenotypic integration. The closed forms for the residual covariance matrix and its determinant and inverse were derived to overcome the computational complexity typical of our high-dimensional model. We used a worked example to validate the utility of this model. The implications of this model for genetic research of evo-devo are discussed.     

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.     

A likelihood approach for mapping growth trajectories using dominant markers in a phase-unknown full-sib family

Dominant markers have been commonly used in mapping quantitative trait loci (QTLs) in outcrossing species, in which not much prior genome information is available. But the dominant nature of these markers may lead to reduced QTL mapping precision and power. A new statistical method is proposed to incorporate growth laws into a QTL mapping framework, under which the use of the efficiency of dominant markers can be increased. This new method can be used to identify specific QTLs affecting differentiation in growth trajectories, and further estimate the timing of a QTL to turn on, or turn off, affecting growth during the entire ontogeny of a species. Using this method based on dominant markers we have successfully mapped a QTL for stem height growth trajectories to a linkage group in a forest tree. The implications of this method for the understanding of the genetic architecture of growth using dominant markers are discussed.Communicated by F. Salamini