A QTL affecting milk yield and composition maps to bovine chromosome 20: a confirmation

The prolactin receptor: A cross-species comparison of gene structure, transcriptional regulation, tissue-specificity, and genetic variation

The conserved and multifaceted functions of prolactin (PRL) are coordinated through varied distribution and expression of its cell-surface receptor (PRLR) across a range of tissues and physiological states. The resultant heterogeneous expression of PRLR mRNA and protein across different organs and cell types supports a wide range of PRL-regulated processes including reproduction, lactation, development, and homeostasis. Genetic variation within the PRLR gene also accounts for several phenotypes impacting agricultural production and human pathology. The goal of this review is to highlight the many elements that control differential expression of the PRLR across tissues, and the various phenotypes that exist across species due to variation in the PRLR gene.

Genome-wide association study for buffalo mammary gland morphology

This research communication describes a genome-wide association study for Italian buffalo mammary gland morphology. Three single nucleotide polymorphisms (AX-85117983, AX-8509475 and AX-85117518) were identified to be significantly associated with buffalo anterior teat length, posterior teat length and distance between anterior and posterior teat, respectively. Two significant signals for buffalo mammary gland morphology were observed in two genomic regions on the chromosome 10, and chromosome 20. One of the regions located on the chromosome 10 has the most likely candidate genes ACTC1 and GJD2, both of which have putative roles in the regulation of mammary gland development. This study provides new insights into the genetic variants of buffalo mammary gland morphology and may be beneficial for understanding of the genetic regulation of mammary growth.

Identification of single nucleotide polymorphisms of PIK3R1 and DUSP1 genes and their genetic associations with milk production traits in dairy cows

Background

Previously, phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1) and dual specificity phosphatase 1 (DUSP1) were identified as promising candidate genes for milk production traits due to their being differentially expressed between the dry period and the peak of lactation in livers of dairy cows. Hence, in this study, the single nucleotide polymorphisms (SNPs) of PIK3R1 and DUSP1 genes were identified and their genetic associations with milk yield, fat yield, fat percentage, protein yield, and protein percentage, were investigated using 1067 Chinese Holstein cows from 40 sire families.

Results

By re-sequencing the entire coding region and 2000 bp of the 5′ and 3′ flanking regions of the two genes, one SNP in the 5′ untranslated region (UTR), three in the 3′ UTR, and two in the 3′ flanking region of PIK3R1 were identified, and one in the 5′ flanking region, one in the 3′ UTR, and two in the 3′ flanking region of DUSP1 were found. Subsequent single-locus association analyses showed that five SNPs in PIK3R1, rs42590258, rs210389799, rs208819656, rs41255622, rs133655926, and rs211408208, and four SNPs in DUSP1, rs207593520, rs208460068, rs209154772, and rs210000760, were significantly associated with milk, fat and protein yields in the first or second lactation (P values ≤ 0.0001 and 0.0461). In addition, by the Haploview 4.2 software, the six and four SNPs in PIK3R1 and DUSP1 respectively formed one haplotype block, and the haplotype-based association analyses showed significant associations between their haplotype combinations and the milk traits in both two lactations (P values ≤ 0.0001 and 0.0364). One SNP, rs207593520(T/G), was predicted to alter the transcription factor binding sites (TFBSs) in the 5′ flanking region of DUSP1. Further, the dual-luciferase assay showed that the transcription activity of allele T in rs207593520 was significantly higher than that of allele G, suggesting the activation of transcriptional activity of DUSP1 gene by allele T of rs207593520. Thus, the rs207593520 SNP was highlighted as a potential causal mutation that should be further verified.

Conclusions

We demonstrated novel and significant genetic effects of the PIK3R1 and DUSP1 genes on milk production traits in dairy cows, and our findings provide information for use in dairy cattle breeding.

Genome-wide association studies for genetic effects that change during lactation in dairy cattle

