Effects on production traits of haplotypes among casein genes in Norwegian goats and evidence for a site of preferential recombination

In goat milk the most abundant proteins are the casein genes, CSN1S1, CSN2, CSN1S2, and CSN3. Mutations have been identified within these genes affecting the level of gene expression, and effects on milk production traits have been reported. The aim of this study was to detect polymorphisms (SNPs) in the casein genes of Norwegian goats, resolve haplotype structures within the loci, and assess the effect of these haplotypes on milk production traits. Four hundred thirty-six Norwegian bucks were genotyped for 39 polymorphic sites across the four loci. The numbers of unique haplotypes present in each locus were 10, 6, 4, and 8 for CSN1S1, CSN2, CSN1S2, and CSN3, respectively. The effects of the CSN1S1 haplotypes on protein percentage and fat kilograms were significant, as were the effects of CSN3 haplotypes on fat percentage and protein percentage. A deletion in exon 12 of CSN1S1, unique to the Norwegian goat population, explained the effects of CSN1S1 haplotypes on fat kilograms, but not protein percentage. Investigation of linkage disequilibrium between all possible pairs of SNPs revealed higher levels of linkage disequilbrium for SNP pairs within casein loci than for SNP pairs between casein loci, likely reflecting low levels of intragenic recombination. Further, there was evidence for a site of preferential recombination between CSN2 and CSN1S2. The value of the haplotypes for haplotype-assisted selection (HAS) is discussed.     

Evidence for a geographical cline of casein haplotypes in Portuguese cattle breeds

Genetic variants of bovine milk proteins have been intensively used to characterize breeds and as markers for population/QTL studies throughout the world. However, a large number of cattle breeds including those found in Portugal, remain unstudied. In this work, we have analysed the genetic variation of six milk protein loci in 10 Portuguese cattle breeds by isoelectric focusing. High genetic diversities were generally found across breeds, with the exception of Mirandesa that showed a trend to fixation of the most common alleles in five loci, as well as of the rarer CSN3B allele. The casein haplotype BA2A was often the most frequent, followed by haplotypes BA2B and BA1A. Remarkably, CA2A was found to be the second most frequent haplotype in Southern breeds, supporting a geographical cline between Central-Northern European breeds and Bos indicus populations. Our data suggest that high genetic similarity among neighbouring Portuguese breeds is mainly caused by gene flow, and that the geographical distribution of particular casein haplotypes may indicate an influence of African cattle.     

Casein SNP in Norwegian goats: additive and dominance effects on milk composition and quality

Background

The four casein proteins in goat milk are encoded by four closely linked casein loci (CSN1S1, CSN2, CSN1S2 and CSN3) within 250 kb on caprine chromosome 6. A deletion in exon 12 of CSN1S1, so far reported only in Norwegian goats, has been found at high frequency (0.73). Such a high frequency is difficult to explain because the national breeding goal selects against the variant''s effect.

Methods

In this study, 575 goats were genotyped for 38 Single Nucleotide Polymorphisms (SNP) located within the four casein genes. Milk production records of these goats were obtained from the Norwegian Dairy Goat Control. Test-day mixed models with additive and dominance fixed effects of single SNP were fitted in a model including polygenic effects.

Results

Significant additive effects of single SNP within CSN1S1 and CSN3 were found for fat % and protein %, milk yield and milk taste. The allele with the deletion showed additive and dominance effects on protein % and fat %, and overdominance effects on milk quantity (kg) and lactose %. At its current frequency, the observed dominance (overdominance) effects of the deletion allele reduced its substitution effect (and additive genetic variance available for selection) in the population substantially.

Conclusions

The selection pressure of conventional breeding on the allele with the deletion is limited due to the observed dominance (overdominance) effects. Inclusion of molecular information in the national breeding scheme will reduce the frequency of this deletion in the population.     

Multiple quantitative trait locus analysis of bovine chromosome 6 in the Israeli Holstein population by a daughter design

