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.     

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.     

A multivariate model for ordinal trait analysis

Many economically important characteristics of agricultural crops are measured as ordinal traits. Statistical analysis of the genetic basis of ordinal traits appears to be quite different from regular quantitative traits. The generalized linear model methodology implemented via the Newton-Raphson algorithm offers improved efficiency in the analysis of such data, but does not take full advantage of the extensive theory developed in the linear model arena. Instead, we develop a multivariate model for ordinal trait analysis and implement an EM algorithm for parameter estimation. We also propose a method for calculating the variance-covariance matrix of the estimated parameters. The EM equations turn out to be extremely similar to formulae seen in standard linear model analysis. Computer simulations are performed to validate the EM algorithm. A real data set is analyzed to demonstrate the application of the method. The advantages of the EM algorithm over other methods are addressed. Application of the method to QTL mapping for ordinal traits is demonstrated using a simulated baclcross (BC) population.     

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.     

Joint mapping of quantitative trait Loci for multiple binary characters

Joint mapping for multiple quantitative traits has shed new light on genetic mapping by pinpointing pleiotropic effects and close linkage. Joint mapping also can improve statistical power of QTL detection. However, such a joint mapping procedure has not been available for discrete traits. Most disease resistance traits are measured as one or more discrete characters. These discrete characters are often correlated. Joint mapping for multiple binary disease traits may provide an opportunity to explore pleiotropic effects and increase the statistical power of detecting disease loci. We develop a maximum-likelihood method for mapping multiple binary traits. We postulate a set of multivariate normal disease liabilities, each contributing to the phenotypic variance of one disease trait. The underlying liabilities are linked to the binary phenotypes through some underlying thresholds. The new method actually maps loci for the variation of multivariate normal liabilities. As a result, we are able to take advantage of existing methods of joint mapping for quantitative traits. We treat the multivariate liabilities as missing values so that an expectation-maximization (EM) algorithm can be applied here. We also extend the method to joint mapping for both discrete and continuous traits. Efficiency of the method is demonstrated using simulated data. We also apply the new method to a set of real data and detect several loci responsible for blast resistance in rice.     

An EM algorithm for mapping binary disease loci: application to fibrosarcoma in a four-way cross mouse family

Many diseases show dichotomous phenotypic variation but do not follow a simple Mendelian pattern of inheritance. Variances of these binary diseases are presumably controlled by multiple loci and environmental variants. A least-squares method has been developed for mapping such complex disease loci by treating the binary phenotypes (0 and 1) as if they were continuous. However, the least-squares method is not recommended because of its ad hoc nature. Maximum Likelihood (ML) and Bayesian methods have also been developed for binary disease mapping by incorporating the discrete nature of the phenotypic distribution. In the ML analysis, the likelihood function is usually maximized using some complicated maximization algorithms (e.g. the Newton-Raphson or the simplex algorithm). Under the threshold model of binary disease, we develop an Expectation Maximization (EM) algorithm to solve for the maximum likelihood estimates (MLEs). The new EM algorithm is developed by treating both the unobserved genotype and the disease liability as missing values. As a result, the EM iteration equations have the same form as the normal equation system in linear regression. The EM algorithm is further modified to take into account sexual dimorphism in the linkage maps. Applying the EM-implemented ML method to a four-way-cross mouse family, we detected two regions on the fourth chromosome that have evidence of QTLs controlling the segregation of fibrosarcoma, a form of connective tissue cancer. The two QTLs explain 50-60% of the variance in the disease liability. We also applied a Bayesian method previously developed (modified to take into account sex-specific maps) to this data set and detected one additional QTL on chromosome 13 that explains another 26% of the variance of the disease liability. All the QTLs detected primarily show dominance effects.     

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.     

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.     

A genetical genomics approach to genome scans increases power for QTL mapping

