Comment on "On the Metropolis-Hastings acceptance probability to add or drop a quantitative trait locus in Markov chain Monte Carlo-based Bayesian analyses"

Using the realized relationship matrix to disentangle confounding factors for the estimation of genetic variance components of complex traits

BACKGROUND

In the analysis of complex traits, genetic effects can be confounded with non-genetic effects, especially when using full-sib families. Dominance and epistatic effects are typically confounded with additive genetic and non-genetic effects. This confounding may cause the estimated genetic variance components to be inaccurate and biased.

METHODS

In this study, we constructed genetic covariance structures from whole-genome marker data, and thus used realized relationship matrices to estimate variance components in a heterogenous population of approximately 2200 mice for which four complex traits were investigated. These mice were genotyped for more than 10,000 single nucleotide polymorphisms (SNP) and the variances due to family, cage and genetic effects were estimated by models based on pedigree information only, aggregate SNP information, and model selection for specific SNP effects.

RESULTS AND CONCLUSIONS

We show that the use of genome-wide SNP information can disentangle confounding factors to estimate genetic variances by separating genetic and non-genetic effects. The estimated variance components using realized relationship were more accurate and less biased, compared to those based on pedigree information only. Models that allow the selection of individual SNP in addition to fitting a relationship matrix are more efficient for traits with a significant dominance variance.

Bayesian QTL mapping for multiple families derived from crossing a set of inbred lines to a reference line

In some crop species, germplasm collections consisting of a large number of accessions that include traditional landraces, modern cultivars and wild species have recently been established. Such collections are regarded as useful stocks of genes for breeding programs. However, to efficiently utilize these collections for plant breeding, understanding genetic variation in agronomic traits at the QTL level between the accessions is indispensable. One effective way to extract the actual QTL information included in these collections is to perform QTL analysis jointly for multiple families derived from crossing some accessions of the collection with a single reference line such as a standard commercial variety. We developed a Bayesian method for jointly analyzing QTL in such interconnected multiple families derived from a set of inbred lines crossed to the reference line, to detect QTL segregating between any of the inbred lines and the reference line. In this study, we considered multiple recombinant inbred lines, each of which was derived from crossing each of the inbred lines to the reference line. The method was evaluated through the use of simulated data sets for its efficiency in detecting QTL and identifying families segregating at each QTL.

Predicting unobserved phenotypes for complex traits from whole-genome SNP data

Genome-wide association studies (GWAS) for quantitative traits and disease in humans and other species have shown that there are many loci that contribute to the observed resemblance between relatives. GWAS to date have mostly focussed on discovery of genes or regulatory regions habouring causative polymorphisms, using single SNP analyses and setting stringent type-I error rates. Genome-wide marker data can also be used to predict genetic values and therefore predict phenotypes. Here, we propose a Bayesian method that utilises all marker data simultaneously to predict phenotypes. We apply the method to three traits: coat colour, %CD8 cells, and mean cell haemoglobin, measured in a heterogeneous stock mouse population. We find that a model that contains both additive and dominance effects, estimated from genome-wide marker data, is successful in predicting unobserved phenotypes and is significantly better than a prediction based upon the phenotypes of close relatives. Correlations between predicted and actual phenotypes were in the range of 0.4 to 0.9 when half of the number of families was used to estimate effects and the other half for prediction. Posterior probabilities of SNPs being associated with coat colour were high for regions that are known to contain loci for this trait. The prediction of phenotypes using large samples, high-density SNP data, and appropriate statistical methodology is feasible and can be applied in human medicine, forensics, or artificial selection programs.

Results from recent genome-wide association studies indicate that for most complex traits, there are many loci that contribute to variation in observed phenotype and that the effect of a single variant (single nucleotide polymorphism, SNP) on a phenotype is small. Here, we propose a method that combines the effects of multiple SNPs to make a prediction of a phenotype that has not been observed. We apply the method to data on mice, using phenotypic and genomic data from some individuals to predict phenotypes in other, either related or unrelated, individuals. We find that correlations between predicted and actual phenotypes are in the range of 0.4 to 0.9. The method also shows that the SNPs used in the prediction appear in regions that are known to contain genes associated with the traits studied. The prediction of unobserved phenotypes from high-density SNP data and appropriate statistical methodology is feasible and can be applied in human medicine, forensics, or artificial breeding programs.

A Bayesian method for simultaneously detecting Mendelian and imprinted quantitative trait loci in experimental crosses of outbred species

Genomic imprinting is interpreted as a phenomenon, in which some genes inherited from one parent are not completely expressed due to modification of the genome caused during gametogenesis. Subsequently, the expression level of an allele at the imprinted gene is changed dependent on the parental origin, which is referred to as the parent-of-origin effect. In livestock, some QTL for reproductive performance and meat productivity have been reported to be imprinted. So far, methods detecting imprinted QTL have been proposed on the basis of interval mapping, where only a single QTL was tested at a time. In this study, we developed a Bayesian method for simultaneously mapping multiple QTL, allowing the inference about expression modes of QTL in an outbred F2 family. The inference about whether a QTL is Mendelian or imprinted was made using Markov chain Monte Carlo estimation by comparing the goodness-of-fits between models, assuming the presence and the absence of parent-of-origin effect at a QTL. We showed by the analyses of simulated data sets that the Bayesian method can effectively detect both Mendelian QTL and imprinted QTL.

