Focusing on the Goat Casein Complex

Discriminant canonical tool for inferring the effect of αS1, αS2, β, and κ casein haplotypes and haplogroups on zoometric/linear appraisal breeding values in Murciano-Granadina goats

Genomic tools have shown promising results in maximizing breeding outcomes, but their impact has not yet been explored. This study aimed to outline the effect of the individual haplotypes of each component of the casein complex (αS1, β, αS2, and κ-casein) on zoometric/linear appraisal breeding values. A discriminant canonical analysis was performed to study the relationship between the predicted breeding value for 17 zoometric/linear appraisal traits and the aforementioned casein gene haplotypic sequences. The analysis considered a total of 41,323 zoometric/linear appraisal records from 22,727 primiparous does, 17,111 multiparous does, and 1,485 bucks registered in the Murciano-Grandina goat breed herdbook. Results suggest that, although a lack of significant differences (p > 0.05) was reported across the predictive breeding values of zoometric/linear appraisal traits for αS1, αS2, and κ casein, significant differences were found for β casein (p < 0.05). The presence of β casein haplotypic sequences GAGACCCC, GGAACCCC, GGAACCTC, GGAATCTC, GGGACCCC, GGGATCTC, and GGGGCCCC, linked to differential combinations of increased quantities of higher quality milk in terms of its composition, may also be connected to increased zoometric/linear appraisal predicted breeding values. Selection must be performed carefully, given the fact that the consideration of apparently desirable animals that present the haplotypic sequence GGGATCCC in the β casein gene, due to their positive predicted breeding values for certain zoometric/linear appraisal traits such as rear insertion height, bone quality, anterior insertion, udder depth, rear legs side view, and rear legs rear view, may lead to an indirect selection against the other zoometric/linear appraisal traits and in turn lead to an inefficient selection toward an optimal dairy morphological type in Murciano-Granadina goats. Contrastingly, the consideration of animals presenting the GGAACCCC haplotypic sequence involves also considering animals that increase the genetic potential for all zoometric/linear appraisal traits, thus making them recommendable as breeding animals. The relevance of this study relies on the fact that the information derived from these analyses will enhance the selection of breeding individuals, in which a desirable dairy type is indirectly sought, through the haplotypic sequences in the β casein locus, which is not currently routinely considered in the Murciano-Granadina goat breeding program.

Molecular Characterization and Phylogenetic Analysis of Casein Gene Family in Camelus ferus

Camel milk is known for its exceptional medical uses. It has been used since ancient times to treat infant diarrhea, hepatitis, insulin-dependent diabetes (IDDM), lactose intolerance, alcohol-induced liver damage, allergies, and autism. It has the power to treat several diseases, with cancer being the most significant. This study investigated the evolutionary relationship, physiochemical characteristics, and comparative genomic analysis of the casein gene family (CSN1S1, CSN2, CSN1S2, and CSN3) in Camelus ferus. Molecular phylogenetics showing the camelid species clustered casein nucleotide sequences into four groups: CSN1S1, CSN2, CSN1S2, and CSN3. The casein proteins from camels were evaluated and found to be unstable, thermostable, and hydrophilic. CSN1S2, CSN2, and CSN3 were acidic, but CSN1S1 was basic. CSN1S1 showed positive selection for one amino acid (Q), CSN1S2 and CSN2 for three (T, K, Q), and CSN3 showed no positive selection. We also compared high-milk-output species such as cattle (Bos Tarus) and low-milk-yield species such as sheep (Ovies Aries) with camels (Camel ferus) and discovered that YY1 sites are more frequent in sheep than in camels and very low in cattle. We concluded that the ratio of YY1 sites in these species may affect milk production.

A New AS-PCR Method to Detect CSN201 Allele, Genotyping at Ca-Sensitive Caseins Loci and Milk Traits Association Studies in Autochthonous Lazio Goats

Calcium-sensitive caseins are the main protein component of milk. In the goat, they are encoded by three genes (CSN1S1, CSN2, and CSN1S2) located on chromosome 6. A high number of alleles has been discovered for these genes in the goat species, responsible for changes in the milk’s qualitative and quantitative characteristics. This study aimed to develop an Allele-Specific PCR (AS-PCR), which allowed us to unequivocally detect goat carriers of the CSN201 allele. Subsequently, the calcium-sensitive casein loci genotype was investigated in three native goat breeds of the Lazio Region (Bianca Monticellana, Capestrina, and Ciociara Grigia). No individuals were carriers of the CSN1S101, CSN1S1E, CSN201, CSN1S2D, and CSN1S20 alleles, while a high frequency of the alleles CSN1S1F and CSN1S1A*,B* was observed. Association analyses between the different genotypes at the CSN1S1 locus and some milk traits, namely the fat and protein yielded and the fat, protein, solids-not-fat, and casein percentages without an effect on the milk yield, were observed.

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.