Genetic effects on milk production traits in dairy cattle might change during lactation. However, most genome-wide association studies (GWAS) for milk production traits assume that genetic effects are constant during lactation. This assumption might lead to missing these quantitative trait loci (QTL) whose effects change during lactation. This study aimed to screen the whole genome specifically for QTL whose effects change during lactation. For this purpose, 4 different GWAS approaches were performed using test-day milk protein content records: (1) separate GWAS for specific lactation stages, (2) GWAS for estimated Wilmink lactation curve parameters, (3) a GWAS using a repeatability model where SNP effects are assumed constant during lactation, and (4) a GWAS for genotype by lactation stage interaction using a repeatability model and accounting for changing genetic effects during lactation. Separate GWAS for specific lactation stages suggested that the detection power greatly differs between lactation stages and that genetic effects of some QTL change during lactation. The GWAS for estimated Wilmink lactation curve parameters detected many chromosomal regions for Wilmink parameter a (protein content level), whereas 2 regions for Wilmink parameter b (decrease in protein content toward nadir) and no regions for Wilmink parameter c (increase in protein content after nadir) were detected. Twenty chromosomal regions were detected with effects on milk protein content; however, there was no evidence that their effects changed during lactation. For 5 chromosomal regions located on chromosomes 3, 9, 10, 14, and 27, significant evidence was observed for a genotype by lactation stage interaction and thus their effects on milk protein content changed during lactation. Three of these 5 regions were only identified using a GWAS for genotype by lactation stage interaction. Our study demonstrated that GWAS for genotype by lactation stage interaction offers new possibilities to identify QTL involved in milk protein content. The performed approaches can be applied to other milk production traits. Identification of QTL whose genetic effects change during lactation will help elucidate the genetic and biological background of milk production.

Association of candidate gene polymorphisms with milk technological traits, yield, composition, and somatic cell score in Italian Holstein-Friesian sires

Advances in DNA-based marker technology have enabled the identification of genomic regions underlying complex phenotypic traits in livestock species. The incorporation of detected quantitative trait loci into genetic evaluation provides great potential to enhance selection accuracies, hence expediting the genetic improvement of economically important traits. The objective of the present study was to investigate 96 single nucleotide polymorphisms (SNP) located in 53 candidate genes previously reported to have effects on milk production and quality traits in a population of highly selected Holstein-Friesian bulls. A total of 423 semen samples were used to genotype the bulls through a custom oligo pool assay. Forty-five SNP in 32 genes were found to be associated with at least 1 of the tested traits. Most significant and favorable SNP trait associations were observed for polymorphisms located in CCL3 and AGPAT6 genes for fat yield (0.037 and 0.033 kg/d, respectively), DGKG gene for milk yield (0.698 kg/d), PPARGC1A, CSN1S1, and AGPAT6 genes for fat percentage (0.127, 0.113, and 0.093%, respectively), GHR gene for protein (0.064%) and casein percentage (0.053%), and TLR4 gene for fat (0.090%), protein (0.066%), and casein percentage (0.050%). Somatic cell score was favorably affected by GHR (-0.095) and POU1F1 (-0.137), and interesting SNP-trait associations were observed for polymorphisms located in CSN2, POU1F1, and AGPAT6 genes for rennet coagulation time (-0.592, -0.558, and -0.462 min, respectively), and GHR and CSN2 genes for curd firmness 30 min after rennet addition (1.264 and 1.183 mm, respectively). In addition to the influence of individual SNP, the effects of composite genotypes constructed by grouping SNP according to their individual effects on traits considered in the analysis were also examined. Favorable and significant effects on milk traits were observed for 2 composite genotypes, one including 10 SNP and the other 4 SNP. The former was associated with an increase of milk (0.075 kg/d), fat (0.097 kg/d), protein (0.083 kg/d), and casein yields (0.065 kg/d), and the latter was associated with an increase of fat (0.244%), protein (0.071%), and casein percentage (0.047%). Although further research is required to validate the identified SNP loci in other populations and breeds, our results can be considered as a preliminary foundation for further replication studies on gene-assisted selection programs.

Transcriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing

Background

Recently, RNA sequencing (RNA-seq) has rapidly emerged as a major transcriptome profiling system. Elucidation of the bovine mammary gland transcriptome by RNA-seq is essential for identifying candidate genes that contribute to milk composition traits in dairy cattle.

Results

