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

Genomic Prediction of Complex Traits in Animal Breeding with Long Breeding History, the Dairy Cattle Case

In accordance with the infinitesimal model for quantitative traits, a very large number of genes affect nearly all economic traits. In only two cases has the causative polymorphism been determined for genes affecting economic traits in dairy cattle. Most current methods for genomic evaluation are based on the "two-step" method. Genetic evaluations are computed by the individual animal model, and functions of the evaluations of progeny-tested sires are the dependent variable for estimation of marker effects. With the adoption of genomic evaluation in 2008, annual rates of genetic gain in the US increased from ∼50-100% for yield traits and from threefold to fourfold for lowly heritable traits, including female fertility, herd-life and somatic cell concentration. Gradual elimination of the progeny test scheme has led to a reduction in the number of sires with daughter records and less genetic ties between years. As genotyping costs decrease, the number of cows genotyped will continue to increase, and these records will become the basic data used to compute genomic evaluations, most likely via application of "single-step" methodologies. Less emphasis in selection goals will be placed on milk production traits, and more on health, reproduction, and efficiency traits and "environmentally friendly" production. Genetic variance for economic traits is maintained by increase in frequency of rare alleles, new mutations, and changes in selection goals and management.

A strategy for QTL detection in half-sib populations

A statistical analysis strategy for the detection of quantitative trait loci (QTLs) in half-sib populations is outlined. The initial exploratory analysis is a multiple regression of the trait score on a subset of markers to allow a rapid identification of possible chromosomal regions of interest. This is followed by multiple marker interval mapping with regression methods within and across families fitting one or two QTLs. Empirical thresholds are determined by experiment-wise permutation tests for different significance levels and empirical confidence intervals for the QTLs' positions are obtained by bootstrapping methods. For traits with evidence for a significant single-QTL effect, an approximate maximum likelihood analysis is performed to obtain estimates of QTL effect and the probability of the QTL genotype for each parent of a group of half-sibs. The strategy is demonstrated in an analysis of previously published data on chromosome 6 and five production traits from a granddaughter design in dairy cattle. The results confirm and extend evidence for QTLs affecting protein percentage. Informativeness of markers limited the possibility of mapping more than one QTL on the same linkage group.

Estimation of quantitative trait loci effects in dairy cattle populations

Standard animal model programs can be modified to include the effect of a quantitative gene, even if only a fraction of the population is genotyped. Five methods to estimate the effect of a diallelic quantitative gene affecting a quantitative trait were compared to a standard animal model (model I) on simulated populations, based on mean squared errors and bias. In models II, III, and IV complete linkage between a single genetic marker and the quantitative trait gene was assumed. In models II and III the elements of the incidence matrix for the gene effect were 0 or 1 for genotyped individuals, and the probabilities of the possible candidate gene genotypes for individuals that were not genotyped. In model III segregation analysis was used to compute these probabilities. If only some of the cows were genotyped, the model III estimates were nearly unbiased, while model II underestimated the simulated effects. When only sires were genotyped, model II overestimated the simulated effect. In models V and VI two markers bracketing the quantitative gene with recombination frequencies of 0.1 and 0.2 with the quantitative gene were simulated, and the algorithm of Whittaker et al. (1996) was used to derive estimates of gene effect and location. In model V marker allele effects were included in the animal model analysis. In model VI, the model I genetic evaluations were analyzed. Model V estimates for both effect and location of the quantitative gene were unbiased, while model VI estimates were only 0.25 of the simulated effect.

A whole genome scan for quantitative trait loci affecting milk protein percentage in Israeli-Holstein cattle, by means of selective milk DNA pooling in a daughter design, using an adjusted false discovery rate criterion

Selective DNA pooling was employed in a daughter design to screen all bovine autosomes for quantitative trait loci (QTL) affecting estimated breeding value for milk protein percentage (EBVP%). Milk pools prepared from high and low daughters of each of seven sires were genotyped for 138 dinucleotide microsatellites. Shadow-corrected estimates of sire allele frequencies were compared between high and low pools. An adjusted false discovery rate (FDR) method was employed to calculate experimentwise significance levels and empirical power. Significant associations with milk protein percentage were found for 61 of the markers (adjusted FDR = 0.10; estimated power, 0.68). The significant markers appear to be linked to 19--28 QTL. Mean allele substitution effects of the putative QTL averaged 0.016 (0.009--0.028) in units of the within-sire family standard deviation of EBVP% and summed to 0.460 EBVP%. Overall QTL heterozygosity was 0.40. The identified QTL appear to account for all of the variation in EBVP% in the population. Through use of selective DNA pooling, 4400 pool data points provided the statistical power of 600,000 individual data points.

