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.

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.

Feeding System Resizes the Effects of DGAT1 Polymorphism on Milk Traits and Fatty Acids Composition in Modicana Cows

Simple Summary

Genetic selection for single-locus polymorphisms could offer suitable opportunities to rapidly improve milk traits in local unselected cattle breeds characterized by low production levels. Since these hardy breeds are generally raised in traditional extensive and semi-intensive systems, which make wide use of grazing resources, the interactive effect between genotype and feeding system is worthy of investigation. In Modicana cattle breed, milk composition and fatty acid profile were influenced by both genetic polymorphisms at the DGAT1 K232A locus and feeding systems. The milk from homozygous AA cows was associated with a more favorable fatty acid composition due to a lower percentage of total saturated fatty acids, saturated to unsaturated ratio, atherogenic index, and a greater presence of oleic acid and total unsaturated fatty acids. Our finding confirmed the important role of pasture feeding on milk composition: the high nutritional and healthy value of milk obtained in extensive systems by pasture-fed cows. The interaction between the two experimental factors also appears to play a role: in our experimental condition, it seems that high pasture feeding can resize the effect of the DGAT1 genotype on milk traits and fatty acid composition in Modicana cows.

Abstract

The interaction between genetic polymorphism and feeding system on milk traits and fatty acid composition was investigated in Modicana cows. Two DGAT1 K232A genotypes (AK and AA) and two feeding regimes, extensive system (EX) with 8 h of grazing without concentrate (EX) and semi-intensive systems (SI) with 2 h of grazing with concentrate, were investigated. DGAT1 genotype did not influence milk yield and composition. The feeding system affected milk composition: protein was significantly higher in SI and lactose in the EX system. A significant genotype × feeding system interaction was observed: the protein and casein levels of AK cows were higher in the SI compared to the EX system. Milk fatty acids profile, total saturated to total unsaturated fatty acids, n-6 to n-3 ratios, and atherogenic index were affected by the feeding system, improving the healthy properties of milk from animals reared in the extensive system. DGAT1 genotype influenced the fatty acid composition: milk from AA cows had a more favorable fatty acid composition due to lower total saturated fatty acids, saturated to unsaturated ratio, atherogenic index, and higher levels of oleic acid and total unsaturated fatty acids. Furthermore, an interaction genotype x feeding system was observed: the AK milk was richer in short-chain FAs (C4:0–C8:0) and C10:0 only in the EX but not in the SI system. Our data suggest that a high amount of green forage in the diet of Modicana cows can resize the effect of the DGAT1 genotype on milk traits and fatty acids composition.

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.

Effect of GH p.L127V Polymorphism and Feeding Systems on Milk Production Traits and Fatty Acid Composition in Modicana Cows

Growth hormone participates in the regulation of lactation and lipid metabolism. A trial study was conducted to evaluate the effects of genetic polymorphism at GH p.L127V and its interaction with feeding system (extensive, EX; semi-intensive, SI) on milk traits and fatty acids composition in Modicana cows. In the semi-intensive farm (SI) diet consisted of hay, concentrate and 2 h of grazing. In the extensive farm (EX) feeding consisted in 8 h of grazing and hay. The frequencies of LL, LV and VV genotypes were, respectively: 0.64, 0.34, 0.02. GH polymorphism and its interaction with feeding system did not influence milk yield and composition. Cows carrying LL genotype produced milk with lower 6:0 and 8:0 and higher 16:1 c9 and 18:1 c9, total UFA and total MUFA. Feeding significantly affected fatty acids: in EX cows lower SFA and higher PUFA and UFA were found, compared to SI cows. The lower, more favorable atherogenic index of milk from EX system was coherent with the improved healthy characteristics of milk from animals fed almost exclusively on pasture. A significant interaction genotype x feeding system was evident for 18:1, higher in the LL cows only in the EX system, but not in the SI system.

Genetic polymorphisms at candidate genes affecting fat content and fatty acid composition in Modicana cows: effects on milk production traits in different feeding systems