Nine Israeli Holstein sire families with 2978 daughters were analyzed for quantitative trait loci effects on chromosome 6 for five milk production traits by a daughter design. All animals were genotyped for 2 markers. The three families with significant effects were genotyped for up to 10 additional markers spanning positions 0-122 cM of BTA6. Two sires were segregating for a locus affecting protein and fat percentage near position 55 cM with an estimated substitution effect of 0.18% protein, which is equivalent to one phenotypic standard deviation. This locus was localized to a confidence interval of 4 cM. One of these sires was also heterozygous for a locus affecting milk, fat, and protein production near the centromere. The hypothesis of two segregating loci was verified by multiple regression analysis. A third sire was heterozygous for a locus affecting milk and protein percentage near the telomeric end of the chromosome. Possible candidates for the major quantitative gene near position 55 cM were determined by comparative mapping. IBSP and SSP1 were used as anchors for the orthologous region on human chromosome 4. Twelve genes were detected within a 2-Mbp sequence. None of these genes have been previously associated with lactogenesis.     

Genetic polymorphism of αs1- and αs2-caseins in Hungarian Milking Goats

The aim of this study was to perform an initial characterization of milk quality and to determine genetic polymorphism at the CSN1S1 and CSN1S2 locus in two herds of local dairy goats (Hungarian Milking). The fat, protein and lactose level of milk samples in Hungarian Milking Goats were compared to other local goat breeds worldwide and it was concluded that the mean milk production of the Hungarian goats should be improved. The presence of the A, B, C + D, E, F and O alleles of CSN1S1 locus and A + B + C + E, D, F and O alleles of CSN1S2 locus were genotyped for by PCR-AS and PCR-RFLP methods in 103 goats. The strong B allele of CSN1S1 is more frequent in the local Hungarian Milking than in the imported Alpine and Saanen goats. The relatively high incidence of the O allele of CSN1S2 gene is also characteristic for the Hungarian Milking Goats and those special allele distribution patterns could be used to develop selection strategies to breed specialised lines of Hungarian local breeds.     

Linkage between amylase-1 locus and a major gene for milk fat content in cattle

Summary
Linkage between the amylase-1 ( Am-1 ) locus and a quantitative trait locus influencing fat content in milk was studied in offspring from heterozygous sires of the Swedish Red and White dairy breed. The effect on bull breeding values for fat content was estimated as interactions between sire and paternal Am-1 allele using a model eliminating the direct effects of sire and Am-1 allele. There were strong indications of linkage, confirming results of previous studies. The interaction was caused by strong associations in 7 out of 14 sire families. A test for within-family variance heterogeneity performed on the whole population of breeding bulls also supported the presence of a major gene for fat content in milk. The results indicate that there is genetic linkage between the Am-1 locus and a locus with large effect on milk fat content.     

Goat casein genotypes are associated with milk production traits in the Sarda breed

The aim of the current work was to analyze, in the Sarda breed goat, genetic polymorphism within the casein genes and to assess their influence on milk traits. Genetic variants at the CSN1S1, CSN2, CSN1S2 and CSN3 gene loci were investigated using PCR‐based methods, cloning and sequencing. Strong alleles prevailed at the CSN1S1 gene locus and defective alleles also were revealed. Null alleles were evidenced at each calcium‐sensitive gene locus. At the CSN3 gene locus, we observed a prevalence of the CSN3 A and B alleles; the occurrence of rare alleles such as CSN3 B'', C, C', D, E and M; and the CSN3 S allele (GenBank KF644565 ) described here for the first time in Capra hircus. Statistical analysis showed that all genes, except CSN3, significantly influenced milk traits. The CSN1S1 BB and AB genotypes were associated with the highest percentages of protein (4.41 and 4.40 respectively) and fat (5.26 and 5.34 respectively) (P < 0.001). A relevant finding was that CSN2 and CSN1S2 genotypes affected milk protein content and yield. The polymorphism of the CSN2 gene affected milk protein percentage with the highest values recorded in the CSN2 AA goats (4.35, at P < 0.001). The CSN1S2 AC goats provided the highest fat (51.02 g/day) and protein (41.42 g/day) (P < 0.01) production. This information can be incorporated into selection schemes for the Sarda breed goat.     

An allele associated with a non-detectable amount of αs2 casein in goat milk

