Effect of κ-Casein Polymorphism on Milk Composition in the Orobica Goat

Colostrum Quality in Different Goat Breeds Reared in Northern Italy

This study aimed to characterize the colostrum quality in three different local goat breeds of Northern Italy (i.e., Frisa Valtellinese, Orobica, and Lariana) and a cosmopolitan one (i.e., Camosciata delle Alpi) (n = 30 per breed), reared under traditional semi-extensive and intensive systems, respectively. Lariana showed the highest percentage of fat (10.18 ± 3.14%) and total solids (30.73 ± 4.89%) but the lowest percentage of lactose (1.87 ± 0.82%; p < 0.05); Orobica had the lowest percentage of fat (7.13 ± 2.48%), total solids (24.11 ± 5.48%), and protein (10.77 ± 4.53%) but the highest percentage of lactose (3.16 ± 0.73%; p < 0.05). This suggests that breeds which have a more pronounced meat aptitude (i.e., Frisa and Lariana) have a higher concentration of components than breeds with more dairy aptitude (i.e., Orobica and Camosciata). Uni- and multivariate analyses showed that IgG is the parameter that best differentiates local breeds from cosmopolitan ones (p < 0.01). Colostrum from Frisa goats showed the highest IgG concentration (100.90 ± 8.11 mg/mL), while the lowest concentration was in the Camosciata breed (74.75 ± 20.16 mg/mL). Finally, the highest lactoferrin concentration was in Frisa (1781.3 ± 892.6 µg/mL) and the lowest in Camosciata and Lariana (763.1 ± 357.9 and 1148.0 ± 858.6 µg/mL, respectively; p < 0.05). Differences between Camosciata and local breeds could be due to the different farming systems, in addition to the genetic characteristics. The higher quality of colostrum produced by some local goats could be an adaptive characteristic that helps the growth and survival of the kids.

Review: Genetic and protein variants of milk caseins in goats

The milk casein genes in goats, are highly polymorphic genes with numerous synonymous and non-synonymous mutations. So far, 20 protein variants have been reported in goats for alpha-S1-casein, eight for beta-casein, 14 for alpha-S2-casein, and 24 for kappa-casein. This review provides a comprehensive overview on identified milk casein protein variants in goat and non-coding DNA sequence variants with some affecting the expression of the casein genes. The high frequency of some casein protein variants in different goat breeds and geographical regions might reflect specific breeding goals with respect to milk processing characteristics, properties for human nutrition and health, or adaptation to the environment. Because protein names, alongside the discovery of protein variants, go through a historical process, we linked old protein names with new ones that reveal more genetic variability. The haplotypes across the cluster of the four genetically linked casein genes are recommended as a valuable genetic tool for discrimination between breeds, managing genetic diversity within and between goat populations, and breeding strategies. The enormous variation in the casein proteins and genes is crucial for producing milk and dairy products with different properties for human health and nutrition, and for genetic improvement depending on local breeding goals.

Does the Acknowledgement of αS1-Casein Genotype Affect the Estimation of Genetic Parameters and Prediction of Breeding Values for Milk Yield and Composition Quality-Related Traits in Murciano-Granadina?

Simple Summary

Genetic evaluations and the selection of breeding animals require the accurate estimation of genetic parameters for economically important traits. As a result, dairy livestock has evolved in response to the needs of producers and consumers. Genetic selection in goats has been mostly based on quantitative traits such as milk yield, fat, protein, and dry matter. However, as reported by the increased heritability values of these parameters after the inclusion of the different allelic variants of αS1 casein in evaluation models, the selection of animals carrying this gene could result in a more efficient genetic selection. High levels of genetic polymorphism (89.58% of polymorphic SNP—as only five out of the 48 SNPs assessed were monomorphic) that are related to greater production of coagulable proteins in milk, a fact that could be associated with a higher yield and improved curd firmness properties.

Abstract

