Complex genetic effects in quantitative trait locus identification: a computationally tractable random model for use in F(2) populations

Methodology for mapping quantitative trait loci (QTL) has focused primarily on treating the QTL as a fixed effect. These methods differ from the usual models of genetic variation that treat genetic effects as random. Computationally expensive methods that allow QTL to be treated as random have been explicitly developed for additive genetic and dominance effects. By extending these methods with a variance component method (VCM), multiple QTL can be mapped. We focused on an F(2) crossbred population derived from inbred lines and estimated effects for each individual and their corresponding marker-derived genetic covariances. We present extensions to pairwise epistatic effects, which are computationally intensive because a great many individual effects must be estimated. But by replacing individual genetic effects with average genetic effects for each marker class, genetic covariances are approximated. This substantially reduces the computational burden by reducing the dimensions of covariance matrices of genetic effects, resulting in a remarkable gain in the speed of estimating the variance components and evaluating the residual log-likelihood. Preliminary results from simulations indicate competitiveness of the reduced model with multiple-interval mapping, regression interval mapping, and VCM with individual genetic effects in its estimated QTL positions and experimental power.     

Covariance between relatives in multibreed populations: additive model

Covariance between relatives in a multibreed population was derived for an additive model with multiple unlinked loci. An efficient algorithm to compute the inverse of the additive genetic covariance matrix is given. For an additive model, the variance for a crossbred individual is a function of the additive variances for the pure breeds, the covariance between parents, and segregation variances. Provided that the variance of a crossbred individual is computed as presented here, the covariance between crossbred relatives can be computed using formulae for purebred populations. For additive traits the inverse of the genotypic covariance matrix given here can be used both to obtain genetic evaluations by best linear unbiased prediction and to estimate genetic parameters by maximum likelihood in multibreed populations. For nonadditive traits, the procedure currently used to analyze multibreed data can be improved using the theory presented here to compute additive covariances together with a suitable approximation for nonadditive covariances.Supported in part by the Illinois Agricultural Experiment Station, Hatch Projects 35-0345 (RLF) and 35-0367 (MG)     

Reducing bias in the measurement of selection

Selection on quantitative characters is commonly mesured in natural populations using regression techniques based on phenotypic covariances between traits and fitness. However, such methods do not give an accutate view of the causal relationship between the phenotype and fitness if enviornmental factors also contribute to covariances between traits and fitness. A recently developed method for estimating selection eliminates the problem of bias resulting from enviormental covariances. This underappreciated method represents a significant addition to the toolbox of evolutionary ecologist.     

Quantitative trait locus analysis of longitudinal quantitative trait data in complex pedigrees

There is currently considerable interest in genetic analysis of quantitative traits such as blood pressure and body mass index. Despite the fact that these traits change throughout life they are commonly analyzed only at a single time point. The genetic basis of such traits can be better understood by collecting and effectively analyzing longitudinal data. Analyses of these data are complicated by the need to incorporate information from complex pedigree structures and genetic markers. We propose conducting longitudinal quantitative trait locus (QTL) analyses on such data sets by using a flexible random regression estimation technique. The relationship between genetic effects at different ages is efficiently modeled using covariance functions (CFs). Using simulated data we show that the change in genetic effects over time can be well characterized using CFs and that including parameters to model the change in effect with age can provide substantial increases in power to detect QTL compared with repeated measure or univariate techniques. The asymptotic distributions of the methods used are investigated and methods for overcoming the practical difficulties in fitting CFs are discussed. The CF-based techniques should allow efficient multivariate analyses of many data sets in human and natural population genetics.     

Pedigree-free animal models: the relatedness matrix reloaded

