A method to optimize selection on multiple identified quantitative trait loci

A mathematical approach was developed to model and optimize selection on multiple known quantitative trait loci (QTL) and polygenic estimated breeding values in order to maximize a weighted sum of responses to selection over multiple generations. The model allows for linkage between QTL with multiple alleles and arbitrary genetic effects, including dominance, epistasis, and gametic imprinting. Gametic phase disequilibrium between the QTL and between the QTL and polygenes is modeled but polygenic variance is assumed constant. Breeding programs with discrete generations, differential selection of males and females and random mating of selected parents are modeled. Polygenic EBV obtained from best linear unbiased prediction models can be accommodated. The problem was formulated as a multiple-stage optimal control problem and an iterative approach was developed for its solution. The method can be used to develop and evaluate optimal strategies for selection on multiple QTL for a wide range of situations and genetic models.     

Mapping quantitative trait loci with DNA microsatellites in a commercial dairy cattle population

Individual loci affecting economically important traits can be located using genetic linkage between quantitative trait loci and genetic markers. In the ‘granddaughter’ experimental design, heterozygous grandsires and their sons are genotyped for the genetic marker, while the quantitative trait records of the granddaughters are used for statistical analysis. Ten DNA microsatellite markers were used to look for associations with quantitative trait loci affecting milk production traits in seven Israeli Holstein grandsire families. At least 60% more grandsires were heterozygous, and 40% fewer individuals were discarded because of unknown paternal allele origin, as compared with diallelic markers. The effects of paternal alleles for locus D21S4 on kg milk and protein were significant (P < 0.025). The allele substitution effects for sire 783 were 283 kg milk and 5.7 kg protein. For both traits, progeny of sire 783 that inherited allele ‘18’ had higher evaluations than progeny that inherited allele ‘21’. These results were verified by genotyping 151 of his daughters. Thus, the rate of genetic gain for protein production can be increased by selecting progeny of sire 783 carrying allele ‘18’ at this locus.     

Spatially and temporally varying selection on intrapopulation quantitative trait loci for a life history trade-off in Mimulus guttatus

Why do populations remain genetically variable despite strong continuous natural selection? Mutation reconstitutes variation eliminated by selection and genetic drift, but theoretical and experimental studies each suggest that mutation‐selection balance insufficient to explain extant genetic variation in most complex traits. The alternative hypothesis of balancing selection, wherein selection maintains genetic variation, is an aggregate of multiple mechanisms (spatial and temporal heterogeneity in selection, frequency‐dependent selection, antagonistic pleiotropy, etc.). Most of these mechanisms have been demonstrated for Mendelian traits, but there is little comparable data for loci affecting quantitative characters. Here, we report a 3‐year field study of selection on intrapopulation quantitative trait loci (QTL) of flower size, a highly polygenic trait in Mimulus guttatus. The QTL exhibit antagonistic pleiotropy: alleles that increase flower size, reduce viability, but increase fecundity. The magnitude and direction of selection fluctuates yearly and on a spatial scale of metres. This study provides direct evidence of balancing selection mechanisms on QTL of an ecologically relevant trait.     

Efficiency of direct selection on quantitative trait loci for a two-trait breeding objective

Selection on known loci affecting quantitative traits (DSQ) was compared to phenotypic selection index for a single and a two-trait selection objective. Two situations were simulated; a single known quantitative locus, and ten identified loci accounting for all the additive genetic variance. Selection efficiency of DSQ relative to traitbased selection was higher for two-trait selection, than was selection on a single trait with the same heritability. The advantage of DSQ was greater when the traits were negatively correlated. Relative selection efficiency (RSE) for a single locus responsible for 0.1 of the genetic variance was 1.11 with heritabilities of 0.45 and 0.2 and zero genetic and phenotypic correlations between the traits. RSE of DSQ for ten known loci was 1.5 to 1.8 in the first 3 generations of selection, but declined in each subsequent generation. With DSQ most loci reached fixation after 7 generations. Response to trait-based selection continued through generation 15 and approached the response obtained with DSQ after 10 generations. The cumulative genetic response after 10 generations of DSQ was only 93% to 97% of the economically optimum genotype because the less favorable allele reached fixation for some loci, generally those with effects in opposite directions on the two traits.     