A total of 2090 lactation records for 710 Murciano-Granadina goats were collected during the years 2005–2016 and analyzed to investigate the influence of the αS1-CN genotype on milk yield and components (protein, fat, and dry matter). Goats were genetically evaluated, including and excluding the αS1-CN genotype, in order to assess its repercussion on the efficiency of breeding models. Despite no significant differences being found for milk yield, fat and dry matter heritabilities, protein production heritability considerably increased after aS1-CN genotype was included in the breeding model (+0.23). Standard errors suggest that the consideration of genotype may improve the model’s efficiency, translating into more accurate genetic parameters and breeding values (PBV). Genetic correlations ranged from −0.15 to −0.01 between protein/dry matter and milk yield/protein and fat content, while phenotypic correlations were −0.02 for milk/protein and −0.01 for milk/fat or protein content. For males, the broadest range for reliability (R_AP) (0.45–0.71) was similar to that of females (0.37–0.86) when the genotype was included. PBV ranges broadened while the maximum remained similar (0.61–0.77) for males and females (0.62–0.81) when the genotype was excluded, respectively. Including the αS1-CN genotype can increase production efficiency, milk profitability, milk yield, fat, protein and dry matter contents in Murciano-Granadina dairy breeding programs.

A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content

The quantity of milk and milk fat and proteins are particularly important traits in dairy livestock. However, little is known about the regions of the genome that influence these traits in goats. We conducted a genome wide association study in French goats and identified 109 regions associated with dairy traits. For a major region on chromosome 14 closely associated with fat content, the Diacylglycerol O-Acyltransferase 1 (DGAT1) gene turned out to be a functional and positional candidate gene. The caprine reference sequence of this gene was completed and 29 polymorphisms were found in the gene sequence, including two novel exonic mutations: R251L and R396W, leading to substitutions in the protein sequence. The R251L mutation was found in the Saanen breed at a frequency of 3.5% and the R396W mutation both in the Saanen and Alpine breeds at a frequencies of 13% and 7% respectively. The R396W mutation explained 46% of the genetic variance of the trait, and the R251L mutation 6%. Both mutations were associated with a notable decrease in milk fat content. Their causality was then demonstrated by a functional test. These results provide new knowledge on the genetic basis of milk synthesis and will help improve the management of the French dairy goat breeding program.

Single nucleotide polymorphism of candidate genes in non-descript local goats of Sri Lanka

In the present study, genetic polymorphism in exon 4 of kappa casein (k-CSN3), exon 2-3 of alpha lactalbumin (LALBA) and exon 1 of gonadotropin releasing hormone receptor (GnRHR) genes were analyzed as candidate genes for milk production, milk quality and prolificacy in non-descriptive local goats in Sri Lanka. Altogether eleven, one and three single nucleotide polymorphisms (SNPs) were identified in k-CSN3, LALBA and GnRHR gene fragments, respectively utilizing the DNA sequencing technique for the first time in Sri Lanka. Seven polymorphic sites out of eleven in k-CSN3 gene fragment and the recorded variable site in LALBA gene fragment were homozygous while all three polymorphic sites in GnRHR gene fragment were heterozygous. Two of the SNPs recorded in the present study are known to be unique for Sri Lankan non-descript goat population at G203T and A730G in k-CSN3 and GnRHR genes, respectively. The study records another two SNPs in GnRHR gene, which are already known to be correlated with higher fecundity in goats (G757A and G891T). Results of the present study will be extremely important in future attempts to develop markers to improve the milk production, milk composition and litter size of non-descript local goats in Sri Lanka.

The Application of Genomic Technologies to Investigate the Inheritance of Economically Important Traits in Goats

Goat genomics has evolved at a low pace because of a lack of molecular tools and sufficient investment. Whilst thousands and hundreds of quantitative trait loci (QTL) have been identified in cattle and sheep, respectively, about nine genome scans have been performed in goats dealing with traits as conformation, growth, fiber quality, resistance to nematodes, and milk yield and composition. In contrast, a great effort has been devoted to the characterization of candidate genes and their association with milk, meat, and reproduction phenotypes. In this regard, causal mutations have been identified in the αS1-casein gene that has a strong effect on milk composition and the PIS locus that is linked to intersexuality and polledness. In recent times, the development of massive parallel sequencing technologies has allowed to build a reference genome for goats as well as to monitor the expression of mRNAs and microRNAs in a broad array of tissues and experimental conditions. Besides, the recent design of a 52K SNP chip is expected to have a broad impact in the analysis of the genetic architecture of traits of economic interest as well as in the study of the population structure of goats at a worldwide scale.

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.

