DNA Sequence Variants and Protein Haplotypes of Casein Genes in German Black Pied Cattle (DSN)

Validation of Bos taurus SNPs for Milk Productivity of Sahiwal Breed (Bos indicus), Pakistan

The aim of the present study was the validation of the already reported Bos taurus SNPs in the Sahiwal breed. A total of nine SNPs of the casein gene were studied. Out of nine, seven Bos taurus SNPs of casein protein genes were found to be significantly associated with milk productivity traits. The genomic DNA was extracted from the mammary alveolar endothelial cells of a flock of 80 purebred Sahiwal lactating dams available at Khizrabad Farm near Sargodha. New allele-specific primers were designed from the NCBI annotated sequence database of Bos taurus to obtain 100 nt-long PCR products. Each dam was tested separately for all the SNPs investigated. Animals with genotype GG for the SNPs rs43703010, rs10500451, and 110323127, respectively, exhibited high milk yield. Similarly, animals with genotype AA for the SNPs rs11079521, rs43703016, and rs43703017 showed high milk yield consistently. For the SNP rs43703015, animals with genotype CC showed high milk productivity. These above-mentioned SNPs have previously been reported to significantly up-regulate casein protein contents in Bos taurus. Our results indicated SNPs that significantly affect the milk protein contents may also significantly increase per capita milk yield. These finding suggest that the above-mentioned reported SNPs can also be used as genetic markers of milk productivity in Sahiwal cattle.

Genetic variants of CSN1S1, CSN2, CSN3, and BLG genes and their association with dairy production traits in Sahiwal cattle and Nili-Ravi buffaloes

Milk protein genes are associated with milk yield and composition in dairy animals. The present study aimed to identify milk protein genes (CSN1S1, CSN2, CSN3, and BLG) genetic variants and their association with milk yield in Sahiwal cattle and Nili-Ravi buffaloes. One hundred animals from each species were selected to collect blood samples and milk production records. Primers were designed for these milk protein genes for PCR amplification. Sequencing of resultant PCR products revealed a higher number of SNPs (13 vs. 7, 5 vs. 1, and 6 vs. 2) in Sahiwal as compared to Nili-Ravi animals in CSN1S1, CSN2, and CSN3 genes, respectively. However, a single SNP was observed in BLG gene of both species. Association analysis revealed that one SNP in BLG gene of Nili-Ravi was associated (p < 0.05) with 305-day milk yield. Two SNPs at CSN1S1 gene in Sahiwal were associated with dry-period. Similarly, one SNP at CSN1S1 and two SNPs at CSN3 gene showed significant association (p < 0.05) with average calving-interval in Sahiwal while two SNPs in CSN1S1 gene were associated (p < 0.05) with this trait in Nili-Ravi. These SNPs could be helpful as candidate variants for marker-assisted selection in cattle and buffaloes for improvement of lactation performance.

Detection of A2A2 genotype of beta casein protein (CSN2) gene in local, exotic and cross bred cattle in Pakistan

Genetic variants of bovine Beta-casein protein (CSN2) gene especially A1 and A2 are the most important variants in dairy cattle. A1 milk protein is considered as risk factor for different disease and milk intolerance which release Beta-Casomorphin-7 during digestion which is a bioactive opioid but not released from A2 milk protein. This opioid is responsible for several human health problems like Coronary Heart disease, type 1 diabetics, milk intolerance and other neurological disorders. In present study, 360 blood sample were collected from Lohani, Achai, jersey, Holstein Friesian, Achai x jersey, Friesian x Sahiwal and Sahiwal x Friesian from different region of Khyber Pakhtunkhwa (KP) province. The polymerase chain reaction (PCR) amplicons were sequenced for the identification of polymorphism in exon 7 of Beta-casein protein (CSN2) gene. Sequencing analysis explored CSN2 genotype in exon 7 using the Genomic sequence from GenBank (X.71104) g.8101 C > A at codon 67. The allelic and genotypic frequencies of CSN2 gene were analyzed and observed that Holstein Friesian cattle exhibited A1A2 33%, A1A1 50% and A2A2 17%, Jersey cattle show 68% A1A1, 18% A1A2 and 14% A2A2, Sahiwal x Friesian 56% A1A1, 26% A1A2 and 18% A2A2, Jersey × Achai 78% A2A2, 15% A1A2 and 7% A1A1, Achai 100% A2A2 Lohani 100% A2A2. This is a preliminary study, conducted with meager resources, therefore, it is very difficult to make conclusion that which particular breed possess harmful alleles and which breed possess useful alleles of beta-casein gene. Therefore, a comprehensive molecular work is needed to be performed with greater number of samples sequencing.

