Invited review: milk protein polymorphisms in cattle: effect on animal breeding and human nutrition

The Aggregated and Micellar Forms of β-Casein Purified from Donkey and Bovine Milk Present Potential as Carriers for Bioactive Nutritional Compounds

In recent years, there has been a growing interest in the pure casein fraction of milk protein, particularly β-casein due to its physicochemical properties as well as its bio- and techno-functional properties. The utilization of self-assembled β-caseins from bovine origin as nanocarriers for the delivery of nutraceutical compounds or drugs has increased dramatically. Concerning β-caseins from other milk sources, the use of hypoallergenic donkey β-caseins as a potential delivery vehicle for nutraceutical hydrophobic compounds is beginning to generate interest. The present review deals with casein micelles models, bovine and donkey β-casein molecular structures, as well as their physical-chemical properties that account for their exploitation in nutraceutics and pharmaceutics. This review work suggests the possibility of developing delivery systems for hydrophobic bioactive compounds using β-casein purified from hypoallergenic donkey milk, highlighting the potential of this protein as an innovative and promising vehicle for enhancing the enrichment and bioavailability of various bioactive substances in food products.

Transcriptomic profiling of gastrointestinal tracts in dairy cattle during lactation reveals molecular adaptations for milk synthesis

During lactation, dairy cattle's digestive tract requires significant adaptations to meet the increased nutrient demands for milk production. As we attempt to improve milk-related traits through selective pressure, it is crucial to understand the biological functions of the epithelia of the rumen, small intestine, and colonic tissues in response to changes in physiological state driven by changes in nutrient demands for milk synthesis. In this study, we obtained a total of 108 transcriptome profiles from three tissues (epithelia of the colon, duodenum, and rumen) of five Holstein cows, spanning eight time points from the early, mid, late lactation periods to the dry period. On average 97.06% of reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We analyzed 27,607 gene expression patterns at multiple periods, enabling direct comparisons within and among tissues during different lactation stages, including early and peak lactation. We identified 1645, 813, and 2187 stage-specific genes in the colon, duodenum, and rumen, respectively, which were enriched for common or specific biological functions among different tissues. Time series analysis categorized the expressed genes within each tissue into four clusters. Furthermore, when the three tissues were analyzed collectively, 36 clusters of similarly expressed genes were identified. By integrating other comprehensive approaches such as gene co-expression analyses, functional enrichment, and cell type deconvolution, we gained profound insights into cattle lactation, revealing tissue-specific characteristics of the gastrointestinal tract and shedding light on the intricate molecular adaptations involved in nutrient absorption, immune regulation, and cellular processes for milk synthesis during lactation.

Effects of A1 Milk, A2 Milk and the Opioid-like Peptide β-Casomorphin-7 on the Proliferation of Human Peripheral Blood Mononuclear Cells

Special attention is given to cow's milk and its variants, with ongoing discussions about health-related impacts primarily focusing on the A1 variant in contrast to the A2 variant. The difference between these variants lies in a single amino acid alteration at position 67 of β-casein. This alteration is presumed to make the A1 variant more susceptible to enzymatic breakdown during milk digestion, leading to an increased release of the peptide β-casomorphin-7 (BCM-7). BCM-7 is hypothesized to interact with µ-opioid receptors on immune cells in humans. Although BCM-7 has demonstrated both immunosuppressive and inflammatory effects, its direct impact on the immune system remains unclear. Thus, we examined the influence of A1 and A2 milk on Concanavalin A (ConA)-stimulated human peripheral blood mononuclear cells (PBMCs), as well as the effect of experimentally digested A1 and A2 milk, containing different amounts of free BCM-7 from β-casein cleavage. Additionally, we evaluated the effects of pure BCM-7 on the proliferation of ConA-stimulated PBMCs and purified CD4+ T cells. Milk fundamentally inhibited PBMC proliferation, independent of the β-casein variant. In contrast, experimentally digested milk of both variants and pure BCM-7 showed no influence on the proliferation of PBMCs or isolated CD4+ T cells. Our results indicate that milk exerts an anti-inflammatory effect on PBMCs, regardless of the A1 or A2 β-casein variant, which is nullified after in vitro digestion. Consequently, we deem BCM-7 unsuitable as a biomarker for food-induced inflammation.

Developmental programming of production and reproduction in dairy cows: IV. Association of maternal milk fat and protein percentage and milk fat to protein ratio with offspring's birth weight, survival, productive and reproductive performance and AMH concentration from birth to the first lactation period

Although the association of maternal milk production with developmental programming of offspring has been investigated, there is limited information available on the relationship of maternal milk components with productive and reproductive performance of the offspring. Therefore, the present study was conducted to analyze the association of maternal milk fat and protein percentage and milk fat to protein ratio with birth weight, survival, productive and reproductive performance and AMH concentration in the offspring. In study I, data of birth weight, milk yield and reproductive variables of offspring born to lactating dams (n = 14,582) and data associated with average maternal milk fat percentage (MFP), protein percentage (MPP) and fat to protein ratio (MFPR) during 305-day lactation were retrieved. Afterwards, offspring were classified in various categories of MFP, MPP and MFPR. In study II, blood samples (n = 339) were collected from offspring in various categories of MFP, MPP and MFPR for measurement of serum AMH. Maternal milk fat percentage was positively associated with birth weight and average percentage of milk fat (APMF) and protein (APMP) and milk fat to protein ratio (FPR) during the first lactation, but negatively associated with culling rate during nulliparity in the offspring (P < 0.05). Maternal milk protein percentage was positively associated with birth weight, APMF, APMP, FPR and culling rate, but negatively associated with milk yield and fertility in the offspring (P < 0.05). Maternal FPR was positively associated with APMF and FPR, but negatively associated with culling rate, APMP and fertility in the offspring (P < 0.05). However, concentration of AMH in the offspring was not associated with MFP, MPP and MFPR (P > 0.05). In conclusion, the present study revealed that maternal milk fat and protein percentage and their ratio were associated with birth weight, survival, production and reproduction of the offspring. Yet it was a preliminary research and further studies are required to elucidate the mechanisms underlying these associations.

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.