Tunisian camel casein gene characterization reveals similarities and differences with Sudanese and Nigerian populations

Milk is a primary protein source that has always played a role in mammalian health. Despite the intensification of research projects on dromedary and the knowledge of the genetic diversity at the casein loci, the genetic structure of the Tunisian camel population still needs exploration. This study sought to determine the genetic diversity of 3 casein gene variants in 5 Tunisian camel ecotypes: c.150G>T at CSN1S1 (α_S1-casein), g.2126A>G at CSN2 (β-casein), and g.1029T>C at CSN3 (κ-casein). The obtained results were compared with data published on Sudanese and Nigerian camels to establish the level of differentiation within and between populations. A total of 159 blood samples were collected from 5 Tunisian camel ecotypes and the extracted DNA was genotyped by PCR-RFLP. A streamlined genotyping protocol was also developed for CSN3. Results indicated that allele T was quite rare (0.06) at CSN1S1 for all ecotypes. Minor allele frequency was found for G (0.462) in CSN2 except for Ardhaoui Medenine ecotype who deviated from the average CSN2 allele frequency of the total population. Allele C showed minor allele frequency of 0.384 in CSN3. Among the Tunisian population, GAT (0.343) was the most represented haplotype in all ecotypes except for Ardhaoui Medenine, where GGC (0.322) was the most frequent one. Significant differences in heterozygosity and local inbreeding were observed across the Tunisian, Sudanese, and Nigerian populations, although the global fixation index indicated that only 2.2% of the genetic variance is related to ecotype differences. Instead, phylogenetic analysis revealed a closer link between the Tunisian and Sudanese populations through a clade subdivision with 3 main branches among the ecotypes. This study represents the first attempt to understand casein gene variability in Tunisian camels; with further study, milk traits and genetic differentiation among populations can be associated with the history of camel domestication.

Genetic Diversity in Casein Gene Cluster in a Dromedary Camel (C. dromedarius) Population from the United Arab Emirates

Genetic polymorphisms, causing variation in casein genes (CSN1S1, CSN1S2, CSN2, and CSN3), have been extensively studied in goats and cows, but there are only few studies reported in camels. Therefore, we aimed to identify alleles with functional roles in the United Arab Emirates dromedary camel (Camelus dromedarius) population to complement previous studies conducted on the same species. Using targeted next-generation sequencing, we sequenced all genes in the casein gene cluster in 93 female camels to identify and characterize novel gene variants. Most variants were found in noncoding introns and upstream sequences, but a few variants showed the possibility of functional impact. CSN2 was found to be most polymorphic, with total 91 different variants, followed by CSN1S1, CSN3 and CSN1S2. CSN1S1, CSN1S2 and CSN2 each had at least two variants while CSN3 had only one functional allele. In future research, the functional impact of these variants should be investigated further.

Authenticity assessment of dairy products by capillary electrophoresis

Milk and derivatives are a very important part in the diet of the world population. Products from goat, buffalo, and sheep species have a greater economic value than the cow ones, therefore, authenticity frauds by improperly adding cow's milk occur frequently: dairy products are among the seven more attractive foods for adulteration. Milk from each of the above-cited animal species has its own definite profile of whey proteins (variants of α-lactalbumin and β-lactoglobulin) and its definite profile of caseins (variants of α_S1 -, α_S2 -, β-, and κ-casein). Such proteins can be usefully exploited as markers of authenticity by using capillary electrophoresis which is the technique of choice for the analysis of proteins. Due to the multiple adjustable parameters that are unknown to other analytical techniques, capillary electrophoresis is able to detect frauds in milk mixtures and cheese with little use of solvents, fast analysis time, and ease of operation. This makes it attractive and competitive for routine checks that are very important to fight the adulteration market. Advantages and limitations are discussed.

Capture Sequencing to Explore and Map Rare Casein Variants in Goats

