Dirk-Jan de Koning

Genomewide Rapid Association Using Mixed Model and Regression: A Fast and Simple Method For Genomewide Pedigree-Based Quantitative Trait Loci Association Analysis

For pedigree-based quantitative trait loci (QTL) association analysis, a range of methods utilizing within-family variation such as transmission-disequilibrium test (TDT)-based methods have been developed. In scenarios where stratification is not a concern, methods exploiting between-family variation in addition to within-family variation, such as the measured genotype (MG) approach, have greater power. Application of MG methods can be computationally demanding (especially for large pedigrees), making genomewide scans practically infeasible. Here we suggest a novel approach for genomewide pedigree-based quantitative trait loci (QTL) association analysis: genomewide rapid association using mixed model and regression (GRAMMAR). The method first obtains residuals adjusted for family effects and subsequently analyzes the association between these residuals and genetic polymorphisms using rapid least-squares methods. At the final step, the selected polymorphisms may be followed up with the full measured genotype (MG) analysis. In a simulation study, we compared type 1 error, power, and operational characteristics of the proposed method with those of MG and TDT-based approaches. For moderately heritable (30%) traits in human pedigrees the power of the GRAMMAR and the MG approaches is similar and is much higher than that of TDT-based approaches. When using tabulated thresholds, the proposed method is less powerful than MG for very high heritabilities and pedigrees including large sibships like those observed in livestock pedigrees. However, there is little or no difference in empirical power of MG and the proposed method. In any scenario, GRAMMAR is much faster than MG and enables rapid analysis of hundreds of thousands of markers.

Fine mapping and imprinting analysis for fatness trait QTLs in pigs

Abstract. Quantitative trait loci (QTL) for fatness traits were reported recently in an experimental Meishan × Large White and Landrace F2 cross. To further investigate the regions on pig Chr 2 (SSC2), SSC4, and SSC7, 25 additional markers from these regions were typed on 800 animals (619 F2 animals, their F1 parents, and F0 grandfathers). Compared with the published maps, a modified order of markers was observed for SSC4 and SSC7. QTL analyses were performed both within the half-sib families as well as across families (line cross). Furthermore, a QTL model accounting for imprinting effects was tested. Information content could be increased considerably on all three chromosomes. Evidence for the backfat thickness QTL on SSC7 was increased, and the location could be reduced to a 33-cM confidence interval. The QTL for intramuscular fat on SSC4 could not be detected in this half-sib analysis, whereas under the line cross model a suggestive QTL on a different position on SSC4 was detected. For SSC2, in the half-sib analysis, a suggestive QTL for backfat thickness was detected with the best position at 26 cM. Imprinting analysis, however, revealed a genome-wise, significant, paternally expressed QTL on SSC2 with the best position at 63 cM. Our results suggest that this QTL is different from the previously reported paternally expressed QTL for muscle mass and fat deposition on the distal tip of SSC2p.

Best Linear Unbiased Prediction of Genomic Breeding Values Using a Trait-Specific Marker-Derived Relationship Matrix

Background With the availability of high density whole-genome single nucleotide polymorphism chips, genomic selection has become a promising method to estimate genetic merit with potentially high accuracy for animal, plant and aquaculture species of economic importance. With markers covering the entire genome, genetic merit of genotyped individuals can be predicted directly within the framework of mixed model equations, by using a matrix of relationships among individuals that is derived from the markers. Here we extend that approach by deriving a marker-based relationship matrix specifically for the trait of interest. Methodology/Principal Findings In the framework of mixed model equations, a new best linear unbiased prediction (BLUP) method including a trait-specific relationship matrix (TA) was presented and termed TABLUP. The TA matrix was constructed on the basis of marker genotypes and their weights in relation to the trait of interest. A simulation study with 1,000 individuals as the training population and five successive generations as candidate population was carried out to validate the proposed method. The proposed TABLUP method outperformed the ridge regression BLUP (RRBLUP) and BLUP with realized relationship matrix (GBLUP). It performed slightly worse than BayesB with an accuracy of 0.79 in the standard scenario. Conclusions/Significance The proposed TABLUP method is an improvement of the RRBLUP and GBLUP method. It might be equivalent to the BayesB method but it has additional benefits like the calculation of accuracies for individual breeding values. The results also showed that the TA-matrix performs better in predicting ability than the classical numerator relationship matrix and the realized relationship matrix which are derived solely from pedigree or markers without regard to the trait. This is because the TA-matrix not only accounts for the Mendelian sampling term, but also puts the greater emphasis on those markers that explain more of the genetic variance in the trait.