We describe a method for integrating gene expression information into genome scans and show that this can substantially increase the statistical power of QTL mapping. The method has three stages. First, standard clustering methods identify small (size 5-20) groups of genes with similar expression patterns. Second, each gene group is tested for a causative genetic locus shared with the clinical trait of interest. This is done using an EM algorithm approach that treats genotype at the putative causative locus as an unobserved variable and combines expression information from all of the genes in the group to infer genotype information at the locus. Finally, expression QTL (eQTL) are mapped for each gene group that shares a causative locus with the clinical trait. Such eQTL are candidates for the causative locus. Simulation results show that this method has far superior power to standard QTL mapping techniques in many circumstances. We applied this method to existing data on mouse obesity. Our method identified 27 putative body weight QTL, whereas standard QTL mapping produced only one. Furthermore, most gene groups with body weight QTL included cis genes, so candidate genes could be immediately identified. Eleven body weight QTL produced 16 candidate genes that have been previously associated with body weight or body weight-related traits, thus validating our method. In addition, 15 of the 16 other loci produced 32 candidate genes that have not been associated with body weight. Thus, this method shows great promise for finding new causative loci for complex traits.     

Progress and promise of genome-wide association studies for human complex trait genetics

Enormous progress in mapping complex traits in humans has been made in the last 5 yr. There has been early success for prevalent diseases with complex phenotypes. These studies have demonstrated clearly that, while complex traits differ in their underlying genetic architectures, for many common disorders the predominant pattern is that of many loci, individually with small effects on phenotype. For some traits, loci of large effect have been identified. For almost all complex traits studied in humans, the sum of the identified genetic effects comprises only a portion, generally less than half, of the estimated trait heritability. A variety of hypotheses have been proposed to explain why this might be the case, including untested rare variants, and gene-gene and gene-environment interaction. Effort is currently being directed toward implementation of novel analytic approaches and testing rare variants for association with complex traits using imputed variants from the publicly available 1000 Genomes Project resequencing data and from direct resequencing of clinical samples. Through integration with annotations and functional genomic data as well as by in vitro and in vivo experimentation, mapping studies continue to characterize functional variants associated with complex traits and address fundamental issues such as epistasis and pleiotropy. This review focuses primarily on the ways in which genome-wide association studies (GWASs) have revolutionized the field of human quantitative genetics.     

Statistical optimization of parametric accelerated failure time model for mapping survival trait loci

Most existing statistical methods for mapping quantitative trait loci (QTL) are not suitable for analyzing survival traits with a skewed distribution and censoring mechanism. As a result, researchers incorporate parametric and semi-parametric models of survival analysis into the framework of the interval mapping for QTL controlling survival traits. In survival analysis, accelerated failure time (AFT) model is considered as a de facto standard and fundamental model for data analysis. Based on AFT model, we propose a parametric approach for mapping survival traits using the EM algorithm to obtain the maximum likelihood estimates of the parameters. Also, with Bayesian information criterion (BIC) as a model selection criterion, an optimal mapping model is constructed by choosing specific error distributions with maximum likelihood and parsimonious parameters. Two real datasets were analyzed by our proposed method for illustration. The results show that among the five commonly used survival distributions, Weibull distribution is the optimal survival function for mapping of heading time in rice, while Log-logistic distribution is the optimal one for hyperoxic acute lung injury.     

远交群体动态性状基因定位的似然分析Ⅰ.理论方法

受动物遗传育种中用来估计动态性状育种值的随机回归测定日模型思想的启发 ,将关于时间 (测定日期 )的Legendre多项式镶嵌在遗传模型的每个遗传效应中 ,以刻画QTL对动态性状变化过程的作用 ,从而建立起动态性状基因定位的数学模型。利用远交设计群体 ,阐述了动态性状基因定位的似然分析原理 ,推导了定位参数似然估计的EM法两步求解过程。结合动态性状遗传分析的特点和普通数量性状基因定位研究进展 ,还提出了有关动态性状基因定位进一步研究的设想     

Generalized Linear Model for Mapping Discrete Trait Loci Implemented with LASSO Algorithm

