Persistence of accuracy of genome-wide breeding values over generations when including a polygenic effect

Background

When estimating marker effects in genomic selection, estimates of marker effects may simply act as a proxy for pedigree, i.e. their effect may partially be attributed to their association with superior parents and not be linked to any causative QTL. Hence, these markers mainly explain polygenic effects rather than QTL effects. However, if a polygenic effect is included in a Bayesian model, it is expected that the estimated effect of these markers will be more persistent over generations without having to re-estimate the marker effects every generation and will result in increased accuracy and reduced bias.

Methods

Genomic selection using the Bayesian method, ''BayesB'' was evaluated for different marker densities when a polygenic effect is included (GWpEBV) and not included (GWEBV) in the model. Linkage disequilibrium and a mutation drift balance were obtained by simulating a population with a Ne of 100 over 1,000 generations.

Results

Accuracy of selection was slightly higher for the model including a polygenic effect than for the model not including a polygenic effect whatever the marker density. The accuracy decreased in later generations, and this reduction was stronger for lower marker densities. However, no significant difference in accuracy was observed between the two models. The linear regression of TBV on GWEBV and GWpEBV was used as a measure of bias. The regression coefficient was more stable over generations when a polygenic effect was included in the model, and was always between 0.98 and 1.00 for the highest marker density. The regression coefficient decreased more quickly with decreasing marker density.

Conclusions

Including a polygenic effect had no impact on the selection accuracy, but showed reduced bias, which is especially important when estimates of genome-wide markers are used to estimate breeding values over more than one generation.     

Sex ratio at birth in India, its relation to birth order, sex of previous children and use of indigenous medicine

Objective

Sex-ratio at birth in families with previous girls is worse than those with a boy. Our aim was to prospectively study in a large maternal and child unit sex-ratio against previous birth sex and use of traditional medicines for sex selection.

Main Outcome Measures

Sex-ratio among mothers in families with a previous girl and in those with a previous boy, prevalence of indigenous medicine use and sex-ratio in those using medicines for sex selection.

Results

Overall there were 806 girls to 1000 boys. The sex-ratio was 720∶1000 if there was one previous girl and 178∶1000 if there were two previous girls. In second children of families with a previous boy 1017 girls were born per 1000 boys. Sex-ratio in those with one previous girl, who were taking traditional medicines for sex selection, was 928∶1000.

Conclusion

Evidence from the second children clearly shows the sex-ratio is being manipulated by human interventions. More mothers with previous girls tend to use traditional medicines for sex selection, in their subsequent pregnancies. Those taking such medication do not seem to be helped according to expectations. They seem to rely on this method and so are less likely use more definitive methods like sex selective abortions. This is the first such prospective investigation of sex ratio in second children looked at against the sex of previous children. More studies are needed to confirm the findings.     

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0068-2) contains supplementary material, which is available to authorized users.     

Selection in action: dissecting the molecular underpinnings of the increasing muscle mass of Belgian Blue Cattle

Background

Belgian Blue cattle are famous for their exceptional muscular development or “double-muscling”. This defining feature emerged following the fixation of a loss-of-function variant in the myostatin gene in the eighties. Since then, sustained selection has further increased muscle mass of Belgian Blue animals to a comparable extent. In the present paper, we study the genetic determinants of this second wave of muscle growth.

Results

A scan for selective sweeps did not reveal the recent fixation of another allele with major effect on muscularity. However, a genome-wide association study identified two genome-wide significant and three suggestive quantitative trait loci (QTL) affecting specific muscle groups and jointly explaining 8-21% of the heritability. The top two QTL are caused by presumably recent mutations on unique haplotypes that have rapidly risen in frequency in the population. While one appears on its way to fixation, the ascent of the other is compromised as the likely underlying MRC2 mutation causes crooked tail syndrome in homozygotes. Genomic prediction models indicate that the residual additive variance is largely polygenic.

Conclusions

Contrary to complex traits in humans which have a near-exclusive polygenic architecture, muscle mass in beef cattle (as other production traits under directional selection), appears to be controlled by (i) a handful of recent mutations with large effect that rapidly sweep through the population, and (ii) a large number of presumably older variants with very small effects that rise slowly in the population (polygenic adaptation).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-796) contains supplementary material, which is available to authorized users.     

Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar)

Background

