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

CSN1S1, CSN3 and LPL: Three Validated Gene Polymorphisms Useful for More Sustainable Dairy Production in the Mediterranean River Buffalo

The search for DNA polymorphisms useful for the genetic improvement of dairy farm animals has spanned more than 40 years, yielding relevant findings in cattle for milk traits, where the best combination of alleles for dairy processing has been found in casein genes and in DGAT1. Nowadays, similar results have not yet been reached in river buffaloes, despite the availability of advanced genomic technologies and accurate phenotype records. The aim of the present study was to investigate and validate the effect of four single nucleotide polymorphisms (SNP) in the CSN1S1, CSN3, SCD and LPL genes on seven milk traits in a larger buffalo population. These SNPs have previously been reported to be associated with, or affect, dairy traits in smaller populations often belonging to one farm. A total of 800 buffaloes were genotyped. The following traits were individually recorded, monthly, throughout each whole lactation period from 2010 to 2021: daily milk yield (dMY, kg), protein yield (dPY, kg) and fat yield (dFY, kg), fat and protein contents (dFP, % and dPP, %), somatic cell count (SCC, 103 cell/mL) and urea (mg/dL). A total of 15,742 individual milk test day records (2496 lactations) were available for 680 buffalo cows, with 3.6 ± 1.7 parities (from 1 to 13) and an average of 6.1 ± 1.2 test day records per lactation. Three out four SNPs in the CSN1S1, CSN3 and LPL genes were associated with at least one of analyzed traits. In particular, the CSN1S1 (AJ005430:c.578C>T) gave favorable associations with all yield traits (dMY, p = 0.022; dPY, p = 0.014; dFY, p = 0.029) and somatic cell score (SCS, p = 0.032). The CSN3 (HQ677596: c.536C>T) was positively associated with SCS (p = 0.005) and milk urea (p = 0.04). Favorable effects on daily milk yield (dMY, p = 0.028), fat (dFP, p = 0.027) and protein (dPP, p = 0.050) percentages were observed for the LPL. Conversely, the SCD did not show any association with milk traits. This is the first example of a confirmation study carried out in the Mediterranean river buffalo for genes of economic interest in the dairy field, and it represents a very important indication for the preselection of young bulls destined for breeding programs aimed at more sustainable dairy production.

Effect of Antioxidant Supplementation on Milk Yield and Quality in Italian Mediterranean Lactating Buffaloes

Buffaloes are raised mainly to obtain milk that is nutritionally very rich. The technological characteristics of buffalo milk are optimal for processing into cheese, and it is mainly used to produce mozzarella cheese. Under stressful conditions, buffaloes, like other animals, produce milk qualitatively poorly. The stressors that can affect the quality of production are, in addition to other factors, deficiencies in nutrients such as vitamins, antioxidants, and minerals. In this study, we evaluated the effect of antioxidant supplementation on the quality of buffalo milk. Sixty-six buffaloes were enrolled and subdivided into two balanced groups of 33 each. The ZnSe group received 0.2 kg/head/day of Bufalo Plus® containing antioxidants and barley meal, CaCO3 and MgCO3 mix; the control group was supplemented with 0.2 kg/head/day of barley meal, CaCO3 and MgCO3 mix. The two groups were fed ad libitum with a total mixed ration (TMR). The amount of diet distributed was recorded daily, and the residue in the trough manger was recorded three times per week. TMR samples were taken every two weeks for each group. Daily milk yield was recorded twice a week. Milk samples were collected every four weeks and analysed for chemical and technological properties. Furthermore, milk total antioxidant capacity was determined. The results obtained showed that the antioxidant supplement had no effect on feed intake, feeding behaviour, and feed efficiency. The treatment positively influenced milk production while it did not affect the chemical characteristics of the milk. In addition, the supplement of antioxidants improved the milk clotting properties (MCP). The supplement did not affect the antioxidant activity of the milk.

Mutation of Signal Transducer and Activator of Transcription 5 (STAT5) Binding Sites Decreases Milk Allergen αS1-Casein Content in Goat Mammary Epithelial Cells