Generalized estimating equation (GEE) algorithm under a heterogeneous residual variance model is an extension of the iteratively reweighted least squares (IRLS) method for continuous traits to discrete traits. In contrast to mixture model-based expectation–maximization (EM) algorithm, the GEE algorithm can well detect quantitative trait locus (QTL), especially large effect QTLs located in large marker intervals in the manner of high computing speed. Based on a single QTL model, however, the GEE algorithm has very limited statistical power to detect multiple QTLs because of ignoring other linked QTLs. In this study, the fast least absolute shrinkage and selection operator (LASSO) is derived for generalized linear model (GLM) with all possible link functions. Under a heterogeneous residual variance model, the LASSO for GLM is used to iteratively estimate the non-zero genetic effects of those loci over entire genome. The iteratively reweighted LASSO is therefore extended to mapping QTL for discrete traits, such as ordinal, binary, and Poisson traits. The simulated and real data analyses are conducted to demonstrate the efficiency of the proposed method to simultaneously identify multiple QTLs for binary and Poisson traits as examples.     

On the generalized poisson regression mixture model for mapping quantitative trait loci with count data

Statistical methods for mapping quantitative trait loci (QTL) have been extensively studied. While most existing methods assume normal distribution of the phenotype, the normality assumption could be easily violated when phenotypes are measured in counts. One natural choice to deal with count traits is to apply the classical Poisson regression model. However, conditional on covariates, the Poisson assumption of mean-variance equality may not be valid when data are potentially under- or overdispersed. In this article, we propose an interval-mapping approach for phenotypes measured in counts. We model the effects of QTL through a generalized Poisson regression model and develop efficient likelihood-based inference procedures. This approach, implemented with the EM algorithm, allows for a genomewide scan for the existence of QTL throughout the entire genome. The performance of the proposed method is evaluated through extensive simulation studies along with comparisons with existing approaches such as the Poisson regression and the generalized estimating equation approach. An application to a rice tiller number data set is given. Our approach provides a standard procedure for mapping QTL involved in the genetic control of complex traits measured in counts.     

A unified DNA sequence and non‐DNA sequence mapping model of complex traits

Increasing evidence shows that quantitative inheritance is based on both DNA sequence and non‐DNA sequence variants. However, how to simultaneously detect these variants from a mapping study has been unexplored, hampering our effort to illustrate the detailed genetic architecture of complex traits. We address this issue by developing a unified model of quantitative trait locus (QTL) mapping based on an open‐pollinated design composed of randomly sampling maternal plants from a natural population and their half‐sib seeds. This design forms a two‐level hierarchical platform for a joint linkage‐linkage disequilibrium analysis of population structure. The EM algorithm was implemented to estimate and test DNA sequence‐based effects and non‐DNA sequence‐based effects of QTLs. We applied this model to analyze genetic mapping data from the OP design of a gymnosperm coniferous species, Torreya grandis, identifying 25 significant DNA sequence and non‐DNA sequence QTLs for seedling height and diameter growth in different years. Results from computer simulation show that the unified model has good statistical properties and is powerful for QTL detection. Our model enables the tests of how a complex trait is affected differently by DNA‐based effects and non‐DNA sequence‐based transgenerational effects, thus allowing a more comprehensive picture of genetic architecture to be charted and quantified.     

The full EM algorithm for the MLEs of QTL effects and positions and their estimated variances in multiple-interval mapping

The advent of complete genetic linkage maps of DNA markers has made systematic studies of mapping quantitative trait loci (QTL) in experimental organisms feasible. The method of multiple-interval mapping provides an appropriate way for mapping QTL using genetic markers. However, efficient algorithms for the computation involved remain to be developed. In this article, a full EM algorithm for the simultaneous computation of the MLEs of QTL effects and positions is developed. EM-based formulas are derived for computing the observed Fisher information matrix. The full EM algorithm is compared with an ECM algorithm developed by Kao and Zeng (1997, Biometrics 53, 653-665). The validity of the inverted observed Fisher information matrix as an estimate of the variance matrix of the MLEs is demonstrated by a simulation study.     

Characterization of genomic imprinting effects and patterns with parametric accelerated failure time model