Flesh colour and growth related traits in salmonids are both commercially important and of great interest from a physiological and evolutionary perspective. The aim of this study was to identify quantitative trait loci (QTL) affecting flesh colour and growth related traits in an F2 population derived from an isolated, landlocked wild population in Norway (Byglands Bleke) and a commercial production population.

Methods

One hundred and twenty-eight informative microsatellite loci distributed across all 29 linkage groups in Atlantic salmon were genotyped in individuals from four F2 families that were selected from the ends of the flesh colour distribution. Genotyping of 23 additional loci and two additional families was performed on a number of linkage groups harbouring putative QTL. QTL analysis was performed using a line-cross model assuming fixation of alternate QTL alleles and a half-sib model with no assumptions about the number and frequency of QTL alleles in the founder populations.

Results

A moderate to strong phenotypic correlation was found between colour, length and weight traits. In total, 13 genome-wide significant QTL were detected for all traits using the line-cross model, including three genome-wide significant QTL for flesh colour (Chr 6, Chr 26 and Chr 4). In addition, 32 suggestive QTL were detected (chromosome-wide P < 0.05). Using the half-sib model, six genome-wide significant QTL were detected for all traits, including two for flesh colour (Chr 26 and Chr 4) and 41 suggestive QTL were detected (chromosome-wide P < 0.05). Based on the half-sib analysis, these two genome-wide significant QTL for flesh colour explained 24% of the phenotypic variance for this trait.

Conclusions

A large number of significant and suggestive QTL for flesh colour and growth traits were found in an F2 population of Atlantic salmon. Chr 26 and Chr 4 presented the strongest evidence for significant QTL affecting flesh colour, while Chr 10, Chr 5, and Chr 4 presented the strongest evidence for significant QTL affecting growth traits (length and weight). These QTL could be strong candidates for use in marker-assisted selection and provide a starting point for further characterisation of the genetic components underlying flesh colour and growth.     

Performance of Polygenic Scores for Predicting Phobic Anxiety

Context

Anxiety disorders are common, with a lifetime prevalence of 20% in the U.S., and are responsible for substantial burdens of disability, missed work days and health care utilization. To date, no causal genetic variants have been identified for anxiety, anxiety disorders, or related traits.

Objective

To investigate whether a phobic anxiety symptom score was associated with 3 alternative polygenic risk scores, derived from external genome-wide association studies of anxiety, an internally estimated agnostic polygenic score, or previously identified candidate genes.

Design

Longitudinal follow-up study. Using linear and logistic regression we investigated whether phobic anxiety was associated with polygenic risk scores derived from internal, leave-one out genome-wide association studies, from 31 candidate genes, and from out-of-sample genome-wide association weights previously shown to predict depression and anxiety in another cohort.

Setting and Participants

Study participants (n = 11,127) were individuals from the Nurses'' Health Study and Health Professionals Follow-up Study.

Main Outcome Measure

Anxiety symptoms were assessed via the 8-item phobic anxiety scale of the Crown Crisp Index at two time points, from which a continuous phenotype score was derived.

Results

We found no genome-wide significant associations with phobic anxiety. Phobic anxiety was also not associated with a polygenic risk score derived from the genome-wide association study beta weights using liberal p-value thresholds; with a previously published genome-wide polygenic score; or with a candidate gene risk score based on 31 genes previously hypothesized to predict anxiety.

Conclusion

There is a substantial gap between twin-study heritability estimates of anxiety disorders ranging between 20–40% and heritability explained by genome-wide association results. New approaches such as improved genome imputations, application of gene expression and biological pathways information, and incorporating social or environmental modifiers of genetic risks may be necessary to identify significant genetic predictors of anxiety.     

Introgression of a major QTL from an inferior into a superior population using genomic selection

Background

Selection schemes aiming at introgressing genetic material from a donor into a recipient line may be performed by backcross-breeding programs combined with selection to preserve the favourable characteristics of the donor population. This stochastic simulation study investigated whether genomic selection can be effective in preserving a major quantitative trait locus (QTL) allele from a donor line during the backcrossing phase.

Methods

In a simulation study, two fish populations were generated: a recipient line selected for a production trait and a donor line characterized by an enhanced level of disease resistance. Both traits were polygenic, but one major QTL affecting disease resistance was segregating only within the donor line. Backcrossing was combined with three types of selection (for total merit index) among the crossbred individuals: classical selection, genomic selection using genome-wide dense marker maps, and gene-assisted genomic selection. It was assumed that production could be observed directly on the selection candidates, while disease resistance had to be inferred from tested sibs of the selection candidates.