Genetic variations in the four casein genes CSN1S1, CSN2, CSN1S2, and CSN3 have obtained substantial attention since they affect the milk protein yield, milk composition, cheese processing properties, and digestibility as well as tolerance in human nutrition. Furthermore, milk protein variants are used for breed characterization, biodiversity, and phylogenetic studies. The current study aimed at the identification of casein protein variants in five domestic goat breeds from Sudan (Nubian, Desert, Nilotic, Taggar, and Saanen) and three wild goat species [Capra aegagrus aegagrus (Bezoar ibex), Capra nubiana (Nubian ibex), and Capra ibex (Alpine ibex)]. High-density capture sequencing of 33 goats identified in total 22 non-synonymous and 13 synonymous single nucleotide polymorphisms (SNPs), of which nine non-synonymous and seven synonymous SNPs are new. In the CSN1S1 gene, the new non-synonymous SNP ss7213522403 segregated in Alpine ibex. In the CSN2 gene, the new non-synonymous SNPs ss7213522526, ss7213522558, and ss7213522487 were found exclusively in Nubian and Alpine ibex. In the CSN1S2 gene, the new non-synonymous SNPs ss7213522477, ss7213522549, and ss7213522575 were found in Nubian ibex only. In the CSN3 gene, the non-synonymous SNPs ss7213522604 and ss7213522610 were found in Alpine ibex. The identified DNA sequence variants led to the detection of nine new casein protein variants. New variants were detected for alpha S1 casein in Saanen goats (CSN1S1^∗C1), Bezoar ibex (CSN1S1^∗J), and Alpine ibex (CSN1S1^∗K), for beta and kappa caseins in Alpine ibex (CSN2^∗F and CSN3^∗X), and for alpha S2 casein in all domesticated and wild goats (CSN1S2^∗H), in Nubian and Desert goats (CSN1S2^∗I), or in Nubian ibex only (CSN1S2^∗J and CSN1S2^∗K). The results show that most novel SNPs and protein variants occur in the critically endangered Nubian ibex. This highlights the importance of the preservation of this endangered breed. Furthermore, we suggest validating and further characterizing the new casein protein variants.

Polymorphism of CSN1S1 (g.12164G>A) and CSN2 (g.8913C>A) genes in pure and cross dairy goats

Casein genes directly control milk protein of animals. CSN1S1 (αS1-Casein) and CSN2 (β-Casein) genes influence on milk protein fractions. Genetic polymorphisms of CSN1S1 gene at g.12164G>A locus and CSN2 gene at g.8913C>A locus were identified by PCR-RFLP technique. Animal samples were pure dairy goats providing PE (5 hds.), Saanen (8 hds.) and their crosses providing Sapera (50% Saanen, 50% PE) (51 hds.) and SaanPE (75% Saanen, 25% PE) (3 hds.) from IRIAP dairy goat station. Allele frequency, genotype frequency, heterozygosity value, and Hardy-Weinberg (H-W) equilibrium value were analyzed by Popgen32 program. CSN1S1_g.12164G>A locus resulted in two alleles, i.e. G allele (192 bp, 145 bp, and 101 bp) and A allele (337 bp and 101 bp). The G allele from the highest frequenciest was successively Saanen (0.625), Sapera (0.578), PE (0.400), and SaanPE (0.333). Most dairy goats were heterozygote (Ho>He) and in H-W equilibrium (q2 count < q2P0.05). Whereas CSN2_g.8913C>A locus was monomorphic for possesing only C allele (233 bp and 162 bp), without A allele (416 bp). The existent g.12164G>A SNP of the CSN1S1 gene of could be a potencial molecular selection marker of milk protein content in dairy goat.

Bayesian Analysis of the Association between Casein Complex Haplotype Variants and Milk Yield, Composition, and Curve Shape Parameters in Murciano-Granadina Goats

Considering casein haplotype variants rather than SNPs may maximize the understanding of heritable mechanisms and their implication on the expression of functional traits related to milk production. Effects of casein complex haplotypes on milk yield, milk composition, and curve shape parameters were used using a Bayesian inference for ANOVA. We identified 48 single nucleotide polymorphisms (SNPs) present in the casein complex of 159 unrelated individuals of diverse ancestry, which were organized into 86 haplotypes. The Ali and Schaeffer model was chosen as the best fitting model for milk yield (Kg), protein, fat, dry matter, and lactose (%), while parabolic yield-density was chosen as the best fitting model for somatic cells count (SCC × 103 sc/mL). Peak and persistence for all traits were computed respectively. Statistically significant differences (p < 0.05) were found for milk yield and components. However, no significant difference was found for any curve shape parameter except for protein percentage peak. Those haplotypes for which higher milk yields were reported were the ones that had higher percentages for protein, fat, dry matter, and lactose, while the opposite trend was described by somatic cells counts. Conclusively, casein complex haplotypes can be considered in selection strategies for economically important traits in dairy goats.

Analysis of Candidate Genes for Growth and Milk Performance Traits in the Egyptian Barki Sheep

The most common sheep breeds of Egypt are Ossimi, Rahmani, and Barki breeds. The latter one is well adapted to the challenging desert environment, characterized by food shortage and a high temperature fluctuation. Growth performance of Barki sheep has an important economic value in terms of minimizing the shortage of mutton meat in Egypt. Further, milk production is of great importance for feeding newborn lambs. Eight candidate genes, recently associated with production traits in different breeds, were used to study the effect of genotype on lamb growth and ewe milk traits. The examined genes were LEP, IGF1, DGAT1, STAT5A, PRL, CSN1S2, GHR, and GHRHR, of which one representative single nucleotide polymorphism (SNP) located in the coding region was selected for genotyping. Data from 251 Barki sheep were used in this study. Association analysis between SNPs and lamb growth traits identified rs420693815 of the LEP gene to be significantly associated with weaning weight and average daily gain. In ewes, significant effects on milk yield and composition have been estimated for LEP (rs420693815), STAT5A (rs161082816), PRL (rs422713690), and GHRHR (rs414991449). The results indicated that these genes might be considered as interesting candidates for further investigations to improve growth and milk performance in Barki sheep.

