Evidence for quantitative trait loci affecting twinning rate in North American Holstein cattle

Twinning in dairy cattle has been associated with many negative health and reproductive events that cause economic loss to the producer. Reports have suggested that twinning rates are increasing and that there may be a positive relationship between milk production and twinning frequency. Putative quantitative trait loci (QTL) for twinning and ovulation rate on bovine chromosomes 5, 7, 19 and 23 have been previously identified in other populations. The objective of this study was to detect and possibly confirm the existence and effects of these QTL in the North American Holstein population. Half-sib families of 20 North American Holstein sires with above average twinning rate predicted transmitting abilities (PTA) comprised the sample population under investigation. Twinning rate PTA values had been estimated from calving data. DNA extracted from semen samples was analysed using 45-61 microsatellite markers across the four chromosomes. Marker heterozygosity of the patriarchs averaged 62%. Evidence of twinning QTL was found in multiple families on chromosomes 5, 7 and 23 and in one family on chromosome 19. Four of the sires formed one three-generation family: one sire and three half-sib sons with sons of their own. This extended family was analysed with additional markers confirming a twinning QTL of significant size on chromosome 5.     

A search for quantitative trait loci for milk production traits on chromosome 6 in Finnish Ayrshire cattle

Cattle chromosome 6 was scanned with 11 markers, ten microsatellites and the casein haplotype, to identify quantitative trait loci (QTLs) affecting the following milk production traits: milk yield, fat percentage, fat yield, protein percentage and protein yield. Twelve Finnish Ayrshire half-sib families with a total of 480 sons were genotyped and used in a grand-daughter design. Interval mapping was performed with a multiple-marker regression approach with a one-QTL and a two-QTL model, and the significance threshold values were determined empirically using a permutation test. Across-family analysis with the one-QTL model revealed an effect on protein percentage (P < 0.05) and on milk yield (P < 0.05). The analysis with the two-QTL model identified significant effects (P < 0.05) on protein percentage, milk yield, and fat yield. Comparing these two cases, the results suggest the existence of two QTLs on chromosome 6 with an effect on milk production traits. One of the QTLs was located around the casein genes. As the other QTL was similar in location and effect to a QTL found previously in Holstein-Friesians, an identity-by-descent approach could be applied to fine map this region.     

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.     

Detection of quantitative trait loci for meat quality traits in cattle

A whole-genome scan was carried out to detect quantitative trait loci (QTL) affecting sensory, organoleptic, physical and chemical properties of meat. The study used phenotypic data from 235 second-generation cross-bred bull calves of a Charolais × Holstein experimental population. Loin muscle samples were evaluated for yield force, intramuscular fat and nitrogen contents, myofibrillar fragmentation index, haem pigment concentration, moisture content and pH at 24 h postmortem. A sensory assessment was performed on grilled loin and roasted silverside joints by trained panellists. A linear regression analysis based on 165 markers revealed 35 QTL at the 5% chromosome-wide significance level (20 for sensory traits and 15 for physical and chemical traits), five of which were highly significant ( F -value: ≥9). The most significant QTL was located on chromosome 6 (with the best likely position at 39 cM) and affected haem pigment concentration. The Holstein allele for this QTL was associated with an increase of 0.53 SD in the haem scores. A QTL for pH_24h was identified on chromosome 14 (at 40 cM) and a QTL for moisture content was identified on chromosome 22 (at 21 cM). Two highly significant QTL were identified for sensory panel-assessed traits: beef odour intensity (grilled sample) on chromosome 10 (at 119 cM), and juiciness (roast sample) on chromosome 16 (at 70 cM). The proportion of phenotypic variance explained by the significant QTL ranged from 3.6% (for nitrogen content on chromosome 10) to 9.5% (for juiciness, roast sample on chromosome 16).     

Identification of quantitative trait loci for growth and carcass composition in cattle

A genomic screening to detect quantitative trait loci (QTL) affecting growth, carcass composition and meat quality traits was pursued. Two hundred nineteen microsatellite markers were genotyped on 176 of 620 (28%) progeny from a Brahman x Angus sire mated to mostly MARC III dams. Selective genotyping, based on retail product yield (%) and fat yield (%), was used to select individuals to be genotyped. Traits included in the study were birth weight (kg), hot carcass weight (kg), retail product yield, fat yield, marbling score (400 = slight00 and 500 = small00), USDA yield grade, and estimated kidney, heart and pelvic fat (%). The QTL were classified as significant when the expected number of false positives (ENFP) was less than 0.05 (F-statistic greater than 17.3), and suggestive when the ENFP was <1 (F-statistic between 10.2 and 17.3). A significant QTL (F = 19; ENFP = 0.02) was detected for marbling score at centimorgan (cM) 54 on chromosome 2. Suggestive QTL were detected for fat yield at 50 cM, for retail product yield at 53 cM, and for USDA yield grade at 63 cM on chromosome 1, for marbling score at 56 cM, for retail product yield at 70 cM, and for estimated kidney, heart and pelvic fat at 79 cM on chromosome 3, for marbling score at 44 cM, for hot carcass weight at 49 cM, and for estimated kidney, heart and pelvic fat at 62 cM on chromosome 16, and for fat yield at 35 cM on chromosome 17. Two suggestive QTL for birth weight were identified, one at 12 cM on chromosome 20 and the other at 56 cM on chromosome 21. An additional suggestive QTL was detected for retail product yield, for fat yield, and for USDA yield grade at 26 cM on chromosome 26. Results presented here represent the initial search for quantitative trait loci in this family. Validation of detected QTL in other populations will be necessary.     