Results

Classical selection was inefficient in preserving the target QTL through the backcrossing phase. In contrast, genomic selection (without specific knowledge of the target QTL) was usually effective in preserving the target QTL, and had higher genetic response to selection, especially for disease resistance. Compared with pure genomic selection, gene-assisted selection had an advantage with respect to disease resistance (28–40% increase in genetic gain) and acted as an extra precaution against loss of the target QTL. However, for total merit index the advantage of gene-assisted genomic selection over genomic selection was lower (4–5% increase in genetic gain).

Conclusion

Substantial differences between introgression programs using classical and genomic selection were observed, and the former was generally inferior with respect to both genetic gain and the ability to preserve the target QTL. Combining genomic selection with gene-assisted selection for the target QTL acted as an extra precaution against loss of the target QTL and gave additional genetic gain for disease resistance. However, the effect on total merit index was limited.     

The importance of identity-by-state information for the accuracy of genomic selection

Background

It is commonly assumed that prediction of genome-wide breeding values in genomic selection is achieved by capitalizing on linkage disequilibrium between markers and QTL but also on genetic relationships. Here, we investigated the reliability of predicting genome-wide breeding values based on population-wide linkage disequilibrium information, based on identity-by-descent relationships within the known pedigree, and to what extent linkage disequilibrium information improves predictions based on identity-by-descent genomic relationship information.

Methods

The study was performed on milk, fat, and protein yield, using genotype data on 35 706 SNP and deregressed proofs of 1086 Italian Brown Swiss bulls. Genome-wide breeding values were predicted using a genomic identity-by-state relationship matrix and a genomic identity-by-descent relationship matrix (averaged over all marker loci). The identity-by-descent matrix was calculated by linkage analysis using one to five generations of pedigree data.

Results

We showed that genome-wide breeding values prediction based only on identity-by-descent genomic relationships within the known pedigree was as or more reliable than that based on identity-by-state, which implicitly also accounts for genomic relationships that occurred before the known pedigree. Furthermore, combining the two matrices did not improve the prediction compared to using identity-by-descent alone. Including different numbers of generations in the pedigree showed that most of the information in genome-wide breeding values prediction comes from animals with known common ancestors less than four generations back in the pedigree.

Conclusions

Our results show that, in pedigreed breeding populations, the accuracy of genome-wide breeding values obtained by identity-by-descent relationships was not improved by identity-by-state information. Although, in principle, genomic selection based on identity-by-state does not require pedigree data, it does use the available pedigree structure. Our findings may explain why the prediction equations derived for one breed may not predict accurate genome-wide breeding values when applied to other breeds, since family structures differ among breeds.     

Different models of genetic variation and their effect on genomic evaluation

Background

The theory of genomic selection is based on the prediction of the effects of quantitative trait loci (QTL) in linkage disequilibrium (LD) with markers. However, there is increasing evidence that genomic selection also relies on "relationships" between individuals to accurately predict genetic values. Therefore, a better understanding of what genomic selection actually predicts is relevant so that appropriate methods of analysis are used in genomic evaluations.

Methods

Simulation was used to compare the performance of estimates of breeding values based on pedigree relationships (Best Linear Unbiased Prediction, BLUP), genomic relationships (gBLUP), and based on a Bayesian variable selection model (Bayes B) to estimate breeding values under a range of different underlying models of genetic variation. The effects of different marker densities and varying animal relationships were also examined.

Results

This study shows that genomic selection methods can predict a proportion of the additive genetic value when genetic variation is controlled by common quantitative trait loci (QTL model), rare loci (rare variant model), all loci (infinitesimal model) and a random association (a polygenic model). The Bayes B method was able to estimate breeding values more accurately than gBLUP under the QTL and rare variant models, for the alternative marker densities and reference populations. The Bayes B and gBLUP methods had similar accuracies under the infinitesimal model.

Conclusions

Our results suggest that Bayes B is superior to gBLUP to estimate breeding values from genomic data. The underlying model of genetic variation greatly affects the predictive ability of genomic selection methods, and the superiority of Bayes B over gBLUP is highly dependent on the presence of large QTL effects. The use of SNP sequence data will outperform the less dense marker panels. However, the size and distribution of QTL effects and the size of reference populations still greatly influence the effectiveness of using sequence data for genomic prediction.     

