Kappa-casein polymorphisms in Indian dairy cattle and buffalo: a new genetic variant in buffalo

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.

Polymorphisms of the kappa casein (CSN3) gene and inference of its variants in water buffalo (Bubalus bubalis)

Kappa casein plays a crucial role in the formation of stable casein micelles and has a key influence on milk-clotting properties. However, current understanding of buffalo CSN3 gene polymorphisms is not sufficient. In this study, the polymorphisms in the complete coding sequence (CDS) of the buffalo CSN3 were detected using PCR product direct sequencing. The CDS of CSN3 for river and swamp buffalo was the same in length, which contained an open reading frame of 573 nucleotides encoding a peptide containing 190 amino acid residues. A total of eight single nucleotide polymorphisms (SNPs) was identified in two types of buffalo. Among them, c.86C>T, c.252G>C, c.445G>A, c.467C>T and c.516A>C were non-synonymous, which leads to p.Pro8Leu, p.Lys63Asn, p.Val128Ile, p.Thr135Ile and p.Glu151Asp substitutions in buffalo kappa casein ( κ -CN), respectively. The substitution of p.Thr135Ile may exert a vital effect on the function of buffalo κ -CN. Eleven haplotypes were defined based on the SNPs found in buffalo, and accordingly, seven protein variants and four synonymous variants of buffalo κ -CN were inferred, called variants A, B, B 1 , C, C 1 , C 2 , D, E, F, F 1 and G. The variants observed in water buffalo did not exist in the Bos genus. In addition, 14 amino acid differential sites of κ -CN between buffalo and the Bos genus were identified, of which 3 were located at glycosylation sites (80S, 96T, 141S) and 4 at phosphorylation sites (19S, 80S, 96T, 141S). It is speculated that they may lead to differences in the physicochemical properties of κ -CN between buffalo and the Bos genus. This study will lay a foundation for exploring the association between the variation in the CSN3 gene and the lactation traits of buffalo.

Analysis of genetic variations across regulatory and coding regions of kappa-casein gene of Indian native cattle (Bos indicus) and buffalo (Bubalus bubalis)

The promoter region of kappa-casein (κ-CN) gene in Indian native cattle and buffalo breeds was sequenced and analyzed for nucleotide variations. Sequence comparison across breeds of Indian cattle revealed a total of 7 variations in the promoter region, of which − 515 G/T, − 427 C/T, − 385 C/T, − 283 A/G and − 251 C/T were located within consensus binding sites for octamer-binding protein (OCT1)/pregnancy specific mammary nuclear factor (PMF), activator protein-2 (AP2), hepatocyte nuclear factor (HNF-1) and GAL4 transcription factors (TFs), respectively. These variations might be involved in gain or loss of potential transcription factor binding sites (TFBSs). Unlike the other 4 variants, the − 283 (A/G) variant located within HNF-1 TFBS was specific to Indian cattle as this change has not been observed in the Bos taurus sequence. Other TFBSs viz., MGF, TBP, NF-1, milk box and C/EBP were conserved across species. For the Indian native buffalo breeds, only 3 changes were identified in the promoter region; − 305 (A/C), − 160 (T/C) and − 141 (A/G) and most of the TFBSs were found to be conserved. However, deletion of two adjacent nucleotides located in and around binding site for C/EBP TF was identified in buffalo when compared with promoter sequence of bovine κ-CN. For κ-CN of Indian native cattle, a strong linkage disequilibrium (LD) was observed for variations 515 G/T, − 427 C/T and − 385 C/T in the promoter region; and for variations at codons 136 and 148 of exon-IV. Further, among intragenic haplotypes, variation − 427 C/T was found to be in LD with variations at codons 136 and 148. The information generated in the present work provides comprehensive characterization of κ-CN gene promoter and coding regions in Indian cattle and buffaloes and reported variations could become important candidates for carrying out further research in dairy traits.

•

The κ-casein gene (κ-CN) and its promoter region was sequence characterized in 15 Indian native cattle and 8 buffalo breeds.

•

Among the identified variations, four were located within TFBSs in Indian cattle, while all TFBSs were conserved in buffaloes.

•

The κ-CN CDS was highly conserved with only 3 and 2 non-synonymous changes for Indian cattle and buffaloes, respectively.

•

Across promoter and CDS, − 427 C/T was found to be in LD with variation at codons 136 and 148.

Effect of CSN1S1-CSN3 (α(S1)-κ-casein) composite genotype on milk production traits and milk coagulation properties in Mediterranean water buffalo

