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.     

Genome screen for twinning rate QTL in four North American Holstein families

The objective of this study was to identify twinning rate quantitative trait loci (QTL) by typing pooled samples in a preliminary screening followed by interval mapping to test QTL effects. Four elite North American Holstein half-sib sire families with high twinning rate predicted transmitting abilities (PTA) were used in this study. Chromosomes 5, 7, 19 and 23 were not genotyped as these chromosomes were scanned for QTL in these families in a previous study. DNA was extracted from phenotypically extreme sons in each sire family. Two pools were prepared from sons of sires in each phenotypic tail, two each for high and low PTA levels for twinning rates. Each pool contained DNA from 4 to 15% of all sons of the sire depending on family. A total of 268 fluorescently labelled microsatellite markers were tested for heterozygosity in sires. About 135--170 informative markers per family were genotyped using pooled DNA samples. Based on the preliminary evidence for potential twinning rate QTL from pooled typing, interval mapping was performed subsequently on 12 chromosomal regions by family combinations. Evidence of QTL for twinning rate was found in one family on BTA 21 and 29 at a chromosome-wide P<0.05 and on BTA 8, 10 and 14 with a chromosome-wide P<0.01.     

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.     

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.     

Detection of quantitative trait loci for meat quality traits in cattle

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

Alternative parameterizations of relatedness in whole genome association analysis of pre-weaning traits of Nelore-Angus calves

Gestation length, birth weight, and weaning weight of F₂ Nelore-Angus calves (n = 737) with designed extensive full-sibling and half-sibling relatedness were evaluated for association with 34,957 SNP markers. In analyses of birth weight, random relatedness was modeled three ways: 1) none, 2) random animal, pedigree-based relationship matrix, or 3) random animal, genomic relationship matrix. Detected birth weight-SNP associations were 1,200, 735, and 31 for those parameterizations respectively; each additional model refinement removed associations that apparently were a result of the built-in stratification by relatedness. Subsequent analyses of gestation length and weaning weight modeled genomic relatedness; there were 40 and 26 trait-marker associations detected for those traits, respectively. Birth weight associations were on BTA14 except for a single marker on BTA5. Gestation length associations included 37 SNP on BTA21, 2 on BTA27 and one on BTA3. Weaning weight associations were on BTA14 except for a single marker on BTA10. Twenty-one SNP markers on BTA14 were detected in both birth and weaning weight analyses.     

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.     

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.     

Genome-wide identification of quantitative trait loci for carcass composition and meat quality in a large-scale White Duroc × Chinese Erhualian resource population

Carcass and meat quality traits are economically important in pigs. In this study, 17 carcass composition traits and 23 meat quality traits were recorded in 1028 F₂ animals from a White Duroc × Erhualian resource population. All pigs in this experimental population were genotyped for 194 informative markers covering the entire porcine genome. Seventy-seven genome-wide significant quantitative trait loci (QTL) for carcass traits and 68 for meat quality were mapped to 34 genomic regions. These results not only confirmed many previously reported QTL but also revealed novel regions associated with the measured traits. For carcass traits, the most prominent QTL was identified for carcass length and head weight at 57 cM on SSC7, which explained up to 50% of the phenotypic variance and had a 95% confidence interval of only 3 cM. Moreover, QTL for kidney and spleen weight and lengths of cervical vertebrae were reported for the first time in pigs. For meat quality traits, two significant QTL on SSC5 and X were identified for both intramuscular fat content and marbling score in the longissimus muscle, while three significant QTL on SSC1 and SSC9 were found exclusively for IMF. Both LM and the semimembranous muscle showed common QTL for colour score on SSC4, 5, 7, 8, 13 and X and discordant QTL on other chromosomes. White Duroc alleles at a majority of QTL detected were favourable for carcass composition, while favourable QTL alleles for meat quality originated from both White Duroc and Erhualian.     

Genome‐scan analysis for genetic mapping of quantitative trait loci underlying birth weight and onset of puberty in doe kids (Capra hircus)

