# Combining different proteomic approaches to resolve complexity of the milk protein fraction of dromedary, Bactrian camels and hybrids, from different regions of Kazakhstan

Alma Ryskaliyeva, Céline Henry, Guy Miranda, Bernard Faye, Gaukhar

Konuspayeva, Patrice Martin

► To cite this version:

Alma Ryskaliyeva, Céline Henry, Guy Miranda, Bernard Faye, Gaukhar Konuspayeva, et al.. Combining different proteomic approaches to resolve complexity of the milk protein fraction of dromedary, Bactrian camels and hybrids, from different regions of Kazakhstan. PLoS ONE, 2018, 13 (5), pp.e0197026. 10.1371/journal.pone.0197026. hal-02621893

HAL Id: hal-02621893

https://hal.inrae.fr/hal-02621893

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

# Combining different proteomic approaches to resolve complexity of the milk protein fraction of dromedary, Bactrian camels and hybrids, from different regions of Kazakhstan

Alma Ryskaliyeva¹*, Céline Henry², Guy Miranda¹, Bernard Faye³, Gaukhar Konuspayeva⁴, Patrice Martin¹

1 INRA, UMR GABI, AgroParisTech, Université Paris-Saclay, Jouy-en-Josas, France, 2 INRA, MICALIS Institute, Plateforme d'Analyse Protéomique Paris Sud-Ouest (PAPPSO), Université Paris-Saclay, Jouy-en-Josas, France, 3 CIRAD, UMR SELMET, France, 4 Al-Farabi Kazakh State National University, Biology department, Almaty, Kazakhstan

* alma.ryskaliyeva@inra.fr

<figure><img src="image_2.png" /></figure>

## Abstract

Nutritional suitability of milk is not only related to gross composition, but is also strongly affected by the microheterogeniety of the protein fraction. Hence, to go further into the evaluation of the potential suitability of non-bovine milks in human/infant nutrition it is necessary to have a detailed characterization of their protein components. Combining proven proteomic approaches (SDS-PAGE, LC-MS/MS and LC-ESI-MS) and cDNA sequencing, we provide here in depth characterization of the milk protein fraction of dromedary and Bactrian camels, and their hybrids, from different regions of Kazakhstan. A total 391 functional groups of proteins were identified from 8 camel milk samples. A detailed characterization of 50 protein molecules, relating to genetic variants and isoforms arising from post-translational modifications and alternative splicing events, belonging to nine protein families (κ-, α<sub>s1</sub>-, α<sub>s2</sub>-, β-; and γ-CN, WAP, α-LAC, PGRP, CSA/LPO) was achieved by LC-ESI-MS. The presence of two unknown proteins UP1 (22,939 Da) and UP2 (23,046 Da) was also reported as well as the existence of a β-CN short isoform (946 Da lighter than the full-length β-CN), arising very likely in both genetic variants (A and B) from proteolysis by plasmin. In addition, we report, for the first time to our knowledge, the occurrence of a α<sub>s2</sub>-CN phosphorylation isoform with 12P groups within two recognition motifs, suggesting thereby the existence of two kinase systems involved in the phosphorylation of caseins in the mammary gland. Finally, we demonstrate that genetic variants, which hitherto seemed to be species-specific (e.g. β-CN A for Bactrian and β-CN B for dromedary), are in fact present both in *Camel dromedarius* and *C. bactrianus*.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Ryskaliyeva A, Henry C, Miranda G, Faye B, Konuspayeva G, Martin P (2018) Combining different proteomic approaches to resolve complexity of the milk protein fraction of dromedary, Bactrian camels and hybrids, from different regions of Kazakhstan. PLoS ONE 13(5): e0197026. https://doi.org/10.1371/journal.pone.0197026

**Editor:** Thierry Rabilloud, Centre National de la Recherche Scientifique, FRANCE

**Received:** December 20, 2017

**Accepted:** April 25, 2018

**Published:** May 10, 2018

**Copyright:** © 2018 Ryskaliyeva et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All data underlying the study are within the paper and its Supporting Information files. Data from LC-MS/MS analyses will be available upon requirement by contacting the corresponding author by email alma.ryskaliyeva@inra.fr.

**Funding:** The study was carried out within the Bolashak International Scholarship of the first author, funded by the JSC «Center for International
Programs» (Kazakhstan). The research was partly supported by a grant from the Ministry of Education and Science of the Republic of Kazakhstan under name "Proteomic investigation of camel milk" #1729/GF4, which is duly appreciated. INRA (Jouy-en-Josas, France) for providing necessary facilities and technical support.

**Competing interests:** The authors have declared that no competing interests exist.

Introduction

According to the most recent statistics, the world camel population is estimated to be about 29 millions [1]. *Camelus dromedarius* is the most frequent and widespread domestic camel species composing 90% of the total camel population [2]. Camels have been domesticated in a number of arid regions, including Northern and Eastern Africa, the Arabian Peninsula and Central and South West Asia. *Camelus bactrianus* forms numerical inferiority, mostly inhabits in Mongolia, China, and Central Asia. Alternatively, there are also crossed camels (hybrids) which are found mainly in Russia, Iran, Turkmenistan, and in Kazakhstan.

Kazakhstan is a specific region where both domesticated species (*C. dromedarius* and *C. bactrianus*) along with wild Bactrian camels (*Camelus ferus*) are maintained in mixed herds [3]. There are about 160,000 camel heads reared in this country for milk production [1]. Camel milk is consumed as fresh milk and as a traditional fermented drink called *shubat*, which is very popular in Central Asia countries. Besides nutritional qualities, camel fresh and fermented milk have been reported to display potential health-promoting properties [4–9] which depend very heavily on its unique protein content.

Advanced improvement in proteomic techniques allow nowadays obtaining a precise image of the protein fraction of milk. Recently, proteomic approaches, based on mass spectrometry [10] and isobaric tag for relative and absolute quantification [11], have been used to analyze the proteome of dromedary camel milk and Bactrian camel milk whey, respectively. These techniques were useful to gain knowledge on the detection, quantification and characterization of camel milk proteins. These studies confirm that camel milk is a rich source of biologically active proteins and peptides [12], [13].

Whey proteins which were reported to display a wide range of bioactivities [14], including immuno-modulating [15], anti-carcinogenic [16], antibacterial, and antifungal activities [17], account for 20% of total camel milk proteins. Pattern-recognition proteins, such as the peptidoglycan recognition protein (PGRP), an intracellular component of neutrophils, modulate anti-inflammatory reaction of the immune response [18]. LTF interacts with lipopolysaccharides of Gram-negative bacteria whereas lysozyme C binds and hydrolyzes peptidoglycans, preferably of Gram-positive bacteria, but with a lower affinity than PGRP [19]. Present at a very low level in ruminant milks [20], PGRP has been detected in mammary secretions of porcine and camel [18] and was shown to participate in granule-mediated killing of gram-positive and negative bacteria [21]. Proteose peptone component 3 (PP3 or Lactophorin or GlyCAM1) plays an important immunological role in the lactating camel, to prevent the occurrence of mastitis, or for its newborn by inhibiting pathogen multiplication in the respiratory and gastrointestinal tracts of the suckling young [22]. Likewise, camel milk contains the whey acidic protein (WAP), also found in rodents and lagomorphs [23]. The biological function of this protein is unknown. However, proteins such as elafin and antileukoproteinase 1, containing WAP domains, are known to function as protease inhibitor involved in the immune defence of multiple epithelia and has been identified as candidate molecular markers for several cancers [24].

As in cow milk, *ca*. 80% of the total protein fraction of camel milk are represented by caseins (CN) that are synthesized under multi-hormonal control in the mammary gland of mammals. Associated with amorphous calcium phosphate nanoclusters they form large and stable colloidal aggregates, the so-called CN micelles, which figure as calcium-transport vehicles. These CN micelles provide neonates with calcium at a very high concentration, which is achieved during their packaging in the secretion pathway [25]. Recently it was reported that $\alpha_{s1}$- and $\alpha_{s2}$-CN display molecular chaperone-like activity inhibiting CN aggregation and triggering micelle structure [26].
However, there is no comprehensive investigation on milk protein variations and variability in composition between individual camels. In addition, proteomic studies did not consider the molecular diversity of each type of protein, arising from genetic polymorphisms (mutations), defects in the processing of primary transcripts and post-translational modifications (PTM) such as phosphorylation, factors that significantly have a pronounced impact on protein structure, and finally on milk properties. Milk protein polymorphism is a unique biological paradigm that could help to understand CN intracellular transport, micelle formation and organization, biodiversity and evolution [27], the release of bioactive peptides with implications in human health [28].

Therefore, to gain an insight into the molecular diversity of camel milk proteins, we design a comprehensive strategy combining classical (SDS-PAGE) and advanced proteomic approaches (LC-MS/MS, LC-ESI-MS), as well as cDNA sequencing. Here we report a complete profiling of the milk protein fraction of Bactrian and dromedary camels from Kazakhstan, including a detailed characterization of camel CN and whey proteins including variants related to genetic polymorphisms, splicing defects, phosphorylation levels. In addition, we introduce a reference point for further investigation in camel milk protein polymorphism.

# Materials and methods

## Ethics statements

All animal studies were carried out in compliance with European Community regulations on animal experimentation (European Communities Council Directive 86/609/EEC) and with the authorization of the Kazakh Ministry of Agriculture. Milk sampling was performed in appropriate conditions supervised by a veterinary accredited by the French Ethics National Committee for Experimentation on Living Animals. No endangered or protected animal species were involved in this study. No specific permissions or approvals were required for this study with the exception of the rules of afore-mentioned European Community regulations on animal experimentation, which were strictly followed.

## Milk samples collection and preparation

In total 179 raw milk samples (Table 1) were collected during morning milking on healthy dairy camels belonging to two camel species: *C. bactrianus* (n = 72) and *C. dromedarius* (n = 65), and their hybrids (n = 42), at different lactation stages, ranging between 30 and 90 days postpartum. Bactrian camels were originating from Kazakh type whereas dromedary camels were from Turkmen Arvana breed. Unfortunately, the information about the nature and the level of hybridization of hybrids was not available. All species are well adapted to the local environment of Kazakhstan.

Camels grazed on four various natural pastures with the distance more than 3,500 kms between the regions at extreme points of Kazakhstan: Almaty (AL) at the foot of Tien Shan Mountain, Shymkent (SH) along deserts Kyzylkum and Betpak-Dala, Kyzylorda (KZ) on the

<table><caption>Table 1. Camel milk samples collected (n = 179) in the 3 species of the 4 regions of Kazakhstan.</caption><thead><tr><td>ID</td><td>Region</td><td>Coding</td><td>Bactrian (B)</td><td>Dromedary (D)</td><td>Hybrid (H)</td><td>Total number of camels for each region</td></tr></thead><tbody><tr><td>1</td><td>Almaty</td><td>AL B/D/H</td><td>13</td><td>20</td><td>1</td><td>34</td></tr><tr><td>2</td><td>Shymkent</td><td>SH B/D/H</td><td>20</td><td>21</td><td>20</td><td>61</td></tr><tr><td>3</td><td>Kyzylorda</td><td>KZ B/D/H</td><td>18</td><td>16</td><td>20</td><td>54</td></tr><tr><td>4</td><td>Atyrau</td><td>ZKO B/D/H</td><td>21</td><td>8</td><td>1</td><td>30</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0197026.t001
edge of the steppe, and Atyrau (ZKO) at the mouth of the Caspian Sea (Fig 1). Whole-milk samples were centrifuged at 2,500 g for 20 min at 4°C (Allegra X-15R, Beckman Coulter, France) to separating fat from skimmed milk. Samples were quickly frozen and stored at -80°C (fat) and -20°C (skimmed milk) until analysis.

## Selection of milk samples for analysis

Of the 179 milk samples collected, 63, including *C. bactrianus* (n = 19), *C. dromedarius* (n = 20), and hybrids (n = 24) from four different regions of Kazakhstan were selected for SDS-PAGE analysis (Fig 2). Each Bactrian and dromedary camel group formed by 5 animals, except Bactrians of Atyrau regions (n = 4). For hybrids, there were 4 groups comprising 10 animals (Kyzylorda and Shymkent regions), whereas there were only 1 and 3 animals for Almaty and Atyrau regions, respectively. This selection was based on lactation stages and number of parities (from 2 to 14) of each camel group composed by the species and grazing regions. It should be emphasized that data available on animals: breed, age, lactation stage and calving number, were estimated by a local veterinarian, since no registration of camels in farms is maintained. Due to the lack of sufficient information, dromedary milk samples (n = 5) from Almaty region were excluded from subsequent analyses. Then, 8 of the 58 remaining milk samples from three different regions (*C. bactrianus*, n = 3, *C. dromedarius*, n = 3, and hybrids, n = 2) exhibiting the most representative SDS-PAGE patterns were analyzed by LC-MS/MS after a tryptic digestion of excised gel bands. Additionally, 30 milk samples (*C. bactrianus*,

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Geographical location of camel milk sampling. Reprinted from http://camelides.cirad.fr/fr/science/pdf/presentation_these_konuspayeva.pdf under a CC BY license, with permission from Konuspayeva Gaukhar, original copyright 2007. https://www.cia.gov/library/publications/the-world-factbook/geos/kz.html https://upload.wikimedia.org/wikipedia/commons/thumb/b/b0/Kazakhstan_on_the_globe_%28Eurasia_centered%29.svg/512px-Kazakhstan_on_the_globe_%28Eurasia_centered%29.svg.png.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0197026.g001
<figure><img src="image_5.png" /><figcaption>Fig 2. Diagram of the experimental scheme designed for quantification and identification of camel milk proteins.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0197026.g002

n = 10; *C. dromedarius*, n = 10; hybrids, n = 10), taken from the 63 milks analyzed by SDS-PAGE, were analyzed by LC-ESI-MS (Bruker Daltonics).

## Coomassie blue (Bradford) protein assay

To estimate the concentration of total protein in a milk sample the Coomassie Blue Protein Assay was used [29]. Absorbance at 590 nm was measured using the UV-Vis spectrophotometer (UVmini-1240, Shimadzu). The reference standard curve was done with commercial bovine serum albumin (BSA) powder dissolved in MilliQ water and diluted to a concentration
of 1 mg/mL. Series of dilutions (0.1, 0.2, 0.4, 0.6, and 0.8 µg/µL) were prepared from the stock solution, in duplicate to ensure the protein concentration is within the range of the assay.

## 1D sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE)