Newly identified mutations at the CSN1S1 gene in Ethiopian goats affect casein content and coagulation properties of their milk

Very high casein content and good coagulation properties previously observed in some Ethiopian goat breeds led to investigating the αs1-casein (CSN1S1) gene in these breeds. Selected regions of the CSN1S1 gene were sequenced in 115 goats from 5 breeds (2 indigenous: Arsi-Bale and Somali, 1 exotic: Boer, and 2 crossbreeds: Boer × Arsi-Bale and Boer × Somali). The DNA analysis resulted in 35 new mutations: 3 in exons, 3 in the 5' untranslated region (UTR), and 29 in the introns. The mutations in exons that resulted in an amino acid shift were then picked to evaluate their influence on individual casein content (αs1-, αs2-, β-, and κ-CN), micellar size, and coagulation properties in the milk from the 5 goat breeds. A mutation at nucleotide 10657 (exon 10) involved a transversion: CAG→CCG, resulting in an amino acid exchange Gln77→Pro77. This mutation was associated with the indigenous breeds only. Two new mutations, at nucleotide 6072 (exon 4) and 12165 (exon 12), revealed synonymous transitions: GTC→GTT in Val15 and AGA→AGG in Arg100 of the mature protein. Transitions G→A and C→T at nucleotides 1374 and 1866, respectively, occurred in the 5' UTR, whereas the third mutation involved a transversion T→G at nucleotide location 1592. The goats were grouped into homozygote new (CC), homozygote reference (AA), and heterozygote (CA) based on the nucleotide that involved the transversion. The content of αs1-CN (15.32g/kg) in milk samples of goats homozygous (CC) for this newly identified mutation, Gln77→Pro77 was significantly higher than in milks of heterozygous (CA; 9.05g/kg) and reference (AA; 7.61g/kg) genotype animals. The αs2-, β-, and κ-CN contents showed a similar pattern. Milk from goats with a homozygous new mutation had significantly lower micellar size. Milk from both homozygote and heterozygote new-mutation goats had significantly shorter coagulation rate and stronger gel than the reference genotype. Except the transversion, the sequence corresponded to allele A and presumably derived from it. Therefore, this allele is denoted by A3. All goats from the reference genotype (AA) were homozygous for the allele at nucleotide position 1374 and 1866, whereas all mutations in the 5' UTR existed in a heterozygous form in both heterozygous (CA) and the new mutation (CC) genotype. The newly identified mutation (CC) detected in some of the goat breeds is, therefore, important in selection for genetic improvement and high-quality milk for the emerging goat cheese-producing industries. The finding will also benefit farmers raising these goat breeds due to the increased selling price of goats. Further studies should investigate the effect of this amino acid exchange on the secondary and tertiary structure of the αs1-CN molecule and on the susceptibility of peptide hydrolysis by digestive enzymes.

Identification of an intronic regulatory mutation at the buffalo αS1-casein gene that triggers the skipping of exon 6

The characterization of casein polymorphism is an essential step in order to understand the genetic basis of milk quality in dairy ruminants. In this work, we report the identification of a regulatory mutation at the buffalo αs1-casein (CSN1S1) gene that alters the normal processing of the primary transcript. Sequencing of CSN1S1 cDNA from individuals harbouring this new variant revealed that its most distinctive feature is the loss of exon 6 that encodes eight amino acids between positions 35-42 of mature protein. In an effort to map the causal mutation, we sequenced a genomic region spanning exons 5-7 of the buffalo CSN1S1 gene. This experiment allowed us to establish that exon 6-skipping is produced by a G to C substitution at the first position of intron 6 that inactivates the donor splice site. This mutation can be typed by PCR-RFLP by using either TaaI or Bpu10I diagnostic restriction enzymes, and it has a frequency of 0.18 in Romanian buffaloes. This exon skipping phenomenon is the first one described in buffalo CSN1S1 locus.

Associations of acetyl-coenzyme A carboxylase α, stearoyl-coenzyme A desaturase, and lipoprotein lipase genes with dairy traits in Alpine goats