Animal models typically require a known genetic pedigree to estimate quantitative genetic parameters. Here we test whether animal models can alternatively be based on estimates of relatedness derived entirely from molecular marker data. Our case study is the morphology of a wild bird population, for which we report estimates of the genetic variance-covariance matrices (G) of six morphological traits using three methods: the traditional animal model; a molecular marker-based approach to estimate heritability based on Ritland's pairwise regression method; and a new approach using a molecular genealogy arranged in a relatedness matrix (R) to replace the pedigree in an animal model. Using the traditional animal model, we found significant genetic variance for all six traits and positive genetic covariance among traits. The pairwise regression method did not return reliable estimates of quantitative genetic parameters in this population, with estimates of genetic variance and covariance typically being very small or negative. In contrast, we found mixed evidence for the use of the pedigree-free animal model. Similar to the pairwise regression method, the pedigree-free approach performed poorly when the full-rank R matrix based on the molecular genealogy was employed. However, performance improved substantially when we reduced the dimensionality of the R matrix in order to maximize the signal to noise ratio. Using reduced-rank R matrices generated estimates of genetic variance that were much closer to those from the traditional model. Nevertheless, this method was less reliable at estimating covariances, which were often estimated to be negative. Taken together, these results suggest that pedigree-free animal models can recover quantitative genetic information, although the signal remains relatively weak. It remains to be determined whether this problem can be overcome by the use of a more powerful battery of molecular markers and improved methods for reconstructing genealogies.     

A phylogenetic model for investigating correlated evolution of substitution rates and continuous phenotypic characters

The comparative approach is routinely used to test for possible correlations between phenotypic or life-history traits. To correct for phylogenetic inertia, the method of independent contrasts assumes that continuous characters evolve along the phylogeny according to a multivariate Brownian process. Brownian diffusion processes have also been used to describe time variations of the parameters of the substitution process, such as the rate of substitution or the ratio of synonymous to nonsynonymous substitutions. Here, we develop a probabilistic framework for testing the coupling between continuous characters and parameters of the molecular substitution process. Rates of substitution and continuous characters are jointly modeled as a multivariate Brownian diffusion process of unknown covariance matrix. The covariance matrix, the divergence times and the phylogenetic variations of substitution rates and continuous characters are all jointly estimated in a Bayesian Monte Carlo framework, imposing on the covariance matrix a prior conjugate to the Brownian process so as to achieve a greater computational efficiency. The coupling between rates and phenotypes is assessed by measuring the posterior probability of positive or negative covariances, whereas divergence dates and phenotypic variations are marginally reconstructed in the context of the joint analysis. As an illustration, we apply the model to a set of 410 mammalian cytochrome b sequences. We observe a negative correlation between the rate of substitution and mass and longevity, which was previously observed. We also find a positive correlation between ω = dN/dS and mass and longevity, which we interpret as an indirect effect of variations of effective population size, thus in partial agreement with the nearly neutral theory. The method can easily be extended to any parameter of the substitution process and to any continuous phenotypic or environmental character.     

The genetic analysis of age-dependent traits: modeling the character process

The extension of classical quantitative genetics to deal with function-valued characters (also called infinite-dimensional characters) such as growth curves, mortality curves, and reaction norms, was begun by Kirkpatrick and co-workers. In this theory, the analogs of variance components for single traits are covariance functions for function-valued traits. In the approach presented here, we employ a variety of parametric models for covariance functions that have a number of desirable properties: the functions (1) are positive definite, (2) can be estimated using procedures like those currently used for single traits, (3) have a small number of parameters, and (4) allow simple hypotheses to be easily tested. The methods are illustrated using data from a large experiment that examined the effects of spontaneous mutations on age-specific mortality rates in Drosophila melanogaster. Our methods are shown to work better than a standard multivariate analysis, which assumes the character value at each age is a distinct character. Advantages over existing methods that model covariance functions as a series of orthogonal polynomials are discussed.     

Accurate genetic and environmental covariance estimation with composite likelihood in genome-wide association studies

Genetic and environmental covariances between pairs of complex traits are important quantitative measurements that characterize their shared genetic and environmental architectures. Accurate estimation of genetic and environmental covariances in genome-wide association studies (GWASs) can help us identify common genetic and environmental factors associated with both traits and facilitate the investigation of their causal relationship. Genetic and environmental covariances are often modeled through multivariate linear mixed models. Existing algorithms for covariance estimation include the traditional restricted maximum likelihood (REML) method and the recent method of moments (MoM). Compared to REML, MoM approaches are computationally efficient and require only GWAS summary statistics. However, MoM approaches can be statistically inefficient, often yielding inaccurate covariance estimates. In addition, existing MoM approaches have so far focused on estimating genetic covariance and have largely ignored environmental covariance estimation. Here we introduce a new computational method, GECKO, for estimating both genetic and environmental covariances, that improves the estimation accuracy of MoM while keeping computation in check. GECKO is based on composite likelihood, relies on only summary statistics for scalable computation, provides accurate genetic and environmental covariance estimates across a range of scenarios, and can accommodate SNP annotation stratified covariance estimation. We illustrate the benefits of GECKO through simulations and applications on analyzing 22 traits from five large-scale GWASs. In the real data applications, GECKO identified 50 significant genetic covariances among analyzed trait pairs, resulting in a twofold power gain compared to the previous MoM method LDSC. In addition, GECKO identified 20 significant environmental covariances. The ability of GECKO to estimate environmental covariance in addition to genetic covariance helps us reveal strong positive correlation between the genetic and environmental covariance estimates across trait pairs, suggesting that common pathways may underlie the shared genetic and environmental architectures between traits.     

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.     