Non-parametric association analysis of additive and dominance effects of casein complex SNPs on milk content and quality in Murciano-Granadina goats

Goat milk casein proteins (αS1, αS2, β and κ) are encoded by four loci (CSN1S1, CSN1S2, CSN2 and CSN3, respectively) clustered within 250 kb in chromosome 6. In this study, 159 Murciano-Granadina goats were genotyped for 48 SNPs within the entire casein region. Phenotypes on milk yield and components were obtained from 2,594 dairy registries. Additive and dominance effects on milk composition and quality were studied using non-parametric tests and principal component analysis to prevent SNPs multicollinearity. Two deletions in exon 4 (CSN1S1 and CSN3), one in exon 7 (CSN2) and one in exon 15 (CSN1S2) have been found at frequencies ranging from 0.12 to 0.50. Bonferroni-corrected significant SNP additive and dominance effects were found for milk yield, fat, protein, dry matter and lactose, and somatic cells. Exons 15 and 7 were significantly associated with milk yield and components except for lactose and somatic cells, while exon 4 was significantly associated with milk yield and components except for protein and dry matter. SNPs' associations with somatic cells were less frequent and weaker than those with milk yield and components. As caseins increase, somatic cells decrease, reducing milk enzymatic activity and consumption suitability. Hence, including molecular information in breeding schemes may promote production efficiency, as selecting against undesirable alleles could prevent the compromises derived from their dominance effects.

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.

Casein Gene Cluster in Camelids: Comparative Genome Analysis and New Findings on Haplotype Variability and Physical Mapping

The structure of casein genes has been fully understood in llamas, whereas in other camelids, this information is still incomplete. In fact, structure and polymorphisms have been identified in three (CSN1S1, αs1-CN; CSN2, β-CN; CSN3, κ-CN) out of four casein genes, whereas controversial information is available for the CSN1S2 (αs2-CN) in terms of structure and genetic diversity. Data from the genome analysis, whose assembly is available for feral camel, Bactrian, dromedary, and alpaca, can contribute to a better knowledge. However, a majority of the scaffolds available in GenBank are still unplaced, and the comparative annotation is often inaccurate or lacking.Therefore, the aims of this study are 1) to perform a comparative genome analysis and synthesize the literature data on camelids casein cluster; 2) to analyze the casein variability in two dromedary populations (Sudanese and Nigerian) using polymorphisms at CSN1S1 (c.150G > T), CSN2 (g.2126A > G), and CSN3 (g.1029T > C); and 3) to physically map the casein cluster in alpaca. Exon structures, gene and intergenic distances, large insertion/deletion events, SNPs, and microsatellites were annotated. In all camelids, the CSN1S2 consists of 17 exons, confirming the structure of llama CSN1S2 gene. The comparative analysis of the complete casein cluster (∼190kb) shows 12,818 polymorphisms. The most polymorphic gene is the CSN1S1 (99 SNPs in Bactrian vs. 248 in dromedary vs. 626 in alpaca). The less polymorphic is the CSN3 in the Bactrian (22 SNPs) and alpaca (301 SNPs), whereas it is the CSN1S2 in dromedary (79 SNPs). In the two investigated dromedary populations, the allele frequencies for the three markers are slightly different: the allele C at CSN1S1 is very rare in Nigerian (0.054) and Sudanese dromedaries (0.094), whereas the frequency of the allele G at CSN2 is almost inverted. Haplotype analysis evidenced GAC as the most frequent (0.288) and TGC as the rarest (0.005). The analysis of R-banding metaphases hybridized with specific probes mapped the casein genes on chromosome 2q21 in alpaca. These data deepen the information on the structure of the casein cluster in camelids and add knowledge on the cytogenetic map and haplotype variability.

Casein haplotype diversity in seven dairy goat breeds