Milk yield and composition are of great economic importance for the dairy goat industry. The identification of genes associated with phenotypic differences for these traits could allow for the implementation of gene-assisted selection programs in goats. Associations between polymorphisms at 3 candidate genes and milk production traits in Alpine goats farmed in Italy were investigated in the present research. Considered genes were acetyl-coenzyme A carboxylase α (ACACA), the major regulatory enzyme of fatty acid biosynthesis; stearoyl-coenzyme A desaturase (SCD), involved in the biosynthesis of monounsaturated fatty acids in the mammary gland; and lipoprotein lipase (LPL), which plays a central role in plasma triglyceride metabolism. An approach somewhat similar to the granddaughter design for detecting quantitative trait loci in dairy cattle was followed. Effects of genotypes of a sample of 59 Alpine bucks on phenotypes of their 946 daughters raised in 75 flocks were investigated. Data comprised 13,331 daily records for milk yields (L/d), fat and protein yields (kg/d), and fat and protein contents (%) of 2,200 lactations. Population genetics parameters were calculated and associations between milk production traits and 10 single nucleotide polymorphisms (SNP) at the 3 genes were tested. Two markers at the ACACA, 1 for the SCD and 1 at the LPL locus, deviated significantly from the Hardy-Weinberg equilibrium, with an observed heterozygosity lower than expected. Flock, age of the goat, kidding season, and stage of lactation affected all traits considered, except fat percentage. Three SNP were found to be significantly associated with milk production traits. The SNP located on the ACACA gene showed an effect on milk yield, with daughters of TT bucks having an average test-day milk yield of about 0.3 to 0.25 L/d lower than the other 2 genotypes. The marker on the LPL locus was highly associated with milk yield, with the largest values for CC daughters (about 0.50L more than GG). The TGT deletion located on the untranslated region of the SCD gene showed significant effects on average milk and protein yields. The homozygote-deleted genotype had values about 0.5 L/d and 16 g/d lower for milk and protein daily yield, respectively, compared with the TGT/TGT genotype. Differences between genotypes were quite constant across most of the lactation. Associations found in the present study, which should be tested in a larger sample, especially for those markers that show rare genotypes, may offer useful indications for the genetic improvement of dairy traits in goats.

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.

Effects of αs1-casein (CSN1S1) and κ-casein (CSN3) genotypes on milk coagulation properties in Murciano-Granadina goats

The effects of the caprine αs1-casein (CSN1S1) polymorphisms on milk quality and cheese yield have been widely studied in French and Italian goat breeds. Much less is known about the consequences of κ-casein (CSN3) genotype on the technological and coagulation properties of goat milk. In the current study, we have performed an association analysis between polymorphisms at the goat CSN1S1 and CSN3 genes and milk coagulation (rennet coagulation time, curdling rate and curd firmness) and technological (time to cutting of curd and cheese yield) properties. In this analysis, we have included 193 records from 74 Murciano-Granadina goats (with genotypes constituted by different combinations of alleles B, E and F of the gene CSN1S1 and alleles A and B of the gene CSN3) distributed in three herds, which were collected bimonthly during a whole lactation. Data analysis, using a linear mixed model for repeated observations, revealed significant associations between CSN1S1 genotypes and the rate of the curdling process. In this way, milk from EE goats had a significantly higher curdling rate than milk from BB individuals (P<0·05). Contrary to previous experiments performed in French breeds, cheese yield was not significantly different in BB, EE and EF goats. Moreover, we have shown that CSN3 genotype has a significant effect on the rennet coagulation time (BB>AB, P<0·05) but not on cheese yield. No interaction between the CSN1S1 and CSN3 genotypes was observed.

Molecular characterization and phylogenetic analysis of a yak (Bos grunniens) κ-casein cDNA from lactating mammary gland