DNA fingerprint bands applied to linkage analysis with quantitative trait loci in chickens

Summary
An efficient approach to detect association between quantitative traits and bands of DNA fingerprint patterns uses intra-family tail analysis, which compares fingerprints of DNA mixes from individuals at the two tails of a phenotypic distribution. In analysis of 67 paternal half-sibs of a meat-type chicken family, of 57 sire bands generated by two probes, one sire-specific band (S6–6) was associated with abdominal fat deposition. The band effect was estimated by a linear model analysis to be 0–88 standard deviations, or about 30% of the family mean. The association between band S6–6 and abdominal fat was further examined by testing progeny of paternal half-sibs of the chickens which were used in the tail analysis, establishing genetic linkage between the DNA marker and a genetic locus affecting abdominal fat deposition.     

The efficiency of designs for fine-mapping of quantitative trait loci using combined linkage disequilibrium and linkage

In a simulation study, different designs were compared for efficiency of fine-mapping of QTL. The variance component method for fine-mapping of QTL was used to estimate QTL position and variance components. The design of many families with small size gave a higher mapping resolution than a design with few families of large size. However, the difference is small in half sib designs. The proportion of replicates with the QTL positioned within 3 cM of the true position is 0.71 in the best design, and 0.68 in the worst design applied to 128 animals with a phenotypic record and a QTL explaining 25% of the phenotypic variance. The design of two half sib families each of size 64 was further investigated for a hypothetical population with effective size of 1000 simulated for 6000 generations with a marker density of 0.25 cM and with marker mutation rate 4 × 10^-4 per generation. In mapping using bi-allelic markers, 42~55% of replicated simulations could position QTL within 0.75 cM of the true position whereas this was higher for multi allelic markers (48~76%). The accuracy was lowest (48%) when mutation age was 100 generations and increased to 68% and 76% for mutation ages of 200 and 500 generations, respectively, after which it was about 70% for mutation ages of 1000 generations and older. When effective size was linearly decreasing in the last 50 generations, the accuracy was decreased (56 to 70%). We show that half sib designs that have often been used for linkage mapping can have sufficient information for fine-mapping of QTL. It is suggested that the same design with the same animals for linkage mapping should be used for fine-mapping so gene mapping can be cost effective in livestock populations.     

Expression quantitative trait loci: present and future

The last few years have seen the development of large efforts for the analysis of genome function, especially in the context of genome variation. One of the most prominent directions has been the extensive set of studies on expression quantitative trait loci (eQTLs), namely, the discovery of genetic variants that explain variation in gene expression levels. Such studies have offered promise not just for the characterization of functional sequence variation but also for the understanding of basic processes of gene regulation and interpretation of genome-wide association studies. In this review, we discuss some of the key directions of eQTL research and its implications.     

Mapping of quantitative trait loci on porcine chromosome 4

A F₂ population derived from a cross between European Large White and Chinese Meishan pigs was established in order to study the genetic basis of breed differences for growth and fat traits. Chromosome 4 was chosen for initial study as previous work had revealed quantitative trait loci (QTLs) on this chromosome affected growth and fat traits in a Wild Boar × Large White cross. Individuals in the F₂ population were typed for nine markers spanning a region of approximately 124 c m . We found evidence for QTLs affecting growth between weaning and the end of test (additive effect: 43·4 g/day) and fat depth measured in the mid-back position (additive effect: 1·82 mm). There was no evidence of interactions between the QTLs and sex, grandparents or F₁ sires, suggesting that the detected QTLs were fixed for alternative alleles in the Meishan and Large White breeds. Comparison of locations suggests that these QTLs could be the same as those found in the Wild Boar × Large White cross.     

Model selection for quantitative trait loci mapping in a full-sib family