Protein deprivation facilitates the independent evolution of behavior and morphology

Ecological conditions such as nutrition can change genetic covariances between traits and accelerate or slow down trait evolution. As adaptive trait correlations can become maladaptive following rapid environmental change, poor or stressful environments are expected to weaken genetic covariances, thereby increasing the opportunity for independent evolution of traits. Here, we demonstrate the differences in genetic covariance among multiple behavioral and morphological traits (exploration, aggression, and body weight) between southern field crickets (Gryllus bimaculatus) raised in favorable (free‐choice) versus stressful (protein‐deprived) nutritional environments. We also quantify the extent to which differences in genetic covariance structures contribute to the potential for the independent evolution of these traits. We demonstrate that protein‐deprived environments tend to increase the potential for traits to evolve independently, which is caused by genetic covariances that are significantly weaker for crickets raised on protein‐deprived versus free‐choice diets. The weakening effects of stressful environments on genetic covariances tended to be stronger in males than in females. The weakening of the genetic covariance between traits under stressful nutritional environments was expected to facilitate the opportunity for adaptive evolution across generations. Therefore, the multivariate gene‐by‐environment interactions revealed here may facilitate behavioral and morphological adaptations to rapid environmental change.     

A comparative method for both discrete and continuous characters using the threshold model

The threshold model developed by Sewall Wright in 1934 can be used to model the evolution of two-state discrete characters along a phylogeny. The model assumes that there is a quantitative character, called liability, that is unobserved and that determines the discrete character according to whether the liability exceeds a threshold value. A Markov chain Monte Carlo algorithm is used to infer the evolutionary covariances of the liabilities for discrete characters, sampling liability values consistent with the phylogeny and with the observed data. The same approach can also be used for continuous characters by assuming that the tip species have values that have been observed. In this way, one can make a comparative-methods analysis that combines both discrete and continuous characters. Simulations are presented showing that the covariances of the liabilities are successfully estimated, although precision can be achieved only by using a large number of species, and we must always worry whether the covariances and the model apply throughout the group. An advantage of the threshold model is that the model can be straightforwardly extended to accommodate within-species phenotypic variation and allows an interface with quantitative-genetics models.     

Evolutionary quantitative genetics of nonlinear developmental systems

In quantitative genetics, the effects of developmental relationships among traits on microevolution are generally represented by the contribution of pleiotropy to additive genetic covariances. Pleiotropic additive genetic covariances arise only from the average effects of alleles on multiple traits, and therefore the evolutionary importance of nonlinearities in development is generally neglected in quantitative genetic views on evolution. However, nonlinearities in relationships among traits at the level of whole organisms are undeniably important to biology in general, and therefore critical to understanding evolution. I outline a system for characterizing key quantitative parameters in nonlinear developmental systems, which yields expressions for quantities such as trait means and phenotypic and genetic covariance matrices. I then develop a system for quantitative prediction of evolution in nonlinear developmental systems. I apply the system to generating a new hypothesis for why direct stabilizing selection is rarely observed. Other uses will include separation of purely correlative from direct and indirect causal effects in studying mechanisms of selection, generation of predictions of medium‐term evolutionary trajectories rather than immediate predictions of evolutionary change over single generation time‐steps, and the development of efficient and biologically motivated models for separating additive from epistatic genetic variances and covariances.     

Developmental interactions and the constituents of quantitative variation