Caseins: Versatility of Their Micellar Organization in Relation to the Functional and Nutritional Properties of Milk

The milk of mammals is a complex fluid mixture of various proteins, minerals, lipids, and other micronutrients that play a critical role in providing nutrition and immunity to newborns. Casein proteins together with calcium phosphate form large colloidal particles, called casein micelles. Caseins and their micelles have received great scientific interest, but their versatility and role in the functional and nutritional properties of milk from different animal species are not fully understood. Caseins belong to a class of proteins that exhibit open and flexible conformations. Here, we discuss the key features that maintain the structures of the protein sequences in four selected animal species: cow, camel, human, and African elephant. The primary sequences of these proteins and their posttranslational modifications (phosphorylation and glycosylation) that determine their secondary structures have distinctively evolved in these different animal species, leading to differences in their structural, functional, and nutritional properties. The variability in the structures of milk caseins influence the properties of their dairy products, such as cheese and yogurt, as well as their digestibility and allergic properties. Such differences are beneficial to the development of different functionally improved casein molecules with variable biological and industrial utilities.

Identification and annotation of milk associated genes from milk somatic cells using expression and RNA-seq data

It is of interest to identify and annotate milk associated genes using expression profiling and RNA-Seq data from milk somatic cells. RNA-Seq data was pre-processed and mapping was done to identify differentially expressed genes (DEG). The functional insights about the up and down regulated genes were gleaned using the protein-protein interaction Network in the STRING database followed by CytoHubba analysis in Cytoscope. Gene ontology, annotation and pathway enrichment was completed using ShinyGO, David tool and QTL analysis. These analysis shows that 21 genes are linked with the secretion of milk.

Demographic pattern of A1/A2 beta casein variants indicates conservation of A2 type haplotype across native cattle breeds (Bos indicus) of India

Genetic variations of the beta casein gene hold importance because of their probable association with human health. Comparative sequence analysis of β-casein gene across Indian native, crossbred and exotic breeds in India revealed 15 SNPs and 4 INDELs corresponding to 14 haplotypes. The frequency of A2 type haplotype was maximum (0.941) across all Indian native breeds. Among the 15 variants reported for taurine breeds, only three (A1, A2 and B) were observed in analysed populations. Allelic profiling of A1/A2 β-casein variants in ~ 4000 animals belonging to three cattle types and breeding bulls also revealed the predominance of A2 allele (0.95) in Indian cattle. The high proportion of A2 allele/haplotype indicates that Indian native cattle are the best suited to meet the demands for A2 milk globally. However, a higher percentage of heterozygous genotype (A1A2) in breeding bulls warrants the need to screen sire lines so as to drift the herd towards A2.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03232-0.

Challenging Sustainable and Innovative Technologies in Cheese Production: A Review

It is well known that cheese yield and quality are affected by animal genetics, milk quality (chemical, physical, and microbiological), production technology, and the type of rennet and dairy cultures used in production. Major differences in the same type of cheese (i.e., hard cheese) are caused by the rennet and dairy cultures, which affect the ripening process. This review aims to explore current technological advancements in animal genetics, methods for the isolation and production of rennet and dairy cultures, along with possible applications of microencapsulation in rennet and dairy culture production, as well as the challenge posed to current dairy technologies by the preservation of biodiversity. Based on the reviewed scientific literature, it can be concluded that innovative approaches and the described techniques can significantly improve cheese production.

Design and performance of a bovine 200 k SNP chip developed for endangered German Black Pied cattle (DSN)