Genomic imprinting, a genetic phenomenon of non-equivalent allele expression that depends on parental origins, has been ubiquitously observed in nature. It does not only control the traits of growth and development but also may be responsible for survival traits. Based on the accelerated failure time model, we construct a general parametric model for mapping the imprinted QTL (iQTL). Within the framework of interval mapping, maximum likelihood estimation of iQTL parameters is implemented via EM algorithm. The imprinting patterns of the detected iQTL are statistically tested according to a series of null hypotheses. BIC model selection criterion is employed to choose an optimal baseline hazard function with maximum likelihood and parsimonious parameters. Simulations are used to validate the proposed mapping procedure. A published dataset from a mouse model system was used to illustrate the proposed framework. Results show that among the five commonly used survival distributions, Log-logistic distribution is the optimal baseline hazard function for mapping QTL of hyperoxic acute lung injury (HALI) survival; under the log-logistic distribution, four QTLs were identified, in which only one QTL was inherited in Mendelian fashion, whereas others were imprinted in different imprinting patterns.     

A model for quantitative trait loci mapping,linkage phase,and segregation pattern estimation for a full-sib progeny

Quantitative trait loci (QTL) mapping is an important approach for the study of the genetic architecture of quantitative traits. For perennial species, inbred lines cannot be obtained due to inbreed depression and a long juvenile period. Instead, linkage mapping can be performed by using a full-sib progeny. This creates a complex scenario because both markers and QTL alleles can have different segregation patterns as well as different linkage phases between them. We present a two-step method for QTL mapping using full-sib progeny based on composite interval mapping (i.e., interval mapping with cofactors), considering an integrated genetic map with markers with different segregation patterns and conditional probabilities obtained by a multipoint approach. The model is based on three orthogonal contrasts to estimate the additive effect (one in each parent) and dominance effect. These estimatives are obtained using the EM algorithm. In the first step, the genome is scanned to detect QTL. After, segregation pattern and linkage phases between QTL and markers are estimated. A simulated example is presented to validate the methodology. In general, the new model is more effective than existing approaches, because it can reveal QTL present in a full-sib progeny that segregates in any pattern present and can also identify dominance effects. Also, the inclusion of cofactors provided more statistical power for QTL mapping.     

Mapping quantitative trait loci for longitudinal traits in line crosses

Quantitative traits whose phenotypic values change over time are called longitudinal traits. Genetic analyses of longitudinal traits can be conducted using any of the following approaches: (1) treating the phenotypic values at different time points as repeated measurements of the same trait and analyzing the trait under the repeated measurements framework, (2) treating the phenotypes measured from different time points as different traits and analyzing the traits jointly on the basis of the theory of multivariate analysis, and (3) fitting a growth curve to the phenotypic values across time points and analyzing the fitted parameters of the growth trajectory under the theory of multivariate analysis. The third approach has been used in QTL mapping for longitudinal traits by fitting the data to a logistic growth trajectory. This approach applies only to the particular S-shaped growth process. In practice, a longitudinal trait may show a trajectory of any shape. We demonstrate that one can describe a longitudinal trait with orthogonal polynomials, which are sufficiently general for fitting any shaped curve. We develop a mixed-model methodology for QTL mapping of longitudinal traits and a maximum-likelihood method for parameter estimation and statistical tests. The expectation-maximization (EM) algorithm is applied to search for the maximum-likelihood estimates of parameters. The method is verified with simulated data and demonstrated with experimental data from a pseudobackcross family of Populus (poplar) trees.     

Evolutionary consequences of many-to-one mapping of jaw morphology to mechanics in labrid fishes

Many physiological traits consist of two hierarchically related levels: physical structures and the emergent functional properties of those structures. Because selection tends to act on the emergent functional traits, the evolution of structural phenotypes will depend on the nature of the form-function relationship. Complex physiological or biomechanical traits are often characterized by many-to-one mapping: numerous structural phenotypes can yield equivalent functions. We suggest that this redundancy can promote the evolution of phenotypic diversity, and we illustrate this effect with a combination of empirical and analytical studies of a complex biomechanical trait, the four-bar linkage found in the jaws of labrid fishes. We show that labrid jaws are subject to many-to-one mapping of form-to-jaw mechanical properties but that some mechanical types have higher levels of morphological redundancy than others. This variation in redundancy has affected the diversity and distribution of labrid jaw shapes: labrid species are disproportionately concentrated around functional traits with higher potential for redundancy. Many-to-one mapping can also mitigate evolutionary constraints imposed by mechanical trade-offs by allowing a species to simultaneously optimize multiple functional properties. Many-to-one mapping may be an important factor in generating the uneven patterns of diversity in physiological traits.