Development is the process by which genotypes are transformed into phenotypes. Consequently, development determines the relationship between allelic and phenotypic variation in a population and, therefore, the patterns of quantitative genetic variation and covariation of traits. Understanding the developmental basis of quantitative traits may lead to insights into the origin and evolution of quantitative genetic variation, the evolutionary fate of populations, and, more generally, the relationship between development and evolution. Herein, we assume a hierarchical, modular structure of trait development and consider how epigenetic interactions among modules during ontogeny affect patterns of phenotypic and genetic variation. We explore two developmental models, one in which the epigenetic interactions between modules result in additive effects on character expression and a second model in which these epigenetic interactions produce nonadditive effects. Using a phenotype landscape approach, we show how changes in the developmental processes underlying phenotypic expression can alter the magnitude and pattern of quantitative genetic variation. Additive epigenetic effects influence genetic variances and covariances, but allow trait means to evolve independently of the genetic variances and covariances, so that phenotypic evolution can proceed without changing the genetic covariance structure that determines future evolutionary response. Nonadditive epigenetic effects, however, can lead to evolution of genetic variances and covariances as the mean phenotype evolves. Our model suggests that an understanding of multivariate evolution can be considerably enriched by knowledge of the mechanistic basis of character development.     

Natural selection, evolvability and bias due to environmental covariance in the field in an annual plant

Estimates of the form and magnitude of natural selection based on phenotypic relationships between traits and fitness measures can be biased when environmental factors influence both relative fitness and phenotypic trait values. I quantified genetic variances and covariances, and estimated linear and quadratic selection coefficients, for seven traits of an annual plant grown in the field. For replicates of 50 paternal half-sib families, coefficients of selection were calculated both for individual phenotypic values of the traits and for half-sib family mean values. The potential for evolutionary response was supported by significant heritability and phenotypic directional selection for several traits but contradicted by the absence of significant genetic variation for fitness estimates and evidence of bias in phenotypic selection coefficients due to environmental covariance for at least two of the traits analysed. Only studies of a much wider range of organisms and traits will reveal the frequency and extent of such bias.     

A general model of functional constraints on phenotypic evolution

A general model of the functional constraints on the rate and direction of phenotypic evolution is developed using a decomposition of the Lande-Arnold model of multivariate phenotypic evolution. The important feature of the model is the F matrix of performance coefficients reflecting the causal relationship between morphophysiological (m-p) and functional performance traits. The structure of F, which reflects the functional architecture of the organism, constrains the shape of the adaptive landscape and thus the rate and direction of m-p trait evolution. The rate of m-p trait evolution is a function of the pattern of coefficients in a row of F. The sums and variances of these rows are related to current concepts of evolvability. The direction of m-p trait evolution through m-p trait space is a function of the functional covariances among m-p traits. The functional covariance between a pair of m-p traits is a measure of how much the traits function together and is computed as the covariance between rows of F. Finally, it is shown that genetic covariances between m-p traits and performance traits are a function of the F matrix, but a G matrix that includes these covariances cannot be used to model functional constraints effectively.     

RESISTANCE OF GENETIC CORRELATION STRUCTURE TO DIRECTIONAL SELECTION IN DROSOPHILA MELANOGASTER

The genetic covariance and correlation matrices for five morphological traits were estimated from four populations of fruit flies, Drosophila melanogaster, to measure the extent of change in genetic covariances as a result of directional selection. Two of the populations were derived from lines that had undergone selection for large or small thorax length over the preceding 23 generations. A third population was constituted using flies from control lines that were maintained with equivalent population sizes as the selected lines. The fourth population contained flies from the original cage population from which the selected and control lines had been started. Tests of the homogeneity of covariance matrices using maximum likelihood techniques revealed significant changes in covariance structure among the selected lines. Prediction of base population trait means from selected line means under the assumption of constant genetic covariances indicated that genetic covariances for the small population differed more from the base population than did the covariances for the large population. The predicted small population means diverged farther from the expected means because the additive genetic variance associated with several traits increased in value and most of the genetic covariances associated with one trait changed in sign. These results illustrate that genetic covariances may remain nearly constant in some situations while changing markedly in others. Possible developmental reasons for the genetic changes are discussed.     

Analysis of cytoplasmic and maternal effects. II. Genetic models for triploid endosperms