Statistical methods for mapping quantitative trait loci (QTLs) in full-sib forest trees, in which the number of alleles and linkage phase can vary from locus to locus, are still not well established. Previous studies assumed that the QTL segregation pattern was fixed throughout the genome in a full-sib family, despite the fact that this pattern can vary among regions of the genome. In this paper, we propose a method for selecting the appropriate model for QTL mapping based on the segregation of different types of markers and QTLs in a full-sib family. The QTL segregation patterns were classified into three types: test cross (1:1 segregation), F₂ cross (1:2:1 segregation) and full cross (1:1:1:1 segregation). Akaike’s information criterion (AIC), the Bayesian information criterion (BIC) and the Laplace-empirical criterion (LEC) were used to select the most likely QTL segregation pattern. Simulations were used to evaluate the power of these criteria and the precision of parameter estimates. A Windows-based software was developed to run the selected QTL mapping method. A real example is presented to illustrate QTL mapping in forest trees based on an integrated linkage map with various segregation markers. The implications of this method for accurate QTL mapping in outbred species are discussed.     

Genetic mapping of quantitative trait loci for aseasonal reproduction in sheep

The productivity and economic prosperity of sheep farming could benefit greatly from more effective methods of selection for year-round lambing. Identification of QTL for aseasonal reproduction in sheep could lead to more accurate selection and faster genetic improvement. One hundred and twenty microsatellite markers were genotyped on 159 backcross ewes from a Dorset × East Friesian crossbred pedigree. Interval mapping was undertaken to map the QTL underlying several traits describing aseasonal reproduction including the number of oestrous cycles, maximum level of progesterone prior to breeding, pregnancy status determined by progesterone level, pregnancy status determined by ultrasound, lambing status and number of lambs born. Seven chromosomes (1, 3, 12, 17, 19, 20 and 24) were identified to harbour putative QTL for one or more component traits used to describe aseasonal reproduction. Ovine chromosomes 12, 17, 19 and 24 harbour QTL significant at the 5% chromosome-wide level, chromosomes 3 and 20 harbour QTL that exceeded the threshold at the 1% chromosome-wide level, while the QTL identified on chromosome 1 exceeded the 1% experiment-wide significance level. These results are a first step towards understanding the genetic mechanism of this complex trait and show that variation in aseasonal reproduction is associated with multiple chromosomal regions.     

An integrative approach for the identification of quantitative trait loci

The genetic dissection of complex traits is one of the most difficult and most important challenges facing science today. We discuss here an integrative approach to quantitative trait loci (QTL) mapping in mice. This approach makes use of the wealth of genetic tools available in mice, as well as the recent advances in genome sequence data already available for a number of inbred mouse strains. We have developed mapping strategies that allow a stepwise narrowing of a QTL mapping interval, prioritizing candidate genes for further analysis with the potential of identifying the most probable candidate gene for the given trait. This approach integrates traditional mapping tools, fine mapping tools, sequence-based analysis, bioinformatics and gene expression.     

Sequential sampling in determining linkage between marker loci and quantitative trait loci

Summary As compared to classical, fixed sample size techniques, simulation studies showed that a proposed sequential sampling procedure can provide a substantial decrease (up to 50%, in some cases) in the mean sample size required for the detection of linkage between marker loci and quantitative trait loci. Sequential sampling with truncation set at the required sample size for the non-sequential test, produced a modest further decrease in mean sample size, accompanied by a modest increase in error probabilities. Sequential sampling with observations taken in groups produced a noticeable increase in mean sample size, with a considerable decrease in error probabilities, as compared to straightforward sequential sampling. It is concluded that sequential sampling has a particularly useful application to experiments aimed at investigating the genetics of differences between lines or strains that differ in some single outstanding trait.     

Revisiting the quantitative trait loci for milk production traits on BTA6

A parallel association study was performed in two independent cattle populations based on 41 validated, targeted single nucleotide polymorphisms (SNPs) and four microsatellite markers to re-evaluate the multiple quantitative trait loci (QTL) architecture for milk performance on bovine chromosome 6 (BTA6). Two distinct QTL located in the vicinity of the middle region of BTA6, but differing unambiguously regarding their effects on milk composition and yield traits were validated in the German Holstein population. A highly significant association of the protein variant ABCG2 p.Tyr581Ser with milk composition traits reconfirmed the causative molecular relevance of the ABCG2 gene in QTL region 1, whereas in QTL region 2, significant and tentative associations between gene variants RW070 and RW023 (located in the promoter region and exon 9 of the PPARGC1A gene for milk yield traits) were detected. For the German Fleckvieh population, only RW023 showed a tentative association with milk yield traits, whereas those loci with significant effects in German Holsteins (ABCG2 p.Tyr581Ser, RW070) showed fixed alleles. Even though our new data highlight two variants in the PPARGC1A gene (RW023, RW070) in QTL region 2, based on the results of our study, currently no unequivocal conclusion about the causal background of this QTL affecting milk yield traits can be drawn. Notably, the German Holstein and Fleckvieh populations, known for their divergent degree of dairy type, differ substantially in their allele frequencies for the growth-associated NCAPG p.Ile442Met locus.     