Effects of feeding whole-cracked rapeseeds, nitrate, and 3-nitrooxypropanol on protein composition, minerals, and vitamin B in milk from Danish Holstein cows

The present study was conducted to assess the individual or combined effects of feeding dietary fat (whole-cracked rapeseed), nitrate, and 3-nitrooxypropanol (3-NOP) on protein profile, mineral composition, B vitamins, and nitrate residues in milk from dairy cows. Forty-eight Danish Holstein cows used in an 8 × 8 incomplete Latin square design were fed 8 factorially arranged diets ((30 or 63 g crude fat/kg DM) × (0 or 10 g nitrate/kg DM) × (0 or 80 mg 3-NOP/kg DM)) over 6 periods of 21 d each. In each period, milk samples were collected from individual cows during the third week by pooling milk obtained from 4 consecutive milkings, and analyzed for protein profile including protein modifications, mineral composition, riboflavin, cobalamin, and presence of nitrate residues. Fat supplementation led to an increase in the phosphorylation degree of αS1-CN by 8.5% due to a decreased relative proportion of αS1-CN 8P and an increased relative proportion of αS1-CN 9P and further to a decrease in the relative proportion of αS2-CN by 2.4%. Additionally, fat supplementation decreased the relative proportions of glycosylated and unglycosylated forms of κ-CN, consequently leading to a 3.6% decrease in total κ-CN. In skim milk, K, Ca, P, and Mg concentrations were altered by individual use of fat, nitrate, and 3-NOP. Feeding nitrate resulted in a 5.4% increase in riboflavin concentration in milk while supplementing 3-NOP increased cobalamin concentration in milk by 21.1%. The nitrate concentration in milk was increased upon feeding nitrate however, this increased concentration was well below the maximum permissible limit of nitrate in milk (<50 mg/L). In conclusion, no major changes were observed in milk protein, and mineral compositions by feeding fat, nitrate, and 3-NOP to dairy cows while the increased riboflavin and cobalamin by nitrate and 3-NOP, respectively, could be of beneficial nutritional value for milk consumers.

The A1/A2 β-casein genotype of cows, but not their horn status, influences peptide generation during simulated digestion of milk

The effect of the horn status of cows on their milk composition and quality is a controversial research topic. In this study, 128 milk samples from 64 horned and 64 disbudded Brown Swiss and Original Braunvieh cows were collected from alpine farms where both horned and disbudded cows were grazing on mountain pastures. The samples were analyzed for their detailed composition and protein digestion in a simulated in vitro digestion (INFOGEST). To exclude probable influences on digestion, the β-casein genotype with its variants A1 and A2 was also included in the study. The effects of horn status and β-casein genotype were investigated in linear mixed models, which included additional influencing random factors such as Original Braunvieh blood proportion, stage of lactation, and farm. Horn status did not have any effect on milk composition or digestion. In contrast, milk from A1A1 cows showed a different protein digestion than milk of A1A2 and A2A2 cows in the gastric phase, including smaller amounts of β-casomorphin(BCM)21-associated peptides and larger amounts of BCM11-associated peptides. Abundances of BCM7 did not differ between β-casein genotypes. At the end of the intestinal phase, the digested milk of A1A1 and A2A2 b-casein genotypes did not differ.

Associations Between BLG, CSN3, DGAT1, GH, PIT1, and PRL Gene Polymorphisms and Milk Production Traits in Holstein Dairy Cows: A Meta-Analysis

The current meta-analysis was employed to combine the results of multiple studies into a single estimate related to B-LG, CSN3, DGAT1, PRL, GH, and PIT1 gene polymorphisms and their effects on milk production traits. The purpose of this meta-analysis was to investigate associations between B-LG, CSN3, DGAT1, GH, PIT1, and PRLgene polymorphisms with milk production traits in Holstein dairy cows. An extensive search was done from 2002 to 2022 year. Statistical analyses were performed by using Stata 11.2 software. Genetic models viz codominant (AA vs. AB, AA vs. BB, and AB vs. BB), dominant (AA + AB vs. BB), recessive (AA vs. AB + BB), and completely over-dominant (AA + BB vs. AB) were applied. The results of meta-analysis of association between B-LG genotypes with milk yield where found a significant (P < 0.05) and with fat and protein contents (P < 0.01). In CSN3 polymorphisms of A/A and A/B genotypes had a significant effect on fat yield (P < 0.05) and protein content (P < 0.01). DGAT1 polymorphisms had a significant effect on milk yield, fat yield, protein yield (P < 0.05), with fat and protein contents showed high effect (P < 0.01). No significant association was detected between GH and milk traits (P > 0.05). PIT1 genotype polymorphisms had a significant effect on milk yield (P < 0.05) and protein content (P < 0.01). PRL genotype polymorphisms were significantly associated with milk yield (P < 0.05), fat content and protein yield (P < 0.01). The B-LG, DGAT1,CSN3 and PRL gene polymorphisms could be utilized as good markers to improve milk production traits in the Holstein cattle breed.

Effects of protein genetic variants on their phosphorylation levels, milk composition, milk proteome, and milk coagulation ability in Chinese Holstein bovine milk

Milk samples were collected from 3625 Chinese Holstein cows to assess the effects of κ-casein (κ-CN) and β-lactoglobulin (β-LG) genetic variants on its milk coagulation properties. The results show that Chinese Holstein cows have a higher frequency of the κ-CN AA and AB variants, and β-LG of the AB and AA variants. Of these, κ-CN B variants, the β-LG AA and BB variants were more frequent in milk showing good coagulation. The effects of the genetic variants on milk composition, milk proteome, and protein phosphorylation sites were studied. The results showed that higher concentrations of protein and dry matter were found in κ-CN BE variant. Moreover, large variations in milk proteome among different κ-CN and β-LG variants were observed. Highly phosphorylated for κ-CN, especially Ser97, was observed in cows with the κ-CN BE variant, but no effect of β-LG variants on phosphorylation site was found. Of the various factors examined, variation of κ-CN phosphorylation sites Ser97 may be the most important in affecting casein structure and milk coagulation ability. Some milk protein contents were found to be negative factors for milk coagulation. In summary, this study showed that κ-CN genetic variants contained different milk compositions and phosphorylation site Ser97 influenced milk coagulation.