The goat CSN1S2 locus is characterized by the presence of three alleles, A, B and C, all associated with about 2.5 g/l of protein per allele. The SDS-PAGE analysis of 441 individual milk samples obtained from goats belonging to a population reared in Southern Italy showed that the milk produced by three goats did not apparently contain alpha s2-casein, whereas milk produced by 37 goats showed a less intense electrophoretic band of this casein fraction (about 50%). These results can be explained by hypothesizing the presence of another allele at this locus, CSN1S2o, associated with a 'null' content of alpha s2-casein. Southern blot, PCR and PCR-RFLP analyses of the DNA region containing the CSN1S2 gene of individuals producing milk with and without alpha s2-casein did not show differences between the two groups. As a consequence, goats producing milk without alpha s2-casein carry an apparently intact gene. The first results obtained by sequencing part of the CSN1S2o allele revealed a G-->A transition at nucleotide 80 of the 11th exon which creates a stop codon and could be responsible for the absence of the alpha s2-casein in goat milk. This mutation eliminates a NcoI restriction site. A test based on this polymorphism has been established in order to identify carriers of the CSN1S2o allele.     

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.     

An MspI polymorphism at the bovine growth hormone (bGH) gene is linked to a locus affecting milk protein percentage

SSCP analysis of the bovine growth hormone (bGH) gene in Israel Holstein dairy cattle uncovered five intragenic haplotypes, denoted A to E. Of these, Haplotype E differed from the others at six fragments; one of which corresponded to the polymorphic MspI site in intron III, at which haplotype E carried the disabled MspI (-) allele. Haplotype E was observed in a single sire only, carrying haplotype A as the second bGH allele. In 523 daughters of this sire genotyped for the MspI polymorphism, heterozygous (+/-) as compared to homozygous (+/+) daughters, showed a significant increasing effect on protein percentage and kg protein per year; and a decreasing effect (P < 0.10) on milk somatic cell counts (MSSC). None of the daughters were homozygous (-/-), indicating that the frequency of this allele in the general population was essentially zero. Calculated skewness (g1) values for the two daughter groups differed significantly with (+/-) daughters showing negative skewness (in the direction of lower protein percentage), and (+/+) daughters positive skewness (in the direction of higher protein percentage). The direction of skewness in each group is indicative of the presence of a QTL having an increasing effect on milk protein percentage in coupling linkage with the MspI (-) allele in this sire, but at some distance from it. Maximum likelihood estimates of the proportion of recombination (r) between the putative QTL and bGH, and the allele substitution effect at the QTL (d), were r = 0.33, a = 0.07% protein, with standard errors 0.058 and 0.009% protein, respectively.     

An SNP in the goat CSN2 promoter region is associated with the absence of beta-casein in milk

So far, at least eight alleles in the goat CSN2 locus have been associated with the level of β -casein expression in milk. Alleles CSN2 ^A , CSN2 ^{A 1}, CSN2 ^B , CSN2 ^C , CSN2 ^D and CSN2 ^E have been associated with normal content (allele effects of about 5 g of β -casein per litre), whereas the CSN2 ⁰ and CSN2 ⁰¹ alleles have been associated with non-detectable levels of β -casein. Most of these alleles have been characterized genetically. Herein, we report the identification of a previously unreported SNP in the goat CSN2 promoter region ( AJ011018 :g.1311T>C), which is associated with the absence of β -casein in the milk. Furthermore, we developed a PCR-based method that allows detection of this mutation.     

Characterization of DGAT1 Allelic Effects in a Sample of North American Holstein Cattle

A putative causative mutation underlying a QTL was identified as a lysine to alanine non-conservative substitution at amino acid 232 of the gene encoding the acylCoA:diacylglycerol acyltransferase (DGAT1) protein. Our goal was to characterize the allelic substitution effects of this DGAT1 mutation in a large sample of Holstein bulls from North America. Statistically significant effects were identified for all of the milk production traits and somatic cell scores. Estimated average effects of substituting the lysine allele for the alanine variant on Holstein bull daughter yield deviations were ?81 kg, 3.7 kg, ?1.1 kg, 0.063%, 0.012%, and ?0.023 units for milk yield, fat yield, protein yield, fat component, protein component, and SCS, respectively. These estimates were largely in agreement with previous studies; however, the magnitudes of the estimates were much smaller in this study. Impacts on economic indices for net merit, cheese merit, and fluid merit were modest. Because of the strong antagonism between fat and protein yield and how those traits influence economic indices, selection for DGAT1 genotypes will likely not find widespread application in the U.S.     