Feeding greatly affects milk yield and composition. The research is highlighting the potential of genetic polymorphism at some loci to affect milk yield and quality traits. These loci can be up/down regulated depending on the production environment; therefore, we hypothesized that milk yield and composition could differ when cows with different genotype at SCD, DGAT1 and ABCG2 loci are reared in different feeding systems. The polymorphisms of SCD, DGAT1 and ABCG2 genes were investigated in Modicana breed. In all, three polymorphic sites, responsible for the genetic variation of quantitative trait loci and therefore defined quantitative trait nucleotides, were genotyped: the transition g.10329C&gt;T in 5th exon determines a substitution p.A293V in the SCD, the dinucleotide mutation g.10433-10434AA&gt;GC in 8th exon responsible for p.K232A substitution in the DGAT1 and the transition g.62569A&gt;C in the 14th exon responsible for p.Y581S substitution in the ABCG2 gene. In the sample of 165 Modicana cows, SCD and DGAT1 genes resulted polymorphic; the alleles g.10329T and g.10433-10434GC were the most frequent in SCD and DGAT1 (0.73 and 0.91) respectively, whereas ABCG2 locus was monomorphic for allele A (p.581Y). Sequencing analysis was carried out on 14 samples with different genotypes to confirm the results of the PCR-RFLP protocols. Based on the genotypes at SCD locus, 47 Modicana cows were selected for the nutritional trial: 24 cows in a semi-intensive farm, with 2 h/day grazing on natural pasture, and 23 cows in an extensive farm, with 8 h/day grazing on natural pasture. Monthly, milk yield and composition were evaluated and individual milk samples were analyzed for fatty acids composition by gas chromatography. No differences in milk yield, fat, protein, lactose, casein and urea were associated to SCD genotype. Feeding systems affected milk yield and composition. No significant genotype×feeding system interaction was observed for milk yield and composition. Fatty acids composition was significantly affected only by the feeding system. Significant interactions were found between SCD genotype and feeding system for six fatty acids: 4:0, 6:0, 8:0, 10:0, 12:0 and t11 18:1. We concluded that the feeding system was the factor that mostly affected milk production and composition; moreover, our results do not confirm what reported in literature as regard the effect of the SCD polymorphism on milk fatty acid composition. The high amount of pasture seemed to have resized the SCD polymorphism effects because of the different fatty acids composition of the diet.

Milk protein polymorphisms and casein haplotypes in Butana cattle

Butana is a Bos indicus dairy cattle breed that is well adapted to the local environment of Sudan. The breed has been gradually declining in number due to breed substitution. Therefore, conservation and improvement strategies are required to maintain this breed. The aim of the present study was to assess genetic variation that is characteristic for Butana cattle in the milk protein genes CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB. In a first step, genomic DNA of five unrelated individuals was comparatively sequenced across all exon and flanking sequences. Ninety-three single nucleotide polymorphisms (SNPs) were identified in Butana cattle compared with the Bos taurus reference sequence at Ensembl. We confirmed the recently identified protein variants CSN2*J, CSN2*L, and LALBA*E. Fifty-two SNPs in non-coding regions are novel. Among the novel SNPs, five are located in promoter regions, three of them are in putative transcription factor binding sites (TFBSs) of the CSN1S2 promoter. Fifteen SNPs potentially affect miRNA target sites. In a second step, 50 unrelated Butana cattle were genotyped. This allowed deriving haplotypes for the casein gene cluster on BTA6. The most frequent haplotype was CSN1S1*C-CSN2*A ² -CSN1S2*A-CSN3*A (C-A ² -A-A, frequency 0.1546). Considering the newly identified CSN1S2 promoter variants, the most frequent haplotype was C-A ² -TTC-A-A (0.1046), with TTC as the promoter variant. The information on protein and promoter variants can be used for the development of conservation and breeding strategies for this local breed.

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.

Molecular characterization of bovine CSN1S2*B and extensive distribution of zebu-specific milk protein alleles in European cattle

The B allele of the bovine alpha (S2)-casein gene (CSN1S2) was characterized at the molecular level and the distribution of zebu-specific milk protein alleles was determined in 26 cattle breeds originating from 3 continents. The CSN1S2*B allele is characterized by a C --> T transition affecting nucleotide 17 of exon 3, which leads to a change in the eighth amino acid of the mature protein, from Ser to Phe (i.e., TCC --> TTC). DNA-based methods were developed to identify carriers of CSN1S2*B and the other alleles (CSN1S2*A, C, and D) at the same locus. CSN1S2*B and other zebu-specific milk protein alleles and casein haplotypes are widely distributed in European cattle breeds, particularly those of southeastern origin. Alleles CSN1S2*B and CSN3*H are important in searching for zebu imprints in European cattle breeds. Diversity estimates at the milk protein loci were highest in the zebus followed by southeastern European taurines. Anatolian Black had the highest number of zebu alleles among European taurines. Common, group, and intergroup haplotypes occurred in the breeds and demonstrated relationships that concurred with developmental histories, genetic makeup, and, in particular, exposed the extent of zebu influence on southeastern European cattle.