We used massive, parallel, high-throughput, RNA-seq to generate the bovine transcriptome from the mammary glands of four lactating Holstein cows with extremely high and low phenotypic values of milk protein and fat percentage. In total, we obtained 48,967,376–75,572,578 uniquely mapped reads that covered 82.25% of the current annotated transcripts, which represented 15549 mRNA transcripts, across all the four mammary gland samples. Among them, 31 differentially expressed genes (p < 0.05, false discovery rate q < 0.05) between the high and low groups of cows were revealed. Gene ontology and pathway analysis demonstrated that the 31 differently expressed genes were enriched in specific biological processes with regard to protein metabolism, fat metabolism, and mammary gland development (p < 0.05). Integrated analysis of differential gene expression, previously reported quantitative trait loci, and genome-wide association studies indicated that TRIB3, SAA (SAA1, SAA3, and M-SAA3.2), VEGFA, PTHLH, and RPL23A were the most promising candidate genes affecting milk protein and fat percentage.

Conclusions

This study investigated the complexity of the mammary gland transcriptome in dairy cattle using RNA-seq. Integrated analysis of differential gene expression and the reported quantitative trait loci and genome-wide association study data permitted the identification of candidate key genes for milk composition traits.

Genome-wide associations for milk production and somatic cell score in Holstein-Friesian cattle in Ireland

BACKGROUND

Contemporary dairy breeding goals have broadened to include, along with milk production traits, a number of non-production-related traits in an effort to improve the overall functionality of the dairy cow. Increased indirect selection for resistance to mastitis, one of the most important production-related diseases in the dairy sector, via selection for reduced somatic cell count has been part of these broadened goals. A number of genome-wide association studies have identified genetic variants associated with milk production traits and mastitis resistance, however the majority of these studies have been based on animals which were predominantly kept in confinement and fed a concentrate-based diet (i.e. high-input production systems). This genome-wide association study aims to detect associations using genotypic and phenotypic data from Irish Holstein-Friesian cattle fed predominantly grazed grass in a pasture-based production system (low-input).

RESULTS

Significant associations were detected for milk yield, fat yield, protein yield, fat percentage, protein percentage and somatic cell score using separate single-locus, frequentist and multi-locus, Bayesian approaches. These associations were detected using two separate populations of Holstein-Friesian sires and cows. In total, 1,529 and 37 associations were detected in the sires using a single SNP regression and a Bayesian method, respectively. There were 103 associations in common between the sires and cows across all the traits. As well as detecting associations within known QTL regions, a number of novel associations were detected; the most notable of these was a region of chromosome 13 associated with milk yield in the population of Holstein-Friesian sires.

CONCLUSIONS

A total of 276 of novel SNPs were detected in the sires using a single SNP regression approach. Although obvious candidate genes may not be initially forthcoming, this study provides a preliminary framework upon which to identify the causal mechanisms underlying the various milk production traits and somatic cell score. Consequently this will deepen our understanding of how these traits are expressed.

A microsatellite-based analysis for the detection of selection on BTA1 and BTA20 in northern Eurasian cattle (Bos taurus) populations

BACKGROUND

Microsatellites surrounding functionally important candidate genes or quantitative trait loci have received attention as proxy measures of polymorphism level at the candidate loci themselves. In cattle, selection for economically important traits is a long-term strategy and it has been reported that microsatellites are linked to these important loci.

METHODS

We have investigated the variation of seven microsatellites on BTA1 (Bos taurus autosome 1) and 16 on BTA20, using bovine populations of typical production types and horn status in northern Eurasia. Genetic variability of these loci and linkage disequilibrium among these loci were compared with those of 28 microsatellites on other bovine chromosomes. Four different tests were applied to detect molecular signatures of selection.

RESULTS

No marked difference in locus variability was found between microsatellites on BTA1, BTA20 and the other chromosomes in terms of different diversity indices. Average D' values of pairwise syntenic markers (0.32 and 0.28 across BTA 1 and BTA20 respectively) were significantly (P < 0.05) higher than for non-syntenic markers (0.15). The Ewens-Watterson test, the Beaumont and Nichol's modified frequentist test and the Bayesian FST-test indicated elevated or decreased genetic differentiation, at SOD1 and AGLA17 markers respectively, deviating significantly (P < 0.05) from neutral expectations. Furthermore, lnRV, lnRH and lnRtheta' statistics were used for the pairwise population comparison tests and were significantly less variable in one population relative to the other, providing additional evidence of selection signatures for two of the 51 loci. Moreover, the three Finnish native populations showed evidence of subpopulation divergence at SOD1 and AGLA17. Our data also indicate significant intergenic linkage disequilibrium around the candidate loci and suggest that hitchhiking selection has played a role in shaping the pattern of observed linkage disequilibrium.

CONCLUSION

