BACKGROUND
large amounts of genomic information generated from single nucleotide polymorphism  microarrays have become available in recent years for many species
<cit> .
this genomic information is used to detect polymorphisms that contribute to variation in economically important traits, such as production traits in farm animals
<cit> .
microarray technology is also used to screen the expression levels of thousands of genes, i.e., the transcriptome
<cit> .
studies have shown that genetic background can have a large impact on differential expression
<cit> .
integrating genome and transcriptome information can help to elucidate the underlying biology of the genotype-phenotype map, using expression quantitative trait locus  mapping
<cit> .
however, in the eqtl approach, associations between snps, transcript level, and phenotypes are analyzed individually.
this is likely to lead to “missing heritability”
<cit> , because corrections for multiple testing lead to a high false negative rate and multiple snps and transcript level that jointly explain the phenotype are ignored
<cit> .
here we propose and demonstrate bayesian models that model all snps and transcript level simultaneously to obtain explained variances by the whole genome and whole transcriptome.
in these models, we identify eqtls as those snps whose effects disappear when transcript level are added to the model.
genomic- and transcriptomic-explained variances are further partitioned by chromosome and genome sections to offer a view of the genetic architecture on different aggregation levels.
the choice of bayesian variable selection  models was due to its features to separate markers with large/moderate or small effects, and to locate the important regions in the genome or transcriptome which serves a better qtl mapping method because it produces clearer signals for qtl
<cit> .
furthermore the prediction based on genomic variables using bvs is more accurate even when the prior is not correct
<cit> .
it is important to say that simpler methods suffer from “missing heritability” too
<cit> .
the aim of this study was to explore the contributions of various sources of variation, such as population structure, snp variants, and gene expression levels, to a set of growth related traits  in mice.
these traits are very important, both in terms of agricultural production and for obesity in humans.
bayesian mixed models and bayesian variable selection models were applied to model pedigree, snps and/or gene expressions and to derive explained variances for these components.
in addition, they were used to partition of snps and gene expression by chromosome and genome sections.
to validate the estimates of explained variances, the predictive ability of these models was studied using cross validation.
data
an m <dig> × icr f <dig> population of  <dig> mice was available with complete records for body weight at 8 weeks  and  <dig> records for feed intake  and feed efficiency , measured during the period 3 weeks to 8 weeks
<cit> .
an additional  <dig> pedigree records were available that described the family structure up to the f <dig> founder lines.
data was obtained in three batches and the sex of the animals was recorded.
at the end of the experiment, the mice were sacrificed and liver tissue was extracted for genome-wide expression profiling.
rna isolation, cdna synthesis, array hybridization, normalization of probe level intensity, and annotation of data were performed as described in detail by
<cit> .
genotypes for  <dig> highly informative single nucleotide polymorphisms  were available for each animal.
these tag-snps were used to trace the genomic variation in this f <dig> population.
density functions of phenotypes are available in additional file
<dig> and the whole data were made publicly available at .
methods
the most complete model used describes phenotypes y  by an intercept μ, environmental effects of batch and sex b, a polygenic effect based on pedigree u, regressions on snp covariates a, regressions on gene expression covariates g, and a model residual e, as:
y=1µ+xb+zu+wa+qg+e
where x is the design matrix for batch and sex effects, z is a design matrix that links polygenic effects to the observed records, w is a matrix with  <dig> snp covariates, and q is a matrix with  <dig>  gene expression covariates.
the snp and gene expression covariates were centered and scaled to unit variance.
based on work of
<cit> , the bayesian mixed model version assigns normal prior to the vectors u, a, g, and e in , i.e.,
u∼n <dig> aσu <dig> a∼n <dig> iσs <dig> g∼n <dig> iσg <dig> e∼n <dig> iσs <dig>  where
σe <dig> is the polygenic variance and a is the numerator relationship matrix based on pedigree information,
σs <dig> is the per-snp explained variance,
σg <dig> is the per-gene expression explained variance, and
σe <dig> is the residual or environmental variance.
these four variances are estimated in the model using flat prior distributions, i.e.,
σu <dig> σs <dig> σg <dig> σe2∼bern.
the remaining parameters in , μ and b, are assigned flat prior distributions, which is the bayesian analog of fitting “fixed effects”  estimates.
a markov chain monte carlo  algorithm was applied in the software bayz
<cit>  to obtain samples from the posterior distribution of the model parameters
fμ,b,u,a,g,σu <dig> σs <dig> σg <dig> σe2|y.
mcmc algorithms for sampling effects and variances in mixed models have been extensively described, for a general overview see
<cit> .
the monte carlo accuracy of the mcmc algorithm was evaluated by correlating repeated estimates for the parameter vectors u, a and g, requiring a correlation > <dig>  from repeated mcmc runs, and by computing the effective sample sizes for the variance components using the r coda package
<cit> .
the explained variance in y from  is var + var + var + var.
to obtain posterior means  and posterior standard deviations  on the explained variances for snps and gene expressions, var and var were evaluated based on the posterior samples for a and g from the mcmc, i.e., as the pm and psd of var values over mcmc cycles, where at is the posterior sample for a from mcmc cycle t. this procedure is not required for the polygenic variance, because z is a design matrix, unlike w and q, which are covariate matrices.
the second model used was a bayesian variable selection model, where the approach of george and mcculloch
<cit>  was followed to fit mixture distributions with small and large variances as the prior distribution for regression coefficients.
in model , such a mixture prior was applied to snps as well as gene expression regression coefficients, with independent parameters and mixture indicators for snps and for gene expressions.
the basic model of george and mcculloch
<cit>  was further extended to incorporate the variances in the mixture distribution as unknown model parameters, which allows the model to learn the relative importance of snps and gene expressions from the data.
this variable selection model thus takes the prior distributions for a and g as follows:
ai∼γain <dig> τa12+1−γain <dig> τa <dig>
gi∼γgin <dig> τg12+1−γgin <dig> τg <dig>
where
τa <dig> and
τa <dig> are the “large” and “small” variances in the mixture distribution for a,
τg <dig> and
τg <dig> are the “large” and “small” variances in the mixture distribution for g, and
γa and
γg are vectors of 0/ <dig> indicator variables for a and g, respectively, indicating whether the ith element in a or g, respectively, comes from the distribution with large or small variance.
the variances
τa <dig> τa <dig> τg <dig> τg <dig> were all estimated from the data using unbounded flat prior distributions.
the constraints
τa12>τa <dig> and
τg12>τg <dig> were applied using a rejection sampler, so that “large” and “small” effects remained identifiable.
the priors for the indicator variables were taken as
γai∼bernπa and
γgi∼bernπg, where
bernπmeans a bernoulli distribution for a 0/ <dig> indicator with a probability π for a  <dig>  the parameters
πa,πgwere taken as known.
the mcmc implementation of this model is relatively straightforward, because conditional on the indicator variables the model remains a mixed model.
the updating of the mixture indicators is described in
<cit> .
this model is also run in the software bayz
<cit> , and the monte carlo accuracy was evaluated in the same way as the mixed model version.
from the posterior samples for a and g in the variable selection model, explained variances were computed and partitioned by chromosome and by genome section.
the variable selection model is more suited to make such a partitioning, because unlike the mixed model version, it allows for different variance contributions per snp.
the explained variances were evaluated in the same way as for the mixed model, by evaluating var and var over mcmc cycles t, except that the a and g samples are obtained under the mixture model prior assumptions.
the same expressions can be straightforwardly evaluated for parts of the snps or gene expressions to obtain explained variances per chromosome and for small windows of snps within chromosomes.
variance within a chromosome was computed using a 5-snp sliding window to obtain a genomic variance profile.
it is difficult to choose an optimal windows size as it depends on extend of ld, marker density and an arbitrary cut-off for what is considered important ld.
in the data analyzed here, average r <dig> between adjacent snps was  <dig> , and average r <dig> between snps two apart was  <dig> , which we considered sufficiently high to warrant computation of variances in a 5-snp window.
to study the relative importance of family structure, snps, and gene expressions, six sub models and the complete model  were used.
these were models that use only pedigree information , only snp data , only gene expression data , snp + gex, ped + gex, ped + snp, and the complete model ped + snp + gex.
these models always included sex and batch effects.
the predictive ability of the models was evaluated using an 11-fold cross-validation.
for body weight,  <dig> records were divided randomly in  <dig> groups, each with  <dig> individuals.
feed intake and feed efficiency, with  <dig> records in total, were randomly divided in  <dig> groups of  <dig> records and one group of  <dig> records.
the complete model, including all variance parameters, was re-estimated on each set of  <dig> folds and predictions were computed for the phenotypes in the remaining 11th fold.
all predictions from the 11-fold cross validation were collected to compute correlations between predicted and actual phenotypes, and regressions of predicted phenotypes on actual phenotypes, using the whole data set.
the slope of the regression lines of predicted phenotypes on actual phenotypes are expected to be  <dig> if the model produces unbiased predictions, which would validate the estimates of explained variances.
the university of nebraska institutional animal care and use committee approved all procedures and protocols.
RESULTS
table
<dig> presents estimates of explained variances for the three traits using the seven models considered.
the results in table
<dig> were obtained using the bayesian mixed model.
we first discuss the models that consider genetic and genomic information, which are the ped, snp and ped + snp models.
the ped model is the classical polygenic model, using family structure to estimate narrow sense heritability, which yielded estimates of 42%, 53%, and 58% for bw, fi, and fe, respectively.
genomic information alone  explained less variance, i.e., 36%, 28%, and 24% for bw, fi, and fe respectively.
it is a common finding that snps explain less variance than the classical heritability estimates
<cit> , which is attributed to causal variants having lower minor allele frequency than the genotyped snps
<cit> , insufficient modeling of identity by descent by snps
<cit> , and incomplete linkage disequilibrium  between causal variants and genotyped snps
<cit> .
combining pedigree and snp data  increased the explained variance above that of using pedigree only, i.e., for bw the ped + snp model obtained an explained variance of 59%, compared to 42% for the ped model.
this phenomenon is particularly common in the analysis of an f <dig> population, where increased genetic variance in the f <dig> can be captured by snps, but not by pedigree.
in the ped + snp model, the part covered by pedigree decreased compared to the ped only model, showing that snps cover part of the family relationships
<cit> .
explained variances are for residuals , polygenic effects , snps , and gene expressions .
the table shows estimates as the posterior mean with posterior standard deviation in parentheses and the proportion of explained variance as percentage of the total.
overall, explained variances increase by adding gene expression information , i.e., in the most complete model  explained variances were 88%, 75%, and 71% for bw, fi, and fe respectively.
this confirms the assumption that gene expressions can explain a larger part of phenotypic variance than genetic or genomic information, by capturing environmental, and possibly non-additive, genetic effects through the gene expressions
<cit> .
information on the genetic architecture of these traits is best judged from the relative contributions of genomic and transcriptomic data in the snp + gex model.
this model shows that, for these traits, the liver transcriptome contributes a larger portion of explained variance.
this is most pronounced for bw, with 18% of explained variance from the genome and 63% from the liver transcriptome.
thus, in this case, the predominant model is that snps regulate gene expressions to exert their effect on the phenotype.
figure
<dig> shows a decomposition of the explained variances at the chromosome and sub-chromosome level for the models using genomic  and genomic with transcriptomic  data for the trait bw.
for the traits fi and fe, see additional file
<dig> and additional file
<dig> respectively.
these results are based on the bayesian variable selection model to better differentiate between genomic regions contributing more and less variance.
the genomic variances at the sub-chromosome level are explained variances in a sliding 5-snp window.
at the chromosome level, chromosome  <dig> particularly stands out, with a relatively large contribution from the snps effects via transcriptome, but only a small contribution from the genome alone in explaining the phenotype.
this does not mean there is no important qtl on this chromosome.
in fact, there is a large qtl on chromosome 10; however, it is an eqtl whose effect can be captured by gene expressions.
figures
1b and c show the details at the sub-chromosome level, with figure
1b showing the explained genomic variances when fitting snps alone , and figure
1c showing the explained genomic variances when adding gene expressions to the model .
the differences between these two graphs show locations of qtls that regulate the liver transcripts and qtls that exert their effect on the phenotype through another route.
for bw, the main qtls are found on chromosomes  <dig>   <dig>   <dig>   <dig>  and  <dig>  and the qtls on  <dig>   <dig>  and  <dig> appear to be eqtls affecting gene expression in the liver.
the qtl on chromosome  <dig> is an intermediate case whose effect is reduced, but does not completely disappear, when adding gene expressions to the model.
thus, this chromosome  <dig> qtl regulates liver transcripts, but must also have effects on bw through other routes, possibly by regulating genes outside the liver.
the chromosome  <dig> qtl is a clear case of a qtl whose effect on bw does not work via the regulation of liver transcripts.
the qtl locations are in agreement with qtls detected for body weight in other studies
<cit> .
the same graphs for traits fi and fe are provided as supplementary material.
these traits show relatively more cases where qtl effects remain after adding liver transcriptome data, which is in agreement with results in table
<dig>
this method/approach is suitable for gene-level resolution.
however, gene-level resolution is highly data dependent, i.e.
it requires high marker density and a study population with ld blocks that span small genomic regions.
in this work we have used f <dig> crosses from outbred lines, which has large ld blocks and this kind of data has limited resolution for fine-mapping of qtl.
one may argue that the most complete model is more interesting to investigate genetic architecture and chromosomal/sub-chromosomal variance but as we have shown snps and pedigree are largely confounded and they explain about the same variance.
this confounded explained variance is getting worse in the case that both pedigree and snps are in one model  which is shown in higher confidence intervals of explained variance by pedigree.
the model with only omics information  is therefore simpler, more accurate and as effective as the model that also uses pedigree information.
this is interesting for future applications of omics technologies, because we expect that pedigree information often will be absent.
figures
<dig> and
<dig> present detailed graphs of the genomic variances  and the distribution of chromosomal genomic values or breeding values
<cit>  of the animals  for chromosomes  <dig> and  <dig>  and for models fitting snp only  or snp + gex .
breeding value is defined as the value of an individual as a parent based on sum of its genes effects
<cit> .
chromosome  <dig> is the case of an apparent eqtl, showing that genomic variance disappears, and that a tri-modal distribution of genomic values collapses, when gene expressions are added to the model.
chromosome  <dig> is the case of a qtl that does not regulate liver transcripts.
the detailed picture of chromosome  <dig> shows that adding gene expressions to the model makes the effects of this qtl clearer: genomic variances outside the qtl region reduce, and a clear tri-modal distribution of chromosomal genomic values is seen in the snp + gex model, but not in the snp-only model.
table
<dig> shows the rank correlations between the predicted values from using pedigree , genomic , or transcriptomic  information.
pedigree and genomic values correlate better than pedigree/genomic values with transcriptomic values.
this confirms that pedigree and genomic information overlap to a reasonable degree, but this is less true for transcriptomic information.
the prediction of phenotypes from these models, using cross-validation, is shown in table
<dig>  showing correlation between predicted and actual phenotypes, and the regression of predicted phenotype on actual phenotype.
the scatter plot of predicted versus actual phenotypes is shown in additional file
<dig>  the results for explained variance and for prediction do not necessarily coincide, because prediction is also affected by the accuracy of the parameter estimates.
the results show that predictions from the snp model are all as good, or better, than from the ped model, while the explained variances from the snp model were lower .
this can be explained by the snp predictions being more accurate than ped predictions.
models including gene expressions show the highest correlations with phenotypes, meaning that models including gene expressions also provide accurate predictions.
the regressions of predicted phenotype on actual phenotype are mostly around  <dig>  indicating that the predictions are unbiased and that the explained variances where correctly assessed.
ρ correlation between true phenotype and predicted value.
β, regression of predicted values on true phenotypes.
CONCLUSIONS
with increased availability of various -omics data, integrative approaches are promising tools for understanding the genetic architecture of complex traits.
we have developed a complementary approach to the univariate “eqtl” mapping, by considering bayesian models that fit all genome-wide snps and transcript abundances in one model, and that estimate and partition explained variances by chromosome and genome segments.
our results show that, using gene expressions, more of the phenotypic variance can be explained and phenotypes can be better predicted.
predictions were also shown to be unbiased, which validates the assessed explained variances.
the improvement of phenotype predictions using gene expression data will be useful for several applications in agriculture and medicine, although it should be assessed on a case-by-case basis as to whether a suitable tissue can be sampled for the gene expression measurements.
partitioning of the explained genomic variance at the level of chromosomes and genome segments showed clear examples of eqtl locations as regions where genomic variance disappears when gene expressions are added to the model.
our study used only gene expressions from the liver, and an obvious further extension is to include expressions from other tissues.
the qtls that did not disappear when transcripts are added to the model may be eqtls that affect gene expression in a tissue other than liver.
the bayesian model is quite efficient for handling large sets of covariates, and extensions to include multiple sets of expressions will be feasible.
we have not provided formal statistical tests in this model, but the bayesian approach lends itself naturally to obtaining confidence intervals for  parameter estimates.
the estimates of total explained variances from the bayesian mixed model can also be obtained by a residual maximum likelihood  approach.
we verified this, and the bayesian and reml estimates generally agree.
however, using reml it is not feasible to utilize mixture priors to better discriminate between snps which contribute more or less variance, and to partition the variances at the sub-chromosome level, which is all straightforward in a bayesian approach.
our approach can easily allow up scaling to higher-density arrays, even to whole-genome sequence data with the variance components analysis as it was for gene expression probes in this study.
abbreviations
bw: body weight; fi: feed intake; fe: feed efficiency; snps: single nucleotide polymorphisms; reml: restricted maximum likelihood; qtl: quantitative trait loci; eqtl: expression quantitative trait loci.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
ae developed the data analysis pipeline, performed statistical analyses, interpreted the results and wrote the manuscript.
ps and lj were involved in project design, statistical analyses, interpretation of results and manuscript editing.
dp and ma prepared the data for the analysis.
all authors have read and approved the final manuscript.
supplementary material
additional file 1
figure s <dig>  distribution of phenotypes of traits body weight including  <dig> animals, feed intake and feed efficiency including  <dig> animals each.
click here for file
additional file 2
figure s <dig>  decomposition of the proportion of variance explained by snps at the level of chromosomes and individual snps in two models: the independent model snp and the conditional model snp+gex for feed intake.
explained variances from snps in snp model  and snp+gex model  in each chromosome.
explained variance by individual snps in snp model and  snp+gex model.
click here for file
additional file 3
figure s <dig>  decomposition of the proportion of variance explained by snps at the level of chromosomes and individual snps in two models: the independent model snp and the conditional model snp+gex for feed efficiency.
explained variances from snps in snp model  and snp+gex model  in each chromosome.
explained variance by individual snps in snp model and  snp+gex model.
click here for file
additional file 4
figure s <dig>  comparison of predicted breeding values versus phenotypes in the models using pedigree information only , snps information only  and gene expression information only  for three traits body weight, feed intake and feed efficiency according to correlation shown in table
<dig>
click here for file
acknowledgement
this research is supported in part by the quantomics research project that has been co-financed by the european commission within the 7th framework programme, contract no.
<dig>  this work is a part of phd project scholarship from the ministry of science, research and technology of iran.

 large amounts of genomic information generated from single nucleotide polymorphism microarrays have become available in recent years for many species. this genomic information is used to detect polymorphisms that contribute to variation in economically important traits, such as production traits in farm animals. microarray technology is also used to screen the expression levels of thousands of genes, i.e., the transcriptome. studies have shown that genetic background can have a large impact on differential expression. integrating genome and transcriptome information can help to elucidate the underlying biology of the genotype-phenotype map, using expression quantitative trait locus mapping. however, in the eqtl approach, associations between snps, transcript level, and phenotypes are analyzed individually. this is likely to lead to “missing heritability”, because corrections for multiple testing lead to a high false negative rate and multiple snps and transcript level that jointly explain the phenotype are ignored. here we propose and demonstrate bayesian models that model all snps and transcript level simultaneously to obtain explained variances by the whole genome and whole transcriptome. in these models, we identify eqtls as those snps whose effects disappear when transcript level are added to the model. genomic- and transcriptomic-explained variances are further partitioned by chromosome and genome sections to offer a view of the genetic architecture on different aggregation levels. the choice of bayesian variable selection models was due to its features to separate markers with large/moderate or small effects, and to locate the important regions in the genome or transcriptome which serves a better qtl mapping method because it produces clearer signals for qtl. furthermore the prediction based on genomic variables using bvs is more accurate even when the prior is not correct. it is important to say that simpler methods suffer from “missing heritability” too. the aim of this study was to explore the contributions of various sources of variation, such as population structure, snp variants, and gene expression levels, to a set of growth related traits in mice. these traits are very important, both in terms of agricultural production and for obesity in humans. bayesian mixed models and bayesian variable selection models were applied to model pedigree, snps and/or gene expressions and to derive explained variances for these components. in addition, they were used to partition of snps and gene expression by chromosome and genome sections. to validate the estimates of explained variances, the predictive ability of these models was studied using cross validation. data an m × icr f population of mice was available with complete records for body weight at 8 weeks and records for feed intake and feed efficiency, measured during the period 3 weeks to 8 weeks. an additional pedigree records were available that described the family structure up to the f founder lines. data was obtained in three batches and the sex of the animals was recorded. at the end of the experiment, the mice were sacrificed and liver tissue was extracted for genome-wide expression profiling. rna isolation, cdna synthesis, array hybridization, normalization of probe level intensity, and annotation of data were performed as described in detail by. genotypes for highly informative single nucleotide polymorphisms were available for each animal. these tag-snps were used to trace the genomic variation in this f population. density functions of phenotypes are available in additional file and the whole data were made publicly available at. methods the most complete model used describes phenotypes y by an intercept μ, environmental effects of batch and sex b, a polygenic effect based on pedigree u, regressions on snp covariates a, regressions on gene expression covariates g, and a model residual e, as: y=1µ+xb+zu+wa+qg+e where x is the design matrix for batch and sex effects, z is a design matrix that links polygenic effects to the observed records, w is a matrix with snp covariates, and q is a matrix with gene expression covariates. the snp and gene expression covariates were centered and scaled to unit variance. based on work of, the bayesian mixed model version assigns normal prior to the vectors u, a, g, and e in, i.e., u∼n aσu a∼n iσs g∼n iσg e∼n iσs where σe is the polygenic variance and a is the numerator relationship matrix based on pedigree information, σs is the per-snp explained variance, σg is the per-gene expression explained variance, and σe is the residual or environmental variance. these four variances are estimated in the model using flat prior distributions, i.e., σu σs σg σe2∼bern. the remaining parameters in, μ and b, are assigned flat prior distributions, which is the bayesian analog of fitting “fixed effects” estimates. a markov chain monte carlo algorithm was applied in the software bayz to obtain samples from the posterior distribution of the model parameters fμ,b,u,a,g,σu σs σg σe2|y. mcmc algorithms for sampling effects and variances in mixed models have been extensively described, for a general overview see. the monte carlo accuracy of the mcmc algorithm was evaluated by correlating repeated estimates for the parameter vectors u, a and g, requiring a correlation > from repeated mcmc runs, and by computing the effective sample sizes for the variance components using the r coda package. the explained variance in y from is var + var + var + var. to obtain posterior means and posterior standard deviations on the explained variances for snps and gene expressions, var and var were evaluated based on the posterior samples for a and g from the mcmc, i.e., as the pm and psd of var values over mcmc cycles, where at is the posterior sample for a from mcmc cycle t. this procedure is not required for the polygenic variance, because z is a design matrix, unlike w and q, which are covariate matrices. the second model used was a bayesian variable selection model, where the approach of george and mcculloch was followed to fit mixture distributions with small and large variances as the prior distribution for regression coefficients. in model, such a mixture prior was applied to snps as well as gene expression regression coefficients, with independent parameters and mixture indicators for snps and for gene expressions. the basic model of george and mcculloch was further extended to incorporate the variances in the mixture distribution as unknown model parameters, which allows the model to learn the relative importance of snps and gene expressions from the data. this variable selection model thus takes the prior distributions for a and g as follows: ai∼γain τa12+1−γain τa gi∼γgin τg12+1−γgin τg where τa and τa are the “large” and “small” variances in the mixture distribution for a, τg and τg are the “large” and “small” variances in the mixture distribution for g, and γa and γg are vectors of 0/ indicator variables for a and g, respectively, indicating whether the ith element in a or g, respectively, comes from the distribution with large or small variance. the variances τa τa τg τg were all estimated from the data using unbounded flat prior distributions. the constraints τa12>τa and τg12>τg were applied using a rejection sampler, so that “large” and “small” effects remained identifiable. the priors for the indicator variables were taken as γai∼bernπa and γgi∼bernπg, where bernπmeans a bernoulli distribution for a 0/ indicator with a probability π for a the parameters πa,πgwere taken as known. the mcmc implementation of this model is relatively straightforward, because conditional on the indicator variables the model remains a mixed model. the updating of the mixture indicators is described in. this model is also run in the software bayz, and the monte carlo accuracy was evaluated in the same way as the mixed model version. from the posterior samples for a and g in the variable selection model, explained variances were computed and partitioned by chromosome and by genome section. the variable selection model is more suited to make such a partitioning, because unlike the mixed model version, it allows for different variance contributions per snp. the explained variances were evaluated in the same way as for the mixed model, by evaluating var and var over mcmc cycles t, except that the a and g samples are obtained under the mixture model prior assumptions. the same expressions can be straightforwardly evaluated for parts of the snps or gene expressions to obtain explained variances per chromosome and for small windows of snps within chromosomes. variance within a chromosome was computed using a 5-snp sliding window to obtain a genomic variance profile. it is difficult to choose an optimal windows size as it depends on extend of ld, marker density and an arbitrary cut-off for what is considered important ld. in the data analyzed here, average r between adjacent snps was, and average r between snps two apart was, which we considered sufficiently high to warrant computation of variances in a 5-snp window. to study the relative importance of family structure, snps, and gene expressions, six sub models and the complete model were used. these were models that use only pedigree information, only snp data, only gene expression data, snp + gex, ped + gex, ped + snp, and the complete model ped + snp + gex. these models always included sex and batch effects. the predictive ability of the models was evaluated using an 11-fold cross-validation. for body weight, records were divided randomly in groups, each with individuals. feed intake and feed efficiency, with records in total, were randomly divided in groups of records and one group of records. the complete model, including all variance parameters, was re-estimated on each set of folds and predictions were computed for the phenotypes in the remaining 11th fold. all predictions from the 11-fold cross validation were collected to compute correlations between predicted and actual phenotypes, and regressions of predicted phenotypes on actual phenotypes, using the whole data set. the slope of the regression lines of predicted phenotypes on actual phenotypes are expected to be if the model produces unbiased predictions, which would validate the estimates of explained variances. the university of nebraska institutional animal care and use committee approved all procedures and protocols.  table  presents estimates of explained variances for the three traits using the seven models considered. the results in table  were obtained using the bayesian mixed model. we first discuss the models that consider genetic and genomic information, which are the ped, snp and ped + snp models. the ped model is the classical polygenic model, using family structure to estimate narrow sense heritability, which yielded estimates of 42%, 53%, and 58% for bw, fi, and fe, respectively. genomic information alone explained less variance, i.e., 36%, 28%, and 24% for bw, fi, and fe respectively. it is a common finding that snps explain less variance than the classical heritability estimates, which is attributed to causal variants having lower minor allele frequency than the genotyped snps, insufficient modeling of identity by descent by snps, and incomplete linkage disequilibrium between causal variants and genotyped snps. combining pedigree and snp data increased the explained variance above that of using pedigree only, i.e., for bw the ped + snp model obtained an explained variance of 59%, compared to 42% for the ped model. this phenomenon is particularly common in the analysis of an f population, where increased genetic variance in the f can be captured by snps, but not by pedigree. in the ped + snp model, the part covered by pedigree decreased compared to the ped only model, showing that snps cover part of the family relationships. explained variances are for residuals, polygenic effects, snps, and gene expressions. the table shows estimates as the posterior mean with posterior standard deviation in parentheses and the proportion of explained variance as percentage of the total. overall, explained variances increase by adding gene expression information, i.e., in the most complete model explained variances were 88%, 75%, and 71% for bw, fi, and fe respectively. this confirms the assumption that gene expressions can explain a larger part of phenotypic variance than genetic or genomic information, by capturing environmental, and possibly non-additive, genetic effects through the gene expressions. information on the genetic architecture of these traits is best judged from the relative contributions of genomic and transcriptomic data in the snp + gex model. this model shows that, for these traits, the liver transcriptome contributes a larger portion of explained variance. this is most pronounced for bw, with 18% of explained variance from the genome and 63% from the liver transcriptome. thus, in this case, the predominant model is that snps regulate gene expressions to exert their effect on the phenotype. figure shows a decomposition of the explained variances at the chromosome and sub-chromosome level for the models using genomic and genomic with transcriptomic data for the trait bw. for the traits fi and fe, see additional file and additional file respectively. these results are based on the bayesian variable selection model to better differentiate between genomic regions contributing more and less variance. the genomic variances at the sub-chromosome level are explained variances in a sliding 5-snp window. at the chromosome level, chromosome particularly stands out, with a relatively large contribution from the snps effects via transcriptome, but only a small contribution from the genome alone in explaining the phenotype. this does not mean there is no important qtl on this chromosome. in fact, there is a large qtl on chromosome 10; however, it is an eqtl whose effect can be captured by gene expressions. figures  1b and c show the details at the sub-chromosome level, with figure 1b showing the explained genomic variances when fitting snps alone, and figure 1c showing the explained genomic variances when adding gene expressions to the model. the differences between these two graphs show locations of qtls that regulate the liver transcripts and qtls that exert their effect on the phenotype through another route. for bw, the main qtls are found on chromosomes and and the qtls on and appear to be eqtls affecting gene expression in the liver. the qtl on chromosome is an intermediate case whose effect is reduced, but does not completely disappear, when adding gene expressions to the model. thus, this chromosome qtl regulates liver transcripts, but must also have effects on bw through other routes, possibly by regulating genes outside the liver. the chromosome qtl is a clear case of a qtl whose effect on bw does not work via the regulation of liver transcripts. the qtl locations are in agreement with qtls detected for body weight in other studies. the same graphs for traits fi and fe are provided as supplementary material. these traits show relatively more cases where qtl effects remain after adding liver transcriptome data, which is in agreement with results in table  this method/approach is suitable for gene-level resolution. however, gene-level resolution is highly data dependent, i.e. it requires high marker density and a study population with ld blocks that span small genomic regions. in this work we have used f crosses from outbred lines, which has large ld blocks and this kind of data has limited resolution for fine-mapping of qtl. one may argue that the most complete model is more interesting to investigate genetic architecture and chromosomal/sub-chromosomal variance but as we have shown snps and pedigree are largely confounded and they explain about the same variance. this confounded explained variance is getting worse in the case that both pedigree and snps are in one model which is shown in higher confidence intervals of explained variance by pedigree. the model with only omics information is therefore simpler, more accurate and as effective as the model that also uses pedigree information. this is interesting for future applications of omics technologies, because we expect that pedigree information often will be absent. figures  and present detailed graphs of the genomic variances and the distribution of chromosomal genomic values or breeding values of the animals for chromosomes and and for models fitting snp only or snp + gex. breeding value is defined as the value of an individual as a parent based on sum of its genes effects. chromosome is the case of an apparent eqtl, showing that genomic variance disappears, and that a tri-modal distribution of genomic values collapses, when gene expressions are added to the model. chromosome is the case of a qtl that does not regulate liver transcripts. the detailed picture of chromosome shows that adding gene expressions to the model makes the effects of this qtl clearer: genomic variances outside the qtl region reduce, and a clear tri-modal distribution of chromosomal genomic values is seen in the snp + gex model, but not in the snp-only model. table  shows the rank correlations between the predicted values from using pedigree, genomic, or transcriptomic information. pedigree and genomic values correlate better than pedigree/genomic values with transcriptomic values. this confirms that pedigree and genomic information overlap to a reasonable degree, but this is less true for transcriptomic information. the prediction of phenotypes from these models, using cross-validation, is shown in table  showing correlation between predicted and actual phenotypes, and the regression of predicted phenotype on actual phenotype. the scatter plot of predicted versus actual phenotypes is shown in additional file the results for explained variance and for prediction do not necessarily coincide, because prediction is also affected by the accuracy of the parameter estimates. the results show that predictions from the snp model are all as good, or better, than from the ped model, while the explained variances from the snp model were lower. this can be explained by the snp predictions being more accurate than ped predictions. models including gene expressions show the highest correlations with phenotypes, meaning that models including gene expressions also provide accurate predictions. the regressions of predicted phenotype on actual phenotype are mostly around indicating that the predictions are unbiased and that the explained variances where correctly assessed. ρ correlation between true phenotype and predicted value. β, regression of predicted values on true phenotypes.  with increased availability of various -omics data, integrative approaches are promising tools for understanding the genetic architecture of complex traits. we have developed a complementary approach to the univariate “eqtl” mapping, by considering bayesian models that fit all genome-wide snps and transcript abundances in one model, and that estimate and partition explained variances by chromosome and genome segments. our results show that, using gene expressions, more of the phenotypic variance can be explained and phenotypes can be better predicted. predictions were also shown to be unbiased, which validates the assessed explained variances. the improvement of phenotype predictions using gene expression data will be useful for several applications in agriculture and medicine, although it should be assessed on a case-by-case basis as to whether a suitable tissue can be sampled for the gene expression measurements. partitioning of the explained genomic variance at the level of chromosomes and genome segments showed clear examples of eqtl locations as regions where genomic variance disappears when gene expressions are added to the model. our study used only gene expressions from the liver, and an obvious further extension is to include expressions from other tissues. the qtls that did not disappear when transcripts are added to the model may be eqtls that affect gene expression in a tissue other than liver. the bayesian model is quite efficient for handling large sets of covariates, and extensions to include multiple sets of expressions will be feasible. we have not provided formal statistical tests in this model, but the bayesian approach lends itself naturally to obtaining confidence intervals for parameter estimates. the estimates of total explained variances from the bayesian mixed model can also be obtained by a residual maximum likelihood approach. we verified this, and the bayesian and reml estimates generally agree. however, using reml it is not feasible to utilize mixture priors to better discriminate between snps which contribute more or less variance, and to partition the variances at the sub-chromosome level, which is all straightforward in a bayesian approach. our approach can easily allow up scaling to higher-density arrays, even to whole-genome sequence data with the variance components analysis as it was for gene expression probes in this study. abbreviations bw: body weight; fi: feed intake; fe: feed efficiency; snps: single nucleotide polymorphisms; reml: restricted maximum likelihood; qtl: quantitative trait loci; eqtl: expression quantitative trait loci. competing interests the authors declare that they have no competing interests. authors’ contributions ae developed the data analysis pipeline, performed statistical analyses, interpreted the results and wrote the manuscript. ps and lj were involved in project design, statistical analyses, interpretation of results and manuscript editing. dp and ma prepared the data for the analysis. all authors have read and approved the final manuscript. supplementary material additional file 1 figure s distribution of phenotypes of traits body weight including animals, feed intake and feed efficiency including animals each. click here for file additional file 2 figure s decomposition of the proportion of variance explained by snps at the level of chromosomes and individual snps in two models: the independent model snp and the conditional model snp+gex for feed intake. explained variances from snps in snp model and snp+gex model in each chromosome. explained variance by individual snps in snp model and snp+gex model. click here for file additional file 3 figure s decomposition of the proportion of variance explained by snps at the level of chromosomes and individual snps in two models: the independent model snp and the conditional model snp+gex for feed efficiency. explained variances from snps in snp model and snp+gex model in each chromosome. explained variance by individual snps in snp model and snp+gex model. click here for file additional file 4 figure s comparison of predicted breeding values versus phenotypes in the models using pedigree information only, snps information only and gene expression information only for three traits body weight, feed intake and feed efficiency according to correlation shown in table click here for file acknowledgement this research is supported in part by the quantomics research project that has been co-financed by the european commission within the 7th framework programme, contract no. this work is a part of phd project scholarship from the ministry of science, research and technology of iran.