Bayesian mapping of quantitative trait loci for multiple complex traits with the use of variance components

Complex traits important for humans are often correlated phenotypically and genetically. Joint mapping of quantitative-trait loci (QTLs) for multiple correlated traits plays an important role in unraveling the genetic architecture of complex traits. Compared with single-trait analysis, joint mapping addresses more questions and has advantages for power of QTL detection and precision of parameter estimation. Some statistical methods have been developed to map QTLs underlying multiple traits, most of which are based on maximum-likelihood methods. We develop here a multivariate version of the Bayes methodology for joint mapping of QTLs, using the Markov chain-Monte Carlo (MCMC) algorithm. We adopt a variance-components method to model complex traits in outbred populations (e.g., humans). The method is robust, can deal with an arbitrary number of alleles with arbitrary patterns of gene actions (such as additive and dominant), and allows for multiple phenotype data of various types in the joint analysis (e.g., multiple continuous traits and mixtures of continuous traits and discrete traits). Under a Bayesian framework, parameters--including the number of QTLs--are estimated on the basis of their marginal posterior samples, which are generated through two samplers, the Gibbs sampler and the reversible-jump MCMC. In addition, we calculate the Bayes factor related to each identified QTL, to test coincident linkage versus pleiotropy. The performance of our method is evaluated in simulations with full-sib families. The results show that our proposed Bayesian joint-mapping method performs well for mapping multiple QTLs in situations of either bivariate continuous traits or mixed data types. Compared with the analysis for each trait separately, Bayesian joint mapping improves statistical power, provides stronger evidence of QTL detection, and increases precision in estimation of parameter and QTL position. We also applied the proposed method to a set of real data and detected a coincident linkage responsible for determining bone mineral density and areal bone size of wrist in humans.

Simultaneous fine mapping of multiple closely linked quantitative trait Loci using combined linkage disequilibrium and linkage with a general pedigree

Within a small region (e.g., <10 cM), there can be multiple quantitative trait loci (QTL) underlying phenotypes of a trait. Simultaneous fine mapping of closely linked QTL needs an efficient tool to remove confounded shade effects among QTL within such a small region. We propose a variance component method using combined linkage disequilibrium (LD) and linkage information and a reversible jump Markov chain Monte Carlo (MCMC) sampling for model selection. QTL identity-by-descent (IBD) coefficients between individuals are estimated by a hybrid MCMC combining the random walk and the meiosis Gibbs sampler. These coefficients are used in a mixed linear model and an empirical Bayesian procedure combines residual maximum likelihood (REML) to estimate QTL effects and a reversible jump MCMC that samples the number of QTL and the posterior QTL intensities across the tested region. Note that two MCMC processes are used, i.e., an (internal) MCMC for IBD estimation and an (external) MCMC for model selection. In a simulation study, the use of the multiple-QTL model clearly removes the shade effects between three closely linked QTL located at 1.125, 3.875, and 7.875 cM across the region of 10 cM, using 40 markers at 0.25-cM intervals. It is shown that the use of combined LD and linkage information gives much more useful information compared to using linkage information alone for both single- and multiple-QTL analyses. When using a lower marker density (11 markers at 1-cM intervals), the signal of the second QTL can disappear. Extreme values of past effective size (resulting in extreme levels of LD) decrease the mapping accuracy.

Bayesian mapping of QTL in outbred F2 families allowing inference about whether F0 grandparents are homozygous or heterozygous at QTL

In this paper, we propose a new Bayesian method for QTL analysis in outbred F2 families based on Markov chain Monte Carlo (MCMC) estimation allowing inference about whether each of F0 founders (grandparents) is homozygous or heterozygous at QTL. This, in turn, allows us to select a model accurately explaining observations of phenotypes for F2 individuals. The proposed method performs the fitting a statistical model of the two possible QTL states in each F0 grandparent, that is, homozygous and heterozygous at QTL, and gives a posterior distribution for the QTL states in each F0 grandparent. We confine ourselves to the discrimination of two QTL states, homozygous or heterozygous, for each of the F0 grandparents without taking into consideration whether common alleles are shared by F0 grandparents. The statistical model includes allelic effects and dominance effects for each QTL. The number of parameters representing allelic effects and dominance effects is therefore changed depending on the QTL states. A Reversible Jump MCMC technique is used for transition between the models of different dimensions. The effectiveness of the proposed method was investigated using simulation experiments. It was practicable to estimate the QTL states of F0 grandparents as well as the number, the locations and the effects of QTL segregating in an outbred F2 family.