Selection, drift, gene flow and breeding have extensively shaped the genomic variability of domestic animals. In goat species, several mutations identified within the casein genes have been shown to affect the level of gene expression of milk production traits. The four casein genes - CSN1S1, CSN2, CSN1S2 and CSN3 - are organized in a cluster of 250 kb located in chromosome 6, and due to tight linkage, their genetic variability is well depicted by haplotypes which are transmitted to the progeny. Thirty single nucleotide polymorphisms (SNPs) located within the casein gene cluster were used to characterize the haplotype variability of six southern Italian goat breeds (Girgentana, Maltese, Rossa Mediterranea, Argentata dell'Etna, Messinese, Capra dell'Aspromonte). A representative sample of the Norwegian dairy goat breed (Norsk melkegeit) has been used as an out-group to obtain a weighted measure of genetic diversity in the metapopulation. A total of 54 haplotypes were detected among the seven breeds: 26, 9, 8 and 11 haplotypes were found at CSN1S1, CSN2, CSN1S2 and CSN3 respectively. The number of haplotypes per breed was 14 (Norwegian), 26 (Messinese), 27 (Rossa Mediterranea and Girgentana) and 31 (Maltese, Argentata dell'Etna and Capra dell'Aspromonte). The Maltese breed showed the highest number of private haplotypes, whereas the Norwegian goat recorded the highest number of shared haplotypes. The linkage disequilibrium analysis showed higher levels of association for the SNP pairs within casein loci than SNP pairs between casein loci, likely reflecting low levels of intra-genic recombination. The highest linkage disequilibrium values were found in CSN1S1 and CSN2 genes in all the breeds, except for Argentata dell'Etna and Rossa Mediterranea. The resolution of the haplotype diversity at the casein cluster can be exploited both for selective and conservative plans.

Genomic analysis of the origins of extant casein variation in goats

The variation in the casein genes has a major impact on the milk composition of goats. Even though many casein polymorphisms have been identified so far, we do not know yet whether they are evolutionarily ancient (i.e., they existed before domestication) or young (i.e., they emerged after domestication). Herewith, we identified casein polymorphisms in a data set of 106 caprine whole-genome sequences corresponding to bezoars (Capra aegagrus, the ancestor of domestic goats) and 4 domestic goat (Capra hircus) populations from Europe, Africa, the Far East, and the Near East. Domestic and wild goat populations shared a substantial number of casein SNP, from 36.1% (CSN2) to 55.1% (CSN1S2). The comparison of casein variation among bezoars and the 4 domestic goat populations demonstrated that more than 50% of the casein SNP are shared by 2 or more populations, and 18 to 44% are shared by all populations. Moreover, the majority of casein alleles reported in domestic goats also segregate in the bezoar, including several alleles displaying significant associations with milk composition (e.g., the A/B alleles of the CSN1S1 and CSN3 genes, the A allele of the CSN2 gene). We conclude that much of the current diversity of the caprine casein genes comes from ancient standing variation segregating in the ancestor of modern domestic goats.

Population structure, genetic diversity and selection signatures within seven indigenous Pakistani goat populations

Goat farming in Pakistan depends on indigenous breeds that have adapted to specific agro-ecological conditions. Pakistan has a rich resource of goat breeds, and the genetic diversity of these goat breeds is largely unknown. In this study, genetic diversity and population structure were characterized from seven indigenous goat breeds using the goat 50K SNP chip. The genetic diversity analysis showed that Bugi toori goats have the highest inbreeding level, consistent with the highest linkage disequilibrium, lowest diversity and long run of heterozygosity segments. This indicates that this breed should be prioritized in future conservation activities. The population structure analysis revealed four fairly distinct clusters (including Bugi toori, Bari, Black Tapri and some Kamori) and three other breeds that are seemingly the results of admixture between these or related groups (some Kamori, Pateri, Tapri and White Tapri). The selection signatures were evaluated in each breed. A total of 2508 putative selection signals were reported. The 26 significant windows were identified in more than four breeds, and selection signatures spanned several genes that directly or indirectly influence traits included coat colour variation (KIT), reproduction (BMPR1B, GNRHR, INSL6, JAK2 and EGR4), body size (SOCS2), ear size (MSRB3) and milk composition (ABCG2, SPP1, CSN1S2, CSN2, CSN3 and PROLACTIN).

A 100-Year Review: Advances in goat milk research

In the century of research chronicled between 1917 and 2017, dairy goats have gone from simply serving as surrogates to cows to serving as transgenic carriers of human enzymes. Goat milk has been an important part of human nutrition for millennia, in part because of the greater similarity of goat milk to human milk, softer curd formation, higher proportion of small milk fat globules, and different allergenic properties compared with cow milk; however, key nutritional deficiencies limit its suitability for infants. Great attention has been given not only to protein differences between goat and cow milk, but also to fat and enzyme differences, and their effect on the physical and sensory properties of goat milk and milk products. Physiological differences between the species necessitate different techniques for analysis of somatic cell counts, which are naturally higher in goat milk. The high value of goat milk throughout the world has generated a need for a variety of techniques to detect adulteration of goat milk products with cow milk. Advances in all of these areas have been largely documented in the Journal of Dairy Science (JDS), and this review summarizes such advances.

Population genetic structure and milk production traits in Girgentana goat breed