BACKGROUND

German Black Pied cattle (DSN) are an endangered dual-purpose breed which was largely replaced by Holstein cattle due to their lower milk yield. DSN cattle are kept as a genetic reserve with a current herd size of around 2500 animals. The ability to track sequence variants specific to DSN could help to support the conservation of DSN's genetic diversity and to provide avenues for genetic improvement.

RESULTS

Whole-genome sequencing data of 304 DSN cattle were used to design a customized DSN200k SNP chip harboring 182,154 variants (173,569 SNPs and 8585 indels) based on ten selection categories. We included variants of interest to DSN such as DSN unique variants and variants from previous association studies in DSN, but also variants of general interest such as variants with predicted consequences of high, moderate, or low impact on the transcripts and SNPs from the Illumina BovineSNP50 BeadChip. Further, the selection of variants based on haplotype blocks ensured that the whole-genome was uniformly covered with an average variant distance of 14.4 kb on autosomes. Using 300 DSN and 162 animals from other cattle breeds including Holstein, endangered local cattle populations, and also a Bos indicus breed, performance of the SNP chip was evaluated. Altogether, 171,978 (94.31%) of the variants were successfully called in at least one of the analyzed breeds. In DSN, the number of successfully called variants was 166,563 (91.44%) while 156,684 (86.02%) were segregating at a minor allele frequency > 1%. The concordance rate between technical replicates was 99.83 ± 0.19%.

CONCLUSION

The DSN200k SNP chip was proved useful for DSN and other Bos taurus as well as one Bos indicus breed. It is suitable for genetic diversity management and marker-assisted selection of DSN animals. Moreover, variants that were segregating in other breeds can be used for the design of breed-specific customized SNP chips. This will be of great value in the application of conservation programs for endangered local populations in the future.

First DNA Sequencing in Beninese Indigenous Cattle Breeds Captures New Milk Protein Variants

This study investigated polymorphisms in the milk protein genes CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB, and casein haplotypes in Beninese indigenous cattle. Considering 67 animals, DNA sequencing of the genes' exons, flanking regions and parts of the 5'-upstream regions identified 1058 genetic variants including 731 previously unknown. In addition, four novel milk protein variants were detected, including CSN3K (p.Ala66Val), LALBAF (p.Arg58Trp), LGBB1 (p.Ala134Val) and LGBK (p.Thr92Asnfs*13). CSN3K is caused by a novel SNP (BTA6:85656526C>T, exon 4) whereas LALBAF and LGBB1 are due to rs714688595C>T (exon 1) and rs109625649C>T (exon 4), respectively. Regarding LGBK, a frameshift insertion of one adenine residue at BTA11:103257980 (exon 3) induces a premature translation termination resulting in a 46% reduction of the reference protein sequence. The casein polymorphisms formed five main CSN1S1-CSN2-CSN1S2-CSN3 haplotypes including B-A1-A-B, B-A1-A-A and C-A2-A-B which are predominant in the investigated cattle breeds. Moreover, in silico analyses of polymorphisms within the 5'- and 3'- untranslated regions of all six milk proteins revealed effects on microRNA and transcription factor binding sites. This study suggests a large genetic variation of milk protein genes in Beninese cattle, which should be investigated in further studies for their effects on milk production, including quality and yield traits.

Casein (CSN) gene variants and parity affect the milk protein traits in crossbred (Bos taurus x Bos indicus) cows in sub-tropical climate

Genotypes at four casein (CSN) loci-A26181G of CSN1S1, C6227T of CSN1S2, A8101C of CSN2, and A13104C of CSN3-along with non-genetic factors were studied for their effects on various milk protein traits in 100 crossbred cows with major inheritance of Holstein Frisian (Bos taurus) and Tharparkar (Bos indicus). Results demonstrated the presence of all CSN genotypes with a predominance of heterozygotes. At CSN2 (A8101C; His67Pro) locus, the A2 allele, desirable for human health, was present in 62% as heterozygous and 29% in homozygous condition. Among non-genetic factors, parity of the cows had a significant influence on the milk protein traits in these crossbreds. The genotypes at CSN1S1, CSN2, and CSN3 loci were found to influence (p<0.05 to 0.01) the casein and whey protein yields and composition traits. The casein index and total milk yield were most influenced by the CSN1S2 locus. The AA (A1 milk) genotype of CSN2 had significantly higher yields and percentages of casein and whey proteins. Positive influence of CC genotype of CSNS3 on milk proteins of was observed similar to Bos taurus cows; however, such influence of AA genotype of CSN2 locus may be distinctive to the crossbred cows maintained in subtropical condition. Overall, the results revealed the diverse effects of CSN genotypes on milk proteins in crossbred cattle.