Combined detection and introgression of QTL in outbred populations

Background

Detecting a QTL is only the first step in genetic improvement programs. When a QTL with desirable characteristics is found, e.g. in a wild or unimproved population, it may be interesting to introgress the detected QTL into the commercial population. One approach to shorten the time needed for introgression is to combine both QTL identification and introgression, into a single step. This combines the strengths of fine mapping and backcrossing and paves the way for introgression of desirable but unknown QTL into recipient animal and plant lines.

Methods

The method consisting in combining QTL mapping and gene introgression has been extended from inbred to outbred populations in which QTL allele frequencies vary both in recipient and donor lines in different scenarios and for which polygenic effects are included in order to model background genes. The effectiveness of the combined QTL detection and introgression procedure was evaluated by simulation through four backcross generations.

Results

The allele substitution effect is underestimated when the favourable QTL allele is not fixed in the donor line. This underestimation is proportional to the frequency differences of the favourable QTL allele between the lines. In most scenarios, the estimates of the QTL location are unbiased and accurate. The retained donor chromosome segment and linkage drag are similar to expected values from other published studies.

Conclusions

In general, our results show that it is possible to combine QTL detection and introgression even in outbred species. Separating QTL mapping and introgression processes is often thought to be longer and more costly. However, using a combined process saves at least one generation. With respect to the linkage drag and obligatory drag, the results of the combined detection and introgression scheme are very similar to those of traditional introgression schemes.     

A gene-derived SNP-based high resolution linkage map of carrot including the location of QTL conditioning root and leaf anthocyanin pigmentation

Background

Purple carrots accumulate large quantities of anthocyanins in their roots and leaves. These flavonoid pigments possess antioxidant activity and are implicated in providing health benefits. Informative, saturated linkage maps associated with well characterized populations segregating for anthocyanin pigmentation have not been developed. To investigate the genetic architecture conditioning anthocyanin pigmentation we scored root color visually, quantified root anthocyanin pigments by high performance liquid chromatography in segregating F₂, F₃ and F₄ generations of a mapping population, mapped quantitative trait loci (QTL) onto a dense gene-derived single nucleotide polymorphism (SNP)-based linkage map, and performed comparative trait mapping with two unrelated populations.

Results

Root pigmentation, scored visually as presence or absence of purple coloration, segregated in a pattern consistent with a two gene model in an F₂, and progeny testing of F₃-F₄ families confirmed the proposed genetic model. Purple petiole pigmentation was conditioned by a single dominant gene that co-segregates with one of the genes conditioning root pigmentation. Root total pigment estimate (RTPE) was scored as the percentage of the root with purple color.All five anthocyanin glycosides previously reported in carrot, as well as RTPE, varied quantitatively in the F₂ population. For the purpose of QTL analysis, a high resolution gene-derived SNP-based linkage map of carrot was constructed with 894 markers covering 635.1 cM with a 1.3 cM map resolution. A total of 15 significant QTL for all anthocyanin pigments and for RTPE mapped to six chromosomes. Eight QTL with the largest phenotypic effects mapped to two regions of chromosome 3 with co-localized QTL for several anthocyanin glycosides and for RTPE. A single dominant gene conditioning anthocyanin acylation was identified and mapped.Comparative mapping with two other carrot populations segregating for purple color indicated that carrot anthocyanin pigmentation is controlled by at least three genes, in contrast to monogenic control reported previously.

Conclusions

This study generated the first high resolution gene-derived SNP-based linkage map in the Apiaceae. Two regions of chromosome 3 with co-localized QTL for all anthocyanin pigments and for RTPE, largely condition anthocyanin accumulation in carrot roots and leaves. Loci controlling root and petiole anthocyanin pigmentation differ across diverse carrot genetic backgrounds.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1118) contains supplementary material, which is available to authorized users.     

A novel sex-determining QTL in Nile tilapia (Oreochromis niloticus)

Background

Fish species often exhibit significant sexual dimorphism for commercially important traits. Accordingly, the control of phenotypic sex, and in particular the production of monosex cultures, is of particular interest to the aquaculture industry. Sex determination in the widely farmed Nile tilapia (Oreochromis niloticus) is complex, involving genomic regions on at least three chromosomes (chromosomes 1, 3 and 23) and interacting in certain cases with elevated early rearing temperature as well. Thus, sex ratios may vary substantially from 50%.

