Whole genome scan to detect quantitative trait loci for bovine milk protein composition

The objective of this study was to perform a whole genome scan to detect quantitative trait loci (QTL) for milk protein composition in 849 Holstein–Friesian cows originating from seven sires. One morning milk sample was analysed for the major milk proteins using capillary zone electrophoresis. A genetic map was constructed with 1341 single nucleotide polymorphisms, covering 2829 centimorgans (cM) and 95% of the cattle genome. The chromosomal regions most significantly related to milk protein composition ( P _genome < 0.05) were found on Bos taurus autosomes (BTA) 6, 11 and 14. The QTL on BTA6 was found at about 80 cM, and affected α_S1-casein, α_S2-casein, β-casein and κ-casein. The QTL on BTA11 was found at 124 cM, and affected β-lactoglobulin, and the QTL on BTA14 was found at 0 cM, and affected protein percentage. The proportion of phenotypic variance explained by the QTL was 3.6% for β-casein and 7.9% for κ-casein on BTA6, 28.3% for β-lactoglobulin on BTA11, and 8.6% for protein percentage on BTA14. The QTL affecting α_S2-casein on BTA6 and 17 showed a significant interaction. We investigated the extent to which the detected QTL affecting milk protein composition could be explained by known polymorphisms in β-casein , κ -casein , β-lactoglobulin and DGAT1 genes. Correction for these polymorphisms decreased the proportion of phenotypic variance explained by the QTL previously found on BTA6, 11 and 14. Thus, several significant QTL affecting milk protein composition were found, of which some QTL could partially be explained by polymorphisms in milk protein genes.     

Conflicting candidates for cattle QTLs

Genome scans have identified quantitative trait loci (QTLs) affecting milk yield and composition in dairy cattle. For one QTL on bovine chromosome 6 (BTA6), previously fine-mapped to a 420-Kb region, mutations in two different genes (OPN and ABCG2) have been proposed as the underlying functional mutation. Comparing the arguments for each gene suggests that both mutations are equally probable. However, functional studies and/or additional populations are required to provide a definite answer.     

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.     

Comparative mapping of bovine chromosome 27 with human chromosome 8 near a dairy form QTL in cattle

In the absence of a complete and annotated bovine genome sequence, detailed human-bovine comparative maps are one of the most effective tools for identification of positional candidate genes contributing to quantitative trait loci (QTL) in cattle. In the present study, eight genes from human chromosome 8 were selected for mapping in cattle to improve breakpoint resolution and confirm gene order on the comparative map near the 40 cM region of the BTA27 linkage map where a QTL affecting dairy form had previously been identified. The resulting map identified ADRB3 as a positional candidate gene for the QTL contributing to the dairy form trait based on its estimated position between 40 and 45 cM on the linkage map. It is also a functional candidate gene due to its role in fat metabolism, and polymorphisms in the ADRB3 gene associated with obesity and metabolic disease in humans, as well as, carcass fat in sheep. Further studies are underway to investigate the existence of polymorphisms in the bovine ADRB3 gene and their association with traits related to fat deposition in cattle.     

Quantitative trait locus-by-environment interaction for milk yield traits on Bos taurus autosome 6

Genotype-by-environment interactions for production traits in dairy cattle have often been observed, while QTL analyses have focused on detecting genes with general effects on production traits. In this study, a QTL search for genes with environmental interaction for the traits milk yield, protein yield, and fat yield were performed on Bos taurus autosome 6 (BTA6), also including information about the previously investigated candidate genes ABCG2 and OPN. The animals in the study were Norwegian Red. Eighteen grandsires and 716 sires were genotyped for 362 markers on BTA6. Every marker bracket was regarded as a putative QTL position. The effects of the candidate genes and the putative QTL were modeled as a regression on an environmental parameter (herd year), which is based on the predicted herd-year effect for the trait. Two QTL were found to have environmentally dependent effects on milk yield. These QTL were located 3.6 cM upstream and 9.1 cM downstream from ABCG2. No environmentally dependent QTL was found to significantly affect protein or fat yield.     