Meta-analysis of results from quantitative trait loci mapping studies on pig chromosome 4

Meta-analysis of results from multiple studies could lead to more precise quantitative trait loci (QTL) position estimates compared to the individual experiments. As the raw data from many different studies are not readily available, the use of results from published articles may be helpful. In this study, we performed a meta-analysis of QTL on chromosome 4 in pig, using data from 25 separate experiments. First, a meta-analysis was performed for individual traits: average daily gain and backfat thickness. Second, a meta-analysis was performed for the QTL of three traits affecting loin yield: loin eye area, carcass length and loin meat weight. Third, 78 QTL were selected from 20 traits that could be assigned to one of three broad categories: carcass, fatness or growth traits. For each analysis, the number of identified meta-QTL was smaller than the number of initial QTL. The reduction in the number of QTL ranged from 71% to 86% compared to the total number before the meta-analysis. In addition, the meta-analysis reduced the QTL confidence intervals by as much as 85% compared to individual QTL estimates. The reduction in the confidence interval was greater when a large number of independent QTL was included in the meta-analysis. Meta-QTL related to growth and fatness were found in the same region as the FAT1 region. Results indicate that the meta-analysis is an efficient strategy to estimate the number and refine the positions of QTL when QTL estimates are available from multiple populations and experiments. This strategy can be used to better target further studies such as the selection of candidate genes related to trait variation.     

Mapping of quantitative trait loci affecting behaviour in swine

Behavioural indices in vertebrates are under genetic control at least to some extent. In spite of significant behavioural problems in farm animals, information on the genetic background of behaviour is sparse. The aim of this study was to map QTL for behavioural indices in swine under healthy conditions and after infection with Sarcocystis miescheriana , as behaviour can be significantly influenced by disease . This well-described parasite model subsequently leads to acute (day 14 p.i.), subclinical (day 28 p.i.) and chronic disease (day 42 p.i.), allowing the study and comparison of the behaviour of pigs under four different states of health or disease. The study was based on a well-described Pietrain/Meishan F₂ family that has recently allowed the detection of QTL for disease resistance. We have mapped six genome-wide significant and 24 chromosome-wide significant QTL for six behavioural indices in swine. Six of these QTL (i.e. 20% of total QTL) showed effects on behavioural traits of the healthy pigs (day 0). Some of them (QTL on SSC11 and 18) lost influence on behavioural activities during disease, while the effects of others (QTL on SSC5, SSC8) partly remained during the whole experiment, although with different effects on the distinct behavioural indices. The disease model has been of high relevance to detect effects of gene loci on behavioural indices. Considering the importance of segregating alleles and environmental conditions that allow the identification of the phenotype, we conclude that there are indeed QTL with interesting effects on behavioural indices in swine.     

Mapping quantitative trait loci regulating chicken body composition traits

Genome scans were conducted on an F₂ resource population derived from intercross of the White Plymouth Rock with the Silkies Fowl to detect QTL affecting chicken body composition traits. The population was genotyped with 129 microsatellite markers and phenotyped for 12 body composition traits on 238 F₂ individuals from 15 full-sib families. In total, 21 genome-wide QTL were found to be responsible for 11 traits, including two newly studied traits of proventriculus weight and shank girth. Three QTL were genome-wide significant: at 499 c m on GGA1 (explained 3.6% of phenotypic variance, P  < 0.01) and 51 c m on GGA5 (explained 3.3% of phenotypic variance, P  < 0.05) for the shank & claw weight and 502 c m on GGA1 (explained 1.4% of phenotypic variance, P  < 0.05) for wing weight. The QTL on GGA1 seemed to have pleiotropic effects, also affecting gizzard weight at 490 c m , shank girth at 489 c m and intestine length at 481 c m . It is suggested that further efforts be made to understand the possible pleiotropic effects of the QTL on GGA1 and that on GGA5 for two shank-related traits.     