Identification of an ovulation rate QTL in cattle on BTA14 using selective DNA pooling and interval mapping

Increased twinning incidence in beef cattle has the potential to improve production efficiency. However, phenotypic selection for twinning rate is difficult because of the trait's low heritability and the long time interval necessary to collect phenotypic records. Therefore, this trait and the correlated trait of ovulation rate are ideal candidates for marker-assisted selection. The objective of this study was to complete a genome-wide search for ovulation rate quantitative trait loci (QTL) in two related sire families. The families (paternal halfsib sires 839802 and 839803) were from a population of cattle selected for ovulation rate at the USDA Meat Animal Research Center, Clay Center, Nebraska. Putative ovulation rate QTL have previously been identified in the 839802 family on chromosomes 7 and 19; however, marker coverage in the original scan was not complete. This study fills the gaps in marker coverage of the earlier study by adding approximately 60 informative microsatellites to each sire family. Each family was genotyped using selective DNA pooling. Sons and daughters were included in either the high or low pool based on their estimated breeding value deviations from the mid-parent average (EBVMD) for ovulation rate. Approximately 40% (839802) and 26% (839803) of available progeny comprised the high and low pools combined. Pooled typing revealed possible associations (nominal P < 0.05) between ovulation rate and marker genotype for 11 and 15 microsatellites in the 839802 and 839803 families, respectively. Subsequent interval mapping strengthened support for the presence of an ovulation rate QTL on BTA14 (chromosome-wise P < 0.02).     

Mapping of bovine ovulation rate QTL; an analytical approach for three generation pedigrees

Interval mapping was conducted for ovulation rate quantitative trait loci (QTL) using data from two related families from the United States Department of Agriculture, US Meat Animal Research Center twinning cattle herd. Both families are extended, three generation pedigrees from which records of sons, daughters and granddaughters were analysed. Both a method of analysis and results from that analysis are reported herein. Results from one of the two families (839802) were previously reported, but reanalysis here including the second, related family (839803) and a revised statistical model lessens support for the previously reported QTL. Results from interval mapping provided evidence for QTL in regions corresponding to those previously suggested for chromosomes 7 (chromosome-wise P < 0.05) and 19 (chromosome-wise P < 0.01) in the 839802 family, although statistical significance was reduced. In contrast to the previous report, evidence for a chromosome 5 QTL in the same family was greatly reduced while support for a QTL on chromosome 10 increased (chromosome-wise P < 0.01). Analysis of data from the related 839803 family failed to replicate evidence of QTL observed on either chromosome 7 or chromosome 19 in the 839802 family.     

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.     

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.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

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.     

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.     

Mapping quantitative trait loci for cytokines in the pig

Increased disease resistance through improved general immune capacity would be beneficial for the welfare and productivity of farm animals. Cytokines are essential diagnostic parameters in veterinary practice. To identify quantitative trait loci (QTL) for cytokine levels in serum in the pig, Interferon‐gamma (IFN‐γ) and Interleukin 10 (IL‐10) levels and the ratio of IFN‐γ to IL‐10 were measured in a composite pig population, before and after challenge with modified live CSF (classical swine fever) vaccine. Through interval mapping using the variance component approach and the permutation test, 11 QTL (five for IFN‐γ, two for IL‐10 and four for the ratio of IFN‐γ to IL‐10) with significance levels of P < 0.10 were identified, of which five were significant at the P < 0.05 level. The most significant QTL (P < 0.01) was found on chromosome 16, with effect on the ratio of IFN‐γ to IL‐10. Within these QTL regions, a number of known genes were revealed and their potential relationships to the studied traits were discussed. Some of these genes may serve as candidate genes for these traits in swine.     

Joint tests for quantitative trait loci in experimental crosses