Identification of a 3.7-Mb region for a marbling QTL on bovine chromosome 4 by identical-by-descent and association analysis

QTL mapping for growth and carcass traits was performed using a paternal half-sib family composed of 325 Japanese Black cattle offspring. Nine QTL were detected at the 1% chromosome-wise significance level at a false discovery rate of less than 0.1. These included two QTL for marbling on BTA 4 and 18, two QTL for carcass weight on BTA 14 and 24, two QTL for longissimus muscle area on BTA 1 and 4, two QTL for subcutaneous fat thickness on BTA 1 and 15 and one QTL for rib thickness on BTA 6. Although the marbling QTL on BTA 4 has been replicated with significant linkages in two Japanese Black cattle sires, the three Q (more marbling) haplotypes, each inherited maternally, were apparently different. To compare the three Q haplotypes in more detail, high-density microsatellite markers for the overlapping regions were developed within the 95% CIs (65 markers in 44–78 cM). A detailed haplotype comparison indicated that a small region (<3.7 Mb) around 46 cM was shared between the Qs of the two sires, whose dams were related. An association of this region with marbling was shown by a regression analysis using the local population, in which the two sires were produced and this was confirmed by an association study using a population collected throughout Japan. These results strongly suggest that the marbling QTL on BTA 4 is located in the 3.7-Mb region at around 46 cM.     

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. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users.     

Mapping of QTL for Body Conformation and Behavior in Cattle

Genome scans for quantitative trait loci (QTL) in farm animals have concentrated on primary production and health traits, and information on QTL for other important traits is rare. We performed a whole genome scan in a granddaughter design to detect QTL affecting body conformation and behavior in dairy cattle. The analysis included 16 paternal half-sib families of the Holstein breed with 872 sons and 264 genetic markers. The markers were distributed across all 29 autosomes and the pseudoautosomal region of the sex chromosomes with average intervals of 13.9 cM and covering an estimated 3155.5 cM. All families were analyzed jointly for 22 traits using multimarker regression and significance thresholds determined empirically by permutation. QTL that exceeded the experiment-wise significance threshold (5% level) were detected on chromosome 6 for foot angle, teat placement, and udder depth, and on chromosome 29 for temperament. QTL approaching experiment-wise significance (10% level) were located on chromosome 6 for general quality of feet and legs and general quality of udder, on chromosome 13 for teat length, on chromosome 23 for general quality of feet and legs, and on chromosome 29 for milking speed. An additional 51 QTL significant at the 5% chromosome-wise level were distributed over 21 chromosomes. This study provides the first evidence for QTL involved in behavior of dairy cattle and identifies QTL for udder conformation on chromosome 6 that could form the basis of recently reported QTL for clinical mastitis.     

Revisiting the quantitative trait loci for milk production traits on BTA6

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

Quantitative trait loci with effects on feed efficiency traits in Hereford × composite double backcross populations

Two half-sib families of backcross progeny were produced by mating F₁ Line 1 Hereford (L1) × composite gene combination (CGC) bulls with L1 and CGC cows. Feed intake and periodic weights were measured for 218 backcross progeny. These progenies were genotyped using 232 microsatellite markers that spanned the 29 BTA. Progeny from L1 and CGC females was analysed separately using composite interval mapping to find quantitative trait loci (QTL) affecting daily dry matter intake (DMI), average daily gain (ADG), feed conversion (FCR) and residual feed intake (RFI). Results from both backcrosses were pooled to find additional QTL. In the backcross to L1, QTL were detected for RFI and DMI on BTA11, FCR on BTA16, and ADG on BTA9. In the backcross to CGC, QTL were detected for RFI on BTA10, FCR on BTA12 and 16 and ADG on BTA15 and 17. After pooling, QTL were detected for RFI on BTA 2, 6, 7, 10, 11, 13 and 16; for FCR on BTA 9, 12, 16, 17 and 21; for ADG on BTA 9, 14, 15, 17; and for DMI on BTA 2, 5, 6, 9, 10, 11, 20 and 23.     