Results

This study focused on mapping sex-determining quantitative trait loci (QTL) in families with skewed sex ratios. These included four families that showed an excess of males (male ratio varied between 64% and 93%) when reared at standard temperature (28°C) and a fifth family in which an excess of males (96%) was observed when fry were reared at 36°C for ten days from first feeding. All the samples used in the current study were genotyped for two single-nucleotide polymorphisms (rs397507167 and rs397507165) located in the expected major sex-determining region in linkage group 1 (LG 1). The only misassigned individuals were phenotypic males with the expected female genotype, suggesting that those offspring had undergone sex-reversal with respect to the major sex-determining locus. We mapped SNPs identified from double digest Restriction-site Associated DNA (ddRAD) sequencing in these five families. Three genetic maps were constructed consisting of 641, 175 and 1,155 SNPs from the three largest families. QTL analyses provided evidence for a novel genome-wide significant QTL in LG 20. Evidence was also found for another sex-determining QTL in the fifth family, in the proximal region of LG 1.

Conclusions

Overall, the results from this study suggest that these previously undetected QTLs are involved in sex determination in the Nile tilapia, causing sex reversal (masculinisation) with respect to the XX genotype at the major sex-determining locus in LG 1.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1383-x) contains supplementary material, which is available to authorized users.     

Genome-wide QTL analysis of meat quality-related traits in a large F2 intercross between Landrace and Korean native pigs

Background

We conducted a genome-wide linkage analysis to identify quantitative trait loci (QTL) that influence meat quality-related traits in a large F₂ intercross between Landrace and Korean native pigs. Thirteen meat quality-related traits of the m. longissimus lumborum et thoracis were measured in more than 830 F₂ progeny. All these animals were genotyped with 173 microsatellite markers located throughout the pig genome, and the GridQTL program based on the least squares regression model was used to perform the QTL analysis.

Results

We identified 23 genome-wide significant QTL in eight chromosome regions (SSC1, 2, 6, 7, 9, 12, 13, and 16) (SSC for Sus Scrofa) and detected 51 suggestive QTL in the 17 chromosome regions. QTL that affect 10 meat quality traits were detected on SSC12 and were highly significant at the genome-wide level. In particular, the QTL with the largest effect affected crude fat percentage and explained 22.5% of the phenotypic variance (F-ratio = 278.0 under the additive model, nominal P = 5.5 × 10⁻⁵⁵). Interestingly, the QTL on SSC12 that influenced meat quality traits showed an obvious trend for co-localization.

Conclusions

Our results confirm several previously reported QTL. In addition, we identified novel QTL for meat quality traits, which together with the associated positional candidate genes improve the knowledge on the genetic structure that underlies genetic variation for meat quality traits in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0080-6) contains supplementary material, which is available to authorized users.     

Estimation of (co)variances for genomic regions of flexible sizes: application to complex infectious udder diseases in dairy cattle

Background

Multi-trait genomic models in a Bayesian context can be used to estimate genomic (co)variances, either for a complete genome or for genomic regions (e.g. per chromosome) for the purpose of multi-trait genomic selection or to gain further insight into the genomic architecture of related traits such as mammary disease traits in dairy cattle.

Methods

Data on progeny means of six traits related to mastitis resistance in dairy cattle (general mastitis resistance and five pathogen-specific mastitis resistance traits) were analyzed using a bivariate Bayesian SNP-based genomic model with a common prior distribution for the marker allele substitution effects and estimation of the hyperparameters in this prior distribution from the progeny means data. From the Markov chain Monte Carlo samples of the allele substitution effects, genomic (co)variances were calculated on a whole-genome level, per chromosome, and in regions of 100 SNP on a chromosome.

Results

Genomic proportions of the total variance differed between traits. Genomic correlations were lower than pedigree-based genetic correlations and they were highest between general mastitis and pathogen-specific traits because of the part-whole relationship between these traits. The chromosome-wise genomic proportions of the total variance differed between traits, with some chromosomes explaining higher or lower values than expected in relation to chromosome size. Few chromosomes showed pleiotropic effects and only chromosome 19 had a clear effect on all traits, indicating the presence of QTL with a general effect on mastitis resistance. The region-wise patterns of genomic variances differed between traits. Peaks indicating QTL were identified but were not very distinctive because a common prior for the marker effects was used. There was a clear difference in the region-wise patterns of genomic correlation among combinations of traits, with distinctive peaks indicating the presence of pleiotropic QTL.

