Detection of quantitative trait loci influencing conformation traits and calving ease in Holstein-Friesian cattle.

An extension of our previous genome scan of a North American Holstein-Friesian population was conducted to identify quantitative trait loci (QTL) affecting conformation traits. Resource families consisted of 1404 sons of 10 elite sires. Genome coverage was estimated to be 2713.5 cM (90%) for 406 markers using a granddaughter design. Regression interval mapping was used to detect QTL affecting 22 conformation traits, including body, udder, feet and legs, and dairy conformation as well as calving ease. Analysis of the families jointly identified 41 chromosome-wise significant QTL influencing conformation traits and 3 significant QTL influencing calving ease on 20 chromosomes. The false discovery rate method was used to account for multiple testing and 3/4 of the suggestive and 5/6 of significant QTL should be real effects. Fourteen of the 44 QTL were significant at the genome-wise level. Comparison of these results with other published reports identifies common QTL affecting conformation traits. Regions on 10 chromosomes appear to affect multiple traits, including conformation, milk production, and somatic cell score, within these particular US Holstein families. Additional work is needed to determine the precise locations of the QTL and select positional candidate genes influencing these traits.

Detection of quantitative trait loci affecting milk production, health, and reproductive traits in Holstein cattle.

We report putative quantitative trait loci affecting female fertility and milk production traits using the merged data from two research groups that conducted independent genome scans in Dairy Bull DNA Repository grandsire families to identify quantitative trait loci (QTL) affecting economically important traits. Six families used by both groups had been genotyped for 367 microsatellite markers covering 2713.5 cM of the cattle genome (90%), with an average spacing of 7.4 cM. Phenotypic traits included PTA for pregnancy rate and daughter deviations for milk, protein and fat yields, protein and fat percentages, somatic cell score, and productive life. Analysis of the merged dataset identified putative quantitative trait loci that were not detected in the separate studies, and the pregnancy rate PTA estimates that recently became available allowed detection of pregnancy rate QTL for the first time. Sixty-one putative significant marker effects were identified within families, and 13 were identified across families. Highly significant effects were found on chromosome 3 affecting fat percentage and protein yield, on chromosome 6 affecting protein and fat percentages, on chromosome 14 affecting fat percentage, on chromosome 18 affecting pregnancy rate, and on chromosome 20 affecting protein percentage. Within-family analysis detected putative QTL associated with pregnancy rate on six chromosomes, with the effect on chromosome 18 being the most significant statistically. These findings may help identify the most useful markers available for QTL detection and, eventually, for marker-assisted selection for improvement of these economically important traits.

Interval and composite interval mapping of somatic cell score, yield, and components of milk in dairy cattle.

Single-marker, interval-mapping (IM) and composite interval mapping (CIM) were used to detect quantitative trait loci (QTL) controlling milk, fat and protein yields, and somatic cell score (SCS). A granddaughter design was used to combine molecular genetic information with predicted transmitting abilities (PTA) and estimated daughter yield deviations (DYD) from eight Dairy Bull DNA Repository Holstein families. Models that included and excluded weights accounting for the uncertainty of the response variable were evaluated in each trait, family and phenotype (DYD and PTA) combination. The genotypic information consisted of 174 microsatellite markers along 29 Bos taurus autosomes. The average number of informative markers per autosome was three and the number of informative sons per family and marker varied between 21 and 173. Within-family results from the least squares single-marker analyses were used in expectation-maximization likelihood IM and CIM implemented with QTL Cartographer. Different CIM model specifications, offering complementary control on the background QTL outside the interval under study, were evaluated. Permutation techniques were used to calculate the genome-wide threshold test statistic values based on 1,000 samples. Results from the DYD and PTA analyses were highly consistent across traits and families. The minor differences in the estimates from the models that accounted for or ignored the uncertainty of the DYD (variance) and PTA (inverse of reliability) may be associated to the elevated and consistent precision of the DYD and PTA among sons. The CIM model best supported by the data had 10 markers controlling for background effects. On autosome (BTA) three, a QTL at 32 cM influencing protein yield was located in family five and a QTL at 74 cM for fat yield was located in family eight. Two map positions associated with SCS were detected on BTA 21, one at 33 cM in family one and the other at 84 cM in family three. A QTL for protein yield was detected between 26 and 36 cM on BTA six, family six, and a QTL for milk yield was detected at 116 cM on BTA seven in family three. The IM and CIM approaches detected a QTL at 3 cM on BTA 14 influencing fat yield in family four. Two map positions on BTA 29 were associated with significant variation of milk (0 cM) and fat yield (14 cM) in family seven. These results suggest the presence of one QTL with pleiotropic effects on multiple traits or multiple QTL within the marker interval. Findings from this study could be used in subsequent fine-mapping work and applied to marker-assisted selection of dairy production and health traits.

Detection of quantitative trait loci influencing dairy traits using a model for longitudinal data.