Genetic models for quantitative traits of triploid endosperms are proposed for the analysis of direct gene effects, cytoplasmic effects, and maternal gene effects. The maternal effect is partitioned into maternal additive and dominance components. In the full genetic model, the direct effect is partitioned into direct additive and dominance components and high-order dominance component, which are the cumulative effects of three-allele interactions. If the high-order dominance effects are of no importance, a reduced genetic model can be used. Monte Carlo simulations were conducted in this study for demonstrating unbiasedness of estimated variance and covariance components from the MINQUE (0/1) procedure, which is a minimum norm quadratic unbiased estimation (MINQUE) method setting 0 for all the prior covariances and 1 for all the prior variances. Robustness of estimating variance and covariance components for the genetic models was tested by simulations. Both full and reduced genetic models are shown to be robust for estimating variance and covariance components under several situations of no specific effects. Efficiency of predicting random genetic effects for the genetic models by the MINQUE (0/1) procedure was compared with the best linear unbiased prediction (BLUP). A worked example is given to illustrate the use of the reduced genetic model for kernel growth characteristics in corn (Zea mays L.).     

How to Explore Morphological Integration in Human Evolution and Development?

Most studies in evolutionary developmental biology focus on large-scale evolutionary processes using experimental or molecular approaches, whereas evolutionary quantitative genetics provides mathematical models of the influence of heritable phenotypic variation on the short-term response to natural selection. Studies of morphological integration typically are situated in-between these two styles of explanation. They are based on the consilience of observed phenotypic covariances with qualitative developmental, functional, or evolutionary models. Here we review different forms of integration along with multiple other sources of phenotypic covariances, such as geometric and spatial dependencies among measurements. We discuss one multivariate method [partial least squares analysis (PLS)] to model phenotypic covariances and demonstrate how it can be applied to study developmental integration using two empirical examples. In the first example we use PLS to study integration between the cranial base and the face in human postnatal development. Because the data are longitudinal, we can model both cross-sectional integration and integration of growth itself, i.e., how cross-sectional variance and covariance is actually generated in the course of ontogeny. We find one factor of developmental integration (connecting facial size and the length of the anterior cranial base) that is highly canalized during postnatal development, leading to decreasing cross-sectional variance and covariance. A second factor (overall cranial length to height ratio) is less canalized and leads to increasing (co)variance. In a second example, we examine the evolutionary significance of these patterns by comparing cranial integration in humans to that in chimpanzees.     

Genetic analysis of crossfostering data with sire and dam records

A method is presented for the analysis of data from crossfostering experiments in which parts of litters are reciprocally interchanged at birth. Observed variances and covariances of differently related individuals are expressed as functions of theoretical causal components of phenotypic variance (additive direct, dominance direct, additive maternal, dominance maternal, direct-maternal covariance, and environmental). Causal components are estimated by weighted least squares analysis of this system of equations, including a ridge-regression procedure to examine consequences of correlation between observed components. Ridge regression suggests that dominance direct genetic variance is generally underestimated, but that narrow-sense heritability estimates are reliable.     

Patterns of phenotypic covariation and correlation in modern humans as viewed from morphological integration

Proportionality of phenotypic and genetic distance is of crucial importance to adequately focus on population history and structure, and it depends on the proportionality of genetic and phenotypic covariance. Constancy of phenotypic covariances is unlikely without constancy of genetic covariation if the latter is a substantial component of the former. If phenotypic patterns are found to be relatively stable, the most probable explanation is that genetic covariance matrices are also stable. Factors like morphological integration account for such stability. Morphological integration can be studied by analyzing the relationships among morphological traits. We present here a comparison of phenotypic correlation and covariance structure among worldwide human populations. Correlation and covariance matrices between 47 cranial traits were obtained for 28 populations, and compared with design matrices representing functional and developmental constraints. Among-population differences in patterns of correlation and covariation were tested for association with matrices of genetic distances (obtained after an examination of 10 Alu-insertions) and with Mahalanobis distances (computed after craniometrical traits). All matrix correlations were estimated by means of Mantel tests. Results indicate that correlation and covariance structure in our species is stable, and that among-group correlation/covariance similarity is not related to genetic or phenotypic distance. Conversely, genetic and morphological distance matrices were highly correlated. Correlation and covariation patterns were largely associated with functional and developmental factors, which probably account for the stability of covariance patterns.