Advertisement
Research| Volume 102, ISSUE 5, P3839-3850, May 2019

Proteomics and microstructure profiling of goat milk protein after homogenization

Open ArchivePublished:March 01, 2019DOI:https://doi.org/10.3168/jds.2018-15363

      ABSTRACT

      This study investigated the protein changes in goat milk during the homogenization process using label-free quantification. We quantified 310 and 315 proteins in the control group (CG) and homogenized group (HG), respectively, and 16 proteins were significantly different between the 2 groups. For HG, the goat milk protein particle sizes were smaller and more evenly distributed and exhibited an increase in the regular arrangement of the secondary structures. Proteomics analysis verified that xanthine dehydrogenase and asparaginase-like 1 expression in CG were higher than in HG, whereas the opposite was observed for fructose-bisphosphate aldolase, κ-casein, and β-casein. Significant changes were found in the homogenization-treated goat milk proteome that were related to goat milk glycolysis/gluconeogenesis metabolism. This work provides updated information on the current proteome characteristics of homogenized goat milk, which may be important for applying the protein component of goat milk to human nutrition and health.

      Key words

      INTRODUCTION

      Homogenization has become a standardized processing method in the dairy industry. In homogenization, mechanical forces reduce fat globule size, alter the structure of milk proteins, and reduce enzyme activity and vitamin levels (
      • Pereda J.
      • Jaramillo D.P.
      • Quevedo J.M.
      • Ferragut V.
      • Guamis B.
      • Trujillo A.J.
      Characterization of volatile compounds in ultra-high-pressure homogenized milk.
      ). Temperature elevation also occurs due to a variety of phenomena, including shear stress, turbulence, cavitation, and the transformation of kinetic energy into heat during the pressure increase (
      • Amador-Espejo G.G.
      • Suàrez-Berencia A.
      • Juan B.
      • Bárcenas M.E.
      • Trujillo A.J.
      Effect of moderate inlet temperatures in ultra-high-pressure homogenization treatments on physicochemical and sensory characteristics of milk.
      ). These modifications generate a synergistic effect with the mechanical changes happening during the homogenization process. Although homogenization reduces Trp fluorescence intensity, it does not significantly change the tertiary or secondary structure of casein (
      • Qi P.X.
      • Ren D.
      • Xiao Y.
      • Tomasula P.M.
      Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
      ).
      In recent years, the relatively new and rapidly developing field of proteomics has offered a key technique to identify, characterize, and quantify proteins. Goat milk proteins are composed of caseins (α-, β-, and κ-CN), whey proteins (WP; α-LA, β-LG, BSA, immunoglobulins, and other minor proteins), proteose peptones, and membrane proteins [mostly milk fat globule membrane protein (MFGMP);
      • Chen D.
      • Zhao X.
      • Li X.
      • Wang J.
      • Wang C.
      Milk compositional changes of Laoshan goat milk from partum up to 261 days postpartum.
      ;
      • Li Y.
      • Chen D.
      • Li J.
      • Zhang X.-X.
      • Wang C.-F.
      • Wang J.-M.
      Changes in superoxide dismutase activity postpartum from Laoshan goat milk and factors influencing its stability during processing.
      ]. Furthermore, 9 αS1-CN, 4 αS2-CN, 12 β-CN, 14 κ-CN, 3 α-LA, and 11 β-LG have been identified in bovine milk protein (
      • Le T.T.
      • Deeth H.C.
      • Larsen L.B.
      Proteomics of major bovine milk proteins: Novel insights.
      ). A variety of milk proteins have been verified by applying quantitative proteomics and bioinformatics to milk proteins using label-free quantification (
      • Anagnostopoulos A.K.
      • Katsafadou A.I.
      • Pierros V.
      • Kontopodis E.
      • Fthenakis G.C.
      • Arsenos G.
      • Karkabounas S.C.
      • Tzora A.
      • Skoufos I.
      • Tsangaris G.T.
      Milk of Greek sheep and goat breeds; characterization by means of proteomics.
      ;
      • Lu J.
      • Liu L.
      • Pang X.
      • Zhang S.
      • Jia Z.
      • Ma C.
      • Zhao L.
      • Lv J.
      Comparative proteomics of milk fat globule membrane in goat colostrum and mature milk.
      ).
      During conventional homogenization, finer emulsions are achieved from smaller fat globules, improving the digestibility and absorption of milk fat (
      • Brick T.
      • Ege M.
      • Boeren S.
      • Böck A.
      • Von E.M.
      • Vervoort J.
      • Hettinga K.
      Effect of processing intensity on immunologically active bovine milk serum proteins.
      ). The natural fat globules (diameter 1–10 μm) are broken into small globules (diameter <1 μm), resulting in the destruction of MFGMP (
      • Qi P.X.
      • Ren D.
      • Xiao Y.
      • Tomasula P.M.
      Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
      ) and inadequate coverage of the newly formed fat surface, even if a portion of native membrane remains associated with the fat droplets (
      • Huppertz T.
      Homogenization of milk: High-pressure homogenizers.
      ). In addition, homogenization to reduce the size of the fat globules also substantially changes casein, WP, and MFGMP by altering protein–protein interactions. To cover the increased fat globule surface area and prevent the fat from coalescing in homogenized milk, surface-active proteins, especially casein, are adsorbed onto newly formed globules. Treated milk is whiter than raw milk due to the increase in the number of fat globules, which diffract light more efficiently. Caseins are preferentially adsorbed over WP at the oil–water interface in homogenized milk (
      • Cano-Ruiz M.E.
      • Richter R.L.
      Effect of homogenization pressure on the milk fat globule membrane proteins.
      ), which could be useful for fine-tuning processing conditions when manufacturing dairy products (
      • Liang Y.
      • Matia-Merino L.
      • Gillies G.
      • Patel H.
      • Ye A.
      • Golding M.
      The heat stability of milk protein-stabilized oil-in-water emulsions: A review.
      ).
      Previous studies focused on the effects of temperature, pressure, time, and number of treatments on fats, proteins, and the activity of enzymes and vitamins in dairy products (
      • Pinho C.R.G.
      • Franchi M.A.
      • Tribst A.A.L.
      • Cristianini M.
      Effect of ultra high pressure homogenization on alkaline phosphatase and lactoperoxidase activity in raw skim milk.
      ). Changes in protein functionality as well as protein and fat absorption after homogenization have also been assessed (
      • Ye A.
      • Cui J.
      • Dalgleish D.
      • Singh H.
      Effect of homogenization and heat treatment on the behavior of protein and fat globules during gastric digestion of milk.
      ); however, no one has studied the proteome of goat milk after homogenization. In the present study, the changes in goat milk protein composition and expression during homogenization processing were identified via label-free technology based on the MaxQuant algorithm for fast and cost-efficient protein quantification. Gene Ontology (GO; http://www.geneontology.org/) annotation and the Kyoto Encyclopedia of Genes and Genomes (KEGG; https://www.genome.jp/kegg/) were used to analyze and classify proteomics data, and morphological analysis of goat milk protein was used to predict microstructural changes in the protein in homogenized goat milk. The effect of homogenization on the glycolytic pathway was then evaluated by comparing changes in homogenized milk protein expression. This work expands our knowledge of the goat milk proteome and is helpful to further clarify the differences between homogenized and fresh goat milk.

      MATERIALS AND METHODS

      Sample Collection and Preparation

      Thirty Laoshan goats aged 30 to 32 mo with no signs of acute mastitis or other clinical diseases were selected from Three Hi Goat Farms in Shandong Province (Zhujiayu Village, Daolang Town, Tai'an City, China). This farm is well organized and uses freestall housing and a TMR silage diet. Feces are cleared manually and milking is performed manually at 1000 h. All goat milk samples were collected at 75 d after lactation and transported to the laboratory with dry ice within 1 h. Fresh mature milk samples from these goats were randomly divided into 3 groups, with 10 in each group, and pooled. These 3 mixed mature milk samples were filtered and divided into 2 groups again for a total of 6 samples. The samples in the first group served as the control group (CG) and the other group was the homogenized group (HG), which was treated at 65 ± 2°C and 20 ± 2 MPa (Figure 1). Whole-milk samples were supplemented with SDT lysis buffer (4% SDS, 100 mmol/L Tris-HCl, 1 mmol/L dithiothreitol, pH 7.6) and boiled for 15 min. After centrifugation at 14,000 × g for 40 min at 4°C (Eppendorf 5430R, Hamburg, Germany), the content of protein in the aqueous phase was determined using the BCA Protein Assay Kit (Bio-Rad, Hercules, CA). The supernatant was used for filter-aided sample preparation as described below (
      • Wisniewski J.R.
      • Zougman A.
      • Nagaraj N.
      • Mann M.
      Universal sample preparation method for proteome analysis.
      ).
      Figure thumbnail gr1
      Figure 1Workflow for collection and preparation of control and homogenized goat milk samples.

      Filter-Aided Sample Preparation

      Protein samples (200 μg) were diluted in SDT buffer, incubated for 5 min in boiling water, and cooled to room temperature. Two centrifugations at 14,000 × g for 15 min at 4°C were performed after adding 2 volumes of 200 μL of UA buffer (8 mol/L urea, 150 mmol/L Tris-HCl, pH 8.0). We then added 100 μL of 100 mmol/L iodoacetamide (163–2109, Bio-Rad) in UA buffer and the samples were incubated in the dark for 30 min. The filters were washed 3 times with 100 μL of UA buffer and twice with 100 μL of 25 mmol/L NH4HCO3 buffer. Trypsin buffer (4 μg of trypsin in 40 μL of 100 mmol/L NH4HCO3) was then added to the samples overnight at 37°C. The peptides from each sample were desalted on C18 cartridges [Empore C18 SPE cartridges (standard density); bed i.d., 7 mm; volume, 3 mL; Sigma-Aldrich, St. Louis, MO] and reconstituted in 40 µL of 0.1% (vol/vol) formic acid (
      • Wisniewski J.R.
      • Zougman A.
      • Nagaraj N.
      • Mann M.
      Universal sample preparation method for proteome analysis.
      ). These samples were prepared for analysis using LC/MS.

      LC/MS

      The peptide mixture was injected onto a reverse-phase trap column (Acclaim PepMap100, 100 µm × 2 cm, nanoViper C18, Thermo Fisher Scientific, Waltham, MA) followed by a C18 reverse-phase analytical column (Easy Column, 10 cm long, 75 μm inner diameter, 3 μm resin, Thermo Fisher Scientific). The mobile LC phase contained buffer A (0.1% formic acid) and B (84% acetonitrile and 0.1% formic acid) and was eluted at a flow rate of 300 nL/min with linear gradients. The MS was programmed using a data-dependent top-10 method. The survey scan was from 300 to 1,800 m/z with a resolution of 70,000 at 200 m/z (
      • Zhu Y.
      • Xu H.
      • Chen H.
      • Xie J.
      • Shi M.
      • Shen B.
      • Deng X.
      • Liu C.
      • Zhan X.
      • Peng C.
      Proteomic analysis of solid pseudopapillary tumor of the pancreas reveals dysfunction of the endoplasmic reticulum protein processing pathway.
      ).

      ELISA

      Two candidate proteins (κ-CN and β-LG) were chosen to further verify the proteomics results because they are abundant and play an important role after homogenization. The content of κ-CN (http://www.jonln.com/) and β-LG (http://www.cloud-clone.com/) in different samples was measured according to the manual provided with the ELISA kits. The ELISA kit for κ-CN (cat. no. JL46016) was obtained from Jianglai Biological (Shanghai, China); the ELISA kit for β-LG (cat. no. SEBO23 Bo) was obtained from Cloud-Clone Corp. (Wuhan, China).

      Data Analysis

      The MS data were analyzed using MaxQuant software version 1.3.0.5 (Max Planck Institute of Biochemistry, Martinsried, Germany) with the parameters in Table 1 (
      • Cox J.
      • Hein M.Y.
      • Luber C.A.
      • Paron I.
      • Nagaraj N.
      • Mann M.
      Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ.
      ). Contaminated protein was discarded during this process.
      Table 1Description of related parameters for data analysis
      ItemValue
      EnzymeTrypsin
      Maximum missed cleavages2
      Main search6 ppm
      First search20 ppm
      MS/MS tolerance20 ppm
      Fixed modificationsCarbamidomethyl (C)
      Variable modificationsOxidation (M), acetyl (protein N-term)
      Databaseuniprot_Caprinae_32203_20170613.fasta
      Database patternReverse
      Peptide false discovery rate (FDR)≤0.01
      Protein FDR≤0.01
      Time window (match between runs)2 min
      Protein quantificationRazor and unique peptides were used for protein quantification
      Label-free quantification (LFQ)True
      LFQ minimum ratio count1

      Cluster Analysis, GO, and KEGG Enrichment Analyses

      Protein expression data were used to perform hierarchical clustering analysis. Cluster 3.0 (http://bonsai.hgc.jp/∼mdehoon/software/cluster/software.htm) and Java Treeview software (http://jtreeview.sourceforge.net) were used. The goat reference database for differentially expressed protein searches was downloaded as FASTA files from the UniProtKB database (http://www.uniprot.org/, June 13, 2017, release;
      • Götz S.
      • Garcia-Gomez J.M.
      • Terol J.
      • Williams T.D.
      • Nagaraj S.H.
      • Nueda M.J.
      • Robles M.
      • Talon M.
      • Dopazo J.
      • Conesa A.
      High-throughput functional annotation and data mining with the Blast2GO suite.
      ). The FASTA protein sequences of differentially expressed proteins were compared with online GO database to retrieve GO and were subsequently mapped to pathways in KEGG UniProt-Caprinae database (
      • Moriya Y.
      • Itoh M.
      • Okuda S.
      • Yoshizawa A.
      • Kanehisa M.
      KAAS: An automatic genome annotation and pathway reconstruction server.
      ).

      Field Emission Scanning Electron Microscope

      A 10% acetic acid solution was used to adjust the pH (4.1) to the casein isoelectric point for casein and WP collection (
      • Zhu Y.
      • Wang J.
      • Wang C.
      Research on the preparation, uniformity and stability of mixed standard substance for rapid detection of goat milk composition.
      ). The WP and casein solution at 100 mg/mL was prepared in PBS (0.1 M, pH 7.0) containing Na2HPO4 (0.577 mol/L) and NaH2PO4 (0.423 mol/L). The protein solution (5 μL) was applied dropwise onto clean chips and air-dried at room temperature. The chips with the attached protein were covered with 2.5% glutaraldehyde and kept overnight at 4°C, after which they were washed thrice with PBS. The samples were passed once through a graded ethanol series (10, 30, 50, 70, 90, 95, and 100%) individually for 15 min each, transferred to a dryer, and incubated for 10 min at 37°C to prepare them for field emission scanning electron microscope analysis. The microstructure of WP and casein was analyzed with a SU8010 (Hitachi, Tokyo, Japan; http://www.hitachi.com/) field emission scanning electron microscope made in Japan.

      Fourier-Transform Infrared Measurements and Analysis

      The Fourier transform infrared (FTIR) spectra of the protein samples were recorded with an FTIR spectrometer (IRPrestige-21, Shimadzu, Kyoto, Japan) in the wavelength range 4,000 to 400 cm−1 with a resolution of 4 cm−1 using 10 scans. Spectral collection, peak identification, and basic infrared spectrogram analysis were performed using Omnic 8.2 (Nicolet, Thermo Fisher Scientific). The second derivative was calculated in the amide I region, 1,700 to 1,600 cm−1. Ten major peaks were resolved and fitted using Peakfit version 4.12 (Seasolve, San Jose, CA). A Gaussian function was used during the fitting routine for all peaks corresponding to α-helix (1,650–1,660 cm−1), β-sheet (1,600–1,640 cm−1 and 1,682–1,700 cm−1), β-turn (1,661–1,681 cm−1), and random coil (1,641–1,650 cm−1;
      • Long G.
      • Ji Y.
      • Pan H.
      • Sun Z.
      • Li Y.
      • Qin G.
      Characterization of thermal denaturation structure and morphology of soy glycinin by FTIR and SEM.
      ). The areas of all component bands assigned to a given secondary structural element were summed and divided by the total integrated area to yield the percentage of the element.

      RESULTS

      Quantification and Hierarchical Clustering of Proteins

      We investigated the protein changes in goat milk after homogenization using label-free quantification. Venn diagrams showed 224 and 259 proteins common to the 3 biological replicates in the CG and HG samples as well as proteins uniquely expressed in goat milk in each group (Figure 2A and B). As seen in Figure 2C, 310 and 315 proteins were quantified in the CG and HG, respectively, and 307 proteins were common between the 2 groups. Sixteen different proteins were also found between CG and HG by a t-test of the label-free quantification intensity data (P < 0.05; Figure 2D). Cluster analysis of the different groups revealed that CG1, CG2, and CG3 clustered and that HG1, HG2, and HG3 were clustered, as depicted in Figure 2D. Furthermore, the cluster graph could be divided into clusters I and II (Figure 2D).
      Figure thumbnail gr2
      Figure 2Venn diagram of quantified proteins in 3 biological duplicates of control group (CG; A) and homogenized group (HG; B) and the number of protein groups commonly present in goat milk from the CG and HG (C). Hierarchical clustering of differentially expressed proteins between the CG and HG (D).

      GO Analysis of Identified Proteins

      Gene Ontology is a standardized system of gene function classification in which gene properties and products are described using the categories biological process (BP), molecular function (MF), and cellular component (CC). We determined the biological activity of proteins with significant differences (P < 0.05) based on MF, localized proteins based on CC, and participating biological pathways based on BP (Figure 3). The main BP were involved in cellular process (12 proteins) and single-organism process (11 proteins). Biological regulation (10 proteins) and regulation of biological process (10 proteins) also contributed to the BP designation. In addition, multicellular organismal process (9 proteins), localization (9 proteins), and response to stimulus (9 proteins) were well represented (Figure 3A).
      Figure thumbnail gr3
      Figure 3Gene Ontology (GO; http://www.geneontology.org/) analysis for biological processes (BP; A), molecular functions (MF; B), and cellular components (CC; C) of identified proteins.
      The most obviously enriched MF for the identified proteins were related to binding (12 proteins) followed by catalytic activity (8 proteins; Figure 3B). The identified proteins were mainly described as from the organelle (13 proteins), cell (13 proteins), cell part (13 proteins), extracellular region part (10 proteins), and extracellular region (10 proteins). In addition, the membrane (9 proteins) was present in goat milk proteins (Figure 3C).

      Analysis of KEGG Pathways

      Homogenization affects protein structure and protein expression level (Figure 2D), which alters the metabolic pathways corresponding to the significantly different proteins (P < 0.05). Many metabolic pathways were tightly linked to glycolysis, as illustrated in Figure 4. These pathways included carbon fixation in photosynthetic organisms, fructose and mannose metabolism, carbon metabolism, biosynthesis of AA, purine metabolism, and caffeine metabolism.
      Figure thumbnail gr4
      Figure 4Significantly different proteins involved in glycolysis and gluconeogenesis (P < 0.05). PRPP = 5-phosphoribosyl 1-pyrophosphate. Numbers represent the EC numbers of the respective enzymes.

      κ-CN and β-LG Proteins in Goat Milk from the CG and HG

      To further validate the LC/MS results, the content of κ-CN and β-LG in the CG and HG were determined by ELISA. As Figure 5A shows, the standard curve coverage of κ-CN ranged from 6.25 to 200 μg/mL. Figure 5B shows that the standard curve of β-LG were 0 to 5 ng/mL. As expected, the variations of β-LG and κ-CN content in CG and HG as determined by ELISA were consistent with the LC/MS results (Figure 5C and D). In addition, Figure 5C and D revealed that κ-CN and β-LG levels in the HG were slightly higher than those in the CG, and these ELISA results also agreed with the LC/MS data. The consistent variation trends determined from the ELISA and LC/MS results indicate that the label-free technology used in this study was robust.
      Figure thumbnail gr5
      Figure 5Standard curves for κ-CN (A) and β-LG (B) by ELISA and comparison of LC/MS and ELISA results for κ-CN (C) and β-LG (D). CG = control group; HG = homogenized group. Error bars represent SD.

      Morphological Analysis of Goat Milk Protein

      The morphology of WP and casein from the CG and HG observed by field emission scanning electron microscopy is shown in Figure 6. The WP in CG (Figure 6A) and HG (Figure 6B) exhibits a regular spherical structure, but the WP in CG shows slight agglomeration. The distribution of WP particles in the homogenized goat milk was more uniform, the particle size was finer, the particle surface was smoother, and the arrangement was more regular, displaying no agglomeration. As depicted in Figure 6C and D, in fresh goat milk, casein particles were attached to each other to form micelles, and these micelles were connected to each other again to form a compact network structure. In the homogenized goat milk, the size distribution of the casein micelles was more uniform and the micelles were more spherical, the agglomeration phenomenon was substantially reduced, micelle size was reduced, and the surface morphology was more regular.
      Figure thumbnail gr6
      Figure 6Field emission scanning electron microscope images of whey protein and casein from the control group and homogenized group. A = whey protein from the control group, B = whey protein from the homogenized group, C = casein micelles from the control group, and D = casein micelles from the homogenized group.

      FTIR Analysis of WP and Casein

      The apparent secondary structure of WP and casein from the CG and HG was analyzed using FTIR spectroscopic methods. Figure 7A, D, B, and E show FTIR spectra and the second derivative of the WP amide band from the CG and HG. The FTIR spectra in the amide I region (1,700–1,600 cm−1) is known for its sensitivity to protein secondary conformational change. The characteristic absorption peaks for raw milk and homogenized milk were 1,654.9 and 1,647.2 cm−1, respectively. As illustrated in Figure 7C and F, the amide I region of the deconvoluted FTIR spectrum of WP from the CG was compared with that of WP from the HG. Figure 7G and H compare the secondary structural content of WP and casein from the CG with that of WP and casein from the HG. It was clear that there were significant differences in α-helix, β-sheet, and random coil for WP between the CG and HG (P < 0.05). However, for casein, there was no significant difference between the 2 groups (P < 0.05).
      Figure thumbnail gr7
      Figure 7Fourier transform infrared spectra of the control group (CG; A) and homogenized group (HG; D); the second derivatives of amide I of the CG (B) and HG (E); Gaussian curve fitting of amide I of the CG (C) and HG (F); and comparison of secondary structural content of WP (G) and CN (H) between the CG and HG. Different letters (a–d) indicate significant differences (P < 0.05) between the CG and HG. Error bars represent SD.

      DISCUSSION

      Effect of Homogenization on Goat Milk Proteins

      We quantified 315 proteins in goat milk treated by homogenization (Figure 2C). Previously observed major proteins, such as β-LG (B2YKY6), κ-CN (A0A0U5A2G6), β-CN (Q5YD57), xanthine dehydrogenase (XDH; W5PMT0), and aldolase (ALDO; W5NQP9), were all identified at high abundance in our study, which indicates the robustness of the methodology (Figure 2D). In previous investigations, 423 proteins were identified and XDH was significantly different in abundance between the colostrum and mature goat milk (
      • Lu J.
      • Wang X.
      • Zhang W.
      • Liu L.
      • Pang X.
      • Zhang S.
      • Lv J.
      Comparative proteomics of milk fat globule membrane in different species reveals variations in lactation and nutrition.
      ). In addition, 53 unique carbonylated peptides (37 carbonylation sites, 15 proteins) were identified in raw milk and differently processed milk products (
      • Milkovska-Stamenova S.
      • Mnatsakanyan R.
      • Hoffmann R.
      Protein carbonylation sites in bovine raw milk and processed milk products.
      ). We found 16 proteins to be significantly different between fresh milk and homogenized milk (P < 0.05). Figure 2D clearly shows that protein expression in cluster I of the CG was higher than that of the HG; in contrast, the protein expression in cluster II showed the opposite trend. The clustering results demonstrated protein degradation and denaturation in goat milk during homogenization and the formation of new products. In the homogenized groups, XDH showed lower expression and ALDO was abundant protein (Figure 2D). This tendency indicated that homogenization destroyed the original structure of fat globules. Analysis of these proteins would increase our understanding and knowledge of goat milk proteins after homogenization treatment.
      Several studies have focused on the effects of homogenization on fat globule microstructure in cow milk used confocal laser scanning microscopy (
      • Yao Y.
      • Zhao G.
      • Zou X.
      • Huang L.
      • Wang X.
      Microstructural and lipid composition changes in milk fat globules during milk powder manufacture.
      ).
      • Cadesky L.
      • Walkling-Ribeiro M.
      • Kriner K.T.
      • Karwe M.V.
      • Moraru C.I.
      Structural changes induced by high-pressure processing in micellar casein and milk protein concentrates.
      observed structural changes induced by high-pressure processing in casein micellar and milk protein concentrates. Previous research has shown that homogenization reduces milk fat globule size and decreases globule size range, resulting in a newly formed membrane (
      • Qi P.X.
      • Ren D.
      • Xiao Y.
      • Tomasula P.M.
      Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
      ). We focused on changes in WP and casein micelles (Figure 6), and morphology results showed that homogenization reduced WP size and increased protein distribution uniformity, with no obvious agglomeration phenomena observed (Figure 6A and B). Figure 7G shows that the α-helix content in the HG was much higher than that in the CG, although the percentage of random coil in the HG was significantly reduced compared with that in the CG (P < 0.05), which further indicates that homogenization stabilized WP. By changing milk protein granule size with homogenization treatment, the physicochemical properties of proteins, such as color, viscosity, pH, and acidity, were also affected (
      • Amador-Espejo G.G.
      • Suàrez-Berencia A.
      • Juan B.
      • Bárcenas M.E.
      • Trujillo A.J.
      Effect of moderate inlet temperatures in ultra-high-pressure homogenization treatments on physicochemical and sensory characteristics of milk.
      ). Ultimately, homogenization reduced the heat coagulation time and the thermal stability of milk (
      • Meena G.S.
      • Singh A.K.
      • Gupta V.K.
      • Borad S.
      • Arora S.
      • Tomar S.K.
      Effect of pH adjustment, homogenization and diafiltration on physicochemical, reconstitution, functional and rheological properties of medium protein milk protein concentrates (MPC70).
      ). The effects of heat coagulation time on microstructure and rheological properties, including solubility, heat stability, water binding, dispersibility, foaming, and emulsifying of hazelnut milks, have been investigated (
      • Gul O.
      • Saricaoglu F.T.
      • Mortas M.
      • Atalar I.
      • Yazici F.
      Effect of high pressure homogenization (HPH) on microstructure and rheological properties of hazelnut milk.
      ).
      Homogenization disrupts fat droplets and causes the globule membrane to rearrange, resulting in increased interactions between WP, casein, and fat globules (
      • Qi P.X.
      • Ren D.
      • Xiao Y.
      • Tomasula P.M.
      Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
      ). Figure 2D shows significant differences in the expression of XDH (W5PMT0) and asparaginase like 1 (ASRGL1; W5PZJ8) between the CG and HG (P < 0.05). The XDH was the most abundant protein in goat milk fat globule membrane (
      • Lu J.
      • Liu L.
      • Pang X.
      • Zhang S.
      • Jia Z.
      • Ma C.
      • Zhao L.
      • Lv J.
      Comparative proteomics of milk fat globule membrane in goat colostrum and mature milk.
      ). In addition, XDH (part of the MFGMP) was partially dissociated on the milk fat globule membrane surface (
      • Fauquant C.
      • Briard V.
      • Leconte N.
      • Michalski M.C.
      Differently sized native milk fat globules separated by microfiltration: Fatty acid composition of the milk fat globule membrane and triglyceride core.
      ). One possible reason for the decrease in XDH and ASRGL1 is the destruction of fat globulin protein during homogenization. When the milk samples were homogenized, XDH was removed from the fresh goat milk. In addition, the shear force and heat might destroy enzyme activity, resulting in decreased ASRGL1 levels.
      • Meena G.S.
      • Singh A.K.
      • Gupta V.K.
      • Borad S.
      • Arora S.
      • Tomar S.K.
      Effect of pH adjustment, homogenization and diafiltration on physicochemical, reconstitution, functional and rheological properties of medium protein milk protein concentrates (MPC70).
      suggested that homogenization disrupted casein aggregates due to shearing and dispersed them uniformly. The treatment enhanced the net negative charge and electrostatic repulsive forces between the casein micelles, preventing casein aggregation. As seen in Figure 6C and D, the original casein micelle structure was destroyed. Additionally, particle size distribution was relatively uniform, the agglomeration phenomenon was substantially reduced, and aggregation was decreased (Figure 6C and D). We noted a possible increase in κ-CN dissociation, which was validated by data showing that κ-CN levels in the HG were higher than those in the CG by ELISA and LC/MS (Figure 5). As expected, the content of κ-CN and β-CN in cluster I of the HG were also higher than those in the CG (Figure 2D). Due to homogenization, original internal β-CN was liberated from the destroyed casein micelles and was therefore significantly higher in the HG than in the CG (P < 0.05; Figure 2D). This result validated earlier reports indicating that the αS1-CN percentage does not change, although the κ-CN percentage increased and the dissociation of β-LG slightly increased (
      • Cano-Ruiz M.E.
      • Richter R.L.
      Effect of homogenization pressure on the milk fat globule membrane proteins.
      ). This was due to κ-CN dissociation from the fat globule surface, regardless of changes to the casein micelle microstructure. However, Figure 7H indicated that homogenization did not significantly change casein secondary structure (P < 0.05). For CG, the fat globules were simply trapped within the gel structure formed by the casein micelles, whereas for HG the fat globules were coated with caseins and actively are incorporated into the gel matrix (
      • Ye A.
      • Cui J.
      • Dalgleish D.
      • Singh H.
      Effect of homogenization and heat treatment on the behavior of protein and fat globules during gastric digestion of milk.
      ). The globules themselves help form the matrix as the clot forms from smaller fat globules in the homogenized milk (
      • Cano-Ruiz M.E.
      • Richter R.L.
      Effect of homogenization pressure on the milk fat globule membrane proteins.
      ). Thus, homogenized goat milk is more conducive to clot formation, and this change has a large influence on yogurt and cheese processing.
      Figure 5D shows that the β-LG content in the HG was higher than that in the CG. In ultra-high-pressure homogenization-treated milk, small amounts of β-LG might denature, which is likely due to the increase in mechanical forces and temperature (
      • Qi P.X.
      • Ren D.
      • Xiao Y.
      • Tomasula P.M.
      Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
      ) resulting from the transformation between sulfhydryl groups (–SH) and disulfide bonds (–S–S–;
      • Hernández A.
      • Harte F.M.
      Manufacture of acid gels from skim milk using high-pressure homogenization.
      ). The secondary structure of WP changed and ultimately increased in stability (
      • Qi P.X.
      • Ren D.
      • Xiao Y.
      • Tomasula P.M.
      Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
      ). Sulfhydryl groups and disulfide bonds are important for protein structure, protein–protein interaction, and function. Unfortunately, the destructive power was much less than heat damage (
      • Pereda J.
      • Ferragut V.
      • Quevedo J.M.
      • Guamis B.
      • Trujillo A.J.
      Heat damage evaluation in ultra-high pressure homogenized milk.
      ), which is why we observed no significant change in β-LG solely due to homogenization, although its level increased slightly (P < 0.05).
      • Qi P.X.
      • Ren D.
      • Xiao Y.
      • Tomasula P.M.
      Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
      showed that conventional homogenization treatment does not appear to affect the levels of free sulfhydryl in soluble WP (e.g., β-LG, α-LA, IgG) in whole milk compared with those in their skimmed counterparts.

      Proteins Primarily Involved in GO

      Gene Ontology results suggested that cellular processes had the largest proportion in BP, and binding was the most important in MF (Figure 3). At the same time, organelle, cell, and cell part were dominant in CC (Figure 3). In a recent paper, MFGMP mainly originated from the plasma membrane (59 proteins), cytosol (32 proteins), and cytoplasmic vesicle (31 proteins;
      • Lu J.
      • Wang X.
      • Zhang W.
      • Liu L.
      • Pang X.
      • Zhang S.
      • Lv J.
      Comparative proteomics of milk fat globule membrane in different species reveals variations in lactation and nutrition.
      ). According to previous research, localization for the major classes of WP from sheep is extracellular, followed by membrane and cytoplasmic (
      • Anagnostopoulos A.K.
      • Katsafadou A.I.
      • Pierros V.
      • Kontopodis E.
      • Fthenakis G.C.
      • Arsenos G.
      • Karkabounas S.C.
      • Tzora A.
      • Skoufos I.
      • Tsangaris G.T.
      Milk of Greek sheep and goat breeds; characterization by means of proteomics.
      ).
      • Yang M.
      • Cao X.
      • Wu R.
      • Liu B.
      • Ye W.
      • Yue X.
      • Wu J.
      Comparative proteomic exploration of whey proteins in human and bovine colostrum and mature milk using iTRAQ-coupled LC-MS/MS.
      found that proteins combined with a partner molecule or other proteins to form protein complexes that expand protein function. For example, specific binding to antigens leads to complement activation and promote phagocytosis, WP binding to proteins and ions, and then association with hydrolase and enzyme inhibitor activity.
      Based on GO functional annotations in other research, quantified bovine milk WP are related to cellular processes, intracellular organelles, and molecular functions (
      • Yang M.
      • Cao X.
      • Wu R.
      • Liu B.
      • Ye W.
      • Yue X.
      • Wu J.
      Comparative proteomic exploration of whey proteins in human and bovine colostrum and mature milk using iTRAQ-coupled LC-MS/MS.
      ).
      • Lu J.
      • Liu L.
      • Pang X.
      • Zhang S.
      • Jia Z.
      • Ma C.
      • Zhao L.
      • Lv J.
      Comparative proteomics of milk fat globule membrane in goat colostrum and mature milk.
      noted that MFGMP were involved in complement and coagulation cascades, whereas adherens junction and galactose metabolism were only observed in colostrum. We focused on the effect of homogenization treatment on ALDO. In our study, ALDO in BP included cellular component organization or biogenesis, developmental process, cellular process, single-organism process, metabolic process, multicellular organismal process, regulation of biological process, biological adhesion, localization, and biological regulation. The MF for ALDO included binding and catalytic activity, as shown in Figure 3; ALDO was also involved in all CC in Figure 3 except for membrane part, synapse part, and synapse. On the basis of these findings, ALDO is involved in glycolysis and gluconeogenesis metabolism.

      KEGG Pathway Analysis

      • Ye A.
      • Cui J.
      • Dalgleish D.
      • Singh H.
      Effect of homogenization and heat treatment on the behavior of protein and fat globules during gastric digestion of milk.
      studied the effect of homogenization and heat treatment on the behavior of protein and fat globules during the gastric digestion of milk. After homogenization, milk fat is easily digested and absorbed because the diameter of the milk fat globule is reduced.
      • Michalski M.C.
      • Januel C.
      Does homogenization affect the human health properties of cow's milk?.
      suggested that changes in homogenized milk could improve milk digestibility and enhance protein bioavailability. The fat and lactose molecules of goat milk shrank during homogenization, and smaller lactose particles are more conducive to digestion and absorption (
      • Korpela R.
      • Paajanen L.
      • Tuure T.
      Homogenization of milk has no effect on the gastrointestinal symptoms of lactose intolerant subjects.
      ). In the human small intestine, lactose is hydrolyzed by lactase into glucose and galactose and then absorbed. Hexokinase, phosphofructokinase-1, and pyruvate kinase are key enzymes in glycolysis, especially phosphofructokinase-1 (
      • Yang J.
      • An D.
      • Zhang P.
      Expression profiling of cassava storage roots reveals an active process of glycolysis/gluconeogenesis.
      ). In our study, we noted a significant difference in ALDO expression after homogenization (P < 0.05), and ALDO also had an effect on glycolysis (Figure 4). This could be because homogenization causes changes in protein expression, which in turn affects glycolysis and gluconeogenesis.
      In general, ALDO showed significant changes in expression between the 2 groups (P < 0.05), with the HG showing higher levels (Figure 2D). Aldolase promotes conversion of glyceraldehyde-3P into sedoheptulose-1,7-bisphosphate through carbon fixation in photosynthetic organisms (Figure 4), which ultimately improves glycerate-3P production. During fructose and mannose metabolism, ALDO converts d-fructose to generate d-fructose-1P and glycerone-P (Figure 4); ALDO also converts fructose-6P to glycerone-P in carbon metabolism, which can be converted to glyceraldehyde-3P (Figure 4). Via glycolysis and gluconeogenesis, ALDO converts α/β-d-glucose generated from β-d-fructose-1,6P2 to glyceraldehyde-3P (Figure 4). Finally, ALDO converts glyceraldehyde-3P to glycerate-3P during glycolysis metabolism (Figure 4). Therefore, the presence of ALDO during homogenization promotes the conversion of glucose sugar to pyruvate. In the role of metabolism, pyruvate mainly (1) reduces to lactic acid, providing energy, and (2) goes into mitochondria to form acetyl CoA for the citrate cycle (
      • Nelson D.L.
      • Cox M.M.
      Lehninger Principles of Biochemistry. W. H.
      ).
      In addition, KEGG analysis indicated that xanthine-transformed glyoxylate also participates in glyoxylate and dicarboxylate metabolism. The XDH expression was significantly higher in the CG than the HG (P < 0.05; Figure 2D). From caffeine metabolism, we found that the HG accumulates more xanthine (Figure 4). In addition, XDH has low levels in the HG, so xanthine was more likely to be converted to glycine than urate in purine metabolism (Figure 4). This conversion reduces the pressure on glycolysis and gluconeogenesis. Based on the metabolic pathways altered by proteins that are significantly different after homogenization (P < 0.05), homogenization improves glycolysis and gluconeogenesis in goat milk; this may promote the production of pyruvic acid. Pyruvate can realize the mutual transformation of sugar, fat, and AA in vivo through acetyl-CoA and tricarboxylic acid cycle (
      • Nelson D.L.
      • Cox M.M.
      Lehninger Principles of Biochemistry. W. H.
      ). Moreover, changes in external conditions (such as pressure or temperature) produce different effects in goat milk (
      • Pinho C.R.G.
      • Franchi M.A.
      • Tribst A.A.L.
      • Cristianinia M.
      Effect of high pressure homogenization process on Bacillus stearothermophilus and Clostridium sporogenes spores in skim milk.
      ). Previous studies found that ultra-high pressure homogenization destroys microorganisms mechanically and with sudden pressure drops (
      • Dumay E.
      • Chevalier-Lucia D.
      • Picart-Palmade L.
      • Benzaria A.
      • Gràcia-Julià A.
      • Blayo C.
      Technological aspects and potential applications of (ultra) high-pressure homogenisation.
      ). Ultra-high-pressure homogenization was also shown to preserve vitamin C and riboflavin better than thermal milk treatment (
      • Sharabi S.
      • Okun Z.
      • Shpigelman A.
      Changes in the shelf life stability of riboflavin, vitamin C and antioxidant properties of milk after (ultra) high pressure homogenization: Direct and indirect effects.
      ). This technology has been used to improve rheological characteristics, physical properties, and milk stability as well as the polysaccharide and protein characteristics (
      • D'Incecco P.
      • Rosi V.
      • Cabassi G.
      • Hogenboom J.A.
      • Pellegrino L.
      Microfiltration and ultra-high-pressure homogenization for extending the shelf-storage stability of UHT milk.
      ). Homogenization and heat treatment can simultaneously improve the production process of dairy products (
      • Li Y.
      • Joyner Melito H.S.
      • Carter B.G.
      • Drake M.A.
      Effects of fat content, pasteurization method, homogenization pressure, and storage time on the mechanical and sensory properties of bovine milk.
      ).

      CONCLUSIONS

      By using microstructure information, we found that the goat milk protein became smaller and more evenly distributed, increasing the regular arrangement of the secondary structures, during the homogenization process. Proteomics analysis verified that the expression of XDH and ASRGL1 was higher in CG than in HG, whereas the opposite was observed for ALDO, κ-CN, and β-CN. Significant changes were found in the homogenized goat milk proteomics, which were related to glycolysis and gluconeogenesis. Based on this experimental study, future research will be dedicated to identifying the proteins involved in the deposits formed under homogenization of goat milk. It is important to study the effects of various homogeneous conditions of proteins through proteomics and then observe the changes in metabolism rather than just performing sensory evaluation in goat milk.

      ACKNOWLEDGMENTS

      This work was supported by the National Natural Science Foundation of China (31501501), Natural Science Foundation of Shandong Province (ZR2014JL019). Technical support by Shanghai Applied Protein Technology (Shanghai, China) is gratefully acknowledged.

      REFERENCES

        • Amador-Espejo G.G.
        • Suàrez-Berencia A.
        • Juan B.
        • Bárcenas M.E.
        • Trujillo A.J.
        Effect of moderate inlet temperatures in ultra-high-pressure homogenization treatments on physicochemical and sensory characteristics of milk.
        J. Dairy Sci. 2014; 97 (24342690): 659-671
        • Anagnostopoulos A.K.
        • Katsafadou A.I.
        • Pierros V.
        • Kontopodis E.
        • Fthenakis G.C.
        • Arsenos G.
        • Karkabounas S.C.
        • Tzora A.
        • Skoufos I.
        • Tsangaris G.T.
        Milk of Greek sheep and goat breeds; characterization by means of proteomics.
        J. Proteomics. 2016; 147 (27102495): 76-84
        • Brick T.
        • Ege M.
        • Boeren S.
        • Böck A.
        • Von E.M.
        • Vervoort J.
        • Hettinga K.
        Effect of processing intensity on immunologically active bovine milk serum proteins.
        Nutrients. 2017; 9 (28858242): 963
        • Cadesky L.
        • Walkling-Ribeiro M.
        • Kriner K.T.
        • Karwe M.V.
        • Moraru C.I.
        Structural changes induced by high-pressure processing in micellar casein and milk protein concentrates.
        J. Dairy Sci. 2017; 100 (28647329): 7055
        • Cano-Ruiz M.E.
        • Richter R.L.
        Effect of homogenization pressure on the milk fat globule membrane proteins.
        J. Dairy Sci. 1997; 80: 2732-2739
        • Cox J.
        • Hein M.Y.
        • Luber C.A.
        • Paron I.
        • Nagaraj N.
        • Mann M.
        Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ.
        Mol. Cell. Proteomics. 2014; 13 (24942700): 2513
        • D'Incecco P.
        • Rosi V.
        • Cabassi G.
        • Hogenboom J.A.
        • Pellegrino L.
        Microfiltration and ultra-high-pressure homogenization for extending the shelf-storage stability of UHT milk.
        Food Res. Int. 2018; 107 (29580509): 477-485
        • Chen D.
        • Zhao X.
        • Li X.
        • Wang J.
        • Wang C.
        Milk compositional changes of Laoshan goat milk from partum up to 261 days postpartum.
        Anim. Sci. J. 2018; 89 (29956422): 1355-1363
        • Dumay E.
        • Chevalier-Lucia D.
        • Picart-Palmade L.
        • Benzaria A.
        • Gràcia-Julià A.
        • Blayo C.
        Technological aspects and potential applications of (ultra) high-pressure homogenisation.
        Trends Food Sci. Technol. 2013; 31: 13-26
        • Fauquant C.
        • Briard V.
        • Leconte N.
        • Michalski M.C.
        Differently sized native milk fat globules separated by microfiltration: Fatty acid composition of the milk fat globule membrane and triglyceride core.
        Eur. J. Lipid Sci. Technol. 2010; 107: 80-86
        • Götz S.
        • Garcia-Gomez J.M.
        • Terol J.
        • Williams T.D.
        • Nagaraj S.H.
        • Nueda M.J.
        • Robles M.
        • Talon M.
        • Dopazo J.
        • Conesa A.
        High-throughput functional annotation and data mining with the Blast2GO suite.
        Nucl. Acids Res. 2008; 36: 3420-3435
        • Gul O.
        • Saricaoglu F.T.
        • Mortas M.
        • Atalar I.
        • Yazici F.
        Effect of high pressure homogenization (HPH) on microstructure and rheological properties of hazelnut milk.
        Innov. Food Sci. Emerg. Technol. 2017; 41: 411-420
        • Hernández A.
        • Harte F.M.
        Manufacture of acid gels from skim milk using high-pressure homogenization.
        J. Dairy Sci. 2008; 91 (18832197): 3761-3767
        • Huppertz T.
        Homogenization of milk: High-pressure homogenizers.
        in: Encyclopedia of Dairy Sciences. Vol. 10. NIZO Food Research, Ede, the Netherlands2011: 755-760
        • Korpela R.
        • Paajanen L.
        • Tuure T.
        Homogenization of milk has no effect on the gastrointestinal symptoms of lactose intolerant subjects.
        Milchwissenschaft. 2005; 60: 3-6
        • Le T.T.
        • Deeth H.C.
        • Larsen L.B.
        Proteomics of major bovine milk proteins: Novel insights.
        Int. Dairy J. 2017; 67: 2-15
        • Li Y.
        • Chen D.
        • Li J.
        • Zhang X.-X.
        • Wang C.-F.
        • Wang J.-M.
        Changes in superoxide dismutase activity postpartum from Laoshan goat milk and factors influencing its stability during processing.
        Ital. J. Anim. Sci. 2018; 17: 835-844
        • Li Y.
        • Joyner Melito H.S.
        • Carter B.G.
        • Drake M.A.
        Effects of fat content, pasteurization method, homogenization pressure, and storage time on the mechanical and sensory properties of bovine milk.
        J. Dairy Sci. 2018; 101 (29398020): 2941-2955
        • Liang Y.
        • Matia-Merino L.
        • Gillies G.
        • Patel H.
        • Ye A.
        • Golding M.
        The heat stability of milk protein-stabilized oil-in-water emulsions: A review.
        Curr. Opin. Colloid Interface Sci. 2017; 28: 63-73
        • Long G.
        • Ji Y.
        • Pan H.
        • Sun Z.
        • Li Y.
        • Qin G.
        Characterization of thermal denaturation structure and morphology of soy glycinin by FTIR and SEM.
        Int. J. Food Prop. 2015; 18: 763-774
        • Lu J.
        • Liu L.
        • Pang X.
        • Zhang S.
        • Jia Z.
        • Ma C.
        • Zhao L.
        • Lv J.
        Comparative proteomics of milk fat globule membrane in goat colostrum and mature milk.
        Food Chem. 2016; 209 (27173528): 10-16
        • Lu J.
        • Wang X.
        • Zhang W.
        • Liu L.
        • Pang X.
        • Zhang S.
        • Lv J.
        Comparative proteomics of milk fat globule membrane in different species reveals variations in lactation and nutrition.
        Food Chem. 2016; 196 (26593540): 665-672
        • Meena G.S.
        • Singh A.K.
        • Gupta V.K.
        • Borad S.
        • Arora S.
        • Tomar S.K.
        Effect of pH adjustment, homogenization and diafiltration on physicochemical, reconstitution, functional and rheological properties of medium protein milk protein concentrates (MPC70).
        J. Food Sci. Technol. 2018; 55 (29606752): 1376-1386
        • Michalski M.C.
        • Januel C.
        Does homogenization affect the human health properties of cow's milk?.
        Trends Food Sci. Technol. 2006; 17: 423-437
        • Milkovska-Stamenova S.
        • Mnatsakanyan R.
        • Hoffmann R.
        Protein carbonylation sites in bovine raw milk and processed milk products.
        Food Chem. 2017; 229 (28372194): 417
        • Moriya Y.
        • Itoh M.
        • Okuda S.
        • Yoshizawa A.
        • Kanehisa M.
        KAAS: An automatic genome annotation and pathway reconstruction server.
        Nucl. Acids Res. 2007; 35: W182-W185
        • Nelson D.L.
        • Cox M.M.
        Lehninger Principles of Biochemistry. W. H.
        Freeman, London, UK2005
        • Pereda J.
        • Ferragut V.
        • Quevedo J.M.
        • Guamis B.
        • Trujillo A.J.
        Heat damage evaluation in ultra-high pressure homogenized milk.
        Food Hydrocoll. 2009; 23: 1974-1979
        • Pereda J.
        • Jaramillo D.P.
        • Quevedo J.M.
        • Ferragut V.
        • Guamis B.
        • Trujillo A.J.
        Characterization of volatile compounds in ultra-high-pressure homogenized milk.
        Int. Dairy J. 2008; 18: 826-834
        • Pinho C.R.G.
        • Franchi M.A.
        • Tribst A.A.L.
        • Cristianini M.
        Effect of ultra high pressure homogenization on alkaline phosphatase and lactoperoxidase activity in raw skim milk.
        Procedia Food Sci. 2011; 1: 874-878
        • Pinho C.R.G.
        • Franchi M.A.
        • Tribst A.A.L.
        • Cristianinia M.
        Effect of high pressure homogenization process on Bacillus stearothermophilus and Clostridium sporogenes spores in skim milk.
        Procedia Food Sci. 2011; 1: 869-873
        • Qi P.X.
        • Ren D.
        • Xiao Y.
        • Tomasula P.M.
        Effect of homogenization and pasteurization on the structure and stability of whey protein in milk.
        J. Dairy Sci. 2015; 98 (25704975): 2884-2897
        • Sharabi S.
        • Okun Z.
        • Shpigelman A.
        Changes in the shelf life stability of riboflavin, vitamin C and antioxidant properties of milk after (ultra) high pressure homogenization: Direct and indirect effects.
        Innov. Food Sci. Emerg. Technol. 2018; 47: 161-169
        • Wisniewski J.R.
        • Zougman A.
        • Nagaraj N.
        • Mann M.
        Universal sample preparation method for proteome analysis.
        Nat. Methods. 2009; 6 (19377485): 359-362
        • Yang M.
        • Cao X.
        • Wu R.
        • Liu B.
        • Ye W.
        • Yue X.
        • Wu J.
        Comparative proteomic exploration of whey proteins in human and bovine colostrum and mature milk using iTRAQ-coupled LC-MS/MS.
        Int. J. Food Sci. Nutr. 2017; 68 (28276902): 671-681
        • Yang J.
        • An D.
        • Zhang P.
        Expression profiling of cassava storage roots reveals an active process of glycolysis/gluconeogenesis.
        J. Integr. Plant Biol. 2011; 53: 193-211
        • Yao Y.
        • Zhao G.
        • Zou X.
        • Huang L.
        • Wang X.
        Microstructural and lipid composition changes in milk fat globules during milk powder manufacture.
        RSC Advances. 2015; 5: 62638-62646
        • Ye A.
        • Cui J.
        • Dalgleish D.
        • Singh H.
        Effect of homogenization and heat treatment on the behavior of protein and fat globules during gastric digestion of milk.
        J. Dairy Sci. 2017; 100 (27837978): 36-47
        • Zhu Y.
        • Wang J.
        • Wang C.
        Research on the preparation, uniformity and stability of mixed standard substance for rapid detection of goat milk composition.
        Anim. Sci. J. 2018; 89 (29349891): 794-801
        • Zhu Y.
        • Xu H.
        • Chen H.
        • Xie J.
        • Shi M.
        • Shen B.
        • Deng X.
        • Liu C.
        • Zhan X.
        • Peng C.
        Proteomic analysis of solid pseudopapillary tumor of the pancreas reveals dysfunction of the endoplasmic reticulum protein processing pathway.
        Mol. Cell. Proteomics. 2014; 13 (24997997): 2593-2603