The aim of this work was to evaluate the genetic status of the Girgentana goat, an endangered breed from Sicily (Italy), using microsatellite markers. Furthermore, as the main purpose of the Girgentana breed is milk production, quantitative milk traits were investigated, including fatty acid profile. Molecular data from CSN1S1, CSN2, CSN1S2, and CSN3 casein genes were also used to infer haplotypes. A total of 264 individuals were collected. Samples of Maltese (n = 41) and Derivata di Siria (n = 33) goat breeds were also used to understand the genetic relationship among breeds. Test-day records for milk production were collected to determine daily milk yield, fat, protein, casein, lactose, and somatic cell count. Individual milk samples were also collected for fatty acid extraction. Wright’s statistics, gene flow, Nei genetic distance, factorial correspondence analysis, and Bayesian assignment test showed the existence of genetic variability and differentiation among breeds. The AMOVA results indicated that 89.96% of the total variance was partitioned within populations. The Girgentana breed appears to have a subdivided population, and has not experienced a recent bottleneck. A high variability in milk yield was observed. Mean morning milk yield was 1448 ± 404 g, with 4.30 ± 0.87% and 3.72 ± 0.44% of fat and protein percentages, respectively. The average somatic cell count found in Girgentana goat milk was higher than the threshold of 1 500 000 cells/mL advised in Europe for fresh milk. Gross milk and fatty acid composition were similar to that reported in the literature for other local goat breeds.

Lama glama αS1-casein: Identification of new polymorphisms in the CSN1S1 gene

South American camelids have been poorly genetically investigated and little information is available in llamas (Lama glama) regarding the diversity of the caseins at the protein and gene level. Exon skipping and duplication events previously reported in the α_S1-casein gene (CSN1S1) led us to investigate the genetic variability at this locus. Seventy-two positive clones for the α_S1-casein transcripts were analyzed and randomly sequenced. The comparative analysis of the sequences revealed 2 transitions, c.366A>G and c.690T>C, at the 10th nucleotide of exon 12 and 94 bp of exon 19, respectively. These SNP are responsible for 2 amino acid changes, Ile→Val in position 86 and Tyr→His in position 194 of the mature protein. Both polymorphisms clarify the genetic events behind the protein variants A and B. This result was confirmed by isoelectric focusing analysis of llama milk samples. Quick methods based on PCR-RFLP and allele-specific PCR were set up for allelic discrimination in a population of 128 animals. Based on genotyping results, 4 haplotypes were observed and the estimated frequencies indicated B as the most common haplotype (0.629) in the investigated population. These data add knowledge to the genetic variability of a species little investigated, and open opportunity for new investigation in the field of milk protein for South American camelids, including the possibility, in the future, to select alleles with favorable characteristics.

Molecular Characterization of the Llamas (Lama glama) Casein Cluster Genes Transcripts (CSN1S1, CSN2, CSN1S2, CSN3) and Regulatory Regions

In the present paper, we report for the first time the characterization of llama (Lama glama) caseins at transcriptomic and genetic level. A total of 288 casein clones transcripts were analysed from two lactating llamas. The most represented mRNA populations were those correctly assembled (85.07%) and they encoded for mature proteins of 215, 217, 187 and 162 amino acids respectively for the CSN1S1, CSN2, CSN1S2 and CSN3 genes. The exonic subdivision evidenced a structure made of 21, 9, 17 and 6 exons for the αs1-, β-, αs2- and κ-casein genes respectively. Exon skipping and duplication events were evidenced. Two variants A and B were identified in the αs1-casein gene as result of the alternative out-splicing of the exon 18. An additional exon coding for a novel esapeptide was found to be cryptic in the κ-casein gene, whereas one extra exon was found in the αs2-casein gene by the comparison with the Camelus dromedaries sequence. A total of 28 putative phosphorylated motifs highlighted a complex heterogeneity and a potential variable degree of post-translational modifications. Ninety-six polymorphic sites were found through the comparison of the lama casein cDNAs with the homologous camel sequences, whereas the first description and characterization of the 5'- and 3'-regulatory regions allowed to identify the main putative consensus sequences involved in the casein genes expression, thus opening the way to new investigations -so far- never achieved in this species.

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.

The β-casein in camels: molecular characterization of the CSN2 gene, promoter analysis and genetic variability