Genetic diversity and population structure in divergent German cattle selection lines on the basis of milk protein polymorphisms

The aim of this study was to analyze the genetic structure of the casein cluster in eight selection lines of the Holstein Friesian (HF), German Simmental (GS) and German Black Pied cattle (“Deutsches Schwarzbuntes Niederungsrind”, DSN) breeds. A total of 2962 milk samples were typed at αs1-casein (αs1-CN), β-casein (β-CN), αs2-casein (αs2-CN) and κ-casein (κ-CN) loci using isoelectric focusing. The number of alleles per locus ranged from one (αs2-CN) to five (β-CN), and the average expected heterozygosity and polymorphic information content of all loci were 0.33 and 0.27, respectively. The unrooted dendrogram revealed that the selection lines of the endangered DSN breed were clearly separated from the HF and GS breeds due to their predominance of the β-CN A1 allele and the comprehensive haplotype BA1A (in the abbreviation of αs1-β-κ-CN). Temporal changes in allele distributions indicated decreasing genetic diversity at the casein loci, explaining the moderate level of genetic differentiation among selection lines (7.1 %). The variability of the casein should be exploited in future using breeding programs to select genetic lines for specific protein production in bovine milk but also to preserve biodiversity.

Genomic Loci Affecting Milk Production in German Black Pied Cattle (DSN)

German Black Pied cattle (DSN) is an endangered population of about 2,550 dual-purpose cattle in Germany. Having a milk yield of about 2,500 kg less than the predominant dairy breed Holstein, the preservation of DSN is supported by the German government and the EU. The identification of the genomic loci affecting milk production in DSN can provide a basis for selection decisions for genetic improvement of DSN in order to increase market chances through the improvement of milk yield. A genome-wide association analysis of 30 milk traits was conducted in different lactation periods and numbers. Association using multiple linear regression models in R was performed on 1,490 DSN cattle genotyped with BovineSNP50 SNP-chip. 41 significant and 20 suggestive SNPs affecting milk production traits in DSN were identified, as well as 15 additional SNPs for protein content which are less reliable due to high inflation. The most significant effects on milk yield in DSN were detected on chromosomes 1, 6, and 20. The region on chromosome 6 was located nearby the casein gene cluster and the corresponding haplotype overlapped the CSN3 gene (casein kappa). Associations for fat and protein yield and content were also detected. High correlation between traits of the same lactation period or number led to some SNPs being significant for multiple investigated traits. Half of all identified SNPs have been reported in other studies, previously. 15 SNPs were associated with the same traits in other breeds. The other associated SNPs have been reported previously for traits such as exterior, health, meat and carcass, production, and reproduction traits. No association could be detected between DGAT1 and other known milk genes with milk production traits despite the close relationship between DSN and Holstein. The results of this study confirmed that many SNPs identified in other breeds as associated with milk traits also affect milk traits in dual-purpose DSN cattle and can be used for further genetic analysis to identify genes and causal variants that affect milk production in DSN cattle.

Comparative Effects of Milk Containing A1 versus A2 β-Casein on Health, Growth and β-Casomorphin-7 Level in Plasma of Neonatal Dairy Calves

Simple Summary