The objective of this study was to locate quantitative trait loci (QTL) causing variation in birth weight and age of puberty of doe kids in a population of Rayini cashmere goats. Four hundred and thirty kids from five half‐sib families were genotyped for 116 microsatellite markers located on the caprine autosomes. The traits recorded were birth weight of the male and female kids, body weight at puberty, average daily gain from birth to age of puberty and age at puberty of the doe kids. QTL analysis was conducted using the least squares interval mapping approach. Linkage analysis indicated significant QTL for birth weight on Capra hircus chromosomes (CHI) 4, 5, 6, 18 and 21. Five QTL located on CHI 5, 14 and 29 were associated with age at puberty. Across‐family analysis revealed evidence for overlapping QTL affecting birth weight (78 cM), body weight at puberty (72 cM), average daily gain from birth to age of puberty (72 cM) and age at puberty (76 cM) on CHI 5 and overlapping QTL controlling body weight at puberty and age at puberty on CHI 14 at 18–19 cM. The proportion of the phenotypic variance explained by the detected QTL ranged between 7.9% and 14.4%. Confirming some of the previously reported results for birth weight and growth QTL in goats, this study identified more QTL for these traits and is the first report of QTL for onset of puberty in doe kids.     

Investigation of the genetic architecture of a bone carcass weight QTL on BTA6

A previous analysis of an F₂/Backcross Charolais × Holstein cross population identified the presence of a highly significant QTL on chromosome 6 (BTA6) affecting the proportion of bone in the carcass. Two closely linked QTL affected birth weight (BW) and body length at birth (BBL). In this report, the marker density around the QTL on BTA6 was increased, adding four additional microsatellite markers across the chromosome and 46 SNPs within the target QTL confidence interval. Of the SNPs, 26 were in positional candidate genes and the remaining 20 provided an even distribution of markers in the target QTL region. As a bone‐related trait, the sum of the bone weight for all the left fore‐ and hindquarter joints of the carcass was analysed. We also studied the BW and BBL. Analyses of the data substantially reduced the QTL confidence interval. No strong evidence was found that the QTL for the three traits studied are different, and we conclude that the results are consistent with a single pleiotropic QTL influencing the three traits, with the largest effects on the proportion of bone in the carcass. The analyses also suggest that none of the SNPs tested is the sole causative variant of the QTL effects. Specifically, the SNP in the NCAPG gene previously reported as a causal mutation for foetal growth and carcass traits in other cattle populations was excluded as the causal mutation for the QTL reported here. Polymorphisms located in other previously identified candidate genes including SPP1, ABCG2, IBSP, MEPE and PPARGC1A were also excluded. The results suggest that SNP51_BTA‐119876 is the polymorphism in strongest linkage disequilibrium with the causal mutation(s). Further research is required to identify the causal variant(s) associated with this bone‐related QTL.     

Detection of QTL for milk production traits in cattle by application of a specifically developed marker map of BTA6

Interval mapping was carried out to identify quantitative trait loci (QTL) for milk production traits in five granddaughter design families of the German Holstein population. Fourteen randomly generated markers spanning the whole of BTA6 and six targeted microsatellite markers from BTA6q21-31 were included in the analysis. In one family a QTL with effects on milk fat yield and milk protein yield was mapped to the interval TGLA37-FBN13 (3 CM proximal to FBN13, lodscore 3.22) in the middle part of the chromosome. Although there are several reports about QTL with effects on milk production traits on BTA6 in the literature, a QTL with effects on milk fat and milk protein yield has not been previously described.     

Sonographic biometrical normograms and estimation of fetal weight in the baboon (Papio anubis)

Background  In order to consider the non-human primate as an adequate model for studying prenatal diagnosis and therapy, comparative data on fetal growth should be available.
Methods  Sixty ultrasound scans were performed in 22 baboons between 14 and 167 days of gestation. Measurements included greatest length, head circumference, biparietal diameter (BPD), transcerebellar diameter, abdominal circumference (AC), femur length (FL), and amniotic fluid index. For all parameters growth curves were established and compared with human curves. In 18 animals, birth weight and placental weight were determined. Different equations described in the literature for estimating the human fetal weight were tested in the baboon.
Results  The fetal and placental growth pattern in the baboon was comparable with humans. The best predictor of fetal weight was the formula presented by Combs: 0.23966 × AC² × FL + 1.623 × BPD³.
Conclusions  A high similarity between baboon and human growth charts is shown. The best equation for estimating the baboon fetal weight is proposed.     