Conclusions

The results show that it is possible to estimate, genome-wide and region-wise genomic (co)variances of mastitis resistance traits in dairy cattle using multivariate genomic models.     

Limits on the reproducibility of marker associations with southern leaf blight resistance in the maize nested association mapping population

Background

A previous study reported a comprehensive quantitative trait locus (QTL) and genome wide association study (GWAS) of southern leaf blight (SLB) resistance in the maize Nested Association Mapping (NAM) panel. Since that time, the genomic resources available for such analyses have improved substantially. An updated NAM genetic linkage map has a nearly six-fold greater marker density than the previous map and the combined SNPs and read-depth variants (RDVs) from maize HapMaps 1 and 2 provided 28.5 M genomic variants for association analysis, 17 fold more than HapMap 1. In addition, phenotypic values of the NAM RILs were re-estimated to account for environment-specific flowering time covariates and a small proportion of lines were dropped due to genotypic data quality problems. Comparisons of original and updated QTL and GWAS results confound the effects of linkage map density, GWAS marker density, population sample size, and phenotype estimates. Therefore, we evaluated the effects of changing each of these parameters individually and in combination to determine their relative impact on marker-trait associations in original and updated analyses.

Results

Of the four parameters varied, map density caused the largest changes in QTL and GWAS results. The updated QTL model had better cross-validation prediction accuracy than the previous model. Whereas joint linkage QTL positions were relatively stable to input changes, the residual values derived from those QTL models (used as inputs to GWAS) were more sensitive, resulting in substantial differences between GWAS results. The updated NAM GWAS identified several candidate genes consistent with previous QTL fine-mapping results.

Conclusions

The highly polygenic nature of resistance to SLB complicates the identification of causal genes. Joint linkage QTL are relatively stable to perturbations of data inputs, but their resolution is generally on the order of tens or more Mbp. GWAS associations have higher resolution, but lower power due to stringent thresholds designed to minimize false positive associations, resulting in variability of detection across studies. The updated higher density linkage map improves QTL estimation and, along with a much denser SNP HapMap, greatly increases the likelihood of detecting SNPs in linkage with causal variants. We recommend use of the updated genetic resources and results but emphasize the limited repeatability of small-effect associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1068) contains supplementary material, which is available to authorized users.     

Linear models for joint association and linkage QTL mapping

Background

Populational linkage disequilibrium and within-family linkage are commonly used for QTL mapping and marker assisted selection. The combination of both results in more robust and accurate locations of the QTL, but models proposed so far have been either single marker, complex in practice or well fit to a particular family structure.

Results

We herein present linear model theory to come up with additive effects of the QTL alleles in any member of a general pedigree, conditional to observed markers and pedigree, accounting for possible linkage disequilibrium among QTLs and markers. The model is based on association analysis in the founders; further, the additive effect of the QTLs transmitted to the descendants is a weighted (by the probabilities of transmission) average of the substitution effects of founders'' haplotypes. The model allows for non-complete linkage disequilibrium QTL-markers in the founders. Two submodels are presented: a simple and easy to implement Haley-Knott type regression for half-sib families, and a general mixed (variance component) model for general pedigrees. The model can use information from all markers. The performance of the regression method is compared by simulation with a more complex IBD method by Meuwissen and Goddard. Numerical examples are provided.

Conclusion

The linear model theory provides a useful framework for QTL mapping with dense marker maps. Results show similar accuracies but a bias of the IBD method towards the center of the region. Computations for the linear regression model are extremely simple, in contrast with IBD methods. Extensions of the model to genomic selection and multi-QTL mapping are straightforward.     

Detection of multiple quantitative trait loci and their pleiotropic effects in outbred pig populations

Background

Simultaneous detection of multiple QTLs (quantitative trait loci) may allow more accurate estimation of genetic effects. We have analyzed outbred commercial pig populations with different single and multiple models to clarify their genetic properties and in addition, we have investigated pleiotropy among growth and obesity traits based on allelic correlation within a gamete.

Methods