Simultaneous mining of linkage and linkage disequilibrium to fine map quantitative trait loci in outbred half-sib pedigrees: revisiting the location of a quantitative trait locus with major effect on milk production on bovine chromosome 14

A maximum-likelihood QTL mapping method that simultaneously exploits linkage and linkage disequilibrium and that is applicable in outbred half-sib pedigrees is described. The method is applied to fine map a QTL with major effect on milk fat content in a 3-cM marker interval on proximal BTA14. This proximal location is confirmed by applying a haplotype-based association method referred to as recombinant ancestral haplotype analysis. The origin of the discrepancy between the QTL position derived in this work and that of a previous analysis is examined and shown to be due to the existence of distinct marker haplotypes associated with QTL alleles having large substitution effects.     

Casein haplotypes and their association with milk production traits in the Finnish Ayrshire cattle

Polymorphism of casein genes was studied in half-sib families of artificial insemination bulls of the Finnish Ayrshire dairy breed. Ten grandsires and 300 of their sons were genotyped for the following polymorphisms: α_S1-casein (B, C), β-casein (A¹, A²), the microsatellite within the K-casein gene (ms⁵, ms⁴) and K-casein (A, B, E). Nine different combinations of these alleles, casein haplotypes, were found. Associations between casein haplotypes and milk production traits (milk and protein yield, fat and protein percentage and milking speed) were studied with ordinary least-squares analysis to find a direct effect of the haplotypes or an association within individual grandsire families using the granddaughter design. Estimated breeding values of sons were obtained from cow evaluation by animal model. No direct effect of the casein haplotypes on the traits was found. Within grandsire families, in one out of four families the chromosomal segment characterized by haplotype 3 (B-A²-ms⁴-A) was associated with an increase in milk yield ( P <0.01) and a decrease in fat percentage ( P < 0.01) when contrasted with haplotype 8 (B-A¹-ms⁴-E). The results provide evidence that in the Finnish Ayrshire breed at least one quantitative trait locus affecting the genetic variation in yields traits is segregating linked to either haplotype 3 (B-A²-ms⁴-A) or 8 (B-A¹-ms⁴-E).     

The CSN1S1?N and F alleles identified by PCR-SSCP and their associations with milk yield and composition in Chinese dairy goats

A method was depicted to identify null allele CSN1S1 N and low allele CSN1S1 F of the CSN1S1 gene of goat using PCR-SSCP (polymerase chain reaction-single strand conformation polymorphism). First, primer A was designed to amplify the exon 9 of CSN1S1 gene which produced three genotypes AA, AB, and BB. Among these three genotypes, only AA and AB individuals had a cytosine deletion at exon 9 after DNA sequencing, which cannot be used to identify the N and F alleles. Therefore, primer B was used to amplify intron 14 of CSN1S1 of described AA and AB individuals. Genotypes FF, FN and NN were detected within AA individuals and genotypes FO and NO were detected in the above AB individuals. The frequencies of F and N alleles in 708 samples from Xinong Saanen (XS) and Guanzhong (GZ) dairy goat breeds were 0.1139, 0.0927, and 0.2376, 0.1193, respectively. In 268 XS samples, the individuals with NN genotype contained a significant lower protein content than that of other genotypes (P<0.01). Individuals of FF genotype had significant higher milk yield than that of NO genotype in the first milk lactation of 202 XS individuals (P<0.05). Therefore, the variability at CSN1S1 locus contains enough genetic diversity to be potentially useful in improving the quality and production of milk in Chinese dairy goat breeds.     

Association of the DGAT1 gene polymorphism in stud bulls with milk productivity in cows