Both major and low-abundant proteins resolved by SDS-PAGE were identified after excision by mass analysis of the tryptic hydrolysate. The method used in the study was based on that from Laemmli [30]. Twenty-five micrograms of each individual skimmed milk sample were loaded into 12.5% acrylamide resolving gel and subjected to electrophoresis. Samples were prepared with Laemmli Lysis-Buffer (Sigma-Aldrich). Separations were performed in a vertical electrophoresis apparatus (Bio-Rad, Marnes-la-Coquette, France). After GelCode Blue Safe Protein staining and gel scanning using Image Scanner iii (Epson Expression™ 10,000 XL, Sweden), resolved bands were excised from the gel and submitted to digestion by trypsin. Thereafter, tryptic peptides were analyzed by LC-MS/MS.

## Identification of proteins by LC-MS/MS analysis

In order to identify the main protein contained in each electrophoretic band, mono dimensional electrophoresis (1D SDS-PAGE) followed by trypsin digestion and by LC-MS/MS analysis, was used essentially as described [31]. Briefly, after a 10 cm migration of samples in such an 1D SDS-PAGE, the 16 main electrophoretic bands (1.5 mm³) were cut on each gel lane, transferred into 96-well microtiter plates (FrameStar, 4titude, 0750/Las). Reduction of disulfide bridges of proteins was carried out by incubating at 37°C for one hour with dithiothreitol (DTT, 10 mM, Sigma), meanwhile the alkylation of free cysteinyl residues with iodoacetamide (IAM, 50 mM, Sigma) at room temperature for 45 min in total obscurity. After gel pieces were washed twice, first, with 100 µL 50% ACN/50 mM NH₄HCO₃ and then with 50 µL ACN, they were finally dried. The hydration was performed at 37°C overnight using digestion buffer 400 ng lys-C protease + trypsin. Hereby, peptides were extracted with 50% ACN/0.5% TFA and then with 100% ACN. Peptide solutions were dried in a concentrator and finally dissolved into 70 µL 2% ACN in 0.08% TFA.

The identification of peptides was obtained using UltiMate™ 3000 RSLCnano System (Thermo Fisher Scientific) coupled either to LTQ Orbitrap XL™ Discovery mass spectrometer or QExactive (Thermo Fischer Scientific).

Four µL of each sample was injected with flow of 20 µL/min on a precolumn cartridge (stationary phase: C18 PepMap 100, 5 µm; column: 300 µm x 5 mm) and desalted with a loading buffer 2% ACN and 0.08% TFA. After 4 min, the precolumn cartridge was connected to the separating RSLC PepMap C18 column (stationary phase: RSLC PepMap 100, 2 µm; column: 75 µm x 150 mm). Elution buffers were A: 2% ACN in 0.1% formic acid (HCOOH) and B: 80% ACN in 0.1% HCOOH. The peptide separation was achieved with a linear gradient from 0 to 35% B for 34 min at 300 nL/min. One run took 42 min, including the regeneration and the equilibration steps at 98% B.

Peptide ions were analyzed using Xcalibur 2.1 with the following machine set up in CID mode: 1) full MS scan in Orbitrap with a resolution of 15 000 (scan range [m/z] = 300–1600) and 2) top 8 in MS/MS using CID (35% collision energy) in Ion Trap. Analyzed charge states were set to 2–3, the dynamic exclusion to 30 s and the intensity threshold was fixed at 5.0 x 10².