Effects of calf breed on milk production and other economic traits of Holstein dams

Production of purebred or crossbred feeder calves for beef, especially HolsteinxJapanese-Black (HxJB) and Japanese-Black (JB), from dairy cattle using artificial insemination or embryo transfer have been used widely in Japan. However, dairy farmers feel uneasy about the effects of calf breed on the economic traits of dams. In this study, those effects were investigated in 798 Holstein heifers bearing Holstein, HxJB, JB or other breed calves. The results of the least-squares ANOVA indicated the effects of calf breed to be significant for gestation length (P<0.01) and calf birth weight (P<0.01) but not for milk yield, fat yield, protein yield, peak yield, day of peak, number of artificial inseminations per pregnancy or days open. Frequency of dystocia was lower in dams bearing HxJB calves than in those bearing Holstein calves (P<0.05). There were no significant differences for frequencies of still birth, retained placenta and subsequent pregnancy. The present data suggest that the effects of calf breed do not place a serious problem on the economic traits of Holstein dams. In conclusion, it is indicated that the production of JB and HxJB feeder calves from Holstein dams does not result in a decrease in dam productivity.     

A genome scan for quantitative trait loci affecting three ear traits in a White Duroc × Chinese Erhualian resource population

Chinese Erhualian pigs have larger and floppier ears compared with White Duroc pigs (small, half- or fully-pricked ears). To identify quantitative trait loci (QTL) for ear weight and area as well as erectness, a genome-wide scan with 194 microsatellites was performed in a White Duroc × Chinese Erhualian resource population (>1000 F₂ animals). Twenty-three genome-wide significant QTL and 12 suggestive QTL were identified. All QTL for ear erectness and size detected in two previous studies, bar two on SSC6 and 9, were confirmed here. The 1% genome-wide significant QTL at 70 cM on SSC5 and at 58 cM on SSC7 have profound and pleiotropic effects on the three ear traits, with Erhualian alleles increasing weight and area but decreasing erectness. Notably, the 95% confidence interval of the QTL for weight and area on SSC7 spanned only 3 cM. New QTL reaching 1% genome-wide significance were found on SSC8 (at 37 cM) for all three ear traits, on SSC4 and 16 for weight and area, and on SSCX for area. Unexpectedly, Erhualian alleles at these loci were associated with lighter and smaller or erect ear. Some new suggestive QTL were also found on other chromosome regions. Almost all the QTL for weight and area had essentially additive effects, while the QTL for erectness on SSC2, 5 and 7 showed not only additive effects but also partial dominance effects of Erhualian alleles. The two most significant QTL on SSC7 and SSC5 could be promising targets for fine mapping and identification of the causative mutations.     

Genome-wide scan for bovine twinning rate QTL using linkage disequilibrium

Twinning is a complex trait with negative impacts on health and reproduction, which cause economic loss in dairy production. Several twinning rate quantitative trait loci (QTL) have been detected in previous studies, but confidence intervals for QTL location are broad and many QTL are unreplicated. To identify genomic regions or genes associated with twinning rate, QTL analysis based on linkage combined with linkage disequilibrium (LLD) and individual marker associations was conducted across the genome using high-throughput single nucleotide polymorphism (SNP) genotypes. A total of 9919 SNP markers were genotyped with 200 sires and sons in 19 half-sib North American Holstein dairy cattle families. After SNPs were genotyped, informative markers were selected for genome-wide association tests and QTL searches. Evidence for twinning rate QTL was found throughout the genome. Thirteen markers significantly associated with twinning rate were detected on chromosomes 2, 5 and 14 ( P  < 2.3 × 10⁻⁵). Twenty-six regions on fourteen chromosomes were identified by LLD analysis at P  < 0.0007. Seven previously reported ovulation or twinning rate QTL were supported by results of single marker association or LLD analyses. Single marker association analysis and LLD mapping were complementary tools for the identification of putative QTL in this genome scan.     

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.     