Hitchhiking due to tight linkage with alleles at candidate genes, e.g. the POLL gene, is a possible explanation for this pattern. The potential impact of selective breeding by man on cattle populations is discussed in the context of selection effects. Our results also suggest that a practical approach to detect loci under selection is to simultaneously apply multiple neutrality tests based on different assumptions and estimations.

Promoter region of the bovine growth hormone receptor gene: single nucleotide polymorphism discovery in cattle and association with performance in Brangus bulls

Expression of the GH receptor (GHR) gene and its binding with GH is essential for growth and fat metabolism. A GT microsatellite exists in the promoter of bovine GHR segregating short (11 bp) and long (16 to 20 bp) allele sequences. To detect SNP and complete an association study of genotype to phenotype, we resequenced a 1,195-bp fragment of DNA including the GT microsatellite and exon 1A. Resequencing was completed in 48 familialy unrelated Holstein, Jersey, Brown Swiss, Simmental, Angus, Brahman, and Brangus cattle. Nine SNP were identified. Phylogeny analyses revealed minor distance (i.e., <5%) in DNA sequence among the 5 Bos taurus breeds; however, sequence from Brahman cattle averaged 27.4 +/- 0.07% divergence from the Bos taurus breeds, whereas divergence of Brangus was intermediate. An association study of genotype to phenotype was completed with data from growing Brangus bulls (n = 553 from 96 sires) and data from 4 of the SNP flanking the GT microsatellite. These SNP were found to be in Hardy-Weinberg equilibrium and in phase based on linkage disequilibrium analyses (r(2) = 0.84 and D'= 0.92). An A/G tag SNP was identified (ss86273136) and was located in exon 1A, which began 88 bp downstream from the GT microsatellite. Minor allele frequency of the tag SNP was greater than 10%, and Mendelian segregation was verified in 3 generation pedigrees. The A allele was derived from Brahman, and the G allele was derived from Angus. This tag SNP genotype was a significant effect in analyses of rib fat data collected with ultrasound when bulls were ~365 d of age. Specifically, bulls of the GG genotype had 6.1% more (P = 0.0204) rib fat than bulls of the AA and AG genotypes, respectively. Tag SNP (ss86273136), located in the promoter of GHR, appears to be associated with a measure of corporal fat in Bos taurus x Bos indicus composite cattle.

Identification of bovine QTL for growth and carcass traits in Japanese Black cattle by replication and identical-by-descent mapping

To map quantitative trait loci (QTL) for growth and carcass traits in a purebred Japanese Black cattle population, we conducted multiple QTL analyses using 15 paternal half-sib families comprising 7860 offspring. We identified 40 QTL with significant linkages at false discovery rates of less than 0.1, which included 12 for intramuscular fat deposition called marbling and 12 for cold carcass weight or body weight. The QTL each explained 2%-13% of the phenotypic variance. These QTL included many replications and shared hypothetical identical-by-descent (IBD) alleles. The QTL for CW on BTA14 was replicated in five families with significant linkages and in two families with a 1% chromosome-wise significance level. The seven sires shared a 1.1-Mb superior Q haplotype as a hypothetical IBD allele that corresponds to the critical region previously refined by linkage disequilibrium mapping. The QTL for marbling on BTA4 was replicated in two families with significant linkages. The QTL for marbling on BTA6, 7, 9, 10, 20, and 21 and the QTL for body weight on BTA6 were replicated with 1% and/or 5% chromosome-wise significance levels. There were shared IBD Q or q haplotypes in the marbling QTL on BTA4, 6, and 10. The allele substitution effect of these haplotypes ranged from 0.7 to 1.2, and an additive effect between the marbling QTL on BTA6 and 10 was observed in the family examined. The abundant and replicated QTL information will enhance the opportunities for positional cloning of causative genes for the quantitative traits and efficient breeding using marker-assisted selection.

The role of the bovine growth hormone receptor and prolactin receptor genes in milk, fat and protein production in Finnish Ayrshire dairy cattle

We herein report new evidence that the QTL effect on chromosome 20 in Finnish Ayrshire can be explained by variation in two distinct genes, growth hormone receptor (GHR) and prolactin receptor (PRLR). In a previous study in Holstein-Friesian dairy cattle an F279Y polymorphism in the transmembrane domain of GHR was found to be associated with an effect on milk yield and composition. The result of our multimarker regression analysis suggests that in Finnish Ayrshire two QTL segregate on the chromosomal region including GHR and PRLR. By sequencing the coding sequences of GHR and PRLR and the sequence of three GHR promoters from the pooled samples of individuals of known QTL genotype, we identified two substitutions that were associated with milk production traits: the previously reported F-to-Y substitution in the transmembrane domain of GHR and an S-to-N substitution in the signal peptide of PRLR. The results provide strong evidence that the effect of PRLR S18N polymorphism is distinct from the GHR F279Y effect. In particular, the GHR F279Y has the highest influence on protein percentage and fat percentage while PRLR S18N markedly influences protein and fat yield. Furthermore, an interaction between the two loci is suggested.

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.

Characterization and RH mapping of bovine microsatellites generated from a microdissected BTA20-specific DNA library

Bovine chromosome 20 (BTA20) is associated with several quantitative trait loci (QTL) for meat tenderness, birth weight, milk yield and composition. Fine mapping of these QTL requires the development of additional informative markers to increase the resolution of the BTA20 genetic and physical maps. A BTA20-specific library was constructed by means of microdissection and microcloning, and screened for dinucleotide repeats with (CA)16 and (GT)16 oligos. A total of 60 new microsatellites (MS) were developed and characterized for polymorphism using the U.S. Department of Agriculture (USDA)/Meat Animal Research Center (MARC) bovine reference family, of which 53 markers were informative in this family. The number of alleles for these loci varied from 1 to 14, with an average of 6.5. Thirty-three of these MSs, together with 105 markers previously mapped to BTA20, were scored on a 7000-rad cattle-hamster whole-genome radiation hybrid panel (SUNbRH), resulting in a high-resolution RH7000 rad map for BTA20.

Detection of different quantitative trait loci for antibody responses to keyhole lympet hemocyanin and Mycobacterium butyricum in two unrelated populations of laying hens

Quantitative trait loci involved in the primary antibody response to keyhole lympet hemocyanin (KLH) and Mycobacterium butyricum were detected in two independent populations of laying hens. The first population was an F2 cross (H/L) of lines divergently selected for either high or low primary antibody responses to SRBC, and the second population was an F2 cross between 2 commercial layer lines displaying differences in feather pecking behavior (FP). Both populations were typed with microsatellite markers widely distributed over the genome with similar intervals between markers. Titers of antibodies binding KLH and M. butyricum were measured for all individuals by ELISA. Two genetic models were applied to detect QTL involved in the humoral immune response: a half-sib model and a line-cross model, both using the regression interval method. In the half-sib analysis, 2 QTL (on GGA14 and GGA27) were detected for the antibody response to KLH for the H/L population, and 2 QTL (on GGA14 and GGA18) were detected for the FP population. Only 1 QTL was detected for M. butyricum on GGA14 in the FP population using the half-sib analysis model. Two QTL were detected for the FP population on GGA2 and GGA3 using the line-cross analysis model. A QTL for the primary antibody response to KLH detected on GGA14 was validated in both populations under the half-sib analysis model. The present data suggest differences in the genetic regulation of antibody responses to two different T-cell dependent antigens.

Quantitative Trait Loci for Health Traits in Finnish Ayrshire Cattle

A whole-genome scan was conducted to search for quantitative trait loci (QTL) affecting health traits in Finnish Ayrshire dairy cattle. The mapping population consisted of 12 bulls and their 491 sons in a granddaughter design. A total of 150 markers were typed covering all 29 autosomes. The traits under study were somatic cell score, mastitis, and a group of other veterinary treatments. Effects of the QTL and positions were estimated with the regression method. When carrying out interval mapping on each chromosome, cofactors were used to adjust for QTL identified at other chromosomes. Empirical P-values were obtained by permutation. Altogether 17 QTL were detected with genomewise significant P-values in the across family analysis. Quantitative trait loci affecting SCS were identified on chromosomes 1, 3, 11, 18, 21, 24, 27, 29, and QTL for mastitis on chromosomes 14, 18. Quantitative trait loci for other veterinary treatments were found on chromosomes 1, 2, 5, 8, 15, 22, and 23. The allele substitution effects were from 0.5 to 1.7 genetic standard deviations. The positions of these health QTL did not overlap with milk QTL detected in previous studies of the same population.

Genetic Manipulation of Mammary Gland Development and Lactation

The mammalian genome is believed to contain some 30,000 to 40,000 different genes. Of these an estimated 42% have no known function. Genetically engineered mouse models (GEMM) have been a powerful tool available for determining gene function in vivo. In the mammary gland, a variety of genetic engineering approaches have been applied successfully to understanding the importance of specific gene products to mammary gland development and lactation. Our own laboratory has applied genetically engineered mice to facilitate understanding of the regulation of mammary gland development and lactation by insulin-like growth factors (IGF) and by the transcription factor, upstream stimulatory factor (USF-2). Our studies on transgenic mice that overexpress IGF-I have demonstrated the importance of IGF-dependent signaling pathways to maintenance of mammary epithelial cells during the declining phase of lactation. Our analysis of early developmental processes in mammary tissue from mice that carry a targeted mutation in the IGF-I receptor gene suggests that IGF-dependent stimulation of cell cycle progression is more important to early mammary gland development than potential antiapoptotic effects. Lastly, our studies on mice that carry a targeted mutation of the Usf2 gene have demonstrated that this gene is necessary for normal lactation and have highlighted the importance of this gene to the maintenance of protein synthesis. These studies, as well as studies of others, have highlighted both the strengths and limitations inherent in the use of GEMM. Limitations serve as the driving force behind development of new experimental strategies and genetic engineering schemes that will allow for a full understanding of gene function within the mammary gland.

Molecular dissection of a quantitative trait locus: a phenylalanine-to-tyrosine substitution in the transmembrane domain of the bovine growth hormone receptor is associated with a major effect on milk yield and composition

We herein report on our efforts to improve the mapping resolution of a QTL with major effect on milk yield and composition that was previously mapped to bovine chromosome 20. By using a denser chromosome 20 marker map and by exploiting linkage disequilibrium using two distinct approaches, we provide strong evidence that a chromosome segment including the gene coding for the growth hormone receptor accounts for at least part of the chromosome 20 QTL effect. By sequencing individuals with known QTL genotype, we identify an F to Y substitution in the transmembrane domain of the growth hormone receptor gene that is associated with a strong effect on milk yield and composition in the general population.

A genome scan for quantitative trait loci affecting milk production in Norwegian dairy cattle

An autosomal genome scan for quantitative trait loci (QTL) affecting milk production traits was carried out on the Norwegian Dairy Cattle population. Six half-sibling families with a total of 285 sons organized according to a granddaughter design were analyzed by a multiple marker regression method. Suggestive QTL for one or several of the five milk traits (milk yield, protein percentage, protein yield, fat percentage and fat yield) were detected on chromosomes 3, 5, 6, 11, 13, 18 and 20. Among these results, the findings on chromosomes 3, 6, and 20 are highly supported by literature. The most convincing result was found close to marker FBN9 on chromosome 6, where a QTL was detected with alleles that cause a marked reduction in both protein and fat percentages and an increase in milk yield. The results for fat and protein percentage were highly significant even after accounting for multiple testing across the genome. Using bootstrapping, a 95% confidence interval for the position of the QTL for the percentage traits on chromosome 6 was estimated to 16 cM.

Genome-wide search for markers associated with bovine spongiform encephalopathy

A genome-wide search for markers associated with BSE incidence was performed by using Transmission-Disequilibrium Tests (TDTs). Significant segregation distortion, i.e., unequal transmission probabilities of alleles within a locus, was found for three marker loci on Chromosomes (Chrs) 5, 10, and 20. Although TDTs are robust to false associations owing to hidden population substructures, it cannot distinguish segregation distortion caused by a true association between a marker and bovine spongiform encephalopathy (BSE) from a population-wide distortion. An interaction test and a segregation distortion analysis in half-sib controls were used to disentangle these two alternative hypotheses. None of the markers showed any significant interaction between allele transmission rates and disease status, and only the marker on Chr 10 showed a significant segregation distortion in control individuals. Nevertheless, the control group may have been a mixture of resistant and susceptible but unchallenged individuals. When new genotypes were generated in the vicinity of these three markers, evidence for an association with BSE was confirmed for the locus on Chr 5.

A genome scan to identify quantitative trait loci affecting economically important traits in a US Holstein population

Quantitative trait loci affecting economically important traits were studied for eight large US Holstein grandsire families by using the granddaughter design. A total of 155 microsatellite markers located throughout the bovine genome were selected for the scan. The data analyzed include genotypes for 50 markers not previously reported. Results analyses of 105 marker genotypes reported previously were updated. Effects of marker alleles were analyzed for 38 traits including traits for milk production, somatic cell score, productive life, conformation, calving ease, and 16 canonical traits derived from conformation and production traits. Permutation tests were used to calculate empirical traitwise error rates. A traitwise critical value of P = 0.1 was used to determine significance. Ten putative quantitative trait loci associated with seven of the new markers were identified within specific families. One marker on chromosome 14 was associated with differences in fat yield, fat percentage, and a canonical production trait in two families. Markers on chromosomes 18 and 22 were associated with differences in rump angle in the same family. Markers were associated with differences in udder depth and fore udder attachment on chromosomes 16 and 20, respectively. One marker on chromosome 27 was associated with a difference in the dairy capacity composite index, and another marker on chromosome 13 was associated with a difference in a canonical conformation trait. These additional markers complete our genome scan to identify quantitative trait loci affecting economically important traits in a selected commercial Holstein population. The quantitative trait loci identified in this genome scan may be useful for marker-assisted selection to increase the rate of genetic improvement on traits such as disease resistance and conformation traits associated with fitness while accelerating genetic improvement for production.

Detection of quantitative trait loci affecting milk production traits on 10 chromosomes in Holstein cattle

Sons (n = 71 to 75) of each of six Holstein sires were genotyped at 69 microsatellite loci covering a total of 676 cM on chromosomes 3, 5, 9, 10, 13, 15, 17, 20, 23, and 26. Estimates of quantitative trait loci (QTL) effect and location were made using a least squares interval mapping approach based on daughter yield deviations of sons for 305 d milk, fat, and protein yield and fat and protein percentage. Thresholds for statistical significance of QTL effects were determined from interval mapping of 10,000 random permutations of the data across the bull sire families and within each sire family separately. Analyses combining data across sires indicated the presence of QTL affecting milk, fat, and protein yield on chromosomes 20 and 26 and a QTL affecting fat and protein percentage on chromosome 3. Analyses within each sire family separately indicated the presence of segregating QTL in at least one family on 7 of the 10 chromosomes. Statistically significant estimates of QTL effects on breeding value ranged from 438 to 658 kg of milk, from 17.4 to 24.9 kg of fat, 13.0 to 17.0 kg of protein, 0.04 to 0.17% fat, and 0.07 to 0.10% protein.

Whole genome scan to detect quantitative trait loci for conformation and functional traits in dairy cattle

A granddaughter design was used to locate quantitative trait loci determining conformation and functional traits in dairy cattle. In this granddaughter design, consisting of 20 Holstein Friesian grandsires and 833 sons, genotypes were determined for 277 microsatellite markers covering the whole genome. Breeding values for 27 traits, regarding conformation (18), fertility (2), birth (4), workability (2), and udder health (1), were evaluated in an across-family analysis using multimarker regression. Significance thresholds were determined using a permutation test. The across-family analysis suggested the presence of 61 quantitative trait loci when 27 (i.e., one for each trait) were expected by chance. The test statistic exceeded the genomewise significance threshold for the following traits and chromosomes: chest width on chromosome 2; gestation length on chromosome 4; stature, body capacity, and size on chromosome 5; dairy character on chromosome 6; angularity on chromosome 12; fore udder attachment on chromosome 13; and fore udder attachment and front teat placement on chromosome 19. The quantitative trait loci for size traits on chromosomes 2, 5, and 6 may also have an effect on calving ease. The quantitative trait loci for udder traits on chromosomes 13 and 19 may also affect somatic cell score and mastitis resistance. If there are no negative effects on other economically important traits, marker assisted selection using markers associated with these quantitative trait loci can be applied.

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.

Candidate gene markers associated with somatotropic axis and milk selection

One of the obstacles to progress in dairy cattle selection is that milk production traits are only expressed after the first calving. However, the use of the quantitative trait loci (QTL) technology will improve the efficiency of dairy industry with a positive image for the consumers. QTL are part of the genome showing a preponderant action and explaining the major part of variation of the trait production. At the present time, the two major strategies developed to detect such QTL are the candidate gene approach and the positional genetics approach. The somatotropic axis contains the most promising candidates in this respect, as it strongly regulates milk production. Then, the identification of favorable QTL associated with the somatotropic axis that are significantly correlated with genetic merits for milk production could lead to more effective selection programs.

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.