Genome wide QTL mapping to identify candidate genes for carcass traits in Hanwoo (Korean Cattle)

Meat quality traits are the most economically important traits affecting the beef industry in Korea. We performed a whole genome quantitative trait locus (QTL) mapping study of carcass data in Hanwoo Korean cattle. Two hundred sixty-six Hanwoo steers from 65 sires were genotyped using a 10K Affymetrix SNP chip. The average SNP interval across the bovine genome was 1.5Mb. Associations between each individual SNP and four carcass traits [carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling (MAR)] were assessed using a linear mixed model of each trait. Combined linkage and linkage disequilibrium analysis (LDLA) detected six potential QTL on BTA04, 06, 13, 16, 17, and 23 at the chromosome-wise level (P<0.05). Two MAR QTL were detected at 52.2 cM of BTA06 and 46.04 cM of BTA17. We identified three genes (ARAP2, LOC539460, and LOC511424) in the QTL region of BTA06 and seven genes (RPS14, SCARB1, LOC782103, BRI3BP, AACS, DHX37, and UBC) in the QTL region of BTA17. One significant QTL for CWT was detected at 100 cM on BTA04 and the corresponding QTL region spanned 1.7 cM from 99.7 to 101.4 cM. For EMA QTL, one significant QTL was detected at 3.9 cM of BTA23 and the most likely QTL interval was 1.4 cM, placing 15 candidate genes in the marker bracket. Finally, two QTL for BFT were identified at 68 cM on BTA13 and 24 cM on BTA16. The LPIN3 gene, which is functionally associated with lipodystrophy in humans, is located in the BFT QTL on BTA13. Thus, two potential candidate genes, acetoacetyl-CoA synthetase (AACS) and lipin (LPIN), were detected in QTL regions on BTA17 for MAR and BTA13 for BFT, respectively. In conclusion, LDLA analysis can be used to detect chromosome regions harboring QTL and candidate genes with a low density SNP panel, yielding relatively narrow confidence intervals regarding location.     

Forward to a detailed sequence map of bovine chromosome 6

Numerous QTL for a variety of phenotypic traits in dairy and beef cattle have been mapped on bovine chromosome 6 (BTA6). The complete and validated information on the molecular genome organization is an essential prerequisite for the conclusive identification of the causative sequence variation underlying the QTL. In our study we describe efforts to improve the genomic sequence map assembly of BTA6 by filling-in gaps and by suggesting sequence contig rearrangements. This is achieved by the generation and in silico mapping of BAC-end sequences (BESs) from clones containing sequences placed on our high-resolution radiation hybrid (RH) map of BTA6 onto the genome sequence map. Linking high-resolution RH mapping with in silico mapping of BESs on BTA6 enabled the detection of discrepancies in chromosomal assignments of genome sequence contigs and improved the resolution of non-conclusive assignments on the genome sequence assembly. Furthermore, 37% of BESs enabled chromosomal assignment of contigs previously unassigned. Anchoring of 66% of BESs onto HSA4 confirmed the synteny of the respective region of BTA6 including the known evolutionary breakpoints. The BESs will play an important role in the ongoing efforts to complete the sequence of the bovine genome and will also provide a source for the identification of new polymorphic sites in the genome sequence to resolve QTL-containing intervals.     

Genome Assembly Anchored QTL Map of Bovine Chromosome 14