Effect of composition, casein genetic variants and glycosylation degree on bovine milk whipping properties

This study investigated the individual and combined effects of ĸ-Casein (ĸ-CN; AA, AB, BB), β-Casein (β-CN; A1A1, A1A2, A2A2) and high and low ratios of glycosylated ĸ-CN to total ĸ-CN, referred to as the glycosylation degree (GD), on bovine cream whipping properties. The genetic variants of individual cows were identified using reversed-phase high-performance liquid chromatography (RP-HPLC) and verified through liquid chromatography-mass spectrometry (LC-MS). A previously discovered relationship between days-in-milk and GD was validated and used to obtain high and low GD milk. Whipped creams were created through the mechanical agitation of fat standardised cream from milk of different ĸ-CN, β-CN, and GD combinations, and whipping properties (the ability to whip, overrun, whipping time and firmness) were evaluated. No significant correlation was measured in whipping properties for cream samples from milks with different ĸ-CN and β-CN genetic variants. However, 80 % of samples exhibiting good whipping properties (i.e., the production of a stiffened peak) were from milk with low GD suggesting a correlation between whipping properties and levels of glycosylation. Moreover, cream separated from skim milk of larger casein micelle size showed superior whipping properties with shorter whipping times (<5 min), and higher firmness and overrun. Milk fat globule (MFG) size, on the other hand, did not affect whipping properties. Results indicate that the GD of κ-CN and casein micelle size may play a role in MFG adsorption at the protein and air interface of air bubbles formed during whipping; hence, they govern the dynamics of fat network formation and influencing whipping properties.

Genetic variability among and within domestic Old and New World camels at the α-lactalbumin gene (LALBA) reveals new alleles and polymorphisms responsible for differential expression

α-Lactalbumin (α-LA), which is encoded by the LALBA gene, is a major whey protein that binds to Ca2+ and facilitates lactose synthesis as a regulatory subunit of the synthase enzyme complex. In addition, it has been shown to play central roles in immune modulation, cell-growth regulation, and antimicrobial activity. In this study, a multitechnical approach was used to fully characterize the LALBA gene and its variants in both coding and regulatory regions for domestic camelids (dromedary, Bactrian camel, alpaca, and llama). The gene analysis revealed a conserved structure among the camelids, but a slight difference in size (2,012 bp on average) due to intronic variations. Promoters were characterized for the transcription factor binding sites (11 found in total). Intraspecies sequence comparison showed 36 SNPs in total (2 in the dromedary, none in the Bactrian camel, 22 in the alpaca, and 12 in the llama), whereas interspecies comparison showed 86 additional polymorphic sites. Eight SNPs were identified as trans-specific polymorphisms, and 2 of them (g.112A>G and g.1229A>G) were particularly interesting in the New World camels. The first creates a new binding site for transcription factor SP1. An enhancing effect of the g.112G variant on the expression was demonstrated by 3 independent pGL3 gene reporter assays. The latter is responsible for the p.78Ile>Val AA replacement and represents novel allelic variants (named LALBA A and B). A link to protein variants has been established by isoelectric focusing (IEF), and bioinformatics analysis revealed that carriers of valine (g.1229G) have a higher glycosylation rate. Genotyping methods based on restriction fragment length polymorphism (PCR-RFLP) were set up for both SNPs. Overall, adenine was more frequent (0.54 and 0.76) at both loci. Four haplotypes were found, and the AA and GA were the most common with a frequency of 0.403 and 0.365, respectively. Conversely, a putative biological gain characterizes the haplotype GG. Therefore, opportunities for rapid directional selection can be realized if this haplotype is associated with favorable milk protein properties. This study adds knowledge at the gene and protein level for α-LA (LALBA) in camelids and importantly contributes to a relatively unexplored research area in these species.

Influence of β-casein genotype on physicochemical properties and functionality of bovine milk

Several studies have been focused on the effect of milk protein genetic variants on milk physicochemical properties and functionality in recent years. β-casein, an important protein related to milk processibility, has been reported to have 2 main genetic variants A1 and A2, for which cows may be homozygous or heterozygous. In this study, several physicochemical properties of milk with β-casein variants A1A1, A1A2, and A2A2 from 3 collection occasions were analyzed. Higher manganese content and lower pH were found to be associated with the A1A1 variant compared with the other 2 genotypes. Better rennet and acid coagulation were found in A1A1 milk compared with A1A2 and A2A2 milk (although P > 0.05), whereas A2A2 milk was more stable to creaming compared with the other 2 genotypes, which may be linked to its smaller fat globule size. Thus, milk from cows with A1A1 genotype could be preferable for cheese making, while that with A2A2 variant can be used in formulations requiring good stability against creaming, and for example, yogurt making, where the softer yogurt texture may be easier to digest.

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.

Sequencing and Characterization of αs2-Casein Gene (CSN1S2) in the Old-World Camels Have Proven Genetic Variations Useful for the Understanding of Species Diversification