Bovine milk generally contains two types of the milk protein β-casein, A1 and A2. Enzymatic digestion of the A1 type yields the opioid peptide β-casomorphin-7, which is suggested to adversely affect human and animal health. This study aimed to compare the effects of milk containing either homozygote A1 or A2 β-casein on health and growth parameters in 47 dairy calves during the first three weeks of life. Additionally, we studied, for the first time, the levels of intact β-casomorphin-7 in plasma of calves fed milk of alternative β-casein genotypes. Milk feeding of “A2-milk” led to a lower milk intake and a looser fecal consistency (higher prevalence of diarrhea) compared to “A1-milk”. Nevertheless, weight gains and end weights of calves of both feeding groups were similar, which might be caused by the associated higher protein content of milk containing the A2 variant. Intact β-casomorphin-7 was detected in plasma after A1- and A2-milk consumption, but was almost 5 times higher for A1-calves. In summary, A2-milk minimized the cleavage of the opioid peptide β-casomorphin-7 and might have advantages in the development of pre-weaned dairy calves.

Abstract

Research has shown that digestion of A1 β-casein (β-CN) affects gastrointestinal motility and opioid activity through the release of the peptide β-casomorphin-7 (β-CM7). In the case of the A2 variant, the cleavage of β-CM7 does not occur or occurs at a very low rate. Therefore, the aim of the study was to compare the effects of milk containing either homozygote A1 or A2 β-CN on health and growth parameters of dairy calves. Forty-seven neonatal calves (24 females, 23 males) of the breeds German Holstein (GH, n = 9), German Simmental (GS, n = 33) and their crossing (GH × GS, n = 5) were used in a 21-day feeding study. Fecal score (FS), respiratory frequency (RF), and rectal body temperature (BT) were recorded daily, whereas body weight was measured at birth and at day 21 to estimate the average daily weight gain (ADG). Additionally, blood was collected from calves three times during the experimental period and, for the first time, the respective plasma samples were analyzed for intact β-CM7. Consumption of A2-milk led to a lower daily milk intake (dMI) (p < 0.05). Furthermore, fecal consistency was softer for calves fed A2-milk (p < 0.05). Although 44% of A2-calves had diarrhea or revealed a tendency towards it (FS ≥ 3), A1-calves had a prevalence of 21%. Calves with a FS of 4 were offered an electrolyte solution and received a dietary food supplement for the stabilization of the fluid and electrolyte balance. Nevertheless, similar ADG and end weights (EW) of calves fed A1- or A2-milk (p > 0.05) indicate that A2-milk may compensate higher diarrhea rates and lower dMI due to the associated higher protein content. This is the first report of intact β-CM7 in plasma of calves fed milk of either A1 or A2 β-CN. Evidence from this study suggests that due to the change in the amino-acid sequence, A2-milk might be able to prevent or, at least, to minimize the cleavage of β-CM7 in calves.

Associations between maternal milk protein genotypes with preweaning calf growth traits in beef cattle

The objectives of this study were to investigate milk casein polymorphisms in dams and to determine the impacts of maternal casein genotypes on growth traits of their sucking calves. Milk samples from 433 dams of the breeds German Angus (GA) and German Simmental (GS) were typed at the milk protein loci α s1-casein (αs1-CN), β-casein (β-CN), α s2-casein (αs2-CN), and κ-casein (κ-CN) via isoelectric focusing. Associations between casein genotypes in maternal milk with growth traits of their 1,872 calves were analyzed until the age of weaning using linear mixed models, considering either genotypes of individual casein loci (model 1) or composite α s1-β-α s2-κ-CN genotypes within the casein cluster (model 2). Besides environmental effects such as sex, age of the dam, and calving year-season, genetic effects (breed group and maternal and paternal effects) were considered in statistical models. The composite casein genotype BBǀA2A2ǀAAǀAB (order of genes on bovine chromosome 6: α s1-ǀβ-ǀα s2-ǀκ-CN) was associated with greater average daily weight gains (ADG) and heavier age-adjusted weaning weights (WW) of calves (P < 0.05). The effects of composite genotypes on birth weight of calves were similar (P > 0.05; model 2). With regard to individual casein loci, greater ADG and WW were observed for calves from dams with the genotypes κ-CN BB and α s1-CN BB, respectively (P < 0.05; model 1). Age-adjusted WW was largest for calves from dams carrying the κ-CN genotype BB (215 kg) compared with calves representing the maternal AB and AA genotypes (both 204 kg). Results from the present study suggested selectable casein genotypes due to their nutritional value of milk (value in terms of offspring performances), offering new perspectives for breeding strategies in beef cattle to improve preweaning calf performance.