A longitudinal-linkage analysis approach was developed and applied to an outbred population. Nonlinear mixed-effects models were used to describe the lactation patterns and were extended to include marker information following single-marker and interval mapping models. Quantitative trait loci (QTL) affecting the shape and scale of lactation curves for production and health traits in dairy cattle were mapped in three U.S. Holstein families (Dairy Bull DNA Repository families one, four, and five) using the granddaughter design. Information on 81 informative markers on six Bos taurus autosomes (BTA) was combined with milk yield, fat, and protein percentage and somatic cell score (SCS) test-day records. Six percent of the single-marker tests surpassed the experiment-wise significance threshold. Marker BL41 on BTA3 was associated with decrease in milk yield during mid-lactation in family one. The scale and shape of the protein percentage lactation curve in family four varied with BMC4203 (BTA6) allele that the son received from the grandsire. Some map locations were associated with variation in the lactation pattern of multiple traits. In family four, the marker HUJI177 (BTA3) was associated with changes in the milk yield and protein percentage curves suggesting a QTL with pleiotropic effects or multiple QTL in the region. The interval mapping model uncovered a QTL on BTA7 associated with variation in milk-yield pattern in family four and a QTL on BTA21 affecting SCS in family five. The developed approach can be extended to random regressions, covariance functions, spline, gametic and variance component models. The results from the longitudinal-QTL approach will help to understand the genetic factors acting at different stages of lactation and will assist in positional candidate gene research. Identified positions can be incorporated into marker-assisted selection decisions to alter the persistency and peak production or the fluctuation of SCS during a lactation.

An ordered comparative map of the cattle and human genomes.

A cattle-human whole-genome comparative map was constructed using parallel radiation hybrid (RH) mapping in conjunction with EST sequencing, database mining for unmapped cattle genes, and a predictive bioinformatics approach (COMPASS) for targeting specific homologous regions. A total of 768 genes were placed on the RH map in addition to 319 microsatellites used as anchor markers. Of these, 638 had human orthologs with mapping data, thus permitting construction of an ordered comparative map. The large number of ordered loci revealed > or =105 conserved segments between the two genomes. The comparative map suggests that 41 translocation events, a minimum of 54 internal rearrangements, and repositioning of all but one centromere can account for the observed organizations of the cattle and human genomes. In addition, the COMPASS in silico mapping tool was shown to be 95% accurate in its ability to predict cattle chromosome location from random sequence data, demonstrating this tool to be valuable for efficient targeting of specific regions for detailed mapping. The comparative map generated will be a cornerstone for elucidating mammalian chromosome phylogeny and the identification of genes of agricultural importance."Ought we, for instance, to begin by discussing each separate species-in virtue of some common element of their nature, and proceed from this as a basis for the consideration of them separately?" from Aristotle, On the Parts of Animals, 350 B.C.E.

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.

Exclusion probabilities of 22 bovine microsatellite markers in fluorescent multiplexes for semiautomated parentage testing.

Six multiplexes developed for semiautomated fluorescence genotyping were evaluated for parentage testing. These multiplexes contained primer pairs for the amplification of 22 microsatellites on 17 bovine autosomes. Exclusion probabilities were determined from genotypes of 1022 Holstein cattle and 311 beef cattle belonging to five breeds. Two cases were considered: case 1, genotypes are known for an alleged parent and an offspring but genotypes of a confirmed parent are unknown; and case 2, genotypes are known for an alleged parent, a confirmed parent and an offspring. If the alleged parent is not the true parent, then the 22 markers will exclude the alleged parent with a probability of > 0.9986 for case 1 and with a probability of > 0.99999 for case 2. On the basis of these exclusion probabilities, the probability that an alleged parent will be falsely included as a parent is in the range of 1/716 to 1/2845 for case 1 and 1/1.2 million to 1/159753 for case 2. In addition to these results, a rapid and efficient non-organic method for extraction of DNA from semen is described.

A male linkage map of the cattle (Bos taurus) genome.

A male linkage map of the cattle (Bos taurus) genome was constructed using nine large half-sib families. The map consists of 269 loci, of which 249 are microsatellites and 20 are structural genes. Among the 249 microsatellites, 140 are markers selected from other maps and 98 are new assignments. Chromosome assignment were established for 35 new markers by somatic cell hybrid analysis, of which 26 were confirmed by linkage analysis. Genome coverage is 1975 cM contained within terminal markers on all 29 autosomes. The average distance between adjacent loci is 9.7 cM, with 72.1% of the map intervals < or = 15 cM and 4.9% of the intervals > or = 25 cM. The inclusion of mapped markers permitted integration and comparisons with other maps, facilitating the identification of discrepancies in chromosome assignment, gene order, and map distance. The inclusion of Type I and blood group markers in the map was useful for comparative mapping, revealing possible blood group orthologies between humans and cattle. The map generated will serve as a useful tool for comparative mapping, mapping of quantitative trait loci and marker assisted selection.

Isolation and characterization of 45 polymorphic microsatellites from the bovine genome.

A small-insert bovine genomic library was constructed in pBluescript II SK(+) and enriched for microsatellites by selective rescue of single-stranded pBluescript DNA carrying (CA)n/(TG)n tandem repeats. Approximately 50% of the clones in the enriched library contained (CA)n repeats or CA-rich sequences. Sequencing of clones selected for (CA)n repeats resulted in the identification and characterization of 45 (CA)n polymorphic microsatellites. Genotyping in 9 large paternal half-sib families indicated that 40 of these microsatellite markers exhibit autosomal Mendelian inheritance. The numbers of alleles range from 2 to 18, with an average of 6.3 per locus. The polymorphic microsatellite markers we have identified and characterized will contribute to the construction of a high-resolution linkage map of bovine genome.