The CSN1S2 gene encodes αs2-casein, the third most abundant protein in camel milk. Despite its importance in foals, human nutrition, and dairy processing, the CSN1S2 gene in camels has received little attention. This study presents the first complete characterization of the CSN1S2 gene sequence in Old-World camels (Camelus bactrianus and Camelus dromedarius). Additionally, the gene promoter, consisting of 752 bp upstream of exon 1, was analyzed. The entire gene comprises 17 exons, ranging in length from 24 bp (exons 4, 8, 11, and 13) to 280 bp (exon 17). Interesting was the identification of the exon 12 in both species. The promoter analysis revealed 24 putative binding sites in the Bactrian camel and 22 in dromedary camel. Most of these sites were typical elements associated with milk protein, such as C/EBP-α, C/EBP-β, Oct-1, and AP1. The SNP discovery showed relatively high genetic diversity compared to other camel casein genes (CSN1S1, CSN2, and CSN3), with a total of 34 polymorphic sites across the two species. Particularly noteworthy is the transition g.311G>A in the CSN1S2 promoter, creating a new putative consensus binding site for a C/EBP-β in the Bactrian camel. At the exon level, two novel variants were found. One was detected in exon 6 of the Bactrian camel (g.3639C>G), resulting in an amino acid replacement, p.36Ile>Met. The second variant was found in noncoding exon 17 of dromedary CSN1S2 (g.1511G>T). Although this mutation occurs in the 3'-UnTranslated Region, it represents the first example of exonic polymorphism in the CSN1S2 for this species. This SNP also affects the binding sites of different microRNAs, including the seed sequence of the miRNA 4662a-3p, highlighting its role as a regulatory factor for CSN1S2 gene. A PCR-RFLP was set up for genotyping a dromedary Tunisian population (n = 157), and the minor allele frequency was found to be 0.27 for the G allele, indicating a potential yield improvement margin. The interspersed elements (INEs) analysis revealed 10 INEs covering 7.34% and 8.14% of the CSN1S2 sequence in the Bactrian and dromedary camels, respectively. Furthermore, six elements (A, B, F, H, I, and L) are shared among cattle and camels and are partially found in other ruminants, suggesting a common ancestral origin of these retrotransposons. Conversely, elements C, D, E, and G are specific to camels.

Difficulties in Establishing the Adverse Effects of β-Casomorphin-7 Released from β-Casein Variants-A Review

β-Casomorphin-7 (BCM-7) is a peptide released through the proteolysis of β-casein (β-CN), which is considered a bioactive peptide displaying evidence of promoting the binding and activation of the μ-opioid receptor located in various body parts, such as the gastrointestinal tract, the immune system and potentially the central nervous system. The possible effects of BCM-7 on health are a theme rising in popularity due to evidence found in several studies on the modulation of gastrointestinal proinflammatory responses that can trigger digestive symptoms, such as abdominal discomfort. With the advancement of studies, the hypothesis that there is a correlation of the possible effects of BCM-7 with the microbiota-gut-brain axis has been established. However, some studies have suggested the possibility that these adverse effects are restricted to a portion of the population, and the topic is controversial due to the small number of in vivo studies, which makes it difficult to obtain more conclusive results. In addition, a threshold of exposure to BCM-7 has not yet been established to clarify the potential of this peptide to trigger physiological responses at gastrointestinal and systemic levels. The proportion of the population that can be considered more susceptible to the effects of BCM-7 are evidenced in the literature review. The challenges of establishing the adverse effects of BCM-7 are discussed, including the importance of quantifying the BCM-7 release in the different β-CN genotypes. In summary, the reviewed literature provides plausible indications of the hypothesis of a relationship between β-CN A1/BCM-7 and adverse health effects; however, there is need for further, especially in vivo studies, to better understand and confirm the physiological effects of this peptide.

CSN1S1 and LALBA Polymorphisms and Other Factors Influencing Yield, Composition, Somatic Cell Score, and Technological Properties of Cow's Milk

We evaluated the influence of CSN1S1 and LALBA polymorphisms on cow's milk yield and quality. The analysis was done on Czech Simmental and Holstein cows. Non-genetic factors were included as well. CSN1S1 did not influence the milk performance in the first lactation. In the second lactation, cows with the BB genotype had significantly higher milk, protein, and fat yields than BC. The differences between LALBA genotypes were non-significant in the first lactation, while in the second lactation, the fat percentage was significantly higher in BB than in AB. The farm significantly influenced milk, protein, and fat yields in both the first and second lactations and fat percentage in the first lactation. The effect of CSN1S1 and LALBA genotypes on the milk technological quality was non-significant. Breed did not influence any of the evaluated technological traits and SCS. The ethanol test was not influenced by farm, season, lactation phase, protein percentage, breed, or non-fat solids percentage. Farm, season, and protein percentage significantly influenced milk fermentation ability, renneting, and SCS. The lactose content is a good indicator of udder health.

Using QTL to Identify Genes and Pathways Underlying the Regulation and Production of Milk Components in Cattle

Milk is a complex liquid, and the concentrations of many of its components are under genetic control. Many genes and pathways are known to regulate milk composition, and the purpose of this review is to highlight how the discoveries of quantitative trait loci (QTL) for milk phenotypes can elucidate these pathways. The main body of this review focuses primarily on QTL discovered in cattle (Bos taurus) as a model species for the biology of lactation, and there are occasional references to sheep genetics. The following section describes a range of techniques that can be used to help identify the causative genes underlying QTL when the underlying mechanism involves the regulation of gene expression. As genotype and phenotype databases continue to grow and diversify, new QTL will continue to be discovered, and although proving the causality of underlying genes and variants remains difficult, these new data sets will further enhance our understanding of the biology of lactation.

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.

Mechanisms and clinical management of eosinophilic oesophagitis: an overview

Since the first description of eosinophilic oesophagitis (EoE) less than three decades ago, we have observed a striking increase in the number of patients diagnosed with EoE and the understanding of its clinical and immunopathogenic background. Nonetheless, a plethora of open questions await elucidation. In this Review, we discuss the current state of knowledge regarding the underlying mechanisms, particularly environmental factors and their interaction with genetic susceptibility. Subsequently, we discuss how to translate these factors into the diagnostic and therapeutic management of this chronic, immune-mediated disorder. Finally, we dissect the still long list of unmet needs, such as reasons for and handling refractory EoE and atypical clinical presentations. These open questions can guide us through future research steps and potentially foster reconsideration of the diagnostic guidelines of EoE.