Raw data were converted to mzXML by MS convert (ProteoWizard version 3.0.4601). UniProtKB Cetartiodactyla database was used (157,113 protein entries, version 2015), in conjunction with contaminant databases were searched by algorithm X!TandemPiledriver (version 2015.04.01.1) with the software X!TandemPipeline (version 3.4) developed by the PAPPSO
platform (http://pappso.inra.fr/bioinfo/). The protein identification was run with a precursor mass tolerance of 10 ppm and a fragment mass tolerance of 0.5 Da. Enzymatic cleavage rules were set to trypsin digestion ("after R and K, unless P follows directly after") and no semi-enzymatic cleavage rules were allowed. The fix modification was set to cysteine carbamido methylation and methionine oxidation was considered as a potential modification. Results were filtered using inbuilt X!TandemParser with peptide *E*-value of 0.05, a protein *E*-value of -2.6, and a minimum of two peptides.

## LC-ESI-MS

Fractionation of camel milk proteins and determination of their molecular masses, performed by coupling RP-HPLC to ESI-MS (micrOTOF™ II focus ESI-TOF mass spectrometer; Bruker Daltonics), were essentially as described [31]. In total 20 µL of skimmed milk samples were first clarified by the addition of 230 µL of clarification solution 0.1 M bis-Tris buffer pH 8.0, containing 8 M urea, 1.3% trisodium citrate, and 0.3% DTT. Clarified milk samples (25 µL) were directly injected onto a Biodiscovery C5 reverse phase column (300 Å pore size, 3 µm, 150 x 2.1 mm; Supelco, France). The mobile phase of the column corresponded to a gradient mixture of Solvent A (H₂O/TFA 100:0.25, v/v) and Solvent B (ACN/TFA 100:0.20, v/v). Elution was achieved using a linear gradient from 5% to 27% B in 20 min, from 27% to 33% B in 0.1 min, from 33% to 34% B in 11.1 min, from 34% to 40% B in 0.1 min, from 40% to 41% B in 14.9 min, and from 41% to 90% B in 0.1 min. This gradient elution was followed by an isocratic elution at 90% B for 4.9 min, and a linear return to 5% B in 0.1 min. The temperature of the column was adjusted to 52°C and the flow rate to 0.2 mL/min. Eluted peaks were detected by UV-absorbance at 214 nm. The liquid effluent was introduced to the mass spectrometer. Positive ion mode was used, and mass scans were acquired over a mass-to-charge ratio (m/z) ranging between 600 and 3000 Da.

The LC/MS system was controlled by the HyStar software (Bruker Daltonics). Peak profiles from UV 214 nm and Extracted Ion Chromatograms (EIC), multicharged ion spectra, deconvoluted spectra and determination of masses were obtained with DataAnalysis Version 4.0 SP1 software (Bruker Daltonics).

## Milk fat globule collection and RNA extraction

Milk was centrifuged at 2,500 g for 20 min to pellet somatic cells (SC) and to separate the upper milk fat globule (MFG) fraction. The MFG fraction was mixed with Trizol LS and heated briefly at 30°C while shaking, to emulsify fat. Total RNA was extracted from milk fat using Trizol (Invitrogen) following the protocol from the manufacturer, as described in Breton *et al.* [32].

## First-strand cDNA synthesis and PCR amplification

First-strand cDNA was synthesized from 5 to 10 ng of total RNA primed with oligo(dT)20 and random primers (3:1, vol/vol) using Superscript III reverse transcriptase (Invitrogen Life Technologies Inc., Carlsbad, CA) according to the manufacturer's instructions. One microliter of 2 U/µL RNase H (Invitrogen Life Technologies) was then added and the reaction mix was incubated for 20 min at 37°C to remove RNA from heteroduplexes. Single-strand cDNA thus obtained was stored at -20°C. cDNA samples covering the entire coding regions of caseins were amplified. PCR was performed in an automated thermocycler GeneAmp® PCR System 2,400 (Perkin-Elmer, Norwalk, USA) with GoTaq® G2 Flexi DNA Polymerase Kit (Promega Corporation, USA). Reactions were carried out with 0.2 mL thin-walled PCR tubes with flat cap strips (Thermo Scientific, UK), in 50 µL volumes containing 5X Green or Colorless
GoTag® Flexi Buffer, MgCl₂ Solution 25 mM, PCR Nucleotide Mix 10 mM each, GoTag® G2 Flexi DNA Polymerase (5 U/µL), 10 mM each oligonucleotide primer, template DNA and nuclease-free water, up to the final volume. Primer pairs, purchased from Eurofins (Eurofins genomics, Germany), were designed using published *Camelus* nucleic acid sequence. Sequencing of PCR fragments was performed with primer pairs used for PCR and sequenced from both strands, according to the Sanger method by Eurofins.

# Results

## Total protein content

Using the Bradford assay for estimating the protein concentration in milk samples, we observed that the highest protein concentration occurred with Bactrian camel milk samples, but the difference was slight comparing with crossed camel species. The total protein value in raw camel milk from Shymkent region was estimated to be ca. 33 g/L (33.15 ± 6.64 g/L) for *C. bactrianus* (n = 5), and 31 g/L (30.83 ± 5.82 g/L) for *C. dromedarius* (n = 7), whereas hybrids (n = 9) displayed an intermediate value 31.5 g/L (31.43 ± 4.56 g/L). On average, Bactrian milk was considered to have a higher total protein content than that of Dromedary [33] and hybrid milks. Our results are in agreement with data reported previously by Konuspayeva et al. [34]. No significant differences were found across species from different geographical locations.

## Identification of main milk proteins from 1D SDS-PAGE by LC-MS/MS

After first adjusting protein concentrations at the same value, 63 individual camel milk samples were separated onto SDS-PAGE. The comparative analysis of whole milk samples by SDS-PAGE displayed rather similar electrophoretic profiles with related migration characteristics and the same apparent molecular weights between individual milk samples of different species and regions. A typical gel pattern from which proteins were identified in individual *C. bactrianus*, *C. dromedarius* and hybrid milk samples of Kyzylorda region is shown in Fig 3.

Sixteen main bands relatively well-resolved were excised from the electrophoretic pattern. The most intense band observed around 26 kDa was identified as β-CN. Quantitative analyses on camel milk proteins carried out before have demonstrated significantly higher amounts of β-CN compared to the homologous bovine CN [35]. The most representative other bands were characterized as being: WAP (12.5 kDa), α-LAC (14.3 kDa), GlyCAM 1 (15.4 kDa and 17.2 kDa), κ-CN (20.3 kDa), PGRP (21.3 kDa), αs2-CN (22.9 kDa), αs1-CN (25.7 kDa), neutrophil gelatinase (28.3 kDa), lipoprotein lipase (46.5 kDa), perilipin-2 (47.2 kDa), butyrophilin (51.0 kDa), amine oxidase (55.3 kDa), lactadherin (56.2 kDa), heat shock protein (70.0 kDa), LTF (77.1 kDa), lactoperoxidase (87.7 kDa), and xanthine oxidase (150 kDa). Masses mentioned above correspond to theoretical masses of proteins identified on the basis of tryptic profiles after LC-MS/MS analysis. Globally, the electrophoretic patterns of Kazakh camel milk samples agree with those reported recently for Israeli and Tunisian camel milk samples [36], [37]. However, surprisingly the prominent fact was the apparent absence in Kazakh milk samples of camel serum albumin (CSA), the major whey protein with a molecular mass equal to 66.0 kDa in camel colostrum [36]. By contrast, this protein has been successfully identified, with the best E-value, in Tunisian fresh milk samples [37].

## Qualitative proteome of camel skimmed milk by LC-MS/MS

We took advantage of LC-MS/MS analysis to identify proteins in electrophoretic bands to go further into the description of the protein fraction of camel milk. Indeed, for each band analyzed by LC-MS/MS, between 10 and 70 different proteins were identified. In such a way,
<figure><img src="image_6.png" /><figcaption>**Fig 3.** 1D SDS-PAGE pattern of *C. bactrianus* (KZB), *C. dromedarius* (KZD) and hybrid (KZH) skimmed milk samples of Kyzylorda (KZ) region. Red frames and black boxes aligned correspond to electrophoretic bands that were excised from the gel and subsequently analyzed for protein identification, after tryptic digestion, by LC-MS/MS. Molecular weight markers from 210 to 8 kDa are indicated at the right of the gel.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0197026.g003

using UniprotKB taxonomy cetartiodactyla (SwissProt + Trembl) database, a total of 391 func-
tional groups of proteins (proteins belonging to a same group share common peptides) were
identified after LC-MS/MS analysis of 8 camel milk samples (S1 Table). A set of 235 proteins
was observed as common to the 8 milk samples. As example, a list of the first 70 common pro-
teins found in milk samples of the three species from Shymkent region is given in Table 2.

Eight proteins were identified as authentically matching with proteins in *C. dromedarius*
protein database, two with *C. bactrianus* protein database, 46 with *C. ferus* protein database,
and the remaining (n = 14) with the other mammalian species such as, *Lama guanicoe*, *Bos tau*-
*rus*, *Sus scrofa* and *Ovis aries* protein databases. Immune-related proteins such as GlyCAM1,
lactadherin (MFG-E8), and LTF, as well as milk fat globule membrane (MFGM)-enriched pro-
teins such as xanthine oxidase (XO), butyrophilin (BTN), actin, ras-related protein Rab-18,
ADP-ribosylation factor 1, tyrosine-protein kinase and GTP-binding protein SAR1b, were
detected. Likewise, proteins originating from blood such as serpin A3-1, apolipoprotein A-1,
α-1-antitrypsin like protein, α-1-acid glycoprotein, β-2-microglobulin, complement C3-like
protein were found in all milk samples analyzed.

Camel milk protein profiling by LC-ESI-MS

Thirty individual milk samples, including *C. bactrianus* (n = 10), *C. dromedarius* (n = 10), and hybrids (n = 10) taken from the 58 milk samples analyzed in SDS-PAGE were submitted to LC-ESI-MS analysis. Milk proteins separated by RP-HPLC were identified based on their molecular mass, arising from ESI-MS. Putative genetic variants and post-translational (glycosylation and phosphorylation) isoforms were determined by deconvoluting multiple charged ion spectra in a real mass scale. Knowing their primary structures, it is possible to determine molecular masses of non post-translationally modified proteins, and then we can precisely know the mass of phosphorylation isoforms resulting from the addition of phosphate groups
<table><caption>Table 2. Top 70 proteins identified by LC-MS/MS from individual *C. bactrianus* (B), *C. dromedarius* (D) and hybrid (H) milk samples of Shymkent region.</caption><thead><tr><th rowspan="2">ID</th><th rowspan="2">Accession</th><th rowspan="2">Description</th><th rowspan="2">Mr</th><th colspan="4">(-) log E-value</th><th colspan="4">Coverage, %</th><th colspan="4">Number of Spectra*</th></tr><tr><th>kDa</th><th>B</th><th>D</th><th>H</th><th>B</th><th>D</th><th>H</th><th>B</th><th>D</th><th>H</th></tr></thead><tbody><tr><td>1</td><td>O97943-2</td><td>Short isoform of Alpha-S1-casein (C. dromedarius)</td><td>25,7</td><td>453,63</td><td>535,3</td><td>585,79</td><td>88</td><td>92</td><td>92</td><td>574</td><td>679</td><td>842</td></tr><tr><td>2</td><td>A0A077SL35</td><td>Beta-casein (C. bactrianus)</td><td>26,1</td><td>299,05</td><td>285,89</td><td>395,72</td><td>77</td><td>79</td><td>79</td><td>542</td><td>507</td><td>745</td></tr><tr><td>3</td><td>O97944</td><td>Alpha-S2-casein (C. dromedarius)</td><td>22,9</td><td>285,16</td><td>320,65</td><td>276,39</td><td>70</td><td>70</td><td>66</td><td>357</td><td>451</td><td>376</td></tr><tr><td>4</td><td>W6GH05</td><td>Lactoferrin (C. dromedarius)</td><td>77,1</td><td>723,25</td><td>557,5</td><td>1174,16</td><td>85</td><td>80</td><td>88</td><td>253</td><td>157</td><td>829</td></tr><tr><td>5</td><td>P15522-2</td><td>Isoform B of Glycosylation-dependent cell adhesion molecule 1 (C. dromedarius)</td><td>15,4</td><td>112,39</td><td>189,45</td><td>192,31</td><td>64</td><td>70</td><td>69</td><td>220</td><td>359</td><td>358</td></tr><tr><td>6</td><td>L0P304</td><td>Kappa-casein (C. bactrianus)</td><td>20,3</td><td>127,64</td><td>149,66</td><td>194,03</td><td>51</td><td>53</td><td>51</td><td>217</td><td>222</td><td>329</td></tr><tr><td>7</td><td>S9WF76</td><td>Lactadherin-like protein (C. ferus)</td><td>45,6</td><td>263,03</td><td>350,07</td><td>501,73</td><td>51</td><td>53</td><td>57</td><td>162</td><td>236</td><td>436</td></tr><tr><td>8</td><td>P00710</td><td>Alpha-lactalbumin (C. dromedarius)</td><td>14,3</td><td>288,2</td><td>277,9</td><td>329,27</td><td>83</td><td>80</td><td>83</td><td>161</td><td>159</td><td>175</td></tr><tr><td>9</td><td>Q9GK12</td><td>Peptidoglycan recognition protein 1 (C. dromedarius)</td><td>21,3</td><td>308,94</td><td>359,49</td><td>336,8</td><td>73</td><td>73</td><td>79</td><td>112</td><td>156</td><td>140</td></tr><tr><td>10</td><td>S9Z0L8</td><td>Amine oxidase [flavin-containing] (C. ferus)</td><td>55,3</td><td>233,97</td><td>231,94</td><td>277,32</td><td>74</td><td>70</td><td>75</td><td>110</td><td>138</td><td>190</td></tr><tr><td>11</td><td>S9Y4T1</td><td>Xanthine dehydrogenase/oxidase (C. ferus)</td><td>150</td><td>270,07</td><td>435,86</td><td>330,26</td><td>40</td><td>45</td><td>41</td><td>79</td><td>138</td><td>121</td></tr><tr><td>12</td><td>S9X3X3</td><td>Butyrophilin subfamily 1 member A1 (C. ferus)</td><td>51,09</td><td>95,37</td><td>140,19</td><td>162,3</td><td>47</td><td>50</td><td>51</td><td>61</td><td>107</td><td>111</td></tr><tr><td>13</td><td>S9X4G0</td><td>Neutrophil gelatinase-associated lipocalin-like protein (C. ferus)</td><td>28,3</td><td>104,88</td><td>98,15</td><td>165,14</td><td>48</td><td>45</td><td>62</td><td>44</td><td>33</td><td>82</td></tr><tr><td>14</td><td>S9X1L5</td><td>Lipoprotein lipase isoform 3 (Fragment) (C. ferus)</td><td>46,5</td><td>107,66</td><td>133,72</td><td>54,73</td><td>43</td><td>55</td><td>27</td><td>43</td><td>56</td><td>29</td></tr><tr><td>15</td><td>S9YK74</td><td>Perilipin (C. ferus)</td><td>47,2</td><td>116,47</td><td>169,78</td><td>167,59</td><td>56</td><td>60</td><td>52</td><td>39</td><td>67</td><td>46</td></tr><tr><td>16</td><td>P09837</td><td>Whey acidic protein (C. dromedarius)</td><td>12,5</td><td>42,97</td><td>59,92</td><td>82,11</td><td>55</td><td>75</td><td>84</td><td>36</td><td>43</td><td>73</td></tr><tr><td>17</td><td>S9X4X6</td><td>Uncharacterized protein (C. ferus)</td><td>43,1</td><td>47,87</td><td>14</td><td>17,5</td><td>30</td><td>13</td><td>16</td><td>36</td><td>9</td><td>16</td></tr><tr><td>18</td><td>S9X7Q1</td><td>Lactoperoxidase isoform 1 preproprotein (C. ferus)</td><td>87,7</td><td>83,1</td><td>62,59</td><td>65,16</td><td>26</td><td>22</td><td>20</td><td>25</td><td>19</td><td>17</td></tr><tr><td>19</td><td>S9XDK9</td><td>Complement C3-like protein (C. ferus)</td><td>262,7</td><td>87,48</td><td>44,5</td><td>292,06</td><td>12</td><td>4</td><td>24</td><td>25</td><td>8</td><td>69</td></tr><tr><td>20</td><td>S9XR87</td><td>Beta-2-microglobulin (C. ferus)</td><td>14,8</td><td>36,87</td><td>36,15</td><td>40,76</td><td>43</td><td>43</td><td>43</td><td>23</td><td>28</td><td>44</td></tr><tr><td>21</td><td>O18831</td><td>Growth/differentiation factor 8 (S. scrofa)</td><td>42,7</td><td>44,62</td><td>49,29</td><td>67,44</td><td>28</td><td>24</td><td>32</td><td>23</td><td>22</td><td>17</td></tr><tr><td>22</td><td>P68103</td><td>Elongation factor 1-alpha 1 (B. taurus)</td><td>50</td><td>58,23</td><td>42,59</td><td>78,3</td><td>31</td><td>27</td><td>31</td><td>20</td><td>12</td><td>24</td></tr><tr><td>23</td><td>S9YCI6</td><td>Peptidyl-prolyl cis-trans isomerase (C. ferus)</td><td>23,8</td><td>44,54</td><td>43,1</td><td>71,24</td><td>54</td><td>54</td><td>57</td><td>20</td><td>21</td><td>22</td></tr><tr><td>24</td><td>S9XP75</td><td>Monocyte differentiation antigen CD14 (C. ferus)</td><td>29,7</td><td>49,75</td><td>44,58</td><td>73,91</td><td>27</td><td>27</td><td>27</td><td>17</td><td>16</td><td>22</td></tr><tr><td>25</td><td>S9WCV2</td><td>Sulfhydryl oxidase (C. ferus)</td><td>72,2</td><td>31,64</td><td>76,49</td><td>14,43</td><td>18</td><td>27</td><td>8</td><td>16</td><td>26</td><td>6</td></tr><tr><td>26</td><td>S9YC53</td><td>Alpha-1-antitrypsin-like protein (C. ferus)</td><td>51,9</td><td>64,76</td><td>64,84</td><td>51,84</td><td>28</td><td>31</td><td>32</td><td>14</td><td>21</td><td>15</td></tr><tr><td>27</td><td>S9YS49</td><td>Putative E3 ubiquitin-protein ligase Roquin (C. ferus)</td><td>158,5</td><td>56,65</td><td>90,28</td><td>8,84</td><td>10</td><td>14</td><td>2</td><td>14</td><td>29</td><td>3</td></tr><tr><td>28</td><td>T0NN97</td><td>Uncharacterized protein (C. ferus)</td><td>151,4</td><td>38,86</td><td>53,95</td><td>36,72</td><td>12</td><td>13</td><td>11</td><td>13</td><td>16</td><td>15</td></tr><tr><td>29</td><td>A0A0F6YEF6</td><td>Anti-HCV NS3/4A serine protease immunoglobulin heavy chain (Fragment) (C. dromedarius)</td><td>13,4</td><td>38,24</td><td>35</td><td>52,95</td><td>24</td><td>24</td><td>24</td><td>12</td><td>10</td><td>26</td></tr><tr><td>30</td><td>S9X9X0</td><td>Vitelline membrane outer layer protein 1-like protein (C. ferus)</td><td>24,2</td><td>44,94</td><td>52,11</td><td>75,95</td><td>51</td><td>50</td><td>53</td><td>12</td><td>14</td><td>31</td></tr><tr><td>31</td><td>S9X358</td><td>Tissue alpha-L-fucosidase (C. ferus)</td><td>42,7</td><td>25,83</td><td>23,27</td><td>24,46</td><td>23</td><td>18</td><td>18</td><td>12</td><td>12</td><td>7</td></tr><tr><td>32</td><td>S9WS72</td><td>Sodium-dependent phosphate transport protein 2B-like protein (C. ferus)</td><td>75,4</td><td>24,99</td><td>15,41</td><td>24,7</td><td>18</td><td>7</td><td>9</td><td>11</td><td>8</td><td>8</td></tr><tr><td>33</td><td>S9X2V0</td><td>UDP-Gal:betaGlcNAc beta 1,4-galactosyltransferase 1, membrane-bound form-like protein (C. ferus)</td><td>37,3</td><td>30,18</td><td>13,09</td><td>38,99</td><td>23</td><td>11</td><td>33</td><td>11</td><td>7</td><td>15</td></tr><tr><td>34</td><td>S9XSQ6</td><td>Vitamin D-binding protein-like protein (C. ferus)</td><td>49,3</td><td>41,66</td><td>109,25</td><td>18,01</td><td>24</td><td>45</td><td>8</td><td>10</td><td>30</td><td>3</td></tr><tr><td>35</td><td>S9WGZ9</td><td>Rab GDP dissociation inhibitor beta isoform 1 (C. ferus)</td><td>35,7</td><td>27,75</td><td>22,32</td><td>55,78</td><td>32</td><td>25</td><td>43</td><td>10</td><td>8</td><td>15</td></tr><tr><td>36</td><td>P19120</td><td>Heat shock cognate 71 kDa protein (B. taurus)</td><td>71,1</td><td>37,14</td><td>60,92</td><td>111,21</td><td>17</td><td>24</td><td>35</td><td>9</td><td>16</td><td>30</td></tr><tr><td>37</td><td>S9XA25</td><td>Ezrin isoform 5-like protein (Fragment) (C. ferus)</td><td>71</td><td>24,32</td><td>28,23</td><td>47,46</td><td>13</td><td>10</td><td>20</td><td>9</td><td>7</td><td>13</td></tr><tr><td>38</td><td>S9XI30</td><td>Uncharacterized protein (C. ferus)</td><td>22</td><td>25,16</td><td>66,23</td><td>20,66</td><td>38</td><td>54</td><td>39</td><td>9</td><td>31</td><td>8</td></tr><tr><td>39</td><td>S9Y2X0</td><td>Platelet glycoprotein 4 (C. ferus)</td><td>38,3</td><td>28,14</td><td>40,27</td><td>55,15</td><td>15</td><td>17</td><td>29</td><td>9</td><td>8</td><td>17</td></tr><tr><td>40</td><td>G9F6X8</td><td>Protein disulfide-isomerase (S. scrofa)</td><td>56,3</td><td>20,09</td><td>20,62</td><td>21,76</td><td>17</td><td>15</td><td>9</td><td>9</td><td>10</td><td>3</td></tr><tr><td>41</td><td>S9WZP7</td><td>Serpin A3-8 (C. ferus)</td><td>74,9</td><td>15,26</td><td>28,06</td><td>72,32</td><td>11</td><td>12</td><td>20</td><td>8</td><td>12</td><td>28</td></tr><tr><td>42</td><td>A7MBJ4</td><td>Receptor-type tyrosine-protein phosphatase F (B. taurus)</td><td>211,1</td><td>23,69</td><td>21,19</td><td>20,54</td><td>6</td><td>6</td><td>4</td><td>8</td><td>8</td><td>8</td></tr><tr><td>43</td><td>S9YFG2</td><td>Complement factor D (C. ferus)</td><td>38,1</td><td>37,99</td><td>33,7</td><td>44,44</td><td>14</td><td>14</td><td>14</td><td>8</td><td>9</td><td>13</td></tr><tr><td>44</td><td>S9WPL9</td><td>Uncharacterized protein (C. ferus)</td><td>84,6</td><td>31,09</td><td>39,23</td><td>57,64</td><td>12</td><td>10</td><td>17</td><td>7</td><td>6</td><td>18</td></tr></tbody></table>

(Continued)
<table><caption>**Table 2.** (Continued)</caption><thead><tr><th rowspan="2">ID</th><th rowspan="2">Accession</th><th rowspan="2">Description</th><th rowspan="2">Mr</th><th colspan="3">(-) log E-value</th><th colspan="3">Coverage, %</th><th colspan="3">Number of Spectra*</th></tr><tr><th>kDa</th><th>B</th><th>D</th><th>H</th><th>B</th><th>D</th><th>H</th><th>B</th><th>D</th><th>H</th></tr></thead><tbody><tr><td>45</td><td>S9XE02</td><td>Beta-2-glycoprotein 1-like protein (C. ferus)</td><td>30,9</td><td>46,22</td><td>64,63</td><td>48,09</td><td>42</td><td>43</td><td>27</td><td>7</td><td>12</td><td>5</td></tr><tr><td>46</td><td>B8XH67</td><td>Na(+)/H(+) exchange regulatory cofactor NHE-RF (S. scrofa)</td><td>39,2</td><td>20,94</td><td>37,33</td><td>32,94</td><td>22</td><td>22</td><td>30</td><td>7</td><td>9</td><td>11</td></tr><tr><td>47</td><td>Q28452</td><td>Quinone oxidoreductase (L. guanicoe)</td><td>35,1</td><td>23,01</td><td>24,97</td><td>51,04</td><td>28</td><td>25</td><td>43</td><td>7</td><td>5</td><td>11</td></tr><tr><td>48</td><td>S9YU13</td><td>Glutathione S-transferase-like protein (C. ferus)</td><td>27,7</td><td>34,68</td><td>8,59</td><td>22</td><td>21</td><td>15</td><td>19</td><td>7</td><td>3</td><td>4</td></tr><tr><td>49</td><td>P00727-2</td><td>Isoform 2 of Cytosol aminopeptidase (B. taurus)</td><td>53,9</td><td>14,58</td><td>13,07</td><td>35,48</td><td>15</td><td>9</td><td>18</td><td>6</td><td>4</td><td>7</td></tr><tr><td>50</td><td>S9WUC8</td><td>Ig kappa chain V-II region RPMI 6410-like protein (C. ferus)</td><td>26,7</td><td>22,72</td><td>12,39</td><td>32,25</td><td>19</td><td>10</td><td>27</td><td>6</td><td>4</td><td>15</td></tr><tr><td>51</td><td>S9WRY0</td><td>L-lactate dehydrogenase B chain isoform 1-like protein (C. ferus)</td><td>30,2</td><td>25,51</td><td>28,95</td><td>37,66</td><td>22</td><td>22</td><td>22</td><td>6</td><td>7</td><td>6</td></tr><tr><td>52</td><td>S9X5V9</td><td>Fc of IgG binding protein (Fragment) (C. ferus)</td><td>254,9</td><td>20,45</td><td>65,48</td><td>10,4</td><td>3</td><td>9</td><td>2</td><td>6</td><td>16</td><td>4</td></tr><tr><td>53</td><td>B5B0D4</td><td>Major allergen beta-lactoglobulin (B. taurus)</td><td>19,8</td><td>9,94</td><td>54,67</td><td>12,3</td><td>21</td><td>59</td><td>25</td><td>6</td><td>27</td><td>7</td></tr><tr><td>54</td><td>S9XE13</td><td>Uncharacterized protein (C. ferus)</td><td>82,6</td><td>12,16</td><td>32,28</td><td>71,26</td><td>5</td><td>10</td><td>15</td><td>5</td><td>17</td><td>22</td></tr><tr><td>55</td><td>S9Y8C6</td><td>Phosphoglucomutase 1 isoform 3-like protein (C. ferus)</td><td>68,3</td><td>11,37</td><td>18,27</td><td>42,46</td><td>9</td><td>11</td><td>22</td><td>5</td><td>6</td><td>11</td></tr><tr><td>56</td><td>S9Y3S5</td><td>Olfactory receptor (C. ferus)</td><td>108,3</td><td>19,83</td><td>37,9</td><td>43,59</td><td>4</td><td>6</td><td>7</td><td>5</td><td>9</td><td>11</td></tr><tr><td>57</td><td>S9XT33</td><td>Lipopolysaccharide-binding protein (C. ferus)</td><td>47,5</td><td>11,05</td><td>9,43</td><td>16,89</td><td>12</td><td>10</td><td>17</td><td>5</td><td>4</td><td>7</td></tr><tr><td>58</td><td>S9YL21</td><td>Apolipoprotein A-I (C. ferus)</td><td>22,5</td><td>19,05</td><td>68,11</td><td>82,98</td><td>27</td><td>49</td><td>57</td><td>5</td><td>12</td><td>14</td></tr><tr><td>59</td><td>S9Y5X2</td><td>Cell death activator CIDE-A-like protein (C. ferus)</td><td>41,2</td><td>16,09</td><td>22,41</td><td>15,11</td><td>11</td><td>14</td><td>10</td><td>4</td><td>11</td><td>4</td></tr><tr><td>60</td><td>S9XC74</td><td>Osteopontin isoform OPN-c (C. ferus)</td><td>34,6</td><td>8,45</td><td>12,77</td><td>34,67</td><td>8</td><td>11</td><td>23</td><td>4</td><td>10</td><td>24</td></tr><tr><td>61</td><td>T0NLV9</td><td>Epoxide hydrolase 1 (C. ferus)</td><td>54,3</td><td>7,1</td><td>11,53</td><td>12,54</td><td>9</td><td>14</td><td>10</td><td>4</td><td>6</td><td>4</td></tr><tr><td>62</td><td>S9WKD1</td><td>Ribonuclease 4 (C. ferus)</td><td>26,8</td><td>11,56</td><td>25,5</td><td>29,63</td><td>12</td><td>22</td><td>23</td><td>3</td><td>7</td><td>9</td></tr><tr><td>63</td><td>S9WDV3</td><td>Fibrinogen gamma chain isoform gamma-B (C. ferus)</td><td>50,5</td><td>7,22</td><td>33,1</td><td>76,88</td><td>6</td><td>49</td><td>37</td><td>3</td><td>12</td><td>33</td></tr><tr><td>64</td><td>S9XLJ3</td><td>Brain-specific serine protease 4-like protein (C. ferus)</td><td>44,6</td><td>7,52</td><td>8,26</td><td>9,3</td><td>10</td><td>10</td><td>12</td><td>3</td><td>3</td><td>4</td></tr><tr><td>65</td><td>S9WY98</td><td>Sodium/glucose cotransporter 1 (C. ferus)</td><td>78,3</td><td>9,01</td><td>5,72</td><td>14,27</td><td>4</td><td>3</td><td>5</td><td>3</td><td>2</td><td>5</td></tr><tr><td>66</td><td>S9WX48</td><td>Alpha-1-acid glycoprotein (C. ferus)</td><td>22,9</td><td>6,65</td><td>7,25</td><td>29,92</td><td>14</td><td>8</td><td>44</td><td>3</td><td>2</td><td>17</td></tr><tr><td>67</td><td>W5P9V5</td><td>Uncharacterized protein (O. aries)</td><td>85,3</td><td>7,14</td><td>7,74</td><td>8,54</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>68</td><td>T0NLF0</td><td>Vitronectin (C. ferus)</td><td>56,2</td><td>11,18</td><td>7,8</td><td>25,19</td><td>7</td><td>4</td><td>10</td><td>3</td><td>2</td><td>4</td></tr><tr><td>69</td><td>Q0IIG8</td><td>Ras-related protein Rab-18 (B. taurus)</td><td>22,9</td><td>3,49</td><td>15,33</td><td>21,19</td><td>10</td><td>28</td><td>24</td><td>2</td><td>5</td><td>4</td></tr><tr><td>70</td><td>S9YNY9</td><td>Nucleobindin-1 (C. ferus)</td><td>53</td><td>4,54</td><td>48,27</td><td>20,28</td><td>5</td><td>41</td><td>21</td><td>2</td><td>16</td><td>9</td></tr></tbody></table>

Molecular masses (M<sub>r</sub>) of proteins are expressed in kDa, E-value in log, coverage in %. Spectra indicates the number of spectra permitting the identification of proteins. Major proteins identified in excised gel bands after SDS-PAGE are given in bold type.

*abundance of each protein was estimated from spectral count. The number of spectra of C. bactrianus (B) classified the table.

https://doi.org/10.1371/journal.pone.0197026.t002

(±79.98 Da). Likewise, masses of isoforms arising from cryptic splice site usage, usually leading to the loss of the first codon (CAG) of an exon specifying a glutaminyl residue (-128 Da), are easily deduced. A camel mass reference database was thus created for the main milk proteins by combining the data available from *C. dromedarius*, *C. bactrianus*, *C. ferus*, and *Lama glama* milk protein sequences published in UniProtKB (ExPASy SIB Bioinformatics Resource Portal) and the National Centre for Biotechnology Information (NCBI).

To illustrate the efficiency of such an approach, a typical protein profile obtained with a milk from a hybrid camel sampled in Kyzylorda region is given in Fig 4. The analysis of molecular isoforms, identified from mass data, are reported in Table 3, in which experimental and theoretical molecular masses of camel milk proteins are given and confronted. The mass accuracy has allowed distinguishing about 50 protein molecules corresponding to isoforms belonging to nine protein families, eluted from the reverse-phase column as 15 peaks.

In peak I, the two molecular masses (21,157 Da and 21,184 Da) found were associated with glycoforms of κ-CN. The molecular mass of 21,157 Da corresponds to mono-phosphorylated
<figure><img src="image_7.png" /><figcaption>**Fig 4.** LC-ESI-MS profile of clarified crossed camel milk of Kyzylorda region. Nine major milk protein fractions were identified in the following order: peak I and II contained glycosylated ant natural isoforms of κ-CN; peak III: WAP; peaks IV, V: αs1-CN; peak VI: α-LAC, αs1-CN and UP1; peak VII: αs2-CN and UP1; peaks VIII, IX, and X αs2-CN along with UP2 in peak X; peak XI: PGRP and UP2; peak XII: CSA/LPO; peaks XIII and XIV: β-CN, and peak XV: γ2-CN.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0197026.g004

variant A of κ-CN with tri-saccharides ((GaN-Ga-SA2) x 3 or (GaN-Ga) + (GaN-Ga-SA3) x 2, or (GaN-Ga-SA) + (GaN-Ga-SA2) + (GaN-Ga-SA3)). The molecular mass of 21,184 Da was expected to be non-phosphorylated variant B of κ-CN with penta-saccharides ((GaN-Ga) x 3 + (GaN-Ga-SA2) x 2, or (GaN-Ga) + (GaN-Ga-SA) x 4, or (GaN-Ga) x 2 + (GaN-Ga-SA) x 2 + (GaN-Ga-SA2), or (GaN-Ga) x 3 + (GaN-Ga-SA) + (GaN-Ga-SA3)). Peak II contained molecules of which the molecular masses (18,210 Da and 18,236 Da) were identified as non-phosphorylated variant B of κ-CN along with the A variant modified at its N-terminal residue to form a pyro-glutamic acid (pyro-E), which is formed spontaneously by cyclization of the N-terminal E residue. The two molecular masses: 12,564 Da and 12,644 Da, detected in peak III, were assigned to the WAP peptide chain without or with one P group, respectively. Peaks IV, V, and VI were shown to contain αs1-CN. The molecular mass of 23,878 Da observed in peak IV was interpreted as being a short isoform (201-residues) of αs1-CN variant A with 4P groups, arising from exons 13' and 16 skipping events in the mature mRNA during the course of primary transcripts splicing, resulting in deleted sequences (residues E112-Q117 and E155-E162). Despite identification of only one splicing isoform with 4P groups (23,878 Da) in this milk sample, isoforms with 3P and 5P, along with cryptic splice site usage were identified in several other milk samples. Peak V consisted of three relative groups of three masses with sequential increments (s.i.) of 80 Da: 24,547 Da—24,707 Da, 24,675 Da—24,835 Da, and 24,689 Da—24,849 Da. The mass difference (128 Da) between the first and the second group (Table 3) corresponds to the loss of glutaminyl residue 83 (ΔQ83), encoded by the first codon (CAG) of exon 11. As reported previously [38], 24,755 Da was identified as the short isoform (207-residues) of the αs1-CN variant A originating from exon 16 skipping during the course of the primary transcript processing. The mass difference (14 Da) between the second (24,675 Da) and
<table><caption>Table 3. Identification of camel milk protein (hybrid from Kyzylorda region) from observed molecular masses using LC-ESI-MS.</caption><thead><tr><th>Peak</th><th>Ret.Time, min</th><th>Observed M<sub>r</sub>, Da</th><th>TheoreticalM<sub>r</sub>, Da</th><th>Protein description</th><th>UniProt/ NCBI GenBank Accession number</th><th>Intensity</th></tr></thead><tbody><tr><td rowspan="2">I</td><td rowspan="2">4.50</td><td>21,157</td><td>21,158</td><td>κ-CN A, 1P, (GaN-Ga-SA2)x3*, pyro-E</td><td rowspan="2"></td><td>1,361</td></tr><tr><td><strong>21,184</strong></td><td><strong>21,182</strong></td><td><strong>κ-CN B, 0P, (GaN-Ga)x3 + (GaN-Ga-SA2)x2**, pyro-E</strong></td><td><strong>5,810</strong></td></tr><tr><td rowspan="2">II</td><td rowspan="2">18.61</td><td>18,210</td><td>18,210</td><td>κ-CN B, 0P ?</td><td>L0P304</td><td>161</td></tr><tr><td>18,236</td><td>18,235</td><td>κ-CN A, 0P, pyro-E</td><td>P79139</td><td>72</td></tr><tr><td rowspan="2">III</td><td rowspan="2">24.32</td><td>12,564</td><td>12,564</td><td>WAP, 0P</td><td>P09837</td><td>1,756</td></tr><tr><td>12,644</td><td>12,644</td><td>WAP, 1P</td><td></td><td>1,575</td></tr><tr><td>IV</td><td>24.97</td><td>23,878</td><td>23,878</td><td>α<sub>s1</sub>-CN A—short isoform (Δex 16 and Δex 13'), 4P</td><td></td><td>242</td></tr><tr><td rowspan="10">V</td><td rowspan="10">26.23</td><td>24,547</td><td>24,547</td><td>α<sub>s1</sub>-CN C—short isoform (Δex 16), 5P, splice variant (ΔQ83)</td><td></td><td>4,885</td></tr><tr><td><strong>24,627</strong></td><td><strong>24,627</strong></td><td><strong>α<sub>s1</sub>-CN C—short isoform (Δex 16), 6P, splice variant (ΔQ83)</strong></td><td><strong>21,606</strong></td></tr><tr><td>24,707</td><td>24,707</td><td>α<sub>s1</sub>-CN C—short isoform (Δex 16), 7P, splice variant (ΔQ83)</td><td>6,990</td></tr><tr><td>24,675</td><td>24,675</td><td>α<sub>s1</sub>-CN C—short isoform (Δex 16), 5P</td><td>9,441</td></tr><tr><td><strong>24,755</strong></td><td><strong>24,755</strong></td><td><strong>α<sub>s1</sub>-CN C—short isoform (Δex 16), 6P</strong></td><td><strong>47,392</strong></td></tr><tr><td>24,835</td><td>24,835</td><td>α<sub>s1</sub>-CN C—short isoform (Δex 16), 7P</td><td>7,046</td></tr><tr><td>24,689</td><td>24,689</td><td>α<sub>s1</sub>-CN A—short isoform (Δex 16), 5P</td><td>9,748</td></tr><tr><td><strong>24,768</strong></td><td><strong>24,769</strong></td><td><strong>α<sub>s1</sub>-CN A—short isoform (Δex 16), 6P</strong></td><td><strong>50,634</strong></td></tr><tr><td>24,849</td><td>24,849</td><td>α<sub>s1</sub>-CN A—short isoform (Δex 16), 7P</td><td>6,909</td></tr><tr><td rowspan="10">VI</td><td rowspan="10">28.53</td><td><strong>14,430</strong></td><td><strong>14,430</strong></td><td><strong>α-LAC</strong></td><td><strong>P00710</strong></td><td><strong>17,797</strong></td></tr><tr><td>22,939</td><td>n/a</td><td>Uncharacterized protein 1 (UP1)</td><td>n/a***</td><td>2,701</td></tr><tr><td>23,020</td><td>n/a</td><td>UP1+80Da</td><td>n/a</td><td>2,489</td></tr><tr><td>23,099</td><td>n/a</td><td>UP1+160Da</td><td>n/a</td><td>1,079</td></tr><tr><td>25,646</td><td>25,645</td><td>α<sub>s1</sub>-CN C, 6P, splice variant (ΔQ83)</td><td>3,501</td></tr><tr><td>25,693</td><td>25,693</td><td>α<sub>s1</sub>-CN C, 5P</td><td>564</td></tr><tr><td><strong>25,773</strong></td><td><strong>25,773</strong></td><td><strong>α<sub>s1</sub>-CN C, 6P</strong></td><td><strong>7,880</strong></td></tr><tr><td>25,787</td><td><strong>25,787</strong></td><td><strong>α<sub>s1</sub>-CN A, 6P</strong></td><td><strong>O97943-1</strong></td><td><strong>3,472</strong></td></tr><tr><td>21,825</td><td>21,826</td><td>α<sub>s2</sub>-CN, 7P</td><td></td><td>552</td></tr><tr><td rowspan="3">VII</td><td rowspan="3">30.05</td><td><strong>21,906</strong></td><td><strong>21,906</strong></td><td><strong>α<sub>s2</sub>-CN, 8P</strong></td><td><strong>O9794</strong></td><td><strong>5,242</strong></td></tr><tr><td>21,984</td><td>21,986</td><td>α<sub>s2</sub>-CN, 9P</td><td></td><td>403</td></tr><tr><td>23,178</td><td>n/a</td><td>UP1+240Da</td><td>n/a</td><td>1,256</td></tr><tr><td rowspan="2">VIII</td><td rowspan="2">31.11</td><td>21,986</td><td>21,986</td><td>α<sub>s2</sub>-CN, 9P</td><td>O97944</td><td>356</td></tr><tr><td><strong>22,066</strong></td><td><strong>22,066</strong></td><td><strong>α<sub>s2</sub>-CN, 10P</strong></td><td></td><td><strong>4,790</strong></td></tr><tr><td rowspan="2">IX</td><td rowspan="2">33.18</td><td>22,066</td><td>22,066</td><td>α<sub>s2</sub>-CN, 10P</td><td></td><td>148</td></tr><tr><td><strong>22,145</strong></td><td><strong>22,146</strong></td><td><strong>α<sub>s2</sub>-CN, 11P</strong></td><td></td><td><strong>1,964</strong></td></tr><tr><td rowspan="3">X</td><td rowspan="3">35.05</td><td>22,226</td><td>22,226</td><td>α<sub>s2</sub>-CN, 12P</td><td></td><td>894</td></tr><tr><td>23,046</td><td>n/a</td><td>Uncharacterized protein 2 (UP2)</td><td>n/a</td><td>231</td></tr><tr><td><strong>19,143</strong></td><td><strong>19,143</strong></td><td><strong>PGRP</strong></td><td><strong>Q9GK12</strong></td><td><strong>7,207</strong></td></tr><tr><td rowspan="3">X</td><td rowspan="3">37.16</td><td>23,206</td><td>n/a</td><td>UP2+160Da</td><td>n/a</td><td>1,592</td></tr><tr><td>23,286</td><td>n/a</td><td>UP2+240Da</td><td>n/a</td><td>735</td></tr><tr><td>66,481</td><td>66,477</td><td>CSA ?</td><td>XP_010981066.1</td><td>1,096</td></tr><tr><td rowspan="3">XII</td><td rowspan="3">38.09</td><td></td><td></td><td>66,491</td><td>Q9GJW6</td><td>2,663</td></tr><tr><td>66,512</td><td>n/a</td><td>CSA ? LPO?</td><td></td><td></td></tr><tr><td>67,342</td><td>n/a</td><td>CSA ? LPO?</td><td></td><td>1,010</td></tr></tbody></table>

(Continued)
<table><caption>**Table 3.** (Continued)</caption><thead><tr><th>Peak</th><th>Ret.Time, min</th><th>Observed M<sub>r</sub>, Da</th><th>TheoreticalM<sub>r</sub>, Da</th><th>Protein description</th><th>UniProt/ NCBI GenBank Accession number</th><th>Intensity</th></tr></thead><tbody><tr><td rowspan="8">XIII</td><td rowspan="8">40.67</td><td>24,746</td><td>24,745</td><td>β-CN A, 3P, splice variant (ΔQ29)</td><td rowspan="8">A0A077SL35</td><td>2,073</td></tr><tr><td>24,793</td><td>24,792</td><td>β-CN A, 2P</td><td>5,469</td></tr><tr><td>24,825</td><td>24,825</td><td>β-CN A, 4P, splice variant (ΔQ29)</td><td>9,586</td></tr><tr><td>24,873</td><td>24,872</td><td>β-CN A, 3P</td><td>10,177</td></tr><tr><td>24,953</td><td>24,953</td><td>β-CN A, 4P</td><td>84,494</td></tr><tr><td>24,842</td><td>24,842</td><td>β-CN B, 4P, splice variant (ΔQ29)</td><td>10,029</td></tr><tr><td>24,891</td><td>24,890</td><td>β-CN B, 3P</td><td>10,365</td></tr><tr><td>24,970</td><td>24,971</td><td>β-CN B, 4P</td><td>87,973</td></tr><tr><td rowspan="6">XIV</td><td rowspan="6">43.71</td><td>23,878</td><td>23,878</td><td>β-CN A-short isoform (Δ946 Da), 4P, splice variant (ΔQ29)</td><td rowspan="6">Q9TVD0</td><td>707</td></tr><tr><td>23,963</td><td>23,958</td><td>β-CN A-short isoform (Δ946 Da), 5P, splice variant (ΔQ29)</td><td>244</td></tr><tr><td>23,929</td><td>23,926</td><td>β-CN A-short isoform (Δ946 Da), 3P</td><td>438</td></tr><tr><td>24,006</td><td>24,006</td><td>β-CN A-short isoform (Δ946 Da), 4P</td><td>9,026</td></tr><tr><td>23,895</td><td>23,896</td><td>β-CN B-short isoform (Δ946 Da), 4P, splice variant (ΔQ29)</td><td>625</td></tr><tr><td>24,024</td><td>24,024</td><td>β-CN B-short isoform (Δ946 Da), 4P</td><td>5,545</td></tr><tr><td rowspan="2">XV</td><td rowspan="2">47.02</td><td>12,357</td><td>12,358</td><td>γ2-CN A, 0P</td><td rowspan="2"></td><td>1,473</td></tr><tr><td>12,376</td><td>12,376</td><td>γ2-CN B, 0P</td><td>1,065</td></tr></tbody></table>

Major proteins within each peak are in bold. Proteins and isoforms previously described are on grey background.

\*(GaN-Ga-SA2) x 3, or (GaN-Ga) + (GaN-Ga-SA3) x 2, or (GaN-Ga-SA)+(GaN-Ga-SA2)+(GaN-Ga-SA3).

\*(GaN-Ga) x 3 + (GaN-Ga-SA2) x 2, or (GaN-Ga) + (GaN-Ga-SA) x 4, or (GaN-Ga) x 2 + (GaN-Ga-SA) x 2 + (GaN-Ga-SA2), or (GaN-Ga) x 3 + (GaN-Ga-SA) + (GaN-Ga-SA3).

\*\*n/a—not applicable.

https://doi.org/10.1371/journal.pone.0197026.t003

the third (24,689 Da) group is due to the aa substitution E30D reported by Shuiep *et al.* [39] characterizing the C variant. Thus, it is concluded that the third mass group gathers α<sub>s1</sub>-CN short isoforms (207-residues) of variant C, with 5P, 6P and 7P, respectively, described in C. *dromedarius*. While cryptic splice site isoforms (ΔQ83) of variant C, with different phosphorylation levels, were not found in the milk sample shown at Fig 4, they were successfully found in several milk samples. Whereas, α<sub>s1</sub>-CN short isoform was systematically present in all camel milk samples with 5, 6 and 7P (Table 3), by contrast, α<sub>s1</sub>-CN short isoforms of variant C occurred in some milk samples with 4P (24,611 Da) and up to 9P (24,929 Da). Herein, α<sub>s1</sub>-CN short isoforms of variants A and C carrying 6P groups are isoforms with the highest mass signal intensity values 50,634 vs. 47,392, respectively.

Peak VI was more complex to interpret. Masses found in this peak belonged to four different molecular mass groups: 14,430 Da (ascribed to α-LAC), 22,939–23,099 Da (s.i. of 80 Da), 25,646 Da and 25,693–25,773 Da (s.i. of 80 Da), and 25,787 Da. Masses around 23 kDa (22,939–23,099 Da), with a mass increment of two P groups (160 Da), were not referenced to any protein in our database. These findings strongly suggest the existence of an additional uncharacterized phosphorylated protein, namely UP1, which remains to be identified. The third mass group, 25,646 Da and 25,693–25,773 Da, corresponds to a mixture of two long isoforms (214 and 215 aa residues, respectively) of α<sub>s1</sub>-CN variant C with 5P and 6P (25,693–25,773 Da) which differs from variant A by an aa substitution (E30D) in the mature protein [40]. The mass of 25,646 Da corresponds to a 214 aa residues isoform of α<sub>s1</sub>-CN variant C (ΔQ83), with 6P. The last molecular mass (25,787 Da) found in this peak was related to the
mature variant A of α<sub>s1</sub>-CN bearing 6P groups, which is by far much less abundant than the short α<sub>s1</sub>-CN A-6P isoform (intensity of the mass signals: 3,472 vs. 50,634).

The four subsequent peaks (VII, VIII, IX, and X) all contained α<sub>s2</sub>-CN molecules, with
phosphorylation levels ranging between 7P (21,825 Da, peak VII) and 12P (22,226 Da, peak
X). Observed molecular masses of 21,825–21,984 Da were in perfect concordance with those
predicted for α<sub>s2</sub>-CN displaying 7P and 9P, whereas α<sub>s2</sub>-CN with 8P (21,906 Da) was the most
frequent isoform. In addition, the mass of 23,179 Da in peak VII probably corresponds to the
UP1 found in fraction VI with one more P group. Masses ranging between 21,986 and 22,226
Da (s.i. of 80 Da) found in peaks VIII, IX, and X were related to α<sub>s2</sub>-CN variant A with 9P to
12P. These results suggest three more potential phosphorylation sites than reported by Kappe-
ler *et al.* [38] who mentioned a maximum of 9 S residues phosphorylated in camel α<sub>s2</sub>-CN.
More recently, Felfoul *et al.* [37] detected two α<sub>s2</sub>-CN isoforms with 10 and 11P groups in
camel milk. Interestingly, peak X contains a second uncharacterized protein (UP2) with a
molecular mass of 23,046 Da, not referring to any mass in our database for camel milk pro-
teins. Such a mass was found in all camel milk samples analyzed so far (n = 30). This suggests
the possible existence of a further phosphoprotein in camel milk, very likely a CN, since two
putative related isoforms with two (23,206 Da) and three (23,286 Da) additional P groups were
detected in peak XI, in which the most abundant mass found (19,143 Da) was attributed to
PGRP.

In the hybrid from Kyzylorda region (Table 3), masses found in peak XII ranged between 66,481 and 67,342 Da. The most abundant masses 66,481 Da and 66,512 Da might be related to CSA of which the theoretical mass (peptide sequence predicted from the *C. dromedarius* genome, NCBI Accession number XP_010981066.1) is 66,477 Da. The mass differences of 4 Da and 35 Da could be attributed to putative genetic polymorphisms. The molecular weight reported by Felfoul *et al.* [37] from fresh camel milk was estimated as 66,600 Da. However, one cannot exclude that such masses could correspond to LPO depending on cleavage sites of the propeptide, when comparing with bovine LPO and human myeloperoxidase [41].

Molecular masses of 24,793–24,953 Da (s.i. of 80 Da) found in peak XIII, were ascribed to β-CN variant A with 2P, 3P and 4P, first described in the *C. bactrianus*. Molecular masses of 24,891–24,970 Da, which differ from β-CN A-3P and 4P by a 18 Da, correspond to β-CN variant B, first described in *C. dromedarius*. The mass difference of 18 Da between variants A and B is due to the M186I substitution. Isoforms of β-CN with 4P predominate whatever the milk sample and the genetic variant were, with equivalent intensity values of the mass signal for variants A and B, exemplified by a heterozygous hybrid camel: 84,494 vs. 87,973, respectively. In addition, the molecular mass of 24,842 Da, observed in peak XIII, corresponds to a splicing variant of β-CN B-4P. Such an isoform, which was so far considered as typical to the dromedary camel, was also found in hybrids and Bactrian camels. It is due to a cryptic splice site usage leading to the loss of the first codon (CAG) of exon 6, encoding residue Q29 in the protein.

Surprisingly, in the next peak (XIV), molecular masses around 24,000 Da (23,878 Da to 24,024 Da) were observed. Given the elution time and the mass range, these masses were very likely relative to the β-CN fraction, especially since a 18 Da mass differential existing between the pair of molecular masses (24,006 Da and 24,024 Da), is consistent with the occurrence of β-CN variants A and B, in both species. The important mass reduction, ~ 946 Da, relatively to the full-length β-CN, is hypothesized to be due to the cleavage by plasmin of the first seven N-terminal residues (1REKEEFK7) of the mature protein, given that this heptapeptide accounts for 947 Da. Furthermore, molecular masses equivalent to 23,878 Da and 23,895 Da are supposed to originate in the cryptic splice site usage (ΔQ29), previously mentioned.
Finally, in the last peak (XV) mass values 12,357 Da and 12,376 Da again with the mass difference in 18 Da were observed. These masses correspond very likely to camel $\gamma_2$-CN A and B (12,357 Da vs. 12,376 Da, respectively), which are degradation products of $\beta$-CN [42].

This extensive analysis shows that mass accuracy provided by LC-ESI-MS was effective to allow protein identification of most of the protein isoforms by comparison of masses observed experimentally to theoretical molecular masses, and sufficiently powerful to recognize post-translational modifications (PTM) such as phosphorylation of CN, as well as genetic variants and long and short isoforms due to splicing inaccuracies.

## Multiple spliced variants of CSN1S1

To confirm the occurrence of CSN1S1 multiple splice variants, we took advantage of the possibility to extract RNA from milk fat globules to sequence PCR fragments of cDNA encoding $\alpha_{s1}$-CN. Three different CSN1S1 transcripts were found in each species and both genetic variants A and C. The nucleotide sequence of the most frequent variant transcript was shown to be deleted of exon 16, encoding the octapeptide EQAYFHLE. Besides, we also observed an isoform displaying the same sequence in which the first codon of exon 11 was lacking. Finally, a full-length transcript including exon 16 and the first codon of exon 11 was also detected, at a lower concentration.

## Discussion

Given the growing interest in camel milk, due to the health potential of its bioactive components [43] and frequently reported high anti-microbial activity [44], over the past 20 years and even more during the last decade, the milk protein fraction of Camelids, from all around the world has been extensively investigated [10], [11], [36], [37], [39], [40], [45–53]. All these studies have explored, with more or less efficient approaches, the composition of the major milk proteins. However, the molecular diversity of these major proteins had not yet been studied. Then, our main objective was i) to provide, if not a comprehensive, at least an in-depth description of the protein fraction of camel milk; ii) to go further into an extensive analysis of the molecular diversity of major milk proteins from Camelids (*C. dromedarius*, *C. bactrianus*, and hybrids) sampled from different regions in Kazakhstan. For these purposes, different proteomic tools and methodological approaches were applied. For short, up to 391 protein species were identified in cumulating LC-MS/MS analyses of 8 individual *Camelus* milk, and the extensive characterization of CN and whey protein polymorphisms, using LC-ESI-MS, revealed a minimum of 50 molecular species.