A genome scan for QTL influencing milk production and health traits in dairy cattle

A genome scan was conducted in the North American Holstein-Friesian population for quantitative trait loci (QTL) affecting production and health traits using the granddaughter design. Resource families consisted of 1,068 sons of eight elite sires. Genome coverage was estimated to be 2,551 cM (85%) for 174 genotyped markers. Each marker was tested for effects on milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score, and productive herd life using analysis of variance. Joint analysis of all families identified marker effects on 11 chromosomes that exceeded the genomewide, suggestive, or nominal significance threshold for QTL effects. Large marker effects on fat percentage were found on chromosomes 3 and 14, and multimarker regression analysis was used to refine the position of these QTL. Half-sibling families from Israeli Holstein dairy herds were used in a daughter design to confirm the presence of the QTL for fat percentage on chromosome 14. The QTL identified in this study may be useful for marker-assisted selection and for selection of a refined set of candidate genes affecting these traits.

A new approach to the problem of multiple comparisons in the genetic dissection of complex traits

Saturated genetic marker maps are being used to map individual genes affecting quantitative traits. Controlling the "experimentwise" type-I error severely lowers power to detect segregating loci. For preliminary genome scans, we propose controlling the "false discovery rate," that is, the expected proportion of true null hypotheses within the class of rejected null hypotheses. Examples are given based on a granddaughter design analysis of dairy cattle and simulated backcross populations. By controlling the false discovery rate, power to detect true effects is not dependent on the number of tests performed. If no detectable genes are segregating, controlling the false discovery rate is equivalent to controlling the experimentwise error rate. If quantitative loci are segregating in the population, statistical power is increased as compared to control of the experimentwise type-I error. The difference between the two criteria increases with the increase in the number of false null hypotheses. The false discovery rate can be controlled at the same level whether the complete genome or only part of it has been analyzed. Additional levels of contrasts, such as multiple traits or pedigrees, can be handled without the necessity of a proportional decrease in the critical test probability.

Detection of putative loci affecting milk production and composition, health, and type traits in a United States Holstein population

Quantitative trait loci affecting milk yield and composition, health, and type traits were studied for seven large grandsire families of US Holstein using the granddaughter design. The families were genotyped at 20 microsatellite markers on 15 chromosomes, and the effects of the marker alleles were analyzed for 28 traits (21 type traits, 5 milk yield and composition traits, somatic cell score, and productive herd life). Markers BM415 on chromosome 6 and BM6425 on chromosome 14 were associated with effects on protein percentage in a single grandsire family. The latter marker had a lower probability of being associated with changes in milk yield and fat percentage in the same family. Increases in productive herd life were associated with an allele at marker BM719 on chromosome 16 in one grandsire family.

Comparative gene mapping: a fine-scale survey of chromosome rearrangements between ruminants and humans

A total of 202 genes were cytogenetically mapped to goat chromosomes, multiplying by five the total number of regional gene localizations in domestic ruminants (255). This map encompasses 249 and 173 common anchor loci regularly spaced along human and murine chromosomes, respectively, which makes it possible to perform a genome-wide comparison between three mammalian orders. Twice as many rearrangements as revealed by ZOO-FISH were observed. The average size of conserved fragments could be estimated at 27 and 8 cM with humans and mice, respectively. The position of evolutionary breakpoints often correspond with human chromosome sites known to be vulnerable to rearrangement in neoplasia. Furthermore, 75 microsatellite markers, 30 of which were isolated from gene-containing bacterial artificial chromosomes (BACs), were added to the previous goat genetic map, achieving 88% genome coverage. Finally, 124 microsatellites were cytogenetically mapped, which made it possible to physically anchor and orient all the linkage groups. We believe that this comprehensive map will speed up positional cloning projects in domestic ruminants and clarify some aspects of mammalian chromosomal evolution.

Mapping quantitative trait loci for milk production and health of dairy cattle in a large outbred pedigree