Do non-bovine domestic animals produce A2 milk?: an in silico analysis

A2 milk is an easily digestible product since it has only A2 beta-casein. In cattle, the A1 and A2 alleles are found in the population and the A2 milk is produced from A2A2 animals. Little is known about these alleles in other domestic dairy species. The present study aims to analyze sequence of genetic material available on public databases and quantify the animals genotyped. Eight domestic species were analyzed. There is strong evidence that domestic non-bovine species only carry A2 beta-casein. The data reported here for goats already confirm it due to the large number of animals genotyped as well as buffaloes. It means that they naturally produce A2 milk and no selection must be done. Thus, the fact that A2 milk is easier to digest can be used to add value to dairy product of these species. It helps to conquer new markets. It also improves people's health and breeder profitability.

Genome-wide association study for milk production and conformation traits in Canadian Alpine and Saanen dairy goats

Increasing the productivity of Canadian dairy goats is critical to the competitiveness of the sector; however, little is known about the underlying genetic architecture of economically important traits in these populations. Consequently, the objectives of this study were as follows: (1) to perform a single-step GWAS for milk production traits (milk, protein, and fat yields, and protein and fat percentages in first and later lactations) and conformation traits (body capacity, dairy character, feet and legs, fore udder, general appearance, rear udder, suspensory ligament, and teats) in the Canadian Alpine and Saanen breeds; and (2) to identify positional and functional candidate genes related to these traits. The data available for analysis included 305-d milk production records for 6,409 Alpine and 3,434 Saanen does in first lactation and 5,827 Alpine and 2,632 Saanen does in later lactations; as well as linear type conformation records for 5,158 Alpine and 2,342 Saanen does. Genotypes were available for 833 Alpine and 874 Saanen animals. Both single-breed and multiple-breed GWAS were performed using single-trait animal models. Positional and functional candidate genes were then identified in downstream analyses. The GWAS identified 189 unique SNP that were significant at the chromosomal level, corresponding to 271 unique positional candidate genes within 50 kb up- and downstream, across breeds and traits. This study provides evidence for the economic importance of several candidate genes (e.g., CSN1S1, CSN2, CSN1S2, CSN3, DGAT1, and ZNF16) in the Canadian Alpine and Saanen populations that have been previously reported in other dairy goat populations. Moreover, several novel positional and functional candidate genes (e.g., RPL8, DCK, and MOB1B) were also identified. Overall, the results of this study have provided greater insight into the genetic architecture of milk production and conformation traits in the Canadian Alpine and Saanen populations. Greater understanding of these traits will help to improve dairy goat breeding programs.

Molecular Characterization and Phylogenetic Analysis of Casein Gene Family in Camelus ferus

Camel milk is known for its exceptional medical uses. It has been used since ancient times to treat infant diarrhea, hepatitis, insulin-dependent diabetes (IDDM), lactose intolerance, alcohol-induced liver damage, allergies, and autism. It has the power to treat several diseases, with cancer being the most significant. This study investigated the evolutionary relationship, physiochemical characteristics, and comparative genomic analysis of the casein gene family (CSN1S1, CSN2, CSN1S2, and CSN3) in Camelus ferus. Molecular phylogenetics showing the camelid species clustered casein nucleotide sequences into four groups: CSN1S1, CSN2, CSN1S2, and CSN3. The casein proteins from camels were evaluated and found to be unstable, thermostable, and hydrophilic. CSN1S2, CSN2, and CSN3 were acidic, but CSN1S1 was basic. CSN1S1 showed positive selection for one amino acid (Q), CSN1S2 and CSN2 for three (T, K, Q), and CSN3 showed no positive selection. We also compared high-milk-output species such as cattle (Bos Tarus) and low-milk-yield species such as sheep (Ovies Aries) with camels (Camel ferus) and discovered that YY1 sites are more frequent in sheep than in camels and very low in cattle. We concluded that the ratio of YY1 sites in these species may affect milk production.

Does a Little Difference Make a Big Difference? Bovine β-Casein A1 and A2 Variants and Human Health-An Update

For over 20 years, bovine beta-casein has been a subject of increasing scientific interest because its genetic A1 variant during gastrointestinal digestion releases opioid-like peptide β-casomorphin-7 (β-CM-7). Since β-CM-7 is involved in the dysregulation of many physiological processes, there is a growing discussion of whether the consumption of the β-casein A1 variant has an influence on human health. In the last decade, the number of papers dealing with this problem has substantially increased. The newest clinical studies on humans showed a negative effect of variant A1 on serum glutathione level, digestive well-being, cognitive performance score in children, and mood score in women. Scientific reports in this field can affect the policies of dairy cattle breeders and the milk industry, leading to the elimination of allele A1 in dairy cattle populations and promoting milk products based on milk from cows with the A2A2 genotype. More scientific proof, especially in well-designed clinical studies, is necessary to determine whether a little difference in the β-casein amino acid sequence negatively affects the health of milk consumers.

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.

Acid and Rennet Coagulation Properties of A2 Milk

This study investigated the acid and rennet milk coagulation properties of A2 milk (β-casein (CN) A2A2 genotype), in comparison to a control milk (blend of A2A1/A1A1/A2A2 genotypes). Acid and rennet coagulation were evaluated using the Optigraph^® system, measuring the coagulation time, aggregation rate, and gel density or curd firmness. The acidification kinetics were monitored using a CINAC^® system, evaluating the time to reach pH 4.6, the acidification rate, the maximum acidification rate, the time required to reach it, and the latency time. The water-holding capacity of acid milk gels and the potential yield, total solids, and syneresis of enzymatic gels were also evaluated. Some variables were highly influenced by the farm factor, showing the importance of the effect of extrinsic parameters. Acid and enzymatic coagulation times were not affected in either milk. The A2 milk presented higher acid gel density and latency time than the control milk. Although the differences in water-holding capacity were not statistically significant, the A2 milk presented lower values, related with the higher gel density. The A2 milk also showed higher rennet aggregation rate and curd firmness than the control milk. Potential yield and syneresis were higher in the A2 milk, which is in accordance with the higher firmness of curd. Coagulation results and gel and curd properties indicate that it is possible to manufacture acid and rennet coagulation dairy products from A2 milk with no major differences when compared with a control milk.