Selective genotyping is common because it can increase the expected correlation between QTL genotype and phenotype and thus increase the statistical power of linkage tests (i.e., regression-based tests). Linkage can also be tested by assessing whether the marginal genotypic distribution conforms to its expectation, a marginal-based test. We developed a class of joint tests that, by constraining intercepts in regression-based analyses, capitalize on the information available in both regression-based and marginal-based tests. We simulated data corresponding to the null hypothesis of no QTL effect and the alternative of some QTL effect at the locus for a backcross and an F2 intercross between inbred strains. Regression-based and marginal-based tests were compared to corresponding joint tests. We studied the effects of random sampling, selective sampling from a single tail of the phenotypic distribution, and selective sampling from both tails of the phenotypic distribution. Joint tests were nearly as powerful as all competing alternatives for random sampling and two-tailed selection under both backcross and F2 intercross situations. Joint tests were generally more powerful for one-tailed selection under both backcross and F2 intercross situations. However, joint tests cannot be recommended for one-tailed selective genotyping if segregation distortion is suspected.     

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.     

Sib mating designs for mapping quantitative trait loci

The power to separate the variance of a quantitative trait locus (QTL) from the polygenic variance is determined by the variability of genes identical by descent (IBD) at the QTL. This variability may increase with inbreeding. Selfing, the most extreme form of inbreeding, increases the variability of the IBD value shared by siblings, and thus has a higher efficiency for QTL mapping than random mating. In self-incompatible organisms, sib mating is the closest form of inbreeding. Similar to selfing, sib mating may also increase the power of QTL detection relative to random mating. In this study, we develop an IBD-based method under sib mating designs for QTL mapping. The efficiency of sib mating is then compared with random mating. Monte Carlo simulations show that sib mating designs notably increase the power for QTL detection. When power is intermediate, the power to detect a QTL using full-sib mating is, on average, 7% higher than under random mating. In addition, the IBD-based method proposed in this paper can be used to combine data from multiple families. As a result, the estimated QTL parameters can be applied to a wide statistical inference space relating to the entire reference population. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mapping of quantitative trait loci for carcass traits in a Japanese Black (Wagyu) cattle population

To detect quantitative trait loci (QTL) that influence economically important traits in a purebred Japanese Black cattle population, we performed a preliminary genome-wide scan using 187 microsatellite markers across a paternal half-sib family composed of 258 offspring. We located six QTL at the 1% chromosome-wise level on bovine chromosomes (BTA) 4, 6, 13, 14 and 21. A second screen of these six QTL regions using 138 additional paternal offspring half-sib from the same sire, provided further support for five QTL: carcass weight on BTA14 (22-39 cM), one for rib thickness on BTA6 (27-58 cM) and three for beef marbling score (BMS) on BTA4 (59-67 cM), BTA6 (68-89 cM) and BTA21 (75-84 cM). The location of QTL for subcutaneous fat thickness on BTA13 was not supported by the second screen (P > 0.05). We determined that the combined contribution of the three QTLs for BMS was 10.1% of the total variance. The combined phenotypic average of these three Q was significantly different (P < 0.001) from those of other allele combinations. Analysis of additional half-sib families will be necessary to confirm these QTL.     

Quantitative trait loci mapping in dairy cattle: review and meta-analysis

From an extensive review of public domain information on dairy cattle quantitative trait loci (QTL), we have prepared a draft online QTL map for dairy production traits. Most publications (45 out of 55 reviewed) reported QTL for the major milk production traits (milk, fat and protein yield, and fat and protein concentration (%)) and somatic cell score. Relatively few QTL studies have been reported for more complex traits such as mastitis, fertility and health. The collated QTL map shows some chromosomal regions with a high density of QTL, as well as a substantial number of QTL at single chromosomal locations. To extract the most information from these published records, a meta-analysis was conducted to obtain consensus on QTL location and allelic substitution effect of these QTL. This required modification and development of statistical methodologies. The meta-analysis indicated a number of consensus regions, the most striking being two distinct regions affecting milk yield on chromosome 6 at 49 cM and 87 cM explaining 4.2 and 3.6 percent of the genetic variance of milk yield, respectively. The first of these regions (near marker BM143) affects five separate milk production traits (protein yield, protein percent, fat yield, fat percent, as well as milk yield).     

Optimal selection on two quantitative trait loci with linkage

A mathematical approach to optimize selection on multiple quantitative trait loci (QTL) and an estimate of residual polygenic effects was applied to selection on two linked or unlinked additive QTL. Strategies to maximize total or cumulative discounted response over ten generations were compared to standard QTL selection on the sum of breeding values for the QTL and an estimated breeding value for polygenes, and to phenotypic selection. Optimal selection resulted in greater response to selection than standard QTL or phenotypic selection. Tight linkage between the QTL (recombination rate 0.05) resulted in a slightly lower response for standard QTL and phenotypic selection but in a greater response for optimal selection. Optimal selection capitalized on linkage by emphasizing selection on favorable haplotypes. When the objective was to maximize total response after ten generations and QTL were unlinked, optimal selection increased QTL frequencies to fixation in a near linear manner. When starting frequencies were equal for the two QTL, equal emphasis was given to each QTL, regardless of the difference in effects of the QTL and regardless of the linkage, but the emphasis given to each of the two QTL was not additive. These results demonstrate the ability of optimal selection to capitalize on information on the complex genetic basis of quantitative traits that is forthcoming.