Quantitative trait loci (QTL) affecting milk production and health of dairy cattle were mapped in a very large Holstein granddaughter design. The analysis included 1794 sons of 14 sires and 206 genetic markers distributed across all 29 autosomes and flanking an estimated 2497 autosomal cM using Kosambi's mapping function. All families were analyzed jointly with least-squares (LS) and variance components (VC) methods. A total of 6 QTL exceeding approximate experiment-wise significance thresholds, 24 QTL exceeding suggestive thresholds, and 34 QTL exceeding chromosome-wise thresholds were identified. Significance thresholds were determined via data permutation (for LS analysis) and chi-square distribution (for VC analysis). The average bootstrap confidence interval for the experiment-wise significant QTL was 48 cM. Some chromosomes harbored QTL affecting several traits, and these were always in coupling phase, defined by consistency with genetic correlations among traits. Chromosome 17 likely harbors 2 QTL affecting milk yield, and some other chromosomes showed some evidence for 2 linked QTL affecting the same trait. In each of these cases, the 2 QTL were in repulsion phase in those families appearing to be heterozygous for both QTL, a finding which supports the build-up of linkage disequilibrium due to selection.

Estimation of candidate gene effects in dairy cattle populations

Standard programs for animal models can be modified to include the effect of a candidate gene, even if only a fraction of the population is genotyped. The elements of the incidence matrix for this effect is 0 or 1 for genotyped individuals and for the probabilities of the candidate gene genotypes for individuals that were not genotyped. The effects of a diallelic candidate gene that were estimated by this method were compared on simulated populations with three alternative estimation methods: analysis of genetic evaluations, yield deviations, and daughter yield deviations, all of which were derived from a standard animal model. The bases of comparison were mean squared error and bias. Four types of experimental designs were considered: genotyping sires only, genotyping cows only, genotyping half of the cows but no sires, and genotyping half of the sires and half of the cows. The estimates that were derived from the three alternative methods all underestimated the simulated effects. The genetic evaluations were more biased than were the yield deviations and the daughter yield deviations. The proposed method was significantly biased only for the design in which half of the daughters, but no sires, were genotyped. Bias and mean squared errors were always lowest by the proposed method.

High-precision genotyping by denaturing capillary electrophoresis