QTL for resistance to summer mortality and OsHV-1 load in the Pacific oyster (Crassostrea gigas).

Summer mortality is a phenomenon severely affecting the aquaculture production of the Pacific oyster (Crassostrea gigas). Although its causal factors are complex, resistance to mortality has been described as a highly heritable trait, and several pathogens including the virus Ostreid Herpes virus type 1 (OsHV-1) have been associated with this phenomenon. A QTL analysis for survival of summer mortality and OsHV-1 load, estimated using real-time PCR, was performed using five F(2) full-sib families resulting from a divergent selection experiment for resistance to summer mortality. A consensus linkage map was built using 29 SNPs and 51 microsatellite markers. Five significant QTL were identified and assigned to linkage groups V, VI, VII and IX. Analysis of single full-sib families revealed differential QTL segregation between families. QTL for the two-recorded traits presented very similar locations, highlighting the interest of further study of their respective genetic controls. These QTL show substantial genetic variation in resistance to summer mortality, and present new opportunities for selection for resistance to OsHV-1.

QTL for body weight, morphometric traits and stress response in European sea bass Dicentrarchus labrax

Natural mating and mass spawning in the European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) complicate genetic studies and the implementation of selective breeding schemes. We utilized a two-step experimental design for detecting QTL in mass-spawning species: 2122 offspring from natural mating between 57 parents (22 males, 34 females and one missing) phenotyped for body weight, eight morphometric traits and cortisol levels, had been previously assigned to parents based on genotypes of 31 DNA microsatellite markers. Five large full-sib families (five sires and two dams) were selected from the offspring (570 animals), which were genotyped with 67 additional markers. A new genetic map was compiled, specific to our population, but based on the previously published map. QTL mapping was performed with two methods: half-sib regression analysis (paternal and maternal) and variance component analysis accounting for all family relationships. Two significant QTL were found for body weight on linkage group 4 and 6, six significant QTL for morphometric traits on linkage groups 1B, 4, 6, 7, 15 and 23 and three suggestive QTL for stress response on linkage groups 3, 14 and 23. The QTL explained between 8% and 38% of phenotypic variance. The results are the first step towards identifying genes involved in economically important traits like body weight and stress response in European sea bass.

Does genomic selection have a future in plant breeding

Plant breeding largely depends on phenotypic selection in plots and only for some, often disease-resistance-related traits, uses genetic markers. The more recently developed concept of genomic selection, using a black box approach with no need of prior knowledge about the effect or function of individual markers, has also been proposed as a great opportunity for plant breeding. Several empirical and theoretical studies have focused on the possibility to implement this as a novel molecular method across various species. Although we do not question the potential of genomic selection in general, in this Opinion, we emphasize that genomic selection approaches from dairy cattle breeding cannot be easily applied to complex plant breeding.

The X chromosome harbors quantitative trait loci for backfat thickness and intramuscular fat content in pigs.

Genetic markers have been used in experimental crosses in pigs todissect genetic variation in quantitative traits (e.g., Andersson et al.1994; Knott et al. 1998; Rohrer and Keele 1998). We reportedrecently on the search for autosomal quantitative trait loci (QTL)for fatness traits in an experimental intercross between the obeseChinese Meishan breed and lean Dutch White production lines (DeKoning et al. 1999; Rattink et al. 2000). In this paper, the analysisto look for QTL on the X Chromosome (Chr) in this cross ispresented.Briefly, 19 Meishan boars were mated to 120 sows of Whitelines from five Dutch breeding companies. From the F