Confirmation of quantitative trait loci affecting fatness in chickens

In this report we describe the analysis of an advanced intercross line (AIL) to confirm the quantitative trait locus (QTL) regions found for fatness traits in a previous study. QTL analysis was performed on chromosomes 1, 3, 4, 15, 18, and 27. The AIL was created by random intercrossing in each generation from generation 2 (G₂) onwards until generation 9 (G₉) was reached. QTL for abdominal fat weight (AFW) and/or percentage abdominal fat (AF%) on chromosomes 1, 3 and 27 were confirmed in the G₉ population. In addition, evidence for QTL for body weight at the age of 5 (BW5) and 7 (BW7) weeks and for the percentage of intramuscular fat (IF%) were found on chromosomes 1, 3, 15, and 27. Significant evidence for QTL was detected on chromosome 1 for BW5 and BW7. Suggestive evidence was found on chromosome 1 for AFW, AF% and IF%, on chromosome 15 for BW5, and on chromosome 27 for AF% and IF%. Furthermore, evidence on the chromosome-wise level was found on chromosome 3 for AFW, AF%, and BW7 and on chromosome 27 for BW5. For chromosomes 4 and 18, test statistics did not exceed the significance threshold.     

Confirmation of two quantitative trait loci regions for nematode resistance in commercial British terminal sire breeds

Sheep internal parasites (nematodes) remain a major health challenge and are costly for pasture-based production systems. Most current breeding programmes for nematode resistance are based on indicator traits such as faecal egg counts (FEC), which are costly and laborious to collect. Hence, genetic markers for resistance would be advantageous. However, although some quantitative trait loci (QTL) have been identified, these QTL are often not consistent across breeds and few breeding strategies for nematode resistance in sheep are currently using molecular information. In this study, QTL for nematode resistance on ovine chromosomes (OAR) 3 and 14, previously identified in the Blackface breed, were explored using commercial Suffolk (n = 336) and Texel lambs (n = 879) sampled from terminal sire breeder flocks in the United Kingdom. FEC were used as the indicator trait for nematode resistance, and these were counted separately for Nematodirus and Strongyles genera. Microsatellite markers were used to map the QTL and the data were analysed using interval mapping regression techniques and variance component analysis. QTL for Nematodirus and Strongyles FEC were found to be segregating on OAR3 at 5% chromosome region-wide significance threshold in both Suffolk and Texel sheep, and Nematodirus FEC QTL were segregating on OAR14 in both breeds. In addition, QTL for growth traits were also found to be segregating at 5% chromosome region-wide on OAR3 and OAR14. The confirmation that FEC QTL segregate in the same position in three widely used breeds widens their potential applicability to purebred Blackface, Suffolk and Texel sheep, with benefits likely to be observed in their commercial crossbred progeny.     

A search for quantitative trait loci for ovulation rate in cattle

Seventy-seven polymorphic microsatellites were analysed in offspring of three elite sires that were part of the foundation of an experimental population selected for twinning rate at the US Meat Animal Research Center, Clay Center, Nebraska. All females were assessed for ovulation rate by rectal palpation of corpora lutea over 8–10 consecutive oestrous cycles from approximately 12 to 18 months of age, and associations between ovulation rate and sire allele were examined in each of the three sire groups. A preliminary analysis was performed using selectively genotyped daughters of each sire. Markers found significant or approaching significance were also genotyped in all daughters, sons and granddaughters of these sires. A test of marker associations limited to the granddaughter data provided an independent confirmation of marker effect and significance relative to the initial test with daughter data. Putative ovulation rate quantitative trait loci were detected on chromosomes 7 and 23. Marker UWCA20 on chromosome 7 was associated with an effect in excess of one phenotypic standard deviation and accounted for approximately 10% of phenotypic variation ovulation rate. Marker CYP21 (steroid 21-hydroxylase) on chromosome 23 was associated with an effect of slightly less than half a phenotypic standard deviation and accounted for approximately 4% of phenotypic variation.