The substitution of lysine for alanine (K232A) in the acyl-CoA-diacylglycerol acyltransferase, which is encoded by the DGAT1 gene, was tested for the significance for breeding evaluation of stud bulls of the holsteinized Black-and-White breed. The breeding value (deviation of milk productivity in daughters compared with cows of the same age) was estimated by the DYDC (daughter yield deviation to contemporary) method with modification. The frequency of allele 232K in the bulls examined was 0.28, lower than in Holstein bulls (0.4-0.6). The greatest effect of the A232K substitution was observed for the percent fat content in milk (1.4 sigma) and milk yield (0.76 sigma), and a lower effect was established for the milk protein yield (0.47 sigma) and percent protein content in milk (0.44 sigma). In the case of milk fat yield, the effect was nonsignificant. A method was proposed for converting the data on fat yield in order to obtain significant results in this case as well. The effect of the A232K substitution was estimated at 154 kg for milk yield, 2.8 kg for protein yield, 0.079% for percent fat content, and 0.015% for percent protein content. The results are discussed in the context of multiple pleiotropic effects of the K232A substitution in the DGAT1 gene. It is proposed that the K232A substitution of the DGAT1 gene be used as a golden standard in comparisons of the effect on milk productivity for the total gene set. This approach will allow a meta-analysis of the gene effects in spite of the differences in dairy cattle breeds and methods used to analyze their breeding value. In view of more than 30-year experience of using sperm of Holstein stud bulls, including those bred in North America, it was noted that the effect of the A232K substitution on milk productivity traits agreed well with the data reported for the North American commercial population of Holstein cattle.     

Additional support for an association between OLR1 and milk fat traits in cattle

The bovine oxidized LDL receptor 1 (OLR1) was chosen as a candidate gene for association tests with milk composition traits. Genotyping of 773 Italian Brown Swiss for a SNP at position 8232 in OLR1 (NW_215807:g.8232C>A) revealed a frequency of 0.95 for the g.8232C allele. The University of Wisconsin Holstein resource population was genotyped for the OLR1 NW_215807:g.[7160C>T; 7161A>G; 8232C>A] SNPs, and four haplotypes were inferred based on the genotypes of sires and their daughters. Oxidized LDL receptor 1 haplotypes were significantly associated with fat percentage (P = 0.0015). Haplotype [C; A; C] was associated with a significant increase in fat percentage when compared with the other haplotypes.     

Phenylalanine hydroxylase gene: novel missense mutation in exon 7 causing severe phenylketonuria.

By direct sequence analysis of 94 mutant phenylalanine hydroxylase alleles using polymerase chain reaction-based techniques, we identified a C to T transition in exon 7 of the human phenylalanine hydroxylase gene that is associated with RFLP haplotypes 1 and 4. A leucine for proline substitution at position 281 can be predicted from the nucleotide sequence of the mutant codon. Expression analysis in cultured mammalian cells after site-directed mutagenesis proved that the base substitution is a disease causing gene lesion. Dot-blot hybridization analysis using allele-specific oligonucleotides revealed that 25% of all mutant haplotype 1 alleles in the German population bear this mutation. In addition, this mutation could be detected on one mutant haplotype 4 allele. The fact that this mutation is associated with only 25% of all mutant haplotype 1 alleles suggests that multiple mutations may be associated with this haplotype. The occurrence of several different mutations would be in agreement with the clinical heterogeneity observed in the group of patients whose PKU alleles belong to haplotype 1.     

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

Individual loci affecting economically important traits can be located using genetic linkage between quantitative trait loci and genetic markers. In the ‘granddaughter’ experimental design, heterozygous grandsires and their sons are genotyped for the genetic marker, while the quantitative trait records of the granddaughters are used for statistical analysis. Ten DNA microsatellite markers were used to look for associations with quantitative trait loci affecting milk production traits in seven Israeli Holstein grandsire families. At least 60% more grandsires were heterozygous, and 40% fewer individuals were discarded because of unknown paternal allele origin, as compared with diallelic markers. The effects of paternal alleles for locus D21S4 on kg milk and protein were significant (P < 0.025). The allele substitution effects for sire 783 were 283 kg milk and 5.7 kg protein. For both traits, progeny of sire 783 that inherited allele ‘18’ had higher evaluations than progeny that inherited allele ‘21’. These results were verified by genotyping 151 of his daughters. Thus, the rate of genetic gain for protein production can be increased by selecting progeny of sire 783 carrying allele ‘18’ at this locus.