Bovine chromosome 14 (BTA14) has been widely explored for quantitative trait loci (QTL) and genes related to economically important traits in both dairy and beef cattle. We reviewed more than 40 investigations and anchored 126 QTL to the current genome assembly (Btau 4_0). Using this anchored QTL map, we observed that, in dairy cattle, the region spanning 0 – 10 Mb on BTA14 has the highest density QTL map with a total of 56 QTL, mainly for milk production traits. It is very likely that both somatic cell score (SCS) and clinical mastitis share some common QTL in two regions: 61.48 Mb - 73.84 Mb and 7.86 Mb – 39.55 Mb, respectively. As well, both ovulation rate and twinning rate might share a common QTL region from 34.16 Mb to 65.38 Mb. However, there are no common QTL locations in three pregnancy related phenotypes: non-return rate, pregnancy rate and daughter pregnancy rate. In beef cattle, the majority of QTL are located in a broad region of 15 Mb – 45 Mb on the chromosome. Functional genes, such as CRH, CYP11B1, DGAT1, FABP4 and TG, as potential candidates for some of these QTL, were also reviewed. Therefore, our review provides a standardized QTL map anchored within the current genome assembly, which would enhance the process of selecting positional and physiological candidate genes for many important traits in cattle.     

Comparative map alignment of BTA27 and HSA4 and 8 to identify conserved segments of genome containing fat deposition QTL

Quantitative trait loci (QTL) associated with fat deposition have been identified on bovine Chromosome 27 (BTA27) in two different cattle populations. To generate more informative markers for verification and refinement of these QTL-containing intervals, we initiated construction of a BTA27 comparative map. Fourteen genes were selected for mapping based on previously identified regions of conservation between the cattle and human genomes. Markers were developed from the bovine orthologs of genes found on human Chromosomes 1 (HSA1), 4, 8, and 14. Twelve genes were mapped on the bovine linkage map by using markers associated with single nucleotide polymorphisms or microsatellites. Seven of these genes were also anchored to the physical map by assignment of fluorescence in situ hybridization probes. The remaining two genes not associated with an identifiable polymorphism were assigned only to the physical map. In all, seven genes were mapped to BTA27. Map information generated from the other seven genes not syntenic with BTA27 refined the breakpoint locations of conserved segments between species and revealed three areas of disagreement with the previous comparative map. Consequently, portions of HSA1 and 14 are not conserved on BTA27, and a previously undefined conserved segment corresponding to HSA8p22 was identified near the pericentromeric region of BTA8. These results show that BTA27 contains two conserved segments corresponding to HSA8p, which are separated by a segment corresponding to HSA4q. Comparative map alignment strongly suggests the conserved segment orthologous to HSA8p21-q11 contains QTL for fat deposition in cattle. Received: 25 February 2000 / Accepted: 30 March 2000     

QTL mapping in outbred half-sib families using Bayesian model selection

In this article, we propose a model selection method, the Bayesian composite model space approach, to map quantitative trait loci (QTL) in a half-sib population for continuous and binary traits. In our method, the identity-by-descent-based variance component model is used. To demonstrate the performance of this model, the method was applied to map QTL underlying production traits on BTA6 in a Chinese half-sib dairy cattle population. A total of four QTLs were detected, whereas only one QTL was identified using the traditional least square (LS) method. We also conducted two simulation experiments to validate the efficiency of our method. The results suggest that the proposed method based on a multiple-QTL model is efficient in mapping multiple QTL for an outbred half-sib population and is more powerful than the LS method based on a single-QTL model.     

A whole genome scan for QTL affecting milk protein percentage in Italian Holstein cattle, applying selective milk DNA pooling and multiple marker mapping in a daughter design