The aim of this study was to estimate effects of CSN1S1-CSN3 (α(S1)-κ-casein) composite genotypes on milk production traits and milk coagulation properties (MCP) in Mediterranean water buffalo. Genotypes at CSN1S1 and CSN3 and coagulation properties [rennet clotting time (RCT), curd firming time (K₂₀), and curd firmness (A₃₀)] were assessed by reversed-phase HPLC and computerized renneting meter analysis, respectively, using single test-day milk samples of 536 animals. Alternative protein variants of α(S1)-CN and κ-CN were detected by HPLC, and identification of the corresponding genetic variants was carried out by DNA analysis. Two genetic variants were detected at CSN1S1 (A and B variants) and 2 at CSN3 (X1 and X2 variants). Statistical inference was based on a linear model including the CSN1S1-CSN3 composite genotype effect (7 genotypes), the effects of herd-test-day (8 levels), and a combined days in milk (DIM)-parity class. Composite genotype AB-X2X2 was associated with decreased test-day milk yield [-0.21 standard deviation (SD) units of the trait] relative to genotype BB-X2X2. Genotypes did not affect milk protein content, but genotype AB-X1X1 was associated with increased fat content compared with genotype BB-X2X2 (+0.28 SD units of the trait) and AB-X1X1 (+0.43 SD units of the trait). For RCT, the largest difference (+1.91 min; i.e., 0.61 SD units of the trait) was observed between genotype AA-X1X2 and AB-X1X1. Direction of genotype effects on K(20) was consistent with that for RCT. The maximum variation in K₂₀ due to genotype effects (between AA-X1X2 and AB-X1X1 genotypes) was almost 0.9 SD units of the trait. Magnitude of genotype effects was smaller for A₃₀ than for RCT and K₂₀, with a maximum difference of 0.5 SD units of the trait between genotype AA-X1X2 and AA-X1X1. The B allele at CSN1S1 was associated with increased RCT and K₂₀ and with weaker curds compared with allele A. Allele X2 at CSN3 exerted opposite effects on MCP relative to CSN1S1 B. Because of linkage disequilibrium, allele B at CSN1S1 and allele X2 at CSN3 tend to be associated and this likely makes their effects cancel each other. This study indicates a role for casein genes in variation of MCP of buffalo milk. Further studies are necessary to estimate the effects of casein genetic variants on variation of cheese yield.

Separation and quantification of water buffalo milk protein fractions and genetic variants by RP-HPLC

A RP-HPLC method, developed for the separation and quantification of the most common genetic variants of bovine milk proteins, was successfully applied to the analysis of water buffalo milk. All the most common buffalo casein and whey proteins fractions, as well as their genetic variants, were detected and separated simultaneously in 40 min. Purified buffalo proteins were used as calibration standards and a total of 536 individual milk samples were analysed for protein composition. α(S1)-, α(S2)-, βγ-, and κ-casein were 32.2%, 15.8%, 36.5%, and 15.5%, respectively, of total casein content, whereas content of β-Lactoglobulin was approximately 1.3 times as high as that of α-Lactalbumin. The existence of a polymorphism of κ-casein was demonstrated in Mediterranean water buffalo and α(S1)- and κ-casein genetic variants were successfully detected by RP-HPLC.

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.

Short communication: new alleles of the bovine kappa-casein gene revealed by resequencing and haplotype inference analysis

We tested the hypothesis that extensive undiscovered genetic diversity exists in important functional genes from domestic and wild cattle species (Bos spp.). We resequenced 483 bp of a key exon (exon IV) from the kappa (kappa)-casein gene (CSN3) for a panel of samples of domestic cattle from 8 countries and a close relative species, the gayal (Bos frontalis). Six single nucleotide polymorphisms were identified. Haplotype inference revealed 12 haplotypes, of which 8 were newly discovered. Among these 8 new haplotypes, 5 differed by one nonsynonymous mutation and 3 differed by one silent mutation from previously well-characterized CSN3 alleles. From those, one was shared by the gayal and Zebu, was different from CSN3*B at position Ile136Thr, and showed a close phylogenetic relationship with the banteng, gaur, and yak. The other 7 new haplotypes were detected in our panel of worldwide local cattle breeds but were absent from previously reported commercial breeds. These results support the hypothesis that genetic diversity at the coding region of CSN3 has been underestimated. This study also highlights how important it is to resequence functionally important genes in worldwide local cattle breeds, many of which are threatened by extinction or replacement by commercial breeds.

A proteomic characterization of water buffalo milk fractions describing PTM of major species and the identification of minor components involved in nutrient delivery and defense against pathogens