Frequencies of milk protein variants and haplotypes estimated from genotypes of more than 1 million bulls and cows of 12 French cattle breeds

Due to their major effects on milk composition and cheese-making properties and their putative effects on human health, there is a great deal of interest in bovine milk protein variants. The objectives of this study were to estimate frequencies of milk protein variants and haplotypes in 12 cattle breeds as well as their trends over time to assess the effect of selection on milk traits. Milk protein variants and haplotypes were identified from SNP genotype data from more than 1 million animals from 12 dairy, beef, or dual-purpose cattle breeds that had been genotyped for genomic selection. We examined a total of 15 loci in the genes that encode β-lactoglobulin (β-LG) and 3 caseins (α_S1-CN, β-CN, and κ-CN); genotypes were directly called from customized SNP chips (50.6%) or imputed (49.4%). Variants A and B of β-LG were frequent in the 12 breeds. For the caseins, we found 3 variants for α_S1-CN (B, C, and D), 6 for β-CN (A1, A2, A3, B, C, and I), and 5 for κ-CN (A, B, C, D, and E). For α_S1-CN, the B variant was the most frequent in all breeds except Jersey. For β-CN, the A2 variant was the most abundant in all breeds except Tarentaise, although in Normande animals, the I variant (30.9%) was almost as common as A2 (39.7%). The C variant was very rare except in the Tarentaise sample (4.8%). The most frequent variant for κ-CN was A in 5 breeds (including Holstein), and B in the 7 other breeds. The B variant was present at a particularly high frequency in Jersey (82.6%) and Normande (85.5%) animals. The C and E variants of κ-CN appeared to be particularly frequent in the Tarentaise (12.7%) and Holstein (9%) breeds, respectively. We found 20 haplotype combinations of α_S1-β-κ CN that were present at a frequency >0.1% in at least one breed; however, only 6 to 9 haplotypes were found in any given breed, demonstrating a strong degree of linkage disequilibrium. The most frequent haplotypes were B-A1-A, B-A2-A, B-A2-B, B-I-B, C-A2-A, and C-A2-B. Some alleles were predominantly found in only one haplotype, such as the E and C variants of κ-CN and the I variant of β-CN, which were mainly found in the B-A1-E, B-A1-C, and B-I-B haplotypes, respectively. We observed changes in the frequency of certain variants over time in several breeds, such as an increase in the frequency of variants A of β-LG, I of β-CN, and B of κ-CN. With these results, we update and complete frequency data that were first estimated 30 to 50 yr ago, and, for the first time in these breeds, we assess the effect of selection on milk protein variants.

Influence of the Casein Composite Genotype on Milk Quality and Coagulation Properties in the Endangered Agerolese Cattle Breed

The aim of this study was the characterization of CSN1S1, CSN2 and CSN3 genetic variability in Agerolese cattle, and the investigation of the effect of casein composite genotypes (CSN1S1, CSN2 and CSN3) on quality and coagulation traits of the corresponding milk. To these purposes, blood and milk from 84 cows were sampled and analysed. Allele frequencies at CSN2 and CSN3 revealed no Hardy-Weinberg equilibrium in the population with a prevalence of allele A² for CSN2 and allele B for CSN3. BBA¹A²AB and BBA²A²AB composite genotypes were the most common in the population. BBA¹A²AB showed a higher total solids and fat content (12.70 ± 0.16 and 3.93 ± 0.10, respectively), while BBA²A²BB showed the best coagulation properties (RCT 12.62 ± 0.81; k₂₀ 5.84 ± 0.37; a₃₀ 23.72 ± 1.10). Interestingly, the A² allele of CSN2 was very widespread in the population; thus, it will be intriguing to verify if A²A² Agerolese cattle milk and the derived cheese may have better nutraceutical characteristics.