Genotyping, as applied to linkage mapping, human identification, or mapping of genetic traits, mandates electrophoretic separation systems that enable a user to identify alleles with high precision to obtain a correct genotype. For 2-bp microsatellites or short tandem repeats (STRs), standard deviations of +/-0.3 nucleotide are required to ensure with 99.7% probability the identity or dissimilarity of tested alleles. A complete system, consisting of commercially available laser-induced fluorescence capillary electrophoresis (ABI PRISM 310) and performance optimized polymer 4 (POP-4), was evaluated for microsatellite separations. POP-4 is a low viscosity polymer for use in uncoated fused microbore silica capillaries. It separates DNA fragments that differ in size by 1 nucleotide up to 250 nucleotides and that differ in size by 2 nucleotides for fragments up to at least 350 nucleotides in length in about 30 min. The presence of denaturants and, more importantly, operation at 60 degrees C was mandatory for high-precision and high-resolution sizing operation. Reproducible separation performance was achieved in excess of 100 injections per capillary with resulting standard deviations in the range of 0.04 to 0.17 nucleotide. Comparative sizing of known CEPH (Centre d'Etudes du Polymorphisme Humaine) samples performed at 22 independent test sites showed the usefulness of the system for genotyping with standard deviations of 0.24 nucleotide, or better.

Multiple marker mapping of quantitative trait loci of Finnish dairy cattle by regression

A total of 453 bulls belonging to 11 half-sib families of Finnish Ayrshires were genotyped for six microsatellite markers on chromosome 9. The data were used in an attempt to map quantitative trait loci applying regression as a multimarker approach. For association analysis with a granddaughter design, the EBV for 12 traits were used: milk yield, protein yield, fat percentage, protein percentage, daughter weight, bull growth, calf mortality, days open, fertility treatments, nonreturn rate, SCC, and clinical mastitis. The empirical values of significance thresholds were determined using a permutation test on the experimental data. Although no significant effects were found, the results indicate some support for the existence of a locus on chromosome 9 that affects milk and protein yields.

Multiple genotype analysis and sexing of IVF bovine embryos

Twenty-one in vitro-fertilized bovine blastocysts were quartered, lysed and subjected to primer elongation preamplification (PEP) procedure, allowing for the analysis of up to 40 genotypes per quarter embryo. The quarter-embryos were sexed by polymerase chain reaction (PCR) using BRY.1, Bov97M and ZFX/ZFY loci, and then genotyped for k-casein, bovine leukocyte adhesion deficiency (BLAD) and microsatellite D9S1. The mitochondrial cytochrome B locus was used as an internal control with a 95% success rate. The PEP procedure amplified genomic fragments in 93% of all cases. The embryos were identified to be 11 males and 10 females. Sexing accuracy was 87% for BRY.1, 97% for ZFX/ZFY and 100% for Bov97M. False genotyping was due mostly to amplification of BRY.1 in the female embryos and to the nonamplification of the ZFY locus in the male embryos. The results indicate that the combined use of Bov97M and ZFX/ZFY loci is a highly accurate procedure for sexing bovine embryos. Genotyping for kappa-casein, D9S1 and BLAD was successful in 94, 99 and 91% of assays, respectively. Sex ratios and allele frequencies of embryos for gk-casein, BLAD and D9S1 were all close to the observed frequencies in the Israeli Holstein population. These results support the conclusion that the genotyping of embryos is as accurate as that of mature animals. Thus, marker-assisted selection can be efficiently applied at the preimplantation embryo level for loci of economic importance.

Marker-assisted selection for a sex-limited character in a nucleus breeding population

The benefits of marker-assisted selection were examined by simulation of an adult multiple ovulation and embryo transfer nucleus breeding scheme. Animals were either typed for two polymorphic marker loci, 20 centimorgans apart, flanking a single biallelic quantitative trait locus and were evaluated using a model accounting for marker information, or animals were not typed but were evaluated by a conventional BLUP animal model. Selection was for a single trait measured on females, and each dam had 4 sons and 4 daughters. Nucleus foundation animals were chosen from a base population in linkage equilibrium. With the favorable allele at an initial frequency of 0.5, marker-assisted selection substantially increased responses at the quantitative trait locus but reduced the polygenic responses. Cumulative genetic gain increased by up to 3, 9, 12, and 6% after one, two, three, and six generations of selection, respectively. If the favorable allele was initially rare (frequency of 0.1), the merits of marker-assisted selection were even more pronounced (genetic gains increased by up to 9, 19, 24 and 15%, respectively). The superiority of marker-assisted selection over conventional BLUP increased when a restriction was placed on selection of full brothers and decreased when variance of the quantitative trait locus used in the evaluation model was overestimated.

Misidentification rate in the Israeli dairy cattle population and its implications for genetic improvement

The DNA microsatellites can be efficiently used to determine incorrect paternity attribution of cattle without genotyping of dams. Allelic frequencies of the population were determined for 12 microsatellites using the maternal alleles of 102 AI sires. The frequency of the most common microsatellite allele ranged from 0.27 to 0.58. Most loci had at least one allele that was present in only a single individual. Paternity of 9 of 173 cows (5.2%) and 3 of 102 bulls (2.9%) was excluded because putative paternal alleles were not present in progeny for at least one locus. For 4 of the 9 cows and all 3 bulls, exclusion was based on at least two loci. Mean probability of exclusion was 0.85 for cows and 0.99 for bulls. With an assumed cost of US $5 per genotype, a misidentification rate of 5%, and a discount rate of 0.05, additional profit for the Israeli-Holstein breeding program from genotyping 100 test daughters of each young sire becomes positive within 10 yr and reaches nearly US $2.4 million after 20 yr.

Prediction of informativeness for microsatellite markers among progeny of sires used for detection of economic trait loci in dairy cattle

Individual loci affecting economic traits can be located using genetic linkage. Application of either daughter or granddaughter design requires determination of allele origin in the progeny. If only the sires and their progeny are genotyped, the paternal allele origin of progeny having the same genotype as the sire cannot be determined. The expected frequency of informative sons can be predicted for each sire and genetic marker from the allele frequencies in the population. The accuracy of a predictor of the frequency of informative progeny was tested on 103 grandsire x microsatellite combinations. Number of sons per grandsire varied from 24 to 129. Allele frequencies in the population were estimated by genotyping seven sires. The regression of the frequency of informative sons on the predicted frequency was 1.04 with a zero intercept model. Thus, considering the large number of genetic markers available for analysis, predicted informative frequency is a useful criterion for selection of genetic markers.

Prospects for the genetic manipulation of dairy cattle: opportunities beyond BST

The dairy industry, with regulatory approvals of recombinant chymosin and bovine somatotropin (BST), has been at the forefront of food and agricultural biotechnology. The commercial fate of these products is one of several factors that may affect the success of future genetic manipulations in dairy cattle and dairy products. Other factors include technical and reproductive constraints in cattle and the cost of producing transgenic cattle. Early applications of genetic manipulation in cattle, for reasons of cost recoupment, may favor production of heterologous proteins in milk for pharmaceutical and medical use. Such applications could benefit genetic modification of milk and milk proteins for food use by providing advance knowledge and experience in mammalian protein expression. Other research opportunity areas that could affect prospects for genetic manipulation of dairy cattle include genome mapping, metabolic pathways, growth and development, and cattle/microbe interactions.