Water buffalo has been studied in relation to the exclusive use of its milk for the manufacture of high-quality dairy products. Buffalo milk presents physicochemical features different from that of other ruminant species, such as a higher content of fatty acids and proteins. We report here a detailed proteomic analysis of buffalo skim milk, whey and milk fat globule membrane fractions. Notwithstanding the poor information available on buffalo genome, identification of protein isoforms corresponding to 72 genes was achieved by a combined approach based on 2-DE/MALDI-TOF PMF and 1-DE/muLC-ESI-IT-MS-MS. Major protein components, i.e. alpha(Sl)-, alpha(S2)-, beta-, kappa-caseins, alpha-lactalbumin and beta-lactoglobulin, were characterized for PTM, providing a scientific basis to coagulation/cheese making processes used in dairy productions. Minor proteins detected emphasized the multiple functions of milk, which besides affording nutrition to the newborn through its major components, also promotes development and digestive tract protection in the neonate, and ensures optimal mammary gland function in the mother. Defense against pathogens is guaranteed by an arsenal of antimicrobial/immunomodulatory proteins, which are directly released in milk or occur on the surface of secreted milk-lipid droplets. Proteins associated with cell signaling or membrane/protein trafficking functions were also identified, providing putative insights into major secretory pathways in mammary epithelial cells.

Sequence analysis of UTR and coding region of kappa-casein gene of Indian riverine buffalo (Bubalus bubalis)

In this study, complete nucleotide as well as derived amino acid sequence characterization of water buffalo (Bubalus bubalis) kappa-casein gene has been presented. Kappa-casein cDNA clones were identified and isolated from a buffalo lactating mammary gland cDNA library. Sequence analysis of kappa-casein cDNA revealed 850 nucleotides with an open reading frame (ORF) of 573 nucleotides, encoding mature peptide of 169 amino acids. The 5' untranslated region (UTR) comprised 71 nucleotides, while 3' UTR was of 206 nucleotides. A total of 11 nucleotide and seven amino acid changes were observed in, buffalo (Bubalus bubalis) as compared to cattle (Bos taurus), sheep (Ovis aries) and goat (Capra hircus). Among these nucleotide changes, eight were unique in buffalo as they were fully conserved in cattle, sheep and goat. Majority of the nucleotide changes and all the amino acid changes; 14 (Asp-Glu), 19(Asp/Ser-Asn), 96(Ala-Thr), 126(Ala-Val), 128(Ala/Gly-Val), 156(Ala/Pro-Val) and 168(Ala/Glu-Val) were limited to exon IV. Three glycosylation sites, Thr 131, Thr 133 and Thr 142 reported in cattle and goat kappa-casein gene were also conserved in buffalo, however, in sheep Thr 142 was replaced by Ala. Chymosin hydrolysis site, between amino acids Phe 105 and Met 106, important for rennet coagulation process, were found to be conserved across four bovid species. Buffalo kappa-casein with the presence of amino acids Thr 136 and Ala 148 seems to be an intermediate of "A" and "B" variants of cattle. Comparison with other livestock species revealed buffalo kappa-casein sharing maximum nucleotide (95.5%) and amino acid (92.6%) similarity with cattle, whereas with pig it showed least sequence similarity of 76.0% and 53.2%, respectively. Phylogenetic analysis based on both nucleotide and amino acid sequence indicated buffalo kappa-casein grouping with cattle, while sheep and goat forming a separate cluster close to them. The non-ruminant species viz. camel, horse and pig were distantly placed, in separate lineages.

SSCP analysis at the bovine CSN3 locus discriminates six alleles corresponding to known protein variants (A, B, C, E, F, G) and three new DNA polymorphisms (H, I, A1)

A high resolution SSCP protocol was developed for simultaneous discrimination of the known CSN3 alleles A, B, C, E, F and G. Furthermore, three new DNA polymorphisms were identified in different Bos taurus and Bos indicus breeds or crosses. Mendelian segregation was shown for two of these polymorphisms (named CSN3*H and 1), and the third (named CSN3*A1) was found in unrelated animals, thus indicating the presence of three additional alleles at the bovine CSN3 locus. DNA sequencing revealed single mutations that led to a Thr/Ile substitution in amino acid position 135 for CSN3*H and to a Ser/Ala substitution in position 104 of the deduced amino acid sequence of CSN3*1 (GenBank accession numbers AF105260 and AF121023) compared to CSN3*A. In CSN3*A1, a silent mutation in the third codon position of Pro150 was found (GenBank accession number AF092513).