κ-Casein is one of the major proteins in the milk of mammals. It plays an important role in determining the size and specific function of milk micelles. We have previously identified and characterized a genetic variant of yak κ-casein by evaluating genomic DNA. Here, we isolate and characterize a yak κ-casein cDNA harboring the full-length open reading frame (ORF) from lactating mammary gland. Total RNA was extracted from mammary tissue of lactating female yak, and the κ-casein cDNA were synthesized by RT-PCR technique, then cloned and sequenced. The obtained cDNA of 660-bp contained an ORF sufficient to encode the entire amino acid sequence of κ-casein precursor protein consisting of 190 amino acids with a signal peptide of 21 amino acids. Yak κ-casein has a predicted molecular mass of 19,006.588 Da with a calculated isoelectric point of 7.245. Compared with the corresponding sequences in GenBank of cattle, buffalo, sheep, goat, Arabian camel, horse, and rabbit, yak κ-casein sequence had identity of 64.76-98.78% in cDNA, and identity of 44.79-98.42% and similarity of 53.65-98.42% in deduced amino acids, revealing a high homology with the other livestock species. Based on κ-casein cDNA sequences, the phylogenetic analysis indicated that yak κ-casein had a close relationship with that of cattle. This work might be useful in the genetic engineering researches for yak κ-casein.

Invited review: milk protein polymorphisms in cattle: effect on animal breeding and human nutrition

The 6 main milk proteins in cattle are encoded by highly polymorphic genes characterized by several nonsynonymous and synonymous mutations, with up to 47 protein variants identified. Such an extensive variation was used for linkage analysis with the description of the casein cluster more than 30 yr ago and has been applied to animal breeding for several years. Casein haplotype effects on productive traits have been investigated considering information on the whole casein complex. Moreover, mutations within the noncoding sequences have been shown to affect the specific protein expression and, as a consequence, milk composition and cheesemaking. Milk protein variants are also a useful tool for breed characterization, diversity, and phylogenetic studies. In addition, they are involved in various aspects of human nutrition. First, the occurrence of alleles associated with a reduced content of different caseins might be exploited for the production of milk with particular nutritional qualities; that is, hypoallergenic milk. On the other hand, the frequency of these alleles can be decreased by selection of sires using simple DNA tests, thereby increasing the casein content in milk used for cheesemaking. Furthermore, the biological activity of peptides released from milk protein digestion can be affected by amino acid exchanges or deletions resulting from gene mutations. Finally, the gene-culture coevolution between cattle milk protein genes and human lactase genes, which has been recently highlighted, is impressive proof of the nonrandom occurrence of milk protein genetic variation over the centuries.

Short communication: casein haplotype variability in sicilian dairy goat breeds

In the Mediterranean region, goat milk production is an important economic activity. In the present study, 4 casein genes were genotyped in 5 Sicilian goat breeds to 1) identify casein haplotypes present in the Argentata dell'Etna, Girgentana, Messinese, Derivata di Siria, and Maltese goat breeds; and 2) describe the structure of the Sicilian goat breeds based on casein haplotypes and allele frequencies. In a sample of 540 dairy goats, 67 different haplotypes with frequency >or=0.01 and 27 with frequency >or=0.03 were observed. The most common CSN1S1-CSN2-CSN1S2-CSN3 haplotype for Derivata di Siria and Maltese was FCFB (0.17 and 0.22, respectively), whereas for Argentata dell'Etna, Girgentana and Messinese was ACAB (0.06, 0.23, and 0.10, respectively). According to the haplotype reconstruction, Argentata dell'Etna, Girgentana, and Messinese breeds presented the most favorable haplotype for cheese production, because the casein concentration in milk of these breeds might be greater than that in Derivata di Siria and Maltese breeds. Based on a cluster analysis, the breeds formed 2 main groups: Derivata di Siria, and Maltese in one group, and Argentata dell'Etna and Messinese in the other; the Girgentana breed was between these groups but closer to the latter.

A critical analysis of production-associated DNA polymorphisms in the genes of cattle, goat, sheep, and pig

Increasing productivity is one of the main objectives in animal production. Traditional breeding methods have led to increased gains in some traits but gains are not easily attainable in traits with low heritabilities. Exploiting the genetic variations underlying desired phenotypes is the goal of today's animal producers. Such positive genetic variants must, however, be known before possible application. Consequently, candidate genes of traits of interest have been searched for possible relationships with such traits or to explain reported quantitative trait loci (QTL) for such traits. DNA variants or polymorphisms have been identified in many such genes and their relationships with production traits determined. However, only a few genes have been evaluated, given the wealth of information on reported QTL for production traits, and in most cases genes are only partially investigated. This review presents available information on DNA variants for production traits and discusses steps that are required for effective utilization of this information for successful marker-assisted selection programs.