α_S1-Casein (encoded by the CSN1S1 gene) is associated with food allergy more than other milk protein components. Milk allergy caused by α_S1-casein is derived from cow milk, goat milk and other ruminant milk. However, little is known about the transcription regulation of α_S1-casein synthesis in dairy goats. This study aimed to investigate the regulatory roles of signal transducer and activator of transcription 5 (STAT5) on α_S1-casein in goat mammary epithelial cells (GMEC). Deletion analysis showed that the core promoter region of CSN1S1 was located at −110 to −18 bp upstream of transcription start site, which contained two putative STAT5 binding sites (gamma-interferon activation site, GAS). Overexpression of STAT5a gene upregulated the mRNA level and the promoter activity of the CSN1S1 gene, and STAT5 inhibitor decreased phosphorylated STAT5 in the nucleus and CSN1S1 transcription activity. Further, GAS site-directed mutagenesis and chromatin immunoprecipitation (ChIP) assays revealed that GAS1 and GAS2 sites in the CSN1S1 promoter core region were binding sites of STAT5. Taken together, STAT5 directly regulates CSN1S1 transcription by GAS1 and GAS2 sites in GMEC, and the mutation of STAT5 binding sites could downregulate CSN1S1 expression and decrease α_S1-casein synthesis, which provide the novel strategy for reducing the allergic potential of goat milk and improving milk quality in ruminants.

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.

Complete CSN1S2 Characterization, Novel Allele Identification and Association With Milk Fatty Acid Composition in River Buffalo

The αs2-casein is one of the phosphoproteins secreted in all ruminants' milk, and it is the most hydrophilic of all caseins. However, this important gene (CSN1S2) has not been characterized in detail in buffaloes with only two alleles detected (reported as alleles A and B), and no association studies with milk traits have been carried out unlike what has been achieved for other species of ruminants. In this study, we sequenced the whole gene of two Mediterranean river buffalo homozygotes for the presence/absence of the nucleotide C (g.7539G>C) realized at the donor splice site of exon 7 and, therefore, responsible for the skipping of the same exon at mRNA level (allele B). A high genetic variability was found all over the two sequenced CSN1S2 alleles. In particular, 74 polymorphic sites were found in introns, six in the promoter, and three SNPs in the coding region (g.11072C>T, g.12803A>T, and g.14067A>G) with two of them responsible for amino acid replacements. Considering this genetic diversity, those found in the database and the SNP at the donor splice site of exon 7, it is possible to deduce at least eight different alleles (CSN1S2 A, B, B1, B2, C, D, E, and F) responsible for seven different possible translations of the buffalo αs2-casein. Haplotype data analysis suggests an evolutionary pathway of buffalo CSN1S2 gene consistent with our proposal that the published allele CSN1S2 A is the ancestral αs2-CN form, and the B2 probably arises from interallelic recombination (single crossing) between the alleles D and B (or B1). The allele CSN1S2 C is of new identification, while CSN1S2 B, B1, and B2 are deleted alleles because all are characterized by the mutation g.7539G>C. Two SNPs (g.7539G>C and g.14067A>G) were genotyped in 747 Italian buffaloes, and major alleles had a relative frequency of 0.83 and 0.51, respectively. An association study between these SNPs and milk traits including fatty acid composition was carried out. The SNP g.14067A>G showed a significant association (P < 0.05) on the content of palmitic acid in buffalo milk, thus suggesting its use in marker-assisted selection programs aiming for the improvement of buffalo milk fatty acid composition.

Genome-Wide Identification of Discriminative Genetic Variations in Beef and Dairy Cattle via an Information-Theoretic Approach

