Short communication: casein haplotype variability in sicilian dairy goat breeds

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.

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.

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.

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 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.

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.

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.

Effect of CSN1S1 gene polymorphism and stage of lactation on milk yield and composition of extensively reared goats

The effect of CSN1S1 genotype and lactation stage on milk yield and composition were investigated in 80 extensively reared goats. Milk yield was recorded in early, mid and late lactation and individual milk samples were collected to determine: fat, protein, lactose and casein content, pH, freezing point, somatic cell count (SCC) and total microbic mesophilic count (TMC). Relative casein composition and amino acid profile were quantified by HPLC. Fatty acid profile was measured by gas-chromatography. Genotype did not affect milk yield, while this trait was significantly affected by lactation stage (P < 0·01). CSN1S1 BB goats produced significantly higher protein and casein percentages (P < 0·05). αs1-casein (CN) was significantly higher in BB and AB goats than AF and BF, showing intermediate values in AA goats (P < 0·01). The protein percentage and the αs1 and αs2-CN fractions were not affected by lactation stage, while the casein content and the β and κ-CN significantly increased throughout lactation (P < 0·01). C4 : 0 and C6 : 0 were not affected by genotype, while C8 : 0 and C10 : 0 were higher in the AA goats than BB; most of the long chain FA were higher in BB than AA goats. MUFA and PUFA increased in late lactation. In addition, BB goats showed higher essential amino acids, resulting in an optimal composition from the nutritional point of view, when compared with AA goats. The increase of MUFA, PUFA, essential and cis-FA in late lactation indicate that the lipid composition of goat's milk, with the progress of lactation, tends to improve its nutritional value.

Genetic polymorphism at the CSN1S1 gene in Girgentana dairy goat breed

The aim of this work was to evaluate the variability of the αs1-casein locus in the endangered Girgentana dairy goat breed in order to define genetic improvement and a conservation program for this breed. The study was performed on 200 dairy goats by means of different PCR protocols. The most frequent alleles were A (0.590) and F (0.290) followed by B (0.065) and N (0.047). CSN1S1 E allele was identified with a very low frequency (0.008). The most common genotype was AF (0.365) followed by AA (0.340). The high frequency of the strong genotypes is associated with the production of milk with high fat and protein content and with optimal technological properties. In Girgentana goat breed, the CSN1S1 genotype information could be utilised in selection strategies for milk protein content and milk yield, in order to select genetic lines for the production of ‘drinking milk’ using weak and null genotypes, and for niche products using strong genotypes.

An initial comparative map of copy number variations in the goat (Capra hircus) genome

BACKGROUND

The goat (Capra hircus) represents one of the most important farm animal species. It is reared in all continents with an estimated world population of about 800 million of animals. Despite its importance, studies on the goat genome are still in their infancy compared to those in other farm animal species. Comparative mapping between cattle and goat showed only a few rearrangements in agreement with the similarity of chromosome banding. We carried out a cross species cattle-goat array comparative genome hybridization (aCGH) experiment in order to identify copy number variations (CNVs) in the goat genome analysing animals of different breeds (Saanen, Camosciata delle Alpi, Girgentana, and Murciano-Granadina) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome.

RESULTS

We identified a total of 161 CNVs (an average of 17.9 CNVs per goat), with the largest number in the Saanen breed and the lowest in the Camosciata delle Alpi goat. By aggregating overlapping CNVs identified in different animals we determined CNV regions (CNVRs): on the whole, we identified 127 CNVRs covering about 11.47 Mb of the virtual goat genome referred to the bovine genome (0.435% of the latter genome). These 127 CNVRs included 86 loss and 41 gain and ranged from about 24 kb to about 1.07 Mb with a mean and median equal to 90,292 bp and 49,530 bp, respectively. To evaluate whether the identified goat CNVRs overlap with those reported in the cattle genome, we compared our results with those obtained in four independent cattle experiments. Overlapping between goat and cattle CNVRs was highly significant (P < 0.0001) suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Genes with environmental functions were over-represented in goat CNVRs as reported in other mammals.

CONCLUSIONS

We describe a first map of goat CNVRs. This provides information on a comparative basis with the cattle genome by identifying putative recurrent interspecies CNVs between these two ruminant species. Several goat CNVs affect genes with important biological functions. Further studies are needed to evaluate the functional relevance of these CNVs and their effects on behavior, production, and disease resistance traits in goats.

Divergence at the casein haplotypes in dairy and meat goat breeds

Casein genes have been proved to have an influence on milk properties, and are in addition appropriate for phylogeny studies. A large number of casein polymorphisms exist in goats, making their analysis quite complex. The four casein loci were analyzed by molecular techniques for genetic polymorphism detection in the two dairy goat breeds Bunte Deutsche Edelziege (BDE; n=96), Weiße Deutsche Edelziege (WDE; n=91), and the meat goat breed Buren (n=75). Of the 35 analyzed alleles, 18 were found in BDE, and 17 in Buren goats and WDE. In addition, a new allele was identified at the CSN1S1 locus in the BDE, showing a frequency of 0·05. This variant, named CSN1S1*A′, is characterized by a t→c transversion in intron 9. Linkage disequilibrium was found at the casein haplotype in all three breeds. A total of 30 haplotypes showed frequencies higher than 0·01. In the Buren breed only one haplotype showed a frequency higher than 0·1. The ancestral haplotype B-A-A-B (in the order: CSN1S1-CSN2-CSN1S2-CSN3) occurred in all three breeds, showing a very high frequency (>0·8) in the Buren.