Amyloid fibril formation by αS1- and β-casein implies that fibril formation is a general property of casein proteins

Caseins are a diverse family of intrinsically disordered proteins present in the milks of all mammals. A property common to two cow paralogues, α_S2- and κ-casein, is their propensity in vitro to form amyloid fibrils, the highly ordered protein aggregates associated with many age-related, including neurological, diseases. In this study, we explored whether amyloid fibril-forming propensity is a general feature of casein proteins by examining the other cow caseins (α_S1 and β) as well as β-caseins from camel and goat. Small-angle X-ray scattering measurements indicated that cow α_S1- and β-casein formed large spherical aggregates at neutral pH and 20°C. Upon incubation at 65°C, α_S1- and β-casein underwent conversion to amyloid fibrils over the course of ten days, as shown by thioflavin T binding, transmission electron microscopy, and X-ray fibre diffraction. At the lower temperature of 37°C where fibril formation was more limited, camel β-casein exhibited a greater fibril-forming propensity than its cow or goat orthologues. Limited proteolysis of cow and camel β-casein fibrils and analysis by mass spectrometry indicated a common amyloidogenic sequence in the proline, glutamine-rich, C-terminal region of β-casein. These findings highlight the persistence of amyloidogenic sequences within caseins, which likely contribute to their functional, heterotypic self-assembly; in all mammalian milks, at least two caseins coalesce to form casein micelles, implying that caseins diversified partly to avoid dysfunctional amyloid fibril formation.

Development of a duplex qPCR assay with locked nucleic acid probes for A, B and E kappa-casein variants detection

Milk proteins determine important milk technological characteristics. Among caseins, Ƙ-casein has been correlated with fat and protein content and cheese yield. Fourteen Ƙ-caseins variants have been described but the alleles A, B and E are the most important ones due to their frequency and/or influence on the technological aptitudes of milk. Therefore, in the present study two different duplex qPCR assays with locked nucleic acid probes (for positions 13104 and 13124 of the Ƙ-casein gene) were developed for the detection of A, B and E variants. Firstly, DNA isolation method from milk somatic cells and hair was optimised. The developed 13124-qPCR assay showed an increased sensitivity reaching up to 6.7 copies DNA copies/reaction at a 95% confidence level with A, B and E alleles reference samples. The 13104-qPCR assay reached up to 6.7 DNA copies/reaction for A allele reference sample and 67 DNA copies/reaction for B and E samples. Intra-assay variation results were below 6%. Applicability was determined using DNA samples from animals with known genotype for Ƙ-casein (AA, AB, BB, BE, AE, EE) and both assays were able to discriminate among the six genotypes with 100% accuracy. Thus, this qPCR method represents a sensitive and rapid option for the detection of Ƙ-casein alleles in both hair and milk samples.

Genetic polymorphism of Pit-1 and CSN3 genes in Holstein calves and its associations with calf birth weight

Abstract. The aim of this study was to examine the polymorphic structures of Pit-1 and CSN3 genes of Holstein calves bred in Gümüşhane province of Türkiye, to determine the distribution of genotype and allele gene frequencies, as well as examine the effects of determined polymorphisms on birth weight of calves. HinfI polymorphisms of Pit-1 and CSN3 genes were identified in DNA isolated from blood samples of 100 Holstein calves used in the study, using the PCR-RFLP method. According to the Hardy–Weinberg genetic equilibrium test, it was observed that the distribution of genotype frequencies of HinfI polymorphisms of Pit-1 genes in the studied population was in equilibrium, but not in equilibrium in terms of CSN3 gene location. The AA, AB, and BB genotype frequencies of the Pit-1 gene in the population were 13.4 %, 40.2 %, and 46.3 %, respectively; the frequency of the A allele was 0.34, while for B allele it was 0.66. The AA, AB, and BB genotype frequencies of the CSN3 gene were found to be 24.5 %, 36.7 %, and 38.8 %, respectively; the frequency of the A allele was 0.43 and the frequency of the B allele was 0.57. According to the Hardy–Weinberg genetic equilibrium test, the distribution of genotype frequencies was in equilibrium in the Pit-1/HinfI polymorphism, but not in the CSN3/HinfI polymorphism. A statistically significant relationship was not found between the genotypes of both polymorphic regions and calf birth weight.

Impact of somatic cell count combined with differential somatic cell count on milk protein fractions in Holstein cattle

Udder health in dairy herds is a very important issue given its implications for animal welfare and the production of high-quality milk. Somatic cell count (SCC) is the most widely used means of assessing udder health status. However, differential somatic cell count (DSCC) has recently been proposed as a new and more effective means of evaluating intramammary infection dynamics. Differential SCC represents the combined percentage of polymorphonuclear neutrophils and lymphocytes (PMN-LYM) in the total SCC, with macrophages (MAC) accounting for the remaining proportion. The aim of this study was to evaluate the association between SCC and DSCC and the detailed milk protein profile in a population of 1,482 Holstein cows. A validated reversed-phase HPLC method was used to quantify 4 caseins (CN), namely α_S1-CN, α_S2-CN, κ-CN, and β-CN, and 3 whey protein fractions, namely β-lactoglobulin, α-lactalbumin, and lactoferrin, which were expressed both quantitatively (g/L) and qualitatively (as a percentage of the total milk nitrogen content, %N). A linear mixed model was fitted to explore the associations between somatic cell score (SCS) combined with DSCC and the protein fractions expressed quantitatively and qualitatively. We ran an additional model that included DSCC expressed as PMN-LYM and MAC counts, obtained by multiplying the percentages of PMN-LYM and MAC by SCC for each cow in the data set. When the protein fractions were expressed as grams per liter, SCS was significantly negatively associated with almost all the casein fractions and positively associated with the whey protein α-lactalbumin, while DSCC was significantly associated with α_S1-CN, β-CN, and α-lactalbumin, but in the opposite direction to SCS. We observed the same pattern with the qualitative data (i.e., %N), confirming opposite effects of SCS and DSCC on milk protein fractions. The PMN-LYM count was only slightly associated with the traits of concern, although the pattern observed was the same as when both SCS and DSCC were included in the model. The MAC count, however, generally had a greater impact on many casein fractions, in particular decreasing both β-CN content (g/L) and proportion (%N), and exhibited the opposite pattern to the PMN-LYM count. Our results show that information obtained from both SCS and DSCC may be useful in assessing milk quality and protein fractions. They also demonstrate the potential of MAC count as a novel udder health trait.

