- Research article
- Open Access
Comparative milk proteome analysis of Kashmiri and Jersey cattle identifies differential expression of key proteins involved in immune system regulation and milk quality
BMC Genomics volume 21, Article number: 161 (2020)
Exploration of the bioactive components of bovine milk has gained global interest due to their potential applications in human nutrition and health promotion. Despite advances in proteomics profiling, limited studies have been carried out to fully characterize the bovine milk proteome. This study explored the milk proteome of Jersey and Kashmiri cattle at day 90 of lactation using high-resolution mass spectrometry based quantitative proteomics nano-scale LC-MS/Q-TOF technique. Data are available via ProteomeXchange with identifier PXD017412.
Proteins from whey were fractionated by precipitation into high and low abundant proteins. A total of 81 high-abundant and 99 low-abundant proteins were significantly differentially expressed between Kashmiri and Jersey cattle, clearly differentiating the two breeds at the proteome level. Among the top differentiating proteins, the Kashmiri cattle milk proteome was characterised by increased concentrations of immune-related proteins (apelin, acid glycoprotein, CD14 antigen), neonatal developmental protein (probetacellulin), xenobiotic metabolising enzyme (flavin monooxygenase 3 (FMO3), GLYCAM1 and HSP90AA1 (chaperone) while the Jersey milk proteome presented higher concentrations of enzyme modulators (SERPINA1, RAC1, serine peptidase inhibitor) and hydrolases (LTF, LPL, CYM, PNLIPRP2). Pathway analysis in Kashmiri cattle revealed enrichment of key pathways involved in the regulation of mammary gland development like Wnt signalling pathway, EGF receptor signalling pathway and FGF signalling pathway while a pathway (T-cell activation pathway) associated with immune system regulation was significantly enriched in Jersey cattle. Most importantly, the high-abundant FMO3 enzyme with an observed 17-fold higher expression in Kashmiri cattle milk seems to be a characteristic feature of the breed. The presence of this (FMO3) bioactive peptide/enzyme in Kashmiri cattle could be economically advantageous for milk products from Kashmiri cattle.
In conclusion, this is the first study to provide insights not only into the milk proteome differences between Kashmiri and Jersey cattle but also provides potential directions for application of specific milk proteins from Kashmiri cattle in special milk preparations like infant formula.
Bovine milk is a valued natural product which delivers a matrix of essential nutrients including growth and immune factors to offspring and a key raw material for human food preparations [1, 2]. Some studies have characterized the bovine milk proteome, its bioactive profile, and the extent of cross reactivity of bovine bioactive milk peptides on various biological functions [3,4,5,6,7]. Milk proteins are generally categorized into three major groups: caseins, whey proteins and milk fat globule membrane proteins [4, 8]. Most of the polypeptides in milk are an essential source of amino acids to neonates  and many resist proteolysis [10, 11]. Milk peptides also facilitate absorption of other nutrients in the gastro-intestinal tract, provide humoral immune responses and support intestinal development . Besides, digestion or fermentation of milk proteins also produces a number of bioactive peptides, which contribute as well to the various functional properties of milk [13, 14]. The major proteins in milk are far outnumbered by numerous other minor proteins which play important roles in a wide range of physiological activities including antioxidant activity, post-natal development of new-borns, maturation of the immune system, establishment of symbiotic microflora, and protection against various pathogens [15, 16].
Several studies have characterised the milk proteome in different species and breeds using different quantitative proteomic techniques [7, 16,17,18,19,20]. The differences in the milk proteome profile have been attributed to genetic, management and disease factors [7, 21]). Although the diverse composition and biological functions of bovine milk has been reported extensively [22,23,24], the comparative abundance of milk proteins in Indian cattle breeds have not been investigated till date. Kashmiri and Jersey cattle are two important milk animals which contribute significantly to the total milk production in the Indian northern state of Kashmir. The Kashmiri cattle is an indigenous breed kept mainly for milk production in the hilly regions of Kashmir. Kashmiri cattle are small, hardy and adapted to the hilly regions of Kashmir. Whereas, Jersey is a well-established dairy breed imported to augment the milk production ability of Kashmiri cattle through cross breeding. We hypothesize that the proteome profile of Kashmiri cattle milk may have special properties or differ from that of the well-established Jersey dairy breed due to its different genetic background and milk producing ability. Therefore, the aim of this study was to study the protein profiles of Kashmiri and Jersey cattle milk which could reveal important protein factors underlying the physiological differences and differences in milk traits between the two breeds.
Proteome profile of bovine milk
Proteins from whey were fractionated by precipitation into high and low abundant proteins. A total of 180 proteins were differentially expressed (DE) (FDR < 0.1) between Kashmiri and Jersey cattle. Specifically, 91 and 89 proteins were significantly upregulated (FDR < 0.1) in Kashmiri and Jersey cattle, respectively (Additional file 2: Table S2a and S2b, Additional file 3). The most upregulated high abundant proteins (fold change (FC) > 2) were CSN2, CD4 and LF, and low abundant proteins were FMO3, GLYCAM1, APLN and BTC in Kashmiri cattle (Table 1, Fig. 1). Whereas, LALBA, ZNF496, CSN3 and LGB were the most upregulated high abundant proteins and RAC1, B2M and SAR1B were the most upregulated minor milk proteins in Jersey cattle (Table 1).
Enriched gene ontology terms of significantly upregulated proteins in Kashmiri and Jersey cattle
Gene ontology (GO) enrichment of significantly upregulated proteins in Kashmiri and Jersey cattle found a total of 4 enriched GO terms in Kashmiri and 4 in Jersey cattle (Table 2). Only extracellular region (GO:0005576) reached significance after FDR correction in both breeds (Table 2).
Protein categories identified through GO annotation
The identified differentially upregulated proteins in Kashmiri and Jersey cattle were categorized according to their GO annotation (Additional file 2: Table S103). Most of the significantly upregulated proteins in both cattle breeds were enzyme modulators (SERPINA3, BTN1A1, SERPINC1, SERPINF2, Serin peptidase inhibitor, RAC1, RRAS, BTN1A1 and uterine milk protein) and hydrolases (GNB2, CTSD, GNB1, PNLIPRP2, CYM) (Fig. 1 a and b). However, proteins belonging to the chaperone classes (HSP90AA1, YWHAB, YWHAZ) were significantly upregulated in Kashmiri cattle only (Fig. 2a and b).
Enriched pathways by significantly upregulated proteins in Kashmiri and Jersey cattle
Significantly upregulated proteins in Kashmiri and Jersey cattle were enriched to 12 and 4 pathways at uncorrected P < 0.05, respectively (Table 3). When FDR correction was applied, 10 and one proteins remained significant (FDR < 0.1) in Kashmiri and Jersey cattle, respectively (Table 3). Of all the pathways, only EGF receptor signalling pathway was enriched at uncorrected P < 0.05 by significantly upregulated proteins in both breeds.
The present study was designed to characterize and compare the milk proteome of Kashmiri and Jersey cattle. Over the past few decades, interest to reveal the dynamics of milk proteome has grown and there have been remarkable developments in the techniques used for fractionation and identification of proteins [25,26,27]. In the present study, a combination of fractionation and mass spectrometry techniques were used to comprehensively characterize the milk proteome profiles of Kashmiri and Jersey cattle breeds.
A total of 180 proteins were found to be differentially expressed between Kashmiri and Jersey cattle. Interestingly, 90 and 89 of the differentially expressed proteins were significantly upregulated in Kashmiri and Jersey cattle, respectively. Enzyme modulators were the major class of up-regulated proteins in both Kashmiri (20.51%) and Jersey cattle (14.28%). Hydrolases represented 12.82 and 14.28% of upregulated proteins in Kashmiri and Jersey cattle, respectively. Interestingly, chaperone class of proteins was only observed in milk of Kashmiri cattle. Chaperones help in the folding of newly synthesized proteins and prevent their premature (mis) folding at least until a domain capable of forming a stable structure is synthesized. As expected and in agreement with earlier studies ([26, 27]), the casein and whey fraction proteins were highly expressed in both breeds. However, a different set of high abundant milk proteins were significantly upregulated in each of the breeds. For example, the abundantly expressed proteins beta-casein, lactoferrin and CD4 were significantly upregulated in Kashmiri while beta-lacto globulin, kappa-casein and alpha-lactalbumin were significantly upregulated in Jersey (Table 1). Interestingly, the low abundant proteins FMO3, GLYCAM1, CD9, APLN, BTC, enterotoxin-binding glycoprotein PP16K, ORM1, serin peptidase inhibitor clade A, adipocyte differentiation-related protein and uterine milk protein were significantly upregulated in Kashmiri while ATP synthase subunit A, RAC1, B2M, SAR1B, TCN2 and MFGE8 were upregulated in Jersey. These results indicate a clear distinction as well as wide differences in the proteome profiles between the breeds which could be explained by high selection pressure for milk production traits in Jersey.
The differences in the expression of high abundant proteins between the breeds might confer differential benefits to their milks. For example, different levels of phosphorylation of beta-casein has been reported to affect the availability of calcium and protein micelle stability of milk , which might have important consequences on the nutrition and technological properties of milk and milk products. Additionally, other key bioactive proteins identified in this study that are well known to exert beneficial effects on human nutrition and health include lactoferrin, GLYCAM1, betacellulin, apelin, LALBA and serine peptidase inhibitor, etc. Iron sequestering properties of lactoferrin (LF), along with blockade of microbial carbohydrate metabolism and destabilisation of the bacterial cell wall [29, 30], has been shown to produce bactericidal and bacteriostatic effects in a wide range of microorganisms, including gram positive and gram negative bacteria, aerobes, anaerobes, yeasts and parasites [31,32,33]. Similarly, GLYCAM1 with a 7.93-fold expression in Kashmiri cattle is known to act as an antimicrobial peptide with ability to protect the intestinal mucosal tract of neonates largely due to its lubricating properties [34, 35]. In addition to these, apelin peptides might be involved in maturation of the gastrointestinal tract [36, 37]. Betacellulin (BTC), a key epidermal growth factor (EGF)  might regulate the development and maturation of the neonatal gut and immune system . EGFs are major growth promoting factors in human milk  but the biological significance of BTC in bovine milk is currently unclear and needs further investigation. However, one plausible explanation for the presence of BTC in bovine milk might be to stimulate the proliferation of the gastrointestinal epithelia in new-borns, as has been proposed for milk-borne EGF and TGF-α (Transforming growth factor alpha) in other species . With respect to Jersey breed, peptides resulting from partial digestion of high abundant proteins such as LALBA, CSN2 and CSN3 in the small intestine may influence gut functions including immune stimulation, mineral and trace element absorption and host defence against infection . Alpha-lactalbumin enhances infant gastrointestinal function , motility and antimicrobial activity . CSN3 is readily hydrolysed in calf’s stomach, allowing the formation of a coagulum that can be readily digested  and also provides heat stability to milk by stabilising the casein micelle . Moreover, CSN3 prevents infection by disrupting the attachment of pathogens to mucosal cells . CSN3 digestion results in the formation of a glycomacropeptide which in turn enhances mineral absorption . Bovine beta 2-microglobulin (B2M) is an antibacterial protein present in milk fat globules. B2M possesses potent antibacterial activities against Gram positive pathogenic bacteria . Bovine milk is an abundant source of bioavailable B12 vitamin wherein when complexed with transcobalamin, a major vitamin B12 binding protein in cows’ milk , stimulates vitamin B12 absorption through intestinal epithelial cells . Lactadherin is secreted by mammary epithelial cells and stored in milk fat globules . Lactadherin, as one of the immune components in bovine milk has been found to prevent rota viral infection in infants by removing the sialic acid from the viral coat [52, 53].
It is worthwhile to note that the low abundant protein, flavin-containing monooxygenase 3 (FMO3) had 16.6 fold expression rate in Kashmiri as compared to Jersey. This is the first report wherein FMO3 has been found to be highly expressed in Kashmiri cattle. Increased presence of FMO3 might be important due to its ability to oxidise trimethylamine (TMA), a compound with fishy odour, to TMAO (Trimethylamine N-oxide), an odourless oxide. Absence of FMO3 leads to fishy flavour in milk due to increased build-up of TMA, and thus might play an important role in maintaining the quality of milk [54,55,56]. Moreover, FMO3 belongs to a drug metabolising enzyme class with ability to oxidize xenobiotics, pesticides and other foreign inhabitants in body fluids including milk and serum [57,58,59,60] and hence presents an efficient defence mechanism in new-borns. The presence of FMO3 at high concentrations in Kashmiri cattle milk can favour utilization of Kashmiri cattle milk in commercial preparations for the promotion of human health and nutritional status. In fact, bio-mining of such bioactive milk protein constituent and marketing it as ingredients may not only serve as a lucrative business for the Indian dairy industry but also in the development of products for consumers with special needs like allergy and milk tolerance.
The GO analysis of significantly up-regulated proteins revealed only one significantly enriched GO term (extracellular region) after FDR correction in both breeds and limited functional overlap was found between the present proteomic data and our earlier transcriptome data  indicating the failure of RNA-based analyses to represent completely protein dynamics .
Pathway analysis helps in biological interpretation of proteomic and other high-throughput data in cells or organisms . Most of the pathways (Wnt signaling pathway, EGF receptor signaling pathway, FGF signaling pathway, PI3 kinase pathway) significantly enriched by the significantly upregulated proteins in Kashmiri cattle are involved in mammary gland development. Wnt signaling pathway regulates mammary development  during various stages of mammary morphogenesis . The proteins enriched in the Wnt signalling pathway were GNB1(G protein subunit beta 1), GNB2 (G protein subunit bBeta 2) and ACTG1(actin gamma 1). ACTG1 plays a critical role in branching and alveolar development of the mammary gland through cytoskeletal remodelling . FGF signalling pathway controls mammary epithelial cell branching and morphogenesis  and activates PI3 kinase pathway through phosphorylation . Epidermal growth factor family plays essential roles in regulating cell proliferation, survival and differentiation of mammary epithelial cells through STAT5A, a key non-tyrosine kinase protein indirectly regulated by JAK2/ELF5, insulin growth factor, estrogen, and progesterone signalling pathways . In Jersey cattle, two significantly (p < 0.05) enriched pathways, blood coagulation/coagulation cascades and T cell activation pathways are associated with immune system regulation . SERPINA1, SERPINC1, SERPINF2 are important proteins in blood coagulation pathway whereas, B2M and RAC1 play critical roles in T cell activation pathway. These proteins play fundamental roles in innate immunity in addition to enhancing adaptive immune responses . Altogether, a wide range of proteins were detected in this study including proteins involved in immune response, host defense and milk quality as well as qualitative and quantitative differences in their milk proteome.
A total of 91 and 89 proteins were significantly upregulated in Kashmiri and Jersey cattle, respectively. A different set of high- abundant and low-abundant proteins were significantly upregulated in Kashmiri and Jersey cattle, clearly differentiating the two breeds at the proteome level. Immune-related proteins (CD4, LF and GLYCAM 1) and drug metabolising enzyme (FMO3) were abundantly expressed in Kashmiri cattle milk. The presence of FMO3 at high concentrations in Kashmiri cattle milk could favour its utilization in commercial preparations for human health promotion and consequently serve as a boost for increased business opportunities for the Indian dairy industry.
Experimental animals and sampling
The ethical clearance was approved by the Institutional Animal Ethics Committee (IAEC) of Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir. A total of three healthy Kashmiri and three Jersey cows in their 3rd lactation from the university dairy farm (Mountain Livestock Research Institute, Share-Kashmir University of Agricultural Sciences and Technology of Kashmir, India) were selected for the study. The animals were kept under similar feeding and management conditions to minimise environmental variation. Fresh milk samples (200 mL) were aseptically collected from all the four quarters (50 mL per quarter) at day 90 in milk (D90), mixed thoroughly, placed on ice and immediately transported to the laboratory for further analysis.
Milk samples were processed differently for high and low abundance protein analysis. For high-abundance protein analysis, 50 mL of milk was immediately placed on ice after collection followed by centrifugation at 4000×g for 10 min at 4 °C within 2 h of collection. The fat layer was removed and skimmed fraction was stored at − 20 °C. Whereas, for low abundance protein analysis, 0.24 mL (100X) mammalian protease inhibitor cocktail (Sigma, Milwaukee, WI, USA) was added to 50 mL of milk followed by centrifugation at 4000×g for 15 min at 4 °C. The cream layer was removed and the skimmed or whey portion was depleted of casein using a previously described method . Briefly, 60 mM CaCl2 was added to skimmed sample and the pH was adjusted to 4.3 using 30% acetic acid (Fisher Scientific, Fair Lawn, NJ, USA). Samples were then centrifuged at 189,000×g at 4 °C for 70 min and the supernatant was collected and stored at − 80 °C.
Enrichment of low abundance proteins
Low abundance minor proteins were enriched using the ProteoMiner Kit (BioRad Laboratories, Hercules, CA, USA) as per manufacturer’s protocol. Whey samples were placed in individual ProteoMiner columns, mixed thoroughly by shaking (gently) followed by incubation at room temperature for 2 h. Subsequently, samples were washed thoroughly using HPLC grade water to remove excess proteins by centrifugation at 7000 g for 5 min. Low abundance proteins were eluted off the beads by addition of 20 μl 4 x Laemmli sample buffer (8% SDS, 40% glycerol, 250 mM Tris, pH 6.8, 400 mM DTT with trace amount of bromophenol blue).
In-solution digestion of proteins and nano-scale LC/MS analysis on QTOF
The pellets after acetone precipitation (high abundant proteins) or TCA (Trichloroacetic acid)-acetone precipitation (low abundant proteins) were dissolved in 50 mM ammonium bicarbonate (dilution 1:3) and 0.1% SDS. 100 μg of the extracted protein was subjected to in solution trypsin digestion with carbamidomethylation at cysteine (fixed) and oxidation at methionine (variable). The dissolved pellet was treated with 10 μl of 100 mM DTT (Dithiothreitol) followed by incubation on a thermo mixer (Eppendorf ThermoMixer® C,) at 95 °C for 1 h. The sample was treated with 18 μl of 250 mM IDA (Iodoacetamide) and then incubated in the dark for 45 min at room temperature. To stop the IDA reaction, 40 μl DTT was added at room temperature and incubated for 10 min. To this solution, 50 mM ammonium bicarbonate and 0.1% SDS was added to make up the volume to 300 μl. For Enzymatic cleavage of the protein, trypsin in the ratio 50:1 (w/v) was added to sample and incubated on the thermo mixer at 37 °C overnight. To stop the trypsin activity, the peptides were then extracted in 0.1% formic acid followed by incubation at 37 °C for 45 min. The extracted mixture was then centrifuged at 13000 g for 10 min and the supernatant was placed in a separate Eppendorf tube. This supernatant was subjected to speed vac at 45 °C. The resulting peptides were then dissolved in 20 μl of 0.1% formic acid and 10 μL of this solution was used on C18 UPLC column for separation of peptides. The mass spectrometer was operated in positive ion mode, and MS spectra were acquired over a range of 375–1500 m/z. For MS and MS/MS scans, the resolution of the orbitrap fusion was set at 120,000 and 50,000 at 200 m/z, respectively. Data-dependent acquisition mode was set as top speed, and ions were fragmented (10 fragment files collected after every full scan) through higher energy collisional dissociation, and cycle time was 3 s with peptide mass tolerance and fragment mass tolerance of 50 ppm and 100 ppm, respectively. The automatic gain control target values for master scan modes and MS/MS were set to 4e5 and 1e5, respectively. Dynamic exclusion duration was 40 s.
Protein identification and differential expression analysis
The individual peptides MSMS spectra were searched against the Swiss-Prot databases using the Mascot Distiller Search engine (v. 2.6.0) for protein identification and expression analysis was performed with PLGS software (Protein Lynx Global Server, Waters, India) by Sandor’s Lifesciences, Hyderabad, India. The results were filtered based on the peptide Benjaminin and Hochberg corrected p-value < 0.1 (FDR < 0.1) or uncorrected p-value < 0.05. Both unique and razor peptides were selected for protein quantification, protein ratios were calculated as the median of only unique or razor peptides of the protein. All peptide ratios were normalized based on the median ratio. The protein species quantification results were statistically analysed by student’s t-test, and the p-value was corrected by the method of Benjamin and Hochberg FDR analysis. An FDR < 0.1 was considered significant due to the low number of samples analysed.
Gene ontology and pathway analysis
Gene ontology (GO) and pathway enrichment analysis of differentially expressed proteins was accomplished with Gene Ontology Consortium data base (http://www.geneontology.org) (Falcon and Gentleman, 2007). GO terms and KEGG pathways (http://www.genome.jp/kegg/) with FDR < 0.1were considered significantly enriched.
Availability of data and materials
The datasets generated and analysed during the current study are available as Additional files.
Epidermal growth factor
Early growth response protein 1
EH domain-containing protein 1
False discovery rate
Fibroblast growth factor
Glycosylation-dependent cell adhesion molecule 1
Heat shock protein90AA1
Liquid chromatography-mass spectrometry/quantitative time of flight
Mammary epithelial cell
Pancreatic lipase related protein 2
Ras-related C3 botulinum toxin substrate 1
Serine protease inhibitor1
Transforming growth factor
Toll like receptor 2
Zinc finger protein 496
Reinhardt TA, Lippolis JD. Bovine milk fat globule membrane proteome. J Dairy Res. 2006;73(4):406–16.
Séverin S, Wenshui X. Milk biologically active components as nutraceuticals. Crit Rev Food Sci Nutr. 2005;45(7–8):645–56.
Dai W, Chen Q, Wang Q, White RR, Liu J, Liu H. Complementary transcriptomic and proteomic analyses reveal regulatory mechanisms of milk protein production in dairy cows consuming different forages. Sci Rep. 2017;7:44234.
Lönnerdal B. Infant formula and infant nutrition: bioactive proteins of human milk and implications for composition of infant formulas. Am J Clin Nutr. 2014;99(3):712S–7S.
Liao Y, Jiang R, Lönnerdal B. Biochemical and molecular impacts of lactoferrin on small intestinal growth and development during early life. Biochem Cell Biol. 2012;90(3):476–84.
Reznikov EA, Comstock SS, Yi C, Contractor N, Donovan SM. Dietary Bovine Lactoferrin Increases Intestinal Cell Proliferation in Neonatal Piglets, 2. J Nutr. 2014;144(9):1401–8.
Ibeagha-Awemu EM, Ibeagha AE, Messier S, Zhao X. Proteomics, genomics, and pathway analyses of Escherichia coli and Staphylococcus aureus infected milk whey reveal molecular pathways and networks involved in mastitis. J Proteome Res. 2010;9(9):4604–19.
Wada Y, Lönnerdal B. Bioactive peptides derived from human milk proteins—mechanisms of action. J Nutr Biochem. 2014;25(5):503–14.
Casado B, Affolter M, Kussmann M. OMICS-rooted studies of milk proteins, oligosaccharides and lipids. J Proteomics. 2009;73(2):196–208.
German JB, Dillard CJ, Ward RE. Bioactive components in milk. Curr Opin Clin Nutr Metab Care. 2002;5(6):653–8.
Lönnerdal B. Nutritional and physiologic significance of human milk proteins. Am J Clin Nutr. 2003;77(6):1537S–43S.
Mohanty DP, Mohapatra S, Misra S, Sahu PS. Milk derived bioactive peptides and their impact on human health–a review. Saudi J Biol Sci. 2016;23(5):577–83.
Meisel H. Biochemical properties of peptides encrypted in bovine milk proteins. Curr Med Chem. 2005;12(16):1905–19.
Ibeagha-Awemu EM, Liu J, Zhao X. Bioactive components in yogurt products. In: Park YW, editor. Bioactive components in Milk and dairy products, 2009; 2009. p. 235–50.
Ward RE, German JB. Understanding milk’s bioactive components: a goal for the genomics toolbox. J Nutr. 2004;134(4):962S–7S.
Tacoma R, Fields J, Ebenstein DB, Lam YW, Greenwood SL. Characterization of the bovine milk proteome in early-lactation Holstein and Jersey breeds of dairy cows. J Proteomics. 2016;130:200–10.
Litwińczuk Z, Król J, Brodziak A, Barłowska J. Changes of protein content and its fractions in bovine milk from different breeds subject to somatic cell count. J Dairy Sci. 2011;94(2):684–91.
D'auria E, Agostoni C, Giovannini M, Riva E, Zetterström R, Fortin R, Roncada P. Proteomic evaluation of milk from different mammalian species as a substitute for breast milk. Acta Paediatr. 2005;94(12):1708–13.
Yang Y, Bu D, Zhao X, Sun P, Wang J, Zhou L. Proteomic analysis of cow, yak, buffalo, goat and camel milk whey proteins: quantitative differential expression patterns. J Proteome Res. 2013;12(4):1660–7.
Wang X, Zhao X, Huang D, Pan X, Qi Y, Yang Y, Cheng G. Proteomic analysis and cross species comparison of casein fractions from the milk of dairy animals. Sci Rep. 2017;7:43020.
Oltenacu PA, Algers B. Selection for increased production and the welfare of dairy cows: are new breeding goals needed? Ambio. 2005;34(4):311–5.
Zhang L, van Dijk AD, Hettinga K. An interactomics overview of the human and bovine milk proteome over lactation. Proteome Sci. 2016;15(1):1.
Mol P, Kannegundla U, Dey G, Gopalakrishnan L, Dammalli M, Kumar M, TSK P. Bovine Milk Comparative Proteome Analysis from Early, Mid, and Late Lactation in the Cattle Breed, Malnad Gidda (Bos indicus). Omics. 2018;22(3):223–35.
Affolter M, Grass L, Vanrobaeys F, Casado B, Kussmann M. Qualitative and quantitative profiling of the bovine milk fat globule membrane proteome. J Proteome. 2010;73(6):1079–88.
Murakami K, Lagarde M, Yuki Y. Identification of minor proteins of human colostrum and mature milk by two-dimensional electrophoresis. Electrophoresis. 1998;19(14):2521–7.
Yamada M, Murakami K, Wallingford JC, Yuki Y. Identification of low-abundance proteins of bovine colostral and mature milk using two-dimensional electrophoresis followed by microsequencing and mass spectrometry. Electrophoresis. 2002;23(7–8):1153–60.
Palmer DJ, Kelly VC, Smit AM, Kuy S, Knight CG, Cooper GJ. Human colostrum: identification of minor proteins in the aqueous phase by proteomics. Proteomics. 2006;6(7):2208–16.
Amigo L, Recio I, Ramos M. Genetic polymorphism of ovine milk proteins: its influence on technological properties of milk—a review. Int Dairy J. 2000;10(3):135–49.
Masson PL, Heremans JF. Lactoferrin in milk from different species. Comp Biochem Physiol. 1971;1:119–29.
Hennart PF, Brasseur DJ, Delogne-Desnoeck JB, Dramaix MM, Robyn CE. Lysozyme, lactoferrin, and secretory immunoglobulin A content in breast milk: influence of duration of lactation, nutrition status, prolactin status, and parity of mother. Am J Clin Nutr. 1991;53(1):32–9.
Weinberg ED. Iron withholding: a defense against infection and neoplasia. Physiol Rev. 1984;64(1):65–102.
Arnold RR, Russell JE, Champion WJ, Brewer M, Gauthier JJ. Bactericidal activity of human lactoferrin: differentiation from the stasis of iron deprivation. Infect Immun. 1982;35(3):792–9.
Ellison RT III, LaForce FM, Giehl TJ, Boose DS, Dunn BE. Lactoferrin and transferrin damage of the gram-negative outer membrane is modulated by Ca2+ and Mg2+. Microbiology. 1990;136(7):1437–46.
Carraway KL, Hull SR. Cell surface mucin-type glycoproteins and mucin-like domains. Glycobiology. 1991;1(2):131–8.
Larocca D, Peterson JA, Walkup G, Urrea R, Ceriani RL. Cloning and sequencing of a complementary DNA encoding a Mr 70,000 human breast epithelial mucin-associated antigen. Cancer Res. 1990;50(18):5925–30.
Wang G, Kundu R, Han S, Qi X, Englander EW, Quertermous T, Greeley GH Jr. Ontogeny of apelin and its receptor in the rodent gastrointestinal tract. Regul Pept. 2009;158(1–3):32–9.
Aydin S. The presence of the peptides apelin, ghrelin and nesfatin-1 in the human breast milk, and the lowering of their levels in patients with gestational diabetes mellitus. Peptides. 2010;31(12):2236–40.
Dunbar AJ, Priebe IK, Belford DA, Goddard C. Identification of betacellulin as a major peptide growth factor in milk: purification, characterization and molecular cloning of bovine betacellulin. Biochem J. 1999;344(Pt 3):713.
Grosvenor CE, Picciano MF, Baumrucker CR. Hormones and growth factors in milk. Endocrine Rev. 1993;14(6):710–28.
Carpenter G. Epidermal growth factor is a major growth-promoting agent in human milk. Science. 1980;210(4466):198–9.
Ceciliani F, Pocacqua V, Provasi E, Comunian C, Bertolini A, Bronzo V, Sartorelli P. Identification of the bovine $\alpha $1-acid glycoprotein in colostrum and milk. Vet Res. 2005;36(5–6):735–46.
Lönnerdal B. Human milk: bioactive proteins/peptides and functional properties. In Protein in Neonatal and Infant Nutrition: Recent Updates,2016, (Vol. 86, pp. 97–107). Karger Publishers.
Berger M, Gray JA, Roth BL. The expanded biology of serotonin. Ann Rev Med. 2009;60:355–66.
Bounous G, Kongshavn PA. Influence of dietary proteins on the immune system of mice. J Nutr. 1982;112(9):1747–55.
Creamer LK, Plowman JE, Liddell MJ, Smith MH, Hill JP. Micelle stability: κ-casein structure and function. J Dairy Sci. 1998;81(11):3004–12.
Strömqvist M, Falk P, Bergström S, Hansson L, Lönnerdal B, Normark S, Hernell O. Human milk kappa-casein and inhibition of Helicobacter pylori adhesion to human gastric mucosa. J Pediatr Gastroenterol Nutr. 1995;21(3):288–96.
Kelleher SL, Chatterton D, Nielsen K, Lönnerdal B. Glycomacropeptide and α-lactalbumin supplementation of infant formula affects growth and nutritional status in infant rhesus monkeys. Am J Clin Nutr. 2003;77(5):1261–8.
Pringnitz DJ, Butler JE. Quantitation of bovine β2-microglobulin: occurrence in body fluids, on milk fat globules and origin in milk. Mol Immunol. 1985;22(7):779–86.
Fedosov SN, Petersen TE, Nexø E. Transcobalamin from cow milk: isolation and physico-chemical properties. Biochim Biophys Acta (BBA)-Protein Struct Mol Enzymol. 1996;1292(1):113–9.
Hine B, Boggs I, Green R, Miller JW, Hovey RC, Humphrey R, Wheeler TT. Transcobalamin derived from bovine milk stimulates apical uptake of vitamin B12 into human intestinal epithelial cells. J Cell Biochem. 2014;115(11):1948–54.
Hvarregaard J, Andersen MH, Berglund L, Rasmussen JT, Petersen TE. Characterization of glycoprotein PAS-6/7 from membranes of bovine milk fat globules. FEBS J. 1996;240(3):628–36.
Yolken RH, Peterson JA, Vonderfecht SL, Fouts ET, Midthun K, Newburg DS. Human milk mucin inhibits rotavirus replication and prevents experimental gastroenteritis. J Clin Invest. 1992;90(5):1984–91.
Becker-Dreps S, Choi WS, Stamper L, Vilchez S, Velasquez DE, Moon SS, Permar SR. Innate immune factors in Mothers’ breast Milk and their lack of association with rotavirus vaccine immunogenicity in Nicaraguan infants. J Pediatr Infect Dis Soc. 2017;6(1):87–90.
Lundén A, Gustafsson V, Imhof M, Gauch R, Bosset JOJ. Dairy Res. 2002;69(3):383–90.
Peter H, Maria K. Studies on the mechanism of hepatic microsomal N-oxide formation. The role of cytochrome P-450 and mixed-function amine oxidase in the N-oxidation of NN-dimethylaniline. Biochem J. 1977;164(3):487.
Spellacy E, Watts RWE, Goolamali SK. Trimethylaminuria. J Inherit Metab Dis. 1979;2:85–8.
Eswaramoorthy S, Bonanno JB, Burley SK, Swaminathan S. Mechanism of action of a flavin-containing monooxygenase. Proc Natl Acad Sci. 2006;103(26):9832–7.
Cashman JR, Zhang J. Human flavin-containing monooxygenases. Annu Rev Pharmacol Toxicol. 2006;46:65–100.
Krueger SK, Williams DE. Mammalian flavin-containing monooxygenases: structure/function, genetic polymorphisms and role in drug metabolism. Pharmacol Ther. 2005;106(3):357–87.
Chen GP, Ziegler DM. Liver microsome and flavin-containing monooxygenase catalyzed oxidation of organic selenium compounds. Arch Biochem Biophys. 1994;312(2):566–72.
Bhat SA, Ahmad SM, Ibeagha-Awemu EM, Bhat BA, Dar MA, Mumtaz PT, Ganai NA. Comparative transcriptome analysis of mammary epithelial cells at different stages of lactation reveals wide differences in gene expression and pathways regulating milk synthesis between Jersey and Kashmiri cattle. PLoS One. 2019;14(2).
Feder ME, Walser JC. The biological limitations of transcriptomics in elucidating stress and stress responses. J Evol Biol. 2005;18(4):901–10.
Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2015;44(D1):D457–62.
Incassati A, Chandramouli A, Eelkema R, Cowin P. Key signaling nodes in mammary gland development and cancer: β-catenin. Breast Cancer Res. 2010;12(6):213.
Yu QC, Verheyen EM, Zeng YA. Mammary development and breast cancer: a Wnt perspective. Cancers. 2016;8(7):65.
Brisken C, Heineman A, Chavarria T, Elenbaas B, Tan J, Dey SK, et al. Essential function of Wnt-4 in mammary gland development downstream of progesterone signaling. Genes Dev. 2000;14(6):650–4.
Zhang X, Martinez D, Koledova Z, Qiao G, Streuli CH, Lu P. FGF ligands of the postnatal mammary stroma regulate distinct aspects of epithelial morphogenesis. Development. 2014;141(17):3352–62.
Turner N, Grose R. Fibroblast growth factor signalling: from development to cancer. Nat Rev Cancer. 2010;10(2):116.
Furth PA, Nakles RE, Millman S, Diaz-Cruz ES, Cabrera MC. Signal transducer and activator of transcription 5 as a key signaling pathway in normal mammary gland developmental biology and breast cancer. Breast Cancer Res. 2011;13(5):220.
Oikonomopoulou K, Ricklin D, Ward PA, Lambris JD. Interactions between coagulation and complement—their role in inflammation. Semin Immunopathol. 2012;341:51–165.
Carroll MC. Complement and humoral immunity. Vaccine. 2008;26:I28–33.
Kunz C, Lönnerdal B. Human-milk proteins: analysis of casein and casein subunits by anion-exchange chromatography, gel electrophoresis, and specific staining methods. Am J Clin Nutr. 1990;51(1):37–46.
Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, Inuganti A, Griss J, Mayer G, Eisenacher M, Pérez E, Uszkoreit J, Pfeuffer J, Sachsenberg T, Yilmaz S, Tiwary S, Cox J, Audain E, Walzer M, Jarnuczak AF, Ternent T, Brazma A, Vizcaíno JA. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47(D1):D442–50 PubMed ID: 30395289.
The authors would like to thank the staff of the University Dairy Farm of Mountain Livestock Research Institute-Manasbal, Srinagar-Kashmir for care and management of the experimental animals. We also acknowledge the DBT-supported bioinformatics Infrastructuralfacility for data analysis.
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE  partner repository with the dataset identifier PXD017412
This work was supported by the Department of Biotechnology, Ministry of Science and Technology, Government of India under Biotechnology Research and Development scheme. The funding body has no role in the study design and data collection, analysis, interpretation of data and in writing the manuscript.
The ethical clearance was approved by the Institutional Animal Ethics Committee (IAEC) of Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Significantly upregulated proteins in Kashmiri cattle (high and low abundant proteins) Table S2. Significantly upregulated proteins in Jersey cattle (high and low abundant proteins). Table S3. Sample chromatograms of individual samples. Table S4-S6. Peptide information of significantly upregulated proteins in Jersey cattle. Table S7-S9. Peptide information of significantly upregulated proteins in Kashmiri cattle.
Classification of significantly upregulated proteins in Kashmiri and Jersey cattle.
Total proteome obtained from different samples.
About this article
Cite this article
Bhat, S.A., Ahmad, S.M., Ibeagha-Awemu, E.M. et al. Comparative milk proteome analysis of Kashmiri and Jersey cattle identifies differential expression of key proteins involved in immune system regulation and milk quality. BMC Genomics 21, 161 (2020). https://doi.org/10.1186/s12864-020-6574-4
- Milk proteome
- FMO3 enzyme