We report on a complete genome scan for quantitative trait loci (QTL) affecting milk protein percentage (PP) in the Italian Holstein-Friesian cattle population, applying a selective DNA pooling strategy in a daughter design. Ten Holstein-Friesian sires were chosen, and for each sire, about 200 daughters, each from the high and low tails of estimated breeding value for PP, were used to construct milk DNA pools. Sires and pools were genotyped for 181 dinucleotide microsatellites covering all cattle autosomes. Sire marker allele frequencies in the pools were obtained by shadow correction of peak height in the electropherograms. After quality control, pool data from eight sires were used for all subsequent analyses. The QTL heterozygosity estimate was lower than that of similar studies in other cattle populations. Multiple marker mapping identified 19 QTL located on 14 chromosomes (BTA1, 2, 3, 4, 5, 6, 8, 9, 12, 14, 17, 20, 23 and 27). The sires were also genotyped for seven polymorphic sites in six candidate genes (ABCG2, SPP1, casein kappa, DGAT1, GHR and PRLR) located within QTL regions of BTA6, 14 and 20 found in this study. The results confirmed or excluded the involvement of some of the analysed markers as the causative polymorphic sites of the identified QTL. The QTL identified, combined with genotype data of these candidate genes, will help to identify other quantitative trait genes and clarify the complex QTL patterns observed for a few chromosomes. Overall, the results are consistent with the Italian Holstein population having been under long-term selection for high PP.     

A refined radiation hybrid map of the telomeric region of bovine chromosome 18q25-q26 compared with human chromosome 19q13

Genome-wide scans have mapped economically important quantitative trait loci (QTL) for mastitis susceptibility in dairy cattle at the telomeric end of bovine chromosome 18 (BTA18). In order to increase the density of markers in this chromosomal region and to improve breakpoint resolution in the human-bovine comparative map, this study describes the chromosomal assignment of seven newly developed gene-associated markers and five microsatellites and eight previously mapped sequence tagged site markers near these QTL. The orientation of KCNJ14, BAX, CD37, NKG7, LIM2, PRKCG, TNNT1, MGC2705, RPL28, EPN1, ZNF582, ZIM2, STK13, ZNF132 and SLC27A5 on the 3000-rad radiation hybrid (RH) map of BTA18 is homologous to the organization found on the corresponding 10 Mbp of human chromosome 19q (HSA19q). The resulting bovine RH map with a length of 20.9 cR spans over about 11 cM on the bovine linkage map. The location of KCNJ14 and SLC27A5 flanking the RH map on BTA18q25-26 has been confirmed by fluorescence in situ hybridization. The data of this refined human-bovine comparative map should improve selection of candidate genes for mastitis susceptibility in dairy cattle.     

Identification of QTL for UV-protective eye area pigmentation in cattle by progeny phenotyping and genome-wide association analysis

Pigmentation patterns allow for the differentiation of cattle breeds. A dominantly inherited white head is characteristic for animals of the Fleckvieh (FV) breed. However, a minority of the FV animals exhibits peculiar pigmentation surrounding the eyes (ambilateral circumocular pigmentation, ACOP). In areas where animals are exposed to increased solar ultraviolet radiation, ACOP is associated with a reduced susceptibility to bovine ocular squamous cell carcinoma (BOSCC, eye cancer). Eye cancer is the most prevalent malignant tumour affecting cattle. Selection for animals with ACOP rapidly reduces the incidence of BOSCC. To identify quantitative trait loci (QTL) underlying ACOP, we performed a genome-wide association study using 658,385 single nucleotide polymorphisms (SNPs). The study population consisted of 3579 bulls of the FV breed with a total of 320,186 progeny with phenotypes for ACOP. The proportion of progeny with ACOP was used as a quantitative trait with high heritability (h(2)?=?0.79). A variance component based approach to account for population stratification uncovered twelve QTL regions on seven chromosomes. The identified QTL point to MCM6, PAX3, ERBB3, KITLG, LEF1, DKK2, KIT, CRIM1, ATRN, GSDMC, MITF and NBEAL2 as underlying genes for eye area pigmentation in cattle. The twelve QTL regions explain 44.96% of the phenotypic variance of the proportion of daughters with ACOP. The chromosomes harbouring significantly associated SNPs account for 54.13% of the phenotypic variance, while another 19.51% of the phenotypic variance is attributable to chromosomes without identified QTL. Thus, the missing heritability amounts to 7% only. Our results support a polygenic inheritance pattern of ACOP in cattle and provide the basis for efficient genomic selection of animals that are less susceptible to serious eye diseases.