## Interspecies in-depth proteomic analysis of camel milk proteins

To our knowledge, the number of proteins identified in this study was relatively higher compared to the numbers reported in previous studies on the camel proteome [10], [11]. The largest camel milk proteome determined so far comprised about 238 proteins including some known camel proteins and heavy-chain immunoglobulins [10]. In this study carried out on *C. dromedarius*, proteins were identified from 2D SDS-PAGE with subsequent matrix-assisted laser desorption/ionization (MALDI) time-of-flight mass spectrometry analysis. However, it should be mentioned that several of the 238 proteins identified matched with the same protein in different species. Hence, at most *ca*. 140 proteins may be considered as unique. By comparison, in the present study a total of 391 unique protein species were determined from LC-MS/MS analyses of *C. bactrianus* (n = 3), *C. dromedarius* (n = 3), and hybrids (n = 2), sampled from three different regions (Atyrau, Shymkent and Kyzylorda). Proteins such as flavin monamine oxidase, perilipin 2, neutrophil gelatinase-associated lipocalin-like protein, brain-
specific serine protease 4-like protein and others, which were not determined previously, were successfully detected. Conversely, about 30 proteins identified by Alhaider and co-workers [10] were not found in our study.

However, as for other mammals, CN represent the major protein fraction of camel milk (80%), among which β-CN is the most abundant [54]. Quantitative analyses performed by Kappeler *et al.* [35] on camel milk CN have demonstrated significant higher amounts of β-CN (15 g/L vs. 10 g/L) compared to the homologous bovine β-CN and significant lower amounts of κ-CN (0.8 g/L vs. 3.5 g/L). Regarding relative proportions, as previously reported [38], α<sub>s1</sub>-, α<sub>s2</sub>-, β- and κ-CN contribute to about 22%, 9.5%, 65%, and 3.5% of total CN, respectively. Taking into account the 30 milk samples analyzed in LC-ESI-MS, relative proportions of individual CN, estimated from the mass signal intensity of each CN family (summing the mass signal of its phosphorylation and splicing isoforms) relatively to the sum of mass signal intensities of all CN families (considering that ionizing properties of caseins and their isoforms are comparable), were 37% α<sub>s1</sub>-CN, 6.1% α<sub>s2</sub>-CN, 53.1% β-CN, and 3.8% κ-CN. These values varied considerably compared to those reported previously by Kappeler *et al.* [38] essentially as far as α<sub>s1</sub>-CN and β-CN are concerned. Whereas α<sub>s1</sub>-CN accounts for 36.1% for *C. bactrianus*, it reaches 37.4% and 37.6% in *C. dromedarius* and hybrids, respectively (Table 4). Percentage of α<sub>s1</sub>-CN calculated in our study was 15% higher than the value reported by Kappeler *et al.* Such an increase is compensated in part by a decrease of 12% of β-CN. The small amount of κ-CN observed is probably underestimated, since most of the highly glycosylated isoforms were not detected. However, this is in agreement with the fact that the size distribution of CN micelles is inversely related to κ-CN content [55], [56], since camel CN micelles are the largest, ranging in size between 280–550 nm [57].