The β-casein is the most abundant protein in camel milk and its encoding gene (CSN2) is considered in other species a 'major' gene for the presence of alleles associated to different level of expression. In the present paper, we report for the first time the characterization of the nucleotide sequence of the whole β-casein-encoding gene (CSN2) plus 2,141 bp at the 5'-flanking region in Camelus dromedarius. The promoter region and the complete cDNA are also provided for the first time in Camelus bactrianus. The gene is spread over 7.8 kb and consists of 9 exons varying in length from 24 bp (exon 5) to 519 bp (exon 7), and 8 introns from 95 bp (intron 5) to 1,950 bp (intron 1). The composite response element (CoRE) region was identified in the promoter, whereas the presence of mature microRNA sequences improves the knowledge on the factors putatively involved in the gene regulation. A total of 46 polymorphic sites have been detected. The transition g.2126A>G falls within the TATA-box of dromedary CSN2 promoter with a putative influence on the transcription factor binding activity. The frequency of the G allele is 0.35 in a population of 180 she-camels belonging to 4 different ecotypes. In the same population, a conservative SNP (g.4175C>A) was found at the codon 7 of the signal peptide, whereas a comparative analysis with a cDNA sequence available in the database evidenced a missense SNP (g.4180T(Leu)>G(Arg)) at exon 2. Four SNPs were found in the bactrian camel. The SNP c.666G>A is responsible for the amino acid change Met(201)→Ile and it represents the first missense allele at the β-casein in camels. Finally, five interspersed repeated elements were identified at intronic level, whereas the presence of putative bio-functional peptides belonging to ACE-inhibitor and anti-oxidative families confirms the potential protective role of the camel milk for the human nutrition.

Characterization of β-casein gene in Indian riverine buffalo

The study aimed at characterization of buffalo β-casein gene and its promoter by PCR-SSCP analysis. Complete β-casein exon VII region analysis revealed two SSCP band patterns, with pattern-I representing predominant allele B (85%) present in homozygous (genotype BB) condition and pattern-II representing a rare allele A1 present in heterozygous condition (genotype A1B). Sequencing of two patterns revealed three nucleotide substitutions at codon 68, 151 and 193 of exon VII. The cDNA sequence of buffalo β-casein gene indicated three further nucleotide substitutions between allele A1 and B at codon 10, 39, and 41. Analysis of β-casein proximal promoter region (-350 upstream to +32) revealed four SSCP band patterns. These SSCP patterns corresponded to nucleotide substitutions at seven locations within 382 bp 5' UTR region of β-casein gene. Haplotype analysis suggested pattern-I of exon VII (wild type) was associated with three types of promoters and pattern-II of exon VII (rare type) corresponded to one exclusive type of promoter. The study suggested two haplotypes of exon VII and four haplotypes of promoter for buffalo β-casein.

Influence of fresh forage-based diets and αs₁-casein (CSN1S1) genotype on nutrient intake and productive, metabolic, and hormonal responses in milking goats

Polymorphism at the αS1-casein locus (CSN1S1) in goats influences several milk production traits. Milk from goats carrying strong alleles, which are associated with high αS1-casein (αS1-CN) synthesis, has higher fat and casein contents, longer coagulation time and higher curd firmness than milk from goats with weak alleles linked to low αS1-CN content. Nutrition also affects these milk properties; therefore, it is important to better understand the interaction between dietary characteristics and the CSN1S1 genotype in goats. This study aimed to investigate the effect of fresh forage based diet or energy supplement on feeding behavior, milk production, and metabolic and hormonal parameters of Girgentana goats with different genotypes at CSN1S1 loci. From a group of goats genotyped by PCR at the DNA level, 12 were selected because they had the same genotype for αS2-CN, β-CN, and κ-CN but a different genotype for αS1-CN: 6 were homozygous for strong alleles at the CSN1S1 loci (AA) and 6 were heterozygous for a weak allele (AF). Goats of each genotype were allocated to 3 subgroups and fed 3 diets ad libitum in a 3×3 Latin square design. The diets were sulla (Hedysarum coronarium L.) fresh forage, sulla fresh forage plus 800 g/d of barley meal (SFB), and mixed hay plus 800 g/d of barley meal (MHB). Diet had a stronger effect than CSN1S1 genotype. The SFB diet led to the highest energy intake, dry matter (DM) digestibility, and milk yield. The fresh forage diets (SFF and SFB) increased DM and crude protein (CP) intake, CP digestibility, and milk CN compared with the MHB diet. The diets supplemented with energy (SFB, MHB) reduced milk fat and urea, improved CP utilization for casein synthesis, and limited body fat mobilization, in accordance with a lower level of nonesterified fatty acids and higher levels of glucose and IGF-1. With regard to CSN1S1 genotype, AA goats showed higher CP digestibility and lower free thyroxine hormone and cholesterol levels than AF goats. Significant diet × genotype interactions indicated how AA goats, compared with AF goats, showed higher DM digestibility and milk yield when fed the SFB diet, which had more energy. A reduction in free triiodothyronine hormone occurred in AF goats fed the MHB diet, whereas no differences were observed in AA goats. These results demonstrate how goats with a higher capacity for αS1-CN synthesis exhibit more efficient energy and protein utilization, evident at the digestive level, and better productive responses to high-nutrition diets.

Molecular characterization and genetic variability at κ-casein gene (CSN3) in camels