Analyzing the associations between genotypic changes and phenotypic traits on a genome-wide scale can contribute to understanding the functional roles of distinct genetic variations during breed development. We performed a whole-genome analysis of Angus and Jersey cattle breeds using conditional mutual information, which is an information-theoretic method estimating the conditional independency among multiple factor variables. The proposed conditional mutual information-based approach allows breed-discriminative genetic variations to be explicitly identified from tens of millions of SNP (single nucleotide polymorphism) positions on a genome-wide scale while minimizing the usage of prior knowledge. Using this data-driven approach, we identified biologically relevant functional genes, including breed-specific variants for cattle traits such as beef and dairy production. The identified lipid-related genes were shown to be significantly associated with lipid and triglyceride metabolism, fat cell differentiation, and muscle development. In addition, we confirmed that milk-related genes are involved in mammary gland development, lactation, and mastitis-associated processes. Our results provide the distinct properties of Angus and Jersey cattle at a genome-wide level. Moreover, this study offers important insights into discovering unrevealed genetic variants for breed-specific traits and the identification of genetic signatures of diverse cattle breeds with respect to target breed-specific properties.

An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats

The casein alpha s1 (CSN1S1) gene encodes α-s1 casein, one of the proteins constituting milk, which affects milk performance, as well as improving the absorption of calcium and bone development in mammals. A previous study found that an 11-bp insertion/deletion (indel) of this gene strongly affected litter size in goats. However, to our knowledge, the relationships between this polymorphism and the milk performance and body measurement traits of goats have not been reported. In this paper, the previously identified indel has been recognized in three Chinese goat breeds, namely the Guanzhong dairy goat (GZDG; n = 235), Shaanbei white cashmere goat (SBWC; n = 1092), and Hainan black goat (HNBG; n = 278), and the following three genotypes have been studied for all of the breeds: insertion/insertion (II), deletion/deletion (DD), and insertion/deletion (ID). The allele frequencies analyzed signified that the frequencies of the "D" allele were higher (47.8%-65.5%), similar to the previous report, which indicates that this polymorphism is genetically stable in different goat breeds. Further analysis showed that this indel was markedly associated with milk fat content, total solids content, solids-not-fat content, freezing point depression, and acidity in GZDG (p < 0.05), and also affected different body measurement traits in all three breeds (p < 0.05). The goats with II genotypes had superior milk performance, compared with the others; however, goats with DD genotypes had better body measurement sizes. Hence, it may be necessary to select goats with an II or DD genotype, based on the desired traits, while breeding. Our study provides information on the potential impact of the 11-bp indel polymorphism of the CSN1S1 gene for improving the milk performance and body measurement traits in goats.

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.

Pathway-based genome-wide association analysis of milk coagulation properties, curd firmness, cheese yield, and curd nutrient recovery in dairy cattle

It is becoming common to complement genome-wide association studies (GWAS) with gene-set enrichment analysis to deepen the understanding of the biological pathways affecting quantitative traits. Our objective was to conduct a gene ontology and pathway-based analysis to identify possible biological mechanisms involved in the regulation of bovine milk technological traits: coagulation properties, curd firmness modeling, individual cheese yield (CY), and milk nutrient recovery into the curd (REC) or whey loss traits. Results from 2 previous GWAS studies using 1,011 cows genotyped for 50k single nucleotide polymorphisms were used. Overall, the phenotypes analyzed consisted of 3 traditional milk coagulation property measures [RCT: rennet coagulation time defined as the time (min) from addition of enzyme to the beginning of coagulation; k₂₀: the interval (min) from RCT to the time at which a curd firmness of 20 mm is attained; a₃₀: a measure of the extent of curd firmness (mm) 30 min after coagulant addition], 6 curd firmness modeling traits [RCT_eq: RCT estimated through the CF equation (min); CF_P: potential asymptotic curd firmness (mm); k_CF: curd-firming rate constant (% × min^-1); k_SR: syneresis rate constant (% × min^-1); CF_max: maximum curd firmness (mm); and t_max: time to CF_max (min)], 3 individual CY-related traits expressing the weight of fresh curd (%CY_CURD), curd solids (%CY_SOLIDS), and curd moisture (%CY_WATER) as a percentage of weight of milk processed and 4 milk nutrient and energy recoveries in the curd (REC_FAT, REC_PROTEIN, REC_SOLIDS, and REC_ENERGY calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk), milk pH, and protein percentage. Each trait was analyzed separately. In total, 13,269 annotated genes were used in the analysis. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway databases were queried for enrichment analyses. Overall, 21 Gene Ontology and 17 Kyoto Encyclopedia of Genes and Genomes categories were significantly associated (false discovery rate at 5%) with 7 traits (RCT, RCT_eq, k_CF, %CY_SOLIDS, REC_FAT, REC_SOLIDS, and REC_ENERGY), with some being in common between traits. The significantly enriched categories included calcium signaling pathway, salivary secretion, metabolic pathways, carbohydrate digestion and absorption, the tight junction and the phosphatidylinositol pathways, as well as pathways related to the bovine mammary gland health status, and contained a total of 150 genes spanning all chromosomes but 9, 20, and 27. This study provided new insights into the regulation of bovine milk coagulation and cheese ability that were not captured by the GWAS.