Even though, there are 2 potential phosphorylation sites in κ-CN (S141 and S159) conserved and phosphorylated in sheep and goats [27] only isoforms with a single or no P group in the first chromatographic peak comprising glycosylated isoforms with 3 or 5 carbohydrate motifs were detected. Five glycosylated isoforms of camel κ-CN ranging in size between *ca.* 24 and 25.9 kDa were found in camel milk using 2D SDS-PAGE [50].

In addition, γ₂-CN, a C-terminal product resulting from a highly specific proteolysis of β-
CN by the natural milk protease (plasmin) was successfully found in the milk samples ana-
lyzed. Previously published data suggested that the proportion of γ-CN in total CN fraction is
highest at the beginning and the end of lactation, and in very low yielding animals [56]. The
molecular masses observed in this study (12,357 Da and 12,376 Da) were lower from those pre-
viously observed by Kappeler [38]: 13.9 kDa, 15.7 kDa and 15.75 kDa.

Immune-related proteins such as GlyCAM1, MFGE8 and LTF were detected in camel milk. GlyCAM1, also named lactophorin or PP3 is a cysteine free protein, which belongs to the family of GlyCAM-type molecules [58]. Two splicing variants A and B were distinguished in camel milk [47]. Variant A encoding 137 aa residues has a M<sub>r</sub> of 15.7 kDa, while variant B encoding 122 aa residues has a M<sub>r</sub> of 13.8 kDa. The primary structure of Variant A reveals 54% identity with a protein isolated from bovine milk [59]. Until late, it has been claimed that camel GlyCAM1 is neither glycosylated nor phosphorylated as bovine GlyCAM1. However, Girardet *et al.* [22] suggested the probable existence of one O-glycosylation site (16TDT18) in variant A of which the apparent M<sub>r</sub> was estimated as 22.5 kDa from SDS-PAGE. Using the same approach, two bands were found, in which we identified GlyCAM1 from LC-MS/MS analysis 22 kDa and 10 kDa, corresponding probably to the glycosylated and putatively phosphorylated isoform of GlyCAM1 observed by Girardet *et al.* [22], and to a product of proteolysis, respectively. Surprisingly, no molecular masses corresponding to camel GlyCAM1 A and B were identified by LC-ESI-MS analysis. Likewise, LC-ESI-MS did not permit to detect LTF, even though, SDS-PAGE and LC-MS/MS data confirm its presence in analyzed camel milk
<table><caption>Table 4. Relative proportion of each CN expressed in %, estimated from the mass signal intensity of each CN family relatively to the sum of mass signal intensities of all CN families in the three camel species.</caption><thead><tr><th rowspan="2"></th><th colspan="2">κ-CN</th><th colspan="2">α<sub>s1</sub>-CN</th><th colspan="2">α<sub>s2</sub>-CN</th><th colspan="2">β-CN</th></tr><tr><th>m</th><th>σ</th><th>m</th><th>σ</th><th>m</th><th>σ</th><th>m</th><th>σ</th></tr></thead><tbody><tr><td>Bactrian</td><td>3,09</td><td>1,89</td><td>36,09</td><td>2,33</td><td>7,13</td><td>1,49</td><td>53,68</td><td>2,08</td></tr><tr><td>Dromedary</td><td>3,63</td><td>2,13</td><td>37,39</td><td>3,89</td><td>5,79</td><td>0,98</td><td>53,19</td><td>3,46</td></tr><tr><td>Hybrid</td><td>4,77</td><td>3,01</td><td>37,57</td><td>3,03</td><td>5,25</td><td>1,56</td><td>52,41</td><td>4,18</td></tr></tbody></table>