Three closed populations, (A) 427 individuals from a Yorkshire and Large White synthetic breed, (B) 547 Large White individuals and (C) 531 Large White individuals, were analyzed using a variance component method with one-QTL and two-QTL models. Six markers on chromosome 4 and five to seven markers on chromosome 7 were used.

Results

Population A displayed a high test statistic for the fat trait when applying the two-QTL model with two positions on two chromosomes. The estimated heritabilities for polygenic effects and for the first and second QTL were 19%, 17% and 21%, respectively. The high correlation of the estimated allelic effect on the same gamete and QTL test statistics suggested that the two separate QTL which were detected on different chromosomes both have pleiotropic effects on the two fat traits. Analysis of population B using the one-QTL model for three fat traits found a similar peak position on chromosome 7. Allelic effects of three fat traits from the same gamete were highly correlated suggesting the presence of a pleiotropic QTL. In population C, three growth traits also displayed similar peak positions on chromosome 7 and allelic effects from the same gamete were correlated.

Conclusion

Detection of the second QTL in a model reduced the polygenic heritability and should improve accuracy of estimated heritabilities for both QTLs.     

A high-density SNP genetic linkage map for the silver-lipped pearl oyster,Pinctada maxima: a valuable resource for gene localisation and marker-assisted selection

Background

The silver-lipped pearl oyster, Pinctada maxima, is an important tropical aquaculture species extensively farmed for the highly sought "South Sea" pearls. Traditional breeding programs have been initiated for this species in order to select for improved pearl quality, but many economic traits under selection are complex, polygenic and confounded with environmental factors, limiting the accuracy of selection. The incorporation of a marker-assisted selection (MAS) breeding approach would greatly benefit pearl breeding programs by allowing the direct selection of genes responsible for pearl quality. However, before MAS can be incorporated, substantial genomic resources such as genetic linkage maps need to be generated. The construction of a high-density genetic linkage map for P. maxima is not only essential for unravelling the genomic architecture of complex pearl quality traits, but also provides indispensable information on the genome structure of pearl oysters.

Results

A total of 1,189 informative genome-wide single nucleotide polymorphisms (SNPs) were incorporated into linkage map construction. The final linkage map consisted of 887 SNPs in 14 linkage groups, spans a total genetic distance of 831.7 centimorgans (cM), and covers an estimated 96% of the P. maxima genome. Assessment of sex-specific recombination across all linkage groups revealed limited overall heterochiasmy between the sexes (i.e. 1.15:1 F/M map length ratio). However, there were pronounced localised differences throughout the linkage groups, whereby male recombination was suppressed near the centromeres compared to female recombination, but inflated towards telomeric regions. Mean values of LD for adjacent SNP pairs suggest that a higher density of markers will be required for powerful genome-wide association studies. Finally, numerous nacre biomineralization genes were localised providing novel positional information for these genes.

Conclusions

This high-density SNP genetic map is the first comprehensive linkage map for any pearl oyster species. It provides an essential genomic tool facilitating studies investigating the genomic architecture of complex trait variation and identifying quantitative trait loci for economically important traits useful in genetic selection programs within the P. maxima pearling industry. Furthermore, this map provides a foundation for further research aiming to improve our understanding of the dynamic process of biomineralization, and pearl oyster evolution and synteny.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-810) contains supplementary material, which is available to authorized users.     

Detecting parent of origin and dominant QTL in a two-generation commercial poultry pedigree using variance component methodology

Introduction

Variance component QTL methodology was used to analyse three candidate regions on chicken chromosomes 1, 4 and 5 for dominant and parent-of-origin QTL effects. Data were available for bodyweight and conformation score measured at 40 days from a two-generation commercial broiler dam line. One hundred dams were nested in 46 sires with phenotypes and genotypes on 2708 offspring. Linear models were constructed to simultaneously estimate fixed, polygenic and QTL effects. Different genetic models were compared using likelihood ratio test statistics derived from the comparison of full with reduced or null models. Empirical thresholds were derived by permutation analysis.

Results

Dominant QTL were found for bodyweight on chicken chromosome 4 and for bodyweight and conformation score on chicken chromosome 5. Suggestive evidence for a maternally expressed QTL for bodyweight and conformation score was found on chromosome 1 in a region corresponding to orthologous imprinted regions in the human and mouse.

Conclusion

Initial results suggest that variance component analysis can be applied within commercial populations for the direct detection of segregating dominant and parent of origin effects.