De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery

The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.

Mediterranean river buffalo CSN1S1 gene: search for polymorphisms and association studies

The aim of the present study was to investigate the variability at CSN1S1 locus of the Italian Mediterranean river buffalo and to study possible allele effects on milk yield and its composition. Effects of parity, calving season and month of production were also evaluated. Three single-nucleotide polymorphisms were detected. The first mutation, located at position 89 of the 17th exon (c.628C>T), is responsible for the amino acid change p.Ser178 (B allele)/Leu178 (A allele). The other two polymorphisms, detected at the positions 144 (c.882G>A) and 239 (c.977A>G) of 19th exon, respectively, are silent (3ʹ UTR, untranslated region). Associations between the CSN1S1 genotypes and milk production traits were investigated using 4122 test day records of 503 lactations from 175 buffalo cows. Milk yield, fat and protein percentages were analysed using a mixed linear model. A significant association between the c.628C>T SNP and the protein percentage was found. In particular, the CC genotype showed an average value ~0.04% higher than the CT and TT genotypes. The allele substitution effect of cytosine into thymine was –0.014, with a quite low (0.3%) protein percentage contribution to total phenotypic variance. A large dominance effect was detected. Characterisation of the CSN1S1 transcripts and a method based on MboI amplification created restriction site PCR for a rapid genotyping of c.628C>T are provided.

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.

Short communication: CSN1S1-CSN3 (α(S1)-κ-casein) composite genotypes affect detailed milk protein composition of Mediterranean water buffalo

The aim of the study was to investigate the effect of composite CSN1S1-CSN3 [α(S1)-κ-casein (CN)] genotype on milk protein composition in Mediterranean water buffalo. Content of α(S1)-CN, α(S2)-CN, β-CN, γ-CN, κ-CN, glycosylated and unglycosylated κ-CN, α-lactalbumin, and β-lactoglobulin was measured by reversed-phase HPLC using 621 individual milk samples. Genotypes at CSN1S1 and CSN3 were also obtained by reversed-phase HPLC. Two alleles were detected at CSN1S1 (corresponding to the A and B variants, O62823: p.Leu193Ser,) and at CSN3 (corresponding to the X1 and X2 variants, CAP12622.1: p.Ile156Thr). Increased proportions of α(S1)-CN in total casein (TCN) were associated with genotypes carrying CSN1S1 A. Genotypes associated with a marked decrease of the proportion of α(S1)-CN in TCN (composite genotypes AB-X1X1 and BB-X1X2) were associated with marked increases in the proportion of α(S2)-CN. In addition, composite genotypes carrying the X1 allele at CSN3 were associated with a greater proportion of α(S2)-CN in TCN relative to those carrying CSN3 X2. Composite genotypes greatly affected also the variability of ratios of κ-CN to TCN, with genotypes carrying the X1 allele at CSN3 being associated with decreased ratios. The decreased content of glycosylated κ-CN associated with CSN3 X1 was responsible for the overall lower content of total κ-CN in milk of X1-carrying animals. Increasing the frequency of specific genotypes might be an effective way to alter milk protein composition, namely the proportion of α(S1)-CN, α(S2)-CN, and κ-CN in TCN, and the degree of glycosylation of κ-CN.