m = mean.

σ = standard deviation.

https://doi.org/10.1371/journal.pone.0197026.t004

samples. On the other hand, molecular masses ranging between 74,338 Da-79,621 Da could be attributed to camel LTF of which the theoretical mass reported by Kappeler *et al.* [60] for the mature protein (689 aa residues long) without PTM is 75,250 Da. Therefore, the mass difference observed is very likely attributable to PTM. In addition, Konuspayeva *et al.* [61] reported that the level of LTF is affected by seasonal variations.

Elsewhere, MFGM-enriched proteins such as XO, BTN, fatty acid synthase, actin, ras-related protein Rab-18, ADP-ribosylation factor 1, tyrosine-protein kinase, GTP-binding protein SAR1b were identified in Kazakh camel milk samples in accordance with previous results obtained with *C. dromedarius* [62] and *C. bactrianus* [11] milk samples. Surprisingly, whereas BTN was present in all milk samples, it seems to be absent in *C. bactrianus* from Atyrau region. This could be due to the way the band in the electrophoresis gel was cut, since BTN was found in the other seven samples analyzed. Regarding proteins originating from blood, such as serpin A3-1, apolipoprotein A-1, α-1-antitrypsin like protein, α-1-acid glycoprotein, β-2-microglobulin, complement C3-like protein, they were also found in Kazakh camel milks, in agreement with findings of Yang *et al.* [11] reported for Bactrian camels from China. By contrast, as mentioned in the Results section, no trace of CSA was found in Kazakh milk samples from LC-MS/MS analyses, whereas its presence is suspected from LC-ESI-MS.