Tunisian camel casein gene characterization reveals similarities and differences with Sudanese and Nigerian populations

Milk is a primary protein source that has always played a role in mammalian health. Despite the intensification of research projects on dromedary and the knowledge of the genetic diversity at the casein loci, the genetic structure of the Tunisian camel population still needs exploration. This study sought to determine the genetic diversity of 3 casein gene variants in 5 Tunisian camel ecotypes: c.150G>T at CSN1S1 (α_S1-casein), g.2126A>G at CSN2 (β-casein), and g.1029T>C at CSN3 (κ-casein). The obtained results were compared with data published on Sudanese and Nigerian camels to establish the level of differentiation within and between populations. A total of 159 blood samples were collected from 5 Tunisian camel ecotypes and the extracted DNA was genotyped by PCR-RFLP. A streamlined genotyping protocol was also developed for CSN3. Results indicated that allele T was quite rare (0.06) at CSN1S1 for all ecotypes. Minor allele frequency was found for G (0.462) in CSN2 except for Ardhaoui Medenine ecotype who deviated from the average CSN2 allele frequency of the total population. Allele C showed minor allele frequency of 0.384 in CSN3. Among the Tunisian population, GAT (0.343) was the most represented haplotype in all ecotypes except for Ardhaoui Medenine, where GGC (0.322) was the most frequent one. Significant differences in heterozygosity and local inbreeding were observed across the Tunisian, Sudanese, and Nigerian populations, although the global fixation index indicated that only 2.2% of the genetic variance is related to ecotype differences. Instead, phylogenetic analysis revealed a closer link between the Tunisian and Sudanese populations through a clade subdivision with 3 main branches among the ecotypes. This study represents the first attempt to understand casein gene variability in Tunisian camels; with further study, milk traits and genetic differentiation among populations can be associated with the history of camel domestication.

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.

Worldwide Research Trends on Milk Containing Only A2 β-Casein: A Bibliometric Study

Simple Summary

A1 β-casein has been correlated with adverse health outcomes, and, as a consequence, milk containing only A2 β-casein has emerged on the market. There has been a relevant increase in publications in this area since 2010. Food Science Technology and Agriculture were the main research areas of this topic. The term β-casomorphin was the most frequently used. The USA, New Zealand, and Australia were the most productive countries, though the most productive research institutions were, in absolute terms, from India, France, and Germany. The majority of the most cited studies that refer to A2 β-casein and health were reviews, and a few clinical trials have also been published.

Abstract

The protein fraction of β-casein may play a key role in the manifestation of a new intolerance: milk protein intolerance. The most common forms of β-casein among dairy cattle breeds are A1 and A2 β-casein. During gastrointestinal digestion of A1 β-casein, an opioid called peptide β-casomorphin-7 (BCM-7) is more frequently released, which can lead to adverse health outcomes. For that reason, novel products labelled as “A2 milk” or “A1-free dairy products” have appeared on the market. In this context, a bibliometric analysis on A2 β-casein research was carried out through the Web of Science (WoS) database. The main objective of this work was to provide an overview of the state of the art in the field of β-casein A2 by analyzing the number of publications per year, trends in thematic content, the most frequently used terms, and the most important institutions and countries in the field. This bibliometric study showed that a greater effort is needed to determine the possible implications of this novel product for human health and the market.

Genetic variants associated with two major bovine milk fatty acids offer opportunities to breed for altered milk fat composition

Background

Although bovine milk is regarded as healthy and nutritious, its high content of saturated fatty acids (FA) may be harmful to cardiovascular health. Palmitic acid (C16:0) is the predominant saturated FA in milk with adverse health effects that could be countered by substituting it with higher levels of unsaturated FA, such as oleic acid (C18:1cis-9). In this work, we performed genome-wide association analyses for milk fatty acids predicted from FTIR spectroscopy data using 1811 Norwegian Red cattle genotyped and imputed to a high-density 777k single nucleotide polymorphism (SNP)-array. In a follow-up analysis, we used imputed whole-genome sequence data to detect genetic variants that are involved in FTIR-predicted levels of C16:0 and C18:1cis-9 and explore the transcript profile and protein level of candidate genes.

Results

Genome-wise significant associations were detected for C16:0 on Bos taurus (BTA) autosomes 11, 16 and 27, and for C18:1cis-9 on BTA5, 13 and 19. Closer examination of a significant locus on BTA11 identified the PAEP gene, which encodes the milk protein β-lactoglobulin, as a particularly attractive positional candidate gene. At this locus, we discovered a tightly linked cluster of genetic variants in coding and regulatory sequences that have opposing effects on the levels of C16:0 and C18:1cis-9. The favourable haplotype, linked to reduced levels of C16:0 and increased levels of C18:1cis-9 was also associated with a marked reduction in PAEP expression and β-lactoglobulin protein levels. β-lactoglobulin is the most abundant whey protein in milk and lower levels are associated with important dairy production parameters such as improved cheese yield.

Conclusions

The genetic variants detected in this study may be used in breeding to produce milk with an improved FA health-profile and enhanced cheese-making properties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00731-9.

Role of CSN2, CSN3, and BLG genes and the polygenic background in the cattle milk protein profile