A QUANTITATIVE TRAIT LOCUS FOR LIVE WEIGHT MAPS TO BOVINE CHROMOSOME 23

Abstract. A multiple-marker mapping approach was used to search for quantitative trait loci (QTLs) affecting production, health, and fertility traits in Finnish Ayrshire dairy cattle. As part of a whole-genome scan, altogether 469 bulls were genotyped for six microsatellite loci in 12 families on Chromosome (Chr) 23. Both multiple-marker interval mapping with regression and maximum-likelihood methods were applied with a granddaughter design. Eighteen traits, belonging to 11 trait groups, were included in the analysis. One QTL exceeded experiment level and one QTL genome level significance thresholds. Across-families analysis provided strong evidence (Pexperiment= 0.0314) for a QTL affecting live weight. The QTL for live weight maps between markers BM1258 and BoLA DRBP1. A QTL significant at genome level (Pgenome= 0.0087) was mapped for veterinary treatment, and the putative QTL probably affects susceptibility to milk fever or ketosis. In addition, three traits exceeded the chromosome 5% significance threshold: protein percentage of milk, calf mortality (sire), and milking speed. In within-family analyses, protein percentage was associated with markers in one family (LOD score = 4.5).

Accuracy of genomic prediction using low-density marker panels.

Genomic selection has been widely implemented in national and international genetic evaluation in the dairy cattle industry, because of its potential advantages over traditional selection methods and the availability of commercial high-density (HD) single nucleotide polymorphism (SNP) panels. However, this method may not be cost-effective for cow selection and for other livestock species, because the cost of HD SNP panels is still relatively high. One possible solution that can enable other species to benefit from this promising method is genomic selection with low-density (LD) SNP panels. In this simulation study, LD SNP panels designed with different strategies and different SNP densities were compared. The effects of number of quantitative trait loci, heritability, and effective population size were evaluated in the framework of genomic selection with LD SNP panels. Methodologies of Bayesian variable selection; BLUP with a trait-specific, marker-derived relationship matrix; and BLUP with a realized relationship matrix were employed to predict genomic estimated breeding values with both HD and LD SNP panels. Up to 95% of accuracy obtained by using an HD panel can be obtained by using only a small proportion of markers. The LD panel with markers selected on the basis of their effects always performs better than the LD panel with evenly spaced markers. Both the genetic architecture of the trait and the effective population size have a significant effect on the performance of the LD panels. We concluded that, to implement genomic selection with LD panels, a training population of sufficient size and genotyped with an HD panel is necessary. The trade-off between the LD panels with evenly spaced markers and selected markers must be considered, which depends on the number of target traits in a breeding program and the genetic architecture of these traits. Genomic selection with LD panels could be feasible and cost-effective, though before implementation, a further detailed genetic and economic analysis is recommended.

Controlling false positives in the mapping of epistatic QTL.

This study addresses the poorly explored issue of the control of false positive rate (FPR) in the mapping of pair-wise epistatic quantitative trait loci (QTL). A nested test framework was developed to (1) allow pre-identified QTL to be used directly to detect epistasis in one-dimensional genome scans, (2) to detect novel epistatic QTL pairs in two-dimensional genome scans and (3) to derive genome-wide thresholds through permutation and handle multiple testing. We used large-scale simulations to evaluate the performance of both the one- and two-dimensional approaches in mapping different forms and levels of epistasis and to generate profiles of FPR, power and accuracy to inform epistasis mapping studies. We showed that the nested test framework and genome-wide thresholds were essential to control FPR at the 5% level. The one-dimensional approach was generally more powerful than the two-dimensional approach in detecting QTL-associated epistasis and identified nearly all epistatic pairs detected from the two-dimensional approach. However, only the two-dimensional approach could detect epistatic QTL with weak main effects. Combining the two approaches allowed effective mapping of different forms of epistasis, whereas using the nested test framework kept the FPR under control. This approach provides a good search engine for high-throughput epistasis analyses.