κ-casein is a glycosilated protein belonging to a family of phosphoproteins (αs1, β, αs2, κ) that represent the major protein component in mammalian milk. κ-casein plays an essential role in the casein micelle stabilization, determining the size and the specific function. In the present paper, we report for the first time the characterization of the nucleotide sequence of the whole κ-casein-encoding gene (CSN3) plus 1045 nucleotides at the 5' flanking region in Camelus dromedarius. The promoter region and the complete cDNA were also provided for the first time in Camelus bactrianus. The gene is spread over 9.3kb and consists of 5 exons varying in length from 33bp (exon 3) to 494bp (exon 4), and 4 introns from 1200bp (intron 3) to 2928bp (intron 2). Highly conserved sequences, located in the 5' flanking region, have been found. The regulatory regions of camels seems to be more related to equids than to other compared species. 17 polymorphic sites have been detected, one of these (g.1029T>C) is responsible for the creation of a new putative consensus sequence for the transcription factor HNF-1. In general, these SNPs are the first reported in camels for casein loci. Finally, seven interspersed repeated elements were also identified at intronic level.

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.

The Garfagnina goat: a zootechnical overview of a local dairy population

Domestic livestock with a limited distribution are increasingly recognized in the action plans of the European Union as a reason for protecting rural land. The preservation and enhancement of the native germplasm and traits selected through the ages in different areas of farming is the first step in increasing typical products at a time when high quality products are increasingly in demand. This is the first time that a zootechnical overview has been performed on the Italian native goat population named "Garfagnina," which is registered on the Tuscan regional repertory of genetic resources at risk of extinction. The aim of the study was to give added value to this population by focusing on particular traits that could be used for promoting typical products. Data on the size of the local goats, zoometric measures, breeding system, milk quality, and genetic polymorphisms were collected to get insight into the current state of the population of this type of goat. The native goat population is reared in Tuscany in central Italy, mostly for its milk. The local goat farms considered in our study are located in the hills and mountains of the northwestern Tuscan Apennine area. For every farm we measured at least 10% of the reproductive females (273), randomly chosen, and all reproductive males (47) for a total of 320 subjects. Regarding the management of the animals and the feeding system, semi-extensive farming is practiced in all the flocks. From a morphological point of view the animals are relatively homogeneous, especially in terms of zoometric data, whereas they show a wider variability regarding coat. Milk gross and fatty acid composition were similar to that reported in the literature for bulk goat milk. Moreover, the average of somatic cell count and standard plate count found in Garfagnina goat milk indicated good hygienic farm management and correct milking practices, although milking is mainly manual. The average number of globules per milliliter found in Garfagnina goat milk was almost double compared with the literature, whereas the average diameter was lower. Milk coagulation properties were scarce, thus indicating poor cheesemaking aptitude of Garfagnina milk. Selecting haplotypes carrying alleles associated with a higher expression of the specific casein could help improve milk cheesemaking aptitude. Moreover, the rather high frequency of the faint CSN1S1*F allele and the occurrence of CSN2*0 might suggest that Garfagnina goat milk could be used, after an appropriate selection, for direct consumption of milk at low casein content for intolerant human subjects.

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.

Characterization of casein gene complex and genetic diversity analysis in Indian goats

Milk protein polymorphism plays an important role in genetic diversity analysis, phylogenetic studies, establishing geographical diversity, conservation decision, and improving breeding goals. Milk protein polymorphism in Indian goat breeds has not been well studied; therefore, an investigation was carried out to analyze the genetic structure of the casein gene and milk protein diversity at six milk protein loci in nine Indian goat breeds/genetic groups from varied agro-climatic zones. Milk protein genotyping was carried out in 1098 individual milk samples by SDS-PAGE at alphaS1-CN (CSN1S1), beta-CN (CSN2), alphaS2-CN (CSN1S2), kappa-CN (CSN3), beta-LG, and alpha-LA loci. Indian goats exhibited alphaS1-casein A allele in higher frequency in the majority of breeds except Ganjam and local goats. The alphaS1-casein A allele frequencies varied from 0.45 to 0.77. A total of 16 casein haplotypes were observed in seven breeds and breed specific haplotypes were observed with respect to geographic region. The average number of alleles was lowest in Ganjam (1.66 +/- 0.81) and highest in Sirohi goats (2.50 +/- 1.05). Expected heterozygosity at six different loci demonstrated genetic diversity and breed fragmentation. Neighbor-Joining tree was built basing on Nei's distance. There was about 16.95% variability due to differences between breeds, indicating a strong subdivision. Principal component analysis was carried out to highlight the relationship among breeds. The variability among goat breeds was contributed by alphaS2-CN, beta-LG and alphaS1-CN. The Indian goats exhibited alphaS1-CN (CSN1S1) A allele in higher frequency in all the breeds indicating the higher casein yield in their milk.

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.