for bw, the main quantitative trait loci are found on chromosomes and and the qtls on and appear to be expression quantitative trait locus affecting gene expression in the liver. chromosome is the case of an apparent eqtl, showing that genomic variance disappears, and that a tri-modal distribution of genomic values collapses, when gene expressions are added to the model. with increased availability of various -omics data, integrative approaches are promising tools for understanding the genetic architecture of complex traits. we analyzed three growth-related traits: body weight, feed intake, and feed efficiency, in an f population of mice. to understand the genetic architecture of complex traits and bridge the genotype-phenotype gap, it is useful to study intermediate -omics data, e.g. the present study introduces a method for simultaneous quantification of the contributions from single nucleotide polymorphisms and transcript abundances in explaining phenotypic variance, using bayesian whole-omics models. the models that include transcripts explained more phenotypic variance and were better at predicting phenotypes than a model using snps alone. bayesian mixed models and variable selection models were used and, based on parameter samples from the model posterior distributions, explained variances were further partitioned at the level of chromosomes and genome segments. the genomic variation was covered by tag snps, and transcript abundances were available from probes measured in the liver. comparison of these models showed that for bw, a large part of the variation explained by snps could be covered by the liver transcript abundances; this was less true for fi and fe. partitioning of explained variances at the chromosome and genome-segment level clearly separated regulatory and structural genomic variation as the areas where snp effects disappeared/remained after adding transcripts to the model. the predictions from these bayesian models are generally unbiased, validating the estimates of explained variances.
 to understand the genetic architecture of complex traits and bridge the genotype-phenotype gap, it is useful to study intermediate -omics data, e.g. the transcriptome. the present study introduces a method for simultaneous quantification of the contributions from single nucleotide polymorphisms and transcript abundances in explaining phenotypic variance, using bayesian whole-omics models. bayesian mixed models and variable selection models were used and, based on parameter samples from the model posterior distributions, explained variances were further partitioned at the level of chromosomes and genome segments.  we analyzed three growth-related traits: body weight, feed intake, and feed efficiency, in an f population of mice. the genomic variation was covered by tag snps, and transcript abundances were available from probes measured in the liver. explained variances were computed for models using pedigree, snps, transcripts, and combinations of these. comparison of these models showed that for bw, a large part of the variation explained by snps could be covered by the liver transcript abundances; this was less true for fi and fe. for bw, the main quantitative trait loci are found on chromosomes and and the qtls on and appear to be expression quantitative trait locus affecting gene expression in the liver. chromosome is the case of an apparent eqtl, showing that genomic variance disappears, and that a tri-modal distribution of genomic values collapses, when gene expressions are added to the model.  with increased availability of various -omics data, integrative approaches are promising tools for understanding the genetic architecture of complex traits. partitioning of explained variances at the chromosome and genome-segment level clearly separated regulatory and structural genomic variation as the areas where snp effects disappeared/remained after adding transcripts to the model. the models that include transcripts explained more phenotypic variance and were better at predicting phenotypes than a model using snps alone. the predictions from these bayesian models are generally unbiased, validating the estimates of explained variances. bayesianbody weightfeed intakegenometranscriptomeeqtlvariance 