A heat shock protein (HSPA6 also called HSP70B') occurred at rank 23 amongst the first third of the most represented proteins in Kazakh camel milks (Table 2). Expression of heat shock proteins, including HSP70 is increased during heat stress and involved in defense against dehydration or thermal stress in arid environments [63], [64]. The entire sequence of this protein has been deduced from the nucleotide sequence of a full-length cDNA in *C. dromedarius* [65]. Comprising 643 aa residues, the camel protein, of which the M<sub>r</sub> is 70,543 Da in agreement with the molecular mass estimated from SDS-PAGE, shares a high similarity (94% identity) with cow and pig HSP70.

Against all expectations, peptides with sequence similarity with bovine β-lactoglobulin, the major allergen in bovine milk, were identified in the 8 camel milk samples (Bactrian, dromedary and hybrids) from Kazakhstan, analyzed by LC-MS/MS. The coverage percentage ranged between 30 and 60% in individual milk samples, and reached 71% cumulating all the peptides found. Five peptides related to bovine β-lactoglobulin were also detected by Alhaider *et al.* [10] in camel milk from Saudi Arabia and the United States. Youcef *et al.* [48] revealed a weak cross reaction between dromedary whey proteins and IgG anti bovine β-lactoglobulin. Such findings disagree with the usually admitted notion that β-lactoglobulin is absent in camel milk [50], [66]. Even though we cannot exclude a possible contamination by bovine milk (unlikely with the 8 camel milk samples analyzed by LC-MS/MS) or the presence in camel milk of a Progesterone Associated Endometrial Protein (PAEP) displaying strong similarities with β-
lactoglobulin. However, significant similarities between human PAEP and the peptides having allowed the identification of β-lactoglobulin in *C. bactrianus* milk, were not found.

## Molecular diversity of camel caseins: Genetic polymorphism and alternative splicing

Regarding camel α<sub>s1</sub>-CN, the situation is particularly confusing. Kappeler *et al.* [38] first described two cDNA (short and long), encoding two protein isoforms of 207 and 215 aa, named A and B variants. The A variant corresponds to the short isoform (207 aa), in which the octapeptide 155EQAYFHLE162 encoded by exon 16 was missing, whereas this octapeptide is present in the 215 aa-long isoforms. In our study, two isoforms long and short showing a 1,018 Da mass difference were found, in which the short isoform was the major component (ca. 90%) of total camel α<sub>s1</sub>-CN. Such an alternative splicing event has been first reported in goats [67], sheep [68], [69] and later in lama [70]. In addition, we observed the existence of two distinct genetic variants called A and C, arising from the E30D aa substitution, as previously reported by Shuiep *et al.* [39]. Since, variants A and B described by Kappeler *et al.* [38] displayed a E aa residue in position 30 of the mature peptide chain, it becomes obvious that Kappeler's A and B variants derived in fact from a single allele, of which the primary transcript is subject to exon 16 skipping during the splicing process. In other words, the B variant is nothing other than a splicing variant of a single allele that we propose to call CSN1S1*A.

Recently, Erhardt *et al.* [40] reported in *C. dromedarius* from different regions of Sudan, the existence of a further variant, called D, clearly displaying a different IEF behavior. Excluding this D variant, which was not precisely characterized, there are α<sub>s1</sub>-CN long and short non-allelic isoforms arising from alternative splicing of a single primary transcript and only two perfectly characterized genetic variants A and C resulting from a single G>T nucleotide substitution in exon 4 and leading to E30D aa substitution. This molecular diversity is becoming more complex due to different phosphorylation levels ranging between 5-8P groups (see thereafter) and due to isoforms arising from cryptic splice site usage [67], [71], [72], leading to the loss of a Q residue corresponding to the first codon of exon 11. Results from cDNA sequencing substantiate this.

Electrophoretic and LC-MS analyses as well as cDNA sequencing confirmed that β-CN occurs as two genetic variants A and B, with the aa substitution M186I (yielding a -18 Da mass difference). The most frequent form of β-CN had 4P groups, one P group more than reported for Somali, Turkana and Pakistani camels by Kappeler *et al.* [38]. Surprisingly, in Kazakh populations, a second series of β-CN components with lower molecular masses (mass difference: -946 Da), relatively to the full-length β-CN were found. This phenomenon, observed with both genetic variants, might be due to the cleavage by plasmin of the first seven N-terminal residues (REKEEFK) of the mature protein. A mass difference of 947 Da was observed between the native full-length protein with 4P (24,953 Da and 24,971 Da for A and B variants, respectively) and the plasmin cleavage product at the same phosphorylation level (24,006 Da and 24,024 Da for A and B variants, respectively). The occurrence of a K residue in position 7 of the mature β-CN does not occur in any other species, of which the N-terminal sequence is known [27]. However, our results strongly suggest that the peptide bond 7K-T8 is sensible to plasmin that is, like trypsin, a serine protease. Indeed, REKEEFK was present amongst tryptic peptides identified in LC-MS/MS analysis.

There is another even less probable possibility, involving the deletion of exon 5 that encodes 8 aa residues (ESITHINK for a mass of 923 Da), since a similar event was previously characterized from mare [73] and donkey [74] milks. However, sequencing of camel β-CN cDNA has not revealed any deletion in the mRNA encoding this protein (results not shown),
consistently with Kappeler *et al.* [38] who only reported a full-length sequence for β-CN, conversely to α<sub>s1</sub>-CN. Since in our study we were not able to provide any further confirmation of the presence of shorter mRNA of camel β-CN in which exon 5 is spliced out, we give preference to the cleavage by plasmin of the first seven N-terminal residues of β-CN rather than an alternative splicing process.

Surprisingly, two so far uncharacterized proteins (UP1 and UP2) with molecular masses around 23,000 Da and different phosphorylation levels were observed, suggesting they are possibly proteins related to CN. However, to prove this hypothesis further research for in depth characterization of these proteins is necessary.

## Post-translational modifications of milk proteins: Phosphorylation of caseins

Among the various approaches developed in proteomics, electrospray ionization (ESI) mass spectrometry (MS) is eminently suitable for studying PTM, including phosphorylation and glycosylation, since the technique provides molecular mass determination of native proteins. Phosphorylation of proteins is one of the most frequent PTM in eukaryotic cells. It has become a common knowledge that phosphorylation of CN occurs at S or T aa residues in tripeptide sequences S/T-X-A where X represents any aa residue and A is an acidic residue [75]. This consensus sequence is recognized by FAM20C, a Golgi CN-kinase, which phosphorylates secreted phosphoproteins, including both CN and members of the small integrin-binding ligand N-linked glycoproteins (SIBLING) protein family, which modulate biomineralization [76]. Each phosphorylation event adds 79.98 Da to the molecular mass of the peptide chain [77]. It was predicted with high confidence 8 probably phosphorylated S residues in α<sub>s1</sub>-CN (S18, S68, S70, S71, S72, S73, S193, and S202), 9 potential phosphorylated S residues in α<sub>s2</sub>-CN (S8, S9, S10, S32, S53, S108, S110, S113, and S121), 4 S residues in β-CN (S15, S17, S18, and S19), and 2 S residues in κ-CN (S141 and S159). However, up to 9P residues per α<sub>s1</sub>-CN molecule were observed whatever the genetic variant is. Theoretically, given the S/T-X-A consensus rule, there are 4 T residues that could be phosphorylated (T55, T80, T153, and T196), leading to a maximum of 12 P groups per molecule. Therefore, we can put forward that at least one of the four T residues is phosphorylated in the α<sub>s1</sub>-CN-9P.

With 11 potentially phosphorylated aa residues matching the S/T-X-A motif (Fig 5), camel α<sub>s2</sub>-CN displays the highest phosphorylation level, in agreement with Felfoul *et al.* who reported recently 11P groups [37]. To reach such a phosphorylation level, besides the nine SerP, two putative ThrP (T118 and T132) have to be phosphorylated. In all the Kazakh milk

<figure><img src="image_8.png" /><figcaption>**Fig 5.** Amino acid sequence of mature camel α<sub>s2</sub>-CN with potential phosphorylation sites. Seryl and Threonyl residues matching the S/T-X-A motif are in red and blue, respectively, and underlined. Threonyl residues matching the S/T-X-X-A motif are in green and underlined.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0197026.g005
samples analyzed in LC-ESI-MS we found α<sub>s2</sub>-CN with 12 P groups, as the molecular mass of 22,226 Da observed corresponds to the mass of the peptide backbone (21,266 Da) increased by 960 Da, a mass increment which coincides with 12 P groups. That means that at least another S/T residue that does not match with the canonic sequence recognized by the mammary kinase(s), is potentially phosphorylated. According to Allende *et al.* [78] the sequence S/T-X-X-A follow-through with the minimum requirements for phosphorylation by the CN-kinase II (CK2). It is critical to highlight in this regard that E or D in this site can be replaced by SerP or ThrP. Two T residues, namely T39 and T129 in the camel α<sub>s2</sub>-CN fully meet the requirements of the above-mentioned motif (Fig 5) and could be phosphorylated. Such an event is the only hypothesis to reach 12P for camel α<sub>s2</sub>-CN. Since these two kinases are very likely secreted, the idea that phosphorylation at T39/T129 may occur in the extracellular environment cannot be excluded. This warrants further investigation. Fam20C, which is very likely the major secretory pathway protein kinase [79], might be responsible for the phosphorylation of S and T residues within S/T-X-A motif, whereas a CK2-type kinase might be responsible for phosphorylation of T residue within an S/T-X-X-A motif. This is in agreement with the hypothesis put forward by Bijl *et al.* [80] and Fang *et al.* [81], who suggest from phenotypic correlations and hierarchical clustering the existence of at least 2 regulatory systems for phosphorylation of α<sub>s</sub>-CN. Elsewhere, bovine milk osteopontin which is a multiphosphorylated glycoprotein also found in bone, was shown to contain 27 SerP and one ThrP [82]. Twenty five SerP and one ThrP were located in S/T-X-E/S(P)/D motifs, whereas two SerP were found in the sequence S-X-X-E/S(P).

**Conclusions**

In this study, six main findings combining proven proteomic and molecular biology approaches are provided. The first one is an enhancing of our knowledge of camel milk protein composition. The second one is deciphering the extreme complexity of camel CN fraction due to PTM (phosphorylation) and splicing events (exon skipping and cryptic splice site usage). The third finding is the detection of two unknown proteins, UP1 and UP2 that remain to be characterized. In addition, we provide results substantiating the possible existence of a camel β-lactoglobulin. However, this result requires further investigation, currently in progress in the laboratory. Afterwards, we report for the first time the presence of α<sub>s2</sub>-CN-12P, and short iso-forms of β-CN probably arising from proteolysis by plasmin, the natural protease of milk. The ultimate finding is the demonstration that genetic variants, which hitherto seemed specific to a species (β-CN A for Bactrian and β-CN B for dromedary), are in fact present in both *dromedarius* and *bactrianus*.

**Supporting information**

S1 Table. Functional groups of proteins identified in camel milk by LC-MS/MS.
(XLSX)

S1 Text. Permission to publish content under CC-BY license.
(PDF)

**Acknowledgments**

The study was carried out within the Bolashak International Scholarship of the first author, funded by the JSC «Center for International Programs» (Kazakhstan). The research was partly supported by a grant from the Ministry of Education and Science of the Republic of Kazakhstan under name "Proteomic investigation of camel milk" #1729/GF4, which is duly
appreciated. The authors thank all Kazakhstani camel milk farms for rendering help in sample
collection, as well as PAPPSO and @BRIDGe teams at INRA (Jouy-en-Josas, France) for pro-
viding necessary facilities and technical support.

**Author Contributions**

**Conceptualization:** Alma Ryskaliyeva, Bernard Faye, Gaukhar Konuspayeva, Patrice Martin.

**Data curation:** Alma Ryskaliyeva, Céline Henry, Guy Miranda.

**Formal analysis:** Alma Ryskaliyeva, Céline Henry, Guy Miranda.

**Funding acquisition:** Gaukhar Konuspayeva.

**Investigation:** Alma Ryskaliyeva, Bernard Faye, Gaukhar Konuspayeva, Patrice Martin.

**Methodology:** Guy Miranda.

**Project administration:** Patrice Martin.

**Resources:** Patrice Martin.

**Supervision:** Bernard Faye, Gaukhar Konuspayeva, Patrice Martin.

**Writing – original draft:** Alma Ryskaliyeva, Patrice Martin.

**Writing – review & editing:** Céline Henry, Guy Miranda, Bernard Faye, Gaukhar Konuspayeva.

**References**

1. FAO, "FAOSTAT," *Food and Agriculture Organization of the United Nations*, 2017. [Online]. Available: http://www.fao.org/faostat/en/#home.

2. Mohandesan E. et al., "Combined hybridization capture and shotgun sequencing for ancient DNA analysis of extinct wild and domestic dromedary camel," *Mol. Ecol. Resour.*, vol. 17, no. 2, pp. 300–313, 2017. https://doi.org/10.1111/1755-0998.12551 PMID: 27289015

3. Nurseitova M., Konuspayeva G., and Jurjanz S., "Comparison of dairy performances between dromedaries, bactrian and crossbred camels in the conditions of South Kazakhstan," *Emirates J. Food Agric.*, vol. 26, no. 4, pp. 366–370, 2014.

4. Agrawal R. P. et al., "Effect of camel milk on glycemic control, risk factors and diabetes quality of life in type-1 diabetes: A randomised prospective controlled study," *J. Camel Pract. Res.*, vol. 10, no. 1, pp. 45–50, 2003.

5. Manaer T., Yu L., Zhang Y., Xiao X. J., and Nabi X. H., "Anti-diabetic effects of shubat in type 2 diabetic rats induced by combination of high-glucose-fat diet and low-dose streptozotocin," *J. Ethnopharmacol.*, vol. 169, pp. 269–274, 2015. https://doi.org/10.1016/j.jep.2015.04.032 PMID: 25922265

6. Sboui A., Khorchani T., Djegham M., Agrebi A., Elhatmi H., and Belhadj O., "Anti-diabetic effect of camel milk in alloxan-induced diabetic dogs: A dose-response experiment," *J. Anim. Physiol. Anim. Nutr. (Berl)*, vol. 94, no. 4, pp. 540–546, 2010.

7. Al-Ayadhi L. Y. and Elamin N. E., "Camel milk as a potential therapy as an antioxidant in Autism Spectrum Disorder (ASD)," *Evid. Based. Complement. Alternat. Med.*, vol. 2013, p. 602834, 2013. https://doi.org/10.1155/2013/602834 PMID: 24069051

8. El-Fakharany E. M. et al., "Influence of camel milk on the hepatitis C virus burden of infected patients," *Exp. Ther. Med.*, vol. 13, no. 4, pp. 1313–1320, 2017. https://doi.org/10.3892/etm.2017.4159 PMID: 28413471

9. Korashy H. M., Maayah Z. H., Abd-Allah A. R., El-Kadi A. O. S., and Alhaider A. a, "Camel milk triggers apoptotic signaling pathways in human hepatoma HepG2 and breast cancer MCF7 cell lines through transcriptional mechanism.," *J. Biomed. Biotechnol.*, vol. 2012, pp. 1–9, 2012. https://doi.org/10.1155/2012/728342 PMID: 21836813

10. Alhaider A. et al., "Through the eye of an electrospray needle: Mass spectrometric identification of the major peptides and proteins in the milk of the one-humped camel (Camelus dromedarius)," *J. Mass Spectrom.*, vol. 48, no. 7, pp. 779–794, 2013. https://doi.org/10.1002/jms.3213 PMID: 23832934
11. Yang Y., Bu D., Zhao X., Sun P., Wang J., and Zhou L., "Proteomic analysis of cow, yak, buffalo, goat and camel milk whey proteins: Quantitative differential expression patterns," *J. Proteome Res.*, vol. 12, no. 4, pp. 1660–1667, 2013. https://doi.org/10.1021/pr301001m PMID: 23464874

12. Hsieh C. C., Hernández-Ledesma B., Fernández-Tomé S., Weinborn V., Barile D., and De Moura Bell J. M. L. N., "Milk proteins, peptides, and oligosaccharides: Effects against the 21st century disorders," *BioMed Research International*, vol. 2015. 2015.

13. Mati A., Senoussi-Ghezali C., Si Ahmed Zennia S., Almi-Sebbane D., El-Hatmi H., and Girardet J. M., "Dromedary camel milk proteins, a source of peptides having biological activities–A review," *International Dairy Journal*, vol. 73. pp. 25–37, 2017.

14. Davoodi S. H. et al., "Health-related aspects of milk proteins," *Iranian Journal of Pharmaceutical Research*, vol. 15, no. 3. pp. 573–591, 2016. PMID: 27980594

15. Legrand D., Elass E., Pierce A., and Mazurier J., "Lactoferrin and host defence: An overview of its immuno-modulating and anti-inflammatory properties," *BioMetals*, vol. 17, no. 3. pp. 225–229, 2004. PMID: 15222469

16. Habib H. M., W. H. Ibrahim, Schneider-Stock R., and Hassan H. M., "Camel milk lactoferrin reduces the proliferation of colorectal cancer cells and exerts antioxidant and DNA damage inhibitory activities," *Food Chem.*, vol. 141, no. 1, pp. 148–152, 2013. https://doi.org/10.1016/j.foodchem.2013.03.039 PMID: 23768340

17. Kanwar J. R. et al., "Multifunctional iron bound lactoferrin and nanomedicinal approaches to enhance its bioactive functions," *Molecules*, vol. 20, no. 6. pp. 9703–9731, 2015. https://doi.org/10.3390/molecules20069703 PMID: 26016555

18. Kappeler S. R., Heuberger C., Farah Z., and Puhan Z., "Expression of the peptidoglycan recognition protein, PGRP, in the lactating mammary gland.," *J. Dairy Sci.*, vol. 87, no. 8, pp. 2660–8, 2004. https://doi.org/10.3168/jds.S0022-0302(04)73392-5 PMID: 15328291

19. Sharma P. et al., "Structural basis of recognition of pathogen-associated molecular patterns and inhibition of proinflammatory cytokines by camel peptidoglycan recognition protein," *J. Biol. Chem.*, vol. 286, no. 18, pp. 16208–16217, 2011. https://doi.org/10.1074/jbc.M111.228163 PMID: 21454594

20. Chace Tydell C., Yount N., Tran D., Yuan J., and Selsted M. E., "Isolation, characterization, and antimicrobial properties of bovine oligosaccharide-binding protein. A microbicidal granule protein of eosinophils and neutrophils," *J. Biol. Chem.*, vol. 277, no. 22, pp. 19658–19664, 2002. https://doi.org/10.1074/jbc.M200659200 PMID: 11880375

21. Dziarski R., Kashyap D. R., and Gupta D., "Mammalian Peptidoglycan Recognition Proteins Kill Bacteria by Activating Two-Component Systems and Modulate Microbiome and Inflammation," *Microb. Drug Resist.*, vol. 18, no. 3, pp. 280–285, 2012. https://doi.org/10.1089/mdr.2012.0002 PMID: 22432705

22. Girardet J. M., Saulnier F., Gaillard J. L., Ramet J. P., and Humbert G., "Camel (camelus dromedarius) milk PP3: evidence for an insertion in the amino-terminal sequence of the camel milk whey protein.," *Biochem. Cell Biol.*, vol. 78, no. 1, pp. 19–26, 2000. PMID: 10735560

23. Hennighausen L. G. and Sippel A. E., "Mouse whey acidic protein is a novel member of the family of 'four-disulfide core' proteins," *Nucleic Acids Res.*, vol. 10, no. 8, pp. 2677–2684, 1982. PMID: 6896234

24. Bouchard D., Morisset D., Bourbonnais Y., and Tremblay G. M., "Proteins with whey-acidic-protein motifs and cancer," *Lancet Oncology*, vol. 7, no. 2. pp. 167–174, 2006. https://doi.org/10.1016/S1470-2045(06)70579-4 PMID: 16455481

25. D. J. McMahon and B. S. Oommen, "Casein micelle structure, functions, and interactions," in *Advanced Dairy Chemistry: Volume 1A: Proteins: Basic Aspects*, 4th Edition, 2013, pp. 185–209.

26. Sakono M., Motomura K., Maruyama T., Kamiya N., and Goto M., "Alpha casein micelles show not only molecular chaperone-like aggregation inhibition properties but also protein refolding activity from the denatured state," *Biochem. Biophys. Res. Commun.*, vol. 404, no. 1, pp. 494–497, 2011. https://doi.org/10.1016/j.bbrc.2010.12.009 PMID: 21144837

27. P. Martin, C. Cebo, and G. Miranda, "Interspecies comparison of milk proteins: Quantitative variability and molecular diversity," in *Advanced Dairy Chemistry: Volume 1A: Proteins: Basic Aspects*, 4th Edition, 2013, pp. 387–429.

28. Balteanu V. A., Carsai T. C., and Vlaic A., "Identification of an intronic regulatory mutation at the buffalo os1-casein gene that triggers the skipping of exon 6," *Mol. Biol. Rep.*, vol. 40, no. 7, pp. 4311–4316, 2013. https://doi.org/10.1007/s11033-013-2518-2 PMID: 23640099

29. Bradford M. M., "A rapid and sensitive method for the quantitation of microgram quantities of protein using the principle of protein dye binding," *Anal. Biochem.*, vol. 72, pp. 248–254, 1976. PMID: 942051

30. Sds-page L. and Sds T., "Laemmli buffer Background," *Nature*, pp. 5–6, 1970.
31. Saadaoui B., Bianchi L., Henry C., Miranda G., Martin P., and Cebo C., "Combining proteomic tools to characterize the protein fraction of llama (Lama glama) milk," *Electrophoresis*, vol. 35, no. 10, pp. 1406-1418, 2014. https://doi.org/10.1002/elp.s.201300383 PMID: 24519815

32. Brenaut P., Bangera R., Bevilacqua C., Rebours E., Cebo C., and Martin P., "Validation of RNA isolated from milk fat globules to profile mammary epithelial cell expression during lactation and transcriptional response to a bacterial infection," *J. Dairy Sci.*, vol. 95, no. 10, pp. 6130-6144, 2012. https://doi.org/10.3168/jds.2012-5604 PMID: 22921620

33. Zhao D. B., Bai Y. H., and Niu Y. W., "Composition and characteristics of Chinese Bactrian camel milk," *Small Ruminant Research*, vol. 127, pp. 58-67, 2015.

34. Konuspayeva G., Faye B., and Loiseau G., "The composition of camel milk: A meta-analysis of the literature data," *Journal of Food Composition and Analysis*, vol. 22, no. 2, pp. 95-101, 2009.

35. Kappeler S. R., Farah Z., and Puhan Z., "5'-Flanking Regions of Camel Milk Genes Are Highly Similar to Homologue Regions of Other Species and Can be Divided into Two Distinct Groups," *J. Dairy Sci.*, vol. 86, no. 2, pp. 498-508, 2003. https://doi.org/10.3168/jds.S0022-0302(03)73628-5 PMID: 12647956

36. Merin U. et al., "A comparative study of milk serum proteins in camel (Camelus dromedarius) and bovine colostrum," *Livest. Prod. Sci.*, vol. 67, no. 3, pp. 297-301, 2001.

37. Felfoul I., Jardin J., Gaucheron F., Attia H., and Ayadi M. A., "Proteomic profiling of camel and cow milk proteins under heat treatment," *Food Chem.*, vol. 216, pp. 161-169, 2017. https://doi.org/10.1016/j.foodchem.2016.08.007 PMID: 27596405

38. Kappeler S., Farah Z., and Puhan Z., "Sequence analysis of Camelus dromedarius milk caseins.," *J. Dairy Res.*, vol. 65, no. 2, pp. 209-222, 1998. PMID: 9627840

39. Shuiep E. T. S., Giambra I. J., El Zubeir I. E. Y. M., and Erhardt G., "Biochemical and molecular characterization of polymorphisms of αs1-casein in Sudanese camel (Camelus dromedarius) milk," *Int. Dairy J.*, vol. 28, no. 2, pp. 88-93, 2013.

40. Erhardt G. et al., "Alpha S1-casein polymorphisms in camel (Camelus dromedarius) and descriptions of biological active peptides and allergenic epitopes," *Trop. Anim. Health Prod.*, vol. 48, no. 5, pp. 879-887, 2016. https://doi.org/10.1007/s11250-016-0997-6 PMID: 26922739

41. Dull T. J., Uyeda C., Strosberg A. D., Nedwin G., and Seilhamer J. J., "Molecular cloning of cDNAs encoding bovine and human lactoperoxidase," *DNA Cell Biol.*, vol. 9, no. 7, pp. 499-509, 1990. https://doi.org/10.1089/dna.1990.9.499 PMID: 2222811

42. Beg O. U., von Bahr-Lindström H., Zaidi Z. H., and Jörnvall H., "Characterization of a camel milk protein rich in proline identifies a new β-casein fragment," *Regul. Pept.*, vol. 15, no. 1, pp. 55-61, 1986. PMID: 3763959

43. Al haj O. A. and Al Kanhal H. A., "Compositional, technological and nutritional aspects of dromedary camel milk," *International Dairy Journal*, vol. 20, no. 12. pp. 811-821, 2010.

44. E. I. El-Agamy, "Bioactive Components in Camel Milk," in *Bioactive Components in Milk and Dairy Products*, 2009, pp. 159-194.

45. El-Agamy E. I., Nawar M., Shamsia S. M., Awad S., and Haenlein G. F. W., "Are camel milk proteins convenient to the nutrition of cow milk allergic children?," *Small Rumin. Res.*, vol. 82, no. 1, pp. 1-6, 2009.

46. Wangoh J., Farah Z., and Puhan Z., "Iso-electric focusing of camel milk proteins," *Int. Dairy J.*, vol. 8, no. 7, pp. 617-621, 1998.

47. Kappeler S., Farah Z., and Puhan Z., "Alternative Splicing of Lactophorin mRNA from Lactating Mammary Gland of the Camel (Camelus dromedarius)," *J Dairy Sci*, vol. 82, no. November 1999, pp. 2084-2093, 1999. https://doi.org/10.3168/jds.S0022-0302(99)75450-0 PMID: 10531593

48. Youcef N. et al., "Cross reactivity between dromedary whey proteins and IgG anti bovine ??-lactalbumin and anti bovine ??-lactoglobulin," *Am. J. Appl. Sci.*, vol. 6, no. 8, pp. 1448-1452, 2009.

49. Ereifej K. I., Alu'datt M. H., Alkhalidy H. A., Alli I., and Rababah T., "Comparison and characterisation of fat and protein composition for camel milk from eight Jordanian locations," *Food Chem.*, vol. 127, no. 1, pp. 282-289, 2011.

50. Hinz K., O'Connor P. M., Huppertz T., Ross R. P., and Kelly A. L., "Comparison of the principal proteins in bovine, caprine, buffalo, equine and camel milk," *J. Dairy Res.*, vol. 79, no. 2, pp. 185-191, 2012. https://doi.org/10.1017/S0022029912000015 PMID: 22365180

51. Salmen S. H., Abu-Tarboush H. M., Al-Saleh A. A., and Metwalli A. A., "Amino acids content and electrophoretic profile of camel milk casein from different camel breeds in Saudi Arabia," *Saudi J. Biol. Sci.*, vol. 19, no. 2, pp. 177-183, 2012. https://doi.org/10.1016/j.sjbs.2011.12.002 PMID: 23961177

52. Ochirkhuyag B., Chobert J. M., Dalgallarrondo M., Choiset Y., and Haertlé T., "Characterization of caseins from Mongolian yak, khainak, and bactrian camel'n," *Lait*, vol. 77, no. 5, pp. 601-613, 1997.
53. Konuspayeva G., Faye B., Loiseau G., and Levieux D., "Lactoferrin and immunoglobulin contents in camel's milk (Camelus bactrianus, Camelus dromedarius, and Hybrids) from Kazakhstan.," *J. Dairy Sci.*, vol. 90, no. 1, pp. 38–46, 2007. https://doi.org/10.3168/jds.S0022-0302(07)72606-1 PMID: 17183073

54. Pauciullo A., Giambra I. J., Iannuzzi L., and Erhardt G., "The β-casein in camels: Molecular characterization of the CSN2 gene, promoter analysis and genetic variability," *Gene*, vol. 547, no. 1, pp. 159–168, 2014. https://doi.org/10.1016/j.gene.2014.06.055 PMID: 24973699

55. Bijl E., de Vries R., van Valenberg H., Huppertz T., and van Hooijdonk T., "Factors influencing casein micelle size in milk of individual cows: Genetic variants and glycosylation of κ-casein," *Int. Dairy J.*, vol. 34, no. 1, pp. 135–141, 2014.

56. Ostersen S., Foldager J., and Hermansen J. E., "Effects of stage of lactation, milk protein genotype and body condition at calving on protein composition and renneting properties of bovine milk.," *J. Dairy Res.*, vol. 64, no. 2, pp. 207–219, 1997. PMID: 9161914

57. Bornaz S., Sahli A., Attalah A., and Attia H., "Physicochemical characteristics and renneting properties of camels' milk: A comparison with goats', ewes' and cows' milks," *Int. J. Dairy Technol.*, vol. 62, no. 4, pp. 505–513, 2009.

58. Beg O. U., von Bahr-Lindström H., Zaidi Z. H., and Jörnvall H., "Characterization of a heterogeneous camel milk whey non-casein protein," *FEBS Lett.*, vol. 216, no. 2, pp. 270–274, 1987. PMID: 3495459

59. Sørensen E. S. and Petersen T. E., "Purification and characterization of three proteins isolated from the proteose peptone fraction of bovine milk," *J. Dairy Res.*, vol. 60, pp. 189–197, 1993. PMID: 8320368

60. Kappeler S. R., Ackermann M., Farah Z., and Puhan Z., "Sequence analysis of camel (Camelus dromedarius) lactoferrin," *Int. Dairy J.*, vol. 9, no. 7, pp. 481–486, 1999.

61. Konuspayeva G., Serikbayeva A., Loiseau G., Narmuratova M., and Faye B., "Lactoferrin of camel milk of Kazakhstan," in *Desertification Combat and Food Safety: the Added Value of Camel Producers*, 2005, vol. 362, pp. 158–167.

62. Saadaoui B., Henry C., Khorchani T., Mars M., Martin P., and Cebo C., "Proteomics of the milk fat globule membrane from Camelus dromedarius," *Proteomics*, vol. 13, no. 7, pp. 1180–1184, 2013. https://doi.org/10.1002/pmic.201200113 PMID: 23349047

63. Sharma S. et al., "Effect of melatonin administration on thyroid hormones, cortisol and expression profile of heat shock proteins in goats (Capra hircus) exposed to heat stress," *Small Rumin. Res.*, vol. 112, no. 1–3, pp. 216–223, 2013.

64. Rhoads R. P., Baumgard L. H., Suagee J. K., and Sanders S. R., "Nutritional Interventions to Alleviate the Negative Consequences of Heat Stress," *Adv. Nutr. An Int. Rev. J.*, vol. 4, no. 3, pp. 267–276, 2013.

65. Elrobh M. S., Alanazi M. S., Khan W., Abduljaleel Z., Al-Amri A., and Bazzi M. D., "Molecular cloning and characterization of cDNA encoding a putative stress-induced heat-shock protein from Camelus dromedarius," *Int. J. Mol. Sci.*, vol. 12, no. 7, pp. 4214–36, 2011. https://doi.org/10.3390/ijms12074214 PMID: 21845074

66. Restani P. et al., "Cross-reactivity between milk proteins from different animal species.," *Clin. Exp. Allergy*, vol. 29, no. 7, pp. 997–1004, 1999. PMID: 10383602

67. Leroux C., Mazure N., and Martin P., "Mutations away from splice site recognition sequences might cismodulate alternative splicing of goat α(s1)-casein transcripts. Structural organization of the relevant gene," *J. Biol. Chem.*, vol. 267, no. 9, pp. 6147–6157, 1992. PMID: 1372900

68. Chianese L., Garro G., Mauriello R., Laezza P., Ferranti P., and Addeo F., "Occurrence of five as1-casein variants in ovine milk," *J. Dairy Res.*, vol. 63, pp. 49–59, 1996. PMID: 8655742

69. Ferranti P. et al., "Primary structure of ovine alpha s1-caseins: localization of phosphorylation sites and characterization of genetic variants A, C and D.," *J. Dairy Res.*, vol. 62, no. 2, pp. 281–296, 1995. PMID: 7601973

70. Pauciullo A. and Erhardt G., "Molecular characterization of the llamas (Lama glama) casein cluster genes transcripts (CSN1S1, CSN2, CSN1S2, CSN3) and regulatory regions," *PLoS One*, vol. 10, no. 4, 2015.

71. Boumahrou N. et al., "Evolution of major milk proteins in Mus musculus and Mus spretus mouse species: a genoproteomic analysis.," *BMC Genomics*, vol. 12, no. 1, p. 80, 2011.

72. Ferranti P., Lilla S., Chianese L., and Addeo F., "Alternative nonallelic deletion is constitutive of ruminant α(s1)-casein," *J. Protein Chem.*, vol. 18, no. 5, pp. 595–602, 1999. PMID: 10524777

73. Miranda G., Mahé M. F., Leroux C., and Martin P., "Proteomic tools characterize the protein fraction of Equidae milk," in *Proteomics*, 2004, vol. 4, no. 8, pp. 2496–2509. https://doi.org/10.1002/pmic.200300765 PMID: 15274143
74. Cunsolo V., Saletti R., Muccilli V., Gallina S., Di Francesco A., and Foti S., "Proteins and bioactive peptides from donkey milk: The molecular basis for its reduced allergenic properties," *Food Research International*, vol. 99. pp. 41–57, 2017. https://doi.org/10.1016/j.foodres.2017.07.002 PMID: 28784499

75. Mercier J. C., "Phosphorylation of caseins, present evidence for an amino acid triplet code posttranslationally recognized by specific kinases," *Biochimie*, vol. 63, no. 1. pp. 1–17, 1981. PMID: 7011421

76. Ishikawa H. O., Xu A., Ogura E., Manning G., and Irvine K. D., "The raine syndrome protein FAM20C is a golgi kinase that phosphorylates bio-mineralization proteins," *PLoS One*, vol. 7, no. 8, 2012.

77. Larsen M. R., Trelle M. B., Thingholm T. E., and Jensen O. N., "Analysis of posttranslational modifications of proteins by tandem mass spectrometry," *BioTechniques*, vol. 40, no. 6. pp. 790–798, 2006. PMID: 16774123

78. Allende J. E. and Allende C. C., "Protein kinases. 4. Protein kinase CK2: an enzyme with multiple substrates and a puzzling regulation," *FASEB J.*, vol. 9, no. 5, pp. 313–323, 1995. PMID: 7896000

79. Tagliabracci V. S. et al., "A Single Kinase Generates the Majority of the Secreted Phosphoproteome," *Cell*, vol. 161, no. 7, pp. 1619–1632, 2015. https://doi.org/10.1016/j.cell.2015.05.028 PMID: 26091039

80. Bijl E., van Valenberg H. J. F., Huppertz T., van Hooijdonk A. C. M., and Bovenhuis H., "Phosphorylation of αS1-casein is regulated by different genes," *J. Dairy Sci.*, vol. 97, no. 11, pp. 7240–7246, 2014. https://doi.org/10.3168/jds.2014-8061 PMID: 25200775

81. Fang Z. H., Visker M. H. P. W., Miranda G., Delacroix-Buchet A., Bovenhuis H., and Martin P., "The relationships among bovine αS-casein phosphorylation isoforms suggest different phosphorylation pathways," *J. Dairy Sci.*, vol. 99, no. 10, pp. 8168–8177, 2016. https://doi.org/10.3168/jds.2016-11250 PMID: 27522420

82. Sørensen E. S., Petersen T. E., and Højrup P., "Posttranslational modifications of bovine osteopontin: Identification of twenty-eight phosphorylation and three O-glycosylation sites," *Protein Sci.*, vol. 4, no. 10, pp. 2040–2049, 1995. https://doi.org/10.1002/pro.5560041009 PMID: 8535240