To devise better selection strategies in dairy cattle breeding programs, a deeper knowledge of the role of the major genes encoding for milk protein fractions is required. The aim of the present study was to assess the effect of the CSN2, CSN3, and BLG genotypes on individual protein fractions (α_S1-CN, α_S2-CN, β-CN, κ-CN, β-LG, α-LA) expressed qualitatively as percentages of total nitrogen content (% N), quantitatively as contents in milk (g/L), and as daily production levels (g/d). Individual milk samples were collected from 1,264 Brown Swiss cows reared in 85 commercial herds in Trento Province (northeast Italy). A total of 989 cows were successfully genotyped using the Illumina Bovine SNP50 v.2 BeadChip (Illumina Inc.), and a genomic relationship matrix was constructed using the 37,519 SNP markers obtained. Milk protein fractions were quantified and the β-CN, κ-CN, and β-LG genetic variants were identified by reversed-phase HPLC (RP-HPLC). All protein fractions were analyzed through a Bayesian multitrait animal model implemented via Gibbs sampling. The effects of days in milk, parity order, and the CSN2, CSN3, and BLG genotypes were assigned flat priors in this model, whereas the effects of herd and animal additive genetic were assigned Gaussian prior distributions, and inverse Wishart distributions were assumed for the respective co-variance matrices. Marginal posterior distributions of the parameters of interest were compared before and after the inclusion of the effects of the 3 major genes in the model. The results showed that a high portion of the genetic variance was controlled by the major genes. This was particularly apparent in the qualitative protein profile, which was found to have a higher heritability than the protein fraction contents in milk and their daily yields. When the genes were included individually in the model, CSN2 was the major gene controlling all the casein fractions except for κ-CN, which was controlled directly by the CSN3 gene. The BLG gene had the most influence on the 2 whey proteins. The genetic correlations showed the major genes had only a small effect on the relationships between the protein fractions, but through comparison of the correlation coefficients of the proteins expressed in different ways they revealed potential mechanisms of regulation and competitive synthesis in the mammary gland. The estimates for the effects of the CSN2 and CSN3 genes on protein profiles showed overexpression of protein synthesis in the presence of the B allele in the genotype. Conversely, the β-LG B variant was associated with a lower concentration of β-LG compared with the β-LG A variant, independently of how the protein fractions were expressed, and it was followed by downregulation (or upregulation in the case of the β-LG B) of all other protein fractions. These results should be borne in mind when seeking to design more efficient selection programs aimed at improving milk quality for the efficiency of dairy industry and the effect of dairy products on human health.

Comparative Analysis of the Polymorphism of the Casein Genes in Camels Bred in Kazakhstan

Caseins play an important role in determining the technological properties and quantitative characteristics of camel milk. To date, only a few studies on the genetic polymorphism of casein genes have been reported in the camel populations of Kazakhstan. Therefore, this work aimed to identify the genetic polymorphism level of casein genes among camel populations of the Almaty region of Kazakhstan. The PCR-RFLP method was used for this purpose and the following genotypes were revealed as a result: CSN3 gene—CC, CT, TT, where the T allele predominated in all populations, with a frequency of 0.60; CSN2 gene—AA, AG, GG, with the predomination of A allele (0.64); and CSN1S1 gene—GG and GT, with the predomination of G allele (0.94). Statistical analysis was carried out using the POPGENE and GenAlEx software. The χ2 values were equal to 12.1 (CSN3), 8.6 (CSN2), and 14.5 (CSN1S1). As a result, three out of 53 animals were designated as the “core” of the population—animals with the desired genotypes: CC genotype for the CSN3 gene and AA genotype for the CSN2 gene. Such animals can be selected for further use with an increase in the number of livestock with high productivity rates.

Screening for phenotypic outliers identifies an unusually low concentration of a β-lactoglobulin B protein isoform in bovine milk caused by a synonymous SNP

Background

Milk samples from 10,641 dairy cattle were screened by a mass spectrometry method for extreme concentrations of the A or B isoforms of the whey protein, β-lactoglobulin (BLG), to identify causative genetic variation driving changes in BLG concentration.

Results

A cohort of cows, from a single sire family, was identified that produced milk containing a low concentration of the BLG B protein isoform. A genome-wide association study (GWAS) of BLG B protein isoform concentration in milk from AB heterozygous cows, detected a group of highly significant single nucleotide polymorphisms (SNPs) within or close to the BLG gene. Among these was a synonymous G/A variation at position + 78 bp in exon 1 of the BLG gene (chr11:103256256G > A). The effect of the A allele of this SNP (which we named B’) on BLG expression was evaluated in a luciferase reporter assay in transfected CHO-K1 and MCF-7 cells. In both cell types, the presence of the B’ allele in a plasmid containing the bovine BLG gene from -922 to + 898 bp (relative to the transcription initiation site) resulted in a 60% relative reduction in mRNA expression, compared to the plasmid containing the wild-type B sequence allele. Examination of a mammary RNAseq dataset (n = 391) identified 14 heterozygous carriers of the B’ allele which were homozygous for the BLG B protein isoform (BB’). The level of expression of the BLG B’ allele was 41.9 ± 1.0% of that of the wild-type BLG B allele. Milk samples from three cows, homozygous for the A allele at chr11:103,256,256 (B’B’), were analysed (HPLC) and showed BLG concentrations of 1.04, 1.26 and 1.83 g/L relative to a mean of 4.84 g/L in milk from 16 herd contemporaries of mixed (A and B) BLG genotypes. The mechanism by which B’ downregulates milk BLG concentration remains to be determined.

Conclusions

High-throughput screening and identification of outliers, enabled the discovery of a synonymous G > A mutation in exon 1 of the B allele of the BLG gene (B’), which reduced the milk concentration of β-lactoglobulin B protein isoform, by more than 50%. Milk from cows carrying the B’ allele is expected to have improved processing characteristics, particularly for cheese-making.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00711-z.

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.