Mapping quantitative trait loci regulating chicken body composition traits

Genome scans were conducted on an F₂ resource population derived from intercross of the White Plymouth Rock with the Silkies Fowl to detect QTL affecting chicken body composition traits. The population was genotyped with 129 microsatellite markers and phenotyped for 12 body composition traits on 238 F₂ individuals from 15 full-sib families. In total, 21 genome-wide QTL were found to be responsible for 11 traits, including two newly studied traits of proventriculus weight and shank girth. Three QTL were genome-wide significant: at 499 c m on GGA1 (explained 3.6% of phenotypic variance, P  < 0.01) and 51 c m on GGA5 (explained 3.3% of phenotypic variance, P  < 0.05) for the shank & claw weight and 502 c m on GGA1 (explained 1.4% of phenotypic variance, P  < 0.05) for wing weight. The QTL on GGA1 seemed to have pleiotropic effects, also affecting gizzard weight at 490 c m , shank girth at 489 c m and intestine length at 481 c m . It is suggested that further efforts be made to understand the possible pleiotropic effects of the QTL on GGA1 and that on GGA5 for two shank-related traits.     

The use of crossbreeding with beef bulls in dairy herds: effects on calving difficulty and gestation length

This study was designed to analyse the evolution in the use of beef bull semen for dairy cattle insemination and, mainly, to assess calving difficulty, gestation length and proportion of stillbirths after breeding pure Holsteins or crossbreeding. Data were collected during 2004 to 2011 for 552 571 Holstein calvings (457 070 Holstein × Holstein, 43 384 Holstein × Limousine, 32 174 Holstein × Belgian Blue and 19 943 Holstein × Galician Blonde). The highest calving difficulty, compared with pure Holsteins was for crosses with Belgian Blue followed by Limousine and Galician Blonde. The Holstein × Limousine and Holstein × Galician Blonde crossbred calves had significantly longer gestation lengths than Holstein × Holstein and Holstein × Belgian Blue calves. Between the latter two, pure Holstein had the shortest gestation length. Calving difficulty and gestation length decreased as the age of the dam advanced. The most difficult calvings were observed in twin calvings, followed by the calvings of male calves and female calves. The gestations leading to the birth of male calves were longer than those leading to female calves and twin calves. Stillbirths were not related to the breed used for mating. Through examining these parameters, sire breed should be considered when selecting a beef breed for the insemination of milk-producing dams.     

Quantitative trait loci (QTL) for lean body mass and body length in MRL/MPJ and SJL/J F(2) mice

Studies on the genetic mechanisms involved in the regulation of lean body mass (LBM) in mammals are minimal, although LBM is associated with a competent immune system and an overall good (healthy) body functional status. In this study, we performed a high-density genome-wide scan using 633 (MRL/MPJ × SJL/J) F₂ intercross to identify the quantitative trait loci (QTL) involved in the regulation of LBM. We hypothesized that additional QTL can be identified using a different mouse cross (MRL/SJL cross). Ten QTL were identified for LBM on chromosomes (chrs) 2, 6, 7, 9,13 and 14. Of those ten, QTL on chrs 6, 7 and 14 were exclusive to LBM, while QTL on chrs 4 and 11 were exclusively body length. LBM QTL on chrs 2 and 9 overlap with those of size. Altogether, the ten LBM QTL explained 41.2% of phenotypic variance in F₂ mice. Five significantly interacting loci that may be involved in the regulation of LBM were identified and accounted for 24.4% of phenotypic variance explained by the QTL. Five epistatic interactions, contributing 22.9% of phenotypic variance, were identified for body length. Interacting loci on chr 2 may influence LBM by regulating body length. Therefore, epistatic interactions as well as single QTL effects play an important role in the regulation of LBM. Electronic Publication