Regression based fast multi-trait genome-wide QTL analysis

Multivariate simple interval mapping (SIM) is one of the most popular approaches for multiple quantitative trait locus (QTL) analysis. Both maximum likelihood (ML) and least squares (LS) multivariate regression (MVR) are widely used methods for multi-trait SIM. ML-based MVR (MVR-ML) is an expectation maximization (EM) algorithm based iterative and complex time-consuming approach. Although the LS-based MVR (MVR-LS) approach is not an iterative process, the calculation of likelihood ratio (LR) statistic in MVR-LS is also a time-consuming complex process. We have introduced a new approach (called FastMtQTL) for multi-trait QTL analysis based on the assumption of multivariate normal distribution of phenotypic observations. Our proposed method can identify almost the same QTL positions as those identified by the existing methods. Moreover, the proposed method takes comparatively less computation time because of the simplicity in the calculation of LR statistic by this method. In the proposed method, LR statistic is calculated only using the sample variance–covariance matrix of phenotypes and the conditional probability of QTL genotype given the marker genotypes. This improvement in computation time is advantageous when the numbers of phenotypes and individuals are larger, and the markers are very dense resulting in a QTL mapping with a bigger dataset.

Statistical models for genotype by environment data: from conventional ANOVA models to eco-physiological QTL models

To study the performance of genotypes under different growing conditions, plant breeders evaluate their germplasm in multi-environment trials. These trials produce genotype × environment data. We present statistical models for the analysis of such data that differ in the extent to which additional genetic, physiological, and environmental information is incorporated into the model formulation. The simplest model in our exposition is the additive 2-way analysis of variance model, without genotype × environment interaction, and with parameters whose interpretation depends strongly on the set of included genotypes and environments. The most complicated model is a synthesis of a multiple quantitative trait locus (QTL) model and an eco-physiological model to describe a collection of genotypic response curves. Between those extremes, we discuss linear-bilinear models, whose parameters can only indirectly be related to genetic and physiological information, and factorial regression models that allow direct incorporation of explicit genetic, physiological, and environmental covariables on the levels of the genotypic and environmental factors. Factorial regression models are also very suitable for the modelling of QTL main effects and QTL × environment interaction. Our conclusion is that statistical and physiological models can be fruitfully combined for the study of genotype × environment interaction.

Multiple Quantitative Trait Locus Analysis of Bovine Chromosome 6 in the Israeli Holstein Population by a Daughter Design

Nine Israeli Holstein sire families with 2978 daughters were analyzed for quantitative trait loci effects on chromosome 6 for five milk production traits by a daughter design. All animals were genotyped for 2 markers. The three families with significant effects were genotyped for up to 10 additional markers spanning positions 0–122 cM of BTA6. Two sires were segregating for a locus affecting protein and fat percentage near position 55 cM with an estimated substitution effect of 0.18% protein, which is equivalent to one phenotypic standard deviation. This locus was localized to a confidence interval of 4 cM. One of these sires was also heterozygous for a locus affecting milk, fat, and protein production near the centromere. The hypothesis of two segregating loci was verified by multiple regression analysis. A third sire was heterozygous for a locus affecting milk and protein percentage near the telomeric end of the chromosome. Possible candidates for the major quantitative gene near position 55 cM were determined by comparative mapping. IBSP and SSP1 were used as anchors for the orthologous region on human chromosome 4. Twelve genes were detected within a 2-Mbp sequence. None of these genes have been previously associated with lactogenesis.